A structure-preserving staggered semi-implicit four-split finite volume scheme for continuum mechanics on unstructured meshes
Abstract
We present a new semi-implicit structure-preserving (SP) finite volume discretization for a unified first-order hyperbolic model of continuum mechanics. This framework provides a common mathematical description for fluids and solids and incorporates transport, viscous effects, heat conduction, and elastic deformations within a single system of hyperbolic partial differential equations. The coexistence of several physical mechanisms gives rise to multiple characteristic wave speeds, resulting in severe time-step restrictions for fully explicit discretizations. To overcome this difficulty, we develop a four-split scheme in which the governing equations are decomposed into convective, temperature, mechanical, and pressure subsystems. The convective subsystem is the only one that is advanced explicitly in time, while the remaining subsystems are treated implicitly in a sequential manner. As a consequence, the time step restriction of the proposed method depends only on the material velocity and is independent of the acoustic, shear, and thermal wave speeds that characterize the model. The scheme is designed to preserve key structural properties of the underlying equations on unstructured grids. In particular, a compatible vertex-staggered spatial discretization on triangles allows the curl-free involutions associated with the distortion field and thermal impulse to be respected whenever the relaxation source terms in the governing PDE are linear or absent. At the same time, the scheme is asymptotic preserving as it is consistent with the low Mach number limit and with the stiff relaxation limits of the GPR model that recover the classical Navier-Stokes-Fourier equations in fluid mechanics. A series of numerical experiments covering both fluids and solids demonstrates the accuracy, robustness, and multi-scale capabilities of the proposed approach.
Related
- A Structure-Preserving and Thermodynamically Compatible Cell-Centered Lagrangian Finite Volume Scheme for Continuum Mechanics
- Virtual Convolved Action for Continuum Solid Mechanics and Electromagnetics
- A Discrete-to-Continuum Variational Theory for Area-Preserving Surface Parameterization
- Geometrization of continuum mechanics and nonequilibrium thermodynamics: Heat conduction
- Uniqueness and structure-preserving numerical simulations for an Ultrahyperbolic Schrödinger equation
BibTeX
@misc{dumbser2026structurepreserving,
title = {{A structure-preserving staggered semi-implicit four-split finite volume scheme for continuum mechanics on unstructured meshes}},
author = {M. Dumbser and W. Boscheri and M. Tavelli and Andrea Thomann},
howpublished = {HAL},
year = {2026},
month = sep,
doi = {10.48550/arxiv.2609.25935},
eprint = {2609.25935},
archivePrefix = {arXiv},
url = {https://hal.science/hal-05835901v1},
note = {CC BY-NC-SA 4.0},
}
RIS
TY - UNPB
TI - A structure-preserving staggered semi-implicit four-split finite volume scheme for continuum mechanics on unstructured meshes
AU - M. Dumbser
AU - W. Boscheri
AU - M. Tavelli
AU - Andrea Thomann
PY - 2026
DA - 2026/09/22
PB - HAL
DO - 10.48550/arxiv.2609.25935
UR - https://hal.science/hal-05835901v1
AB - We present a new semi-implicit structure-preserving (SP) finite volume discretization for a unified first-order hyperbolic model of continuum mechanics. This framework provides a common mathematical description for fluids and solids and incorporates transport, viscous effects, heat conduction, and elastic deformations within a single system of hyperbolic partial differential equations. The coexistence of several physical mechanisms gives rise to multiple characteristic wave speeds, resulting in severe time-step restrictions for fully explicit discretizations. To overcome this difficulty, we develop a four-split scheme in which the governing equations are decomposed into convective, temperature, mechanical, and pressure subsystems. The convective subsystem is the only one that is advanced explicitly in time, while the remaining subsystems are treated implicitly in a sequential manner. As a consequence, the time step restriction of the proposed method depends only on the material velocity and is independent of the acoustic, shear, and thermal wave speeds that characterize the model. The scheme is designed to preserve key structural properties of the underlying equations on unstructured grids. In particular, a compatible vertex-staggered spatial discretization on triangles allows the curl-free involutions associated with the distortion field and thermal impulse to be respected whenever the relaxation source terms in the governing PDE are linear or absent. At the same time, the scheme is asymptotic preserving as it is consistent with the low Mach number limit and with the stiff relaxation limits of the GPR model that recover the classical Navier-Stokes-Fourier equations in fluid mechanics. A series of numerical experiments covering both fluids and solids demonstrates the accuracy, robustness, and multi-scale capabilities of the proposed approach.
ER -