Vamos a estudiar un modelo de Tierra que consta de un núcleo de radio r_c=3490 km y densidad constante ρ_c=11.0 g/cm³ y un manto de radio interior r_c y radio exterior R=6371 km de densidad constante ρ_m=4.437 g/cm³.

Naturalmente, ρ_c, ρ_m y la densidad media ρ=5.517 g/cm³ están relacionadas mediante la ecuación

${\displaystyle \rho R^3={\rho }_c{r}_c^3+{\rho }_m(R^3-r_c^3)}$

Variación de la aceleración de la gravedad g con el radio r

En la página titulada “Aceleración de la gravedad” demostramos que la aceleración de la gravedad g en un punto P en el interior de una distribución esférica de masa, se debe solamente a la masa contenida en la esfera de radio r<R, y se calcula mediante la fórmula

${\displaystyle g={\displaystyle \underset{0}{\overset{r}{\int }}\frac{G\rho (x)4\pi x^2}{r^2}}dx}$

• La aceleración de la gravedad en un punto situado en el interior del núcleo r<r_c

${\displaystyle g={\displaystyle \underset{0}{\overset{r}{\int }}\frac{G{\rho }_{c}4\pi x^2}{r^2}dx}=\frac{4}{3}\pi G{\rho }_{c}rr<r_c}$

g=19.577 (r/R) m/s²

• La aceleración de la gravedad en un punto situado en el interior del manto r_c<r<R

${\displaystyle \begin{array}{l}g={\displaystyle \underset{0}{\overset{r_c}{\int }}\frac{G{\rho }_{c}4\pi x^{2}dx}{r^2}}+{\displaystyle \underset{r_c}{\overset{r}{\int }}\frac{G{\rho }_{m}4\pi x^{2}dx}{r^2}}=\frac{4\pi }{r^2}G{\rho }_c\frac{r_c^3}{3}+\frac{4\pi }{r^2}G{\rho }_m\frac{r^3-r_c^3}{3}=\\ g=\frac{4\pi }{3}G\left\{({\rho }_c-{\rho }_m)\frac{r_c^3}{r^2}+{\rho }_{m}r\right\}{r}_c<r<R\end{array}}$

g=1.920(R/r)²+7.898(r/R) m/s²

• La aceleración de la gravedad en un punto situado en el exterior de la Tierra r>R

${\displaystyle g=G\frac{M}{r^2}=\frac{4\pi }{3}G\rho \frac{R^3}{r^2}}$

g=9.82(R/r)² m/s²

En la figura, se muestra la dependencia del módulo de g con el cociente r/R (en color rojo) y se compara con el valor de g en función de (r/R) suponiendo la Tierra homogénea (en color azul), de densidad constante e igual a la densidad media

La aceleración de la gravedad crece linealmente con el radio r, en el núcleo para r<r_c. En el manto g varía relativamente poco, disminuyendo al principio y aumentando al final, tal como vemos con más detalle en la figura inferior.

Calculamos el valor mínimo de la aceleración de la gravedad en el manto, derivando g con respecto de r e igualando a cero.

$$\begin{gathered} {\displaystyle \frac{dg}{dr}=g=\frac{4\pi }{3}G\left\{-2({\rho }_c-{\rho }_m)\frac{r_c^3}{r^3}+{\rho }_m\right\}=0 \\ {r}_{\min }=r_c\sqrt[3]{\frac{2({\rho }_c-{\rho }_m)}{{\rho }_m}}=0.786·R} \end{gathered}$$

En una primera aproximación, tomamos la aceleración de la gravedad en el manto como constante e igual al siguiente valor medio (recta de color azul).

${\displaystyle g_m=\frac{g(r_{\min })+g(R)}{2}=\frac{9.32+9.82}{2}=9.57{\text{m/s}}^{\text{2}}}$

Energía potencial

Una partícula de masa m experimenta una fuerza F=mg dirigida radialmente hacia el centro de la Tierra. La energía potencial E_p(r)  correspondiente a esta fuerza conservativa se calcula del siguiente modo

