AeroSpike 3.1
| MAIN PAGE | RESUME |


 AeroSpike Nozzle Method of Characteristics Program
Copyright © 1999-2022 John Cipolla/AeroRocket


Main start-up screen for the Aerospike Nozzle and Minimum Length Nozzle (MLN) analyses


Aerospike Rocket Motor Designed Using AeroSpike 3.1 NEW


PART 1: AEROSPIKE NOZZLE DESIGN
ANNULAR (3-D) & LINEAR (2-D) CONTOURS
BACK TO TOP

Summary of Features
1. Determine the shape of an annular or linear aerospike nozzle given the thruster exit area ratio (Aei/At), projected area expansion ratio (Ae/At), pressure ratio (Pc/Pa), thruster internal radius (Ra), radius to lip of cowl (Re), total nozzle length from origin (Lnozzle), chamber temperature (Tc), chamber pressure (Pc), ratio of specific heats (
g) and gas constant.
2. Click the UpDown command button to move a locator to one of seven points in the flow field.
3. All important flow properties are displayed in real time as the locator moves from point to point in the flow field described by the characteristic mesh of the aerospike.
4. Generate color contour plots of Mach number (Mn), Pressure (P/Pc), Temperature (T/Tc) and density (R/Rc) with a single click.
5. Plot Mn, P/Pc, T/Tc, R/Rc, CF, CFvac, ISP, ISPvac as a function of aerospike nozzle axial location at a particular Pc/Pa.
6. Plot CF, Thrust, CFvac, ISP, and ISPvac verses pressure ratio (Pc/Pa) on a semi-log scale.
7. Units include, MKS (meter-newton-sec), CGS (centimeter-dyne-second), FPS (foot-pound-second) and IPS (inch-pound-second).
8. Graphically display the outer flow boundary for under expanded flow, over expanded flow and the angle of the outer boundary flow.
9. Graphically display the initial shock wave formed at the lip of the cowl for over expanded flow (Pa/Pc > Pe/Pc) and the shock angle from the lip.
10. Define gas properties for inert gases, liquid propellant gases and solid fuel propellant gases or insert your own values.
11. Define the analysis for annular (3-D) or linear (2-D) aerospike nozzles.
12. Define the angle the sonic section of the thruster makes with the axis of the aerospike nozzle.
13. Added a hybrid rocket motor propellant having the following fuel and oxidizer to the list of combustion gases: 85% Nitrous Oxide, 15% HTPB.
14. Added the ability to save F(x) verses PR (Pressure Ratio) and F(x) verses x to a CSV file for use with Notepad or Excel.
15.
In the Aerospike Nozzle Data section added a display of Truncation as percent of total aerospike length.
16.
In the Aerospike Nozzle Data section added a display of Distance from throat (origin) to end of thruster.
17.
In the Aerospike Nozzle Data section added a display of Distance from end of thruster to end of ramp.
18.
NEW! Added the ability to include base thrust for truncated aerospike nozzles.

Propellant Gases Available

Inert Gases
Dry Air Hydrogen Helium Water Vapor Argon Carbon Dioxide
Carbon Monoxide Nitrogen Oxygen Nitrogen Monoxide Nitrous Oxide Chlorine
Methane          
 
Liquid Fuel Propellant Gases
Oxygen, 75% Ethyl Alcohol(1.43) Oxygen, Hydrazine(.09) Oxygen, Hydrogen(4.02)
Oxygen, RP-1(2.56) Oxygen, UDMH(1.65) Fluorine, Hydrazine(2.3)
Fluorine, Hydrogen(7.60) Nitrogen Tetroxide, Hydrazine(1.34) Nitrogen Tetroxide, 50% UDMH, 50% Hydrazine(2.0)
Nitric Acid, RP-1(4.8) Nitric Acid, 50% UDMH, 50% Hydrazine(2.20) Liquid Oxygen, Liquid Methane (2.70, 2.80, 2.90)

Solid Fuel Propellant Gases
Ammonium Nitrate, 11% Binder, 4-20% Mg Ammonium Perchlorate, 18% Binder, 4-20% Al Ammonium Perchlorate, 12% Binder, 4-20% Al

Hybrid Rocket Motor Propellant Gases
85% Nitrous Oxide, 15% HTPB    

User-Defined Gases
Specify custom or user-defined gases by inserting Ratio of specific heats for exhaust (g) and Gas constant of exhaust (Rgas) in the Aerospike Nozzle Data section.

General Discussion
Aerospike NozzleAeroSpike performs an expansion-wave analysis from the throat of the thruster nozzle, where Mn = 1.0, to the thruster nozzle internal-exit as a series of simple wave expansions. Then, for the external ramp AeroSpike performs a series of Prandtl-Meyer expansions from the lip of the cowl, where R=Re, to the entire length of the external ramp of the aerospike nozzle. The ideal contour or shape of the external ramp of the aerospike nozzle is determined using isentropic supersonic flow theory. Then, depending on whether the flow is underexpanded or if the flow is overexpanded AeroSpike performs either a Prandtl-Meyer expansion analysis or an oblique shock wave analysis to determine the angle of the outer flow boundary from the lip of the cowl. As a by product of the oblique shock wave analysis AeroSpike determines the shock wave angle for overexpanded flow and plots both the outer boundary contour and the initial shock wave from the lip of the cowl. If the Check to Include base thrust check box is not checked then base pressure is assumed equal to atmospheric pressure (Pb = Patm) which means base thrust is zero and only the centerbody and thrusters contribute to total aerospike thrust. However, if the Check to Include base thrust check box is checked then base pressure and atmospheric pressure are non-equal resulting in the following aerospike nozzle total thrust equations, Ftotal = Fcenterbody + Fbase + Fthruster. and CF = Ftotal  / (At * Pc).

Procedure
From the menu on the top of the main start-up screen, select units (MNS, CGS, FPS or IPS) from the Units menu and then the propellant gas from the Gases menu. A number of inert gases, liquid fuel propellants and solid fuel propellants are available. The value for the ratio of specific heats (
g) are determined from the Units and Gases menus and are passed on to the Aerospike Nozzle program after clicking the Aerospike Nozzle command button on the main start-up screen. The ratio of specific heats, gas constant (Rgas), chamber pressure (Pc), and pressure ratio (PR) are required for the Aerospike Nozzle analysis. Additionally, the pressure ratio (PR) represents the maximum value for Pc/Pa that AeroSpike will use to plot CF, Thrust, CFvac , ISP and ISPvac as a function of PR. The chamber pressure is computed based on the atmospheric pressure (Pa) and the pressure ratio (Pc/Pa). These values are automatically passed to the Aerospike Nozzle analysis when the command button is clicked. However, the user can over-ride any input value by inserting new data directly into each data entry box on the Aerospike Nozzle Design screen. Each time the user changes any data entry the results are automatically updated and displayed. The user only needs to click the Plot button to see a new contour plot of the results or the UpDown button to see flow results at any of the characteristic mesh points.

Toolbar Operations

1. Click [X] to switch between the main data entry area (Figure-2) and the secondary data entry area (Figure-3). The main data entry area is displayed by default. Specify either annular aerospike geometry or linear aerospike geometry by clicking one of two option buttons in the secondary data entry area. In addition, the thruster sonic-section angle (60 degrees to 120 degrees) is located in the secondary data entry area. The thruster sonic-section angle is measured from the axis of the aerospike nozzle to the section that defines the throat of the thruster (where Mach number = 1). Default = 90 degrees. Finally, check the Check to Include base thrust check box to include base thrust for the determination of total thrust and thrust coefficient (CF) for truncated aerospike nozzles.
2. Send all flow properties (X, Y, Mn etc.) at each characteristic mesh point to the printer.
3. Send an image of the screen to the printer.
4. Save all flow properties (X,Y, Mn etc) at each characteristic mesh point to a data file.
5. Read the nozzle description file from a previous session.
6. Save the nozzle description file from a previous session.
7. Refresh the displayed analysis to the default analysis seen during start-up.
8. Return to the main start-up screen.

