Saint Venant PDF

Tecnura
ISSN: 0123-921X
[email protected]
Universidad Distrital Francisco Jos de
Caldas
Colombia
Amars Castro, Gloria Estefany; Guerrero Barbosa, Thomas Edison; Snchez Ortiz, Edgar
Antonio
Comportamiento de las ecuaciones de Saint-Venant en 1D y aproximaciones para
diferentes condiciones en rgimen permanente y variable
Tecnura, vol. 19, nm. 45, julio-septiembre, 2015, pp. 75-87
Universidad Distrital Francisco Jos de Caldas
Bogot, Colombia
Disponible en: http://www.redalyc.org/articulo.oa?id=257040047007
Cmo citar el artculo

Nmero completo
Sistema de Informacin Cientfica
Ms informacin del artculo Red de Revistas Cientficas de Amrica Latina, el Caribe, Espaa y Portugal
Pgina de la revista en redalyc.org Proyecto acadmico sin fines de lucro, desarrollado bajo la iniciativa de acceso abierto
Tecnura
http://revistas.udistrital.edu.co/ojs/index.php/Tecnura/issue/view/680
DOI: http://dx.doi.org/10.14483/udistrital.jour.tecnura.2015.3.a06
Investigacin
Comportamiento de las ecuaciones de Saint-Venant en 1D y aproximaciones

para diferentes condiciones en rgimen permanente y variable
Performance of the Saint Venant equations in 1D and approaches to different
conditions in steady and variable state
Gloria Estefany Amars Castro*, Thomas Edison Guerrero Barbosa**, Edgar Antonio Snchez Ortiz***
Fecha de recepcin: 6 de septiembre de 2014 Fecha de aceptacin: 6 de abril de 2015
Como citar: Amars Castro, G. E., Guerrero Barbosa, T. E., & Snchez Ortiz, E. A. (2015). Comportamiento de las
ecuaciones de Saint-Venant en 1D y aproximaciones para diferentes condiciones en rgimen permanente y va-
riable. Revista Tecnura, 19(45), 75-87. doi: 10.14483/udistrital.jour.tecnura.2015.3.a06
RESUMEN simple que produzca soluciones tiles a las necesi-

La importancia del comportamiento de las ecuacio- dades de la ingeniera, y se encuentra que dentro
nes de Saint-Venant y sus implicaciones en el anlisis del anlisis de redes de alcantarillado hay muchos
del flujo aplicado a modelos de trnsito hidrulico factores que pueden afectar la solucin de las ecua-
es el tema principal de esta investigacin. Partien- ciones de Saint-Venant.
do de las aplicaciones y usos de dichos modelos Palabras clave: ecuaciones de Saint-Venant, mode-
matemticos, se determina y se evala el compor- los matemticos, onda dinmica, pendiente subcr-
tamiento de estas ecuaciones en 1D y sus aproxima- tica, pendiente supercrtica.
ciones para diferentes condiciones posibles de una
situacin real mediante la experimentacin num- ABSTRACT
rica computacional para un tramo de un cauce con The importance of the behavior of the Saint Ve-
una seccin hidrulica prismtica, en casos en los nant equations and its implications for flow analy-
cuales la pendiente del canal presenta tres tipos de sis applied to models of hydraulic transit is the main
condiciones: supercrtica, subcrtica y pendiente ho- theme of this research, where the basis of the appli-
rizontal, para un caudal de referencia. El objetivo cations and uses of these mathematical models is
principal de esta investigacin es encontrar una de- determined and the performance is evaluated these
mostracin de la aplicabilidad de las ecuaciones de equations in 1D and its possible approaches for di-
Saint-Venant para diferentes condiciones de acerca- fferent conditions of a real situation using compu-
miento a esta problemtica desde el punto de vista tational numerical experimentation for a section of
de la solucin de las ecuaciones de Saint-Venant en a prismatic channel with a hydraulic section, in ca-
una dimensin. Se propone un esquema de solucin ses where the channel slope presents three types of
* Ingeniera civil, estudiante de la Maestra en Ingeniera Civil con nfasis en Recursos Hdricos, Universidad del Norte, Departamento de
Ingeniera Civil y Ambiental, Barranquilla, Colombia. Contacto: [email protected]
** Ingeniero civil, magster en Ingeniera Civil. Docente Auxiliar, Universidad Francisco de Paula Santander Ocaa, Departamento de Ingeni-
era Civil, Ocaa, Colombia. Contacto: [email protected]
*** Licenciado en Fsica y Matemticas, magster en Educacin. Director seccional, Universidad Francisco de Paula Santander Ocaa, Ocaa,
Colombia. Contacto: [email protected]
Tecnura p-ISSN: 0123-921X e-ISSN: 2248-7638 Vol. 19 No. 45 Julio - Septiembre 2015 pp. 75-87
[ 75 ]
Comportamiento de las ecuaciones de Saint-Venant en 1D y aproximaciones para diferentes condiciones en rgimen permanente y variable
Amars Castro, G. E., Guerrero Barbosa, T. E., & Snchez Ortiz, E. A.
conditions: supercritical, subcritical and horizontal useful solutions to the needs of engineering, finding
slope for a design rate. The main objective of this re- that within the sewer network analysis there are
search is to find a demonstration of the applicability many factors that can affect the solution of the Saint
of the Saint Venant equations for different conditions Venant equations.
of approaching this problem from the point of view Keywords: mathematical modeling, Saint-Venant
of the solution of the Saint Venant equations in one equations, slope subcritical, supercritical slope,
dimension. Scheme simple solution that produces wave dynamics.
INTRODUCCIN inestables debido a la alta acumulacin de erro-

