这段代码是YEqn方程的第一相的代码,我的算例中第一相有五种组分,现在我想针对这五种组分分别改写他们对应的输运方程,应该怎么去修改?
forALL(Y1, i)
{
tmp<fvScalarMatrix>Y1iEqn(phase1.YiEqn(Y1[i]));
if (Y1iEqn.valid())
Y1iEqn =
(
Y1iEqn
==
*massTransfer[Y1[i].name()]
+ fvOptions(alpha1, rho1, Y1[i])
);
Y1iEqn->relax();
Y1iEqn->solve(mesh.solver(Yi))
}