Rocket Motor Nozzle

Download as pdf or txt
Download as pdf or txt
You are on page 1of 33

The University of Akron

IdeaExchange@UAkron

Williams Honors College, Honors Research The Dr. Gary B. and Pamela S. Williams Honors
Projects College

Spring 2020

Rocket Motor Nozzle


Corey Hillegass
[email protected]

Follow this and additional works at: https://ideaexchange.uakron.edu/honors_research_projects

Part of the Aerodynamics and Fluid Mechanics Commons, Aviation and Space Education Commons,
Computer-Aided Engineering and Design Commons, Heat Transfer, Combustion Commons,
Manufacturing Commons, Propulsion and Power Commons, and the Structures and Materials Commons
Please take a moment to share how this work helps you through this survey. Your feedback will
be important as we plan further development of our repository.

Recommended Citation
Hillegass, Corey, "Rocket Motor Nozzle" (2020). Williams Honors College, Honors Research Projects.
1135.
https://ideaexchange.uakron.edu/honors_research_projects/1135

This Dissertation/Thesis is brought to you for free and open access by The Dr. Gary B. and Pamela
S. Williams Honors College at IdeaExchange@UAkron, the institutional repository of The University
of Akron in Akron, Ohio, USA. It has been accepted for inclusion in Williams Honors College,
Honors Research Projects by an authorized administrator of IdeaExchange@UAkron. For more
information, please contact [email protected], [email protected].
Senior Design Final Report:

Rocket Motor Nozzle


Corey Hillegass

Anthony Smith

Joshua Meyers

College of Engineering - The University of Akron

Spring 2020
Executive Summary
A rocket nozzle was designed for the Akronauts Rocket Design Team specifically for a
new liquid propellant motor that is being tested. Compressible fluid flow calculations were
made to create an inner nozzle geometry to successfully accelerate hot gases from combustion to
generate a sufficient amount of thrust. The nozzle was then designed around this geometry in
order to recreate heritage nozzle assembly protocol dealing with solid fuel motors. Once a first
iteration of design was made, Matlab and Ansys Workbench were used to analyze and compare
computation results. After initial analysis, new iterations to the nozzle insert and carrier design
were made to combat thermal and structural effects of the high temperature and pressure
environment it was exposed to. Through this project, the team was able to develop the
groundwork necessary for self-sustainable rocket nozzle manufacturing for the design team.

1
Table of Contents

Executive Summary 1

Table of Contents 2

Introduction 3

Conceptual Design 3
Theory & Geometry 3
Calculating Throat and Exit Areas 4
Evolution of Design Concepts 8

Embodiment Design 12
Materials Research 12

Detail Design 13
Geometric Design and Material Selection 13
Estimated Costs 17
Simulation 18
Prototype 24

Conclusions 25

Works Cited 26

Appendix 27
Matlab Code 27

Signature Page 31

2
Introduction
The senior design team aimed to research, analyze and manufacture a rocket motor
nozzle for the Akronauts rocket design team. In the following pages the team will be showing the
development and process for development of this rocket nozzle. The intent for this project is to
be utilized on an experimental research and development project new to the Akronaut rocket
design team. This project is developing a liquid fueled motor to be used in a future competition.
The team has only used solid fuel for rocket propulsion, so this project is assisting in deployment
of new technologies to the team. The rocket motor nozzle previously had been an afterthought, it
was assumed to have been manufactured correctly and performed efficiently. Upon further
analysis, the team decided to tackle the project of designing and testing a new nozzle that would
work in conjunction with the liquid motor development to utilize techniques learned from the
classroom into actual application. The goal for this project is to understand the rocket motor
nozzle, design an efficient and manufacturable nozzle, and further the progression of the team.

Conceptual Design

Theory & Geometry


Nozzle dimensions are determined by the temperatures and pressures that are created by
the fuel used in the combustion chamber. The Akronauts have been using commercial and
student built solid fuel motors since the creation of the team but have recently been putting time
into researching the idea of building a liquid motor propulsion system. Data simulation software
such as Rocket Propulsion Analysis (for liquid fuel) were used to calculate the chamber
temperatures and pressures so that throat and exit parameters (pressures, temperatures, and areas)
can be calculated so that a nozzle can be optimized to atmospheric pressure at launch. The
reason that the decision was made to have a perfectly expanded nozzle upon launch is because
the rocket needs as much thrust at lift-off as possible. A higher velocity off the launch rail will

