Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “kinetic distributions”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

A moment-conserving discontinuous Galerkin representation of the relativistic Maxwellian distribution

Kinetic simulations of relativistic gases and plasmas are critical for understanding diverse astrophysical and terrestrial systems, but the accurate construction of the relativistic Maxwellian, the Maxwell–Jüttner distribution, on a discrete simulation grid is challenging. Difficulties arise from the finite velocity bounds of the domain, which may not capture the entire distribution function, as well as errors introduced by projecting the function onto a discrete grid. Here, we present a novel scheme for iteratively correcting the moments of the projected distribution applicable to all grid-based discretizations of the relativistic kinetic equation. In addition, we describe how to compute the needed nonlinear quantities, such as Lorentz boost factors, in a discontinuous Galerkin scheme through a combination of numerical quadrature and weak operations. The resulting method accurately captures the distribution function and ensures that the moments match the desired values to machine precision.

astrophysical plasmas↗
DOE OSTI2025

On High-Order/Low-Order and Micro-Macro Methods for Implicit Time-Stepping of the BGK Model

In this paper, a high-order/low-order (HOLO) method is combined with a micro-macro (MM) decomposition to accelerate iterative solvers in fully implicit time-stepping of the Bhatnagar–Gross–Krook (BGK) equation for gas dynamics. The MM formulation represents a kinetic distribution as the sum of a local Maxwellian and a perturbation. In highly collisional regimes, the perturbation away from initial and boundary layers is small and can be compressed to reduce the overall storage cost of the distribution. The convergence behavior of the MM methods, the usual HOLO method, and the standard source iteration method is analyzed on a linear BGK model. Both the HOLO and MM methods are implemented using a discontinuous Galerkin (DG) discretization in phase space, which naturally preserves the consistency between high- and low-order models required by the HOLO approach. Furthermore, the accuracy and performance of these methods are compared on the Sod shock tube problem and a sudden wall heating boundary layer problem. Overall, the results demonstrate the robustness of the MM and HOLO approaches and illustrate the compression benefits enabled by the MM formulation when the kinetic distribution is near equilibrium.

BGK model↗
DOE OSTI2025

Verification of the PERSENT Software

Ongoing commercial design activities require a thorough verification of the Argonne Reactor Computation codes be performed. DIF3D is central to this system and substantial work has been done to verify its accuracy on several identified commercial needs. This manuscript details the verification work done on PERSENT which relies upon the DIF3D code for its forward and adjoint flux solution. Previous work identified the PERSENT features required to be verified to support commercial design activities, features of which are generally applicable to hexagonal-Z fast reactor designs. The scope of this verification effort includes verifying PERSENT’s ability to correctly calculate four key quantities: perturbation worth distributions, kinetics parameters, sensitivity coefficients, and cross section uncertainty quantification. This manuscript provides the verification tasks and their results with respect to these quantities needed for commercial design activities. For the perturbation worth distributions, hand calculations are deployed to verify the PERSENT calculated results. Similarly, hand calculation of the PERSENT computed kinetics parameters is also used to verify the PERSENT results. In both of these, the input to PERSENT is manipulated to ensure the hand calculation exactly matches the equations PERSENT is calculating. The sensitivity coefficients involve calculating the derivatives of a parameter (such as reactivity worth), with respect to the cross section data. Direct finite difference calculations with DIF3D are used to verify the PERSENT calculated results. For the uncertainty quantification, manufactured input to PERSENT is used to allow an exact hand calculation to reproduce the PERSENT calculated results. The work detailed in this report verified that significant issues were identified for earlier versions of PERSENT for sensitivity coefficients which were corrected in this work and thus version 12.1.0 of PERSENT must be used to reproduce all of the verified work in this report.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2026

Stochastic Modeling of the Joint Neutron Number-Cumulative Fission Fragment Kinetic Energy Deposition Distribution and its Statistical Moments [Slides]

We investigate the joint distribution of the neutron number and cumulative fission-fragment kinetic energy (FKE) deposition, with a specific focus on low-order statistical moments: the mean, variance, and correlation. Starting from a point-kinetic framework, we derive a forward Master equation (FME) for the joint distribution and develop the corresponding moment equations.

42 ENGINEERING↗
DOE OSTI2024

Convective and Orographic Origins of the Mesoscale Kinetic Energy Spectrum

