[fsch1]

Hello,
I'm running the reactingPrcelFilmFoam with rain onto a sphere.
But I get a reproducable error:

mpirun --np 4 --hostfile AllMachines reactingParcelFilmFoam -parallel > log_reacting2
[2] #0  Foam::error::printStack(Foam::Ostream&) at ??:?
[2] #1  Foam::sigFpe::sigHandler(int) at ??:?
[2] #2  ? in "/lib/x86_64-linux-gnu/libc.so.6"
[2] #3  log in "/lib/x86_64-linux-gnu/libm.so.6"
[2] #4  Foam::H2O::mu(double, double) const at ??:?
[2] #5  Foam::regionModels::surfaceFilmModels::liquidFilmThermo::mu() const at ??:?
[2] #6  Foam::regionModels::surfaceFilmModels::liquidViscosity::correct(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) at ??:?
[2] #7  Foam::regionModels::surfaceFilmModels::thermoSingleLayer::solveEnergy() at ??:?
[2] #8  Foam::regionModels::surfaceFilmModels::thermoSingleLayer::evolveRegion() at ??:?
[2] #9  Foam::regionModels::regionModel::evolve() at ??:?
[2] #10  ? at ??:?
[2] #11  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
[2] #12  ? at ??:?

Does anybody got an idea what to change? In serial it's running perfectly...


Here the last two timesteps in the log-file:

Courant Number mean: 0.00336716026 max: 0.2003539134
Film max Courant number: 0.118061767
deltaT = 2.073198415e-05
Time = 1.12667651782929123527310366625897586345672607421875


Solving 3-D cloud reactingCloud1

Cloud: reactingCloud1 injector: model1
Added 1 new parcels

Cloud: reactingCloud1
Current number of parcels       = 1336
Current mass in system          = 0.002272204635
Linear momentum                 = (0.01426717076 -9.194179476e-06 -0.008108808571)
|Linear momentum|                = 0.01641051561
Linear kinetic energy           = 0.07473454549
model1:
number of parcels added     = 22423
mass introduced             = 0.02816691295
Parcels absorbed into film      = 708
New film detached parcels       = 359
New film splash parcels         = 932
Parcel fate (number, mass)
- escape                      = 575, 0.0007306381428
- stick                       = 0, 0
Temperature min/max             = 299.65107, 300.0447904
Mass transfer phase change      = 0


Evolving thermoSingleLayer for region wallFilmRegion
diagonal:  Solving for deltaf*rhof, Initial residual = 0, Final residual = 0, No Iterations 0
smoothSolver:  Solving for Ufx, Initial residual = 0.0004691227802, Final residual = 8.353019671e-13, No Iterations 3
smoothSolver:  Solving for Ufy, Initial residual = 0.0004440606956, Final residual = 9.299441697e-13, No Iterations 3
smoothSolver:  Solving for Ufz, Initial residual = 0.0003639826031, Final residual = 6.334200633e-13, No Iterations 3
smoothSolver:  Solving for hf, Initial residual = 0.0007227416005, Final residual = 8.001247771e-13, No Iterations 3
smoothSolver:  Solving for deltaf, Initial residual = 0.002681936015, Final residual = 1.417259447e-11, No Iterations 3
diagonal:  Solving for deltaf*rhof, Initial residual = 0, Final residual = 0, No Iterations 0
smoothSolver:  Solving for deltaf, Initial residual = 2.679234213e-05, Final residual = 2.915438171e-12, No Iterations 2
diagonal:  Solving for deltaf*rhof, Initial residual = 0, Final residual = 0, No Iterations 0
smoothSolver:  Solving for deltaf, Initial residual = 2.660526029e-05, Final residual = 1.048142881e-11, No Iterations 2
diagonal:  Solving for deltaf*rhof, Initial residual = 0, Final residual = 0, No Iterations 0
smoothSolver:  Solving for deltaf, Initial residual = 4.02713407e-07, Final residual = 5.588853589e-14, No Iterations 2
diagonal:  Solving for deltaf*rhof, Initial residual = 0, Final residual = 0, No Iterations 0
smoothSolver:  Solving for Ufx, Initial residual = 0.001348005764, Final residual = 4.203500098e-13, No Iterations 3
smoothSolver:  Solving for Ufy, Initial residual = 0.0004765992185, Final residual = 3.690060877e-13, No Iterations 3
smoothSolver:  Solving for Ufz, Initial residual = 0.0002895922165, Final residual = 3.100166034e-13, No Iterations 3
smoothSolver:  Solving for hf, Initial residual = 0.05942934433, Final residual = 6.175152751e-14, No Iterations 3
smoothSolver:  Solving for deltaf, Initial residual = 7.196314673e-06, Final residual = 1.313969295e-12, No Iterations 2
diagonal:  Solving for deltaf*rhof, Initial residual = 0, Final residual = 0, No Iterations 0
smoothSolver:  Solving for deltaf, Initial residual = 3.230331151e-10, Final residual = 3.150961891e-13, No Iterations 1
diagonal:  Solving for deltaf*rhof, Initial residual = 0, Final residual = 0, No Iterations 0
smoothSolver:  Solving for deltaf, Initial residual = 4.637893188e-13, Final residual = 4.637893188e-13, No Iterations 0
diagonal:  Solving for deltaf*rhof, Initial residual = 0, Final residual = 0, No Iterations 0
smoothSolver:  Solving for deltaf, Initial residual = 4.637779949e-13, Final residual = 4.637779949e-13, No Iterations 0
diagonal:  Solving for deltaf*rhof, Initial residual = 0, Final residual = 0, No Iterations 0

