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 1,387 records · Page 77
All sourcesNASA NTRSDOE OSTI
DOE OSTI2017

Stability of Electrodeposition at Solid-Solid Interfaces and Implications for Metal Anodes

We generalize the conditions for stable electrodeposition at isotropic solid-solid interfaces using a kinetic model which incorporates the effects of stresses and surface tension at the interface. We develop a stability diagram that shows two regimes of stability: a previously known pressure-driven mechanism and a new density-driven stability mechanism that is governed by the relative density of metal in the two phases. We show that inorganic solids and solid polymers generally do not lead to stable electrodeposition, and provide design guidelines for achieving stable electrodeposition.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2018

Nonequilibrium Pair Breaking in Ba ( Fe 1 - x Co x ) 2 As 2 Superconductors: Evidence for Formation of a Photoinduced Excitonic State

Ultrafast terahertz (THz) pump-probe spectroscopy reveals an unusual out-of-equilibrium Cooper pair nonlinear dynamics and a nonequilibrium state driven by femtosecond (fs) photoexcitation of superconductivity (SC) in iron pnictides. Following fast SC quench via hot-phonon scattering, a second, abnormally slow (many hundreds of picoseconds), SC quench regime is observed prior to any recovery. Importantly, a nonlinear pump fluence dependence is identified for this remarkably long prebottleneck dynamics that are sensitive to both doping and temperature. Here, using quantum kinetic modeling we argue that the buildup of excitonic interpocket correlation between electron-hole (e-h) quasiparticles (QP) quenches SC after fs photoexcitation leading to a long-lived, many-QP excitonic state.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Terahertz Second-Harmonic Generation from Lightwave Acceleration of Symmetry-Breaking Nonlinear Supercurrents

We report terahertz (THz) light-induced second harmonic generation, in superconductors with inversion symmetry that forbid even-order nonlinearities. The THz second harmonic emission vanishes above the superconductor critical temperature and arises from precession of twisted Anderson pseudospins at a multicycle, THz driving frequency that is not allowed by equilibrium symmetry. We explain the microscopic physics by a dynamical symmetry breaking principle at sub-THz-cycle by using quantum kinetic modeling of the interplay between strong THz-lightwave nonlinearity and pulse propagation. In this work, the resulting nonzero integrated pulse area inside the superconductor leads to light-induced nonlinear supercurrents due to subcycle Cooper pair acceleration, in contrast to dc-biased superconductors, which can be controlled by the band structure and THz driving field below the superconducting gap.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Phonon-Assisted Ballistic Current from First-Principles Calculations

The bulk photovoltaic effect (BPVE) refers to current generation due to illumination by light in a homogeneous bulk material lacking inversion symmetry. In addition to the intensively studied shift current, the ballistic current, which originates from asymmetric carrier generation due to scattering processes, also constitutes an important contribution to the overall kinetic model of the BPVE. In this Letter, we use a perturbative approach to derive a formula for the ballistic current resulting from the intrinsic electron-phonon scattering in a form amenable to first-principles calculation. We then implement the theory and calculate the ballistic current of the prototypical BPVE material BaTiO 3 using quantum-mechanical density functional theory. The magnitude of the ballistic current is comparable to that of the shift current, and the total spectrum (shift plus ballistic) agrees well with the experimentally measured photocurrents. Furthermore, we show that the ballistic current is sensitive to structural change, which could benefit future photovoltaic materials design.

14 SOLAR ENERGY↗
DOE OSTI2017

Elucidation of Iron Gettering Mechanisms in Boron-Implanted Silicon Solar Cells

