Search NASA⌕ Search

SEARCH · Search NASA

Results for “kinetic modeling”

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

Analysis of Nitrogen Shockwave Radiation at Velocity Around 7km/s

In the last decade, the Electric Arc Shock Tube (EAST) facility has been widely used to provide experimental data aiming to validate radiative heating models used for atmospheric entry, mostly at velocities above 9km/s in air. In more recent years, another experimental campaign aimed to look at velocity below 9 km/s and its analysis resulted in a significant improvement of the air kinetic model using new rates found in the literature but discrepancies were still observed, especially for the N2 second positive system. Despite these previous studies, the overall uncertainties related to surface heating for lunar return remain around 50% To reduce these uncertainties, the last experimental campaign at EAST, named test 62, was conducted in pure nitrogen at velocity ranging from 6 to 11 km/s at a freestream pressure of 0.2 Torr to provide more fundamental data on nitrogen. A previous study analyzed the data obtained at a velocity of 10.3 km/s, where atomic radiation is dominant, and provided spatially resolved density measurements of different states of nitrogen. This study will focus on experimental data obtained at shock velocity of 6.88 and 7.15 km/s where molecular nitrogen is expected to be the dominant source of radiation.

Augustin Tibere-Inglesse↗

Introducing the VRT gas turbine combustor

An innovative annular combustor configuration is being developed for aircraft and other gas turbine engines. This design has the potential of permitting higher turbine inlet temperatures by reducing the pattern factor and providing a major reduction in NO(x) emission. The design concept is based on a Variable Residence Time (VRT) technique which allows large fuel particles adequate time to completely burn in the circumferentially mixed primary zone. High durability of the combustor is achieved by dual function use of the incoming air. The feasibility of the concept was demonstrated by water analogue tests and 3-D computer modeling. The computer model predicted a 50 percent reduction in pattern factor when compared to a state of the art conventional combustor. The VRT combustor uses only half the number of fuel nozzles of the conventional configuration. The results of the chemical kinetics model require further investigation, as the NO(x) predictions did not correlate with the available experimental and analytical data base.

Melconian, Jerry O.↗

Introducing the VRT gas turbine combustor

An innovative annular combustor configuration is being developed for aircraft and other gas turbine engines. This design has the potential of permitting higher turbine inlet temperatures by reducing the pattern factor and providing a major reduction in NO(x) emission. The design concept is based on a Variable Residence Time (VRT) technique which allows large fuel particles adequate time to completely burn in the circumferentially mixed primary zone. High durability of the combustor is achieved by dual function use of the incoming air. The feasibility of the concept was demonstrated by water analogue tests and 3-D computer modeling. The computer model predicted a 50 percent reduction in pattern factor when compared to a state of the art conventional combustor. The VRT combustor uses only half the number of fuel nozzles of the conventional configuration. The results of the chemical kinetics model require further investigation, as the NO(x) predictions did not correlate with the available experimental and analytical data base.

Melconian, Jerry O.↗

Explicit least squares system parameter identification for exact differential input/output models

The equation error for a class of systems modeled by input/output differential operator equations has the potential to be integrated exactly, given the input/output data on a finite time interval, thereby opening up the possibility of using an explicit least squares estimation technique for system parameter identification. The paper delineates the class of models for which this is possible and shows how the explicit least squares cost function can be obtained in a way that obviates dealing with unknown initial and boundary conditions. The approach is illustrated by two examples: a second order chemical kinetics model and a third order system of Lorenz equations.

Pearson, A. E.↗

Induced nucleation of carbon dust in red giant stars

This study quantitatively tests the proposed model of induced nucleation of carbonaceous grains in carbon-rich red giant stars. Induced nucleation is the process of grain growth initiated by the presence of reactive surfaces provided by seed particles. The numerical study was performed using a deailed chemical kinetic model of carbon deposition, grain coagulation, and homogeneous nucleation of polycyclic aromatic hydrocarbons (PAHs). The model uses a method of moments to keep track of developing grain population in the forming dust shell. We test the efficiency of grain formation for large ranges of dust shell parameters typical for carbon stars. Our model is capable of producing a range of optically thick and thin dust shells in carbon stars. Results are in accord with (IRAS) spectral classes of carbon stars. The resulting composite grains produced are consistent with those recently found in ancient meteorites. This model also provides a realistic explanation for high abundances of (PAHs) in the interstellar medium and some planetary nebulae.

