Search NASA⌕ Search

SEARCH · Search NASA

Results for “reaction kinetics and dynamics”

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 55 records · Page 3

On Kinetics Modeling of Vibrational Energy Transfer

Two models of vibrational energy exchange are compared at equilibrium to the elementary vibrational exchange reaction for a binary mixture. The first model, non-linear in the species vibrational energies, was derived by Schwartz, Slawsky, and Herzfeld (SSH) by considering the detailed kinetics of vibrational energy levels. This model recovers the result demanded at equilibrium by the elementary reaction. The second model is more recent, and is gaining use in certain areas of computational fluid dynamics. This model, linear in the species vibrational energies, is shown not to recover the required equilibrium result. Further, this more recent model is inconsistent with its suggested rate constants in that those rate constants were inferred from measurements by using the SSH model to reduce the data. The non-linear versus linear nature of these two models can lead to significant differences in vibrational energy coupling. Use of the contemporary model may lead to significant misconceptions, especially when integrated in computer codes considering multiple energy coupling mechanisms.

Gilmore, John O.↗

Dynamic mechanical characterization of cure of a polyimide-graphite fiber composite (PMR 15/Celion 6000)

Dynamic mechanical studies of thermosetting PMR 15 polyimide/Celion 6000 composites were used to characterize cure behavior. Variation in composite shear modulus with cure time for laminates partially cured under 3.45 MPa pressure was compared with that of laminates fully cured under pressure to select a partially cured system whose behavior corresponded with that of the pressure-cured composites. An empirical kinetic model was developed that relates relative composite shear moduli to cure time and temperature. This model, coupled with standard statistical techniques, was used in the determination of an overall activation energy, E = 145 kj/mole (35 kcal/mole) and order of reaction, n, for the cure reaction(s). The empirical relationship defining dynamic modulus as a function of time and temperature remained linear well into the glassy region. Changes in dynamic glass-transition temperature and in the breadth of the transition peaks with cure time are discussed. Shear modulus and damping were found to be more sensitive measures of relative extent of cure than Fourier transform infrared spectroscopy.

Hurwitz, F. I.↗

A Role of the Reaction Kernel in Propagation and Stabilization of Edge Diffusion Flames of C1-C3 Hydrocarbons

Diffusion flame stabilization is of essential importance in both Earth-bound combustion systems and spacecraft fire safety. Local extinction, re-ignition, and propagation processes may occur as a result of interactions between the flame zone and vortices or fire-extinguishing agents. By using a computational fluid dynamics code with a detailed chemistry model for methane combustion, the authors have revealed the chemical kinetic structure of the stabilizing region of both jet and flat-plate diffusion flames, predicted the flame stability limit, and proposed diffusion flame attachment and detachment mechanisms in normal and microgravity. Because of the unique geometry of the edge of diffusion flames, radical back-diffusion against the oxygen-rich entrainment dramatically enhanced chain reactions, thus forming a peak reactivity spot, i.e., reaction kernel, responsible for flame holding. The new results have been obtained for the edge diffusion flame propagation and attached flame structure using various C1-C3 hydrocarbons.

Takahashi, Fumiaki↗

Numerical study of shock-induced combustion in methane-air mixtures

The shock-induced combustion of methane-air mixtures in hypersonic flows is investigated using a new reaction mechanism consisting of 19 reacting species and 52 elementary reactions. This reduced model is derived from a full kinetic mechanism via the Detailed Reduction technique. Zero-dimensional computations of several shock-tube experiments are presented first. The reaction mechanism is then combined with a fully implicit Navier-Stokes computational fluid dynamics (CFD) code to conduct numerical simulations of two-dimensional and axisymmetric shock-induced combustion experiments of stoichiometric methane-air mixtures at a Mach number of M = 6.61. Applications to the ram accelerator concept are also presented.

Yungster, Shaye↗

Incorporation of a Chemical Equilibrium Equation of State into LOCI-Chem

