Rocket Motor Nozzle
Rocket Motor Nozzle
Rocket Motor Nozzle
IdeaExchange@UAkron
Williams Honors College, Honors Research The Dr. Gary B. and Pamela S. Williams Honors
Projects College
Spring 2020
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:
Anthony Smith
Joshua Meyers
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
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.
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:
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
6
Figure 1: Cone and Bell Shape Design
7
Figure 3: Solidworks Sketch of Nozzle Geometry
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.
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.
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.
11
Embodiment Design
Materials Research
Anticipated Temperatures at Nozzle Locations
Temperatures (K) at: Liquid
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
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.
14
Figure 11: Carrier Drawing
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
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.
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.
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.
22
The equivalent stress calculation correlates to the Ansys calculations. Now onto the throat insert.
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.
24
Figure 18: Second Printed 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
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
%% 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)
30
Signature Page
Senior Design Project
STUDENT(S)
Name: E-mail:
ADVISOR(S)
Name: E-mail:
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