Numerical Methods & Engineering Math

FE Environmental Study Guide Exam

Track your progress as you read

Finished sections get a check mark.

1 Introduction, Mathematics & Ethics

Exam Format, NCEES CBT Logistics & Blueprint Map Calculus & Differential Equations Probability, Statistics & Regression Numerical Methods & Engineering Math Mathematics Applications for Environmental FE

2 Ethics, Economics & Materials

NCEES Model Rules & Engineering Ethics Professional Practice, Liability & Public Safety Engineering Economics: Time Value of Money Benefit-Cost, Replacement & Inflation Analysis Materials: Concrete, Corrosion & Polymers

3 Environmental Science & Chemistry

Geosynthetics & Environmental Materials Environmental Chemistry & Stoichiometry Water Chemistry: pH, Alkalinity & Hardness Microbiology, BOD Kinetics & Monod Growth Henry's Law & Phase Equilibrium for Air Quality

4 Risk Assessment & Water Resources

Environmental Science Applications EPA Four-Step Risk Assessment Paradigm Exposure, Dose, Hazard Quotient & Cancer Risk Hydrology & Rainfall-Runoff Analysis Open-Channel Flow & Manning Equation

5 Groundwater & Water Treatment

Groundwater: Darcy, Theis & Cooper-Jacob Streeter-Phelps & Water-Quality Modeling Water Treatment: Sedimentation & Filtration Disinfection, CT Rule & Chlorination Wastewater: Activated Sludge & BOD Removal

6 Advanced Treatment & Air Quality

Advanced Treatment: MBR, GAC & RO Sludge Handling & Water Reuse NAAQS, Criteria Pollutants, HAPs & AQI Gaussian Plume, Pasquill Stability & Effective Stack Height AP-42 Factors & Air Pollution Control Technologies

7 Air Permitting & Solid/Hazardous Waste

PSD, NSR, BACT, LAER, Title V & Offsets RCRA Lists, Characteristics & Generator Categories Subtitle D Landfill Design, Liners & Leachate CERCLA, UST, TCLP & Land Disposal Restrictions MSW Collection, Transfer, Composting, WTE & Recycling

8 Sustainability, EHS & Exam Day

LCA Stages, Carbon Footprint & Scope 1/2/3 Emissions LEED Credits, ISO 14001 EMS & Triple Bottom Line PELs, TLVs, Hierarchy of Controls & Confined Spaces EPCRA Reporting, TRI Form R, DOT Hazmat & RCRA Manifests Exam Format, NCEES CBT Logistics & Blueprint Map

Numerical Methods & Engineering Math

Key Takeaways Explain This Topic AI

Trapezoidal rule integrates tabulated data.

Newton-Raphson solves implicit Manning depth.

Linear algebra solves mass balances.

Bisection brackets roots.

Keep extra sig figs until the end.

Last updated: July 2026

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 ∫ a b f ( x ) d x ≈ ∑ i = 1 n − 1 f ( x i ) + f ( x i + 1 ) 2 ( x i + 1 − x i ) \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) ∫ a b f ( x ) d x ≈ i = 1 ∑ n − 1 2 f ( x i ) + f ( x i + 1 ) ( x i + 1 − x i )

Worked example: Flows (m³/h) at 0,1,2,3 h: 90, 130, 110, 85. V ≈ 90 + 130 2 ( 1 ) + 130 + 110 2 ( 1 ) + 110 + 85 2 ( 1 ) = 110 + 120 + 97.5 = 327.5 m 3 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 V ≈ 2 90 + 130 ( 1 ) + 2 130 + 110 ( 1 ) + 2 110 + 85 ( 1 ) = 110 + 120 + 97.5 = 327.5 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 ) x_{n+1} = x_n - f(x_n)/f'(x_n) x n + 1 = x n − f ( x n ) / f ′ ( x n ) .

Manning example: Q = 1.49 n A R h 2 / 3 S 1 / 2 Q = \frac{1.49}{n} A R_h^{2/3} S^{1/2} Q = n 1.49 A R h 2/3 S 1/2 with A and R_h as functions of depth y. Define f ( y ) = Q c a l c ( y ) − Q g i v e n f(y)=Q_{calc}(y)-Q_{given} f ( y ) = Q c a l c ( y ) − Q g i v e n ; 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 o u t Q_1+Q_2+Q_3=Q_{out} Q 1 + Q 2 + Q 3 = Q o u t and mass balances. Use substitution or calculator matrix solve.

Error and
Significant
Figures

Issue Mitigation
Coarse Δx Finer intervals reduce integration error
Bad initial guess Try bisection first
Unit mix Convert before numerics
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

Identify tabulated vs. implicit equation.

Select trapezoidal, Newton, or linear algebra.

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: ∫ a b f d x ≈ h 3 [ f 0 + 4 f 1 + 2 f 2 + 4 f 3 + ⋯ + f n ] \int_a^b f\,dx \approx \frac{h}{3}\left[f_0 + 4f_1 + 2f_2 + 4f_3 + \cdots + f_n\right] ∫ a b f d x ≈ 3 h [ f 0 + 4 f 1 + 2 f 2 + 4 f 3 + ⋯ + f n ]

Worked example: Flow readings every 2 h: 80, 120, 100, 140 m³/h (h = 2). V ≈ 2 3 [ 80 + 4 ( 120 ) + 2 ( 100 ) + 4 ( 140 ) ] = 2 3 [ 80 + 480 + 200 + 560 ] = 880 m 3 V \approx \frac{2}{3}[80 + 4(120) + 2(100) + 4(140)] = \frac{2}{3}[80 + 480 + 200 + 560] = 880\text{ m}^3 V ≈ 3 2 [ 80 + 4 ( 120 ) + 2 ( 100 ) + 4 ( 140 )] = 3 2 [ 80 + 480 + 200 + 560 ] = 880 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 c a l c ( y ) − Q g i v e n f(y) = Q_{calc}(y) - Q_{given} f ( y ) = Q c a l c ( y ) − Q g i v e n , 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:

Source Q (MGD) Cl (mg/L)
A Q_A 25
B Q_B 180
C 10 - Q_A - Q_B 8

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 6(25) + Q_B(180) + (4-Q_B)(8) = 500 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 + Δ t ( d C d t ) n C_{n+1} = C_n + \Delta t \left(\frac{dC}{dt}\right)_n C n + 1 = C n + Δ t ( d t d C ) 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 ≈ ∑ f i + f i + 1 2 Δ x V \approx \sum \frac{f_i+f_{i+1}}{2}\Delta x V ≈ ∑ 2 f i + f i + 1 Δ x

x n + 1 = x n − f ( x n ) / f ′ ( x n ) x_{n+1}=x_n-f(x_n)/f'(x_n) 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.