Search NASA⌕ Search

SEARCH · Search NASA

Results for “kinetic methods”

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 361 records · Page 20
All sourcesNASA NTRSDOE OSTI
NASA NTRS1979

Matrix perturbation for analytical model improvement

A matrix perturbation method is proposed to calculate the Jacobian matrix and to compute the new eigendata for the parameter estimation procedure. The advantages of the method are the applicability to large complex structures without knowing the analytical expressions for the mass and stiffness matrices, and a cost effective approach for the re-computation of the eigendata. This method also allows the use of other measurements such as modal forces, kinetic energy distribution, and strain energy distributions in the estimation procedure. A realistic sample problem is presented to demonstrate the effectiveness of the proposed method.

Chen, J. C.↗
NASA NTRS1982

Kinetics and mechanism of the oxidation process of two-component Fe-Al alloys

The oxidation process of two-component Fe-Al alloys containing up to 7.2% Al and from 18 to 30% Al was studied. Kinetic measurements were conducted using the isothermal gravimetric method in the range of 1073-1223 K and 1073-1373 K for 50 hours. The methods used in studies of the mechanism of oxidation included: X-ray microanalysis, X-ray structural analysis, metallographic analysis and marker tests.

Przewlocka, H.↗
NASA NTRS2002

Semiclassical Calculation of Reaction Rate Constants for Homolytical Dissociations

There is growing interest in extending organometallic chemical vapor deposition (OMCVD) to III-V materials that exhibit large thermal decomposition at their optimum growth temperature, such as indium nitride. The group III nitrides are candidate materials for light-emitting diodes and semiconductor lasers operating into the blue and ultraviolet regions. To overcome decomposition of the deposited compound, the reaction must be conducted at high pressures, which causes problems of uniformity. Microgravity may provide the venue for maintaining conditions of laminar flow under high pressure. Since the selection of optimized parameters becomes crucial when performing experiments in microgravity, efforts are presently geared to the development of computational OMCVD models that will couple the reactor fluid dynamics with its chemical kinetics. In the present study, we developed a method to calculate reaction rate constants for the homolytic dissociation of III-V compounds for modeling OMCVD. The method is validated by comparing calculations with experimental reaction rate constants.

Beatriz H. Cardelino↗
NASA NTRS2006

BLIMPK/Streamline Surface Catalytic Heating Predictions on the Space Shuttle Orbiter

This paper describes the results of an analysis of localized catalytic heating effects to the U.S. Space Shuttle Orbiter Thermal Protection System (TPS). The analysis applies to the High-temperature Reusable Surface Insulation (HRSI) on the lower fuselage and wing acreage, as well as the critical Reinforced Carbon-Carbon on the nose cap, chin panel and the wing leading edge. The object of the analysis was to use a modified two-layer approach to predict the catalytic heating effects on the Orbiter windward HRSI tile acreage, nose cap, and wing leading edge assuming localized highly catalytic or fully catalytic surfaces. The method incorporated the Boundary Layer Integral Matrix Procedure Kinetic (BLIMPK) code with streamline inputs from viscous Navier-Stokes solutions to produce heating rates for localized fully catalytic and highly catalytic surfaces as well as for nominal partially catalytic surfaces (either Reinforced Carbon-Carbon or Reaction Cured Glass) with temperature-dependent recombination coefficients. The highly catalytic heating results showed very good correlation with Orbiter Experiments STS-2, -3, and -5 centerline and STS-5 wing flight data for the HRSI tiles. Recommended catalytic heating factors were generated for use in future Shuttle missions in the event of quick-time analysis of damaged or repaired TPS areas during atmospheric reentry. The catalytic factors are presented along the streamlines as well as a function of stagnation enthalpy so they can be used for arbitrary trajectories.

Marichalar, Jeremiah J.↗
NASA NTRS2014

Evaluation of Human and Anthropomorphic Test Device Finite Element Models under Spaceflight Loading Conditions

