Ecuaciones de aguas poco profundas

Las ecuaciones de aguas poco profundas se derivan de las ecuaciones de conservación de la masa y conservación del momento lineal (las ecuaciones de Navier-Stokes), que se mantienen incluso cuando se rompen los supuestos de aguas poco profundas, como a través de un salto hidráulico. En el caso de un lecho horizontal, con unas fuerzas de Coriolis, de fricción y viscosidad despreciables, las ecuaciones de aguas someras son: $${\displaystyle {\begin{aligned}{\frac {\partial (\rho \eta )}{\partial t}}&+{\frac {\partial (\rho \eta u)}{\partial x}}+{\frac {\partial (\rho \eta v)}{\partial y}}=0,\\[3pt]{\frac {\partial (\rho \eta u)}{\partial t}}&+{\frac {\partial }{\partial x}}\left(\rho \eta u^{2}+{\frac {1}{2}}\rho g\eta ^{2}\right)+{\frac {\partial (\rho \eta uv)}{\partial y}}=0,\\[3pt]{\frac {\partial (\rho \eta v)}{\partial t}}&+{\frac {\partial (\rho \eta uv)}{\partial x}}+{\frac {\partial }{\partial y}}\left(\rho \eta v^{2}+{\frac {1}{2}}\rho g\eta ^{2}\right)=0.\end{aligned}}}$$ Aquí η es la altura total de la columna de fluido (profundidad instantánea del fluido en función de x, y y t), y el vector 2D (u, v) es la velocidad de flujo horizontal del fluido, promediada a través de la columna vertical. Además g es la aceleración debida a la gravedad y ρ es la densidad del fluido. La primera ecuación se deriva de la conservación de la masa, las dos segundas de la conservación del momento.^[3]

Expandiendo las derivadas en lo anterior usando la regla del producto, se obtiene la forma no conservativa de las ecuaciones de aguas someras. Como las velocidades no están sujetas a una ecuación de conservación fundamental, las formas no conservativas no se mantienen a través de un choque o salto hidráulico. También se incluyen los términos apropiados para las fuerzas de Coriolis, fricción y viscosidad, para obtener (para una densidad de fluido constante): $${\displaystyle {\begin{aligned}{\frac {\partial h}{\partial t}}&+{\frac {\partial }{\partial x}}{\Bigl (}(H+h)u{\Bigr )}+{\frac {\partial }{\partial y}}{\Bigl (}(H+h)v{\Bigr )}=0,\\[3pt]{\frac {\partial u}{\partial t}}&+u{\frac {\partial u}{\partial x}}+v{\frac {\partial u}{\partial y}}-fv=-g{\frac {\partial h}{\partial x}}-ku+\nu \left({\frac {\partial ^{2}u}{\partial x^{2}}}+{\frac {\partial ^{2}u}{\partial y^{2}}}\right),\\[3pt]{\frac {\partial v}{\partial t}}&+u{\frac {\partial v}{\partial x}}+v{\frac {\partial v}{\partial y}}+fu=-g{\frac {\partial h}{\partial y}}-kv+\nu \left({\frac {\partial ^{2}v}{\partial x^{2}}}+{\frac {\partial ^{2}v}{\partial y^{2}}}\right),\end{aligned}}}$$ donde:

u	es la velocidad en ladirección x, o velocidad zonal
v	es la velocidad en ladirección y , o velocidad meridional
H	es la altura media de la superficie de presión horizontal
h	es la desviación de la altura de la superficie de presión horizontal con respecto a su altura media, donde h: η(x, y, t) = H(x, y) + h(x, y, t)
b	es la altura topográfica desde una referencia D, donde b: H(x, y) = D + b(x,y)
g	es la aceleración debida a la gravedad
f	es el coeficiente de Coriolis asociado con la fuerza de Coriolis. En la Tierra, f es igual a 2Ω sin(φ), donde Ω es la velocidad de rotación angular de la Tierra (π/12 radians/hora), y φ es la latitud
k	es el coeficiente de arrastre por viscosidad
ν	es la viscosidad cinemática

A menudo ocurre que los términos cuadráticos en u y v, que representan el efecto de la advección a granel, son pequeños en comparación con los otros términos. Esto se llama equilibrio geostrófico, y equivale a decir que el número de Rossby es pequeño. Suponiendo también que la altura de la ola es muy pequeña comparada con la altura media (h ≪ H), se tiene (sin fuerzas viscosas laterales): $${\displaystyle {\begin{aligned}{\frac {\partial h}{\partial t}}&+H\left({\frac {\partial u}{\partial x}}+{\frac {\partial v}{\partial y}}\right)=0,\\[3pt]{\frac {\partial u}{\partial t}}&-fv=-g{\frac {\partial h}{\partial x}}-bu,\\[3pt]{\frac {\partial v}{\partial t}}&+fu=-g{\frac {\partial h}{\partial y}}-bv.\end{aligned}}}$$

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 del canal.

