Monotonicity-preserving finite element schemes based on differentiable nonlinear stabilization
Visualitza/Obre
Cita com:
hdl:2117/100080
Tipus de documentArticle
Data publicació2017-01
Condicions d'accésAccés obert
Llevat que s'hi indiqui el contrari, els
continguts d'aquesta obra estan subjectes a la llicència de Creative Commons
:
Reconeixement-NoComercial-SenseObraDerivada 3.0 Espanya
ProjecteNUMEXAS - NUMERICAL METHODS AND TOOLS FOR KEY EXASCALE COMPUTING CHALLENGES IN ENGINEERING AND APPLIED SCIENCES (EC-FP7-611636)
COMFUS - Computational Methods for Fusion Technology (EC-FP7-258443)
COMFUS - Computational Methods for Fusion Technology (EC-FP7-258443)
Abstract
In this work, we propose a nonlinear stabilization technique for scalar conservation laws with implicit time stepping. The method relies on an artificial diffusion method, based on a graph-Laplacian operator. It is nonlinear, since it depends on a shock detector. Further, the resulting method is linearity preserving. The same shock detector is used to gradually lump the mass matrix. The resulting method is LED, positivity preserving, and also satisfies a global DMP. Lipschitz continuity has also been proved. However, the resulting scheme is highly nonlinear, leading to very poor nonlinear convergence rates. We propose a smooth version of the scheme, which leads to twice differentiable nonlinear stabilization schemes. It allows one to straightforwardly use Newton’s method and obtain quadratic convergence. In the numerical experiments, steady and transient linear transport, and transient Burgers’ equation have been considered in 2D. Using the Newton method with a smooth version of the scheme we can reduce 10 to 20 times the number of iterations of Anderson acceleration with the original non-smooth scheme. In any case, these properties are only true for the converged solution, but not for iterates. In this sense, we have also proposed the concept of projected nonlinear solvers, where a projection step is performed at the end of every nonlinear iterations onto a FE space of admissible solutions. The space of admissible solutions is the one that satisfies the desired monotonic properties (maximum principle or positivity).
CitacióBadia, S., Bonilla, J. Monotonicity-preserving finite element schemes based on differentiable nonlinear stabilization. "Computer methods in applied mechanics and engineering", Gener 2017, vol. 313, p. 133-158.
ISSN0045-7825
Versió de l'editorhttp://www.sciencedirect.com/science/article/pii/S0045782516306405
Fitxers | Descripció | Mida | Format | Visualitza |
---|---|---|---|---|
art018-lssc.pdf | 1,000Mb | Visualitza/Obre |