移植sprayFoam液相功能

[钟余]

@钟余 如果时间精度为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 影响没那么大。可能还是 @流体中二少年 提到的其他几个原因。