3
help the rocket stabilize itself and fly straight as compared to a slower rocket, which has more of
a possibility to ‘tip over’ and fly at an angle that is not perpendicular to the ground.
For nozzle shape, the project was centered around a bell nozzle design, which is different
from the student made nozzles by the Akronauts in the past. Past nozzle designs are based off of
a conical design with a half angle of 15°. This half angle is generally used in conically designed
nozzles because a higher or lower angle will lose thrust. The bell design is based on the conical
design when it comes to the throat and exit area, but has a parabolic contour from the throat to
the exit, as the name proposes. This difference being the bell is proven to minimize losses in
thrust that a conical contour cannot. The bell shape can have a lower exit angle from a two
dimensional standpoint which lowers the loss of thrust that is not going along the center axis of
the rocket. Other minor losses the bell shape mitigates the energy that flow expansion waves
have throughout the nozzle. An increasing, yet diminishing, rate of area that an exhaust gas
experiences from the throat to exit through a method of characteristics is to help minimize this
loss in energy.

Calculating Throat and Exit Areas


With the pressure, temperature, mass flow rate, and specific heat ratio taken from
Burnsim analysing the solid fuel, it is now possible to calculate throat area. The relationship
between the chamber pressure and the throat pressure is

Equation 1: Stagnation Pressure


where p0 is the chamber pressure, p is the throat pressure, and M is equal to 1 because mach 1 is
the velocity of exhaust at the throat. K is the specific heat ratio. Temperature at the throat is
found by using this equation

Equation 2: Stagnation Temperature

4
where T0 is the chamber pressure, T is the throat temperature and M is equal to one again. These
throat pressure and temperature values are then used the find throat area:

Equation 3: Throat Area


where p1 is throat pressure m dot is mass flow rate, R is the gas constant for the fuel used, and T1
is the throat temperature. The exit velocity can also be calculated using this formula:

Equation 4: Exit Velocity


where T1 is the chamber temperature, P1 is the chamber pressure and P2 is the exit pressure. The
exit pressure for an ideally expanded nozzle is equal to the ambient temperature the nozzle is
built for, which in this case, its sea-level conditions. In order to find the expansion ratio, Mach
number is needed. And in order to find mach number, the exit temperature needs to be
calculated:

Equation 5: Exit Temperature

Equation 6: Mach Number

Equation 7: Expansion Ratio


Multiplying the throat area to the expansion ratio value gives you the exit area. The final
dimension needed to design the nozzle contour is the length from throat to exit. This dimension
is based on the length of the conical nozzle. A 100% bell nozzle is the same length as a conical
nozzle, which is calculated by

Equation 8: Nozzle Length

5
where r2 is the exit radius, rt is the throat radius, and alpha is the half angle of the nozzle (ideal
half angle for a conical nozzle is 15 degrees). Usually, the bell nozzle is constructed to around
80% of the length of the cone length because you hardly lose any efficiency doing so. In
addition to decreasing length, you are decreasing weight and material cost. In this application,
the choice was made because the team wanted the most efficiency possible in this section of the
design, therefore a 100% bell nozzle is the configuration moving forward.
Table 1 shows parameters calculated by Matlab. The code with input parameters is in the
appendix.

Liquid

Expansion Ratio 3.94

Throat Diameter (in) 0.776

Exit Diameter (in) 1.5376

Length (in) 1.4235


Table 1: Nozzle Dimensions
For the shape of the nozzle, the Rao nozzle approximation method was utilized to design the
nozzle. The initial inlet contour of the nozzle is a radius of 1.5 times the throat diameter until you
reach the actual throat, then the contour is an arc with a radius of 0.382 times the radius of the
throat. The contour becomes a parabolic spline up to the exit. 𝛳i and 𝛳n are the initial and final
parabola angles in the nozzle design and can be found in this graph from the Rocket Propulsion
Elements textbook which is based off of expansion ratio and nozzle length.

6
Figure 1: Cone and Bell Shape Design

Figure 2: RAO Nozzle angles vs. Expansion Ratio Plot

7
Figure 3: Solidworks Sketch of Nozzle Geometry

Evolution of Design Concepts


With past solid rocket motor designs, the interface of the nozzle to the motor casing was
relatively simple. The nozzle would be inserted in a machined housing and held in place by
retaining rings shown in Figure 4. Due to the design of this nozzle being for a liquid propelled
motor, the design had to be reconstructed into two subsystems: a combustion chamber and a
nozzle.

Figure 4: Nozzle Post Burn (Solid Fuel)

8
There were guidelines regarding flow geometry based on legacy experience and proven
concepts. However, for this particular purpose, there isn’t a “standard” design that is used
regarding how to bring a nozzle to life. This must be created nearly from scratch and adapted as
efficiently as possible to a given application. Small scale liquid rocket motors,of course, have
been developed and successfully used in the past, so the team was able to use these past designs
by both colleges and hobbyists alike as input to develop the official Akronauts design. As seen
in Figure 5, which is along the lines of one of the original design ideas, the portion of the nozzle
which is located downflow of the throat has exposed geometry, with a carrier ring to hold it in
place. Based on extensive material research, it was deemed that the green portion of the nozzle in
the figure below would likely pose a risk in the area of structural integrity if graphite is used as
the insert/flowpath material.

