Skip to content

Bug liquidFilmStepWithSprinklers

Summary

tutorials/lagrangian/reactingParcelFoam/liquidFilmStepWithSprinklers breaks after about 1.9s . endTime in controlDict is 1, but particles are injected at 5 and it is suggested to increase endTime to 10. Same error in serial and parallel. Tested on 2 different machines.

Steps to reproduce

Increase endTime from 1 to 5s;./Allrun; At 1.9375s it will stop. Error below.

Relevant logs and/or images

Solving3-D cloud reactingCloud1
Cloud: reactingCloud1
    Current number of parcels       = 0
    Current mass in system          = 0
    Linear momentum                 = (0 0 0)
   |Linear momentum|                = 0
    Linear kinetic energy           = 0
    Average particle per parcel     = 0
    Injector model1:
      - parcels added               = 0
      - mass introduced             = 0
    Surface film:
      - parcels absorbed            = 0
      - mass absorbed               = 0
      - parcels ejected             = 0
      - new splash parcels          = 0
    Parcel fate: system (number, mass)
      - escape                      = 0, 0
    Parcel fate: patch (number, mass) base
      - escape                      = 0, 0
      - stick                       = 0, 0
    Parcel fate: patch (number, mass) outlet
      - escape                      = 0, 0
      - stick                       = 0, 0
    Parcel fate: patch (number, mass) inlet
      - escape                      = 0, 0
      - stick                       = 0, 0
    Parcel fate: patch (number, mass) sides
      - escape                      = 0, 0
      - stick                       = 0, 0
    Parcel fate: patch (number, mass) top
      - escape                      = 0, 0
      - stick                       = 0, 0
    Temperature min/max             = 0, 0
    Mass transfer phase change      = 0
    Mass transfer devolatilisation  = 0
    Mass transfer surface reaction  = 0

diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
PIMPLE: iteration 1

Evolving kinematicThinFilm for region region0
DILUPBiCGStab:  Solving for Uf_filmx, Initial residual = 0.0118513, Final residual = 9.60929e-15, No Iterations 1
DILUPBiCGStab:  Solving for Uf_filmy, Initial residual = 0.0109228, Final residual = 6.0994e-11, No Iterations 1
DILUPBiCGStab:  Solving for Uf_filmz, Initial residual = 0.000543556, Final residual = 8.40472e-10, No Iterations 1
DILUPBiCGStab:  Solving for hf_film, Initial residual = 0.00211817, Final residual = 4.50239e-12, No Iterations 1
DILUPBiCGStab:  Solving for Uf_filmx, Initial residual = 0.0122752, Final residual = 4.82225e-15, No Iterations 1
DILUPBiCGStab:  Solving for Uf_filmy, Initial residual = 0.00786743, Final residual = 4.88628e-11, No Iterations 1
DILUPBiCGStab:  Solving for Uf_filmz, Initial residual = 0.000531824, Final residual = 8.05986e-10, No Iterations 1
DILUPBiCGStab:  Solving for hf_film, Initial residual = 0.000296377, Final residual = 5.39824e-13, No Iterations 1
DILUPBiCGStab:  Solving for Uf_filmx, Initial residual = 0.00201142, Final residual = 4.85637e-15, No Iterations 1
DILUPBiCGStab:  Solving for Uf_filmy, Initial residual = 0.000200432, Final residual = 3.70026e-13, No Iterations 1
DILUPBiCGStab:  Solving for Uf_filmz, Initial residual = 0.000302661, Final residual = 3.90512e-10, No Iterations 1
DILUPBiCGStab:  Solving for hf_film, Initial residual = 5.02967e-05, Final residual = 1.05152e-13, No Iterations 1
DILUPBiCGStab:  Solving for Uf_filmx, Initial residual = 0.000488423, Final residual = 2.39743e-15, No Iterations 1
DILUPBiCGStab:  Solving for Uf_filmy, Initial residual = 9.81888e-05, Final residual = 2.3487e-13, No Iterations 1
DILUPBiCGStab:  Solving for Uf_filmz, Initial residual = 0.000156307, Final residual = 1.99035e-10, No Iterations 1
DILUPBiCGStab:  Solving for hf_film, Initial residual = 1.40484e-05, Final residual = 6.25537e-14, No Iterations 1
Film h min/max   = (4.55048e-06 0.00165964)
Film mag(U) min/max   = (0 2.51884)
smoothSolver:  Solving for Ux, Initial residual = 0.00197335, Final residual = 2.55965e-07, No Iterations 5
smoothSolver:  Solving for Uy, Initial residual = 0.00928135, Final residual = 8.40566e-07, No Iterations 5
smoothSolver:  Solving for Uz, Initial residual = 0.00411528, Final residual = 5.61892e-07, No Iterations 5
DILUPBiCGStab:  Solving for h, Initial residual = 0.00497396, Final residual = 1.94128e-07, No Iterations 1
T gas min/max   = 299.53, 300.008
DICPBiCGStab:  Solving for p_rgh, Initial residual = 0.0032935, Final residual = 0.000297896, No Iterations 11
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 6.43459e-07, global = -3.59919e-08, cumulative = -0.000505954
DICPBiCGStab:  Solving for p_rgh, Initial residual = 0.00119447, Final residual = 0.000117951, No Iterations 10
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 2.54758e-07, global = -3.68368e-08, cumulative = -0.000505991
DICPBiCGStab:  Solving for p_rgh, Initial residual = 0.000190021, Final residual = 8.92924e-07, No Iterations 69
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 1.92866e-09, global = 3.00806e-10, cumulative = -0.000505991
DILUPBiCGStab:  Solving for epsilon, Initial residual = 0.00229694, Final residual = 4.45827e-07, No Iterations 1
DILUPBiCGStab:  Solving for k, Initial residual = 0.00441566, Final residual = 7.895e-07, No Iterations 1
ExecutionTime = 57.35 s  ClockTime = 57 s
Courant Number mean: 0.466217 max: 0.868391
Max film Courant Number: 2.49983 Film velocity magnitude: 4.99966
deltaT = 0.0125
Time = 1.9375


