rhoCentralDyMFoam是否可以添加通量修正

[zxyf1]

各位CFDer你们好，

由于motion function library的限制，我在使用foam-extend-3.2中的动网格来进行建模，但是由于涉及到激波问题，需要使用rhoCentralFoam来进行求解，然而extend版本中并不包含动网格求解器rhoCentralDyMFoam, 因此需要我进行自定义。

目前我遇到的问题的是，我参照icoDyMFoam和openFoam-4.1中的rhoCentralDyMFoam来进行自定义，可是每次在动网格更新后温度都会升到一个极大值，之后求解器奔溃。现在我想参照pimpleDyMFoam来添加correctPhi.H来绕过这一极值问题，可是correctPhi似乎是作用于不可压流求解器的，对于这一部分我有所疑虑。在此向牛人请教，请问如此是否可行？

[CFD中文网]

correctPhi并不是动网格必须的条件，不是你必须的，我觉得最简单的是把你的运动库从extend移植到官方版本。

然后关于温度错误，你确定是correctPhi引起的？你完全可以关闭correctPhi求解。

[zxyf1]

感谢你的回复，把extend的运动库移植到官方版本，这个方法我不是没有想过，但是事实证明几乎是不可行的，我所使用的是layerAdditionRemoval，如果你能成功移植，务必请告知我，万分感谢。
关于correctPhi，目前我没有使用这段通量修正，我觉得温度错误极有可能是通量问题造成的，目前我想要确定的是，correctPhi能否用于高速可压缩流，特别是有激波间断的情况？

[CFD中文网]

当然可以移植:cheeky: ，你不需要把官方所有的热物理库移植过去，如果你明白求解的方程，只需要自己硬植入就可以了。
correctPhi用于在比如新开始一个时间步计算的情况下，计算守恒的phi场，动网格的算法主要在fvc::makeRelative()；correctPhi当然可以用于可压缩流，compressibleInterDyMFoam就植入了。
我查看了一下，rhoCentralDyMFoam默认并没有使用correctPhi，你试过网格不运动，结果正确否？

[zxyf1]

确实，默认的rhoCentralDyMFoam没有使用correctPhi, 而且在其他可压流求解器中有使用correctPhi，不过需要注意的是，由于我在这里涉及到的是激波模拟，因此有压力间断，但是correctPhi似乎是通过压力进行通量修正，我不敢确定这是否会对压力间断产生影响，因此目前暂时还没有使用。我尝试过网格不运动，即直接使用rhoCentralFoam情况似乎是正确的。