Figure 5: Nozzle/Combustion Chamber Design Concept


Looking forward, the design was then geared towards having a carrier that fully
encapsulated the graphite portion of the nozzle. The purpose of using graphite compared to other
options is for reasons of cost, heat tolerance, weight, and machinability. As shown in Figures 6 &
7, the design was geared towards a bimaterial flowpath. This would allow the highest
temperatures to be dealt with by the graphite while downstream the flow would be much cooler

9
(relatively speaking). However, even though a metallic material would likely be able to
withstand the cooler portion of the post-throat burn, there are anticipated issues that would arise
regarding thermal expansion. If the 2 materials expand and contract at different rates, the
interfacing geometry where the materials change on the flow path may become uneven, which
would pose unnecessary issues. On top of this, since metallic materials are inherently more
conductive, exposing the lower portion to the burn would likely heat up the entire assembly,
which could eventually cause structural integrity issues resulting from repeated thermal shock.

Figure 6: Original Concept

Figure 7: Improved Concept


In the above figures, it can be seen that sealing solutions are starting to be incorporated
into the nozzle. To ensure pressure & thrust is only dispersed where intended, o-rings were added
between the diametral interface of the insert and carrier. A face seal was also added on the flange
of the assembly to ensure there will be no leakage from the combustion chamber to the nozzle.
Moving forward, for the reasons stated above, the goal swiftly moved to isolating the metallic
material from as much heat as possible. With this goal in mind, the resulting design eliminated a
bimaterial flowpath, and utilized the existing geometry of the nozzle to extend graphite insert as

10
the entire flowpath. As seen in Figure 8, the insert is held captive by the metallic carrier. Based
on Ansys simulations, this proved extremely effective in lowering the overall temperature of the
nozzle.

Figure 8: Graphite Insert Flowpath Design.


After the above design change was incorporated into Ansys simulations, it was
discovered that temperatures were still marginal around the flange, flowing downward
throughout the nozzle. To remedy the high temperatures, a phenolic insulator was added to the
upward part of the nozzle to essentially eliminate a conductive path between the graphite and
metallic carrier. One o-ring was removed from the design, leaving only 2 necessary o-rings in the
entire design. These changes can be seen in Figure 9. This configuration was the final resulting
design of many months of research, and the interface geometry between the two materials was
able to be manipulated to reduce temperatures to a safe range, while maintaining large factors of
safety throughout the design.

Figure 9: Final Design with Phenolic Liner

11
Embodiment Design

Materials Research
Anticipated Temperatures at Nozzle Locations
Temperatures (K) at: Liquid

Combustion Chamber 2964

Throat 2711

Exit 1631
Table 2: Anticipated Temperatures at Nozzle Locations
Considering the extremely high temperatures that are experienced during a burn, our
choices for material are limited and crucial. In the past, graphite has been used for nozzle
material. An on-going issue with this is that each nozzle is only good for one launch due to
excessive wear, which altars the flow geometry. Even though liquid fuels reduce wear, research
and legacy experience indicates that wear will still occur in graphite nozzles during each burn.
Graphite may be used for testing, but once a design is optimized, it will likely turn into a more
costly option, especially if a heat exchanging system is incorporated into the design.
Ideally, 3D printing could be used to produce a nozzle, and the possibility of printing
inconel was heavily considered. However, the 3D printing process involves melting a material so
that it can be manipulated into a desired geometry. Being that our nozzle will be exposed to
temperatures in excess of 3000 °K, it is unlikely to find a material to withstand required
temperatures while being affordably printed. Thus, machining a selected material utilizing CNC
technology has become the preferred option.
OFHC (Oxygen Free High Conductivity Copper), stainless steel, carbon fiber/other
composites, ceramics were looked into as options. One of the bigger challenges that the team is
running into is finding a material that will be able to withstand the required temperatures. It is
undetermined as to how beneficial preheating the oxidizer/fuel will be. Utilizing regenerative
cooling with a high conductivity material such as steel or copper will be a riskier path to take

