Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “gradient methods”

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 1,513 records · Page 84
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

A compact x-ray spectrometer for measurements of electron temperature distributions in inertial confinement fusion implosions at OMEGA

The Wedge Range Filter (WRF), commonly used for proton spectroscopy at the OMEGA Laser Facility and National Ignition Facility, is adapted to measure the x-ray continuum spectrum through transmission measurement using a continuous-gradient filter. Continuum x rays emitted from the hotspot of an implosion contain information about the plasma composition and electron temperature. The WRF data are leveraged to probe this distribution, specifically the electron temperature distribution. In this work, the data recorded with the WRF are forward modeled using a temperature distribution model folded with the WRF response function. An uncertainty analysis is conducted through a Bayesian regression algorithm using a Hamiltonian Monte Carlo sampler. This analysis enables the uncertainties in the instrument response to be folded into the uncertainty estimation of the electron temperature and absolute x-ray emission. Data analysis for a series of OMEGA implosions is presented and compared with radiation hydrodynamic simulations.

Lasers↗
DOE OSTI2025

Examination of coal combustion management sites for microbiological and chemical signatures of groundwater impacts

Coal combustion accounts for 40% of the world’s electricity and generates more than a billion tons of coal combustion products (CCP) annually, half of which end up in landfills and impoundments. CCP contain mixtures of chemicals that can be mobile in the environment and impact the quality of surface water and potable groundwater. In this investigation, water samples from 14 coal combustion management sites across 4 physiographic regions in the United States, paired with background and down-gradient groundwater samples, were analyzed for water chemistry and microbiology. The objective was to determine if microbiology data alone, or supported by chemistry data, could reliably differentiate source waters and identify sites where CCP is known or expected to be influencing groundwater. Two percent of the total amplicons showed genus level conservation across CCP management sites, regions, and sample types; corresponding to ubiquitous, facultatively aerobic proteobacterial taxa that are generally recognized for the potential to respire using different terminal electron acceptors. Ordination plots did not reveal significant differences ( p > 0.05) in 16S rRNA gene amplicon diversity by CCP management site, water sample types, or physiographic regions. Contrastingly, chemistry distinguished sample types by standard water quality metrics (total dissolved solids, Ca:SO 4 ratio), alkali earth metals (K, Na, Li), selenium, boron, and fluoride. A focused evaluation of 16S rRNA gene amplicons for a subset of CCP management sites revealed microbiological features and chemical drivers (F, Ca, temperature) that positively identified the single CCP management site confirmed to have groundwater impacted by CCP leachate. At this site, 9 genera (>0.5% relative abundance) were exclusive to CCP porewater and downgradient groundwater. Inferred metabolisms for these taxa indicates potential for N and S biogeochemical transformations and 1-C metabolism that are consistent with a reducing environment, as evidenced by low ORP and depleted SO 4 2− . This research contributes to a growing understanding of conditions where these data types, analyses, and interpretation methods could be applied for distinguishing influence from CCP on the surrounding environment, as well as practical limitations.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2026

Soil nitrogen mineralization rates, nutrient stocks, stable isotopes, and water volumetric measurements across terrestrial-aquatic interfaces from three wetlands at the Tanglewood Biological Station, Alabama

