33 lines
1.1 KiB
C
33 lines
1.1 KiB
C
Switch Equilibrium = SLFMdict.lookup("Equilibrium");
|
|
word NstFolder = SLFMdict.lookup("NstFolder");
|
|
scalarField NstList = SLFMdict.lookup("NstList");
|
|
scalarField NstIndex(NstList.size(),0);
|
|
|
|
forAll(NstList, n)
|
|
{
|
|
NstIndex[n] = n;
|
|
}
|
|
|
|
typedef List<scalarField> scalarFieldArray1d;
|
|
typedef List<scalarFieldArray1d> scalarFieldArray2d;
|
|
typedef List<scalarFieldArray2d> scalarFieldArray3d;
|
|
|
|
SLFM::FlameStructure slfmLibrary(runTime.constant()/NstFolder, NstList, etaGrid, composition.species());
|
|
|
|
scalarFieldArray2d& Y_SLFM (slfmLibrary.Y());
|
|
scalarFieldArray2d& W_SLFM (slfmLibrary.W());
|
|
scalarFieldArray1d& T_SLFM (slfmLibrary.T());
|
|
scalarFieldArray1d& Q_SLFM (slfmLibrary.Q());
|
|
scalarFieldArray1d& rho_SLFM (slfmLibrary.rho());
|
|
scalarFieldArray1d& h_SLFM (slfmLibrary.h());
|
|
scalarFieldArray1d& Rgas_SLFM (slfmLibrary.Rgas());
|
|
|
|
for(label n=0 ; n<NstList.size() ; n++)
|
|
{
|
|
for(label i=0 ; i<Y.size() ; i++)
|
|
{
|
|
Rgas_SLFM[n] += Y_SLFM[i][n]/composition.W(i); // = 1/meanMW
|
|
}
|
|
Rgas_SLFM[n] = 8314.46*Rgas_SLFM[n]; //universal gas constant = 8314.46 [J/K.Kmole]
|
|
|
|
}
|