Search NASA⌕ Search

SEARCH · Search NASA

Results for “Kinetic Modelling”

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 775 records · Page 43

Fundamental Studies of Bifunctional Catalysts for Tandem Reactions

During the past few years of our funding from DOE Basic Energy Sciences, we have focused on the catalytic conversion of biomass-derived 5-hydroxymethyl furfural (HMF) as a platform molecule for sustainable synthesis of value-added chemicals. The combination of the diene and aldehyde functionalities in HMF enables catalytic production of acetalized HMF derivatives with diol or epoxy reactants to allow reversible synthesis of norcantharimide derivatives upon Diels-Alder reaction with maleimides. Whereas the electron-withdrawing nature of aldehyde group inhibits the reactivity of diene group in HMF to undergo Diels-Alder reactions with dienophiles, by converting the aldehyde group in HMF to an electron-donating hydroxyl group, HMF can append maleimide-based chemicals through Diels-Alder coupling. Acetalization of HMF not only produces an electron-rich diene for Diels-Alder reaction, but it also alters the reversibility of acetal formation and thereby allows for the controlled release of chemicals by hydrolysis. Therefore, we synthesized various norcantharimide derivatives from biomass-derived HMF by acetalization over an acid catalyst, Amberlyst-15, followed by Diels-Alder reaction with maleimide. The norcantharimides release the starting materials by retro Diels-Alder reaction that is triggered by acetal hydrolysis under acidic (≤pH 3) conditions. In subsequent work, we explored the synthesis and properties of functional polyurethanes and polyesters with tunable properties from biomass-derived (HMF)-Acetone-HMF (HAH) monomers. HAH can be selectively hydrogenated over Cu and Ru catalysts to produce partially-hydrogenated HAH (PHAH) and fully-hydrogenated HAH (FHAH). HAH contains functional groups, including hydroxyl, furan, enone, and ketone functionalities, that can be exploited to further tune the polymer properties. The π-electron conjugation between the enone and furan groups in HAH monomer stabilizes the furan group to prevent the uncontrolled degradation of furans during reactions, such as etherification. Therefore, HAH is a renewable monomer that can be used to synthesize polymers, having high molecular weight, symmetric functionalities, an inexpensive production price, and being derived from renewable resources without competition with food resources. We demonstrated the synthesis of functional polyurethanes and a renewable polyester from HAH-derived monomers, and 4,4’-methylenebis(phenyl isocyanate) (MDI). We showed that it is essential to be able to selective hydrogenate C=C double bonds in HAH to be able to achieve effective Diels Alder reactions of the furan ring. Therefore, we studied the hydrogenation at temperatures from 313 – 393 K of HAH over Pd, Ru, and Cu based catalysts. HAH was selectively hydrogenated to produce partially-hydrogenated monomers (PHAH) over Cu and Ru catalysts and to fully-hydrogenated HAH monomers (FHAH) over the Ru catalyst. Pd based catalysts yielded a mixture of partially and fully hydrogenated monomers. Reaction kinetics models were employed to quantify the kinetic behavior for hydrogenation over Ru, Cu, and Pd catalysts. A 5-step pathway exhibited over Pd and Ru catalysts consists of both series and parallel reaction steps, where HAH is both converted to fully hydrogenated products sequentially via series reactions of partially hydrogenated intermediates, as well as converted directly in parallel reactions to form the fully hydrogenated products. In contrast, a 3-step pathway over the Cu catalyst consists only of the consecutive reaction steps, where the final product was formed via series reactions of intermediate products. Additionally, reaction over the Cu catalyst did not hydrogenate the furan rings of the HAH molecule and yielded a different final product than those hydrogenation over Pd and Ru catalysts. Using the results from our detailed reaction kinetics studies, we developed universal conditions for the maximum production of the various hydrogenated products for both batch and plug flow systems. In subsequent work, we employed reaction schemes to optimize the yields of products from selective hydrogenations of HAH in isopropanol solvent and in the presence of liquid water. Reaction schemes consisting of 7, 9, and 11 steps were examined to describe the rates of formation of the observed products and reaction intermediates for hydrogenation of HAH over Ru and Pd catalysts, and a 3-step scheme was studied over Cu catalysts. Rate constants and activation energies were calculated using these reaction schemes, and we then apply these reaction schemes to explore the effects of water addition on the hydrogenation pathways. The effects of water addition to isopropanol (IPA) solvents on the hydrogenation of HAH were markedly different over Pd, Ru, and Cu catalysts. Over the Pd catalyst, the addition of water to IPA increased hydrogenation rates and promoted hydrogenation of furan rings. The addition of water to IPA yielded significant carbon losses over the Ru catalyst, and slowed hydrogenation steps over Cu, while significantly inhibiting hydrogenation of the ketone group. This behavior opened routes toward increased production rates of PHAH=O (a partially hydrogenated form of HAH containing a C=O bond), a product in which the diene groups of the furan rings were not hydrogenated. Addition of water also allowed increased feed concentrations of HAH that were previously not possible in pure IPA solvents.

