Hugoniot

Introducción Propiedades de la onda de choque Aplicaciones Propagación de perturbaciones en un gas en movimiento
Introducción a las Ondas de Choque
Luis Moraga
Centro de Física Experimental, Facultad de Ciencias, Universidad de Chile
Curso de Pre- y Postgrado ONDAS DE CHOQUE, 2008
Ondas de Choque
Asunto:
Introducción
La naturaleza de las ondas de choque:
Continuidad en las superficies de discontinuidad:
Adiabática de choque
Propiedades de la onda de choque
Ondas de choque en gases politrópicos
Ondas de choque fuertes
Atenuación de las ondas de choque
Aplicaciones
Frenado hipersónico de un objeto romo
Propagación de perturbaciones en un gas en movimiento
Características
Problemas resueltos
Ondas de Choque
Figura: Un tsunami.
Ondas de Choque
Figura: Ejemplo de una onda de choque
Ondas de Choque
Figura: Hugoniot y compañeros (1870)
Ondas de Choque
Figura: Formación de una onda de choque
Ondas de Choque
Figura: Pasaje de una onda sonora y de una onda de choque.
Ondas de Choque
Al paso de una onda sonora (1), el sistema oscila a lo largo de

una adiabática (curva de entropía constante). Por el contrario,
al paso de una onda de choque, el sistema salta de una
manera irreversible a lo largo de una Hugoniot hasta otro
estado (2) con mayor presión, densidad y entropía. Pasada la
onda de choque, el sistema relaja a lo largo de una nueva
adiabática hasta un estado (3) con la presión original, pero
menor densidad y mayor entropía y temperatura.
Ondas de Choque
¿ Qué cambia de un modo continuo a través de la superficie de

discontinuidad?
I El flujo de masa.
I El flujo de momentum.
I El flujo de energía.
Ondas de Choque
Onda de marea propagándose a lo largo de un canal rectangular:
v es la velocidad del agua detrás del frente de discontinuidad. u

es la velocidad del frente de discontinuidad.
Conservación de masa: h2 v∆t = (h2 − h1 )u∆t.

h1
v =u 1− .
h2
Ondas de Choque
Masa adicional que adquiere rapidez v = ρh2 u∆t − ρh2 v∆t

Fuerza neta hacia la derecha = 12 ρgh22 − 12 ρgh21
(ρh2 u∆t − ρh2 v∆t)v = ( 21 ρgh22 − 12 ρgh21 )∆t

h2
u2 = 21 g(h2 + h1 )
h1
√
Si h1 ' h2 , entonces u ' gh.
Ondas de Choque
Condiciones de continuidad:
Continuidad de la corriente de masa:
ρ1 v1 = ρ2 v2 = j. (1)
j = corriente de masa = densidad de momentum.

Ondas de Choque
Continuidad de la corriente de momentum:
p1 + ρ1 v12 = p2 + ρ2 v22 . (2)
Continuidad de la corriente de energía:
w1 + 21 v12 = w2 + 21 v22 . (3)
w = + pV
w: densidad de entalpía. : densidad de energía interna.
Ondas de Choque
Ecuaciones fundamentales:
p2 − p1
j2 = (Línea de Rayleigh.)
V1 − V2
1 − 2 + 12 (V1 − V2 )(p1 + p2 ) = 0
(Adiabática de choque o Hugoniot.)
V1 = 1/ρ1 , V2 = 1/ρ2 : Volúmenes específicos.
Ondas de Choque
Figura: p1 , V1 estados iniciales. p2 , V2 estados finales.
Ondas de Choque
Gases politrópicos
Un gas politrópico es uno en el cual la presión es inversamente
proporcional a una potencia (entera, fraccional o irracional) del
volumen:
Constante
p= .
Vγ
En este gas la curva de entropía constante (adiabática de
Poisson) es
pV γ = Otra constante.
La energía interna es
pV c2
= cV T = = .
γ−1 γ(γ − 1)
La entalpía es
γpV c2
w = cp T = = .
γ−1 γ−1
Ondas de Choque
En un gas politrópico la adiabática de choque (o de Hugoniot)
w1 − w2 + 12 (V1 + V2 )(p2 − p1 ) = 0
resulta ser
ρ1 V2 (γ + 1)p1 + (γ − 1)p2
= =
ρ2 V1 (γ − 1)p1 + (γ + 1)p2
(γ + 1)M12
=
(γ − 1)M12 + 2
donde M1 = v1 /c1 es el número de Mach.
Ondas de Choque
Adiabática de Hugoniot en un gas politrópico
Ondas de Choque
Relación entre las temperaturas a uno u otro lado de la

discontinuidad:
T2 p2 (γ + 1)p1 + (γ − 1)p2
=
T1 p1 (γ − 1)p1 + (γ + 1)p2
[2γM12 − (γ − 1)][(γ − 1)M12 + 2]
= .
(γ + 1)2 M1
Ondas de Choque
Relación entre las presiones a uno u otro lado de la

discontinuidad:
p2 2γM12 γ − 1
= − .
p1 γ+1 γ+1
Ondas de Choque
Relación entre los números de Mach (o las velocidades) a uno

