cfdemSolverIB - particle deletion

Submitted by Tonyxz on Fri, 12/16/2016 - 16:13

Dear all,
I modified the twoSpheresGlowonskiMPI Tutorial for the cfdemSolverIB such that the inlet and outlet are not fixed walls but patches with fixed velocity and fixed pressure, respectively. The two particles are inserted with flow velocity and move towards the outlet until the first particle is about to leave the region and be deleted. Then the error message below occurs.
I attached the case and would be very happy about any suggestions on how to fix this.

Thank you,
Tony

-----------------------------------------------------------------------------------------------------------
Time = 0.0507

Selected 1 cells for refinement out of 2480.
Refined from 2480 to 2487 cells.
Selected 0 split points out of a possible 179.
Courant Number mean: 0.00989863 max: 0.316167
- evolve()

timeStepFraction() = 1
Starting up LIGGGHTS
Executing command: 'run 10 '
run 10
Setting up run ...
Memory usage per processor = 40.4865 Mbytes
Step Atoms KinEng rke Volume
5061 2 31.579261 0 2
CFD Coupling established at step 5070
5071 2 35.469635 0 2
Loop time of 5.20003e-05 on 4 procs for 10 steps with 2 atoms

Pair time (%) = 1.25007e-06 (2.40396)
Neigh time (%) = 0 (0)
Comm time (%) = 1.80005e-05 (34.6161)
Outpt time (%) = 5.99993e-06 (11.5383)
Other time (%) = 2.67498e-05 (51.4416)

Nlocal: 0.5 ave 1 max 0 min
Histogram: 2 0 0 0 0 0 0 0 0 2
Nghost: 1.5 ave 2 max 1 min
Histogram: 2 0 0 0 0 0 0 0 0 2
Neighs: 0.5 ave 1 max 0 min
Histogram: 2 0 0 0 0 0 0 0 0 2

Total # of neighbors = 2
Ave neighs/atom = 1
Neighbor list builds = 0
Dangerous builds = 0
LIGGGHTS finished
nr particles = 2
SetInterFace
SetInterFace done
evolve done.
DILUPBiCG: Solving for Ux, Initial residual = 0.0148187, Final residual = 6.7233e-07, No Iterations 2
DILUPBiCG: Solving for Uy, Initial residual = 0.0186391, Final residual = 9.30683e-07, No Iterations 2
DILUPBiCG: Solving for Uz, Initial residual = 0.0508122, Final residual = 1.13935e-07, No Iterations 2
DICPCG: Solving for p, Initial residual = 0.952282, Final residual = 5.26765e-07, No Iterations 37
time step continuity errors : sum local = 7.71607e-10, global = -8.84254e-13, cumulative = 1.35398e-10
DICPCG: Solving for p, Initial residual = 0.481984, Final residual = 8.66238e-07, No Iterations 34
time step continuity errors : sum local = 2.45525e-10, global = 3.99138e-12, cumulative = 1.39389e-10
DICPCG: Solving for p, Initial residual = 0.0986404, Final residual = 9.99391e-07, No Iterations 31
time step continuity errors : sum local = 1.24615e-10, global = 7.84079e-13, cumulative = 1.40173e-10
DICPCG: Solving for p, Initial residual = 0.0110188, Final residual = 9.92681e-07, No Iterations 27
time step continuity errors : sum local = 1.14316e-10, global = 4.39973e-12, cumulative = 1.44573e-10
particleCloud.calcVelocityCorrection()
DICPCG: Solving for phiIB, Initial residual = 0.0928716, Final residual = 8.68297e-07, No Iterations 29
ExecutionTime = 22.95 s ClockTime = 23 s

Time = 0.0508

Selected 0 cells for refinement out of 2487.
Selected 30 split points out of a possible 179.
Unrefined from 2487 to 2277 cells.
Courant Number mean: 0.00973243 max: 0.336885
- evolve()

timeStepFraction() = 1
Starting up LIGGGHTS
Executing command: 'run 10 '
run 10
Setting up run ...
Memory usage per processor = 40.4865 Mbytes
Step Atoms KinEng rke Volume
5071 1 0.14694091 0 2
CFD Coupling established at step 5080
5081 1 0.14823158 0 2
Loop time of 4.52498e-05 on 4 procs for 10 steps with 1 atoms

