Search NASA⌕ Search

SEARCH · Search NASA

Results for “Kinetics 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 847 records · Page 47

Kinetic corrections from analytic non-Maxwellian distribution functions in magnetized plasmas

In magnetized plasma physics, almost all developed analytic theories assume a Maxwellian distribution function (MDF) and in some cases small deviations are described using the perturbation theory. The deviations with respect to the Maxwellian equilibrium, called kinetic effects, are required to be taken into account especially for fusion reactor plasmas. Generally, because the perturbation theory is not consistent with observed steady-state non-Maxwellians, these kinetic effects are numerically evaluated by very central processing unit (CPU)-expensive codes, avoiding the analytic complexity of velocity phase space integrals. We develop here a new method based on analytic non-Maxwellian distribution functions constructed from non-orthogonal basis sets in order to (i) use as few parameters as possible, (ii) increase the efficiency to model numerical and experimental non-Maxwellians, (iii) help to understand unsolved problems such as diagnostics discrepancies from the physical interpretation of the parameters, and (iv) obtain analytic corrections due to kinetic effects given by a small number of terms and removing the numerical error of the evaluation of velocity phase space integrals. This work does not attempt to derive new physical effects even if it could be possible to discover one from the better understandings of some unsolved problems, but here we focus on the analytic prediction of kinetic corrections from analytic non-Maxwellians. As applications, examples of analytic kinetic corrections are shown for the secondary electron emission, the Langmuir probe characteristic curve, and the entropy. This is done by using three analytic representations of the distribution function: the Kappa distribution function, the bi-modal or a new interpreted non-Maxwellian distribution function (INMDF). The existence of INMDFs is proved by new understandings of the experimental discrepancy of the measured electron temperature between two diagnostics in JET. As main results, it is shown that (i) the empirical formula for the secondary electron emission is not consistent with a MDF due to the presence of super-thermal particles, (ii) the super-thermal particles can replace a diffusion parameter in the Langmuir probe current formula, and (iii) the entropy can explicitly decrease in presence of sources only for the introduced INMDF without violating the second law of thermodynamics. Moreover, the first order entropy of an infinite number of super-thermal tails stays the same as the entropy of a MDF. The latter demystifies the Maxwell's demon by statistically describing non-isolated systems.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Estimating parameters from rotating ring disc electrode measurements

Rotating ring disc electrode (RRDE) experiments are a classic tool for investigating kinetics of electrochemical reactions. Several standardized methods exist for extracting transport parameters and reaction rate constants using RRDE measurements. Here in this work, we compare some approximate solutions to the convective diffusion used popularly in the literature to a rigorous numerical solution of the Nernst-Planck equations coupled to the three dimensional flow problem. In light of these computational advancements, we explore design aspects of the RRDE that will help improve sensitivity of our parameter estimation procedure to experimental data. We use the oxygen reduction in acidic media involving three charge transfer reactions and a chemical reaction as an example, and identify ways to isolate reaction currents for the individual processes in order to accurately estimate the exchange current densities.

30 DIRECT ENERGY CONVERSION↗

Critical Parameter Identification of Fuel-Cell Models Using Sensitivity Analysis

Numerical modeling has been a vital tool in proton-exchange-membrane fuel-cell (PEMFC) analysis; however, the predictive capabilities of these models depend on the input physical parameters, several of which are either not experimentally measured or have large scatter in measured values. This article presents an uncertainty propagation-based sensitivity analysis to identify the model parameters that impact the model predictions most. A comprehensive 2-D membrane electrode assembly (MEA) model is used to perform local sensitivity analysis at multiple operating conditions, which encompass the range of environments and operating conditions a cell can encounter. While at lower humidities, cathode kinetics and membrane-ohmic-loss related parameters are crucial, gas transport and porous-media saturation behavior are more important at humidified conditions. Several of these findings are different from previous studies presented in literature. Identifying the crucial parameters helps focus future material and cell optimization studies as well as experimental studies to quantify these parameters with higher accuracy.

25 ENERGY STORAGE↗

An Integrated Simulation of Multiple-Pass U-10Mo Alloy Hot Rolling and Static Recrystallization

