Skip to article frontmatterSkip to article content
Site not loading correctly?

This may be due to an incorrect BASE_URL configuration. See the MyST Documentation for reference.

6.4 Geometric Choking

Chalmers University of Technology
Department of Mechanical Engineering
Division of Fluid Dynamics

For steady-state nozzle flow, the massflow is obtained as

m˙=ρuA=const\dot{m}=\rho uA=const

Eqn. (6.32) can be evaluated at any location inside the nozzle and if evaluated at sonic conditions we get

m˙=ρuA\dot{m}=\rho^* u^* A^*

By definition u=au^*=a^* and thus

m˙=ρaA\dot{m}=\rho^* a^* A^*

ρ\rho^* and aa^* can be obtained using the ratios ρ/ρo\rho^*/\rho_o and a/aoa^*/a_o

ρ=(ρρo)ρo=poRTo(2γ+1)1/(γ1)\rho^*=\left(\frac{\rho^*}{\rho_o}\right)\rho_o=\frac{p_o}{RT_o}\left(\frac{2}{\gamma+1}\right)^{1/(\gamma-1)}
a=(aao)ao=ao(2γ+1)1/2=γRTo(2γ+1)1/2a^*=\left(\frac{a^*}{a_o}\right)a_o=a_o\left(\frac{2}{\gamma+1}\right)^{1/2}=\sqrt{\gamma RT_o}\left(\frac{2}{\gamma+1}\right)^{1/2}

Eqns. (6.36) and (6.35) in Eqn. (6.34) gives

m˙=poRTo(2γ+1)1/(γ1)γRTo(2γ+1)1/2A\dot{m}=\frac{p_o}{RT_o}\left(\frac{2}{\gamma+1}\right)^{1/(\gamma-1)} \sqrt{\gamma RT_o}\left(\frac{2}{\gamma+1}\right)^{1/2} A^*

which can be rewritten as

m˙=poAToγR(2γ+1)(γ+1)/(γ1)\dot{m}=\frac{p_o A^*}{\sqrt{T_o}}\sqrt{\frac{\gamma}{R}\left(\frac{2}{\gamma+1}\right)^{(\gamma+1)/(\gamma-1)}}

Eqn. (6.38) valid for:

It should be noted that the choked massflow can be calculated using Eqn. (6.38) even for cases with shocks downstream of the throat.