CFD Online Logo CFD Online URL
Home > Forums > Software User Forums > OpenFOAM

reactingFoam - turbulent reacting flow

Register Blogs Members List Search Today's Posts Mark Forums Read

LinkBack Thread Tools Search this Thread Display Modes
Old   December 7, 2009, 13:57
Default reactingFoam - turbulent reacting flow
New Member
Ilja Sabelfeld
Join Date: Nov 2009
Posts: 22
Rep Power: 15
hamburgFoam is on a distinguished road
Hello FoamFriends,

I am a freshman at OpenFoam and trying to impliment a turbulent reacting flow. my idea was to use the reactingFoam solver.

My case:
i have zylinder with a fuel (CH4) inlet, an air inlet and a coflow inlet. at the other end of the zylinder is an outlet.

first of all i implimented a simpleFoam case just with a velocity field (without the reaction). that works fine. so, now i am trying to add the reaction.

i was trying to run the reactingFoam tutorial ( i updated the path of the "chem.inp" and "therm.dat" files as describe and run the case. after the there was a error message.

"keyword psiChemistryModel is undefined in dictionary"

so, i added the keyword as in the dieselFoam tutorial.

"psiChemistryModel ODEChemistryModel<gasThermoPhysics>;"

another error message came up.

Reading chemistry properties

Reading g

Reading thermophysicalProperties
Selecting psiChemistryModel ODEChemistryModel<gasThermoPhysics>
Selecting thermodynamics package hPsiMixtureThermo<reactingMixture<gasThermoPhysics >>
Selecting chemistryReader chemkinReader
#0 Foam::error:rintStack(Foam::Ostream&) in "/home/openfoam/OpenFOAM/OpenFOAM-1.6/lib/linuxGccDPOpt/"
#1 Foam::sigFpe::sigFpeHandler(int) in "/home/openfoam/OpenFOAM/OpenFOAM-1.6/lib/linuxGccDPOpt/"
#2 Uninterpreted:
#3 Foam:imensionedField<double, Foam::volMesh>:perator/=(Foam:imensionedField<double, Foam::volMesh> const&) in "/home/openfoam/OpenFOAM/OpenFOAM-1.6/lib/linuxGccDPOpt/"
#4 Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>:
perator/=(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) in "/home/openfoam/OpenFOAM/OpenFOAM-1.6/lib/linuxGccDPOpt/"
#5 Foam::multiComponentMixture<Foam::sutherlandTransp ort<Foam::specieThermo<Foam::janafThermo<Foam:er fectGas> > > >::correctMassFractions() in "/home/openfoam/OpenFOAM/OpenFOAM-1.6/lib/linuxGccDPOpt/"
#6 Foam::multiComponentMixture<Foam::sutherlandTransp ort<Foam::specieThermo<Foam::janafThermo<Foam:er fectGas> > > >::multiComponentMixture(Foam::dictionary const&, Foam::List<Foam::word> const&, Foam::HashPtrTable<Foam::sutherlandTransport<Foam: :specieThermo<Foam::janafThermo<Foam:erfectGas> > >, Foam::word, Foam::string::hash> const&, Foam::fvMesh const&) in "/home/openfoam/OpenFOAM/OpenFOAM-1.6/lib/linuxGccDPOpt/"
#7 Foam::reactingMixture<Foam::sutherlandTransport<Fo am::specieThermo<Foam::janafThermo<Foam:erfectGa s> > > >::reactingMixture(Foam::dictionary const&, Foam::fvMesh const&) in "/home/openfoam/OpenFOAM/OpenFOAM-1.6/lib/linuxGccDPOpt/"
#8 Foam::hPsiMixtureThermo<Foam::reactingMixture<Foam ::sutherlandTransport<Foam::specieThermo<Foam::jan afThermo<Foam:erfectGas> > > > >::hPsiMixtureThermo(Foam::fvMesh const&) in "/home/openfoam/OpenFOAM/OpenFOAM-1.6/lib/linuxGccDPOpt/"
#9 Foam::hCombustionThermo::addfvMeshConstructorToTab le<Foam::hPsiMixtureThermo<Foam::reactingMixture<F oam::sutherlandTransport<Foam::specieThermo<Foam:: janafThermo<Foam:erfectGas> > > > > >::New(Foam::fvMesh const&) in "/home/openfoam/OpenFOAM/OpenFOAM-1.6/lib/linuxGccDPOpt/"
#10 Foam::hCombustionThermo::NewType(Foam::fvMesh const&, Foam::word const&) in "/home/openfoam/OpenFOAM/OpenFOAM-1.6/lib/linuxGccDPOpt/"
#11 Foam:siChemistryModel:siChemistryModel(Foam::fvMes h const&, Foam::word const&) in "/home/openfoam/OpenFOAM/OpenFOAM-1.6/lib/linuxGccDPOpt/"
#12 Foam::ODEChemistryModel<Foam:siChemistryModel, Foam::sutherlandTransport<Foam::specieThermo<Foam: :janafThermo<Foam:erfectGas> > > >::ODEChemistryModel(Foam::fvMesh const&, Foam::word const&, Foam::word const&) in "/home/openfoam/OpenFOAM/OpenFOAM-1.6/lib/linuxGccDPOpt/"
#13 Foam:siChemistryModel::addfvMeshConstructorToTable <Foam::ODEChemistryModel<Foam:siChemistryModel, Foam::sutherlandTransport<Foam::specieThermo<Foam: :janafThermo<Foam:
erfectGas> > > > >::New(Foam::fvMesh const&, Foam::word const&, Foam::word const&) in "/home/openfoam/OpenFOAM/OpenFOAM-1.6/lib/linuxGccDPOpt/"
#14 Foam:siChemistryModel::New(Foam::fvMesh const&) in "/home/openfoam/OpenFOAM/OpenFOAM-1.6/lib/linuxGccDPOpt/"
#15 main in "/home/openfoam/OpenFOAM/OpenFOAM-1.6/applications/bin/linuxGccDPOpt/reactingFoam"
#16 __libc_start_main in "/lib/"
#17 _start at /usr/src/packages/BUILD/glibc-2.9/csu/../sysdeps/i386/elf/start.S:122

so, i don't know what that means. has anyone an helpfull advice for me how to fix my problem??? ...or can tell me how to fix the reactingFoam tutorial from the link above? ...or can share another reactingFoam tutorial with me?

Thanks in advance.


hamburgFoam is offline   Reply With Quote


Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are Off
Pingbacks are On
Refbacks are On

Similar Threads
Thread Thread Starter Forum Replies Last Post
modelling laminar and turbulent flow in the same Ol FLUENT 2 November 25, 2005 03:52
Validation case for turbulent flow Ratan Main CFD Forum 0 October 4, 2005 04:02
PISO for reacting flow newcomer_ to_ CFD Main CFD Forum 0 September 25, 2005 19:49
regarding reacting flow with coupled solver kiran FLUENT 2 April 15, 2005 05:57
Turbulent flow over flat plate Narayanan Main CFD Forum 1 May 8, 2000 08:03

All times are GMT -4. The time now is 09:23.