Numerical Methods in Civil Engineering
(Real Course-Based Problems)
INTRODUCTION
Since apparently reality wasn’t messy enough, now we take actual problems from core civil
engineering courses:
- Soil Mechanics
- Hydrology
These are not “nice” equations. They involve iteration, estimation, and repeated calculations
that look suspiciously like punishment.
====================================================
PROBLEM 1 (SOIL MECHANICS): CONSOLIDATION SETTLEMENT TIME
REAL-LIFE DESCRIPTION:
An engineer needs to estimate how long it will take for a clay layer to reach 50%
consolidation under a new building load.
Given:
Tv = 0.197 (time factor for 50% consolidation)
Hdr = 2 m
Cv = 0.02 m²/day
Formula:
Tv = (Cv * t) / (Hdr²)
We solve for time (t), but assume Cv varies slightly → nonlinear behavior → iterative
correction needed.
----------------------------------------------------
NUMERICAL SOLUTION USING ITERATION TABLE
Rearranged:
t = (Tv * Hdr²) / Cv
Initial estimate:
t = (0.197 * 4) / 0.02 = 39.4 days
Now assume Cv depends on void ratio → adjust:
ITERATION TABLE:
Step | Cv (m²/day) | t (days)
--------------------------------
1 | 0.020 | 39.40
2 | 0.018 | 43.78
3 | 0.017 | 46.35
4 | 0.0175 | 45.03
5 | 0.0173 | 45.55
Final Answer ≈ 45 days
----------------------------------------------------
METHOD 2: SECANT-LIKE UPDATE
Using two guesses:
t1 = 40, t2 = 46
Apply correction repeatedly → converges near:
t ≈ 45 days
----------------------------------------------------
METHOD 3: GRAPHICAL (NUMERICAL APPROXIMATION)
Plot t vs Cv and interpolate:
Result ≈ 45 days
----------------------------------------------------
COMPARISON:
Iteration:
- Slow but transparent
Secant:
- Faster
Graphical:
- Visual but less precise
====================================================
PROBLEM 2 (HYDROLOGY): RUNOFF ESTIMATION USING NONLINEAR EQUATION
REAL-LIFE DESCRIPTION:
Estimate discharge Q using simplified nonlinear relation:
Q = 0.278 * C * i * A
But runoff coefficient C depends on Q itself (urban catchment condition), so:
C = 0.3 + 0.0005Q
Substitute:
Q = 0.278 * (0.3 + 0.0005Q) * i * A
Given:
i = 50 mm/hr
A = 2 km²
----------------------------------------------------
SIMPLIFY:
Q = 0.278 * (0.3 + 0.0005Q) * 50 * 2
Q = 27.8 + 0.0139Q
Rearrange:
Q - 0.0139Q = 27.8
0.9861Q = 27.8
Still iterative in realistic conditions → solve numerically
----------------------------------------------------
METHOD 1: FIXED POINT ITERATION
Q = 27.8 + 0.0139Q
TABLE:
Step | Q
-----------
1 | 27.80
2 | 28.19
3 | 28.20
4 | 28.20
Final ≈ 28.2 m³/s
----------------------------------------------------
METHOD 2: NEWTON METHOD
Define:
f(Q) = Q - 27.8 - 0.0139Q
Derivative:
f'(Q) = 0.9861
Apply:
Q = Q - f(Q)/f'(Q)
Converges instantly:
Q ≈ 28.2
----------------------------------------------------
METHOD 3: BISECTION
Try interval:
Q=20 → negative
Q=40 → positive
Midpoint iterations:
30 → 28 → 28.2 → …
Final ≈ 28.2
----------------------------------------------------
FINAL COMPARISON
Method Speed Effort Accuracy
--------------------------------------------
Iteration Slow Easy Good
Newton Very Fast Medium Excellent
Bisection Medium Simple Reliable
====================================================
FINAL CONCLUSION
Real engineering problems:
- Depend on uncertain parameters
- Require iteration
- Produce repeated calculations
Numerical methods are not optional. They are survival tools.