I was following up to build a liquid droplet interacting with solid wall using many-body dissipative particle dynamics.
I could made a script to simulate a liquid droplet on static position but I was stuck how to set different attraction forces depending on the coordinates of droplets by the following correlation.
Asl (attraction forces) – x is the position of the liquid droplets
-8.662 x < -40
1.03 x 10-3x2 - 0.16 – 16.71 -40 ≤ x ≤ 40
-21.462 x > 40
My question is how to read the position of liquid droplets and how to update the attraction forces depending on the position into the script in LAMMPS.
I attached my script to let you make sense my problem.