This work focuses on the formulation of a four-equation model for simulating unsteady two-phase mixtures with phase transition and strong discontinuities. The main assumption consists in a homogeneous temperature, pressure and velocity fields between the two phases. Specifically, we present the extension of a residual distribution scheme to solve a four-equation two-phase system with phase transition written in a non-conservative form, i.e. in terms of internal energy instead of the classical total energy approach. This non-conservative formulation allows avoiding the classical oscillations obtained by many approaches, that might appear for the pressure profile across contact discontinuities. The proposed method relies on a Finite Element based Residual Distribution scheme which is designed for an explicit second-order time stepping. We test the non-conservative Residual Distribution scheme on several benchmark problems and assess the results via a cross-validation with the approximated solution obtained via a conservative approach, based on a HLLC scheme. Furthermore, we check both methods for mesh convergence and show the effective robustness on very severe test cases, that involve both problems with and without phase transition.
翻译:这项工作的重点是为模拟具有阶段过渡和强烈不连续状态的不稳定的两阶段混合物制定一个四等分模型,主要假设是两个阶段之间的同质温度、压力和速度字段。具体地说,我们提出一个剩余分配办法的延伸,以解决四等分两阶段的系统,以非保守的形式,即以内部能源而不是传统总能源法的形式,进行分阶段过渡,即以内部能源而不是传统总能源法的形式,对四等分制进行计算。这种非保守的提法可以避免许多方法获得的典型振动,这些方法可能出现在接触不连续状态的压力剖面上。提议的方法依赖于一个基于精密元素的残余分配办法,该办法的设计是明确的第二顺序时间步。我们根据几个基准问题测试非保守剩余分配办法,并通过一种保守办法,以HLLC办法为基础,对结果进行交叉校准。此外,我们检查两种组合方法,并显示在非常严重的试验案例中的有效稳健性,这既涉及有阶段过渡的问题,也不涉及阶段过渡的问题。