@pengdi 多相密度基方法可参考Saurel的文章:doi:10.1017/S0022112008002061,基于迎风格式应对多可压缩流动问题处理效果较好,求解器是自己开发的。
yhdthu
帖子
-
多相流超声速计算出现负温度 -
用较小的时间步长,结果反而出问题了动量插值里引入的压力梯度三阶导项前面的系数乘了deltaT,如果deltaT太小压力光顺效果会变差。
-
多相流超声速计算出现负温度用CompressibleIsoInterFoam算超声速射流有两个本质上无法解决的问题:
- 压力基Rhie-Chow插值属于中心插值,可压缩出现shock时会振荡,如果开启了choi修正加入了耗散会把激波抹平(当然看你关注什么,如果只是关注推力之类的统计量可能没啥影响,但是研究流场结构就够呛了)
- 如果你认真推一遍公式,会发现它的能量(温度)方程压根就是错的,算出负温度不足为奇。另外,我记得thermo里的mixture函数在每个时间步计算时都会initial一下,温度会重新算,可以按“Negative initial”关键字搜索定位一下代码。
总之吧,这个求解器勉强能用,如果你换套网格、换个格式(比如snGrad用limited)或许到这个时刻就能算通了,然而隐患一直都在,不知啥时候又会出现。
如果真想把浸没射流算好,建议用多相的density-based solver,我们算的效果长这样:
-
求问,有没有两相可压缩等熵流的验证算例?@yhdthu 大佬别闹,求赐教
-
求问,有没有两相可压缩等熵流的验证算例?@heike256 东岳哥好,我是把两相密度基算法(p-relax, 5方程)用压力基的求解套路实现一下,受制于Rhie-Chow插值,处理压缩性较强的问题不是很给力(会振荡),所以想找个相对简单点的算例测试一下,暂时没找到
-
求问,有没有两相可压缩等熵流的验证算例?如题,最近写了个新的求解器,想验证一下算法,有哪位大佬指教一下
-
关于MRF的疑问@李东岳 哈哈哈,物理隔绝了
-
请教:第三代涡识别方法Liutex的代码在linux中运行的问题@yfleung jhd上的代码有bug,且效率低,我自己重写可用。
https://cloud.tsinghua.edu.cn/d/05b0149e37904e7ca61c/ -
关于MRF的疑问@izumi 我理解是这样:SRF给边界是相对速度,所以跟着坐标系一起转。MRF是绝对速度,所以固定不动。
-
复制之前算例的场信息到新算例,如何改变边界条件?? -
openfoam添加湍流入口方法?我看到DFSEM入口有bug?
讨论地址:https://www.cfd-online.com/Forums/openfoam-solving/177711-turbulentdfseminlet.html
有问题的是11式,本来是个无量纲参数,根号出来其实有量纲,不知有人测试过没?
-
网格生成与合成的问题模型还是很简单的,就是一个柱体+方形,但是我用icem画必须要用两个O grid,网格过渡很不好(至少我看起来很难受),不知有什么其他思路?
-
网格生成与合成的问题由于对网格质量和控制的要求,我把流体域分成了两部分进行网格划分,显然二者有个重合的界面。
我先用blockMesh生成了核心区的网格。之后我想用snappyhexmesh生成非核心区网格。
我想请问,是否可以将已经有了的核心区的界面导入,相当于固定了面网格,去生成非核心区网格呢?我目前没查到相关方法,我试了mergeMesh,但效果很不好,不知各位有何高见?
-
并行效率疑问所以最终结论是PCG+DIC效率最高是么?
PCG+GAMG是不是稳定性好一些? -
将多相流模型完全双曲化都有什么好处呢?如题。
-
压力后处理问题@xpqiu 好的,谢谢邱老师,我再试试😄
-
压力后处理问题@xpqiu 好的,我的网格在观测点处太粗了,时间步长选择的和探针的采样频率一致,是不是需要比这个值还小呢?
-
压力后处理问题@xpqiu 是的,定性的可以对上。计算中波动被磨平了很多,您认为是怎么回事呢?
-
压力后处理问题邱老师好,我做了监测点的PSD。请问对比时域的绝对值是不是没意义呢?
-
压力后处理问题水听器是没于水中测声压的探头。我有一组实验数据,我想把模拟的检测结果和实验对比。我不知道后处理的时候应该直接用监测点的瞬时值呢?还是应该用压力脉动值?请大家给些建议