Wave or Physics-Appropriate Multidimensional Upwinding Approach for Compressible Multiphase Flows
arXiv:2501.02416 · doi:10.1016/j.jcp.2025.114157
Abstract
This paper introduces multidimensional algorithms for simulating multiphase flows, leveraging the wave structure of the Euler equations in characteristic space and the physical properties of variables in physical space. The algorithm applies different reconstruction schemes to acoustic, vorticity, and entropy waves in characteristic space to enhance accuracy and minimize numerical artifacts. In characteristic space, upwind schemes are used for acoustic waves, central schemes for vorticity and entropy waves, and Tangent of Hyper-bola for INterface Capturing (THINC) reconstruction for material interfaces and contact discontinuities (a subset of entropy waves). This approach prevents spurious vortices in periodic shear layers, accurately captures vortical structures in gas-gas and gas-liquid interactions, and improves the accuracy of shock-entropy wave interactions. In physical space, phasic densities are computed using THINC in regions of contact discontinuities and material interfaces, while tangential velocities are calculated with central schemes to improve vortical structures. An adaptive reconstruction technique is also introduced to mitigate oscillations near shocks, which arise from primitive variable reconstruction, by combining primitive and characteristic variable reconstructions with the liquid phase being identified using the stiffened gas parameter. The proposed multidimensional upwinding approach outperforms traditional schemes, demonstrating superior accuracy in capturing physical phenomena, reducing numerical artifacts, and better matching experimental results across complex test cases.
References in corpus (9)
- A positivity-preserving high-order weighted compact nonlinear scheme for compressible gas-liquid flows
- MFC: An open-source high-order multi-component, multi-phase, and multi-scale compressible flow solver
- High-Order Central-Upwind shock capturing scheme using a Boundary Variation Diminishing (BVD) Algorithm
- Gradient Based Reconstruction: Inviscid and viscous flux discretizations, shock capturing, and its application to single and multicomponent flows
- A low-dissipation shock-capturing framework with flexible nonlinear dissipation control
- First order hyperbolic approach for Anisotropic Diffusion equation
- On the role of spectral properties of viscous flux discretization for flow simulations on marginally resolved grids
- Implicit gradients based conservative numerical scheme for compressible flows
- Efficient high-order Gradient-based Reconstruction for compressible flows