Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “density distribution”

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

Passivation‐Induced Species Dynamics and Microstructural Evolution in Solid‐State Lithium–Sulfur Cathodes

Solid-state lithium–sulfur (SSLS) batteries offer high theoretical energy density, yet their practical viability is hindered by poor sulfur utilization and limited rechargeability. At the core of this challenge lies the passivating nature of Li 2 S, which restricts ionic and electronic transport, suppresses interfacial activity, and severely impedes the reversibility of electrochemical reactions. In this study, we elucidate the mechanistic origins of these limitations by resolving how charge and discharge species form, grow, and spatially evolve within the cathode microstructure under varied current densities and electrode compositions. By resolving the species distribution at the particle scale and coupling it with Raman spectroscopy and X-ray diffraction, we demonstrate how Li 2 S formation induces localized surface passivation that progressively limits electrochemical accessibility within the cathode microstructure. Sulfur utilization is found to be strongly governed by the interplay between sulfur loading, residual porosity, and interfacial architecture. High sulfur contents result in buried, electrochemically isolated domains due to poor solid electrolyte (SE) percolation, while low sulfur contents trigger SE degradation via parasitic reactions. The resulting sulfur-porosity maps delineate the mechanistic boundaries between reversible and transport-limited regimes, offering actionable design guidance for SSLS cathodes with enhanced sulfur utilization.

electrode microstructure↗
DOE OSTI2022

Laser Powder Bed Fusion of ODS 14YWT from Gas Atomization Reaction Synthesis Precursor Powders

Abstract Laser powder bed fusion (LPBF) additive manufacturing (AM) is a promising route for the fabrication of oxide dispersion strengthened (ODS) steels. In this study, 14YWT ferritic steel powders were produced by gas atomization reaction synthesis (GARS). The rapid solidification resulted in the formation of stable, Y-containing intermetallic Y 2 Fe 17 on the interior of the powder and a stable Cr-rich oxide surface. The GARS powders were consolidated with LPBF. Process parameter maps identified a stable process window resulting in a relative density of 99.8%. Transmission electron microscopy and high-energy x-ray diffraction demonstrated that during LPBF, the stable phases in the powder dissociated in the liquid melt pool and reacted to form a high density (1.7 × 10 20 /m 3 ) of homogeneously distributed Ti 2 Y 2 O 7 pyrochlore dispersoids ranging from 17 to 57 nm. The use of GARS powder bypasses the mechanical alloying step typically required to produce ODS feedstock. Preliminary mechanical tests demonstrated an ultimate tensile and yield strength of 474 MPa and 312 MPa, respectively.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Axially-deformed solution of the Skyrme-Hartree-Fock-Bogoliubov equations using the transformed harmonic oscillator basis (IV) HFBTHO (v4.0): A new version of the program

We describe the new version 4.0 of the code HFBTHO that solves the nuclear Hartree-Fock-Bogoliubov problem by using the deformed harmonic oscillator basis in cylindrical coordinates. In the new version, we have implemented the restoration of rotational, particle number, and reflection symmetry for even-even nuclei. The restoration of rotational symmetry does not require using bases closed under rotation. Furthermore, we added the SeaLL1 functional and improved the calculation of the Coulomb potential. Finally, we refactored the code to facilitate maintenance and future developments.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Surprisingly high irradiation-induced defect mobility in Fe 3 O 4 as revealed through in situ transmission electron microscopy

Radiation-induced defects are expected to change the corrosion kinetics when a material is exposed to a corrosive environment. Characterization of the defects forming under irradiation and their ability to move/diffuse is therefore key to understanding how they impact the corrosion kinetics and possibly mechanisms. For that matter, Fe/Fe 3 O 4 interfaces were irradiated using 1 MeV Kr 2+ at temperatures of 273 K (25 °C), 573 K (300 °C) and 773 K (500 °C) to a maximum dose of 10 dpa (as measured in the Fe matrix) in situ in a transmission electron microscope. The evolution of the microstructure was followed as a function of the dose and temperature. The formation of defects was followed in the metal and the oxide as well as at the metal/oxide interface and characterized in terms of size and density. Additionally, the dynamic observation of the spatial distribution of the defects allowed for a qualitative assessment of the mobility of defects in the oxide as a function of temperature, revealing evidence for enhanced mobility even more so than in the Fe matrix. Mechanisms of induced defect mobility in the spinel are discussed to account for the observations.

36 MATERIALS SCIENCE↗
DOE OSTI2022

EDGE2D-EIRENE and ERO2.0 predictions of nitrogen molecular break-up and transport in the divertor of JET low-confinement mode plasmas