In an effort to develop occupant protection standards for future multipurpose crew vehicles, the National Aeronautics and Space Administration (NASA) has looked to evaluate the test device for human occupant restraint with the modification kit (THOR‐K) anthropomorphic test device (ATD) in relevant impact test scenarios. With the allowance and support of the National Highway Traffic Safety Administration, NASA has performed a series of sled impact tests on the latest developed THOR‐K ATD. These tests were performed to match test conditions from human volunteer data previously collected by the U.S. Air Force. The objective of this study was to evaluate the THOR‐K finite element (FE) model and the Total HUman Model for Safety (THUMS) FE model with respect to the tests performed. These models were evaluated in spinal and frontal impacts against kinematic and kinetic data recorded in ATD and human testing. Methods: The FE simulations were developed based on recorded pretest ATD/human position and sled acceleration pulses measured during testing. Predicted responses by both human and ATD models were compared to test data recorded under the same impact conditions. The kinematic responses of the models were quantitatively evaluated using the ISO‐metric curve rating system. In addition, ATD injury criteria and human stress/strain data were calculated to evaluate the risk of injury predicted by the ATD and human model, respectively. Results: Preliminary results show well‐correlated response between both FE models and their physical counterparts. In addition, predicted ATD injury criteria and human model stress/strain values are shown to positively relate. Kinematic comparison between human and ATD models indicates promising biofidelic response, although a slightly stiffer response is observed within the ATD. Conclusion: As a compliment to ATD testing, numerical simulation provides efficient means to assess vehicle safety throughout the design process and further improve the design of physical ATDs. The assessment of the THOR‐K and THUMS FE models in a spaceflight testing condition is an essential first step to implementing these models in the computational evaluation of spacecraft occupant safety. Promising results suggest future use of these models in the aerospace field.

Putnam, Jacob P.↗
NASA NTRSDate not supplied

Direct Simulation Monte Carlo Studies of the Gas Sampling for the VATMOS-SR Mission Concept

VATMOS-SR (Venus ATMOSpheric - Sample Return) is a small spacecraft mission concept which would return a gas sample from the upper atmosphere of Venus to Earth for scientific analysis. This could be the first sample return mission for an extra-terrestrial atmosphere, and potentially the first sample return from an Earth-sized planet. The VATMOS-SR mission concept consists of a SmallSat atmospheric sampling probe (45 deg. sphere/cone geometry, <1 m diameter) that is designed to skim through the Venus upper atmosphere and acquire gas samples below the homopause altitude (around ~110 km altitude), where the different atmospheric gases are mixed. The velocity of the spacecraft where sampling would occur is expected to be between ~10.5 km/s and ~13.1 km/s, depending on the trajectory chosen. This presentation will discuss hypervelocity sampling in the upper atmosphere of Venus, with respect to the VATMOS-SR mission concept. VATMOS-SR would enable critical atmospheric measurements to form a full picture of how, why, and when Venus evolved to be so different from Earth and Mars. The abundances and isotopic compositions of volatile elements (such as N, C, S, O and the noble gases) in planetary atmospheres record volatile delivery during accretion, outgassing from planetary interiors, and atmospheric loss to space. Precise and accurate determinations of atmospheric volatile signatures are the key to understanding the origins and geodynamical evolution of Venus compared to the other terrestrial planets. Hypersonic sampling poses unique technical and scientific challenges. To ensure it is possible to relate the composition of the sampled gases to the free stream atmospheric composition, large-scale numerical simulations are employed to model the flow through the VATMOS-SR sampling system. In particular, an emphasis is placed on quantifying noble gas isotopic fractionation that occurs during the sample acquisition and transfer process in order to determine how measured isotopic ratios of noble gases in the sample compare to actual isotopic ratios in the Venusian atmosphere. The Direct Simulation Monte Carlo (DSMC) code SPARTA, an open source software package developed by Sandia National Laboratories, is used in this work. SPARTA, based on Bird’s DSMC method, is a molecular-level gas-kinetic technique. As SPARTA is able to model hypervelocity reacting flows in strong chemical and thermal non-equilibrium, this software package is well suited to determine relevant flow properties for the VATMOS-SR mission concept, and to numerically quantify the expected level of elemental and/or isotopic fractionation in the sample acquired by VATMOS-SR. This presentation will show results from 3D simulations correlating the noble gas isotopic fractionation in the gas acquired at hypervelocity speeds to its ambient atmosphere value. In particular, emphasis will be placed at Xenon isotopes of masses 128 and 130, as precise measurements of that ratio would yield a comparison to Earth’s atmosphere. Additionally, sensitivity studies that quantify the uncertainties due to the freestream parameters as well as the modeling parameters will be performed.

direct simulation Monte Carlo↗
NASA NTRSDate not supplied

The VATMOS-SR Mission Concept: DSMC Studies of the Gas Sampling

