用的求解器是意大利一哥们Davide Modesti写的是AUSM格式(上传文件),里面前台阶算例中用的无粘流,入口压力、速度、温度给的是1、1、3,热物性模型用的是const+hconst。相同热物性配置,我把他入口静温、静压和速度分别改成16095、130和572就会发散。经过debug,发现是AUSM格式里面c12的原因会导致数值爆炸。改用Sutherland+janaf模型后发现人工粘性打开就容易发生,温度出现负值。
有没有大佬能知道这是为什么呢,是AUSM格式本身鲁棒性不好么。
[rhoEnergyFoam-main.zip](Request Entity Too Large)