H2 injection temperature as IC under janaf limit (modifying janafThermoI.H)
Hi everybody,
I'm newby in OpenFOAM and in C++ programming language and I try to simulate combustion air-H2 inside a scramjet combustion chamber. For my thesis I need to set a injection temperature of hydrogen around 133K but it is already out of temperature range. I have calculated coefficients for the lower temperature using Chemkin Fitdata function, but I don't know how to modify OpenFoam to read this new file (i.e therm_lowrange.dat). Looking around the code I found the file where temperature is checked, janafThermoI.H: template<class equationOfState> inline void Foam::janafThermo<equationOfState>::checkT(const scalar T) const { if (T < Tlow_ || T > Thigh_) { FatalErrorIn ( "janafThermo<equationOfState>::checkT(const scalar T) const" ) << "attempt to use janafThermo<equationOfState>" " out of temperature range " << Tlow_ << " -> " << Thigh_ << "; T = " << T << abort(FatalError); } } I would like to bypass this check for example say OpenFoam to read something like that from thermophysicalProperties file: CHEMKINThermoFile_lowrange "$FOAM_CASE/chemkin/therm_lowrange.dat" before abort (if doesn't exist). What is the easy way to do it? Do you know if foamChemistryReader can do this? Thanks for any suggestions! Luigi |
This is my new janafThemoI:
Code:
template<class equationOfState> Code:
Solving chemistry |
All times are GMT -4. The time now is 13:06. |