Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “kinetic parameters”

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,117 records · Page 62
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Origin of Rapid Delithiation In Secondary Particles Of LiNi 0.8 Co 0.15 Al 0.05 O 2 and LiNi y Mn z Co 1− y − z O 2 Cathodes

Abstract Most research on the electrochemical dynamics in materials for high‐energy Li‐ion batteries has focused on the global behavior of the electrode. This approach is susceptible to misleading analyses resulting from idiosyncratic kinetic conditions, such as surface impurities inducing an apparent two‐phase transformation within LiNi 0.8 Co 0.15 Al 0.05 O 2 . Here, nano‐focused X‐ray probes are used to measure delithiation operando at the scale of secondary particle agglomerates in layered cathode materials during charge. After an initial latent phase, individual secondary particles undergo rapid, stochastic, and largely uniform delithiation, which is in contrast with the gradual increase in cell potential. This behavior reproduces across several layered oxides. Operando X‐ray microdiffraction (‐XRD) leverages the relationship between Li content and lattice parameter to further reveal that rate acceleration occurs between Li‐site fraction ( x Li ) ≈0.9 and ≈0.5 for LiNi 0.8 Co 0.15 Al 0.05 O 2 . Physics‐based modeling shows that, to reproduce the experimental results, the exchange current density ( i 0 ) must depend on x Li , and that i 0 should increase rapidly over three orders of magnitude at the transition point. The specifics and implications of this jump in i 0 are crucial to understanding the charge‐storage reaction of Li‐ion battery cathodes.

25 ENERGY STORAGE↗
DOE OSTI2020

Inflationary stimulated Raman scattering in shock-ignition plasmas

In the shock-ignition inertial confinement fusion scheme, high-intensity lasers propagate through an inhomogeneous coronal plasma, driving a shock designed to cause fuel ignition. During the high-intensity ignitor laser pulse, in the long scale length coronal plasma, back-scattered stimulated Raman scattering (SRS) is likely to be in the kinetic regime. In this work, we use one-dimensional particle-in-cell simulations to show that there is a non-linear frequency shift caused by kinetic effects, resulting in the growth of SRS in an inhomogeneous plasma far exceeding the predictions of fluid theory, the so-called inflationary SRS or iSRS. We find that iSRS occurs over a wide range of density scale lengths relevant to shock-ignition and other directly-driven inertial confinement fusion schemes. The presence of iSRS in shock-ignition plasmas has implications for the theoretical gains from shock-ignition inertial confinement fusion. Here, we quantify the intensity threshold for the onset of iSRS for shock-ignition relevant parameters.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Cross-beam energy transfer in direct-drive ICF. I. Nonlinear and kinetic effects

Results are presented from a series of simulations examining the susceptibility of the cross-beam energy transfer (CBET) instability to nonlinear processes in the context of direct-drive inertial confinement fusion experiments on the OMEGA laser facility. These form the basis for the second paper of this series [A. G. Seaton, L. Yin, R. Follett, B. J. Albright, and A. Le, “Cross-beam energy transfer in direct-drive ICF. II. Theory and simulation of mitigation through increased laser bandwidth,” Phys. Plasmas 29, 042707 (2022)], where we examine the efficacy of increases in laser bandwidth at suppressing CBET. We choose laser and plasma conditions for the simulations that are favorable to CBET and promote nonlinearity. Through a comparison of outputs from the particle-in-cell code vector particle in cell (VPIC) and the linearized fluid code laser-plasma simulation environment (LPSE), a series of nonlinear effects have been identified in the kinetic simulations that include particle trapping, the two-ion wave decay, and ion-acoustic wave self-focusing. These effects produce time-dependent energy transfer, in contrast to the linearized fluid simulations in which a steady state is reached after an initial transient. Ion trapping is shown to allow for increased energy transfer relative to fluid simulations, with the remaining nonlinear processes acting to reduce the energy transfer. Nonlinear dynamics is contrasted for low- and high-intensity beams as well as between speckled and planar beams. For the parameters under consideration, beam profile has a significant effect on nonlinear dynamics, though the greatest sensitivity is to beam intensity.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

