Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “KINETICS”

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 37 records · Page 2

Energy dynamics and current sheet structure in fluid and kinetic simulations of decaying magnetohydrodynamic turbulence

Simulations of decaying magnetohydrodynamic (MHD) turbulence are performed with a fluid and a kinetic code. The initial condition is an ensemble of long-wavelength, counter-propagating, shear-Alfvén waves, which interact and rapidly generate strong MHD turbulence. The total energy is conserved and the rate of turbulent energy decay is very similar in both codes, although the fluid code has numerical dissipation, whereas the kinetic code has kinetic dissipation. The inertial range power spectrum index is similar in both the codes. The fluid code shows a perpendicular wavenumber spectral slope of k{sub ⊥}{sup −1.3}. The kinetic code shows a spectral slope of k{sub ⊥}{sup −1.5} for smaller simulation domain, and k{sub ⊥}{sup −1.3} for larger domain. We estimate that collisionless damping mechanisms in the kinetic code can account for the dissipation of the observed nonlinear energy cascade. Current sheets are geometrically characterized. Their lengths and widths are in good agreement between the two codes. The length scales linearly with the driving scale of the turbulence. In the fluid code, their thickness is determined by the grid resolution as there is no explicit diffusivity. In the kinetic code, their thickness is very close to the skin-depth, irrespective of the grid resolution. This work shows that kinetic codes can reproduce the MHD inertial range dynamics at large scales, while at the same time capturing important kinetic physics at small scales.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Asymmetric “melting” and “freezing” kinetics of the magnetostructural phase transition in B2-ordered FeRh epilayers

Synchrotron X-ray diffraction was used to study the phase transformation processes during the magnetostructural transition in a B2-ordered FeRh (001)-oriented epilayer grown on MgO by sputtering. Out-of-plane lattice constant measurements within the hysteretic regime of the transition reveal a microstructure consistent with the coexistence of lattice-expanded and contracted phases in spatially distinct regions. It was found that the phase separation is more pronounced during cooling than heating. Furthermore, whilst lattice-expanded domains that span the height of the film can be undercooled by several kelvins, there is no equivalent superheating. This asymmetry between the cooling and heating processes in FeRh is consistent with the difference in the kinetics of generic freezing and melting transitions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Coupled Kinetics of Ferrihydrite Transformation and As(V) Sequestration under the Effect of Humic Acids: A Mechanistic and Quantitative Study

In natural environments, kinetics of As(V) sequestration/release is usually coupled with the dynamic iron mineral transformation, which is further influenced by the presence of the natural organic matter (NOM). However, both quantitative and mechanistic understanding of this process is still lacking. In this study, we investigated the effect of humic acids (HA), a major component of NOM, on the coupled kinetics of ferrihydrite transformation into hematite/goethite and sequestration of As(V) in iron minerals. Time-resolved As(V) and HA interactions with iron minerals were studied in stirred-flow kinetic experiments, using aberration-corrected scanning transmission electron microscopy (STEM), and X-ray adsorption spectroscopy. Based on the experimental results, we developed a mechanistic kinetics model for As(V) reactions with iron minerals with varying iron mineral compositions and solution chemistry. Our results showed that HA slowed ferrihydrite transformation and the rates of As(V) speciation changes within iron minerals, while HA did not incorporate into iron minerals. Our kinetics model is able to account for the variations of iron mineral compositions, solution chemistry, and As(V) speciation, which has significant environmental implications on predicting the As(V) behavior in the environment.

ferrihydrite, humic acid↗

Nonlinear theory of drift-cyclotron kinetics and the possible breakdown of gyro-kinetics

A nonlinear theory of drift-cyclotron kinetics (termed cyclo-kinetics here) is formulated to test the breakdown of the gyro-kinetic approximations. Six dimensional cyclo-kinetics can be regarded as an extension of five dimensional gyro-kinetics to include high-frequency cyclotron waves, which can interrupt the low-frequency gyro-averaging in the (sixth velocity grid) gyro-phase angle. Nonlinear cyclo-kinetics has no limit on the amplitude of the perturbations. Formally, there is no gyro-averaging when all cyclotron (gyro-phase angle) harmonics of the perturbed distribution function (delta-f) are retained. Retaining only the (low frequency) zeroth cyclotron harmonic in cyclo-kinetics recovers both linear and nonlinear gyro-kinetics. Simple recipes are given for converting continuum nonlinear delta-f gyro-kinetic transport simulation codes to cyclo-kinetics codes by retaining (at least some) higher cyclotron harmonics.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Assessment of ion kinetic effects in shock-driven inertial confinement fusion implosions using fusion burn imaging

