Search NASA⌕ Search

SEARCH · Search NASA

Results for “rocket equation”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 235 records · Page 13

Hypergolic bipropellant spray combustion and flow modelling in rocket engines

A predictive tool for hypergolic bipropellant spray combustion and flow evolution in small rocket combustion chambers is described. It encompasses a computational technique for the gas-phase governing equations, a discrete particle method for liquid bipropellant sprays, and constitutive models for combustion chemistry, interphase exchanges, and unlike impinging hypergolic spray interactions. Emphasis is placed on the phenomenological modeling of the hypergolic liquid bipropellant gasification processes. Sample computations with the N2H4-N2O4 propellant system are given in order to show some of the capabilities and inadequacies of this tool.

Larosiliere, Louis M.↗

Similarity solution of the boundary-layer equations for laser heated flows

The laminar boundary-layer equations with local similarity approximation are solved in order to estimate the heat transfer rate to the wall in a laser-heated rocket thruster using pure hydrogen, where the average temperature of the hot plasma core is about 14,000 K and the core Reynolds number based on nozzle throat diameter is about 2000. Under these conditions, the density-viscosity product at the wall can be 10 times the freestream value, and the hydrogen is completely dissociated. Hence the boundary layer equations with variable transport properties are solved by a quasi-linearization technique, and the equilibrium properties of hydrogen are used in the calculations. Velocity profiles are obtained, and the wall shear and wall stagnation enthalpy gradient are plotted against the pressure gradient parameter.

Wu, P. K.↗

Solid rocket combustion simulator technology using the hybrid rocket for simulation

The hybrid rocket is reexamined in light of several important unanswered questions regarding its performance. The well-known heat transfer limited burning rate equation is quoted, and its limitations are pointed out. Several inconsistencies in the burning rate determination through fuel depolymerization are explicitly discussed. The resolution appears to be through the postulate of (surface) oxidative degradation of the fuel. Experiments are initiated to study the fuel degradation in mixtures of nitrogen/oxygen in the 99.9 percent/0.1 percent to 98 percent/2 percent range. The overall hybrid combustion behavior is studied in a 2 in-diameter rocket motor, where a PMMA tube is used as the fuel. The results include detailed, real-time infrared video images of the combustion zone. Space- and time-averaged images give a broad indication of the temperature reached in the gases. A brief outline is shown of future work, which will specifically concentrate on the exploration of the role of the oxidizer transport to the fuel surface, and the role of the unburned fuel that is reported to escape below the classical time-averaged boundary layer flame.

Ramohalli, Kumar↗

A Computational and Experimental Investigation of Shear Coaxial Jet Atomization

The instability and subsequent atomization of a viscous liquid jet emanated into a high-pressure gaseous surrounding is studied both computationally and experimentally. Liquid water issued into nitrogen gas at elevated pressures is used to simulate the flow conditions in a coaxial shear injector element relevant to liquid propellant rocket engines. The theoretical analysis is based on a simplified mathematical formulation of the continuity and momentum equations in their conservative form. Numerical solutions of the governing equations subject to appropriate initial and boundary conditions are obtained via a robust finite difference scheme. The computations yield real-time evolution and subsequent breakup characteristics of the liquid jet. The experimental investigation utilizes a digital imaging technique to measure resultant drop sizes. Data were collected for liquid Reynolds number between 2,500 and 25,000, aerodynamic Weber number range of 50-500 and ambient gas pressures from 150 to 1200 psia. Comparison of the model predictions and experimental data for drop sizes at gas pressures of 150 and 300 psia reveal satisfactory agreement particularly for lower values of investigated Weber number. The present model is intended as a component of a practical tool to facilitate design and optimization of coaxial shear atomizers.

Ibrahim, Essam A.↗

Investigation of Navier-Stokes Code Verification and Design Optimization

