Isentropic nozzle flow

Last updated October 05, 2024

In fluid mechanics, isentropic nozzle flow describes the movement of a fluid through a narrow opening without an increase in entropy (an isentropic process).

Overview

Whenever a gas is forced through a tube, the gaseous molecules are deflected by the tube's walls. If the speed of the gas is much less than the speed of sound, the gas density will remain constant and the velocity of the flow will increase. However, as the speed of the flow approaches the speed of sound, compressibility effects on the gas are to be considered. The density of the gas becomes position dependent. While considering flow through a tube, if the flow is very gradually compressed (i.e. area decreases) and then gradually expanded (i.e. area increases), the flow conditions are restored (i.e. return to its initial position). So, such a process is a reversible process. According to the Second Law of Thermodynamics, whenever there is a reversible and adiabatic flow, constant value of entropy is maintained. Engineers classify this type of flow as an isentropic flow of fluids. Isentropic is the combination of the Greek word "iso" (which means - same) and entropy.

When the change in flow variables is small and gradual, isentropic flows occur. The generation of sound waves is an isentropic process. A supersonic flow that is turned while there is an increase in flow area is also isentropic. Since there is an increase in area, therefore we call this an isentropic expansion. If a supersonic flow is turned abruptly and the flow area decreases, the flow is irreversible due to the generation of shock waves. The isentropic relations are no longer valid and the flow is governed by the oblique or normal shock relations.

Set of Equations

Below are nine equations commonly used when evaluating isentropic flow conditions.^[1] These assume the gas is calorically perfect; i.e. the ratio of specific heats is a constant across the temperature range. In typical cases the actual variation is only slight.

$M={\frac {v}{c_{s}}}$
$c_{s}={\sqrt {\frac {\gamma p}{\rho }}}={\sqrt {\gamma RT}}$
${\frac {p}{\rho ^{\gamma }}}={\frac {p_{t}}{\rho _{t}^{\gamma }}}$
${\frac {p}{p_{t}}}=\left({\frac {\rho }{\rho _{t}}}\right)^{\gamma }=\left({\frac {T}{T_{t}}}\right)^{\frac {\gamma }{\gamma -1}}$
$q={\frac {1}{2}}\rho v^{2}={\frac {1}{2}}\gamma pM^{2}$
${\frac {p}{p_{t}}}=\left(1+{\frac {\gamma -1}{2}}M^{2}\right)^{\frac {-\gamma }{\gamma -1}}$
${\frac {T}{T_{t}}}=\left(1+{\frac {\gamma -1}{2}}M^{2}\right)^{-1}$
${\frac {\rho }{\rho _{t}}}=\left(1+{\frac {\gamma -1}{2}}M^{2}\right)^{\frac {-1}{\gamma -1}}$
${\frac {A}{A^{*}}}={\frac {1}{M}}\left({\frac {\gamma -1}{2}}\right)^{\frac {-(\gamma +1)}{2(\gamma -1)}}\left(1+{\frac {\gamma -1}{2}}M^{2}\right)^{\frac {-\gamma }{\gamma -1}}$

Properties without a subscript are evaluated at the point of interest (this point may be chosen anywhere along the length of the nozzle, but once chosen, all properties in a calculation must be evaluated at the same point)
Subscript $t$ denotes a property at total/stagnation conditions. In a rocket or jet engine, this means the conditions inside the combustion chamber. For example, $p_{t}$ is total pressure/stagnation pressure/chamber pressure (all equivalent).
$M$ is the local Mach number of the gas
$v$ is the speed of the gas (m/s)
$c_{s}$ is the local speed of sound through the gas (m/s)
$\gamma$ is the ratio of specific heats of the gas
$p$ is the pressure of the gas (Pa)
$\rho$ is the density of the gas (kg/m³)
$T$ is the temperature of the gas (K)
$A$ is the cross sectional area of the nozzle at the point of interest (m²)
$A^{*}$ is the cross sectional area of the nozzle at the sonic point, or the point where gas velocity is Mach 1 (m²). Ideally this will occur at the nozzle throat.

Stagnation properties

