Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Density gradient function”

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 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Nucleon localization function in rotating nuclei

An electron localization function was originally introduced to visualize in positional space bond structures in molecules. It became a useful tool to describe electron configurations in atoms, molecules, and solids. In nuclear physics, a nucleon localization function (NLF) has been used to characterize cluster structures in light nuclei, formation of fragments in fission, and pasta phases appearing in the inner crust of neutron stars. In the case of a deformed rotating nucleus, several NLFs can be introduced, depending on the definition of the spin-quantization axis, direction of the total angular momentum, and self-consistent symmetries of the system. Contributions to the NLF from the current density, spin-current tensor density, and density gradient terms are negligible in the nuclear interior. The oscillating pattern of the simplified NLF can be explained in terms of a constructive interference between kinetic-energy and particle densities. The characteristic nodal pattern seen in the NLF in the direction of major axis of a rotating nucleus comes from single-particle orbits carrying large aligned angular momentum. The variation of the NLF along the minor axis of the nucleus can be traced back to deformation-aligned orbits.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Investigating trends in actinide covalency and magnetism with 35/37 Cl SSNMR spectroscopy and first-principles calculations

Solid-state NMR (SSNMR) spectroscopy is a powerful technique for studying actinide chemistry but has been significantly limited due to the complex paramagnetism, and radiological hazards presented by these materials. Lanthanide and actinide salts often feature magnetic ordering and can be paramagnetic, ferromagnetic, or antiferromagnetic depending on temperature and electronic structure. Paramagnetic interactions can manifest in SSNMR both as secular spectral shifts and/or couplings as well as contributions from non-secular relaxation. Both effects can be directly measured with NMR and used to extrapolate rich chemical information such as coordination environments, bonding characteristics, local molecular dynamics, and correlation times. Typically, these studies are carried out on high-γ and highly abundant NMR-active isotopes (e.g., 1 H, 6/7 Li, 19 F, 23 Na, etc.) or on enriched rare isotopes (e.g., 2 H and 17 O), which can be expensive. Herein, we present a facile methodology to measure the 35/37 Cl electric-field gradient (EFG) and paramagnetic shift anisotropy (SA) tensor components using static wideline SSNMR measurements of LaCl 3 , NdCl 3 , UCl 3 , and UCl 4 . The static powder spectra were measured with both 35 Cl and 37 Cl SSNMR to increase the fidelity of the extracted tensor parameters. Variable temperature NMR of a select case confirms the Curie–Weiss paramagnetism. Relaxation measurements of both nuclei further corroborate observations owing to the paramagnetic relaxation enhancement and reveal simultaneous quadrupolar relaxation mechanisms. Density functional theory (DFT) calculations using Hubbard U corrections to the uranium valence orbitals show excellent agreement with experimental EFG tensor parameters and help describe the bonding characteristics in these lanthanide and actinide systems.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

Analysis of Hydrocyclone Geometry via Rapid Optimization Based on Computational Fluid Dynamics

Abstract Hydrocyclones exploit density gradients for the centrifugal separation of dispersions in a continuous liquid. Selection of the geometrics for optimal separation is case specific, like the media characteristics. The existing optimization method based on computational fluid dynamics (CFD) provides a powerful analytical tool but requires long computational times. The most common praxis for CFD optimization is via the single‐factor optimization method (SFOM). In this study, a novel approach is presented as an improved rapid optimization method that implements a dynamic‐mesh and user‐defined function optimization method (DUOM). The DUOM adapts the dynamic‐mesh approach from other applications to the optimization analysis of hydrocyclones. The DUOM reduced the computational time by 31.1 %, compared to the SFOM.

Liu, Lin↗
DOE OSTI2022

Oxygen–Chlorine Chemisorption Scaling for Seawater Electrolysis on Transition Metals: The Role of Redox

