CFD中文网

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

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

目前也不确定是openMPI的问题还是OpenFOAM太挑剔，所以各位使用集群计算的老师同学们，能不能看看你们集群的hostname是不是有下划线的，能让我定位一下问题的位置。

在学校的集群上装好OpenFOAM后发现只能单核计算，不能并行计算，报错如下：
While trying to create a regular expression of the node names
used in this application, the regex parser has detected the
presence of an illegal character in the following node name:

node: hpc_cn09

Node names must be composed of a combination of ascii letters,
digits, dots, and the hyphen ('-') character. See the following
for an explanation:

https://en.wikipedia.org/wiki/Hostname

Please correct the error and try again.
估计是计算节点名称“ hpc_cn09”不能有下划线吧？但是要改集群的名称太不现实了，请问有能不改名称的解决方法吗？

那可能是你的组分里根本就没有定义N2，你把inertSpecie N2换成inertSpecie H2O试试。

@CarmeloSun 增加一个温度标量方程，并与自带的能量方程通过一个能量源项来耦合

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

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

图1

图2

图3

@东岳 可能是有这个问题，但是我这几天把能用的格式的试了一遍，基本都是这个结果，我试着把动量源项加到压力方程里也是这个结果。似乎U方程算出来就是这个结果，真是奇怪

并且我的压力、密度在这些位置都是正常的，没有这样的突变，只有Ux是这种现象。

OpenFOAM版本：8
求解器：reactionFoam
我在reactingFoam的基础上做了一个两段多孔介质燃烧器求解器，需要在reactingFoam基础上增加一个温度标量方程，并与自带的能量方程通过一个能量源项来耦合，另外还需要一个fvOptions多孔介质阻力模型动量源项。所有的工作我都已经完成了，但是在测试的时候发现这个多孔介质动量源项对边界处，和两段多孔介质材料交界处的速度有不合理的影响。
如图：初始化均匀温度场300k，气体入口也是300k，速度0.7m/s，但是在入口处的第一个网格位置速度会低于0.7，同理，在中部两段材料交界的位置左右两个网格也会出现一边高一边低的现象。请问这个问题的原因是什么？有没有处理的方法？