This dataset supports a broader study investigating wetland hydrologic and biogeochemical responses to inundation events. Soil samples were collected across four sampling events along terrestrial-aquatic gradients at three wetland sites located within the Tanglewood Biological Station in Alabama from April 2024 to June 2025. The contents in this data package include soil in-situ nitrogen mineralization rates (measurements of net nitrification, net ammonification, and net mineralization), nutrient stocks (total carbon, total nitrogen, and organic matter), stable isotopes (carbon and nitrogen), and water volumetric measurements (water-filled pore space). Water level data related to each wetland location can be found at https://data.ess-dive.lbl.gov/view/doi:10.15485/2530253 (Kirker et al., 2024), related water geochemistry data can be found at https://data.ess-dive.lbl.gov/datasets/doi:10.15485/3001967 (Forbes et al., 2025), and related surface water sediment chemistry data can be found at https://data.ess-dive.lbl.gov/datasets/doi:10.15485/3377325 (Molina Serpas et al., 2026). In addition to this readme, this data package also includes a file-level metadata (FLMD) file that describes each file and a data dictionary (DD) that describes all column/row headers and variable definitions. This dataset is comprised of (1) file-level metadata; (2) data dictionary; (3) field metadata and international generic sample numbers (IGSNs); (4) readme; (5) the field protocol; and (6) a subfolder with sample data. The sample data subfolder contains (1) net nitrification rate, (2) net ammonification rate, (3) areal net mineralization rate, (4) percent organic matter, (5) water-filled pore space, (6) total carbon content, (7) total nitrogen content, (8) stable carbon isotope (delta carbon-13), and (9) stable nitrogen isotope (delta nitrogen-15), and (10) methods codes. All files are .csv or .pdf.

13-C↗
DOE OSTI2025

Baltimore Social-Environmental Collaborative (BSEC) Doppler Lidar & Derived Products

This repository contains all processed Doppler‐lidar outputs from the PSU lidar deployed for the Baltimore Social‐Environmental Collaborative (BSEC) project. Vertical Stare Scans (fixed‐beam, vertical profiling): 1 Hz backscatter intensity (m⁻¹ sr⁻¹), signal‐to‐noise ratio (unitless), and Doppler vertical‐velocity (m s⁻¹) on ~30 m range gates, stored as CF-compliant NetCDF. Wind Profiles (horizontal‐wind retrieval): daily NetCDF outputs of retrieved horizontal wind speed (m s⁻¹) and direction (degrees), computed from the angled‐scan returns. Profile Statistics (summary statistics on the vertical velocity): 15 min windows (default) of mean, variance, skewness, kurtosis, high-frequency variance, etc., as a function of height; saved as CF-compliant NetCDF files. Boundary Layer Height (BLH) (fuzzy-logic output): 15 min BLH estimates (m), with lower/upper fuzzy bounds (m) and a quality flag (0–4) indicating data status (e.g., no data, good, below range, ran out of signal, cloud-topped). Cloud Base Height (Haar-gradient detection): 15 min estimates of cloud-base height (m) with a cloud-detection quality flag (0–3: none, low, moderate, high). All five product streams are organized by year and date under their own top-level folders (01_Vertical_Stare_Scans/ through 05_Cloud_Height/). Each folder contains a data_ /YYYY/ subdirectory with daily CF-compliant NetCDF outputs (96 windows per day at 15 min intervals). Global attributes in each file include creation history, version (2.0.0), institution, and source. Instrument & MeasurementsThe PSU Doppler Lidar samples aerosol backscatter (m⁻¹ sr⁻¹), signal-to-noise ratio, and radial velocity at ~1 Hz. Vertical stare scans point the beam straight up; after collecting angled scans through multiple elevation angles, the "Wind Profiles" product contains the fully retrieved horizontal wind speed and direction. Data were collected continuously at ~30 m range resolution, with a typical height ceiling of ~12 km. How to Use Open any NetCDF with Python's xarray, MATLAB, or similar CF-compliant tools. Stare scans and angled-scan retrievals (Wind Profiles) are CF-compliant daily NetCDF files. Profile-Statistics, BLH, and Cloud Height files are daily 15 min summaries (96 time steps per file). Inspect the included variables (e.g., vertical_velocity_variance, wind_speed, BLH, cloud_base_height) for your analyses. Use the quality flags (BLH_flag, cloud_flag) to filter out poor-quality retrievals. For more information or questions about processing methods, please contact:Nicholas E. Prince ⟨nec5299@psu.edu⟩Penn State Department of Meteorology & Atmospheric Science

Air Quality↗
DOE OSTI2026

Baltimore Social-Environmental Collaborative (BSEC) Doppler Lidar & Derived Products

