# Rayleigh flow

Rayleigh flow

Rayleigh flow refers to diabatic flow through a constant area duct where the effect of heat addition or rejection 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. Rayleigh flow is named after John Strutt, 3rd Baron Rayleigh.

Theory

The Rayleigh flow model begins with a differential equation that relates the change in Mach number with the change in stagnation temperature, T0. The differential equation is shown below.

$frac\left\{dM^2\right\}\left\{M^2\right\} = frac\left\{1 + gamma M^2\right\}\left\{1 - M^2\right\}left\left(1 + frac\left\{gamma - 1\right\}\left\{2\right\}M^2 ight\right)frac\left\{dT_0\right\}\left\{T_0\right\}$

Solving the differential equation leads to the relation shown below, where T0* is the stagnation temperature at the throat location of the duct which is required for thermally choking the flow.

$frac\left\{T_0\right\}\left\{T_0^*\right\} = frac\left\{2left\left(gamma + 1 ight\right)M^2\right\}\left\{left\left(1 + gamma M^2 ight\right)^2\right\}left\left(1 + frac\left\{gamma - 1\right\}\left\{2\right\}M^2 ight\right)$

These values are significant in the design of combustion systems. For example, if a turbojet combustion chamber has a maximum temperature of T0* = 2000 K, T0 and M at the entrance to the combustion chamber must be selected so thermal choking does not occur, which will limit the mass flow rate of air into the engine and decrease thrust.

For the Rayleigh flow model, the dimensionless change in entropy relation is shown below.

$Delta S = frac\left\{Delta s\right\}\left\{c_p\right\} = lnleft \left[M^2left\left(frac\left\{gamma + 1\right\}\left\{1 + gamma M^2\right\} ight\right)^frac\left\{gamma + 1\right\}\left\{gamma\right\} ight\right]$

The above equation can be used to plot the Rayleigh line on a Mach number versus ΔS graph, but the dimensionless enthalpy, H, versus ΔS diagram is more often used. The dimensionless enthalpy equation is shown below with an equation relating the static temperature with its value at the choke location for a calorically perfect gas where the heat capacity at constant pressure, cp, remains constant.

$H = frac\left\{h\right\}\left\{h^*\right\} = frac\left\{c_pT\right\}\left\{c_pT^*\right\} = frac\left\{T\right\}\left\{T^*\right\}$

$frac\left\{T\right\}\left\{T^*\right\} = left \left[frac\left\{left\left(gamma + 1 ight\right)M\right\}\left\{1 + gamma M^2\right\} ight\right] ^2$

The above equation can be manipulated to solve for M as a function of H. However, due to the form of the T/T* equation, a complicated multi-root relation is formed for M = M(T/T*). Instead, M can be chosen as an independent variable where ΔS and H can be matched up in a chart as shown in Figure 1. Figure 1 shows that heating will increase an upstream, subsonic Mach number until M = 1 and the flow chokes. Conversely, adding heat to a duct with an upstream, supersonic Mach number will cause the Mach number to decrease until the flow chokes. Cooling produces the opposite result for each of those two cases. The Rayleigh flow model reaches maximum entropy at M = 1. For subsonic flow, the maximum value of H occurs at M = .845. This indicates that cooling, instead of heating, causes the Mach number to move from .845 to 1. This is not necessarily correct as the stagnation temperature always increases to move the flow from a subsonic Mach number to M = 1, but from M = .845 to M = 1 the flow accelerates faster than heat is added to it. Therefore, this is a situation where heat is added but T/T* decreases in that region.

Additional Rayleigh Flow Relations

The area and mass flow rate are held constant for Rayleigh flow. Unlike Fanno flow, the Fanning friction factor, "f", remains constant. These relations are shown below with the * symbol representing the throat location where choking can occur.

$A = A^* = constant$

$dot\left\{m\right\} = dot\left\{m\right\}^* = constant$

Differential equations can also be developed and solved to describe Rayleigh flow property ratios with respect to the values at the choking location. The ratios for the pressure, density, static temperature, velocity and stagnation pressure are shown below, respectively. They are represented graphically along with the stagnation temperature ratio equation from the previous section. A stagnation property contains a '0' subscript.

