Cap tulo X: Flujos Compresibles Estacionarios: Segunda Parte

INTRODUCCION Repaso de conceptos de la clase anterior. OBJETIVOS DE LA CLASE DESARROLLO

Flujo adiab atico con fricci on en un conducto de secci on constante. Diagramas de Fanno. Aplicaciones. Ondas de choque Normal: Relaciones. Ondas de choque oblicuas.

CONCLUSIONES

()

1 / 1

Objetivos:

Extender el an alisis de ujos compresibles a ujos en un conducto de secci on constante, donde la fricci on no puede ser despreciada. Presentar las relaciones que vinculan entalp a especica con entrop a especica en ujos compresibles. Breve an alisis de las ondas de choque.

()

2 / 1

Flujos en conductos de secci on constante.

Flujo adiab atico con fricci on Analicemos el ujo en un conducto tan largo que no podamos despreciar la fricci on del gas contra la pared. Supongamos adem as que no existe la posibilidad de intercambio de calor entre el gas y las paredes del conducto. Si la velocidad en el conducto es del orden de la velocidad del sonido o superiores, el ujo ser a turbulento. Sin embargo, consideramos un s olo aspecto de la turbulencia: supondremos que los perles de velocidades son constantes, u = cte en una secci on. Luego, caudal A u equivale al caudal real en la secci on.

()

3 / 1

Flujos en conductos de secci on constante.

Flujo adiab atico con fricci on Como el caudal es constante a lo largo del conducto, y estudiamos secciones constantes, entonces: j = u = constante = j (x) donde x es la coordenada sobre el eje del conducto. Como no hay intercambio de calor con el medio, la ecuaci on de la energ a se escribe seg un: h1 + u2 j2 u2 1 = h2 + 2 = const = h + 2 2 2 2 (1)

Consideremos la entrop a del gas. Debido a la fricci on interna, no permanece constante sino que crece con la evoluci on del gas en el conducto. Luego,
()

ds >0 dx
4 / 1

Flujos en conductos de secci on constante.

Recordando la relaci on de Gibbs T ds = dh d (2)

reemplazando sobre (??) y derivando respecto a x: T si desarrollamos 1 dp j 2 d ds + + dx dx 3 dx dp + dx s ds dx (3)

d = dx p

(4)

Reemplazando en (??) tendremos: T j 2 3 s ds 1 = dx 1 j 2 2 p dp dx (5)

()

5 / 1

Flujos en conductos de secci on constante.

Conforme a las relaciones termodin amicas s =

p

<0 p Cp

p :

factor t ermico is obaro, coeciente de compresibilidad t ermica.

Recordando que ds/dx aumenta, la expresi on del primer miembro en (??) resulta una cantidad positiva. Por otro lado, en el segundo dp miembro, el signo de debe coincidir con el de dx 1 Entonces si M a < 1 = M a > 1 =
()

j 2 2 p

= (1
p

u2 ) = M a2 1 c2

dp <0 dx dp >0 dx

(6)
6 / 1

Flujos en conductos de secci on constante.

An alogamente, como j = u = cte. en las secciones, si derivamos respecto a x: d du u + dx dx d du Como y u son cantidades positivas, el signo de es contrario al de . dx dx ds Introduciendo la expresi on de que obtuvimos en (??) en (??) dx d 1 dp ds = 2 = 1 dx c dx p Cp dx Resumiendo, podemos concluir que si: du >0 M a < 1 = dx du M a > 1 = <0 dx (M a2 1)
T p Cp c2

+ M a2 d <0 dx d >0 dx

1 dp c2 dx

Es decir que la velocidad crece aguas abajo en movimiento subs onico y decrece en movimiento supers onico, a la vez que la densidad decrece y crece respectivamente.
() 7 / 1

Flujos en conductos de secci on constante.

Estudiaremos la forma de las curvas de la entrop a en funci on de la presi on. A partir de (??) puede deducirse una relaci on para la variaci on de la entrop a con la presi on. 1 M a2 1 ds = (7) dp T j 2 3
s p

de aqu , vemos que para M a = 1, la entrop a presenta un extremo.Las curvas s(p) tienen una forma:

Podemos estudiar las evoluciones a partir de la Figura.

() 8 / 1

Flujos en conductos de secci on constante.

Sea la velocidad a la entrada del conducto u, subs onica. Como debe cumplirse que dp ds > 0 y la relaci on deducida en (??), < 0, se dene una regi on donde u < c. dx dx ds La curva no puede recorrerse en su totalidad pues a partir del m aximo, < 0. dx Luego, el ujo que ingresa a un conducto de secci on constante subs onico, permanece subs onico.

()

9 / 1

Flujos en conductos de secci on constante.