09 BIOMASS FUELS↗

Lower length scale model for palladium attack of silicon carbide in TRISO fuel

TRistructural ISOtropic (TRISO) particle fuels rely on silicon carbide (SiC) as the primary barrier for metallic fission product (FP) release. Palladium (Pd) generated by fission degrades the SiC layer, resulting in the formation of lamellar layers of palladium silicides (PdxSi) and carbon (C) perpendicular to the direction of attack. The Pd attack has been hypothesized to be responsible for failure of the SiC layer and enhance FP release. To better understand and quantify Pd attack of SiC in TRISO particles, a multiscale, mechanistic model of Pd transport is being developed by the NEAMS program. Previous work provided an initial hypothesis for modeling lamellar microstructure formation in SiC due to Pd attack using a phase-field model. The work described in this report builds on the previous model by using molecular dynamics (MD) simulations to parameterize the phase-field model kinetics, and build a reduced order model in BISON using the improved mesoscale Pd penetration model.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

A Test of the Interstellar Boundary EXplorer Ribbon Formation in the Outer Heliosheath

NASA’s Interstellar Boundary EXplorer ( IBEX ) mission is imaging energetic neutral atoms (ENAs) propagating to Earth from the outer heliosphere and local interstellar medium (LISM). A dominant feature in all ENA maps is a ribbon of enhanced fluxes that was not predicted before IBEX . While more than a dozen models of the ribbon formation have been proposed, consensus has gathered around the so-called secondary ENA model. Two classes of secondary ENA models have been proposed; the first class assumes weak scattering of the energetic pickup protons in the LISM, and the second class assumes strong but spatially localized scattering. Here we present a numerical test of the “weak scattering” version of the secondary ENA model using our gyro-averaged kinetic model for the evolution of the phase-space distribution of protons in the outer heliosheath. As input for our test, we use distributions of the primary ENAs from our MHD-plasma/kinetic-neutral model of the heliosphere-LISM interaction. The magnetic field spectrum for the large-scale interstellar turbulence and an upper limit for the amplitude of small-scale local turbulence (SSLT) generated by protons are taken from observations by Voyager 1 in the LISM. The hybrid simulations of energetic protons are also used to set the bounding wavenumbers for the spectrum of SSLT. Our test supports the “weak scattering” version. This makes an additional solid step on the way to understanding the origin and formation of the IBEX ribbon and thus to improving our understanding of the interaction between the heliosphere and the LISM.

79 ASTRONOMY AND ASTROPHYSICS↗

A Two-Fluid Study of Oblique Tearing Modes and 3-D Magnetic Reconnection

