I am dealing with steady shear of polymer blends use fix/form command and DPD thermostat. If got a DPD thermostat then the fix/sllod cannot be used (am I right ?) which brings the troubles e.g. the temp/deform is not correct to dump the system temperature under high rate. I used the temp/profile to obtain the real temperature. However I am quite doubt whether the flow kinetic components are correctly excluded when compute the system pressure and stress tensors following the “compute press all pressure mytemp”, though it seems the energy tensors have been dealt in the temperature compute command stated in the manual, considering this is not a fix/sllod condition.
Following is the main script for the Shear and Temperature & Pressure Computation.
/*~~~~~~~~~~~~~
thermo_style custom step temp cpu etotal pxx pyy pzz lx ly lz xy xz yz
thermo 100
Thanks for you reply. In fact I want to know whether the temp/profile used in the pressure calculation will subtract out the flow velocity as in the temp/profile calculation.