I want to apply time varying body force on the particles in poiseuille using ‘fix addforce’ command. The force can be sinusoidal or square function. I think the former is difficult to implement but there should be a way to implement the square function. But I am not sure about how to do this. Any inputs will be valuable.