Numerical Methods & Engineering Math

Quick Answer: FE numerical methods mean trapezoidal integration, root-finding (Newton-Raphson), and solving linear systems — not programming. Use them when data are tabulated or when Manning and similar equations cannot be solved algebraically.

Environmental data rarely arrive as neat continuous functions. Hourly flows, discrete sampling, and implicit hydraulic equations require numerical techniques the Handbook supports conceptually.

Trapezoidal

Rule

[ \int_a^b f(x),dx \approx \sum_{i=1}^{n-1} \frac{f(x_i)+f(x_{i+1})}{2}(x_{i+1}-x_i) ]

Worked example: Flows (m³/h) at 0,1,2,3 h: 90, 130, 110, 85.

[ V \approx \frac{90+130}{2}(1)+\frac{130+110}{2}(1)+\frac{110+85}{2}(1)=110+120+97.5=327.5\text{ m}^3 ]

Apply to pollutant mass by multiplying each interval's average flow by concentration.

Simpson's

Rule (Conceptual)

Parabolic segments through three points give higher accuracy for smooth functions. FE items more often specify trapezoidal — read the stem.

Newton-Raphson

Solve ( f(x)=0 ) iteratively: ( x_{n+1} = x_n - f(x_n)/f'(x_n) ).

Manning example: ( Q = \frac{1.49}{n} A R_h^{2/3} S^{1/2} ) with A and R_h as functions of depth y. Define ( f(y)=Q_{calc}(y)-Q_{given} ); iterate until |f| is small.

Bisection

If ( f(a) ) and ( f(b) ) have opposite signs, halve the interval — robust when Newton diverges.

Linear

Systems

Three-stream blending with unknown flows yields equations like ( Q_1+Q_2+Q_3=Q_{out} ) and mass balances. Use substitution or calculator matrix solve.

Error and

Significant

Figures

Exam

Traps

• Trapezoidal with unequal Δx still uses each actual interval width.
• Newton needs a derivative — for Manning, differentiate numerically if unsure.
• Extraneous roots — pick physical depth (positive, below banks).

FE

Strategy

1. Identify tabulated vs. implicit equation.
2. Select trapezoidal, Newton, or linear algebra.
3. Verify answer magnitude (flow, volume, depth).

Tip: Store intermediate results in calculator memory — retyping causes most numerical errors under time pressure.

FE Exam

Integration

Environmental FE items on this topic often combine regulatory classification with a quantitative step. Read the stem for the governing law (CWA, CAA, RCRA, OSHA) before selecting equations. Flag multi-step problems and return after your first pass — average time is under three minutes per question across 110 items.

Practice locating handbook relationships by keyword during timed drills. Confirm units on every constant: mg/L versus μg/m³, ft³/s versus MGD, and days versus seconds in decay and pumping problems are frequent error sources.

FE Exam Integration

Simpson's Rule — Worked Tabular Example

For even number of intervals with uniform spacing h:

[ \int_a^b f,dx \approx \frac{h}{3}\left[f_0 + 4f_1 + 2f_2 + 4f_3 + \cdots + f_n\right] ]

Worked example: Flow readings every 2 h: 80, 120, 100, 140 m³/h (h = 2).

[ V \approx \frac{2}{3}[80 + 4(120) + 2(100) + 4(140)] = \frac{2}{3}[80 + 480 + 200 + 560] = 880\text{ m}^3 ]

Trapezoidal on same data gives 840 m³ — Simpson captures curvature when the hydrograph peaks mid-interval.

Secant and Bisection for Manning Depth

When (f(y) = Q_{calc}(y) - Q_{given}), bisection between y_low = 1 ft and y_high = 6 ft guarantees convergence if signs differ.

Worked example: Q_given = 150 cfs; at y = 3 ft, Q_calc = 130 cfs; at y = 4 ft, Q_calc = 175 cfs. Bracket [3, 4]. Midpoint y = 3.5 ft → interpolate linearly: y ≈ 3 + (150-130)/(175-130) × 1 ≈ 3.44 ft.

Matrix Solution for Blending

Three-source blend: sources A, B, C with flows and chloride. Two unknown flows with total Q fixed:

Target blend 50 mg/L at 10 MGD → one chloride balance equation. If Q_A = 6 MGD, solve Q_B from (6(25) + Q_B(180) + (4-Q_B)(8) = 500) → Q_B ≈ 1.85 MGD.

Round-Off and Significant Figures

Environmental permits quote 2–3 significant figures. Carry one extra digit through multi-step CT, Manning, or dose calculations; round only at the end. Reporting 3.456789 mg/L when choices are 3.4, 3.5, 4.1 signals a unit error, not precision.

Euler's Method (First-Order ODE)

[ C_{n+1} = C_n + \Delta t \left(\frac{dC}{dt}\right)_n ]

If (dC/dt = -0.1C) with C₀ = 100 mg/L and Δt = 1 day: C₁ = 100 - 10 = 90; C₂ = 81 — matches exponential decay for small Δt. FE may ask whether explicit Euler over- or under-predicts decay (over-predicts remaining concentration for convex decay curves).

Error Propagation (Conceptual)

When result R = A × B, relative errors add in quadrature approximately. If flow is ±5% and concentration ±10%, mass load uncertainty ≈ √(5² + 10²) ≈ 11% — choose conservative design factors accordingly.

Trapezoidal Rule

[ V \approx \sum \frac{f_i+f_{i+1}}{2}\Delta x ]

[ x_{n+1}=x_n-f(x_n)/f'(x_n) ] for Manning depth when Q(y) is implicit.

Romberg and Richardson (Recognition)

FE items may state that Richardson extrapolation improves trapezoidal estimates when interval is halved. If error scales as h², combining coarse and fine estimates cancels leading error — conceptual only; follow stem formulas.

Fixed-Point Iteration

Rewrite f(x)=0 as x=g(x); iterate x_{n+1}=g(x_n). Converges if |g′(x)|<1 near root. Slower than Newton but needs no derivative — useful when f′ is messy.