Finite Element Approximation of the Diffusion Operator on Tetrahedra

Authors:
Mario Putti;Christian Cordes
Affiliations:
-;-
Venue:
SIAM Journal on Scientific Computing
Year:
1998

Citing 0
Cited 7

Unstructured grid optimization for improved monotonicity of discrete solutions of elliptic equations with highly anisotropic coefficients

Journal of Computational Physics
Mesh locking effects in the finite volume solution of 2-D anisotropic diffusion equations

Journal of Computational Physics
Coupling geological and numerical models to simulate groundwater flow and contaminant transport in fractured media

Computers & Geosciences
Comparison of 3D flow fields arising in mixed and standard unstructured finite elements

ICCS'03 Proceedings of the 1st international conference on Computational science: PartI
Guaranteed and Fully Robust a posteriori Error Estimates for Conforming Discretizations of Diffusion Problems with Discontinuous Coefficients

Journal of Scientific Computing
The MAST-edge centred lumped scheme for the flow simulation in variably saturated heterogeneous porous media

Journal of Computational Physics
Monotonic solution of heterogeneous anisotropic diffusion problems

Journal of Computational Physics

Quantified Score

Hi-index	0.02

Visualization

Abstract

Linear Galerkin finite element discretizations of the Laplace operator produce nonpositive stiffness coefficients for internal element edges of two-dimensional Delaunay triangulations. This property, also called the positive transmissibility (PT) condition, is a prerequisite for the existence of an M-matrix and ensures that nonphysical local extrema are not present in the solution. For tetrahedral elements, it has already been shown that the linear Galerkin approach does not in general satisfy the PT condition. We propose a modification of the three-dimensional Galerkin scheme that, if a Delaunay triangulation is used, satisfies the PT condition for internal edges and, if further conditions on the boundary are specified, yields an $M$-matrix. The proposed approach can also be extended to the general diffusion operator with nonconstant or anisotropic coefficients.