CFD Online Discussion Forums

CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)
-   OpenFOAM (https://www.cfd-online.com/Forums/openfoam/)
-   -   help with reactingFoam (https://www.cfd-online.com/Forums/openfoam/71494-help-reactingfoam.html)

piccinini January 5, 2010 06:46
help with reactingFoam
 
Does any one know what could be causing this error?
I'm running reactingFoam solver.

Thanks.

Code:
/*---------------------------------------------------------------------------*\
| =========                |                                                |
| \\      /  F ield        | OpenFOAM: The Open Source CFD Toolbox          |
|  \\    /  O peration    | Version:  1.6                                  |
|  \\  /    A nd          | Web:      www.OpenFOAM.org                      |
|    \\/    M anipulation  |                                                |
\*---------------------------------------------------------------------------*/
Build  : 1.6-f802ff2d6c5a
Exec  : reactingFoam
Date  : Jan 05 2010
Time  : 09:40:58
Host  : mtrs0030
PID    : 31721
Case  : /home/piccinini/work/diffflame/rans
nProcs : 1
SigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Create mesh for time = 0

Reading chemistry properties


Reading g

Reading thermophysicalProperties
Selecting psiChemistryModel ODEChemistryModel<gasThermoPhysics>
Selecting thermodynamics package hPsiMixtureThermo<reactingMixture<gasThermoPhysics>>
Selecting chemistryReader foamChemistryReader
#0  Foam::error::printStack(Foam::Ostream&) in "/cluster/apps/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libOpenFOAM.so"
#1  Foam::sigFpe::sigFpeHandler(int) in "/cluster/apps/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libOpenFOAM.so"
#2  ?? in "/lib/libc.so.6"
#3  Foam::fvPatchField<double>::operator/=(Foam::fvPatchField<double> const&) in "/cluster/apps/OpenFOAM/OpenFOAM-1.6/applications/bin/linux64GccDPOpt/reactingFoam"
#4  Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>::operator/=(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) in "/cluster/apps/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libreactionThermophysicalModels.so"
#5  Foam::multiComponentMixture<Foam::sutherlandTransport<Foam::specieThermo<Foam::janafThermo<Foam::perfectGas> > > >::correctMassFractions() in "/cluster/apps/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libreactionThermophysicalModels.so"
#6  Foam::multiComponentMixture<Foam::sutherlandTransport<Foam::specieThermo<Foam::janafThermo<Foam::perfectGas> > > >::multiComponentMixture(Foam::dictionary const&, Foam::List<Foam::word> const&, Foam::HashPtrTable<Foam::sutherlandTransport<Foam::specieThermo<Foam::janafThermo<Foam::perfectGas> > >, Foam::word, Foam::string::hash> const&, Foam::fvMesh const&) in "/cluster/apps/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libreactionThermophysicalModels.so"
#7  Foam::reactingMixture<Foam::sutherlandTransport<Foam::specieThermo<Foam::janafThermo<Foam::perfectGas> > > >::reactingMixture(Foam::dictionary const&, Foam::fvMesh const&) in "/cluster/apps/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libreactionThermophysicalModels.so"
#8  Foam::hPsiMixtureThermo<Foam::reactingMixture<Foam::sutherlandTransport<Foam::specieThermo<Foam::janafThermo<Foam::perfectGas> > > > >::hPsiMixtureThermo(Foam::fvMesh const&) in "/cluster/apps/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libreactionThermophysicalModels.so"
#9  Foam::hCombustionThermo::addfvMeshConstructorToTable<Foam::hPsiMixtureThermo<Foam::reactingMixture<Foam::sutherlandTransport<Foam::specieThermo<Foam::janafThermo<Foam::perfectGas> > > > > >::New(Foam::fvMesh const&) in "/cluster/apps/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libreactionThermophysicalModels.so"
#10  Foam::hCombustionThermo::NewType(Foam::fvMesh const&, Foam::word const&) in "/cluster/apps/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libreactionThermophysicalModels.so"
#11  Foam::psiChemistryModel::psiChemistryModel(Foam::fvMesh const&, Foam::word const&) in "/cluster/apps/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libchemistryModel.so"
#12  Foam::ODEChemistryModel<Foam::psiChemistryModel, Foam::sutherlandTransport<Foam::specieThermo<Foam::janafThermo<Foam::perfectGas> > > >::ODEChemistryModel(Foam::fvMesh const&, Foam::word const&, Foam::word const&) in "/cluster/apps/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libchemistryModel.so"
#13  Foam::psiChemistryModel::addfvMeshConstructorToTable<Foam::ODEChemistryModel<Foam::psiChemistryModel, Foam::sutherlandTransport<Foam::specieThermo<Foam::janafThermo<Foam::perfectGas> > > > >::New(Foam::fvMesh const&, Foam::word const&, Foam::word const&) in "/cluster/apps/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libchemistryModel.so"
#14  Foam::psiChemistryModel::New(Foam::fvMesh const&) in "/cluster/apps/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libchemistryModel.so"
#15  main in "/cluster/apps/OpenFOAM/OpenFOAM-1.6/applications/bin/linux64GccDPOpt/reactingFoam"
#16  __libc_start_main in "/lib/libc.so.6"
#17  _start at /usr/src/packages/BUILD/glibc-2.9/csu/../sysdeps/x86_64/elf/start.S:116
Floating exception

naveen January 6, 2010 05:15
reacting foam problem
 
hi

it is the problem of boundary conditions....

hamburgFoam January 6, 2010 17:50
hey,

I was able to reproduce a similar error when I either set all the initial mass fractions to 0 or all the mass fractions were set to 0 at any boundary.
look at the Ydefault or other spezies files you set up.

Ilja

shajitah May 4, 2010 21:04
Dear guys,

Do you know how I can download the tutorials of version 1.6 of OpenFOAM solvers e.g. reactingFoam? I mean, when you install OpenFOAM it automatically provides some tutorials for some solvers, but for some of the other solvers there's no tutorial, so I just wanna know how I can get to those ones specifically for version 1.6.

I look forward to hearing from you at your earliest convenient time.

Thanks,
Sina

naveen May 5, 2010 01:16
help with reactingFoam
 
Dear shajitah,

The reactingFoam solver is already inbuilt in the OpenFOAM 1.6 version.

Regards

Naveen

shajitah May 5, 2010 12:25
Dear Naveen,

Thanks for your quick reply. I know the solver reactingFoam exists, but I'm looking for a tutorial which is directly corresponding to reactingFoam. I have found a reactingFoam tutorial for version 1.0 on internet, but it doesn't work on version 1.6. I'm roughly sure that there would be a tutorial for most of the solvers, but can't find anything for reactingFoam!

Regards,
Sina

ddovizio March 9, 2011 13:08
1 Attachment(s)
Dear guys,
I have the same piccinini's problem.
I am using OpenFoam 1.7.1.
I have a simple methane-air non premixed combustor and used a one step mechanism for CH4 (as the counteFlowFflame2D tutorial using reactinFoam) and so I specified in my 0 directory CH4, O2, N2 and Ydefault.
In these files the mass fractions are not set all to 0.
Could you explain me what should be wrong?
I attached my 0 folder.

Regards,

Daniele

mturcios777 March 9, 2011 19:46
There is a tutorial for reactingFoam available on the wiki:
http://openfoamwiki.net/index.php/Tu..._firstTutorial

with a write up available at the following link:
http://www.tfd.chalmers.se/~hani/kur...actingFoam.pdf

The tutorial uses a single step reaction for heptane while the write-up uses a 2 step reaction mechanism for some fictitious species (they way it is done sheds light on how CHEMKIN works).

yash.aesi October 4, 2013 05:08
mass fraction of products
 
Greetings oll,
I am using OF-2.2 and tried to solve my case using reactingFoam solver with PasR model . i am getting the temperature profile almost as required but i am not getting the mass fraction of the products as desired . its coming almost double of what i required.

i checked my BC's many times and i don't think there is any mistake . i am attaching the files plz have a look and tell me then where i am doing wrong ??

Thanks in advance :)


Regards ,
Sonu


All times are GMT -4. The time now is 15:05.