The mesoscale spectrum describes the distribution of kinetic energy in the Earth's atmosphere between length scales of 10 and 400 km. Since the first observations, the origins of this spectrum have been controversial. At synoptic scales, the spectrum follows a –3 spectral slope, consistent with two-dimensional turbulence theory, but a shallower –5/3 slope was observed at the shorter mesoscales. The cause of the shallower slope remains obscure, illustrating our lack of understanding. Through a novel coarse-graining methodology, we are able to present a spatio-temporal climatology of the spectral slope. We find convection and orography have a shallowing effect and can quantify this using “conditioned spectra.” These are typical spectra for a meteorological condition, obtained by aggregating spectra where the condition holds. This allows the investigation of new relationships, such as that between energy flux and spectral slope. Potential future applications of our methodology include predictability research and model validation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Uniform secondary current distribution at disk electrodes under Tafel kinetics enabled by concentric current-shielding rings

Non-uniform secondary current distribution at rotating disk electrodes (RDE) is a common problem when using resistive electrolyte media or large applied currents. In a recent publication, we have shown that auxiliary electrodes such as the ring of a rotating ring disk electrode (RRDE) can help suppress current non-uniformities at the disk enabling reliable electroanalytical measurements. However, this previous work considered linear kinetics where current distribution non-uniformities were moderate. In the present contribution, we consider current distribution non-uniformities encountered under Tafel kinetics. We show, for the case of 2,5-dihydroxy-1,4-benzoquinone (DHBQ) reduction, that optimally chosen ring conditions serve to provide effective shielding at the disk edge rendering the overall disk current distribution to be uniform. Numerical modeling and scaling analysis (using the Wagner number) are presented to aid a user in determining the optimal ring current density for achieving uniform disk current distribution under Tafel kinetics. This approach is especially useful when studying soluble-soluble redox transitions for which, unlike deposit distribution in electrodeposition, the current distribution non-uniformity is not visually apparent.

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

Relativistic approach to manipulating angular distribution of charged particles via kinetic equations

Deflection angles of charged particles interacting with materials play a critical role in various plasma applications. The development of a mathematically well-posed kinetic collision operator that accounts for deflection angles of strong Coulomb interactions remains a fundamental open problem. This paper presents a relativistic method for modifying the electromagnetic field in an anisotropic and adjustable manner to manipulate a system of charged particles, specifically by the transfer of angular momentum from a superluminal wave source to particles at specific times and locations. The method provides a mechanism to influence the scattering outcomes of strong interactions by manipulating the angular distribution of particles, and thus the deflection angles of their interactions with a solid surface, without requiring detailed knowledge of the kinetic collision operator. To this end, we demonstrate how a specific type of singularity, generated by Maxwell's equations for a superluminal wave source at the boundary of the plasma, can modify the electromagnetic field in a highly directional manner. The proposed method can lead to the development of novel approaches for controlling interactions of charged particles with a material in plasma systems. Published by the American Physical Society 2025

Moini, Nima (ORCID:0009000929568824)↗
DOE OSTI2025

Low-energy electron driven reactions in 2-bromo-5-nitrothiazole

Thiazole derivatives are biologically relevant molecules, used also in pharmaceutical applications. Herein, we report results for electron attachment to 2-bromo-5-nitrothiazole (BNT) in the gas phase. Employing two crossed electron–molecule beam experiments, we determined the efficiency curves of various fragment anions as a function of the initial electron energy between about 0 and 10 eV as well as the emission angle and kinetic energy distributions of Br− and NO2− ions formed from a resonance near 4 eV. The experiments were supported by quantum chemical calculations, exploring possible dissociation pathways along with their reaction energies. We also compared the electron attachment characteristics of BNT with those of the native thiazole molecule by performing electron attachment experiments and calculations for this molecule as well. Compared to thiazole, which is primarily degraded only by electrons with kinetic energies between about 5 and 10 eV, BNT is susceptible to low-energy electrons near 0 eV with enhanced cross section for (dissociative) electron attachment. However, although BNT offers two localization sites with high electron affinity (Br and NO2 moieties), we do not find the corresponding anions as the dominant negatively charged species formed upon electron attachment. Instead, the reaction channels with an abstraction of Br and NO2 as neutral radicals prevail, accompanied by the opening of the thiazole ring due to the relatively weak C–S bond.

Chen, Jiakuan (ORCID:0009000729608641)↗
DOE OSTI2025

Tropical Cyclone Super Resolution using conditional diffusion denoising probabilistic model from mesoscale simulation to LES

Accurate modeling of tropical cyclone wind fields is essential for the design, risk assessment, and operational planning of offshore energy infrastructure. While mesoscale simulations are widely used thanks to their computational efficiency, they lack the necessary resolution to capture key features such as wind shear and veer profiles as well as the distribution turbulent kinetic energy (TKE). High-fidelity large-eddy simulation (LES) models on the other hand, can resolve turbulent structures and provide a more accurate representation of the complex wind field, albeit at a higher computational cost. To address this modeling gap, we introduce a two-part generative framework to enhance the resolution and physics-capturing ability of mesoscale simulations. First, a reduced-order model based on Karhunen–Loève (KL) decomposition is used to extract dominant spatial modes from one-dimensional mean wind profiles. A multilayer perceptron (MLP) is trained to map mesoscale mode weights to their LES counterparts, enabling accurate reconstruction of vertical velocity profiles. Second, a conditional Diffusion Denoising Probabilistic Model (DDPM) is developed to super-resolve coarse and low-fidelity mesoscale velocity fields, recovering fine-scale turbulence structures and stress distributions. The framework is evaluated across different tropical cyclone intensity categories defined by the Saffir–Simpson scale and demonstrates strong performance in both interpolation and extrapolation tasks. The generated fields accurately reproduce spatial coherence, stress distributions, and spectral energy characteristics observed in LES data. By bridging the fidelity gap between mesoscale and LES outputs, this approach offers a scalable, data-driven solution for enhancing the representation of tropical cyclone wind fields, enabling more robust offshore energy infrastructure systems design in tropical-cyclone-prone areas.

17 WIND ENERGY↗
DOE OSTI2025

Approximating the particle distribution in rotating and tandem mirror traps

Steady-state distribution functions can be used to calculate stability conditions for modes, radiation energy losses and particle loss rates. Heuristic analytic approximations to these distributions can capture key behaviors of the true distributions such as the relative speeds of different transport processes while possessing computational advantages over their numerical counterparts. In this paper, we motivate and present a closed-form analytic model for a distribution of particles in a centrifugal or tandem mirror. We find that our model outperforms other known models in approximating numerical steady-state simulations outside of a narrow range of low confining potentials. We demonstrate the model’s suitability in the high confining potential regime for applications such as loss-cone stability thresholds, fusion yields and available energy.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Elucidating the Mechanism of Helium Evaporation from Liquid Water

Here, we investigate the evaporation of trace amounts of helium solvated in liquid water using molecular dynamics simulations and theory. Consistent with experimental observations, we find a super-Maxwellian distribution of kinetic energies of evaporated helium. This excess of kinetic energy over typical thermal expectations is explained by an effective continuum theory of evaporation based on a Fokker–Planck equation, parametrized molecularly by a potential of mean force and position-dependent friction. Using this description, we find that helium evaporation is strongly influenced by the friction near the interface, which is anomalously small near the Gibbs dividing surface due to the ability of the liquid–vapor interface to deform around the gas particle. Our reduced description provides a mechanistic interpretation of trace gas evaporation as the motion of an underdamped particle in a potential subject to a viscous environment that varies rapidly across the air–water interface. From it we predict the temperature dependence of the excess kinetic energy of evaporation, which is yet to be measured.

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

Toward Efficient Entropic Recycling by Mastering Ring–Chain Kinetics

Traditional chemical recycling approaches for condensation polymers suffer compounding energy losses and CO 2 emissions across multiple polymerization and depolymerization cycles. Entropic recycling can address these energy losses by entrapping free energy within the deconstruction products. Entropic recycling involves depolymerization to macrocyclic monomers, but such processes have not been feasible due to the high dilutions typically required to generate macrocyclic compounds. Here, we leverage selective catalysis to allow entropic recycling at concentrations 20–2000× higher than typical for macrocyclization reactions. We find that Ru-based olefin metathesis catalysts containing bulky iodine ligands significantly bias the ring–chain kinetic product distribution during ring-closing metathesis (RCM) toward the formation of oligomeric cycloalkenes. Further improvements in reaction concentration and macrocycle yield are obtained by using high catalyst loadings and by predisposing the alkene substrates to undergo favorable macrocyclization. These RCM optimizations translate effectively to cyclodepolymerization (CDP) of an olefin-containing polymer, with RCM and CDP affording similar macrocycle product distributions under identical reaction conditions. Macrocycle polymerization by entropy-driven ring-opening metathesis provides much higher molecular weight polymers than condensation polymerization of linear analogues, reducing the time to achieve high molecular weight from hours to minutes and enabling polymerization at room temperature. Finally, our findings re-emphasize the importance of energy consumption during a polymer’s lifecycle and provide a framework for the design of efficient entropic recycling systems.

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

Michel Electron Selection with SPINE for DUNE Far Detector Simulation

Michel electrons are a valuable input for particle detector calibration due to their consistent kinetic energy distribution. This report details the evaluation of a Michel electron identification method's application to simulated data from the DUNE (Deep Underground Neutrino Experiment) far detector. This method, which relies on the neural network-based particle classification software SPINE (Scalable Particle Imaging with Neural Embeddings), was developed and calibrated using simulated data for the SBND (Short-Baseline Neutrino Detector) experiment before being applied to simulated DUNE data from a 1x2x6 subset of far detector modules.

Wilson, Dante [Colorado State U.]↗
DOE OSTI2025

Coupling of Lipid Peroxidation and Criegee Intermediate Mediated Autoxidation in the Heterogeneous Oxidation of Linoleic Acid Aerosols

Autoxidation is an established mechanism for the degradation of organic molecules, which is relevant in the atmosphere, combustion processes, the environment, and the rancidification of lipids (commonly known as lipid peroxidation). Autoxidation proceeds via radical chain reactions involving hydroxyl (•OH), peroxy (RO2•), and alkoxy radicals, which are also prominent oxidants in the atmosphere. Recent reports have provided evidence for an alternative autoxidation mechanism driven instead by Criegee intermediates (CIs), which are produced from the reaction of β-hydroxy peroxy radicals (β-OH-RO2•). This work evaluates the contributions of these two mechanisms in the •OH initiated heterogeneous oxidation of linoleic acid (LA) aerosols. Reaction kinetics and product distributions are monitored using a vacuum ultraviolet photoionization aerosol time-of-flight mass spectrometer. To explain the observed kinetics, a kinetic model is developed that incorporates both the conventional peroxidation and alternative CI-mediated autoxidation mechanisms. We observe that the CI-mediated autoxidation pathways enhance the heterogeneous autoxidation rate, while the peroxidation reactions, although present, contributes less to the overall oxidation rate. α-Acyloxyalkyl hydroperoxides (AAHPs) are identified as key indicators for bimolecular reactions of CI with LA, highlighting the role of LA as a CI scavenger. Moreover, the measured functionalized LA products with hydroxyl or carbonyl group(s), serve as markers for the peroxidation reactions. In summary, this work presents a quantitative framework to understand the coupled reaction network of •OH, RO2•, β-OH-RO2• radicals, and CI in driving heterogeneous autoxidation, which is crucial for understanding degradation mechanisms of organic molecules in the environment and atmosphere.

Criegee intermediate↗
DOE OSTI2025

Kinetics Measurements in Resistive Electrolytes Using Ring-Disk Electrode: Ring as Current “Shield” Enables Uniform Disk Current Distribution

Rotating disk electrodes are commonly used for electrochemical kinetics measurements. A major disadvantage of these types of electrodes is their nonuniform secondary current distribution, especially when performing electroanalytical measurements in resistive electrolytes. Such nonuniform current distribution can render the values of kinetics constants extracted from the disk electrode to be highly inaccurate. Furthermore, one emerging class of electrolytes that suffer from low ionic conductivities is deep eutectic solvents (DES). DES are a promising class of electrolytes for various emerging applications; however, due to their resistive nature, the secondary current distribution when using them is typically highly nonuniform. For example, the Wagner number when measuring Cu²⁺/Cu⁺ kinetics in choline chloride–ethylene glycol DES (1:4 molar ratio of ChCl:EG) is very low (<0.1), indicating highly nonuniform current distribution over the disk electrode. We show here that the Cu²⁺/Cu⁺ exchange current density measured using disk electrodes is very inaccurate due to the aforementioned nonuniform current distribution. To obtain uniform disk current distribution, we employ here a coplanar concentric rotating ring-disk electrode (RRDE), where the ring serves the function of a current “shield.” Specifically, we show using modeling that the ring minimizes the current distribution nonuniformity at the disk by effectively shielding the disk against current spikes near the disk edge. This enables improved precision in electrode kinetics measurements for the Cu²⁺/Cu⁺ couple in resistive DES. To enable broad applicability of this technique, an analytical expression based on the Wagner number is integrated into an iterative algorithm to help users identify ring conditions to achieve uniform current distribution and thus improved electroanalytics at the disk.

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

Sulfur-functionalized solid-phase materials for the selective separation of arsenic and selenium

Radioactive arsenic (As) isotopes are of growing interest for applications in nuclear medicine, national security, and environmental research. Recent efforts at the Facility for Rare Isotope Beams (FRIB) have focused on aqueous harvesting of selenium-72,73 ( 72,73 Se) and their daughter isotopes, arsenic-72,73 ( 72,73 As), which are particularly valuable for medical applications and nuclear data studies, respectively. Both conventional isotope production and harvesting methods require chemical separations to purify radioactive As from parent and co-produced Se radioisotopes. While several solid-phase separation methods for As and Se exist, many depend on complex oxidation state control or highly acidic conditions. This study presents results for sulfur-based solid-phase materials selected to enable uptake at lower acidity and eliminate the need for intricate redox chemistry. Specifically, the performance of three covalently bound sulfur-based ligands were evaluated: (1) thiophenol-polystyrene, (2) propanethiol-silica, and (3) thiourea-silica. Uptake characteristics—including distribution coefficients (Dw), kinetics, and column separation behavior—were assessed using 75 Se and 73 As in hydrochloric (HCl) acid and nitric (HNO 3 ) solutions. The resins demonstrated high-yield (>95%) and high-purity As recovery across a range of HCl concentrations. Comparable results in HNO 3 were achieved when combined with anion exchange chromatography. Furthermore, the potential application of these materials for medical isotope generators was also investigated through ligand stability and repeated elution studies. Overall, sulfur leaching from the resins was negligible at the concentrations relevant for these separations but increased with higher acid concentrations.

Arsenic↗
DOE OSTI2025

The role of the Pt-group dehydrogenation catalyst in alkane metathesis for polyolefin deconstruction

Recent proposed approaches in the depolymerization of waste plastics employ an olefin intermediate to produce alkanes or alkenes using olefin metathesis in tandem chemistry. Here, in this study, we investigated the role of the dehydrogenation catalyst on reaction rate, kinetics, and product distribution in heterogeneous tandem dehydrogenation and olefin metathesis (alkane metathesis) of three different alkane reactants, including polyethylene. We found that many properties to which alkane dehydrogenation rates were sensitive-including metal composition, nanoparticle size, and surface doping of Re species also controlled activity in Tandem D/OM. When comparing Pd, Pt, and Pt 3 Sn 1 , supported Pd in tandem with a Re 2 O 7 olefin metathesis catalyst showed four-fold higher activity (surface area basis) compared to Pt or Pt 3 Sn 1 catalysts on the same support, mainly due to differences in the rate of hydrogenation. Catalyst preparation resulted in metal nanoparticles partially covered by ReO x , as seen from elemental mapping. Co-location of Re 2 O 7 and Pd correlated with increased rates of hydrogenation (i.e., an increase in the rate of alkane formation and simultaneous lowering of the rate of alkene formation), with a reaction order in catalyst study that further supported this conclusion. The Pd and Re 2 O 7 system displayed marked improvement compared to Pt or Pt 3 Sn 1 with Re 2 O 7 , and previous work, in the depolymerization rate of a linear polyethylene feedstock, with over 94 % reduction in polymer molecular weight in 15 h at 190 °C using less catalyst and increased reactant loadings, while keeping solvent to polymer consumption below 2.5.

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

Quantification of [MTBD][beti] loading and its effects on Pt/HSC electrode performance in hydrogen fuel cells

High surface area carbon-supported platinum catalysts (Pt/HSC) are widely used in polymer electrolyte membrane fuel cells but often suffer from limited proton and oxygen transport within porous domains. To address these challenges, we integrate the ionic liquid [MTBD][beti] into Pt/HSC catalyst layers using two deposition methods-one-pot and sequential deposition-to tune IL distribution within micropores and mesopores. A combination of ex situ and in operando techniques were employed to elucidate electrode structure-performance relationships across a range of IL loadings. Sequential deposition achieved more efficient pore filling and higher IL retention at lower IL:C ratios, enabling enhanced proton conductivity and protection of active sites. Compared to the IL-free Pt/HSC, IL-modified electrodes demonstrated up to 28% improvement in mass activity and enhanced high-current-density performance under low relative humidity, while maintaining comparable performance under humidified conditions. Electrochemical impedance spectroscopy and CO displacement experiments reveal that improvements are linked to reduced ionic resistance and lower sulfonate adsorption on Pt sites, rather than changes in electrochemical surface area. However, excessive IL loading leads to mass transport losses. These results highlight the importance of selective pore filling and efficient IL distribution in achieving kinetic gains without compromising ionic and gas transport.

08 HYDROGEN↗