I want to adopt pair_dpd from Lammps to implement Poiseuille flow. Since the DPD has its own way of controlling the temperature, I just use “addforce” to drive the fluid.
In my simulation, the velocity of the fluid keeps going up and shows no sign of stopping, what could be the reason for this? Below is my code. Thank you for your help.