This repository contains all processed Doppler‐lidar outputs from the PSU lidar deployed for the Baltimore Social‐Environmental Collaborative (BSEC) project. Vertical Stare Scans (fixed‐beam, vertical profiling): 1 Hz backscatter intensity (m⁻¹ sr⁻¹), signal‐to‐noise ratio (unitless), and Doppler vertical‐velocity (m s⁻¹) on 30 m range gates (and 3 m range gates), stored as CF-compliant NetCDF. Wind Profiles (horizontal‐wind retrieval): daily NetCDF outputs of retrieved horizontal wind speed (m s⁻¹) and direction (degrees), computed from the angled‐scan returns. Profile Statistics (summary statistics on the vertical velocity): 15 min windows (default) of mean, variance, skewness, kurtosis, high-frequency variance, etc., as a function of height; saved as CF-compliant NetCDF files. Boundary Layer Height (BLH) (fuzzy-logic output): 15 min BLH estimates (m), with lower/upper fuzzy bounds (m) and a quality flag (0–4) indicating data status (e.g., no data, good, ran out of signal, below range, cloud-topped). Cloud Base Height (Haar-gradient detection): 10 min estimates of cloud-base height (m). All five product streams are organized by year and date under their own top-level folders (01_Vertical_Stare_Scans/ through 05_Cloud_Height/). Each folder contains a data_ /YYYY/ subdirectory with daily CF-compliant NetCDF outputs (96 windows per day at 15 min intervals). Global attributes in each file include creation history, version (3.0.0), institution, and source. Instrument & MeasurementsThe PSU Doppler Lidar samples aerosol backscatter (m⁻¹ sr⁻¹), signal-to-noise ratio, and radial velocity at ~1 Hz. Vertical stare scans point the beam straight up; after collecting angled scans through multiple elevation angles, the "Wind Profiles" product contains the fully retrieved horizontal wind speed and direction. Data were collected continuously at ~30 m range resolution (and 3 m for the year of 2025), with a typical height ceiling of ~12 km. How to Use Open any NetCDF with Python's xarray, MATLAB, or similar CF-compliant tools. Stare scans and angled-scan retrievals (Wind Profiles) are CF-compliant daily NetCDF files. Profile-Statistics, BLH, and Cloud Height files are daily 15 min (10 min for Cloud Heights) summaries (96 time steps per file). Inspect the included variables (e.g., vertical_velocity_variance, wind_speed, BLH, cloud_base_height) for your analyses. Use the quality flags (BLH_flag, cloud_flag) to filter out poor-quality retrievals. For more information or questions about processing methods, please contact:Nicholas E. Prince ⟨nec5299@psu.edu⟩Penn State Department of Meteorology & Atmospheric Science

Air Quality↗
DOE OSTI2026

Baltimore Social-Environmental Collaborative (BSEC) Doppler Lidar & Derived Products

