Characterization of the Impingement Dynamics of Pulsed Rocket Plumes with the Ground at Low Ambient Pressure

Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/76629/1/AIAA-2007-5707-681.pd

Micro-Nozzle Simulation and Test for an Electrothermal Plasma Thruster

With an increased demand in Cube Satellite (CubeSat) development for low cost science and exploration missions, a push for the development of micro-propulsion technology has emerged, which seeks to increase CubeSat capabilities for novel mission concepts. One type of micro-propulsion system currently under development, known as Pocket Rocket, is an electrothermal plasma micro-thruster. Pocket Rocket uses a capacitively coupled plasma, generated by radio-frequency, in order to provide neutral gas heating via ion-neutral collisions within a gas discharge tube. When compared to a cold-gas thruster of similar size, this gas heating mechanism allows Pocket Rocket to increase the exit thermal velocity of its gaseous propellant for increased thrust. Previous experimental work has only investigated use of the gas discharge tube\u27s orifice for propellant expansion into vacuum. This thesis aims to answer if Pocket Rocket may see an increase in thrust with the addition of a micro-nozzle, placed at the end of the gas discharge tube. With the addition of a conical ε = 10, α = 30° micro-nozzle, performance increases of up to 6% during plasma operation, and 25% during cold gas operation, have been observed. Propellant heating has also been observed to increase by up to 60 K within the gas discharge tube

Underexpanded Supersonic Plume Surface Interactions: Applications for Spacecraft Landings on Planetary Bodies

