Navigation

    CFD中文网

    CFD中文网

    • Login
    • Search
    • 最新
    1. Home
    2. Chen_hao
    C
    • Profile
    • Following 1
    • Followers 0
    • Topics 1
    • Posts 11
    • Groups 0

    Chen_hao

    @Chen_hao

    2
    Profile views
    11
    Posts
    0
    Followers
    1
    Following
    Joined Last Online

    Chen_hao Unfollow Follow

    Latest posts made by Chen_hao

    • RE: 关于修改湍流模型中非线性雷诺应力项的问题

      @李东岳 @
      谢谢老师指导,我最后按照方法二,完全按照他的格式来修改的,修改到最后我发现一个问题,如下

       public:
      
          //typedef typename BasicTurbulenceModel::alphaField alphaField;
          //typedef typename BasicTurbulenceModel::rhoField rhoField;
          //typedef typename BasicTurbulenceModel::transportModel transportModel;
      
      
          //- Runtime type information
          TypeName("MHDkOmega");
      
      
          // Constructors
      
              //- Construct from components
              MHDkOmega
              (
                  //const alphaField& alpha,   ????
                 //const rhoField& rho,
                  const geometricOneField& alpha,  
                  const geometricOneField& rho,
                  const volVectorField& U,
                  const surfaceScalarField& alphaRhoPhi,
                  const surfaceScalarField& phi,
                  const transportModel& transport,
                  const word& propertiesName = turbulenceModel::propertiesName,
                  const word& type = typeName
              );
      

      其中 alpha是定义为geometricOneField也就是定义为1,可我的模型是三层液体,采用了multiphaseInterFoam求解器又基于VOF,我才意识到是不是湍流模型一般不能适用于多相的情况,这个alpha能不能修改来适用于这个求解器,担心自己做了无用功,希望老师能帮忙解答:xinlei:

      posted in OpenFOAM
      C
      Chen_hao
    • RE: 关于修改湍流模型中非线性雷诺应力项的问题

      @xpqiu @李东岳

      QQ图片20220701101118.png

      老师,改程序中有个关于语法的问题不太确定,ShiQuadraticKE.C 中程序没有采用模板化的,但是KOmaga.C是采用模板化的,在我更改KOmaga.H文件的继承类时,我上图这个写法对不对,直接在nonlinearEddyViscosityincompressible::RASModel 中类比eddyViscosity<RASModel<BasicTurbulenceModel>>。

      posted in OpenFOAM
      C
      Chen_hao
    • RE: 关于修改湍流模型中非线性雷诺应力项的问题

      @李东岳
      我会采用老师提供的第一种思路进行修改,感激不尽~

      posted in OpenFOAM
      C
      Chen_hao
    • RE: 关于修改湍流模型中非线性雷诺应力项的问题

      @xpqiu 谢谢老师的回答,ShihQuadraticKE模型是一个很好的参考!

      这类模型除了修改 νt ,还需要把雷诺应力的非线性部分加上。

      您所指的加上雷诺应力非线性部分就是UEqn.H中turbulence->divDevRhoReff(rho, U)项的修改对吗?

      posted in OpenFOAM
      C
      Chen_hao
    • RE: 关于修改湍流模型中非线性雷诺应力项的问题

      @李东岳
      理解了,谢谢老师的耐心回答:146:

      posted in OpenFOAM
      C
      Chen_hao
    • RE: 关于修改湍流模型中非线性雷诺应力项的问题

      @李东岳
      谢谢老师夸奖。您的意思是,只需要更换KOmega中的G,目的是求解更精准的$ν_t$,从而作用到动量方程的U求解,故不需要再修改动量方程中的雷诺应力项,就可以达到效果。不知道理解的对不对,我一直在思考既然修改的话,需要两个地方统一变换,才能达到效果,希望老师能解答这个迷惑。

      posted in OpenFOAM
      C
      Chen_hao
    • RE: 关于修改湍流模型中非线性雷诺应力项的问题

      @李东岳
      45cef8c3-7e3b-4142-888f-5fac91b8244a-IVRS%`)91~)KLI}NLFKA2DI.png

      关于雷诺应力张量形式,上面有错误,正确版如下:
      431a9ef8-7ab7-4486-909b-7e8b4d3848ad-L[OCC5NV(U5VNABTKSCMH]P.png

      不知道有没有错误,我也在写其代码,跟老师讨论。

      posted in OpenFOAM
      C
      Chen_hao
    • RE: 关于修改湍流模型中非线性雷诺应力项的问题

      @李东岳
      十分感谢老师的耐心回答!回复不及时请见谅,明白老师的意思了,方程转换好了(在编辑中稍等一会给您发过来),就您说的源项部分我单纯定义了一个标量场(语言基础差点),自己在尝试,请老师见谅。

      c87f5636-4cb4-493e-bd12-df498f8e1f46-image.png

      posted in OpenFOAM
      C
      Chen_hao
    • RE: 关于修改湍流模型中非线性雷诺应力项的问题

      @李东岳
      1d4c6455-e7c6-4745-adca-6daf38318f29-H[]@OX2Z%XH)C}(8LJMOR1Q.png

      posted in OpenFOAM
      C
      Chen_hao
    • 关于修改湍流模型中非线性雷诺应力项的问题

      @李东岳 谢谢老师对我问题的指导,理解了针对“动量方程”雷诺应力项的修改就是通过修改linearViscousStress.C。
      我先介绍下我的研究方向:模拟三层液态金属组成的电池的不稳定性,随着电压增加时三层流体会发生不稳定性即湍流,目前组里采用DNS计算消耗计算资源太大,所以我需要采用湍流模型来尝试一下这个发生过程,我采用了基于磁流体修改后的K-omega模型,其需要加入磁动能源项与各项异性的雷诺应力,以下是我需要修改的非线性雷诺应力项:
      QQ截图20220617161225.png
      我的问题是:雷诺应力也包含在kOmega.H中如下:

      tmp<volTensorField> tgradU = fvc::grad(U);
          volScalarField G  
          (
              this->GName(),
              nut*(tgradU() && dev(twoSymm(tgradU())))
          );
          tgradU.clear();
      
       tmp<fvScalarMatrix> omegaEqn	
          (
              fvm::ddt(alpha, rho, omega_) 
            + fvm::div(alphaRhoPhi, omega_)
            - fvm::laplacian(alpha*rho*DomegaEff(), omega_)
           ==
              gamma_*alpha*rho*G*omega_/k_
            - fvm::SuSp(((2.0/3.0)*gamma_)*alpha*rho*divU, omega_)
            - fvm::Sp(beta_*alpha*rho*omega_, omega_)
            + fvOptions(alpha, rho, omega_)  
            + fvm::SuSp(Comg_*S/k_, omega_) 
          );
      

      其中G中包含线性雷诺应力,针对这个非线性雷诺应力是不是我需要同时:
      1、修改动量方程中的,即通过修改linearViscousStress.C;
      2、修改kOmega模型中的,即也要修改上面的G。

      初学openfoam见识浅薄,希望各位老师和同学能不吝赐教。

      posted in OpenFOAM
      C
      Chen_hao