Tidal disruption of dissipative planetesimals

A self-consistent numerical model is developed for the tidal disruption of a solid planetesimal. The planetesimal is treated as a highly viscous, slightly compressible fluid whose disturbed parts are an inviscid, pressureless fluid undergoing distortion and disruption. The distortions were constrained to being symmetrical above and below the equatorial plane. The tidal potential is expanded in terms of Legendre polynomials, which eliminates the center of mass acceleration effects, permitting definition of equations of motion in a noninertial frame. Consideration is given to viscous dissipation and to characteristics of the solid-atmosphere boundary. The model is applied to sample cases in one, two and three dimensions.