This repository contains all processed Doppler‐lidar outputs from the PSU lidar deployed for the Baltimore Social‐Environmental Collaborative (BSEC) project. Vertical Stare Scans (fixed‐beam, vertical profiling): 1 Hz backscatter intensity (m⁻¹ sr⁻¹), signal‐to‐noise ratio (unitless), and Doppler vertical‐velocity (m s⁻¹) on 30 m range gates (and 3 m range gates), stored as CF-compliant NetCDF. Wind Profiles (horizontal‐wind retrieval): daily NetCDF outputs of retrieved horizontal wind speed (m s⁻¹) and direction (degrees), computed from the angled‐scan returns. Profile Statistics (summary statistics on the vertical velocity): 15 min windows (default) of mean, variance, skewness, kurtosis, high-frequency variance, etc., as a function of height; saved as CF-compliant NetCDF files. Boundary Layer Height (BLH) (fuzzy-logic output): 15 min BLH estimates (m), with lower/upper fuzzy bounds (m) and a quality flag (0–4) indicating data status (e.g., no data, good, ran out of signal, below range, cloud-topped). Cloud Base Height (Haar-gradient detection): 10 min estimates of cloud-base height (m). All five product streams are organized by year and date under their own top-level folders (01_Vertical_Stare_Scans/ through 05_Cloud_Height/). Each folder contains a data_ /YYYY/ subdirectory with daily CF-compliant NetCDF outputs (96 windows per day at 15 min intervals). Global attributes in each file include creation history, version (3.0.1), institution, and source. Instrument & MeasurementsThe PSU Doppler Lidar samples aerosol backscatter (m⁻¹ sr⁻¹), signal-to-noise ratio, and radial velocity at ~1 Hz. Vertical stare scans point the beam straight up; after collecting angled scans through multiple elevation angles, the "Wind Profiles" product contains the fully retrieved horizontal wind speed and direction. Data were collected continuously at ~30 m range resolution (and 3 m for the year of 2025), with a typical height ceiling of ~12 km. How to Use Open any NetCDF with Python's xarray, MATLAB, or similar CF-compliant tools. Stare scans and angled-scan retrievals (Wind Profiles) are CF-compliant daily NetCDF files. Profile-Statistics, BLH, and Cloud Height files are daily 15 min (10 min for Cloud Heights) summaries (96 time steps per file). Inspect the included variables (e.g., vertical_velocity_variance, wind_speed, BLH, cloud_base_height) for your analyses. Use the quality flags (BLH_flag, cloud_flag) to filter out poor-quality retrievals. For more information or questions about processing methods, please contact:Nicholas E. Prince ⟨nec5299@psu.edu⟩Penn State Department of Meteorology & Atmospheric Science

Air Quality↗
DOE OSTI2019

Comment on “Superconductivity at low density near a ferroelectric quantum critical point: Doped SrTiO 3 ”

Wölfle and Balatsky have proposed a microscopic pairing mechanism for doped SrTiO 3 (STO) based on the gradient coupling of electronic density to the soft transverse optical (TO) phonon mode. Since this coupling to TO phonons is usually weak, this conclusion is surprising, especially for a low-density superconductor such as STO, where the density of states is small. A crucial step in the argument made by Wölfle and Balatsky is that the displacement vector of the TO mode is not strictly perpendicular to the momentum vector, making a deformation coupling possible. We show that they have made a mistake in computing the eigenvector and have grossly overestimated this lack of orthogonality. When corrected, the coupling is negligible. We also use transport data to put upper bounds on the coupling constant which are much smaller than the estimate by Wölfle and Balatsky. Lastly, we also object to their use of the Eliashberg equation when the phonon frequency is larger than the Fermi energy.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2019

Variational control forces for enhanced sampling of nonequilibrium molecular dynamics simulations

We introduce a variational algorithm to estimate the likelihood of a rare event within a nonequilibrium molecular dynamics simulation through the evaluation of an optimal control force. Optimization of a control force within a chosen basis is made possible by explicit forms for the gradients of a cost function in terms of the susceptibility of driven trajectories to changes in variational parameters. We consider probabilities of time-integrated dynamical observables as characterized by their large deviation functions and find that in many cases, the variational estimate is quantitatively accurate. Additionally, we provide expressions to exactly correct the variational estimate that can be evaluated directly. We benchmark this algorithm against the numerically exact solution of a model of a driven particle in a periodic potential, where the control force can be represented with a complete basis. Additionally, we then demonstrate the utility of the algorithm in a model of repulsive particles on a line, which undergo a dynamical phase transition, resulting in singular changes to the form of the optimal control force. In both systems, we find fast convergence and are able to evaluate large deviation functions with significant increases in statistical efficiency over alternative Monte Carlo approaches.

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

Vegetation classification map and covariates associated with NEON AOP survey, East River, CO 2018