Por otro lado, si el ujo ingresa supers onico, el desplazamiento a lo largo de la rama de la izquierda del gr aco describe el ujo. Supongamos que el m aximo de entrop a se alcanza para x = lK . Si el largo del conducto l es menor a lk , el ujo permanece supers onico en todo el largo. Si el largo l > lk , nuevamente, no es posible continuar la evoluci on de la curva. En este caso, se desarrolla un onda de choque que permite que el ujo pase a la regi on ds > 0. subs onica y que cumpla dx El ujo alcanzar a nalmente en el extremo del conducto la presi on exterior correspondiente.

()

10 / 1

Flujos en conductos de secci on constante.

Para analizar la relaci on entre la entalp a y la entrop a, recordemos las relaciones: u2 2 + h = C1 = cte j = u = C2 Luego, 1 2 C2 2 + h = C3

Es decir que para una densidad de ujo dada j , se tiene h = h(). Como s = s(, h) podemos construir para distintas j diagramas h s que son llamados diagramas de Fanno.

Al igual que con las toberas, el caudal m asico se encuentra acotado por el fen omeno de bloqueo o atoramiento correspondiente a la condici on M a = 1. () 11 / 1

Diagramas de Fanno en toberas convergentes.

()

12 / 1

Diagramas de Fanno en toberas de Laval.

()

13 / 1

Onda de choque Normal. Analicemos la onda de choque de la gura

Considerando la ecuaci on de continuidad: m = 1 A1 V1 = 2 A2 V2 La ecuaci on de conservaci on de la cantidad de movimiento nos permite escribir: X Fx = p1 A1 p2 A2 = m (V2 V1 ) Como el espesor de la onda es muy peque no, A2 reemplazando obtenemos: p2 p1 = 1 V12 2 V22 que es el salto de presi on a trav es de la onda.
() 14 / 1

A1 y entonces (8)

Onda de choque Normal.

Considerando que el proceso es adiab atico, hemos visto que 2k p1 p2 V12 V22 = k 1 1 2 Resolviendo, surgen las relaciones de la onda de Choque. p2 p1 T2 T1 V2 V1 = = = 2kM a2 1 (k 1) k+1 2 [2kM a2 1 (k 1)][2 + (k 1)M a1 ] 2 2 (k + 1) M a1 (k 1)M a2 1 1 +2 = 2 (k + 1)M a2 1

(9)

()

15 / 1

Onda de choque Normal. Para un gas perfecto diat omico, podemos entonces sintetizar las ecuaciones en el siguiente gr aco:

()

16 / 1

Onda de choque Normal.

Recordando la relaci on de Gibbs (??) y las relaciones term odin amicas de gases perfectos: k p dh = d k1 dT 1 ds = Rd + Cv T Surge al integrar que: s2 s1 = Cv ln Para p2 /p1 p2 p1 2 1 k ! = Cv ln p2 p1 (k 1)p2 /p1 + k + 1 (k + 1)p2 /p1 + k + 1

1 + con peque no, esta expresi on tiene como l mite: 3 s2 s1 k 2 1 p2 p1 = 2 Cv 12k p1

(10)

Que indica que las ondas de choque son de compresi on, necesariamente p2 es mayor que p1 .
() 17 / 1

Onda de choque Normal. Considerando la ecuaci on de conservaci on (??), la reescribimos por comodidad 1 expresando vol umenes espec cos v = en vez de densidad. j2 = Como j>0 = p < 0 p > 0 v < 0 v > 0 p 2 p1 p = v1 v2 v (11)

Dado que p es positivo, necesariamente la segunda posibilidad es la u nica que puede tomar lugar. La ecuaci on de conservaci on de la energ a en estos sistemas, escrita en t erminos de v (jv1 )2 (jv1 )2 = h2 + h1 + 2 2 Reemplazando (??) en la anterior, se obtiene la llamada relaci on de Hugoniot: 1 (v1 + v2 )(p2 p1 ) = 0 2 Por otro lado, en t erminos de energ a interna e = h pv h1 h2 + e1 e2 +
()

(12)

1 (v1 v2 )(p1 + p2 ) = 0 2

(13)
18 / 1

Onda de choque Normal.

Si el gas es perfecto la relaci on de Hugoniot conduce a: (k + 1)(v 1/v2 ) (k 1) p2 = (14) p1 (k + 1) (k 1)(v1 /v2 )

()

19 / 1

Conductos de secci on constante aislados t ermicamente.

Ecuaci on de continuidad jA = m Ecuaci on de conservaci on de la cantidad de movimiento. d(u2 /2) 1 dp u2 + = dx dx 8h Ecuaci on de la energ a h+ u2 k p u2 = + = h0 2 k1 2

