Dear:
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.
Bests，
Shuai Zhang
Dear:
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.
Adding a force is adding an acceleration. So indeed you should expect an acceleration until the added force is equivalent to the dissipated energy.
Also, unless I miss my guess, your input looks more like a Couette flow than a Poisseuille flow to me. Please compare it to the inputs in the examples/flow folder.