Logo Studenta

guias-de-astro

¡Este material tiene más páginas!

Vista previa del material en texto

Guías de problemas de Astrofísica
6 de diciembre de 2018
Los resueltos de
ALF
N
o
 h
ay
 p
ro
bl
em
a Wi
lly!
Sobre estos apuntes Estos apuntes/resueltos que usted está viendo fueron creados por un alumno mientras cursaba la
materia. Es por ello que podrían haber errores de tipeo, errores conceptuales, de interpretación en los resultados, etc. Use
estos apuntes con precaución. Estos apuntes no son oficiales de ninguna cátedra. Lea atentamente el prospecto. En caso de
notar algún efecto adverso suspenda inmediatamente su uso y consulte con su profesor de cabecera.
El alumno autor de estos apuntes cursó la materia durante el segundo cuatrimestre de 2018, este link conduce a la
página oficial del curso.
Encontrá más resueltos de Alf en este link.
Box 1 - ¿Cómo se hacen estos apuntes?
Estos apuntes están hechos usando un programa llamado Lyxa. Para hacer los dibujos se usó Inkscape y después se
insertó las imágenes en formato svgb directamente en Lyx.
En este repositorio de GitHub se encuentra la plantilla (template) que Alf usa actualmente, con todo lo necesario
para compilarla y empezar a divertirse.
aLyx es una interfaz gráfica para Latex que hace que la escritura se vuelva extremadamente fluida y veloz (al punto de poderse
tomar apuntes en vivo durante una clase).
bsvg es el formato nativo de Inkscape.
1
https://losresueltosdealf.wordpress.com/
http://materias.df.uba.ar/astrofa2018c2/
https://losresueltosdealf.wordpress.com/
https://www.lyx.org/
https://inkscape.org/es/
https://github.com/SengerM/lyx
ÍNDICE ÍNDICE
Índice
Guía 1 - El Universo 4
Problema 1 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 4
Problema 3 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 4
Problema 4 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 4
Problema 5 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 5
Problema 6 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 5
Problema 8 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 5
Problema 9 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 6
Problema 13 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 6
Problema 15 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 6
Problema 17 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 6
Problema 18 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 7
Problema 19 FALTA . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 7
Problema 21 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 7
Guía 2 - Fluidos clásicos 9
Problema 2 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 9
Ítem b . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 9
Ítem c . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 9
Problema 4 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 10
Ítem a . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 10
Problema 5 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 11
Problema 7 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 11
Problema 8 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 12
Ítem a . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 12
Ítem b . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 12
Problema 9 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 13
Ítem a . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 13
Ítem b . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 14
Guía 3 - Relatividad 15
Problema 2 (No entiendo enunciado) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 15
Ítem a . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 15
Problema 3 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 15
Ítem a . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 15
Ítem b . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 15
Problema 4 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 15
Ítem a . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 15
Problema 5 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 16
Ítem a . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 16
Problema 6 (averiguar si esto está bien) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 17
Problema 7 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 18
Guía 4 - Atmósferas estelares 19
Problema 1 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 19
Ítem a (falta una cosa!) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 19
Ítem c . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 19
Problema 3 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 19
Ítem a . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 19
Problema 6 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 20
Ítem a . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 20
Problema 7 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 21
Ítem a . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 21
Problema 8 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 22
Ítem a . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 22
Ítem b . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 22
Ítem c . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .. . . . . . . . . . . . . . 22
Problema 9 FALTA . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 23
Los resueltos de
ALF
N
o
 h
ay
 p
ro
bl
em
a Wi
lly!
←click me!
2 Si ves un error, avisale a Alf , 201812062058 CC0
https://losresueltosdealf.wordpress.com/
https://losresueltosdealf.wordpress.com/contacto/
https://en.wikipedia.org/wiki/Public_domain
ÍNDICE ÍNDICE
Guía 5 - Medio interestelar 24
Problema 2 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 24
Problema 3 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 24
Ítem a . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 24
Ítem b PREGUNTAR . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 25
Guía 6 - Estructura y evolución estelar 27
Problema 1 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 27
Ítem b . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 27
Problema 2 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 27
Problema 5 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 28
Ítem a . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 28
Ítem b . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 28
Ítem c . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 29
Guía 7 - Campos magnéticos 30
Problema 1 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 30
Ítem a . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 30
Problema 3 - MHD . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 30
Ítem b . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 31
Problema 4 - Reconexión de Sweet-Parker . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 32
Ítem a . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 32
Ítem b . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 32
Problema 5 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 34
Ítem a . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 34
Ítem b . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 34
Guía 8 - Galaxias 36
Problema 1 (LSR) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 36
Problema 2 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 36
Problema 3 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 38
Ítem a . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 38
Ítem b . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 39
Guía 9 - Cosmología 40
Problema 2 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 40
Ítem a . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 40
Ítem b . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 40
Ítem c . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 41
Ítem d . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 41
Ítem e . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 42
Ítem f . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 42
Ítem g . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 42
Formulas 44
PREGUNTAS
1. PREGUNTA 1 - Escala hidrodinámica . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 27
2. PREGUNTA 2 - Energía potencial gravitatoria? . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 29
3. PREGUNTA 3 - Qué es el teorema del virial? . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 29
Los resueltos de
ALF
N
o
 h
ay
 p
ro
bl
em
a Wi
lly!
←click me!
3 Si ves un error, avisale a Alf , 201812062058 CC0
https://losresueltosdealf.wordpress.com/
https://losresueltosdealf.wordpress.com/contacto/
https://en.wikipedia.org/wiki/Public_domain
GUÍA 1 - EL UNIVERSO
Guía 1 - El Universo
Problema 1
Los portadores de información que se me ocurren son los siguientes:
Radiación electromagnética. Esto abarca todas las partes del espectro: radio, micro ondas, infrarojo, visible, ultravioleta,
rayos X y rayos γ.
Rayos cósmicos. Distintos tipos de partículas que provienen del espacio, por ejemplo protones, neutrinos, etc.
Radiación gravitatoria. Esto es lo más reciente, la primera onda gravitatoria se midió en 2015 en LIGO.
Problema 3
La radiación de cuerpo negro tiene una densidad espectral dada por la ley de Planck que es
I =

2hν3
c2
1
exp
(
hν
kT
)
− 1
→ En términos de ν
2hc2
λ5
1
exp
(
hc
λkT − 1
) → En términos de λ
Para saber en qué zona del espectro se encuentra la mayor radiación simplemente hay que calcular el máximo de esto. Es
decir, derivar e igualar a cero. El problema es que queda una ecuación trascendente
λ0 =
ch
5kT
1
1− exp
(
− chλ0kT
)
No se me ocurre forma no-numérica de resolver este problema.
Problema 4
La frecuencia máxima es
νmáx =
1,6× 10−19 C× 1× 10−8 T
4π × 9,3× 10−31 kg × 3× 108 m s−1
(
1,6× 10−10 J
9,3× 10−31 kg × (3× 108 m s−1)2
)
= 1,6 Hz
Este valor de frecuencia está en el espectro ELF (extremely low frequency). Para detectar frecuencias tan bajas se utilizan
instalaciones especiales, por ejemplo la siguiente ground dipole antenna
La Nebulosa del Cangrejo de ve así:
Los resueltos de
ALF
N
o
 h
ay
 p
ro
bl
em
a Wi
lly!
←click me!
4 Si ves un error, avisale a Alf , 201812062058 CC0
https://es.wikipedia.org/wiki/Ley_de_Planck
https://en.wikipedia.org/wiki/Extremely_low_frequency
https://losresueltosdealf.wordpress.com/
https://losresueltosdealf.wordpress.com/contacto/
https://en.wikipedia.org/wiki/Public_domain
Problema 5 GUÍA 1 - EL UNIVERSO
La imagen la tomé de este link.
Problema 5
Fácil:
ν = E
h
= 1,45× 109 Hz
y en longitud de onda
λ = c
ν
= 20,6 cm
Esta línea de emisión, conocida como “la línea de 21 cm”, se usa en radiostronomía para estudiar las nubes de polvo, según
Wikipedia.
Problema 6
Ver apuntes de la práctica.
Problema 8
La luminosidad del sol será
L� = 1,38
kW
m2 × Superficie de esfera con radio 1 AU
= 36× 1022 kW
Para calcular la temperatura del sol utilizamos la ley de Stephan-Boltzmann E = σT 4 con lo cual
T� efectiva =
(
L�
4πR2�σ
) 1
4
≈ 5700 K
Los resueltos de
ALF
N
o
 h
ay
 p
ro
bl
em
a Wi
lly!
←click me!
5 Si ves un error, avisale a Alf , 201812062058 CC0
http://hubblesite.org/image/4028/newshttps://es.wikipedia.org/wiki/L%C3%ADnea_de_hidr%C3%B3geno
https://es.wikipedia.org/wiki/Ley_de_Stefan-Boltzmann
https://losresueltosdealf.wordpress.com/
https://losresueltosdealf.wordpress.com/contacto/
https://en.wikipedia.org/wiki/Public_domain
Problema 9 GUÍA 1 - EL UNIVERSO
Problema 9
De acuerdo con la tercera ley de Kepler se tiene que
T 2 = 4π
2
G (m1 +m2)
a3→ Kepler
por lo tanto se puede despejar M = m1 +m2 que es la masa del sistema.
Problema 13
El paralaje anual se explica mediante el siguiente dibujito:
θ
Usando trigonometría se tiene que tanα = 1 AUd de donde se puede despejar
d = 1 AUtanα
α ≈ 0→ ≈ 1 AU
α
Para α = 0,75′′ = 0,0002181662 rad se obtiene
d ≈ 2,77× 105 AU ≈ 4,38 ly
De acuerdo con Wikipedia esto es correcto.
Problema 15
Es una combinación de problemas previos. El dibujito es así
Tierra →
← Sistema binario
1.2''
12 pc
a
Entonces el semieje mayor de la elipse que describe la órbita es a = 12 pc × tan
(
1,2′′
2
)
. Conociendo este semieje y el
período orbital (consigna) se puede calcular la masa del sistema utilizando la tercera ley de Kepler T 2 = 4π
2
G(m1+m2)a
3.
Problema 17
Se puede obtener la velocidad del astro con respecto al laboratorio utilizando el corrimiento Doppler. Si ν es la frecuencia
de la fuente en su propio frame y ν′ es la frecuencia observada en un frame que se mueve con velocidad β en un ángulo θ
según el siguiente dibujito
θ
entonces
ν′ = ν (1− β cos θ) γ→ Doppler
Debido a que λobservada > λemitida entonces νobservada < νemitida lo cual implica que el astro se está alejando.
Los resueltos de
ALF
N
o
 h
ay
 p
ro
bl
em
a Wi
lly!
←click me!
6 Si ves un error, avisale a Alf , 201812062058 CC0
https://es.wikipedia.org/wiki/Leyes_de_Kepler
https://es.wikipedia.org/wiki/Paralaje_anual
https://es.wikipedia.org/wiki/Alfa_Centauri
https://losresueltosdealf.wordpress.com/
https://losresueltosdealf.wordpress.com/contacto/
https://en.wikipedia.org/wiki/Public_domain
Problema 18 GUÍA 1 - EL UNIVERSO
Problema 18
El dibujito del problema es así:
Vega
Cúmulo NGC104
dTierra-NGC104
dTierra-Vega=dSol-Vega
La luminosidad de la estrella Vega es
LVega =
˛
FVega · ds
=
˛
r=dTierra-Vega
dFVegacTierra
= 4πd2Tierra-Vega dFVegacTierra
Por otro lado la luminosidad de la estrella RR-Lyrae es
LRR-Lyrae = 1,05 LVega
4πd2Tierra-NGC104 dFRR-LyraecTierra =
2,29× 10−6 × 4πd2Tierra-NGC104 dFVegacTierra = ← dFRR-LyraecTierra = 2,29× 10
−6 dFVegacTierra
Juntando las dos ecuaciones tenemos que
2,29× 10−6
1,05 × d
2
Tierra-NGC104 = d2Tierra-Vega
por lo tanto
dTierra-NGC104 =
√
1,05
2,29× 10−6 dTierra-Vega
= 5250 pc
Problema 19 FALTA
Qué es el “flujo”?
Problema 21
De acuerdo a la tercera ley de Kepler T 2 = 4π
2a3
G(mVía lactea+m�) donde T es el período orbital y a el semieje mayor de la
elipse. Considerando mVía Lactea � m� se obtiene que
mVía Lactea =
4π2a3
GT 2
Por otro lado, para una órbita circular como propone la consigna, se tiene que a = R0 por lo que
mVía Lactea =
4π2R30
GT 2
Los resueltos de
ALF
N
o
 h
