Hi, everyone:
I want to bulid a sandwich structure with two solids on both sides and liquid in the moddile. The two solid walls need to be thermostated as 100 K and 140 K, and I want to see the tempureture distribution of the liquid region. I wonder which thermostatting method is best.
- Firstly, Itried temp/rescale, and the enegry is output. However, the linear temperature distribution can not be obtained after 3000000 steps running. What is more, unger this situation, the momentum conversion can be remained or not ?
fix 1 all nve
run 200000
fix 2 cold temp/rescale 10 100 100 0.5 1.0
fix_modify 2 energy yes
fix 3 hot temp/rescale 10 140 140 0.5 1.0
fix_modify 3 energy yes
thermo_style custom atoms step temp pe ke etotal press vol f_2 f_3
- I tried
fix 2 cold nvt temp 100 100 100 drag 0.2
fix_modify 2 energy yes
fix 3 hot nvt temp 140 140 100 drag 0.2
fix_modify 3 energy yes
thermo_style custom atoms step temp pe ke etotal press vol f_2 f_3