mirror of
https://github.com/ParticulateFlow/CFDEMcoupling-PFM.git
synced 2025-12-08 06:37:44 +00:00
65 lines
2.2 KiB
Plaintext
65 lines
2.2 KiB
Plaintext
# Pour granular particles into a cylinder, then induce flow
|
|
echo both
|
|
atom_style granular
|
|
atom_modify map array
|
|
communicate single vel yes
|
|
|
|
boundary m m m
|
|
newton off
|
|
|
|
units cgs
|
|
|
|
region reg block -1.5 1.5 -1.5 1.5 -0.1 5.54 units box
|
|
create_box 1 reg
|
|
|
|
neighbor 0.1 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
|
|
pair_coeff * *
|
|
|
|
# timestep, gravity
|
|
timestep 0.00001
|
|
fix gravi all gravity 981 vector 0.0 0.0 -1.0
|
|
|
|
# walls
|
|
fix zwalls1 all wall/gran model hertz tangential history primitive type 1 zplane 0.0
|
|
fix zwalls2 all wall/gran model hertz tangential history primitive type 1 zplane 5.53
|
|
fix cylwalls all wall/gran model hertz tangential history primitive type 1 zcylinder 1.385 0. 0.
|
|
|
|
# particle distributions and insertion
|
|
region bc cylinder z 0.0 0.0 1.2 0. 5.5 units box
|
|
fix pts1 all particletemplate/sphere 1 atom_type 1 density constant 2. radius constant 0.05
|
|
fix pdd1 all particledistribution/discrete 1 1 pts1 1.0
|
|
|
|
fix ins all insert/pack seed 100001 distributiontemplate pdd1 vel constant 0. 0. -100. insert_every once overlapcheck yes all_in yes particles_in_region 10000 region bc
|
|
|
|
# apply nve integration to all particles that are inserted as single particles
|
|
fix integr all nve/sphere
|
|
|
|
# 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
|
|
|
|
# insert the first particles so that dump is not empty
|
|
run 1
|
|
dump dmp all custom 5000 post/dump.liggghts_init id type x y z 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 10000 upto
|
|
write_restart post/restart/liggghts.restart
|