ay
 p
ro
bl
em
a Wi
lly!
←click me!
7 Si ves un error, avisale a Alf , 201812062058 CC0
https://losresueltosdealf.wordpress.com/
https://losresueltosdealf.wordpress.com/contacto/
https://en.wikipedia.org/wiki/Public_domain
Problema 21 GUÍA 1 - EL UNIVERSO
Finalmente, el período orbital T se puede conocer según
T = 2π
ω
ω = r × v → = 2π
R0v
siendo v la velocidad orbital. Entonces
mVía Lactea =
4π2R30
G
(
2π
R0v
)2
Ahora sólo hay que reemplazar los números y hacer la cuenta.
Los resueltos de
ALF
N
o
 h
ay
 p
ro
bl
em
a Wi
lly!
←click me!
8 Si ves un error, avisale a Alf , 201812062058 CC0
https://losresueltosdealf.wordpress.com/
https://losresueltosdealf.wordpress.com/contacto/
https://en.wikipedia.org/wiki/Public_domain
GUÍA 2 - FLUIDOS CLÁSICOS
Guía 2 - Fluidos clásicos
Problema 2
Ítem b
Para u = 0 la distribución de Maxwell-Boltzmann es
f (r,p) = n (r)
(2πmkT (r))3/2
exp
(
− p
2
2mkT (r)
)
Como se puede ver tiene simetría esférica en p, es decir que f (r,p) = f (r, |p |) = f (r, p). Entonces el número de partículas
con momento entre p y p+ dp se puede expresar según
n (r, p) =
2πˆ
θ=0
π̂
φ=0
f (r,p) d3r d3p
d3p = p2 sinφdθ dφ dp =
2πˆ
θ=0
π̂
φ=0
f (r,p) p2 sinφdθ dφ d3r dp
f = f (r, p)→ = p2f (r, p) d3r
2πˆ
θ=0
dθ
π̂
φ=0
sinφdφ dp
= 4πf (r, p) p2 d3r dp
= 4πn (r)
(2πmkT (r))3/2
p2 exp
(
− p
2
2mkT (r)
)
d3r dp
= f (r, p) dr dp
Por otro lado se sabe que el número de partículas con momento entre p1 y p2 viene dado por
n (r, p ∈ (p1, p2)) =
p2ˆ
p1
fp
(
r, p′
)
dp′ d3r
siendo fp (r, p) = f (r, p) la distribución hallada previamente. Usando la relación p = mv lo anterior puede expresarse como
n (r, v ∈ (v1, v2)) =
mv2ˆ
mv1
fp
(
r, p′
)
dp′ d3r
=
v2ˆ
v1
fp
(
r,mv′
) dp′
dv′︸ ︷︷ ︸
fv(r,v′)
dv′ d3r
por lo tanto
fv (r, v) = fp (r,mv)
dp
dv
fp (r, p) =
4πn (r)
(2πmkT (r))3/2
p2e−
p2
2mkT (r)
p = mv
→ = 4π
(
m
2πkT (r)
)3/2
v2 exp
(
− mv
2
2kT (r)
)
Ítem c
Cada una de estas cantidades es 
〈v〉 (r) def=
∞̂
0
vf (r, v) dv =
√
8kT (r)
πm
n (r)
〈
v2
〉
(r) def=
∞̂
0
v2f (r, v) dv = 3kT (r)n (r)
m
Los resueltos de
ALF
N
o
 h
ay
 p
ro
bl
em
a Wi
lly!
←click me!
9 Si ves un error, avisale a Alf , 201812062058 CC0
https://losresueltosdealf.wordpress.com/
https://losresueltosdealf.wordpress.com/contacto/
https://en.wikipedia.org/wiki/Public_domain
Problema 4 GUÍA 2 - FLUIDOS CLÁSICOS
Problema 4
Ítem a
La fotósfera del sol está compuesta en un 75 % de hidrógeno. Una línea espectral de frecuencia νemitida emitida por un
átomo que se mueve con velocidad β = vc
νobservada = νemitida (1− β · r̂) γ→ Doppler shift
= νemitida (1− β cos θ) γ
siendo r̂ el versor que apunta en la dirección de observación y θ el ángulo entre r̂ y β. Debido a que la distribución de
velocidades es isótropa entonces θ tiene una distribución uniforme entre 0 y 2π. En cambio β sigue la distribución de
Maxwell-Boltzmann.
Por algún motivo que no termino de entender hay que considerar que θ = 0 o θ = π. Esto me parece medio trucho pero
es la forma en que usualmente se trata el problema (ver por ejemplo este link). Eligiendo θ = π queda
νobservada = νemitida (1 + β) γ
Si la velocidad de las partículas es v � c entonces el factor γ ≈ 1 y se puede simplificar
νobservada = νemitida
(
1 + v
c
)
De aquí se puede despejar
v = c
(
νobservada
νemitida
− 1
)
→ Cambio de variable
Considérese ahora la distribución de Maxwell-Boltzmann
fv (v) =
( m
2πkT
)3/2
exp
(
−mv
2
2kT
)
que en una única dimensión (digamos vx) es
fv (v) =
√
m
2πkT exp
(
−mv
2
2kT
)
→ En una dimensión
Se puede hacer ahora un cambio de variable para ponerla como función de la frecuencia
fν (ν) = fv (v (ν))
dv
dν
ν = νobservada → = 4π
√
m
2πkT
c3
ν0
(
ν
ν0
− 1
)2
exp
(
−mc
2
2kT
(
ν
ν0
− 1
)2)
Ahora se puede calcular el ensanchamiento Doppler según
Ensanchamiento→σν =
√
〈ν2〉 − 〈ν〉2
...
=
√
kT
mc2
ν0 �
Las cuentas las hice con el programa wxMaxima, a continuación dejo el código y los resultados:
Los resueltos de
ALF
N
o
 h
ay
 p
ro
bl
em
a Wi
lly!
←click me!
10 Si ves un error, avisale a Alf , 201812062058 CC0
https://en.wikipedia.org/wiki/Photosphere#Composition_of_the_Sun
https://es.wikipedia.org/wiki/Ensanchamiento_Doppler
https://es.wikipedia.org/wiki/Ensanchamiento_Doppler#Derivaci%C3%B3n
https://wxmaxima-developers.github.io/wxmaxima/
https://losresueltosdealf.wordpress.com/
https://losresueltosdealf.wordpress.com/contacto/
https://en.wikipedia.org/wiki/Public_domain
Problema 5 GUÍA 2 - FLUIDOS CLÁSICOS
(%i5) assume(m>0)$ assume(ν_0>0)$ assume(T>0)$ assume(c>0)$ assume(k>0)$
(%i6) f_v(v) := sqrt(m/(2*%pi*k*T))*exp(-m*v^2/(2*k*T)); /* Maxwell-Boltzmann */
(%o6) fv( )v :=
m
2 π k T
exp
( )−m v2
2 k T
(%i7) v(ν) := c*(ν/ν_0 - 1); /* Doppler */
(%o7) v( )ν :=c
ν
ν0
−1
(%i8) f_ν(ν) := f_v(v(ν))*diff(v(ν), ν, 1); /* Cambio de variable */
(%o8) fν( )ν :=fv( )v( )ν
d
d ν
v( )ν
(%i9) σ_ν: sqrt( integrate(f_ν(ν)*ν^2, ν, -inf, inf)- integrate(f_ν(ν)*ν, ν, -inf, inf)^2 );
(σ_ν)
c m
23 2/ π T k ν0
3
c m
+
23 2/ π T3 2/ k3 2/ ν0
3
c3 m3 2/
23 2/ π T k ν0
−ν0
2
(%i10) factor(σ_ν);
(%o10)
T k ν0
c m
Problema 5
De acuerdo a lo que se vio en la clase teórica el tensor de presiones es
Pij (r) =
1
m
ˆ
pipjf (r,p) d3p→ Tensor de presiones
Si hay equilibrio termodinámico local entonces f es la distribución de Maxwell-Boltzmann fp (r,p) = n(r)(2πkT (r))3/2 exp
(
− p
2
2mkT (r)
)
.
Haciendo las cuentas (me las hizo la compu) se puede verificar que
Pij (r) =
{
n (r) kT (r)m3/2 para i = j
0 para i 6= j
Entonces
p = 13 tr (Pij)
= mnkT
Me está sobrando la masa m, no sé... Como se puede ver se satisface que para un fluido ideal en equilibrio termodinámico
Pij = pδij→ Para fluido ideal en eq.
Problema 7
El fenómeno será mas intenso en el ecuador. Si la estrella tiene un radio R entonces una partícula de masa m en el ecuador
sentirá las fuerzas F gravitatoria = −
mMG
R2
r̂
F centrífuga = mω2Rr̂
En el caso límite ambas fuerzas son iguales con lo cual se obtiene
ω0 =
√
MG
R3
Utilizando datos de Wikipedia obtengo
ωpara el sol0 ∼ 6,6× 10−4
lo que corresponde a un período de
T para el sol0 ∼ 0,11 dı́as
El período de rotación real del sol es
Tsol ∼ 24,5 dı́as
lo cual está bastante lejos del valor límite.
Los resueltos de
ALF
N
o
 h
ay
 p
ro
bl
em
a Wi
lly!
←click me!
11 Si ves un error, avisale a Alf , 201812062058 CC0
https://losresueltosdealf.wordpress.com/
https://losresueltosdealf.wordpress.com/contacto/
https://en.wikipedia.org/wiki/Public_domain
Problema 8 GUÍA 2 - FLUIDOS CLÁSICOS
Problema 8
Ítem a
Las únicas fuerzas a las que estarán sometidos los elementos de fluido, vistos desde un sistema de referencia inercial (o
sea que no rota con el planeta) serán la fuerza gravitatoria y la fuerza centrífuga.
La densidad de fuerza centrífuga viene dada por
f centrífuga = −ρω × (ω × r) → Densidad de fuerza centrífuga
Para ω = ωẑ la fuerza centrífuga será (coordenadas cilíndricas)
f centrífuga = −ρω2ẑ × (ẑ × r)
ẑ × r =
∣∣∣∣∣∣
r̂ θ̂ ẑ
0 0 1
r 0 z
∣∣∣∣∣∣→ = −ρω2rẑ × θ̂
ẑ × θ̂ = −r̂ → = ρω2rr̂de cilíndricas
Como se puede ver
f centrífuga =
∂
∂r
(
ρr2ω2
2
)
r̂
= −∇
(
−ρω
2
2 r
2
)
→ Coordenadas cilíndricas
por lo tanto
φcentrífugo = −
ρω2
2 r
2→ Coordenadas cilíndricas
Por otro lado la fuerza gravitatoria se obtiene a partir del potencial gravitatorio. Si uno quisiera hacer bien este problema
debería escribir
φgravitatorio = −ρG
˚
Vdel planeta
d3r′∣∣ r − r′ ∣∣
y, como Vdel planeta es función de la fuerza centrífuga y la fuerza centrífuga es función de la deformación del planeta, se
obtendría un sistema acoplado e imposible de resolver manualmente. CREO que hay que asumir que el planeta no se deforma
demasiado, entonces aún se lo puede considerar como una esfera en lo que respecta al φgravitatorio. Para una distribución
esférica el potencial gravitatorio en el interior es (lo saqué de este link)
φgravitatorio ≈ dφgravitatoriocel interior de una esfera
= −2π3 ρ
2Gr2→ Coordenadas esféricas
Entonces
φ = φcentrífugo + φgravitatorio
= −ρω
2
2 r
2
de cilíndricas −
2π
3 ρ
2Gr2
rde cilíncricas = r sinφ→ = −
ρω2
2 r
2 sin2 φ− 2π3 ρ
2Gr2
donde r ≡ rde esféricas.
Ítem b
Si la superficie del planeta es una isobara entonces
d∇pcSuprficie del planeta = 0
Debido a que el planeta está compuesto por un fluido ideal entonces la dinámica de este queda descripta por la ecuación de
Navier-Stokes
∂tu+ (u ·∇)u = −
∇p
ρ
+ f→ Navier-Stokes
Los resueltos de
ALF
N
o
 h
ay
 p
ro
bl
em
a Wi
lly!
←click me!
12 Si ves un error, avisale a Alf , 201812062058 CC0
https://es.wikipedia.org/wiki/Fuerza_centr%C3%ADfuga
https://en.wikipedia.org/wiki/Gravitational_potential#Spherical_symmetry
https://losresueltosdealf.wordpress.com/
https://losresueltosdealf.wordpress.com/contacto/
https://en.wikipedia.org/wiki/Public_domain
Problema 9 GUÍA 2 - FLUIDOS CLÁSICOS
En el caso estacionario se tiene que ∂tu = 0. Por otro lado el término (u ·∇)u se anula debido a la simetría del problema:
es de esperar que
u = u (r, φ) θ̂→ Por la simetría
Entonces
(u ·∇)u = u (r, θ) 1
r sinφ
∂
∂θ︸ ︷︷ ︸
θ̂·∇
u (r, φ) θ̂ = 0
En consecuencia la ecuación de Navier-Stokes queda
∇p = ρf
f = −∇φ→ = −ρ∇φ
Si d∇pcSuprficie del planeta = 0 entonces
dφcSuperficie del planeta = φ0→ Constante
−ρω
2
2 R
2 sin2 φ− 2π3 ρ
2GR2 =
de donde se puede despejar
R =
√
− φ0
ρω2
2 sin
2 φ+ 2π3 ρ2G
Para dar un valor de φ0 se puede considerar un punto de referencia como por ejemplo un polo. En este caso
dφcEn un polo = φ (R0, φ = 0)
= −2π3 ρ
2GR20
Entonces
R (φ) = R0
√
1
3ω2
4πρG sin
2 φ+ 1
Como se puede ver R (θ) ≥ R0 = dRcEn el polo lo cual tiene sentido X.
Problema 9
Ítem a
La ecuación de Navier-Stokes para un fluido ideal es(
∂
∂t
+ u ·∇
)
u = −∇p
ρ
+ f→ Navier-Stokes
Creo que hay que asumir que el fluido está en reposo, esto es
u ≡ 0→ Fluido en reposo
(que es distinto del caso estacionario ∂u∂t = 0). La densidad de fuerza por unidad de masa f será producto del potencial
gravitatorio, es decir que
f = −∇φ
con φ = MGr siendo M la masa del planeta. La ecuación de Navier-Stokes queda entonces
∇p
ρ
+ ∇φ = 0
Considerando a la atmósfera como formada por un gas ideal entonces vale la ecuación de estado del gas ideal
p = nkT
= ρ
m
kT
Introduciendo esto en la ecuación de Navier-Stokes, y usando la hipótesis de atmósfera isotérmica,
kT
m
∇ρ
ρ
+ ∇φ = 0
Los resueltos de
ALF
N
o
 h
ay
 p
ro
bl
em
a Wi
lly!
←click me!
13 Si ves un error, avisale a Alf , 201812062058 CC0
https://losresueltosdealf.wordpress.com/
https://losresueltosdealf.wordpress.com/contacto/
https://en.wikipedia.org/wiki/Public_domain
Problema 9 GUÍA 2 - FLUIDOS CLÁSICOS
Reemplazando ∇ en coordenadas esféricas se obtienen las siguientes tres ecuaciones
r̂) kT
m
1
ρ
∂ρ
∂r
− MG
r2
= 0
θ̂) ∂ρ
∂θ
= 0
φ̂) ∂ρ
∂φ
= 0
Evidentemente ρ sólo depende de r, como era de esperarse dada la simetría del problema. Para la ecuación de la componente
radial se tiene la ecuación diferencial
∂ρ
∂r
− mMG
kTr2
ρ = 0
cuya solución es
ρ (r) = C exp
(
−mMG
kTr
)
Asumiendo que ρ (R) = ρR (la densidad a nivel del suelo) entonces
ρ (r) = ρR exp
(
−mMG
kT
[
1
R
− 1
r
])
En cuanto a la presión, ésta se obtiene simplemente reemplazando ρ en la ecuación de estado p = ρmkT .
Para r ≈ R podemos considerar una aproximación de Taylor para el exponente, es decir
1
r
− 1
R
≈ −r −R
R2
+ . . .
y si además definimos z def= r −R entonces
ρ (z) ≈ ρR exp
(
−zmMG
kTR2
)
MG
R2
def= g → = ρR exp
(
−zmg
kT
)
H
def= kT
mg
→ = ρR exp
(
− z
H
)
Obsérvese que también podríamos haber aproximado a ρ del siguiente modo
ρ (r ∼ R) ≈ ρR
(
1− mGM
R2kT
(r −R) + . . .
)
g
def= GM
R2
z
def= r −R
→ = ρR (1− mgkT z + . . .)
= ρR
(
1− z
H
+ . . .
)
Ítem b
Es hacer la integralita
´∞
0 ρRe
− zH dz = ρRH.
Los resueltos de
ALF
N
o
 h
