How to move particles?

Submitted by hulli on Thu, 08/28/2014 - 18:57

Hi all,
In order to learn CDFDEM, I wanted to adapt the cylinder tutorial for my needs. Basically, I want to have an inflow boundary condition at the top face and an outflow boundary condition at the bottom face. Since I am a bit more familiar with OpenFoam I set up the case in OpenFoam first and let it run. In OpenFoam the flow evolves nicely from top to bottom the when using the pisoFoam solver. However, when it comes to CFDEM, I do not understand why there is a “dummy” in the “0” folder? At least in OpenFoam I need U, p…, to get the simulation running. So, I copied all files (U,p…) in the “0” folder of CFDEM. Anyways, I am not able to change the motion of the particles calculated by the DEM routine. Nothing changes, the particles always fall at the bottom of the cylinder, even if I choose extremely high velocity values.

I appreciate your help…

Forums:

j-kerbl's picture

j-kerbl | Mon, 09/01/2014 - 14:04

Hi,
the dummy doesn't really do anything, except protecting an otherwise possibly empty folder from deletion during e.g. compression.
I'm not sure I really get what you're trying to accomplish, but as a start with CFDEMcoupling for OpenFOAMers it is always good to get to know LIGGGHTS, the DEM side of the coupling. Maybe the LIGGGHTS tutorials and/or help pages may introduce it to you a little bit more than the in-files in the DEM folders of the CFDEMcoupling cases.

I hope that helps you out a little bit!

Cheers
Josef

hulli | Mon, 09/01/2014 - 20:17

Hi Josef,
thank you very much for your answer... I stated getting familiar with liggghts.

However, I mostly get the ERROR: Invalid pair style.

I found out that some package is missing. Someone suggested the granular package.
I am a bit confused, because I installed CFDEM as suggested on the web page.

I hope anybody can help me out!!!

Thank you very much!!!

j-kerbl's picture

j-kerbl | Tue, 09/02/2014 - 09:48

Hi,

to check your liggghts installation you can try to cd into ~/LIGGGHTS/LIGGGHTS-Public/examples/LIGGGHTS/meshGran and run 'cfdemLiggghts in.meshGran' there.

What is cfdemSysTest reporting?

Cheers
Josef

hulli | Tue, 09/02/2014 - 20:00

Hi Josef,
Thank you so much for your suggestion to use 'cfdemLiggghts in.meshGran' command. It works.
In order to understand what I am doing, I intended to start the simulation manually as suggested here: http://www.cfdem.com/forums/setting-cfd-dem-coupling.
However, when I started liggghts by the command: liggghts < in.liggghts_init, I always got the error message. But when I tried your suggestion, the tutorial worked. I am somehow still struggeling to couple it manually to a fow which frives my particles.
But thanky you sooooooooooo much!!!

j-kerbl's picture

j-kerbl | Wed, 09/03/2014 - 13:38

Hi again,

I did some short research and I think I found out what the problem might be!
Sorry I overlooked that before, but is the tutorial you started with the fillCylinder from the cfdemPostProc folder?
If so, this is not a coupled case itself, more a post-processing tool.
To bring some light into it, look at the CFD/constant/couplingProperties file.
There the particles have for example 'noDrag' selected as a forceModel.

A better starting point would be cfdemSolverPiso tutorials. Especially the ErgunTestMPI and settlingTestMPI.

I hope thatt does the trick! :-)
Cheers
Josef

hulli | Wed, 09/03/2014 - 20:23

Hi Josef,
thank you so much for your help!!!!!
I did exactly what you suggested. Basically I took the settlingTestMPI couplingProperties file and tried to understand why the filCylinder was not running or coupling. In order to fully understand what it is doing it tried to couple it manually, which worked out thanks to your advice using the cfdemLiggghts in.*** command. Then I could see the particles using that pyton script to rewrite the dump files and to get .vtk files.
In the CFD folder I am a bit struggling to get .vtk files. Basically, I use the command foamToVTK, after the calculation finished. However, even if I have lots of DEM.vtk files I only get a single 0.vtk file from the CFD part, even if I took an incredible small time step or write interval. Do I have to change the write interval somewhere else? Does the write interval from the openFoam part has to be similar than the on for the DEM dump files.
Thanks so much for your help again!!!!

j-kerbl's picture

j-kerbl | Thu, 09/04/2014 - 11:29

Hi,
to have a look at the OpenFOAM data you can do this as usual. Open the controldict.foam in the CFD/system folder with paraview. Remember that CFDEMcoupling runs in parallel, therefore you have to select decomposed in paraview or reconstruct the case before. The CFD writeinterval is set in system/controlDict and DEM dump-interval in the in-file and they do not have to be the same.

If you are interested in more details, there are courses held several times a year.

Cheers
Josef

hulli | Thu, 09/04/2014 - 19:15

Hi Josef,
Thank you so much for your help! Yes that was exactly what I was looking for. Well, basically a course would be very helpful but I have to get into it myself, I have to give a lecture on state of the art numerical simulation techniques in two weeks.... Aaah time is running haha ...
Anyways, if you are still keen to help me out …after some steps calculated … I always end up with that error in the Ergun tutorial… Is this related to the times steps chosen in the DEM / CFD model ?

[3] [1] [2] [0] #0 Foam::error::printStack(Foam::Ostream&)##0 Foam::error::printStack(Foam::Ostream&)0 Foam::error::printStack(Foam::Ostream&)#0 Foam::error::printStack(Foam::Ostream&) in "/opt/openfoam230/platforms/linuxGccDPOpt/lib/libOpenFOAM.so"
[1] #1 Foam::sigFpe::sigHandler(int) in "/opt/openfoam230/platforms/linuxGccDPOpt/lib/libOpenFOAM.so"
[3] #1 Foam::sigFpe::sigHandler(int) in "/opt/openfoam230/platforms/linuxGccDPOpt/lib/libOpenFOAM.so"
[0] #1 Foam::sigFpe::sigHandler(int) in "/opt/openfoam230/platforms/linuxGccDPOpt/lib/libOpenFOAM.so"
[2] #1 Foam::sigFpe::sigHandler(int) in "/opt/openfoam230/platforms/linuxGccDPOpt/lib/libOpenFOAM.so"
[0] #2 Uninterpreted:
[0] #3 in "/opt/openfoam230/platforms/linuxGccDPOpt/lib/libOpenFOAM.so"
[1] #2 Uninterpreted:
[1] #3 in "/opt/openfoam230/platforms/linuxGccDPOpt/lib/libOpenFOAM.so"
[3] #2 Uninterpreted:
[3] #3 in "/usr/lib/libmpi.so.0"
[0] #4 in "/opt/openfoam230/platforms/linuxGccDPOpt/lib/libOpenFOAM.so"
[2] #2 Uninterpreted:
[2] #3 in "/usr/lib/libmpi.so.0"
[3] #4 in "/usr/lib/libmpi.so.0"
[1] #4 in "/usr/lib/openmpi/lib/openmpi/mca_coll_tuned.so"
[0] #5 in "/usr/lib/openmpi/lib/openmpi/mca_coll_tuned.so"
[3] #5 in "/usr/lib/libmpi.so.0"
[2] #4 in "/usr/lib/openmpi/lib/openmpi/mca_coll_tuned.so"
[1] #5 in "/usr/lib/openmpi/lib/openmpi/mca_coll_tuned.so"
[0] #6 PMPI_Allreduce in "/usr/lib/openmpi/lib/openmpi/mca_coll_tuned.so"
[3] #6 PMPI_Allreduce in "/usr/lib/openmpi/lib/openmpi/mca_coll_tuned.so"
[1] #6 PMPI_Allreduce in "/usr/lib/openmpi/lib/openmpi/mca_coll_tuned.so"
[2] #5 in "/usr/lib/libmpi.so.0"
[0] #7 Foam::reduce(double&, Foam::sumOp const&, int, int) in "/usr/lib/libmpi.so.0"
[3] #7 Foam::reduce(double&, Foam::sumOp const&, int, int) in "/usr/lib/libmpi.so.0"
[1] #7 Foam::reduce(double&, Foam::sumOp const&, int, int) in "/opt/openfoam230/platforms/linuxGccDPOpt/lib/openmpi-system/libPstream.so"
[0] #8 Foam::PCG::solve(Foam::Field&, Foam::Field const&, unsigned char) const in "/usr/lib/openmpi/lib/openmpi/mca_coll_tuned.so"
[2] #6 PMPI_Allreduce in "/opt/openfoam230/platforms/linuxGccDPOpt/lib/openmpi-system/libPstream.so"
[3] #8 Foam::PCG::solve(Foam::Field&, Foam::Field const&, unsigned char) const in "/opt/openfoam230/platforms/linuxGccDPOpt/lib/openmpi-system/libPstream.so"
[1] #8 Foam::PCG::solve(Foam::Field&, Foam::Field const&, unsigned char) const in "/usr/lib/libmpi.so.0"
[2] #7 Foam::reduce(double&, Foam::sumOp const&, int, int) in "/opt/openfoam230/platforms/linuxGccDPOpt/lib/libOpenFOAM.so"
[3] #9 Foam::fvMatrix::solveSegregated(Foam::dictionary const&) in "/opt/openfoam230/platforms/linuxGccDPOpt/lib/libOpenFOAM.so"
[0] #9 Foam::fvMatrix::solveSegregated(Foam::dictionary const&) in "/opt/openfoam230/platforms/linuxGccDPOpt/lib/openmpi-system/libPstream.so"
[2] #8 Foam::PCG::solve(Foam::Field&, Foam::Field const&, unsigned char) const in "/opt/openfoam230/platforms/linuxGccDPOpt/lib/libOpenFOAM.so"
[1] #9 Foam::fvMatrix::solveSegregated(Foam::dictionary const&) in "/opt/openfoam230/platforms/linuxGccDPOpt/lib/libOpenFOAM.so"
[2] #9 Foam::fvMatrix::solveSegregated(Foam::dictionary const&) in "/opt/openfoam230/platforms/linuxGccDPOpt/lib/libfiniteVolume.so"
[0] #10 in "/opt/openfoam230/platforms/linuxGccDPOpt/lib/libfiniteVolume.so"
[3] #10

