Burt Rutan's White Knight and SpaceShip One, Photo Courtesy of Scaled Composites
Thermodynamics and Propulsion

6.7 Examples of Lost Work in Engineering Processes

  1. Lost work in Adiabatic Throttling: Entropy and Stagnation Pressure Changes
    Figure 6.8: Adiabatic Throttling
    Image fig3AdiabaticThrottling_web

    A process we have encountered before is adiabatic throttling of a gas, by a valve or other device as shown in Figure 6.8. The velocity is denoted by $ c$ . There is no shaft work and no heat transfer and the flow is steady. Under these conditions we can use the first law for a control volume (the Steady Flow Energy Equation) to make a statement about the conditions upstream and downstream of the valve:

    $\displaystyle h_1 +\frac{c_1^2}{2} =h_2 +\frac{c_2^2}{2} =h_t,$

    where $ h_t$ is the stagnation enthalpy, corresponding to a (possibly fictitious) state with zero velocity. The stagnation enthalpy is the same at stations 1 and 2 if $ Q=W=0$ , even if the flow processes are not reversible.

    For a perfect gas with constant specific heats, $ h = c_p T$ and $ h_t
=c_p T_t$ . The relation between the static and stagnation temperatures is:

    $\displaystyle \frac{Tt}{T}$ $\displaystyle = 1 + \frac{c^2}{2c_pT} = 1 + \frac{(\gamma-1)c^2}{2\underbrace{\gamma RT}_{a^2}} =1+\frac{(\gamma-1)c^2}{2a^2},$    
    $\displaystyle \frac{T_t}{T}$ $\displaystyle =1+\left(\frac{\gamma-1}{2}\right)M^2,$    

    where $ a$ is the speed of sound and $ M$ is the Mach number, $ M = c / a$ . In deriving this result, use has only been made of the first law, the equation of state, the speed of sound, and the definition of the Mach number. Nothing has yet been specified about whether the process of stagnating the fluid is reversible or irreversible.

    When we define the stagnation pressure, however, we do it with respect to isentropic deceleration to the zero velocity state. For an isentropic process

    $\displaystyle \frac{P_2}{P_1}=\left(\frac{T_2}{T_1}\right)^{\gamma/(\gamma-1)}.$

    The relation between static and stagnation pressures is

    $\displaystyle \frac{P_t}{P}=\left(\frac{T_t}{T}\right)^{\gamma/(\gamma-1)}.$

    Figure 6.9: Static and stagnation pressures and temperatures
    Image fig3StaticAndStagnationTemperaturesTS_web

    The stagnation state is defined by $ P_t$ , $ T_t$ . In addition, $ s_\textrm{stagnation state} =s_\textrm{static state}$ . The static and stagnation states are shown in $ T$ -$ s$ coordinates in Figure 6.9.

    Stagnation pressure is a key variable in propulsion and power systems. To see why, we examine the relation between stagnation pressure, stagnation temperature, and entropy. The form of the combined first and second law that uses enthalpy is

    $\displaystyle Tds =dh - \frac{1}{\rho}dP.$ (6..8)

    Figure 6.10: Stagnation and static states
    Image fig3StagnationAndStaticStates_web

    This holds for small changes between any thermodynamic states and we can apply it to a situation in which we consider differences between stagnation states, say one state having properties $ (T_t,\;P_t)$ and the other having properties $ (T_t+dT_t,\;P_t +dP_t)$ (see Figure 6.10). The corresponding static states are also indicated. Because the entropy is the same at static and stagnation conditions, $ ds$ needs no subscript. Writing (6.8) in terms of stagnation conditions yields

    $\displaystyle ds = \frac{c_pdT_t}{T_t} - \frac{1}{\rho_tT_t}dP_t = \frac{c_p dT_t}{T_t}-\frac{R}{P_t}dP_t.$

    Both sides of the above are perfect differentials and can be integrated as

    $\displaystyle \frac{\Delta s}{R}=\frac{\gamma}{\gamma-1}\ln\left(\frac{T_{t2}}{T_{t1}}\right)
