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,279 records · Page 71

MimicGAN: Robust Projection onto Image Manifolds with Corruption Mimicking

In the past few years, Generative Adversarial Networks (GANs) have dramatically advanced our ability to represent and parameterize high-dimensional, non-linear image manifolds. As a result, they have been widely adopted across a variety of applications, ranging from challenging inverse problems like image completion, to problems such as anomaly detection and adversarial defense. A recurring theme in many of these applications is the notion of projecting an image observation onto the manifold that is inferred by the generator. In this context, Projected Gradient Descent (PGD) has been the most popular approach, which essentially optimizes for a latent vector that minimizes the discrepancy between a generated image and the given observation. However, PGD is a brittle optimization technique that fails to identify the right projection (or latent vector) when the observation is corrupted, or perturbed even by a small amount. Such corruptions are common in the real world, for example images in the wild come with unknown crops, rotations, missing pixels, or other kinds of non-linear distributional shifts which break current encoding methods, rendering downstream applications unusable. To address this, we propose corruption mimicking—a new robust projection technique, that utilizes a surrogate network to approximate the unknown corruption directly at test time, without the need for additional supervision or data augmentation. The proposed method is significantly more robust than PGD and other competing methods under a wide variety of corruptions, thereby enabling a more effective use of GANs in real-world applications. Finally, more importantly, we show that our approach produces state-of-the-art performance in several GAN-based applications—anomaly detection, domain adaptation, and adversarial defense, that benefit from an accurate projection.

97 MATHEMATICS AND COMPUTING↗

Comparison of EBSD, DIC, AFM, and ECCI for active slip system identification in deformed Ti-7Al

The use of SEM-DIC, AFM, ECCI, and HR-EBSD to characterize slip-system activity was assessed on the same material volume of Ti-7Al. This study presents a robust comparison of the various methods for the first time, including an assessment of their advantages and disadvantages, and how they can be used effectively in a complementary fashion. The analysis of the different approaches was carried out in a blind, round-robin manner at three different universities. A Ti-7Al specimen was deformed in uniaxial tension to approximately 3% axial strain, and the active slip systems were independently identified using (i) trace analysis; (ii) in-SEM digital image correlation, (iii) observations of residual dislocations from ECCI, and (iv) long-range rotation gradients through HR-EBSD, with consistent trace identification in all cases. Displacement data from AFM was used to augment SEM-DIC displacement data by providing complementary out-of-plane displacement information. Furthermore, short-range dislocation gradients (measured by DIC) provided insight into the residual geometrically necessary dislocation (GND) content, and was consistent with the GND content extracted from EBSD data and ECCI images, confirming the presence of residual GNDs on the dominant slip systems resulting in visible slip bands. Furthermore, these approaches can be used in tandem to provide multi-modal information on slip band identification, strain and orientation gradients, out-of-plane displacements, and the presence of GNDs and SSDs, all of which can be used to inform and validate the development of dislocation-based crystal plasticity and strain gradient models.

36 MATERIALS SCIENCE↗

Joint Inversion of Body Wave Arrival Times, Surface Wave Dispersion Data and Receiver Functions: Method and Application to South China

In this work, we have developed a new joint inversion method that incorporates body wave arrival times, surface wave dispersion and receiver functions to simultaneously update earthquake locations and constrain three-dimensional P-wave (Vp) and S-wave velocity (Vs) models. Due to complementary sensitivities of the three types of data, the proposed joint inversion algorithm can reduce the intrinsic non-uniqueness of inversions using fewer types of data and better determine smooth velocity variations and velocity discontinuities. Synthetic tests demonstrate the advantages of this new joint inversion algorithm in resolving velocity structures, especially in constraining velocity gradients across the Moho interface. We have applied the proposed joint inversion algorithm to image the lithosphere velocity structure of south China. The inverted Vp and Vs models fit body wave arrival times, surface wave dispersion and receiver functions well. We further analyzed the distribution of Vs gradients across the Moho interface in detail, which helps us better understand tectonics in south China.

58 GEOSCIENCES↗

Throughfall-Reduction Drying Effects on Soil Respiration, Moisture, and Temperature in Four Distinct Lowland Panamanian Forests from 2015-2022