In fluid dynamics, a stagnation point is a point in a flow field where the local velocity of the fluid is zero. The isentropic stagnation state is the state a flowing fluid would attain if it underwent a reversible adiabatic deceleration to zero velocity. There are both actual and the isentropic stagnation states for a typical gas or vapor. Sometimes it is advantageous to make a distinction between the actual and the isentropic stagnation states. The actual stagnation state is the state achieved after an actual deceleration to zero velocity (as at the nose of a body placed in a fluid stream), and there may be irreversibility associated with the deceleration process. Therefore, the term "stagnation property" is sometimes reserved for the properties associated with the actual state, and the term total property is used for the isentropic stagnation state. The enthalpy is the same for both the actual and isentropic stagnation states (assuming that the actual process is adiabatic). Therefore, for an ideal gas, the actual stagnation temperature is the same as the isentropic stagnation temperature. However, the actual stagnation pressure may be less than the isentropic stagnation pressure. For this reason the term "total pressure" (meaning isentropic stagnation pressure) has particular meaning compared to the actual stagnation pressure.

Flow analysis

The isentropic efficiency is ${\frac {h_{1}-h_{2a}}{h_{1}-h_{2}}}$ . The variation of fluid density for compressible flows requires attention to density and other fluid property relationships. The fluid equation of state, often unimportant for incompressible flows, is vital in the analysis of compressible flows. Also, temperature variations for compressible flows are usually significant and thus the energy equation is important. Curious phenomena can occur with compressible flows.

For simplicity, the gas is assumed to be an ideal gas.
The gas flow is isentropic.
The gas flow is constant.
The gas flow is along a straight line from gas inlet to exhaust gas exit.
The gas flow behavior is compressible.

There are numerous applications where a steady, uniform, isentropic flow is a good approximation to the flow in conduits. These include the flow through a jet engine, through the nozzle of a rocket, from a broken gas line, and past the blades of a turbine.

m

= Mach number

V

= velocity

R

= universal gas constant

p

= pressure

k

= specific heat ratio

T

= temperature

* = sonic conditions

ρ

= density

A

= area

M m

= molar mass

Energy equation for the steady flow: $q_{net}+h+{\frac {V^{2}}{2}}=w_{net}+h_{0}+{\frac {V_{0}^{2}}{2}}$

To model such situations, consider the control volume in the changing area of the conduit of Fig. The continuity equation between two sections an infinitesimal distance $dx$ apart is $\rho AV=(\rho +d\rho )(A+dA)(V+dV)$

If only the first-order terms in a differential quantity are retained, continuity takes the form ${\frac {dV}{V}}+{\frac {dA}{A}}+{\frac {d\rho }{\rho }}=0$

The energy equation is: ${\frac {V^{2}}{2}}+{\frac {k}{k-1}}\cdot {\frac {p}{\rho }}={\frac {\left(V+dV\right)^{2}}{2}}+{\frac {k}{k-1}}\cdot {\frac {p+dp}{\rho +d\rho }}$

This simplifies to, neglecting higher-order terms: $VdV+{\frac {k}{k-1}}\cdot {\frac {\rho dp-pd\rho }{\rho ^{2}}}=0$ Assuming an isentropic flow, the energy equation becomes: $VdV+{\frac {k\cdot p}{\rho ^{2}}}\cdot d\rho =0$ Substitute from the continuity equation to obtain ${\frac {dV}{V}}\cdot \left({\frac {\rho V^{2}}{kp}}-1\right)={\frac {dA}{A}}$ or, in terms of the Mach number: ${\frac {dV}{V}}\cdot \left(M^{2}-1\right)={\frac {dA}{A}}$

This equation applies to a steady, uniform, isentropic flow. There are several observations that can be made from an analysis of Eq. (9.26). They are:

For a subsonic flow in an expanding conduit ( $M < 1$ and $dA > 0$ ), the flow is decelerating ( $dV < 0$ ).
For a subsonic flow in a converging conduit ( $M < 1$ and $dA < 0$ ), the flow is accelerating ( $dV > 0$ ).
For a supersonic flow in an expanding conduit ( $M > 1$ and $dA > 0$ ), the flow is accelerating ( $dV > 0$ ).
For a supersonic flow in a converging conduit ( $M > 1$ and $dA < 0$ ), the flow is decelerating ( $dV < 0$ ).
At a throat where $dA = 0$ , either $M = 1$ or $dV = 0$ (the flow could be accelerating through $M = 1$ , or it may reach a velocity such that $dV = 0$ ).

Supersonic flow

A nozzle for a supersonic flow must increase in area in the flow direction, and a diffuser must decrease in area, opposite to a nozzle and diffuser for a subsonic flow. So, for a supersonic flow to develop from a reservoir where the velocity is zero, the subsonic flow must first accelerate through a converging area to a throat, followed by continued acceleration through an enlarging area.

