Articles

Limite de Roche

a distância limite para a qual um satélite pode se aproximar sem quebrar depende da rigidez do satélite. Num extremo, um satélite completamente rígido manterá a sua forma até que as forças de maré o quebrem. No outro extremo, um satélite altamente fluido gradualmente se deforma levando a um aumento das forças de maré, fazendo com que o satélite alongasse, aumentando ainda mais as forças de maré e fazendo com que se rompesse mais rapidamente.

a maioria dos satélites reais estaria em algum lugar entre estes dois extremos, com a resistência à tração tornando o satélite nem perfeitamente rígido nem perfeitamente fluido. Por exemplo, um asteroide empilhado de escombros se comportará mais como um fluido do que um sólido rochoso; um corpo gelado se comportará rigidamente no início, mas se tornará mais fluido à medida que o aquecimento das marés se acumula e seus gelos começam a derreter.

mas note que, como definido acima, o limite de Roche refere-se a um corpo mantido unido apenas pelas forças gravitacionais que causam partículas de outra forma não relacionadas a coalesce, formando assim o corpo em questão. O limite de Roche também é normalmente calculado para o caso de uma órbita circular, embora seja direto modificar o cálculo para se aplicar ao caso (por exemplo) de um corpo passando o primário em uma trajetória parabólica ou hiperbólica.

cálculo por satélite rígido

o limite de Roche de corpo rígido é um cálculo simplificado para um satélite esférico. Formas irregulares, como as de deformação de maré no corpo ou as órbitas primárias são negligenciadas. Supõe-se que esteja em equilíbrio hidrostático. Estes pressupostos, embora irrealistas, simplificam grandemente os cálculos.

O limite de Roche para uma rígida esférica satélite é a distância d {\displaystyle d}

d

, da primária, em que a força gravitacional em um teste de massa na superfície do objeto é exatamente igual à força de maré puxando a massa para longe do objeto: d = R M ( 2 ρ M ρ m ) 1 3 {\displaystyle d=R_{M}\left(2{\frac {\rho _{M}}{\rho _{m}}}\right)^{\frac {1}{3}}}

{\displaystyle d=R_{M}\left(2{\frac {\rho _{M}}{\rho _{m}}}\right)^{\frac {1}{3}}}

, onde R M {\displaystyle R_{M}}

R_M

é o raio do primário, ρ M {\displaystyle \rho _{M}}

\rho_M

é a densidade do primário, e ρ m {\displaystyle \rho _{m}}

\rho_m

é a densidade do satélite. Isso pode ser equivalentemente escrita como d = R m ( 2 M M M m ) 1 3 {\displaystyle d=R_{m}\left(2{\frac {M_{M}}{M_{m}}}\right)^{\frac {1}{3}}}

{\displaystyle d=R_{m}\left(2{\frac {M_{M}}{M_{m}}}\right)^{\frac {1}{3}}}

, onde R m {\displaystyle R_{m}}

R_m

é o raio do secundário, M M {\displaystyle M_{M}}

M_M

é a massa do primário, e M m {\displaystyle M_{m}}

M_{m}

é a massa do secundário.

isto não depende do tamanho dos objetos, mas da razão de densidades. Esta é a distância orbital dentro da qual material solto (por exemplo, regolito) na superfície do satélite mais próximo do primário seria puxado para longe, e do mesmo modo o material do lado oposto ao primário também irá para longe, em vez de para, o satélite.

Note que este é um resultado aproximado como força de inércia e estrutura rígida são ignorados em sua derivação.

O período orbital depende então apenas da densidade do sistema:

P = 2 π d 3 G M M ) 1 / 2 = 2 π ( d, 3 ( 4 / 3 ) π G R M 3 ρ M ) 1 / 2 = 6 π G ρ m {\displaystyle P=2\pi \left({\frac {d^{3}}{GM_{M}}}\right)^{1/2}=2\pi \left({\frac {d^{3}}{(4/3)\pi GR_{M}^{3}\rho _{M}}}\right)^{1/2}={\sqrt {\frac {6\pi }{G\rho _{m}}}}}

{\displaystyle P=2\pi \left({\frac {d^{3}}{GM_{M}}}\right)^{1/2}=2\pi \left({\frac {d^{3}}{(4/3)\pi GR_{M}^{3}\rho _{M}}}\right)^{1/2}={\sqrt {\frac {6\pi }{G\rho _{m}}}}}

, onde G é a constante gravitacional. Por exemplo, uma densidade de 3.346 g / cc (a densidade da nossa lua) corresponde a um período orbital de 2.552 horas.

Derivação da formulaEdit

Derivação do limite de Roche

a fim de determinar o limite de Roche, considere uma pequena massa de u {\displaystyle u}

u

na superfície do satélite mais próximo do primário. Existem duas forças nesta massa u {\displaystyle u}

u

: a atração gravitacional em direção ao satélite e a atração gravitacional em direção ao primário. Suponha que o satélite está em queda livre em torno do primário e que a força de maré é o único termo relevante da atração gravitacional do primário. Esta suposição é uma simplificação, pois a queda livre só se aplica verdadeiramente ao centro planetário, mas será suficiente para esta derivação.

A força gravitacional F G {\displaystyle F_{\text{G}}}

F_{{\text{G}}}

sobre a massa u {\displaystyle u}

u

para o satélite de massa m {\displaystyle m}

m

e raio r {\displaystyle r}

r

pode ser expressa de acordo com a lei de Newton da gravitação. F G = G m u r 2 {\displaystyle F_{\text{G}}={\frac {Gmu}{r^{2}}}}

F_{{\text{G}}}={\frac {Gmu}{r^{2}}}

a força de maré F T {\displaystyle F_{\text{T}}}

F_{{\text{T}}}

sobre a massa u {\displaystyle u}

u

para o primário, com raio R {\displaystyle R}

R

e massa M {\displaystyle M}

M

, a uma distância d {\displaystyle d}

d

entre os centros dos dois corpos, pode ser expressa, aproximadamente, como F T = 2 G M u d r o 3 {\displaystyle F_{\text{T}}={\frac {2GMur}{d^{3}}}}

F_{{\text{T}}}={\frac {2GMur}{d^{3}}}

. para obter esta aproximação, encontre a diferença na força gravitacional do primário no centro do satélite e na borda do satélite mais próxima do primário.: F T = G M u ( d − r ) 2 − G M u d 2 {\displaystyle F_{\text{T}}={\frac {GMu}{(d-r)^{2}}}-{\frac {GMu}{d^{2}}}}

F_{{\text{T}}}={\frac {GMu}{(d-r)^{2}}}-{\frac {GMu}{d^{2}}}

F T = G M u d 2 − ( d − r ) 2 d 2 ( d − r ) 2 {\displaystyle F_{\text{T}}=GMu{\frac {d^{2}-(d-r)^{2}}{d^{2}(d-r)^{2}}}}

F_{{\text{T}}}=GMu{\frac {d^{2}-(d-r)^{2}}{d^{2}(d-r)^{2}}}

F T = G M u 2 d r − r 2 d 4 − 2 d 3 r + r 2 d 2 {\displaystyle F_{\text{T}}=GMu{\frac {2dr-r^{2}}{d^{4}-2d^{3}r+r^{2}d^{2}}}}

F_{{\text{T}}}=GMu{\frac {2dr-r^{2}}{d^{4}-2d^{3}r+r^{2}d^{2}}}

Na aproximação onde r ≪ R {\displaystyle r\ll R}

r\ll R

e R < d {\displaystyle R<d}

Rd

, pode-se dizer que o r 2 {\displaystyle r^{2}}

r^{2}

no numerador e a cada período com r {\displaystyle r}

r

no denominador vai para zero, o que nos dá: F T = G M u 2 d r d 4 {\displaystyle F_{\text{T}}=GMu{\frac {2dr}{d^{4}}}}

