Unstable velocity verlet algorithm

[franziss]

I implemented a velocity verlet algorithm, and use it to simulate movement of particles. Depending on the parameters, sometimes the simulation of the movment of the particles is stable, while other times, the simulation may go haywire and the position of the particles becomes exponential.

Is there a way to prevent unstable simulation? Is there any verlet algorithm that guarantees stability? E.g. a regularization parameter in the verlet algo? Or should I try constraint algorithm? I have a look at constraint algorithm, it seems that the distance between two particles is fixed to a constraint. Does constraint algorithm allows a maximum threshold, i.e. the distance between two particles cannot exceed a threshold?

Thank you for your help.

 

[LightHero]

The velocity Verlet algorithm does not guarantee stability, but there are ways you can improve the stability of your simulations. For example, you can reduce the time step or increase the damping factor to make the simulation more stable. You can also use a regularization parameter in the Verlet algorithm, such as adding a small amount of damping or friction to the system. Alternatively, you can use a constraint algorithm that constrains the distance between particles to a maximum threshold. This will ensure that the particles don't move too far from their initial positions and will help to keep the system stable.