in "/opt/openfoam230/platforms/linuxGccDPOpt/lib/libfiniteVolume.so"
[1] #10 [0] at cfdemSolverPiso.C:0
[0] #11
[3] at cfdemSolverPiso.C:0
[3] #11 in "/opt/openfoam230/platforms/linuxGccDPOpt/lib/libfiniteVolume.so"
[2] #10 [1] at cfdemSolverPiso.C:0
[1] #11

[0] in "/home/poseidon/OpenFOAM/poseidon-2.3.0/platforms/linuxGccDPOpt/bin/cfdemSolverPiso"
[0] #12 __libc_start_main
[3] in "/home/poseidon/OpenFOAM/poseidon-2.3.0/platforms/linuxGccDPOpt/bin/cfdemSolverPiso"
[3] #12 __libc_start_main[1] in "/home/poseidon/OpenFOAM/poseidon-2.3.0/platforms/linuxGccDPOpt/bin/cfdemSolverPiso"
[1] #12 __libc_start_main[2] at cfdemSolverPiso.C:0
[2] #11 in "/lib/i386-linux-gnu/libc.so.6"
[0] #13 in "/lib/i386-linux-gnu/libc.so.6"
[3] #13
in "/lib/i386-linux-gnu/libc.so.6"
[1] #13

[2] in "/home/poseidon/OpenFOAM/poseidon-2.3.0/platforms/linuxGccDPOpt/bin/cfdemSolverPiso"
[2] #12 __libc_start_main[3] in "/home/poseidon/OpenFOAM/poseidon-2.3.0/platforms/linuxGccDPOpt/bin/cfdemSolverPiso"
[poseidon:13091] *** Process received signal ***
[poseidon:13091] Signal: Floating point exception (8)
[poseidon:13091] Signal code: (-6)
[poseidon:13091] Failing at address: 0x3323
[poseidon:13091] [ 0] [0xb777e40c]
[poseidon:13091] [ 1] [0xb777e424]
[poseidon:13091] [ 2] /lib/i386-linux-gnu/libc.so.6(gsignal+0x4f) [0xb4b831df]
[poseidon:13091] [ 3] /opt/openfoam230/platforms/linuxGccDPOpt/lib/libOpenFOAM.so(_ZN4Foam6sigFpe10sigHandlerEi+0x60) [0xb53cb230]
[poseidon:13091] [ 4] [0xb777e400]
[poseidon:13091] [ 5] /usr/lib/libmpi.so.0(+0x26488) [0xb455b488]
[poseidon:13091] [ 6] /usr/lib/openmpi/lib/openmpi/mca_coll_tuned.so(+0x3ed3) [0xb40e6ed3]
[poseidon:13091] [ 7] /usr/lib/openmpi/lib/openmpi/mca_coll_tuned.so(+0x17ec) [0xb40e47ec]
[poseidon:13091] [ 8] /usr/lib/libmpi.so.0(PMPI_Allreduce+0x185) [0xb4578fd5]
[poseidon:13091] [ 9] /opt/openfoam230/platforms/linuxGccDPOpt/lib/openmpi-system/libPstream.so(_ZN4Foam6reduceERdRKNS_5sumOpIdEEii+0x230) [0xb4add550]
[poseidon:13091] [10] /opt/openfoam230/platforms/linuxGccDPOpt/lib/libOpenFOAM.so(_ZNK4Foam3PCG5solveERNS_5FieldIdEERKS2_h+0x7f1) [0xb5162641]
[poseidon:13091] [11] /opt/openfoam230/platforms/linuxGccDPOpt/lib/libfiniteVolume.so(_ZN4Foam8fvMatrixIdE15solveSegregatedERKNS_10dictionaryE+0x16b) [0xb67e62db]
[poseidon:13091] [12] cfdemSolverPiso() [0x8067e10]
[poseidon:13091] [13] cfdemSolverPiso() [0x8062c1a]
[poseidon:13091] [14] /lib/i386-linux-gnu/libc.so.6(__libc_start_main+0xf3) [0xb4b6e4d3]
[poseidon:13091] [15] cfdemSolverPiso() [0x8066511]
[poseidon:13091] *** End of error message ***
[1] in "/home/poseidon/OpenFOAM/poseidon-2.3.0/platforms/linuxGccDPOpt/bin/cfdemSolverPiso"
[poseidon:13089] *** Process received signal ***
[poseidon:13089] Signal: Floating point exception (8)
[poseidon:13089] Signal code: (-6)
[poseidon:13089] Failing at address: 0x3321
[poseidon:13089] [ 0] [0xb775040c]
[poseidon:13089] [ 1] [0xb7750424]
[poseidon:13089] [ 2] /lib/i386-linux-gnu/libc.so.6(gsignal+0x4f) [0xb4b551df]
[poseidon:13089] [ 3] /opt/openfoam230/platforms/linuxGccDPOpt/lib/libOpenFOAM.so(_ZN4Foam6sigFpe10sigHandlerEi+0x60) [0xb539d230]
[poseidon:13089] [ 4] [0xb7750400]
[poseidon:13089] [ 5] /usr/lib/libmpi.so.0(+0x26488) [0xb452d488]
[poseidon:13089] [ 6] /usr/lib/openmpi/lib/openmpi/mca_coll_tuned.so(+0x3ed3) [0xb40b8ed3]
[poseidon:13089] [ 7] /usr/lib/openmpi/lib/openmpi/mca_coll_tuned.so(+0x17ec) [0xb40b67ec]
[poseidon:13089] [ 8] /usr/lib/libmpi.so.0(PMPI_Allreduce+0x185) [0xb454afd5]
[poseidon:13089] [ 9] /opt/openfoam230/platforms/linuxGccDPOpt/lib/openmpi-system/libPstream.so(_ZN4Foam6reduceERdRKNS_5sumOpIdEEii+0x230) [0xb4aaf550]
[poseidon:13089] [10] /opt/openfoam230/platforms/linuxGccDPOpt/lib/libOpenFOAM.so(_ZNK4Foam3PCG5solveERNS_5FieldIdEERKS2_h+0x7f1) [0xb5134641]
[poseidon:13089] [11] /opt/openfoam230/platforms/linuxGccDPOpt/lib/libfiniteVolume.so(_ZN4Foam8fvMatrixIdE15solveSegregatedERKNS_10dictionaryE+0x16b) [0xb67b82db]
[poseidon:13089] [12] cfdemSolverPiso() [0x8067e10]
[poseidon:13089] [13] cfdemSolverPiso() [0x8062c1a]
[poseidon:13089] [14] /lib/i386-linux-gnu/libc.so.6(__libc_start_main+0xf3) [0xb4b404d3]
[poseidon:13089] [15] cfdemSolverPiso() [0x8066511]
[poseidon:13089] *** End of error message ***
[0] in "/home/poseidon/OpenFOAM/poseidon-2.3.0/platforms/linuxGccDPOpt/bin/cfdemSolverPiso"
[poseidon:13088] *** Process received signal ***
[poseidon:13088] Signal: Floating point exception (8)
[poseidon:13088] Signal code: (-6)
[poseidon:13088] Failing at address: 0x3320
[poseidon:13088] [ 0] [0xb77b340c]
[poseidon:13088] [ 1] [0xb77b3424]
[poseidon:13088] [ 2] /lib/i386-linux-gnu/libc.so.6(gsignal+0x4f) [0xb4bb81df]
[poseidon:13088] [ 3] /opt/openfoam230/platforms/linuxGccDPOpt/lib/libOpenFOAM.so(_ZN4Foam6sigFpe10sigHandlerEi+0x60) [0xb5400230]
[poseidon:13088] [ 4] [0xb77b3400]
[poseidon:13088] [ 5] /usr/lib/libmpi.so.0(+0x26488) [0xb4590488]
[poseidon:13088] [ 6] /usr/lib/openmpi/lib/openmpi/mca_coll_tuned.so(+0x3e18) [0xb411be18]
[poseidon:13088] [ 7] /usr/lib/openmpi/lib/openmpi/mca_coll_tuned.so(+0x17ec) [0xb41197ec]
[poseidon:13088] [ 8] /usr/lib/libmpi.so.0(PMPI_Allreduce+0x185) [0xb45adfd5]
[poseidon:13088] [ 9] /opt/openfoam230/platforms/linuxGccDPOpt/lib/openmpi-system/libPstream.so(_ZN4Foam6reduceERdRKNS_5sumOpIdEEii+0x230) [0xb4b12550]
[poseidon:13088] [10] /opt/openfoam230/platforms/linuxGccDPOpt/lib/libOpenFOAM.so(_ZNK4Foam3PCG5solveERNS_5FieldIdEERKS2_h+0x7f1) [0xb5197641]
[poseidon:13088] [11] /opt/openfoam230/platforms/linuxGccDPOpt/lib/libfiniteVolume.so(_ZN4Foam8fvMatrixIdE15solveSegregatedERKNS_10dictionaryE+0x16b) [0xb681b2db]
[poseidon:13088] [12] cfdemSolverPiso() [0x8067e10]
[poseidon:13088] [13] cfdemSolverPiso() [0x8062c1a]
[poseidon:13088] [14] /lib/i386-linux-gnu/libc.so.6(__libc_start_main+0xf3) [0xb4ba34d3]
[poseidon:13088] [15] cfdemSolverPiso() [0x8066511]
[poseidon:13088] *** End of error message ***
in "/lib/i386-linux-gnu/libc.so.6"
[2] #13
--------------------------------------------------------------------------
mpirun noticed that process rank 1 with PID 13089 on node poseidon exited on signal 8 (Floating point exception).
….so I can not simulate simulations that last longer than 0.07 s, but they run fine when testing the hydrodynamics only with pisoFoam solver.
Thank you very much I appreciate all help !!!