To clarify what controls species oxidation selectivity in seawater electrolysis, density functional theory (DFT) is used to identify chemisorption enthalpy trends and scaling relations for the simplest relevant adsorbates (O, Cl, and H) on relevant surfaces of 3d transition metals, as well as Pd and Pt, in face-centered-cubic and, if different, their ground-state crystal structures. Approximations are tested for electron exchange-correlation (XC) and van der Waals interactions to assess their ability to reproduce experimental adsorption enthalpies of H and O on Pt(111). The vdW-uncorrected generalized gradient approximation to XC of Perdew, Burke, and Ernzerhof (PBE) agrees most closely with experiments. Using DFT-PBE thereafter, it is determined that the O chemisorption enthalpy on this wide range of transition-metal surfaces is proportional to the sum of first and second atomic ionization energies, akin to a Born–Haber cycle for a redox reaction, indicating that metal redox activity controls O chemisorption strength. Then it is shown that the O and Cl chemisorption enthalpies are strongly correlated, suggesting that the transition metals considered will oxidize unselectively water and Cl – . This strong correlation appears also for crystal reduction potentials of binary oxides and chlorides, indicating a fundamental challenge for future seawater electrode materials design.

08 HYDROGEN↗
DOE OSTI2024

Development of an Optical Library for Coevaporated CdSe x Te 1– x

The conversion efficiency of CdTe solar cells may be improved by bandgap engineering, i.e., changing the bandgap value through the addition of Se in the absorber. The Se alloying enables a short-circuit current density improvement, as it leads to a bandgap energy value decrease. Furthermore, it has been associated with increased minority carrier lifetimes, assuring high open-circuit voltage values. An Se gradient profile control can further optimize the solar cell performance. Thus, an optical model baseline of the CdSe x Te 1–x (CST) compound was developed. Spectroscopic ellipsometry measurements were conducted to accurately extract the optical constants of ten CST layers deposited through coevaporation with x varying from 0 to 1. Using the measured dielectric function spectra from the discrete CST layers with varying x, and considering the composition-induced shift in the critical point energies, an energy-shift model was employed to develop the accurate optical library for the CST compound for any x value to provide data for future modeling and optimization. Furthermore, the library accuracy was validated through optical simulations of the quantum efficiency of a graded CST solar cell using the finite-difference time-domain method by replicating the Se profile in the absorber layer measured through secondary ion mass spectrometry.

14 SOLAR ENERGY↗
DOE OSTI2025

Beyond Magic Barrels: Digital manufacturing for crystallization, process development and optimization of explosive materials: Part II Resveratrol Exemplar

This SAND report summarizes work supported by an Engineering Sciences Research Foundation (ESRF) Lab Directed Research and Development (LDRD) project entitled “Beyond Magic Barrels: Digital manufacturing for crystallization, process development and optimization of explosive materials.” This SAND report is written in two parts with Part 1 discusses recrystallization of our explosive exemplar and Part 2 summarizing our work with recrystallization of resveratrol. We have studied resveratrol recrystallization with a multiscale approach combining experiments, modeling and simulation. At the single crystal scale, microscopy experiments illuminate crystal time-dependent growth rates using advanced image analysis. Bench scale experiments were carried out to look at growth of multiple particles in a small reactor creating thousands of particles and analyzing the results with microscopy and μCT. For the modeling we combine kinetic Monte Carlo (kMC) models with subscale information from density functional theory (DFT) or molecular dynamics. This work is discussed in Part 1 and can also be found in a paper from the project discussing a coarse-grained kMC model specifically developed for resveratrol. For well-mixed systems, we have population balance equations (PBE) linked with species mass conservation forming a set of ordinary differential equations that can be solved quickly. For more complicated geometries, such as the vat crystallization used throughout the complex, a coupled computational fluid dynamic (CFD)/PBE method was developed to account for gradients in temperature and concentration and differences in crystallization rates throughout the domain. These simulations are more complex and require high performance computing. We present results for two cases: 5% seed fast cool with parameters fit to the well-mixed case and 5% seed slow cool using the same parameters. We show reasonable agreement with experiments though are particles are significantly larger than the experiments.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Ten-moment fluid model for low-temperature magnetized plasmas

