Applications of the Conservation Laws

J. Gordon Leishman

doi:https://doi.org/10.15394/eaglepub.2022.1066.n19

23 Applications of the Conservation Laws

Introduction

The practical use of the conservation equations of fluid dynamics is best learned by studying and understanding exemplar problems. Remember that exemplars of the field are “Key examples chosen so as to be typical of designated levels of quality of competence.”^[1] Such exemplars in fluid dynamics include applications to flow through pipes and ducts, venturimeters used for flow measurement, Pitot tubes and Pitot-static systems for flow speed and airspeed measurement, understanding the essential performance and efficiency of propulsion systems such as jet engines, and the forces on airfoil sections and other bodies. After working through a few exemplars, the systematic solution process for applying the conservation equations will become apparent so new and more ambitious problems can be tackled more confidently.

Generally speaking, all three conservation laws must be used to solve practical problems in fluid dynamics, i.e., using the governing equations to conserve mass, momentum, and energy. However, using all three equations will not always be required. Conservation of mass is generally always required in problem-solving. Energy principles will inevitably be needed whenever pressure is involved as an unknown in the problem. The momentum equation will be necessary whenever forces are involved. An auxiliary equation, the equation of state, may also be necessary for problems involving compressibility effects. The Bernoulli equation can be used instead of the formal energy equation, but only for steady, incompressible, inviscid flows without work or energy addition, and this restriction on its use should always be remembered.

Learning Objectives

Calculate the flow properties through a venturi and the contraction and test section of a wind tunnel using the fundamental conservation laws of aerodynamics.
Understand the principles associated with the operation of a Pitot tube and how to measure static and dynamic pressure,
Use the conservation laws in the integral form to find the drag on a body.
Analyze the essential performance of a propulsive system such as a jet engine.

Flow in a Venturi

Consider the steady, incompressible flow (i.e., $\varrho$ is a constant) through a convergent-divergent nozzle, as shown below, which is called a venturi. If the flow is assumed to be incompressible, then the fluid could be a liquid, or it would have to be a gas such as air but flowing at low velocities, i.e., all flow velocities need to be below a Mach number of 0.3, which may not be known a priori.

The advantages of a pressure drop in the throat of a venturi are used in many practical engineering applications.

The intake to a venturi is called the mouth, and the narrowest part is called the throat. Following Bernoulli’s principle, there is a pressure drop at the throat of a venturi where the flow speeds up, which can be used to advantage in several practical applications such as flow meters, carburetors, and wind tunnels.

The first step in analyzing the flow through a venturi is to draw a control surface around the venturi with $\vec{n} \, dS = d\vec{S}$ pointing out of the control volume by convention, as shown below. Let the inlet conditions be at section 1, the throat at section 2, and the outlet at section 3. The cross-sectional areas in each case are assumed to be known, such as by measurement. The flow can be considered steady and radially axisymmetric, i.e., an appropriate one-dimensional flow assumption.

Flow model for the flow through a venturi.

The steady form of the continuity equation in integral form is

(1) $\begin{equation*} \oiint_S \varrho \, \vec{V} \bigcdot d\vec{S} = 0 \end{equation*}$

In words: “What mass flow comes into the control volume per unit of time then leaves the control volume in the same time, i.e., no fluid mass accumulates inside the control volume.”

There is no flow over the walls of the venturi, so what fluid mass flow comes into the mouth of the venturi per unit of time also leaves simultaneously at the venturi’s exit. If the cross-sectional area variation is relatively moderate, the problem can be significantly simplified with a one-dimensional flow assumption. However, further consideration is also needed as to why a flow through a rapidly converging or expanding venturi (or duct) may not be as readily justifiable as a one-dimensional problem.

The flow enters the mouth of the venturi of area $A_1$ with an average velocity $V_1$ , has a velocity $V_2$ at the throat of area $A_2$ , and velocity $V_3$ at the exit of area $A_3$ . The fluid mass flow coming in is

(2) $\begin{equation*} \iint_{\rm intake} \varrho \, \vec{V} \bigcdot d\vec{S} = -\varrho_1 A_1 V_1 \end{equation*}$

The minus sign appears because the flow is opposite to the direction of $d\vec{S}$ at the entrance face. At the exit face, then

(3) $\begin{equation*} \iint_{\rm exit} \varrho \, \vec{V} \bigcdot d\vec{S} = \varrho_3 A_3 V_3 \end{equation*}$

Therefore, by using the continuity equation, then

(4) $\begin{equation*} \oiint_S \varrho \, \vec{V} \bigcdot d\vec{S} = \iint_{\rm intake} \varrho \, \vec{V} \bigcdot d\vec{S} + \iint_{\rm exit} \varrho \, \vec{V} \bigcdot d\vec{S} = 0 \end{equation*}$

which, for a one-dimensional flow, becomes

(5) $\begin{equation*} -\varrho_1 V_1 A_1 + \varrho_3 V_3 A_3 = 0 \end{equation*}$

If the flow is incompressible, $\varrho_1 = \varrho_2 = \varrho_3 = \varrho$ . Under these circumstances, the continuity equation can now be written as

(6) $\begin{equation*} -\varrho V_1 A_1 + \varrho V_3 A_3 = 0 \end{equation*}$

or

(7) $\begin{equation*} \varrho V_1 A_1 = \varrho V_3 A_3 = \overbigdot{m} \end{equation*}$

Also, it will be apparent that

(8) $\begin{equation*} \varrho V_1 A_1 = \varrho V_2 A_2 = \overbigdot{m} = \varrho \, Q \end{equation*}$

Canceling through the density gives

(9) $\begin{equation*} V_1 A_1 = V_2 A_2 = V_3 A_3 = Q \end{equation*}$

This equation states that the volume flow rate $Q$ through the venturi is constant.

It can also be seen from this latter equation that if the throat area decreases, then the flow velocity must increase there so that the continuity of the flow is satisfied, i.e., the mass and volume flow rate are constant at any one cross-section. Conversely, if the area of the throat increases, the flow velocity must decrease. Moreover, from the previous discussion of Bernoulli’s equation, remember that when the velocity increases, the pressure decreases; conversely, when the velocity decreases, the pressure increases, i.e.,

(10) $\begin{equation*} p + \frac{1}{2} \varrho V^2 = \mbox{\small constant} \end{equation*}$

or

(11) $\begin{equation*} p_1 + \frac{1}{2} \varrho V_1^2 = p_2 + \frac{1}{2} \varrho V_2^2 \end{equation*}$

(12) $\begin{equation*} p_1 - p_2 = \frac{1}{2} \varrho \left( V_2^2 - V_1^2 \right) \end{equation*}$

The pressure difference $p_1 - p_2$ can be measured across the mouth and throat of the venturi, which can be accomplished by drilling two small holes^[2] (i.e., pressure taps) perpendicular to the venturi walls at locations 1 and 2. The pressures are then measured by connecting tubes from these holes to a differential pressure gauge or the two sides (legs) of a U-tube manometer; see here for a simple experimental demonstration of this process. This pressure difference can then be related to the unknown velocity (or flow rate) by means of calculation or direct calibration using measurements.

From the continuity equation, then

(13) $\begin{equation*} V_{2} = V_{1}\left(\frac{A_1}{A_2}\right) \end{equation*}$

and from Bernoulli’s equation

(14) $\begin{equation*} p_{1} + \frac{1}{2}\varrho {V_{1}}^{2} = p_{2} + \frac{1}{2}\varrho {V_{2}}^{2} \end{equation*}$

Substituting gives

(15) $\begin{eqnarray*} p_{1} + \frac{1}{2} \varrho V_{1}^{2} & = & p_{2} + \frac{1}{2}\varrho V_{1}^{2} \left( \frac{A_1}{A_2} \right)^2 \\ p_{1} - p_{2} & = & \frac{1}{2} \varrho V_{1}^{2} \bigg( \left( \frac{A_1}{A_{2}} \right)^{2} - 1 \bigg) \end{eqnarray*}$

and so

(16) $\begin{equation*} Q = A_1 \sqrt{ \frac{2( p_{1} - p_{2})}{\varrho \bigg( \left( \displaystyle{\frac{A_{1}}{A_{2}} } \right)^{2} - 1 \bigg) } } \end{equation*}$

Therefore, knowing the pressure difference $p_{1} - p_{2}$ and the geometry of the venturi (the area ratio), then the flow rate $Q$ can be determined. Finally, notice that the change in pressure can be related to a change in the hydrostatic head of the fluid $\Delta h$ , using

(17) $\begin{equation*} p_{1} - p_{2} = \varrho_l \, g \, \Delta h \end{equation*}$

where $\varrho_l$ is the density of the liquid in the U-tube manometer.

Applications of a Venturi

A venturi finds many applications in practical engineering applications. Its primary characteristic is that the flow velocity is higher in the venturi’s throat. As a result, the pressure at the throat is lower than the atmospheric pressure, and this pressure difference is helpful for various purposes.

Flowmeter (Venturimeter)

As the previous analysis suggests, a venturi forms the basis for a device for measuring volume (or mass) flow rates. Such flowmeters, or venturimeters as they are often called, are commercially available in various sizes and capacities. They are widely used in the water, chemical, pharmaceutical, and oil and gas industries to measure fluid flow rates along ducts and pipes, an example being shown in the figure below. While venturimeters operate fundamentally on Bernoulli’s principle that the pressure decreases in the throat of the venturi, they are usually calibrated at the factory to account for any losses so that volumetric flow rates can be more accurately related to the measured pressure drop across the throat. There are usually four pressure taps at each section, which are pneumatically averaged together to give a more consistent and accurate reading of the static pressure drop between the entrance and the throat.

A classic venturimeter measures the mass or volumetric flow rate of a fluid through a pipe.

The volumetric flow rate can be obtained from the ideal fluid assumption derived previously with a correction factor to account for any viscous losses, i.e.,

(18) $\begin{equation*} Q = K_{\rm cal} \, A_1 \sqrt{ \frac{2( p_{1} - p_{2})}{\varrho \bigg( \left( \displaystyle{\frac{A_{1}}{A_{2}} } \right)^{2} - 1 \bigg) } } \end{equation*}$

where $K_{\rm cal}$ is a calibration factor. The value of $K_{\rm cal}$ is usually close to 1, but it will depend on the specific venturi.

Suction Source for Flight Instruments

In early aircraft, a venturi (or sometimes a pair of venturi tubes) was mounted on the side of an aircraft’s fuselage, as shown in the photograph below. The venturi provided low (suction) pressure for air-driven gyroscopic instruments such as the artificial horizon and directional gyro. Modern light aircraft, however, are fitted with mechanical engine-driven vacuum pumps to provide the needed suction pressure, so venturis are no longer required.

A pair of venturi tubes mounted on the side of an early aircraft provided low (suction) pressure to drive gyroscopic flight instruments.

From the continuity equation, then the velocity at the throat, $V_t$ , for a given airspeed, $V_{\infty}$ , is

(19) $\begin{equation*} V_{t} = V_{\infty}\left(\frac{A_i}{A_t}\right) \end{equation*}$

where $A_i$ is the inlet area, and $A_t$ is the area of the throat. From the Bernoulli equation, then