This package includes geospatial data layers developed to investigate how environmental gradients—specifically topography and near-surface soil properties—drive the spatial arrangement of dominant plant communities in mountainous watersheds. The geospatial products, which support the analysis of these ecological relationships, are derived from airborne hyperspectral and LiDAR datasets acquired by the National Ecological Observatory Network (NEON) Airborne Observation Platform (AOP), in conjunction with an extensive ground field campaign conducted in summer 2018. This work is part of the DOE Watershed Function Science Focus Area (SFA) and features geospatial datasets developed based on observations and ground data collected at East River, Colorado, in collaboration with the National Ecological Observatory Network (NEON) Airborne Observation Platform (AOP) survey in June 2018. Classification Map: - Classification Map (PNG, GeoTIFF): Derived from hyperspectral and LiDAR airborne data using a machine learning approach. - Class Code Mapper (CSV): Associates pixel values with corresponding vegetation/non-vegetation classes. - Classification Reference Data (CSV): Reference data used in the machine learning procedure. LiDAR-Derived Products: - Topographical Metrics (GeoTIFFs): Elevation, slope, curvature, TWI, TPI, solar insolation, and canopy height model (CHM), smoothed with a 5x5 pixel window. Vegetation Indices: - GeoTIFFs of NDVI, NDNI, NDWI: Vegetation indices derived from hyperspectral data. Urban Masks: - Urban Mask (GeoTIFF): Applied to the mapping to convert bare soil classes to urban classes. Software Compatibility: GeoTIFFs: Can be visualized with GIS software or libraries that support GeoTIFF images. CSV Files: Can be opened with any software that handles comma-separated values. The FLMD file provides details and links to the source datasets used to derive the products. The manuscript (in the Method session) provides details on how each product was derived. This work was supported by the Watershed Function Science Focus Area at Lawrence Berkeley National Laboratory funded by the US Department of Energy, Office of Science, Biological and Environmental Research under Contract No. DE-AC02-05CH11231. Update on 2026-03-25: Since the original dataset publication date of 02/28/2020, this package has a new classification map derived by an improved methodology. This update also includes additional ground data that improved the representation of some of the communities. See the methods for further details on what has changed between versions.

2018 NEON and 2025 CHESS Campaigns↗
DOE OSTI2024

The rise and fall of stretched bond errors: Extending the analysis of Perdew–Zunger self-interaction corrections of reaction barrier heights beyond the LSDA

Incorporating self-interaction corrections (SIC) significantly improves chemical reaction barrier height predictions made using density functional theory methods. We present a detailed orbital-by-orbital analysis of these corrections for three semi-local density functional approximations (DFAs) situated on the three lowest rungs of Jacob’s ladder of approximations. The analysis is based on Fermi–Löwdin Orbital Self-Interaction Correction (FLOSIC) calculations performed at several steps along the reaction pathway from the reactants (R) to the transition state (TS) to the products (P) for four representative reactions selected from the BH76 benchmark set. For all three functionals, the major contribution to self-interaction corrections of the barrier heights can be traced to stretched bond orbitals that develop near the TS configuration. The magnitude of the ratio of the self-exchange–correlation energy to the self-Hartree energy (XC/H) for a given orbital is introduced as an indicator of one-electron self-interaction error. XC/H = 1.0 implies that an orbital’s self-exchange–correlation energy exactly cancels its self-Hartree energy and that the orbital, therefore, makes no contribution to the SIC in the FLOSIC scheme. For the practical DFAs studied here, XC/H spans a range of values. Here, the largest values are obtained for stretched or strongly lobed orbitals. We show that significant differences in XC/H for corresponding orbitals in the R, TS, and P configurations can be used to identify the major contributors to the SIC of barrier heights and reaction energies. Based on such comparisons, we suggest that barrier height predictions made using the strongly constrained and appropriately normed meta-generalized gradient approximation may have attained the best accuracy possible for a semi-local functional using the Perdew–Zunger SIC approach.

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

Quantifying electron temperature distributions from time-integrated x-ray emission spectra