The goals and highlights of this research are to: 1) Undertake two-fluid simulations of oblique tearing modes with $m_{i}/m_{e}$ = 100 and 1836, and $b_g$ ≤ 10 with further bridging the fluid and kinetic models and linear growth rates agree with kinetic theory for $b_g$ ≤ 1; 2) 3-D simulations of a doubly-periodic Harris Sheet with $m_{i}/m_{e}$ = 100 and $b_g$ = 2:5 that features current layer splits into multiple ~ de-scale diffusion regions, an amplification of current by 2x, and similar features to kinetic simulations; 3) Influence of boundary conditions: periodic vs. line-tied; and 4) Address how much of the kinetic nature of MR can be captured by a two-fluid model. (A fluid model requires dissipation.)

58 GEOSCIENCES↗

Development of Kinetics and Mathematical Models for High-Pressure Gasification of Lignite-Switchgrass Blends: Cooperative Research and Development Final Report, CRADA Number CRD-11-447

NREL will work with Participant as a subtier partner under DE-FOA-0000240 titled "Co-Production of Power, Fuels, and Chemicals via Coal/Biomass Mixtures." The goal of the project is to determine the gasification characteristics of switchgrass and lignite mixtures and develop kinetic models. NREL will utilize a pressurized thermogravimetric analyzer to measure the reactivity of chars generated in a pressurized entrained-flow reactor at Participant's facilities and to determine the evolution of gaseous species during pyrolysis of switchgrass-lignite mixtures. Mass spectrometry and Fourier-transform infrared analysis will be used to identify and quantify the gaseous species. The results of the project will aid in defining key reactive properties of mixed coal biomass fuels.

09 BIOMASS FUELS↗

Modeling of feature profile evolution for ion etching

A kinetic model is presented to investigate the profile evolution during ion etching. The effects of ion sputtering, redeposition, and diffusion processes are all taken into consideration in the formation mechanism of surface profile. The dominant factors accounting for the surface smoothening and roughening during ion etching are well explained in this study. Under high ion flux or ion energy, the sputtering effect plays a controlling role in roughening the surface profile with a high etching rate. While decreasing ion flux or ion energy, the surface profile is smoothened by the diffusion mechanism with a long time ion irradiation. For a low temperature, the characteristic length of nanostructures decreases with a sputtered feature profile due to the low mobility. Our simulation results are consistent well with many experimental observations. This theoretical model provides an efficient numerical approach to fully understand the mechanism for the formation of surface profile allowing for designing of appropriate experiments to form specific nanostructures through ion-beam technology.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Parameter study of a Kinetic Phase Transition model

We investigate the ability of the two model parameters in the Kinetic Phase Transition (KPT) model of Carl Greeff to describe phase transitions of various types when combined with the tin Equation of states for the phases, SESAME 2162. While the two parameters are correlated, one combination can be associated with the nucleation process giving rise to the hysteresis seen in phase transition paths in pressure-temperature space and another can be associated with the overall timescale of the transition. We acknowledge the desire to have a model where the parameters are more independent of each other for easier model parameter determination.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

A review of low-rank methods for time-dependent kinetic simulations

Time-dependent kinetic models are ubiquitous in computational science and engineering. The underlying integro-differential equations in these models are high-dimensional, comprised of a six–dimensional phase space, making simulations of such phenomena extremely expensive. In this article we demonstrate that in many situations, the solution to kinetics problems lives on a low dimensional manifold that can be described by a low-rank matrix or tensor approximation. We then review the recent development of so-called low-rank methods that evolve the solution on this manifold. The two classes of methods we review are the dynamical low-rank (DLR) method, which derives differential equations for the low-rank factors, and a Step-and-Truncate (SAT) approach, which projects the solution onto the low-rank representation after each time step. Thorough discussions of time integrators, tensor decompositions, and method properties such as structure preservation and computational efficiency are included. We further show examples of low-rank methods as applied to particle transport and plasma dynamics.

97 MATHEMATICS AND COMPUTING↗

Integrative mapping reveals molecular features underlying the mechanism of nucleocytoplasmic transport