A classical, non-singular, bouncing universe

We present a model for a classical, non-singular bouncing cosmology without violation of the null energy condition (NEC). The field content is General Relativity plus a real scalar field with a canonical kinetic term and only renormalizable, polynomial-type self-interactions for the scalar field in the Jordan frame. The universe begins vacuum-energy dominated and is contracting at t=-∞. We consider a closed universe with a positive spatial curvature, which is responsible for the universe bouncing without any NEC violation. An Rφ 2 coupling between the Ricci scalar and the scalar field drives the scalar field from the initial false vacuum to the true vacuum during the bounce. The model is sub-Planckian throughout its evolution and every dimensionful parameter is below the effective-field-theory scale M P , so we expect no ghost-type or tachyonic instabilities. This model solves the horizon problem and extends co-moving particle geodesics to past infinity, resulting in a geodesically complete universe without singularities. We solve the Friedman equations and the scalar-field equation of motion numerically, and analytically under certain approximations.

Astronomy & Astrophysics↗
DOE OSTI2021

Fermi acceleration in relativistic collisionless plasma shocks correlates with anisotropic energy gains

Collisionless shocks generated by two colliding relativistic electron-positron plasma shells are studied using particle-in-cell simulations. Shocks are mediated by the Weibel instability (WI), and the kinetic energy of the fastest accelerated particles is found to be anisotropically modified by WI-induced electric fields. Specifically, we show that all particles interacting with the shock bifurcate into two groups based on their final relativistic Lorentz factor γ: slow (γ < γ bf ) and fast (γ < γ bf ), where γ bf is the bifurcation Lorentz factor that was found to be approximately twice the initial (upstream) Lorentz factor γ0. We have found that the kinetic energies of the slow particles are equally affected by the longitudinal and transverse components of the shock electric field, whereas the fast particles are primarily accelerated by the transverse field component.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

MOSCATO Development and Integration in Fiscal Year 2022

During FY21, we conducted ongoing development work for the MOSCATO (Molten Salt Chemistry and Transport) solver. The code development work primarily consisted of transitioning capabilities from the original version of the solver, which was written in OpenFOAM, into Nek5000. In doing so, a fast, highly parallelizable solver was created that is capable of complex chemistry and corrosion simulations for engineering-scale molten salt systems. The Nek5000 version of MOSCATO is now fully featured and capable of higher-fidelity simulations than were previously possible. Demonstration cases including a thermal convection loop have been simulated to test these new capabilities. We built upon the work for FY22 and improved the code from several different perspectives. First, we improved the user interface by adding a new component to the official Nek5000 input file (.par). This new part contains documents parameters like, salt properties (density, viscosity, Cp, thermal conductivity), diffusion coefficients, standard potential, etc. Second, we built a conversion script to extract salt properties from the MSTDB-TP salt database and write to MOSCATO input file. Third, we migrated the code to NekRS, which is the GPU branch of Nek5000 and suitable for next generation supercomputers. Verification and Validation (V&V) work was also continued in FY22. Two tasks were performed. The first V&V task involved the validation of the Poisson-Nernst-Planck equation solver and Butler-Volmer electrode kinetics, by comparing with numerical and experimental data about thermoelectric cells. The second task involved the comparisons to corrosion results from a thermal convection loop run during the MSRE era. Satisfactory agreement was obtained from both tasks.

Yuan, Haomin↗
DOE OSTI2024

Material Discovery and Design Principles of Perovskite Oxides for Reversible Solid Oxide Cells (R-SOC)

