Many interesting applications of hyperbolic systems of equations are stiff, and require the time step to satisfy restrictive stability conditions. One way to avoid small time steps is to use implicit time integration. Implicit integration is quite straightforward for first-order schemes. High order schemes instead also need to control spurious oscillations, which requires limiting in space and time also in the linear case. We propose a framework to simplify considerably the application of high order non-oscillatory schemes through the introduction of a low order implicit predictor, which is used both to set up the nonlinear weights of a standard high order space reconstruction, and to achieve limiting in time. In this preliminary work, we concentrate on the case of a third-order scheme, based on diagonally implicit Runge Kutta (DIRK) integration in time and central weighted essentially non-oscillatory (CWENO) reconstruction in space. The numerical tests involve linear and nonlinear scalar conservation laws.
Quinpi: Integrating Conservation Laws with CWENO Implicit Methods / Puppo, Gabriella; Semplice, Matteo; Visconti, Giuseppe. - In: COMMUNICATIONS ON APPLIED MATHEMATICS AND COMPUTATION. - ISSN 2096-6385. - (2023). [10.1007/s42967-021-00171-0]
Quinpi: Integrating Conservation Laws with CWENO Implicit Methods
Gabriella Puppo
;Matteo Semplice;Giuseppe Visconti
2023
Abstract
Many interesting applications of hyperbolic systems of equations are stiff, and require the time step to satisfy restrictive stability conditions. One way to avoid small time steps is to use implicit time integration. Implicit integration is quite straightforward for first-order schemes. High order schemes instead also need to control spurious oscillations, which requires limiting in space and time also in the linear case. We propose a framework to simplify considerably the application of high order non-oscillatory schemes through the introduction of a low order implicit predictor, which is used both to set up the nonlinear weights of a standard high order space reconstruction, and to achieve limiting in time. In this preliminary work, we concentrate on the case of a third-order scheme, based on diagonally implicit Runge Kutta (DIRK) integration in time and central weighted essentially non-oscillatory (CWENO) reconstruction in space. The numerical tests involve linear and nonlinear scalar conservation laws.File | Dimensione | Formato | |
---|---|---|---|
Puppo_Quinpi_ Integrating_2022.pdf
accesso aperto
Tipologia:
Versione editoriale (versione pubblicata con il layout dell'editore)
Licenza:
Tutti i diritti riservati (All rights reserved)
Dimensione
3.44 MB
Formato
Adobe PDF
|
3.44 MB | Adobe PDF |
I documenti in IRIS sono protetti da copyright e tutti i diritti sono riservati, salvo diversa indicazione.