plane poiseuille flow with LIGGGHTS

Hi,

I want guidance related to setting up plane poiseuille flow with LIGGGHTS.

I've added "fix 9 all addforce 0.1 0.0 0.0"
Can anyone give me some guidance to proceed furthur ?

I'm attaching my input file here

"
atom_style granular
atom_modify map array
communicate single vel yes

boundary p f f
newton off

units si

region reg block 0 0.4 0 0.1 0 0.1 units box
create_box 1 reg

neighbor 0.001 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
#mm

# walls
#fix xwalls1 all wall/gran model hertz tangential history primitive type 1 xplane 0.0
#fix xwalls2 all wall/gran model hertz tangential history primitive type 1 xplane 0.4
fix ywalls1 all wall/gran model hertz tangential history primitive type 1 yplane 0.0
fix ywalls2 all wall/gran model hertz tangential history primitive type 1 yplane 0.1
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 0.1

# particle distributions and insertion
fix pts1 all particletemplate/sphere 15485863 atom_type 1 density constant 2000 radius constant 0.0005
fix pdd1 all particledistribution/discrete 15485867 1 pts1 1.0

fix ins all insert/pack seed 32452843 distributiontemplate pdd1 vel constant 0.0 0.0 0.0 insert_every once overlapcheck yes all_in yes particles_in_region 1000 region reg

#poiseuille flow
fix 9 all addforce 0.1 0.0 0.0

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

# screen output
compute 1 all erotate/sphere
thermo_style custom step atoms ke c_1 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/vtk 800 post/newModels_*.vtk id type x y z ix iy iz vx vy vz fx fy fz omegax omegay omegaz radius
#dump dmp all custom 5000 post/dump.liggghts_init id type type x y z ix iy iz vx vy vz fx fy fz omegax omegay omegaz radius

run 10000 upto
write_restart post/restart/liggghts.restart
"