************************************ using userdefined CFDEM_SRC_DIR = /public1/home/scb0030/cfdem_xia/openfoam/CFDEM/CFDEMcoupling-PUBLIC-5.0/src using userdefined CFDEM_SOLVER_DIR = /public1/home/scb0030/cfdem_xia/openfoam/CFDEM/CFDEMcoupling-PUBLIC-5.0/applications/solvers using userdefined CFDEM_DOC_DIR = /public1/home/scb0030/cfdem_xia/openfoam/CFDEM/CFDEMcoupling-PUBLIC-5.0/doc using userdefined CFDEM_UT_DIR = /public1/home/scb0030/cfdem_xia/openfoam/CFDEM/CFDEMcoupling-PUBLIC-5.0/applications/utilities using userdefined CFDEM_TUT_DIR = /public1/home/scb0030/cfdem_xia/openfoam/CFDEM/CFDEMcoupling-PUBLIC-5.0/tutorials using userdefined CFDEM_LAMMPS_LIB_DIR = /public1/home/scb0030/cfdem_xia/openfoam/LIGGGHTS/LIGGGHTS-PUBLIC/src/../lib using userdefined CFDEM_LIGGGHTS_LIB_PATH = /public1/home/scb0030/cfdem_xia/openfoam/LIGGGHTS/LIGGGHTS-PUBLIC/src using userdefined CFDEM_ADD_LIBS_DIR = /public1/home/scb0030/cfdem_xia/openfoam/CFDEM/CFDEMcoupling-PUBLIC-5.0/src/lagrangian/cfdemParticle/etc/addLibs_universal using userdefined CFDEM_ADD_LIBS_NAME = additionalLibs_5.0 using userdefined CFDEM_LIB_DIR = /public1/home/scb0030/cfdem_xia/openfoam/CFDEM/CFDEMcoupling-PUBLIC-5.0/platforms/linux64IccDPInt32Opt/lib using userdefined CFDEM_USER_LIB_DIR = /public1/home/scb0030/cfdem_xia/openfoam/CFDEM/scb0030-PUBLIC-5.0/platforms/linux64IccDPInt32Opt/lib using userdefined CFDEM_APP_DIR = /public1/home/scb0030/cfdem_xia/openfoam/CFDEM/CFDEMcoupling-PUBLIC-5.0/platforms/linux64IccDPInt32Opt/bin using userdefined CFDEM_USER_APP_DIR = /public1/home/scb0030/cfdem_xia/openfoam/CFDEM/scb0030-PUBLIC-5.0/platforms/linux64IccDPInt32Opt/bin ************************************ /*---------------------------------------------------------------------------*\ | ========= | | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | | \\ / O peration | Version: 5.0 | | \\ / A nd | Web: www.OpenFOAM.org | | \\/ M anipulation | | \*---------------------------------------------------------------------------*/ Build : 5.0 Exec : blockMesh Date : Jun 25 2021 Time : 09:54:48 Host : "gb0807.para.bscc" PID : 91785 I/O : uncollated Case : /public1/home/scb0030/Example/ErgunTestMPI/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 "/public1/home/scb0030/Example/ErgunTestMPI/CFD/constant/polyMesh" Creating block mesh from "/public1/home/scb0030/Example/ErgunTestMPI/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 : 8 Number of boundary faces : 14 Number of defined boundary faces : 14 Number of undefined boundary faces : 0 Checking patch -> block consistency Creating block offsets Creating merge list . Creating polyMesh from blockMesh Creating points with scale 1 Block 0 cell size : i : 0.001725 .. 0.001725 j : 0.001725 .. 0.001725 k : 0.0023041667 .. 0.0023041667 Block 1 cell size : i : 0.001725 .. 0.001725 j : 0.0010104816 .. 0.0010104816 k : 0.0023041667 .. 0.0023041667 Block 2 cell size : i : 0.0010104816 .. 0.0010104816 j : 0.0027052731 .. 0.0027052731 k : 0.0023041667 .. 0.0023041667 Block 3 cell size : i : 0.0027052731 .. 0.0027052731 j : 0.0010104816 .. 0.0010104816 k : 0.0023041667 .. 0.0023041667 Block 4 cell size : i : 0.0010104816 .. 0.0010104816 j : 0.001725 .. 0.001725 k : 0.0023041667 .. 0.0023041667 Creating cells Creating patches Writing polyMesh ---------------- Mesh Information ---------------- boundingBox: (-0.0138 -0.0138 0) (0.0138 0.0138 0.0553) nPoints: 5225 nCells: 4608 nFaces: 14400 nInternalFaces: 13248 ---------------- Patches ---------------- patch 0 (start: 13248 size: 192) name: inlet patch 1 (start: 13440 size: 192) name: outlet patch 2 (start: 13632 size: 768) name: wall End /*---------------------------------------------------------------------------*\ | ========= | | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | | \\ / O peration | Version: 5.0 | | \\ / A nd | Web: www.OpenFOAM.org | | \\/ M anipulation | | \*---------------------------------------------------------------------------*/ Build : 5.0 Exec : decomposePar Date : Jun 25 2021 Time : 09:54:48 Host : "gb0807.para.bscc" PID : 91792 I/O : uncollated Case : /public1/home/scb0030/Example/ErgunTestMPI/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 = 1152 Number of faces shared with processor 1 = 192 Number of faces shared with processor 2 = 192 Number of processor patches = 2 Number of processor faces = 384 Number of boundary faces = 288 Processor 1 Number of cells = 1152 Number of faces shared with processor 0 = 192 Number of faces shared with processor 3 = 192 Number of processor patches = 2 Number of processor faces = 384 Number of boundary faces = 288 Processor 2 Number of cells = 1152 Number of faces shared with processor 0 = 192 Number of faces shared with processor 3 = 192 Number of processor patches = 2 Number of processor faces = 384 Number of boundary faces = 288 Processor 3 Number of cells = 1152 Number of faces shared with processor 1 = 192 Number of faces shared with processor 2 = 192 Number of processor patches = 2 Number of processor faces = 384 Number of boundary faces = 288 Number of processor faces = 768 Max number of cells = 1152 (0% above average 1152) Max number of processor patches = 2 (0% above average 2) Max number of faces between processors = 384 (0% above average 384) Time = 0 Processor 0: field transfer Processor 1: field transfer Processor 2: field transfer Processor 3: field transfer End /*---------------------------------------------------------------------------*\ | ========= | | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | | \\ / O peration | Version: 5.0 | | \\ / A nd | Web: www.OpenFOAM.org | | \\/ M anipulation | | \*---------------------------------------------------------------------------*/ Build : 5.0 Exec : /public1/home/scb0030/cfdem_xia/openfoam/CFDEM/CFDEMcoupling-PUBLIC-5.0/platforms/bin/cfdemSolverPiso -parallel Date : Jun 25 2021 Time : 09:54:50 Host : "gb0807.para.bscc" PID : 91831 I/O : uncollated Case : /public1/home/scb0030/Example/ErgunTestMPI/CFD nProcs : 4 Slaves : 3 ( "gb0807.para.bscc.91832" "gb0807.para.bscc.91833" "gb0807.para.bscc.91834" ) 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 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 momentum exchange field Ksl Reading voidfraction field voidfraction = (Vgas/Vparticle) Creating density field rho Reading particle velocity field Us Reading/calculating face flux field phi Generating interstitial face flux field phiByVoidfraction Selecting incompressible transport model Newtonian Selecting turbulence model type laminar Selecting laminar stress model Stokes No finite volume options present Reading g Selecting locateModel engine 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 2021-06-24-14:42:24 by scb0030, git commit unknown) log ../DEM/log.liggghts thermo_log ../DEM/post/thermo.txt atom_style granular atom_modify map array communicate single vel yes boundary m m m newton off units si processors 2 2 1 # read the restart file read_restart ../DEM/post/restart/liggghts.restart Reading restart file ... WARNING: Restart file version does not match LIGGGHTS version (../read_restart.cpp:507) --> restart file = Version LIGGGHTS-PUBLIC 3.8.0, compiled 2021-06-24-14:33:37 by scb0030, git commit unknown --> LIGGGHTS = Version LIGGGHTS-PUBLIC 3.8.0, compiled 2021-06-24-14:42:24 by scb0030, git commit unknown version 3 8 version 3 8 version 3 8 version 3 8 WARNING: Restart file used different # of processors (../read_restart.cpp:574) WARNING: Restart file used different 3d processor grid (../read_restart.cpp:587) orthogonal box = (-0.015 -0.015 -0.001) to (0.015 0.015 0.0554) 2 by 2 by 1 MPI processor grid 9996 atoms neighbor 0.0005 bin neigh_modify delay 0 # Material properties required for granular pair styles fix m1 all property/global youngsModulus peratomtype 5.e6 fix m2 all property/global poissonsRatio peratomtype 0.45 fix m3 all property/global coefficientRestitution peratomtypepair 1 0.3 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.00001 fix gravi all gravity 9.81 vector 0.0 0.0 -1.0 # walls fix zwalls1 all wall/gran model hertz tangential history primitive type 1 zplane 0.0 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 fix zwalls2 all wall/gran model hertz tangential history primitive type 1 zplane 0.0553 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 fix cylwalls all wall/gran model hertz tangential history primitive type 1 zcylinder 0.01385 0. 0. 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 # cfd coupling fix cfd all couple/cfd couple_every 100 mpi couple_every as specified in LIGGGHTS is overriden by calling external program (../fix_cfd_coupling.cpp:114) fix cfd2 all couple/cfd/force/implicit #fix cfd2 all couple/cfd/force/accumulator RongDrag 10 1.5e-3 #fix cfd2 all couple/cfd/force/implicit/accumulated #CrankNicolson 0.5 # apply nve integration to all particles that are inserted as single particles fix integr all nve/sphere # center of mass compute centerOfMass all com # compute explicit dragforce compute explDrag all reduce update_on_run_end yes sum f_dragforce[1] f_dragforce[2] f_dragforce[3] # sum of explicit and implicit drag force given from CFD to DEM variable totalDragX equal f_cfd2[1] variable totalDragY equal f_cfd2[2] variable totalDragZ equal f_cfd2[3] # explicit drag force given from CFD to DEM variable explicitDragX equal c_explDrag[1] variable explicitDragY equal c_explDrag[2] variable explicitDragZ equal c_explDrag[3] variable time equal step*dt fix extra all print 10 "${time} ${explicitDragX} ${explicitDragY} ${explicitDragZ} ${totalDragX} ${totalDragY} ${totalDragZ}" file ../DEM/post/forces.txt title "# time expDrag(X Y Z) totDrag(X Y Z)" # screen output compute rke all erotate/sphere thermo_style custom step atoms ke c_rke vol c_centerOfMass[3] c_explDrag[1] c_explDrag[2] c_explDrag[3] f_cfd2[1] f_cfd2[2] f_cfd2[3] thermo 10 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) dump dmp all custom 5000 ../DEM/post/dump*.liggghts_run id type x y z vx vy vz fx fy fz f_dragforce[1] f_dragforce[2] f_dragforce[3] radius #dump dmp all custom/vtk 5000 ../DEM/post/dump_liggghts_run.*.vtk id type x y z vx vy vz fx fy fz f_dragforce[1] f_dragforce[2] f_dragforce[3] radius run 0 Resetting global state of Fix contacthistory Style contacthistory from restart file info Resetting per-atom state of Fix contacthistory Style contacthistory from restart file info Setting up run at Fri Jun 25 09:54:50 2021 Memory usage per processor = 11.7605 Mbytes Step Atoms KinEng rke Volume centerOf explDrag explDrag explDrag cfd2[1] cfd2[2] cfd2[3] 10000 9996 3.6856019e-08 3.1131695e-08 5.076e-05 0.0075988038 0 0 0 0 0 0 Loop time of 5.48363e-06 on 4 procs for 0 steps with 9996 atoms, finish time Fri Jun 25 09:54:50 2021 Pair time (%) = 0 (0) Neigh time (%) = 0 (0) Comm time (%) = 0 (0) Outpt time (%) = 0 (0) Other time (%) = 5.48363e-06 (100) Nlocal: 2499 ave 2508 max 2492 min Histogram: 1 0 0 2 0 0 0 0 0 1 Nghost: 743.25 ave 750 max 736 min Histogram: 1 0 0 0 0 2 0 0 0 1 Neighs: 17413 ave 17521 max 17335 min Histogram: 1 0 1 0 1 0 0 0 0 1 Total # of neighbors = 69652 Ave neighs/atom = 6.96799 Neighbor list builds = 0 Dangerous builds = 0 Selecting IOModel basicIO Selecting probeModel off Selecting voidFractionModel divided Selecting averagingModel dense Selecting clockModel off start clock measurement at t >0.001 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.0 OF build: 5.0 CFDEM build: [1] #0 Foam::error::printStack(Foam::Ostream&) at ??:? [3] #0 Foam::error::printStack(Foam::Ostream&) at ??:? [2] #0 Foam::error::printStack(Foam::Ostream&) at ??:? [0] #0 Foam::error::printStack(Foam::Ostream&) at ??:? [3] #1 Foam::sigFpe::sigHandler(int) at ??:? [1] #1 Foam::sigFpe::sigHandler(int) at ??:? [0] #1 Foam::sigFpe::sigHandler(int) at ??:? [2] #1 Foam::sigFpe::sigHandler(int) at ??:? [3] #2 ? in "/lib64/libc.so.6" [1] #2 ? in "/lib64/libc.so.6" [0] #2 ? in "/lib64/libc.so.6" [2] #2 ? in "/lib64/libc.so.6" [3] #3 Foam::twoWayMPI::setCG() const at ??:? [1] #3 Foam::twoWayMPI::setCG() const at ??:? [0] #3 Foam::twoWayMPI::setCG() const at ??:? [2] #3 Foam::twoWayMPI::setCG() const at ??:? [3] #4 Foam::cfdemCloud::cfdemCloud(Foam::fvMesh const&) at ??:? [1] #4 Foam::cfdemCloud::cfdemCloud(Foam::fvMesh const&) at ??:? [0] #4 Foam::cfdemCloud::cfdemCloud(Foam::fvMesh const&) at ??:? [3] #5 ? at ??:? [2] #4 Foam::cfdemCloud::cfdemCloud(Foam::fvMesh const&) at ??:? [1] #5 ? at ??:? [0] #5 ? at ??:? [3] #6 __libc_start_main in "/lib64/libc.so.6" [2] #5 ? at ??:? [1] #6 __libc_start_main in "/lib64/libc.so.6" [0] #6 __libc_start_main in "/lib64/libc.so.6" [3] #7 ? at ??:? [2] #6 __libc_start_main in "/lib64/libc.so.6" [1] #7 ? at ??:? [0] #7 ? at ??:? [2] #7 ? at ??:? srun: error: gb0807: tasks 0-3: Floating point exception /*---------------------------------------------------------------------------*\ | ========= | | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | | \\ / O peration | Version: 5.0 | | \\ / A nd | Web: www.OpenFOAM.org | | \\/ M anipulation | | \*---------------------------------------------------------------------------*/ Build : 5.0 Exec : reconstructParMesh Date : Jun 25 2021 Time : 09:54:50 Host : "gb0807.para.bscc" PID : 91878 I/O : uncollated Case : /public1/home/scb0030/Example/ErgunTestMPI/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 This is an experimental tool which tries to merge individual processor meshes back into one master mesh. Use it if the original master mesh has been deleted or if the processor meshes have been modified (topology change). This tool will write the resulting mesh to a new time step and construct xxxxProcAddressing files in the processor meshes so reconstructPar can be used to regenerate the fields on the master mesh. Not well tested & use at your own risk! Merge tolerance : 1e-07 Write tolerance : 1e-08 Doing geometric matching on correct procBoundaries only. This assumes a correct decomposition. Found 4 processor directories Reading database "CFD/processor0" Reading database "CFD/processor1" Reading database "CFD/processor2" Reading database "CFD/processor3" End.