Problem During Running Custom Code

Submitted by Yeasir Mohammad Akib on Sun, 06/19/2022 - 11:22


I am running the following code and facing some problem.

#Simple chute wear test

atom_style granular
atom_modify map array
boundary m m m
newton off

communicate single vel yes

units si

region domain block -0.5 0.5 -0.5 0.5 -0.5 0.5 units box
create_box 1 domain

neighbor 0.2 bin
neigh_modify delay 0

#Material properties required for new pair styles

fix m1 all property/global youngsModulus peratomtype 386e6
fix m2 all property/global poissonsRatio peratomtype 0.3
fix m3 all property/global coefficientRestitution peratomtypepair 1 0.64
fix m4 all property/global coefficientFriction peratomtypepair 1 0.6

#New pair style
pair_style gran model hertz tangential history #Hertzian without cohesion
pair_coeff * *

timestep 0.00001

fix gravi all gravity 9.81 vector 0.0 0.0 -1.0

#the chute
fix box all mesh/surface file meshes/box.stl type 1
fix granwalls all wall/gran model hertz tangential history mesh n_meshes 1 meshes box

#distributions for insertion
fix pts1 all particletemplate/sphere 15485863 atom_type 1 density constant 2500 radius constant 0.0015
fix pts2 all particletemplate/sphere 15485867 atom_type 1 density constant 2500 radius constant 0.0025
fix pdd1 all particledistribution/discrete 32452843 2 pts1 0.3 pts2 0.7

##region and insertion
#group nve_group region domain
#region bc block -0.015 0.015 -0.01 0.01 0.01 0.05 units box

fix ins_mesh all mesh/surface/planar file meshes/insertionface.stl type 1

#particle insertion
fix ins all insert/stream seed 32452867 distributiontemplate pdd1 &
nparticles 6000 massrate 0.1 overlapcheck yes all_in yes vel constant 0.0 0.0 -1.0 &
insertion_face ins_mesh extrude_length 0.25

#apply nve integration to all particles that are inserted as single particles
fix integr all nve/sphere #Changing nve_group- all

#output settings, include total thermal energy
compute 1 all erotate/sphere
thermo_style custom step atoms ke c_1 vol
thermo 1000
thermo_modify lost ignore norm no

#insert the first particles so that dump is not empty
run 1
dump dmp all custom/vtk 200 post/chute_*.vtk id type type x y z ix iy iz vx vy vz fx fy fz omegax omegay omegaz radius
dump dumpstress all mesh/gran/VTK 200 post/mesh_*.vtk stress wear cad

#insert particles
run 100000 upto
unfix ins


Error message is:


LIGGGHTS (Version LIGGGHTS-PUBLIC 3.8.0, compiled 2022-06-13-16:05:55 by slm_lab_utrgv, git commit 86544c3f1215bc9744d7929e3eacf53ba086af9e)
Created orthogonal box = (-0.5 -0.5 -0.5) to (0.5 0.5 0.5)
1 by 2 by 2 MPI processor grid

Reading STL file 'meshes/box.stl' (mesh processing step 1/3)
Fix particledistribution/discrete (id pdd1): distribution based on mass%:
pts1: d=3.000000e-03 (max. bounding sphere) mass%=30.000000%
pts2: d=5.000000e-03 (max. bounding sphere) mass%=70.000000%
Fix particledistribution/discrete (id pdd1): distribution based on number%:
pts1: d=3.000000e-03 (max. bounding sphere) number%=66.489362%
pts2: d=5.000000e-03 (max. bounding sphere) number%=33.510638%

Reading STL file 'meshes/insertionface.stl' (mesh processing step 1/3)
Setting up run at Sun Jun 19 04:21:26 2022

Import and parallelization of mesh box containing 2 triangle(s) successful
Import and parallelization of mesh ins_mesh containing 2 triangle(s) successful
INFO: Particle insertion ins: 319.158713 particles every 25000 steps - particle rate 1276.634850 (mass rate 1.000000e-01)
6000 particles (mass 4.699856e-01) within 450000 steps
Memory usage per processor = 89.5836 Mbytes
Step Atoms KinEng 1 Volume
0 0 0 0 83440.821
WARNING: Fix insert/stream: Extruded insertion face extends outside domain, may not insert all particles correctly (../fix_insert_stream.cpp:561)
ERROR: Fix insert/stream (id ins): insertion volume could not be distributed properly in parallel. Bad decomposition or insertion face extrusion is too small or outside domain (../fix_insert_stream.cpp:552)
`���}U (../fix_insert_stream.cpp:552)
Primary job terminated normally, but 1 process returned
a non-zero exit code. Per user-direction, the job has been aborted.
mpirun detected that one or more processes exited with non-zero status, thus causing
the job to be terminated. The first process to do so was:

Process name: [[43806,1],0]
Exit code: 1

Can someone please help?