One-dimensional Saint-Venant equation

In fluid dynamics, the one-dimensional (1-D) Saint-Venant equation was derived by Adhémar Jean Claude Barré de Saint-Venant, and is commonly used to model transient open-channel flow and surface runoff. It is a simplification of the two-dimensional (2-D) shallow water equations, which are also known as the two-dimensional Saint-Venant equations. The 1-D simplification is used exclusively in models including HEC-RAS, SWMM5, InfoWorks, ISIS, Flood Modeller, MIKE 11, and MIKE SHE because it is significantly easier to solve than the full shallow water equations. Common applications of the 1-D Saint-Venant Equation include dam break analyses, storm pulses in an open channel, as well as storm runoff in overland flow.

Derivation from Navier–Stokes equations

The 1-D Saint-Venant equation can be derived from the Navier–Stokes equations that describe fluid motion. The x-component of the Navier–Stokes equations – when expressed in Cartesian coordinates in the x-direction – can be written as:

\frac{\partial u}{\partial t} + u \frac{\partial u}{\partial x} + v \frac{\partial u}{\partial y}+ w \frac{\partial u}{\partial z}= -\frac{\partial p}{\partial x} \frac{1}{\rho} + \nu \left(\frac{\partial^2 u}{\partial x^2} + \frac{\partial^2 u}{\partial y^2} + \frac{\partial^2 u}{\partial z^2}\right)+ f_x,

where u is the velocity in the x-direction, v is the velocity in the y-direction, w is the velocity in the z-direction, t is time, p is the pressure, ρ is the density of water, ν is the kinematic viscosity, and f_x is the body force in the x-direction.

1.	If it is assumed that friction is taken into account as a body force, then ν can be assumed as zero so: $\nu \left(\frac{\partial^2 u}{\partial x^2} + \frac{\partial^2 u}{\partial y^2} + \frac{\partial^2 u}{\partial z^2}\right)= 0.$
2.	Assuming one-dimensional flow in the x-direction it follows that: $v\frac{\partial u}{\partial y}+ w \frac{\partial u}{\partial z} = 0$
3.	Assuming also that the pressure distribution is approximately hydrostatic it follows that: $p = \rho g h$ or in differential form: $\partial p = \rho g \left(\partial h \right).$ And when these assumptions are applied to the x-component of the Navier–Stokes equations: $-\frac{\partial p}{\partial x} \frac{1}{\rho} = -\frac{1}{\rho}\frac{\rho g \left(\partial h \right)}{\partial x} = -g \frac{\partial h}{\partial x}.$
4.	There are 2 body forces acting on the channel fluid, gravity, and friction: $f_x =f_{x,g} + f_{x,f}$ where f_x,g is the body force due to gravity and f_x,f is the body force due to friction.
5.	f_x,g can be calculated using basic physics and trigonometry: $F_{g} = (\sin\theta) gM$ where F_g is the force of gravity in the x-direction, θ is the angle, and M is the mass. Figure 1: Diagram of block moving down an inclined plane. The expression for sin θ can be simplified using trigonometry as: $\sin\theta = \frac{opp}{hyp}.$ For small θ (reasonable for almost all streams) it can be assumed that: $\sin\theta = \tan\theta = \frac{opp}{adj} = S$ and given that f_x represents a force per unit mass, the expression becomes: $f_{x,g} = gS.$
6.	Assuming the energy grade line is not the same as the channel slope, and for a reach of consistent slope there is a consistent friction loss, it follows that: $f_{x,f} = S_f g.$
7.	All of these assumptions combined arrives at the 1-dimensional Saint-Venant equation in the x-direction: $\frac{\partial u}{\partial t} + u \frac{\partial u}{\partial x} + g \frac{\partial h}{\partial x} + g (S - S_f) = 0$ $(a)\quad \ \ (b)\quad \ \ \ (c) \qquad (d) \ \ \ (e)$ where (a) is the local acceleration term, (b) is the convective acceleration term, (c) is the pressure gradient term, (d) is the gravity term, and (e) is the friction term.^[1]

Terms

The local acceleration (a) can also be thought of as the “unsteady term” as this describes some change in velocity over time. The convective acceleration (b) is an acceleration caused by some change in velocity over position, for example the speeding up or slowing down of a fluid entering a constriction or an opening, respectively. Both these terms make up the inertia terms of the 1-dimensional Saint-Venant equation.

The pressure gradient term (c) describes how pressure changes with position, and since the pressure is assumed hydrostatic, this is the change in head over position. The gravity term (d) is the acceleration due to slope, while the friction term (e) accounts for losses in energy due to friction

Common simplifications

Dynamic wave

The dynamic wave is the term used to describe the full 1-dimensional Saint-Venant equation. It is numerically challenging to solve, but is valid for all channel flow scenarios. The dynamic wave is used for modeling transient storms in modeling programs including HEC-RAS,^[2] InfoWorks,^[3] MIKE 11,^[4] Wash 123d ^[5] and SWMM.

Kinematic wave