In this paper, a one-dimensional 10-moment multi-fluid plasma model is developed and applied to low-temperature magnetized plasmas. The 10-moment model solves for six anisotropic pressure terms, in addition to density and three components of fluid momentum, which allows the model to capture finite kinetic effects. The results are benchmarked against a 5-moment model, which assumes that the gas constituents follow a Maxwellian velocity distribution function (VDF), and a particle-in-cell Monte Carlo collision model that allows for arbitrary non-Maxwellian VDFs. The models are compared in a one-dimensional, low-temperature, partially magnetized plasma test case. The 10-moment results accurately reproduce the anisotropic temperature profile in low-temperature magnetized plasmas, where shear gradients exist due to the E×B drift. We discuss the mechanisms by which the anisotropic pressure can be generated in low-temperature magnetized plasmas. In addition, the importance of a self-consistent heat flux closure to the 10-moment model is studied, showing consistency with other models only when the assumptions of the underlying model are met. The 10-moment model allows for study of electron inertia effects and non-Maxwellian VDFs without the need for kinetic methods that are more computationally expensive.

Kuldinow, Derek Amur (ORCID:0000000319730196)↗
DOE OSTI2023

Ab Initio Evaluation of Plutonium Dioxide $S(α,β)$ and Thermal Neutron Cross Sections

Plutonium Dioxide (PuO 2 ) is an advanced fuel-cycle nuclear fuel material for thermal and fast reactors, with a melting point at 3261.8 K and thermal conductivity above 3 Wm -1 K -1 for temperatures below 1500 K. Plutonium from spent UO 2 is reprocessed to form chemically stable mixed oxide (MOX) fuel. In contrast to metallic fuels, the low PuO 2 conductivity results in high thermal gradients between fuel center and surface, resulting in structural annealing and efficient fission product migration over burnup. Crystal binding affects the doppler broadening of epithermal resonances for nuclear fuels and scattering behavior of thermal neutrons. Currently, there is no ENDF/B Thermal Scattering Law, i.e., $S(α,β)$, evaluation for PuO 2 , only free atom cross sections for Plutonium and Oxygen, respectively, which neglect lattice contributions to total cross section. Herein, ab initio lattice dynamics (AILD) techniques are employed to calculate the phonon density of states (DOS) using spin-orbit-coupling density functional theory (DFT) to predict the PuO2 paramagnetic ground-state structure. These results highlight the evaluation of $S(α,β)$ for Plutonium and Oxygen, respectively in PuO 2 and consequential generation of thermal neutron scattering cross sections for high fidelity criticality safety analysis and reactor calculations.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Evaluating the gradients of localized diabatic state energies and couplings at minimum cost

We calculate the gradients of Boys diabatic state energies and couplings when the electronic vector space is generated by configuration interaction singles. Our approach follows the Lagrangian approach of Paz and Glover (rather than direct differentiation of the adiabatic-to-diabatic approaches that have been published previously). The result is that we achieve a dramatic increase in savings over previous approaches, and the present approach should be immediately useful to scientists focused on electronic relaxation, especially chemists studying electron transfer who wish to go beyond the Condon approximation. Here, a future extension to time-dependent density functional theory in the Tamm–Dancoff approximation is clear.

Chemical physics↗
DOE OSTI2020

Ground-state properties of the Hubbard model in one and two dimensions from the Gutzwiller conjugate gradient minimization theory

We introduce Gutzwiller conjugate gradient minimization (GCGM) theory, an ab initio quantum many-body theory for computing the ground-state properties of infinite systems. GCGM uses the Gutzwiller wave function but does not use the commonly adopted Gutzwiller approximation (GA), which is a major source of inaccuracy. Instead, the theory uses an approximation that is based on the occupation probability of the on-site configurations, rather than approximations that decouple the site-site correlations as used in the GA. We test the theory in the one-dimensional and two-dimensional Hubbard models at various electron densities and find that GCGM reproduces energies and double occupancies in reasonable agreement with benchmark data at a very small computational cost.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Focusing ions at atmospheric pressure using nonlinear DC voltage sequences applied to a stacked ring ion guide

