[OMN]

Hello

I am writing a new thermo class to model the mixture of two gases, for example air (itself a mixture) and methane.

I think I have understood how to define it, but I have trouble understanding how to link the scalar used to transport the methane specie, called Y in the 0 directory, to the properties of the methane specie defined in the thermophysicalProperties file.

What I have done is the following:
1) Added the following
makeBasicMixture
(
multiComponentMixture,
constTransport,
hConstThermo,
perfectGas
);
to the file basicMixtures.C

2) Added the following
Y_
(
IOobject
(
“Y”,
mesh.time().timeName(),
mesh,
IOobject::MUST_READ,
IOobject::AUTO_WRITE
),
mesh
),
to the file basicThermo.C

3) Created a new perfectGas directory, called speciesPerfectGas, and changed the calculation of psi to the following:
inline scalar speciesPerfectGas:si(scalar, scalar T, scalar Y) const
{
return 1.0/((R()*(1-Y) + R_Y*Y)*T);
}


I have tried to find where R() is calculated, to know how to define R_Y in the previous code, but I can't find it. Would someone know that please?

Is there an example of multiComponentMixture somewhere?

Cheers

Olivier

[Chris Lucas]

Hi,

have a look at the class specie in the folder "thermophysicalModels/specie/specie"

Regards,
Christian

[OMN]

Thank you Christian!

Cheers

Olivier

[carowjp]

Hello Olivier,

I am working on a similar feature. Would you be willing to share your implementation?

thank you and best regards,

James

[OMN]

Yes of course James.

Have you got an email?

regards

Olivier

[carowjp]

My email is:

james dot p dot carow at wmich dot edu

Thank you Olivier and have a good weekend!

James