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

Chen_hao

@李东岳
理解了，谢谢老师的耐心回答

xpqiu

@chen_hao
你可以看看 src/TurbulenceModels/incompressible/turbulentTransportModels/RAS/ShihQuadraticKE 这个湍流模型
这是一个非线性的湍流模型而且也只考虑了二阶非线性项，跟你要实现的那个类似。
这类模型除了修改 $\nu_t$ ，还需要把雷诺应力的非线性部分加上。

李东岳

对。

你要是动kOmega模型，可以把$G$换成你的$\tau$，在linearStress.C那面还要处理一下devTau函数，要把非线性的贡献加上（$\tau$减去线性部分）。
你要是动非线性类模型如ShihQuadraticKE，把$G$换成你的$\tau$。这样不需要特为处理devTau函数了。但要处理下omega方程（如果没有的话）。

Chen_hao

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

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

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

xpqiu

@chen_hao 在关于修改湍流模型中非线性雷诺应力项的问题中说：

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

不是，不需要修改这个函数，你仔细看看 ShihQuadraticKE 这个模型的继承关系，它继承的是 nonLinearEddyViscousity 类。这个类计算雷诺应力的方式不一样，除了线性项还有一个 nonLinearStress 项。你只要按照这个框架去写，把你的非线性项赋值给 nonLinearStress 就行了。

Chen_hao

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

Chen_hao

@xpqiu @李东岳

QQ图片20220701101118.png

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

李东岳

我觉得你不用这么写。你动继承肯定会有一些问题，目前也看不出来，不是很strightworrd。比较简单的是：

kOmega.H头文件不用动，你先去把tau植入到你的kOmega.C里面
计算G
去devTau函数里面把tau的非线性部分加上

你要是按照ShiQuadraticKE套路来，

把k方程改成你的omega方程，不用动.H文件
非线性项按照你的tau去写

第一种简单一些，但有一些硬植入。就是devTau函数这块
第二种更灵活。做完了就是非线性komega模型。

Chen_hao

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

 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能不能修改来适用于这个求解器，担心自己做了无用功，希望老师能帮忙解答

李东岳

从这个alpha来看，确实不能用于多相

Chen_hao

@李东岳

老师，像这种形式的alpha的定义，能用于多相的情况吗？

typedef typename BasicTurbulenceModel::alphaField alphaField
const alphaField& alpha

BasicTurbulenceModel这个通用类是代指的是KomegaSST？（就本例）

Chen_hao

@李东岳
请老师看看G这部分代码可不可以这样写，
矢量场Snl有点搞不清=和（）赋值，tgradU后面用不用加括号~

tmp<volTensorField> tgradU = fvc::grad(U);//tmp reload to Field
    volScalarField S2(2*magSqr(symm(tgradU())));
    
    tensor I(1,0,0,0,1,0,0,0,1);  
    volSymmTensorField S(twoSymm(tgradU));
    volTensorField W(twoSkew(tgradU));
    volTensorField Snl =
        1/omega_
       *(
           S&S-(1/3)*magSqr(W)*I
         + twoSymm(S&W)
         + (2.5)*(W&W-(1/3)*magSqr(W)*I)
        );
    
    volScalarField::Internal GbyNu(dev(twoSymm(tgradU()())+Snl) && tgradU()());
    volScalarField::Internal G(this->GName(), nut()*GbyNu);
    tgradU.clear();

李东岳

可以。我记得湍流方程右边都是只考虑网格点的作用，不需要考虑边界。因此你的G声明了内部场，没有问题。

CFD中文网

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