res (Tsai, 2003). En estos casos generalmente se su-
En Colombia los modelos de trnsito permanente giere el uso de un esquema de solucin mediante
son de uso comn en las modelaciones en ingenie- la formulacin y desarrollo de la onda cinemti-
ra. Software de uso frecuente como el Sistema de ca, de tal forma que considere la no linealidad en
Anlisis de Ros, desarrollado por el Centro de In- flujos bajo estas condiciones. Estas modelaciones
geniera Hidrolgica de California, EUA (HEC-RAS, matemticas son consideradas como una aproxi-
por su sigla en ingls), permite determinar la altura macin a la realidad, por lo cual es primordial es-
requerida para construir estructuras de control como tablecer las condiciones de frontera que permitan
puentes, diques, vertederos y obras de proteccin en determinar las condiciones aproximadas del com-
aguas bajas, redes hdricas y flujo libre; sin embargo, portamiento del flujo y del canal, ya que cuanto
cuando la lmina del flujo crece rpidamente, los ms se adapten las ecuaciones de partida y la pro-
errores de los esquemas no linealizados aumentan. gramacin numrica para su resolucin al fenme-
Las ecuaciones de Saint-Venant (1871) han no real que se est estudiando, mejores resultados
demostrado ser un modelo matemtico til con se obtendrn.
aplicaciones en varias ramas de la ingeniera hi- Generalmente en los estudios hidrulicos, como
drulica y costera; estos modelos se basan en el es el caso de los efectos de la propagacin de ave-
desarrollo de ecuaciones diferenciales encargadas nidas en ros, clculos de velocidades y niveles de
de modelar variacin de caudal y nivel de un flu- agua, son usados mtodos de abstracciones uni-
jo en un espacio unidimensional y el tiempo de dimensionales en rgimen permanente gradual-
manera no permanente en canales abiertos. En ge- mente variado y fondo fijo, pues se les considera
neral estas modelaciones se presentan en tres di- como una herramienta suficiente para determinar
mensiones, pero las aplicaciones prcticas suelen el comportamiento del ro y se deja a un lado la
utilizar una aproximacin unidimensional en la di- evolucin temporal en el proceso. Aplicaciones de
reccin de mayor relevancia, que usualmente sue- partculas hidrodinmicas con modelos matemti-
le ser a lo largo del canal en la direccin de flujo. cos basados en las ecuaciones de Saint-Venant en
En la revisin bibliogrfica se encuentran evi- canales abiertos son mostradas por (Chang, Kao,
dencias de otros estudios que se encargan de pro- Chang, & Hsu, 2011).
poner diferentes escenarios con el fin de explicar La onda dinmica se considera como el mode-
el comportamiento del flujo mediante las ecuacio- lo de trnsito distribuido ms completo existente,
nes de Saint-Venant y establecer su aplicabilidad ya que tiene en cuenta todas las fuerzas que inte-
para determinar el comportamiento del flujo. ractan en el trnsito de la creciente (inerciales,
En las modelaciones de flujo en canales de presin, gravitacional, friccin) y tiene en cuenta
alta pendiente (So > 5%) se generan soluciones los efectos de remanso en el clculo; en su forma
[ 76 ]
completa, estas ecuaciones se consideran de solu- mediante condiciones de piscina nivelada, pen-
cin compleja aun con el uso de las ms modernas diente subcrtica y condicin aguas abajo profun-
herramientas computacionales (Cea, 2008). didad crtica, pendiente supercrtica y condicin
A pesar de la complejidad de estas ecuaciones, aguas abajo profundidad normal.
son estudiadas con el fin de entenderlas y con esto
mejorar las ideas en el campo de la mecnica de METODOLOGA
los fluidos en las soluciones de flujos particulares y
para examinar leyes de modelamiento matemtico Teniendo en cuenta el concepto de onda cine-
(White, 1991). mtica como un modelo simplificado del flujo
Otros estudios aplicativos referentes a la imple- no permanente en canales abiertos y conociendo
mentacin del esquema numrico en volmenes que las ondas dinmicas resuelven las ecuacio-
finitos para las ecuaciones de flujo Saint-Venant nes completas, se concluira que las ondas din-
han permitido conocer el perfil de flujo superficial micas resuelven las ecuaciones completas, pues
y la infiltracin en el suelo en sistemas de riego, se esperara que la onda dinmica presentara me-
como es el caso de Pino, Meja y Chvarry (2012), jor comportamiento que la onda cinemtica de-
quienes determinaron, tras un proceso de simula- bido a que la onda dinmica contiene trminos
cin de varios casos, que la pendiente longitudinal inerciales (aceleracin local y aceleracin con-
tiene un papel importante en el flujo para las mel- vectiva), fuerza gravitacional, fuerza de friccin;
gas; segn su modelo, la pendiente debe ser me- por el contrario, la onda cinemtica nicamente
nor a 0,001 m/m. Las pendientes mayores generan contiene los trminos de fuerza gravitacional y
flujos rpidos o supercrticos. de fuerza de friccin pero segn estudios realiza-
Al aplicar el anlisis de flujo de rgimen varia- dos los ltimos 100 cien aos indican lo contra-
ble es imprescindible considerar las restricciones rio. Lighthill y Whitham (1955) afirman que bajo
que permiten el desarrollo de las ecuaciones de las condiciones apropiadas para el flujo de ave-
Saint-Venant, por lo tanto es necesario utilizar un nida, las ondas dinmicas se atenan muy rpi-
software que permita modelar las situaciones en damente, y son las ondas cinemticas, viajando
que sea aplicable la suposicin de onda cinemti- a una velocidad menor, las que asumen el papel
ca como la onda dinmica. preponderante.
Aplicaciones de las ecuaciones de Saint-Venant
contrarias a los casos ya conocidos, como en redes Tipos de trnsito hidrulico
fijas convencionales o mtodos de partculas en mo-
vimiento, se evidencian en Yao-Hsin (2013), quien Las ecuaciones de Saint-Venant tienen varias for-
muestra una formulacin a partir de la reasignacin mas simplificadas, cada una de las cuales define
de partculas computacionales a lo largo de las cur- un modelo de trnsito distribuido unidimensional
vas caractersticas. En la prediccin de las variables (Zapata Gaviria, 2012). Aunque las fuerzas expre-
de flujo para diferentes condiciones tiene gran im- sadas en la ecuacin de momentum siempre estn
portancia desarrollar modelos y simulaciones nu- presentes en cualquier problema de hidrodinmi-
mricas hidrodinmicas con precisin adecuada ca, se puede dar el caso de que no todas sean re-
que permitan la correcta planificacin de obras. levantes. Esto hace que sea posible tener varios
El objetivo principal de esta investigacin es tipos de trnsito hidrulico de crecientes, como se
evaluar el comportamiento de las ecuaciones de muestra en la ecuacin (1).
Saint-Venant para condiciones en rgimen perma-
1 1
nente y variable, bajo diferentes escenarios de es- (1)

