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

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

@xpqiu @李东岳

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

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

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

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

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

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

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

@李东岳

关于雷诺应力张量形式，上面有错误，正确版如下：

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

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

@李东岳

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

我的问题是：雷诺应力也包含在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见识浅薄，希望各位老师和同学能不吝赐教。