interPhaseChangeFoam空化求解,负绝对压力问题
-
最近在做一个2D翼形的模拟,求解器用的是interPhaseChangeFoam
边界条件:
入口:总压,用的是prghTotalPressure;速度pressureInletVelocity 根据phi计算流速为5.33m/s;
出口:zeroGradient velocity & p_rgh湍流模型使用的是schnerrsauer,选择合适的n值进行计算
turbulence model用的是DDES模型
为了简化问题,暂时关闭了重力项和表面张力
在求解压力泊松方程时候,出现了压力负值的情况(我用的是绝对压力),经观察,出现负值的情况为alpha出现尖点的位置(for instance cavitation collapse)
传质源项的隐式处理,增强了压力泊松方程的对角占优特性,这是不可压流体求解器,理论上说,压力泊松是椭圆形方程,如果场内出现强扰动会瞬时传遍全场(空泡破碎产生高压),但是相方程是如何对泊松方程影响的我还不清楚,如何解决这个问题呢?
-
相方程是如何对泊松方程影响的我还不清楚
这个问题之前的不清楚原因。
多相流的压力方程来源于两相的不可压缩连续性方程如:
\begin{equation}
\frac{\partial \alpha}{\partial t}+\nabla\cdot(\alpha\mathbf{U}_ \alpha)=0
\end{equation}
\begin{equation}
\frac{\partial \beta}{\partial t}+\nabla\cdot(\beta\mathbf{U}_ \beta)=0
\end{equation}
两个加起来就是divU=0,要依据这个组建压力方程。
对于你的有源项的相连续性方程,两个加起来之后就不再是divU=0,而是divU=S,基于这个方程出来的就是你的压力方程。所以说:相方程左边组建速度散度,相方程右边组建压力方程源项。大体是这个关系。