Surface film: thermoSingleLayer
added mass         = 0.0005961575051
current mass       = 0.002751873343
min/max(mag(U))    = 4.552324176e-06, 1.670292839
min/max(delta)     = 3.262762823e-10, 0.002147979312
coverage           = 0.08537182179
injected mass      = 0.002257751105
- patch: region0_to_wallFilmRegion_Spiegel_gesamt: 0.002257751105
min/mean/max(T)    = -118.8767876, 299.9606451, 300.0737498
mass phase change  = 0
vapourisation rate = 0

diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
PIMPLE: iteration 1
DILUPBiCG:  Solving for Ux, Initial residual = 4.774473578e-06, Final residual = 1.613095847e-09, No Iterations 1
DILUPBiCG:  Solving for Uy, Initial residual = 3.756932482e-05, Final residual = 2.595369948e-08, No Iterations 1
DILUPBiCG:  Solving for Uz, Initial residual = 2.989558412e-05, Final residual = 6.391515742e-09, No Iterations 1
DILUPBiCG:  Solving for O2, Initial residual = 0.0004161546431, Final residual = 3.672645244e-08, No Iterations 1
DILUPBiCG:  Solving for H2O, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCG:  Solving for h, Initial residual = 0.003491197827, Final residual = 1.481118942e-09, No Iterations 2
T gas min/max   = 299.9735727, 300.0360395
GAMG:  Solving for p_rgh, Initial residual = 0.001156861795, Final residual = 6.238410992e-08, No Iterations 3
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 7.874500283e-12, global = -3.265306044e-12, cumulative = 4.919632552e-10
GAMG:  Solving for p_rgh, Initial residual = 4.300866817e-07, Final residual = 7.61214236e-10, No Iterations 3
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 9.638893995e-14, global = -1.261879254e-15, cumulative = 4.919619933e-10
GAMG:  Solving for p_rgh, Initial residual = 2.33478321e-09, Final residual = 1.448029858e-11, No Iterations 3
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 1.996194293e-15, global = 2.584716104e-16, cumulative = 4.919622518e-10
GAMG:  Solving for p_rgh, Initial residual = 6.908149188e-11, Final residual = 5.805426344e-13, No Iterations 3
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 2.481609042e-16, global = -4.516502538e-18, cumulative = 4.919622472e-10
PIMPLE: iteration 2
DILUPBiCG:  Solving for Ux, Initial residual = 2.361148179e-06, Final residual = 2.203516898e-09, No Iterations 1
DILUPBiCG:  Solving for Uy, Initial residual = 1.890741739e-05, Final residual = 6.802702732e-09, No Iterations 1
DILUPBiCG:  Solving for Uz, Initial residual = 1.500670702e-05, Final residual = 5.349820067e-09, No Iterations 1
DILUPBiCG:  Solving for O2, Initial residual = 0.0004763365249, Final residual = 4.006099592e-08, No Iterations 1
DILUPBiCG:  Solving for H2O, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCG:  Solving for h, Initial residual = 0.001538415813, Final residual = 7.699140697e-10, No Iterations 2
T gas min/max   = 299.9735703, 300.0332232
GAMG:  Solving for p_rgh, Initial residual = 7.872105478e-05, Final residual = 9.410877388e-09, No Iterations 3
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 1.189758131e-12, global = -8.658627291e-13, cumulative = 4.910963845e-10
GAMG:  Solving for p_rgh, Initial residual = 4.866550957e-08, Final residual = 7.69534777e-11, No Iterations 3
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 9.898296426e-15, global = -4.272769698e-17, cumulative = 4.910963418e-10
GAMG:  Solving for p_rgh, Initial residual = 6.472856085e-10, Final residual = 2.913561058e-12, No Iterations 3
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 5.340546246e-16, global = -9.911746922e-17, cumulative = 4.910962427e-10
GAMG:  Solving for p_rgh, Initial residual = 1.912666082e-11, Final residual = 2.763374625e-13, No Iterations 3
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 1.988935974e-16, global = -4.849715864e-19, cumulative = 4.910962422e-10
DILUPBiCG:  Solving for epsilon, Initial residual = 4.980080585e-07, Final residual = 3.237954098e-10, No Iterations 1
DILUPBiCG:  Solving for k, Initial residual = 3.570996542e-05, Final residual = 1.728293291e-09, No Iterations 1
ExecutionTime = 1213.99 s  ClockTime = 1218 s