To achieve a desired microstructure and minimize the thickness variation in rolled foils, researchers must understand the effects of foil fabrication process variables on microstructure evolution. We developed an integrated simulation of deformation and recrystallization that employs the finite element method (FEM) and the kinetic Monte Carlo (KMC) Potts model, respectively, to investigate microstructure evolution during multiple-pass hot rolling and heat treatment in polycrystalline U-10Mo fuel. Scanning electron microscopy and electron backscatter diffraction images of microstructures were directly used as input in FEM calculation of deformation, and the calculated strains were used to determine the driving force of nucleation and growth of recrystallized grains in the Potts model. Grain structures predicted by the Potts model were used to update the grain structure and material properties for FEM. Simulation alternated between FEM and the Potts model to simulate grain structure evolution during multiple rolling and heat treatments. The initial model parameters were determined by benchmarking the recrystallization kinetics against experimental data. Then, the model was applied to predict the grain structure evolution. Results showed that our model can capture the coupling between deformation and recrystallization and can quantitatively reproduce the observed U-10Mo recrystallization and grain growth kinetics. The simulation results demonstrated that the developed model can predict U-10Mo grain structures as a function of initial microstructure and foil fabrication parameters.

36 MATERIALS SCIENCE↗

Kinetics and Thermochemistry of the Cl((sup 2)P(sub J)) + C2Cl4 Association Reaction

A laser flash photolysis-resonance fluorescence technique has been employed to study the kinetics of the Cl(sup 2)P(sub j) + C2Cl4 association reaction as a function of temperature (231-390 K) and pressure (3-700 Torr) in nitrogen buffer gas. The reaction is found to be in the falloff regime between third and second order over the range of conditions investigated, although the second-order limit is approached at the highest pressures and lowest temperatures. At temperatures below 300 K, the association reaction is found to be irreversible on the experimental time scale of approximately 20 m-s. The kinetic data at T is less than 300 K have been employed to obtain falloff parameters in a convenient format for atmospheric modeling. At temperatures above 330 K, reversible addition is observed, thus allowing equilibrium constants for C2Cl5 formation and dissociation to be determined. Second- and third-law analyses of the equilibrium data lead to the following thermochemical parameters for the association reaction: Delta-H(298) = -18.1 +/- 1.3 kcal/mol, Delta-H(0) = -17.6 +/- 1.3 kcal/mol, and Delta-S(298) = -27.7 +/- 3.0 cal/mol.K. In conjunction with the well-known heats of formation of Cl((sup 2)P(sub j)) and C2Cl4 the above Delta-H values lead to the following heats of formation for C2Cl5, at 298 and 0 K: Delta-H(f,298) = 8.0 +/- 1.3 kcal/mol and Delta-H(f,0) = 8.1 +/- 1.5 kcal/mol. The kinetic and thermochemical parameters reported above are compared with other reported values, and the significance of reported association rate coefficients for understanding tropospheric chlorine chemistry is discussed.

Nicovich, J. M.↗

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↗

Approach to Evaluating Reorganization Energies of Interfacial Electrochemical Reactions

Reaction rate coefficients for electron-transfer processes at the electrode–electrolyte interface are commonly estimated by using the Butler–Volmer equation, but their values are inaccurate beyond a few tenths of volts of overpotential. The Marcus–Hush–Chidsey (MHC) formalism yields correct asymptotic behavior of the rate coefficients vs applied overpotential but has complex dependencies on the redox system’s intrinsic parameters, which can be difficult to model or measure. In this work, we bridge the two kinetics formalisms to estimate the reorganization energy, one of the important parameters for the MHC formalism, and investigate its dependence on other intrinsic parameters such as activation barriers, electronic coupling strength, and the density of states of the electrode surface. We examine the sensitivity of the reorganization energy to these parameters, establish some general relationships for accurately predicting rate coefficients using the MHC formalism over a wide range of applied overpotentials, and compare this approach to calculating MHC rate constants with other empirical approaches for the mechanisms of CO 2 reduction on different metal electrode surfaces.

Butler−Volmer↗

Ghosts- and tachyon-free regions of the Randall-Sundrum model parameter space

Model building within the Randall-Sundrum (RS) framework generally involves placing the Standard Model fields in the bulk. Such fields may possess non-zero values for their associated brane-localized kinetic terms (BLKTs) in addition to possible bulk mass parameters. Here in this paper we clearly identify the regions of the RS model parameter space where the presence of bulk mass terms and BLKTs yield a setup which is free from both ghost and tachyon instabilities. Such physically acceptable parameter space regions can then be used to construct realistic and phenomenologically viable RS models.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Direct Numerical Simulation of Involute Channel Turbulence