Input Variable Definitions
1. Thruster exit area ratio (Aei/At): Ratio of thruster internal exit area (Aei) to thruster throat area (At). Equation 1 is inverted to find Pc/Pei from Aei/At.
2. Thruster pressure ratio (Pc/Pei): Ratio of chamber pressure to thruster exit pressure. Found by interation of Equation 1 and displayed in lower data region.
3. Aerospike expansion ratio (Ae/At): The projected area of the aerospike nozzle (Ae =
p * Re^2) divided by the total thruster throat area.
4. Ratio of specific heats (
g): Selected from a pull-down menu or user-defined.
5. Gas constant of exhaust (Rgas): Selected from a pull-down menu or user-defined.
6. Aerospike pressure ratio (Pc/Pa): Ratio of the chamber pressure (Pc) to the atmospheric pressure (Pa).
7. Thruster internal circular radius (Ra): Radius of the internal portion of the thruster duct from point 1 (throat) to point 2 (Ra).
8. Radius to lip of cowl (Re): Radius that defines the projected area of the aerospike nozzle (Re).
9. Aerospike length from origin (Lnozzle): Total length of the aerospike nozzle from the origin (throat) of the thruster to the end of the ramp.
10. Chamber temperature (Tc): Chamber temperature in either degrees Rankine or degrees Kelvin depending on the units selected.
11. Chamber pressure (Pc): Chamber pressure whose units depend on the units selected.
12. Width of ramp for linear aerospike nozzles (Lramp).
.
Equation 1: Thruster Cross-Sectional Area and Pressure Ratio Relationship.


Figure 1. Aerospike Nozzle Displaying Basic Geometry and the External Expansion Fan.


AeroSpike Validation-1
EXTERNAL FLOW FIELD FOR AN ANNULAR AEROSPIKE NOZZLE


Figure 2. Aerospike Nozzle - Optimum Expansion (PR = 71.5) and 20% plug nozzle configuration


Figure 3. Aerospike Nozzle - Secondary Input Data Entry Area for Annular/Linear nozzle and thruster angle inputs.

Truncated Aerospike Nozzle Base Pressure, Total Thrust and Thrust Coefficient
To determine truncated aerospike nozzle base pressure (Pb) for computing total thrust and pressure coefficient (CF) simply check the Check to Include base thrust check box in the Annular or Linear Ramp Selection and Thruster Angle data entry area. When checked this check box activates base pressure computation for truncated aerospike nozzles where base pressure is included for determination of total aerospike thrust and thrust coefficient. Truncated aerospike nozzle base thrust is determined using two curve-fit relationships for computing base thrust. The first relationship is atmospheric pressure (Patm) verses pressure ratio (PR=Pc/Pa) and the second relationship is base pressure verses percent truncation. Where, X% truncation refers to an aerospike nozzle where (100-X)% of the expansion ramp has been removed leaving a blunt base region. The plots in Figure 4 and Figure 5 illustrate the relationship between Patm verses PR and Pb verses Percent Truncation for the computation of base thrust, total thrust and thrust coefficient. The curve-fit for base pressure verses percent truncation displayed in Figure 5 was developed using several Computational Fluid Dynamics (CFD) analyses using aerospike nozzles having 20%, 30%, 40% and 50% truncation.

Aerospike nozzle total thrust is computed using the following equations knowing PR and percent truncation.
Fbase= (Pb - Patm) * Abase, Ftotal = Fcenterbody + Fbase + Fthruster and CF = Ftotal  / (At * Pc) where Patm = fn(PR) and Pb = fn(% Truncation)


Figure 4. Atmospheric pressure (Pa) verses pressure ratio (PR=Pc/Patm)

 
Figure 5. Base pressure (Pb) verses percent truncation curve-fit using CFD analyses of aerospike nozzles having 20%, 30%, 40% and 50% truncation.
 
CF vs. PR Validation   Flow Field Validation


Figure 6. CF verses Pressure Ratio (Pc/Pa)  - Semi-Log plot, Maximum PR = 1000. AeroSpike CF verses PR compared to 20% (80% truncated) plug nozzle  Base Flow CFD analysis includes base pressure coefficient
. Reference: AIAA 2001-1051, T. Ito, K. Fujii, Flow Field Analysis of the Base Region of Axisymmetric Aerospike Nozzles.



"I used (AeroSpike) to design several types of nozzle(s) and found your software is really useful".
Takashi Ito, JAXA

Reference: "Aerospike Nozzle Flow Fields", AIAA 2001-1051, by Takashi Ito. Contour plots used with permission from reference.

 

Figure 7. AeroSpike program nozzle flow field results compared to Base Flow CFD analyses for PR=9, PR=71 and PR=500.
 

Aerospike Rocket Motor Design Example    

AeroCFD 5.2 Analysis by John Cipolla

Annular aerospike designed to produce 225,000 pounds of thrust.

 

AeroCFD analysis of an aerospike rocket. M = 1.5.

