Así que hace un tiempo hice un pequeño proyecto en el que cogí un "modelo solar estándar" de este documento que me da información útil para hacer una estimación. (Como es lógico, el enlace que se da para descargar los datos ha cambiado en los últimos diez años; no he investigado para ver si los datos siguen estando disponibles públicamente).
-
Sólo alrededor del 1,5% de la masa del sol es algo distinto al hidrógeno y al helio-4. Esto es así desde el núcleo hasta la superficie. Supondremos que el sol sólo contiene hidrógeno y helio-4.
-
Toda la masa del sol, excepto el 0,2% más externo (hasta el 90% del radio del sol) está a una temperatura $kT>54\,\mathrm{eV}$ , que es la energía necesaria para girar $\mathrm{He^+}$ en $\mathrm{He^{2+}}$ . (Esta energía es cuatro veces la energía de Rydberg.) Así que en algún lugar por encima del 99% de la masa del sol está completamente ionizada.
-
La temperatura central $kT\approx 1300\,\mathrm{eV}$ es mucho menor que la masa del electrón, por lo que la materia del núcleo no es relativista.
-
Voy a suponer que los electrones no están degenerados; esta herramienta (vía esta pregunta ) me hace pensar que es una suposición bastante segura para la materia en el núcleo con densidad $\rho \approx 150\,\mathrm{g/cm^3}$ y la temperatura $T \approx 10^7\,\mathrm K$ .
En ese caso podemos tratar el núcleo del sol como una mezcla de tres gases ideales que no interactúan, $\mathrm H^+$ , $\mathrm{He}^{2+}$ y $\mathrm e^-$ . Como dice George Herold, cada partícula de gas ideal tiene una energía cinética media $\frac32 kT$ por lo que querremos las densidades numéricas. La densidad numérica del hidrógeno $n_\mathrm{H}$ es $$ n_\mathrm{H} = \rho f_\mathrm{H}/{\mu_\mathrm{H} } $$ donde $\rho$ es la densidad de masa, $f_\mathrm{H}$ es la fracción de masa de hidrógeno, y $\mu_\mathrm{H} = 1\,\mathrm{gram/mole}$ es la masa atómica del hidrógeno. Se tiene una expresión similar para el helio (con $\mu_\mathrm{He} = 4\,\mathrm{gram/mole}$ ). La densidad del número de electrones, gracias a la ionización completa, es sólo $$ n_\mathrm{e} = n_\mathrm{H} + 2n_\mathrm{He}. $$ Aquí hay una figura que muestra la temperatura, la densidad de masa y la composición de mi fuente anterior y la densidad numérica calculada aquí:
Tenga en cuenta que la escala horizontal (radio) está ponderada por la masa: se encuentra aproximadamente la mitad de la masa del sol entre 0,1 y 0,3 radios solares, por lo que ese intervalo ocupa aproximadamente la mitad del eje horizontal. Esto es puramente una técnica de visualización, para que tu ojo no se distraiga con las capas externas (relativamente) frías y difusas del sol.
Para encontrar el total densidad de energía térmica, tenemos que integrar Encontramos la densidad de energía térmica $$ \epsilon = (n_\mathrm{H} + n_\mathrm{He} + n_\mathrm{e})\frac32 kT $$ y el volumen de una cáscara delgada de radio $r$ es $$ dV = 4\pi r^2 dr $$ Esta integral $\int\epsilon\, dV$ me da una energía cinética total almacenada $E=3.09\times10^{41}\,\mathrm{J}$ de los cuales aproximadamente el 95% está contenido en la mitad del radio del sol.
Ahora bien, si el sol tuviera una densidad uniforme se podría estimar su energía potencial gravitatoria la energía que se liberó cuando todas las piezas cayeron juntas, como $$ U_\text{uniform sphere} = -\frac35 \frac{GM_\text{sphere}^2}{R_\text{sphere}} = - 2.3\times10^{41}\,\mathrm J \text{ (uniformly dense sun)}. $$ Eso está muy cerca de nuestro calor almacenado. Podemos hacerlo un poco mejor, ya que conocemos el perfil de densidad del sol, encontrando la energía potencial liberada a medida que se deposita cada cáscara esférica, $$ U = - \int_0^{M_\text{sun}} \frac{G M_\text{enclosed}(r)}{r} dM = -6.15\times10^{41}\,\mathrm{J}. $$ Esta autoenergía gravitacional es aproximadamente el doble de la energía cinética almacenada --- que un verdadero astrónomo habría predicho como consecuencia de la teorema del virial .
0 votos
Se considera que las estrellas tienen capacidades térmicas negativas: adsabs.harvard.edu/full/1977MNRAS.181..405L
0 votos
@DavePhD, supongo que es porque cuando una estrella pierde energía, se calienta más al contraerse más. Pero aun así, debe haber una fuente de energía (energía potencial en este caso) para realizar el calentamiento.
0 votos
@AbanobEbrahim, como has dicho, el calor específico es una función de la temperatura. Debido a esto, tu suposición de que el calor específico a 250 K es válido para una temperatura de 15 millones de K es una extrapolación extrema que definitivamente producirá un gran error en tu respuesta. Y, añado, "no", no sé cómo estimar buenos datos para este problema.