Skip to main content

Derive kinematic wave approximation

Short answer

Neglect local acceleration, convective acceleration and pressure gradient terms in the momentum equation, giving S_f = S_0 and locally uniform flow Q = α A^(5/3). Substituting this into continuity yields the kinematic wave equation ∂A/∂t + c ∂A/∂x = q with celerity c = dQ/dA ≈ (5/3)V, describing flood wave propagation without attenuation.

Derive the kinematic wave approximation for the movement of fluid.

Starting point

The momentum equation of the Saint-Venant system is

∂Q∂t⏟1+∂∂x(Q2A)⏟2+gA∂y∂x⏟3−gA(S0−Sf)=0\underbrace{\frac{\partial Q}{\partial t}}_{1}+\underbrace{\frac{\partial}{\partial x}\left(\frac{Q^2}{A}\right)}_{2}+\underbrace{gA\frac{\partial y}{\partial x}}_{3}-gA(S_0-S_f)=0

Dividing by gAgA gives the form with slopes: Sf=S0−∂y∂x−Vg∂V∂x−1g∂V∂tS_f=S_0-\dfrac{\partial y}{\partial x}-\dfrac{V}{g}\dfrac{\partial V}{\partial x}-\dfrac1g\dfrac{\partial V}{\partial t}, in which terms 1, 2 and 3 are the local acceleration, convective acceleration and pressure terms.

Approximation

For flood waves on steep channels these three terms are very small compared with S0S_0 and SfS_f. Neglecting them:

S0−Sf=0 ⇒ Sf=S0S_0-S_f=0\ \Rightarrow\ S_f=S_0

The flow is then locally uniform: the discharge depends only on the area at the same section. Using Manning's equation Q=1nA5/3P2/3S01/2Q=\dfrac{1}{n}\dfrac{A^{5/3}}{P^{2/3}}S_0^{1/2},

Q=αkAβ,β=53,  αk=S01/2nP2/3Q=\alpha_kA^{\beta},\qquad \beta=\frac53,\ \ \alpha_k=\frac{S_0^{1/2}}{nP^{2/3}}

Kinematic wave equation

Continuity: ∂A∂t+∂Q∂x=q\dfrac{\partial A}{\partial t}+\dfrac{\partial Q}{\partial x}=q. From the relation above, ∂Q∂x=dQdA∂A∂x\dfrac{\partial Q}{\partial x}=\dfrac{dQ}{dA}\dfrac{\partial A}{\partial x}, hence

∂A∂t+c∂A∂x=q,or∂Q∂t+c∂Q∂x=c q,c=dQdA=αkβAβ−1=βV\frac{\partial A}{\partial t}+c\frac{\partial A}{\partial x}=q,\qquad\text{or}\qquad \frac{\partial Q}{\partial t}+c\frac{\partial Q}{\partial x}=c\,q,\qquad c=\frac{dQ}{dA}=\alpha_k\beta A^{\beta-1}=\beta V

A flood wave therefore moves downstream at the kinematic celerity c≈53Vc\approx\frac53V without changing shape (no attenuation or backwater effects).

✓ Answer checked by Claude Opus 5.5 (high effort)

All Computational Techniques in Civil Engineering chapter 5 questions →