Theory
Isentropic flow relations:
$$\frac{T}{T_0}=\left(1+\frac{\gamma-1}{2}M^2\right)^{-1}$$ $$\frac{A}{A^*}=\frac{1}{M}\left[\frac{2}{\gamma+1}\left(1+\frac{\gamma-1}{2}M^2\right)\right]^{\frac{\gamma+1}{2(\gamma-1)}}$$Set exit Mach number, throat area, and stagnation conditions to instantly compute isentropic flow distributions of Mach, pressure, and temperature along the nozzle axis. Optional normal shock visualization included.
Isentropic flow relations:
$$\frac{T}{T_0}=\left(1+\frac{\gamma-1}{2}M^2\right)^{-1}$$ $$\frac{A}{A^*}=\frac{1}{M}\left[\frac{2}{\gamma+1}\left(1+\frac{\gamma-1}{2}M^2\right)\right]^{\frac{\gamma+1}{2(\gamma-1)}}$$The core of the design relies on isentropic (constant entropy) flow relations. The first key equation relates the local static temperature (T) to the stagnation temperature (T₀), which is the temperature if the flow were brought to rest adiabatically.
$$\frac{T}{T_0}=\left(1+\frac{\gamma-1}{2}M^2\right)^{-1}$$Here, M is the local Mach number, and γ is the specific heat ratio (cp/cv). As M increases, T drops significantly. This is why nozzles get very cold internally.
The most critical equation for shape design is the area-Mach number relation. It defines the nozzle's cross-sectional area (A) at any point relative to the throat area (A*), where M=1.
$$\frac{A}{A^*}=\frac{1}{M}\left[\frac{2}{\gamma+1}\left(1+\frac{\gamma-1}{2}M^2\right)\right]^{\frac{\gamma+1}{2(\gamma-1)}}$$For a given exit Mach number (M_exit) and γ, this equation gives the expansion area ratio. This is the value you are directly calculating with the simulator's input parameters. A/A* is minimum (1) at the throat and grows for both subsonic and supersonic Mach numbers.
Rocket Engine Nozzles: This is the most iconic application. The bell-shaped nozzle on a rocket is a De Laval nozzle. It maximizes thrust by efficiently expanding hot, high-pressure combustion gases to supersonic speeds. The simulator's stagnation conditions (P₀, T₀) directly model the conditions in a rocket's combustion chamber.
Supersonic Wind Tunnels: To test aircraft and missile models at supersonic speeds, engineers use wind tunnels with a De Laval nozzle in the test section. By carefully designing the area ratio, they can set a precise Mach number for testing, just like you can set the exit Mach number in this tool.
Steam Turbines (Historical & Modern): Gustaf de Laval originally invented this nozzle for steam turbines. It allowed steam to reach supersonic speeds, dramatically increasing the efficiency and power output of early turbines. The principle is still used in some high-power applications today.
Jet Engine Afterburners & Nozzles: Military fighter jets with afterburners often use variable-geometry convergent-divergent nozzles. During afterburner use, the nozzle opens up to act as a De Laval nozzle, allowing the higher-temperature exhaust to expand to supersonic speeds for extra thrust.
When you start using this tool, there are a few common pitfalls to watch out for. First, make sure you understand the assumption that "the inlet condition comes from an infinitely large tank." The simulator calculates starting from a state where the flow velocity at the inlet is nearly zero (Mach number ≈ 0). However, an actual engine's combustion chamber is finite in size and has its own flow. If the inlet Mach number is 0.1 or higher, be aware that the simple area ratio formula used by this tool alone won't determine the accurate shape.
Next, be careful of errors in setting the "specific heat ratio γ." A common mistake is designing a rocket nozzle for air (γ=1.4). Rocket combustion gases contain a lot of water vapor and carbon dioxide, so γ can be around 1.2. For example, simply changing γ from 1.4 to 1.25 increases the required area ratio to achieve the same exit Mach number of 5.0 by about 40. Always double-check the gas you're using.
Finally, there's the "gap between ideal and reality." This calculation is based on "one-dimensional isentropic flow," which completely ignores wall friction, heat losses, and two-dimensional flow effects. In actual nozzles, especially small ones, the boundary layer can reduce the effective flow area, preventing the design performance from being achieved. Remember, even if you achieve an "optimal design" in simulation, the real design work begins from there.
The principle of the de Laval nozzle is fundamental to all fields dealing with supersonic flow. First, consider "turbomachinery." The stator vanes in jet engine turbines and compressors function as nozzles and diffusers. While the flow passages in turbine cascades have complex shapes, unlike the smooth contraction and expansion handled by this tool, the fundamental principle of converting energy by accelerating and decelerating the flow is the same.
Next, it's also applied in "HVAC (Heating, Ventilation, and Air Conditioning) duct design." Especially for high-velocity air supply or cleanroom air outlets, nozzle design principles are utilized to prevent flow separation and achieve uniform flow, even in the subsonic regime. Conversely, for exhaust ducts, understanding them as "diffusers," where the flow slows down in an expanding section to recover pressure (reducing pressure loss), becomes crucial.
Furthermore, there are interesting applications in the field of "MEMS (Micro-Electro-Mechanical Systems)." In micro- and nano-scale gas flow channels, gas behavior deviates from continuum fluid mechanics (rarefied gas effects). However, supersonic micro-nozzles are researched for micro-thrusters and analytical devices, and their basic design still starts with this area ratio equation. The fundamental principles remain universal, even when the scale changes.
Once you're comfortable with this tool's calculations, your next step is to grasp the mathematical background of "why it works that way." The key lies in the "Saint-Venant and Wantzel equations," derived by combining the conservation laws of mass, momentum, and energy with the equation of state. In particular, how flow velocity changes with cross-sectional area $A$ can be understood from the following differential relation:
$$ \frac{dA}{A} = (M^2 -1)\frac{dV}{V} $$
This single equation tells the whole story. If $M<1$ (subsonic), $dA$ and $dV$ have opposite signs (acceleration in a converging duct). If $M>1$ (supersonic), they have the same sign (acceleration in a diverging duct). When $M=1$, $dA$ must be zero (the throat). From this one equation, the necessary shape of the de Laval nozzle becomes clear.
To get closer to practical application, your next steps should be learning about "behavior under off-design conditions" and "two-dimensional and three-dimensional effects." Learn how to quantitatively predict why the "shock waves" you experimented with in the tool form and how they move, based on "pressure ratios." Furthermore, real nozzle flow is axisymmetric and two-dimensional, where the wall "boundary layer" interacts with the central "core flow." Using CFD (Computational Fluid Dynamics) to visualize these phenomena will let you experience the difference from the ideal case and should significantly deepen your understanding.