(20) $\begin{equation*} p_{\infty} + \frac{1}{2}\varrho {V_{\infty}}^{2} = p_{t} + \frac{1}{2}\varrho {V_{t}}^{2} = p_{t} + \frac{1}{2}\varrho V_{\infty}^2\left(\frac{A_1}{A_2}\right)^2 \end{equation*}$

Therefore, the suction pressure available at the throat relative to static pressure is

(21) $\begin{equation*} p_{t} - p_{\infty} = \frac{1}{2}\varrho {V_{\infty}}^{2} \left( 1 - \left(\frac{A_1}{A_2}\right)^2 \right) \end{equation*}$

which shows that it decreases with the square of the airspeed at a given density altitude.

Carburetor

A carburetor provides an air-fuel mixture to a piston (reciprocating) engine. A carburetor has a venturi through which incoming air is mixed with fuel. Fuel from a tank is derived to a fuel reservoir where the fuel level is maintained using a float. The fuel delivery line opens into the venturi in the carburetor’s throat, which first atomizes the fuel. The lower pressure in this region helps to draw fuel into the airstream, mix it with the air, and vaporize it downstream of the throat before the resulting fuel-air mixture goes into the cylinders to be burned. Metering the fuel needs fine adjustment to ensure that the correct air-fuel mixture arrives at each of the cylinders, which will maximize power from the engine and also minimize hydrocarbon emissions.

Wind-Tunnel

A low-speed wind tunnel is a large venturi where the airflow is driven by a fan connected to a motor drive. The wind tunnel fan is a large propeller, and it is specifically designed to move the air through the test section efficiently, as shown in the figure below. This is an open-return or Eiffel type of wind tunnel, the other type being a closed-return or loop type.

The airflow enters the mouth of area $A_1$ at a flow velocity $V_1$ with pressure $p_1$ . The wind tunnel then contracts to a smaller area $A_2$ at the test section where the velocity has increased to $V_2$ . The velocity in the test section must increase if continuity is satisfied. The model (such as a wing or complete airplane model) is placed in the test section, where its aerodynamic characteristics are measured. The flow then passes downstream into a diverging duct called a diffuser, where just before the fan, the area is $A_3$ , the velocity reduces to $V_3$ , and the pressure increases to $p_3$ .

From the continuity equation, the air velocity in the test section is

(22) $\begin{equation*} V_2 = \left( \frac{A_1}{A_2} \right) V_1 \end{equation*}$

In turn, the velocity at the exit of the diffuser before the fan is

(23) $\begin{equation*} V_3 = \left( \frac{A_2}{A_3} \right) V_2 \end{equation*}$

The pressure at various locations in the wind tunnel is related to the velocity by Bernoulli’s equation

(24) $\begin{equation*} p_1 + \frac{1}{2} \varrho V_1^2 = p_2 + \frac{1}{2} \varrho V_2^2 = p_3 + \frac{1}{2} \varrho V_3^2 \end{equation*}$

The velocity in the working or test section (the most important quantity) can be related to the pressure drop across sections 1 and 2, i.e., the pressure drop between the mouth of the contraction section and the test section. From the Bernoulli equation, then

(25) $\begin{equation*} V_2^2 = \frac{2}{\varrho} \left( p_2 - p_1 \right) + V_1^2 \end{equation*}$

This means that

(26) $\begin{equation*} V_2^2 = \frac{2}{\varrho} \left( p_1 - p_2 \right) + \left( \frac{A_2}{A_1} \right)^2 V_2^2 \end{equation*}$

Solving for $V_2$ (the flow velocity in the test section) gives

(27) $\begin{equation*} V_{2} = \sqrt{ \frac{2\left(p_{1} - p_{2}\right)}{\varrho \bigg( 1 - \left( \displaystyle {\frac{A_{2}}{A_{1}} }\right)^{2} \bigg) } } \end{equation*}$

The area ratio $A_2/A_1$ is fixed for a given wind tunnel. Recall that the density is a constant for an incompressible flow. Its value can be obtained from measurements of static pressure and temperature in conjunction with the equation of state.

Therefore, the above equation can be used to determine the velocity in the test section by measuring the pressure drop from the intake to the contraction and the test section. This latter technique is used in most low-speed wind tunnels to measure the flow speed in the test section, i.e., the drop in pressure between some point on the contraction and the test section is measured. This pressure drop is then related to the flow velocity in the test section using Bernoulli’s equation and confirmed by calibration. In the calibration, a Pitot probe (see discussion below) is placed in the test section, and the pressure drop is measured for a range of flow speeds; any discrepancy leads to a calibration factor, $K_{\rm cal}$ , that can be used to determine the flow speed more accurately, i.e., using

(28) $\begin{equation*} V_{2} = \sqrt{ \frac{2 \, K_{\rm cal} \left(p_{1} - p_{2}\right)}{\varrho \bigg( 1 - \left( \displaystyle {\frac{A_{2}}{A_{1}} }\right)^{2} \bigg) } } \end{equation*}$

In practice, the calibration factor for a low-speed wind tunnel is close to 1.

Cavitating Venturi

Cavitating venturis can be used to control or limit the flow rate of a liquid through a pipe. The throat of a cavitating venturi is sized such that for a given differential pressure drop between the inlet and the throat, the pressure at the throat is reduced to its vapor pressure point, $p_v$ , as illustrated in the figure below. The resulting vapor and bubbles that form, called cavitation, which is like boiling without adding heat, will then limit the mass flow rate through the venturi to prevent any additional increase in flow rate for a given upstream inlet pressure. As the bubbles travel downstream into the throat, they eventually burst, and the static pressure recovers to the pre-cavitation values. The mass flow rate through the pipe is not dependent on the downstream pressure, so the mass flow can then be controlled by changing the inlet pressure and, hence, the amount of cavitation in the throat. Cavitating venturimeters are often used as propellant mass flow limiters in rocket engines.

