Flujo Compresible Capítulo VI: Ondas de choque

CAPITULO VI ONDAS DE CHOQUE

6.1 PROPAGACION DE LAS PERTURBACIONES FINITAS

Anteriormente se mostró que las pequeñas perturbaciones elásticas se propagan en un medio

dado (líquido o gas) con una velocidad igual a la velocidad local del sonido. En mecánica de fluidos,
así mismo; en otros capítulos de Física, además de las pequeñas perturbaciones, tiene a lugar otros
tipos de perturbaciones que no pueden ser consideradas como pequeñas, las cuales generalmente
son llamadas perturbaciones finitas, por ejemplo; fuertes explosiones, ondas de detonación, ondas de
choque formadas en los bordes de ataque de un proyectil o avión que se mueve con velocidad
supersónica, etc.

En general, cuando ocurren pequeñas perturbaciones, el valor del incremento de cualquiera de

los parámetros es pequeño en comparación con su valor hasta antes de ocurrir la perturbación. Por
ejemplo, una pequeña perturbación de presión en un fluido cualquiera puede ser caracterizada por
una magnitud adimensional del tipo:

p − p0 ∆p
∆p = =
1 (6.1)
p0 p0

donde p0 es la presión del fluido hasta antes de aparecer la perturbación, p es el valor máximo de
la presión después de ocurrida la pequeña perturbación. En las perturbaciones finitas, el valor de ∆p
puede ser del orden de la unidad y significativamente mayor.

La diferencia fundamental entre estas dos formas de perturbaciones consiste en el

comportamiento de los parámetros del flujo (donde ocurre la perturbación) y que determinan el
movimiento del flujo. En las pequeñas perturbaciones, todos los parámetros del flujo son funciones
continuas de las coordenadas y del tiempo; pero, en las perturbaciones finitas los parámetros del flujo
(velocidad, presión, temperatura y densidad) sufren un salto finito; siendo esto lo que diferencia a las
pequeñas perturbaciones de las perturbaciones finitas.

Las perturbaciones que ocurren en flujos compresibles, así mismo; la propagación del sonido,
pueden ser perturbaciones pequeñas o finitas dependiendo de las condiciones de cómo ocurren. Es
sabido que en condiciones normales las perturbaciones acústicas son pequeñas perturbaciones y se
propagan con la velocidad del sonido, pero en fuertes explosiones estas serán perturbaciones finitas
y su velocidad de propagación puede superar significativamente a la velocidad del sonido.

Tanto el movimiento de las pequeñas perturbaciones como el movimiento de las

perturbaciones finitas, son descritas por ecuaciones completamente diferentes. La primera puede
describirse a través de una ecuación diferencial lineal en derivadas parciales, llamada ecuación de
onda y que generalmente tiene la forma de:

∂2 f ∂2 f
2
= u2 2
∂t ∂x
En tanto que el movimiento en perturbaciones finitas, se describe con una ecuación diferencial
no lineal y que tiene una solución compleja, dependiendo de las condiciones de ocurrencia.

6.2 GENERACION DE LAS ONDAS FINITAS Y ONDAS DE CHOQUE

La formación de una perturbación finita en un gas puede analizarse a través de un cuadro

cualitativo, considerando la compresión del gas; observándose la Fig. 6.1 siguiente

Ing. César A. Quispe Gonzáles, M Sc. Página 104

 Flujo Compresible Capítulo VI: Ondas de choque

Figura 6.1 - Esquema de la generación de una onda finita de compresión.

Sea que inicialmente en un conducto finito, dentro del cual se tiene un gas estancado con
parámetros p0 , T0 y ρ 0 . Para un tiempo t1 cualquiera, un pistón que se encuentra dentro del
conducto, comienza a moverse instantáneamente, comprimiendo la capa de gas adyacente a su
superficie y generando una perturbación que viaja hacia la derecha con la velocidad del sonido
a1 = kRT1 . La onda se propaga en un plano y delante de ella, los parámetros quedan invariables,
ya que estos no sufren la influencia de la onda, pero detrás de ella, los parámetros sufren una
variación por efecto de la generación de la perturbación y la presión sufre un incremento siendo igual
a p + dp .

Para el instante de tiempo t 2 , se origina una nueva perturbación que también se moverá con la
velocidad local del sonido a2 = kRT2 , pero como T2 > T1 debido a la compresión del gas, la onda
a2 se propagará con mayor velocidad que la onda a1 , alcanzándola en algún instante, de tal manera
que los planos de las ondas 2 y 1 se superponen creando un nuevo frente que divide el flujo en no
perturbado delante de éste frente y perturbado atrás de él.

Si se supone que el pistón es acelerado del reposo, variando su velocidad en ∆V . La

velocidad final puede ser idealizada como una composición de un cierto número de incrementos de
velocidad que es alcanzada debido a los incrementos infinitesimales de velocidad, cada uno de ellos
con magnitud dV como se observa en la Fig. 6.2. Así mismo, las ondas más próximas al pistón
tienden a alcanzar a aquellas que se encuentran más distantes dentro del tubo, sumándose una tras
otra, generando un solo frente en el cual la diferencia de de presión a ambos lados de este frente va
aumentando a medida que se suman más ondas, generando una fuerte onda llamada onda de
choque.

Figura 6.2 - Incrementos de velocidad e incrementos de presión en la onda de compresión

En base a todo lo analizado, se puede decir que la onda de choque es una onda poderosa que
se forma como consecuencia de la superposición de otras perturbaciones, de espesura finita y muy

Ing. César A. Quispe Gonzáles, M Sc. Página 105

 Flujo Compresible Capítulo VI: Ondas de choque

pequeña, comparable con el libre camino medio de las moléculas λ y que a condiciones
atmosféricas se puede medir en Amstrong ( 
A=10 -10
)
m .

Tras el pistón, en el dominio de caída de presión, el cuadro será inverso, ya que las ondas se
propagaran en un medio rarificado (con la presión disminuyendo), la temperatura disminuye al igual
que la velocidad del sonido, lo que hace que las ondas iniciales se propaguen a mayor velocidad que
las ondas posteriores.

Si el sistema de coordenada está unido al conducto; entonces el movimiento de la onda de

choque es no estacionario y todos los parámetros del flujo parámetros del flujo dependerán del
tiempo en cualquier punto. Para transformar el problema en estacionario, es necesario introducir un
sistema de coordenadas ligado al plano de la onda de choque, entonces la onda puede considerarse
estacionaria y los parámetros adelante y atrás de la onda serán valores constantes que varían en la
onda de choque.

En esta estrecha región, donde se desarrolla la onda de choque; se producen drásticas

variaciones de los parámetros del fluido, las gradientes dentro de la onda de choque son elevadas y
por tanto, los efectos viscosos (efectos de disipación) tienen una determinante influencia en el
proceso, resultando que dentro de la onda de choque, el flujo sea no isoentrópico (Fig. 6.3)

Para una onda estacionaria, se notará que tanto la presión, la temperatura y la densidad sufren
un fuerte incremento al pasar por la onda de choque, mientras que la velocidad sufrirá una
disminución drástica de su valor. El flujo dentro de la onda será no isoentrópico, pero fuera de ella
sigue siendo isoentrópico y adiabático.

Figura 6.3 - Variación de los parámetros del flujo en una onda de choque

La formación de la onda de choque se da en flujos supersónicos, cuando en su movimiento

encuentra el borde de ataque de cualquier cuerpo que generan perturbaciones en el flujo (Fig. 6.4).

Figura 6.4 - Formación de la onda de choque

Todo flujo está conformado por una gran cantidad de moléculas, las cuales inicialmente
impactan sobre el borde de ataque del cuerpo que está delante del flujo. Debido a este impacto hay
una variación de la energía molecular y de cantidad de movimiento.

Ing. César A. Quispe Gonzáles, M Sc. Página 106

 Flujo Compresible Capítulo VI: Ondas de choque

Para el caso de un flujo subsónico V∞ < a∞ , en este caso, las perturbaciones provocadas por
la presencia de un cuerpo, se propagan en ambas direcciones, hacia adelante con velocidad V∞ + a∞
y hacia atrás del flujo, a∞ − V∞ llevando información sobre la presencia del cuerpo y haciendo que las
líneas de corriente se ajusten a la presencia de él.

Para el caso de un flujo supersónico V∞ > a∞ , las perturbaciones se propagan hacia adelante
con velocidad V∞ + a∞ , pero no pueden volver hacia atrás, ya que las perturbaciones son arrastradas
por el flujo, debido a su mayor velocidad; no habiendo información alguna sobre la presencia del
cuerpo. Estas ondas de perturbación, al no poder propagarse hacia atrás, generan un proceso de
coalescencia (aglutinamiento) a una corta distancia del borde de ataque del cuerpo, generando un
frente de onda llamada onda de choque.

Las ondas de choque pueden clasificarse como ondas de choque normales, cuando el plano de
la onda es perpendicular a la línea de corriente; onda de choque oblicua, cuando el plano de la onda
forma algún ángulo β con las líneas de corriente y onda de choque curvilínea, cuando la onda se
desarrolla sobre una superficie curvilínea.

6.3 CONSIDERACIONES PARA EL ESTUDIO DE LAS ONDAS DE CHOQUE NORMALES

