CALTIULY 10 FLUJO DE.GaSES Y VAPORES 10.1 Velocidad gel sonido y ntimero ‘de jiech iia velocidad de propagacién de una onda acuisti- ca en un medio eldstico tridimensional, se puede ex presar @) - : Vs “fF (10.1) done = ape en donde a (10.2) es:el mééulo de slesticidad volumétrico de la sustan cia que consti el medio en el. que se produce le propagacién, yp es su densidad.Ap es el pulso de presién aplicado, yAV/V la variacién relativa del volumen correspondiente. El signo: menos sé debe @ que-él pulso de presién que actta desde efuera pro- duce la disminucién Gel volumen del sistema. Para expreser la velocidad del ‘soriido en un me- dio geseoso en funcién de propiedades termodindmi- cas ¢5 neceserio determinar la relecién entre K y p> es decir el tipo de proceso que se verifica du- rante la propegacién de la onda actistica. ‘ on wton, de acuerdo con'las nociones de Boyle, planteéd le hipétesis de que el fenémeno fuese iso- eles sobre la onda aciistica ver la Dindmica de los Gases" del av- erior, 1976, Secciones 1.7 y 3.6; (@) Para "Introduccién , tor, Ed. del 219térmico. Se debe a Laplace la interpretacién correc ta de que se trata de un procéso isoentrépico, es decir adiabdtico reversible, lo cual. se. justifica teniendo en cuenta que la propagacién de la onda a~ ctstica se produce cuando la .perturbacién que provo cael pulso de presién es sumamenté débdil. El médulo de elasticidad volumétrico, ec. (10. 2), se puede expresar ; 4 ap gay dv siendo v el volumen especf{fico del gas. De acuerdo con la ley isoentrépica, ec. (2.18), p v8= cite, se deduce pee why SB = tc yh =-tove vt a=} pvi - av y remplazando K=3pv'v=tp sustituyendo K en la expresién (10.1) y teniendo en cuenta la ecuacién de estado del ges perfecto, se obtiene =/t? =Jtas (10.3) en donde se pone en evidencia que para un gas deter minado la velocidad del sonido es funcién solamente de la temperatura. Por ejemplo: para el aire se ob- tiene : s=¥1,4 x 287,JT, =.20V7 fi/s] Se llama "nimero de Mach" a la relacién entre la velocidad que se considere y la velocidad del so nido en el-lugar de la medicién, es decir = V/s : (10.4) 220La definicién' de este nimcro adimersional impli. ca que cuando la velocidad del flujo‘es la del soni, do resulta M = 1; y por-lo tanto M<1 corresponde a flujo subsénico;y M>1 a supersénico. Es evidente que los nimeros de Mach no son proporcionales ‘a las velocidades, ya.que las vélocidades del sonido no son constantes sino que varfan con la temperatura. Ejemplo 10.1: Calculer le velocidad del sonido en 2: el-.aire y én el helio a 15°C de temperatura. a) aire: Vs,= 20 YT = 209288 = 340 m/s 2 ad = 207736 ka/e ke " ~ pb) helio: “BR onoatémico) Gi ¥s =/1,66 x 2077,6 x 288 = 997 m/s_ Ejemplo 10.2:Un proyectil se desplaza a la veloci~ dad de 1000 m/s. Calcular el ntmero de Hach si la-temperaturae del aire es 15°C y -53°C. foe Nee LOCOh a) KH = Ve = BOS DeE = 3724 b) T = -'53 + 273 '= 220 K n Ma UT = a0 foo = 3937 10.2 Eeuacién de le energfa en futicién del numero de_ Hach En la secciérd 3i5:¢ se tonsideré el caso de flu “21jo estacionario de wun fluido compresible en un con- ducto adiabdético que no réaliza trabajo sobre el me dio exterior. bn tales condiciones ecuacién de 1 energia se reduce @ la expresién (3.17) y se verifi ca que el salto de entalpia es igual a la variacién de la-energfa cinética. Es importante destacar que este resultado se verifica con o sin friccién (ver secciones 3.6°y 5-9)- Si-por un conducto se descarga un gas desde un tanyue en donde la presién Pg-se mantiene constante (Fig. 3.7), la écuacién de la energia se pucde ex- presar yr : : : mtt7* i) (3.19) en donde Cp es el calor especifico a presién cons- gante del gas, T su temperatura estdtica o del flu- jo en movimiento, To su temperetura de estancamien- to 7 V la velocidad de la corriente. Por medio de le ecuacién de Mayer Cp - Cv = R, dividienéo por Cv se deduce R Cn 5 = = a ‘ ae le ¥-1 ade donde y como .¥= Cp/Cv resulta Dividiendo la ee. (3-19) por Cp 7, se obtiens remplazando Cp: se puede expresar Bo _ tele Vie 6 aoe taty © en. cuenta la definicién del -nimero ge Hach, ec. (10.4), se ovtiene Finalnente (10.5) expresién adimensional que se-conoce.como "ecuacién de le energfa en funcién del nimero de Mach", y de acuerdo con lo sefialedo. anteriormente, se: verifica en el flujo adiabdético con o sin friecién. Es: decir que. es vélida cuando se mentienen-las condiciones isoentrépicas y también. cuando se producen-en el flujo pérdidas irreversibles por efectos frieciona- les con eumento de entropfa: f For medio de la ec. (2.26) paran =} y es decir cuando se verifica ld ley isoentrépica de Poisson, ec. (2.18) pv'=cte, se suede expresar ¥ Bae (1 4 ER yt jr (20.6) Esta relacién entre la presién de impacto o de estancanie: Po y la est&tica p, se cumple solamen te. si el flujo es isoentrépico, condicién que se po ne en evidencia en el exponente. Esto significa que ta relacién (10.6) no representa a la ecuacién de ls erergia, sino gue su validez es mds restririgida ¥ se verifica solamente en el flujo isoentrépico, es decir ‘en el caso adiabdtico y ademdés ceversibdl 10.3 Descarga de gases por tcberas dad obtener la ma- iujo para una determi Una tobera tiene como fi yor velocidad de ‘se2 223nada diferencia de presién, es decir convertir la energia disponible a'le entreda en cnergie cinética ala salida. Lo contrario es un difusor, de modo que produciré la dasaceleracién de la corriente. E- xisten dos tipos bdsicos de toberas: la convergente 0 subsénica (Fig. 3,6,a) y le convergente-divergen te, llamada también tobera Laval o supersénice (Fig. 3.6eb). A mediados del siglo B= pAV XIX Pliegner en Alemania verificé en forma expéri mental que en la descar- 8a de gases por medio de toberas, que disminuyen: do la pregitn de salida él caudal en masa aumen~ ta al principio, pero cu endo la presién de des+ Figura 10.1 Carga es eproximadamente la mitad de la del tanqueyel caudal alcanza un maxi mo y se mantiene constante aunque le presién siga disminuyendo y tienda a cero. En le figura 10.1 se ‘representa el caudal mésico descargedo en-funcién dé la relacién p/py, siendo-p la présién de salida ¥. Po le del tanque, o sea la presién de estancamien to. Este valor -de B/Po para el que se alcanza el ca udal m&ximo se denomina "relacién de presidén eriti~ ca". Para determinarlo partimos de la eécuacién de la .energia para flujo adiabético, ec. (3:19), de-la cual se obtiene Vi ! I ! é P/Po =0 = 2 Cp (fo - 7) = 2 Cp To (1 -F Si no se tiene en cuenta la friccién viscosa - la que de acuerdo con la teorfa. de la capa limite 224ndtl queda reduciéa a un peguefio esposor a lo go de las paredes del conducto - la entropia del flujo en la tobera.se-considera constante y por con siguiente se puede expresar T/¥o en funcién de p/Po por. medio de relaciones isoentrénicas, es decir Voe2 . to [: -(2) Intresa obtener el producto PV,.0 sea el caudal por unidad de 4 rea o "intensidad de . caudal", en funcién de ; ° OPderig + la velacién p/pg. En ocrit la figura 10.2.se ob- Figura 10.2 serva que para baja ve locidad'el producto. pV comienza aumentando debido al incremento dela velocidad, manteniéndose prdc- ticamente le densidad constante. Si la presién de’ descarga se enula, en el case limite de una expansién al vecfo, la densidad tam- bién tiende a cero. Como la velocidad méxima posi- ble en un conducto tiene un limite, que esté deter- minedo por la eritalpfa dé estancamiento, 1o cual se pone en evidencia en-la ec. (3.19) Vmdx =/2 Cp fo, se concluye gue el producto pV también tiende a ce ro. Se observa en la figura 10.2 que, en consecuen- cia, la funcién representada deberd pasar por un md ximo,al que le corresponderd la méxima intensidad de caudal. Evidentemente en un.conducto de seccién varia- ble el imo caudal nor unicad de drea se verifica en le minima seccién transversel, es decir en la “ garganta". Pera obtener el valer P/By que le corres.ponde,.o sea la relécién de presidn critita, se de- be igualar a cero la defivade de la funcién pV con respecto a la varigble P/Po: con iveda segunda negativa. Teniendo en cuenta relaciones isoentrépicas, se deduce 1 . 2 ee 2 no a p= 02 (= ¥ por lo tento = & oe eas Po, y remplazando en le anterior, se cbtiene iz Hi 2 cp to ph, @) “ G)"| (20.7) que es la funcién buscada y se denomina "ecuacién es. el producto p mente lo mis~ tener en cu de Saint-Venant". Es evidente qui V pasa por un maximo, correspondi mo ocurre con el radicando. Luego, ch be © ite las conét €s que depend condicio-~ nés en el tanque, serd necesd acero la derivada de la funcién Bel ‘ -xF en donde. x = p/Do es decir ay. 2 ral _ ee abe 0 ax ¥ x e x. y realizando-operaciones, sé deduce 8 e 2(2 -(4)F TOs a (Pacts YHL ( ) con lo que se pone en evidencia que la relecién: de presién crftica se obtiene en la gargante del con- ducto y d spente solezente del expenente isoentrépi-(¥= 1,4) dicha relecién vale U,528.. barge veyor ve agua saturedo o sobrecalentudo aproximadamente 0,54 y para vapor hdmedo 0,577. Vemos ahora,a denostrar que ‘en correpondencia de;la'relacién de présién crftica, es decir en la gerganta, se obtiene la velocidad del sonido o sea nimero.de Hach igual a-uno. Partimos de la rele- cién isdentrépica (10.6) y hacemos Kk = 1, con lo gue se deduce ota (ta A lo éual verifica el valor obtenido en la expresién (10.8), que es lo gue yuerfamos demostrar. 10.4 Ley de las dreas Una tobere convergente-divergente es un conduc to con variacién de su ‘seccién transversal, con un -@reea minima que sé denomina "gargenta". (Fig. 3. 6.b). El flujo en le zona convergente.es subsénico y como méximo se puede obtener la velocidad del so nido en la gargania. Fara interpretar 1 la posibilidad:de aleanzar ma yor velocidad en la parte divergente, es decir flu jo supersénico, consideremos la ectiacién de conti- nuidad (3.12): pa pva en donde el “asterisco significa condicién er{tica, 0 sea yue A* es el drea dela garganta, V* 1a velo- cided del sonido yp* la densidad del fluido en la; : seecién critica. Por medio de la ec. (10.7) rempla- zando se obtiene . ays _ (ye du lai ao Ge : (3 conoce coo "2 y de las dreas" y en la figura 10.3. Se.observa gue tie 2 intéticas; come en la rganta es M la dereche corresponde a flujo subsénico izquierda a supersénico. Es evidente que Dae 4/a® representa la forma de una tobera convergente-divergente, y que los puntos A y B co- rresponden 2 velocidad supersénica y subsénica, reg ivenente. A AAS Figura 10.3 Se aclera que en la figura se han utilizado dos ejes de abscisas P/Po ¥ Hy -ye que ambos pardémetros estén univocamente vinculados por medio de le rela- cién iscentrénica (10.6). Ia ley de las dreas se Puede enunciar de le siguie: ay 226"bara una miswa.relecién de dreas se obtienép dos posibles condiciones en el- flujo isoentrépico: una a velocidad subsénica y otra supersénica™. 10.5 Zobera Laval supersénica Para que se cumpla le ley de-las dreas en un conducto convergente-divergente.se deben cumplir determinadas condiciones: en primer lugar sé: tie~ ne que-alcanzar en la garganta la velocidad del so nido, es decir le relacién de presién critica; si ‘esto no se verifica'la vélocidad.aumenta hacia la gargenta hesta un. valor méximo ‘subsénico y luego er flujo se desacelera ectuando la zona divergente como difusor subsénico.. Ademds la présién de descarga deberd ser lo su ficientemente baja como-para permitir le expansién supersénica. En la figura 10.3 se verifica que pa~ ra obtener un nimero de Hach supersénico determina do se deberé adoptar una tobera Laval cuya drea de salida tiene gue cumplir la-relacién A/A* corres- pondiente, y la pregidén dé descarga deberd ser me- nor o igual que el valor determinedo en abscisas. El disefio de la: iparte divergente es sumemente importante y agui né.lo travemos: por ser un proble ma de diridmica de lés’ gases, que se.resuelve por el método de “caracteristicas hodogréficas". Es in teresante a ar que se dese evitar el despegue de la cape li a pera lo-cuel el: dngulo de diver- cane no puede ser meyor de 7 a 8 grados, por lo cual la zone -divergente es de mayor longitud que la convergente. En ‘cambio la zona subsénica es SsU- ficiente gue sea ntinua, generalmenteun arco de circunferencia o dc elipse. Es evidente que el caudal mésico m =pa V es stante en todas las:secciones del conducto. in la zona divergent aunenta el d4rea, pero como la presién tiene que disminuir corisiderablemente,se produce una reduccién correspondiente Gc densidad, de modo tal que para que el producto de los tres factores se 26, constante, la velocided debe- rd euméntar por encima de la locidud del sonido. Se aclara que si en lugar Ge disminuirse le presién de descarga manteniendo le de estencamien to Do constante, se aumentara la po del-tanque dé jando constante le presién de salida, el caudal ya no serd constante como se ha representado en la figura 10.1,.sino que aumenta (como se he in dicado en linea de trezos) debido @l eumento de densidad del fluido. El criteric de seleccién en el salto de presién di - Si la relecién entre la presién de a la de estancemiento es : D> Po Poferit corresponde uta tobera simplemente convergente - Si la relecién es Pig (= erit se requiere una tobera Laval supersénica. La velocidad de salida se obtiene en funcién de la variacién de entalpfa, tal como se he obte- nido en la seccidén 3.5.c, ec. (3.18). “ericrdo on 230eniento ho, se pucde ex de donde - h) (10.10) Esta ecvacién da la velocidad del flujo en cualquier seccién de la tobera en donde: la entalpfa sea h,. su- ponierido .condiciones isoentrépicas. En realided se produce en la corriente aumento de entropfa por fric cién viscosa y el salto de entaltia disponible resul ta menor que el isoentrépico. En la figura 10.4 se ha’ representedo la zona de trabajo en el, diegrama de Nollier: se observa que ~ las condiciones de es- tencaniento (o de entra da) corresponden en es- te caso a vapor sobreca lentado a la presién po y temperatura ty, y en la seccién de selida a la presién pg. se trata de vapor himedo. Bl pun to s representa el-esta do del vevor 2 la sali- da para flujo isoentrépico, y s' en el caso real i- rreversible. Se llama “eficiencia dele tobera" a la relacién 7 ho'= hs! 304 X = De Ry ho - hs oe es decir que es el Gociente entre el salto entélpi- co real observado y la variacién de entalpfa supo- niendo condiciones isoentrépicas. En consecuencia, la velocidad real que se alcenza a la sli se puede expresar Vsr.=/2(ho - hs!) =¥ 20;(ho - hs) 231se ha representado también el ¢ ceién erftica o gar+ ita; es evidente que le rresidn Bg se obtiene. de le relacién critica Ubicado el punto g en el diegrama, se puede obtener le velocidad en la gargenta, o sea la del sonido Vg =f 2 (ho + hg) y también su temperatura Tg (o su titulo si estuvie ra.eni la zona de’ vapor htimedo). El caudal en masa se -obtiene de la ecuacién de continuidad (3.12) en funcién de la densidad o del -volumen especffico hepa v=Al fke/s] 10.6 Ondas de chogque . La formacién de ondas de choque es caracteristi ca del flujo supersénico; debido a la elevada com- presibilidad del gas se producen a través de un.es- pesor sumamente pequefio una disminucién brusca dé velocidad y un aumento correspondiente de presién, de densidad y de temperatura. la onda acisticea constituye un caso particular de propagacién de pulso de presién sumamente débil; si la intensidad de la perturhbacién es fuerte, como en el caso de un cuerpo que se desplaza a velocidad 232supersénica, se fornan ondas’ de choque que se mani~ fiestan como una explosién. la éxplicacién de’este fenémeno.— que se puede observar en fotograffas por medio. de sistemas espe- ciales de ‘visualizacién de flujo - se basa en los procesOs sumamente irreversibles debidos a la.visco sidad y a la Condueti bt Lt dad térmica del fluido con presible. la.dinémica de los gases estudia.el comportami- ento de la corriente delante y detrds de la onda de choque, considerdndola mea discontinuidad a través de la.cual se produce un aumento de entropia por ca usa interna debido a. desorden molecular. El fenéme= no fisico en sus detalles es muy complejo, sobreto~ do por el pequefio espesor de la onda que.es del or- den del libre camino medio molecular. Si-un cuerpo se desplaza a velocidad supersénita, por ejemplo un proyectil (Fig.10 +5), todo-el sistema de ondas que se forma en el medio gase oso circundente.se-desplaza a la misma. velocidad. $i se im- pone e todo el sistema igual Figura 10.5 velocided pero en sentido con trario, serd. posible considerar a la onda como un proceso estacionario. : En la figura 10.6 se ha representado una onda estacionaria formada en'un conducto adiabdético,nor @ Hayores detalles sobre ondas de choqué en "Intro- duccién 'e le Dindémica dé log Gases", del autor, Ed. dél Interior, Caps.5, 1976. 922Corriente, que se denomi ne "onda de choaue recta". nor seré el f s, ye través de las ondas de choque recto el flujo siempre cambia de supersdénico a sub sénico. En rel las ondas @e chogue pueden ocu-~ rrir con una cierta inclinacién con respecto direccién de la corriente (como se han représenta- do en la figura 10.5) y las designa “ondas de choque oblicuas".- £1 fh fs de la onda obli- cua puede ser subsé menor ve locided. oD A través de le onda de choqve las siguientes condiciones: e cyumplen V = cte a) Ecuacién de continuidad (3.12) - La masa que fluye por unidad de tiempo es constante ya que no puede haber generacién, sug traccién o ecumulacién de masa és de la onda. Ad aunque el conducto sea convergente . © diver, secciones transversales delan— te y det de la onda pueden considerarse cons- tantes debido' a su pequefio espesor, es: decir Aa = Ao. En consecuencia, le ecuacién de continui- dad a través de la onda se puede escribir P= Ve (10.12) 34oc z ievecién de la energia (3.19) + + Cp f= Cp To 4n las secciones 3.6 y.5.9 se insistié en que es te ecuacién se verifica en ‘el flujo adiabdtico con o sin friccién, es decir que incluye la posi bilidad-de eumento de entropfa, y por consiguien te la ecuacién de la energia es vélida para las condiciones delante y detrds de la orida de choe que, de modo que se puede expresar 2 2 +Cp -= + Cp T, =Cp fo (10.13) Adem&és se pone de manifiesto que la entaipfa de estencamiento se mantiene constante a través de la onda, es decir Cp To, = Cp Tog y si se considera que el calor especifico no va~ via, se verifice to. = Tog es cecir que a través de la onda de choque la ten peratura de estancamiento permanece constanie. c)Bcuaecién de la cantidad de movimiento 10 Por medio de la ecuacién de la cantidad de mo viniento apliceda a una superficie de control que encierre a la onda, y teniendo en cuenta que As = Ag pon ser e] espesor sumamente pequefio, se obticne™ Dt MW =p, +R Ve =cte (10.14) +7 Relaciones de Fanno y de Rayleigh Una manera“de determinar las condiciones del ujo a través de la onda de choque recta, es por dio de la interseceién-de las denominades curvas. @ Ls deduccién puede consultarse en "Introduccién Dindaica de los Gasés", del eutor. ra)2 Hayleigh, las. que pueden representarse aas de coordenadas. de Fanno y en varios in la seccién anterior se mencioné que a. través onda de chogue se verifican las siguientes e- ones Relecién de Fanno Relacién. de Reyleigh continuidad pV = cte ~ de le cantidad de moviniento pt+pV -= cte ta solucién analitica qué sé obtiene por medio de las ecuaciones de la energfa y. de continuidad se denonina "relacién de Fanno". Como la ecuacién de le regia en forna expreada es. vdlida para la co- rriente adiabética con o sin friccién, se concluye que la relacién de Fanno se puede aplicar en los ca sos de flujos con friccién, pero con temperatura de estancamiento fo constante, és decir sin intercambio de calor con. el exterior. Ia solucién que se obtiene por medio-de las ecug ciones de la centidad de movimiento y de continuidad se designa “relacién de Réyleigh". Como en las ecua- encionadas no hey condicién To constante, se ue esta relacién tiene significado fisico en en un conducto con intercambio de calor con rior. : . scindiendo de la interpretacién de las rela- « Penne y de Rayleigh por. separado, se tra> solucién sinulténea que, evidentene de- be verificar las condiciones delentc > cetrds de la onda de chogue recta. Pera ld condicién delante de ra qué el tstado del gas es tnico, 236es decir que en cualquier sistema de coordenadas~se representard ‘por un solo punto. Iuego; por medio de las relaciones de Fanno 0 de Rayleigh, el estado dél gas va a seguir un camino di- ferente. Es evidente que si las curvas correspondiens tes a estas funciones se interceptan en otroppunto, en el mismo se cumplirdn de. nuevo simultdéneamente dix chas relaciones, y en el caso que en el estado inicial xistan condiciones supersénicas, la segunda intercep- cidn representard el flujo subsénico detrds de la onda. En la figura 10.7 . se han representado » | ! Rama las curvas de Fanno subsénicas y de Rayleigh; se observa que los pun tos de interseccién de terminan:un aumen to de entropia S,-' para pasar del est: do l antes de la on da, al punto 2 aque i corresponde 21 flujo 1 ! subsénico detrds del ; 1 choque. Se demuestra - us) = 8 que en los puntos de Figura 10: adxima entropia de ambas curvas es if=1; en consécuene cia, las ramas superiores corresponderdén al flujo sub sénico y las inferiores al supersénico. Un pasaje brusco del estado 2 al 1 resulta imposi ble pues estarfa en contradiccién con el segundo prin cipio, ya que significaria una disminucién de entropia en un sistema adiabdtico (Ver Secc.. 5.6) supersénicas 237: Sevdescarga vapor de agua por una to- entrada p,= 10 atm, p 2 jiento fs = 500°C, y la presién a la salida es p,= 1 atm. Se pide: a) Qué tipo de tobera corresponde. b) La velocidad isoentrépica de salida. c) La velocidad de.salida si el rend@imiento de la tobera es, = 0,9. d) El caudal de vapor descérgado si el didmetro de la ach de 4 at : (aoe acre corresponde jobera Laval (supersénica) b) Ia velocidad tsoeiirépica © @ de salida es V =/ 2(hi~he) En el diagrama de Hollier se deterninan las entalpias i} con lo cual se obtiene : v =f2(3480 - 2845) 10°, = 1127 m/s ¢) Ta velocidad para % = 0,9 v= 2, (iyha) = (2x0, 9( 3480-2845) 10° = anes a) El caudal m = A’'V/v en donde dichos Be on corresponden ‘a la seccién. minima”o garganta. . Ag= ae. 0; 785x410 cn = 7,85xi0 nt = = hg} are obtener hg en el diagrama de Mollier se tiene que determinar previamente la presién del vapor en la*seccién de garganta, es decir la-relacién de présién criciee ny w &Fy = 0,54 p= 0,54 x 10 = 5,4 atm = 0,54 liPa yen el diagrana se obtiene hy = 3270 kJ/kg Remplazando se. obtiene vg = $2(3480 ~ 3270)10 = 648 m/s Para calcular el. volumen especifico del vapor sobrecalentado se utiliza la expresién (8.20) y la Tabla 9 3 ~ 4616x673 -_ ' 544) _ %> Se ~ 0,00348 ~ 0,000486 (24) = 0,572 m/kg Renplazando 7 me PPE a = 1185x107 x648 _ 8,89 ka/e 2 0,572 e) El numero de Hach de selida pea 7 ee, 1 = 7s 7 a8 = 1,74 (isoentrépico)