The significance and nature of ion kinetic effects in D{sup 3}He-filled, shock-driven inertial confinement fusion implosions are assessed through measurements of fusion burn profiles. Over this series of experiments, the ratio of ion-ion mean free path to minimum shell radius (the Knudsen number, N{sub K}) was varied from 0.3 to 9 in order to probe hydrodynamic-like to strongly kinetic plasma conditions; as the Knudsen number increased, hydrodynamic models increasingly failed to match measured yields, while an empirically-tuned, first-step model of ion kinetic effects better captured the observed yield trends [Rosenberg et al., Phys. Rev. Lett. 112, 185001 (2014)]. Here, spatially resolved measurements of the fusion burn are used to examine kinetic ion transport effects in greater detail, adding an additional dimension of understanding that goes beyond zero-dimensional integrated quantities to one-dimensional profiles. In agreement with the previous findings, a comparison of measured and simulated burn profiles shows that models including ion transport effects are able to better match the experimental results. In implosions characterized by large Knudsen numbers (N{sub K} ∼ 3), the fusion burn profiles predicted by hydrodynamics simulations that exclude ion mean free path effects are peaked far from the origin, in stark disagreement with the experimentally observed profiles, which are centrally peaked. In contrast, a hydrodynamics simulation that includes a model of ion diffusion is able to qualitatively match the measured profile shapes. Therefore, ion diffusion or diffusion-like processes are identified as a plausible explanation of the observed trends, though further refinement of the models is needed for a more complete and quantitative understanding of ion kinetic effects.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Rate theory of solvent exchange and kinetics of Li{sup +} − BF{sub 4}{sup −}/PF{sub 6}{sup −} ion pairs in acetonitrile

In this paper, we describe our efforts to apply rate theories in studies of solvent exchange around Li{sup +} and the kinetics of ion pairings in lithium-ion batteries (LIBs). We report one of the first computer simulations of the exchange dynamics around solvated Li{sup +} in acetonitrile (ACN), which is a common solvent used in LIBs. We also provide details of the ion-pairing kinetics of Li{sup +}-[BF{sub 4}] and Li{sup +}-[PF{sub 6}] in ACN. Using our polarizable force-field models and employing classical rate theories of chemical reactions, we examine the ACN exchange process between the first and second solvation shells around Li{sup +}. We calculate exchange rates using transition state theory and weighted them with the transmission coefficients determined by the reactive flux, Impey, Madden, and McDonald approaches, and Grote-Hynes theory. We found the relaxation times changed from 180 ps to 4600 ps and from 30 ps to 280 ps for Li{sup +}-[BF{sub 4}] and Li{sup +}-[PF{sub 6}] ion pairs, respectively. These results confirm that the solvent response to the kinetics of ion pairing is significant. Our results also show that, in addition to affecting the free energy of solvation into ACN, the anion type also should significantly influence the kinetics of ion pairing. These results will increase our understanding of the thermodynamic and kinetic properties of LIB systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Growth Kinetics and Size Distribution Dynamics of Viscous Secondary Organic Aerosol

Low bulk diffusivity inside viscous semisolid atmospheric secondary organic aerosol (SOA) can prolong equilibration time scale, but its broader impacts on aerosol growth and size distribution dynamics are poorly understood. In this article, we present quantitative insights into the effects of bulk diffusivity on the growth and evaporation kinetics of SOA formed under dry conditions from photooxidation of isoprene in the presence of a bimodal aerosol consisting of Aitken (ammonium sulfate) and accumulation (isoprene or α-pinene SOA) mode particles. Aerosol composition measurements and evaporation kinetics indicate that isoprene SOA is composed of several semivolatile organic compounds (SVOCs), with some reversibly reacting to form oligomers. Model analysis shows that liquid-like bulk diffusivities can be used to fit the observed evaporation kinetics of accumulation mode particles but fail to explain the growth kinetics of bimodal aerosol by significantly under-predicting the evolution of the Aitken mode. In contrast, the semisolid scenario successfully reproduces both evaporation and growth kinetics, with the interpretation that hindered partitioning of SVOCs into large viscous particles effectively promotes the growth of smaller particles that have shorter diffusion time scales. This effect has important implications for the growth of atmospheric ultrafine particles to climatically active sizes.