The nozzles on a rocket designed to place satellites in orbit are constructed using such converging-diverging geometry. The energy and continuity equations can take on particularly helpful forms for the steady, uniform, isentropic flow through the nozzle. Apply the energy equation with $Q, W S = 0$ between the reservoir and some location in the nozzle to obtain $c_{p}\cdot T_{0}={\frac {V^{2}}{2}}+c_{p}\cdot T$

Any quantity with a zero subscript refers to a stagnation point where the velocity is zero, such as in the reservoir. Using several thermodynamic relations, equations can be put in the forms:

${\begin{aligned}{\frac {T_{0}}{T}}&=1+{\frac {k-1}{2}}M^{2}\\{\frac {p_{0}}{p}}&=\left(1+{\frac {k-1}{2}}M^{2}\right)^{\frac {k}{k-1}}\\{\frac {\rho _{0}}{\rho }}&=\left(1+{\frac {k-1}{2}}M^{2}\right)^{\frac {1}{k-1}}\end{aligned}}$

If the above equations are applied at the throat (the critical area signified by an asterisk (*) superscript, where $M = 1$ ), the energy equation takes the forms

${\begin{aligned}{\frac {T^{*}}{T_{0}}}&={\frac {2}{k+1}}\\{\frac {p^{*}}{p_{0}}}&=\left({\frac {2}{k+1}}\right)^{\frac {k}{k-1}}\\{\frac {\rho ^{*}}{\rho _{0}}}&=\left({\frac {2}{k+1}}\right)^{\frac {1}{k-1}}\end{aligned}}$

The critical area is often referenced even though a throat does not exist. For air with $k = 1.4$ , the equations above provide ${\begin{aligned}T^{*}&=0.833333\cdot T_{0}\\p^{*}&=0.528282\cdot p_{0}\\\rho ^{*}&=0.633938\cdot \rho _{0}\end{aligned}}$

The mass flux through the nozzle is of interest and is given by: ${\dot {m}}=\rho AV={\frac {p}{RT}}\cdot A\cdot M{\sqrt {kR_{specific}T}}=p{\sqrt {kM_{m} \over RT}}AM$

With the use of Eq. (9.28), the mass flux, after applying some algebra, can be expressed as ${\dot {m}}=p_{0}MA{\sqrt {\frac {kM_{m}}{RT_{0}}}}\left(1+{\frac {k-1}{2}}M^{2}\right)^{\frac {k+1}{-2(k-1)}}$

If the critical area is selected where $M = 1$ , this takes the form ${\dot {m}}=p_{0}A^{*}{\sqrt {\frac {kM_{m}}{RT_{0}}}}\left(1+{\frac {k-1}{2}}\right)^{\frac {k+1}{-2(k-1)}}$ which, when combined with previous it provides: ${\frac {A}{A^{*}}}={\frac {1}{M}}\left({\frac {k+1}{2+(k-1)\cdot M^{2}}}\right)^{\frac {k+1}{-2(k-1)}}$

Converging nozzle

Consider a converging nozzle connecting a reservoir with a receiver. If the reservoir pressure is held constant and the receiver pressure reduced, the Mach number at the exit of the nozzle will increase until $M e = 1$ is reached, indicated by the left curve in figure 2. After $M e = 1$ is reached at the nozzle exit for $p r = 0.5283 p 0$ , the condition of choked flow occurs and the velocity throughout the nozzle cannot change with further decreases in $p r$ . This is due to the fact that pressure changes downstream of the exit cannot travel upstream to cause changes in the flow conditions. The right curve of figure 2. represents the case when the reservoir pressure is increased and the receiver pressure is held constant. When $M e = 1$ , the condition of choked flow also occurs; but Eq indicates that the mass flux will continue to increase as $p 0$ is increased. This is the case when a gas line ruptures.