${\displaystyle \begin{array}{l}F=\frac{-d{E}_p}{dr}\\ {\displaystyle \underset{r}{\overset{\infty }{\int }}d{E}_p}=-{\displaystyle \underset{r}{\overset{\infty }{\int }}F(r)·dr}\\ E_p(r)={\displaystyle \underset{r}{\overset{\infty }{\int }}F(r)·dr}\end{array}}$

Tomamos como nivel cero de energía potencial E_p(∞)=0

• Energía potencial de la partícula para r>R

Para r>R, la expresión de la fuerza sobre la partícula es

${\displaystyle F(r)=-mG\frac{4\pi }{3}\rho \frac{R^3}{r^2}}$

La expresión de la energía potencial en este intervalo es

${\displaystyle E_p(r)={\displaystyle \underset{r}{\overset{\infty }{\int }}-G\frac{4\pi }{3}\rho R^3\frac{m}{r^2}dr=}-\frac{4\pi }{3}mG\rho \frac{R^3}{r}}$

• Energía potencial de la partícula para r_c<r<R

La expresión de la fuerza sobre la partícula en esta región es

${\displaystyle F(r)=-\frac{4\pi }{3}mG\left\{({\rho }_c-{\rho }_m)\frac{r_c^3}{r^2}+{\rho }_{m}r\right\}}$

La expresión de la energía potencial en este intervalo es

${\displaystyle \begin{array}{l}{E}_p(r)={\displaystyle \underset{r}{\overset{\infty }{\int }}F(r)dr=}{\displaystyle \underset{r}{\overset{R}{\int }}-\frac{4\pi }{3}mG\left\{({\rho }_c-{\rho }_m)\frac{r_c^3}{r^2}+{\rho }_{m}r\right\}dr+}{\displaystyle \underset{R}{\overset{\infty }{\int }}-G\frac{4\pi }{3}\rho R^3\frac{m}{r^2}}\\ =-\frac{4\pi }{3}mG\left(-\frac{1}{2}{\rho }_m{r}^2+({\rho }_c-{\rho }_m)\frac{r_c^3}{r}+\frac{3}{2}{\rho }_m{R}^2\right)\end{array}}$

Para obtener esta expresión, se ha empleado la relación entre las densidades del manto ρ_m, núcleo ρ_c y media ρ.

${\displaystyle \rho R^3={\rho }_c{r}_c^3+{\rho }_m(R^3-r_c^3)}$

• Energía potencial de la partícula para 0<r<r_c

La expresión de la fuerza sobre la partícula en esta región es

${\displaystyle F(r)=-\frac{4\pi }{3}mG{\rho }_{c}r}$

La expresión de la energía potencial en este intervalo es

${\displaystyle \begin{array}{l}{E}_p(r)=\underset{r}{\overset{\infty }{{\displaystyle \int }}}F(r)dr=\underset{r}{\overset{r_c}{{\displaystyle \int }}}\left(-\frac{4\pi }{3}mG{\rho }_{c}r\right)dr+\underset{r_c}{\overset{R}{{\displaystyle \int }}}-\frac{4\pi }{3}mG\left(({\rho }_c-{\rho }_m)\frac{r_c^3}{r^2}+{\rho }_{m}r\right)dr\\ +\text{}\underset{R}{\overset{\infty }{{\displaystyle \int }}}-G\frac{4\pi }{3}\rho R^3\frac{m}{r^2}=-\frac{4\pi }{3}mG\left(-\frac{1}{2}{\rho }_c{r}^2+\frac{3}{2}({\rho }_c-{\rho }_m)r_c^2+\frac{3}{2}{\rho }_m{R}^2\right)\end{array}}$

En la figura, se representa la energía potencial E_p(r) (en color rojo) en función del cociente (r/R), y se compara con la energía potencial correspondiente al modelo de Tierra homogénea (en color azul).

El principio de conservación de la energía

${\displaystyle \frac{1}{2}m{v}^2+E_p(r)=\text{cte}}$

Conocida la posición r₀ inicial y la velocidad inicial v₀ de la partícula, calculamos la velocidad v de la partícula para una determinada posición r.