Nuclear pore complexes (NPCs) enable rapid, selective, and robust nucleocytoplasmic transport. To explain how transport emerges from the system components and their interactions, we used experimental data and theoretical information to construct an integrative Brownian dynamics model of transport through an NPC, coupled to a kinetic model of transport in the cell. The model recapitulates key aspects of transport for a wide range of molecular cargoes, including preribosomes and viral capsids. Our model quantifies how flexible phenylalanine-glycine (FG) repeat proteins create an entropic barrier to passive diffusion and how this barrier is selectively lowered in facilitated diffusion by the many transient interactions of nuclear transport receptors with the FG repeats. Selective transport is enhanced by “fuzzy” multivalent interactions, redundant FG repeat mass, coupling to the energy-dependent RanGTP concentration gradient, and exponential dependence of transport kinetics on the transport barrier. Our model will facilitate rational modulation of the NPC and its artificial mimics.

Raveh, Barak [School of Computer Science and Engin↗

Direct implicit and explicit energy-conserving particle-in-cell methods for modeling of capacitively coupled plasma devices

Achieving large-scale kinetic modeling is a crucial task for the development and optimization of modern plasma devices. With the trend of decreasing pressure in applications, such as plasma etching, kinetic simulations are necessary to self-consistently capture the particle dynamics. The standard, explicit, electrostatic, momentum-conserving particle-in-cell method suffers from restrictive stability constraints on spatial cell size and temporal time step, requiring resolution of the electron Debye length and electron plasma period, respectively. This results in a very high computational cost, making the technique prohibitive for large volume device modeling. We investigate the direct implicit algorithm and the explicit energy conserving algorithm as alternatives to the standard approach, both of which can reduce computational cost with a minimal (or controllable) impact on results. These algorithms are implemented into the well-tested EDIPIC-2D and LTP-PIC codes, and their performance is evaluated via 2D capacitively coupled plasma discharge simulations. The investigation reveals that both approaches enable the utilization of cell sizes larger than the Debye length, resulting in a reduced runtime, while incurring only minor inaccuracies in plasma parameters. The direct implicit method also allows for time steps larger than the electron plasma period; however, care must be taken to avoid numerical heating or cooling. It is demonstrated that by appropriately adjusting the ratio of cell size to time step, it is possible to mitigate this effect to an acceptable level.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

JP-10 combustion studied with shock tube experiments and modeled with automatic reaction mechanism generation

Our research presents shock tube experiments and kinetic modeling efforts on the pyrolysis and combustion of JP-10. The experiments were performed at 6–8 atm using 2000 ppm of JP-10 over a temperature range of 1000–1600 K for pyrolysis and oxidation equivalence ratios from 0.14 to 1.0. This work distinguishes itself from previous studies as GC/MS was used to identify and quantify the products within the shocked samples, enabling the tracking of product yield dependence on equivalence ratio as well as identifying several new intermediates that form during JP-10’s decomposition. A detailed, comprehensive model of JP-10’s combustion and pyrolysis kinetics was constructed with the help of RMG, an open-source reaction mechanism generation software package. The resulting model, which includes 691 species reacting in 15,518 reactions, was extensively validated against the shock tube experimental dataset as well as newly published flow tube pyrolysis data from Ghent. Most of the significant rate coefficients were computed using quantum chemistry. The model succeeds in identifying all major pyrolysis and combustion products and captures key trends in the product distribution. Simulated ignition delays agree within a factor of 4 with most experimental ignition delay data gathered from literature. The proposed experimental work and modeling efforts yield new insights on JP-10’s complex decomposition and oxidation chemistry and identify key pathways towards aromatics formation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

X-ray heating and electron temperature of laboratory photoionized plasmas

We discuss experimental and modeling results for the x-ray heating and temperature of laboratory photoionized plasmas. A method is used to extract electron temperature based on the analysis of transmission spectroscopy data that is independent of atomic kinetics modeling. Boltzmann kinetic and radiation-hydrodynamics simulations emphasized the critical interplay between atomic physics and heating in photoionized plasmas. They also demonstrated the dramatic impact of photoexcitation on excited state populations, line emissivity, and radiation cooling. Temperature calculations performed with astrophysical codes significantly overestimated the measurement.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

First steps toward predicting corrosion behavior of structural materials in molten salts

To address the need for physics-based models to predict corrosion behavior of materials in molten salts, here we propose potential methods to aid in selection of optimum materials for structural components in molten salt powered technologies. In the present work, the role of alloy thermodynamics and kinetics on governing corrosion rates of Ni-based alloys will be discussed combining experimental and computational methods. A few strategies are presented to quantify corrosion rates of Ni-based materials isothermally exposed in purified KCl-MgCl 2 salts at 600° C-800° C. The influence of capsule material, potential corrosion products and role of alloy composition on the observed corrosion rates was discussed with coupled thermodynamic-kinetic models. Larger depths of corrosion attack were observed in alloy 230 specimens compared to other alloys under similar conditions was attributed to the much higher chemical activity of Cr in the alloy that results in a larger Cr chemical potential gradient.

36 MATERIALS SCIENCE↗

Fast microflow kinetics and acid catalyst deactivation in glucose conversion to 5-hydroxymethylfurfural

5-Hydroxymethyl furfural (HMF) is an important platform chemical because it can be upgraded to various drop-in and performance-advantaged products. The cascade reaction of HMF production from glucose over a Lewis acid (CrCl 3 ) and a Brønsted acid (HCl) catalyst in aqueous media is investigated in a microreactor at short residence times and high temperatures. We study the formation of various chromium species using UV-visible spectrophotometry and elucidate the Cr(III) speciation. The catalyst reactivity increases sharply at short residence times, and then drops at long times. This indicates that the catalyst treatment plays a vital role in getting optimal reactivity, and recording the catalyst history is necessary. We develop a kinetic model to describe the catalyst speciation as well as the Lewis and Brønsted acid-catalyzed reaction kinetics using a hierarchical approach. Here, the model is in good agreement with experiments. We demonstrate the benefits of tandem Lewis-external added Brønsted acid catalysis in processing time, productivity, and catalyst stability. We apply this model to optimize the HMF yield and obtain ~36% yield at 200 °C in 7 min and report the highest productivity of ~15% yield per min, demonstrating the opportunity of reaching high productivity at short residence times.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The Importance of Decarbonylation Mechanisms in the Atomic Layer Deposition of High-Quality Ru Films by Zero-Oxidation State Ru(DMBD)(CO) 3

Achieving facile nucleation of noble metal films through atomic layer deposition (ALD) is extremely challenging. To this end, η 4 -2,3-dimethylbutadiene ruthenium tricarbonyl (Ru(DMBD)(CO) 3 ), a zero-valent complex, has recently been reported to achieve good nucleation by ALD at relatively low temperatures and mild reaction conditions. Here, we study the growth mechanism of this precursor by in situ quartz-crystal microbalance and quadrupole mass spectrometry during Ru ALD, complemented by ex situ film characterization and kinetic modeling. These studies reveal that Ru(DMBD)(CO) 3 produces high-quality Ru films with excellent nucleation properties. This results in smooth, coalesced films even at low film thicknesses, all important traits for device applications. However, Ru deposition follows a kinetically limited decarbonylation reaction scheme, akin to typical CVD processes, with a strong dependence on both temperature and reaction timescale. The non-self-limiting nature of the kinetically driven mechanism presents both challenges for ALD implementation and opportunities for process tuning. By surveying reports of similar precursors, we suggest that the findings can be generalized to the broader class of zero-oxidation state carbonyl-based precursors used in thermal ALD, with insight into the design of effective saturation studies.

chemical vapor deposition↗

Shock Tube and Flame Speed Measurements of 2,4,4-Trimethyl-1-Pentene: A Co-Optima Biofuel

Abstract The combustion of 2,4,4-trimethyl-1-pentene (diisobutylene, C8H16), which is a biofuel and a component of surrogate fuels, is examined in this work. Carbon monoxide time–histories and ignition delay times are collected behind reflected shock waves utilizing a shock tube and mid-infrared laser absorption spectroscopy. Measurements were obtained near 10 atm pressure during stoichiometric oxidation of 0.15%C8H16/O2/Ar. Simulated results from chemical kinetic models are provided, and sensitivity analyses are used to discuss differences between models for both ignition delay times and carbon monoxide formation. In addition, laminar burning speeds are obtained at 1 atm, 428 K, and equivalence ratios, phi, between 0.91 and 1.52 inside a spherical chamber facility. Measured burning speeds are found to be less than that of ethanol over the equivalence ratio span. Burning speed measurements are compared to predictions of chemical kinetic mechanisms and are in agreement for the richest conditions; however, at lean conditions, the model predicts a far slower-burning speed. The maximum burning speed occurs at an equivalence ratio of 1.08 with a magnitude of 0.70 m/s. The current work provides the crucial experimental data needed for assessing the feasibility of this biofuel and for the development of future combustion chemical kinetics models.

Energy & Fuels↗

Development Of Thermodynamic and Kinetic Simulation Tools and Testing Procedures for Enhanced Durability of Concrete Containing Industrial By-Products

This project developed screening tools that enable evaluation of alternative cementitious binders that create concretes to significantly reduce energy and emissions while remaining cost competitive on both initial and long-term costs. The team began with the viewpoint that acceptance of new cementitious binder products has a substantially greater chance of successful implementation when capital investment is not excessive and the end product has customers that have experience and resources to use this. As such, a binder system that is based on portland cement with blended using industrial by-products (alternative cementitious materials) has the potential for dramatic and meaningful impact. The team has focused on developing implementable solutions in specifications and current practice. This however requires three main factors: 1) ability to screen byproducts and alternative materials for success, 2) ability to ‘treat’ materials chemically to enhance kinetics, and 3) ability to provide predictions of performance of both binders and concrete from first principles. The project developed/refined a state of the art and scientifically based screening test for SCM called the pozzolanic reactivity test. The team developed kinetic models to simulate these materials as well as experimental approaches to alter selected reactions. Simulation tools were developed that enable the performance of concrete to be predicted based on the chemistry and reactivity of the cement and alternative SCM. Specifically, this project: • enhanced the kinetic reactivity models for use in multi-scale computational programs that use thermodynamics to predict reaction products. 2 • developed scaling models to extend thermodynamic modeling to link these models with pore structure. This enables strength, transport property, and coupled transport prediction. • developed tools to predict performance of cementitious materials using the pozzolanic reactivity test and chemical composition. The predicted properties are consistent with AASHTO R101 and the CEB-FIP model code and can be measured using associated test procedures. • demonstrated mechanical and fracture based modeling tools that thermodynamic predictions and inputs to predict concrete service life. These results have been used to demonstrate the value of enabling specifications to include ASTM C 595 cement as well ASTM C150 cement. In addition, these products are being used to expedite the evaluation of alternative SCM to aide in determining which materials have potential value and what ‘compositions’ of blended cements merit further investment.

42 ENGINEERING↗

Pre-steady-state kinetics of nanocrystal:molybdenum nitrogenase biohybrids reveals hole-scavenging efficiency is critical to N 2 reduction

Molybdenum (Mo) nitrogenase is a two-component enzyme complex that catalyzes the reduction of dinitrogen to ammonia and protons to hydrogen gas. We have shown that electrons for dinitrogen reduction can be delivered photochemically to the catalytic MoFe protein component by cadmium sulfide (CdS) nanocrystals. In this study, we used electron paramagnetic resonance spectroscopy to measure the transient populations of catalytic intermediates. We fit the populations with a pre-steady-state kinetic model, which allowed us to distinguish between productive and non-productive reaction pathways and extract the rate constants for the reaction. Our results demonstrated that the rate of catalytic electron delivery into MoFe protein increased with the concentration of the sacrificial electron donor. This enabled electron delivery to exceed the rate of hydride protonation, a relaxation pathway that competes with N 2 binding. Thus, managing the balance between electron transfer and hole transfer reactions is required to achieve a kinetic regime that favors N 2 reduction.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