Anchored HLLD: A Riemann solver for a hypoelastic diffuse-interface model of compressible fluid-solid interaction
Abstract
We present an HLLD Riemann solver for a multicomponent hypoelastic diffuse-interface model of compressible fluid-solid interaction, with no separately computed interface and a deviatoric stress advanced by an objective hypoelastic rate law outside conservation-law form. The difficulty specific to this model is that the stress evolution is nonconservative and, unlike the familiar volume-fraction terms of Kapila-type models, cannot be left outside the Riemann solve once the shear waves are resolved. The shear coefficient ρ(G+τxx) sets the impedance of the resolved shear waves and so belongs to their eigenstructure. Anchored HLLD carries it inside the solve by evaluating the HLLD fan twice at each face, anchored in turn to each adjacent cell-centered state, with every cell advancing from the flux anchored to its own state. HLL, HLLC, and anchored HLLD form a three-tier hierarchy of increasing wave resolution in one framework. We prove local structural properties: consistency at the reconstruction order, exact species-mass conservation, and tangential-momentum and energy drift vanishing under refinement at the reconstruction's face-jump order on smooth fields, with one-dimensional tests illustrating these structural properties and demonstrating first-order drift convergence on discontinuous fields. Across one-dimensional tests and two-dimensional production simulations in the open-source solver MFC, anchored HLLD resolves both the material contact and the shear waves, HLLC the contact alone, and HLL neither sharply. Where refinement studies are shown, anchored HLLD converges toward exact or semi-analytic solutions. A three-dimensional shock-ellipsoid interaction demonstrates the solver with all six stress components active.
Related
- A conservative coupling method of sharp-interface and multi-species model for compressible reacting gas-liquid flows with phase change
- A Momentum-weighted Interface Stabilization Approach for Fluid-Structure Interaction with Strong Added-Mass Effects
- A Finite Element Model for Fluid--Poroelastic Structure Interaction with Dissipative Interface
- Anchored Exploration in Product Evolution: A Phylogenetic Network Framework for Semiconductor Designs
- Probing Local Charm Hadronization Dynamics with Hadron-Anchored Energy Correlators
BibTeX
@misc{zhang2026anchored,
title = {{Anchored HLLD: A Riemann solver for a hypoelastic diffuse-interface model of compressible fluid-solid interaction}},
author = {Yanjun Zhang and Tim Colonius},
howpublished = {arXiv},
year = {2026},
month = oct,
doi = {10.48550/arxiv.2610.31958},
eprint = {2610.31958},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2610.31958},
}
RIS
TY - UNPB
TI - Anchored HLLD: A Riemann solver for a hypoelastic diffuse-interface model of compressible fluid-solid interaction
AU - Yanjun Zhang
AU - Tim Colonius
PY - 2026
DA - 2026/10/22
PB - arXiv
DO - 10.48550/arxiv.2610.31958
UR - https://arxiv.org/abs/2610.31958
AB - We present an HLLD Riemann solver for a multicomponent hypoelastic diffuse-interface model of compressible fluid-solid interaction, with no separately computed interface and a deviatoric stress advanced by an objective hypoelastic rate law outside conservation-law form. The difficulty specific to this model is that the stress evolution is nonconservative and, unlike the familiar volume-fraction terms of Kapila-type models, cannot be left outside the Riemann solve once the shear waves are resolved. The shear coefficient ρ(G+τ_xx) sets the impedance of the resolved shear waves and so belongs to their eigenstructure. Anchored HLLD carries it inside the solve by evaluating the HLLD fan twice at each face, anchored in turn to each adjacent cell-centered state, with every cell advancing from the flux anchored to its own state. HLL, HLLC, and anchored HLLD form a three-tier hierarchy of increasing wave resolution in one framework. We prove local structural properties: consistency at the reconstruction order, exact species-mass conservation, and tangential-momentum and energy drift vanishing under refinement at the reconstruction's face-jump order on smooth fields, with one-dimensional tests illustrating these structural properties and demonstrating first-order drift convergence on discontinuous fields. Across one-dimensional tests and two-dimensional production simulations in the open-source solver MFC, anchored HLLD resolves both the material contact and the shear waves, HLLC the contact alone, and HLL neither sharply. Where refinement studies are shown, anchored HLLD converges toward exact or semi-analytic solutions. A three-dimensional shock-ellipsoid interaction demonstrates the solver with all six stress components active.
ER -