[AleDR]

Hi guys,

I'm really new in OpenFOAM, but I have some little experience in combustion modelling with Fluent and burner laboratory testing.
I really appreciate the efforts of the OF community, so I would make some contributions if possible!

Reading through some posts here, I found that combustion modelling (a hard topic for any CFD engineer!) is somewhat not well established in OpenFOAM. But this is the main objective for the community, collaborating for common knowledge... So why don't start a sub-topic in this forum? Hope that I will find others interested in it!


Anyway before starting I have a question:
- Is radiation modelling possible for combustion equipment simulation in OpenFOAM ?

I think about standard workflow in Fluent:
Cold flow solution -> Reacting flow solution -> Radiation model activated --> Solution!
I have not understood if this is possible in OF too!


Another topic is the stability of combustion model in reactingFoam...
I have gone through the tutorial "Tut reactingFoam firstTutorial" in OpenFOAM Wiki adapting it to 1.7.1 version and I get this error:

--> FOAM FATAL ERROR:
attempt to use janafThermo<equationOfState> out of temperature range 200 -> 5000; T = 5000.01

From function janafThermo<equationOfState>::checkT(const scalar T) const
in file /home/opencfd/OpenFOAM/OpenFOAM-1.7.x/src/thermophysicalModels/specie/lnInclude/janafThermoI.H at line 63.

FOAM aborting

#0 Foam::error:: PrintStack(Foam::Ostream&) in "/opt/openfoam171/lib/linuxGccDPOpt/libOpenFOAM.so"
#1 Foam::error::abort() in "/opt/openfoam171/lib/linuxGccDPOpt/libOpenFOAM.so"
#2 Foam::specieThermo<Foam::janafThermo<Foam: :PerfectGas> >::H(double) const in "/opt/openfoam171/lib/linuxGccDPOpt/libreactionThermophysicalModels.so"
#3 Foam::ODEChemistryModel<Foam: :PsiChemistryModel, Foam::sutherlandTransport<Foam::specieThermo<Foam: :janafThermo<Foam: :PerfectGas> > > >::solve(double, double) in "/opt/openfoam171/lib/linuxGccDPOpt/libchemistryModel.so"
#4
in "/opt/openfoam171/applications/bin/linuxGccDPOpt/reactingFoam"
#5 __libc_start_main in "/lib/tls/i686/cmov/libc.so.6"
#6
in "/opt/openfoam171/applications/bin/linuxGccDPOpt/reactingFoam"

It seems to me (please, tell me if I'm wrong...) that the solver does not limit the temperature rise during iterations. In Fluent this is done automatically and you get a warning for this, but the solver continues to iterate and often this problem is overcome as the solution starts to settle down. The maxTemperature is a setting of the software and can be modified by the user. Maybe this should be implemented in OpenFOAM too, but I have to be a little skilled with reactingFoam before!!!

In other words, the Tut reactionFoam firstTutorial does not work for me!!!

Here are the files if you want to have a deeper look into my modifications (actually the tutorial files I have downloaded featured a controlDict with a wrong line for the solver... dieselFoam instead of reactingFoam ! I made some other minor modifications too and I think it should be necessary to add a good and physically sound /0 too...).
After that we could -somehow!- correct the tutorial on the OpenFOAM Wiki, in order to give our contribution and build a good tutorial section for the entire community (it seems to me very poor!!!).

I'll be waiting for collaborations, thank you!!!

.Alex.