Next Article in Journal
Flexural Properties of Wet-Laid Hybrid Nonwoven Recycled Carbon and Flax Fibre Composites in Poly-Lactic Acid Matrix
Next Article in Special Issue
Fabrication and Testing of the Cold Gas Propulsion System Flight Unit for the Adelis-SAMSON Nano-Satellites
Previous Article in Journal
Crack Growth in a Range of Additively Manufactured Aerospace Structural Materials
Previous Article in Special Issue
A Dual Mode Propulsion System for Small Satellite Applications
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Single-Use Microthruster Concept for Small Satellite Attitude Control in Formation-Flying Applications

1
Department of Mechanical Engineering, University of Vermont, Burlington, VT 05405, USA
2
Benchmark Space Systems Inc., South Burlington, VT 05403, USA
*
Author to whom correspondence should be addressed.
Aerospace 2018, 5(4), 119; https://doi.org/10.3390/aerospace5040119
Submission received: 5 June 2018 / Revised: 29 October 2018 / Accepted: 5 November 2018 / Published: 14 November 2018

Abstract

:
In recent years, the maturation of small satellite technology has led to their adoption for a variety of space missions. The next generation of small satellite missions, however, will likely have the satellites operating in formations or “constellations” to perform missions that are not currently possible. A key enabling technology for constellation-based missions is a miniaturized propulsion system that is capable of delivering the extremely low impulse levels required for maintaining precise relative position and orientation. Existing propulsion solutions for this regime suffer from compromises on power, safety, and cost that have limited their adoption. In this work, we describe a new, low-power micropropulsion concept based on the thermal decomposition of an inert chemical blowing agent (CBA) as the propellant. A meso-scale prototype device is designed, fabricated, and tested. The experimental results indicate that this concept, when appropriately scaled, is capable of providing thrust levels (∼1 μ N) and impulse-bits (∼0.1 μ N·s) that are commensurate with the intended application.

Graphical Abstract

1. Introduction

As the capabilities of small satellites—defined by NASA as having masses less than 180 kg—have matured, mission designers have begun to consider their use in formation flying mission architectures that would be impossible and/or cost-prohibitive to perform with larger satellites. Two immediate examples are multi-point Earth observations or synthetic aperture arrays for deep space exploration; an extensive compilation of formation flying missions is found in the recent review by Bandyopadhyay et al. [1] Small satellites are uniquely suited for these types of missions, as all of the elements can be launched simultaneously and for a fraction of the cost of traditional missions. There has been particular interest in missions designed around the CubeSat platform since the supporting infrastructure for launching and deploying satellites built to this standard is already in place. Indeed, a number of CubeSat-based missions are in development that demonstrate the potential of the platform and the demand for continued improvement of the supporting technologies [2,3,4].
Much of the formation flying literature has focused, and rightly so, on the aspects of dynamics and control. However, the implementation of such control strategies hinges upon the existence and availability of the associated satellite technology and hardware. Although current small satellites designs are considerably more capable than those in previous generations, the hardware and system technology still lag mission concepts; indeed, additional satellite subsystem development is required to support new mission concepts. A recent review of the state of small satellite technology is given in the 2015 NASA Ames Small Spacecraft Technology State of the Art technical report [5]. Perhaps the most critical of these subsystems is that of propulsion; this subsystem must be capable of providing extremely low level impulses needed for the relative position and orientation control necessary to enable precise formation flying [6,7,8]. Small satellite attitude control thrusters are particularly challenging as they must provide reliable, low impulse-bit operation while conforming to the size, weight, power, and cost constraints of the form factor. In addition to the technical challenges these propulsion systems must address, there are often regulatory hurdles that impede their adoption for small satellites [9,10]. These regulations require propulsion systems to meet range safety and secondary payload requirements that limit the storage tank pressurization, amount of stored chemical energy, and toxicity of the propellant. Such requirements immediately eliminate many chemical-based propulsion options, and limit the efficacy of others.
We note in passing that electric propulsion concepts have also been proposed in recent years for CubeSat platforms [11,12,13,14]. These include electrospray thrusters, pulsed-plasma thrusters (PPT) and field-emission electric propulsion thrusters (FEEP) among others. Readers interested in non-chemical options may wish to consult the cited references for greater details. Based on its inherently low thrust levels and high specific impulse, electric propulsion would seem to be a natural consideration for formation-flying applications. However, these systems pose challenges in terms of system mass and power requirements that raise questions as to their viability for the CubeSat platform. As such, we restrict our focus here to chemical-based options.
A new small satellite propulsion attitude control concept is currently in development at the University of Vermont with the above limitations and considerations in mind. The proposed system consists of a physical array of discrete, single-use microthrusters, each containing a non-toxic solid propellant based on a chemical blowing agent (CBA). The operation of the thruster is triggered by the thermal decomposition of the CBA into hot, pressurized gas. The packaging into an array permits multiple firings by the system. Based on the goal of providing thrust levels (∼1 μ N) and impulse-bits (∼0.1 μ N·s) necessary for small attitude adjustments, our design calculations show that this array would require MEMS-manufacturing techniques. This paper describes the design concept for the system, the development of meso-scale proof-of-concept thruster, and its performance characteristics under testing. The experimental results are supplemented by numerical simulations of the thruster firing to better understand the observed performance within the context of the transient gas dynamics within the device. The results of these efforts point to the viability of a system that can be mass-produced at a low-cost for future small satellite formation flying missions.