Cadwell, Brian J.↗

Physical and numerical sources of computational inefficiency in integration of chemical kinetic rate equations: Etiology, treatment and prognosis

The design of a very fast, automatic black-box code for homogeneous, gas-phase chemical kinetics problems requires an understanding of the physical and numerical sources of computational inefficiency. Some major sources reviewed in this report are stiffness of the governing ordinary differential equations (ODE's) and its detection, choice of appropriate method (i.e., integration algorithm plus step-size control strategy), nonphysical initial conditions, and too frequent evaluation of thermochemical and kinetic properties. Specific techniques are recommended (and some advised against) for improving or overcoming the identified problem areas. It is argued that, because reactive species increase exponentially with time during induction, and all species exhibit asymptotic, exponential decay with time during equilibration, exponential-fitted integration algorithms are inherently more accurate for kinetics modeling than classical, polynomial-interpolant methods for the same computational work. But current codes using the exponential-fitted method lack the sophisticated stepsize-control logic of existing black-box ODE solver codes, such as EPISODE and LSODE. The ultimate chemical kinetics code does not exist yet, but the general characteristics of such a code are becoming apparent.

Pratt, D. T.↗

Production, Outflow Velocity, and Radial Distribution of H2O and OH in the Coma of Comet C/1995 O1 [Hale-Bopp] from Wide Field Imaging of OH

Observations of OH are a useful proxy of the water production rate (Q(sub H2O)) and outflow velocity (V(sub out)) in comets. From wide field images taken on 03/28/1997 and 04/08/1997 that capture the entire scale length of the OH coma of comet C/1995 O1 (Hale-Bopp), we obtain Q(sub H2O) from the model-independent method of aperture summation. With an adaptive ring summation algorithm, we extract the radial brightness distribution of OH 0-0 band emission out to cometocentric distances of up to 10(exp 6) km, both as azimuthal averages and in quadrants covering different position angles relative to the comet-Sun line. These profiles are fit using both fixed and variable velocity 2-component spherical expansion models to estimate V(sub OH) with increasing distance from the nucleus. The OH coma of Hale-Bopp was more spatially extended than previous comets, and this extension is best matched by a variable acceleration of H2O and OH that acted across the entire coma, but was strongest within 1-2 x 10(exp 4) km from the nucleus. Our models indicate that V(sub OH) at the edge of our detectable field of view (10(exp 6) km) was approx. 2-3 times greater in Hale-Bopp than for a 1P/Halley-class comet at 1 AU, which is consistent with the results of more sophisticated gas-kinetic models, extrapolation from previous observations of OH in comets with Q(sub H2O) greater than 10(exp 29)/s , and direct radio measurements of the outer coma Hale-Bopp OH velocity. The most probable source of this acceleration is thermalization of the excess energy of dissociation of H2O and OH over an extended collisional coma. When the coma is broken down by quadrants in position angle, we find an azimuthal asymmetry in the radial distribution that is characterized by an increase in the spatial extent of OH in the region between the orbit-trailing and anti-sunward directions. Model fits specific to this area and comparison with radio OH measurements suggest greater acceleration here, with V(sub OH) approx. 1.5 times greater at a 10(exp 6) km cometocentric distance than elsewhere in the coma. We discuss several mechanisms that may have acted within the coma to produce the observed effect.

Harris, Walter M.↗

An ODE-Fitting Approach to Estimate Critical Aircraft Performance Parameters for Trajectory Prediction

Ground-based decision support tools (DST) in air traffic management (ATM) typically perform trajectory prediction based on aircraft performance model (APM) parameters, but some or all of these parameters might not be readily available. In particular, the three critical parameters required for trajectory prediction are the thrust setting, drag coefficients, and takeoff weight of the aircraft. Unfortunately, these parameters are coupled and appear together in physics-based kinetic models. Past approaches utilize data from a specific phase of flight (climb, level flight or descent), where one or more of the parameters are assumed to be known and estimate the remaining unknown parameters. This approach introduces bias/errors and also does not extend to scenarios where all of the above parameters are not known with sufficient accuracy. This paper is the first of its kind to propose a generalized framework for simultaneous estimation of all three critical APM parameters (thrust, drag, and mass). The proposed approach utilizes data from both the climb and descent phases and fits the ordinary differential equation for altitude in each phase using historical trajectory data available from radar tracks or ADS-B. The approach yields a set of optimized APM parameters that are best suited to fit each historical flight record. The methodology is applied on on sample flights from three different aircraft types, and the results demonstrate low fit error and consequently will yield a high level of prediction accuracy.

trajectory prediction, aircraft performance model,↗

An ODE-Fitting Approach to Estimate Critical Aircraft Performance Parameters for Trajectory Prediction

Ground-based decision support tools (DST) in air traffic management (ATM) typically perform trajectory prediction based on aircraft performance model (APM) parameters, but some or all of these parameters might not be readily available. In particular, the three critical parameters required for trajectory prediction are the thrust setting, drag coefficients, and takeoff weight of the aircraft. Unfortunately, these parameters are coupled and appear together in physics-based kinetic models. Past approaches utilize data from a specific phase of flight (climb, level flight or descent), where one or more of the parameters are assumed to be known and estimate the remaining unknown parameters. This approach introduces bias/errors and also does not extend to scenarios where all of the above parameters are not known with sufficient accuracy. This paper is the first of its kind to propose a generalized framework for simultaneous estimation of all three critical APM parameters (thrust, drag, and mass). The proposed approach utilizes data from both the climb and descent phases and fits the ordinary differential equation for altitude in each phase using historical trajectory data available from radar tracks or ADS-B. The approach yields a set of optimized APM parameters that are best suited to fit each historical flight record. The methodology is applied on on sample flights from three different aircraft types, and the results demonstrate low fit error and consequently will yield a high level of prediction accuracy.

trajectory prediction↗

Modeling of hypersonic reacting flows

The effect of nonequilibrium thermochemical processes on the aerodynamic and heat transfer characteristics of a hypersonic vehicle is illustrated with examples. It is shown that the conventional method of predicting chemical reactions always predict that the flow is closer to equilibrium than it actually is, leading to incorrect predictions of aerodynamic characteristics of a vehicle. A method is presented for predicting chemical processes using a two-temperature kinetic model, which is developed on the basis of combining conservation equations for vibrational energy and for electron electronic energy, which are also derived. The validity of the two-temperature model is proven using two types of experimental data: radiation and shock shapes.

Park, Chul↗

The Effects of Consistent Chemical Kinetics Calculations on the Pressure-Temperature Profiles and Emission Spectra of Hot Jupiters

In this work we investigate the impact of calculating non-equilibrium chemical abundances consistently with the temperature structure for the atmospheres of highly-irradiated, close-in gas giant exoplanets. Chemical kinetics models have been widely used in the literature to investigate the chemical compositions of hot Jupiter atmospheres which are expected to be driven away from chemical equilibrium via processes such as vertical mixing and photochemistry. All of these models have so far used pressure-temperature (P-T) profiles as fixed model input. This results in a decoupling of the chemistry from the radiative and thermal properties of the atmosphere, despite the fact that in nature they are intricately linked. We use a one-dimensional radiative-convective equilibrium model, ATMO, which includes a sophisticated chemistry scheme to calculate P-T profiles which are fully consistent with non-equilibrium chemical abundances, including vertical mixing and photochemistry. Our primary conclusion is that, in cases of strong chemical disequilibrium, consistent calculations can lead to differences in the P-T profile of up to 100 K compared to the P-T profile derived assuming chemical equilibrium. This temperature change can, in turn, have important consequences for the chemical abundances themselves as well as for the simulated emission spectra. In particular, we find that performing the chemical kinetics calculation consistently can reduce the overall impact of non-equilibrium chemistry on the observable emission spectrum of hot Jupiters. Simulated observations derived from non-consistent models could thus yield the wrong interpretation. We show that this behaviour is due to the non-consistent models violating the energy budget balance of the atmosphere.

planets and satellites↗

Magnetogasdynamic Power Extraction and Flow Conditioning for a Gas Turbine

An extension of the Russian AJAX concept to a turbojet is being explored. This magnetohydrodynamic (MHD) energy bypass engine cycle incorporating conventional gas turbine technology has MHD flow conditioning at the inlet to electromagnetically extract part of the inlet air kinetic energy. The electrical power generated can be used for various on-board vehicle requirements including plasma flow control around the vehicle or it may be used for augmenting the expanding flow in the high speed nozzle by MHD forces to generate more thrust. In order to achieve this interaction, the air needs to be ionized by an external means even up to fairly high flight speeds, and the leading candidates may be classified as electrical discharge devices. The present kinetic modeling calculations suggest that the use of electron beams with characteristics close to the commercially available e-beam systems (electron energy approx. 60 keV, beam current approx. 0.2 mA/sq cm) to sustain ionization in intermediate pressure, low-temperature (P = 0.1 atm, T = 300 K) supersonic air flows allows considerable reduction of the flow kinetic energy (up to 10 to 20 percent in M = 3 flows). The calculations also suggest that this can be achieved at a reasonable electron beam efficiency (eta approx. 5), even if the e-beam window losses are taken into account. At these conditions, the exit NO and O atom concentrations due to e-beam initiated chemical reactions do not exceed 30 ppm. Increasing the beam current up to approx. 2 mA/sq cm, which corresponds to a maximum electrical conductivity of sigma(sub max) approx. 0.8 mho/m at the loading parameter of K = 0.5, would result in a much greater reduction of the flow kinetic energy (up to 30 to 40 percent). The MHD channel efficiency at these conditions would be greatly reduced (to eta approx. 1) due to increased electron recombination losses in the channel. At these conditions, partial energy conversion from kinetic energy to heat would result in a significant total pressure loss (P(sub 0)/P(sub 0i) approx. 0.3). The total pressure loss can be reduced operating at the loading parameter closer to unity, at the expense of the reduced electrical power output. Raising the beam current would also result in the increase of the exit O atom concentrations (up to 600 ppm) and NO (up to 150 ppm).

Adamovich, Igor V.↗

Numerical Simulation of Vitiation Effects on a Hydrogen-Fueled Dual-Mode Scramjet

The Wind-US computational fluid dynamics (CFD) flow solver was used to simulate dual-mode direct-connect ramjet/scramjet engine flowpath tests conducted in the University of Virginia (UVa) Supersonic Combustion Facility (SCF). The objective was to develop a computational capability within Wind-US to aid current hypersonic research and provide insight to flow as well as chemistry details that are not resolved by instruments available. Computational results are compared with experimental data to validate the accuracy of the numerical modeling. These results include two fuel-off non-reacting and eight fuel-on reacting cases with different equivalence ratios, split between one set with a clean (non-vitiated) air supply and the other set with a vitiated air supply (12 percent H2O vapor). The Peters and Rogg hydrogen-air chemical kinetics model was selected for the scramjet simulations. A limited sensitivity study was done to investigate the choice of turbulence model and inviscid flux scheme and led to the selection of the k-epsilon model and Harten, Lax and van Leer (for contact waves) (HLLC) scheme for general use. Simulation results show reasonably good agreement with experimental data and the overall vitiation effects were captured.

Vyas, Manan A.↗

Integrated Physics-based Modeling and Experiments for Improved Prediction of Combustion Dynamics in Low-Emission Systems

Concurrent simulation and experiment was undertaken to assess the ability of a hybrid RANS-LES model to predict combustion dynamics in a single-element lean direct-inject (LDI) combustor showing self-excited instabilities. High frequency pressure modes produced by Fourier and modal decomposition analysis were compared quantitatively, and trends with equivalence ratio and inlet temperature were compared qualitatively. High frequency OH PLIF and PIV measurements were also taken. Submodels for chemical kinetics and primary and secondary atomization were also tested against the measured behavior. For a point-wise comparison, the amplitudes matched within a factor of two. The dependence on equivalence ratio was matched. Preliminary results from simulation using an 18-reaction kinetics model indicated instability amplitudes closer to measurement. Analysis of the simulations suggested a band of modes around 1400 Hz were due to a vortex bubble breakdown and a band of modes around 6 kHz were due to a precessing vortex core hydrodynamic instability. The primary needs are directly coupled and validated ab initio models of the atomizer free surface flow and the primary atomization processes, and more detailed study of the coupling between the 3D swirling flow and the local thermoacoustics in the diverging venturi section.

combustion↗

Seventh Annual V. M. Goldschmidt Conference

Topics considered include: Subduction of the Aseismic Cocos Ridge Displaced Magma Sources Beneath the Cordillera de Talamanca, Costa Rica; Topography of Transition Zone Discontinuities: A Measure of 'Olivine' Content and Evidence for Deep Cratonic Roots; Uranium Enrichment in Lithospheric Mantle: Case Studies from French Massif Central; Rare-Earth-Element Anomalies in the Decollement Zone of the nankai Accretionary Prism, Japan: Evidence of Fluid Flow?; Rare Earth Elements in Japanese Mudrocks: The Influence of Provenance; The Evolution of Seawater Strontium Isotopes in the Last Hundred Million Years: Reinterpretation and Consequences for Erosion and Climate Models; From Pat to Tats: The Lead Isotope Legacy in the Studies of the Continental Crust-Upper Mantle System; Geochronology of the Jack Hills Detrital Zircons by Precise Uranium-Lead Isotope-Dilution Analysis of Crystal Fragments; Iridium in the Oceans; The Helium-Heat-Lead Paradox; Control of Distribution Patterns of Heavy Metals in Ganga Plain Around Kanpur Region, India, by Fluvial Geomorphic Domains; Geochemical and Isotopic Features of Ferrar Magmatic Provience (Victoria Land, Antarctica); Rare Earth Elements in Marine Fine-Grained Sediments from the Northwestern Portuguese Shelf (Atlantic); Aspects of Arc Fluxes; General Kinetic Model for Dolomite Precipitation Rate with Application to the Secular History of Seawater Composition; High-Precision Uranium-series Chronology from Speleothems; Trace-Element Modeling of Aqueous Fluid-Peridotite Interaction in the Mantle Wedge of Subduction Zones; Rainfall Variations in Southeastern Australia over the Last 500,000 Years from Speleothem Deposition; The Role of Water in High-Pressure Fluids; The Kinetic Conditions of Metamorphic Minearogenesis: Evidence from Minerals and Assemblages.

Source record↗

The excitation and distribution of CO (J = 6-5) emission in the Orion Nebula

Observations of the 691-GHz (J = 6-5) transition of CO in the BN/KL region of Orion obtained in February 1981, at the IR Telescope Facility at Mauna Kea are reported. The system employs a heterodyne receiver with an overall noise temperature of 3900 K DSB at 432 microns, 64 5-MHz IF-filter-bank channels, and a chopping secondary with 120-arcsec excursion, and has 35-arcsec resolution. Sample data are presented graphically and analyzed using a rate equation and a kinetic model. A 35 x 45-arcsec core with hot broad plateau emission (antenna temperature 180 + or - 36 K) surrounded by an area of strong narrow-line 120-K emission is found. It is suggested that the plateau emission originates in a wide thin sheet after a shock wave, with H2 density at least 10 to the 6th/cu cm and gas kinetic temperature greater than 500 K.

Buhl, D.↗

A DSMC Surface Chemistry Model for Carbon-Based Ablators

A detailed molecular surface chemistry model for the DSMC (Direct Simulation Monte Carlo) method is proposed and implemented into the SPARTA (Stochastic PArallel Rarefied-gas Time-accurate Analyzer) DSMC solver. Molchanova et al. constructed a molecular model for surface recombination in DSMC that includes different surface processes (adsorption, desoprtion, Eley-Rideal and Langmuir-Hinshelwood). All surface processes can be divided into two groups: surface mechanisms, which involve only the particle adsorbed by the surface (desorption and Langmuir-Hinshelwood), and impact mechanisms, which also involve gas-phase particles (adsorption, Eley-Rideal). Using a similar approach, the 14-reaction kinetic model of oxygen-carbon interaction suggested by Zhlukhtov and Abe, as well as more recent models by Alba et al., Poovathinghal et al., and a new model developed in the scope of this work, are implemented in SPARTA. The computational results for the different oxidation models are compared with experimental results from Murray et al. (oxidation of a vitreous carbon surface due to a hyperthermal beam of O and O2), with a particular focus on fluxes, angular and Time-Of-Flight distributions of scattered particles.

oxidation↗

Study of Ring Current Dynamics During Geomagnetic Storms

This research program considered modeling the dynamical evolution of the ring current during several geomagnetic storms. The first year (6/01/1997-5/31/1998) of this successful collaborative research between the University of New Hampshire (UNH) and the University of California Los Angeles (UCLA) was supported by NASA grant NAG5-4680. The second and third years (6/01/1998-5/31/2000) were funded at UNH under NASA grant NAG5-7368. Research work at UNH concentrated on further development of a kinetic model to treat all of the important physical processes that affect the ring current ion population during storm conditions. This model was applied to simulate ring current development during several International Solar-Terrestrial Physics (ISTP) events, and the results were directly compared to satellite observations. A brief description of our major accomplishments and a list of the publications and presentations resulting from this effort are given.

Jordanova, Vania K.↗