CFD中文网

@dzw05 你遇到的极值局部点距离给定压力条件的边界接近吗？我今天下午算的时候发现时间步过大时，距离fixedValue p_rgh边界最近的cell中心也会出现温度极值....调小了时间步后就没有了。我算的这段时间一直感觉温度和压强的联系太密切了...

@AppleKiller 不可压缩不可以

@李东岳 在 Simple算法中通量phi的计算 中说：

也就是说在边界处phi是不正确的。不知道是我没有找到phi边界修正，

边界处的phi不需要计算，是用户给的。加入速度1m/s，乘以网格边界面积，就是phi啊，并且不变。

如果是pressureInletVelocity 边界呢？ 代码里可以看出pressureInletVelocity 需要用边界上的phi来更新出新的边界速度。如果是这种边界边界上的phi又应该如何确认呢？

在初始化场的时候，比如这种

volScalarField test ( IOobject ( ... ) mesh //如果不写边界条件类型，除了特殊的边界条件，都为calculated，如果计算，就有值。一般情况下，边界条件很重要，要谨慎 scalar(0) );

Have no idea :chigua:

@王金成 我用的1812，可以编译SWAK，不过需要自己做些配置，不然会报错，过段时间我写个教程。

我最近看了下，

替代文字

上面这个方程应该是稳态算法的方程。里面MRF.DDt(U)用于处理$\Omega\times\bfU_I$。

96bfbd01-3d8c-469a-a303-a79d75b4d80c-image.png

但是上面这个方程不是OpenFOAM植入的方程。在瞬态情况下应该是下面这个：
$$
\frac{\p \bfU_I}{\p t}+\nabla\cdot(\bfU_I\bfU_R)+\Omega\times\bfU_I=...
$$
下面的瞬态代码分别对应：

tmp<fvVectorMatrix> tUEqn ( fvm::ddt(U) //第1项 + fvm::div(phi, U) //第2项 + MRF.DDt(U) //第3项 + turbulence->divDevReff(U) == fvOptions(U) );

:chitang:

@东岳 非常感谢！那我试着推推看这里:jingya:

@李东岳 额，从controlDict里的时间来看，多相计算的物理时间要远大于单相，但还是没能得到合理结果。我再算下去看看。