Pair time (%) = 9.99891e-07 (2.20971)
Neigh time (%) = 0 (0)
Comm time (%) = 1.14999e-05 (25.4141)
Outpt time (%) = 5.99993e-06 (13.2596)
Other time (%) = 2.67501e-05 (59.1166)

Nlocal: 0.25 ave 1 max 0 min
Histogram: 3 0 0 0 0 0 0 0 0 1
Nghost: 0.75 ave 1 max 0 min
Histogram: 1 0 0 0 0 0 0 0 0 3
Neighs: 0 ave 0 max 0 min
Histogram: 4 0 0 0 0 0 0 0 0 0

Total # of neighbors = 0
Ave neighs/atom = 0
Neighbor list builds = 0
Dangerous builds = 0
LIGGGHTS finished
nr particles = 2
SetInterFace
[2] #0 Foam::error::printStack(Foam::Ostream&)[0] #0 Foam::error::printStack(Foam::Ostream&)[1] #0 Foam::error::printStack(Foam::Ostream&)[3] #0 Foam::error::printStack(Foam::Ostream&) at ??:?
[0] #1 Foam::sigFpe::sigHandler(int) at ??:?
at ??:?
[2] #1 Foam::sigFpe::sigHandler(int)[1] #1 Foam::sigFpe::sigHandler(int) at ??:?
[3] #1 Foam::sigFpe::sigHandler(int) at ??:?
[0] #2 ? at ??:?
[1] #2 ? at ??:?
[2] #2 ? at ??:?
[3] #2 ? in "/usr/lib/libc.so.6"
[0] #3 Foam::voidFractionModel::pointInParticle(int, Foam::Vector, Foam::Vector, double) const in "/usr/lib/libc.so.6"
[2] #3 Foam::voidFractionModel::pointInParticle(int, Foam::Vector, Foam::Vector, double) const in "/usr/lib/libc.so.6"
[1] #3 Foam::voidFractionModel::pointInParticle(int, Foam::Vector, Foam::Vector, double) const in "/usr/lib/libc.so.6"
[3] #3 Foam::voidFractionModel::pointInParticle(int, Foam::Vector, Foam::Vector, double) const at voidFractionModel.C:?
[0] #4 Foam::cfdemCloudIB::setInterFace(Foam::GeometricField&) at voidFractionModel.C:?
[2] #4 Foam::cfdemCloudIB::setInterFace(Foam::GeometricField&) at voidFractionModel.C:?
[1] #4 Foam::cfdemCloudIB::setInterFace(Foam::GeometricField&) at voidFractionModel.C:?
[3] #4 Foam::cfdemCloudIB::setInterFace(Foam::GeometricField&) at ??:?
[0] #5 Foam::cfdemCloudIB::evolve(Foam::GeometricField&, Foam::GeometricField&) at ??:?
at ??:?
[2] #5 at ??:?
Foam::cfdemCloudIB::evolve(Foam::GeometricField&, Foam::GeometricField&)[1] #5 Foam::cfdemCloudIB::evolve(Foam::GeometricField&, Foam::GeometricField&)[3] #5 Foam::cfdemCloudIB::evolve(Foam::GeometricField&, Foam::GeometricField&) at ??:?
[3] #6 at ??:?
[0] #6 at ??:?
[1] #6 at ??:?
[2] #6 ???? at ??:?
[3] #7 __libc_start_main at ??:?
[1] #7 __libc_start_main at ??:?
[0] #7 __libc_start_main at ??:?
[2] #7 __libc_start_main in "/usr/lib/libc.so.6"
[1] #8 in "/usr/lib/libc.so.6"
[3] #8 in "/usr/lib/libc.so.6"
[0] #8 ?? in "/usr/lib/libc.so.6"
[2] #8 ?? at ??:?
[boreas:17748] *** Process received signal ***
[boreas:17748] Signal: Floating point exception (8)
[boreas:17748] Signal code: (-6)
[boreas:17748] Failing at address: 0x3ea00004554
[boreas:17748] [ 0] /usr/lib/libc.so.6(+0x330b0)[0x7f44766f60b0]
[boreas:17748] [ 1] /usr/lib/libc.so.6(gsignal+0xcf)[0x7f44766f604f]
[boreas:17748] [ 2] /usr/lib/libc.so.6(+0x330b0)[0x7f44766f60b0]
[boreas:17748] [ 3] /home/laura/OpenFOAM/laura-3.0.0/platforms/linux64GccDPInt32Opt/lib/liblagrangianCFDEM-PUBLIC-3.0.0.so(_ZNK4Foam17voidFractionModel15pointInParticleEiNS_6VectorIdEES2_d+0x56)[0x7f4478a97716]
[boreas:17748] [ 4] /home/laura/OpenFOAM/laura-3.0.0/platforms/linux64GccDPInt32Opt/lib/liblagrangianCFDEM-PUBLIC-3.0.0.so(_ZN4Foam12cfdemCloudIB12setInterFaceERNS_14GeometricFieldIdNS_12fvPatchFieldENS_7volMeshEEE+0x3a2)[0x7f4478a043a2]
[boreas:17748] [ 5] at ??:?
[boreas:17750] *** Process received signal ***
[boreas:17750] Signal: Floating point exception (8)
[boreas:17750] Signal code: (-6)
[boreas:17750] Failing at address: 0x3ea00004556
[boreas:17750] [ 0] /usr/lib/libc.so.6(+0x330b0)[0x7f1accbb00b0]
[boreas:17750] [ 1] /usr/lib/libc.so.6(gsignal+0xcf)[0x7f1accbb004f]
[boreas:17750] [ 2] /usr/lib/libc.so.6(+0x330b0)[0x7f1accbb00b0]
[boreas:17750] [ 3] /home/laura/OpenFOAM/laura-3.0.0/platforms/linux64GccDPInt32Opt/lib/liblagrangianCFDEM-PUBLIC-3.0.0.so(_ZN4Foam12cfdemCloudIB6evolveERNS_14GeometricFieldIdNS_12fvPatchFieldENS_7volMeshEEES5_+0x26e)[0x7f4478a04ace]
[boreas:17748] [ 6] cfdemSolverIB[0x41b9cf]
[boreas:17748] [ 7] /usr/lib/libc.so.6(__libc_start_main+0xf1)[0x7f44766e3291]
[boreas:17748] [ 8] cfdemSolverIB[0x41eb9a]
[boreas:17748] *** End of error message ***
/home/laura/OpenFOAM/laura-3.0.0/platforms/linux64GccDPInt32Opt/lib/liblagrangianCFDEM-PUBLIC-3.0.0.so(_ZNK4Foam17voidFractionModel15pointInParticleEiNS_6VectorIdEES2_d+0x56)[0x7f1acef51716]
[boreas:17750] [ 4] /home/laura/OpenFOAM/laura-3.0.0/platforms/linux64GccDPInt32Opt/lib/liblagrangianCFDEM-PUBLIC-3.0.0.so(_ZN4Foam12cfdemCloudIB12setInterFaceERNS_14GeometricFieldIdNS_12fvPatchFieldENS_7volMeshEEE+0x3a2)[0x7f1aceebe3a2]
[boreas:17750] [ 5] /home/laura/OpenFOAM/laura-3.0.0/platforms/linux64GccDPInt32Opt/lib/liblagrangianCFDEM-PUBLIC-3.0.0.so(_ZN4Foam12cfdemCloudIB6evolveERNS_14GeometricFieldIdNS_12fvPatchFieldENS_7volMeshEEES5_+0x26e)[0x7f1aceebeace]
[boreas:17750] [ 6] cfdemSolverIB[0x41b9cf]
[boreas:17750] [ 7] /usr/lib/libc.so.6(__libc_start_main+0xf1)[0x7f1accb9d291]
[boreas:17750] [ 8] cfdemSolverIB[0x41eb9a]
[boreas:17750] *** End of error message ***
at ??:?
[boreas:17749] *** Process received signal ***
[boreas:17749] Signal: Floating point exception (8)
[boreas:17749] Signal code: (-6)
[boreas:17749] Failing at address: 0x3ea00004555
[boreas:17749] [ 0] /usr/lib/libc.so.6(+0x330b0)[0x7fe7841390b0]
[boreas:17749] [ 1] /usr/lib/libc.so.6(gsignal+0xcf)[0x7fe78413904f]
[boreas:17749] [ 2] /usr/lib/libc.so.6(+0x330b0)[0x7fe7841390b0]
[boreas:17749] [ 3] at ??:?
[boreas:17747] *** Process received signal ***
[boreas:17747] Signal: Floating point exception (8)
[boreas:17747] Signal code: (-6)
[boreas:17747] Failing at address: 0x3ea00004553
[boreas:17747] [ 0] /usr/lib/libc.so.6(+0x330b0)[0x7fabc25460b0]
[boreas:17747] [ 1] /usr/lib/libc.so.6(gsignal+0xcf)[0x7fabc254604f]
[boreas:17747] [ 2] /usr/lib/libc.so.6(+0x330b0)[0x7fabc25460b0]
[boreas:17747] [ 3] /home/laura/OpenFOAM/laura-3.0.0/platforms/linux64GccDPInt32Opt/lib/liblagrangianCFDEM-PUBLIC-3.0.0.so(_ZNK4Foam17voidFractionModel15pointInParticleEiNS_6VectorIdEES2_d+0x56)[0x7fe7864da716]
[boreas:17749] [ 4] /home/laura/OpenFOAM/laura-3.0.0/platforms/linux64GccDPInt32Opt/lib/liblagrangianCFDEM-PUBLIC-3.0.0.so(_ZNK4Foam17voidFractionModel15pointInParticleEiNS_6VectorIdEES2_d+0x56)[0x7fabc48e7716]
[boreas:17747] [ 4] /home/laura/OpenFOAM/laura-3.0.0/platforms/linux64GccDPInt32Opt/lib/liblagrangianCFDEM-PUBLIC-3.0.0.so(_ZN4Foam12cfdemCloudIB12setInterFaceERNS_14GeometricFieldIdNS_12fvPatchFieldENS_7volMeshEEE+0x3a2)[0x7fe7864473a2]
[boreas:17749] [ 5] /home/laura/OpenFOAM/laura-3.0.0/platforms/linux64GccDPInt32Opt/lib/liblagrangianCFDEM-PUBLIC-3.0.0.so(_ZN4Foam12cfdemCloudIB12setInterFaceERNS_14GeometricFieldIdNS_12fvPatchFieldENS_7volMeshEEE+0x3a2)[0x7fabc48543a2]
[boreas:17747] [ 5] /home/laura/OpenFOAM/laura-3.0.0/platforms/linux64GccDPInt32Opt/lib/liblagrangianCFDEM-PUBLIC-3.0.0.so(_ZN4Foam12cfdemCloudIB6evolveERNS_14GeometricFieldIdNS_12fvPatchFieldENS_7volMeshEEES5_+0x26e)[0x7fe786447ace]
[boreas:17749] [ 6] cfdemSolverIB[0x41b9cf]
[boreas:17749] [ 7] /usr/lib/libc.so.6(__libc_start_main+0xf1)[0x7fe784126291]
[boreas:17749] [ 8] cfdemSolverIB[0x41eb9a]
[boreas:17749] *** End of error message ***
/home/laura/OpenFOAM/laura-3.0.0/platforms/linux64GccDPInt32Opt/lib/liblagrangianCFDEM-PUBLIC-3.0.0.so(_ZN4Foam12cfdemCloudIB6evolveERNS_14GeometricFieldIdNS_12fvPatchFieldENS_7volMeshEEES5_+0x26e)[0x7fabc4854ace]
[boreas:17747] [ 6] cfdemSolverIB[0x41b9cf]
[boreas:17747] [ 7] /usr/lib/libc.so.6(__libc_start_main+0xf1)[0x7fabc2533291]
[boreas:17747] [ 8] cfdemSolverIB[0x41eb9a]
[boreas:17747] *** End of error message ***
--------------------------------------------------------------------------
mpirun noticed that process rank 1 with PID 17748 on node boreas exited on signal 8 (Floating point exception).

AttachmentSize
Binary Data twospheresglowinskimpi.tar_.gz82.75 KB

JohLoh | Tue, 12/20/2016 - 15:48

hello Tony,
I get the same error message. If I use periodic boundaries in DEM, the error somehow disappears. Looks like the particle number "nr particles = 2" isn't updated ... Has anyone experience with this error?
Best
Johannes

cgoniva's picture

cgoniva | Fri, 12/30/2016 - 13:05

Hello,

the IB solver could not handle the particle leaving the domain (switching to periodic domain in DEM works (mostly) as the particle would not be deleted).

Thank you for bringing up the issue and providing the testcase!

In the next release there will be a change in the code to overcome this problem and allow for running your testcase.

Best regards,

Christoph

Tonyxz | Wed, 02/01/2017 - 15:48

Hello,
thank you for the bugfix in the new release!
I installed the new LIGGGHTS and CFDEM version and use it with OpenFOAM 3.0.1.
Unfortunately, I still get the same error when running the case I uploaded with my first post.
Do I have to make other changes to overcome the problem?

Best regards
Tony