Many modern ion mobility (IM) and mass spectrometers (MS) operate under low pressure (= 10 Torr) and employ high voltage radiofrequencies (RF) to provide ion confinement. Unfortunately, RF effectiveness drastically decreases as pressure increases, and few techniques for focusing ions at elevated pressures exist. Here we demonstrate a new approach for focusing ions at atmospheric pressure (AP) by applying nonlinear DC voltage sequences following quadratic and power (exponential) functions to a stacked ring ion guide. In this work, we used ion trajectory simulations to rigorously explore how ions react to nonlinear electric fields and validate the simulations with a set of ion current measurements performed at AP. Ion trajectory simulations show that ions initially defocus near the entrance of the device but then become intensely focused as they travel through the device. Contour plots for both nonlinear voltage sequences show electric field lines that increasingly curve inwards as a function of distance, resulting in spatial ion focusing. Experimental ion current and spot size measurements were performed at AP using a 10-cm stacked ring ion guide and a segmented Faraday cup detector. Quadratic sequences produced ~5% smaller spot sizes (~22.8 mm) and ~25% higher ion current compared to a linear voltage sequence (~24.0 mm). Alternatively, power sequences produced ~64% smaller spot sizes (~8.7 mm), albeit with ~10x lower ion current. However, both nonlinear voltage sequences produced similar ion currents at the center of the Faraday cup detector, indicating that higher ion densities are achieved when using nonlinear voltage gradients. These results demonstrate a new way to focus ions at AP, and the capabilities demonstrated here provide fundamental insights on how to keep ions inside analytical devices at elevated pressures without RF.

47 OTHER INSTRUMENTATION↗
DOE OSTI2022

Focusing Ions at Atmospheric Pressure using Nonlinear DC Voltage Sequences Applied to a Stacked Ring Ion Guide

Many modern ion mobility (IM) and mass spectrometers (MS) operate under low pressure (= 10 Torr) and employ high voltage radiofrequencies (RF) to provide ion confinement. Unfortunately, RF effectiveness drastically decreases as pressure increases, and few techniques for focusing ions at elevated pressures exist. Here we demonstrate a new approach for focusing ions at atmospheric pressure (AP) by applying nonlinear DC voltage sequences following quadratic and power (exponential) functions to a stacked ring ion guide. We used ion trajectory simulations to rigorously explore how ions react to nonlinear electric fields and validate the simulations with a set of ion current measurements performed at AP. Ion trajectory simulations show that ions initially defocus near the entrance of the device but then become intensely focused as they travel through the device. Contour plots for both nonlinear voltage sequences show electric field lines that increasingly curve inwards as a function of distance, resulting in spatial ion focusing. Experimental ion current and spot size measurements were performed at AP using a 10-cm stacked ring ion guide and a segmented Faraday cup detector. Quadratic sequences produced ~5% smaller spot sizes (~22.8 mm) and ~25% higher ion current compared to a linear voltage sequence (~24.0 mm). Alternatively, power sequences produced ~64% smaller spot sizes (~8.7 mm), albeit with ~10x lower ion current. However, both nonlinear voltage sequences produced similar ion currents at the center of the Faraday cup detector, indicating that higher ion densities are achieved when using nonlinear voltage gradients. These results demonstrate a new way to focus ions at AP, and the capabilities demonstrated here provide fundamental insights on how to keep ions inside analytical devices at elevated pressures without RF.

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

Assessing the Limitations of Self-Interaction-Corrected Functionals for Describing the Hydrated Electron