2. Design Concept

The proposed system is intended to offer relative position and attitude control of a CubeSat-class or smaller spacecraft operating in a constellation architecture. The target performance characteristics for the design are shown in Table 1; these performance parameters are selected to provide the satellite with the ability to maintain relative orientation and position in the face of disturbance torques and/or to allow for re-orientation with milli-arcsecond resolution.
To address the regulatory challenges typically associated with chemical propellants, this system is designed to use a chemical blowing agent as the propellant. Blowing agents are a family of chemical compounds that produce a significant volume of gas as a result of thermally-induced decomposition. They are commonly used as an additive in polymer manufacturing and food production to reduce the density of the surrounding medium [15,16]. They are attractive as a potential propellant since they are non-toxic, inert at temperatures below the activation temperature. Moreover, CBAs are low in cost due to their extensive use in other industries.
To achieve the minimum thrust and impulse bits required, our thruster concept envisions an individual thruster “element” consisting of a 32 × 32 × 10 μ m propellant chamber/plenum, paired with an electro-thermal heating element and a converging-diverging nozzle. The concept further envisions that the individual thrusters would be packaged as an array of microthrusters with individual decomposition chambers, with each array being placed at the corners of the CubeSat to offer control (Figure 1). Based on current MEMS-manufacturing capabilities, it is deemed feasible to create a 7 mm square array of these cells, with spacing for structural and thermal integrity, consisting of 19,600 individual cells, each designed to provide ∼1 μ N of thrust.
The operation of the microthruster array is straightforward. Every element in the array can be addressed and ignited individually, simultaneously or in controlled sequences to produce the required vector thrust for the specific application. During the manufacturing process, a diaphragm is inserted to cover the propellant chamber that is designed to burst at a predetermined plenum pressure. Potential shell coating materials with MEMS applications include ceramic polymer compositions, polymer composites, and polymer metal composites. This design offers high production principals due to its simplicity and very few steps of assembly, which could be made in strip or sheet forms. Further, the array system can be directly integrated into other interceptor parts reducing even further the assembly, manufacturing, cost and weight.
The division of the propellant into small quantities also results in safer operation, and minimizes the possibility of hazard occurring in the storage of the propellant. The propellant powder can be transported safely if it is packed in small containers or capsules and then packaged in each separate cell. These capsules are safe even in external fire and can be extinguished easily by water spray. The absence of toxic products in the decomposition reaction makes the propellant and the resulting microthruster array environmentally friendly in every stage from the manufacturing of the propellant to array assembly.
While this system is conceptually straightforward, there are outstanding questions regarding the performance of the blowing agent as a propellant and the impact of the diaphragm-rupture event on the performance of the thruster. The following sections describe the development of a meso-scale thruster prototype that is used to study the performance of the blowing agent, and the development of a numerical model that was used to study internal flow dynamics after the diaphragm rupture.

3. Proof of Concept Thruster Design

3.1. Design Overview

The design and analysis of an individual thruster is the natural starting point for the characterization and optimization of this propulsion concept, and is the focus of the remainder of this article. A single solid propellant thruster with nozzle exit at ambient pressure was thus designed and tested. Figure 2 shows a schematic diagram of the complete propulsion system with key components labeled.
A ubiquitous component of the thruster design is the converging-diverging supersonic nozzle which produces thrust by efficiently converting the thermal energy of the inlet gases into kinetic energy at the outlet. In this meso-scale prototype, the nozzle was designed with 2 mm throat and an expander half angle of 15 based on design considerations for supersonic micronozzles presented by Louisos et al. [17]. The nozzle geometry, shown in Figure 3, is symmetric about the throat and has an outlet area expansion ratio of 36. Both the propellant chamber and the nozzle were fabricated using additive manufacturing and thermally resistant plastic based filament (i.e., polylactic-acid/PLA) with melting temperature 180–220 C.
In the current thruster design, the propellant decomposition chamber houses a nichrome heating element that holds an aluminum foil packet containing a mass of the CBA powder (see Figure 4a). The heater is designed for maximum surface area coverage of the powder using a coiled geometry. Additionally, hollow cylinders are employed to prevent direct contact between the nichrome wire and the decomposition chamber; such direct contact was found to thermally damage the plastic chamber after a few tests. The electrical leads of the nichrome heating element are passed out of the chamber via ferrule elements to provide external triggering.
The rupture disc is a single-use pressure relief device used to separate nozzle and chamber components until the maximum decomposition pressure condition is achieved. With no moving parts or release times, the rupture disc is a simple, reliable, and fast acting mechanism to release the pressure (Figure 4b,c). In this study, aluminum 8011 alloy was used for its combination of light weight and corrosion-resistance. The burst pressure is adjusted by varying the number of layers; in this study, 4–8 layers were considered.
A key consideration of this design is that of modularity, as the chamber, nozzle, and mounting cart are all capable of being re-configured. Consequently, longer chambers and arbitrary nozzle dimensions can be fabricated according to the application, making the fabrication process effective and efficient. Additionally, for this design, the rupture disc and the igniter can each be located at different positions and in different configurations. For the rupture disc, obvious locations include the entrance to the nozzle or at the nozzle throat. The position of the rupture disc leads to different flow behavior within the nozzle; this is discussed in detail in the numerical simulations described in Section 5. The igniter/propellant assembly can also be packaged in different configurations. The experiments have shown that the configuration impacts the decomposition rate of the propellant and the completeness of the decomposition, and hence the overall microthruster performance.

