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 649 records · Page 36

Quantitative kinetic rules for plastic strain-induced α - ω phase transformation in Zr under high pressure

Plastic strain-induced phase transformations (PTs) and chemical reactions under high pressure are broadly spread in modern technologies, friction and wear, geophysics, and astrogeology. However, because of very heterogeneous fields of plastic strain $E$ p and stress σ tensors and volume fraction c of phases in a sample compressed in a diamond anvil cell (DAC) and impossibility of measurements of σ and $E$ p , there are no strict kinetic equations for them. Here, we develop a kinetic model, finite element method (FEM) approach, and combined FEM-experimental approaches to determine all fields in strongly plastically predeformed Zr compressed in DAC, and specific kinetic equation for α-ω PT consistent with experimental data for the entire sample. Since all fields in the sample are very heterogeneous, data are obtained for numerous complex 7D paths in the space of 3 components of the plastic strain tensor and 4 components of the stress tensor. Kinetic equation depends on accumulated plastic strain (instead of time) and pressure and is independent of plastic strain and deviatoric stress tensors, i.e., it can be applied for various above processes. Our results initiate kinetic studies of strain-induced PTs and provide efforts toward more comprehensive understanding of material behavior in extreme conditions.

36 MATERIALS SCIENCE↗

CO 2 methanation reaction pathways over unpromoted and NaNO 3 -promoted Ru/Al 2 O 3 catalysts

Catalytic CO 2 sorbents, materials that adsorb and pre-concentrate CO 2 on the catalyst surface prior to subsequent conversion, are becoming important materials in CO 2 capture and utilization. In this work, a prototypical CO 2 methanation catalyst – Ru/Al 2 O 3 – and a related catalytic sorbent – NaNO 3 /Ru/Al 2 O 3 – are used for CO 2 methanation in flowing hydrogen in a fixed bed reactor at temperatures ranging from 220 to 280 °C. Activation energies for the NaNO 3 /Ru/Al 2 O 3 material are slightly higher than unpromoted Ru/Al 2 O 3 catalysts, and the reaction orders vary more significantly. In situ IR spectroscopy and steady-state isotopic kinetic analysis (SSITKA) using in situ IR/MS spectroscopy show that bicarbonate and linear carbonyl species are the likely reaction intermediates over unpromoted Ru/Al 2 O 3 , while bidentate carbonate, formate and linear carbonyl species are among likely reaction intermediates over NaNO 3 /Ru/Al 2 O 3 . Rate laws consistent with the obtained experimental data are proposed after kinetic modeling of multiple plausible reaction pathways. In conclusion, results suggest that the pathway over the NaNO 3 /Ru/Al 2 O 3 catalyst likely has an additional kinetically relevant irreversible step in the CO 2 methanation reaction pathway.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

High-pressure ignition delay time measurements of a four-component gasoline surrogate and its high-level blends with ethanol and methyl acetate

Ignition delay times of a gasoline surrogate (iso-octane/toluene/n-heptane/1-hexene at 55/25/15/5% by liquid volume, developed in McCormick et al. 2017) with high-level bioblendstock-gasoline surrogate blends (50% and 85% biofuel by liquid volume of each ethanol and methyl acetate) were collected behind reflected shock waves. Post-reflected-shock temperatures ranged from 968 to 1361 K at pressures of about 4, 10, 25, and 50 atm. Data were collected for real fuel–air mixtures at fuel lean and stoichiometric conditions (φ = 0.5, 1.0) with focus on the more-practical, fuel-lean conditions. Ignition delay times were measured from OH* chemiluminescence around 307 nm using an endwall diagnostic. The RON and MON of the surrogate were 90.3 and 84.7, respectively. Ethanol was chosen due to its wide use in flex-fuel vehicles and availability, while also providing an increased octane rating. Methyl acetate was chosen for its especially high octane rating to investigate the effect of an extreme case. To validate the gasoline surrogate, the data are compared to real gasoline (RD387) and several gasoline surrogate experiments from the literature. Using the wide range of pressures studied, the gasoline surrogate’s pressure dependence was quantified to account for test-to-test variations using regression analysis. Similarly, a global correlation for gasoline and its surrogates was developed using all available data from the literature. Two modern chemical kinetics models targeting gasoline and its surrogates are compared to the ignition delay time measurements. These new high-pressure, high-bioblendstock concentration tests provide required chemical kinetic data for optimizing fuel and engine design.