Para el estudio de las ondas de choque normales, el flujo puede considerarse como un flujo
unidimensional, donde los parámetros del flujo varían solo en una dirección, que corresponde a la
dirección del flujo. Además, se debe considerar las siguientes hipótesis:
i) La onda de choque es perpendicular a las líneas de corriente.
ii) Flujo a régimen permanente, la derivada de cualquier parámetro P respecto al tiempo es
nula ( ∂P dt = 0 ) .
iii) Las líneas de corriente del flujo son paralelas, generando áreas constantes entre ellas.
iv) Las fuerzas actuantes sobre el cuerpo son despreciadas.
v) El proceso de todo el sistema es adiabático (sin transferencia de calor dq = 0 )
vi) No se realiza ningún trabajo externo.

6.4 RELACION FUNDAMENTAL DE LA ONDA DE CHOQUE NORMAL

Para encontrar las relaciones que existen entre los parámetros del flujo, entre las secciones
delante de ella (corriente abajo) respecto a los parámetros atrás de la onda (corriente arriba) se
utilizan las ecuaciones de conservación de mas, de cantidad de movimiento y de energía. El esquema
de cálculo de la onda de choque normal y estacionaria se muestra en la Fig. 6.5

Figura 6.5 - Esquema de cálculo de la onda de choque normal

Aplicando la ecuación de continuidad, considerando que por ser la onda de espesura muy
ínfima, la sección de paso permanece constante se tiene:

ρ1 V1 = ρ 2 V2 (6.2)

Aplicando la ecuación de cantidad de movimiento, considerando que no existen fuerzas

volumétricas y superficiales, de tiene:

Ing. César A. Quispe Gonzáles, M Sc. Página 107

 Flujo Compresible Capítulo VI: Ondas de choque

p1 + ρ1 V12 = p2 + ρ2 V22 (6.3)

La ecuación de energía para la onda de choque normal, puede expresarse mediante la

relación:

V12 k p1 V22 k p2
+ = + (6.4)
2 k − 1 ρ1 2 k − 1 ρ2

Para hallar la relación entre las velocidades delante y detrás de la onda, se divide la Ec. (6.3)
entre la Ec. (6.2) obteniéndose:

p1 + ρ1 V12 p2 + ρ 2 V22 p1 p
= ⇒ + V1 = 2 + V2 (6.4)
ρ1 V1 ρ 2 V2 ρ1 V1 ρ 2 V2

La ecuación de energía, antes y después de la onda, también puede escribirse bajo la forma:

V12 k p1 a2 k + 1 *2 p1 k + 1 *2 ( k − 1) 2
+ = 0 = a ⇒ = a − V1
2 k − 1 ρ1 k − 1 2 ( k − 1) ρ1 2k 2k
V22 k p2 a2 k + 1 *2 p2 k + 1 *2 ( k − 1) 2
+ = 0 = a ⇒ = a − V2
2 k − 1 ρ 2 k − 1 2 ( k − 1) ρ2 2k 2k

Reemplazando estos resultados en la Ec. (6.4) y realizando algunas transformaciones, se

tiene:

1  k + 1 *2 ( k − 1) 2  1  k + 1 *2 ( k − 1) 2 
 a − V1  + V1 =  a − V2  + V2
V1  2k 2k  V2  2k 2k 
k + 1 a*2  ( k − 1)  k + 1 a*2  ( k − 1) 
− V1  − 1 = − V2  − 1
2k V1  2k  2k V2  2k 
k + 1 *2  1 1   k +1 k + 1 *2  V2 − V1   k +1 
a  −  = (V2 − V1 )   ⇒ a 
 V1V2
 = (V2 − V1 ) 
 
2k  V1 V2   2k  2k    2k 

lo que permite encontrar una relación fundamental y muy importante para la onda de choque normal:

V1 V2 a*2
V1 V2 = a*2 ⇒ = ⇒ λ1 λ2 = 1 (6.5)
a* a* a*2
Quiere decir que el producto de la velocidad del flujo delante de ella por el valor de la velocidad
detrás de ella es igual al cuadrado de la velocidad crítica del flujo; ó también, el producto de los
coeficientes de velocidad antes y después de la onda es igual a la unidad.

Como en la onda de choque normal V1 > V2 , entonces se tiene que λ1 > λ2 y para que se
cumpla la condición de la Ec. (6.5) es necesario que λ1 > 1 y λ2 < 1 . Estas últimas desigualdades
denotan que cuando el flujo atraviesa una onda de choque normal, pasa de un régimen supersónico a
un régimen subsónico.

Ing. César A. Quispe Gonzáles, M Sc. Página 108

 Flujo Compresible Capítulo VI: Ondas de choque

6.5 VARIACION DE LOS PARAMETROS DEL FLUJO EN LA ONDA DE CHOQUE

La variación de los parámetros en una onda de choque normal puede ser calculada de manera
cuantitativa si es que se asume que los parámetros del flujo incidente (antes de la onda de choque)
son conocidos.

Para calcular la variación de la velocidad ∆V = V1 − V2 , se utiliza la Ec. (6.5) de donde se tiene

que V2 = a*2 V1 , obteniéndose:

a*2  a*2   1 
∆V = V1 − V2 = V1 − = V1 1 − 2  ⇒ ∆V = V1 − V2 = V1 1 − 2  (6.6)
V1  V   λ 
 1   1 

La variación de presión en la onda de choque normal se halla a partir de la Ec. (6.3) de donde:

∆p = p2 − p1 = ρ1 V12 − ρ 2 V22 = ρ1 V1 (V1 − V2 ) = ρ1 V1 ∆V

Sustituyendo en esta ecuación, la expresión hallada en la Ec. (5.6) se tiene:

 1 
∆p = p2 − p1 = ρ1 V12 1 − 2  (6.7)
 λ1 

Para encontrar la diferencia de la densidad, se utiliza la Ec. (6.2) de donde ρ 2 = ρ1V1 V2 , y de

*2
la Ec. (6.5) se obtiene V2 = a V1 , entonces:

V1 V2  V2 
∆ρ = ρ 2 − ρ1 = ρ1 − ρ1 = ρ1  1*2 − 1 = ρ1  1*2 − 1 ⇒ ∆ρ = ρ 2 − ρ1 = ρ1 ( λ12 − 1) (6.8)
V2 a  a 

El salto de temperatura se halla a partir de la ecuación de energía, escrita bajo la siguiente

forma:

V12 k V2 k k −1 2
+
2 k −1
RT1 = 2 +
2 k −1
RT2 ⇒ ∆T = T2 − T1 =
2k R
(V1 − V22 )

El valor de V2 = a*2 V1 se halla de la Ec. (6.5), por lo que al sustituir en la ecuación anterior se
obtiene:

k − 1  2 a*4  k − 1 2  a*4  k −1 2  1 
∆T =  V1 − 2  = V1 1 − 4  ⇒ ∆T = T2 − T1 = V1 1 − 4  (6.9)
2k R  V1  2 k R  V1  2 k R  λ1 

Las expresiones anteriormente halladas, expresan la variación de los parámetros del flujo al
pasara a través de una onda de choque, en función del coeficiente de velocidad λ1 . Para expresar
estas variaciones en función del número de Mach M1 , se debe tener en cuenta:

1 k −1 2 1
2
= +
λ 1 k + 1 k + 1 M12

Por lo que, al sustituir el valor de λ por su correspondiente valor de M, se obtiene:

Ing. César A. Quispe Gonzáles, M Sc. Página 109

 Flujo Compresible Capítulo VI: Ondas de choque

 k −1 2 1  2  1 
∆V = V1 1 − + 2 
⇒ ∆V = V1 1 − 2  (6.10)
 k + 1 k + 1 M1  k + 1  M1 
2  1 
∆p = ρ1V12 1 − 2  (6.11)
k +1  M1 
M12 − 1
∆ρ = ρ1 (6.12)
k −1 2
1+ M1
2
k −1   k −1 2 1  
2

∆T = T1M1 1 − 
2
+ 2 
 (6.13)
2   k + 1 k + 1 M1  

En todas las fórmulas deducidas anteriormente, puede observarse que la variación de

parámetros solo es posible para el caso en que λ1 > 1 o que M 1 > 1 . Si se considera que
λ1 = M1 = 1 , entonces; las variaciones de parámetros serían nulas y se tendría una función continua
y para λ1 < 1 ó M 1 < 1 no se podría satisfacer la condición expresada por la Ec. (6.5) por lo que las
ondas de choque normales solo son posibles en flujos supersónicos y no es posible en flujos sónicos
o subsónicos.

Los valores de de variación relativa de presión ∆p p1 y de variación relativa de densidad

∆ρ ρ1 pueden ser hallados con las fórmulas deducidas anteriormente. Para el caso de la variación
relativa de la densidad, se tiene de la Ec. (6.12):

 1   1 
M12 1 − 2  1 − 2 
∆ρ  M1  =  M1 
=
ρ1  1 k −1  1 k −1 
M12  2 +   2+ 
 M1 2   M1 2 

De esta ecuación, en el límite cuando M 1 → ∞ , la variación relativa de densidad será:

∆ρ 2 ρ k +1 ρ
= ó también 2 = y para el caso del aire, k = 1.4 ; obteniéndose 2 = 6
ρ1 k −1 ρ1 k − 1 ρ1

También se debe indicar que las relaciones anteriormente halladas no consideran los procesos
de ionización y disociación del gas a altas temperaturas. Cuando se tiene en cuenta procesos de
ionización y disociación térmica, la razón de densidades puede ser mayor que 6. Tal es el caso de las
comprobaciones cuando T1 = 2260 º C se tiene ρ 2 ρ1 = 6.04 ; cuando T1 = 14000 º C se tiene
ρ 2 ρ1 = 11.1 y cuando T1 = 50000 º C se tiene ρ 2 ρ1 = 8.97

6.6 RELACIONES PARA LA ONDA DE CHOQUE NORMAL

