(P7) Corriente sobre álabes (Examen mayo 2018)
Una corriente bidimensional y uniforme de un gas (presión P1 ,t¡ temperatura T1 y velocidad V1 según el eje x) incide
sobre una cascada de álabes que distan entre sí una distancia H, ejerciendo sobre cada uno de ellos una fuerza (por
unidad de longitud) F = Fx ex +Fz ez . Aplicando las ecuaciones de conservación en forma integral, calcular las condiciones
aguas abajo de la cascada, P2 , T2 y V2 (módulo V2 y ángulo α) cuando dicha corriente se ha uniformizado. Supongan gas
perfecto, movimiento estacionario, fuerzas másicas despreciables y álabes aislados térmicamente. Utilizar el volumen de
control estacionario dela figura (8.11), donde la superficie entre álabes es de simetría.
1
Conservación de la masa
Z Z
d +
ρdV ρ(v − vc ) · ndS = 0
dt∀c
Sc
Como en la mayoría de los problemas de este tipo, el flujo es estacionario (primera integral nula por no depender del
tiempo) y el volumen de control está en reposo (vc = 0). Las superficies de control a considerar serán:
S1 y S2 , superficies de entrada y salida. Hay que tener especial cuidad con la de salida, dado que la velocidad
formará un ángulo con la misma.
SL envolvente lateral del fluido. Sus normales serán perpendiculares a la velocidad en todos sus puntos.
Ss superficie del álabe. Impermeable, por lo que no hay intercambios de masa ni cantidad de movimiento a través
de ella.
Es decir, que sólo tendremos en cuenta la entrada y la salida en este apartado. Además, como vamos a trabajar en
cartesianas por comodidad, obviaremos la dimensión de profundidad (no tenemos datos para usarla, además de resultar
irrelevante), por lo que las superficies las sustituiremos por las alturas L dadas.
Nos quedará
Z
−1
ρ1 V1 0 dS = −ρ1 V1 H
s1 0
Z
1
ρ2 V2 cosα V2 senα dS = ρ2 V2 cosα · H
s2 0
Uniendo ambas expresiones
(−ρ1 V1 + ρ2 V2 cosα) H = 0 (1)
Conservación de la cantidad de movimiento
Z Z Z Z Z
d
ρvdV
+ ρv(v − v c ) · ndS = − ∆P dS + τ ·ndS + ρf
mdV
dt
∀c Sc Sc Sc ∀c
Donde hemos despreciado además las fuerzas másicas por tratarse de un gas.
Además, las fuerzas viscosas son despreciables en la interfaz fluido-fluido, y en la superficie de contacto con el ala las
vamos a añadir a las provocadas por las diferencias de presión en la misma y y agruparlas en el término (−Fx , Fy )1 de
fuerza neta resultante. Los términos debidos a la diferencia de presión en la envolvente lateral también los consideraremos
nulos por ser P = Patm en los mismos.
Tendremos entonces
A la izquierda del igual
Z Z
−1 1
ρ1 V1 0 V1 0 dS + ρ2 V2 cosα V2 senα V2 cosα V2 senα dS
S1 0 S2 0
A la derecha del igual
Z Z
− (P1 − Patm ) · −1 0 dS − (P2 − Patm ) · 1 0 dS
S1 S1
Integrando ambas expresiones y agrupando para los ejes OX y OZ de acuerdo a los criterios antes expuestos, nos
quedará:
Eje OX: −Fx = −ρ1 V1 H + ρ2 V2 cos2 α − (P1 − P2 ) H (2)
Eje OZ: Fz = ρ2 V2 senα · cosα · H (3)
1 Hemos cambiado el signo de F y a positivo porque tendrá sentido opuesto al indicado en los ejes dados en el enunciado.
2
Conservación de la energía
v 2 v2
Z Z
d
(v − vc ) · ndS =
ρ
e+
dV + ρ e+
dt∀
c
2 Sc 2
Z Z Z Z Z
− P v · ndS + ρf m· vdV
+ v·τ
·ndS − q · ndS
+ QdV
∀c
Sc S
c S
c
∀c
donde hemos tenido en cuenta que la pared de los álabes es adiabática y no se genera calor interno. Además, la superficie
del álabe no permite el intercambio de propiedades que incluyan la velocidad, luego sólo integraremos en las superficies
de entrada y salida.
Nos quedaría
V2 V2
Z Z Z Z
ρ1 e + 1 (−V1 )dS + ρ2 e + 2 (V2 cosα) · ndS = − P1 (−V1 ) dS − P2 (V2 cosα) dS
S1 2 S2 2 S1 S2
Integrando y sustituyendo e por su valor e = Cv · T
V2 V2
−ρ1 V1 Cv · T1 + 1 + ρ2 V2 cosα Cv · T2 + 2 = P1 V1 − P2 V2 cosα
2 2
1 P
Usando ahora la relación entre los calores específicos, la presión, la densidad y el coeficiente γ Cv T = γ−1 ρ
V12 V22
P1 P2
−ρ1 V1 + + ρ2 V2 cosα + = P1 V1 − P2 V2 cosα (4)
ρ1 (γ − 1) 2 ρ2 (γ − 1) 2
Solución final
Suponiendo que la densidad es conocida y aproximadamente igual a la entrada y la salida del volumen de control,
podemos eliminarla como incógnita y nos quedarán cuatro ecuaciones ((1), (2), (3) y (4)) y cuatro incógnitas (P2 , V2 ,
senα y cosα).
La temperatura podemos obtenerla de la relación con γ ya dicha , eliminando Cv dado que conocemos T1
1 P1 1 P1
Cv T1 = ⇒ Cv =
γ−1 ρ γ − 1 ρT1
1 P2 P2
T2 = ⇒ T2 = T1
γ − 1 ρCv P1