It is interesting that the exit pressure $p e$ is able to be greater than the receiver pressure $p r$ . Nature allows this by providing the streamlines of a gas the ability to make a sudden change of direction at the exit and expand to a much greater area resulting in a reduction of the pressure from $p e$ to $p r$ . The case of a converging-diverging nozzle allows a supersonic flow to occur, providing the receiver pressure is sufficiently low. This is shown in figure 3 assuming a constant reservoir pressure with a decreasing receiver pressure. If the receiver pressure is equal to the reservoir pressure, no flow occurs, represented by curve A. If $p r$ is slightly less than $p 0$ , the flow is subsonic throughout, with a minimum pressure at the throat, represented by curve B. As the pressure is reduced still further, a pressure is reached that result in $M = 1$ at the throat with subsonic flow throughout the remainder of the nozzle. There is another receiver pressure substantially below that of curve C that also results in isentropic flow throughout the nozzle, represented by curve D; after the throat the flow is supersonic. Pressures in the receiver in between those of curve C and curve D result in non-isentropic flow (a shock wave occurs in the flow). If $p r$ is below that of curve D, the exit pressure $p e$ is greater than $p r$ . Once again, for receiver pressures below that of curve C, the mass flux remains constant since the conditions at the throat remain unchanged. It may appear that the supersonic flow will tend to separate from the nozzle, but just the opposite is true. A supersonic flow can turn very sharp angles, since nature provides expansion fans that do not exist in subsonic flows. To avoid separation in subsonic nozzles, the expansion angle should not exceed 10°. For larger angles, vanes are used so that the angle between the vanes does not exceed 10°.

Related Research Articles

The Mach number, often only Mach, is a dimensionless quantity in fluid dynamics representing the ratio of flow velocity past a boundary to the local speed of sound. It is named after the Austrian physicist and philosopher Ernst Mach.

Compressible flow is the branch of fluid mechanics that deals with flows having significant changes in fluid density. While all flows are compressible, flows are usually treated as being incompressible when the Mach number is smaller than 0.3. The study of compressible flow is relevant to high-speed aircraft, jet engines, rocket motors, high-speed entry into a planetary atmosphere, gas pipelines, commercial applications such as abrasive blasting, and many other fields.

An isentropic process is an idealized thermodynamic process that is both adiabatic and reversible. The work transfers of the system are frictionless, and there is no net transfer of heat or matter. Such an idealized process is useful in engineering as a model of and basis of comparison for real processes. This process is idealized because reversible processes do not occur in reality; thinking of a process as both adiabatic and reversible would show that the initial and final entropies are the same, thus, the reason it is called isentropic. Thermodynamic processes are named based on the effect they would have on the system. Even though in reality it is not necessarily possible to carry out an isentropic process, some may be approximated as such.

In fluid dynamics, the Euler equations are a set of partial differential equations governing adiabatic and inviscid flow. They are named after Leonhard Euler. In particular, they correspond to the Navier–Stokes equations with zero viscosity and zero thermal conductivity.

The Rankine–Hugoniot conditions, also referred to as Rankine–Hugoniot jump conditions or Rankine–Hugoniot relations, describe the relationship between the states on both sides of a shock wave or a combustion wave in a one-dimensional flow in fluids or a one-dimensional deformation in solids. They are named in recognition of the work carried out by Scottish engineer and physicist William John Macquorn Rankine and French engineer Pierre Henri Hugoniot.

A de Laval nozzle is a tube which is pinched in the middle, making a carefully balanced, asymmetric hourglass shape. It is used to accelerate a compressible fluid to supersonic speeds in the axial (thrust) direction, by converting the thermal energy of the flow into kinetic energy. De Laval nozzles are widely used in some types of steam turbines and rocket engine nozzles. It also sees use in supersonic jet engines.

In fluid dynamics, the pressure coefficient is a dimensionless number which describes the relative pressures throughout a flow field. The pressure coefficient is used in aerodynamics and hydrodynamics. Every point in a fluid flow field has its own unique pressure coefficient, $C p$ .

An orifice plate is a device used for measuring flow rate, for reducing pressure or for restricting flow.

An oblique shock wave is a shock wave that, unlike a normal shock, is inclined with respect to the direction of incoming air. It occurs when a supersonic flow encounters a corner that effectively turns the flow into itself and compresses. The upstream streamlines are uniformly deflected after the shock wave. The most common way to produce an oblique shock wave is to place a wedge into supersonic, compressible flow. Similar to a normal shock wave, the oblique shock wave consists of a very thin region across which nearly discontinuous changes in the thermodynamic properties of a gas occur. While the upstream and downstream flow directions are unchanged across a normal shock, they are different for flow across an oblique shock wave.

In fluid dynamics, Rayleigh flow refers to frictionless, non-adiabatic fluid flow through a constant-area duct where the effect of heat transfer is considered. Compressibility effects often come into consideration, although the Rayleigh flow model certainly also applies to incompressible flow. For this model, the duct area remains constant and no mass is added within the duct. Therefore, unlike Fanno flow, the stagnation temperature is a variable. The heat addition causes a decrease in stagnation pressure, which is known as the Rayleigh effect and is critical in the design of combustion systems. Heat addition will cause both supersonic and subsonic Mach numbers to approach Mach 1, resulting in choked flow. Conversely, heat rejection decreases a subsonic Mach number and increases a supersonic Mach number along the duct. It can be shown that for calorically perfect flows the maximum entropy occurs at $M = 1$ .

