05.04.2016 Views

Modern Engineering Thermodynamics

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

16.6 Choked Flow 665<br />

T a<br />

a<br />

p a = constant<br />

T<br />

T b<br />

T a′<br />

b<br />

a ′<br />

p exit = p b = constant<br />

FIGURE 16.14<br />

Isentropic expansion of a vapor near its saturation at state a through a nozzle to a nonequilibrium supersaturated state at a’. This is<br />

followed by an irreversible nucleation and condensation process to equilibrium state b.<br />

s<br />

16.6 CHOKED FLOW<br />

The mass flow rate per unit area in an isentropic nozzle can be determined for an ideal gas from Eqs. (16.11),<br />

(16.12), and (16.13) and the ideal gas formula, ρ = p/RT, as<br />

_m /A = ρV = ðp/RTÞðMcÞ<br />

pffiffiffiffiffiffiffiffiffiffiffiffi<br />

= pMð kg c RTÞðT os /T os Þ 1/2 ðp os /p os Þ/RT<br />

= p os Mðkg c /RT os Þ 1/2 ðT os /TÞ 1/2 ðp/p os Þ<br />

(16.21)<br />

<br />

= p os M kg 1/2 <br />

c<br />

1 + k − 1 ðk+1Þ/2ð1−kÞ<br />

RT os 2 M2<br />

Figure 16.15 is a schematic of the variation in _m /A with the ratio of the back pressure p B to the upstream stagnation<br />

pressure p os .WhenM throat = M = 1:0, the nozzle is passing the maximum possible flow _m max :<br />

Since we always have M throat ≤ 1:0, clearly the maximum nozzle mass flow rate occurs when the throat velocity<br />

is sonic. Then _m = _m max ,M throat = M = 1:0, A throat = A , and Eq. (16.21) becomes<br />

1/2<br />

_m max /A kg<br />

<br />

= p c k + 1 ðk+1Þ/2ð1−kÞ<br />

os<br />

(16.22a)<br />

2<br />

RT os<br />

Note that the value of _m max /A depends only on the<br />

upstream isentropic stagnation properties and is completely<br />

independent of the downstream conditions. For air,<br />

k = 1.40 and, in <strong>Engineering</strong> English units, this equation<br />

reduces to<br />

<br />

ð _m max /A Þ air = 0:532 lbm pffiffiffi<br />

<br />

<br />

. R<br />

os<br />

pffiffiffiffiffiffi<br />

(16.22b)<br />

lbf .s T os<br />

where _m max must be in lbm/s, A must be in in 2 , p os<br />

must be in lbf/in 2 ,andT os must be in R. This is called<br />

Fliegner’s (or sometimes Zeuner’s) formula and was<br />

experimentally discovered in the 1870s. In metric<br />

SI units, Eq. (16.22a) reduces to<br />

<br />

ð _m max /A Þ air<br />

= 0:0404 kg pffiffiffi<br />

<br />

<br />

. K<br />

os<br />

pffiffiffiffiffiffi<br />

(16.22c)<br />

N.s T os<br />

where _m max must be in kg/s, A must be in m 2 , p os must<br />

be in N/m 2 ,andT os must be in K.<br />

m/A<br />

m max /A*<br />

Choked flow<br />

region<br />

M throat = 1<br />

0 p B /p os<br />

M throat < 1<br />

FIGURE 16.15<br />

The relative variation in dimensionless isentropic converging nozzle<br />

air mass flow rate with increasing back pressure to upstream<br />

stagnation pressure ratio. Since the Mach number in a converging<br />

nozzle cannot exceed 1.0, the maximum flow rate through the<br />

nozzle occurs when the Mach number at the throat is 1.0.<br />

1.0

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!