Sharp-interface Simulations of Energetic Multiphase Flows with Large Density and Viscosity Ratios
Abstract
Flows with high density ratios, such as wave breaking and air entrainment in maritime applications, remain challenging to simulate due to their energetic and strongly nonlinear nature. In such regimes, maintaining numerical robustness is difficult when using the commonly adopted velocity-based formulation. The Consistent Mass-Momentum (CMOM) transport framework improves numerical robustness by enforcing fundamental physical properties, most notably momentum conservation and semi-discrete energy-conserving. However, CMOM replaces the advection of a continuous velocity field with that of a discontinuous momentum field. When combined with sharp interface methods, this leads to severe momentum shocks, for which conventional shock-capturing schemes are ineffective. To reconcile physical fidelity with numerical robustness, this work proposes a Synchronized Donor-Region of Momentum fluxes (SynDRoM) that enforces monotonicity of the transported velocity field. The resulting algorithm effectively eliminates spurious velocity oscillations without sacrificing physical fidelity, as demonstrated through scalar transport and interfacial shear instability test cases. Beyond difficulties from large density ratio, improper estimation of viscosity in the vicinity of the interface can introduce numerical instabilities at finite time steps, thereby undermining overall robustness. To address this issue, a viscosity limiter based on the bounded kinetic viscosity concept is introduced and validated using a gravity-driven plane shear flow. Finally, a breaking wave simulation is performed to assess the combined performance of the proposed physics-preserving numerical schemes for multiphase flows.
Related
- Neural-Network-based Viscosity Closure for Non-Newtonian Multiphase Flows
- Constraining Effective Viscosity in the Intracluster Medium via the Thermal Sunyaev-Zeldovich Effect -- Predictions from the SLOW Constrained Coma Cluster Simulations
- A Robust All-Mach Six-Equation Diffuse-Interface Method for Multiphase Flows with Surface Tension
- Variational formulation of the constitutive relaxation at finite rate of the mechanical state in multiphase flows with the diffuse-interface method
- Energetic particles, Shafranov shift and finite β effects on TAE, KBM and ITG instabilities in global electromagnetic gyrokinetic simulations
BibTeX
@misc{huang2026sharpinterface,
title = {{Sharp-interface Simulations of Energetic Multiphase Flows with Large Density and Viscosity Ratios}},
author = {Tzu-Yao Huang and Nicolas Valle and Artur K. Lidtke and Kelli Hendrickson and Gabriel D. Weymouth},
howpublished = {arXiv},
year = {2026},
month = oct,
doi = {10.48550/arxiv.2606.02467},
eprint = {2606.02467},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2606.02467},
note = {CC BY 4.0},
}
RIS
TY - UNPB
TI - Sharp-interface Simulations of Energetic Multiphase Flows with Large Density and Viscosity Ratios
AU - Tzu-Yao Huang
AU - Nicolas Valle
AU - Artur K. Lidtke
AU - Kelli Hendrickson
AU - Gabriel D. Weymouth
PY - 2026
DA - 2026/10/02
PB - arXiv
DO - 10.48550/arxiv.2606.02467
UR - https://arxiv.org/abs/2606.02467
AB - Flows with high density ratios, such as wave breaking and air entrainment in maritime applications, remain challenging to simulate due to their energetic and strongly nonlinear nature. In such regimes, maintaining numerical robustness is difficult when using the commonly adopted velocity-based formulation. The Consistent Mass-Momentum (CMOM) transport framework improves numerical robustness by enforcing fundamental physical properties, most notably momentum conservation and semi-discrete energy-conserving. However, CMOM replaces the advection of a continuous velocity field with that of a discontinuous momentum field. When combined with sharp interface methods, this leads to severe momentum shocks, for which conventional shock-capturing schemes are ineffective. To reconcile physical fidelity with numerical robustness, this work proposes a Synchronized Donor-Region of Momentum fluxes (SynDRoM) that enforces monotonicity of the transported velocity field. The resulting algorithm effectively eliminates spurious velocity oscillations without sacrificing physical fidelity, as demonstrated through scalar transport and interfacial shear instability test cases. Beyond difficulties from large density ratio, improper estimation of viscosity in the vicinity of the interface can introduce numerical instabilities at finite time steps, thereby undermining overall robustness. To address this issue, a viscosity limiter based on the bounded kinetic viscosity concept is introduced and validated using a gravity-driven plane shear flow. Finally, a breaking wave simulation is performed to assess the combined performance of the proposed physics-preserving numerical schemes for multiphase flows.
ER -