东岳流体对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.注释掉的两行是和等号右边的语句是一个意思吗?在编写的时候是不是哪种表述都可以?
麻烦有经验的老师进行讲解~感谢~