EDGE2D-EIRENE simulations of nitrogen-seeded partially detached JET L-mode plasmas show that the divertor N I to N V radiation distributions are highly sensitive to the upstream electron density, and less sensitive to changes in the assumed cross-field particle diffusivity for nitrogen ions. The EDGE2D-EIRENE simulations reproduce the peak intensities of N I to N V as measured by vertically viewing divertor filterscopes to within 50 % for a narrow range of the upstream electron density within the experimental uncertainties, while the predicted profiles are narrower than the measured ones. Including nitrogen atoms only in the ERO2.0 simulations implies lower N III and N IV peak intensities by one and two thirds, respectively, compared to EDGE2D-EIRENE. If nitrogen is assumed to recycle exclusively as molecules instead of atoms, ERO2.0 predicts that the N III and N IV intensities in the divertor increase by up to a factor of two and that the time-averaged, volume-integrated number of N 5+ to N 7+ ions in the plasma also increases by approximately a factor of 2.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Charge transfer effect on local lattice distortion in a HfNbTiZr high entropy alloy

We report it is often assumed that atoms are hard spheres in the estimation of local lattice distortion (LLD) in high-entropy alloys (HEAs). However, our study demonstrates that the hard sphere model misses the key effect, charge transfer among atoms with different electronegativities, in the understanding of the severely distorted local structure in the body-centered cubic (BCC) HfNbTiZr HEA. Through the characterization and simulations of the local structure of the HfNbTiZr HEA, we found that the LLD regions in HfNbTiZr have larger average size compared with face-centered cubic (FCC) HEAs, but the magnitude of LLD significantly varies from one atomic site to another because of electron transfer from larger atoms, Hf and Zr, to smaller ones, Nb and Ti. Our finding may form the basis for the design of high performance HEAs through tuning their site-to-site chemical heterogeneity.

36 MATERIALS SCIENCE↗
DOE OSTI2026

A novel conditional formulation of the Vlasov–Ampère equations: a conservative, positivity, asymptotic and Gauss law preserving scheme

We propose a novel reformulation of the Vlasov–Ampère equations for plasmas that reveals discrete symmetries that enables simultaneous conservation of mass, momentum and energy; preservation of Gauss’s law; positivity of the distribution function; and consistency with quasi-neutral asymptotics. The approach employs variable and coordinate transformations to yield a coupled system comprising a modified Vlasov equation and associated moment–field equations. The modified Vlasov equation advances a conditional distribution function that excludes mass, momentum and energy densities, which are instead evolved through moment equations enforcing the relevant symmetries, conservation laws and involution constraints. This reformulation aligns naturally with a recent slow-manifold reduction technique, which separates fast electron time scales and simplifies the treatment of the quasi-neutral limit within the reduced moment–field subsystem. Using this framework, we develop a numerical method for the reduced 1D1V subsystem that, for the first time in the literature, satisfies all key physical constraints while maintaining a quasi-neutral asymptotic behaviour. The advantages of the method are demonstrated on canonical electrostatic test problems, including the multiscale ion acoustic shock wave.

1D1V↗
DOE OSTI2021

From Femtoseconds to Gigaseconds: The SolDeg Platform for the Performance Degradation Analysis of Silicon Heterojunction Solar Cells

Heterojunction Si solar cells exhibit notable performance degradation. Here, we modeled this degradation by electronic defects getting generated by thermal activation across energy barriers over time. To analyze the physics of this degradation, we developed the SolDeg platform to simulate the dynamics of electronic defect generation. First, femtosecond molecular dynamics simulations were performed to create a-Si/c-Si stacks, using the machine learning-based Gaussian approximation potential. Second, we created shocked clusters by a cluster blaster method. Third, the shocked clusters were analyzed to identify which of them supported electronic defects. Fourth, the distribution of energy barriers that control the generation of these electronic defects was determined. Fifth, an accelerated Monte Carlo method was developed to simulate the thermally activated time-dependent defect generation across the barriers. Our main conclusions are as follows. (1) The degradation of a-Si/c-Si heterojunction solar cells via defect generation is controlled by a broad distribution of energy barriers. (2) We developed the SolDeg platform to track the microscopic dynamics of defect generation across this wide barrier distribution and determined the time-dependent defect density N(t) from femtoseconds to gigaseconds, over 24 orders of magnitude in time. (3) We have shown that a stretched exponential analytical form can successfully describe the defect generation N(t) over at least 10 orders of magnitude in time. (4) We found that in relative terms, Voc degrades at a rate of 0.2%/year over the first year, slowing with advancing time. (5) We developed the time correspondence curve to calibrate and validate the accelerated testing of solar cells. We found a compellingly simple scaling relationship between accelerated and normal times t normal ∝ t accel T(accel)/T(normal) . (6) We also carried out experimental studies of defect generation in a-Si:H/c-Si stacks. We found a relatively high degradation rate at early times that slowed considerably at longer time scales.

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