tudio y modificando las condiciones de frontera
[ 77 ]
Modelo de onda cinemtica una de continuidad y otra de momentum. La solu-

cin de las ecuaciones requiere de mtodos num-
El modelo de onda cinemtica es la mejor alterna- ricos directos, los cuales convierten las ecuaciones
tiva en la modelacin del trnsito hidrulico en re- diferenciales en ecuaciones de diferencias finitas
des debido a la simplicidad de su estructura, buen para la conservacin de la masa y del momentum
comportamiento y representatividad ante una mo- lineal, las cuales pueden ser lineales o no linea-
delacin, sin embargo el que desprecie ciertas les. Los clculos de las variables desconocidas se
condiciones como trminos de aceleracin local, hacen simultneamente en funcin del espacio y
aceleracin convectiva y fuerzas de presin, limita del tiempo. En el caso del trnsito hidrulico de
en cierta medida al modelo para ciertas circunstan- crecientes, las variables desconocidas son: Q y A.
cias de escenarios dependiendo del caso especifico Si se tiene una funcin f(x) continua en un inter-
(condiciones de pendiente y velocidad entre otras). valo (a,b), el valor de f(x+x) puede ser calculado
En este caso la ecuacin de Saint Venant de utilizando el valor de f(x), que se supone conoci-
conservacin de momentum se reduce a la ecua- do, y las derivadas de la funcin f en el punto x co-
cin (2) y la ecuacin (3). nocido. Utilizando series de Taylor y resolviendo
para f(x), se obtiene la ecuacin (4).
(2)