With rapid progress made in employing computational techniques for various complex Navier-Stokes fluid flow problems, design optimization problems traditionally based on empirical formulations and experiments are now being addressed with the aid of computational fluid dynamics (CFD). To be able to carry out an effective CFD-based optimization study, it is essential that the uncertainty and appropriate confidence limits of the CFD solutions be quantified over the chosen design space. The present dissertation investigates the issues related to code verification, surrogate model-based optimization and sensitivity evaluation. For Navier-Stokes (NS) CFD code verification a least square extrapolation (LSE) method is assessed. This method projects numerically computed NS solutions from multiple, coarser base grids onto a freer grid and improves solution accuracy by minimizing the residual of the discretized NS equations over the projected grid. In this dissertation, the finite volume (FV) formulation is focused on. The interplay between the xi concepts and the outcome of LSE, and the effects of solution gradients and singularities, nonlinear physics, and coupling of flow variables on the effectiveness of LSE are investigated. A CFD-based design optimization of a single element liquid rocket injector is conducted with surrogate models developed using response surface methodology (RSM) based on CFD solutions. The computational model consists of the NS equations, finite rate chemistry, and the k-6 turbulence closure. With the aid of these surrogate models, sensitivity and trade-off analyses are carried out for the injector design whose geometry (hydrogen flow angle, hydrogen and oxygen flow areas and oxygen post tip thickness) is optimized to attain desirable goals in performance (combustion length) and life/survivability (the maximum temperatures on the oxidizer post tip and injector face and a combustion chamber wall temperature). A preliminary multi-objective optimization study is carried out using a geometric mean approach. Following this, sensitivity analyses with the aid of variance-based non-parametric approach and partial correlation coefficients are conducted using data available from surrogate models of the objectives and the multi-objective optima to identify the contribution of the design variables to the objective variability and to analyze the variability of the design variables and the objectives. In summary the present dissertation offers insight into an improved coarse to fine grid extrapolation technique for Navier-Stokes computations and also suggests tools for a designer to conduct design optimization study and related sensitivity analyses for a given design problem.

Vaidyanathan, Rajkumar↗

Numerical Investigation of Slag Behavior for RSRM

It is known that the flow field of the redesigned solid rocket motor (RSRM) is very complicated due to the complex characteristics of turbulent multi-phase flow, chemical reaction, particle combustion, evaporation, breakup and agglomeration etc. It requires multi-phase calculations, chemical reaction simulation, and particle combustion, evaporation, and breakup models to obtain a better understanding of thermophysics for the RSRM design using numerical methods. Also, the slag buildup due to the molten particles is another factor affecting the performance of the RSRM. To achieve this goal, the volume of fluid (VOF) method is used to capture the free surface motion so as to simulate the accumulation of the molten particles (slag) of the RSRM. A finite rate chemistry model is used to simulate the chemical reaction effects. For multi-phase calculations, the Hermsen combustion model is used for the aluminum particle combustion analysis and the Taylor Analogy Breakup (TAB) model is used for the particle breakup analysis. An interphase mas-exchange model introduced by Spalding is used for the evaporation calculation. The particle trajectories are calculated using a one-step implicit method for several groups of particle sizes by which the drag forces and heat fluxes are then coupled with the gas phase equations. The preliminary results predicted a reasonable physical simulation of the particle effects using a simple two dimensional solid rocket motor configuration. It shows that the AL/AL2O3 particle sizes are reduced due to the combustion, evaporation, and breakup. The flow field is disturbed by the particles. Mach number distributions in the nozzle are deformed due to the effect of particle concentrations away from the center line.

Liaw, P.↗

A theoretical evaluation of laser-sustained plasma thruster performance

An extensive numerical experiment has been conducted to evaluate rocket thruster performance using a laser-sustained hydrogen plasma as the propellant. The plasma was sustained using a 30 kW CO2 laser beam operated at 10.6 microns focused inside the thruster. The steady-state Navier-Stokes equations coupled with the laser power absorption process have been solved numerically. A pressure based Navier-Stokes solver using body-fitted coordinate was used to calculate the laser-supported rocket flow which included both recirculating and transonic flow regions. The local thermodynamic equilibrium (LTE) assumption was used for the plasma thermophysical and optical properties. Geometric ray tracing was adopted to describe the laser beam. Several different throat size thrusters operated at 150 and 300 kPa chamber stagnation pressure were studied. It was found that the thruster performance (vacuum specific impulse) was highly dependent on the operating conditions, and a properly designed laser supported thruster can attain a specific impulse around 1500 secs. The heat loading on the thruster wall was also estimated and was in the range of that for a conventional chemical rocket.

Jeng, San-Mou↗

A Generalized Eulerian-Lagrangian Analysis, with Application to Liquid Flows with Vapor Bubbles