Multiscale Mechanisms of Twisted Carbon Nanotube Yarns Probed In Situ by Soft X-rays during Tensile Loading

Piecing together carbon nanotubes (CNTs) into assemblies has so far failed to achieve the same elite strength performance metrics as individual CNTs, highlighting a critical deficiency in understanding the effects that the processing of individual nanostructures have on the performance of their derived macroscale assemblies, thereby hindering the development of a process-structure-performance map for these materials. Here, in this work, we propose a method to decouple the distribution orientation of nanoscale tortuosity and the microscale twist of CNT dry-spun yarns under applied loads via in situ soft X-ray probing at high energy (1200 eV) and low energy (280 eV), respectively. With this decoupling enabled by in situ soft X-ray scattering, we acquired a deeper understanding of the deformation mechanisms of these yarns. We found that for untreated yarns, the twist angle of collective CNT bundles at the macroscale is more sensitive to applied stress than the nanoscale alignment distribution. We also found that increasing nominal twist densities of yarns as well as increased strengthening via plasma treatments and polymer infiltration act to decrease the yarns’ sensitivity to realignment at the nanoscale and prevent failure by the slip mechanism.

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

Material Characterization-Based Wear Mechanism Investigation for Biomass Hammer Mills

Biomass, as harvested, is composed of inorganic compounds both intrinsically and extrinsically and can be abrasive. The present study investigates the wear modes and mechanisms of two types of blades of hammer mills used in biomass size reduction (impacting the particle size and distribution) and densification (impacting the size, shape, and density). The dominant wear modes for the stage 1 steel blades are determined to be erosive and polishing wear. For the stage 2 blades with a carbide weld overlay, the main wear mechanisms are erosion and fracture. Partial replacement of Co by Fe in the outer layer of the carbide grits, likely induced by diffusion during high-temperature welding, has been correlated to the observed microcracking. Finally, the microcracking is believed to weaken the grit strength and fracture toughness to make the overlay prone to fracture and erosion due to repetitive contact with the inorganic contents in chopping biomass.

42 ENGINEERING↗
DOE OSTI2024

Correlated Anion Disorder in Heteroanionic Cubic TiOF 2

Resolving anion configurations in heteroanionic materials is crucial for understanding and controlling their properties. For anion-disordered oxyfluorides, conventional Bragg diffraction cannot fully resolve the anionic structure, necessitating alternative structure determination methods. We have investigated the anionic structure of anion-disordered cubic (ReO 3 -type) TiOF 2 using X-ray pair distribution function (PDF), 19 F MAS NMR analysis, density functional theory (DFT), cluster expansion modeling, and genetic-algorithm structure prediction. Our computational data predict short-range anion ordering in TiOF 2 , characterized by predominant cis-[O 2 F 4 ] titanium coordination, resulting in correlated anion disorder at longer ranges. To validate our predictions, we generated partially disordered supercells using genetic-algorithm structure prediction and computed simulated X-ray PDF data and 19 F MAS NMR spectra, which we compared directly to experimental data. To construct our simulated 19 F NMR spectra, we derived new transformation functions for mapping calculated magnetic shieldings to predicted magnetic chemical shifts in titanium (oxy)fluorides, obtained by fitting DFT-calculated magnetic shieldings to previously published experimental chemical shift data for TiF 4 . We find good agreement between our simulated and experimental data, which supports our computationally predicted structural model and demonstrates the effectiveness of complementary experimental and computational techniques in resolving anionic structure in anion-disordered oxyfluorides. From additional DFT calculations, we predict that increasing anion disorder makes lithium intercalation more favorable by, on average, up to 2 eV, highlighting the significant effect of variations in short-range order on the intercalation properties of anion-disordered materials.

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

When van der Waals Met Kagome: A 2D Antimonide with a Vanadium-Kagome Network