Reversible solid oxide cells (R-SOCs) are highly efficient devices for energy conversion and storage, capable of operating for both hydrogen utilization and production. In fuel cell mode, an R-SOC consumes hydrogen or natural gas to generate electricity, while in electrolysis mode, it produces hydrogen from steam. The discover of new materials with rapid oxygen surface exchange kinetics and enduring stability is crucial for the economically viable commercialization of R-SOCs. To facilitate this pursuit, we conducted extensive Density Functional Theory (DFT) calculations and developed Machine Learning (ML) models to predict critical catalytic properties essential for R-SOCs, such as oxygen surface exchange/diffusivity, and area-specific resistance (ASR). BaCoxFeyZrzO3-d(BFCZ)(x+y+z=1) emerged as a promising family of electrode materials with high activity and stability, validated through systematic experimental study. Moreover, a robust numerical multiphysics model was developed to optimize materials and microstructure parameters, providing the ability to predict the performance of functional R-SOCs.

Liu, Jian↗
DOE OSTI2025

Optimizing ensemble NV − spin properties of fluorescent diamond microparticles by systematic low pressure high temperature annealing

Low pressure high temperature annealing is a means for driving nitrogen and defect diffusion in diamond to reduce internal lattice damage without the need for technically complicated high-pressure cells. Herein, we perform a systematic time (5, 15, and 30 min) and temperature (1200 °C–1800 °C) study of effects of low-pressure high temperature annealing on photoluminescence, spin concentrations, and spin relaxation properties of NV centers in ca. 3 μm synthetic type 1b diamond particles. Annealing in the temperature range of ca. 1400 °C–1700 °C for even 5 min leads to a higher optically detected magnetic resonance contrast as compared to standard annealing at 900 °C for 2 h. Particles annealed at 1700 °C for 5 min exhibit a contrast close to about 13% as compared to about 9% for those annealed at 900 °C for 2 h. A reduction in the zero-field splitting strain parameter from E ≈ 4.5 MHz to ≈ 2.5 MHz and spectral linewidth from Δν ≈ 7 MHz to ≈ 4 MHz are observed even after 5 min annealing at 1700 °C. Improvements in these spectral parameters resulted in a roughly 2-fold reduction in the noise level of temperature monitoring experiment utilizing an ensemble of NV centers in the particles. Annealing in the temperature range of 1600 °C for 15 or 30 min or 1700 °C for 5 min resulted in NV T 1 relaxation times approaching ca. 5 ms typically observed for bulk diamond. Quantitative electron paramagnetic resonance (EPR) allowed for estimations of thermal activation energies of paramagnetic center annihilation. Monitoring the primary defect concentration (P1 and other defects with half integer spins) and utilizing second order kinetic modeling, an activation energy of 3.63 ± 0.28 eV was estimated. Alternatively, using the NV half field EPR signal and first order kinetic modeling, a similar activation energy 3.89 ± 0.29 eV was estimated.

NV centers↗
DOE OSTI2022

Non-Idealities in Lab-Scale Kinetic Testing: A Theoretical Study of a Modular Temkin Reactor

The Temkin reactor can be applied for industrial relevant catalyst testing with unmodified catalyst particles. It was assumed in the literature that this reactor behaves as a cascade of continuously stirred tank reactors (CSTR). However, this assumption was based only on outlet gas composition or inert residence time distribution measurements. The present work theoretically investigates the catalytic CO2 methanation as a test case on different catalyst geometries, a sphere, and a ring, inside a single Temkin reaction chamber under isothermal conditions. Axial gas-phase species profiles from detailed computational fluid dynamics (CFD) are compared with a CSTR and 1D plug-flow reactor (PFR) model using a sophisticated microkinetic model. In addition, a 1D chemical reactor network (CRN) model was developed, and model parameters were adjusted based on the CFD simulations. Whereas the ideal reactor models overpredict the axial product concentrations, the CRN model results agree well with the CFD simulations, especially under low to medium flow rates. This study shows that complex flow patterns greatly influence species fields inside the Temkin reactor. Although residence time measurements suggest CSTR-like behavior, the reactive flow cannot be described by either a CSTR or PFR model but with the developed CRN model.

Wehinger, Gregor D. (ORCID:0000000217743391)↗
DOE OSTI2020

Validation of the two-region Rossi-alpha model for reflected assemblies

