I was trying to simulate a Couette flow & oscillatory Couette flow and wanted to know the thermal temperature at different regions(layers along the height of channel).

1) I was using temp/region to do it but I wasnt explicitly sure if it removed the advection velocity while providing the temperature.

2) Also I was assuming the temperature provided by default thermo command and using thermo_style temp would be the same but it comes out to be different.

Would appreciate any comments provided.

Thank you.

BR/Joseph

ps: input script for reference

.

region fluid block INF INF 2.39 52.5 INF INF units box #<<<<< CHNAGE DIMENSION CORRSPONDING TO PENETRATION DEPTH

region btm_wall block INF INF INF 0.95 INF INF units box

region top_wall block INF INF 53.86 INF INF INF units box #<<<<< CHNAGE DIMENSION CORRSPONDING TO PENETRATION DEPTH

.

# initial velocities

compute mobile flow temp

velocity flow create 1.10 482748 temp mobile units box

thermo_modify temp mobile

# thermostating

compute mobile_thermo all temp/partial 0 0 1

fix 3 all langevin 1.10 1.10 1.0 699483

fix 1 flow nve

fix 2 lower nve

fix_modify 3 temp mobile_thermo

thermo_modify temp mobile_thermo

# Couette flow

velocity lower set 0.0 0.0 0.0 units box

fix 8 upper move wiggle 10.0 0.0 0.0 20.0 units box #<<<<< CHANGE AMPLITUDE AND TIMEPERIOD

fix 4 boundary setforce 0.0 0.0 0.0

#variable pxy equal pxy

region wall_fluid_1 block INF INF 53.135 53.385 INF INF units box

#group gwall_fluid_1 region wall_fluid_1

#compute temp_tst_1 gwall_fluid_1 temp/com

compute temp_tst_1 all temp/region wall_fluid_1

.

# Run

timestep 0.002

thermo 1000

#thermo_style custom step temp evdwl ecoul epair ebond edihed eimp eangle emol etotal

#thermo_style custom step temp etotal c_temp_tst_1 c_temp_tst_2 c_temp_tst_3 c_temp_tst_4 c_temp_tst_5 c_temp_tst_6 c_temp_tst_7 c_temp_tst_8 c_temp_tst_9

.