Under a NASA MSFC SBIR Phase 2 effort an analysis has been developed for liquid flows with vapor bubbles such as those in liquid rocket engine components. The analysis is based on a combined Eulerian-Lagrangian technique, in which Eulerian conservation equations are solved for the liquid phase, while Lagrangian equations of motion are integrated in computational coordinates for the vapor phase. The novel aspect of the Lagrangian analysis developed under this effort is that it combines features of the so-called particle distribution approach with those of the so-called particle trajectory approach and can, in fact, be considered as a generalization of both of those traditional methods. The result of this generalization is a reduction in CPU time and memory requirements. Particle time step (stability) limitations have been eliminated by semi-implicit integration of the particle equations of motion (and, for certain applications, the particle temperature equation), although practical limitations remain in effect for reasons of accuracy. The analysis has been applied to the simulation of cavitating flow through a single-bladed section of a labyrinth seal. Models for the simulation of bubble formation and growth have been included, as well as models for bubble drag and heat transfer. The results indicate that bubble formation is more or less 'explosive'. for a given flow field, the number density of bubble nucleation sites is very sensitive to the vapor properties and the surface tension. The bubble motion, on the other hand, is much less sensitive to the properties, but is affected strongly by the local pressure gradients in the flow field. In situations where either the material properties or the flow field are not known with sufficient accuracy, parametric studies can be carried out rapidly to assess the effect of the important variables. Future work will include application of the analysis to cavitation in inducer flow fields.

Dejong, Frederik J.↗

Nighttime ion composition measurements at the geomagnetic equator

Two ion composition profiles, representative of the nighttime equatorial ionosphere between 90 km and 300 km, are presented. These profiles were obtained by two rocket-borne ion mass spectrometers on a single night for solar zenith angles of 112 deg and 165 deg. For both flights, the principal ion above 200 km is O(+). The downward drift of the atomic ions O(+) and N(+), coinciding with the postsunset lowering of the F2 peak, is observed through an enhancement of the density of O(+) at altitudes above 200 km and N(+) above 240 km. Below the drift region, O(+) and N(+) are observed in concentrations larger than expected. The NO(+) altitude distribution retains its shape throughout the night, and below 210 km, is the principal ion. The behavior of O2(+) can be explained by the O(+), electron density and theoretical neutral nitric oxide concentrations. Light metallic ions, including Mg(+), Na(+), and possibly Si(+), are observed to altitudes approaching 300 km and are affected by vertical drift.

Goldberg, R. A.↗

Circumferential flow analysis at the aft field joint of the Space Shuttle solid rocket motor

Flow analyses have been performed to determine the nature of the three-dimensional flow field in the vicinity of the aft-most field joint of the Space Shuttle Redesigned Solid Rocket Motor (RSRM). Specific objectives included the quantification of the circumferential pressure and velocity gradients at the joint location which might be caused by the non-uniform erosion of the rubber inhibitor which protrudes from the wall into the flow field. Three-dimensional Navier-Stokes equations have been solved in conjunction with the conservation equation for the turbulence energy and the dissipation rate. The numerical predictions have been compared with the measurements from a 7.5 percent scale cold flow model of the redesigned solid rocket motor.

Majumdar, Alok K.↗

Hybrid 3-D rocket trajectory program. Part 1: Formulation and analysis. Part 2: Computer programming and user's instruction

Models utilizing various sub-sets of the six degrees of freedom are used in trajectory simulation. A 3-D model with only linear degrees of freedom is especially attractive, since the coefficients for the angular degrees of freedom are the most difficult to determine and the angular equations are the most time consuming for the computer to evaluate. A computer program is developed that uses three separate subsections to predict trajectories. A launch rail subsection is used until the rocket has left its launcher. The program then switches to a special 3-D section which computes motions in two linear and one angular degrees of freedom. When the rocket trims out, the program switches to the standard, three linear degrees of freedom model.

Huang, L. C. P.↗

Numerical study of high-area-ratio H2/O2 rocket nozzles

Calculations were made for high-area-ratio H2/O2 rocket nozzles by using the RPLUS code which employs an implicit finite volume, lower-upper symmetric successive over-relaxation scheme to solve the Navier-Stokes equations and the species equations in a coupled manner. The combustion processes of hydrogen and oxygen are modeled by an 8-species and 18-step reaction mechanism, and the turbulence is simulated by the Baldwin-Lomax turbulence model with a pressure gradient correction on the van Driest's damping constant in the inner-layer model. The calculated results for two different rocket nozzles are presented as contours, profiles, and performance values. The predicted specific impulse and thrust agree well with the experimental data and those from the industry standard code. It is demonstrated that the RPLUS code can be used for analysis of rocket nozzles.

Kim, S. C.↗

Prediction of propellant behavior in spinning flow of a space motor

