关于of-V中interCondensatingEvaporatingFoam处理密度疑问

[沐风]

各位Foamer大家好，我目前在尝试开发固相和液相密度不同时熔化/凝固问题的模拟（melting or solidification with density change）。基于这一问题，openfoam中可以参考的求解器有可压缩+传热的buoyantPimpleFoam求解器和openfoam-v版本中interCondensatingEvaporatingFoam求解器。buoyantPimpleFoam求解器中涉及到关于密度和压力的状态方程，关于该求解器中的thermo.rho我目前还不能理解。因此目前我主要关注同样和我问题相关的interCondensatingEvaporatingFoam求解器，关于该求解器，我有如下三点疑问：

1. 为什么在求解器主程序中pimple.loop中每次循环中要新声明一个rhophi面标量场，我曾将该声明放到creatfields文件中，但发现程序运行结果存在不同；
2. 为什么在速度方程中添加了fvm::Sp(fvc::ddt(rho) + fvc::div(rhoPhi), U)这一项，我看CFD-online中有人给出的解释为：这项是一stabilisation项，如果连续性方程满足该项为0，说实话对此我不是特别理解
3. 为什么在构建interCondensatingEvaporatingFoam求解器在构建压力方程中，没有密度对时间偏导这一项，是基于以上问题1和问题2的处理，将密度对时间偏导这一项可以省略吗？而我看buoyantPimpleFoam求解器中有密度对时间偏导这一项。

最后，希望中文网的各位大佬指点迷津，All suggeations are appreciated!

[沐风]

@东岳

[沐风]

经过我的研究，第二个问题涉及到守恒和非守恒形式的动量方程转变，第三个问题中密度对时间的偏微分可以通过连续性方程转换到全微分中，以上两个问题均可以从东岳老师写的‘无痛苦N-S方程’中找到答案，现在就差第一个问题了

[李东岳]

@沐风 在 关于of-V中interCondensatingEvaporatingFoam处理密度疑问 中说：

为什么在求解器主程序中pimple.loop中每次循环中要新声明一个rhophi面标量场，我曾将该声明放到creatfields文件中，但发现程序运行结果存在不同；

你可以放在createFields中，但是要在循环中置0

[沐风]

@东岳 经过昨天阅读输出主程序rhophi后的代码，发现这个rhophi声明并赋值为0的目的是为了在子循环求解相连续性方程并修正质量流率。东岳老师之前写的‘VOF 模型，interFoam 计算流程’提供了很好的参考。PS:发现提的三个问题好弱，吾辈还需多读书！

[李东岳]

多谢老铁对我资料的认可

[fireztw]

@沐风 您好，我最近也在研究可以计算凝固问题的模型，初步的模拟计划就是水滴在低温平板上逐渐凝固成冰，最好考虑水和冰的密度差异，请问楼主现在可以计算了吗？有点问题想请教一下，非常感谢！