|
[Sponsors] |
April 11, 2012, 07:15 |
sprayFoam crashes
|
#1 |
New Member
Lukas Fischer
Join Date: Mar 2009
Location: Innsbruck, Austria
Posts: 15
Rep Power: 17 |
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). |
|
May 10, 2012, 09:34 |
|
#2 | |
New Member
Oliver Pasqual
Join Date: May 2011
Posts: 13
Rep Power: 15 |
Quote:
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:rintStack(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:(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 |
||
July 13, 2013, 00:14 |
|
#3 |
New Member
Joyal
Join Date: Jun 2013
Posts: 7
Rep Power: 13 |
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 |
|
July 14, 2013, 12:08 |
|
#4 | |
New Member
Lukas Fischer
Join Date: Mar 2009
Location: Innsbruck, Austria
Posts: 15
Rep Power: 17 |
Quote:
|
||
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
reactingFoam crashes mysteriously | jose_rodrig | OpenFOAM Running, Solving & CFD | 9 | August 4, 2015 11:18 |
modified code crashes in parallel | mabinty | OpenFOAM Programming & Development | 2 | January 11, 2013 06:15 |
checkMesh and decomposePar crashes with salome IDEAS mesh | pajofego | OpenFOAM | 0 | June 23, 2011 09:43 |
[OpenFOAM] paraFoam crashes with channelFoam example | mgdenno | ParaView | 2 | February 5, 2011 08:46 |
TwoPhaseEulerFoam bed tutorial case crashes in 1.6. Stable in 1.5 | hemph | OpenFOAM Bugs | 2 | February 17, 2010 09:10 |