Objectives: Climatic drying is predicted for many tropical forests, yet effects on soil properties across moisture and soil gradients within tropical forests remain poorly characterized, hampering predictions of forest-climate feedbacks. We hypothesized that throughfall reduction (i.e. chronic drying) would suppress soil CO2 fluxes (i.e., respiration) in already-drier tropical forests by further reductions in soil moisture, but increases CO2 fluxes in wetter tropical forests by alleviating anaerobiosis and soil saturation. We conducted a field study in four distinct lowland Panamanian forests that vary in background rainfall and soil fertility. In each forest we constructed throughfall reduction structures that reduce moisture inputs to plots by ~50% (see methods). In drying and paired control plots we measured soil CO2 fluxes, soil moisture (volume/volume), soil temperature (°C), air temperature (°C), and forest floor biomass during wet-dry cycles (2015 – 2022). We also collected litterfall biomass and surveyed all tree species for diameter at breast height (dbh), and identified tree species in 2018 and 2019.Results:We found that soil moisture peaked in the wet season and declined in the dry season. Measured soil CO2 fluxes declined in the dry season and peaked in the early wet season ahead of peak soil moisture, resulting in a lower soil moisture optimum for respiration than previously modeled. Chronic throughfall exclusion also suppressed soil moisture across the four forests, and also initially suppressed soil CO2 fluxes across forests. There was sustained suppression of soil CO2 fluxes after four years in the wettest forest only (-28 ± 4% during the dry season), but elevated soil CO2 fluxes in a fertile forest after four years (+75 ± 28% during the late wet season). The unexpected negative drying effect in the wettest, most infertile forest could have resulted from reduced vertical flushing of nutrients into soils, as the drying effect increased with time. Including hydro-nutrient interactions in ecosystem models could improve predictions of tropical forest-climate feedbacks. Datasets included: Datasets included here include .csv and .xls files for soil respiration (carbon dioxide flux) collected using the LiCor (see methods), air temperature °C , soil moisture (vl/vl), and soil temperature °C. Soil moisture and soil temperature were collected at 0-10cm depths using hand-held probes (see methods). Data are on ~3 week timescales from 2015-2022 with some gaps. There is also a .kml file that includes coordinates for all 32 plots included in the study of four forests (n = 4 throughfall reduction and n = 4 control plots per site). No special software is needed to open these files.

54 ENVIRONMENTAL SCIENCES↗

Quantitative Insight to Fission Gas Bubble Distribution in Irradiated Annular U-10Zr Metallic Fuel Using Deep Learning

U-10Zr Metal fuel is a promising nuclear fuel candidate for next-generation sodium-cooled fast spectrum reactors. Porosity is one of the key facts to impact the performance of metallic fuel. Additionally, a mechanical understanding of fission gas bubbles evolution behavior is a prerequisite for fuel development and qualification. Previous study of fission gas bubbles relied on a simple threshold method working on low resolution optical microscopy images, which has challenges in recognizing bubble boundaries, and caused inaccurate statistics of bubble properties. In this paper, a pre-trained deep learning model on Scanning Electron Microscopy (SEM) images from an annular U-10Zr fuel (AF1), was applied to another U-10Zr annular fuel (AF2). More accurate fission gas bubble segmentation results were generated, which leads to more precise qualitative analysis on the morphology, size, density, and orientation of bubbles. Furthermore, we investigated the lanthanide movement along the radial temperature gradient and obtained conclusive findings.

36 MATERIALS SCIENCE↗

Prediction and uncertainty quantification of shale well performance using multifidelity Monte Carlo