La ecuación de energía para una línea de corriente que pasa por una onda de choque, se
escribe como:

V12 k V2 k
+ RT1 = 2 + RT2 (6.14)
2 k −1 2 k −1
A su vez, la temperatura de estancamiento a ambos lados de la onda será:

Ing. César A. Quispe Gonzáles, M Sc. Página 110

 Flujo Compresible Capítulo VI: Ondas de choque

V12 k k 
+ RT1 = RT01 
2 k −1 k −1  V12 k V12 k
 ⇒ + RT1 = + RT2 ⇒ T01 = T02 (6.15)
V12
k k  2 k − 1 2 k −1
+ RT2 = RT02
2 k −1 k −1 

Se observa que de acuerdo a la Ec. (6.14), las partes izquierdas de las dos ecuaciones
anteriores son iguales, por lo que las partes derechas también deberán ser iguales, desprendiéndose
que en una onda de choque normal, la temperatura de estancamiento permanece constante.

Las razones de los parámetros del flujo en la onda de choque normal pueden hallarse en
función del coeficiente de velocidad λ . Si se escribe la temperatura de estancamiento como una
función del coeficiente de velocidad se tendrá:

T1 k −1 2 
= 1− λ1
T01 k + 1  T2 ( k + 1) − ( k − 1) λ2
2

 ⇒ =
T2 k −1 2  T1 ( k + 1) − ( k − 1) λ12
= 1− λ2
T02 k + 1 

El valor de λ1 se puede hallar de la Ec. (6.5) por lo que se puede escribir:

T2 ( k + 1) − ( k − 1) (1 λ1 )
2

= ⇒
T2
= 2
( k + 1) λ12 − ( k − 1) (6.16)
T1 ( k + 1) − ( k − 1) λ12 T1 λ1 ( k + 1) − ( k − 1) λ12 

La relación de densidades y velocidades se puede hallar a partir de la ecuación de continuidad

y con ayuda de la Ec. (6.5):

ρ 2 V1 V1 a* λ1 ρ 2 V1
ρ1 V1 = ρ 2 V2 ⇒ = = = ⇒ = = λ12 (6.17)
ρ1 V2 V2 a* λ2 ρ1 V2

La relación de presiones en la onda de choque normal, puede ser establecida de la ecuación

de estado:

p2 = ρ 2 RT2  p2 ρ 2 T2
 ⇒ =
p1 = ρ1 RT1  p1 ρ1 T1

Sustituyendo en esta expresión, las Ecs. (6.17) y (6.16) que expresan las relaciones de
densidad y de temperatura, permite obtener:

p2 2
= λ1 2
( k + 1) λ12 − ( k − 1)
⇒
p2 ( k + 1) λ1 − ( k − 1)
=
2

(6.18)
p1 λ1 ( k + 1) − ( k − 1) λ12  p1 ( k + 1) − ( k − 1) λ12

Para encontrar la relación de presiones de estancamiento en la onda de choque, se sigue el

siguiente razonamiento:

p02  p02   p2   p1 
=     
p01  p2  ISO  p1 O.C .  p01  ISO

Sustituyendo en esta ecuación, la relación hallada en la Ec. (6.18) se tendrá:

Ing. César A. Quispe Gonzáles, M Sc. Página 111

 Flujo Compresible Capítulo VI: Ondas de choque

k k
p02  k − 1 2 
−
k −1  ( k + 1) λ12 − ( k − 1)   k − 1 2  k −1
= 1 − λ2   2 
1− λ1 
p01  k + 1   ( k + 1) − ( k − 1) λ1   k + 1 

Como λ2 = 1 λ1 se tiene:

k k
−
p02  ( k + 1) λ1 − ( k − 1)   ( k + 1) λ12 − ( k − 1)   ( k + 1) − ( k − 1) λ12  k −1
2 k −1
=
p01  ( k + 1) λ12   2  
  ( k + 1) − ( k − 1) λ1   k +1 
k k
−
p02 ( k + 1) λ1 − ( k − 1)  k −1 ( k + 1) λ12 − ( k − 1)  ( k + 1) − ( k − 1) λ12  k −1
2

=  k 
p01 k 2 −  ( k + 1) − ( k − 1) λ12  k
−
( k + 1) λ1  k −1 
k − 1 ( )
k + 1 k −1
1
2 k  k + 1 − k − 1 λ 2  k −1
p02
= ( λ1 ) k −1 
( ) ( ) 1
 (6.19)
 ( k + 1) λ1 − ( k − 1) 
2
p01

Estas mismas razones, halladas en las Ecs. (6.16), (6.17), (6.18) y (6.19) también pueden
expresarse por medio del número de Mach M . Escribiendo las temperaturas de estancamiento en
función de M y recordando que ellas son iguales a ambos lados de la onda, se tiene:

T01 k −1 2 
= 1+ M1 
T1 2 T2 2 + ( k − 1) M1
2

 ⇒ = (6.20)
T02 k − 1 2  T1 2 + ( k − 1) M 22
= 1+ M2 
T2 2 

Asumiendo que las condiciones iniciales son conocidas, entonces es preferible expresar la Ec.
(6.20) en función del número M1 , para lo cual hace falta hallar una relación que enlace los números
M1 y M 2 . De la relación fundamental de la onda de choque normal expresada por la Ec. (5.5) se
puede obtener:

2
V V a* a * a* a0 a* a0  a*  a0 a0
V1V2 = a *2
⇒ 1 2 = ⇒ M1M 2 = = 
a1 a2 a1 a2 a0 a1 a0 a2  a0  a1 a2

Recordando que a = kRT , así como la relación de parámetros críticos y de estancamiento,

la expresión anterior puede reescribirse como:

2
 T*  T0 T0  2  k −1 2 k −1 2
M1M 2 =   =  1+ M1 1 + M2
  T1 T2  k + 1  2 2
 T0 

Esta ecuación se puede resolver para M 2 de la forma siguiente:

2 2
 2   k − 1 2  k − 1 2   2   1 k − 1  1 k −1 
2
M M =
1
2
2  1 + M1 1 + M2  ⇒ 1 =    2+  2 + 
 k +1  2  2   k + 1   M1 2  M 2 2 

Ing. César A. Quispe Gonzáles, M Sc. Página 112

 Flujo Compresible Capítulo VI: Ondas de choque

2
 k +1
2  
 k + 1   1 k −1   1 k −1  1  2  − k −1
  = 2 +  2 +  ⇒ =
1 k −1
 2   M1 2   M2 2  M 22 + 2
M12 2
2
 k +1  k −1  1 k −1  k −1 1 1 
kM12 −
k −1 
  −  2+  k− 2  
1  2  2  M1 2  2 M1 2 M1  2 
= = =
M 22 1 k −1 1 k −1 1  k −1 2 
2
+ 2
+ 1 + M1 
M1 2 M1 2 M12  2 
k −1 2
1+ M1 2 + ( k − 1) M12
M 22 = 2 ⇒ M2 = (6.21)
k M12 −
k −1 2 k M12 − ( k − 1)
2
Luego, sustituyendo la Ec. (5.21) en la Ec. (5.20) se podrá obtener:

T2  2 k M1 − ( k − 1)   2 + ( k − 1) M1 
2 2
T2 2 + ( k − 1) M12
= ⇒ = (6.22)
T1 2 + ( k − 1) M12 T1 2
( k + 1) M12
2 + ( k − 1)
2 k M12 − ( k − 1)

La relación de presiones en la onda puede hallarse de la ecuación de cantidad de movimiento,

de donde se tiene p1 + ρ1 V12 = p2 + ρ2 V22 .

El segundo término puede reescribirse como ρ V 2 = ρ M2a2 = ρ M2k ( p ρ ) = k p M2 ,

entonces se tendrá:

p2 1 + k M12
p1 + k p1M12 = p2 + k p2 M 22 ⇒ p1 (1 + k M12 ) = p2 (1 + k M 22 ) ⇒ = (6.23)
p1 1 + k M 22

Sustituyendo en esta ecuación, la expresión hallada para M 2 , se obtiene:

p2
=
1 + k M12
=
(1 + k M12 )  2k M12 − ( k − 1) 
p1 2 + ( k − 1) M12 2k M12 − ( k − 1) + 2 k + k ( k − 1) M12
1+ k
2k M12 − ( k − 1)

La reducción de esta expresión permite obtener:

p2 2k M12 + k M12 − k + 1 + 2 k 2 M14 − k 2 M12 2k M1 ( k M1 + 1) − k ( k M1 − 1) + ( k M1 + 1)

2 2 2 2

= =
p1 k 2 M12 + k M12 + 1 + k ( k M12 + 1) ( k + 1)
p2 2kM1 − ( k − 1)
2
p2 2k k −1
= ⇒ = M12 − (6.24)
p1 k +1 p1 k + 1 k +1

La relación de densidades en la onda puede hallarse con la ecuación de estado para el gas
perfecto y la ecuación de continuidad:

Ing. César A. Quispe Gonzáles, M Sc. Página 113

 Flujo Compresible Capítulo VI: Ondas de choque

p2 = ρ 2 RT2  p2 ρ 2 T2 ρ 2 V1 p2 T1
 ⇒ = ⇒ = =
p1 = ρ1 RT1  p1 ρ1 T1 ρ1 V2 p1 T2

Sustituyendo en esta expresión, las relaciones halladas en las Ecs. (6.22) y (6.24), se halla:

ρ 2 V1  2kM1 − ( k − 1) 
2 2

= =
( k + 1) M12
ρ1 V2  k +1   2kM12 − ( k − 1)   2 + ( k − 1) M12 
    