Renewed interest in development of advanced high-speed transport, reentry vehicles and propulsion systems has led to a resurgence of research into high speed aerodynamics. As this flow regime is typically dominated by hot reacting gaseous flow, efficient models for the characteristic chemical activity are necessary for accurate and cost effective analysis and design of aerodynamic vehicles that transit this regime. The LOCI-Chem code recently developed by Ed Luke at Mississippi State University for NASA/MSFC and used by NASA/MSFC and SSC represents an important step in providing an accurate, efficient computational tool for the simulation of reacting flows through the use of finite-rate kinetics [3]. Finite rate chemistry however, requires the solution of an additional N-1 species mass conservation equations with source terms involving reaction kinetics that are not fully understood. In the equilibrium limit, where the reaction rates approach infinity, these equations become very stiff. Through the use of the assumption of local chemical equilibrium the set of governing equations is reduced back to the usual gas dynamic equations, and thus requires less computation, while still allowing for the inclusion of reacting flow phenomenology. The incorporation of a chemical equilibrium equation of state module into the LOCI-Chem code was the primary objective of the current research. The major goals of the project were: (1) the development of a chemical equilibrium composition solver, and (2) the incorporation of chemical equilibrium solver into LOCI-Chem. Due to time and resource constraints, code optimization was not considered unless it was important to the proper functioning of the code.

Cox, Carey F.↗

A Multi-Physics Study on High-Specific Power Li-O2 Batteries for Electric Aircraft

Commercialization of lithium-air batteries faces many challenges, such as electrolyte decomposition, short cycle life, low energy and power density, etc. However, commercialization of Li-O2 batteries for aeronautics is much more challenging due to additional safety constraints on cyclability and performance (high specific power and specific energy). For this presentation, we will discuss inter-related aspects of physics-based modeling of a pack: cell and battery model calibration. In addition, we will evaluate and present optimal battery designs for high discharge current density, high discharge time, and low battery mass using simulation-based optimization.The Finite Element Model (FEM) used to simulate a Li-O2 cell is based on the work of Bevara [1]. The different aspects of the model are based on: porous electrode theory and concentrated electrolyte theory; quantum tunneling model for the resistance of conformal layer of discharge product (Li2O2) [1]; Butler-Volmer kinetics for electrochemical reaction; Fick's diffusion for oxygen transport; and an oxygen dissolution model is applied at the air/electrolyte interface [2]. The electrolyte properties such as ion conductivity, ion diffusion, oxygen diffusion, and mass density of the electrolyte were taken from Molecular Dynamics (MD) simulations [3]; while the other model parameters, which includes mass of cell components, were calibrated to match experiments at high discharge current densities. The cell mass includes the anode, cathode, separator, electrolyte, and other components (such as current collector). This calibrated model is used to perform parametric studies on cathode thickness, porosity, tortuosity, carbon particle size, electrolyte transport and material properties, partial pressure of oxygen, discharge time, and discharge current density to study optimal designs for high specific power and energy. References:1. Bevara, V. & Andrei, P. (2014), J. Electrochem. Soc. 161 (14), A2068-A2079.2.Mehta, M. & Andrei, P. (2015), J. Power Sources. 286, 299-308.3.Liyana-Arachchi, T.; Haskins, J.; Burke, C.; Diederichsen, K.; McCloskey, B.; & Lawson, J. (2018), J. Phys. Chem. B. 122 (36), 8548 - 8559.4.Choi, W.; Kikumoto, H.; Choudhary, R. & Ooka, R. (2018), Applied Energy, 209, 306-321.

Mehta, Mohit↗

Surrogate Fuel Chemistry from ReaxFF Molecular Dynamics

Reactive molecular dynamics with ReaxFF simulations are used to study the chemistry of four component jet fuel mixtures of n-dodecane, iso-octane, n-propylbenzene, and mesitylene. ReaxFF does not require knowledge of the chemistry a priori providing insight into reactions otherwise not seen in established chemical kinetic mechanisms. The results show that a difference in the local chemical environment can change the fuel decomposition pathways across the surrogate mixtures where the same fuel species are used. These simulations highlight the need to consider cross-chemical interactions in chemical kinetic mechanism development and the importance of resonance stabilized polyaromatic hydrocarbon reactions for soot growth and formation. A method for predicting the threshold sooting index from the carbon species formed at the end of the simulations is established. The results of the simulations are consistent with measurements of the threshold sooting index with the predicted values within 20% of the measured data

Surrogate Fuels↗

Computation of thermo-chemical nonequilibrium Martian atmospheric entry flows

The aerothermodynamics of a Martian atmospheric entry trajectory are studied. A chemical kinetics model for the CO2-N2 mixture found in the Martian atmosphere is implemented in a two-dimensional computational fluid dynamics method. The gas mixture is allowed to be in thermo-chemical nonequilibrium and to be composed of eight species characterized by two temperatures. The CO2 reaction model is validated by comparing computations to experimental results. Results for typical Martian entry conditions are presented and limitations of the current thermophysical model are discussed.

