Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Kinetic”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

Comparison of hydrodynamic and semi-kinetic treatments for plasma flow along closed field lines

Hydrodynamic and semi-kinetic treatments of plasma flow along closed geomagnetic field lines are compared. The hydrodynamic treatment is based on a simplified 16-moment set of transport equations as the equations for the heat flows are not solved; the heat flows are treated heuristically. The semi-kinetic treatment is based on a particle code. The comparison deals with the distributions of the plasma density, flow velocity, and parallel and perpendicular temperatures as obtained from the two treatments during the various stages of the flow. In the kinetic treatment, the appropriate boundary condition is the prescription of the velocity distribution functions for the particles entering the flux tubes at the ionospheric boundaries; those particles leaving the system are determined by the processes occurring in the flux tube. The prescribed distributions are half-Maxwellian with temperature T(sub 0) and density n(sub 0). In the hydrodynamic model, the prescribed boundary conditions are on density (n(sub 0)), flow velocity (V(sub 0)) and temperature (T(sub 0). It was found that results from the hydrodynamic treatment critically depend on V(sub 0); for early stages of the flow this treatment yields results in good agreement with those from the kinetic treatment, when V(sub 0) = square root of (kT(sub 0)/2 (pi)m), which is the average velocity of particles moving in a given direction for a Maxwellian distribution. During this early stage, the flows developing form the conjugate ionospheres show some distinct transitions. For the first hour or so, the flows are highly supersonic and penetrate deep into the opposite hemispheres, and both hydrodynamics and kinetic treatments yield almost similar features. It is found that during this period heatflow effects are negligibly small. When a flow penetrates deep into the opposite hemisphere, the kinetic treatment predicts reflection and setting up of counterstreaming. In contrast, the hydrodynamic treatment yields a shock in the flow. The reasons for this difference in the two treatments is discussed, showing that in view of the relatively warm ions, the coupling of ion beams and the consequent shock formation in the offequatorial region are not likely due to the enhancements in the beam temperatures. The counterstreaming in the kinetic treatment and the shock in the hydrodynamic treatment first advance upward to the equator and then downward to the ionospheric boundary from where the flow originated. The transit time for this advancement is found to be about 1 hour for the respective models. After 2 hours or so, both models predict that the flows from the ionospheric boundaries are generally subsonic with respect to the local ion-sound speed. At late stages of the flow, when a substantial fraction of ions entering the flux tube begin to return back in the kinetic treatment, the hydrodynamic treatment with the boundary condition V(sub 0) = square root of (kT(sub 0)/2(pi)m) yields an over-refilling, and the choice of V(sub 0) becomes uncertain.

Singh, Nagendra↗
NASA NTRS2012

Polarization and Compressibility of Oblique Kinetic Alfven Waves

Even though solar wind, as a collisionless plasma, is properly described by the kineticMaxwell-Vlasov description, it can be argued that much of our understanding of solar wind observational data comes from an interpretation and numerical modeling which is based on a fluid description of magnetohydrodynamics. In recent years, there has been a significant interest in better understanding the importance of kinetic effects, i.e. the differences between the kinetic and usual fluid descriptions. Here we concentrate on physical properties of oblique kinetic Alfvn waves (KAWs), which are often recognized as one of the key ingredients in the solar wind turbulence cascade. We use three different fluid models with various degrees of complexity and calculate polarization and magnetic compressibility of oblique KAWs (propagation angle q = 88), which we compare to solutions derived from linear kinetic theory. We explore a wide range of possible proton plasma b = [0.1,10.0] and a wide range of length scales krL = [0.001,10.0]. It is shown that the classical isotropic two-fluid model is very compressible in comparison with kinetic theory and that the largest discrepancy occurs at scales larger than the proton gyroscale. We also show that the two-fluid model contains a large error in the polarization of electric field, even at scales krL 1. Furthermore, to understand these discrepancies between the two-fluid model and the kinetic theory, we employ two versions of the Landau fluid model that incorporate linear low-frequency kinetic effects such as Landau damping and finite Larmor radius (FLR) corrections into the fluid description. It is shown that Landau damping significantly reduces the magnetic compressibility and that FLR corrections (i.e. nongyrotropic contributions) are required to correctly capture the polarization.We also show that, in addition to Landau damping, FLR corrections are necessary to accurately describe the damping rate of KAWs. We conclude that kinetic effects are important even at scales which are significantly larger than the proton gyroscale krL 1.

Turbulence↗
NASA NTRS2017

Kinetic Scale Structure of Low-Frequency Waves and Fluctuations

The dissipation of solar wind turbulence at kinetic scales is believed to be important for the heating of the corona and for accelerating the wind. The linear Vlasov kinetic theory is a useful tool for identifying various wave modes, including kinetic Alfven, fast magnetosonic/whistler, and ion-acoustic (or kinetic slow), and their possible roles in the dissipation. However, the kinetic mode structure in the vicinity of ion-cyclotron modes is not clearly understood. The present paper aims to further elucidate the structure of these low-frequency waves by introducing discrete particle effects through hybrid simulations and Klimontovich formalism of spontaneous emission theory. The theory and simulation of spontaneously emitted low-frequency fluctuations are employed to identify and distinguish the detailed mode structures associated with ion-Bernstein modes versus quasi-modes. The spontaneous emission theory and simulation also confirm the findings of the Vlasov theory in that the kinetic Alfven waves can be defined over a wide range of frequencies, including the proton cyclotron frequency and its harmonics, especially for high-beta plasmas. This implies that these low-frequency modes may play predominant roles even in the fully kinetic description of kinetic scale turbulence and dissipation despite the fact that cyclotron harmonic and Bernstein modes may also play important roles in wave-particle interactions.

waves↗
DOE OSTI2025

Catalytic Resonance Theory: Experimental and Kinetic Interpretation of Programmable Catalysis

The complexity of programmable catalysts that change in the physical or electronic state on the time scale of a catalytic turnover results in uncertainty in the kinetic outcome of common measurements with dynamic catalyst experiments. To provide an interpretation of utility to experimentalists, a general programmable catalytic reaction exhibiting high turnover efficiency was simulated to understand the kinetic behavior of catalysts with varying frequency, binding energy, oscillation amplitude, and minimum oscillation binding energy at different temperatures and reactant gas pressures. The time-averaged catalytic rate for varying temperature in the form of an Arrhenius plot exhibited three distinct kinetic regimes, with the intermediate temperatures or applied frequencies exhibiting zero slope, indicative of barrier-free kinetic control. Transitions in experimentally measurable kinetic regimes with both temperature and applied catalyst oscillation frequency were associated with transitions in the sensitivity of reaction rate constants and degrees of rate control. Furthermore, these distinct kinetic macroscopic features specific to programmable catalysts were identified for observation in experimental dynamic kinetic measurements.

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

Oxidation kinetic modelling of Fe-based oxygen carriers for chemical looping applications: Impact of the topochemical effect

Chemical looping is a promising technology for fossil fuel utilization due to its high fuel conversion efficiency with in-situ CO 2 capture capability. Metal oxide are used as oxygen carriers (OCs) and circulate between a fuel reactor and an air reactor to perform reduction and oxidation reactions, respectively. In general, OC exiting the fuel reactor is not reduced fully to its metallic state due to many factors including carbon deposition and OC deactivation. Therefore, the effect of the initial reduction state on the OC oxidation in the air reactor is a significant parameter for consideration in developing the oxidation kinetic model. The objective of this work is to develop a physically significant kinetic model that applies to the oxidation of both initially fully and partially reduced OC with air. For this study, 1.5 mm Fe-based OC particle supported with TiO 2 was used as the model OC particle due to its complex multistep reaction nature. The oxidation kinetics were experimentally investigated in a thermogravimetric analyzer (TGA). Results indicate a significant difference in the oxidation rate profile for the OCs when oxidized from an initially fully reduced compared to an initially partially reduced state. Elemental mapping via energy-dispersive X-ray spectroscopy (EDS) reveals a shrinking-core type topochemical pattern across the OC particle, which was identified to be the cause of the dependency of kinetics on the initial reduction state. A generalized kinetic model was developed based on the observed shrinking-core behavior without presuming any rate-determining steps and experimentally validated over a broad range of temperatures (800–1000 °C) and oxygen concentrations (5, 7, 10, and 15 mol%). Impacts of particle porosity, size, and core-shell structure on the OC oxidation kinetics were analyzed in the developed oxidation kinetic model to suggest methods of improving the oxidation rate of the OC without modifying the chemical composition.

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

Product-specific kinetics reveal effect of epoxy crystallization on thermoset thermal degradation

Crystallization is a common problem for epoxy resins, which are ubiquitous in industrial and commercial use. Integration of crystallized epoxy monomers into cured thermosets has been shown to alter the thermosets’ final mechanical properties. However, no studies have investigated the impact of these crystals on the thermal stability of the thermoset. Here we investigate the degradation kinetics of a bisphenol F–based epoxy thermoset with and without crystallized monomers using Product-Specific Kinetic (PSK) analysis coupled with Evolved Gas Analysis-Mass Spectrometry (EGA-MS). PSK analysis revealed significant differences in evolved product ion kinetics, suggesting changes in the degradation kinetics between thermoset configurations. Further, it was concluded that early stages of degradation are influenced most by crystal presence due to the high concentration of unreacted epoxy monomers and lower cross-linking density of the cured network. After post-cure annealing, significant changes are observed in the degradation kinetics of thermosets without crystal inclusions. Conversely, post-cure annealing procedures of crystal integrated thermosets showed little change in the thermoset degradation kinetics across all conversion extents. These findings suggest that post-cure annealing of thermosets with crystals present at the onset does not alter the cross-linking density of the polymer network enough to significantly change the degradation kinetics. We hypothesize this is because the excess monomers from the melted crystals are unable to find suitable reaction sites for complete binding into the polymer network, which has direct implications for the material properties and final thermal stability of the thermoset.

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

High-Temperature 2D Optical Relaxation Visualizes Enhanced Oxygen Exchange Kinetics at Metal-Mixed Conducting Oxide Interfaces

Solid-state heterointerfaces are of interest for emergent local behavior that is distinct from either bulk parent compound. One technologically relevant example is the case of mixed ionic/electronic conductor (MIEC)–metal interfaces, which play an important role in electrochemistry. Metal–MIEC composite electrodes can demonstrate improved catalytic activity vs single-phase MIECs, improving fuel cell efficiency. Similarly, MIEC surface reaction kinetics are often evaluated using techniques that place metal current collectors in contact with the surface under evaluation, potentially altering the response vs the native surface. Here, techniques enabling direct and local in situ observation of the behavior at and around such heterointerfaces are needed. Here, we develop a spatially resolved optical transmission relaxation (2D-OTR) method providing continuous evaluation of local, high-temperature, controlled atmosphere defect kinetics across a ~1 cm2 sample area simultaneously in a contact-free manner. We apply it to observe the spatial variance of oxygen incorporation and evolution rates at ~525–620 °C, in response to step changes in oxygen partial pressure, on MIEC SrTi 0.65 Fe 0.35 O 3–x films as a function of distance from porous Pt and Au layers. Using this model geometry, we find significant enhancements in kinetics adjacent to the metals that decay over a few millimeter distance. To extract kinetic parameters, we fit the short-term optical data (initial portion of relaxations) with an exponential decay function appropriate for surface-exchange-limited kinetics, yielding apparent surface exchange coefficients (k chem ) with spatial resolution, decreasing with distance from the metal. To understand the kinetic processes governing the complete (long-term) optical relaxations, we performed COMSOL simulations, which demonstrated that a combination of laterally varying k chem and in-plane diffusion controls the observed kinetics over the full time range. Further support for spatially varying kchem comes from demonstrations of changing surface and bulk chemistry vs distance from the metal–MIEC interface, by X-ray photoelectron and optical absorption spectroscopies, respectively. Although microporous Pt and Au are not excellent electrodes in isolation, both metals exert a synergistic effect on the oxygen surface exchange rate in the presence of the mixed conducting film.

36 MATERIALS SCIENCE↗
DOE OSTI2024

The kinetic analog of the pressure–strain interaction

Energy transport in weakly collisional plasma systems is often studied with fluid models and diagnostics. However, the applicability of fluid models is limited when collisions are weak or absent, and using a fluid approach can obscure kinetic processes that provide key insights into the physics of energy transport. Kinetic diagnostics retain all of the information in 3D-3V phase space and thereby reach beyond the insights of fluid models to elucidate the mechanisms responsible for collisionless energy transport. In this work, we derive the Kinetic Pressure–Strain (KPS): a kinetic analog of the pressure–strain interaction, which is the channel between flow energy density and internal energy density in fluid models. Through two case studies of electron Landau damping, we demonstrate that the KPS diagnostic can elucidate kinetic mechanisms that are responsible for energy transport in this channel, just as the related field–particle correlation is known to identify kinetic mechanisms of transport between electromagnetic field energy density and kinetic energy density in particle flows. In addition, we show that resonant electrons play a major role in transferring energy between fluid flows and internal energy during the process of Landau damping.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Early stage growth of amorphous thin film: Average kinetics, nanoscale dynamics, and pressure dependence

We used the coherent grazing incidence small angle x-ray scattering technique to study the average kinetics and nanoscale dynamics during early-stage a-WSi 2 sputter deposition. The kinetic and dynamic properties are examined as a function of pressure, which is known to be a critical factor in determining final surface roughness. Surface growth kinetics and dynamics are characterized by time parameters extracted from the height-height structure factor and correlation functions. The roughness at a given length scale reaches a maximum before relaxing down to a steady state. The lateral length scale dependence and pressure dependence are then compared among the measured kinetics and dynamics time parameters. Surfaces grown at lower pressures are smoother, leading to longer correlation times. The time parameters to reach a dynamic steady state and a kinetic steady state show contrasting pressure dependence. A dynamic steady state is reached earlier than the kinetic steady state at high pressure. A more random deposition direction and lower kinetic energy at higher pressures can explain these phenomena, along with the hypothesis that larger nanoclusters form in vapor before arriving at the surface. Here, a continuum model is applied to simulate the overall behavior with mixed success.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Mass Transport Limitations and Kinetic Consequences of Corn Stover Deacetylation

Alkaline pretreatment of herbaceous feedstocks such as corn stover prior to mechanical refining and enzymatic saccharification improves downstream sugar yields by removing acetyl moieties from hemicellulose. However, the relationship between transport phenomena and deacetylation kinetics is virtually unknown for such feedstocks and this pretreatment process. Here, we report the development of an experimentally validated reaction–diffusion model for the deacetylation of corn stover. A tissue-specific transport model is used to estimate transport-independent kinetic rate constants for the reactive extraction of acetate, hemicellulose and lignin from corn stover under representative alkaline conditions (5–7 g L -1 NaOH, 10 wt% solids loadings) and at low to mild temperatures (4–70°C) selected to attenuate individual component extraction rates under differential kinetic regimes. The underlying transport model is based on microstructural characteristics of corn stover derived from statistically meaningful geometric particle and pore measurements. These physical descriptors are incorporated into distinct particle models of the three major anatomical fractions (cobs, husks and stalks) alongside an unsorted, aggregate corn stover particle, capturing average Feret lengths of 917–1239 μm and length-to-width aspect ratios of 1.8–2.9 for this highly heterogeneous feedstock. Individual reaction–diffusion models and their resulting particle model ensembles are used to validate and predict anatomically-specific and bulk feedstock performance under kinetic-controlled vs. diffusion-controlled regimes. In general, deacetylation kinetics and mass transfer processes are predicted to compete on similar time and length scales, emphasizing the significance of intraparticle transport phenomena. Critically, we predict that typical corn stover particles as small as ~2.3 mm in length are entirely diffusion-limited for acetate extraction, with experimental effectiveness factors calculated to be 0.50 for such processes. Debilitatingly low effectiveness factors of 0.021–0.054 are uncovered for cobs—implying that intraparticle mass transfer resistances may impair observable kinetic measurements of this anatomical fraction by up to 98%. These first-reported quantitative maps of reaction vs. diffusion control link fundamental insights into corn stover anatomy, biopolymer composition, practical size reduction thresholds and their kinetic consequences. These results offer a guidepost for industrial deacetylation reactor design, scale-up and feedstock selection, further establishing deacetylation as a viable biorefinery pretreatment for the conversion of lignocellulosics into value-added fuels and chemicals.

09 BIOMASS FUELS↗
DOE OSTI2026

Catalytic resonance theory for the kinetic signatures of multiple wavelength photocatalysis

Multi-step surface chemistry promoted with one or two wavelengths of light exhibited distinct kinetic signatures indicative of the light-adsorbate interactions and the number of photon-sensitive elementary steps. In this work, kinetic Monte Carlo simulations identified the kinetic response of general surface mechanisms to variation in the per-site photon flux of one or two wavelengths of incident light. Photocatalytic rates were described via a non-dimensional photon flux, identifying multiple kinetic regimes with unique degrees of rate control. Under photon-controlled kinetic conditions, maximum quantum yield and turnover frequency were obtained under constant illumination, while pulsing of one or more light sources was shown to exhibit slower rates and less efficient photocatalytic promotion due to the inherent dynamic nature of light, which comprises a stream of photons. Simulations provided distinct kinetic signatures in catalytic rates and Arrhenius kinetics for specific light-surface-adsorbate interactions corresponding to photocatalytic promotion of specific steps in surface chemistry.

catalytic resonance↗
DOE OSTI2020

Intrinsic Kinetic Limitations in Substituted Lithium-Layered Transition-Metal Oxide Electrodes

Substituted Li-layered transition-metal oxide (LTMO) electrodes such as Li x Ni y Mn z Co 1-y-z O 2 (NMC) and Li x Ni y Co 1-y-z Al z O 2 (NCA) show reduced first cycle Coulombic efficiency (90-87% under standard cycling conditions) in comparison with the archetypal Li x CoO 2 (LCO; similar to 98% efficiency). Focusing on Li x Ni 0.8 Co 0.15 Al 0.0.5 O 2 as a model compound, we use operando synchrotron X-ray diffraction (XRD) and nuclear magnetic resonance (NMR) spectroscopy to demonstrate that the apparent first-cycle capacity loss is a kinetic effect linked to limited Li mobility at x > 0.88, with near full capacity recovered during a potentiostatic hold following the galvanostatic charge- discharge cycle. This kinetic capacity loss, unlike many capacity losses in LTMOs, is independent of the cutoff voltage during delithiation and it is a reversible process. The kinetic limitation manifests not only as the kinetic capacity loss during discharge but as a subtle bimodal compositional distribution early in charge and, also, a dramatic increase of the charge-discharge voltage hysteresis at x > 0.88. 7 Li NMR measurements indicate that the kinetic limitation reflects limited Li transport at x > 0.86. Electrochemical measurements on a wider range of LTMOs including Li x (Ni,Fe) y Co 1-y O 2 suggest that 5% substitution is sufficient to induce the kinetic limitation and that the effect is not limited to Ni substitution. In this paper, we outline how, in addition to a reduction in the number of Li vacancies and shrinkage of the Li-layer size, the intrinsic charge storage mechanism (two-phase vs solid-solution) and localization of charge give rise to additional kinetic barriers in NCA and nonmetallic LTMOs in general.

25 ENERGY STORAGE↗
DOE OSTI2020

Two-fluid and kinetic transport physics of Kelvin–Helmholtz instabilities in nonuniform low-beta plasmas

Hall-magnetohydrodynamic (Hall-MHD) theory, two-fluid simulations, and kinetic simulations are used in this work to investigate the cross-field transport properties of Kelvin–Helmholtz instabilities in nonuniform low-beta collisionless plasmas. Hall-MHD analysis shows how the linear properties of the instability are modified by density gradients and magnetization. High-order accurate two-fluid and kinetic simulations, with complete dynamics of finite-mass electrons and ions, are applied to a suite of parameter cases to systematically assess the effects of diamagnetic drift, magnetization, charge separation, and finite Larmor motion. Initialization of exact two-species kinetic equilibria facilitates the study of isolated physical effects and enables detailed cross-comparisons between two-fluid and kinetic simulations, including for cases where ion gyroradii are comparable to gradient scale lengths. For nonuniform plasmas with significant space charge, the results of two-fluid and kinetic simulations are found to disagree with Hall-MHD predictions. Kelvin–Helmholtz instability growth rates, per unit shear, are shown to be smaller when ion diamagnetic drift and $E \times B$ drift are parallel and larger when the two drifts are antiparallel. The effect is attributed to polarization drift in the shear layer, which leads to redistribution of charge, alters the electric field that drives plasma advection, and consequently modifies growth rates. Instability-induced mass transport for different parameters is characterized in terms of the flux across the shear layer and a simplified diffusion model. Distribution functions from kinetic simulations are shown to deviate substantially from Maxwellian reconstructions, indicating the importance of kinetic physics during the nonlinear phase of the instability.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Accurate parameterization of the kinetic energy functional for calculations using exact-exchange

In this report, electronic structure calculations based on Kohn-Sham density functional theory (KSDFT) that incorporate exact exchange or hybrid functionals are associated with large computational expense, a consequence of the inherent cubic scaling bottleneck and large associated prefactor, which limits the length and time scales that can be accessed. Though orbital-free density functional theory (OFDFT) calculations scale linearly with system size and are associated with significantly smaller prefactor, they are limited by the absence of accurate density-dependent kinetic energy functionals. Therefore, the development of accurate density-dependent kinetic energy functionals is important for OFDFT calculations of large realistic systems. To this end, we propose a method to train kinetic energy functional models at the exact-exchange level of theory by using a dictionary of physically relevant terms that have been proposed in the literature in conjunction with linear or nonlinear regression methods to obtain the fitting coefficients. For our dictionary, we use gradient expansion of the kinetic energy, nonlocal models proposed in the literature and their nonlinear combinations, such as a model that incorporates spatial correlations between higher order derivatives of electron density at two points. The predictive capabilities of these models are assessed by using a variety of model one-dimensional systems that exhibit diverse bonding characteristics, such as a chain of eight hydrogens, LiF, LiH, C 4 H 2 , C 4 N 2 and C 3 O 2 . We show that by using data from model one-dimensional KSDFT calculations performed using the exact-exchange functional for only a few neutral structures, it is possible to generate models with high accuracy for charged systems as well as electron and kinetic energy densities during self-consistent field iterations. In addition, we show that it possible to learn both the orbital dependent terms, i.e., the kinetic energy and the exact-exchange energy, and models that incorporate additional nonlinearities in spatial correlations, such as a quadratic model, are needed to capture subtle features of the kinetic energy density that are present in exact-exchange-based KSDFT calculations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

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

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

42 ENGINEERING↗
DOE OSTI2021

Methods—Analyzing Electrochemical Kinetic Parameters in Deep Eutectic Solvents Using an Extended Butler-Volmer Equation

Deep eutectic solvents (DESs) are promising electrolytes for electrochemical redox reactions, which can be used in redox flow batteries (RFBs). However, in some systems like the Fe 2+/3+ redox reaction in ethaline, traditional Tafel-based kinetic analysis generates unreasonable kinetic parameters (i.e., large anodic/cathodic charge transfer coefficients ( α a , α c ) along with low exchange current densities ( i 0 )). This hinders a comprehensive kinetic and kinetic mechanism study. Here, we perform a detailed investigation of the Tafel analysis using a series of synthetic rotating disk electrode (RDE) data. We find the Tafel analysis only works well when i 0 < 0.57 i lim (limiting current density) in our scenario and leads to abnormal kinetic values once i 0 exceeds this limiting value. Thus, we propose an extended Butler-Volmer (ex-BV) analysis based on modern non-linear fitting techniques to obtain the actual kinetic parameters for such systems. The results show that this method fits the RDE data closely and generates reliable α a , α c and i 0 values, demonstrating that it is a good replacement for traditional Tafel analysis for kinetic studies in high-viscosity electrolytes such as DES systems.

Electrochemistry↗
DOE OSTI2025

Capability in Theory, Modeling, and Validation for a Range of Innovative Fusion Concepts using High-Fidelity Moment-Kinetic Models

A computational modeling capability is created and available to the fusion community to understand and design lower-cost and innovative fusion concepts. The approach uses high- fidelity kinetic, moment-kinetic, and moment models and includes sophisticated plasma- boundary interactions. A majority of fusion-relevant simulations are performed with magnetohydrodynamic models and hybrid particle-in-cell codes, with limited-fidelity electron and kinetic physics. However, in fusion configurations like Z-pinches, field-reversed- configurations, plasma jet magneto-inertial fusion, spinning mirrors, and others, kinetic effects (both electron and ions) are critical to understand the physics and design scaling into the highly kinetic regime of a burning fusion plasma. Furthermore, as present fusion machines move towards a burning plasma regime, liquid-metal blankets are needed to handle first-wall heat- flux, reduce erosion, and eventually for energy conversion and fuel breeding. The work performed under this ARPA-E BETHE Capability Team advances the state-of-the-art in modeling and understanding plasma dynamics in fusion devices and its coupling with liquid-metal dynamics. These are critical areas of research for fusion energy to become realizable. To address these complex problems, we have leveraged and extended computational capabilities through the code, Gkeyll (developed jointly with Princeton Plasma Physics Laboratory and academic partners), for kinetic and moment modeling of fusion plasmas. The Concept Teams supported by this Capability Team include the Wisconsin High-field Axisymmetric Mirror (WHAM), Centrifugal Mirror Experiment (CFME), Plasma-Jet Magneto- Inertial Fusion (PJMIF), and solid and liquid wall plasma-material interaction studies relevant to a number of fusion concepts including Zap Energy’s Z-pinch. This software is open-source and available to the fusion community as a high-fidelity tool for the design of lower-cost fusion experiments. 3D gyrokinetic simulations of WHAM are now possible for long enough time scales to understand the evolution of interchange instabilities. 3D multi-fluid simulations of CMFE at higher Mach numbers are now possible for detailed design iterations with the goal of stability. The state-of-the-art in understanding shock formation and shock mitigation regimes in merging liners for PJMIF have been furthered by our kinetic simulations. Our novel models and frameworks studying plasma-material interaction by incorporating wall emission for various solid wall materials of relevance to pulsed and steady fusion concepts have advanced the state-of-the-art in our understanding of particle fluxes, heat fluxes, and other quantities at cathodes and anodes. The results from this work may explain discrepancies between experimental and theoretical predictions of achieved current densities in pulsed concepts such as Z-pinches. Another significant contribution of this Capability Team is the development and deployment of a novel experimental platform, LEX (Liquid Electrode eXperiment), at Virginia Tech to understand liquid metal free-surface response to electromagnetic pulses. The novel experiments along with model validation quantified the effect of different materials and sizes of liquid metal droplets on the radiative power balance of fusion plasmas for pulsed concepts. Furthermore, these experiments provided mitigation strategies for violent liquid metal response for high current pulses as would be expected in fusion regimes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1983

Spectral study of wintertime kinetic energy of the Northern Hemisphere in the troposphere

Characteristics of the kinetic energy of wind fields at various pressure levels were analyzed, and significant wavenumbers in the wavenumber-frequency domain were identified. The nonlinear interaction terms of the kinetic energy equation were examined, and the distribution of the kinetic energy at the 850 mb, 500 mb, and 200 mb levels was calculated. A 5 deg latitude-longitude square grid was used, with NMC data for the 1975-1976 winter in the 20-60 deg N at 500 mb and 20-85 deg N for the 200 mb and 850 mb levels. The kinetic energy distribution was determined to be geography-dependent, with wavenumbers 6-9 westerly waves in the midfrequency range contributing significantly to kinetic energy maxima over the North Pacific and the east coast of North America. The contribution of the nonlinear interactions of these waves, which correspond to the longitudinal convergence of the kinetic energy flux, was found to be larger than the meridional convergence of the kinetic energy flux, and to occur mainly between 30-50 deg N. The nonlinear interactions were a negative contribution over the North Pacific at the 200 mb level.

Lee, H. N.↗