Development of accurate, robust liquid equations of state for multi-phase CFD simulations with a modified AUSM+-up scheme

• Published in: Computers & fluids 2013-04, Vol.77, p.166-180
• Main Authors: Robbins, D.J., Cant, R.S., Gladden, L.F.
• Format: Article
• Language: English
• Subjects: AUSM+ scheme; Compressible flow; Compressible liquid; Computational fluid dynamics; Computer simulation; Correlation; Equations of state; Fluid flow; Flux; Mathematical models; Multiphase flow; Solvers
• ISSN: 0045-7930; 1879-0747
• DOI: 10.1016/j.compfluid.2013.01.031

► Using EOS to reform the primitive variables for compressible multi-phase solver. ► Both phases assumed compressible. ► Eight equations of state compared for accuracy and robustness. ► Novel enthalpy coupling of cubic EOSs for robustness and speed. ► Modification of source term discretisation for smooth, robust solution profiles. Eight equations of state (EOS) have been evaluated for the simulation of compressible liquid water properties, based on empirical correlations, the principle of corresponding states and thermodynamic relations. The IAPWS-IF97 EOS for water was employed as the reference case. These EOSs were coupled to a modified AUSM+-up convective flux solver to determine flow profiles for three test cases of differing flow conditions. The impact of the non-viscous interaction term discretisation scheme, interfacial pressure method and selection of low-Mach number diffusion were also compared. It was shown that a consistent discretisation scheme using the AUSM+-up solver for both the convective flux and the non-viscous interfacial term demonstrated both robustness and accuracy whilst facilitating a computationally cheaper solution than discretisation of the interfacial term independently by a central scheme. The simple empirical correlations gave excellent results in comparison to the reference IAPWS-IF97 EOS and were recommended for developmental work involving water as a cheaper and more accurate EOS than the more commonly used stiffened-gas model. The correlations based on the principles of corresponding-states and the modified Peng–Robinson cubic EOS also demonstrated a high degree of accuracy, which is promising for future work with generic fluids. Further work will encompass extension of the solver to multiple dimensions and to account for other source terms such as surface tension, along with the incorporation of phase changes.