[zhixuan]

hi all,

As a newbie of OF, I'm trying to simulate flame in a combustion chamber. As a start, I'm thinking to modify the LTSReactingFoam tutorial, which is a case for 2D laminar counter flow. I only imported my own geometry and changed the boundary condition accordingly. The pre-process runs very well and I can read the geometry and initial condition from paraView. However, it fails to start the first iteration and here is the output when I run LTSReactingFoam:

Create time

Create mesh for time = 0


PIMPLE: Operating solver in PISO mode


Reading g
Creating reaction model

Selecting combustion model laminar<psiChemistryCombustion>
Selecting chemistry type 
{
    chemistrySolver EulerImplicit;
    chemistryThermo psi;
}

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

Selecting chemistryReader foamChemistryReader
chemistryModel: Number of species = 5 and reactions = 1
    using integrated reaction rate
Reading field U

Reading/calculating face flux field phi

Creating turbulence model.

Selecting turbulence model type laminar
Creating field dpdt

Creating field kinetic energy K

No finite volume options present

#0  Foam::error::printStack(Foam::Ostream&) in "/opt/openfoam230/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#1  Foam::sigFpe::sigHandler(int) in "/opt/openfoam230/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#2   in "/lib/x86_64-linux-gnu/libc.so.6"
#3  Foam::divide(Foam::Field<double>&, Foam::UList<double> const&, Foam::UList<double> const&) in "/opt/openfoam230/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#4  Foam::operator/(Foam::UList<double> const&, Foam::UList<double> const&) in "/opt/openfoam230/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#5  
 at /home/opencfd/OpenFOAM/OpenFOAM-2.3.0/src/finiteVolume/lnInclude/compressibleCourantNo.H:41
#6  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#7  
 in "/opt/openfoam230/platforms/linux64GccDPOpt/bin/LTSReactingFoam"
Floating point exception (core dumped)

I also repeated the same processes in ReactingFoam tutorial. It reports the same thing. So I'm wondering if there is any hints on how to debug this problem.

Thank you a lot!

[zhixuan]

if the code is needed, here is the link for download.
https://www.dropbox.com/s/k6a89q9oz0...oam_lam.tar.gz

Thank you