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{matrix} (1) & 1 / ρ = C_{0} + C_{1} T + C_{2} T^{2} - C_{3} p - C_{4} p T \end{matrix}$

shrine

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

多项式状态方程

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

李东岳

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

unravelkaka

前辈你好，我最经也遇到这个问题，请问有时间交流一下吗？

CFD中文网

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