The annular aerospike rocket motor pictured above (left) is a design developed and rendered by Richard Caldwell (Rocket Nut) intended to produce 225,000 pounds of thrust at sea level. This design concept was developed using AeroSpike 2.6 software to specify the internal thruster and external ramp geometries for efficient operation from sea level to orbital altitude. Rocket Nut's design uses an untruncated Annular Aerospike ramp as illustrated in Figure-3. Please click here for more information.

AeroSpike Validation-2
X-33 XRS-2200 AEROSPIKE ROCKET ENGINE ANALYSIS EXAMPLE
X-33 spaceplaneIn this example the Linear Ramp or 2-D option is used to determine sea level and vacuum thrust for the X-33 XRS-2200 aerospike rocket engine. Click the Hide or Show Aerospike Nozzle Data (X in the toolbar) to select the Linear Aerospike option and then insert 90 for the 2-D Ramp width in the space provided. The data from Figure-9 were used as input for the aerospike analysis illustrated in Figure-10 where vacuum  thrust and specific impulse (Isp) are predicted to be 264,600 lbf and 454.7 sec where Pc/Pa = 100000. Dimensional data for this analysis are based on Boeing's results for the XRS-2200 linear aerospike rocket engine where vacuum thrust and Isp are 266,230 lbf and 436.5 sec for a 0.6% variation in thrust and 4.2% variation in Isp. Due to unavailable thruster dimensions and conflicting dimensional and performance information provided by NASA and relevant technical papers this aerospike rocket engine analysis is an approximation.

NASA Description of the X-33 Spaceplane: The X-33 was to have been a wedged-shaped subscale technology demonstrator prototype of a potential future Reusable Launch Vehicle (RLV) that Lockheed Martin dubbed VentureStar. The company hoped to develop VentureStar early this century. Through demonstration flight and ground research, NASA's X-33 program was to have provide the information needed for industry representatives such as Lockheed Martin to decide whether to proceed with the development of a full-scale, commercial RLV program. The X-33 design was based on a lifting body shape with two revolutionary linear aerospike rocket engines and a rugged metallic thermal protection system. The vehicle also was to have had lightweight components and fuel tanks built to conform to the vehicle's outer shape. Time between X-33 flights was planned to normally be seven days, but the program hoped to demonstrate a two-day turnaround between flights during the flight-test phase of the program. The X-33 was to have been an unpiloted vehicle that took off vertically like a rocket and landed horizontally like an airplane. It was planned to reach altitudes of up to 50 miles and high hypersonic speeds. The X-33 Program was managed by the Marshall Space flight Center and was planned to have been launched at a special launch site on Edwards Air Force Base. Technical problems with the composite liquid hydrogen tank resulted in the program being cancelled in February 2001.

XRS-2200 Linear Aerospike Engine
XRS-2200 Engine 5K ft Vacuum
Thrust, lbf 204,420 266,230
Specific Impulse, sec 339 436.5
Propellants Oxygen, Hydrogen
Mixture Ratio (O/H) 5.5
Chamber Pressure, psia 857
Cycle Gas Generator
Area Ratio (Ae/At) 58
Throttling, Percent Thrust 50 - 100
Dimensions, inches
Forward End
Aft End
Forward to Aft
----
134
wide x 90 long
42 wide x 90 long
90
 
Figure-8 aerospike engine side view.

Figure-9, XRS-2200 aerospike rocket engine description data repeated in Figure-10.

XRS-2200 Linear Aerospike Engine Results
Figure-10, XRS-2200 vacuum (Pc/Pa=100,000) analysis. Input data for the results in Figure-10 based on data from Figure-9 using AeroSpike version 2.6.0.5.

Aerospike Engine Results Compared 5K ft (Pc/Pa = 70.0621) Vacuum (Pc/Pa = 100,000)
Thrust, lbf Isp, sec Thrust, lbf Isp, sec
AeroSpike 3.1 155,500 267 264,600 454.7
Boeing XRS-2200 Results 204,420 339 266,230 436.5

BOUNDARY SHAPE: The outer boundary angle from the lip of the thruster cowl to the end of the external ramp increases as altitude increases. The pressure ratio (Pc/Pa) defines the extent to which the outer boundary expands as altitude and pressure ratio increase. For pressure ratio, Pc is the thruster chamber pressure and Pa is the local atmospheric pressure. This section compares the XRS-2200 expansion boundary shapes at 5K feet (Pc/Pa = 70.0621) and 50K feet (Pc/Pa = 509.21) determined by Navier Stokes and AeroSpike.

Aerospike outer boundary shape

 

