Abstract
The Laplace Tidal Equations (LTEs) play a basic role for investigating the dissipation of long-term orbital energy as heat in surface oceans of icy satellites. Here we address solving LTEs by means of finite elements using PDE2D, a general-purpose finite element program for solving multidimensional PDEs. We use PDE2D to solve the equations for eccentricity and obliquity gravitational potentials for a shallow water approximation, and calculate ocean dissipation in three icy satellites, Enceladus, Europa and Titan. Compared with previous modeling procedures, PDE2D solves for unknowns that are approximated by piecewise cubic polynomials, thus attaining a higher O(dφ4) + O(dθ4) numerical precision. Our modeling results show, within few percent errors, a good agreement with earlier semi-analytical models and numerical solutions. We provide an open access to our PDE2D codes which are intended to be simple and easy to adjust for a wide range of tidal heating related problems. Our modeling results show that PDE2D can be a useful tool for researchers studying tidal heating related problems.
Sign-in/Register to access full text options 
Free) 
Talk to us
Join us for a 30 min session where you can share your feedback and ask us any queries you have
Schedule a call
