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/)
-   -   problem with solving lagrange reaction cloud (https://www.cfd-online.com/Forums/openfoam-solving/134489-problem-solving-lagrange-reaction-cloud.html)

Polli April 30, 2014 07:53
problem with solving lagrange reaction cloud
 
Hello to all!
I have the folowing problem:
I simulate on OpenFOAM 2.2.2 an lagrance cloud which consists of ash, C, and C16H34 as a liquid. The evaporation i simulate with the liquid evapouration model. The simulation deals with particles which fly trough a hot gas and start to burn. The simulation runs perfect for the fist 10 000 interations but then i always get the following error message. I tried to reduce the timestep but this also not work. And one more strange thing, after the breakdown of the simulation, wenn i restart it from the last timestep i also get an error message. Below i poste the console output with the two last normal timesteps and the error message. After the error message i restartet the simulation and then there was the next error.

Thank you for your help, yours Polli

CONSOLE OUTPUT:

Courant Number mean: 0.0008352037363 max: 0.02073954834
deltaT = 4.994994995e-08
Time = 0.00184316


Solving 3-D cloud coalCloud1
Cloud: coalCloud1
Current number of parcels = 200
Current mass in system = 3.536393974e-07
Linear momentum = (-2.329083518e-08 2.185374263e-08 -6.694384904e-08)
|Linear momentum| = 7.41722859e-08
Linear kinetic energy = 5.237042721e-07
Rotational kinetic energy = 0
model1:
number of parcels added = 200
mass introduced = 3.690682587e-07
Parcels absorbed into film = 0
New film detached parcels = 0
Parcel fate (number, mass)
- escape = 0, 0
- stick = 0, 0
Temperature min/max = 309.6698539, 1321.440261
Mass transfer phase change = 1.542162181e-08
Mass transfer devolatilisation = 0
Mass transfer surface reaction = 7.23953612e-12

diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCG: Solving for Ux, Initial residual = 0.001104156049, Final residual = 9.873439027e-26, No Iterations 39
DILUPBiCG: Solving for Uy, Initial residual = 0.0005135342004, Final residual = 1.822588789e-26, No Iterations 40
DILUPBiCG: Solving for Uz, Initial residual = 0.0005317415987, Final residual = 9.042617802e-26, No Iterations 40
DILUPBiCG: Solving for CH4, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCG: Solving for C16H34, Initial residual = 1.11068148e-05, Final residual = 4.869725666e-09, No Iterations 1
DILUPBiCG: Solving for O2, Initial residual = 1.886436085e-05, Final residual = 5.272595733e-09, No Iterations 1
DILUPBiCG: Solving for H2, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCG: Solving for CO2, Initial residual = 2.05067378e-05, Final residual = 5.533931152e-09, No Iterations 1
DILUPBiCG: Solving for H2O, Initial residual = 2.049604594e-05, Final residual = 5.535477216e-09, No Iterations 1
DILUPBiCG: Solving for h, Initial residual = 1.869263163e-05, Final residual = 2.909299214e-26, No Iterations 39
T gas min/max = 312.9801108, 2822.970886
T gas min/max = 534.7538798, 948.6851826
T gas min/max = 312.9821183, 2822.878886
GAMG: Solving for p, Initial residual = 0.008959627278, Final residual = 9.390391906e-08, No Iterations 1
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 8.316632968e-10, global = -4.769657762e-10, cumulative = -9.890697957e-08
GAMG: Solving for p, Initial residual = 1.105018814e-07, Final residual = 4.992748749e-09, No Iterations 1
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 2.016364421e-11, global = -9.559535127e-13, cumulative = -9.890793552e-08
GAMG: Solving for p, Initial residual = 5.365688439e-09, Final residual = 5.365688439e-09, No Iterations 0
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 2.139834542e-11, global = -6.946863034e-13, cumulative = -9.890863021e-08
DILUPBiCG: Solving for epsilon, Initial residual = 5.823533933e-05, Final residual = 1.163355584e-26, No Iterations 38
DILUPBiCG: Solving for k, Initial residual = 0.0001720544673, Final residual = 9.004543255e-26, No Iterations 38
ExecutionTime = 22044.09 s ClockTime = 22890 s

Courant Number mean: 0.0008349075938 max: 0.02088308343
deltaT = 4.994994995e-08
Time = 0.00184321


Solving 3-D cloud coalCloud1
[5] #0 Foam::error::printStack(Foam::Ostream&) at ??:?
[5] #1 Foam::sigFpe::sigHandler(int) at ??:?
[5] #2 in "/lib64/libc.so.6"
[5] #3 Foam::liquidMixtureProperties::X(Foam::Field<doubl e> const&) const at ??:?
[5] #4 Foam::LiquidEvaporation<Foam::ReactingCloud<Foam:: ThermoCloud<Foam::KinematicCloud<Foam::Cloud<Foam: :ReactingMultiphaseParcel<Foam::ReactingParcel<Foa m::ThermoParcel<Foam::KinematicParcel<Foam:article > > > > > > > > >::Tvap(Foam::Field<double> const&) const at ??:?
[5] #5 void Foam::ReactingParcel<Foam::ThermoParcel<Foam::Kine maticParcel<Foam:article> > >::calcPhaseChange<Foam::ReactingParcel<Foam::Ther moParcel<Foam::KinematicParcel<Foam:article> > >::TrackingData<Foam::ReactingMultiphaseCloud<Foam ::ReactingCloud<Foam::ThermoCloud<Foam::KinematicC loud<Foam::Cloud<Foam::ReactingMultiphaseParcel<Fo am::ReactingParcel<Foam::ThermoParcel<Foam::Kinema ticParcel<Foam:article> > > > > > > > > > >(Foam::ReactingParcel<Foam::ThermoParcel<Foam::Ki nematicParcel<Foam:article> > >::TrackingData<Foam::ReactingMultiphaseCloud<Foam ::ReactingCloud<Foam::ThermoCloud<Foam::KinematicC loud<Foam::Cloud<Foam::ReactingMultiphaseParcel<Fo am::ReactingParcel<Foam::ThermoParcel<Foam::Kinema ticParcel<Foam:article> > > > > > > > > >&, double, int, double, double, double, double, double, double, double, int, double, Foam::Field<double> const&, Foam::Field<double>&, double&, double&, double&, Foam::Field<double>&) at ??:?
[5] #6 void Foam::ReactingMultiphaseParcel<Foam::ReactingParce l<Foam::ThermoParcel<Foam::KinematicParcel<Foam: article> > > >::calc<Foam::ReactingParcel<Foam::ThermoParcel<Fo am::KinematicParcel<Foam: article> > >::TrackingData<Foam::ReactingMultiphaseCloud<Foam ::ReactingCloud<Foam::ThermoCloud<Foam::KinematicC loud<Foam::Cloud<Foam::ReactingMultiphaseParcel<Fo am::ReactingParcel<Foam::ThermoParcel<Foam::Kinema ticParcel<Foam: article> > > > > > > > > > >(Foam::ReactingParcel<Foam::ThermoParcel<Foam::Ki nematicParcel<Foam: article> > >::TrackingData<Foam::ReactingMultiphaseCloud<Foam ::ReactingCloud<Foam::ThermoCloud<Foam::KinematicC loud<Foam::Cloud<Foam::ReactingMultiphaseParcel<Fo am::ReactingParcel<Foam::ThermoParcel<Foam::Kinema ticParcel<Foam: article> > > > > > > > > >&, double, int) at ??:?
[5] #7 bool Foam::KinematicParcel<Foam: article>::move<Foam::ReactingParcel<Foam::ThermoPa rcel<Foam::KinematicParcel<Foam: article> > >::TrackingData<Foam::ReactingMultiphaseCloud<Foam ::ReactingCloud<Foam::ThermoCloud<Foam::KinematicC loud<Foam::Cloud<Foam::ReactingMultiphaseParcel<Fo am::ReactingParcel<Foam::ThermoParcel<Foam::Kinema ticParcel<Foam: article> > > > > > > > > > >(Foam::ReactingParcel<Foam::ThermoParcel<Foam::Ki nematicParcel<Foam: article> > >::TrackingData<Foam::ReactingMultiphaseCloud<Foam ::ReactingCloud<Foam::ThermoCloud<Foam::KinematicC loud<Foam::Cloud<Foam::ReactingMultiphaseParcel<Fo am::ReactingParcel<Foam::ThermoParcel<Foam::Kinema ticParcel<Foam: article> > > > > > > > > >&, double) at ??:?
[5] #8 void Foam::Cloud<Foam::ReactingMultiphaseParcel<Foam::R eactingParcel<Foam::ThermoParcel<Foam::KinematicPa rcel<Foam: article> > > > >::move<Foam::ReactingParcel<Foam::ThermoParcel<Fo am::KinematicParcel<Foam: article> > >::TrackingData<Foam::ReactingMultiphaseCloud<Foam ::ReactingCloud<Foam::ThermoCloud<Foam::KinematicC loud<Foam::Cloud<Foam::ReactingMultiphaseParcel<Fo am::ReactingParcel<Foam::ThermoParcel<Foam::Kinema ticParcel<Foam: article> > > > > > > > > > >(Foam::ReactingParcel<Foam::ThermoParcel<Foam::Ki nematicParcel<Foam: article> > >::TrackingData<Foam::ReactingMultiphaseCloud<Foam ::ReactingCloud<Foam::ThermoCloud<Foam::KinematicC loud<Foam::Cloud<Foam::ReactingMultiphaseParcel<Fo am::ReactingParcel<Foam::ThermoParcel<Foam::Kinema ticParcel<Foam: article> > > > > > > > > >&, double) at ??:?
[5] #9 void Foam::KinematicCloud<Foam::Cloud<Foam::ReactingMul tiphaseParcel<Foam::ReactingParcel<Foam::ThermoPar cel<Foam::KinematicParcel<Foam: article> > > > > >::evolveCloud<Foam::ReactingParcel<Foam::ThermoPa rcel<Foam::KinematicParcel<Foam: article> > >::TrackingData<Foam::ReactingMultiphaseCloud<Foam ::ReactingCloud<Foam::ThermoCloud<Foam::KinematicC loud<Foam::Cloud<Foam::ReactingMultiphaseParcel<Fo am::ReactingParcel<Foam::ThermoParcel<Foam::Kinema ticParcel<Foam: article> > > > > > > > > > >(Foam::ReactingParcel<Foam::ThermoParcel<Foam::Ki nematicParcel<Foam: article> > >::TrackingData<Foam::ReactingMultiphaseCloud<Foam ::ReactingCloud<Foam::ThermoCloud<Foam::KinematicC loud<Foam::Cloud<Foam::ReactingMultiphaseParcel<Fo am::ReactingParcel<Foam::ThermoParcel<Foam::Kinema ticParcel<Foam: article> > > > > > > > > >&) at ??:?
[5] #10 void Foam::KinematicCloud<Foam::Cloud<Foam::ReactingMul tiphaseParcel<Foam::ReactingParcel<Foam::ThermoPar cel<Foam::KinematicParcel<Foam: article> > > > > >::solve<Foam::ReactingParcel<Foam::ThermoParcel<F oam::KinematicParcel<Foam: article> > >::TrackingData<Foam::ReactingMultiphaseCloud<Foam ::ReactingCloud<Foam::ThermoCloud<Foam::KinematicC loud<Foam::Cloud<Foam::ReactingMultiphaseParcel<Fo am::ReactingParcel<Foam::ThermoParcel<Foam::Kinema ticParcel<Foam: article> > > > > > > > > > >(Foam::ReactingParcel<Foam::ThermoParcel<Foam::Ki nematicParcel<Foam: article> > >::TrackingData<Foam::ReactingMultiphaseCloud<Foam ::ReactingCloud<Foam::ThermoCloud<Foam::KinematicC loud<Foam::Cloud<Foam::ReactingMultiphaseParcel<Fo am::ReactingParcel<Foam::ThermoParcel<Foam::Kinema ticParcel<Foam: article> > > > > > > > > >&) at ??:?
[5] #11
[5] at ??:?
[5] #12 __libc_start_main in "/lib64/libc.so.6"
[5] #13
[5] at /home/abuild/rpmbuild/BUILD/glibc-2.18/csu/../sysdeps/x86_64/start.S:125
[linux-finq:07753] *** Process received signal ***
[linux-finq:07753] Signal: Floating point exception (8)
[linux-finq:07753] Signal code: (-6)
[linux-finq:07753] Failing at address: 0x3e800001e49
[linux-finq:07753] [ 0] /lib64/libc.so.6(+0x358d0) [0x7f75efdfe8d0]
[linux-finq:07753] [ 1] /lib64/libc.so.6(gsignal+0x39) [0x7f75efdfe849]
[linux-finq:07753] [ 2] /lib64/libc.so.6(+0x358d0) [0x7f75efdfe8d0]
[linux-finq:07753] [ 3] /opt/OpenFOAM-2.2.2/platforms/linux64GccDPOpt/lib/libliquidMixtureProperties.so(_ZNK4Foam23liquidMix tureProperties1XERKNS_5FieldIdEE+0x243) [0x7f75f3fd9c63]
[linux-finq:07753] [ 4] /opt/OpenFOAM-2.2.2/platforms/linux64GccDPOpt/lib/liblagrangianIntermediate.so(_ZNK4Foam17LiquidEvap orationINS_13ReactingCloudINS_11ThermoCloudINS_14K inematicCloudINS_5CloudINS_24ReactingMultiphasePar celINS_14ReactingParcelINS_12ThermoParcelINS_15Kin ematicParcelINS_8particleEEEEEEEEEEEEEEEEEE4TvapER KNS_5FieldIdEE+0x1c) [0x7f75f55036bc]
[linux-finq:07753] [ 5] MycoalChemistryFoam4(_ZN4Foam14ReactingParcelINS_1 2ThermoParcelINS_15KinematicParcelINS_8particleEEE EEE15calcPhaseChangeINS6_12TrackingDataINS_23React ingMultiphaseCloudINS_13ReactingCloudINS_11ThermoC loudINS_14KinematicCloudINS_5CloudINS_24ReactingMu ltiphaseParcelIS6_EEEEEEEEEEEEEEEEvRT_didddddddidR KNS_5FieldIdEERSP_RdST_ST_SS_+0x82) [0x49f472]
[linux-finq:07753] [ 6] MycoalChemistryFoam4(_ZN4Foam24ReactingMultiphaseP arcelINS_14ReactingParcelINS_12ThermoParcelINS_15K inematicParcelINS_8particleEEEEEEEE4calcINS7_12Tra ckingDataINS_23ReactingMultiphaseCloudINS_13Reacti ngCloudINS_11ThermoCloudINS_14KinematicCloudINS_5C loudIS8_EEEEEEEEEEEEEEvRT_di+0x399) [0x4a2179]
[linux-finq:07753] [ 7] MycoalChemistryFoam4(_ZN4Foam15KinematicParcelINS_ 8particleEE4moveINS_14ReactingParcelINS_12ThermoPa rcelIS2_EEE12TrackingDataINS_23ReactingMultiphaseC loudINS_13ReactingCloudINS_11ThermoCloudINS_14Kine maticCloudINS_5CloudINS_24ReactingMultiphaseParcel IS7_EEEEEEEEEEEEEEEEbRT_d+0x3d0) [0x4d9b50]
[linux-finq:07753] [ 8] MycoalChemistryFoam4(_ZN4Foam5CloudINS_24ReactingM ultiphaseParcelINS_14ReactingParcelINS_12ThermoPar celINS_15KinematicParcelINS_8particleEEEEEEEEEE4mo veINS8_12TrackingDataINS_23ReactingMultiphaseCloud INS_13ReactingCloudINS_11ThermoCloudINS_14Kinemati cCloudISA_EEEEEEEEEEEEvRT_d+0x365) [0x4df385]
[linux-finq:07753] [ 9] MycoalChemistryFoam4(_ZN4Foam14KinematicCloudINS_5 CloudINS_24ReactingMultiphaseParcelINS_14ReactingP arcelINS_12ThermoParcelINS_15KinematicParcelINS_8p articleEEEEEEEEEEEE11evolveCloudINS9_12TrackingDat aINS_23ReactingMultiphaseCloudINS_13ReactingCloudI NS_11ThermoCloudISC_EEEEEEEEEEvRT_+0x276) [0x4e0396]
[linux-finq:07753] [10] MycoalChemistryFoam4(_ZN4Foam14KinematicCloudINS_5 CloudINS_24ReactingMultiphaseParcelINS_14ReactingP arcelINS_12ThermoParcelINS_15KinematicParcelINS_8p articleEEEEEEEEEEEE5solveINS9_12TrackingDataINS_23 ReactingMultiphaseCloudINS_13ReactingCloudINS_11Th ermoCloudISC_EEEEEEEEEEvRT_+0x2bd) [0x4e06fd]
[linux-finq:07753] [11] MycoalChemistryFoam4() [0x4428a5]
[linux-finq:07753] [12] /lib64/libc.so.6(__libc_start_main+0xf5) [0x7f75efdeabe5]
[linux-finq:07753] [13] MycoalChemistryFoam4() [0x44ce2d]
[linux-finq:07753] *** End of error message ***
--------------------------------------------------------------------------
mpirun noticed that process rank 5 with PID 7753 on node linux-finq exited on signal 8 (Floating point exception).
--------------------------------------------------------------------------
TPT-Sim7@linux-finq:~/OpenFOAM/TPT-Sim7-2.2.2/run/Funkenberechnung1/0_Datumsordner/04_29_40mJ_0,4mmS4_constcp,vielStaub_fastReaction_ neuePyrolyse> mpirun -np 12 MycoalChemistryFoam4 -parallel
/*---------------------------------------------------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.2.2 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
Build : 2.2.2
Exec : MycoalChemistryFoam4 -parallel
Date : Apr 30 2014
Time : 10:03:53
Host : "linux-finq"
PID : 9957
Case : /home/TPT-Sim7/OpenFOAM/TPT-Sim7-2.2.2/run/Funkenberechnung1/0_Datumsordner/04_29_40mJ_0,4mmS4_constcp,vielStaub_fastReaction_ neuePyrolyse
nProcs : 12
Slaves :
11
(
"linux-finq.9958"
"linux-finq.9959"
"linux-finq.9960"
"linux-finq.9961"
"linux-finq.9962"
"linux-finq.9963"
"linux-finq.9964"
"linux-finq.9965"
"linux-finq.9966"
"linux-finq.9967"
"linux-finq.9968"
)

Pstream initialized with:
floatTransfer : 0
nProcsSimpleSum : 0
commsType : nonBlocking
polling iterations : 0
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.0018


Reading g
Creating combustion model

Selecting combustion model infinitelyFastChemistry<psiThermoCombustion,gasHTh ermoPhysics>
Selecting thermodynamics package
{
type hePsiThermo;
mixture singleStepReactingMixture;
transport sutherland;
thermo janaf;
energy sensibleEnthalpy;
equationOfState perfectGas;
specie specie;
}

Selecting chemistryReader foamChemistryReader
Fuel heat of combustion :44634083.44
stoichiometric air-fuel ratio :3.532718617
stoichiometric oxygen-fuel ratio :3.532718617
Maximum products mass concentrations:
CO2: 0.6968982505
H2O: 0.3031017495
Combustion mode: explicit
Creating component thermo properties:
multi-component carrier - 7 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

No finite volume options present


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 3-D injection
Choosing nParticle to be a fixed value, massTotal variable now does not determine anything.
Selecting distribution model RosinRammler
Selecting dispersion model stochasticDispersionRAS
Selecting patch interaction model standardWallInteraction
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:
C16H34
Selecting devolatilisation model constantRateDevolatilisation
Participating volatile species:
CH4: particle mass fraction = 0
H2: particle mass fraction = 0
CO2: particle mass fraction = 0
Selecting surface reaction model COxidationKineticDiffusionLimitedRate
C(s): particle mass fraction = 0.431848

Constructing limestone cloud
employing parcel parcelTypeId of 2
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 3-D injection
Selecting distribution model RosinRammler
Selecting dispersion model stochasticDispersionRAS
Selecting patch interaction model standardWallInteraction
Selecting surface film model none
Selecting U integration scheme Euler
employing parcel parcelTypeId of 2
Selecting heat transfer model RanzMarshall
Selecting T integration scheme analytical
Selecting radiationModel none
Courant Number mean: 0.0002599751084 max: 0.01511134008

PIMPLE: Operating solver in PISO mode


Starting time loop

Courant Number mean: 0.0002599751084 max: 0.01511134008
deltaT = 1.99920032e-08
Time = 0.00180002


Solving 3-D cloud coalCloud1
[3] #0 Foam::error::printStack(Foam::Ostream&)[5] #0 [7] Foam::error::printStack(Foam::Ostream&)#0 Foam::error::printStack(Foam::Ostream&) at ??:?
[5] #1 Foam::sigFpe::sigHandler(int) at ??:?
[7] #1 Foam::sigFpe::sigHandler(int) at ??:?
[3] #1 Foam::sigFpe::sigHandler(int) at ??:?
[5] #2 at ??:?
[3] #2 at ??:?
[7] #2 in "/lib64/libc.so.6"
...


All times are GMT -4. The time now is 13:19.