In fluid dynamics, stagnation pressure is the static pressure at a stagnation point in a fluid flow. At a stagnation point the fluid velocity is zero. In an incompressible flow, stagnation pressure is equal to the sum of the free-stream static pressure and the free-stream dynamic pressure.

Choked flow is a compressible flow effect. The parameter that becomes "choked" or "limited" is the fluid velocity.

Subsonic wind tunnels are used for operations at low Mach numbers, with speeds in the test section up to 480 km/h. They may be of open-return type or closed-return flow. These tunnels use large axial fans to move air and increase dynamic pressure, overcoming viscous losses. The design principles of subsonic wind tunnels are based on the continuity equation and Bernoulli's principle, which allow for the calculation of important parameters such as the tunnel's contraction ratio.

Accidental release source terms are the mathematical equations that quantify the flow rate at which accidental releases of liquid or gaseous pollutants into the ambient environment which can occur at industrial facilities such as petroleum refineries, petrochemical plants, natural gas processing plants, oil and gas transportation pipelines, chemical plants, and many other industrial activities. Governmental regulations in many countries require that the probability of such accidental releases be analyzed and their quantitative impact upon the environment and human health be determined so that mitigating steps can be planned and implemented.

In fluid dynamics, dynamic pressure is the quantity defined by:

A supersonic expansion fan, technically known as Prandtl–Meyer expansion fan, a two-dimensional simple wave, is a centered expansion process that occurs when a supersonic flow turns around a convex corner. The fan consists of an infinite number of Mach waves, diverging from a sharp corner. When a flow turns around a smooth and circular corner, these waves can be extended backwards to meet at a point.

In fluid dynamics, Fanno flow is the adiabatic flow through a constant area duct where the effect of friction is considered. Compressibility effects often come into consideration, although the Fanno flow model certainly also applies to incompressible flow. For this model, the duct area remains constant, the flow is assumed to be steady and one-dimensional, and no mass is added within the duct. The Fanno flow model is considered an irreversible process due to viscous effects. The viscous friction causes the flow properties to change along the duct. The frictional effect is modeled as a shear stress at the wall acting on the fluid with uniform properties over any cross section of the duct.

The Cauchy number (Ca) is a dimensionless number in continuum mechanics used in the study of compressible flows. It is named after the French mathematician Augustin Louis Cauchy. When the compressibility is important the elastic forces must be considered along with inertial forces for dynamic similarity. Thus, the Cauchy Number is defined as the ratio between inertial and the compressibility force in a flow and can be expressed as

Non ideal compressible fluid dynamics (NICFD), or non ideal gas dynamics, is a branch of fluid mechanics studying the dynamic behavior of fluids not obeying ideal-gas thermodynamics. It is for example the case of dense vapors, supercritical flows and compressible two-phase flows. With the term dense vapors, we indicate all fluids in the gaseous state characterized by thermodynamic conditions close to saturation and the critical point. Supercritical fluids feature instead values of pressure and temperature larger than their critical values, whereas two-phase flows are characterized by the simultaneous presence of both liquid and gas phases.

Taylor–Maccoll flow refers to the steady flow behind a conical shock wave that is attached to a solid cone. The flow is named after G. I. Taylor and J. W. Maccoll, whom described the flow in 1933, guided by an earlier work of Theodore von Kármán.

References

↑ "Isentropic Flow Equations". www.grc.nasa.gov. Retrieved 2023-08-31.

Colbert, Elton J. Isentropic Flow Through Nozzles . University of Nevada, Reno. 3 May 2001. Accessed 15 July 2014.
Benson, Tom. "Isentropic Flow". NASA.gov. National Aeronautics and Space Administration. 21 June 2014. Accessed 15 July 2014.
Bar-Meir, Genick. "Isenotropic Flow". Archived 2015-02-26 at the Wayback Machine Potto.org. Potto Project. 21 November 2007. Accessed 15 July 2014.

This page is based on this Wikipedia article
Text is available under the CC BY-SA 4.0 license; additional terms may apply.
Images, videos and audio are available under their respective licenses.

[1] "Isentropic Flow Equations". www.grc.nasa.gov. Retrieved 2023-08-31.

[1]