Courant Number mean: 0.003355201296 max: 0.1996423283
Film max Courant number: 0.1171148411
deltaT = 2.08059609e-05
Time = 1.126697323790192939441112685017287731170654296875


Solving 3-D cloud reactingCloud1
Cloud: reactingCloud1
Current number of parcels       = 1337
Current mass in system          = 0.002278235033
Linear momentum                 = (0.01426757853 -9.2759301e-06 -0.008108441306)
|Linear momentum|                = 0.0164106887
Linear kinetic energy           = 0.07473325128
model1:
number of parcels added     = 22423
mass introduced             = 0.02816691295
Parcels absorbed into film      = 708
New film detached parcels       = 360
New film splash parcels         = 932
Parcel fate (number, mass)
- escape                      = 575, 0.0007306381428
- stick                       = 0, 0
Temperature min/max             = 299.65107, 300.0447904
Mass transfer phase change      = 0


Evolving thermoSingleLayer for region wallFilmRegion
diagonal:  Solving for deltaf*rhof, Initial residual = 0, Final residual = 0, No Iterations 0
smoothSolver:  Solving for Ufx, Initial residual = 0.0005563004946, Final residual = 8.347499125e-13, No Iterations 3
smoothSolver:  Solving for Ufy, Initial residual = 0.0004842313087, Final residual = 7.369156009e-13, No Iterations 3
smoothSolver:  Solving for Ufz, Initial residual = 0.0003696127878, Final residual = 6.177397222e-13, No Iterations 3
smoothSolver:  Solving for hf, Initial residual = 0.03199865499, Final residual = 4.415082567e-12, No Iterations 3

[fsch1]

Hello,
is there really nobody who has an idea where to look for the error?
It happens always after a parcel detaches..

Or do you need more information?

[fsch1]

Does nobody can help me?

Is the error caused by rho (or mu in the post below)?

[2] #0  Foam::error::printStack(Foam::Ostream&) at ??:?
[2] #1  Foam::sigFpe::sigHandler(int) at ??:?
[2] #2  ? in "/lib/x86_64-linux-gnu/libc.so.6"
[2] #3  ? in "/lib/x86_64-linux-gnu/libm.so.6"
[2] #4  pow in "/lib/x86_64-linux-gnu/libm.so.6"
[2] #5  Foam::H2O::rho(double, double) const at ??:?
[2] #6  Foam::regionModels::surfaceFilmModels::liquidFilmThermo::rho() const at ??:?
[2] #7  Foam::regionModels::surfaceFilmModels::thermoSingleLayer::correctThermoFields() at ??:?
[2] #8  Foam::regionModels::surfaceFilmModels::kinematicSingleLayer::preEvolveRegion() at ??:?
[2] #9  Foam::regionModels::surfaceFilmModels::thermoSingleLayer::preEvolveRegion() at ??:?
[2] #10  Foam::regionModels::regionModel::evolve() at ??:?
[2] #11  ? at ??:?
[2] #12  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
[2] #13  ? at ??:?

[cyw]

Have you got and idea about the reason?