AeroSpike Rocket Motor Designed Using AeroSpike 2.6 NEW

Aerospike hot fire jet and inset illustrating AeroSpike 2.6 results superimposed on aerospike.

This section for the design and test of an actual aerospike rocket motor/thruster displays the geometry, fabrication methods and photographic results for a model rocket scale aerospike rocket motor. Click here to see a 10-second 400KB QuickTime movie of an aerospike rocket motor designed using AeroSpike 3.1. Requires QuickTime from Apple Computer. This design continues to evolve so updated information will appear as results are available.

Aerospike Input Data, 1st column. Some Results, 2nd column
Thruster exit area ratio (Aei/At) 2.0   Chamber temperature 1047.0
Aerospike expansion ratio (Ae/At) 10.5   Chamber pressure 154.0
Ratio of specific heats for exhaust 1.2   Aerospike thrust 2.802
Gas constant of exhaust (Rgas) 247139   Ramp base radius (Rbase) 0.0
Aerospike pressure ratio (Pc/Pa) 10.5   CF - Thrust coefficient 0.888
Thruster internal circular radius (Ra) 0.2   CF - Vacuum thrust coefficient 1.748
Radius to lip of cowl (Re) 0.25   Isp - Specific Impulse 51.3
Aerospike length from origin (Lnozzle) 0.96   Isp - Vacuum specific impulse 101.0
Propellant weight (lbs) 0.0077   Estimated burn time (sec) 0.172


Aerospike generated using AeroSpike 3.1. This plot properly scaled was used as a template.

Aerospike Ramp Fabrication   Aerospike Rocket Motor: Exploded View
 
Aerospike fabricated from a solid bar of 6061-T6 aluminum using a lathe to machine the article using a template generated by AeroSpike 2.6.
 
  Exploded view of the aerospike rocket motor. Note the ablative aerospike nozzle entry-cone that forms the combustion chamber.
Aerospike Rocket Motor Test (Tburn = 0.0 sec)   Aerospike Rocket Motor Test (Tburn = 0.172 sec)
 
Aerospike rocket motor primed and ready for hot fire testing.   Aerospike image clipped from QuickTime movie. Requires QuickTime from Apple.

PART 2: 2-D & 3-D MINIMUM LENGTH NOZZLE DESIGN USING
THE METHOD OF CHARACTERISTICS (FREE BONUS ADDITION)
BACK TO TOP


Summary of Features

1. Determine shapes and flow properties of 2-D Minimum Length Nozzles (MLN) given exit Mach number (Mdesign) and throat diameter (Dt).
2. Determine shapes and flow properties of 3-D Minimum Length Nozzles using an approximation procedure based on 2-D results.
3. Click the UpDown command button to move a locator from point to point in the flow field.
4. All important flow properties are displayed in real time as the locator moves from point to point in the flow field described by the characteristic mesh.
5. Generate color contour plots of Mach number (Mn), Pressure (P/Pc), Temperature (T/Tc) and density (R/Rc) with a single click.
6. Units include, MKS (meter-newton-second), CGS (centimeter-dyne-second), FPS (foot-pound-second) and IPS (inch-pound-second).
7. Define gas properties for inert gases, liquid propellant gases and solid fuel propellant gases or insert your own values.
8. Output all flow variables to the printer or text file for use with spreadsheet applications.

Propellant Gases Available

Inert Gases
Dry Air Hydrogen Helium Water Vapor Argon Carbon Dioxide
Carbon Monoxide Nitrogen Oxygen Nitrogen Monoxide Nitrous Oxide Chlorine
Methane          
 
Liquid Fuel Propellant Gases
Oxygen, 75% Ethyl Alcohol(1.43) Oxygen, Hydrazine(.09) Oxygen, Hydrogen(4.02)
Oxygen, RP-1(2.56) Oxygen, UDMH(1.65) Fluorine, Hydrazine(2.3)
Fluorine, Hydrogen(7.60) Nitrogen Tetroxide, Hydrazine(1.34) Nitrogen Tetroxide, 50% UDMH, 50% Hydrazine(2.0)
Nitric Acid, RP-1(4.8) Nitric Acid, 50% UDMH, 50% Hydrazine(2.20)  

Solid Fuel Propellant Gases
Ammonium Nitrate, 11% Binder, 4-20% Mg Ammonium Perchlorate, 18% Binder, 4-20% Al Ammonium Perchlorate, 12% Binder, 4-20% Al