2D materials showcase unconventional properties emerging from quantum confinement effects. In this work, a “soft chemical” route allows for the deintercalation of K + from the layered antimonide KV 6 Sb 6 , resulting in the discovery of a new metastable 2D-Kagome antimonide K 0.1(1) V 6 Sb 6 with a van der Waals gap of 3.2 Å. The structure of K 0.1(1) V 6 Sb 6 was determined via the synergistic techniques, including X-ray pair distribution function analysis, advanced transmission electron microscopy, and density functional theory calculations. The K 0.1(1) V 6 Sb 6 compound crystallizes in the monoclinic space group C2/m (a = 9.57(2) Å, b = 5.502(8) Å, c = 10.23(2) Å, β = 97.6(2)°, Z = 2). The [V 6 Sb 6 ] layers in K 0.1(1) V 6 Sb 6 are retained upon deintercalation and closely resemble the layers in the parent compound, yet deintercalation results in a relative shift of the adjacent [V 6 Sb 6 ] layers. The magnetic properties of the K 0.1(1) V 6 Sb 6 phase in the 2–300 K range are comparable to those of KV 6 Sb 6 and another Kagome antimonide KV 3 Sb 5 , consistent with nearly temperature-independent paramagnetism. Electronic band structure calculation suggests a nontrivial band topology with flat bands and opening of band crossing afforded by deintercalation. Transport property measurements reveal a metallic nature for K 0.1(1) V 6 Sb 6 and a low thermal conductivity of 0.6 W K –1 m –1 at 300 K. Additionally, ion exchange in KV 6 Sb 6 via a solvothermal route leads to a successful partial exchange of K + with A + (A = Na, Rb, and Cs). Here, this study highlights the tunability of the layered structure of the KV 6 Sb 6 compound, providing a rich playground for the realization of new 2D materials.

2D↗
DOE OSTI2020

A Drift‐Kinetic Method for Obtaining Gradients in Plasma Properties From Single‐Point Distribution Function Data

Abstract In this paper, we derive a new drift‐kinetic method for estimating gradients in the plasma properties through a velocity space distribution at a single point. The gradients are intrinsically related to agyrotropic features of the distribution function. This method predicts the gradients in the magnetized distribution function and can predict gradients of arbitrary moments of the gyrotropic background distribution function. The method allows for estimates on density and pressure gradients on the scale of a Larmor radius, proving to resolve smaller scales than any method currently available to spacecraft. The model is verified with a set of fully kinetic VPIC particle‐in‐cell simulations.

Wetherton, B. A.↗
DOE OSTI2023

Subseasonal Clustering of Atmospheric Rivers Over the Western United States

Abstract The serial occurrence of atmospheric rivers (ARs) along the US West Coast can lead to prolonged and exacerbated hydrologic impacts, threatening flood‐control and water‐supply infrastructure due to soil saturation and diminished recovery time between storms. Here a statistical approach for quantifying subseasonal temporal clustering among extreme events is applied to a 41‐year (1979–2019) wintertime AR catalog across the western United States (US). Observed AR occurrence, compared against a randomly distributed AR timeseries with the same average event density, reveals temporal clustering at a greater‐than‐random rate across the western US with a distinct geographical pattern. Compared to the Pacific Northwest, significant AR clusters over the northern Coastal Range of California and Sierra Nevada are more frequent and occur over longer time periods. Clusters along the California Coastal Range typically persist for 2 weeks, are composed of 4–5 ARs per cluster, and account for over 85% of total AR occurrence. Across the northwest Coast‐Cascade Ranges, clusters account for ∼50% of total AR occurrence, typically last 8–10 days, and contain 3–4 individual AR events. Based on precipitation data from a high‐resolution dynamical downscaling of reanalysis, the fractions of total and extreme hourly precipitation attributable to AR clusters are largest along the northern California coast and in the Sierra Nevada. Interannual variability among clusters highlights their importance for determining whether a particular water year is anomalously wet or dry. The mechanisms behind this unusual clustering are unclear and require further research.

Meteorology & Atmospheric Sciences↗
DOE OSTI2024

A Comprehensive Northern Hemisphere Particle Microphysics Data Set From the Precipitation Imaging Package

Microphysical observations of precipitating particles are critical data sources for numerical weather prediction models and remote sensing retrieval algorithms. However, obtaining coherent data sets of particle microphysics is challenging as they are often unindexed, distributed across disparate institutions, and have not undergone a uniform quality control process. This work introduces a unified, comprehensive Northern Hemisphere particle microphysical data set from the National Aeronautics and Space Administration precipitation imaging package (PIP), accessible in a standardized data format and stored in a centralized, public repository. Data is collected from 10 measurement sites spanning 34° latitude (37°N–71°N) over 10 years (2014–2023), which comprise a set of 1,070,000 precipitating minutes. The provided data set includes measurements of a suite of microphysical attributes for both rain and snow, including distributions of particle size, vertical velocity, and effective density, along with higher-order products including an approximation of volume-weighted equivalent particle densities, liquid equivalent snowfall, and rainfall rate estimates. The data underwent a rigorous standardization and quality assurance process to filter out erroneous observations to produce a self-describing, scalable, and achievable data set. Case study analyses demonstrate the capabilities of the data set in identifying physical processes like precipitation phase-changes at high temporal resolution. Bulk precipitation characteristics from a multi-site intercomparison also highlight distinct microphysical properties unique to each location. This curated PIP data set is a robust database of high-quality particle microphysical observations for constraining future precipitation retrieval algorithms, and offers new insights toward better understanding regional and seasonal differences in bulk precipitation characteristics.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Liquid state properties of SEI components in dimethoxyethane

