En Física se diseñan y utilizan modelos geométricos para simplificar la complejidad de los elementos físicos reales. Se usan densidades homogéneas, radios y superficies perfectas, para como elementos de 1 ,2 o 3 dimensiones, todo esto es para poder obtener resultados matemáticos que se pueden analizar y extrapolar a otros modelos más complejos. Esto se utiliza en toda la física, en Mecánica de rotación, en electromagnetismo, en campos gravitatorios, etc. Esto es el día a día de la física y se utiliza desde la época de Newton quien descubrió el cálculo diferencial y el cálculo integral.
Como ejemplo esto se puede observar y estudiar en el siguiente enlace, donde vienen modelos de masas lineales, masas en forma de circunferencias lineales, masas en superficies, para calcular el campo de elementos de complejidad superior.
http://hyperphysics.phy-astr.gsu.edu/hbasees/electric/elelin.html
Por lo tanto, esto no es un problema para que salgan cosas raras en las esferas huecas de espesor infinitesimal.
Totalmente de acuerdo hasta ahí.
Ahora con un ejemplo voy a tratar de explicar la potencia de esto que he dicho anteriormente, para el cálculo de la gravedad que crea un casquete esférico comprendidos entre los radios Rint y Rext.
a) Vamos a suponer que r es un punto exterior a la esfera r>Rext, para este caso:
\( dg=-\frac{G dm}{r^2} \) es la gravedad que genera una esfera hueca de espesor infinitesimal en el punto r que dista del centro de la esfera.
\( g\left ( r \right )=\int_{Rint}^{Rext}dg=\int_{Rint}^{Rext}-\frac{G dm}{r^2}=-\frac{G}{r^2}\int_{Rint}^{Rext} dm=-\frac{G M}{r^2},, M=M_{total esfera} \)
El cálculo es enormemente rápido y sencillo, no hay que hacer integrales triples y el resultado es el mismo. Aquí queda demostrado la potencia de utilizar este tipo de elementos matemáticos ideales.
Ciertamente es rápido y sencillo. Además, yo añadiría que el cálculo no tiene ningún sentido, ni físico ni matemático, pero, como el resultado final es correcto, no merece la pena discutir tu cálculo, al menos de momento.
b) Vamos a suponer que r es un punto donde r= Rext y esto sirva como ejemplo para cualquier punto \( r\geq{Rint} \) y \( r\leq{Rext} \):
Aquí tenemos que:\( dg=-\frac{G dm}{2R_ext^2} \) para la esfera elemental r=Rext
\( dg=-\frac{G dm}{r^2} \) para el resto de las esferas elementales \( Rint\leq{r}< Rext \)
Como vemos en la primera esfera, en r=Rext, dg tiene el valor la mitad que el resto de esferas elementales que están comprendidas en Rint y Rext. Por lo tanto:
\( \color{red}g\left ( r \right )=\int_{Rint}^{Rext}dg=-\frac{GM}{R_{ext}^2}+ \frac{G dm_{esferaR_{ext}}}{2R{ext}^2}, \)
donde \( dm_{esferaR_{ext}}=\rho4\pi R_{ext}^2 dR \)
Sustituyendo:\( g\left ( r \right )=-\frac{GM}{R_{ext}^2}+ \frac{G \rho4\pi R_{ext}^2 dR}{2R{ext}^2}=-\frac{GM}{R_{ext}^2}+ 2\pi G \rho dR \)
\( \color{blue}\rho dR=\sigma_{\rm esfera Rext} \)
sustituyendo:
\( g\left ( r \right )=-\frac{GM}{R_{ext}^2}+ 2\pi G \sigma_{esferaRext} \)
Como saben en matemáticas 2 cosas son iguales si son totalmente iguales y por ello vemos que el resultado no es el oficial pues hay un término adicional con lo cual me ratifico que los resultados de la integración de esferas de espesores finitos oficiales tienen un error también.
Es difícil discutir esto, porque tus cuentas no tienen ni pies ni cabeza, pero, afortunadamente, el error esencial en el que incurres no depende de la inexactitud de todo lo demás, así que te lo puedo señalar igualmente.
Imagina que quieres calcular una integral sencilla, como la de la función \( g_1(r)= r^2 \) en el intervalo \( [1, 5] \). Es fácil:
\( \displaystyle \int_1^5g_1(r)\,dr = \int_1^5 r^2\,dr = \left[\frac{r^3}3\right]_1^5 = \frac{125}3-\frac13 = 41.33\ldots \)
Pero entonces vienes y me dices,
¡eh!, que no es ésa la función que hay que integrar. La correcta es:\( \displaystyle g_2(r) = \cases{r^2&si $r\neq 4$,\cr 8&si $r=4$,} \)
de modo que para \( r=4 \) vale justo la mitad de lo que has considerado tú. Por lo tanto la integral está mal.Y yo te respondo ¡no! ¡Da igual! También se cumple que
\( \displaystyle \int_1^5g_2(r)\,dr = \int_1^5 r^2\,dr = \left[\frac{r^3}3\right]_1^5 = \frac{125}3-\frac13 = 41.33\ldots \)
porque un principio fundamental del cálculo integral dice que:
Si modificamos un integrando en un número finito de puntos, la integral no se ve alterada.Por lo tanto, da igual, a la hora de calcular la integral, que la función \( g(r) \) valga \( 8 \) o \( 16 \) para \( r=4 \), porque el resultado será el mismo.
Lo mismo sucede con tu ejemplo. Que el valor de \( dg \) sea distinto en \( r= R_{\rm ext} \) no altera el valor de una integral en el intervalo \( [R_{\rm int}, R_{\rm ext}] \).
En particular, la igualdad que te he puesto en rojo está mal (entre otras muchas cosas), porque no puedes traducir que un integrando tome un valor distinto en un punto en que la variación en ese punto añada un término nuevo al resultado. Variar el valor de un integrando en un punto no suma nada a una integral, no la modifica.
Aunque todos tus cálculos están llenos de incoherencias, la más relevante en lo que nos ocupa es lo que te he puesto en azul: a la izquierda tienes una expresión "infinitesimal", con un \( dR \), mientras que a la derecha tienes un valor finito. Necesitas hacer una pirueta así para deshacerte del sumando "infinitesimal" que te ha quedado después de haber calculado ya la integral, pero es es una trampa muy gorda, que no se nota mucho porque es un sinsentido en medio de un mar de sinsentidos.
No te puedo concretar más tu error porque eso requiere entrar en los cálculos que haces que, como digo, no tienen sentido alguno. Lo máximo que puedo hacer es hacer bien tus cálculos, es decir, calcular el campo determinado por un cascarón esférico a partir del campo de una esfera sin grosor que das por conocido. Richard R Richard ha hecho el cálculo en un mensaje anterior partiendo de cero, pero si prefieres hacerlo a partir del campo de una esfera sin grosor, también se puede hacer como a ti te gusta, y sale lo mismo, pero para ello hay que ser cuidadoso al combinar la distribución tridimensional de masa del cascarón con las distribuciones bidimensionales de las esferas en que pretendes descomponerlo (cosa que tú no haces). Vamos a ello:
El resultado bidimensional del que quieres partir, enunciado correctamente, es éste:
\( \displaystyle \int_{S_R}\frac{Gdm}{\|\bar x-\bar p\|^3}(\bar x-\bar p)=-\frac{GM}{r^2}\bar u \),
donde \( S_R \) es la esfera de radio \( R \) con centro en el origen de coordenadas, \( \bar p \) es un punto situado a una distancia \( r>R \) del centro de la esfera y \( \bar u=\bar p/\|\bar p\| \) es el vector radial unitario. Si \( r<R \) la integral da \( 0 \) y si \( r=R \) la integral da lo mismo salvo por que hay que poner un \( 2 \) adicional en el denominador, aunque eso no hace falta para nada.
Lo primero que necesitamos aclarar, para que se entienda qué es esa integral, es que \( dm = \rho_R\,d\bar x \), donde \( \rho_R \) es la densidad superficial de la esfera de radio \( R \), que suponemos constante, y \( d\bar x \) es el elemento de área. Similarmente, la masa \( M \) la podemos expresar como \( M = \rho_R V_R \), donde \( V_R \) es la superficie de la esfera de radio \( R \). Explícitamente:
\( \displaystyle \int_{S_R}\frac{G\rho_Rd\bar x}{\|\bar x-\bar p\|^3}(\bar x-\bar p)=-\frac{G\rho_R}{r^2}4\pi R^2\bar u \),
y, simplificando:
\( \displaystyle \int_{S_R}\frac{d\bar x}{\|\bar x-\bar p\|^3}(\bar x-\bar p)=-\frac{1}{r^2}4\pi R^2\bar u \).
Con esta integral podemos calcular el campo generado por el cascarón \( C \) de centro el origen de coordenadas y radios \( R_1<R_2 \):
\( \displaystyle \int_C\frac{G\,dm}{\|\bar x-\bar p\|^3}(\bar x-\bar p) \),
donde ahora \( dm \) es un elemento de masa tridimensional, que podemos expresar como \( dm = \rho\,d\bar x \), donde \( \rho \) es la densidad (tridimensional) del cascarón, que suponemos constante, y ahora \( d\bar x \) es el elemento de volumen. Así nos queda
\( \displaystyle \int_C\frac{G\rho\,d\bar x}{\|\bar x-\bar p\|^3}(\bar x-\bar p) \).
Ahora usamos una versión del teorema de Fubini para variedades diferenciales por la que la integral sobre un cascarón se puede descomponer así:
\( \displaystyle \int_C\frac{G\rho\,d\bar x}{\|\bar x-\bar p\|^3}(\bar x-\bar p)= G\rho\int_{R_1}^{R_2}\int_{S_R}\frac{d\bar x}{\|\bar x-\bar p\|^3}(\bar x-\bar p)\,dR \).
Aquí \( R \) varía en el intervalo \( [R_1, R_2] \) y, para cada \( R \), la \( \bar x \) varía ahora en la esfera de radio \( R \), de modo que el integrando es una de las integrales de las que hemos hablado antes. Si el punto \( \bar p \) está a una distancia \( r \) del centro del cascarón y \( R_1\leq r\leq R \), tenemos que el integrando vale \( 0 \) siempre que \( r<R \), por lo que queda:
\( \displaystyle \int_C\frac{G\rho\,d\bar x}{\|\bar x-\bar p\|^3}(\bar x-\bar p)= G\rho\int_{R_1}^r\int_{S_R}\frac{d\bar x}{\|\bar x-\bar p\|^3}(\bar x-\bar p)\,dR \).
Y si usamos el valor de la integral interior que hemos obtenido antes, queda:
\( \displaystyle \int_C\frac{G\rho\,d\bar x}{\|\bar x-\bar p\|^3}(\bar x-\bar p)= -G\rho\int_{R_1}^r\frac{1}{r^2}4\pi R^2\bar u\,dR = -\frac{G\rho 4\pi}{r^2}\int_{R_1}^r R^2\,dR\, \bar u \).
Y éste es el punto donde tú tendrías que decir ¡No! ¡Está mal!, porque cuando has sustituido la integral interior por su valor, lo que haces es correcto si \( R<r \), pero para \( R=r \) te falta un \( 2 \) en el denominador. Y aquí es donde te digo que eso da igual: si queremos calcular la integral en el intervalo \( [R_1,r] \) de la función
\( g(R)=\frac{1}{r^2}4\pi R^2\bar u \),
da igual el valor que asignemos a \( g(r) \), porque el valor de una integral no se altera si modificamos el integrando en un único punto \( R= r \). En particular, no falta ningún "término adicional".
Por simplificar, voy a considerar, concretamente, que el punto \( \bar p \) está en el borde exterior del cascarón, es decir, que \( r= R_2 \), que es el caso que tú considerabas. Entonces, el módulo del campo gravitatorio queda
\( \displaystyle \frac{G\rho 4\pi}{r^2}\int_{R_1}^{R_2} R^2\,dR= \frac{G\rho 4\pi}{r^2}\left[\frac{R^3}3\right]_{R_1}^{R_2}=\frac{G\rho}{r^2}\frac43\pi(R_2^3-R_1^3)= \frac{GM}{r^2} \), que es el resultado que puedes ver en cualquier parte, sin ninguna presunta alteración, en contra de lo que pretendes.
Si haces el cálculo para \( R_1<r<R_2 \) sale la expresión dada por AlbertR en un mensaje anterior:
\( \displaystyle \frac{G\rho 4\pi}{r^2}\int_{R_1}^r R^2\,dR = \frac{G\rho 4\pi}{r^2}\frac{r^3-R_1^3}3 = \frac{G\rho}{r^2}\frac{4\pi}{3}(R_2^3-R_1^3)\frac{r^3-R_1^3}{R_2^3-R_1^3}=\frac{GM}{R_2^3-R_1^3} \frac{r^3-R_1^3}{r^2} \)
y si hacer tender \( R_2 \) y \( R_1 \) a un valor \( r \) común, el resultado que obtienes es distinto según cómo lo hagas. Por ejemplo, si haces \( R_1 = r-a\epsilon \) y \( R_2 = r+b\epsilon \), la intensidad del campo gravitatorio del cascarón en un punto que dista \( r \) del centro, cuando \( \epsilon \) tiende a \( 0 \), se aproxima a
\( \displaystyle\frac{GM}{\color{red}r^2}\frac a{{\color{red} a+}b} \),
Corregidas DOS erratas.que puede ser cualquier cosa, según como elijas \( a, b \). Por eso no tiene sentido físico atribuir una intensidad al campo gravitatorio de una esfera sin grosor en un punto de su superficie, porque en cualquier caso real la esfera tendrá un grosor, por pequeño que sea, y el resultado será muy distinto si el punto de la "superficie" que tomas está en su cara interna (ahí será \( 0 \)) en su cara externa (donde será \( GM/r^2 \)) o bien en un punto intermedio, donde podrá tomar cualquier valor intermedio.
¿Qué pasa si hacemos que Rint lo hacemos tender a Rext, es decir, que se hace el espesor menor (pero con la misma masa M de la esfera)? Lo que ocurre es que cada circunferencia elemental va ganando en masa, debido a que la masa se está concentrando y esto provoca que el término adicional vaya ganando en importancia. de tal forma que cuando Rint=Rext lo que pasa es que toda la masa está toda incluida en la esfera de radio Rext y espesor infinitesimal, con una superficie \( S=4\pi R_{ext}^2 \) y una \( \sigma=\frac{M{totalesfera}}{4\pi R_{ext}^2} \)
Con ello en el límite Rext=Rint
\( g\left ( r \right )=-\frac{GM}{R_{ext}^2}+ \frac{GM}{2R_{ext}^2}=- \frac{GM}{2R_{ext}^2} \)
La masa no ha desaparecido, se ha reconcentrado en la superficie de la esfera hueca de espesor infinitesimal y el resultado de la gravedad en su superficie sale lo que ya llevo diciendo durante estos días.
No, nada de eso tiene sentido. No puedes jugar tan alegremente con densidades tridimensionales y densidades superficiales. El "término adicional" directamente no existe y, si bien puedes, efectivamente, considerar cascarones cada vez más finos siempre con la misma masa, el valor del campo gravitatorio no tenderá a un valor fijo, sino que puedes hacer que tienda a cualquier cosa según la proporción con la que se acercan a \( r \) los dos radios.
Por último por favor, coger la página dónde indiqué el cálculo de la gravedad y supongan desde el principio que el punto de cálculo está en R. Verán que la integral existe, está acotada y sale lo que yo os estoy diciendo.
He realizado esta cálculo en coordenadas cartesianas y ocurre lo mismo siempre la integral se calcula normalmente y el valor es el de la gravedad que genera un plano infinito de espesor infinitesimal.
La integral es correcta y da lo que tú dices, pero:
1) Ese valor es irrelevante para el cálculo del campo de un cascarón (por las razones que te he explicado), como lo es para cualquier otro cálculo.
y
2) Ese valor no tiene una interpretación física, porque sólo puedes considerar un modelo bidimensional cuando éste aproxima a los modelos tridimensionales de grosor despreciable, pero en este caso no es así, sino que al sustituir la esfera bidimensional por un cascarón, el valor del campo en lo que era un punto de la esfera dependerá sustancialmente de su posición en el cascarón.