Search NASA⌕ Search

SEARCH · Search NASA

Results for “quasi-one-dimensional”

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 199 records · Page 11

THRSTER: A THRee-STream Ejector Ramjet Analysis and Design Tool

An engineering tool for analyzing ejectors in rocket based combined cycle (RBCC) engines has been developed. A key technology for multi-cycle RBCC propulsion systems is the ejector which functions as the compression stage of the ejector ramjet cycle. The THRee STream Ejector Ramjet analysis tool was developed to analyze the complex aerothermodynamic and combustion processes that occur in this device. The formulated model consists of three quasi-one-dimensional streams, one each for the ejector primary flow, the secondary flow, and the mixed region. The model space marches through the mixer, combustor, and nozzle to evaluate the solution along the engine. In its present form, the model is intended for an analysis mode in which the diffusion rates of the primary and secondary into the mixed stream are stipulated. The model offers the ability to analyze the highly two-dimensional ejector flowfield while still benefits from the simplicity and speed of an engineering tool. To validate the developed code, wall static pressure measurements from the Penn-State and NASA-ART RBCC experiments were used to compare with the results generated by the code. The calculated solutions were generally found to have satisfactory agreement with the pressure measurements along the engines, although further modeling effort may be required when a strong shock train is formed at the rocket exhaust. The range of parameters in which the code would generate valid results are presented and discussed.

Chue, R. S.↗

Solid Fuel Burning in Steady, Strained, Premixed Flow Fields: The Graphite/Air/Methane System

A detailed numerical investigation was conducted on the simultaneous burning of laminar premixed CH4/air flames and solid graphite in a stagnation flow configuration. The graphite and methane were chosen for this model, given that they are practical fuels and their chemical kinetics are considered as the most reliable ones among solid and hydrocarbon fuels, respectively. The simulation was performed by solving the quasi-one-dimensional equations of mass, momentum, energy, and species. The GRI 2.1 scheme was used for the gas-phase kinetics, while the heterogeneous kinetics were described by a six-step mechanism including stable and radical species. The effects of the graphite surface temperature, the gas-phase equivalence ratio, and the aerodynamic strain rate on the graphite burning rate and NO, production and destruction mechanisms were assessed. Results indicate that as the graphite temperature increases, its burning rate as well as the NO, concentration increase. Furthermore, it was found that by increasing the strain rate, the graphite burning rate increases as a result of the augmented supply of the gas-phase reactants towards the surface, while the NO, concentration decreases as a result of the reduced residence time. The effect of the equivalence ratio on both the graphite burning rate and NO, concentration was found to be non-monotonic and strongly dependent on the graphite temperature. Comparisons between results obtained for a graphite and a chemically inert surface revealed that the chemical activity of the graphite surface can result to the reduction of NO through reactions of the CH3, CH2, CH, and N radicals with NO.

Egolfopoulos, Fokion N.↗

The NASA Ames 16-Inch Shock Tunnel Nozzle Simulations and Experimental Comparison

The 16-Inch Shock Tunnel at NASA Ames Research Center is a unique test facility used for hypersonic propulsion testing. To provide information necessary to understand the hypersonic testing of the combustor model, computational simulations of the facility nozzle were performed and results are compared with available experimental data, namely static pressure along the nozzle walls and pitot pressure at the exit of the nozzle section. Both quasi-one-dimensional and axisymmetric approaches were used to study the numerous modeling issues involved. The facility nozzle flow was examined for three hypersonic test conditions, and the computational results are presented in detail. The effects of variations in reservoir conditions, boundary layer growth, and parameters of numerical modeling are explored.

TokarcikPolsky, S.↗

Progress in the Development of a Nozzle Design Methodology for Pulsed Detonation Engines

The Johns Hopkins University Applied Physics Laboratory (JHU/APL), in support of the NASA Glenn Research Center (NASA GRC), is investigating performance methodologies and system integration issues related to Pulsed Detonation Engine (PDE) nozzles. The primary goal of this ongoing effort is to develop design and performance assessment methodologies applicable to PDE exit nozzle(s). APL is currently focusing its efforts on a common plenum chamber design that collects the exhaust products from multiple PDE tubes prior to expansion in a single converging-diverging exit nozzle. To accomplish this goal, a time-dependent, quasi-one-dimensional analysis for determining the flow properties in and through a single plenum and exhaust nozzle is underway. In support of these design activities, parallel modeling efforts using commercial Computational Fluid Dynamics (CFD) software are on-going. These efforts include both two and three-dimensional as well as steady and time-dependent computations to assess the flow in and through these devices. This paper discusses the progress in developing this nozzle design methodology.

