icoFoam求解器解析中问题。

OF九月九

大家好，在icoFoam解析中
···

        volScalarField rAU(1.0/UEqn.A());
        volVectorField HbyA("HbyA", U);
        HbyA = rAU*UEqn.H();//公式(15)

···UEqn.A()、UEqn.H()是什么意思？

CFD中文网

@OF九月九在 [icoFoam求解器解析中问题。]

\begin{equation}
UEqn.A() = \frac{1}{{{A_\mathrm{p}}}}
\end{equation}

\begin{equation}
UEqn.H() = \left( { - \sum {{A_\mathrm{n}}u_n^r} + E_\mathrm{p}^r} \right)
\end{equation}

Clear?

:sunglasses:

OF九月九

我知道是这个，我以为.A()和.H()包含什么重要东西呢。谢了！:happy:

张某人

@cfd中文网你好，不过似乎UEqn.A()应该等于Ap吧。按你给出的定义，rAU*UEqn.H()就等于Ap*UEqn.H()了，这和HbyA的定义不一致了。

东岳

是的，你说的正确，UEqn.A()等于Ap。

maoyanjun_dut

if (nonOrth == nNonOrthCorr)
{
phi = phiHbyA - pEqn.flux();
//使用求解的压力场修正通量场，在最后一次修正的时候通量守恒，Issa指出，大约需要2-3次内循环步。
//对应方程(26)，pEqn.flux()返回方程(26)方程右边第二项，也为fvc::interpolate(rUA)*fvc::snGrad(p)*mag(mesh.Sf())。
//某些可压缩求解器其中的pEqn.flux()可能为+号，即为phi = phiHbyA + pEqn.flux()。这是因为pEqn中的laplacian项为−号
}
@东岳做个勘误，注释中的 fvc::interpolate(rUA)*fvc::snGrad§*mag(mesh.Sf()) 似乎应该是 rAU

还行麻烦解释一下

surfaceScalarField phiHbyA
            (
                "phiHbyA",
                (fvc::interpolate(HbyA) & mesh.Sf())
                //此处依据Rhie-chow插值原理，HbyA使用线性插值得到，
                //即需要在算例中设定interpolate(HbyA)的格式
              + fvc::interpolate(rAU)*fvc::ddtCorr(U, phi)
            );

里面的 fvc::interpolate(rAU)*fvc::ddtCorr(U, phi) 的意思。谢谢啦。

张某人

@东岳李老师，您好。个人认为icoFoam解析中公式(22)->(23)推导过程中，(22)式等号右侧高斯积分中的1/Ap,f在(23)式中应该变成散度项里的1/Ap。

东岳

@张某人多谢，已更新
@maoyanjun_dut 有空我研究研究