Project 2

PDF Created with deskPDF PDF Writer - Trial :: http://www.docudesk.com

Build Unsteady Model of “Pike” .. Use Integrator of your choice

Calculate:
Chamber pressure profile
Regression rate profile
Massflow rate (compare to choking massflow)
Mass depletion vs time
plot Thrust profile
plot Total Impulse profile
Effective Mean Specific Impulse
Allow:
St. Robert’s Parameter Input
Variable Step Size
Variable Thermodynamic Properties (as inputs to the problem)
Erosive burn model for cylindrical port (Not Bates grain)

1
MAE 5540 - Propulsion Systems
 Project 2 (2)
PDF Created with deskPDF PDF Writer - Trial :: http://www.docudesk.com

Part 1, cylindrical port

Fuel Grain Geometry
L0 = 35 cm Animal WorksTM,
D0= 7.6 cm L700 Motor Geometry
D0=3 cm
ρ propellant = 1260 kg/M3

Single propellant
Nozzle Geometry segment
A* = 1.887 cm 2
Aexit/A* = 4.0
θexit = 20 deg.

Assume ends are burn inhibited

2
MAE 5540 - Propulsion Systems
 Project 2 (3)
PDF Created with deskPDF PDF Writer - Trial :: http://www.docudesk.com

Combustion Gas Properties

γ = 1.18
MW = 23 kg/kg-mol
T0= 2900 K

Burn Parameters
a= 0.12 cm/(sec-kPa n)
n=0.16
M crit = 0.3
k = 0.2
(cylindrical port only)

3
MAE 5540 - Propulsion Systems
 Project 2 (4)
PDF Created with deskPDF PDF Writer - Trial :: http://www.docudesk.com

Part 1, cylindrical port

Fuel Grain Geometry
L0 = 35 cm Animal WorksTM,
D0= 7.6 cm L700 Motor Geometry
D0=3 cm
ρ propellant = 1260 kg/M3
Repeat results
Nozzle Geometry Using bates grain
A* = 1.887 cm 2 With 3 segments
Aexit/A* = 4.0
θexit = 20 deg.

Assume ends are not! burn inhibited

4
MAE 5540 - Propulsion Systems
 Project 2 (5)
PDF Created with deskPDF PDF Writer - Trial :: http://www.docudesk.com

Combustion Gas Properties

γ = 1.18
MW = 23 kg/kg-mol
T0= 2900 K

Burn Parameters
a= 0.12 cm/(sec-kPa n)
n=0.16
M crit = 0.3
Set to Zero for
k = 0.2
Bates grain
(cylindrical port only)
Assume no erosive’
burning
5
MAE 5540 - Propulsion Systems
 Project 2 (6)
PDF Created with deskPDF PDF Writer - Trial :: http://www.docudesk.com

Examine sensitivity of calculations to burn rate parameters,

Critical Mach number (for erosion) … cylindrical port
Only, Assume bates grain does not burn erosivley

What is the effect of Flame temperature (T0)

Plot Regression rate versus Chamber pressure

Prepare report stating your results and conclusions

6
MAE 5540 - Propulsion Systems
 Cylindrical Port: State
PDF Created with deskPDF PDF Writer - Trial :: http://www.docudesk.com

equations to program
• Use Trapezoidal rule or Runge-Kutta to integrate
 γ +1 
∂P0 Aburn aPo n  A *
 2  (γ −1) 
=  ρ p RgT0 − P0  − P0  γ RgT0 
∂t Vc Vc  γ + 1 
 
 
• Recursive propagation of chamber diameter
(k +1)∆t •
(k )∆t •
(k +1)∆t •
Rburn k +1 = Ri initial + ∫ rdt = Ri initial + ∫ rdt + ∫ rdt →
0 0 (k )∆t
(k +1)∆t • •
Rburn k +1 = Rburn k + ∫ rdt ≈ Rburn k + r ∆t = Rburn k + aPo n k ∆t
(k )∆t
7
MAE 5540 - Propulsion Systems
 8
Cylindrical Port Solution

MAE 5540 - Propulsion Systems

PDF Created with deskPDF PDF Writer - Trial :: http://www.docudesk.com
 Cylindrical Port Solution (2)
PDF Created with deskPDF PDF Writer - Trial :: http://www.docudesk.com

Chamber Pressure
Linear time
scale

Log time
scale

9
MAE 5540 - Propulsion Systems
 Cylindrical Port Solution (3)
PDF Created with deskPDF PDF Writer - Trial :: http://www.docudesk.com

Thrust/Impulse
Linear time
scale

10
MAE 5540 - Propulsion Systems
 Cylindrical Port Solution (4)
PDF Created with deskPDF PDF Writer - Trial :: http://www.docudesk.com

Massflow/regression rate
Linear time
scale

11
MAE 5540 - Propulsion Systems
 Cylindrical Port Solution (5)
PDF Created with deskPDF PDF Writer - Trial :: http://www.docudesk.com

Burn surface area,

Linear time Chamber Mach number
scale

12
MAE 5540 - Propulsion Systems
 Bates grain Port: State
PDF Created with deskPDF PDF Writer - Trial :: http://www.docudesk.com

equations to program
• Use Trapezoidal rule or Runge-Kutta to integrate
 γ +1 
∂P0 Aburn aPo n  A *
 2  (γ −1) 
=  ρ p RgT0 − P0  − P0  γ RgT0 
∂t Vc Vc  γ + 1 
 
•  
r = a ⋅ Po
n

• →
sregression = ∫ r ⋅ dt
t

( Aburn )total = N ⋅ π ⋅ 
(
 D 2 − ( d + 2 ⋅ s )2
0 0 ) 
+ ( L0 − 2 ⋅ s ) ⋅ ( d 0 + 2 ⋅ s ) 
2 
 
N ⋅π 
( ol )total = ( + ⋅ ) ⋅ ( − ⋅ ) + ⋅ ( 2 ⋅ s )
2 2
V d 2 s L 2 s D
4  
0 0 0

13
MAE 5540 - Propulsion Systems
 14
Bates Grain Solution

MAE 5540 - Propulsion Systems

PDF Created with deskPDF PDF Writer - Trial :: http://www.docudesk.com
 Bates Grain Solution (2)
PDF Created with deskPDF PDF Writer - Trial :: http://www.docudesk.com

Chamber Pressure

Linear time
scale

Log time
scale

15
MAE 5540 - Propulsion Systems
 16
Thrust/Impulse
Bates Grain Solution (3)

MAE 5540 - Propulsion Systems

PDF Created with deskPDF PDF Writer - Trial :: http://www.docudesk.com
 Bates Grain Solution (4)
PDF Created with deskPDF PDF Writer - Trial :: http://www.docudesk.com

Massflow/regression rate

17
MAE 5540 - Propulsion Systems
 Bates Grain Solution (5)
PDF Created with deskPDF PDF Writer - Trial :: http://www.docudesk.com

Massflow/regression rate
Burn surface area

18
MAE 5540 - Propulsion Systems