Add dpd thermostat

I am using the lastest LAMMPS to simulate a 3D couette flow.

How do I change the code here if I use dpd thermostat in Couette flow simulation? I’m going to finally look at the correlation between speed and viscosity.

Even if this is not an example, we would like to simply conclude the relationship between speed-viscosity with and without dpd thermostat. Please let me know if you have any other codes. Thank you.

LJ potentials

pair_style lj/cut 2.5
pair_coeff * * 1.0 1.0 2.5

define groups

region 1 block INF INF INF INF INF 1
group lower region 1
region 2 block INF INF INF INF 11 INF
group upper region 2
group boundary union lower upper
group flow subtract all boundary

set group lower type 1
set group upper type 3

initial velocities

#compute mobile flow temp
compute mobile flow temp/partial 1 1 1
velocity flow create 1.0 482748 temp mobile
fix 1 all nve
fix 2 flow temp/rescale 100 1.0 1.0 0.02 1.0
fix_modify 2 temp mobile

Couette flow

thermo 100
thermo_modify temp mobile

velocity lower set 0.0 0.0 0.0
velocity upper set 0.0 0.0 0.0
fix 3 boundary setforce 0.0 0.0 0.0

Run

dump 1 all atom 5 dump.flow
run 5000

timestep 0.005
velocity upper set 2.0 0.0 0.0
velocity lower set -2.0 0.0 0.0
run 60000
fix 4 flow ave/spatial 100 10 2000 z 0.0 4.0 vx ave running file vel.profile
run 40000

Please study the documentation: pair_style dpd command — LAMMPS documentation