A Computational Modeling for Semi-Coupled Multiphase Flow in Atmospheric Icing Conditions

ABSTRACT: A two-dimensional second-order positivity-preserving finite volume upwind scheme is developed for a semi-coupled algorithm involving the air and droplet flow fields in the Eulerian frame, which shares the grid for each phase. Special emphasis is placed on the computational modeling, which is induced from a strongly-coupled algorithm that satisfies the strict hyperbolicity and its numerical scheme based on the Harten-Lax-van Leer-Contact solver preserving the positivity to handle multiphase flow in the Eulerian frame. The proposed modeling associated with the semi-coupled algorithm including the Navier-Stokes and droplet equations takes into account different boundary conditions on the solid surface for each phase. The verification and validation studies show that the new scheme can solve the air and droplet flow fields in fairly good agreement with the exact analytical solutions and experimental data. In particular, it accurately predicted the maximum value of the droplet impingement intensity near the stagnation region and the droplet impingement area.

摘要：针对欧拉框架下空气与液滴流场的半耦合算法（两相共享同一套计算网格），本文开发了二维二阶保正有限体积迎风格式。研究重点聚焦于计算建模部分，该建模源自满足严格双曲性的强耦合算法，且其数值格式基于保正的哈滕-拉克斯-范里尔-接触（Harten-Lax-van Leer-Contact, HLLC）求解器，用于处理欧拉框架下的多相流问题。本文所提出的建模方法结合了包含纳维-斯托克斯（Navier-Stokes）方程与液滴运动方程的半耦合算法，并考虑了各相在固体表面的不同边界条件。验证与确认研究表明，该新格式求解得到的空气与液滴流场结果与精确解析解及实验数据吻合度良好。尤为关键的是，该格式可精准预测驻点区域附近的液滴撞击强度最大值以及液滴撞击区域。