Computational Techniques in Water Resources ENCE 370 Chapter 3: Method of Characteristics

Fundamentals of the Method of Characteristics in Water Resources

In hydraulic engineering, fluid transients caused by rapid valve closures, pump trips, or sudden surge operations generate high-frequency pressure waves known as water hammer. Similarly, dam breaks and gate operations initiate unsteady wave propagation in open channels. Direct numerical integration of hyperbolic partial differential equations describing these events can introduce numerical dispersion and instability. The Method of Characteristics solves this challenge by reformulating partial derivatives into total derivatives along pathlines called characteristic curves.

1. Concept, Solution Requirements, Advantages, and Limitations of the Method of Characteristics

The core concept behind the Method of Characteristics involves finding trajectories along the space-time ($x-t$) plane along which information about physical variables (such as piezometric head $H$ and discharge $Q$, or water depth $y$ and velocity $v$) propagates at specific characteristic wave speeds ($a$ or $c$).

A. Solution Requirements

• The governing partial differential equations must be strictly hyperbolic in nature.
• Accurate physical wave speed propagation velocities (e.g., speed of sound in water through elastic pipes $a$, or shallow water wave celerity $c = \sqrt{gy}$) must be determined.
• Proper spatial discretization ($\Delta x$) and temporal discretization ($\Delta t$) satisfying the Courant criterion must be chosen.
• Well-defined boundary conditions at all physical terminals (reservoirs, valves, junctions, dead ends) must be established.

B. Advantages of the Method of Characteristics

• High Accuracy: Accurately tracks steep wave fronts and sharp pressure peaks without introducing excessive numerical smoothing or artificial numerical diffusion.
• Physical Clarity: Preserves the physical concept of wave propagation along directional lines ($C^+$ and $C^-$).
• Flexibility with Boundary Conditions: Simplifies complex boundary evaluations by isolating upstream and downstream wave contributions.

C. Limitations of the Method of Characteristics

• Restricted primarily to 1D and simple 2D hyperbolic transient systems.
• Requires strict adherence to grid spacing constraints; handling variable pipe diameters or complex networks requires spatial interpolation, which may introduce interpolation errors.

2. Governing Equations and Characteristics Formulation for Pipe Flow

Unsteady flow in elastic closed conduits is governed by 1D hyperbolic continuity and momentum equations:

Where $H$ is the piezometric head, $Q$ is discharge, $a$ is wave speed, $g$ is acceleration due to gravity, $D$ is pipe diameter, $A$ is cross-sectional area, and $f$ is the Darcy-Weisbach friction factor.

By applying linear combination multipliers ($\lambda$) using the Method of Characteristics, these equations transform into two distinct sets of ordinary differential equations valid along the positive ($C^+$) and negative ($C^-$) characteristic lines:

3. Finite Difference Solution and Rectangular Grid Algorithm for Water Hammer

To implement the Method of Characteristics on a digital computer, the $x-t$ continuum is discretized into a rectangular grid. Let node $P$ be an unknown interior grid point at time $t + \Delta t$, calculated from adjacent known grid points $A$ (upstream) and $B$ (downstream) at time $t$.

Integrating the differential equations along $C^+$ ($AP$) and $C^-$ ($BP$) using first-order finite differences yields algebraic compatibility equations:

Where the constants $C_P$, $C_M$, and impedance factor $B_P$ are defined as:

Solving simultaneously for an interior node $P$ gives the explicit solutions for head ($H_P$) and discharge ($Q_P$):

4. Stability Analysis and Courant Condition in the Method of Characteristics

The stability of the numerical solution generated by the Method of Characteristics requires that the numerical domain of dependence completely encloses the physical domain of dependence. This requires satisfying the **Courant-Friedrichs-Lewy (CFL)** condition:

When $C_r = 1$, characteristic lines pass exactly through grid intersections ($A$ and $B$). If $C_r < 1$, spatial grid interpolations are necessary, which can introduce numerical damping. If $C_r > 1$, the numerical solution becomes unstable and experiences catastrophic error growth.

5. Boundary Conditions for Unsteady Pipe Flow Problems

At boundaries (e.g., reservoirs, downstream control valves, series pipe junctions), only one characteristic line enters from the internal grid domain. Thus, an auxiliary boundary relation is required:

• Constant Head Reservoir Boundary (Upstream): $H_P = H_{res}$. The unknown discharge $Q_P$ is computed directly using the incoming $C^-$ equation: $Q_P = \frac{H_{res} – C_M}{B_P}$.
• Downstream Valve Boundary: Governed by the orifice equation $Q_P = \tau \cdot Q_0 \sqrt{\frac{H_P}{H_0}}$, where $\tau$ is the non-dimensional valve opening ratio. Combining this with the incoming $C^+$ equation yields a quadratic equation in $Q_P$.
• Dead End / Closed Valve: $Q_P = 0$, giving $H_P = C_P$.

6. Method of Characteristics for 1D Unsteady Open Channel Flow

For open channel flow governed by the full 1D Saint-Venant equations, the characteristic directions are defined by the flow velocity $v$ and shallow water celerity $c = \sqrt{gy}$:

Integrating the compatibility equations along these curved $C^+$ and $C^-$ lines on an $x-t$ grid yields the time-marching solution for water depth $y_P$ and velocity $v_P$, making the Method of Characteristics a versatile tool for routing dam-break waves and hydrographs in rivers.