Para determinar la importancia de la fricci on en el ujo, comparamos el t ermino inercial frente al de fricci on: u2 u2 : 2L 8h =
L 4h L 4h

1 1

= =

fricci on despreciable fricci on dominante

()

20 / 1

Conductos de secci on constante aislados t ermicamente. Combinando apropiadamente, se llega a: k1 m h0 p2 = k A M a2 [1 + (k 1)M a2 /2] h0 T = 1 + (k 1)M a2 /2 Por lo que la presi on y la temperatura en cada posici on a lo largo del eje del conducto dependen del valor local del n umero de Mach. La distribuci on del numero de Mach surge de: 1 d 1 du 1 dA + + =0 dx u dx A dx k 1 dp k p d du +u =0 k 1 dx k 1 2 dx dx h0 = h[1 + (k 1)M a2 /2] Luego, 1 d(M a2 ) 2 du 1 dp 1 d = + M a2 dx u dx p dx dx
() 21 / 1

Conductos de secci on constante aislados t ermicamente. Se llega a la siguiente ecuaci on diferencial 2 + (k 1)M a2 kM a2 1 dx = 2 Ma 1 M a2 8h que integrada resulta: [2 + (k 1)M a2 (0)] M a2 1 kx 1 k+1 = + ln M a2 M a2 (0) 2 [2 + k + 1)M a2 ] M a2 (0) 4h Si la entrada al conducto es subs onica, p(L) = pa (15)

()

22 / 1

Conductos de secci on constante aislados t ermicamente.

Si la longitud del conducto es tal que se produce la onda de choque (condici on s onica a la salida), se supone que la onda de choque est a situada a una distancia L1 > L de la entrada, entonces: [2 + (k 1)M a2 (0)] M a2 kL1 1 1 k+1 = + ln M a2 (L1 ) M a2 (0) 2 [2 + k + 1)M a(L1 )2 ] M a2 (0) 4h [2 + (k 1)M a2 (L1 )] M a2 1 kL 1 k+1 = + ln 2 2 M a (L) M a (L1 ) 2 [2 + k + 1)M a(L)2 ] M a2 (L1 ) 4h A lo largo de la onda de choque, el salto de M a est a dado por M a2 2 (L1) = 2 + (k 1)M a2 1 (L1 ) 2kM a2 1 (L1 ) + 1 k (18)

(16)

(17)

Y el problema se cierra imponiendo la condici on de descarga a presi on ambiente. 2 h0 (0) k1 m (19) p(L) = pa = k A M a2 (L)[1 + (k 1)M a2 (L)/2]

()

23 / 1

Conductos de secci on constante aislados t ermicamente. En un conducto puede ingresar ujo supers onico si, por ejemplo, este proviene desde una tobera acoplada al mismo. Luego, M a(0) y m son determinados por las condiciones de salida en la tobera. k+1 2 2(k 1)0 c0 Ag k+1 donde Ag es el a rea de la garganta. k+1 k1 m = M a(0)0 c0 A 1 + M a2 (0) 2(k 1) 2 m = Conocidas estas condiciones, se pueden determinar las condiciones a la salida, M a(L) y p(L).

()

24 / 1

Conductos de secci on constante aislados t ermicamente.

La presencia de una onda de choque tiene asociada el crecimiento de entrop a, a un en el caso de movimientos que pueden ser considerados en todo el espacio como ujos perfectos (ausencia de efectos viscosos y termoconducci on). El incremento de entrop a implica la presencia de un mecanismo de disipaci on de energ a. Por lo que la paradoja de DAlembert en el caso existan discontinuidades como las ondas de choque no tiene m as validez ya que al movimiento de un cuerpo en este caso el uido opone una resistencia. El aumento de entrop a ejerce tambi en otra inuencia sobre el movimiento. En general, si el movimiento del gas antes del choque es potencial. Luego del mismo, se convierte en rotacional.

()

25 / 1

Onda de choque oblicua. Hab amos mencionado que una perturbaci on puntual sobre un ujo supers onico tiene consecuencias solamente aguas abajo de la misma. El ujo alrededor de un cuerpo avanzando a velocidad s onica presenta caracter sticas muy distintas respecto a uno que lo hace a velocidad subs onica. Recordando el a ngulo de Mach que se forma a partir de la consideraci on de la velocidad del ujo u y la de la perturbaci on c, se observan ondas de choque oblicuas respecto a la direcci on principal del escurrimiento.

()

26 / 1

