移植sprayFoam液相功能

[钟余]

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

[钟余]

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

[流体中二少年]

@钟余 原因比较复杂，大概有：网格尺寸(一般诱导区长度或者半反应区长度需要有10个以上的网格来解析)，化学机理(适合高温高压)，还有就是detonationFoam的rhoY并不没有进行通量计算，当然我不知这个影响大不大(这个你可以参考deepflame的通量格式)，还有就是和点火个数有关(这个你得自己试试换下点火个数)。我没跑过20kpa的算例，你可以试试经典Oran得6670pa算例。

[钟余]

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

[流体中二少年]

@钟余 你是想映射一段ZND爆轰作为初始流场？这个可以用sdtoolbox算一段一维爆轰数据，然后set到流场里面，这个不难。其次，这个两步反应模型，我没具体了解过，不过你可以参考下这个https://github.com/risinyoung/rhoCentral1StepFoam。

[钟余]

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

[李东岳]

优质讨论！！

[尚善若水]

我的计算经验，rhoYi 不进行 flux splitting 影响没那么大。可能还是 @流体中二少年 提到的其他几个原因。

[钟余]

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

[尚善若水]

@钟余 你有用oran机理测试吗？你看这篇文章，机理影响还挺大的

[钟余]

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

[钟余]

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

[尚善若水]

@钟余 oran的机理我也没找到，你可以试下detonationFoam用的Li的机理，detonationFoam算的挺好的

[钟余]

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

[尚善若水]

@钟余 如果有结果可以分享下，我也很好奇😂

[钟余]

@尚善若水 我又尝试了Marinov的二十五步反应机理，想复刻Oran的6670Pa的算例，结果还是不行

[尚善若水]

@钟余 好像和你前面发的差别也挺大？

[钟余]

@尚善若水 我又询问了detonationFoam的发布者孙杰，他也给了我解释：这个是当时的胞格文件，这个算例非常难算，因为他的诱导区长度很大，反应物的活性很差，很难起爆，起爆能量高了又会受初始起爆区的影响很久胞格数目和初始的高压区数目一样，没法过渡到自身特征的胞格尺寸，当时我也卡了很久，最后的经验是计算域首先一定要长，另外就是初始受起爆区影响的时候（看胞格数目是不是和高压区的数目一致），把自适应网格层数减小利于拜托初始设置的影响，等摆脱这个阶段后再改回去保证网格收敛。
但我还是不知道有些文献是如何设置的，在非常短的时间内就可以形成稳定胞格结构。

[尚善若水]

@钟余 可能oran的机理更好起爆？或者插值一维计算结果，你看这个BYCFoam: An Improved Solver for Rotating Detonation Engines Based on OpenFOAM。

[钟余]

@尚善若水 我看了下，他引入一维的结果后，好像确实很快就形成了54*30mm的稳定胞格，这种应该可以去除初始高压区的影响或许我也可以去尝试一下。他是先运用一维网格计算得到的一维CJ爆轰波，再利用mapfields映射的吗，我对这方面的操作不是很了解。