Para un determinado valor de la energía total, una partícula que atraviese la Tierra a través de uno de sus diámetros alcanzará una velocidad mayor, para una determinada distancia al centro, en el modelo de dos capas, que en el modelo de Tierra homogénea.

Presión en el centro de la Tierra

La ecuación fundamental de la hidrostática es

${\displaystyle dp=-\rho gdr}$

En la página titulada "La aceleración de la gravedad" calculamos la presión en el centro de la Tierra suponiendo que la densidad de la Tierra es constante e igual a la densidad media. La presión en el modelo más realista de núcleo más manto ambos de densidad constante es

${\displaystyle \begin{array}{l}{p}_0-p=-\underset{0}{\overset{R}{{\displaystyle \int }}}\rho (r)·g(r)dr=-\underset{r_c}{\overset{R}{{\displaystyle \int }}}{\rho }_m\frac{4\pi }{3}G\left(({\rho }_m-{\rho }_c)\frac{r_c^3}{r^2}-{\rho }_{m}r\right)dr\\ -\text{}\underset{0}{\overset{r_c}{{\displaystyle \int }}}{\rho }_c\frac{4\pi }{3}G{\rho }_{c}r·dr=\frac{4\pi }{3}G\left\{\left(\frac{3}{2}{\rho }_m^2-{\rho }_m{\rho }_c-\frac{1}{2}{\rho }_c^2\right)r_c^2+({\rho }_m{\rho }_c-{\rho }_m^2)\frac{r_c^3}{R}-{\rho }_m^2\frac{R^2}{2}\right\}\end{array}}$

donde p₀=1.013·10⁵ Pa es la presión para r=R es decir, la presión atmosférica, que es despreciable frente a la presión creada en el centro por el manto y el núcleo.

p=3.29·10¹¹ Pa

Aproximación

Si consideramos la aceleración de la gravedad constante en el manto e igual a su valor medio g_m. La expresión para la presión p en el centro de la Tierra es mucho más simple.

${\displaystyle \begin{array}{l}{p}_0-p=-{\displaystyle \underset{0}{\overset{R}{\int }}\rho (r)·g(r)dr=}-{\displaystyle \underset{r_c}{\overset{R}{\int }}{\rho }_m{g}_{m}dr}-{\displaystyle \underset{0}{\overset{r_c}{\int }}{\rho }_c\frac{4\pi }{3}G{\rho }_{c}r·dr}=\\ -{\rho }_m{g}_m(R-r_c)-\frac{2\pi }{3}G{\rho }_c^2{r}_c^2\end{array}}$

p=3.28·10¹¹ Pa

Momento de inercia

El momento de inercia es un parámetro muy importante en el estudio del movimiento de la Tierra.

Consideramos la distribución esférica de masa M y radio R, dividida en capas esféricas de radio x y de espesor dx. Cada capa esférica a su vez, la dividimos en anillos de radio variable x·sinθ. La masa contenida en el anillo es dm=ρ·(2πxsinθ)·(x·dθ)·dx.

El momento de inercia del anillo respecto del eje de rotación Z es

(xsinθ)²·dm

Integramos con respecto del ángulo θ para calcular el momento de inercia de la capa esférica de radio x y de espesor dx, respecto del eje de rotación Z.

${\displaystyle dI={\displaystyle \underset{0}{\overset{\pi }{\int }}{(x·\sin \theta )}^2·dm}=2\pi \rho x^4·dx{\displaystyle \underset{0}{\overset{\pi }{\int }}{\sin }^3\theta }d\theta =\frac{8\pi }{3}\rho x^4·dx}$

El momento de inercia del núcleo con respeto al eje Z es

${\displaystyle I_c={\displaystyle \underset{0}{\overset{r_c}{\int }}\frac{8\pi }{3}{\rho }_c{x}^{4}dx}=\frac{8\pi }{15}{\rho }_c{r}_c^5}$

El momento de inercia del manto respecto del eje Z es

${\displaystyle I_m={\displaystyle \underset{r_c}{\overset{R}{\int }}\frac{8\pi }{3}{\rho }_m{x}^{4}dx}=\frac{8\pi }{15}{\rho }_m(R^5-r_c^5)}$

