Flux-corrected transport algorithms preserving the eigenvalue range of symmetric tensor quantities

JOURNAL OF COMPUTATIONAL PHYSICS

卷:350

Article

Lohmann, Christoph¹ 

[1] TU Dortmund Univ, Inst Appl Math LS 3, Vogelpothsweg 87, D-44227 Dortmund, Germany

关键词: Tensor quantity;    Continuous Galerkin method;    Flux-corrected transport;    Artificial diffusion;    Local discrete maximum principles;   

DOI  :  10.1016/j.jcp.2017.09.009

来源: Elsevier

【 摘 要 】

This paper presents a new approach to constraining the eigenvalue range of symmetric tensors in numerical advection schemes based on the flux-corrected transport (FCT) algorithm and a continuous finite element discretization. In the context of element-based FEM-FCT schemes for scalar conservation laws, the numerical solution is evolved using local extremum diminishing (LED) antidiffusive corrections of a low order approximation which is assumed to satisfy the relevant inequality constraints. The application of a limiter to antidiffusive element contributions guarantees that the corrected solution remains bounded by the local maxima and minima of the low order predictor. The FCT algorithm to be presented in this paper guarantees the LED property for the maximal and minimal eigenvalues of the transported tensor at the low order evolution step. At the antidiffusive correction step, this property is preserved by limiting the antidiffusive element contributions to all components of the tensor in a synchronized manner. The definition of the element-based correction factors for FCT is based on perturbation bounds for auxiliary tensors which are constrained to be positive semidefinite to enforce the generalized LED condition. The derivation of sharp bounds involves calculating the roots of polynomials of degree up to 3. As inexpensive and numerically stable alternatives, limiting techniques based on appropriate estimates are considered. The ability of the new limiters to enforce local bounds for the eigenvalue range is confirmed by numerical results for 2D advection problems. (C) 2017 Elsevier Inc. All rights reserved.

【 授权许可】

Free