12
since any failure of this system could result in melting, so unless deemed necessary, this path will
likely not be taken. Carbon fiber and other composites appear to be impractical and costly, so
those are cut from the list of possible options.
Due to high temperature endurance and manufacturability, graphite and ceramics were
explored as the most practical options. In the automotive industry, ceramic coatings over steel
exhaust components are becoming more popular. Combining the machinability and cost
effectiveness of steel with the high temperature endurance of certain ceramics could deal with
the heat dilemma. This would involve machining the nozzle out of steel and then coating the
critical flow surfaces with a ceramic compound.
While it may be a simpler and “safer” option to make the entire nozzle out of a ceramic
such as TaC (tantalum carbide), this may not be practical due to cost. One of the primary topics
of research that has been discussed is the primary structure of the nozzle being machined out of
steel or aluminum, while a thin layer of UHTC (Ultra High Temperature Ceramics) could be
formed and used as an insert so that the entire critical flow surface would be made of the
ceramic. This would provide temperature resistance in the necessary areas while maintaining
machinability and reduce cost with the nozzle “shell”.

Detail Design

Geometric Design and Material Selection


Considering factors such as cost, manufacturability, and reliability, the decided design of
the nozzle is a 2 piece system (2 main pieces) comprised of a grade GT-62B compression molded
babbit impregnated carbon-graphite insert, a 310 stainless steel carrier, a phenolic insulating
liner, and viton o-rings for sealing. The specific grade of graphite was picked based on legacy
experience and the 310 stainless steel was chosen over aluminum primarily due to its superior
high temperature tolerance under working conditions. Both aluminum or 310 stainless are

13
relatively cost effective options, however at this stage of the liquid rocket motor program, the
stainless will yield a much higher factor of safety compared to aluminum. It is, however, possible
for a switch to aluminum to be made once reliable test data is obtained and the design is
reanalyzed. Aluminum is much lighter, which is ideal for flight. Figures 10, 11, and 12 show
mechanical drawings of the insert, carrier, and over all nozzle assembly. Table 3 shows the bill
of materials for assembly.

Figure 10: Insert Drawing

14
Figure 11: Carrier Drawing

Table 3: Bill of Materials

15
Figure 12: Assembly Drawing
Observing the assembly drawing in Figure 12 and the insert drawing in Figure 10, it can
be seen that the flow path on the final design is 100% graphite. Compared to earlier designs, this
is much simpler and more reliable due to the fact that a bi-material flow path would likely
expand at dissimilar raters and any potential heat spike or failure could cause the metal to melt.
The current design utilizes the nozzles geometry to hold the entire design captive with 6 grade 8
¼-20 steel alloy bolts. As seen in the assembly drawing, the phenolic liner isolates the “upper”
portion of the nozzle, which is where thermal simulations showed to be the hottest part of the
geometry (inside the flowpath). Since the phenolic liner has such a low thermal conductivity, it
insulated the hottest part of the graphite from the stainless carrier; thus, the overall temperature
of the carrier was reduced since such a large conductive path for heat was essentially eliminated.
Due to the tight tolerancing on the “lower” part of the nozzle, the insert should. remain centered
in the carrier with no structural issues in the hoop direction. The selected o-rings are made of
viton, which has proven itself via legacy experience for this application. Since this material also

16
has a very low thermal conductivity, there should be no issues with the integrity of the o-rings.
They will, however, be changed after each burn since they will be exposed to such a severe
thermal shock.
Regarding the drafting principles used in the technical drawings, they were designed with
manufactuability, along with machining capability in mind. Since these components could be
machined on a restricted time schedule and most college students don’t have a fluent knowledge
of GD&T (Geometric Dimensioning and Tolerancing), it was not heavily incorporated into the
drawings. That being said, it wouldn’t be hard, and in a high volume industrial situation, true
position, circular runout, and total runout tolerances (with respect to the appropriate datums)
would have likely been incorporated into the drawings. This would hold tighter and truer
tolerances on the finished product. However, another factor that the team had to consider from a
practicality stand point was inspection. For certain geometric tolerancing details, the design team
does not have the means of inspecting things (accurately) to such a level of precision. That being
said, the drawings were drafted with the current manufacturing and inspection capabilities of the
team in mind.

Estimated Costs
Based on our design, the following table outlines the approximate material costs. The
actual costs vary based on machine availability and if new manufacturing machinery needed to
be acquired. Utilizing the University of Akron machining shop for most construction, the
Akronauts rocket design team can lower these costs and assume the proper machinability is able
to be achieved or achieved through company connections. For the purposes of this project, it is
assumed that the hours used for the construction and development of the design are being done
by the team members. For advanced machining, an outside resource with a cost to be determined
would be utilized, possibly in return for sponsorship of the project. For this, the costs can be
whittled down to material acquisition. The cost associated with the phenolic liner is tasked to the
combustion chamber design group as this part is originating from the combustion chamber and
not from the rocket nozzle.

17
Part Quantity Cost Store

310 Stainless Steel Rod 4”x6” 1 ~$410 McMaster-Carr