F_{{\text{T}}}=GMu{\frac {2dr}{d^{4}}}

F T = 2 G M u r d 3 {\displaystyle F_{\text{T}}={\frac {2GMur}{d^{3}}}}

F_{{\text{T}}}={\frac {2GMur}{d^{3}}}

The Roche limit is reached when the gravitational force and the tidal force balance each other out.

F G = F T {\displaystyle F_{\text{G}}=F_{\text{T}}\;}

F_{{\text{G}}}=F_{{\text{T}}}\;

ou

G m u r 2 = 2 G M u d r o 3 {\displaystyle {\frac {Gmu}{r^{2}}}={\frac {2GMur}{d^{3}}}}

\frac{Gmu}{r^2} = \frac{2GMur}{d^3}

,

o que dá o limite de Roche, d {\displaystyle d}

d

, como d = r ( 2 M ) 1 3 {\displaystyle d=r\left(2\,{\frac {M}{m}}\right)^{\frac {1}{3}}}

d=r\left(2\,{\frac {M}{m}}\right)^{{{\frac {1}{3}}}}

O raio do satélite não deverá aparecer na expressão para o limite, então é re-escrito em termos de densidades.

Para uma esfera de massa M {\displaystyle M}

M

pode ser escrito como M = 4 π ρ M R 3 3 {\displaystyle M={\frac {4\pi \rho _{M}R^{3}}{3}}}

M = \frac{4\pi\rho_M R^3}{3}

onde R {\displaystyle R}

R

é o raio do primário.

E também

m = 4 π ρ m r 3 3 {\displaystyle m={\frac {4\pi \rho _{m}r^{3}}{3}}}

m = \frac{4\pi\rho_m r^3}{3}

onde r {\displaystyle r}

r

é o raio do satélite.

Substituindo para as massas, a equação para o limite de Roche, e cancelando os 4 π / 3 {\displaystyle 4\pi /3}

4\pi /3

dá-d = r ( 2 ρ M R 3 ρ m r 3 ) 1 / 3 {\displaystyle d=r\left({\frac {2\rho _{M}R^{3}}{\rho _{m}r^{3}}}\right)^{1/3}}

d = r \left( \frac{ 2 \rho_M R^3 }{ \rho_m r^3 } \right)^{1/3}

,

o que pode ser simplificado para o seguinte limite de Roche:

d = R ( 2 ρ M ρ m ) 1 3 ≈ 1.26 R ( ρ M ρ m ) 1 3 {\displaystyle d=R\left(2\,{\frac {\rho _{M}}{\rho _{m}}}\right)^{\frac {1}{3}}\approx de 1,26 R\left({\frac {\rho _{M}}{\rho _{m}}}\right)^{\frac {1}{3}}}

d=R\left(2\,{\frac {\rho _{M}}{\rho _{m}}}\right)^{{{\frac {1}{3}}}}\cerca de 1,26 R\left({\frac {\rho _{M}}{\rho _{m}}}\right)^{{{\frac {1}{3}}}}

. uma fórmula mais precisa

Uma vez que um satélite próximo provavelmente estará orbitando em uma órbita quase circular com rotação síncrona, considere como a força centrífuga da rotação afetará os resultados. Que força é

F C = ω 2 u r = G M u d r o 3 {\displaystyle F_{C}=\omega ^{2}ur={\frac {GMur}{d^{3}}}}

F_C = \omega^2 ur = \frac{GMur}{d^3}

e é adicionada a FT. Fazendo a força, o equilíbrio cálculo produz esse resultado para o limite de Roche:

d = R M ( 3 ρ M ρ m ) 1 3 ≈ 1.442 R M ( ρ M ρ m ) 1 3 {\displaystyle d=R_{M}\left(3\;{\frac {\rho _{M}}{\rho _{m}}}\right)^{\frac {1}{3}}\approx 1.442R_{M}\left({\frac {\rho _{M}}{\rho _{m}}}\right)^{\frac {1}{3}}}