Hybrid Rocket Motor Propellant Gases
85% Nitrous Oxide, 15% HTPB    

User-Defined Gases
Specify custom or user-defined gases by inserting Ratio of specific heats (g) in the Minimum Length Nozzle Data section.

General Discussion
The Minimum Length Nozzle routine performs a minimum length nozzle (MLN) design using the method of characteristics. A minimum length nozzle has the smallest possible throat-to-exit length that is still capable of maintaining uniform supersonic flow at the exit. Strictly speaking a minimum length nozzle requires a sharp corner at the throat. However, sometimes a sharp corner at the throat may be impractical. A nearly minimum length nozzle may be generated by specifying a very small but finite radius of curvature at the throat with the inflection point of the throat-curve just downstream of the throat. For a nearly minimum length nozzle simply specify the streamline from the throat-curve so the curvature lines up with the nozzle wall shape generated by MLN.

A straight sonic line is assumed to occur at the throat of the minimum length nozzle. For the example presented in Figure 2 and Figure 3, where the exit Mach number is 2.4, the first characteristic (C_ ) propagating from the corner of the throat is inclined by a small amount (
q = 0.375 deg) from the normal sonic line. The slope of the first characteristic is dy/dx = (q - m) = -73.725 deg. See Figure 1 below. The remaining expansion fan is divided into six increments. The Mach number at each point in the flow is determined from the Prandtl-Meyer function using the Newton-Raphson iteration method and the unit processes dictated by the method of characteristics. The nozzle contour is drawn by starting at the throat corner where the maximum expansion angle of the wall, qw_max is equal to one-half the Prandtl-Meyer function, n(Mn) / 2, at the design exit Mach number. For a minimum length nozzle the maximum expansion angle is equal to one-half the Prandtl-Meyer function for the design exit Mach number. For other nozzles the maximum expansion angle must be less than n(Mn=Mdesign) / 2. For a detailed discussion of the method of characteristics please refer to the reference, Modern Compressible Flow, With Historical Perspective, by John D. Anderson, pages 260 to 282.

PLEASE NOTE: For the 2-D Minimum Length Nozzle selection the "X" and "Y" coordinates of the nozzle contour represent the horizontal and vertical dimensions that define the 2-D characteristic mesh. Therefore, the Exit Area Ratio (Aexit/Athroat) = [2*Yexit*WIDTH] / [2*Ythroat*WIDTH] = Yexit/Ythroat because the flow is 2-Dimensional. Likewise, for the 3-D Minimum Length Nozzle selection the "X" and "Y" coordinates of the contour represent the horizontal and radial dimensions that define the 3-D axisymmetric mesh. Therefore, the Exit Area Ratio (Aexit/Athroat) = [
p*Yexit^2] / [p*Ythroat^2] = (Yexit/Ythroat)^2 because the flow is 3-Dimensional and not 2-Dimensional. Finally, MLN Project files generated by previous versions of MLN must be updated by adding "1" for 2-D flow or "2"  for 3-D flow at the bottom of the MLN file. Do not forget to save each Project file with an MLN extension when updating older Project files for use with AeroSpike 2.5 or higher.

Procedure
From the menu on the top of the main start-up screen, select units (MNS, CGS, FPS or IPS) from the Units menu and then the propellant gas from the Gases menu. A number of inert gases, liquid propellants and solid fuel propellants are available. The value for the ratio of specific heats (
g) are determined from the Units and Gases menus and are passed on to the MLN program after clicking the Minimum Length Nozzle command button on the main start-up screen. Only the ratio of specific heats are required for the MLN analysis. The other values including the gas constant (Rgas), chamber pressure (Pc) and pressure ratio (PR) are not required for the MLN analysis. When performing an MLN analysis the only values required are the Inclination angle from the sonic line (see above), Design Mach number (Mdesign), and throat diameter (Dt). The ratio of specific heats has already been specified from the main screen selection. However, the user can over-ride the inserted ratio of specific heats by simply inserting his own ratio of specific heats in the data entry box. Each time the user changes any data entry the results are automatically updated and displayed. The user only needs to click the Plot button to see a new contour plot of the results or the UpDown button to see flow results at any of the characteristic mesh points.

Toolbar Operations