K-shell x-ray emission spectroscopy is a standard tool used to diagnose the plasma conditions created in high-energy-density physics experiments. In the simplest approach, the emissivity-weighted average temperature of the plasma can be extracted by fitting an emission spectrum to a single temperature condition. It is known, however, that a range of plasma conditions can contribute to the measured spectra due to a combination of the evolution of the sample and spatial gradients. In this work, we define a parameterized model of the temperature distribution and use Markov Chain Monte Carlo sampling of the input parameters, yielding uncertainties in the fit parameters to assess the uniqueness of the inferred temperature distribution. Here we present the analysis of time-integrated S and Fe x-ray spectroscopic data from the Orion laser facility and demonstrate that while fitting each spectral region to a single temperature yields two different temperatures, both spectra can be fit simultaneously with a single temperature distribution. We find that fitting both spectral regions together requires a maximum temperature of $1310^{+90}_{-70}$ eV with significant contributions from temperatures down to 200 eV.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Plasma gradient effect on direct laser acceleration

The transfer of a high-intensity laser pulse energy to a high-energy electron beam via the direct laser acceleration mechanism is shown to be significantly enhanced through control of the plasma density gradient. Experiments performed using the OMEGA EP facility's high-intensity beams altered the plasma density and gradients by changing the Mach number and the angle of the gas-jet nozzle to the laser axis. When a long density gradient at the rear of the target is used, the total high-energy electron number measured was enhanced by 4.5 times compared to a shorter rear gradient. Complementary two-dimensional simulations, which follow the laser field evolution and the corresponding electron dynamics, strongly support the key trends observed in the experiment. The effect is twofold, the long density gradient provides the longest acceleration distance while it minimizes the formation of the sheath field as the electron beam exits into the vacuum. This study shows the importance of tailoring the plasma density.

Laser plasma interactions↗
DOE OSTI2023

Joint Inversion of Surface Electrical Resistivity Tomography and Seismic Refraction Data between the 200 Areas

Geologic stratigraphy on the Hanford Site influences groundwater and contaminant migration through the aquifer system and the vadose zone. The current geologic framework model (GFM) relies heavily on a sparse distribution of borehole data in some locations to map geologic contacts and hydrologic properties in the subsurface. Non-invasive geophysical methods such as electrical resistivity tomography (ERT), transient electromagnetic surveying, and seismic imaging are being used at Hanford to map subsurface structure in areas with limited well observations. This is to develop and mature the capability of geophysical methods to aid in GFM refinement, to identify regions of subsurface complexity, and for optimal well siting. A joint inversion of co-located seismic refraction and ERT data was carried out for data collected on a ~2.3-km profile between the 200 Areas on the Hanford Site. While ERT and seismic refraction images have sensitivity to overlapping physical properties (porosity, moisture content, lithology), the resolution and physics used to acquire each of these datasets are different and therefore information can be different or mutually complementary. Performing a joint inversion provides a reasonable option for a coherent, coupled interpretation for mutually complementary datasets. Between the 200 Areas, there are few boreholes to interpret the geologic framework model, and these data sets were obtained to provide a first line of evidence toward identifying stratigraphic structure. The seismic refraction and ERT data were independently inverted during fiscal year 2022 and broadly showed a two-layer structure with a trough-like feature that is ~1 km wide and upwards of 150 m deep. The depth of the trough feature was greater in the ERT image compared to the seismic image, which indicated a maximum depth of approximately 110 m. The objective of the joint inversion described in this report was to invert the seismic refraction and ERT data together while constraining the ERT image to be structurally similar to the seismic refraction image. The approach was applied using the geophysical inverse modeling program E4D, which has the capability to invert first-arrival times from seismic refraction data and ERT resistances using a “cross-gradient” constraint. The application of cross-gradient constraints with different weights produces ERT models that show a high degree of similarity within the upper 100 m (above ~120 m elevation). None of the ERT models show an improved structural similarity to the seismic result; therefore, it is recommended that further attempts to jointly interpret these models focus on petrophysics and image resolution. Petrophysical measurements of core samples would improve knowledge of what drives the ERT response in this region and, along with downhole geophysical measurements, could be used to “ground truth” the surface-based geophysical results. Image resolution studies would provide insight into which regions of the inverted images are reliable and which regions are poorly constrained.