d=R_{M}\left(3\;{\frac {\rho _{M}}{\rho _{m}}}\right)^{{{\frac {1}{3}}}}\cerca de 1.442R_{M}\left({\frac {\rho _{M}}{\rho _{m}}}\right)^{{{\frac {1}{3}}}}

………. (1)

ou: d = R m ( 3 M M M m ) 1 3 ≈ 1.442 R m ( M M ) 1 3 {\displaystyle d=R_{m}\left(3\;{\frac {M_{M}}{M_{m}}}\right)^{\frac {1}{3}}\approx 1.442R_{m}\left({\frac {M_{M}}{M_{m}}}\right)^{\frac {1}{3}}}

{\displaystyle d=R_{m}\left(3\;{\frac {M_{M}}{M_{m}}}\right)^{\frac {1}{3}}\approx 1.442R_{m}\left({\frac {M_{M}}{M_{m}}}\right)^{\frac {1}{3}}}

………. (2)

Use m = 4 π ρ m r 3 3 {\displaystyle m={\frac {4\pi \rho _{m}r^{3}}{3}}}

m = \frac{4\pi\rho_m r^3}{3}

(onde r {\displaystyle r}

r

é o raio do satélite) para substituir ρ m {\displaystyle \rho _{m}}

\rho_m

na fórmula(1), podemos ter uma terceira fórmula:
d = ( 9 M M 4 π ρ m ) 1 3 ≈ 0.8947 ( M M ρ m ) 1 3 {\displaystyle d=\left({\frac {9M_{M}}{4\pi \rho _{m}}}\right)^{\frac {1}{3}}\approx 0.8947\left({\frac {M_{M}}{\rho _{m}}}\right)^{\frac {1}{3}}}

d=\left({\frac {9M_{M}}{4\pi \rho _{m}}}\right)^{{{\frac {1}{3}}}}\cerca de 0.8947\left({\frac {M_{M}}{\rho _{m}}}\right)^{{{\frac {1}{3}}}}

………. Assim, é suficiente observar a massa da estrela (planeta) e estimar a densidade do planeta (satélite) para calcular o limite de Roche do planeta (satélite) no sistema estelar (planetário).

limite de Roche, Hill e esfera de raio da planetEdit

a Comparação da Colina esferas e Roche limites do Sol-sistema Terra-Lua (sem escala) com as regiões sombreadas denotando estável órbitas dos satélites de cada corpo

Considere a possibilidade de um planeta com uma densidade de ρ m {\displaystyle \rho _{m}}

\rho_m

e um raio de r {\displaystyle r}

r

, em órbita de uma estrela withis é o significado físico do limite de Roche, lóbulo de Roche e Hill sphere.

Fórmula(2) pode ser descrito como: R Roche = R Hill 3 M m 3 = R secundária 3 M m 3 {\displaystyle R_{\text{Roche}}=R_{\text{Hill}}{\sqrt{\frac {3 M}{m}}}=R_{\text{secundário}}{\sqrt{\frac {3 M}{m}}}}

R_{{{\text{Roche}}}}=R_{{{\text{Hill}}}}{\sqrt{{\frac {3 M}{m}}}}=R_{{\text{secundário}}}{\sqrt{{\frac {3 M}{m}}}}

, um matemático perfeito simetria.esta é a significância astronômica do limite de Roche e da esfera de Hill. Nota: O limite de Roche e a esfera de Hill são completamente diferentes uns dos outros, mas são ambos trabalhos de Édouard Roche.a esfera de Colina de um corpo astronômico é a região na qual domina a atração de satélites, enquanto o limite de Roche é a distância mínima para a qual um satélite pode se aproximar de seu corpo primário sem força de maré superando a gravidade interna que mantém o satélite Unido.

Nota : limite de Roche e esfera de Hill são completamente diferentes um do outro, mas são ambos trabalho de Édouard Roche.esfera de Colina de um corpo astronômico é a região em que domina a atração de satélites, enquanto limite de Roche é a distância mínima para a qual um satélite pode se aproximar de seu corpo primário sem força de maré superando a gravidade interna que mantém o satélite Unido.

