@Samuel-Tu 雷诺数400 入口速度1m/s 粘度0.00125 直径0.5 那个算例和对比论文已经发在下面了

大神们好，我按照别人论文想模拟一个差不多贴底的圆柱绕流，用的是icofoam，只是把downwall变成粘性底层改成noslip条件，不知道为什么我把圆柱周围网格和底层网格越加密，升力系数会变得很大对不上，阻力系数出现很多凸起不平! 算例链接和对比文献如下：链接：https://pan.baidu.com/s/1KMeXCWMSQIES7VshrJYhzg 提取码：gr0b

@东岳 初步判定可能是of2.3.1版本对应的mpi环境配置问题，当输入mpirun命令后会提示：It looks like opal_init failed for some reason; your parallel process is likely to abort. There are many reasons that a parallel process can fail during opal_init; some of which are due to configuration or environment problems. This failure appears to be an internal failure; here's some additional information (which may only be relevant to an Open MPI developer):
然后具体解决办法在摸索中

各位Foamer大家好，最近我在CFD-onlie上下载了fabian_roesler等人贡献的基于焓-多孔介质方法模拟熔化过程的并行计算solver（https://www.cfd-online.com/Forums/openfoam-solving/93620-melting-problem-looking-appropriate-solvers-7.html ），根据solver的中的内容，我配置了对应的case。发现采用该solver单核串行执行case没有问题，但是采用并行执行case时，会提示：shift: can't shift that many。我搞了很久也没找到解决方法，恳请各位大佬帮我看下问题出现的原因，不胜感激！

对应的solver和case在附件中，of版本用的是2.3.1 parallel-convMeltFoam.zip case.zip

@BlookCFD 谢谢，这个问题已解决，这个主要涉及到interFoam系列算法子循环求解相连续性方程并修正质量流率，可以将surfaceScalarField rhoPhi在循环外声明创建，但是循环内部每次循环初始值要赋0.

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

同学，你好。很高兴发现你也从事interPhaseChangeFoam求解器相关研究工作。我想问一下，在interPhaseChangeFoam求解器主程序中，在pimple.loop()每次循环中中为什么要声明创建surfaceScalarField rhoPhi场，这样做的目的是什么？

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

@东岳

各位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!