3.2. Performance Estimation

The propellant used in the experiments was the commercial blowing agent AccuCell-196 provided by KibbeChem Inc., Elkhart, IN, USA; the CBA composition details are proprietary to the manufacturer. This CBA has a gas yield of 218 cm 3 /g and decomposition temperature in the range of 180–205 C according to the manufacturer. AccuCell-196 is a brand of azodicarbonamide, C 2 H 4 N 4 O 2 , which decomposes into a mixture of nitrogen (N 2 ), carbon monoxide (CO) and carbon dioxide (CO 2 ) in a ratio of 65:32:3. The residual solids are made up of urazole, biurea, cyanuric acid, urea, and ammonia salt. Using a mixture-based approach to obtain physical properties of the gas, along with quasi-1D nozzle theory, it is possible to calculate a thrust value for the prototype. It was assumed that the nozzles perfectly expand the flow such that the pressure at the exit is equal to the ambient pressure that each nozzle was designed for. The initial conditions were estimated to be a temperature of 500 K and a pressure of five times that of atmospheric. The throat diameter of 2 mm was conservatively chosen based on the resolution of the 3-D printer used in fabrication. The prototype was designed for eventual operation at rarefied ambient pressures, although initial testing would be at atmospheric pressures (Table 2). Under these idealized design assumptions, a thrust for this meso-scale device operating under steady-state conditions was estimated to be 7.0 N.

4. Experimental Results

4.1. Propellant Characterization

A first step in the CBA propellant characterization was that of obtaining an empirical correlation between the CBA mass and the corresponding adiabatic decomposition temperature. Experiments were thus performed in a high-pressure vessel apparatus (Parr Instruments, Moline, IL, USA) with a volume of 25 mL volume. The vessel configuration is that of a cylinder and a band heater, with the latter connected to a PID controller to increase temperature gradually. Masses between 0.3 and 1.0 ± 0.0001 g of the chemical blowing agent were enclosed inside the vessel and heated up at a constant rate of 15 C per minute until the thermal decomposition was initiated; the heating rate used was based upon the heating element adapted to the test apparatus. Once initiated, the decomposition is self-sustaining and proceeds until completion. The experimental arrangement used is shown in Figure 5.
Figure 6 presents a collection of the temperature histories within the vessel corresponding to different masses of the CBA (azodicarbonamide). The reaction was observed to be exothermic as expected. When a critical temperature is reached, the vessel goes into rapid “thermal runaway” and the chemical reaction produces gas and heat, with the exothermic peak in the range of 98–119 C. For these experiments, one expects the decomposition temperature characteristics to be independent of the propellant mass. This is essentially what is observed, with any variability observed being attributed to the limitations of the testing apparatus and measurement capabilities. Further, the similarity in the trends shows the repeatability of the decomposition process. It is important to note that the extended amount of time required for the decomposition to occur is a direct result of the experimental test chamber used to conduct the testing. The test apparatus contained a small mass of propellant within a relatively large volume that was thermally insulated and the heating mechanism was not optimized for its performance.
According to data reported in the literature [18], the range for the decomposition of pure azodicarbonamide initiates at a significantly higher temperature (∼190 C). This suggests that azodicarbonamide samples used in these experiments were mixed with other compounds to effectively reduce its decomposition temperature and presumably to improve the decomposition performance. Possible additives include metal salts, such as Pb, Zn, and Cd compounds, some of which are used as stabilizers. Therefore, the degree to which the decomposition temperature is reduced depends on which additive is chosen and the amount used [19]. Furthermore, Sims et al. [20] found that the decomposition temperature of chemical blowing agent systems based on azodicarbonamide strongly depends on the powder particle size and the imposed heating rate. This is especially noteworthy when considering possible deleterious impacts upon thrust production due to the presence of small particulates resulting from incomplete decomposition.
The corresponding measurements of the maximum decomposition pressure occurring in the vessel appear in Figure 7. These data confirm a linear relationship between the peak pressure and the mass of CBA. A curious result that was persistent throughout all of the experiments, however, was the fact that the observed linear relation did not intercept through the origin. The experimental resolution of the apparatus prevented accurate measurements for masses below 0.3 g, and so it was impossible to resolve this question. Lastly, from the combined pressure and temperature profiles, it was observed that the peak exotherms at 98–119 C correspond approximately to the maximum rate of gas yield from the chemical decomposition process.

4.2. Test-Firing and Thrust Measurements