$frac\left\{p\right\}\left\{p^*\right\} = frac\left\{gamma + 1\right\}\left\{1 + gamma M^2\right\}$

$frac\left\{ ho\right\}\left\{ ho^*\right\} = frac\left\{1 + gamma M^2\right\}\left\{left\left(gamma + 1 ight\right)M^2\right\}$

$frac\left\{T\right\}\left\{T^*\right\} = frac\left\{left\left(gamma + 1 ight\right)^2M^2\right\}\left\{left\left(1 + gamma M^2 ight\right)^2\right\}$

$frac\left\{V\right\}\left\{V^*\right\} = frac\left\{left\left(gamma + 1 ight\right)M^2\right\}\left\{1 + gamma M^2\right\}$

$frac\left\{p_0\right\}\left\{p_0^*\right\} = frac\left\{gamma + 1\right\}\left\{1 + gamma M^2\right\}left \left[left\left(frac\left\{2\right\}\left\{gamma + 1\right\} ight\right)left\left(1 + frac\left\{gamma - 1\right\}\left\{2\right\}M^2 ight\right) ight\right] ^frac\left\{gamma\right\}\left\{gamma - 1\right\}$

Applications

The Rayleigh flow model has many analytical uses, most notably involving aircraft engines. For instance, the combustion chambers inside turbojet engines usually have a constant area and the fuel mass addition is negligible. These properties make the Rayleigh flow model applicable for heat addition to the flow through combustion, assuming the heat addition does not result in dissociation of the air-fuel mixture. Producing a shock wave inside the combustion chamber of an engine due to thermal choking is very undesirable due to the decrease in mass flow rate and thrust. Therefore, the Rayleigh flow model is critical for an initial design of the duct geometry and combustion temperature for an engine.

The Rayleigh flow model is also used extensively with the Fanno flow model. These two models intersect at points on the enthalpy-entropy and Mach number-entropy diagrams, which is meaningful for many applications. However, the entropy values for each model are not equal at the sonic state. The change in entropy is 0 at M = 1 for each model, but the previous statement means the change in entropy from the same arbitrary point to the sonic point is different for the Fanno and Rayleigh flow models. If initial values of si and Mi are defined, a new equation for dimensionless entropy versus Mach number can be defined for each model. These equations are shown below for Fanno and Rayleigh flow, respectively.

$Delta S_F = frac\left\{s - s_i\right\}\left\{c_p\right\} = lnleft \left[left\left(frac\left\{M\right\}\left\{M_i\right\} ight\right)^frac\left\{gamma - 1\right\}\left\{gamma\right\}left\left(frac\left\{1 + frac\left\{gamma - 1\right\}\left\{2\right\}M_i^2\right\}\left\{1 + frac\left\{gamma - 1\right\}\left\{2\right\}M^2\right\} ight\right)^frac\left\{gamma + 1\right\}\left\{2gamma\right\} ight\right]$

$Delta S_R = frac\left\{s - s_i\right\}\left\{c_p\right\} = lnleft \left[left\left(frac\left\{M\right\}\left\{M_i\right\} ight\right)^2left\left(frac\left\{1 + gamma M_i^2\right\}\left\{1 + gamma M^2\right\} ight\right)^frac\left\{gamma + 1\right\}\left\{gamma\right\} ight\right]$

Figure 5 shows the Rayleigh and Fanno lines intersecting with each other for initial conditions of si = 0 and Mi = 3. The intersection points are calculated by equating the new dimensionless entropy equations with each other, resulting in the relation below.

$left\left(1 + frac\left\{gamma - 1\right\}\left\{2\right\}M_i^2 ight\right)left \left[frac\left\{M_i^2\right\}\left\{left\left(1 + gamma M_i^2 ight\right)^2\right\} ight\right] = left\left(1 + frac\left\{gamma - 1\right\}\left\{2\right\}M^2 ight\right)left \left[frac\left\{M^2\right\}\left\{left\left(1 + gamma M^2 ight\right)^2\right\} ight\right]$

