一些constrainPressure的疑问



  • 先贴icofoam中的一段代码

                adjustPhi(phiHbyA, U, p);
    
                // Update the pressure BCs to ensure flux consistency
                constrainPressure(p, U, phiHbyA, rAU);
    

    疑问是constrainPressure是做什么用的?有的求解器前面跟了adjustPhi,像icofoam,adjustPhi是修正压力边界使连续性方程有解,有的求解器像buoyantBoussinesqPimpleFoam没有adjustPhi只有constrainPressure。
    再来看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])
                );
            }
        }
    }
    

    所以以上代码的意思是压力的边界条件如果是fixedFluxPressure,那么压力的边界面法向梯度等于

    (边界流量 - 边界流量)/网格面积*密度/Uequ.A() ????????这不成0了

    完全一脸茫然??????这干嘛用的,有啥意义呢?



  • 。。。。。。自问自答 ,汗

    phiHbyABf不是边界流量。。。。

    所以以下代码应该是符合动量方程,猜测可能是进入压力PISO循环之前修改下边界条件的系数

                    (
                        phiHbyABf[patchi]
                      - rho.boundaryField()[patchi]
                       *MRF.relative(SfBf[patchi] & UBf[patchi], patchi)
                    )
                   /(magSfBf[patchi]*rhorAUBf[patchi])
    

Log in to reply
 

CFD中文网 2016 - 2020 | 京ICP备15017992号-2