El momento de inercia de la Tierra respecto del eje Z es

I=I_c+I_m=9.543·10³⁶+7.412·10³⁷=8.366·10³⁷ kgm²=0.345·MR²

Siendo M=5.98·10²⁴ kg la masa de la Tierra, y R=6.371·10⁶ m el radio de la Tierra

El momento de inercia de una distribución esférica y homogénea de masa es I=2MR²/5=0.40MR²

Túnel que atraviesa la Tierra por un diámetro

Se construye un túnel que atraviesa la Tierra por uno de sus diámetros. Vamos a estudiar el movimiento de una partícula de masa m que se suelta en un extremo del túnel. Calcularemos el tiempo que tarda en completar una oscilación.

La fuerza que ejerce la Tierra sobre una partícula de masa m es $m\overrightarrow{g}$ , dirigida radialmente hacia el centro de la Tierra, su aceleración es $\overrightarrow{g}$ .

• La ecuación del movimiento de la partícula cuando se encuentra en el manto r_c<x<R es

${\displaystyle \frac{d^{2}x}{d{t}^2}=-\frac{4\pi }{3}G\left\{({\rho }_c-{\rho }_m)\frac{r_c^3}{x^2}+{\rho }_{m}x\right\}}$

La partícula parte de la posición x=R con velocidad nula dx/dt=0, llega al núcleo, x=r_c en el instante t_m con velocidad -v_m.

Para calcular v_m aplicamos el principio de conservación de la energía.  La partícula parte del reposo desde un punto que dista r=R del centro de la Tierra y llega a un punto en el interior del manto que dista r centro de la Tierra, con velocidad v.

${\displaystyle -\frac{4\pi }{3}mG\rho R^2=\frac{1}{2}m{v}^2-\frac{4\pi }{3}mG\left(-\frac{1}{2}{\rho }_m{r}^2+({\rho }_c-{\rho }_m)\frac{r_c^3}{r}+\frac{3}{2}{\rho }_m{R}^2\right)}$

Cuando llega al límite entre el manto y el núcleo, r=r_c,

${\displaystyle v_m^2=\frac{8\pi }{3}G\left\{({\rho }_c-\frac{3}{2}{\rho }_m)r_c^2-({\rho }_c-{\rho }_m)\frac{r_c^3}{R}+{\rho }_m\frac{R^2}{2}\right\}}$

Para llegar a esta expresión se ha empleado la relación entre las densidades del manto ρ_m, núcleo ρ_c y media ρ.

${\displaystyle \rho R^3={\rho }_c{r}_c^3+{\rho }_m(R^3-r_c^3)}$

Con los datos:

• G=6.67·10^-11 N·m²/kg²,

• Densidad del núcleo ρ_c=11000 kg/m³,

• Densidad del manto ρ_m=4437 kg/m³,

• Radio del núcleo r_c=3.490·10⁶ m

• Radio de la Tierra R=6.371·10⁶ m

Se obtiene v_m=7444.3 m/s

Para obtener el tiempo t_m, resolvemos la ecuación diferencial de segundo orden por procedimientos numéricos

• La ecuación del movimiento en el núcleo es

${\displaystyle \frac{d^{2}x}{d{t}^2}=-\frac{4}{3}\pi G{\rho }_{c}x}$

En el núcleo la partícula describe un MAS cuya frecuencia angular es

${\displaystyle {\omega }^2=\frac{4}{3}\pi G{\rho }_c}$

ω=0.00175 rad/s

La solución de la ecuación diferencial del movimiento es

x=Asin(ωt)+Bcos(ωt)
dx/dt=Aωcos(ωt)-Bωsin(ωt)

Los coeficientes A y B se determina a partir de las condiciones iniciales

Ponemos el contador de tiempo parcial a cero t=0, x=r_c, dx/dt=-v_m.

${\displaystyle x=\frac{-v_m}{\omega }\text{sin}(\omega t)+r_c\cos (\omega t)}$

El tiempo que tarda en llegar al centro de la Tierra x=0 es

${\displaystyle t_c=\frac{1}{\omega }\arctan \frac{r_c\omega }{v_m}}$