To quantify the thrust characteristics of the prototype, a specially designed test rig was constructed. The test rig consists of the thruster assembly mounted a 3D printed cart designed to operate on a linear air track. All experiments were performed with the nozzle exit at ambient pressure. A DC power supply (5 A and 3.5 V) was used to deliver current to the nichrome heating element via electrical contacts extending from the propellant chamber. This power level was chosen to produce the fastest possible heating and decomposition of the propellant. A photograph of the testing apparatus prior to firing appears in Figure 8. During an experiment, electrical contact is made and the nichrome elements rapidly heats up and decomposes the CBA, thus pressurizing the chamber. At the maximum designed chamber pressure, the rupture disc fails and thrust is produced, thereby accelerating the test rig on the linear air track. The finite duration of the thrust production as the chamber depressurizes corresponds to a total impulse (or impulse bit) for the firing.
Experiments were recorded using a high-speed CCD camera to analyze the position of the test rig on the air track as a function of time. From this data, the velocity and acceleration of the cart were obtained and used to determine instantaneous thrust production. Thrust profiles appear in Figure 9 for different rupture disc thicknesses (i.e., burst pressures). Each dataset shown represents the ensemble average over four trials. At t = 0.0 s, the nichrome igniter is energized and the CBA is heated until the decomposition process is initiated. In the testing, there was a lag observed between the igniter activation and the generation of thrust. This lag was typically on the order of 1–2 s, depending on the level of the burst pressure and the efficiency of the heater/CBA packaging. In all cases, the thrust rapidly reaches a peak value and then decays. The peak value of the thrust occurs repeatably at t 0.02 s for all burst pressures; however, the decay process is lengthier at higher burst pressures. Conservatively, the thrust production is completed within approximately 0.2 s in all cases, and much sooner at lower pressures. The presence of thrust oscillations that can become negative is also observed at the highest burst pressures (Figure 9c,d). This curious result suggests the result of significant flow oscillations occurring within the microthruster including backflow at the exit. This phenomenon is further addressed in Section 5.
The impulse bit delivered by the microthruster during its firing is a key parameter in characterizing the propulsion performance. From the terminal velocity of the test rig on the linear air track, the total impulse can be determined from the test rig’s change in momentum. The results are shown in Figure 10. The error bars in the figure are based upon the variations in total impulse over four trials at each burst pressure (number of foil layers). It was found that a nearly linear relation exists between the impulse bit for this microthruster and the burst pressure, with a range of 4 mN·s (4 layers) to 15 mN·s (8 layers).
It is clear that the thrust levels observed experimentally are far below the estimated values in the preliminary design discussed in Section 3. Although numerous simplifying assumptions were made in that analysis, this cannot account for magnitude of the discrepancy. It is clear that some phenomenon is present in the experimental system that was not accurately presented in the design model. To investigate this further, it was deemed appropriate to examine the internal unsteady flow dynamics using numerical simulations.

5. Numerical Simulations of Transient Nozzle Flow

A Computational Fluid Dynamics (CFD) model was constructed to investigate the transient flow structure resulting from disk rupture and the effects on thruster performance. The explicit, finite volume, compressible flow solver rhoCentralFoam [21], part of the OpenFOAM v4.1 CFD suite, was used to conduct fully-transient laminar simulations. The central upwind scheme of Kurganov et al. [22] was used with van Leer-limiting [23].
The computational domain was developed based on the prototype geometry and consists of the plenum, converging-diverging nozzle, and a region downstream that extend 10 times the nozzle exit diameter in the axial and transverse directions. Symmetry is leveraged to reduce the computational expenditure by simulating a 2 wedge of the full axisymmetric nozzle. To model the effect of the rupture disc, mesh elements upstream of the disc are initialized at the rupture pressure of the disc ( P 0 = 544.5 kPa) and the decomposition temperature of the propellant ( T 0 = 400 K), while elements downstream are initialized at ambient conditions ( P = 101.3 kPa and T = 300 K). As a result, a discontinuity (shock) lies between these two regions initially. The blowing agent is assumed to be completely decomposed and the effects of solid residue are not considered. The plenum and nozzle walls are treated as no-slip and adiabatic. A far-field pressure relaxed condition is imposed at the outlet of the extend downstream region with p = 101.3 kPa and T = 300 K.
The computational mesh was generated using Pointwise® grid generation software. The final static mesh consisted of four structured domains with a total of ≈125,000 cells. Octree style dynamic mesh refinement based on the magnitude of the density gradient is implemented for better resolution of the transient shock structure. A value of | | ρ | | = 1000 kg/m 4 is used as a threshold for cell splitting (refinement) and a value of | | ρ | | = 800 kg/m 4 is used for recombination of split cells (coarsening). Cells are allowed to be subdivided a maximum of two times producing 16 sub-cells per original cell in refinement regions. The number of cells in the dynamically refined mesh varied throughout the course of the simulation and primary stayed within the range of 400,000–1,000,000 total cells. Good agreement was found between the time-dependent thrust and mass flow rate estimated by the dynamically refined mesh and static mesh indicating the achievement of grid independence.
Performance results from the transient simulation are presented in Figure 11 with the instantaneous mass-flow rate and thrust plotted in Figure 11a,b, respectively. The mass flow is assessed at the nozzle exit plane according to
m ˙ = A ρ ( U · n ) d A
and the thrust is assessed by applying
F = A ( ( p p ) · n ) d A + A ( τ · n ) d A A U ( ρ U · n ) d A
to the nozzle and plenum walls. During the firing, it is observed that the nozzle experiences periods of both outflow and backflow with corresponding periods of positive and negative thrust. This is consistent with the experimental observations in Figure 9. Plots of the cumulative mass flow exiting the thruster and the total impulse delivered appear in Figure 11c,d.
The oscillatory thrust and mass flow appears to be caused by acoustic resonance in the plenum and converging portion of the nozzle. When the rupture disc breaks, a shock front forms and propagates into the converging portion of the nozzle. The shock then interacts with the converging side-walls and eventually the steep density gradient at the nozzle throat. When the original shock front is reflected off the nozzle throat, it creates a region of rarefaction upstream of the throat. The rarefaction region breaks the choked condition and subsequently draws flow back into the nozzle. Pseudo-Schlieren contours and the axial velocity component are plotted for t = 1.00 ms, the first spike in thrust/mass-flow, and t = 1.50 ms the first trough in thrust/mass-flow in Figure 12. From the Schlieren contours, the complex structure of the flow is evident. At t = 1.00 ms, there are oblique shocks and a Mach-disc forming inside the nozzle, characteristic of the over-expanded condition of the test-firing. A strong shock can also be seen in converging nozzle section. At t = 1.50 ms, the shock structure in the expander is replaced by vortical inflow.
Additional simulations were performed to investigate whether the resonant behavior could be mitigated by moving the rupture disc to the throat from its original upstream position; all other conditions were held constant. The pressure distributions used to initialize the two simulation cases are shown in Figure 13a. A snapshot of the instantaneous velocity field depicting the resultant flow structure of the two rupture disc locations is shown in Figure 13b. When the rupture disk located at the throat a coherent jet structure forms while in the case of the original configuration it does not.
Thrust and impulse profiles for two rupture disc locations were computed and are featured in Figure 14. Simulation results suggest that relocating the rupture disc to the throat yields a smaller amplitude thrust profile but a larger total impulse. This appears to be caused by two phenomena working in tandem. By moving the rupture disc to the throat, the original shock exits out the expander instead of reflecting around the plenum and converging portion of the nozzle. This reduces plenum pressure oscillations and makes the nozzle behavior closer to quasi-stead-state. Moreover, placing the rupture disc at the throat give the nozzle a larger high-pressure reservoir allowing for a longer burn at higher pressures. Together, these two effects yield a smoother thrust profile and larger total impulse.