u otro lado de la discontinuidad:
2 + (γ − 1)M12
M22 = ;
2γM12 − (γ − 1)
o (más simétricamente)
γM12 M22 − 12 (γ − 1)(M12 + M22 ) = 1.
Ondas de Choque
(γ − 1)p2 (γ + 1)p1
V2 ρ1 γ−1
= = .
V1 ρ2 γ+1
Note que el valor máximo de ρ2 es ρ2 = 4ρ1 para un gas
monoatómico y ρ2 = 6ρ1 para uno diatómico.
T2 γ − 1 p2
= .
T1 γ + 1 p1
En cambio, T2 puede ser arbitrariamente mayor que T1 .
Ondas de Choque
¿Dónde se produce la discontinuidad?
Figura: Formación de una onda de choque.
Ondas de Choque
Figura: Las áreas achuradas son iguales.
Z
(x − xs )dv = 0.
abcd
Ondas de Choque
En primera aproximación, hay una relación lineal entre la

rapidez del fluido y la del frente de choque u ' c0 + αv, con
c4 ∂ 2 V

α =
2v 2 ∂p2 S
1
= 2 (γ + 1)
para un gas politrópico.
Ondas de Choque
Área ABC = Área A0 DE
Ondas de Choque
Distancia entre B 0 y C 0 = α(∆v)1 t.

