I’m developing a simulation software for flight simulation.

I have problems with solver stability.

I have a rigidbody subject to gravity and aerodynamic forces. My timestep is fixed to 0.01s and lift and drag forces are computed at timestep n using body linear and angular velocity. (I compute AoA and relative speed body/fluid in wing coordinate system) and apply to the actor with simply functions: AddForceAtLocalPosition for example.

after few timesteps (200/300) the simulation diverge and velocities goes to infinite. If I lower the timestep obviously the instability come up at a later moment.

I’ve tried to use function as:

ComputeVelocityDeltaFromImpulse and use velocity at timestep n+1 to compute forces in order to try an “implicit method” but this not improve the simulation stability.

I understood that Physx solver use semi implicit integrator method but i don’t if i can modify some parameters about the solver.

Any method to improve solver stability without reducing the timestep too much in this case? There a lot of simulator and i can’t figure out how they deal such stability problems for realtime simulations?

Thank you,