skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: The numerical solution of diffusion problems in strongly heterogeneous non-isotropic materials

Abstract

A new second-order finite-difference algorithm for the numerical solution of diffusion problems in strongly heterogeneous and non-isotropic 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 second-order methods. The new algorithm is formulated for logically rectangular grids and is derived using the support-operators 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 first-order operators in the support-operators 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 » 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.« less

Authors:
;  [1];  [2]
  1. Los Alamos National Lab., NM (United States)
  2. Univ. of New Mexico, Albuquerque, NM (United States)
Publication Date:
OSTI Identifier:
535127
DOE Contract Number:  
W-7405-ENG-36
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 non-isotropic 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 non-isotropic 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 non-isotropic materials". United States. doi:10.1006/jcph.1996.5633.
@article{osti_535127,
title = {The numerical solution of diffusion problems in strongly heterogeneous non-isotropic materials},
author = {Hyman, J and Shashkov, M and Steinberg, S},
abstractNote = {A new second-order finite-difference algorithm for the numerical solution of diffusion problems in strongly heterogeneous and non-isotropic 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 second-order methods. The new algorithm is formulated for logically rectangular grids and is derived using the support-operators 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 first-order operators in the support-operators 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}
}