Generalization of Strang-Fix remark for discrete maximum principles on tetrahedral meshes

1 Introduction

• Maximum principles are qualitative properties of elliptic equations that numerical approximations are expected to preserve through discrete analogues.The cited discussion presents them as properties that adequately imitate the physical behavior of modeled phenomena.Source: §1 · p. 2
• Prior work established angle-based DMP guarantees for simplicial meshes, including nonobtuse tetrahedral meshes and higher-dimensional settings.The cited literature also considers nonlinear elliptic models and mesh generation and refinement.Source: §1 · p. 2
• Nonobtuse dihedral angles guarantee nonpositive contributions to off-diagonal finite-element entries, supporting the usual M-matrix basis for DMPs.The standard entry calculation uses facet areas and element volumes.Source: §1 · p. 2
• In two dimensions, adjacent-angle cotangents express off-diagonal entries, and Strang and Fix’s trigonometric observation yields a weaker sign condition than elementwise nonobtuseness.That condition remains within the M-matrix framework and is closely related to Delaunay triangulations.Source: §1 · pp. 2–3
• The paper generalizes the edge-length and dihedral-angle-cotangent representation of off-diagonal stiffness entries to tetrahedral meshes and gives a new proof of that representation.This representation supports proofs of DMPs for some structured tetrahedral meshes containing obtuse dihedral angles.Source: §1 · p. 3

2 Model problem and maximum principles

• The model is a Poisson test problem on a bounded polytopic domain with Lipschitz boundary and sufficiently smooth forcing.The classical maximum principle supplies the continuous sign-preservation property used to motivate its discrete analogue.Source: §2 · pp. 4–5
• The section uses continuous piecewise-linear basis functions on a simplicial mesh, with interior-node functions spanning the finite-dimensional approximation space.The finite-element approximation is defined in that space.Source: §2 · pp. 4–5
• For the model problem, nonpositivity of the forcing implies nonpositive load-vector entries, and a monotone stiffness matrix then yields nonpositive finite-element coefficients.The basis functions are nonnegative, so nonpositive coefficients give a nonpositive approximation.Source: §2 · p. 5
• The stiffness matrix is symmetric positive definite, so the DMP criterion reduces to ensuring that its off-diagonal entries are nonpositive.Symmetric positive-definite matrices with nonpositive off-diagonal entries are Stieltjes matrices, a class of monotone matrices.Source: §2 · p. 5
• The section notes that nonzero Dirichlet and mixed boundary conditions can be handled with small technical modifications.It refers to prior works for details.Source: §2 · p. 5

• The section introduces tetrahedral geometric identities and Cayley–Menger determinants and cofactors before presenting the local-entry result.Source: §3 · pp. 6–7
• Theorem 2 gives a geometric formula for each tetrahedron’s contribution to an off-diagonal stiffness entry.Source: §3 · pp. 5–7
• Corollary 1 obtains a global entry by summing local contributions over tetrahedra sharing its edge.Source: §3 · p. 8
• The representation can avoid computing mesh-element face areas and volumes, an implementation benefit.Source: §3 · p. 8
• When corresponding edge lengths are equal, negative contributions can offset positive contributions arising from obtuse dihedral angles.Source: §3 · p. 8
• Allowing slightly obtuse dihedral angles is practically attractive because nonobtuseness is hard to preserve during global and especially local mesh refinements.Source: §3 · p. 8

4 Examples and final comments

• In Example 1, four tetrahedra surround B1B2; symmetry makes two dihedral angles sum to π, canceling their contributions to a12.Moving B1 and B2 along opposite arcs can make the other two relevant angles nonobtuse, so their contributions are nonpositive.Source: §4 · pp. 8–9
• Example 1 permits a dihedral angle arbitrarily close to π while retaining a12 ≤ 0.The vertices B3–B6 form a rectangle, and B1 and B2 are mirror images on a circle perpendicular to its plane.Source: §4 · p. 9
• In Example 2, one dihedral angle is obtuse, but the angles along edges ℓ56 and ℓ34 sum to less than π, while those along ℓ46 and ℓ35 are acute.Together, these angle relations make a12 strictly negative.Source: §4 · pp. 9–10
• Clusters around every interior edge as in the examples yield a valid DMP for the model problem despite some obtuse dihedral angles.The conclusion is scoped to meshes with this cluster structure.Source: §4 · p. 10