Title | ||
---|---|---|
A constrained finite element method satisfying the discrete maximum principle for anisotropic diffusion problems |
Abstract | ||
---|---|---|
Nonlinear constrained finite element approximations to anisotropic diffusion problems are considered. Starting with a standard (linear or bilinear) Galerkin discretization, the entries of the stiffness matrix are adjusted so as to enforce sufficient conditions of the discrete maximum principle (DMP). An algebraic splitting is employed to separate the contributions of negative and positive off-diagonal coefficients which are associated with diffusive and antidiffusive numerical fluxes, respectively. In order to prevent the formation of spurious undershoots and overshoots, a symmetric slope limiter is designed for the antidiffusive part. The corresponding upper and lower bounds are defined using an estimate of the steepest gradient in terms of the maximum and minimum solution values at surrounding nodes. The recovery of nodal gradients is performed by means of a lumped-mass L"2 projection. The proposed slope limiting strategy preserves the consistency of the underlying discrete problem and the structure of the stiffness matrix (symmetry, zero row and column sums). A positivity-preserving defect correction scheme is devised for the nonlinear algebraic system to be solved. Numerical results and a grid convergence study are presented for a number of anisotropic diffusion problems in two space dimensions. |
Year | DOI | Venue |
---|---|---|
2009 | 10.1016/j.jcp.2009.01.031 | J. Comput. Physics |
Keywords | Field | DocType |
finite element method,35j25,antidiffusive numerical flux,diffusion problem,gradient recovery,proposed slope,nonnegativity constraints,anisotropic diffusion,antidiffusive part,anisotropic diusion,discrete maximum principle,35b50,stiffness matrix,nonlinear algebraic system,algebraic splitting,numerical result,slope limiting,65n30,anisotropic diffusion problem,65n22,upper and lower bounds,satisfiability | Anisotropic diffusion,Discretization,Mathematical optimization,Maximum principle,Nonlinear system,Mathematical analysis,Galerkin method,Finite element method,Stiffness matrix,Flux limiter,Mathematics | Journal |
Volume | Issue | ISSN |
228 | 9 | Journal of Computational Physics |
Citations | PageRank | References |
18 | 0.98 | 4 |
Authors | ||
3 |
Name | Order | Citations | PageRank |
---|---|---|---|
Dmitri Kuzmin | 1 | 167 | 23.90 |
M. J. Shashkov | 2 | 33 | 3.57 |
Daniil Svyatskiy | 3 | 170 | 12.13 |