Nuclear criticality safety, nonproliferation and safeguards, emergency response, and stockpile stewardship utilize estimates of the k eff multiplication factor. The value of k eff cannot be directly measured, but it can be inferred from the prompt neutron period. One modality of measuring the prompt neutron period is the Rossi-alpha method, which is predicated on fitting a histogram of detection time differences due to the nonrandom temporal distribution of same-fission-chain neutrons. Recent works have developed the motivation and theory to expand traditional one-region point kinetic models resulting in one-exponential histogram fits to two-region models that result in two-exponential fits. This paper validates the new two-region model using organic scintillator measurements of copper-reflected weapons-grade plutonium (0:83 ≤ k eff ≤ 0:94) and high-density- polyethylene-reflected highly enriched uranium (0:73 ≤ k eff ≤ 0:95). Furthermore, the results show that more thermal systems have shorter prompt neutron periods in the core region due to increased induced-fission probabilities for moderated neutrons. Here, a new parameter introduced by the two- region model is also shown to be correlated to the amount of reflection, and may be used to infer assembly properties such as type and amount of reflector.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2022

Accurate quantification of lattice temperature dynamics from ultrafast electron diffraction of single-crystal films using dynamical scattering simulations

In ultrafast electron diffraction (UED) experiments, accurate retrieval of time-resolved structural parameters, such as atomic coordinates and thermal displacement parameters, requires an accurate scattering model. Unfortunately, kinematical models are often inaccurate even for relativistic electron probes, especially for dense, oriented single crystals where strong channeling and multiple scattering effects are present. This article introduces and demonstrates dynamical scattering models tailored for quantitative analysis of UED experiments performed on single-crystal films. As a case study, we examine ultrafast laser heating of single-crystal gold films. Comparison of kinematical and dynamical models reveals the strong effects of dynamical scattering within nm-scale films and their dependence on sample topography and probe kinetic energy. Applying to UED experiments on an 11 nm thick film using 750 keV electron probe pulses, the dynamical models provide a tenfold improvement over a comparable kinematical model in matching the measured UED patterns. Also, the retrieved lattice temperature rise is in very good agreement with predictions based on previously measured optical constants of gold, whereas fitting the Debye–Waller factor retrieves values that are more than three times lower. Altogether, these results show the importance of a dynamical scattering theory for quantitative analysis of UED and demonstrate models that can be practically applied to single-crystal materials and heterostructures.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Self-similar Reynolds-averaged mechanical–scalar turbulence models for Rayleigh–Taylor mixing induced by power-law accelerations in the small Atwood number limit

Analytical self-similar solutions to two-, three-, and four-equation Reynolds-averaged mechanical–scalar turbulence models describing turbulent Rayleigh–Taylor mixing driven by a temporal power-law acceleration are derived in the small Atwood number (Boussinesq) limit. The solutions generalize those previously derived for constant acceleration Rayleigh–Taylor mixing for models based on the turbulent kinetic energy K and its dissipation rate ε, together with the scalar variance S and its dissipation rate χ [O. Schilling, “Self-similar Reynolds-averaged mechanical–scalar turbulence models for Rayleigh–Taylor, Richtmyer–Meshkov, and Kelvin–Helmholtz instability-induced mixing in the small Atwood number limit,” Phys. Fluids 33, 085129 (2021)]. The turbulent fields are expressed in terms of the model coefficients and power-law exponent, with their temporal power-law scalings obtained by requiring that the self-similar equations are explicitly time-independent. Mixing layer growth parameters and other physical observables are obtained explicitly as functions of the model coefficients and parameterized by the exponent of the power-law acceleration. Values for physical observables in the constant acceleration case are used to calibrate the two-, three-, and four-equation models, such that the self-similar solutions are consistent with experimental and numerical simulation data corresponding to a canonical (i.e., constant acceleration) Rayleigh–Taylor turbulent flow. The calibrated four-equation model is then used to numerically reconstruct the mean and turbulent fields, and turbulent equation budgets across the mixing layer for several values of the power-law exponent. Finally, the reference solutions derived here can be used to understand the model predictions for strongly accelerated or decelerated Rayleigh–Taylor mixing in the large Reynolds number limit.

42 ENGINEERING↗
DOE OSTI2025