La velocidad que alcanza al llegar al centro de la Tierra es

${\displaystyle \begin{array}{l}v=-v_m\cos (\omega t_c)-r_c\omega \text{sin}(\omega t_c)=-\sqrt{v_m^2+{(r_c\omega )}^2}\\ v^2=\frac{8\pi }{3}G\left(\frac{3}{2}({\rho }_c-{\rho }_m)r_c^2+\frac{1}{2}{\rho }_m{R}^2-({\rho }_c-{\rho }_m)\frac{r_c^3}{R}\right)\end{array}}$     

Para obtener la expresión anterior se ha empleado las relaciones trigonométricas

${\displaystyle \text{sin}\theta =\frac{\tan \theta }{\sqrt{1+{\tan }^2\theta }}\cos \theta =\frac{1}{\sqrt{1+{\tan }^2\theta }}}$

Obtenemos esta velocidad aplicando el principio de conservación de la energía. La partícula parte del reposo desde un punto que dista r=R del centro de la Tierra y llega a un punto en el interior del núcleo que dista r centro de la Tierra, con velocidad v.

${\displaystyle -\frac{4\pi }{3}mG\rho R^2=\frac{1}{2}m{v}^2-\frac{4\pi }{3}mG\left(-\frac{1}{2}{\rho }_c{r}^2+\frac{3}{2}({\rho }_c-{\rho }_m)r_c^2+\frac{3}{2}{\rho }_m{R}^2\right)}$

La velocidad cuando pasa por el centro de la Tierra, para r=0 es

${\displaystyle v^2=\frac{8\pi }{3}G\left(\frac{3}{2}({\rho }_c-{\rho }_m)r_c^2+\frac{1}{2}{\rho }_m{R}^2-({\rho }_c-{\rho }_m)\frac{r_c^3}{R}\right)}$

Con los datos del modelo de Tierra: t_c=392.4 s, v=9635.9 m/s, que en el modelo de Tierra homogénea da el valor v=7913.0 m/s

Completamos un cuarto de oscilación en el tiempo t_m+t_c. El periodo de unas oscilación es P=4(t_c+t_m)

Aproximación para el movimiento de la partícula en el manto

Se describe de forma aproximada el movimiento de la partícula, suponiendo que la aceleración de la gravedad en el manto es constante e igual al valor medio g_m.

En el manto el movimiento de la partícula es rectilíneo uniformemente acelerado

${\displaystyle v=-g_{m}tx=R-\frac{1}{2}{g}_m{t}^2}$

Llega al núcleo x=r_c en el instante t_m con velocidad -v_m.

${\displaystyle t_m=\sqrt{\frac{2(R-r_c)}{g_m}}{v}_m^2=2g_m(R-r_c)}$

Sabiendo que la aceleración de la gravedad media en el manto es g_m=9.57 m/s², la velocidad de la partícula cuando alcanza el núcleo es v_m=7424.5 m/s, un valor próximo al calculado mediante el principio de conservación de la energía o resolviendo la ecuación diferencial del movimiento por procedimientos numéricos.

• El tiempo que tarda en llegar al núcleo desde que se suelta en el extremo del túnel es t_m=775.9 s

• El tiempo que tarda en llegar al origen a través del núcleo es t_c=393.2 s

A partir de aquí, la partícula describe un MAS, hasta que llega a la frontera entre el núcleo y el manto x=-r_c, con velocidad –v_m. En el manto se frena la partícula hasta que llega a x=-R en reposo, completando medio ciclo de oscilación. El periodo P=4(t_c+t_m)=4676.4 s=77.9 min

Túnel a lo largo de una cuerda

Se construye un túnel a lo largo de una cuerda que dista d del centro de la Tierra

La fuerza que mueve la partícula a lo largo del túnel es la componente (en color azul) mg·cosθ=mg·x/r. Siendo r la distancia entre el centro de la Tierra y la partícula r²=x²+d².

• La ecuación del movimiento de la partícula en el manto es