ρ 2 V1
= =
( k + 1) M1 2

(6.25)
ρ1 V2 2 + ( k − 1) M12

Para la relación de presiones de estancamiento se tendrá:

k k
−
p02  p02   p2   p1   k − 1 2  k −1  2k k −1   k −1 2  k −1
=      = 1 + M2   M12 −  1 + M1 
p01  p2  ISO  p1 O.C .  p01  ISO  2   k +1 k +1 2 

Al sustituir en esta expresión, la relación hallada para M 2 en la Ec. (6.21) y reduciendo, se

obtiene:
k 1
−
p02  ( k + 1) M1  k −1  2kM12 k − 1 
2 k −1
=   −  (6.26)
p01  2 + (k − 1)M12   k + 1 k + 1 

5.7 RELACION DE RANKINE-HUGONIOT

Esta relación permite establecer la dependencia entre la presión y la densidad en una onda de
choque normal. Primero fue Rankine en Inglaterra (1851) y después Hugoniot en Francia (1887)
quienes observaron la formación de ondas con cambios abruptos de presión, que hoy se conocen
como ondas de choque.

Estas ecuaciones son deducidas a partir de las tres ecuaciones de conservación, continuidad,
cantidad de movimiento y energía; las cuales fueron utilizadas en la deducción de la Ec. (5.18). Al
sustituir en la razón de presiones en la onda de choque, expresada por la Ec. (5.18) el valor de λ
obtenido en la Ec. (5.17) que expresa la razón de densidades, se tiene:

ρ2
( k + 1) − ( k − 1)
p2 ( k + 1) λ − ( k − 1)
2
p2 ρ1
= 1
⇒ = (6.27)
p1 ( k + 1) − ( k − 1) λ12 p1 ρ
( k + 1) − ( k − 1) 2
ρ1

Al resolver la Ec. (6.27) para hallar la razón de densidades en función de la razón de presiones
se encuentra:

p2
( k + 1) + ( k − 1)
ρ2 p1
= (6.28)
ρ1 p
( k + 1) + ( k − 1) 2
p1

Esta expresión se llama adiabata de Rankine-Hugoniot y establece la dependencia de la razón

de presión respecto a la razón de densidad en la onda de choque, que difiere de la expresión
1k
adiabática e isoentrópica dada por ρ 2 ρ1 = ( p2 p1 ) .

Ing. César A. Quispe Gonzáles, M Sc. Página 114

 Flujo Compresible Capítulo VI: Ondas de choque

Si en la Ec. (6.27) se utiliza la ecuación de estado para expresar la relación de densidades

mediante la relación de presiones y de temperaturas y después se resuelve para hallar T2 T1 se
obtiene:

p2 ρ 2 T2 T2 p2 ρ1  p2   ( k + 1) + ( k − 1)( p2 p1 ) 
= ⇒ = =   
p1 ρ1 T1 T1 p1 ρ 2  p1   ( k + 1)( p2 p1 ) + ( k − 1) 
2
p p 
( k + 1) 2 + ( k − 1)  2 
T2 p1  p1 
= (6.29)
T1 p
( k + 1) 2 + ( k − 1)
p1

6.8 VARIACION DE LA ENTROPIA EN UNA ONDA DE CHOQUE

La variación de entropía puede ser expresada por:

p2 v p ρ
s2 − s1 = cv ln + c p ln 2 = cv ln 2 + c p ln 1
p1 v1 p1 ρ2

Sustituyendo en esta ecuación los valores de los calores específicos, se puede obtener:

1 k 1 k
R p k R ρ1 s2 − s1  p  k −1  ρ  k −1  p  k −1  ρ  k −1
s2 − s1 = ln 2 + ln ⇒ = ln  2  + ln  1  = ln  2   1 
k − 1 p1 k − 1 ρ 2 R  p1   ρ2   p1   ρ 2 

Al sustituir la razón de presiones hallada en la Ec. (5.18) y de densidades en la Ec. (6.17) se

tiene:

1 k 1
−
s2 − s1  ( k + 1) λ12 − ( k − 1)  k −1  1  k −1  ( k + 1) − ( k − 1) λ12  k −1 −
2k
= ln  2   2 = ln   λ1 k −1

 ( k + 1) − ( k − 1) λ1   λ1   ( k + 1) λ1 − ( k − 1) 
2
R

El término bajo el logaritmo natural es similar a la razón de presiones de estancamiento hallada

en la Ec. (6.19), por lo que se puede escribir:

s2 − s1 p
= − ln 02 (6.30)
R p01

Similar situación se tiene si se utiliza el número M. Sustituyendo la razón de presiones dada en

la Ec. (6.24) y la razón de densidades dada en la Ec. (6.25); se obtiene:

1 k

s2 − s1  2kM12 − ( k − 1)  k −1  2 + ( k − 1) M12  k −1
= ln    
  ( k + 1) M1 
2
R  k +1
1 k
−
s2 − s1  2kM12 − ( k − 1)  k −1  ( k + 1) M12  k −1
= ln    2 
R  k +1   2 + ( k − 1) M1 

Cambiando de signo la parte derecha de la ecuación se tiene:

Ing. César A. Quispe Gonzáles, M Sc. Página 115

 Flujo Compresible Capítulo VI: Ondas de choque

1 k
−
s2 − s1  2kM12 − ( k − 1)  k −1  ( k + 1) M12  k −1
= − ln    2
(6.31)
R  k +1   2 + ( k − 1) M1 

La parte derecha de la ecuación, dentro de las llaves; concuerda con la Ec. (6.26), entonces:

s2 − s1 p
= − ln 02
R p01

Cuando M 1 > 1 , la Ec. (6.31) muestra que cuando el flujo atraviese la onda de choque se
tendrá un incremento en la entropía, tal como está estipulado por la segunda ley de la termodinámica.
Cuando M 1 < 1 esta ecuación da como resultado una variación de entropía negativa, lo cual
físicamente no puede ser posible en los flujos adiabáticos. Por lo tanto, esta ecuación certifica que la
onda de choque normal solo será posible si el flujo incidente es supersónico.

6.9 LINEAS DE FANNO

Varios procesos que ocurren en flujos gaseosos se pueden entender mejor si son analizados
por medio de diagramas, como el diagrama h − s .

Los flujos adiabáticos en los cuales se tienen efectos de fricción, se pueden mostrar mediante
un diagrama T − s o su equivalente h − s . En este caso, los parámetros del flujo se podrán
determinar mediante la línea de Fanno, que estada basada en las ecuaciones de continuidad, de
energía y de estado y que se pueden expresar bajo la forma s = s ( h, ρ )

Considerando que para los gases perfectos se tiene:

1 dv dρ
v= ⇒ ln v = C − ln ρ ⇒ =−
ρ v ρ
el cambio de entropía puede expresarse como:

dT dv dT dρ
ds = cv +R ⇒ ds = cv −R (6.32)
T v T ρ

También de la ecuación de continuidad se halla que ρV = const ⇒ d ρ ρ = − dV V

De la ecuación de energía se puede hallar:

V2 k k 2k dV d ( T0 − T )
+ RT = RT0 ⇒ V= R (T0 − T ) ⇒ =
2 k −1 k −1 k −1 V 2 (T0 − T )

Sustituyendo estas dos últimas expresiones en la Ec. (6.32) se obtiene:

dT d (T0 − T ) ds dT R d (T0 − T )
ds = cv +R ⇒ = +
T 2 (T0 − T ) cv T cv 2 (T0 − T )

Recordando la Ecuación de Mayer, R cv = k − 1 , por lo que la expresión anterior queda:

ds dT k − 1 d (T0 − T )
= + (6.33)
cv T 2 (T0 − T )

Ing. César A. Quispe Gonzáles, M Sc. Página 116

 Flujo Compresible Capítulo VI: Ondas de choque

Al integrar la Ec. (6.33) se obtiene el cambio de entropía en función de la temperatura,

considerando que la temperatura de estancamiento permanece constante. Al graficar esta relación en
un sistema termodinámico T − s se obtiene una curva que se conoce como línea de Fanno y que se
muestra en la Fig. 6.6. Ya que el cambio de entalpia es proporcional al cambio de temperatura, la
representación del diagrama T − s será similar al diagrama h − s .

Figura 6.6 Líneas de Fanno en el diagrama h − s (T − s )

La línea de Fanno también puede ser hallada en función del número M, a partir de otra
expresión para la variación de entropía, recordando que dh = c p dT :

dT dp dh dp
ds = c p −R ⇒ ds = −R
T p T p

La cual puede reescribirse como:

 1 R dp  dh 1
ds = dh  −  ⇒ = (6.34)
 T p dh  ds 1 − R dp
T p dh

En la ecuación anterior, es necesario encontrar una expresión para dp dh , la cual puede ser
hallada de la siguiente manera:

dρ
dV
De la ecuación de continuidad: +
= 0 y de la ecuación de energía se tiene
ρ V
dh + VdV = 0 por lo que al combinar estas dos expresiones se halla: dh = V 2 d ρ ρ .

dρ dp dT
Para un gas perfecto: = − , por lo que la ecuación anterior se vuelve:
ρ p T
 dp dT 
dh = V 2  − 
 p T 

Ing. César A. Quispe Gonzáles, M Sc. Página 117

 Flujo Compresible Capítulo VI: Ondas de choque

dh dT dh
Así también, de h = c pT se obtiene = = que al sustituir en la expresión anterior
h T c pT
 dp dh 
dh = V 2  −
 p c pT 
da
 
