reactingFoam计算速度的问题



  • OpenFOAM版本:8
    求解器:reactingFoam
    发现一个很头疼的问题,在使用reactingFoam计算甲烷-空气预混燃烧时,计算的速度似乎跟我设置的当量比有很大的关系,当我把当量比设为0.6时,计算会越算越慢,从最开始1s能算几个步长到后面10几秒才能算一个步长,当计算工况到10s以后基本就卡住不动了。但是当量比设为1的话计算就很正常,每个步长计算的时间基本都是一样。
    另外一个就是关于constant/chemistryProperties文件里的这个设置:equilibriumRateLimiter on/off,(图1)我没有找到关于这个开关的资料,但是代码说明里的建议是保持关闭。但是当我把这个开关打开时发现计算速度就会变得很快,不会随当量比发生明显的变化。只是计算的结果看上去很难受,火焰锋面处的体积放热率、组分浓度是锯齿状的(图2)。而正常情况因是图3的情况。

    d63462a7-247b-49f4-8e66-d2b8ce43e35c-图片.png
    图1

    eca3fdd1-6575-406b-aa01-96e8cfa5447b-图片.png
    图2

    3962fe71-0428-4710-b89f-c7dd67a00117-图片.png
    图3



  • 目前还没有找到关于计算速度变慢的原因,但是经过这段时间的调试发现了一些规律并尝试着把问题解决了。
    规律:当计算变慢时火焰峰值温度都在2500k左右,而实际上合理的温度应该在1800附近。只要当计算随着迭代到达了2500k速度马上就降下来了,最后基本就算不动了。
    解决办法:调小时间步长,比如用0.1的库朗数去计算,温度就不会到这么高,随着迭代温度回到1800k后可以适当增大库朗数。


Log in to reply