Dear CFDEM users,
I have been running some IB-solver cases with a high demand on computational power. For that reason I had to start sometimes with a relatively small timestep before increasing the timestep of the simulation after some time.
To reload the DEM particles in my restarted simulation I use the read_dump command. Now, however, when implementing more than 10 particles with the read_dump command I obtain the following error within my LOG-file :
--- --- --- --- --- --- --- --- --- --- --- --- --- --- --- --- --- --- --- --- --- --- --- --- --- --- --- --- --- ---
Reading restart file ...
WARNING: Restart file version does not match LIGGGHTS version (../read_restart.cpp:502)
--> restart file = Version LIGGGHTS-PUBLIC 3.3.1, compiled 2016-05-20-09:49:08 by bnaar, git commit 99507f5217ad4fd4be5107d6262bd7274c9a5c03 based on LAMMPS 23 Nov 2013
--> LIGGGHTS = Version LIGGGHTS-PUBLIC 3.3.1, compiled 2016-05-20-09:57:43 by bnaar, git commit 99507f5217ad4fd4be5107d6262bd7274c9a5c03 based on LAMMPS 23 Nov 2013
orthogonal box = (-0.0024 -0.0024 0) to (0.0024 0.0024 0.0074)
2 by 2 by 4 MPI processor grid
46 atoms
read_dump ../DEM/post/dump.liggghts_init 0 x y z vx vy vz box yes replace yes
Scanning dump file ...
Reading snapshot from dump file ...
orthogonal box = (-0.0024 -0.0024 0) to (0.0024 0.0024 0.0074)
46 atoms before read
46 atoms in snapshot
0 atoms purged
46 atoms replaced
0 atoms trimmed
0 atoms added
46 atoms after read
neighbor 0.00008 bin
neigh_modify delay 0
# Material properties required for new pair styles
fix m1 all property/global youngsModulus peratomtype 5.e5
fix m2 all property/global poissonsRatio peratomtype 0.5
fix m3 all property/global coefficientRestitution peratomtypepair 1 1
fix m4 all property/global coefficientFriction peratomtypepair 1 0.25
# pair style
pair_style gran model hertz tangential history #Hertzian without cohesion
pair_coeff * *
# timestep, gravity
timestep 0.000005
# walls
fix xwalls1 all wall/gran model hertz tangential history primitive type 1 zcylinder 0.0024 0 0
Resetting global state of Fix history_xwalls1 Style property/atom from restart file info
Resetting per-atom state of Fix history_xwalls1 Style property/atom from restart file info
# change the particles density
set group all density 1100 # RBCs
Setting atom values ...
[12] #0 Foam::error::printStack(Foam::Ostream&)[10] #0 Foam::error::printStack(Foam::Ostream&) at ??:?
[12] #1 Foam::sigSegv::sigHandler(int) at ??:?
[10] #1 Foam::sigSegv::sigHandler(int) at ??:?
[12] #2 at ??:?
[10] #2 in "/lib64/libc.so.6"
[12] #3 LAMMPS_NS::FixPropertyAtom::unpack_restart(int, int) in "/lib64/libc.so.6"
[10] #3 LAMMPS_NS::FixPropertyAtom::unpack_restart(int, int) at ??:?
[12] #4 LAMMPS_NS::Modify::add_fix(int, char**, char*) at ??:?
[10] #4 LAMMPS_NS::Modify::add_fix(int, char**, char*) at ??:?
[12] #5 LAMMPS_NS::FixWallGran::post_create() at ??:?
[12] #6 LAMMPS_NS::Modify::add_fix(int, char**, char*) at ??:?
[10] #5 LAMMPS_NS::FixWallGran::post_create() at ??:?
[12] #7 LAMMPS_NS::Input::execute_command() at ??:?
[10] #6 LAMMPS_NS::Modify::add_fix(int, char**, char*) at ??:?
[12] #8 LAMMPS_NS::Input::file() at ??:?
[10] #7 LAMMPS_NS::Input::execute_command() at ??:?
[12] #9 LAMMPS_NS::Input::file(char const*) at ??:?
[10] #8 LAMMPS_NS::Input::file() at ??:?
[12] #10 Foam::twoWayMPI::twoWayMPI(Foam::dictionary const&, Foam::cfdemCloud&) at ??:?
[12] #11 Foam::dataExchangeModel::adddictionaryConstructorToTable::New(Foam::dictionary const&, Foam::cfdemCloud&) at ??:?
[10] #9 LAMMPS_NS::Input::file(char const*) at ??:?
[12] #12 Foam::dataExchangeModel::New(Foam::dictionary const&, Foam::cfdemCloud&) at ??:?
[10] #10 Foam::twoWayMPI::twoWayMPI(Foam::dictionary const&, Foam::cfdemCloud&) at ??:?
[12] #13 Foam::cfdemCloud::cfdemCloud(Foam::fvMesh const&) at ??:?
[10] #11 Foam::dataExchangeModel::adddictionaryConstructorToTable::New(Foam::dictionary const&, Foam::cfdemCloud&) at ??:?
[12] #14 Foam::cfdemCloudIB::cfdemCloudIB(Foam::fvMesh const&) at ??:?
[12] #15 at ??:?
[10] #12 Foam::dataExchangeModel::New(Foam::dictionary const&, Foam::cfdemCloud&)
[12] at ??:?
[12] #16 __libc_start_main at ??:?
[10] #13 Foam::cfdemCloud::cfdemCloud(Foam::fvMesh const&) in "/lib64/libc.so.6"
[12] #17
at ??:?
[10] #14 Foam::cfdemCloudIB::cfdemCloudIB(Foam::fvMesh const&)[12] at ??:?
[n11-41:61400] *** Process received signal ***
[n11-41:61400] Signal: Segmentation fault (11)
[n11-41:61400] Signal code: (-6)
AND MORE ...
--- --- --- --- --- --- --- --- --- --- --- --- --- --- --- --- --- --- --- --- --- --- --- --- --- --- --- --- --- ---
Does anyone know what is going on? Because the restarting of the files is working for cases with less than 10 particles.
Best regards,
Boris