Converging-Diverging Nozzle Dynamics

The converging-diverging (CD) nozzle, conceived by Swedish engineer Carl Gustaf Patrik de Laval in 1888 for steam impulse turbines, stands as one of the quintessential aerodynamic geometries in aerospace propulsion. Without a convergent-divergent area variation, accelerating a continuous gaseous flow from stagnation conditions to supersonic velocities ($M > 1$) is physically impossible in steady 1D flow.

In modern rocketry, ramjets, and supersonic jet exhausts, de Laval nozzles convert high-temperature, high-pressure thermal stagnation enthalpy ($h_0$) stored inside a combustion reservoir into high-velocity directed kinetic energy ($\frac{1}{2}V^2$), maximizing thrust ($F = \dot{m}V_e + (p_e - p_a)A_e$) and specific impulse ($I_{sp}$).

A de Laval nozzle consists of three distinct thermodynamic regions as illustrated in the vector schematic below:

• Stagnation Reservoir / Chamber ($p_0, T_0$): A large volume where flow speed is negligible ($V \approx 0$). Pressure and temperature equal total stagnation properties.
• Convergent Section ($dA/dx < 0$): Channels gas into a narrowing duct. Under subsonic conditions ($M < 1$), decreasing cross-sectional area forces flow acceleration ($dV > 0$).
• Throat ($A = A^* = A_{min}$): The location of minimum area. According to the Area-Velocity relation, Mach 1 ($M = 1$) can only occur precisely at this geometric station.
• Divergent Section ($dA/dx > 0$): An expanding area duct. In compressible supersonic flow ($M > 1$), expanding area accelerates the fluid further, increasing velocity at the expense of temperature and pressure.

Consider an experiment where a high-pressure reservoir ($p_0$) exhausts through a CD nozzle into an ambient receiver with variable back pressure ($p_b$). As $p_b$ is progressively reduced from $p_0$:

1. The mass flow rate ($\dot{m}$) initially increases proportionally to the pressure differential.
2. At a critical back pressure threshold ($p_b = p_{choke}$), acoustic velocity ($V = a^*$) is attained at the minimum area station ($M_{throat} = 1$).
3. Choked Condition: Any further reduction in back pressure ($p_b \to 0$) fails to transmit pressure disturbance waves upstream past the sonic throat, because acoustic disturbances propagate at speed $a$, whereas the fluid at the throat moves at speed $a$ in the opposite direction ($u - a = 0$).

Hence, once choked, the mass flow rate through the nozzle remains strictly constant and frozen, irrespective of whether the nozzle exhausts into atmospheric sea-level air or absolute hard vacuum in space.

Depending on the ratio of ambient back pressure to stagnation pressure ($p_b/p_0$), the nozzle traverses seven distinct aerodynamic operational states:

Plotting the normalized pressure ratio ($p(x)/p_0$) along the normalized nozzle length ($x/L$) provides a graphical representation of compressible internal aerodynamics:

The computational core of AeroFlow is driven by closed-form analytical relationships coupled with bisection numerical inverses:

The AeroFlow Simulator translates the analytical formulations above into a visual interactive laboratory environment. Developed at the Department of Aerospace Engineering, MIT Campus, it features: