Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “density evolution”

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 289 records · Page 16

Enhanced Activity in Layered Metal-Oxide-Based Oxygen Evolution Catalysts by Layer-by-Layer Modulation of Metal-Ion Identity

Few-layered potassium nickel and cobalt oxides show drastic differences in catalytic activity based on metal ion preorganization. Uniform compositions [(CoO 2 /K) 6 or (NiO 2 /K) 6 ] show limited activity, while homogeneously mixed-metal cobalt/nickel oxides [(Co n Ni (1–n) O 2 /K) 6 ] display moderate improvement. However, a layer-by-layer arrangement of alternating cobalt and nickel oxide sheets [e.g., (CoO 2 /K/NiO 2 /K)] provides superior catalytic performance, reducing the oxygen evolution overpotential by ∼200–400 mV. Density functional theory simulations provide an illustration of the electronic properties (density of states and localization of orbitals) that promote catalysis in the layer-segregated materials over those of homogeneous composition. This study reveals that atomic preorganization of metal ions within layered catalysts plays a more crucial role than the overall metal composition in enhancing catalytic efficiency for oxygen evolution.

catalysts↗

A theoretical study of the lifetime and transport of large ionospheric density structures

A three-dimensional time-dependent ionospheric model was used to study the spatial and temporal evolution and transport of large-scale high-density ionospheric structures for a range of solar cycle, seasonal, and IMF conditions. Both density depletions and enhancements were considered. It was found that, depending on the IMF, horizontal plasma convection can cause an initial structure to break up into multiple structures of various sizes, to become stretched into elongated segments, or to remain as a single distorted structure. The lifetime of an F-region density structure depends on several factors, including its magnitude, the initial location where it was formed, the season, the solar cycle, and the convection pattern. For example, in summer, the effects of a large density structure can disappear in a few hours or last as long as nine hours, while in winter the effects can persist for 24 hours. The passage of perturbed plasma flux tubes through sunlit and auroral regions can significantly increase the lifetime of plasma enhancements.

Schunk, R. W.↗

Phase-field modeling of stored-energy-driven grain growth with intra-granular variation in dislocation density

Abstract We present a phase-field (PF) model to simulate the microstructure evolution occurring in polycrystalline materials with a variation in the intra-granular dislocation density. The model accounts for two mechanisms that lead to the grain boundary migration: the driving force due to capillarity and that due to the stored energy arising from a spatially varying dislocation density. In addition to the order parameters that distinguish regions occupied by different grains, we introduce dislocation density fields that describe spatial variation of the dislocation density. We assume that the dislocation density decays as a function of the distance the grain boundary has migrated. To demonstrate and parameterize the model, we simulate microstructure evolution in two dimensions, for which the initial microstructure is based on real-time experimental data. Additionally, we applied the model to study the effect of a cyclic heat treatment (CHT) on the microstructure evolution. Specifically, we simulated stored-energy-driven grain growth during three thermal cycles, as well as grain growth without stored energy that serves as a baseline for comparison. We showed that the microstructure evolution proceeded much faster when the stored energy was considered. A non-self-similar evolution was observed in this case, while a nearly self-similar evolution was found when the microstructure evolution is driven solely by capillarity. These results suggest a possible mechanism for the initiation of abnormal grain growth during CHT. Finally, we demonstrate an integrated experimental-computational workflow that utilizes the experimental measurements to inform the PF model and its parameterization, which provides a foundation for the development of future simulation tools capable of quantitative prediction of microstructure evolution during non-isothermal heat treatment.

Materials Science↗

Particle tracing in the magnetosphere - New algorithms and results

A fast and efficient method is employed to trace charged particles through realistic magnetospheric electric and magnetic fields, greatly reducing computer simulation times. The method is applicable for particles having arbitrary charge, energy, or pitch angle but which conserve the first two adiabatic invariants. An efficient method is used to classify drift orbits, which greatly facilitates identification of the last closed drift path or other drift boundaries. The time-independent evolution of the bounce-averaged phase space density along convective drift orbits is calculated. These three tools can be used to quantitatively describe convective evolution of the particle distribution from the tail, an essential step in understanding the production of unstable distributions in the magnetosphere.