Uncertainty quantification is an integral component of reservoir management, especially considering the inherent uncertainty in subsurface systems. While a standard practice to estimate the uncertainty, Monte Carlo (MC) simulation is computationally intense when the sampling population comprises high-fidelity simulations. Alternatively, the Multi-fidelity Monte Carlo (MFMC) simulation overcomes this computational intensity by integrating low- and high-fidelity simulations. Our goal is to minimize the number of expensive high-fidelity simulations while maintaining accuracy and using numerous fast and cheap low-fidelity simulations to efficiently sample to input parameter space of interest. We selected gas production from unconventional wells to demonstrate the potential speedups and accuracy of the MFMC approach. The model fidelity usually determines the trade-off between accuracy and efficiency. While the high-fidelity model is more accurate, the low-fidelity model is more efficient. Our high-fidelity simulation includes reservoir simulations of a hydraulically fractured well. On the other hand, our low-fidelity model comprises the parallel-plate flow model. We used differential programming to efficiently solve the 1D flow model, where automatic differentiation is used to efficiently compute the gradients. We matched the production profile of high-fidelity simulations with our low-fidelity simulations. Then, we used a support vector regression to map the high- and low-fidelity input parameters. The mapping function is essential to tune the low-dimensional parameter space of the low-fidelity model to the high-dimensional parameter space of the high-fidelity model. We found that we can use a combination of 9 high fidelity and 10,000 low fidelity simulations to efficiently and accurately simulate pressure management. This method is at least two orders of magnitude faster than only using high-fidelity simulations. Finally, from a broader perspective, MFMC could efficiently estimate the uncertainty of various systems and models, integrating low- and high-fidelity models.

04 OIL SHALES AND TAR SANDS↗

Soil gas probes for monitoring trace gas messengers of microbial activity

Abstract Soil microbes vigorously produce and consume gases that reflect active soil biogeochemical processes. Soil gas measurements are therefore a powerful tool to monitor microbial activity. Yet, the majority of soil gases lack non-disruptive subsurface measurement methods at spatiotemporal scales relevant to microbial processes and soil structure. To address this need, we developed a soil gas sampling system that uses novel diffusive soil probes and sample transfer approaches for high-resolution sampling from discrete subsurface regions. Probe sampling requires transferring soil gas samples to above-ground gas analyzers where concentrations and isotopologues are measured. Obtaining representative soil gas samples has historically required balancing disruption to soil gas composition with measurement frequency and analyzer volume demand. These considerations have limited attempts to quantify trace gas spatial concentration gradients and heterogeneity at scales relevant to the soil microbiome. Here, we describe our new flexible diffusive probe sampling system integrated with a modified, reduced volume trace gas analyzer and demonstrate its application for subsurface monitoring of biogeochemical cycling of nitrous oxide (N 2 O) and its site-specific isotopologues, methane, carbon dioxide, and nitric oxide in controlled soil columns. The sampling system observed reproducible responses of soil gas concentrations to manipulations of soil nutrients and redox state, providing a new window into the microbial response to these key environmental forcings. Using site-specific N 2 O isotopologues as indicators of microbial processes, we constrain the dynamics of in situ microbial activity. Unlocking trace gas messengers of microbial activity will complement -omics approaches, challenge subsurface models, and improve understanding of soil heterogeneity to disentangle interactive processes in the subsurface biome.

04 OIL SHALES AND TAR SANDS↗

Augmenting a Simulation Campaign for Hybrid Computer Model and Field Data Experiments

The Kennedy and O’Hagan (KOH) calibration framework uses coupled Gaussian processes (GPs) to meta-model an expensive simulator (first GP), tune its “knobs” (calibration inputs) to best match observations from a real physical/field experiment and correct for any modeling bias (second GP) when predicting under new field conditions (design inputs). There are well-established methods for placement of design inputs for data-efficient planning of a simulation campaign in isolation, that is, without field data: space-filling, or via criterion like minimum integrated mean-squared prediction error (IMSPE). Analogues within the coupled GP KOH framework are mostly absent from the literature. Here, in this study, we derive a closed form IMSPE criterion for sequentially acquiring new simulator data for KOH. We illustrate how acquisitions space-fill in design space, but concentrate in calibration space. Closed form IMSPE precipitates a closed-form gradient for efficient numerical optimization. We demonstrate that our KOH-IMSPE strategy leads to a more efficient simulation campaign on benchmark problems, and conclude with a showcase on an application to equilibrium concentrations of rare earth elements for a liquid–liquid extraction reaction.

97 MATHEMATICS AND COMPUTING↗

Virtual Agents-Based Attack-Resilient Distributed Control for Islanded AC Microgrid