To facilitate cost-effective manufacturing of boron-implanted silicon solar cells as an alternative to BBr 3 diffusion, we performed a quantitative test of the gettering induced by solar-typical boron-implants with the potential for low saturation current density emitters (< 50 fA/cm 2 ). We show that depending on the contamination level and the gettering anneal chosen, such boron-implanted emitters can induce more than a 99.9% reduction in bulk iron point defect concentration. The iron point defect results as well as synchrotron-based Nano-X-ray-fluorescence investigations of iron precipitates formed in the implanted layer imply that, with the chosen experimental parameters, iron precipitation is the dominant gettering mechanism, with segregation-based gettering playing a smaller role. We reproduce the measured iron point defect and precipitate distributions via kinetics modeling. First, we simulate the structural defect distribution created by the implantation process, and then we model these structural defects as heterogeneous precipitation sites for iron. Unlike previous theoretical work on gettering via boron- or phosphorus-implantation, our model is free of adjustable simulation parameters. The close agreement between the model and experimental results indicates that the model successfully captures the necessary physics to describe the iron gettering mechanisms operating in boron-implanted silicon. Furthermore, this modeling capability allows high-performance, cost-effective implanted silicon solar cells to be designed.

14 SOLAR ENERGY↗
DOE OSTI2024

Using in situ UO 2 bicrystal sintering to understand grain boundary dislocation nucleation kinetics and creep

Capillary evolution at bicrystal UO 2 grain boundaries is characterized using in situ transmission electron microscopy. The discontinuous nature of the densification process, both particle rotation and axial strain, along with the large activation stress for densification support a hypothesis that grain boundary strain in UO 2 follows nucleation rate limited kinetics at low to intermediate stresses, that is, less than ≈ 10 8 PA. Further, the temperature dependence of the average activation stress for sintering agrees well with analysis of bulk sintering data and creep data reported within the literature when analyzed in the context of a grain boundary dislocation nucleation rate limited kinetic model.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Quantifications of growth rate and particle size of magnesium silicate hydrate cement using small‐angle X‐ray scattering

The focus of recent research on cement binders has shifted to alternatives such as magnesium silicate hydrate (M-S-H) to mitigate the environmental impact of conventional concrete production. However, a detailed understanding of the morphology and growth rate of M-S-H, which informs its applicability as a cementing agent, is lacking. These properties are important as they influence cement hydration and consequently the development of macroscopic properties. Here, in this study, the particle size distribution and growth rate of M-S-H synthesized by mixing aqueous stock solutions were analyzed using synchrotron-based small-angle X-ray scattering. The reaction solutions contained dissolved magnesium and silicon with concentrations ranging from 2 to 500 mM, corresponding to saturation indices with respect to M-S-H endmembers, M 0.75 -S-H and M 1.50 -S-H, between 5.8 and 10.3. The precipitation rates of M-S-H align with affinity-based kinetic models and can be described by the following equations: For M 0.75 -S-H: 𝑟 M 0.75 −S−H =10 −9.54±0.70 ⁢(Ω$^{0.59±0.08}_{M_{0.75}⁢\textrm{SH}}$−1), and for M 1.50 -S-H: 𝑟 M 1.50 −S−H =10 −7.56±0.37 ⁢(Ω$^{0.34±0.04}_{M_{1.50}⁢\textrm{SH}}$−1), where 𝑟 is the rate in mol/L/s, and Ω represents the saturation ratio with respect to M-S-H. Moreover, we observed that changing the [Mg]/[Si] ratio affects particle size, a trend that cannot be fully explained by changes in the saturation index alone. These insights are critical for improving the predictability and control of cementation in concrete that utilizes M-S-H-based binders.

M-S-H↗
DOE OSTI2022

High-Fuel Loading Ignition Delay Time Characterization of Hydrogen/Natural Gas/Ammonia at Gas Turbine-Relevant Conditions Inside a High-Pressure Shock Tube

Abstract This study explores the combustion characterization of high-fuel percentage, air-diluted mixtures of H2 mixed with natural gas (NG) as well as mixtures of H2 and NH3 at temperatures and pressures relevant to turbine operating conditions (20–30 bar, 1000–1500 K). Lower temperatures (below 1070 K) exhibit preignition characteristics due to non-homogeneity. An attempt to mitigate these occurrences at high pressures is investigated using the constrained reaction volume (CRV) stage-filling technique. Due to the need to further refine the facility CRV stage-filling uncertainty, only higher temperature data will be interpreted at this time. The test conditions in this study closely replicate the temperatures, pressures, and mixtures that would be seen in hydrogen-powered gas turbines, making it the first to explore such conditions. The experimental IDTs were compared against the current state-of-the-art chemical kinetic models for mechanism validation. The current work will advance H2-powered turbines and aims to determine the optimum molecular ratio of H2 when mixed with natural gas.