(4)
(3)

Esta ecuacin puede ser utilizada para solucio-
Lo anterior implica que el flujo puede ser trata- nar la ecuacin diferencial parcial, en donde la
do como uniforme no permanente. Esto significa derivada es funcin de los valores conocidos de la
que no hay efectos transientes de curvas de reman- funcin inicial.
so. En este caso existe un balance entre el efecto
de la aceleracin debido a la gravedad y el efecto Solucin por diferencias finitas
de retraso de las fuerzas de friccin.
La celeridad de onda es la velocidad con la El esquema de diferencias finitas debe solucionar
cual se mueve el frente de onda con respecto a los las ecuaciones de Saint Venant. La variable depen-
bordes del conducto. En el caso de la onda cine- diente nicamente es Q. La ecuacin de continui-
mtica, los trminos de aceleracin y de presin dad, por ser una ecuacin diferencial parcial, se
son muy pequeos comparados con los gravitacio- puede aplicar al mtodo de las diferencias finitas,
nales y los de friccin; es decir, el movimiento de con el fin de conocer Q para cualquier lugar x del
la onda est descrito ms por la ecuacin de con- conducto en cualquier instante de tiempo.
tinuidad que por la de momentum. La otra variable dependiente (y o A) se calcu-
la utilizando Q(x,t) y una ecuacin de resistencia
Esquema numrico fluida, tal como la ecuacin de Gauckler Man-
ning. Luego de un sencillo procedimiento mate-
Las ecuaciones de Saint Venant para la onda cine- mtico se llega a la ecuacin (5).
mtica son dos ecuaciones diferenciales parciales,

2 2

(5)

2
[ 78 ]
Donde los valores a y b dependen de la geome- tales como la propagacin de la onda sobre una pen-
tra del conducto, y los trminos de x y t deben diente horizontal, subcrtica o supercrtica. En cuan-
estar relacionados mediante la condicin de Cou- to al clculo del trnsito distribuido de la onda sobre
rant correspondiente a la ecuacin (6). el canal evaluado, fue necesario implementar una
rutina (algoritmo) que lograra la resolucin numrica
(6) de las ecuaciones diferenciales de momento y conti-

