移植sprayFoam液相功能

[钟余]

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

[尚善若水]

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

[钟余]

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

[尚善若水]

@钟余 Cheng也用的Oran机理

[钟余]

@尚善若水 是的，我也看到了，但目前也还没找到Oran的机理

[流体中二少年]

@钟余 不好意思，最近比较忙，没想到你已经做了这么多的工作。是这样的，计算域长点，到后面会稳定两个胞格，但是我用了取巧的方法，你稍微给大点宽度，现在是60mm，你给个61.8mm试试，然后给4个点火源。这是我最近用5阶weno算的结果： ，和openfoam算的差不多。

[尚善若水]

@流体中二少年 这是用PeleC算的吗？我觉得机理可能影响也很大，请问能分享下Oran机理文件吗？cfd@ik.me 谢谢指教。

[流体中二少年]

@尚善若水 对，用的PeleC。很抱歉，这个机理也是从别人那里来的，我没办法直接分享给你，不过我的经验是无论Li的机理还是Burke的机理，都能算出这个效果。

[尚善若水]

@流体中二少年 好的，多谢啦

[钟余]

@流体中二少年 好的，我现在去尝试一下。有好的结果我再分享出来。感谢大佬能为我们解惑