58 GEOSCIENCES↗
DOE OSTI2021

Reduced order models for Lagrangian hydrodynamics

It is reported, as a mathematical model of high-speed flow and shock wave propagation in a complex multimaterial setting, Lagrangian hydrodynamics is characterized by moving meshes, advection-dominated solutions, and moving shock fronts with sharp gradients. These challenges hinder the existing projection-based model reduction schemes from being practical. We develop several variations of projection-based reduced order model techniques for Lagrangian hydrodynamics by introducing three different reduced bases for position, velocity, and energy fields. A time-windowing approach is also developed to address the challenge imposed by the advection-dominated solutions. Lagrangian hydrodynamics is formulated as a nonlinear problem, which requires a proper hyper-reduction technique. Therefore, we apply the over-sampling DEIM and SNS approaches to reduce the complexity due to the nonlinear terms. Finally, we also present both a posteriori and a priori error bounds associated with our reduced order model. We compare the performance of the spatial and time-windowing reduced order modeling approaches in terms of accuracy and speed-up with respect to the corresponding full order model for several numerical examples, namely Sedov blast, Gresho vortices, Taylor-Green vortices, and triple-point problems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Intrinsic Anomalous Hall Conductivity in a Nonuniform Electric Field

We study how the intrinsic anomalous Hall conductivity is modified in two-dimensional crystals with broken time-reversal symmetry due to weak inhomogeneity of the applied electric field. Focusing on a clean noninteracting two-band system without band crossings, we derive the general expression for the Hall conductivity at small finite wave vector q to order q 2 , which governs the Hall response to the second gradient of the electric field. Using the Kubo formula, we show that the answer can be expressed through the Berry curvature, Fubini-Study quantum metric, and the rank-3 symmetric tensor which is related to the quantum geometric connection and physically corresponds to the gauge-invariant part of the third cumulant of the position operator. We further compare our results with the predictions made within the semiclassical approach. By deriving the semiclassical equations of motion, we reproduce the result obtained from the Kubo formula in some limits. We also find, however, that the conventional semiclassical description in terms of the definite position and momentum of the electron is not fully consistent because of singular terms originating from the Heisenberg uncertainty principle. We thus present a clear example of a case when the semiclassical approach inherently suffers from the uncertainty principle, implying that it should be applied to systems in nonuniform fields with extra care.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Many-Body Effects in FeN4 Center Embedded in Graphene

We introduce a computational approach to study porphyrin-like transition metal complexes, bridging density functional theory and exact many-body techniques, such as the density matrix renormalization group (DMRG). We first derive a multi-orbital Anderson impurity Hamiltonian starting from first principles considerations that qualitatively reproduce generalized gradient approximation (GGA)+U results when ignoring inter-orbital Coulomb repulsion U ′ and Hund exchange J. An exact canonical transformation is used to reduce the dimensionality of the problem and make it amenable to DMRG calculations, including all many-body terms (both intra- and inter-orbital), which are treated in a numerically exact way. We apply this technique to FeN 4 centers in graphene and show that the inclusion of these terms has dramatic effects: as the iron orbitals become single occupied due to the Coulomb repulsion, the inter-orbital interaction further reduces the occupation, yielding a non-monotonic behavior of the magnetic moment as a function of the interactions, with maximum polarization only in a small window at intermediate values of the parameters. Furthermore, U ′ changes the relative position of the peaks in the density of states, particularly on the iron d z 2 orbital, which is expected to affect the binding of ligands greatly.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Improving deep learning model performance under parametric constraints for materials informatics applications