VATMOS-SR (Venus ATMOSpheric - Sample Return) is a small spacecraft mission concept that would return a gas sample from the upper atmosphere of Venus to Earth for scientific analysis. This could be the first sample return mission for an extra-terrestrial atmosphere and potentially the first sample return from an Earth-sized planet. The VATMOS-SR mission concept consists of a SmallSat atmospheric sampling probe (45 deg. sphere/cone geometry, <1 m diameter) that is designed to skim through the Venus upper atmosphere and acquire gas samples below the homopause altitude (around ~110 km altitude), where the different atmospheric gases are mixed. The velocity of the spacecraft where sampling would occur is expected to be between ~10.5 km/s and ~13.1 km/s, depending on the trajectory chosen. This presentation will discuss hypervelocity sampling in the upper atmosphere of Venus with respect to the VATMOS-SR mission concept. VATMOS-SR would enable critical atmospheric measurements to form a complete picture of how, why, and when Venus evolved to be so different from Earth and Mars. The abundances and isotopic compositions of volatile elements (such as N, C, S, O, and the noble gases) in planetary atmospheres record volatile delivery during accretion, outgassing from planetary interiors, and atmospheric loss to space. Precise and accurate determinations of volatile atmospheric signatures are the key to understanding the origins and geodynamical evolution of Venus compared to the other terrestrial planets. Hypersonic sampling poses unique technical and scientific challenges. To ensure it is possible to relate the composition of the sampled gases to the free stream atmospheric composition, large-scale numerical simulations are employed to model the flow through the VATMOS-SR sampling system. In particular, an emphasis is placed on quantifying noble gas isotopic fractionation that occurs during the sample acquisition and transfer process to determine how measured isotopic ratios of noble gases in the sample compared to actual isotopic ratios in the Venusian atmosphere. The Direct Simulation Monte Carlo (DSMC) code SPARTA, an open-source software package developed by Sandia National Laboratories, is used in this work. SPARTA, based on Bird’s DSMC method, is a molecular-level gas-kinetic technique. As SPARTA is able to model hypervelocity reacting flows in strong chemical and thermal non-equilibrium, this software package is well suited to determine relevant flow properties for the VATMOS-SR mission concept and to numerically quantify the expected level of elemental and/or isotopic fractionation in the sample acquired by VATMOS-SR. This presentation will show results from 3D simulations correlating the noble gas isotopic fractionation in the gas acquired at hypervelocity speeds to its ambient atmosphere value. In particular, emphasis will be placed on Xenon isotopes of masses 128 and 130, as precise measurements of that ratio would yield comparison to Earth’s atmosphere. Additionally, sensitivity studies that quantify the uncertainties due to the freestream parameters, as well as the modeling parameters, will be performed.

direct simulation Monte Carlo↗
DOE OSTI2024

Accelerating particle-in-cell kinetic plasma simulations via reduced-order modeling of space-charge dynamics using dynamic mode decomposition

We present a data-driven reduced-order modeling of the space-charge dynamics for electromagnetic particle-in-cell (EMPIC) plasma simulations based on dynamic mode decomposition (DMD). The dynamics of the charged particles in kinetic plasma simulations such as EMPIC is manifested through the plasma current density defined along the edges of the spatial mesh. We showcase the efficacy of DMD in modeling the time evolution of current density through a low-dimensional feature space. Not only do such DMD based predictive reduced-order models help accelerate EMPIC simulations, they also have the potential to facilitate investigative analysis and control applications. Here, we demonstrate the proposed DMD-EMPIC scheme for reduced-order modeling of current density and speedup in EMPIC simulations involving electron beam under the influence of magnetic field, virtual cathode oscillations, and backward wave oscillator.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1981

A study of the viscous and nonadiabatic flow in radial turbines

A method for analyzing the viscous nonadiabatic flow within turbomachine rotors is presented. The field analysis is based upon the numerical integration of the incompressible Navier-Stokes equations together with the energy equation over the rotors blade-to-blade stream channels. The numerical code used to solve the governing equations employs a nonorthogonal boundary fitted coordinate system that suits the most complicated blade geometries. Effects of turbulence are modeled with two equations; one expressing the development of the turbulence kinetic energy and the other its dissipation rate. The method of analysis is applied to a radial inflow turbine. The solution obtained indicates the severity of the complex interaction mechanism that occurs between different flow regimes (i.e., boundary layers, recirculating eddies, separation zones, etc.). Comparison with nonviscous flow solutions tend to justify strongly the inadequacy of using the latter with standard boundary layer techniques to obtain viscous flow details within turbomachine rotors. Capabilities and limitations of the present method of analysis are discussed.

