Civil Engineering A path through The Default Library ======================================== How the bridges, buildings, roads, dams, pipes and water a society runs on are designed to stand up, last and keep people safe. Step 1. Calculus I: Slopes, Totals and the Bridge Between Them -------------------------------------------------------------- Why: Loads, flows and settlements change across space and time: calculus describes them. Your car's speedometer is a derivative. Its odometer is an integral. Everything a first calculus course teaches is a way of moving between the two. Look at a dashboard. The odometer says how far you have gone in total. The speedometer says how fast the total is changing right now. Press the accelerator and the speedometer moves first; the odometer follows, a little behind, like a patient accountant adding up what the speedometer reported. Nobody gave the speedometer a ruler or a clock to look at the whole trip. It reads one instant. Yet it reports a rate, which is a statement about change, from a single moment. How can a single moment say anything about change? That puzzle is what the course exists to answer, and the answer turned out to run most of modern engineering and computing. This is the second time the library has met calculus. Calculus: Derivatives, Integrals, and the Language of Motion (12/11) is the easier first read, with the tangent line and the rectangle sum. Here we go a step further: what the course trains you to see, how a computer gets it wrong, and where the idea is used afterwards. The course in one box Who takes it: first-year students in Mechanical & Aerospace Engineering and in Computer Science, and almost every other science or engineering degree. What it is for: turning "how much?" into "how fast?" and back, so that change can be calculated instead of guessed. Comes after: nothing in this series. Algebra and a little trigonometry are all it assumes. Leads to: Multivariable Calculus (34/02), Differential Equations (34/03), Probability & Statistics (34/05), Physics I: Mechanics (34/06), Thermodynamics (34/14), and from them nearly everything else on both majors. Honest line: you will not be able to differentiate a hard formula by hand after this. You will see why a flat spot on a curve is where the best answer hides, and why the area under a rate is a total. A limit is where a calculation settles The speedometer puzzle is this. Speed is distance divided by time, but in an instant no time passes and no distance is covered, so you seem to be dividing zero by zero. The way out is to stop asking what happens at zero and ask what a sequence of better and better answers is settling towards. Measure over one second, then a tenth, then a hundredth, and watch the average speed home in on a number. That number is the limit. speed now = limit of (distance change ÷ time change) as the time change shrinks to 0 A graduate does not carry the formal definition around, but carries the habit: if a quantity looks like 0 ÷ 0, ask what it is approaching. Take sin(x) ÷ x as x goes to zero, with x in radians. At x = 0.1 it is 0.99833; at 0.01, 0.99998; at 0.001, 0.9999998. It never "is" 1 at zero, since it is undefined there, but it settles on 1. That one limit is why the derivative of sine turns out to be cosine. Interactive · Shrink the gap until the computer gives up We estimate the slope of ex at x = 1 by (e1+h − e) ÷ h. The true slope is e = 2.718281828459045. Make the gap h smaller. Gap h = 10−k The error falls as the gap shrinks, as the limit promises, until about h = 10−8, where it bottoms out near 7 × 10−9. Shrink further and it gets worse: the two values being subtracted agree in almost all of their roughly 16 stored digits, so the subtraction wipes out the information. At 10−16 the answer is exactly 0. This is a computer-science lesson hiding in a calculus course: the limit is a mathematical idea, and a machine can only approach it so far. The derivative is a slope, and a slope is a rate Zoom in on any smooth curve and it straightens out. The steepness of that straight piece is the derivative. In plain words: how much the output moves per unit nudge of the input. If the curve is distance against time, the derivative is speed. If it is temperature against position, it is how sharply it warms as you walk. Rate and slope are one thing seen two ways. f′(x) = limit of [f(x+h) − f(x)] ÷ h as h → 0 The graduate's real asset is not computing derivatives but reading them. Where the curve is rising, the derivative is positive. Where it is flat, the derivative is zero. Where the curve is steepest, the derivative is at its peak. Try it: the lower graph is the derivative of the upper one, drawn point by point. Interactive · A curve and its slope, side by side x³ − 3x sin x bell: e^(−x²) Position x Check the three curves. Every time the top curve peaks or dips, the bottom one crosses zero. Every time the top is climbing, the bottom is above the axis. For x³ − 3x the slope is 3x² − 3, which is zero at x = −1 and x = 1: the hill top and the valley floor. The chain rule: rates multiply down a chain Rates often come in chains. A bicycle's pedals turn the chainring, the chain turns the rear wheel, the wheel moves the bike. If the wheel turns 3 times per pedal turn, and the bike moves 2 metres per wheel turn, the bike moves 6 metres per pedal turn. You multiply the ratios. That is the chain rule, and it is the most used rule in the course. d(outer)/d(input) = d(outer)/d(middle) × d(middle)/d(input) The classic worked case is a balloon. Its volume depends on its radius, and the radius depends on time. At radius r the surface has area 4πr², so one extra centimetre of radius adds a shell of volume about 4πr² cm³: that is dV/dr. Multiply by how fast the radius is growing and you have how fast the volume is growing. At r = 10 cm, growing 2 cm each second, volume grows at 4π × 100 × 2 ≈ 2,513 cm³ per second. Pull the sliders and confirm it. Interactive · Inflate a balloon Radius now Radius is growing at The picture shows why it works: the shaded shell is what one second adds, a skin of area 4πr² and thickness equal to the radius growth. The same rule runs inside every neural network, where the "chain" is dozens of layers long and the multiplication is called backpropagation (Machine Learning, 34/32; Multivariable Calculus, 34/02). Where the slope is zero, something is best At the top of a hill you stop going up and have not yet started going down. The slope is zero. So to find the biggest or smallest value of anything smooth (strongest, cheapest, lightest, fastest) look for where its derivative is zero. That is optimisation, and it is why engineers love the derivative. Take a 30 cm by 20 cm sheet of card. Cut a square of side x from each corner and fold up the sides to make an open box. Its volume is V = x(30 − 2x)(20 − 2x). Cut too little and the box is shallow. Cut too much and it is narrow. In between there is a best x, and the slope tells you where: at x ≈ 3.92 cm the volume peaks at about 1,056 cm³. Interactive · The best box Corner cut x Jump to the flat spot Notice what the slope does as you drag. To the left of the peak the slope is positive, so a bigger cut helps. To the right it is negative, so a smaller cut helps. At the peak it is zero and nothing nearby is better. Gradient descent, the algorithm that trains machine-learning models, is exactly this reasoning repeated: look at the slope, step downhill, repeat. The integral adds up a rate into a total The odometer is the speedometer added up over time. Draw speed against time; the area between the curve and the time axis is distance. That is the integral. Anywhere a rate is known and a total is wanted (litres per minute into a tank, watts into a battery, newtons over a distance) the answer is an area. Try it before you read on. A car accelerates steadily from 0 to 20 m/s in 10 seconds, cruises at 20 m/s for 30 seconds, then brakes steadily to a stop in 10 seconds. How far has it gone? Interactive · Predict, then reveal 600 m 800 m 1,000 m 1,200 m Choose an answer. The shape is two triangles and a rectangle: ½ × 10 × 20 + 30 × 20 + ½ × 10 × 20 = 100 + 600 + 100 = 800 m. The common wrong guess, 1,000 m, uses the top speed for all 50 seconds and forgets that the car spends time slower than that. The Fundamental Theorem: two jobs, one idea Here is the surprise that makes the course a single subject. Build a running total from a rate: the amount in the tank after time T. Then ask how fast that total is growing at T. The answer is just the rate at T. Accumulating and differentiating undo each other, the way the odometer and speedometer do. Newton and Leibniz, working separately in the late 1600s, each saw this. if A(T) = area under the rate up to T, then A′(T) = the rate at T Its practical power is that you can find a total without adding anything: find a function whose slope is your rate and subtract its values at the two ends. Below, water flows into a tank at a rate that wobbles: 4 + 3 sin(t/2) litres per minute. Drag the time and watch the right-hand graph, the tank's total, steepen exactly when the rate on the left is high. Interactive · The tank and its rate Time T The readout checks the theorem two ways: area by adding 600 thin slices, and the exact formula 4T + 6(1 − cos(T/2)). At T = 12 minutes both give 48.24 litres. The slope of the total at that moment is 4 + 3 sin 6 ≈ 3.16 litres per minute, the rate at that moment. Every smooth thing is nearly a straight line The last habit is the most useful in daily engineering. Close to any point, a smooth curve is almost its tangent line, so you can replace it with a straight line and be nearly right. The small-angle rule sin x ≈ x is why a pendulum swings at a steady rhythm. The rule √(1 + x) ≈ 1 + x/2 gives mental square roots: √1.1 ≈ 1.05, true value 1.0488. f(a + x) ≈ f(a) + f′(a) · x   (good when x is small) How small is small? The error shrinks faster than the step. Halve the step and the error falls about four times for √(1 + x) and ex, and about eight times for sine. Interactive · How good is the straight line? sin x ≈ x √(1+x) ≈ 1 + x/2 eˣ ≈ 1 + x Step x from the point What a graduate just sees A rate and a slope are the same thing, and "how fast?" always means a derivative. A total is an area under a rate, and adding up thin slices is a calculation, not a mystery. Rates multiply along a chain: gear ratios, units, layers of a network. The best value of a smooth quantity sits where its slope is zero. Differentiating and accumulating are inverse jobs, so a total can often be found without adding. Up close, everything smooth is a straight line, and engineers exploit that constantly. Subtracting two nearly equal numbers on a computer destroys precision, so a naive slope estimate gets worse when the gap gets too small. Where it shows up next Physics I (34/06) and Dynamics (34/12): velocity is the derivative of position, acceleration of velocity, and Newton's law ties force to the second one. Differential Equations (34/03): when a law is written as a rate, solving it means integrating, which is the whole next course. Thermodynamics (34/14) and Heat Transfer (34/16): heat added is an integral of power; temperature profiles are read from slopes. Probability & Statistics (34/05): the chance of a range of outcomes is the area under a density curve. Algorithms (34/24) and Machine Learning (34/32): growth rates of running time are compared with limits, and models are trained by following slopes downhill. If you want to go deeper 📘 Textbook: James Stewart, Daniel Clegg and Saleem Watson, Calculus: Early Transcendentals (9th ed., Cengage). The most-assigned text; best for worked examples and a huge problem bank. 📘 Textbook: Hass, Heil, Bogacki and Weir, Thomas' Calculus (15th ed., Pearson). Strong on applications, popular in engineering departments. ▶️ Gentler: Grant Sanderson's Essence of Calculus series on 3Blue1Brown, a visual tour of the same ideas. Also Steven Strogatz, Infinite Powers, for the story. 🌐 Free course: MIT OpenCourseWare 18.01SC Single Variable Calculus, with video lectures, notes and problems. Khan Academy also covers it. The takeaway Calculus I teaches two moves and one bridge. The derivative turns a total into a rate (a slope); the integral turns a rate into a total (an area); and the Fundamental Theorem says they undo each other. From that come the habits of mind: chain rates together, look for zero slope to find the best, and treat smooth things as lines up close. Our slope estimate for ex was best near a gap of 10−8 and useless at 10−16, a reminder that the ideas are exact and the machines are not. Try this: next time you are in a car or on a bike, watch the speedometer and ask what the area under its history means. Then on a hill, find the spot where the road is flat. That is a derivative equal to zero. Sources & Further Reading Sources. Every number in the interactives was computed in node: the forward-difference error of ex at x = 1 for h = 10−1 to 10−16 (best 6.6 × 10−9 at 10−8, exactly 0 at 10−16, in IEEE double precision); sin(x)/x at 0.1, 0.01, 0.001; the box V(x) = 4x³ − 100x² + 600x, zero slope at x = (50 − √700)/6 = 3.9237, V = 1,056.3 cm³; the balloon 4π × 10² × 2 = 2,513.27 cm³/s; the car, 800 m; the tank, A(12) = 48.239 L; small-step errors of the three approximations. Book editions checked against the publishers' pages (Cengage, Pearson) and the course numbers against MIT OpenCourseWare. 🌐 MIT OpenCourseWare 18.01SC Single Variable Calculus ▶️ 3Blue1Brown: Essence of Calculus 🌐 Khan Academy: Calculus 1 📘 Stewart, Clegg and Watson, Calculus: Early Transcendentals, 9th ed. (Cengage) 📘 Hass, Heil, Bogacki and Weir, Thomas' Calculus, 15th ed. (Pearson) 📘 Steven Strogatz, Infinite Powers, a popular history of calculus Step 2. Physics I: Mechanics — Three Ways to Keep the Books ----------------------------------------------------------- Why: Every structure is Newton's laws standing still. Why doubling your speed quadruples your stopping distance, why a skater spins up when she pulls her arms in, and the three questions an engineer asks before touching a calculator. A car doing 50 km/h stops in about 14 metres on a dry road. Most people guess that at 100 km/h it takes about 28. It takes about 56. The extra speed doesn't just cost twice the distance, it costs four times, and the same fact is why a crash at 100 km/h is not "twice as bad" as one at 50. Nobody needs a formula to be afraid of that, but a person who has done a first mechanics course sees it: the brakes have to get rid of the car's energy, energy grows with the square of speed, and the road can only take it away at a steady rate per metre. That one chain of thought is most of the course. Physics I is less a list of topics than a way of keeping accounts. There are three ledgers, force, energy and momentum, and a good physicist picks the one in which the problem is easiest. If you want an easier first read, the library has The Four Forces (13/01), Energy: The One Currency (13/02) and Orbits (13/11). This piece goes further: it makes you use the ledgers. The course in one box Who takes it: first-year students of Mechanical & Aerospace Engineering, and in many programmes Computer Science too, as the science requirement. What it's for: turning "what will this object do?" into something you can calculate. Comes after: Calculus I (34/01). Leads to: Physics II: Electricity & Magnetism (34/07), Statics (34/10), Dynamics (34/12), Materials Science & Manufacturing (34/13), Thermodynamics (34/14). Honest line: you won't be able to design a suspension after this. You will see why a car's tyres, not its engine, decide how fast it can take a corner. Newton's laws are one accounting system Isaac Newton set out his three laws in the Principia in 1687. Read as a bookkeeper would, they say: momentum (mass times velocity) is the thing being tracked; a force is a deposit or withdrawal; and every deposit has a matching withdrawal somewhere else. The first law says that with no net force the account doesn't change: things keep their velocity. The second says that a net force changes momentum at a steady rate, which for a fixed mass is the famous rule below. The third says forces come in pairs between two bodies: if you push a wall, the wall pushes you, and the pair never acts on the same body, which is why they don't cancel. Fnet = m · athe net force on a body is its mass times its acceleration Worked number: a 1,500 kg car getting from 0 to 100 km/h (27.8 m/s) in 8 seconds accelerates at 3.47 m/s², so the road must push it forward with about 5.2 kN. Two ice skaters of 50 kg and 70 kg who push apart from rest leave with equal and opposite momentum: if the heavier moves at 1.0 m/s the lighter goes at 1.4 m/s the other way. Total momentum stays at zero. Try the classic: you on a scale in an elevator. The scale doesn't read your weight, it reads the force pushing up on you. Interactive · The elevator scale (a 70 kg rider) Speeding up, going up Cruising Speeding up, going down Cable cut Elevator acceleration (up is positive) Accelerating upward at 2 m/s² the scale reads about 827 N, like 84 kg. At constant speed it reads your true weight, whether the elevator is moving or still. In free fall it reads zero. That is exactly the "weightless" astronaut of Orbits. Free-body diagrams, and friction as a force that only pushes as hard as it must The graduate's first reflex on any problem is to draw one object, alone, and every push on it: weight, the surface pushing back, friction, tension. Nothing else. That is a free-body diagram. Add the forces as vectors, split them along the directions that matter, and F = ma does the rest. Friction is the force beginners get wrong. Static friction is not a fixed number: it pushes exactly as hard as needed to stop sliding, up to a ceiling of μs times the normal force (how hard the surface pushes back). Only once the ceiling is passed does the block slide, and then friction drops to μk × normal force. f ≤ μs · Nfriction is at most a coefficient times how hard the surfaces press together On a ramp the weight splits into a part along the slope, mg sin θ, and a part pressing into it, mg cos θ. Sliding begins when the first beats μs times the second, which gives tan θ = μs. The mass cancels. Check that in the tool: change the mass and nothing about when it slips changes. (The coefficients are round textbook values; real surfaces vary a lot.) Interactive · A block on a ramp Wood on wood Rubber on concrete Steel on ice Ramp angle Mass Arrows are drawn as fractions of the weight. With μs = 0.5 the block holds until 26.6° and then slides. At 40° it accelerates down at about 4.05 m/s² (with μk = 0.3). This is why a pile of sand has a steepest slope, and why designers of ramps, screws and wedges care so much about one number. Energy: the balance that doesn't care about the path Force problems get painful when the path is complicated. Energy is the escape. The work a force does is force times distance moved along it, and the work–energy theorem says the net work on a body equals its change in kinetic energy. When only gravity or springs are involved, the total of kinetic plus potential energy is simply conserved, whatever the route. KE = ½ · m · v²speed counts squared; this is why speed is so costly A bobsled dropping 30 m ends at √(2 × 9.81 × 30) = 24.3 m/s (87 km/h) with no friction, however many bends the track has. Now back to the opening: the brakes do work μ × m × g × d against the car's ½mv², so the stopping distance is v² ÷ (2μg). Double v, quadruple d. Interactive · Brake, or hit Dry (μ 0.7) Wet (0.5) Snow (0.2) Ice (0.1) Speed when you brake Obstacle ahead Braking only, no reaction time. μ values are round, typical figures for tyres. The sharpest result: brake from 100 km/h with an obstacle 14 m away, the distance that stops a 50 km/h car on a dry road. You still hit it at about 87 km/h, with 75% of the energy you started with: 14 metres of hard braking has removed only a quarter of it. Momentum: what survives every crash Energy has a weakness: in a collision much of it turns to heat, sound and bent metal, so the energy ledger is no use. Momentum, though, is conserved in any collision, because the two forces are equal and opposite for the same short time (the third law again). The total before equals the total after, with direction as a sign. How bouncy the collision is gets one number, the restitution e: 1 for perfectly bouncy, 0 for stuck together. Momentum doesn't care about e. Kinetic energy does. Interactive · Two carts collide Equal masses, bouncy Stick together Light hits heavy Blue mass Blue velocity Gold mass Gold velocity Bounciness e Velocity to the right is positive; arrows show velocity. With equal masses and e = 1, the moving cart stops dead and the other leaves with its speed: the toy known as Newton's cradle. Stuck together, two 1,000 kg cars (one at 10 m/s, one parked) leave at 5 m/s and half the kinetic energy, 25 kJ, becomes damage. That lost half is the reason cars have crumple zones: they spend the energy deforming metal instead of people. Rotation: the same laws with new names Everything above has a spinning twin. Force becomes torque, force times the perpendicular distance from the pivot: 80 N on a 0.25 m wrench is 20 N·m, and a longer wrench gives more for the same push. Mass becomes moment of inertia I, which counts not just how much mass there is but how far from the axis it sits. Momentum becomes angular momentum L = I × ω. τ = I · α   ·   L = I · ωwith no outside torque, L is conserved Mass far from the axis counts for much more, because I grows with the square of distance. That's why a solid disc beats a hoop of the same mass and radius down a ramp: the hoop keeps all its mass at the rim, so it spends more of its energy on spinning (the disc's acceleration is two thirds of g sin θ, the hoop's one half, whatever the size or mass). And it's why a skater spins faster when she pulls her arms in. Interactive · The skater's weights Arms out Halfway Tucked in Weights' distance from the axis A model: body I = 0.9 kg·m², plus two 2 kg weights; start at one turn per second with arms out. From arms out to tucked in, I falls from 3.15 to 0.99 kg·m², so the spin rises from 1 to about 3.2 turns per second. The energy also rises, from 62 J to 198 J, and the difference is the work her muscles do pulling the weights inward. Circular motion and gravity: a force that points at the middle Going round a bend at constant speed is still accelerating, because the velocity keeps changing direction. The acceleration points to the centre and equals v²/r, so something must push sideways with force mv²/r. On a flat road that something is friction. No more friction, no more turn: the car continues in a straight line, which feels, from inside, like being thrown outward. vmax = √(μ · g · r)fastest speed a flat bend of radius r allows Notice the engine's power isn't in the formula. On a 50 m bend with μ = 0.5 (wet) the limit is 56 km/h; dry (0.7), 67 km/h; ice (0.1), 25 km/h. Interactive · Take the bend Dry (0.7) Wet (0.5) Ice (0.1) Speed Radius of the bend Replace friction with gravity and you get orbits. Set v²/r equal to GM/r² and the period T satisfies T² ∝ r³, Kepler's third law. Check it: the Space Station's lap of 92.6 minutes at 6,779 km from Earth's centre predicts, for the Moon at 384,400 km, a period of 27.5 days. The real month is 27.3. One law, from a bend in a road to the Moon. Which ledger? Choosing force, energy or momentum This is the real skill. A hard problem in one ledger is a one-liner in another. Three rough rules: if the question asks for a force or a time, use F = ma; if it asks for a speed between two positions, and the path is messy, use energy; if there's a collision or explosion, use momentum. If something spins, use angular momentum. Try five. Interactive · Pick the first tool Force Energy Momentum Angular momentum Next question What a graduate just sees Stopping distance goes with speed squared: double the speed, four times the distance. A scale reads the push on you, not your mass: it changes in an accelerating elevator. Static friction is only as large as it has to be, up to a limit that depends on μ and the pressing force. Energy ignores the path; use it when the route is complicated. Momentum survives every collision; kinetic energy rarely does. Mass far from the axis counts for much more in rotation, because distance is squared. A turn needs a sideways force; remove the grip and you go straight. Weightless means falling together, not "no gravity". Where it shows up next Statics (34/10) is F = ma with a = 0: the free-body diagram becomes the bridge or crane that stays put. Dynamics (34/12) and Vibrations (34/17) take rotation and the second law into machines and springs, usually as differential equations (Differential Equations, 34/03). Thermodynamics (34/14) is the energy ledger with heat added: the lost half of that crash is where it starts. Physics II (34/07) redoes the same accounting for charges, with Coulomb's force and voltage as energy per charge. In Programming Fundamentals (34/09), you'll turn F = ma into a simulation loop; in Computer Science game engines and robotics code, it's the physics step. If you want to go deeper 📘 Standard textbooks. Halliday, Resnick & Walker, Fundamentals of Physics (12th ed., Wiley, 2021), or Young & Freedman, University Physics (16th ed., Pearson); the usual first-year choices, with worked examples and thousands of problems. For a sharper, more mathematical route, Kleppner & Kolenkow, An Introduction to Mechanics (2nd ed., Cambridge). 📘 Gentler, free. The Feynman Lectures on Physics, Volume I, free online at feynmanlectures.caltech.edu; chapters 4 (energy) and 10 (momentum) are a good start. 🌐 Free course. MIT OpenCourseWare 8.01SC Classical Mechanics (Fall 2016): lectures, problem sets and solutions. The takeaway Mechanics is three ledgers kept over the same world. Force tells you what changes motion and when; energy lets you skip the path; momentum survives collisions; and each has a spinning twin. Doubling speed quadruples stopping distance because energy goes with v². The graduate's skill is choosing which ledger makes this problem easy. Try this: next time you're a passenger, estimate how far the car needs to stop at the speed you're doing, using "14 m at 50 km/h, four times that at 100". Then notice how far ahead the next car really is. Sources & Further Reading Sources. Newton's laws first published in Philosophiae Naturalis Principia Mathematica (1687). Textbook physics throughout. Friction coefficients are rounded from published tables, which vary widely between sources (for example Engineering Toolbox's "Friction coefficients for common materials", which gives rubber on dry concrete 0.6–0.85, on wet concrete 0.45–0.75). Every number in the interactives was computed in node and checked by hand: 14.05 m and 56.2 m braking distances at 50 and 100 km/h (μ = 0.7); 86.6 km/h residual speed; 5.21 kN; 24.3 m/s; skater I = 3.15 and 0.99 kg·m², 3.18 turns per second, 62 J and 198 J; 56.4, 66.7 and 25.2 km/h bend limits; Kepler scaling from the Space Station (92.6 min at 6,779 km) to the Moon (27.46 days from 384,400 km). 🌐 MIT OpenCourseWare: 8.01SC Classical Mechanics 📘 The Feynman Lectures on Physics, Vol. I (free online, feynmanlectures.caltech.edu) 📘 Kleppner & Kolenkow, An Introduction to Mechanics, 2nd ed. (Cambridge University Press) 🌐 Engineering Toolbox: friction coefficients Step 3. Linear Algebra: Machines That Move Space ------------------------------------------------ Why: A building frame has thousands of joints; matrices solve them together. A grid of numbers is really a machine that moves all of space at once. Learn to see it, and circuits, bridges, search engines and AI all start to look alike. You pinch two fingers on a photo and it rotates and zooms smoothly under them. A game character turns to face you and every corner of its face swings round together. A search for "best pizza near me" ranks millions of pages in a fraction of a second. A bridge designer presses "solve" and a screen of colour shows how every beam is being squeezed. Underneath all four is the same object: a grid of numbers, and one habit of thought about what a grid of numbers does. A graduate of a linear algebra course no longer sees a table of numbers. They see a machine that grabs the whole plane, or the whole of space, and moves it. This piece builds that picture with seven things you can drag. You will not come out able to compute a determinant of a 5×5 by hand (nobody does that on purpose). You will come out seeing why a system of equations sometimes has no answer, why a bridge hums at particular pitches, and why the chips that run AI are built for one operation only. The course in one box Who takes it: nearly every engineering and computer science student, usually in the first or second year. It is the single most reused maths course in both majors. What it's for: handling many unknowns and many dimensions at once, with one set of tools. Comes after: only school algebra and a little geometry. (Calculus helps but is not required.) Leads to: Statics (34/10), Control Systems (34/20) and Machine Learning (34/32) list it as a prerequisite; Vibrations (34/17), Computer Architecture (34/25) and Multivariable Calculus (34/02) lean on it constantly. Honest line: you won't be able to run a finite element program after this piece. You will see why it spends nearly all its time solving one equation, Ax = b. If you want an easier first read, try the library's Vectors, Matrices and Linear Transformation (12/14); this one goes further into solving, eigenvectors, fitting and why hardware cares. A matrix moves all of space at once Start with the picture a graduate keeps. A 2×2 matrix is not four numbers. It is two arrows. The first column says where the arrow pointing one step right (call it i) ends up. The second column says where the arrow pointing one step up (j) ends up. Once you know where those two land, you know where everything lands, because every point is some amount of i plus some amount of j, and the machine is straight and even: grid lines stay straight, and evenly spaced lines stay evenly spaced. Ax = x1(column 1) + x2(column 2) Read that in words: the matrix times a point is "that many of the first column, plus that many of the second". Try it. The grey grid is the original plane; the blue grid is where the matrix sends it. Interactive · Move the plane Do nothing Rotate 30° Stretch Shear Mirror Column 1, right part (a) Column 2, right part (b) Column 1, up part (c) Column 2, up part (d) Without JavaScript: the rotate-30° matrix sends the arrow (1, 0) to (0.87, 0.5) and (0, 1) to (−0.5, 0.87). Every grid line turns 30° and the unit square keeps area 1. Notice what the shaded unit square does. In "Rotate" it keeps its size. In "Stretch" it doubles. In "Shear" it leans over but keeps its area. That area change is the next idea. The determinant is how much area gets scaled The determinant is one number that says how much the machine scales area. If the unit square becomes a shape with area 5, the determinant is 5. If it becomes a mirror image, the determinant is negative: area is kept (or scaled) but orientation is flipped. For a 2×2 matrix with columns (a, c) and (b, d) the rule is short: det = ad − bc Worked number: the matrix with columns (2, 1) and (1, 3) has det = 2×3 − 1×1 = 5, so it stretches every area in the plane by a factor of five. In three dimensions the determinant scales volume. The interesting case is zero. A determinant of zero means the machine flattens space: the plane is squashed onto a line (or a point), and the area becomes nothing. Squashing cannot be undone, because many different starting points land on the same spot. Drag the slider until the determinant reaches zero. Interactive · Squash the plane flat Tilt the second arrow: s Without JavaScript: the matrix with columns (1, 0.5) and (s, 1) has determinant 1 − 0.5s. At s = 2 the determinant is 0, the whole plane lands on one line, and the points (0, 0) and (2, −1) both land on the origin. The orange dot is a probe at (2, −1). At s = 2 it lands exactly on the origin, the same place as the point (0, 0). Two different inputs, one output: no machine could tell you which one you started from. That is the whole meaning of "this matrix has no inverse", and it is why a zero determinant is the alarm bell in every solver you will ever meet. Solving many equations at once: Ax = b Now use the machine backwards. Given the output b, which input x did the machine start from? That is the equation Ax = b, and it is the most common problem in applied mathematics. Every time you have several unknowns that are tied together by several straight-line relationships, you have one. A circuit is the cleanest example. Two loops share a resistor. Kirchhoff's voltage law says the voltage gained around each loop equals the voltage lost, which gives one equation per loop. With a 10 V battery, R1 = 2 Ω, R2 = 3 Ω (shared) and R3 = 4 Ω: 5I1 − 3I2 = 10−3I1 + 7I2 = 0 The determinant is 5×7 − 3×3 = 26, so there is exactly one answer: I1 = 2.69 A and I2 = 1.15 A. Change the resistors below and watch the matrix and the answer change together. Interactive · Solve a circuit with a matrix Battery V R1 (loop 1 only) R2 (shared) R3 (loop 2 only) Without JavaScript: with 10 V, 2 Ω, 3 Ω and 4 Ω the equations are 5I1 − 3I2 = 10 and −3I1 + 7I2 = 0, its determinant is 26, and the loop currents are 2.69 A and 1.15 A. The reason an engineer cares: nothing about this method knows it is a circuit. A truss is the same. At each joint of a bridge, the horizontal forces must add to zero and so must the vertical ones: two equations per joint. A simple triangle truss has three joints and so six equations, for three member forces and three support reactions, and the answer is one solve of Ax = b. A real structure has thousands of joints. Finite element programs routinely solve systems of many thousands of equations, and research codes have reported solving over a million unknowns at once. The method is the same one you just used; only the size has changed. (Statics, 34/10, is this idea applied to forces.) A zero determinant in a structure's equations is a warning: the structure may be able to move without any member pushing back. Eigenvectors: the directions a transformation only stretches Most arrows are turned by a matrix as well as stretched. But for many matrices there are special directions where the arrow comes out pointing the same way (or exactly backwards), only longer or shorter. Those are the eigenvectors, and the stretch factor along each is its eigenvalue. (Eigen is German for "own": the matrix's own directions.) Av = λv In words: feed the matrix the special arrow v and it hands back the same arrow multiplied by a plain number λ. For the matrix with columns (2, 1) and (1, 2) the special directions are 45° (stretched by 3) and 135° (stretched by 1). Every other arrow gets bent towards the stretchy one. Turn the dial and look for the moments the two arrows line up. Interactive · Find the directions that only stretch Stretch & mix Shear Rotation Two springs Direction of the input arrow Without JavaScript: for the matrix with columns (2, 1) and (1, 2), an arrow at 45° comes out at 45° three times longer, and an arrow at 135° comes out at 135° unchanged. A pure rotation has no such direction at all. Three things to notice. The shear has only one special direction (along the x-axis, where nothing happens at all). The rotation has none, because every arrow turns. And "Two springs" has two special directions with different stretch factors, 1 and 3; hold that thought. Eigenvectors at work: bridges that hum, pages that rank Vibration. Picture two equal masses on springs between two walls. Push them any way you like and the motion looks messy. But there are exactly two clean patterns: both swing together, or they swing in opposite directions. These are the modes, and they are the eigenvectors of the "two springs" matrix above, (1, 1) and (1, −1). The eigenvalues 1 and 3 fix how fast each mode vibrates: the squared frequency is the eigenvalue times k/m, so the opposite-motion mode vibrates √3 ≈ 1.73 times faster. Any messy motion is just a mix of the modes. That is how engineers find which pitches will shake a wing, a bridge deck or a drive shaft (Vibrations, 34/17). Ranking pages. Google's founders modelled a web surfer who follows a random link from each page, and now and then jumps to a random page instead. The long-run share of time spent on each page is its PageRank, and it is the eigenvector (with eigenvalue 1) of the matrix describing one step of that surfer. In their 1998 paper Brin and Page used a "damping" probability of 0.85 for following links. You can find the eigenvector without solving anything fancy: start with equal scores and apply the machine over and over. Press Step. Interactive · A four-page web finds its eigenvector Step onceRun 20 stepsReset Without JavaScript: A links to B and C, B links to C, C links to A, D links to C. After many steps with damping 0.85 the ranks settle to A 0.373, B 0.196, C 0.394, D 0.038. D, which nobody links to, keeps only the base share. After about twenty steps the numbers stop changing: the vector has become an eigenvector, because applying the machine again returns the same vector. Page C wins because two pages point at it and one of them (A) is itself well-ranked. Least squares: the best answer when there is no exact one Real data is never perfect. Six measurements that "should" lie on a straight line will not, so the equations Ax = b have more equations than unknowns and no exact solution. Do you give up? No. You choose the line whose total mistake is smallest. Measure each miss vertically, square it, and add up the squares. The least squares line makes that total as small as possible. In pictures: each squared miss is literally a square. Drag the line and watch the squares shrink and grow. Interactive · Make the squares as small as possible Slope Height at x = 0 Snap to the least-squares line Without JavaScript: for the six points (1, 1.2), (2, 2.9), (3, 3.1), (4, 4.8), (5, 5.2), (6, 7.1) the best line is y = 1.089x + 0.24, with total squared miss 0.798. The guess y = x gives 2.75. The graduate's picture is geometric: the right-hand side b is a point outside the flat sheet of everything the matrix can reach, and the best answer is the point on the sheet closest to it, found by dropping a perpendicular. That perpendicular condition gives the one equation to remember: ATA x̂ = ATb In words: multiply both sides by the matrix flipped on its side and you get a square system you can solve. The history is a priority fight: Legendre published the method in 1805, Gauss said he had used it since 1795 and published in 1809, and Gauss's use of it to predict where the asteroid Ceres would reappear in 1801 made it famous. Today it sits under every trend line, every calibration of a sensor, and the first model most machine-learning students ever train. Why GPUs and AI are matrix multiplication Multiplying two matrices is just many dot products: every cell of the answer is one row of the first matrix paired with one column of the second, multiplied and added. Click a cell below and see the row and column that make it. Interactive · Build a product cell by cell Tap a cell of C = A × B. Matrix size n (n × n times n × n) 1 billion multiplies per second 100 billion per second 10 trillion per second The speeds are round illustrative numbers, not any particular chip. Without JavaScript: multiplying two n×n matrices takes n×n×n multiplications. For n = 1,000 that is 1,000,000,000, and the 1,000,000 cells of the answer can all be computed independently. Two facts follow. First, the work is huge: n×n matrices need n³ multiplications, so 1,000×1,000 is a billion. Second, the cells are independent, so a thousand workers (or a million) can each take cells at the same time. That is exactly the shape of work a graphics chip was built for, because drawing a 3D scene means pushing every corner of every triangle through a matrix that rotates, moves and projects it. Modern neural networks turned out to be the same shape. A layer of a network is a matrix of learned numbers times a vector (or a matrix) of inputs. In their scaling-law paper, Kaplan and colleagues estimate that a transformer language model with N parameters needs roughly 2N arithmetic operations to process one token. GPT-3, with 175 billion parameters, therefore needs about 3.5×1011 operations for each token, nearly all of them inside matrix products. NVIDIA's Volta GPUs (2017) put dedicated "tensor core" circuits on the chip that do nothing but a small matrix multiply-and-add (4×4) per clock tick. When people say AI runs on GPUs, they mean it runs on hardware shaped like linear algebra. (Computer Architecture, 34/25, and Machine Learning, 34/32.) What a graduate just sees A matrix is a machine that moves all of space; its columns say where the basis arrows go. The determinant is the area (or volume) scale factor, and zero means the machine flattens space and cannot be undone. "Solve Ax = b" is the same question for a circuit, a truss, a fitted line and a game physics step. If there are as many independent equations as unknowns there is one answer; fewer, many answers; contradictory, none. Eigenvectors are the directions a machine only stretches; they simplify any problem the machine drives. A vibrating structure is a mix of modes, one per eigenvector; the eigenvalues are the pitches. When there is no exact answer, take the closest one: project onto what you can reach. Matrix multiplication is many independent dot products, which is why it parallelises and why hardware is built around it. Where it shows up next Statics (34/10) and Mechanics of Materials (34/11): joint and member forces are a linear system; stiffness matrices tie loads to deflections. Vibrations (34/17) and Dynamics (34/12): mode shapes and natural frequencies are eigenvectors and eigenvalues of mass and stiffness matrices. Control Systems (34/20): the eigenvalues of a system's matrix say whether a drone's wobble dies away or grows. Machine Learning (34/32): data is a matrix, a layer is a matrix, and training adjusts matrices; least squares is the model you meet first. Computer Architecture (34/25) and graphics: GPUs and tensor cores exist to multiply matrices fast. If you want to go deeper 📘 Gilbert Strang, Introduction to Linear Algebra (Wellesley-Cambridge Press; 6th edition, 2023). The standard first course for engineers; best for the four fundamental subspaces and for application. 📘 David C. Lay, Steven R. Lay and Judi J. McDonald, Linear Algebra and Its Applications (Pearson; 6th edition, 2021). Gentle, example-heavy, widely assigned. For proofs, Sheldon Axler's Linear Algebra Done Right (Springer; 4th edition, free to read online at linear.axler.net) is the elegant theoretical route. ▶️ 3Blue1Brown, Essence of Linear Algebra: the best visual series for learning alone; it builds exactly the "matrix as moving space" picture used here. 🌐 MIT OpenCourseWare 18.06, Linear Algebra (Gilbert Strang): full video lectures, free. The takeaway Linear algebra is the study of machines that move all of space at once. Read a matrix by its columns, read its determinant as the area scale, solve Ax = b when you need to run the machine backwards, look for the eigenvector directions when you want it to be simple, and project onto what you can reach when there is no exact answer. Try this: next time an app zooms, rotates or recolours a photo, picture the two column arrows moving. And for the fun of it, in the circuit widget raise R2 to 10: the shared resistor starts to decide what both loops do. Sources & Further Reading Sources. Strang, Introduction to Linear Algebra (6th ed., 2023) and MIT 18.06 for the standard course content. Kirchhoff loop equations and the spring-mass modes are textbook derivations; every number above was recomputed in code. Legendre's 1805 publication and Gauss's priority claim and Ceres prediction: Wikipedia, "Least squares"; MacTutor biography of Legendre. PageRank and the 0.85 damping factor: Brin and Page, "The Anatomy of a Large-Scale Hypertextual Web Search Engine" (1998). Volta tensor cores (4×4 matrix multiply-accumulate): Markidis et al., "NVIDIA Tensor Core Programmability, Performance & Precision" (2018). About 2N operations per token: Kaplan et al., "Scaling Laws for Neural Language Models" (2020); GPT-3 at 175 billion parameters: Brown et al. (2020). Million-unknown structural solves: NASA Technical Reports Server and OSTI conference papers on parallel sparse solvers. Go deeper. The books and courses above. Step 4. Programming Fundamentals: Instructions, State and Taking Things Apart ----------------------------------------------------------------------------- Why: Civil engineers model in code and spreadsheets long before they pour concrete. A computer can't guess what you meant. Everything in a first programming course is a way of saying it exactly enough, and of working out why you didn't. Tap "sort by price" and a list of a thousand items rearranges before you've lifted your finger. Somewhere, a person wrote down steps so exact that a machine with no judgment, no memory of yesterday and no idea what a "price" is could do it perfectly, a thousand times over. That is the whole trade. A programmer isn't someone who knows a language, they're someone who can state a task so precisely that nothing is left to guess. And since people are bad at being precise on the first try, the other half of the skill is finding out what you actually said, as opposed to what you meant. Mechanical and aerospace engineers use the same habits to simulate a wing before it's built; computer scientists use them to build every other subject on the degree. Python and MATLAB are the common first languages for engineers, and this piece uses Python-style code. The ideas carry to every language. The course in one box Who takes it: first-year students of Computer Science (usually their very first course) and Mechanical & Aerospace Engineering (often in MATLAB or Python). What it's for: telling a machine what to do exactly, and finding out when you didn't. Comes after: nothing; it has no prerequisites. Leads to: Data Structures (34/23), Computer Architecture (34/25) and Software Engineering (34/31) on the computing side; Differential Equations (34/03), Dynamics (34/12) and Control Systems (34/20) on the engineering side, all of which are solved on a computer in practice. Honest line: you won't be able to ship an app after this. You will be able to read ten lines of code, predict what they do, and make a physics model in a loop. A program is instructions plus state Strip away the syntax and a program is two things: a list of instructions, and the state, the current values of all its named boxes. Each instruction reads the state and changes it a little. Running the program is just that, one step after another, and the answer is whatever state is left (or printed) at the end. The graduate habit is to run the program in their head: keep a small table of variable values and update it line by line. Below is a tiny interpreter that does it visibly. Press Step. Then edit the program: change the 6 in the Collatz example, or write your own with =, +=, if/else, while, for … in range(…) and print. Interactive · Step through a program Add 1 to 5 Collatz steps Count down by 3 The program (edit me) StepRun to the endStart again A real mini-interpreter, written for this page. It stops after 400 steps, in case you write a loop that never ends. The first program adds 1 + 2 + 3 + 4 + 5 and prints 15. The Collatz one halves even numbers and triples-plus-one odd ones until it reaches 1: starting from 6 it takes 8 steps. Nobody knows if every starting number reaches 1, and yet the program is five lines long. That's what "simple instructions, deep behaviour" feels like. Variables and types: what a number is made of A variable is a name for a box in memory. A type says what kind of thing is in the box, and therefore what you may do with it: you can add two numbers, join two pieces of text, but adding "3" and 4 is an error or a surprise. Types exist because the machine stores everything as patterns of bits, and the pattern means nothing without the type. The box has a size too. A signed 16-bit integer holds −32,768 to 32,767, and nothing more. Decimals are stored as floating point, binary fractions of limited precision, so 0.1 can't be stored exactly and 0.1 + 0.2 gives 0.30000000000000004. On 4 June 1996 the first Ariane 5 rocket veered off course and broke up about 37 seconds after lift-off, a loss of more than US$370 million. The investigation traced it to the guidance software, reused from Ariane 4: a 64-bit floating-point value (the horizontal bias) was converted to a 16-bit signed integer, and on the faster rocket it didn't fit. The failure shut down both inertial units. A number too big for its box was enough. Interactive · Does it fit? 8 bits 16 bits 32 bits Number to store Ariane's code reportedly raised an error on overflow; some systems wrap silently instead, which is shown here. Conditionals and loops: choosing and repeating With state and instructions, only two more ideas are needed to compute anything: choose (if something is true, do this, else that) and repeat (while something is true, or for each value in a range). Together they let ten lines do work that would take ten million lines written out. The classic trap is in the edges. In Python's range(1, 10) the loop starts at 1 and stops before 10, so it runs 9 times, not 10. Off-by-one errors are the "fencepost" mistake: a fence of 10 sections needs 11 posts. Try the bounds yourself. Interactive · Edit the loop Add up i Add up i² Just count Start Stop (never reached) Step With the defaults, range(1, 11) runs 10 times and adds up to 55. Make the step negative and the start bigger than the stop to count down. Set the start above the stop with a positive step and the body never runs at all, which is perfectly legal and a common source of silent bugs. Functions and recursion: naming a piece, then calling it again A function packages some instructions under a name, takes inputs, and gives back an output. That is decomposition, the central skill of engineering of any kind: break a big problem into pieces small enough to understand and test alone, give each a name, then build on the names. The tracer above could be one function; so could braking distance. Recursion is a function that calls itself on a smaller version of its problem, until it reaches a case it can answer directly. The Fibonacci numbers (0, 1, 1, 2, 3, 5, 8, …, each the sum of the previous two) have a one-line recursive definition. Run it naively and the machine re-solves the same sub-problems again and again. Interactive · Count the calls Naive Remember answers (memo) Compute fib(n), n = Time it on this device Naively, fib(5) makes 15 calls, fib(20) makes 21,891, and fib(30), whose answer is just 832,040, makes 2,692,537. Remembering answers brings fib(30) down to 59 calls. Same answer, a different algorithm, and the course after this one (Algorithms, 34/24) is the study of exactly that difference. Data: lists and dictionaries Most programs aren't about single numbers but collections. A list is an ordered row of values, found by position (the third item). A dictionary maps a key to a value (a word to its count, a student ID to a name) and finds an entry by key without scanning. Choosing the right container is the first lesson of Data Structures (34/23). The difference is speed. To find a word in a list of n you may have to look at all n. A dictionary (usually built as a hash table) jumps almost straight to it, in roughly constant time on average, however large it is. Paste some text below and look a word up. Interactive · Count the words Textthe cat sat on the mat and the dog sat on the cat Look up a word Debugging is the scientific method Programs don't do what you meant; they do what you wrote. Finding the gap is science: look at the result, form a hypothesis ("the loop is skipping the last item"), make a prediction ("if I print the count here it will be one short"), run the experiment (print it), and revise. Guessing and changing things at random is the opposite, and it's how you add new bugs. The strongest move is bisection. If the code is a chain of steps and the result is wrong, check the value halfway along. If it's right, the fault is in the second half; if wrong, the first. Each check halves the suspects. In 1947 operators of Harvard's Mark II computer taped a moth found stuck in a relay into the logbook with the note "First actual case of bug being found" — a joke that only works because engineers already called faults "bugs". Interactive · Find the bug in as few probes as you can Start with x = 5. Slide to a step and probe the value there against what you'd get by hand. Look at x Probe hereBlame this stepNew bug Twelve suspect lines need at most four probes when you halve each time, because 2⁴ = 16 is at least 12. Checking from the top, one line at a time, can take twelve. Modelling: a physics simulation in a dozen lines Engineers model long before they build. Most physical laws give a rate of change (acceleration from forces, as in Physics I), and a computer can follow a rate forward in small steps: update the position using the velocity, update the velocity using the acceleration, repeat. That is the Euler method, and the loop is the whole simulation. x, y = 0.0, 0.0 vx, vy = 21.2, 21.2 # 30 m/s at 45 degrees dt = 0.1 # time step in seconds while y >= 0: x = x + vx * dt y = y + vy * dt vy = vy - 9.81 * dt # gravity changes vy only print(x) # about 93.9 metres Smaller time steps are more accurate and cost more. For the cannonball at 30 m/s and 45°, the exact range without air is 91.7 m. With a step of 0.5 s the simulation says 102.1 m, 0.1 s gives 93.9 m, and 0.001 s gives 91.8. Now add air drag, for which no tidy formula exists: the same loop still works, which is why engineers simulate. Interactive · A cannonball in code No air With air drag Launch speed Launch angle Time step Drag is modelled as a quadratic term with a baseball-like coefficient (0.007 per metre); a model, not a measurement. Readable code is a letter to a person A program is run by a machine once and read by people many times: the next engineer, the reviewer, you in six months. The well-known line from the preface of Structure and Interpretation of Computer Programs is: "Programs must be written for people to read, and only incidentally for machines to execute." The machine will accept either version below. Only one of them can be checked by a human. def f(a, b): return (a / 3.6) ** 2 / (2 * b * 9.81) def braking_distance(speed_kmh, friction): speed = speed_kmh / 3.6 # convert to metres per second return speed ** 2 / (2 * friction * 9.81) Both give 14.05 m for 50 km/h and a friction of 0.7, the braking distance of Physics I. Only the second tells you what it means, and so lets a reader spot a wrong constant. Good names, small functions and comments that say why are the cheapest debugging tools there are. People disagree on how much structure is enough: some teams favour heavy up-front design, others short cycles of build, test and revise; Software Engineering (34/31) covers both camps. What a graduate just sees A program is instructions plus state; to understand it, track the variables line by line. Every value has a type and a size, and a number too large for its box is a real way to lose a rocket. Loops run until their stop value, not through it: watch the edges. Break a problem into named functions you can test one at a time. A function that calls itself can be elegant and wildly wasteful; counting calls shows which. The right container (list or dictionary) can be the difference between instant and never. Debugging means hypotheses and experiments; halve the suspects each time. Simulation is a rate of change plus a loop, and smaller steps cost more for less error. Code is read far more often than it is written. Where it shows up next Data Structures (34/23) and Algorithms (34/24): the list-versus-dictionary and naive-versus-memo contrasts above, measured properly. Computer Architecture (34/25): the fixed-size boxes and overflow are the hardware's registers; instructions plus state is literally how a processor works. Software Engineering (34/31): decomposition, testing and readable code, scaled to teams and years. Differential Equations (34/03), Dynamics (34/12), Control Systems (34/20): the Euler loop is the starting point of numerical solvers, vehicle simulators and the code that runs a controller. Machine Learning (34/32): training is a loop that nudges numbers until the error falls. If you want to go deeper 📘 Standard textbook. John Guttag, Introduction to Computation and Programming Using Python (3rd ed., MIT Press, 2021), the text behind MIT's introductory course, written for engineers and scientists, with modelling and data. 📘 Gentler, free. Allen Downey, Think Python (3rd ed., Green Tea Press / O'Reilly), free to read online at greenteapress.com. For the deeper "why", Structure and Interpretation of Computer Programs by Abelson and Sussman is free online from MIT Press. 🌐 Free course. MIT OpenCourseWare 6.0001 Introduction to Computer Science and Programming in Python (Fall 2016), or Harvard's CS50x (starts in C, then Python and SQL) and CS50P (Python only), free on OpenCourseWare. The takeaway A program is instructions acting on state. Variables, types, choices, loops and functions are how you say it exactly; lists and dictionaries are how you hold data; debugging is hypothesis and experiment; and a simulation is a rate of change inside a loop. The graduate's real skill is reading code, predicting it, and taking a big problem apart. Try this: take any small thing you do by hand, like splitting a bill or converting recipe units, and write it as exact steps for a friend to follow without any judgment. Then watch where they get stuck: those are your bugs. Sources & Further Reading Sources. Ariane 5 flight 501 (4 June 1996; 64-bit floating-point to 16-bit signed integer conversion overflow; about 37 seconds; more than US$370 million): Wikipedia, "Ariane flight V88", which summarises the Inquiry Board's report. The Mark II moth (9 September 1947, logbook note "First actual case of bug being found", now at the Smithsonian) from the Computer History Museum and JSTOR Daily, "The bug in the computer bug story", which notes the term was already in use. The SICP quotation from the preface of Abelson & Sussman, Structure and Interpretation of Computer Programs (MIT Press, 1985). CS50x language order from cs50.harvard.edu/x; Think Python 3rd ed. from greenteapress.com. Computed and checked in node, including by running the page's own interpreter on its example programs: 15; 8 Collatz steps; −2; 40,000 wraps to −25,536 in 16 bits; fib call counts 15, 177, 21,891 and 2,692,537 (and 2n−1 with a memo); projectile ranges 91.74 m exact and 102.1, 93.9, 91.8 m for steps of 0.5, 0.1 and 0.001 s; 14.05 m braking distance. 🌐 MIT OpenCourseWare: 6.0001 🌐 Harvard CS50x and CS50P 📘 Think Python, 3rd edition (free online) 📘 Guttag, Introduction to Computation and Programming Using Python, 3rd ed. (MIT Press) 🌐 JSTOR Daily: the bug in the computer bug story Step 5. Engineering Drawing & CAD: The Contract Between Designer and Maker -------------------------------------------------------------------------- Why: Drawings are the contract between engineer and builder. A drawing is not a picture of a part. It is a promise about what will be made, how exactly, and who is to blame when it isn't. Buy a shelf bracket, and the screw holes line up with the studs in your wall. Buy a replacement bearing for a bicycle, a car, a washing machine, and it fits a shaft that was made in a different factory, on a different continent, years earlier. Nobody phoned anyone. What made all that work is a sheet of paper, or these days a file, that says: make it exactly like this, to within this much. That sheet is an engineering drawing, and it is the quietest important idea in mechanical engineering. A graduate looks at one and reads it the way you read a sentence: the views, the numbers, the little symbols. The skill sounds dull until you see how many ways a part can come back wrong when a drawing is vague. The course in one box Who takes it: mechanical, aerospace and civil engineers, and anyone in product design, usually in year 1. What it's for: saying precisely what shape a part is, how exact it must be, and then turning that into a model that can be analysed and made. Comes after: nothing; it needs only geometry. Leads to: Machine Design (34/19), where every part you design is communicated by a drawing; it also supports Materials Science & Manufacturing (34/13), Mechanics of Materials (34/11) and the Capstone (34/22). Honest line: you won't learn a CAD package here (those change every decade); you will learn what every package is asking you to say. A drawing is a contract A drawing is not an illustration. It is the document the designer hands to the maker, and it is what the buyer pays against. If the maker builds exactly what is drawn and it doesn't work, the question "whose fault is that?" is settled by reading the drawing. If the maker builds something different and it fails, that is settled by the same document. This is why engineers insist on one rule: if it matters, it is on the drawing. Not in an email, not in a phone call, not "everyone knows". A change that isn't on the drawing doesn't exist, and a change that is on it has to be checked by someone who understands what it does. The sad proof is the Hyatt Regency walkways in Kansas City, which collapsed in 1981 and killed 114 people (some sources say 113). The contract drawings showed the hanger rods holding two walkways as single long rods. The steel fabricator's shop drawings changed that to two sets of shorter rods, which roughly doubled the load on the critical connection. The US National Bureau of Standards investigation found two problems: the original connection was already inadequate (it would have reached only about 60% of the capacity the building code expected), and the change made it far worse. The lesson isn't that someone was careless with a pencil. It is that a drawing carries a decision, and an unchecked change to a drawing is an unchecked change to a building. A drawing says what to make, how exactly, and who is responsible when it isn't. Three flat views tell you a solid A 3D object has to be described on flat paper. The method used for roughly two centuries is orthographic projection: imagine the object inside a glass box, and look straight at each face from far away. Projection lines are parallel, there is no perspective, and so a length on the drawing is a true length: parallel edges stay parallel, and a hole that is 12 mm across is drawn 12 mm across in the view where you look down its axis. Usually three views are enough: front, top and right side. The front and top share the same widths, the front and side share the same heights, and from any two you can find the third. Edges you can't see are drawn dashed. If the inside matters, you imagine slicing the part and drawing the cut face with hatching, a section view. Several parts together appear in an assembly drawing, often with an exploded view and a numbered parts list. There are two conventions for where the views go on the sheet. Third angle (used in the US and Canada) puts the top view above the front. First angle (used in most of Europe and Asia) puts it below. Same views, different layout; a small symbol in the title block says which one you are reading. Try it below: the three views come from a real 3D model, projected live, with hidden lines worked out by ray-casting. Interactive · Rotate the part, read its three views Hidden lines Third angle Turn Tilt The part is an L-bracket, 60 × 40 × 40 mm, with a 12 mm vertical hole in the base, a 10 mm horizontal hole through the upright and a notch in one corner. Solid lines are visible edges and dashed lines are hidden ones. Drag the picture to rotate it. Every dimension is a range: tolerance and stack-up No machine makes anything exactly 25.000 mm, and nothing needs to. What matters is how far off is still acceptable. A tolerance is that allowed range: 25 ± 0.1 mm means anything from 24.9 to 25.1 is good. Tight tolerances cost money (more careful machines, more inspection, more scrap), so the designer's job is to give each dimension the loosest tolerance that still works. The trouble is that tolerances add up. Put three parts in a housing: the gap left over is the housing size minus all three parts, so every ± feeds into the gap. Adding the extremes gives the worst case, which is safe but costly. But it is very unlikely that every part lands at its extreme in the same direction. If you treat each ± as three standard deviations of a bell curve, the combined spread is the square root of the sum of the squares, the statistical (root-sum-square) stack: worst case = Σ ti    statistical = √(Σ ti²) Worked number: a 100 mm housing holds parts of 30, 40 and 28 mm, leaving a 2 mm gap. With tolerances ±0.30, ±0.20, ±0.20 and ±0.15 mm, worst case is ±0.85 mm (gap 1.15 to 2.85 mm). The statistical range is ±0.44 mm. That 2:1 difference is why statistical stacks let designers loosen tolerances; it is also a bet that your manufacturing really is centred and random. Interactive · Tolerance stack-up calculator Housing (100 mm) Part 1 (30 mm) Part 2 (40 mm) Part 3 (28 mm) The statistical figure assumes each tolerance is three standard deviations of a normal distribution and the parts are independent. Check this against real inspection data before relying on it. One more everyday consequence: how you dimension decides how tolerances add. Dimension a row of holes as a chain (each from the last) and the errors accumulate along the chain. Dimension them all from one datum edge and each hole only inherits its own tolerance. Good drawings dimension from the surfaces that matter, and never make the machinist do arithmetic. Fits, and GD&T in plain words When a shaft goes into a hole, what matters is the gap between them. A fit is a standard pairing of tolerances chosen for the job. ISO 286 names them with letters and numbers: a capital letter for the hole, a lowercase one for the shaft, and a number for the tolerance grade. For a 25 mm size, H7 for the hole means 25.000 to 25.021 mm. Pair it with g6 and the shaft runs freely; with k6 you may get either a small gap or a small squeeze (a transition fit); with p6 you must press the shaft in (an interference fit), and the parts then grip by friction alone. Interactive · Pick a fit H7/g6 H7/h6 H7/k6 H7/p6 Deviations are for the 18 to 30 mm size band from ISO 286 tables (H7: 0 to +21 µm, g6: −20 to −7, h6: −13 to 0, k6: +2 to +15, p6: +22 to +35). A micrometre is a thousandth of a millimetre; a human hair is very roughly 50 to 100 of them. A ± size tolerance only controls the size at each point. It says nothing about whether the hole is straight, round, square to the face, or in the right place. Geometric dimensioning and tolerancing (GD&T) is a language of symbols for exactly that: flatness, perpendicularity, roundness, and above all position. Each rule reads like a short sentence: "this hole's axis must lie within a 0.2 mm circle, measured from datum surfaces A, B and C." The datums are the three flat faces that sit the part on a table, against one wall, and against another; they are how a maker and an inspector agree on where "measure from" is. A circular position zone also holds more good parts than a square ± box with the same corner reach: a ±0.1 mm square has area 0.040 mm², and the circle through its corners (0.283 mm across) 0.063 mm², about 57% more. In the US the rulebook is ASME Y14.5 (the 2018 edition, reaffirmed in 2024); ISO has its own set. Parametric CAD: you draw the intent, not the shape Computer-aided design began as an electronic drafting board. Ivan Sutherland's Sketchpad, in the early 1960s, let a person draw with a light pen and have the computer keep the lines connected. Modern CAD went further: since the late 1980s (PTC's Pro/ENGINEER shipped in 1988), mainstream CAD records how you built the shape. You draw a 2D sketch, add constraints (this line is horizontal, that circle is 12 mm, these two are equal), then extrude or cut it into a solid. Every step goes into a feature tree and every number is a parameter, so changing one rebuilds the part. The point is design intent: the rules you chose say what must stay true when something changes. Widen this plate: should the hole stay 30 mm from the wall, or stay in the middle? Both are legitimate, and they are different parts. Try it. Interactive · A sketch that keeps its promises Hole 30 mm from the left Hole stays centred Plate width How many constraints does a sketch need? Tap constraints on and off. The plate sketch is a rectangle (8 coordinates) and a circle (3 numbers) = 11 degrees of freedom. Each constraint removes some. From the model to simulation and the shop floor Once a part exists as a precise digital solid, two things become cheap. First, you can simulate it. The standard tool is the finite element method (FEA): chop the part into thousands of little elements, write the equations of stiffness for each one, and solve all of them together for how the part deforms and where the stress is. It is a numerical answer, and the discipline is to ask how wrong is it? The way to find out is to refine the mesh and see whether the answer settles. The bar below is a real, tiny finite-element model: a tapered aluminium bar, pulled with 10 kN at the thin end. Interactive · Refine the mesh Number of elements Bar: 100 mm long, cross-section tapering from 200 mm² at the wall to 50 mm² at the free end, E = 70 GPa, 10 kN axial load. Each element has linear displacement, so its stress is one constant number. The exact answers come from integrating the load over the taper. Look at the two errors. Movement settles quickly; stress, which depends on the slope of the movement, converges much more slowly, and the worst stress is usually where it matters most. A graduate never trusts a coloured stress plot with one mesh. Second, the same model feeds manufacturing. CAM (computer-aided manufacturing) turns the solid into toolpaths: a list of moves for a CNC machine (often G-code) to cut the part from stock. For 3D printing the model is sliced into layers and each layer into a toolpath. The drawing doesn't disappear: the tolerances, fits and datums are still what the inspector checks, even if the model carries them. What a graduate just sees A drawing is a contract: if it matters, it is on the drawing, and a change needs a check. Any object can be read from two or three orthographic views; hidden features are dashed. Every dimension is a range, and tolerances cost money. Tolerances stack: worst case is the sum, statistical is the root-sum-square, and datum dimensioning stops errors accumulating. A fit is a pairing of two tolerances, and "clearance, transition, interference" is a ladder from sliding to pressed. A size tolerance doesn't say the part is straight or in the right place; GD&T does. A CAD model is a record of intent, with sketch, constraints and a feature tree; what you constrain is what stays true. A simulation is a number with an error bar, and you check it by refining the mesh. Where it shows up next Machine Design (34/19): every shaft, gear and bolted joint you design ends up as a dimensioned drawing, with fits for the bearing seats. Materials Science & Manufacturing (34/13): the process (casting, machining, printing) sets the tolerances you can realistically ask for. Mechanics of Materials (34/11): FEA is the numerical cousin of the beam and stress formulas. Systems Engineering & the Capstone (34/22): interface drawings are how subsystems built by different teams agree to fit. In the real world: product design, a car's body panels fitting with an even gap, an aircraft's thousands of parts matching across suppliers. If you want to go deeper 📘 Standard texts: Giesecke and co-authors, Technical Drawing with Engineering Graphics (Pearson, 15th edition; a 16th is listed), the usual first-year text. Gene Cogorno, Geometric Dimensioning and Tolerancing for Mechanical Design (McGraw Hill, 3rd edition), for GD&T. The standard itself is ASME Y14.5-2018. 📘 Gentler: the free MIT OpenCourseWare short course How to CAD Almost Anything (RES.16-002), a hands-on introduction to parametric modelling through reverse engineering. 🌐 Free course: NPTEL: Engineering Graphics and Design (IIT Delhi, Profs. Naresh Datla and S. R. Kale). The takeaway A drawing is the contract between a designer and a maker: views to say the shape, dimensions and tolerances to say how exactly, fits and GD&T to say how parts meet. Tolerances add up, so you stack them deliberately. CAD records your intent so changes rebuild correctly, and simulation and manufacturing both read the same model, which is why you check each by asking how wrong it could be. Try this: pick any small object with a few parts (a hinge, a jar lid) and sketch its front, top and side views. Then guess which two dimensions matter most for it to work, and what a ±1 mm error on each would do. Sources & Further Reading Sources. Hyatt Regency collapse: US National Bureau of Standards, NBSIR 82-2465 (1982) and its NIST-hosted summaries; the Linda Hall Library notes (114 deaths; NIST pages and UPI report 113). Third/first angle: ISO 128 and ASME Y14.3 as summarised by Xometry and Peachpit. ISO 286 fits for 18 to 30 mm: ISO 286-2 values as reproduced by cadforum.cz (H7, g6, h6, k6, p6). ASME Y14.5-2018 (R2024): ASME and ANSI webstore listings. Pro/ENGINEER commercial shipments from January 1988: PTC history summaries. Sketchpad (Sutherland, early 1960s): PTC/NovEdge history summaries. Numbers computed in node: the three views and hidden lines by ray-casting; stack-up (worst case 0.85, statistical 0.439); the p6 and g6 limits; position-zone areas (0.040 and 0.063 mm²); and the tapered-bar finite-element model checked against its exact integral (exact tip movement 0.13203 mm). 🌐 NPTEL: Engineering Graphics and Design 🌐 MIT OpenCourseWare: How to CAD Almost Anything 📘 Giesecke et al., Technical Drawing with Engineering Graphics 📘 Cogorno, Geometric Dimensioning and Tolerancing for Mechanical Design Step 6. Multivariable Calculus: Calculus on a Landscape ------------------------------------------------------- Why: Wind, water and soil pressure vary over surfaces and volumes. A raindrop on a hillside picks the steepest way down without thinking. The arrow it follows is the same one that trains every neural network and drives heat through a metal plate. Stand on a hillside with a map. Height depends on two things at once: how far east you are and how far north. Single-variable calculus handles one input at a time. The hillside needs a new question: of all the directions I could step, which one rises fastest? Almost nothing real depends on one number. Temperature in an engine block depends on position in three dimensions. The error of a language model depends on billions of adjustable weights. Air pressure depends on place and time. This course takes the slope-and-area ideas of Calculus I (34/01) and lets them live in more than one dimension. Surprisingly little new is needed: one new arrow, the gradient, and one new accounting rule, the idea that what flows out of a region is what is made inside it. The course in one box Who takes it: second-year Mechanical & Aerospace Engineering students; Computer Science students meet it in year 4 when machine learning needs it. What it is for: calculating change and totals when many things vary at once. Comes after: Calculus I (34/01). Leads to: Physics II: Electricity & Magnetism (34/07), Fluid Mechanics (34/15), Machine Learning (34/32), and through them Heat Transfer (34/16) and Propulsion (34/21). Honest line: you will not be able to grind through a surface integral after this. You will see why a gradient points uphill, why training an AI is "just" rolling downhill, and what the big theorems say in plain words. Two inputs make a landscape, and a partial derivative is one slope A function of two inputs is a surface: height above each point of a map. Engineers draw it the way hikers do, with contour lines: closely spaced lines mean steep ground, widely spaced lines mean gentle ground. To ask "how steep?" you must say in which direction. Walk due east and note the slope: that is the partial derivative with respect to x, written ∂f/∂x. Hold east still and walk north: that is ∂f/∂y. The recipe is simple: treat every other variable as frozen and differentiate as before. ∂f/∂x = slope of f as only x changes, with y held still The landscape below has a tall hill near (1, 1) and a lower one near (−1.2, −0.8). Drag the point on the map (or use the sliders). The cyan arrow shows the slope eastward and the green arrow the slope northward; a leftward or downward arrow means the ground falls in that direction. Interactive · Walk the map East–west position x North–south position y At the starting point (−0.2, 0.3) the height is 0.684, the eastward slope +0.621 and the northward slope +0.226: the ground rises both ways, more steeply towards the east. The gradient points straight uphill Put the two partial slopes together as one arrow, east component first: that is the gradient, written ∇f. It has a beautiful property: it points in the direction of steepest ascent, and its length is how steep that is. Step in any other direction and the slope is the gradient's length times the cosine of the angle between your step and the arrow. Straight along it, cosine is 1: the most you can get. Sideways to it, cosine is 0: you walk along a contour and stay level. slope in direction u = ∇f · u = |∇f| cos(angle) Two consequences follow. The gradient is always perpendicular to the contour line through that point, which is why water runs across contours, not along them. And the same arrow governs heat: in Fourier's law, heat flows down the temperature gradient, from hot to cold, at a rate proportional to its steepness. Many diffusion laws have this shape: flow follows the negative gradient. Interactive · Which way is uphill? Direction you step (0° = east, 90° = north) Turn to the gradient Turn the pink direction until the slope is largest: it lines up with the orange arrow. Turn it 90° from the arrow and the slope is zero. At the very top of a hill the gradient shrinks to nothing, which is the multivariable version of "zero slope means a peak" from Calculus I. Rolling downhill is how a network learns Now flip the sign. A machine-learning model has adjustable numbers (weights), and a loss measures how wrong it is on the training data. The loss is a landscape over the weights, with one dimension per weight. To learn, take a small step against the gradient, so downhill, then recompute and repeat. This is gradient descent, which Augustin-Louis Cauchy proposed in 1847 for solving equations, long before computers. The chain rule from Calculus I is what lets the gradient be computed layer by layer, in the method called backpropagation, published by Rumelhart, Hinton and Williams in Nature in 1986. A model like GPT-3 (2020) has 175 billion weights, so its landscape has 175 billion dimensions, yet the rule is the same arrow. new position = old position − (step size) × ∇(loss) The step size matters, and you can feel why. In the bowl below, one direction is twelve times steeper than the other. Too small a step crawls; just right is fast; too large and each step overshoots the valley floor by more than it started and the loss explodes. For this bowl, steps above 1/6 (about 0.167) diverge. Interactive · Learning rate Too small Good Zig-zag Too big Step size (learning rate) Steps taken Real training uses refinements (momentum, adaptive step sizes, random mini-batches), but all of them are ways of coping with exactly this picture: a long narrow valley where the steepest direction is not the direction of the destination. Double and triple integrals add up a landscape's worth Integrals extend the same way. A single integral adds thin slices along a line. A double integral adds tiny patches over an area: if each patch has a density (kilograms per square metre), the sum is the plate's mass. A triple integral does the same over a solid, with density in kilograms per cubic metre. The recipe is slice, multiply, add, repeated once per dimension. mass = sum over all pieces of (density × piece size)  =  ∬ ρ dA The same sums also give the centroid, the balance point: weight each position by its piece's mass and divide by the total. Engineers use it constantly: the centre of mass of a rocket, a car or an aircraft decides whether it is stable. Try a plate 2 m long and 1 m wide whose density rises along the length. Interactive · Where does the plate balance? Density rises along the length: a (kg per m² per m) Density rises across the width: b The density here is 1 + ax + by kg/m². With a = 3 and b = 0, the plate's mass is 2 + 2×3 = 8 kg and it balances at x = 1.25 m, shifted from the middle (1 m) towards the heavy end. A vector field puts an arrow on every point Wind on a weather map, water in a pipe, the pull of a magnet: each attaches an arrow (a direction and strength) to every point in space. That is a vector field. Two numbers computed from the field's partial derivatives tell you most of what it is doing locally. Divergence: is the field spreading out from here, or converging? A tap under water is a source; a drain is a sink. A positive divergence means a small blob of fluid is expanding. Curl: is the field spinning? Put a tiny paddle wheel in the flow. If it turns, the curl is not zero, and the sign tells you which way. divergence = ∂vx/∂x + ∂vy/∂y    curl = ∂vy/∂x − ∂vx/∂y Pick a field and let a round blob of fluid (grey dashes show where it started) be carried along. You will see that a source inflates the blob without turning it, a whirlpool turns it without changing its area, and a saddle stretches it one way and squeezes it the other while keeping its area. Interactive · Drift a blob through a field Source Whirlpool Shear Saddle Time elapsed Divergence and curl are the two words in Maxwell's equations of electromagnetism (Physics II, 34/07) and in the equations of fluid flow (Fluid Mechanics, 34/15). Spreading and spinning are the two things a field can do locally, besides moving along as a whole. What flows out of a region is what is made inside Here is the idea behind the three big theorems, in the language of a bathtub. If water pours out through the boundary of a region, then some taps inside must be producing it (or the level inside must be falling). Add up every tiny source and sink inside and you get the net flow out through the boundary. Nothing more is needed. net flow out through the boundary = total divergence inside That is the divergence theorem, the three-dimensional version of the statement. Green's theorem is the flat version (the name honours George Green, whose Essay of 1828 contained the related surface-to-volume formula), and Stokes's theorem relates the spinning of a field inside a surface to the flow around its edge. All three say the same kind of thing: a sum over a boundary equals a sum of local changes inside. It is why one can measure a rotating machine from outside and why conservation laws can be written either as "this region's total changes by what crosses its walls" (integral form) or as an equation about each point (differential form). Interactive · Count the flow two ways Source Cubic Source in, sink out Steady wind Circle radius Circle centre, left or right The first number adds up the field's push through the circle itself. The second never looks at the circle's edge: it adds the divergence in every patch inside. They agree for every field, radius and position. For the plain source (x, y), the divergence is 2 everywhere, so a unit circle encloses 2 × π = 6.283 units of production, and exactly 6.283 flows out through its rim. What a graduate just sees A slope needs a direction; a partial derivative is the slope with every other input held still. The gradient points straight uphill, its length is the steepness, and it is perpendicular to the contour. "Downhill" is the whole idea behind training a model: step against the gradient, and mind the step size. Heat, and many other things, flow down a gradient from high to low. A double or triple integral is slice, multiply, add, once per dimension; mass and balance point come from it. A field can only spread (divergence) or spin (curl), locally. What leaves a region through its wall is exactly what its insides produce. Where it shows up next Physics II: Electricity & Magnetism (34/07): Maxwell's equations are statements about the divergence and curl of the electric and magnetic fields. Fluid Mechanics (34/15): conservation of mass is a divergence equation, and the divergence theorem turns "what crosses the walls of this control volume" into "what happens inside it". Heat Transfer (34/16): conduction is heat flowing down the temperature gradient; the heat equation is the gradient and divergence together. Machine Learning (34/32): training is gradient descent on a loss landscape with millions or billions of dimensions. Mechanics of Materials (34/11) and Dynamics (34/12): centroids and moments of inertia are integrals of exactly the kind you added up above. If you want to go deeper 📘 Textbook: Stewart, Clegg and Watson, Calculus: Early Transcendentals (9th ed.), chapters 14 to 16: the standard multivariable text, rich in examples. 📘 Textbook: Jerrold Marsden and Anthony Tromba, Vector Calculus (6th ed., W. H. Freeman): the classic for the geometry of divergence, curl and the big theorems. ▶️ Gentler: Khan Academy's Multivariable Calculus course, which has videos with a strong visual feel, and Grant Sanderson's multivariable lessons. 🌐 Free course: MIT OpenCourseWare 18.02SC Multivariable Calculus, with lectures, notes and exercises. The takeaway With many inputs, a slope needs a direction. The partial derivatives are the slopes along the axes; the gradient joins them into one arrow that points uphill and measures how steep it is. Walk against it and you are doing gradient descent, which trains models with billions of weights and also describes how heat flows. Integrals extend to areas and volumes, giving mass and balance points, and for flows the whole subject rests on one bookkeeping rule: what leaves a region through its boundary is what its interior produces. Try this: on any hill, a park slope or a street, find the line along which a ball would roll: that is the negative gradient. Then notice that the ball's path crosses the contours at right angles. Sources & Further Reading Sources. Gradient descent: A.-L. Cauchy, "Méthode générale pour la résolution des systèmes d'équations simultanées", Comptes rendus 25 (1847), pp. 536–538; Claude Lemaréchal, "Cauchy and the gradient method", Documenta Mathematica, Extra Volume ISMP (2012). Backpropagation: Rumelhart, Hinton and Williams, Nature 323 (1986), pp. 533–536. GPT-3's 175 billion parameters: Brown et al., "Language Models are Few-Shot Learners" (2020, arXiv 2005.14165). Green's Essay of 1828 containing the theorem named after him: MacTutor biography and the essay itself. Computed in node: the landscape f = 2e−((x−1)²+(y−1)²)/1.5 + 1.2e−((x+1.2)²+(y+0.8)²)/1.0 and its analytic gradient (checked against finite differences, agreement to 3 × 10−10); gradient descent on ½(x² + 12y²), stable for step sizes below 2/12; plate mass and centroid by 3,200-cell sums against the exact formulas; blob area against ediv·t; circle flux against the divergence integral for four fields (agreement to 4 decimals). 🌐 MIT OpenCourseWare 18.02SC Multivariable Calculus 🌐 Khan Academy: Multivariable calculus 🌐 Documenta Mathematica, ISMP volume, which contains Lemaréchal's short history of gradient descent 📘 Stewart, Clegg and Watson, Calculus: Early Transcendentals, 9th ed. (Cengage) 📘 Marsden and Tromba, Vector Calculus, 6th ed. Step 7. Differential Equations: Laws Written as Rates ----------------------------------------------------- Why: Settlement, flow and cooling are differential equations. A cup of coffee cools, a capacitor charges, a car bounces over a bump, a fish stock collapses. Each follows a rule that says how fast things change, not where they will be. Solving it is the whole course. Pour coffee and watch it cool. It loses heat fast at first, when it is far hotter than the room, and ever more slowly as it approaches room temperature. Nobody wrote down the temperature at every minute in advance. All you need is a rule: the hotter the cup is above the room, the faster it cools. That sentence, turned into an equation, is a differential equation, and it predicts the whole curve. This is the secret structure of physics. Nature rarely says "the position will be this". It says "the rate of change is that, given the state right now". Newton's second law is such a rule. So are the equations of electrical circuits, chemical reactions, populations, heat, and the flow of air over a wing. A differential equations course teaches how to read these rules and how to turn them into predictions, either exactly or step by step on a computer. The course in one box Who takes it: second-year Mechanical & Aerospace Engineering students. It is not on the Computer Science route in this series, but CS students meet the same ideas whenever they simulate, animate or train anything that evolves in time. What it is for: predicting how a system evolves when you only know the rule for its rate of change. Comes after: Calculus I (34/01). Leads to: Dynamics (34/12), Fluid Mechanics (34/15), Heat Transfer (34/16), Vibrations (34/17), Control Systems (34/20). Honest line: you will not learn the solution tricks that fill the textbook. You will see the three behaviours that recur everywhere (settling, oscillating, tipping) and why a simulation can lie to you. An equation about a rate hands you a direction at every point An ordinary differential equation, an ODE, is an equation that mentions a quantity and how fast it is changing. The simplest useful kind says the rate depends on how far the quantity is from a target: dy/dt = k (target − y) Read it aloud: the further I am from the target, the faster I move towards it. Coffee towards the room's temperature, a capacitor towards the supply voltage, a medicine level towards zero. The picture is a slope field: at every point of the plane the equation fixes the direction a solution must go through, drawn as a tiny segment. A solution is any curve that follows those segments the whole way, like a leaf carried along a stream. Different starting points give different curves, all flowing toward the target. Interactive · Follow the arrows Pull strength k Target value Starting value Notice that the segments get flatter as you approach the dashed line: the rate shrinks with the distance. That is why the curve approaches the target but never quite arrives. The exact solution is the target plus the starting gap times a decaying exponential, e−kt. The time constant: how fast a first-order system forgets The rule dy/dt = k(target − y) has one number that sums it up: the time constant, τ = 1/k. In one τ the system covers about 63.2% of the remaining gap; after about 4.6 τ it has covered 99%. Engineers read response speed from that single number. fraction of the gap covered after time t = 1 − e−t/τ In an electrical circuit, a resistor R in series with a capacitor C charges towards the supply with τ = R × C. A 10 kΩ resistor and a 100 µF capacitor give τ = 10,000 Ω × 0.0001 F = 1 second: after 1 s the capacitor holds 63.2% of the supply (3.16 V of 5 V), after 4.6 s, 99%. A cooling coffee follows the same curve; only the story and the units change. Interactive · Charge a capacitor Resistance R Capacitance C Notice what the sliders do: double either one and the curve stretches twice as wide, but keeps exactly the same shape. The same shape is why a thermometer takes a while to read a fever, and why a control engineer (Control Systems, 34/20) asks first: "what is the time constant?" Solving it step by step: Euler's method Most differential equations have no neat formula. The fall-back is to walk along the slope field in small steps. At the current point read the slope, move a short time h in that direction, and repeat. This is Euler's method, set out by Leonhard Euler in his Institutionum calculi integralis (1768), and it is the ancestor of every physics engine, weather model and circuit simulator. ynext = y + h × (the rate at y) Test it on a rule we can also solve exactly, dy/dt = −y, starting from 1. Each Euler step multiplies the value by (1 − h). Small steps follow the true curve e−t closely. But raise h and something stranger happens. Interactive · Euler versus the exact answer Step size h At h = 1 the first step lands exactly on zero. Beyond 1 the steps overshoot, and past h = 2 each step makes the number bigger in size, so the "answer" explodes even though the real one decays. That is numerical instability: a stepping rule has its own stability limit, and a simulation that ignores it will happily print nonsense. At h = 0.5 Euler gets 0.000244 at t = 6, where the true value is 0.00248: right shape, ten times too small. Smaller steps, or smarter methods (the Runge–Kutta family), fix that. Second order: the spring, the mass and the damper Newton's second law says force equals mass times acceleration, and acceleration is the rate of change of the rate of change of position. So a mass on a spring obeys a second-order equation: the second derivative of position x depends on position and velocity. m x″ + c x′ + k x = 0 The three terms are inertia (mass m), a damper that resists velocity (c) and a spring that pulls back with force proportional to displacement (k). Two numbers decide everything: the natural frequency, ω = √(k/m), and the damping ratio, ζ = c ÷ (2√(km)). ζ below 1, underdamped: it rings, swinging past rest and fading. ζ = 1, critically damped: the fastest return without overshoot. ζ above 1, overdamped: no swing, but a sluggish crawl home. Engineers tune for it. A door closer is aimed near critical damping so the door shuts without banging or bouncing. A shock absorber is a compromise: enough damping to stop the bouncing, not so much that the wheel can't follow the road. Interactive · Tune the damper No damping Light Critical Heavy Mass m Spring stiffness k Damping c With m = 1 kg and k = 25 N/m the natural frequency is 5 rad/s, or about 0.8 Hz. Setting c = 10 N·s/m gives ζ = 1, critical. With c = 2, ζ = 0.2 and the first swing past rest reaches 52.7% of the starting displacement on the other side. This same equation, with different symbols, describes the RLC circuit, a building in an earthquake and a robot arm: that is the next courses (Vibrations, 34/17, and Control Systems, 34/20). Equilibria: where things settle, and where they tip Often you do not need the whole curve, only the places where change stops: the equilibria, where the rate equals zero. Two questions follow. Where are they? And is each one stable (nudge it and it comes back) or unstable (nudge it and it runs away)? The answer is the slope of the rate curve at that point, from Calculus I: negative slope means stable, positive means unstable. A toy model shows why this matters. Fish regrow according to dy/dt = y(1 − y/100) per year (a standard logistic model), and a fleet takes a fixed catch H every year. For H = 20 there are two equilibria, at about 27.6 and 72.4 fish: the upper is stable, the lower unstable. Start just above 27.6 and the stock recovers to 72.4. Start just below it and the stock collapses to zero in a few years. Raise the catch above 25, the most the stock can regrow, and no equilibrium remains at all. Interactive · A fishery on the edge Catch 20, start 30 Catch 20, start 25 Catch 26 Fixed catch H per year Starting population The lesson generalises: a system can look healthy while it sits near an unstable tipping point, and a slow change in a parameter (here the catch) can make the stable state vanish suddenly. Bridges, power grids, aircraft and ecosystems are all analysed this way. The model is a classroom toy, not a description of any real fishery. Why physics is written in differential equations Why do laws come as rates? Because change is local. What a system does next depends only on what is happening here, now: the forces, the gradients, the current state. Newton's second law (second order in position), Maxwell's equations for electric and magnetic fields, the equations of heat flow and the equations of fluid motion all have this shape. Give the rule and the starting conditions and the future follows. Two kinds appear. An ordinary equation (ODE) has one changing variable, usually time: a pendulum, a circuit, a rocket's trajectory. A partial differential equation (PDE) has several, such as time and position: the temperature along a rod, the pressure in the air around a wing. The Navier–Stokes equations of fluid flow (Fluid Mechanics, 34/15) are PDEs so hard that the Clay Mathematics Institute named a question about them one of its seven Millennium Prize problems in 2000. Most real ones are solved on computers using the stepping idea above. A glimpse of PDEs: the heat equation Joseph Fourier wrote the heat equation in his Théorie analytique de la chaleur (1822). It says: the temperature at a point changes at a rate proportional to how much its surroundings, on average, are hotter or cooler than it is, which is the curvature of the temperature curve. ∂u/∂t = α ∂²u/∂x² A peak has negative curvature, so it cools; a dip has positive curvature, so it warms. Heat flows from hot to cold and the profile smooths out. The rod below is insulated at both ends, so the average temperature, the total heat, never changes, only its distribution. Fourier's great discovery (see the library's piece on the Fourier transform, 12/12) is that smooth wiggles of different sharpness die away at different rates: a wiggle five times sharper dies 25 times faster, because it decays like e−(n²π²)t. Interactive · Heat spreading along a rod Hot spot Half hot Gentle wave Sharp wave Time Time step, as a fraction of the stability limit Try the two waves at the same time. The gentle wave 1 + 0.8 cos(πx) has faded to 0.298 high at t = 0.1, matching the exact 0.8e−π²t. The sharp wave 1 + 0.8 cos(5πx) is gone almost at once, 25 times faster in the exponent. Then drag the time-step slider above 0.5 and watch the simulation explode: Euler's lesson again, now in a PDE. Real heat solvers (Heat Transfer, 34/16) must respect the same limit or use methods built to avoid it. What a graduate just sees A law about a rate is a direction at every point; a solution is a curve that follows those directions. First-order systems approach their target by a time constant: 63% in one τ, 99% in about 4.6 τ. Cooling, charging and decay are the same curve with different labels. A step-by-step simulation has a stability limit: too big a step creates growth that is not in the physics. Mass, spring and damper give three behaviours: ring, crawl, or the sweet spot at ζ = 1. Equilibria matter, and slope tells stable from unstable. Systems tip. Sharp features smooth out first in diffusion, and total heat is conserved. Where it shows up next Dynamics (34/12): every motion problem becomes an equation of the form above. Vibrations (34/17): the damped spring is the building block of everything that shakes, from engines to aircraft wings. Control Systems (34/20): a thermostat, cruise control or autopilot is a differential equation you shape with feedback; time constants and damping ratios are its vocabulary. Heat Transfer (34/16) and Fluid Mechanics (34/15): the heat equation and Navier–Stokes are PDEs; the stepping and stability ideas carry straight over. Computer science: physics engines in games, climate and weather models, and circuit simulators are loops of Euler-like steps; neural ODEs even treat a network as a flow. If you want to go deeper 📘 Textbook: William Boyce, Richard DiPrima and Douglas Meade, Elementary Differential Equations and Boundary Value Problems (12th ed., Wiley): the classic course text. 📘 Textbook: Steven Strogatz, Nonlinear Dynamics and Chaos (3rd ed., CRC Press): an unusually readable treatment of equilibria, stability and tipping. ▶️ Gentler: 3Blue1Brown's differential equations series, and Strogatz's Infinite Powers for the story of why calculus describes nature. 🌐 Free course: MIT OpenCourseWare 18.03SC Differential Equations, with video lectures, notes and problem sets. The takeaway Physics speaks in rates. A differential equation gives a direction at every point and a solution follows it. First-order systems approach a target with a time constant (63% in one τ). Second-order systems, from springs and masses, ring, crawl or sit at the sweet spot of critical damping. Equilibria are where change stops and slope tells stable from unstable. When no formula exists, you step along on a computer, but only small enough steps tell the truth: past a limit, a simulation blows up even though nature does not. Try this: take a hot drink and a thermometer (or just your hand) and note how long it takes to lose half its excess heat over the room. It will take about the same time again to lose half of what is left. That constant "half-life" is the time constant in disguise. Sources & Further Reading Sources. Euler's method: Institutionum calculi integralis, vol. 1 (1768); Encyclopedia of Mathematics, "Euler method". Newton's law of cooling originates in Newton's anonymous 1701 paper "Scala graduum caloris"; the modern differential form is not in the original (Wikipedia, "Newton's law of cooling", and history-of-science references). Fourier's Théorie analytique de la chaleur (1822). Millennium Prize problems announced by the Clay Mathematics Institute in 2000, Navier–Stokes among them. Computed in node: RC charging (τ = RC = 1 s for 10 kΩ and 100 µF, 63.2% at τ, 99.0% at 4.605τ); Euler at h = 0.5 gives 2.44 × 10−4 at t = 6 against the exact 2.48 × 10−3; spring closed-form solutions for the three regimes checked against a step-by-step integration (agreement to 5 decimals) and the first-swing overshoot e−πζ/√(1−ζ²) = 52.7% for ζ = 0.2; harvested logistic equilibria 50(1 ± √(1 − H/25)) = 27.64 and 72.36 for H = 20, with collapse after 4.3 years from a start of 25; the explicit heat scheme (100 cells, stable below 0.5) against the exact decay 0.8e−π²t = 0.2982 (numerical 0.2981) at t = 0.1. 🌐 MIT OpenCourseWare 18.03SC Differential Equations ▶️ 3Blue1Brown: differential equations 📘 Boyce, DiPrima and Meade, Elementary Differential Equations and Boundary Value Problems, 12th ed. (Wiley) 📘 Steven Strogatz, Nonlinear Dynamics and Chaos, 3rd ed. (CRC Press) 🌐 Newton's law of cooling (with the history of the law) Step 8. Statics: Why Things Stand Still --------------------------------------- Why: Nothing that stands still has unbalanced forces; that one rule designs bridges. Every bridge, shelf and crane is a bargain between forces that cancel exactly. Learn to see the bargain and you can read a structure at a glance. A shelf sags a little under its books, and then it stops. It does not creep lower all night, and it does not suddenly give up. Somewhere inside it, a quiet bargain has been struck: the wall pushes up on the bracket with exactly as much force as the books push down. If the books were heavier by one paperback, the bargain would be re-struck with a slightly bigger push and a slightly bigger sag. That bargain is the whole of statics. A bridge, a crane, a chair, a bone, the wing spar of an airliner on the ground: each is a collection of pushes and pulls that cancel exactly. When engineers say a structure is "in equilibrium," they are saying something stronger than "it isn't moving." They mean every piece of it, however small, has forces on it that add to zero. Once you can see that cancellation, a surprising amount of the built world becomes legible. You can look at a crane and tell which arm is being stretched and which squashed. You can see why a tall narrow thing falls over before it slides, and why a seesaw does not care how big the children are, only where they sit. The course in one box Who takes it: mechanical, aerospace and civil engineers, usually in the second year, straight after first-year physics. What it is for: working out the forces inside and on a thing that is standing still, so you know what each part has to carry. Comes after: Physics I: Mechanics (34/06) and Linear Algebra (34/04). Leads to: Mechanics of Materials (34/11), which asks how a part responds to those forces, and from there Machine Design (34/19). What this piece will and won't give you: you won't be able to size a real bridge member after this. You will see why a truss is made of triangles, and why a tall crate tips before a squat one. An easier first read on the same ground: Everyday Physics: The Answers Hiding in Plain Sight (Physics 5) and The Four Forces That Run Everything (Physics 1). Nothing is accelerating, so everything adds to zero Newton's second law says that a net force makes a body accelerate. Statics is the special case where nothing accelerates, so the net force is zero. That is one sentence, and an entire first course grows from it, because it can be applied to any piece of a structure you care to isolate. The habit that makes it work is the free-body diagram. Draw a closed outline around the thing you are studying (the whole bridge, one beam, one bolt) and cut every connection to the outside world. At each cut, draw the force that connection supplied: gravity, a push from a support, a pull from a cable. If you have drawn them all, the arrows must balance. If they don't, you forgot one. Most of a statics course is learning to draw that outline honestly. Two kinds of balance are needed for something that can turn as well as slide. The forces must add to zero, and the turning effects (moments, next section) must add to zero. In symbols: ΣF = 0    and    ΣM = 0 Read it as: add up every force, add up every turning effect, and both totals come out as nothing. In a flat problem that is three equations: horizontal forces, vertical forces, and moments. Three equations can find three unknowns, which is why a simple beam with one pin support and one roller support (three unknown reactions) can be solved on the back of an envelope. Supports are described by what they refuse to let the beam do. A roller blocks only up-and-down motion, so it supplies one force. A pin blocks sliding in both directions, so it supplies two. A fixed wall blocks sliding and turning, so it supplies two forces and a moment. The forces a support supplies are called reactions, and they are exactly as large as they need to be. Interactive · Where does the load go? A beam is held up at both ends. Drag the red load along it, or use the slider. Position of the load Size of the load Put 1,000 N on a 10 m beam 6 m from the left end. The right support carries 1,000 × 6 ÷ 10 = 600 N and the left carries 400 N. The nearer support always carries more. Slide the load to the middle and they split it equally; slide it over a support and that support takes everything. A force's turning power is force times distance Push a door near its hinge and it barely moves. Push at the handle and it swings. The force is the same; the distance from the hinge is not. The turning effect of a force about a point is the force multiplied by the shortest distance from the point to the force's line of action. Engineers call it a moment (physicists say torque). M = F × d    (newton-metres) This is why a long spanner loosens a stubborn nut: 100 N pushed 0.2 m from the nut gives 20 N·m, and the same push 0.4 m out gives 40 N·m. And it is why the second balance condition is so useful. If you add up moments about a point where an unknown force acts, that force has zero distance and drops out of the equation. Choosing the pivot cleverly is half the craft. A seesaw is the cleanest case. It balances when the turning effect on one side equals the turning effect on the other. The weights matter, and the distances matter, but only their product counts. Interactive · The seesaw Left mass Left distance from pivot Right mass Right distance from pivot Start with 30 kg on the left at 2 m and 20 kg on the right at 2 m. The left side wins. Now move the lighter child out to 3 m: 30 × 2 = 60 and 20 × 3 = 60, and the seesaw sits level. (The reading shows newton-metres, which is the same product times 9.81.) Gravity cancels from both sides, so the balance holds on the Moon too. A truss is a machine for sorting pushes from pulls Take a rectangle of four sticks pinned at the corners and it folds flat with a finger. Add one diagonal and it cannot. A triangle is the only polygon that is rigid with pinned corners, so bridges, roof frames, cranes and pylons are built of them. In an ideal truss every member is pinned at both ends and loads act only at the joints, so each member carries a force along its own length and nothing else: either a pull (tension) or a push (compression). To find which, use the method of joints. Isolate one joint, draw the arrows of the members meeting there, and make them balance. Start at a joint with only two unknown members, solve it, and move to the next. For a real truss with dozens of joints this becomes a big set of linear equations, which is exactly what Linear Algebra (34/04) is for. The tool below does it live: seven member forces and three support reactions, ten equations solved at once. Interactive · Which members pull, which push? A 6 m span truss. Red members are in tension, blue in compression; thicker means more force. Load at bottom middle Load at top left Load at top right Load Truss depth With 10 kN at the bottom middle and a depth of 3 m, each support carries 5 kN. The top member is squashed (compression, 5 kN), the two bottom members are stretched (2.5 kN of tension each), and each of the four sloping members carries 5.59 kN, alternating between push and pull. Now squash the truss to 1.5 m deep: the top member's force doubles to 10 kN. A shallow truss works its members much harder than a deep one, for the same reason a thin plank sags more than a deep beam. You will meet that reason again in Mechanics of Materials (34/11). Switch the load to a top joint and watch members change sex. A member that was pulling can start pushing. Engineers care, because tension and compression fail differently: a rod in tension can be thin, but a long member in compression may bow sideways long before it is crushed. Friction pushes back only as hard as it has to Friction is the most misread force in school physics. It is not a fixed number. If you push a heavy crate gently and it does not move, friction is pushing back with exactly your push, no more. Push harder and friction rises to match, until it hits a ceiling. Past the ceiling, the crate slides. For dry solid surfaces, the ceiling is roughly proportional to how hard the surfaces are pressed together (the normal force), and, to a good approximation, does not depend on the contact area. This is the Amontons–Coulomb law: Ffriction ≤ μ × N The coefficient μ (mu) depends on the pair of materials. A sliding object usually feels a slightly smaller friction than the ceiling that held it still, so there are two coefficients, static and kinetic. The numbers below are made up for illustration (0.5 static, 0.4 kinetic), not the values for any particular pair of materials. Interactive · Push the crate Mass of crate Your push A 50 kg crate presses on the floor with 490.5 N. At μ = 0.5, friction can resist up to 245 N. Push with 100 N and friction pushes back with exactly 100 N and nothing happens. Push with 300 N and it breaks free; kinetic friction is then 0.4 × 490.5 = 196 N, leaving 104 N to accelerate 50 kg at about 2.1 m/s². That sudden drop is also why things "stick-slip": a squeaking door or a violin bow. The centre of gravity decides what tips A body can be treated, for balance, as if all its weight acted at a single point: the centre of gravity. For a uniform object this is its centroid, the geometric centre of its shape. A rectangle's centroid is at the middle; a triangle's is one third of the way up from the base. For awkward shapes, engineers find the centroid by splitting them into simple pieces and averaging their positions, weighted by area. The rule for tipping is then simple and a little startling: an object stays up as long as the vertical line through its centre of gravity falls inside its base of support. The moment that line passes outside the edge it is resting on, gravity's turning effect about that edge has nothing to oppose it, and the object goes over. That gives a clean comparison between tipping and sliding on a tilting surface. A block tips when the tilt reaches the angle whose tangent is its width divided by its height. It slides when the tilt reaches the angle whose tangent is μ. Whichever threshold comes first wins. Interactive · Tilt the ramp Ramp angle Block height : width Friction coefficient μ A crate 1 m tall and 0.5 m wide tips when tan θ reaches 0.5, which is 26.6°. With μ = 0.4 it would slide at 21.8° first. Raise μ to 0.8 (rough rubber matting, say) and it slides only at 38.7°, so now it tips at 26.6° instead. Same crate, different outcome, and the reason is entirely in the comparison of two numbers. A low, wide object, or one with a heavy bottom, has a larger tipping angle. That is why lorries are loaded heavy-side-down and why a laden roof rack makes a car feel twitchy. What a graduate just sees Anything standing still has forces that add to zero, and turning effects that add to zero. If a piece of the picture does not balance, a force is missing. Before solving anything, they cut the thing out of the world and draw what each cut supplies. A support exerts exactly the reaction needed, and the nearer support carries more of a load. A long lever arm is a force multiplier; the pivot you pick for adding up moments can make an unknown vanish. Triangles are rigid. Every member of a truss is secretly either pulling or pushing along its length, and the pushed ones are the ones to worry about. A deep truss works its members less than a shallow one for the same load. Friction is a limit, not a number: it matches your push until the limit, then lets go. Tall and narrow tips, low and wide slides: compare the width-to-height ratio with μ. Where it shows up next Mechanics of Materials (34/11): statics gives the force in each member or cross-section; that course asks how big a bar must be so that force doesn't break or bend it. Dynamics (34/12): everything here is Newton's second law with the acceleration set to zero. Let it be non-zero and the same free-body diagram describes a moving machine. Machine Design (34/19): every shaft, gear tooth and bolt is sized from the reactions found by statics. Linear Algebra (34/04): real structures have hundreds of joints, and the equations are solved as matrices, just as the truss tool does. In the world: crane booms, bridge trusses, aircraft wing ribs, scaffolding, bike frames, the load rating on a ladder, the rollover rating of a vehicle. If you want to go deeper 📘 Standard textbook: R. C. Hibbeler, Engineering Mechanics: Statics (15th edition, Pearson). The most widely assigned text, with a great many worked problems. 📘 Standard textbook: J. L. Meriam and L. G. Kraige, Engineering Mechanics: Statics (Wiley). Another common choice, careful on free-body diagrams. 📘 Gentler: J. E. Gordon, Structures: Or Why Things Don't Fall Down (Penguin 1978; Da Capo reprint). No equations needed to follow it; written by a materials scientist who loves bridges, boats and bones. 🌐 Free book: Engineering Statics: Open and Interactive by Daniel W. Baker and William Haynes, a free online textbook with interactive figures. 🌐 Free course: MIT OpenCourseWare 2.001 Mechanics & Materials I, which treats statics as the first step into materials. The takeaway If a thing is standing still, every force on it, and every turning effect, adds to zero. Draw the thing alone, draw what each cut supplies, and the unknown pushes and pulls fall out of three equations. That is how a beam shares its load between supports, how a truss sorts its members into pushers and pullers, and why a tall crate tips while a squat one slides. Try this: this week, look at a bookshelf, a bridge or a bike stand and ask three questions. Where are the supports? Which parts are being pulled and which pushed? What would make it tip? Sources & Further Reading Sources. The equilibrium equations, support types, method of joints, Amontons–Coulomb friction, centroids and tipping are standard textbook statics (see the books below). Every number in the interactives was computed in node: beam reactions (RA = W(L−x)/L), seesaw moments, the ten-equation truss solve (checked by hand for the 3 m deep case: chords 5 kN and 2.5 kN, diagonals 5/sin 63.4° = 5.59 kN), friction (N = 50 × 9.81 = 490.5 N; a 300 N push gives 2.08 m/s²) and tipping (arctan 0.5 = 26.6°; arctan 0.4 = 21.8°; arctan 0.8 = 38.7°). The friction coefficients 0.5 and 0.4 are illustrative. The truss is the ideal, pin-jointed model, as in textbook problems; real joints are stiffer. Go deeper. Hibbeler's Statics, 15th edition, is listed by Pearson (ISBN 9780137514663). Gordon's Structures was first published by Penguin in 1978. 🌐 Engineering Statics: Open and Interactive 🌐 MIT OpenCourseWare 2.001 Mechanics & Materials I 📘 J. E. Gordon, Structures: Or Why Things Don't Fall Down Step 9. Mechanics of Materials: Why Things Bend, Buckle and Break ----------------------------------------------------------------- Why: Why beams bend, columns buckle and parts crack, and how to size them so they don't. A bridge cable, a bicycle frame and a paperclip are all springs with limits. Learn the few ideas that say how far they stretch, and what finally breaks them. Bend a paperclip once and it holds its new shape. Bend it back and forth a dozen times and it snaps in your fingers, even though you never once pulled as hard as it can be pulled. A ruler squeezed end to end skips the squashing and bows sideways with a sudden bang. A diving board flexes under your weight and springs back, but a bit more weight would leave it permanently bent. None of that is mysterious once you have the picture a mechanical engineer carries. Every solid is a very stiff spring. Every load makes it stretch or squash or twist by an amount you can predict, and the way it fails depends on what kind of load, how often, and where the stress piles up. Statics told you how big the forces are. This course is about what they do to the stuff. The course in one box Who takes it: mechanical, aerospace, civil and materials engineers, usually in the second year (sometimes called Strength of Materials or Solid Mechanics). What it is for: predicting how much a part deforms and when it will fail, so you can choose its shape and size. Comes after: Statics (34/10), and through it Physics I: Mechanics (34/06) and Linear Algebra (34/04). Leads to: Machine Design (34/19), which turns it into shafts, bolts and gears. It also underlies Materials Science & Manufacturing (34/13) and Vibrations (34/17). What this piece will and won't give you: you won't be able to certify a beam. You will see why a beam is shaped like an I, why a long thin column fails differently from a short fat one, and why a small hole can matter more than a big load. Easier first read: Everyday Physics: The Answers Hiding in Plain Sight (Physics 5). Read Statics: Why Things Stand Still (34/10) first. Stress and strain: divide out the size A thick rod and a thin rod made of the same steel and pulled with the same force behave very differently, and a long rod stretches more than a short one. To compare materials rather than particular objects, engineers divide the size out. Stress is force per unit area: how hard the material is working. Strain is stretch per unit length: how much it has deformed, as a fraction. A rod that grows from 1,000 mm to 1,001 mm has strain 0.001, or 0.1%. σ = F ÷ A     ε = ΔL ÷ L For small stretches, stress and strain are proportional. Robert Hooke published the idea in 1678 as ut tensio, sic vis, "as the extension, so the force". The constant of proportionality for a material is its Young's modulus E, named after Thomas Young. It is the material's stiffness, independent of shape: σ = E × ε Steel's E is about 200 GPa, aluminium's about 69 GPa, and pine's roughly 10 GPa (values vary by alloy and species). Read it as: for the same stress, aluminium stretches about three times as much as steel. Try it. Interactive · Pull a rod Steel Aluminium 6061-T6 Pulling force Rod diameter With 10 kN on a 10 mm steel rod 1 m long: area 78.5 mm², stress 127 MPa, strain 0.000637, stretch 0.64 mm. Halve the diameter and the stress quadruples, because the area falls with the square of the diameter. Swap to aluminium and the stretch nearly triples (200 ÷ 69) while the stress stays put. That is why stiffness-critical parts (not just strength-critical ones) are chosen by E. The stress–strain curve is a material's personality Pull a bar in a test machine, record stress against strain, and you get a curve that is a fingerprint of the material. It has a few landmarks worth knowing by heart. First comes a steep straight line: the elastic range. Let go and the bar returns exactly to its original length. The slope is E. At the yield point the line gives way: atoms begin to slip past each other, and release no longer restores the old shape, so some stretch is permanent. In a ductile metal like structural steel the curve then rises further (the metal work-hardens) to the ultimate strength, the highest stress it ever carries. After that the bar thins in one place (necking) and finally breaks. A brittle material such as glass or cast iron stays elastic almost to the end and snaps with little or no warning. Engineers prefer ductile materials for structures because ductility gives notice: a steel beam sags visibly and absorbs energy long before it breaks. The yield point, not the ultimate strength, is usually the number a designer works from. Interactive · Stretch it until it gives Illustrative curves: a steel-like metal (E = 200 GPa, yield 250 MPa, ultimate 400 MPa) and a glass-like brittle solid (E = 70 GPa, snaps at 50 MPa). They are idealised in shape, not measured. Ductile, steel-like Brittle, glass-like Whole curve Zoom in on the first 0.4% Stretch the bar to The elastic range of the steel-like metal ends at strain 0.00125, since 250 MPa ÷ 200,000 MPa = 0.00125. That is just 1.25 mm on a metre of bar, a sliver of the total 25% it can stretch before breaking. Use the zoom button to see that sliver. Nearly everything a structure does in service happens inside it, which is why it is designed to stay there. Beams: stiffness grows with depth cubed When a beam sags, its top edge is squashed and its bottom edge is stretched. Somewhere between them is a layer that does neither, the neutral axis. Material near the middle is barely working. Material at the top and bottom is working hardest and is also the farthest from the middle, which gives it the most leverage against bending. The shape's resistance to bending is captured in one number, the second moment of area I. For a rectangle of width b and depth d, I = b·d³ ÷ 12. The cube is the whole story. For a beam resting on two supports with a load in the middle: sag = P L³ ÷ (48 E I) Double the depth of a beam, keeping the width, and its stiffness rises 2³ = 8 times. Compare a plank on edge with the same plank laid flat: the edge-up plank is twice as deep (a factor of 8) but half as wide (a factor of ½), so it is 4 times stiffer. Doubling the span makes the sag 2³ = 8 times bigger. Depth and span, not weight, rule the stiffness. That also explains the I-beam. If material near the middle isn't doing much, move it to the flanges and leave a thin web to hold them apart. The shape gives up some stiffness for a lot less steel. Interactive · Shape the beam Steel Aluminium Pine Solid rectangle I-section Lay it flat Load at the middle A solid steel 50 × 100 mm beam on edge, 2 m span, 5 kN at the middle, sags 1.0 mm. Lay it flat and it sags 4 mm. Make the same outline an I-section (flanges and web each a tenth of the size) and it uses 28% of the steel but keeps 54% of the stiffness, sagging 1.85 mm. The beam is shaped like an I because that puts the steel where the bending lives. Buckling: slender columns bow before they crush You would expect a column to fail when the compression exceeds the material's strength. A short, stout one does. A long thin one never gets there. Press on the ends of a plastic ruler and at some load it suddenly bows out sideways, far below what it could carry if it stayed straight. That is buckling: an instability, not a strength failure. Leonhard Euler worked out the load at which a perfectly straight slender column bows, in the 1750s: Pcr = π² E I ÷ (K L)² The load depends on stiffness (E I), not strength, and falls with the square of the length: double the length and the load a column can take falls to a quarter. The factor K depends on how the ends are held: 1 for pinned ends, 0.5 if both are clamped (four times stronger), 2 for a flagpole fixed at the bottom only (a quarter as strong). A column's crushing load, on the other hand, is just yield stress times area and doesn't care how long it is. Interactive · Crush or bow? Pinned both ends Clamped both ends Clamped and pinned Flagpole Rod diameter Length A 20 mm steel rod, 1 m long and pinned at both ends, buckles at 15.5 kN, though it would take 78.5 kN to crush it. At 2 m it buckles at only 3.9 kN. At 0.2 m the buckling load is far above the crushing load, so it crushes instead. The formula is for ideal columns; real ones have imperfections and do worse, and real designs use empirical corrections with generous margins. In structures, compression members are the ones engineers watch for exactly this reason, and why a long pushed member is deeper or braced. Torsion: a twisted shaft is a bar in shear Stretching is not the only way to load a bar. Push the top of a deck of cards sideways and the cards slide over each other: that is shear, a sliding stress (force parallel to the surface, divided by the area). A bolt holding two plates that try to slide past each other is in shear. The shear modulus G plays the role for sliding that E plays for stretching; for steel it is about 80 GPa. Twist a shaft and every thin slice slides a little relative to its neighbour, and the sliding is biggest at the surface. The surface stress is torque times radius divided by the shaft's polar second moment J (for a solid round shaft J = πD⁴ ÷ 32): τ = T r ÷ J     twist = T L ÷ (G J) Because the centre of a shaft barely strains, you can bore it out and keep most of its torsional stiffness. That is why drive shafts, bicycle cranks and aircraft tubes are hollow. Interactive · Twist a shaft Outer diameter Hollowness Torque A solid 20 mm steel shaft carrying 100 N·m has a surface shear stress of 63.7 MPa and twists 4.6° over a metre. Bore it out to 16 mm (80% of the outer diameter) and it keeps 59% of the twisting stiffness with 36% of the steel. The surface stress rises to 108 MPa, so you do pay something, which is why designers also check it. The yield strength in shear is roughly 0.58 times the tensile yield (about 144 MPa for 250 MPa steel). Holes, cracks and repetition: why parts break below their strength Stress flows through a part like water through a pipe. Where the pipe narrows or bends, the flow bunches up. A hole in a plate in tension makes the stress lines crowd around it, so the stress at the hole's edge is higher than the average, by a stress concentration factor K. For a small round hole in a wide plate, K = 3. For an elliptical hole whose width across the pull is a times its height, K = 1 + 2a. A sharper, longer notch across the pull means a bigger K, and a crack is the sharpest notch of all. Interactive · A hole in a loaded plate Average stress in the plate Shape of the hole At 80 MPa average stress (well under yield) a round hole brings the edge to 240 MPa, a factor of safety of about 1.04. At 100 MPa it goes to 300 MPa and the metal there yields. A 5 : 1 slot puts K at 11. The factor of safety is just capacity divided by what you expect to load, with a deliberate margin for the unknowns. Aircraft designers typically work with slimmer margins than builders of bridges, because every kilogram of margin is paid for in fuel, so they have less room for a surprise like a concentration they missed. Repetition makes it worse. A load applied again and again, far below the yield stress, can start a tiny crack at a stress concentration and grow it a little each cycle until what remains cannot carry the load. That is fatigue, and it is the paperclip's story at engineering scale. The de Havilland Comet, an early jet airliner, suffered two mid-air break-ups over the Mediterranean in 1954 (10 January off Elba, 35 people killed; 8 April near Naples, 21 killed). At Farnborough, investigators sealed a whole fuselage in a water tank and pressurised it again and again to imitate flights. It burst open after 3,057 pressure cycles (1,221 real flights plus 1,836 simulated), from a crack at the corner of the forward escape-hatch window. On the Elba wreckage the fatigue crack began at a rivet hole by the opening for an aerial. The story is often told as "square windows"; sharp corners were stress concentrators, but they were one of several in that airframe. The lesson that stayed: repeated pressurisation plus a stress concentration can break even a strong skin, so pressurised airframes must be designed and tested for fatigue. What a graduate just sees Strength is stress (force over area), not force. Halving the diameter quadruples the stress. Every solid is a spring, and Young's modulus is its stiffness: steel stretches a third as much as aluminium for the same stress. Most of the working life of a structure is spent in the first fraction of a percent of strain, deliberately. Ductile means warning; brittle means none. A deep beam is stiff because stiffness grows with the cube of the depth. Put material where it is farthest from the middle. A long thin thing in compression is a buckling problem, not a crushing problem, and doubling its length cuts its capacity to a quarter. Hollow tubes are nearly as stiff in twist as solid bars at a fraction of the weight. Sharp corners, holes and cracks multiply stress; round them. Repetition finds them. Where it shows up next Statics (34/10) feeds it: the reaction or member force found there becomes the F and T here. Machine Design (34/19): sizing a shaft, a bolt or a gear tooth is these formulas with factors of safety and fatigue data added. Materials Science & Manufacturing (34/13): why steel yields where it does, why glass is brittle, and how heat treatment moves the curve you played with. Vibrations (34/17): a beam is a spring, and its stiffness sets the pitch at which it rings. In the world: aircraft skins and wing spars, bridge girders, bicycle frames, tall chimneys, engine crankshafts, the reason window corners on ships and planes are rounded. If you want to go deeper 📘 Standard textbook: R. C. Hibbeler, Mechanics of Materials (Pearson). Widely assigned, clear on method, with a great many worked problems. 📘 Standard textbook: F. P. Beer, E. R. Johnston Jr., J. T. DeWolf and D. F. Mazurek, Mechanics of Materials (McGraw Hill). The other common choice. 📘 Standard textbook: J. M. Gere and B. J. Goodno, Mechanics of Materials (Cengage). Known for patient explanations. 📘 Gentler: J. E. Gordon, Structures: Or Why Things Don't Fall Down, and its companion The New Science of Strong Materials. Written for non-engineers, and mostly equation-free. 🌐 Free course: MIT OpenCourseWare 2.001 Mechanics & Materials I, with lecture notes and problem sets. The takeaway Stress is force over area, strain is stretch over length, and Young's modulus is the exchange rate between them. A beam's stiffness rises with the cube of its depth, so it is shaped like an I. A slender column bows rather than crushes, at a load that falls with the square of its length. A twisted shaft is a bar in shear, and a hollow one is almost as good as a solid one. And every hole, corner and crack concentrates stress, which repeated loads then turn into fatigue. Try this: bend a paperclip back and forth and count the bends until it breaks. Then look at any bridge, crane or bicycle and find where the metal is placed far from the middle, and where the corners are rounded. Sources & Further Reading Sources. The formulas (σ = F/A, ε = ΔL/L, Hooke's law, I = bd³/12, centre-load sag PL³/48EI, Euler's Pcr = π²EI/(KL)², τ = Tr/J, twist TL/GJ, the Inglis ellipse K = 1 + 2a/b) are standard textbook results. Hooke published the anagram's solution ut tensio, sic vis in 1678 (Wikipedia, "Hooke's law"). Typical E values: Wikipedia, "Young's modulus" (steel about 200 GPa, aluminium about 68–69 GPa, wood roughly 9–12 GPa). Euler's column formula dates from 1757 (Wikipedia, "Buckling"). Comet: BOAC Flight 781, 10 January 1954, 35 dead; South African Airways Flight 201, 8 April 1954, 21 dead; water-tank test fuselage G-ALYU failed on 24 June 1954 after 3,057 cycles (1,221 actual flights plus 1,836 simulated), from Wikipedia, "De Havilland Comet"; the fatigue origin on G-ALYP at a rivet hole by the aerial (ADF) window, and the square-window caveat, from Aerossurance, "Common Comet Misconceptions". Every number in the interactives was computed in node (steel rod: 127.3 MPa, 0.637 mm; beam: 1.00, 4.00 and 1.85 mm; column: 15.5, 3.88 and 78.5 kN; shaft: 63.7 and 107.8 MPa, 4.56°; notch: K = 3, 11). The two stress–strain curves are idealised and the G, yield and ultimate values are round textbook figures. Go deeper. Gordon's two books were first published by Penguin in 1978 (Structures) and earlier (The New Science of Strong Materials). 🌐 MIT OpenCourseWare 2.001 Mechanics & Materials I 🌐 Aerossurance: Common Comet Misconceptions 📘 J. E. Gordon, Structures: Or Why Things Don't Fall Down Step 10. Materials Science & Manufacturing: Why Things Are Made of What They Are -------------------------------------------------------------------------------- Why: Concrete, steel, timber and asphalt: why each is used where it is. A paperclip bends, a mug shatters, a rubber band springs back. Same push, three different answers, and the whole course is the reason why. Take a steel paperclip and bend it open. It stays bent. Bend that same kink back and forth, and it gets stubborn at the bend, then warm, then it snaps, usually after a handful of bends. Now try the same with a ceramic mug handle: it doesn't bend at all, it just breaks. Try it with a rubber band: it stretches a long way and comes back every time. Three materials, three completely different answers to the same push. A materials scientist doesn't see three mysteries. They see three arrangements of atoms, and one question that explains the lot: what happens to the atoms when you load them? That question, plus a second one (how was this thing made, and what did making it do to the atoms?), is the whole course. It decides why a bicycle frame is aluminium, why a surgeon's scalpel is hard and a car spring is not, why an aeroplane wing is turning from metal to carbon fibre, and why two parts with the same shape and the same alloy can behave differently. The course in one box Who takes it: mechanical, aerospace, civil and materials engineers, usually in year 2, beside Mechanics of Materials. What it's for: choosing what a part is made of, and from that, how to make it. Comes after: Physics I: Mechanics (34/06). Leads to: Machine Design (34/19), which needs it to pick shaft steels and bolt grades; it also feeds Mechanics of Materials (34/11), Propulsion (34/21) and the Capstone (34/22). Honest line: you won't be able to pick a steel grade from a handbook after this; you will see why a blacksmith plunges a blade into water, and why a bridge is not made of glass. Atoms are springs: where stiffness comes from A solid is atoms held at a comfortable spacing by bonds. Push two atoms closer than that spacing and they repel strongly. Pull them apart and they attract. At the right distance, they sit in an energy valley. Near the bottom of any smooth valley, the shape is almost a parabola, which is exactly the shape of a spring. So, for small stretches, every bond is a spring, and a solid is a huge three-dimensional mattress of them. Engineers measure how stiff the mattress is by stress (force per area) against strain (fractional stretch). Their ratio is Young's modulus, E: stress = E × strain    (σ = E ε) In words: the stiffer the bonds and the tighter the packing, the more stress it takes to stretch the material by the same fraction. Diamond sits at roughly 1,000 GPa. Steel is about 200 GPa, aluminium about 70, nylon about 3, and rubber 0.01 to 0.1, because rubber's long molecular chains are coiled like tangled string, so stretching it mostly uncoils the chains rather than stretching bonds. Worked number: pull a steel bar 1 m long with a 10 mm × 10 mm section using 10 kN. Stress = 10,000 N ÷ 100 mm² = 100 MPa. Strain = 100 MPa ÷ 200,000 MPa = 0.0005, so the bar grows 0.5 mm. The atomic spacing, about 0.25 nm, has changed by about 0.0001 nm. Now pull with a gentler 1 kN on the same section: steel stretches 0.05 mm, aluminium 0.14 mm, nylon 3.4 mm. Stiffness is a property of the bonds. Heat-treating a steel makes it far stronger, and does almost nothing to its stiffness. That last point is the surprise that a first-year loses and a graduate keeps: a soft mild steel and a hardened spring steel have almost the same E. A thin spring steel wire does not sag less than a mild steel one. It simply springs back from a much bigger load. Four words engineers keep apart: stiff, strong, tough, hard In everyday speech they blur together. In engineering they are different questions, and a material can say yes to one and no to another. Stiff: how little it stretches elastically under load (E). Strong: how much stress it carries before it yields (starts to stretch permanently) or breaks. Tough: how much energy it absorbs before it breaks, which is the area under its stress-strain curve. Hard: how well its surface resists being dented or scratched. Glass is stiff and, in a clean sheet, decently strong, but not tough: it breaks with no warning because cracks run through it. A rubber band is not stiff or strong, but is tough. Mild steel is the compromise that built the modern world: stiff, strong enough, and above all tough, because it yields, warns you by bending, and soaks up energy. Try the pull test below. Interactive · Pull test Mild steel Hardened spring steel Aluminium 6061-T6 Nylon 66 Glass Pull until it breaks (left to right) Simplified curves built from typical handbook values (modulus, yield, tensile strength, elongation); real batches vary. "Energy to fracture" is the area under the whole curve. The hardened spring steel row is illustrative. Each plot is scaled to its own material, so compare the table, not the picture. Compare steel and glass in the table. Their stiffness is within a factor of three, and yet steel absorbs roughly 4,500 times more energy per volume before it fails. Hardened spring steel has the same stiffness again, over five times the yield stress, and about the same area, but loses the long warning ductility gives. Strength, stiffness and toughness really are three separate dials. (Cracks have their own measure, fracture toughness; the area under the curve is the version for an uncracked bar.) Why metals are soft: crystals, grains and dislocations Metals are crystals: atoms stacked in neat repeating layers. A simple calculation of how much shear it takes to slide one whole layer over the next gives a stress of roughly one-tenth of the shear modulus G. For steel, G is about 80 GPa, so theory promises something like 8 GPa. Real mild steel yields at about 0.25 GPa, some thirty times lower, and for many pure metals the gap is hundreds of times. Something is wrong with the picture. In 1934, three researchers (Taylor, Orowan and Polanyi, working separately) gave the answer: a real crystal is not perfect. It contains dislocations, lines where an extra half-plane of atoms is wedged in. Slip doesn't need a whole layer to slide at once. It needs the dislocation to creep along, shifting a few atoms at a time, the way a ruck in a carpet moves far more easily than dragging the whole carpet. So everything a metallurgist does to make a metal stronger is a way of getting in the way of dislocations: Smaller grains. A metal part is made of many small crystals (grains) with misaligned layers. A dislocation can't cross the boundary easily, so more boundaries mean more strength. Alloying. Foreign atoms snag dislocations. Work hardening. Deforming a metal makes more dislocations that tangle each other. That's your paperclip getting stubborn at the bend. Precipitates. Tiny hard particles in the way (how aircraft aluminium is strengthened). The grain-size effect has a famous simple rule, the Hall–Petch relation: yield strength = σ0 + k ÷ √d Yield strength is a base strength plus a term that grows as the grain diameter d shrinks, with the square root making it gentler than you'd guess. For one measured low-carbon steel, σ₀ ≈ 100 MPa and k ≈ 600 MPa·µm½. Then grains of 25 µm give about 220 MPa; 10 µm give about 290; 3 µm give about 450. Interactive · Shrink the grains Average grain diameter The grains are a random pattern at the chosen average size, not a micrograph. Strength uses the Hall–Petch fit for one low-carbon steel; other metals have other constants, and the rule fails at nanometre grain sizes. Heat treatment: cooking steel Steel is iron with a little carbon (up to about 2%). Heat it red-hot and its atoms rearrange into a form that dissolves the carbon. Cool it slowly and the carbon has time to leave, leaving soft, tough mixtures of tiny crystals. Quench it (plunge it into water or oil) and the carbon has no time to escape. It is trapped, and the lattice is forced into a strained, distorted shape called martensite. That shape is extremely hard, and dislocations barely move. It is also brittle, and often has internal stresses big enough to crack the part. The fix is tempering: reheat the quenched steel to somewhere between about 150 and 600 °C. The trapped carbon can move a little, the strain relaxes, and you trade hardness for toughness. The higher the temper, the softer and tougher. A blacksmith reads the temperature from the colour of the thin oxide film on the bright steel: faint yellow, then straw, brown, purple, blue. Choose a tool below and see the colour it is tempered to. Interactive · Temper it Razor Rock drill Cold chisel Screwdriver Spring Tempering temperature Colours and uses are the traditional blacksmith's chart for carbon steel; the exact temperatures depend on the steel and how it was heated. The hardness and toughness trend is shown in words, because real numbers depend on the alloy. Tools that need an edge, such as razors and drills, are tempered low: hard, a little brittle. Things that must flex or take shocks are tempered higher. Notice what is not changing: the composition. Same iron, same carbon, different heat history, a different part. The same logic, with a different mechanism, is how aircraft aluminium is "aged" to full strength. Choosing a material: the Ashby trade-off Engineers rarely want "the strongest" or "the stiffest". They want the lightest part that is stiff enough, or the cheapest that is strong enough. The materials scientist Michael Ashby turned this into a method: plot every material on a chart of two properties (log scales, because properties span many powers of ten), and use a performance index to rank them for a given job. For a part that must stay below a given deflection at the lowest weight, the index depends on the shape of the job. For a tie rod or cable in tension it is E ÷ ρ (stiffness over density). For a beam whose cross-section you may scale, it is √E ÷ ρ. For a flat panel, ∛E ÷ ρ. Why the roots? Because in a beam you can use a stiffer material or make it deeper, and depth is worth a lot (it counts by the cube). Light, bulky materials win. Interactive · Which material makes the lightest part? Tie rod Beam Panel Anchor the guide line on: Mass at equal stiffness is 100% for steel; smaller is lighter. Data are typical values (carbon fibre depends strongly on layup; the "all fibres aligned" figure applies only along the fibres). Stiffness-limited design only: strength, toughness, cost and fatigue change the ranking. Two things stand out. First, for a tie rod, steel, aluminium, titanium and magnesium all score about the same (the metals are stiff in proportion to their weight). Aluminium wins only when the part can be made bulkier, in a beam or a panel, where a beam of aluminium is about 41% lighter than a steel one of equal stiffness. That is the real reason aeroplanes are aluminium: skin and spars are thin-walled, bulky shapes. Second, carbon fibre composite beats everything metallic, but only along its fibres; that's why a composite wing is built by laying plies in chosen directions. Boeing has said that about half of the 787's primary structure by weight is composite. Corrosion: a battery you didn't mean to build Rust is not "iron getting old". It is an electrochemical cell. Where iron meets water and oxygen, some patches give up electrons and dissolve (the anode), and the electrons run through the metal to other patches where oxygen picks them up (the cathode). Any water film is the electrolyte. The same thing happens on a larger scale when two different metals touch while wet. The more active one becomes the anode and corrodes faster than it would alone; the nobler one is protected. Galvanising uses this deliberately: zinc sacrifices itself to protect steel, and ships carry blocks of zinc or aluminium for the same job. Aluminium and stainless steel resist corrosion because of a thin, tough, self-healing oxide film (passivation) that blocks the reaction. Area matters too. The cathode reaction needs surface to run on, so a small anode with a big cathode concentrates the attack on a small area. Steel rivets in a big copper sheet is the bad case; a copper rivet in a steel plate, the forgiving one. Interactive · Who corrodes? Metal one Metal two Surface area, metal one : metal two Order follows the standard seawater galvanic series (most active first). It is a guide, not a prediction: real rates depend on the electrolyte, temperature and surface films. The area rule is a first-order approximation. How you make it changes what you made A drawing says "steel, this shape". It doesn't say how the steel was formed, and that changes the part: Casting (pour liquid metal into a mould) makes complex shapes cheaply, but the metal can trap gas pores and shrinks as it freezes, so cast parts tend to be weaker than wrought ones. Forging (hammer or press hot metal into shape) squeezes out pores and makes the grains follow the shape, like grain in wood. Crankshafts and connecting rods are often forged. Machining (cutting from solid stock) is accurate and keeps the stock's properties, but wastes material and leaves tool marks that can start fatigue cracks. Additive manufacturing (3D printing, building in layers) needs no tooling and can make internal channels no other process can, but parts are often weaker across layers than along them, and surfaces are rough. Injection moulding (molten polymer injected into a steel mould) makes millions of identical parts cheaply, after paying for the mould. Plastic chains line up in the flow direction, and where two flows meet they weld weakly. Cost shapes the choice as much as physics. Each process has a one-off cost (tooling, setup) and a cost per part. Divide the first over the quantity and add the second, and each process wins in its own range. Interactive · Which process, at what quantity? How many parts? The dollar figures are invented round numbers to show the shape of the trade-off, not quotes. Real tooling costs vary hugely (a simple mould can cost thousands, a car-door die millions), and the processes make parts in different materials with different properties. What a graduate just sees Stiffness belongs to the bonds, strength belongs to the microstructure: heat treatment changes the second and barely touches the first. Stiff, strong, tough and hard are four different questions, and brittle failure comes with no warning. Metals are weaker than their atoms promise because of dislocations, and every strengthening trick is a way to block them. A bent paperclip is harder at the bend: work hardening, then fatigue. Same alloy, different heat treatment or process, different part. "Lightest for a given stiffness" has an answer that depends on the shape of the job, and bulky, light materials win in bending. Two metals touching in water make a battery; the active one pays. Tooling cost against per-part cost decides how things are made. Where it shows up next Mechanics of Materials (34/11) uses E and yield strength as the inputs for sizing beams and shafts. Machine Design (34/19) is where you choose the steel for a gear or shaft and the grade of a bolt, using the heat-treatment and strength trade-offs from here. Propulsion (34/21): turbine blades live in gas hotter than their own melting point, using single-crystal alloys and cooling channels. Materials set the engine's temperature limit. Engineering Drawing & CAD (34/18): the drawing's material note and surface finish are how this course reaches the shop floor. In the real world: an aircraft's choice of aluminium or carbon, a car's crumple zone, a hip implant (titanium), the case of your phone. If you want to go deeper 📘 Standard texts: William Callister and David Rethwisch, Materials Science and Engineering: An Introduction (Wiley, 10th edition 2020), the most widely assigned survey. Michael Ashby, Materials Selection in Mechanical Design (Butterworth-Heinemann), for the property-chart method used above. For manufacturing, Kalpakjian and Schmid, Manufacturing Engineering and Technology (Pearson). 📘 Gentler: Mark Miodownik, Stuff Matters (2013), a prize-winning tour of everyday materials. J. E. Gordon, The New Science of Strong Materials, a classic with no equations required. 🌐 Free: MIT OpenCourseWare 3.091, Introduction to Solid-State Chemistry (bonding and crystals), and the University of Cambridge's DoITPoMS teaching packages, with micrographs and simulations. The takeaway Atoms are springs, so stiffness comes from the bonds and hardly changes with processing. Strength, hardness and toughness come from the microstructure (grains, dislocations, precipitates, trapped carbon), which heat treatment and manufacturing control. To choose a material, ask which property limits the part and use the right index, because the best material for a rod, a beam and a panel differ. And wherever metals meet water, expect a battery. Try this: bend a paperclip once, then bend the kink back and forth until it snaps, counting the bends. Feel the kink warm up. Then find two different metals touching in your home (a steel screw in a brass or aluminium fitting) and look for white powder or rust at the join. Sources & Further Reading Sources. Young's moduli: Wikipedia, "Young's modulus" (approximate values for relative comparison); aluminium and nylon properties from typical data sheets (MakeItFrom). Ideal shear strength about G/10 and the 1934 dislocation papers of Taylor, Orowan and Polanyi: standard textbook account (e.g. Wikipedia, "Dislocation"). Hall–Petch constants for low-carbon steel (σ₀ ≈ 100 MPa, k ≈ 600 MPa·µm½): Gotoh et al., Kanazawa University repository record on ECAP low-carbon steel. Tempering colours and uses: Wikipedia, "Tempering (metallurgy)". Performance indices √E/ρ, E/ρ, ∛E/ρ: Ashby, Materials Selection in Mechanical Design, summarised in NPTEL and Ansys teaching notes. Galvanic series: Wikipedia, "Galvanic series". Boeing 787 composite share: SAMPE/NASAMPE summary of Boeing figures. Numbers computed in node: bar stretch (0.5 mm), comparison stretches, energies to fracture, mass ratios, Hall–Petch values, and process switch points. 🌐 MIT OpenCourseWare: 3.091 Introduction to Solid-State Chemistry 🌐 DoITPoMS, University of Cambridge 📘 Mark Miodownik, Stuff Matters (2013) 📘 Callister and Rethwisch, Materials Science and Engineering: An Introduction Step 11. Probability & Statistics: Reasoning with Uncertainty ------------------------------------------------------------- Why: Loads, floods and material strengths are random; codes are built on statistics. Nothing you measure is exact, and no part comes out the same twice. Statistics is the discipline of knowing how much to trust a number. A machinist turns a hundred steel shafts to "10 millimetres". Measure them with a good micrometer and no two are the same: 9.981, 10.012, 10.004. Nothing is wrong. That scatter is a fact of every process, every material and every sensor. Now stack five of those parts end to end in a gearbox. Will the stack still fit the housing? A pollster calls a thousand people and says 52% will vote yes, give or take three. A drug trial says "significant". A forecast says 30% chance of rain. Every one of these is the same question: given scatter, how much should I trust this number? A graduate of probability and statistics doesn't see a single number any more. They see the shape of the scatter behind it. This piece builds that sight with eight things you can run: dice, tolerances, polls, a drug-trial trap, a belief that updates, and a gearbox stack. The course in one box Who takes it: almost every engineering and science student, usually in the second year; computer scientists meet it again before machine learning. What it's for: deciding what to believe, and how far to trust a measurement, when outcomes are uncertain. Comes after: Calculus I (34/01), which helps with continuous distributions, though the ideas here need little more than arithmetic. Leads to: Machine Learning (34/32) needs it directly. It also shows up in Materials Science & Manufacturing (34/13: process control and defects), Systems Engineering (34/22: reliability) and Computer Security (34/33: randomness and guessing). Honest line: you won't be able to design a clinical trial after this. You will see why "significant" doesn't mean "important", why a bigger sample beats a cleverer guess, and why most real quantities look like a bell curve. For easier first reads, see the library's Probability: How to Not Fool Yourself (12/03), Statistics for Reading the World (12/05) and Bayes' Theorem (12/17); this piece goes further into the engineering side. Random variables: a shape with a centre and a spread A random variable is a number that comes out of a process you can't predict one run at a time. The key move is to stop asking "what will it be?" and ask "what is the shape of what it could be?" Roll two dice and add them: 7 is six times as likely as 2 or 12, so the shape is a tent. Roll many times and the tent appears. Interactive · Roll dice, grow a shape Dice added together Roll 10 timesRoll 1,000 timesClear Without JavaScript: two dice add to 7 with probability 6/36 and to 2 or 12 with probability 1/36 each. The mean total is 7 and the standard deviation is about 2.42. Two numbers summarise a shape. The mean is its balance point: 7 for two dice. The standard deviation (the "spread") says how far values typically stray from the mean: about 2.4 for two dice, which means totals usually land between roughly 5 and 9. The variance is the standard deviation squared, and it is the one that adds up neatly when you combine independent things. standard deviation = √( average of (value − mean)² ) Now the shape that matters most. Measured quantities (shaft diameters, heights, sensor noise) tend to form the normal curve, the bell. Its great convenience is a fixed rule: about 68.3% of values fall within one standard deviation of the mean, 95.4% within two, and 99.7% within three. Engineers use it every day. Say the shaft's spec is 10.00 mm ± 0.05 mm. Slide the process around and see what share of parts get thrown away. Interactive · How many parts fail the spec? Average diameter Spread (standard deviation) Without JavaScript: a centred process with standard deviation 0.025 mm (the spec limits are exactly ±2 standard deviations) scraps about 4.55% of parts. Cut the spread to 0.01 mm and about 0.00006% are scrapped. Shift the average 0.02 mm off-centre with spread 0.02 mm and about 6.7% are scrapped. Two lessons graduates keep. First, spread, not average, decides quality: a process centred perfectly but wide still scraps parts. Second, being off-centre costs a lot even when the spread is fine. That is why factories track both. The central limit theorem: why the bell keeps appearing Why should so many things be bell-shaped? The answer is one of the most useful theorems in mathematics. Take almost any random process, however lopsided, and average a few dozen independent draws. The averages pile up in a bell, centred on the true mean, with spread shrinking like 1/√n. Early forms are due to de Moivre in 1733 (for coin flips) and Laplace, whose general version dates from about 1810. Pick a lopsided source and drag the sample size. Interactive · Average anything, get a bell Flat (0 to 1) Very lopsided Coin (0 or 1) Draws averaged together: n Without JavaScript: for the lopsided source (mean 1, standard deviation 1) the average of n = 30 draws has standard deviation 1/√30 = 0.183, and the 3,000 averages form a close-to-normal bell around 1. At n = 1 you see the raw source. By n = 30 even the lopsided one is a near-bell. The prediction (the curve) uses only two numbers from the source, its mean and its standard deviation. The shape of the source is forgotten. That is why a measurement made of many small independent errors is normal, and why so much of statistics can pretend things are bell-shaped. standard error of an average = σ / √n In words: the average of n draws wobbles less than one draw, by the square root of n. To halve the wobble you need four times the data, not twice. It is the quietest expensive fact in all of experimental science. Sampling and confidence intervals: how sure is a poll? You can't ask everyone, so you ask n people and treat the share saying yes as an estimate. How far off might it be? The standard error from the last section gives the answer. For a share p the wobble is √(p(1−p)/n), and a 95% confidence interval is the estimate plus and minus about two of those. For n = 1,000 and a share near half, that is ±3.1 percentage points: the familiar "margin of error". The wording is subtle, and the interactive makes it concrete. "95% confidence" is a statement about the method, not about one interval: if you ran the procedure again and again, about 95 intervals in 100 would capture the true value. Here the truth is secretly 52%. Run 100 polls. Interactive · 100 polls, one truth People asked in each poll Run 100 polls Without JavaScript: with 1,000 people per poll the margin of error is about 3.1 points; run 100 polls and about 95 of the intervals cover the true value. With 4,000 people the margin halves to about 1.5 points. Watch two things. Roughly five in a hundred intervals miss the truth even when everything is done correctly. And quadrupling the sample from 1,000 to 4,000 only halves the interval. A poll of a country and a poll of a town need the same sample size for the same precision, as long as the sample is random. The real danger is not the margin of error; it is a sample that isn't random at all, which no formula fixes. P-values and Bayes: two ways to weigh evidence A hypothesis test asks whether data are surprising if nothing is going on. The p-value is the probability, assuming nothing is going on, of seeing a result at least as extreme as yours. Flip a coin 100 times and get 60 heads: for a fair coin the chance of a split at least that lopsided in either direction is about 5.7%. It is not the probability that the coin is fair, and it is not the probability that the result was "due to chance". The American Statistical Association issued a rare formal statement on this in 2016, with six principles. Two of them: p-values do not measure the probability that a hypothesis is true, and a p-value does not measure the size or importance of an effect. The classic trap is testing many things. At the usual cut-off of 0.05, a test of a true "nothing" comes out "significant" one time in 20 by luck. Test twenty unrelated things and the chance that at least one lights up is 1 − 0.9520, about 64%. Try it: each press tests twenty effects that are, by construction, nothing. Interactive · Test twenty things that do nothing Run 20 testsReset the tally Without JavaScript: twenty tests of pure noise at the 0.05 level give at least one false "discovery" about 64% of the time, and one in 20 tests on average. If you report only the one that lit up, you have a headline and no finding. This is part of why many published results fail to repeat: in the Open Science Collaboration's 2015 attempt to replicate 100 psychology studies, 97% of the originals had statistically significant results but only 36% of the replications did. Bayes' rule answers a different question: given the data, how should my belief move? Start with a prior belief, multiply by how well each possibility explains the data, and renormalise. Say you test parts and want the defect rate. Know nothing at first (any rate from 0% to 100% equally likely) and add parts one at a time. Interactive · Learn a defect rate from parts Add a good partAdd 10 good partsAdd a defective partReset Without JavaScript: after 20 good parts and no defects, a flat starting belief gives a most-plausible defect rate near 0%, an average of 4.5%, and a 90% credible range of roughly a quarter of a percent to 13%. Twenty clean parts cannot rule out a defect rate near one in eight. Frequentists and Bayesians disagree about the philosophy: one treats probability as long-run frequency and never assigns it to a hypothesis, the other as degrees of belief that data update. In practice both are used, often on the same project, and on big data sets they usually agree. The Bayesian version has one practical advantage here: after twenty clean parts it honestly says "we still can't rule out one in eight". Correlation, regression and the trap of causation Correlation (called r) measures how tightly two quantities move together along a straight line, from −1 to +1. Regression fits the best line through the cloud (by least squares, as in Linear Algebra, 34/04). Both are summaries, and any summary can hide the story. In 1973 the statistician Francis Anscombe made four data sets with, to two decimals, the same mean, the same spread, the same r (0.82) and the same line (y = 3 + 0.5x). Look at each. Interactive · Four data sets, one set of numbers Set 1Set 2Set 3Set 4 Without JavaScript: all four Anscombe sets have mean x = 9, mean y = 7.50, correlation 0.816 and line y = 3.00 + 0.50x. Set 1 is a noisy line, set 2 a curve, set 3 a line with one outlier, and set 4 a single far point that creates the whole slope. The numbers agree and the stories do not: one noisy line, one curve, one line with a rogue point, one vertical stack held up by a single point. Always plot first. The other trap is the famous one. Ice-cream sales and drownings rise together, because summer heat drives both. A confounder (a hidden third factor) makes two things move together without either causing the other. The remedy graduates know is randomised assignment: if a coin decides who gets the treatment, nothing else can be systematically linked to it. Where you can't randomise, causal claims need much more than a line through a cloud. Engineering: tolerance stacks, reliability and Monte Carlo Five parts, each 10 mm ± 0.1 mm, line up in an assembly. What is the worst the total can be out? Worst-case stacking adds the tolerances: ±0.5 mm. It is safe and expensive, because it assumes every part is at its extreme in the same direction at once. If the parts vary independently, variances add, so the standard deviations add in quadrature: root-sum-square, √5 × 0.1 = ±0.224 mm. For five parts that is under half the worst-case figure. When the formulas get hard, engineers use Monte Carlo: let the computer build the assembly ten thousand times with random parts and count how often it fails. Here each part's tolerance is treated as three standard deviations. Choose how many parts are stacked. Interactive · Monte Carlo tolerance stack Parts in the stack Build 10,000 assemblies again Without JavaScript: five parts at ±0.1 mm (three standard deviations) give a worst-case limit of ±0.5 mm, a root-sum-square limit of ±0.224 mm, and in simulation only about 0.3% of assemblies exceed the root-sum-square limit; none come close to 0.5 mm. The honest caveat is that root-sum-square assumes the parts vary independently and are centred. A tool wearing in the same direction on every part breaks that assumption, and so does one supplier's batch. Engineers disagree about when to trust statistical stacking over worst-case; safety-critical designs often insist on worst-case. Reliability follows the same arithmetic. A system that needs all of its parts is a series system, and probabilities multiply: ten components that each work 99% of the time give 0.9910 ≈ 90.4% for the whole. Add a backup in parallel and failure needs both to fail: two units that each work 90% of the time give 1 − 0.1² = 99%. The two rules drive why aircraft duplicate systems and why complicated machines need very reliable components. What a graduate just sees Every number has a spread behind it; the first question is "how scattered?", not "what is it?". Spread, not the average, decides how many parts are scrapped, and centring matters too. Averages of many independent things are bell-shaped, whatever the source looked like. Halving the uncertainty costs four times the data. A 95% interval is a claim about a method: it misses one time in twenty by design. A p-value is "how surprising is this if nothing is going on", not "how likely is it true". Test twenty things and one will look significant; ask how many were tried. Plot the data; correlation is not causation, and only randomising breaks a confounder. Independent errors add in quadrature; series reliabilities multiply. Where it shows up next Machine Learning (34/32): training is fitting a model to noisy data; ideas of overfitting, error bars and Bayes are this course made bigger. Materials Science & Manufacturing (34/13): process capability and scrap rates are the normal curve against a spec limit. Mechanics of Materials (34/11) and Machine Design (34/19): safety factors exist because loads and strengths are distributions, not numbers. Systems Engineering (34/22): reliability budgets multiply series probabilities and add redundancy to push them back up. Computer Security (34/33) and Algorithms (34/24): good random numbers, randomised algorithms, and the odds of guessing a secret. If you want to go deeper 📘 Douglas Montgomery and George Runger, Applied Statistics and Probability for Engineers (Wiley; 7th edition). The standard engineering course text, with plenty of tolerance, quality and reliability problems. 📘 Joseph Blitzstein and Jessica Hwang, Introduction to Probability (CRC Press; 2nd edition, 2019). Best for deep intuition about probability itself, from Harvard's Stat 110. 📖 David Spiegelhalter, The Art of Statistics (2019): the most readable route for learning alone, built around real cases and honest about uncertainty. 🌐 MIT OpenCourseWare 18.05, Introduction to Probability and Statistics (Spring 2022): free notes and problem sets, with Bayesian and classical statistics side by side. The takeaway Nothing you measure is exact, so every number is a shape. Summarise a shape by its centre and its spread, trust averages more as the square root of the sample, read a confidence interval as a promise about a method, treat a p-value as a measure of surprise, and plot before you believe a correlation. Try this: next time you see a poll or a "significant" headline, find the sample size and ask how many things were tested. Then, in the dice widget above, add four dice and watch the tent turn into a bell. Sources & Further Reading Sources. The normal-curve rule (68.3%, 95.4%, 99.7%), the central limit theorem, confidence intervals, the binomial test, the Beta posterior, root-sum-square stacking and series/parallel reliability are standard textbook results (Montgomery and Runger; MIT 18.05); every number shown was recomputed in code. Wasserstein and Lazar, "The ASA's Statement on p-Values: Context, Process, and Purpose", The American Statistician (2016). Open Science Collaboration, "Estimating the reproducibility of psychological science", Science (2015): 97% of originals significant against 36% of replications. Anscombe, "Graphs in Statistical Analysis", The American Statistician (1973). Early forms of the central limit theorem: de Moivre (1733) and Laplace (c. 1810), as surveyed in recent histories such as Ulyanov, "From classical to modern central limit theorems" (arXiv 2405.19828). Go deeper. The books and course above; Wasserstein and Lazar's statement is free and short (doi.org/10.1080/00031305.2016.1154108). Step 12. Fluid Mechanics: Why Air and Water Push Back ----------------------------------------------------- Why: Water pushes, flows and erodes; fluid mechanics is the rulebook. Water presses, air lifts, and the equations that govern both are worth a million dollars to anyone who can prove they behave. What a graduate still sees in every flow. Stand at the bottom of a swimming pool's deep end and your ears tell you something about the water above. Watch a heavy airliner climb off a runway. Tilt your hand out of a moving car's window and feel it rise. Squeeze a garden hose and watch the jet go further. All of these are one subject: a fluid, which is just something that flows, pushing on things. The strange part is that the full rulebook, the Navier–Stokes equations, has been known for about 180 years, fits on one line, and nobody can prove that it always behaves. Engineers use it anyway, every day. A fluid mechanics graduate learns which simple pictures work (and when), and which ones are traps. The course in one box Who takes it: mechanical, aerospace, civil and chemical engineers, usually in year 3. What it's for: predicting the forces, flows and losses when air or water moves past, through or around something. Comes after: Multivariable Calculus (34/02), Differential Equations (34/03), Thermodynamics (34/14). Leads to: Heat Transfer (34/16), Propulsion (34/21), and everything that flies or pumps. It also meets Vibrations (34/17) in wind-driven failures and Control Systems (34/20) in flight control. Honest line: you won't be able to size a pump or design a wing after this. You will see why a dam is thick at the bottom, why a wing needs speed, and why the Navier–Stokes problem is worth a million dollars. Pressure: the weight of everything above you Pressure is force per area. In a still fluid it comes from one thing: the weight of the column of fluid above. Dive down and you carry more water, so pressure rises in a straight line with depth. It pushes equally in every direction, which is why a diver's ears feel it from both sides. p = ρ · g · hpressure (extra, above the air) = density × gravity × depth For water that is about 9.8 kilopascals per metre, so 10 m down you carry about one extra atmosphere. Now the dam. The pressure grows with depth, so the push on the wall is a triangle, tiny at the top and largest at the bottom. Add up the triangle and the total push per metre of wall is ½·ρ·g·H². Double the water depth and the push quadruples. That is why a dam is thick at the bottom, and why a reservoir twice as deep is a very different structure, not a slightly bigger one. Interactive · The dam Water depth behind the wall At 10 m deep: base pressure 98 kPa, total push 491 kN per metre of wall. Fresh water, 1,000 kg/m³, g = 9.81 m/s². Try 10 m, then 20 m: the push goes from 491 to 1,962 kN per metre. Four times, not two. Buoyancy: the weight of the water you pushed aside The pressure under an object is bigger than the pressure on top, because it is deeper. That difference is a net upward push. Archimedes' principle says exactly how big: the weight of the fluid the object displaced. A block floats when it can displace its own weight before it is fully under. So it floats low or high depending on one ratio, its density against the fluid's. A steel ship is the proof that shape matters: steel is eight times denser than water, but a hull encloses so much air that the average density of the whole ship is below water's. An iceberg is ice (about 917 kg/m³) in seawater (about 1,025): 917 ÷ 1,025 = 0.895, so about 89% of it is below the surface. Interactive · Float or sink Fresh waterSeawaterMercury Density of the block Set to ice (917) A 600 kg/m³ block floats with 60% below the surface. Rough values: fresh water 1,000, seawater 1,025, mercury 13,534 kg/m³. Try a 7,800 block (steel) in mercury: it floats. Continuity and Bernoulli: faster flow, lower pressure Mass can't vanish. If water flows through a pipe that narrows, the same amount must pass each cross-section every second, so it must speed up. That is the continuity equation: area times speed stays constant. Squeeze the hose end with your thumb and the jet speeds up by exactly the area ratio. A₁v₁ = A₂v₂area × speed = the same all along the pipe Now the clever step. Speeding water up takes energy, and in smooth flow that energy comes from pressure. Daniel Bernoulli's relation is just energy conservation along a streamline: pressure plus the kinetic energy per volume plus the height term stays constant. p + ½ρv² + ρgh = constantpressure + motion energy + height energy, per cubic metre So where the water speeds up, the pressure drops. Worked example: a 10 cm pipe narrows to 5 cm. The area ratio is 4, so water at 2 m/s becomes 8 m/s, and the pressure falls by ½·1,000·(8² − 2²) = 30,000 Pa = 30 kPa. A carburettor, a perfume spray and a Venturi flow meter all run on this. Where it breaks. Bernoulli assumes steady flow, no friction, and no change in density. Friction in a long pipe steals pressure. Air above roughly Mach 0.3 (about a third of the speed of sound) starts to compress. And if the pressure falls far enough, water boils at room temperature (cavitation), which wrecks propellers and pumps. Interactive · The Venturi Speed in the wide pipe (10 cm) Throat diameter Pipe 10 cm to 5 cm at 2 m/s: 8 m/s in the throat, pressure down by 30 kPa. Water, 1,000 kg/m³, ideal flow, 150 kPa (absolute) before the throat. Squeeze the throat hard enough and the pressure falls to the boiling pressure of water. Lift: why "equal transit time" is wrong You may have been told that air going over the curved top of a wing has further to go, so it must go faster to "catch up" with air going underneath, and Bernoulli does the rest. The last step is right: the pressure on top is lower. The first steps are wrong. There is no rule that two neighbouring bits of air must meet again at the back edge, and they don't: in measurements and simulations the air over the top gets there first. Worse, the "longer path" story predicts far too little lift, and it can't explain a flat plate, a symmetric wing or a plane flying upside down, all of which fly. NASA Glenn's education pages call it "one of the most widely circulated, incorrect explanations". The better picture: a wing, tilted into the wind (the angle of attack), turns the air downward. Pushing air down pushes the wing up, which is Newton's third law. The pressure difference between top and bottom is how that push reaches the wing's skin, so Bernoulli and Newton are two bookkeepings of one event, not rivals. Why does the air go faster on top? Because the wing sets up a circulation of flow around itself, fixed by the sharp trailing edge, not by any race. L = ½ · ρ · v² · S · CLlift = dynamic pressure × wing area × lift coefficient For small angles the coefficient CL grows in proportion to the angle, until the flow can no longer follow the top of the wing. Then it separates and lift collapses: a stall. Worked case: a 1,100 kg light plane with 16 m² of wing at 50 m/s (97 knots) needs CL = 0.44, about 5° above its zero-lift angle. Slow to 30 m/s and the same plane needs 14°, almost at the stall. This is why landing speeds are what they are. Interactive · Hold the plane up Airspeed Angle of attack (from zero lift) Find the level-flight angle At 50 m/s and 5.1° the lift of 10.8 kN balances the weight of 10.8 kN. The stall is at about 28 m/s. A model, not a real aircraft: wing area 16 m², aspect ratio 7.5, sea-level air, lift slope about 5 per radian (thin-airfoil theory with a finite-wing correction), stall assumed at 16°. Viscosity and the Reynolds number: sticky against sluggish Viscosity is a fluid's internal stickiness: how hard layers resist sliding over each other. Honey is far more viscous than water; water is about 55 times more viscous than air. Whenever a fluid moves, two things compete. Inertia is the tendency of the moving fluid to keep going, and viscosity is the tendency to smooth the motion out. The ratio of the two is one number. Re = ρ · v · L / μReynolds number = density × speed × size, divided by viscosity Osborne Reynolds showed in 1883, by injecting a thread of dye into water flowing in a glass pipe, that this one number decides how a flow behaves. Low Re: the dye thread stays a straight line (laminar flow). High Re: it breaks into chaotic swirls (turbulent flow). In a pipe the change happens roughly between Re 2,000 and 4,000; it is a range, not a switch. The same number explains why a swimming bacterium can't coast. For a 2 micrometre bacterium swimming at 30 micrometres a second, Re is about 6 × 10⁻⁵, so the water behaves like thick syrup and it stops the instant it stops pushing (Purcell's classic 1977 essay "Life at low Reynolds number"). A human swimmer is near 4 million. Interactive · What is the Reynolds number of...? BacteriumSlow water in a 2 cm pipeSame pipe, 10× fasterWalking personThrown baseballSwimmer WaterAir Size L Speed v Water at 0.1 m/s in a 2 cm pipe: Re ≈ 2,000, right at the edge of laminar flow. Ten times faster is about 20,000: turbulent. Water 20 °C (998 kg/m³, 1.0 mPa·s); air 20 °C (1.2 kg/m³, 0.018 mPa·s). The labels follow pipe flow; around other bodies the thresholds differ. Boundary layers and drag: the thin skin that decides everything At the surface of any solid, fluid sticks to it (the no-slip condition). So speed rises from zero at the wall to the full stream speed over a thin skin called the boundary layer. On a flat plate, 0.3 m from the front edge in a 10 m/s breeze, that skin is only about 3 mm thick. Almost all the friction, and a lot of the drama, lives there. Drag comes in two kinds. Friction drag is the skin rubbing. Pressure drag appears when the boundary layer peels off the body behind its widest point (separation), leaving a low-pressure wake. Streamlining means keeping the flow attached so the wake stays thin. The whole story fits one formula. F = ½ · ρ · v² · Cd · Adrag = dynamic pressure × drag coefficient × frontal area Because it goes with the square of speed, the power to overcome it (force times speed) goes with the cube. A car with Cd = 0.30 and 2.2 m² of frontal area needs about 8.5 kW to push air aside at 100 km/h and 18.7 kW at 130 km/h: 2.2 times the power for 30% more speed. A golf ball shows the other trick. Its dimples trip the boundary layer into turbulence early; a turbulent layer clings to the ball longer, the wake shrinks, and drag falls (measured drops of more than a factor of two are reported for dimpled over smooth spheres at the same Re). A rougher surface can mean less drag. Interactive · Drag, force and power Flat plate face-on (1.2)Sphere (0.47)Car (0.30)Streamlined teardrop (0.05) Speed Car at 100 km/h: drag about 307 N, about 8.5 kW. At 30% more speed: 1.69× the force, 2.20× the power. Frontal area 2.2 m², air 1.204 kg/m³. The coefficients are typical round values and vary with the exact shape and with Re. Navier–Stokes: Newton's second law for a blob of fluid Take a tiny parcel of fluid and apply F = ma. That's it. Claude-Louis Navier wrote a version with viscosity in 1822, and George Stokes re-derived it on firmer ground in 1845. For water or slow air (incompressible flow), in words and then in symbols: ρ (∂u/∂t + u·∇u) = −∇p + μ∇²u + ρgmass × acceleration = pressure push + viscous smearing + gravity ρ(∂u/∂t + u·∇u): mass times acceleration. The parcel accelerates for two reasons: the flow at its spot changes in time (∂u/∂t), and it is carried to a place where the flow is different (u·∇u). −∇p: the parcel is pushed from high pressure toward low. μ∇²u: viscosity averages the velocity with the neighbours, so fast fluid drags slow fluid along. It smooths and dissipates, turning motion into heat. ρg: gravity (or any body force). Plus a second rule, ∇·u = 0, which is continuity. The troublemaker is u·∇u: the velocity appears multiplied by its own slope. This nonlinear term lets big eddies feed small ones, and it is why turbulence is so hard. Its job shows up already in a one-dimensional cousin, Burgers' equation, below. Try it: the nonlinear term steepens a smooth wave into a front, and viscosity decides how sharp the front can get. Interactive · Steepening against smoothing Viscosity ν Time Burgers' equation (fast parts catch slow parts, viscosity smooths). Without viscosity the front would become infinitely sharp at t ≈ 0.16. Computed live (finite differences, 400 points) and checked against the exact solution: errors under 1%. It is a 1-D toy, with no pressure and no turbulence. The million-dollar question. In 2000 the Clay Mathematics Institute named Navier–Stokes as one of seven Millennium Prize Problems, with $1,000,000 for each. The question is not whether the equations are useful (they are, endlessly). It is whether, in three dimensions, smooth starting flows always stay smooth, or whether a solution can ever "blow up" into infinite velocities in finite time. Nobody knows. As of this writing no proof has been accepted by the Clay Institute; claims appear regularly and are checked slowly. Dimensional analysis, wind tunnels and CFD You can't solve Navier–Stokes by hand for a car. So engineers do two things. First, test a small model. The trick, from dimensional analysis, is that two flows behave alike when their dimensionless numbers match, mostly the Reynolds number. A 1:10 model car must therefore run ten times faster in air, which would be 300 m/s, almost Mach 0.9: not the same physics at all. The fix is to change the fluid: water has about 15 times less kinematic viscosity than air, so the model only needs about 20 m/s underwater. Interactive · The scale model Air tunnelWater tunnel Model scale Real car at 30 m/s; a 1:10 model in water must run at 19.9 m/s to match the Reynolds number. Kinematic viscosity: air 1.516 × 10⁻⁵ m²/s, water 1.004 × 10⁻⁶ m²/s (20 °C). Speed of sound: air 343 m/s, water 1,480 m/s. Second, computational fluid dynamics (CFD): slice space into millions of cells and march the equations forward on a computer. Resolving every eddy (direct numerical simulation) costs roughly the cube of the Reynolds number, which is why aircraft and cars are simulated with turbulence models that average the small eddies away. CFD is numerical analysis and parallel computing as much as fluid mechanics, and it is where Computer Science meets this course. Modern work also trains machine-learning models on CFD results (Machine Learning, 34/32) to get approximate answers fast. In practice teams use CFD, tunnel tests and flight tests together, because each checks the others. What a graduate just sees Pressure is weight above you, so deep water pushes hard: depth, not volume, loads a dam. Squeeze a flow and it speeds up and its pressure drops. Check whether the "ideal" assumption still holds. A wing works by turning air downwards. "Equal transit time" is a myth. Lift and drag grow with speed squared; the power to overcome drag grows with speed cubed. Before anything else, ask for the Reynolds number: it tells you whether viscosity or inertia is in charge. Thin boundary layers carry the friction, and separation of the layer causes most drag. A rougher surface can have less drag (golf ball), and slower water can match faster air (model testing). The Navier–Stokes equations are Newton's second law with pressure, viscosity and gravity, and nobody can yet prove they always stay smooth. Where it shows up next Heat Transfer (34/16): a moving fluid carries heat. The boundary layer you met here is also the thermal layer that decides how fast a radiator or a chip cools. Propulsion (34/21): jets and rockets are pipes, nozzles and compressors where Bernoulli, compressibility and thermodynamics meet. Vibrations (34/17) and Control Systems (34/20): wind makes structures flutter, and a flight controller exists to keep a wing at the right angle of attack. The real world: aircraft and car design, wind turbines, pumps, pipelines, blood-flow models, weather forecasting. Computer Science: CFD codes are large parallel programs, and fluid solvers drive the water in films and games. If you want to go deeper 📘 Textbook: Frank M. White, Fluid Mechanics (McGraw-Hill): a standard undergraduate text, strong on worked examples. 📘 Textbook: Kundu, Cohen and Dowling, Fluid Mechanics (Academic Press; 7th edition, 2024): broader, with more on the mathematics and CFD. 📘 Gentler: Steven Vogel, Life in Moving Fluids (Princeton University Press): fluid mechanics through animals and plants, nearly no equations. ▶️ Free course: MIT OpenCourseWare 2.06 Fluid Dynamics (Spring 2013), an undergraduate course with problem sets and exams. 🌐 Easier first read in this library: Flying Machines (29/05) and Chaos Theory and the Butterfly Effect (12/07). The takeaway A fluid is a pile of tiny parcels obeying F = ma with three forces: pressure, viscosity and gravity. Pressure is the weight above, so dams are thick at the bottom. Continuity and Bernoulli say squeezed flow speeds up and loses pressure. Lift comes from turning air downward, not from equal transit times. The Reynolds number says whether a flow is smooth or turbulent, and drag grows with the square of speed. The full equations are the best rulebook we have, and whether they always behave is still an open million-dollar question. Try this: next time you drive or ride, notice how much harder the engine works at motorway speed than at 50 km/h. That is the cube law. Then hold a spoon's back under a thin tap stream and feel the water pulled around the curve. Sources & Further Reading Sources. Hydrostatics, continuity, Bernoulli, Navier–Stokes and dimensional analysis are standard textbook material (White; Kundu et al.). Navier 1822 and Stokes 1845: the 2024 review "200 Years of the Navier–Stokes Equation" (arXiv 2401.13669). Millennium Prize: Clay Mathematics Institute (claymath.org), official problem by C. Fefferman. Equal-transit myth and the downward-turning explanation: NASA Glenn Research Center, "Incorrect Lift Theory" and Tom Benson's Beginner's Guide to Aeronautics. Reynolds's dye experiment (1883) and the pipe transition range of about 2,000–4,000. Low-Re swimming: E. M. Purcell, "Life at low Reynolds number", American Journal of Physics 45 (1977). Golf-ball drag reduction: Balaras et al., APS Division of Fluid Dynamics 2016. Computed in node: pressure, force, Venturi, lift, Reynolds numbers, drag power and scale-model speeds; the Burgers solver was checked against the exact Cole–Hopf solution. 🌐 MIT OpenCourseWare 2.06 Fluid Dynamics 🌐 Clay Mathematics Institute: Navier–Stokes equation 🌐 Purcell, Life at Low Reynolds Number (a reprint) 📘 Steven Vogel, Life in Moving Fluids Step 13. Structural Analysis: How Forces Find Their Way to the Ground --------------------------------------------------------------------- Why: How forces find their way down through a structure to the ground. A structure is a stack of paths from the roof to the soil. Learn to follow the load and you can see which parts are working, which are spare, and what happens when one goes. On the 80th floor of a tall building on a windy day, a glass of water on a desk can ripple even though nothing is wrong. The building is swaying, slowly and by a few centimetres, the way a tree does. Nobody designed it to be still. Someone decided how far it would be allowed to move, and checked that the people on the 80th floor would not be seasick, that the lifts would still run, and that the cladding would not crack. That decision comes from structural analysis: given a structure and the loads on it, find every force, every bending moment and every movement inside it. Statics (34/10) gave you the forces in a structure that is simple enough to solve by balance alone. Real structures are not that simple. They are continuous frames with more members than the minimum, loaded by wind, earthquakes and lorries as well as their own weight, and the answer has to say not just whether they stay up but how much they move. The course in one box Who takes it: civil and structural engineers, usually in the third year, as a first course (determinate structures) and a second (indeterminate structures and matrix methods); aerospace engineers meet the same ideas for airframes. Programmes that follow the ASCE Civil Engineering Body of Knowledge put it ahead of structural design.What it is for: working out the forces, bending moments and movements in a structure so that a designer knows what each part has to carry.Comes after: Statics (34/10), Mechanics of Materials (34/11) and Linear Algebra (34/04).Leads to: Structural Design: Concrete & Steel (34/65), which sizes the members this course gives forces for; and the earthquake side of Vibrations (34/17).What this piece will and won't give you: you won't be able to analyse a real building after this; nobody does it by hand any more. You will see why a loaded structure sends its weight down a path, why spare members save lives, and what the computer does when it solves a frame. Easier first reads: Statics: Why Things Stand Still (34/10) and Mechanics of Materials (34/11). The loads, spans and steel section here are illustrative round numbers, not design values. Loads find a path to the ground, and they never arrive alone Every force on a building has to reach the ground by a continuous route, called the load path. A floor slab spreads weight to its beams; the beams pass it to girders; the girders pass it to columns; the columns pass it to footings; the footings pass it to the soil (34/66). Break the chain anywhere and the load looks for another route or the structure falls. Engineers read a building by tracing these paths, and count how much floor each column is responsible for, its tributary area. An interior column in a 6 m by 6 m grid under ten floors of 5 kPa carries 10 × 36 × 5 = 1,800 kN at its base, about 180 tonnes. The loads have different characters. Dead load is the structure itself, permanent and well known. Live load is people, furniture and stored goods, which come and go. Wind pushes sideways and its effect grows with height. Earthquake does not push at all: the ground moves, and the building's own mass resists the movement, so the force is mass times acceleration and a heavier building attracts a bigger force. That is why a heavy building is often governed by earthquake and a light, tall one by wind. Interactive · Gravity, wind and earthquake on one tower Add storeys and see which sideways load takes over. The plan is 30 m by 30 m; the earthquake force is a fraction of the weight. StoreysBay sizeFloor loadWind pressureEarthquake fraction of weight These loads never act alone and never all peak together. Codes give load combinations: a recipe of factors that multiplies each load by how uncertain it is and how likely it is to coincide with the others. Two examples of the shape: ASCE 7 in the United States uses 1.2 D + 1.6 L, and the Eurocodes use 1.35 G + 1.5 Q (G is permanent load, Q variable). Here is the surprise: another combination takes only 0.9 of the dead load and subtracts the wind. When wind lifts a light roof, the dead weight is helping, so you must assume you have less of it. The factors above are examples, not design values. Interactive · Which combination governs? Try a heavy dead load, then a light roof with a big wind. Switch between two code styles. ASCE 7 styleEurocode styleDead load effectLive load effectWind effect Redundancy: spare load paths forgive mistakes A structure is statically determinate when the equilibrium equations alone give every force, with no member to spare: remove any one and it becomes a mechanism. It is indeterminate when it has more members or supports than that. For a flat pin-jointed truss, count members m, support reactions r and joints j. If m + r = 2j it is determinate. The truss in Statics has 7 members, 3 reactions and 5 joints: 7 + 3 = 10 = 2 × 5. Add one extra diagonal and the sum exceeds 10 by one: one degree of indeterminacy, one spare load path. Spare paths are what you want when something goes wrong, and something does. On 15 December 1967 the Silver Bridge over the Ohio River, between Point Pleasant, West Virginia and Gallipolis, Ohio, collapsed in rush-hour traffic and 46 people died. A single eyebar in one suspension chain failed: a crack about 2.5 mm deep, grown by stress corrosion, broke it. Each link of the chain had only two eyebars, so when one went the other was overloaded. Investigators concluded that with three or more eyebars per link, one failure might not have been fatal; the crack was beyond what 1960s inspection could find. The collapse led to the first national bridge inspection programme in the Federal-Aid Highway Act of 1968. The price of redundancy is that indeterminate structures are sensitive to things statics does not see: if one support settles, or one member heats up, the forces redistribute (the soil enters the problem again, 34/66). Interactive · Cut a hanger from a bridge A deck hangs from identical hangers. Cut some, and see whether the rest can take it. Static sharingSudden loss, load doubledHangersHangers cutDeck load Beams and frames: the moment diagram is the picture of the structure The bending moment at each point of a beam is the quantity that decides how big the beam must be (34/11), so the central drawing of structural analysis is the moment diagram: the moment plotted along the member. For a simply supported beam under a uniform load w, the moment is a parabola with a peak of wL²÷8 at midspan. For a beam fixed at both ends, the ends hold it: the peak moment at each end is wL²÷12 and at midspan only wL²÷24. In numbers, for 15 kN/m on 6 m: 67.5 kN·m simply supported against 45 at the ends and 22.5 at the middle if fixed. Same beam, same load, one third less peak moment, and much less deflection, because indeterminate beams are stiffer. δ = 5wL⁴ ÷ 384EI   (simply supported, uniform load) Read it as: sag grows with the fourth power of the span, so doubling the span makes the sag sixteen times larger unless the beam is made stiffer. A frame is beams and columns joined by stiff joints; the joint passes moment between them, which is why a beam's end moments land in the columns and why the column sizes depend on the floors above. Codes limit deflection as a fraction of the span, with values such as span ÷ 360 under live load common for floors; the actual limit depends on the code and what the floor carries. The steel beam here sags 15.1 mm over 6 m, which is span ÷ 396. Interactive · Same beam, different supports Change how the ends are held and watch the moment diagram and the sag change. SimpleFixed-fixedFixed + rollerCantileverSpanUniform load Influence lines: ask what one load does, wherever it stands A bridge carries a moving load, so no single picture of forces is enough. The trick is to turn the question round. Put a single unit load at every position along the span and plot what it does to one quantity, such as the moment at midspan. That plot is an influence line. For a simply supported span the midspan moment is a triangle: zero with the load at either support, largest with it in the middle. Then for any set of loads, the effect is the sum of each load times the height of the line under it. Take a truck with three axles of 50, 150 and 150 kN spaced 4 m apart, crossing a 20 m span. The worst midspan moment is 1,350 kN·m, reached when the lead axle is 14 m along and a heavy axle sits at midspan. Parking all 350 kN at midspan would give 1,750 kN·m, so a truck is kinder than its total weight because its axles spread out. Real codes use design vehicles and lane loads of this kind (the AASHTO HL-93 loading in the United States or Load Model 1 in the Eurocodes, for example), and the maximum is found with influence lines. The worst place for the truck is not where you would guess. Interactive · Drive a truck across the bridge Drag the position of the truck. The green curve is the midspan moment as it crosses. SpanLead axle position The stiffness method: stiff things attract force Indeterminate structures were once solved by hand with clever iterative tricks, such as Hardy Cross's moment distribution (1930). Computers need something more mechanical, and the answer is to treat every member as a spring. A bar of area A, modulus E and length L is a spring with stiffness k = EA÷L. A beam is a spring for bending, stiffness proportional to EI over length. Join the springs at nodes, write the force balance at each node in terms of unknown displacements, and you get one equation: K u = P Read it as: the stiffness matrix times the displacements gives the applied loads. Solve for the displacements, then every member force follows from its own stretch times its own stiffness. A stiffer member moves the same distance but carries more force, so stiff things attract load. The finite element method, FEA, is exactly this with the structure cut into thousands or millions of small elements. Only the size of K changes. Try three bars between two walls: stiffnesses 100, 50 and 100 kN/mm, with a 100 kN pull on the first node. Equilibrium alone cannot say how the load splits, but K = [150, −50; −50, 150] does. The nodes move 0.75 mm and 0.25 mm. Bar 1 is pulled with 75 kN, bars 2 and 3 are pushed with 25 kN each. Interactive · Solve a structure with a matrix Change the stiffness of each bar and the loads on the nodes. The 2 by 2 matrix is solved live. Bar 1 stiffnessBar 2 stiffnessBar 3 stiffnessLoad on node 1Load on node 2 Tall buildings: a vertical cantilever fighting the wind A tall building is a cantilever sticking out of the ground, loaded sideways by wind. For a uniform sideways load w, the top moves by wH⁴÷8EI: fourth power of the height again. So the whole art of tall-building structure is to raise I, and I grows with the square of how far material sits from the centre line. Material in the middle is nearly useless for sway; material at the edge is worth far more. That gives the three systems. A core, a hollow tube of concrete or steel walls round the lifts and stairs, is the spine. A tube puts the stiff material on the outside wall: closely spaced columns joined by deep beams, as in the framed, braced and bundled tubes developed by Fazlur Khan in the 1960s and 70s (examples are the John Hancock Center, Chicago, and the Willis Tower). Outriggers are stiff arms or trusses that tie the core to the perimeter columns, usually at a few floors, so that the outer columns pull and push like the flanges of a giant I-beam when the core tries to bend. Where the structure alone cannot calm the sway, a tuned mass damper (a heavy pendulum, 660 tonnes in Taipei 101) takes the energy out. Interactive · Stiffen a tower Start with the core alone. Then engage the perimeter columns, as outriggers or a tube would, and watch the top. HeightWidthPerimeter engaged What a graduate just sees A structure is a set of paths that carry load to the ground; tracing the path is how you read it. Loads come in families, and codes combine them with factors. The worst case is not always the heaviest one: wind uplift beats dead weight. Count members, reactions and joints. A determinate structure has none to spare; an indeterminate one has spare paths. The moment diagram is the beam's story. Fixing the ends lowers the peak moment and the sag. Sag grows with the fourth power of span and falls with EI. An influence line turns "where is the worst place for the load?" into reading a graph. Everything is a spring. Stiff members attract force, and the computer solves K u = P. To stop a tall building swaying, put the stiff material as far from the centre as you can. Where it shows up next Structural Design: Concrete & Steel (34/65): this course gives the forces and moments; design sizes the section and the bars for them, using the load combinations shown here. Geotechnical Engineering (34/66): the supports in every model are soil, and a support that settles changes the forces in an indeterminate frame. Vibrations (34/17): earthquake and wind response of a building is its mass and stiffness as a spring-and-mass problem. Linear Algebra (34/04): K u = P with a million unknowns is a sparse linear system. In the world: the structural model behind every skyscraper, bridge rating, stadium roof and aircraft wing, and the FEA package an engineer runs before anything is built. If you want to go deeper 📘 Standard textbook: R. C. Hibbeler, Structural Analysis (Pearson). Widely assigned, with many worked problems on beams, frames, influence lines and matrix methods. 📘 Standard textbook: Aslam Kassimali, Structural Analysis (Cengage). Another common course text, clear on the stiffness method. 📘 Gentler: Mario Salvadori, Why Buildings Stand Up: The Strength of Architecture (W. W. Norton, 1980). Load paths, arches, frames and skyscrapers for the general reader. And J. E. Gordon, Structures: Or Why Things Don't Fall Down. 🌐 Free course: MIT OpenCourseWare 1.050 Solid Mechanics, which covers determinate planar structures and statically indeterminate systems. For a graduate step, 1.571 Structural Analysis and Control. The takeaway A structure is a set of load paths to the ground. Loads arrive in families and combine under code factors; a structure with spare paths forgives mistakes that a determinate one does not; the moment diagram shows what each member must resist; an influence line finds the worst place for a moving load; and a computer solves any of it by treating each member as a spring, so that stiff things attract force. A tall building is then a cantilever, and the cure for sway is to move the stiff material outwards. Try this: next time you cross a bridge or stand in a doorway, find the load path. Where does the weight go, which members are in tension, and could the structure survive losing one of them? Sources & Further Reading Sources. Equilibrium, the determinacy count (m + r = 2j), the moment diagrams and deflection formulas for a uniformly loaded beam (5wL⁴/384EI, wL²/8, wL²/12, wL²/24, wL⁴/185EI, wL⁴/8EI), influence lines, the direct stiffness method and cantilever drift (wH⁴/8EI) are standard textbook structural analysis. Every number in the interactives was computed in node: column load 10 × 36 × 5 = 1,800 kN; base shears for 10 storeys (earthquake 4,500 kN at 0.1 of weight, wind 1,050 kN); the combination examples (1.2D + 1.6L = 400 kN for D = 200, L = 100); the beam moments and deflections (EI = 16,720 kN·m², E = 200 GPa and I = 8.36×10⁻⁵ m⁴ about an IPE 300 section: 15.1 mm, span ÷ 396); the truck (1,350 kN·m at lead axle 14 m; 1,750 for a point load at midspan); the three-bar chain (u = 0.75 and 0.25 mm; forces 75, −25, −25 kN); and the tower (Core I = 417 m⁴ for a 12 m core with 0.4 m walls, tube I = 10,468 m⁴ with 0.25 m equivalent wall, E = 30 GPa and 1.5 kPa wind, all illustrative; real tubes lose part of their stiffness to shear lag). Silver Bridge facts (15 December 1967, 46 dead, eyebar 330, crack about 2.5 mm, two eyebars per link, 1968 inspection programme) come from the Wikipedia article on the Silver Bridge. The factors 1.2/1.6 and 1.35/1.5, the HL-93 and Load Model 1 names, the span ÷ 360 limit, Hardy Cross 1930, Fazlur Khan's tubes and the 660-tonne Taipei 101 damper are standard engineering knowledge not re-checked against a primary source in this session. Go deeper. Hibbeler, Kassimali, Salvadori and Gordon are standard titles; check the current edition. The MIT OpenCourseWare pages were confirmed to exist. 🌐 MIT OpenCourseWare 1.050 Solid Mechanics 📘 Mario Salvadori, Why Buildings Stand Up 📘 J. E. Gordon, Structures: Or Why Things Don't Fall Down Step 14. Geotechnical Engineering: Everything Stands on Soil ------------------------------------------------------------ Why: Everything stands on soil, and soil is the least predictable material there is. Soil is the one structural material nobody manufactured. Learn how water, weight and shaking decide whether the ground holds. At the beach, walk on the wet sand just above the waves and it is firm enough to run on. Stand still on it and bounce a little, and a shine of water creeps up around your feet and the ground goes soft. Dry sand a metre higher is loose and slow. Sand under water is different again. It is the same grains every time; what changed is where the water is, and how fast you shook it. That is the whole problem of geotechnical engineering. Steel is made in a mill to a specification. Concrete is mixed to a recipe. Soil is whatever the last glacier, river or volcano left, it differs from one corner of a site to the next, and every building, bridge and dam in the world finally rests on it. Sometimes the ground says so loudly. The Tower of Pisa began to sink while it was still being built, in 1178, because its shallow foundation, about three metres deep, sat in weak subsoil. Work stopped for most of a century, partly through war, and that pause let the soil settle and gain strength, which is a large part of why the tower is still standing. It is the most famous lesson in this subject, and we will come back to it. The course in one box Who takes it: civil engineers, usually in the third year, as a first course in soil mechanics followed by a second in foundation engineering; geological, mining and environmental engineers take it too. Programmes that follow the ASCE Civil Engineering Body of Knowledge treat geotechnical engineering as one of the core areas of the degree.What it is for: working out how much weight soil can carry, how far it will sink, and when a slope or a bank will move.Comes after: Mechanics of Materials (34/11) and Fluid Mechanics (34/15).Leads to: Structural Design: Concrete & Steel (34/65), whose footings and piles rest on this; Hydraulics & Water Resources (34/67), for dams and seepage; and Construction Management (34/70), where excavations and temporary works live.What this piece will and won't give you: you won't be able to design a foundation after this, and no engineer would trust a number from a toy. You will see why wet soil is weaker than dry, why clay sinks for decades, and why an earthquake can turn a street into soup. Easier first reads on the same ground: Statics: Why Things Stand Still (34/10), Mechanics of Materials (34/11) and Everyday Physics: The Answers Hiding in Plain Sight (Physics 5). All soil values below are illustrative round numbers, not design values. Soil is grains, water and air, and the water carries part of the load A lump of soil is three things at once: solid grains, water, and air. The grains touch at their corners and the gaps between them, the voids, are filled with water, air or both. Engineers start every problem by counting the three. Take a soil with a void ratio of 0.7 (voids are 70% of the volume of the solids). Four-tenths of its volume is empty space. Fill the voids with water and it weighs about 19.3 kN per cubic metre instead of 15.3 when dry, a quarter more. unit weight γ = (Gs + S·e) ÷ (1 + e) × γw Read it as: the weight of the grains plus the weight of the water, divided by the total volume. Here e is the void ratio, S the fraction of voids holding water, Gs about 2.65 for quartz grains, and γw = 9.81 kN/m³ for water. Interactive · Count the three phases Change how loosely the grains are packed and how much of the gap is filled with water. Void ratio eSaturation S The key discovery of the subject, by Karl Terzaghi in the 1920s, is what the water does under load. Put weight on soil and two things can carry it: the grain-to-grain contacts and the pressure of the water in the voids. Only the contacts give friction, and friction is what soil strength is. Terzaghi's principle says that soil behaves according to the part carried by the grains, the effective stress: σ′ = σ − u Read it as: effective stress is the total stress minus the water pressure. At 10 m depth with the water table at 2 m, the total stress is about 196 kPa and the water pressure 78 kPa, so the grains carry only 118 kPa. Raise the water table and the grains carry less. Push water upward through sand fast enough and the effective stress falls to zero, the sand floats, and a cofferdam floor boils. This is "quicksand", and it needs flow, not a special kind of sand. Interactive · Who carries the load, the grains or the water? Move the point deeper, raise the water table, or push water upward through the soil. Depth of the pointDepth of the water tableUpward flow gradient How fast water moves is its own measurement. Darcy's law says the velocity is the permeability times the hydraulic gradient, and permeability differs by many orders of magnitude between clean gravel and clay. That one fact, how slowly water leaves, is what the next idea is about. Clay sinks for decades: consolidation Squeeze a wet sponge and the water leaves at once. Squeeze a sponge whose holes are a thousand times smaller and it takes a long time. A clay layer under a new load behaves like the second sponge. At first the water takes the extra load, because it cannot escape. As it seeps out, the load moves to the grains, which pack closer, and the ground settles. That slow transfer is consolidation. Terzaghi gave it two equations. The amount of final settlement depends on how compressible the clay is and how much the load is. The time it takes depends on how far the water has to travel: t = Tv · Hdr² ÷ cv Read it as: time grows with the square of the drainage path. Twice the distance means four times the wait. A 10 m clay layer that can drain only upward needs Tv = 0.848 for 90% of its settlement; with cv = 2 m²/yr that is 42 years. If sand lets it drain both up and down, the path is 5 m and the wait is 10.6 years. Interactive · Watch a clay layer settle Pick the load, the clay, and whether water can leave from one side or both. Drag the time. Drains one wayDrains both waysClay thicknessAdded load ΔσConsolidation coefficient cvTime This is Pisa. Construction stopped around 1178 and resumed in 1233 and 1272, and the long pauses let the soil under it consolidate. The tilt reached 5.5 degrees by 1990, and the tower was closed on 7 January 1990. Remedial work from 1993 to 2001 used about 870 tonnes of lead counterweights and removed soil from beneath the high side, about 70 tonnes in all. The tilt came down to 3.97 degrees, by about 45 cm at the top, back to where it had been in 1838. The tower reopened on 15 December 2001. Notice that engineers did not strengthen the tower; they fixed the ground, and they did it by taking soil away. Soil strength is friction, and friction needs squeezing Pile dry sand and it stops at a steady angle, the angle of repose, whatever the size of the pile. That angle is a picture of friction between grains. Soil shear strength follows the same law as a block on a ramp (34/10), with one change: the squeeze that counts is the effective stress. τf = c′ + σ′ tan φ′ Read it as: strength is a little cohesion plus the effective squeeze times a friction factor. Sand has no cohesion to speak of, so for a thin layer sliding on a slope the factor of safety is simply tan φ′ ÷ tan β. A dry sand slope with φ′ = 32° stands at any angle under 32°. At 25° the factor of safety is 1.34. Now let rain fill the layer. The water pressure cancels part of the squeeze, and the factor of safety drops to 0.65. The slope that stood now slides, though no one touched it. Interactive · Steepen a slope, then add rain Raise the angle until it fails. Then fill the soil with water and see how far the safe angle falls. Slope angleFriction angle φ′Water in the slip layerCohesion c′ A foundation can only push as hard as the ground can push back A footing spreads a column's weight over a patch of soil. Load it more and more and, in the end, the soil under it shears along a bowl-shaped surface and the footing plunges. The pressure at that point is the ultimate bearing capacity. The bearing-capacity equation for a long strip footing, in Terzaghi's form, has three terms, one for each way the soil helps: qult = c Nc + q Nq + ½ γ B Nγ Read it as: cohesion, plus the weight of the soil beside the footing (the deeper you bury it, the more it helps), plus the weight of the soil under the footing (a wider footing mobilises more of it). The N values depend only on the friction angle (here the classical Prandtl–Reissner values, and Vesic's for Nγ). Engineers then divide by a safety factor, commonly 3, and check that the actual pressure is below it. A footing 2 m wide, 1 m deep in loose sand (φ = 30°) has an ultimate capacity near 690 kPa. Flood the site, and it falls to about 290. If the good ground is deep, no spread footing will do. Engineers drive or drill piles through the weak layer so that they carry the load by friction along their sides, by bearing at their tips on something firm, or both. Interactive · Size a footing Pick the ground, the footing, and the load per metre of footing length. Then flood the site. Soft clayLoose sandDense sandFloodedFooting width BBurial depthColumn load Shaking can turn sand into liquid Loose, water-filled sand wants to pack tighter when shaken. If the water cannot leave in time, it takes the load instead, the water pressure climbs, and the effective stress falls. Strength is effective stress times friction, so strength falls with it. When the pressure of the water has risen to equal the squeeze the grains were under, the effective stress is zero and the sand has no strength: this is liquefaction. It happens in loose sandy and silty ground below the water table; clay and dense gravel mostly resist. The Christchurch earthquake of 22 February 2011 was a magnitude 6.2 earthquake, centred about 6.7 km south-east of the city centre, and 185 people died. Liquefaction caused extensive damage to homes in the eastern suburbs, covering streets and gardens with silt and sand; the figures for how much silt was removed vary between sources (several hundred thousand tonnes). Some of the worst-hit areas were designated a residential red zone, judged not feasible to rebuild. The lesson for engineers is that a site can be perfectly strong on a calm day and have almost no strength for thirty seconds. Interactive · Raise the water pressure Slide from calm ground to full liquefaction in sand at 5 m depth (water table at 1 m). Excess water pressure Soft ground amplifies the shaking On 19 September 1985 a magnitude 8.0 earthquake struck off the Pacific coast of Michoacán, more than 350 km from Mexico City. By the time the waves reached the city they were far weaker than at the coast, yet the damage was terrible: 412 buildings collapsed in the city and thousands more were seriously damaged. Estimates of the dead range widely, from about 5,000 to far more, with about 10,000 the figure most often given. The cause was the ground. Much of the city sits on the dried bed of Lake Texcoco, soft, water-saturated volcanic clay 7 to 37 m deep, and a soft layer on stiff rock acts like a tuning fork. T = 4 H ÷ Vs Read it as: the layer's natural period is four times its thickness divided by the speed of shear waves in it. For 37 m of clay with Vs = 74 m/s that is 2 s. The lakebed's own period is commonly put at about 2 to 2.5 s, which matched the natural period of buildings of 6 to 15 storeys, and those suffered most. Resonance (Vibrations, 34/17) did the rest. A simple rule gives a building a period of about 0.1 s per storey, shorter than the ground's, but a cracked and yielding building gets more flexible, and its period stretches towards the ground's. Interactive · Tune a building to the ground Press the lakebed button, then try buildings of different heights. Slide the damage control to see a cracked building drift into resonance. Lakebed-like groundFirm groundThickness of soft layerShear-wave speedBuilding heightPeriod lengthening What a graduate just sees Soil is grains, water and air, and the water carries part of the load. Wet is weaker than dry because only the grains give friction. Effective stress, σ′ = σ − u, is the stress that matters. Anything that changes the water pressure changes the strength. Clay is a sponge with tiny holes: it settles for years, and the wait grows as the square of the drainage path. Sand has friction but no glue, so a slope's safe angle is its friction angle, and rain lowers it a lot. A footing's capacity rises with width, depth and friction, and falls when the ground floods. Loose saturated sand shaken fast has no time to drain. It loses strength and flows, and a buried tank floats up. Soft layers have a natural period, 4H/Vs, and shaking at that period is magnified. The usual fix for bad ground is to remove the water, spread the load, or go through it to something better. Where it shows up next Structural Design: Concrete & Steel (34/65): every footing and pile cap in that course carries the load found here, and the settlement a building can tolerate decides how it is built. Hydraulics & Water Resources (34/67): earth dams, levees and seepage under them are Darcy's law plus effective stress. Structural Analysis (34/64): a structure's supports are never perfectly rigid, and a support that settles changes the forces in an indeterminate frame. Vibrations (34/17): site amplification is resonance, with the soil as the spring. In the world: tunnels, excavations, retaining walls, offshore wind-turbine foundations, landslide hazard maps, and the hidden half of every building project. If you want to go deeper 📘 Standard textbook: Braja M. Das, Principles of Geotechnical Engineering (Cengage). A common first course text with many worked problems. 📘 Standard textbook: Muni Budhu, Soil Mechanics and Foundations (Wiley). Another widely assigned text, with an emphasis on worked examples. 📘 The classic: Karl Terzaghi, Ralph Peck and Gholamreza Mesri, Soil Mechanics in Engineering Practice (Wiley). The practitioner's book by the founder of the subject, for after the first course. 📘 Gentler: Henry Petroski, To Engineer Is Human: The Role of Failure in Successful Design (St. Martin's Press, 1985). Case histories, no equations. 🌐 Free course: NPTEL: Soil Mechanics/Geotechnical Engineering I, Prof. Dilip Kumar Baidya, IIT Kharagpur, free video lectures. The takeaway Soil is grains, water and air, and the part the grains carry, the effective stress, is what gives it strength. Change the water and you change the ground: rain weakens a slope, a flood weakens a footing, an earthquake can take a sand's strength away for half a minute, and a soft layer can amplify the shaking above it. Clay is slow, so what you build today is still settling in forty years. Try this: next time you walk across wet sand, stamp and watch the water rise around your foot. That is effective stress changing under your feet. Then notice how the shoreline slope relates to the angle the sand will stand at. Sources & Further Reading Sources. The three-phase relations, Terzaghi's effective stress principle, one-dimensional consolidation, the infinite-slope model, the strip-footing bearing capacity equation and the quarter-wavelength resonance of a soil layer are standard textbook geotechnics (see the books above). Every number in the interactives was computed in node: the void-ratio example (e = 0.7: 19.3 and 15.3 kN/m³), effective stress at 10 m (196 − 78.5 = 117.5 kPa), the consolidation example (sc = 716 mm; t90 = 42.4 yr one-way, 10.6 yr two-way; U checked against Tv = 0.197 and 0.848), the slope example (tan 32° ÷ tan 25° = 1.34 dry, 0.65 wet), bearing capacity (Nq = 18.40, Nc = 30.14 and Vesic's Nγ = 22.40 at 30°; qult about 694 kPa dry and 293 kPa flooded) and the soil period (4 × 37 ÷ 74 = 2.0 s). All soil properties, the 5%-damped building and the 0.1 s-per-storey rule are illustrative simplifications. Facts about Pisa (construction dates, 5.5 and 3.97 degrees, 1990 and 2001, 870 tonnes of lead, 70 tonnes of soil), Christchurch (22 February 2011, Mw 6.2, 6.7 km, 185 deaths, red zone) and Mexico City (19 September 1985, Mw 8.0, over 350 km, 412 collapsed, 7–37 m of lakebed clay, 6–15 storeys) are taken from the Wikipedia articles on those events; the death-toll range for Mexico City varies widely by source. Go deeper. Das, Budhu, Terzaghi–Peck–Mesri and Petroski are standard titles; check the current edition with your bookshop. The NPTEL course listing was confirmed on the NPTEL site. 🌐 NPTEL Soil Mechanics/Geotechnical Engineering I (IIT Kharagpur) 📘 Karl Terzaghi, Ralph Peck and Gholamreza Mesri, Soil Mechanics in Engineering Practice 📘 Henry Petroski, To Engineer Is Human Step 15. Hydraulics & Water Resources: Moving Water Where People Need It ------------------------------------------------------------------------ Why: Rivers, floods, pipes and dams: moving water where people need it, and away from where they don't. Pipes, pumps, rivers, floods, dams and levees: the part of civil engineering that decides where water goes, and what happens on the day it goes somewhere else. Turn on a tap on the tenth floor and water comes out. Somebody had to lift it there, and somebody had to decide how wide the pipe should be so it would not come out as a trickle. Now look at a river after a storm. It was ankle-deep last week and now it is brown, fast and over the car park. Same water, same laws, a different job: in the pipe, the engineer is delivering it; in the river, the engineer is surviving it. This is the course where fluid mechanics goes to work on whole landscapes. The graduate does not solve Navier–Stokes. They keep a short list of habits: count water in and water out, think in "head" (energy as a height), ask how often the bad day comes, and remember that every barrier has a way to fail. The course in one box Who takes it: civil and environmental engineers (and agricultural and some mechanical engineers), usually in year 3, as one course in hydraulics and often a second in hydrology. Civil programmes shaped by the ASCE Body of Knowledge all include it; geoscience and public-policy programmes teach the hydrology half. What it's for: sizing the pipes, pumps, channels, reservoirs and flood defences that move and store water, and judging the risk that they are overwhelmed. Comes after: Fluid Mechanics (34/15), Probability & Statistics (34/05). Leads to: Environmental Engineering (34/68) for the water once it is dirty, Geotechnical Engineering (34/66) for dams and levees as soil structures, and Corporate Finance (34/51), the engineering-economics course, for deciding what protection is worth paying for. Honest line: you won't be able to design a dam or a drainage network after this. You will see why narrowing a pipe is so costly, why a "100-year flood" can strike twice in a decade, and what went wrong at New Orleans, Teton and Vajont. This piece picks up where Fluid Mechanics (34/15) stops. Pressure, Bernoulli and the Reynolds number are assumed; if they are rusty, read 34/15 first. Pipes and pumps: where two curves cross Engineers measure the energy in water as a height: head, in metres. Bernoulli's three terms (pressure, speed, elevation) each become a height, and they add up to the total head. In a real pipe the total slowly drops along the way, because friction turns it into heat. That drop is the head loss. hf = f · (L / D) · v² / 2ghead lost = friction factor × (length ÷ diameter) × speed-head This is the Darcy–Weisbach relation. The friction factor f comes from the Reynolds number and the roughness of the pipe wall (the Moody chart, or a formula such as Swamee–Jain). The part to keep in your head is the diameter. For a fixed flow, the speed goes as 1/D², so v² goes as 1/D⁴, and the extra 1/D makes it about 1/D⁵. Halve the diameter and the loss is roughly 32 times bigger. In the model below, 30 litres a second through 1,500 m of 200 mm pipe loses about 7 m of head; through 100 mm it would lose about 233 m, which no pump would supply. Pipe diameter is the cheapest dial in the whole system. A pump adds head, but less of it the more water you ask for. Plot the pump's head against flow and the line falls. Plot what the pipe needs (the lift in height, plus the friction loss) and the line rises. The pump runs where the two curves cross. Everything else in a pumping design is moving one curve to put that crossing where you want it. Whole networks follow the same logic as an electric circuit: head is voltage, flow is current, and at every junction the flows add up to zero (see Circuit Analysis, 34/34). The classic hand method for loops of pipe, by Hardy Cross (1936), is just that. Interactive · Find the operating point Pipe diameter Lift (static head) A made-up pump (50 m of head at zero flow, 35 m at 60 L/s) feeding 1,500 m of pipe, roughness 0.1 mm, plus 5 velocity heads of fittings. Water at 20 °C. Try 250 mm and 15 m, then shrink the pipe to 150 mm and watch the flow fall. Open channels: when water outruns its own waves In a pipe, pressure pushes the water. In a river, a canal or a storm drain flowing part-full, only gravity does, and the water has a free surface that can wave. That surface decides everything. A small ripple on shallow water travels at √(g·y), where y is the depth. Compare the water's speed with the ripple speed and you get the Froude number. Fr = v / √(g · y)flow speed ÷ speed of a surface wave Below 1 the flow is subcritical: calm, deep, and ripples can travel upstream, so a downstream gate or bridge pier "tells" the water far behind it. Above 1 the flow is supercritical: fast, shallow, and ripples are swept away downstream, so nothing downstream can warn it. For uniform flow in a long channel the working rule is Manning's equation (from the Irish engineer Robert Manning, 1890s): speed grows with the slope's square root and with the channel's "hydraulic radius" to the power 2/3, and falls with roughness. When fast shallow water meets slow deep water, it cannot change gently. It jumps: a sudden rise in surface, a roaring, churning roller, and a large share of the energy burned off as turbulence. The depths on the two sides are linked by one formula, and the energy lost is a large fraction of what the flow carried. That is useful: below a spillway, a concrete stilling basin forces a jump so the water leaves the dam slowly instead of scouring the river bed. You can see the small version in your kitchen: the circular ring around where a tap's stream hits the sink is a hydraulic jump. Interactive · The hydraulic jump Depth before the jump Speed before the jump Rectangular channel, horizontal bed, no friction on the bed. Sequent depth y₂ = (y₁/2)(√(1 + 8Fr₁²) − 1); energy lost = (y₂ − y₁)³ / (4y₁y₂). The jump is drawn to a made-up scale. A catchment: rain in, flood out Hydrology starts with an accounting identity: rain that falls on a catchment either soaks in, evaporates, or runs off. Only the runoff makes the river rise, and how fast it arrives decides how high. Dry, bushy ground soaks up the first part of a storm. A roof, a road and a car park do not. A standard classroom model (the Curve Number method from the US Soil Conservation Service, now NRCS, with a triangular "unit hydrograph" for timing) captures this in a few lines. Land cover sets a curve number from about 30 (very absorbent) to 100 (a paved surface). The more paved, the bigger the runoff for the same rain, and the water reaches the river sooner, because pipes and gutters carry it faster than meadow does. Both effects stack on the peak. This is the reason developers are now asked for detention ponds and permeable paving: to give back some of the storage the pavement took. Interactive · Pave the catchment Share of land paved Storm: rain in 3 hours Illustrative 10 km² catchment: rural land CN 70, paved CN 98, area-weighted. Lag falls from 3 h (unpaved) to 1.2 h (fully paved). 80 mm in 3 h is a very heavy storm. Dashed line: the same storm on unpaved land. Return periods: why a 100-year flood isn't once a century Engineers cannot say "it will flood in 2031". They say how likely a given flood is in any one year. A "100-year flood" is the flow that has a 1-in-100 chance, or 1%, of being equalled or beaten in any year. That is all it means. It is a probability per year, not a schedule, and not a promise that 99 quiet years must follow. P(at least one in n years) = 1 − (1 − 1/T)nT = return period in years Over a 30-year mortgage the chance of meeting a 100-year flood at least once is 1 − 0.99³⁰ = 26%. Over a 100-year life it is 63%. That also means there is a 37% chance of none in a century, and a 26% chance of two or more. The independence is the key: the river does not remember last year's flood. Floods in successive years are rare but entirely possible. Two cautions belong in every graduate's head. First, the 100-year flow is estimated, usually by fitting a distribution to a record of perhaps 50 to 100 years, then extrapolating. The tail is uncertain. Second, the method assumes the climate and the catchment are the same as in the record. Paving, deforestation and a warming climate all break that assumption, which is why many agencies now say "1% annual exceedance probability" and revise their maps. Interactive · The flood odds calculator Return period T Years you are exposed Roll one random period Independent years, constant chance. The strip shows one random run with a fixed starting seed, then a new one each press. Red cell: a flood at least this big. Dams and reservoirs: moving water through time Rivers flood in spring and run dry in late summer, while people need water all year. A dam is a machine for shifting water through time: store the surplus, release it when it is short. Dams also make electricity, hold back floods and raise a water level for navigation or irrigation. Their costs are real too: sediment slowly fills reservoirs, evaporation loses water, and a valley, its people and its fish are changed for ever. The engineering of storage is bookkeeping, a mass balance for the whole reservoir: storage today = storage yesterday + inflow − release − losses. The surprise is the limit. A reservoir cannot supply more than the river brings in on average, however big. Below that limit, bigger storage buys reliability, but with sharply diminishing returns. Play with it below. Interactive · Size the reservoir Storage capacity Yearly demand One hundred invented years of river flow (mean 100 units, varying by about ±40% and with some dry runs), not a real river. No evaporation, no sediment, no environmental release. When barriers fail: New Orleans, Teton and Vajont A levee is a wall of earth or concrete that keeps a river or the sea off land it would otherwise cover. It does one job: hold a water level up to a design height. It moves risk rather than removing it: land behind it develops, and the day the design is exceeded or the structure fails the losses are larger. The 100-year idea matters here: a levee built to the 100-year level is overtopped, on average, once in a hundred years. Hurricane Katrina struck on 29 August 2005. The storm surge overtopped some levees and floodwalls, notably on the eastern side of the city, and that part was a question of design height. The disaster came from where defences failed below their design level. The 17th Street and London Avenue canals and part of the Industrial Canal broke while water was below the top of the walls. By 31 August about 80% of the city was flooded, some areas under up to 15 feet of water. Published death counts for Louisiana differ, roughly 1,000 to 1,500. The US Army Corps of Engineers, which had designed and built the system, acknowledged in a June 2006 report faulty design, incomplete sections and substandard construction. A 2007 American Society of Civil Engineers review called it the worst engineering catastrophe in US history, and concluded about two-thirds of the flooding could have been avoided if the levees had held. The lesson is the one every graduate keeps: a defence is a system of links, and what matters is the weakest one (here soil under the walls and the walls' foundations; see Geotechnical Engineering, 34/66). It also shows how a "protected" city can be misled into thinking it is safe. Two dam failures make the same point from other angles, and both are often told wrong. Teton Dam, Idaho, 5 June 1976. It was a 93 m earth-fill dam, finished in November 1975, failing during its first filling. Seepage found a path through or beneath the dam and began carrying soil away: internal erosion, called piping. The leak grew from a trickle to a breach within hours; the reservoir emptied and 11 people died. An independent panel named internal erosion the most probable cause, but could not say whether it began through fissured rock under the grout seal or through cracks in the dam's core. Neither a flood nor an earthquake overwhelmed it. A dam can be strong and still lose to the water seeping through and under it. Water in soil, and what it can carry off, is the subject of Geotechnical Engineering (34/66). Vajont, Italy, 9 October 1963. The concrete arch dam, at 262 m among the tallest in the world, did not fail. It stood. About 260 million cubic metres of mountainside slid into the reservoir, the displaced water went over the crest, and the flood killed about 2,000 people in the valley below (the official count is 1,917; estimates run up to roughly 2,500). The warning signs on the slope had been disregarded, and a handful of engineers were convicted in a trial moved to L'Aquila. Experts still debate how much rain, reservoir level and earthquakes each contributed. The lesson: the engineer designs the dam, but the dam sits in a landscape, and the landscape must be understood too. Water scarcity: the same equations, worldwide On the planet as a whole, water is not scarce; fresh water is badly placed and timed. The hydrology graduate thinks in balances, region by region and season by season. Where demand approaches the renewable supply, extra demand is met by drawing down stored water: aquifers, glaciers, reservoirs. Irrigated agriculture is the largest user worldwide, by most estimates around two thirds to 70% of withdrawals, so most of the water problem is a food problem. The Joint Monitoring Programme of WHO and UNICEF reported in August 2025 that, in 2024, about 2.1 billion people (roughly one in four) still lacked safely managed drinking water, including 106 million who drank directly from untreated surface water. Those figures describe services, not rainfall. Scarcity is as often a pipe, a pump, a price and a government as it is a drought. Engineers also disagree about the cure: big dams and transfers versus recycling, efficiency, demand management and small local systems. The right mix depends on the place, and civil engineers inherit the argument. What a graduate just sees Head is energy as a height; pipes lose it to friction, and the loss is roughly 1/D⁵: a skinny pipe is expensive to push through. A pump works where its falling curve crosses the pipe's rising one. Change either curve and the flow moves. Fast shallow water outruns its own ripples; where it meets slow deep water it jumps, and the jump is a free energy-burner. Pavement raises the runoff and shortens the delay, so the flood peak grows twice over. A 100-year flood is a 1% chance each year. Over a mortgage, one in four. A reservoir can never yield more than the average inflow, and each extra bit of storage buys less. Defences fail at their weakest link, often below their design load: soil, seepage, foundations. Failures happen to dams and levees as parts of a landscape: Teton leaked, Vajont slid, New Orleans' walls gave way. Where it shows up next Environmental Engineering (34/68): the same pipes and channels now carry sewage and treated water; the river's flow sets how much pollution it can take. Geotechnical Engineering (34/66): seepage, slope stability and foundations are what decided Teton, Vajont and New Orleans. Probability & Statistics (34/05) and Engineering Economics (34/51): choosing the return period is a decision about money and lives. The real world: city drainage, irrigation schemes, hydropower, flood maps and insurance, port and canal design, and the water utility in every town. If you want to go deeper 📘 Textbook: Larry W. Mays, Water Resources Engineering (Wiley): the broad standard, covering hydraulics and hydrology together. 📘 Textbook: David A. Chin, Water-Resources Engineering (Pearson): strong on worked design problems. 📘 Classic for the river half: Ven Te Chow, Open-Channel Hydraulics (1959), and Chow, Maidment and Mays, Applied Hydrology (McGraw-Hill). 📘 Gentler: John McPhee, The Control of Nature (1989): essays, including a long one on the Army Corps' fight to hold the Mississippi in its channel. ▶️ Free course: MIT OpenCourseWare 1.72 Groundwater Hydrology (Fall 2005), a graduate course for the water-underground side; for the surface side, the fluid mechanics course linked in 34/15 is the base. 🌐 Easier first read in this library: Fluid Mechanics (34/15). The takeaway Water engineering is accounting plus risk. In a pipe, head is lost to friction, the pump sits where two curves cross, and diameter is the strongest lever. In a river, speed against wave speed decides whether the flow is calm or jumps. A catchment turns rain into a flood, and paving makes that flood bigger and quicker. A 100-year flood is a 1% chance a year, not a schedule. Dams shift water in time, can never beat the average inflow, and fail through seepage, slopes and weak links in the chain, as at Teton, Vajont and New Orleans. Try this: after the next heavy rain, find where the water on your street goes and what it runs over. Then look up the flood map for your town and see what its "1%" line is. Is your house inside it? Sources & Further Reading Sources. Darcy–Weisbach, Swamee–Jain, Froude number, Manning's equation, the hydraulic jump relations, the SCS Curve Number method and the return-period formula are standard textbook material (Mays; Chin; Chow; NRCS National Engineering Handbook Part 630). All interactive numbers were computed in node. Teton Dam: US Bureau of Reclamation fact sheet and the 1976 Independent Panel report (hosted at damfailures.org): 11 deaths, internal erosion the most probable cause, exact initiation not determined; construction and failure dates from the Teton Dam article on Wikipedia. Vajont: Wikipedia article (260 million m³, 262 m dam, official toll 1,917, range up to about 2,500, trial at L'Aquila). New Orleans 2005: Wikipedia, Effects of Hurricane Katrina in New Orleans (80% flooded by 31 August, failures below design height, June 2006 Corps acknowledgement, ASCE review). Water access: WHO/UNICEF Joint Monitoring Programme, report of August 2025 (2.1 billion lack safely managed drinking water in 2024; 106 million use untreated surface water). Hardy Cross (1936), the Dutch 1953 flood and the irrigation share of withdrawals are cited from standard texts and FAO figures, not re-checked for this piece. 🌐 WHO/UNICEF JMP: household water, sanitation and hygiene, 2025 update 🌐 Dam Failures: Teton Dam (Idaho, 1976) 🌐 MIT OpenCourseWare 1.72 Groundwater Hydrology 📘 John McPhee, The Control of Nature Step 16. Vibrations: Everything Is a Spring and a Mass ------------------------------------------------------ Why: Earthquakes and wind shake structures; resonance is the danger. Why a guitar string has a pitch, a washing machine walks, and a bridge can wobble itself apart: one equation, and the number that decides how badly things shake. Pluck a guitar string and it rings at one pitch, no matter how hard you pluck it. Push a child on a swing in time with the swing and a gentle shove builds a big arc. Drive over a washboard road at the wrong speed and the whole car shudders. Different things, one idea: almost every object has a frequency it prefers to move at, and it responds wildly to being pushed at that rhythm. A vibrations graduate sees this everywhere. A phone buzzing on a table, a washing machine walking across the floor, a skyscraper swaying in a typhoon, a tall bridge that wobbles under a crowd: each is a mass, a spring and some friction, and each has a number that says how badly it will misbehave. The course in one box Who takes it: mechanical, aerospace and civil engineers, usually in year 3. What it's for: predicting how a structure or machine moves when it is pushed, and how to stop that motion from becoming a failure. Comes after: Dynamics (34/12) and Differential Equations (34/03). Leads to: Control Systems (34/20), which uses the very same equation, plus Machine Design (34/19) for shafts, gears and bearings, and Propulsion (34/21) for blades that must not flutter. Fluid Mechanics (34/15) joins it whenever wind or water drives the motion. Honest line: you won't be able to design a damper for a tower after this. You will see why a bridge can wobble, why soft mounts protect a machine, and what Taipei 101's giant golden ball is for. Everything is a spring and a mass Take anything that springs back when displaced: a diving board, a building, a guitar string, a car on its tyres. Close up, it behaves like a mass on a spring. Newton's law says that the spring pushes back in proportion to how far you stretch it, so the acceleration is proportional to the displacement and points the other way. That gives a motion that repeats, at a frequency that depends on only two things. ωn = √(k / m)natural frequency (rad/s) = square root of stiffness ÷ mass Stiffer means faster; heavier means slower. Divide by 2π to get hertz. Worked example: a 1 kg mass on a 100 N/m spring has ω = 10 rad/s, so f = 1.59 Hz, one swing every 0.63 s. A car body of about 400 kg per wheel on a 20,000 N/m spring comes to 1.1 Hz, roughly a walking pace, which is why a car feels like it is floating. Four times the stiffness, or a quarter of the mass, doubles the frequency, because of the square root. The size of the pluck doesn't appear in the formula. The frequency is a property of the object itself, which is what makes it so useful. Interactive · Mass and spring 1 kg, 100 N/m4 kg, 400 N/m16 kg, 400 N/m Mass Spring stiffness Pause ω = √(k/m) = 10 rad/s → 1.59 Hz, one swing every 0.63 s. Ideal, undamped, computed live. Note that 4 kg on 400 N/m has the same frequency as 1 kg on 100 N/m: only the ratio k/m counts. Quadruple the mass and the frequency halves. Damping: where the energy goes Real things stop ringing because friction, air and bending metal turn motion into heat. That is damping. Engineers measure it with a single dimensionless number, the damping ratio ζ (zeta): the actual damping divided by the critical damping, the exact amount that just prevents any overshoot. ζ = c / (2√(k·m))damping ratio = damper strength ÷ the critical value ζ below 1 (underdamped): it overshoots and rings, each swing smaller. A guitar string has tiny ζ and rings for seconds; a car's shock absorbers are set well below 1 so the wheels still follow the road. ζ = 1 (critical): the quickest return to rest with no overshoot. A good door closer aims here. ζ above 1 (overdamped): no overshoot, but a slow, syrupy creep. For large ζ it behaves like a first-order system, a single exponential with a time constant, which is exactly what the next subject uses to describe cruise control. For ζ below 1 each swing is smaller than the last by a fixed factor. At ζ = 0.05 each swing is 73% of the one before, so it rings for dozens of cycles. At ζ = 0.2 it is 28%: three swings and it's gone. Interactive · Let go and watch 0.020.20.71 (critical)2 (overdamped) Damping ratio Underdamped: each swing is 28% of the last. Exact solution of m·x″ + c·x′ + k·x = 0, released from rest, checked against a step-by-step numerical solution. Forced vibration and resonance Now push the mass back and forth at a steady frequency. At first the mass just follows the push. If you push very fast, it barely moves, because there isn't time. But push at the natural frequency and each push arrives just as the mass is moving the right way, so the energy adds up cycle after cycle until damping drains energy as fast as you put it in. That is resonance. The amplitude, as a multiple of what a steady push of the same size would do, is the famous curve: amplitude ÷ static = 1 / √[(1 − r²)² + (2ζr)²]r = driving frequency ÷ natural frequency At r = 1 the first bracket vanishes and the amplification is 1/(2ζ). A lightly damped structure with ζ = 0.02 moves 25 times further than a static push; with ζ = 0.5 barely more than a static push. So damping, not the size of the push, decides how bad resonance is. The motion also shifts from moving with the push (below resonance) to moving against it (above), with a 90° lag exactly at resonance. Interactive · The resonance curve Slow pushAt resonanceFast push Damping ratio ζ Driving frequency At resonance with ζ = 0.1 the motion is 5× a static push. Steady-state response of a damped mass-spring to a sinusoidal force. For ζ = 0.05 the peak is 10×; for ζ = 0.5 it shrinks to 1.15× and moves to lower frequency. Isolation: soft mounts and the √2 rule Sometimes the push is something you want to keep away from a structure: an engine shaking a car body, a washing machine's drum shaking a floor, a seismometer ignoring the ground. The trick is to put a spring between them. The fraction of shaking that gets through is called transmissibility, and it has a surprising shape: below r = √2 (about 1.41) the mount makes things worse, and only above it does isolation begin. Above √2, more damping lets more through at high frequency. So the rule: make the mount's natural frequency well below the frequency you want to block. To cut a 30 Hz engine's force to 10%, with light damping you need r ≈ 3.3, so a mount with a natural frequency near 9 Hz. In practice that means a soft spring that sags a few millimetres under the load. The price is that every time the machine starts or stops it sweeps through the mount's natural frequency, so you want some damping to survive the crossing. Interactive · An engine on mounts Mount natural frequency Mount damping ratio A 9 Hz mount under a 30 Hz engine lets about 13% of the force through. The curve is the standard transmissibility of a damped mass-spring. Try moving the mount to 25 Hz and watch it amplify. Modes: a structure is a choir of springs A single mass has one natural frequency. A string, beam, plate or aircraft wing is a continuous chain of masses and springs, so it has many: a mode is one pattern of motion in which every part moves in step at one shared frequency. The lowest mode is the whole string bowing out. The second has a still point in the middle (a node), the third has two, and so on. Pluck the string and you excite several at once; that blend is the timbre of the instrument. For an ideal string the modes are whole-number multiples of the lowest, a harmonic series. A guitar's low E at 82.41 Hz also rings at 164.8, 247.2 and so on. A beam fixed at one end (a diving board, a wing, a skyscraper) is different: its second mode is 6.27 times the first and its third 17.5 times, so the overtones are not neat multiples. Engineers find these modes first, then keep forces away from them. Interactive · Mode shapes Guitar stringCantilever beam Mode 1Mode 2Mode 3Mode 4Mode 5Mode 6 Pause Mode 2 of a string tuned to 82.41 Hz: 164.8 Hz, with 1 still point (node). Exact mode shapes (sine for the string; the standard cosh/cos formula for the beam, using βL = 1.875, 4.694, 7.855, 10.996, 14.137, 17.279). Shown in slow motion, with the amplitude exaggerated. Orange dots are nodes. When it goes wrong: Tacoma Narrows and the wobbly bridge Tacoma Narrows (1940). The usual textbook story says the wind shook the bridge at its natural frequency and resonance tore it down. That version is wrong. The bridge opened on 1 July 1940 and fell on 7 November, in a wind of roughly 65 km/h (reports range from 35 to 42 mph). The deck was slender and flexible in twisting, and on the day it went from bouncing to twisting about once every four seconds. In 1991 K. Y. Billah and R. H. Scanlan, in the American Journal of Physics, argued that the cause was aeroelastic flutter: the twisting itself changed the wind's forces on the deck so that they fed energy into the twist. This is a self-excited motion, with no steady external beat to resonate with. Once the wind passed a threshold speed, the aerodynamic effect acted like negative damping. Some details are still argued over, but "simple resonance" is not the accepted account. London Millennium Bridge (2000). The footbridge opened on 10 June 2000 and was closed two days later after it swayed sideways under crowds. Walking pushes sideways on the deck about once per two steps. If the deck sways, people adjust their stride, and on this bridge their combined pushes ended up in phase with the sway, adding energy instead of removing it. The effect only began once the crowd passed a critical size. The Arup engineers' study showed that any bridge with a lateral frequency below about 1.3 Hz could suffer the same if enough people were on it. The fix was not to stiffen it but to add dampers (viscous and tuned mass dampers) and it reopened in early 2002. Engineers still debate how much of the effect is people falling into step with each other and how much is each person balancing against the sway. Both stories share one mechanism: something that pushes in a way that reduces the effective damping. Once the damping reaches zero, any disturbance grows. Try it. Interactive · Crowd, damping and the tipping point As openedWith added dampers People on the deck With 100 people the sway still dies out. The tipping point is about 157 people. An illustrative model, not the real bridge's data: 100 tonnes of moving mass, 0.5 Hz, 0.5% natural damping, each person adding 20 N·s/m of negative damping, and +5% damping for the "added dampers". It shows the mechanism: the threshold is where the crowd's push equals the structure's own damping. Tuned mass dampers: a second spring that eats the vibration You can fight resonance by adding damping, or by adding a second mass-spring tuned to the first. A tuned mass damper is a small mass on a spring and damper, attached to the structure and tuned to its natural frequency. When the structure resonates, the small mass swings against it and absorbs energy. The best tuning, worked out by J. P. Den Hartog and still the starting point, is simple: for a mass ratio μ, tune the damper to 1/(1 + μ) of the building's frequency, with damping ratio √(3μ/(8(1+μ)³)). Taipei 101 is the famous case: a 660-tonne steel sphere about 5.5 m across (728 short tons) hangs like a pendulum in the upper floors, and in Typhoon Soudelor in 2015 it swung about a metre. The makers claim it can reduce sway by up to 40% (a company figure). In the simulator below, a damper that weighs only 2% of the structure cuts the peak response from 50× to about 9×. Interactive · Tune the damper Damper mass Damper tuning (× building frequency) Damper damping ratio Optimal settingsMistuned (15% off) Peak response without the damper: 50×. With it: about 8.7×. Two-mass model computed with complex arithmetic; the structure has 1% damping. The damper's mass is the slider's value in tenths of a percent, so 20 means 2%. At 100 (10%) you can see the response split into two lower humps. What a graduate just sees Everything with mass and springiness has a natural frequency, √(k/m), set by the object and not by how hard it is hit. Stiffer goes up, heavier goes down, and by the square root: four times the stiffness doubles the frequency. Damping ratio, not the force, decides how bad resonance is: the peak is about 1/(2ζ). Below √2 times its natural frequency a soft mount amplifies; above it, it protects. Any structure has a whole family of modes; engineers keep forces away from the low ones. Tacoma Narrows was flutter, a self-excited instability, not simple resonance. When something fails by growing vibration, look for what has made the damping effectively negative. If you can't remove the push, add damping (or a tuned damper) rather than making it stiffer. Where it shows up next Control Systems (34/20): the equation of a mass with damping and a spring is the standard second-order system, with the same ζ and ωn. Feedback is how you add damping by software: a flight controller does for a drone what a viscous damper does for a bridge. Machine Design (34/19): shafts and gearboxes have critical speeds at which they whip; turbines and rotors are balanced so their push stays small. Propulsion (34/21) and Fluid Mechanics (34/15): compressor and turbine blades and wings can flutter, the Tacoma mechanism at jet speeds. The real world: earthquake-resistant towers, car and aircraft ride, hard-drive heads, noise and comfort engineering, musical instrument design. Computer Science: analysing vibrations is Fourier analysis (see Fourier Transform, 12/12), and software simulating springs underlies physics engines, animation and robotics. If you want to go deeper 📘 Textbook: Singiresu S. Rao, Mechanical Vibrations (Pearson): the most widely assigned undergraduate text, with many worked examples. 📘 Textbook: Daniel J. Inman, Engineering Vibration (Pearson): clearer on modelling and on MATLAB use. J. P. Den Hartog's Mechanical Vibrations (Dover) is the short classic that gave us the tuned damper. 📘 Gentler: The Feynman Lectures on Physics, Vol. I, ch. 21 "The Harmonic Oscillator" and ch. 23 "Resonance", free online. ▶️ Free course: MIT OpenCourseWare 2.003SC Engineering Dynamics (Fall 2011), an introduction to vibrations and dynamics of lumped mechanical models, with recitations on frequency response. 🌐 Easier first read in this library: Waves and Sound (13/13) and Fourier Transform (12/12). The takeaway Everything springy is a mass and a spring with natural frequency √(k/m). Damping (ζ) decides whether it rings, creeps or returns cleanly, and sets the peak of resonance at about 1/(2ζ). A soft mount protects only above √2 times its natural frequency. Structures have many modes, and the dramatic failures, such as Tacoma Narrows and the Millennium Bridge, came from forces that made the effective damping negative, not from a simple beat. A tuned mass damper is a second spring placed where it can steal the energy. Try this: hold a ruler overhanging a desk and twang it. Shorten the overhang and hear the pitch rise. Then stick a lump of putty on the tip and hear it drop: stiffness up, frequency up; mass up, frequency down. Sources & Further Reading Sources. Single-degree-of-freedom theory, modes, isolation and the Den Hartog absorber are standard textbook material (Rao; Inman; Den Hartog). Tacoma Narrows: K. Y. Billah and R. H. Scanlan, "Resonance, Tacoma Narrows bridge failure, and undergraduate physics textbooks", American Journal of Physics 59 (1991) 118–124; MIT 18.03 notes (Appendix B); dates and wind reports vary by source (35–42 mph). Millennium Bridge: P. Dallard et al., The Structural Engineer 79(22), 2001; Ingenia, "Stabilising the London Millennium Bridge"; Wikipedia, "Synchronous lateral excitation". Taipei 101: Britannica and the Guinness record page (660 metric tons; floors given as 87th to 92nd by different sources; the 40% is a manufacturer's claim). Computed in node: all frequencies, ratios, peaks and the damper optimum (checked against a numerical scan of the two-mass model); the release response was checked against a step-by-step numerical solution. 🌐 MIT OpenCourseWare 2.003SC Engineering Dynamics 🌐 Billah and Scanlan (1991), scanned paper and Ingenia: Stabilising the London Millennium Bridge 📘 Rao, Mechanical Vibrations · Inman, Engineering Vibration · Den Hartog, Mechanical Vibrations Step 17. Structural Design: Concrete & Steel That Bend Before They Break ------------------------------------------------------------------------ Why: Turning forces into real beams, columns and slabs, under the codes that keep them safe. Design is deciding how rare a failure must be, then putting each material where it works best, so that a structure that fails does so slowly and loudly. A beam that held its load on test day has told you very little. What matters is whether it will hold on the worst day of fifty years: the day a crowd stands in the wrong place, the concrete in that corner was a bit weak, and a storm is blowing. Engineers cannot know that day in advance. What they can do is decide how unlikely they will let it be. Structural design is the art of doing that with steel and concrete. Structural analysis (34/64) hands you the forces. Design turns them into a real thing: how deep the beam, how many bars of reinforcement and where, how many bolts, and how to arrange it so that if the structure is ever overloaded it gives warning rather than letting go. It is a small number of ideas applied with great care. The course in one box Who takes it: civil and structural engineers, usually in the third and fourth years, as two courses: reinforced concrete design and steel design. Architects meet a lighter version. Programmes that follow the ASCE Civil Engineering Body of Knowledge require design in at least one structural material.What it is for: choosing the size and detail of beams, columns and slabs so that they safely carry the forces from structural analysis, under the rules of a design code.Comes after: Structural Analysis (34/64) and Materials Science & Manufacturing (34/13).Leads to: Construction Management (34/70), which has to build what is drawn, and the capstone Systems Engineering (34/22). Foundations in this design rest on Geotechnical Engineering (34/66).What this piece will and won't give you: you won't be able to design a beam: that needs a code, and every value here is an example, not a design value. You will see why the steel in a concrete beam goes where it does, why codes use safety factors, and why a structure that bends first is safer. Easier first reads: Mechanics of Materials (34/11), Materials Science & Manufacturing (34/13) and Statics (34/10). The Hyatt Regency walkway collapse is told in Engineering Drawing & CAD (34/18) and not repeated here. Limit states: design for the bad day, with factors drawn from statistics Modern codes design for limit states: conditions beyond which a structure no longer does its job. The ultimate limit state is collapse or loss of stability. The serviceability limit state is the quieter failure: a floor too bouncy to live on, a crack too wide, a beam that sags until the doors jam. For the ultimate limit state, the rule in words is: the design strength of the member must be at least as big as the factored loads. φ Rn ≥ Σ γi Qi Read it as: take the nominal strength, shrink it by a factor φ below 1, and make sure it still beats the loads, each stretched by a factor γ above 1. The factors are not superstition. Loads vary (a live load is far more uncertain than a dead load, which is why it gets the bigger factor) and strength varies (one batch of concrete, one steel heat). Treat both as bell curves. The structure fails where the load curve overlaps the strength curve, and the factors are chosen so the overlap is acceptably small. Eurocode EN 1990 puts the target for an ordinary building at a reliability index β of 3.8 over 50 years, which is a probability of failure of about 7 in 100,000. Notice what the overlap depends on. A member with a mean strength 1.6 times the mean load, with 12% scatter in strength and 20% in load, has β of only 2.2 and a failure chance of about 1.5%. Raise the ratio to 2.2 and β reaches 3.6, about 1 in 7,000. Calibrating φ and γ so that ordinary members land near the target is what code committees spend years on. Older "allowable stress" design instead divided strength by one overall factor and did not know how safe the result was. Interactive · Why safety factors exist Move the average strength above the average load, and change how much each one scatters. The overlap is the chance of failure. Mean strength ÷ mean loadScatter in strengthScatter in load Reinforced concrete: put the steel where the concrete would crack Concrete is very good at being squeezed and poor at being stretched: its tensile strength is only about a tenth of its compressive strength. Steel is good at both. Reinforced concrete puts each where it is needed. Load a simply supported beam and it sags: the top squeezes and the bottom stretches, and between them is a line, the neutral axis, with no stress at all. So the steel goes at the bottom, as far from the neutral axis as the cover allows, where its pull works with the longest lever arm. The two materials expand with temperature at nearly the same rate, and ridges on the bars grip the concrete, which is why the combination works at all. In a cantilever the top is the stretched side, so the bars go in the top. Mn ≈ As fy (d − a/2) Read it as: the pull in the steel times the lever arm between the pull and the push in the concrete. Four 20 mm bars (1,257 mm²) of 500 MPa steel, 450 mm below the top of a 300 by 500 mm beam of 30 MPa concrete, make a push zone 82 mm deep and a capacity of 257 kN·m. It also sets how the beam will fail. If the steel is modest, it yields first and the beam bends a long way with wide cracks, which is a warning. If there is too much steel, the concrete crushes first, suddenly. Codes force the first kind by limiting the steel and by cutting φ for the second. Interactive · Where do you put the rebar? Slide the bars from the bottom to the top of the beam and watch the capacity collapse. Add bars until the failure mode changes. Number of 20 mm barsBar depth below the topFactored load on 6 m span Prestressing: squeeze the beam before you load it Cracks let in water and corrosion, and a cracked beam is less stiff. Prestressing avoids them. The engineer stretches steel strands, bolts them against the ends of the beam, and lets them squeeze the concrete before any load arrives. When the load comes and tries to stretch the bottom, the squeeze has to be undone first, and the concrete never goes into tension. The strands are placed low, so that the squeeze itself bends the beam upward and fights the sag. σ = −P/A ∓ P e c/I ± M c/I Read it as: the even squeeze, plus the bending the squeeze causes because it is off-centre, plus the bending from the load. For a 300 by 600 mm section over 10 m, 20 kN/m of load would put 13.9 MPa of tension on the bottom, four times what concrete can take, so it would crack. With 1,000 kN of prestress 200 mm below the centre, the bottom is in 2.8 MPa of compression. The prestressing also lets a floor be thinner for the same span. Eugène Freyssinet was the pioneer of the method in the 1920s and 1930s. The thing to watch is the unloaded stage: the beam may crack on top, because the squeeze pushes it up. Interactive · Squeeze the beam in advance Without prestress, the beam cracks. Add force and lift the strands away from the centre. Prestress forceStrand height below centreTotal load Steel: strong, slim, and only as good as its connections Steel is strong, stiff and made in a factory to tight tolerances, so steel members are slim: I-beams, tubes and angles with a high stiffness for their weight (34/11). Slimness is its weakness. A slim column buckles; a slim flange wrinkles (local buckling); a deep beam can twist sideways (lateral-torsional buckling). Much of steel design is deciding how to prevent these, by bracing, by choosing compact sections, and by checking each failure mode with its own formula. The real work is at the connections. Members are joined with bolts or welds, and the joint must be strong enough in every way it could fail: the bolts shearing, the plate crushing against the bolt, the end tearing out, the plate rupturing through the bolt holes. The idea underneath is a chain. A chain is as strong as its weakest link, so the designer chooses which link is the weakest, and makes it the one that stretches. A tension plate with holes shows it: the whole plate yields at 0.9 × Ag × fy, but through the holes only the net area carries load, and failure there is rupture, at 0.75 × An × fu. A plate 200 by 12 mm with two 22 mm holes yields at 540 kN against rupture at 562 kN, so it yields first, and shows it. With three holes the rupture value is 482 kN. It now breaks without stretching. (The factors shown follow AISC 360; Eurocode 3 uses different partial factors for the same two checks.) A famous connection failure is told in Engineering Drawing & CAD (34/18). Interactive · Drill holes in a tension plate Add holes and make them bigger until the weak link changes from yielding to rupture. Number of holesHole diameter Ductility: the structure should bend before it breaks Strength is not the only thing a structure needs. Ductility is the ability to deform a long way past the first yield while still carrying load. A ductile structure announces its distress with visible cracks and sag, and gives people time to leave. A brittle one fails without warning. Mild steel can stretch by a fifth or more of its length before it breaks; concrete crushes at a strain of 0.3%. That is why codes limit the steel in a beam and demand detailing, such as closely spaced ties, that keeps the concrete confined. Ductility has a second, larger use: earthquakes. Making a building strong enough to stay elastic in the worst earthquake would be absurdly costly. Instead, codes let it yield, because yielding absorbs energy. Seismic codes divide the elastic force by a response modification factor (called R in ASCE 7, with values between about 1.5 and 8, or a behaviour factor q in Eurocode 8) in return for demanding detailing that guarantees the structure can move that far without failing. The rule behind it, the "equal displacement" approximation, says that a yielding structure moves about as far as an elastic one would, so the ductility it needs is about R. If the structure is brittle, the discount is not available. Interactive · Trade strength for ductility Divide the earthquake force by a factor R, then give the structure the ductility to match. Force reduction factor RDuctility the structure has Codes around the world, and one collapse that rewrote them The physics is the same everywhere. The codes are not, because each country writes its own, shaped by its climate, earthquakes, materials, history and legal system. Nearly all use limit-state design with factors. Here are the best-known names. These are examples; always use the code in force where the building stands. RegionConcreteSteelLoads and quakes Europe (EU and many others)Eurocode 2 (EN 1992)Eurocode 3 (EN 1993)EN 1990, 1991, Eurocode 8 (EN 1998) United StatesACI 318AISC 360ASCE 7 MexicoNTC for concreteNTC for steelNTC for earthquake design (Mexico City) ChinaGB 50010GB 50017GB 50011 (seismic) IndiaIS 456IS 800IS 1893 (seismic) JapanAIJ standardsAIJ standardsBuilding Standard Law AustraliaAS 3600AS 4100AS 1170 CanadaCSA A23.3CSA S16National Building Code BrazilABNT NBR 6118ABNT NBR 8800ABNT NBR 6123 (wind) Typical factors show the family resemblance: ASCE 7 uses 1.2 D + 1.6 L and ACI 318 a strength factor of 0.90 for a ductile beam; Eurocodes use 1.35 G + 1.5 Q and partial factors of 1.5 on concrete and 1.15 on reinforcing steel. Different numbers, same method. These are examples, not design values. The code that changed the most was written after a gas leak. At about 5:45 am on 16 May 1968, a resident on the 18th floor of Ronan Point, a 22-storey tower of precast concrete panels in Canning Town, East London, lit a match to light the stove. Leaking gas exploded, and the blast blew out a load-bearing wall panel in the corner of the building. The floors above collapsed and the failure ran down the whole corner of the building. Four people died and 17 were hurt; the toll was low only because most residents were still in their bedrooms. The tower had been completed two months earlier. The Griffiths inquiry found that the design had met the regulations of the day. The regulations simply did not yet ask for what was missing: the floor-to-wall connections were not designed for sideways loads such as an explosion or wind, nothing gave the load another route when one panel was lost, and later investigation found missing bolts and holes stuffed with rubbish. It recommended that new tall blocks be designed to resist progressive collapse. In 1970, the Fifth Amendment to the UK Building Regulations introduced requirements against disproportionate collapse; many other jurisdictions have since adopted similar rules. The same two ideas run through this whole course: ties that hold the structure together, and spare paths for the load when one is lost (34/64). What a graduate just sees Safety factors are not guesses: they are chosen so that the overlap between load and strength is acceptably small, which depends on how much each scatters. Concrete is poor in tension, so the steel goes where the beam stretches: the bottom of a sagging beam, the top of a cantilever. Too little steel in a concrete beam gives warning; too much gives a sudden crush. Prestress puts concrete into compression before the load comes, so that it never cracks. In a steel joint, choose which part is the weakest link, and make it the one that stretches. A structure that can bend a long way before it breaks is worth more than one that is slightly stronger and brittle. Seismic codes trade strength for ductility, which only works if the ductility is really there. Codes differ by country; the physics and the method do not. Where it shows up next Construction Management (34/70): the bars, ties and bolts you specify have to be placed, inspected and paid for on site. Geotechnical Engineering (34/66): footings and piles are reinforced concrete, designed with these same rules. Structural Analysis (34/64): the load combinations and moment diagrams that feed every check here. Vibrations (34/17): ductility is how structures survive an earthquake that exceeds the design force. In the world: every building permit, every bridge assessment, every retrofit that adds ties or confinement to an old building. If you want to go deeper 📘 Standard textbook (concrete): James K. Wight and James G. MacGregor, Reinforced Concrete: Mechanics and Design (Pearson). Careful on why the rules exist. 📘 Standard textbook (concrete): Darwin, Dolan and Nilson, Design of Concrete Structures (McGraw-Hill). Another common text. 📘 Standard textbook (steel): Jack C. McCormac and Stephen F. Csernak, Structural Steel Design (Pearson), or Charles G. Salmon and John E. Johnson's Steel Structures: Design and Behavior. 📘 Gentler: Mario Salvadori, Why Buildings Stand Up (W. W. Norton, 1980), and Henry Petroski, To Engineer Is Human (1985), for how design is shaped by failures. 🌐 Free course: MIT OpenCourseWare 1.054 Mechanics and Design of Concrete Structures, an undergraduate course on reinforced concrete behaviour and design codes. The takeaway A design code is a way of saying how rare a failure must be, turned into factors on loads and strengths. Concrete carries the squeeze and steel carries the stretch, placed where each belongs. Prestress keeps the concrete from stretching at all. In steel, the connection decides, and the weakest link should be the one that yields. Above all, a structure should bend before it breaks, because that is what gives people time, and what lets an earthquake design trade strength for movement. Try this: look at a concrete beam or bridge soffit and ask where the cracks, if any, would run, and which side the steel must be on. Then look at a steel joint and ask which part would give way first. Sources & Further Reading Sources. Limit-state design, the reliability index (β = (μR − μS) ÷ √(σR² + σS²), the normal-distribution model used for illustration), the rectangular stress block for reinforced concrete, the prestress stress formula, the net-section tension check and the equal-displacement rule are standard textbook material. Every number in the interactives was computed in node: β = 2.16 (about 1.5%) at 1.6 × / 12% / 20% and 3.62 (about 1 in 7,000) at 2.2 ×; Φ(−3.8) = 7.2×10⁻⁵; Mn = 257 kN·m for four 20 mm bars at d = 450 mm (a = 82 mm, steel strain 1.07%, f′c 30 MPa, fy 500 MPa, β1 = 0.836 from the ACI formula); the prestress example (P/A = 5.56 MPa, Pec/I = 11.11 MPa, Mc/I = 13.89 MPa: bottom −2.78 MPa, top −8.33 MPa); the plate (540 kN yield against 562 kN rupture with two 22 mm holes, 482 kN with three); and the R = 4, μ = 3 example. The strength factors 0.90 and 0.75 follow AISC 360 for yielding and rupture in tension; the flexure factor and the concrete strength relations follow ACI 318 in simplified form, and the moment-capacity routine solves the section by strain compatibility, so a bar near the compression face gives a small, brittle capacity. The reliability target (EN 1990 reliability class RC2: β = 3.8 for 50 years, failure probability about 7×10⁻⁵) was confirmed from a web search of EN 1990 summaries. Ronan Point facts (16 May 1968, about 5:45 am, 18th floor, 22-storey, Canning Town, 4 dead and 17 injured, Griffiths inquiry, 1970 Fifth Amendment on disproportionate collapse, missing bolts found in later investigation) come from the Wikipedia article on Ronan Point. The code names in the table, the partial factors 1.5 and 1.15, the response modification factor range 1.5 to 8, tensile strength about a tenth of compressive strength, and Freyssinet as pioneer of prestressing are standard engineering knowledge not re-checked against a primary source in this session; the code table lists best-known names only and codes are revised often. The ASCE 7 target reliability was not verified. Go deeper. Wight & MacGregor, Darwin–Dolan–Nilson, McCormac & Csernak, Salmon & Johnson, Salvadori and Petroski are standard titles; check the current edition. The MIT OpenCourseWare page was confirmed. 🌐 MIT OpenCourseWare 1.054 Mechanics and Design of Concrete Structures 📘 Mario Salvadori, Why Buildings Stand Up 📘 Henry Petroski, To Engineer Is Human Step 18. Environmental Engineering: Clean Water In, Clean Water Out ------------------------------------------------------------------- Why: Clean water in, clean water out: the engineering behind public health. The engineering behind the tap, the toilet, the stack and the bin: one accounting rule, a few microbes, and the discovery that most public health happens in pipes. You drink a glass of water from a tap and never think about it. That glass was a river or an aquifer a day or two ago, carrying mud, microbes and farm runoff. A plant you have never seen coagulated, filtered and disinfected it, and a pipe network kept it clean for a long walk to your kitchen. You flush the toilet, and the same thing happens in reverse. Taken together, this is the quiet machine that made cities safe to live in. Before it, cities were where people went to die of cholera and typhoid. This course teaches the one tool that runs the whole machine, and then walks through its main parts: water, wastewater, air and waste. The course in one box Who takes it: civil and environmental engineers (and chemical engineers, in a related form), usually year 3 or 4. Civil programmes shaped by the ASCE Body of Knowledge teach it as the "environmental" leg alongside structural, geotechnical, water and transportation; dedicated environmental programmes build a whole degree on it. What it's for: keeping people healthy by controlling what gets into water, air and soil, and cleaning what already has. Comes after: Fluid Mechanics (34/15), Differential Equations (34/03). Leads to: Hydraulics & Water Resources (34/67) for the same water seen as a resource, and Corporate Finance (34/51, engineering economics) and the Systems Engineering capstone (34/22) for choosing what to build. Honest line: you won't be able to design a treatment plant after this. You will see why every treatment plant is a chain of barriers, why a river can "sag" for oxygen, and what went wrong in Flint. Not medical or regulatory advice: this piece explains how engineering works. Local drinking-water rules differ by country. Mass balance: the one tool that runs everything Matter is not created or destroyed, so for any region you can draw a box around, what builds up inside equals what comes in, minus what goes out, plus what is made, minus what is destroyed. That is a mass balance. A lake, a treatment tank, a room, a river reach, a landfill and the atmosphere over a city are all boxes. The only questions are which flows and reactions you include. V · dC/dt = Q · Cin − Q · C − k · V · Cbuild-up = inflow − outflow − decay This is the "well-mixed tank" version. Two numbers describe it. The residence time is V/Q, the time to flush one tank-volume through. And at steady state, when build-up stops, the concentration settles at Cin · Q / (Q + kV): it is lower if the water is flushed fast or the pollutant decays fast. The same equation shows why a small lake with a slow outflow is easily poisoned, and why it takes so long to recover: the clean-up follows the same exponential curve down, at the same rate. Interactive · A lake with a polluted stream Flow through the lake Natural decay of the pollutant Source switched onSource switched off An invented lake of 5 million m³ and an incoming stream at 100 mg/L of a pollutant. The tank is assumed perfectly mixed. Drinking water: stack the barriers River water is a hazard soup: clay, algae, organic matter, bacteria, viruses, protozoa such as Giardia and Cryptosporidium, and chemicals. No single step removes all of it, so a plant is a chain of barriers, each catching what the others miss. The classic chain is five steps. Coagulation. Tiny particles carry a charge that makes them repel each other and stay suspended for ever. Add an aluminium or iron salt and the charge is neutralised. Flocculation. Gentle stirring lets them collide and glue into flakes (flocs). Sedimentation. In a quiet tank the heavy flocs settle out. Filtration. What is left passes through sand or other media, which catches the finest particles, including most protozoa. Disinfection. Chlorine, chloramine, ozone or ultraviolet light inactivate what remains, and a small chlorine residual guards the pipes all the way to your tap. Success is counted in logs. One "log" removal means 90% gone; two logs, 99%; four logs, 99.99%. Regulators set targets in logs for each organism class (the US surface-water rules, for example, require 3-log for Giardia and 4-log for viruses overall), and every barrier gets credit for its share. Disinfection follows a simple rule of thumb (Chick–Watson): kill, in logs, grows with the concentration of disinfectant times the contact time, called CT. The catch is that organisms differ by orders of magnitude. Viruses fall to a little chlorine; Cryptosporidium shrugs off any practical dose, so for it the plant relies on filtration and UV. And chlorine reacts with natural organic matter to form by-products, so the dose is a balance of two risks, with the harm from leaving microbes in the water much larger. That balance is why the system uses several barriers rather than one big dose. Interactive · Barriers in series VirusGiardiaCryptosporidium Coagulation + filtration Chlorine dose Contact time Illustrative orders of magnitude, not design values: chlorine kill = k · C · t, with k = 0.3 (virus), 0.02 (Giardia), 0.0002 (Cryptosporidium) per mg/L per minute. Real values depend on temperature, pH and the strain. Total capped at 12 logs. John Snow, 1854, and the quiet machine In the 1850s most doctors held that cholera came from miasma, foul air from rot and filth. In late August 1854 an outbreak began in Soho, London, around Broad Street. The physician John Snow, who suspected the water, went door to door, marked the deaths on a map, and found them clustered around one hand pump. The Reverend Henry Whitehead, a local curate who had doubted him, helped trace the source to a household where a baby had cholera and a leaking cesspool near the well. The council removed the pump handle on 8 September. Here the popular story needs two corrections. The outbreak was already declining before the handle came off, and Snow himself said so. And the authorities did not accept his theory at the time: the Board of Health blamed miasma. His stronger evidence was a "natural experiment" from the same year. Two water companies supplied the same streets, house by house; one, Southwark and Vauxhall, drew its water from the sewage-laden lower Thames, while the other, Lambeth, had moved its intake upstream. Households supplied by the first suffered many times more cholera than those by the second, with the poor and the rich on both. Over the following decades London built its sewers and filtered its water, and cholera retreated. Hundreds died in the Broad Street outbreak itself (the number depends on the area and dates counted; one standard summary gives 420). The same story is still being written. The WHO/UNICEF Joint Monitoring Programme reported in August 2025 that about 2.1 billion people, roughly one in four, lacked safely managed drinking water in 2024, including 106 million drinking directly from untreated surface water. It also noted that 2.8 billion people gained safely managed sanitation between 2000 and 2024. The engineering for what remains is known. The hard parts are money, maintenance, institutions and trust. Wastewater and the river's oxygen sag What goes down the drain is mostly water, with organic matter that bacteria will eat. Treating it is a copy of what a healthy river does, in a tank, faster. Primary treatment lets solids settle. Secondary treatment (usually the activated-sludge process) bubbles air through a tank of microbes that consume the dissolved organics. Tertiary steps remove nitrogen and phosphorus and disinfect. The measure of how much "food" the water carries is the BOD (biochemical oxygen demand): the oxygen that microbes will use up, in mg/L, usually measured over 5 days at 20 °C. Why does that matter? Because the microbes breathe the river's oxygen. Put a polluted outfall into a river and two things happen at once. Bacteria consume oxygen as they eat, so dissolved oxygen (DO) falls. The air tops the river back up through its surface at a rate proportional to the deficit from saturation. At first the eating wins and DO sags; downstream the food runs out and the air wins, so DO recovers. That curve is the oxygen sag, described by H. W. Streeter and Earle Phelps in 1925 from Ohio River data collected in 1914–15. dD/dt = k1L − k2Ddeficit grows by decay of waste, shrinks by re-aeration Fish begin to struggle somewhere below 4 to 5 mg/L. Temperature hits this twice. Warm water holds less oxygen at saturation (about 9.1 mg/L at 20 °C and 7.6 at 30 °C, in fresh water at sea level), and warm bacteria eat faster. The same sewage that a river handled in cool weather can push a warm river to zero. Interactive · The oxygen sag Waste load (BOD just below the outfall) Water temperature River speed Re-aeration rate at 20 °C Streeter–Phelps with an initial deficit of 0.5 mg/L and a decay rate k₁ of 0.3 per day at 20 °C. Temperature corrections: k₁ × 1.048(T−20), k₂ × 1.024(T−20). Dashed line: 4 mg/L, an illustrative stress level. Flint, 2014: when the chemistry was changed and the pipes were not On 25 April 2014 the city of Flint, Michigan, as a cost-saving step, switched its drinking water from Detroit's supply to the Flint River. The river water was more corrosive, and the city did not add the corrosion inhibitor (a phosphate) that would have helped the old pipes keep their protective coating. Lead pipes and solder in older homes began to leach lead. Treatment plants do not only clean water; they set its chemistry so that it does not eat the pipes it travels through. What followed is a case study in how a system fails. Residents complained of brown, foul water. In 2015 a Virginia Tech team led by Marc Edwards found that at least a quarter of sampled households exceeded the federal action level for lead of 15 parts per billion, and some homes far more. The paediatrician Mona Hanna-Attisha reported that the share of children with blood lead above 5 µg/dL had risen from 2.4% to 4.9% (in the highest-risk areas, 4% to 10.6%); a later 2018 follow-up found a smaller rise that was not statistically significant after correction, so the size of the effect is itself debated, and this is a good example of a small-numbers problem. There was also an outbreak of Legionnaires' disease in the county, with 87 cases and 12 deaths between June 2014 and November 2015; the state never confirmed a definitive link to the water, though later studies pointed to low chlorine levels. Flint returned to Detroit's water on 16 October 2015. An EPA report in 2017 placed most of the blame for the oversight failures on the state's environmental agency. The last lead pipe in the city's replacement programme was replaced on 1 July 2025, as reported by the Wikipedia article on the crisis. Take three lessons from it. A chemistry change is an engineering change. Regulation that only tests a few samples can miss a problem that families can smell. And data from people who live there should be taken seriously early. Air pollution: a plume spreads sideways and up Smoke from a chimney is a mass balance in motion. The wind carries it downwind and turbulence spreads it sideways and upwards, so the concentration at any point falls the further it goes, and a bell-shaped (Gaussian) profile fits both directions. Ground-level concentration behind a stack depends on the emission rate, the wind speed (more wind, more dilution), the stack's effective height, and above all on how turbulent the air is. In a sunny, unstable afternoon the plume spreads fast; on a calm clear night in a stable layer it stays thin and travels far. C = Q / (π · u · σy · σz) · e−H²/2σz²ground-level concentration on the centre line Raise the stack and the worst ground-level concentration falls steeply, since the plume has to spread much more before it reaches the ground. In the toy below, doubling the stack from 60 m to 120 m in neutral air at 5 m/s cuts the peak by a factor of about 5. But the lesson cuts both ways: a tall stack sends the problem further away rather than removing it. That is why modern rules also limit emissions at the source and why acid rain crossed borders. "Dilution is the solution to pollution" has been abandoned for a reason. Interactive · The chimney plume UnstableSlightly unstableNeutralStable night Effective stack height Wind speed Emission 100 g/s, flat open country, Briggs spreading formulas for Pasquill classes A, C, D and F. A classroom model: real regulatory models (such as AERMOD) add terrain, buildings and weather data. Distance on a log scale. Solid waste: what we throw away still exists Mass balance applies to the bin too. Everything that is bought is eventually thrown away, burned, recycled or composted, and nothing disappears. The standard ladder, the waste hierarchy, ranks the options: prevent waste, reuse, recycle, recover energy, and only then dispose. A modern engineered landfill is a lined, drained, capped box whose rotting waste makes methane, so the engineer collects the gas (some plants burn it for power) and treats the liquid that seeps out, the leachate. The World Bank estimated global municipal solid waste at 2.56 billion tonnes in 2022 and projected 3.86 billion by 2050 if nothing changes, a rise of about 50%. In low-income countries and across Sub-Saharan Africa and South Asia, the great majority is openly dumped, with collection rates in some regions as low as 31% to 67%. The world's waste problem is largely a collection and funding problem, then a design problem. Recycling rates differ widely, and what "recycling" counts depends on the country's rules, so check any number's definition before comparing. What climate changes Environmental engineers design for a "design storm", a "design drought", a "design river temperature". Climate change moves each of them. Warmer water holds less oxygen and speeds the bacteria, so the sag deepens. A warmer atmosphere holds more moisture, so heavy downpours are heavier, which can overwhelm the combined sewers that carry rain and sewage in the same pipe in many older cities, spilling untreated sewage into rivers. Drought leaves less water to dilute a discharge. Rising seas push salt into rivers and aquifers. Treatment itself has a footprint: aeration is energy-hungry, and wastewater and landfills release methane and nitrous oxide. The graduate recognises that the past record no longer sets the design for the future, the same point as in Hydraulics & Water Resources (34/67). What a graduate just sees Draw a box; ask what comes in, what goes out, what is made, what is destroyed. Everything else follows. Residence time V/Q tells you how fast a lake, a tank or a room clears itself, and how slowly it recovers. No single barrier cleans water: coagulation, filtration and disinfection each catch what the last one missed. Treatment is measured in logs, and organisms differ enormously: what kills a virus may do nothing to Cryptosporidium. A river is a reactor: waste eats oxygen, air brings it back, and the lowest point is downstream of the outfall, not at it. Warm water holds less oxygen and bacteria work faster: summer is the dangerous season. Water treatment sets chemistry as well as cleanliness, and Flint shows what happens when you ignore the pipes. A tall stack reduces nearby pollution but only moves it; the fix is at the source. Where it shows up next Hydraulics & Water Resources (34/67): the supply, the river's flow and flood risk that decide how much dilution a river offers. Corporate Finance (34/51, engineering economics) and Systems Engineering (34/22): treatment is a choice between costs, risks and who pays. Chemistry and biology sit under every step here, and a graduate often moves toward public health, ecology or policy. The real world: municipal water and sewer utilities, industrial wastewater, remediation of contaminated sites, air-permit consulting, and humanitarian water supply after disasters. If you want to go deeper 📘 Textbook: Mackenzie L. Davis and David A. Cornwell, Introduction to Environmental Engineering (McGraw-Hill): the standard survey of water, air and waste. 📘 Textbook: James R. Mihelcic and Julie B. Zimmerman, Environmental Engineering: Fundamentals, Sustainability, Design (Wiley): mass balance first, with sustainability and global examples. 📘 Reference for the plant: Metcalf & Eddy, Wastewater Engineering: Treatment and Resource Recovery (McGraw-Hill). 📘 Gentler: Steven Johnson, The Ghost Map (2006): the Broad Street outbreak as a gripping story. ▶️ Free course: MIT OpenCourseWare 1.061 Transport Processes in the Environment (Fall 2008), an undergraduate course on how substances move through rivers and lakes (mass balance, diffusion and exchange). For treatment itself, 1.85 Water and Wastewater Treatment Engineering is a graduate course. 🌐 Easier first read in this library: Hydraulics & Water Resources (34/67) for the water side. The takeaway Environmental engineering runs on one rule: what goes into a box equals what comes out plus what builds up, adjusted for what is made or destroyed. With it you see why a slow lake stays dirty, why drinking water needs a chain of barriers counted in logs, why a river dips in oxygen downstream of an outfall, and why a tall chimney only moves smoke. John Snow showed that water carries disease; Flint showed that the chemistry of water is part of the pipe's design. Climate shifts the numbers the designs were built on. Try this: find your water utility's annual water-quality report (most publish one) and read what they test, the source, and how it is disinfected. Then find where your wastewater goes, and which river or sea receives it. Sources & Further Reading Sources. Mass balance, Chick–Watson, BOD, Streeter–Phelps and the Gaussian plume are standard textbook material (Davis and Cornwell; Mihelcic and Zimmerman). Streeter–Phelps: derived in 1925 from Ohio River data (US Public Health Service, 1914–15), with typical rate constants and temperature factors as in the Wikipedia article. Dissolved-oxygen saturation values are approximate fresh-water values at sea level. John Snow and Broad Street: Wikipedia, 1854 Broad Street cholera outbreak (31 August start, handle removed 8 September, outbreak may already have been declining, Whitehead, the two water companies). Flint: Wikipedia, Flint water crisis (25 April 2014 switch, no corrosion inhibitor, Virginia Tech study, Hanna-Attisha figures and the 2018 follow-up, Legionnaires' 87 cases and 12 deaths, switch back on 16 October 2015, EPA 2017 report, last lead line replaced 1 July 2025). Access figures: WHO/UNICEF JMP, August 2025 (2.1 billion without safely managed drinking water in 2024; 106 million on untreated surface water; 2.8 billion gained safely managed sanitation 2000–2024). Waste: World Bank, What a Waste update (2.56 billion tonnes in 2022; 3.86 billion projected for 2050; collection rates). Log targets for US surface-water rules, the Briggs spreading coefficients and the Chick–Watson constants are textbook values, not re-checked against a primary page for this piece; the constants are labelled illustrative. All interactive numbers were computed in node. 🌐 WHO/UNICEF JMP: household water, sanitation and hygiene, 2025 update 🌐 MIT OpenCourseWare 1.061 Transport Processes in the Environment 🌐 World Bank: What a Waste 📘 Steven Johnson, The Ghost Map Step 19. Transportation Engineering: Why Traffic Jams Form Out of Nothing ------------------------------------------------------------------------- Why: How roads, junctions and transit move people, and why new lanes fill up. Roads behave like a fluid, a signal like a budget of seconds, a curve like a promise to the driver. And adding a lane does not always do what you hope. You are driving along a motorway at 100 km/h. Without any accident, any roadwork or any visible reason, the brake lights ahead come on. Everyone slows to a crawl, then ten minutes later the road opens and you speed up again. You pass the place where it happened and there is nothing there. Where did the jam come from? The answer is the first lesson of this course: a jam is not a thing, it is a wave in the traffic itself. And traffic engineering is the craft of keeping the numbers away from the point where those waves begin, and keeping people alive when they are wrong about them. The course in one box Who takes it: civil engineers (and urban planners and some industrial engineers), usually in year 3 or 4, as one course in highway and traffic engineering. Programmes shaped by the ASCE Body of Knowledge all teach it; planning schools teach the demand and land-use half. What it's for: moving people and goods safely and efficiently through roads, junctions and transit, and understanding how travellers respond to what we build. Comes after: Probability & Statistics (34/05), Differential Equations (34/03). Leads to: Construction Management (34/70) to build it, Corporate Finance (34/51, engineering economics) for judging whether it is worth it, and the Systems Engineering capstone (34/22). Honest line: you won't be able to design an interchange after this. You will see why a ring road makes phantom jams, why green time is a budget, why 30 km/h is a different world from 50, and what the evidence says about "just build more lanes". Speed, density and flow: the fundamental diagram Look at one lane from above. Three numbers describe it. Density (k) is how many cars per kilometre. Speed (v) is how fast they go. Flow (q) is how many cars pass a point per hour. They are tied by pure bookkeeping. q = k · vflow = density × speed Now add one fact about drivers: they slow down as the road fills. At low density, cars go at the free-flow speed and the flow rises with density. But at some point the cars are too close to keep going fast, speed drops, and flow falls even though density keeps rising. The plot of flow against density is the fundamental diagram, and it rises to a peak: the lane's capacity. Past the peak the road is congested: more cars on it means fewer cars getting through. The picture below uses a simple triangle (a free-flow line up to a peak, then a congested line down to jam density) with round, illustrative numbers. Capacity is not an absolute: a motorway lane tops out at very roughly 2,000 to 2,400 vehicles an hour in good conditions (the Highway Capacity Manual in the US publishes the standard values). Engineers grade a road by level of service, from A (free-flowing) to F (breakdown), using density or delay, and the key danger is the knee near capacity, where a small disturbance turns into a long jam. Interactive · The fundamental diagram Cars per kilometre in one lane Illustrative triangular diagram: free-flow speed 110 km/h, peak at 22 cars/km (2,420 cars/h), congestion wave moving backward at 20 km/h, jam density 143 cars/km. Real roads vary. The "state" labels are mine, not Highway Capacity Manual grades. Phantom jams: a jam with no cause Drivers do not react instantly or exactly. One driver brakes slightly harder than needed; the next, with a delay, brakes a bit harder still; and a few cars back the road has come to a stop. That stop-start region moves backwards through the traffic like a wave, while the cars move forwards through it. Near the capacity peak, the traffic is unstable: small hesitations grow, not fade. In 2008 a Japanese team (Sugiyama and colleagues) put vehicles on a circular track, asked the drivers to hold the same speed, and watched. Starting from uniform flow, a jam cluster appeared on its own and travelled backward at about the speed seen on real motorways. Their conclusion: a bottleneck is only a trigger, and the real origin of a jam is the collective behaviour of the vehicles once density passes a critical value. You can see it in a model. The Nagel–Schreckenberg model (1992) is a ring of cells with cars that obey four rules each step: speed up (to a limit), slow down to avoid the car ahead, occasionally brake at random for no reason, and move. Without the random rule, nothing happens below a critical density. With it, jams are born, travel backward, and cost the road a fraction of its flow. Interactive · The ring road Cars on the ring (share of cells) Random hesitation per driver per second New random run Each row is one second, the road runs left to right and wraps around (200 cells of 7.5 m). Red cars are stopped, blue are moving fast. Jams show up as red streaks leaning backward. Maximum speed is 5 cells per second. Signals: green time is a budget A signal gives different streams the junction in turn. Every change of phase wastes a few seconds (the first car to start up is slow, and a safety clearance is needed), so a longer cycle wastes proportionally less. But a long cycle means long waits, with long queues. There is a best cycle length in between. Frank Webster, at Britain's Road Research Laboratory, worked this out in 1958 and gave two classics. The best cycle length is roughly C = (1.5 L + 5) / (1 − Y), where L is the lost time per cycle and Y is the sum of the "critical flow ratios" (the flow divided by the saturation flow, for each phase). The closer Y gets to 1, the faster the best cycle grows, because the junction is nearly full. And green time is split in proportion to those flow ratios. Webster's delay formula, used below, has a bowl shape: too short a cycle wastes time in change-overs, too long makes cars wait. Near the bottom the curve is flat, so the engineer has slack. Interactive · Time the signal Traffic on the first road (per lane) Traffic on the crossing road (per lane) Cycle length Two phases, one lane each, saturation flow 1,800 cars/h, lost time 6 s per phase. Webster's delay formula for random arrivals. Red marker: your cycle. Green line: Webster's best. When a road is full (above about 97% of its capacity) delay is unbounded and the curve is cut off. Geometric design: seeing far enough to stop A road is designed around what the driver can see and do. The master idea is stopping sight distance: if an obstacle appears at the limit of vision, can the driver stop before it? It adds two distances: how far the car travels while the driver notices and reacts, and how far it needs to brake. SSD = 0.278·V·t + V² / (254 · (a/9.81 + G))V in km/h, t reaction time (s), a braking deceleration, G grade The usual North American design values are a reaction time of 2.5 s and a deceleration of 3.4 m/s² (gentle, so that wet or worn tyres still work). At 100 km/h, that gives about 183 m. The first term grows with speed, the second with the square of speed, so at high speed braking dominates. The same square law sets the tightest curve: the radius must satisfy R = V² / (127(e + f)), where e is the banking (superelevation) and f is the sideways friction the tyres can give while keeping passengers comfortable. Codes differ by country, and the numbers here are examples, not design values. Interactive · Can you stop in time? Speed Reaction time Road grade (negative is downhill) Banking e of the curve Braking deceleration 3.4 m/s². The curve uses side friction f = 0.12, an illustrative mid value (design tables use lower values at higher speeds). Try 4 s: it is what you might allow for a driver looking at a phone. Safety: speed, and the Vision Zero idea Roads kill. The WHO estimates about 1.16 million people die each year in road crashes, 92% of them in low- and middle-income countries, which have about 60% of the world's vehicles; more than half are pedestrians, cyclists or motorcyclists. The rate per 100,000 people differs by a factor of ten between countries. Using the WHO estimates compiled on Wikipedia (the years differ): Norway 2.1 (2022), Sweden 2.1 (2022), Japan 2.1 (2023), India 12.6 (2023), Mexico 12.0 (2021), the United States 14.2 (2021), Brazil 15.7 (2021) and Nigeria 20.7 (2019). Countries do not count the same way, and WHO adjusts for that, so national reports can differ. The physical reason speed matters is energy: the kinetic energy goes as v², and the stopping distance as v². The WHO says each 1% rise in mean speed raises the fatal-crash risk about 4%, and the risk of a pedestrian dying when hit by a car front rises about 4.5 times from 50 km/h to 65 km/h. A 2011 AAA Foundation study (Tefft) put the average death risk for a pedestrian hit by a car at about 10% at 23 mph (37 km/h), 25% at 32 mph, 50% at 42 mph, 75% at 50 mph and 90% at 58 mph (93 km/h). Vision Zero begins from this. Sweden's parliament adopted it in 1997: nobody should be killed or seriously injured, because "life and health can never be exchanged for other benefits". The roads should be designed for human error, with the designer sharing the responsibility, and speeds matched to what bodies survive (30 km/h where cars meet pedestrians). Sweden's road deaths fell from 541 in 1997 to 192 in 2021, though progress has stalled since about 2013. Elsewhere results are mixed: Washington, D.C. saw deaths rise in most years after adopting it in 2015, so the label alone changes nothing without road changes. Interactive · A child steps out Your speed Distance when they appear Your reaction time Dry road, braking at 7 m/s² after your reaction time. Death risk interpolated from the Tefft (2011) points above; at low speeds the risk is lower than 10% but not zero. An illustration, not a prediction for any one crash. Induced demand and Braess's paradox The natural response to a congested road is to widen it. Sometimes this works for a while. But the evidence shows that new capacity is largely taken up by new traffic: people drive more, at more times, to farther places, and others move out. This is induced demand. The best-known study, by Gilles Duranton and Matthew Turner (American Economic Review, 2011), found that on US interstate highways in metropolitan areas, vehicle-kilometres travelled rose proportionally with lane-kilometres: an elasticity close to 1. In this case the new lanes fill to about the old level of traffic per lane and congestion is back. Not every study finds that much. A 2004 meta-analysis put the response to a 10% capacity increase at about 4% more driving at once and about 10% within a few years. Mokhtarian and colleagues (2002) found little effect, and critics such as Polzin argue that much extra driving is just population growth or otherwise benign. The reverse is also seen: when roads lose capacity, some trips simply disappear. The safe summary: expect much of the new capacity to fill, more over years than months, and more in growing metropolitan areas; the exact amount depends on the place. Interactive · Widen the road Lane-kilometres added No new traffic (0)Short run (0.4)Long run (1.0) Traffic after = (1 + added share)^ε, where ε is the demand elasticity: 0 means no new trips, 1 means traffic grows exactly with capacity (the Duranton–Turner estimate for US interstates). Traffic per lane is a rough guide to congestion, not a full model. Braess's paradox is a different effect, and often misreported. Dietrich Braess showed in 1968 that adding a road to a network can make everyone's trip longer, if each driver picks the quickest route for themselves. In the textbook example below, 4,000 drivers go from Start to End, two roads take T/100 minutes (T is the number of cars on them) and two take a fixed 45 minutes. Without the shortcut, 2,000 take each route and the trip is 65 minutes. Add a free, instant link and everyone takes the new route: 80 minutes. Nobody can do better by switching on their own, yet everyone is worse off. The paradox is a maths result under particular assumptions: selfish route choice, delays that depend on flow, a certain network shape. Whether real cities contain such links is a separate question. There are reported cases (the Cheonggyecheon expressway in Seoul, Stuttgart in the 1960s, 42nd Street in New York in 1990), but they are mostly anecdotes with other things changing at once. Youn, Gastner and Jeong (2008) found roads in Boston, New York and London where models suggest the effect could occur, and laboratory experiments with real people exist (Rapoport and colleagues, 2009). The safest reading: the paradox is real in theory, plausible in particular places, and no reason to think closing any given road will help. Interactive · Braess's network Drivers No shortcutAdd the shortcut Equilibrium: no driver can save time by switching. Start to A and B to End take T/100 minutes; Start to B and A to End take 45; the shortcut A to B takes 0. With few drivers (under 3,000) the shortcut helps; above that it hurts, until about 9,000, when traffic is so heavy it stops mattering. Transit and land use A lane of cars cannot move many people: at roughly 2,000 vehicles an hour and about 1.2 people per car, it is about 2,400 people an hour. A bus running every two minutes with 80 passengers is similar. A bi-articulated bus every minute with 160 passengers carries about 9,600 an hour, and a metro train of 1,200 passengers every two and a half minutes about 28,800. (All round, illustrative figures: capacity is vehicle size times frequency.) That is the geometric advantage of transit, and it only shows up if the vehicles are frequent and well filled. An empty bus loses to a car on every measure. Filling them is a land-use question. Transit works when many people live and work near stops, walk to them safely, and can reach many places by it; low-density, separated land uses make it expensive. Land use and transport shape each other over decades: roads create sprawl, and sprawl creates traffic. Planning research has examined how density, mix of uses, street design, distance to destinations and distance to transit each affect travel; individual effects are modest, and combinations matter more (Ewing and Cervero's 2010 review is the standard citation). What a graduate just sees Flow = density × speed. Flow peaks at capacity, then falls: a fuller road can move fewer cars. Jams are waves, moving backward, and need no accident: close to capacity, a tiny hesitation grows. Green time is a budget; the best cycle is not the shortest nor the longest, and it lengthens fast as a junction fills. Stopping distance has two parts, and the braking part goes as speed squared; 4 s of distraction is a long road. Energy and stopping distance go as v², so 30 km/h and 50 km/h are different worlds for a pedestrian. Vision Zero designs the road for human error; the label alone does not save lives. New lanes fill with new trips (elasticity near 1 in the best-known study), more slowly than they appear. Braess: in theory a shortcut can slow everyone; in practice be careful about claims. Where it shows up next Construction Management (34/70): road and transit projects are built under traffic, and staging is a traffic problem. Corporate Finance (34/51, engineering economics): the value of time, accident costs and the case for a new line are the whole business case. Probability & Statistics (34/05): arrivals at a junction are random, which is why delay curves look the way they do. The real world: traffic signals and smart corridors, ramp metering, bus rapid transit, school-zone and speed-limit decisions, crash analysis and city street redesigns. If you want to go deeper 📘 Textbook: Fred L. Mannering and Scott S. Washburn, Principles of Highway Engineering and Traffic Analysis (Wiley): the standard undergraduate text. 📘 Textbook: Nicholas J. Garber and Lester A. Hoel, Traffic and Highway Engineering (Cengage): broad, with design standards. 📘 Reference: the Highway Capacity Manual (Transportation Research Board), the US source for capacity and level-of-service methods. 📘 Gentler: Tom Vanderbilt, Traffic: Why We Drive the Way We Do (and What It Says About Us) (2008): a readable tour of jams, signs and driver behaviour. ▶️ Free course: MIT OpenCourseWare 1.201J Transportation Systems Analysis: Demand and Economics (Fall 2008), a graduate course on demand, pricing and project evaluation. 🌐 Easier first read in this library: Probability & Statistics (34/05) for the randomness that makes queues. The takeaway Traffic is a fluid with a peak: flow rises with density to a capacity, then falls, and near that peak jams form from nothing as waves that travel backward. Signals trade lost time against queues, and Webster's best cycle lengthens as a junction fills. Roads are designed around stopping distance, which grows with the square of speed, and the same square law decides who survives a collision. Vision Zero designs for human error. More lanes tend to fill with new trips, and Braess's paradox shows that adding a link can hurt everyone, though real examples are less clear than the maths. Try this: on your next trip, notice when traffic slows with no visible cause, and count how far back the stop-start region reaches. Then on a street near you, pace out the distance a car at 50 km/h needs to stop (about 35 m with a 1.5 s reaction time) and see how much of the street that is. Sources & Further Reading Sources. Flow, density and speed relations, Webster's cycle and delay formulas, stopping sight distance and curve radius are standard textbook material (Mannering and Washburn; Garber and Hoel; AASHTO Green Book). Stopping distance with t = 2.5 s and a = 3.4 m/s² (about 183 m at 100 km/h): Wikipedia, Stopping sight distance. Nagel–Schreckenberg model (1992): Wikipedia. Phantom jams: Sugiyama et al., New Journal of Physics 10, 033001 (2008). Road deaths worldwide, speed and pedestrian risk: WHO fact sheet Road traffic injuries (page dated 20 July 2026); rates by country: WHO estimates compiled on Wikipedia, List of countries by traffic-related death rate (years vary). Pedestrian death risk by speed: B. C. Tefft, Impact Speed and a Pedestrian's Risk of Severe Injury or Death, AAA Foundation for Traffic Safety (2011). Vision Zero, Sweden 541 deaths (1997) and 192 (2021) and the critiques: Wikipedia, Vision Zero. Induced demand: Duranton and Turner, American Economic Review 101(6), 2011; other studies as summarised in Wikipedia, Induced demand. Braess (1968) example, Seoul/Stuttgart/42nd Street reports, Youn et al. (2008) and Rapoport et al. (2009): Wikipedia, Braess's paradox. Not re-checked for this piece: the Highway Capacity Manual capacity range, Ewing and Cervero (2010), Newman and Kenworthy (1989) and the transit capacity round numbers (labelled illustrative). All interactive numbers were computed in node. 🌐 Duranton and Turner, The Fundamental Law of Road Congestion (AER 2011) 🌐 Sugiyama et al., Traffic jams without bottlenecks (New J. Phys. 2008) 🌐 WHO: Road traffic injuries 🌐 MIT OpenCourseWare 1.201J Transportation Systems Analysis 📘 Tom Vanderbilt, Traffic Step 20. Construction Management: Why Big Projects Run Late ----------------------------------------------------------- Why: Scheduling, cost, contracts and safety: how a design actually gets built. A building is a network of tasks, a price with a risk range and a contract that decides who pays. Learn the critical path, the contingency and the overrun data that explain why megaprojects keep running late. A new hospital wing is announced: finished in three years, for a stated price. Four years later it opens, and the bill has grown by a third. Nobody was dishonest in the headline sense, and the engineering was sound. The delay was not a mystery either: a delivery came late, the wrong crew was waiting for a space that was not ready, and the soil under one corner was not what the survey said. Building a structure is the easy half of civil engineering. Getting it built, on time and on budget, by hundreds of people who do not work for one organisation, is the other half, and statistically the half that goes wrong more often. A graduate of construction management does not see a building rise. They see a network of tasks, a crew chart, a price with a risk range, a contract that decides who pays when the soil surprises, and a site where people could be hurt. The course in one box Who takes itCivil engineering students, usually in the final year, and dedicated construction management or construction engineering degrees, which teach it throughout. What it is forPlanning, pricing, contracting and running a construction project so that the design becomes a finished structure without anyone being hurt or ruined. Comes afterEngineering Drawing & CAD (34/18), because drawings are the contract, and Operations Management (34/52), because a site is a process with a bottleneck. Leads toThe capstone, Systems Engineering (34/22), where many parts and stakeholders must agree; and engineering economics, Corporate Finance (34/51), for the question "should we build it at all?". It is the last link after Structural Analysis (34/64), Geotechnical Engineering (34/66), Hydraulics (34/67) and Transportation Engineering (34/69): they say what to build, this says how. What this piece gives youYou will not be able to run a site. You will see why a delay on one task delays everything and another delays nothing, why contingency is a number and not a feeling, and why megaprojects keep overrunning. Easier first read: Judging a Proposal: The Outside View (04/10). The project triangle, and cutting a job into pieces Every project balances three things: scope (what you build), cost and time. Change one and another moves. The picture of a triangle with these corners, and quality in the middle, goes back at least to Martin Barnes, who taught a version of it in 1969. The slogan "good, fast, cheap: pick two" is its folk version. Real projects are subtler than the slogan, because you can trade, but not freely: adding people to a late project does not make it proportionally quicker. The first tool for taming a project is the work breakdown structure (WBS): cut the whole job into pieces until every piece has one owner, one budget and a way of knowing when it is done. A house extension might split into site preparation, foundations, frame, roof, services and finishes. Each in turn splits further. The tree makes the scope visible, so nothing is missed and nothing is built twice, and it is the skeleton for the schedule, the estimate and the contracts that follow. Here is the trade-off in miniature. Tick the pieces of scope; then choose how many people to put on site. More workers finish sooner, but each new worker adds less, because people get in each other's way. The time-related cost (wages and daily site overheads) falls at first as the job shortens, then rises again. Interactive · Scope, crew and the triangle Site prep Foundations Frame Roof Services Finishes Crew on site Crowding Illustrative: each person's productivity is 1 ÷ (1 + crowding × (crew − 1)), wage $300 a day, site overhead $1,500 a day, labour-days and material costs invented for the example (the same logic as Fred Brooks's 1975 observation about adding people to late software projects). Scheduling with the critical path method Once the job is in pieces, each piece has a duration and some must wait for others: you cannot raise the frame before the foundations. Draw that as a network and you have what developers at DuPont and Remington Rand worked out in the late 1950s as the critical path method (CPM), around the same time as the US Navy and Booz Allen Hamilton developed the related PERT. The method is a pair of sweeps. The forward pass finds the earliest each task can start and finish, assuming every predecessor finishes as early as possible. The backward pass finds the latest each can finish without making the whole project later. The difference between them is the task's float. float = latest start − earliest starttasks with zero float form the critical path: the project cannot finish before their total length The critical path is the longest chain of tasks from start to finish. Its length is the project's duration, and a day lost on any task along it is a day lost for the project. Tasks off the path have float: they can slip that many days without consequence. Below is a house in nine tasks, including a window order that takes twelve days to make. Pick a task and change its duration. Shorten a critical task and the whole project shortens (until another path becomes critical). Shorten a non-critical task and nothing happens. Delay "Services" by a few days and watch where float runs out. Interactive · Edit the network ABCDEFGHI Duration of the chosen task Reset Windows arrive late (25 d) Fast-track the roof (3 d) Services slip (12 d) A=Excavate, B=Foundations, C=Frame, D=Roof, E=Services, F=Cladding and windows (needs D and I), G=Interior finishes (needs E and F), H=Handover, I=Windows made (no predecessor). Finish-to-start links only; durations in days; invented for illustration. Resource levelling: the peak is the cost CPM assumes you have as many workers as you need. You never do. Draw the number of workers on site on each day of the schedule as a bar chart (a resource histogram) and you will find peaks and valleys. Peaks cost: the extra crane, the extra bunkhouse, workers with nothing to do the day after. Resource levelling shifts tasks that have float, so that crews overlap less. The finish date does not move, because float is free. The peak does. Here the same nine tasks have crews attached (for instance, six people on the frame, five on services). Left alone and started as early as possible, the roof and services overlap and the site peaks at 11. Services have five days of float, so slide them back and the peak falls. Change a duration in the network above and both charts redraw: the schedule you edit is the one levelled. Interactive · Flatten the peak Crew per task: A 3, B 4, C 6, D 6, E 5, F 4, G 6, H 1, I 0 (made off site). The algorithm tries each start between earliest and latest, keeping the arrangement that minimises the sum of squared daily crews. Illustrative. Real levelling is harder (software uses rules of thumb because the exact problem is famously hard), but the principle is what a graduate keeps: float is a resource you can spend. Cost estimating and contingency An estimate is a prediction made when you know least. The common practice, set out by AACE International, classifies estimates by how much of the design is done: an early concept-stage estimate can easily be out by tens of percent either way, while one made near tender, when drawings are nearly complete, is within roughly ten to fifteen percent. The honest way to state an estimate is as a range, not a single number. The second trap is that the sum of the most likely costs of each package is not the most likely total. Costs are skewed: a package is more likely to overrun by a lot than to come in far below. The sum of the most-likely values is lower than the average outcome. The solution used on serious projects is a quantitative risk analysis: give each package a low, most likely and high value, simulate thousands of combinations (a Monte Carlo simulation), and read off the cost you can be 50% or 80% sure of. The gap above the base is the contingency. Interactive · How much contingency? Size of the upside (overrun) risk How much the packages move together Seven invented work packages with low / most-likely / high costs (in $ thousands) such as Structure 270 / 300 / 420, triangular distributions, 4,000 trials with a fixed random seed. "Move together" means that in that share of trials all packages hit the same percentile (bad weather, a bad market). Illustrative, not a pricing method. Two lessons come out of the simulation. Contingency is a calculated number, not a pad. And correlation matters more than people expect: if everything goes wrong together (a strike, a price spike), the risk does not average out, and the contingency you need grows. Contracts and delivery methods A contract answers the question every project eventually asks: who pays when something goes wrong? How the work is organised determines the answer. Several families are in use around the world: MethodHow it worksTrade-off Design-bid-buildOwner hires a designer; a finished design is tendered to contractors, usually on price.Clear roles, competitive price; slow, and designer and builder do not talk until late. Design-buildOne firm is responsible for design and construction.Faster, fewer design disputes; the owner has less control of detail. Construction manager at riskA builder advises during design, then guarantees a maximum price.Early buildability advice; depends on the quality of the guarantee. Alliance or integrated deliveryOwner, designer and builder share risk and reward under one agreement.Aligns incentives; needs trust and open books. Public-private partnershipA private consortium designs, builds, finances and often operates for decades.Private money and long-term upkeep; the public pays over decades. Standard forms differ by country. FIDIC, the international federation of consulting engineers, publishes widely used international forms (its colour-coded Red, Yellow and Silver books, first issued in 1999 with second editions in 2017). In the UK the NEC suite and JCT forms are common, while US projects often use AIA or ConsensusDocs forms. Public work is tendered under national or regional procurement law, such as EU procurement directives. Price terms vary too: a lump sum fixes the price and puts the overrun on the contractor; cost-plus reimburses actual cost plus a fee and puts it on the owner; a target cost shares the gap. Slide through those three. Interactive · Who carries the overrun? Lump sum Cost plus fee Target cost, shared What the work really costs Contractor's share of the gap (target cost only) Target $10.0 M, fee $0.5 M. Lump sum: owner pays target + fee. Cost plus: owner pays actual + fee. Target cost: owner pays actual + fee + share × (target − actual). Simplified; real contracts add caps, change orders and claims. Illustrative. No form is best. A lump sum rewards efficiency but invites cutting corners and fighting over changes; cost-plus is flexible and gives the builder no reason to save; sharing the gap, as in NEC-style target cost contracts, tries to keep both sides pulling together. Megaprojects and the habit of overrunning The most reliable fact about large projects is that they run over. Bent Flyvbjerg, with Mette Skamris Holm and Søren Buhl, studied 258 transport infrastructure projects worth about $90 billion (published in 2002). Costs were underestimated in about 9 of 10. Averages: rail projects ended 44.7% above their estimate, bridges and tunnels 33.8%, roads 20.4%. In Flyvbjerg's 2014 overview of megaprojects (projects of a billion dollars or more), overruns of up to 50% were common and over 50% not uncommon: the Sydney Opera House 1,400%, the UK NHS IT system 400–700%, Boston's Big Dig 220%, Denver International Airport 200%, and the Channel Tunnel 80%, all in real terms. Overruns had stayed high for the 70 years for which comparable data exist, and all countries with data suffered them. In later work with a database of more than 16,000 projects, he reports that only about 0.5% meet cost, schedule and benefit targets together (secondary reports of his 2023 book). Why? Flyvbjerg's answer is two-part: optimism bias (planners imagine this project going well, a psychological effect Daniel Kahneman and Amos Tversky named the planning fallacy) and strategic misrepresentation (the projects that win approval are those that look cheapest, so there is a selection pressure for under-estimates). The remedy he advocates is the "outside view": ask how projects of this class actually finished, and uplift your estimate accordingly. The idea is not undisputed. Peter Love and Chukwuka Ahiaga-Dagbui argued in 2018 that the study's methods and its "error or lie" framing are flawed. Try the reference-class idea on a budget of your own. Interactive · What would an average project do? Rail Bridge or tunnel Road Budget at the decision to build Contingency you have set aside Averages (rail 44.7%, fixed links 33.8%, roads 20.4%) are from Flyvbjerg, Holm and Buhl (2002), 258 projects, 1927 onward, and describe typical outcomes, not any single project's forecast. Safety: the part with no second chance Construction is among the most dangerous peacetime occupations. The numbers below are official counts, and they are not directly comparable: the US counts by occupation, the EU by accidents at work, Britain by incidents reportable to the regulator, and workforces differ in size. WhereYearWhat was countedSource United States20245,070 fatal work injuries in all jobs; 1,032 of them to construction and extraction workers, about one in fiveBLS Census of Fatal Occupational Injuries European Union2023792 fatal accidents at work in construction, 24.0% of all; 6.3 per 100,000 people employed, third highest among sectors after mining and quarrying (10.8) and agricultureEurostat Great Britain2023/2451 construction workers killed (47 the year before); provisional 2025/26 figures show 25 of 126 worker deathsHealth and Safety Executive In the US, falls, slips and trips are the largest single cause: in 2023 they accounted for 421 of the 1,075 deaths in the construction industry, about 39% (BLS). The International Labour Organization lists falls, electrocution, crush injuries and being caught between objects as the leading hazards worldwide. The standard way to think about them is the hierarchy of controls: first eliminate the hazard (design out the work at height), then substitute, then engineer a control (guardrails, scaffolds), then change procedures, and only last rely on personal protective equipment. Personal protective equipment such as a harness is the weakest layer; designing the hazard out is the strongest. Deaths are a lagging indicator. Good sites watch leading ones: near-miss reports, stop-work authority and whether the newest labourer will tell the foreman something is wrong. That is psychological safety (34/53) with literal consequences. For the engineering side of how designs fail, see the Hyatt Regency case in Engineering Drawing & CAD (34/18). BIM: building it twice, once in the computer Building information modelling (BIM) replaces a set of 2D drawings with one shared 3D model whose elements carry data: this beam, this duct, this supplier, this cost. The first payoff is clash detection: finding in the model that the ventilation duct runs through a structural beam, which is cheaper than finding it on site with the crane booked. Later the same model feeds estimates, the construction sequence and the facility manager. ISO 19650, published from 2018, sets an international framework for managing information on a BIM project, building on earlier British standards. Governments have pushed it: the UK set out a strategy in 2011 aimed at collaborative 3D BIM on its projects by 2016, and Hong Kong made it mandatory for government projects above HK$30 million from 1 January 2018. A 2004 NIST study estimated that poor interoperability cost the US capital facilities industry $15.8 billion a year. Cautions are fair, too: adoption outside rich countries has been slow, a 2019 claim of savings up to 20% was an official estimate and not a measured result, and the industry warns about "BIM-washing", a label on work that does not use it properly. The model is only as good as the information put in. What a graduate just sees A project is scope, cost and time, held in tension; change one and another moves. Only the critical path sets the finish date. Everything else has float, and float can be spent. The peak crew, not the average, drives the cost of labour and space. The most likely cost of each part does not add up to the most likely total. Contingency is a computed number. A contract is a risk allocation: ask who pays when the soil surprises. Big projects overrun on average. The right question is "how did projects like this end?" A fall is the likeliest killer on a site; remove the hazard before you issue the harness. The drawing is a model; the more information it carries, the earlier clashes show up. Where it shows up next Systems Engineering & the Capstone (34/22): a project is many parts and stakeholders agreeing; the schedule, risk register and interface list are the same tools. Operations Management (34/52): bottlenecks, queues and buffers on a site are the same laws; the critical path is a bottleneck in time. Structural Analysis (34/64), Geotechnical (34/66), Hydraulics (34/67), Transportation (34/69): they supply the design; this decides whether it can be built. Corporate Finance (34/51): discounting, risk and the choice between building and not. In the world: general contractors, owner's representatives, quantity surveyors, highway agencies, developers and every public project that has made the news for being late. If you want to go deeper 📘 Textbook. Chris Hendrickson, Project Management for Construction (Carnegie Mellon, available free online): scheduling, estimating, contracts and safety, written for civil engineering students. 📘 Textbook. Daniel Halpin and Bolivar Senior, Construction Management (Wiley): the standard undergraduate text, with worked CPM and cost examples. 📘 Gentler. Bent Flyvbjerg and Dan Gardner, How Big Things Get Done (Currency, 2023): readable, data-driven account of why projects fail and what the successful ones do. 🌐 Free course. MIT OpenCourseWare 1.040 Project Management (Spring 2009, Moavenzadeh): an undergraduate civil engineering course on managing construction and infrastructure projects. The takeaway Break the job into pieces, find the critical path, and spend the float. Level the crews so the peak does not set the bill. Price with a range, and hold a contingency that a simulation, not hope, supports. Choose a contract that puts risk with whoever can best control it. Look at how similar projects actually ended before believing yours will differ. And design out the fall before you hand out the harness. Try this: the next time you pass a building site, pick one activity you can see and ask what must happen before it can start and what is waiting on it. You have just drawn a bit of a network. Sources & Further Reading Sources. Project triangle: Martin Barnes, 1968 thesis and 1969 course, as described in the Wikipedia article on the project management triangle. CPM and PERT: late 1950s, Walker and Kelley (DuPont, Remington Rand); Booz Allen Hamilton and the US Navy. Cost estimate classes: AACE International (ranges are described in rough terms here). Overruns: Flyvbjerg, Holm & Buhl, "Underestimating Costs in Public Works Projects: Error or Lie?", Journal of the American Planning Association 68(3), 2002 (258 projects, about US$90 billion; rail 44.7%, fixed links 33.8%, roads 20.4%); Flyvbjerg, "What You Should Know About Megaprojects, and Why: An Overview", Project Management Journal (2014); Flyvbjerg & Gardner, How Big Things Get Done (2023), figures for the 16,000-project database taken from secondary reports; critique: Love & Ahiaga-Dagbui, Transportation Research Part A 113 (2018). Safety: US Bureau of Labor Statistics, Census of Fatal Occupational Injuries (2024 release; 2023 data for the falls figure); Eurostat, accidents at work by economic activity (2023); UK Health and Safety Executive, work-related fatal injuries in Great Britain. BIM: ISO 19650; NIST GCR 04-867 (2004) as cited in the Wikipedia article on BIM; national mandates as described there. FIDIC, NEC, JCT, AIA and ConsensusDocs forms are named from general knowledge. Interactive figures were computed and checked in code; datasets are illustrative. 🌐 MIT OCW 1.040 Project Management 📘 Hendrickson, Project Management for Construction 📘 Halpin & Senior, Construction Management 📘 Flyvbjerg & Gardner, How Big Things Get Done Step 21. Corporate Finance: The Price of Time and Risk ------------------------------------------------------ Why: Should we build it? Engineering economics is finance with a long horizon. A pound today is worth more than a pound next year, and a risky pound is worth less than a safe one. Corporate finance is the craft of putting a number on both, so that a company can decide which projects are worth their money. Someone offers you $1,000 today or $1,100 in a year. Almost nobody needs a course to say "it depends": on whether you trust them, on what else you could do with the money, on whether you need it now. Those three words, trust, alternatives and need, are the whole subject. Corporate finance is what happens when you ask the question about a factory instead of a loan, with a hundred million at stake and a board waiting for the answer. Every company faces the same three decisions. Which projects should we spend on? Where does the money come from, borrowed or from owners? And what is the whole business worth? A graduate has one tool for all three: a way of turning cash that arrives at different times, with different risks, into a single number today. The course in one box Who takes itBusiness, economics and engineering students, usually in year 2 or 3; every MBA takes it in the first year, and it sits in Wharton's undergraduate business core and in Harvard Business School's Required Curriculum. Finance majors then split it into investments, derivatives and banking. What it is forDeciding whether a project, a company or a funding plan is worth more than it costs, once time and risk are priced in. Comes afterAccounting (34/50) for the cash flows, and Probability & Statistics (34/05) for risk. Leads toEntrepreneurship (34/56). Civil and other engineers meet it as "engineering economics" when they ask whether a bridge is worth building. What this piece gives youYou will not be able to value a real company. You will be able to compute a net present value, see why a headline return can mislead, and know what a discount rate is made of. Easier first reads in the library: Compound Interest Will Blow Your Mind (03/03) and Investing Without Gambling (03/10). This piece goes further. It describes how finance works and what economists argue about; it is not investment advice. Money has a time stamp: compounding and present value If a bank pays 6% a year, $1,000 today becomes $1,060 next year and $1,123.60 the year after, because the interest earns interest. That is compounding. Run it backwards and you get its twin, present value: what a future sum is worth now. At 6%, a promise of $1,000 in 20 years is worth only about $312 today, because $312 invested at 6% would grow to $1,000 by then. FV = PV × (1 + r)n   ·   PV = FV ÷ (1 + r)nr = interest rate per period, n = number of periods. The rate that turns future money into today's money is called the discount rate. Interactive · Grow it, then discount it Amount Interest or discount rate Years At the default, $1,000 grows to $3,207 in 20 years, and $1,000 due in 20 years is worth $312 now. The rule of 72 gives a shortcut: money doubles in roughly 72 divided by the percentage rate, so about 12 years at 6% (the exact answer is 11.9). The same arithmetic is why waiting costs: a rate is the price of time. Net present value: the one rule for saying yes A project is a stream of cash: an outlay now, receipts later. To judge it, discount each receipt back to today and add them up. The sum, minus the outlay, is the net present value (NPV). If NPV is positive, the project earns more than the discount rate demands, and doing it makes the owners richer by exactly that amount. If it is negative, they would do better lending the money at the discount rate. NPV = −Outlay + CF1/(1+r) + CF2/(1+r)2 + … + CFn/(1+r)naccept when NPV > 0. The discount rate is your opportunity cost for projects of this risk. Worked: spend $1,000 now for $300 a year over five years. At an 8% discount rate the five receipts are worth $1,198 today, so NPV is +$198. At 15% the receipts are worth $1,006, and NPV is nearly zero. The rate at which NPV hits zero is the project's internal rate of return (IRR), here 15.2%. Slide the rate and watch the line cross zero. Interactive · NPV against the discount rate Outlay today Cash per year Years Discount rate When IRR misleads IRR is popular because it is one percentage, easy to compare with a hurdle rate. It is also the part of finance where the textbooks warn you most. Two classic failures. Several IRRs. If the cash flow changes sign more than once, for example a mine that needs a big cleanup bill at the end, the NPV curve can cross zero twice. A project with flows of −100, +230 and −132 (in millions) has IRRs of both 10% and 20%, and it is worth doing only when the discount rate lies between them. A single IRR number hides that. Scale. Project A turns 100 into 150 (IRR 50%). Project B turns 1,000 into 1,300 (IRR 30%). If you can do only one and the cost of capital is 10%, A adds 36 and B adds 182. Higher IRR, much lower value. The two NPV curves cross at about 27.8%, and below that rate B wins. The rule graduates keep: when IRR and NPV disagree, trust NPV, because it measures money, not percent. Interactive · Two traps Cleanup bill at the end Small and fast vs big and slow Discount rate Illustrative cash flows (millions), invented to show the traps. Risk and return, and why diversification is nearly free So far the future was certain. It is not. A share can fall 40% in a year, and investors want to be paid for taking that chance. The central observation of modern finance, made by Harry Markowitz in 1952, is that you are paid for the risk you cannot avoid, and that you can avoid a lot of it for nothing. Hold many assets that do not move in lockstep and their ups and downs partly cancel. Risk here is measured as volatility: the standard deviation of returns. For a portfolio of N assets, each with volatility σ and an average pairwise correlation ρ, the portfolio's volatility is: σp = σ × √( 1/N + (1 − 1/N) × ρ )diversification removes the 1/N part; the ρ part, risk shared by all assets, is what is left Take stocks with 30% volatility each and a correlation of 0.3. One stock: 30%. Ten stocks: 18.2%. Fifty: 16.8%. The floor is 30% × √0.3 = 16.4%, no matter how many you hold. That floor is market risk, which cannot be diversified away; the rest is specific risk, which can. Interactive · How many stocks is enough? Each stock's volatility Correlation between stocks Number of stocks held Assumes equal volatility and equal pairwise correlation, a simplification. Real correlations vary and rise in crises. CAPM: the idea, and its critics If investors diversify, the only risk worth a reward is the market risk left over. The Capital Asset Pricing Model, developed in the early 1960s by Jack Treynor, William Sharpe, John Lintner and Jan Mossin (Sharpe's version appeared in 1964), turns that into a price. A stock's beta measures how strongly it moves with the market: beta 1 moves with it, 0.5 half as much, 1.5 one and a half times. The model says the required return is the safe rate plus beta times the market's extra reward: Required return = rf + β × (market return − rf)the market return minus rf is called the equity risk premium Interactive · The security market line Risk-free rate Equity risk premium Beta Inputs are illustrative, not forecasts. Estimates of the premium range widely. With a 3% safe rate, a 5% premium and a beta of 1.2, the required return is 9%. That number is the discount rate you use in NPV. The model is elegant, which is why it is taught everywhere, and why it is argued about. Critics: Fama and French (1992) found that beta alone did a poor job of explaining average returns across US stocks, while company size and the ratio of book value to price did better, and Fama and French later added those as extra factors. Richard Roll (1977) argued the model cannot truly be tested, since the "market" should hold every asset, including property and human skills. Many practitioners use CAPM anyway as a rough, transparent starting point. Treat the number as an estimate with a wide error bar. Cost of capital and valuation: discounted cash flow A company's cost of capital is the return its funders demand: lenders want interest, owners want a CAPM-style return. Blend the two, weighted by how much each provides, and adjust the debt for tax, and you get the weighted average cost of capital (WACC). It is the discount rate for the firm's cash flows. A discounted cash flow (DCF) valuation then does what NPV did for projects, but for the whole firm. Forecast the free cash flow for some years, add a terminal value for everything after, and discount it all. With this firm's cash of $100 million, 8% growth for five years, 2% afterwards and a 9% WACC, the business is worth about $1.88 billion, and 74% of that comes from the terminal value. That is the first lesson of DCF: most of the answer is the distant future, so the answer is only as good as the growth and discount-rate assumptions. Move the terminal growth up to 4% and the value jumps by almost a third. Interactive · Value a company Growth for the next 5 years Growth forever after Discount rate (WACC) A company with $100 million of free cash flow now. If growth forever would reach the discount rate, the tool caps it one point below. The other common method uses multiples: look at what similar companies sell for relative to a measure like earnings or cash flow, and apply the same ratio. The tool above would call this firm's value "18.8 times cash flow". Multiples are quick and market-based, but they only move with the market: if the whole sector is overpriced, so is yours. Careful analysts use both and ask why they disagree. Capital structure: Modigliani and Miller in words Should a firm fund itself with debt or equity? In 1958 Franco Modigliani and Merton Miller proved something startling: in a world with no taxes, no bankruptcy costs and no information gaps, it does not matter. The firm's value comes from what its assets earn, not from how the claims on them are sliced. An investor could borrow personally to copy a leveraged firm, so the market will not pay extra for the firm doing it. More debt makes the equity riskier, so shareholders demand a higher return that exactly cancels the cheapness of debt. They received Nobel prizes for this (Modigliani 1985, Miller 1990). The theorem is a benchmark: it tells you where to look for real effects. The first is tax. In most countries interest is deductible and dividends are not, so debt carries a tax shield and a levered firm's WACC falls. The second is distress: more debt raises the chance of bankruptcy, which is costly. The textbook "trade-off theory" balances the two; pecking order theorists argue managers instead prefer internal cash, then debt, then new shares, because of information gaps. Both are live in the literature. Interactive · Debt, tax and the cost of capital Debt as a share of firm value Corporate tax rate Assumes an all-equity return of 10% and debt at 5%. Distress costs are not modelled, which is why the tax benefit looks endless here. Real firms hit a limit. With no tax, the WACC stays at 10% at any debt level while the cost of equity climbs (13.3% at 40% debt). With a 25% tax rate, 40% debt cuts WACC to 9%. Pushed to 90% the tool says 7.75%, which no lender would permit: that is the limit the distress cost supplies. Efficient markets: the debate that never ended If prices already reflect what everyone knows, nobody can reliably beat the market, and an NPV calculation is just reading the price. Eugene Fama's 1970 review gave this efficient market hypothesis its structure: weak (past prices carry no news), semi-strong (public information is already in the price) and strong (even private information is). The evidence is mixed in a way that keeps both camps in business. For: most actively managed funds trail a plain index over long periods, which the S&P SPIVA scorecards have reported year after year across many markets. Against: Robert Shiller showed in 1981 that stock prices move far more than later dividends justify, and bubbles and crashes seem hard to explain with rational pricing. Behavioural economists, such as those in Behavioral Economics (34/46), build on this. A middle path: Sanford Grossman and Joseph Stiglitz argued in 1980 that markets cannot be perfectly efficient, because then nobody would pay to gather the information that makes them efficient. Fama, Shiller and Lars Peter Hansen shared the 2013 Nobel prize in economics: a prize for one side's theory, the other side's critique and the statistical methods both used. Luck matters more than people think. Try it: give a thousand "fund managers" a coin each and see how many beat the market ten years in a row by chance alone. Interactive · Skill or coin flips? Managers Years in a row Flip the coins Each manager has a 50% chance of beating the market each year, independent of the others. The simulation uses your browser's random numbers. For 1,000 managers and ten years the expected number with a perfect record is about one (1,000 ÷ 1,024), and there is a 62% chance of at least one. That one looks like a genius. Real streaks are not purely luck, but the arithmetic is why careful readers ask how many managers started. What a graduate just sees A sum of money has a date on it. Never add cash flows from different years without discounting. A discount rate is an opportunity cost: what else the money could earn at the same risk. NPV measures value created in dollars; IRR is only a percentage and can mislead, especially with scale or sign changes. You are paid for risk you cannot diversify away, not for risk you could have removed for free. A valuation is mostly about the distant future, so ask which assumption carries it. Debt is cheaper because it is safer for lenders, which makes equity riskier. The free lunch is in the tax code, not in the maths. Beating the market is hard; many apparent winners are lucky survivors. Every number has a model behind it. Know which, and what it assumes. Where it shows up next Entrepreneurship (34/56): equity, dilution and the valuation of a company that does not yet earn. Economics and Macroeconomics (34/42, 34/43): the interest rate is set by savers, firms and central banks; every discount rate in this piece sits on top of it. Behavioral Economics (34/46) and Game Theory (34/45): where the efficient-market story strains, and how markets with strategic players behave. In the world: a mortgage, a pension, a merger announcement, a startup funding round, a new power plant, a government deciding whether a railway pays back. If you want to go deeper 📘 Textbook. Richard Brealey, Stewart Myers, Franklin Allen and Alex Edmans, Principles of Corporate Finance (McGraw Hill, 14th edition, 2022). Strong on intuition and real cases. 📘 Textbook. Jonathan Berk and Peter DeMarzo, Corporate Finance (Pearson, current edition). More systematic, built around no-arbitrage reasoning. 📘 Gentler. Peter Bernstein, Against the Gods: The Remarkable Story of Risk (Wiley, 1996) for the history of risk; Burton Malkiel, A Random Walk Down Wall Street (Norton) for the efficient-markets argument in plain language. 🌐 Free course. MIT OpenCourseWare 15.401 Finance Theory I (Fall 2008, Andrew Lo): lectures and notes on present value, risk, portfolio theory and valuation. The takeaway Corporate finance turns cash that arrives at different times and different risks into a single number today. Discount at a rate that reflects what the money could earn elsewhere at the same risk; trust NPV over IRR; remember that diversification removes risk for free and the rest is what you are paid for. Every valuation is a model, and its biggest input is usually a guess about the far future. Try this: take any offer with delayed money (a "0% for 12 months" deal, a pension quote, a subscription that saves 20% if you pay a year ahead). Pick a discount rate you believe, and compute the present value of both sides. Sources & Further Reading Sources. Time value, NPV, IRR and the diversification formula are standard textbook finance (Brealey et al.; Berk & DeMarzo). Markowitz, "Portfolio Selection", Journal of Finance (1952); Sharpe, "Capital Asset Prices", Journal of Finance (1964), with Treynor, Lintner and Mossin independently; Roll, "A Critique of the Asset Pricing Theory's Tests" (1977); Fama and French, "The Cross-Section of Expected Stock Returns", Journal of Finance (1992); Modigliani and Miller, "The Cost of Capital, Corporation Finance and the Theory of Investment", American Economic Review 48(3) (1958); Fama, "Efficient Capital Markets: A Review of Theory and Empirical Work" (1970); Shiller (1981) on excess volatility; Grossman and Stiglitz (1980). Nobel prizes: Markowitz, Miller and Sharpe 1990; Modigliani 1985; Fama, Hansen and Shiller 2013. The SPIVA claim refers to S&P Dow Jones Indices' published scorecards; the share varies by year and market, so no percentage is quoted. Interactive numbers were computed in node: $1,000 at 6% for 20 years = $3,207.14, present value $311.80, doubling 11.90 years; NPV of 5 × $300 on $1,000 at 8% = $197.81, IRR 15.24%; the −100/+230/−132 flows have IRRs of exactly 10% and 20%; A versus B at 10% = $36.36 against $181.82, crossover 27.8%; diversification 30.0%, 18.2%, 16.8% and floor 16.4%; CAPM 3% + 1.2 × 5% = 9%; DCF enterprise value $1,877.9 million with terminal 74.1%, 2,472.7 at 4% terminal growth; WACC 9.0% and cost of equity 12.5% at 40% debt and 25% tax; 1,000/1,024 = 0.98 managers, 62.4% chance of at least one. 🌐 MIT OCW 15.401 Finance Theory I 📘 Brealey, Myers, Allen & Edmans, Principles of Corporate Finance 📘 Bernstein, Against the Gods 📘 Malkiel, A Random Walk Down Wall Street Step 22. Systems Engineering & the Capstone: Many Parts, One Product -------------------------------------------------------------------- Why: A project is many subsystems and stakeholders agreeing; that's the capstone. A product is many subsystems agreeing with each other. Most famous failures happened on the lines between them. A phone in your pocket does several things at once. The radio talks to a tower, the screen glows, the processor thinks, the battery pushes current through all of it, and the whole slab has to stay cool enough to hold. Every one of those parts was designed by a different team. Yet when you press a button, they agree. Now imagine one of them quietly disagrees. One team counts in one unit, another in a different one. Nothing is wrong with any single part. Everything is wrong with the product. This subject is about that gap: not how to design a part, but how to make a hundred parts mean the same thing. Engineers who have taken it carry a habit for life. When someone proudly shows a clever component, they ask three questions before admiring it: what does it plug into, what does it assume about its neighbours, and what happens when it is wrong? The course in one box Who takes it: final-year mechanical and aerospace engineers, usually as a capstone (a year-long team project), and in many universities every engineering major has a version of it. What it is for: turning a pile of well-designed parts into one product that meets what someone actually needed, on a budget and a schedule. Comes after: Machine Design (34/19), Control Systems (34/20) and Software Engineering (34/31), which is also a Computer Science subject. Leads to: nothing on the map; it is the last stop of the Mechanical & Aerospace path. What it leads to is the job, where most of the work is this. What this piece gives you, and what it doesn't: you won't be able to run a real programme or write a requirements document an auditor would accept. You will see why a spacecraft is weighed more often than it is admired, and why the famous failures were almost all about neighbours disagreeing. A product is many subsystems agreeing A system is a set of parts whose combined behaviour is something none of them can do alone. A car is not an engine plus wheels. It is an engine that must fit in the space the body left, make heat the radiator can reject, and shake at frequencies the mounts can absorb. Each sentence in that list is an agreement between two subsystems. Systems engineers write those agreements down in two forms. A requirement says what the whole thing, or one part, must do, in a sentence a test can settle: "The device shall have a mass of no more than 12 kg." A good requirement can be checked by measuring; "lightweight" cannot. An interface says exactly what passes between two parts: the bolt pattern, the voltage and connector, the message format, the units. A document that holds these agreements is called an interface control document, and it is the most boring and most valuable thing in the project. Most of a system's trouble is not inside the parts. It is on the lines between them. The craft is a loop: break the whole into subsystems, give each a budget and an interface, let teams work, then bring the results back and see whether the whole still works. The V: build down, prove up The classic picture is a letter V. Down the left side you go from need to detail: what the customer wants, then system requirements, then the design of each subsystem, then the parts. At the bottom you build. Up the right side you test, and the key idea is that each level of testing answers to the level on the left at the same height. You test the part against the part's design, the assembled subsystem against its requirements, the whole product against the system requirements, and finally the product in real use against the original need. That gives the two words every systems engineer uses precisely. Verification asks whether you built it right: does it meet the requirements as written? Validation asks whether you built the right thing: does it actually solve the problem the person had? A product can pass every verification and still fail validation, if the requirements were wrong. The V-model has several roots: it appeared in the late 1980s in Germany (for defence procurement) and in the United States, and it is a teaching picture, not a law. Real projects loop back and forth across it, and agile software teams (see Software Engineering, 34/31) spiral through small versions of it. Interactive · Walk the V Now classify. Is each of these verification or validation? Worked example: a requirement says the part "shall have a mass of no more than 12 kg". You put it on a scale and read 11.6 kg. That is verification, and the requirement is met with 0.4 kg (about 3%) to spare. Whether anyone can lift a 12 kg part into place is a validation question, and the scale cannot answer it. Trade studies: choosing without a perfect answer Early in a project there are always several ways to do a job, and none is best at everything. A trade study scores each option against the criteria that matter, weights the criteria by how much they matter, and adds it up. The sum is not the point. The point is that the weights make the argument visible: two people who disagree about the answer are usually disagreeing about the weights, and now they can say so. score = (w1s1 + w2s2 + …) ÷ (w1 + w2 + …) Interactive · Wheels, tracks or legs? A small delivery robot. The scores (1 = poor, 10 = excellent) are invented for the exercise; the point is what the weights do. Drag a weight and watch the winner change. Try it: with the starting weights, plain wheels win comfortably. Put nearly all the weight on performance over rough ground and legs edge ahead. Nobody's arithmetic changed, only what they cared about. A mature team also checks sensitivity: how far can a weight move before the winner flips? If a tiny change flips it, the study has not decided anything yet. Budgets and margins: mass, power and heat Hardware obeys sums. A spacecraft has a mass budget (everything added up against what the rocket can lift), a power budget (everything it draws against what the panels and battery supply) and a thermal budget (everything it dissipates against what it can radiate away). Each subsystem gets an allowance. The system engineer's job is to add them up, often, and to be suspicious of the total. The suspicion is earned. Designs get heavier. NASA's systems engineering handbook, in an example, says spacecraft dry mass can grow by as much as 25 to 30 percent between detailed design and launch. So teams hold a margin: they add a deliberate percentage to the current best estimate, larger when the design is young and shrinking as parts are actually built and weighed. A margin is not padding. It is an honest admission of what you don't know yet. Interactive · A mass budget with margin An invented small satellite with a 500 kg launch limit. Its best estimate is the sum of the parts below. Choose how much margin you hold, then how much the design really grows. Margin you hold Growth the design really shows The best estimate here is 415 kg. With 20% margin the budget reads 498 kg, comfortably under 500. Then the design grows by 25%: 415 × 1.25 = 518.75 kg, and the satellite no longer fits. The margin was not too generous; the growth was bigger than the margin. That is how real programmes learn their margins were thin, usually late. The budgets are coupled, which is the deeper lesson. Take mass in a rocket. The rocket equation says that the ratio of liftoff mass to final mass is e raised to (Δv divided by exhaust speed). One extra kilogram of dead weight at the end therefore costs several kilograms of propellant at the start. mstart ÷ mend = eΔv ÷ (Isp × g0) Interactive · The snowball A stage with 1,000 kg of structure and payload that must change its speed by 4,000 m/s. The engine's specific impulse is 320 s (a good chemical engine). Add dead weight and see what it costs at the start. Extra dead weight Velocity change needed At 4,000 m/s the ratio is 3.58, so each extra kilogram at the end costs 2.58 kg more propellant, and 3.58 kg of extra liftoff mass in all. Ten kilograms of "just a bracket" becomes about 36 kg on the pad. This simple model ignores that the bigger tanks weigh more too, so reality is worse. This is why aerospace engineers fight for grams. Heat has a vacuum-sized version of the same coupling. In space a spacecraft can only get rid of heat by radiating it. The hotter and larger its radiating surface, the more it sheds, following the fourth power of temperature: P = ε σ A T4  →  A = P ÷ (ε σ T4) Interactive · Where does the heat go? Every watt the electronics use comes out as heat. With a radiator at 300 K (27 °C), emissivity 0.85, and ignoring sunlight, how much area must it have? Heat to reject Radiator temperature At 300 K, 100 W needs about 0.26 m² of radiator. So a power decision, a few more watts for a faster computer, becomes a thermal decision (more area or a hotter radiator) and a mass decision (more structure to carry it). One number changed and three subsystems moved. That chain is what "system" means. Interfaces are where the failures live Two famous losses show this better than any lecture, and both are now textbook cases. On 23 September 1999, NASA's Mars Climate Orbiter lost contact during its orbit-insertion burn. The investigation found the root cause was a failure to use metric units in a piece of ground software. A ground file that described the effect of the spacecraft's small thruster firings reported it in pound-force seconds. The navigation software that used them expected newton-seconds, as the specification said. One pound-force is 4.448 newtons, so every number was read about 4.45 times too small. The board noted that the project's own checks did not catch it: it cited a lack of complete end-to-end verification of the navigation software and models. It was meant to pass Mars at about 226 km. Calculations afterwards put it as low as about 57 km, below the roughly 80 km thought survivable. On 4 June 1996, the first Ariane 5 rocket, flight 501, veered and broke up about 37 seconds after the launch sequence began, roughly 30 seconds after liftoff. The inquiry, chaired by Jacques-Louis Lions, traced it to the inertial reference system, which had reused software from Ariane 4. A number describing horizontal motion, held as a 64-bit floating-point value, was converted to a 16-bit signed integer. Ariane 5 flew a faster early trajectory than Ariane 4, the value no longer fitted, and the conversion raised an error that was not handled. The two redundant units had identical hardware and software, so both shut down together, and the rocket's computer then read diagnostic data as if it were real guidance. The board's lesson was not "a programmer blundered". It was that an assumption valid for the old rocket was never tested against the new one: the board reported that the new rocket's trajectory data was left out of the unit's specification, and the units were not part of the full closed-loop simulation. Three of the seven vulnerable conversions had been left unguarded on the belief that the values were physically limited, and that belief proved wrong for this one. Interactive · The two interface lab 1. Units. Slide the real impulse. The ground file gives a number in pound-force seconds; navigation reads it as newton-seconds. Real impulse of a burst 2. Range. A 16-bit signed integer holds -32,768 to 32,767. Slide the value being converted, and choose whether the conversion is guarded. Value to convert Both stories share a shape. Every part worked for the conditions it was designed for. The failure sat in an assumption between parts. Systems engineers answer this with the same few habits: write the interface down, put units and ranges in it, test the parts together as early as possible, and test with real data rather than the data you expect. FMEA: asking how it breaks before it does Failure mode and effects analysis, or FMEA, is a structured way of asking "how could this break, and what happens then?" It was first described in a US military procedure in 1949, and by the early 1960s NASA contractors were using versions of it on programmes such as Apollo. For each part, a team lists the ways it can fail, the effect of each, and how it could be noticed. Then they rate three things from 1 to 10: how severe the effect is, how often it will happen (occurrence), and how hard it is to detect before it hurts someone. Multiply the three and you get a risk priority number (RPN), used to rank what to fix first. Interactive · Rank the risks An electric scooter. The ratings are invented for the exercise. Pick a failure and move its three ratings (for detection, 10 means hardest to detect). Severity Occurrence Detection (10 = hardest) The starting numbers show the catch. Battery overheating while charging has the worst possible consequence (severity 10), yet its RPN of 60 ranks it third, below a loosening bolt that scores 192. That is not an arithmetic error: multiplying scales hides a very bad but rare event behind a mild but common one. Critics have said exactly this, and the 2019 joint automotive handbook of the AIAG and VDA replaced the RPN with an "action priority" table that looks at severity first. The honest use of FMEA is the conversation it forces, not the number. Integration, test and what the capstone teaches The right-hand side of the V is where the plan meets the world. Integration means putting the parts together a little at a time, not all at once on the last week, so a surprise is found when it is cheap. Test is the evidence: vibration tables for launch, thermal-vacuum chambers for space, wind tunnels, fatigue rigs, and for software, automated suites (see 34/31). The rule teams repeat: test as you fly, and fly as you tested. The senior capstone is the one place most students get to feel all of this at small scale. In the US, accreditation criteria from ABET require engineering programmes to include a culminating major design experience that draws on earlier coursework and works within multiple realistic constraints, which is why almost every programme has one. A typical capstone gives a team of four to six students a real or realistic customer, a budget of a few thousand dollars, and two semesters. The lessons are rarely the equations. They are these: The customer's first sentence is not the requirement; asking the second question is. The part you assumed would be easy is the one that is late. The interface you did not write down is the one that does not fit. Your design is only done when somebody else can build, test and run it from your documents. What a graduate just sees The hard part of a product is on the lines between its parts, so read the interfaces first. A requirement you cannot test is a wish; every requirement should say how it will be verified. Verification is "built it right"; validation is "built the right thing". They fail separately. Everything gets heavier, hotter and later than the first estimate, so hold margin and shrink it only with evidence. In aerospace a kilogram late in the stack costs several at the start, and a watt costs heat, area and mass. Redundancy that runs the same software fails the same way twice. Ask "how could this break?" before it does, and rank by consequence, not only by a product of scores. Disagreements about a design are usually disagreements about weights; make them explicit. Where it shows up next Machine Design (34/19) gives the parts and their loads; this is how they become a product. Mechanics of Materials (34/11) supplies the mass you can afford to remove. Control Systems (34/20): a controller is a subsystem whose interface is a sensor, an actuator and a rate. A wrong sign or a wrong unit there is the Mars Climate Orbiter story with a motor. Heat Transfer (34/16) and Thermodynamics (34/14) are the thermal budget. Vibrations (34/17) is why launch loads decide what survives. Software Engineering (34/31) is the other half of nearly every modern product; the Ariane story sits exactly between them. Computer Security (34/33) is the same discipline applied to an adversary. In the real world: spacecraft, aircraft, cars, medical devices and phones are built this way, and in companies the job title is often "systems engineer" or "integration and test engineer". If you want to go deeper 📘 Systems Engineering: Principles and Practice, Kossiakoff, Seymour, Flanigan and Biemer (Wiley, 3rd ed., 2020): the standard university textbook for the whole subject. 📘 Product Design and Development, Ulrich, Eppinger and Yang (McGraw-Hill, 7th ed., 2020): a common capstone text on turning a need into a product. 🌐 NASA Systems Engineering Handbook: free to download from NASA, and the best-known practical reference. 📘 To Engineer Is Human: The Role of Failure in Successful Design, Henry Petroski (1985): a gentle, readable book on why things fail and what designers learn. 🌐 MIT OpenCourseWare 16.842, Fundamentals of Systems Engineering (Prof. Olivier de Weck): free and complete, though a graduate-level course, so expect it to be dense. The takeaway A product is many subsystems agreeing. Systems engineering is the discipline of writing the agreements down (requirements and interfaces), proving them level by level (the V, verification and validation), spending scarce mass, power and heat honestly (budgets and margins), and hunting the ways it can break (FMEA). The big failures were not exotic physics. They were two neighbours with different assumptions. Try this: pick something in your home with two or more parts that must work together, such as a phone and its charger, or a door and its lock. Write down three things the one part assumes about the other. Which assumption would be most expensive if it were wrong? Sources & Further Reading Sources. Mars Climate Orbiter: NASA press release of 10 November 1999 on the Mishap Investigation Board's Phase I report; Wikipedia, "Mars Climate Orbiter". Ariane 5 flight 501: the Inquiry Board report of 19 July 1996 (chair Jacques-Louis Lions); Wikipedia, "Ariane flight V88". V-model and verification versus validation: Wikipedia, "V-model". FMEA history and the move to action priority: Wikipedia, "Failure mode and effects analysis". Dry-mass growth of 25 to 30 percent: NASA Systems Engineering Handbook SP-6105 (1995), p. 62, as quoted in a project document. Accreditation: ABET Criterion 5. Computed in node: rocket-equation ratio 3.577, radiator area 0.256 m², the mass-budget sums, trade-study scores and RPN products. 🌐 NASA Systems Engineering Handbook 🌐 MIT OpenCourseWare 16.842: Fundamentals of Systems Engineering 🌐 NASA: Mars Climate Orbiter Phase I report press release 🌐 Ariane 501 Inquiry Board report (mirror) 📘 To Engineer Is Human, Henry Petroski