Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “gradient sampling”

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

Coarse-Grained Simulations of Polyrotaxane Hydrogels under Quiescent and Shear Conditions

Cyclodextrin-based polyrotaxanes (PR) form hydrogels in water through cyclodextrin (CD) aggregation and crystallization. These networks can break under shear flows, making them versatile platforms for extrusion-based 3D printing. To optimize the material properties of 3D-printed PR gels, a microscopic understanding of the structural evolution during 3D printing is necessary. Here, we employ coarse-grained (CG) simulations to reveal the PR assembly process at the molecular level. Our simulations reproduce the experimental crystal morphologies of PR at varying concentrations and chain lengths under quiescent conditions. Using nonequilibrium simulations, we show shear flow ruptures crystalline domain connectivity in PR gels and stacks the lamellae in the gradient direction, allowing the materials to flow during 3D printing. After the cessation of flow, the anisotropic crystal alignment and absence of available dangling PR diminish the intercrystal connectivity. The printed materials are therefore mechanically weaker than the pristine hydrogels, in agreement with experimental results. Nonetheless, by relating the microscopic structural evolution with viscoelastic properties of PR gels and solutions, we elucidate how flow conditions and sample composition affect the 3D printing performance of PR hydrogels.

Smith, Cameron D. [Dartmouth College, Hanover, NH ↗
DOE OSTI2020

Coupling High-Throughput Experiments and Regression Algorithms to Optimize PGM-Free ORR Electrocatalyst Synthesis

Over the past decades, significant improvement has been achieved in the performance of platinum group metal-free (PGM-free) materials as an alternative to Pt-based electrocatalysts for oxygen reduction reaction (ORR). However, further progress in ORR activity requires evaluation of precursors and synthesis approaches. In response to this challenge, we generated a first of its kind experimental data set of 36 samples using high-throughput synthesis and activity measurements. Several control parameters (e.g., Fe precursor identity, the precursor content, and pyrolysis temperature) were varied. We then developed several state-of-the-art machine learning (ML) based regression models to predict ORR activity, dependent on selected synthesis variables. Through an iterative algorithm, higher prediction accuracy (smaller root-mean-square error) was achieved. We identified that gradient boosting regression (GBR) and support vector regression (SVR), among several methods, work best for this data set. Aided by our ML-based surrogate models, we decided to alter catalyst synthesis conditions, which resulted in a 36% increase in measured ORR activity in comparison to the maximum ORR mass activity value of 21.9 A/g catalyst in the original data set. Overall, this combined experiment and machine learning approach represents a promising path forward toward developing highly efficient next-generation ORR electrocatalysts and, more generally, functional materials.

25 ENERGY STORAGE↗
DOE OSTI2022

Process-based climate change assessment for European winds using EURO-CORDEX and global models

Wind energy is an important pillar of decarbonization strategies and potentially vulnerable to climate change. Existing wind climate change assessments rely on climate models but a systematic investigation of the global-to-regional climate modeling chain is missing. In this study, I highlight key limitations, namely (a) the differing representation of land use change in global and regional climate models which compromises comparability, and (b) the consistency of large-scale features along the global-to-regional climate modeling chain. To this end, I analyze the large European Coordinated Downscaling Experiment (EURO-CORDEX) ensemble (rcp85: N = 49; rcp45: N= 18; rcp26: N = 22) along with the driving global models (rcp85: N = 7; rcp45: N = 5; rcp26: N = 7), finding evidence that climate change reduces mean wind speeds by up to -0.8 ms -1 (offshore) and -0.3 ms -1 (onshore). I provide physical explanations for these changes by identifying two key drivers. First, onshore wind speeds drop in the driving global models in regions and scenarios with strong land use change but show no drop in EURO-CORDEX where land use is held constant. Second, offshore wind reductions follow decreases in the equator-to-pole temperature gradient remarkably well with correlations reaching around 0.9 in resource-rich European countries like Ireland, the United Kingdom and Norway, implying that arctic amplification is a severe risk for European offshore wind energy. My results suggest that earlier conclusions of negligible climate change impacts on wind energy might be premature if either land use changes strongly or polar amplification is at or above the range sampled in global climate models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

HfTa_BCC_SolidSolution_128atoms_VASP6

We performed density functional theory (DFT) calculations for body-centered-cubic (BCC) structures with 128 lattices sites of solid solution binary alloys hafnium-tantalum (Hf-Ta). The electronic structures of alloys have been calculated using Vienna Ab initio Simulation Package (VASP). Within this package the DFT approach is used to reduce many-body Schrodinger equation to set of single particle Kohn-Sham (KS) equations. The generalized electronic exchange-correlation functional is described by generalized gradient approximation with the Perdew-Burke-Ernzerhof parametrization. The electron-ion interactions is described by pseudopotentials developed within the plane-wave basis projector augmented-wave (PAW) approach. These pseudopotentials are available at the VASP portal (http://cms.mpi.univie.ac.at/vasp/). Our calculations have been run with the pseudopotentials treating s and p semi-core states as valence in case for the elements Hf and Ta. The electronic densities and potentials are expanded over plane-waves with energy cutoff of 350 eV. 2x2x2 k-mesh and normal precision were used. The alloys were modeled by supercell containing 128 randomly distributed atoms. At initial step the atoms occupy perfect BCC lattice sites. This initial structure was optimized until energy changes less than 1e-6 eV, while forces acting on atoms don't exceed 1e-2 eV/angstrom. The electron-ion interaction is described by PAW pseudopotentials. The calculations have been collected by sampling chemical compositions across the entire compositional range. The chemical compositions have been sampled by progressively changing the number of atoms per constituent by 4. For each chemical composition of binaries and ternaries, the first-principle calculations have been run for 100 randomized arrangements of the constituents on the BCC lattice sites. We collected data for a total of 3,040 randomized atomic structures over 31 chemical compositions. Further methodological and structural information is contained in the dataset README.txt file.

36 MATERIALS SCIENCE↗
DOE OSTI2022

REGAL International Program: Analysis of experimental data for depletion code validation

The Rod-Extremity and Gadolinia AnaLysis (REGAL) Program is a joint international effort to expand the nuclide inventory experimental data for irradiated nuclear fuel, with a specific focus on addressing two challenging needs associated with the characterization of modern, high duty, nuclear fuel. The first challenge is filling the gaps in experimental nuclide inventory data for gadolinia (UO 2 –Gd 2 O 3 ) fuel rods. The huge absorption cross sections of Gd-155 and Gd-157 in the Gd dopant in these rods lead to atypical spatial self-shielding patterns and have an impact on the neutronic environment within the fuel assembly compared to regular UO 2 fuel rods. Additionally, the second challenge is investigating the impact of burnup gradients at rod extremities on fuel composition and neutron leakage, to provide relevant experimental data for assessing computational capabilities to model such impact. A benchmark has been defined as a first step in the development of best-estimate models in the preliminary phase of the experimental data evaluation. Comparison of experimental results obtained in Phase I of the program for two measured pressurized water reactor (PWR) samples, one UO 2 and one UO 2 –Gd 2 O 3 sample, with calculated results obtained with different computational tools based on the defined benchmark are presented and discussed.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

Macroscopic trends of linear tearing stability in cylindrical current profiles

Abstract The likelihood of realising tokamak power-plants will be greatly improved by the discovery of high-gain equilibria that resist the formation of small islands and hence avoid the disruptive neoclassical tearing mode. We propose a series of studies to understand how simple tokamak design can leverage aspects of tearing onset physics to maximise passive resistance to island formation. Here we investigate the variation that current profiles can bring about in preventing tearing onset through the cylindrical linear tearing stability parameter Δ ′ . A database of 159148 realistic pilot-plant current profiles was generated with Monte Carlo sampling, and the distribution of Δ ′ values was linked with interpretable profile characteristics. In agreement with prior theoretical and experimental studies, Δ ′ was found to be strongly correlated with the existence and steepness of a local toroidal current well or hill, with the former destabilising and the latter stabilising. In the absence of these two cases, the remaining Δ ′ values were linearly bounded by the toroidal current gradient at the rational surface.

Physics↗
DOE OSTI2025

CO2 and CH4 leaf-level fluxes and soil porewater concentrations from common vegetation patches in Louisiana’s coastal wetlands

This dataset contains leaf-level flux and soil porewater concentration measurements of carbon dioxide (CO2) and methane (CH4 ) in plots in the footprint of Ameriflux sites US-LA2 and US-LA3. Leaf fluxes in US-LA2 were measured on patches dominated by Sagittaria lancifolia and co-dominated by Sagittaria lancifolia and Typha latifolia. In US-LA3, fluxes were measured from distinct Juncus roemerianus and Spartina alterniflora patches. The porewater concentrations were collected across a vertical profile (~50 cm depth) at centric locations within 25 m2 plots where we measured the leaf fluxes. US-LA3 included an additional set of measurements in open water spots. We aimed to evaluate differences in leaf fluxes and porewater pools of CO2 and CH4 of representative ecohydrological patches across a salinity gradient. We also used this dataset to help develop ELM-Wet, a more realistic representation of wetland carbon biogeochemical processes within the U.S. Department of Energy’s Energy Exascale Earth System Model (E3SM) Land Model version 1 (ELM v.1). The files can be opened with regular text editors or spreadsheet programs. Version 2.0 (8/26/2025): This is the latest version of this dataset. The update includes additional samples of soil porewater CH4/CO2 concentrations from June-2021 to November-2022, as well as minor adjustments made to V1 samples via changing Henry's solubility to account for porewater salinity. Additionally leaf-level measurments of spectral indices, PSRI, NDVI, and PRI have been added to complement Leaf-level flux measurements. All V1 data sets have been integrated into V2 sheets, ensuring data from the previous version is contained with the additional samples and consistent with V2 metadata.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Physics-Driven Convolutional Autoencoder Approach for CFD Data Compressions: Preprint

With the growing size and complexity of turbulent flow models, data compression approaches are of the utmost importance to analyze, visualize, or restart the simulations. Recently, in-situ autoencoder-based compression approaches have been proposed and shown to be effective at producing reduced representations of turbulent flow data. However, these approaches focus solely on training the model using point-wise sample reconstruction losses that do not take advantage of the physical properties of turbulent flows. In this paper, we show that training autoencoders with additional physics-informed regularizations, e.g., enforcing incompressibility and preserving enstrophy, improves the compression model in three ways: (i) the compressed data better conform to known physics for homogeneous isotropic turbulence without negatively impacting point-wise reconstruction quality, (ii) inspection of the gradients of the trained model uncovers changes to the learned compression mapping that can facilitate the use of explainability techniques, and (iii) as a performance byproduct, training losses are shown to converge up to 12x faster than the baseline model.

auto-encoders↗
DOE OSTI2020

Hard, transparent, sp3-containing 2D phase formed from few-layer graphene under compression

Despite several theoretically proposed two-dimensional (2D) diamond structures, experimental efforts to obtain such structures are in initial stage. Recent high-pressure experiments provided significant advancements in the field, however, expected properties of a 2D-like diamond such as sp 3 content, transparency and hardness, have not been observed together in a compressed graphene system. Here, we compress few-layer graphene samples on SiO 2 /Si substrate in water and provide experimental evidence for the formation of a quenchable hard, transparent, sp 3 -containing 2D phase. Our Raman spectroscopy data indicates phase transition and a surprisingly similar critical pressure for two-, five-layer graphene and graphite in the 4e6 GPa range, as evidenced by changes in several Raman features, combined with a lack of evidence of significant pressure gradients or local non-hydrostatic stress components of the pressure medium up to ≈ 8 GPa. The new phase is transparent and hard, as evidenced from indentation marks on the SiO 2 substrate, a material considerably harder than graphene systems. Furthermore, we report the lowest critical pressure (≈ 4 GPa) in graphite, which we attribute to the role of water in facilitating the phase transition. Overall, theoretical calculations and experimental data indicate a novel, surface to- bulk phase transition mechanism that gives hint of diamondene formation.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Active matter shepherding and clustering in inhomogeneous environments

Here, we consider a mixture of active and passive run-and-tumble disks in an inhomogeneous environment where only half of the sample contains quenched disorder or pinning. The disks are initialized in a fully mixed state of uniform density. We identify several distinct dynamical phases as a function of motor force and pinning density. At high pinning densities and high motor forces, there is a two-step process initiated by a rapid accumulation of both active and passive disks in the pinned region, which produces a large density gradient in the system. This is followed by a slower species phase separation process where the inactive disks are shepherded by the active disks into the pin-free region, forming a nonclustered fluid and producing a more uniform density with species phase separation. For higher pinning densities and low motor forces, the dynamics becomes very slow and the system maintains a strong density gradient. For weaker pinning and large motor forces, a floating clustered state appears, and the time-averaged density of the system is uniform. We illustrate the appearance of these phases in a dynamic phase diagram.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

PyTREES

PyTREES (Python tool for Training/Testing Robust Explainable Ensembles on Spectra) is software that implements a data-driven approach to predicting the amount of specific oxides present in materials samples of laser-induced breakdown spectroscopy (LIBS); such as from the ChemCam instrument suite onboard the NASA Curiosity rover. PyTREES is designed to input LIBS data in the format provided by the ChemCam team [1]. PyTREES then applies appropriate pre-processing to this data [2], and implements several regression methods for predicting oxides from spectra. The regression methods include: ensemble methods (random forest, extra trees, and gradient boosting regression) and blended submodels using the “double blending” technique. PyTREES additionally implements methods for quantifying the importance of features in regression model: (1) mean decrease in impurity (MDI) and (2) permutation importance to investigate the wavelengths used by the regression methods. [1] Gasda et al. (2021). Spectrochim Acta B, 181, 106223. [2] Clegg et al. (2017). Spectrochim Acta B , 129, 64–85.

Oyen, Diane↗
DOE OSTI2021

Molecular-level insights into structure and dynamics in ionic liquids and polymer gel electrolytes

We report designing new electrolytes requires a better understanding of the correlation between their transport properties and their molecular structure. In this work, we present a detailed study of ionic liquids and polymer gel electrolytes probing their structure and dynamics by nuclear magnetic resonance (NMR) spectroscopy. In particular, ammonium- and phosphonium-based ionic liquids (ILs) are combined with different LiTFSI concentrations, and then with different ratios of poly(methylmethacrylate) polymer. The temperature dependence of self-diffusion coefficients of mobile species, D Li+ , D TFSI- and D P4441+ or D N4441+ , measured by pulsed field gradient (PFG) NMR spectroscopy obeys the Arrhenius equation. Diffusivity of [P4441][TFSI] ILs is found to be greater than those of the [N4441][TFSI] ILs in both liquid and gel electrolytes. Solid state NMR experiments including 13 C and 19 F MAS, 13 C{ 19 F}, 13 C{ 1 H} CPMAS probe the local structure and molecular-level interactions between ions and polymer in gel electrolytes, particularly for samples with higher PMMA content (≥25 wt%). Finally, the fast field cycling relaxometry has been used to unveil the rotational and translational dynamics of P4441 + or N4441 + and TFSI - by measuring 1 H and 19 F R 1 relaxation rate profiles at different temperatures. A comprehensive NMR analysis including relaxation studies at low magnetic field provides decisive new insights regarding the formation of ionic clusters and the interaction of ions with the polymer chain in the case of the gel electrolytes.

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

Structure optimization with stochastic density functional theory

Linear-scaling techniques for Kohn–Sham density functional theory are essential to describe the ground state properties of extended systems. Still, these techniques often rely on the localization of the density matrix or accurate embedding approaches, limiting their applicability. In contrast, stochastic density functional theory (sDFT) achieves linear- and sub-linear scaling by statistically sampling the ground state density without relying on embedding or imposing localization. In return, ground state observables, such as the forces on the nuclei, fluctuate in sDFT, making optimizing the nuclear structure a highly non-trivial problem. In this work, we combine the most recent noise-reduction schemes for sDFT with stochastic optimization algorithms to perform structure optimization within sDFT. We compare the performance of the stochastic gradient descent approach and its variations (stochastic gradient descent with momentum) with stochastic optimization techniques that rely on the Hessian, such as the stochastic Broyden–Fletcher–Goldfarb–Shanno algorithm. In conclusion, we further provide a detailed assessment of the computational efficiency and its dependence on the optimization parameters of each method for determining the ground state structure of bulk silicon with varying supercell dimensions.

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

Search for Time-Dependent C P Violation in D 0 → π + π − π 0 Decays

A measurement of time-dependent C P violation in D 0 → π + π − π 0 decays using a p p collision data sample collected by the LHCb experiment in 2012 and from 2015 to 2018, corresponding to an integrated luminosity of 7.7 fb − 1 , is presented. The initial flavor of each D 0 candidate is determined from the charge of the pion produced in the D * ( 2010 ) + → D 0 π + decay. The decay D 0 → K − π + π 0 is used as a control channel to validate the measurement procedure. The gradient of the time-dependent C P asymmetry Δ Y in D 0 → π + π − π 0 decays is measured to be Δ Y = ( − 1.3 ± 6.3 ± 2.4 ) × 10 − 4 , where the first uncertainty is statistical and the second is systematic, which is compatible with C P conservation. © 2024 CERN, for the LHCb Collaboration 2024 CERN

Aaij, R. (ORCID:0000000305331952)↗
DOE OSTI2021

Trends in Spitzer Secondary Eclipses

It is well established that the magnitude of the incident stellar flux is the single most important factor in determining the day–night temperature gradients and atmospheric chemistries of short-period gas giant planets. However, it is likely that other factors, such as planet-to-planet variations in atmospheric metallicity, C/O ratio, and cloud properties, also contribute to the observed diversity of infrared spectra for this population of planets. In this study, we present new 3.6 and 4.5 μm secondary eclipse measurements for five transiting gas giant planets: HAT-P-5b, HAT-P-38b, WASP-7b, WASP-72b, and WASP-127b. We detect eclipses in at least one bandpass for all five planets and confirm circular orbits for all planets except WASP-7b, which shows evidence for a nonzero eccentricity. Building on the work of Garhart et al., we place these new planets into a broader context by comparing them with the sample of all planets with measured Spitzer secondary eclipses. We find that incident flux is the single most important factor for determining the atmospheric chemistry and circulation patterns of short-period gas giant planets. Although we might also expect surface gravity and host star metallicity to play a secondary role, we find no evidence for correlations with either of these two variables.

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

HfTi_BCC_SolidSolution_128atoms_VASP6

We performed density functional theory (DFT) calculations for body-centered-cubic (BCC) structures with 128 lattices sites of solid solution binary alloys hafnium-titanium (Hf-Ti). The electronic structures of alloys have been calculated using Vienna Ab initio Simulation Package (VASP). Within this package the DFT approach is used to reduce many-body Schrodinger equation to set of single particle Kohn-Sham (KS) equations. The generalized electronic exchange-correlation functional is described by generalized gradient approximation with the Perdew-Burke-Ernzerhof parametrization. The electron-ion interactions is described by pseudopotentials developed within the plane-wave basis projector augmented-wave (PAW) approach. These pseudopotentials are available at the VASP portal (http://cms.mpi.univie.ac.at/vasp/). Our calculations have been run with the pseudopotentials treating s and p semi-core states as valence in case for the elements Hf and Ti. The electronic densities and potentials are expanded over plane-waves with energy cutoff of 350 eV. 2x2x2 k-mesh and normal precision were used. The alloys were modeled by supercell containing 128 randomly distributed atoms. At initial step the atoms occupy perfect bcc lattice cites. This initial structure was optimized until energy changes less than 1e-6 eV, while forces acting on atoms don't exceed 1e-2 eV/angstrom. The electron-ion interaction is described by PAW pseudopotentials. The calculations have been collected by sampling chemical compositions across the entire compositional range. The chemical compositions have been sampled by progressively changing the number of atoms per constituent by 4. For each chemical composition of binaries and ternaries, the first-principle calculations have been run for 100 randomized arrangements of the constituents on the BCC lattice sites. We collected data for a total of 3,100 randomized atomic structures over 31 chemical compositions. The calculations have been collected on Air Force HPC11 cluster using the VASP 6.5.1. Additional methodology and file structure information is available in the dataset README.txt file.

36 MATERIALS SCIENCE↗
DOE OSTI2026

TiV_BCC_SolidSolution_128atoms_VASP6

We performed density functional theory (DFT) calculations for body-centered-cubic (BCC) structures with 128 lattices sites of solid solution binary alloys titanium-vanadium (Ti-V). The electronic structures of alloys have been calculated using Vienna Ab initio Simulation Package (VASP). Within this package the DFT approach is used to reduce many-body Schrodinger equation to set of single particle Kohn-Sham (KS) equations. The generalized electronic exchange-correlation functional is described by generalized gradient approximation with the Perdew-Burke-Ernzerhof parametrization. The electron-ion interactions is described by pseudopotentials developed within the plane-wave basis projector augmented-wave (PAW) approach. These pseudopotentials are available at the VASP portal (http://cms.mpi.univie.ac.at/vasp/). Our calculations have been run with the pseudopotentials treating s and p semi-core states as valence in case for the elements Ti and V. The electronic densities and potentials are expanded over plane-waves with energy cutoff of 350 eV. 2x2x2 k-mesh and normal precision were used. The alloys were modeled by supercell containing 128 randomly distributed atoms. At initial step the atoms occupy perfect bcc lattice cites. This initial structure was optimized until energy changes less than 1e-6 eV, while forces acting on atoms don't exceed 1e-2 eV/angstrom. The electron-ion interaction is described by PAW pseudopotentials. The calculations have been collected by sampling chemical compositions across the entire compositional range. The chemical compositions have been sampled by progressively changing the number of atoms per constituent by 4. For each chemical composition of binaries and ternaries, the first-principle calculations have been run for 100 randomized arrangements of the constituents on the BCC lattice sites. We collected data for a total of 3,100 randomized atomic structures over 31 chemical compositions. The calculations have been collected on Air Force HPC11 cluster using the VASP 6.5.1. Additional methodology and file structure information is available in the dataset README.txt file.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Unraveling Strain Gradient Induced Electromechanical Coupling in Twisted Double Bilayer Graphene Moiré Superlattices

Moiré superlattices of 2D materials with a small twist angle are thought to exhibit appreciable flexoelectric effect, though unambiguous confirmation of their flexoelectricity is challenging due to artifacts associated with commonly used piezoresponse force microscopy (PFM). For example, unexpectedly small phase contrast (≈8°) between opposite flexoelectric polarizations is reported in twisted bilayer graphene (tBG), though theoretically predicted value is 180°. We report a methodology is developed to extract intrinsic moiré flexoelectricity using twisted double bilayer graphene (tDBG) as a model system, probed by lateral PFM. For small twist angle samples, it is found that a vectorial decomposition is essential to recover the small intrinsic flexoelectric response at domain walls from a large background signal. The obtained threefold symmetry of commensurate domains with significant flexoelectric response at domain walls is fully consistent with the theoretical calculations. Incommensurate domains in tDBG with relatively large twist angles can also be observed by this technique. A general strategy is provided here for unraveling intrinsic flexoelectricity in van der Waals moiré superlattices while providing insights into engineered symmetry breaking in centrosymmetric materials.

36 MATERIALS SCIENCE↗