关于rhoCentralFoam计算高超声速问题时温度出现负值的原因



  • 问题描述:
    以openfoam自带算例斜劈(或其他算例均可)为例,修改如下:
    M=6或10等,热物理参数使用pefectgas, 入口温度为60K,静压1270pa,模拟试验状态,因为总温有限,用等熵公式算出静温一般较低;边界条件不变
    问题是,不管Dt和CFL数怎么调,很容易在计算若干步出现初始T0为负值的情况;
    个人理解:M较大时,流场参数的梯度较大,过低的初始温度和静压加剧了计算不稳定,若改成自然来流,则能正常计算,通过fvOptions对温度进行限制根本不起作用,请大家指导一二



  • @longming 你好,我在使用rhoPimpleFoam时遇到了同样的问题,请问你的问题现在解决了吗?



  • 你可以试试修改flux limiter。