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.

154 records · Page 9

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↗

SUPIN: A Tool for the Aerodynamic Design and Analysis of Supersonic Inlets

The Supersonic Inlet Design and Analysis Tool (SUPIN) is a computational program that performs geometric modeling and aerodynamic design and analysis of a set of supersonic inlets whose types include the axisymmetric pitot, three-dimensional pitot, axisymmetric, two-dimensional, two-dimensional bifurcated-duct, and streamline-traced inlets. The inlets are modeled by dividing the inlet flowpath into components that start at the freestream and ends at the engine face. The primary components of the inlet include the external supersonic diffuser, throat section, and subsonic diffuser. Each component is characterized by a set of inputs that include geometric and aerodynamic factors. Examples of the geometric factors are angles, lengths, area ratios, and coordinates that describe the geometry of the components. Examples of the aerodynamic factors are Mach numbers, flow ratios, and pressure ratios that set constraints or desired aerodynamic properties of the component. These input factors are specified within a text-based input data file that is read into SUPIN. SUPIN is a Fortran 95 program that performs the inlet design and analysis as a batch or non-interactive process. SUPIN writes output data files describing the inlet geometry and aerodynamic performance. SUPIN uses analytic, empirical, and numerical methods to perform the inlet design and analysis. The inlet geometry is constructed using planar and three-dimensional geometry constructs and is represented by coordinates, angles, areas, profiles, and surfaces. Surfaces are output in the form of a Plot3D surface grid file. SUPIN does have some capability to automatically generate planar and three-dimensional, multi-block, structured grids for computational fluid dynamics (CFD) analysis. The aerodynamic analysis performed by SUPIN computes quasi-one-dimensional properties, such as Mach numbers, pressures, and temperatures, through the inlet flowpath. Planar aerodynamic properties, such as those for shock and Mach waves, are computed for the supersonic compression of the external supersonic diffuser and throat section. The inlet performance is characterized by the inlet flow rates, total pressure recovery, and inlet drag. The primary intent of this document is to describe the usage of SUPIN. The methods used within SUPIN are documented and may serve as a useful reference for understanding the aerodynamic design and analysis of supersonic inlets.

supersonic inlets↗