XiFoam计算气体爆炸增加障碍物后发散

[xuqiming]

我前期基于OF7 XiFoam进行了氢气爆炸模拟，与实验结果一致性较好（case1所示），收敛性很好。现在加入障碍物后（case2），计算在进行到5ms左右发散。发散时能看见障碍附近出现很大的湍动能。我感觉是k-epsilon不适合这种工况，修改为SST仍然发散。请问各位大佬有什么好的建议吗。
使用的基本设置如下

• k-epsilon
• 0.4 mm网格，（测试过无关性）
• k-1e-3，zerogradient
• epsilon-0.5，epsilonWallFunction
• combustion-Gulder,Xi-transport
  

[李东岳]

我不太熟悉这个算法，所以不能给出针对性的意见。

不过这个求解器用的人非常少，这么多年的演变也很少，以后可能会删除了。不知道reactingFoam能不能实现你要的功能。

[xuqiming]

@李东岳 感谢李老师的回复。reactingFoam求解器之前实现了无障碍下的计算和计算精度。但是由于反应机制的刚性和反应步数问题，导致计算量过大，同时reactingFoam在进行计算时压力震荡剧烈，一直不清楚是离散格式设置的问题还是求解器本身的问题，所以最近在尝试XiFoam这个似乎更简单一些的求解器。

[xuqiming]

经过查看该求解器的原始文献发现，当Xi模型选择为transport模型时，使用如下的方程进行计算，经过尝试发现的确是该方程影响了计算的收敛性。个人认为可能是由于该方程中没有扩散项导致的不稳定性。将湍流方程修改为RNG KEPSILON和relizable KEPSILON时能够进行计算，但是计算结果欠佳。希望有大佬能够予以指点。

[李东岳]

这个应该是模型里面的变量b吧？这个模型确实没有扩散项。就跟相方程一样。

[好名字]

@xuqiming 您好打扰您了，我在用XiFoam做氢-空气预混燃烧的时候发现点火强度strength对燃烧过程影响很大，但是我没有找到strength的详细定义，希望能得到您的指导
ignitionSites
(
{
location (0.055 0.041 0.0005);
diameter 0.003;
start 0;
duration 0.001;
strength 2;
}
);

[xuqiming]

@好名字 根据我查阅的资料，strength 是点火能量的表征，影响的是初始时刻的点火效果，也就是b值和燃烧产物的质量分数之类的物理量。

[好名字]

@xuqiming 谢谢您的解答