CFDEM Coupling blows up : Help !

Submitted by chaitanya_wadi on Thu, 07/25/2013 - 16:45

Hi I am a regular user of CFDEM softwares OpnFOAM and LIGGGHTS.
I want to simulate channel flow for which I am coupling channel flow tutorial in openfoam with liggghts ( after modifying some parameters)
I am getting following error which I don't understand why.

Create time

Create mesh for time = 0

Reading field p

Reading physical velocity field U
Note: only if voidfraction at boundary is 1, U is superficial velocity!!!

Reading momentum exchange field Ksl

Reading voidfraction field voidfraction = (Vgas/Vparticle)

Creating dummy density field rho

Reading particle velocity field Us

Reading/calculating face flux field phi

Selecting incompressible transport model Newtonian
Selecting turbulence model type RASModel
Selecting RAS turbulence model laminar

Reading g
Selecting locateModel engine
Selecting dataExchangeModel twoWayMPI
Starting up LIGGGHTS for first time execution
Executing input script '../DEM/in.foumeterpaticles'
LIGGGHTS (Version LIGGGHTS-PUBLIC 2.3.2, compiled 2013-04-12-14:14:24 by gauss1 based on LAMMPS 20 Apr 2012)
Created orthogonal box = (-0.003 -0.002 -0.001) to (4.1 2.04 2.077)
2 by 1 by 1 MPI processor grid
nevery as specified in LIGGGHTS is overriden by calling external program (cfd_datacoupling_mpi.cpp:57)
0 atoms in group nve_group
Fix particledistribution/discrete (id pdd1): distribution based on mass%:
pts1: d=5.000000e-02 (max. bounding sphere) mass%=100.000000%
Fix particledistribution/discrete (id pdd1): distribution based on number%:
pts1: d=5.000000e-02 (max. bounding sphere) number%=100.000000%
Setting up run ...
INFO: Resetting random generator for region bc
INFO: Particle insertion ins: 0.000000 particles every 0 steps - particle rate 0.000000 (mass rate 0.000000)
Memory usage per processor = 12.6764 Mbytes
Step Atoms KinEng 1 ts[1] ts[2] Volume
0 0 -0 0 0 0 17.410161
INFO: Particle insertion ins: inserted 296 particle templates (mass 38.746309) at step 1
- a total of 296 particle templates (mass 38.746309) inserted so far.
1 296 7.5361509e-12 0 0 0 17.410161
Loop time of 0.000897527 on 2 procs for 1 steps with 296 atoms

Pair time (%) = 1.43051e-06 (0.159384)
Neigh time (%) = 0.000635505 (70.8062)
Comm time (%) = 1.00136e-05 (1.11569)
Outpt time (%) = 2.34842e-05 (2.61655)
Other time (%) = 0.000227094 (25.3022)

Nlocal: 148 ave 156 max 140 min
Histogram: 1 0 0 0 0 0 0 0 0 1
Nghost: 1 ave 2 max 0 min
Histogram: 1 0 0 0 0 0 0 0 0 1
Neighs: 0.5 ave 1 max 0 min
Histogram: 1 0 0 0 0 0 0 0 0 1

Total # of neighbors = 1
Ave neighs/atom = 0.00337838
Neighbor list builds = 1
Dangerous builds = 0
Selecting IOModel off
Selecting voidFractionModel centre
centreVoidFraction constructor done
Selecting averagingModel dense
Selecting clockModel standardClock
start clock measurement at t >1e-05
Selecting meshMotionModel noMeshMotion

CFDEMcoupling version: cfdem-2.5.3

, compatible to LIGGGHTS version: 2.3.2

If BC are important, please provide volScalarFields -imp/expParticleForces-
ignoring ddt(voidfraction)
Selecting forceModel KochHillDrag
Selecting forceModel Archimedes
accounting for Archimedes on DEM and CFD side!
Selecting momCoupleModel implicitCouple

Selecting liggghtsCommandModel execute ,provide dicts, numbered from 0 to n
liggghtsCommand run 100
firstCouplingStep = 1
lastCouplingStep = 1000000000
couplingStepInterval = 1
cg is set to: 1
solving volume averaged Navier Stokes equations of type B