nuidad de una forma rpida y estable. Esta rutina fue
elaborada en Matlab 2011 y consider dentro de su
APLICACIN PARA EVALUAR LAS estructura la posibilidad de evaluar diferentes con-
ECUACIONES DE SAINT-VENANT diciones de frontera, tales como: hidrograma y lim-
nigrama en la frontera aguas arriba e hidrogramas,
Mediante el uso de las ecuaciones de Saint-Venant limnigramas, condiciones de flujo normal y flujo cr-
y los mtodos de solucin tanto implcito como tico sobre la frontera aguas abajo.
explcito, se procedi a desarrollar los siguientes El anlisis del escenario propuesto indica que el
casos de estudios y para cada uno se obtuvo re- caudal pico mostrado en la grfica del hidrograma
sultados en una situacin en particular. La seccin de entrada es un tiempo de 3,65 hr (condiciones t-
transversal del canal de irrigacin se muestra en la picas de la ciudad de Barranquilla); para la progra-
figura 1, as como otras caractersticas de esta se macin es necesario que el tiempo de modelacin
presentan a continuacin: sea superior por lo menos al caudal pico del hi-
drograma es por esto que se decidi tomar 5 h de
Tipo de cauce: canal prefabricado. tal forma que supere el caudal pico mostrado en el
Uso: irrigacin. hidrograma; adicionalmente se utiliz un diferen-
Forma de la seccin transversal: trapezoidal. cial de tiempo (t) de 0,1 hrs y un diferencial de
Relacin de taludes laterales: 2H:1V. posicin (x) de 305 m. Los casos sobre los cuales
Base: 6,1 m. se evalu la ecuacin de Saint Venant se presentan
Longitud: 2743,2 m. a continuacin.
Cota inicial: 29,1 m.s.n.m.
Rugosidad (n de Manning): 0,025. Caso 1: flujo horizontal
La pendiente longitudinal utilizada vari segn Para este escenario de modelacin se plan-
las condiciones hidrulicas que se deseaban evaluar, te una pendiente de fondo de S0=0,0000005,
Figura 1. Seccin Transversal del Canal de Irrigacin.
Fuente: elaboracin propia.
[ 79 ]
suficientemente baja, simulando las condiciones la simulacin de flujo horizontal correspondiente

de un embalse. Las condiciones de frontera para a las figuras 2 y 3.
simulacin fueron: Al igual que las simulaciones de los escenarios
anteriores, esta se ejecut con el modelo de onda di-
Aguas arriba: hidrograma de caudales. nmica, seleccionando como valor de =0,75. Para
Aguas abajo: profundidad normal. este caso las condiciones iniciales definidas son:
A continuacin se presentan las condiciones de Caudal inicial: 10 m3/s.

frontera aguas arriba y aguas abajo definidas para Profundidad inicial: 7,119 m.
Figura 2. Condicin de frontera aguas arriba, simulacin de flujo horizontal.
Figura 3. Simulacin de flujo horizontal.
[ 80 ]
Caso 2: flujo subcrtico Limnigrama de niveles aguas arriba.

Limnigrama de niveles aguas abajo.
Se ha considerado un flujo subcrtico, con pen-
diente de fondo de S0=0,0005; las condiciones de Las figuras 4 y 5 muestran las condiciones de
frontera aplicadas en esta simulacin fueron: frontera aguas arriba y aguas abajo utilizadas para
esta simulacin.
Figura 4. Condicin de frontera aguas arriba, simulacin de flujo subcrtico.
Figura 5. Condicin de frontera aguas arriba, simulacin de flujo subcrtico.
[ 81 ]
Caso 3: flujo supercrtico Hidrograma de caudales aguas abajo.
Se ha considerado un flujo supercrtico con pen- Las figuras 6 y 7 muestran las condiciones de
diente de fondo de S0=0,05; las condiciones de frontera aguas arriba y aguas abajo utilizadas para
frontera aplicadas en esta simulacin fueron: esta simulacin.
Hidrograma de caudales aguas arriba.
Figura 6. Condicin de frontera aguas arriba, simulacin de flujo supercrtico.
Fuente: elaboracin propia
Figura 7. Condicin de frontera aguas abajo, simulacin de flujo supercrtico.
[ 82 ]
RESULTADOS la entrada y a las 5 horas el caudal en la salida es de