54 ENVIRONMENTAL SCIENCES↗

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↗

Multiple highly expressed phospho enol pyruvate carboxylase genes have divergent enzyme kinetic properties in two C4 grasses

Abstract Background and Aims Phosphoenolpyruvate (PEP) carboxylase (PEPC) catalyses the irreversible carboxylation of PEP with bicarbonate to produce oxaloacetate. This reaction powers the carbon-concentrating mechanism (CCM) in plants that perform C4 photosynthesis. This CCM is generally driven by a single PEPC gene product that is highly expressed in the cytosol of mesophyll cells. We found two C4 grasses, Panicum miliaceum and Echinochloa colona, that each have two highly expressed PEPC genes. We characterized the kinetic properties of the two most abundant PEPCs in E. colona and P. miliaceum to better understand how the enzyme’s amino acid structure influences its function. Methods Coding sequences of the two most abundant PEPC proteins in E. colona and P. miliaceum were synthesized by GenScript and were inserted into bacteria expression plasmids. Point mutations resulting in substitutions at conserved amino acid residues (e.g. N-terminal serine and residue 890) were created via site-directed PCR mutagenesis. The kinetic properties of semi-purified plant PEPCs from Escherichia coli were analysed using membrane-inlet mass spectrometry and a spectrophotometric enzyme-coupled reaction. Key Results The two most abundant P. miliaceum PEPCs (PmPPC1 and PmPPC2) have similar sequence identities (>95 %), and as a result had similar kinetic properties. The two most abundant E. colona PEPCs (EcPPC1 and EcPPC2) had identities of ~78 % and had significantly different kinetic properties. The PmPPCs and EcPPCs had different responses to allosteric inhibitors and activators, and substitutions at the conserved N-terminal serine and residue 890 resulted in significantly altered responses to allosteric regulators. Conclusions The two, significantly expressed C4Ppc genes in P. miliaceum were probably the result of genomes combining from two closely related C4Panicum species. We found natural variation in PEPC’s sensitivity to allosteric inhibition that seems to bypass the conserved 890 residue, suggesting alternative evolutionary pathways for increased malate tolerance and other kinetic properties.

DiMario, Robert J. (ORCID:0000000250566868)↗

Rigorous Kinetic Modeling and Optimization Study of a Modified Claus Unit for an Integrated Gasification Combined Cycle (IGCC) Power Plant with CO{sub 2} Capture