A direct numerical simulation (DNS) study was performed on turbulent flow in the high flux isotope reactor involute channel geometry to develop a numerical database and determine the differences compared with a flat parallel channel. The varying channel curvature along the walls was studied for differences in mean profiles. Parameters of interest include streamwise velocity, turbulent kinetic energy (TKE), and turbulence dissipation rate, as well as Reynolds stresses and turbulence transport terms. Profile sampling was carried out at 10 locations along the span of the involute. Additional DNS studies were performed on smaller domains of comparable curvature to the involute domain: a high curvature channel (high circular), a low curvature channel (low circular), and a flat channel (flat). Here, each of these four cases was compared against each other and to other DNS studies performed on parallel flows. The results indicate that the bulk involute channel flow does not differ significantly from a flat parallel channel flow and that the curvature of the walls does not significantly alter the mean flow parameters. However, the regions of the involute channel near the side walls exhibit relatively low magnitude twin recirculation structures driven toward the side walls from the centerline of the channel, which warrants further study.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

An adaptive least-squares global sensitivity method and application to a plasma-coupled combustion prediction with parametric correlation

We introduce an efficient non-intrusive surrogate-based methodology for global sensitivity analysis and uncertainty quantification. Modified covariance-based sensitivity indices (mCov-SI) are defined for outputs that reflect correlated effects. The overall approach is applied to simulations of a complex plasma-coupled combustion system with disparate uncertain parameters in sub-models for chemical kinetics and a laser-induced breakdown ignition seed. The surrogate is based on an Analysis of Variance (ANOVA) expansion, such as widely used in statistics, with orthogonal polynomials representing the ANOVA subspaces and a polynomial dimensional decomposition (PDD) representing its multi-dimensional components. The coefficients of the PDD expansion are obtained using a least-squares regression, which both avoids the direct computation of high-dimensional integrals and affords an attractive flexibility in choosing sampling points. This facilitates importance sampling using a Bayesian calibrated posterior distribution, which is fast and thus particularly advantageous in common practical cases, such as our large-scale demonstration, for which the asymptotic convergence properties of polynomial expansions cannot be realized due to computation expense. Effort, instead, is focused on efficient finite-resolution sampling. Standard covariance-based sensitivity indices (Cov-SI) are employed to account for correlation of the uncertain parameters. Magnitude of Cov-SI is unfortunately unbounded, which can produce extremely large indices that limit their utility. Alternatively, mCov-SI are then proposed in order to bound this magnitude $\in$ [0, 1]. The polynomial expansion is coupled with an adaptive ANOVA strategy to provide an accurate surrogate as the union of several low-dimensional spaces, avoiding the typical computational cost of a high-dimensional expansion. It is also adaptively simplified according to the relative contribution of the different polynomials to the total variance. Finally, the approach is demonstrated for a laser-induced turbulent combustion simulation model, which includes parameters with correlated effects.

42 ENGINEERING↗

Epitaxial growth and atomic arrangement in Fe2CrO4 on crystal symmetry matched (001) MgAl2O4

We report on the chemical, structural, and optical properties of molecular beam epitaxy (MBE) synthesized thin films of multifunctional Fe2CrO4 on (001)-oriented MgAl2O4 (MAO). Substrate temperatures near 500 °C are required to obtain smooth films with the expected out-of-plane (OOP) lattice parameter. Mg diffusion from the MAO substrate is kinetically suppressed up to 500 °C. We discuss anti-phase boundaries in symmetry matched epitaxial systems. This research provides new insight into the epitaxial growth and crystalline properties of crystal-symmetry-matched Fe2CrO4/MAO(001) heterostructures.

Scafetta, Mark↗

Deposition kinetics and characterization of stable ionomers from hexamethyldisiloxane and methacrylic acid by plasma enhanced chemical vapor deposition

A novel ionomer of hexamethyldisiloxane and methacrylic acid was synthesized by plasma enhanced chemical vapor deposition (PECVD). The PECVD process, being solventless, allows mixing of monomers with very different solubilities, and for polymers formed at high deposition rates and with high structural stability (due to the high number of cross-links and covalent bonding to the substrate) to be obtained. A kinetic study over a large set of parameters was run with the aim of determining the optimal conditions for high stability and proton conductivity of the polymer layer. Copolymers with good stability over 6 months' time in air and water were obtained, as demonstrated by ellipsometry, X-Ray reflectivity, and FT-IR spectroscopy. Stable coatings showed also proton conductivity as high as 1.1 ± 0.1 mS cm{sup −1}. Chemical analysis showed that due to the high molecular weight of the chosen precursors, it was possible to keep the plasma energy-input-per-mass low. This allowed limited precursor fragmentation and the functional groups of both monomers to be retained during the plasma polymerization.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Energy-dependent optimization of the prompt fission neutron spectrum with CGMF