42 ENGINEERING↗

Pressure-dependent kinetics of cyclopentene consumption

Recent experiments have uncovered deficiencies in the ability of current kinetic models for cyclopentane combustion to predict species profiles as a function of O2 concentration. One likely source of these discrepancies flows from the limited availability of relevant kinetic data for key functionalized intermediates, requiring the use of distant analogies to estimate their rates. This work fills this gap for cyclopentene, the most prominent C5 intermediate in low-temperature cyclopentane combustion. Ab initio transition-state-theory based master equation methods are used to calculate pressure-dependent rate constants across four potential energy surfaces (C5H8 + OH, C5H8 + HO2, C5H7, and C5H7 + O2). Potential energy surfaces are characterized at the CCSD(T)-F12/cc-pVDZ-F12//revDSD-PBEP86-D3BJ/def2-TZVPP level of theory. Pressure-dependent rate constants and branching fractions are determined via 1D master equation calculations. The results are used to analyze the competition for radicals between cyclopentene and cyclopentane in the intermediate temperature region, and to identify the dominant pathways flowing from these initiation channels. On the OH initiation surface, we find good agreement with recently published experimental rate constants at high temperature, and we elucidate the complex temperature dependence below 800 K, where 𝜋𝜋 addition becomes prominent. On the HO2 initiation surface, we find an interesting alternative route to bicyclic ether formation via well-skipping over the HO2 adduct well. On the cyclopentenyl radical surface, the major product is cyclopentadiene plus H. On the cyclopentenyl + O2 surface, we find that resonance stabilization largely prevents the allylic radical from undergoing oxidation, whereas the alkylic radical has a somewhat elevated propensity for oxidation pathways relative to cyclopentyl.

ab initio calculations↗

Large thermal hysteresis in a single-phase NiTiNb shape memory alloy

Large thermal hysteresis ($T$ hys ) is favorable in shape memory alloys (SMAs) because it renders wide service and storage temperature range of the SMA-made parts. In this letter, we show that large intrinsic $T$ hys up to 100 K is achieved in a low Nb content (2 at.%), β-Nb-free and single-phase NiTiNb SMA fabricated by casting, forging, wire-drawing and annealing. Such a large $T$ hys is comparable to that of the widely used dual-phase NiTiNb SMAs in which high Nb content and β-Nb are required. By exploiting the grain-size-dependence of $T$ hys in our single-phase NiTiNb SMA and that in an equiatomic NiTi SMA, we attribute the large $T$ hys to the enhanced kinetic resistance of transformation by the small grains. In conclusion, the connection between the $T$ hys and the kinetic resistance is demonstrated by a modified dislocation-based kinetic model.

36 MATERIALS SCIENCE↗

Resolving Dynamics in the Ensemble: Finding Paths through Intermediate States and Disordered Protein Structures

Proteins have been found to inhabit a diverse set of three-dimensional structures. The dynamics that govern protein interconversion between structures happen over a wide range of time scales-picoseconds to seconds. Our understanding of protein functions and dynamics is largely reliant upon our ability to elucidate physically populated structures. From an experimental structural characterization perspective, we are often limited to measuring the ensemble-averaged structure both in the steady-state and time-resolved regimes. Generating kinetic models and understanding protein structure–function relationships require atomistic knowledge of the populated states in the ensemble. Here, we present ensemble refinement methodologies that integrate time-resolved experimental signals with molecular dynamics models. We first discuss integration of experimental structural restraints to molecular models in disordered protein systems that adhere to the principle of maximum entropy for creating a complete set of ensemble structures. We then propose strategies to find kinetic pathways between the refined structures, using time-resolved inputs to guide molecular dynamics trajectories and the use of inference to generate tailored stimuli to prepare a desired ensemble of protein states.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Stiff neural ordinary differential equations

