Lumping Techniques for Mixed Finite Element Diffusion Discretizations
- Lawrence Livermore National Lab. (LLNL), Livermore, CA (United States)
We are interested in positivity preserving spatial discretizations applicable to arbitrary order mixed finite element spatial discretizations of the radiative diffusion equations. Though not guaranteed to yield strict positivity, matrix lumping is appealing as a method of increasing solution positivity because unlike ad hoc fixups or strictly non-negative, non-linear solution representations, matrix lumping keeps the mixed finite element diffusion equations as a linear system of equations. Self-lumping schemes have previously been shown to be a generalization of more traditional matrix lumping techniques. Self-lumping schemes restrict integration quadrature points to the interpolation points used to define interpolatory finite element basis functions. As shown previously, the positivity and order of accuracy of self-lumping schemes are directly related to the choice of finite element interpolation point. Here, we develop self-lumping schemes for arbitrary order, mixed finite element spatial discretizations of the diffusion equation in Cartesian geometry. After introducing several possible self-lumping schemes appropriate for mixed finite elements, we test the ability of each scheme to increase solution positivity relative to an unlumped discretization using a simple test problem designed to induce negative solutions, evaluating efficacy as a function of finite element order and finite element interpolation point type. The method of manufactured solutions is then used to evaluate order of convergence for promising self-lumping schemes.
- Research Organization:
- Lawrence Livermore National Laboratory (LLNL), Livermore, CA (United States)
- Sponsoring Organization:
- USDOE National Nuclear Security Administration (NNSA)
- Grant/Contract Number:
- AC52-07NA27344
- OSTI ID:
- 1469467
- Report Number(s):
- LLNL-JRNL--744788; 898234
- Journal Information:
- Journal of Computational and Theoretical Transport, Journal Name: Journal of Computational and Theoretical Transport Journal Issue: na Vol. na; ISSN 2332-4309
- Publisher:
- Taylor and FrancisCopyright Statement
- Country of Publication:
- United States
- Language:
- English
Similar Records
Characterization of high order spatial discretizations and lumping techniques for discontinuous finite element SN transport
High-order solution methods for grey discrete ordinates thermal radiative transfer