Molded Carbon-Graphite Rod 1 $229.22 graphitestore.com


[~$57 per nozzle]
[2.5” x 12” good for 4 nozzles]

Viton O-ring 330 10-pack 1 $15.75 McMaster-Carr


[$1.58 per ring]

Viton O-Ring 233 10-pack 1 $14.39 McMaster-Carr


[$1.44 per ring]
Table 3: Bill of Materials Cost

Simulation
Engineering tools like Ansys Workbench were used in order to simulate the environment
that the nozzle would be exposed and asked to withstand. Ansys Fluent was used specifically to
output data to compare to the isentropic flow calculations that are in textbooks like Rocket
Propulsion Elements and Introduction to Compressible Fluid Flow. Once that simulation was
run and checked, the next step was to model and simulate the stresses on the carrier and graphite
throat insert. This was done using the Static Structural tool in Ansys. Once that was effectively
modeled, the model was then put into Thermal Transient to analyze how the system would react
to a five second burn time.
For Fluent analysis, a 2D model and a 3D model were made to analyze the model. The
job of the 2D model was more for observing the reaction of flow through a converging-diverging
nozzle. The model also had an absent area behind the nozzle exit in order for Ansys to simulate
the exhaust plume that would be generated in real life.
For the fluent simulation, a surface was modeled in Solidworks. The surface plane was
modeled as only half of the nozzle along the flowpath of the exhaust. If you look closely at the
Figure 13, you can see a centerline running though the plane surface. Once the flow was
simulated, symmetry was used to create a simpler simulation. Figure 13 is a contour of velocity
through a nozzle.

18
The velocity reacts as it should, with the velocity accelerating to the throat, the flow
being choked at the throat and diverging section further accelerating the flow to increase thrust.
The transition from a teal to bright green in the throat is where the flow experiences mach 1
velocities. Further down the stream in the exhaust, you can observe boundaries where velocities
suddenly decrease. These are oblique shock waves that are typical in compressible flow exiting
a nozzle. The simulation included a Spalart-Allmaras turbulence instead of inviscid
characteristics.

Figure 13: 2D Simulation with Plume model


For the 3D model, simulations were conducted: an inviscid model, and a
Spalart-Allmaras model. This time only simulating the flow within the nozzle as opposed to the
nozzle flow and plume. With the given chamber pressure and temperature, the simulation can be
used to output important data in the throat and exit to confirm calculations in Matlab. The
3-dimensional model was sliced into thirteen planes with a distance of 5 mm in between each
plane placed along the flow path. Using the surface integral tool in the Results section of Fluent,
data like velocity, pressure, temperature, and mach number can be averaged at these
cross-sectional planes and recorded. The data was then compared to values calculated in Matlab
using isentropic flow equations in the plots below.

19
The simulations behave like the isentropic flow equations and like the charts in Rocket
Propulsion Elements. The data from Matlab and the two simulations are in Table 4 below.
Values ending in t are throat values and exit values end in e.

Table 4: Values of important Throat and Exit Parameters

The next simulation was done in Thermal Transient, where the nozzle assembly is tested
to see if the material could withstand a 5 second burn. The Fluent data that was taken and
applied to the same sliced planes described earlier with the slicing of the Structural model. There
are 14 varying temperatures applied to these slices in the along the flowpath walls to simulate the
stagnation temperatures the nozzle will experience during a burn. Then and ambient temperature
of 293.15 °K was applied to the surfaces experiencing ambient temperature throughout the
burntime

20
This simulation led to some design changes to further increase the robustness of the
nozzle. The original design of the part had no phenolic liner extruded into the carrier but after
simulations, it was evident that the design would melt and fail if the Akronauts would follow
through with it and run a test. The inner ridge that is exposed to the temperatures of the
combustion chamber would have experienced temperatures of greater than 2000 °K. This led to
the decision to extrude the phenolic liner one inch into the nozzle assembly that was originally
going to line just the walls of the combustion chamber. This lead to also shortening the inlet,
which would thicken the thinnest part of the nozzle insert, where the combustion chamber meets
the nozzle insert. After the changes were made, the carrier would be shielded from the large
amount of heat from the hottest part of the engine. The most at-risk part of the nozzle assembly
became the exit retention ring of the carrier. The temperature contour is shown in Figure 14.

.
Figure 14: Temperature Contour of Nozzle Assembly
According to simulation, the maximum temperature experienced by the carrier will be 1229 °K
at the point shown in Figure 14. With this being the case, this data led to the decision that 310
stainless steel would make for the best carrier based on its thermal properties like a melting point
ranging around 1700 °K.
In a Static Structural analysis both the throat insert and carrier were analyzed in separate
simulations. The carrier was subjected to a 350 psi pressure load to the internal surfaces and was
treated as a pressure vessel. The entire carrier wouldn’t be subjected to the pressures as high as

