ENGINEERING HYDROLOGY:  CHAPTER 096 - INTRODUCTION TO DYNAMIC WAVES



1. RATIONALE


1.01
A kinematic wave simplifies the momentum equation to a statement of steady uniform flow.


1.02
A diffusion wave simplifies the momentum equation to a statement of steady nonuniform flow.


1.03
A third type of open channel flow wave, the dynamic wave, considers the complete momentum equation, including its inertial components.


1.04
As such, the dynamic wave contains more physical information that either kinematic or diffusion waves.


1.05
However, dynamic wave solutions are more complex.


1.06
In a dynamic wave equation, the equations of mass and momentum conservation are solved by a numerical procedure, typically using finite differences.


1.07
The partial differential equations governing the motion are discretized following a chosen numerical scheme.


1.08
Among the various schemes that have been used for this purpose, the Preissmann scheme is perhaps the most popular.


1.09
This is a four-point scheme, centered in the temporal derivatives, and slightly off-centered in the spatial derivatives.


1.10
The off-centering of the spatial derivatives produces a small amount of numerical diffusion, which is necessary to control the numerical stability of the nonlinear scheme.


1.11
This produces a workable yet sufficiently accurate scheme.



2. METHODOLOGY


2.01
The stream channel reach is divided into several subreaches for computational purposes.


2.02
The application of the Preissmann scheme to the governing equations for the various subreaches results in a matrix solution requiring a double-sweep algorithm.


2.03
A double-sweep algorithm is a type of Gauss elimination technique which accounts only for the nonzero entries of the coefficient matrix, which are located within a narrow band surrounding the main diagonal.


2.04
The upstream boundary condition is a discharge hydrograph.


2.05
The downstream boundary condition is a stage hydrograph.


2.06
The solution of the set of hyperbolic equations marches in time until a specified number of time intervals is completed.


2.07
In practice, a dynamic wave solution represents an order-of-magnitude increase in complexity and associated data requirements, when compared to either kinematic or diffusion wave solutions.


2.08
Its use is recommended in situations where neither kinematic or diffusion waves are likely to adequately represent the physical phenomena.


2.09
In particular, dynamic wave solutions are applicable to flow over very mild slopes, flow into large reservoirs, strong backwater conditions, and flow reversals.


2.10
In general, the dynamic wave solution is recommended for cases warranting a precise calculation of the unsteady variation of river stages.



3. RELEVANCE TO ENGINEERING HYDROLOGY


3.01
Dynamic wave solutions are often referred to as hydraulic river routing.


3.02
As such, they have the capability to calculate unsteady discharges and stages when presented with the appropriate geometric channel data and initial and boundary conditions.


3.04
Kinematic wave solutions calculate unsteady discharges; the corresponding stages are subsequently calculated from the appropriate rating curve.


3.05
An equilibrium, that is, a steady uniform, rating curve is used for this purpose.


3.06
Diffusion wave may or may not use equilibrium ratings to calculate stages.


3.07
Dynamic waves rely on the physics of the phenomena, as built into the governing equations, to generate their own unsteady rating curves.


3.08
A looped rating curve is produced at every cross section, as shown here.


3.09
For any given stage, the discharge is higher in the rising limb of the hydrograph, and lower in the receding limb.


3.10
This loop is due to hydrodynamic reasons, and should not be confused with other loops, which may be due to erosion, sedimentation, or changes in bed configuration.


3.11
The width of the loop is a measure of the flow unsteadiness, with wider loops corresponding to highly unsteady flows, that is, dynamic wave flow.


3.12
If the loop is narrow, it implies that the flow is mildly unsteady, perhaps a diffusion wave.


3.13
If the loop is practically nonexistent, the flow can be approximated as a kinematic wave.


3.14
In fact, the basic assumption of kinematic flow is that momentum can be simulated as steady flow, that is, that the rating curve is single valued, with no loop.


3.15
The preceding observations lead to the conclusion that the relevance of dynamic waves to engineering hydrology is directly related to the extent of flow unsteadiness and the size of the associated loop in the rating curve.


3.16
For highly unsteady flows such as dam-breach flood waves, it may well be the only way to account for the looped rating.


3.17
For other less unsteady flows, kinematic and diffusion waves remain a viable alternative, provided their applicability can be clearly demonstrated.


3.18
The applicability of kinematic and diffusion waves is a function of the product of hydrograph time-to-peak and bottom slope.


3.19
A source of complexity in dynamic wave solutions is the need to specify a downstream boundary condition, typically a discharge-stage rating.


3.20
This discharge-stage rating is usually specified as a single-valued rating, that is, a kinematic rating.


3.21
This contradicts the solution at the downstream boundary, because the dynamic wave solution generates its own looped ratings at the interior points.


3.22
A finely balanced model will usually go unstable in this situation, but models with heavy numerical diffusion may survive, only to give erroneous results.


3.23


Narrator: Victor M. Ponce

Music: Fernando Oñate

Editor: Flor Pérez


Copyright © 2011

Visualab Productions

All rights reserved