TidalPy
TidalPy copied to clipboard
RadialSolver: Dynamic Liquid Layer Stability Issues
Describe the bug This is a thread to track improvements / issues related to integration stability when using dynamic liquid layers in RadialSolver.
- Causes
- Stability problems appear when a dynamic liquid layer is present in a world with a moderately large forcing period.
- The value of the forcing period depends on the density of the liquid. A low density liquid tends to be more stable at higher forcing periods than a high density layer.
Note that "forcing period" refers to the tidal mode that is being calculated. A planet with a low orbital period may still require larger forcing period modes depending on its eccentricity, obliquity, and rotation rate.