Leary, B. A.↗

Uncertainty Propagation for Turbulent, Compressible Flow in a Quasi-1D Nozzle Using Stochastic Methods

This paper describes a fully spectral, Polynomial Chaos method for the propagation of uncertainty in numerical simulations of compressible, turbulent flow, as well as a novel stochastic collocation algorithm for the same application. The stochastic collocation method is key to the efficient use of stochastic methods on problems with complex nonlinearities, such as those associated with the turbulence model equations in compressible flow and for CFD schemes requiring solution of a Riemann problem. Both methods are applied to compressible flow in a quasi-one-dimensional nozzle. The stochastic collocation method is roughly an order of magnitude faster than the fully Galerkin Polynomial Chaos method on the inviscid problem.

Zang, Thomas A.↗

Optimal Area Profiles for Ideal Single Nozzle Air-Breathing Pulse Detonation Engines

The effects of cross-sectional area variation on idealized Pulse Detonation Engine performance are examined numerically. A quasi-one-dimensional, reacting, numerical code is used as the kernel of an algorithm that iteratively determines the correct sequencing of inlet air, inlet fuel, detonation initiation, and cycle time to achieve a limit cycle with specified fuel fraction, and volumetric purge fraction. The algorithm is exercised on a tube with a cross sectional area profile containing two degrees of freedom: overall exit-to-inlet area ratio, and the distance along the tube at which continuous transition from inlet to exit area begins. These two parameters are varied over three flight conditions (defined by inlet total temperature, inlet total pressure and ambient static pressure) and the performance is compared to a straight tube. It is shown that compared to straight tubes, increases of 20 to 35 percent in specific impulse and specific thrust are obtained with tubes of relatively modest area change. The iterative algorithm is described, and its limitations are noted and discussed. Optimized results are presented showing performance measurements, wave diagrams, and area profiles. Suggestions for future investigation are also discussed.

Paxson, Daniel E.↗

Combustion Instabilities Modeled

NASA Lewis Research Center's Advanced Controls and Dynamics Technology Branch is investigating active control strategies to mitigate or eliminate the combustion instabilities prevalent in lean-burning, low-emission combustors. These instabilities result from coupling between the heat-release mechanisms of the burning process and the acoustic flow field of the combustor. Control design and implementation require a simulation capability that is both fast and accurate. It must capture the essential physics of the system, yet be as simple as possible. A quasi-one-dimensional, computational fluid dynamics (CFD) based simulation has been developed which may meet these requirements. The Euler equations of mass, momentum, and energy have been used, along with a single reactive species transport equation to simulate coupled thermoacoustic oscillations. A very simple numerical integration scheme was chosen to reduce computing time. Robust boundary condition procedures were incorporated to simulate various flow conditions (e.g., valves, open ends, and choked inflow) as well as to accommodate flow reversals that may arise during large flow-field oscillations. The accompanying figure shows a sample simulation result. A combustor with an open inlet, a choked outlet, and a large constriction approximately two thirds of the way down the length is shown. The middle plot shows normalized, time-averaged distributions of the relevant flow quantities, and the bottom plot illustrates the acoustic mode shape of the resulting thermoacoustic oscillation. For this simulation, the limit cycle peak-to-peak pressure fluctuations were 13 percent of the mean. The simulation used 100 numerical cells. The total normalized simulation time was 50 units (approximately 15 oscillations), which took 26 sec on a Sun Ultra2.

Paxson, Daniel E.↗

Pulse Detonation Engine Modeled

