Distribución lateral de velocidades en cauces naturales

INGENIERÍA DEL AGUA · VOL. 12 · Nº 3 SEPTIEMBRE 2005 1 1 Laboratorio de Hidráulica, Depto. de Ing. Civil, Univ. Tecnológica Nacional. – Fac. Reg. Córdoba. E-mail: [email protected] 2 INA (Instituto Nacional del Agua) y Fac. de Ingeniería - Universidad de Buenos Aires. E-mail: [email protected] 3 CRSA - INA (Instituto Nacional del Agua). E-mail: [email protected].edu Artículo recibido el 3 de octubre de 2003, recibido en forma revisada el 26 de noviembre de 2004 y aceptado para su publicación el 14 de diciembre de 2004 . Pueden ser remitidas discusiones sobre el artículo hasta seis meses después de la publicación del mismo siguiendo lo indicado en las “Instrucciones para autores”. En el caso de ser aceptadas, éstas serán publicadas conjuntamente con la respuesta de los autores. INTRODUCCIÓN Los estudios de ingeniería hidráulica relacionados con obras fluviales requieren usualmente conocer la distribución lateral de velocidades en cauces naturales, de geometría irregular y en varias ocasiones constituidos por un cauce principal y planicies de inundación. Estas velocidades, que se asumen promediadas en la vertical, responden en términos generales al hecho de ser mayores donde la profundidad de flujo también lo es, como se representa en la Figura 1. Los modelos numéricos hidrodinámicos bidimensionales (como por ejemplo el programa RMA2, del U.S. Army Corps of Engineers) son considerados como modelos teóricos apropiados para resolver este problema. De todos modos, como en la mayoría de las aplicaciones la escala longitudinal de variación del flujo es mucho mayor a la escala lateral (dada por el ancho de flujo), las variaciones longitudinales y laterales de flujo pueden ser calculadas independientemente (Menéndez 2003). Las primeras, a través de los modelos hidrodinámicos unidimensionales basados en las ecuaciones de Saint – Vénant; en tanto que las segundas, a través de los modelos 1D de distribución lateral. Figura 1. Perfil transversal de velocidad típico de una sección compuesta Resumen: En este trabajo se presentan los resultados de la aplicación de dos modelos para la predicción de la distribución lateral de velocidades en cinco casos de cauces naturales. Los modelos considerados son el Método del Canal Dividido (DCM) a través de su implementación en el programa HEC-RAS, y el Método de la Distribución Lateral (LDM) a través de la solución analítica propuesta por Shiono y Knight. Los casos analizados corresponden a cinco secciones ubicadas sobre diversos ríos: el río Severn en Montford (UK), y las estaciones de medición Eldorado y Puerto Libertad sobre el Alto río Paraná, el río de La Suela, en la provincia de Córdoba, y el río Colorado en Buta Ranquil (Neuquén), todas estas en la República Argentina. Se presentan los valores ajustados de los parámetros de los modelos, haciendo hincapié en la interpretación física de los mismos. Se compara el desempeño de los modelos aplicados, y se analizan sus ventajas y desventajas relativas. Se concluye que, a pesar de que el modelo SKM es matemáticamente más complejo, la ganancia obtenida para estimar la distribución lateral de velocidades en cauces naturales es lo suficientemente significativa como para recomendar su utilización. Palabras clave: Hidrodinámica, distribución de velocidades, cauces naturales. DISTRIBUCION LATERAL DE VELOCIDADES EN CAUCES NATURALES Juan F. Weber1, Ángel N. Menéndez2, Leticia Tarrab3 INGENIERÍA DEL AGUA · VOL. 12 · Nº 3 SEPTIEMBRE 2005 2 Juan F. Weber, Ángel N. Menéndez y Leticia Tarrab INGENIERÍA DEL AGUA · VOL. 12 · Nº 3 SEPTIEMBRE 2005 3 DISTRIBUCION LATERAL DE VELOCIDADES EN CAUCES NATURALES Uno de los modelos más difundidos para la predicción de la distribución lateral de velocidades es el denominado Método del Canal Dividido (DCM por sus siglas en inglés) basado en la hipótesis del método de Lotter (1933) para el cálculo del coeficiente de Manning equivalente. Este método, de sencilla aplicación, no contempla la redistribución lateral de velocidades impuesta por la difusión turbulenta. Para considerarla, es necesario partir de la ecuación general de movimiento (ecuación de Navier-Stokes) y llegar al denominado Método de la Distribución Lateral (LDM, por sus siglas en inglés) que presenta la velocidad promediada en la vertical como una función de la coordenada lateral, a través de una ecuación diferencial ordinaria no lineal de segundo orden, que es posible resolver analíticamente asumiendo una función sencilla para representar la geometría del lecho. La solución analítica del LDM da origen al denominado Método de Shiono y Knight (SKM, por sus siglas en inglés) que es el segundo modelo considerado en este trabajo. MODELOS CONSIDERADOS Método del canal dividido (DCM) Lotter (1933) y posteriormente Einstein y Banks (1950) trabajaron sobre la hipótesis que el caudal total en la sección puede obtenerse como la suma de los caudales parciales fluyentes por cada subsección, asumiendo que la pendiente de la línea de energía es constante para todas las subsecciones y que los perímetros mojados parciales correspondientes a cada subárea son los dados por el contorno fijo únicamente. Esto implica que no existe difusión lateral de cantidad de movimiento; las pérdidas de energía son sólo debidas a la fricción. El caudal parcial Qi de cada subsección se puede obtener como: (1) donde: S0 pendiente longitudinal Ki conductividad hidráulica de cada subsección, dada por (2) en esta ecuación: Ai Area de la subsección i Ri Radio hidráulico de la subsección i ni Coeficiente de rugosidad de Manning de la subsección i La distribución lateral de velocidades puede ser estimada a partir de los caudales Qi calculados con la ecuación (1). La velocidad media en una subsección, Vi, se puede estimar como: (3) Este modelo está implementado en el programa HEC-RAS en su comando Flow Distribution. Modelo HEC-RAS Este modelo unidimensional tiene como objetivo la resolución de la ecuación de Bernoulli, para el flujo permanente gradualmente variado: (4) (5) (6) Donde: Yi Tirante o profundidad de flujo en la sección i Zi Cota o elevación del fondo en la sección i Vi Velocidad media en la sección i ai Coeficiente de distribución de velocidades o de Coriolis correspondiente a la sección i he Pérdidas de energía en el tramo entre las secciones 1 y 2 L Longitud del tramo entre las secciones 1 y 2 Sfi Gradiente longitudinal de pérdidas por fricción en la sección i C Coeficiente por pérdidas por contracción – expansión. El programa utilizado es HEC – RAS (River Analysis System) del U.S. Army Corps of Engineers en su versión 3.0 (HEC, 2001). Para el desarrollo del presente trabajo resulta de interés la forma en que el modelo unidimensional estima una distribución lateral de velocidades. El usuario debe indicar al programa que desea obtener el perfil transversal de velocidades mediante la opción Flow Distribution. En ella, debe indicar el número de puntos en los que desea obtener los valores de la velocidad promediada en el área de influencia de cada uno de ellos. INGENIERÍA DEL AGUA · VOL. 12 · Nº 3 SEPTIEMBRE 2005 2 Juan F. Weber, Ángel N. Menéndez y Leticia Tarrab INGENIERÍA DEL AGUA · VOL. 12 · Nº 3 SEPTIEMBRE 2005 3 DISTRIBUCION LATERAL DE VELOCIDADES EN CAUCES NATURALES Suponiendo conocidos los valores de Sf y el nivel del agua en una sección determinada (los cuales fueron calculados al resolver el perfil de flujo en el problema dado), los pasos que sigue el programa son los siguientes: 1. Divide el área de flujo en la cantidad de partes definida por el usuario (Figura 2), y luego calcula el área, perímetro mojado, radio hidráulico y profundidad hidráulica para cada parte. Figura 2. Subdivisión de la sección transversal en HECRAS 2. Con los valores calculados previamente, y el valor del n de Manning correspondiente, calcula la conductividad y el porcentaje de flujo en cada parte. 3. Suma las conductividades de cada parte. En general, esta suma será distinta a la conductividad de la sección completa. Para corregir esta diferencia, calcula la relación entre la conductividad de la sección completa y la suma de las conductividades parciales. Luego multiplica la conductividad de cada sector por esa relación. Con esa conductividad determina el caudal en cada sector. 4. Por último, calcula la velocidad media en cada sector como el cociente del caudal parcial antes determinado y el área correspondiente. Método de la Distribución Lateral (LDM) Una de las principales limitaciones del modelo DCM es la imposibilidad que tiene de considerar las pérdidas energéticas debido a la difusión lateral de cantidad de movimiento. Para poder incluirlas, es necesario plantear un modelo basado en las ecuaciones generales de movimiento o ecuaciones de Navier – Stokes, con las simplificaciones correspondientes. En este caso, es posible obtener una ecuación diferencial que permita predecir la distribución lateral de velocidades en la sección, asumiendo las condiciones de flujo permanente y uniforme y una cierta profundidad de flujo dada. Esta ecuación da lugar al denominado Lateral Distribution Method (LDM). La solución de la ecuación diferencial del LDM puede hacerse en forma numérica o analítica. El primer procedimiento fue propuesto por primera vez por Wark et al.(1990). En el segundo caso, el desarrollo original es debido a Shiono y Knight (1988, 1991). La integración de la ecuación de cantidad de movimiento en la dirección del flujo (ecuación de Navier – Stokes) en combinación con la ecuación de continuidad, para un flujo permanente y uniforme asumiendo como modelo de turbulencia el de la viscosidad de torbellino o de Boussinesq, da lugar a la ecuación diferencial del denominado método de la distribución lateral (LDM): (7) en la cual: r: Densidad del fluido g: Aceleración de la gravedad, S0: Pendiente longitudinal del lecho Y: Profundidad de flujo V: Velocidad del Flujo y: Coordenada horizontal medida transversal a la corriente s: talud lateral f: Coeficiente de fricción de Darcy l: Viscosidad adimensional de torbellino de Boussinesq. (rUV)d Tensión de Reynolds debido a las corrientes secundarias Shiono y Knight (1988, 1991) obtienen las siguientes soluciones analíticas de la ecuación (7) : Para fondo plano (horizontal) (s à •): (8) Para fondo lineal (inclinado) (0 < s< •): (9) En estas ecuaciones: (10) . (11). (12) INGENIERÍA DEL AGUA · VOL. 12 · Nº 3 SEPTIEMBRE 2005 4 Juan F. Weber, Ángel N. Menéndez y Leticia Tarrab INGENIERÍA DEL AGUA · VOL. 12 · Nº 3 SEPTIEMBRE 2005 5 DISTRIBUCION LATERAL DE VELOCIDADES EN CAUCES NATURALES (13) (14) (15) En las ecuaciones (8) y (9), los coeficientes ai y bi se obtienen al plantear las condiciones de contorno que corresponden a la de no deslizamiento en los contornos rígidos y las de continuidad y derivabilidad de la función velocidad en cada cambio de pendiente transversal del fondo (Ecuaciones (16) y (17)), llegando de ese modo a un sistema lineal de ecuaciones que permite calcular los valores de estos coeficientes (Weber y Menéndez, 2003). (16) (17) A pesar de que Abril (2003) ha demostrado que la ecuación (17) debería reemplazarse por la condición de continuidad de la tensión de corte, su utilización no produce errores notables. El parámetro G de la Ecuación (15) representa el gradiente lateral de la fuerza unitaria aparente debida de la tensión de Reynolds originada por la corrientes secundarias que se desarrollan al pie del talud del canal principal. Este parámetro se determina experimentalmente. Figura 3. Subdivisión típica de un canal de sección compuesta. A partir de la ecuación (7), o de las ecuaciones (8) y (9), es posible ver que para un tramo determinado de la sección transversal, el modelo analítico a considerado tiene dos parámetros, uno que pondera la fricción con el contorno fijo (parámetro f) y otro que considera la difusión lateral de momentum (parámetro l). Estrictamente hablando, el modelo tiene 2n parámetros, siendo n el número de tramos rectos en que es posible discretizar la sección transversal. De hecho, se deberán definir, por cada tramo recto, un par de valores l y f. CASOS DE ESTUDIO Se presentan a continuación los casos considerados en el presente trabajo. Todos se corresponden con mediciones de campo. El primer conjunto de datos proviene del río Severn, en Montford (UK) cuya aplicación del LDM ya fuera publicada por Knight et al (1989). Se muestran los resultados en particular por el interés que demuestra la comparación con el DCM. Los otros cuatro conjuntos de datos provienen de mediciones efectuadas en ríos argentinos. Dos secciones sobre el Alto Río Paraná (provincia de Misiones), correspondientes a un río de grandes dimensiones, una sección sobre el río Colorado (provincia de Neuquén) y otra sección sobre el río La Suela (provincia de Córdoba) representativo de un pequeño río de montaña. Río Severn en Montford Esta serie de datos está constituida por 3 perfiles transversales de velocidades medidos en el río Severn en Montford, Reino Unido. Estos datos fueron obtenidos del trabajo de Knight et al (1989). El río Severn es el río más largo de Inglaterra, con más de 100 puentes que lo atraviesan. Uno de ellos está localizado en Montford, donde se ubica la estación de medición considerada. La cuenca de aporte en esa sección es de 2025 km²; el módulo del río es de 43 m³/s. La precipitación media anual en la cuenca es de 1147 mm. En la Figura 4 puede verse la localización geográfica de la sección en estudio. La sección del río consta de dos planicies de inundación, con diferentes anchos y rugosidades. En la Figura 5 puede verse la sección transversal considerada. INGENIERÍA DEL AGUA · VOL. 12 · Nº 3 SEPTIEMBRE 2005 4 Juan F. Weber, Ángel N. Menéndez y Leticia Tarrab INGENIERÍA DEL AGUA · VOL. 12 · Nº 3 SEPTIEMBRE 2005 5 DISTRIBUCION LATERAL DE VELOCIDADES EN CAUCES NATURALES Figura 4. Localización de la sección de análisis – río Severn en Montford A continuación, en la Tabla 1 se resumen las características principales de los tres casos analizados. Tabla 1. Características de los casos de estudio – río Severn en Montford Caudal (m³/s) Cota (m.s.n.m.) Profundidad (m) Caso 1 330.8 58.092 7.81 Caso 2 220.6 57.202 6.92 Caso 3 188.8 56.732 6.15 Figura 5. Sección transversal – río Severn en Montford Alto río Paraná El río Paraná se forma en territorio brasileño con la unión de los ríos Paranaiba y Grande y tras recorrer aproximadamente 4000 km. confluye con el río Uruguay formando el río de La Plata. Se distinguen dos sectores principales: Alto Paraná y Paraná. El Río Alto Paraná sirve de apoyo al límite internacional desde el punto tripartito en que recibe al río Iguazú hasta su propia confluencia con el río Paraguay. Involucra a dos provincias argentinas: Misiones y Corrientes. El mismo, desarrolla una cuenca de 933000 km2 hasta el sitio de Posadas, que se localiza a 1597 km de su desembocadura en el Estuario del Plata, en el extremo de un trayecto de 350 km a partir de la presa de Itaipú. El tramo del Alto Paraná se encuentra muy encajonado, las profundidades de flujo en las secciones de estudio oscilan entre los 20 y 60 m, mientras que el ancho del río en el nivel de aguas medias varía entre 200 y 2000 m. Las pendientes de las aguas medias varían entre 0.045 y 0.10 m/km y las velocidades medias superficiales de las mismas entre 4.2 y 6 km/h (MOP, 1973).El caudal medio anual medido en Posadas, para una serie de 98 años del período 1901-1998 es de 12313 m3/s, con un caudal máximo diario de 53227 m3/s que corresponde al año 1904-1905 y un caudal mínimo diario de 3906 m3/s registrado en el año 1944-1945. Las estaciones de aforos estudiadas corresponden a las ubicadas en Puerto Libertad y en Eldorado (Figura 6). Los dos perfiles de aforo considerados fueron ubicados y relevados entre los días 17 y 25 de marzo de 1971 con sondas ecógrafas y relacionados altimétricamente con los ceros de las escalas hidrométricas En el anexo se muestran los perfiles relevados. Se realizaron 102 aforos en Puerto Libertad y 106 en Eldorado. Las velocidades se midieron en varias verticales en el cauce principal y a 0.2, 0.6 y 0.8 de profundidad, mientras que en las cercanías de las márgenes se midió a 0.6 de profundidad. En la Tabla 1 se exponen las cantidades de aforos líquidos que se han seleccionado en las estaciones consideradas y a su vez discretizados en estacionarios y no estacionarios. Se consideró como aforos estacionarios (o cuasi-estacionarios) a aquellos en los cuales la variación del nivel de agua durante el aforo era menor que 0.05 m. Además se descartaron algunos aforos en los cuales se indicaba que podría haber ocurrido alguna falla de instrumental. Es de destacar la escasez de aforos realizados en régimen estrictamente estacionario, los cuales son justamente los de mayor interés a la hora de aplicar la metodología descripta. Tabla 1 Aforos efectuados en los distintos regímenes AGUAS Estacionarias No estacionarias TOTAL Pto Libertad 21 8167 102 El dorado 39 57.202 106 Río Severn en Montford 50.0 51.0 52.0 53.0 54.0 55.0 56.0 57.0 58.0 59.0 60.0 0 20 40 60 80 100 120 Progresivas (m) Cota (m) Figura 6. Secciones transversales – Alto río Paraná INGENIERÍA DEL AGUA · VOL. 12 · Nº 3 SEPTIEMBRE 2005 6 Juan F. Weber, Ángel N. Menéndez y Leticia Tarrab INGENIERÍA DEL AGUA · VOL. 12 · Nº 3 SEPTIEMBRE 2005 7 DISTRIBUCION LATERAL DE VELOCIDADES EN CAUCES NATURALES Río Colorado en Buta Ranquil El río Colorado, perteneciente al grupo de los sistemas patagónicos de vertiente atlántica, está formado por la confluencia de los ríos Grande y Barrancas a 835 metros sobre el nivel del mar. Desde sus orígenes en la Cordillera de los Andes, hasta su desembocadura en el Océano Atlántico, presenta una extensión de 1200 kilómetros con dirección NO-SE, de los cuales 920 corresponden al Colorado propiamente dicho. Sus aguas son compartidas por las Provincias de Mendoza, Neuquén, La Pampa, Río Negro y Buenos Aires, que lo convierten en una cuenca hídrica interprovincial. El área de la cuenca imbrífera es de aproximadamente 15300 km2, correspondiente al río Colorado aguas arriba de la estación de aforos de Buta Ranquil (esta estación se encuentra a unos 25 km de la confluencia de los ríos Grande y Barrancas). De régimen nival, su período de creciente se inicia en octubre, alcanza su máximo en diciembre y retoma en marzo niveles menores a la media anual (en Buta Ranquil –alta cuenca-: 148 m3/s; en Pichi Mahuida –curso medio-: 132 m3/s). El caudal medio anual medido en Buta Ranquil, para una serie de 58 años del período 1940-1997 es de 147 m3/s, con un una caudal máximo diario de 1053 m3/s ocurridos en el año rico 1982-1983 y un caudal mínimo diario de 130 m3/s registrados en el año 1968-1969. La estación de aforos Buta Ranquil está ubicada a 850 m.s.n.m en los 37º 06’ de latitud sur y 69º 44’ de longitud oeste. El cauce en esa sección se encuentra encajonado por taludes rocosos de roca basáltica casi verticales lo cual le confiere una gran estabilidad a la sección a lo largo del tiempo. La geometría de la sección con la subdivisión de tramos adoptada para la modelación se muestra en la Figura 7. Figura 7. Sección transversal en la estación de aforos – río Colorado en Buta Ranquil Para la calibración del modelo se utilizaron datos provistos por EVARSA correspondientes al período Enero-Julio del 2003 (período A). Los mismos correspondían 10 aforos en donde se midieron velocidades transversales, relevamientos de fondo y pendiente de las aguas. Además se disponía de 11 valores de caudal asociados a una profundidad correspondientes al período de Mayo-Diciembre del 2002 (período B), los cuales se utilizaron con fines comparativos Río de La Suela La cuenca del Río La Suela, con un área de 132 km2, se encuentra ubicada entre los 64º35’ y 64º46’ de longitud Oeste; y 31º31’ y 31º39’ de latitud Sur, en la provincia de Córdoba, Argentina. El cierre de la cuenca es la estación de aforos La Suela, ubicada a 896 m.s.n.m. El río La Suela, tiene una longitud total en su curso principal de 31 km. Su comportamiento responde a un régimen torrencial, siendo el caudal medio en el período 1971 y 1986 de 1.56 m3/s, mientras que la crecida máxima en el año 1975 alcanzó un pico de 500 m3/s. El cauce en la sección de aforos es de lecho arenoso con carencia de vegetación. La margen derecha está formada por un talud rocoso, y la margen izquierda por una extensa planicie de inundación, por lo que la misma presenta mayor rugosidad superficial. La pendiente media del fondo en la sección de aforos para el año 1982 era de 0.3% (Dasso, 1977). Respecto a los datos de aforos, esta estación cuenta con más de 1000 mediciones tanto para aguas bajas como para crecientes entre los años 1971 y 1984, mientras que se relevaron perfiles transversales en los años 1971, 1975, 1979, 1982. Para profundidades menores a 1,2 m las mediciones se efectuaron a vadeo con un velocímetro Siap y las correspondientes a tirantes superiores mediante un teleférico operado desde la orilla, con un molinete Ott Kempten. Las velocidades en las verticales se midieron a 0.2 y 0.8 de profundidad y a 0.6 para profundidades bajas. Las crecidas se aforaban en forma continua hasta el paso total de la onda. Se utilizó la sección topográfica relevada en el año 1982 y se seleccionaron aquellos valores de velocidad correspondientes a los aforos del período 1980-1983, en el cual se alcanzó cierto equilibrio morfológico en la geometría transversal. Otra condición para la selección de los mismos fue que la profundidad máxima superara los 0.5 m, de modo de reducir la incertidumbre asociada al procedimiento de medición, ya que para niveles 5 m 1 m 1243 5 6 INGENIERÍA DEL AGUA · VOL. 12 · Nº 3 SEPTIEMBRE 2005 6 Juan F. Weber, Ángel N. Menéndez y Leticia Tarrab INGENIERÍA DEL AGUA · VOL. 12 · Nº 3 SEPTIEMBRE 2005 7 DISTRIBUCION LATERAL DE VELOCIDADES EN CAUCES NATURALES de agua tan bajos pequeñas variaciones en el nivel de fondo influyen las velocidades del flujo. En la Figura 8 se muestra el perfil transversal, con la subdivisión de tramos adoptada para la aplicación de los modelos. En total se utilizaron 32 aforos, de los cuales se eligieron 10 para la calibración de los perfiles laterales de velocidad. Figura 8. Geometría de la sección de aforos – río La Suela APLICACIÓN DE LOS MODELOS Se aplicaron los modelos descritos a los casos considerados. El método DCM se aplicó a través del programa HEC-RAS del U.S. Army Corps of Engineers (HEC, 2001) mientras que el SKM se aplicó a través del programa LATERAL, especialmente diseñado para la aplicación de este método a secciones irregulares cuyo algoritmo ha sido detallado en trabajos previos (Weber y Menéndez, 2003a). Río Severn en Montford Se aplicó el modelo HEC-RAS a los 3 casos analizados, con el objeto de obtener una estimación de la distribución lateral de velocidades, considerando los valores para el n de Manning variables entre 0,030 a 0,032 en el cauce principal y 0,019 a 0,100 en las planicies de inundación, calculados a partir de los valores de f obtenidos del trabajo de Knight et al (1989). Para la instrumentación del programa HEC-RAS, se consideró una geometría prismática de 10 km de longitud, con una pendiente longitudinal de 2 x 10-4, la cual es también un dato obtenido de la misma referencia. En los tres casos, se consideró conocida la profundidad de flujo y el caudal (que fueron establecidas como condición de borde del problema), por lo cual se debió verificar la condición de flujo uniforme, la cual no fue posible cumplir. De todos modos, la variación de tirantes y velocidades a lo largo de los 10000 m simulados es pequeña (me- nor al 10 % en el caso de los tirantes y menor al 15 % en el caso de las velocidades medias). Se implementó la solución analítica de Shiono y Knight (SK) descripta para los 3 casos de estudio mencionados. Para ello, y en cada caso de estudio, se utilizó el programa LATERAL (Weber y Menéndez, 2003a) considerando para la rugosidad (n de Manning) los valores anteriormente mencionados y para la viscosidad adimensional de torbellino l los siguientes valores: 0,07 para el cauce principal, 0,2 para los taludes y 3,0 para las planicies de inundación (Knight et al, 1989). En los gráficos siguientes (Figura 9 a Figura 11) se presentan los perfiles transversales de velocidad calculados por ambos modelos junto con los datos de campaña para los tres casos considerados. Figura 9. Resultados caso 1 – río Severn en Montford Figura 10. Resultados caso 2 – río Severn en Montford Figura 11. Resultados caso 3 – río Severn en Montford 1 2 3 4 5 6 7 8 9 10 1 s5 m Caso 1 - 0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 1.8 0 20 40 60 80 100 120 Progresivas (m) Velocidad (m/s) campaña SKM fondo DCM Caso 2 - 0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 0 20 40 60 80 100 120 Progresivas (m) Velocidad (m/s) campaña SKM fondo DCM Caso 3 - 0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 0 20 40 60 80 100 120 Progresivas (m) Velocidad (m/s) campaña SKM fondo DCM INGENIERÍA DEL AGUA · VOL. 12 · Nº 3 SEPTIEMBRE 2005 8 Juan F. Weber, Ángel N. Menéndez y Leticia Tarrab INGENIERÍA DEL AGUA · VOL. 12 · Nº 3 SEPTIEMBRE 2005 9 DISTRIBUCION LATERAL DE VELOCIDADES EN CAUCES NATURALES Como puede verse en las figuras anteriores, el modelo DCM representa razonablemente la distribución transversal de velocidades, presentando dificultades en la estimación del pico de velocidades que se produce en el cauce principal. En efecto, se presentan errores por defecto de 9,6%, 13,1% y 13,5% para los casos 1, 2 y 3 respectivamente; lo que muestra que el error en la estimación del pico de velocidades es mayor cuanto menor la profundidad de flujo. En general puede verse un buen ajuste de la solución analítica a los datos de campaña, tal como ya lo habían presentado Knight et al (1989). Río Alto Paraná en Eldorado En la primera etapa se calibraron los coeficientes del modelo SK con los datos de los aforos considerados estacionarios. Luego de la calibración, se aplicaron los modelos en 10 casos seleccionados, con la finalidad de mostrar el desempeño de los mismos en distintas condiciones de flujo. La viscosidad adimensional de torbellino, l, se estimó con la ecuación (18), la cual fue establecida en base a una formulación empírica. .(18) Donde: B Ancho del canal principal H Profundidad del canal principal Dr Profundidad relativa de las planicies: (H-h)/H Esta ecuación, de acuerdo a Tarrab (2004) y Tarrab y Weber (2004) ha permitido predecir la viscosidad adimensional de torbellino adecuadamente no sólo en ambas secciones del río Paraná sino también en los ríos Colorado y La Suela. El coeficiente de rugosidad en el cauce principal adoptado (tramos 5 y 6) responde a: (19) Donde Ym es la profundidad media del tramo. De los 10 casos estudiados, sólo se presentan aquí los resultados correspondientes a tres de ellos, por razones de brevedad. En ese sentido, se han seleccionado los casos 1, 7 y 10 asociados a condiciones de aguas bajas, medias y altas, respectivamente. En las figuras siguientes (Figura 12 a Figura 14) se muestran los perfiles laterales de velocidad predichos tanto por el modelo SKM como por el DCM (HEC-RAS) a los casos 1, 7 y 10. Figura 12. Resultados caso 1, río Alto Paraná en Eldorado Figura 13. Resultados caso 7 – río Alto Paraná en Eldorado Figura 14. Resultados caso 10, río Alto Paraná en Eldorado Se aprecia una buena estimación de la distribución lateral de velocidades por el modelo SKM, tanto para aguas bajas, como para aguas altas, a pesar de que en el cauce principal, específicamente en el tramo 10, se observa un comportamiento atípico de las velocidades medias. El valor de la velocidad decrece abruptamente en este tramo, manteniéndose con un valor similar al del tramo anterior (aumentando en algunos casos) en el tramo correspondiente al talud derecho. Caso 1 0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 -315 -215 -115 -15 85 185 Progresiva (m) Velocidad (m/s) SKM DCM campaña Fondo Caso 7 0 0.5 1 1.5 2 2.5 3 -315 -215 -115 -15 85 185 Progresiva (m) Velocidad (m/s) SKM DCM campaña Fondo Caso 10 0 1 2 3 4 5 -315 -215 -115 -15 85 185 Progresiva (m) Velocidad (m/s) SKM DCM campaña Fondo INGENIERÍA DEL AGUA · VOL. 12 · Nº 3 SEPTIEMBRE 2005 8 Juan F. Weber, Ángel N. Menéndez y Leticia Tarrab INGENIERÍA DEL AGUA · VOL. 12 · Nº 3 SEPTIEMBRE 2005 9 DISTRIBUCION LATERAL DE VELOCIDADES EN CAUCES NATURALES Este comportamiento se podría atribuir a la forma del perfil en ese tramo, y a la posible existencia de corrientes secundarias ocasionadas por este fuerte cambio de pendiente lateral en el cauce principal. De todas maneras, la calibración de dicho coeficiente en ese tramo (aumentándolo con un factor de 2.1 con respecto al del cauce principal) resultó efectiva en el momento de predecir el perfil de velocidades. En cuanto al modelo DCM, el ajuste no ha sido tan exitoso. En aguas bajas, el modelo produce una buena estimación, aunque sobreestima en pequeña magnitud a las velocidades en el cauce principal, incrementándose esta diferencia a medida que los caudales aumentan. En el caso 10 (Figura 14) las velocidades en el cauce principal son aproximadamente el doble del valor de las medidas en campo. Río Alto Paraná en Puerto Libertad En una primera instancia, se calibró el modelo SKM con los aforos estacionarios (o cuasi-esta- cionarios). Los valores de n en el cauce principal se establecieron como una función de la profundidad de flujo local de la forma: (20) . En los demás tramos (taludes, márgenes y planicies) se adoptó un n de Manning igual a 0.06. Luego de la calibración con los datos seleccionados para tal finalidad, se aplicaron los modelos en 10 casos de manera de determinar y graficar la distribución lateral de velocidades en algunas situaciones representativas de los aforos totales. Los casos de aplicación se seleccionaron de manera de que los mismos fueran en lo posible estacionarios y que a su vez cubrieran el rango experimental de la muestra total. La viscosidad adimensional de torbellino, l, se estimó con la ecuación (18). La variación del valor de l en los tramos que conforman el cauce principal (cauce principal y taludes) no es significativa en el rango de profundidades analizado, sin embargo, en los demás tramos, las variaciones son mayores. De los 10 casos estudiados, sólo se presentan aquí los resultados correspondientes a tres de ellos, por razones de brevedad. En ese sentido, se han seleccionado los casos 1, 7 y 10 asociados a condiciones de aguas bajas, medias y altas, respectivamente. En las figuras siguientes (Figura 15 a Figura 17) se muestran los perfiles laterales de velocidad predichos tanto por el modelo SKM como por el DCM (HEC-RAS) a los casos 1, 7 y 10. Figura 15. Resultados caso 1, río Alto Paraná en Puerto Libertad Figura 16. Resultados caso 7, río Alto Paraná en Puerto Libertad Figura 17. Resultados caso 10, río Alto Paraná en Puerto Libertad Se puede observar un buen desempeño del modelo SKM en todos los casos, inclusive en los casos de inestacionariedad del flujo y eventos extraordinarios (crecidas) como ocurre en el caso de aplicación número 10. Los quiebres en el perfil de velocidades en los tramos cercanos a las márgenes, se deben al abrupto incremen- to de la viscosidad adimensional de torbellino. Caso 1 0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 -200 -150 -100 -50 0 50 100 150 Progresiva (m) Velocidad (m/s) SKM DCM campaña Fondo Caso 7 0 0.5 1 1.5 2 2.5 3 3.5 4 -200 -150 -100 -50 0 50 100 150 Progresiva (m) Velocidad (m/s) SKM DCM campaña Fondo Caso 1 0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 -200 -150 -100 -50 0 50 100 150 Progresiva (m) Velocidad (m/s) SKM DCM campaña Fondo