Skip to content
  • 最新
  • 版块
  • 东岳流体
  • 随机看[请狂点我]
皮肤
  • Light
  • Cerulean
  • Cosmo
  • Flatly
  • Journal
  • Litera
  • Lumen
  • Lux
  • Materia
  • Minty
  • Morph
  • Pulse
  • Sandstone
  • Simplex
  • Sketchy
  • Spacelab
  • United
  • Yeti
  • Zephyr
  • Dark
  • Cyborg
  • Darkly
  • Quartz
  • Slate
  • Solar
  • Superhero
  • Vapor

  • 默认(不使用皮肤)
  • 不使用皮肤
折叠
CFD中文网

CFD中文网

C

cavitazione.club

@cavitazione.club
关于
帖子
2
主题
2
群组
0
粉丝
0
关注
0

帖子

最新

  • 粘度场异常
    C cavitazione.club
    1. 速度场
      U.png
    2. 压力场
      p.png
    3. 粘度场
      nut.png

    粘度边界:进出口calculated,机翼表面nutUSpaldingWallFunction,上下边界symmetry简化问题

    计算收敛,domainminandmax监测稳定无震荡

    求解器:default interPhaseChangeFoam,ofv1912

    湍流模型:kEqn LES

    粘度场异常原因?


  • Rhie-Chow interplation
    C cavitazione.club

    Original RC interplation
    No time dependency, No relaxation of the velocity field:

    // Rhie-Chow interplation
    phi = (U_avg_f & mesh.Sf()) 
        - ( (DUf*( gradp_f - gradp_avg_f)) & mesh.Sf() );
    

    Time dependency, relaxation of the velocity field:

    // Rhie-Chow interplation
    phi = (U_avg_f & mesh.Sf()) 
        - ( (DUf*( gradp_f - gradp_avg_f)) & mesh.Sf() )
        + (scalar(1) - URFU)*(phi.prevIter() - (U_avg_prevIter_f & mesh.Sf()))
        + DTf*( phi_old - (U_old_f& mesh.Sf()));
    

    OpenFOAM
    No time dependency, with or without relaxation of the velocity field (check simpleFoam solver):

    surfaceScalarField phiHbyA("phiHbyA", fvc::interpolate(HbyA) & mesh.Sf());
    

    Questions: Where is the second term of the original RC? Where is the third term of the original RC if taking the relaxation of the velocity field into consideration?

    Time dependency, relaxation of the velocity field (check pisoFoam solver):

    surfaceScalarField phiHbyA
     (
         "phiHbyA",
         (fvc::interpolate(HbyA) & mesh.Sf())
       + fvc::ddtPhiCorr(rAU, U, phi)
    );
    

    Questions: ddtPhiCorr takes the function of (the third term+the fourth term of the original RC)? I guess NO. Let's recall the definition of ddtPhiCorr:

    tmp<fluxFieldType> ddtPhiCorr
    (
        new fluxFieldType
        (
            ddtIOobject,
            rDeltaT
          * this->fvcDdtPhiCoeff(U.oldTime(), phiAbs.oldTime())
          * (
                 fvc::interpolate(rA*rho.oldTime())*phiAbs.oldTime()
               - (
                        fvc::interpolate(rA*rho.oldTime()*U.oldTime())
                      & mesh().Sf()
                  )
             )
        )
    );
    

    Look like ddtPhiCorr returns the following codes:

    fvcDdtPhiCoeff*DTf*( phi_old - (U_old_f& mesh.Sf()));
    

    Last question: How to take the effect of relaxation of the velocity field into ddtPhiCorr?
    谢谢

  • 登录

  • 登录或注册以进行搜索。
  • 第一个帖子
    最后一个帖子
0
  • 最新
  • 版块
  • 东岳流体
  • 随机看[请狂点我]