@xuluxi18 才看见,rhoCentralFoam能跑通,rhoPimpleFoam跑不通么
教授
帖子
-
OpenFOAM膨胀喷管边界条件设置 -
buoyantfoam计算发散求助是的,越复杂的求解器越容易发散。
-
如何根据stl模型定义cellZone?你需要看cellZone,不是看stl。paraFoam可以打开cellZone
-
buoyantBoussinesqSimpleFoam的p_rgh收敛问题压力是没下去。原因不知道了,需要详细debug
-
buoyantBoussinesqSimpleFoam的p_rgh收敛问题压力残差确实没下去。你的残差有几条线甚至在0,看起来没有求解。
-
buoyantfoam计算发散求助0.6可以了。还发散需要debug算例了。
-
OpenFOAM求解sinq传输方程搞出来了,一个小bug卡2小时
-
OpenFOAM求解sinq传输方程$$
\frac{\p q}{\p t}+\frac{\p \mathrm{sin}(q)}{\p x}=0
$$$$
-1.5<x<1.5
$$$$
if x<0, q = 0.78539
$$$$
if x>0, q = 11.7809
$$方程看起来简单,没想到还搞了2小时,下面是我处理到最好的解了。
算例在这里(写好求解器后一键出图),感兴趣可以研究研究
Leveque那面也有更多的信息:https://www.clawpack.org/riemann_book/html/Nonconvex_scalar.html#Sinusoidal-flux
-
OpenFOAM计算斜激波,波后压力振荡一说到这个库朗数,我看了一下,这面是用最大波速来计算的CFL。我用compressibleCourant对比了下。二者一个0.17,一个0.14。还有点差距。波速是求解之前预估的,phi是求解之后。我估计是差距的来源。不过没有细看代码,自己想的。
我今天在准备这个算例的时候发现,http://dyfluid.com/class.html#hyperbolic 特征速度和传输速度有2倍的差异。针对这种偏hyperbolic的方程可以用俩种方法计算库朗数:
- 特征速度
- 传输速度
所以,rhoCentralFoam应该就是用的特征速度计算的库朗数,而不是传输速度。如果想看传输度的库朗数,需要用
compressibleCourant.H
来计算。 -
buoyantfoam计算发散求助Quality是0.354
网格质量比较差。一般是看综合的Quality