UNIVERSIDAD NACIONAL MAYOR DE SAN

MARCOS
(Universidad del Perú, Decana de América)

Facultad de Ciencias Físicas

E.A.P. INGENIERÍA MECÁNICA DE FLUIDOS

Curso: Aerodinámica I

Sección: II

Semestre: 2017-1

Profesor: Dr. Miguel Ernesto Arámbulo Manrique

Trabajo: Onda de choque en toberas convergente -

divergente

Alumno:
 ALVARADO HUAMÁN, JULIO 13130034

2017
 Ondas de Choque en Toberas
Concepto:
En la mecánica de fluidos, una onda de choque es una onda de presión abrupta producida por un objeto que
viaja más rápido que la velocidad del sonido en dicho medio, que a través de diversos fenómenos produce
diferencias de presión extremas y aumento de la temperatura (si bien la temperatura de remanso
permanece constante de acuerdo con los modelos más simplificados). La onda de presión se desplaza como
una onda de frente por el medio.

Una de sus características es que el aumento de presión en el medio se percibe como explosiones.

También se aplica el término para designar a cualquier tipo de propagación ondulatoria, y que transporta,
por tanto energía a través de un medio continuo o el vacío, de tal manera que su frente de onda comporta
un cambio abrupto de las propiedades del medio.

Ejemplos
 Explosiones, como por ejemplo bombas cuyas ondas son las responsables de mover objetos y
destruirlos. Para esas ondas de detonación existen modelos matemáticos empíricos y teóricos
exactos.
 Los aviones supersónicos provocan ondas de choque al volar por encima de régimen transónico (M
> 0,8) pues aparecen zonas donde el aire supera la velocidad del sonido localmente, por ejemplo
sobre el perfil del ala, aunque el propio avión no viaje a M > 1.
 Meteoritos que entran en la atmósfera producen ondas de choque. El aumento de temperatura
producido por la onda de choque es la responsable de que se vean los meteoros.
 En los alrededores del canal del relámpago hay un aire muy caliente que, con ondas de choque,
produce el trueno en tormentas. Es decir que es como una explosión a lo largo del camino que
recorre el relámpago. Debido a las fluctuaciones irregulares que influyen el camino de las ondas, no
solo se oye un golpe sino una serie de más o menos golpes fuertes en una distancia lejana.
 En los propulsores de los cohetes pueden aparecer ondas de choque si han sido mal diseñados.
Esas ondas pueden causar la destrucción del cohete, por lo que deben ser amortiguadas.
Características de una tobera convergente-divergente cuando se presentan las
ondas de choque.

Figura 1

Punto d: PE = PB y se ha mantenido el flujo isentrópico.

Punto f: La presión PE en el plano de salida, no se ve afectada por el aumento de presión desde el
punto d hasta este punto. El aumento de PE a PB ocurre fuera de la tobera.
Punto g: Su valor de contrapresión es justo el necesario para originar un choque normal permanente
en el plano de salida de la tobera. PE en el plano de salida (flujo abajo del choque), es igual a la
contrapresión PB y M < 1 a la salida de la tobera.
Punto h: Al subir la contrapresión del punto g al punto h, el choque normal se mueve dentro de la
tobera.

ECUACIONES A UTILIZAR:
Para flujo isentrópico de un gas ideal:

ṁ = 𝜌𝐴𝑉 (1)

𝑃 = 𝜌𝑅𝑇 (2)

𝐶 = √𝐾𝑅𝑇 (3)

𝑉
𝑀= (4)
𝐶

Donde:
ṁ : flujo másico
𝜌 : densidad
A : área
V : velocidad
R : constante del gas
P : presión
T : temperatura
C : velocidad del sonido
K : razón de los calores específicos
M : número de mach.
Expresiones para las propiedades de estancamiento isentrópico local de un gas ideal.

𝐾
𝑃0 𝐾−1 𝐾−1
= [1 + 𝑀2 ] (5)
𝑃 2

