@李东岳 就算仿照sprayfoam进行添加的，createFields.H里还添加了下面部分，twoPhaseMixtureThermo里将RhoThermo改为了rhoReactionThermo，应该就这些。
rhoReactionThermo& thermo2 = mixture.thermo2();
basicSpecieMixture& composition = thermo2.composition();
PtrList<volScalarField>& Y = composition.Y();
const word inertSpecie(thermo2.lookup("inertSpecie"));
if (!composition.species().found(inertSpecie))
{
FatalIOErrorIn(args.executable().c_str(), thermo2)
<< "Inert specie " << inertSpecie << " not found in available species "
<< composition.species()
<< exit(FatalIOError);
}
const label inertIndex(composition.species()[inertSpecie]);
multivariateSurfaceInterpolationScheme<scalar>::fieldTable fields;
forAll(Y, i)
{
fields.add(Y[i]);
}
surfaceScalarField alpharho2
(
IOobject
(
"alphaRho2Phi",
runTime.timeName(),
mesh,
IOobject::READ_IF_PRESENT,
IOobject::AUTO_WRITE
),
alpha2*rho2
);
surfaceScalarField alphaRho2Phi
(
IOobject
(
"alphaRho2Phi",
runTime.timeName(),
mesh,
IOobject::READ_IF_PRESENT,
IOobject::AUTO_WRITE
),
fvc::interpolate(alpha2*rho2)*phi
);