23,3m3/s y en la entrada el caudal es de 5,66 m3/s.
Caso 1 Mientras tanto, en el limnigrama notamos que la va-
riacin de la profundidad en la entrada es muy simi-
A continuacin se presenta el trnsito del hidrograma lar a la variacin de la profundidad en la entrada.
de caudales y del limnigrama de profundidades a lo Anteriormente observamos las variaciones de
largo del canal para las condiciones iniciales y fina- los caudales y las profundidades respecto al tiem-
les, como se ve en la figura 8. Respecto a los hidro- po; ahora bien, en la figura 9 se presenta la varia-
gramas observamos que en el tiempo 3,65 h el caudal cin de los caudales y las profundidades respecto
pico de la salida (26,82 m3/s) cruza el hidrograma de al tiempo y al espacio.
Figura 8. Hidrogramas y limnigramas de entrada y salida para simulacin de flujo horizontal.
Figura 9. Resultados de la simulacin de flujo horizontal con modelo de onda dinmica.
[ 83 ]
Caso 2 El hidrograma de salida alcanza el mismo caudal

pico que el hidrograma de entrada (62,3 m3/s) pico
La simulacin de este escenario se llev a cabo con el despus de 4 horas de simulacin. Se aprecia que
modelo de onda dinmica. En la figura 10 se presen- el modelo de onda dinmica simula correctamente
ta el hidrograma de entrada, el hidrograma de sali- el escenario de flujo subcrtico unidimensional para
da, las profundidades de entrada y las profundidades el escenario de limnigramas como condiciones de
de salida. En la figura 11 se muestra cmo transitan frontera. Los resultados no tienen inconsistencias en
el hidrograma y el limnigrama a lo largo del canal. el hidrograma de salida con un valor de = 0,75.
Figura 10. Hidrogramas y limnigramas de entrada y salida para simulacin de flujo subcrtico.
Figura 11. Resultados de la simulacin de flujo subcrtico con modelo de onda dinmica.
[ 84 ]
Caso 3 largo del canal. El hidrograma de salida alcanza