Pierro, Michael↗
DOE OSTI2023

Thixotropy, antithixotropy, and viscoelasticity in hysteresis

Thixotropy, antithixotropy, and viscoelasticity are three types of time-dependent dynamics that involve fundamentally different underlying physical processes. Here, we show that the three dynamics exhibit different signatures in hysteresis by examining the fingerprints of the simplest thixotropic kinetic model, a new antithixotropic model that we introduce here, and the Giesekus model. We start by showing that a consistent protocol to generate hysteresis loops is a discrete shear-rate controlled ramp that begins and ends at high shear rates, rather than at low shear rates. Using this protocol, we identify two distinguishing features in the resulting stress versus shear rate loops. The first is the direction of the hysteresis loops: clockwise for thixotropy, but counterclockwise for viscoelasticity and antithixotropy. A second feature is achieved at high ramping rates where all responses lose hysteresis: the viscoelastic response shows a stress plateau at low shear rates due to lack of stress relaxation, whereas the thixotropic and antithixotropic responses are purely viscous with minimal shear thinning or thickening. We establish further evidence for these signatures by experimentally measuring the hysteresis of Laponite suspensions, carbon black suspensions, and poly(ethylene oxide) solutions, each representing a historically accepted example of each class of material behavior. The signatures measured in experiments are consistent with those predicted by the three models. This study reveals different fingerprints in hysteresis loops associated with thixotropy, antithixotropy, and viscoelasticity, which may be helpful in distinguishing the three time-dependent responses.

Mechanics↗
DOE OSTI2020

Ultimate Suppression of Thermal Transport in Amorphous Silicon Nitride by Phononic Nanostructure.

Engineering the thermal conductivity of amorphous materials is highly essential for the thermal management of future electronic devices. Here, we demonstrate the impact of ultrafine nanostructuring on the thermal conductivity reduction of amorphous silicon nitride (a-Si3N4) thin films, in which the thermal transport is inherently impeded by the atomic disorders. Ultrafine nanostructuring with feature sizes below 20 nm allows us to fully suppress contribution of the propagating vibrational modes (propagons), leaving only the diffusive vibrational modes (diffusons) to contribute to thermal transport in a-Si3N4. A combination of the phonon-gas kinetics model and the Allen-Feldmann theory reproduced the measured results without any fitting parameters. The thermal conductivity reduction was explained as extremely strong diffusive boundary scattering of both propagons and diffusIons. These findings give rise to substantial tunability of thermal conductivity of amorphous materials, which enables us to provide better thermal solutions in microelectronic devices.

Tambo, Naoki↗
DOE OSTI2017

Recipient-Biased Competition for an Intracellularly Generated Cross-Fed Nutrient Is Required for Coexistence of Microbial Mutualists

ABSTRACT Many mutualistic microbial relationships are based on nutrient cross-feeding. Traditionally, cross-feeding is viewed as being unidirectional, from the producer to the recipient. This is likely true when a producer’s waste, such as a fermentation product, has value only for a recipient. However, in some cases the cross-fed nutrient holds value for both the producer and the recipient. In such cases, there is potential for nutrient reacquisition by producer cells in a population, leading to competition against recipients. Here, we investigated the consequences of interpartner competition for cross-fed nutrients on mutualism dynamics by using an anaerobic coculture pairing fermentative Escherichia coli and phototrophic Rhodopseudomonas palustris . In this coculture, E. coli excretes waste organic acids that provide a carbon source for R. palustris . In return, R. palustris cross-feeds E. coli ammonium (NH 4 + ), a compound that both species value. To explore the potential for interpartner competition, we first used a kinetic model to simulate cocultures with varied affinities for NH 4 + in each species. The model predicted that interpartner competition for NH 4 + could profoundly impact population dynamics. We then experimentally tested the predictions by culturing mutants lacking NH 4 + transporters in both NH 4 + competition assays and mutualistic cocultures. Both theoretical and experimental results indicated that the recipient must have a competitive advantage in acquiring cross-fed NH 4 + to sustain the mutualism. This recipient-biased competitive advantage is predicted to be crucial, particularly when the communally valuable nutrient is generated intracellularly. Thus, the very metabolites that form the basis for mutualistic cross-feeding can also be subject to competition between mutualistic partners. IMPORTANCE Mutualistic relationships, particularly those based on nutrient cross-feeding, promote stability of diverse ecosystems and drive global biogeochemical cycles. Cross-fed nutrients within these systems can be either waste products valued by only one partner or nutrients valued by both partners. Here, we explored how interpartner competition for a communally valuable cross-fed nutrient impacts mutualism dynamics. We discovered that mutualism stability necessitates that the recipient have a competitive advantage against the producer in obtaining the cross-fed nutrient, provided that the nutrient is generated intracellularly. We propose that the requirement for recipient-biased competition is a general rule for mutualistic coexistence based on the transfer of intracellularly generated, communally valuable resources.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2017