KIPM Detector Characterization in QUIET

Kinetic Inductance Phonon‑Mediated (KIPM) detectors are superconducting microwave resonators that sense energy depositions in the substrate through a transient shift in their resonant frequency. In this work, I installed new devices in the QUIET cryogenic facility and used a network analyzer to drive them across a range of power levels. At each setting, I measured the transmitted signal and fit it to a simple notch‑filter model, allowing us to extract key performance metrics. Despite the increased noise at low power, the model consistently captured the resonance behavior. I first measured the resonance parameters as a function of RF drive power, then repeated the measurement across a range of temperatures to test whether the device behaved like an MKID. While the notch-filter fits captured the resonance reliably, the observed temperature dependence did not match the expected MKID signature — indicating that these first in-house fabricated chips are not yet functioning as true MKIDs. Overall, this reframes our work as a study of resonance parameter extraction and identification, pointing towards future iterations of device design and fabrication to get the true MKID response.

Savitala, Akash [U. Washington, Seattle (main); Fe↗
DOE OSTI2021

Early- and late-time dynamics of laser-produced plasmas by combining emission and absorption spectroscopy

Optical emission spectroscopy and laser absorption spectroscopy are two highly complementary techniques for diagnosing plasmas, specifically in the determination of important parameters such as temperature and density. In this work, we characterized a laser-produced plasma by combining emission and tunable absorption spectroscopy. Plasmas are produced by focusing pulses from a frequency quadrupled Nd:YAG laser onto a glass target containing a minor concentration of aluminum in an ambient atmosphere of ~100 Torr N 2 . Here, emission spectroscopy is used to investigate the line behavior of Al and Ca, electron density, and excitation temperature at times = 5 μs, while absorption spectroscopy is used to investigate the line behavior of Al and Ca, species density, and kinetic temperature at times = 5-100 μs. A comparison is made between the excitation temperature measured from emission using the Boltzmann plot method and the kinetic temperature measured from absorption; the implications of this comparison on local thermodynamic equilibrium and line-of-sight averaging are discussed. By combining the two methods of optical emission and laser absorption spectroscopy, a comprehensive characterization of the plasma over the entire lifetime of its evolution is possible, leading to a greater understanding of the plasma behavior including particle condensation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Thermodynamics and collisionality in firehose-susceptible high- β plasmas

We study the evolution of collisionless plasmas that, due to their macroscopic evolution, are susceptible to the firehose instability, using both analytic theory and hybrid-kinetic particle-in-cell simulations. We establish that, depending on the relative magnitude of the plasma β, the characteristic time scale of macroscopic evolution and the ion-Larmor frequency, the saturation of the firehose instability in high-β plasmas can result in three qualitatively distinct thermodynamic (and electromagnetic) states. By contrast with the previously identified ‘ultra-high-beta’ and ‘Alfvén-inhibiting’ states, the newly identified ‘Alfvén-enabling’ state, which is realised when the macroscopic evolution time τ exceeds the ion-Larmor frequency by a β-dependent critical parameter, can support linear Alfvén waves and Alfvénic turbulence because the magnetic tension associated with the plasma’s macroscopic magnetic field is never completely negated by anisotropic pressure forces. We characterise these states in detail, including their saturated magnetic-energy spectra. The effective collision operator associated with the firehose fluctuations is also described; we find it to be well approximated in the Alfvén-enabling state by a simple quasi-linear pitch-angle scattering operator. The box-averaged collision frequency is ν eff ∼ β/τ, in agreement with previous results, but certain subpopulations of particles scatter at a much larger (or smaller) rate depending on their velocity in the direction parallel to the magnetic field. Our findings are essential for understanding low-collisionality astrophysical plasmas including the solar wind, the intracluster medium of galaxy clusters and black hole accretion flows. We show that all three of these plasmas are in the Alfvén-enabling regime of firehose saturation and discuss the implications of this result.

astrophysical plasmas↗
DOE OSTI2023

Changing Directions: Influence of Ligand Electronics on the Directionality and Kinetics of Photoinduced Charge Transfer in Cu(I)Diimine Complexes

A key challenge to the effective utilization of solar energy is to promote efficient photoinduced charge transfer, specifically avoiding unproductive, circuitous electron-transfer pathways and optimizing the kinetics of charge separation and recombination. We hypothesize that one way to address this challenge is to develop a fundamental understanding of how to initiate and control directional photoinduced charge transfer, particularly for earth-abundant first-row transition-metal coordination complexes, which typically suffer from relatively short excited-state lifetimes. Here, in this paper, we report a series of functionalized heteroleptic copper(I)bis(phenanthroline) complexes, which have allowed us to investigate the directionality of intramolecular photoinduced metal-to-ligand charge transfer (MLCT) as a function of the substituent Hammett parameter. Ultrafast transient absorption suggests a complicated interplay of MLCT localization and solvent interaction with the Cu(II) center of the MLCT state. This work provides a set of design principles for directional charge transfer in earth-abundant complexes and can be used to efficiently design pathways for connecting the molecular modules to catalysts or electrodes and integration into systems for light-driven catalysis.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

A theoretical investigation of the hydrolysis of uranium hexafluoride: the initiation mechanism and vibrational spectroscopy

Depleted uranium hexafluoride (UF 6 ), a stockpiled byproduct of the nuclear fuel cycle, reacts readily with atmospheric humidity, but the mechanism is poorly understood. Here we compare several potential initiation steps at a consistent level of theory, generating underlying structures and vibrational modes using hybrid density functional theory (DFT) and computing relative energies of stationary points with double-hybrid (DH) DFT. A benchmark comparison is performed to assess the quality of DH-DFT data using reference energy differences obtained using a complete-basis-limit coupled-cluster (CC) composite method. The associated large-basis CC computations were enabled by a new general-purpose pseudopotential capability implemented as part of this work. Dispersion-corrected parameter-free DH-DFT methods, namely PBE0-DH-D3(BJ) and PBE-QIDH-D3(BJ), provided mean unsigned errors within chemical accuracy (1 kcal mol -1 ) for a set of barrier heights corresponding to the most energetically favorable initiation steps. The hydrolysis mechanism is found to proceed via intermolecular hydrogen transfer within van der Waals complexes involving UF 6 , UF 5 OH, and UOF 4 , in agreement with previous studies, followed by the formation of a previously unappreciated dihydroxide intermediate, UF 4 (OH) 2 . The dihydroxide is predicted to form under both kinetic and thermodynamic control, and, unlike the alternate pathway leading to the UO 2 F 2 monomer, its reaction energy is exothermic, in agreement with observation. Finally, harmonic and anharmonic vibrational simulations are performed to reinterpret literature infrared spectroscopy in light of this newly identified species.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Computational Analysis of the Effect of Structured Packing Design on Absorption Column Hydrodynamics for Post-Combustion Carbon Capture Applications

Solvent based post-combustion carbon capture technologies have a potential for reducing carbon emissions from fossil-fuel-fired power plants and industrial sources where CO2 emissions are inherently harder to mitigate, such as steel or cement industries. A prominent technology to achieve this is by retrofitting absorption columns to the existing infrastructure. While these systems have been among the less costly alternatives for carbon capture, they still impose a considerable energy penalty to the operation of power plants or industrial facilities. The optimization of CO2 capture rate in solvent-based absorption process is complex as it depends on several factors including CO2 solubility, solvent reaction kinetics and temperature effects on the solubility, reaction rates, surface tension, and thermophysical properties of the solvent and the flue gas. The overall heat and mass transfer also depends on the hydrodynamics, which in turn, is affected by the packing geometry. In the current work, we systematically quantify the effects the design of the structured packing has on the column hydrodynamics, by performing detailed CFD simulations for different geometrical configurations and operating conditions. We then obtain relationships between the key hydrodynamic metrics, such as liquid holdup, interfacial area, wetted area, and pressure drop to the parameters defining the packing geometries and identify new more effective packing designs for the given operating conditions.

Shah, Yash Girish↗