Khalil, I.↗
NASA NTRS1986

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.↗
DOE OSTI2024

On the prediction of pressure effects for the combination kinetics of two alkyl radicals with the geometric mean rule

Here, alkyl combination reactions, along with their reverse alkane unimolecular decomposition reactions, are widely present in pyrolysis and combustion systems. These reactions occur mainly at temperatures at which their rate coefficients depend on pressure, and only a few studies have explored this dependence for large size alkyls. The accurate calculation of these kinetic data is limited by the size of the alkyl fragments involved, since the theoretical methods used are computationally expensive. Therefore, little reliable data is available for the modeling of real fuels alkane decomposition kinetics. In this paper, we propose an alternative approach that allows an accurate estimation of pressure-dependent kinetic rates for the combination of long-chain alkyl radicals, based on tabulated reaction models. The rates of combination of linear alkyls radicals forming ethane to pentane were calculated using the CASPT2-F12/cc-pVDZ-F12 level within VRC-TST theory combined with master-equation simulations. The computed high-pressure limit rate coefficients showed that n-propyl could be a minimum-size representative structure for larger n-alkyl chain, but this approach fails to capture the pressure dependence of the reaction that remains size-dependent for larger alkyls. The Geometric Mean Rule (GMR), a simple rule that links the self and cross-reactions of three alkyl radicals, is used to address this issue. Our theoretical calculations show for the first time that the GMR applies not only at the high-pressure limit, but also for every pressure explored. The GMR is shown to be accurate, with a maximum 22 % deviation in the falloff region for the estimation of rates for small alkyl combinations. Using this method, we propose a tabulated reaction rate rule for the combination of alkyl radicals up to C5+C5. The use of computed rate coefficients with the GMR approximation opens ways to develop pressure-dependent reaction rate rules for large n-alkyl combinations, based on accurate and computationally expensive theoretical calculations on smaller alkyl fragments.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2026

A Perspective on Traditional and Data Driven Electrochemical Modeling and Analysis

To understand the behavior of electrochemical systems, we need to reduce the dimensionality of the measured current-voltage-time (I-V-t) data by fitting models, thus enabling us to analyze and compare the governing physics. Traditionally, the process for this is an 'expert first' approach: defining the model and its explicit assumptions based on inductive reasoning or empirical observation, fitting small portions of the I-V-t data where assumptions are most valid or carefully designing experiments to enforce key assumptions, and then interpreting the model parameters. However, modern data-driven methods enable a new paradigm: a 'data first' approach, where the latent behaviors governing the system's measured response are identified directly using machine-learning models that optimize both model structure and parameters from the I-V-t data, guaranteeing that the learned model explains as much of the observed system response as possible. After model identification, the model can then be interrogated by an expert to connect observed behaviors with underlying physics. This talk will review several different types of electrochemical analysis (electrochemical impedance, differential voltage-capacity, electrochemical kinetics) and compare the traditional and data-driven methods for analyzing the data.

42 ENGINEERING↗
NASA NTRS1990

Reduced and simplified chemical kinetics for air dissociation using Computational Singular Perturbation

The Computational Singular Perturbation CSP methods is employed (1) in the modeling of a homogeneous isothermal reacting system and (2) in the numerical simulation of the chemical reactions in a hypersonic flowfield. Reduced and simplified mechanisms are constructed. The solutions obtained on the basis of these approximate mechanisms are shown to be in very good agreement with the exact solution based on the full mechanism. Physically meaningful approximations are derived. It is demonstrated that the deduction of these approximations from CSP is independent of the complexity of the problem and requires no intuition or experience in chemical kinetics.

Goussis, D. A.↗
NASA NTRS2005

Summary of Simplified Two Time Step Method for Calculating Combustion Rates and Nitrogen Oxide Emissions for Hydrogen/Air and Hydrogen/Oxygen