Simulating the hydrated electron using density functional theory is challenging due to the prevalence of self-interaction error in standard functionals. Hybrid functionals like PBEh(40) can reasonably describe the chemistry of an excess electron in water and partially mitigate self-interaction error by incorporating exact Hartree–Fock exchange, but they are computationally expensive making them impractical for large-scale and long-time ab initio molecular dynamics simulations. Explicit self-interaction correction schemes that are applied on an orbital-by-orbital basis offer a potential alternative when the correction is limited to the singly occupied molecular orbital obtained with a generalized gradient approximation functional. Here, we examine whether the Perdew–Zunger self-interaction correction scheme applied to the revPBE functional can provide a computationally efficient and physically sensible alternative to PBEh(40) for the hydrated electron. We find that functionals incorporating a self-interaction correction scheme should be viewed with caution when applied to the hydrated electron and its reactivity. Furthermore, we show that it is critical to consider extensive sampling and diverse chemical environments when validating their performance.

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

Collisionality driven turbulent particle transport changes in DIII-D H-mode plasmas

The results of the experimental dimensionless scan in this paper confirm that there is an increase in density peaking towards lower collisionality and that this can be partly linked to a shift in the turbulence regime from ITG towards TEM. However at the lowest collisionality, the changes in turbulence and transport are much more pronounced than expected from direct collisionality effect on the turbulence. In this paper, the collisionality, ν * is varied by a factor 5, while keeping ρ * , q , β , M , fixed. Additionally, a 3 Hz gas puff modulation is applied to modulate the electron density profile and extract the perturbed transport coefficients using two diagnostics. The transport analysis shows that the increase in density peaking at low ν * is linked to an increase in the inward particle pinch and not an increase in core fueling. These observations are not only in agreement with prior modeling scans of how turbulence changes as a function of collisionality and its impact upon the particle fluxes, but also with the multi-machine database (Fable E. et al 2010 Plasma Phys. Control. Fusion 52 015007) (Angioni C. et al 2003 Phys. Rev. Lett. 90 205003). The changes in turbulence across the collisionality scan were captured at large scale by the BES and at smaller scale by the DBS. A comparison with gradient-driven GENE simulations showed similar trends at both scales. Moreover, the changes observed in overall transport are in agreement with gradient-driven TGLF particle flux simulations. This indicates that TGLF/GENE when given the gradients as input, are able to reproduce the experimentally observed turbulence changes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Tree-inspired lignin microrods-based composite heterogeneous nanochannels for ion transport and osmotic energy harvesting

One of the key processes of tree lignification is that lignin penetrates into the cell wall and fills in the cell wall framework, thereby increasing the hardness and hydrophobic of the tree channels, which is beneficial to consolidate and support the tree cell wall and water transport. Inspired by this natural process, we demonstrated a lignin-based nanofluidic heterogeneous membrane that closely mimics the channels in tree, which can realize ion transport function and effectively capture reverse electrodialysis. The membrane was synthesized by heating dealkaline lignin and PVA at 200 °C and this formed fusiform microrods and a closed-packed membrane. Simultaneously, this membrane composited with anodized aluminum (AAO) channels membrane at 200 °C to form asymmetric heterogeneous nanochannels membrane, which can transport counter-ions and harvest osmotic energy. This membrane implements ion current rectification in 0.1 M KCl electrolyte solution at pH 3 due to the confinement of pores and opposite surface charges in lignin-based heterogeneous nanochannel. An output power density of 0.97 W m- 2 is obtained under a 50-fold salinity gradient, which can be further improved to 1.19 W m -2 by increasing the salinity gradient from 50-fold to 500-fold. Accordingly, this nanofluidic membranes were prepared by using lignin, the key component in tree, which not only mimicked a crucial process of the water and ionic transport process of channels in tree, but also had the prospect in the field of osmotic energy harvesting.

09 BIOMASS FUELS↗
DOE OSTI2024

Quantification of grain boundary effects on the geometrically necessary dislocation density evolution and strain hardening of polycrystalline Mg 4Al using in situ tensile testing in scanning electron microscope and HR-EBSD