Abstract Modern machine learning (ML) and deep learning (DL) techniques using high-dimensional data representations have helped accelerate the materials discovery process by efficiently detecting hidden patterns in existing datasets and linking input representations to output properties for a better understanding of the scientific phenomenon. While a deep neural network comprised of fully connected layers has been widely used for materials property prediction, simply creating a deeper model with a large number of layers often faces with vanishing gradient problem, causing a degradation in the performance, thereby limiting usage. In this paper, we study and propose architectural principles to address the question of improving the performance of model training and inference under fixed parametric constraints. Here, we present a general deep-learning framework based on branched residual learning (BRNet) with fully connected layers that can work with any numerical vector-based representation as input to build accurate models to predict materials properties. We perform model training for materials properties using numerical vectors representing different composition-based attributes of the respective materials and compare the performance of the proposed models against traditional ML and existing DL architectures. We find that the proposed models are significantly more accurate than the ML/DL models for all data sizes by using different composition-based attributes as input. Further, branched learning requires fewer parameters and results in faster model training due to better convergence during the training phase than existing neural networks, thereby efficiently building accurate models for predicting materials properties.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Evaluation of the Effect of Burnup Acceleration on UO 2 Microstructure Evolution

Accelerated fuel qualification (AFQ) is a methodology by which new nuclear fuels are developed in an accelerated time frame compared with historical fuel qualification approaches. AFQ generally relies on high-fidelity physics-based modeling and simulation tools to adequately describe fuel performance as well as on revolutionary methods to accelerate burnup accumulation and collect relevant data more quickly. This report summarizes the use of advanced fuel modeling and simulation tools to evaluate microstructures from commercially irradiated fuel and microstructures from proposed MiniFuel irradiations, in which burnup accumulation is accelerated while prototypic temperature conditions are maintained. In this milestone, we used the mesoscale fuel performance code MARMOT to model the evolution of irradiated UO 2 microstructures and their potential restructuring at high burnup. The simulation conditions were informed by BISON models of both commercially irradiated fuel and MiniFuel. A first set of simulations investigated the recrystallization behavior of fully dense microstructures and showed full recrystallization at burnups as low as 52 MWd/kgU at 950°C. However, these simulations did not account for the presence of fission gas bubbles (predicted by BISON). Therefore, a second set of simulations including fission gas bubbles was performed and indicated that at the lowest temperature considered (650°C), the porous UO 2 microstructures have the highest total Gibbs free energies and are likely to recrystallize earlier than higher temperature cases (800 and 950°C), which agrees with high-burnup fuel characterization data. The results also showed that at lower temperature (650°C), the total free energies of the PWR fuel and MiniFuel microstructures are not significantly different. However, at the highest temperature (950°C), MiniFuel microstructures have a lower free energy than that of the PWR fuel microstructure. The competing effects between the temperature-dependent grain nucleation rate and the reduction of the free energy of the microstructure at higher temperature as a result of diffusion indicated that restructuring may occur at even higher temperatures than those considered in this study. In addition to the microstructure evolution modeling efforts, the burnup gradient across a single fuel specimen was also considered. This evaluation was for the VXF-15 position of the High Flux Isotope Reactor (HFIR) using the code suite HFIRCON, which was developed to automate the workflow for evaluating targets and fuel as they are irradiated in HFIR. The burnup gradient evaluation showed a dependence on both the axial and radial locations within the specimen, with a maximum difference of 1.7 between the inner and outermost radial layers. This relationship was further supported by considering the fission product speciation with respect to location within the specimen, which showed a higher concentration of 239 Pu, 240 Pu, and 241 Pu on the outer radial locations of the specimen than the center. The findings of the burnup and speciation evaluation show that some amount of self-shielding is occurring in the specimen when irradiated in the high-flux environment of HFIR; however, this impact is more pronounced for natural uranium when compared to 6% enrichment due to the higher ratio of 238 U in the specimen. Further analyses are required to understand the sensitivity of this gradient to spatial mesh and enrichment of the specimen.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