${\displaystyle \begin{array}{l}\frac{d^{2}x}{d{t}^2}=-\frac{4\pi }{3}G\left\{({\rho }_c-{\rho }_m)\frac{r_c^3}{r^2}+{\rho }_{m}r\right\}\frac{x}{r}\\ \frac{d^{2}x}{d{t}^2}=-\frac{4\pi }{3}G\left\{({\rho }_c-{\rho }_m)\frac{r_c^3}{{(x^2+d^2)}^{3/2}}+{\rho }_{m}x\right\}\end{array}}$

La partícula parte de la posición $x=\sqrt{R^2-d^2}$ con velocidad nula dx/dt=0, llega al núcleo, $x=\sqrt{r_c^2-d^2}$ en el instante t_m con velocidad -v_m, que se calculan resolviendo la ecuación diferencial de segundo orden por procedimientos numéricos.

Aplicando el principio de conservación de la energía obtenemos el mismo valor para v_m que el calculado para d=0, es decir, v_m=7444.3 m/s, ya que la partícula parte del reposo desde una posición que dista R del centro de la Tierra y llega a una posición que dista r_c del centro de la Tierra.

• La ecuación del movimiento en el núcleo es

${\displaystyle \begin{array}{l}\frac{d^{2}x}{d{t}^2}=\left(-\frac{4}{3}\pi G{\rho }_{c}r\right)\frac{x}{r}\\ \frac{d^{2}x}{d{t}^2}=-\frac{4}{3}\pi G{\rho }_{c}x\end{array}}$

En el núcleo, la partícula describe un MAS cuya frecuencia angular es la misma que la calculada en el apartado precedente

${\displaystyle {\omega }^2=\frac{4}{3}\pi G{\rho }_c}$

La solución de la ecuación diferencial del movimiento es

x=Asin(ωt)+Bcos(ωt)
dx/dt=Aωcos(ωt)-Bωsin(ωt)

Los coeficientes A y B se determina a partir de las condiciones iniciales

Ponemos el contador de tiempo parcial a cero t=0, $x=\sqrt{r_c^2-d^2}$ , dx/dt=-v_m.

${\displaystyle x=\frac{-v_m}{\omega }\text{sin}(\omega t)+\sqrt{r_c^2-d^2}\cos (\omega t)}$

El tiempo que tarda en llegar al origen x=0 es

${\displaystyle t_c=\frac{1}{\omega }\arctan \frac{\omega \sqrt{r_c^2-d^2}}{v_m}}$

El tiempo que pasa la partícula en el interior del núcleo 4·t_c, va disminuyendo a medida que aumenta d.

Aplicando el principio de conservación de la energía calculamos la velocidad con la que llega la partícula al origen. La partícula parte del reposo desde una posición que dista R del centro de la Tierra y llega a una posición que dista d del centro de la Tierra.

${\displaystyle \begin{array}{l}-\frac{4\pi }{3}mG\rho R^2=\frac{1}{2}m{v}^2-\frac{4\pi }{3}mG\left(-\frac{1}{2}{\rho }_c{d}^2+\frac{3}{2}({\rho }_c-{\rho }_m)r_c^2+\frac{3}{2}{\rho }_m{R}^2\right)\\ v^2=\frac{8\pi }{3}G\left(-\frac{1}{2}{\rho }_c{d}^2+\frac{3}{2}({\rho }_c-{\rho }_m)r_c^2+\frac{1}{2}{\rho }_m{R}^2-({\rho }_c-{\rho }_m)\frac{r_c^3}{R}\right)\end{array}}$

Por ejemplo si d=0.2R=1.274·10⁶ m, v=9373.4 m/s

A partir de aquí, la partícula describe un MAS, hasta que llega a la frontera entre el núcleo y el manto $x=-\sqrt{r_c^2-d^2}$ , con velocidad –v_m. En el manto, se frena la partícula hasta que llega a $x=-\sqrt{R^2-d^2}$ en reposo, completando medio ciclo de oscilación.  El periodo es P=4(t_c+t_m)

Aproximación

Si tomamos la aceleración de la gravedad en el manto constante g_m, la componente de la fuerza sobre la partícula a lo largo del túnel -g_m·x/r deja de ser constante. Por lo que la ecuación del movimiento se ha de integrar aplicando procedimientos numéricos.

