Quote:
Originally Posted by lukasfischer
(Post 354146)
hi forum users,
i am trying to model a nebulizer by injecting a spray source in a gas flow. my simulation crashes randomly while solving the cloud, i think the cause is that single parcels become colder than the temperature set as TMin=... in constantProperties in sprayCloudProperties.
when i restart the simulation it will continue, probably because of the small initial timestep in controlDict.
do you think that this could be the right explaination, if yes, is there a way to control the time step like with the maxCo number, but for the termal calculation of the cloud? Or could i set dummy thermal properties for the H2O
to do nothing or at least to not cause a FPE?
attached is a plot of the last timestep, the parcels are colored by temp, one of the spheres is at 150K which is the value i set as TMin. the slice is colored by pressure.
thanks, Lukas
Code:
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 9.48149e10, global = 1.20385e11, cumulative = 1.17464e08
DILUPBiCG: Solving for epsilon, Initial residual = 1.06239e05, Final residual = 4.85898e07, No Iterations 1
DILUPBiCG: Solving for k, Initial residual = 1.54945e05, Final residual = 1.11891e06, No Iterations 1
ExecutionTime = 806.61 s ClockTime = 916 s
Courant Number mean: 0.0209411 max: 0.599903
deltaT = 0.000481764
Time = 2.85017
Solving cloud sprayCloud
> Cloud: sprayCloud
Added 10 new parcels
Cloud: sprayCloud
Current number of parcels = 687
Current mass in system = 3.27202e09
Linear momentum = (1.2832e12 3.82908e12 1.3838e10)
Linear momentum = 1.38439e10
Linear kinetic energy = 2.35283e11
Rotational kinetic energy = 0
Total number of parcels added = 66196
Total mass introduced = 6.15151e07
Parcel fates:
 escape = 0, 0
 stick = 0, 0
