Skip to main content

First-order kinematic wave, flood peak at community

Short answer

Apply the first-order implicit scheme Q_{i+1}^{n+1} = (Q_{i+1}^n + C Q_i^{n+1})/(1+C) with c = 1.2 m/s, Δt = 1 h, Δx = 7.2 km, giving C = 0.6. The downstream peak discharge is (100 + 0.6×150)/1.6 = 118.75 m³/s.

A flood of 150 m3^3/s peak discharge passed a gauging station at 12:00 noon on a river. There is a community adjacent to the river 7.2 km downstream. What will be the value of peak discharge at that community at 12:00 noon if the velocity of flow is 1.2 m/s2^2 [as printed] and the peak discharge at that community at 9:00 A.M. is 100 m3^3/s? Assume width of river as [?] and use the first order accurate numerical scheme of the kinematic wave equation. Take Δx=7.2\Delta x = 7.2 km and Δt=1\Delta t = 1 hr.

The first-order accurate (implicit) scheme of the linear kinematic wave equation ∂Q/∂t+c ∂Q/∂x=0\partial Q/\partial t+c\,\partial Q/\partial x=0 is

Qi+1n+1−Qi+1nΔt+c Qi+1n+1−Qin+1Δx=0 ⇒ Qi+1n+1=Qi+1n+C Qin+11+C,C=cΔtΔx\frac{Q_{i+1}^{n+1}-Q_{i+1}^n}{\Delta t}+c\,\frac{Q_{i+1}^{n+1}-Q_i^{n+1}}{\Delta x}=0 \ \Rightarrow\ Q_{i+1}^{n+1}=\frac{Q_{i+1}^n+C\,Q_i^{n+1}}{1+C},\qquad C=\frac{c\Delta t}{\Delta x}

Interpretation of the data

  • The gauging station is point ii; the community (Δx=7.2\Delta x=7.2 km downstream) is i+1i+1.
  • Qin+1=150Q_i^{n+1}=150 m3^3/s (peak at the station at 12:00 noon).
  • Qi+1n=100Q_{i+1}^n=100 m3^3/s is the flow already at the community at the earlier level (taken as the 9:00 A.M. value, with Δt=1\Delta t=1 h per step).
  • The velocity given as 1.2 (m/s) is used as the wave celerity cc; the channel width is then not needed, because the equation is written in terms of QQ.

Calculation

C=c ΔtΔx=1.2×36007200=0.6C=\frac{c\,\Delta t}{\Delta x}=\frac{1.2\times3600}{7200}=0.6 Qi+1n+1=100+0.6×1501+0.6=1901.6=118.75 m3/sQ_{i+1}^{n+1}=\frac{100+0.6\times150}{1+0.6}=\frac{190}{1.6}=118.75\ \text{m}^3/\text{s}

Answer: the discharge at the community at 12:00 noon is about 118.75 m3^3/s (the peak is delayed and attenuated numerically). The data are incomplete in the question (width missing, unit of velocity printed as m/s2^2), so this reading is assumed.

✓ Answer checked by Claude Opus 5.5 (high effort)

All Computational Techniques in Civil Engineering chapter 5 questions →