CFD中文网

    CFD中文网

    • 登录
    • 搜索
    • 最新

    fvVectorMatrix Ueqn

    OpenFOAM
    2
    5
    970
    正在加载更多帖子
    • 从旧到新
    • 从新到旧
    • 最多赞同
    回复
    • 在新帖中回复
    登录后回复
    此主题已被删除。只有拥有主题管理权限的用户可以查看。
    • Q
      qingdong_wang 最后由 编辑

      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()
                  )
      
      
      1 条回复 最后回复 回复 引用
      • 李东岳
        李东岳 管理员 最后由 编辑

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

        CFD高性能服务器 http://dyfluid.com/servers.html
        2023年,线下CFD课,预热一下 http://dyfluid.com/class.html

        1 条回复 最后回复 回复 引用
        • Q
          qingdong_wang 最后由 编辑

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

          1 条回复 最后回复 回复 引用
          • Q
            qingdong_wang 最后由 编辑

            然后还有一点我不太明白的地方, 就是这个地方的这个地方的粒子力为啥先乘了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
            + ....
            
            
            1 条回复 最后回复 回复 引用
            • 李东岳
              李东岳 管理员 最后由 编辑

              pEqn里面用到了rAUc

              CFD高性能服务器 http://dyfluid.com/servers.html
              2023年,线下CFD课,预热一下 http://dyfluid.com/class.html

              1 条回复 最后回复 回复 引用
              • First post
                Last post