1. Show or hide the main data window from view or from being printed.
2. Send all flow propeties (X, Y, Mn etc) at each characteristic mesh point to the printer.
3. Send an image of the screen to the printer.
4. Save all flow properties (X,Y, Mn etc) at each characteristic mesh point to a data file.
5. Read the nozzle description file from a previous session.
6. Save the nozzle description file from a previous session.
7. Refresh the displayed analysis to the default analysis seen during start-up.
8. Return to the main start-up screen.

Input Variable Definitions
1. Ratio of specific heats (
g): Selected from pull-down menu or user-defined.
2. Inclination angle from sonic line: This angle is used to compute the slope of the first characteristic from the edge of the throat.
3. Design Mach number (Mdesign): The Mach number at the exit of the nozzle where the flow is uniform.
4. Throat diameter (Dt): The entrance to the minimum length nozzle where Mn = 1.0
5. Area ratio (Aexit/Athroat): The resulting exit area ratio of the nozzle determined by the method of characteristics.
6. Specify whether the nozzle is 2-D or 3-D by clicking either the 2-D characteristics or 3-D approximation option buttons.


Figure 1. Description of the inclination angle (q) from the sonic line (at throat) where Mn =1.0


Minimum Length Nozzle Validation-1
Example 11.1, on page 282 from Modern Compressible Flow, With Historical Perspective, by John D. Anderson

Figure 2. Example 11.1, on page 282 from Modern Compressible Flow, With Historical Perspective, by John D. Anderson

NOTE ABOUT MLN ANALYSIS ACCURACY: The Minimum Length Nozzle (MLN) analysis illustrated in Figure-2 uses exact input data from Modern Compressible Flow With Historical Perspective. For maximum accuracy simply insert 0.0 degrees in the Inclination angle from sonic line input block. When this simple modification is performed the exit Area ratio (AR) becomes 2.43 which compares to AR = 2.403 for exact 2-D isentropic flow and represents a 1.124% difference from isentropic theory.


Figure 3. MLN Characteristic Mesh For exit Mach number of 2.4 where Inclination angle from sonic line = 0.0 degrees produces AR = 2.43.


Minimum Length Nozzle Validation-2
Figure 17.5, Gasdynamics: Theory and Applications, 2-D and Approximate 3-D MLN Validation at M = 3.0
The following table compares 2-D and 3-D MLN results with data scaled from Figure 17.5, Gasdynamics: Theory and Applications. Two wall-points, one from the center and one at the end of the nozzle contour have been selected for comparison. Notice that 3-D Minimum Length Nozzles are substantially shorter than equivalent 2-D Minimum Length Nozzles that have identical
Area Ratio (Aexit/Athroat). For comparison purposes all results are referenced to the curved sonic line analysis for 2-D and 3-D axisymmetric nozzles. Please reference Gasdynamics: Theory and Applications* page 325, Figure 17.15 where g = 1.4 and Mexit = 3. Finally, please note that MLN uses the straight sonic line method of characteristics analysis.

2-D MLN Analysis X_Coordinate   Y_Coordinate Difference X_Coordinate Difference Y_Coordinate Difference
AeroSpike 3.1 6.522   3.331 -10.7% 17.43 4.79% 4.354 -6.8%
Straight Sonic Line* 6.522   3.30 -11.5% 16.83 -0.53% 4.198 -10.1%
Curved Sonic Line* 6.522   3.73 - 16.92 - 4.670 -
                 
3-D MLN Analysis X_Coordinate   Y_Coordinate Difference X_Coordinate Difference Y_Coordinate Difference
AeroSpike 3.1 3.126   1.603 -0.06% 8.353 -2.68% 2.087 2.9%
Axisymmetric
Curved Sonic Line*
3.126   1.604 - 8.59 - 2.028 -


Figure 4. 2-D MLN Characteristic Mesh when
g =1.4 and exit Mach number = 3


Figure 5. Approximate 3-D Axisymmetric MLN Characteristic Mesh when
g =1.4 and exit Mach number = 3
BACK TO TOP


AeroSpike System Requirements
(1) Screen resolution: 800 X 600
(2) System: Windows 98, XP, Vista, Windows 7, Windows 10 (32 bit and 64 bit), NT or Mac with emulation
(3) Processor Speed: Pentium 3 or 4
(4) Memory: 64 MB RAM
(5) English (United States) Language
(6) 256 colors

Please note this web page requires your browser to have
Symbol fonts to properly display Greek letters (
a, m, p, and w)