Candler, Graham V.↗

Time-dependent Computational Studies of Premixed Flames in Microgravity

This report describes the research performed at the Center for Reactive Flow and Dynamical Systems in the Laboratory for Computational Physics and Fluid Dynamics, at the Naval Research Laboratory, in support of NASA Microgravity Science and Applications Program. The primary focus of this research is on investigating fundamental questions concerning the propagation and extinction of premixed flames in earth gravity and in microgravity environments. Our approach is to use detailed time-dependent, multispecies, numerical models as tools to simulate flames in different gravity environments. The models include a detailed chemical kinetics mechanism consisting of elementary reactions among the eight reactive species involved in hydrogen combustion, coupled to algorithms for convection, thermal conduction, viscosity, molecular and thermal diffusion, and external forces. The external force, gravity, can be put in any direction relative to flame propagation and can have a range of values. Recently more advanced wall boundary conditions such as isothermal and no-slip have been added to the model. This enables the simulation of flames propagating in more practical systems than before. We have used the numerical simulations to investigate the effects of heat losses and buoyancy forces on the structure and stability of flames, to help resolve fundamental questions on the existence of flammability limits when there are no external losses or buoyancy forces in the system, to understand the interaction between the various processes leading to flame instabilities and extinguishment, and to study the dynamics of cell formation and splitting. Our studies have been able to bring out the differences between upward- and downward-propagating flames and predict the zero-gravity behavior of these flames. The simulations have also highlighted the dominant role of wall heat losses in the case of downward-propagating flames. The simulations have been able to qualitatively predict the formation of multiple cells and the cessation of cell-splitting. Our studies have also shown that some flames in a microgravity environment can be extinguished due to a chemical instability and without any external losses. However, further simulations are needed to more completely understand upward-propagating and zero-gravity flames as well as to understand the potential effect of radiative heat losses.

Kailasanath, K.↗

A Theoretical Investigation of Oxidation Efficiency of a Volatile Removal Assembly Reactor Under Microgravity Conditions

Volatile Removal Assembly (VRA) is a subsystem of the Closed Environment Life Support System (CELSS) installed in the International Space Station. It is used for removing contaminants (volatile organics) in the wastewater produced by the space station crews. The major contaminants are formic acid, ethanol, and propylene glycol. The VRA contains a slim packbed reactor (3.5 cm diameter and four 28 cm long tubes in series) to perform catalyst oxidation of wastewater at elevated pressure and temperature under microgravity conditions. In the reactor, the contaminants are burned with oxygen gas (O2) to form water and carbon dioxide (CO2) that dissolves in the water stream. Optimal design of the reactor requires a thorough understanding about how the reactor performs under microgravity conditions. The objective of this study was to develop a mathematical model to interpret experimental data obtained from normal and microgravity conditions, and to predict the performance of VRA reactor under microgravity conditions. Catalyst oxidation kinetics and the total oxygen-water contact area control the efficiency of catalyst oxidation for mass transfer, which depends on oxygen gas holdup and distribution in the reactor. The process involves bubbly flow in porous media with chemical reactions in microgravity environment. This presents a unique problem in fluid dynamics that has not been studied. Guo et al. (2004) developed a mathematical model that predicts oxygen holdup in the VRA reactor. No mathematical model has been found in the literature that can be used to predict the efficiency of catalyst oxidation under microgravity conditions.

Guo, Boyun↗

Predictive Modeling of Fast-Curing Thermosets in Nozzle-Based Extrusion

This work presents an approach to modeling the dynamic spreading and curing behavior of thermosets in nozzle-based extrusions. Thermosets cover a wide range of materials, some of which permit low-temperature processing with subsequent high-temperature and high-strength working properties. Extruding thermosets may overcome the limited working temperatures and strengths of conventional thermoplastic materials used in additive manufacturing. This project aims to produce technology for the fabrication of thermoset-based structures leveraging advances made in nozzle-based extrusion, such as fused deposition modeling (FDM), material jetting, and direct writing. Understanding the synergistic interactions between spreading and fast curing of extruded thermosetting materials will provide essential insights for applications that require accurate dimensional controls, such as additive manufacturing [1], [2] and centrifugal coating/forming [3]. Two types of thermally curing thermosets -- one being a soft silicone (Ecoflex 0050) and the other being a toughened epoxy (G/Flex) -- served as the test materials in this work to obtain models for cure kinetics and viscosity. The developed models align with extensive measurements made with differential scanning calorimetry (DSC) and rheology. DSC monitors the change in the heat of reaction, which reflects the rate and degree of cure at different crosslinking stages. Rheology measures the change in complex viscosity, shear moduli, yield stress, and other properties dictated by chemical composition. By combining DSC and rheological measurements, it is possible to establish a set of models profiling the cure kinetics and chemorheology without prior knowledge of chemical composition, which is usually necessary for sophisticated mechanistic modeling. In this work, we conducted both isothermal and dynamic measurements with both DSC and rheology. With the developed models, numerical simulations yielded predictions of diameter and height of droplets, along with width and height of extruded lines cured at varied temperatures. Experimental results carried out on a goniometric platform and a nozzle-based 3D printer showed agreement with the numerical simulations. Finally, this presentation will show how the models are adaptable to the planning of tool paths and designs in additive manufacturing.