-\ln\left(\frac{P_{t2}}{P_{t1}}\right).$

    For a process with $ Q=W=0$ , the stagnation enthalpy, and hence the stagnation temperature, is constant. In this situation, the stagnation pressure is related directly to the entropy as,

    $\displaystyle \frac{\Delta s}{R} = -\ln\left(\frac{P_{t2}}{P_{t1}}\right).$ (6..9)

    Figure 6.11: Losses reflected in changes in stagnation pressure when $ T_t = \textrm {const.}$
    Image fig3StagnationPressureAndDeltaS_web

    Figure 6.11 shows this relation on a $ T$ -$ s$ diagram. We have seen that the entropy is related to the loss, or irreversibility. The stagnation pressure plays the role of an indicator of loss if the stagnation temperature is constant. The utility is that it is the stagnation pressure (and temperature) which are directly measured, not the entropy. The throttling process is a representation of flow through inlets, nozzles, stationary turbomachinery blades, and the use of stagnation pressure as a measure of loss is a practice that has widespread application.

    Equation (6.9) can be put in several useful approximate forms. First, we note that for aerospace applications we are (hopefully!) concerned with low loss devices, so that the stagnation pressure change is small compared to the inlet level of stagnation pressure,

    $\displaystyle \frac{\Delta P_t}{P_{t1}} = \frac{(P_{t1} -P_{t2})}{P_{t1}}\ll 1.$

    Expanding the logarithm (using $ \ln(1-x) \cong -x + \ldots$ ),

    $\displaystyle \ln\left(\frac{P_{t2}}{P_{t1}}\right) =\ln\left(1-\frac{\Delta P_t}{P_{t1}}\right)
\approx -\frac{\Delta P_t}{P_{t1}},$

    or

    $\displaystyle \frac{\Delta s}{R} \approx \frac{\Delta P_t}{P_{t1}}.$

    Another useful form is obtained by dividing both sides by $ c^2/2$ and taking the limiting forms of the expression for stagnation pressure in the limit of low Mach number ($ M\ll 1$ ). Doing this, we find:

    $\displaystyle \frac{T\Delta s}{c^2/2} \cong \frac{\Delta P_t}{\rho c^2/2}.$

    The quantity on the right can be interpreted as the change in the ``Bernoulli constant'' for incompressible (low Mach number) flow. The quantity on the left is a non-dimensional entropy change parameter, with the term $ T \Delta s$ now representing the loss of mechanical energy associated with the change in stagnation pressure.

    To summarize:

    1. For many applications the stagnation temperature is constant and the change in stagnation pressure is a direct measure of the entropy increase.
    2. Stagnation pressure is the quantity that is actually measured so that linking it to entropy (which is not measured) is useful.
    3. We can regard the throttling process as a ``free expansion'' at constant temperature $ T_{t1}$ from the initial stagnation pressure to the final stagnation pressure. We thus know that, for the process, the work we need to do to bring the gas back to the initial state is $ T_t \Delta s$ , which is the ``lost work'' per unit mass.





    Muddy Points

    Why do we find stagnation enthalpy if the velocity never equals zero in the flow? (MP 6.13)

    Why does $ T_t$ remain constant for throttling? (MP 6.14)

  2. Adiabatic Efficiency of a Propulsion System Component (Turbine)

    Figure 6.12: Schematic of turbine and associated thermodynamic representation in $ h$ -$ s$ coordinates
    Image fig3SchematicOfTurbine_web

    A schematic of a turbine and the accompanying thermodynamic diagram are given in Figure 6.12. There is a pressure and temperature drop through the turbine and it produces work. There is no heat transfer so the expressions that describe the overall shaft work and the shaft work per unit mass are


    If the difference in the kinetic energy at inlet and outlet can be neglected, Equation (6.11) reduces to

    $\displaystyle (h_2 - h_1)= w_\textrm{shaft}.$

    The adiabatic efficiency of the turbine is defined as

    $\displaystyle \eta_\textrm{ad} = \left[\frac{\textrm{actual work}}
{\textrm{ideal work $(\Delta s=0)$}}\right]_\textrm{for a given
pressure ratio}.$

    The performance of the turbine can be represented in an $ h$ -$ s$ plane (similar to a $ T$ -$ s$ plane for a perfect gas with constant specific heats) as shown in Figure 6.12. From the figure the adiabatic efficiency is

    $\displaystyle \eta_\textrm{ad} = \frac{h_1 - h_2}{h_1-h_{2s}} = \frac{h_1-h_{2s}-(h_2-h_{2s})}{h_1-h_{2s}}.$

    The adiabatic efficiency can therefore be written as

    $\displaystyle \eta_\textrm{ad} = 1 - \left(\frac{\Delta h}{\textrm{Ideal work}}\right).$

    The non-dimensional term $ (\Delta h/\textrm{Ideal work})$ represents the departure from isentropic (reversible) processes and hence a loss. The quantity $ \Delta h$ is the enthalpy difference for two states along a constant pressure line (see diagram). From the combined first and second laws, for a constant pressure process, small changes in enthalpy are related to the entropy change by $ Tds