Agrupando términos se obtiene:

dp dh  V2  dp dp p p
dh = V 2 −V 2 ⇒ dh 1 +  =V 2 ⇒ = 2+
p c pT  c pT  p dh V c pT
 
Al sustituir en la Ec. (5.34) se obtiene:

dh 1 dh 1 1
= ⇒ = =
ds 1 R  p p  ds 1
−
R
−
R 1 R ( k − 1) R
−  2+  − −
T p  V c pT T V 2 c pT T M 2 k R T k RT

Permitiendo obtener finalmente:

dh 1 kTM 2 dh kTM 2
= 2 = ⇒ = (6.35)
ds M k − 1 − M 2 ( k − 1) M 2 k − 1 − M 2 k + M 2 ds M 2 − 1
M 2 kT
Cuando se integra la Ec. (6.35) se obtiene la ecuación de la entalpia en función de la
temperatura, del número M y de la entropía. Esta ecuación puede ser graficada en un sistema
h − s obteniéndose la línea de Fanno, mostrada en la Fig. 6.6.

Se observa que la función tiene un punto de máxima entropía, el cual puede ser hallado de la
condición:

ds M 2 − 1
= = 0 ⇒ M2 −1 = 0 ⇒ M = 1 (6.36)
dh kTM 2
Por lo tanto, la máxima entropía ocurre a condiciones sónicas. Para la parte superior de la
2 2
curva dh ds < 0 y d h ds < 0 , correspondiendo a un flujo subsónico por lo que de acuerdo a la
ecuación de energía, se tiene que a menores velocidades corresponden mayores entalpías. En la
2 2
parte inferior de la curva, dh ds > 0 y d h ds > 0 ; lo que corresponde al flujo supersónico. Ya que
el flujo es adiabático, entonces de acuerdo a al segundo principio de la termodinámica la entropía
debe aumentar, por lo tanto, el flujo solo podrá moverse de izquierda a derecha en el diagrama h − s
(T − s ) . Se debe indicar también que la curva fue obtenida bajo la condición de
m A = ρV = const y que para cada gasto másico, corresponderá una curva en particular.

De acuerdo a la ecuación de energía, cuando la velocidad aumenta, la temperatura (o entalpia)

disminuye y viceversa. Ya que el flujo es adiabático, con rozamiento, entonces la entropía aumenta,
disminuyendo la entalpia, lo que origina un incremento de la velocidad. Por eso, en el flujo subsónico;
los efectos de fricción hacen que la densidad disminuya y parte de la energía de fricción es
recuperada por el flujo aumentando su velocidad, a fin de mantener el gasto másico constante,
debido a lo cual también aumenta el número M . En el caso del flujo supersónico, se observa un
fenómeno inverso, los efectos de fricción hacen que la densidad aumente a costa de la disminución
de la velocidad, por lo que el número de Mach disminuye. Cuando el flujo alcanza M = 1 , entonces el

Ing. César A. Quispe Gonzáles, M Sc. Página 118

 Flujo Compresible Capítulo VI: Ondas de choque

flujo se estrangula (se ahoga) y se mantendrá sónico a la salida a costa de la disminución del gasto
másico.

En el régimen subsónico también se observa que al disminuir el número M hacia valores muy
pequeños, la línea de Fanno se acerca asintóticamente a la línea de entalpía de estancamiento,
volviéndose casi horizontal en el extremo subsónico, lo que corresponde a un proceso con
estrangulación, en la que la velocidad permanece casi constante, pero se observan incrementos de
entropía.

Se debe indicar también que para deducir las ecuaciones de la Línea de Fanno solo fueron
utilizadas las ecuaciones de continuidad y de energía, no utilizándose la ecuación de cantidad de
movimiento, por lo que la Línea de Fanno es sólo una solución particular de un caso general.

6.10 LINEAS DE RAYLEIGH

Las curvas alineas de Rayleigh se derivan de los principios de conservación de la masa y los
principios de momento lineal. Las curvas de Rayleigh corresponden a un flujo que sufre una variación
de calor, como consecuencia de esto la entalpía de estancamiento y la temperatura de estancamiento
son variables, por lo que las ecuaciones presentadas para flujo adiabático no pueden ser aplicadas
en este caso. Para analizar este caso, se debe considerar un proceso de calentamiento o
enfriamiento simple, con la finalidad de poder despreciar los efectos de la fricción.

Anteriormente fue demostrado que para gases perfectos dv v = − d ρ ρ , por lo que el cambio
de entropía puede expresarse como:

dp dv  dp c p d ρ   dp dρ 
ds = cV + cp ⇒ ds = cV  −  = cV  −k
ρ 
(6.32)
p v  p cv ρ   p

De la ecuación de continuidad se halla que ρ V = m = const ⇒ d ρ ρ = − dV V y

V = m ρ . Con estas relaciones, la ecuación de la cantidad de movimiento se puede escribir como:

ρ 2V 2 m2 m2
p + ρV 2 = p + = p+ = const = C ⇒ p=C− (6.33)
ρ ρ ρ
Al diferenciar esta ecuación se obtiene:

 ρ C − m2   ρ C − ρ C + m  d ρ
2

d 
  
ρ ρ 2  m2 dρ
=  = 
dp dp
⇒ =
p ρ C − m2 ρ C − m2 2
p ρC −m ρ
ρ ρ

por lo que la Ec. (5.32) se puede reescribir como:

 m2 dρ dρ   m2  dρ
ds = cV  2
− k  ⇒ ds = cV  2
−k (6.34)
 ρC −m ρ ρ   ρC −m  ρ

La ecuación (6.34) expresa a la entropía como una función de la densidad s = f ( ρ ) . La

entalpia también puede ser hallada como una función de la densidad, a partir de

dh = c p dT .

La diferencial dT puede ser hallada de la ecuación de estado, que al ser diferenciada da:

Ing. César A. Quispe Gonzáles, M Sc. Página 119

 Flujo Compresible Capítulo VI: Ondas de choque

dp d ρ dT dT dp d ρ  m 2  dρ p  m2  dρ
= + ⇒ = − = 2
− 1 ⇒ dT =  2
− 1
p ρ T T p ρ  ρC −m  ρ Rρ  ρC −m  ρ

Sustituyendo el valor de p obtenido en la Ec. (6.33), la ecuación anterior queda:

1  m2   m 2  dρ 1  Cm 2 m2 m2 m2  d ρ
dT =  C −  − 1  =  − C − + 
Rρ  ρ   ρ C − m2  ρ R ρ  ρ C − m2 ρ ρ C − m2 ρ  ρ
1  ρ Cm − C ρ ( ρ C − m ) − m + m ( ρ C − m )  d ρ
2 2 4 2 2

dT =  
Rρ 
 ρ ( ρ C − m2 )  ρ

1  m ( ρ C − m ) − ( ρ C − m )( ρ C − m )  d ρ
2 2 2 2
1  2m 2 − ρ C  d ρ
dT =   ⇒ dT =  
Rρ 
 ρ ( ρ C − m 2
) 

ρ Rρ  ρ  ρ

lo que permite reescribir la ecuación de la entalpia como una función de la densidad h = f (ρ),
obteniéndose:

c p  2m 2 − ρ C  d ρ
dh =   (6.35)
Rρ  ρ  ρ

La representación gráfica de las ecuaciones que determinan la entropía s y la entalpia h como

función de la densidad ρ en un diagrama h − s genera una curva que se conoce como línea de
Rayleigh y que se muestra en la Fig. 6.7.

Figura 5.7 - Línea de Rayleigh en el sistema h − s (T − s ) .

Al igual que la línea de Fanno, la línea de Rayleigh posee dos valores de entalpia para cada
valor de entropía y la rama superior de la curva representa los parámetros del flujo en el régimen
subsónico, mientras que la rama inferior representa a los parámetros del flujo en el régimen
supersónico. Al realizar una transferencia de calor al gas (calentamiento del gas), la entropía
aumenta acercándose al punto de mayor entropía; moviéndose por la línea de Rayleigh de izquierda
a derecha, mientras que al retirar calor del flujo (enfriamiento el flujo) la entropía tiende a disminuir y
las propiedades del flujo se trasladarán de derecha a izquierda.

Ing. César A. Quispe Gonzáles, M Sc. Página 120

 Flujo Compresible Capítulo VI: Ondas de choque

De la gráfica se observa que la línea de Rayleigh tiene un valor de entropía máxima, la cual
puede ser calculada al dividir las expresiones halladas en las Ecs. (6.34) y (6.35), de donde se
obtiene:

 m2  dρ  m2 
cV  2
− k   2
−k
ρC − m  ρ = cV R ρ  ρ C − m
=   =0
ds
dh c p  2m − ρ C  d ρ c p
2
 2m − ρ C 
2

   
Rρ  ρ  ρ  ρ 

luego:

m2 ρ 2V 2
−k = 0 ⇒ =k
ρ C − m2  m2  2 2
ρ p+ − ρ V
 ρ 
ρ 2V 2 p
⇒ = k ⇒ V 2 = k = a2
2 2 2 2
ρ p+ ρ V − ρ V ρ

Se observa que el flujo alcanza su máxima entropía cuando la velocidad del flujo iguala a la
velocidad del sónico, es decir a condiciones sónicas, donde V = a; ( M = 1) .

Así mismo, se observa que la curva de Rayleigh tiene un punto donde la temperatura del flujo
alcanza un máximo valor. Para hallar este punto se procede de la siguiente manera. De la ecuación
de movimiento unidimensional (Ecuación de Euler) se tiene:

