Ecuaciones unidimensionales de Saint-Venant
Contenido
Las ecuaciones unidimensionales (1-D) de Saint-Venant fueron derivadas por Adhémar Jean Claude Barré de Saint-Venant, y se utilizan comúnmente para modelar el flujo en canal abierto transitorio y la escorrentía superficial. Pueden considerarse como una contracción de las ecuaciones bidimensionales (2-D) de aguas poco profundas, que también se conocen como ecuaciones bidimensionales de Saint-Venant. Las ecuaciones de Saint-Venant 1-D contienen hasta cierto punto las principales características de la forma de la sección transversal "Sección (geometría)") del canal.
Las ecuaciones 1-D se utilizan ampliamente en modelos informáticos como TUFLOW"), Mascaret&action=edit&redlink=1 "MASCARET (Software) (aún no redactado)") (EDF), lang=es SIC (Irstea), HEC-RAS,[5] SWMM5, ISIS,[5] InfoWorks,[5] Flood Modeller, SOBEK 1DFlow MIKE 11"),[5] y MIKE SHE") porque son significativamente más fáciles de resolver que las ecuaciones completas de aguas poco profundas. Las aplicaciones comunes de las ecuaciones 1-D de Saint-Venant incluyen enrutamiento de inundaciones a lo largo de los ríos (incluyendo la evaluación de las medidas para reducir los riesgos de inundación), el análisis de la rotura de presas, los pulsos de tormenta en un canal abierto, así como la escorrentía de la tormenta en el flujo terrestre.
Ecuaciones
El sistema de ecuaciones diferenciales parciales que describen el flujo incompresible 1-D en un canal abierto de sección transversal arbitraria -tal como fue derivado y planteado por Saint-Venant en su artículo de 1871 (ecuaciones 19 y 20)- es:[6].
and.
donde x es la coordenada espacial a lo largo del eje del canal, t denota el tiempo, A(x,t) es el área de la sección transversal del flujo en la ubicación x, u(x,t) es la velocidad de flujo, ζ(x,t) es la elevación de la superficie libre y τ(x,t) es el esfuerzo cortante de la pared a lo largo del perímetro mojado P(x,t) de la sección transversal en x. Además ρ es la densidad del fluido (constante) y g es la aceleración gravitatoria.
La solución del sistema de ecuaciones hiperbólico (1)–(2) se obtiene a partir de la geometría de las secciones transversales, proporcionando una relación funcional entre el área de la sección transversal A y la elevación de la superficie ζ en cada posición x. Por ejemplo,Por ejemplo, para una sección transversal rectangular, con anchura de canal constante B y elevación del lecho del canal z , el área de la sección transversal es: . La profundidad instantánea del agua es con z(x) el nivel del lecho (es decir, la elevación del punto más bajo del lecho por encima del fecha, véase la figura de sección transversal). Para paredes de canales sin movimiento, el área de la sección transversal A en la ecuación (1) puede escribirse como:
con b(x,h) la anchura efectiva de la sección transversal del canal en la ubicación x cuando la profundidad del fluido es h - así que para canales rectangulares.[7].
El esfuerzo cortante de la pared τ depende de la velocidad del flujo u, pueden relacionarse utilizando, por ejemplo, la ecuación de Darcy-Weisbach, la fórmula de Manning o la fórmula de Chézy.
Además, la ecuación (1) es la ecuación de continuidad, que expresa la conservación del volumen de agua para este fluido homogéneo incompresible. La ecuación (2) es la ecuación de momentum, que da el equilibrio entre las fuerzas y las tasas de cambio del momentum.
La pendiente del lecho S(x), la pendiente de fricción S(x, t) y el radio hidráulico R(x, t) se definen como:
y.
En consecuencia, la ecuación del momento (2) puede escribirse como:[7].
Conservación del momento
La ecuación del momento (3) también puede ser planteada en la llamada forma de conservación"), o forma euleriana, mediante algunas manipulaciones algebraicas sobre las ecuaciones de Saint-Venant, (1) y (3). En términos del caudal "Caudal (fluido)") :[8].
donde A, I y I son funciones de la geometría del canal, descrito en términos de la anchura del canal B(σ,x). Aquí, σ es la altura sobre el punto más bajo de la sección transversal en la ubicación x. Entonces σ es la altura sobre el nivel del lecho z(x) (del punto más bajo en la sección transversal):.
Arriba - en la ecuación del momento (4) en forma de conservación -A, I y I es evaluada como . El término describe la fuerza hidrostática en una determinada sección transversal. Y, para un canal no prismático "Prisma (geometría)"), da los efectos de las variaciones de la geometría a lo largo del eje x del canal.
En las aplicaciones, dependiendo del problema en cuestión, a menudo se prefiere utilizar la ecuación del momento en forma no conservativa, (2) or (3), o la forma de conservación (4). Por ejemplo, en el caso de la descripción de los saltos hidráulicos, se prefiere la forma de conservación ya que el flujo de momento es continuo a través del salto.
Características
Las ecuaciones de Saint-Venant (1)–(2) puede analizarse mediante el método de características").[9][10][11][12] Las dos aceleraciones dx/dt en las curvas características son:[8]
with.
El Número de Froude Fr= |u|/c determina si el flujo es subcrítico () o supercrítico ().
Para un canal rectangular y prismático de anchura constante B, es decir, con y los invariantes de Riemann") son[9] y
por lo que las ecuaciones en forma característica son::[9]
Los invariantes de Riemann y el método de las características para un canal prismático de sección transversal arbitraria son descritos por Didenkulova y Pelinovsky (2011).[12].
Las características y los invariantes de Riemann proporcionan información importante sobre el comportamiento del flujo, además de que pueden utilizarse en el proceso de obtención de soluciones (analíticas o numéricas).[13][14][15][16].
Modelación derivada
La onda dinámica es la ecuación unidimensional completa de Saint-Venant. Es numéricamente difícil de resolver, pero es válida para todos los escenarios de flujo del canal. La onda dinámica se utiliza para modelar tormentas transitorias en programas de modelización como Mascaret&action=edit&redlink=1 "MASCARET (Software) (aún no redactado)") (EDF), SIC (Irstea), HEC-RAS,[17] InfoWorks_ICM Archivado el 25 de octubre de 2016 en Wayback Machine.,[18] MIKE 11"),[19] Wash 123d[20] y SWMM5").
En el orden de las simplificaciones crecientes, al eliminar algunos términos de las ecuaciones completas de Saint-Venant 1D (también conocida como ecuación de onda dinámica), obtenemos la también clásica ecuación de onda difusiva y la ecuación de onda cinemática.
Para la onda difusiva se supone que los términos de inercia son menores que los términos de gravedad, fricción y presión. Por lo tanto, la onda difusiva puede describirse con mayor precisión como una onda no inercial, y se escribe como:
La onda difusiva es válida cuando la aceleración inercial es mucho menor que todas las demás formas de aceleración o, en otras palabras, cuando hay principalmente un flujo subcrítico, con valores de Froude bajos. Los modelos que utilizan la hipótesis de la onda difusiva incluyen MIKE SHE")[21] y LISFLOOD-FP.[22] In the software SIC (Irstea) estas opciones también están disponibles, ya que los 2 términos de inercia (o cualquiera de ellos) se pueden eliminar en opción desde la interfaz.
Para la onda cinemática se supone que el flujo es uniforme, y que la pendiente de fricción es aproximadamente igual a la pendiente del canal. Esto simplifica la ecuación completa de Saint-Venant a la onda cinemática:
La onda cinemática es válida cuando el cambio en la altura de la onda a lo largo de la distancia y la velocidad a lo largo de la distancia y el tiempo es insignificante en relación con la pendiente del lecho, por ejemplo, para flujos poco profundos sobre pendientes pronunciadas.[23] La onda cinemática se usa en HEC-HMS").[24].