La forma diferencial de las ecuaciones básicas del flujo unidimensional en canales


Save this PDF as:
 WORD  PNG  TXT  JPG

Tamaño: px
Comenzar la demostración a partir de la página:

Download "La forma diferencial de las ecuaciones básicas del flujo unidimensional en canales"

Transcripción

1 Ingeniería Investigación y Tecnología, ISSN , III.2, 53-61, 2002 (artículo arbitrado) DOI: La forma diferencial de las ecuaciones básicas del flujo unidimensional en canales G. Sotelo-Ávila y A. Cafaggi-Félix División de Ingeniería Civil, Topográfica y Geodésica Facultad de Ingeniería, UNAM s: y (recibido: octubre de 2001; aceptado: diciembre de 2001) Resumen Las tres ecuaciones diferenciales del flujo unidimensional fueron obtenidas en su forma general y presentadas en la Revista Ingeniería. Investigación y Tecnología (Vol. Il-No.3-Iulio--Septiembre-2001). En este artículo se aplican al caso particular del canal y se deducen las ecuaciones para los casos de flujo permanente, gradual y espacialmente variado. Se muestra cómo los principios de energía y de cantidad de movimiento no describen el mismo flujo con idénticos grados de precisión, lo cual sólo se logra cuando se hacen simplificaciones importantes en ambas ecuaciones. Descriptores: flujo en canales, ecuaciones básicas en la hidráulica. Abstract The particular eguations for open channel flow are obtained from the general form of tfie diferential eguations of continuity, energy and momentum for unsteady, incompressible and one dimensional flow presented in a previous numberofthis publication. The cases ofsteady, gradually andspatially varied flow are determined. Keywords: open channel flow, hydraulic basic eguations. Antecedentes Existe confusión en la literatura al utilizar las ecuaciones de energía y cantidad de movimiento en el análisis del flujo unidimensional en canales. Los trabajos de Keulegan (1942), Strelkoff (1969), Yen y Wenzel (1970), Yen (1973) y Field et al., (1988) han explicado muchos aspectos de análisis, pero aún persiste la confusión en algunas nociones fundamentales, en particular, cuando se confunde el coeficiente a de Coriolis en la ecuación de la energía con p de Boussinesq en la ecuación de cantidad de movimiento, además, es incorrecto el uso de la ecuación de energía mezclada con la pendiente de fricción, cuando de hecho la pendiente de la línea de energía es más apropiado. El análisis deyen (1970 y 1973) es en muchos aspectos más general que el que se presenta en este artículo, pero también es menos accesible para los lectores, además de que usa únicamente la ecuación de cantidad de movimiento como básica en sus desarrollos y no discute la variación que tienen a y p con el tirante, que en canales de sección compuesta es muy importante, como puede verse en Sotelo (1998). El propósito de este trabajo es el de aplicar las tres ecuaciones diferenciales del flujo unidimensional obtenidas en su forma general y