${\displaystyle \begin{array}{l}\frac{d^{2}x}{d{t}^2}=-g_m\frac{x}{r}\\ \frac{d^{2}x}{d{t}^2}=-g_m\frac{x}{\sqrt{d^2+x^2}}\end{array}}$

Caso particular d>r_c

La partícula no toca el núcleo, su movimiento transcurre en el manto. La ecuación del movimiento de la partícula en el manto es

${\displaystyle \begin{array}{l}\frac{d^{2}x}{d{t}^2}=-\frac{4\pi }{3}G\left\{({\rho }_c-{\rho }_m)\frac{r_c^3}{r^2}+{\rho }_{m}r\right\}\frac{x}{r}\\ \frac{d^{2}x}{d{t}^2}=-\frac{4\pi }{3}G\left\{({\rho }_c-{\rho }_m)\frac{r_c^3}{{(x^2+d^2)}^{3/2}}+{\rho }_{m}x\right\}\end{array}}$

La partícula parte de la posición $x=\sqrt{R^2-d^2}$ con velocidad nula dx/dt=0, y llega al origen  x=0, en el instante t_m con velocidad -v_m, que se calculan resolviendo la ecuación diferencial de segundo orden por procedimientos numéricos. El periodo es P=4·t_m

Aplicando el principio de conservación de la energía calculamos la velocidad con la que llega la partícula al origen. La partícula parte del reposo desde una posición que dista R del centro de la Tierra y llega a una posición que dista d del centro de la Tierra en el interior del manto.

${\displaystyle \begin{array}{l}-\frac{4\pi }{3}mG\rho R^2=\frac{1}{2}m{v}^2-\frac{4\pi }{3}mG\left(-\frac{1}{2}{\rho }_m{d}^2+({\rho }_c-{\rho }_m)\frac{r_c^3}{d}+\frac{3}{2}{\rho }_m{R}^2\right)\\ v^2=\frac{8\pi }{3}G\left(-\frac{1}{2}{\rho }_m{d}^2+({\rho }_c-{\rho }_m)\frac{r_c^3}{d}+\frac{1}{2}{\rho }_m{R}^2-({\rho }_c-{\rho }_m)\frac{r_c^3}{R}\right)\end{array}}$

Para llegar a esta expresión se ha empleado la relación entre las densidades del manto ρ_m, núcleo ρ_c y media ρ.

${\displaystyle \rho R^3={\rho }_c{r}_c^3+{\rho }_m(R^3-r_c^3)}$

Por ejemplo si d=0.6R=3.823·10⁶ m, v=6965.3 m/s

Actividades

Se introduce

La distancia d, del túnel al centro de la Tierra, en el control titulado Distancia al centro

Se pulsa el botón titulado Nuevo

Se observa el movimiento de la partícula en el túnel que atraviesa la Tierra, cuando parte del reposo en uno de los extremos del túnel. Una flecha indica la componente de fuerza de atracción sobre la partícula a lo largo del túnel.

Se representa la función energía potencial E_p(r) en el núcleo y en el manto.

• La energía total una recta horizontal de color negro, la energía total es negativa y es el segmento AB
• La energía potencial un segmento vertical AC de color rojo,
• La energía cinética un segmento vertical BC de color azul

La pendiente de la curva cambiada de signo, nos proporciona el valor y sentido de la fuerza de atracción que ejerce la Tierra sobre la partícula, una flecha horizontal de color azul.

Se proporciona los datos del

• Tiempo t, en minutos
• La posición x/R a lo largo del túnel, una fracción del radio R de la Tierra
• La velocidad v de la partícula en m/s

Los datos del modelo del interior de la Tierra son:

• G=6.67·10^-11 N·m²/kg²,
• Densidad del núcleo ρ_c=11000 kg/m³,
• Densidad del manto ρ_m=4437 kg/m³,
• Radio del núcleo r_c=3.490·10⁶ m
• Radio de la Tierra R=6.371·10⁶ m

Referencias

Snyder R. Two-density model of the Earth. Am. J. Phys. 54 (6) June 1986, pp. 511-513