Due to its dependence on a communication network, distributed secondary control of microgrids is susceptible to denial-of-service (DoS) attacks in channel shutdown mode, which may negatively impact the network connectivity and thus deteriorate the coordination and power sharing among distributed generators (DGs). Honeypot is a common method for cyber deception by introducing fake targets. However, in the context of microgrid, the misleading information spread by honeypots will also impact the system performance. This paper proposes an attack-resilient distributed control for AC microgrids utilizing virtual agents (VAs) to counteract both DoS edge and node attacks. The VAs are designed to not impact the system’s steady state during normal operation but to share information among neighboring real agents and serve as dummy targets for DoS attacks. The control with VAs is implemented by a primal-dual gradient based distributed algorithm to efficiently obtain a practical solution for voltage/frequency regulation and power sharing. The simulation results on a 4-DG test system and a modified IEEE 34-bus system show that 1) VAs do not impact the normal functionality of the test system, and 2) deploying VAs can enhance the resilience of the microgrid control against DoS edge and node attacks.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Controlling Oxidation of Nb in Oxygen Abundant Environments

Modern particle accelerators depend on Superconducting Radio Frequency (SRF) cavities made from high-purity niobium (Nb) to achieve optimal performance, including high quality factors and strong accelerating gradients. However, when exposed to air, niobium naturally forms a complex oxide layer that can introduce surface imperfections and carbon contamination. This work examines an alternative oxidation strategy under oxygen-rich conditions to better regulate the oxide formation process. The ultimate objective is to improve surface uniformity and cleanliness, thereby reducing defect density and enhancing performance. We used Confocal Microscopy, Scanning Electron Microscopy (SEM), and X-Ray Photoelectron Spectroscopy (XPS) to analyze surface changes. We used standard metrics, like Arithmetic Average Roughness (Ra) and Root Mean Square (Rq). Three oxidation methods were applied: short, extended, and a repeated HF and H₂O₂ oxidation process. Preliminary results show promising reductions in carbon contamination and surface defects.

Romero, Juan [Fermilab]↗

A reaction–diffusion model for grayscale digital light processing 3D printing

We report that digital light processing (DLP) 3D printing is an additive manufacturing process that utilizes light patterns to photopolymerize a liquid resin into a solid. Due to the accuracy of modern digital micromirror devices (DMD) and recent advances in resin chemistry, it is now possible to create functionally graded structures using different light intensity values, also known as grayscale DLP (g-DLP). Different intensities of light lead to differences in the polymer crosslinking density after curing, which ultimately produces a part with gradients of material properties. However, g-DLP is a complicated process. First, the DLP printing is a highly coupled chemical and physical process that involves light propagation, chemical reactions, species diffusion, heat transfer, volume shrinkage, and changes in mechanical behaviors of the curing resin. Second, in g-DLP, light gradients create strong in plane gradients of chemical species concentrations in the curing liquid resin due to the strong dependence of light intensity on the rate of monomer crosslinking. Furthermore, light gradients through the depth create concentration gradients due to the degree of cure dependent light absorption and the use of photoabsorbers. These complex physical features of the printing process must be understood in order to properly control printing parameters such as light exposure time, printing speed, and grayscale variations to achieve accurate mechanical properties. In this paper, a photopolymerization reaction–diffusion model is developed and used in conjunction with experiments to investigate the coupled effects of light propagation, chemical reaction rates, and species diffusion during g-DLP 3D printing. The model is implemented numerically utilizing the finite difference method and simulation results are compared to experimental findings of simple printed structures. The agreement between experimental and model predictions of simple quantities of interest, such as geometric feature sizes, shows that the model can capture the overcure due to free-radical and other species diffusion during printing when grayscale patterns are employed. This model lays the groundwork for future extensions that can incorporate more complex coupled physics such as heat transfer, volume shrinkage, and material property evolution, which are critically important in utilizing g-DLP 3D printing for the fabrication of high-performance parts which excellent geometric and material property tolerances.

36 MATERIALS SCIENCE↗

Comparative Density Functional Theory Study of Magnetic Exchange Couplings in Dinuclear Transition-Metal Complexes