Mass transfer phase change = 3.88596e07
D32 (mu) = 9.88594
Liquid penetration 95% mass (m) = 0.0168679
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCG: Solving for Ux, Initial residual = 0.000911625, Final residual = 2.26302e06, No Iterations 3
DILUPBiCG: Solving for Uy, Initial residual = 0.000594203, Final residual = 5.52574e06, No Iterations 2
DILUPBiCG: Solving for Uz, Initial residual = 0.000514736, Final residual = 2.54549e06, No Iterations 3
DILUPBiCG: Solving for C7H16, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCG: Solving for O2, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCG: Solving for CO2, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCG: Solving for H2O, Initial residual = 0.000472812, Final residual = 6.91721e06, No Iterations 2
DILUPBiCG: Solving for hs, Initial residual = 0.000381923, Final residual = 3.39677e06, No Iterations 3
T gas min/max = 322.537, 370
DICPCG: Solving for p, Initial residual = 0.00858937, Final residual = 0.000689481, No Iterations 74
DICPCG: Solving for p, Initial residual = 0.00152945, Final residual = 0.000148412, No Iterations 7
DICPCG: Solving for p, Initial residual = 0.000221837, Final residual = 2.20728e05, No Iterations 10
DICPCG: Solving for p, Initial residual = 3.29432e05, Final residual = 2.72098e06, No Iterations 70
DICPCG: Solving for p, Initial residual = 4.25008e06, Final residual = 9.68375e07, No Iterations 4
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 1.49365e09, global = 1.20547e11, cumulative = 1.17344e08
DICPCG: Solving for p, Initial residual = 0.00196189, Final residual = 0.000190001, No Iterations 26
DICPCG: Solving for p, Initial residual = 0.00049205, Final residual = 4.17605e05, No Iterations 59
DICPCG: Solving for p, Initial residual = 7.16986e05, Final residual = 6.93027e06, No Iterations 19
DICPCG: Solving for p, Initial residual = 1.03981e05, Final residual = 1.03773e06, No Iterations 14
DICPCG: Solving for p, Initial residual = 1.55307e06, Final residual = 8.79028e07, No Iterations 1
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 1.35458e09, global = 6.85588e10, cumulative = 1.10488e08
DICPCG: Solving for p, Initial residual = 0.000217007, Final residual = 2.06744e05, No Iterations 12
DICPCG: Solving for p, Initial residual = 3.76131e05, Final residual = 2.93503e06, No Iterations 71
DICPCG: Solving for p, Initial residual = 4.70404e06, Final residual = 9.41462e07, No Iterations 4
DICPCG: Solving for p, Initial residual = 1.03837e06, Final residual = 7.77936e07, No Iterations 1
DICPCG: Solving for p, Initial residual = 7.8544e07, Final residual = 7.8544e07, No Iterations 0
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 1.21031e09, global = 1.17959e11, cumulative = 1.1037e08
DICPCG: Solving for p, Initial residual = 4.69634e05, Final residual = 4.51442e06, No Iterations 17
DICPCG: Solving for p, Initial residual = 7.84907e06, Final residual = 9.82477e07, No Iterations 65
DICPCG: Solving for p, Initial residual = 1.31751e06, Final residual = 6.55612e07, No Iterations 1
DICPCG: Solving for p, Initial residual = 6.63309e07, Final residual = 6.63309e07, No Iterations 0
DICPCG: Solving for p, Initial residual = 6.63309e07, Final residual = 6.63309e07, No Iterations 0
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 1.02209e09, global = 2.36106e11, cumulative = 1.10606e08
DILUPBiCG: Solving for epsilon, Initial residual = 1.15936e05, Final residual = 4.86978e07, No Iterations 1
DILUPBiCG: Solving for k, Initial residual = 1.76057e05, Final residual = 1.14782e06, No Iterations 1
ExecutionTime = 806.99 s ClockTime = 916 s
Courant Number mean: 0.0209503 max: 0.600233
deltaT = 0.000481576
Time = 2.85065
Solving cloud sprayCloud
> Cloud: sprayCloud
Added 16 new parcels
[3] #0 Foam::error::printStack(Foam::Ostream&) in "/opt/openfoam210/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[3] #1 Foam::sigFpe::sigHandler(int) in "/opt/openfoam210/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[3] #2 in "/lib/x86_64linuxgnu/libc.so.6"
[3] #3 Foam::H2O::D(double, double, double) const in "/opt/openfoam210/platforms/linux64GccDPOpt/lib/libliquidProperties.so"
[3] #4
[3] in "/opt/openfoam210/platforms/linux64GccDPOpt/bin/sprayFoam"
[3] #5
[3] in "/opt/openfoam210/platforms/linux64GccDPOpt/bin/sprayFoam"
[3] #6
[3] in "/opt/openfoam210/platforms/linux64GccDPOpt/bin/sprayFoam"
[3] #7
[3] in "/opt/openfoam210/platforms/linux64GccDPOpt/bin/sprayFoam"
[3] #8
[3] in "/opt/openfoam210/platforms/linux64GccDPOpt/bin/sprayFoam"
[3] #9
[3] in "/opt/openfoam210/platforms/linux64GccDPOpt/bin/sprayFoam"
[3] #10
[3] in "/opt/openfoam210/platforms/linux64GccDPOpt/bin/sprayFoam"
[3] #11 __libc_start_main in "/lib/x86_64linuxgnu/libc.so.6"
[3] #12
[3] in "/opt/openfoam210/platforms/linux64GccDPOpt/bin/sprayFoam"
[lukaslaptop:02279] *** Process received signal ***
[lukaslaptop:02279] Signal: Floating point exception (8)
[lukaslaptop:02279] Signal code: (6)
[lukaslaptop:02279] Failing at address: 0x3e8000008e7
[lukaslaptop:02279] [ 0] /lib/x86_64linuxgnu/libc.so.6(+0x33d80) [0x7f2c47a0cd80]
[lukaslaptop:02279] [ 1] /lib/x86_64linuxgnu/libc.so.6(gsignal+0x35) [0x7f2c47a0cd05]
[lukaslaptop:02279] [ 2] /lib/x86_64linuxgnu/libc.so.6(+0x33d80) [0x7f2c47a0cd80]
[lukaslaptop:02279] [ 3] /opt/openfoam210/platforms/linux64GccDPOpt/lib/libliquidProperties.so(_ZNK4Foam3H2O1DEddd+0x2a) [0x7f2c4ad9c7ea]
[lukaslaptop:02279] [ 4] sprayFoam() [0x479f4d]
[lukaslaptop:02279] [ 5] sprayFoam() [0x4a582e]
[lukaslaptop:02279] [ 6] sprayFoam() [0x4a5f30]
[lukaslaptop:02279] [ 7] sprayFoam() [0x4a710f]
[lukaslaptop:02279] [ 8] sprayFoam() [0x4a7de1]
[lukaslaptop:02279] [ 9] sprayFoam() [0x4a89a3]
[lukaslaptop:02279] [10] sprayFoam() [0x431558]
[lukaslaptop:02279] [11] /lib/x86_64linuxgnu/libc.so.6(__libc_start_main+0xff) [0x7f2c479f7eff]
[lukaslaptop:02279] [12] sprayFoam() [0x424d89]
[lukaslaptop:02279] *** End of error message ***

mpirun noticed that process rank 3 with PID 2279 on node lukaslaptop exited on signal 8 (Floating point exception).

Hi Lukas,
have you solved your problem?
The similar problem occured to me.
The error message is shown as:

Starting time loop
Courant Number mean: 0 max: 0
Time = 1e05
Solving cloud coalCloud1
> Cloud: coalCloud1
Added 10 new parcels
#0 Foam::error::printStack(Foam::Ostream&) in "/home/oliver/OpenFOAM/OpenFOAM2.1.0/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#1 Foam::sigFpe::sigHandler(int) in "/home/oliver/OpenFOAM/OpenFOAM2.1.0/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#2
at sigaction.c:0
#3 Foam::H2O::D(double, double, double) const in "/home/oliver/OpenFOAM/OpenFOAM2.1.0/platforms/linux64GccDPOpt/lib/libliquidProperties.so"
#4

The case is modified from the original one, and only mesh was changed.....
so it will be highly appreciated if you can give me some suggestions on this problem.
Best regards.
oliver