A simplified single rate expression for hydrogen combustion and nitrogen oxide production was developed. Detailed kinetics are predicted for the chemical kinetic times using the complete chemical mechanism over the entire operating space. These times are then correlated to the reactor conditions using an exponential fit. Simple first order reaction expressions are then used to find the conversion in the reactor. The method uses a two time step kinetic scheme. The first time averaged step is used at the initial times with smaller water concentrations. This gives the average chemical kinetic time as a function of initial overall fuel air ratio, temperature, and pressure. The second instantaneous step is used at higher water concentrations (greater than l x 10(exp -20)) moles per cc) in the mixture which gives the chemical kinetic time as a function of the instantaneous fuel and water mole concentrations, pressure and temperature (T(sub 4)). The simple correlations are then compared to the turbulent mixing times to determine the limiting properties of the reaction. The NASA Glenn GLSENS kinetics code calculates the reaction rates and rate constants for each species in a kinetic scheme for finite kinetic rates. These reaction rates are used to calculate the necessary chemical kinetic times. This time is regressed over the complete initial conditions using the Excel regression routine. Chemical kinetic time equations for H2 and NOx are obtained for H2/Air fuel and for H2/O2. A similar correlation is also developed using data from NASA's Chemical Equilibrium Applications (CEA) code to determine the equilibrium temperature (T(sub 4)) as a function of overall fuel/air ratio, pressure and initial temperature (T(sub 3)). High values of the regression coefficient R squared are obtained.

Marek, C. John↗
DOE OSTI2024

Factors that Impact the Redox Kinetics Associated with Sr1-xAxFeO3 (A = Ca, Ba) in Air Separations Technologies

Perovskite oxides, including SrFeO3-δ, have garnered significant interest for their impressive oxygen storage thermodynamics and kinetics. While traditional oxygen production is performed through cryogenic separation methods at the industrial scale, oxygen carrier-based separations have offered an attractive alternative at the modular scale. Perovskites are considered ideal oxygen carriers, with their structural flexibility minimizing both the thermodynamic barrier and the structural strain related to the oxygen uptake/release cycle. Our latest work has identified a few major factors that directly impact the oxygen release thermodynamics and kinetics for Sr1-xAxFeO3-δ (A = Ca, Ba) materials including elemental composition, crystallographic symmetry, surface area, and pretreatment conditions. This presentation will discuss these factors using Sr1¬-xCaxFeO3 and Ba1-xSrxFeO3 to illustrate their individual and collective impacts on the design of high-performance oxygen carriers when operating at 300-700 °C. Additionally, we will discuss how each of these factors can be leveraged for the numerous catalytic applications also dependent on the redox properties of mixed metal oxides.

Popczun, Eric↗
NASA NTRS1978

Pulsational stabilities of a star in thermal imbalance - Comparison between the methods

The stability coefficients for quasi-adiabatic pulsations for a model in thermal imbalance are evaluated using the dynamical energy (DE) approach, the total (kinetic plus potential) energy (TE) approach, and the small amplitude (SA) approaches. From a comparison among the methods, it is found that there can exist two distinct stability coefficients under conditions of thermal imbalance as pointed out by Demaret. It is shown that both the TE approaches lead to one stability coefficient, while both the SA approaches lead to another coefficient. The coefficient obtained through the energy approaches is identified as the one which determines the stability of the velocity amplitudes. For a prenova model with a thin hydrogen-burning shell in thermal imbalance, several radial modes are found to be unstable both for radial displacements and for velocity amplitudes. However, a new kind of pulsational instability also appears, viz., while the radial displacements are unstable, the velocity amplitudes may be stabilized through the thermal imbalance terms.

Vemury, S. K.↗
DOE OSTI2025

On High-Order/Low-Order and Micro-Macro Methods for Implicit Time-Stepping of the BGK Model

In this paper, a high-order/low-order (HOLO) method is combined with a micro-macro (MM) decomposition to accelerate iterative solvers in fully implicit time-stepping of the Bhatnagar–Gross–Krook (BGK) equation for gas dynamics. The MM formulation represents a kinetic distribution as the sum of a local Maxwellian and a perturbation. In highly collisional regimes, the perturbation away from initial and boundary layers is small and can be compressed to reduce the overall storage cost of the distribution. The convergence behavior of the MM methods, the usual HOLO method, and the standard source iteration method is analyzed on a linear BGK model. Both the HOLO and MM methods are implemented using a discontinuous Galerkin (DG) discretization in phase space, which naturally preserves the consistency between high- and low-order models required by the HOLO approach. Furthermore, the accuracy and performance of these methods are compared on the Sod shock tube problem and a sudden wall heating boundary layer problem. Overall, the results demonstrate the robustness of the MM and HOLO approaches and illustrate the compression benefits enabled by the MM formulation when the kinetic distribution is near equilibrium.

BGK model↗