Xie, Jingjin↗

Data Parallel Line Relaxation (DPLR) Code User Manual: Acadia - Version 4.01.1

Data-Parallel Line Relaxation (DPLR) code is a computational fluid dynamic (CFD) solver that was developed at NASA Ames Research Center to help mission support teams generate high-value predictive solutions for hypersonic flow field problems. The DPLR Code Package is an MPI-based, parallel, full three-dimensional Navier-Stokes CFD solver with generalized models for finite-rate reaction kinetics, thermal and chemical non-equilibrium, accurate high-temperature transport coefficients, and ionized flow physics incorporated into the code. DPLR also includes a large selection of generalized realistic surface boundary conditions and links to enable loose coupling with external thermal protection system (TPS) material response and shock layer radiation codes.

Wright, Michael J.↗

Considerations for Femtosecond Laser Electronic Excitation Tagging in High-Speed Flows

Femtosecond laser electronic excitation tagging (FLEET) is an unseeded method for molecular tagging which offers valuable opportunities for measurement of high-speed (transonic, supersonic or hypersonic) flows. The unique nature of high-speed testing demands certain performance from FLEET such as satisfactory signal-to-noise ratio (SNR) at depressed static conditions (i.e., low temperatures, pressures and densities), wide dynamic range for velocity determination (especially single-shot), and measurements with acceptable accuracy and precision. This dissertation strives to evaluate FLEET in those regards and provide strategies to maximize the method's capabilities. A zero-dimensional kinetics model in nitrogen explains FLEET signal changes with pressure/density and/or temperature in terms of plasma-chemical reactions. Poorly known rate coefficients are tuned by comparing model output to measurements, with temporal agreement up to several hundred nanoseconds. Modeling reveals that initial signal peaks at reduced density because of slowed temporal evolution (and decay) of excited populations. Low temperatures enhance signal by enlarging cluster ion populations which contribute to excited species via electron-ion dissociative recombination. A purpose-built free jet facility provides experimental validation of the kinetics model and assesses FLEET velocimetry in low temperature and pressure/density conditions. Signal, lifetime, accuracy and precision results are obtained from unheated subsonic through Mach 4.0 operation of the facility, with best results noted. FLEET measurements of a sweeping jet (SWJ) actuator in compressible operation showcase its advantages in a highly unsteady jet containing subsonic through supersonic velocities. FLEET velocimetry is performed in the device's internal and external flow fields, with the latter compared to hot-wire anemometry. Internal measurements reveal the absence of shockwaves theorized to occur at high pressure ratios. Simultaneous qualitative measurements of compressible jet mixing are shown as a proof-of-concept. Overall, the work demonstrates that previous understanding of SWJ incompressible operation readily extends into the compressible realm. Practical aspects of performing FLEET velocimetry are detailed, along with strategies for improving measurement quality. Determination of a fundamental precision in nitrogen and air is attempted. Experiments show that increasing time delay and/or SNR improves velocimetry precision. A comparison of five camera systems indicates sensors with larger pixels capture higher SNR data and produce more precise results.

FLEET↗

Analytical ultrasonics for characterization of metallurgical microstructures and transformations

