[Polli]

I found the following bug. If you have an particle which consists of ash, C, and burnable liquids like C16H34, sometimes the simulation breaks down when all the C16H34 or other liquids are away. When you stop the simulation it is not possible to restart it from that point when the liquids are away. But also in time of simulation sometimes it breaks down on any timestep and it is not possible to restart it. Do anybody know how to solve this problem? I think it is someting in the devolatiliaton part. Here the console output:

Solving 3-D cloud coalCloud1
Cloud: coalCloud1
Current number of parcels = 200
Current mass in system = 3.557672736e-07
Linear momentum = (-7.131925902e-08 -2.35527338e-08 5.219148202e-09)
|Linear momentum| = 7.528882709e-08
Linear kinetic energy = 4.973265222e-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 = 305.1890576, 1414.751503
Mass transfer phase change = 1.328956196e-08
Mass transfer devolatilisation = 0
Mass transfer surface reaction = 1.14231585e-11


Solving 3-D cloud limestoneCloud1
Cloud: limestoneCloud1
Current number of parcels = 0
Current mass in system = 0
Linear momentum = (0 0 0)
|Linear momentum| = 0
Linear kinetic energy = 0
Rotational kinetic energy = 0
model1:
number of parcels added = 0
mass introduced = 0
Parcels absorbed into film = 0
New film detached parcels = 0
Parcel fate (number, mass)
- escape = 0, 0
- stick = 0, 0
Temperature min/max = 0, 0

diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCG: Solving for Ux, Initial residual = 0.004174570485, Final residual = 6.55497394e-26, No Iterations 48
DILUPBiCG: Solving for Uy, Initial residual = 0.004820024763, Final residual = 7.073102926e-26, No Iterations 49
DILUPBiCG: Solving for Uz, Initial residual = 0.007326765268, Final residual = 7.541125317e-26, No Iterations 49
DILUPBiCG: Solving for CH4, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCG: Solving for C16H34, Initial residual = 4.72855249e-05, Final residual = 3.048644737e-10, No Iterations 2
DILUPBiCG: Solving for O2, Initial residual = 2.705945236e-05, Final residual = 2.948377816e-09, No Iterations 3
DILUPBiCG: Solving for H2, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCG: Solving for CO2, Initial residual = 2.872733169e-05, Final residual = 9.596927472e-09, No Iterations 2
DILUPBiCG: Solving for H2O, Initial residual = 2.870491424e-05, Final residual = 9.555075994e-09, No Iterations 2
DILUPBiCG: Solving for h, Initial residual = 3.50509694e-05, Final residual = 9.366565726e-26, No Iterations 45
T gas min/max = 311.1077118, 3039.43868
T gas min/max = 529.5399963, 937.6982209
T gas min/max = 311.1114847, 3039.661303
GAMG: Solving for p, Initial residual = 0.03405915517, Final residual = 1.281626863e-05, No Iterations 2
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 4.096748315e-08, global = 2.88424725e-10, cumulative = -9.818412466e-06
GAMG: Solving for p, Initial residual = 1.63088632e-05, Final residual = 9.331245638e-08, No Iterations 15
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 3.015852833e-10, global = 2.112726392e-12, cumulative = -9.818410353e-06
GAMG: Solving for p, Initial residual = 7.457479388e-06, Final residual = 9.592142917e-09, No Iterations 24
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 3.087218428e-11, global = -2.364163373e-13, cumulative = -9.81841059e-06
DILUPBiCG: Solving for epsilon, Initial residual = 4.928435687e-05, Final residual = 3.629117024e-26, No Iterations 48
DILUPBiCG: Solving for k, Initial residual = 0.0003198610339, Final residual = 9.268242123e-26, No Iterations 48
ExecutionTime = 20543.52 s ClockTime = 24489 s

Courant Number mean: 0.0008461850078 max: 0.1804167747
deltaT = 1.664381602e-07
Time = 0.00188419