j-kerbl's picture

j-kerbl | Mon, 09/08/2014 - 12:07

Hi,

it is rather hard to see from that error message what went wrong. Could you attach a log-file of the output cfdemSolverPiso produces?
Is it the unmodified tutorial case or did you change it already to fit your needs?

Cheers
Josef

hulli | Mon, 09/08/2014 - 20:26

Hi Josef,
please find below the log file. This is one of my backup files where I at least was able to get some steps calculated. However I do not see any particle movement. I do not understand why the simulation breaks down when I simply change the simulation times in my runs. In the meantime I benchmarked nearly all parameters and got the impression that it must be related to the timing of the coupling, because I am able to run the DEM or CFD independently and they run for ages without breaking down. Thank you so much I really appreciate it !

// run_parallel_cfdemSolverPiso_ErgunTestMPI_CFDDEM //

/home/poseidon/CFDEM/poseidon-PUBLIC-2.3.0/run/ErgunTestMPI/CFD

/*---------------------------------------------------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.3.0 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
Build : 2.3.0-f5222ca19ce6
Exec : cfdemSolverPiso -parallel
Date : Sep 03 2014
Time : 16:41:33
Host : "poseidon"
PID : 30181
Case : /home/poseidon/CFDEM/poseidon-PUBLIC-2.3.0/run/ErgunTestMPI/CFD
nProcs : 4
Slaves :
3
(
"poseidon.30182"
"poseidon.30183"
"poseidon.30184"
)

Pstream initialized with:
floatTransfer : 0
nProcsSimpleSum : 0
commsType : nonBlocking
polling iterations : 0
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

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.liggghts_resume'
LIGGGHTS (Version LIGGGHTS-PUBLIC 3.0.2, compiled 2014-07-02-05:52:57 by poseidon based on LAMMPS 23 Nov 2013)
Reading restart file ...
WARNING: Restart file version does not match LIGGGHTS version (../read_restart.cpp:470)
--> restart file = Version LIGGGHTS-PUBLIC 3.0.2, compiled 2014-09-01-19:53:49 by poseidon based on LAMMPS 23 Nov 2013
--> LIGGGHTS = Version LIGGGHTS-PUBLIC 3.0.2, compiled 2014-07-02-05:52:57 by poseidon based on LAMMPS 23 Nov 2013
WARNING: Restart file used different # of processors (../read_restart.cpp:517)
WARNING: Restart file used different 3d processor grid (../read_restart.cpp:530)
orthogonal box = (-0.015 -0.015 -0.001) to (0.015 0.015 0.0554)
2 by 2 by 1 MPI processor grid
10000 atoms
Resetting global state of Fix history_zwalls1 Style property/atom from restart file info
Resetting per-atom state of Fix history_zwalls1 Style property/atom from restart file info
Resetting global state of Fix history_zwalls2 Style property/atom from restart file info
Resetting per-atom state of Fix history_zwalls2 Style property/atom from restart file info
Resetting global state of Fix history_cylwalls Style property/atom from restart file info
Resetting per-atom state of Fix history_cylwalls Style property/atom from restart file info
Setting atom values ...
10000 settings made for density
nevery as specified in LIGGGHTS is overriden by calling external program (../cfd_datacoupling_mpi.cpp:57)
WARNING: Dump mesh/stl cannot find any fix of type 'mesh/surface' to dump (../dump_mesh_stl.cpp:123)
Resetting global state of Fix contacthistory Style contacthistory from restart file info
WARNING: Dump mesh/stl cannot find any fix of type 'mesh/surface' to dump (../dump_mesh_stl.cpp:123)
Resetting per-atom state of Fix contacthistory Style contacthistory from restart file info
WARNING: Dump mesh/stl cannot find any fix of type 'mesh/surface' to dump (../dump_mesh_stl.cpp:123)
WARNING: Dump mesh/stl cannot find any fix of type 'mesh/surface' to dump (../dump_mesh_stl.cpp:123)
Setting up run ...
Memory usage per processor = 10.1338 Mbytes
Step Atoms KinEng 1 Volume centerOf dragtota dragtota dragtota
10000 10000 9.4102837e-08 8.1266333e-08 5.076e-05 0.0075829221 0 0 0
10001 10000 9.4018037e-08 8.1217291e-08 5.076e-05 0.0075829201 0 0 0
Loop time of 0.00827169 on 4 procs for 1 steps with 10000 atoms

Pair time (%) = 0.000979304 (11.8392)
Neigh time (%) = 0 (0)
Comm time (%) = 0.0019486 (23.5574)
Outpt time (%) = 0.00494307 (59.7589)
Other time (%) = 0.000400722 (4.8445)

Nlocal: 2500 ave 2522 max 2490 min
Histogram: 2 0 1 0 0 0 0 0 0 1
Nghost: 734.5 ave 756 max 709 min
Histogram: 1 0 0 0 1 0 0 1 0 1
Neighs: 17447.8 ave 17591 max 17363 min
Histogram: 1 1 0 1 0 0 0 0 0 1

Total # of neighbors = 69791
Ave neighs/atom = 6.9791
Neighbor list builds = 0
Dangerous builds = 0
Selecting IOModel basicIO
Selecting probeModel off
Selecting voidFractionModel divided
Selecting averagingModel dense
Selecting clockModel standardClock
start clock measurement at t >0.0005
Selecting smoothingModel off
Selecting meshMotionModel noMeshMotion

CFDEMcoupling version: cfdem-2.7.1

, compatible to LIGGGHTS version: 3.0.2

, compatible to OF version: 2.3.x-commit-4d6f4a3115ff76ec4154c580eb041bc95ba4ec09

If BC are important, please provide volScalarFields -imp/expParticleForces-
ignoring ddt(voidfraction)
Selecting forceModel KochHillDrag
Selecting forceModel gradPForce
gradPForce is applied only to DEM side
using interpolated value of pressure gradient.
Selecting forceModel viscForce
viscForce is applied only to DEM side
using interpolated value of pressure gradient.
Selecting momCoupleModel implicitCouple

implicit momentum exchange field calculate if alphaP larger than : 1
Selecting liggghtsCommandModel runLiggghts ,provide dicts, numbered from 0 to n
cg is set to: 1
solving volume averaged Navier Stokes equations of type A
check if Archimedes lift is necessary (depends on pressure formulation)

Starting time loop

Time = 0.0005

Courant Number mean: 0.000904159 max: 0.0216998

timeStepFraction() = 0.25
update Ksl.internalField()
DILUPBiCG: Solving for Ux, Initial residual = 0.0667276, Final residual = 5.21943e-06, No Iterations 2
DILUPBiCG: Solving for Uy, Initial residual = 0.0757376, Final residual = 7.33718e-06, No Iterations 2
DILUPBiCG: Solving for Uz, Initial residual = 1, Final residual = 1.51776e-06, No Iterations 3
suppressing ddt(voidfraction)
DICPCG: Solving for p, Initial residual = 1, Final residual = 0.0948392, No Iterations 26
suppressing ddt(voidfraction)
DICPCG: Solving for p, Initial residual = 0.0018859, Final residual = 0.000152842, No Iterations 7
time step continuity errors : sum local = 1.39912e-05, global = -6.08989e-07, cumulative = -6.08989e-07
suppressing ddt(voidfraction)
DICPCG: Solving for p, Initial residual = 0.00266956, Final residual = 0.000193059, No Iterations 26
suppressing ddt(voidfraction)
DICPCG: Solving for p, Initial residual = 0.000235493, Final residual = 1.63875e-05, No Iterations 7
time step continuity errors : sum local = 1.32659e-06, global = 5.48844e-08, cumulative = -5.54105e-07
suppressing ddt(voidfraction)
DICPCG: Solving for p, Initial residual = 9.76065e-05, Final residual = 9.40296e-06, No Iterations 25
suppressing ddt(voidfraction)
DICPCG: Solving for p, Initial residual = 1.14044e-05, Final residual = 1.13792e-06, No Iterations 5
time step continuity errors : sum local = 9.20556e-08, global = 5.18716e-09, cumulative = -5.48918e-07
suppressing ddt(voidfraction)
DICPCG: Solving for p, Initial residual = 5.05308e-06, Final residual = 9.54519e-07, No Iterations 24
suppressing ddt(voidfraction)
DICPCG: Solving for p, Initial residual = 1.07054e-06, Final residual = 5.16025e-07, No Iterations 1
time step continuity errors : sum local = 4.17431e-08, global = 4.80791e-09, cumulative = -5.4411e-07
suppressing ddt(voidfraction)
DICPCG: Solving for p, Initial residual = 7.81263e-07, Final residual = 7.81263e-07, No Iterations 0
suppressing ddt(voidfraction)
DICPCG: Solving for p, Initial residual = 7.81263e-07, Final residual = 7.81263e-07, No Iterations 0
time step continuity errors : sum local = 6.31992e-08, global = 1.19325e-08, cumulative = -5.32177e-07
suppressing ddt(voidfraction)
DICPCG: Solving for p, Initial residual = 8.30306e-07, Final residual = 8.30306e-07, No Iterations 0
suppressing ddt(voidfraction)
DICPCG: Solving for p, Initial residual = 8.30306e-07, Final residual = 8.30306e-07, No Iterations 0
time step continuity errors : sum local = 6.71665e-08, global = 1.36465e-08, cumulative = -5.18531e-07
suppressing ddt(voidfraction)
DICPCG: Solving for p, Initial residual = 8.37179e-07, Final residual = 8.37179e-07, No Iterations 0
suppressing ddt(voidfraction)
DICPCG: Solving for p, Initial residual = 8.37179e-07, Final residual = 8.37179e-07, No Iterations 0
time step continuity errors : sum local = 6.77225e-08, global = 1.39453e-08, cumulative = -5.04585e-07
ExecutionTime = 0.11 s ClockTime = 1 s

Starting time loop

Time = 0.001

Courant Number mean: 0.00583949 max: 0.025538

timeStepFraction() = 0.5
update Ksl.internalField()
DILUPBiCG: Solving for Ux, Initial residual = 0.0247312, Final residual = 8.04009e-07, No Iterations 2
DILUPBiCG: Solving for Uy, Initial residual = 0.0247434, Final residual = 8.01159e-07, No Iterations 2
DILUPBiCG: Solving for Uz, Initial residual = 0.0628047, Final residual = 6.74471e-06, No Iterations 2
suppressing ddt(voidfraction)
DICPCG: Solving for p, Initial residual = 0.0204451, Final residual = 0.00168776, No Iterations 26
suppressing ddt(voidfraction)
DICPCG: Solving for p, Initial residual = 0.0963279, Final residual = 0.00796892, No Iterations 7
time step continuity errors : sum local = 1.12115e-05, global = 4.13178e-07, cumulative = -9.14073e-08
suppressing ddt(voidfraction)
DICPCG: Solving for p, Initial residual = 0.0356656, Final residual = 0.00327296, No Iterations 7
suppressing ddt(voidfraction)
DICPCG: Solving for p, Initial residual = 0.00371321, Final residual = 0.000301265, No Iterations 24
time step continuity errors : sum local = 3.74038e-07, global = -2.73887e-09, cumulative = -9.41462e-08
suppressing ddt(voidfraction)
DICPCG: Solving for p, Initial residual = 0.00190124, Final residual = 0.00014793, No Iterations 6
suppressing ddt(voidfraction)
DICPCG: Solving for p, Initial residual = 0.000178641, Final residual = 1.39498e-05, No Iterations 13
time step continuity errors : sum local = 1.7102e-08, global = -3.98603e-09, cumulative = -9.81322e-08
suppressing ddt(voidfraction)
DICPCG: Solving for p, Initial residual = 9.68504e-05, Final residual = 9.11995e-06, No Iterations 4
suppressing ddt(voidfraction)
DICPCG: Solving for p, Initial residual = 1.06908e-05, Final residual = 8.61115e-07, No Iterations 26
time step continuity errors : sum local = 1.05548e-09, global = -7.84144e-11, cumulative = -9.82106e-08
suppressing ddt(voidfraction)
DICPCG: Solving for p, Initial residual = 6.2111e-06, Final residual = 8.5677e-07, No Iterations 3
suppressing ddt(voidfraction)
DICPCG: Solving for p, Initial residual = 9.32623e-07, Final residual = 9.32623e-07, No Iterations 0
time step continuity errors : sum local = 1.14324e-09, global = 3.62607e-10, cumulative = -9.7848e-08
suppressing ddt(voidfraction)
DICPCG: Solving for p, Initial residual = 1.13108e-06, Final residual = 6.46992e-07, No Iterations 1
suppressing ddt(voidfraction)
DICPCG: Solving for p, Initial residual = 6.46451e-07, Final residual = 6.46451e-07, No Iterations 0
time step continuity errors : sum local = 7.92445e-10, global = 4.07323e-10, cumulative = -9.74407e-08
ExecutionTime = 0.15 s ClockTime = 1 s

Starting time loop

Time = 0.0015

Courant Number mean: 0.0431084 max: 0.0434284

timeStepFraction() = 0.75
update Ksl.internalField()
DILUPBiCG: Solving for Ux, Initial residual = 0.073771, Final residual = 2.96044e-06, No Iterations 2
DILUPBiCG: Solving for Uy, Initial residual = 0.0737709, Final residual = 2.9719e-06, No Iterations 2
DILUPBiCG: Solving for Uz, Initial residual = 0.0783965, Final residual = 8.60726e-06, No Iterations 2
suppressing ddt(voidfraction)
DICPCG: Solving for p, Initial residual = 0.288544, Final residual = 0.0272918, No Iterations 6
suppressing ddt(voidfraction)
DICPCG: Solving for p, Initial residual = 0.0334503, Final residual = 0.00267758, No Iterations 15
time step continuity errors : sum local = 2.13726e-06, global = 8.58242e-07, cumulative = 7.60801e-07
suppressing ddt(voidfraction)
DICPCG: Solving for p, Initial residual = 0.0205541, Final residual = 0.00161566, No Iterations 20
suppressing ddt(voidfraction)
DICPCG: Solving for p, Initial residual = 0.00176087, Final residual = 0.000151493, No Iterations 10
time step continuity errors : sum local = 1.11357e-07, global = 1.85896e-08, cumulative = 7.79391e-07
suppressing ddt(voidfraction)
DICPCG: Solving for p, Initial residual = 0.00123283, Final residual = 8.78163e-05, No Iterations 23
suppressing ddt(voidfraction)
DICPCG: Solving for p, Initial residual = 9.43201e-05, Final residual = 8.69712e-06, No Iterations 8
time step continuity errors : sum local = 6.40606e-09, global = -1.01465e-10, cumulative = 7.79289e-07
suppressing ddt(voidfraction)
DICPCG: Solving for p, Initial residual = 6.76292e-05, Final residual = 5.60665e-06, No Iterations 18
suppressing ddt(voidfraction)
DICPCG: Solving for p, Initial residual = 6.2621e-06, Final residual = 9.11889e-07, No Iterations 5
time step continuity errors : sum local = 6.71716e-10, global = -1.33258e-10, cumulative = 7.79156e-07
suppressing ddt(voidfraction)
DICPCG: Solving for p, Initial residual = 4.43325e-06, Final residual = 9.00162e-07, No Iterations 3
suppressing ddt(voidfraction)
DICPCG: Solving for p, Initial residual = 9.44394e-07, Final residual = 9.44394e-07, No Iterations 0
time step continuity errors : sum local = 6.9562e-10, global = -1.51335e-10, cumulative = 7.79005e-07
ExecutionTime = 0.19 s ClockTime = 1 s

Starting time loop

Time = 0.002

Courant Number mean: 0.0433244 max: 0.044053
Starting up LIGGGHTS
Executing command: 'run 200 '
Setting up run ...
Memory usage per processor = 10.1341 Mbytes
Step Atoms KinEng 1 Volume centerOf dragtota dragtota dragtota
10001 10000 9.4018037e-08 8.1217291e-08 5.076e-05 0.0075829201 0 0 0
10010 10000 9.4176164e-08 8.1256577e-08 5.076e-05 0.0075829027 0 0 0
10020 10000 9.4610531e-08 8.1328049e-08 5.076e-05 0.0075828838 0 0 0
10030 10000 9.4463571e-08 8.1305261e-08 5.076e-05 0.0075828654 0 0 0
10040 10000 9.4615162e-08 8.0814247e-08 5.076e-05 0.0075828475 0 0 0
10050 10000 9.5022239e-08 8.0361739e-08 5.076e-05 0.0075828302 0 0 0
10060 10000 9.5195197e-08 8.0117333e-08 5.076e-05 0.0075828134 0 0 0
10070 10000 9.559187e-08 7.9978117e-08 5.076e-05 0.0075827971 0 0 0
10080 10000 9.5622062e-08 7.9418016e-08 5.076e-05 0.0075827814 0 0 0
10090 10000 9.5641847e-08 7.9111929e-08 5.076e-05 0.0075827663 0 0 0
CFD Coupling established at step 10100
10100 10000 9.5907089e-08 7.8898709e-08 5.076e-05 0.0075827516 0 0 0
10110 10000 9.6130655e-08 7.8687735e-08 5.076e-05 0.0075827374 0 0 0
10120 10000 9.6278054e-08 7.8257395e-08 5.076e-05 0.0075827237 0 0 0
10130 10000 9.6309129e-08 7.8030132e-08 5.076e-05 0.0075827104 0 0 0
10140 10000 9.6735191e-08 7.7876652e-08 5.076e-05 0.0075826975 0 0 0
10150 10000 9.7331625e-08 7.7753438e-08 5.076e-05 0.007582685 0 0 0
10160 10000 9.7890128e-08 7.7710797e-08 5.076e-05 0.0075826728 0 0 0
10170 10000 9.8437661e-08 7.757888e-08 5.076e-05 0.007582661 0 0 0
10180 10000 9.8197666e-08 7.6656342e-08 5.076e-05 0.0075826496 0 0 0
10190 10000 9.8481757e-08 7.6506127e-08 5.076e-05 0.0075826386 0 0 0
CFD Coupling established at step 10200
10200 10000 9.833396e-08 7.6005531e-08 5.076e-05 0.0075826281 0 0 0
10201 10000 9.8360816e-08 7.5976959e-08 5.076e-05 0.0075826271 0 0 0
Loop time of 1.28183 on 4 procs for 200 steps with 10000 atoms

Pair time (%) = 0.215467 (16.8093)
Neigh time (%) = 0.0159923 (1.24761)
Comm time (%) = 0.505157 (39.4091)
Outpt time (%) = 0.0692928 (5.40578)
Other time (%) = 0.475919 (37.1282)

Nlocal: 2500 ave 2522 max 2490 min
Histogram: 2 0 1 0 0 0 0 0 0 1
Nghost: 734.5 ave 756 max 709 min
Histogram: 1 0 0 0 1 0 0 1 0 1
Neighs: 17450.8 ave 17595 max 17366 min
Histogram: 1 1 0 1 0 0 0 0 0 1

Total # of neighbors = 69803
Ave neighs/atom = 6.9803
Neighbor list builds = 2
Dangerous builds = 0
LIGGGHTS finished

Coupling...

timeStepFraction() = 1
update Ksl.internalField()
DILUPBiCG: Solving for Ux, Initial residual = 0.126159, Final residual = 6.93304e-06, No Iterations 2
DILUPBiCG: Solving for Uy, Initial residual = 0.126576, Final residual = 6.69229e-06, No Iterations 2
DILUPBiCG: Solving for Uz, Initial residual = 0.131515, Final residual = 3.11751e-07, No Iterations 3
suppressing ddt(voidfraction)
DICPCG: Solving for p, Initial residual = 0.923966, Final residual = 0.0923256, No Iterations 18
suppressing ddt(voidfraction)
DICPCG: Solving for p, Initial residual = 0.024663, Final residual = 0.0018473, No Iterations 23
time step continuity errors : sum local = 1.71327e-05, global = -2.31868e-07, cumulative = 5.47137e-07
suppressing ddt(voidfraction)
DICPCG: Solving for p, Initial residual = 0.0295681, Final residual = 0.0016789, No Iterations 25
suppressing ddt(voidfraction)
DICPCG: Solving for p, Initial residual = 0.00203084, Final residual = 0.000197223, No Iterations 6
time step continuity errors : sum local = 2.26485e-06, global = 2.84716e-07, cumulative = 8.31854e-07
suppressing ddt(voidfraction)
DICPCG: Solving for p, Initial residual = 0.00117286, Final residual = 0.000113826, No Iterations 16
suppressing ddt(voidfraction)
DICPCG: Solving for p, Initial residual = 0.000127847, Final residual = 9.23984e-06, No Iterations 20
time step continuity errors : sum local = 1.05859e-07, global = -3.16994e-09, cumulative = 8.28684e-07
suppressing ddt(voidfraction)
DICPCG: Solving for p, Initial residual = 7.27518e-05, Final residual = 6.60729e-06, No Iterations 12
suppressing ddt(voidfraction)
DICPCG: Solving for p, Initial residual = 8.33729e-06, Final residual = 8.20613e-07, No Iterations 23
time step continuity errors : sum local = 9.39093e-09, global = -2.16363e-11, cumulative = 8.28662e-07
ExecutionTime = 0.63 s ClockTime = 3 s

Starting time loop

Time = 0.0025

Courant Number mean: 0.0380539 max: 0.0439754

timeStepFraction() = 0.25
update Ksl.internalField()
DILUPBiCG: Solving for Ux, Initial residual = 0.552074, Final residual = 7.94001e-06, No Iterations 2
DILUPBiCG: Solving for Uy, Initial residual = 0.567565, Final residual = 8.41684e-06, No Iterations 2
DILUPBiCG: Solving for Uz, Initial residual = 0.762603, Final residual = 6.27591e-06, No Iterations 2
suppressing ddt(voidfraction)
DICPCG: Solving for p, Initial residual = 0.0962798, Final residual = 0.00846693, No Iterations 18
suppressing ddt(voidfraction)
DICPCG: Solving for p, Initial residual = 0.00640221, Final residual = 0.000621111, No Iterations 22
time step continuity errors : sum local = 9.05222e-06, global = -4.65951e-07, cumulative = 3.62711e-07
suppressing ddt(voidfraction)
DICPCG: Solving for p, Initial residual = 0.00492205, Final residual = 0.000300325, No Iterations 25
suppressing ddt(voidfraction)
DICPCG: Solving for p, Initial residual = 0.000340946, Final residual = 1.94393e-05, No Iterations 7
time step continuity errors : sum local = 3.0284e-07, global = 7.41712e-08, cumulative = 4.36882e-07
suppressing ddt(voidfraction)
DICPCG: Solving for p, Initial residual = 0.000145077, Final residual = 1.30596e-05, No Iterations 9
suppressing ddt(voidfraction)
DICPCG: Solving for p, Initial residual = 1.50116e-05, Final residual = 1.36102e-06, No Iterations 23
time step continuity errors : sum local = 2.11984e-08, global = 6.52626e-10, cumulative = 4.37535e-07
suppressing ddt(voidfraction)
DICPCG: Solving for p, Initial residual = 6.28261e-06, Final residual = 8.30612e-07, No Iterations 19
suppressing ddt(voidfraction)
DICPCG: Solving for p, Initial residual = 9.83973e-07, Final residual = 9.83973e-07, No Iterations 0
time step continuity errors : sum local = 1.5324e-08, global = 1.5101e-09, cumulative = 4.39045e-07
suppressing ddt(voidfraction)
DICPCG: Solving for p, Initial residual = 1.11819e-06, Final residual = 6.65748e-07, No Iterations 1
suppressing ddt(voidfraction)
DICPCG: Solving for p, Initial residual = 6.67881e-07, Final residual = 6.67881e-07, No Iterations 0
time step continuity errors : sum local = 1.04013e-08, global = 1.32153e-09, cumulative = 4.40366e-07
suppressing ddt(voidfraction)
DICPCG: Solving for p, Initial residual = 6.76206e-07, Final residual = 6.76206e-07, No Iterations 0
suppressing ddt(voidfraction)
DICPCG: Solving for p, Initial residual = 6.76206e-07, Final residual = 6.76206e-07, No Iterations 0
time step continuity errors : sum local = 1.0531e-08, global = 1.29751e-09, cumulative = 4.41664e-07
suppressing ddt(voidfraction)
DICPCG: Solving for p, Initial residual = 6.77494e-07, Final residual = 6.77494e-07, No Iterations 0
suppressing ddt(voidfraction)
DICPCG: Solving for p, Initial residual = 6.77494e-07, Final residual = 6.77494e-07, No Iterations 0
time step continuity errors : sum local = 1.0551e-08, global = 1.29374e-09, cumulative = 4.42958e-07
ExecutionTime = 0.7 s ClockTime = 3 s

Starting time loop

Time = 0.003

Courant Number mean: 0.0384664 max: 0.0479352

timeStepFraction() = 0.5
update Ksl.internalField()
DILUPBiCG: Solving for Ux, Initial residual = 0.699066, Final residual = 2.79663e-06, No Iterations 2
DILUPBiCG: Solving for Uy, Initial residual = 0.705729, Final residual = 3.19124e-06, No Iterations 2
DILUPBiCG: Solving for Uz, Initial residual = 0.930427, Final residual = 4.03897e-06, No Iterations 2
suppressing ddt(voidfraction)
DICPCG: Solving for p, Initial residual = 0.163906, Final residual = 0.0110152, No Iterations 25
suppressing ddt(voidfraction)
DICPCG: Solving for p, Initial residual = 0.00321408, Final residual = 0.000238419, No Iterations 7
time step continuity errors : sum local = 1.37535e-05, global = -3.83531e-06, cumulative = -3.39235e-06
suppressing ddt(voidfraction)
DICPCG: Solving for p, Initial residual = 0.0032055, Final residual = 0.000293893, No Iterations 14
suppressing ddt(voidfraction)
DICPCG: Solving for p, Initial residual = 0.000306236, Final residual = 2.69609e-05, No Iterations 18
time step continuity errors : sum local = 1.53125e-06, global = 1.5102e-07, cumulative = -3.24133e-06
suppressing ddt(voidfraction)
DICPCG: Solving for p, Initial residual = 0.000107608, Final residual = 9.71187e-06, No Iterations 19
suppressing ddt(voidfraction)
DICPCG: Solving for p, Initial residual = 1.10265e-05, Final residual = 7.99136e-07, No Iterations 23
time step continuity errors : sum local = 4.52823e-08, global = 3.7192e-09, cumulative = -3.23761e-06
suppressing ddt(voidfraction)
DICPCG: Solving for p, Initial residual = 5.34748e-06, Final residual = 9.78494e-07, No Iterations 5
suppressing ddt(voidfraction)
DICPCG: Solving for p, Initial residual = 1.05932e-06, Final residual = 6.16257e-07, No Iterations 1
time step continuity errors : sum local = 3.4918e-08, global = -8.52179e-09, cumulative = -3.24613e-06
suppressing ddt(voidfraction)
DICPCG: Solving for p, Initial residual = 7.75233e-07, Final residual = 7.75233e-07, No Iterations 0
suppressing ddt(voidfraction)
DICPCG: Solving for p, Initial residual = 7.75233e-07, Final residual = 7.75233e-07, No Iterations 0
time step continuity errors : sum local = 4.39257e-08, global = -1.02239e-08, cumulative = -3.25636e-06
suppressing ddt(voidfraction)
DICPCG: Solving for p, Initial residual = 7.96042e-07, Final residual = 7.96042e-07, No Iterations 0
suppressing ddt(voidfraction)
DICPCG: Solving for p, Initial residual = 7.96042e-07, Final residual = 7.96042e-07, No Iterations 0
time step continuity errors : sum local = 4.51049e-08, global = -1.04429e-08, cumulative = -3.2668e-06
ExecutionTime = 0.76 s ClockTime = 4 s

Starting time loop

Time = 0.0035

Courant Number mean: 0.0338913 max: 0.0484831

timeStepFraction() = 0.75
update Ksl.internalField()
DILUPBiCG: Solving for Ux, Initial residual = 0.79977, Final residual = 1.21753e-06, No Iterations 2
DILUPBiCG: Solving for Uy, Initial residual = 0.807072, Final residual = 1.38985e-06, No Iterations 2
DILUPBiCG: Solving for Uz, Initial residual = 0.915679, Final residual = 2.07722e-06, No Iterations 2
suppressing ddt(voidfraction)
DICPCG: Solving for p, Initial residual = 0.0401306, Final residual = 0.00257249, No Iterations 25
suppressing ddt(voidfraction)
DICPCG: Solving for p, Initial residual = 0.00127819, Final residual = 0.000100209, No Iterations 7
time step continuity errors : sum local = 1.22596e-05, global = -4.46137e-06, cumulative = -7.72816e-06
suppressing ddt(voidfraction)
DICPCG: Solving for p, Initial residual = 0.00131329, Final residual = 9.8475e-05, No Iterations 13
suppressing ddt(voidfraction)
DICPCG: Solving for p, Initial residual = 0.000103265, Final residual = 9.96352e-06, No Iterations 9
time step continuity errors : sum local = 1.20864e-06, global = 8.58029e-08, cumulative = -7.64236e-06
suppressing ddt(voidfraction)
DICPCG: Solving for p, Initial residual = 4.43452e-05, Final residual = 4.20103e-06, No Iterations 18
suppressing ddt(voidfraction)
DICPCG: Solving for p, Initial residual = 4.59264e-06, Final residual = 9.4796e-07, No Iterations 6
time step continuity errors : sum local = 1.14968e-07, global = 2.02349e-08, cumulative = -7.62213e-06
suppressing ddt(voidfraction)
DICPCG: Solving for p, Initial residual = 2.40082e-06, Final residual = 8.68508e-07, No Iterations 2
suppressing ddt(voidfraction)
DICPCG: Solving for p, Initial residual = 8.75651e-07, Final residual = 8.75651e-07, No Iterations 0
time step continuity errors : sum local = 1.06199e-07, global = 1.48856e-09, cumulative = -7.62064e-06
suppressing ddt(voidfraction)
DICPCG: Solving for p, Initial residual = 8.82991e-07, Final residual = 8.82991e-07, No Iterations 0
suppressing ddt(voidfraction)
DICPCG: Solving for p, Initial residual = 8.82991e-07, Final residual = 8.82991e-07, No Iterations 0
time step continuity errors : sum local = 1.07089e-07, global = -6.91649e-10, cumulative = -7.62133e-06
ExecutionTime = 0.81 s ClockTime = 4 s

Starting time loop

Time = 0.004

Courant Number mean: 0.0334389 max: 0.0553343
Starting up LIGGGHTS
Executing command: 'run 200 pre no post no'
Step Atoms KinEng 1 Volume centerOf dragtota dragtota dragtota
10201 10000 9.8360816e-08 7.5976959e-08 5.076e-05 0.0075826271 9.3304136e-05 -6.6351777e-05 0.34857644
10210 10000 1.4040492e-07 7.5963471e-08 5.076e-05 0.007582736 9.3304136e-05 -6.6351777e-05 0.34857644
10220 10000 2.9135778e-07 7.7173332e-08 5.076e-05 0.0075831586 9.3304136e-05 -6.6351777e-05 0.34857644
10230 10000 5.3728733e-07 8.0197319e-08 5.076e-05 0.0075838819 9.3304136e-05 -6.6351777e-05 0.34857644
10240 10000 8.6611509e-07 8.5391866e-08 5.076e-05 0.00758489 9.3304136e-05 -6.6351777e-05 0.34857644
10250 10000 1.243949e-06 7.8893965e-08 5.076e-05 0.0075861691 9.3304136e-05 -6.6351777e-05 0.34857644
10260 10000 1.6991422e-06 8.6145517e-08 5.076e-05 0.007587707 9.3304136e-05 -6.6351777e-05 0.34857644
10270 10000 2.2185131e-06 9.8654202e-08 5.076e-05 0.0075894935 9.3304136e-05 -6.6351777e-05 0.34857644
10280 10000 2.7921873e-06 1.1398369e-07 5.076e-05 0.0075915197 9.3304136e-05 -6.6351777e-05 0.34857644
10290 10000 3.4178861e-06 1.3123593e-07 5.076e-05 0.0075937787 9.3304136e-05 -6.6351777e-05 0.34857644
CFD Coupling established at step 10300
10300 10000 4.092993e-06 1.5025003e-07 5.076e-05 0.0075962647 9.3304136e-05 -6.6351777e-05 0.34857644
10310 10000 4.8167129e-06 1.6978662e-07 5.076e-05 0.0075989736 9.3304136e-05 -6.6351777e-05 0.34857644
10320 10000 5.5889693e-06 1.8915659e-07 5.076e-05 0.0076019022 9.3304136e-05 -6.6351777e-05 0.34857644
10330 10000 6.4092489e-06 2.0759449e-07 5.076e-05 0.0076050486 9.3304136e-05 -6.6351777e-05 0.34857644
10340 10000 7.277966e-06 2.2514771e-07 5.076e-05 0.0076084116 9.3304136e-05 -6.6351777e-05 0.34857644
10350 10000 8.1978522e-06 2.4381778e-07 5.076e-05 0.0076119907 9.3304136e-05 -6.6351777e-05 0.34857644
10360 10000 9.1688381e-06 2.6225563e-07 5.076e-05 0.007615786 9.3304136e-05 -6.6351777e-05 0.34857644
10370 10000 1.0191898e-05 2.8039326e-07 5.076e-05 0.0076197979 9.3304136e-05 -6.6351777e-05 0.34857644
10380 10000 1.1269978e-05 2.9878294e-07 5.076e-05 0.0076240268 9.3304136e-05 -6.6351777e-05 0.34857644
10390 10000 1.240344e-05 3.1733062e-07 5.076e-05 0.0076284737 9.3304136e-05 -6.6351777e-05 0.34857644
CFD Coupling established at step 10400
10400 10000 1.3593407e-05 3.3509534e-07 5.076e-05 0.0076331394 9.3304136e-05 -6.6351777e-05 0.34857644
10401 10000 1.371544e-05 3.3678412e-07 5.076e-05 0.0076336181 9.3304136e-05 -6.6351777e-05 0.34857644
Loop time of 1.24455 on 4 procs for 200 steps with 10000 atoms
LIGGGHTS finished

Coupling...

timeStepFraction() = 1
update Ksl.internalField()
DILUPBiCG: Solving for Ux, Initial residual = 0.843951, Final residual = 5.89195e-07, No Iterations 2
DILUPBiCG: Solving for Uy, Initial residual = 0.851278, Final residual = 6.4595e-07, No Iterations 2
DILUPBiCG: Solving for Uz, Initial residual = 0.927633, Final residual = 1.19818e-06, No Iterations 2
suppressing ddt(voidfraction)
DICPCG: Solving for p, Initial residual = 0.0184113, Final residual = 0.00131549, No Iterations 25
suppressing ddt(voidfraction)
DICPCG: Solving for p, Initial residual = 0.000661022, Final residual = 6.132e-05, No Iterations 6
time step continuity errors : sum local = 0.0398332, global = -0.0398245, cumulative = -0.0398322
suppressing ddt(voidfraction)
DICPCG: Solving for p, Initial residual = 0.00057953, Final residual = 5.36262e-05, No Iterations 23
suppressing ddt(voidfraction)
DICPCG: Solving for p, Initial residual = 5.50728e-05, Final residual = 5.07195e-06, No Iterations 7
time step continuity errors : sum local = 0.039821, global = -0.0398195, cumulative = -0.0796517
suppressing ddt(voidfraction)
DICPCG: Solving for p, Initial residual = 2.10759e-05, Final residual = 1.35184e-06, No Iterations 19
suppressing ddt(voidfraction)
DICPCG: Solving for p, Initial residual = 1.58732e-06, Final residual = 9.22573e-07, No Iterations 1
time step continuity errors : sum local = 0.0398201, global = -0.0398198, cumulative = -0.119472
suppressing ddt(voidfraction)
DICPCG: Solving for p, Initial residual = 1.40219e-06, Final residual = 9.58386e-07, No Iterations 1
suppressing ddt(voidfraction)
DICPCG: Solving for p, Initial residual = 9.60766e-07, Final residual = 9.60766e-07, No Iterations 0
time step continuity errors : sum local = 0.0398201, global = -0.0398199, cumulative = -0.159291
ExecutionTime = 1.2 s ClockTime = 6 s

Starting time loop

Time = 0.0045

Courant Number mean: 0.0336723 max: 0.0609544

timeStepFraction() = 0.25
update Ksl.internalField()
DILUPBiCG: Solving for Ux, Initial residual = 0.880891, Final residual = 5.23238e-07, No Iterations 2
DILUPBiCG: Solving for Uy, Initial residual = 0.886507, Final residual = 5.47076e-07, No Iterations 2
DILUPBiCG: Solving for Uz, Initial residual = 0.927454, Final residual = 9.36911e-07, No Iterations 2
suppressing ddt(voidfraction)
DICPCG: Solving for p, Initial residual = 0.00300589, Final residual = 0.000288095, No Iterations 23
suppressing ddt(voidfraction)
DICPCG: Solving for p, Initial residual = 0.00029603, Final residual = 2.19779e-05, No Iterations 7
time step continuity errors : sum local = 0.000969545, global = -6.28741e-08, cumulative = -0.159291
suppressing ddt(voidfraction)
DICPCG: Solving for p, Initial residual = 0.000542116, Final residual = 3.07743e-05, No Iterations 23
suppressing ddt(voidfraction)
DICPCG: Solving for p, Initial residual = 3.21364e-05, Final residual = 2.65512e-06, No Iterations 6
time step continuity errors : sum local = 0.000965097, global = -6.58992e-09, cumulative = -0.159291
suppressing ddt(voidfraction)
DICPCG: Solving for p, Initial residual = 1.96895e-05, Final residual = 1.83815e-06, No Iterations 10
suppressing ddt(voidfraction)
DICPCG: Solving for p, Initial residual = 1.93694e-06, Final residual = 8.74695e-07, No Iterations 13
time step continuity errors : sum local = 0.00096468, global = 3.29739e-08, cumulative = -0.159291
suppressing ddt(voidfraction)
DICPCG: Solving for p, Initial residual = 1.2776e-06, Final residual = 6.02596e-07, No Iterations 1
suppressing ddt(voidfraction)
DICPCG: Solving for p, Initial residual = 6.01597e-07, Final residual = 6.01597e-07, No Iterations 0
time step continuity errors : sum local = 0.000964613, global = 3.01445e-08, cumulative = -0.159291
suppressing ddt(voidfraction)
DICPCG: Solving for p, Initial residual = 6.2454e-07, Final residual = 6.2454e-07, No Iterations 0
suppressing ddt(voidfraction)
DICPCG: Solving for p, Initial residual = 6.2454e-07, Final residual = 6.2454e-07, No Iterations 0
time step continuity errors : sum local = 0.000964617, global = 2.93621e-08, cumulative = -0.159291
suppressing ddt(voidfraction)
DICPCG: Solving for p, Initial residual = 6.26617e-07, Final residual = 6.26617e-07, No Iterations 0
suppressing ddt(voidfraction)
DICPCG: Solving for p, Initial residual = 6.26617e-07, Final residual = 6.26617e-07, No Iterations 0
time step continuity errors : sum local = 0.000964618, global = 2.92049e-08, cumulative = -0.159291
suppressing ddt(voidfraction)
DICPCG: Solving for p, Initial residual = 6.27017e-07, Final residual = 6.27017e-07, No Iterations 0
suppressing ddt(voidfraction)
DICPCG: Solving for p, Initial residual = 6.27017e-07, Final residual = 6.27017e-07, No Iterations 0
time step continuity errors : sum local = 0.000964618, global = 2.91773e-08, cumulative = -0.159291
ExecutionTime = 1.27 s ClockTime = 6 s

Starting time loop

Time = 0.005

Courant Number mean: 0.0357893 max: 0.06601

timeStepFraction() = 0.5
update Ksl.internalField()
DILUPBiCG: Solving for Ux, Initial residual = 0.877061, Final residual = 5.43118e-07, No Iterations 2
DILUPBiCG: Solving for Uy, Initial residual = 0.887549, Final residual = 5.6082e-07, No Iterations 2
DILUPBiCG: Solving for Uz, Initial residual = 0.920607, Final residual = 9.76777e-07, No Iterations 2
suppressing ddt(voidfraction)
DICPCG: Solving for p, Initial residual = 0.00163281, Final residual = 0.00014086, No Iterations 20
suppressing ddt(voidfraction)
DICPCG: Solving for p, Initial residual = 0.000148121, Final residual = 1.45277e-05, No Iterations 10
time step continuity errors : sum local = 0.000967783, global = -1.4135e-06, cumulative = -0.159293
suppressing ddt(voidfraction)
DICPCG: Solving for p, Initial residual = 0.000514841, Final residual = 4.71435e-05, No Iterations 14
suppressing ddt(voidfraction)
DICPCG: Solving for p, Initial residual = 4.81931e-05, Final residual = 3.84456e-06, No Iterations 19
time step continuity errors : sum local = 0.000965355, global = -3.0548e-08, cumulative = -0.159293
suppressing ddt(voidfraction)
DICPCG: Solving for p, Initial residual = 1.92389e-05, Final residual = 1.88791e-06, No Iterations 9
suppressing ddt(voidfraction)
DICPCG: Solving for p, Initial residual = 1.98101e-06, Final residual = 8.92532e-07, No Iterations 3
time step continuity errors : sum local = 0.000964672, global = 1.22823e-07, cumulative = -0.159293
suppressing ddt(voidfraction)
DICPCG: Solving for p, Initial residual = 1.32861e-06, Final residual = 7.34657e-07, No Iterations 1
suppressing ddt(voidfraction)
DICPCG: Solving for p, Initial residual = 7.33645e-07, Final residual = 7.33645e-07, No Iterations 0
time step continuity errors : sum local = 0.000964642, global = 1.40617e-07, cumulative = -0.159293
suppressing ddt(voidfraction)
DICPCG: Solving for p, Initial residual = 7.5058e-07, Final residual = 7.5058e-07, No Iterations 0
suppressing ddt(voidfraction)
DICPCG: Solving for p, Initial residual = 7.5058e-07, Final residual = 7.5058e-07, No Iterations 0
time step continuity errors : sum local = 0.000964645, global = 1.43236e-07, cumulative = -0.159292
suppressing ddt(voidfraction)
DICPCG: Solving for p, Initial residual = 7.52815e-07, Final residual = 7.52815e-07, No Iterations 0
suppressing ddt(voidfraction)
DICPCG: Solving for p, Initial residual = 7.52815e-07, Final residual = 7.52815e-07, No Iterations 0
time step continuity errors : sum local = 0.000964646, global = 1.43613e-07, cumulative = -0.159292
ExecutionTime = 1.34 s ClockTime = 6 s

Starting time loop

Time = 0.0055

Courant Number mean: 0.0341785 max: 0.0674852

timeStepFraction() = 0.75
update Ksl.internalField()
DILUPBiCG: Solving for Ux, Initial residual = 0.878075, Final residual = 5.53829e-07, No Iterations 2
DILUPBiCG: Solving for Uy, Initial residual = 0.886317, Final residual = 5.70215e-07, No Iterations 2
DILUPBiCG: Solving for Uz, Initial residual = 0.917164, Final residual = 9.6365e-07, No Iterations 2
suppressing ddt(voidfraction)
DICPCG: Solving for p, Initial residual = 0.00148635, Final residual = 0.000130099, No Iterations 10
suppressing ddt(voidfraction)
DICPCG: Solving for p, Initial residual = 0.000138597, Final residual = 1.09383e-05, No Iterations 10
time step continuity errors : sum local = 0.000966758, global = -1.194e-06, cumulative = -0.159293
suppressing ddt(voidfraction)
DICPCG: Solving for p, Initial residual = 0.000512476, Final residual = 4.46237e-05, No Iterations 15
suppressing ddt(voidfraction)
DICPCG: Solving for p, Initial residual = 4.5975e-05, Final residual = 4.06425e-06, No Iterations 20
time step continuity errors : sum local = 0.000965394, global = 1.05576e-07, cumulative = -0.159293
suppressing ddt(voidfraction)
DICPCG: Solving for p, Initial residual = 1.91072e-05, Final residual = 1.62248e-06, No Iterations 9
suppressing ddt(voidfraction)
DICPCG: Solving for p, Initial residual = 1.72103e-06, Final residual = 9.7355e-07, No Iterations 1
time step continuity errors : sum local = 0.000964705, global = 1.09346e-07, cumulative = -0.159293
suppressing ddt(voidfraction)
DICPCG: Solving for p, Initial residual = 1.46251e-06, Final residual = 9.84669e-07, No Iterations 1
suppressing ddt(voidfraction)
DICPCG: Solving for p, Initial residual = 9.86174e-07, Final residual = 9.86174e-07, No Iterations 0
time step continuity errors : sum local = 0.000964711, global = 1.12688e-07, cumulative = -0.159293
suppressing ddt(voidfraction)
DICPCG: Solving for p, Initial residual = 9.92469e-07, Final residual = 9.92469e-07, No Iterations 0
suppressing ddt(voidfraction)
DICPCG: Solving for p, Initial residual = 9.92469e-07, Final residual = 9.92469e-07, No Iterations 0
time step continuity errors : sum local = 0.000964712, global = 1.13443e-07, cumulative = -0.159293
ExecutionTime = 1.39 s ClockTime = 7 s

Starting time loop

Time = 0.006

Courant Number mean: 0.0348439 max: 0.0695409
Starting up LIGGGHTS
Executing command: 'run 200 pre no post no'
Step Atoms KinEng 1 Volume centerOf dragtota dragtota dragtota
10401 10000 1.371544e-05 3.3678412e-07 5.076e-05 0.0076336181 -0.0054974844 0.00012212879 0.7773872
10410 10000 1.6812824e-05 3.5539783e-07 5.076e-05 0.007638172 -0.0054974844 0.00012212879 0.7773872
10420 10000 2.0815632e-05 3.7825318e-07 5.076e-05 0.0076438286 -0.0054974844 0.00012212879 0.7773872
10430 10000 2.5279728e-05 4.016518e-07 5.076e-05 0.0076501133 -0.0054974844 0.00012212879 0.7773872
10440 10000 3.0204316e-05 4.260652e-07 5.076e-05 0.0076570266 -0.0054974844 0.00012212879 0.7773872
10450 10000 3.5583287e-05 4.4952769e-07 5.076e-05 0.0076645683 -0.0054974844 0.00012212879 0.7773872
10460 10000 4.1403788e-05 4.6798961e-07 5.076e-05 0.0076727384 -0.0054974844 0.00012212879 0.7773872
10470 10000 4.7665931e-05 4.8840897e-07 5.076e-05 0.0076815366 -0.0054974844 0.00012212879 0.7773872
10480 10000 5.4363941e-05 5.1418936e-07 5.076e-05 0.0076909621 -0.0054974844 0.00012212879 0.7773872
10490 10000 6.1475866e-05 5.4061475e-07 5.076e-05 0.0077010135 -0.0054974844 0.00012212879 0.7773872
CFD Coupling established at step 10500
10500 10000 6.9012057e-05 5.746662e-07 5.076e-05 0.0077116889 -0.0054974844 0.00012212879 0.7773872
10510 10000 7.6965257e-05 6.1442407e-07 5.076e-05 0.0077229866 -0.0054974844 0.00012212879 0.7773872
10520 10000 8.5327543e-05 6.6297443e-07 5.076e-05 0.0077349049 -0.0054974844 0.00012212879 0.7773872
10530 10000 9.4088216e-05 7.1545442e-07 5.076e-05 0.0077474417 -0.0054974844 0.00012212879 0.7773872
10540 10000 0.00010324438 7.7705855e-07 5.076e-05 0.007760595 -0.0054974844 0.00012212879 0.7773872
[3] #0 Foam::error::printStack(Foam::Ostream&) in "/opt/openfoam230/platforms/linuxGccDPOpt/lib/libOpenFOAM.so"
[3] #1 Foam::sigFpe::sigHandler(int) in "/opt/openfoam230/platforms/linuxGccDPOpt/lib/libOpenFOAM.so"
[3] #2 Uninterpreted:
[3] #3 LIGGGHTS::ContactModels::ContactModel >::collision(LIGGGHTS::ContactModels::CollisionData&, LIGGGHTS::ContactModels::ForceData&, LIGGGHTS::ContactModels::ForceData&) in "/home/poseidon/OpenFOAM/poseidon-2.3.0/platforms/linuxGccDPOpt/lib/liblagrangianCFDEM-PUBLIC-2.3.0.so"
[3] #4 LIGGGHTS::Walls::Granular > >::compute_force(LAMMPS_NS::FixWallGran*, LIGGGHTS::ContactModels::CollisionData&, double*) in "/home/poseidon/OpenFOAM/poseidon-2.3.0/platforms/linuxGccDPOpt/lib/liblagrangianCFDEM-PUBLIC-2.3.0.so"
[3] #5 LAMMPS_NS::FixWallGran::post_force_primitive(int) in "/home/poseidon/OpenFOAM/poseidon-2.3.0/platforms/linuxGccDPOpt/lib/liblagrangianCFDEM-PUBLIC-2.3.0.so"
[3] #6 LAMMPS_NS::FixWallGran::post_force_wall(int) in "/home/poseidon/OpenFOAM/poseidon-2.3.0/platforms/linuxGccDPOpt/lib/liblagrangianCFDEM-PUBLIC-2.3.0.so"
[3] #7 LAMMPS_NS::FixWallGran::post_force(int) in "/home/poseidon/OpenFOAM/poseidon-2.3.0/platforms/linuxGccDPOpt/lib/liblagrangianCFDEM-PUBLIC-2.3.0.so"
[3] #8 LAMMPS_NS::Modify::post_force(int) in "/home/poseidon/OpenFOAM/poseidon-2.3.0/platforms/linuxGccDPOpt/lib/liblagrangianCFDEM-PUBLIC-2.3.0.so"
[3] #9 LAMMPS_NS::Verlet::run(int) in "/home/poseidon/OpenFOAM/poseidon-2.3.0/platforms/linuxGccDPOpt/lib/liblagrangianCFDEM-PUBLIC-2.3.0.so"
[3] #10 LAMMPS_NS::Run::command(int, char**) in "/home/poseidon/OpenFOAM/poseidon-2.3.0/platforms/linuxGccDPOpt/lib/liblagrangianCFDEM-PUBLIC-2.3.0.so"
[3] #11 void LAMMPS_NS::Input::command_creator(LAMMPS_NS::LAMMPS*, int, char**) in "/home/poseidon/OpenFOAM/poseidon-2.3.0/platforms/linuxGccDPOpt/lib/liblagrangianCFDEM-PUBLIC-2.3.0.so"
[3] #12 LAMMPS_NS::Input::execute_command() in "/home/poseidon/OpenFOAM/poseidon-2.3.0/platforms/linuxGccDPOpt/lib/liblagrangianCFDEM-PUBLIC-2.3.0.so"
[3] #13 LAMMPS_NS::Input::one(char const*) in "/home/poseidon/OpenFOAM/poseidon-2.3.0/platforms/linuxGccDPOpt/lib/liblagrangianCFDEM-PUBLIC-2.3.0.so"
[3] #14 Foam::twoWayMPI::couple() const in "/home/poseidon/OpenFOAM/poseidon-2.3.0/platforms/linuxGccDPOpt/lib/liblagrangianCFDEM-PUBLIC-2.3.0.so"
[3] #15 Foam::cfdemCloud::evolve(Foam::GeometricField&, Foam::GeometricField, Foam::fvPatchField, Foam::volMesh>&, Foam::GeometricField, Foam::fvPatchField, Foam::volMesh>&) in "/home/poseidon/OpenFOAM/poseidon-2.3.0/platforms/linuxGccDPOpt/lib/liblagrangianCFDEM-PUBLIC-2.3.0.so"
[3] #16
[3] in "/home/poseidon/OpenFOAM/poseidon-2.3.0/platforms/linuxGccDPOpt/bin/cfdemSolverPiso"
[3] #17 __libc_start_main in "/lib/i386-linux-gnu/libc.so.6"
[3] #18
[3] in "/home/poseidon/OpenFOAM/poseidon-2.3.0/platforms/linuxGccDPOpt/bin/cfdemSolverPiso"
[poseidon:30184] *** Process received signal ***
[poseidon:30184] Signal: Floating point exception (8)
[poseidon:30184] Signal code: (-6)
[poseidon:30184] Failing at address: 0x75e8
[poseidon:30184] [ 0] [0xb77cf40c]
[poseidon:30184] [ 1] [0xb77cf424]
[poseidon:30184] [ 2] /lib/i386-linux-gnu/libc.so.6(gsignal+0x4f) [0xb4bd41df]
[poseidon:30184] [ 3] /opt/openfoam230/platforms/linuxGccDPOpt/lib/libOpenFOAM.so(_ZN4Foam6sigFpe10sigHandlerEi+0x60) [0xb541c230]
[poseidon:30184] [ 4] [0xb77cf400]
[poseidon:30184] [ 5] /home/poseidon/OpenFOAM/poseidon-2.3.0/platforms/linuxGccDPOpt/lib/liblagrangianCFDEM-PUBLIC-2.3.0.so(_ZN8LIGGGHTS13ContactModels12ContactModelINS0_9GranStyleILi3ELi1ELi0ELi0ELi0EEEE9collisionERNS0_13CollisionDataERNS0_9ForceDataES8_+0xf1) [0xb5cca9c1]
[poseidon:30184] [ 6] /home/poseidon/OpenFOAM/poseidon-2.3.0/platforms/linuxGccDPOpt/lib/liblagrangianCFDEM-PUBLIC-2.3.0.so(_ZN8LIGGGHTS5Walls8GranularINS_13ContactModels12ContactModelINS2_9GranStyleILi3ELi1ELi0ELi0ELi0EEEEEE13compute_forceEPN9LAMMPS_NS11FixWallGranERNS2_13CollisionDataEPd+0x183) [0xb5ccb1a3]
[poseidon:30184] [ 7] /home/poseidon/OpenFOAM/poseidon-2.3.0/platforms/linuxGccDPOpt/lib/liblagrangianCFDEM-PUBLIC-2.3.0.so(_ZN9LAMMPS_NS11FixWallGran20post_force_primitiveEi+0x406) [0xb5b02a86]
[poseidon:30184] [ 8] /home/poseidon/OpenFOAM/poseidon-2.3.0/platforms/linuxGccDPOpt/lib/liblagrangianCFDEM-PUBLIC-2.3.0.so(_ZN9LAMMPS_NS11FixWallGran15post_force_wallEi+0x218) [0xb5b00ee8]
[poseidon:30184] [ 9] /home/poseidon/OpenFOAM/poseidon-2.3.0/platforms/linuxGccDPOpt/lib/liblagrangianCFDEM-PUBLIC-2.3.0.so(_ZN9LAMMPS_NS11FixWallGran10post_forceEi+0x51) [0xb5b00c71]
[poseidon:30184] [10] /home/poseidon/OpenFOAM/poseidon-2.3.0/platforms/linuxGccDPOpt/lib/liblagrangianCFDEM-PUBLIC-2.3.0.so(_ZN9LAMMPS_NS6Modify10post_forceEi+0x5a) [0xb5d25a7a]
[poseidon:30184] [11] /home/poseidon/OpenFOAM/poseidon-2.3.0/platforms/linuxGccDPOpt/lib/liblagrangianCFDEM-PUBLIC-2.3.0.so(_ZN9LAMMPS_NS6Verlet3runEi+0x3e4) [0xb5e6ab54]
[poseidon:30184] [12] /home/poseidon/OpenFOAM/poseidon-2.3.0/platforms/linuxGccDPOpt/lib/liblagrangianCFDEM-PUBLIC-2.3.0.so(_ZN9LAMMPS_NS3Run7commandEiPPc+0x3c6) [0xb5e17fe6]
[poseidon:30184] [13] /home/poseidon/OpenFOAM/poseidon-2.3.0/platforms/linuxGccDPOpt/lib/liblagrangianCFDEM-PUBLIC-2.3.0.so(_ZN9LAMMPS_NS5Input15command_creatorINS_3RunEEEvPNS_6LAMMPSEiPPc+0x4b) [0xb5cef5ab]
[poseidon:30184] [14] /home/poseidon/OpenFOAM/poseidon-2.3.0/platforms/linuxGccDPOpt/lib/liblagrangianCFDEM-PUBLIC-2.3.0.so(_ZN9LAMMPS_NS5Input15execute_commandEv+0x8b6) [0xb5ce9de6]
[poseidon:30184] [15] /home/poseidon/OpenFOAM/poseidon-2.3.0/platforms/linuxGccDPOpt/lib/liblagrangianCFDEM-PUBLIC-2.3.0.so(_ZN9LAMMPS_NS5Input3oneEPKc+0x98) [0xb5cea498]
[poseidon:30184] [16] /home/poseidon/OpenFOAM/poseidon-2.3.0/platforms/linuxGccDPOpt/lib/liblagrangianCFDEM-PUBLIC-2.3.0.so(_ZNK4Foam9twoWayMPI6coupleEv+0xa26) [0xb58682d6]
[poseidon:30184] [17] /home/poseidon/OpenFOAM/poseidon-2.3.0/platforms/linuxGccDPOpt/lib/liblagrangianCFDEM-PUBLIC-2.3.0.so(_ZN4Foam10cfdemCloud6evolveERNS_14GeometricFieldIdNS_12fvPatchFieldENS_7volMeshEEERNS1_INS_6VectorIdEES2_S3_EES9_+0x68) [0xb579f328]
[poseidon:30184] [18] cfdemSolverPiso() [0x80606ef]
[poseidon:30184] [19] /lib/i386-linux-gnu/libc.so.6(__libc_start_main+0xf3) [0xb4bbf4d3]
[poseidon:30184] [20] cfdemSolverPiso() [0x8066511]
[poseidon:30184] *** End of error message ***
--------------------------------------------------------------------------
mpirun noticed that process rank 3 with PID 30184 on node poseidon exited on signal 8 (Floating point exception).
--------------------------------------------------------------------------

j-kerbl's picture

j-kerbl | Tue, 09/09/2014 - 14:15

Hi,

it is hard to say what really happened in your case.
I see that you changed some timing settings, but I'm not sure if that is a problem here.

My suggestion on how to make a nice video to show would be the following modifications to the ErgunTestMPI:
CFD/system/controlDict: endTime 2.0
CFD/0/U:
type uniformFixedValue;
uniformValue table
(
(0.000 (0 0 0.002))
(0.500 (0 0 0.100))
);

This should run without a problem with the latest PUBLIC version.

If it doesn't, it would be nice if I could have your case files to see if I can reproduce the error.

Cheers
Josef