rhoCentralDyMFoam是否可以添加通量修正



  • 各位CFDer你们好,

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

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



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

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



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



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



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


Log in to reply