Surface flow formula (Water Overlay)

The printable version is no longer supported and may have rendering errors. Please update your browser bookmarks and please use the default browser print function instead.

Imbalances in water levels across the grid drive the flow of water until a state of equilibrium is reached in terms of h (the height of the water column) and flux. Behavior of the flow is described by a second-order semi-discrete central-upwind scheme produced by Kurganov and Petrova (2007)^[1], which is based on the 2-D Saint-Venant equations (a.k.a. shallow water equations):

{\begin{aligned}{\frac {\partial h}{\partial t}}&+{\frac {\partial (hu)}{\partial x}}+{\frac {\partial (hv)}{\partial y}}=0,\\[3pt]{\frac {\partial (hu)}{\partial t}}&+{\frac {\partial }{\partial x}}\left(hu^{2}+{\frac {1}{2}}gh^{2}\right)+{\frac {\partial (huv)}{\partial y}}=-gh{\frac {\partial B}{\partial x}}-ghn^{2}u{\sqrt {u^{2}+v^{2}}}h^{-{\frac {4}{3}}},\\[3pt]{\frac {\partial (hv)}{\partial t}}&+{\frac {\partial (huv)}{\partial x}}+{\frac {\partial }{\partial y}}\left(hv^{2}+{\frac {1}{2}}gh^{2}\right)=-gh{\frac {\partial B}{\partial y}}-ghn^{2}v{\sqrt {u^{2}+v^{2}}}h^{-{\frac {4}{3}}},\end{aligned}}

where

u	is the velocity in the x-direction
v	is the velocity in the y-direction
h	is the water depth
B	is the bottom elevation
g	is the acceleration due to gravity
n	is the Gauckler–Manning coefficient

References

↑ Kurganov A, Petrova G (2007) ∙ A Second-Order Well-Balanced Positivity Preserving Central-Upwind Scheme for the Saint-Venant System ∙ found at: http://www.math.tamu.edu/~gpetrova/KPSV.pdf (last visited 2019-04-11)