Fluid satellitesEdit

A more accurate approach for calculating the Roche limit takes the deformation of the satellite into account. Um exemplo extremo seria um satélite líquido que orbitasse um planeta, onde qualquer força que atuasse sobre o satélite iria deformá-lo em um esferóide prolato.

o cálculo é complexo e seu resultado não pode ser representado numa fórmula algébrica exata. Roche-se derivou a seguinte solução aproximada para o limite de Roche:

d ≈ 2.44 R ( ρ M ρ m ) 1 / 3 {\displaystyle d\approx de 2,44 R\left({\frac {\rho _{M}}{\rho _{m}}}\right)^{1/3}}

d \a cerca de 2.44R\left( \frac {\rho_M} {\rho_m} \right)^{1/3}

no Entanto, uma melhor aproximação que leva em conta o principal oblateness e o satélite de massa é:

d ≈ 2.423 R ( ρ M ρ m ) 1 / 3 ( ( 1 + m 3 M ) + c 3 R ( 1 + m ) 1 − c / R ) 1 / 3 {\displaystyle d\approx 2.423 R\left({\frac {\rho _{M}}{\rho _{m}}}\right)^{1/3}\left({\frac {(1+{\frac {m}{3M}})+{\frac {c}{3}}(1+{\frac {m}{M}})}{1-c/R}}\right)^{1/3}}

d \a cerca de 2.423 R\left( \frac {\rho_M} {\rho_m} \right)^{1/3} \left( \frac{(1+\frac{m}{3M})+\frac{c}{3}(1+\frac{m}{M})}{1-c/R} \right)^{1/3}

em que c / R {\displaystyle c/R}

c/R

é o oblateness do primário. O factor numérico é calculado com a ajuda de um computador. a solução de fluido é apropriada para corpos que são apenas vagamente mantidos juntos, como um cometa. Por exemplo, a órbita em decomposição do cometa Shoemaker–Levy 9 em torno de Júpiter passou dentro de seu limite de Roche em julho de 1992, fazendo com que se fragmentasse em um número de peças menores. Em sua próxima aproximação em 1994, os fragmentos caíram no planeta. Shoemaker-Levy 9 foi observado pela primeira vez em 1993, mas sua órbita indicou que havia sido capturada por Júpiter algumas décadas antes.como o caso do satélite fluido é mais delicado do que o rígido, o satélite é descrito com algumas suposições simplificadoras. Primeiro, suponha que o objeto consiste incompressível de fluidos que tem densidade constante ρ m {\displaystyle \rho _{m}}

\rho_m

e o volume V {\displaystyle V}

V

que não dependem externa ou interna de forças. em segundo lugar, assuma que o satélite se move em uma órbita circular e ele permanece em rotação síncrona. Isto significa que a velocidade angular ω {\displaystyle \omega }

\omega

em que gira em torno do seu centro de massa é a mesma que a velocidade angular em que se move em torno do baricentro do sistema global.

a velocidade angular ω {\displaystyle \omega }

\omega

é dada pela terceira lei de Kepler: ω 2 = G m + m D 3 . {\displaystyle \omega ^{2}=G\, {\frac {M+m}{D^{3}}}}.}

\omega^2 = G \, \frac{m + m}{D^3}.

quando M é muito maior que m, esta será próxima de

ω 2 = G M D 3 . {\displaystyle \omega ^{2}=G\, {\frac {m}{D^{3}}}}.}

\omega^2 = G \, \frac{m}{D^3}.

a rotação síncrona implica que o líquido não se move e o problema pode ser considerado estático. Portanto, a viscosidade e o atrito do líquido neste modelo não desempenham um papel, uma vez que estas quantidades desempenhariam um papel apenas para um fluido em movimento.dados estes pressupostos, devem ser tidas em conta as seguintes forças::a força gravitacional devida ao corpo principal; a força centrífuga no sistema rotativo de referência; e o campo de auto-gravitação do satélite.