=dh$ , or approximately,

    $\displaystyle T_2 \Delta s = \Delta h.$

    The adiabatic efficiency can thus be approximated as

    $\displaystyle \eta_\textrm{ad} = 1 - \frac{T_2\Delta s}{h_1-h_{2s}}=
1 - \left(\frac{\textrm{Lost work}}{\textrm{Ideal work}}\right).$

    The quantity $ T \Delta s$ represents a useful figure of merit for fluid machinery inefficiency due to irreversibility.





    Muddy Points

    How do you tell the difference between shaft work/power and flow work in a turbine, both conceptually and mathematically? (MP 6.15)

  3. Isothermal Expansion with Friction

    Figure 6.13: Isothermal expansion with friction
    Image fig3IsothermalExpansionWithFriction_web

    In a more general look at the isothermal expansion, we now drop the restriction to frictionless processes. As seen in Figure 6.13, work is done to overcome friction. If the kinetic energy of the piston is negligible, a balance of forces tells us that

    $\displaystyle W_\textrm{system on piston} =W_\textrm{done by friction} +W_\textrm{received}.$

    During the expansion, the piston and the walls of the container will heat up because of the friction. The heat will be (eventually) transferred to the atmosphere; all frictional work ends up as heat transferred to the surrounding atmosphere.

    $\displaystyle W_\textrm{friction} =Q_\textrm{friction}.$

    The amount of heat transferred to the atmosphere due to the frictional work only is thus,

    $\displaystyle Q_\textrm{friction} =\underbrace{W_\textrm{system on piston}}_\textrm{Work produced}
-\underbrace{W_\textrm{received}}_\textrm{work received}.$

    The entropy change of the atmosphere (considered as a heat reservoir) due to the frictional work is Image fig3IrreversibleCycle_web

    The engine operates in a cycle and the entropy change for the complete cycle is zero (because entropy is a state variable). Therefore,

    $\displaystyle \Delta S = 0 + \underbrace{\Delta S_\textrm{heat source} +
\Delta S_\textrm{heat sink}}_{\Delta S_\textrm{surroundings}}.$

    The total entropy change is

    $\displaystyle \Delta S_\textrm{total} =\Delta S_\textrm{heat source} +\Delta S_\textrm{heat sink} =
-\frac{Q}{T} + \frac{Q_0}{T_0}.$

    Suppose we had an ideal reversible engine working between these same two temperatures, which extracted the same amount of heat, $ Q$ , from the high temperature reservoir, and rejected heat of magnitude $ Q_\textrm{0,rev}$ to the low temperature reservoir. The work done by this reversible engine is

    $\displaystyle W_\textrm{rev} =Q-Q_\textrm{0,rev}.$

    For the reversible engine the total entropy change over a cycle is

    $\displaystyle \Delta S_\textrm{total} =\Delta S_\textrm{heat source} +\Delta S_\textrm{heat sink} =
-\frac{Q}{T} + \frac{Q_\textrm{0,rev}}{T_0} = 0.$

    Combining the expressions for work and for the entropy changes,

    $\displaystyle Q_0=Q_\textrm{0,rev} +W_\textrm{rev} -W.$

    The entropy change for the irreversible cycle can therefore be written as

    $\displaystyle \Delta S_\textrm{total} = \underbrace{-\frac{Q}{T} + \frac{Q_\textrm{0,rev}}{T_0}}_{=0}
