[kenyi89]

Dear all,

Without altering anything I'm trying to run the 'column' tutorial for the MPPICFoam solver in OpenFOAM 2.3.1. I run blockMesh followed by MPPICFoam in the case directory and get:

#0 Foam::error:rintStack(Foam::Ostream&) at ~/dpm-mpc/src/OSspecific/POSIX/printStack.C:221
#1 Foam::sigFpe::sigHandler(int) at ~/dpm-mpc/src/OSspecific/POSIX/signals/sigFpe.C:108
#2 in "/lib/x86_64-linux-gnu/libc.so.6"
#3
at ~/dpm-mpc/src/lagrangian/intermediate/lnInclude/MPPICCloud.C:292 (discriminator 1)
#4
at ~/dpm-mpc/src/lagrangian/intermediate/lnInclude/KinematicCloud.C:120
#5
at ~/dpm-mpc/src/lagrangian/intermediate/lnInclude/MPPICCloud.C:179 (discriminator 1)
#6
at ~/dpm-mpc/applications/solvers/lagrangian/DPMFoam/MPPICFoam/../DPMFoam.C:85
#7 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#8
at ??:?

I haven't bothered finding the problem yet because I was assuming the tutorials to run successfully. Can anyone either confirm the crash (in which case I will dive into it) or have a look at what I could have done wrong?

Thanks.

[thiagopl]

Hi,

I did the same here and it started with no errors message.
OF 2.3.1.

PIMPLE: Operating solver in PISO mode

Starting time loop

Courant Number mean: 0 max: 0
Time = 0.0002

Evolving kinematicCloud

Solving 3-D cloud kinematicCloud

--> Cloud: kinematicCloud injector: model1
    Added 2873 new parcels

smoothSolver:  Solving for kinematicCloud:alpha, Initial residual = 2.1334886e-05, Final residual = 1.2261187e-14, No Iterations 1
Cloud: kinematicCloud
    Current number of parcels       = 2873
    Current mass in system          = 0.076947165
    Linear momentum                 = (0 0 -0.00015113283)
   |Linear momentum|                = 0.00015113283
    Linear kinetic energy           = 1.4846822e-07
    model1:
        number of parcels added     = 2873
        mass introduced             = 0.076947165
    Parcels absorbed into film      = 0
    New film detached parcels       = 0
    Parcel fate (number, mass)      : patch sides
      - escape                      = 0, 0
      - stick                       = 0, 0
    Parcel fate (number, mass)      : patch top
      - escape                      = 0, 0
      - stick                       = 0, 0
    Parcel fate (number, mass)      : patch bottom
      - escape                      = 0, 0
      - stick                       = 0, 0
    Min cell volume fraction        = 0.24636485
    Max cell volume fraction        = 0.36328377
    Min dense number of parcels     = 87

smoothSolver:  Solving for U.airx, Initial residual = 0, Final residual = 0, No Iterations 0
smoothSolver:  Solving for U.airy, Initial residual = 0, Final residual = 0, No Iterations 0
smoothSolver:  Solving for U.airz, Initial residual = 1, Final residual = 6.8420209e-10, No Iterations 1
GAMG:  Solving for p, Initial residual = 1, Final residual = 0.0057238502, No Iterations 8
time step continuity errors : sum local = 1.7069106e-08, global = -2.908067e-22, cumulative = -2.908067e-22
GAMG:  Solving for p, Initial residual = 0.00036321186, Final residual = 5.4531259e-07, No Iterations 8
time step continuity errors : sum local = 2.5696931e-11, global = -2.9074366e-22, cumulative = -5.8155037e-22
ExecutionTime = 0.08 s  ClockTime = 0 s

Courant Number mean: 1.5346575e-10 max: 2.2167832e-10
Time = 0.0004...

[kenyi89]

The program starts executing fine, it crashes after 0.041 seconds of simulation (a wallclock time of about 40 seconds on my PC), see also the log file I attached in the previous post.

After I reported I was told a bugfix causing the crash is already implemented in 2.3.x, so the issue is solved.