CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Running, Solving & CFD

problem with solving lagrange reaction cloud

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Display Modes
Old   April 30, 2014, 07:53
Default problem with solving lagrange reaction cloud
  #1
New Member
 
Werner
Join Date: Apr 2014
Posts: 11
Rep Power: 3
Polli is on a distinguished road
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:rintStack(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:rintStack(Foam::Ostream&)[5] #0 [7] Foam::error:rintStack(Foam::Ostream&)#0 Foam::error:rintStack(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"
...
Polli is offline   Reply With Quote

Reply

Tags
cloud, combustion, lagrange particle

Thread Tools
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Cannot run the code properly: very large time step continuity error crst15 OpenFOAM Running, Solving & CFD 9 December 14, 2014 19:17
Simulation seems to converge but crashes suddenly xxxx OpenFOAM 16 September 12, 2014 08:07
transsonic nozzle with rhoSimpleFoam Unseen OpenFOAM Running, Solving & CFD 7 April 16, 2014 03:38
Extrusion with OpenFoam problem No. Iterations 0 Lord Kelvin OpenFOAM 6 April 12, 2011 11:24
Low Mach number Compressible jet flow using LES ankgupta8um OpenFOAM Running, Solving & CFD 7 January 15, 2011 14:38


All times are GMT -4. The time now is 22:04.