Neural Ordinary Differential Equations (ODEs) are a promising approach to learn dynamical models from time-series data in science and engineering applications. This work aims at learning neural ODEs for stiff systems, which are usually raised from chemical kinetic modeling in chemical and biological systems. We first show the challenges of learning neural ODEs in the classical stiff ODE systems of Robertson’s problem and propose techniques to mitigate the challenges associated with scale separations in stiff systems. We then present successful demonstrations in stiff systems of Robertson’s problem and an air pollution problem. The demonstrations show that the usage of deep networks with rectified activations, proper scaling of the network outputs as well as loss functions, and stabilized gradient calculations are the key techniques enabling the learning of stiff neural ODEs. The success of learning stiff neural ODEs opens up possibilities of using neural ODEs in applications with widely varying time-scales, such as chemical dynamics in energy conversion, environmental engineering, and life sciences.

97 MATHEMATICS AND COMPUTING↗

Measurement techniques for analysis of fission fragment excited gases

Spectroscopic analysis of fission fragment excited He, Ar, Xe, N2, Ne, Ar-N2, and Ne-N2 have been conducted. Boltzmann plot analysis of He, Ar and Xe have indicated a nonequilibrium, recombining plasma, and population inversions have been found in these gases. The observed radiating species in helium have been adequately described by a simple kinetic model. A more extensive model for argon, nitrogen and Ar-N2 mixtures was developed which adequately describes the energy flow in the system and compares favorably with experimental measurements. The kinetic processes involved in these systems are discussed.

Schneider, R. T.↗

Investigation of the hydrochlorination of SiCl4

The hyrochlorination of SiC14 and m.g. silicon metal to produce SiHC13, was investigated. Reaction kinetic measurements were carried out to collect additional rate data at 525 C and 550 C. A theoretical study was carried out to provide a kinetic model and a rate equation for the hydrochlorination reaction. Results of this preliminary study show that the rate of formation of SiHC13 follows a pseudo first order kinetics. The rate constants were measured at three temperatures, 550 C, 500 C and 450 C, respectively. The activation energy was determined from the Arrhenius plot to give a value of 13.2 Kcal/mole. The design of a quartz reactor to measure reaction rates and equilibrium conversion of SiHC13 at reaction temperature up to 650 C was completed.

Mui, J. Y. P.↗

Departures of the electron energy distribution from a Maxwellian in hydrogen. I - Formulation and solution of the electron kinetic equation. II - Consequences

The problem of calculating the steady-state free-electron energy distribution in a hydrogen gas is considered in order to study departures of that distribution from a Maxwellian at sufficiently low degrees of ionization. A model kinetic equation is formulated and solved analytically for the one-particle electron distribution function in a steady-state partially ionized hydrogen gas, and it is shown that the formal solution can be accurately approximated by using the WKB method. The solutions obtained indicate that the high-energy tail of the distribution is susceptible to distortion by imbalanced inelastic collisions for ionization fractions not exceeding about 0.1 and that such departures from a Maxwellian can lead to significant changes in the collisional excitation and ionization rates of ground-state hydrogen atoms. Expressions for the electron-hydrogen collision rates are derived which explicitly display their dependence on the hydrogen departure coefficients. The results are applied in order to compare self-consistent predictions with those based on the a priori assumption of a Maxwellian distribution for models of the thermal ionization equilibrium of hydrogen in the optically thin limit, spectral-line formation by a gas consisting of two-level atoms, and radiative transfer in finite slabs by a gas of four-level hydrogen atoms.

Shoub, E. C.↗

Robe Development for Electrical Conductivity Analysis in an Electron Gun Produced Helium Plasma