Throughout the course of FY21, significant effort was put into investigating models within the LANL developed Hauser-Feshbach fission fragment decay code, CGMF, to understand and potentially solve the long-standing challenge of a too-soft prompt fission neutron spectrum, PFNS. Several inputs and models to CGMF were investigated, including the discrete nuclear levels, the optical model potential, level densities, and the fission fragment initial conditions. Some of the global models within CGMF led to a slight hardening of the neutron spectrum—particularly the likely incomplete discrete levels in through which γ-rays decay—but none of the changes where large enough for the tail of the PFNS to reproduce experimental data. A significant hardening of the spectrum tail was observed when the fission fragment initial conditions were optimized based on their sensitivities to the PFNS data for thermal incident neutrons. In this way, the parameters for the CGMF mass and total kinetic energy distributions, along with the spin cutoff factor were adjusted to better reproduce the experimental PFNS measurements. This optimization hardened the tail of the PFNS slightly but led to unphysical mass distributions for the fission fragments before neutron emission. It was clear from the above that we do not expect to be able to produce an evaluation-quality PFNS with CGMF in the near future. Challenges at thermal will persist–and possibly worsen–with increasing incident energy, where more models are needed to completely describe the fission. Basic-science research funding exceeding the amount available and scope of our NCSP funds would be needed to tackle this decade-long challenge impacting many fission-fragment event generator. And, in fact, Amy Lovell won LDRD ECR funding to do so over the next few years. Therefore, we focused in FY22 on extending evaluation capabilities beyond thermal incident neutrons, to take into account the incident energy dependence of the PFNS and fission fragment initial condition distributions in CGMF. We chose to set up the evaluation methodology to perform PFNS evaluations with CGMF across incident-neutron energies, in order to have it readily available for future NCSP evaluations when the PFNS from CGMF has improved. In this report, we outline the evaluation methodology, along with the results of the optimization, including full model calculations with CGMF using the evaluated parameters.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Kinetic simulation of magnetic field generation and collisionless shock formation in expanding laboratory plasmas

Recent laboratory experiments with laser-produced plasmas have observed and studied a number of fundamental physical processes relevant to magnetized astrophysical plasmas, including magnetic reconnection, collisionless shocks, and magnetic field generation by Weibel instability, opening up new experimental platforms for laboratory astrophysics. We develop a fully kinetic simulation model for first-principles simulation of these systems including the dynamics of magnetic fields—magnetic field generation by the Biermann battery effect or Weibel instability; advection by the ion flow, Hall effect, and Nernst effect; and destruction of the field by dissipative mechanisms. Key dimensionless parameters describing the system are derived for scaling between kinetic simulation, recent experiments, and astrophysical plasmas. First, simulations are presented which model Biermann battery magnetic field generation in plasmas expanding from a thin target. Ablation of two neighboring plumes leads to the formation of a current sheet as the opposing Biermann-generated fields collide, modeling recent laser-driven magnetic reconnection experiments. Second, we simulate recent experiments on collisionless magnetized shock generation, by expanding a piston plasma into a pre-magnetized ambient plasma. For parameters considered, the Biermann effect generates additional magnetic fields in the curved shock front and thereby increases shock particle reflection. In conclusion, both cases show the importance of kinetic processes in the interaction of plasmas with magnetic fields and open opportunities to benchmark these important processes through comparison of theory and experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Oxidation kinetics of biochar from woody and herbaceous biomass

The goal of this study was to determine oxidation kinetics for biochar produced from fast pyrolysis of various biomass feedstocks. In particular, the role of inherent ash content on the oxidation rate was evaluated. Thermogravimetric analysis (TGA) and fluidized bed combustion experiments were used to explore oxidation kinetics of six fast pyrolysis produced biochars with diverse ash content. Here, the reaction rates varied by a factor of three under chemical kinetic-limited conditions, demonstrating inorganic content impacts oxidation rate. Chemical kinetic rate coefficients were proposed as a function of compositional parameters to determine overall fit and impact. Potassium content was found to have a positive correlation, best describing the differences in the oxidation kinetic rate coefficients. Additionally, feedstocks were subjected to a 1 M HCl acid wash mitigating the catalytic activity of the metals. Acid washed biochars had lower oxidation kinetic rates compared to their unwashed counterparts, indicating the removal of catalytically active metals reduced oxidation rate. Gas composition (CO/CO 2 ) was measured during fluidized bed experiments for both acid-washed and unwashed biochars, which varied between the six biochars. Formation of CO 2 was greatly affected by catalytic metals, finding potassium content to correlate well with a higher percentage of CO 2 formation as compared to CO. Comparison of oxidation rates were made between the two experimental apparatuses to measure the effect of attrition on biochar oxidation.

