[Flujo Uniforme]   [Fórmula de Chézy]   [Fórmula de Manning]   [Rugosidad de Manning]   [Cálculo del Flujo Uniforme]   [Cálculo del Caudal de Avenida]   [Flujo Superficial Uniforme]   [Preguntas]   [Problemas]   [Bibliografía]   •   CAPÍTULO 5:  FLUJO UNIFORME 5.1  FLUJO UNIFORME [Fórmula de Chézy]   [Fórmula de Manning]   [Rugosidad de Manning]   [Cálculo del Flujo Uniforme]   [Cálculo del Caudal de Avenida]   [Flujo Superficial Uniforme]   [Preguntas]   [Problemas]   [Bibliografía]   •   [Arriba]   El flujo uniforme se presenta únicamente en canales prismáticos. En el flujo uniforme, la profundidad de flujo, el área de flujo, la velocidad media, y la descarga (o caudal) son constantes a lo largo del canal. El término flujo de equilibrio se utiliza para describir la condición de flujo en canales no prismáticos (naturales) de sección transversal no uniforme, el cual es equivalente al flujo uniforme en canales prismáticos. En el flujo uniforme, todas las pendientes, la pendiente de fricción Sf, la pendiente de energía Se, la pendiente de la superficie del agua Sw, y la pendiente de fondo So, son constantes e iguales a una pendiente: S. Sf  =  Se  =  Sw  =  So  =  S (5-1) El flujo uniforme no permanente no existe; si el flujo es no permanente, entonces no es uniforme. Sin embargo, para el número de Vedernikov V = 1, el flujo uniforme se convierte en neutralmente estable, lo cual conduce a las ondas de rollo (Sección 1.3). Esta condición es la "inestabilidad del flujo uniforme" descrita por Chow (1959). Cuando V < 1, las perturbaciones del flujo se atenúan, y las olas de rodillo no se forman. Establecimiento del flujo uniforme Desde un punto de vista mecánico, el flujo uniforme ocurre en un volumen de control cuando la fuerza de fricción es igual a la fuerza gravitacional. En ausencia de controles de sección (Sección 4.3), todos los flujos en canales tienden a ser uniformes. En principio, a la naturaleza "le gusta el flujo uniforme." En el flujo uniforme, la característica única de la curva de gasto (es decir, la singularidad de la relación profundidad-descarga) lo califica como control. Por lo tanto, el flujo uniforme crítico es un tipo de control muy fuerte. La profundidad del flujo uniforme se conoce como profundidad normal. La Figura 5-1 muestra la formación de flujo uniforme en un canal relativamente largo. La figura superior representa el flujo normal subcrítico, con las secciones de control aguas arriba y aguas abajo. La figura central representa el flujo crítico, con las secciones de control aguas arriba y aguas abajo. La figura inferior representa el flujo normal supercrítico, con la sección de control aguas arriba únicamente. Fig. 5-1  Formación del flujo uniforme (Chow, 1959). Velocidad del flujo uniforme En general, la velocidad media del flujo uniforme se describe por la siguiente fórmula: V  =  C R x S y (5-2) en la cual C = coeficiente de fricción, R = radio hidráulico, definido como R = A/P, y x e y son los exponentes de R y S, respectivamente. Los exponentes varían de acuerdo al tipo de rugosidad (laminar, turbulenta, de transición, o mixta laminar-turbulenta) y la forma de la sección transversal (arbitraria, hidráulicamente ancha, rectangular, trapezoidal, triangular, o inherentemente estable). En la práctica, hay dos fórmulas clásicas para el flujo uniforme: (1) la fórmula de Chézy, y (2) la fórmula de Manning. Actualmente se usan también variaciones de estas fórmulas, como son las fórmulas de Chézy adimensional y Manning-Strickler. 5.2  LA FÓRMULA DE CHÉZY [Fórmula de Manning ]   [Rugosidad de Manning]   [Cálculo del Flujo Uniforme]   [Cálculo del Caudal de Avenida]   [Flujo Superficial Uniforme]   [Preguntas]   [Problemas]   [Bibliografía]   •   [Arriba]   [Flujo Uniforme]   Para derivar la fórmula de Chézy, el esfuerzo cortante desarrollado a lo largo del fondo del canal se modela con la llamada ley de fricción cuadrática: τb  =  ρ f V 2 (5-3) en la cual ρ = densidad de la masa, f = un tipo de factor de fricción (coeficiente de resistencia), y V = velocidad media. Esta ecuación es adimensional, por lo tanto, tiene una base teórica. El esfuerzo cortante desarrollado a lo largo del perímetro mojado de un volumen de control de longitud L es (Fig. 5-2): Fs  =  τb PL  =  ρ f V 2 PL (5-4) Fig. 5-2  Volumen de control para el flujo uniforme (Chow, 1959). El peso del agua contenido en el volumen de control es W. Esta fuerza gravitacional se resuelve a lo largo de la dirección del movimiento de la siguiente manera: Fg  =  W sen θ (5-5) En un canal de pequeña pendiente:  sen θ ≅ tan θ = S. Por lo tanto: Fg  =  W tan θ  =  W S  =  γ A L S (5-6) Igualando las fuerzas de fricción (Eq. 5-4) y gravitacional (Eq. 5-6): ρ f V 2 P  =  γ A S (5-7) lo cual se reduce a: f V 2  =  g (A /P ) S  =  g R S (5-8) en la que R = radio hidráulico. Resolviendo para V: V  =  (g/f )1/2 (R S ) 1/2 (5-9) V  =  C (R S )1/2 (5-10) en la cual C = coeficiente de Chézy, se define como sigue: C  =  (g/f )1/2 (5-11) Por lo tanto, el factor de fricción f en la Ec. 5-3 es:            g f  =  ______           C 2 (5-12) La Ecuación 5-10 es la fórmula de Chézy. De la Ec. 5-8 se puede derivar una variación de la fórmula de Chézy:               V 2 S  =  f   _____                 g R (5-13) lo cual equivale a:                D       V 2 S  =  f   _____   _____                  R       g D (5-14) Dada la Ec. 4-6, la Ec. 5-14 se reduce a:                D        S  =  f   _____   F 2                R        (5-15) La Ecuación 5-15 es básicamente igual que la Ec. 4-5, la cual se deriva de la ecuación de Darcy-Weisbach aplicada para el flujo en canales abiertos. Por lo tanto, la ecuación de Chézy adimensional (Ec. 5-15) y la ecuación de Darcy-Weisbach modificada (Ec. 4-5) son esencialmente las mismas. Para un canal hidráulicamente ancho, D ≅ R, y la Ec. 5-15 se reduce a:                        S  =  f F 2                        (5-16) lo cual es igual a la Ec. 4-8. La Tabla 5-1 muestra los valores correspondientes de f, el factor de fricción de Darcy-Weisbach f, y los coeficientes de Chézy. Ver también Video de Laboratorio:   Establecimiento del flujo uniforme . Tabla 5-1  Valores correspondientes de f, f, y C. Factor de fricción f f de Darcy-Weisbach Chézy C Unidades SI Unidades acostumbradas en EE.UU. 0.002 0.016 70.02 126.83 0.003 0.024 57.17 103.55 0.004 0.032 49.51 89.68 0.005 0.040 44.29 80.21 Historia de la Fórmula de Chézy Antoine Chézy nació el 1 de Septiembre de 1718 en Chalon-sur-Marne, Francia, y murió el 4 de Octubre de 1798. En 1749, trabajando en Amsterdam, Cornelius Velsen declaró: "La velocidad debe ser proporcional a la raíz cuadrada de la pendiente." En 1757, en Hannover, Alemania, Albert Brahms escribió: "La acción desaceleradora de la fricción de fondo en un flujo uniforme no sólo es igual a la acción aceleradora de la gravedad sino también es proporcional al cuadrado de la velocidad." Velsen y Brahms estaban trabajando en las leyes generales y teorías de Torricelli y Bernoulli. Chézy usó algunas de estas ideas para desarrollar su fórmula. A Chézy se le dio la tarea de determinar la sección transversal y la descarga para un canal en el río Yvette, cerca de París, pero a una mayor elevación. Desde 1769, fue recopilando datos experimentales del canal de Courpalet y del río Sena. Sus estudios y conclusiones figuran en un informe al Sr. Perronet con fecha del 21 de octubre de 1775. El documento original, escrito en Francés, se titula "Tesis sobre la velocidad del flujo en un canal," y está firmado por el Sr. Chézy, Inspector General de la Escuela de Puentes y Caminos. El documento se encuentra en el expediente No. 847, manuscrito No. 1915, de la colección que obra en la biblioteca de la escuela. En 1776, Chézy escribió otro artículo, titulado: "La fórmula para encontrar la velocidad uniforme que el agua tendrá en un canal de pendiente conocida." Este documento, el cual reside en el mismo archivo [No. 847, Ms. No. 1915], contiene la famosa fórmula de Chézy: V = 272 (ah/p)1/2 en la cual h es la pendiente, a es el área, y p es el perímetro mojado. El coeficiente 272 es aplicable al canal de Courpalet en un sistema de unidades antiguo. En el sistema métrico, el valor equivalente es: V = 31 (ah/p)1/2 Para el río Sena, el valor del coeficiente es 44. Clemens tradujo los dos documentos de Chézy al Inglés. Riche de Prony, uno de los ex-alumnos de Chézy, fue el primero en utilizar su fórmula. Más tarde, en 1801, en Alemania, Eytelwein utilizó las ideas de Chézy y De Prony para popularizar el uso de la fórmula. 5.3  LA FÓRMULA DE MANNING [Rugosidad de Manning]   [Cálculo del Flujo Uniforme]   [Cálculo del Caudal de Avenida]   [Flujo Superficial Uniforme]   [Preguntas]   [Problemas]   [Bibliografía]   •   [Arriba]   [Flujo Uniforme]   [Fórmula de Chézy]   La fórmula de Manning, en unidades SI, es:            1 V  =  ____  R 2/3 S 1/2            n (5-17) en la cual n = coeficiente de fricción de Manning, factor de fricción, o simplemente n de Manning. En el sistema de unidades acostumbradas en EE.UU. la fórmula de Manning es:            1.486 V  =  ________  R 2/3 S 1/2               n (5-18) El valor 1.486 es un factor de conversión que se deriva de la equivalencia (3.28 piés/m)1/3 = 1.48592 ≅ 1.486. Este factor es requerido para expresar la ecuación original de Manning (Eq. 5-17) en el sistema de unidades acostumbradas de EE.UU. A fin de comparar con la fórmula de Chézy, la ecuación de Manning se expresa como sigue:            1.486 V  =  ________  R 1/6 R 1/2 S 1/2               n (5-19) La relación entre los coeficientes de Manning y Chézy se encuentra comparando las Ecs. 5-10 y 5-19:            1.486 C  =  ________  R 1/6               n (5-20) La Ecuación 5-20 implica que mientras C varía con el radio hidráulico, el valor de n es constante. Esto puede ser aproximadamente cierto para los canales prismáticos (artificiales), pero por lo general no lo es para los canales naturales (Barnes, 1967). En los canales naturales, el valor de n puede variar con el nivel de la superficie del agua y la profundidad del flujo. Esto es atribuible a: Las variaciones naturales en la rugosidad del canal con el aumento del nivel del agua, incluyendo el efecto de los flujos de desborde (Fig. 2-15), o Los cambios morfológicos en la fricción total de fondo, la cual se compone de fricción de los granos y fricción de forma, a medida que el flujo varía desde un nivel bajo, a través de un nivel intermedio, hasta un nivel alto (Simons y Richardson, 1966). Fórmulas empíricas para el n de Manning Se han desarrollado varias correlaciones entre el n de Manning y el tamaño de la partícula (diámetro del grano). Williamson (1951) ha correlacionado el factor de fricción de Darcy-Weisbach f con la rugosidad relativa, produciendo la siguiente relación (Henderson, 1966):                      ks f  =  0.113 ( ____ ) 1/3                      R (5-21) en la cual ks = rugosidad de grano, en unidades de longitud, y R = radio hidráulico. Dado que f = 8 (g /C 2), la Ec. 5-21 se reduce a:                8g                 R C  =  ( ________ )1/2 ( _____ ) 1/6              0.113              ks (5-22) En unidades acostumbradas en EE.UU., la Ec. 5-22 puede ser reducida a:               1.486 R1/6 C  =  ________________             0.0311 ks1/6 (5-23) Al comparar la Ec. 5-23 con la Ec. 5-20, se puede expresar n en términos de la rugosidad de fondo de la siguiente manera (ks en pies): n  =  0.0311 ks1/6 (5-24) En términos de rugosidad relativa y rugosidad absoluta, la expresión general para el n de Manning es (Chow, 1959): n  =  [f (R/ks)] ks1/6 (5-25) lo cual implica que en la Ec. 5-24 la rugosidad relativa es una constante (0.0311). Suponiendo que la rugosidad límite pueda ser representada por el tamaño d84 de la partícula, es decir, aquél para el cual el 84% de los granos (por peso) son más finos, la Ec. 5-24 se convierte en: n  =  0.0311 d841/6 (5-26) Strickler utilizó una constante (0.0342) para la rugosidad relativa f(R/ks), y el tamaño medio de la partícula d50 como el diámetro representativo de grano, para producir la ecuación: n  =  0.0342 d501/6 (5-27) Dado que d84 > d50, se observa que las ecuaciones de Strickler y Williamson son mutuamente consistentes. La Tabla 5-2 muestra los valores de n de Manning calculados con la fórmula de Strickler (Ec. 5-27). Tabla 5-2  Valores n de Manning calculados con la fórmula de Strickler (Ec. 5-27). Tamaño medio de la partícula d50 (pies) n de Manning 0.0001 0.007 0.001 0.011 0.01 0.016 0.1 0.023 1 0.034 Historia de la Fórmula de Manning Robert Manning nació en Normandía, Francia, en 1816, y murió en 1897. En 1826, se mudó a Waterford, Irlanda, y trabajó como contador. En 1846, durante el año de la gran carestía, Manning fue reclutado por la División de Drenaje Arterial de la Oficina Irlandesa de Obras Públicas. Después de trabajar como dibujante, fue ascendido a ingeniero asistente. En 1848, llegó a ser ingeniero distrital, cargo que ocupó hasta 1855. Como ingeniero distrital, Manning leyó "Traité d'Hydraulique" de d'Aubisson des Voissons, desarrollando luego un gran interés por la hidráulica. De 1855 a 1869, Manning trabajó para el Marqués de Downshire en la supervisión de la construcción de la Bahía del Puerto de Dundrum, en Irlanda, y el diseño de un sistema de abastecimiento de agua en Belfast. Después de la muerte del marqués en 1869, Manning regresó a la Oficina Irlandesa de Obras Públicas como ingeniero asistente. En 1874 Manning ascendió a ingeniero en jefe, cargo que ocupó hasta su jubilación en 1891. Manning no recibió ningún tipo de educación o capacitación formal en la mecánica de fluidos. Su experiencia en la contabilidad y su naturaleza pragmática influyeron en su trabajo y lo llevaron a simplificar problemas. Manning comparó las siete fórmulas más conocidas de la época: Du Buat (1786), Eytelwein (1814), Weisbach (1845), St. Venant (1851), Neville (1860), Darcy y Bazin (1865), y Ganguillet y Kutter (1869). Calculó la velocidad obtenida con cada fórmula para una pendiente dada, y el radio hidráulico variando entre 0.25 m y 30 m. Para cada condición, Manning encontró el valor medio de las siete velocidades y desarrolló la fórmula de mejor ajuste. La fórmula original fue la siguiente: V = 32 [RS (1 + R1/3)]1/2 la cual más tarde simplificó a: V = C Rx S1/2 En 1885, Manning asignó al exponente x el valor de 2/3 y escribió la fórmula de la siguiente manera: V = C R2/3 S1/2 En una carta a Flamant, Manning declaró: "El recíproco de C corresponde aproximadamente al valor de n, según Ganguillet y Kutter, siendo C y n constantes para el mismo canal." El 4 de Diciembre de 1889, a la edad de 73 años, Manning propuso por primera vez su fórmula a la Institución de Ingenieros Civiles de Irlanda. La fórmula salió a la luz en 1891, en un artículo escrito por él mismo titulado "Sobre el flujo de agua en canales abiertos y tuberías," publicado en las Actas de la Institución de Ingenieros Civiles de Irlanda. A Manning no le gustaba su propia ecuación por dos razones: En primer lugar, en aquellos tiempos resultaba difícil calcular la raíz cúbica de un número y además elevarlo al cuadrado para llegar a la potencia 2/3. En segundo lugar, la ecuación era dimensionalmente incorrecta. Por tales razones, Manning desarrolló la siguiente ecuación, la cual es adimensional: V = C (gS)1/2 [R1/2 + (0.22/m1/2 )(R - 0.15 m)] en la cual m = "altura de una columna de mercurio que equilibra la atmósfera," y C = un número adimensional "que varía con la naturaleza de la superficie." Sin embargo, en los libros de texto de finales del siglo XIX, la fórmula de Manning fue escrita como sigue: V = (1/n) R2/3 S1/2 A través de su "Manual de Hidráulica," King (1918) popularizó la fórmula de Manning en la forma como se conoce hoy en día, así como la aceptación de que el coeficiente C de Manning fuera el recíproco de la n de Kutter. En los Estados Unidos, al valor n se le conoce como el factor de fricción de Manning, o simplemente como la n de Manning. En Europa, la K de Strickler es la misma que la C de Manning, es decir, el recíproco de n. Cuando se usa K en lugar de n, la ecuación de Manning se conoce como la ecuación de Manning-Strickler, o simplemente Strickler. 5.4  LA RUGOSIDAD DE MANNING [Cálculo del Flujo Uniforme]   [Cálculo del Canal de Avenida]   [Flujo Superficial Uniforme]   [Preguntas]   [Problemas]   [Bibliografía]   •   [Arriba]   [Flujo Uniforme]   [Fórmula de Chézy ]   [Fórmula de Manning]   Dada la Ecuación 5-17 (o 5-18), una vez que son conocidas tres variables, la cuarta se puede calcular. Comúnmente, R y S son conocidas, y n se estima, a partir de lo cual se puede calcular V. Éste es el método directo, el cual es la forma más común de usar la ecuación de Manning. Cuando se requiere mayor precisión, o cuando el valor de n no puede ser estimado con certeza, se puede calcular n midiendo la velocidad V, el radio hidráulico R, y la pendiente de fondo S. A este procedimiento se le conoce como el método inverso, o el método de calibración. En la práctica, la mayoría de las aplicaciones de la fórmula de Manning usan el método directo. Estimación del n de Manning No existe un procedimiento exacto para estimar el n de Manning. A continuación se presentan algunas recomendaciones probadas en la práctica. Recomendaciones para la estimación del n de Manning Entender los factores que afectan al n de Manning. Consultar una tabla de valores típicos y estimar el valor de acuerdo a experiencia. Consultar varias colecciones de fotos en las cuales el valor de n ha sido documentado con precisión. Familiarizarse con la apariencia de canales típicos para los cuales los valores de n sean conocidos. Chow (1959) ha presentado una colección ilustrada de veinticuatro (24) canales para los cuales los valores de n han sido establecidos. El rango de valores documentados por Chow es de n = 0.012 (un canal revestido con placas de concreto), a n = 0.150 (un río natural en suelo de arcilla arenosa, con taludes y fondo irregulares). Chow documentó valores tan bajos como n = 0.008 (plástico acrílico) y tan altos como n = 0.200 (llanuras de inundación de corrientes naturales, con vegetación densa en el verano) (Tabla 5-4). Estos valores son aplicables al flujo turbulento. Barnes (1967) ha presentado una colección ilustrada a colores de cincuenta (50) cauces fluviales para los cuales el n de Manning ha sido calculado. La colección de Barnes se puede consultar en Características de Rugosidad de Canales Naturales. El valor más bajo de n documentado por Barnes es n = 0.024, en el Río Columbia en Vernita, Washington (Fig. 5-3). El valor más alto es n = 0.075, en Rock Creek cerca de Darby, Montana (Fig. 5-4). Fig. 5-3  El Río Columbia en Vernita, Washington. Fig. 5-4  Rock Creek cerca de Darby, Montana. Arcement y Schneider (1989) han presentado una colección ilustrada a colores de quince (15) llanuras de inundación en el sureste de los Estados Unidos, para las cuales el valor de n ha sido calculado. La colección de Arcement y Schneider se puede consultar en Coeficientes de Rugosidad de Manning para Canales Naturales y Llanuras de Inundación. El valor más pequeño de n documentado por Arcement y Schneider es n = 0.100, correspondiente a Cypress Creek cerca de Downsville, Louisiana (Fig. 5-5). El valor más alto es n = 0.200, correspondiente a Thompson Creek cerca de Clara, Mississippi (Fig. 5.6). Fig. 5-5  Cypress Creek cerca de Downsville, Louisiana. Fig. 5-5  Thompson Creek cerca de Clara, Mississippi. Los factores que afectan al n de Manning En la práctica, el valor de n de Manning varía en forma considerable. En los cauces naturales puede variar desde aproximadamente 0.020 para algunos ríos muy grandes con fondo liso (Fig. 5-7), hasta un poco más de 0.200 en ríos pequeños de montaña, con pendiente pronunciada (Fig. 5-8). La Tabla 5-3 describe los diversos factores que afectan el coeficiente de rugosidad de Manning. Fig. 5-7  Río Paraguay en Fuerte Coimbra, Mato Groso del Sur, Brasil. Fig. 5-8  Arroyo Rachichuela, en la cuenca del río La Leche, Lambayeque, Perú. Tabla 5-3  Factores que afectan el coeficiente de rugosidad de Manning. Factor Descripción Rugosidad de la superficie Las granulometrías finas conducen a valores bajos, mientras que las granulometrías gruesas conducen a valores altos. Vegetación El tipo, altura, densidad, y distribución espacial de la vegetación afectan a la velocidad de flujo. Los valores de n en canales con vegetación pueden exceder 0.25, y en algunos casos, llegar a más de 0.4. Irregularidades del canal Los bancos de arena, promontorios, depresiones, huecos, y montículos en el lecho del canal crean rugosidad adicional, lo cual produce pérdidas de energía local. Alineación del Canal Generalmente, un canal recto contará con un valor menor de n, mientras que un canal sinuoso tendrá un valor mayor de n. La sinuosidad puede aumentar la rugosidad del canal hasta en un 30% (Chow, 1959). Agradación (depósito) y degradación (erosión) Los cambios en la morfología del canal aumentarán o disminuirán la rugosidad en forma impredecible. El efecto dependerá del tipo de material que forma el lecho, la relación de ancho a profundidad, y la cantidad de sedimentos transportados (carga de sedimentos). Obstrucciones del canal Las obstrucciones de troncos, pilares de puentes y otros obstáculos tienden a aumentar la rugosidad del canal. El efecto dependerá del tipo de obstrucción, tamaño, forma, cantidad y distribución espacial. Tamaño y forma del canal Generalmente, los canales pequeños tienen mayor rugosidad, mientras que los canales grandes tienen menor rugosidad (compare la Fig. 5.7 con la Fig. 5.8). La relación de aspecto es usualmente mayor en canales grandes, lo cual disminuye la rugosidad. Descarga y nivel de la superficie del agua La rugosidad varía en forma impredecible con la descarga y el nivel de la superficie del agua. Las velocidades medias varían de niveles muy bajos a niveles muy altos, y la variación es generalmente compleja. La Figura 2-15 muestra un esquema típico. Estación del año En canales con vegetación, la rugosidad de la superficie aumenta durante el período de crecimiento y disminuye durante el estado latente. Este comportamiento está sujeto a un efecto latitudinal. Carga suspendida y carga de fondo El transporte de sedimentos, ya sea como carga suspendida o carga de fondo, consumirá una energía adicional, lo cual lleva a un incremento de rugosidad. Cowan (1956) ha desarrollado una metodología racional para la estimación del n de Manning. La ecuación de Cowan es:   n  =  (no  +  n1  +  n2  +  n3  +  n4 ) m5   (5-28) en la cual: no = valor básico de n para un canal recto, uniforme, y liso n1 = valor adicional para tomar en cuenta irregularidades de la superficie n2 = valor adicional para tomar en cuenta las variaciones en tamaño y forma de la sección transversal n3 = valor adicional para tomar en cuenta las obstrucciones n4 = valor adicional para tomar en cuenta el efecto de la vegetación m5 = factor para tomar en cuenta la sinuosidad del canal (meandros). La Tabla 5-3 muestra los valores que deben ser usados en la Ec. 5-28. Tabla 5-3  Correcciones para el n de Manning (Ec. 5-28). Condiciones del Flujo Valores Tipo de material en el perímetro del canal Tierra Arena, limo, y arcilla no 0.020 Corte en roca Afloramiento de la roca 0.025 Grava fina Grava de hasta 8 mm de diámetro 0.024 Grava gruesa Grava de más de 8 mm de diámetro 0.028 Grado de irregularidades en la superficie Liso La mejor condición n1 0.000 Menor Canales bien dragados, con pendientes laterales poco erosionadas 0.005 Moderado Canales poco dragados, con pendientes laterales moderadamente erosionadas 0.010 Severo Canales muy erosionados, con superficies altamente irregulares y/o excavadas en roca gruesa 0.020 Variaciones en la forma y tamaño de la sección transversal Gradual Variaciones suaves o pequeñas n2 0.000 Alternadas ocasionalmente Secciones grandes y pequeñas alternadas ocasionalmente, con desplazamiento lateral del flujo 0.005 Alternadas con frecuencia Secciones grandes y pequeñas alternadas con frecuencia, con desplazamiento lateral frecuente del flujo 0.010-0.015 Efecto de obstrucciones Despreciable (a) La medida en que las obstrucciones ocupan o reducen el área de flujo, (b) El carácter de las obstrucciones (bordes agudos u objetos angulares, los cuales inducen una mayor turbulencia que los objetos de superficie curva o lisa), y (c) El lugar y espaciamiento de las obstrucciones, transversales y longitudinales, a lo largo del canal n3 0.000 Menor 0.010-0.015 Apreciable 0.015-0.030 Mayor 0.030-0.060 Efecto de la vegetación Bajo Césped o maleza, en la cual la profundidad de flujo es de 2 a 3 veces la altura de la vegetación n4 0.005-0.010 Medio Césped o maleza, en la cual la profundidad del flujo es de 1 a 2 veces la altura de la vegetación 0.010-0.025 Alto Césped o maleza, en la cual la profundidad de flujo es casi igual a la altura de la vegetación 0.025-0.050 Muy alto Césped o maleza, en la cual la profundidad del flujo es menos de la mitad (1/2) de la altura de la vegetación 0.050-0.100 Sinuosidad del canal Baja Sinuosidad menor a 1.2 m5 1.00 Media Sinuosidad entre 1.2 y 1.5 1.15 Alta Sinuosidad mayor a 1.5 1.30 La Tabla 5-4 muestra los valores de n de Manning para canales de diversos tipos, documentados por Chow (1959). Para cada tipo de canal, se muestran los valores mínimos, normales, y máximos de n. Los valores recomendados para el diseño se indican en negrita. Tabla 5-4  Rango recomendado para los valores de n de Manning. 1 1 2 3 4 Tipo de canal y descripción Mínimo Normal Máximo A Conductos cerrados fluyendo parcialmente llenos   A-1 Metal     a. Latón, liso 0.009 0.010 0.013     b. Acero       1. Soldado 0.010 0.012 0.014       2. Riveteado 0.013 0.016 0.017     c. Hierro fundido       1. Revestido 0.010 0.013 0.014       2. Sin revestir 0.011 0.014 0.016     d. Hierro forjado       1. Negro 0.012 0.014 0.015       2. Galvanizado 0.013 0.016 0.017     e. Metal corrugado       1. Drenaje subsuperficial 0.017 0.019 0.021       2. Drenaje pluvial 0.021 0.024 0.030 1 Chow, V. T. 1959. Open-channel hydraulics. McGraw-Hill, New York.  Haga click -aquí- para mostrar completa la Tabla  5-4. 5.5  CÁLCULO DEL FLUJO UNIFORME [Cálculo del Caudal de Avenida]   [Flujo Superficial Uniforme]   [Preguntas]   [Problemas]   [Bibliografía]   •   [Arriba]   [Flujo Uniforme]   [Fórmula de Chézy]   [Fórmula de Manning]   [Rugosidad de Manning]   La descarga del flujo en canales es (Ec. 5-2): Q  =  V A  =  C R x S y A (5-29) La Ecuación 5-29 se puede expresar como sigue: Q  =  K S y  =  K S 1/2 (5-30) en la cual K = conducción, definida así: K  =  C R x A (5-31) O, alternativamente:             Q K  =  ______           S 1/2 (5-32) De acuerdo a la ecuación de Chézy: K  =  C A R 1/2 (5-33) De acuerdo a la ecuación de Manning, en unidades SI:            1 K  =  ____  A R 2/3            n (5-34) En el sistema de unidades acostumbradas en EE.UU:            1.486 K  =  ________  A R 2/3               n (5-35) La conducción K contiene información sobre la rugosidad y el tamaño y forma de la sección transversal, y no es función de la pendiente del canal. Canales de rugosidad compuesta Por lo general, un canal que desborda sus bancos tiene más de un valor de n de Manning: uno dentro del banco, otro fuera del banco izquierdo, y otro fuera del banco derecho (Fig. 5-9). Se puede calcular un valor compuesto de n bajo el supuesto de que las velocidades en las tres subsecciones son iguales. Esta suposición es conveniente, pero evade la situación de la posible falta de uniformidad del flujo en una sección transversal compuesta. Supóngase un canal de rugosidad variable a lo largo de su perímetro mojado, siendo N el número de subsecciones. Los perímetros húmedos son: P1, P2, P3, ..., PN. Los valores correspondientes de rugosidad son: n1, n2, n3, ..., nN. Asumiendo que todas las velocidades son iguales: V1 = V2 = V3 = VN = V (5-36) Para cualquier subsección i :            1 Vi  =  ____  Ri 2/3 S 1/2            ni (5-37)            1 Vi  =  ____  (Ai / Pi ) 2/3 S 1/2            ni (5-38) El área de flujo para la subsección i es:              Vi 3/2 ni 3/2 Pi Ai  =  _________________                    S 3/4 (5-39) El área total de flujo es:              V 3/2 n 3/2 P A  =  _________________                    S 3/4 (5-40) El área total de flujo es igual a la suma de las subáreas. Por lo tanto: V 3/2 n 3/2 P  =  ∑ (Vi 3/2 ni 3/2 Pi )                                          N (5-41) De acuerdo a la Ec. 5-36, todas las velocidades son iguales. Por lo tanto, la Ec. 5-41 se reduce a:  n 3/2 P  =  ∑ (ni 3/2 Pi )                   N (5-42) Por lo tanto, el valor del n de Manning para un canal de sección transversal compuesta es:                     ∑ (ni 3/2 Pi )                 N n  =  [ _________________ ] 2/3                         P   (5-43) Fig. 5-9 Un canal de sección transversal compuesta. Canales de sección transversal compuesta La sección transversal de un canal puede estar compuesta de varias subsecciones distintas, cada una con rugosidad diferente. Por ejemplo, un canal aluvial sujeto a inundaciones estacionales generalmente consiste de un canal principal y dos canales laterales (Fig. 5-9). Los canales laterales son generalmente más rugosos que el canal principal. Por lo tanto, la velocidad media en el canal principal es generalmente mayor que las de los canales laterales. La ecuación de Manning puede aplicarse por separado a cada subsección, y la descarga total es igual a la suma de las descargas de las subsecciones. Para el canal en su conjunto, la velocidad media es igual a la descarga total dividida entre el área total. El coeficiente de distribución de velocidades aplicable a todo el canal es diferente de los coeficientes de distribución de velocidades aplicables a cada subsección. El coeficiente total de distribución de velocidad puede calcularse como se detalla a continuación. Asumir un número total de subsecciones N y varias subsecciones i, en la cuali varía de 1 a N. De acuerdo a la ecuación de continuidad (Ec. 2-4): Qi  =  Vi Ai (5-44) en la cual Qi = conducción a través de la subsección i, y Vi = velocidad media a través la subsección i, con área de flujo Ai. Además, a partir de la Ec. 5-30: Qi  =  Ki S 1/2 (5-45) en la cual Ki = conducción a través de la subsección i (Ec. 5-33 o 5-34). Combinando las Ec. 5-44 y 5-45: Vi  =  (Ki / Ai) S 1/2 (5-46) La descarga total es:                N Q  =  ∑ Q i         i = 1 (5-47) Usando la Ec. 5-45 en la Ec. 5-47:                      N Q  =  ( Σ Ki )  S 1/2             i = 1 (5-48) De la Ec. 2-4, Q = V A. Por lo tanto:                     N V  =  ( Σ Ki ) ( S 1/2 / A )            i = 1 (5-49) De la Ec. 2-24, el coeficiente de energía se define como sigue:                      N             Σ Vi 3 Ai            i = 1 α  =  _____________                                 V 3A (5-50) Asimismo, a partir de la Ec. 2-31, el coeficiente de cantidad de movimiento es:                      N             Σ Vi 2 Ai            i = 1 β  =  _____________                                 V 2A (5-51) Sustituyendo las Ecs. 5-46 y 5-49 en la Ec. 5-50, el coeficiente compuesto de energía α, aplicable a toda la sección transversal es:                                     N              A2 [  Σ (αi Ki 3 / Ai 2 ) ]                      i = 1 α  =  ____________________________                                                          N                      ( Σ Ki ) 3                         i = 1 (5-52) Asimismo, sustituyendo las Ecs. 5-46 y 5-49 en la Ec. 5-51, el valor compuesto de β es:                                  N              A [  Σ (βi Ki 2 / Ai ) ]                    i = 1 β  =  ____________________________                                                          N                      ( Σ Ki ) 2                         i = 1 (5-53)  Ejemplo 5-1. Calcular los coeficientes α y β correspondientes al caudal pico en una corriente natural (río) de sección compuesta (ver Fig. 5-9), dados los siguientes datos geométricos e hidráulicos (Adaptado de Chow, 1959): (1) (2) (3) (4) (5) (6) (7) i Descripción Área de flujo A (m2) Perímetro mojado P (m) n de Manning α β 1 Canal principal 5360 225 0.035 1.10 1.04 2 Llanura de inundación izquierda 5710 405 0.040 1.11 1.04 3 Llanura de inundación derecha 6030 480 0.045 1.13 1.05 Solución. Para cada componente i, use la Ec. 5-34 para calcular la conducción K, y con ella, los componentes de las Ecs. 5-52 y 5-53. La siguiente tabla muestra un resumen de los resultados. (1) (2) (3) (4) (5) (6) (7) (8) (9) (10) i Descripción Área de flujo A (m2) Perímetro mojado P (m) Radio hidráulico R (m) R 4/3 n K βK 2/A αK 3/A2 1 Canal principal 5,360 225 23.8 8.28 0.035 1,884,282 688,906,605 256,153,800,700 2 Llanura de inundación izquierda 5,710 405 14.10 5.84 0.040 1,238,819 297,772,279 64,725,339,200 3 Llanura de inundación derecha 6,030 480 12.56 5.40 0.045 1,075,270 201,329,329 38,636,378,030 Suma - 17,100 - - - - 4,198,371 1,188,008,213 359,515,517,900 Aplique las sumas de las Columnas 3, 8, y 10 de esta tabla en la Ec. 5-52 para calcular α = 1.42. De igual manera, aplique las sumas de las Columnas 3, 8, y 9 en la Ec. 5-53 para calcular β = 1.15. Cálculo del flujo uniforme De acuerdo a la Fig. 5-10, se deduce la siguiente proporción: z /1 = x /y. Por lo tanto, el ancho de superficie T es:              T  =  b  +  2x  =  b  +  2zy              (5-54) El área de flujo A es: A  =  (b  + x ) y  =  (b  + zy ) y              (5-55) Fig. 5-10  Una sección transversal trapezoidal típica. El perímetro mojado P es:              P  =  b  +  2 (y 2  +  z 2y 2 )1/2   (5-56) Simplificando:              P  =  b  +  2 y ( 1  +  z 2 )1/2              (5-57) A partir de la ecuación de Manning, la descarga Q es:             k Q  =  _____ A  R 2/3 S 1/2             n (5-58) en la cual k = 1 en unidades SI, y k = 1.486 en el sistema acostumbradas en EE.UU. Dado que R = A /P, la Ec. 5-58 se reduce a:     Q n             A 5/3 _________  =  _______     k S 1/2          P 2/3 (5-59) Sustituyendo las Ecs. 5-55 y 5-57 en la Ec. 5-59:     Q n                      [ (b  + zy ) y ] 5/3 _________  =  _____________________________     k S 1/2           [b  +  2 y ( 1  +  z 2 )1/2] 2/3 (5-60) Simplificando:                                       Q n                                                      Q n     [ (b  + zy ) y ] 5/2  -  ( ________ ) 3/2 [ 2 y ( 1  +  z 2 )1/2 ]  -  ( ________ ) 3/2 b  =  0                                    k S 1/2                                                   k S 1/2 (5-61) La profundidad normal y se calcula por medio de la Ec. 5-61 con los siguientes datos de entrada: (1) descarga Q, (2) ancho de fondo b, (3) talud z [z: H a 1: V, Fig. 5-10], (4) pendiente de fondo S, y (5) n de Manning. Por lo tanto, con las Ecs. 5-54 y 5-55:             A D  =  ______             T (5-62)             Q V  =  ______             A (5-63) La Ecuación 5-61 es la fórmula general para el flujo uniforme o normal, aplicable a los canales prismáticos de sección transversal trapezoidal. Para un canal rectangular: z = 0. Asimismo, para un canal triangular de sección simétrica: b = 0. Para resolver la Ec. 5-61 se expresa:                                                   Q n                                                      Q n     f (y)  =  [ (b  + zy ) y ] 5/2  -  ( ________ ) 3/2 [ 2 y ( 1  +  z 2 )1/2 ]  -  ( ________ ) 3/2 b                                                 k S 1/2                                                  k S 1/2 (5-64) Sustituyendo la variable x = y:                                                   Q n                                                      Q n     f (x)  =  [ (b  + zx ) x ] 5/2  -  ( ________ ) 3/2 [ 2 x ( 1  +  z 2 )1/2 ]  -  ( ________ ) 3/2 b                                                 k S 1/2                                                  k S 1/2 (5-65) La Ecuación 5-65 se resuelve mediante un procedimiento de prueba y error. A continuación se describe un algoritmo basado en el valor de la función. La Figura 5-11 es un ejemplo del flujo normal. Algoritmo para calcular la profundidad normal basado en el valor de la función Asumir un valor inicial de x = 0. Así: f (0) = - [ ( Qn ) / (k So ) ] 3/2 b lo cual resulta en un número negativo grande. Por lo tanto se confirma que el valor inicial de la función es menor que cero. Asumir un valor inicial apropiado del intervalo Δx, por ejemplo: Δx = 1. Establecer x = x + Δx Calcular f (x) Terminar cuando Δx < Δx TOL. Un valor típico de Δx TOL es 0.0001. Si f (x) < 0, regresar al Paso 3. Si f (x) > 0, establecer:  Δx = 0.1 Δx Establecer:  x = x - 9 Δx Regresar al Paso 4. Fig. 5-11  Flujo normal en el canal Wellton-Mohawk, Wellton, Arizona. La aproximación de Newton El algoritmo anteriormente descrito utiliza únicamente el valor de la función para la aproximación a la raíz. Como alternativa, la aproximación de Newton basada en la tangente converge más rápidamente. Sin embargo, para que el método de Newton funcione correctamente, es necesario primero acercarse a la raíz usando la iteración basada en el valor de la función. De lo contrario, es probable que el método de Newton no converja. Con referencia a la Fig. 5-12, el valor de la tangente en xo es:                    f(xo) f '(xo)  =  _________                  xo - xr (5-66) en la cual xo = valor de prueba de x, f(xo) = valor de la función en xo, xr = nuevo valor de x, que se aproxima más a la raíz. Fig. 5-12  Esquema ilustrativo de la iteración de Newton. Resolviendo para xr  de la Ec. 5-66:                     f (xo) xr  =  xo  -  ________                     f '(xo) (5-67) Como se muestra en la Fig. 5-12, cuando f (xo) aumenta con xo (como es el caso de la Ec. 5-65) al pasar por la raíz, el valor de la función y el valor de la tangente son positivos; por lo tanto, el denominador de la Ec. 5-66 también es positivo, y xr  se encuentra a la izquierda de xo. Con cada iteración subsecuente, la raíz es aproximada en forma de zig-zag, hasta que la tolerancia especificada sea satisfecha. Se puede mostrar fácilmente que la Ec. 5-67 también se aplica cuando f (xo) disminuye a medida que xo aumenta, como en el caso del flujo crítico, ver la Sección 4.2. El valor de f '(x) es:                                                                                                         Q n    f ' (x)  =  x 5/2 (5/2) (b  + zx ) 3/2 z  +  (b  + zx ) 5/2 (5/2) x 3/2  -  ( ________ ) 3/2 [ 2 ( 1  +  z 2 )1/2 ]                                                                                                        k S 1/2 (5-68) Simplificando la Ec. 5-68:                                                                                               Q n    f ' (x)  =  5 z x 5/2 (b  + zx ) 3/2  +  5 x 3/2 (b  + zx ) 5/2  -  ( ________ ) 3/2 [ 4 ( 1  +  z 2 )1/2 ]                                                                                             k S 1/2 (5-69) El procedimiento para la aproximación de Newton se describe a continuación. Algoritmo para calcular la profundidad normal:  Aproximación de Newton Asumir un valor inicial de xo = 0. Asumir un valor inicial apropiado del intervalo Δx, por ejemplo: Δx = 1. Establecer xo = xo + Δx Calcular f (xo) Si f (xo) < 0, regresar al Paso 3. Si f (xo) > 0, calcular la raíz xr  usando las Ecs. 5-65 y 5-69:                     f (xo) xr  =  xo  -  ________                     f '(xo) Terminar cuando la diferencia (xr  -  xo) es demasiado pequeña. Un valor típico de la diferencia es 0.0001. De lo contrario, establecer xo = xr  y regresar al Paso 6.  Ejemplo 5-2. Usando CANAL EN LÍNEA 01, calcular la profundidad normal para las siguientes condiciones: Q = 3 m3/s, b = 5 m, z = 1, S = 0.001, n = 0.015.   CÁLCULO EN LÍNEA. Usando CANAL EN LÍNEA 01, la profundidad normal es yn = 0.473 m, la velocidad normal es vn = 1.16 m/s, y el número de Froude normal es Fn = 0.562.  Ejemplo 5-3. Usando CANAL EN LÍNEA 01, calcular la profundidad normal para las siguientes condiciones: Q = 20 pies3/s, b = 13 pies, z = 2, S = 0.0008, n = 0.013.   CÁLCULO EN LÍNEA. Usando CANAL EN LÍNEA 01, la profundidad normal es yn = 0.631 pies, la velocidad normal es vn = 2.221 pies/s, y el número de Froude normal es Fn = 0.514. 5.6  CÁLCULO DEL CAUDAL DE AVENIDA [Flujo Superficial Uniforme]   [Preguntas]   [Problemas]   [Bibliografía]   •   [Arriba]   [Flujo Uniforme]   [Fórmula de Chézy ]   [Fórmula de Manning]   [Rugosidad de Manning]   [Cálculo del Flujo Uniforme]   La combinación de velocidades y profundidades altas que predominan durante una inundación aumenta el riesgo de accidentes y lesiones corporales (Fig. 5-13). Por lo tanto, generalmente no es posible medir la descarga durante una inundación. Sin embargo, se puede obtener un estimado de la descarga máxima de forma indirecta, por medio del uso de fórmulas de flujo en canales. Éste es el principio del método de pendiente-área. Fig. 5-13  Una inundación severa en un río tropical. Para aplicar el método de pendiente-área en un tramo de río, se requieren los siguientes datos: La longitud del tramo; La caída, es decir, el cambio de elevación en la superficie del agua a través del tramo; El área de flujo, el perímetro mojado, y los coeficientes de carga de velocidad en las secciones aguas arriba y aguas abajo; y El valor medio del n de Manning para el tramo. La selección de un tramo adecuado está basada en los siguientes criterios: Las marcas que indican la elevación del agua se deben ser fácilmente identificables (Fig. 5-14). El tramo debe ser lo suficientemente largo para que la caída pueda ser medida con precisión. La forma de la sección transversal y las dimensiones del canal deben ser relativamente constantes. El tramo debe ser relativamente recto, aunque es preferible un tramo con contracción a un tramo con expansión. Deben evitarse los puentes, curvas, cascadas y otras características que provoquen una falta de uniformidad del flujo. Fig. 5-14  Un campesino mostrando el nivel de agua alcanzado por las inundaciones, Karnataka, India (1991). La precisión del método de pendiente-área mejora a medida que aumenta la longitud del tramo (Fig. 5-15). Un tramo adecuado debe satisfacer uno o más de los siguientes criterios: La relación de la longitud del tramo a la profundidad hidráulica debe ser mayor que 75; La caída debe ser mayor que o igual a 0.15 m: F ≥ 0.15 m; y La caída debe ser mayor que cualquiera de las cargas de velocidad calculadas en las secciones aguas arriba y aguas abajo. Fig. 5-15  Ilustración del método de pendiente-área. El procedimiento consta de los siguientes pasos: Calcular la conducción K en las secciones aguas arriba y aguas abajo:            1 K1 = ( __ ) A1 R1 2/3            n (5-70a)            1 K2 = ( __ ) A2 R2 2/3            n (5-70b) en las cuales K = conducción, A = área de flujo, R = radio hidráulico, n = valor de Manning promedio en el tramo, y 1 y 2 indican los secciones aguas arriba y aguas abajo, respectivamente (La Ec. 5-70 esta dada en unidades SI). Calcular la conducción equivalente del tramo, igual a la media geométrica de las conducciones aguas arriba y aguas abajo: K = ( K1 K2 )1/2 (5-71) en la cual K = conducción del tramo. Calcular la primera aproximación a la pendiente de energía:            F S  =   ___            L (5-72) en la cual S = primera aproximación a la pendiente de energía, F = caída, y L = longitud del tramo. Calcular la primera aproximación a la descarga pico: Qi  =  K S 1/2 (5-73) en la cual Qi  = primera aproximación a la descarga pico. Calcular las cargas de velocidad:            α1 ( Qi /A1 ) 2 hv1 =  ______________                    2g (5-74a)            α2 ( Qi /A2 ) 2 hv2 =  ______________                    2g (5-74b) en las cuales hv1 y hv2 = cargas de velocidad en las secciones aguas arriba y aguas abajo respectivamente, α1 y α2 = coeficientes de carga de velocidad aguas arriba y aguas abajo respectivamente, y g = aceleración de la gravedad. Calcular un valor actualizado de la pendiente de energía:            F + k ( hv1 - hv2 ) Si =  ___________________                      L (5-75) en la cual Si = valor actualizado de la pendiente de energía, y k = coeficiente de pérdida de carga. Para el caso de expansión de flujo, es decir, A2 > A1, k = 0.5; para la contracción de flujo, es decir, A1 > A2, k = 1. Calcular un valor actualizado de la descarga pico Qi  =  K Si 1/2 (5-76) Regresar al paso 5 y repetir los pasos 5 a 7. El procedimiento termina cuando la diferencia entre dos valores sucesivos de la descarga pico Q obtenida en el paso 7 es despreciable. En la práctica, generalmente se requieren por lo menos dos iteraciones.  Ejemplo 5-4. Calcular la descarga pico (el caudal de avenida) por el método de pendiente-área, dados los siguientes datos: Longitud del tramo L = 1,450 m; caída F = 2.2 m; coeficiente de Manning en el tramo n = 0.035; área de flujo aguas arriba A1 = 432 m2; perímetro mojado aguas arriba P1 = 85 m; coeficiente de carga de velocidad aguas arriba α1 = 1.15; área de flujo aguas abajo A2 = 455 m2; perímetro mojado aguas abajo P2 = 92 m; y coeficiente de carga de velocidad aguas abajo α2 = 1.12. Solución. Use la Ec. 5-70 para calcular las conducciones aguas arriba: K1 = 36,486; y aguas abajo: K2 = 37,736. Use la Ec. 5-71 para calcular la conducción promedio en el tramo: K = 37,106. Use la Ec. 5-72 para calcular la primera aproximación (i = 0) a la pendiente de la línea de energía. Use la Ec. 5-73 para calcular la primera aproximación (i = 0) de la descarga pico. Proceda en forma recursiva. Los resultados se muestran en la tabla siguiente. (1) (2) (3) (4) (5) i hv1 (m) hv2 (m) Pendiente de energía S Descarga pico Q (m3/s) 0 - - 0.001517 1,445.3 1 0.656 0.576 0.001545 1,458.4 2 0.668 0.586 0.001546 1,458.7 La Columna 5 muestra la descarga pico calculada después de cada paso por la Ec. 5-73. El resultado final se muestra en la Fila i = 2, Columna 5 (segunda y última iteración).   CÁLCULO EN LÍNEA. Usando el calculador PENDIENTE-AREA, la descarga pico es: Q = 1,458.6 m3/s. 5.7  FLUJO SUPERFICIAL UNIFORME [Preguntas]   [Problemas]   [Bibliografía]   •   [Arriba]   [Flujo Uniforme]   [Fórmula de Chézy]   [[Fórmula de Manning]   [Rugosidad de Manning]   [Cálculo del Flujo Uniforme]   [Cálculo del Caudal de Avenida]   El flujo sobre una cuenca o superficie de captación se conoce como flujo superficial. Generalmente en el flujo superficial, la profundidad de flujo es muy pequeña comparada con el ancho. Bajo estas condiciones, el flujo puede ser laminar o turbulento, dependiendo de las rugosidades absoluta y relativa. Si las velocidades y profundidades de flujo son suficientemente pequeñas, el flujo puede ser laminar; de lo contrario, el flujo puede ser transitorio o turbulento, dependiendo del número de Reynolds (Sección 1.4). En el flujo superficial, suele presentarse un régimen mixto laminar-turbulento. Este tipo de flujo se caracteriza por cambios de régimen bajo las condiciones de flujo espacialmente variado que normalmente se encuentran en el flujo de superficie. Cuando el flujo es laminar, el exponente de la curva de gasto es β = 3. En condiciones de flujo turbulento, el exponente es β = 3/2 para la fricción de Chézy, y β = 5/3 para la fricción de Manning. Los regímenes mixtos de flujo laminar-turbulento presentan valores de β que varían entre laminar y turbulento. Flujo superficial laminar Con referencia a la Fig. 5-16, el esfuerzo de corte que actúa en el nivel P es: τa  =  γ (ym - y ) S (5-77) De acuerdo a la ley de viscosidad de Newton, el esfuerzo de corte en P es proporcional al gradiente vertical de velocidades:                 dv τr  =  μ   _____                 dy (5-78) en la cual μ = una constante de proporcionalidad conocida como viscosidad dinámica. Fig. 5-16  Ilustración del flujo superficial uniforme. Igualando los esfuerzos actuantes y resistentes:        dv μ   _____  =  γ (ym - y ) S        dy (5-79) En forma diferencial:         μ dv    =  γ (ym - y ) S dy         (5-80) El peso específico  γ = ρg, y la viscosidad dinámica μ = ρν, en la cual ν = viscosidad cinemática. Por lo tanto, la Ec. 5-80 se reduce a:             gS dv  =  _____  (ym - y ) dy              ν (5-81) Integrando la Ec. 5-81:             gS v  =  ∫ _____  (ym - y ) dy              ν (5-82)           gS                    y 2 v  =  _____  [ ym y  -  _____ ]  +  C            ν                      2 (5-83) en la cual C es una constante de integración. Para v = 0, y = 0, por lo tanto: C = 0, y la relación velocidad media/ profundidad de flujo es:           gS                    y 2 v  =  _____  [ ym y  -  _____ ]            ν                      2 (5-84) La Ecuación 5-84 indica que el perfil de velocidades en el flujo superficial uniforme tiene una distribución parabólica. La curva de descarga-profundidad se obtiene integrando la Ec. 5-84 entre los límites de 0 y ym, es decir, desde el fondo hasta la superficie, lo cual resulta en:                          gS                         y 2 q  =  ∫ v dy  =  _____   ∫  [ ym y  -  _____ ] dy                           ν                           2 (5-85)           gS       ym 2       ym 3 q  =  _____  [ _____  -  _____ ]            ν          2            6 (5-86) la cual se reduce a:            gS               q  =  ______   ym 3            3ν           (5-87)                           q  =  CL ym 3                       (5-88) en la cual CL = coeficiente de la curva de gasto correspondiente al flujo laminar, definido como sigue:              gS               CL  =  ______                3ν           (5-89) En el flujo laminar, el exponente de la curva de gasto es β = 3 (Ec. 5-88), y la curva de gasto es una función de la fricción interna, o viscosidad interna, representada por la viscosidad cinemática ν. Por lo tanto, el flujo laminar es una función de la temperatura. Dada la Ec. 5-77, la velocidad media en el flujo laminar, v = q /ym, es:                           v  =  CL ym 2                       (5-90) La curva de gasto correspondiente al flujo turbulento de acuerdo a Chézy es:                             q  =  C S 1/2 ym 3/2                       (5-91) La curva de gasto correspondiente al flujo turbulento de acuerdo a Manning, en unidades SI es:                             q  =  (1/n) S 1/2 ym 5/3                       (5-92) De igual manera, en las unidades acostumbradas en EE.UU.:                             q  =  (1.486 / n) S 1/2 ym 5/3                       (5-93) Puede observarse que el exponente de la curva de gasto varía de β = 3 para el flujo laminar (Ec. 5-88), a β = 3/2 para la fricción turbulenta de acuerdo a Chézy (Ec. 5-91), o a β = 5/3 para la fricción turbulenta de acuerdo a Manning (Ec. 5-92). En el flujo superficial uniforme, los valores de β comprendidos en el intervalo entre flujo laminar y turbulento representan la condición mixta laminar-turbulenta (Sección 1.3). El número de Vedernikov es:              (β - 1) v               V  =  ____________                (g y )1/2           (5-94) En el caso de V  =  1, el flujo es neutralmente estable, lo cual lleva al desarrollo de ondas de rollo (Fig. 5-17). A continuación se describe la relación entre el exponente β y el número de Vedernikov V . La relación entre el exponente β y el número de Vedernikov En el flujo laminar:  β = 3. Por lo tanto, V  =  2 F, y el flujo llega a ser inestable cuando F  =  0.5. En el flujo turbulento con la fricción de Chézy:  β = 1.5. Por lo tanto, V  =  0.5 F, y el flujo llega a ser inestable cuando F  =  2. En el flujo turbulento con la fricción de Manning:  β = 5/3. Por lo tanto, V  =  (2/3) F, y el flujo llega a ser inestable cuando F  =  3/2, es decir, F  =  1.5.  Ejemplo 5-5. Calcular el número de Vedernikov, dados velocidad de flujo v = 1 m/s; profundidad de flujo y = 1 m; y exponente de la curva de gasto β = 1.63. Solución. Usando la Ec. 5-94: V = (β - 1) v / (gy) = 0.201.   CÁLCULO EN LÍNEA. Usando el calculador ENLINEA VEDERNIKOV, el número de Vedernikov es: V = 0.201. Fig. 5-17  Ondas de rollo en el aliviadero del Reservorio Turner, Condado de San Diego, California. PREGUNTAS [Problemas]   [Bibliografía]   •   [Arriba]   [Flujo Uniforme]   [Fórmula de Chézy]   [Fórmula de Manning]   [Rugosidad de Manning]   [Cálculo del Flujo Uniforme]   [Cálculo del Caudal de Avenida]   [Flujo Superficial Uniforme]   ¿Cuándo se convierte en inestable el flujo uniforme? ¿En qué se basa la fórmula de Chézy? ¿Cuál es la diferencia entre las fórmulas de Manning y Chézy? ¿Cuál es el valor mínimo del n de Manning que se puede lograr en la práctica? ¿Cuál es el rango de valores del n de Manning medido por Barnes? ¿Cuál es el rango de valores del n de Manning medido por Arcement y Schneider en llanuras de inundación? ¿Por qué el cálculo de rugosidad compuesta utilizando la Ec. 5-43 es sólo una aproximación? ¿Cuáles son las cinco variables que se utilizan en el cálculo del flujo uniforme en un canal trapezoidal? ¿Por qué es mejor usar la aproximación de Newton para calcular la profundidad normal, en lugar de usar únicamente la aproximación que usa el valor de la función? ¿Cuál es la relación mínima de la longitud del tramo a la profundidad hidráulica en el método de pendiente-área? ¿Cuál es el exponente de la curva de gasto en condiciones de flujo laminar? ¿Cuál es el exponente de la curva de gasto en condiciones de flujo turbulento de acuerdo a Chézy en un canal hidráulicamente ancho? ¿Cuál es el exponente de la curva de gasto en condiciones de flujo turbulento de acuerdo a Manning en un canal hidráulicamente ancho? ¿Bajo qué valor del número de Froude es probable que el flujo llegue a ser inestable en condiciones de flujo laminar? PROBLEMAS [Bibliografía]   •   [Arriba]   [Flujo Uniforme]   [Fórmula de Chézy]   [Fórmula de Manning]   [Rugosidad de Manning]   [Cálculo del Flujo Uniforme]   [Cálculo del Caudal de Avenida]   [Flujo Superficial Uniforme]   [Preguntas]   Demostrar que el factor de fricción de Darcy-Weisbach y la n de Manning se relacionan por medio de la siguiente ecuación: fD = 8 g n 2 / (k 2 R 1/3) en la cual el factor de fricción fD = Darcy-Weisbach, g = aceleración de la gravedad, R = radio hidráulico, y k = una constante para el sistema de unidades, igual a 1 en unidades SI y 1.486 en unidades acostumbradas en EE.UU. Expresar la relación en unidades SI y en unidades acostumbradas en EE.UU . Calcular la descarga Q usando la ecuación de Manning, con los siguientes datos: área de flujo A = 23.5 pies2, radio hidráulico R = 5.6 pies, pendiente de fondo S = 0.0025, y n de Manning = 0.035. Calcular la descarga Q usando la ecuación de Manning, con los siguientes datos: área del flujo A = 45 m2; radio hidráulico R = 6 m; pendiente de fondo S = 0.003; n de Manning = 0.04. Dado f = 0.0025, calcular la descarga Q para un área de flujo A = 12.4 m2, radio hidráulico R = 2.1 m; y pendiente de fondo S = 0.0015. Dado f = 0.0035, calcular la descarga Q para un área de flujo A = 18 pies2, radio hidráulico R = 4.5 pies, y la pendiente de fondo S = 0.0018. Usar CANAL EN LÍNEA 01 para calcular la profundidad normal, la velocidad y el número de Froude con los siguientes datos: Q = 150 m3/s, b = 10 m, z = 2, So = 0.0005, n = 0.025. Usar el CANAL EN LÍNEA 01 para calcular la profundidad normal, la velocidad y número de Froude con los siguientes datos: Q = 250 pies cúbicos por segundo, b = 20 pies, z = 1, So = 0.001, n = 0.030. Usando CANAL EN LÍNEA 15, calcular la descarga para un canal prismático, con: b = 20 pies, y = 3 pies, z = 2, n = 0.025, S = 0.0016. Fig. 5-18  Esquema para un canal trapezoidal. Usando el CANAL EN LÍNEA 15, calcular la descarga para un canal prismático con b = 6 m, y = 1 m, z = 1.5, n = 0.015, S = 0.0002. Una inundación en Clearwater Creek ha dejado marcas de agua observables en un determinado tramo del río. Los datos hidráulicos se han medido en dos secciones transversales A y B, a una distancia de 1,850 pies. La caída del tramo entre las secciones transversales es de 9.1 pies y el n de Manning promedio es 0.035. Aguas arriba, el área de flujo, perímetro mojado, y coeficiente α de Coriolis son 550 pies2, 55 pies, y 1.17, respectivamente. Aguas abajo, el área de flujo, perímetro mojado, y coeficiente α son 620 pies2, 52 pies, y 1.10, respectivamente. Usar PENDIENTE-ÁREA para calcular la descarga de avenida. Calcular la descarga por unidad de ancho en un plano de flujo superficial, bajo flujo laminar y con una profundidad media de 1.5 cm y pendiente de fondo de 0.001. Asumir que la temperatura del agua es T = 20oC. Expresar la descarga en L/s/m. Usando el cálculo diferencial, calcular la profundidad relativa y/D correspondiente a la descarga máxima en una alcantarilla circular de diámetro D. Explicar la razón por la cual la descarga máxima no occurre para y = D. BIBLIOGRAFÍA •   [Arriba]   [Flujo Uniforme]   [Fórmula de Chézy]   [Fórmula de Manning]   [Rugosidad de Manning]   [Cálculo del Flujo Uniforme]   [Cálculo del Caudal de Avenida]   [Flujo Superficial Uniforme]   [Preguntas]   [Problemas]   Arcement, G. J. y V. R. Schneider. 1989. Guide for selecting Manning's roughness coefficients for natural channels and flood plains. U.S. Geological Survey Water-Supply Paper 2339, Washington, D.C. Barnes, H. A. 1967. Roughness characteristics of natural channels. U.S. Geological Survey Water-Supply Paper 1849, Washington, D.C. Cowan, W. L. 1956. Estimating hydraulic roughness coefficients. Agricultural Engineering, Vol. 37, No. 7, pp. 473-475, Julio. Chow, V. T. 1959. Open-channel hydraulics. McGraw-Hill, Nueva York. Henderson, M. H. 1966. Open-channel flow. Macmillan, Nueva York. Simons, D. B., y E. V. Richardson. 1966. Resistance to flow in alluvial channels. U.S. Geological Survey Professional Paper 422-J, Washington, D.C. Williamson, J. 1951. The laws of flow in rough pipes. La Houille Blanche, Vol. 6, No. 5, Septiembre-Octubre, p. 738. http://hidraulicadecanales.sdsu.edu 201102 06:35