fvVectorMatrix Ueqn



  • fvVectorMatrix UEqn
        (
            fvm::ddt(alphacRho, U)
          + MRF.DDt(alphacRho, U)
          + fvm::ddt(rho, U) + fvm::div(rhoPhi, U)
          + MRF.DDt(rho, U)
          + turbulence->divDevRhoReff(rho, U)
         ==
            fvOptions(rho, U)
              + cloudSU
        );
    
    
    

    MPPICInterfoam 中的UEqn == 后面是源项。 但是我不太明白的是, 为啥这个matrix中间有==呢?因为简单的icoFoam中没有后面的源项,就会感觉比较容易理解,但是加上了等号,我就不能理解了。
    还有就是这个cloudSU也是源项吗还是ö粒子力?这个方程在后面不是已经加上particle forces了吗?

    UEqn
                     ==
                            fvc::reconstruct
                            (
                                    phicForces/rAUcf // is added as particle force
                              +
                                    (
                                            fvc::interpolate
                                            (
                                            mixture.sigmaK()
                                            )*fvc::snGrad(alpha1) // can be supplanted by the mixture.surfaceTensionForce()
                                            - ghf*fvc::snGrad(rho)
                                            - fvc::snGrad(p_rgh)
                                    ) * mesh.magSf()
                )
    
    

  • 管理员

    矩阵求解是Ax=b,等号后面的是b
    源项处理一部分隐性,一部分显性,phicForce是显性的部分



  • 矩阵是Ax = b, 然后solve(UEqn = fvc:xxx)是求解的Ax=b? 那么这样的话A 是前面的UEqn (fvVectorMatrix)?那这个里面的matrix为啥还有等号呢



  • 然后还有一点我不太明白的地方, 就是这个地方的这个地方的粒子力为啥先乘了rAUc然后又除了rAUc?好像是关于离散吗?

    volScalarField rAUc(1.0/UEqn.A());
            surfaceScalarField rAUcf(fvc::interpolate(rAUc));
    
    
            surfaceScalarField phicForces
            (
                    (fvc::interpolate(rAUc*cloudVolSUSu) & mesh.Sf())
            );
    
        if (pimple.momentumPredictor())
        {
            solve
            (
                            UEqn
                     ==
                            fvc::reconstruct
                            (
                                    phicForces/rAUcf // is added as particle force
    + ....
    
    

  • 管理员

    pEqn里面用到了rAUc


Log in to reply
 


CFD中文网 | 东岳流体学术 | 东岳流体商业 | 吉ICP备20003622号-1