𝑇0 𝐾−1
=1+ 𝑀2 (6)
𝑇 2

1
𝜌0 𝐾−1 𝐾−1
= [1 + 𝑀2 ] (7)
𝜌 2

Donde:
𝑃0 : Presión de estancamiento.
𝑇0 : Temperatura de estancamiento.
𝜌0 : Densidad de estancamiento.

Para la condición crítica M = 1

𝐾+1
𝐾−1
2 2(𝐾−1)
𝐴 1 1+ 2 𝑀
= [ 𝐾−1 ] (8)
𝐴∗ 𝑀 1+
2

Ecuaciones para onda de choque normal en el flujo de un gas ideal que fluye a través de una tobera:

Figura 2

𝑇0𝑥 = 𝑇0𝑦 (9)

Nota: los subíndices x y y se usan para indicar las condiciones flujo arriba y flujo abajo de la onda de choque,
respectivamente.

𝐾−1
𝑇𝑦 1+ 2 𝑀2 𝑥
= 𝐾−1 2 (10)
𝑇𝑥 1+ 𝑀 𝑦
2

𝑇𝑦 𝑃 2 𝑀𝑦 2
= ( 𝑦) ( ) (11)
𝑇𝑥 𝑃𝑥 𝑀𝑥

𝐾−1
𝑃𝑦 𝑀𝑥 √1+ 2 𝑀2 𝑥
= (12)
𝑃𝑥 𝐾−1 2
𝑀𝑦 √1+ 𝑀 𝑦
2
 𝑃𝑦 1+𝐾𝑀2 𝑥
= (13)
𝑃𝑥 1+𝐾𝑀2 𝑦

2
𝑀2 𝑥 +𝐾−1
𝑀2 𝑦 = 2𝐾 2 (14)
𝑀 𝑥 −1
𝐾−1

𝐾
𝑃0𝑦 𝑃 𝑇 𝐾−1
= ( 𝑦) ( 𝑦 ) (15)
𝑃0𝑥 𝑃𝑥 𝑇𝑥

PROBLEMA

a)Una tobera convergente - divergente tiene una relación de área de salida / área de la garganta de 2.
Entra aire a la garganta con una presión de estancamiento de 360K El área de la garganta es 500
mm2:
Determine:
1. La velocidad y la presión del aire a la salida, si la parte divergente actúa como tobera y el proceso es
isentrópico (Pd).
2. La velocidad y la presión del aire a la salida, si la parte divergente actúa como difusor y el proceso
es isoentrópico (Pc).
3. La velocidad y la presión del aire a la salida, si se produce una onda de choque justo en la salida
(Pg).

DATOS
AE = A* = 2:0
A* = 500mm2 = 0.0005m2
P0 = 1000kPa
T0 = 360K

CALCULOS

1.-
PE =? VE =?
M > 1 (flujo supersónico) Representa el punto d en Figura 1
Flujo isoentrópico.

Despejando el valor de M de la ecuación (8) con K = 1:4 y A/A* = 2 se obtiene:

ME = 2.197 (16)
ME = 0.3079 (17)

Como en este caso la sección divergente está actuando como tobera, M > 1 y se selecciona el valor de
ME = 2.197.
Resolviendo para el valor de P 0/PE de la ecuación (5) con K = 1.4 y ME = 2.197 se obtiene:

𝑃0
= 10.64 (18)
𝑃𝐸

𝑃0 100 𝐾𝑃𝑎
𝑃𝐸 = =  𝑃𝐸 = 93.98 𝐾𝑃𝑎 (19)
10.64 10.64

Resolviendo para el valor de T0/TE de la ecuación (6) con K = 1.4 y ME = 2.197 se obtiene:
𝑇0
= 1.9654 (20)
𝑇𝐸

𝑇0 360 𝐾
𝑇𝐸 = =  𝑇𝐸 = 183.17 𝐾 (21)
1.9654 1.9654