A model for 2-phase flow dynamics in a spinning spherical rocket motor, developed for the Ofeq satellite program, is based on the sectional approach to solution of the flow equations. An experimental model was built, to enable the motion of the aluminum/aluminum-oxide particles resulting from combustion within a solid-fuel motor to be simulated by injected paraffin droplets. The injected droplets included under-5 micrometer droplets which move with the gas stream and larger droplets, averaging 20 micrometer diameter, which simulate the motion of aluminum particles. The test chamber comprised a pair of cylindrical pyrex tubes with a sharp contraction in diameter, rotated at various speeds by a frequency-controllable motor. An optical system, based on a 5 W argon ion laser with a beam splitter and frequency shifter, mounted on a movable table, facilitated sectional measurements of the three velocity components and determination of size-velocity relationships. Preliminary results indicate that the effect of rotation on axial velocity is negligible, while its effect on tangential velocity approximates to solid-body rotation.

Gany, A.↗

Modeling and Validation of a Propellant Mixer for Controller Design

A mixing chamber used in rocket engine testing at the NASA Stennis Space Center is modelled by a system of two nonlinear ordinary differential equations. The mixer is used to condition the thermodynamic properties of cryogenic liquid propellant by controlled injection of the same substance in the gaseous phase. The three inputs of the mixer are the positions of the valves regulating the liquid and gas flows at the inlets, and the position of the exit valve regulating the flow of conditioned propellant. Mixer operation during a test requires the regulation of its internal pressure, exit mass flow, and exit temperature. A mathematical model is developed to facilitate subsequent controller designs. The model must be simple enough to lend itself to subsequent feedback controller design, yet its accuracy must be tested against real data. For this reason, the model includes function calls to thermodynamic property data. Some structural properties of the resulting model that pertain to controller design, such as uniqueness of the equilibrium point, feedback linearizability and local stability are shown to hold under conditions having direct physical interpretation. The existence of fixed valve positions that attain a desired operating condition is also shown. Validation of the model against real data is likewise provided.

Richter, Hanz↗

Probable Rotation States of Rocket Bodies in Low Earth Orbit

In order for Active Debris Removal to be accomplished, it is critically important to understand the probable rotation states of orbiting, spent rocket bodies. As compared to the question of characterizing small unresolved debris, in this problem there are several advantages: (1) objects are of known size, mass, shape and color, (2) they have typically been in orbit for a known period of time, (3) they are large enough that resolved images may be obtainable for verification of predicted orientation, and (4) the dynamical problem is simplified to first order by largely cylindrical symmetry. It is also nearly certain for realistic rocket bodies that internal friction is appreciable in the case where residual liquid or, to a lesser degree, unconsolidated solid fuels exist. Equations of motion have been developed for this problem in which internal friction as well as torques due to solar radiation, magnetic induction, and gravitational gradient are included. In the case of pure cylindrical symmetry, the results are compared to analytical predictions patterned after the standard approach for analysis of symmetrical tops. This is possible because solar radiation and gravitational torques may be treated as conservative. Agreement between results of both methods ensures their mutual validity. For monotone symmetric cylinders, solar radiation torque vanishes if the center of mass resides at the geometric center of the object. Results indicate that in the absence of solar radiation effects, rotation states tend toward an equilibrium configuration in which rotation is about the axis of maximum inertia, with the axis of minimum inertia directed toward the center of the earth. Solar radiation torque introduces a modification to this orientation. The equilibrium state is asymptotically approached within a characteristic timescale given by a simple ratio of relevant characterizing parameters for the body in question. Light curves are simulated for the expected asymptotic final rotation states of model objects, and these are compared to data derived from physical models of the same objects, tested in the Optical Measurements Center at JSC. Comparison to relevant light curves from actual orbiting rocket bodies are also performed, and diagnostic features of such curves are examined.

Ojakangas, Gregory W.↗

A numerical model for coupling between atomization and spray dynamics in liquid rocket thrust chambers

This paper describes a novel method of coupling the atomization and spray combustion processes encountered in coaxial injection elements of liquid rocket engine thrust chambers. This method is based on the Jet-Embedding technique in which the liquid jet core equations and the gas phase equations are solved separately. The liquid and gas phase solutions, however, are coupled through the boundary conditions at the interface between the phases. The computational grid for the gas phase calculations are adapted to the shape of the liquid jet core. The axial variation of droplet sizes are calculated using a stability analysis appropriate for the atomization regime of liquid jet breakup. The predictions of this method have been validated with experimental data on low speed water jets. Using this method, calculations are performed for the SSME fuel preburner single injector flow field. The results obtained are in good agreement with the predictions of the volume-of-fluid method.

Giridharan, M. G.↗