reactingFoam计算速度的问题

[Wust.liao]

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

图1

图2

图3

[Wust.liao]

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

[香柏树]

@wust-liao 你好，同学。我刚开始用reactingFoam计算甲烷预混燃烧，一直出现mpirun报错，请问你遇到过么，知道怎么解决么？万分感谢~

[Wust.liao]

@香柏树 你有试过自带的案例吗，先排除软件自身的问题

[香柏树]

@wust-liao 自带sandia的案例算了没问题，我是把sandia的rans改为les，换了一个复杂的燃烧室计算的。不知道哪块出问题了，导致只能断断续续地计算

[Wust.liao]

@香柏树看上去就是发散了，换个简单的燃烧室算算，模型一个一个往上加。不容易一步到位的

[yuky]

@香柏树 请问你用ReactingFoam ras改les成功了吗，想和您请教一下