Las ecuaciones 1-D se utilizan ampliamente en modelos informáticos como TUFLOW, Mascaret (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.

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]

$${\displaystyle {\frac {\partial A}{\partial t}}+{\frac {\partial \left(Au\right)}{\partial x}}=0}$$

(1)

and

$${\displaystyle {\frac {\partial u}{\partial t}}+u\,{\frac {\partial u}{\partial x}}+g\,{\frac {\partial \zeta }{\partial x}}=-{\frac {P}{A}}\,{\frac {\tau }{\rho }},}$$

(2)

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 _b, el área de la sección transversal es: A = B (ζ - z_b) = B h. La profundidad instantánea del agua es h(x,t) = ζ(x,t) − z_b(x) con z_b(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: $${\displaystyle A(x,t)=\int _{0}^{h(x,t)}b(x,h')\,dh',}$$ 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 b(x, h) = B(x) 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_f(x, t) y el radio hidráulico R(x, t) se definen como: $${\displaystyle S=-{\frac {\mathrm {d} z_{\mathrm {b} }}{\mathrm {d} x}},}$$ $${\displaystyle S_{\mathrm {f} }={\frac {\tau }{\rho gR}}}$$ y $${\displaystyle R={\frac {A}{P}}.}$$

En consecuencia, la ecuación del momento (2) puede escribirse como:^[7]

$${\displaystyle {\frac {\partial u}{\partial t}}+u\,{\frac {\partial u}{\partial x}}+g\,{\frac {\partial h}{\partial x}}+g\,\left(S_{\mathrm {f} }-S\right)=0.}$$

( 3)

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 Q = Au:^[8]

$${\displaystyle {\frac {\partial Q}{\partial t}}+{\frac {\partial }{\partial x}}\left({\frac {Q^{2}}{A}}+g\,I_{1}\right)+g\,A\,\left(S_{f}-S\right)-g\,I_{2}=0,}$$

( 4)

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_b(x) (del punto más bajo en la sección transversal): $${\displaystyle {\begin{aligned}A(\sigma ,x)&=\int _{0}^{\sigma }B(\sigma ',x)\;\mathrm {d} \sigma ',\\I_{1}(\sigma ,x)&=\int _{0}^{\sigma }(\sigma -\sigma ')\,B(\sigma ^{\prime },x)\;\mathrm {d} \sigma '\qquad {\text{and}}\\I_{2}(\sigma ,x)&=\int _{0}^{\sigma }(\sigma -\sigma ')\,{\frac {\partial B(\sigma ',x)}{\partial x}}\;\mathrm {d} \sigma '.\end{aligned}}}$$

Arriba - en la ecuación del momento (4) en forma de conservación -A, I₁ y I₂ es evaluada como σ = h(x,t). El término g I₁ describe la fuerza hidrostática en una determinada sección transversal. Y, para un canal no prismático, g I ₂ 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.

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] $${\displaystyle {\frac {\mathrm {d} x}{\mathrm {d} t}}=u\pm c,}$$ with $${\displaystyle c={\sqrt {\frac {gA}{B}}}.}$$

El Número de Froude Fr= |u|/c determina si el flujo es subcrítico (Fr < 1) o supercrítico (Fr > 1).

Para un canal rectangular y prismático de anchura constante B, es decir, con A = B h y c = √gh los invariantes de Riemann son^[9] $${\displaystyle r_{+}=u+2{\sqrt {gh}}}$$ y $${\displaystyle r_{-}=u-2{\sqrt {gh}},}$$ por lo que las ecuaciones en forma característica son::^[9] $${\displaystyle {\begin{aligned}&{\frac {\mathrm {d} }{\mathrm {d} t}}\left(u+2{\sqrt {gh}}\right)=g\left(S-S_{f}\right)&&{\text{along}}\quad {\frac {\mathrm {d} x}{\mathrm {d} t}}=u+{\sqrt {gh}}\quad {\text{and}}\\&{\frac {\mathrm {d} }{\mathrm {d} t}}\left(u-2{\sqrt {gh}}\right)=g\left(S-S_{f}\right)&&{\text{along}}\quad {\frac {\mathrm {d} x}{\mathrm {d} t}}=u-{\sqrt {gh}}.\end{aligned}}}$$ 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]}

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 (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: $${\displaystyle g{\frac {\partial h}{\partial x}}+g(S_{f}-S)=0.}$$ 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: $${\displaystyle S_{f}-S=0.}$$ 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]