21
the chamber all throughout the nozzle length was subjected to it anyways to ensure that the
design could withstand the maximum pressure. The results of the analysis are in Figure 15.

Figure 15: Von Mises Equivalent Stress Contour of Carrier


The contour above is a Von Mises equivalent stress contour that displays the distribution
of stress on the part. As predicted, the most stress is concentrated in the mid-section of the
carrier between the lip of the exit end and the bolting flange of the nozzle. The maximum stress
experienced in the part will be 17.8 MPa along the inside wall of the part according to the
program. On the outside wall, the part will experience stresses of 13 MPa. In order to confirm
this, a simple hoop stress calculation can be done to ensure the validity of the Ansys model. The
equation for hoop stress for a cylindrical pressure vessel is
PR
σ= t
Equation 9: Equivalent Stress for a Cylindrical Pressure Vessel
Where P is the internal pressure of the vessel, R is the effective radius of the vessel, and t is the
thickness of the vessel wall. With P = 350 psi, R = 1.4 in, and t = 0.2 in, the equivalent stress of
the carrier wall becomes
350*1.4
σ= 0.2
= 2450 psi
2450 psi = 16.9 M P a
13 M P a < 16.9 M P a < 17.8 M P a

22
The equivalent stress calculation correlates to the Ansys calculations. Now onto the throat insert.

Figure 16: Von Mises Stress Contour of Nozzle Insert


The insert was subjected to varying pressures throughout the length of the nozzle. The part was
sliced like the flow geometry in Fluent, and those pressures at each plane were applied to the
flowpath walls. The highest stress along the throat is right at the base of the fillet on the inner
rim of the nozzle opening. After evaluation, this makes sense because it is the thinnest part of
the throat design while being exposed to the chamber pressure. The hoop stress was evaluated at
this point.
350*.965
σ= 0.3196
= 1057 psi
1057 psi = 7.29 M P a
Comparing this to the simulation, 7.29 MPa is a very close estimation to the maximum stress
outputted by the model, which is 7.59 MPa.
The factor of safety can now be calculated for each part. For the carrier, the stress
experienced by the part is estimated to be around 2450 psi. If the yield strength of 310 stainless
steel is 33000 psi, then the factor of safety for the carrier would be around 13.5. The carrier
design is more than enough to withstand the pressure experienced during a burn. For the throat
insert, the maximum stress experienced is estimated to be 1057 psi. If the yield strength of
graphite is 7500 psi, then the factor of safety will be around 7. The throat design will be safe to
withstand a burn under these conditions.

23
Prototype
For the purposes of both the presentation and to have a physical model to look and touch, the
team decided to undertake rapid prototyping measures for the rocket motor nozzle. The simplest way to
achieve having a prototype constructed was through 3D printing. In total, 3 models were 3D printed and
assembled. These models, pictured below, are all created using 3D modeling software of SolidWorks.
They were saved out to .stl format files and input into a 3D slicer program Cura. Once sliced and
formatted, they were printed using 1.75mm gold, white, and blue PLA by a Creality Ender 3 Pro Printer.
This was done by one of the team members and without utilization of campus materials due to the
University shut down. The 1st model was created early on and helped the team reconsider the design,
proving to be a worthwhile undertaking. By looking at the design of Figure 18, the team redesigned the
insert and lengthened the carrier. For the final design in Figure 19, the phenolic insert was also included to
complete the design process.

Figure 17: Conceptual Design

24
Figure 18: Second Printed Design

Figure 19: Final Design

Conclusions
Further design reviews by the Akronauts Rocket Design Team would be the next step in
the design process of the nozzle. The steps that follow would be to submit these engineering
drawings to a machine shop and initiate the manufacturing process. Based on data analysis done
in Rocket Propellant Analysis, Matlab, and Ansys Workbench, the chosen fuel/Oxidizer mix of
ethanol and nitrous oxide was analyzed, the nozzle geometry was created, the nozzle assembly
was modeled, and the nozzle was confirmed and simulated. The overall design was invented
with legacy experience and new techniques learned in order to expand the groups engineering
skill sets. Based on the research and analysis of the design, the group feels that the design can
move forward and be trusted to run safely and solve the problem of a new, creative nozzle
design.

25
Works Cited

