CFD中文网

@尚善若水 好的，感谢大佬，他给出的算例是Li的21步H2，O2，N2的化学反应，我看下他的文献中有没有H2+O2+Ar的机理

@尚善若水 我现在去看一下，我尝试的两种化学反应机理，差距确实也挺大

@尚善若水 你好，这个我还没有，不过我用了2种基元反应模型，一种是17步反应的，一种是18步反应的，效果都不太行

@流体中二少年 大佬你好，我最近也尝试了Oran的6670Pa的算例，算出来却是这个效果，不知道为什么会这样，我把算例贴出来7Ar.zip

@李东岳 好的，感谢东岳老师

@尚善若水 感谢大佬，我看的是李健的气相爆轰数值模拟教程

在数值模拟中，记录爆轰波的胞格结构有两种方法，一种是记录每一个网格点上的最大历史压力，maxp，一种则是记录网格点上的最大涡量。但前者在规则反射的情况下，爆轰波的反射波强度较高，与三波点压力相近，会导致胞格结构模糊，因此，我想尝试第二种方法，即记录网格点上的最大涡量，不知如何实现。

@流体中二少年 好的我这就去看，感谢大佬

@流体中二少年 感谢大佬，我也询问过发布detonationFoam的大佬孙杰，他主要说过网格尺度不够对胞格大小有影响，以及实验胞格在20倍诱导区以上，模拟在10倍诱导区。我确实没有想过化学机理，rhoY及点火个数的影响。我研究的是马赫反射，胞格作为特征尺度很重要，我想询问一下，在Openfoam中有没有办法实现ZND爆轰波，或者2步反应模型，以自己调整胞格大小，类似这种。

@钟余 如果时间精度为1阶，则胞格网格更小

@流体中二少年 前辈你好，我现在主要是用开源的detonationFoam，但得到的结果感觉不太理想，例如胞格大小较实验小得多（采用基元反应模型，20Kpa，2H2+O2+2AR），并且胞格形状也不太稳定 ，我想咨询一下，这是因为什么原因造成的

@李东岳 好的，感谢东岳老师

@李东岳 李博，我看主程序，及rhoYEqn，不知道它是如何实现采用算子分裂分步法，将刚醒反应性的时间演化与对流项和扩散项的时间演化分开处理

Openfoam中如何实现算子分裂方法

我采用的是一个爆轰求解器，里面采用了自适应网格加密，如果想把计算结果作为另外一个计算模型的初始场，可不可以

由于detonationFoam是基于OpenFOAM开发的，因此直接使用了OpenFOAM提供的数值方案。对流通量采用带有HLLC-P近似黎曼求解器[36]的二阶MUSCL方案[46]进行计算，而扩散通量则采用二阶中心差分方案。采用算子分裂分步法[18]将刚性反应项的时间演化与对流项和扩散项的时间演化分开处理。在非反应流的时间推进中，采用了一阶隐式欧拉方案。我在文献中看见该求解器采用如上的数值格式设置，但不知道如何在Openfoam里进行设置，比如二阶MUSCL方案，是包含在HLLCP里面的吗

好的，感谢李博

@wangfei9088 好的，感谢大佬