This function computes the divergent section of an axisymmetric minimum length nozzle for a perfect gas (inviscid, compressible, irrotational flow). The kernel and transition regions are computed using the axisymmetric method of the characteristics. A normal sonic line is assumed at the throat.
The function requires as input:
- the exit Mach number (>1, supersonic);
- the heat capacity ratio (gamma=cp/cv);
- the number of characteristics (integer >0) in the linear kernel region;
- the number of characteristics (integer >0) in the compressed region downstream of the throat (x=0). When equal to one, the compression region is not resolved;
- the exponent of the power-law that defines the characteristics in the compressed region;
- maximum tolerance (residuals)
- aspect ratio of the mesh in the transition region (usually close to 1)
The function outputs a non-dimensional plot. Lengths are normalized to the throat radius. The code stores the coordinates of the contour points in a .txt file.
The code is based on the following references:
 J.D. Anderson. Modern Compressible Flow: With Historical Perspective. Aero-
nautical and Aerospace Engineering Series. McGraw-Hill Education, 2003 (Library of Congress CN: 2002067852)
 B. M. Argrow and G. Emanuel. Comparison of Minimum Length Nozzles. Journal of Fluids Engineering, 110(3):283-288, 09 1988.
 Kuno Foelsch. The analytical design of an axially symmetric Laval nozzle for a
parallel and uniform jet. Journal of the Aeronautical Sciences, 16(3):161-166,
 Ying-Nien Yu. A summary of design techniques for axisymmetric hypersonic
wind tunnels. Technical report, NATO - Science and Technology Organization,
Ludovico Fossà (2021). Axisymmetric Minimum Length Nozzle Design Tool (https://www.mathworks.com/matlabcentral/fileexchange/85253-axisymmetric-minimum-length-nozzle-design-tool), MATLAB Central File Exchange. Retrieved .
Find the treasures in MATLAB Central and discover how the community can help you!Start Hunting!