Sheldon, R. B.↗

A three-dimensional model of corotating streams in the solar wind. II - Hydrodynamic streams

Theoretical aspects of corotating solar wind dynamics on a global scale are examined by means of numerical simulations executed with a nonlinear, inviscid, adiabatic, single fluid, three-dimensional hydrodynamic formulation. The factors most strongly affecting stream evolution, namely, the intrinsic correlations among density, temperature, and velocity existing near the source, the amplitude, longitudinal breadth and latitudinal breadth of the stream, and the heliographic latitude of the centroid of the stream, are considered. Calculations indicate that even in the presence of substantial meridional gradients, the latitudinal transport of mass, energy, and momentum by streams appears to have little impact on the overall evolution of the structure.

Pizzo, V. J.↗

Elucidating Gas Reduction Effects of Organosilicon Additives in Lithium-Ion Batteries

Lithium-ion batteries (LIBs) with nonaqueous liquid electrolytes are prone to gas generation at elevated voltages and temperatures, degrading battery performance and posing serious safety risks. Organosilicon (OS) additives are an emerging candidate solution for gassing problems in LIBs, but a detailed understanding of their functional mechanisms remains elusive. In this work, we present a combined computational and experimental study to elucidate the gasreducing effects of OS additives. Cell volume measurements and gas chromatography–mass spectrometry reveal that OS additives can substantially reduce gas evolution in LIBs, particularly CO 2 regardless of source. Through density functional theory calculations, we identify multiple plausible pathways for CO 2 evolution, including (1) nucleophile-induced ring-opening of ethylene carbonate (EC) and the subsequent electro-oxidation and (2) direct electro-oxidation of lithium carbonate (Li 2 CO 3 ). Correspondingly, we find that OS additives function via two primary mechanisms: (1) scavenging of nucleophiles such as superoxide (O 2 •– ), peroxide (O 2 2– ), and carbonate ion (CO 3 2– ); (2) oligomerization with ethylene carbonate oxide ion and ethylene dicarbonate ion. Moreover, we discover that OS additives possess strong lithium coordination affinity, which helps further reduce the nucleophilic reaction energies and hence increases their nucleophile-scavenging efficiency. Finally, we provide a mechanistic interpretation for the enhanced gasreduction effects observed with fluorinated OS compounds, corroborated by surface analysis results from X-ray photoelectron spectroscopy. Our study offers the first molecular-level insights into how OS additives contribute to reduced gas formation in LIBs, paving the way for improved safety and performance of LIBs.

25 ENERGY STORAGE↗

Unveiling the Electrocatalytic Hydrogen Evolution Reaction Pathway on RuP 2 through Ab Initio Grand Canonical Monte Carlo

In this study, the high catalytic reactivity of ruthenium phosphide (RuP 2 ) has been identified by first-principles density functional theory (DFT) calculations for the electrocatalytic hydrogen evolution reaction (HER). Complex surface reconstructions are considered by applying the ab initio grand canonical Monte Carlo (ai-GCMC) algorithm, efficiently providing a sufficient phase-space exploration of possible surfaces. Combined with surface-phase Pourbaix diagrams, we are able to identify the actual surfaces that obtained under specific experimental environments, thus leading to a more accurate understanding of the nature of the active sites and the binding strength of adsorbates. Specifically, through hundreds of surface reconstructions and hydrogenation states generated with ai-GCMC, we identify the most favorable surface phases of RuP 2 under aqueous acidic conditions. We discover that the HER activity is determined by multiple surfaces with different stoichiometries within a narrow electrode potential window. Low HER overpotential (η) has been found for each of the identified surfaces, as low as 0.04 V. High H-coverage reconstructed surfaces have been discovered under acidic conditions, and the surface Ru sites introduced by additional Ru adatoms or exposed by P-vacancies serve as the active sites for HER based on their nearly reversible H binding. Furthermore, this work provides atomistic insights into the origin of high HER activity on RuP 2 by exploring the dynamic surface phases of electrocatalysts and features a generalizable method to explore the reconstructed/hydrogenated surface space as a function of experimental conditions.

