A novel hybrid approach for accurate simulation of compressible multi-component flows across all-Mach number

TL;DR

This paper introduces a hybrid numerical method combining Godunov-type schemes and projection techniques to accurately simulate compressible multi-component flows across all Mach numbers, effectively handling shock waves and interfaces.

Contribution

The work presents a novel hybrid approach that integrates AUSM-based advection with a projection method, improving accuracy across all flow regimes including low Mach and high-speed flows.

Findings

Accurately simulates flows from incompressible to supersonic regimes.

Effectively captures shock waves and material interfaces without oscillations.

Demonstrates robustness in complex multiphase flow simulations.

Abstract

Numerical simulation of multi-component flow systems characterized by the simultaneous presence of pressure-velocity coupling and pressure-density coupling dominated regions remains a significant challenge in computational fluid dynamics. Thus, this work presents a novel approach that combines the Godunov-type scheme for high-speed flows with the projection solution procedure for incompressible flows to address this challenge. The proposed hybrid approach begins by splitting the inviscid flux into the advection part and the pressure part. The solution variables are first updated to their intermediate states by solving the advection part with the all-speed AUSM (Advection Upwind Splitting Method) Riemann solver. The advection flux in AUSM is modified to eliminate the pressure flux term that deteriorates the accuracy at the low Mach region. To prevent the advection flux from causing spurious velocities when surface tension is present, the pressure-velocity coupling term is modified to ensure it vanishes at material interfaces. Then, we derive the pressure Helmholtz equation to solve the final pressure and update the intermediate states to the solution variables at the next time step. The proposed hybrid approach retains the upwind property of the AUSM scheme for high Mach numbers while recovering central schemes and the standard projection solution for low Mach limits. To accurately resolve the complex flow structures including shock waves and material interfaces without numerical oscillations, a newly proposed homogenous ROUND (Reconstruction Operator on Unified Normalised-variable Diagram) reconstruction strategy is employed in this work. By simulating high-speed compressible multiphase flows and incompressible multiphase flows, this study demonstrates the ability of the proposed method to accurately handle flow regimes across all Mach numbers.