Ionic composition of a humid air plasma under ionizing radiation

A kinetic model is proposed for ion–molecular processes involving charged particles of a humid air plasma produced by a fast electron beam. The model includes more than 600 processes involving electrons and 41 positive and 14 negative ions, including hydrated ions H{sub 3}O{sup +} (H{sub 2}O){sub n} and O{sub 2}{sup −}(H{sub 2}O){sub n} with n = 1, 2, …, 12. The energy costs of production of electron–ion pairs and electronic and vibrational (for water molecules, also rotational) excitation of molecules are calculated in nitrogen, oxygen, water vapor, air, and humid air. A method is proposed for calculating the energy costs in mixtures by the calculation data in pure gases. The evolution of the plasma composition is studied by the numerical solution of a system of 56 time-dependent balance equations for the number of charged particles of plasma by the fourth-order Runge–Kutta method. The steady-state composition of plasma is determined by solving nonlinear steady-state balance equations for the ionization rates of humid air from 10 to 10{sup 16} cm{sup –3}/s and the fraction of water molecules from 10{sup –3}% to 1.5%. It is established that, for water vapor content (the ratio of the number density of water molecules to the total number density of air molecules) of 0.015–1.5% in air at atmospheric pressure and room temperature, the main ion species are two types of positive ions H{sub 3}O{sup +} (H{sub 2}O){sub n} with the number of water molecules n = 5, 6 and three species of negative ions O{sub 2}{sup −}(H{sub 2}O){sub n} with n = 5, 8, 9.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

Analysis of the Range of Applicability of Thermodynamic Calculations in the Engineering of Nitride Fuel Elements

The domains of applicability of thermodynamic calculations in the engineering of nitride fuel are analyzed. Characteristic values of the following parameters, which affect directly the concentration equilibration time, are estimated: nuclide production rate; characteristic times to local equilibrium in the considered temperature range; characteristic time needed for a stationary temperature profile to be established; characteristic time needed for a quasi-stationary concentration field to be established on a scale comparable to the size of a fuel pellet. It is demonstrated that equilibrium thermodynamic calculations are suitable for estimating the chemical and phase composition of fuel. However, a two-layer kinetic model should be developed in order to characterize the transport processes in condensed and gaseous phases. The process of diffusive transport needs to be taken into account in order to determine the composition in the hot region at the center of a fuel element.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2018

Heating and Nonequilibrium Distributions of Ions in a Reverse Shock Wave of the SN 1987A Remnant

A hydrodynamical description of supernova remnants is based on the approximation of locally equilibrium particle distributions. Shock waves in supernova remnants at various stages of ejecta propagation are collisionless and form nonequilibrium particle distributions that relax slowly to quasiequilibrium distributions within a time longer than the hydrodynamic time. A kinetic model of the heating of ions behind the front of a shock wave in the SN 1987A remnant is considered with allowance for a complex chemical composition of the ejecta, and nonequilibrium distributions of ions in the vicinity of this shock wave are calculated. In addition to the quasi-Maxwellian peak, which determines the effective temperature of a given charge state of an ion, nonequilibriumdistributions of ions contain, in some cases, a suprethermal component, which may describe to the injection of ions in the process of cosmic-ray acceleration.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

Experimental and theoretical study of the near IR emission of xenon excited by a fast electron beam