In situ tensile testing in a scanning electron microscope (SEM) in conjunction with high-resolution electron backscatter diffraction (HR-EBSD) under load was used to characterize the evolution of geometrically necessary dislocation (GND) densities at individual grain boundaries as a function of applied strain in a polycrystalline Mg—4Al alloy. The increase in GND density was investigated at plastic strains of 0 %, 0.6 %, 2.2 %, 3.3 % from the area including 76 grains and correlated with (i) geometric compatibility between slip systems across grain boundaries, and (ii) plastic incompatibility. We develop expressions for the grain boundary GND density evolution as a function of plastic strain and plastic incompatibility, from which uniaxial tensile stress-strain response of polycrystalline Mg—4Al are computed and compared with experimental measurement. The findings in this study contribute to understanding the mechanisms governing the strain hardening response of single-phase polycrystalline alloys and more reliable prediction of mechanical behaviors in diverse microstructures.

Geometrically necessary dislocations↗
DOE OSTI2025

Score-Based Physics-Informed Neural Networks for High-Dimensional Fokker–Planck Equations

The Fokker-Planck (FP) equation is a foundational partial differential equation (PDE) in stochastic processes involving Brownian motions. However, the curse of dimensionality (CoD) poses a formidable challenge when dealing with high-dimensional FP equations. Although Monte Carlo simulation and (vanilla) Physics-Informed Neural Networks (PINNs) have shown the potential to tackle CoD, both methods exhibit significant numerical errors in high dimensions when dealing with the probability density function (PDF) associated with Brownian motion. The point-wise PDF values tend to decrease exponentially as dimensionality increases, surpassing the precision of numerical simulations and resulting in substantial errors. In addition, due to its massive sampling, Monte Carlo fails to offer fast sampling. Modeling the logarithm likelihood (LL) via vanilla PINNs transforms the FP equation into a notoriously difficult Hamilton-Jacobi-Bellman (HJB) equation, which is impractical for PINN learning, whose error grows rapidly with dimension. To this end, we propose a novel approach utilizing a score-based solver to fit the score function in stochastic differential equations (SDEs). The score function, defined as the gradient of the LL, plays a fundamental role in inferring LL and PDF and enables fast SDE sampling, offering an effective means to overcome the CoD. Three fitting methods, Score Matching (SM), Sliced Score Matching (SSM), and Score-PINN, are introduced, each contributing unique advantages in computational complexity, accuracy, and generality. The proposed score-based SDE solver operates in two stages: first, employing score matching or Score-PINN to acquire the score function; and second, solving the LL via an ordinary differential equation (ODE) using the obtained score function. Comparative evaluations across these methods showcase varying trade-offs. The proposed methodology is evaluated across diverse SDEs, including anisotropic Ornstein-Uhlenbeck processes, geometric Brownian motion, and Brownian motion with varying eigenspace. We also test various distributions, including Gaussian, Log-normal, Laplace, and Cauchy distributions. The numerical results demonstrate the score-based SDE solver’s stability, speed, and performance across different experimental settings, solidifying its potential as a solution to CoD for high-dimensional FP equations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Atomic Gradient Structure Alters Electronic Structure in 3D across the Bulk and Enhances Photoactivity

Abstract Understanding structure–function relationships enables the design of materials with tailored functionalities. The long‐standing challenge is to design materials with high active volume to improve efficiency. Tailored grain boundaries and lattice defects are traditionally used to tune the electronic structure near interfaces or defects to promote electron and hole separation. However, the active volume of point defect sites or interfaces in these traditional photocatalysts is extremely low. This study reports a structure with continuous atomic positional deformation across the bulk, altering the electronic structure in 3D and creating a significantly higher active volume. Such a structure in anatase is obtained and tuned by phase transformation during the heating process. Transmission electron microscopy and density functional theory results reveal that atomic deformations result in continuous band bending across the particles, facilitating electron–hole separation, inhibiting their recombination, and inducing dramatically enhanced photoactivity. These findings enable a different materials design paradigm that can potentially be harnessed for a broad range of applications.

Ren, Peng↗