66 lines
1.4 KiB
C
66 lines
1.4 KiB
C
scalarFieldArray1d C1table //AMC C1 coefficient table
|
|
(
|
|
etamax+1,
|
|
scalarField(NVar+1, 0)
|
|
);
|
|
|
|
scalarFieldArray2d Ttable
|
|
(
|
|
etamax+1,
|
|
scalarFieldArray1d
|
|
(
|
|
NVar+1,
|
|
scalarField(NstList.size(), 0)
|
|
)
|
|
);
|
|
|
|
scalarFieldArray2d Rtable
|
|
(
|
|
etamax+1,
|
|
scalarFieldArray1d
|
|
(
|
|
NVar+1,
|
|
scalarField(NstList.size(), 0)
|
|
)
|
|
);
|
|
|
|
scalarFieldArray3d Ytable
|
|
(
|
|
postSpecieNames.size(),
|
|
scalarFieldArray2d
|
|
(
|
|
etamax+1,
|
|
scalarFieldArray1d
|
|
(
|
|
NVar+1,
|
|
scalarField(NstList.size(), 0)
|
|
)
|
|
)
|
|
);
|
|
|
|
Info<<"Construct C1 coefficient table"<<endl;
|
|
for(label j=1 ; j<etamax; j++)
|
|
{
|
|
bpdf.C1coeff(etaValue[j], varValue, C1table[j]);
|
|
}
|
|
|
|
Info<<"Construct Temp and Rgas table"<<endl;
|
|
for(label j=0 ; j<=etamax ; j++)
|
|
{
|
|
scalar mf = etaValue[j];
|
|
for(label v=0 ; v<=NVar ; v++)
|
|
{
|
|
scalar mfVar = varValue[v]*mf*(1.0-mf);
|
|
for(label n=0 ; n<NstList.size() ; n++)
|
|
{
|
|
Ttable[j][v][n] = bpdf.evaluate(mf, mfVar, etaValue, Y_SLFM[n][Ysize]);
|
|
Rtable[j][v][n] = bpdf.evaluate(mf, mfVar, etaValue, Rgas_SLFM[n]);
|
|
forAll(postSpecieIndices, yi)
|
|
{
|
|
const label y = postSpecieIndices[yi];
|
|
Ytable[yi][j][v][n] = bpdf.evaluate(mf, mfVar, etaValue, Y_SLFM[n][y]);
|
|
}
|
|
}
|
|
}
|
|
Info<<j/etamax*100<<"%"<<endl;
|
|
}
|