72 int main(
int argc,
char *argv[])
81 IFstream controlFile(controlFileName);
84 if (!controlFile.good())
87 <<
"Cannot read file " << controlFileName
95 mixture rMix(control.lookup(
"reactants"));
96 mixture pMix(control.lookup(
"products"));
99 Info<<
nl <<
"Reading Burcat data dictionary" <<
endl;
104 IFstream BurcatCpDataFile(BurcatCpDataFileName);
107 if (!BurcatCpDataFile.good())
110 <<
"Cannot read file " << BurcatCpDataFileName
119 rMix[0].volFrac()*
thermo(CpData.lookup(rMix[0].name()))
122 for (label i = 1; i < rMix.size(); i++)
124 reactants = reactants
125 + rMix[i].volFrac()*
thermo(CpData.lookup(rMix[i].name()));
131 2*pMix[0].volFrac()*
thermo(CpData.lookup(pMix[0].name()))
134 for (label i = 1; i < pMix.size(); i++)
137 + 2*pMix[i].volFrac()*
thermo(CpData.lookup(pMix[i].name()));
140 Info <<
"Adiabatic flame temperature of mixture " << rMix.name() <<
" = "
141 << products.TH(reactants.H(T0), 1000.0) <<
" K" <<
endl;