double getCriticalVelocitySquared(TimeStep step) { final velocity = getCriticalVelocity(step); return velocity * velocity; }