Latest / Mechanical Engineering Made Simple / (#118) Introduction to Fluid Mechanics - Lesson 3
Transcript
- 0:00Welcome back, design junkies, to Mechanical Engineering Made
- 0:02Simple. If you enjoy the podcast, show
- 0:04your support with a like, a follow or a comment.
- 0:07Thanks for having me back. Oh, absolutely.
- 0:09So imagine for a second that you just, you know, turn off your
- 0:12kitchen faucet and the pipes behind your wall just suddenly
- 0:15explode. Like with the force of a hand
- 0:18grenade. Which is a terrifying thought,
- 0:19honestly. Right.
- 0:20Or picture yourself flying a jet so incredibly fast that the
- 0:26invisible thin air in front of you just instantaneously turns
- 0:30into a solid, impenetrable wall of force.
- 0:34A wall that can literally RIP metal apart.
- 0:36Yeah, that's the sound barrier for you.
- 0:38Exactly. So today we are bringing you
- 0:40into a very special master class.
- 0:42We are talking directly to you, the student of the physical
- 0:45world. And our mission today is
- 0:47ambitious, but I promise it is going to be incredibly
- 0:50rewarding, going to be a fun one.
- 0:51We are going to build a comprehensive mental model of
- 0:54fluid dynamics from the absolute ground up.
- 0:57We're talking compressible gases at supersonic speeds, then
- 1:02moving into the erratic nature of unsteady liquid flows right,
- 1:06and finally arriving at the cutting edge computer
- 1:09simulations that allow us to actually see and predict all
- 1:13this chaos. It really is a profound journey.
- 1:16You know, when you look at the history of how we understand
- 1:19fluids, which by the way means anything that flows liquids or
- 1:23gases, you realize that for a long time we were just kind of
- 1:26guessing. Just trial and error.
- 1:28Exactly. And we're going to trace the
- 1:30path of how we stopped guessing and started manipulating the
- 1:34very fabric of pressure and temperature to do things like
- 1:37we'll go to the moon. Yeah, you really can't
- 1:39appreciate why we need to build multi $1,000,000 supercomputers
- 1:43to simulate like a simple splash of water until you first
- 1:46understand the sheer mathematical headache of
- 1:49supersonic gas. Oh, it's a massive headache.
- 1:51I love that framing. So let's start at the very
- 1:53beginning. Let's talk about the fundamental
- 1:55nature of the stuff we were moving around.
- 1:57The fluids themselves. Right.
- 1:59When we talk about fluids in engineering, we're usually
- 2:02categorizing them based on one major physical property, which
- 2:05is how they handle being squeezed.
- 2:08I always thought of water and air as fundamentally different
- 2:10things, but mechanically they're both just fluids, right?
- 2:14That is the perfect starting point.
- 2:16Yes, both are fluids, but they sit on opposite sides of a very
- 2:19important line, which is compressibility.
- 2:22We define fluids broadly as either compressible or
- 2:25incompressible. Let's look at liquids first.
- 2:29Like water or hydraulic oil, for almost all engineering purposes
- 2:33liquids are treated as completely incompressible.
- 2:36So if I take like a plastic medical syringe, fill it
- 2:40entirely with water, put my thumb over the opening and push
- 2:43down on the plunger as hard as I can.
- 2:45The plunger doesn't move. It doesn't budget.
- 2:47At all, right? I mean, it barely yields on a
- 2:48microscopic level, but practically speaking it does not
- 2:51move. The density of that water, the
- 2:53amount of mass packed into that specific volume, remains
- 2:56essentially constant regardless of how much pressure your thumb
- 2:59is applying. Because the molecules are
- 3:01already packed close closely together.
- 3:03Exactly. Gases, however, are a completely
- 3:06different story. They are highly compressible
- 3:08because the molecules have vast amounts of empty space between
- 3:11them. So if you empty that same
- 3:13syringe, fill it with air, cap it and push.
- 3:16You can easily compress that volume of air down to a fraction
- 3:19of its original size. You're forcing those air
- 3:21molecules closer together, increasing the density.
- 3:24Like squeezing a massive volume of atmospheric air into a thick
- 3:29aluminum scuba tank. Precisely.
- 3:31But there is a massive caveat to this rule, and it's actually one
- 3:35of the most useful shortcuts in all the fluid mechanics.
- 3:38I love a good shortcut. What is it?
- 3:40Even though gas is highly compressible, we can actually
- 3:42treat it mathematically as if it were an incompressible liquid,
- 3:46provided the conditions are right.
- 3:47Wait, really? How can we just ignore the fact
- 3:49that air squishes? It all depends on how fast the
- 3:52air is moving relative to the speed of sound.
- 3:55If a gas is moving slowly, specifically below Mach 0.3,
- 4:00which is roughly 230 mph, the change in its density as it
- 4:04flows over an object is remarkably tiny.
- 4:06How tiny are we talking? It's actually less than 5%.
- 4:09Oh wow. And in the grand scheme of
- 4:11designing, say, an air conditioning duct, a 5% change
- 4:15in density just doesn't matter enough to complicate the math.
- 4:19Exactly. For practical everyday
- 4:21calculations, a 5% density change is totally negligible.
- 4:25So the air blowing out of the desk fan in your office behaves,
- 4:30mathematically speaking, almost exactly like water flowing
- 4:33through a wide pipe. That's wild.
- 4:35The equations are literally the same.
- 4:37They are. It is only when we start pushing
- 4:40gas to extreme high velocities or dealing with massive altitude
- 4:44and pressure changes that we cross a threshold.
- 4:46And suddenly that compressibility cannot be
- 4:49ignored. Right.
- 4:50The air stops acting like water and starts acting like a highly
- 4:52volatile spring. OK, so when we cross that
- 4:54threshold, when we start dealing with air moving at hundreds of
- 4:57miles an hour, we have to rely on foundational thermodynamics,
- 5:00right? We do.
- 5:01And the big one, the classic rule that governs all of this,
- 5:03is the perfect gas law. Yes, the absolute cornerstone of
- 5:07thermodynamics, the relationship where PV equals RT.
- 5:11Right. Pressure times volume equals the
- 5:14gas constant times temperature. Now I want to make sure I'm
- 5:17visualizing this correctly rather than just memorizing a
- 5:20formula. What is this law physically
- 5:22telling us about the error? It tells us that for a perfect
- 5:26gas, these core properties, pressure, volume and
- 5:31temperature, are inextricably LinkedIn this delicate balancing
- 5:34act. Let's break it down.
- 5:36Sure. Pressure is just the physical
- 5:38force of gas molecules bouncing off the walls of a container.
- 5:41Temperature is a measure of how fast those molecules are
- 5:43vibrating and zipping around, and volume is simply the space
- 5:47they have to move in. So if I take a sealed box of air
- 5:50and I crush it down, decrease in the volume, I am forcing the
- 5:54same number of molecules into a tighter space.
- 5:56Right. So they're going to hit the
- 5:57walls more often. Which means the pressure goes
- 5:59up. Exactly.
- 6:00But the wall also dictates that if you compress it and hold the
- 6:03pressure steady, somehow the temperature has to change.
- 6:06OK, how so? If you do work on a gas by
- 6:09compressing it quickly, you're adding energy to the system and
- 6:12the temperature spikes. This is why you know the bottom
- 6:15of a bicycle pump gets hot when you're vigorously feeling a
- 6:18tire. Oh, yeah, I've totally felt
- 6:19that. That makes total sense.
- 6:21Now this brings up something really fascinating about how
- 6:24gases actually store heat, Right.
- 6:26Yeah. The specific heats.
- 6:28Right. Depending on how you hold a gas
- 6:31while you heat it up, it behaves differently.
- 6:34There's specific heat at a constant volume and specific
- 6:37heat at a constant pressure. Why does the container matter so
- 6:40much? Let's use a thought experiment.
- 6:42Imagine you have a rigid heavy steel box filled with air.
- 6:47It is sealed completely shut. OK, so the volume cannot change.
- 6:51Exactly. Now you put a blowtorch to the
- 6:53bottom of that box. You are adding heat energy
- 6:56because the box can't expand. Every single ounce of that heat
- 6:59energy goes directly into raising the internal energy of
- 7:02the gas. It makes the molecules vibrate
- 7:04and bounce around faster and faster.
- 7:06Right, so the temperature rises rapidly.
- 7:08This is heating at a constant volume.
- 7:10We denote that specific heat as C sub V.
- 7:12OK, that is straightforward. Heat goes in, temperature goes
- 7:15up. But now imagine a different
- 7:17container. Imagine a hollow cylinder with a
- 7:20movable piston sliding inside it, kind of like the cylinder in
- 7:23your car's engine. OK, the gas is trapped under the
- 7:26piston. Yep, if you hold a steady weight
- 7:29on top of the piston, the pressure inside remains
- 7:31constant. Now apply the blowtorch again.
- 7:34Well, as the gas heats up, the molecules get energized, they
- 7:37hit the piston harder, and because the piston can move, the
- 7:40gas expands, it pushes the piston upward.
- 7:43Right. The gas is physically doing
- 7:45mechanical work to push that heavy piston against gravity.
- 7:49It's expending energy to expand. Oh, I see where this is going.
- 7:53Because it is spending some of that incoming blowtorch heat to
- 7:56do the physical work of lifting the piston, there is less energy
- 8:00leftover to simply raise the temperature of the gas itself.
- 8:02Oh wow. So if I want to raise the
- 8:04temperature of the gas in the cylinder by 10°, I have to apply
- 8:08the blowtorch for a much longer time than I did for the rigid
- 8:11steel box. You have to add more heat
- 8:14Exactly. The specific heat at a constant
- 8:17pressure, which is C sub P, is always higher than the specific
- 8:22heat at a constant volume C sub V.
- 8:25And I saw in the reading that the mathematical difference
- 8:27between those two heat values is really elegant.
- 8:30It is C sub P -. C sub V is exactly equal to R,
- 8:33the gas constant. That difference in heat energy
- 8:37is literally a direct measurement of the mechanical
- 8:40work the gas performed during its expansion.
- 8:42That is an amazing way to visualize the math.
- 8:44So we have pressure, volume, temperature and heat all dancing
- 8:48together, trading energy back and forth.
- 8:50The thermodynamic dance, Yeah. But how does this dance dictate
- 8:53something like the speed of sound?
- 8:54I mean, when I think of sound, I just think of noise.
- 8:56Someone talking, a dog barking. Physically, what is it?
- 8:59If we strip away biology and just look at the physics, a
- 9:02sound wave is nothing more than a microscopic disturbance in a
- 9:05fluid. It is a tiny temporary pressure
- 9:09discontinuity. Like a ripple, yeah.
- 9:11When you speak, your vocal cords push a layer of air molecules.
- 9:15Those molecules bump into their neighbors, transferring their
- 9:18momentum, and then bounce back. Those neighbors bump into the
- 9:20next layer, and so on. It's like a microscopic domino
- 9:23effect. Exactly.
- 9:24It is a longitudinal wave of compression propagating through
- 9:27the fluid. The actual air molecules don't
- 9:30travel from my mouth to your ear, they just wiggle back and
- 9:33forth in place. But the wave of pressure
- 9:35travels. Right, and the velocity at which
- 9:38that specific pressure wave propagates forward is what we
- 9:41call the speed of sound. And what governs how fast that
- 9:44domino effect can happen? Because I was looking at the
- 9:47equations for this and it totally upended my assumptions.
- 9:49What did you assume? I figured that the speed of
- 9:52sound would depend on how thick or pressurized the air is, like
- 9:55denser air faster sound, but it doesn't.
- 9:58It does not. The propagation velocity of a
- 10:01sound wave, usually denoted as A, is calculated by taking sqrt
- 10:06K * r * t. OK, let's break that down.
- 10:09K is a specific ratio of heats, R is the gas constant, and T is
- 10:14the absolute temperature. Now for a given gas, like the
- 10:18air we breathe, K&R are fixed numbers.
- 10:21They just don't change. Which means the only variable
- 10:23left in the equation is T temperature.
- 10:26That is the beautiful, isolating truth of it.
- 10:29The speed of sound in a perfect gas is strictly entirely
- 10:33proportional to the square root of its absolute temperature.
- 10:36Nothing else. Matters nothing else matters.
- 10:38Not the ambient pressure, not the density, just the
- 10:40temperature. Let me make sure I grasp this.
- 10:43If I am in a room at sea level and it's a comfortable 16°C, the
- 10:48speed of sound is roughly 340 meters per second.
- 10:51Correct. If I take a highly pressurized
- 10:53chamber at sea level, pump it up to 10 times atmospheric
- 10:57pressure, but keep the temperature exactly at 16°C, the
- 11:01speed of sound is exactly the same.
- 11:03Exactly the same, because while you increase the pressure, the
- 11:06density increase proportionally, and those two effects cancel
- 11:09each other out perfectly. The only thing that determines
- 11:12how fast those molecules can pass the bump of a sound wave to
- 11:15their neighbor is how fast they're already vibrating, and
- 11:18temperature is just a measure of that vibration.
- 11:20Hotter molecules are zipping around faster, so they transmit
- 11:24the wave faster. Colder molecules are sluggish,
- 11:26so the wave travels slower. So when a commercial airliner
- 11:29climbs up to 30,000 feet where the air pressure is incredibly
- 11:33thin, it's not the lack of air pressure that changes the speed
- 11:36of sound. No, it's the fact that it is
- 11:39freezing cold. Up there.
- 11:41In the freezing upper atmosphere, the speed of sound
- 11:43drops significantly. This is why engineers don't
- 11:46typically measure high speed flight in purely mph.
- 11:49Right, they use the Mach number named after Ernst Mock.
- 11:52Exactly, the Mach number is incredibly simple.
- 11:55It's a ratio m = U / a. You take the objects actual
- 12:01velocity U and you divide it by the local speed of sound a in
- 12:06the exact environment it is flying through.
- 12:08So if the object's velocity is equal to the speed of sound, the
- 12:11Mach number is exactly 1. So a plane flying at a constant
- 12:14500 mph on the deck right above the ocean might be flying at
- 12:18Mach 0.6 because the air is warm and sound travels fast, right?
- 12:24But if that exact same plane maintains 500 mph and climbs to
- 12:2840,000 feet where it's 30 below 0.
- 12:31Because the local speed of sound has dropped due to the cold,
- 12:34that same 500 mile per hour speed might now be Mach 0.8.
- 12:40The plane is much closer to breaking the sound barrier, even
- 12:43though it hasn't actually accelerated.
- 12:45OK, this brings us to a phenomenon that is absolutely
- 12:47crucial to understand before we can talk about building jet
- 12:51engines or rockets. The sound barrier itself.
- 12:53Right, What physically happens to the air as an object
- 12:56accelerates toward Mach one I want to visualize.
- 12:58This Let's use the Doppler effect to visualize it.
- 13:01Imagine a stationary object like a speaker resting on the ground,
- 13:04emitting A continuous tone. OK, it is radiating sound waves
- 13:08outward in perfect concentric spheres, like dropping a Pebble
- 13:12in a perfectly still pond. The ripple spread out evenly in
- 13:16all directions. I picture the ripples.
- 13:18Now put that speaker on a cart and start driving it forward,
- 13:22but drive it slower than the speed of sound so the Mach
- 13:25number is less than one. The speaker emits a wave, moves
- 13:28forward a bit, emits another wave, moves forward and so on.
- 13:31Exactly. It's essentially chasing its own
- 13:33sound waves. The waves in front of the cart
- 13:36get brunched up, squeeze closer together.
- 13:39The waves behind the card get stretched out.
- 13:41That bunching up is the Doppler effect.
- 13:44Because the waves in front are squeezed tighter, the frequency
- 13:47increases. Right, which is why when an
- 13:50ambulance is driving towards you, the siren sounds higher
- 13:52pitched. As it passes you and drives
- 13:55away, the waves stretch out and the pitch drops noticeably.
- 13:59I have a specific way I've always liked to think about
- 14:01this, to ground it in reality. I think of it like a traffic jam
- 14:05on a highway. I like this, walk me through it.
- 14:07Imagine those sound waves radiating forward from the
- 14:10object are like the red brake lights of cars.
- 14:13If a car hits its brakes, the red lights come on, the driver
- 14:17behind them sees it and hits their brakes, and their lights
- 14:19come on. It creates a chain reaction, a
- 14:22warning signal traveling backward up the highway, telling
- 14:25all the cars, hey, there's an obstacle ahead, slow down.
- 14:28Adjust your path. That's a great visual.
- 14:30When an airplane is flying subsonic, slower than sound, the
- 14:34pressure waves it emits act exactly like those brake lights.
- 14:37They travel forward faster than the plane, warning the
- 14:40undisturbed air molecules miles ahead.
- 14:43They are essentially communicating with the air.
- 14:46They are sending a physical signal that says a plane is
- 14:48coming. And because the air gets that
- 14:50advance note, notice the molecules have time to smoothly
- 14:53part and flow gracefully around the shape of the wings and the
- 14:56fuselage. The air is prepared, but now
- 14:59let's accelerate our cart. What happens when the Mach
- 15:02number equals exactly 1? The object is now moving at the
- 15:05exact propagation velocity of its own sound.
- 15:07Waves. It's keeping pace perfectly, so
- 15:10it emits a wave, but it doesn't fall behind.
- 15:12It rides right on top of it. Yes, the waves can no longer
- 15:16escape out in front. Every single pressure wave the
- 15:19object emits forward stacks up perfectly on top of the previous
- 15:22one, right at the leading edge of the object.
- 15:24It creates A formidable, intense wall of compressed air.
- 15:28Returning to the traffic analogy, if you are driving your
- 15:31car faster than the brake light warning can physically travel up
- 15:35the highway, you outrun your own warning signals.
- 15:38You slam into the stop traffic ahead of you without giving them
- 15:41any prior notice to move out of the way.
- 15:43You violently crash into the undisturbed air, and if you push
- 15:47even faster, if your Mach number goes greater than one, you
- 15:49completely leave your waves behind.
- 15:52You are constantly plunging into completely undisturbed,
- 15:55completely unwarned air. The air molecules don't have
- 15:58time to gently part, they are violently shoved aside.
- 16:00And this abrupt, violent collision is what creates a
- 16:03shockwave. Exactly.
- 16:05And what happens to all the sound the object is making?
- 16:07Well, it can't go forward. Right, because you are out
- 16:10running your own waves. All those overlapping pressure
- 16:13spheres form a strict geometric shape behind you.
- 16:16It creates a cone. We call this the mock cone.
- 16:19The mock cone. The tip of the cone is pinned to
- 16:21the nose of the aircraft, and the cone flares out behind it.
- 16:25All of the sound, all of the pressure disturbances are
- 16:28completely confined to the inside volume of this cone.
- 16:31That is wild. So if I am standing in a quiet
- 16:34desert and a supersonic fighter jet flies directly over my head
- 16:38at Mach 2, I literally will not hear it coming.
- 16:41You will hear absolute silence as it approaches, and even as it
- 16:44is directly overhead, it isn't until the edge of that invisible
- 16:48mock cone physically drags across your body on the ground
- 16:52that the sound finally reaches you.
- 16:53And. Because all those compressed
- 16:55sound waves are piled up on the edge of that cone.
- 16:57When it sweeps over you, you don't just hear the engines, you
- 16:59experience a massive explosive Sonic boom.
- 17:03Exactly. And the faster the jet flies,
- 17:05the sharper the angle of that cone becomes, right?
- 17:07It sweeps back tighter against the jet.
- 17:09Right, there's math for that too.
- 17:11The angle of the Mach cone, which we call alpha, is defined
- 17:14mathematically by sine. Alpha equals 1 / m so at exactly
- 17:19Mach 1, the cone is flat. It's a wall.
- 17:22But at Mach 2 it sweeps back to a 30° angle.
- 17:25And at Mach 5, it's an incredibly tight needle shape.
- 17:28So let's pivot here. If pushing an object through the
- 17:32air at high speeds creates these massive, violent pressure
- 17:35differences in shock waves, how on earth do aerospace engineers
- 17:41design things like rockets? That's the $1,000,000 question.
- 17:45How do you take hot chaotic gas and intentionally force it past
- 17:49that Sonic limit without it just blowing the engine apart?
- 17:52That is one of the great triumphs of mechanical
- 17:54engineering. To understand it, we have to
- 17:56leave the comfort of standard, everyday fluid dynamics behind.
- 18:00We have to look at how pipes and nozzles handle supersonic flow,
- 18:03because the rules are entirely reversed.
- 18:06OK, let's explore this. We are moving into the bizarre,
- 18:08paradoxical geometry of supersonic nozzles.
- 18:11To set the baseline, we have to imagine a theoretically perfect
- 18:15scenario. We assume the flow of gas
- 18:17through our nozzle is isentropic.
- 18:19What does that mean practically? It means we are assuming a
- 18:22reversible adiabatic process. Adiabatic means no heat is being
- 18:26transferred in or out of the pipe walls.
- 18:28The gas isn't losing heat to the metal.
- 18:30And reversible. Reversible means there is 0
- 18:32friction. The entropy of the system, the
- 18:35fundamental thermodynamic measure of disorder and wasted
- 18:37energy remains completely constant.
- 18:40The change in entropy is 0 OK. A magical, frictionless, perfect
- 18:45pipe. In this perfect system, we have
- 18:47an unbreakable rule called the continuity equation.
- 18:50It simply states that the mass flow rate must remain constant.
- 18:53Matter cannot magically appear or disappear.
- 18:56So whatever amount of mass of gas enters the front of the
- 18:58pipe, that exact same amount of mass must exit the back of the
- 19:02pipe every single second. Exactly.
- 19:04The mass flow rate M is calculated by multiplying 3
- 19:07things, the density of the gas rho, the velocity it's moving U,
- 19:11and the cross-sectional area of the pipe.
- 19:13Rho times U * a right? If you change one of those 3
- 19:16variables, the other two must automatically adjust to keep the
- 19:19total mass flow rate identical. OK, I can visualize this using
- 19:22normal everyday subsonic flow, like water coming out of a
- 19:25garden hose. Let's use that.
- 19:26The water is moving relatively slowly, so the Mach number is
- 19:29far less than one. It behaves just how your
- 19:32intuition expects. If I have water flowing out of
- 19:35the hose and I want it to spray faster to reach the other side
- 19:37of the yard, I put my thumb partially over the opening.
- 19:40You are decreasing the cross-sectional area of the
- 19:43pipe. Right, and because water is
- 19:45incompressible, its density doesn't change.
- 19:47So to force the exact same amount of mass through a much
- 19:50smaller hole every second, the water has no choice but to
- 19:53dramatically increase its velocity.
- 19:55Perfect for subsonic flow. A decrease in area causes an
- 20:00increase in velocity. It is highly intuitive, but here
- 20:05is the massive counterintuitive twist.
- 20:08If the gas inside that pipe is already moving at supersonic
- 20:11speeds. If the Mach number is greater
- 20:13than one, the physical rules of the universe flip completely
- 20:17upside down in a supersonic flow.
- 20:20If you put your thumb over the end of the hose.
- 20:22If you decrease the area of the pipe.
- 20:24The velocity of the gas actually decreases, it slows down.
- 20:27Let me pause you there. That makes 0 logical sense.
- 20:29Why would making the hole smaller slowest speeding gas
- 20:32down? It all comes back to our first
- 20:34topic, extreme compressibility. When a gas is moving that
- 20:38incredibly fast, faster than the speed of sound, and you suddenly
- 20:42force it through a tighter bottleneck, the gas doesn't
- 20:44smoothly speed up. It compresses.
- 20:46Exactly. The intense collision with the
- 20:48narrowing walls causes the gas to compress violently.
- 20:52The molecules just slam into each other.
- 20:54So the pressure and the density absolutely skyrocket.
- 20:57In fact, the math dictates that in a supersonic flow, the
- 21:00density of the gas increases at a much faster rate than the
- 21:03physical area of the pipe decreases.
- 21:05Oh wait, let me think about the continuity equation.
- 21:08Mass flow equals density times velocity times area, right?
- 21:13If the area shrinks a little bit, but the density absolutely
- 21:16explodes to a huge number, then to keep the overall mass flow
- 21:20equation balanced, the velocity variable is mathematically
- 21:24forced to drop. You got it.
- 21:26The immense spike in density essentially chokes the speed.
- 21:29Squeezing A supersonic flow slows it down.
- 21:31Which means if shrinking the pipe slows A supersonic gas
- 21:35down, to make a supersonic gas go even faster, you have to do
- 21:39the exact opposite. That is the astonishing reality.
- 21:42To accelerate supersonic flow, to make it go faster, the pipe
- 21:46must actually widen. You have to increase the
- 21:49cross-sectional area. That is so deeply
- 21:51counterintuitive to my brain. If I take my thumb off the
- 21:54garden hose, the water slows down to a trickle.
- 21:57But for a supersonic gas, making the exit wider is the
- 22:01accelerator pedal. It's wild, isn't it?
- 22:04How did engineers ever figure out how to harness this?
- 22:07It was figured out in the late 19th century by an engineer
- 22:10named Gustav de Laval working on steam turbines.
- 22:13He invented what we now call the convergent divergent nozzle.
- 22:17The nozzle. Right.
- 22:19It is arguably one of the most elegant, pivotal pieces of
- 22:21engineering in human history. It is the iconic hourglass bell
- 22:26shaped noggle. You see on the back of it every
- 22:28rocket engine that has ever gone to space.
- 22:30OK, I want to walk through exactly how this bell shape
- 22:32works from front to back, step by step, because the geometry
- 22:35feels like a magistrate. We start at the top right at the
- 22:38combustion Chamber of the rocket.
- 22:39Inside that chamber you have incredibly hot, highly
- 22:43pressurized gas from burning rocket fuel.
- 22:47But crucially, inside the massive chamber, the gas is
- 22:51moving relatively slowly. It is subsonic.
- 22:54OK, so as it flows out of the chamber, it enters the first
- 22:57half of the hourglass shape. The pipe starts to narrow.
- 23:00This is the convergence section. And because the gas is currently
- 23:03subsonic, it behaves normally just like your garden hose.
- 23:07The narrowing area causes the gas to accelerate.
- 23:09It gets faster and faster and faster as the pipe squeezes
- 23:13tighter. Until it hits the absolute
- 23:15narrowest point in the middle of the hourglass Engineers call
- 23:17this the throat. Right.
- 23:19And if the pressure inside the combustion chamber is high
- 23:22enough, something incredible happens precisely at that
- 23:24microscopic slice of the. Throat.
- 23:26An absolute physical limit is reached, yes.
- 23:29The gas accelerates to exactly Mach one, the local speed of
- 23:32sound. We call this the critical
- 23:34velocity. And when a gas hits Mach 1,
- 23:36specifically at the narrowest throat of a pipe, it triggers a
- 23:39condition engineers call choked flow.
- 23:42What does it mean for a pipe to be choked?
- 23:44It means the flow has mathematically maxed out.
- 23:47You cannot push any more mass through that hole at the throat.
- 23:51The pressure of the gas drops to a very specific fraction of the
- 23:55original chamber pressure. We call this the critical
- 23:57pressure or P star, right? Yes, for a standard gas like
- 24:00air, that critical pressure ratio P star divided by P
- 24:03naughty is a fixed constant 0.528.
- 24:08So the pressure at the throat is always just over half of
- 24:11whatever the pressure is in the main chamber.
- 24:13Exactly. Once you hit that critical
- 24:15pressure ratio and the gas achieves Mach one at the throat,
- 24:18the nozzle is completely choked. And here here's the crazy part.
- 24:22Once it's choked, you could alter the environment completely
- 24:25and the mass flow rate won't change.
- 24:27What do you mean? You can take the exit of that
- 24:28nozzle and expose it to the vacuum of deep space.
- 24:31You could drop the outside pressure to absolute zero and
- 24:34the amount of gas moving through that throat will not increase by
- 24:37a single solitary gram. I am struggling to understand
- 24:40why if I hook a vacuum cleaner up to a pipe, it sucks more air
- 24:45through it. Why doesn't the immense infinite
- 24:48vacuum of space physically suck the gas out of the rocket
- 24:51faster? It goes back to how fluid knows
- 24:54what the environment is doing, the information that there is a
- 24:57vacuum downstream, that suction is physically transmitted
- 25:02backward through the gas in the form of a pressure waste, a
- 25:06sound wave. OK, so the vacuum tries to send
- 25:08a signal upstream saying hey, it's empty down here, flow
- 25:11faster, yes. But that signal can only travel
- 25:14through the gas at the local speed of sound.
- 25:16Meanwhile, the gas itself is physically flowing downstream
- 25:19through the throat at the exact speed of sound.
- 25:22Oh my God, it's like a salmon trying to swim upstream against
- 25:24a waterfall. If the waterfall is falling at
- 25:2720 miles an hour and the salmon can only swim at 20 miles an
- 25:29hour, the salmon stays perfectly stationary relative to the
- 25:33shore. It is exactly like that.
- 25:35The pressure signal trying to report the vacuum conditions can
- 25:38never physically swim past the throat.
- 25:40It is swept downstream as fast as that tries to propagate
- 25:43upstream. That's brilliant.
- 25:45Therefore, the combustion chamber is completely blindly
- 25:48isolated from whatever is happening past the throat.
- 25:51The chamber has no idea the vacuum is there.
- 25:54The throat limits the flow absolutely.
- 25:56That is a phenomenal piece of physics.
- 25:58But here's the problem. A rocket needs its exhaust to
- 26:02leave the engine much, much faster than Mach 1 to generate
- 26:05enough thrust to overcome gravity.
- 26:07If the throat limits it to Mach one, how do we get more speed?
- 26:11This is where Gustav de la Valle's genius comes in.
- 26:14Immediately after the narrow throat, the nozzle begins to
- 26:17gently flare outward. The pipe lightens.
- 26:20This is the divergent section, the bottom half of the
- 26:22hourglass, the big bell shape you see on a space shuttle.
- 26:25And because the gas is already traveling at exactly Mach 1 when
- 26:28it enters this widening section, our paradoxical supersonic rule
- 26:32kicks. In exactly the gas is now
- 26:33supersonic, so the gradually widening area of the bell
- 26:37doesn't slow the gas down like it would in a garden hose.
- 26:39Instead, the widening area allows the gas to rapidly,
- 26:42violently expand. The density plummets, and to
- 26:46satisfy the continuity equation, as the density drops, the
- 26:49velocity is forced to absolutely skyrocket.
- 26:51The gas expands and accelerates to Mach 2, Mach 3, Mach 10,
- 26:55shooting out the back and generating massive world shaking
- 26:59thrust. When you step back and look at
- 27:01the bigger picture, you have to admire the profound elegance of
- 27:04this. We are not using complex moving
- 27:07mechanical pumps to force gas past the sound barrier.
- 27:11There are no moving parts in the nozzle.
- 27:13Just geometry. We are simply manipulating the
- 27:16stationary physical geometry of a piece of metal to trick a hot
- 27:19gas into breaking the laws of intuition.
- 27:22We are exploiting the mathematical quirks of density
- 27:25and pressure ratios to send humanity into orbit.
- 27:28It is absolute design perfection, but nature usually
- 27:31doesn't like it when we push things to the absolute limit.
- 27:34No it doesn't. What happens if our perfect
- 27:36geometry meets the messy real world?
- 27:38Like a rocket engine is designed to expand gas to a specific
- 27:42pressure, what if the atmospheric pressure outside the
- 27:44rocket doesn't perfectly match what the nozzle is expecting?
- 27:47That leads to highly violent consequences.
- 27:49The gas hits a wall and this is where we have to explore what
- 27:52happens when supersonic flows fail and how shock waves are
- 27:55formed inside machinery. Let's dive into that.
- 27:58We've got this beautifully smooth supersonic acceleration
- 28:01happening in the divergent bell of our nozzle.
- 28:04But what if the back pressure, the actual pressure of the
- 28:07atmosphere pushing back against the exit of the engine, is too
- 28:09high? Say the rocket is still on the
- 28:12launchpad at sea level where the air is thick.
- 28:14If the back pressure from the atmosphere is too high, the gas
- 28:17inside the nozzle overextends. It drops its pressure too low
- 28:21during its acceleration, and suddenly it finds itself pushing
- 28:25against the wall of thick atmospheric air.
- 28:27So it hits a literal wall of air.
- 28:29Yeah, the supersonic velocity simply cannot be maintained
- 28:32against that heavy resistance. So it has to slow down.
- 28:35Yes, but fluid dynamics dictates that a supersonic flow cannot
- 28:40just gently, smoothly decelerate back down to subsonic speeds.
- 28:44It can't just ease off the gas pedal, it has to slam on the
- 28:47brakes. Oh boy.
- 28:48And the way a fluid slams on the brakes is by spontaneously
- 28:51generating a shockwave right inside the divergent bell of the
- 28:55nozzle. Wait, it generates a Sonic boom
- 28:57inside the metal engine itself? It does.
- 28:59The gas is cruising along at Mach 2.
- 29:02It hits this invisible threshold of atmospheric resistance and
- 29:05instantaneously, across a microscopic distance, it drops
- 29:09to subsonic speeds. That sounds incredibly
- 29:11destructive. It can be.
- 29:13As a rocket climbs through the atmosphere and the outside air
- 29:15pressure changes, this internal shockwave will actually migrate
- 29:19up and down the inside of the bell.
- 29:21Let's clearly define what this shockwave is physically, because
- 29:24it's not just a loud noise. Earlier we said a normal sound
- 29:28wave is a tiny gentle pressure disturbance.
- 29:32What is the shockwave? The shockwave is a macroscopic
- 29:35pressure discontinuity. It is a massive, incredibly
- 29:38violent step change in the state of the gas.
- 29:40Imagine a wall roughly a few micrometers thick.
- 29:43OK, picture a wall. On the upstream side of that
- 29:45invisible wall, the gas is moving at supersonic speed with
- 29:48low pressure and relatively low temperature.
- 29:50The instant the gas molecules pass through that micrometer
- 29:53thin wall, the velocity plummets to subsonic, the pressure spikes
- 29:57violently, and the temperature absolutely skyrockets.
- 29:59It's like the gas is running full speed into a brick wall of
- 30:02its own making. And does a shockwave travel at
- 30:05the speed of sound like a normal noise?
- 30:08No, this is another key difference.
- 30:09A normal sound wave always travels at exactly the local
- 30:13seed of sound, regardless of how loud it is.
- 30:15But a shockwave is different. A shockwave actually travels
- 30:18faster the larger the pressure rises across it.
- 30:22So bigger explosion creates a faster shockwave.
- 30:25Exactly. If a bomb goes off, the initial
- 30:27shockwave rips outward much, much faster than Mach one, but
- 30:31as it expands, it loses energy, the pressure difference drops,
- 30:35the wave slows down, and eventually it decays into a
- 30:38regular Sonic sound wave that just sounds like a distant boom.
- 30:41I was looking at the math that governs this brutal transition.
- 30:45There is a set of relationships called the rank and Hugoniot
- 30:48equations. Now, I am not going to pretend I
- 30:49can solve nonlinear differential equations in my head, but
- 30:52conceptually, what do these equations fundamentally prove
- 30:55about how shockwaves behave? The rank and Hugoniot equations
- 30:59are the definitive map of a shockwave.
- 31:02They precisely correlate exactly how pressure, density and
- 31:05temperature on the front side of the shock relate to the
- 31:07backside, right? But the most profound revelation
- 31:10hidden inside these equations comes when you look at the
- 31:13entropy calculations. Entropy being the measure of
- 31:16disorder or wasted heat in a system.
- 31:18Right. The math lays out an absolute
- 31:20ironclad rule of physics. A shockwave can only develop
- 31:24when the upstream flow of the gas is supersonic.
- 31:27It is physically mathematically impossible to generate a
- 31:30shockwave in a subsonic flow. Why is it impossible?
- 31:33Because the equations show that if a shockwave were to form in a
- 31:36subsonic flow, it would cause the total entropy of the system
- 31:39to decrease. And according to the second law
- 31:42of thermodynamics, in an isolated system entropy can only
- 31:45increase or remain constant, it can never go backward.
- 31:48The universe forbids order from spontaneously forming from chaos
- 31:52without adding outside energy. Therefore, the universe forbids
- 31:55subsonic shock waves. That is mind bending.
- 31:58The physical geometry of our universe requires you to be
- 32:02breaking the sound barrier to create this phenomenon.
- 32:05Now these shock waves aren't just flat walls, right?
- 32:08They take on different shapes depending on what they are
- 32:11crashing into. Yes, the reference is
- 32:13distinguished between normal and oblique shock waves.
- 32:18Let's talk about that. A normal shockwave is the
- 32:20simplest type. It stands at exactly 90°
- 32:23perpendicular to the flow direction.
- 32:25It's a flat wall. You typically only see normal
- 32:28shocks contain inside pipes or hovering right in the throat of
- 32:31a jet engine intake. OK, what about an oblique?
- 32:34Shock An oblique shockwave happens at an angle.
- 32:37Like the wake behind a boat, but in the air.
- 32:39That's a great visual. Imagine a supersonic fighter
- 32:42jet. The nose of the jet is usually
- 32:44shaped like a sharp wedge or a needle.
- 32:47As that sharp point slices into the supersonic air, it deflects
- 32:50the air at an angle. It can't compress straight back.
- 32:53Exactly. The air compresses diagonally,
- 32:55creating an oblique shockwave that sweeps back from the nose.
- 32:58The shape of the object doing the piercing really matters
- 33:01here, doesn't it? If the wedge angle of The Jets
- 33:03nose is sharp enough, the shockwave stays perfectly
- 33:06attached to the very tip of the metal.
- 33:08It slices the air clean. It does.
- 33:11But what if you try to push a blunt fat object to supersonic
- 33:15speeds? If the physical angle of the
- 33:17object is too blunt. If the metal curves away too
- 33:20steeply, the air simply cannot physically turn fast enough to
- 33:24get out of the way. The mathematical limit is
- 33:26exceeded. What happens then?
- 33:28When that happens, the shockwave actually detaches from the
- 33:31vehicle entirely. It detaches.
- 33:32Where does it go? It stands off in front of the
- 33:34object. It's called a bow shock, where a
- 33:36detached shockwave. It hovers a few inches or
- 33:39sometimes feet directly in front of the blunt nose, moving along
- 33:43with the craft. Wow.
- 33:44This creates a massive, violent zone of highly pressurized,
- 33:47extremely hot subsonic air trapped right between the
- 33:51shockwave and the front of the vehicle.
- 33:53Which creates incredible aerodynamic drag.
- 33:55This is why supersonic planes like the Concorde or fighter
- 33:58jets look like sharp needles. They want to keep the shockwave
- 34:01cleanly attached to minimize drag.
- 34:04Yes, but on the other hand, look at the Apollo space capsules
- 34:07returning from the moon. They hit the atmosphere at Mach
- 34:0925 and their heat Shields were completely flat and blunt.
- 34:13And that was a deliberate genius engineering choice.
- 34:17Because they wanted to stop. Right, The Apollo capsule wanted
- 34:20to slow down. They wanted massive drag.
- 34:23But more importantly, by using a blunt shape to intentionally
- 34:26create a detached bow shock, that hovering layer of
- 34:30incredibly hot compressed air acted like a thermal buffer.
- 34:34Oh, that's smart. The most intense heat of re
- 34:36entry was concentrated in the shock wave standing a foot in
- 34:39front of the capsule rather than directly touching the metal
- 34:42surface of the heat shield. They used a wall of invisible
- 34:45compressed air to shield themselves from burning up.
- 34:48That is brilliant. OK, so we've established that a
- 34:50wedge or a concave curve compresses supersonic air and
- 34:55creates a shockwave. But what happens if a supersonic
- 34:58flow goes around a convex curve? What do you mean?
- 35:01Like what if the air is flowing fast over a flat surface and the
- 35:04surface suddenly drops away like a ledge?
- 35:06Does it create a reverse shockwave?
- 35:08It creates something entirely different and quite beautiful,
- 35:12called a Prantlemeyer Expansion fan.
- 35:14An expansion fan. Remember, a shockwave is a
- 35:16sudden violent compression. But when supersonic flow goes
- 35:20around a convex corner when it has more room to expand, it
- 35:24doesn't violently snap. Instead, it smoothly spreads
- 35:28out. How does a fan work?
- 35:30From the very point of the corner, the flow generates a fan
- 35:33made of countless infinitely thin pressure waves called Mach
- 35:37lines radiating outward. OK, so it curves around.
- 35:40As the supersonic air turns around the curve and passes
- 35:43through this fan, it's smoothly and continuously accelerates,
- 35:46expands and drops in pressure. There is no violent shock, no
- 35:50sudden wall, no loss of total energy, just a beautiful
- 35:54continuous acceleration around the bend.
- 35:56Wow. OK.
- 35:57So far we've been talking about these flows in a highly
- 35:59theoretical way. We assumed earlier that the flow
- 36:02in our nozzle was isentropic, perfect, frictionless and
- 36:06retaining all its heat. Which is never true in the real
- 36:08world. Exactly.
- 36:09Real pipes are made of rough metal that causes friction.
- 36:12Real jet engines have intense fire inside them.
- 36:15How do we model reality? That is where engineers rely on
- 36:192 highly specific models, Fanoflow and really flow.
- 36:23Let's start with Fano. Fanoflow is the mathematical
- 36:26model for a pipe where the flow is still adiabatic, meaning we
- 36:31assume no heat is passing through the walls, but we
- 36:34introduced the irreversible reality of friction.
- 36:37The gas is physically rubbing against the rough interior of
- 36:40the pipe, creating drag. OK, so Fano flow is the math of
- 36:43friction. What is Rayleigh Flow?
- 36:45Rayleigh flow approaches the problem from the opposite side.
- 36:48Rayleigh flow assumes the pipe is perfectly smooth and
- 36:50frictionless, but we add intense, deliberate heat
- 36:53transfer. Like in an engine.
- 36:55Yes, this is the model used to understand what happens inside
- 36:59the combustion Chamber of a jet engine, where burning jet fuel
- 37:02is constantly dumping massive amounts of thermal energy
- 37:05directly into the flowing gas. Friction and heat and the most
- 37:08mind blowing thing about these two models is what happens when
- 37:11you plot them on a graph. Oh yes, the entropy diagram.
- 37:14If you chart both Fanno flow and Rayleigh flow on a graph
- 37:18measuring entropy, they both do something completely bizarre.
- 37:21They do. Whether you were dealing with
- 37:23subsonic flow or supersonic flow, whether you are applying
- 37:26friction or applying heat, both the Fanno line and the Rayleigh
- 37:30line naturally inevitably curve and converge toward exactly Mach
- 37:351. The speed of sound, yes.
- 37:38If you have a subsonic flow in a pipe and you add friction, the
- 37:41friction causes the pressure to drop and the flow to accelerate
- 37:44toward Mach. One, and if it's supersonic.
- 37:46If the flow is supersonic, friction causes it to decelerate
- 37:49toward Mach 1. The same is true for heat.
- 37:52In Rayleigh flow, adding heat drives the gas toward the speed
- 37:55of sound. Mach 1 represents the exact
- 37:58state of maximum entropy for that specific pipe.
- 38:01It is like the thermodynamic destiny of the system.
- 38:04No matter what you do, the universe is trying to drive the
- 38:06flow toward the speed of sound, and once it hits that bottleneck
- 38:09of maximum entropy, it physically cannot go any further
- 38:12without changing the pipe itself.
- 38:14It is a profound bottleneck in nature.
- 38:16The physics simply converge. I want to pivot drastically here
- 38:20because this raises a major question for me.
- 38:22We have spent this entire time talking about gases.
- 38:25Compressible, squishy, high speed gases.
- 38:27But what happens when we look at liquids?
- 38:29We established the very beginning that liquids, like
- 38:32Water Act completely incompressible.
- 38:35They don't stretch, they don't squeeze.
- 38:36They are stubborn. So what happens if a heavy,
- 38:39stubborn, incompressible liquid suddenly undergoes an extreme
- 38:42pressure change? When you subject an
- 38:45incompressible liquid to extreme dynamic forces, the liquid stops
- 38:49acting so polite it becomes highly destructive.
- 38:53This brings us to the chaotic, often violent world of unsteady
- 38:57flow and a phenomenon known as the water hammer.
- 39:00Let's get into it. We are shifting our perspective
- 39:03from steady state flow, which is smooth, constant predictable
- 39:06movement, to unsteady flow where velocities and pressures are
- 39:10rapidly changing over time. Where do engineers actually deal
- 39:13with this in the real world? Unsteady flow is everywhere.
- 39:16It is the core of any system that transmits power or provides
- 39:19automatic precise control using liquids hydraulics.
- 39:23Like construction equipment. The heavy machinery on a
- 39:25construction site lifting tons of dirt, the flight control
- 39:28surfaces, moving the flaps on a commercial airliner, even the
- 39:31hydraulic brakes in your car. All of these rely on unsteady
- 39:34fluid dynamics. The foundation of this usually
- 39:36starts with a very simple, almost peaceful thought
- 39:39experiment. Imagine a clear plastic tube
- 39:43bent into AU shape and partially filled with liquid.
- 39:47Right, a simple U2 manometer. The water rests evenly in both
- 39:51vertical columns, perfectly balanced by gravity.
- 39:54Now imagine you blow gently into one side of the tube, pushing
- 39:58the liquid column down on the left, which forces the liquid on
- 40:01the right side to rise. And then you let go.
- 40:03You stop blowing. What happens?
- 40:05The liquid column immediately falls back down, but its
- 40:08momentum carries it past the center resting point.
- 40:11It pushes the left side back up. Exactly.
- 40:13The liquid will oscillate back and forth, sloshing from side to
- 40:17side, exactly like a pendulum swinging on a grandfather clock.
- 40:20And what's amazing is that you can calculate the exact timing
- 40:23of that slosh using the same physics as a pendulum.
- 40:27The cycle time t = 2 * π * sqrt l / 2 * G It is pure harmonic
- 40:35motion. It is, however, unlike a perfect
- 40:37pendulum in a vacuum, the water doesn't slosh forever.
- 40:40The liquid is rubbing against the walls of the plastic tube.
- 40:43Right friction. This laminar frictional
- 40:45resistance steals a tiny bit of kinetic energy with every swing.
- 40:49It acts as a dampener. Eventually the oscillations get
- 40:52smaller and smaller until the fluid settles perfectly back to
- 40:55rest. OK, so that's a gentle
- 40:57oscillation, but let's scale that up to industrial sizes.
- 40:59Let's look at what happens in a straight pipe.
- 41:01Imagine a massively long straight pipe connected to a
- 41:05huge water reservoir. There's a closed valve at the
- 41:08end. When you suddenly crank that
- 41:10valve open, the water doesn't just instantly blast out at its
- 41:14absolute maximum speed. You know it takes time.
- 41:17In the very first fraction of a second after opening the valve,
- 41:21the entire head of pressure, the massive weight of the water
- 41:25sitting in the reservoir, is dedicated entirely to overcoming
- 41:28inertia. Overcoming inertia.
- 41:30It is applying force to accelerate a very heavy
- 41:33stationary mass of water residing in the pipe.
- 41:35So it starts accelerating, yes. But as the water speeds up, the
- 41:39dynamic changes. Fluid friction against the inner
- 41:41walls of the pipe kicks in. The faster the water goes, the
- 41:44harder the friction fights back, stealing energy from the driving
- 41:47pressure. It creates this smooth
- 41:50mathematical curve on a graph, a logarithmic curve.
- 41:53Yes, extremely fast acceleration at first.
- 41:57When the water is slow and friction is low, then the
- 41:59acceleration tapers off gently. It curves until the driving
- 42:03pressure from the reservoir perfectly balances the friction
- 42:06loss of the pipe and the water achieves a final steady
- 42:09velocity. OK, so opening a valve is
- 42:12mathematically quite elegant. It is a smooth transition to a
- 42:15steady state. Closing a valve, however, is an
- 42:18absolute nightmare. The water hammer.
- 42:20The water hammer. I think almost everyone has
- 42:23experienced this in an old house or an apartment building.
- 42:26You wash your hands, you turn the tap off quickly, and
- 42:28suddenly the pipes behind the drywall violently bang like
- 42:31someone hit them with a literal hammer.
- 42:33Yeah. It's startling.
- 42:34I want to break down the exact mechanics of a water hammer.
- 42:37Let's go millisecond by millisecond.
- 42:39You've got water flowing steadily down a pipe and someone
- 42:42violently slams a lever valve shut.
- 42:45What exactly happens in that first instant right at the valve
- 42:48face? In that very first millisecond,
- 42:51the physical momentum of the fluid right at the valve face
- 42:53drops to absolute zero. It hits a solid wall.
- 42:58But water is heavy, it has mass, and it was moving with velocity.
- 43:03It possessed kinetic energy. That energy cannot just vanish.
- 43:06And because water is essentially incompressible, it cannot squish
- 43:10like a pillow to absorb the blow.
- 43:12So where does the energy go? It converts instantaneously into
- 43:16a massive localized pressure spike right against the metal
- 43:19face of the closed valve. Exactly.
- 43:22The formula for this spike is terrifyingly simple and linear.
- 43:26DP equals rho times V * A. So the change in pressure equals
- 43:31the density of the fluid multiplied by its original
- 43:34flowing velocity, multiplied by the speed of sound in that
- 43:37specific fluid. And here is where the chaos
- 43:39begins. That pressure spike doesn't just
- 43:41sit there. The water directly behind the
- 43:43valve has stopped, but the water 100 feet back up the pipe has no
- 43:48idea the valve is closed yet. Because the signal hasn't
- 43:50reached. It it is still rushing forward
- 43:52at full speed. Yeah.
- 43:53So that trailing water violently crashes into the stopped water.
- 43:56It compresses the water by a microscopic fraction of a
- 43:59millimeter and instantly stops itself.
- 44:01It's a microscopic high speed pileup.
- 44:03Car after car crashing into the back of the jamb.
- 44:06Exactly, and as each layer of water crashes and stops, the
- 44:10high pressure zone expands backward.
- 44:12This creates a high pressure shockwave that travels
- 44:14physically upstream, moving backward against the original
- 44:17flow direction at the speed of sound through water.
- 44:19OK, so I am picturing this invisible wall of intense
- 44:23destructive pressure traveling at thousands of miles an hour
- 44:27all the way back up the pipe. What happens when it finally
- 44:30reaches the open reservoir at the other end?
- 44:33When the shockwave reaches the open tank, the highly
- 44:36pressurized water inside the pipe suddenly finds itself open
- 44:39to the much lower normal resting pressure of the reservoir.
- 44:42Nature hates an imbalance. Oh man, so the highly
- 44:45pressurized pipe literally vomits the water backward into
- 44:48the tank? The flow actually reverses
- 44:50direction. It reverses, but this sudden
- 44:53release of water drops the pressure inside the pipe
- 44:55significantly below normal. Now you have a low pressure
- 44:58wave, and that low pressure wave travels all the way back
- 45:01downstream, back toward the closed valve.
- 45:03And when that low pressure wave finally hits the closed valve,
- 45:06the water is trying to pull away from the valve face.
- 45:09It creates a suction. A powerful suction.
- 45:11That suction then reaches back to the tank, pulls fresh water
- 45:14from the tank back into the pipe, and the entire high
- 45:17pressure slam cycle starts all over again.
- 45:20It bounces high pressure traveling upstream, low pressure
- 45:23traveling downstream, over and over again, bouncing back and
- 45:27forth. That loud, violent banging noise
- 45:30in your walls is the physical shockwave ricocheting off the
- 45:32ends of the plumbing, vibrating the metal until friction finally
- 45:36dampens the energy away. Now, you mentioned earlier that
- 45:38the speed of this destructive wave, the Sonic velocity,
- 45:42depends on the fluid. But we have to factor in the
- 45:45real world here. Water isn't flowing through a
- 45:48magical rigid vacuum, it's flowing inside a real, physical
- 45:52pipe. Right, a pipe made of steel or
- 45:54copper or concrete. And under that immense violent
- 45:57pressure spike of a water hammer, even solid steel is not
- 46:00perfectly rigid. The metal actually stretches.
- 46:02It does. When that shockwave hits, the
- 46:04massive internal pressure causes the pipe to expand radially.
- 46:07It literally bulges outward slightly as the wave passes.
- 46:11Right. And because the walls of the
- 46:13pipe yield and give way under the pressure, the volume of the
- 46:16pipe increases slightly. This makes the incompressible
- 46:19fluid appear to be compressible. It suddenly has extra room to
- 46:23give. So the environment itself
- 46:25changes the fundamental math of the fluid.
- 46:27We can't just use the standard stiffness or bulk modulus of
- 46:30water anymore. We have to calculate a modified
- 46:34bulk modulus K prime. Right AK prime.
- 46:37That actually includes the elasticity, the Young's modulus,
- 46:39of the solid pipe material itself.
- 46:41We have to treat the water in the pipe as a single combined
- 46:44flexible system. Correct.
- 46:46And because the pipe absorbs a significant portion of the shock
- 46:49energy by physically expanding, the velocity of the shockwave
- 46:52actually drops. If you have an ordinary steel or
- 46:55concrete water pipe, the wave velocity drops from the pure
- 46:58speed of sound and water down to roughly 600 to 1200 meters per
- 47:02second. Which is still incredibly fast,
- 47:04over 2000 miles an hour. So how do civil engineers
- 47:07prevent this phenomenon from just destroying municipal water
- 47:10grids and industrial plants? The simplest, most effective way
- 47:13is to just not slam the door. You use slow closing valves.
- 47:18Engineers use a complex formula called Olivais equation to
- 47:22calculate the exact maximum pressure spike based on the
- 47:25closure time of the valve. The goal is to stretch the
- 47:28physical closure time out to be longer than the time it takes
- 47:31the shockwave to make a full round trip up to the reservoir
- 47:34and back down the pipe. Exactly.
- 47:36So if you close it slowly, the low pressure reflection wave
- 47:40returning from the reservoir meets the high pressure wave
- 47:42still being generated by the closing valve and they
- 47:45destructively interfere. Yes, the waves essentially
- 47:48cancel each other out, bleeding off the kinetic energy safely
- 47:51and gradually, rather than in one catastrophic spike.
- 47:55I want to make sure we understand the real world stakes
- 47:58of this. This isn't just about avoiding
- 47:59an annoying banging noise in your apartment walls.
- 48:02Think about a massive industrial hydropower station built inside
- 48:06a mountain. OK, picturing it.
- 48:08The stakes are absolute life and death.
- 48:10A large hydropower station uses a giant, heavily reinforced pipe
- 48:15called a penstock. This pipe funnels thousands of
- 48:18tons of water down a massive vertical drop directly into a
- 48:22spinning turbine. Wow, OK.
- 48:24Now imagine a worst case scenario.
- 48:26The electrical grid suddenly drops its load.
- 48:29Maybe a substation blows out. The turbine suddenly has no
- 48:33resistance and it starts over speeding.
- 48:35If it spins too fast, it will tear itself apart.
- 48:38So the safety systems have to immediately choke off the water
- 48:41supply to save the turbine. But if those massive multi ton
- 48:45penstock valves at the bottom of the mountain are commanded to
- 48:47close too fast, the resulting water hammer wouldn't just make
- 48:51a noise. It would be a bomb.
- 48:52The kinetic energy of thousands of tons of falling water
- 48:55instantly converting to pressure would easily exceed the yield
- 48:58strength of the steel. The pressure spike would
- 49:00literally explode the Penn stock right out of the side of the
- 49:03mountain, flooding the facility and causing a catastrophic
- 49:06failure. That is genuinely terrifying.
- 49:09The invisible forces hidden in moving water are immense.
- 49:13OK, I want to take a step back and look at where we are.
- 49:15We have equations for perfect gases in perfect nozzles.
- 49:19We have equations for single shock waves traveling up
- 49:22straight pipes. Nimple shapes.
- 49:24But real world fluid dynamics is almost never a straight pipe.
- 49:27It's turbulent, chaotic air tumbling over a swept jet wing,
- 49:32or chaotic water violently spinning in three dimensions
- 49:35through a curved turbine blade. It's messy.
- 49:38Trying to calculate valve leakage, irregular pipe friction
- 49:42and three-dimensional interacting shock waves by hand.
- 49:45It's just impossible, right? It is fundamentally
- 49:48mathematically impossible able to solve by hand using the
- 49:51governing physics equations, which are known as the Navier
- 49:53Stokes equations. Why are they unsolvable?
- 49:55What makes them so uniquely difficult?
- 49:58Because the Navier Stokes equations are deeply non linear,
- 50:00in almost all practical intermediate flow regimes you
- 50:03have two competing forces. You have the inertia term, which
- 50:06is the forward momentum of the fluid, it's desire to keep
- 50:10moving in chaotic directions. OK, inertia and you have the
- 50:13viscosity term, which is the internal friction of the fluid,
- 50:16the stickiness that tries to smooth the flow out.
- 50:19Because inertia and viscosity fight each other in incredibly
- 50:22complex non linear ways, a change in one intimately changes
- 50:27the other in unpredictable loops.
- 50:29There are no analytical pencil and paper solutions for complex
- 50:32shapes. Which perfectly sets up the
- 50:34final act of our master class, teaching computers to see the
- 50:37flow. Because we could not solve the
- 50:40math by hand, humanity had to invent an entirely new field of
- 50:43engineering. We had to invent computational
- 50:45fluid dynamics, or CFD. Since we cannot calculate a
- 50:49single elegant equation that solves the entire volume of
- 50:52fluid holistically at once, we had to find a workaround.
- 50:55And what's the? Workaround The solution is
- 50:56something called the finite difference method.
- 50:59Instead of looking at the hole, we discretize the fluid.
- 51:01We break the continuous flowing mass of air or water down into a
- 51:05highly structured mathematical grid or a mesh made of millions
- 51:08of tiny individual boxes. I love comparing this to the
- 51:11evolution of photography. Back in the day, theoretical
- 51:15fluid dynamics was like an analog photograph.
- 51:18You were trying to capture the entire continuous image in one
- 51:22single holistic chemical exposure, one giant math
- 51:26equation. That is a highly accurate
- 51:28analogy. But CFD is like a modern digital
- 51:31camera sensor. Instead of looking at the whole
- 51:33picture, CFD breaks the universe down into millions of tiny
- 51:37independent pixels or cells. And the magic happens inside
- 51:41each individual pixel. To figure out what the fluid is
- 51:44doing inside one specific box, we use a mathematical trick
- 51:47called Taylor series expansion. Explain that simply, if
- 51:50possible. Essentially, Taylor series
- 51:52expansion allows us to take those impossible curvy non
- 51:55linear partial differential equations and approximate them
- 51:58into simpler straight line algebraic equations that a
- 52:01computer's processor can actually solve.
- 52:04We calculate the basic velocity and pressure in one single cell,
- 52:07and then we pass that result across the boundary to the
- 52:10neighboring cell. The computer solves every cell
- 52:13and then refreshes the entire frame thousands of times a
- 52:16second to simulate the passage of time.
- 52:18And there are different mathematical strategies for how
- 52:20you pass that information between the cells, right?
- 52:23The literature mentions forward difference, backward difference
- 52:26and central difference. Yes, let's look at central
- 52:29difference. If a computer wants to know the
- 52:31rate of change inside a target cell, Central difference looks
- 52:34at the cell immediately behind it, in the cell immediately in
- 52:37front of it, and averages them. It looks at both neighbors.
- 52:40That seems like a smart way to do it.
- 52:42It provides what we call 2nd order accuracy, meaning it is
- 52:46highly precise and captures fine details.
- 52:49However, it has a fatal flaw, which is if you're simulating
- 52:52flows with very high Reynolds numbers, meaning very fast,
- 52:56highly turbulent, chaotic flow, using central difference can
- 53:00actually cause the entire simulation to become numerically
- 53:03unstable. Unstable how?
- 53:05The tiny rounding errors between the cells amplify, bounce back
- 53:08and forth, and the math literally explodes to Infinity.
- 53:12The simulation crashes. So how do software engineers fix
- 53:15a crashing simulation? They switched to a technique
- 53:18called the upwind difference method.
- 53:20It relies on a very simple physical truth in a highly fast
- 53:24rushing flow. Almost all the relevant physical
- 53:27information is coming from upstream, blowing down toward
- 53:31you. What happened is downstream
- 53:32doesn't really affect you. Exactly so.
- 53:34The mathematical formula intentionally weights the
- 53:37upstream cells much heavier than the downstream cells.
- 53:41By favoring the upwind data, it suppresses those chaotic
- 53:44mathematical oscillations and keeps the simulation perfectly
- 53:47stable even at extreme speeds. OK, so that is how we build the
- 53:51grid and keep it stable. Now let's look at how we
- 53:54actually simulate physical events inside that digital
- 53:57space, starting with incompressible fluids like
- 54:00water. There is a historic landmark
- 54:03breakthrough in CFD called the MAS method or Marker and Sell
- 54:07method. Why was it so revolutionary?
- 54:09The Mac method was developed specifically to solve a problem
- 54:12that plagued early computers simulating free surfaces.
- 54:15Imagine trying to simulate water sloshing in 1/2 empty tank or a
- 54:19wave crashing on a beach. That sounds incredibly
- 54:21complicated. It is the boundary where the
- 54:24water meets the empty air is incredibly complex to track
- 54:27mathematically. In the Mac method, the computer
- 54:31generates the background mesh, but then it drops 10s of
- 54:34thousands of weightless marker particles into the simulated
- 54:37fluid. What do you mean by weightless?
- 54:39They don't have mass, they don't affect the physics equations,
- 54:42they're essentially digital dye. They're just dropped into the
- 54:45mathematical current, and they simply track exactly where the
- 54:49fluid physically exists within the grid as the velocity and
- 54:52pressure equations push them around over time.
- 54:54OK so if a cell has markers in it, it's water.
- 54:56If it's empty it's air, right? I remember seeing the visual
- 54:59they used to prove this method worked and it is iconic.
- 55:02They use used the Mac method to simulate a single liquid drop
- 55:06falling into a shallow puddle. We have all seen that famous
- 55:10high speed macro photograph where a milk drop hits a surface
- 55:14and the splash violently forms a perfect symmetrical little crown
- 55:17of liquid spikes shooting upward.
- 55:19Yes, the coronet splash. It is a highly complex chaotic 3
- 55:23dimensional physical event. The MSC method, using a
- 55:26primitive computer from decades ago, was able to perfectly
- 55:30simulate the physics of that crown forming spike by spike,
- 55:33droplet by droplet, purely by tracking those weightless marker
- 55:37particles through a grid. It proved that computers could
- 55:40replicate the chaotic beauty of fluid nature.
- 55:42It was a triumph. But as impressive as the MSC
- 55:46method is, what is truly fascinating is how much more
- 55:49difficult the math becomes when we switch back to our first
- 55:52topic, Compressible fluids. Supersonic air.
- 55:56You can't just track velocity and pressure with digital die
- 55:59anymore. This is the shockwaves.
- 56:00Exactly. Shock waves are violent
- 56:02discontinuities. To simulate them, the computer
- 56:05is forced to solve the equation of state and the energy equation
- 56:08simultaneously with mass and momentum.
- 56:11The temperature and density are wildly, violently fluctuating
- 56:14from cell to cell. So how do supercomputers compute
- 56:17a shockwave without crashing? We use what are called time
- 56:19marching methods. A prominent example is the
- 56:21McCormack method. It is incredibly useful because
- 56:25it uses conservation forms for mass, momentum, and energy.
- 56:29But to save massive amounts of computational power, it
- 56:32temporarily ignores the discussity of the air.
- 56:35It treats the air as frictionless.
- 56:37So it's a simplification. Yes, but this method is
- 56:40brilliant and highly efficient for predicting macro level
- 56:43aerodynamic events. Like calculating exactly where a
- 56:47detached bow shockwave will form a foot in front of the nose cone
- 56:50of a supersonic rocket. It paints the big picture.
- 56:53But what if the big picture isn't enough?
- 56:55What if you need to calculate the viscosity?
- 56:58What if the microscopic friction of the air rubbing against the
- 57:01metal is the whole point of the engineering problem?
- 57:03When you need absolute granular reality, you have to step up to
- 57:07the heavy artillery. You use something like the IAF
- 57:10method, the Implicit Approximate Factorization method, which is
- 57:13also commonly known as the Beam warming method.
- 57:16OK, the beam warming method. This is heavy duty mathematics.
- 57:18It does not take shortcuts. It solves the full complete
- 57:22viscous compressible Naver Stokes equations.
- 57:25Where would an engineer ever need that much punishing
- 57:27mathematical detail? Look inside a modern commercial
- 57:30jet engine. A turbofan.
- 57:33Imagine the massive titanium fan blades at the front of that
- 57:36engine. They are rotating at thousands
- 57:38of revolutions per minute. At those rotational speeds, the
- 57:41tips of the blades are actually breaking the sound barrier,
- 57:44while the base of the blades near the hub are moving
- 57:47subsonic. So the air flowing through the
- 57:49fan is this nightmare chaotic mixture of subsonic and
- 57:53supersonic flow at the exact same time.
- 57:56Transonic flow. Precisely, and in that transonic
- 57:59environment, microscopic shock waves are constantly forming,
- 58:03detaching and interacting with the turbulent friction boundary
- 58:07layer right on the curved surface of the metal plates.
- 58:10That is intensely complicated. It is the beam warming method
- 58:14running on parallel supercomputers can physically
- 58:16map that incredibly complex density distribution.
- 58:20It allows engineers to literally see the friction in the
- 58:22microscopic shock waves rubbing against the titanium, allowing
- 58:25them to tweak the curve of the blade by a fraction of a
- 58:28millimeter to save millions of dollars in jet fuel.
- 58:30It is almost God like the ability to take the invisible,
- 58:34chaotic, destructive forces of nature, forces that would RIP a
- 58:38building apart and map them out and perfect digital color on a
- 58:41monitor. It's an incredible time to be an
- 58:43engineer. OK, as we wrap up this master
- 58:45class deep dive, I want you to give us something to chew on.
- 58:48We have covered a massive amount of ground.
- 58:51What is the ultimate destination of all this physics and computer
- 58:55science? I think you just have to look
- 58:57back at the incredible journey we just took.
- 58:59We started with scientists in the 1800s playing with sealed
- 59:02boxes, realizing that pressure times volume equals temperature,
- 59:07just the basic rudimentary relationship.
- 59:10And we traced that intellectual lineage through supersonic
- 59:13rockets, through exploding water pipes, all the way to
- 59:16supercomputers calculating 3 dimensional friction on a
- 59:19rotating jet turbine. Here's the provocative thought
- 59:22for the future. As computing power marches
- 59:25relentlessly toward quantum levels, the mathematical pixels
- 59:28in our CFD meshes will become infinitely small.
- 59:31We might soon reach a point where we can run fluid dynamic
- 59:35models so perfectly precise, so indistinguishable from reality,
- 59:40that physical wind tunnels and metal prototypes become entirely
- 59:44historical relics. You mean never building a test
- 59:46model again? Imagine an aerospace engineer
- 59:49designing A next generation supersonic aircraft.
- 59:52They can test every single physical parameter, stress the
- 59:55metal against virtual shock waves, and fly the plane for
- 59:58thousands of simulated hours in every weather condition
- 1:00:01imaginable, entirely inside a silicon chip, before a single
- 1:00:04piece of actual metal is ever cut on a factory floor.
- 1:00:08We will transition from reacting to the physical world to
- 1:00:11generating it perfectly in code. That is profound, designing and
- 1:00:15perfecting reality in the digital realm before it ever
- 1:00:17physically exists. Well designed junkies.
- 1:00:19That brings us to the end of our master class on the invisible
- 1:00:21forces around us. If you haven't already, give us
- 1:00:23a like and follow. And please leave us a comment if
- 1:00:25you enjoy the podcast. We love to hear from the folks
- 1:00:27who make this possible.