CFD Online Discussion Forums

CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)
-   OpenFOAM Running, Solving & CFD (https://www.cfd-online.com/Forums/openfoam-solving/)
-   -   coalChemistryFoam: "Foam::error::printStack(Foam::Ostream&) at ??:?" (https://www.cfd-online.com/Forums/openfoam-solving/148303-coalchemistryfoam-foam-error-printstack-foam-ostream.html)

musabai February 9, 2015 02:24

coalChemistryFoam: "Foam::error::printStack(Foam::Ostream&) at ??:?"
 
Hi everyone, I am new to OF.

I want to run a case using my own geometry. I copy everything from coalChemistryFoam, but change geometry to my own. I get this error while running this. Don't know how to solve it. Perhaps someone help me? Thank you very much. Mesh is ok.

Solving 2-D cloud coalCloud1
#0 Foam::error::printStack(Foam::Ostream&) at ??:?
#1 Foam::sigFpe::sigHandler(int) at ??:?
#2 in "/lib/x86_64-linux-gnu/libc.so.6"

PHP Code:

/*---------------------------------------------------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  2.3.1                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
Build  2.3.1-bcfaaa7b8660
Exec   
coalChemistryFoam
Date   
Feb 08 2015
Time   
23:06:43
Host   
"baosheng"
PID    9163
Case   : /home/baosheng/OpenFOAM/baosheng-2.3.1/run/researchcoalCombustion/coalChemistryFoam/simplifiedSiwek
nProcs 
1
sigFpe 
Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking Monitoring run-time modified files using timeStampMaster
allowSystemOperations 
Allowing user-supplied system call operations

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

Create mesh 
for time 0


Reading g
Creating combustion 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
chemistryModel
Number of species and reactions 2
Selecting ODE solver seulex
    using integrated reaction rate
Creating component thermo properties
:
    
multi-component carrier 6 species
    liquids 
1 components
    solids 
2 components

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

Creating finite volume options from 
"system/fvOptions"

Selecting finite volume options model type fixedTemperatureConstraint
    Source
source1
    
applying source at time 0.1 for duration 0.4
    
selecting cells using cellSet ignitionCells
    
selected 0 cell(swith volume 0


Constructing coal cloud
Constructing particle forces
    Selecting particle force sphereDrag
    Selecting particle force gravity
Constructing cloud functions
    none
Constructing particle injection models
Creating injector
model1
Selecting injection model manualInjection
    Constructing 2
-D injection
Selecting distribution model RosinRammler
Selecting dispersion model stochasticDispersionRAS
Selecting patch interaction model standardWallInteraction
Selecting stochastic collision model none
Selecting surface film model none
Selecting U integration scheme Euler
Selecting heat transfer model RanzMarshall
Selecting T integration scheme analytical
Selecting composition model singleMixtureFraction
--> FOAM Warning 
    
From function phaseProperties::initialiseGlobalIds(...)
    
in file phaseProperties/phaseProperties/phaseProperties.C at line 231
    Assuming no mapping between solid 
and carrier species
Selecting phase change model liquidEvaporation
Participating liquid species
:
    
H2O
Selecting devolatilisation model constantRateDevolatilisation
Participating volatile species
:
    
CH4particle mass fraction 0.127444
    H2
particle mass fraction 0.020889
    CO2
particle mass fraction 0.062667
Selecting surface reaction model COxidationKineticDiffusionLimitedRate
    C
(s): particle mass fraction 0.659000048

Constructing limestone cloud
Constructing particle forces
    Selecting particle force sphereDrag
    Selecting particle force gravity
Constructing cloud functions
    none
Constructing particle injection models
Creating injector
model1
Selecting injection model manualInjection
    Constructing 2
-D injection
Selecting distribution model RosinRammler
Selecting dispersion model stochasticDispersionRAS
Selecting patch interaction model standardWallInteraction
Selecting stochastic collision model none
Selecting surface film model none
Selecting U integration scheme Euler
Selecting heat transfer model RanzMarshall
Selecting T integration scheme analytical
Selecting radiationModel P1
Selecting absorptionEmissionModel binaryAbsorptionEmission
Selecting absorptionEmissionModel constantAbsorptionEmission
Selecting absorptionEmissionModel cloudAbsorptionEmission
Selecting scatterModel cloudScatter
Courant Number mean
0 max0

PIMPLE
Operating solver in PISO mode


Starting time loop

Courant Number mean
0 max0
deltaT 
0.000119047619
Time 
0.000119048


Solving 2
-D cloud coalCloud1
#0  Foam::error::printStack(Foam::Ostream&) at ??:?
#1  Foam::sigFpe::sigHandler(int) at ??:?
#2   in "/lib/x86_64-linux-gnu/libc.so.6"
#3  void Foam::ReactingMultiphaseParcel<Foam::ReactingParcel<Foam::ThermoParcel<Foam::KinematicParcel<Foam::particle> > > >::calc<Foam::ReactingParcel<Foam::ThermoParcel<Foam::KinematicParcel<Foam::particle> > >::TrackingData<Foam::ReactingMultiphaseCloud<Foam::ReactingCloud<Foam::ThermoCloud<Foam::KinematicCloud<Foam::Cloud<Foam::ReactingMultiphaseParcel<Foam::ReactingParcel<Foam::ThermoParcel<Foam::KinematicParcel<Foam::particle> > > > > > > > > > >(Foam::ReactingParcel<Foam::ThermoParcel<Foam::KinematicParcel<Foam::particle> > >::TrackingData<Foam::ReactingMultiphaseCloud<Foam::ReactingCloud<Foam::ThermoCloud<Foam::KinematicCloud<Foam::Cloud<Foam::ReactingMultiphaseParcel<Foam::ReactingParcel<Foam::ThermoParcel<Foam::KinematicParcel<Foam::particle> > > > > > > > > >&, double, int) at ??:?
#4  bool Foam::KinematicParcel<Foam::particle>::move<Foam::ReactingParcel<Foam::ThermoParcel<Foam::KinematicParcel<Foam::particle> > >::TrackingData<Foam::ReactingMultiphaseCloud<Foam::ReactingCloud<Foam::ThermoCloud<Foam::KinematicCloud<Foam::Cloud<Foam::ReactingMultiphaseParcel<Foam::ReactingParcel<Foam::ThermoParcel<Foam::KinematicParcel<Foam::particle> > > > > > > > > > >(Foam::ReactingParcel<Foam::ThermoParcel<Foam::KinematicParcel<Foam::particle> > >::TrackingData<Foam::ReactingMultiphaseCloud<Foam::ReactingCloud<Foam::ThermoCloud<Foam::KinematicCloud<Foam::Cloud<Foam::ReactingMultiphaseParcel<Foam::ReactingParcel<Foam::ThermoParcel<Foam::KinematicParcel<Foam::particle> > > > > > > > > >&, double) at ??:?
#5  void Foam::InjectionModel<Foam::KinematicCloud<Foam::Cloud<Foam::ReactingMultiphaseParcel<Foam::ReactingParcel<Foam::ThermoParcel<Foam::KinematicParcel<Foam::particle> > > > > > >::inject<Foam::ReactingParcel<Foam::ThermoParcel<Foam::KinematicParcel<Foam::particle> > >::TrackingData<Foam::ReactingMultiphaseCloud<Foam::ReactingCloud<Foam::ThermoCloud<Foam::KinematicCloud<Foam::Cloud<Foam::ReactingMultiphaseParcel<Foam::ReactingParcel<Foam::ThermoParcel<Foam::KinematicParcel<Foam::particle> > > > > > > > > > >(Foam::ReactingParcel<Foam::ThermoParcel<Foam::KinematicParcel<Foam::particle> > >::TrackingData<Foam::ReactingMultiphaseCloud<Foam::ReactingCloud<Foam::ThermoCloud<Foam::KinematicCloud<Foam::Cloud<Foam::ReactingMultiphaseParcel<Foam::ReactingParcel<Foam::ThermoParcel<Foam::KinematicParcel<Foam::particle> > > > > > > > > >&) at ??:?
#6  void Foam::KinematicCloud<Foam::Cloud<Foam::ReactingMultiphaseParcel<Foam::ReactingParcel<Foam::ThermoParcel<Foam::KinematicParcel<Foam::particle> > > > > >::evolveCloud<Foam::ReactingParcel<Foam::ThermoParcel<Foam::KinematicParcel<Foam::particle> > >::TrackingData<Foam::ReactingMultiphaseCloud<Foam::ReactingCloud<Foam::ThermoCloud<Foam::KinematicCloud<Foam::Cloud<Foam::ReactingMultiphaseParcel<Foam::ReactingParcel<Foam::ThermoParcel<Foam::KinematicParcel<Foam::particle> > > > > > > > > > >(Foam::ReactingParcel<Foam::ThermoParcel<Foam::KinematicParcel<Foam::particle> > >::TrackingData<Foam::ReactingMultiphaseCloud<Foam::ReactingCloud<Foam::ThermoCloud<Foam::KinematicCloud<Foam::Cloud<Foam::ReactingMultiphaseParcel<Foam::ReactingParcel<Foam::ThermoParcel<Foam::KinematicParcel<Foam::particle> > > > > > > > > >&) at ??:?
#7  void Foam::KinematicCloud<Foam::Cloud<Foam::ReactingMultiphaseParcel<Foam::ReactingParcel<Foam::ThermoParcel<Foam::KinematicParcel<Foam::particle> > > > > >::solve<Foam::ReactingParcel<Foam::ThermoParcel<Foam::KinematicParcel<Foam::particle> > >::TrackingData<Foam::ReactingMultiphaseCloud<Foam::ReactingCloud<Foam::ThermoCloud<Foam::KinematicCloud<Foam::Cloud<Foam::ReactingMultiphaseParcel<Foam::ReactingParcel<Foam::ThermoParcel<Foam::KinematicParcel<Foam::particle> > > > > > > > > > >(Foam::ReactingParcel<Foam::ThermoParcel<Foam::KinematicParcel<Foam::particle> > >::TrackingData<Foam::ReactingMultiphaseCloud<Foam::ReactingCloud<Foam::ThermoCloud<Foam::KinematicCloud<Foam::Cloud<Foam::ReactingMultiphaseParcel<Foam::ReactingParcel<Foam::ThermoParcel<Foam::KinematicParcel<Foam::particle> > > > > > > > > >&) at ??:?
#8  Foam::ReactingMultiphaseCloud<Foam::ReactingCloud<Foam::ThermoCloud<Foam::KinematicCloud<Foam::Cloud<Foam::ReactingMultiphaseParcel<Foam::ReactingParcel<Foam::ThermoParcel<Foam::KinematicParcel<Foam::particle> > > > > > > > >::evolve() at ??:?
#9  
 
at ??:?
#10  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#11  
 
at ??:?
Floating point exception (core dumped


ssss February 9, 2015 04:42

That is an error because there is a division by zero. Check your initial conditions and boundaries.

Be also careful because you are operating a PIMPLE in Piso mode, maybe you want to change it:

PIMPLE: Operating solver in PISO mode

Check in openfoamwiki.net the Pimple guide

musabai February 20, 2015 14:07

Hi SSSS,

I changed my geometry. At first in the z-direction the thickness is just 0.001m. It had this problem. Then I changed to 0.01 m and it had this problem too. but when the thickness is 0.1 m, there is no problem.


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