Resolviendo para el valor de CE de la ecuación (3) se obtiene:

𝐽
𝐶𝐸 = √(1.4)(287 )(183.17 𝐾) = 271.29 𝑚/𝑠 (22)
𝐾𝑔.𝑠

Despejando para el valor de VE de la ecuación (4) se obtiene:

𝑚
𝑉𝐸 = 𝑀𝐸 𝐶𝐸 → 𝑉𝐸 = (2.197) (271.29 ) → 𝑉𝐸 = 596.02 𝑚/𝑠 (23)
𝑠

2.-
PE =? VE =?
M < 1 (flujo subsónico) Representa el punto c en Figura 1
Flujo isoentrópico.

Despejando el valor de M de la ecuación (8) con K = 1.4 y A/A* = 2 se obtiene:

ME = 2.197 (24)
ME = 0.3079 (25)

Como en este caso la sección divergente está actuando como difusor, M < 1 y se selecciona el valor de
ME = 0.3079 Resolviendo para el valor de PE/P0 de la ecuación (5) se obtiene:

𝑃𝐸
= 0.9360 (26)
𝑃0

𝑃𝐸 = 0.9360 = 0.9360 (1000 𝐾𝑃𝑎) → 𝑃𝐸 = 936 𝐾𝑃𝑎 (27)

Resolviendo para el valor de TE/T0 de la ecuación (6) se obtiene:

𝑇𝐸
= 0.9813 (28)
𝑇0

𝑇𝐸 = 0.9813𝑇0 = 0.9813(360𝐾) → 𝑇𝐸 = 353.268 𝐾 (29)

Resolviendo para el valor de CE de la ecuación (3) se obtiene:

𝐽
𝐶𝐸 = √(1.4)(287 )(353.268 𝐾) = 376.75 𝑚/𝑠 (30)
𝐾𝑔.𝑠

Despejando para el valor de VE de la ecuación (4) se obtiene:

𝑚
𝑉𝐸 = 𝑀𝐸 𝐶𝐸 → 𝑉𝐸 = (0.3079) (376.75 ) → 𝑉𝐸 = 116.00 𝑚/𝑠 (31)
𝑠
3.-
PE =? VE =?
Si se produce onda de choque justo en la salida Representa el punto g en Figura 1

Las condiciones en x, que es un instante antes de que ocurra la onda de choque, se conocen de la parte 1 y
estos resultados son:

Mx = 2.197 (32)
Px = 93.98kPa (33)
Tx = 183.17K (34)
P0x = 93.98kPa (35)

Resolviendo para el valor de My de la ecuación (14) se obtiene:

2
2.1972 +
𝑀𝑦 = √ 1.4 − 1 → 𝑀𝑦 = 0.54746
2(1.4)
∗ 2.1972 − 1
1.4 − 1

Resolviendo para el valor de Py/Px de la ecuación (13) se obtiene:

𝑃𝑌 1 + (1.4)(2.197)2 𝑃𝑌
= → = 5.46459
𝑃𝑋 1 + (1.4)(0.54746)2 𝑃𝑋

Resolviendo para el valor de Ty/Tx de la ecuación (11) se obtiene:

𝑇𝑌 0.54746 2 𝑇𝑌
= 5.464592 ( ) → = 1.854
𝑇𝑋 2.197 𝑇𝑋
Resolviendo para el valor de P 0y/P0x de la ecuación (15) se obtiene:

𝑃0𝑌 1.4 𝑃0𝑌

= 5.46459(1.854)1−1.4 → = 0.63
𝑃0𝑋 𝑃0𝑋

Entonces:

𝑃𝑦 = 5.46459𝑃𝑥 = 5.46459 (93.98𝑘𝑃𝑎) → 𝑃𝑦 = 𝑃𝐸 = 513.562𝑘𝑃𝑎