+\frac{W_\textrm{rev}-W}{T_0}.$

    The difference in work that the two cycles produce is proportional to the entropy that is generated during the cycle:

    $\displaystyle T_0\Delta S_\textrm{total} =W_\textrm{rev} -W.$

    The second law states that the total entropy generated is greater than zero for an irreversible process, so that the reversible work is greater than the actual work of the irreversible cycle.

    An ``engine effectiveness,'' $ E_\textrm{engine}$ , can be defined as the ratio of the actual work obtained divided by the work that would have been delivered by a reversible engine operating between the two temperatures $ T$ and $ T_0$ :

    $\displaystyle E_\textrm{engine}$ $\displaystyle =\frac{\eta_\textrm{engine}}{\eta_{\textrm{reversible engine}}}$    
      $\displaystyle =\frac{W}{W_\textrm{rev}} =\frac{\textrm{Actual work obtained}}{\...
...{Work that \emph{would be delivered} by a reversible cycle between $T$, $T_0$}}$    
    $\displaystyle E_\textrm{engine}$ $\displaystyle =\frac{W_\textrm{rev}-T_0\Delta S_\textrm{total}}{W_\textrm{rev}}= 1-\frac{T_0\Delta S_\textrm{total}}{W_\textrm{rev}}.$    

    The departure from a reversible process is directly reflected in the entropy change and the decrease in engine effectiveness.





    Muddy Points

    Why does $ \Delta S_\textrm{irrev}=\Delta S_\textrm{total}$ ? (MP 6.17)

    In discussing the terms ``closed system'' and ``isolated system,'' can you assume that you are discussing a cycle or not? (MP 6.18)

    Does a cycle process have to have $ \Delta S=0$ ? (MP 6.19)

    In a real heat engine, with friction and losses, why is $ \Delta S$ still 0 if $ TdS=dQ+d\Phi$ ? (MP 6.20)

  4. Propulsive Power and Entropy Flux

    The final example in this section combines a number of ideas presented in this subject and in Unified in the development of a relation between entropy generation and power needed to propel a vehicle.

    Figure 6.15 shows an aerodynamic shape (airfoil) moving through the atmosphere at a constant velocity. A coordinate system fixed to the vehicle has been adopted so that we see the airfoil as fixed and the air far away from the airfoil moving at a velocity $ c_0$ . Streamlines of the flow have been sketched, as has the velocity distribution at station ``0'' far upstream and station ``d'' far downstream.

    The airfoil has a wake, which mixes with the surrounding air and grows in the downstream direction. The extent of the wake is also indicated. Because of the lower velocity in the wake the area between the stream surfaces is larger downstream than upstream.

    Figure 6.15: Airfoil with wake and control volume for analysis of propulsive power requirement
    Image fig3AirfoilWithWake

    We use a control volume description and take the control surface to be defined by the two stream surfaces and two planes at station 0 and station d. This is useful in simplifying the analysis because there is no flow across the stream surfaces. The area of the downstream plane control surface is broken into $ A_1$ , which is the area outside the wake and $ A_2$ , which is the area occupied by wake fluid, i.e., fluid that has suffered viscous losses. The control surface is also taken far enough away from the vehicle so that the static pressure can be considered uniform. For fluid which is not in the wake (no viscous forces), the momentum equation is

    $\displaystyle cdc =-dP/\rho.$

    Uniform static pressure therefore implies uniform velocity, so that on $ A_1$ the velocity is equal to the upstream value, $ c_0$ . The downstream velocity profile is actually continuous, as indicated. It is approximated in the analysis as a step change to make the algebra a bit simpler. (The conclusions apply to the more general velocity profile as well and we would just need to use integrals over the wake instead of the algebraic expressions below.)

    The equation expressing mass conservation for the control volume is

    $\displaystyle \rho_0A_0c_0 = \rho_0 A_1 c_0 + \rho_2 A_2 c_2.$ (6..12)

    The vertical face of the control surface is far downstream of the body. By this station, the wake fluid has had much time to mix and the velocity in the wake is close to the free stream value, $ c_0$ . We can thus write,

    $\displaystyle \textrm{wake velocity }= c_2 =(c_0-\Delta c);\qquad \Delta c/c_0\ll 1.$ (6..13)

    (We chose our control surface so the condition $ \Delta c/c_0 \ll 1$ was upheld.)

    The integral momentum equation (control volume form of the momentum equation) can be used to find the drag on the vehicle:

    $\displaystyle \rho_0 A_0 c_0^2 = -\textrm{Drag} +\rho_0 A_1 c_0^2 +\rho_2 A_2 c_2^2.$ (6..14)

    There is no pressure contribution in Eq. (6.14) because the static pressure on the control surface is uniform. Using the form given for the wake velocity and expanding the terms in the momentum equation we obtain

    $\displaystyle \rho_0 A_0 c_0^2 = -\textrm{Drag} +\rho_0 A_1 c_0^2 +\rho_2 A_2[c_0^2 - 2c_0\Delta c + (\Delta c)^2].$ (6..15)

    The last term in the right hand side of the momentum equation, $ \rho_2 A_2 (\Delta c)^2$ , is small by virtue of the choice of control surface and we can neglect it. Doing this and grouping the terms on the right hand side of Eq. (6.15) in a different manner, we have

    $\displaystyle c_0[\rho_0 A_0 c_0] = c_0[\rho_0 A_1 c_0 +\rho_2 A_2 (c_0-\Delta c)] +\{-\textrm{Drag}
-\rho_2 A_2 c_0 \Delta c\}.$

    The terms in the square brackets on both sides of this equation are the continuity equation multiplied by $ c_0$ . They thus sum to zero leaving the curly bracketed terms as

    $\displaystyle \textrm{Drag} = -\rho_2 A_2 c_0 \Delta c.$ (6..16)

    The wake mass flow is $ \rho_2 A_2 c_2 = \rho_2 A_2 (c_0 - \Delta
c)$ . All this flow has a velocity defect (compared to the free stream) of $ \Delta c$ , so that the defect in flux of momentum (the mass flow in the wake times the velocity defect) is, to first order in $ \Delta c$ ,

    $\displaystyle \textrm{Momentum defect in wake }= -\rho_2 A_2 c_0 \Delta c =