The modified Claus process is one of the most common technologies for sulfur recovery from acid gas streams. Important design criteria for the Claus unit, when part of an Integrated Gasification Combined Cycle (IGCC) power plant, are the ability to destroy ammonia completely and the ability to recover sulfur thoroughly from a relatively low purity acid gas stream without sacrificing flame stability. Because of these criteria, modifications to the conventional process are often required, resulting in a modified Claus process. For the studies discussed here, these modifications include the use of a 95% pure oxygen stream as the oxidant, a split flow configuration, and the preheating of the feeds with the intermediate pressure steam generated in the waste heat boiler (WHB). In the future, for IGCC plants with CO{sub 2} capture, the Claus unit must satisfy emission standards without sacrificing the plant efficiency in the face of typical disturbances of an IGCC plant, such as rapid change in the feed flow rates due to load-following and wide changes in the feed composition because of changes in the coal feed to the gasifier. The Claus unit should be adequately designed and efficiently operated to satisfy these objectives. Even though the Claus process has been commercialized for decades, most papers concerned with the modeling of the Claus process treat the key reactions as equilibrium reactions. Such models are validated by manipulating the temperature approach to equilibrium for a set of steady-state operating data, but they are of limited use for dynamic studies. One of the objectives of this study is to develop a model that can be used for dynamic studies. In a Claus process, especially in the furnace and the WHB, many reactions may take place. In this work, a set of linearly independent reactions has been identified, and kinetic models of the furnace flame and anoxic zones, WHB, and catalytic reactors have been developed. To facilitate the modeling of the Claus furnace, a four-stage method was devised so as to determine which set of linearly independent reactions would best describe the product distributions from available plant data. Various approaches are taken to derive the kinetic rate expressions, which are either missing in the open literature or found to be inconsistent. A set of plant data is used for optimal estimation of the kinetic parameters. The final model agrees well with the published plant data. Using the developed kinetics models of the Claus reaction furnace, WHB, and catalytic stages, two optimization studies are carried out. The first study shows that there exists an optimal steam pressure generated in the WHB that balances hydrogen yield, oxygen demand, and power generation. In the second study, it is shown that an optimal H{sub 2}S/SO{sub 2} ratio exists that balances single-pass conversion, hydrogen yield, oxygen demand, and power generation. In addition, an operability study has been carried out to examine the operating envelope in which both the H{sub 2}S/SO{sub 2} ratio and the adiabatic flame temperature can be controlled in the face of disturbances typical for the operation of an IGCC power plant with CO{sub 2} capture. Impact of CO{sub 2} capture on the Claus process has also been discussed.

20 FOSSIL-FUELED POWER PLANTS↗

Covalent bonds are created by the drive of electron waves to lower their kinetic energy through expansion

An analysis based on the variation principle shows that in the molecules H{sub 2}{sup +}, H{sub 2}, B{sub 2}, C{sub 2}, N{sub 2}, O{sub 2}, F{sub 2}, covalent bonding is driven by the attenuation of the kinetic energy that results from the delocalization of the electronic wave function. For molecular geometries around the equilibrium distance, two features of the wave function contribute to this delocalization: (i) Superposition of atomic orbitals extends the electronic wave function from one atom to two or more atoms; (ii) intra-atomic contraction of the atomic orbitals further increases the inter-atomic delocalization. The inter-atomic kinetic energy lowering that (perhaps counter-intuitively) is a consequence of the intra-atomic contractions drives these contractions (which per se would increase the energy). Since the contractions necessarily encompass both, the intra-atomic kinetic and potential energy changes (which add to a positive total), the fact that the intra-atomic potential energy change renders the total potential binding energy negative does not alter the fact that it is the kinetic delocalization energy that drives the bond formation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Explicit integration with GPU acceleration for large kinetic networks

We demonstrate the first implementation of recently-developed fast explicit kinetic integration algorithms on modern graphics processing unit (GPU) accelerators. Taking as a generic test case a Type Ia supernova explosion with an extremely stiff thermonuclear network having 150 isotopic species and 1604 reactions coupled to hydrodynamics using operator splitting, we demonstrate the capability to solve of order 100 realistic kinetic networks in parallel in the same time that standard implicit methods can solve a single such network on a CPU. This orders-of-magnitude decrease in computation time for solving systems of realistic kinetic networks implies that important coupled, multiphysics problems in various scientific and technical fields that were intractable, or could be simulated only with highly schematic kinetic networks, are now computationally feasible.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Hybrid simulations of magnetic reconnection with kinetic ions and fluid electron pressure anisotropy

We present the first hybrid simulations with kinetic ions and recently developed equations of state for the electron fluid appropriate for reconnection with a guide field. The equations of state account for the main anisotropy of the electron pressure tensor.Magnetic reconnection is studied in two systems, an initially force-free current sheet and a Harris sheet. The hybrid model with the equations of state is compared to two other models, hybrid simulations with isothermal electrons and fully kinetic simulations. Including the anisotropicequations of state in the hybrid model provides a better match to the fully kinetic model. In agreement with fully kinetic results, the main feature captured is the formation of an electron current sheet that extends several ion inertial lengths. This electron current sheet modifies the Hall magnetic field structure near the X-line, and it is not observed in the standard hybrid model with isotropic electrons. The saturated reconnection rate in this regime nevertheless remains similar in all three models. Here, implications for global modeling are discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