For the kinematic wave it is assumed that the flow is uniform, and that the friction slope is approximately equal slope of the channel. This simplifies the full Saint-Venant equation to the kinematic wave:

(S - S_f) = 0.

The kinematic wave is valid when the change in wave height over distance and velocity over distance and time is negligible relative to the bed slope, e.g. for shallow flows over steep slopes.^[6] The kinematic wave is used in HEC-HMS.^[7]

Diffusive wave

For the diffusive wave it is assumed that the inertial terms are less than the gravity, friction, and pressure terms. The diffusive wave can therefore be more accurately described as a non-inertia wave, and is written as:

g \frac{\partial h}{\partial x} + g (S - S_f) = 0.

The diffusive wave is valid when the inertial acceleration is much smaller than all other forms of acceleration, or in other words when there is primarily subcritical flow. Models that use the diffusive wave assumption include MIKE SHE^[8] and LISFLOOD-FP.^[9]

Example

Following is an example solution of the kinematic wave using an explicit finite difference algorithm.

Setup

Multiplying the kinematic wave equation by the area of the channel, it follows that:

\frac{\partial Q}{\partial t} + V_w \frac{\partial Q}{\partial x} = 0

where Q is the channel flow, and V_w is the velocity of wave propagation.

V_w is the celerity, and is equal to the derivative of the discharge with respect to the depth, normalized by the channel topwidth (T):

V_w = \frac{\frac{dQ}{dh}}{T} = \frac{1}{T} \frac{d}{dh} \left(\frac{1.0}{n}AR^{2/3}S^{1/2} \right)

where n is Manning’s roughness coefficient, and R is the hydraulic radius.

Assuming a wide channel, the hydraulic radius is approximately the same as the channel depth, simplifying the equation above to:

V_w = \frac{5}{3}V.

To describe the stability of the finite difference scheme the Courant number (C) is defined as:

C = V_w \frac{\Delta t}{\Delta x}

where Δt is the time discretization, and Δx is the space discretization.

Solution

Table 1: Flow prediction (Q, m³/s) using an explicit finite difference formulation of the kinematic wave with a Courant number of 0.9.

Suppose V_w is 1 m/s, the flow of interest is at 200 m intervals over a 1000 m reach, the initial flow for the entire reach is 3 m³/s, and an estimate is used for the hydrograph at location "0". Then the kinematic wave can be used to estimate the flow over time. If C is less than 1 there is some numerical dispersion but the solution is stable, but if C is greater than 1 the solution is numerically unstable. The case of C = 0.9 yields a Δt of 180 seconds, or 0.05 hours.

The notation used for solving the kinematic wave is:

Q^{j+1}_i

where $Q^j_i$ is the flow at location i at time j. Note that the j-term is not an exponent.

An explicit solution algorithm is used, with all of the flow values known for the initial conditions (t = 0), and at the upstream boundary (x = 0) for all values of t. The solution is explicit because known flows upstream at the previous time step (i-1, j) and at the same location at the previous timestep (i,j) are used to calculate the flow for the next timestep (j+1).

If the forward difference approximation is used then:

Q^{j+1}_i = \left(1-C\right) Q^j_i + C\left(Q^j_{i-1}\right).

Using this methodology a storm can be routed through the system, as displayed in Table 1.

The hydrographs at the upstream boundary and the downstream boundary of the routing channel are compared in figure 3.

To illustrate how the forward difference algorithm works, the flow at the location 400 m downstream of the upstream boundary after 0.5 hours will be illustrated (see highlighted entries in Table 1).

From table 1 it follows that:

Q^j_i = 9.6 \text{ m}^3/\text{s}

Q^j_{i-1} = 10.3 \text{ m}^3/\text{s}

Therefore:

Q^{j+1}_i = \left(1-0.9\right) 9.6 + 0.9\left(10.3\right) = 10.2 \text{ m}^3/\text{s}.

Figure 3: Flow at the upstream (x = 0 m) end and downstream (x = 1000 m) end of the example routing channel.

Courant number analysis

To demonstrate the effect of the Courant number (C) on the predicted hydrograph, the routing is repeated for various C values and the hydrograph at the downstream boundary ( x = 1000 m) is displayed. The results for a C of 0.7, 0.8, 0.9 and 1.0 are displayed in figure 4, and the results for a C of 1.0, 1.1, 1.2, and 1.3 are displayed in figure 5. When using an explicit algorithm, if C is less than 1 numerical dispersion occurs; however, the results are stable. This is illustrated by the hydrographs varying more as C decreases. When using an explicit algorithm, if C is greater than 1 numerical instabilities occur and the results become unstable. This is illustrated by the increasingly oscillatory behavior in the downstream hydrograph as C increases.

Figure 4: Flow at the downstream (x = 1000 m) end of the example routing channel for Courant Numbers of 1 or less (a), and a zoom in at the hydrograph peaks (b).

Figure 5: Flow at the downstream (x = 1000 m) end of the example routing channel for Courant numbers of 1 or greater (a), and a zoom in at the hydrograph peaks (b).

References