2 DOI: presentadas en la Revista Ingeniería. Investigación y Tecnología (2001), al caso particular del canal, sin tener que deducirlas. Se aplican los tres principios básicos a un volumen de control diferencial, si bien en la práctica es común utilizar únicamente dos de ellas, la de continuidad y una de movimiento, dependiendo del flujo particular de que se trate. Se hace por ello una comparación entre las ecuaciones de energía y cantidad de movimiento toda vez que son diferentes entre sí. El volumen de control VC en un canal es deformable, limitado lateralmente por la superficie de control SC, coincidente con la pared del conducto y con una superficie libre, además de dos secciones transversales al flujo separadas la distancia ds, donde s representa la coordenada curvilínea que sigue el eje de la conducción. La otra variable independiente es el tiempo t. El volumen de líquido contenido en el VC en cualquier instante es du = Ads, donde el área A es función de (s, t) y ds es independiente de ambos. La densidad p permanece constante. La velocidad media V también es función de (s, t). El tirante y se mide normal a la plantilla. Ecuaciones diferenciales de continuidad, energía y cantidad de movimiento La forma general de las ecuaciones diferenciales del flujo unidimensional presentadas por Sotelo y Cafaggi (2001) se exponen a continuación. a) Ecuación de continuidad d(pva) 5(pA) _ Q (1) ds dt El coeficiente de Coriolis a sirve para corregir el cálculo del flujo de energía cinética, y el coeficiente de Boussinesq p se usa para corregir el cálculo de la cantidad de movimiento, ambos empleando la velocidad media V. La pendiente de energía SEse define en la forma i K = s ds donde dfir es la pérdida total de energía que ocurre en la distancia ds a lo largo del conducto en cada instante. c) Ecuación de cantidad de movimiento d_ ( + \A ap z+ -^ + p A ds v P 9 2 g 2 g ds (3) 2(P -1) 2g Q dt g dt ds (4) La irreversibilidad en la ecuación de la cantidad de movimiento se explica como la energía dh consumida por el trabajo que tienen que realizar las fuerzas de fricción que se oponen al flujo en la distancia ds, es decir, representa la pendiente de fricción. dh, =S, ds 1 gpa Donde x0 es el esfuerzo tangencial medio en el perímetro de contacto P entre líquido y pared de la sección. (5) donde p, A, y V son funciones de s y f. b) Ecuación de la energía. E cuaciones d iferenciales del flujo a su p erficie libre d f P V2 ^ z + + a ds { 99 20; r v V2 ^ \_da P +(a-p) VP0 2 g Q dt + P ^ + ^ = o g dt ds ( 2 ) Ecuación de continuidad Cuando el canal es de fondo plano, x representa la coordenada rectilínea que sigue el fondo de la conducción. Dado que la densidad p es constante, la ecuación de continuidad en un canal de fondo

3 DOI: plano, donde las secciones transversales están separadas la distancia, (Figura 1) es: o d(va) da dt d Q = da dt (6) Debido a que en el planteamiento de esta ecuación de energía se utilizó la de continuidad en algunos pasos de su desarrollo, es válido sustituir la ecuación (6a) en algunos términos de dicha ecuación. Por tanto, la ecuación (9) también se (6a) escribe como sigue Ecuación de la energía ( V2> y cos0 + a yc eos 0 da Q st En la figura 1 se observa que la carga de posición hasta el centroide de cualquier sección es zg=z +(y -yg)cos0 donde z es el desnivel entre el fondo de la sección y el plano de referencia. La carga de presión también en el centroide (considerada hidrostática) es =ygcos0 pg Por tanto, la carga piezométrica en cualquier punto de la sección transversal se puede calcular como la suma de ambas +a v l» +1 ylaa+pav+ss=o (l0) Q 2g Q 2g dt g dt El término yc A representa el momento estático del área con respecto de la superficie libre. Al ocurrir el incremento dy del tirante y, se presenta un incremento de dicho momento que, de acuerdo con el teorema de ejes paralelos, vale d(yc A )-d(yc + dy) + (dy)2 - Ayc =Ady Por ser (dy)2 un término de orden superior, se puede despreciar, por tanto zg+ - z +ycos0 pg (7) o (Ayc ) = A (11) dy La pendiente longitudinal del fondo se define como S o =senq (8) donde el signo negativo se debe a que z disminuye cuando x aumenta y siempre se interpreta como sen 0, cualquiera que sea su magnitud. Sustituyendo las ecuaciones (3), (7) y (8) en la (2), la ecuación de la energía para un canal de eje recto y fondo plano se puede escribir en la forma -S0 d d + y eos 0 + a V 2 # > y es válido entonces el siguiente desarrollo d_ dt (yca) > > Í h A f Pero también es valido el desarrollo que sigue d_ dt (yca)=y da. dyg c + A-^A dt dt da. y G A G dt dy c dy dy dt y para que amkjos desarrollos sean iguales, se debe da cumplir que dyg/d y = 1, y que yg =0. Por tanto, la dt ecuación (10) se reduce a

4 DOI: c d a d Sn + y cos0 + d x f V2A c d 1 d a -Sn + UCOS g Q f V2 ^ a Q 2 g a V2 dq p V2 5A p5v _ SF =0 Q 2g Q 2g dt g dt : 12) + (pav)+se=o 2gQ dt : 15a) Por separado se realiza el siguiente desarrollo O bien 1 d 2gQ f r< i \ a 91 \ A2 1 d 2gQ (av2q) = -S0 + ycos0 1 2gQ ^ i \ / 2 \ \ji dq Q (av ) +av 2gQ ', Q2 ^ dt a Q 3 \ + SE=0 : 15bi 1 d f v2 X d r av2 "j a V2 dq a Q = Q l 2g J l 24 J Q 2 g Otro desarrollo separado es el siguiente [13) La ecuación (1 5b) es una forma general de la ecuación de la energía en canales de cualquier pendiente y es idéntica a la ecuación de energía obtenida por Field et al., (1998), para canales de pendiente pequeña si se hace eos 0=1. 1 d ípq2] _ 1 P Q2 da! 2Q dq 2gQ dt lp A J 2g Q A2 dt A dt Pj_5Q P_Q_5A g A dt 2g A2 dt donde se ha considerado que P es independiente de t. Además: dq/dt = dva/dt = A(dV/dt)+V(dA/dt). Q/A2=V/A, por tanto Ecuación de cantidad de movimiento Se sustituyen las ecuaciones (5), (6a), (7) y (8) en la (4) y se obtiene c d d f V2A -Sn + y eos g V2 dp p V2 5Q V 5A 1av _ n S, =0 2g Q 2g ga dt g dt (16) 1 d _p i r fp 2] A a 5V + V W9A 2gQ dt g A _ dt dt _P Ü JL W_dA ili 2g A dt Pero O d 2 g ) ~ o x l 29 ) viap 2 g 1 d 2gQ dt (pqv) = por tanto v^_sp_ a_ 2 g fnv2 o 3 l P^r 2<u 12<U + = + (14) g dt g A dt 2g A dt g dt Q 2g dt Sustituyendo las ecuaciones (13) y (14) en la (12), una forma más compacta de la ecuación de la energía en canales es de manera que la suma de dos de los términos en la ecuación (16) resulta d f vm V2 5p _ d ( v X r V2 a P P - P l 20 J 2 g l 0 / l 2^ J

5 DOI: Por otra parte De la ecuación (1 1) se puede escribir d rrv n _ d rrc 2 i ) p p l 9 ) { ga a ) 1 d... 1 A5x( yc)_a dy _ dy ga A2 A d x { ga y Sustituyendo la ecuación (18) en la (17) d_ V 2 "\ + vlap PT g J3 V^_dA + _1_ d_ A g A o d V, P l 2g) (18) (19) Pero Por tanto, la ecuación (22) se puede escribir -S0 c d A 9A dy 1 dq _ S, =0 ga dt 1 dq _ 1 d(va) _ V da 1dV ga dt ga dt ga dt g dt y la ecuación (22) se presenta de la forma Además, otro de los términos de la ecuación (16) se desarrolla como sigue p V2 dq _ P V2 d(va) _ Q g Q g P X l f v ^ + A ^ V ^ - ^ + P - ^ Q g t J A g g -s0 3 + i A A 2 '\ V da dv H hS, 0 ga ga dt g dt Que es idéntica a la ecuación (21) cuando ésta se emplea en canales de pendiente pequeña (cos0= 1). Casos p articu lares l xl so - yi 9a + p5_r Q g A g ^ 2g y (20) a) Flujo espacialmente variado no permanente sustituyendo las ecuaciones (19) y (20) en la (16) y simplificando se obtiene c ^ Q 1 d -S n+ y eos A P l 9A) V da 15V c S, =0 ga dt g dt (21) que es otra forma de presentar la ecuación de cantidad de movimiento en canales. Field et al.. (1998) presentan la ecuación de cantidad de movimiento en la forma Es común presentar la ecuación de continuidad (6a) en la forma..da. dv da. V + A + = 0 dt (23) La ecuación de energía en su versión (1 2) o (1 5b) y la de cantidad de movimiento en su versión (16) o (21) permanecen las mismas, pero es más común utilizar la segunda. b) Flujo gradualmente variado no permanente -S0 _L d_ A + A ya 1 dq c n S, =0 ga dt ( 22) Cualquier forma de sección Hipótesis

6 DOI: 1. No existe aportación o salida del agua en toda su longitud dq/ = Utiliza las ecuaciones de continuidad y cantidad de movimiento. 3. El valor de (3 es constante. Al sustituir la ecuación de continuidad (6a) en la de cantidad de movimiento (16), se tiene G asto creciente Hipótesis 1. dq/ es el gasto unitario q* que se aporta lateralmente en dirección perpendicular al eje del canal. 2. Utiliza los principios de continuidad y cantidad de movimiento. 3. Se considera P = a. Sn + ucosg +,V 2 ' 2# S, =0 g dt 1 Sección rectangular (Ecuaciones de Saint-Venant) Hipótesis (24) Como sólo se tiene variación con respecto a x, de la ecuación (21) y las hipótesis mencionadas, se puede escribir donde -S0 +cos0 dy a qh d CQ2 A + S, =0 (27) 1. El canal es de pendiente pequeña (eos 0= 1) y de sección transversal rectangular. A=by, donde b es el ancho del canal. 2. Utiliza las ecuaciones de continuidad y cantidad de movimiento. 3. El valor de p es 1. En un canal rectangular se define el gasto unitario q = QIb, y se sustituye en la ecuación (6a) obteniendo la ecuación de continuidad * 3 = 0 dt La ecuación (24) se puede escribir en la forma (25) además _a_ d _ ( CP" ga ^ A, da _ otq2 da ga3 dy 2aQ dq ; ga2 donde T es el ancho de la superficie libre, por lo que a d ' Q2'" ga { A, El número de Froude es av2 dy 2aQ dq g( A / T) ga2 c dy d 6V3 -S \ 29 j 5V h S c 0 g dt 1 c d y V SV 15V. S q S, =0 g g dt (26) que es la ecuación para flujo no permanente y gradualmente variado en un canal de sección rectangular. Las ecuaciones (25) y (26) también se conocen como ecuaciones de Saint - Venant. c) Flujo espacialm ente variado perm anente y se puede escribir a d Í Q2\ ga v Ay 2 dy 2aQ dq ga2 (28) Con la ecuación (28) en la (27) y despejando resulta 2aQ dq Sn - S, S, - dy ga2 (29) cos0 - F2

7 DOI: Cuando la pendiente es pequeña, se hace eos 0= 1, en la ecuación (29). G asto decreciente Hipótesis 1. dq/ es el gasto unitario q* que sale lateralmente del canal. 2. Utiliza los principios de continuidad y de energía. 3. La pérdida de energía más importante es por fricción, por lo que SE= S,. 4. Se considera p = a. De la ecuación (6a) en la ecuación de la energía (1 2), y desarrollando uno de los términos, se puede escribir / y V V2) cos0 +a 2g) a V2 dq a V2 5Q aa5q aq da _ n S, =0 (30) Q 2 g Q 2g g dt ga dt como el flujo es permanente 3Q/St = 0, de la ecuación de continuidad (6a) en la ecuación (30), y al tener solo variación en x s e tiene f c á V_ 2 '\ dq S + y cos0 +a + -aq + S, =0 (31] v Yg ga2 La energía específica se define como por lo que la ecuación (3 1) se puede escribir c, y s a c2\ aq u dq du S0 +-^-(cos0 - F ) S, =0 UX^ ^ ga Av despejando dy/, se tiene q _c _ ctq dq d y _ < ~ ^ d x co s0-f 2 (32) Cuando el canal es de pendiente pequeña, se hace eos 0=1 en la ecuación (32). Gradualm ente variado permanente Hipótesis 1. Utiliza los principios de continuidad y de energía. 2. La pérdida de energía más importante es por fricción, por lo que SE= Sf. 3. Se considera p = a. Cuando dq/ = 0, la ecuación (32) es dy _ S0 S i cos0-f2 (33) que es la ecuación dinámica. Idéntico resultado se obtiene a partir de la ecuación (29) que emplea el principio de cantidad de movimiento; cuando el canal es de pendiente pequeña, se hace eos 0= 1. V2 E =y cos0 + a ; 2 g además = (co s0 - F2), dy de _ dy de dy

8 DOI: T F igura!. Corte longitudinal y sección transversal de un canal de eje recto y fondo plano C onclusiones La forma general de cada una de las ecuaciones diferenciales del flujo unidimensional (continuidad, energía y cantidad del movimiento, previamente obtenidas), permiten deducir las respectivas ecuaciones aplicadas al canal. Con las ecuaciones generales del canal se obtienen las respectivas ecuaciones diferenciales para cada caso particular sin necesidad de volver a plantear cada vez un nuevo desarrollo. El análisis muestra, en cada caso, cuando la ecuación de la energía o de la cantidad de movimiento describe con mayor grado de precisión al flujo que se trate. R eferen cias Field W.G., Lambert M.F. y Williams B.J. (1998). Energy and Momentum in One Dimensional Open Channel Flow. Journal ofhydraulic Research, Vol. 36, No. 1. Keulegan G.H. (1942). Equation of Motion for the Steady Mean Flow of Water in Open Channels. Research Paper RP1488, National Bureau of Standards, US. Department of Commerce, Sotelo G. (1998). Algoritmo del método de Blalock y Sturm para determinar los tirantes críticos múltiples en canales compuestos. Ingeniería Hidráulica en México, Vol. XIII, No. 1, Sotelo G. y Cafaggi A. (2001). La forma diferencial de las ecuaciones básicas del flujo unidimensional. Ingeniería. Investigación y Tecnología, Vol. II, No.3, Strelkof T. (1969). One Dimensional Equations of Open Chanel Flow. Journal of the Hydraulics División ASCE, Vol. 95, No. HY3, Yen B.C. y Wenzel H.G. (1970). Dynamic Equations for Spatially Varied Flow. Journal of the Hydraulics División ASCE, Vol. 96, No. HY3, Yen B.C. (1973). Open Channel Flow Equations Revisited. Journal of the Hydraulics División ASCE, Vol. 99, No. EM5, ,