25 ENERGY STORAGE↗

Testing the parametric model for self-interacting dark matter using matched halos in cosmological simulations

Here, we systemically evaluate the performance of the self-interacting dark matter (SIDM) halo model proposed in Ref.Yang et al. (2023) with matched halos from high-resolution cosmological CDM and SIDM simulations. The model incorporates SIDM effects along mass evolution histories of CDM halos and it is applicable to both isolated halos and subhalos. We focus on the accuracy of the model in predicting halo density profiles at z = 0 and the evolution of maximum circular velocity. We find the model predictions agree with the simulations within 10%–50% for most of the simulated (sub)halos, 50%–100% for extreme cases. This indicates that the model effectively captures the gravothermal evolution of the halos with very strong, velocity-dependent self-interactions. For an example application, we apply the model to study the impact of various SIDM scenarios on strong lensing perturber systems, demonstrating its utility in predicting SIDM effects for small-scale structure analyses. Our findings confirm that the model is an effective tool for mapping CDM halos into their SIDM counterparts.

Yang, Daneng [Chinese Academy of Sciences, Nanjing↗

Second mode interactions in supersonic boundary layers

The nonlinear evolution of a two-dimensional second mode unstable wave in a Mach 4.5 wall-bounded flow is computed by solving the full time-dependent compressible Navier-Stokes equations. A highly accurate solution is obtained using spectral collocation methods. It is shown that departure from linearity first occurs in the critical layer due to the cubic nonlinearities in the momentum equation. This is a direct result of the large density perturbations in this regime. Time evolution studies of the growth rate as a function of normal distance from the plate suggests that the mode is evolving toward a nonlinear saturated state, and that this problem is possibly amenable to standard weakly nonlinear perturbation methods.

Erlebacher, Gordon↗

Early time evolution of a chemically produced electron depletion

The early time evolution of an ionospheric electron depletion produced by a radially expanding electron attachment chemical release is studied with a two-dimensional simulation model. The model includes electron attachment chemistry, incorporates fluid electrons, particle ions and neutrals, and considers the evolution in a plane perpendicular to the geomagnetic field for a low beta plasma. Timescales considered are of the order of or less than the cyclotron period of the negative ions that result as a by-product of the electron attacment reaction. This corresponds to time periods of tenths of seconds during recent experiemts. Simulation results show that a highly sheared azimuthal electron flow velocity develops in the radially expanding depletion boundary. This sheared electron flow velocity and the steep density gradients in the boundary give rise to small-scale irregulatities in the form of electron density cavities and spikes. The nonlinear evolution of these irregularities results in trapping and ultimately turbulent heating of the negative ions.

Scales, W. A.↗

The Cosmic Evolution of C IV Absorbers at 1.4 < z < 4.5: Insights from 100,000 Systems in DESI Quasars

We present the largest catalog to date of triply ionized carbon (C IV ) absorbers detected in quasar spectra from the Dark Energy Spectroscopic Instrument. Using an automated matched-kernel convolution method with adaptive signal-to-noise thresholds, we identify 101,487 C IV systems in the redshift range 1.4 < z < 4.5 from 300,637 quasar spectra. Completeness is estimated via Monte Carlo simulations, and the catalog is 50% complete at EW C IV ≥ 0.4 Å. The differential equivalent width frequency distribution declines exponentially and shows weak redshift evolution. The absorber incidence per unit comoving path increases by a factor of 2–5 from z ≈ 4.5 to z ≈ 1.4, with stronger redshift evolution for strong systems. Using column densities derived from the apparent optical depth method, we constrain the cosmic mass density of C IV , Ω C IV , which increases by a factor of ∼3.8 from (0.82 ± 0.05) × 10 −8 at z ≈ 4.5 to (3.16 ± 0.2) × 10 −8 at z ≈ 1.4. From Ω C IV , we estimate a lower limit on intergalactic medium metallicity ${\mathrm{log}}({Z}_{{\rm{IGM}}}/{Z}_{\odot })\gtrsim -3.25$ at z ∼ 2.3, with a smooth decline at higher redshifts. These trends trace the cosmic star formation history and He II photoheating rate, suggesting a link between C IV enrichment, star formation, and UV background over ∼3 Gyr. The catalog also provides a critical resource for future studies connecting circumgalactic metals to galaxy evolution, especially near cosmic noon.

79 ASTRONOMY AND ASTROPHYSICS↗

Collapse and relaxation of rotating stellar systems

The evolution of initially spherical, uniform-density, and uniform-velocity-dispersion stellar systems is investigated by means of a three-dimensional N-body computer model. Some 100,000 simulation stars are used to follow the collapse and relaxation of the system for various amounts of angular momentum in solid-body rotation. For initially low values of the angular momentum satisfying the Ostriker-Pebbles stability criterion the systems quickly relax to an axisymmetric shape and resemble elliptical galaxies in appearance. The maximum flattening for these systems is equivalent only to an E2 system. For larger values of the initial angular momentum bars develop, and the systems undergo a much more drastic evolution. The apparent rotational and random velocities of the barred systems are very sensitive to the viewing direction. An additional complication is the frequent misalignment of the apparent major axis with the direction that reflects the maximum rotation.

Hohl, F.↗

Experimental and simulation study of target biasing effects on plasma transport in linear plasma device MPS-LD

Linear plasma devices (LPDs) are important experimental platforms for investigating plasma–material interactions (PMI). In PMI experiments, it has been found that applying a target bias not only effectively modifies the incident ion energy, but also induces significant changes in the electron density and electron temperature, whereby the evolution of these plasma parameters is primarily governed by plasma transport processes. However, at present, the physical process and mechanism underlying such bias-induced variations remain unclear. In this work, biasing experiments under argon plasma discharge conditions were first carried out on the MPS-LD device. For the corresponding experiments, an electric potential model was newly developed based on the BOUT++ LPD module, enabling self-consistent simulations of plasma transport under biased conditions. Numerical simulations were then performed to reproduce the experimental results and to validate the accuracy of the proposed model. Finally, by combining experimental measurements with numerical simulations, a bias-voltage scan was performed to investigate how the electron density and electron temperature vary with the bias voltage (U bias ). The results show that applying negative bias decreases the target electron density (n e,T ) while increasing the target electron temperature (T e,T ). In contrast, positive bias increases both n e,T and T e,T ; however, at high positive bias, n e,T first reaches a maximum and subsequently decreases with further increases in U bias . The underlying physical mechanisms are analyzed using particle flux, momentum, and energy conservation. It indicates that the applied bias regulates the parallel electric field, thereby changing ion and electron velocities, and consequently affecting the electron density. At high positive bias, the ion velocity is further influenced by ion viscosity, leading to the reversal in n e,T . Meanwhile, the enhanced parallel electric field drives stronger currents, significantly increasing ion–electron frictional work and converting the input bias power into electron energy, which raises the electron temperature. In conclusion, these results contribute to a deeper understanding of the effects and mechanisms of biasing on plasma transport in the MPS-LD device.

BOUT++ simulation↗

Numerical calculations of the collapse of nonrotating, magnetic gas clouds

Results of the first self-consistent numerical calculations of the dynamic collapse of a magnetized protostellar gas cloud are presented. Symmetry about an axis parallel to the initial magnetic field direction has been assumed, so that the calculations could be performed on a two-dimensional grid. Also, the cloud was taken to be nonrotating and isothermal, and the magnetic field was assumed to remain frozen in to the gas. As starting models for the calculations, gas spheres with uniform density and magnetic field were used. The time evolution of the clouds has been calculated for roughly two initial free-fall times, at which point the central density has increased by a factor of approximately 10,000 to 1,000,000. Several such calculations have been performed for different values of the cloud's initial thermal, magnetic, and gravitational energies. In virtually all cases it is found that, once a flattened core forms in the cloud, the central magnetic field strength, B, varies with gas density, rho, according to (d log B/d log rho) = 1/2. This behavior is independent of the initial energy ratios mentioned above. It is also found that the magnetic field is able to prevent completely the collapse of part of the outer envelope of the cloud.

Scott, E. H.↗

Shock propagation in aerogel and TPP foams for inertial fusion energy target design

Achieving practical inertial fusion energy (IFE) requires the development of target designs with well-characterized microstructure and compression response. We measured shock dynamics in low-density (17.5–500 mg/cm 3 ) aerogel and two-photon polymerization (TPP) foams using x-ray phase contrast imaging (XPCI) methods and the Velocity Interferometer System for Any Reflector. By analyzing shock front evolution, we examined how target type and density influence shock propagation and energy dissipation. Talbot-XPCI shows that aerogels support a smooth, bowed shock front due to their homogeneous nanometer-scale pore network. In contrast, TPP foams exhibit irregular, stepwise propagation driven by interactions with their periodic micrometer-scale lattice. Shock velocity follows a power-law relation: aerogels deviate from classical ρ −1/2 scaling due to pore-collapse dissipation, while TPP foams follow the trend with larger uncertainties from density variations. Comparisons with xRAGE simulations reveal systematic underestimation of shock speeds. These results provide the first experimental constraints on shock propagation in TPP foams over a wide density range and highlight the influence of internal structure on anisotropic shock behavior. Our findings support improved benchmarking of EOS and hydrodynamic models and inform the design of foam architectures that promote implosion symmetry in IFE capsules.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Coupled neutronics, thermochemistry, corrosion modeling and sensitivity analyses for isotopic evolution in molten salt reactors

This study presents a computational methodology for analyzing isotopic evolution and associated uncertainties in molten salt reactors (MSRs), focusing on both fluoride- and chloride-based fuel salts. The primary goal is to enhance the understanding of isotopic behavior in MSRs and provide data to support future experimental efforts. The methodology integrates transport-coupled depletion calculations using OpenMC, equilibrium thermodynamics modeling with Thermochimica, and a corrosion model. Sensitivity analyses are performed to evaluate the impact of power density, air ingress, and humidity content on isotopic evolution in MSR concepts. This study examines representative F- and Cl-based MSR designs, highlighting the dominant influence of power density on isotopic composition, which significantly affects isotope production and depletion rates, accounting for approximately 76% of the observed variance in element concentration. Air ingress and humidity content also affect the redox potential, solubility of heavier elements, and corrosion rates, thereby altering the expected isotopic evolution in the reactor. On average, air ingress accounts for around 17% of the variance in element concentrations, while humidity explains the remaining 7%. These variances differ significantly from element to element, depending on the element’s role in depletion, redox potential evolution, and galvanic corrosion. The findings indicate that power density, air ingress, and humidity content are all critical factors for optimizing reactor design and operational strategies. Furthermore, the study provides expected ranges for key impurities in the fuel salt, which are crucial for guiding future experimental studies and refining MSR designs. Finally, this study demonstrates the importance of modeling depletion coupled with the evolution of redox potential and chemical interactions in MSR fuel salts.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

The evolution of cocoons surrounding light, extragalactic jets

If the mass density of supersonic, collimated material is less than that of the surrounding medium, a so-called light jet will be enveloped by a cocoon of overpressured shocked gas. Hydrodynamical simulations are used to understand the evolution of the cocoon. The cocoon's evolution is also compared to a simple analytic theory. To reconcile the theory with the simulations, the growth of the jet head must be taken into account. The overpressured cocoon stage exists for a relatively short astronomical time, after which only the region of the cocoon near the jet head remains overpressured. The spatial distribution of the optical emission often observed in distant extragalactic jet systems can be explained with this improved understanding of cocoon evolution.

Cioffi, Denis F.↗