Please check the following script. It is running fine with atom_style sph but with hybrid sph peri it gives Segmentation fault error.
Many thanks
echo screen
processors 2 2 2
atom_style hybrid sph peri
atom_modify map array
newton off
boundary f f f
units si
variable h equal 0.0002
variable c equal 14.80 # soundspeed for Tait’s EOS
variable dt equal 0.0000001*{h}/{c} # CFL criterion
region box block 0 0.04 0 0.04 0 0.02 units box
create_box 2 box
region drop sphere 0.02 0.02 0.0065 0.00135
region water block 0.015 0.025 0.015 0.025 0.0046 0.005
lattice bcc 0.00015
create_atoms 1 region drop
lattice sc 0.00015
create_atoms 2 region water
#lattice sc 1.0
set region drop type 1
set region water type 2
group drop type 1
group water type 2
mass 1 1.75972385E-9 #1.0767E-07
mass 2 1.76714587E-9 #1.09721E-07
set group drop sph/rho 995.8
set group water sph/rho 1000
pair_style hybrid/overlay sph/rhosum 1 sph/taitwater peri/pmb
pair_coeff * * sph/taitwater 995.8 {c} 1 {h}
pair_coeff 1 1 sph/rhosum {h}
pair_coeff 2 2 peri/pmb 6.2696893e24 0.00045001 0.00015 0.0
fix gfix water gravity 9.81 vector 0 0 -1
fix gfix drop gravity 9.81 vector 0 0 -1
velocity drop ramp vz -2.67 -2.67 z 0.0185 0.0115
velocity water set 0.0 0.0 0.0 sum yes units box
compute rho_peratom all sph/rho/atom
compute e_peratom all sph/e/atom
compute esph all reduce sum c_e_peratom
compute ke all ke
variable etot equal c_esph+c_ke+f_gfix
fix dtfix all dt/reset 1 NULL {dt} {h} units box
fix integrate_drop_fix drop sph
fix integrate_bc_fix water sph/stationary
fix F1 water nve
compute C1 water damage/atom
thermo_style custom step time atoms c_ke c_esph v_etot f_gfix press f_dtfix
thermo 10
thermo_modify lost ignore norm no
dump D1 all custom 100 …/perisph/post/dump.drop id type xs ys zs vx vy vz fx fy fz
variable skin equal 0.3*{h}
neighbor {skin} bin
neigh_modify every 1 delay 0 check no
run 5000