The cavitation number is usually expressed as

(29) $\begin{equation*} C_{\rm cav} = \dfrac{p - p_v}{\frac{1}{2} \varrho V_t^2} \end{equation*}$

where $V_t$ is the flow velocity at the throat. This is a form of pressure coefficient, so cavitation will first occur at $C_{\rm cav} = 1$ . The volumetric flow rate through a venturi is given by

(30) $\begin{equation*} Q = A_1 \sqrt{ \frac{2( p_{e} - p_{t})}{\varrho \bigg( \left( \displaystyle{\frac{A_{e}}{A_{t}} } \right)^{2} - 1 \bigg) } } \end{equation*}$

where $p_e$ is the entrance pressure to the venturi, and $p_t$ is the pressure in the throat. The entrance area is $A_t$ , and the throat area is $A_t$ , both of which can be measured. If the pressure in the throat equals the vapor pressure, i.e., $p_t = p_v$ and $C_{\rm cav} = 1$ , then the flow rate will become cavitation-limited to a maximum value of

(31) $\begin{equation*} Q_{\rm max} = A_e \sqrt{ \frac{2( p_{e} - p_{v})}{\varrho \bigg( \left( \displaystyle{\frac{A_{e}}{A_{t}} } \right)^{2} - 1 \bigg) } } \end{equation*}$

After reaching this point, further attempts to increase the upstream pressure will be ineffective in increasing the flow rate. Vapor pressures, $p_v$ , for different liquids at different temperatures can be found in various online resources and in the table below. Knowing the vapor pressure also allows the venturi to be designed to limit volumetric or mass flow rate. For example, for given values of $Q_{\rm max}$ , $p_e$ , and $A_1$ , the throat area $A_t$ can be determined to limit the flow rate from the onset of cavitation.

**Vapor pressures of different liquids at a temperature of 25°C (77°F)**
Liquid	Vapor pressure (Pa)	Vapor pressure (lb/ft²)
Water	3,166.6	0.0666
JET-A	~ 2,000	~ 0.04206
100LL gasoline	~ 40,000	~ 0.8413
Ethanol	5,853.7	0.1227
Acetone	24,266.4	0.5085
Diethyl ether	53,465.3	1.1195
Methanol	16,931.5	0.356316
Hexane	19,995	0.4201
Chloroform	26,384.9	0.5543
Carbon tetrachloride	12,125.7	0.2555
Benzine	12,661.2	0.2673

Siphon (Syphon)

A siphon, also spelled syphon, is a device used to transfer liquids from one level to another. It consists of a tube or hose bent in a U-shape, with one end placed in some liquid, such as water, and the other set at a lower level, as shown in the figure below. This arrangement allows the water to flow from the tank to a lower level, which relies on the principle of gravity and hydrostatic pressure. To start the flow, the siphon tube must be filled with water, which can be done by suctioning the air out of the tube or by filling it with water and submerging one end in the tank. Once the flow starts, it will continue until the tank is emptied.

The principle of a siphon is used to transfer a liquid from a container placed at a higher level to a lower level.

It is often claimed that a siphon works on the principle of air pressure, but this is incorrect. To see why, the Bernoulli equation applied in the fluid between the tank and the outlet from the siphon gives

(32) $\begin{equation*} p_1 + \frac{1}{2} \varrho V_1^2 + \varrho \, g \, z_1 = p_2 + \frac{1}{2} \varrho V_2^2 + \varrho \, g \, z_2 \end{equation*}$

The change in air pressure between the levels at $z_1$ and $z_2$ is minimal, so it can be assumed that $p_1 = p_2$ . Also, if the fluid in the tank is assumed to be of much greater volume than the volume flow rate out of the siphon, then $z_2$ = constant. Therefore, under these quasi-steady assumptions, the Bernoulli equation becomes

(33) $\begin{equation*} p_1 + \frac{1}{2} \varrho V_1^2 + \varrho \, g \, z_1 = p_1 + 0 + \varrho \, g \, z_2 \end{equation*}$

or

(34) $\begin{equation*} \frac{1}{2} V_1^2 + g \, z_1 = g \, z_2 \end{equation*}$

Rearranging gives the outlet flow velocity as

(35) $\begin{equation*} V_1 = \sqrt{ 2 \, g \left( z_2 - z_1 \right)} = \sqrt{ 2 \, g \, h} \end{equation*}$

which is called Torricelli’s law.^[3] Notice that the flow velocity depends on the difference in the height of the water, $h$ , and acceleration under gravity, so the higher the height difference, the faster the water will flow. There is no significant change in atmospheric pressure between the two levels. Therefore, the siphon effect is produced by the hydrostatic pressure difference within the liquid, not the (slight) difference in atmospheric air pressure.

However, considering viscous effects (and hence losses) will change this result because there will be some small pressure drop along the length of the siphon tube. Calculating such internal pressure losses and flow velocity requires a viscous flow theory, which will be considered in a later chapter of this ebook.

Total, Static, & Dynamic Pressures

It must now be explained more precisely what a fluid’s total and static pressures mean. Consider a flow moving with velocity $V$ at a pressure $p_{s}$ . If moving with the velocity of the airflow, then the pressure that will be felt is $p_{s}$ . This latter pressure is called static pressure and measures the effects produced by the random motion of the fluid molecules. Now, if the boundary or “wall” of a flow is considered (e.g., the wall of a duct, pipe, or wind tunnel), and a small hole in the wall is drilled in perpendicular to the flow, then the pressure that is measured there would be the static pressure $p_{s}$ .