Mehta, R C, and K Suresh. Thermal Stress Analysis of a Solid Rocket Nozzle Throat Insert Using
Finite Element Method. Indian Journal of Engineering & Materials Sciences, 1998, pp.
272–277, Thermal Stress Analysis of a Solid Rocket Nozzle Throat Insert Using Finite
Element Method.
“Molded Carbon-Graphite Rod, Babbit, 2.5‘DIA x 12‘L.” GraphiteStore,
www.graphitestore.com/molded-carbon-graphite-rod-babbit-filled-mp001145.
Oosthuizen, P. H., and William E. Carscallen. Introduction to Compressible Fluid Flow. CRC
Press/Taylor & Francis Group, 2014.
Simonis, Rokas. “Home.” Rokas Simonis Engineering Portfolio, 2019,
simonis.coventry.domains/cfd/axisymmetric-simulation-of-a-rocket-mortar-exhaust-in-an
sys-fluent-dissertation-project-2/.
Sutton, George Paul, and Oscar Biblarz. Rocket Propulsion Elements. Wiley., 2017.

Rao, G. V. R. “Recent Developments in Rocket Nozzle Configurations.” ARS Journal, vol. 31,
no. 11, 1961, pp. 1488–1494., doi:10.2514/8.5837.
“What Is Von Mises Stress?.” What Is Von Mises Stress? - SimScale Documentation,
www.simscale.com/docs/content/simwiki/fea/what-is-von-mises-stress.html.
“Tantalum Carbide Nanopowder/ Nanoparticles” SSNano.com,
<https://www.ssnano.com/inc/sdetail/tantalum_carbide_nanopowder__nanoparticles___ta
c__99_5____1um_/23487>
Sutton, George Paul, and Oscar Biblarz. Rocket Propulsion Elements. Wiley., 2017.

26
Appendix

Matlab Code
clear
clc
close all

m0=24; % total Rocket mass (kg) (70 lbs)


g0=9.81; % gravity (m/s^2)
theta=90*pi/180; % launch angle (radians)
cd=.3;
rho=1.225; % avg density (kg/m^3)
A=pi*0.15621^2;

%% Liquid motor specs


mo=.4; %Mass flow rate Nitrous Oxide (kg/s)
mf=mo/4; %mass flow rate Ethanol (kg/s)
m=mo+mf; %Total Mass flow rate (kg/s)

mp=5; %Propellent mass (kg)


r=4; %Fuel ratio
M=.025076; % (Mf+r*Mo)/(r+1); %Molecular weight of
Propellent (kg/mol)
%
% acp=(cpf+r*cpo)/(r+1); %average Specific Heat (Cp) (KJ/kg-k)
% acv=(cvf+r*cvo)/(r+1); %average specific heat (Cv) (KJ/kg-k)
% k=acp/acv; %average specific heat ratio (Cp/Cv)
kc=1.175; %Specific Heat Ratio Combustion Chamber
kt=1.187; %Specific Heat Ratio Throat
ke=1.2425; %Specific Heat Ratio Exit
k=(kt+ke)/2;%
R= 8.314/M; %Universal gas constant KJ/kg-mol-k
%k1w=(kc+kt+ke)/2
Pc=2.413e+6; % Chamber pressure (pa) (350 psi)
Tc=2964.7; % Chamber Pressure K
Pa=101325; %Ambient Pressure (Pa)
Pe=Pa; %Exit Pressure (Pa)
%% Throat Parameters
Pt=Pc*(2/(kt+1))^(kt/(kt-1)); %Throat Pressure
Tt=2*Tc/(kt+1); %Throat Temperature
% Dt=.0201; %Throat Diameter (m)
% At=pi*(Dt/2)^2; %Throat Area (m^2)
Mt=sqrt(Tt*R*kt)
%% Fuel parameters

27
At=m/(Pc*kt*sqrt((2/(kt+1))^((kt+1)/(kt-1)))/sqrt(kt*R*Tt));%(m/Pt)*sq
rt((R*Tt)/(kt*(2/(kt+1))^((kt+1)/(kt-1)))); %Mass Flow Rate (kg/s)
cr=9; %Chamber Area Ratio (Ac/At)
Ac=cr*At; %Chamber Area
Dc=2*sqrt(Ac/pi); %Exit Diameter (m)
cp=kc*R/(kc-1);
%% Exit Conditions
Te=Tt*(Pe/Pt)^((ke-1)/ke); %Exit Temperature (K)
Ve=sqrt(((2*ke*R*Tc)/(ke-1))*(1-((Pe/Pc)^((ke-1)/ke)))); %Exhaust Gas
Velocity
a=sqrt(ke*R*Te); %Speed of sound at exit
Nm=Ve/a; %Exit Mach Number
%% Exit Parameters
Mo1=1:.01:Nm;
E=(1/Nm)*((1+(k-1)/2*Nm^2)/((k+1)/2))^((k+1)/(2*(k-1))); %Expansion
Ratio
Ae=E*At; %Exit Area
De=(2*sqrt(Ae/pi))*39.37; %Exit Diameter
Dt=(2*sqrt(At/pi))*39.37; %Throat Diameter
Lc=((De/2-Dt/2)/(tan(15*pi/180))); % Nozzle Length (m)

