compressibleInterFoam求解器怎么加入空化模型

影川风

http://cfd-china.com/topic/95/compressibletwophasemixture库编写问题链接是一年前有前辈在论坛里问过的问题。我现在也面临着相似的问题，怎么编写compressibleTwoPhaseMixture库来调用空化模型？

影川风

我想在compressibleInterFoam求解器的基础上加入空化模型，控制方程的编写基本上没啥问题，不过参数调用这块我不怎么会，我请教过别人也看过曾经在这论坛里问过相似问题的前辈问的问题，不过还是不怎么理解。

李东岳

Okay。我看了一下，interPhaseChangeFoam只能用与不可压缩，所以你们打算添加可压缩VOF的方程；

一种方式是基于compressibleInterFoam附加空化模型;
一种方式是基于interPhaseChangeFoam更改速度压力方程适用于可压缩模型;

怎么编写compressibleTwoPhaseMixture库来调用空化模型？

这是一个非常大的问题，恐怕需要你自己完成。如果你选择第二种方法，你需要进行2步：

确认空化模型是否需要改动来可压缩适配，若不需要修改可跳过，若需要修改可修改phaseChangeTwoPhaseMixture.C公式内的函数;
incompressibleTwoPhaseMixture.C这个文件只是计算相界面的粘度、运动粘度等，个人认为不需要修改，可跳过; 不过最好确认下；
修改VOF模型中的速度和压力方程，可参考compressibleInterFoam中的处理方式，主要在压力方程中附加源项;

在做改动前，最后要非常熟悉interPhaseChangeFoam的流程，然后逐步的添加自己的模型。

影川风

@李东岳非常感谢东岳老师耐心的解答，我最近尝试着解决下这个问题

yuan_neu

第二种方案可行性更强
可以参考这个

Numerical Modelling of Evaporation and
Condensation Phenomena

论文里已经有了现成的compressible 的公式推导

小考拉

@影川风前辈您好，我目前也尝试在compressibleInterFoam里面添加空化模型，可以和前辈交流一下吗？

我爱梨涡

@小考拉你好，加个好友？1195143370，我也是interfoam加空化

上级

不知道这篇论文能否给点帮助（同学的论文）：
"叶秉晟. 针对非定常空化流动的多相可压缩求解器开发及应用[D]. 中国科学院大学，2019."

我爱梨涡

@上级老哥，俺下载不了啊如有时间，麻烦老哥发一份可否？1195143370@qq.com,拜谢！！！

shrine

@影川风我也在学习这个模型，请教：

看到tutorials里面只有层流的例子，这个不能计算紊流吗？
thermophysicalProperties三个文件中

phases (water air);

  pMin 10000;

sigma 0.07;

这两个变量不知道什么意思

还有水相的文件中

equationOfState
{
C (0.001278 -2.1055e-06 3.9689e-09 4.3772e-13 -2.0225e-16);
}

这个也不知道是啥。

李东岳

水相前面写的是：

equationOfState rPolynomial;

所以是多项式状态方程，在rPolynomial.H里面写的是：

Reciprocal polynomial equation of state for liquids and solids

    \f[
        1/\rho = C_0 + C_1 T + C_2 T^2 - C_3 p - C_4 p T
    \f]

    This polynomial for the reciprocal of the density provides a much better fit
    than the equivalent polynomial for the density and has the advantage that it
    support coefficient mixing to support liquid and solid mixtures in an
    efficient manner.

所以就是：
\begin{equation}
1/\rho = C_0 + C_1 T + C_2 T^2 - C_3 p - C_4 p T
\end{equation}

shrine

@李东岳在 compressibleInterFoam求解器怎么加入空化模型中说：

多项式状态方程

这是固液两项？如果算气液两相这么设？

李东岳

@shrine 什么都可以这么设定。主要是密度算的对就行。