CFD Online Discussion Forums

CFD Online Discussion Forums (http://www.cfd-online.com/Forums/)
-   OpenFOAM Running, Solving & CFD (http://www.cfd-online.com/Forums/openfoam-solving/)
-   -   sprayFoam crashes (http://www.cfd-online.com/Forums/openfoam-solving/99766-sprayfoam-crashes.html)

lukasfischer April 11, 2012 06:15

sprayFoam crashes
 
1 Attachment(s)
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.48149e-10, global = 1.20385e-11, cumulative = 1.17464e-08
DILUPBiCG:  Solving for epsilon, Initial residual = 1.06239e-05, Final residual = 4.85898e-07, No Iterations 1
DILUPBiCG:  Solving for k, Initial residual = 1.54945e-05, Final residual = 1.11891e-06, 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.27202e-09
    Linear momentum                = (1.2832e-12 3.82908e-12 1.3838e-10)
  |Linear momentum|                = 1.38439e-10
    Linear kinetic energy          = 2.35283e-11
    Rotational kinetic energy      = 0
    Total number of parcels added  = 66196
    Total mass introduced          = 6.15151e-07
    Parcel fates:
      - escape                      = 0, 0
      - stick                      = 0, 0
    Mass transfer phase change      = 3.88596e-07
    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.26302e-06, No Iterations 3
DILUPBiCG:  Solving for Uy, Initial residual = 0.000594203, Final residual = 5.52574e-06, No Iterations 2
DILUPBiCG:  Solving for Uz, Initial residual = 0.000514736, Final residual = 2.54549e-06, 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.91721e-06, No Iterations 2
DILUPBiCG:  Solving for hs, Initial residual = 0.000381923, Final residual = 3.39677e-06, 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.20728e-05, No Iterations 10
DICPCG:  Solving for p, Initial residual = 3.29432e-05, Final residual = 2.72098e-06, No Iterations 70
DICPCG:  Solving for p, Initial residual = 4.25008e-06, Final residual = 9.68375e-07, No Iterations 4
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 1.49365e-09, global = -1.20547e-11, cumulative = 1.17344e-08
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.17605e-05, No Iterations 59
DICPCG:  Solving for p, Initial residual = 7.16986e-05, Final residual = 6.93027e-06, No Iterations 19
DICPCG:  Solving for p, Initial residual = 1.03981e-05, Final residual = 1.03773e-06, No Iterations 14
DICPCG:  Solving for p, Initial residual = 1.55307e-06, Final residual = 8.79028e-07, No Iterations 1
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 1.35458e-09, global = -6.85588e-10, cumulative = 1.10488e-08
DICPCG:  Solving for p, Initial residual = 0.000217007, Final residual = 2.06744e-05, No Iterations 12
DICPCG:  Solving for p, Initial residual = 3.76131e-05, Final residual = 2.93503e-06, No Iterations 71
DICPCG:  Solving for p, Initial residual = 4.70404e-06, Final residual = 9.41462e-07, No Iterations 4
DICPCG:  Solving for p, Initial residual = 1.03837e-06, Final residual = 7.77936e-07, No Iterations 1
DICPCG:  Solving for p, Initial residual = 7.8544e-07, Final residual = 7.8544e-07, No Iterations 0
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 1.21031e-09, global = -1.17959e-11, cumulative = 1.1037e-08
DICPCG:  Solving for p, Initial residual = 4.69634e-05, Final residual = 4.51442e-06, No Iterations 17
DICPCG:  Solving for p, Initial residual = 7.84907e-06, Final residual = 9.82477e-07, No Iterations 65
DICPCG:  Solving for p, Initial residual = 1.31751e-06, Final residual = 6.55612e-07, No Iterations 1
DICPCG:  Solving for p, Initial residual = 6.63309e-07, Final residual = 6.63309e-07, No Iterations 0
DICPCG:  Solving for p, Initial residual = 6.63309e-07, Final residual = 6.63309e-07, No Iterations 0
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 1.02209e-09, global = 2.36106e-11, cumulative = 1.10606e-08
DILUPBiCG:  Solving for epsilon, Initial residual = 1.15936e-05, Final residual = 4.86978e-07, No Iterations 1
DILUPBiCG:  Solving for k, Initial residual = 1.76057e-05, Final residual = 1.14782e-06, 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_64-linux-gnu/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_64-linux-gnu/libc.so.6"
[3] #12 
[3]  in "/opt/openfoam210/platforms/linux64GccDPOpt/bin/sprayFoam"
[lukas-laptop:02279] *** Process received signal ***
[lukas-laptop:02279] Signal: Floating point exception (8)
[lukas-laptop:02279] Signal code:  (-6)
[lukas-laptop:02279] Failing at address: 0x3e8000008e7
[lukas-laptop:02279] [ 0] /lib/x86_64-linux-gnu/libc.so.6(+0x33d80) [0x7f2c47a0cd80]
[lukas-laptop:02279] [ 1] /lib/x86_64-linux-gnu/libc.so.6(gsignal+0x35) [0x7f2c47a0cd05]
[lukas-laptop:02279] [ 2] /lib/x86_64-linux-gnu/libc.so.6(+0x33d80) [0x7f2c47a0cd80]
[lukas-laptop:02279] [ 3] /opt/openfoam210/platforms/linux64GccDPOpt/lib/libliquidProperties.so(_ZNK4Foam3H2O1DEddd+0x2a) [0x7f2c4ad9c7ea]
[lukas-laptop:02279] [ 4] sprayFoam() [0x479f4d]
[lukas-laptop:02279] [ 5] sprayFoam() [0x4a582e]
[lukas-laptop:02279] [ 6] sprayFoam() [0x4a5f30]
[lukas-laptop:02279] [ 7] sprayFoam() [0x4a710f]
[lukas-laptop:02279] [ 8] sprayFoam() [0x4a7de1]
[lukas-laptop:02279] [ 9] sprayFoam() [0x4a89a3]
[lukas-laptop:02279] [10] sprayFoam() [0x431558]
[lukas-laptop:02279] [11] /lib/x86_64-linux-gnu/libc.so.6(__libc_start_main+0xff) [0x7f2c479f7eff]
[lukas-laptop:02279] [12] sprayFoam() [0x424d89]
[lukas-laptop:02279] *** End of error message ***
--------------------------------------------------------------------------
mpirun noticed that process rank 3 with PID 2279 on node lukas-laptop exited on signal 8 (Floating point exception).


MPJ May 10, 2012 08:34

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.48149e-10, global = 1.20385e-11, cumulative = 1.17464e-08
DILUPBiCG:  Solving for epsilon, Initial residual = 1.06239e-05, Final residual = 4.85898e-07, No Iterations 1
DILUPBiCG:  Solving for k, Initial residual = 1.54945e-05, Final residual = 1.11891e-06, 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.27202e-09
    Linear momentum                = (1.2832e-12 3.82908e-12 1.3838e-10)
  |Linear momentum|                = 1.38439e-10
    Linear kinetic energy          = 2.35283e-11
    Rotational kinetic energy      = 0
    Total number of parcels added  = 66196
    Total mass introduced          = 6.15151e-07
    Parcel fates:
      - escape                      = 0, 0
      - stick                      = 0, 0
    Mass transfer phase change      = 3.88596e-07
    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.26302e-06, No Iterations 3
DILUPBiCG:  Solving for Uy, Initial residual = 0.000594203, Final residual = 5.52574e-06, No Iterations 2
DILUPBiCG:  Solving for Uz, Initial residual = 0.000514736, Final residual = 2.54549e-06, 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.91721e-06, No Iterations 2
DILUPBiCG:  Solving for hs, Initial residual = 0.000381923, Final residual = 3.39677e-06, 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.20728e-05, No Iterations 10
DICPCG:  Solving for p, Initial residual = 3.29432e-05, Final residual = 2.72098e-06, No Iterations 70
DICPCG:  Solving for p, Initial residual = 4.25008e-06, Final residual = 9.68375e-07, No Iterations 4
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 1.49365e-09, global = -1.20547e-11, cumulative = 1.17344e-08
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.17605e-05, No Iterations 59
DICPCG:  Solving for p, Initial residual = 7.16986e-05, Final residual = 6.93027e-06, No Iterations 19
DICPCG:  Solving for p, Initial residual = 1.03981e-05, Final residual = 1.03773e-06, No Iterations 14
DICPCG:  Solving for p, Initial residual = 1.55307e-06, Final residual = 8.79028e-07, No Iterations 1
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 1.35458e-09, global = -6.85588e-10, cumulative = 1.10488e-08
DICPCG:  Solving for p, Initial residual = 0.000217007, Final residual = 2.06744e-05, No Iterations 12
DICPCG:  Solving for p, Initial residual = 3.76131e-05, Final residual = 2.93503e-06, No Iterations 71
DICPCG:  Solving for p, Initial residual = 4.70404e-06, Final residual = 9.41462e-07, No Iterations 4
DICPCG:  Solving for p, Initial residual = 1.03837e-06, Final residual = 7.77936e-07, No Iterations 1
DICPCG:  Solving for p, Initial residual = 7.8544e-07, Final residual = 7.8544e-07, No Iterations 0
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 1.21031e-09, global = -1.17959e-11, cumulative = 1.1037e-08
DICPCG:  Solving for p, Initial residual = 4.69634e-05, Final residual = 4.51442e-06, No Iterations 17
DICPCG:  Solving for p, Initial residual = 7.84907e-06, Final residual = 9.82477e-07, No Iterations 65
DICPCG:  Solving for p, Initial residual = 1.31751e-06, Final residual = 6.55612e-07, No Iterations 1
DICPCG:  Solving for p, Initial residual = 6.63309e-07, Final residual = 6.63309e-07, No Iterations 0
DICPCG:  Solving for p, Initial residual = 6.63309e-07, Final residual = 6.63309e-07, No Iterations 0
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 1.02209e-09, global = 2.36106e-11, cumulative = 1.10606e-08
DILUPBiCG:  Solving for epsilon, Initial residual = 1.15936e-05, Final residual = 4.86978e-07, No Iterations 1
DILUPBiCG:  Solving for k, Initial residual = 1.76057e-05, Final residual = 1.14782e-06, 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_64-linux-gnu/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_64-linux-gnu/libc.so.6"
[3] #12 
[3]  in "/opt/openfoam210/platforms/linux64GccDPOpt/bin/sprayFoam"
[lukas-laptop:02279] *** Process received signal ***
[lukas-laptop:02279] Signal: Floating point exception (8)
[lukas-laptop:02279] Signal code:  (-6)
[lukas-laptop:02279] Failing at address: 0x3e8000008e7
[lukas-laptop:02279] [ 0] /lib/x86_64-linux-gnu/libc.so.6(+0x33d80) [0x7f2c47a0cd80]
[lukas-laptop:02279] [ 1] /lib/x86_64-linux-gnu/libc.so.6(gsignal+0x35) [0x7f2c47a0cd05]
[lukas-laptop:02279] [ 2] /lib/x86_64-linux-gnu/libc.so.6(+0x33d80) [0x7f2c47a0cd80]
[lukas-laptop:02279] [ 3] /opt/openfoam210/platforms/linux64GccDPOpt/lib/libliquidProperties.so(_ZNK4Foam3H2O1DEddd+0x2a) [0x7f2c4ad9c7ea]
[lukas-laptop:02279] [ 4] sprayFoam() [0x479f4d]
[lukas-laptop:02279] [ 5] sprayFoam() [0x4a582e]
[lukas-laptop:02279] [ 6] sprayFoam() [0x4a5f30]
[lukas-laptop:02279] [ 7] sprayFoam() [0x4a710f]
[lukas-laptop:02279] [ 8] sprayFoam() [0x4a7de1]
[lukas-laptop:02279] [ 9] sprayFoam() [0x4a89a3]
[lukas-laptop:02279] [10] sprayFoam() [0x431558]
[lukas-laptop:02279] [11] /lib/x86_64-linux-gnu/libc.so.6(__libc_start_main+0xff) [0x7f2c479f7eff]
[lukas-laptop:02279] [12] sprayFoam() [0x424d89]
[lukas-laptop:02279] *** End of error message ***
--------------------------------------------------------------------------
mpirun noticed that process rank 3 with PID 2279 on node lukas-laptop 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 = 1e-05


Solving cloud coalCloud1

--> Cloud: coalCloud1
Added 10 new parcels

#0 Foam::error::printStack(Foam::Ostream&) in "/home/oliver/OpenFOAM/OpenFOAM-2.1.0/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#1 Foam::sigFpe::sigHandler(int) in "/home/oliver/OpenFOAM/OpenFOAM-2.1.0/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#2
at sigaction.c:0
#3 Foam::H2O::D(double, double, double) const in "/home/oliver/OpenFOAM/OpenFOAM-2.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

s1300044 July 12, 2013 23:14

Hi, Lukas,
I would like to simulate a system like nebulier. Could you tell me some comments to setup the model and phycisic model in StarCCM+?

Please provide any information to me.

Thanks.

Joyal

lukasfischer July 14, 2013 11:08

Quote:

Originally Posted by s1300044 (Post 439478)
Hi, Lukas,
I would like to simulate a system like nebulier. Could you tell me some comments to setup the model and phycisic model in StarCCM+?

Please provide any information to me.

Thanks.

Joyal

Sorry, but i only used OpenFOAM to model this. btw, a video of the final simulation is at http://www.youtube.com/watch?v=Ilme9QHcKFc


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