ReactingFoam with dynamic mesh
 

Hi everyone,

I am trying to implement reactingfoam (OpenFoam-2.2.1) for dynamic mesh (layered engineMesh Type) i.e engines.. I tried to follow the same steps as it has been done to incorporate XiFoam in engines to create engineFoam.
The solver complies fine without any error.but When I try running it i get the following error:


Reading g
Creating reaction model

Selecting combustion model PaSR<psiChemistryCombustion>
Selecting chemistry type
{
chemistrySolver ode;
chemistryThermo psi;
}

Selecting thermodynamics package
{
type hePsiThermo;
mixture reactingMixture;
transport sutherland;
thermo janaf;
energy sensibleEnthalpy;
equationOfState perfectGas;
specie specie;
}

Selecting chemistryReader foamChemistryReader

#0 Foam-error-printStack(Foam-Ostream&) in "/usr/apps1/openfoam-2.2.1/OpenFOAM-2.2.1/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#1 Foam-sigFpe-sigHandler(int) in "/usr/apps1/openfoam-2.2.1/OpenFOAM-2.2.1/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#2 __restore_rt at sigaction.c:0
#3 Foam:heThermo<Foam: psiReactionThermo, Foam-SpecieMixture<Foam-reactingMixture<Foam-sutherlandTransport<Foam-species-thermo<Foam-janafThermo<Foam-perfectGas<Foam-specie> >, Foam-sensibleEnthalpy> > > > >-he(Foam-Field<double> const&, Foam-Field<double> const&, int) const in "/usr/apps1/openfoam-2.2.1/OpenFOAM-2.2.1/platforms/linux64GccDPOpt/lib/libreactionThermophysicalModels.so"
#4 Foam:heThermo<Foam: psiReactionThermo, Foam-SpecieMixture<Foam-reactingMixture<Foam-sutherlandTransport<Foam-species-thermo<Foam-janafThermo<Foam-perfectGas<Foam-specie> >, Foam-sensibleEnthalpy> > > > >-init() in "/usr/apps1/openfoam-2.2.1/OpenFOAM-2.2.1/platforms/linux64GccDPOpt/lib/libreactionThermophysicalModels.so"
#5 Foam:heThermo<Foam: psiReactionThermo, Foam-SpecieMixture<Foam-reactingMixture<Foam-sutherlandTransport<Foam-species-thermo<Foam-janafThermo<Foam-perfectGas<Foam-specie> >, Foam-sensibleEnthalpy> > > > >-heThermo(Foam-fvMesh const&, Foam-word const&) in "/usr/apps1/openfoam-2.2.1/OpenFOAM-2.2.1/platforms/linux64GccDPOpt/lib/libreactionThermophysicalModels.so"
#6 Foam: psiReactionThermo-addfvMeshConstructorToTable<Foam-hePsiThermo<Foam-psiReactionThermo, Foam-SpecieMixture<Foam-reactingMixture<Foam-sutherlandTransport<Foam-species-thermo<Foam-janafThermo<Foam-perfectGas<Foam-specie> >, Foam-sensibleEnthalpy> > > > > >-New(Foam-fvMesh const&, Foam-word const&) in "/usr/apps1/openfoam-2.2.1/OpenFOAM-2.2.1/platforms/linux64GccDPOpt/lib/libreactionThermophysicalModels.so"
#7 Foam-autoPtr<Foam-psiReactionThermo> Foam-basicThermo-New<Foam-psiReactionThermo>(Foam-fvMesh const&, Foam-word const&) in "/usr/apps1/openfoam-2.2.1/OpenFOAM-2.2.1/platforms/linux64GccDPOpt/lib/libreactionThermophysicalModels.so"
#8 Foam: psiReactionThermo-New(Foam-fvMesh const&, Foam-word const&) in "/usr/apps1/openfoam-2.2.1/OpenFOAM-2.2.1/platforms/linux64GccDPOpt/lib/libreactionThermophysicalModels.so"
#9 Foam- psiChemistryModel-psiChemistryModel(Foam-fvMesh const&) in "/usr/apps1/openfoam-2.2.1/OpenFOAM-2.2.1/platforms/linux64GccDPOpt/lib/libchemistryModel.so"
#10 Foam-chemistryModel<Foam-psiChemistryModel, Foam-sutherlandTransport<Foam-species-thermo<Foam-janafThermo<Foam-perfectGas<Foam-specie> >, Foam-sensibleEnthalpy> > >-chemistryModel(Foam-fvMesh const&) in "/usr/apps1/openfoam-2.2.1/OpenFOAM-2.2.1/platforms/linux64GccDPOpt/lib/libchemistryModel.so"
#11 Foam-ode<Foam-chemistryModel<Foam-psiChemistryModel, Foam-sutherlandTransport<Foam-species-thermo<Foam-janafThermo<Foam-perfectGas<Foam-specie> >, Foam-sensibleEnthalpy> > > >-ode(Foam-fvMesh const&) in "/usr/apps1/openfoam-2.2.1/OpenFOAM-2.2.1/platforms/linux64GccDPOpt/lib/libchemistryModel.so"
#12 Foam- psiChemistryModel-addfvMeshConstructorToTable<Foam-ode<Foam-chemistryModel<Foam-psiChemistryModel, Foam-sutherlandTransport<Foam-species-thermo<Foam-janafThermo<Foam-perfectGas<Foam-specie> >, Foam-sensibleEnthalpy> > > > >-New(Foam-fvMesh const&) in "/usr/apps1/openfoam-2.2.1/OpenFOAM-2.2.1/platforms/linux64GccDPOpt/lib/libchemistryModel.so"
#13 Foam-autoPtr<Foam-psiChemistryModel> Foam-basicChemistryModel-New<Foam-psiChemistryModel>(Foam-fvMesh const&) in "/usr/apps1/openfoam-2.2.1/OpenFOAM-2.2.1/platforms/linux64GccDPOpt/lib/libchemistryModel.so"
#14 Foam- psiChemistryModel-New(Foam-fvMesh const&) in "/usr/apps1/openfoam-2.2.1/OpenFOAM-2.2.1/platforms/linux64GccDPOpt/lib/libchemistryModel.so"
#15 Foam-combustionModels- psiChemistryCombustion-psiChemistryCombustion(Foam-word const&, Foam-fvMesh const&) in "/usr/apps1/openfoam-2.2.1/OpenFOAM-2.2.1/platforms/linux64GccDPOpt/lib/libcombustionModels.so"
#16 Foam-combustionModels- aSR<Foam-combustionModels- psiChemistryCombustion>: : PaSR(Foam-word const&, Foam-fvMesh const&) in "/usr/apps1/openfoam-2.2.1/OpenFOAM-2.2.1/platforms/linux64GccDPOpt/lib/libcombustionModels.so"
#17 Foam-combustionModels- psiCombustionModel-adddictionaryConstructorToTable<Foam-combustionModels-PaSR<Foam-combustionModels-psiChemistryCombustion> >-New(Foam-word const&, Foam-fvMesh const&) in "/usr/apps1/openfoam-2.2.1/OpenFOAM-2.2.1/platforms/linux64GccDPOpt/lib/libcombustionModels.so"
#18 Foam-combustionModels- siCombustionModel-New(Foam-fvMesh const&) in "/usr/apps1/openfoam-2.2.1/OpenFOAM-2.2.1/platforms/linux64GccDPOpt/lib/libcombustionModels.so"
#19 main in "/usr/erc/people/krish042/OpenFOAM/krish042-2.2.1/platforms/linux64GccDPOpt/bin/my_reactingFoam"
#20 __libc_start_main in "/lib64/libc.so.6"
#21 __gxx_personality_v0 in "/usr/erc/people/krish042/OpenFOAM/krish042-2.2.1/platforms/linux64GccDPOpt/bin/my_reactingFoam"
Floating exception.

Has anyone else faced such problems. I am still new to C++ and any help would be appreciated. I have attached my source code in case you need to take a look.

Thanks

Hi krishna;

Currently I'm working with reactingFoam without reactions (combustion off) and I'm having this error. Did you solve this?.


thank you :)

If you read the error logs carefully, you can see that the problem is related to the creation of the chemistry model. You are using the Foam chemistry model (as opposed to CHEMKIN); do you have an properly formatted Foam chemistry file?

Hi mTurcios, thank you for your quick reply.

I've been modifying counterFlowFlame2D tutorial (reactingFoam for RANS) cause I need it for LES with a 3D configuration. With LES,the solver runs properly (that's why I think that the problem is not my chemistry file), but when I impose my 3D configuration I'm having the following error:

Attachment 28529

I'm think that my problem is related with my boundary conditions. ¿What do you think?

These are my boundary conditions:

Attachment 28524

I would be so grateful if you could help me. :)

Finally I realized what was happening. One of my species in my 0/ files had a fixed value "0", for instance, if openfoam was calculating the concentration of H2 in a boundary with this condition it would calculate Yi=ni/nT but if n1 is equal to zero the concentration would be infinite.