zhyang-dev/openfoam
1
0
Fork 0
mirror of https://develop.openfoam.com/Development/openfoam.git synced 2025-11-28 03:28:01 +00:00
Code Issues Actions Packages Projects Releases Wiki Activity

thermophysical utilities: Add environmental pressure as additional parameter

Browse Source
This commit is contained in:
Henry
2012-06-06 11:21:43 +01:00
parent b27e0dfe26
commit 18bfaf18a5
4 changed files with 9 additions and 8 deletions
Download Patch File Download Diff File Expand all files Collapse all files

3
applications/utilities/thermophysical/adiabaticFlameT/adiabaticFlameT.C
View File

@ -69,6 +69,7 @@ int main(int argc, char *argv[])
dictionary control(controlFile);
scalar P(readScalar(control.lookup("P")));
scalar T0(readScalar(control.lookup("T0")));
const word fuelName(control.lookup("fuel"));
scalar n(readScalar(control.lookup("n")));
@ -179,7 +180,7 @@ int main(int argc, char *argv[])
Info<< "products " << (1/products.nMoles())*products << ';' << endl;
scalar Tad = products.THa(reactants.Ha(T0), 1000.0);
scalar Tad = products.THa(reactants.Ha(P, T0), P, 1000.0);
Info<< "Tad = " << Tad << nl << endl;
}

9
applications/utilities/thermophysical/equilibriumCO/equilibriumCO.C
View File

@ -50,9 +50,8 @@ typedef specieThermo<janafThermo<perfectGas>, absoluteEnthalpy> thermo;
int main(int argc, char *argv[])
{
# include "setRootCase.H"
# include "createTime.H"
#include "setRootCase.H"
#include "createTime.H"
Info<< nl << "Reading Burcat data IOdictionary" << endl;
@ -71,6 +70,7 @@ int main(int argc, char *argv[])
scalar P = 1e5;
scalar T = 3000.0;
SLPtrList<thermo> EQreactions;
@ -121,7 +121,7 @@ int main(int argc, char *argv[])
forAllConstIter(SLPtrList<thermo>, EQreactions, iter)
{
Info<< "Kc(EQreactions) = " << iter().Kc(T) << endl;
Info<< "Kc(EQreactions) = " << iter().Kc(P, T) << endl;
}
Info<< nl << "end" << endl;
@ -131,4 +131,3 @@ int main(int argc, char *argv[])
// ************************************************************************* //

2
applications/utilities/thermophysical/equilibriumFlameT/equilibriumFlameT.C
View File

@ -242,7 +242,7 @@ int main(int argc, char *argv[])
scalar equilibriumFlameTemperatureNew =
products.THa(reactants.Ha(T0), adiabaticFlameTemperature);
products.THa(reactants.Ha(P, T0), P, adiabaticFlameTemperature);
if (j==0)
{

3
applications/utilities/thermophysical/mixtureAdiabaticFlameT/mixtureAdiabaticFlameT.C
View File

@ -69,6 +69,7 @@ int main(int argc, char *argv[])
dictionary control(controlFile);
scalar P(readScalar(control.lookup("P")));
scalar T0(readScalar(control.lookup("T0")));
mixture rMix(control.lookup("reactants"));
mixture pMix(control.lookup("products"));
@ -116,7 +117,7 @@ int main(int argc, char *argv[])
}
Info<< "Adiabatic flame temperature of mixture " << rMix.name() << " = "
<< products.THa(reactants.Ha(T0), 1000.0) << " K" << endl;
<< products.THa(reactants.Ha(P, T0), P, 1000.0) << " K" << endl;
return 0;
}