administrators

看起来像是不可能的神奇错误，类似“煮大米饭一开锅竟然没有饭”。

原因未知。看起来像是不可能的神奇错误

我最近看了下，

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

但是上面这个方程不是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)
);

@hongjiewang

    IOList<scalar> utau
    (
        IOobject
        (
            "utau",
            runTime.constant(),
            "../postProcessing",
            mesh,
            IOobject::NO_READ,
            IOobject::AUTO_WRITE 
        ),
        1
    );
    
    utau[0] = b.weightedAverage(mesh().V()).value();

多谢！已更新

拉格朗日只能用paraFoam，paraFoam需要编译，所以感觉windows弄不了

这个就是发散了。具体原因需要更深层次的debug

@freeman 外挂求解器一般只能用于特定版本

这里面phi是严格守恒的，U不存在守恒的概念。所以你用phi = fvc::flux(U)，也可以，但是没有原来那样consistent

如果你的S在动量方程处理，你就不需要放在压力方程了

你这个最好用wedge网格来计算就不会出现这个问题了