The use of magnetohydrodynamic (MHD) power conversion systems, potentially coupled with a fission power source, is currently being investigated as a driver for an advanced propulsion system, such as a plasma thruster. The efficiency of a MHD generator is strongly dependent on the electrical conductivity of the fluid that passes through the generator; power density increases as fluid conductivity increases. Although traditional MHD flows depend on thermal ionization to enhance the electrical conductivity, ionization due to nuclear interactions may achieve a comparable or improved conductivity enhancement while avoiding many of the limitations inherent to thermal ionization. Calculations suggest that nuclear-enhanced electrical conductivity increases as the neutron flux increases; conductivity of pure He-3 greater than 10 mho/m may be achievable if exposed to a flux greater than 10(exp 12) neutrons/cm2/s.) However, this remains to be demonstrated experimentally. An experimental facility has been constructed at the Propulsion Research Center at the NASA Marshall Space Flight Center, using helium as the test fluid. High energy electrons will be used to simulate the effects of neutron-induced ionization of helium gas to produce a plasma. These experiments will be focused on diagnosis of the plasma in a virtually static system; results will be applied to future tests with a MHD system. Initial experiments will utilize a 50 keV electron gun that can operate at up to a current of 200 micro A. Spreading the electron beam over a four inch diameter window results in an electron flux of 1.5x 10(exp 13) e/sq cm/s. The equivalent neutron flux that would produce the same ionization fraction in helium is 1x10(exp 12) n/sq cm/s. Experiments will simulate the neutron generated plasma modeled by Bitteker, which takes into account the products of thermal neutron absorption in He-3, and includes various ion species in estimating the conductivity of the resulting plasma. Several different probes will be designed and implemented to verify the plasma kinetics model. System parameters and estimated operating ranges are summarized. The predicted ionization fraction, electron density, and conductivity levels are provided in for an equivalent neutron flux of 1x10(exp 12) n/cm2/s. Understanding the complex plasma kinetics throughout a MHD channel is necessary to design an optimal power conversion system for space propulsion applications. The proposed experiments seek to fully characterize the helium plasma and to determine the reliability of each measurement technique, such that they may be applied to more advanced MHD studies. The expected value of each plasma parameter determined from theoretical models will be verified experimentally by several independent techniques to determine the most reliable method of obtaining each parameter. The results of these experiments will be presented in the final paper.

Bragg-Sitton, Shannon M.↗

Wire directed energy deposition of steel-aluminum structures using cold metal transfer process

In this study, a sharp transition from 316L stainless steel to 4043 aluminum alloy was fabricated using wire directed energy deposition (DED) via the cold metal transfer (CMT) process. The CMT process with its inherently low heat input, led to a significant reduction in intermetallic thickness at the bi-metallic interface compared to blown powder DED technique reported in the literature resulting in superior properties when compared to those of dissimilar steel-aluminum welds. Thermo-kinetic modeling confirmed that the intermetallic formation is through a classical nucleation and growth mechanism, and the fraction and thickness can be controlled by adjusting CMT process parameters to kinetically arrest or minimize the intermetallic formation. These findings underscore the efficacy of CMT-based wire DED for fabrication of steel-aluminum bi-metallic structures.

36 MATERIALS SCIENCE↗

Effect of the abolition of intersubunit salt bridges on allosteric protein structural dynamics

A salt bridge, one of the representative structural factors established by non-covalent interactions, plays a crucial role in stabilizing the structure and regulating the protein function, but its role in dynamic processes has been elusive. Here, to scrutinize the structural and functional roles of the salt bridge in the process of performing the protein function, we investigated the effects of salt bridges on the allosteric structural transition of homodimeric hemoglobin (HbI) by applying time-resolved X-ray solution scattering (TRXSS) to the K30D mutant, in which the interfacial salt bridges of the wild type (WT) are abolished. The TRXSS data of K30D are consistent with the kinetic model that requires one monomer intermediate in addition to three structurally distinct dimer intermediates (I 1 , I 2 , and I 3 ) observed in WT and other mutants. The kinetic and structural analyses show that K30D has an accelerated biphasic transition from I 2 to I 3 by more than nine times compared to WT and lacks significant structural changes in the transition from R-like I 2 to T-like I 3 observed in WT, unveiling that the loss of the salt bridges interrupts the R–T allosteric transition of HbI. Besides, the correlation between the bimolecular CO recombination rates in K30D, WT, and other mutants reveals that the bimolecular CO recombination is abnormally decelerated in K30D, indicating that the salt bridges also affect the cooperative ligand binding in HbI. These comparisons of the structural dynamics and kinetics of K30D and WT show that the interfacial salt bridges not only assist the physical connection of two subunits but also play a critical role in the global structural signal transduction of one subunit to the other subunit via a series of well-organized structural transitions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Temperature and density dependence of Kr L-shell spectrum in hot dense plasmas

Kr L-shell spectroscopy modeling results are discussed in this paper, focusing on the n = 4 to n = 2 line transitions of Be- and Li-like Kr ions. Collisional radiative atomic kinetic and Stark-broadened spectral line shape calculations show electron temperature Te and density ne sensitivity in the spectrum. The combination of the Te dependence due to the relative intensity of Be-like to Li-like line emissions in the range from 1.5 to 3 keV and the ne sensitivity from the Stark broadening effect on the line shapes in the range from 5×1023 to 2×1024/ cc results in a spectrum that can be employed to diagnose Te and ne. Two different collisional radiative atomic kinetic models i.e., Prismspect [J. J. MacFarlane, et al., Int. Fusion Sci. Appl. Conf. Proc. 457 (2003)] and ABAKO [Florido, et al., PRE, 80, 056402 (2009)] produce similar results in level populations and spectra. In x-ray spectroscopy of implosion cores, this Kr L-shell spectrum may prove useful in an intermediate Te range in which Ar is too ionized for its K-shell to be of diagnostic value and Kr is not ionized enough for its K-shell emission to be useful.

Physics↗

A Structural and Molecular Approach for the Study Biomarkers

Investigation of the nucleation and growth of crystals in both abiotic and biotic systems is critical to seemingly diverse disciplines of geology, biology, environmental science, and astrobiology. While there are abundant studies devoted to the determination of the structure and composition of inorganic crystals, as well as to the development of thermodynamic and kinetic models, it is only recently that research efforts have been directed towards understanding mineralization in biological systems (i.e., biomineralization). Biomineralization refers to the processes by which living organisms form inorganic solids. Studies of the processes of biomineralization under low temperature aqueous conditions have focused primarily on magnetite forming bacteria and shell forming marine organisms. Many of the biological building materials consist of inorganic minerals (calcium carbonate, calcium phosphate, silica or iron oxide) intricately combined with organic polymers (like proteins). More recently, efforts have been undertaken to explore the nature of biological activities in ancient rocks. In the absence of well-preserved microorganisms or genetic material required for the polmerase chain reaction (PCR) method in molecular phylogenetic studies, using biominerals as biomarkers offers an alternative approach for the recognition of biogenic activity in both terrestrial and extraterrestrial environments. The primary driving force in biomineralization is the interaction between organic and inorganic phases. Thus, the investigation of the ultrastructure and the nature of reactions at the molecular level occurring at the interface between inorganic and organic phases is essential to understanding the processes leading to the nucleation and growth of crystals. It is recognized that crystal surfaces can serve as the substrate for the organization of organic molecules that lead to the formation of polymers and other complex organic molecules, and in discussions of the origins of life, is referred to as organic synthesis on mineral surfaces. Furthermore, it is suggested that the interaction between mineral surfaces and simple organic molecules resulted in the formation of amino acids, RNA, and perhaps other more complex molecules such as proteins. On the other hand, in natural systems, it is recognized that functional groups on cell walls or membranes of microorganisms serve as sites of nucleation and crystallization. The precise replication of biominerals with controlled structure, morphology, size and texture is not confined to higher organisms as it also occurs in primitive prokaryotic cells such as magnetotactic bacteria and cyanobacteria. This suggests that the principal strategies of biomineralization were established early on in the evolutionary history of organisms. It is critical, therefore, to search for common mechanisms within diverse biological systems. One such common factor is the capability for organization and self-assembly. Organic macromolecules such as proteins and lipids can aggregate and polymerize forming membranes or extracellular matrix. At the organic-inorganic interface, several factors such as lattice geometry, polarity, stereochemistry and topography may act in concert to control nucleation and growth of crystals. Although several models have been proposed that discuss the significance of these factors for biomineralization, no comprehensive experimental data are available. In contrast to crystallization in exclusively inorganic systems, the kinetics of reaction and structural relationships between organic and inorganic phases in biominerals or biomimetic material is poorly understood. For example, it is not clear if the concept of epitactic growth (geometrical matching of unit cells at the interface of a secondary crystal growing on a primary crystal) applies to organic-inorganic systems. In contrast to inorganic templates that often have a smooth and rigid surface that promotes epitactic growth, biological substrates are usually rough and result in a large degree of mismatch. It is apparent that factors controlling the reaction at the crystal-matrix interface are strongly dependent upon the nature of the substrate. Therefore, characterization of the assembled organic surface and surface structure of the inorganic phase is crucial to understanding the processes of biomineralization. The focus of our research is the investigation of the processes leading to the nucleation and growth of crystals on both natural and synthetic systems through an interdisciplinary approach that integrates molecular biology, morphology and mineralogy using advanced preparation and analytical techniques. We have studied run-products, particularly magnetite, siderite and other carbonates, that resulted from extracellular biomineralization by extremophiles isolated from a variety of extreme environments ranging from permafrost to hydrothermal vent systems. The results of this study are critical to recognizing biomarkers in terrestrial and extraterrestrial environments.