dp dp
V dV + =0 ⇒ = −V dV ⇒ dp = − ρ V dV
ρ ρ
Así también, de la primera y segunda ley de la termodinámica se puede hallar:
T ds = dh − v dp = dh − dp ρ

Si de las dos ecuaciones anteriores se elimina el término común dp ρ se obtiene:

Tds = dh − V dV .

Dado que para un gas perfecto, la entalpia es función de la temperatura y del calor específico,
entonces:

ds dV dT T
Tds = c p dT − VdV ⇒ T = cp − V ⇒ = (6.36)
dT dT ds c + V dV
p
dT

dp d ρ dT
Para hallara dV dT se diferencia la ecuación de estado, lo que da: = +
p ρ T

De la ecuación de movimiento de Euler dp = ρ V dV y del la ecuación de continuidad

d ρ ρ = − dV V , entonces se puede escribir la ecuación anterior como:

ρ V dV dV dT  V 1 dT  T V  dV dV RV
− =− + ⇒ − +  dV = ⇒  −  =1 ⇒ =
ρ RT V T  RT V  T  V R  dT dT RT − V 2

Ing. César A. Quispe Gonzáles, M Sc. Página 121

 Flujo Compresible Capítulo VI: Ondas de choque

Sustituyendo esta relación en la Ec. (6.36) se puede escribir:

dT T T T T T 1 − k M2
= = = = =
ds  RV  RV 2 RM 2 k RT k R M2 c p (1 − M 2 )
cp + V  2 
cp + 2 cp + c p + 2
 RT − V  RT − V RT − M 2 k RT 1− k M

Al igualar a cero se obtiene el máximo valor de temperatura:

1 − k M 2 = 0 ⇒ M=1 k (6.37)

La Ec. (6.37) representa la pendiente dT ds de la curva de Rayleigh en el diagrama T − s .

Cuando el número de Mach está en el intervalo 0 < M < 1 k , la pendiente es positiva y cuando
1 k < M < 1 , la pendiente será negativa, alcanzando su máximo valor para M = 1 k . Cuando
M = 1 , la pendiente dT ds es infinita. Cuando M>1 , entonces dT ds es positivo, tal como se
puede observar en la Fig. 6.7.

Cuando el flujo es subsónico, el calentamiento del gas acelera al flujo y lo lleva hacia las
proximidades de M = 1 , mientras que el enfriamiento del gas disminuye la entropía y desacelera el
flujo. Cuando el flujo es supersónico, ocurre un caso inverso; el calentamiento disminuye el número M
mientras que el enfriamiento lo eleva. Dado que el calentamiento del gas aumenta la entropía, tanto
en el flujo subsónico como en el supersónico; el flujo tiende a un mismo valor de entropía, pero sólo
en una dirección; por lo que el calentamiento no puede transformar el flujo subsónico en supersónico
y viceversa.

Para un gasto másico dado, existirá una cantidad de calor necesario para que la entropía
alcance un valor máximo y si este calor supera al valor necesario, se dice que el flujo está
térmicamente bloqueado. En el caso del flujo subsónico si la cantidad de calor excede a la cantidad
necesaria, causará la disminución del gasto másico por o que el número de Mach inicial se verá
reducido; en el diagrama h − s la curva pasará a otra línea de Rayleigh, las condiciones sónicas
seguirán existiendo en la salida sólo que la velocidad será mayor.

En el caso del flujo supersónico ocurrirá una onda de choque, pero dado que este choque es
adiabático no ocurrirá una variación en la temperatura de estancamiento T0 en la sección donde se
forma la onda, solo que el flujo el flujo se volverá subsónico y el gasto másico permanecerá
invariable.

Se debe indicar que, al igual que la línea de Fanno, la línea de Rayleigh es una solución parcial
ya que sólo fueron empleadas dos de las ecuaciones de conservación: la de continuidad y cantidad
de movimiento.

6.11 SOLUCION PARA LA ONDA DE CHOQUE MEDIANTE LAS LINEAS DE FANNO Y

RAYLEIGH

Las propiedades del flujo a través de un volumen de control que engloba una onda de choque,
puede ser descrita por las ecuaciones de conservación (de continuidad, cantidad de movimiento y
energía) así como la ecuación de estado. Dado que la línea de Fanno solo utiliza la ecuación de
continuidad, de energía y de estado, entonces es una solución parcial, que se complementará con la
línea de Rayleigh, que fue desarrollada sobre la base de la ecuación de continuidad, cantidad de
movimiento y ecuación de estado.

Los dos puntos donde la línea de Fanno interseca a la línea de Rayleigh representan a los dos
puntos límites de un proceso adiabático, que corresponderán al estado que tiene el flujo, antes de la
onda y después de ella. Cuando el flujo atraviesa la onda de choque, lo hace mediante un proceso
sumamente rápido, por lo que se puede considerar que es un proceso adiabático.

Ing. César A. Quispe Gonzáles, M Sc. Página 122

Figura 6.8 Gráfico de la solución de la

onda de choque
 Flujo Compresible Capítulo VI: Ondas de choque

De acuerdo a la segunda ley de la termodinámica, la entropía o bien aumenta o bien

permanece y dado que el punto 1 representa el estado del flujo antes de pasar la onda de choque,
correspondiendo a un régimen supersónico, entonces este punto 1 deberá encontrarse en la
intersección de las ramas inferiores de las líneas de Fanno y Rayleigh.

Por el contrario, el punto 2 representa el estado del flujo cuando deja la onda de choque, por lo
que este punto deberá encontrarse en la rama superior de las líneas de Fanno y de Rayleigh, tal
como se muestra en la Fig. 6.8.

Figura 6.8 - Gráfico de la solución de la onda de choque

6.12 ONDA DE CHOQUE NORMAL EN TOBERAS CONVERGENTE-DIVERGENTE

El flujo a través de una tobera será isoentrópico solo si la relación de la contrapresión p con la
presión de estancamiento p0 , alcanza un valor especial; en cualquier otro caso, el flujo generará
pérdidas irreversibles. Un flujo será isoentrópico en una tobera convergente-divergente si la tobera
está diseñada de acuerdo a las leyes isoentrópicas; además de que la presión en la sección de salida
es igual a la contrapresión.

Para analizar el efecto de la contrapresión, se analiza un flujo que atraviesa una tobera
convergente- divergente, tal como se muestra en la Fig. 6.9. Cuando el valor de la contrapresión se
encuentra entre los valores p2 y p4 se generará una onda de choque normal estacionaria en la
parte divergente de la tobera. En este caso el flujo se comportará como isoentrópico dentro de la
tobera sólo hasta el punto en que ocurre una súbita transición del régimen supersónico al régimen
subsónico.

Pasada la onda, el flujo nuevamente se comporta como isoentrópico a velocidad subsónica,

desacelerándose hasta que la presión en el plano de salida sea igual a la contrapresión. La
localización de la onda de choque dentro de la tobera dependerá de la contrapresión, a mayor
diferencia entre la contrapresión y la presión de diseño más cerca a la garganta se localizará la onda
de choque normal.

Como también se muestra en la Fig.6.9, si la contrapresión es igual a p4 la onda de choque se

formará en la en el plano determinado por la sección de salida de la tobera, dentro de la tobera el flujo
tendrá una presión igual a la contrapresión y pasando la onda se nivelará con la presión en la sección
de salida.

Si la contrapresión es menor que la presión p4 el flujo en la tobera se mantiene sin

experimentar cambio alguno y el ajuste hacia el valor de la contrapresión ocurre fuera de la tobera.
Cuando la contrapresión se encuentra entre los valores p4 y p6 , el flujo es supersónico a la salida
de la tobera y en este caso; se dice que la tobera está sobre-expandida. En este caso, fuera de la
tobera se forma una sucesión de ondas de choque oblicuas bidimensionales, acompañadas de ondas
Ing. César A. Quispe Gonzáles, M Sc. Página 123
 Flujo Compresible Capítulo VI: Ondas de choque

de expansión. Las ondas de choque oblicuas son más débiles que las ondas normales y el ángulo
entre la onda de choque y la dirección del flujo está dado por la potencia de la onda. Este ángulo y
esta potencia decrece a medida que la contrapresión a medida que se acerca a la presión de diseño
p6 .

Figura 5.9 - Efecto de la contrapresión sobre el flujo en una tobera convergente-divergente

Cuando la contrapresión es igual a la presión de diseño p6 la onda de choque oblicua

desaparece y el flujo se vuelve isoentrópico. Cuando la contrapresión es menor que la contrapresión
de diseño p6 se dice que la tobera está sub-expandida y una sucesión de ondas bidimensionales
de expansión y ondas de compresión ocurren fuera de la tobera.
Ing. César A. Quispe Gonzáles, M Sc. Página 124
 Flujo Compresible Capítulo VI: Ondas de choque

Por el contrario, cuando la presión es mayor que p2 , que es la máxima presión de salida que
corresponde a las condiciones sónicas en la garganta (sección de menor área), el flujo a través de la
tobera es completamente subsónico y el gasto másico del flujo dependerá de la contrapresión.

Así también, como se muestra en la Fig. 6.9, el gasto másico del flujo se incrementará a
medida que la contrapresión decrece, pero esto ocurrirá sólo hasta las condiciones sónicas (hasta
que la contrapresión es igual a la presión crítica), condición que es alcanzada en la garganta de la
tobera. Si la contrapresión continua disminuyendo más allá de la presión crítica el gasto másico
permanecerá constante e igual al máximo posible.