Solving3-D cloud reactingCloud1
[stack trace]
=============
#1  Foam::sigSegv::sigHandler(int) in ~/OpenFOAM/OpenFOAM-v2306/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so
#2  ? in /lib/x86_64-linux-gnu/libc.so.6
#3  Foam::ThermoSurfaceFilm<Foam::KinematicCloud<Foam::Cloud<Foam::ReactingMultiphaseParcel<Foam::ReactingParcel<Foam::ThermoParcel<Foam::KinematicParcel<Foam::particle> > > > > > >::setParcelProperties(Foam::ReactingMultiphaseParcel<Foam::ReactingParcel<Foam::ThermoParcel<Foam::KinematicParcel<Foam::particle> > > >&, int) const in ~/OpenFOAM/ceschia-v2306/platforms/linux64GccDPInt32Opt/lib/liblagrangianIntermediate.so
#4  ? in ~/OpenFOAM/OpenFOAM-v2306/platforms/linux64GccDPInt32Opt/bin/reactingParcelFoam
#5  ? in ~/OpenFOAM/OpenFOAM-v2306/platforms/linux64GccDPInt32Opt/bin/reactingParcelFoam
#6  ? in ~/OpenFOAM/OpenFOAM-v2306/platforms/linux64GccDPInt32Opt/bin/reactingParcelFoam
#7  ? in ~/OpenFOAM/OpenFOAM-v2306/platforms/linux64GccDPInt32Opt/bin/reactingParcelFoam
#8  ? in ~/OpenFOAM/OpenFOAM-v2306/platforms/linux64GccDPInt32Opt/bin/reactingParcelFoam
#9  ? in ~/OpenFOAM/OpenFOAM-v2306/platforms/linux64GccDPInt32Opt/bin/reactingParcelFoam
#10  ? in /lib/x86_64-linux-gnu/libc.so.6
#11  __libc_start_main in /lib/x86_64-linux-gnu/libc.so.6
#12  ? in ~/OpenFOAM/OpenFOAM-v2306/platforms/linux64GccDPInt32Opt/bin/reactingParcelFoam
=============
Segmentation fault (core dumped)

Environment information

  • OpenFOAM version : 2306
  • Operating system : ubuntu 22 and 20
Edited by Alberto ceschin