Flux vector splittings and semi-implicit finite volume schemes for an isentropic two-phase flow model
Abstract
We present two new semi-implicit finite volume schemes for the 1D hyperbolic and thermodynamically compatible (HTC) isentropic two-phase flow model of Romenski et al., written in terms of mixture density, volume fraction, partial gas density, mixture momentum and the relative phase velocity. In gas-liquid mixtures, the sound speed of the weakly compressible liquid is much larger than the gas sound speed and flow velocity, so the time step of fully explicit discretizations is dictated by liquid acoustic waves, even when the interest lies in the dynamics of the gas. To overcome this difficulty, the model's flux is split at the continuous level into an explicit subsystem (gas dynamics) and an implicit subsystem (liquid acoustic waves), leading to a nonlinear system for liquid pressure solved with the Newton method. In the first splitting the entire liquid phase is treated implicitly, while in the second the liquid convective terms are moved to the explicit subsystem, so only liquid acoustic waves are treated implicitly and nested Picard iterations are avoided. Thus, the time step of both schemes is restricted only by flow velocities and gas sound speed, independent of liquid sound speed. We prove and numerically confirm that both schemes exactly preserve the Abgrall condition (uniform velocity and pressure across material interfaces). For well-prepared initial data, both schemes are consistent with the volume balance of the liquid's low Mach number limit. Riemann problems, pressure equilibrium tests and acoustic tests show that in non-stiff regimes the semi-implicit schemes reproduce the results of a fully explicit scheme, while in the stiff regime they allow time steps that are up to 3-4 orders of magnitude larger, avoiding the excessive numerical dissipation of explicit schemes whose numerical viscosity is dictated by liquid sound speed.
Related
- Well-balanced semi-implicit relaxed finite volume schemes for hyperbolic balance laws
- Are implicit schemes needed in adaptive time stepping for PDEs?
- Neural Finite Volume and Mixture-of-Experts Physics-Informed Networks for the Isentropic Cargo–LeRoux Model: An Exact-Solution-Grounded Comparative Study
- GPU-Accelerated Nested Schur-Complement Solvers for Atmospheric Flow Simulations Using Horizontally-Explicit, Vertically-Implicit Schemes and Discontinuous Galerkin Methods
- Tangent bundles of Quot schemes of vector bundles on curves
BibTeX
@misc{dumbser2026flux,
title = {{Flux vector splittings and semi-implicit finite volume schemes for an isentropic two-phase flow model}},
author = {Michael Dumbser and Gabriella Puppo and Alessandra Spilimbergo and Tommaso Tenna},
howpublished = {arXiv},
year = {2026},
month = nov,
doi = {10.48550/arxiv.2611.00781},
eprint = {2611.00781},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2611.00781},
note = {CC BY-NC-SA 4.0},
}
RIS
TY - UNPB
TI - Flux vector splittings and semi-implicit finite volume schemes for an isentropic two-phase flow model
AU - Michael Dumbser
AU - Gabriella Puppo
AU - Alessandra Spilimbergo
AU - Tommaso Tenna
PY - 2026
DA - 2026/11/01
PB - arXiv
DO - 10.48550/arxiv.2611.00781
UR - https://arxiv.org/abs/2611.00781
AB - We present two new semi-implicit finite volume schemes for the 1D hyperbolic and thermodynamically compatible (HTC) isentropic two-phase flow model of Romenski et al., written in terms of mixture density, volume fraction, partial gas density, mixture momentum and the relative phase velocity. In gas-liquid mixtures, the sound speed of the weakly compressible liquid is much larger than the gas sound speed and flow velocity, so the time step of fully explicit discretizations is dictated by liquid acoustic waves, even when the interest lies in the dynamics of the gas. To overcome this difficulty, the model's flux is split at the continuous level into an explicit subsystem (gas dynamics) and an implicit subsystem (liquid acoustic waves), leading to a nonlinear system for liquid pressure solved with the Newton method. In the first splitting the entire liquid phase is treated implicitly, while in the second the liquid convective terms are moved to the explicit subsystem, so only liquid acoustic waves are treated implicitly and nested Picard iterations are avoided. Thus, the time step of both schemes is restricted only by flow velocities and gas sound speed, independent of liquid sound speed. We prove and numerically confirm that both schemes exactly preserve the Abgrall condition (uniform velocity and pressure across material interfaces). For well-prepared initial data, both schemes are consistent with the volume balance of the liquid's low Mach number limit. Riemann problems, pressure equilibrium tests and acoustic tests show that in non-stiff regimes the semi-implicit schemes reproduce the results of a fully explicit scheme, while in the stiff regime they allow time steps that are up to 3-4 orders of magnitude larger, avoiding the excessive numerical dissipation of explicit schemes whose numerical viscosity is dictated by liquid sound speed.
ER -