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

Comparison of laser-produced plasma spatio-temporal electron density evolution measured using interferometry with simulation results

Due to the difficulties associated with experimental measurements of laser-produced plasma (LPP) properties during the earliest stages of plasma evolution, radiation hydrodynamic codes are often used. However, although these codes have been extensively validated in the higher intensity regimes, validation at low to moderate intensities has been limited. In this study, the spatio-temporal electron density evolution in an LPP generated at moderate laser intensities and at various laser wavelengths was validated against the FLASH code for times up to 20 ns. The LPP was generated by focusing the fundamental and various harmonics radiation (1064, 532, and 266 nm) from a 6 ns full width half maximum Nd:YAG laser, at a laser intensity of 10 GW cm -2 , onto a copper target. The spatio-temporal density evolution of the expanding plasma was analyzed using Nomarski interferometry. Experimental measurements were found to be consistent with FLASH simulations, and the dependence of electron densities on wavelength was found to be in agreement with analytical models, varying as n e ∝ λ -0.7 . However, slight differences were noted in the widths and shape functions of the experimental and simulated electron density profiles.

2D fast imaging↗
DOE OSTI2026

Time-resolved interferometric measurements of plasma density evolution in laser-driven capacitor-coil targets

Laser-driven capacitor-coil targets provide a compact platform for generating strong magnetic fields and are widely used in magnetized high-energy-density (HED) plasma experiments. In addition to magnetic-field generation, these targets also produce plasma in the coil region, which can influence the subject physical processes and interact with secondary targets or external plasmas in their applications. However, direct, time-resolved measurements of the plasma density surrounding the coil remain limited. Here, we report interferometric measurements of the plasma density evolution in laser-driven capacitor-coil targets irradiated by the University of Osaka LFEX laser. Two-dimensional electron density maps reveal two distinct plasma sources loading the coil region: plasma generated in the coil itself and plasma produced by laser ablation of the target plates. These results provide quantitative information on plasma loading and evolution in capacitor-coil targets and are directly relevant to the design and modeling of magnetized HED plasma experiments.