9 DOI: Semblanza de los autores Gilberto Sotelo-Ávila. Obtuvo los grados de ingeniero civil, maestro en ingeniería y doctor en ingeniería en la Facultad de Ingeniería de la UNAM. Inició su actividad universitaria en 1958 como investigador en el Instituto de Ingeniería y como profesor en licenciatura y posgrado dentro de la misma Facultad, actividades que funge actualmente. Fue jefe de la División de Ingeniería Civil, Topográfica y Geodésica, consejero técnico de la Facultad de Ingeniería y consejero académico de área en la UNAM. Es investigador nacional nivel 2 del Sistema Nacional de Investigadores del CONACYT y ha recibido los premios Universidad Nacional 1998 y Docencia en Ciencias Exactas en la UNAM, asimismo; el premio nacional Miguel A. Urquijo 1999, otorgado por el Colegio de Ingenieros Civiles de México al mejor artículo técnico publicado durante 1997 y El Consejo Técnico de la Facultad de Ingeniería de la UNAM le otorgó las cátedras especiales Javier Barros Sierra en 1999 y Mariano Hernández Barrenchea en el año Es miembro de la Academia Mexicana de Ciencias de la Asociación Internacional de Investigaciones Hidráulicas del Colegio de Ingenieros Civiles de México y del Consejo Académico del CENEVAL para el examen de calidad profesional de los ingenieros civiles. Adriana Cafaggi-Félix. Ingeniera civil egresada en 1980 de la Fcultad de Ingeniería UNAM, obtuvo el grado de maestra en ingeniería hidráulica en 1985 en la misma institución. Fue jefa del laboratorio y del Departamento de Hidráulica; actualmente es profesora de carrera en el Departamento de Ingeniería Hidráulica. Ha participado en congresos nacionales e internacionales y es autora de varios artículos e instructivos de prácticas para el Laboratorio de Hidráulica.

10 DOI: