In this chapter, we introduce the tools required to model the dynamics of fluids. This will allow us to model how objects can float, how water flows through a pipe, and how airplane wings create lift. We will start by introducing the concept of pressure and modelling static fluids (hydrostatics) before developing models for fluids that flow (hydrodynamics). Fluids are generally defined as the phase of matter in which atoms (or molecules) are only loosely bound to each other, such as in gases or liquids. Most of the formalism that we develop will apply to any fluid (gas, liquid, plasma), although we will often restrict ourselves to modelling the most simple situations (e.g. laminar flow of an incompressible liquid).
where F⊥ is the component of the force perpendicular to the surface. The SI unit for pressure is the Pascal (Pa). Pressure is related to the area, A, over which a force is exerted, and can be thought of as a measure of how concentrated that force is. For example, a force of 10N exerted through a needle (a small area) will result in a much larger pressure than if that force was exerted by a flat hand (a larger area).
When a force is exerted on a fluid, it creates pressure that we model as being everywhere in the fluid.
For each element in the fluid, the pressure from the surrounding fluid exerts an inwards force on the element from all directions (see Figure 2). In reaction, the element exerts an outwards force in all directions, and these forces act on neighbouring elements.
This is somewhat analogous to the tension that exists everywhere in a rope, where each element of the rope experiences forces from the neighbouring elements in rope that try to “pull it apart”. Pressure can be thought of as a “negative” tension, in that the material under pressure is experiencing forces trying to collapse the element onto itself, rather than trying to pull it apart. To create a tension in a rope, one would exert an outwards force on the rope (in order to stretch it), so that the rope exerts an inwards force in reaction. In order to create pressure in a fluid, one must exert an inwards force on the fluid, which then exerts an outwards force in reaction.
If we consider a small cubic volume of fluid, as depicted in the centre of Figure 2, that element of fluid will experience inwards forces in all directions from the pressure in the surrounding fluid, as illustrated by the arrows. If the forces from the pressure result in no net force on the fluid element, then we say that the fluid is in hydrostatic equilibrium, and the element of fluid will be at rest in an inertial frame of reference.
Consider, instead, an element of fluid that at the edge of a container for the fluid (e.g. a cup of water), as depicted in Figure 3.
In this case, there is no fluid on the right-hand side of the fluid element to exert a force towards the left. If the fluid element is in equilibrium, it must then be the container that exerts that force, Fcontainer, on the fluid. By Newton’s Third Law, the element of fluid exerts an outwards force on the container. This is true at all points on the surface of container, which will all experience an outwards force from the pressure of the fluid. If the pressure is constant over a surface, the magnitude of the outwards force on the surface will be equal to the pressure of the fluid multiplied by the area of that surface.
If you place an empty sealed tin can under water, the water will exert a pressure on all of the surfaces of the tin can that leads to a net inwards force on all surfaces of the tin can. If the water pressure is high enough, the tin can will get crushed. If, on the other hand, the tin can is allowed to fill with water, it will not get crushed, as the water inside the tin can will have the same pressure as the water outside the tin can and will exert an equal net outwards force on all surfaces of the tin can. The net force on each surface of the can will be zero, and the tin can will not get crushed, no matter how high the water pressure is.
In general, if there is an interface with fluid on either side of it a different pressures, it is the difference in pressure on either side of the interface that determines the net force exerted on the interface, rather than the absolute pressure.
When discussing Figure 2, we argued that the fluid exerts an equal force, from all directions, on the fluid element, so that the net force on the fluid element is zero. This is not quite correct in the presence of gravity, where the fluid element will have a weight. Thus, if the fluid element is to be in equilibrium, the upwards force (and pressure) from the fluid below must be higher than that from the fluid above the fluid element.
Figure 4 shows an element of fluid that has a height h and a surface area A in the horizontal plane. The pressure, P2, in the fluid below the fluid element must be higher than the pressure, P1, above the fluid element, if the fluid element is in equilibrium.
The element of fluid has a total mass, m, given by:
where, V=Ah, is the volume of the fluid, and, ρ, its density.
The net (horizontal) force exerted by the external fluid on the fluid element is zero along the vertical surfaces. Let P1 be the pressure in the fluid above the fluid element, and P2 be the pressure below the fluid element. If we choose a y axis that is positive upwards and the fluid element does not accelerate in the vertical direction, then the y component of Newton’s Second Law, written for the fluid element, is:
where we used the fact that the force resulting from a pressure is given by the pressure multiplied by the area over which it is exerted. We thus find that the difference in pressure due to gravity in a fluid between two positions, y2 and y1, is given by:
where the y axis is defined to increase in the upwards direction. Since the pressure in the fluid depends on the location in the fluid, we say that there is a “pressure gradient” in the fluid.
We have assumed that the density of the fluid, ρ, is constant, and that the fluid cannot be compressed. This is a very good approximation for a liquid such as water, but not for a gas, whose density will depend on its pressure. If the fluid were a gas (e.g. a column of air in our atmosphere), both the density and the pressure will change as a function of height. We can easily take this into account in our model, if we consider the fluid element to have a very small height, dy, instead of the finite height, h, as in the derivation above. A fluid element with an infinitesimal height, dy, is illustrated in Figure 6.
In the very small height, dy, the density of the fluid, ρ, can be taken to be constant, and the infinitesimal element of fluid will have a mass dm:
We can model the pressure exerted by the fluid above the fluid element as P+dP, and the pressure exerted by the fluid below as P, where dP is a small (negative) change in pressure[38]. The y component of Newton’s Second Law written for the infinitesimal fluid element is thus:
This tells us that the rate of change of pressure with increasing y is negative; in other words, the pressure decreases as the elevation increases, as we had already concluded. We can integrate the equation to obtain the change in pressure in going from y1 to y2:
If the density, ρ, is constant, then this leads to (4). Note that, thus far, we have only modelled how pressure in a fluid changes with height, but we have not determined the absolute pressure in a fluid.
Pascal’s Principle states that if an external pressure is exerted on a fluid, the pressure everywhere in the fluid increases by that amount. For example, if a fluid is contained in a piston with a cross-section area, A, and a force, F, is exerted on the piston (Figure 7), then the pressure everywhere in the fluid increase by F/A.
If we wish to determine the absolute pressure in the water at some depth, h, in the ocean, we need to include the fact that the Earth’s atmosphere exerts a net downwards force on the surface of the ocean in addition to the fact that the pressure changes with depth due to gravity. The pressure from the air in the Earth’s atmosphere is called “atmospheric pressure”, and depends on a variety of conditions, such as the weather. The average pressure from the atmosphere is P0=1.013×105Pa. If the atmospheric pressure is P0 at the surface of the ocean, then the pressure at some depth, h, is given by:
In this section, we describe how one can design instruments to measure pressure. The most straightforward device is a manometer, which is constructed using a U-shaped tube filled with a fluid of known density, ρ, as shown in Figure 10.
A manometer can be used to measure a pressure P relative to atmospheric pressure, P0. One end of the tube is open to atmospheric pressure, and the other is connected to the fluid (e.g. a gas) for which we want to measure the pressure. If the pressure being measured is larger than atmospheric pressure, the fluid in the manometer will experience a greater downwards force on the side of the pressure to be measured than on the side open to atmospheric pressure, as shown in Figure 10. There will be a difference, h, in the level of the fluid on each side of the tube, which is directly proportional to the difference in pressure between the two sides of the tube.
Consider the point in the fluid at location B in Figure 10, where the pressure is PB=P, the pressure to be measured. The point in the fluid at location A, which is at the same height in the fluid, must have the same pressure as point B. We can write the pressure at point A, PA, as the sum of the atmospheric pressure and the pressure from the column of water of height, h:
The difference between a pressure and the atmospheric pressure is called “gauge pressure”, and is all that we can measure if we do not know the absolute value of the atmospheric pressure. Using a manometer, the gauge pressure is given by ρgh, whereas the “absolute pressure”, P, is given by adding the atmospheric pressure to the gauge pressure, P=P0+ρgh. Most pressure measuring devices (“pressure gauges”), measure pressure relative to atmospheric pressure, using a similar mechanism.
The atmospheric pressure at a location on Earth varies based on the weather. A barometer is an instrument designed to measure the atmospheric pressure. A simple barometer can be built from a manometer, with one end closed, as illustrated in Figure 11.
One end of the manometer is sealed on a day where the atmospheric pressure is, say, P0, while the other end of the tube is left open. The height difference, h, between the fluid in either side of the tube is a measure of how different the current atmospheric pressure, P1, is relative to the pressure, P0, when the manometer was sealed. In Figure 11, the barometer is shown on a day where the atmospheric pressure is lower than on the day the manometer was sealed. The difference in pressure is given by:
if we define h to be positive when the side with the pressure P0 is higher (so h is negative in Figure 11 and P1 is less than P0).
We can also measure the absolute atmospheric pressure if we evacuate the air out of the sealed end of the tube, so that P0=0. When doing so, the difference in height between the fluid on either side of the manometer is a measure of the absolute atmospheric pressure.
In this section, we examine how the pressure gradient in a fluid leads to a force of buoyancy on an object that is immersed in the fluid.
In the left panel of Figure 14, we show a hemi-spherical element of fluid with a volume V. The weight of the element of fluid, Fg, is supported by the net upwards force, FB, exerted by the pressure of the fluid surrounding the fluid element. The mass, M, of the element of fluid is given by:
Now, suppose that the fluid element is “displaced” and replaced by the hull of a boat, as shown in the right panel of Figure 14. The net upwards force from the pressure of the fluid must remain the same, FB, but that force is now exerted on the hull of the boat. We call that force the force of “buoyancy”, which is the reason that a boat can float and the reason that you feel lighter when walking in a swimming pool than on land.
Thus, if an object displaces a volume, V, of a fluid with density ρ, when immersed in the fluid, that object will experience an upwards force of buoyancy, FB, with magnitude:
This “principle” was originally discovered by Archimedes, who stated that the force of buoyancy is equal to the weight of the displaced fluid. Note that we drew the fluid element at the surface of the fluid, but this is not required, and a force of buoyancy will be present if the object is completely immersed in the liquid. If you refer back to Figure 4, you will recall that the net upwards force on an element of fluid must be equal to its weight, even if the fluid element is completely immersed.
In the previous sections we developed “hydrostatic” models for fluids when those fluids are at rest (in some inertial reference frame). In this section, we develop “hydrodynamic” models to discuss what happens when fluids flow. We will restrict our models to fluids that flow in a “laminar” fashion, rather than a “turbulent” fashion.
Laminar flow is the flow of a fluid when each particle in the fluid follows a path that can be represented by a line (a “streamline”). Turbulent flow is the flow of a fluid where particles can follow rather complex paths, usually involving “Eddy currents” (little whirlpools). The two types of flow are illustrated in Figure 15.
Consider the laminar flow of a fluid through a pipe whose cross-sectional area narrows from A1 to A2 in the direction of flow, as illustrated in Figure 16.
The particles that make up the fluid have a speed v1 at the wide end of the pipe and speed v2 at the narrow end. The equation of continuity is based on the premise that the fluid that enters the pipe must exit the pipe, as there is nowhere else for the fluid to go. That is, if during a period of time, Δt, a mass, Δm, of fluid enters the wide end of the pipe, then during that same period of time, the same mass of fluid must exit the narrow end of the pipe.
During a period of time, Δt, the fluid at the wide end of the pipe will travel a distance l1=v1Δt. Thus, a volume of fluid, ΔV1, will enter the wide end of the pipe:
If the fluid is compressible, its density can change. Let ρ1 be the density of the fluid at the wide end of the pipe and ρ2 be the density of the fluid at the narrow end. The mass of fluid, Δm, entering the wide end of the pipe is given by:
where A is the cross-sectional area of the surface through which a fluid with speed, v, flows[39]. Q has the dimension of volume per time, and corresponds to the volume of fluid passing through the cross section A per unit time. For an incompressible fluid, the equation of continuity is thus equivalent to stating that the volumetric flow, Q, of the fluid is a constant.
In this section, we examine how the pressure and speed of a fluid change as the fluid flows. We will restrict ourselves to discussing the laminar flow of an incompressible fluid with no friction. Bernoulli was the first to quantitatively describe the flow of incompressible fluids, and we will show in this section how to derive “Bernoulli’s Principle”.
Consider the laminar flow of an incompressible fluid through a pipe that changes height, from y1 to y2, as well as cross-sectional area, from A1 to A2, as shown in Figure 18. The figure shows an element of fluid, in blue, as it moves through the pipe. The top panel corresponds to the location of the fluid element at time t=0, whereas the bottom panel shows the location of the element of fluid at time t=Δt.
To model how the fluid moves through this pipe, we can use energy and the Work-Energy Theorem. We start by considering the amount of work done on the element of fluid as it moves from the position in the top panel to the position in the bottom panel.
The fluid that is to the left of the element of fluid exerts a pressure, P1, on the fluid element that leads to a net force, F1, towards the right. Similarly, the fluid to the right of the element of fluid exerts a net force F2 in the opposite direction, due to the pressure P2 on that side of the fluid element.
In a period of time, Δt, the left part of the fluid element will move a distance l1=v1Δt, while the right part of the fluid element will move a distance l2=v2Δt. We can calculate the work done by each force, defining positive work to be in the direction of motion:
Gravity will also do (negative) work on the fluid as it changes height. In a period of time, Δt, a mass of fluid, Δm, will move from position y=y1 to position y=y2. The mass of fluid that changes height is given by the part of the fluid that moves a distance, l1, on the right side of the pipe:
Since the locations 1 and 2 that we chose are arbitrary, we can state that, for laminar incompressible flow, the following quantity evaluated at any position is a constant:
This statement is what we call Bernoulli’s Equation, and is equivalent to conservation of energy for the fluid. If the fluid is not flowing (v1=v2=0), then this is equivalent to the statement of hydrostatic equilibrium that we derived in (4):
If a fluid is flowing at constant height such that v2>v1 (as in Figure 16), then P2<P1; that is, the pressure in the fluid is lower if the fluid is flowing faster. Note that P is the pressure inside the fluid and is not related to the force that would be exerted by the fluid if it were to collide with an object. It makes sense that the fluid has a lower pressure where it is moving faster, because the net force exerted on the fluid is related to the difference in pressure on either side of the fluid. The fluid will accelerate in the direction where pressure decreases, thus it will be moving faster when it is in a region of low pressure.
Bernoulli’s principle can be used to describe many phenomena. For example, an airplane wing (technically, an “airfoil”) creates lift because the pressure of the air above the wing is lower than the pressure of the air below the wing. This is illustrated in Figure 19, which shows that the laminar flow of the air creates a low pressure area above the wing. As the stream lines of air encounter the wing, those that are above the wing get compressed together, which leads to a faster speed of the air above the wing (equation of continuity). The resulting difference in air pressure above and below the wing results in a net upwards force on the wing.
Bernoulli’s principle also describes why the roof can be lifted off of a house in high winds (Figure 20, left panel). It is not the force of the wind against the roof that blows the roof off of a house; it is the difference in air pressure in the house (normal) and the pressure above the roof (low, due to the flowing wind), that results in a net upwards force on the roof. Bernoulli’s principle is also used to construct atomizers which allow liquid in a bottle to be sprayed (Figure 20, right panel). For example, perfume bottles often have a bulb connected to a tube/spout. When you squeeze the bulb, it causes the air in the tube to flow quickly, creating a low pressure in the vertical segment of the spout. The liquid is forced up by the pressure in the bottle; once the liquid arrives in the fast flowing air, it is sprayed out along with the air.
The following examples illustrate how to apply Bernoulli’s principle.
So far, we have assumed that fluids flow with no friction. In reality, the particles moving in a fluid exert internal friction on each other called “viscosity”. This can be modelled as the friction between different layers of fluid in a laminar flow. For example, you may notice that the water that flows in a wide river flows much faster in the middle of the river than near the river banks, where the water is almost stationary, as shown in Figure 22.
One can model the banks of the river as exerting a frictional force on the layer of water that is in contact with the banks. That layer then exerts a frictional force on the next layer closer to the centre of the river, and so on.
One can define a viscosity coefficient, η, based on measuring the force required to pull a plate past another plate when there is a fluid between the plates. Consider two plates that have an area, A, that are a distance l apart, and contain the fluid of interest between them, as illustrated in Figure 23.
The viscosity of the fluid is defined based on the force that is required to pull the top plate while the bottom plate remains immobile. The layer of fluid directly below the moving plate will move with the plate at a speed, v, while the layer of fluid immediately in contact will the stationary plate will also be stationary. Moving one plate will thus lead to a gradient (a change) in the speed of the fluid as a function of the position between the two plates. The magnitude of the force, F, required to move one plate with speed, v, was empirically determined to be proportional to the area of the plates, A, and the speed, v, while being inversely proportional to the distance, l, between the two plates:
If the viscosity of the fluid is zero, then no force is required to pull the plate. The more viscous the fluid, the more difficult it is to pull the top plate. You can experiment with this by comparing the force required to move a small piece of paper across the top of a puddle of water and across the top of honey.
The presence of viscosity means that any fluid that flows will lose mechanical energy due to internal friction (which will heat up the fluid). Thus, Bernoulli’s equation is not correct if the fluid has viscosity, as a fluid cannot flow through a horizontal pipe without a change in pressure to overcome the losses due to friction.
For the flow of an incompressible viscous fluid through a pipe, one can postulate that the flow rate, Q, is proportional to the change in pressure, ΔP, across the pipe:
where ΔP is taken as the positive difference between the pressure at either end of the pipe. The fluid flows from high pressure to low pressure. We can introduce a constant of proportionality, R, to be the “resistance of the pipe”, so that we can write:
where we wrote the constant of proportionality as 1/R, so that a larger value of R corresponds to a pipe with a higher resistance to flow. That is, for a given pressure difference, as one increases the resistance of the pipe, one decreases the flow rate through that pipe. The relationship above can be used to empirically determine the resistance of a pipe.
The flow through a pipe with a given resistance will be zero if there is no pressure gradient in the fluid along the pipe. Conversely, if there is no flow of fluid in the pipe, the pressure is the same everywhere in the pipe. We can thus also view a drop in pressure in a pipe to be the result of flow of liquid through the pipe. The pressure cannot drop across a horizontal pipe if there is no flow.
When you close the tap on your kitchen faucet, the pressure inside the faucet is close to the pressure in the main water line that supplies your house. As soon as you open the tap and allow water to flow, the pressure in your faucet drops to atmospheric pressure, and the resulting pressure gradient from the main supply forces water to flow out of the faucet. If you try to plug your kitchen faucet with your thumb and stop the flow of water, you will need to exert a force large enough to overcome the pressure that exists in the main water supply. You will find that it is practically impossible to stop the flow of water with your thumb, as the pressure in the main supply needs to be high enough to overcome the resistance of the pipes and still result in a usable flow of water.
Poiseuille first developed a model for the laminar flow of a liquid through a uniform horizontal cylindrical pipe of length, L, with a circular cross-section with radius r. He found that the resistance of such a pipe to a fluid of viscosity, η, is given by:
This makes some intuitive sense, as we expect more resistance (more impedance to flow), if the pipe is longer and if the fluid is more viscous (the resistance is zero if there is no viscosity). We further expect less resistance if the pipe has a larger radius. The resistance found by Poiseuille goes down as the fourth power of the radius. Thus, a pipe that is twice as wide will have a volume flow that is 24=16 times larger because of the reduced resistance.
The laminar flow rate, Q, of a viscous fluid through a pipe of length L and radius R, when there is a pressure difference ΔP, is given by:
where F⊥ is the component of the force perpendicular to the surface.
If a force is exerted on the particles in a fluid (e.g. gravity), a pressure will exist everywhere in the fluid. If the fluid is placed in a container, that pressure leads to an external force on all surfaces of the container.
If two fluids at different pressures exist on either side of an interface/object, the net force on that interface/object from the pressures of the fluids will be proportional to the difference in pressure of the fluids on either side.
A fluid is in hydrostatic equilibrium if the sum of the forces on any fluid element is zero. In the presence of gravity, this always leads to a vertical pressure gradient
Pascal’s Principle states that if an external pressure, P, is applied to one location in a fluid, then the pressure everywhere in the fluid increases by P.
If an object is immersed in a fluid, it will experience a force of buoyancy that is in the opposite direction to the gravitational field in that fluid. The magnitude of the buoyancy force is given by Archimedes’ Principle:
where, ρ, is the density of the fluid and, V, is the volume of the fluid displaced by the object (i.e. the volume of the part of the object that is immersed in the fluid).
We can distinguish between laminar and turbulent flow of fluids. In laminar flow, individual particles in the fluid follow well-defined streamlines. In turbulent flow, individual particle follow complicated paths that usually involve Eddy currents. In general, it is much easier to model the laminar flow of fluids.
The equation of continuity states that the mass flow rate of a fluid through a closed system must be the same everywhere in the system (no fluid can appear or disappear). For laminar flow of a fluid with density, ρ, flowing at speed, v, through a pipe with cross section, A, the mass flow rate is a constant:
A fluid is said to be incompressible if it has constant density. For a fluid of constant density, the volume flow rate, Q, must be constant everywhere in a closed system:
for the laminar flow of a fluid with no viscosity. P is the internal pressure of the fluid, v its speed, and y the height of the fluid relative to a fixed coordinate system. In particular, Bernoulli’s Principle implies that, for a constant height, the internal pressure of a fluid must decrease if its speed increases.
Viscosity, η, for the laminar flow of a fluid can be modelled as the result of the internal friction force between layers of the fluid. Because of viscosity, a fluid cannot flow in a horizontal pipe unless there is a difference in pressure across the pipe. Similarly, there will be no horizontal pressure gradient through a fluid unless the fluid is flowing. In general, the volume flow rate, Q, of an incompressible fluid through a pipe with resistance, R, is given by:
For the laminar flow of a fluid with viscosity, η, through a horizontal cylindrical pipe of length, L, and radius, r, the flow rate is given by Poiseuille’s equation:
a. The forces acting on the door are the force of buoyancy, the door’s weight, and the weights of Rebecca and Ryan, as shown in Figure 27.
We can combine the weight of the door and the weight of the people into the total weight, Fg. We choose the y axis to be positive upwards. The sum of the forces on the door in the y direction is given by:
For the door to float, the net force on the door must be greater than or equal to zero. We want to find the minimum buoyant force for them to float, so we set the net force equal to zero:
where the weight includes the mass of Rebecca (mR), Ryan (mr) and the door (md). We added the subscript W to the right side of the equation to remind ourselves that the buoyant force depends on the density and volume of the displaced water. We want to find the maximum density of the wood in order for Rebcca and Ryan to stay above the water’s surface. This means that the maximum volume of water that can be displaced is the volume of the door, Vw=Vd (so that the surface of the door is level with the surface of the water, as in Figure 25). We can rewrite the mass of the door in terms of its volume and density, and apply our condition that Vw=Vd :
The maximum density of the wood that would allow them to both float is 527kg/m3. Balsa wood has a density that is about 150kg/m3, so would allow them to survive. However, it is unlikely that a random floating door is made of balsa wood (although one would choose lighter materials when constructing a ship).
b. No, they could not both stay on the door because the density of oak is greater than the maximum density of 527kg/m3. We can find the amount of mass that can be added to the door (mA) in order for the person on it to stay above water:
The door can support an additional mass of 66kg, so either Rebecca or Ryan can survive if the other does not get on the door.
Solution 14.2
a. Given that the pressure in the patient’s veins is 80mmHg above atmospheric pressure, we want to find the pressure required at the other end of the needle so that we get the desired flow rate through the needle. We model the needle as a horizontal cylindrical pipe and assume that the saline solution exhibits laminar flow. We can therefore use Poiseuille’s equation:
The pressure at the exit of the needle, P2, is just the blood pressure (80mmHg+1atm). The radius of the needle is 0.60mm/2=0.30mm. The flow rate has to be in units of m3/s. The flow rate in the appropriate units is thus:
Note that, in the first line, we converted 80mmHg into Pascals.
b. We can easily determine the height of the IV bag that is required to give the desired pressure. We choose a coordinate system with a y axis that is vertical (positive upwards) with the origin at the location of the needle (Figure 28).
At the top of the solution in the IV bag, y0, the solution has a speed of zero and is at atmospheric pressure, P0=1atm. The velocity at the needle is 0, and the pressure is 13457Pa+1atm. Bernoulli’s principle states: