多组分气体扩散问题求解

13562349205

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();

ShaneHEEE

可以的

CFD中文网

多组分气体扩散问题求解