Reading/calculating field UMean

Reading/calculating field pMean

Reading/calculating field UPrime2Mean

Reading/calculating field pPrime2Mean

fieldAverage: starting averaging at time 0

Starting time loop

Time = 1e-05

Courant Number mean: 1.51391e-05 max: 2.63183e-05
Starting up LIGGGHTS
firstCouplingStep = 1
lastCouplingStep = 1000000000
couplingStepInterval = 1
Executing command: 'run 100 '
Setting up run ...
Memory usage per processor = 13.5168 Mbytes
Step Atoms KinEng 1 ts[1] ts[2] Volume
1 296 7.5361509e-12 0 0 0 17.410161
CFD Coupling established at step 100
101 296 3.0386214e-07 6.8793777e-15 0 0 17.410161
Loop time of 0.00129139 on 2 procs for 100 steps with 296 atoms

Pair time (%) = 7.71284e-05 (5.97249)
Neigh time (%) = 0 (0)
Comm time (%) = 4.50611e-05 (3.48934)
Outpt time (%) = 1.23978e-05 (0.96003)
Other time (%) = 0.00115681 (89.5781)

Nlocal: 148 ave 156 max 140 min
Histogram: 1 0 0 0 0 0 0 0 0 1
Nghost: 1 ave 2 max 0 min
Histogram: 1 0 0 0 0 0 0 0 0 1
Neighs: 0.5 ave 1 max 0 min
Histogram: 1 0 0 0 0 0 0 0 0 1

Total # of neighbors = 1
Ave neighs/atom = 0.00337838
Neighbor list builds = 0
Dangerous builds = 0
LIGGGHTS finished

Coupling...

timeStepFraction() = 1
calculating ddt(voidfraction) based on couplingTime
update Ksl.internalField()
[0] #0 Foam::error::printStack(Foam::Ostream&)[1] #0 Foam::error::printStack(Foam::Ostream&) in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[1] #1 Foam::sigFpe::sigHandler(int) in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[0] #1 Foam::sigFpe::sigHandler(int) in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[1] #2 in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[0] #2 in "/lib/x86_64-linux-gnu/libc.so.6"
[1] #3 Foam::divide(Foam::Field&, Foam::UList const&, Foam::UList const&) in "/lib/x86_64-linux-gnu/libc.so.6"
[0] #3 Foam::divide(Foam::Field&, Foam::UList const&, Foam::UList const&) in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[1] #4 in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[0] #4 void Foam::divide(Foam::GeometricField&, Foam::GeometricField const&, Foam::GeometricField const&)void Foam::divide(Foam::GeometricField&, Foam::GeometricField const&, Foam::GeometricField const&) in "/home/gauss1/OpenFOAM/gauss1-2.1.1/platforms/linux64GccDPOpt/bin/cfdemSolverPiso"
[0] #5 in "/home/gauss1/OpenFOAM/gauss1-2.1.1/platforms/linux64GccDPOpt/bin/cfdemSolverPiso"
[1] #5

[1] at cfdemSolverPiso.C:0
[1] #6 [0] at cfdemSolverPiso.C:0
[0] #6

[0] in "/home/gauss1/OpenFOAM/gauss1-2.1.1/platforms/linux64GccDPOpt/bin/cfdemSolverPiso"
[0] #7 __libc_start_main[1] in "/home/gauss1/OpenFOAM/gauss1-2.1.1/platforms/linux64GccDPOpt/bin/cfdemSolverPiso"
[1] #7 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
[0] #8 in "/lib/x86_64-linux-gnu/libc.so.6"
[1] #8

