Solving diffusion equations with rough coefficients in rough grids
- Los Alamos National Lab., NM (United States)
- Univ. of New Mexico, Albuquerque, NM (United States)
A finite-difference algorithm for the numerical solution of diffusion problems in strongly heterogeneous and nonisotropic media is constructed for logically rectangular grids. The performance of this algorithm is comparable to other algorithms for problems with smooth coefficients and regular grids, and it is superior for problems with rough coefficients and/or skewed grids. The algorithm is derived using the support-operators method, which constructs discrete analogs of the divergence and flux operator that satisfy discrete analogs of the important integral identities relating the continuum operators. This paper gives the first application of this method to the solution of diffusion problems in heterogeneous an nonisotropic media. The support-operators method forces the discrete analog of the flux operator to be the negative adjoint of the discrete divergence in an inner product weighted by the conductivity, as in the differential case. Once this is accomplished, many other important properties follow; for example, the scheme is conservative and the discrete analog of the variable material Laplacian is symmetric and negative definite. In addition, on any grid, the discrete divergence is zero on constant vectors and the discrete divergence is zero on constant vectors and the discrete divergence is zero on constant vectors and the discrete flux operator is exact for linear functions in case when K is piecewise constant. Moreover, the discrete gradient`s null space is the constant functions, just as in the continuum. Because the algorithm is flux based, it has twice as many unknowns as more standard algorithms. However, the matrices that need to be inverted are symmetric and positive definite, so the most powerful linear solvers can be applied. Also, the scheme is second-order accurate so, all things considered, it is efficient.
- DOE Contract Number:
- W-7405-ENG-36
- OSTI ID:
- 530648
- Journal Information:
- Journal of Computational Physics, Vol. 129, Issue 2; Other Information: PBD: Dec 1996
- Country of Publication:
- United States
- Language:
- English
Similar Records
Support Operators Method for the Diffusion Equation in Multiple Materials
Mimetic difference approximations of partial differential equations