The solid-electrolyte interphase (SEI) layer is a critical constituent of battery technology, which incorporates the use of lithium metals. Since the formation of the SEI is difficult to avoid, the engineering and harnessing of the SEI are absolutely critical to advancing energy storage. One problem is that much fundamental information about SEI properties is lacking due to the difficulty in probing a chemically complex interfacial system. One such property that is currently unknown is the dissolution of the SEI. This process can have significant effects on the stability of the SEI, which is critical to battery performance but is difficult to probe experimentally. Here, we report the use of ab initio computational chemistry simulations to probe the solution state properties of SEI components LiF, Li 2 O, LiOH, and Li 2 CO 3 in order to study their dissolution and other solution-based characteristics. Ab initio molecular dynamics was used to study the solvation structures of the SEI with a combination of radial distribution functions, discrete solvation structure maps, and vibrational density of states, which allows for the determination of free energies. From the change in free energy of dissolution, we determined that LiOH is the most likely component to dissolve in the electrolyte followed by LiF, Li 2 CO 3 , and Li 2 O although none were favored thermodynamically. Furthermore, this indicates that dissolution is not probable, but Li 2 O would make the most stable SEI with regard to dissolution in the electrolyte.

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

Multidimensional modeling of non-equilibrium plasma generated by a radio-frequency corona discharge

Low-temperature plasma (LTP) ignition concepts rely on the production of radical and charged species to speed up the onset of combustion in spark-ignition engines. These features are responsible for the superior performance of LTP igniters under extremely dilute combustion operation that is not achievable by conventional spark igniters. Additionally, LTP discharges extend the lifetime of the igniters, due to the avoidance of spark processes. For these reasons, the engine research community and the automotive industry have shown growing interest in this technology in the recent years. As of today, computational fluid-dynamics (CFD) codes typically used by the multi-dimensional engine modeling community do not have reliable models to describe LTP ignition processes. One key missing piece of information is the physical and chemical properties of the plasma and their effect on combustion ignition. Most non-equilibrium plasma simulations reported in literature are based on simplified, canonical geometries, with simple discharge excitation schemes. Here we conduct multi-dimensional modeling of the non-equilibrium plasma generated by an application-relevant radio-frequency (RF) corona discharge in air. Three test cases are simulated, characterized by different environmental pressure levels and peak electrode voltage values at room temperature. Streamer penetration, electron number density, atomic oxygen production, and bulk gas temperature distribution in the first 10 sinusoidal pulses are presented and discussed. This model can be used as a key tool for an in-depth understanding of RF-corona discharge for automotive applications and provides the basis for future implementations of dedicated LTP ignition models in CFD codes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

The effect of saturated thermal conduction on clouds in a hot plasma

We numerically investigate the internal evolution of multiphase clouds, which are at rest with respect to an ambient, highly ionized medium (HIM) representing the hot component of the circumgalactic medium. Time-dependent saturated thermal conduction and its implications like condensation rates and mixing efficiency are assessed in multiphase clouds. Our simulations are carried out by using the adaptive mesh refinement code FLASH. The model clouds are initially in both hydrostatic and thermal equilibrium and are in pressure balance with the HIM. Thus, they have steep gradients in both temperature and density at the interface to HIM leading to non-negligible thermal conduction. Several physical processes are considered numerically or semi-analytically: thermal conduction, radiative cooling and external heating of gas, self-gravity, mass diffusion, and dissociation of molecules and ionization of atoms. It turns out that saturated thermal conduction triggers a continuous condensation irrespective of cloud mass. Dynamical interactions with ambient HIM all relate to the radial density gradient in the clouds: (1) mass flux due to condensation is the higher the more homogeneous the clouds are; (2) mixing of condensed gas with cloud gas is easier in low-mass clouds, because of their shallower radial density gradient; and thus (3) accreted gas is distributed more efficiently. A distinct and sub-structured transition zone forms at the interface between cloud and HIM, which starts at smaller radii and is much narrower as deduced from analytical theory.

79 ASTRONOMY AND ASTROPHYSICS↗