The application of contact (piezoelectric) and noncontact (laser generation and detection) ultrasonic techniques for dynamic investigation of precipitation hardening processes in aluminum alloys, as well as crystallization and phase transformation in rapidly solidified amorphous and microcrystalline alloys is discussed. From the variations of the sound velocity and attenuation the precipitation mechanism and kinetics were determined. In addition, a correlation was established between the observed changes in the velocity and attenuation and the mechanical properties of age-hardenable aluminum alloys. The behavior of the elastic moduli, determined ultrasonically, were found to be sensitive to relaxation, crystallization and phase decomposition phenomena in rapidly solidified metallic glasses. Analytical ultrasonics enables determination of the activation energies and growth parameters of the reactions. Therefrom theoretical models can be constructed to explain the changes in mechanical and physical properties upon heat treatment of glassy alloys. The composition dependence of the elastic moduli in amorphous Cu-Zr alloys was found to be related to the glass transition temperature, and consequently to the glass forming ability of these alloys. Dynamic ultrasonic analysis was found to be feasible for on-line, real-time, monitoring of metallurgical processes.

Rosen, M.↗

Eddies, streams, and convergence zones in turbulent flows

Recent studies of turbulent shear flows have shown that many of their important kinematical and dynamical properties can be more clearly understood by describing the flows in terms of individual events or streamline patterns. These events or flow regions are studied because they are associated with relatively large contributions to certain average properties of the flow, for example kinetic energy, Reynolds stress, or to particular processes in the flow, such as mixing and chemical reactions, which may be concentrated at locations where streamlines converge for fast chemical reactions (referred to as convergence or C regions), or in recirculating eddying regions for slow chemical reactions. The aim of this project was to use the numerical simulations to develop suitable criteria for defining these eddying or vortical zones. The C and streaming (S) zones were defined in order to define the whole flow field. It is concluded that homogeneous and sheared turbulent flow fields are made up of characteristic flow zones: eddy, C, and S zones. A set of objective criteria were found which describe regions in which the streamlines circulate, converge or diverge, and form high streams of high velocity flow.

Hunt, J. C. R.↗

Low-Temperature Oxidation Reactions and Cool Flames at Earth and Reduced Gravity

Non-isothermal studies of cool flames and low temperature oxidation reactions in unstirred closed vessels are complicated by the perturbing effects of natural convection at earth gravity. Buoyant convection due to self-heating during the course of slow reaction produces spatio-temporal variations in the thermal and thus specie concentration fields due to the Arrhenius temperature dependence of the reaction rates. Such complexities have never been quantitatively modeled and were the primary impetus for the development of CSTR's (continuously stirred tank reactors) 30 years ago. While CSTR's have been widely adopted since they offer the advantage of spatial uniformity in temperature and concentration, all gradients are necessarily destroyed along with any structure that may otherwise develop. Microgravity offers a unique environment where buoyant convection can be effectively minimized and the need for stirring eliminated. Moreover, eliminating buoyancy and the need for stirring eliminates complications associated with the induced hydrodynamic field whose influence on heat transport and hot spot formation, hence explosion limits, is not fully realized. The objective of this research is to quantitatively determine and understand the fundamental mechanisms that control the onset and evolution of low temperature reactions and cool flames in both static and flow reactors. Microgravity experiments will be conducted to obtain benchmark data on the structure (spatio-temporal temperature, concentration, flow fields), the dynamics of the chemical fronts, and the ignition diagrams (pressure vs. temperature). Ground-based experiments will be conducted to ascertain the role of buoyancy. Numerical simulations including detailed kinetics will be conducted and compared to experiment.

Pearlman, Howard↗

Evaporation kinetics of Mg2SiO4 crystals and melts from molecular dynamics simulations

Computer simulations based on the molecular dynamics (MD) technique were used to study the mechanisms and kinetics of free evaporation from crystalline and molten forsterite (i.e., Mg2SiO4) on an atomic level. The interatomic potential employed for these simulations reproduces the energetics of bonding in forsterite and in gas-phase MgO and SiO2 reasonably accurately. Results of the simulation include predicted evaporation rates, diffusion rates, and reaction mechanisms for Mg2SiO4(s or l) yields 2Mg(g) + 20(g) + SiO2(g).

Kubicki, J. D.↗

The topside ionosphere - A region of dynamic transition

The review article concentrates on dynamic processes at work in the topside ionosphere (between the F2 peak and about 3000 km) where the H ion dominates and ionic reactions can be neglected. The history of ionosphere and plasmasphere research using radio waves is reviewed. Low-speed and high-speed multispecies plasma ion flow is studied with various models (13-moment approximation, 5-moment approximation, kinetic models of the polar wind). Experimental observations of the plasmapause, results of vertical soundings of the topside, and global pole-to-pole distributions of ion composition, plasma temperature, and electron density are reviewed.

Banks, P. M.↗