A diffuse-interface method for reducing spurious pressure oscillations in multicomponent transcritical flow simulations

Abstract Spurious pressure oscillations are a numerical instability that is commonly observed in multiphase and multicomponent flow problems near the critical point. A diffuse-interface model has been created to simulate transcritical mixing in multispecies and multiphase systems where spurious pressure oscillations pose a serious challenge in achieving convergence. To reduce the spurious pressure oscillations, three methods have been proposed and implemented in the present paper: (1) artificially thickening the interface between different species, (2) reconstruction of the primitive variables in the characteristics space, and (3) developing a hybrid method (HY) that switches between quasi-conservative (QC) double-flux (DF) and the classical fully-conservative numerical procedures based on the changes in the effective specific heat ratio and the effective reference internal energy. The cases considered include (1) a transcritical shock tube problem between nitrogen and n-dodecane, (2) a near-critical shock-droplet interaction, and (3) the transcritical Spray A benchmark case of Engine Combustion Network (ECN) that represents the injection of n-dodecane into a high-pressure and high-temperature nitrogen environment. Characteristic-wise reconstruction of the primitive variables in our proposed fully-conservative method tends to reduce the spurious pressure oscillations compared to primitive-wise reconstruction including the reconstruction of the speed of sound. The HY model is found to effectively reduce the magnitude of the spurious pressure oscillations, the loss of energy conservation, and provides more accurate results for all tested cases. The artificial thickening of the interface between different species is found to reduce spurious pressure oscillations and reduce energy conservation loss when used in conjunction with a QC method (HY or DF).