1.reactingFoam是否将chemistryProperties文件中的(chemistry on)改为off,并且热物理模型如下设置thermoType
{
type hePsiThermo;
mixture multiComponentMixture;
transport sutherland;
thermo janaf;
energy sensibleEnthalpy;
equationOfState perfectGas;
specie specie;
},就可以求解多组分气体扩散问题呢?
2.观察到chtMultiRegionFoam求解器中,关于流体域的求解代码中也求解了组分方程,那么是不是这个chtMultiRegionFoam求解器也完全适合流体域中多组分扩散的求解呢?以下是求解器源代码中含的YEqn.H文件```
code_text
{
mvConvection = tmp<fv::convectionScheme<scalar>>
(
fv::convectionScheme<scalar>::New
(
mesh,
fields,
phi,
mesh.schemes().div("div(phi,Yi_h)")
)
);
}
reaction.correct();
forAll(Y, i)
{
if (composition.solve(i))
{
volScalarField& Yi = Y[i];
fvScalarMatrix YiEqn
(
fvm::ddt(rho, Yi) + mvConvection->fvmDiv(phi, Yi)
+ thermophysicalTransport.divj(Yi)
==
reaction.R(Yi)
+ fvModels.source(rho, Yi)
);
YiEqn.relax();
fvConstraints.constrain(YiEqn);
YiEqn.solve("Yi");
fvConstraints.constrain(Yi);
}
}
composition.normalise();