6. Discussion

This study has demonstrated the viability of a chemical blowing agent as a propellant in a meso-scale prototype thruster, yet it is clear that the system will require further refinement to be efficacious for use in small satellites. The primary challenges lie in: (1) the reduction in size; and (2) the corresponding ignition method, with the latter being manual at present.
In terms of scale reduction, at small scales, it is expected that heat loss through the chamber walls will emerge as a concern for the ignition and decomposition of the propellant. The increased surface area to volume ratio that comes with the small size of the miniaturization can cause chemical and thermal quenching [24,25]. This aspect is further compounded by the relatively slow decomposition process for a blowing agent propellant in comparison to more energetic solid propellants. The increased heat loss has the potential to inhibit or stifle the chemical decomposition, thus negating control over the reaction. The outcome could be a non-sustainable or unstable reaction process leading to non-reproducible results. It is clear that thermal insulation of the reaction chamber will be a key consideration in the scaling down of this design. Future approaches include low thermal conductivity materials such as Foturan (photosensitive glass) as a potential solution in order to minimize the thermal cross-talk between neighboring cells.
The ignition system, in contrast, will likely benefit from the miniaturization of the geometry. With the unique capabilities available in MEMS-manufacturing, it will be possible to better integrate the ignition heater with the blowing agent in the packaging. Examination of the propellant chamber of the current meso-scale prototype has revealed that the current heater/CBA packaging results in unreacted propellant after firing. An illustration of this is shown in Figure 15. With miniaturization, the increase in contact area between the propellant and the heater will make for a more consistent and repeatable ignition process that should also be faster. Successfully tested micro-igniter designs such as polysilicon resistors can be found in the literature [26,27].
A third and final area for the improvement lies in the re-design of the chamber and nozzle geometry based upon the specific and transient nature of the flow in this device. The numerical simulations have shown complex transient shock patterns can form in devices that use a rupture disc. Further, the simulations suggest that the shocks can interact with the nozzle geometry in such a way as to produce a periodic loss of the choked condition at the throat. Without choking present, supersonic flow conditions cannot be achieved and thus the traditional de Laval nozzle design becomes inappropriate. To remedy this, future directions include moving the rupture disc to the throat from the original upstream position and a new design of decomposition chambers with different volumes and/or geometries to promote flow choking. If this proves unworkable, a redesign or complete elimination of the nozzle may be necessary.

7. Conclusions

In this paper, we describe the developmental efforts towards a new miniaturized thruster propulsion concept capable of delivering the low levels of thrust and impulse necessary for attitude and control adjustments inherent to small satellite formation flying. The propulsion system uses the thermally-driven decomposition of a chemical blowing agent as a propellant. The use of the chemical blowing agent is well-suited for small satellite operations, especially CubeSat platforms, where current launch payload restrictions place constraints on storage tank pressurization, amount of stored chemical energy, and toxicity of the propellant.
A meso-scale proof-of-concept thruster using the CBA azodicarbonamide as a propellant was designed and constructed. Experimental testing demonstrated the viability of the concept, but showed that the performance levels were well below those estimated during the design stage. Subsequent numerical simulations of the gas dynamic flow in the device during firing indicated acoustic resonant phenomena in the converging nozzle and plenum attributed to the degradation of performance relative to the idealized steady-state case. The simulations also suggested that substantial performance gains might be achieved by properly designing the converging nozzle section and plenum to attenuate the strong shock that occurs during burst disc rupture.
The transition to a micro-scale device is expected to offer improved performance based primarily on more efficient heating and decomposition of the CBA. Aside from lower power input at the smaller scale, more efficient coupling of the heating element and the CBA is possible; coupled with the greatly reduced thermal mass (volume), much faster activation times are anticipated. Furthermore, studies have shown that the decomposition reaction occurs much more rapidly under pressure. On the micro-scale, it is easier to fabricate a stronger rupture disc than for the meso-scale design. The combination of reduced volume and stronger rupture disc will naturally lead to elevated pressures and improved performance.
Overall, the outcomes of this study indicate that the concept merits continued development. Particular emphases moving forward lie in the reduction of the device size, the ignition/propellant packaging system, and redesign of the plenum and nozzle geometry to improve the transient flow characteristics.