Numerical and experimental investigations of both far-field and near-field supersonic steady jet interactions with a flat surface at various atmospheric pressures are presented in this paper. These studies were done in assessing the landing hazards of both the NASA Mars Science Laboratory and Phoenix Mars spacecrafts. Temporal and spatial ground pressure measurements in conjunction with numerical solutions at altitudes of approx.35 nozzle exit diameters and jet expansion ratios (e) between 0.02 and 100 are used. Data from steady nitrogen jets are compared to both pulsed jets and rocket exhaust plumes at Mach approx.5. Due to engine cycling, overpressures and the plate shock dynamics are different between pulsed and steady supersonic impinging jets. In contrast to highly over-expanded (e 5 (lunar atmospheric regime), the ground pressure is minimal due to the development of a highly expansive shock structure. We show this is dependent on the stability of the plate shock, the length of the supersonic core and plume decay due to shear layer instability which are all a function of the jet expansion ratio. Asymmetry and large gradients in the spatial ground pressure profile and large transient overpressures are predominantly linked to the dynamics of the plate shock. More importantly, this study shows that thruster plumes exhausting into martian environments possess the largest surface pressure loads and can occur at high spacecraft altitudes in contrast to the jet interactions at terrestrial and lunar atmospheres. Theoretical and analytical results also show that subscale supersonic cold gas jets adequately simulate the flow field and loads due to rocket plume impingement provided important scaling parameters are in agreement. These studies indicate the critical importance of testing and modeling plume-surface interactions for descent and ascent of spacecraft and launch vehicles

Transient Three-Dimensional Side Load Analysis of Out-of-Round Film Cooled Nozzles

The objective of this study is to investigate the effect of nozzle out-of-roundness on the transient startup side loads. The out-of-roundness could be the result of asymmetric loads induced by hardware attached to the nozzle, asymmetric internal stresses induced by previous tests and/or deformation, such as creep, from previous tests. The rocket engine studied encompasses a regeneratively cooled thrust chamber and a film cooled nozzle extension with film coolant distributed from a turbine exhaust manifold. The computational methodology is based on an unstructured-grid, pressure-based computational fluid dynamics formulation, and a transient inlet history based on an engine system simulation. Transient startup computations were performed with the out-of-roundness achieved by four degrees of ovalization of the nozzle: one perfectly round, one slightly out-of-round, one more out-of-round, and one significantly out-of-round. The computed side load physics caused by the nozzle out-of-roundness and its effect on nozzle side load are reported and discussed

Transient Three-Dimensional Startup Side Load Analysis of a Regeneratively Cooled Nozzle

The objective of this effort is to develop a computational methodology to capture the startup side load physics and to anchor the computed aerodynamic side loads with the available data from a regeneratively cooled, high-aspect-ratio nozzle, hot-fired at sea level. The computational methodology is based on an unstructured-grid, pressure-based, reacting flow computational fluid dynamics and heat transfer formulation, a transient 5 s inlet history based on an engine system simulation, and a wall temperature distribution to reflect the effect of regenerative cooling. To understand the effect of regenerative wall cooling, two transient computations were performed using the boundary conditions of adiabatic and cooled walls, respectively. The results show that three types of shock evolution are responsible for side loads: generation of combustion wave; transitions among free-shock separation, restricted-shock separation, and simultaneous free-shock and restricted shock separations; along with the pulsation of shocks across the lip, although the combustion wave is commonly eliminated with the sparklers during actual test. The test measured two side load events: a secondary and lower side load, followed by a primary and peak side load. Results from both wall boundary conditions captured the free-shock separation to restricted-shock separation transition with computed side loads matching the measured secondary side load. For the primary side load, the cooled wall transient produced restricted-shock pulsation across the nozzle lip with peak side load matching that of the test, while the adiabatic wall transient captured shock transitions and free-shock pulsation across the lip with computed peak side load 50% lower than that of the measurement. The computed dominant pulsation frequency of the cooled wall nozzle agrees with that of a separate test, while that of the adiabatic wall nozzle is more than 50% lower than that of the measurement. The computed teepee-like formation and the tangential motion of the shocks during lip pulsation also qualitatively agree with those of test observations. Moreover, a third transient computation was performed with a proportionately shortened 1 s sequence, and lower side loads were obtained with the higher ramp rate

Navier/Stokes/Direct Simulation Monte Carlo Modeling of Small Cold Gas Thruster Nozzle and Plume Flows

This study involves the modeling of small cold-gas (N2) thrusters nozzle and plume flows, their interactions with spacecraft surfaces and the induced pressure environment. These small cold-gas thrusters were used for pitch, yaw and roll control and were mounted on the bottom of the conical Environmental Monitor Payload (EMP) suborbital spacecraft. The pitch and yaw thrusters had 0.906 mm throat diameter and 4.826 mm exit diameter, while the roll thrusters had 1.6 mm throat diameter and 5.882 mm exit diameter. During thruster firing, at altitudes between 670 km and 1200 km, pressure measurements exhibited non-periodic pulses (Gatsonis et al., 1999). The pressure sensor was located inside the EMP and was connected to it\u27s sidewall with a 0.1-m long, 0.022-m diameter tube and the pressure pulses appeared instantaneously with the firings for thrusters without a direct line-of-sight with the sensor entrance. Preliminary analysis showed that the plume of these small EMP thrusters undergoes transition from continuous to rarefied. Therefore, nozzle and plume simulations are performed using a combination of Navier-Stokes and Direct Simulation Monte Carlo codes. This study presents first a validation of the Navier-Stokes code Rampant used for the continuous EMP nozzle and plume simulations. The first Rampant validation example involves a two-dimensional axisymetric freejet expansion and is used to demonstrate the use of Bird\u27s breakdown parameter. Results are compared favorably with those of Bird (1980) obtained through the method of characteristics. The second validation example involves three-dimensional plume simulations of a NASA thruster. This nitrogen nozzle has a throat diameter of 3.18 mm, an exit diameter of 31.8 mm, half-angle of 20 degrees, stagnation temperature of 699 K, stagnation pressure of 6,400 Pa. Simulation results are compared favorably with previous Navier-Stokes and Direct Simulation Monte Carlo numerical work. The third validation example involves three-dimensional simulations of Rothe\u27s (1970) nozzle that has a throat diameter of 2.5 mm, an exit diameter of 20.3 mm, half-angle of 20 degrees, operating at stagnation temperature of 300 K and pressure of 1975 Pa. Numerical results also compared favorably to experimental data. The combined Navier-Stokes/DSMC approach and the EMP simulation results are presented and discussed. The continuous part of the EMP nozzle and plume flow is modeled using the three-dimensional Navier-Stokes Rampant code. The Navier-Stokes domain includes the geometry of the nozzle and the EMP base until transition of the continuous flow established by Bird\u27s breakdown parameter. The rarefied part of the plume flow is modeled using the Direct Simulation Monte Carlo code DAC. Flowfield data obtained inside the breakdown surface from the Navier-Stokes simulation are used as inputs to the DSMC simulations. The DSMC domain includes the input surface and the EMP spacecraft geometry. The combined Navier-Stokes/DSMC simulations show the complex structure of the plume flow as it expands over the EMP surfaces. Plume reflection and backflow are demonstrated. The study also summarizes findings presented by Gatsonis et al. (2000), where the DSMC predictions at the entrance of the pressure sensor are used as inputs to a semi-analytical model to predict the pressure inside the sensor. It is shown that the pressure predictions for the pitch/yaw thrusters are close to the measurements. The plume of a pitch or yaw thruster reaches the pressure sensor after expanding on the EMP base. The pressure predicted for the roll thruster is larger that the measured. This is attributed to the uncertainty in the roll thruster location on the EMP base resulting, in the simulation, in a component of direct flow to the sensor

Supersonic Constricted Plasma Flows

The Pocket Rocket electrothermal microthruster is a miniaturised electric propulsion system designed for nanosatellites operating in space. A weakly ionised capacitively coupled plasma is ignited in the flowing Ar gas propellant within a constricted discharge chamber at 1 Torr using less than 10 W of radiofrequency power. The discharge can operate either continuously or in rapid pulsed mode since plasma breakdown initiates almost instantaneously on a μs time scale. The propellant is heated to temperatures approaching 1000 K and is expanded through a converging-diverging nozzle into vacuum at supersonic velocities. Thrust on the order of 1 mN is generated as a reactionary force to the linear momentum of the expelled neutral gas propellant. This thesis presents a comprehensive model of Pocket Rocket developed with computational fluid dynamics and plasma simulations. Boundary layer effects are significant in the rarefied flow within the constricted discharge chamber. A slip boundary condition with the appropriate tangential momentum and thermal accommodation coefficients must be used to produce results that precisely match experimental measurements. The problem of including vacuum regions within a fluid simulation domain is unconventionally circumvented by taking advantage of the flow velocity choking. The computed sonic surface, thrust force, and specific impulse are in good agreement with theoretical predictions. Volumetric plasma-induced heating of the background neutral gas is primarily due to ion-neutral charge exchange collisions, with very little contribution from electron-neutral elastic collisions. The propellant temperature is described by two local models based on the different ion transport behaviour in the plasma bulk and plasma sheath. The most dominant process is surface bombardment by ions accelerated through the plasma sheath, which heats the discharge chamber wall and is responsible for the creation of secondary electrons that sustain the gamma mode discharge. The geometrical area asymmetry of the grounded and powered electrodes results in a self-bias that manifests as a spatially nonuniform negative charging within the dielectric discharge chamber wall. In the thin sheath regime, the self-biased waveform has a diminished trailing edge at each positive peak, and asymmetrically displaced negative peaks due to the extraneous impedance of the dielectric wall. This leads to a redefinition of the self-bias voltage that uses the maxima envelope of the self-biased waveform instead of the mean, which maintains consistency with different extraneous impedances. The performance of Pocket Rocket is improved by optimising the physical and electrical geometry for thrust and boundary layer effects, and plasma confinement is achieved through the formation of a conical plasma sheath at the nozzle throat. Enhanced recombination in the supersonic expanding plume creates a neutral exhaust, thereby avoiding contamination of externally mounted solar panels and interference with sensitive instruments. Most importantly, the combination of flow velocity choking and plasma confinement results in a convergent plasma simulation that accurately models plasma expansion into vacuum. The computational fluid dynamics and plasma modelling technique and analysis presented in this thesis are not restricted to the Pocket Rocket discharge and may be adapted for other discharges at different pressure regimes and physical scales

Bibliography of Lewis Research Center technical publications announced in 1992

This compilation of abstracts describes and indexes the technical reporting that resulted from the scientific and engineering work performed and managed by the Lewis Research Center in 1992. All the publications were announced in the 1992 issues of STAR (Scientific and Technical Aerospace Reports) and/or IAA (International Aerospace Abstracts). Included are research reports, journal articles, conference presentations, patents and patent applications, and theses

Thirteenth Workshop for Computational Fluid Dynamic Applications in Rocket Propulsion and Launch Vehicle Technology

This conference publication includes various abstracts and presentations given at the 13th Workshop for Computational Fluid Dynamic Applications in Rocket Propulsion and Launch Vehicle Technology held at the George C. Marshall Space Flight Center April 25-27 1995. The purpose of the workshop was to discuss experimental and computational fluid dynamic activities in rocket propulsion and launch vehicles. The workshop was an open meeting for government, industry, and academia. A broad number of topics were discussed including computational fluid dynamic methodology, liquid and solid rocket propulsion, turbomachinery, combustion, heat transfer, and grid generation

Experimental and analytical comparison of flowfields in a 110 N (25 lbf) H2/O2 rocket

A gaseous hydrogen/gaseous oxygen 110 N (25 lbf) rocket was examined through the RPLUS code using the full Navier-Stokes equations with finite rate chemistry. Performance tests were conducted on the rocket in an altitude test facility. Preliminary parametric analyses were performed for a range of mixture ratios and fuel film cooling pcts. It is shown that the computed values of specific impulse and characteristic exhaust velocity follow the trend of the experimental data. Specific impulse computed by the code is lower than the comparable test values by about two to three percent. The computed characteristic exhaust velocity values are lower than the comparable test values by three to four pct. Thrust coefficients computed by the code are found to be within two pct. of the measured values. It is concluded that the discrepancy between computed and experimental performance values could not be attributed to experimental uncertainty