Search NASA⌕ Search

SEARCH · Search NASA

Results for “kinetic theory”

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 487 records · Page 27
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Effect of Strength of Soluble and Solid Acid Catalysts on Selectivity in Aldol Reactions

Aldol reactions can follow the classical condensation path and produce an α,β-unsaturated carbonyl compound and water, however in some cases, aldol chemistry may be steered towards a fission that gives an olefin and a carboxylic acid. This investigation examines how the strength of Brønsted acid sites—both in homogeneous and heterogeneous catalysis—controls these two pathways and their kinetics. The cross-aldol reaction between benzaldehyde and 3-pentanone served as test case. Batch reactions, conducted in toluene as solvent at a temperature of 140 °C under autogenous pressure, were analyzed by GC and in situ ATR-FTIR spectroscopy to determine product distributions and rate constants for condensation and fission pathways. A series of soluble acids including a family of sulfonic acids mostly favored the condensation pathway, with formic acid as the weakest in the series by its pKa being inactive for aldol chemistry. Significant amounts of fission products were rare except for the known selectivity of phosphoric acid, and the circumstances were difficult to delineate. For the sulfonic acid family, the logarithm of the first-order condensation rate constants scaled only roughly with pKa (water) values, whereas a good correlation was obtained with calculated deprotonation Gibbs energies in toluene. A series of H-forms of isomorphously substituted beta zeolites, HESiBEA with E = Al, Ga, Fe, or B, favored the fission pathway. Site density and site strength were characterized by calorimetric measurements of the heats of adsorption of isopropylamine, which decreased in the order Al > Ga, Fe > B. The logarithm of the site-normalized first-order fission rate constants scaled roughly with the heats of adsorption and correlated excellently with reported deprotonation energies. In conclusion, acid strength mainly affects activity and can be seen as a pre-requisite for either aldol condensation or fission chemistry, whereas additional, yet to be full clarified catalyst properties and reaction conditions are required to steer aldol chemistry towards fission selectivity.

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

Multiple Roles of Alkanethiolate‐Ligands in Direct Formation of H 2 O 2 over Pd Nanoparticles

Abstract Coadsorbed organic species including thiolates can promote direct synthesis of hydrogen peroxide from H 2 and O 2 over Pd particles. Here, density functional theory based kinetic modeling, augmented with activity measurements and vibrational spectroscopy are used to provide atomistic understanding of direct H 2 O 2 formation over alkylthiolate(RS) Pd. We find that the RS species are oxidized during reaction conditions yielding RSO 2 as the effective ligand. The RSO 2 ligand shows superior ability for proton transfer to the intermediate surface species OOH, which accelerates the formation of H 2 O 2 . The ligands promote the selectivity also by blocking sites for unselective water formation and by modifying the electronic structure of Pd. The work rationalizes observations of enhanced selectivity of direct H 2 O 2 formation over ligand‐funtionalized Pd nanoparticles and shows that engineering of organic surface modifiers can be used to promote desired hydrogen transfer routes.

Chen, Lin↗
DOE OSTI2022

Multiple Roles of Alkanethiolate‐Ligands in Direct Formation of H 2 O 2 over Pd Nanoparticles

Abstract Coadsorbed organic species including thiolates can promote direct synthesis of hydrogen peroxide from H 2 and O 2 over Pd particles. Here, density functional theory based kinetic modeling, augmented with activity measurements and vibrational spectroscopy are used to provide atomistic understanding of direct H 2 O 2 formation over alkylthiolate(RS) Pd. We find that the RS species are oxidized during reaction conditions yielding RSO 2 as the effective ligand. The RSO 2 ligand shows superior ability for proton transfer to the intermediate surface species OOH, which accelerates the formation of H 2 O 2 . The ligands promote the selectivity also by blocking sites for unselective water formation and by modifying the electronic structure of Pd. The work rationalizes observations of enhanced selectivity of direct H 2 O 2 formation over ligand‐funtionalized Pd nanoparticles and shows that engineering of organic surface modifiers can be used to promote desired hydrogen transfer routes.

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

Accurate numerical, integral methods for computing drift-kinetic Trubnikov-Rosenbluth potentials

A novel numerical method is employed to compute the integral form of the axi-symmetric Trubnikov-Rosenbluth potentials. Two methods for quadrature in pitch-angle are described and their convergence properties are studied. Careful attention is given to quadrature over a singular Green's function. Here it is shown that an infinite series representation of the Green's function can be used more efficiently than its closed form involving complete elliptic integrals. Then a collocation method in speed, with its associated quadrature scheme, is laid out and its convergence properties are studied. Using the proposed scheme, accurate low-order moments of the field collision operator are obtained using relatively few velocity space degrees of freedom. The scheme is showcased by solving for the equilibrium, axi-symmetric bootstrap current in tokamaks. A C 0 Gauss-Lobatto-Legendre finite element pitch-angle basis with vertex nodes at the trapped/passing boundary is shown, in the context of the integral methods used, to be much more efficient than the more common Legendre polynomial expansion.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Structure-preserving neural networks for the regularized entropy-based closure of a linear, kinetic, radiative transport equation