Author Contributions

Propulsion scheme concept, D.L.H., M.R.M., P.C.L. and S.R.-D.; System design, S.R.-D., D.L.H., M.R.M., and P.C.L.; Experimental design, S.R.-D. and L.H.; Experimental testing and analysis, S.R.-D. and L.H.; Computational modeling, J.M.P.; and Manuscript preparation and editing, D.L.H., M.R.M. and J.M.P.

Funding

This research was partially funded by the Vermont Space Grant Consortium under NASA Cooperative Agreement # NNX15AP86H.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Bandyopadhyay, S.; Subramanian, G.P.; Foust, R.; Morgan, D.; Chung, S.J.; Hadaegh, F. A Review of Impending Small Satellite Formation Flying Missions. In Proceedings of the 53rd AIAA Aerospace Sciences Meeting, Kissimmee, FL, USA, 5–9 January 2015. [Google Scholar]
  2. Ruan, W.; Milstein, A.B.; Blackwell, W.; Miller, E.L. A Probabilistic Analysis of Positional Errors on Satellite Remote Sensing Data Using Scattered Interpolation. IEEE Geosci. Remote Sens. Lett. 2017, 14, 261–265. [Google Scholar] [CrossRef]
  3. Rahmat-Samii, Y.; Manohar, V.; Kovitz, J.M. For Satellites, Think Small, Dream Big: A review of recent antenna developments for CubeSats. IEEE Antennas Propag. Mag. 2017, 59, 22–30. [Google Scholar] [CrossRef]
  4. Clark, P.; Malphrus, B.; Reuter, D.; MacDowall, R.; Folta, D.; Hurford, T.; Brambora, C.; Farrell, W. The First Deep Space Cubesat Broadband IR Spectrometer, Lunarcubes, and the Search for Lunar Volatiles. In NASA Technical Reports Server; NASA Goddard Space Flight Center: Greenbelt, MD, USA, 2017. [Google Scholar]
  5. NASA Ames Research Center. Small Spacecraft Technology State of the Art; Technical Report; Mission Design Division, TP-2015-216648; NASA Ames Research Center: Mountain View, CA, USA, 2015.
  6. Reichbach, J.; Sedwick, R.; Martinez-Sanchez, M. Micropropulsion system selection for precision formation flying satellites. In Proceedings of the 37th Joint Propulsion Conference and Exhibit, Salt Lake City, UT, USA, 8–11 July 2001; p. 3646. [Google Scholar]
  7. Cervone, A.; Deeb, A.; van Wees, T.; Jansen, E.; Sundaramoorthy, P.; Chu, J.; Zandbergen, B. A micro-propulsion subsystem to enable formation flying on the DelFFi mission. In Proceedings of the 8th International Workshop on Satellite Constellations and Formation Flying, Delft, The Netherlands, 8–10 June 2015. [Google Scholar]
  8. Persson, S.; Bodin, P.; Gill, E.; Harr, J.; Jörgensen, J. PRISMA-an autonomous formation flying mission. In Proceedings of the ESA Small Satellite Systems and Services Symposium (4S), Sardinia, Italy, 25–29 September 2006. [Google Scholar]
  9. Lee, S.; Hutputanasin, A.; Toorian, A.; Lan, W.; Munakata, R. CubeSat design specification. CubeSat Program 2009, 8651, 22. [Google Scholar]
  10. Stratigentas, D. Small Satellites: Regulatory Requirements and Challenges. Ph.D. Thesis, Martinus Nijhoff Publishers, Leyden, The Netherlands, 2015. [Google Scholar]
  11. Mazouffre, S. Electric propulsion for satellites and spacecraft: Established technologies and novel approaches. Plasma Sources Sci. Technol. 2016, 25, 033002. [Google Scholar] [CrossRef]
  12. Lemmer, K. Propulsion for cubesats. Acta Astronaut. 2017, 134, 231–243. [Google Scholar] [CrossRef]
  13. Levchenko, I.; Xu, S.; Teel, G.; Mariotti, D.; Walker, M.; Keidar, M. Recent progress and perspectives of space electric propulsion systems based on smart nanomaterials. Nat. Commun. 2018, 9, 879. [Google Scholar] [CrossRef] [PubMed]
  14. Levchenko, I.; Bazaka, K.; Ding, Y.; Raitses, Y.; Mazouffre, S.; Henning, T.; Klar, P.J.; Shinohara, S.; Schein, J.; Garrigues, L.; et al. Space micropropulsion systems for Cubesats and small satellites: From proximate targets to furthermost frontiers. Appl. Phys. Rev. 2018, 5, 011104. [Google Scholar] [CrossRef] [Green Version]
  15. Zhang, Y.; Kontopoulou, M.; Ansari, M.; Hatzikiriakos, S.; Park, C.B. Effect of molecular structure and rheology on the compression foam molding of ethylene-α-olefin copolymers. Polym. Eng. Sci. 2011, 51, 1145–1154. [Google Scholar] [CrossRef]
  16. Reglero Ruiz, J.A.; Vincent, M.; Agassant, J.F.; Sadik, T.; Pillon, C.; Carrot, C. Polymer foaming with chemical blowing agents: Experiment and modeling. Polym. Eng. Sci. 2015, 55, 2018–2029. [Google Scholar] [CrossRef] [Green Version]
  17. Louisos, W.F.; Alexeenko, A.A.; Hitt, D.L.; Zilic, A. Design considerations for supersonic micronozzles. Int. J. Manuf. Res. 2008, 3, 80–113. [Google Scholar] [CrossRef]
  18. Bhatti, A.S.; Dollimore, D.; Goddard, R.; O’Donnell, G. The effects of additives on the thermal decomposition of azodicarbonamide. Thermochim. Acta 1984, 76, 273–286. [Google Scholar] [CrossRef]
  19. Bhatti, A.S.; Dollimore, D.; Goddard, R.; O’Donnell, G. The thermal decomposition of azodicarbonamide. Thermochim. Acta 1984, 76, 63–77. [Google Scholar] [CrossRef]
  20. Sims, G.; Jaafar, H. A chemical blowing agent system (CBAS) based on azodicarbonamide. J. Cell. Plast. 1994, 30, 175–188. [Google Scholar] [CrossRef]
  21. Greenshields, C.; Weller, H.; Gasparini, L.; Reese, J. Implementaiton of semi-discrete, non-staggered central schemes in a colocated, polyhedral, finite volume framework, for high-speed viscous flows. Int. J. Numer. Methods Fluids 2010, 63, 1–21. [Google Scholar]
  22. Kurganov, A.; Noelle, S.; Petrova, G. Semi-discrete Central-Upwind Schemes for Hyperbolic Conservation Laws and Hamilton-Jacobi Equations. SIAM J. Sci. Comput. 2001, 23, 707–740. [Google Scholar] [CrossRef]
  23. van Leer, B. Towards the ultimate conservative difference scheme, II: Monotonicity and conservation combined in second order scheme. J. Comput. Phys. 1974, 17, 361–370. [Google Scholar] [CrossRef]
  24. Zhang, K.; Chou, S.; Ang, S.S. Development of a solid propellant microthruster with chamber and nozzle etched on a wafer surface. J. Micromech. Microeng. 2004, 14, 785. [Google Scholar] [CrossRef]
  25. Lee, D.H.; Kwon, S. Heat transfer and quenching analysis of combustion in a micro combustion vessel. J. Micromech. Microeng. 2002, 12, 670. [Google Scholar] [CrossRef]
  26. Chaalane, A.; Chemam, R.; Houabes, M.; Yahiaoui, R.; Metatla, A.; Ouari, B.; Metatla, N.; Mahi, D.; Dkhissi, A.; Esteve, D. A MEMS-based solid propellant microthruster array for space and military applications. J. Phys. Conf. Ser. 2015, 660, 012137. [Google Scholar] [CrossRef] [Green Version]
  27. Wu, X.; Dong, P.; Li, Z.; Li, S.; Liu, Q.; Xu, C.; Wan, H. Design, fabrication and characterization of a solid propellant micro-thruster. In Proceedings of the 2009 4th IEEE International Conference on Nano/Micro Engineered and Molecular Systems, Shenzhen, China, 5–8 January 2009; pp. 476–479. [Google Scholar]