\textrm{Drag}.$

    The combined first and second law gives us a means of relating the entropy and velocity:

    $\displaystyle Tds = dh - dP/\rho.$

    The pressure is uniform ($ dP=0$ ) at the downstream station. There is no net shaft work or heat transfer to the wake so that the mass flux of stagnation enthalpy is constant. We can also approximate that the condition of constant stagnation enthalpy holds locally on all streamlines. Applying both of these to the combined first and second law yields

    $\displaystyle Tds = dh_t -cdc. $

    For the present situation, $ dh_t = 0$ ; $ cdc = c_0 \Delta c$ , so that

    $\displaystyle T_0\Delta s =-c_0 \Delta c.$ (6..17)

    In Equation (6.17) the upstream temperature is used because differences between wake quantities and upstream quantities are small at the downstream control station. The entropy can be related to the drag as

    $\displaystyle \textrm{Drag} = \rho_2 A_2 T_0 \Delta s.$ (6..18)

    The quantity $ \rho_2 A_2 c_0 \Delta s$ is the entropy flux (mass flux times the entropy increase per unit mass; in the general case we would express this by an integral over the locally varying wake velocity and density). The power needed to propel the vehicle is the product of $ \textrm{drag} \times \textrm{flight speed}$ , $ \textrm{Drag} \times c_0$ . From Eq. (6.18), this can be related to the entropy flux in the wake to yield a compact expression for the propulsive power needed in terms of the wake entropy flux:

    $\displaystyle \textrm{Propulsive power needed} =T_0 (\rho_2 A_2 c_0 \Delta s)= T_0 \times \textrm{Entropy flux in wake}.$ (6..19)

    This amount of work is dissipated per unit time in connection with sustaining the vehicle motion. Equation (6.19) is another demonstration of the relation between lost work and entropy generation, in this case manifested as power that needs to be supplied because of dissipation in the wake.





Muddy Points

Is it safe to say that entropy is the tendency for a system to go into disorder? (MP 6.21)

Douglas Quattrochi 2006-08-06