CFD中文网

李

一个2D的算例更好。可以调节进口速度，来测试是否是正常的。
1D算例由于通量守恒，不好测试

李

估计你下载的不是v11的版本。你要卸载互相适配的版本。

李

感谢分享！！

李

这个求解器不适用于openfoam11

李

屌爆了老铁，聊天记录已收藏，明天就发JFM

李

测试这个边界条件，最好使用一个基础的求解器比如pisoFoam，使用自带的2D pitzdaily算例来单核运行测试，测试通过后将其衍生到你的算例。我在近期可以顺着我的路子弄一个。不过如果你可以弄一个上传的话，可以加快这个进程。

李

你有测试pisoFoam这种单相流么

李

https://cfd-china.com/topic/8018

算例发上来我看看

李

constrainPressure主要是更新压力第二类边界条件，公式如下：

\begin{equation}
\left( \nabla p_{rgh} \right)_f \cdot\bfn_f=
\frac{\left(\mathbf{HbyA}_f^{*}- \frac{1}{{{A^n_{\mathrm{P},f}}}}(\bfg\cdot\bfh\nabla\rho)_f - \mathbf{U}_f \right)\cdot\bfS_f}
{
|\bfS_f|
\frac{1}{{{A^n_{\mathrm{P},f}}}}
}
\end{equation}
上述公式与代码并不一致。在OpenFOAM中，constrainPressure为

forAll(pBf, patchi)
    {
        if (isA<fixedFluxPressureFvPatchScalarField>(pBf[patchi]))
        {
            refCast<fixedFluxPressureFvPatchScalarField>
            (
                pBf[patchi]
            ).updateCoeffs
            (
                (
                    phiHbyABf[patchi]
                  - rho.boundaryField()[patchi]
                   *MRF.relative(SfBf[patchi] & UBf[patchi], patchi)
                )
               /(magSfBf[patchi]*rhorAUBf[patchi])
            );
        }
    }

多乘了一个密度。应该改为：

forAll(pBf, patchi)
    {
        if (isA<fixedFluxPressureFvPatchScalarField>(pBf[patchi]))
        {
            refCast<fixedFluxPressureFvPatchScalarField>
            (
                pBf[patchi]
            ).updateCoeffs
            (
                (
                    phiHbyABf[patchi]
                  - MRF.relative(SfBf[patchi] & UBf[patchi], patchi)
                )
               /(magSfBf[patchi]*rhorAUBf[patchi]/rho.boundaryField()[patchi])
            );
        }
    }

在非常老的OpenFOAM版本中，看起来是正确的，与公式一致

setSnGrad<fixedFluxPressureFvPatchScalarField>
    (
        p_rgh.boundaryField(),
        (
            phiHbyA.boundaryField()
          - fvOptions.relative(mesh.Sf().boundaryField() & U.boundaryField())
        )/(mesh.magSf().boundaryField()*rAUf.boundaryField())
    );

李

activeBaffleVelocity

在我那个笔记里面有写

CFD中文网

administrators

帖子

成员列表