Simple算法中通量phi的计算

李东岳

给了速度，通量补救顺其自然的确定了么

dzw05

@东岳代码中用 phi = phiHbyA - pEqn.flux()求通量，而不是通过U求通量。在边界条件中并没有找到显式通过U修正边界通量的语句。

李东岳

phi = phiHbyA - pEqn.flux()是对phi内部面的赋值，phi的边界面给定U，phi边界面自动确定

dzw05

@东岳额，我现在的问题，就是不知道Phi的边界面在哪里自动确定的：（，在进行哪一步操作的时候确定的呢？

wwzhao

@dzw05 phi是surfaceScalarField，是没有boundary face的。

李东岳

surfaceScalarField phi
(
    IOobject
    (
        "phi",
        runTime.timeName(),
        mesh,
        IOobject::READ_IF_PRESENT,
        IOobject::AUTO_WRITE
    ),
    linearInterpolate(U) & mesh.Sf()
);

https://coding.net/u/dyfluid/p/OpenFOAM-1.6/git/blob/master/src/finiteVolume/cfdTools/incompressible/createPhi.H

dzw05

@wwzhao 其实也是有boundaryface的，在phi文件中，是会出现boundaryField的。

wwzhao

@dzw05 surfaceScalarField的boundary field总是calculated或empty，所以相当于没有。此外，surfaceScalarField的嵌套类Boundary没有initEvaluate和evaluate两个函数，所以调用surfaceScalarField的correctBoundaryConditions时会编译报错不通过。

dzw05

@wwzhao 对于boundary field这里明白了，多谢。但回到我帖子的问题，也就是说phi的boundary field没有进行显式的更新，那么以icoFoam为例，压力方程求解过程中更新phi的语句为 phi = phiHbyA - pEqn.flux()，这里求得的通量phi的边界是怎么保证正确的呢？尤其是对于wall边界条件，通量应该是零。

wwzhao

@dzw05

对流项的离散需要用到通量，internal face上的通量即为phi，而boundary face上的通量则通过U计算得到。

用PISO/SIMPLE算法得到的速度场是天然守恒的，所以boundary face上的通量肯定没有问题。

李东岳

dabiu

我同样有个问题就是通量是用phi = phiHbyA - pEqn.flux()修正，为什么不直接在下面这条语句的后面直接用phi = fvc::flux(U)来更新呢，这样直接顺带把边界也更新了？

U = HbyA - rAU*fvc::grad(p);
U.correctBoundaryConditions();

@dzw05 @李东岳 @wwzhao

qingdong_wang

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

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

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

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

CFD中文网

Simple算法中通量phi的计算