Uma vez que todas estas forças são conservadoras, elas podem ser expressas por meio de um potencial. Além disso, a superfície do satélite é equipotencial. Caso contrário, as diferenças de potencial daria origem a forças e movimentos de algumas partes do líquido na superfície, o que contradiz a suposição do modelo estático. Dada a distância do corpo principal, deve determinar-se a forma da superfície que satisfaz a condição equipotencial.

distância Radial de um ponto sobre a superfície do elipsóide para o centro de massa

Como a órbita tem sido assumido circular, o total da força gravitacional e orbital força centrífuga agindo sobre o corpo principal cancelar. Isso deixa duas forças: a força de maré e a força centrífuga de rotação. A força de maré depende da posição em relação ao centro de massa, já considerado no modelo rígido. Para corpos pequenos, a distância das partículas líquidas do centro do corpo é pequena em relação à distância d para o corpo principal. Assim, a força de maré pode ser linearizada, resultando na mesma fórmula para FT como dado acima.enquanto esta força no modelo rígido depende apenas do raio r do satélite, no caso do fluido, todos os pontos da superfície devem ser considerados, e a força de maré depende da distância Δd do centro da massa a uma dada partícula projetada na linha que une o satélite ao corpo principal. Chamamos Δd de distância radial. Desde que a força de maré é linear, Δd, o potencial é proporcional ao quadrado da variável e para m ≪ M {\displaystyle m\ll M}

m\ll M

temos V T = − 3 G M 2 d 3 D d 2 {\displaystyle V_{T}=-{\frac {3GM}{2d^{3}}}\Delta d^{2}\,}

V_T = - \frac{3 G M }{2 d^3}\Delta d^2 \,

da Mesma forma, a força centrífuga tem um potencial

V C = − 1 2 ω 2 ∆ d 2 = − G M 2 d 3 D d 2 {\displaystyle V_{C}=-{\frac {1}{2}}\omega ^{2}\Delta d^{2}=-{\frac {GM}{2d^{3}}}\Delta d^{2}\,}

V_C = - \frac{1}{2} \omega^2 \Delta d^2 = - \frac{G }{2 d^3}\Delta d^2 \,

para velocidade angular de rotação ω {\displaystyle \omega }

\omega

.

queremos determinar a forma do satélite para o qual a soma do potencial de auto-gravitação e VT + VC é constante na superfície do corpo. Em geral, tal problema é muito difícil de resolver, mas neste caso particular, pode ser resolvido por um palpite hábil devido à dependência quadrada do potencial de maré na distância radial Δd para uma primeira aproximação, podemos ignorar o potencial centrífugo VC e considerar apenas o potencial de maré VT.uma vez que o potencial VT muda apenas em uma direção, ou seja, a direção em direção ao corpo principal, o satélite pode ser esperado para tomar uma forma axial simétrica. Mais precisamente, podemos partir do princípio de que se trata de um sólido de revolução. O auto-potencial na superfície de tal sólido de revolução só pode depender da distância radial para o centro de massa. Na verdade, a intersecção do satélite com um plano perpendicular à linha que une os corpos é um disco cujo limite, pelas nossas suposições, é um círculo de potencial constante. Se a diferença entre o potencial de auto-gravitação e o VT for constante, ambos os potenciais devem depender da mesma forma de Δd. Em outras palavras, o auto-potencial tem que ser proporcional ao quadrado de Δd. Então pode-se mostrar que a solução equipotencial é um elipsoide de revolução. Dada uma densidade e volume constantes, o auto-potencial de tal corpo depende apenas da excentricidade ε do elipsóide.:

V s = V s 0 + G π ρ m ⋅ f ( ε ) ⋅ ∆ d 2 , {\displaystyle V_{s}=V_{s_{0}}+G\pi \rho _{m}\cdot f(\varepsilon )\cdot \Delta d^{2},}

