Research reports

SECOND-ORDER OPTIMALLY STABLE IMEX (PSEUDO-)STAGGERED GALERKIN DISCRETIZATION: APPLICATION TO LAVA FLOW MODELING

by F. Gatti and G. Orlando

(Report number 2025-26)

Abstract
We present second-order optimally stable Implicit-Explicit (IMEX) Runge-Kutta (RK) schemes with application to a modified set of shallow water equations that can be used to model the dynamics of lava flows. The schemes are optimally stable in the sense that they satisfy, at the space-time discretization level, a condition analogous to the L-stability of Runge-Kutta methods for ordinary differential equations. A novel (pseudo-)staggered Galerkin scheme is introduced, which can be interpreted as an extension of the classical two-step Taylor-Galerkin (TG2) scheme. The method is derived by combining a von Neumann stability analysis with a Lax-Wendroff procedure. For the discretization of the non-conservative terms that characterize the lava flow model, we employ the Path-Conservative (PC) method. The proposed scheme is evaluated on a number of relevant test cases, demonstrating accuracy, robustness, and well-balancing properties for the lava flow model.

Keywords: Shallow water equations, Well-balanced schemes, IMEX schemes, von Neumann analysis, Lax-Wendroff procedure, (pseudo-)staggered methods

BibTeX
@Techreport{GO25_1147,
  author = {F. Gatti and G. Orlando},
  title = {SECOND-ORDER OPTIMALLY STABLE IMEX (PSEUDO-)STAGGERED GALERKIN DISCRETIZATION: APPLICATION TO LAVA FLOW MODELING},
  institution = {Seminar for Applied Mathematics, ETH Z{\"u}rich},
  number = {2025-26},
  address = {Switzerland},
  url = {https://www.sam.math.ethz.ch/sam_reports/reports_final/reports2025/2025-26.pdf },
  year = {2025}
}

Disclaimer
© Copyright for documents on this server remains with the authors. Copies of these documents made by electronic or mechanical means including information storage and retrieval systems, may only be employed for personal use. The administrators respectfully request that authors inform them when any paper is published to avoid copyright infringement. Note that unauthorised copying of copyright material is illegal and may lead to prosecution. Neither the administrators nor the Seminar for Applied Mathematics (SAM) accept any liability in this respect. The most recent version of a SAM report may differ in formatting and style from published journal version. Do reference the published version if possible (see SAM Publications).

JavaScript has been disabled in your browser