Ecuación de Emden-Chandrasekhar - Emden–Chandrasekhar equation

En astrofísica , la ecuación de Emden-Chandrasekhar es una forma adimensional de la ecuación de Poisson para la distribución de densidad de una esfera de gas isotérmica esféricamente simétrica sujeta a su propia fuerza gravitacional, llamada así por Robert Emden y Subrahmanyan Chandrasekhar . La ecuación fue introducida por primera vez por Robert Emden en 1907. La ecuación dice

${\ Displaystyle {\ frac {1} {\ xi ^ {2}}} {\ frac {d} {d \ xi}} \ left (\ xi ^ {2} {\ frac {d \ psi} {d \ xi}} \ right) = e ^ {- \ psi}}$

donde es el radio adimensional y está relacionado con la densidad de la esfera de gas como , donde es la densidad del gas en el centro. La ecuación no tiene una solución explícita conocida. Si se utiliza un fluido politrópico en lugar de un fluido isotérmico, se obtiene la ecuación de Lane-Emden . La suposición isotérmica generalmente se modela para describir el núcleo de una estrella. La ecuación se resuelve con las condiciones iniciales, ${\ Displaystyle \ xi}$${\ Displaystyle \ psi}$${\ Displaystyle \ rho = \ rho _ {c} e ^ {- \ psi}}$${\ Displaystyle \ rho _ {c}}$

${\ Displaystyle \ psi = 0, \ quad {\ frac {d \ psi} {d \ xi}} = 0 \ quad {\ text {at}} \ quad \ xi = 0.}$

La ecuación también aparece en otras ramas de la física, por ejemplo, la misma ecuación aparece en la teoría de la explosión de Frank-Kamenetskii para una vasija esférica. La versión relativista de este modelo isotérmico de simetría esférica fue estudiada por Subrahmanyan Chandrasekhar en 1972.

Derivación

Para una estrella gaseosa isotérmica , la presión se debe a la presión cinética y la presión de radiación.${\ Displaystyle p}$

${\ Displaystyle p = \ rho {\ frac {k_ {B}} {WH}} T + {\ frac {4 \ sigma} {3c}} T ^ {4}}$

dónde

• ${\ Displaystyle \ rho}$ es la densidad
• ${\ Displaystyle k_ {B}}$es la constante de Boltzmann
• ${\ Displaystyle W}$es el peso molecular medio
• ${\ Displaystyle H}$ es la masa del protón
• ${\ Displaystyle T}$ es la temperatura de la estrella
• ${\ Displaystyle \ sigma}$es la constante de Stefan-Boltzmann
• ${\ Displaystyle c}$ es la velocidad de la luz

La ecuación para el equilibrio de la estrella requiere un equilibrio entre la fuerza de presión y la fuerza gravitacional.

${\ Displaystyle {\ frac {1} {r ^ {2}}} {\ frac {d} {dr}} \ left ({\ frac {r ^ {2}} {\ rho}} {\ frac {dp } {dr}} \ right) = - 4 \ pi G \ rho}$

donde es el radio medido desde el centro y es la constante gravitacional . La ecuación se reescribe como ${\ Displaystyle r}$${\ Displaystyle G}$

${\ Displaystyle {\ frac {k_ {B} T} {WH}} {\ frac {1} {r ^ {2}}} {\ frac {d} {dr}} \ left (r ^ {2} { \ frac {d \ ln \ rho} {dr}} \ right) = - 4 \ pi G \ rho}$

Introduciendo la transformación

${\ Displaystyle \ psi = \ ln {\ frac {\ rho _ {c}} {\ rho}}, \ quad \ xi = r \ left ({\ frac {4 \ pi G \ rho _ {c} WH} {k_ {B} T}} \ right) ^ {1/2}}$

donde está la densidad central de la estrella, conduce a ${\ Displaystyle \ rho _ {c}}$

${\ Displaystyle {\ frac {1} {\ xi ^ {2}}} {\ frac {d} {d \ xi}} \ left (\ xi ^ {2} {\ frac {d \ psi} {d \ xi}} \ right) = e ^ {- \ psi}}$

Las condiciones de contorno son

${\ Displaystyle \ psi = 0, \ quad {\ frac {d \ psi} {d \ xi}} = 0 \ quad {\ text {at}} \ quad \ xi = 0}$

Porque , la solución es como ${\ Displaystyle \ xi \ ll 1}$

${\ Displaystyle \ psi = {\ frac {\ xi ^ {2}} {6}} - {\ frac {\ xi ^ {4}} {120}} + {\ frac {\ xi ^ {6}} { 1890}} + \ cdots}$

Limitaciones del modelo

Asumir que la esfera isotérmica tiene algunas desventajas. Aunque la densidad obtenida como solución de esta esfera de gas isotérmica disminuye desde el centro, disminuye demasiado lentamente para dar una superficie bien definida y una masa finita para la esfera. Se puede demostrar que, como , ${\ Displaystyle \ xi \ gg 1}$

${\ Displaystyle {\ frac {\ rho} {\ rho _ {c}}} = e ^ {- \ psi} = {\ frac {2} {\ xi ^ {2}}} \ left [1 + {\ frac {A} {\ xi ^ {1/2}}} \ cos \ left ({\ frac {\ sqrt {7}} {2}} \ ln \ xi + \ delta \ right) + O (\ xi ^ {-1}) \ derecha]}$

donde y son constantes que se obtendrán con solución numérica. Este comportamiento de la densidad da lugar a un aumento de masa con un aumento de radio. Así, el modelo suele ser válido para describir el núcleo de la estrella, donde la temperatura es aproximadamente constante. ${\ Displaystyle A}$${\ Displaystyle \ delta}$

Solución singular

La introducción de la transformación transforma la ecuación a ${\ Displaystyle x = 1 / \ xi}$

${\ Displaystyle x ^ {4} {\ frac {d ^ {2} \ psi} {dx ^ {2}}} = e ^ {- \ psi}}$

La ecuación tiene una solución singular dada por

${\ Displaystyle e ^ {- \ psi _ {s}} = 2x ^ {2}, \ quad {\ text {o}} \ quad - \ psi _ {s} = 2 \ ln x + \ ln 2}$

Por lo tanto, se puede introducir una nueva variable como , donde se puede derivar la ecuación para , ${\ Displaystyle - \ psi = 2 \ ln x + z}$${\ Displaystyle z}$

${\ Displaystyle {\ frac {d ^ {2} z} {dt ^ {2}}} - {\ frac {dz} {dt}} + e ^ {z} -2 = 0, \ quad {\ text { donde}} \ quad t = \ ln x}$

Esta ecuación se puede reducir a primer orden introduciendo

${\ Displaystyle y = {\ frac {dz} {dt}} = \ xi {\ frac {d \ psi} {d \ xi}} - 2}$

entonces tenemos

${\ Displaystyle y {\ frac {dy} {dz}} - y + e ^ {z} -2 = 0}$

Reducción

Hay otra reducción debida a Edward Arthur Milne . Definamos

${\ Displaystyle u = {\ frac {\ xi e ^ {- \ psi}} {d \ psi / d \ xi}}, \ quad v = \ xi {\ frac {d \ psi} {d \ xi}} }$

luego

${\ Displaystyle {\ frac {u} {v}} {\ frac {dv} {du}} = - {\ frac {u-1} {u + v-3}}}$

Propiedades

• Si es una solución de la ecuación de Emden-Chandrasekhar, entonces también es una solución de la ecuación, donde es una constante arbitraria.${\ Displaystyle \ psi (\ xi)}$${\ Displaystyle \ psi (A \ xi) -2 \ ln A}$${\ Displaystyle A}$
• Las soluciones de la ecuación de Emden-Chandrasekhar que son finitas en el origen tienen necesariamente en${\ Displaystyle d \ psi / d \ xi = 0}$${\ Displaystyle \ xi = 0}$

Ver también

• Ecuación de Lane-Emden
• Teoría de Frank-Kamenetskii
• Ecuación de la enana blanca de Chandrasekhar

Referencias