◀ Back to Note

Tony's Intertwingled Memex • Reference Manual

Open Drainage Channel Hydraulic Cross-Section Optimization (Environmental & Civil Engineering) Solution

Formulas & Definitions

A water resources and civil environmental engineer at the UNLV Howard R. Hughes College of Engineering is designing an open, rectangular reinforced-concrete flood control drainage channel for the Clark County Regional Flood Control District in the Las Vegas valley.

Hydrologic peak runoff calculations dictate that the channel must maintain a fixed cross-sectional water flow area of:

A(x,y)=xy=18 m2A(x, y) = xy = 18\text{ m}^2A(x,y)=xy=18 m2
where xxx is the bottom channel width (in meters, m\text{m}m) and yyy is the vertical water depth (in meters, m\text{m}m).

To minimize excavation costs, frictional head loss along the fluid boundary, and the expensive shotcrete concrete liner material across the bottom bed and both vertical sidewalls, the engineer must minimize the wetted perimeter:

P(x,y)=x+2yP(x, y) = x + 2yP(x,y)=x+2y

  1. Formulate the objective function P(x,y)P(x, y)P(x,y) and the constraint equation g(x,y)=kg(x, y) = kg(x,y)=k.
  2. Compute the gradient vectors ∇P(x,y)\nabla P(x, y)∇P(x,y) and ∇g(x,y)\nabla g(x, y)∇g(x,y).
  3. Formulate the Lagrange multiplier vector equation ∇P=λ∇g\nabla P = \lambda \nabla g∇P=λ∇g and equate the individual scalar components.
  4. Eliminate λ\lambdaλ to establish the optimal geometric aspect ratio between channel bottom width xxx and fluid depth yyy.
  5. Substitute into the flow area constraint to calculate the optimal channel dimensions (xopt,yopt)(x_{\text{opt}}, y_{\text{opt}})(xopt​,yopt​) and the minimum wetted perimeter Pmin⁡P_{\min}Pmin​.
  6. Verify why this optimal hydraulic ratio (x=2yx = 2yx=2y, where width equals twice the depth) represents the classic "most efficient rectangular open channel" in civil fluid mechanics.

Step-by-Step Solution:

  1. Formulate Functions:

    • Objective function to minimize: Wetted perimeter P(x,y)=x+2yP(x, y) = x + 2yP(x,y)=x+2y (base plus two vertical wet sidewalls).
    • Constraint equation: Fixed wet flow area g(x,y)=xy=18g(x, y) = xy = 18g(x,y)=xy=18, where x>0x > 0x>0 and y>0y > 0y>0.
  2. Compute Gradient Vectors: State the Two-Variable Gradient Formula (defined in The Gradient Vector):

    ∇P(x,y)=⟨∂P∂x,∂P∂y⟩=⟨1,2⟩=i^+2j^\nabla P(x, y) = \left\langle \frac{\partial P}{\partial x}, \frac{\partial P}{\partial y} \right\rangle = \langle 1, 2 \rangle = \hat{i} + 2\hat{j}∇P(x,y)=⟨∂x∂P​,∂y∂P​⟩=⟨1,2⟩=i^+2j^​
    ∇g(x,y)=⟨∂g∂x,∂g∂y⟩=⟨y,x⟩=yi^+xj^\nabla g(x, y) = \left\langle \frac{\partial g}{\partial x}, \frac{\partial g}{\partial y} \right\rangle = \langle y, x \rangle = y\hat{i} + x\hat{j}∇g(x,y)=⟨∂x∂g​,∂y∂g​⟩=⟨y,x⟩=yi^+xj^​

  3. Formulate Lagrange Multiplier System: State the Lagrange Multiplier Vector System (defined in Lagrange Multipliers):

    ∇P=λ∇gandxy=18\nabla P = \lambda \nabla g \quad \text{and} \quad xy = 18∇P=λ∇gandxy=18
    Equate scalar components:
    i^:1=λy\hat{i}: 1 = \lambda yi^:1=λy
    j^:2=λx\hat{j}: 2 = \lambda xj^​:2=λx
    Constraint:xy=18\text{Constraint}: xy = 18Constraint:xy=18

  4. Eliminate λ\lambdaλ and Determine Geometric Aspect Ratio: Because physical dimensions require depth y>0y > 0y>0 and width x>0x > 0x>0, we can safely solve for λ\lambdaλ without dividing by zero: From i^\hat{i}i^:

    λ=1y\lambda = \frac{1}{y}λ=y1​
    Substitute λ=1y\lambda = \frac{1}{y}λ=y1​ into j^\hat{j}j^​:
    2=(1y)x  ⟹  2=xy  ⟹  x=2y2 = \left(\frac{1}{y}\right)x \implies 2 = \frac{x}{y} \implies x = 2y2=(y1​)x⟹2=yx​⟹x=2y
    The optimal bottom width must be exactly twice the water depth (x=2yx = 2yx=2y).

  5. Solve for Optimal Channel Dimensions: Substitute x=2yx = 2yx=2y into the flow area constraint xy=18xy = 18xy=18:

    (2y)y=18  ⟹  2y2=18  ⟹  y2=9  ⟹  y=3 m(2y)y = 18 \implies 2y^2 = 18 \implies y^2 = 9 \implies y = 3\text{ m}(2y)y=18⟹2y2=18⟹y2=9⟹y=3 m
    (Discarding negative root y=−3y = -3y=−3 since physical depth must be positive). Now compute the channel bottom width:
    x=2(3)=6 mx = 2(3) = 6\text{ m}x=2(3)=6 m
    The minimum wetted perimeter is:
    Pmin⁡=x+2y=6+2(3)=6+6=12 mP_{\min} = x + 2y = 6 + 2(3) = 6 + 6 = 12\text{ m}Pmin​=x+2y=6+2(3)=6+6=12 m

  6. Hydraulic Engineering Interpretation: In civil fluid mechanics, Manning's equation states that fluid discharge velocity increases with hydraulic radius Rh=A/PR_h = A / PRh​=A/P. For a fixed cross-sectional area AAA, minimizing the wetted perimeter PPP directly maximizes the hydraulic radius RhR_hRh​, maximizing volumetric conveyance capacity while minimizing frictional shear resistance against the rough concrete liner walls. The aspect ratio x=2yx = 2yx=2y (a half-square profile) is the universal theoretically optimal cross-section for open rectangular channels.