CFD中文网

帖子

• ---

  OpenFOAM膨胀喷管边界条件设置
  李 李东岳

  @xuluxi18 才看见，rhoCentralFoam能跑通，rhoPimpleFoam跑不通么

• ---

  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小时，下面是我处理到最好的解了。

  算例在这里（写好求解器后一键出图），感兴趣可以研究研究

  sinq.tar.xz

  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的方程可以用俩种方法计算库朗数：

  1. 特征速度
  2. 传输速度

  所以，rhoCentralFoam应该就是用的特征速度计算的库朗数，而不是传输速度。如果想看传输度的库朗数，需要用compressibleCourant.H来计算。

• ---

  buoyantfoam计算发散求助
  李 李东岳

  Quality是0.354

  网格质量比较差。一般是看综合的Quality