The numerical solution of diffusion problems in strongly heterogeneous nonisotropic materials
Abstract
A new secondorder finitedifference algorithm for the numerical solution of diffusion problems in strongly heterogeneous and nonisotropic media is constructed. On problems with rough coefficients or highly nonuniform grids, the new algorithm is superior to all other algorithms we have compared it with. For problems with smooth coefficients on smooth grids, the method is comparable with other secondorder methods. The new algorithm is formulated for logically rectangular grids and is derived using the supportoperators method. A key idea in deriving the method was to replace the usual inner product of vector functions by an inner product weighted by the inverse of the material properties tensor and to use the flux operator, defined as the material properties tensor times the gradient, rather than the gradient, as one of the basic firstorder operators in the supportoperators method. The discrete analog of the flux operator must also be the negative adjoint of the discrete divergence, in an inner product that is a discrete analog of the continuum inner product. The resulting method is conservative and the discrete analog of the variable coefficient Laplacian is symmetric and negative definite on nonuniform grids. In addition, on any grid, the discrete divergence is zero on constantmore »
 Authors:

 Los Alamos National Lab., NM (United States)
 Univ. of New Mexico, Albuquerque, NM (United States)
 Publication Date:
 OSTI Identifier:
 535127
 DOE Contract Number:
 W7405ENG36
 Resource Type:
 Journal Article
 Journal Name:
 Journal of Computational Physics
 Additional Journal Information:
 Journal Volume: 132; Journal Issue: 1; Other Information: PBD: 15 Mar 1997
 Country of Publication:
 United States
 Language:
 English
 Subject:
 66 PHYSICS; 42 ENGINEERING NOT INCLUDED IN OTHER CATEGORIES; 99 MATHEMATICS, COMPUTERS, INFORMATION SCIENCE, MANAGEMENT, LAW, MISCELLANEOUS; DIFFUSION; PARTIAL DIFFERENTIAL EQUATIONS; FINITE DIFFERENCE METHOD; ALGORITHMS
Citation Formats
Hyman, J, Shashkov, M, and Steinberg, S. The numerical solution of diffusion problems in strongly heterogeneous nonisotropic materials. United States: N. p., 1997.
Web. doi:10.1006/jcph.1996.5633.
Hyman, J, Shashkov, M, & Steinberg, S. The numerical solution of diffusion problems in strongly heterogeneous nonisotropic materials. United States. doi:10.1006/jcph.1996.5633.
Hyman, J, Shashkov, M, and Steinberg, S. Sat .
"The numerical solution of diffusion problems in strongly heterogeneous nonisotropic materials". United States. doi:10.1006/jcph.1996.5633.
@article{osti_535127,
title = {The numerical solution of diffusion problems in strongly heterogeneous nonisotropic materials},
author = {Hyman, J and Shashkov, M and Steinberg, S},
abstractNote = {A new secondorder finitedifference algorithm for the numerical solution of diffusion problems in strongly heterogeneous and nonisotropic media is constructed. On problems with rough coefficients or highly nonuniform grids, the new algorithm is superior to all other algorithms we have compared it with. For problems with smooth coefficients on smooth grids, the method is comparable with other secondorder methods. The new algorithm is formulated for logically rectangular grids and is derived using the supportoperators method. A key idea in deriving the method was to replace the usual inner product of vector functions by an inner product weighted by the inverse of the material properties tensor and to use the flux operator, defined as the material properties tensor times the gradient, rather than the gradient, as one of the basic firstorder operators in the supportoperators method. The discrete analog of the flux operator must also be the negative adjoint of the discrete divergence, in an inner product that is a discrete analog of the continuum inner product. The resulting method is conservative and the discrete analog of the variable coefficient Laplacian is symmetric and negative definite on nonuniform grids. In addition, on any grid, the discrete divergence is zero on constant vectors, the null space for the gradient is the constant functions, and, when the material properties are piecewise constant, the discrete flux operator is exact for piecewise linear functions. We compare the methods on some of the most difficult examples to be found in the literature. 21 refs., 12 figs., 5 tabs.},
doi = {10.1006/jcph.1996.5633},
journal = {Journal of Computational Physics},
number = 1,
volume = 132,
place = {United States},
year = {1997},
month = {3}
}