Tangente del ángulo DA0 E
(∆v)1 ∆v
= = .
`1 + αt(∆v)1 `
Ondas de Choque
Igualdad de áreas:
1 `2 (∆v)1
2 `1 (∆v)1 = 12 `(∆v) = .
2[`1 + αt(∆v)1 ]
Entonces: p
` = `1 1 + αt(∆v)1 /`1 (4)
Ondas de Choque
Además:
(∆v)1
∆v = p , (5)
1 + αt(∆v)1 /`1
Z
E = 21 ρ v 2 dx = 61 (∆v)2 `,
esto es,
E0
E=p . (6)
1 + αt(∆v)1 /`1
Ondas de Choque
Ondas planas
En conclusión, cuando t → ∞,
I La discontinuidad
√ en velocidad disminuye asintóticamente
√
como 1/ t (o, equivalentemente, 1/ x),
I La energía de una onda de
√ choque disminuye
asintóticamente como 1/ t.
√
I La longitud del pulso aumenta asintóticamente como t.
I La pendiente del perfil ∆v/` es, asintóticamente,
independiente del tiempo.
Ondas de Choque
Ondas cilíndricas.
Una onda cilíndrica es unapsolución de la ecuación de onda

que depende sólo de r = x2 + y 2 y t.
Ondas de Choque
La función ψ(x, t) = una constante es una posible onda plana.

La función ψ(r, t) = una constante no es una onda cilíndrica
posible.
Por lo tanto, al paso de una onda de choque cilíndrica, la zona
de compresión debe estar necesariamente seguida de otra de
rarefacción.
Figura: Perfil de una onda de choque cilíndrica.
Ondas de Choque
Ondas esféricas.
Una onda esférica esp

una solución de la ecuación de onda que
depende sólo de r = x2 + y 2 + z 2 y t.
Ondas de Choque
Atenuación de ondas de choque cilíndricas y esféricas.
Se concluye que, asintt́oticamente cuando r → ∞,

I Para una onda cilíndrica,
1
∆v ∝ .
r3/4
I Para una onda esférica,
1
∆v ∝ .
r log(r/a)
Ondas de Choque
Explorador de Marte
Ondas de Choque
Figura: Objeto romo en movimiento supersónico

Ondas de Choque
2γM12 − (γ − 1)
p2 = p1 ;
γ+1
2 + (γ − 1)M12
v2 = c1 ;
2 + (γ + 1)M1
(γ + 1)M12
ρ2 = ρ1 .
2 + (γ − 1)M12
Ondas de Choque
A lo largo de la línea de corriente AO, por Bernoulli:
w + 12 v 2 = w0 .
Para un gas politrópico

1/(γ−1) γ
c2

T ρ0
w= ; ρ = ρ0 ; p = p0 .
γ−1 T0 ρ
Ondas de Choque
T0
= 1 + 12 (γ − 1)M22 ;
T2
ρ0 1/(γ−1)
= 1 + 12 (γ − 1)M22 ;
ρ2
p0 γ/(γ−1)
= 1 + 12 (γ − 1)M22 ;
p2
Entonces,
(γ+1)/(γ−1)
M12

γ+1
p0 = p1 1/(γ−1) .
2
γ − 12 (γ − 1)/M12
Ondas de Choque
La presión en el vértice del objeto romo,

(γ+1)/(γ−1)
M12

γ+1
p0 = p1 1/(γ−1) ,
2
γ − 12 (γ − 1)/M12
es menor (para M1 > 1) que la que se obtendría por pura

retardación adiabática:
γ + 1 (γ+1)/(γ−1) M12

p0 = p1 .
2 γ 1/(γ−1)
Ondas de Choque
Características
Ecuaciones fundamentales
∂ρ
+ ∇ · j = 0; (Ecuación de continuidad.)
∂t
∂ρ ∂v ∂ρ
+ρ +v = 0;
∂t ∂x ∂x

∂p ∂p ∂v ∂p
+v + ρc2 = 0; c = 2
(7)
∂t ∂x ∂x ∂ρ S
∂v ∂v 1 ∂p
+v + = 0; (Ecuación de Euler.) (8)
∂t ∂x ρ ∂x
Ondas de Choque
Características
Características
Ondas de Choque
Características
Propagación de los valores iniciales
Ondas de Choque
Problemas resueltos
Pistón descomprimiendo un cilindro
Como el gas vecino al cilindro se mueve junto con el, las

características parten formando un ángulo fijo con su
trayectoria t − x.
Ondas de Choque
Problemas resueltos
Perfil v contra x
Ondas de Choque
Problemas resueltos
Pistón comprimiendo un cilindro
Transcurrido un tiempo breve (en verde) se tiene una solución

univaluada
Ondas de Choque
Problemas resueltos
Perfil v contra x (tiempos breves)
Ondas de Choque
Problemas resueltos
Cáustica
La familia de rayos reflejados por un espejo esférico tiene una

envolvente. De acuerdo con la óptica geométrica, la intensidad
de la luz es infinita en la envolvente.
Ondas de Choque
Problemas resueltos
Perfil v contra x (tiempos largos)
Para tiempos largos (en rojo) se obtiene una solución

multivaluada.
Ondas de Choque
Problemas resueltos
Flujo a lo largo de un terreno convexo
El flujo a la izquierda de la característica OA es uniforme. A su

derecha, es cada vez más rarificado. No se forma ondas de
choque.
Ondas de Choque
Problemas resueltos
Flujo a lo largo de un terreno cóncavo
A lo largo de OB se forma un frente de choque. Por debajo de

O se tiene un flujo regular con compresión creciente.
Ondas de Choque

Hugoniot

Cargado por

Copyright:

Formatos disponibles

Hugoniot

Cargado por

Información del documento

Descripción original:

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Hugoniot

Cargado por

Copyright:

Formatos disponibles

Introducción Propiedades de la onda de choque Aplicaciones Propagación de perturbaciones en un gas en movimiento

Introducción a las Ondas de Choque

Centro de Física Experimental, Facultad de Ciencias, Universidad de Chile

Curso de Pre- y Postgrado ONDAS DE CHOQUE, 2008

La naturaleza de las ondas de choque:

La naturaleza de las ondas de choque:

Figura: Ejemplo de una onda de choque

La naturaleza de las ondas de choque:

Figura: Hugoniot y compañeros (1870)

La naturaleza de las ondas de choque:

Figura: Formación de una onda de choque

La naturaleza de las ondas de choque:

Figura: Pasaje de una onda sonora y de una onda de choque.

La naturaleza de las ondas de choque:

Al paso de una onda sonora (1), el sistema oscila a lo largo de

Continuidad en las superficies de discontinuidad:

¿ Qué cambia de un modo continuo a través de la superficie de

Continuidad en las superficies de discontinuidad:

Onda de marea propagándose a lo largo de un canal rectangular:

v es la velocidad del agua detrás del frente de discontinuidad. u

Continuidad en las superficies de discontinuidad:

Masa adicional que adquiere rapidez v = ρh2 u∆t − ρh2 v∆t

(ρh2 u∆t − ρh2 v∆t)v = ( 21 ρgh22 − 12 ρgh21 )∆t

Continuidad de la corriente de masa:

j = corriente de masa = densidad de momentum.

Continuidad de la corriente de momentum:

p1 + ρ1 v12 = p2 + ρ2 v22 . (2)

Continuidad de la corriente de energía:

w1 + 21 v12 = w2 + 21 v22 . (3)

Figura: p1 , V1 estados iniciales. p2 , V2 estados finales.

Ondas de choque en gases politrópicos

Ondas de choque en gases politrópicos

En un gas politrópico la adiabática de choque (o de Hugoniot)

donde M1 = v1 /c1 es el número de Mach.

Ondas de choque en gases politrópicos

Adiabática de Hugoniot en un gas politrópico

Ondas de choque en gases politrópicos

Relación entre las temperaturas a uno u otro lado de la

Ondas de choque en gases politrópicos

Relación entre las presiones a uno u otro lado de la

Ondas de choque en gases politrópicos

Relación entre los números de Mach (o las velocidades) a uno

γM12 M22 − 12 (γ − 1)(M12 + M22 ) = 1.

Ondas de choque fuertes

Ondas de choque fuertes

En cambio, T2 puede ser arbitrariamente mayor que T1 .

Atenuación de las ondas de choque

¿Dónde se produce la discontinuidad?

Figura: Formación de una onda de choque.

Atenuación de las ondas de choque

Figura: Las áreas achuradas son iguales.

Atenuación de las ondas de choque

En primera aproximación, hay una relación lineal entre la

para un gas politrópico.

Atenuación de las ondas de choque

Área ABC = Área A0 DE

Atenuación de las ondas de choque

Distancia entre B 0 y C 0 = α(∆v)1 t.

Atenuación de las ondas de choque