Interestingly, the intersection points occur at the given initial Mach number and its post-normal shock value. For Figure 5, these values are M = 3 and .4752, which can be found the normal shock tables listed in most compressible flow textbooks. A given flow with a constant duct area can switch between the Rayleigh and Fanno models at these points.

*Fanno flow
*Isentropic process
*Gas dynamics
*Compressible flow
*Choked flow
*Enthalpy
*Entropy

References

* cite book
last = Zucker | first = Robert D. | coauthors = Biblarz O.
title = Fundamentals of Gas Dynamics | origyear = 2002
publisher = John Wiley & Sons
id = ISBN 0471059676

* cite book
last = Shapiro | first = Ascher H.
title = The Dynamics and Thermodynamics of Compressible Fluid Flow, Volume 1 | origyear = 1953
publisher = Ronald Press
id = ISBN 978-0-471-06691-0

* cite book
last = Hodge | first = B. K. | coauthors = Koenig K.
title = Compressible Fluid Dynamics with Personal Computer Applications | origyear = 1995
publisher = Prentice Hall
id = ISBN 013308552X

* cite book
last = Emanuel | first = G.
title = Gasdynamics: Theory and Applications | origyear = 1986
publisher = AIAA
id = ISBN 0930403126

* [http://meweb.ecn.purdue.edu/~meapplet/java/comp_calculator/rayleigh_calculator.html Purdue University Rayleigh flow calculator]
* [http://www.engr.uky.edu/~jdjacob/me530/rayleigh.html University of Kentucky Rayleigh flow Webcalculator]

Wikimedia Foundation. 2010.

Поможем решить контрольную работу

### Look at other dictionaries:

• Rayleigh–Bénard convection — Bénard cells. Rayleigh–Bénard convection is a type of natural convection, occurring in a plane of fluid heated from below, in which the fluid develops a regular pattern of convection cells known as Bénard cells. Rayleigh–Bénard convection is one… …   Wikipedia

• Rayleigh–Taylor instability — The Rayleigh–Taylor instability, or RT instability (after Lord Rayleigh and G. I. Taylor), is an instability of an interface between two fluids of different densities, which occurs when the lighter fluid is pushing the heavier fluid. [citation… …   Wikipedia

• Rayleigh number — In fluid mechanics, the Rayleigh number for a fluid is a dimensionless number associated with buoyancy driven flow (also known as free convection or natural convection). When the Rayleigh number is below the critical value for that fluid, heat… …   Wikipedia

• Rayleigh-Benard-Experiment — Schematische Darstellung der Bénard Zellen in Seitenansicht Beim Bénard Experiment (Bénard Effekt, Bénard System, Bénard Zellen, Rayleigh Benard Experiment) wird eine dünne, homogene Flüssigkeitsschicht die sich in einem festen Behältnis und in… …   Deutsch Wikipedia

• Rayleigh Medal — The Rayleigh Medal is a prize awarded annually by the Institute of Acoustics for outstanding contributions to acoustics . The prize is named after John William Strutt or (Lord Rayleigh).Recipients of this award include* 1970 David E. Weston... *… …   Wikipedia

• Fanno flow — refers to 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,… …   Wikipedia

• Compressible flow — is the area of fluid mechanics that deals with fluids in which the fluid density varies significantly in response to a change in pressure. Compressibility effects are typically considered significant if the Mach number (the ratio of the flow… …   Wikipedia

• D'Alembert's paradox — In fluid dynamics, d Alembert s paradox (or the hydrodynamic paradox) is a contradiction reached in 1752 by French mathematician Jean le Rond d Alembert.[1] D Alembert proved that – for incompressible and inviscid potential flow – the drag force… …   Wikipedia

• analysis — /euh nal euh sis/, n., pl. analyses / seez /. 1. the separating of any material or abstract entity into its constituent elements (opposed to synthesis). 2. this process as a method of studying the nature of something or of determining its… …   Universalium

• Multiangle light scattering — (MALS) is a technique for determining, independently, the absolute molar mass and the average size of particles in solution, by detecting how they scatter light. Collimated light from a laser source is most often used, in which case the technique …   Wikipedia