el mismo caudal pico que el hidrograma de en-
La simulacin de este escenario se llev a cabo trada (62,3 m3/s) pico en menos de 0,3 horas de
con el modelo de onda dinmica. En la figura 12 simulacin debido a la gran velocidad en la que
se presenta el hidrograma de entrada, el hidrogra- se mueve la onda. Se pudo observar que el tiempo
ma de salida, las profundidades de entrada y las de respuesta del sistema a la onda fue tal que los
profundidades de salida. En la figura 13 se muestra hidrogramas de niveles y caudales de las fronte-
cmo transitan el hidrograma y el limnigrama a lo ras estn muy cercanos entre s. Se aprecia que el
Figura 12. Hidrogramas y limnigramas de entrada y salida para simulacin de flujo supercrtico.
Figura 13. Resultados de la simulacin de flujo supercrtico con modelo de onda dinmica.
[ 85 ]
modelo de onda dinmica simula correctamente (1988), Singh (1996), Tsai (2003), y Macchione y
el escenario de flujo supercrtico para las condi- Morelli (2003).
ciones de frontera dadas. Los resultados presentan El modelo basado en la simplificacin de la
una ligera inestabilidad en el arranque de la simu- onda dinmica permite simular los frentes de on-
lacin, y es debido a que el modelo para esos in- das manteniendo la continuidad en la masa y adi-
tervalos de tiempo no convergi en una solucin cionalmente disminuyendo la atenuacin en el
adecuada, lo cual era de esperarse para el caso de caudal pico. Sin embargo, para condiciones de
flujo evaluado. flujo en canales, donde se pueda generar un fren-
te de onda, las ecuaciones de Saint-Venant no son
CONCLUSIONES vlidas en las vecindades de la discontinuidad.
De acuerdo con los resultados obtenidos para
El modelo de la onda dinmica demuestra ser apli- cada aplicacin, se debe cuestionar el uso de los
cable a todos los casos, a menos que sea posible softwares de manera que se tenga en cuenta en las
demostrar que los modelos o aproximaciones de la modelaciones todos los parmetros que afectan las
onda difusa y la onda cinemtica son adecuados ecuaciones de Saint-Venant.
para representar el flujo. La evaluacin y modelacin del comportamien-
Las condiciones para las cuales son aplicables to del flujo es importante para la aplicacin de la
las ecuaciones de Saint-Venant que se basan en ingeniera hidrulica y costera, es por esto que las
la solucin de la ecuacin de continuidad acom- modelaciones que se realizan deben ajustarse a las
paada por la ecuacin de momentum lineal son: situaciones reales encontradas en campo para de-
para condiciones de pendiente menores al 10% y finir las condiciones de frontera de la estructura
situaciones hidrulicas que cuenten con pendien- hidrulica y sus limitaciones, con el fin de que las
tes pronunciadas, se puede hacer uso del mtodo modelaciones simulen una situacin cercana a la
de la onda cinemtica;, si por el contrario las con- real presentada.
diciones a modelar cuentan con pendientes sua- Los resultados sugieren que los efectos de los
ves (terrenos planos con pendiente inferiores al errores de los parmetros de entrada en los resulta-
0,05%), se debe hacer uso del mtodo de la onda dos de salida son ms significativos en situaciones
dinmica, de tal manera que ajuste de forma ade- especiales, como los valores ms bajos de coefi-
cuada a las condiciones del fenmeno que se de- ciente de rugosidad de Manning y una pendiente
sea modelar. base ms pronunciada en las caractersticas de un
Durante el desarrollo de la investigacin se de- hidrograma de diseo, tal como lo sugiere Barti,
muestra que existen factores que afectan la solu- Sajjad y Akbari (2012).
cin de las ecuaciones, entre los cuales estn el
x (variacin de la distancia) que se escoja para el FINANCIAMIENTO
clculo, las condiciones de frontera, la geometra
y el tipo de estructura hidrulica. La investigacin es resultado del financiamiento
Para verificar el efecto del x, se model con entregado por la Universidad Francisco de Paula
el programa Matlab un tramo con una longitud de Santander Ocaa (UFPSO) a travs de la Divisin
2743,2 m y se variaron los x para el clculo. El de Investigacin y Extensin (DIE) del proyecto ti-
uso de las ecuaciones simplificadas de la onda ci- tulado Estudio del sistema de drenaje del campus
nemtica, en el trnsito de crecientes para cana- de la Universidad Francisco de Paula Santander
les de alta pendiente, se muestran bastante tiles, Ocaa y diseo de estructuras hidroeconmicas
tal y como lo reportan Chow, Maidment y Mays mediante el principio de incertidumbre.
[ 86 ]
REFERENCIAS Macchione, F., & Morelli, M. (2003). Pactical aspect