Suppose a U-tube manometer is connected, as shown in the figure below, and the open end outside the flow (i.e., the reference end) of the manometer is at lower static pressure than in the flow. In that case, a difference in height $h_s$ will be shown on the manometer. The pressure is the static pressure between the flow and the ambient reference pressure, called the gauge pressure. If the static pressures inside and outside are equal, then $h_s = 0$ .

Concepts of the measurement of static pressure, total pressure, and dynamic pressure.

Now consider a Pitot tube (after the hydraulic engineer, Henri Pitot) inserted into the flow, which is an open tube with its end facing into the flow. The fluid cannot flow through the tube so it will stagnate in the tube, i.e., the flow velocity will go to zero at the entrance to the tube because there is no place for the fluid to flow. Hence, the streamline that infringes directly on the mouth of the Pitot tube will have zero velocity; this location is called a stagnation point. In this case, the pressure that is measured is called the total pressure, $p_T$ , which is the sum of the dynamic pressure and the static pressure, i.e.,

(36) $\begin{equation*} p_T = p_s + \frac{1}{2} \varrho V^2 \end{equation*}$

where $\varrho$ is the density of the flowing fluid, with an equivalent height of the hydrostatic head on the manometer of $h_t$ .

Finally, consider the case where the reference pressure is now the static pressure in the flow. In this case, the pressure would be the dynamic pressure of the flow, i.e., the difference between the total pressure and the static pressure

(37) $\begin{equation*} p_s + \frac{1}{2} \varrho V^2 - p_s = \frac{1}{2} \varrho V^2 \end{equation*}$

with an equivalent hydrostatic head of $h_d$ on the manometer.

Pitot Tube & Pitot-Static Tube

Pitot probes (or Pitot tubes), static probes, and Pitot-static probes are often used to measure flow velocities. In each case, the principle is the same: pressure measurements can be used with the Bernoulli equation to calculate the flow velocity. The proviso is that the density of the flow where the probe is inserted is known. These three types of pressure probes are shown in the figure below.

When a Pitot tube and static tube are combined into a single unit, it is called a Pitot-static tube (or probe). A “static tube” *per se* is rarely used to measure static pressure in a flow.

For the analysis of a Pitot tube, consider point 1 upstream of the Pitot tube and point 2 at the entrance to the tube, as shown in the figure below. From Bernoulli’s equation it is known that $p+1/2\varrho V^{2} =$ constant along any given streamline, so that

(38) $\begin{equation*} p_{1} + \frac{1}{2} \varrho_{\infty} V_{1}^{2} = p_{2} + \frac{1}{2} \varrho V_{2}^{2} = p_T = p_{\infty} + \frac{1}{2} \varrho_{\infty} V_{\infty}^{2} \end{equation*}$

The difference in pressure between the total pressure and the static pressure can be measured on a suitably calibrated pressure gauge.

But at point 2 at the entrance to the Pitot tube, the fluid is brought to rest so $V_{2} = 0$ (i.e., it is a stagnation point), and so

(39) $\begin{equation*} p_1 + \frac{1}{2} \varrho_{\infty} V_{1}^2 = p_{2} = p_T =\mbox{\small Total pressure at point 2} \end{equation*}$

The static pressure $p_1 = p_{s} = p_{\infty}$ must be measured to obtain the flow velocity from this latter expression. Like the venturi problem, this static pressure can be measured with a separate static vent or with static ports on the outer side of a concentric Pitot-static tube. Solving for $V_{1}$ gives

(40) $\begin{equation*} V_{1} = \sqrt{\frac{ 2\left(p_{T} - p_{s}\right)}{\varrho_{\infty}}} = \sqrt{\frac{ 2\left(p_{T} - p_{\infty}\right)}{\varrho_{\infty}}} \end{equation*}$

Therefore, the upstream flow velocity can be measured by measuring the difference between a flow’s total and static pressure. Again, this pressure difference can be measured using a differential gauge, manometer, or pressure transducer.

The preceding principles are used in airspeed measurement using an Air Speed Indicator (ASI), which is part of the pneumatic system used on an aircraft for not only airspeed measurement but also altitude and rate of climb, as shown in the figure below. An airspeed indicator is a suitably calibrated differential pressure gauge. The dynamic pressure source needed for an airspeed indicator is either a Pitot tube or a Pitot-static tube. Notice, in this case, the use of a Pitot probe to measure the total pressure probe and a separate static vent (usually placed somewhere on the side of the fuselage) to measure the reference pressure. A heater prevents the Pitot probe from malfunctioning in icing conditions by preventing ice accumulation on the exposed probe.

A Pitot-static system on an airplane provides the reference pressures needed for the airspeed indicator and other pneumatic instruments. Notice, in this case, the use of a Pitot probe to measure total pressure and a separate static port to measure the static pressure.

More often than not, airplanes use Pitot tubes (to measure total pressure $p_T$ ) with static taps (to measure $p_s$ ) placed at some other point on the airframe, as previously shown. In this case, then the airspeed $V_{\infty}$ will be

(41) $\begin{equation*} V_{\infty} = \sqrt{\frac{ 2\left(p_{T} - p_{s}\right)}{\varrho_{\infty}}} \end{equation*}$

In either case, the purpose is the same: to obtain a measurement of the aircraft’s airspeed. This information is then provided to the pilot on an airspeed indicator. Therefore, an airspeed indicator is just a dynamic pressure gauge calibrated in units of speed; usually, units of nautical miles per hour (kts) are used, but sometimes miles per hour (mph). An online simulator of the Pitot-static system helps understand how it works and what happens under various failure scenarios.

Flow Speed in a Wind Tunnel

While it has been previously described how the flow speed in the test section of a wind tunnel can be obtained by measuring the static pressure drop between the inlet and the test section, another way is to use a Pitot probe, as shown in the figure below. The Pitot probe measures the total pressure, $p_T$ , in the upstream section or any upstream section that is convenient to do so. This approach tends to give a more accurate way of measuring dynamic pressure because it is the difference between a higher and a lower pressure rather than two lower pressures of similar magnitude.