ay
 p
ro
bl
em
a Wi
lly!
←click me!
14 Si ves un error, avisale a Alf , 201812062058 CC0
https://losresueltosdealf.wordpress.com/
https://losresueltosdealf.wordpress.com/contacto/
https://en.wikipedia.org/wiki/Public_domain
GUÍA 3 - RELATIVIDAD
Guía 3 - Relatividad
Problema 2 (No entiendo enunciado)
Ítem a
No termino de entender qué está pasando...
θ
β βap
Si fuese como en el dibujito que hice entonces βap = β sin θ y listo...
Problema 3
Ítem a
La diferencia más notable que encuentro es que para θ = π2 en el caso clásico no hay efecto Doppler mientras que en el
caso relativista sí hay.
Ítem b
A partir de
{
ν′ = ν (1− β · r̂) γ
λν = c
se obtiene
λ′ = λ(1− β · r̂) γ
y entonces
z =
λ
(1−β·r̂)γ − λ
λ
= 1− (1− β · r̂) γ(1− β · r̂) γ
Reemplazando β · r̂ = vc cos θ y γ =
1√
1− v2
c2
se puede despejar
v
c
=
(z + 1)2 cos θ ±
√
(z + 1)2 cos2 θ − z2 − 2z
(z + 1)2 cos2 θ + 1
Para θ = 0 lo anterior se reduce a
v
c
= z (z + 2)
(z + 1)2 + 1
Dada la definición de z def= λ
′−λ0
λ0
se tiene que z ∈ (−1,∞). En consecuencia −1 < vc < 1 lo cual demuestra que | v | < c.
Si existieran velocidadessuperlumínicas lo que ocurriría es que z se volvería imaginario, como así también la frecuencia
observada y muchas otras cantidades. No sé cómo se interpreta esto.
Problema 4
Ítem a
La ley que se aplica para sacar todo es la conservación del cuadrimomento
dpµcπ+ = dp
µcµ+ + νµ
= dpµcµ+ + dp
µcνµ
Para cada partícula se tiene 
dpµcπ+ =
[
Eπ+/c
0
]
→ Pión en reposo
dpµcµ+ =
[
Eµ/c
pµ+
]
dpµcνµ =
[
Eν/c
pνµ
]
Los resueltos de
ALF
N
o
 h
ay
 p
ro
bl
em
a Wi
lly!
←click me!
15 Si ves un error, avisale a Alf , 201812062058 CC0
https://losresueltosdealf.wordpress.com/
https://losresueltosdealf.wordpress.com/contacto/
https://en.wikipedia.org/wiki/Public_domain
Problema 5 GUÍA 3 - RELATIVIDAD
por lo tanto {
Eπ+ = Eµ+ + Eνµ
pµ+ = −pνµ ≡ p
Utilizando la relación
E2 = c2p2 +m20c4
que se obtiene de considerar el invariante relativista pµpµ del siguiente modo
dpµpµcFrame cualquiera = dp
µpµcFrame en el que p=0 y γ=1
pµ =
[
E
c
p
]
→ E
2
c2
− p2 = m2c2← pµ = mγ
[
c
p
]
se encuentra que 
Eπ+ = mπ+c2
Eµ+ =
√
m2µ+c
4 + c2p2
Eνµ = cp
Reemplazando en Eπ+ = Eµ+ + Eνµ se obtiene
mπ+c
2 =
√
m2µ+c
4 + c2p2 + cp
de donde se puede despejar
p = c
m2π+ −m2µ+
2mπ+
Entonces
Eπ+ = Eµ+ + Eνµ
mπ+c
2 =
√
m2µ+c
4 + c2p2 + Eνµ
y reemplazando p = c
m2
π+−m
2
µ+
2mπ+
se obtiene finalmente
Eνµ = c2
m2µ+ + 3m2π+
2mπ+
Problema 5
Ítem a
En dibujitos lo que está pasando es lo siguiente
Electrón θ
ϕ
Situación inicial Situación final
Fotón
Los momentos iniciales y finales son
(
pinicialγ
)
µ =
[
hνi
c
hνi
c x̂
]
(
piniciale
)
µ =
[
mc2
0
]

(
pfinalγ
)
µ =
[
hνf
c
hνf
c (x̂ cos θ + ŷ sin θ)
]
(
pfinale
)
µ =
[ √
p2 +mc2
p (x̂ cosφ− ŷ sinφ)
]
donde se ha usado la relación E2 = c2p2 +m2c4 para las componentes temporales de los cuadrimomentos. Por conservación
de cuadrimomento vale que
Conservación→

hνi
c
+mc2 = hνf
c
+
√
p2 +m2c4 → Componente temporal
hνi
c
= hνf
c
cos θ + p cosφ → Componente en x̂
0 = hνf
c
sin θ − p sinφ → Componente en ŷ
Los resueltos de
ALF
N
o
 h
ay
 p
ro
bl
em
a Wi
lly!
←click me!
16 Si ves un error, avisale a Alf , 201812062058 CC0
https://losresueltosdealf.wordpress.com/
https://losresueltosdealf.wordpress.com/contacto/
https://en.wikipedia.org/wiki/Public_domain
Problema 6 (averiguar si esto está bien) GUÍA 3 - RELATIVIDAD
De las segundas dos ecuaciones se puede despejar

p cosφ = h
c
(νi − νf cos θ)
p sinφ = h
c
νf sin θ
. Elevando cada una al cuadrado y sumándolas
se obtiene
p2 = h
2
c2
(
ν2i + ν2f − 2νiνf cos θ
)
De la primera ecuación (conservación de componente temporal) se puede despejar
p2 =
[
h
c
(νi − νf ) +mc2
]2
−m2c4
= h
2
c2
(νi − νf )2 + 2
h
c
(νi − νf )mc2 +���m2c4 −���m2c4
= h
2
c2
(
ν2i + ν2f − 2νiνf
)
+ 2h
c
(νi − νf )mc2
Igualando estas dos últimas expresiones supongo que se despeja el resultado...
h�2
c�2
(
���
�ν2i + ν2f − 2νiνf cos θ
)
= h
�2
c�2
(
���
�ν2i + ν2f − 2νiνf
)
+ 2
�
��h
c
(νi − νf )mc2
−2h
c
νiνf cos θ = −2
h
c
νiνf + 2 (νi − νf )mc2
h
mc3
(1− cos θ) = νi − νf
νiνf
Listo, a menos de un c y un signo menos este es el resultado X.
Problema 6 (averiguar si esto está bien)
La ecuación de estado de un gas ideal es
PV = nRT → Gases ideales
siendo R = kNA la constante universal de los gases ideales (k es la constante de Boltzmann y NA el número de Avogadro).
Ahora es simplemente hacer la cuenta. En particular, haciendo las cuentas con los datos de la consigna, se tiene que
P =
{
c1n
5/3 Caso no relativista
c2n
4/3 Caso relativista
con c1 = 3
2/3h2
20π2/3m y c2 =
31/3ch
8π1/3 . Entonces las ecuaciones de estado son
Eqs. de estado→
{
c1n
4/3V = RT Caso no relativista
c2nV = RT Caso relativista
Para el que esté interesado en las cuentas (hechas con el programa wxMaxima):
(%i2) ε_clasico: p^2/2/m;
ε_relativista: p·c;
(ε_clasico)
p2
2 m
(ε_relativista) c p
(%i4) P_clasica: 8·%pi/3/h^3·integrate(p^3·diff(ε_clasico,p,1), p, 0, p_F);
P_relativista: 8·%pi/3/h^3·integrate(p^3·diff(ε_relativista,p,1), p, 0, p_F);
(P_clasica)
8 π pF
5
15 h3 m
(P_relativista)
2 π c pF
4
3 h3
(%i6) P_clasica: subst((3·h^3/8/%pi·n)^(1/3), p_F, P_clasica);
P_relativista: subst((3·h^3/8/%pi·n)^(1/3), p_F, P_relativista);
(P_clasica)
32 3/ h2 n5 3/
20 π2 3/ m
(P_relativista)
31 3/ c h n4 3/
8 π1 3/
Los resueltos de
ALF
N
o
 h
ay
 p
ro
bl
em
a Wi
lly!
←click me!
17 Si ves un error, avisale a Alf , 201812062058 CC0
http://andrejv.github.io/wxmaxima/
https://losresueltosdealf.wordpress.com/
https://losresueltosdealf.wordpress.com/contacto/
https://en.wikipedia.org/wiki/Public_domain
Problema 7 GUÍA 3 - RELATIVIDAD
Problema 7
Si un fotón fue emitido con una frecuencia νemitido y fue detectado con una frecuencia νobservado entonces el redshift se
calcula según
zgravitatorio
def= νemitido − νobservado
νobservado
→ Redshift
ν ∝ E → = ∆E
E
Redshift gravitatorio→ = −∆φ
c2
Ahora se considerará que 
dφcTierra = 0
dφcDonde se emite el fotón = −
MG
R
La justificación para considerar que dφcTierra = 0 radica en el hecho de que
MTierra
RTierra
� MR , es decir que considero que el fotón
se emitió en un objeto muy masivo y compacto. Entonces
zgravitatorio =
MG
Rc2
Por otro lado el redshift Doppler es tal que νobservado = νemitido (1− β · r̂) γ. Asumiendo que β = βr̂ entonces
zDoppler = 1− γ + βγ
Igualando ambos redshifts y despejando para β (recordar que γ = 1√
1−β2
) se obtiene
β =

GM
(
2Rc2 −GM
)
Rc2 (2Rc2 − 2GM) +G2M2
1
Como β = 1 implica que el objeto se mueve a la velocidad de la luz, descartamos esta solución y nos quedamos con la otra.
Los resueltos de
ALF
N
o
 h
ay
 p
ro
bl
em
a Wi
lly!
←click me!
18 Si ves un error, avisale a Alf , 201812062058 CC0
https://en.wikipedia.org/wiki/Redshift#Measurement,_characterization,_and_interpretation
https://losresueltosdealf.wordpress.com/
https://losresueltosdealf.wordpress.com/contacto/
https://en.wikipedia.org/wiki/Public_domain
GUÍA 4 - ATMÓSFERAS ESTELARES
Guía 4 - Atmósferas estelares
Problema 1
Ítem a (falta una cosa!)
El número de partículas es por definición
nν =
ˆ
dr
ˆ
dpfν (r,p)
En cuanto a la energía, si hν es la energía de un fotón con frecuencia ν entonces
�ν = hνnν
es la energía de n fotones de frecuencia ν.
En astronomía se suele definir la intensidad de radiación Iν de modo tal que
dEnergía = Iν dΩ dν dS dt
Es decir que
Iν ∼
Energía
Ángulo sólido× Frecuencia× Superficie× Tiempo
De acuerdo a lo que vimos en la práctica
Iν (r) =
h4ν3
c2
ˆ
fν (r,p) d3p
Está bien lo anterior?
La intensidad media Jν es es el momento de orden cero
Jν (r, t)
def= 14π
‹
Iν (r, t) dΩ
mientras que el flujo es el momento de primer orden
qν
def= 1
c
‹
uIν dΩ
Ítem c
La componente ij del tensor de presión es
(pν)ij =
1
c3
‹
Jνuiuj dΩ→ Tensor de presión
Problema 3
Ítem a
Por definición
qν
def= 1
c
‹
uIν dΩ
Usando que todo es ideal y que la radiación se “propaga hacia afuera”, esto es u = cr̂
u
entonces
qν = 0
para cualquier superficie esférica que rodee a la estrella... Algo no anda bien. Pregunté en clase y la respuesta fue difusa.
La definición de luminosidad de una estrella es
Lν
def=
‹
qν · ds
Los resueltos de
ALF
N
o
 h
ay
 p
ro
bl
em
a Wi
lly!
←click me!
19 Si ves un error, avisale a Alf , 201812062058 CC0
https://losresueltosdealf.wordpress.com/
https://losresueltosdealf.wordpress.com/contacto/
https://en.wikipedia.org/wiki/Public_domain
Problema 6 GUÍA 4 - ATMÓSFERAS ESTELARES
Asumiendo que
qν = qν (r) r̂→ Asumo esto
entonces Lν se puede calcular en forma trivial sobre cualquier superficie esférica de radio r centrada en la estrella y da
Lν = 4πr2qν (r)
Debido a que r puede ser cualquiera, en particular se puede elegir r = R (o sea la superficie de la estrella) con lo cual
Lν =
{
4πr2qν (r)
4πR2qν (R) ← Evaluado en la superficie
y entonces
qν (r) = qν (R)
R2
r2
Problema 6
Ítem a
Dado un fluido (gas de fotones) la fuerza que éste ejerce sobre una superficie ds = n̂ ds es
dFi = σijdsj→ Lo vi en E1
siendo σij el tensor de esfuerzos. En el caso de un fluido ideal
σij = −p1→ Fluido ideal
con lo cual
dFi = −pδij dsj
= −p dsi
y ésta será la fuerza que se ejerce sobre un elemento de fluido de hidrógeno ionizado.
Para un gas de fotones se tiene que
pν
def= 13 tr
(
(pν)ij
)
(pν)ij =
1
c3
‹
Jνuiuj dΩ→ =
1
3 tr
(
1
c3
‹
Jνuiuj dΩ
)
Isotropía por consigna⇒ Jν = Iν → =
Iν
3c3
‹
tr (uiuj) dΩ
= Iν3c3
‹
|u |2 dΩ
u = cr̂ → = 4πIν3c
Ahora usamos las siguientes formulitas{
qν (R) = 4πIν → Flujo en superficie estrella
Lν = 4πR2qν (R) → Luminosidad
con lo cual
pν =
Lν
12πR2c
Integrando en todo el espectro
p =
ˆ
pν dν
= L12πR2c
Entonces
dF = − L12πR2c ds
Para concluir se debe integrar esta expresión en todo el contorno de un elemento de H II, el cual consideraré como un
objeto sólido con uno contorno del siguiente modo
Los resueltos de
ALF
N
o
 h
ay
 p
ro
bl
em
a Wi
lly!
←click me!
20 Si ves un error, avisale a Alf , 201812062058 CC0
https://losresueltosdealf.wordpress.com/
https://losresueltosdealf.wordpress.com/contacto/
https://en.wikipedia.org/wiki/Public_domain
Problema 7 GUÍA 4 - ATMÓSFERAS ESTELARES
R
R+dr
p(R)
p(R+dR)Elemento de H II
ds S1 S2
S3
S4
R
dθ
Entonces
F =
˛
S
dF
Se cancelan por simetría→ =
�
�
��
ˆ
S1
dF +
ˆ
S2
dF +
�
�
��
ˆ
S3
dF +
ˆ
S4
dF
= L12πc
ˆ
S2
ds
R2
+
ˆ
S4
ds
R2

≈ L12πc
(
S2
R2
+ S4
(R+ dr)2
)
S2 = −S2r̂
S4 = S4r̂
}
→ = Lr̂12πc
(
S4
(R+ dr)2
− S2
R2
)
Taylor a orden cero en dr → ≈ Lr̂12πc
(
S4
R2
− S2
R2
)
= (S4 − S2)L12πcR2 r̂
Con un poco de fe decimos ahora que S4 − S2 = σ que es la sección eficaz y entonces
F = σL12πcR2 r̂
Obsérvese que como S4 > S2 entonces σ > 0 y la fuerza es, efectivamente, “hacia afuera” de la estrella.
Problema 7
Ítem a
La ecuación de transporte de radiación es
∂tIν + u ·∇Iν = (Sν − Iν) ρκν

Sν =
eν + σνJν
κν
→ Función fuente
κν = kν + σν → Opacidad
kν → Coeficiente de absorción
σν → Coeficiente de scattering
ρ → Densidad de absorbentes
eν → Emisividad
Cuando se está en equilibrio termodinámico entonces
En equilibrio→

Iν = Bν → Distribución de Maxwell-Boltzmann
∂tIν = 0 → Estacionario
∇Iν = 0 → Homogeneidad
Jν = Iν → Isotropía
con lo cual la ecuación de transporte se convierte en
0 =
(
eν + σνBν
κν
−Bν
)
ρκν
Los resueltos de
ALF
N
o
 h
ay
 p
ro
bl
em
a Wi
lly!
←click me!
21 Si ves un error, avisale a Alf , 201812062058 CC0
https://losresueltosdealf.wordpress.com/
https://losresueltosdealf.wordpress.com/contacto/
https://en.wikipedia.org/wiki/Public_domain
Problema 8 GUÍA 4 - ATMÓSFERAS ESTELARES
por lo tanto eν+σνBνκν −Bν = 0 de donde se obtiene
eν = (κν − σν)Bν
κ
def= k + σ → = kνBν X
Problema 8
Ítem a
Para una atmósfera estacionaria, plana y estratificada en la que z es la dirección vertical vale que
∂
∂t
≡ 0 → Estacionaria
∂
∂x
= ∂
∂y
≡ 0 → Plana
? → Estratificada
con lo cual la ecuación de transporte ∂tIν + u ·∇Iν = (Sν − Iν) ρκν resulta
c cos θ
ρκν
∂Iν
∂z
= Sν − Iν
donde se usó que u = cû con û apuntando en la dirección del observador. Definiendo

µ
def= cos θ
τν
def=
∞̂
z
ρκν
c
dz
finalmente se obtiene
µ
∂Iν
∂τν
= Iν − Sν
Ítem b
De la ecuación de transporte se puede despejar
Sν = µ
∂Iν
∂τν
− Iν
Ahora se multiplica a ambos lados por e−
τν
µ con lo cual
e−
τν
µ Sν =
(
µ
∂Iν
∂τν
− Iν
)
e−
τν
µ
= µ ∂
∂τν
(
Iνe
− τνµ
)
Multiplicando por ∂τν a ambos lados e integrando se obtiene
El infinito (?)ˆ
El observador
e−
τν
µ Sν dτν = µ
El infinito (?)ˆ
El observador
d
(
Iνe
− τνµ
)
= µ
⌈
Iνe
− τνµ
⌋El infinito (?)
En el observador
El observador ≡ τν = 0→ = µIν (τν = 0)
por lo tanto
Iν (τν = 0) =
1
µ
∞̂
0
e−
τν
µ Sν dτν
No sé por qué el límite de integración superior es ∞, así lo hicimos en clase.
Ítem c
Asumiendo que Sν ≈ aτν + b entonces (sólo es hacer la cuentita)
Iν (τν = 0) = µa+ b
= a cos θ + b
Como a, b > 0 entonces esta es una función que decrece a medida que θ → π2 que es justamente lo que ocurre en el limbo
solar.
Los resueltos de
ALF
N
o
 h
ay
 p
ro
bl
em
a Wi
lly!
←click me!
22 Si ves un error, avisale a Alf , 201812062058 CC0
https://losresueltosdealf.wordpress.com/
https://losresueltosdealf.wordpress.com/contacto/
https://en.wikipedia.org/wiki/Public_domain
Problema 9 FALTA GUÍA 4 - ATMÓSFERAS ESTELARES
Problema 9 FALTA
Hay que modelarlo como un problema de atmósfera plano-paralela pero “al revés” que en el caso del sol. Es decir, tenemos
lo siguiente:
θ
Atmósfera
← Mr. Eyes
← Estrella en el infinito
Intensidad de radiación
justo afuera de la atmósfera, m0
z
u
La propagación en la atmósfera está modelada por la ecuación de transporte
∂tIν + u ·∇Iν = (Sν − Iν) ρκν→ Transport equation
Al igual que en el problema 8, consideramos las hipótesis de atmósfera estacionaria, plana y estratificada con lo cual
∂
∂t
= 0
∂
∂x
= ∂
∂y
= 0
Considerando además que
u = cû = c (x̂ sin θ − ẑ cos θ)
la ecuación de transporte queda
c cos θ∂Iν
∂z
= (Sν − Iν) ρκν
o, en términos de la profundidad óptica τν
def=
´∞
z
ρκν
c dz,
∂Iν
∂τν
cos θ = Iν − Sν
.............. DANDOAS NDKAOS PDNKASOP DNASKOP DNAPOD
Los resueltos de
ALF
N
o
 h
ay
 p
ro
bl
em
a Wi
lly!
←click me!
23 Si ves un error, avisale a Alf , 201812062058 CC0
https://losresueltosdealf.wordpress.com/
https://losresueltosdealf.wordpress.com/contacto/
https://en.wikipedia.org/wiki/Public_domain
GUÍA 5 - MEDIO INTERESTELAR
Guía 5 - Medio interestelar
Problema 2
La energía de una partícula de polvo puede expresarse según
E = Ede traslación + Ede rotación
= m2 v
2 + I2ω
2
donde I es el momento de inercia y ω la velocidad angular. En consecuencia el valor medio de la energía será
〈E〉 = m2
〈
v2
〉
+ I2
〈
ω2
〉
= m2
(〈
v2x
〉
+
〈
v2y
〉
+
〈
v2z
〉)
+ I2
(〈
ω2x
〉
+
〈
ω2y
〉
+
〈
ω2z
〉)
= 〈Tx〉+ 〈Ty〉+ 〈Tz〉+ 〈Rx〉+ 〈Ry〉+ 〈Rz〉
donde Ti y Ri son las energías de traslación y de rotación en cada eje, respectivamente.
Por otro lado, el teorema de equipartición de la energía dice que si cada partícula de un gas ideal contiene i ∈ {1, . . . , N}
grados de libertad, entonces a cada grado de libertad le corresponde una energía media dada por
〈�i〉 =
kT
2 → Equipartición
En el caso de los granos de polvo esto quiere decir que
〈Ti〉 =
kT
2
〈Ri〉 =
kT
2
De aquí se puede despejar 
〈
v2i
〉
= kT
m〈
ω2i
〉
= kT
I
y entonces 
〈
v2
〉
= 3kT
m〈
ω2
〉
= 3kT
I
Problema 3
Ítem a
Si se consideran dos diferenciales de masa dentro de esta esfera de fluido, la energía potencial entre ellos será
dE = − dm1 dm2G
| r1 − r2 |
dm = ρ d3r → = −ρ (r1) ρ (r2)G
| r1 − r2 |
d3r1 d
3r2
La energía potencial total para toda la distribución será
E = −G
ˆ
d3r1
ˆ
d3r2
ρ (r1) ρ (r2)
| r1 − r2 |
Multiplicando y dividiendo por los parámetros apropiados se obtiene
E = ΘGM
2
R
con
Θ = − R
M2
ˆ
d3r1
ˆ
d3r2
ρ (r1) ρ (r2)
| r1 − r2 |
Los resueltos de
ALF
N
o
 h
ay
 p
ro
bl
em
a Wi
lly!
←click me!
24 Si ves un error, avisale a Alf , 201812062058 CC0
https://losresueltosdealf.wordpress.com/
https://losresueltosdealf.wordpress.com/contacto/
https://en.wikipedia.org/wiki/Public_domain
Problema 3 GUÍA 5 - MEDIO INTERESTELAR
Si la distribución es homogénea entonces ρ (r) = ρ = constante con lo cual
Θ = −Rρ
2
M2
ˆ ˆ
d3r1 d
3r2
| r1 − r2 |
→ Para distribución homogénea
r
def= r1 − r2
d3r = d3r1
}
→ = −Rρ
2
M2
ˆ
d3r2
ˆ
d3r
r
d3r en esféricas→ = −Rρ
2
M2
ˆ
d3r2 4π
R̂
0
r dr
= −2πR
3ρ2
M2
ˆ
d3r2︸ ︷︷ ︸
4π
3 R
3
= −8π
2
3
R6ρ2
M2
M = 4π3 ρR
3 → = −32
En este caso la energía es simplemente
E = −32
GM2
R
(En clase quedó 52 ...)
Ítem b PREGUNTAR
La ecuación de Navier-Stokes, luego de aplicar todas las simplificaciones de siempre, es
∇p+ ρ∇φ = 0→ Navier-Stokes
Ahora hay que usar un truco que es multiplicar por ·r e integrar. Con eso sale. El primer término es
ˆ
∇p · r d3r =
Identidad vectorial︷ ︸︸ ︷ˆ
∇ · (pr) d3r −
ˆ
p∇ · r d3r
Teorema de la div.
∇ · r = 3
}
→ =
˛
pr · ds− 3
ˆ
p d3r
dpc¸ ≡ p0 → = p0
˛
r · ds−3
ˆ
p d3r
= 4πR3p0 − 4
ˆ
p d3r
Ahora asumimos equilibrio termodinámico y entonces p = nkT . Además, se asume todo homogéneo menos n (r). Entonces
ˆ
∇p · r d3r = 4πR3p0 − 4kT
ˆ
n (r) d3r︸ ︷︷ ︸
N
= 4πR3p0 − 4kTN
cV (r) =
3
2
k
m
para gas ideal monoat.→ = 4πR3p0 −
8
3cV Tsocmsocmsocm
donde NA es el número de Avogadro. EN CLASE QUEDÓ DISTINTO! Pero no sé qué se usó...
El otro término de Navier-Stokes es
ˆ
∇φ · r d3r = −ΘGM2
ˆ
d3r
r2
← φ = ΘGM
2
R
⇒∇φ = −ΘGM
2
R2
Θ = −32ˆ
d3r
r2
= 4πR
→ = 6πGM2R
Los resueltos de
ALF
N
o
 h
ay
 p
ro
bl
em
a Wi
lly!
←click me!
25 Si ves un error, avisale a Alf , 201812062058 CC0
https://losresueltosdealf.wordpress.com/
https://losresueltosdealf.wordpress.com/contacto/
https://en.wikipedia.org/wiki/Public_domain
Problema 3 GUÍA 5 - MEDIO INTERESTELAR
Juntando los dos términos
4πR3p0 −
8
3cV T
N
NA
+ 6πρGM2R = 0
reemplazando ρ = 4π3
M
R3 , y despejando p0 se obtiene algo que no es lo de la consigna... Debe haber alguna cuenta mal en
algún lugar, ya fue.
Los resueltos de
ALF
N
o
 h
ay
 p
ro
bl
em
a Wi
lly!
←click me!
26 Si ves un error, avisale a Alf , 201812062058 CC0
https://losresueltosdealf.wordpress.com/
https://losresueltosdealf.wordpress.com/contacto/
https://en.wikipedia.org/wiki/Public_domain
GUÍA 6 - ESTRUCTURA Y EVOLUCIÓN ESTELAR
Guía 6 - Estructura y evolución estelar
Problema 1
Ítem b
PREGUNTA 1 - Escala hidrodinámica
No la estoy encontrando en la bibliografía.
Escala de Kelvin-Helmholtz [Carroll and Ostlie, 2017, p. 386] La escala de Kelvin-Helmholtz representa el tiempo de vida
de una estrella asumiendo que su única fuente de energía es la gravedad. Es decir, pensemos en una estrella que mediante
algún proceso mágico convierte energía gravitatoria en radiación, y que su única fuente de energía es la gravitatoria. En este
caso el tiempo de vida de la estrella sería, asumiendo además que su luminosidad es constante,
τKelvin-Helmholtz =
Egravitatoria
L
= GM
2
LR
.
Escala nuclear [Carroll and Ostlie, 2017, p. 388] La escala nuclear representa el tiempo de vida de una estrella asumiendo
que su única fuente de energía es la conversión de materia de los núcleos mediante E = mc2. Entonces, asumiendo luminosidad
constante a lo largo de su vida,
τnuclear ∼
Mc2
L
.
Lo anterior sería válido si toda la masa de la estrella se pudiera convertir en radiación. Sin embargo las reacciones nucleares
permiten que esto ocurra sólo para una fracción de la masa. Es por ello que se añade un factor Q que tiene en cuenta el
“rendimiento” de este proceso (es decir Q ∝ mproductosmreactivos ) y además tiene en cuenta que no toda la estrella tiene temperatura
suficiente para realizar reacciones nucleares (por ejemplo sólo el 10 % de la masa en el núcleo está sometida a condiciones
tales para que se den las reacciones [Carroll and Ostlie, 2017, example 3.2, p. 389]). Entonces
τnuclear =
QMc2
L
.
Problema 2
La expresión
B = c2 (Zmp + (A− Z)mn −m) → Energía de ligadura del núcleo
es la energía de ligadura del núcleo. La energía de ligadura por nucleón es, entonces,
f = B
A
.→ Energía de ligadura por nucleón
El gráfico de f como función de A es una cosa así [Choudhuri, 2010, fig. 4.1], [Carroll and Ostlie, 2017, p. 404]
Los resueltos de
ALF
N
o
 h
ay
 p
ro
bl
em
a Wi
lly!
←click me!
27 Si ves un error, avisale a Alf , 201812062058 CC0
https://losresueltosdealf.wordpress.com/
https://losresueltosdealf.wordpress.com/contacto/
https://en.wikipedia.org/wiki/Public_domain
Problema 5 GUÍA 6 - ESTRUCTURA Y EVOLUCIÓN ESTELAR
y como se puede ver el elemento que está en la cima es el 56Fe. Esto quiere decir que si se tienen, por ejemplo, dos átomos
de hidrógeno y se los combina en un helio, éste último será más estable que los dos hidrógenos. Si se continúa fusionando
núcleos éstos son cada vez más estables hasta llegar al hierro a partir del cual comienzan a hacerse más inestables. Entonces,
en forma natural, lo que ocurrirá es lo siguiente:
Fusión es más
estable
Fisión es más
estable
Problema 5
Ítem a
Para encontrar la ecuación de equilibrio hidrostático de una estrella, se considera la ecuación de Navier-Stokes para un
fluido con viscosidad nula:
∂tu+ (u ·∇)u = −
∇p
ρ
+ f .→ Navier-Stokes
En la ecuación de equilibrio hidrostático se asumen condiciones estáticas (obvio, ja) por lo tanto
u = 0.→ Estático
Dada la simetría esférica del problema es de esperar que p = p (r) por lo que ∇p = ∂p∂r r̂. En cuanto a la fuerza, ésta es
fgravitatoria = −MrGr2 r̂ donde Mr es la masa encerrada por la esfera de radio r. Entonces [Carroll and Ostlie, 2017, p. 377]
∂p
∂r
= ρMrG
r2
.
Ítem b
Si ui es la energía interna por unidad de masa entonces
Ei =
ˆ
ui dM
es la energía interna total. Por otro lado se tiene la relación
M (r) = 4π
rˆ
0
ρ
(
r′
)
r′2 dr′
por lo tanto
dM = 4πρ (r) r2 dr
y entonces, realizando el cambio de variable en la integral de Ei,
Ei = 4π
R̂
0
uρr2 dr
Según consigna u = 3p
ρn
→ = 12π
n
R̂
0
pr2 dr.
Los resueltos de
ALF
N
o
 h
ay
 p
ro
bl
em
a Wi
lly!
←click me!
28 Si ves un error, avisale a Alf , 201812062058 CC0
https://losresueltosdealf.wordpress.com/
https://losresueltosdealf.wordpress.com/contacto/
https://en.wikipedia.org/wiki/Public_domain
Problema 5 GUÍA 6 - ESTRUCTURA Y EVOLUCIÓN ESTELAR
PREGUNTA 2 - Energía potencial gravitatoria?
Cómo hago la cuenta de Eg?
Ítem c
En equilibrio hidrodinámico vale el teorema del virial (?) por lo tanto
Etotal = Einterna + Egravitatoria
Teorema del virial→ = Einterna − nEinterna
Gas ideal⇒ Einterna =
3
2kT → =
3
2kT (1− n) .
Ahora el calor específico es
cg =
dE
dT
= 32k (1− n) .
En ningún lugar se dice nada sobre n más que “es una constante”. Supongo que n > 1, de lo contrario el resultado no es el
esperado.
PREGUNTA 3 - Qué es el teorema del virial?
Por qué vale en equilibrio hidrodinámico?
Los resueltos de
ALF
N
o
 h
ay
 p
ro
bl
em
a Wi
lly!
←click me!
29 Si ves un error, avisale a Alf , 201812062058 CC0
https://losresueltosdealf.wordpress.com/
https://losresueltosdealf.wordpress.com/contacto/
https://en.wikipedia.org/wiki/Public_domain
GUÍA 7 - CAMPOS MAGNÉTICOS
Guía 7 - Campos magnéticos
Problema 1
Ítem a
La radiación de un dipolo alineado con el eje ẑ es (Física 3 o Teórica 1)
B (r, θ) = µ04π
[
− 1
r3
ẑ + 3 cos θ
r3
r̂
]
|m |
donde m es el momento dipolar magnético. De acuerdo con la consigna éste se relaciona con otras cantidades del problema
según
|m | =
√
Ė
12πc2
ω4
√
ε0
µ0
→ Consigna
por lo tanto el campo magnético será
B (r, θ) = µ04π
[
− 1
r3
ẑ + 3 cos θ
r3
r̂
]√
Ė
12πc2
ω4
√
ε0
µ0
.
Ahora hay que encontrar Ė en función de ω y listo. Para un cuerpo rígido la energía de rotación es E = I2ω
2 por lo que
Ė = Iωω̇.
Reemplazando este resultado en la expresión del campo magnético se encuentra que
B (r, θ) = µ04π
[
− 1
r3
ẑ + 3 cos θ
r3
r̂
]√
12πIc2 ω̇
ω3
√
ε0
µ0
.
Finalmente, usando que
ω = 2π
P
→ Consigna
y que ω̇ = − 2πP 2 Ṗ entonces
ω̇
ω3 = −
ṖP
4π4 y el campo magnético es
B (r, θ) = µ04π
[
− 1
r3
ẑ + 3 cos θ
r3
r̂
]√
− 62π Ic
2ṖP
√
ε0
µ0
.
Me da imaginario, con lo cual tiene que estar mal. Cuando lo hicieron en clase dijeron “no ponemos los signos negativos” y
eso fue todo... No sé si pifié en las cuentas o si hubo magia negra en clase.
Problema 3 - MHD
Las ecuaciones MHD son, de acuerdo a lo que vimos en la teórica y en la práctica,
Eqs. MHD→

∂tρ+ ∇ · (ρu) = 0 → Continuidad
ρ (∂t + u ·∇)u = −∇p+
1
c
J ×B + ρcargaE → Navier-Stokes
E + 1
c
u×B = 1
σ
J → Ley de Ohm
∂tB = ∇× (u×B) + η∇2B → Eq. de inducción
c∇×B = 4πJ + ∂tE → La agregaron en la práctica
Por algún motivo que ahora desconozco hay que asumir que
Hay que asumir esto→

| ∂tE |
|J |
� 1
Rm
def= |∇× (u×B) |
| η∇2B |
� 1
.
Los resueltos de
ALF
N
o
 h
ay
 p
ro
bl
em
a Wi
lly!
←click me!
30 Si ves un error, avisale a Alf , 201812062058 CC0
https://losresueltosdealf.wordpress.com/
https://losresueltosdealf.wordpress.com/contacto/
https://en.wikipedia.org/wiki/Public_domainProblema 3 - MHD GUÍA 7 - CAMPOS MAGNÉTICOS
No sé de dónde sale eso, pero en clase lo asumieron. Entonces esto implica que en las ecuaciones MHD se cancela lo que está
indicado en rojo: 
∂tρ+ ∇ · (ρu) = 0 → Continuidad
ρ (∂t + u ·∇)u = −��∇p+
1
c
J ×B +����ρcargaE → Navier-Stokes
E + 1
c
u×B = 1
σ
J → Ley de Ohm
∂tB = ∇× (u×B) +���η∇2B → Eq. de inducción
c∇×B = 4πJ +��∂tE → La agregaron en la práctica
y lo que se canceló en azul es un completo misterio para mí. Es lo que hicieron en clase. Aparentemente el término con
ρcargaE se cancela pues | ρcargaE ||J×B | � 1, pero ni idea por qué.
Por otro lado la consigna dice que u0 = 0, el fluido es incompresible y homogéneo, el campo magnético es uniforme, y
hay algo que es ideal (no sé qué). Entonces
Hipótesis de la consigna→

u = δu
ρ = ρ0
B0 = B0ẑ
p = p0 + δp
J = J0 + δJ
E = E0 + δE
p
ργ
= p0
ργ0
→ Es por eso de ideal
con ∇ρ0 = 0, ∇×B0 = 0. Al meter esto en las ecuaciones MHD se obtiene
Eqs. MHD→

∇ · δu = 0
ρ0∂tδu =
1
c
(J0 × δB + δJ ×B0ẑ)
δE +
1
c
δu ×B0ẑ =
1
σ
δJ
∂tδB = (B0 ·∇) δu −B0 (∇ · δu)
0 = δJ
.
Sabiendo que ∇ · δu = 0 y que δJ = 0 se puede simplificar lo anterior tal que
∇ · δu = 0
ρ0∂tδu =
1
c
J0 × δB
δE +
1
c
δu ×B0ẑ = 0
∂tδB = (B0 ·∇) δu
.
Proponiendo δ ∼ ei(k·r−ωt) lo anterior es 
k · δu = 0
− iωρ0δu −
1
c
J0 × δB = 0
δE +
1
c
δu ×B0ẑ = 0
ωδB + (B0ẑ · k) δu = 0
Ya fue, esto jamás va a dar lo que pide la consigna. Ni siquiera tengo bien en claro cuáles son las ecuaciones MHD, en la
teórica son distintas que en la práctica. En clase fue pura magia negra la resolución de este problema.
Ítem b
Debido a que ω ∝ kz entonces las ondas de Alfvén son no dispersivas. Además como
{
k · δu = 0
k · δB = 0
entonces son ondas
transversales.
Los resueltos de
ALF
N
o
 h
ay
 p
ro
bl
em
a Wi
lly!
←click me!
31 Si ves un error, avisale a Alf , 201812062058 CC0
https://losresueltosdealf.wordpress.com/
https://losresueltosdealf.wordpress.com/contacto/
https://en.wikipedia.org/wiki/Public_domain
Problema 4 - Reconexión de Sweet-Parker GUÍA 7 - CAMPOS MAGNÉTICOS
Problema 4 - Reconexión de Sweet-Parker
Ítem a
En el modelo de Sweet-Parker se asume
1. Una condición estacionaria, es decir que ∂t ≡ 0.
2. Fluido homogéneo ⇒∇ρ = 0.
3. Se asume una geometría idealizada como sigue
uin
uin
uoutuout L
Lminúscula
4. Se asume que L� `.
5. Se considera que uin � vAlfvén ⇒ ρinu2in �
B2in
4π .
6. ρoutu2out �
B2out
4π
?⇒ uout � vAlfvén.
En [Choudhuri, 2010, fig. 8.9] se ilustra este fenómeno:
Ítem b
Entiendo que el problema se resuelve con las ecuaciones MHD. Voy a seguir lo que hicimos en la teórica.
Si la situación es estacionaria entonces la ecuación de continuidad es ∇·(ρu) = 0. Si además se asume fluido incompresible
entonces ∇ · u = 0 e integrando sobre el contorno rectangular idealizado que se dibujó previamente se encuentra fácilmente
que
uin
uout
= `
L
.
Consideremos ahora la ecuación de Navier-Stokes
ρ (∂t + u ·∇)u = −∇p+
1
c
J ×B + ρcargaE.→ Navier-Stokes
En el caso estacionario se anula el término ∂tu y además vale que J = c4π∇ × B pues una de las ecuaciones MHD es
c∇ ×B = 4πJ + ���*
Estacionario
∂tE . Además el término ρcargaE lo despreciamos, no sé por qué. Entonces la ecuación de Navier-Stokes
queda
ρ (u ·∇)u = −∇p+ 14π (∇×B)×B
= −∇p− 14πB × (∇×B)
A× (B ×C) ≡ B (A ·C)− (A ·B)C → = −∇p− 14π (∇ (B ·B)− (B ·∇)B)
= −∇p+ 14π (B ·∇)B −∇
(
B2
4π
)
.
Los resueltos de
ALF
N
o
 h
ay
 p
ro
bl
em
a Wi
lly!
←click me!
32 Si ves un error, avisale a Alf , 201812062058 CC0
https://losresueltosdealf.wordpress.com/
https://losresueltosdealf.wordpress.com/contacto/
https://en.wikipedia.org/wiki/Public_domain
Problema 4 - Reconexión de Sweet-Parker GUÍA 7 - CAMPOS MAGNÉTICOS
Del lado izquierdo se puede aplicar la misma identidad pero en “sentido inverso” de modo tal que (u ·∇)u = ∇u2 − u ×
(∇× u) y entonces
ρ (∇× u)× u+ ρ∇u2 = −∇p+ 14π (B ·∇)B −∇
(
B2
4π
)
.
A continuación se sacará provecho de la simetría de reflexión que tiene el problema en torno a los siguientes ejes
uin
uin
uoutuout
sobre los cuales vale, por simetría, lo siguiente
Sobre el eje rojo→

u ∝ ŷ
B ∝ x̂
∇× u = 0
∂xB = 0
Sobre el eje azul→

u ∝ x̂
B ∝ ŷ
∇× u = 0
∂yB = 0
por lo tanto la componente ŷ en el eje rojo y la componente x̂ en el eje azul son
dComponente ŷ de la ecuaciónceje rojo →
∂
∂y
(
p+ ρu2 + B
2
4π
)
= 0
dComponente x̂ de la ecuaciónceje azul →
∂
∂x
(
p+ ρu2 + B
2
4π
)
= 0
lo cual implica que 
⌈
p+ ρu2 + B
2
4π
⌋
eje rojo
= constante⌈
p+ ρu2 + B
2
4π
⌋
eje azul
= otra constante
.
Evalúese ahora cada una de estas expresiones en los puntos indicados a continuación:
uin
uin
uoutuout
A
B
C
Obtiénese ⌈
p+ ρu2 + B
2
4π
⌋
Punto A
=
⌈
p+ ρu2 + B
2
4π
⌋
Punto B
=
⌈
p+ ρu2 + B
2
4π
⌋
Punto C
y, nuevamente por simetría, {
ducPunto A = 0
dBcPunto A = 0
por lo tanto
dpcA = dpcB + ρu
2
in +
⌈
B2
⌋
B
4π = dpcC + ρu
2
out +
⌈
B2
⌋
C
4π .
Los resueltos de
ALF
N
o
 h
ay
 p
ro
bl
em
a Wi
lly!
←click me!
33 Si ves un error, avisale a Alf , 201812062058 CC0
https://losresueltosdealf.wordpress.com/
https://losresueltosdealf.wordpress.com/contacto/
https://en.wikipedia.org/wiki/Public_domain
Problema 5 GUÍA 7 - CAMPOS MAGNÉTICOS
(Recordar que ∇ρ = 0.) Ahora vamos a considerar que
Asumimos esto→uin � vAlfvén =
dBcB√
4πρ
por lo tanto ρu2in �
dB2cB
4π . En cuanto al punto C vamos a asumir “lo opuesto”, vamos a pensar que
Asumimos esto→ ρu2out �
⌈
B2
⌋
C
4π .
Entonces la ecuación anterior es
dpcA = pin +
B2in
4π = pout + ρu
2
out.
Por algún motivo que ahora no estoy viendo lo anterior implica que
uout =
Bin√
4πρ
.
O sea, no sé qué pasa con las presiones pero parece que son iguales. Aquí hemos encontrado que
uout = dvAcin
y del principio de todo ya sabíamos que uinuout =
`
L lo cual implica que
uin
dvAcin
= `
L
tal como pedía la consigna.
Problema 5
Ítem a
Recordemos que la densidad de energía electromagnética es [Jackson, 1999, eq. (6.106)]
u = 12 (E ·D +B ·H) .
En el vacío
{
D = E
H = B
y considerando además que E ≈ 0 entonces
u = B
2
2 .
La energía total será
Etotal = 2Ede cada tubo
= 2
ˆ
tubo
B2
2 dV
Se asume todo uniforme→ = πR2LB20 .
Ítem b
La “energía procesada por unidad de tiempo” es la potencia del proceso
P = Etotal
τ
.
La estimación de τ es bien cabeza: consiste en asumir dos cilindros que se superponen
τ ∼ 2R
Uin
.
Entonces
P ∼ πR
2LB20Uin
2R .
Los resueltos de
ALF
N
o
 h
ay
 p
ro
bl
em
a Wi
lly!
←click me!
34 Si ves un error, avisale a Alf , 201812062058 CC0
https://losresueltosdealf.wordpress.com/
https://losresueltosdealf.wordpress.com/contacto/
https://en.wikipedia.org/wiki/Public_domain
Problema 5 GUÍA 7 - CAMPOS MAGNÉTICOS
El número de Lundquist es
S−
1
2 = Uin
vA
→ Lundquist
vA =
B0√
4πρ0
→ = Uin
√
4πρ0
B0
donde vA es la velocidad de Alfvén.
Los resueltos de
ALF
N
o
 h
ay
 p
ro
bl
em
a Wi
lly!
←click me!
35 Si ves un error, avisale a Alf , 201812062058 CC0
https://losresueltosdealf.wordpress.com/
https://losresueltosdealf.wordpress.com/contacto/
https://en.wikipedia.org/wiki/Public_domain
GUÍA 8 - GALAXIAS
Guía 8 - Galaxias
Problema 1 (LSR)
El estándar de reposo local (LSR) se define como aquel sistema de referencia (no inercial) en el que el “elemento de fluido
en el que está el sol” está en reposo. Nótese que esto no significa que el sol está en reposo en el LSR. Las coordenadas del
LSR se definen en forma estandarizada según el siguiente esquema:
Posición
Velocidad
La velocidad de una estrella i con respecto al sol (o sea, lo que nosotros medimos) es
v i
Sol
= v i
LSR
− v LSR
Sol
donde v i
Sol
es la velocidad de la estrella medida con respecto al sol, v i
LSR
es la velocidad medida con respecto al LSR y v LSR
Sol
es la velocidad del LSR con respecto al sol. (Voy a usar esta notación v objeto
respecto a
, que podría parecer exagerada,por el único
motivo de que toda la dificultad de este tema es ser prolijo con estos subíndices. En el fondo no es más que un problema de
cinemática, lo que aprendimos en las primeras dos semanas de F1.) Véase que la relación previa es equivalente a
v Sol
LSR
= v i
Sol
− v i
LSR
.
Ahora se considera un promedio sobre muchas estrellas cercanas al sol tal que
v Sol
LSR
= 1
| {Estrellas vecinas al sol} |
∑
i∈{Estrellas vecinas al sol}
(
v i
Sol
− v i
LSR
)
=
〈
v i
Sol
〉
−
〈
v i
LSR
〉
.
Por definición de LSR se tiene que
1
N
N∑
i=1
v i
Cualquier punto
= v LSR
cualquier punto
→ Definición de LSR
lo cual implica que
〈
v i
LSR
〉
= 1
N
N∑
i=1
v i
LSR
Definición de LSR→ = v LSR
LSR
≡ 0
y entonces
v Sol
LSR
=
〈
v i
Sol
〉
.→ Velocidad peculiar del Sol
Es decir, la velocidad del sol con respecto al LSR (a veces llamada velocidad peculiar del sol) es el promedio de la velocidad
de las estrellas vecinas con respecto al sol.
Para medir v i
Sol
se puede utilizar un sistema de coordenadas esféricas centrado en el sol. La componente radial se mide
por efecto Doppler en tanto que las componentes angulares se miden por simple observación de la posición de la estrella y
conocimiento de su distancia respecto al sol.
Problema 2
Consideremos el siguiente dibujito (sacado de este link)
Los resueltos de
ALF
N
o
 h
ay
 p
ro
bl
em
a Wi
lly!
←click me!
36 Si ves un error, avisale a Alf , 201812062058 CC0
https://en.wikipedia.org/wiki/Local_standard_of_rest
https://en.wikipedia.org/wiki/Peculiar_velocity
https://en.wikipedia.org/wiki/Oort_constants
https://losresueltosdealf.wordpress.com/
https://losresueltosdealf.wordpress.com/contacto/
https://en.wikipedia.org/wiki/Public_domain
Problema 2 GUÍA 8 - GALAXIAS
Π0^
Θ0^
r
Θ
Π̂
Θ̂
r̂
Lo que se pide analizar son las componentes en φ̂ y r̂ de las estrellas cercanas medidas con respecto al LSR, es decir
mostrar que 〈
v i
LSR
〉
= r̂ rA sin (2l)∓ φ̂ r (B +A cos (2l)) .
donde el promedio debe entenderse como que se se realiza en una región acotada, por ejemplo como sigue:
r
← Acá se realiza el promedio
Recordando que v LSR
cualquier punto
def= 1N
∑N
i=1 v iCualquier punto véase que lo que se pide no es otra cosa que la velocidad del LSR
en un punto cercano al Sol medida con respecto al LSR solar, es decir〈
v i
LSR
〉
≡ v LSR(R,Θ)
LSR solar
.
Teniendo en mente lo anterior, procedamos ahora sí a resolver el problema. Si la galaxia es un disco con rotación diferencial
entonces
v LSR(R,Θ)
CG
= Ω (R)×R
Ω = Ω
(
−Ẑ
)
R = RΠ̂
→ = Ω (R)R(−Ẑ)× Π̂
= Ω (R)RΘ̂.
Ahora considérese la relación
v LSR(R,Θ)
CG
= v LSR(R,Θ)
LSR solar
+ v LSR solar
CG
por lo que
v LSR(R,Θ)
LSR solar
= v LSR(R,Θ)
CG
− v LSR solar
CG
v LSR(R,Θ)
CG
= Ω (R)RΘ̂
v LSR solar
CG
= Ω (R0)R0Θ̂0
→ = Ω (R)RΘ̂− Ω (R0)R0Θ̂0
Proyecto en r̂ y φ̂→ = Ω (R)R
(
r̂ cosα+ φ̂ sinα
)
− Ω (R0)R0
(
r̂ cos (π/2− `) + φ̂ sin (π/2− `)
)
= (Ω (R)R cosα− Ω (R0)R0 sin `) r̂ + (Ω (R)R sinα− Ω (R0)R0 cos `) φ̂.
Los resueltos de
ALF
N
o
 h
ay
 p
ro
bl
em
a Wi
lly!
←click me!
37 Si ves un error, avisale a Alf , 201812062058 CC0
https://losresueltosdealf.wordpress.com/
https://losresueltosdealf.wordpress.com/contacto/
https://en.wikipedia.org/wiki/Public_domain
Problema 3 GUÍA 8 - GALAXIAS
Ahora hay que hacer un poco de geometría para encontrar relaciones entre las cantidades de la ecuación anterior. En particular
véase que
Geometría pura→
{
R cosα = R0 sin `
R sinα = R0 cos `− r
por lo tanto
v LSR(R,Θ)
LSR solar
= (Ω (R)− Ω (R0))R0 sin ` r̂ + (Ω (R)− Ω (R0))R0 cos ` φ̂− Ω (R) r φ̂
= (Ω (R)− Ω (R0))R0
(
sin ` r̂ + cos ` φ̂
)
− Ω (R) r φ̂.
Ahora hay que hacer unas aproximaciones asumiendo
R0, R� r→ Asumimos esto
lo cual implica que
Ω (R)− Ω (R0) ≈
⌈
dΩ
dR
⌋
R0
(R−R0)︸ ︷︷ ︸
δR
Ω = V
R
→ =
⌈
− V
R2
+ 1
R
dV
dR
⌋
R=R0
(R−R0)
=
(
− V0
R20
+ 1
R0
⌈
dV
dR
⌋
R=R0
)
(R−R0)
R−R0 ≡ −r cos `→ =
(
V0
R20
− 1
R0
⌈
dV
dR
⌋
R=R0
)
r cos `
donde V y V0 son las velocidades que figuran en el primer diagramita de la resolución. Reemplazando esto en v LSR(R,Θ)
LSR solar
se
obtiene que
v LSR(R,Θ)
LSR solar
=
(
V0
R �20
− 1
��R0
⌈
dV
dR
⌋
R=R0
)
��R0
(
r cos ` sin ` r̂ + r cos ` cos ` φ̂
)
− Ω (R) r φ̂
Id. trig.→ =
(
V0
R0
−
⌈
dV
dR
⌋
R=R0
)
︸ ︷︷ ︸
2A
(
r
sin (2`)
2 r̂ + r
cos (2`) + 1
2 φ̂
)
− Ω (R) r φ̂
= Ar sin (2`) r̂ +Ar cos (2`) φ̂+ 12
(
V0
R0
−
⌈
dV
dR
⌋
R=R0
)
r φ̂− Ω (R) r φ̂.
Para la componente radial ya se tiene lo que se estaba buscando X. Para la componente tangencial (o sea φ̂) se puede
aproximar, usando que R,R0 � r,
Ω (R) ≈ Ω (R0) =
V0
R0
y entonces
v LSR(R,Θ)
LSR solar
≈ Ar sin (2`) r̂ +Ar cos (2`) φ̂+ 12
(
V0
R0
−
⌈
dV
dR
⌋
R=R0
)
r φ̂− V0
R0
r φ̂
= Ar sin (2`) r̂ +Ar cos (2`) φ̂− 12
(
V0
R0
+
⌈
dV
dR
⌋
R=R0
)
︸ ︷︷ ︸
B
r φ̂
= Ar sin (2`) r̂ + [A cos (2`)−B] r φ̂ X.
Problema 3
Ítem a
El lagrangiano de una partícula en un potencial central es
L = m2 ṙ
2 + mMG
r2
Los resueltos de
ALF
N
o
 h
ay
 p
ro
bl
em
a Wi
lly!
←click me!
38 Si ves un error, avisale a Alf , 201812062058 CC0
https://losresueltosdealf.wordpress.com/
https://losresueltosdealf.wordpress.com/contacto/
https://en.wikipedia.org/wiki/Public_domain
Problema 3 GUÍA 8 - GALAXIAS
donde r es el vector que une a la partícula y a la masa que produce el campo gravitatorio. En coordenadas polares
r = rr̂ ⇒ ṙ = ṙr̂ + rθ̇θ̂ ⇒ ṙ2 = ṙ2 + r2θ̇2
esto es
L = m2
(
ṙ2 + r2θ̇2
)
+ mMG
r
.
Las ecuaciones de Euler-Lagrange ddt
(
dL
dq̇
)
= dLdq para cada coordenada son
Euler-Lagrange→

mr̈ = d
dr
(
m
2 r
2θ̇2 + mMG
r
)
d
dt
(
r2θ̇
)
= 0
.
De la segunda ecuación tenemos que
L = r2θ̇→ Se conserva
es una cantidad conservada. Reemplazando en la primera
r̈ = d
dr
(
L2
2r2 +
MG
r
)
donde L
2
2r2 +
MG
r = φe es el potencial efectivo.
Ítem b
Se busca una solución de la forma
r (t) = r0 + r1 (t)
donde r0 es una solución estacionaria, i.e. corresponde a un punto de equilibrio del potencial efectivo.
El potencial efectivo se puede aproximar según
dφe
dr
=
⌈
dφe
dr
⌋
r0
+
⌈
d2φe
dr2
⌋
r0
(r − r0) + . . .
y debido a que r0 es un punto de equilibrio entonces ⌈
dφe
dr
⌋
r0
= 0.
Por otro lado, por definición, r − r0 = r1. Entonces la ecuación de movimiento para r1 es
r̈1 +
⌈
d2φe
dr2
⌋
r0
r1 = 0
(en algún lugar le chingué a un signo pues lo tuve que forzar el +) y entonces la frecuencia de oscilación es
ω2 =
⌈
d2φe
dr2
⌋
r0
.
Los resueltos de
ALF
N
o
 h
ay
 p
ro
bl
em
a Wi
lly!
←click me!
39 Si ves un error, avisale a Alf , 201812062058 CC0
https://losresueltosdealf.wordpress.com/
https://losresueltosdealf.wordpress.com/contacto/
https://en.wikipedia.org/wiki/Public_domain
GUÍA 9 - COSMOLOGÍA
Guía 9 - Cosmología
Problema 2
Ítem a
Aclaración La ley de Hubble no se deduce de la ecuación de continuidad, es un hecho observacional. A lo sumo se deduce
de la relatividad general y algún modelo cosmológico, pero eso excede ampliamente este curso. Lo que la consigna quiso decir
es “meta la ley de Hubble en la ecuación de continuidad y vea qué pasa”.
Box 2 - Ley de Hubble-Lemaître
La ley de Hubble-Lemaître (es un hecho observacional que) establece que
1. A escalas extragalácticas (10 Mpc o más) los objetos presentan un redshift que se interpreta como una
velocidad de recesión (o sea, que se alejan) respecto a La Tierra.
2. La velocidad de recesión obtenida mediante el redshift observado satisface
v ≈ H0r
donde v es la velocidad del objeto en cuestión con respecto a La Tierra y r es el vector que apunta desde La
Tierra hasta el objeto en cuestión.
Box 3 - El principio cosmológico
El principio cosmológico establece que a gran escala la densidad del universo es uniforme, es decir
∇ρ = 0 a gran escala.
La ecuación de continuidad es la de siempre
∂ρ
∂t
+ ∇ · (ρv) = 0.→ Continuidad para un fluido
Desarrollando ∇ · (ρv) = ∇ρ · v + ρ∇ · v y asumiendo la validez del principio cosmológico
∇ρ = 0→ Principio cosmológico
entonces la ecuación de continuidad es
∂ρ
∂t
+ ρ∇ · v = 0.
Siahora usamos la ley de Hubble-Lemaître
v = H (t) r→ Ley de Hubble-Lemaître
y usamos que ∇ · r ≡ 3 obtenemos finalmente
H = − 13ρ
∂ρ
∂t
.→ Esto es lo que pedía la consigna
Ítem b
Reemplazando H = ṘR en la ecuación hallada en el ítem previo se obtiene
1
R
∂R
∂t
= − 13ρ
∂ρ
∂t
⇒ ∂
∂t
(
ρR3
)
= 0 ⇒ ρR3 = constante.
Escribiendo a la constante como ρ0R30 se obtiene que
ρR3 = ρ0R30
y la interpretación de esto es que si R (t) representa un volumen entonces
M = ρ (t)R3 (t) → Masa no depende de t
es la masa contenida en dicho volumen, y como ρ y R evolucionan de modo tal que M es independiente de t, se dice que R
es el radio de una esfera que “se mueve con la materia”, o algo así. Como que R fuese el radio del universo es la idea.
Los resueltos de
ALF
N
o
 h
ay
 p
ro
bl
em
a Wi
lly!
←click me!
40 Si ves un error, avisale a Alf , 201812062058 CC0
https://en.wikipedia.org/wiki/Hubble%27s_law
https://en.wikipedia.org/wiki/Hubble%27s_law
https://en.wikipedia.org/wiki/Cosmological_principle
https://losresueltosdealf.wordpress.com/
https://losresueltosdealf.wordpress.com/contacto/
https://en.wikipedia.org/wiki/Public_domain
Problema 2 GUÍA 9 - COSMOLOGÍA
Ítem c
La ecuación de Euler es
∂v
∂t
+ (v ·∇)v = −
�
�
�7
Ppio. cosmológico
∇p
ρ
−∇φ→ Euler
donde φ es el potencial gravitatorio que satisface, al igual que el potencial electrostático, la ecuación de Poisson
∇2φ = 4πGρ.→ Poisson
Reemplazando la ley de Hubble v = H (t) r en la ecuación de Euler se obtiene
r
∂H
∂t
+H2 (r ·∇) r = −∇φ.
r
(
Ḣ +H2
)
= ← (r ·∇) r = r
Multiplicando por ∇· a ambos lados de la ecuación esto es
∇ · r
(
Ḣ +H2
)
= −∇2φ
y ahora usando ∇ · r ≡ 3 y ∇2φ = 4πGρ (Poisson) encontramos que
Ḣ +H2 = −4π3 Gρ.
Reemplazando H ≡ ṘR el término de la izquierda es
Ḣ +H2 = R̈
R
−
�
��Ṙ
2
R2︸ ︷︷ ︸
Ḣ
+
�
��Ṙ
2
R2
por lo tanto
R̈ = −4πG3 ρR.
Finalmente, usando ρ = ρ0R
3
0
R3 ,
R̈ = −4πGρ0R
3
0
3R2 . X
Ítem d
Para integrar se usa el siguiente truco
R̈
def= dṘ
dt
Chain rule→ = dṘ
dR
dR
dt
= dṘ
dR
Ṙ
= d
dR
Ṙ2
2 .
Reemplazando esto en la ecuación del ítem previo se obtiene
dṘ2
dR
= −8πGρ0R
3
0
3R2 ⇒ dṘ
2 = −8πGρ0R
3
0
3R2 dR ⇒ Ṙ
2 = 8πGρ0R
3
0
3R + c
o bien, definiendo la constante de integración negativa,
Ṙ2 = 8πGρ0R
3
0
3R − k.
El motivo por el cual sólo nos interesa k ∈ {0, 1,−1} es porque R se define a través de H = ṘR y como es evidente está
definido a menos de una constante multiplicativa. Esto hace que siempre se pueda redefinir R de modo tal que k ∈ {0, 1,−1}.
Los resueltos de
ALF
N
o
 h
ay
 p
ro
bl
em
a Wi
lly!
←click me!
41 Si ves un error, avisale a Alf , 201812062058 CC0
https://losresueltosdealf.wordpress.com/
https://losresueltosdealf.wordpress.com/contacto/
https://en.wikipedia.org/wiki/Public_domain
Problema 2 GUÍA 9 - COSMOLOGÍA
Ítem e
Del ítem previo se puede despejar ∣∣ Ṙ ∣∣ = √8πGρ0R303R − k.
Aquí se pueden distinguir dos regímenes:
1. Para k = − | k | (o sea k ∈ (−∞, 0]) se tiene que
∣∣ Ṙ ∣∣ = √8πGρ0R303R + | k |.
Acá se puede ver que la expansión (o contracción) del universo es eterna ya que si Ṙ > 0 en algún momento, siempre
lo será sin importar el valor de R.
2. Para k| k | = 1 (o sea k > 0) la ecuación anterior es
∣∣ Ṙ ∣∣ = √8πGρ0R303R − | k |.
En este caso existe un valor de R = Rmáx tal que Ṙ = 0, y es
Rmáx =
8πGρ0R30
3 | k | .
Para saber qué ocurre después, i.e. si se empieza a contraer, si se queda estanco o si se vuelve a expandir, hay que calcular
derivadas temporales de orden superior. Reemplazando el valor de R0 en la expresión encontrada con anterioridad para
R̈ se obtiene ⌈
R̈
⌋
Rmáx
= − 3 | k |
2
16πGρ0R30
y como
⌈
R̈
⌋
Rmáx
< 0 entonces efectivamente sigue una contracción.
Ítem f
Se define
Ω def= ρ3H2
8πG
.→ Parámetro de densidad
Usando ρ = ρ0R
3
0
R3
Ω = 8πGρ0R
3
0
3R
1
R2H2
Ṙ2 = 8πGρ0R
3
0
3R − k → =
Ṙ2 + k
R2H2
de donde se puede despejar
k = R2H2Ω− Ṙ2
H ≡ Ṙ
R
→ = Ṙ2 (Ω− 1) .
Como Ṙ2 > 0 entonces el signo de k depende de Ω ≷ 1.
Ítem g
De acuerdo a lo que vimos en clase la ecuación del calor es
3
2
k
m
ρ
(
∂T
∂t
+ v ·∇T
)
+ p∇ · v = 0→ Eq. del calor
donde k no estoy seguro si es la constante de Boltzmann u otra cosa, y m no sé qué es. El principio cosmológico (ver box 3)
establece que el universo es homogéneo, en consecuencia
∇T = 0.→ Cosmological principle
Los resueltos de
ALF
N
o
 h
ay
 p
ro
bl
em
a Wi
lly!
←click me!
42 Si ves un error, avisale a Alf , 201812062058 CC0
https://losresueltosdealf.wordpress.com/
https://losresueltosdealf.wordpress.com/contacto/
https://en.wikipedia.org/wiki/Public_domain
Problema 2 GUÍA 9 - COSMOLOGÍA
Además la ley de Hubble (ver box 2) dice que v = Hr. Entonces
1
2
k
m
ρ
∂T
∂t
+ pH = 0
donde se usó ∇ · r = 3. Ahora usamos la presión de radiación y la del gas idealpr =
4σ
3c T
4 → Radiation
pm = nkT → Matter
y además H ≡ ṘR y ρ = ρ0
(
R0
R
)3 por lo tanto la ecuación previa es
kρ0
2m
R30
R3
∂T
∂t
+ ξT ζ 1
R
∂R
∂t
= 0 ⇒ kρ0R
3
0
2mξ T
−ζ dT = −R2 dR
donde ξ y ζ son las constantes correspondientes a radiación o a materia. Integrando lo anterior se obtiene
kρ0R
3
0
2mξ
T 1−ζ
1− ζ = −
R3
3 if ζ 6= 1
kρ0R
3
0
2mξ lnT = −
R3
3 if ζ = 1
de donde se despeja
T ∼
{
R
3
1−ζ if ζ 6= 1
eR
3
if ζ = 1
.
En el caso de la presión de radiación ζ = 4 por lo tanto
T ∼ R−1→ Para radiación X
y para la presión del gas ideal ζ = 1 por lo tanto T ∼ eR3 jaja.
Los resueltos de
ALF
N
o
 h
ay
 p
ro
bl
em
a Wi
lly!
←click me!
43 Si ves un error, avisale a Alf , 201812062058 CC0
https://en.wikipedia.org/wiki/Radiation_pressure#Compression_in_a_uniform_radiation_field
https://losresueltosdealf.wordpress.com/
https://losresueltosdealf.wordpress.com/contacto/
https://en.wikipedia.org/wiki/Public_domain
FORMULAS
Formulas
Box 4 - Fórmulas matemáticas
∞̂
0
e−x dx = 1
∞̂
0
√
xe−x dx =
√
π
2
fx (x) distribución
x = g (y) cambio de variable
}
⇒ fy (y) = fx (g (y))
dg
dy
Box 5 - Fórmulas de la guía 1
Black body radiation→ I =

2hν3
c2
1
exp
(
hν
kT
)
− 1
→ En términos de ν
2hc2
λ5
1
exp
(
hc
λkT − 1
) → En términos de λ
E =
ˆ
I dν = σT 4→ Stephan-Boltzmann
L =
˛
F · ds→ Luminosidad en función del flujo
Leyes de Kepler Las leyes de Kepler son:
1. Las órbitas son elípticas, estando el sol en uno de los ejes.
2. El radio vector que une al planeta con el sol recorre áreas iguales en tiempos iguales.
3. Se satisface T 2 = 4π
2
G(m1+m2)a
3 donde T es el período orbital y a el semieje mayor de la elipse.
Método de paralaje Determina distancia de “objetos cercanos” midiendo el ángulo del siguiente dibujito:
θ
d
. Entonces tan θ = 1 AUd .
Box 6 - Fórmulas de fluidos clásicos
Maxwell-Boltzmann→

fp (r,p) =
n (r)
(2πkT (r))3/2
exp
(
− (p−mu (r))
2
2mkT (r)
)
fv (v) =
( m
2πkT
)3/2
exp
(
−mv
2
2kT
)
fv (r, |v |) = 4π
(
m
2πkT (r)
)3/2
v2 exp
(
− mv
2
2kT (r)
)
{
f centrífuga = −ρω × (ω ×R)
fCoriolis = −2ρω × u
∂tu+ (u ·∇)u = −
∇p
ρ
+ f→ Navier-Stokes p = nkT → Gas ideal
Pij (r) =
1
m
ˆ
pipjf (r,p) d3p→ Tensor de presiones
Los resueltos de
ALF
N
o
 h
ay
 p
ro
bl
em
a Wi
lly!
←click me!
44 Si ves un error, avisale a Alf , 201812062058 CC0
https://losresueltosdealf.wordpress.com/
https://losresueltosdealf.wordpress.com/contacto/
https://en.wikipedia.org/wiki/Public_domain
FORMULAS
Box 7 - Fórmulas relatividad especial
pµ = mγ
[
c
v
]
=
[
E/c
p
]
E2
c2
= p2 +m2c2
β =
v
c
γ =
(
1− β2
)−1/2
Redshift→ z =

λobservado − λemitido
λemitido
νemitido − νobservado
νobservado
νobservada = νemitida (1− β · r̂) γ→ Doppler
Box 8 - Fórmulas de transporte de radiación
dEnergía = Iν dΩ dν dS dt Iν (r) =
h4ν3
c2
ˆ
fν (r,p) d3p Jν (r, t)
def= 14π
‹
Iν dΩ qν
def=
?
1
c
‹
uIν dΩ
qν = qλ
∣∣∣∣ dλdν
∣∣∣∣ Lν def= ‹ qν · ds (pν)ij = 1c3
‹
Jνuiuj dΩ pν =
1
3 tr
(
(pν)ij
)
m =

m0 − 2,5 log
(
q
q0
)
m0 − 5 log
(r0
r
) → Magnitud aparente M =

m− 2,5 log
(dqc10 pc
q
)
m− 5 log
(
r
10 pc
) → Magnitud

Más contenidos de este tema