Figure 1. Schematic of the conceptual microthruster arrays mounted on a 3U cubesat. The cubesat dimensions are 10 × 10 × 30 cm and the envisioned size of the discrete microthruster array is less than 1 cm.
Figure 1. Schematic of the conceptual microthruster arrays mounted on a 3U cubesat. The cubesat dimensions are 10 × 10 × 30 cm and the envisioned size of the discrete microthruster array is less than 1 cm.
Aerospace 05 00119 g001
Figure 2. The basic thruster design consists of: (1) decomposition chamber; (2) heating element; (3) rupture disc; and (4) o-ring ( × 2 ). The converging/diverging nozzle has been outlined for clarity.
Figure 2. The basic thruster design consists of: (1) decomposition chamber; (2) heating element; (3) rupture disc; and (4) o-ring ( × 2 ). The converging/diverging nozzle has been outlined for clarity.
Aerospace 05 00119 g002
Figure 3. Mechanical drawing showing side view of the thruster design geometry. Dimensions are in mm.
Figure 3. Mechanical drawing showing side view of the thruster design geometry. Dimensions are in mm.
Aerospace 05 00119 g003
Figure 4. (a) Inside the chamber, the connected heating element holding the aluminum foil cup before enclosure; (b) the aluminum 8011 alloy rupture disc before firing; and (c) the aluminum 8011 alloy rupture disc after firing.
Figure 4. (a) Inside the chamber, the connected heating element holding the aluminum foil cup before enclosure; (b) the aluminum 8011 alloy rupture disc before firing; and (c) the aluminum 8011 alloy rupture disc after firing.
Aerospace 05 00119 g004
Figure 5. Schematic diagram of pressure vessel apparatus.
Figure 5. Schematic diagram of pressure vessel apparatus.
Aerospace 05 00119 g005
Figure 6. Temperature histories for pressure vessel experiments for different masses of the CBA. The onset of the decomposition of the CBA at the activation temperature is accompanied by a sudden increase in temperature and pressure. It is seen that the temperature histories are mostly insensitive to the initial CBA mass, as would be expected.
Figure 6. Temperature histories for pressure vessel experiments for different masses of the CBA. The onset of the decomposition of the CBA at the activation temperature is accompanied by a sudden increase in temperature and pressure. It is seen that the temperature histories are mostly insensitive to the initial CBA mass, as would be expected.
Aerospace 05 00119 g006
Figure 7. Maximum decomposition pressure, as a function of mass of blowing agent. Data taken from the pressure vessel test. The straight line represents the best fit for the testing samples with mass less than 1 g, with a coefficient of determination R 2 = 0.99 .
Figure 7. Maximum decomposition pressure, as a function of mass of blowing agent. Data taken from the pressure vessel test. The straight line represents the best fit for the testing samples with mass less than 1 g, with a coefficient of determination R 2 = 0.99 .
Aerospace 05 00119 g007
Figure 8. A photograph of the prototype thruster in testing configuration on an air track.
Figure 8. A photograph of the prototype thruster in testing configuration on an air track.
Aerospace 05 00119 g008
Figure 9. Thrust profiles obtained from cart velocity measurements after the disc rupture for 5–8 aluminum layers, corresponding to different burst pressures. The aluminum layers have an average thickness of 0.016 mm. Each curve represents the ensemble average of four trials. The results show there is a sharp spike in thrust followed by a rapid decrease thereafter. Note also the presence of negative thrust values; this is addressed in Section 5.
Figure 9. Thrust profiles obtained from cart velocity measurements after the disc rupture for 5–8 aluminum layers, corresponding to different burst pressures. The aluminum layers have an average thickness of 0.016 mm. Each curve represents the ensemble average of four trials. The results show there is a sharp spike in thrust followed by a rapid decrease thereafter. Note also the presence of negative thrust values; this is addressed in Section 5.
Aerospace 05 00119 g009
Figure 10. Average total impulse from four test trials as a function of the rupture of the number of aluminum layers (i.e., burst pressure). The aluminum layers have an average thickness of 0.016 mm.
Figure 10. Average total impulse from four test trials as a function of the rupture of the number of aluminum layers (i.e., burst pressure). The aluminum layers have an average thickness of 0.016 mm.
Aerospace 05 00119 g010
Figure 11. Simulated nozzle performance as a function of time.
Figure 11. Simulated nozzle performance as a function of time.
Aerospace 05 00119 g011
Figure 12. Simulated flow state at the first maxima in thrust (t = 1.00 ms) and the first minima (t = 1.50 ms).
Figure 12. Simulated flow state at the first maxima in thrust (t = 1.00 ms) and the first minima (t = 1.50 ms).
Aerospace 05 00119 g012
Figure 13. Simulation results for the flow field for two different rupture disc locations.
Figure 13. Simulation results for the flow field for two different rupture disc locations.
Aerospace 05 00119 g013
Figure 14. Simulated nozzle performance as a function of time.
Figure 14. Simulated nozzle performance as a function of time.
Aerospace 05 00119 g014aAerospace 05 00119 g014b
Figure 15. Photograph inside the propellant chamber after test firing showing unreacted CBA residue. This test was with a rupture disc consisting of eight layers. The heater coil is visible in the center as is the aluminum foil packaging of the CBA, surrounding the coil.
Figure 15. Photograph inside the propellant chamber after test firing showing unreacted CBA residue. This test was with a rupture disc consisting of eight layers. The heater coil is visible in the center as is the aluminum foil packaging of the CBA, surrounding the coil.
Aerospace 05 00119 g015
Table 1. Target Operating Characteristics.
Table 1. Target Operating Characteristics.
System ParameterTarget Value or Range
Thrust level1–1000 μ N
Impulse bit0.1–100 μ N·s
Mass<0.01 kg
Power consumption<1 W
Table 2. Nozzle operating parameters and thrust prediction for the prototype under assumed steady-state conditions.
Table 2. Nozzle operating parameters and thrust prediction for the prototype under assumed steady-state conditions.
Nozzle ParameterValue
Ambient Pressure (kPa)0.5445
Exit Area Ratio A e / A * 36.0
Exit Mach number5.5
Thrust Estimate (N)7.0

