at one point the position/force vectors get filled with nan. Detect the divergence before that, clue: probably best to act directly inside the solver or create a new component in the scene which watch the solver residual or state vectors.