Thomas-Keprta, Kathie↗

Flow diagnostics, chemically reacting turbulent flow modeling and turbulent mixing: A summary of project studies in combustion

This report includes copies of relevant publications, including conference proceedings, resulting from studies funded by NAG3-1010. Topics include flow visualization, LDV (laser doppler velocimeter), straight and curved simulations, and kinetics modeling of nitrogen oxides, including a turbulent mixing model for O2 and a maximum mixedness model.

Ghorashi, Bahman↗

Reversible Reactions, Mesh Size, and Segmental Dynamics Control Penetrant Diffusion in Ethylene Vitrimers

The diffusion of two aromatic dyes with nearly identical sizes was measured in ethylene vitrimers with precise linker lengths and borate ester cross-links using fluorescence recovery after photobleaching (FRAP). One dye possessed a reactive hydroxyl group, while the second was inert. The reaction of the hydroxyl group with the network is slow relative to the hopping times of the dye, resulting in a large slowdown by a factor of 50 for a reactive probe molecule. A kinetic model was fit to the fluorescence intensity data to determine rate constants for the reversible reaction of the dye from the network, which confirms the role of slow reaction kinetics. A second network cross-linker was also investigated with a substituted boronic ester showing ∼10,000 times faster exchange kinetics. In this system, the two dyes show the same diffusion coefficient, as the reaction is no longer the rate-limiting step. The role of dense meshes on small and large dyes is also discussed in the context of the existing theories. Finally, these results highlight the potential of dynamic networks to control penetrant transport through synergistic effects of the mesh size, dynamic bond kinetics, and penetrant–network interactions.

confinement↗

ECUT (Energy Conversion and Utilization Technologies) program: Biocatalysis project

The Annual Report presents the fiscal year (FY) 1988 research activities and accomplishments, for the Biocatalysis Project of the U.S. Department of Energy, Energy Conversion and Utilization Technologies (ECUT) Division. The ECUT Biocatalysis Project is managed by the Jet Propulsion Laboratory, California Institute of Technology. The Biocatalysis Project is a mission-oriented, applied research and exploratory development activity directed toward resolution of the major generic technical barriers that impede the development of biologically catalyzed commercial chemical production. The approach toward achieving project objectives involves an integrated participation of universities, industrial companies and government research laboratories. The Project's technical activities were organized into three work elements: (1) The Molecular Modeling and Applied Genetics work element includes research on modeling of biological systems, developing rigorous methods for the prediction of three-dimensional (tertiary) protein structure from the amino acid sequence (primary structure) for designing new biocatalysis, defining kinetic models of biocatalyst reactivity, and developing genetically engineered solutions to the generic technical barriers that preclude widespread application of biocatalysis. (2) The Bioprocess Engineering work element supports efforts in novel bioreactor concepts that are likely to lead to substantially higher levels of reactor productivity, product yields and lower separation energetics. Results of work within this work element will be used to establish the technical feasibility of critical bioprocess monitoring and control subsystems. (3) The Bioprocess Design and Assessment work element attempts to develop procedures (via user-friendly computer software) for assessing the energy-economics of biocatalyzed chemical production processes, and initiation of technology transfer for advanced bioprocesses.

Baresi, Larry↗