Innovative Solid Fuels for Hybrid Rockets
Innovative Solid Fuels for Hybrid Rockets
1
Abstract
Hybrid propulsion display several advantages over solid and liquid propulsion.
Among these, the most important are safety, thrust throttability, possibility of
engine re-start, low costs and low environmental impact. The main drawback of
hybrid technology is its fuel low regression rate. The purpose of this work is to
review the available literature on hybrid propulsion. Moreover, a comparison among
different kinds of solid fuels was performed, in terms of average regression rate. The
tested fuels are HTPB and paraffin wax. As expected from a literature review,
paraffin-based solid fuels display higher regression rate when compared to that
typical of HTPB-based fuels. This is due to the entrainment phenomenon that
characterizes liquefying fuels, which form a liquid layer on their burning surface.
2
Table of Contents
1 Motivations and Objectives………………………………………………………………................4
1.1 Motivations………………………………………………………………………………………………………..4
1.2 Objectives………………………………………………………………………………………………………….4
1.3 Outline.............................................……………………………………………………………………4
2.9 Derivation of Hybrid Fuel Regression Rate Equation for Classical Fuels…………….19
2.10 Derivation of Hybrid Fuel Regression Rate Equation for Advanced Fuels………..22
3 Experimental Setup…………………………………………………………………………................25
Nomenclature…………………………………………………………………………………….................31
References………………………………………………………………………………………….................34
3
Chapter 1
1.1 Motivations
Hybrid rocket propulsion displays several advantages when compared to solid and
liquid propulsion. The main advantages of the hybrid technology are the possibility
of thrust modulation, restart and throttling capabilities, safe transportation and
operation, low development and operational costs, and low environmental impact.
On the other hand, low regression rate is the main drawback that has limited so far
the full development of hybrid rocket technology.
The solid fuel regression rate is the rate at which the combustion front penetrates
the solid grain and travels in the direction perpendicular to the grain surface. Many
attempts have been made to increase the regression rate of the hybrid fuel,
including for example addition of energetic ingredients (metal powders, metal
hydrides), increase of convective thermal exchange coefficient through turbulence
enhancement (swirled oxidizer injection), and use of liquefying fuels and
exploitation of the entrainment phenomenon, i.e. the entrainment of droplets from
the unstable liquid-gas interface, which can substantially increase fuel mass transfer
leading to higher regression rates.
1.2 Objectives
1.3 Outline
Chapter 2 discusses the state of the art of hybrid rocket propulsion. This chapter
deals with the historical background of the hybrid rocket propulsion, the main
4
features the hybrid rocket, the combustion process and combustion instability and
the most common methods used to increase the regression rate of hybrid fuels.
In Chapter 5 the conclusion of this work are presented and future work is suggested.
5
Chapter 2
The early history of hybrid rocket development dates back to the early 1930s, at the
same time when liquid and solid rockets development began.
On the other hand, solid propellant rockets offer the advantage of simplicity, but
give lower specific impulse and problems of safety and environmental impact. Solid
propellant rocket doesn’t require the complex and expensive machinery of liquid
systems. However, in the solid propellant the solid fuel and oxidizer are mixed
together producing an explosive fuel. Upon ignition the solid fuel burns
6
uninterrupted until all the fuel is exhausted. The liquid and solid propellant rockets
are represented in figure 1.
The classical hybrid rocket configuration includes a solid fuel grain cast into the
combustion chamber and a separate tank containing a liquid oxidizer. The oxidizer
can be fed into the ports of the fuel grain in the combustion chamber either by gas
pressure or by a pump system. The inverse hybrid motor includes a liquid fuel and a
solid oxidizer. In the gas-generator type of hybrid rocket, the solid fuel is loaded
with a small amount of solid oxidizer, forming a fuel-rich solid-propellant grain.
Oxidizer is then injected into the after burner section to mix and burn with the fuel-
rich gases generated by the solid grain. The combined type has both aspects of the
classical and gas-generator types of hybrid motors because the oxidizer is injected
into both the head end of the fuel grain and the aft mixing chamber.
Classical materials used as fuels in hybrid rockets are polymers such as hydroxyl
terminated polybutadiene (HTPB). Higher regression rate fuels are obtained using
liquefying fuels, such as paraffin and polyethylene waxes, which form a liquid layer
on their burning surface. In this type of fuels the regression rate is generated by the
vaporization of the liquid and the entrainment of liquid droplets into the gas stream.
The entrainment mass transfer increases with deceasing viscosity and surface
tension of the melt layer and increasing mass flux of the oxidizer. Typical oxidizers
7
are liquid oxygen, hydrogen peroxide, nitrogen tetroxide, and nitrous oxide. In a
typical hybrid rocket motor the liquid oxidizer is injected axially into a pre-
combustion chamber where it vaporizes then it passes through the ports of the fuel
grain and reacts with fuel vapor. A post-combustion chamber ensures that all fuel
and oxidizer are burnt before exiting the nozzle. Different types of hybrid rockets
are shown in figure 3.
1. Safety: since fuel and oxidizer are separated by distance and phase, hybrids
have almost no explosion hazard and very few failure modes. Moreover, solid
fuels are not hazardous for storage and transportation, unlike solid
propellants and volatile liquid fuels such as hydrogen.
2. Insensitivity: unlike solid-propellant grains, where fuel and oxidizer are
intimately mixed, hybrid-fuel grains are insensitive to cracks and
imperfections. In conventional hybrids, heterogeneous reactions due to
oxidizer attacking the fuel surface generally do not come into play because
the diffusion flame zone shields the fuel surface from the oxidizer-rich core
flow. Therefore, potential cracks and fuel imperfections that may increase the
fuel surface area do not have a significant effect on the internal ballistics of
the fuel grain.
8
3. Reliability: since only the oxidizer is stored in liquid form, hybrid rockets
require only half the feed-system hardware needed by liquid rockets. This
provides advantages in terms of improved reliability, lower feed system
weight, and less complex mechanical design.
4. Energy management: hybrid rockets can be throttled for thrust control,
maneuvering, motor shutdown, and restart by adjusting only the oxidizer flow
rate while avoiding the necessity of matching hydraulic characteristics with
the fuel, as must be done for liquid propellant rockets.
5. Fuel versatility: the solid-phase fuel provides a convenient matrix for
introducing a variety of additives such as metal particles for high-energy
missions.
6. Design flexibility: though most hybrid development efforts have focused on
large boosters using the conventional configuration, hybrid motors can also
be used for a variety of missions including in-space propulsion, satellite
maneuvering, orbit maintenance, waste generators .
7. Environmental friendliness : hybrid rockets using typical propellants such as
liquid oxygen (LOX) and rubber-base fuel such as HTPB, have environmentally
clean exhaust without hydrogen chloride or aluminum oxide.
8. Low cost: hybrid rockets are very economical to both manufacture and launch
because of their inherent safety and minimal failure modes.
Hybrid rockets also display some disadvantages with respect to liquid and solid
rockets, as described below.
1. Slow regression rate: polymeric hybrid fuels, such as HTPB, regress rather
slowly, generally at least an order of magnitude slower than solid propellants.
To produce the necessary mass flow rate of pyrolyzed vapor from the fuel
grain consistent with a desired thrust level, multiport grains with large wetted
surface areas must be employed. Such grains require proportionally large
pressure cases and can display poor volumetric loading.
2. Low volumetric loading: in addition to slow regression rates leading to poor
volumetric loading, the use of aft combustion chambers downstream the fuel
grain to complete the mixing and combustion of fuel and oxidizer further
exacerbates the mass fraction disadvantage compared to liquid and solid
systems.
3. Fuel residuals: conventional hybrid fuels with multiple ports cannot be burnt
to completion because portions of the fuel web between the ports would
9
dislodge from the host grain and potentially block the nozzle. Therefore, a few
percent of the fuel is left intact at the end of the mission.
4. Mixture ratio shift: because of the strong coupling between the oxidizer flow
through the fuel port and the grain ballistic behavior, the overall mixture ratio
during the combustion displays a time variation due to the increase in the
port size. This variation affects the propulsion system performance. However,
use of secondary oxidizer injectors and innovative grain designs can eliminate
the mixture ratio shift.
5. Mixing/Combustion inefficiencies: hybrids have fundamentally different
mixing and combustion processes than either liquid or solid propellant
systems. In liquid engines, propellant mixing occurs on a droplet-size scale,
whereas in solids the fuel and oxidizer are intimately mixed during the grain
casting process. In hybrids, however, propellant mixing and combustion occur
in a macroscopic diffusion flame zone that has a length scale of the same
order as the fuel grain length. This model of mixing and combustion may
result in slightly lower overall combustion efficiency than competing chemical
systems.
The combustion process occurring in a hybrid engine depends on the type of solid
fuel used.
Figure 4 represents the flame structure of a hybrid rocket motor in the case of
standard fuel. The oxidizer is injected in the axial direction while the fuel pyrolisis in
the trasversal direction induces the so-called blowing effect.
The combustion takes place after ignition, when the gasified solid fuel and the
injected liquid oxidizer mix and reach approximately the stoichiometric ratio. The
flame is located within the dynamic boundary layer. This flame is diffusive and is
controlled more by fluid-dynamics rather than by chemical kinetics. The combustion
is sustained by the thermal energy feedback, by convection and radiation, from the
diffusive flame to the burning surface of the solid grain. The combustion products
are transported downstream by the convective flow and are expanded in the nozzle
for thrust generation. The temperature grows from the wall temperature to the free
stream temperature and has a maximum value in correspondence of the flame,
whereas the speed grows monotonically from wall value to the free stream value.
10
Figure 4 Hybrid propulsion flame structure.
11
Figure 5 Physical processes involved in hybrid rocket combustion.
One of the physical phenomena that limits the burning rate in a hybrid motor is the
so-called blocking effect that is caused by the high velocity injection of the
vaporizing fuel into the gas stream. This difference in the combustion scheme of a
hybrid motor significantly alters the burning rate characteristics compared to a solid
rocket. Blocking can be explained as follows: increasing the heat transfer to the fuel
causes the evaporative mass transfer from the liquid-gas interface to increase, but
the increased blowing from the surface reduces the temperature and the velocity
gradient at the surface thus reducing the convective heat transfer. The blowing also
thickens the boundary layer and displaces the flame sheet further from the fuel
surface leading to a further reduction in convective heat transfer. The position of
the flame sheet and the shape of the thermal and velocity boundary layer is the
result of a complex chemical and fluid-mechanical balance among the oxidizer flow
entering the port, the fuel flow produced by evaporation and the flow of
combustion products. As a result, the burning rate is limited in a fundamental way
which is difficult to overcome by either increasing heat transfer to the fuel or by a
reduction in the fuel heat of gasification. Although radiative heat transfer from
12
flame does not suffer from the blocking effect it is usually small compared to the
convective heat transfer.
There are a lot of methods used to increase the regression rate of solid fuels. These
methods include:
Hybrid motors have two basic types of instabilities: oxidizer feed system-induced
instability (non-acoustic), and flame holding instability (acoustic). Oxidizer feed
system instability is essentially a chugging type and arises when the feed system is
sufficiently soft. In cryogenic systems, this implies a high level of compressibility
from sources such as vapor cavities or two-phase flow in feed lines combined with
insufficient isolation from motor combustion processes. Flame-holding instabilities
arise due to inadequate flame stabilization in the boundary layer and are not
associated with feed system flow perturbations.
Axial injection of oxidizer results in a strong hot gas flow recirculation zone at the
fuel grain leading edge, producing a stable combustion. Whereas conical injection of
oxidizer can produce a weak or nonexistent hot gas flow recirculation zone at the
fuel grain leading edge, resulting in unstable combustion.
where is the inlet mass flow rate and is the outlet mass flow rate through
a supersonic gasdynamic nozzle. By making explicit all terms, one finds:
Where is the mass flow rate of the liquid oxidizer and is the mass flow rate
of the gasified gaseous fuel.
15
Where increase in time and decrease in time.
From these equations it can be concluded that the instantaneous pressure pc and
the instantaneous fuel mass flow rate increase while the instantaneous mass
mixture ratio O/F decreases in time for n<0.5; vice versa for n>0.5.
HTPB is a classical fuel used in hybrid propulsion because of its good physical
properties and good performance. It is a polymer of butadiene terminated at each
end with a hydroxyl functional group. It is a viscous liquid that become an elastic
16
solid when it reacts with another organic material (Isophorone diisocyanate, IPDI).
HTPB has the following formula:
HO-CH2-CH=CH-CH2-(CH2-CH=CH-CH2)n-CH2-CH=CH-CH2-OH
where n is the number of butadiene monomers presents in the polymer. The HTPB
used in this work is HTPB R-45, composed of 45 monomers of butadiene and has the
following properties:
These fuels such as paraffin and polyethylene waxes display a very high regression
rate with respect to the traditional fuels. They produce a very thin, low viscosity, low
surface tension liquid layer on their burning surface. The instability of this layer is
driven by the oxidizer gas flow in the port and leads to entrainment of droplets into
the gas stream greatly increasing the overall fuel mass transfer rate. This mechanism
acts like a continuous spray injection system distributed along the port. Since
droplet entrainment is not limited by diffusive heat transfer to the fuel from the
combustion zone, this mechanism can lead to much higher surface regression rates
than those typical of conventional polymeric fuels. In figure 7 the regression rate of
paraffin-based fuels is compared to that of HTPB.
Figure 7 Regression rates vs. oxidizer flux for HTPB and paraffin-based fuels.
17
The entrainment mass transfer depends on the operational parameters (pressure,
oxidizer flux) and on the material properties of the solid fuel (viscosity, surface
tension). The entrainment mechanism is represented in figure 8.
The experimental data suggested that the dynamic pressure exponent is in the
range of 1-1.5. For example Gater and L’Ecuyer scaling for large mass fluxes
indicates that α is approximately 1.5 and β is equal to 2. The viscosity and surface
tension exponents are both predicted to be 1.
18
2.9 Derivation of Hybrid Fuel Regression Rate Equation for Classical fuels
Figure 9 represents the energy balance at the fuel grain surface. The general steady-
state energy balance can be written as follows:
Energy input fuel surface = Energy output of fuel surface
Qconvection+ Qradiation in = Qconduction in + Qphase change+ Qradiation out
Neglecting radiation and in depth conduction in the fuel mass, the steady state
surface energy balance becomes:
where is the energy transferred to the fuel surface by convection, and , , and
are respectively the solid fuel density, surface regression rate, and overall fuel
heat of vaporization or decomposition. At the fuel surface the heat transferred by
convection equals that transferred by conduction so that:
evaluated at the fuel surface. The central problem in determining the hybrid fuel
regression rate is thereby reduced to determining the basic aerothermal properties
19
of the boundary layer. The heat transfer coefficient at the wall is related to the skin
friction coefficient via the following relationship:
where Cf is the skin friction coefficient with blowing, Ch is the Stanton number, and
Pr is the Prandtl number. Furthermore, the Stanton number can be written in terms
of the heat flux to the fuel surface as:
where ∆h is the enthalpy difference between the flame zone and the fuel surface,
and , are the density and velocity of oxidizer at the edge of the boundary layer.
Combining the equations, the regression rate of the fuel surface can be written as:
From boundary layer theory, one can show that the skin friction coefficient without
blowing ( ) is related to the local Reynolds number by the relation:
( )
Experiments have shown that is related to by the following relation:
( )
Figure 10 represents the blocking factor Cf/Cf0 in function of the blowing coefficient
β.
In a turbulent boundary layer, the prandtl number is very nearly equal to 1. It can be
shown that for Pr=1, β as defined in the previous equation, is also equal to ∆h/hv.
Noting that is the definition of oxidizer mass velocity (G), the equation that
express the regression rate can be written in the final form as:
20
The coefficient 0.036 is applied when the quantities are expressed in the English
Engineering system of units.
In some hybrid motors, radiation may be a significant contributor to the total fuel
surface heat flux. Such motors include those with metal additives to the fuel grain
(such as aluminum) or motors in which soot may be present in significant
concentrations in the combustion chamber. In these instances the regression rate
equation must be modified to account for heat flux from a radiating particle cloud.
The radiative contribution affects surface blowing, and hence the convective heat
flux as well. In this case the total heat flux to the fuel surface is expressed by:
The radiation heat flux has been hypothesized to have the following form:
Where and are respectively the regression rate with and without radiation.
21
2.10 Derivation of Hybrid Fuel Regression Rate Equation for Advanced Fuels
The formation of liquid layer instabilities and entrainment of liquid droplets require
three major modifications in the classical hybrid combustion theory:
1. The ratio of the enthalpy difference to the effective heat of gasification (∆h/hv)
that appears in the thermal blowing parameter expression is altered. The
effective heat of gasification is reduced because the evaporation energy required
for the fuel mass transfer from the surface is partly avoided by the mechanical
entrainment of the liquid, whereas the enthalpy difference between the flame
and the surface is also reduced because some of the reactants are now in liquid
phase.
2. The blocking factor CH/CH0 that modifies the convective heat flux to the surface is
also altered as a result of the presence of the two-phase flow. The blocking factor
can be expressed as a function of evaporation blowing parameter:
The evaporation blowing parameter Bg includes only the gaseous phase mass
transfer from the fuel surface.
3. The ripples formed on the liquid layer surface increase the surface roughness and
the heat transfer from the flame front to the surface.
In general, the total regression rate of a hybrid motor can be written as a sum of the
evaporation regression rate that is generated by the vaporization of the liquid into
the gas stream and the entrainment regression rate that is related to the mass
transfer mechanically extracted from the liquid surface.
For an arbitrary combination of the entrainment and evaporative mass transfer, the
energy balance at liquid gas interface is:
where
22
The nondimensional energy parameters for entrainment Rhe and vaporization Rhv are
introduced because the material that is extracted through the entrainment
mechanism possesses different heating histories (that is, no heat of vaporization is
required for entrainment). The assumption is made that the effective heating in the
liquid phase required for fuel material, which is going through the entrainment mass
transfer mechanism, reduces linearly as the vaporization component of the
regression rate decreases.
The roughness parameter Fr is introduced in the energy equation to account for the
increased heat transfer by wrinkling of the liquid surface. It has been argued by
Gater and L’Ecuyer that the surface roughness decreases with increasing dynamic
pressure of the gas flow. The empirical formula for the roughness correction
parameter suggested by Gater and L’Ecuyer can be expressed in terms of the
operational parameters of the motor as:
23
Figure 11 Blowing correction according to Marxman’s formula.
The entrainment regression rate can be written as in terms of the mass flux in the
port and the total regression rate as:
where aent is a function of the properties of the selected fuel and average gas
density. For simplicity the assumption is made that this coefficient is constant for a
given fuel.
By resolving this set of non linear algebraic equations the total regression rate as a
function of the axial location and local mass flux can be obtained.
24
Chapter 3
Experimental Setup
The experimental set up used for this work includes a combustion chamber, an
oxidizer inlet system, a nitrogen inlet system and a pyrotechnic ignition device.
Nitrogen is used as a cooling flow after combustion shutdown, with the aim to
extinguish the oxidation reactions in the combustion chamber.
25
Figure 13 Single slab grain. Sample sizes: 50 x10 x5 mm.
The basic principle used to measure the flow rate is that when a fluid stream is
restricted, its pressure decreases by an amount which depends on the rate of flow
through the restriction. Therefore, the pressure difference between points before
and after the restriction can be used to indicate flow rate. Several devices are used
to measure the flow rate, the most common devices are the venturi tube, the flow
nozzle, the orifice, and the flow tube. The derivation of the relationship between the
pressure difference and the volume flow rate is the same regardless of which type of
device is used. An orifice plate is used to measure the flow rate.
An orifice plate is a restriction with an opening smaller than the pipe diameter which
is inserted in the pipe; the typical orifice plat has a concentric, sharp edged opening,
as shown in Figure 14. Because of the smaller area the fluid velocity increases,
causing a corresponding decrease in pressure. The flow rate can be calculated from
the measured pressure drop across the orifice plate. The orifice plate is the most
commonly used flow sensor, but it creates a rather large non-recoverable pressure
due to the turbulence around the plate, leading to high energy consumption .The
pressure difference before and after the orifice plate is measured by a differential
manometer.
26
The Reynolds number is determined by this formula:
Where V is the velocity of the flow, D is the diameter of the pipe, and μ is the fluid
dynamic viscosity.
The oxidizer mass flow rate can be determined multiplying the volume flow rate by
the density of the fluid ( ).
The average regression rate is determined as the ratio between the mass of the fuel
consumed during combustion and the product of burning time, fuel density and fuel
combustion area:
,
where ∆m is the burned fuel mass, ߩf is the density of the solid fuel, tb is the burning
time and Ab is the burning surface.
27
Chapter 4
Firing tests were performed in double slab fuel configuration, using pure oxygen as
oxidizer, with operating pressure of 1.5 bar and oxygen mass flux ranging from 150
kg/m2s to 350 kg/m2s. Fuels tested include pure HTPB and paraffin.
The solid fuel regression rate can be expressed as a function of the oxidizer mass
flux as:
where:
The constants a and n are obtained fitting the experimental data with a power law.
The values of these constants obtained for pure HTPB fuel are a=0.05 and n= 0.44.
For paraffin fuel, the coefficients a and n obtained are higher (a=0.061 and
n=0.5437) than those obtained for HTPB. The relationship between the regression
rate and the oxidizer mass flux for the paraffin wax is :
Figure 15 shows a comparison between the regression rate of HTPB and paraffin
wax in function of oxygen mass flux.
28
Figure 15 Regression rate of HTPB and paraffin wax in function of oxygen mass flux.
It can be observed that paraffin regression rate is higher than that of HTPB (+100%
at 200 kg/m2s). This is because in the case of paraffin wax the regression rate is the
sum of the entrainment regression rate and the evaporation regression rate, while
in the case of HTPB the regression rate depends only on evaporation. The
entrainment of droplets by the gaseous oxidizer stream flowing in the port above
the surface is the main reason for this increase of the regression rate.
29
Chapter 5
In this work a survey of the available literature on hybrid propulsion was performed,
focusing on historical background, main features of a hybrid rocket engine, methods
for regression rate increase, combustion instability. Attention was also focused on
the differences between standard and innovative fuels, and their respective
combustion processes. Some firing tests were then performed in order to compare
HTPB and paraffin fuels in terms of average regression rate. Paraffin fuel use results
in a regression rate increase of 100% with respect to HTPB, at the selected oxygen
mass flux. The higher regression rate values displayed by paraffin are due to its
higher tendency to entrainment phenomenon.
30
Nomenclature
= regression rate, mm/s
B = blowing parameter
Bg = evaporation blowing parameter
Fr = heat-transfer correction factor for surface roughness
h = melt layer thickness, mm
G = mass flux of oxidizer, kg/m2-s
z = axial distance along the port, m
σ = surface tension, milliN/m
Pd = dynamic pressure in the port, Pa
tb = burning time, s
μ = viscosity, milliPa-s
CH = Stanton number with blowing
CH0 = Stanton number without blowing
CB1, CB2 = blowing correction coefficients
Cf, Cf0 = skin friction coefficients with and without blowing
∆h = enthalpy difference between the flame and the surfaces, kJ
k, ratio of specific heats (cp/cv)
R = the gas constant [J/Kg K]
Pr = Prandtl number
Re = Reynlods number
hm, he = total heats of melting and entrainment, kJ/Kg
31
Lm, Lv = latent heat of melting and vaporization, kJ/kg
Tg = average gas phase temperature, K
Tm, Tv = melting and vaporization temperature, K
∆T1 = temperature difference, Tv-Tm, K
Rh = ratio of effective heats of gasification for entrainment and vaporization
C = specific heat kJ/kg-s
c* = characteristic velocity, m/s
η c* = efficiency of the characteristic velocity
Cf = skin-friction coefficient
α,β = dynamic pressure and thickness exponents
At = nozzle throat area, m2
Pc = combustion chamber pressure, pa
Tc = combustion chamber temperature, K
Mc = average molar mass of combustion products, kg/kmol
= oxidizer mass flow rate, kg/s
Subscripts
ent = entrainment
g = gas
32
l = liquid
s= solid
v = vaporization
f = fuel
33
References
[1] G.P. Sutton, Rocket Propulsion Elements, Wiley, New York, USA, 6th edition,
1992.
34