The main challenge of large-scale numerical simulation of radiation transport is the high memory and computation time requirements of discretization methods for kinetic equations. In this work, we derive and investigate a neural network-based approximation to the entropy-based closure method to accurately compute the solution of the multi-dimensional moment system with a low memory footprint and competitive computational time. We extend methods developed for the standard entropy-based closure to the regularized entropy-based closures. The main idea is to interpret structure-preserving neural network approximations of the regularized entropy-based closure as a two-stage approximation to the original entropy-based closure. We conduct a numerical analysis of this approximation and investigate optimal parameter choices. Our numerical experiments demonstrate that the method has a much lower memory footprint than traditional methods with competitive computation times and simulation accuracy. The code and all trained networks are provided on GitHub.

entropy closure↗
DOE OSTI2020

Local spin polarizations in relativistic heavy-ion collisions

Here, based on a generalized side-jump formalism for massless chiral fermions, which naturally takes into account the spin-orbit coupling in the scattering of two chiral fermions and the chiral vortical effect in a rotating chiral fermion matter, we have developed a covariant and total-angular-momentum-conserved chiral transport model to study both the global and local polarizations of this matter. For a system of massless quarks of random spin orientations and finite vorticity in a box, we have demonstrated that the model can exactly conserve the total angular momentum of the system and dynamically generate the quark spin polarization expected from a thermally equilibrated quark matter. Using this model to study the spin polarization in relativistic heavy-ion collision, we have found that the local quark spin polarizations depend strongly on the reference frame where they are evaluated as a result of the nontrivial axial charge distribution caused by the chiral vortical effect. We have further shown that because of the anomalous orbital or side-jump contribution to the quark spin polarization, the local quark polarizations calculated in the medium rest frame are qualitatively consistent with the local polarizations of Lambda hyperons measured in experiments.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Integrated modeling methodology for ash agglomeration in poly-disperse fluidized beds using particle population framework

This article discusses a unique agglomeration modeling methodology developed based on binary collisions to combine the effects of heterogeneity in ash chemistry and granular physics. A simple population balance is defined to find changes in the particle size distribution (PSD) of a fluidized bed. Thermodynamic equilibrium calculations and a computational fluid dynamics (CFD) code, are used to obtain hydrodynamic parameters. A method to calculate and use a distribution of collision frequencies in fluidized beds, corresponding to a poly-disperse particle size distribution, was developed in order to incorporate the particle level heterogeneities. A distribution of collision frequencies obtained for poly-dispersed particles showed a three orders of magnitude higher collision frequency amongst the smaller particles than the coarser ones, at the initiation of agglomeration. Ash agglomeration occurred when the slag amount (binder) was less than 10 wt% at temperatures <850 °C for Pittsburgh No. 8 coal.

42 ENGINEERING↗
DOE OSTI2021

Distinct Mechanisms and Hydricities of Cp*Ir-Based CO 2 Hydrogenation Catalysts in Basic Water

Transition metal-catalyzed reversible hydrogenation of CO 2 to formate in aqueous solutions under ambient conditions is an attractive environmentally friendly strategy for the storage and transportation of H 2 in a liquid chemical form. Here, mechanistic details for the CO 2 hydrogenation by a series of pentamethylcyclopentadienyl iridium(III) (Cp*Ir) complexes with picolinamidate ligands are investigated through density functional theory and kinetic isotope effects (KIEs) calculations. Additionally, adopting a speciation approach at a pH of 8.2, CO 2 hydrogenation mechanisms of 37 distinct protonation states and conformers of Cp*Ir type complexes are investigated. This imposes a challenging test for the computational modeling in terms of providing a consistent correlation with experimental kinetics and KIE studies of multiple catalysts instead of a single complex. Overall, H 2 heterolysis to generate an iridium hydride intermediate was demonstrated to be the rate determining step, which proceeds via distinct pathways depending on the stereoelectronic properties of the ligand. A peculiar iridium dihydride route was also uncovered, which could be optimized to accelerate the H 2 heterolysis. We have further computed the thermodynamic and kinetic hydricities of 74 complexes with [Cp*Ir(L)(H)] q ( q = -1, 0, 1 depending on the ligand charge by deprotonation or protonation events) general formula. Hydricity values do not show any noticeable correlation with the thermodynamics (Δ G ) of [Cp*Ir(L)(HCO 2 – )] formation but exhibit linear correlation with the kinetics (Δ G ‡ ) of electrophilic CO 2 attack to iridium hydride species, especially when the charges and local structural effects of the metal hydrides are considered. Insights gained in this study on (i) the factors determining the preferred pathway for the rate-limiting H 2 heterolysis step, (ii) the correlation between thermodynamic hydricity and the kinetics of CO 2 insertion to iridium hydride species, and (iii) mechanistic analysis via combination of free energy profiles and KIE studies for a series of iridium catalysts will provide guidelines for the design of next-generation catalysts for the reversible H 2 storage through CO 2 utilization.

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

Biogas sequestration to carbon nanofibers via tandem catalytic strategies

Upgrading decentralized biogas represents a sustainable route to produce valuable products while mitigating two potent greenhouse gases, namely, methane (CH 4 ) and carbon dioxide (CO 2 ). Conventional dry reforming of CH 4 with CO 2 yields syngas with low H 2 /CO ratios (≤1) and requires high temperatures (>800 °C) to overcome equilibrium constraints and abate coke deposition, which limits commercial implementation. Here we demonstrate the conversion of biogas into value-added carbon nanofibers via reaction integration in tandem reactors, while reducing the reaction temperature, shifting equilibrium limits and yielding H 2 -enriched syngas (H 2 /CO = 2–3) as a byproduct. Experimental and theoretical insights reveal that potassium (K) modification enhances carbon nanofiber formation due to synergistic effects via a balanced interplay between KO x -induced cobalt facets and cobalt carbide species. In conclusion, the energy cost and CO 2 footprint analyses highlight the potential advantages of tandem processes for the sustainable upgrading of biogas into valuable solid carbon products.

09 BIOMASS FUELS↗
DOE OSTI2020

Runaway electron transport in stochastic toroidal magnetic fields

Here, we study the transport and confinement properties of runaway electrons (RE) in the presence of magnetic fields with perturbations producing different levels of stochasticity. We use Kinetic Orbit Runaway Electron Code (KORC) [Carbajal et al. , Phys. Plasmas 24, 042512 (2017) and del-Castillo-Negrete et al. , Phys. Plasmas 25, 056104 (2018)] for simulating the full-orbit (FO) and guiding-center (GC) dynamics of RE in perturbed magnetic fields that exhibit magnetic islands. We extend previous works on this problem [Wingen et al. , Nucl. Fusion 46, 941 (2006); Izzo et al. , Nucl. Fusion 51, 063032 (2011); Papp et al. , Nucl. Fusion 51, 043004 (2011); V. Izzo and P. Parks, Phys. Plasmas 24, 060705 (2017); and Sommariva et al. , Nucl. Fusion 58, 016043 (2018)] by studying in detail full-orbit effects on the RE dynamics. We quantify FO effects on RE transport by performing one-to-one comparisons between FO and GC simulations. It is found that, for the magnetic field configurations considered, GC simulations predict twice the RE losses of FO simulations for 1 MeV and four times the RE losses of FO simulations for 25 MeV. Similarly, we show how different GC and FO dynamics of RE moving around magnetic islands can be, especially in the scenario where the RE Larmor radius is on the order of the size of the magnetic island. We also study the role of rotation of the magnetic islands on RE confinement, and we find that low-frequency toroidal rotation has no observable effect on RE transport in the cases considered. These results shed some light into the potential of avoidance or mitigation mechanisms based on magnetic perturbations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Generalized slab universal instability and its appearance in pair plasma

In this work, a generalized linear dispersion relation of electromagnetic slab universal modes is derived, taking into account arbitrary ion charge state, electron finite Larmor radius (FLR) effects, and Debye shielding effects. As a consequence, it provides more accurate predictions and is applicable to a wider range of plasmas. We find that electron FLR effects have a weakly stabilizing effect on the slab universal instability in electron–ion plasma, while Debye shielding strongly stabilizes this instability when λ D approaches ρ i (λ D is the Debye length and ρi is the ion gyroradius). In particular, we examine the stability criterion for this instability in electron–positron pair plasmas and find that the instability persists in this simplest plasma system as long as the pair plasma number density exceeds the critical value n c = B 2 /(8πm e c 2 ).

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Analytic quasi-steady evolution of marginally unstable kinetic plasma instabilities in the presence of drag and scattering

The 1D bump-on-tail problem is studied in order to determine the influence of drag on quasi-steady solutions near marginal stability (1-γ d /γ L << 1) when effective collisions are much larger than the instability growth rate (ν>>γ). In this common tokamak regime, it is rigorously shown that the paradigmatic Berk–Breizman cubic equation for the nonlinear mode evolution reduces to a much simpler differential equation, dubbed the time-local cubic equation, which can be solved directly. It is found that in addition to increasing the saturation amplitude, drag introduces a shift in the apparent oscillation frequency by modulating the saturated wave envelope. Furthermore, excellent agreement is found between the analytic solution for the mode evolution and both the numerically integrated Berk–Breizman cubic equation and fully nonlinear 1D Vlasov simulations. Experimentally isolating the contribution of drag to the saturated mode amplitude for verification purposes is explored but complicated by the reality that the amount of drag cannot be varied independently of other key parameters in realistic scenarios. While the influence of drag is modest when the ratio of effective drag to effective scattering α/ν is very small, it can become substantial when α/ν ≳ 0.5, suggesting that drag should be accounted for in quantitative models of fast-ion-driven instabilities in fusion plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Computing the Paschen curve for argon with speed-limited particle-in-cell simulation

Upon inclusion of collisions, the speed-limited particle-in-cell (SLPIC) simulation method successfully computed the Paschen curve for argon. The 1D3V simulations modeled an electron cascade across an argon-filled capacitor, including electron-neutral ionization, electron-neutral elastic collisions, electron-neutral excitation, and ion-induced secondary electron emission. In electrical breakdown, the timescale difference between ion and electron motion makes traditional PIC methods computationally slow. To decrease this timescale difference and speed up computation, here we used SLPIC, a time-domain algorithm that limits the speed of the fastest electrons in the simulation. The SLPIC algorithm facilitates a straightforward, fully kinetic treatment of dynamics and collisions. SLPIC was as accurate as PIC, but ran up to 200 times faster. SLPIC accurately computed the Paschen curve for argon over three orders of magnitude in pressure.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

dc electrical conductivity in strongly magnetized plasmas

A generalized Ohm's law is derived to treat strongly magnetized plasmas in which the electron gyrofrequency significantly exceeds the electron plasma frequency. Here, strong magnetization of electrons causes the frictional drag between electrons and ions due to Coulomb collisions to shift, producing an additional transverse resistivity term in the generalized Ohm's law that is perpendicular to both the current ($J$) and the Hall ($JxB$) direction. In the limit of very strong magnetization the parallel resistivity is found to increase by a factor of 3/2, and the perpendicular resistivity by a factor of $\frac{3}{4} \ln (m_i/m_e)$, where $m_i$ and $m_e$ are the ion and electron masses. These results suggest that strong magnetization significantly changes the magnetohydrodynamic evolution of a plasma.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Kinetic Landau-fluid closures of non-Maxwellian distributions

New kinetic Landau-fluid closures, based on the cutoff Maxwellian distribution, are derived. A special static case is considered (the frequency ω=0). In the strongly collisional regime, our model reduces to Braginskii's heat flux model, and the transport is local. In the weak collisional regime, our model indicates that the heat flux is non-local and recovers the Hammett–Perkins model while the value of the cutoff velocity approaches to infinity. We compare the thermal transport coefficient χ of Maxwellian, cutoff Maxwellian and super-Gaussian distribution. The results show that the reduction of the high-speed tail particles leads to the corresponding reduction of the thermal transport coefficient χ across the entire range of collisionality, more reduction of the free streaming transport toward the weak collisional regime. In the collisionless limit, χ approaches to zero for the cutoff Maxwellian and the super-Gaussian distribution but remains finite for Maxwellian distribution. χ is complex if the cutoff Maxwellian distribution is asymmetric. The Im(χ) approaches to different convergent values in both collisionless and strongly collisional limit, respectively. It yields an additional streaming heat flux in comparison with the symmetric cutoff Maxwellian distribution. Furthermore, due to the asymmetric distribution, there is a background heat flux q0 though there is no perturbation. Furthermore, the derived Landau-fluid closures are general for fluid moment models, and applicable for the cutoff Maxwellian distribution in an open magnetic field line region, such as the scape-off-layer of Tokamak plasmas, in the thermal quench plasmas during a tokamak disruption, and the super-Gaussian electron distribution function due to inverse bremsstrahlung heating in laser-plasma studies.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Symbolic regression development of empirical equations for diffusion in Lennard-Jones fluids

Symbolic regression (SR) with a multi-gene genetic program has been used to elucidate new empirical equations describing diffusion in Lennard-Jones (LJ) fluids. Some examples include equations to predict self-diffusion in pure LJ fluids and equations describing the finite-size correction for self-diffusion in binary LJ fluids. The performance of the SR-obtained equations was compared to that of both the existing empirical equations in the literature and to the results from artificial neural net (ANN) models recently reported. It is found that the SR equations have improved predictive performance in comparison to the existing empirical equations, even though employing a smaller number of adjustable parameters, but show an overall reduced performance in comparison to more extensive ANNs.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