Zhang, Yang [Princeton University, NJ (United Stat↗
DOE OSTI2025

Rules for the crystallite size and dislocation density evolution in phases during α-ω transformation in Zr under high-pressure and severe plastic flow

The first in-situ X-ray diffraction (XRD) study of the evolution of the crystallite size and dislocation density in phases during plastic strain-induced phase transformation (PT) is performed utilizing α-ω PT in strongly pre-deformed commercially pure Zr as an important example. Rough diamond anvils (rough-DA) are introduced to intensify all occurring processes during heterogeneous compression of Zr in a diamond anvil cell (DAC). The main rule is found that during α-ω PT the crystallite size and dislocation density in ω-Zr depend solely on the volume fraction of ω-Zr and are independent of pressure, plastic strain tensor, its path, and initial nanostructure. Crystallite size in ω-Zr increases from 10 to 60 nm during the PT, while dislocation density reduces from 1.83×10 15 /m 2 to 0.6×10 15 /m 2 . Rough-DA produce a steady nanostructure in α-Zr before PT with smaller crystallite size and larger dislocation density than smooth-DA, leading to a reduction of the minimum pressure for α-ω PT to a record value 0.67 GPa, 9 times smaller than under hydrostatic loading and 5.1 times lower than the phase equilibrium pressure. In addition to strain, the kinetics of strain-induced PT unexpectedly depends on time. Also, strain-controlled part of kinetics is zero order, in contrast to the first-order kinetics with smooth-DA. The obtained results open a new window for understanding the mutual effects of nanostructure evolution and PT during severe plastic flow in various technological and natural processes. In conclusion, they may bring up economic strategies of producing nanocomposites and single-phase nanostructured materials with optimal properties.

Crystalline size↗
DOE OSTI2025

Rapid synthesis of phase-engineered tungsten carbide electrocatalysts via flash joule heating for high-current-density hydrogen evolution

Fabricating durable and high-performance electrocatalysts operating at high current densities for industrial acidic hydrogen evolution remains a daunting challenge. Tailoring the phase composition of electrocatalysts is a promising strategy to harness synergistic effects and improve charge transfer, thereby optimizing their performance. Here, this work presents a fast, green method based on flash joule heating (FJH) to synthesize phase-engineered tungsten carbide electrocatalysts for the acidic hydrogen evolution reaction (HER) at high current densities. Tungsten carbide electrodes with varying FJH treatment durations (3, 10, 30, and 60 s) are fabricated to fine-tune the mixture of tungsten monocarbide (WC) and tungsten semicarbide (W 2 C) phases. Results show that samples with a 30-s treatment (TC-3) exhibit an optimal balance between these phases, leading to a low overpotential of 180.97 and 387 mV at current densities of 10 mA/cm 2 and 4 A/cm 2 , respectively. TC-3 exhibits significantly lower charge transfer resistance compared to the other electrocatalysts, which can be attributed to its optimal phase ratio. Notably, the TC-3 electrocatalyst remains stable for over 9 days at 4 A/cm 2 due to their controlled phases and excellent corrosion-resistant properties. This work highlights a new method to fabricate cost-effective, high-performance tungsten carbide electrocatalysts with well-controlled phase compositions.

Acidic hydrogen evolution reaction↗
DOE OSTI2024

A FFT-based mesoscale continuum dislocation mechanics with defect energy: Applications to composites and polycrystals

A crystal plasticity elastoviscoplastic FFT (fast Fourier transform) formulation with a mesoscale continuum field dislocation mechanics model is presented, which incorporates a defect energy density that depends on GND densities and an associated material length scale. This allows to thermodynamically derive internal length scale dependent intra-crystalline backstress and Peach–Koehler force acting on GND densities. The model considers GND density evolution through a filtered numerical spectral approach, which is coupled with stress equilibrium through the elastoviscoplastic FFT algorithm. The discrete Fourier transform (DFT) method together with finite difference (FD) schemes is applied to solve both the backstress tensor and the Fourier–Green operator. Numerical results are first reported for two-phase laminate composites with plastic single crystal channels and elastic precipitates for shear loadings. Channel size effects are simulated and analyzed on the overall and local hardening behaviors during monotonous loadings. In addition, the evolutions of GND densities and the role of their associated backstress on size effects are examined during reversible shear loading. In a second part, the role of the defect energy internal length scale on polycrystal’s hardening during tension–compression is discussed. The results are compared to those obtained using FFT-based continuum field dislocation mechanics without defect energy.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Luminosity Function of Quasars at 1.0 < z < 3.5 from SDSS and DESI

We present a study of the evolution of type 1 quasars at 1.0 < z < 3.5, covering the peak epoch of quasar activity. The quasar evolution has been extensively explored by a variety of previous works and the derived quasar luminosity functions (QLFs) are not very consistent with each other, presumably due to the complexities introduced by different quasar selection techniques and associated completeness corrections. We use a new strategy to construct QLFs based on a library of all known quasars. We focus on a wide region of ∼1700 deg 2 and a deep field of ∼265 deg 2 that have rich spectroscopic data primarily from the Sloan Digital Sky Survey and DESI. We then apply traditional color cuts in the rest-frame UV–optical to select quasar candidates and use the quasar library to identify them. Our final sample consists of 62,426 quasars at 1.0 < z < 3.5, with a high completeness (∼96%) and a high purity (∼93%) in the color selection. Simple color cuts can potentially minimize selection biases for the study of quasar evolution. We derive binned QLFs and characterize them using a double power-law model. Sample incompleteness and contamination are considered as part of the uncertainties in the calculation. Compared to previous results, our QLFs are slightly higher at the faint end, and also higher at the bright end at 2.5 < z < 3.5. The QLFs suggest that the quasar evolution at 1.0 < z < 2.5 can be well described by the pure luminosity evolution model, while at 2.5 < z < 3.5, it can be described by either the pure luminosity evolution or the pure density evolution model.

Luminosity function↗
DOE OSTI2024

A magnetic analog of pressure–strain interaction

We study the evolution equation for magnetic energy density for a non-relativistic magnetized plasma in the (Lagrangian) reference frame comoving with the electron bulk velocity. Analyzing the terms that arise due to the ideal electric field, namely, perpendicular electron compression and magnetic field line bending, we recast them to reveal a quantity with a functional form analogous to the often-studied pressure–strain interaction term that describes one piece of internal energy density evolution of the species in a plasma, except with the species pressure tensor replaced by the magnetic stress tensor. We dub it the “magnetic stress–strain interaction.” We discuss decompositions of the magnetic stress–strain interaction analogous to those used for pressure–strain interaction. These analogies facilitate the interpretation of the evolution of the various forms of energy in magnetized plasmas and should be useful for a wide array of applications, including magnetic reconnection, turbulence, collisionless shocks, and wave–particle interactions. We display and analyze all the terms that can change magnetic energy density in the Lagrangian reference frame of the electrons using a particle-in-cell simulation of magnetic reconnection.

Barbhuiya, M. Hasan (ORCID:0000000163301650)↗
DOE OSTI2026

Velocity-space Origins of the Pressure–Strain Interaction in Multipopulation Distributions and Its Application to Magnetic Reconnection

A forefront research question is how energy evolves in weakly collisional plasmas for which departures from local thermodynamic equilibrium (LTE) are significant. The standard approach is studying the terms in the non-LTE energy evolution equation derived by taking the second moment of the Boltzmann equation, but the resultant fluid metrics do not retain information about which particles at which velocities drive energy evolution. A widely studied channel for internal energy density evolution is the pressure–strain interaction. Here, we employ the kinetic pressure–strain, a phase-space diagnostic whose velocity-space integral recovers the pressure–strain interaction to disambiguate the contributions to the pressure–strain interaction from disparate particle populations in composite phase-space densities. We develop phase-space analogs of the pressure–strain interaction decompositions to provide the phase-space origins of normal versus sheared flow. We introduce the “kinetic strain-rate” tensor, the phase-space analog of the strain-rate tensor, which we argue is needed to interpret the phase-space origins of the pressure–strain interaction. To demonstrate the utility of these quantities, we investigate them for composite electron distributions near the electron diffusion region in two-dimensional particle-in-cell simulations of antiparallel symmetric magnetic reconnection. We find that the phase-space-based diagnostics isolate the roles of distinct populations. These results contribute to a growing body of work providing new methods for quantifying phase-space energy evolution for a broad array of processes, from magnetic reconnection to collisionless shocks and turbulence, opening new pathways for answering longstanding problems of particle energization in weakly collisional plasmas.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Accelerating charge estimation in molecular dynamics simulations using physics-informed neural networks: corrosion applications

Molecular Dynamics (MD) simulations are used to understand the effects of corrosion on metallic materials in salt brine. Reactive force fields in classical MD enable accurate modeling of bond formation and breakage in the aqueous medium and at the metal-electrolyte interface, while also facilitating dynamic partial charge equilibration. However, MD simulations are computationally intensive and unsuitable for modeling the long time scales characteristic of corrosive phenomena. To address this, we develop reduced-order machine learning models that provide accurate and efficient predictions of charge density in corrosive environments. Specifically, we use Long Short-Term Memory (LSTM) networks to forecast charge density evolution based on atomic environments represented by Smooth Overlap of Atomic Positions (SOAP) descriptors. A physics-informed loss function enforces charge neutrality and electronegativity equivalence. The atomic charges predicted by the deep learning model trained on this work were obtained two orders of magnitude faster than those from molecular dynamics (MD) simulations, with an error of less than 3% compared to the MD-obtained charges, even in extrapolative scenarios, while adhering to physical constraints. This demonstrates the excellent accuracy, computational efficiency, and validity of the developed model. Lastly, even though developed for corrosion, these protocols are formulated in a phenomenon-agnostic manner, allowing application to various variable-charge interatomic potentials and related fields.

Atomistic models↗
DOE OSTI2026

Evolution of charge density wave soft phonon modes in Pd 𝑥 ⁢ErTe 3

Here, we investigated the lattice dynamics of quasi-two-dimensional Pd-intercalated ErTe 3 in relation to its charge-density-wave (CDW) transitions by means of x-ray diffuse and meV-resolution inelastic x-ray scattering. In pristine ErTe 3 , CDW order develops at orthogonal in-plane wave vectors q$^𝑐_1$ =(0,0,0.29) (the 𝑐−CDW ) and q$^𝑎_2 $= (0.31,0,0) (the 𝑎−CDW ), with transition temperatures $𝑇^𝑐_1$ = 270K and $𝑇^𝑎_2$ = 160K, respectively. Remarkably, we observe diffuse x-ray scattering already near the higher transition temperature $𝑇^𝑐_1$ along 𝑎−CDW but at a slightly different wave vector q$^𝑎_1$ = (0.29,0,0). Inelastic x-ray scattering for Pd 0.01 ⁢ErTe 3 shows that a partial phonon softening at q$^𝑎_1$, underscoring the strong competition between ordering tendencies along the nearly equivalent in-plane axes of the orthorhombic lattice. For intercalation levels 𝑥 ≥ 0.02, the 𝑎−CDW state is suppressed. Nevertheless, a similar correlation between phonon softening and diffuse scattering persists along the [100] direction, again observed at q$^𝑎_1$ = (0.29,0,0) and $𝑇^𝑐_1$. These findings suggest that the 𝑎−CDW is fully suppressed for 𝑥 ≥ 0.02, and that the residual diffuse scattering at originates from the partial phonon softening associated with the 𝑐−CDW, reflected by the near equality of the absolute size of q$^𝑐_1$ and q$^𝑎_1$.

Maity, Avishek [Technische Univ. München, Garching↗
DOE OSTI2024

Determining the N -Representability of a Reduced Density Matrix via Unitary Evolution and Stochastic Sampling

The N-representability problem consists in determining whether, for a given p-body matrix, there exists at least one N-body density matrix from which the p-body matrix can be obtained by contraction, that is, if the given matrix is a p-body reduced density matrix (p-RDM). The knowledge of all necessary and sufficient conditions for a p-body matrix to be N-representable allows the constrained minimization of a many-body Hamiltonian expectation value with respect to the p-body density matrix and, thus, the determination of its exact ground state. However, the number of constraints that complete the N-representability conditions grows exponentially with system size, and hence, the procedure quickly becomes intractable for practical applications. This work introduces a hybrid quantum-stochastic algorithm to effectively replace the N-representability conditions. The algorithm consists of applying to an initial N-body density matrix a sequence of unitary evolution operators constructed from a stochastic process that successively approaches the reduced state of the density matrix on a p-body subsystem, represented by a p-RDM, to a target p-body matrix, potentially a p-RDM. The generators of the evolution operators follow the well-known adaptive derivative-assembled pseudo-Trotter method (ADAPT), while the stochastic component is implemented by using a simulated annealing process. The resulting algorithm is independent of any underlying Hamiltonian, and it can be used to decide whether a given p-body matrix is N-representable, establishing a criterion to determine its quality and correcting it. We apply the proposed hybrid ADAPT algorithm to alleged reduced density matrices from a quantum chemistry electronic Hamiltonian, from the reduced Bardeen–Cooper–Schrieffer model with constant pairing, and from the Heisenberg XXZ spin model. In all cases, the proposed method behaves as expected for 1-RDMs and 2-RDMs, evolving the initial matrices toward different targets.

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

Predicting the evolution of biomass bulk density through feedstock preprocessing: Discrete element modeling, regression analysis, and pilot-scale validation

Bulk density is an important material property of biomass feedstocks, influencing handling, storage, transport costs, and conversion efficiency. In this study, predictive regression models for loose and tapped bulk densities of Alamo and Cave-in-Rock switchgrass are developed using a comprehensive dataset generated via calibrated bonded-sphere discrete element method (DEM) simulations. Here, a key contribution of this study is the use of a DEM-based approach, which correlates density with moisture content and particle size distribution parameters and enables analysis across a continuous particle size range, overcoming limitations of purely experimental data. For comparison, regression models are also developed using only experimental data from pilot-scale runs at the Biomass Feedstock National User Facility at Idaho National Laboratory. Validation against pilot-scale data showed reasonable prediction accuracy for both model types, particularly for smaller particle sizes (post-secondary grinding). While the experimental model showed slightly better performance matching the validation data in some cases, the DEM-based model benefits from a much larger dataset, reduced predictor multicollinearity, and continuous parameter coverage, highlighting the utility of validated simulation models for developing robust predictive tools for biomass preprocessing applications.

09 - BIOMASS FUELS↗
DOE OSTI2025

Measurement of interfacial thermal resistance in high-energy-density matter

Heat transport across interfaces is a ubiquitous phenomenon with many unresolved aspects. In particular, it is unknown if an interfacial thermal resistance (ITR) occurs in matter with high-energy-density where free electrons dominate the heat conduction. Here, we report on the first experimental evidence that a significant heat barrier is present between two different regions of high-energy-density matter: a strongly heated tungsten wire and a surrounding plastic layer that stays relatively cold. We use diffraction-enhanced imaging to track the time evolution of density discontinuities and reconstruct the temperature evolution in the quasi-stationary stage. The clear signatures of a temperature jump demonstrate the importance of the ITR for strongly heated systems with far-reaching implications for interpreting experiments and applications like inertial confinement fusion.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Cation Crossover Limits Accessible Current Densities for Zero-Gap Alkaline CO2 Reduction to Ethylene

Traditional CO2 reduction systems often fail in an alkaline environment due to the interaction of CO2 with a high-pH electrolyte, where carbonate and bicarbonate ion formation results in potassium-containing salt precipitation. The presence of the salt crystals causes a reduction in the selectivity of the electrolyzer toward CO2 conversion. Here, the critical operational variables, which elicit the salting out process, are investigated (i.e., ion transport). When the electrolyzer exceeds a critical current density, H2 evolution dominates CO2 reduction due to salt formation, which is confirmed by postmortem cross-sectional SEM-EDS of the electrode. The critical current density decreases with an increasing membrane thickness or anolyte ionic strength. Cathode salt formation is mediated by the unmitigated crossover of cations from the anolyte to the cathode across an anion exchange membrane, through which cations are imperfectly excluded. It is likely that electric field-driven migration promotes an increase in concentration of potassium across the membrane, until, at the critical current density for that electrolyzer arrangement, the concentration of potassium and bicarbonate ions exceeds the solubility limit of KHCO3, leading to salt precipitation.

CO2 reduction↗
DOE OSTI2026

Post-puff SOL broadening on MAST-U under high-recycling conditions: evidence consistent with cross-field transport changes

Transient broadening of the scrape-off layer (SOL) density profile can modify main-chamber first-wall particle fluxes and divertor loading, yet its control parameters remain debated between divertor-regime transitions, neutral dynamics and changes in cross-field transport. We investigate fueling-driven SOL density-profile evolution and post-fueling relaxation on MAST-Upgrade (MAST-U) in ohmic L-mode, using two otherwise similar double-null Conventional Divertor discharges (I p = 450 kA, B T = 0.33 T) with identical 50 ms low-field side gas puffs; the only intentional difference is the puff start time. Upstream Thomson scattering shows that both discharges develop a transient far-SOL density profile modification, expressed as an increased SOL-width metric λ n e and a far-SOL enhancement consistent with a shoulder-like signature. In the earlier-puff case, the SOL-width metric remains elevated after puff termination and the post-puff decay time scales are systematically longer across the analysed radii. Outer-target Langmuir probes indicate high-recycling conditions during the analysed window (few-eV T e with radially peaked j sat,∥ and q ∥ profiles, without signatures of deep detachment). The outer-divertor collisionality proxy Λ div is elevated in both cases and does not discriminate between the different post-puff persistence. The target-integrated ion flux proxy ∫ J sat dA evolves nearly identically in both discharges when aligned to the puff start. Taken together, these observations suggest that the late/post-puff upstream SOL evolution is not set by parallel exhaust to the outer targetalone and is more consistent with upstream cross-field redistribution, with a possible role for plasma–neutral coupling and fueling geometry.

MAST-U↗
DOE OSTI2025

Revealing Coexisting Cu0–Cu+ Sites in Cu3N Nanoensembles for Selective CC Coupling of CO2 Under Low Overpotential

To address a long-existing debate on what copper species are responsible for efficient CC coupling, especially ethanol formation, in electrochemical CO2 reduction reaction, herein, a comprehensive study using Cu3N nanocubes with a uniform size and shape, alongside a single crystalline phase is reported. The Cu3N nanoensemble electrode has a remarkable Faradaic efficiency (FE) of 64% for ethanol production at a relatively low potential of -0.6 V versus reversible hydrogen electrode. Through in-operando X-ray absorption spectroscopy study, a dynamic phase evolution that directly correlates with changes in FE across varying applied potentials is observed. Notably, the nanoensemble with a composition of ≈71% Cu+ and 29% Cu0 is identified as being selective for ethanol formation at the low overpotential. Conversely, a predominantly metallic Cu phase formed at potentials more negative than -0.6 V favors the hydrogen evolution reaction. Density functional theory calculations at the Cu3N-Cu interface substantiate that the coexistence of Cu0-Cu+ not only energetically favors the ethanol reaction pathway but also destabilizes the intermediates for ethylene pathway.

*CH-CHO as an intermediate for ethanol production↗