Your browser does not seem to support JavaScript. As a result, your viewing experience will be diminished, and you have been placed in read-only mode.
Please download a browser that supports JavaScript, or enable it if it's disabled (i.e. NoScript).
@李东岳 明白了,谢谢解答。李老师课程增加这个内容,感觉挺不错的。现在见到的确实基本都是rhoCentralFoam+reactingFoam(纯气)/+sprayFoam(气液)耦合。其实OpenFOAM本身的超声速求解能力感觉并不强大,好多功能。比如做成全隐式减少对小库朗数的依赖,增加高精度通量计算格式,高精度重构格式等都没有,但是其重心不在这边也可以理解。
胞格结构?李老师,我也不是太懂爆轰。我用这个求解器是想实现超声速流动下的反应横向射流。这两天忙着给老板写本子,有时间我看看这个算例
@尚善若水 是的,两个OpenFOAM公司重心都不在超音速那面了,所以也没开发也没包含。大佬玩哪个方向的,我看你好像做可压缩的比较多
胞格结构
对了,就是这个。我看他们爆轰都验证这个。我对爆轰理论也不太理解,不过对算法比较有兴趣。有机会再讨论 :-)
@李东岳 李老师,我主要做超声速燃烧
@xuluxi18 我现在不搞两相爆轰了,但是据我所知现在用of算爆轰燃烧已经非常成熟了,github上面有很多开源的求解器,例如Deepflame,blastFoam,detonationFoam,这些求解器有显式RK时间推进,耗散更小的黎曼求解器以及支持二维网格的AMR等等,比这个rhoReactingCentralFoam功能强大很多,你可以试试,祝好!
@李东岳 在 移植sprayFoam液相功能 中说:
rhoReactingCentralFoam
感谢李老师的回复,之前用过rhoReactingCentralFoam算爆轰胞格结构,效果还可以,但是这个求解器不支持二维网格的AMR并且KNP格式耗散大了点,后面我用一个HLLC-type的黎曼求解器和用上AMR,对惰性/反应性多组分激波管,旋转爆轰,斜爆轰进行验证,效果也非常好,链接文本
感谢分享,全面!
@流体中二少年 前辈你好,我现在主要是用开源的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 影响没那么大。可能还是 @流体中二少年 提到的其他几个原因。