Multicenter transition-metal complexes (MCTMs) with magnetically interacting ions have been proposed as components for information-processing devices and storage units. For any practical application of MCTMs as magnetic units, it is crucial to characterize their magnetic behavior, and in particular, the isotropic magnetic exchange coupling, J, between its magnetic centers. Due to the large size of typical MCTMs, density functional theory is the only practical electronic structure method for evaluating the J coupling. Here, we assess the accuracy of different density functional approximations for predicting the magnetic couplings of eight dinuclear transition-metal complexes, including five dimanganese, two dicopper, and one divanadium with known reliable experimental J couplings spanning from ferromagnetic to strong antiferromagnetic. The density functionals considered include global hybrid functionals which mix semilocal density functional approximations and exact exchange with a fixed admixing parameter, six local hybrid functionals where the admixing parameters are extended to be spatially dependent, the SCAN and r 2 SCAN meta-generalized gradient approximations (GGAs), and two widely used GGAs. We found that global hybrids tested in this work have a tendency to over-correct the error in magnetic coupling parameters from the Perdew–Burke–Ernzerhof (PBE) GGA as seen for manganese complexes. The performance of local hybrid density functionals shows no improvement in terms of bias and is scattered without a clear trend, suggesting that more efforts are needed for the extension from global to local hybrid density functionals for this particular property. The SCAN and r 2 SCAN meta-GGAs are found to perform as well as benchmark global hybrids on most tested complexes. We further analyze the charge density redistribution of meta-GGAs as well as global and local hybrid density functionals with respect to that of PBE, in connection to the self-interaction error or delocalization error.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Fabrication of THz corrugated wakefield structure and its high power test

We present overall process for developing terahertz (THz) corrugated structure and its beam-based measurement results. 0.2-THz corrugated structures were fabricated by die stamping method as the first step demonstration towards GW THz radiation source and GV/m THz wakefield accelerator. 150-μm thick disks were produced from an OFHC (C10100) foil by stamping. Two types of disks were stacked alternately to form 46 mm structure with ~ 170 corrugations. Custom assembly was designed to provide diffusion bonding with a high precision alignment of disks. The compliance of the fabricated structure have been verified through beam-based wakefield measurement at Argonne Wakefield Accelerator Facility. Both measured longitudinal and transverse wakefield showed good agreement with simulated wakefields. Measured peak gradients, 9.4 MV/m/nC for a long single bunch and 35.4 MV/m/nC for a four bunch trains, showed good agreement with the simulation.

43 PARTICLE ACCELERATORS↗

Efficiently resolving the terrestrial-aquatic interface in E3SM with sub-grid methods to improve coastal simulations (Final Technical Report)

The broad goal of this project was to extend the capabilities of MPAS‐Ocean, the ocean component of DOE’s Energy Exascale Earth System Model (E3SM), to better resolve the hydrodynamics of the terrestrial‐aquatic interface without significantly increasing computing resource requirements. Specifically, our goal was to better resolve the small-scale features across the terrestrial‐aquatic interface which includes the low lying coastal floodplain and a dendritic network of hydraulically efficient channels connected to estuarine systems and back bays, barrier islands and high energy inlets. The processes driving flow across the terrestrial‐aquatic interface include upland and floodplain hydrology, tides, wind and atmospheric pressure (including coastally generated storm surge across continental shelves and within estuaries and bays), wind waves and water temperature and salinity gradients.

58 GEOSCIENCES↗

Machine learning in nuclear materials research