𝑃0𝑌 = 0.630𝑃0𝑋 = 0: 630 (1000𝑘𝑃𝑎) → 𝑃0𝑌 = 630.00𝑘𝑃𝑎

𝑇𝑦 = 1.854𝑇𝑥 = 1.854 (183.17𝐾) → 𝑇𝑦 = 𝑇𝐸 = 339.60𝐾

Resolviendo para el valor de Cy de la ecuación (3) se obtiene:

𝐽
𝐶𝑦 = √(1.4) (287 ) (339.60𝐾) = 369.39 𝑚/𝑠
𝐾𝑔. 𝑠

Despejando para el valor de Vy de la ecuación (4) se obtiene:

𝑚
𝑉𝐸 = 𝑉𝑌 = 𝑀𝑌 𝐶𝑌 → 𝑉𝐸 = (0.54746) (369.39 ) → 𝑉𝐸 = 202.23 𝑚/𝑠
𝑠
 b) En un turborreactor, los gases de combustión se descargan en una tobera convergente-divergente;
si la sección de la garganta tiene área crítica, el flujo que pasa por la garganta es sónico y puede
evolucionar a supersónico por la geometría divergente. El chorro supersónico se descarga a presión
atmosférica, con lo que la presión de salida de la tobera será la atmosférica; en función del valor de
la presión de salida con respecto a la de estancamiento (p0/ps) se pueden tener ondas de choque
en la zona divergente lo que impide que el flujo de salida sea supersónico. Para un valor concreto
de la relación de presiones, no se produce ninguna onda de choque, y se tiene una tobera adaptada
que genera un chorro de salida supersónico. A partir de los datos:

DETERMINE:
1. Condiciones de remanso.
2. Condiciones críticas.

DATOS: Sección de entrada de la tobera: AE = 0,6 m2; MaE = 0,4; TE = 2000 K. Sección en donde se produce la
onda de choque: A1=A2=(A* +AS)/2. Ambiente: presión atmosférica a una altura de vuelo de 10000 m,
patm = 26,42 kPa

1. CONDICIONES DE REMANSO: conociendo la temperatura y el número de Mach, se tiene la

temperatura de remanso

𝑇0
= 1 + 0.2𝑀𝑎𝐸2
𝑇𝐸

𝑇0 = 𝑇𝐸 (1 + 0.2𝑀𝑎𝐸2 ) = 2000(1 + 0.2 ∗ 0.42 ) = 2064 𝐾

La presión de remanso se puede obtener a partir del caudal másico, y de las condiciones de entrada
(E)

𝜌0 𝑀𝑎
ṁ=𝐴 √1.4
√𝑅𝑇0 (1 + 0.2𝑀𝑎2 )3

Despejando 𝜌0 :

39 √287 ∗ 2064 (1 + 0.2 ∗ 0.42 )3

𝜌0 = = 116.18 𝑘𝑃𝑎
0.6 √1.4 0.4

2. CONDICIONES CRÍTICAS: el área crítica se obtiene a partir de las condiciones de la sección de

entrada:

𝐴 1 (1 + 0.2𝑀𝑎2 )³
=
𝐴∗ 𝑀 1.728

Despejando 𝐴∗ :

1.728
𝐴∗ = 0.6 ∗ 0.4 = 0.3773 𝑚²
(1 + 0.2 ∗ 0.42 )3

La temperatura y la presión en condiciones críticas, se obtienen a partir de las de remanso:

𝑇0 2064
𝑇∗ = = = 1720 𝐾
1.2 1.2
∗ 𝑃0 116.18
𝑃 = = = 81.38 𝑘𝑃𝑎
1.23.5 1.23.5
Como la presión crítica es mayor que la de salida (atmosférica = 26,42 kPa), a partir de la garganta (en una
tobera adaptada de área igual a la crítica) el flujo será supersónico hasta la descarga en toberas adaptadas;
en caso contrario entre la garganta y la salida aparecen ondas de choque, que llevan el flujo de salida a
subsónico.