[0] in "/home/gauss1/OpenFOAM/gauss1-2.1.1/platforms/linux64GccDPOpt/bin/cfdemSolverPiso"
[gauss1:11420] *** Process received signal ***
[gauss1:11420] Signal: Floating point exception (8)
[gauss1:11420] Signal code: (-6)
[gauss1:11420] Failing at address: 0x3e800002c9c
[1] in "/home/gauss1/OpenFOAM/gauss1-2.1.1/platforms/linux64GccDPOpt/bin/cfdemSolverPiso"
[gauss1:11421] *** Process received signal ***
[gauss1:11421] Signal: Floating point exception (8)
[gauss1:11421] Signal code: (-6)
[gauss1:11421] Failing at address: 0x3e800002c9d
[gauss1:11420] [ 0] /lib/x86_64-linux-gnu/libc.so.6(+0x364a0) [0x7f5489c3a4a0]
[gauss1:11420] [ 1] /lib/x86_64-linux-gnu/libc.so.6(gsignal+0x35) [0x7f5489c3a425]
[gauss1:11420] [ 2] /lib/x86_64-linux-gnu/libc.so.6(+0x364a0) [0x7f5489c3a4a0]
[gauss1:11420] [ 3] /opt/openfoam211/platforms/linux64GccDPOpt/lib/libOpenFOAM.so(_ZN4Foam6divideERNS_5FieldIdEERKNS_5UListIdEES6_+0x9f) [0x7f548ad8d6bf]
[gauss1:11420] [ 4] cfdemSolverPiso(_ZN4Foam6divideINS_12fvPatchFieldENS_7volMeshEEEvRNS_14GeometricFieldIdT_T0_EERKS6_S9_+0x8f) [0x4307bf]
[gauss1:11420] [ 5] cfdemSolverPiso() [0x421e00]
[gauss1:11420] [ 6] cfdemSolverPiso() [0x418d9b]
[gauss1:11420] [ 7] /lib/x86_64-linux-gnu/libc.so.6(__libc_start_main+0xed) [0x7f5489c2576d]
[gauss1:11420] [ 8] cfdemSolverPiso() [0x41eafd]
[gauss1:11420] *** End of error message ***
[gauss1:11421] [ 0] /lib/x86_64-linux-gnu/libc.so.6(+0x364a0) [0x7fa2b95944a0]
[gauss1:11421] [ 1] /lib/x86_64-linux-gnu/libc.so.6(gsignal+0x35) [0x7fa2b9594425]
[gauss1:11421] [ 2] /lib/x86_64-linux-gnu/libc.so.6(+0x364a0) [0x7fa2b95944a0]
[gauss1:11421] [ 3] /opt/openfoam211/platforms/linux64GccDPOpt/lib/libOpenFOAM.so(_ZN4Foam6divideERNS_5FieldIdEERKNS_5UListIdEES6_+0x9f) [0x7fa2ba6e76bf]
[gauss1:11421] [ 4] cfdemSolverPiso(_ZN4Foam6divideINS_12fvPatchFieldENS_7volMeshEEEvRNS_14GeometricFieldIdT_T0_EERKS6_S9_+0x8f) [0x4307bf]
[gauss1:11421] [ 5] cfdemSolverPiso() [0x421e00]
[gauss1:11421] [ 6] cfdemSolverPiso() [0x418d9b]
[gauss1:11421] [ 7] /lib/x86_64-linux-gnu/libc.so.6(__libc_start_main+0xed) [0x7fa2b957f76d]
[gauss1:11421] [ 8] cfdemSolverPiso() [0x41eafd]
[gauss1:11421] *** End of error message ***
--------------------------------------------------------------------------
mpirun noticed that process rank 0 with PID 11420 on node gauss1 exited on signal 8 (Floating point exception).
--------------------------------------------------------------------------
/*---------------------------------------------------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.1.1 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
Build : 2.1.1-221db2718bbb
Exec : reconstructPar
Date : Jul 25 2013
Time : 20:06:06
Host : "gauss1"
PID : 11492
Case : /home/gauss1/CFDEM/CFDEMcoupling-PUBLIC/tutorials/cfdemSolverPiso/cfdem/CFDEM/ddp/channelflow/CFD
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster
allowSystemOperations : Disallowing user-supplied system call operations

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

--> FOAM FATAL ERROR:
No times selected

From function reconstructPar
in file reconstructPar.C at line 139.

FOAM exiting