************************************ using userdefined CFDEM_SRC_DIR = /home/project/XY/CFDEM/CFDEMcoupling-PUBLIC-5.x/src using userdefined CFDEM_SOLVER_DIR = /home/project/XY/CFDEM/CFDEMcoupling-PUBLIC-5.x/applications/solvers using userdefined CFDEM_DOC_DIR = /home/project/XY/CFDEM/CFDEMcoupling-PUBLIC-5.x/doc using userdefined CFDEM_UT_DIR = /home/project/XY/CFDEM/CFDEMcoupling-PUBLIC-5.x/applications/utilities using userdefined CFDEM_TUT_DIR = /home/project/XY/CFDEM/CFDEMcoupling-PUBLIC-5.x/tutorials using default CFDEM_LAMMPS_LIB_DIR = /home/project/XY/LIGGGHTS/LIGGGHTS-PUBLIC/src/../lib using default CFDEM_LIGGGHTS_LIB_PATH = /home/project/XY/LIGGGHTS/LIGGGHTS-PUBLIC/src using default CFDEM_ADD_LIBS_DIR = /home/project/XY/CFDEM/CFDEMcoupling-PUBLIC-5.x/src/lagrangian/cfdemParticle/etc/addLibs_universal using default CFDEM_ADD_LIBS_NAME = additionalLibs_5.x using default CFDEM_LIB_DIR = /home/project/XY/CFDEM/CFDEMcoupling-PUBLIC-5.x/platforms/linux64GccDPInt32Opt/lib using default CFDEM_USER_LIB_DIR = /home/project/XY/CFDEM/XY-PUBLIC-5.x/platforms/linux64GccDPInt32Opt/lib using default CFDEM_APP_DIR = /home/project/XY/CFDEM/CFDEMcoupling-PUBLIC-5.x/platforms/linux64GccDPInt32Opt/bin using default CFDEM_USER_APP_DIR = /home/project/XY/CFDEM/XY-PUBLIC-5.x/platforms/linux64GccDPInt32Opt/bin ************************************ mesh needs to be built /*---------------------------------------------------------------------------*\ | ========= | | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | | \\ / O peration | Version: 5.x | | \\ / A nd | Web: www.OpenFOAM.org | | \\/ M anipulation | | \*---------------------------------------------------------------------------*/ Build : 5.x-7f7d351b741b Exec : blockMesh -case CFD Date : Jan 04 2021 Time : 21:27:33 Host : "fh2n0480.localdomain" PID : 4085355 I/O : uncollated Case : ./CFD nProcs : 1 sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE). fileModificationChecking : Monitoring run-time modified files using timeStampMaster (fileModificationSkew 10) allowSystemOperations : Allowing user-supplied system call operations // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // Create time Deleting polyMesh directory "./CFD/constant/polyMesh" Creating block mesh from "./CFD/system/blockMeshDict" Creating block edges No non-planar block faces defined Creating topology blocks Creating topology patches Reading patches section Creating block mesh topology Reading physicalType from existing boundary file Default patch type set to empty Check topology Basic statistics Number of internal faces : 0 Number of boundary faces : 6 Number of defined boundary faces : 6 Number of undefined boundary faces : 0 Checking patch -> block consistency Creating block offsets Creating merge list . Creating polyMesh from blockMesh Creating patches Creating cells Creating points with scale 1 Block 0 cell size : i : 0.1 .. 0.1 j : 0.1 .. 0.1 k : 0.1 .. 0.1 Writing polyMesh ---------------- Mesh Information ---------------- boundingBox: (0 0 0) (1 1 4) nPoints: 4961 nCells: 4000 nFaces: 12900 nInternalFaces: 11100 ---------------- Patches ---------------- patch 0 (start: 11100 size: 100) name: inlet patch 1 (start: 11200 size: 1600) name: wall patch 2 (start: 12800 size: 100) name: outlet End starting DEM run in parallel... /home/project/XY/CFDEM/XY-PUBLIC-5.x/run/twoSpheresGlowinskiMPI_Try/DEM // run_liggghts_DEM // /home/project/XY/CFDEM/XY-PUBLIC-5.x/run/twoSpheresGlowinskiMPI_Try/DEM [fh2n0480.localdomain:4085367] MCW rank 0 bound to socket 0[core 0[hwt 0-1]]: [BB/../../../../../../../../..][../../../../../../../../../..] [fh2n0480.localdomain:4085367] MCW rank 1 bound to socket 0[core 1[hwt 0-1]]: [../BB/../../../../../../../..][../../../../../../../../../..] [fh2n0480.localdomain:4085367] MCW rank 2 bound to socket 0[core 2[hwt 0-1]]: [../../BB/../../../../../../..][../../../../../../../../../..] [fh2n0480.localdomain:4085367] MCW rank 3 bound to socket 0[core 3[hwt 0-1]]: [../../../BB/../../../../../..][../../../../../../../../../..] [fh2n0511.localdomain:2379442] MCW rank 4 bound to socket 0[core 0[hwt 0-1]]: [BB/../../../../../../../../..][../../../../../../../../../..] [fh2n0511.localdomain:2379442] MCW rank 5 bound to socket 0[core 1[hwt 0-1]]: [../BB/../../../../../../../..][../../../../../../../../../..] [fh2n0511.localdomain:2379442] MCW rank 6 bound to socket 0[core 2[hwt 0-1]]: [../../BB/../../../../../../..][../../../../../../../../../..] [fh2n0511.localdomain:2379442] MCW rank 7 bound to socket 0[core 3[hwt 0-1]]: [../../../BB/../../../../../..][../../../../../../../../../..] [fh2n0520.localdomain:2925702] MCW rank 12 bound to socket 0[core 0[hwt 0-1]]: [BB/../../../../../../../../..][../../../../../../../../../..] [fh2n0520.localdomain:2925702] MCW rank 13 bound to socket 0[core 1[hwt 0-1]]: [../BB/../../../../../../../..][../../../../../../../../../..] [fh2n0520.localdomain:2925702] MCW rank 14 bound to socket 0[core 2[hwt 0-1]]: [../../BB/../../../../../../..][../../../../../../../../../..] [fh2n0520.localdomain:2925702] MCW rank 15 bound to socket 0[core 3[hwt 0-1]]: [../../../BB/../../../../../..][../../../../../../../../../..] [fh2n0519.localdomain:2727304] MCW rank 8 bound to socket 0[core 0[hwt 0-1]]: [BB/../../../../../../../../..][../../../../../../../../../..] [fh2n0519.localdomain:2727304] MCW rank 9 bound to socket 0[core 1[hwt 0-1]]: [../BB/../../../../../../../..][../../../../../../../../../..] [fh2n0519.localdomain:2727304] MCW rank 10 bound to socket 0[core 2[hwt 0-1]]: [../../BB/../../../../../../..][../../../../../../../../../..] [fh2n0519.localdomain:2727304] MCW rank 11 bound to socket 0[core 3[hwt 0-1]]: [../../../BB/../../../../../..][../../../../../../../../../..] LIGGGHTS (Version LIGGGHTS-PUBLIC 3.8.0, compiled 2020-12-30-10:53:31 by XY, git commit 6e87343906b56807d7e6e102f94210a8f69e2fe8) log ../DEM/log.liggghts thermo_log ../DEM/post/thermo.txt atom_style granular atom_modify map array communicate single vel yes boundary f f f newton off units si processors 4 2 2 region reg block 0. 1. 0. 1. 0. 4. units box create_box 1 reg Created orthogonal box = (0 0 0) to (1 1 4) 4 by 2 by 2 MPI processor grid neighbor 0.3 bin neigh_modify delay 0 binsize 0.01 # Material properties required for new pair styles fix m1 all property/global youngsModulus peratomtype 5.e7 fix m2 all property/global poissonsRatio peratomtype 0.45 fix m3 all property/global coefficientRestitution peratomtypepair 1 0.9 fix m4 all property/global coefficientFriction peratomtypepair 1 0.5 # pair style pair_style gran model hertz tangential history #Hertzian without cohesion pair_coeff * * # timestep, gravity timestep 0.00003 fix gravi all gravity 981 vector 0.0 0.0 -1.0 # walls fix xwalls1 all wall/gran model hertz tangential history primitive type 1 xplane 0. fix xwalls2 all wall/gran model hertz tangential history primitive type 1 xplane 1. fix ywalls1 all wall/gran model hertz tangential history primitive type 1 yplane 0. fix ywalls2 all wall/gran model hertz tangential history primitive type 1 yplane 1. fix zwalls1 all wall/gran model hertz tangential history primitive type 1 zplane 0. fix zwalls2 all wall/gran model hertz tangential history primitive type 1 zplane 4. # cfd coupling fix cfd all couple/cfd couple_every 10 mpi couple_every as specified in LIGGGHTS is overriden by calling external program (../fix_cfd_coupling.cpp:114) fix cfd2 all couple/cfd/force # create single partciles create_atoms 1 single .5 .5 3.5 units box Created 1 atoms create_atoms 1 single .5 .5 3.16 units box Created 1 atoms set atom 1 diameter 0.167 density 1.5 vx 0 vy 0 vz 0 Setting atom values ... 1 settings made for diameter 1 settings made for density 1 settings made for vx 1 settings made for vy 1 settings made for vz set atom 2 diameter 0.167 density 1.5 vx 0 vy 0 vz 0 Setting atom values ... 1 settings made for diameter 1 settings made for density 1 settings made for vx 1 settings made for vy 1 settings made for vz variable vx1 equal vx[1] variable vy1 equal vy[1] variable vz1 equal vz[1] variable vx2 equal vx[2] variable vy2 equal vy[2] variable vz2 equal vz[2] variable x1 equal x[1] variable y1 equal y[1] variable z1 equal z[1] variable x2 equal x[2] variable y2 equal y[2] variable z2 equal z[2] variable time equal step*dt fix extra1 all print 100 "${time} ${vx1} ${vy1} ${vz1}" file ../DEM/post/velocity_particle_1.txt title "%" screen no fix extra2 all print 100 "${time} ${vx2} ${vy2} ${vz2}" file ../DEM/post/velocity_particle_2.txt title "%" screen no fix extra3 all print 100 "${time} ${x1} ${y1} ${z1}" file ../DEM/post/position_particle_1.txt title "%" screen no fix extra4 all print 100 "${time} ${x2} ${y2} ${z2}" file ../DEM/post/position_particle_2.txt title "%" screen no # apply nve integration to all particles that are inserted as single particles fix integr all nve/sphere #wenn das ausgeblendet, dann kein vel update # screen output compute rke all erotate/sphere thermo_style custom step atoms ke c_rke vol thermo 1000 thermo_modify lost ignore norm no compute_modify thermo_temp dynamic yes WARNING: Changing thermo_temp compute object. This object is deprecated and will be removed in the future. (../compute.cpp:161) # insert the first particles so that dump is not empty dump dmp all custom 100 ../DEM/post/dump.liggghts_run id type x y z ix iy iz vx vy vz fx fy fz omegax omegay omegaz radius #force : f_couple_cfd[0] f_couple_cfd[1] f_couple_cfd[2] #node : f_couple_cfd[6] #cell id : f_couple_cfd[7] run 1 Setting up run at Mon Jan 4 21:27:34 2021 Memory usage per processor = 27.4725 Mbytes Step Atoms KinEng rke Volume 0 2 0 0 4 1 2 3.1682527e-06 0 4 Loop time of 6.88393e-05 on 16 procs for 1 steps with 2 atoms, finish time Mon Jan 4 21:27:35 2021 Pair time (%) = 1.37673e-06 (1.99991) Neigh time (%) = 0 (0) Comm time (%) = 1.39623e-05 (20.2824) Outpt time (%) = 2.94335e-05 (42.7567) Other time (%) = 2.40669e-05 (34.9609) Nlocal: 0.125 ave 2 max 0 min Histogram: 15 0 0 0 0 0 0 0 0 1 Nghost: 0.875 ave 2 max 0 min Histogram: 9 0 0 0 0 0 0 0 0 7 Neighs: 0.0625 ave 1 max 0 min Histogram: 15 0 0 0 0 0 0 0 0 1 Total # of neighbors = 1 Ave neighs/atom = 0.5 Neighbor list builds = 0 Dangerous builds = 0 starting CFDDEM run in parallel... /home/project/XY/CFDEM/XY-PUBLIC-5.x/run/twoSpheresGlowinskiMPI_Try/CFD Decomposing case. /*---------------------------------------------------------------------------*\ | ========= | | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | | \\ / O peration | Version: 5.x | | \\ / A nd | Web: www.OpenFOAM.org | | \\/ M anipulation | | \*---------------------------------------------------------------------------*/ Build : 5.x-7f7d351b741b Exec : decomposePar -force Date : Jan 04 2021 Time : 21:27:35 Host : "fh2n0480.localdomain" PID : 4085450 I/O : uncollated Case : /pfs/data4/project/project/XY/CFDEM/XY-PUBLIC-5.x/run/twoSpheresGlowinskiMPI_Try/CFD nProcs : 1 sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE). fileModificationChecking : Monitoring run-time modified files using timeStampMaster (fileModificationSkew 10) allowSystemOperations : Allowing user-supplied system call operations // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // Create time Decomposing mesh region0 Create mesh Calculating distribution of cells Selecting decompositionMethod simple Finished decomposition in 0 s Calculating original mesh data Distributing cells to processors Distributing faces to processors Distributing points to processors Constructing processor meshes Processor 0 Number of cells = 290 Number of faces shared with processor 1 = 108 Number of faces shared with processor 4 = 60 Number of faces shared with processor 8 = 10 Number of faces shared with processor 9 = 1 Number of processor patches = 4 Number of processor faces = 179 Number of boundary faces = 171 Processor 1 Number of cells = 210 Number of faces shared with processor 0 = 108 Number of faces shared with processor 2 = 100 Number of faces shared with processor 5 = 40 Number of faces shared with processor 9 = 14 Number of processor patches = 4 Number of processor faces = 262 Number of boundary faces = 54 Processor 2 Number of cells = 290 Number of faces shared with processor 1 = 100 Number of faces shared with processor 3 = 108 Number of faces shared with processor 6 = 60 Number of faces shared with processor 10 = 10 Number of faces shared with processor 11 = 1 Number of processor patches = 5 Number of processor faces = 279 Number of boundary faces = 71 Processor 3 Number of cells = 210 Number of faces shared with processor 2 = 108 Number of faces shared with processor 7 = 40 Number of faces shared with processor 11 = 14 Number of processor patches = 3 Number of processor faces = 162 Number of boundary faces = 154 Processor 4 Number of cells = 300 Number of faces shared with processor 0 = 60 Number of faces shared with processor 5 = 100 Number of faces shared with processor 12 = 14 Number of faces shared with processor 13 = 1 Number of processor patches = 4 Number of processor faces = 175 Number of boundary faces = 175 Processor 5 Number of cells = 200 Number of faces shared with processor 1 = 40 Number of faces shared with processor 4 = 100 Number of faces shared with processor 6 = 100 Number of faces shared with processor 13 = 10 Number of processor patches = 4 Number of processor faces = 250 Number of boundary faces = 50 Processor 6 Number of cells = 300 Number of faces shared with processor 2 = 60 Number of faces shared with processor 5 = 100 Number of faces shared with processor 7 = 100 Number of faces shared with processor 14 = 14 Number of faces shared with processor 15 = 1 Number of processor patches = 5 Number of processor faces = 275 Number of boundary faces = 75 Processor 7 Number of cells = 200 Number of faces shared with processor 3 = 40 Number of faces shared with processor 6 = 100 Number of faces shared with processor 15 = 10 Number of processor patches = 3 Number of processor faces = 150 Number of boundary faces = 150 Processor 8 Number of cells = 200 Number of faces shared with processor 0 = 10 Number of faces shared with processor 9 = 100 Number of faces shared with processor 12 = 40 Number of processor patches = 3 Number of processor faces = 150 Number of boundary faces = 150 Processor 9 Number of cells = 300 Number of faces shared with processor 0 = 1 Number of faces shared with processor 1 = 14 Number of faces shared with processor 8 = 100 Number of faces shared with processor 10 = 100 Number of faces shared with processor 13 = 60 Number of processor patches = 5 Number of processor faces = 275 Number of boundary faces = 75 Processor 10 Number of cells = 200 Number of faces shared with processor 2 = 10 Number of faces shared with processor 9 = 100 Number of faces shared with processor 11 = 100 Number of faces shared with processor 14 = 40 Number of processor patches = 4 Number of processor faces = 250 Number of boundary faces = 50 Processor 11 Number of cells = 300 Number of faces shared with processor 2 = 1 Number of faces shared with processor 3 = 14 Number of faces shared with processor 10 = 100 Number of faces shared with processor 15 = 60 Number of processor patches = 4 Number of processor faces = 175 Number of boundary faces = 175 Processor 12 Number of cells = 210 Number of faces shared with processor 4 = 14 Number of faces shared with processor 8 = 40 Number of faces shared with processor 13 = 108 Number of processor patches = 3 Number of processor faces = 162 Number of boundary faces = 154 Processor 13 Number of cells = 290 Number of faces shared with processor 4 = 1 Number of faces shared with processor 5 = 10 Number of faces shared with processor 9 = 60 Number of faces shared with processor 12 = 108 Number of faces shared with processor 14 = 100 Number of processor patches = 5 Number of processor faces = 279 Number of boundary faces = 71 Processor 14 Number of cells = 210 Number of faces shared with processor 6 = 14 Number of faces shared with processor 10 = 40 Number of faces shared with processor 13 = 100 Number of faces shared with processor 15 = 108 Number of processor patches = 4 Number of processor faces = 262 Number of boundary faces = 54 Processor 15 Number of cells = 290 Number of faces shared with processor 6 = 1 Number of faces shared with processor 7 = 10 Number of faces shared with processor 11 = 60 Number of faces shared with processor 14 = 108 Number of processor patches = 4 Number of processor faces = 179 Number of boundary faces = 171 Number of processor faces = 1732 Max number of cells = 300 (20% above average 250) Max number of processor patches = 5 (25% above average 4) Max number of faces between processors = 279 (28.8684% above average 216.5) Time = 0 Processor 0: field transfer Processor 1: field transfer Processor 2: field transfer Processor 3: field transfer Processor 4: field transfer Processor 5: field transfer Processor 6: field transfer Processor 7: field transfer Processor 8: field transfer Processor 9: field transfer Processor 10: field transfer Processor 11: field transfer Processor 12: field transfer Processor 13: field transfer Processor 14: field transfer Processor 15: field transfer End do nothing. // run_parallel_cfdemSolverIB_twoSpheresGlowinskiMPI_CFDDEM // /home/project/XY/CFDEM/XY-PUBLIC-5.x/run/twoSpheresGlowinskiMPI_Try/CFD [fh2n0480.localdomain:4085495] MCW rank 0 bound to socket 0[core 0[hwt 0-1]]: [BB/../../../../../../../../..][../../../../../../../../../..] [fh2n0480.localdomain:4085495] MCW rank 1 bound to socket 0[core 1[hwt 0-1]]: [../BB/../../../../../../../..][../../../../../../../../../..] [fh2n0480.localdomain:4085495] MCW rank 2 bound to socket 0[core 2[hwt 0-1]]: [../../BB/../../../../../../..][../../../../../../../../../..] [fh2n0480.localdomain:4085495] MCW rank 3 bound to socket 0[core 3[hwt 0-1]]: [../../../BB/../../../../../..][../../../../../../../../../..] [fh2n0511.localdomain:2379525] MCW rank 4 bound to socket 0[core 0[hwt 0-1]]: [BB/../../../../../../../../..][../../../../../../../../../..] [fh2n0511.localdomain:2379525] MCW rank 5 bound to socket 0[core 1[hwt 0-1]]: [../BB/../../../../../../../..][../../../../../../../../../..] [fh2n0511.localdomain:2379525] MCW rank 6 bound to socket 0[core 2[hwt 0-1]]: [../../BB/../../../../../../..][../../../../../../../../../..] [fh2n0511.localdomain:2379525] MCW rank 7 bound to socket 0[core 3[hwt 0-1]]: [../../../BB/../../../../../..][../../../../../../../../../..] [fh2n0520.localdomain:2925781] MCW rank 12 bound to socket 0[core 0[hwt 0-1]]: [BB/../../../../../../../../..][../../../../../../../../../..] [fh2n0520.localdomain:2925781] MCW rank 13 bound to socket 0[core 1[hwt 0-1]]: [../BB/../../../../../../../..][../../../../../../../../../..] [fh2n0520.localdomain:2925781] MCW rank 14 bound to socket 0[core 2[hwt 0-1]]: [../../BB/../../../../../../..][../../../../../../../../../..] [fh2n0520.localdomain:2925781] MCW rank 15 bound to socket 0[core 3[hwt 0-1]]: [../../../BB/../../../../../..][../../../../../../../../../..] [fh2n0519.localdomain:2727384] MCW rank 8 bound to socket 0[core 0[hwt 0-1]]: [BB/../../../../../../../../..][../../../../../../../../../..] [fh2n0519.localdomain:2727384] MCW rank 9 bound to socket 0[core 1[hwt 0-1]]: [../BB/../../../../../../../..][../../../../../../../../../..] [fh2n0519.localdomain:2727384] MCW rank 10 bound to socket 0[core 2[hwt 0-1]]: [../../BB/../../../../../../..][../../../../../../../../../..] [fh2n0519.localdomain:2727384] MCW rank 11 bound to socket 0[core 3[hwt 0-1]]: [../../../BB/../../../../../..][../../../../../../../../../..] /*---------------------------------------------------------------------------*\ | ========= | | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | | \\ / O peration | Version: 5.x | | \\ / A nd | Web: www.OpenFOAM.org | | \\/ M anipulation | | \*---------------------------------------------------------------------------*/ Build : 5.x-7f7d351b741b Exec : cfdemSolverIB -parallel Date : Jan 04 2021 Time : 21:27:40 Host : "fh2n0480.localdomain" PID : 4085511 I/O : uncollated Case : /pfs/data4/project/project/XY/CFDEM/XY-PUBLIC-5.x/run/twoSpheresGlowinskiMPI_Try/CFD nProcs : 16 Slaves : 15 ( "fh2n0480.localdomain.4085512" "fh2n0480.localdomain.4085513" "fh2n0480.localdomain.4085514" "fh2n0511.localdomain.2379531" "fh2n0511.localdomain.2379532" "fh2n0511.localdomain.2379533" "fh2n0511.localdomain.2379534" "fh2n0519.localdomain.2727390" "fh2n0519.localdomain.2727391" "fh2n0519.localdomain.2727392" "fh2n0519.localdomain.2727393" "fh2n0520.localdomain.2925787" "fh2n0520.localdomain.2925788" "fh2n0520.localdomain.2925789" "fh2n0520.localdomain.2925790" ) 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 (fileModificationSkew 10) allowSystemOperations : Allowing user-supplied system call operations // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // Create time Create mesh for time = 0 Selecting dynamicFvMesh dynamicRefineFvMesh PISO: Operating solver in PISO mode Reading field p Reading physical velocity field U Note: only if voidfraction at boundary is 1, U is superficial velocity!!! Reading physical velocity field U Note: only if voidfraction at boundary is 1, U is superficial velocity!!! Creating dummy density field rho = 1 Reading field phiIB Reading field phiIB Reading/calculating face flux field phi Selecting incompressible transport model Newtonian Selecting turbulence model type laminar Selecting laminar stress model Stokes Reading field interFace Reading g Selecting locateModel engineIB Selecting dataExchangeModel twoWayMPI Starting up LIGGGHTS for first time execution Executing input script '../DEM/in.liggghts_run' LIGGGHTS (Version LIGGGHTS-PUBLIC 3.8.0, compiled 2020-12-30-11:01:52 by XY, git commit 6e87343906b56807d7e6e102f94210a8f69e2fe8) log ../DEM/log.liggghts thermo_log ../DEM/post/thermo.txt atom_style granular atom_modify map array communicate single vel yes boundary f f f newton off units si processors 4 2 2 region reg block 0. 1. 0. 1. 0. 4. units box create_box 1 reg Created orthogonal box = (0 0 0) to (1 1 4) 4 by 2 by 2 MPI processor grid neighbor 0.3 bin neigh_modify delay 0 binsize 0.01 # Material properties required for new pair styles fix m1 all property/global youngsModulus peratomtype 5.e7 fix m2 all property/global poissonsRatio peratomtype 0.45 fix m3 all property/global coefficientRestitution peratomtypepair 1 0.9 fix m4 all property/global coefficientFriction peratomtypepair 1 0.5 # pair style pair_style gran model hertz tangential history #Hertzian without cohesion pair_coeff * * # timestep, gravity timestep 0.00003 fix gravi all gravity 981 vector 0.0 0.0 -1.0 # walls fix xwalls1 all wall/gran model hertz tangential history primitive type 1 xplane 0. fix xwalls2 all wall/gran model hertz tangential history primitive type 1 xplane 1. fix ywalls1 all wall/gran model hertz tangential history primitive type 1 yplane 0. fix ywalls2 all wall/gran model hertz tangential history primitive type 1 yplane 1. fix zwalls1 all wall/gran model hertz tangential history primitive type 1 zplane 0. fix zwalls2 all wall/gran model hertz tangential history primitive type 1 zplane 4. # cfd coupling fix cfd all couple/cfd couple_every 10 mpi couple_every as specified in LIGGGHTS is overriden by calling external program (../fix_cfd_coupling.cpp:114) fix cfd2 all couple/cfd/force # create single partciles create_atoms 1 single .5 .5 3.5 units box Created 1 atoms create_atoms 1 single .5 .5 3.16 units box Created 1 atoms set atom 1 diameter 0.167 density 1.5 vx 0 vy 0 vz 0 Setting atom values ... 1 settings made for diameter 1 settings made for density 1 settings made for vx 1 settings made for vy 1 settings made for vz set atom 2 diameter 0.167 density 1.5 vx 0 vy 0 vz 0 Setting atom values ... 1 settings made for diameter 1 settings made for density 1 settings made for vx 1 settings made for vy 1 settings made for vz variable vx1 equal vx[1] variable vy1 equal vy[1] variable vz1 equal vz[1] variable vx2 equal vx[2] variable vy2 equal vy[2] variable vz2 equal vz[2] variable x1 equal x[1] variable y1 equal y[1] variable z1 equal z[1] variable x2 equal x[2] variable y2 equal y[2] variable z2 equal z[2] variable time equal step*dt fix extra1 all print 100 "${time} ${vx1} ${vy1} ${vz1}" file ../DEM/post/velocity_particle_1.txt title "%" screen no fix extra2 all print 100 "${time} ${vx2} ${vy2} ${vz2}" file ../DEM/post/velocity_particle_2.txt title "%" screen no fix extra3 all print 100 "${time} ${x1} ${y1} ${z1}" file ../DEM/post/position_particle_1.txt title "%" screen no fix extra4 all print 100 "${time} ${x2} ${y2} ${z2}" file ../DEM/post/position_particle_2.txt title "%" screen no # apply nve integration to all particles that are inserted as single particles fix integr all nve/sphere #wenn das ausgeblendet, dann kein vel update # screen output compute rke all erotate/sphere thermo_style custom step atoms ke c_rke vol thermo 1000 thermo_modify lost ignore norm no compute_modify thermo_temp dynamic yes WARNING: Changing thermo_temp compute object. This object is deprecated and will be removed in the future. (../compute.cpp:161) # insert the first particles so that dump is not empty dump dmp all custom 100 ../DEM/post/dump.liggghts_run id type x y z ix iy iz vx vy vz fx fy fz omegax omegay omegaz radius #force : f_couple_cfd[0] f_couple_cfd[1] f_couple_cfd[2] #node : f_couple_cfd[6] #cell id : f_couple_cfd[7] run 1 Setting up run at Mon Jan 4 21:27:41 2021 Memory usage per processor = 27.4725 Mbytes Step Atoms KinEng rke Volume 0 2 0 0 4 1 2 3.1682527e-06 0 4 Loop time of 0.000152375 on 16 procs for 1 steps with 2 atoms, finish time Mon Jan 4 21:27:41 2021 Pair time (%) = 1.20472e-06 (0.790629) Neigh time (%) = 0 (0) Comm time (%) = 2.67074e-05 (17.5274) Outpt time (%) = 4.47649e-05 (29.378) Other time (%) = 7.96984e-05 (52.304) Nlocal: 0.125 ave 2 max 0 min Histogram: 15 0 0 0 0 0 0 0 0 1 Nghost: 0.875 ave 2 max 0 min Histogram: 9 0 0 0 0 0 0 0 0 7 Neighs: 0.0625 ave 1 max 0 min Histogram: 15 0 0 0 0 0 0 0 0 1 Total # of neighbors = 1 Ave neighs/atom = 0.5 Neighbor list builds = 0 Dangerous builds = 0 Selecting IOModel basicIO Selecting probeModel off Selecting voidFractionModel IB W A R N I N G - do not use in combination with differentialRegion model! Selecting averagingModel dilute Selecting clockModel off start clock measurement at t >0.0006 Selecting smoothingModel off Selecting meshMotionModel noMeshMotion CFDEMcoupling version: cfdem-3.8.1 , compatible to LIGGGHTS version: 3.8.0 , compatible to OF version and build: 5.x-commit-538044ac05c4672b37c7df607dca1116fa88df88 You are currently using: OF version: 5.x OF build: 5.x-7f7d351b741b CFDEM build: 3.8.0-3-gb237264-dirty cg is set to: 1 LIGGGHTS hosts 1 atom types. If BC are important, please provide volScalarFields -imp/expParticleForces- ignoring ddt(voidfraction) Selecting forceModel ShirgaonkarIB Selecting forceSubModel ImEx reading switches: looking for treatForceExplicit ... not found in dict, using default. treatForceExplicit = false looking for verbose ... not found in dict, using default. verbose = false Selecting forceModel ArchimedesIB Selecting forceSubModel ImEx reading switches: looking for treatForceExplicit ... not found in dict, using default. treatForceExplicit = true accounting for Archimedes only on DEM side! Selecting liggghtsCommandModel runLiggghts ,provide dicts, numbered from 0 to n Starting time loop --> FOAM Warning : From function void Foam::timeControl::read(const Foam::dictionary&) in file db/functionObjects/timeControl/timeControl.C at line 89 Reading "/pfs/data4/project/project/XY/CFDEM/XY-PUBLIC-5.x/run/twoSpheresGlowinskiMPI_Try/CFD/system/controlDict.functions.probes" from line 62 to line 80 Using deprecated 'outputControl' Please use 'writeControl' with 'writeInterval' Time = 0.0003 Selected 0 cells for refinement out of 4000. Selected 0 split points out of a possible 0. Courant Number mean: 0 max: 0 - evolve() timeStepFraction() = 2 Starting up LIGGGHTS Executing command: 'run 10 ' run 10 Setting up run at Mon Jan 4 21:27:41 2021 Memory usage per processor = 27.4725 Mbytes Step Atoms KinEng rke Volume 1 2 3.1682527e-06 0 4 CFD Coupling established at step 10 11 2 0.00038335858 0 4 Loop time of 0.000405411 on 16 procs for 10 steps with 2 atoms, finish time Mon Jan 4 21:27:41 2021 Pair time (%) = 2.99083e-06 (0.737728) Neigh time (%) = 0 (0) Comm time (%) = 8.36516e-05 (20.6338) Outpt time (%) = 2.43671e-05 (6.01047) Other time (%) = 0.000294401 (72.618) Nlocal: 0.125 ave 2 max 0 min Histogram: 15 0 0 0 0 0 0 0 0 1 Nghost: 0.875 ave 2 max 0 min Histogram: 9 0 0 0 0 0 0 0 0 7 Neighs: 0.0625 ave 1 max 0 min Histogram: 15 0 0 0 0 0 0 0 0 1 Total # of neighbors = 1 Ave neighs/atom = 0.5 Neighbor list builds = 0 Dangerous builds = 0 LIGGGHTS finished Foam::cfdemCloudIB::reAllocArrays() nr particles = 2 evolve done. DILUPBiCG: Solving for Ux, Initial residual = 1, Final residual = 1.17206e-07, No Iterations 2 DILUPBiCG: Solving for Uy, Initial residual = 1, Final residual = 1.29815e-07, No Iterations 2 DILUPBiCG: Solving for Uz, Initial residual = 1, Final residual = 1.62662e-07, No Iterations 2 DICPCG: Solving for p, Initial residual = 0.612887, Final residual = 9.92843e-07, No Iterations 70 time step continuity errors : sum local = 1.30506e-15, global = -1.9984e-19, cumulative = -1.9984e-19 DICPCG: Solving for p, Initial residual = 0.695098, Final residual = 8.50898e-07, No Iterations 69 time step continuity errors : sum local = 1.29661e-15, global = -3.33067e-20, cumulative = -2.33147e-19 DICPCG: Solving for p, Initial residual = 0.759346, Final residual = 6.20971e-07, No Iterations 69 time step continuity errors : sum local = 1.10768e-15, global = -2.9976e-19, cumulative = -5.32907e-19 DICPCG: Solving for p, Initial residual = 0.713144, Final residual = 9.19272e-07, No Iterations 69 time step continuity errors : sum local = 1.08679e-15, global = -4.32987e-19, cumulative = -9.65894e-19 particleCloud.calcVelocityCorrection() DICPCG: Solving for phiIB, Initial residual = 1, Final residual = 7.08889e-07, No Iterations 56 ExecutionTime = 0.35 s ClockTime = 0 s Time = 0.0006 Selected 124 cells for refinement out of 4000. Refined from 4000 to 4868 cells. Selected 0 split points out of a possible 124. Courant Number mean: 3.13168e-06 max: 0.00069752 - evolve() timeStepFraction() = 2 Starting up LIGGGHTS Executing command: 'run 10 ' run 10 Setting up run at Mon Jan 4 21:27:41 2021 Memory usage per processor = 27.4725 Mbytes Step Atoms KinEng rke Volume 11 2 0.00038335858 0 4 CFD Coupling established at step 20 21 2 0.00083081303 0 4 Loop time of 0.000161851 on 16 procs for 10 steps with 2 atoms, finish time Mon Jan 4 21:27:41 2021 Pair time (%) = 2.79484e-06 (1.7268) Neigh time (%) = 0 (0) Comm time (%) = 2.22154e-05 (13.7258) Outpt time (%) = 2.78402e-05 (17.2011) Other time (%) = 0.000109001 (67.3463) Nlocal: 0.125 ave 2 max 0 min Histogram: 15 0 0 0 0 0 0 0 0 1 Nghost: 0.875 ave 2 max 0 min Histogram: 9 0 0 0 0 0 0 0 0 7 Neighs: 0.0625 ave 1 max 0 min Histogram: 15 0 0 0 0 0 0 0 0 1 Total # of neighbors = 1 Ave neighs/atom = 0.5 Neighbor list builds = 0 Dangerous builds = 0 LIGGGHTS finished nr particles = 2 evolve done. DILUPBiCG: Solving for Ux, Initial residual = 0.334417, Final residual = 4.75447e-08, No Iterations 2 DILUPBiCG: Solving for Uy, Initial residual = 0.334417, Final residual = 4.61869e-08, No Iterations 2 DILUPBiCG: Solving for Uz, Initial residual = 0.164322, Final residual = 6.68457e-08, No Iterations 2 DICPCG: Solving for p, Initial residual = 0.765863, Final residual = 8.49725e-07, No Iterations 73 time step continuity errors : sum local = 4.33393e-12, global = 3.3317e-20, cumulative = -9.32577e-19 DICPCG: Solving for p, Initial residual = 0.0500405, Final residual = 9.06754e-07, No Iterations 63 time step continuity errors : sum local = 2.75013e-12, global = 6.24508e-19, cumulative = -3.08069e-19 DICPCG: Solving for p, Initial residual = 0.00197264, Final residual = 9.17166e-07, No Iterations 36 time step continuity errors : sum local = 2.93093e-12, global = -2.41462e-19, cumulative = -5.49531e-19 DICPCG: Solving for p, Initial residual = 0.000120208, Final residual = 9.12865e-07, No Iterations 11 time step continuity errors : sum local = 2.91541e-12, global = 4.74619e-19, cumulative = -7.49127e-20 particleCloud.calcVelocityCorrection() DICPCG: Solving for phiIB, Initial residual = 0.573966, Final residual = 9.67989e-07, No Iterations 54 ExecutionTime = 0.39 s ClockTime = 0 s Time = 0.0009 Selected 368 cells for refinement out of 4868. Refined from 4868 to 7444 cells. Selected 16 split points out of a possible 396. Unrefined from 7444 to 7332 cells. Courant Number mean: 1.8483e-05 max: 0.00534861 - evolve() timeStepFraction() = 2 Starting up LIGGGHTS Executing command: 'run 10 ' run 10 Setting up run at Mon Jan 4 21:27:41 2021 Memory usage per processor = 27.4725 Mbytes Step Atoms KinEng rke Volume 21 2 0.00083081303 0 4 CFD Coupling established at step 30 31 2 0.001042339 0 4 Loop time of 0.000168143 on 16 procs for 10 steps with 2 atoms, finish time Mon Jan 4 21:27:41 2021 Pair time (%) = 3.08122e-06 (1.8325) Neigh time (%) = 0 (0) Comm time (%) = 2.49866e-05 (14.8604) Outpt time (%) = 2.63098e-05 (15.6473) Other time (%) = 0.000113765 (67.6598) Nlocal: 0.125 ave 2 max 0 min Histogram: 15 0 0 0 0 0 0 0 0 1 Nghost: 0.875 ave 2 max 0 min Histogram: 9 0 0 0 0 0 0 0 0 7 Neighs: 0.0625 ave 1 max 0 min Histogram: 15 0 0 0 0 0 0 0 0 1 Total # of neighbors = 1 Ave neighs/atom = 0.5 Neighbor list builds = 0 Dangerous builds = 0 LIGGGHTS finished nr particles = 2 evolve done. DILUPBiCG: Solving for Ux, Initial residual = 0.275794, Final residual = 2.98659e-07, No Iterations 2 DILUPBiCG: Solving for Uy, Initial residual = 0.276106, Final residual = 2.45299e-07, No Iterations 2 DILUPBiCG: Solving for Uz, Initial residual = 0.0256583, Final residual = 1.40342e-07, No Iterations 2 DICPCG: Solving for p, Initial residual = 0.852825, Final residual = 9.58436e-07, No Iterations 76 time step continuity errors : sum local = 1.73977e-11, global = -7.17564e-19, cumulative = -7.92477e-19 DICPCG: Solving for p, Initial residual = 0.112769, Final residual = 8.96923e-07, No Iterations 70 time step continuity errors : sum local = 5.97634e-12, global = -9.59038e-19, cumulative = -1.75151e-18 DICPCG: Solving for p, Initial residual = 0.00624313, Final residual = 8.02321e-07, No Iterations 64 time step continuity errors : sum local = 5.07118e-12, global = 8.17662e-20, cumulative = -1.66975e-18 DICPCG: Solving for p, Initial residual = 0.000526106, Final residual = 8.90144e-07, No Iterations 34 time step continuity errors : sum local = 5.92509e-12, global = 7.34238e-20, cumulative = -1.59632e-18 particleCloud.calcVelocityCorrection() DICPCG: Solving for phiIB, Initial residual = 0.620332, Final residual = 7.5111e-07, No Iterations 57 ExecutionTime = 0.49 s ClockTime = 0 s Time = 0.0012 Selected 0 cells for refinement out of 7332. Selected 144 split points out of a possible 380. Unrefined from 7332 to 6324 cells. Courant Number mean: 2.17962e-05 max: 0.00487998 - evolve() timeStepFraction() = 2 Starting up LIGGGHTS Executing command: 'run 10 ' run 10 Setting up run at Mon Jan 4 21:27:41 2021 Memory usage per processor = 27.4725 Mbytes Step Atoms KinEng rke Volume 31 2 0.001042339 0 4 CFD Coupling established at step 40 41 2 0.0010440444 0 4 Loop time of 0.000163609 on 16 procs for 10 steps with 2 atoms, finish time Mon Jan 4 21:27:41 2021 Pair time (%) = 2.82115e-06 (1.72433) Neigh time (%) = 0 (0) Comm time (%) = 2.39418e-05 (14.6336) Outpt time (%) = 2.47994e-05 (15.1577) Other time (%) = 0.000112046 (68.4844) Nlocal: 0.125 ave 2 max 0 min Histogram: 15 0 0 0 0 0 0 0 0 1 Nghost: 0.875 ave 2 max 0 min Histogram: 9 0 0 0 0 0 0 0 0 7 Neighs: 0.0625 ave 1 max 0 min Histogram: 15 0 0 0 0 0 0 0 0 1 Total # of neighbors = 1 Ave neighs/atom = 0.5 Neighbor list builds = 0 Dangerous builds = 0 LIGGGHTS finished nr particles = 2 evolve done. DILUPBiCG: Solving for Ux, Initial residual = 0.164201, Final residual = 3.76574e-07, No Iterations 2 DILUPBiCG: Solving for Uy, Initial residual = 0.164243, Final residual = 3.59466e-07, No Iterations 2 DILUPBiCG: Solving for Uz, Initial residual = 0.0128872, Final residual = 8.21797e-08, No Iterations 2 DICPCG: Solving for p, Initial residual = 0.564446, Final residual = 5.87881e-07, No Iterations 77 time step continuity errors : sum local = 4.75593e-12, global = 1.38839e-19, cumulative = -1.45749e-18 DICPCG: Solving for p, Initial residual = 0.12397, Final residual = 8.13664e-07, No Iterations 69 time step continuity errors : sum local = 3.85652e-12, global = 2.47082e-19, cumulative = -1.2104e-18 DICPCG: Solving for p, Initial residual = 0.00592188, Final residual = 7.58782e-07, No Iterations 46 time step continuity errors : sum local = 4.45166e-12, global = -7.02179e-19, cumulative = -1.91258e-18 DICPCG: Solving for p, Initial residual = 0.000504464, Final residual = 9.3165e-07, No Iterations 38 time step continuity errors : sum local = 5.29615e-12, global = -4.93942e-19, cumulative = -2.40652e-18 particleCloud.calcVelocityCorrection() DICPCG: Solving for phiIB, Initial residual = 0.265759, Final residual = 7.46798e-07, No Iterations 59 ExecutionTime = 0.55 s ClockTime = 0 s Time = 0.0015 Selected 0 cells for refinement out of 6324. Selected 8 split points out of a possible 272. Unrefined from 6324 to 6268 cells. Courant Number mean: 2.77333e-05 max: 0.00588432 - evolve() timeStepFraction() = 2 Starting up LIGGGHTS Executing command: 'run 10 ' run 10 Setting up run at Mon Jan 4 21:27:41 2021 Memory usage per processor = 27.4725 Mbytes Step Atoms KinEng rke Volume 41 2 0.0010440444 0 4 CFD Coupling established at step 50 51 2 0.0011026021 0 4 Loop time of 0.000164933 on 16 procs for 10 steps with 2 atoms, finish time Mon Jan 4 21:27:41 2021 Pair time (%) = 2.85229e-06 (1.72936) Neigh time (%) = 0 (0) Comm time (%) = 2.79931e-05 (16.9724) Outpt time (%) = 2.44626e-05 (14.8318) Other time (%) = 0.000109625 (66.4664) Nlocal: 0.125 ave 1 max 0 min Histogram: 14 0 0 0 0 0 0 0 0 2 Nghost: 0.875 ave 2 max 0 min Histogram: 8 0 0 0 0 2 0 0 0 6 Neighs: 0.125 ave 1 max 0 min Histogram: 14 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 evolve done. DILUPBiCG: Solving for Ux, Initial residual = 0.0774, Final residual = 1.73435e-07, No Iterations 2 DILUPBiCG: Solving for Uy, Initial residual = 0.0774711, Final residual = 1.66837e-07, No Iterations 2 DILUPBiCG: Solving for Uz, Initial residual = 0.00540484, Final residual = 3.02225e-08, No Iterations 2 DICPCG: Solving for p, Initial residual = 0.683482, Final residual = 7.67603e-07, No Iterations 76 time step continuity errors : sum local = 3.50159e-12, global = -6.85906e-19, cumulative = -3.09243e-18 DICPCG: Solving for p, Initial residual = 0.158635, Final residual = 6.26763e-07, No Iterations 71 time step continuity errors : sum local = 1.10716e-12, global = -4.52701e-19, cumulative = -3.54513e-18 DICPCG: Solving for p, Initial residual = 0.00775263, Final residual = 8.5457e-07, No Iterations 46 time step continuity errors : sum local = 1.78028e-12, global = -7.52564e-19, cumulative = -4.2977e-18 DICPCG: Solving for p, Initial residual = 0.000685754, Final residual = 8.86278e-07, No Iterations 39 time step continuity errors : sum local = 1.85106e-12, global = 1.55092e-19, cumulative = -4.1426e-18 particleCloud.calcVelocityCorrection() DICPCG: Solving for phiIB, Initial residual = 0.232006, Final residual = 9.97434e-07, No Iterations 58 ExecutionTime = 0.61 s ClockTime = 0 s Time = 0.0018 Selected 0 cells for refinement out of 6268. Selected 0 split points out of a possible 264. Courant Number mean: 2.98648e-05 max: 0.00609245 - evolve() timeStepFraction() = 2 Starting up LIGGGHTS Executing command: 'run 10 ' run 10 Setting up run at Mon Jan 4 21:27:41 2021 Memory usage per processor = 27.4725 Mbytes Step Atoms KinEng rke Volume 51 2 0.0011026021 0 4 CFD Coupling established at step 60 61 2 0.0013665062 0 4 Loop time of 0.0001689 on 16 procs for 10 steps with 2 atoms, finish time Mon Jan 4 21:27:41 2021 Pair time (%) = 2.84938e-06 (1.68702) Neigh time (%) = 0 (0) Comm time (%) = 2.2641e-05 (13.405) Outpt time (%) = 2.47519e-05 (14.6548) Other time (%) = 0.000118658 (70.2533) Nlocal: 0.125 ave 2 max 0 min Histogram: 15 0 0 0 0 0 0 0 0 1 Nghost: 0.875 ave 2 max 0 min Histogram: 9 0 0 0 0 0 0 0 0 7 Neighs: 0.0625 ave 1 max 0 min Histogram: 15 0 0 0 0 0 0 0 0 1 Total # of neighbors = 1 Ave neighs/atom = 0.5 Neighbor list builds = 0 Dangerous builds = 0 LIGGGHTS finished nr particles = 2 evolve done. DILUPBiCG: Solving for Ux, Initial residual = 0.0407574, Final residual = 1.34213e-07, No Iterations 2 DILUPBiCG: Solving for Uy, Initial residual = 0.0407573, Final residual = 1.33765e-07, No Iterations 2 DILUPBiCG: Solving for Uz, Initial residual = 0.00212805, Final residual = 6.11183e-06, No Iterations 1 DICPCG: Solving for p, Initial residual = 0.590104, Final residual = 8.74187e-07, No Iterations 75 time step continuity errors : sum local = 2.41605e-12, global = 8.65567e-20, cumulative = -4.05605e-18 DICPCG: Solving for p, Initial residual = 0.125614, Final residual = 6.84886e-07, No Iterations 71 time step continuity errors : sum local = 8.37925e-13, global = -9.04246e-19, cumulative = -4.96029e-18 DICPCG: Solving for p, Initial residual = 0.00556585, Final residual = 8.71515e-07, No Iterations 40 time step continuity errors : sum local = 1.30144e-12, global = -3.0468e-19, cumulative = -5.26497e-18 DICPCG: Solving for p, Initial residual = 0.000493687, Final residual = 9.33962e-07, No Iterations 17 time step continuity errors : sum local = 1.36828e-12, global = -1.32178e-20, cumulative = -5.27819e-18 particleCloud.calcVelocityCorrection() DICPCG: Solving for phiIB, Initial residual = 0.14273, Final residual = 6.8454e-07, No Iterations 55 ExecutionTime = 0.64 s ClockTime = 0 s Time = 0.0021 Selected 0 cells for refinement out of 6268. Selected 0 split points out of a possible 264. Courant Number mean: 3.203e-05 max: 0.00645705 - evolve() timeStepFraction() = 2 Starting up LIGGGHTS Executing command: 'run 10 ' run 10 Setting up run at Mon Jan 4 21:27:41 2021 Memory usage per processor = 27.4725 Mbytes Step Atoms KinEng rke Volume 61 2 0.0013665062 0 4 CFD Coupling established at step 70 71 2 0.0017310924 0 4 Loop time of 0.000164298 on 16 procs for 10 steps with 2 atoms, finish time Mon Jan 4 21:27:41 2021 Pair time (%) = 2.84851e-06 (1.73374) Neigh time (%) = 0 (0) Comm time (%) = 2.85114e-05 (17.3534) Outpt time (%) = 2.4959e-05 (15.1912) Other time (%) = 0.00010798 (65.7216) Nlocal: 0.125 ave 1 max 0 min Histogram: 14 0 0 0 0 0 0 0 0 2 Nghost: 0.875 ave 2 max 0 min Histogram: 8 0 0 0 0 2 0 0 0 6 Neighs: 0.125 ave 1 max 0 min Histogram: 14 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 evolve done. DILUPBiCG: Solving for Ux, Initial residual = 0.0467523, Final residual = 1.61162e-07, No Iterations 2 DILUPBiCG: Solving for Uy, Initial residual = 0.0467522, Final residual = 1.58637e-07, No Iterations 2 DILUPBiCG: Solving for Uz, Initial residual = 0.00276303, Final residual = 8.17083e-06, No Iterations 1 DICPCG: Solving for p, Initial residual = 0.547815, Final residual = 7.49155e-07, No Iterations 76 time step continuity errors : sum local = 1.97094e-12, global = 1.86012e-19, cumulative = -5.09218e-18 DICPCG: Solving for p, Initial residual = 0.145234, Final residual = 7.14162e-07, No Iterations 71 time step continuity errors : sum local = 9.58344e-13, global = -7.04838e-19, cumulative = -5.79702e-18 DICPCG: Solving for p, Initial residual = 0.00643663, Final residual = 8.63424e-07, No Iterations 46 time step continuity errors : sum local = 1.36709e-12, global = -3.71912e-19, cumulative = -6.16893e-18 DICPCG: Solving for p, Initial residual = 0.000618496, Final residual = 9.52735e-07, No Iterations 26 time step continuity errors : sum local = 1.46845e-12, global = 2.44453e-19, cumulative = -5.92448e-18 particleCloud.calcVelocityCorrection() DICPCG: Solving for phiIB, Initial residual = 0.0901749, Final residual = 7.2405e-07, No Iterations 55 ExecutionTime = 0.67 s ClockTime = 0 s Time = 0.0024 Selected 0 cells for refinement out of 6268. Selected 0 split points out of a possible 264. Courant Number mean: 3.48095e-05 max: 0.00688244 - evolve() timeStepFraction() = 2 Starting up LIGGGHTS Executing command: 'run 10 ' run 10 Setting up run at Mon Jan 4 21:27:41 2021 Memory usage per processor = 27.4725 Mbytes Step Atoms KinEng rke Volume 71 2 0.0017310924 0 4 CFD Coupling established at step 80 81 2 0.002103895 0 4 Loop time of 0.000178583 on 16 procs for 10 steps with 2 atoms, finish time Mon Jan 4 21:27:41 2021 Pair time (%) = 2.83995e-06 (1.59027) Neigh time (%) = 0 (0) Comm time (%) = 2.90994e-05 (16.2946) Outpt time (%) = 3.0449e-05 (17.0503) Other time (%) = 0.000116195 (65.0648) Nlocal: 0.125 ave 1 max 0 min Histogram: 14 0 0 0 0 0 0 0 0 2 Nghost: 0.875 ave 2 max 0 min Histogram: 8 0 0 0 0 2 0 0 0 6 Neighs: 0.125 ave 1 max 0 min Histogram: 14 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 evolve done. DILUPBiCG: Solving for Ux, Initial residual = 0.0502708, Final residual = 1.71699e-07, No Iterations 2 DILUPBiCG: Solving for Uy, Initial residual = 0.0502707, Final residual = 1.68515e-07, No Iterations 2 DILUPBiCG: Solving for Uz, Initial residual = 0.00295856, Final residual = 8.76299e-06, No Iterations 1 DICPCG: Solving for p, Initial residual = 0.516727, Final residual = 8.16115e-07, No Iterations 76 time step continuity errors : sum local = 2.42532e-12, global = 2.42121e-19, cumulative = -5.68236e-18 DICPCG: Solving for p, Initial residual = 0.151486, Final residual = 7.18681e-07, No Iterations 71 time step continuity errors : sum local = 1.12843e-12, global = -7.23792e-19, cumulative = -6.40615e-18 DICPCG: Solving for p, Initial residual = 0.0066911, Final residual = 8.57773e-07, No Iterations 46 time step continuity errors : sum local = 1.55998e-12, global = -6.32346e-19, cumulative = -7.03849e-18 DICPCG: Solving for p, Initial residual = 0.000646209, Final residual = 9.37445e-07, No Iterations 39 time step continuity errors : sum local = 1.59589e-12, global = -1.82629e-19, cumulative = -7.22112e-18 particleCloud.calcVelocityCorrection() DICPCG: Solving for phiIB, Initial residual = 0.0476911, Final residual = 5.89507e-07, No Iterations 53 ExecutionTime = 0.69 s ClockTime = 0 s Time = 0.0027 Selected 0 cells for refinement out of 6268. Selected 0 split points out of a possible 264. Courant Number mean: 3.83218e-05 max: 0.00748651 - evolve() timeStepFraction() = 2 Starting up LIGGGHTS Executing command: 'run 10 ' run 10 Setting up run at Mon Jan 4 21:27:41 2021 Memory usage per processor = 27.4725 Mbytes Step Atoms KinEng rke Volume 81 2 0.002103895 0 4 CFD Coupling established at step 90 91 2 0.0024776697 0 4 Loop time of 0.00016368 on 16 procs for 10 steps with 2 atoms, finish time Mon Jan 4 21:27:41 2021 Pair time (%) = 2.85689e-06 (1.74541) Neigh time (%) = 0 (0) Comm time (%) = 2.83186e-05 (17.3012) Outpt time (%) = 2.49669e-05 (15.2535) Other time (%) = 0.000107538 (65.7) Nlocal: 0.125 ave 1 max 0 min Histogram: 14 0 0 0 0 0 0 0 0 2 Nghost: 0.875 ave 2 max 0 min Histogram: 8 0 0 0 0 2 0 0 0 6 Neighs: 0.125 ave 1 max 0 min Histogram: 14 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 evolve done. DILUPBiCG: Solving for Ux, Initial residual = 0.0479488, Final residual = 1.68129e-07, No Iterations 2 DILUPBiCG: Solving for Uy, Initial residual = 0.0479487, Final residual = 1.65139e-07, No Iterations 2 DILUPBiCG: Solving for Uz, Initial residual = 0.00268615, Final residual = 7.92939e-06, No Iterations 1 DICPCG: Solving for p, Initial residual = 0.509783, Final residual = 7.98607e-07, No Iterations 76 time step continuity errors : sum local = 2.59014e-12, global = -2.60685e-19, cumulative = -7.48181e-18 DICPCG: Solving for p, Initial residual = 0.151148, Final residual = 7.118e-07, No Iterations 71 time step continuity errors : sum local = 1.17674e-12, global = -2.77235e-19, cumulative = -7.75904e-18 DICPCG: Solving for p, Initial residual = 0.00668352, Final residual = 8.30888e-07, No Iterations 46 time step continuity errors : sum local = 1.57827e-12, global = 4.13821e-19, cumulative = -7.34522e-18 DICPCG: Solving for p, Initial residual = 0.00063374, Final residual = 9.86068e-07, No Iterations 20 time step continuity errors : sum local = 1.8075e-12, global = 7.24474e-20, cumulative = -7.27277e-18 particleCloud.calcVelocityCorrection() DICPCG: Solving for phiIB, Initial residual = 0.0362045, Final residual = 9.53632e-07, No Iterations 50 ExecutionTime = 0.73 s ClockTime = 0 s Time = 0.003 Selected 0 cells for refinement out of 6268. Selected 0 split points out of a possible 264. Courant Number mean: 4.21584e-05 max: 0.00822633 - evolve() timeStepFraction() = 2 Starting up LIGGGHTS Executing command: 'run 10 ' run 10 Setting up run at Mon Jan 4 21:27:41 2021 Memory usage per processor = 27.4725 Mbytes Step Atoms KinEng rke Volume 91 2 0.0024776697 0 4 CFD Coupling established at step 100 101 2 0.0028728033 0 4 Loop time of 0.000896791 on 16 procs for 10 steps with 2 atoms, finish time Mon Jan 4 21:27:42 2021 Pair time (%) = 3.5938e-06 (0.40074) Neigh time (%) = 0 (0) Comm time (%) = 3.05233e-05 (3.40361) Outpt time (%) = 0.000169469 (18.8973) Other time (%) = 0.000693205 (77.2984) Nlocal: 0.125 ave 1 max 0 min Histogram: 14 0 0 0 0 0 0 0 0 2 Nghost: 0.875 ave 2 max 0 min Histogram: 8 0 0 0 0 2 0 0 0 6 Neighs: 0.125 ave 1 max 0 min Histogram: 14 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 evolve done. DILUPBiCG: Solving for Ux, Initial residual = 0.0437948, Final residual = 1.60816e-07, No Iterations 2 DILUPBiCG: Solving for Uy, Initial residual = 0.0437948, Final residual = 1.58213e-07, No Iterations 2 DILUPBiCG: Solving for Uz, Initial residual = 0.00243064, Final residual = 7.17975e-06, No Iterations 1 DICPCG: Solving for p, Initial residual = 0.517604, Final residual = 7.75355e-07, No Iterations 76 time step continuity errors : sum local = 2.5931e-12, global = -8.17824e-21, cumulative = -7.28095e-18 DICPCG: Solving for p, Initial residual = 0.148878, Final residual = 6.97016e-07, No Iterations 71 time step continuity errors : sum local = 1.15656e-12, global = -7.49346e-20, cumulative = -7.35589e-18 DICPCG: Solving for p, Initial residual = 0.00660915, Final residual = 8.11699e-07, No Iterations 46 time step continuity errors : sum local = 1.54715e-12, global = 9.9911e-20, cumulative = -7.25598e-18 DICPCG: Solving for p, Initial residual = 0.000617411, Final residual = 9.75026e-07, No Iterations 20 time step continuity errors : sum local = 1.78935e-12, global = 1.8335e-19, cumulative = -7.07263e-18 particleCloud.calcVelocityCorrection() DICPCG: Solving for phiIB, Initial residual = 0.0325763, Final residual = 6.45413e-07, No Iterations 50 ExecutionTime = 0.77 s ClockTime = 1 s [...] Time = 0.0744 Selected 6 cells for refinement out of 6394. [13] [13] [13] --> FOAM FATAL IO ERROR: [13] error in IOstream "IOstream" for operation operator>>(Istream&, List&) : reading first token [13] [13] file: IOstream at line 0. [13] [13] From function void Foam::IOstream::fatalCheck(const char*) const [13] in file db/IOstreams/IOstreams/IOstream.C at line 109. [13] FOAM parallel run exiting [13] -------------------------------------------------------------------------- MPI_ABORT was invoked on rank 13 in communicator MPI_COMM_WORLD with errorcode 1. NOTE: invoking MPI_ABORT causes Open MPI to kill all MPI processes. You may or may not see output from other processes, depending on exactly when Open MPI kills them. -------------------------------------------------------------------------- finished CFDDEM run in parallel...