Hey all,

I’m a bit confused as to modeling the following situation.

suppose I’m dropping a particle of mass M from height H. The forces acting upon it include gravity and drag.

now, suppose i know the rate at which this particle ALSO loses its mass during that fall, say, dM/dt = -f(t).

how do I insert that into MD lammps simulation ?

let’s say I begin with F=dP/dt=m*dv/dt + v*dm/dt --> mg + Fdrag = m*dv/dt + v*dm/dt

so apparently it seems reasonable to work with this formulation (variable mass problem).

and yet, not sure as to how actually implement this into lammps ? can lammps work with variable particle mass ?

Thanks in advance !!

p.s i’m new into MD and just starting to mingle in this world