A thermodynamically consistent and robust four-equation model for multi-phase multi-component compressible flows using ENO-type schemes including interface regularization

In this work, a concise and robust computational framework is proposed to simulate compressible multi-phase multi-component flows. To handle both shocks and material interfaces, a positivity-preserving ENO-type scheme is coupled with multi-phase interface regularization terms. The positivity-preserving limiter is conservative and is applied locally for minimal degradation of the baseline ENO-type scheme. The interface regularization terms are extended from the conservative diffuse interface (CDI) model to accommodate multi-phase, multi-component flows. The ENO-type scheme is designed to be consistent with the thermodynamic equilibrium assumptions of the four-equation multi-phase model, naturally enforcing the interface equilibrium condition - preventing oscillations in pressure, velocity, and temperature around isothermal material interfaces - without requiring additional equations for volume fraction or mixture equation of state parameters, as is commonly done for the five-equation model. Additionally, non-dilute species diffusion models are extended to the multi-phase, multi-component setting. We show that this consistent framework is equally applicable for regimes ranging from single-phase to multi-phase multi-component flows. The proposed models and numerical schemes are implemented in the highly parallel Hypersonic Task based Research (HTR) Solver, and high-resolution simulations are performed using both CPUs and GPUs.