Observatoire de Paris
<-   Calculation of the Roche limit   ->
Equilibrium or rupture, under the effect of the self-gravitational field, which must ensure cohesion, and of the gradient terms of the gravitational field, which tears the satellite apart (in the frame of reference of the center of mass of the satellite).
Copyright : ASM

Roche's argument is as follows :  although the satellite is spherical, let us suppose that it is made up of two spheres of radius r and mass m. Let us think about two snow balls of radius r, held together by the gravitation force that each sphere exerts on the other. That force, F_att, is given by Newton's relation :
F_att= Gm^(2)/((2*r))^(2)
Now, let us consider that each satellite is at a distance D of a planet of mass M and radius R. The attraction force F, between the planet and the nearest snow ball, will be greater than the force F' between the planet and the more distant ball. The force F is given by the following relation  :
And force F' is given by  :
The two balls will feel the resulting force F_mar, which tends to separate them. This force is equivalent to the difference between the forces F and F'. Therefore we have  : F_mar=F-F' And, because D>>r : F_mar = -4*G*M*m*r/D^(3)
The two masses will be separated if the force F_mar is greater than the force F_att, that is if :
2^(4)*M/D^(3) > m/r^(3)
Let us now replace the mass M by rho_P*(4/3)*pi*R^(3), where rho_Pis the density of the planet, and the mass m by rho*(4/3)*pi*r^(3), where rho is the density of the satellite.
Therefore, the two masses are separated if the distance D is less than  2^(4/3)*R * (rho_P/rho)^(1/3)
. This is a quite satisfactory approximation since 2^(4/3) is equal to 2.51, whereas the exact value is 2,456.