Emission of xenon excited by a 120-keV electron beam at gas pressures of 100, 200, 500, and 760 Torr nm was studied experimentally and theoretically. More than 30 spectral lines were identified in the wavelength range of 750–1000 nm. A self-consistent kinetic model is developed to calculate the emission intensity of xenon atoms in the near IR range. The model includes balance equations for the number densities of electrons, ions and excimer molecules; equations for the populations of electron levels; and the Boltzmann equation for the low-energy part of the electron energy distribution function with a source of slow electrons. Excitation and ionization rates of xenon by the beam electrons and the energy spectrum of slow electrons are calculated by the Monte Carlo method. It is shown that, under these conditions, the main mechanism of xenon atom excitation is dissociative recombination of Xe{sub 3}{sup +} ions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2018

Simulation of the Plasma Afterglow in the Discharge Gap of a Subnanosecond Switch Based on an Open Discharge in Helium

The phenomenon of subnanosecond electrical breakdown in a strong electric field observed in an open discharge in helium at pressures of 6–20 Torr can be used to create ultrafast plasma switches triggering into a conducting state for a time shorter than 1 ns. To evaluate the possible repetition rate of such a subnanosecond switch, it is interesting to study the decay dynamics of the plasma remaining in the discharge gap after ultrafast breakdown. In this paper, a kinetic model based on the particle-in-cell Monte Carlo collision method is used to study the dynamics of the plasma afterglow in the discharge gap of a subnanosecond switch operating with helium at a pressure of 6 Torr. The simulation results show that the radiative, collisional-radiative, and three-body collision recombination mechanisms significantly contribute to the afterglow decay only while the plasma density remains higher than 10{sup 12} cm{sup −3}; the main mechanism of the further plasma decay is diffusion of plasma particles onto the wall. Therefore, the effect of recombination in the plasma bulk is observed only during the first 10–20 μs of the afterglow. Over nearly the same time, plasma electrons become thermalized. The afterglow time can be substantially reduced by applying a positive voltage U{sub c} to the cathode. Since diffusive losses are limited by the ion mobility, the additional ion drift toward the wall significantly accelerates plasma decay. As U{sub c} increases from 0 to +500 V, the characteristic time of plasma decay is reduced from 35 to 10 μs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Electron Transfer from Semiconductor Nanocrystals to Redox Enzymes

This review summarizes progress in understanding electron transfer from photoexcited nanocrystals to redox enzymes. The combination of the light-harvesting properties of nanocrystals and the catalytic properties of redox enzymes has emerged as a versatile platform to drive a variety of enzyme-catalyzed reactions with light. Transfer of a photoexcited charge from a nanocrystal to an enzyme is a critical first step for these reactions. This process has been studied in depth in systems that combine Cd-chalcogenide nanocrystals with hydrogenases. The two components can be assembled in close proximity to enable direct interfacial electron transfer or integrated with redox mediators to transport charges. Time-resolved spectroscopy and kinetic modeling have been used to measure the rates and efficiencies of the electron transfer. Electron transfer has been described within the framework of Marcus theory, providing insights into the factors that can be used to control the photochemical activity of these biohybrid systems. The range of potential applications and reactions that can be achieved using nanocrystal–enzyme systems is expanding, and numerous fundamental and practical questions remain to be addressed.

Chemistry↗
DOE OSTI2023

Platinum Surface Oxide and Oxygen Reduction Reaction Kinetics during Transient Fuel Cell Operation

Pt surface oxide coverage (POC) and Oxygen Reduction Reaction (ORR) were measured under a range of operating conditions in a Proton Exchange Membrane Fuel Cell (PEMFC). Sensitivity analysis suggests that adsorption of sulfonate groups of the ionomer on the Pt surface plays an important role in POC and ORR kinetics. Although both Pt oxide growth and ORR activity decay follow a logarithm of time behavior, it is found that ORR kinetics are affected by at least two different types of Pt oxide. A semi-empirical ORR kinetic model is proposed taking into account the effect of ionomer and Pt oxide types. The model is capable of providing a quantitative prediction of POC and ORR activity over a range of potential, temperature, relative humidity, and time that is relevant to normal PEMFC operation.

Electrochemistry↗