La Fig. 6.9 también muestra que la presión de salida pS variará con la contrapresión pC .
Cuando el flujo es subsónico en la salida, la presión de salida es igual a la contrapresión. Cuando la
tobera está sobre-expandida la presión de salida es menor que la contrapresión y cuando la tobera
está sub-expandida, la presión de salida es mayor que la contrapresión. En una tobera perfectamente
expandida (sea convergente o convergente-divergente libre de ondas de choque, o sea; trabajando a
régimen de diseño) la presión de salida es exactamente igual la contrapresión.

PROBLEMA 1.- En un flujo de aire, se conoce que los parámetros antes de una onda de choque son
V = 500 m s y temperatura T = 300 K . Determinar la velocidad y la temperatura del flujo después de
pasar la onda.

Solución.- Inicialmente, se determina el número de Mach antes de la onda:

V1 V1 500
M1 = = = = 1.4401
a1 kRT1 1.4 × 287 × 300

La temperatura de estancamiento será:

T01 k −1 2  1.4 − 1 
=1+ M1 ⇒ T01 = 300 1 + 1.44012  = 424.4 K
T1 2  2 

2 2
La velocidad crítica del flujo será: a* = kRT01 = 1.4 × 287 × 424.4 = 377 m s
k +1 1.4 + 1

Ya que en la onda de choque: V1V2 = a

*2
⇒ V2 = a*2 V1 = 377 2 500 = 284.23 m s

Además: λ2 = V2 a* = 284.23 377 = 0.754

Luego, la temperatura del flujo después de pasar la onda será:

T2 k −1 2  1.4 − 1 
=1− λ2 ⇒ T2 = 424.4 ×  1 − 0.7542  = 384.2 K
T02 k +1  1.4 + 1 

PROBLEMA 2.- Un gas tiene una razón de calores específicos k = 1.3 y antes de la onda de choque,
su velocidad es el doble de la velocidad que tiene al pasar la onda. Encontrar los números de Mach
antes y después de la onda, así como la relación de temperaturas T2 T1

Solución.- De acuerdo a la relación fundamental de la onda de choque:

V2 1 1
V1V2 = a*2 ⇒ 2V2V2 = a*2 *
= λ2 = ⇒ λ1 = = 2
a 2 λ2
Transformando los coeficientes de velocidad a sus correspondientes valores de Mach se tiene

Ing. César A. Quispe Gonzáles, M Sc. Página 125

 Flujo Compresible Capítulo VI: Ondas de choque

2λ 2 ( k + 1)
M= ⇒ M1 = 2.3529 y M 2 = 0.682
1 − ( k − 1) λ 2 ( k + 1)
Como en la onda de choque, las temperaturas de estancamiento son iguales T01 = T02 , entonces

T2 k −1 2
T02
1− λ2 T2 ( k + 1) − ( k − 1) λ2 1.3 + 1 − (1.3 − 1) 0.7071
2 2

= k +1 ⇒ = = = 1.265
T1
1−
k −1 2
λ1 T1 ( k + 1) − ( k − 1) λ1 1.3 + 1 − (1.3 − 1)1.4142
2 2

T01 k +1

PROBLEMA 3.- Antes de pasar una onda de choque, un flujo de aire tiene una velocidad
V = 800 m s , una presión p = 100 kPa y una temperatura T = 300 K . Encontrar los parámetros del
flujo al pasar la onda, así como los parámetros de estancamiento antes y después de la onda.

Solución.- El número de Mach antes de la onda será:

V1 V1 800
M1 = = = = 2.3042
a1 kRT1 1.4 × 287 × 300

Luego los parámetros de estancamiento antes de la onda serán:

T01 k −1 2  1.4 − 1 
=1+ M1 ⇒ T01 = 300  1 + 2.30422  = 618.6 K
T1 2  2 
k 1.4
p01  k − 1 2  k −1  1.4 − 1 1.4 −1
= 1 + M1  ⇒ p01 = 100 1 + 2.3042 2  = 1258.72 kPa
p1  2   2 

La velocidad crítica del flujo e halla de:

2 2
a* = kRT01 = 1.4 × 287 × 618.5 = 455.1 m s
k +1 1.4 + 1

Hallando los coeficientes de velocidad:

V1 800 1 1
λ1 = = = 1.7578 y λ2 = = = 0.5689
a* 455.1 λ1 1.7578

En la onda de choque: T01 = T02 = 618.6 K

La presión pasando la onda se puede hallar de:

p2  2k k −1   2 × 1.4 1.4 − 1 
= M12 −  ⇒ p2 = 100  2.30422 −  = 602.77 kPa
p1  k + 1 k +1  1.4 + 1 1.4 + 1 

La presión de estancamiento pasando la onda se halla de:

k 1.4
−
p2  k − 1 2  k −1  1.4 − 1  1.4 −1
= 1 − λ2  ⇒ p02 = 602.77 × 1 − 0.56892  = 731.86 kPa
p02  k + 1   1.4 + 1 

Ya que el proceso de compresión en la onda de choque no es isoentrópico, entonces:

Ing. César A. Quispe Gonzáles, M Sc. Página 126

 Flujo Compresible Capítulo VI: Ondas de choque

ρ 2 ( k + 1)( p2 p1 ) + ( k − 1) (1.4 + 1) 6.0277 + (1.4 − 1)

= = = 3.09
ρ1 ( k + 1) + ( k − 1)( p2 p1 ) (1.4 + 1) + (1.4 − 1) 6.0277

PROBLEMA 4.- Fluye aire a través de dos toberas Laval, ubicadas consecutivamente una tras de
otra. Para la primera tobera se conoce que la relación del área de salida respecto al área mínima es
AS 1 AMIN 1 = 2.5 . En algún régimen de trabajo, se observa que en la sección mínima de la primera y
segunda tobera se tiene la velocidad del sonido, debido a que una onda de choque se forma entre
estas dos toberas.. Determinar la relación de áreas mínimas AMIN 2 AMIN 1 de estas toberas.

Solución.- Las toberas son de tipo Laval (convergente-divergente) y la velocidad del sonido es
observada en las áreas mínimas, entonces AMIN 1 = A1* y AMIN 2 = A2* . La onda de choque se forma
entre ambas toberas, entonces para la primera tobera se puede escribir:

A1* AMIN 1 1
qS 1 = = = = 0.4 (*)
AS 1 AS 1 2.5

Para esta sección, también se cumple:

1 1
A*  k + 1  k −1  k − 1 2  k −1
qS 1 = 1 =   λS 1 1 − λS 1  = 0.4
AS 1  2   k +1 

Resolviendo esta ecuación, se obtiene λ1 = 1.8068 ; luego:

λ2 = 1 λ1 = 1 1.8068 = 0.5535

Para la sección pasando la onda, se puede escribir:

1 1 1 1
A* AMIN 2  k + 1  k −1  k − 1 2  k −1  2.4  0.4  0.4  0.4
qE 2 = 2 = =  λ2 1 − λ2  =   × 0.5535  1 − 0.5535  = 0.7659 (**)
AE 2 AE 2  2   k +1   2   2.4 

Dividiendo la relación (**) entre la relación (*) se obtiene:

qE 2 AMIN 2 AE 2 AMIN 2 0.7659

= ⇒ = = 1.915
qS 1 AMIN 1 AS 1 AMIN 1 0.4

PROBLEMA 5.- Un flujo de aire fluye a través de una tobera cuya área de sección varía de acuerdo a
la ley:

A 2
= 1.5 − 0.5e − x , donde −1 ≤ x ≤ 1
AMIN

La presión de estancamiento a la entrada de la tobera es p0 = 250 kPa y la temperatura de

estancamiento es T0 = 300 K . En alguna sección de la tobera, se produce una onda de choque, en la
cual la relación de densidad ρ 2 ρ 02 = 0.75 . Calcular los parámetros del flujo antes y después de
pasar la onda, la posición de la sección donde ella se forma y la variación de entropía en la onda de
choque.

Solución.- De la ley de área variable se obtiene:

x -1.00 -0.75 -0.50 -0.25 0.00 0.25 0.50 0.75 1.00

A/A MIN 1.316 1.215 1.111 1.030 1.000 1.030 1.111 1.215 1.316

Ing. César A. Quispe Gonzáles, M Sc. Página 127

 Flujo Compresible Capítulo VI: Ondas de choque

De la tabla se observa que la tobera es convergente-divergente, por lo tanto, para que se genere una
*
onda de choque, el flujo tuvo que pasar de subsónico a supersónico, en tal caso AMIN = A .
De la relación de densidades:

1
k + 1   ρ2  
k −1
ρ 2  k − 1 2  k −1 1.4 + 1 
1 − ( 0.75 )  = 0.8076
1.4 −1
= 1− λ2 ⇒ λ2 = 1 −   =
ρ 02  k + 1  k − 1   ρ02   1.4 − 1  


Entonces: λ1 = 1 λ2 = 1 0.8076 = 1.2383

El coeficiente de gasto antes de la onda será:

1 1 1 1
 k + 1  k −1  k − 1 2  k −1  2.4  0.4  0.4  0.4
q1 =   λ1  1 − λ1  =   × 1.2383  1 − 1.2383  = 0.934
 2   k +1   2   2.4 

A* AMIN 1
También: q1 = = = 2 = 0.934
, resolviendo esta ecuación se obtiene x = 0.39
A1 A1 1.5 − 0.5e − x

Por lo que la onda se originará en la sección con x = 0.39 .

Los parámetros antes de la onda serán:

k 1.4
p1  k − 1 2  k −1  1.4 − 1 1.4 −1
= 1 − λ1  ⇒ p1 = 250 × 1 − 1.23832  = 89 kPa
p01  k + 1   1.4 + 1 
T1 k −1 2  1.4 − 1 
=1− λ1 ⇒ T1 = 300 × 1 − 1.23832  = 371.7 K
T01 k +1  1.4 + 1 