09 BIOMASS FUELS↗

An ODE-Fitting Approach to Estimate Critical Aircraft Performance Parameters for Trajectory Prediction

Ground-based decision support tools (DST) in air traffic management (ATM) typically perform trajectory prediction based on aircraft performance model (APM) parameters, but some or all of these parameters might not be readily available. In particular, the three critical parameters required for trajectory prediction are the thrust setting, drag coefficients, and takeoff weight of the aircraft. Unfortunately, these parameters are coupled and appear together in physics-based kinetic models. Past approaches utilize data from a specific phase of flight (climb, level flight or descent), where one or more of the parameters are assumed to be known and estimate the remaining unknown parameters. This approach introduces bias/errors and also does not extend to scenarios where all of the above parameters are not known with sufficient accuracy. This paper is the first of its kind to propose a generalized framework for simultaneous estimation of all three critical APM parameters (thrust, drag, and mass). The proposed approach utilizes data from both the climb and descent phases and fits the ordinary differential equation for altitude in each phase using historical trajectory data available from radar tracks or ADS-B. The approach yields a set of optimized APM parameters that are best suited to fit each historical flight record. The methodology is applied on on sample flights from three different aircraft types, and the results demonstrate low fit error and consequently will yield a high level of prediction accuracy.

trajectory prediction, aircraft performance model,↗

An ODE-Fitting Approach to Estimate Critical Aircraft Performance Parameters for Trajectory Prediction

Ground-based decision support tools (DST) in air traffic management (ATM) typically perform trajectory prediction based on aircraft performance model (APM) parameters, but some or all of these parameters might not be readily available. In particular, the three critical parameters required for trajectory prediction are the thrust setting, drag coefficients, and takeoff weight of the aircraft. Unfortunately, these parameters are coupled and appear together in physics-based kinetic models. Past approaches utilize data from a specific phase of flight (climb, level flight or descent), where one or more of the parameters are assumed to be known and estimate the remaining unknown parameters. This approach introduces bias/errors and also does not extend to scenarios where all of the above parameters are not known with sufficient accuracy. This paper is the first of its kind to propose a generalized framework for simultaneous estimation of all three critical APM parameters (thrust, drag, and mass). The proposed approach utilizes data from both the climb and descent phases and fits the ordinary differential equation for altitude in each phase using historical trajectory data available from radar tracks or ADS-B. The approach yields a set of optimized APM parameters that are best suited to fit each historical flight record. The methodology is applied on on sample flights from three different aircraft types, and the results demonstrate low fit error and consequently will yield a high level of prediction accuracy.

trajectory prediction↗

Selective sodium-ion diffusion channels in Na 2 - x Fe 3 (PO 4 ) 3 positive electrode for Na-ion batteries

The electrochemical community around the world is focusing heavily on sodium-ion batteries R&D due to the huge abundance of sodium and need for large scale deployment of safe and inexpensive batteries for the electrical grid energy storage. Here, we observe the characterization of model parameters such as ionic diffusivity and interfacial kinetics, for a potential sodium-ion battery cathode material, Na 2 - x Fe 3 (PO 4 ) 3 , using experimental and computational investigations. The Na 2 - x Fe 3 (PO 4 ) 3 structure is characterized by two distinct one-dimensional (1-D) Na diffusion channels. Density Functional Theory (DFT) calculation reveals that Na is first fully removed from one of the channels followed by the subsequent removal of Na from the second channel. The experimental exploration reveals that, in the beginning of Na removal at 0 ≤ x ≤ 0.5; the sodium ion diffusivity is of the order of ~10 –11 cm 2 s –1 and then slightly decreases. This is further confirmed by comparing the calculated sodium diffusion barriers along the individual 1-D channels as well as between the channels. Nevertheless, the exchange current density slowly increases up to x=1.0 and remains quite constant thereafter. The magnitude of exchange current density is very low suggesting that the interfacial kinetics are the rate limiting factor. The obtained results indicate that Na 2 - x Fe 3 (PO 4 ) 3 could achieve better rate performance with long cycling stability through engineering of the particle morphology and microstructure.

25 ENERGY STORAGE↗