[Tommy Floessner]

Dear CFD-community,

i tried to simulate a reactive flow of fluids in OF230.
As a base i used the tutorial-case for reactingFoam and applied
several changes that would describe my problem.
The most important change is the main-Solver to rhoReactingFoam.
I'll post the log-file for further Informations about
the used thermophysical-, transport-, combustion-, etc. models.

The Problem in the set up case is a segmentation fault during the initialization
of the field. (see log-file)

I know that the first thing to check, in case of segmentation faults, is
the geometry and boundary conditions.
This i tested sucessfully with a methan-combustion in the reactingFoam tutorial-case.
But i don't want to exclude it as the root of the error by default.


I hope somebody can give me one or two hints, that could help solving the problem.
If any further Information is needed to get rid of the question tags, just ask.


Thanks

Tommy


Log-file:

/*---------------------------------------------------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  2.3.0                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
Build  : 2.3.0-f5222ca19ce6
Exec   : rhoReactingFoam
Date   : Sep 15 2014
Time   : 14:06:53
Host   : "dread-MS-7235"
PID    : 9818
Case   : /home/dread/OpenFOAM/run/rrf_002_poly
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster
allowSystemOperations : Disallowing user-supplied system call operations

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

Create mesh for time = 0


Reading g
Creating reaction model

Selecting combustion model PaSR<rhoChemistryCombustion>
Selecting chemistry type 
{
    chemistrySolver ode;
    chemistryThermo rho;
}

Selecting thermodynamics package 
{
    type            heRhoThermo;
    mixture         reactingMixture;
    transport       polynomial;
    thermo          hPolynomial;
    energy          sensibleEnthalpy;
    equationOfState icoPolynomial;
    specie          specie;
}

AMI: Creating addressing and weights between 500 source faces and 500 target faces
AMI: Patch source sum(weights) min/max/average = 0.999915, 1, 0.999992
AMI: Patch target sum(weights) min/max/average = 1, 1, 1
Selecting chemistryReader foamChemistryReader

// end Log-file

// additional information  

#0  Foam::error::printStack(Foam::Ostream&) in "/home/dread/OpenFOAM/OpenFOAM-2.3.0/platforms/linuxGccDPOpt/lib/libOpenFOAM.so"
#1  Foam::sigFpe::sigHandler(int) in "/home/dread/OpenFOAM/OpenFOAM-2.3.0/platforms/linuxGccDPOpt/lib/libOpenFOAM.so"
#2  Uninterpreted: 
#3   at rhoReactionThermos.C:0
#4  Foam::multiComponentMixture<Foam::polynomialTransport<Foam::species::thermo<Foam::hPolynomialThermo<Foam::icoPolynomial<Foam::specie, 8>, 8>, Foam::sensibleEnthalpy>, 8> >::correctMassFractions() in "/home/dread/OpenFOAM/OpenFOAM-2.3.0/platforms/linuxGccDPOpt/lib/libreactionThermophysicalModels.so"
#5  Foam::multiComponentMixture<Foam::polynomialTransport<Foam::species::thermo<Foam::hPolynomialThermo<Foam::icoPolynomial<Foam::specie, 8>, 8>, Foam::sensibleEnthalpy>, 8> >::multiComponentMixture(Foam::dictionary const&, Foam::List<Foam::word> const&, Foam::HashPtrTable<Foam::polynomialTransport<Foam::species::thermo<Foam::hPolynomialThermo<Foam::icoPolynomial<Foam::specie, 8>, 8>, Foam::sensibleEnthalpy>, 8>, Foam::word, Foam::string::hash> const&, Foam::fvMesh const&) in "/home/dread/OpenFOAM/OpenFOAM-2.3.0/platforms/linuxGccDPOpt/lib/libreactionThermophysicalModels.so"
#6  Foam::reactingMixture<Foam::polynomialTransport<Foam::species::thermo<Foam::hPolynomialThermo<Foam::icoPolynomial<Foam::specie, 8>, 8>, Foam::sensibleEnthalpy>, 8> >::reactingMixture(Foam::dictionary const&, Foam::fvMesh const&) in "/home/dread/OpenFOAM/OpenFOAM-2.3.0/platforms/linuxGccDPOpt/lib/libreactionThermophysicalModels.so"
#7  Foam::heThermo<Foam::rhoReactionThermo, Foam::SpecieMixture<Foam::reactingMixture<Foam::polynomialTransport<Foam::species::thermo<Foam::hPolynomialThermo<Foam::icoPolynomial<Foam::specie, 8>, 8>, Foam::sensibleEnthalpy>, 8> > > >::heThermo(Foam::fvMesh const&, Foam::word const&) in "/home/dread/OpenFOAM/OpenFOAM-2.3.0/platforms/linuxGccDPOpt/lib/libreactionThermophysicalModels.so"
#8  Foam::rhoReactionThermo::addfvMeshConstructorToTable<Foam::heRhoThermo<Foam::rhoReactionThermo, Foam::SpecieMixture<Foam::reactingMixture<Foam::polynomialTransport<Foam::species::thermo<Foam::hPolynomialThermo<Foam::icoPolynomial<Foam::specie, 8>, 8>, Foam::sensibleEnthalpy>, 8> > > > >::New(Foam::fvMesh const&, Foam::word const&) in "/home/dread/OpenFOAM/OpenFOAM-2.3.0/platforms/linuxGccDPOpt/lib/libreactionThermophysicalModels.so"
#9  Foam::autoPtr<Foam::rhoReactionThermo> Foam::basicThermo::New<Foam::rhoReactionThermo>(Foam::fvMesh const&, Foam::word const&) in "/home/dread/OpenFOAM/OpenFOAM-2.3.0/platforms/linuxGccDPOpt/lib/libreactionThermophysicalModels.so"
#10  Foam::rhoReactionThermo::New(Foam::fvMesh const&, Foam::word const&) in "/home/dread/OpenFOAM/OpenFOAM-2.3.0/platforms/linuxGccDPOpt/lib/libreactionThermophysicalModels.so"
#11  Foam::rhoChemistryModel::rhoChemistryModel(Foam::fvMesh const&) in "/home/dread/OpenFOAM/OpenFOAM-2.3.0/platforms/linuxGccDPOpt/lib/libchemistryModel.so"
#12  Foam::chemistryModel<Foam::rhoChemistryModel, Foam::polynomialTransport<Foam::species::thermo<Foam::hPolynomialThermo<Foam::icoPolynomial<Foam::specie, 8>, 8>, Foam::sensibleEnthalpy>, 8> >::chemistryModel(Foam::fvMesh const&) in "/home/dread/OpenFOAM/OpenFOAM-2.3.0/platforms/linuxGccDPOpt/lib/libchemistryModel.so"
#13  Foam::ode<Foam::chemistryModel<Foam::rhoChemistryModel, Foam::polynomialTransport<Foam::species::thermo<Foam::hPolynomialThermo<Foam::icoPolynomial<Foam::specie, 8>, 8>, Foam::sensibleEnthalpy>, 8> > >::ode(Foam::fvMesh const&) in "/home/dread/OpenFOAM/OpenFOAM-2.3.0/platforms/linuxGccDPOpt/lib/libchemistryModel.so"
#14  Foam::rhoChemistryModel::addfvMeshConstructorToTable<Foam::ode<Foam::chemistryModel<Foam::rhoChemistryModel, Foam::polynomialTransport<Foam::species::thermo<Foam::hPolynomialThermo<Foam::icoPolynomial<Foam::specie, 8>, 8>, Foam::sensibleEnthalpy>, 8> > > >::New(Foam::fvMesh const&) in "/home/dread/OpenFOAM/OpenFOAM-2.3.0/platforms/linuxGccDPOpt/lib/libchemistryModel.so"
#15  Foam::autoPtr<Foam::rhoChemistryModel> Foam::basicChemistryModel::New<Foam::rhoChemistryModel>(Foam::fvMesh const&) in "/home/dread/OpenFOAM/OpenFOAM-2.3.0/platforms/linuxGccDPOpt/lib/libchemistryModel.so"
#16  Foam::rhoChemistryModel::New(Foam::fvMesh const&) in "/home/dread/OpenFOAM/OpenFOAM-2.3.0/platforms/linuxGccDPOpt/lib/libchemistryModel.so"
#17  Foam::combustionModels::rhoChemistryCombustion::rhoChemistryCombustion(Foam::word const&, Foam::fvMesh const&) in "/home/dread/OpenFOAM/OpenFOAM-2.3.0/platforms/linuxGccDPOpt/lib/libcombustionModels.so"
#18  Foam::combustionModels::laminar<Foam::combustionModels::rhoChemistryCombustion>::laminar(Foam::word const&, Foam::fvMesh const&) in "/home/dread/OpenFOAM/OpenFOAM-2.3.0/platforms/linuxGccDPOpt/lib/libcombustionModels.so"
#19  Foam::combustionModels::PaSR<Foam::combustionModels::rhoChemistryCombustion>::PaSR(Foam::word const&, Foam::fvMesh const&) in "/home/dread/OpenFOAM/OpenFOAM-2.3.0/platforms/linuxGccDPOpt/lib/libcombustionModels.so"
#20  Foam::combustionModels::rhoCombustionModel::adddictionaryConstructorToTable<Foam::combustionModels::PaSR<Foam::combustionModels::rhoChemistryCombustion> >::New(Foam::word const&, Foam::fvMesh const&) in "/home/dread/OpenFOAM/OpenFOAM-2.3.0/platforms/linuxGccDPOpt/lib/libcombustionModels.so"
#21  Foam::combustionModels::rhoCombustionModel::New(Foam::fvMesh const&) in "/home/dread/OpenFOAM/OpenFOAM-2.3.0/platforms/linuxGccDPOpt/lib/libcombustionModels.so"
#22  
 in "/home/dread/OpenFOAM/OpenFOAM-2.3.0/platforms/linuxGccDPOpt/bin/rhoReactingFoam"
#23  __libc_start_main in "/lib/i386-linux-gnu/libc.so.6"
#24  
 in "/home/dread/OpenFOAM/OpenFOAM-2.3.0/platforms/linuxGccDPOpt/bin/rhoReactingFoam"

[C_Cherry]

Hello!
I use the rhoReactingFoam to solver a chemistry problem, the thermoType is samilar to yours.
I do not know why I can not run my case.I am confused.
Can you help you? Or you can upload your case.
Thanks you!

[wyldckat]

Greetings to all!

@Tommy Floessner: Sorry that no one was able to get to you on this topic. My guess is that there was a division by zero, possibly because the initial temperature for one of the fluids was set to 0 Kelvin. I say this because of these two lines of the output:

#1  Foam::sigFpe::sigHandler(int) in "/home/dread/OpenFOAM/OpenFOAM-2.3.0/platforms/linuxGccDPOpt/lib/libOpenFOAM.so"

#4  Foam::multiComponentMixture<Foam::polynomialTransport<Foam::species::thermo<Foam::hPolynomialThermo<Foam::icoPolynomial<Foam::specie, 8>, 8>, Foam::sensibleEnthalpy>, 8> >::correctMassFractions() in "/home/dread/OpenFOAM/OpenFOAM-2.3.0/platforms/linuxGccDPOpt/lib/libreactionThermophysicalModels.so"

The first line refers to the arithmetic problem, namely "sigFpe" refers to "signal Floating Point Error", i.e. a division by zero or attempting to calculate the logarithm of 0.
The second line indicates that the error occurred when updating/calculating/correcting the mass fractions.


@C_Cherry: Please provide the complete error message that you have gotten, the same way as Tommy did, because your error message might be similar, but could be due to an entirely different reason. Either way, my answer above might help.

Best regards,
Bruno

[C_Cherry]

Thank you for your reply@wyldckat
In my case,the fluid flows in a pipe, the inlet temperature is 473K, the wall temperature is 873K, and the outlet is zeroGradient. The inlet pressure is zeroGradient, and the outlet pressure is 5MPa. The inlet velocity is 0.04244m/s, and the outlet velocity is zeroGradient.
I use the rhoReactingFoam to solve a chemistry problem.
I hope you can give me some help. I appreciate you very much!

Here is the log file.

Create time

Create mesh for time = 0


Reading g
Creating reaction model

Selecting combustion model PaSR<rhoChemistryCombustion>
Selecting chemistry type 
{
    chemistrySolver EulerImplicit;
    chemistryThermo rho;
}

Selecting thermodynamics package 
{
    type            heRhoThermo;
    mixture         reactingMixture;
    transport       polynomial;
    thermo          hPolynomial;
    equationOfState icoPolynomial;
    specie          specie;
    energy          sensibleEnthalpy;
}

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

Reading/calculating face flux field phi

Creating turbulence model.

Selecting turbulence model type RASModel
Selecting RAS turbulence model kEpsilon
kEpsilonCoeffs
{
    Cmu             0.09;
    C1              1.44;
    C2              1.92;
    C3              -0.33;
    sigmak          1;
    sigmaEps        1.3;
    Prt             1;
}

Creating field dpdt

Creating field kinetic energy K

No finite volume options present

Courant Number mean: 4.77042e-06 max: 0.000715562

PIMPLE: Operating solver in PISO mode


Starting time loop

Courant Number mean: 4.77042e-06 max: 0.000715562
deltaT = 1.19904e-05
Time = 1.19904e-05

diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
PIMPLE: iteration 1
DILUPBiCG:  Solving for C10H22, Initial residual = 1, Final residual = 7.07715e-11, No Iterations 2
DILUPBiCG:  Solving for H2, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCG:  Solving for CH4, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCG:  Solving for C2H4, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCG:  Solving for C2H6, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCG:  Solving for C3H6, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCG:  Solving for C3H8, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCG:  Solving for C4H8, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCG:  Solving for C4H10, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCG:  Solving for C5H10, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCG:  Solving for C5H12, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCG:  Solving for C6H12, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCG:  Solving for C6H14, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCG:  Solving for C7H14, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCG:  Solving for C7H16, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCG:  Solving for C8H16, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCG:  Solving for C8H18, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCG:  Solving for C9H18, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCG:  Solving for C9H20, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCG:  Solving for h, Initial residual = 1, Final residual = 3.83483e-08, No Iterations 1
min/max(T) = 473, 873
DICPCG:  Solving for p, Initial residual = 1, Final residual = 56.6857, No Iterations 1001
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 0.00423238, global = -7.44309e-05, cumulative = -7.44309e-05
DICPCG:  Solving for p, Initial residual = 0.179419, Final residual = 0.010874, No Iterations 1001
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 0.000257316, global = -4.964e-06, cumulative = -7.93949e-05
DILUPBiCG:  Solving for epsilon, Initial residual = 0.00233739, Final residual = 8.03855e-08, No Iterations 1
DILUPBiCG:  Solving for k, Initial residual = 1, Final residual = 1.85515e-09, No Iterations 2
ExecutionTime = 4.3 s  ClockTime = 4 s

Courant Number mean: 0.0108272 max: 0.0156534
deltaT = 1.43712e-05
Time = 2.63616e-05

diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
PIMPLE: iteration 1
DILUPBiCG:  Solving for C10H22, Initial residual = 0.227262, Final residual = 1.50724e-09, No Iterations 2
DILUPBiCG:  Solving for H2, Initial residual = 1.01627e-07, Final residual = 1.01627e-07, No Iterations 0
DILUPBiCG:  Solving for CH4, Initial residual = 7.65904e-07, Final residual = 7.65904e-07, No Iterations 0
DILUPBiCG:  Solving for C2H4, Initial residual = 2.39767e-06, Final residual = 1.93074e-10, No Iterations 1
DILUPBiCG:  Solving for C2H6, Initial residual = 1.2649e-06, Final residual = 1.01857e-10, No Iterations 1
DILUPBiCG:  Solving for C3H6, Initial residual = 3.23123e-06, Final residual = 2.60197e-10, No Iterations 1
DILUPBiCG:  Solving for C3H8, Initial residual = 2.64994e-06, Final residual = 2.13387e-10, No Iterations 1
DILUPBiCG:  Solving for C4H8, Initial residual = 3.67143e-06, Final residual = 2.95644e-10, No Iterations 1
DILUPBiCG:  Solving for C4H10, Initial residual = 1.97929e-06, Final residual = 1.59384e-10, No Iterations 1
DILUPBiCG:  Solving for C5H10, Initial residual = 4.00392e-06, Final residual = 3.22418e-10, No Iterations 1
DILUPBiCG:  Solving for C5H12, Initial residual = 2.98688e-06, Final residual = 2.4052e-10, No Iterations 1
DILUPBiCG:  Solving for C6H12, Initial residual = 5.47903e-06, Final residual = 4.41202e-10, No Iterations 1
DILUPBiCG:  Solving for C6H14, Initial residual = 2.56059e-06, Final residual = 2.06193e-10, No Iterations 1
DILUPBiCG:  Solving for C7H14, Initial residual = 5.53991e-06, Final residual = 4.46104e-10, No Iterations 1
DILUPBiCG:  Solving for C7H16, Initial residual = 2.40867e-06, Final residual = 1.9396e-10, No Iterations 1
DILUPBiCG:  Solving for C8H16, Initial residual = 5.69446e-06, Final residual = 4.5855e-10, No Iterations 1
DILUPBiCG:  Solving for C8H18, Initial residual = 4.57642e-07, Final residual = 4.57642e-07, No Iterations 0
DILUPBiCG:  Solving for C9H18, Initial residual = 2.23377e-06, Final residual = 1.79876e-10, No Iterations 1
DILUPBiCG:  Solving for C9H20, Initial residual = 1.28459e-07, Final residual = 1.28459e-07, No Iterations 0
DILUPBiCG:  Solving for h, Initial residual = 0.70559, Final residual = 1.41858e-07, No Iterations 2
min/max(T) = 472.992, 873
DICPCG:  Solving for p, Initial residual = 0.0220879, Final residual = 0.813472, No Iterations 1001
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 0.0237036, global = 9.19515e-05, cumulative = 1.25566e-05
DICPCG:  Solving for p, Initial residual = 0.786122, Final residual = 0.0045141, No Iterations 1001
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 0.000137576, global = 7.97722e-06, cumulative = 2.05338e-05
DILUPBiCG:  Solving for epsilon, Initial residual = 0.00735973, Final residual = 4.25711e-07, No Iterations 1
DILUPBiCG:  Solving for k, Initial residual = 0.723139, Final residual = 9.64343e-07, No Iterations 2
ExecutionTime = 8.52 s  ClockTime = 9 s

Courant Number mean: 0.0152289 max: 0.0186182
deltaT = 1.72256e-05
Time = 4.35873e-05

diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
PIMPLE: iteration 1
DILUPBiCG:  Solving for C10H22, Initial residual = 0.137911, Final residual = 4.1808e-10, No Iterations 2
DILUPBiCG:  Solving for H2, Initial residual = 2.2177e-07, Final residual = 2.2177e-07, No Iterations 0
DILUPBiCG:  Solving for CH4, Initial residual = 1.67136e-06, Final residual = 7.60635e-11, No Iterations 1
DILUPBiCG:  Solving for C2H4, Initial residual = 5.19659e-06, Final residual = 2.36497e-10, No Iterations 1
DILUPBiCG:  Solving for C2H6, Initial residual = 2.7415e-06, Final residual = 1.24766e-10, No Iterations 1
DILUPBiCG:  Solving for C3H6, Initial residual = 7.00318e-06, Final residual = 3.18716e-10, No Iterations 1
DILUPBiCG:  Solving for C3H8, Initial residual = 5.74333e-06, Final residual = 2.6138e-10, No Iterations 1
DILUPBiCG:  Solving for C4H8, Initial residual = 7.95722e-06, Final residual = 3.62134e-10, No Iterations 1
DILUPBiCG:  Solving for C4H10, Initial residual = 4.28982e-06, Final residual = 1.9523e-10, No Iterations 1
DILUPBiCG:  Solving for C5H10, Initial residual = 8.67782e-06, Final residual = 3.94929e-10, No Iterations 1
DILUPBiCG:  Solving for C5H12, Initial residual = 6.4736e-06, Final residual = 2.94614e-10, No Iterations 1
DILUPBiCG:  Solving for C6H12, Initial residual = 1.18748e-05, Final residual = 5.40424e-10, No Iterations 1
DILUPBiCG:  Solving for C6H14, Initial residual = 5.54969e-06, Final residual = 2.52567e-10, No Iterations 1
DILUPBiCG:  Solving for C7H14, Initial residual = 1.20067e-05, Final residual = 5.46429e-10, No Iterations 1
DILUPBiCG:  Solving for C7H16, Initial residual = 5.22043e-06, Final residual = 2.37583e-10, No Iterations 1
DILUPBiCG:  Solving for C8H16, Initial residual = 1.23417e-05, Final residual = 5.61672e-10, No Iterations 1
DILUPBiCG:  Solving for C8H18, Initial residual = 9.98669e-07, Final residual = 9.98669e-07, No Iterations 0
DILUPBiCG:  Solving for C9H18, Initial residual = 4.84136e-06, Final residual = 2.20331e-10, No Iterations 1
DILUPBiCG:  Solving for C9H20, Initial residual = 2.80325e-07, Final residual = 2.80325e-07, No Iterations 0
DILUPBiCG:  Solving for h, Initial residual = 0.348893, Final residual = 2.51228e-07, No Iterations 2
min/max(T) = 472.999, 873
DICPCG:  Solving for p, Initial residual = 0.0961288, Final residual = 0.520595, No Iterations 1001
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 0.00396067, global = 2.93884e-05, cumulative = 4.99222e-05
DICPCG:  Solving for p, Initial residual = 0.356093, Final residual = 0.00296496, No Iterations 1001
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 3.3259e-05, global = 7.44155e-07, cumulative = 5.06663e-05
DILUPBiCG:  Solving for epsilon, Initial residual = 0.00521535, Final residual = 4.15496e-07, No Iterations 1
DILUPBiCG:  Solving for k, Initial residual = 0.470671, Final residual = 1.24424e-07, No Iterations 2
ExecutionTime = 13.51 s  ClockTime = 14 s

Courant Number mean: 0.014602 max: 0.0181641
deltaT = 2.06565e-05
Time = 6.42437e-05
......
......
......
......
......
Courant Number mean: 0.0876726 max: 0.288901
deltaT = 0.002
Time = 0.246

diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
PIMPLE: iteration 1
DILUPBiCG:  Solving for C10H22, Initial residual = 0.00637298, Final residual = 5.46549e-07, No Iterations 4
DILUPBiCG:  Solving for H2, Initial residual = 0.0670608, Final residual = 1.84323e-07, No Iterations 6
DILUPBiCG:  Solving for CH4, Initial residual = 0.108365, Final residual = 2.97852e-07, No Iterations 6
DILUPBiCG:  Solving for C2H4, Initial residual = 0.11579, Final residual = 3.18261e-07, No Iterations 6
DILUPBiCG:  Solving for C2H6, Initial residual = 0.112548, Final residual = 3.09352e-07, No Iterations 6
DILUPBiCG:  Solving for C3H6, Initial residual = 0.116759, Final residual = 3.20924e-07, No Iterations 6
DILUPBiCG:  Solving for C3H8, Initial residual = 0.116146, Final residual = 3.19239e-07, No Iterations 6
DILUPBiCG:  Solving for C4H8, Initial residual = 0.117097, Final residual = 3.21853e-07, No Iterations 6
DILUPBiCG:  Solving for C4H10, Initial residual = 0.115008, Final residual = 3.16112e-07, No Iterations 6
DILUPBiCG:  Solving for C5H10, Initial residual = 0.117303, Final residual = 3.22421e-07, No Iterations 6
DILUPBiCG:  Solving for C5H12, Initial residual = 0.116529, Final residual = 3.20293e-07, No Iterations 6
DILUPBiCG:  Solving for C6H12, Initial residual = 0.117923, Final residual = 3.24124e-07, No Iterations 6
DILUPBiCG:  Solving for C6H14, Initial residual = 0.116027, Final residual = 3.18914e-07, No Iterations 6
DILUPBiCG:  Solving for C7H14, Initial residual = 0.117941, Final residual = 3.24175e-07, No Iterations 6
DILUPBiCG:  Solving for C7H16, Initial residual = 0.115807, Final residual = 3.18308e-07, No Iterations 6
DILUPBiCG:  Solving for C8H16, Initial residual = 0.117987, Final residual = 3.243e-07, No Iterations 6
DILUPBiCG:  Solving for C8H18, Initial residual = 0.101897, Final residual = 2.80074e-07, No Iterations 6
DILUPBiCG:  Solving for C9H18, Initial residual = 0.115517, Final residual = 3.17512e-07, No Iterations 6
DILUPBiCG:  Solving for C9H20, Initial residual = 0.0738386, Final residual = 2.02953e-07, No Iterations 6
DILUPBiCG:  Solving for h, Initial residual = 0.0075258, Final residual = 4.35917e-08, No Iterations 4
min/max(T) = 452.174, 873
DICPCG:  Solving for p, Initial residual = 0.21668, Final residual = 0.195115, No Iterations 1001
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 0.00845023, global = 0.00029595, cumulative = 0.000458058
DICPCG:  Solving for p, Initial residual = 0.125125, Final residual = 0.0586924, No Iterations 1001
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 0.00519283, global = -3.89853e-05, cumulative = 0.000419072
DILUPBiCG:  Solving for epsilon, Initial residual = 0.0448897, Final residual = 4.74853e-08, No Iterations 4
DILUPBiCG:  Solving for k, Initial residual = 0.321046, Final residual = 8.4238e-07, No Iterations 4
ExecutionTime = 958.48 s  ClockTime = 958 s

Courant Number mean: 0.152591 max: 0.715227
deltaT = 0.001
Time = 0.247

diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
PIMPLE: iteration 1
DILUPBiCG:  Solving for C10H22, Initial residual = 0.00430238, Final residual = 4.93826e-07, No Iterations 3
DILUPBiCG:  Solving for H2, Initial residual = 0.085043, Final residual = 2.71822e-07, No Iterations 5
DILUPBiCG:  Solving for CH4, Initial residual = 0.107737, Final residual = 3.44359e-07, No Iterations 5
DILUPBiCG:  Solving for C2H4, Initial residual = 0.110815, Final residual = 3.54199e-07, No Iterations 5
DILUPBiCG:  Solving for C2H6, Initial residual = 0.1095, Final residual = 3.49996e-07, No Iterations 5
DILUPBiCG:  Solving for C3H6, Initial residual = 0.1112, Final residual = 3.55429e-07, No Iterations 5
DILUPBiCG:  Solving for C3H8, Initial residual = 0.110957, Final residual = 3.54652e-07, No Iterations 5
DILUPBiCG:  Solving for C4H8, Initial residual = 0.111334, Final residual = 3.55855e-07, No Iterations 5
DILUPBiCG:  Solving for C4H10, Initial residual = 0.110502, Final residual = 3.53198e-07, No Iterations 5
DILUPBiCG:  Solving for C5H10, Initial residual = 0.111415, Final residual = 3.56116e-07, No Iterations 5
DILUPBiCG:  Solving for C5H12, Initial residual = 0.111109, Final residual = 3.55139e-07, No Iterations 5
DILUPBiCG:  Solving for C6H12, Initial residual = 0.111658, Final residual = 3.56892e-07, No Iterations 5
DILUPBiCG:  Solving for C6H14, Initial residual = 0.11091, Final residual = 3.54502e-07, No Iterations 5
DILUPBiCG:  Solving for C7H14, Initial residual = 0.111665, Final residual = 3.56915e-07, No Iterations 5
DILUPBiCG:  Solving for C7H16, Initial residual = 0.110822, Final residual = 3.54221e-07, No Iterations 5
DILUPBiCG:  Solving for C8H16, Initial residual = 0.111683, Final residual = 3.56971e-07, No Iterations 5
DILUPBiCG:  Solving for C8H18, Initial residual = 0.104853, Final residual = 3.35141e-07, No Iterations 5
DILUPBiCG:  Solving for C9H18, Initial residual = 0.110706, Final residual = 3.53851e-07, No Iterations 5
DILUPBiCG:  Solving for C9H20, Initial residual = 0.0895874, Final residual = 2.86347e-07, No Iterations 5
DILUPBiCG:  Solving for h, Initial residual = 0.00516729, Final residual = 9.07852e-07, No Iterations 3
min/max(T) = 452.165, 873
DICPCG:  Solving for p, Initial residual = 0.140282, Final residual = 2.44093, No Iterations 1001
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 0.083028, global = 0.000464274, cumulative = 0.000883346
DICPCG:  Solving for p, Initial residual = 0.468191, Final residual = 0.0195225, No Iterations 1001
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 0.00396294, global = -5.42544e-05, cumulative = 0.000829091
DILUPBiCG:  Solving for epsilon, Initial residual = 0.0475024, Final residual = 3.84935e-07, No Iterations 3
DILUPBiCG:  Solving for k, Initial residual = 0.389679, Final residual = 5.24729e-07, No Iterations 4
ExecutionTime = 962.36 s  ClockTime = 962 s

Courant Number mean: 0.172735 max: 0.896032
deltaT = 0.000428571
Time = 0.247429

diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
PIMPLE: iteration 1
#0  Foam::error::printStack(Foam::Ostream&) at ??:?
#1  Foam::sigFpe::sigHandler(int) at ??:?
#2  ? in "/home/wangchang/OpenFOAM/OpenFOAM-2.4.0/platforms/linuxGccDPOpt/lib/libOpenFOAM.so"
#3  Foam::EulerImplicit<Foam::chemistryModel<Foam::rhoChemistryModel, Foam::polynomialTransport<Foam::species::thermo<Foam::hPolynomialThermo<Foam::icoPolynomial<Foam::specie, 8>, 8>, Foam::sensibleEnthalpy>, 8> > >::solve(Foam::Field<double>&, double&, double&, double&, double&) const at ??:?
#4  double Foam::chemistryModel<Foam::rhoChemistryModel, Foam::polynomialTransport<Foam::species::thermo<Foam::hPolynomialThermo<Foam::icoPolynomial<Foam::specie, 8>, 8>, Foam::sensibleEnthalpy>, 8> >::solve<Foam::UniformField<double> >(Foam::UniformField<double> const&) at ??:?
#5  Foam::chemistryModel<Foam::rhoChemistryModel, Foam::polynomialTransport<Foam::species::thermo<Foam::hPolynomialThermo<Foam::icoPolynomial<Foam::specie, 8>, 8>, Foam::sensibleEnthalpy>, 8> >::solve(double) at ??:?
#6  Foam::combustionModels::laminar<Foam::combustionModels::rhoChemistryCombustion>::correct() at ??:?
#7  Foam::combustionModels::PaSR<Foam::combustionModels::rhoChemistryCombustion>::correct() at ??:?
#8  ? at ??:?
#9  __libc_start_main in "/lib/i386-linux-gnu/libc.so.6"
#10  ? at ??:?

[wyldckat]

Greetings C_Cherry,

My apologies, it's being nearly a year, and I never managed to find the time to look into the problem you were having. Hopefully you managed to solve the issue.

If you haven't, the crash seems to have occurred when it was trying to calculate the equation for C10H22, due to not being able to solve one of the polynomial expressions (resulted in an invalid floating point operation, hence the sigFpe)... not sure which polynomial it was exactly, if it was receiving pressure and returning temperature or vice-versa.

My guess is that the initial pressure inside the domain was set to 0, instead of 1 or 5 MPa, which would somewhat explain the sudden crash due to not being able to calculate things properly, e.g. using a very high pressure as one of the input variables into the polynomial, could not render a valid solution.

Best regards,
Bruno