Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 1,135 records · Page 63

Regression Analysis of Top of Descent Location for Idle-thrust Descents

In this paper, multiple regression analysis is used to model the top of descent (TOD) location of user-preferred descent trajectories computed by the flight management system (FMS) on over 1000 commercial flights into Melbourne, Australia. The independent variables cruise altitude, final altitude, cruise Mach, descent speed, wind, and engine type were also recorded or computed post-operations. Both first-order and second-order models are considered, where cross-validation, hypothesis testing, and additional analysis are used to compare models. This identifies the models that should give the smallest errors if used to predict TOD location for new data in the future. A model that is linear in TOD altitude, final altitude, descent speed, and wind gives an estimated standard deviation of 3.9 nmi for TOD location given the trajec- tory parameters, which means about 80% of predictions would have error less than 5 nmi in absolute value. This accuracy is better than demonstrated by other ground automation predictions using kinetic models. Furthermore, this approach would enable online learning of the model. Additional data or further knowl- edge of algorithms is necessary to conclude definitively that no second-order terms are appropriate. Possible applications of the linear model are described, including enabling arriving aircraft to fly optimized descents computed by the FMS even in congested airspace. In particular, a model for TOD location that is linear in the independent variables would enable decision support tool human-machine interfaces for which a kinetic approach would be computationally too slow.

trajectory prediction↗

Photoionization of seeded combustion products as a method of enhancing the efficiency of magnetohydrodynamic power generators

Here, in this study, we performed an experimental and computational investigation into the feasibility of utilizing photoionization to enhance the electrical conductivity of seeded oxy-fuel combustion products and improve the performance of magnetohydrodynamic (MHD) power generators. We applied a variety of optical and microwave diagnostics to study the ionization and recombination processes of potassium excited by an excimer laser in a high-velocity oxy-fuel free jet. Computational fluid dynamic (CFD) simulations were performed to model the thermophysical properties and species densities of the free jet. The CFD results were validated with position-dependent potassium concentration measurements. Electron recombination exponential lifetimes were measured through time-resolved microwave transmission. The experimental electron lifetimes were compared with lifetimes calculated from CFD-predicted species densities and literature recombination rates. It was determined that K + or O 2 are the most likely recombination partners for photoionized electrons. Time-resolved fluorescence measurements provided evidence of an ionization pathway involving a two-photon ionization of KOH . Finally, a zero-dimensional chemical kinetic model was developed to assess the fundamental viability of inducing a non-equilibrium electron population to provide a net energy return in combustion-driven MHD power generators. We determined that a high energy return is feasible for targeting electrode boundary layers with ultraviolet photoionization. We also found that photoionization could potentially lower the required temperature of the bulk gas flow.

20 FOSSIL-FUELED POWER PLANTS↗

Stochastic methods for aerosol chemistry: a compact molecular description of functionalization and fragmentation in the heterogeneous oxidation of squalane aerosol by OH radicals

The heterogeneous oxidation of organic aerosol by hydroxyl radicals (OH) can proceed through two general pathways: functionalization, in which oxygen functional groups are added to the carbon skeleton, and fragmentation, in which carbon–carbon bonds are broken, producing higher volatility, lower molecular weight products. An ongoing challenge is to develop a quantitative molecular description of these pathways that connects the oxidative evolution of the average aerosol properties (e.g. size and hygroscopicity) to the transformation of free radical intermediates. In order to investigate the underlying molecular mechanism of aerosol oxidation, a relatively compact kinetics model is developed for the heterogeneous oxidation of squalane particles by OH using free radical intermediates that convert reactive hydrogen sites into oxygen functional groups. Stochastic simulation techniques are used to compare calculated system properties over ten oxidation lifetimes with the same properties measured in experiment. The time-dependent average squalane aerosol mass, volume, density, carbon number distribution of scission products, and the average elemental composition are predicted using known rate coefficients. For functionalization, the calculations reveal that the distribution of alcohol and carbonyl groups is controlled primarily by the initial OH abstraction rate and to lesser extent by the branching ratio between secondary peroxy radical product channels. For fragmentation, the calculations reveal that the formation of activated alkoxy radicals with neighboring functional groups controls the molecular decomposition, particularly at high O/C ratios. Finally, this kinetic scheme provides a framework for understanding the oxidation chemistry of a model organic aerosol and informs parameterizations of more complex systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A granular flow model for dense planetary rings

In the present study of the viscosity of a differentially rotating particle disk, in the limiting case where the particles are densely packed and their collective behavior resembles that of a liquid, the pressure tensor is derived from both the equations of hydrodynamics and a simple kinetic model of collisions due to Haff (1983). Density waves and narrow circular rings are unstable if the liquid approximation applies, and the consequent nonlinear perturbations may generate 'splashing' of the ring material in the vertical direction. These results are pertinent to the origin of the ellipticities of ringlets, the nonaxisymmetric features near the outer edge of the Saturn B ring, and unexplained residuals in kinematic models of the Saturn and Uranus rings.

Borderies, N.↗

Modeling Solar-Wind Heavy-Ions' Potential Sputtering of Lunar KREEP Surface

Recent laboratory data suggest that potential sputtering may be an important weathering mechanism that can affect the composition of both the lunar surface and its tenuous exosphere; its role and implications, however, remain unclear. Using a relatively simple kinetic model, we will demonstrate that solar-wind heavy ions induced sputtering of KREEP surfaces is critical in establishing the timescale of the overall solar-wind sputtering process of the lunar surface. We will also also show that potential sputtering leads to a more pronounced and significant differentiation between depleted and enriched surface elements. We briefly discuss the impacts of enhanced sputtering on the composition of the regolith and the exosphere, as well as of solar-wind sputtering as a source of hydrogen and water on the moon.

Barghouty, A. F.↗

Numerical estimation of phase transformations in solid state during Yb:YAG laser heating of steel sheets

This work concerns the numerical modeling of heat transfer and phase transformations in solid state occurring during the Yb:YAG laser beam heating process. The temperature field is obtained by the numerical solution into transient heat transfer equation with convective term. The laser beam heat source model is developed using the Kriging interpolation method with experimental measurements of Yb:YAG laser beam profile taken into account. Phase transformations are calculated on the basis of Johnson - Mehl - Avrami (JMA) and Koistinen - Marburger (KM) kinetics models as well as continuous heating transformation (CHT) and continuous cooling transformation (CCT) diagrams for S355 steel. On the basis of developed numerical algorithms 3D computer simulations are performed in order to predict temperature history and phase transformations in Yb:YAG laser heating process.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

A Numerical Model for the Dynamics of Pickup Ions Outside the Heliopause and IBEX “Ribbon” Observation

The leading mechanism for the origin of the Interstellar Boundary Explorer (IBEX) “ribbon” of enhanced energetic neutral atoms (ENAs) from the outer heliosphere is the so-called secondary ENA process. In this mechanism, neutralized ions escape the heliosphere, and after several times of charge exchange collision,s some propagate back toward Earth in directions nearly perpendicular to the local interstellar magnetic field. However, the physical processes governing the distribution of the pickup ions (PUIs) outside the heliopause are still unclear. In this study, we build a new global model of the IBEX ribbon, where a key component is to calculate the dynamics of the PUIs outside the heliopause by solving the gyrophase-averaged focused transport equations on top of the background magnetohydrodynamic-kinetic model for the heliosphere-interstellar medium interaction. We discuss how the properties of the simulated ribbon change with different scattering parameters and show simulation results of the ENA sky maps, ribbon centers, ribbon profiles at several azimuthal slices, the ribbon ENA source region, and the velocity distribution of PUIs in the ENA source region for different scattering frequencies. Our results show that the model can provide reasonable ribbon sky maps comparable to the observed ribbon when the PUI mean free path is on the order of several thousand astronomical units.

79 ASTRONOMY AND ASTROPHYSICS↗

On analysis and stochastic modeling of the particle kinetic energy equation in particle-laden isotropic turbulent flows

In this study, we analyze three-dimensional particle-laden, isotropic turbulence to develop an understanding of inertial particle dynamics from a kinetic energy perspective. Data trends implying inhomogeneous sampling of the flow by particles are identified and used to support a proposed particle behavior: particles appear to accumulate in regions of low flow kinetic energy over time because they lose kinetic energy and slow down in such regions, ultimately causing them to spend more time there. To elucidate this behavior, we derive a particle kinetic energy equation from the particle momentum equation, which incorporates inertial effects through the Schiller–Naumann drag correlation. Upon extracting fundamental physics from this equation, hypotheses regarding the role of the Stokes number in the temporal change of particle kinetic energy and the previously proposed particle behavior are evaluated using simulation data considering three Stokes numbers. Finally, a Fokker–Planck equation is used to derive the steady-state probability density function of the particle kinetic energy. The model fits the simulation data well and provides a tool for further investigation into understanding preferential concentration, as well as a reduced order model for predicting particle kinetic energy in turbulent flows.

42 ENGINEERING↗

Consistent Chemical Mechanism from Collaborative Data Processing

Numerical tool of Process Informatics Model (PrIMe) is mathematically rigorous and numerically efficient approach for analysis and optimization of chemical systems. It handles heterogeneous data and is scalable to a large number of parameters. The Boundto-Bound Data Collaboration module of the automated data-centric infrastructure of PrIMe was used for the systematic uncertainty and data consistency analyses of the H 2 /CO reaction model (73/17) and 94 experimental targets (ignition delay times). The empirical rule for evaluation of the shock tube experimental data is proposed. The initial results demonstrate clear benefits of the PrIMe methods for an evaluation of the kinetic data quality and data consistency and for developing predictive kinetic models.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Towards Co-evolution of Membrane Transport and Metabolism

Protocellular boundaries were inextricably connected to the metabolism they encapsulated: to be inheritable, early metabolism must have led to an increased rate of growth and division of vesicles and, similarly, transport through vesicle boundaries must have supported the evolution of metabolism. Even though explaining how this coupling emerged and evolved in the absence of the complex machinery of modern cells is one of the key issues in studies on the origin of life, little is known about the biochemical and biophysical processes that might have been involved. This gap in our knowledge is a major impediment in efforts to construct scenarios for the origin of life and laboratory models of protocells. A combination of experimental and computational studies carried out by us and our collaborators is aimed at helping to close this gap. Properties of membranes might have contributed to the selection of RNA as an early biopolymer. A kinetic mechanism was proposed (Sacerdote & Szostak, 2005) in which ribose was supplied more quickly than other aldopentoses to primordial cells for preferential incorporation of ribonucleotides into nucleic acids. This proposal is based on a finding that ribose permeates membranes an order of magnitude faster than its diastereomers, arabinose and xylose. Our computer simulations, which yield permeation rates in excellent agreement with experiment, and kinetic modeling explain this phenomenon in terms of inter- and intramolecular interactions involving exocyclic hydroxyl groups attached to carbon atoms of the pyranose ring (Wei and Pohorille, 2009). They also constrain scenarios for the formation of the earliest nucleic acids (Wei and Pohorille, 2013). In one scenario, sugars permeate protocellular walls and subsequently are used to synthesize nucleic acids inside protocells. As long as this process proceeds at the rate faster than 6x10(exp -3)/s, ribose derivatives will be available for synthesis easier than their diastereomers. If nucleosides or their activated derivatives are synthesized outside protocells and subsequently transported across protocellular membranes the kinetic mechanism does not apply because all diastereomers, which have their sugars in the furanose rather than pyranose form, permeate the membrane at approximately the same rate. Properties of membranes might have been also coupled to metabolism involving peptides. Recently, Adamala and Szostak (2013) have shown that a dipeptide inside fatty-acid vesicles catalyzes the formation of another dipeptide that binds to vesicle walls and, by doing so, promotes their growth at the expense of other vesicles. This coupling of metabolism, permeability of vesicles and their growth is the first demonstration of evolutionary advantage imparted by small, membrane-bound peptides. Building on this work we have calculated the rate at which different blocked amino acids are delivered to a protocell for synthesis of dipeptides. We have further shown that the dipeptides are located at the water-membrane interface rather than in the center of the bilayer. On these basis it is anticipated that other dipeptides containing aromatic, but not necessarily hydrophobic amino acids (e.g. tyrosine) could have the same catalytic effects. Insight from these studies allows for estimating the rate of vesicle growth and the rates of dipeptide synthesis required to keep the system in balance. These results, in combination with our earlier studies, lead to a general scenario for evolution from membrane-bound dipeptides to ion channels in the origin of life.

Metabolism↗

Analysis and prediction of reaction kinetics using the degree of rate control

“Degree of rate control” (DRC) analysis provides a quantitative approach for analysing the kinetics of multi-step reaction mechanisms that has been widely applied to both heterogeneous and homogeneous catalysis research, as well as electrocatalysis. The DRC of any given transition state or intermediate is defined as a partial derivative such that it approximately equals the fractional increase in net rate to the product of interest per differential decrease in its standard-state free energy for that species (÷RT), holding constant the standard-state free energies of all other transition states and intermediates. Even very complex mechanisms usually have only a few species with non-zero DRCs and are thus the species whose interactions with the catalyst most strongly affect the net rate. These key DRC values thus offer a simple and intuitive route to optimize catalyst materials, especially with the assistance of computational methods like density functional theory (DFT). These high-DRC species are also the species whose energetics must be most accurately measured or calculated to achieve an accurate kinetic model for any reaction mechanism. In simple cases with a single “rate-determining step”, the DRC for its transition state (TS) is + 1. Catalyst-bound intermediates, on the other hand, often have negative DRCs equal to a small integer times their fractional occupancy of catalyst sites. The apparent activation energy equals a weighted average of the standard-state enthalpies (relative to reactants) of all the species (intermediates, transition states and products) in the reaction mechanism, each weighted by its DRC (+RT). It has been shown that the apparent transfer coefficient in electrocatalysis, an inverted form of the Tafel slope, is a weighted average of the number of electrons transferred to generate each intermediate or product species in the mechanism, weighted again by the DRC. Quantitative analysis of kinetic isotope effects (KIEs, or the ratio of net rates for different reactant isotopes) in complex mechanisms has shown that the logarithm of the KIE equals the weighted average over all species in the mechanism of the difference between the two isotopes in their standard-state free energies (÷RT), again weighted by the DRC. The reaction orders with respect to fluid-phase concentrations of reactants, products and intermediates have also been proven to be directly related to DRCs. Thus, there are numerous experimental observables which equate to short linear combinations of DRCs, so that combinations of experimental measurements might provide access to DRC values. Since its invention for transition states in 1994 and its generalization to include intermediates in 2009, the DRC has thus far mainly been calculated for microkinetic models based either entirely on DFT or on DFT with the key energies (i.e., those for high-DRC species) being fine-tuned to match experiments. The relationships summarized in this work provide new opportunities for using experiments earlier in the development and optimization of microkinetic models that require input from computational catalysis.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Experimental Investigation of Uranium Volatility during Vapor Condensation

The predictive models that describe the fate and transport of radioactive materials in the atmosphere following a nuclear incident (explosion or reactor accident) assume that uranium-bearing particulates would attain chemical equilibrium during vapor condensation. Here, we show that kinetically driven processes in a system of rapidly decreasing temperature can result in substantial deviations from chemical equilibrium. This can cause uranium to condense out in oxidation states (e.g., UO 3 vs UO 2 ) that have different vapor pressures, significantly affecting uranium transport. To demonstrate this, we synthesized uranium oxide nanoparticles using a flow reactor under controlled conditions of temperature, pressure, and oxygen concentration. The atomized chemical reactants passing through an inductively coupled plasma cool from ~5000 to 1000 K within milliseconds and form nanoparticles inside a flow reactor. The ex situ analysis of particulates by transmission electron microscopy revealed 2–10 nm crystallites of fcc-UO 2 or α-UO 3 depending on the amount of oxygen in the system. α-UO 3 is the least thermodynamically preferred polymorph of UO 3 . The absence of stable uranium oxides with intermediate stoichiometries (e.g., U 3 O 8 ) and sensitivity of the uranium oxidation states to local redox conditions highlight the importance of in situ measurements at high temperatures. Therefore, we developed a laser-based diagnostic to detect uranium oxide particles as they are formed inside the flow reactor. Our in situ measurements allowed us to quantify the changes in the number densities of the uranium oxide nanoparticles (e.g., UO 3 ) as a function of oxygen gas concentration. Our results indicate that uranium can prefer to be in metastable crystal forms (i.e., α-UO 3 ) that have higher vapor pressures than the refractory form (i.e., UO 2 ) depending on the oxygen abundance in the surrounding environment. This demonstrates that the equilibrium processes may not dominate during rapid condensation processes, and thus kinetic models are required to fully describe uranium transport subsequent to nuclear incidents.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Effect of fulvic acid surface coatings on plutonium sorption and desorption kinetics on goethite

The rates and extent of plutonium (Pu) sorption and desorption onto mineral surfaces are important parameters for predicting Pu mobility in subsurface environments. The presence of natural organic matter, such as fulvic acid (FA), may influence these parameters. We investigated the effects of FA on Pu(IV) sorption/desorption onto goethite in two scenarios: when FA was (1) initially present in solution or (2) found as organic coatings on the mineral surface. A low pH was used to maximize FA coatings on goethite. Experiments were combined with kinetic modeling and speciation calculations to interpret variations in Pu sorption rates in the presence of FA. Our results indicate that FA can change the rates and extent of Pu sorption onto goethite at pH 4. Differences in the kinetics of Pu sorption were observed as a function of the concentration and initial form of FA. The fraction of desorbed Pu decreased in the presence of FA, indicating that organic matter can stabilize sorbed Pu on goethite. These results suggest that ternary Pu–FA–mineral complexes could enhance colloid-facilitated Pu transport. In conclusion, more representative natural conditions need to be investigated to quantify the relevance of these findings.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗

Collisionless electrostatic interchange instabilities

The linear Vlasov dispersion equation for electrostatic plasma instabilities driven by gravity and weak density gradients perpendicular to a uniform magnetic field is derived and solved numerically. Two interchange instabilities emerge: the well-known fluid mode at long wavelengths and a kinetic model at wavelengths short compared with the ion gyroradius. The properties of both instabilities are studied, as well as the effects of gravity on the universal and lower-hybrid density drift instabilities. The results show that the kinetic interchange generally has a larger growth rate than the fluid interchange instability, indicating that, whenever the latter is present in a collisionless plasma, the former may also be found.

Gray, S. P.↗

CFD Analysis of Emissions for a Candidate N+3 Combustor

An effort was undertaken to analyze the performance of a model Lean-Direct Injection (LDI) combustor designed to meet emissions and performance goals for NASA's N+3 program. Computational predictions of Emissions Index (EINOx) and combustor exit temperature were obtained for operation at typical power conditions expected of a small-core, high pressure-ratio (greater than 50), high T3 inlet temperature (greater than 950K) N+3 combustor. Reacting-flow computations were performed with the National Combustion Code (NCC) for a model N+3 LDI combustor, which consisted of a nine-element LDI flame-tube derived from a previous generation (N+2) thirteen-element LDI design. A consistent approach to mesh-optimization, spraymodeling and kinetics-modeling was used, in order to leverage the lessons learned from previous N+2 flame-tube analysis with the NCC. The NCC predictions for the current, non-optimized N+3 combustor operating indicated a 74% increase in NOx emissions as compared to that of the emissions-optimized, parent N+2 LDI combustor.

nitrogen oxides↗

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↗

Cloud based tool for analysis of chemical kinetic mechanisms

The size and complexity of chemical kinetic models used to simulate combustion of hydrocarbon fuels continue to increase with our improved understanding of the underlying physical mechanisms, along with the improved ability to utilize such models as the result of increased computational power and efficiency. As mechanisms grow beyond thousands of species and tens of thousands of reactions, it becomes increasingly difficult to manually check for errors and physical inconsistencies. In this work, we present several automated methods to check for such issues in the specification of chemical reaction models. First, we demonstrate how discontinuities in thermodynamic data can cause simulation difficulties manifested as long wall-clock times and convergence failures. To correct this type of problem, we describe an automated method for refitting thermodynamic parameters. We also outline several methods to check the timescales of reaction rate coefficients to ensure physical consistency. All the methods are made available through a web-based tool (https://combustiontools.llnl.gov) intended to aid mechanism developers and users to improve the accuracy and performance of fuel models used by the combustion community.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Self-consistent quantum kinetic theory of diatomic molecule formation

A quantum kinetic theory of molecule formation is presented which includes three-body recombination and radiative association for a thermodynamically closed system which may or may not exchange energy with its surrounding at a constant temperature. The theory uses a Sturmian representation of a two-body continuum to achieve a steady-state solution of a governing master equation which is self-consistent in the sense that detailed balance between all bound and unbound states is rigorously enforced. The role of quasibound states in catalyzing the molecule formation is analyzed in complete detail. The theory is used to make three predictions which differ from conventional kinetic models. These predictions suggest significant modifications may be needed to phenomenological rate constants which are currently in wide use. Implications for models of low and high density systems are discussed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