compressibleInterFoam alphaEqn.H

hongjiewang

东岳流体对compressibleInterFoam解析中，对dgdt的规定为
\begin{equation}
\alpha_{1}\alpha_{2}(\frac{1}{\rho _{2}}\frac{D\rho _{2}}{Dt}-\frac{1}{\rho _{1}}\frac{D\rho _{1}}{Dt})
\end{equation}
http://dyfluid.com/compressibleInterFoam.html
但是在code中，

         fvScalarMatrix alpha1Eqn
         (
             (
                 LTS
               ? fv::localEulerDdtScheme<scalar>(mesh).fvmDdt(alpha1)
               : fv::EulerDdtScheme<scalar>(mesh).fvmDdt(alpha1)
             )
           + fv::gaussConvectionScheme<scalar>
             (
                 mesh,
                 phiCN,
                 upwind<scalar>(mesh, phiCN)
             ).fvmDiv(phiCN, alpha1)
        // - fvm::Sp(fvc::ddt(dimensionedScalar(dimless, 1), mesh)
        //           + fvc::div(phiCN), alpha1)
          ==
             Su + fvm::Sp(Sp + divU, alpha1)
         );

 forAll(dgdt, celli)
 {
     if (dgdt[celli] > 0.0)
     {
         Sp[celli] -= dgdt[celli]/max(1.0 - alpha1[celli], 1e-4);
         Su[celli] += dgdt[celli]/max(1.0 - alpha1[celli], 1e-4);
     }
     else if (dgdt[celli] < 0.0)
     {
         Sp[celli] += dgdt[celli]/max(alpha1[celli], 1e-4);
     }
 }

volScalarField dgdt(alpha1*fvc::div(phi));

有两个问题想要请教：
1.dgdt在createFields.H中的定义为什么和解析中的不同，还是说后续有对dgdt进行计算~
2.注释掉的两行是和等号右边的语句是一个意思吗？在编写的时候是不是哪种表述都可以？
麻烦有经验的老师进行讲解~感谢~

ShaneHEEE

creatField.H中的相当于初始化，后续会通过求解方程更新cell中的值

2.“==”就相当于“-”，注释掉的两行和等号右边是一个意思

李东岳

@hongjiewang 后续dgdt有计算，我那个没写，工作量太大了。我这面有个英文sci投出去写得很详细，但是一审被居高。也没空翻译成中文

hongjiewang

@shaneheee 嗯嗯好滴谢谢老师的回复

hongjiewang

@李东岳谢谢老师~希望以后有机会能读到老师对这部分的分析

CFD中文网

compressibleInterFoam alphaEqn.H