Another method of measuring the flow velocity in a wind tunnel is to use measurements of total pressure and static pressure.

From the Bernoulli equation, then

(42) $\begin{equation*} P_T = p_2 + \dfrac{1}{2} \varrho V_2^2 = p_{\infty} + \dfrac{1}{2} \varrho V_{\infty}^2 \end{equation*}$

so that

(43) $\begin{equation*} V_2 = V_{\infty} = \sqrt{ \frac{2 \left( p_T - p_{\infty} \right)}{\varrho_{\infty}} } \end{equation*}$

Again, the value of density, $\varrho_{\infty}$ , can be obtained from static pressure and temperature measurements in conjunction with the equation of state. The calibration would be verified by placing a Pitot probe in the test section to obtain any calibration factor, $K_{\rm cal}$ , i.e.,

(44) $\begin{equation*} V_{\infty} = K_{\rm cal} \sqrt{ \frac{2 \left( p_T - p_{\infty} \right)}{\varrho_{\infty}} } \end{equation*}$

where $K_{\rm cal}$ will be very close to unity.

Jet Engine Performance

The thrust of a turbojet engine can be determined using conservation principles applied to a control volume surrounding the engine. There is a mass flow into the engine from the air and an additional mass flow from the fuel. Fuel is much denser than air, so although the volume flow of fuel may be relatively low, its mass flow is still significant and needs to be accounted for.

A control volume approach for analyzing a turbojet engine, which works on the air to increase its momentum in the downstream direction. According to Newton’s third law, the thrust is directed in the upstream direction.

The mass flow of air into the engine will be

(45) $\begin{equation*} \overbigdot{m}_{\rm air} = \varrho_{\infty} V_{\infty} A_i \end{equation*}$

where $A_i$ is the inlet area, and the mass flow rate of fuel is $\overbigdot{m}_{\rm fuel}$ . Therefore, the thrust is

(46) $\begin{equation*} T = (\overbigdot{m}_{\rm air} + \overbigdot{m}_{\rm fuel} ) V_e - \overbigdot{m}_{\rm air} V_{\infty} + ( p_e - p_{\infty} ) A_e ) \end{equation*}$

with $A_e$ as the exit area and $V_e$ as the exit or “jet” velocity. The pressure term in Eq. 46 is relatively small compared to the change in momentum of the flow and so may be neglected, i.e.,

(47) $\begin{equation*} T = (\overbigdot{m}_{\rm air} + \overbigdot{m}_{\rm fuel} ) V_e - \overbigdot{m}_{\rm air} V_{\infty} \end{equation*}$

If the mass flow of the fuel (small) is neglected in comparison to the mass flow rate of air (high), then

(48) $\begin{equation*} T = \overbigdot{m}_{\rm air} V_j - \overbigdot{m}_{\rm air} V_{\infty} = \overbigdot{m} \left(V_j - V_{\infty} \right) \end{equation*}$

Using the principle of conservation of energy, the propulsive efficiency can be defined as

(49) $\begin{equation*} \eta_p = \frac{\mbox{\small Useful power produced}}{\mbox{\small Total power expended }} = \frac{T \, V_{\infty}}{ T \, V_{\infty} + \frac{1}{2} \overbigdot{m} \left( V_j - V_{\infty} \right)^2} \end{equation*}$

which shows that losses appear as a gain in kinetic energy in the downstream jet flow. Substituting for the thrust gives

(50) $\begin{equation*} \eta_p= \frac{\overbigdot{m} \left( V_j - V_{\infty} \right) V_{\infty}}{ \overbigdot{m} \left( V_j - V_{\infty} \right) V_{\infty} + \frac{1}{2} \overbigdot{m} \left( V_j - V_{\infty} \right)^2} = \frac{2}{1 + \left( \displaystyle{\frac{V_j}{V_{\infty} }}\right) } \end{equation*}$

where to be meaningful $V_j > V_{\infty}$ . Notice that a low value of $V_j$ for a given thrust and, hence, higher propulsive efficiency can be achieved by having a significantly high mass flow rate, $\overbigdot{m}$ , through the engine. Therefore, even with this relatively simple analysis, it can be concluded that it is more efficient to create thrust by accelerating a large air volume (mass flow rate) at a lower $V_j$ versus a lower air volume (mass flow rate) at a high $V_j$ .

Forces on a Two-Dimensional body

The conservation laws can be used to determine the drag of a two-dimensional body in terms of off-surface fluid properties, e.g., the force on a body in a wind tunnel. Consider the physical problem of a body immersed in a fluid flow, as shown in the figure below. The fluid will generally exert pressure and viscous forces on the body, which are the sources of the resultant force, $R$ , and hence the lift, $L$ , and the drag, $D$ . Physically, because of these effects, a wake with a lower flow velocity will form downstream of the body. The objective is to use the equations of motion of the fluid (in this case, continuity and momentum) to calculate the force on the body.

The upstream and downstream sections labeled 1 and 2 and the top and bottom sides bound the control volume. The top and bottom sides of the control volume are considered far enough away from the body that they are in the free stream and, therefore, are streamlines to the flow, i.e., there will be no flow over these boundaries. These boundaries, for example, couldm be the floor and ceiling of a wind tunnel. Remember that $d\vec{S}$ points out of the control volume by convention. The upstream or free-stream (undisturbed) velocity is $V_1$ = constant = $V_{\infty}$ . Because of viscosity, a wake of reduced velocity or diminished flow exists behind the body; this gives a wake region a velocity profile in which $V_2 < V_\infty$ .

Consider the forces on the fluid in the control volume. They stem from two sources:

The net effect of the pressure distribution over the control volume, i.e., over the surface ABCD, then
(51) $\begin{equation*} \mbox{\small Net pressure force} = -\iint_{\rm ABCD} \, p \, d\vec{S} \end{equation*}$
The reaction of the surface forces on the body, i.e., those forces produced over the inner part of the control volume as created by the body, i.e., $-\vec{R}$ , in accordance with Newton’s 3rd law.

Remember that the flow will exert pressure and shear forces on the body, resulting in a resultant force, $\vec{R}$ . This means that by Newton’s third law, the airfoil will exert a force $-\vec{R}$ on the fluid within the control volume. Therefore, the total force on the fluid volume is the sum of the pressure forces on the control surface and the reaction force of the body on the fluid, i.e.,

(52) $\begin{equation*} \mbox{\small Net force} = -\iint_{\rm ABCD} p \, d\vec{S} - \vec{R} \end{equation*}$

and from the conservation of momentum equation, then

(53) $\begin{equation*} -\iint_{\rm ABCD} p \, d\vec{S} - \vec{R} = \frac{d}{dt}\oiiint_{\rm C.V.} \varrho\vec{V} d {\cal{V}} + \iint_{\rm ABCD} \left(\varrho \, \vec{V} \bigcdot d\vec{S} \right) \vec{V} \end{equation*}$

For steady flow $\partial ({\rm everything})/ \partial t = 0$ , so that this latter equation is simplified to

(54) $\begin{equation*} \vec{R} = -\iint_{\rm ABCD} \left(\varrho\vec{V}\bigcdot d\vec{S}\right) \vec{V} -\iint_{\rm ABCD} p \, d\vec{S} \end{equation*}$

Notice that this latter equation is a vector equation.

Drag Force

To find the drag, the $x$ -component of this latter equation is needed, noting that $V_{\infty}$ and $V_2$ are the inflow and outflow velocities in the $x$ -direction and that the $x$ -component of $R$ is the drag on the section, i.e., $D$ per unit span. The pressures upstream and downstream are equal, i.e., there is no longitudinal pressure gradient, so that

(55) $\begin{equation*} \iint_{\rm AD} p \, d\vec{S} = \iint_{\rm CB} p \, d\vec{S} \end{equation*}$

Therefore, the force on the fluid and drag on the body can be expressed only in terms of the ${u}$ velocity profiles upstream and downstream, i.e.,

(56) $\begin{equation*} D = -\iint_{\rm ABCD} \left(\varrho\vec{V}\bigcdot d\vec{S} \right) u \end{equation*}$

and so this integral becomes a one-dimensional integral (noting carefully the signs on the components), i.e.,

(57) $\begin{equation*} D = -\bigg( - \int_{D}^{A} \varrho_{_{\infty}} V_{\infty}^{2} dy + \int_{C}^{B}\varrho_{2}V_{2}^{2}dy \bigg) = \int_{D}^{A}\varrho_{_{\infty}}V_{\infty}^{2}dy - \int_{C}^{B}\varrho_{2} V_{2}^{2}dy \end{equation*}$

The continuity equation can also be applied to the control volume, i.e., based on the principle of conservation of mass, then

(58) $\begin{equation*} \int_{D}^{A}\varrho_{_{\infty}}V_{1}dy = \int_{C}^{B} \varrho_{2} V_{2} dy \end{equation*}$

Multiplying by $V_{\infty}$ in this case (which is a constant) gives

(59) $\begin{equation*} \int_{D}^{A}\varrho_{_{\infty}}V_{\infty}^2 dy = \int_{C}^{B} \varrho_{2} V_{2} V_{\infty} dy \end{equation*}$

Substituting back into the expression for $D$ gives

(60) $\begin{eqnarray*} D & = & \int_{D}^{A}\varrho_{_{\infty}}V_{\infty}^{2}dy - \int_{C}^{B}\varrho_{2} V_{2}^{2}dy = \int_{C}^{B} \varrho_{2} V_{2} V_{\infty}dy - \int_{C}^{B}\varrho_{2}V_{2}^{2}dy \\ [6pt] & = & \int_{C}^{B}\varrho_{2}V_{2} \left(V_{\infty} - V_{2}\right)dy \end{eqnarray*}$

This latter expression gives the drag of a two-dimensional body in terms of $V_{\infty}$ and the flow field properties $\varrho_2$ and $V_2$ across a vertical dimension placed downstream of the body.

Notice that $V_{\infty} - V_2$ is the velocity deficiency or decrement at a given station (or section) in the downstream wake. Also, $\varrho_{2}V_{2}$ is the mass flux, so the product $\varrho_{2}V_{2}(V_{\infty}-V_{2})$ gives the decrement in momentum. The integral of this expression leads to the total decrement in momentum behind the body and, hence, the drag of the section. If $\varrho$ = constant, then things can be simplified even further to get

(61) $\begin{equation*} D =\varrho\int_{C}^{B}V_{2}\left(V_{\infty} -V_{2}\right) dy =\varrho_{\infty} \int_{C}^{B}V_{2}\left(V_{\infty}-V_{2}\right) dy \end{equation*}$

However, some caution is appropriate because this approach works well only without significant flow separation and turbulence from the body. If there are significant losses in the wake that manifest as viscous losses from turbulent eddies, the drag will be underestimated.

This latter technique is frequently employed in wind tunnel testing for measuring forces on relatively streamlined bodies; the velocities downstream of the body or airfoil are measured using what is known as a wake rake. This rake is, in fact, an array of Pitot tubes, which are used to measure the dynamic pressure and the flow velocities in the wake. The momentum deficit is measured at discrete points in the wake, so that the drag per unit span is

