@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插值,处理压缩性较强的问题不是很给力(会振荡),所以想找个相对简单点的算例测试一下,暂时没找到
-
如题,最近写了个新的求解器,想验证一下算法,有哪位大佬指教一下
-
@李东岳 哈哈哈,物理隔绝了
-
@yfleung jhd上的代码有bug,且效率低,我自己重写可用。
https://cloud.tsinghua.edu.cn/d/05b0149e37904e7ca61c/ -
@izumi 我理解是这样:SRF给边界是相对速度,所以跟着坐标系一起转。MRF是绝对速度,所以固定不动。
-
我看到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。请问对比时域的绝对值是不是没意义呢?
-
水听器是没于水中测声压的探头。我有一组实验数据,我想把模拟的检测结果和实验对比。我不知道后处理的时候应该直接用监测点的瞬时值呢?还是应该用压力脉动值?请大家给些建议
-
如题,比如相体积分数,有人尝试过么?
-
如题,目前的方法是blockMesh先生成方形单层mesh,再利用extrudeMesh生成360的3D网格。
但是出现个问题,里面的面不会自动的成为interior,不太清楚怎么处理。以下是extrudeMeshDict
constructFrom patch; sourceCase "$FOAM_CASE"; sourcePatches (side1); exposedPatchName side2; extrudeModel sector; nLayers 10;//90; sectorCoeffs { point (0 0 0); axis (0 1 0); angle 20;//360; } flipNormals false; mergeFaces false;
-
@南风 他这个意思是群泡传质速率,不是单泡的。另外这个值与你其他参数取值也有关系
-
通过求解ode方程,得到了一系列采样点,想对其进行fft分析,计算出主频对应的幅值,不知道怎么用of里的fft
-
看到个帖子不错,记录一下。关键词:householder transform
-
有人写过此边界么?
-
前辈好,这个边界目的是为了强制使声偶极子在此处压力为0。我不确定加了会不会导致计算发散。不知有什么办法定义一个域内面积为0的面,然后加边界不?
-
利用baffleDict生成了一对baffle面,上下两个面的边界可以给,但是侧面也需要边界条件。由于是无限薄,侧面其实就是一个edge。我查了貌似没有办法加edge的边界。如果以加入源项的方式是否能等效实现边界条件的加入呢?
-
@东岳 是的,这个operator[]不会每个时间步都去读,感觉这样改起来是个无底洞呢。不知是否能设置个贴体坐标,所有的点都相对其运动起来
-
@wwzhao 不知是否能设置贴体坐标,所有的probes都相对其移动呢?
-
@东岳 请问前辈,我如何能把vectorField转化成Ostream呢?
-
多谢前辈回复,我是想这么弄,直接在controlDict里改probeLocations可否?这些坐标点应该是每一个时间步结束后都会读一遍吧?
-
动网格问题,我想设置一系列点,坐标在相对移动的面一定距离,也就是说监测点probe是移动的,有什么好的方法能实现不?
-
ESI版本1906,有个bug,在compressibleInterFoam下,将turbulence models换成twoPhaseTransport时,除了laminar外都不能算,直接报错。bug提交在8个月之前,但是并没有人回应,不知其他人有没有研究。
-
@东岳 前辈好,最近在算les空化模拟,也卡在计算效率太低上了TAT
-
@random_ran 请问找到进一步提高效率的方法了么?
-
用interPhaseChangeFoam求解器是不能算出所谓“压力波”,因为其求解的压力方程是椭圆形pde,你看到的pressure wiggles是由于算法导致的,是非物理的压力震荡。
若想研究空泡中的shock wave,你得换个求解器cavitatingFoam。空泡群是用的EoS模型来描述的 -
U是体积平均速度
U = U_l\alpha_l + U_v\alpha_v
其中U_l和U_v都是filtered单相流体的速度
-
@random_ran “数学公式”选项,码命令即可
-
好!
但为啥不直接用pyFoam?😂 -
@wwzhao 老师你好~这算是我自己给自己结个题😂
看起来,这文章讲的还是松弛因子的无关性修正,这问题确实有点意思,我这个结果是按照那个提出比较早的方案做的,效果还可以。 -
看样子是压缩过度了,如果界面太sharp,CSF模型计算曲率会出问题,代数重构方法的通病
-
植入后的cavity测试,与时间步无关
-
如题,MMIM即modified momentum interpolation method
原始文献:DISCUSSION ON MOMENTUM INTERPOLATION METHOD FOR COLLOCATED GRIDS OF INCOMPRESSIBLE FLOW
OF中的动量插值方法为OMIM,其中还加入了一项ddtCorr,该项可以消除时间步的影响。
但是在MMIM文献的分析中,表示ddtCorr并不能消除时间步长的影响。而且这个ddtCorr还会带来额外的耗散,如果做LES模拟会降低精度。
最近我尝试实现MMIM,按着文章中所说,发现实现不了。问题主要是更新phi和U的时候二者相差甚远,会发散
有没有有兴趣的同学一起来搞一下?
-
@东岳 如此说来,除了加密网格,没有好的解决方案啦
多相流超声速计算出现负温度
用较小的时间步长,结果反而出问题了
多相流超声速计算出现负温度
求问,有没有两相可压缩等熵流的验证算例?
求问,有没有两相可压缩等熵流的验证算例?
求问,有没有两相可压缩等熵流的验证算例?
关于MRF的疑问
请教:第三代涡识别方法Liutex的代码在linux中运行的问题
关于MRF的疑问
复制之前算例的场信息到新算例,如何改变边界条件??
openfoam添加湍流入口方法?
网格生成与合成的问题
网格生成与合成的问题
并行效率疑问
将多相流模型完全双曲化都有什么好处呢?
压力后处理问题
压力后处理问题
压力后处理问题
压力后处理问题
压力后处理问题
如何利用setField设置非均匀分布的初场?
如何用extrudeMesh生成360度的axisymmetric 3D网格
关于interPhaseChangeFoam里面的一些疑问
有人用过of里的fft吗
anti-symmetry边界?
anti-symmetry边界?
如何在domain内加边界条件?
如何在domain内加边界条件?
如何将计算区域内部的某个面指定为边界,然后加边界条件
如何设置移动的probe点?
如何设置移动的probe点?
如何设置移动的probe点?
如何设置移动的probe点?
如何设置移动的probe点?
compressibleInterFoam + twoPhaseTransport
并行效率疑问
并行效率疑问
并行效率疑问
interPhaseChangeFoam能否求解出空泡溃灭时的压力波
VOF动量方程中的速度指的是哪一相?
求OpenFOAM中MULES算法的解析
有人尝试过在OF中实现MMIM方法么?
【分享】 Residual 动态监测
有人尝试过在OF中实现MMIM方法么?
有人尝试过在OF中实现MMIM方法么?
两相流VOF自由面处总会出现锯齿状
关于MULES的疑问
有人尝试过在OF中实现MMIM方法么?
有人尝试过在OF中实现MMIM方法么?
关于bubble flow中的通量计算问题