Nuclear materials are often demanded to function for extended time in extreme environments, including high radiation fluxes with associated transmutations, high temperature and temperature gradients, mechanical stresses, and corrosive coolants. They also have a wide range of microstructural and chemical makeups, resulting in multifaceted and often out-of-equilibrium interactions. Machine learning (ML) is increasingly being used to tackle these complex time-dependent interactions and aid researchers in developing models and making predictions, sometimes with better accuracy than traditional modeling that focuses on one or two parameters at a time. Conventional practices of acquiring new experimental data in nuclear materials research are often slow and expensive, limiting the opportunity for data-centric ML, but new methods are changing that paradigm. Here we review high-throughput computational and experimental data approaches, especially robotic experimentation and active learning that is based on Gaussian process and Bayesian optimization. We show ML examples in structural materials (e.g., reactor pressure vessel (RPV) alloys and radiation detecting scintillating materials) and highlight new techniques of high-throughput sample preparation and characterizations, and automated radiation/environmental exposures and real-time online diagnostics. Herein, this review suggests that ML models of material constitutive relations in plasticity, damage, and even electronic and optical responses to radiation are likely to become powerful tools as they develop. Finally, we speculate on how the recent trends of using natural language processing (NLP) to aid the collection and analysis of literature data, interpretable artificial intelligence (AI), and the use of streamlined scripting, database, workflow management, and cloud computing platforms that will soon make the utilization of ML techniques as commonplace as the spreadsheet curve-fitting practices of today.

36 MATERIALS SCIENCE↗

First steps toward predicting corrosion behavior of structural materials in molten salts

To address the need for physics-based models to predict corrosion behavior of materials in molten salts, here we propose potential methods to aid in selection of optimum materials for structural components in molten salt powered technologies. In the present work, the role of alloy thermodynamics and kinetics on governing corrosion rates of Ni-based alloys will be discussed combining experimental and computational methods. A few strategies are presented to quantify corrosion rates of Ni-based materials isothermally exposed in purified KCl-MgCl 2 salts at 600° C-800° C. The influence of capsule material, potential corrosion products and role of alloy composition on the observed corrosion rates was discussed with coupled thermodynamic-kinetic models. Larger depths of corrosion attack were observed in alloy 230 specimens compared to other alloys under similar conditions was attributed to the much higher chemical activity of Cr in the alloy that results in a larger Cr chemical potential gradient.

36 MATERIALS SCIENCE↗

Computation of Direct Sensitivities of Spatial Multibody Systems With Joint Friction

Abstract Friction exists in most mechanical systems and may have a major influence on the dynamic performance of the system. The incorporation of friction in dynamic systems has been a subject of active research for several years owing to its high nonlinearity and its dependence on several parameters. Consequently, optimization of dynamic systems with friction becomes a challenging task. Gradient-based optimization of dynamical systems is a prominent technique for optimal design and requires the computation of model sensitivities with respect to the design parameters. The novel contribution of this paper is the derivation of the analytical methodology for the computation of direct sensitivities for smooth multibody systems with joint friction using the Lagrangian index-1 formulation. System dynamics have been computed using two different friction models; the Brown and McPhee, and the Gonthier et al. model. The methodology proposed to obtain model sensitivities has also been validated using the complex finite difference method. A case study has been conducted on a spatial multibody system to observe the effect of friction on the dynamics and model sensitivities, compare sensitivities with respect to different parameters and demonstrate the numerical and validation aspects. Since design parameters can have very different magnitudes and units, the sensitivities have been scaled with the parameters for comparison. Finally, a discussion has been presented on the interpretation of the case study results. Due to the incorporation of joint friction, ‘jumps’ or discontinuities are observed in the model sensitivities akin to those observed for hybrid dynamical systems.

Engineering↗

A Novel LDPP-MADDPG Approach for Distributed Power Allocation in mmWave Cellular Networks

This paper considers the problem of distributed beam scheduling and power allocation problem in millimeter- Wave (mmWave) cellular networks, in which multiple Base Stations (BSs) operate as individual operators over a shared spectrum. We propose a novel learning-aided approach that integrates the Lyapunov Drift-Plus-Penalty (LDPP) framework and Multi-agent Deep Deterministic Policy Gradient (MADDPG) reinforcement learning algorithms. This offers a powerful approach to learning stable and constraint-aware policies, reaping the joint benefit of both LDPP and MADDPG, in complex multiagent environments. The major challenge for this approach is to integrate these two approaches in a meaningful and effective manner. The key idea to solve this problem is to introduce a novel feature of local observation that incorporates potential negative value of the reward function due to the stochastic constraints introduced by the LDPP framework. Empirical results demonstrate that our proposed scheme outperforms the baseline methods under various conditions.

99 - GENERAL AND MISCELLANEOUS↗