Onda de choque oblicua. En el estudio de las ondas de choque normales, determinamos que estas ajustan las condiciones externas (presi on) respecto a las condiciones de salida te oricas, isoentr opicas en conductos. Por el contrario, las ondas de choque oblicuas aparecen como consecuencia de cambios abruptos en la geometr a del ujo. Las ondas de choque oblicuas son similares a las ondas normales, en tanto que para modelizar ambas consideramos discontinuidades en magnitudes como la presi on, densidad, etc. , que tienen asociadas un aumento de la entrop a. En el siguiente esquema se estudia un caso general:
: Angulo del choque : Angulo de deecci on.

Vn1 = V1 sin( ) Vt1 = V1 cos( )

Vn2 = V2 sin( ) Vt2 = V2 cos( )

Se pueden denir las variaciones del n umero de Mach:

Mn1 = M1 sin() Mt1 = M1 cos() Mn2 = M2 sin( ) Mt2 = M2 cos( )

()

27 / 1

Onda de choque oblicua.

Las ecuaciones de conservaci on considerando un volumen de control paralelo al choque son: Conservaci on de la masa Z n V ds = 0 = 1 Vn1 = 2 Vn2
S

(20)

Conservaci on de la cantidad de movimiento. Z Z (V n V )ds =

S S

ds t

En la direcci on tangente: Vt1 (1 Vn1 S1 ) = Vt2 (2 Vn2 S2 ) La que considerando (??) conduce a Vt1 = Vt2 En la direcci on normal: Conservaci on de la energ a h1 + h1 +
2 V1

(21)
2

(p1 + 1 Vn1 )S1 = (p2 + 2 Vn2 )S2

2 V2

= h2 + = h2 +

2
2 2 Vn 1 + V t1

2 2 Vn 1 + Vt1

2
2 Vn 1

recuperando lo hallado en (??) se obtiene: V2 = h2 + n2 2 2 Entonces, las 3 ecuaciones de conservaci on guardan igual forma que las ecuaciones deducidas para ondas h1 + normales. Luego es posible analizar una onda de choque oblicua como normal siempre que consideremos que V Vn y M a M an . () 28 / 1

Onda de choque oblicua. Dado que en la mayor a de los problemas, conocemos los datos del ujo V1 y el a ngulo , resulta m as conveniente expresar en funci on de M1 y . Considerando que Vt1 r = Vt2 = Mt1 kRT1 = Mt2 kRT2 . T1 De donde, Mt2 = Mt1 Luego, como Mn2 = Mt2 tg ( ) entonces: T2 r r T1 T1 Mn2 = Mt1 tg ( ) = M1 cos() tg ( ) (22) T2 T2 Tomando en cuenta la relaci on de temperatura de choques normal, aplicada a choques oblicuos: ! 2 1 1 + k Mn T1 2 2 = 1 2 T2 1 + k Mn 1 2 Reemplazando, se obtiene:
2 Mn 2

(23)

2 M1

cos ()

1+ 1+

2 k1 Mn 2 2 k1 2 Mn 1 2

! tg2 ( ) (24)

De la relaci on para onda de choque normal para el n umero de Mach:

2 Mn 2 = 2 2 + (k 1)Mn 1 2 (1 k) + 2kMn 1 29 / 1

()

Onda de choque oblicua.

Se llega a:

" = arctg

1 sin( )

k1 k+1

sin ( ) +

!#

2 k + 1 M1

La soluci on de esta ecuaci on se representa en la siguiente gura:

Vemos que existen para y M a dados 2 soluciones: una que correponde a la onda de choque fuerte y otra a la onda de choque d ebil, de angulos peque nos y menores relaciones de presi on p2 /p1 . En la pr actica, la primera soluci on se da para p muy elevados que son improbables, luego, se considera solamente la soluci on d ebil. El gr aco muestra que hay soluci on para angulos menores que un cierto l mite. M as all a de ese valor, no es posible encontrar su correspondiente . Se dice en estos caso, que la onda del choque se separa del cuerpo. () 30 / 1

Onda de choque oblicua.

Cerca de la nariz del cuerpo, existe una zona subs onica y la onda de choque se separa del perl. Detr as de la onda de choque, s olo la componente normal de la velocidad debe ser inferior a la velocidad s onica.

()

31 / 1

Conclusiones.

Hemos estudiado la evoluci on de la entrop a para ujos compresibles en conductos largos. Apoy andonos en el segundo principio de la termodin amica, vimos que en un conducto aparece una dependencia de las caracter sticas del r egimen respecto del n umero de M a en la entrada. Bas andonos en las ecuaciones de estado, y en las leyes de conservaci on, presentamos los diagramas de Fanno que vinculan en distintos ujos la entrop a con la entalp a. Hemos cerrado formulando la relaci on que existe entre las magnitudes termodin amicas, cuando existe una onda de choque normal. La descripci on acabada de esta onda, as como la de una onda oblicua, excede los alcances del curso.

()

32 / 1