↑ Saint-Venant, A. (1871), Theorie du mouvement non permanent des eaux, avec application aux crues des rivieres et a l’introduction de marees dans leurs lits. Comptes rendus des seances de l’Academie des Sciences.
↑ Brunner, G. W. (1995), HEC-RAS River Analysis System. Hydraulic Reference Manual. Version 1.0 Rep., DTIC Document.
↑ Searby, D.; Dean, A.; Margetts J. (1998), Christchurch harbour Hydroworks modelling., Proceedings of the WAPUG Autumn meeting, Blackpool, UK.
↑ Havnø, K., M. Madsen, J. Dørge, and V. Singh (1995), MIKE 11-a generalized river modelling package, Computer models of watershed hydrology., 733–782.
↑ Yeh, G.; Cheng, J.; Lin, J.; Martin, W. (1995), A numerical model simulating water flow and contaminant and sediment transport in watershed systems of 1-D stream-river network, 2-D overland regime, and 3-D subsurface media . Computer models of watershed hydrology, 733–782.
↑ Novak, P., et al., Hydraulic Modelling – An Introduction: Principles, Methods and Applications. 2010: CRC Press.
↑ Scharffenberg, W. A., and M. J. Fleming (2006), Hydrologic Modeling System HEC-HMS: User's Manual, US Army Corps of Engineers, Hydrologic Engineering Center.
↑ DHI (Danish Hydraulic Institute) (2011), MIKE SHE User Manual Volume 2: Reference Guide, edited.
↑ Bates, P., T. Fewtrell, M. Trigg, and J. Neal (2008), LISFLOOD-FP user manual and technical note, code release 4.3. 6, University of Bristol.

Physical oceanography

Waves	Airy wave theory Ballantine scale Benjamin–Feir instability Boussinesq approximation Breaking wave Clapotis Cnoidal wave Cross sea Dispersion Edge wave Equatorial waves Fetch Gravity wave Infragravity wave Internal wave Kelvin wave Kinematic wave Longshore drift Luke's variational principle Mild-slope equation Radiation stress Rogue wave Rossby wave Rossby-gravity waves Sea state Seiche Significant wave height Sneaker wave Soliton Stokes boundary layer Stokes drift Stokes wave Swell Trochoidal wave Tsunami megatsunami Undertow Ursell number Wave action Wave base Wave height Wave power Wave radar Wave setup Wave shoaling Wave turbulence Wave–current interaction Waves and shallow water one-dimensional Saint-Venant equation shallow water equations Wind wave model

Circulation	Atmospheric circulation Baroclinity Boundary current Coriolis effect Downwelling Eddy Ekman layer Ekman spiral Ekman transport El Niño Southern Oscillation General Circulation Model Geostrophic current Global Ocean Data Analysis Project Gulf Stream Halothermal circulation Humboldt Current Hydrothermal circulation Langmuir circulation Longshore drift Loop Current Maelstrom Modular ocean model Ocean dynamics Ocean gyre Princeton ocean model Rip current Subsurface currents Sverdrup balance Thermohaline circulation shutdown Upwelling Whirlpool World Ocean Circulation Experiment

Tides	Amphidromic point Earth tide Head of tide Internal tide Lunitidal interval Perigean spring tide Rip tide Rule of twelfths Slack water Tidal bore Tidal force Tidal power Tidal race Tidal range Tidal resonance Tide gauge Tideline

Landforms	Abyssal fan Abyssal plain Atoll Bathymetric chart Coastal geography Cold seep Continental margin Continental rise Continental shelf Contourite Guyot Hydrography Oceanic basin Oceanic plateau Oceanic trench Passive margin Seabed Seamount Submarine canyon

Plate tectonics	Convergent boundary Divergent boundary Fracture zone Hydrothermal vent Marine geology Mid-ocean ridge Mohorovičić discontinuity Vine–Matthews–Morley hypothesis Oceanic crust Outer trench swell Ridge push Seafloor spreading Slab pull Slab suction Slab window Subduction Transform fault Volcanic arc

Ocean zones	Benthic Deep ocean water Deep sea Littoral Mesopelagic Oceanic Pelagic Photic Surf

Sea level	Deep-ocean Assessment and Reporting of Tsunamis Future sea level Global Sea Level Observing System North West Shelf Operational Oceanographic System Ocean Surface Topography Mission Sea-level curve Sea level rise World Geodetic System

Acoustics	Deep scattering layer Hydroacoustics Ocean acoustic tomography Sofar bomb SOFAR channel Underwater acoustics

Related	Argo Benthic lander Color of water DSV Alvin Marginal sea Marine energy Marine pollution Mooring National Oceanographic Data Center Ocean Ocean exploration Ocean observations Ocean reanalysis Ocean surface topography Ocean thermal energy conversion Oceanography Pelagic sediment Sea surface microlayer Sea surface temperature Seawater Science On a Sphere Thermocline Underwater glider Water column World Ocean Atlas

Category Commons

This article is issued from Wikipedia - version of the Wednesday, September 23, 2015. The text is available under the Creative Commons Attribution/Share Alike but additional terms may apply for the media files.