Solving 3-D cloud coalCloud1
[3] #0 Foam::error::printStack(Foam::Ostream&) at ??:?
[3] #1 Foam::sigFpe::sigHandler(int) at ??:?
[3] #2 in "/lib64/libc.so.6"
[3] #3 Foam::liquidMixtureProperties::X(Foam::Field<doubl e> const&) const at ??:?
[3] #4 Foam::LiquidEvaporation<Foam::ReactingCloud<Foam:: ThermoCloud<Foam::KinematicCloud<Foam::Cloud<Foam: :ReactingMultiphaseParcel<Foam::ReactingParcel<Foa m::ThermoParcel<Foam::KinematicParcel<Foam::partic le> > > > > > > > >::Tvap(Foam::Field<double> const&) const at ??:?
[3] #5 void Foam::ReactingParcel<Foam::ThermoParcel<Foam::Kine maticParcel<Foam::particle> > >::calcPhaseChange<Foam::ReactingParcel<Foam::Ther moParcel<Foam::KinematicParcel<Foam::particle> > >::TrackingData<Foam::ReactingMultiphaseCloud<Foam ::ReactingCloud<Foam::ThermoCloud<Foam::KinematicC loud<Foam::Cloud<Foam::ReactingMultiphaseParcel<Fo am::ReactingParcel<Foam::ThermoParcel<Foam::Kinema ticParcel<Foam::particle> > > > > > > > > > >(Foam::ReactingParcel<Foam::ThermoParcel<Foam::Ki nematicParcel<Foam::particle> > >::TrackingData<Foam::ReactingMultiphaseCloud<Foam ::ReactingCloud<Foam::ThermoCloud<Foam::KinematicC loud<Foam::Cloud<Foam::ReactingMultiphaseParcel<Fo am::ReactingParcel<Foam::ThermoParcel<Foam::Kinema ticParcel<Foam::particle> > > > > > > > > >&, 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 ??:?
[3] #6 void Foam::ReactingMultiphaseParcel<Foam::ReactingParce l<Foam::ThermoParcel<Foam::KinematicParcel<Foam::p article> > > >::calc<Foam::ReactingParcel<Foam::ThermoParcel<Fo am::KinematicParcel<Foam::particle> > >::TrackingData<Foam::ReactingMultiphaseCloud<Foam ::ReactingCloud<Foam::ThermoCloud<Foam::KinematicC loud<Foam::Cloud<Foam::ReactingMultiphaseParcel<Fo am::ReactingParcel<Foam::ThermoParcel<Foam::Kinema ticParcel<Foam::particle> > > > > > > > > > >(Foam::ReactingParcel<Foam::ThermoParcel<Foam::Ki nematicParcel<Foam::particle> > >::TrackingData<Foam::ReactingMultiphaseCloud<Foam ::ReactingCloud<Foam::ThermoCloud<Foam::KinematicC loud<Foam::Cloud<Foam::ReactingMultiphaseParcel<Fo am::ReactingParcel<Foam::ThermoParcel<Foam::Kinema ticParcel<Foam::particle> > > > > > > > > >&, double, int) at ??:?
[3] #7 bool Foam::KinematicParcel<Foam::particle>::move<Foam:: ReactingParcel<Foam::ThermoParcel<Foam::KinematicP arcel<Foam::particle> > >::TrackingData<Foam::ReactingMultiphaseCloud<Foam ::ReactingCloud<Foam::ThermoCloud<Foam::KinematicC loud<Foam::Cloud<Foam::ReactingMultiphaseParcel<Fo am::ReactingParcel<Foam::ThermoParcel<Foam::Kinema ticParcel<Foam::particle> > > > > > > > > > >(Foam::ReactingParcel<Foam::ThermoParcel<Foam::Ki nematicParcel<Foam::particle> > >::TrackingData<Foam::ReactingMultiphaseCloud<Foam ::ReactingCloud<Foam::ThermoCloud<Foam::KinematicC loud<Foam::Cloud<Foam::ReactingMultiphaseParcel<Fo am::ReactingParcel<Foam::ThermoParcel<Foam::Kinema ticParcel<Foam::particle> > > > > > > > > >&, double) at ??:?
[3] #8 void Foam::Cloud<Foam::ReactingMultiphaseParcel<Foam::R eactingParcel<Foam::ThermoParcel<Foam::KinematicPa rcel<Foam::particle> > > > >::move<Foam::ReactingParcel<Foam::ThermoParcel<Fo am::KinematicParcel<Foam::particle> > >::TrackingData<Foam::ReactingMultiphaseCloud<Foam ::ReactingCloud<Foam::ThermoCloud<Foam::KinematicC loud<Foam::Cloud<Foam::ReactingMultiphaseParcel<Fo am::ReactingParcel<Foam::ThermoParcel<Foam::Kinema ticParcel<Foam::particle> > > > > > > > > > >(Foam::ReactingParcel<Foam::ThermoParcel<Foam::Ki nematicParcel<Foam::particle> > >::TrackingData<Foam::ReactingMultiphaseCloud<Foam ::ReactingCloud<Foam::ThermoCloud<Foam::KinematicC loud<Foam::Cloud<Foam::ReactingMultiphaseParcel<Fo am::ReactingParcel<Foam::ThermoParcel<Foam::Kinema ticParcel<Foam::particle> > > > > > > > > >&, double) at ??:?
[3] #9 void Foam::KinematicCloud<Foam::Cloud<Foam::ReactingMul tiphaseParcel<Foam::ReactingParcel<Foam::ThermoPar cel<Foam::KinematicParcel<Foam::particle> > > > > >::evolveCloud<Foam::ReactingParcel<Foam::ThermoPa rcel<Foam::KinematicParcel<Foam::particle> > >::TrackingData<Foam::ReactingMultiphaseCloud<Foam ::ReactingCloud<Foam::ThermoCloud<Foam::KinematicC loud<Foam::Cloud<Foam::ReactingMultiphaseParcel<Fo am::ReactingParcel<Foam::ThermoParcel<Foam::Kinema ticParcel<Foam::particle> > > > > > > > > > >(Foam::ReactingParcel<Foam::ThermoParcel<Foam::Ki nematicParcel<Foam::particle> > >::TrackingData<Foam::ReactingMultiphaseCloud<Foam ::ReactingCloud<Foam::ThermoCloud<Foam::KinematicC loud<Foam::Cloud<Foam::ReactingMultiphaseParcel<Fo am::ReactingParcel<Foam::ThermoParcel<Foam::Kinema ticParcel<Foam::particle> > > > > > > > > >&) at ??:?
[3] #10 void Foam::KinematicCloud<Foam::Cloud<Foam::ReactingMul tiphaseParcel<Foam::ReactingParcel<Foam::ThermoPar cel<Foam::KinematicParcel<Foam::particle> > > > > >::solve<Foam::ReactingParcel<Foam::ThermoParcel<F oam::KinematicParcel<Foam::particle> > >::TrackingData<Foam::ReactingMultiphaseCloud<Foam ::ReactingCloud<Foam::ThermoCloud<Foam::KinematicC loud<Foam::Cloud<Foam::ReactingMultiphaseParcel<Fo am::ReactingParcel<Foam::ThermoParcel<Foam::Kinema ticParcel<Foam::particle> > > > > > > > > > >(Foam::ReactingParcel<Foam::ThermoParcel<Foam::Ki nematicParcel<Foam::particle> > >::TrackingData<Foam::ReactingMultiphaseCloud<Foam ::ReactingCloud<Foam::ThermoCloud<Foam::KinematicC loud<Foam::Cloud<Foam::ReactingMultiphaseParcel<Fo am::ReactingParcel<Foam::ThermoParcel<Foam::Kinema ticParcel<Foam::particle> > > > > > > > > >&) at ??:?
[3] #11
[3] at ??:?
[3] #12 __libc_start_main in "/lib64/libc.so.6"
[3] #13
[3] at /home/abuild/rpmbuild/BUILD/glibc-2.18/csu/../sysdeps/x86_64/start.S:125
[linux-finq:22313] *** Process received signal ***
[linux-finq:22313] Signal: Floating point exception (8)
[linux-finq:22313] Signal code: (-6)
[linux-finq:22313] Failing at address: 0x3e800005729
[linux-finq:22313] [ 0] /lib64/libc.so.6(+0x358d0) [0x7f459e1258d0]
[linux-finq:22313] [ 1] /lib64/libc.so.6(gsignal+0x39) [0x7f459e125849]
[linux-finq:22313] [ 2] /lib64/libc.so.6(+0x358d0) [0x7f459e1258d0]
[linux-finq:22313] [ 3] /opt/OpenFOAM-2.2.2/platforms/linux64GccDPOpt/lib/libliquidMixtureProperties.so(_ZNK4Foam23liquidMix tureProperties1XERKNS_5FieldIdEE+0x243) [0x7f45a2300c63]
[linux-finq:22313] [ 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) [0x7f45a382a6bc]
[linux-finq:22313] [ 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:22313] [ 6] MycoalChemistryFoam4(_ZN4Foam24ReactingMultiphaseP arcelINS_14ReactingParcelINS_12ThermoParcelINS_15K inematicParcelINS_8particleEEEEEEEE4calcINS7_12Tra ckingDataINS_23ReactingMultiphaseCloudINS_13Reacti ngCloudINS_11ThermoCloudINS_14KinematicCloudINS_5C loudIS8_EEEEEEEEEEEEEEvRT_di+0x399) [0x4a2179]
[linux-finq:22313] [ 7] MycoalChemistryFoam4(_ZN4Foam15KinematicParcelINS_ 8particleEE4moveINS_14ReactingParcelINS_12ThermoPa rcelIS2_EEE12TrackingDataINS_23ReactingMultiphaseC loudINS_13ReactingCloudINS_11ThermoCloudINS_14Kine maticCloudINS_5CloudINS_24ReactingMultiphaseParcel IS7_EEEEEEEEEEEEEEEEbRT_d+0x3d0) [0x4d9b50]
[linux-finq:22313] [ 8] MycoalChemistryFoam4(_ZN4Foam5CloudINS_24ReactingM ultiphaseParcelINS_14ReactingParcelINS_12ThermoPar celINS_15KinematicParcelINS_8particleEEEEEEEEEE4mo veINS8_12TrackingDataINS_23ReactingMultiphaseCloud INS_13ReactingCloudINS_11ThermoCloudINS_14Kinemati cCloudISA_EEEEEEEEEEEEvRT_d+0x365) [0x4df385]
[linux-finq:22313] [ 9] MycoalChemistryFoam4(_ZN4Foam14KinematicCloudINS_5 CloudINS_24ReactingMultiphaseParcelINS_14ReactingP arcelINS_12ThermoParcelINS_15KinematicParcelINS_8p articleEEEEEEEEEEEE11evolveCloudINS9_12TrackingDat aINS_23ReactingMultiphaseCloudINS_13ReactingCloudI NS_11ThermoCloudISC_EEEEEEEEEEvRT_+0x276) [0x4e0396]
[linux-finq:22313] [10] MycoalChemistryFoam4(_ZN4Foam14KinematicCloudINS_5 CloudINS_24ReactingMultiphaseParcelINS_14ReactingP arcelINS_12ThermoParcelINS_15KinematicParcelINS_8p articleEEEEEEEEEEEE5solveINS9_12TrackingDataINS_23 ReactingMultiphaseCloudINS_13ReactingCloudINS_11Th ermoCloudISC_EEEEEEEEEEvRT_+0x2bd) [0x4e06fd]
[linux-finq:22313] [11] MycoalChemistryFoam4() [0x4428a5]
[linux-finq:22313] [12] /lib64/libc.so.6(__libc_start_main+0xf5) [0x7f459e111be5]
[linux-finq:22313] [13] MycoalChemistryFoam4() [0x44ce2d]
[linux-finq:22313] *** End of error message ***
--------------------------------------------------------------------------
mpirun noticed that process rank 3 with PID 22313 on node linux-finq exited on signal 8 (Floating point exception).