{\displaystyle V_{s}=V_{s_{0}}+G\pi \rho _{m}\cdot f(\varepsilon )\cdot \Delta d^{2},}

, onde V s 0 {\displaystyle V_{s_{0}}}

V_{s_0}

é a constante de auto-potencial na intersecção da aresta circular do corpo e a simetria central avião dada pela equação Δd=0.

O adimensional função f está para ser determinada a partir do precisas solução para o potencial do elipsóide

f ( ε ) = 1 − ε 2 ε 3 ⋅ {\displaystyle f(\varepsilon )={\frac {1-\varepsilon ^{2}}{\varepsilon ^{3}}}\cdot \left}

{\displaystyle f(\varepsilon )={\frac {1-\varepsilon ^{2}}{\varepsilon ^{3}}}\cdot \left}

e, surpreendentemente, não depende do volume do satélite.

O gráfico de uma função adimensional f, que indica como a força das marés de potencial depende da excentricidade ε do elipsóide.

embora a forma explícita da função f pareça complicada, é claro que podemos e podemos escolher o valor de ε Para que o potencial VT seja igual a VS mais uma constante independente da variável Δd. Por inspeção, isso ocorre quando a

2 G π ρ M R 3 d 3 = G π ρ m f ( ε ) {\displaystyle {\frac {2G\pi \rho _{M}R^{3}}{d^{3}}}=G\pi \rho _{m}f(\varepsilon )}

{\displaystyle {\frac {2G\pi \rho _{M}R^{3}}{d^{3}}}=G\pi \rho _{m}f(\varepsilon )}

Esta equação pode ser resolvida numericamente. O grafo indica que existem duas soluções e, assim, a menor representa a forma de equilíbrio estável (o elipsoide com a excentricidade menor). Esta solução determina a excentricidade do elipsoide de maré em função da distância ao corpo principal. A derivada da função f tem um zero onde a excentricidade máxima é atingida. Isto corresponde ao limite de Roche.

o derivado de f determina a excentricidade máxima. Isto dá o limite de Roche.

Mais precisamente, o limite de Roche é determinado pelo fato de que a função f, que pode ser considerada como uma medida não linear da força comprimindo o elipsoide em direção a uma forma esférica, é limitada de modo que há uma excentricidade na qual esta força de contratação se torna máxima. Uma vez que a força de maré aumenta quando o satélite se aproxima do corpo principal, é claro que há uma distância crítica em que o elipsoide é rasgado.

a excentricidade máxima pode ser calculada numericamente como o zero da derivada de f’. Obtém-se

ε max ≈ 0 . 86 {\\displaystyle \varepsilon _{\text {max}}\approx 0{.} 86}

{\displaystyle \varepsilon _{\text{max}} \ approx 0{.}86}

que corresponde à razão dos eixos elipsóides 1:1, 95. Inserindo isto na fórmula para a função f pode-se determinar a distância mínima a que o elipsoide existe. Este é o limite de Roche,

D ≈ 2 . 423 R r ρ ρ M ρ m 3 . {\displaystyle D \ approx 2{.}423\cdot R\cdot {\sqrt{\frac {\rho _{m}} {\rho _{m}}}}}}}}\,.

d \ approx 2{.}423 \cdot R \cdot \sqrt{ \frac {\rho_M} {\rho_m} } \,.

surpreendentemente, incluindo o potencial centrífugo faz notavelmente pouca diferença, embora o objeto se torne um elipsoide de Roche, um elipsoide triaxial geral com todos os eixos com comprimentos diferentes. O potencial torna-se uma função muito mais complicada dos comprimentos dos eixos, exigindo funções elípticas. No entanto, a solução prossegue tanto quanto no caso de maré-apenas, e encontramos

D ≈ 2 . 455 R r ρ ρ M ρ m 3 . {\displaystyle D \ approx 2{.}455\cdot R\cdot {\sqrt{\frac {\rho _{m}} {\rho _{m}}}}}}}}\,.