(62) $\begin{equation*} D \approx \varrho_{_{\infty}} \sum_{i=1}^{N} V_2(x, y_i) \bigg( V_\infty - V_2(x, y_i) \bigg) \Delta y \end{equation*}$

where $V_2(x, y_i)$ is the velocity at the $i$ -th point in the wake, and $\Delta y$ is the spacing between the measurement points; $\Delta y$ does not have to be uniform. The values would be measured for each test point, e.g., different flow speeds and angles of attack.

Lift Force

Consider now how to measure the lift on the body, which is the force in the $y$ direction. The lift on the C.V. is calculated as the net vertical force resulting from the pressure difference between the suface AB (the ceiling) and CD (the floor), i.e.,

(63) $\begin{equation*} L = \iint_{_{\text{fl}}} p \, dS - \iint_{_{\text{ce}}} p \, dS = \iint_{\text{CD}} p \, dS - \iint_{\text{AB}} p \, dS \end{equation*}$

as illustrated in the figure below, where $p_{_{\text{ce}}}(x)$ is the measured pressure distribution (using pressure taps) along the ceiling, and $p_{_{\text{fl}}}(x)$ is the pressure distribution along the floor. Again, while the floor and ceiling are streamlines to the flow there can be a pressure distribution along their lengths from the presence of the body.

The lift on the body can be determined from the pressure distributions on the floor and ceiling of the wind tunnel.

The sectional lift per unit span is obtained using a one-dimensional pressure integral along the length of the control volume, i.e.,

(64) $\begin{equation*} L = \int_{x_1}^{x_2} \bigg( p_{_{\text{fl}}}(x) - p_{_{\text{ce}}}(x) \bigg) dx \end{equation*}$

$x_1$ would be located at the start of the test section and $x_2$ will be at or near the end. Each of the pressure taps would be connected to a pressure measurement system. Because the pressure will be measured at discrete points along the wind tunnel walls, the lift per unit span is given as

(65) $\begin{equation*} L \approx \sum_{i=1}^{N} \bigg( p_{_{\text{fl}}}(x_i) - p_{_{\text{ce}}}(x_i) \bigg) \Delta x \end{equation*}$

where $x_i$ represents the measurement locations along the wind tunnel, and $\Delta x$ is the spacing between measurement points; $\Delta x$ can be non-nuniform. This sum is just the area under the $\Delta$ pressure curve between floor and ceiling along the length of the test section.

Historical Context: NACA and Abbott & Von Doenhoff

The integral method was extensively utilized by NACA researchers, notably Ira Abbot and Albert von Doenhoff, whose influential work in the 1940s was published as “Theory of Wing Sections & Summary of Airfoil Data.” Their approach inferred forces from test section wall pressure distributions and wake surveys, eliminating the need for direct surface measurements like pressures on every airfoil. This method allowed over 200 airfoil sections to be tested and compared, significantly benefiting the aeronautics industry by providing a comprehensive catalog of airfoil shapes that could be selected to meet specific aerodynamic requirements.

Summary & Closure

Many fluid flow problems can be analyzed using the continuity and momentum equations in the integral form, as well as the energy equation in the surrogate form of the Bernoulli equation. All three conservation principles will be needed in most practical problems, along with the complete energy equation and the equation of state if compressibility effects are involved. The momentum equation will always be required whenever forces need to be calculated. Problems that can be assumed to be steady, inviscid, and incompressible are the easiest to understand and predict, and to this end, some exemplars have been discussed and solved in this chapter. The experience gained in solving these exemplary fluid flow problems increases confidence in analyzing more complex problems, such as those involving unsteadiness and compressibility effects in a flow.

5-Question Self-Assessment Quickquiz

For Further Thought or Discussion

Why is the flow through a rapidly converging or expanding venturi (or duct) more difficult to justify as one-dimensional?
What happens when an aircraft flies at higher Mach numbers? Does Bernoulli’s equation apply?
Who was Henri Pitot? What did he do besides invent the Pitot tube?
Using a drinking straw and a ruler, explain how you would measure the flow velocity in a water channel.
What happens in the throat of a venturi as the flow speed becomes subsonic and then supersonic?
Think about some aeronautical problems that cannot be tackled using the conservation principles in integral form.
It is proposed that the drag of a circular cylinder is to be measured in a wind tunnel using the momentum deficiency technique. What is the concern here, and why?
Why does a balloon filled with helium deflate quicker than a balloon filled with air?
If the drag on a body can be measured using a momentum integral method, how can the lift on the body be measured?

Additional Online Resources

Explore some of these additional online resources to help understand the application of the conservation laws when applied to simple fluid flow problems:

Explore this interactive simulation of flow through a pipe.
Video of worked example problems using the conservation laws.
View this video for a simple experimental demonstration of a U-tube manometer.
This Pitot Static System Simulator allows you to visualize the Pitot static system on an airplane under varying atmospheric conditions and what happens when system parts become blocked.

Sadler, D. R. (2005). "Interpretations of Criteria-Based Assessment and Grading in Higher Education." Assessment and Evaluation in Higher Education, 30 (2), 175–194. ↵
The diameter of such holes should be less than 0.5 mm or 20 thousandths of an inch. ↵
Evangelista Torricelli's original derivation can be found in the second book "De motu aquarum" of his "Opera Geometrica." See also: A. Malcherek, "History of the Torricelli Principle and a New Outflow Theory," Journal of Hydraulic Engineering, 142 (11), 1–7, 2016. ↵

License

Icon for the Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International License

Introduction to Aerospace Flight Vehicles Copyright © 2022–2024 by J. Gordon Leishman is licensed under a Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International License, except where otherwise noted.

Digital Object Identifier (DOI)

https://doi.org/https://doi.org/10.15394/eaglepub.2022.1066.n19