E1=(1./Mo1).*((1+(k-1)/2.*Mo1.^2)/((k+1)/2)).^((k+1)/(2*(k-1)));
figure(1)
plot(Mo1,E1)
title('Mach to Area ratio')
grid on
%% Thrust
Cf=
sqrt(((2*k^2)/(k-1))*((2/(k+1))^((k+1)/(k-1)))*(1-(Pe/Pc)^((k-1)/k
)))+((Pe-Pa)/Pc)*E;% Thrust Coefficient (dimensionless)
F=Cf*At*Pc; %Thrust
Is=F/(m*g0); %Specific Impulse (sec)
t=mp/m; % Burn time (s)
MR=(m0-mp)/m0;
F1=m*Ve

%% Launch Rail Velocity


a0=g0*((F/(g0*m0))-1); %Rocket initial acceleration (m/s^2)
vlr=sqrt(2*a0*3.6576); %Launch Rail Velocity (m/s) (22.86m/s = 75ft/s)

%% Altitude
upy0=Ve*log(1/MR)*sin(theta)-g0*t; %Velocity in y direction no drag
(m/s)
upx=Ve*log(1/MR)*cos(theta); %Velocity in x direction n
thetap=atan(upy0/upx)*180/pi;
B=.5*rho*upy0^2/MR;
upy=upy0-B*cd*A/m0;

28
yp=Ve*t*(1-log(1/MR)/(1/MR-1))*sin(theta)-.5*g0*t^2;
xp=Ve*t*(1-log(1/MR)/(1/MR-1))*cos(theta);

yz0=.5*(upy^2)/g0;

yz=yp+yz0;
% Altitude calculation for cesaroni L1350 (1349N thrust, .688kg/s mfr)
is
% 11% greater than OpenRocket's alt calculation
%% Flow Temperature
Mo=linspace(0,Nm);
T=Tc./(1+((ke-1)/2)*Mo.^2);
figure(2)
plot(Mo,T)
grid on
xlabel('Mach Number');ylabel('Temp (K)')
%% Flow Pressure
P=Pc./(1+((ke-1)/2)*Mo.^2).^(ke/(ke-1));
figure(3)
plot(Mo,P)
ylabel('pressure')
%% Printed Variables
%No descriptions will be given, just print functions
%Values Required for Design of Motor
fprintf('Design Values \n')
%fprintf('Pressure at Throat in Pa =
%4.4f\n',Pt)
fprintf('Temperature at Throat in K =
%4.4f\n',Tt)
fprintf('Throat Diameter (in) =
%4.4f\n',Dt)
fprintf('Exit Diameter (in) =
%4.4f\n',De)
fprintf('Nozzle Length (in) =
%4.4f\n',Lc)
fprintf('Expansion Ratio (Ae/At) = %4.2f\n',E)
fprintf('Thrust (N) = %4.2f\n',F)
%fprintf('Altitude (m)=
%4.2f\n',yz)
fprintf('Burn Time (s) = %4.2f\n',t)
fprintf('Specific Impulse (s) =
%4.2f\n',Is)
%fprintf('Launch Acceleration (m/s^2) =
%4.2f\n',a0)
fprintf('Launch Rail Exit Velocity(m/s) =
%4.2f\n',vlr)
%Values Used For Reference and Debugging
fprintf('Calculated Values \n')

29
fprintf('R in J/kgK %4.4f\n',R)
% fprintf('Cp in J/molK = %s\n',Cp)
% fprintf('Cv in J/kgK = %s\n',Cv)

fprintf('Temperature at Exit (K) = %f\n',Te)


fprintf('Area of Throat (m^2) = %f\n',At)
fprintf('Area of Exit (m^2) = %f\n',Ae)
fprintf('Local Sonic Velocity at Exit (m/s) = %f\n',a)
fprintf('Velocity at Exit (m/s) =
%4.4f\n',Ve)
fprintf('Mach Number at Exit =
%4.4f\n',Nm)

30
Signature Page
Senior Design Project

STUDENT(S)

Name: E-mail:

Anthony Smith [email protected]

Corey Hillegass [email protected]

Joshua Meyers [email protected]

ADVISOR(S)

Name: E-mail:

Dr. Mani Kannan [email protected]

Dr. Scott Sawyer [email protected]

FINAL REPORT

The final report and project are submitted to the advisor with following comments.

Comments:

Thank you for your assistance in the completion of the teams’ senior design project.

Advisor Signature:

Date:5/3/2020

31

You might also like