Los parámetros después de la onda, considerando T02 = T01 serán:

T2 k −1 2  1.4 − 1 
=1− λ2 ⇒ T1 = 300 × 1 − 0.80762  = 267.4 K
T01 k +1  1.4 + 1 
p2 ( k + 1) λ1 − ( k − 1) (1.4 + 1)1.23832 − (1.4 − 1)
2

= ⇒ p2 = 89 × = 163.37 kPa
p1 ( k + 1) − ( k − 1) λ12 (1.4 + 1) − (1.4 − 1)1.23832
k 1.4
−
p2  k − 1 2  k −1  1.4 − 1  1.4 −1
= 1 − λ2  ⇒ p02 = 163.37 × 1 − 0.80762  = 244.4 kPa
p02  k + 1   1.4 + 1 

La variación de entropía que ocurre en la onda de choque será:

s2 − s1 p  244.4  J
= − ln 02 ⇒ ∆s = 287  − ln  = 6.51
R p01  250  kg·k

PROBLEMA 6.- El proceso de compresión de un gas en una onda de choque no es isoentrópico. En

ondas de choque suficientemente débiles, este proceso es cercano a un proceso isoentrópico. Para
que valores de número M, antes de la onda de choque normal, la relación de densidades ρ 2 ρ1 se
puede determinar mediante la relación isoentrópica, sin incurrir en un error superior al 1%,
considerando que la razón de calores específicos del gas es k = 1.3

Solución.- De acuerdo a la relación de Hugoniot para la onda de choque normal, se tiene:

Ing. César A. Quispe Gonzáles, M Sc. Página 128

 Flujo Compresible Capítulo VI: Ondas de choque

ρ 2 ( k + 1)( p2 p1 ) + ( k − 1)
=
ρ1 ( k + 1) + ( k − 1)( p2 p1 )
1
ρ   p k
Del proceso isoentrópico se puede obtener:  2  =  2 
 ρ1  ISO  p1 

La relación de densidades ρ 2 ρ1 en la onda de choque, siempre es menor que en el proceso

isoentrópico y considerando la condición del problema, se puede obtener la siguiente relación:

1
ρ2 ρ 
= 0.99  2  ⇒
( k + 1)( p2 p1 ) + ( k − 1) = 0.99  p2  k
 
ρ1  ρ1  ISO ( k + 1) + ( k − 1)( p2 p1 )  p1 

Resolviendo esta expresión para p2 p1 se obtiene que p2 p1 = 2.07

Como para la onda de choque se tiene:

p2  2 k k −1  2 × 1.3 2 1.3 − 1 

= M12 −  ⇒  M1 −  = 2.07 ⇒ M1 = 1.395
p1  k + 1 k +1  1.3 + 1 1.3 + 1 

PROBLEMA 7.- A través de una tobera convergente-divergente, fluye aire con presión de
estancamiento 500 kPa y temperatura de estancamiento de 450 K. El área de la garganta de la
tobera es 5 cm2 y en la sección de salida es 12 cm2. Si la contrapresión es de 170 kPa y en alguna
sección de la parte divergente ocurre una onda de choque. Determinar: a) el número M en la sección
de salida, b) la variación de la presión de estancamiento que ocurre en la onda, c) los números M
antes y después de la onda, d) el área de la sección donde ocurre la onda de choque, e) la
contrapresión si el flujo fuera totalmente isoentrópico.

Solución.- a) Para flujos adiabáticos, a

ambos lados de la onda y considerando que
q* = 1 , se puede establecer la siguiente
relación:

p01 A1*q1* = p02 A2* q2* ⇒ p01 A1* = p02 A2*

Si se divide ambos lados de la ecuación por

pS AS
* *
p01 A p A pS As p A
1
= 02 2
⇒ = S s
pS AS pS AS p02 A2* p01 A1*

pS As 250 12
Sustituyendo valores se obtiene: = × = 1.2
p02 A2* 500 5

k 1 1
− −
p A p  k − 1 2  k −1  k + 1  k −1 1  k −1 2  k −1
De otro lado, se tiene: S s* = S qS−1 = 1 − λS    1− λ = 1.2
p02 A2 p02  k +1   2  λS  k + 1 S 

Desarrollando esta expresión, se tiene:

1
 k + 1  k + 1  k −1  k +1
 = 0 ⇒ λS + 11.357575λS − 6 = 0
2 2
λ +S   1.2λS − 
 k − 1  2   k −1 

Ing. César A. Quispe Gonzáles, M Sc. Página 129

 Flujo Compresible Capítulo VI: Ondas de choque

Resolviendo esta ecuación y tomando el valor positivo se obtiene: λS = 0.5057 o también

M S = 0.4719 .

b) La presión de estancamiento pasando la onda será:

k 1.4
p02  k − 1 2  k −1  0.4  0.4
= 1 + MS  ⇒ p02 = 250 1 − 0.47192  = 291.18 kPa
pS  2   2 

La diferencia de presiones de estancamiento será:

∆p0 = p01 − p02 = 500 − 291.18 = 208.82 kPa

c) Para hallar los números de Mach antes y después de la onda se hace:

1
2 k  k + 1 − k − 1 λ 2  k −1
p02
= ( λ1 ) k −1 
( ) ( ) 1
 = 0.5824
 ( k + 1) λ1 − ( k − 1) 
2
p01

Resolviendo esta ecuación, sabiendo que el flujo antes de la onda es supersónico se obtiene
λ1 = 1.7571 y en consecuencia λ2 = 1 λ1 = 1 1.7571 = 0.5691

Transformando a sus valores de coeficiente de velocidad se obtiene: M1 = 2.3022 y M 2 = 0.5341

d) Para hallar el área donde ocurre la onda de choque se hace:

1 1 1 1
A*  k + 1  k −1  k − 1 2  k −1  2.4  0.4  0.4  0.4
q1 = 1 =   λ1  1 − λ1  =   1.75711 − 1.75712  = 0.4551
A1  2   k +1   2   2.4 

Luego, el área será: A1 = A1* 0.4551 = 5 0.4551 = 10.99 cm 2

e) Para el caso de un flujo totalmente isoentrópico:

1 1
A1* 5  k + 1  k −1  k − 1 2  k −1
= =  λS 1 − λS 
AS 12  2   k +1 

Sabiendo que λ1 antes de la onda es supersónico se obtiene λ1 = 1.7917 permitiendo hallar:

k 1.4
p SISO  k − 1 2  k −1  0.4  0.4
= 1 − λS  ⇒ p ISO
S = 500 1 − 1.7917 2  = 34.27 kPa
p01  k +1   2.4 

PROBLEMA 8.- Un flujo de aire con parámetros de estancamiento p0 = 300 kPa y T0 = 400 K pasa
por una tobera convergente divergente, que en la garganta tiene AMIN = 10 cm 2 y en la salida
AS = 20 cm2 . En la sección donde A = 15 cm 2 se forma una onda de choque. Encontrar los
parámetros del flujo en la sección de salida y la variación de entropía en la onda de choque.

Ing. César A. Quispe Gonzáles, M Sc. Página 130

 Flujo Compresible Capítulo VI: Ondas de choque

Solución.- En la garganta el flujo es sónico, E

*
entonces A1 = AMIN AT = 10 cm2 1 2

Para la sección antes de la onda:

p01=300 kPa
*
A 10
1 T01=400 K AS = 20 cm2
q1 = = = 0.6666
A1 15

También para esta sección se tiene: A = 15 cm2

1 1
 k + 1  k −1  k − 1 2  k −1
q1 =   λ1  1 − λ1  = 0.6666
 2   k +1 

Resolviendo esta ecuación se obtiene: λ1 = 1.5635 y λ2 = 1 λ1 = 1 1.5635 ⇒ λ2 = 0.6396

Para la sección después de la onda se tiene:

1 1 1 1
 k + 1  k −1  k − 1 2  k −1  2.4  0.4  0.4  k −1
q2 =   λ2 1 − λ2  =   0.6396 1 − 0.63962  = 0.8456
 2   k +1   2   2.4 

A2*
Así también q2 = = 0.8456 ⇒ A2* = 15 × 0.8456 = 12.68 cm 2
A2
1 1
A* 12.68  k + 1  k −1  k − 1 2  k −1
Para la sección de salida: qS = 2 = = 0.6342 =   λS 1 − λS 
AS 20  2   k +1 

Resolviendo esta ecuación, se obtiene: λS = 0.4357

Como T02 = T01 = 400 K y para la presión de estancamiento:

1 1
2 k  k + 1 − k − 1 λ 2  k −1
p02 ( ) ( ) 1 2.8 
2.4 − 0.4 × 1.56352  0.4
= ( λ1 ) k −1   ⇒ p02 = 300 (1.5635 ) 0.4   = 236.51 kPa
 ( k + 1) λ1 − ( k − 1) 
2 2
p01  2.4 × 1.5635 − 0.4 

La presión estática será:

k 1.4
pS  k − 1 2  k −1  0.4  0.4
= 1 − λS  ⇒ p2 = 236.51 ×  1 − 0.4357 2  = 211.34 kPa
p02  k + 1   2.4 

La variación de entropía que ocurre en la onda de choque será:

s2 − s1 p  236.51  J
= − ln 02 ⇒ ∆s = 287  − ln  = 68.25
R p01  300  kg·k

Ing. César A. Quispe Gonzáles, M Sc. Página 131