ADDITIONAL REQUIREMENT: Input data for all AeroRocket programs must use a period (.) and not a comma (,) and the computer must be set to the English (United States) language. For example, gas constant should be written as  Rgas = 355.4 (J / kg*K = m^2 / sec^2*K) and not Rgas = 355,4. The English (United States) language is set in the Control Panel by clicking Date, Time, Language and Regional Options then Regional and Language Options and finally by selecting English (United States). If periods are not used in all inputs and outputs the results will not be correct.


AEROSPIKE REVISIONS
AeroSpike 2.3 Features and Error Fixes

1. Fixed plot resolution problem that occurred for some high aspect ratio aerospike nozzles.
2. Fixed the ratio of specific heats (
g) manual entry error that would not accept gamma (g) = 2 and some other minor errors.
3. Fixed the incorrect gas constant (Rgas) value for hydrogen. For hydrogen, Rgas = 4122.11 m^2/(sec^2*K).
4. Added the ability to specify thruster sonic-section (throat) angle. Throat angle can vary from 60 to 120 degrees, the default is 90 degrees.

AeroSpike 2.4.1 Features and Error Fixes

1. Added a hybrid rocket motor propellant having the following fuel and oxidizer to the list of combustion gases: 85% Nitrous Oxide, 15% HTPB.
2.
Added the ability to save F(x) verses PR (Pressure Ratio) and F(x) verses x to a CSV file for use with Notepad or Excel.
3.
In the Aerospike Nozzle Data section added a display of Truncation as percent of total aerospike length.
4. In the Aerospike Nozzle Data section added a display of Distance from throat (origin) to end of thruster.
5. In the Aerospike Nozzle Data section added a display of Distance from
end of thruster to end of ramp.
6. Corrected a few Status Bar display errors for plots of F(x) verses x.

AeroSpike 2.4.2 Error Fix (11/26/2006)
1) The gas Nitrogen Dioxide in the Gases pull-down menu should be labeled Nitrous Oxide (N2O). (Fixed)

AeroSpike 2.5.0 Features (01/23/2007)
1. Added ability to determine shapes and flow properties of 3-D Minimum Length Nozzles using an approximation procedure based on 2-D results.

AeroSpike 2.6.0.1 Features (09/21/2008)
1. Added the ability to include base thrust of truncated aerospike nozzles to determine total thrust and thrust coefficient.
2. To modify or change units in previous versions of AeroSpike the user needed to close the main aerospike nozzle analysis screen and then redefine units in the start-up screen. However, starting with this new version the user goes directly to the start-up screen without closing the aerospike nozzle analysis screen to modify pressure ratio, units, gases and altitude to instantly included those changes on the main aerospike nozzle analysis screen.
3. Darkened all data display boxes to prevent confusion with white data entry boxes for the MLN and aerospike analyses.

AeroSpike 2.6.0.2 Features (09/14/2009)
1) For AeroSpike, fixed all input data text boxes for 32 bit and 64 bit Windows Vista. When operating earlier versions of AeroSpike in Windows Vista the input data text boxes failed to show their borders making it difficult to separate each input data field from adjacent input data fields.
This simple change did not alter any computational result.

AeroSpike 2.6.0.3 Features (04/27/2011)
1) For the Minimum Length Nozzle routine corrected the 3D axisymmetric mesh-data display.
This simple change did not alter any computational result.

AeroSpike 2.6.0.4 Features (09/21/2011)
1) For AeroSpike the color contour plots and aerospike nozzle shape data
are now displayed to scale. In previous versions aerospike color contour plots and aerospike nozzle shape displays were not to scale to allow full use of the available display area. However, recent work indicates it is more useful to display scaled aerospike geometry than to fill the entire plot area. The factor 0.5390 was used to properly scale the X-coordinates of the X, R plot data. This simple change did not alter any computational result.

AeroSpike 2.6.0.5 Features (09/25/2011)
1) For AeroSpike the Truncation as percent of total aerospike length
data field displayed incorrect values intermittently based on UNITs selection. This simple change did not alter any other computational result.

AeroSpike 3.1.01 Features (01/21/2021)
1) Increased numerical analysis speed for aerospike and minimum length nozzle (MLN) computations.
2) In the Gases pull down menu added properties for Liquid Oxygen and Liquid Methane (LOX/LMH4) gases for mass ratios; 2.70, 2.80 and 2.90.


For more information about
AeroSpike please contact AeroRocket.

BACK TO TOP


| MAIN PAGE | RESUME |
Copyright © 1999-2022 John Cipolla/AeroRocket