comparing shock-capturing schemes for dam break
Barati, R., Sajjad, R., & Akbari, G. (2012). Analysis of problems. Journal hydraulic engineering , 129 (3),
dynamic wave model for flood routing in natu- 187-195.
ral rivers. Water Science and Engineering , 5 (3), Pino Vargas, E., Meja Marcacuzo, J. A., & Chvarry Ve-
243-258. larde, E. (2012). Un Modelo Numrico 1D en Vo-
Cea, L. (2008). Modelacin matemtica en lecho fijo del lmenes Finitos para la Solucin de las Ecuaciones
flujo en ros. Modelos 1D y 2D en rgimen perma- de Flujo e Infiltracin del Riego por Gravedad en
nente y variable. Jornadas Tcnicas sobre Hidruli- Melgas. Revista ECIPer , 9 (1), 22-30.
ca Fluvial, CEDEXMinisterio de Fomento. Madrid. Singh, V. P. (1996). Kinematic wave modeling in water
Chang, T. J., Kao, H. M., Chang, K. H., & Hsu, M. H. resources: Sourface water hydrology. New York:
(2011). Numerical simulation of shallow water John Wiley.
dam break flows in open channels using smoothed Tsai, C. W. (2003). Applicability of kinematic, noni-
particle hydrodynamics. Journal Hydrology , 408, nertia, and quasi-steady dynamic wave models for
7890. bend-flow simulation. Journal of Hydraulic Engi-
Chow, V. T., Maidment, D. R., & Mays, L. (1988). Applied neering , 613-627.
Hydrology. Estados Unidos: Mac Graw Hill. White, M. F. (1991). Viscous Fluid Flow. Estados Unidos:
De Saint-Venant, B. (1871). Theorie du mouvement non University of Rhode Island.
permanent des eaux, avec application aux crues de Yao-Hsin, H. (2013). A characteristic particle method
rivieras at a lintroduction des marces dans leur lit. for the Saint-Venant equations. Computers & Fluids,
Compt Rend Acad Sci , 147154. 76, 5872.
Lighthill, M. J., & Whitham, G. B. (1955). On kinetic Zapata Gaviria, J. F. (2012). Propuesta metodolgica
wave II: a theory of traffic flow on crowded roads. para la conceptualizacin de sistemas de drenaje
Proceedings of the Royal Society of London, Series urbano. Revista EPM , 6, 20-31.
A, 229, pp. 317- 345.
[ 87 ]

Saint Venant PDF

Cargado por

Copyright:

Formatos disponibles

Saint Venant PDF

Cargado por

Información del documento

Título original

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Saint Venant PDF

Cargado por

Copyright:

Formatos disponibles

Tecnura

Disponible en: http://www.redalyc.org/articulo.oa?id=257040047007

Cmo citar el artculo

Comportamiento de las ecuaciones de Saint-Venant en 1D y aproximaciones

Fecha de recepcin: 6 de septiembre de 2014 Fecha de aceptacin: 6 de abril de 2015

RESUMEN simple que produzca soluciones tiles a las necesi-

INTRODUCCIN inestables debido a la alta acumulacin de erro-

Modelo de onda cinemtica una de continuidad y otra de momentum. La solu-

Figura 1. Seccin Transversal del Canal de Irrigacin.

Fuente: elaboracin propia.

suficientemente baja, simulando las condiciones la simulacin de flujo horizontal correspondiente

A continuacin se presentan las condiciones de Caudal inicial: 10 m3/s.

Figura 2. Condicin de frontera aguas arriba, simulacin de flujo horizontal.

Fuente: elaboracin propia.

Figura 3. Simulacin de flujo horizontal.

Fuente: elaboracin propia.

Caso 2: flujo subcrtico Limnigrama de niveles aguas arriba.

Figura 4. Condicin de frontera aguas arriba, simulacin de flujo subcrtico.

Fuente: elaboracin propia.

Figura 5. Condicin de frontera aguas arriba, simulacin de flujo subcrtico.

Fuente: elaboracin propia.

Caso 3: flujo supercrtico Hidrograma de caudales aguas abajo.

Figura 6. Condicin de frontera aguas arriba, simulacin de flujo supercrtico.

Fuente: elaboracin propia

Figura 7. Condicin de frontera aguas abajo, simulacin de flujo supercrtico.

Fuente: elaboracin propia.

RESULTADOS la entrada y a las 5 horas el caudal en la salida es de

Figura 8. Hidrogramas y limnigramas de entrada y salida para simulacin de flujo horizontal.

Fuente: elaboracin propia.

Figura 9. Resultados de la simulacin de flujo horizontal con modelo de onda dinmica.

Fuente: elaboracin propia.

Caso 2 El hidrograma de salida alcanza el mismo caudal

Fuente: elaboracin propia.

Fuente: elaboracin propia.

Caso 3 largo del canal. El hidrograma de salida alcanza

Fuente: elaboracin propia.

Fuente: elaboracin propia.

REFERENCIAS Macchione, F., & Morelli, M. (2003). Pactical aspect

También podría gustarte