# Approaches to Low Fuel Regression Rate in Hybrid Rocket Engines

Authors: Dario Pastrone
Journal: International Journal of Aerospace Engineering
 (2012)

---

## Abstract

Hybrid rocket engines are promising propulsion systems which present appealing features such as safety, low cost, and environmental friendliness. On the other hand, certain issues hamper the development hoped for. The present paper discusses approaches addressing improvements to one of the most important among these issues: low fuel regression rate. To highlight the consequence of such an issue and to better understand the concepts proposed, fundamentals are summarized. Two approaches are presented (multiport grain and high mixture ratio) which aim at reducing negative effects without enhancing regression rate. Furthermore, fuel material changes and nonconventional geometries of grain and/or injector are presented as methods to increase fuel regression rate. Although most of these approaches are still at the laboratory or concept scale, many of them are promising.


---

## Body

## 1. Introduction

Hybrid rocket engines (HREs) are chemical rockets which present interesting advantages over liquid rocket engines (LREs) and solid rocket motors (SRMs) and can provide a safe and affordable option for many applications. They have recently come to the fore as they have been chosen to power the second stage of the Tier One, the winner of the $10 million  X-Prize [1]. The Tier One (http://www.scaled.com/projects/tierone/) was the first privately-developed reusable vehicle able to perform a manned suborbital flight over 100 km. The first stage of the Tier One is a twin-turbojet carrier-launch aircraft, called the White Knight, which is able to bring the second stage to 15 km altitude. The second stage, the SpaceShipOne, is a reusable three-place manned space plane, powered by an N2O/HTPB hybrid rocket engine. This HRE has a burning time of about 80 seconds and it is able to produce an average thrust of 75 kN with a vacuum specific impulse of about 250 s. The Tier One was retired after winning the X-Prize, but it remains proof that HREs may be the key propulsion system of the emerging space market [2].In HREs, oxidizer and fuel are separated and stored in two different physical phases. In the most common configuration, a liquid or gaseous oxidizer and a solid-fuel grain are employed (direct HRE). In the conventional configuration the grain is cylindrical with a circular port. Due to this propellant storage, HREs have a peculiar combustion process as shown Figure 1. A boundary layer is formed when the oxidizer is injected into the port of the fuel grain. After ignition a diffusion flame is formed inside this layer. The fuel deriving from the solid grain is gasified by the heat coming from the flame, which in turn is fed by the fuel. The gasified fuel blows from the surface and modifies the boundary layer, blocking the heat transfer.Figure 1
HRE combustion.The aforementioned propellant storage and combustion process determine both positive and negative effects. Many authors have discussed the advantages and the drawbacks of HREs. Characteristic hybrid rocket features are briefly summarized below. Benefits include the following.(i)
Performance: HREs may deliver a higher specific impulse than SRMs. Due to the high-density solid they also may have higher density specific impulse than LREs. This latter benefit may be wasted by sliver or low grain volumetric efficiency.(ii)
Safety: HREs are inherently safe and low cost mainly due to the use of a solid fuel grain which is classically inert. Since the fuel and oxidizer are separated by distance and phase, hybrids have almost no explosion hazard and very few failure modes. The fluid oxidizer is usually depleted in the diffusion flame region. It can reach and attack the fuel surface only when gaseous chemical kinetics are slow (low pressures). No deflagration to detonation transition can be caused by pores, cracks, and imperfections, and inadvertent ignitions are avoided. Last, but not least, defining a maximum expected operating pressure is not a concern as HREs do not have the high temperature sensitivity which characterizes SRMs.(iii)
Reliability and simplicity: feeding system hardware is reduced as only the oxidizer is liquid.(iv)
Fuel versatility: additives for different purposes can be embedded in the fuel grain.(v)
Oxidizer control: liquid propellant control allows for throttling, motor shutdown, and reignition.(vi)
Environmental friendliness: compared to SRMs, oxidizers are chlorine free. Compared to LREs, storable propellants are available, which do not have noxious features such as MMH and nitrogen tetroxide.(vii)
Low cost: hybrid rockets pose almost no explosion hazard during manufacture, transport, ground test, and storage. Low recurring costs are foreseen because of high levels of safety and minimal failure modes.Classical HREs also present disadvantages which include the following.(i)
Performance: HREs cannot reach the high specific impulse of cryogenic bipropellant LREs and have a lower density specific impulse compared to SRMs.(ii)
Low fuel regression rate: the regression rate of conventional binders such as HTPB is typically an order of magnitude lower than solid propellants and hence a large fuel surface is needed to produce the required thrust level, as discussed later in Section2. This limit is usually set by the physical phenomenon of heat transfer from the diffusion flame to the fuel surface. As a consequence, HREs may have poor fuel loading, low thrust densities, and large length-to-diameter ratios.(iii)
Mixture ratio shifting: the regression rate primarily depends on the mass-flux, that is, on the ratio of mass-flow and port areaAp. During combustion of classical side burning grains, Ap increases determining a regression rate reduction. On the other hand, the burning area Ab increases as well. These two competing effects are, in general, not balanced, leading to a change in the fuel mass flow, even when the oxidizer mass-flow is kept constant (see Section 2). The mixture ratio shifts and performance are worsened. Nonconventional injection systems may fix this issue but introduce complexity.(iv)
Mixing inefficiencies: part of the fuel under the flame at the grain port exit may not mix with any oxidizer and thus exit the nozzle before releasing chemical energy. HREs have a lower overall combustion efficiency than LREs and SRMs. A mixer between the grain aft end and the nozzle inlet can enhance mixing, but the dry weight of the system is larger.(v)
Slow transient/response to throttling: due to the thermal lag in the solid fuel, ignition and response to throttling is slow in comparison to LREs. Also, the chamber volume may be larger than in LREs, with large tail-off time.One of the most important issues remains the very low regression rate of the fuel grain. Various methods for enhancing regression rate or, at least, for reducing negative effects of low fuel regression have been suggested in the past. The goal of this paper is to survey broadly some of the more significant approaches proposed. After summarizing some fundamentals (Section2), two approaches are presented (Section 3) which aim at reducing the negative effects of regression rate without enhancing regression rate itself. Approaches which directly try to enhance regression rate are then presented: the modification of fuel is considered in Section 4, while approaches based on unconventional geometries for the grain and/or the injector are presented in Section 5.


## 2. Fundamentals

In this section the fundamentals are summarized in order to better understand the negative effects of low regression rate and in turn to allow one to consider approaches to mitigate this issue. An exhaustive treatment of this subject is beyond the scope of this work. Rudiments of regression rates behavior and models are first discussed. A ballistic model is then used to highlight that cylindrical single-port side-burning grain may have an unacceptable shape if large thrust levels are required. For further information consult [14].


### 2.1. Regression Rate Behavior and Modeling

Typical fuel regression rate behavior is depicted in Figure2 as a function of mass flux G. Three different regions can be identified. For the medium G range the regression rate is diffusion dominated and is a function of the mass flux [15–17]. When G assumes high values, the combustion appears to be controlled by chemical kinetics and not  by diffusion  [18]. As a consequence, regression rate becomes pressure dependent. An upper bound for the mass flux exists, called the flooding limit, which depends on pressure level and propellant combination. When approaching this limit, the mass flux may become so high that the flame is extinguished due to small Damköhler numbers and/or very oxidizer-rich conditions. Finally, at low values of mass flux, the convective heat transfer diminishes and radiation from gaseous species may play a significant role. Consequently, the regression rate is enhanced depending on the partial pressure of emitting components and port diameter. In this regime a lower bound of the mass flux is also present. The regression rate is small, the fuel remains in the thermal layer of the grain a long time, and the solid fuel may be cooked/melt underneath the grain surface. Therefore, if the mass flux is too low (e.g., at the end of HRE operation with blowdown feed systems), chuffing instability may be produced by the perpetual repeating mechanical removal of this soft layer.Figure 2
Typical regression rate behavior.The well-known theory developed by Marxman and colleagues [15–17] set the basis for understanding and describing fuel regression rate. They assumed that the fuel regression rate is controlled by heat transfer to the grain and applied a flux balance at the fuel surface, finding (1)rρF=qΔHv,eff,
where r is the fuel regression rate, ρF is the fuel density, q is the total heat flux received by the fuel grain, and ΔHv,eff is the thermal energy required to change a unit mass of solid fuel into gas.According to their studies, the regression rate is primarily governed by convective heat transfer, that is, by the local mass fluxG. Modeling the heat transfer in a turbulent layer, they proposed a simplified regression rate expression for combustion with negligible radiation.   Taking into account a refit of the data proposed by Altman  Humble  [19] this expression is (2)rρF∝B0.32G0.8x-0.2,5<B<100.
The blowing parameter B is the ratio of the thermal energy of the main stream relative to the surface ΔHfw and ΔHv,eff.  The weak negative dependence on axial position x reflects the effects of boundary layer growth on heat transfer. On the contrary, the total local mass flux increases with axial position along the fuel grain. These two competing effects usually determine a location of minimum regression rate along the grain axis. In any case, these difference are usually small. For these reasons the regression rate is assumed to be constant along the grain axis and semiempirical correlations based on the mass flow entering the port area, that is, GO=ṁO/Ap are widely used, that is, (3)r=aGOn,
where a and n are obtained by experimental data for a given GO range and depend on engine dimensions (scale effect), port/grain geometry, injector geometry, and flow features. Correlation like the one presented in (3) well-describe the regression rate of conventional HREs in the intermediate range of GO, where behavior is dominated by turbulent heat transfer.Radiation from gas-phase products is usually small compared to convection. Nevertheless, convection itself becomes small at low values ofG and radiation plays a role. Radiation effects are also important when the combustion products contain condensed matter, for example, soot and metal/oxides particles. Marxman and coworkers [16, 17] proposed a correction factor to take into account the radiation heat transfer. The total heat flux q in (1) becomes q=qc[(qr/qc)+exp(-qr/qc)]. This expression takes into account that there is a coupling effect between radiant flux qr and convective heat qc. In fact, the radiant flux enhances the blocking effect due to blowing thus depressing convective heat.Other phenomena that affect fuel regression rate have been considered in later works [20–28]. Beside radiation, the most important mechanisms taken into account include variable fluid and transport properties across the boundary layer, and both gas-phase and heterogenous reaction kinetics.The regression rate correlation of (3) should be modified consequently. As an example, Chiaverini et al. [28] proposed a modified relation for cases where kinetics effects can be disregarded. Variable fluid and transport properties into the boundary layer and radiation from soot are considered. The resulting model is (4)ρFrG=0.0155(DhLg)0.3θ0.6ReD-0.2f[B,qrqc],
where Dh is the port hydraulic diameter, Lg the grain length, and θ a temperature factor, representing the ratio of average flame to surface temperature [29].Although correlations like (4) allow for a better fit of experimental data, (3) is often used. The reader should be aware that the values of a and n may be dependent on many factors such as GO range, fuel formulation, grain production method, engine dimensions (scale effect), port/grain geometry, injector geometry, and flow features. For example, if (3) is used to fit experimental data, r may be found to vary with GO to the power of 0.6 instead of the theoretical 0.8 for turbulent flow over a flat plate used in (2). These findings imply that nonconvective processes, such as radiant heat flux and/or finite rate chemical kinetics, may influence the solid-fuel regression rates in this case.


### 2.2. Low Fuel Regression Rate Effects in Classical HREs

When the propellant combination is given, the characteristic velocityc* is a function of the mixture ratio α=ṁO/ṁF (i.e., oxidizer to fuel mass flow ratio) [3]. Chamber pressure affects c* to a lesser extent. The thrust coefficient CF can be evaluated when the following data are provided: combustion gas properties (specific heat ratio), ambient pressure pa, and expansion ratio ϵ=Ae/At (exit to throat area ratio) or pe/pc (exit to chamber pressure ratio). The corresponding effective exhaust velocity c=c*CF or the specific impulse Is=c/g0 is consequently evaluated.If a thrust level is assigned, the needed propellant mass flowmp is (5)ṁp=(1+α)ṁF=Fc*CF,
and, using (3), the fuel mass flow mF is (6)ṁF=aρF(ṁOAp)nAb,
where ρF is the solid fuel density.The mixture ratio is(7)α=ṁOṁF=ṁO(1-n)aρFApnAb.
This relation shows that α changes in proportion to Ab/Apn. In most cases Ab/Apn is a function of time and α changes even if ṁO is kept constant. This characteristic behavior of HREs is called mixture-ratio shifting. It is worthwhile to observe that mixture ratio shifting does not occur if a circular port is used and n=0.5.From (5) and (6): (8)Ab=1aρF(Fc*CF)(1-n)Apn(1+α)(1-n)αn.
Equation (8) shows the relation of geometry (grain/nozzle) engine operation (mixture ratio) and mission requirements (thrust). Due to the low regression rate, the length-to-diameter ratio may be of concern. In fact, large values of Ab are required, while Ap must keep GO in the appropriate range. Assuming a single-circular port, for given values of F, α, and c, the length-to-diameter ratio of the grain becomes (Lg/Dp)∝Dp(2n-1), where Lg is the grain length, and Dp the port diameter. The outer diameter of the grain is Dg=Dp+2w, where w is the web thickness. If n>0.5 the length-to-diameter ratio of the grain may be reduced with a port diameter reduction. On the other hand, the overall length-to-diameter ratio L/D must be considered. The rocket length L includes the oxidizer tank length. The reduction of the port area makes the ratio L/D worse if the oxidizer tank and grain have the same diameter [30]. On the contrary, if the diameter of the oxidizer tank may be larger that grain diameter, both L is reduced and D augmented thus improving the ratio L/D. Therefore, the reduction of Ap may be favorable [31]. However, this approach is limited as Ap has a lower bound due to flooding.As a consequence, HREs are generally restrained to high length-to-diameter ratios and thus are long and skinny. The values ofL/D become unacceptable when requiring high performance (e.g., upper stages and boosters), especially when compared to LREs and SRMs that are currently on the market. Solutions to this issue become necessary.


## 2.1. Regression Rate Behavior and Modeling

Typical fuel regression rate behavior is depicted in Figure2 as a function of mass flux G. Three different regions can be identified. For the medium G range the regression rate is diffusion dominated and is a function of the mass flux [15–17]. When G assumes high values, the combustion appears to be controlled by chemical kinetics and not  by diffusion  [18]. As a consequence, regression rate becomes pressure dependent. An upper bound for the mass flux exists, called the flooding limit, which depends on pressure level and propellant combination. When approaching this limit, the mass flux may become so high that the flame is extinguished due to small Damköhler numbers and/or very oxidizer-rich conditions. Finally, at low values of mass flux, the convective heat transfer diminishes and radiation from gaseous species may play a significant role. Consequently, the regression rate is enhanced depending on the partial pressure of emitting components and port diameter. In this regime a lower bound of the mass flux is also present. The regression rate is small, the fuel remains in the thermal layer of the grain a long time, and the solid fuel may be cooked/melt underneath the grain surface. Therefore, if the mass flux is too low (e.g., at the end of HRE operation with blowdown feed systems), chuffing instability may be produced by the perpetual repeating mechanical removal of this soft layer.Figure 2
Typical regression rate behavior.The well-known theory developed by Marxman and colleagues [15–17] set the basis for understanding and describing fuel regression rate. They assumed that the fuel regression rate is controlled by heat transfer to the grain and applied a flux balance at the fuel surface, finding (1)rρF=qΔHv,eff,
where r is the fuel regression rate, ρF is the fuel density, q is the total heat flux received by the fuel grain, and ΔHv,eff is the thermal energy required to change a unit mass of solid fuel into gas.According to their studies, the regression rate is primarily governed by convective heat transfer, that is, by the local mass fluxG. Modeling the heat transfer in a turbulent layer, they proposed a simplified regression rate expression for combustion with negligible radiation.   Taking into account a refit of the data proposed by Altman  Humble  [19] this expression is (2)rρF∝B0.32G0.8x-0.2,5<B<100.
The blowing parameter B is the ratio of the thermal energy of the main stream relative to the surface ΔHfw and ΔHv,eff.  The weak negative dependence on axial position x reflects the effects of boundary layer growth on heat transfer. On the contrary, the total local mass flux increases with axial position along the fuel grain. These two competing effects usually determine a location of minimum regression rate along the grain axis. In any case, these difference are usually small. For these reasons the regression rate is assumed to be constant along the grain axis and semiempirical correlations based on the mass flow entering the port area, that is, GO=ṁO/Ap are widely used, that is, (3)r=aGOn,
where a and n are obtained by experimental data for a given GO range and depend on engine dimensions (scale effect), port/grain geometry, injector geometry, and flow features. Correlation like the one presented in (3) well-describe the regression rate of conventional HREs in the intermediate range of GO, where behavior is dominated by turbulent heat transfer.Radiation from gas-phase products is usually small compared to convection. Nevertheless, convection itself becomes small at low values ofG and radiation plays a role. Radiation effects are also important when the combustion products contain condensed matter, for example, soot and metal/oxides particles. Marxman and coworkers [16, 17] proposed a correction factor to take into account the radiation heat transfer. The total heat flux q in (1) becomes q=qc[(qr/qc)+exp(-qr/qc)]. This expression takes into account that there is a coupling effect between radiant flux qr and convective heat qc. In fact, the radiant flux enhances the blocking effect due to blowing thus depressing convective heat.Other phenomena that affect fuel regression rate have been considered in later works [20–28]. Beside radiation, the most important mechanisms taken into account include variable fluid and transport properties across the boundary layer, and both gas-phase and heterogenous reaction kinetics.The regression rate correlation of (3) should be modified consequently. As an example, Chiaverini et al. [28] proposed a modified relation for cases where kinetics effects can be disregarded. Variable fluid and transport properties into the boundary layer and radiation from soot are considered. The resulting model is (4)ρFrG=0.0155(DhLg)0.3θ0.6ReD-0.2f[B,qrqc],
where Dh is the port hydraulic diameter, Lg the grain length, and θ a temperature factor, representing the ratio of average flame to surface temperature [29].Although correlations like (4) allow for a better fit of experimental data, (3) is often used. The reader should be aware that the values of a and n may be dependent on many factors such as GO range, fuel formulation, grain production method, engine dimensions (scale effect), port/grain geometry, injector geometry, and flow features. For example, if (3) is used to fit experimental data, r may be found to vary with GO to the power of 0.6 instead of the theoretical 0.8 for turbulent flow over a flat plate used in (2). These findings imply that nonconvective processes, such as radiant heat flux and/or finite rate chemical kinetics, may influence the solid-fuel regression rates in this case.


## 2.2. Low Fuel Regression Rate Effects in Classical HREs

When the propellant combination is given, the characteristic velocityc* is a function of the mixture ratio α=ṁO/ṁF (i.e., oxidizer to fuel mass flow ratio) [3]. Chamber pressure affects c* to a lesser extent. The thrust coefficient CF can be evaluated when the following data are provided: combustion gas properties (specific heat ratio), ambient pressure pa, and expansion ratio ϵ=Ae/At (exit to throat area ratio) or pe/pc (exit to chamber pressure ratio). The corresponding effective exhaust velocity c=c*CF or the specific impulse Is=c/g0 is consequently evaluated.If a thrust level is assigned, the needed propellant mass flowmp is (5)ṁp=(1+α)ṁF=Fc*CF,
and, using (3), the fuel mass flow mF is (6)ṁF=aρF(ṁOAp)nAb,
where ρF is the solid fuel density.The mixture ratio is(7)α=ṁOṁF=ṁO(1-n)aρFApnAb.
This relation shows that α changes in proportion to Ab/Apn. In most cases Ab/Apn is a function of time and α changes even if ṁO is kept constant. This characteristic behavior of HREs is called mixture-ratio shifting. It is worthwhile to observe that mixture ratio shifting does not occur if a circular port is used and n=0.5.From (5) and (6): (8)Ab=1aρF(Fc*CF)(1-n)Apn(1+α)(1-n)αn.
Equation (8) shows the relation of geometry (grain/nozzle) engine operation (mixture ratio) and mission requirements (thrust). Due to the low regression rate, the length-to-diameter ratio may be of concern. In fact, large values of Ab are required, while Ap must keep GO in the appropriate range. Assuming a single-circular port, for given values of F, α, and c, the length-to-diameter ratio of the grain becomes (Lg/Dp)∝Dp(2n-1), where Lg is the grain length, and Dp the port diameter. The outer diameter of the grain is Dg=Dp+2w, where w is the web thickness. If n>0.5 the length-to-diameter ratio of the grain may be reduced with a port diameter reduction. On the other hand, the overall length-to-diameter ratio L/D must be considered. The rocket length L includes the oxidizer tank length. The reduction of the port area makes the ratio L/D worse if the oxidizer tank and grain have the same diameter [30]. On the contrary, if the diameter of the oxidizer tank may be larger that grain diameter, both L is reduced and D augmented thus improving the ratio L/D. Therefore, the reduction of Ap may be favorable [31]. However, this approach is limited as Ap has a lower bound due to flooding.As a consequence, HREs are generally restrained to high length-to-diameter ratios and thus are long and skinny. The values ofL/D become unacceptable when requiring high performance (e.g., upper stages and boosters), especially when compared to LREs and SRMs that are currently on the market. Solutions to this issue become necessary.


## 3. Mitigation Approaches

In this section two solutions which attempt to mitigate the negative effects of low regression rate are presented. Multiport grains allow for reducingL/D values with an almost constant fuel mass, while high mixture ratio aims at reducing the mass of the fuel grain.


### 3.1. Multiport Grains

Assuming that bothAb and Ap are constant, the grain length Lg=Ab/P may be reduced by increasing the perimeter P. This result can be obtained using a multiport grain. The larger the number of ports, the more pronounced the length reduction. Also, volumetric efficiency is increased as the port volume is proportional to Lg. As an example, two possible multiport geometries are shown in Figure 5. In Figure 5(a) a triangular port [32] is shown. Sliver is a drawback which can be reduced using a wagon-wheel grain (Figure 5(b)) where a central circular port with initial radius Rci is surrounded by a row of N quadrangular ports [33]. Grains with two or three rows of quadrangular ports have been tested [30]. Even with these more complicated geometries the unburned mass fraction is high. Furthermore, other drawbacks exist. The individual ports may behave in a different way and dedicated injector or a large prechamber would be required, compromising weight and simplicity. In addition, design and fabrication are complex and grain structural integrity may become an issue, especially towards the end of the burn.Radial-flow/end-burning grains are other interesting geometry alternatives which may improve fuel section volumetric efficiency. Combination of end-burning with side-burning grains is possible. In this case the flow structure is modified and the regression rate may be enhanced too. Approaches which involve end burning grain or radial flows are left for the next sections which discuss regression rate enhancement.Also more complicated geometries can be used which allow for a better use of volume. Rapid prototyping has been identified as a manufacturing technique for hybrid rocket fuel grains that enables the production of complex 3D grain shapes, which may improve hybrid rocket motor performance [34].


### 3.2. High Mixture-Ratio Values

If high mixture-ratio values are used, the regression rate is of minor concern, as the fuel grain must give a minor contribution to the overall mass flow. Unfortunately, the propellant performance is a function of mixture ratio as well. Benefits may be obtained by adopting propellant combinations which present high values ofc* at high mixture rations. Figure 3 shows c* for three different propellant combinations: HTPB with liquid oxygen (LOX) or N2O, and polyethylene (PE) with an 85% concentration solution of hydrogen peroxide (HP). The thrust coefficient CF is primarily a function of the expansion ratio and is only slightly affected by α, reaching a maximum value in the stoichiometric region. As a result the maximum value of c is somehow shifted to the right of α which determines the maximum of c*.Figure 3
Theoretical values of characteristic velocity for different propellant combinationsm [3].It is worthwhile to note here that proper tools are needed to perform a comparison [30, 31, 35–37]. In many applications thrust has a major influence since it affects both the propulsion system design and the trajectory performance; a compromise must be sought, as greater thrust levels reduce the gravitational losses but increase structural mass. Moreover, HREs are characterized by their peculiar combustion process and the consequent link between thrust level and mixture ratio. In order to efficiently perform a coupling of the HRE design parameter optimization with the trajectory optimization, a unique procedure has been developed at the Politecnico di Torino which uses a direct/indirect nested method.As an example, the performance of a microgravity platform for three different propellant combination are compared [31]. Given rocket payload (100 kg) and initial mass (500 kg, comprising payload, fixed masses, propulsion system, and propellant), the time spent above 100 km (tμg) is the performance index to be maximized. The optimization aims at finding the optimal mass split between propellant and propulsion system (i.e., tanks, combustion chamber, nozzle), the optimal grain geometry, and the corresponding optimal trajectory [31]. Single-port grains are considered and different propellant options are compared in Table 2. A simple blowdown feed system is adopted.Results show the superior performance of the HP/PE combination. Propellant consumption and propulsion system dry-mass are competing parameters. The LOX/HTPB propellant combination presents higher maximum values ofc*, but this maximum is reached for small values of α: large grains are required and the propulsion system mass penalty is higher than propellant savings. Also, LOX and HP have larger densities than HTPB, and higher values of density-specific impulse may be obtained with the HP/PE combination which adopts high values of α. On the other hand, N2O is penalized by the low performance in terms of specific impulse. However, N2O presents some interesting features, such as the reduced rocket length L and diameter D, as the large mixture ratio reduces the grain length. Also note that the initial required thrust for the same mission is lower.


## 3.1. Multiport Grains

Assuming that bothAb and Ap are constant, the grain length Lg=Ab/P may be reduced by increasing the perimeter P. This result can be obtained using a multiport grain. The larger the number of ports, the more pronounced the length reduction. Also, volumetric efficiency is increased as the port volume is proportional to Lg. As an example, two possible multiport geometries are shown in Figure 5. In Figure 5(a) a triangular port [32] is shown. Sliver is a drawback which can be reduced using a wagon-wheel grain (Figure 5(b)) where a central circular port with initial radius Rci is surrounded by a row of N quadrangular ports [33]. Grains with two or three rows of quadrangular ports have been tested [30]. Even with these more complicated geometries the unburned mass fraction is high. Furthermore, other drawbacks exist. The individual ports may behave in a different way and dedicated injector or a large prechamber would be required, compromising weight and simplicity. In addition, design and fabrication are complex and grain structural integrity may become an issue, especially towards the end of the burn.Radial-flow/end-burning grains are other interesting geometry alternatives which may improve fuel section volumetric efficiency. Combination of end-burning with side-burning grains is possible. In this case the flow structure is modified and the regression rate may be enhanced too. Approaches which involve end burning grain or radial flows are left for the next sections which discuss regression rate enhancement.Also more complicated geometries can be used which allow for a better use of volume. Rapid prototyping has been identified as a manufacturing technique for hybrid rocket fuel grains that enables the production of complex 3D grain shapes, which may improve hybrid rocket motor performance [34].


## 3.2. High Mixture-Ratio Values

If high mixture-ratio values are used, the regression rate is of minor concern, as the fuel grain must give a minor contribution to the overall mass flow. Unfortunately, the propellant performance is a function of mixture ratio as well. Benefits may be obtained by adopting propellant combinations which present high values ofc* at high mixture rations. Figure 3 shows c* for three different propellant combinations: HTPB with liquid oxygen (LOX) or N2O, and polyethylene (PE) with an 85% concentration solution of hydrogen peroxide (HP). The thrust coefficient CF is primarily a function of the expansion ratio and is only slightly affected by α, reaching a maximum value in the stoichiometric region. As a result the maximum value of c is somehow shifted to the right of α which determines the maximum of c*.Figure 3
Theoretical values of characteristic velocity for different propellant combinationsm [3].It is worthwhile to note here that proper tools are needed to perform a comparison [30, 31, 35–37]. In many applications thrust has a major influence since it affects both the propulsion system design and the trajectory performance; a compromise must be sought, as greater thrust levels reduce the gravitational losses but increase structural mass. Moreover, HREs are characterized by their peculiar combustion process and the consequent link between thrust level and mixture ratio. In order to efficiently perform a coupling of the HRE design parameter optimization with the trajectory optimization, a unique procedure has been developed at the Politecnico di Torino which uses a direct/indirect nested method.As an example, the performance of a microgravity platform for three different propellant combination are compared [31]. Given rocket payload (100 kg) and initial mass (500 kg, comprising payload, fixed masses, propulsion system, and propellant), the time spent above 100 km (tμg) is the performance index to be maximized. The optimization aims at finding the optimal mass split between propellant and propulsion system (i.e., tanks, combustion chamber, nozzle), the optimal grain geometry, and the corresponding optimal trajectory [31]. Single-port grains are considered and different propellant options are compared in Table 2. A simple blowdown feed system is adopted.Results show the superior performance of the HP/PE combination. Propellant consumption and propulsion system dry-mass are competing parameters. The LOX/HTPB propellant combination presents higher maximum values ofc*, but this maximum is reached for small values of α: large grains are required and the propulsion system mass penalty is higher than propellant savings. Also, LOX and HP have larger densities than HTPB, and higher values of density-specific impulse may be obtained with the HP/PE combination which adopts high values of α. On the other hand, N2O is penalized by the low performance in terms of specific impulse. However, N2O presents some interesting features, such as the reduced rocket length L and diameter D, as the large mixture ratio reduces the grain length. Also note that the initial required thrust for the same mission is lower.


## 4. Fuels with Enhanced Regression Rate

In this section approaches are considered based on the modification of the fuel. Possible approaches can be divided down into three categories:(1) adding energetic particles, (2) using energetic polymers (and/or plasticisers) instead of conventional inert components such as HTPB, (3) using fuels such as cryogenic fuels or paraffin-based fuels which exhibit new mass-transfer mechanisms. A combination of the aforementioned options is also possible.


### 4.1. Particle Additives in Solid Fuels

Several additives have been tested. Most researches considered metal additives in polymeric-based solid fuels [6, 8, 9, 27, 41–45]. More recently, metal additives in paraffin-based fuels have also been considered [8, 46–50].Combustion of metals has inherent advantages as they help to increase both specific impulse and density specific impulse. In fact, even if their condensed combustion products determine multiphase losses, they have high reaction heat and high density. Being usually in powdered form they can easily be mixed in solid fuels. In the case of HREs it has been proven that metal additives may enhance regression rate. The addition of various types of metals (e.g., Al, LiH, B, W, Mg) into solid fuels has been tested. Most investigations considered aluminum, lithium, and boron.Regression rate is influenced by component characteristics (e.g., reactivity, heat of oxidation, density) and particle features (e.g., shape, dimension, coating). Aluminum may be considered the most researched additive. It has a high heat of oxidation, high density and ease of ignition. Lithium-based particles are highly reactive, but the heat of oxidation is low compared to more common particles such aluminum. Finally, boron has a very high heat of combustion but is not so easy to ignite, which is of concern in limited volume propulsion systems. As far as dimension is concerned, microsized particles have been extensively tested. More recently nanosized particles have been used. Nanosized particles have a very high surface to volume ratio. Compared to microsized particles, nanosized particles present reduced ignition and burning time, present a more complete combustion (highc* efficiency),  and determine a larger increment od regression rate, but they are more expensive and complicate the grain production process.Nonvolatile fuel particles (e.g., aluminum) in the grain have several effects, influencing both physical properties of the solid and the combustion process. Reasons behind some test results are currently under investigation. The main known effects due to nonvolatile particles are as follows.(i)
Particles such as metals have higher density than fuel binders, and the grain mean density changes with the percentage of weight addedk. It is useful to consider both the regression rate r and the mass burning rate ρFr improvement. For example, if a weight fraction of 20% of ultrafine aluminum (UFAL) is added in HTPB, the regression rate is increased by about 40%, whereas the mass burning rate increases by almost 70% over that of pure HTPB [27].(ii)
The blocking effect on surface heat transfer depends only on the gas blowing rate. Referring only to the volatile fuel mass flow, the densityρv=(1-k)ρF should be used instead of ρF in the previous equations (1) and (2). It follows that, disregarding other effects, the regression rate should be nearly inversely proportional to (1-k).(iii)
The value ofΔHv,eff (here evaluated on a volatile binder mass basis) is higher  than the value of the volatile binder alone ΔHv, as further energy is needed to heat the nonvolatile particle to the surface temperature. Nevertheless, due to the relatively low heat capacity of most metal additives, ΔHv,eff does not differ much from ΔHv. The new value of ΔHv,eff should be used to evaluate B, but the resulting effect is mitigated by the low exponent of B in (2). (Note that a different point of view is often adopted in the literature: reference is made to the whole mass of the fuel grain. In this case  ΔHv,eff is decreased as the nonvolatile particle should be only heated to surface temperature and not gasified.)(iv)
The radiant energy flux coming from condensed particles is large. This effect is mitigated by the coupling effect between radiant and convective heat fluxes.(v)
Additives enhance flame temperature andΔHfw increases. If particle size is large (some μm) the nonvolatile particles perform distributed combustion releasing some heat far from the fuel surface. On the contrary, nanometric particles burn more rapidly and closer to the fuel surface. Also microexplosions may be experienced. Moreover particle coating affects combustion. The new value of ΔHfw should be used to evaluate B, but the resulting effect is mitigated by the low exponent of B in (2).(vi)
The removal of solid from the grain surface may become the rate-limiting process for fuel rate regression [51]. The presence of flakes/agglomerates has been observed on the surface and in the combustion zone. Metals and their oxide may accumulate on the surface and leave it as agglomerates. At high mass fluxes the decomposition of fuel binders such as HTPB occurs very quickly and soot is formed, possibly due to desorption process of HTPB fragments at the grain surface.(vii)
Some of the oxidizer is consumed by metal particles burning and converted to solid oxides, thus reducing local gaseous mass-flux.Beside metals, also oxidizer particles such as ammonium perchlorate (AP) have been used [52–55]. The presence of an oxidizer in the solid grain determines near-surface/heterogeneous reactions. In addition the flame in the boundary layer approaches the surface, thus further enhancing heat transfer. Consequently, the regression rate increases. Nevertheless, the inherently safe inert behavior of the solid grain is compromised.


### 4.2. Energetic Compounds

Energetic formulations can be used in solid fuels to replace low-energy conventional binders such as HTPB. Both polymers and/or plasticizers which contribute to the overall energy of the composition can be used.High nitrogen ingredients (HiN) show great promise [56, 57]. Despite low decomposition temperatures, energetic materials such as triaminoguanidinium azotetrazolate (TAGzT) present positive heat of formation. Containing no oxygen, the dominant chemistry is more centered on condensed phase reactions and much faster decomposition rates are obtained when compared to HMX. A 25% increment in regression rate with respect to a similar HTPB grain has been observed using a 25% wt addition of TAGzT [58]. If nanoaluminum is added in addition to HiN ingredients, Al reacts with the high-temperature nitrogen liberated from the decomposition of HiN material to form AlN near the surface where oxygen is absent.Azido-substituted polymers such as glycidyl azide polymer (GAP) [59], 3,3-bis-(azidomethyl)oxetane (BAMO) and 3-azidomethyl-3-methyl oxetane (AMMO) are effective as the polymer matrix. However, these compounds begin to transition the inert fuel to a fuel-rich solid propellant and the safe nature of classical HREs may be compromised.


### 4.3. Paraffin Fuels

Researchers at Stanford University [60] discovered that paraffin-based fuels exist which have regression rates that are 3-4 times than those of conventional hybrid fuels. This is mainly due to the production of a thin liquid layer on the fuel grain surface which becomes instable: If this layer is characterized by low viscosity μ and low surface tension σ, instability driven by the oxidizer flow may arise and liquid fuel droplets are injected into the boundary layer [5, 61]. This mass-transfer mechanism is similar to a spray injection and does not depend on heat transfer. It enhances fuel mass flow without the blocking effect typical of gaseous fuel blowing. The entrained mass flow ṁentr depends on dynamic pressure pdyn, layer thickness h, and liquid properties (i.e., μ and σ) at the temperature of the melt layer (9)ṁentr∝pdynβhδσθμλ.
As shown in (9), entrainment occurs only when the liquid in the melt layer has low values of μ and σ. For example, high-density polyethylene (HDPE) does form a melt layer, but μ is four orders of magnitude larger than paraffin and droplet entrainment is not significant. In contrast, fuels such as paraffin waxes and polyethylene waxes have low viscosity. These fuels also present other advantages, as they are low-cost, easy to process, and environmentally benign. They may give a specific impulse similar to kerosene, but with higher density and without the combustion product toxicity typical of SRM grains. A blend of different components can be used to mitigate this issue and to tailor regression rate for a given mission. For example, ethylene-vinyl acetate copolymers have been mixed with paraffin-wax to change the melted-fuel viscosity [62]. According to experimental data, regression rate increases almost in proportion with μ-1/6. Beside all these advantages, one issue is the low mechanical strength. Additives may be used to improve this quality, but they may affect viscosity and surface tension adversely.Additives to further enhance the regression rate have also been used. As an example high regression rate has been experienced adding Al particles in paraffin-based fuels: 13% wt of aluminum particles (Silberline) determine a 30% regression rate enhancement [8]. It is believed that aluminum particles coming from the surface may be encapsulated by the wax, thus helping ignition. Regression rate usually increases with loaded weight of additive. This is not the case of some tests with ammonia borane [48], where regression rate was enhanced with small weight fraction, while decreased for higher weight fractions.


### 4.4. Cryogenic Fuels

Many fuels such as methane, ethylene, pentane, and RP1 have been tested as solid cryogenic grains [63–65]. They proved to regress 2 to 10 times faster than storable fuels in the same operating conditions. Solid methane/gaseous oxygen (GOX) proved to regress at nearly 8 mm/s with an oxidizer mass flux of 20 kg/m2/s. Also, solid oxidizers such as CO and O2 have been tested (reverse HRE configuration) [66]. The high regression rate obtained is explained by the theories developed to describe high regression rates in paraffin fuels [5, 61]. Results are encouraging, but complications arise due to the low temperatures, and the production of large grain is difficult.


## 4.1. Particle Additives in Solid Fuels

Several additives have been tested. Most researches considered metal additives in polymeric-based solid fuels [6, 8, 9, 27, 41–45]. More recently, metal additives in paraffin-based fuels have also been considered [8, 46–50].Combustion of metals has inherent advantages as they help to increase both specific impulse and density specific impulse. In fact, even if their condensed combustion products determine multiphase losses, they have high reaction heat and high density. Being usually in powdered form they can easily be mixed in solid fuels. In the case of HREs it has been proven that metal additives may enhance regression rate. The addition of various types of metals (e.g., Al, LiH, B, W, Mg) into solid fuels has been tested. Most investigations considered aluminum, lithium, and boron.Regression rate is influenced by component characteristics (e.g., reactivity, heat of oxidation, density) and particle features (e.g., shape, dimension, coating). Aluminum may be considered the most researched additive. It has a high heat of oxidation, high density and ease of ignition. Lithium-based particles are highly reactive, but the heat of oxidation is low compared to more common particles such aluminum. Finally, boron has a very high heat of combustion but is not so easy to ignite, which is of concern in limited volume propulsion systems. As far as dimension is concerned, microsized particles have been extensively tested. More recently nanosized particles have been used. Nanosized particles have a very high surface to volume ratio. Compared to microsized particles, nanosized particles present reduced ignition and burning time, present a more complete combustion (highc* efficiency),  and determine a larger increment od regression rate, but they are more expensive and complicate the grain production process.Nonvolatile fuel particles (e.g., aluminum) in the grain have several effects, influencing both physical properties of the solid and the combustion process. Reasons behind some test results are currently under investigation. The main known effects due to nonvolatile particles are as follows.(i)
Particles such as metals have higher density than fuel binders, and the grain mean density changes with the percentage of weight addedk. It is useful to consider both the regression rate r and the mass burning rate ρFr improvement. For example, if a weight fraction of 20% of ultrafine aluminum (UFAL) is added in HTPB, the regression rate is increased by about 40%, whereas the mass burning rate increases by almost 70% over that of pure HTPB [27].(ii)
The blocking effect on surface heat transfer depends only on the gas blowing rate. Referring only to the volatile fuel mass flow, the densityρv=(1-k)ρF should be used instead of ρF in the previous equations (1) and (2). It follows that, disregarding other effects, the regression rate should be nearly inversely proportional to (1-k).(iii)
The value ofΔHv,eff (here evaluated on a volatile binder mass basis) is higher  than the value of the volatile binder alone ΔHv, as further energy is needed to heat the nonvolatile particle to the surface temperature. Nevertheless, due to the relatively low heat capacity of most metal additives, ΔHv,eff does not differ much from ΔHv. The new value of ΔHv,eff should be used to evaluate B, but the resulting effect is mitigated by the low exponent of B in (2). (Note that a different point of view is often adopted in the literature: reference is made to the whole mass of the fuel grain. In this case  ΔHv,eff is decreased as the nonvolatile particle should be only heated to surface temperature and not gasified.)(iv)
The radiant energy flux coming from condensed particles is large. This effect is mitigated by the coupling effect between radiant and convective heat fluxes.(v)
Additives enhance flame temperature andΔHfw increases. If particle size is large (some μm) the nonvolatile particles perform distributed combustion releasing some heat far from the fuel surface. On the contrary, nanometric particles burn more rapidly and closer to the fuel surface. Also microexplosions may be experienced. Moreover particle coating affects combustion. The new value of ΔHfw should be used to evaluate B, but the resulting effect is mitigated by the low exponent of B in (2).(vi)
The removal of solid from the grain surface may become the rate-limiting process for fuel rate regression [51]. The presence of flakes/agglomerates has been observed on the surface and in the combustion zone. Metals and their oxide may accumulate on the surface and leave it as agglomerates. At high mass fluxes the decomposition of fuel binders such as HTPB occurs very quickly and soot is formed, possibly due to desorption process of HTPB fragments at the grain surface.(vii)
Some of the oxidizer is consumed by metal particles burning and converted to solid oxides, thus reducing local gaseous mass-flux.Beside metals, also oxidizer particles such as ammonium perchlorate (AP) have been used [52–55]. The presence of an oxidizer in the solid grain determines near-surface/heterogeneous reactions. In addition the flame in the boundary layer approaches the surface, thus further enhancing heat transfer. Consequently, the regression rate increases. Nevertheless, the inherently safe inert behavior of the solid grain is compromised.


## 4.2. Energetic Compounds

Energetic formulations can be used in solid fuels to replace low-energy conventional binders such as HTPB. Both polymers and/or plasticizers which contribute to the overall energy of the composition can be used.High nitrogen ingredients (HiN) show great promise [56, 57]. Despite low decomposition temperatures, energetic materials such as triaminoguanidinium azotetrazolate (TAGzT) present positive heat of formation. Containing no oxygen, the dominant chemistry is more centered on condensed phase reactions and much faster decomposition rates are obtained when compared to HMX. A 25% increment in regression rate with respect to a similar HTPB grain has been observed using a 25% wt addition of TAGzT [58]. If nanoaluminum is added in addition to HiN ingredients, Al reacts with the high-temperature nitrogen liberated from the decomposition of HiN material to form AlN near the surface where oxygen is absent.Azido-substituted polymers such as glycidyl azide polymer (GAP) [59], 3,3-bis-(azidomethyl)oxetane (BAMO) and 3-azidomethyl-3-methyl oxetane (AMMO) are effective as the polymer matrix. However, these compounds begin to transition the inert fuel to a fuel-rich solid propellant and the safe nature of classical HREs may be compromised.


## 4.3. Paraffin Fuels

Researchers at Stanford University [60] discovered that paraffin-based fuels exist which have regression rates that are 3-4 times than those of conventional hybrid fuels. This is mainly due to the production of a thin liquid layer on the fuel grain surface which becomes instable: If this layer is characterized by low viscosity μ and low surface tension σ, instability driven by the oxidizer flow may arise and liquid fuel droplets are injected into the boundary layer [5, 61]. This mass-transfer mechanism is similar to a spray injection and does not depend on heat transfer. It enhances fuel mass flow without the blocking effect typical of gaseous fuel blowing. The entrained mass flow ṁentr depends on dynamic pressure pdyn, layer thickness h, and liquid properties (i.e., μ and σ) at the temperature of the melt layer (9)ṁentr∝pdynβhδσθμλ.
As shown in (9), entrainment occurs only when the liquid in the melt layer has low values of μ and σ. For example, high-density polyethylene (HDPE) does form a melt layer, but μ is four orders of magnitude larger than paraffin and droplet entrainment is not significant. In contrast, fuels such as paraffin waxes and polyethylene waxes have low viscosity. These fuels also present other advantages, as they are low-cost, easy to process, and environmentally benign. They may give a specific impulse similar to kerosene, but with higher density and without the combustion product toxicity typical of SRM grains. A blend of different components can be used to mitigate this issue and to tailor regression rate for a given mission. For example, ethylene-vinyl acetate copolymers have been mixed with paraffin-wax to change the melted-fuel viscosity [62]. According to experimental data, regression rate increases almost in proportion with μ-1/6. Beside all these advantages, one issue is the low mechanical strength. Additives may be used to improve this quality, but they may affect viscosity and surface tension adversely.Additives to further enhance the regression rate have also been used. As an example high regression rate has been experienced adding Al particles in paraffin-based fuels: 13% wt of aluminum particles (Silberline) determine a 30% regression rate enhancement [8]. It is believed that aluminum particles coming from the surface may be encapsulated by the wax, thus helping ignition. Regression rate usually increases with loaded weight of additive. This is not the case of some tests with ammonia borane [48], where regression rate was enhanced with small weight fraction, while decreased for higher weight fractions.


## 4.4. Cryogenic Fuels

Many fuels such as methane, ethylene, pentane, and RP1 have been tested as solid cryogenic grains [63–65]. They proved to regress 2 to 10 times faster than storable fuels in the same operating conditions. Solid methane/gaseous oxygen (GOX) proved to regress at nearly 8 mm/s with an oxidizer mass flux of 20 kg/m2/s. Also, solid oxidizers such as CO and O2 have been tested (reverse HRE configuration) [66]. The high regression rate obtained is explained by the theories developed to describe high regression rates in paraffin fuels [5, 61]. Results are encouraging, but complications arise due to the low temperatures, and the production of large grain is difficult.


## 5. Nonconventional Grain Geometries and Injectors

Methods can be used that modify classical flow structure in order to improve heat transfer. These approaches may be based on nonconventional grain geometry, nonconventional injector design, or both. Nonconventional geometries, for example, radial and/or end burning of the grain, may determine 3D flows with detachments and recirculation. Swirl injectors are able to reduce both thickness and growth of the boundary layer, thus enhancing heat transfer. The heat transfer variation is reduced and regression rate is more uniform.


### 5.1. Head-End Swirl in Side-Burning Grains

Several researchers tested swirl injector at the fore end of conventional cylindrical grains [11, 67]. For instance, Yuasa et al.  [11] reports results for different grain length, geometric swirl factor, and GO. Swirl strength and oxidizer mass flux were varied independently. An axial injector was also used to obtain a baseline correlation. A regression rate 2.7 times greater than that for the baseline case was experienced using GOX and PMMA. Using the correlation r=aGOn it has been found that the exponent n is not affected significantly by the swirl number, but changes with the range of mass flux (n~0.8 for large values of GO according to the theoretical turbulent layer heat transfer, while n~0.6 at lower GO values, probably due to radiation effects). On the contrary, the value of a increases almost linearly with swirl factor.Some researchers also modified the grain geometry [67, 68]. Lee et al. [67] considered grooves in the grain and swirl injectors. Separated and combined effects where analyzed. Swirl injectors alone are able to have double-regression rate, while grooves have a minor effect if used alone. Grooves are supposed to enhance the turbulence in the near-surface region and allow for an increased burning surface of more than 100%, but the complexity and the cost associated with grain manufacture must be taken into account. In [68] a tapered grain is proposed to avoid swirl decay along the grain axis. Conservation of angular momentum accelerates the tangential velocity because of section reduction and the swirl angle may be maintained.It is to be pointed out that some effects of swirling must be addressed, including torque and effects of nonaxial flow in the nozzle (effective throat area and divergence losses). This holds for any of the concepts employing swirl oxidizer injection.


### 5.2. Aft-End Swirl in Side-Burning Grains: The Vortex Hybrid

The key characteristic of vortex hybrid [10, 69, 70] is a unique coaxial, coswirling, counterflowing vortex pair which has been found to improve regression rate. A classical cylindrical grain with circular port section is used. The vortex pair is obtained by injecting the oxidizer through a swirl injector located between the aft end of the fuel grain and the inlet of the exit nozzle. The oxidizer is prevented from directly flowing out of the nozzle by centrifugal forces and pressure gradients. An outer vortex is formed which spirals toward the engine head end. At the engine head end, the outer vortex turns inward and transforms into an inner vortex that spirals toward the nozzle. This flow structure presents many advantages: the residence time is increased, mixing is enhanced, both thickness and growth of the boundary layer are reduced, thus augmenting the heat transfer to the fuel surface and making the regression rate very uniform. Tests have been carried out using GOX with HTPB and other proprietary fuels. HTPB shows a regression rate up to 6 times faster that those in HREs with classical head end injectors. Using the classical correlation of (3), it has been found [14]: (10)r=0.193GO0.54.
The empirical power of 0.54 on mass flux, in conjunction with the single circular port and the uniform regression rate along the axis, guarantees a reduced mixture ratio shifting. On the other hand, the regression rate depends on the injection velocity of the oxidizer. Maximum mass flux tested [10] was about 100 kg/m2 s. Scale effects must be investigated.


### 5.3. Radial-Flow HREs

An alternative, which may improve fuel section volumetric efficiency, is the radial flow hybrid rocket. Multiple fuel plates stacked on top of each other can provide very large burn surface areas and high volumetric efficiency if there is a little gap in between the plates. This configuration can be a good option when a smallL/D is important as in launcher upper stages. It can be also used for space propulsion [71]. The flow can be injected from a central hole in a disk, or radially. The measured regression rates are greater than those observed in axial port designs at the same flux and chamber pressure levels [12]. However, detailed regression on both upper and lower fuel disks is very complex, being influenced by flow impingement, flow separation, and 3D flows that cause nonuniform regression rate, determining the presence of sliver. In addition, structural integrity may be of concern if multiple disks are used.Swirl have been used in radial flow hybrids. One example is the so called vortex flow pancake (VFP) [72] where the swirl injectors are between two end-burning disks of fuel. It has been found that the regression rate is highly dependent on flux. The regression rate appears to be one order of magnitude larger that the one obtained in classical HRE for similar mass fluxes. Nevertheless, the mass flux used for these tests are small with respect to typical values in classical HRE tests: Scale effects must be evaluated.


### 5.4. Side/End-Burning Grains: The Cascaded Multistage Impinging-Jet

The cascaded multistage impinging-jet (CAMUI) is a method proposed to improve thrust level of hybrid rockets without changing fuel material or introducing nonconventional injectors [13, 73]. This concept is based on a unique geometry of the fuel grain which consists of several stages of cylindrical blocks with two axial ports. The jet of the combustion gas collides with the burning surfaces repeatedly, resulting in the high regression rate of the solid fuel. Both side burning and lateral end burning is used, thus reducing grain length for a given Ab. In fact, each block has three burning surfaces, namely, the port surfaces and the two lateral ends (upstream and downstream end faces). To reduce sliver and avoid grain collapsing, the ports must reach the external grain diameter when the upstream face meets the downstream end face. The three surfaces have different regression and predicting grain geometry evolution and fuel flow rate is not an easy task. The development of accurate regression rate correlation is needed. Tests have been performed using GOX/PE propellant combination. The regression rate for port surface is similar to (3): (11)r=a(α′)Gpn,
but in this case a is a function of the local mixture ratio α′ and the propellant mass flux Gp is used. The regression rate of forward end faces depends on the spacing between blocks H and the port area diameter Dp of the upstream block: (12)r(H/Dp)m=a(α′)Gpn.
The downstream end face has a more complicated behavior, presenting non uniform regression rates. Two different regression rate correlations are introduced as wall jets enhance the regression in the central part of the face. Due to the grain configuration low values of L/D can be obtained, but volumetric efficiency and sliver may be issues.


## 5.1. Head-End Swirl in Side-Burning Grains

Several researchers tested swirl injector at the fore end of conventional cylindrical grains [11, 67]. For instance, Yuasa et al.  [11] reports results for different grain length, geometric swirl factor, and GO. Swirl strength and oxidizer mass flux were varied independently. An axial injector was also used to obtain a baseline correlation. A regression rate 2.7 times greater than that for the baseline case was experienced using GOX and PMMA. Using the correlation r=aGOn it has been found that the exponent n is not affected significantly by the swirl number, but changes with the range of mass flux (n~0.8 for large values of GO according to the theoretical turbulent layer heat transfer, while n~0.6 at lower GO values, probably due to radiation effects). On the contrary, the value of a increases almost linearly with swirl factor.Some researchers also modified the grain geometry [67, 68]. Lee et al. [67] considered grooves in the grain and swirl injectors. Separated and combined effects where analyzed. Swirl injectors alone are able to have double-regression rate, while grooves have a minor effect if used alone. Grooves are supposed to enhance the turbulence in the near-surface region and allow for an increased burning surface of more than 100%, but the complexity and the cost associated with grain manufacture must be taken into account. In [68] a tapered grain is proposed to avoid swirl decay along the grain axis. Conservation of angular momentum accelerates the tangential velocity because of section reduction and the swirl angle may be maintained.It is to be pointed out that some effects of swirling must be addressed, including torque and effects of nonaxial flow in the nozzle (effective throat area and divergence losses). This holds for any of the concepts employing swirl oxidizer injection.


## 5.2. Aft-End Swirl in Side-Burning Grains: The Vortex Hybrid

The key characteristic of vortex hybrid [10, 69, 70] is a unique coaxial, coswirling, counterflowing vortex pair which has been found to improve regression rate. A classical cylindrical grain with circular port section is used. The vortex pair is obtained by injecting the oxidizer through a swirl injector located between the aft end of the fuel grain and the inlet of the exit nozzle. The oxidizer is prevented from directly flowing out of the nozzle by centrifugal forces and pressure gradients. An outer vortex is formed which spirals toward the engine head end. At the engine head end, the outer vortex turns inward and transforms into an inner vortex that spirals toward the nozzle. This flow structure presents many advantages: the residence time is increased, mixing is enhanced, both thickness and growth of the boundary layer are reduced, thus augmenting the heat transfer to the fuel surface and making the regression rate very uniform. Tests have been carried out using GOX with HTPB and other proprietary fuels. HTPB shows a regression rate up to 6 times faster that those in HREs with classical head end injectors. Using the classical correlation of (3), it has been found [14]: (10)r=0.193GO0.54.
The empirical power of 0.54 on mass flux, in conjunction with the single circular port and the uniform regression rate along the axis, guarantees a reduced mixture ratio shifting. On the other hand, the regression rate depends on the injection velocity of the oxidizer. Maximum mass flux tested [10] was about 100 kg/m2 s. Scale effects must be investigated.


## 5.3. Radial-Flow HREs

An alternative, which may improve fuel section volumetric efficiency, is the radial flow hybrid rocket. Multiple fuel plates stacked on top of each other can provide very large burn surface areas and high volumetric efficiency if there is a little gap in between the plates. This configuration can be a good option when a smallL/D is important as in launcher upper stages. It can be also used for space propulsion [71]. The flow can be injected from a central hole in a disk, or radially. The measured regression rates are greater than those observed in axial port designs at the same flux and chamber pressure levels [12]. However, detailed regression on both upper and lower fuel disks is very complex, being influenced by flow impingement, flow separation, and 3D flows that cause nonuniform regression rate, determining the presence of sliver. In addition, structural integrity may be of concern if multiple disks are used.Swirl have been used in radial flow hybrids. One example is the so called vortex flow pancake (VFP) [72] where the swirl injectors are between two end-burning disks of fuel. It has been found that the regression rate is highly dependent on flux. The regression rate appears to be one order of magnitude larger that the one obtained in classical HRE for similar mass fluxes. Nevertheless, the mass flux used for these tests are small with respect to typical values in classical HRE tests: Scale effects must be evaluated.


## 5.4. Side/End-Burning Grains: The Cascaded Multistage Impinging-Jet

The cascaded multistage impinging-jet (CAMUI) is a method proposed to improve thrust level of hybrid rockets without changing fuel material or introducing nonconventional injectors [13, 73]. This concept is based on a unique geometry of the fuel grain which consists of several stages of cylindrical blocks with two axial ports. The jet of the combustion gas collides with the burning surfaces repeatedly, resulting in the high regression rate of the solid fuel. Both side burning and lateral end burning is used, thus reducing grain length for a given Ab. In fact, each block has three burning surfaces, namely, the port surfaces and the two lateral ends (upstream and downstream end faces). To reduce sliver and avoid grain collapsing, the ports must reach the external grain diameter when the upstream face meets the downstream end face. The three surfaces have different regression and predicting grain geometry evolution and fuel flow rate is not an easy task. The development of accurate regression rate correlation is needed. Tests have been performed using GOX/PE propellant combination. The regression rate for port surface is similar to (3): (11)r=a(α′)Gpn,
but in this case a is a function of the local mixture ratio α′ and the propellant mass flux Gp is used. The regression rate of forward end faces depends on the spacing between blocks H and the port area diameter Dp of the upstream block: (12)r(H/Dp)m=a(α′)Gpn.
The downstream end face has a more complicated behavior, presenting non uniform regression rates. Two different regression rate correlations are introduced as wall jets enhance the regression in the central part of the face. Due to the grain configuration low values of L/D can be obtained, but volumetric efficiency and sliver may be issues.


## 6. Final Remarks

The low regression rate of fuel grain is one of the most challenging issues for hybrid rockets. Many solutions have been presented in the literature, demonstrating the vital research activity in hybrid rocket propulsion. The most important principles and investigated solutions have been reviewed here. The challenge is to eliminate or, at least, mitigate the negative effects of the low regression rate, without compromising the appealing features of HREs (e.g., safety, low cost, simplicity, environmental friendliness).Propellant combinations which require high mixture ratios reduce the contribution grain fuels make. In this case the choice of the oxidizer plays a significant role. Proper tools which couple engine design and trajectory optimization are needed to find the best propellant combination to be used. In order to have acceptable rocket geometries, multiport grains have to be adopted for missions requiring high performance. However, the drawbacks of multiport grains, such as sliver and structural concerns, drive the researchers to investigate means which enhance the regression rate.Approaches which increase regression rate have been presented in two sections of this paper. A first section has been devoted to advanced fuels, considering both additives (mainly metal particles) and nonconventional fuels, that is, paraffin-based fuels and cryogenic grains. A second section considered methods which aim at increasing the heat transfer to the grain surface, modifying grain and/or injector geometry. Figure4 and Table 1 compare the performance of approaches reviewed here. Dashed lines shows conventional HRE regression rate for comparison. Higher regression rates are usually obtained with respect to classical HREs at similar operating conditions. Most of the solutions presented remain at a laboratory or concept scale and need further research efforts. Nevertheless, the data presented in the literature show that promising approaches do exist and may be adopted/combined to give HREs the opportunity to become a very competitive propulsion system.Table 1
Values ofa and n to be used in (3)  (r in m/s and GO in kg/m2 s).
LegendSystemPropellantsa(m1+2nkg-nsn-1)nNoteReferenceGO (kg/m2/s)Fuel/additivesBPure HTPBGOX/HTPB2.85·10-50.681—[4]35–280WParaffinGOX/Wax9.10·10-50.690—[5]20–120WAlParaffin/13% SilbalGOX/fuel9.40·10-50.766—[6]150–300CCryoGOX/CH44.14·10-50.830—[7]3–30S1Pure HTPBGOX/fuel——baseline, [8, Figure 11][8]80–150S12HTPB/13% Al325GOX/fuel——[8, Figure 11][8]80–120S13HTPB/13% ALEXGOX/fuel——coated, [8, Figure 11][8]70–120HbPure HTPBGOX/HTPB8.7·10-50.530baseline[9]50–400HAlHTPB/AlGOX/fuel1.4·10-50.930—[9]50–400HAPHTPB/APGOX/fuel3.8·10-50.710—[9]50–400HAPAlHTPB/Al/APGOX/fuel1.2·10-50.97—[9]50–400Swirl/grain geometryVXVortexGOX/HTPB1.93·10-40.540—[10]60–110VEnd-burning + SwirlGOX/PMMA3.45·10-50.778—[11]40–60VbEnd-burning + SwirlGOX/PMMA1.45·10-50.749no swirl[11]40–60VEnd-burning + SwirlGOX/PMMA5.96·10-50.641—[11]10–20VbEnd-burning + SwirlGOX/PMMA2.76·10-50.581no swirl[11]10–20RURadialGOX/HTPB9.20·10-50.570upper disk[12]30–70RLRadialGOX/HTPB1.00·10-40.700lower disk[12]30–70CACAMUIGOX/PE2.40·10-50.800(11), α′=3[13]200–700Table 2
Optimal design and performance of sounding rocket for different propellant combinations.
Propellantsa (m1+2nkg-n sn-1)nReferenceFi (kN)αimp (kg)D (m)L (m)tμg (s)HP/PE7.00·10-60.800[38]24.88.573390.425.01299LOX/HTPB9.24·10-60.852[39]23.53.203280.465.52219N2O/HTPB1.87·10-40.347[40]16.210.923400.404.75177Figure 4
Regression rate for different systems and GOX/fuel combinations (see legend in Table1).Multiport geometries.
(a)
Triangular-port geometry(b)
Quadrilateral-port geometry



---

*Source: 649753-2012-06-21.xml*