移植sprayFoam液相功能

[流体中二少年]

新手改求解器，由于老板要求引入液相功能，最近尝试把sprayFoam的液相功能移植到rhoReactingCentralFoam（rhoCentralFoam添加了reactingFoam的反应功能）上，编译求解器时，报的错误意思是方程数据类型不一致，我也大概知道了问题出在哪里，

solve
(
    fvm::ddt(rhoU)
  + fvc::div(phiUp)
  - rho()*g
  - parcels.SU(U)
);

这是动量方程的代码，我的理解是里面既有标量场，又有矢量场，所以数据类型不一致，希望老师们指点以下，给点思路，其实在能量方程那里遇到同样的问题

[Tens]

@流体中二少年 报错信息贴出来才好看

[流体中二少年]

@tens 不好意思，隔这么久回复你，求解器已经改好了，问题在于矩阵的形式不同，接下来准备验证求解器的可靠性

[xuluxi18]

大佬，最近我尝试在of7上使用rhoReactingCentralFoam，但编译完成后作者的example运行不了，报错如下，不知大佬遇到过这个问题没有

[李东岳]

上图就是发散了

[xuluxi18]

是的，应该是动网格有点问题，我把动网格关闭后能够成功运行。

[李东岳]

之前还没发现过这个rhoReactingCentralFoam，挺有意思。这个你们验证过爆轰么。你把自适应网格关了之后，那个高精度还能解析出来么，结果咋样。

[尚善若水]

@李东岳 李老师，基于rhoCentralFoam的求解器是不是从算法上就不太适合非规则网格和二阶格式。

[李东岳]

就不太适合非规则网格和二阶格式

什么意思，没懂。
非规则网格是什么，非结构网格么

[李东岳]

@xuluxi18 在 移植sprayFoam液相功能 中说：

是的，应该是动网格有点问题，我把动网格关闭后能够成功运行。

怎么样，后来的那种菱形的流场结构出来没（专业术语叫什么？）

[尚善若水]

@李东岳 是的呀李老师，截图是我看到的您写的rhoCentralFoam解析，俄罗斯hybridCentralSolvers的开发者也说，rhoCentralFoam不适合非结构网格，他的更好🤭

[李东岳]

@尚善若水 所有2阶格式在非结构网格拓展的时候，都会有一个潜在的问题，就是上游网格的继续上游的网格点的判断。对于超音速流动，格式模板的确定应该需要更加的精确，对于低音速流动因为方程并不是严格的双曲，因此限定性不会那么强。从理论上倒是可以解释。就像他们说的爆轰，如果不适用upwind类的格式，比如就是普通的reactingFoam，应该精度不够出不来这么细节的流场结果。所以他们就得用reactingFoam+upwind类格式才行。

不过我没实测过rhoCentralFoam在非结构网格上的计算效果。

但是帖子中讨论的rhoCentralFoam + 反应流算爆轰我觉得挺有意思，我在想要不要自己研究一下当做以后得课程资料。

[尚善若水]

@李东岳 明白了，谢谢解答。李老师课程增加这个内容，感觉挺不错的。现在见到的确实基本都是rhoCentralFoam+reactingFoam（纯气）/+sprayFoam（气液）耦合。其实OpenFOAM本身的超声速求解能力感觉并不强大，好多功能。比如做成全隐式减少对小库朗数的依赖，增加高精度通量计算格式，高精度重构格式等都没有，但是其重心不在这边也可以理解。

[xuluxi18]

胞格结构？李老师，我也不是太懂爆轰。我用这个求解器是想实现超声速流动下的反应横向射流。这两天忙着给老板写本子，有时间我看看这个算例

[李东岳]

@尚善若水 是的，两个OpenFOAM公司重心都不在超音速那面了，所以也没开发也没包含。大佬玩哪个方向的，我看你好像做可压缩的比较多

胞格结构

对了，就是这个。我看他们爆轰都验证这个。我对爆轰理论也不太理解，不过对算法比较有兴趣。有机会再讨论 :-)

[尚善若水]

@李东岳 李老师，我主要做超声速燃烧

[流体中二少年]

@xuluxi18 我现在不搞两相爆轰了，但是据我所知现在用of算爆轰燃烧已经非常成熟了，github上面有很多开源的求解器，例如Deepflame，blastFoam，detonationFoam，这些求解器有显式RK时间推进，耗散更小的黎曼求解器以及支持二维网格的AMR等等，比这个rhoReactingCentralFoam功能强大很多，你可以试试，祝好！

[流体中二少年]

@李东岳 在 移植sprayFoam液相功能 中说：

rhoReactingCentralFoam

感谢李老师的回复，之前用过rhoReactingCentralFoam算爆轰胞格结构，效果还可以，但是这个求解器不支持二维网格的AMR并且KNP格式耗散大了点，后面我用一个HLLC-type的黎曼求解器和用上AMR，对惰性/反应性多组分激波管，旋转爆轰，斜爆轰进行验证，效果也非常好，链接文本

[李东岳]

感谢分享，全面！