Pulse Detonation Engine Technology is currently being investigated at Glenn for both airbreathing and rocket propulsion applications. The potential for both mechanical simplicity and high efficiency due to the inherent near-constant-volume combustion process, may make Pulse Detonation Engines (PDE's) well suited for a number of mission profiles. Assessment of PDE cycles requires a simulation capability that is both fast and accurate. It should capture the essential physics of the system, yet run at speeds that allow parametric analysis. A quasi-one-dimensional, computational-fluid-dynamics-based simulation has been developed that may meet these requirements. The Euler equations of mass, momentum, and energy have been used along with a single reactive species transport equation, and submodels to account for dominant loss mechanisms (e.g., viscous losses, heat transfer, and valving) to successfully simulate PDE cycles. A high-resolution numerical integration scheme was chosen to capture the discontinuities associated with detonation, and robust boundary condition procedures were incorporated to accommodate flow reversals that may arise during a given cycle. The accompanying graphs compare experimentally measured and computed performance over a range of operating conditions for a particular PDE. Experimental data were supplied by Fred Schauer and Jeff Stutrud from the Air Force Research Laboratory at Wright-Patterson AFB and by Royce Bradley from Innovative Scientific Solutions, Inc. The left graph shows thrust and specific impulse, Isp, as functions of equivalence ratio for a PDE cycle in which the tube is completely filled with a detonable hydrogen/air mixture. The right graph shows thrust and specific impulse as functions of the fraction of the tube that is filled with a stoichiometric mixture of hydrogen and air. For both figures, the operating frequency was 16 Hz. The agreement between measured and computed values is quite good, both in terms of trend and magnitude. The error is under 10 percent everywhere except for the thrust value at an equivalence ratio of 0.8 in the left figure, where it is 14 percent. The simulation results shown were made using 200 numerical cells. Each cycle of the engine, approximately 0.06 sec, required 2.0 min of CPU time on a Sun Ultra2. The simulation is currently being used to analyze existing experiments, design new experiments, and predict performance in propulsion concepts where the PDE is a component (e.g., hybrid engines and combined cycles).

Paxson, Daniel E.↗

Baseline Computational Fluid Dynamics Methodology for Longitudinal-Mode Liquid-Propellant Rocket Combustion Instability

A computational method for the analysis of longitudinal-mode liquid rocket combustion instability has been developed based on the unsteady, quasi-one-dimensional Euler equations where the combustion process source terms were introduced through the incorporation of a two-zone, linearized representation: (1) A two-parameter collapsed combustion zone at the injector face, and (2) a two-parameter distributed combustion zone based on a Lagrangian treatment of the propellant spray. The unsteady Euler equations in inhomogeneous form retain full hyperbolicity and are integrated implicitly in time using second-order, high-resolution, characteristic-based, flux-differencing spatial discretization with Roe-averaging of the Jacobian matrix. This method was initially validated against an analytical solution for nonreacting, isentropic duct acoustics with specified admittances at the inflow and outflow boundaries. For small amplitude perturbations, numerical predictions for the amplification coefficient and oscillation period were found to compare favorably with predictions from linearized small-disturbance theory as long as the grid exceeded a critical density (100 nodes/wavelength). The numerical methodology was then exercised on a generic combustor configuration using both collapsed and distributed combustion zone models with a short nozzle admittance approximation for the outflow boundary. In these cases, the response parameters were varied to determine stability limits defining resonant coupling onset.

Litchford, R. J.↗

Updates to Multi-Dimensional Flux Reconstruction for Hypersonic Simulations on Tetrahedral Grids

The quality of simulated hypersonic stagnation region heating with tetrahedral meshes is investigated by using an updated three-dimensional, upwind reconstruction algorithm for the inviscid flux vector. An earlier implementation of this algorithm provided improved symmetry characteristics on tetrahedral grids compared to conventional reconstruction methods. The original formulation however displayed quantitative differences in heating and shear that were as large as 25% compared to a benchmark, structured-grid solution. The primary cause of this discrepancy is found to be an inherent inconsistency in the formulation of the flux limiter. The inconsistency is removed by employing a Green-Gauss formulation of primitive gradients at nodes to replace the previous Gram-Schmidt algorithm. Current results are now in good agreement with benchmark solutions for two challenge problems: (1) hypersonic flow over a three-dimensional cylindrical section with special attention to the uniformity of the solution in the spanwise direction and (2) hypersonic flow over a three-dimensional sphere. The tetrahedral cells used in the simulation are derived from a structured grid where cell faces are bisected across the diagonal resulting in a consistent pattern of diagonals running in a biased direction across the otherwise symmetric domain. This grid is known to accentuate problems in both shock capturing and stagnation region heating encountered with conventional, quasi-one-dimensional inviscid flux reconstruction algorithms. Therefore the test problems provide a sensitive indicator for algorithmic effects on heating. Additional simulations on a sharp, double cone and the shuttle orbiter are then presented to demonstrate the capabilities of the new algorithm on more geometrically complex flows with tetrahedral grids. These results provide the first indication that pure tetrahedral elements utilizing the updated, three-dimensional, upwind reconstruction algorithm may be used for the simulation of heating and shear in hypersonic flows in upwind, finite volume formulations.

Gnoffo, Peter A.↗

Control Volume Analysis of Boundary Layer Ingesting Propulsion Systems With or Without Shock Wave Ahead of the Inlet

The performance benefit of boundary layer or wake ingestion on marine and air vehicles has been well documented and explored. In this article, a quasi-one-dimensional boundary layer ingestion (BLI) benefit analysis for subsonic and transonic propulsion systems is performed using a control volume of a ducted propulsion system that ingests the boundary layer developed by the external airframe surface. To illustrate the BLI benefit, a relationship between the amount of BLI and the net thrust is established and analyzed for two propulsor types. One propulsor is an electric fan, and the other is a pure turbojet. These engines can be modeled as a turbofan with an infinite bypass ratio for the electric fan, and with a zero bypass ratio for the pure turbojet. The analysis considers two flow processes: a boundary layer being ingested by an aircraft inlet and a shock wave sitting in front of the inlet. Though the two processes are completely unrelated, both represent a loss of total pressure and velocity. In real applications, it is possible to have both processes occurring in front of the inlet of a transonic vehicle. Preliminary analysis indicates that the electrically driven propulsion system benefits most from the boundary layer ingestion and the presence of transonic shock waves, whereas the benefit for the turbojet engine is near zero or negative depending on the amount of total temperature rise across the engine.

Kim, Hyun Dae↗

Analysis of Water Recovery Rate from the Heat Melt Compactor

Human space missions generate trash with a substantial amount of plastic (20% or greater by mass). The trash also contains water trapped in food residue and paper products and other trash items. The Heat Melt Compactor (HMC) under development by NASA Ames Research Center (ARC) compresses the waste, dries it to recover water and melts the plastic to encapsulate the compressed trash. The resulting waste disk or puck represents an approximately ten-fold reduction in the volume of the initial trash loaded into the HMC. In the current design concept being pursued, the trash is compressed by a piston after it is loaded into the trash chamber. The piston face, the side walls of the waste processing chamber and the end surface in contact with the waste can be heated to evaporate the water and to melt the plastic. Water is recovered by the HMC in two phases. The first is a pre-process compaction without heat or with the heaters initially turned on but before the waste heats up. Tests have shown that during this step some liquid water may be expelled from the chamber. This water is believed to be free water (i.e., not bound with or absorbed in other waste constituents) that is present in the trash. This phase is herein termed Phase A of the water recovery process. During HMC operations, it is desired that liquid water recovery in Phase A be eliminated or minimized so that water-vapor processing equipment (e.g., condensers) downstream of the HMC are not fouled by liquid water and its constituents (i.e., suspended or dissolved matter) exiting the HMC. The primary water recovery process takes place next where the trash is further compacted while the heated surfaces reach their set temperatures for this step. This step will be referred to herein as Phase B of the water recovery process. During this step the waste chamber may be exposed to different selected pressures such as ambient, low pressure (e.g., 0.2 atm), or vacuum. The objective for this step is to remove both bound and any remaining free water in the trash by evaporation. The temperature settings of the heated surfaces are usually kept above the saturation temperature of water but below the melting temperature of the plastic in the waste during this step to avoid any encapsulation of wet trash which would reduce the amount of recovered water by blocking the vapor escape. In this paper, we analyze the water recovery rate during Phase B where the trash is heated and water leaves the waste chamber as vapor, for operation of the HMC in reduced gravity. We pursue a quasi-one-dimensional model with and without sidewall heating to determine the water recovery rate and the trash drying time. The influences of the trash thermal properties, the amount of water loading, and the distribution of the water in the trash on the water recovery rates are determined.

Heat Melt Compactor↗

Hall-Effect Thruster Simulations with 2-D Electron Transport and Hydrodynamic Ions

A computational approach that has been used extensively in the last two decades for Hall thruster simulations is to solve a diffusion equation and energy conservation law for the electrons in a direction that is perpendicular to the magnetic field, and use discrete-particle methods for the heavy species. This "hybrid" approach has allowed for the capture of bulk plasma phenomena inside these thrusters within reasonable computational times. Regions of the thruster with complex magnetic field arrangements (such as those near eroded walls and magnets) and/or reduced Hall parameter (such as those near the anode and the cathode plume) challenge the validity of the quasi-one-dimensional assumption for the electrons. This paper reports on the development of a computer code that solves numerically the 2-D axisymmetric vector form of Ohm's law, with no assumptions regarding the rate of electron transport in the parallel and perpendicular directions. The numerical challenges related to the large disparity of the transport coefficients in the two directions are met by solving the equations in a computational mesh that is aligned with the magnetic field. The fully-2D approach allows for a large physical domain that extends more than five times the thruster channel length in the axial direction, and encompasses the cathode boundary. Ions are treated as an isothermal, cold (relative to the electrons) fluid, accounting for charge-exchange and multiple-ionization collisions in the momentum equations. A first series of simulations of two Hall thrusters, namely the BPT-4000 and a 6-kW laboratory thruster, quantifies the significance of ion diffusion in the anode region and the importance of the extended physical domain on studies related to the impact of the transport coefficients on the electron flow field.

plasma↗

Passive Rocket Diffuser Theory: A Re-Examination of Minimum Second Throat Size

Second-throat diffusers serve to isolate rocket engines from the effects of ambient back pressure during testing without using active control systems. Among the most critical design parameters is the relative area of the diffuser throat to that of the nozzle throat. A smaller second throat is generally desirable because it decreases the stagnation-to-ambient pressure ratio the diffuser requires for nominal operation. There is a limit, however. Below a certain size, the second throat can cause pressure buildup within the diffuser and prevent it from reaching the start condition that protects the nozzle from side-load damage. This paper presents a method for improved estimation of the minimum second throat area which enables diffuser start. The new 3-zone model uses traditional quasi-one-dimensional compressible flow theory to approximate the structure of two distinct diffuser flow fields observed in Computational Fluid Dynamics (CFD) simulations and combines them to provide a less-conservative estimate of the second throat size limit. It is unique among second throat sizing methods in that it accounts for all major conical nozzle and second throat diffuser design parameters within its limits of application. The performance of the 3-zone method is compared to the historical normal shock and force balance methods, and verified against a large number of CFD simulations at specific heat ratios of 1.4 and 1.25. Validation is left as future work, and the model is currently intended to function only as a first-order design tool.

Jones, Daniel R.↗

Hybrid Reduced Order Model for N2-N2 Interactions for Application to Dissociation and Energy Transfer Processes

Recent work in the aerothermodynamics community has focused on the development of reduced-order models for thermo-chemical non-equilibrium which avoid the restrictive assumptions of multi-temperature models and the prohibitive cost associated with State-to-State (StS) models. In the present work, this is accomplished by lumping energy states together and assuming groups of states are roughly in equilibrium. As a result, the non-equilibrium behavior of a gas can be captured at a reduced computational cost from a full StS simulation. In this work, we present a hybrid grouping model for studying energy transfer and dissociation in a mixture of nitrogen molecules due to N2-N2 reactions. This is accomplished by making use of a grouping strategy informed by data from the N2-N StS kinetic data. However, due to the massive computational cost associated with constructing StS data for the N2-N2 system, the kinetic data for the hybrid grouping model are calculated using the quasi-classical trajectory (QCT) method by sampling states for trajectory within the groups. This general framework is called the Maximum-Entropy Quasi-Classical Trajectory (ME-QCT) method. The primary challenge associated with this method is that rates for reverse grouped reactions cannot be obtained through detailed balance at a group level, due to the variation of group internal temperatures. To construct the full model for N2-N2 grouped kinetics using the ME-QCT method, detailed balance is invoked at the microscopic level, allowing for the calculation of the full kinetic data from QCT. Results will be presented using the full ME-QCT model for the N2-N2 system in an isothermal and isochoric reactor simulation. In addition, simple CFD test cases for a one-dimensional standing shock and for a quasi-one-dimensional nozzle will be used for demonstration of the ME-QCT method. This method allows for the calculation of non-equilibrium behavior for the N2-N2 system without the prohibitive cost of a full StS simulation. Moreover, it enables the construction of a unified model for the dissociating and recombining non-equilibrium flows.

nonequilibrium modeling↗

Physical Simulation of Rocket Exhaust Aerodynamics Using Heated Ethane: Conceptual Foundations

Heated ethane (C2H6) has been proposed as an alternative to inert gases for use as a motive fluid in the experimental simulation of rocket exhaust plumes. By adjusting stagnation temperature, the isentropic exponent of ethane can be tuned to approximate those produced by common rocket propellants including hydrogen, hypergols, alcohols, and hydrocarbons. As a result, ethane can be made to follow a nozzle expansion process which is nearly identical to realistic rocket engine flow fields. Additionally, its high auto-ignition temperature and resistance to condensation enable the testing of expansion ratios much larger than conventional inertgas testing. NASA SSC has performed quasi-one-dimensional analyses using the Chemical Equilibrium with Applications (CEA) code as a preliminary means to compare flow fields produced by non-reacting ethane to those of reacting combustion products. A LO2/LH2 rocket engine operating at a chamber pressure of 5.0 MPa and a mixture ratio of 6.1 was used as an example case to demonstrate ethane’s efficacy as a simulant. Errors for key similarity parameters were compared to legacy cold-flow test methods. Additional errors induced by machining tolerances and chemical impurities were also examined. Results suggest that at a 3% geometric scale and ~500 K ethane stagnation temperature, an error of less than 2.5% throughout the flow field is realistically achievable along the dimensions of Mach number, Reynolds number, pressure ratio, and isentropic exponent. The development of an experimental test bed for validation of this configuration is currently underway.

Jones, Daniel R.↗

LGGUN User's Manual

LGGUN is a quasi-one-dimensional Godunov code1 which is second order accurate in time and third order accurate in space, uses realistic equations of state for all media, includes wall friction and heat transfer for all gas and dense media zones and includes a simple non-equilibrium model for gas phase turbulence. This manual details the theory and operation of the code.

Ballistic Range↗

Correlated Spatio-temporal Evolution of Extreme-Ultraviolet Ribbons and Hard X-Rays in a Solar Flare

We analyze the structure and evolution of ribbons from the M7.3 SOL2014-04-18T13 flare using ultraviolet images from the Interface Region Imaging Spectrograph and the Solar Dynamics Observatory (SDO)/Atmospheric Imaging Assembly (AIA), magnetic data from the SDO/Helioseismic and Magnetic Imager, hard X-ray (HXR) images from the Reuven Ramaty High Energy Solar Spectroscopic Imager, and light curves from the Fermi/Gamma-ray Burst Monitor, in order to infer properties of coronal magnetic reconnection. As the event progresses, two flare ribbons spread away from the magnetic polarity inversion line. The width of the newly brightened front along the extension of the ribbon is highly intermittent in both space and time, presumably reflecting nonuniformities in the structure and/or dynamics of the flare current sheet. Furthermore, the ribbon width grows most rapidly in regions exhibiting concentrated nonthermal HXR emission, with sharp increases slightly preceding the HXR bursts. The light curve of the ultraviolet emission matches the HXR light curve at photon energies above 25 keV. In other regions the ribbon-width evolution and light curves do not temporally correlate with the HXR emission. This indicates that the production of nonthermal electrons is highly nonuniform within the flare current sheet. Our results suggest a strong connection between the production of nonthermal electrons and the locally enhanced perpendicular extent of flare ribbon fronts, which in turn reflects the inhomogeneous structure and/or reconnection dynamics of the current sheet. Despite this variability, the ribbon fronts remain nearly continuous, quasi-one-dimensional features. Thus, although the reconnecting coronal current sheets are highly structured, they remain quasi-two-dimensional and the magnetic energy release occurs systematically, rather than stochastically, through the volume of the reconnecting magnetic flux.

S. J. Naus↗