Share and Cite

MDPI and ACS Style

Romero-Diez, S.; Hantsche, L.; Pearl, J.M.; Hitt, D.L.; McDevitt, M.R.; Lee, P.C. A Single-Use Microthruster Concept for Small Satellite Attitude Control in Formation-Flying Applications. Aerospace 2018, 5, 119. https://doi.org/10.3390/aerospace5040119

AMA Style

Romero-Diez S, Hantsche L, Pearl JM, Hitt DL, McDevitt MR, Lee PC. A Single-Use Microthruster Concept for Small Satellite Attitude Control in Formation-Flying Applications. Aerospace. 2018; 5(4):119. https://doi.org/10.3390/aerospace5040119

Chicago/Turabian Style

Romero-Diez, Sandra, Lydia Hantsche, Jason M. Pearl, Darren L. Hitt, M. Ryan McDevitt, and Patrick C. Lee. 2018. "A Single-Use Microthruster Concept for Small Satellite Attitude Control in Formation-Flying Applications" Aerospace 5, no. 4: 119. https://doi.org/10.3390/aerospace5040119

APA Style

Romero-Diez, S., Hantsche, L., Pearl, J. M., Hitt, D. L., McDevitt, M. R., & Lee, P. C. (2018). A Single-Use Microthruster Concept for Small Satellite Attitude Control in Formation-Flying Applications. Aerospace, 5(4), 119. https://doi.org/10.3390/aerospace5040119

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop