Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “data reduction 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 883 records · Page 49
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Improving the astrometric solution of the Hyper Suprime-Cam with anisotropic Gaussian processes

Context. We study astrometric residuals from a simultaneous fit of Hyper Suprime-Cam images. Aims. We aim to characterize these residuals and study the extent to which they are dominated by atmospheric contributions for bright sources. Methods. We used Gaussian process interpolation with a correlation function (kernel) measured from the data to smooth and correct the observed astrometric residual field. Results. We find that a Gaussian process interpolation with a von Kármán kernel allows us to reduce the covariances of astrometric residuals for nearby sources by about one order of magnitude, from 30 mas 2 to 3 mas 2 at angular scales of ~1 arcmin. This also allows us to halve the rms residuals. Those reductions using Gaussian process interpolation are similar to recent result published with the Dark Energy Survey dataset. We are then able to detect the small static astrometric residuals due to the Hyper Suprime-Cam sensors effects. We discuss how the Gaussian process interpolation of astrometric residuals impacts galaxy shape measurements, particularly in the context of cosmic shear analyses at the Rubin Observatory Legacy Survey of Space and Time.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Latent Twins

Over the past decade, scientific machine learning has transformed the development of mathematical and computational frameworks for analyzing, modeling, and predicting complex systems. From inverse problems to numerical partial differential equations (PDEs), dynamical systems, and model reduction, these advances have pushed the boundaries of what can be simulated. Yet they have often progressed in parallel, with representation learning and algorithmic solution methods evolving largely as separate pipelines. With Latent Twins, we propose a unifying mathematical framework that creates a hidden surrogate in latent space for the underlying equations. Whereas digital twins mirror physical systems in the digital world, Latent Twins mirror mathematical systems in a learned latent space governed by operators. Through this lens, classical modeling, inversion, model reduction, and operator approximation all emerge as special cases of a single principle. We establish the fundamental approximation properties of Latent Twins for both ordinary differential equations (ODEs) and PDEs and demonstrate the framework across three representative settings: (i) canonical ODEs, capturing diverse dynamical regimes; (ii) a PDE benchmark using the shallow-water equations, contrasting Latent Twin simulations with deep operator network and forecasts with a four-dimensional variational method baseline; and (iii) a challenging real-data geopotential reanalysis dataset, reconstructing and forecasting from sparse, noisy observations. Latent Twins provide a compact, interpretable surrogate for solution operators that evaluate across arbitrary time gaps in a single-shot, while remaining compatible with scientific pipelines such as assimilation, control, and uncertainty quantification. Looking forward, this framework offers scalable, theory-grounded surrogates that bridge data-driven representation learning and classical scientific modeling across disciplines.

Latent Twins↗
DOE OSTI2017

Estimating parameters from rotating ring disc electrode measurements

Rotating ring disc electrode (RRDE) experiments are a classic tool for investigating kinetics of electrochemical reactions. Several standardized methods exist for extracting transport parameters and reaction rate constants using RRDE measurements. Here in this work, we compare some approximate solutions to the convective diffusion used popularly in the literature to a rigorous numerical solution of the Nernst-Planck equations coupled to the three dimensional flow problem. In light of these computational advancements, we explore design aspects of the RRDE that will help improve sensitivity of our parameter estimation procedure to experimental data. We use the oxygen reduction in acidic media involving three charge transfer reactions and a chemical reaction as an example, and identify ways to isolate reaction currents for the individual processes in order to accurately estimate the exchange current densities.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2017

Immersion Cooling of Electronics in DoD Installations

A considerable amount of energy is consumed to cool electronic equipment in data centers. A method for substantially reducing the energy needed for this cooling was demonstrated. The method involves immersing electronic equipment in a non-conductive liquid that changes phase from a liquid to a gas. The liquid used was 3M Novec 649. Two-phase immersion cooling using this liquid is not viable at this time. The primary obstacles are IT equipment failures and costs. However, the demonstrated technology met the performance objectives for energy efficiency and greenhouse gas reduction. Before commercialization of this technology can occur, a root cause analysis of the failures should be completed, and the design changes proven.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2024

TCO Analysis Approach and Regional Analysis of dWPT for Class 8 Tractors

Dynamic Wireless Power Transfer (dWPT) is a method by which battery electric vehicles (BEVs) can charge their battery while traveling on the road without the need for a physical conductive connection to the power source. dWPT has been proposed as a strategy to enable a reduction in vehicle battery capacity and associated mass and cost. In this slide deck presented at the EVs@Scale Consortium - Wireless Power Transfer Pillar Deep-Dive Meeting on November 11th, 2023, NREL provides results from an evaluation of dWPT using data from Class 8 tractors driving in the Atlanta Metro Area. NREL selected data for archetypal days representing local, regional, and long-haul trips, defined according to trip length, that included travel on primary roadways. EVI-InMotion (Electric Vehicle Infrastructure - InMotion), a systems planning and optimization tool developed at NREL, was used to evaluate dWPT performance assuming dWPT charging on 120 road segments for a total roadway lane distance of 2,365 miles. The EVI-InMotion results and representative day drive cycles were analyzed with NREL's T3CO (Transportation Technology Total Cost of Ownership) tool to estimate the total cost of ownership (TCO) for scenarios comprising two model years - 2030 and 2040 - and two technology progress cases. TCO was calculated for diesel, fuel cell electric, BEVs with batteries sized assuming no dWPT capabilities, and 200kWh BEVs with dWPT installed. This analysis finds that en-route stationary charging frequency and downtime when not on electrified roadways are the main contributors to TCO for the dWPT vehicles and that these vehicles can achieve cost parity with FCEVs at low electricity costs. Based on the scenario assumptions used here, low electricity costs would further help the cost parity with diesel vehicles in regional and long-haul cases due to stationary fueling downtime. This presentation also concludes that key factors affecting the parity potential of dWPT-capable vehicles include more extensive dWPT road coverage, higher en-route charging power, less expensive power batteries, and higher hydrogen or diesel fuel costs.

ADVANCED PROPULSION SYSTEMS,ENERGY PLANNING, POLIC↗
DOE OSTI2023

Studying the hadron structure with PANDA and CLAS using machine learning techniques

The hadron spectroscopy and structure are currently very active fields of research to study the non-perturbative regime of quantum chronodynamics. The first one studies the complex structure of excited hadrons by looking at their decay products, while the latter uses lepton scattering on nucleons. Both methods require reconstruction algorithms with great efficiency and good particle identification and background rejection rates. This work aims to provide these by either improving the existing methods or developing new ones. The first part of this document presents a feasibility study of a predicted hybrid charmonium state for the PANDA experiment. Lattice QCD calculations predict the ground state hybrid charmonium to be a spin exotic with quantum numbers of JP C = 1?+ at a mass of around 4.3 GeV with a width to be around 20 MeV. A machine learning based data analysis scheme is proposed to further improve the signal efficiency and the background reduction, alongside with improvements of the analysis software (PandaRoot), that are vital for this study. These improvements include a reworked clustering algorithm for the electromagnetic calorimeter (EMC) and an optimized monte carlo matching for neutral particles. The second part of this document is about studying the proton structure. A multidimensional study of the structure function ratio Fsin(?)LU /FUU has been performed for K±, based on the measurement of beam-spin asymmetries. It uses the high statistics data recorded with the CLAS12 spectrometer at Jefferson Laboratory. Fsin(?)LU is a twist-3 quantity that provides information about the quark gluon correlations in the proton. This document will present for the first time a simultaneous analysis of two kaon channels over a large kinematic range of z, xB , PT and Q2 with virtualities Q2 ranging from 1 GeV2 up to 8 GeV2 using machine learning techniques for improved particle identification.

Kripko, Aron↗
DOE OSTI2021

Updated systematics of intermediate-energy single-nucleon removal cross sections

The body of experimental measurements of intermediate-energy reactions that remove a single nucleon from a secondary beam of neutron- or proton-rich nuclei continues to grow. These data have been analyzed consistently using an approximate, eikonal-model treatment of the reaction dynamics combined with appropriate shell-model descriptions of the projectile initial state, the bound final states spectrum of the reaction residue, and single-particle removal strengths computed from their wave-function overlaps. Here, the systematics of the ratio R s of the measured inclusive cross section to all bound final states and the calculated cross section to bound shell-model states—in different regions of the nuclear chart and involving both very weakly bound and strongly bound valence nucleons—is important in relating the empirically deduced orbital occupancies to those from the best available shell-model predictions. Importantly, several new higher-energy measurements, for which the sudden-approximation aspect of the dynamical description is placed on an even stronger footing, now supplement the previously analyzed measurements. These additional data sets are discussed. Their R s values are shown to conform to and reinforce the earlier-observed systematics, with no indication that the approximately linear reduction in R s with increasing nucleon separation energy is a consequence of a breakdown of the sudden approximation.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2013

Magnetic inhibition of convection and the fundamental properties of low-mass stars. I. Stars with a radiative core

Magnetic fields are hypothesized to inflate the radii of low-mass stars—defined as less massive than 0.8 M {sub ☉}—in detached eclipsing binaries (DEBs). We investigate this hypothesis using the recently introduced magnetic Dartmouth stellar evolution code. In particular, we focus on stars thought to have a radiative core and convective outer envelope by studying in detail three individual DEBs: UV Psc, YY Gem, and CU Cnc. Our results suggest that the stabilization of thermal convection by a magnetic field is a plausible explanation for the observed model-radius discrepancies. However, surface magnetic field strengths required by the models are significantly stronger than those estimated from observed coronal X-ray emission. Agreement between model predicted surface magnetic field strengths and those inferred from X-ray observations can be found by assuming that the magnetic field sources its energy from convection. This approach makes the transport of heat by convection less efficient and is akin to reduced convective mixing length methods used in other studies. Predictions for the metallicity and magnetic field strengths of the aforementioned systems are reported. We also develop an expression relating a reduction in the convective mixing length to a magnetic field strength in units of the equipartition value. Our results are compared with those from previous investigations to incorporate magnetic fields to explain the low-mass DEB radius inflation. Finally, we explore how the effects of magnetic fields might affect mass determinations using asteroseismic data and the implication of magnetic fields on exoplanet studies.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

SU-F-J-161: Prostate Contouring in Patients with Bilateral Hip Prostheses: Impact of Using Artifact-Reduced CT Images and MRI

Purpose: To investigate the usefulness of metal artifact reduction in CT images of patients with bilateral hip prostheses (BHP) for contouring the prostate and determine if the inclusion of MR images provides additional benefits. Methods: Five patients with BHP were CT scanned using our clinical protocol (140kV, 300mAs, 3mm slices, 1.5mm increment, Philips Medical Systems, OH). Images were reconstructed with the orthopaedic metal artifact reduction (O-MAR) algorithm. MRI scanning was then performed (1.5T, GE Healthcare, WI) with a flat table-top (T{sub 2}-weighted, inherent body coil, FRFSE, 3mm slices with 0mm gap). All images were transferred to Pinnacle (Version 9.2, Philips Medical Systems). For each patient, two data sets were produced: one containing the O-MAR-corrected CT images and another containing fused MRI and O-MAR-corrected CT images. Four genito-urinary radiation oncologists contoured the prostate of each patient on the O-MAR-corrected CT data. Two weeks later, they contoured the prostate on the fused data set, blinded to all other contours. During each contouring session, the oncologists reported their confidence in the contours (1=very confident, 3=not confident) and the contouring difficulty that they experienced (1=really easy, 4=very challenging). Prostate volumes were computed from the contours and the conformity index was used to evaluate inter-observer variability. Results: Larger prostate volumes were found on the O-MAR-corrected CT set than on the fused set (p< 0.05, median=36.9cm{sup 3} vs. 26.63 cm{sup 3}). No significant differences were noted in the inter-observer variability between the two data sets (p=0.3). Contouring difficulty decreased with the addition of MRI (p<0.05) while the radiation oncologists reported more confidence in their contours when MRI was fused with the O-MAR-corrected CT data (p<0.05). Conclusion: This preliminary work demonstrated that, while O-MAR correction to CT images improves visualization of anatomy, the addition of MRI enhanced the oncologists’ confidence in contouring the prostate in patients with BHP.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2019

Fission Gamma-Ray Spectra Relevant to Post-Detonation Debris Analysis

A terrorist-generated nuclear attack would invoke unprecedented pressure for a rapid and accurate nuclear forensics characterization of the event. In-situ fission product gamma-ray spectroscopy offers the potential for a significant reduction in analysis time while still providing high fidelity identification of the primary composition of the detonated device. The in-situ spectroscopic approach we are proposing can rapidly identify the species that underwent fission as well as the fission-inducing neutron energy spectrum. To demonstrate this goal, we have used the gammaray history from our previously obtained LLNL-LANL-TUNL fission product-yield (FPY) data to provide insights for post-detonation debris analysis. The proposed method relies on the ratio of photopeak areas for selected gamma-rays representing the FPYs of sixteen fission products generated using mono-energetic neutrons at energies between 0.5 MeV to 14.8 MeV on highlyenriched 239 Pu, 235 U, and 238 U sealed targets. In addition to the fission product ratio measurements we have applied a new regression analysis based on our experimental FPY data and validated the approach using relevant Godiva FPY data and (surrogate) photofission product yields of 240 Pu. The main advantage of the regression analysis is that it investigates the FPY data set as a whole, rather than focusing on individual FPY ratios.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2014

Optimizing mini-ridge filter thickness to reduce proton treatment times in a spot-scanning synchrotron system

Purpose: Study the contributors to treatment time as a function of Mini-Ridge Filter (MRF) thickness to determine the optimal choice for breath-hold treatment of lung tumors in a synchrotron-based spot-scanning proton machine. Methods: Five different spot-scanning nozzles were simulated in TOPAS: four with MRFs of varying maximal thicknesses (6.15–24.6 mm) and one with no MRF. The MRFs were designed with ridges aligned along orthogonal directions transverse to the beam, with the number of ridges (4–16) increasing with MRF thickness. The material thickness given by these ridges approximately followed a Gaussian distribution. Using these simulations, Monte Carlo data were generated for treatment planning commissioning. For each nozzle, standard and stereotactic (SR) lung phantom treatment plans were created and assessed for delivery time and plan quality. Results: Use of a MRF resulted in a reduction of the number of energy layers needed in treatment plans, decreasing the number of synchrotron spills needed and hence the treatment time. For standard plans, the treatment time per field without a MRF was 67.0 ± 0.1 s, whereas three of the four MRF plans had treatment times of less than 20 s per field; considered sufficiently low for a single breath-hold. For SR plans, the shortest treatment time achieved was 57.7 ± 1.9 s per field, compared to 95.5 ± 0.5 s without a MRF. There were diminishing gains in time reduction as the MRF thickness increased. Dose uniformity of the PTV was comparable across all plans; however, when the plans were normalized to have the same coverage, dose conformality decreased with MRF thickness, as measured by the lung V20%. Conclusions: Single breath-hold treatment times for plans with standard fractionation can be achieved through the use of a MRF, making this a viable option for motion mitigation in lung tumors. For stereotactic plans, while a MRF can reduce treatment times, multiple breath-holds would still be necessary due to the limit imposed by the proton extraction time. To balance treatment time and normal tissue dose, the ideal MRF choice was shown to be the thinnest option that is able to achieve the desired breath-hold timing.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

Technical Note: Phantom study to evaluate the dose and image quality effects of a computed tomography organ-based tube current modulation technique

Purpose: This technical note quantifies the dose and image quality performance of a clinically available organ-dose-based tube current modulation (ODM) technique, using experimental and simulation phantom studies. The investigated ODM implementation reduces the tube current for the anterior source positions, without increasing current for posterior positions, although such an approach was also evaluated for comparison. Methods: Axial CT scans at 120 kV were performed on head and chest phantoms on an ODM-equipped scanner (Optima CT660, GE Healthcare, Chalfont St. Giles, England). Dosimeters quantified dose to breast, lung, heart, spine, eye lens, and brain regions for ODM and 3D-modulation (SmartmA) settings. Monte Carlo simulations, validated with experimental data, were performed on 28 voxelized head phantoms and 10 chest phantoms to quantify organ dose and noise standard deviation. The dose and noise effects of increasing the posterior tube current were also investigated. Results: ODM reduced the dose for all experimental dosimeters with respect to SmartmA, with average dose reductions across dosimeters of 31% (breast), 21% (lung), 24% (heart), 6% (spine), 19% (eye lens), and 11% (brain), with similar results for the simulation validation study. In the phantom library study, the average dose reduction across all phantoms was 34% (breast), 20% (lung), 8% (spine), 20% (eye lens), and 8% (brain). ODM increased the noise standard deviation in reconstructed images by 6%–20%, with generally greater noise increases in anterior regions. Increasing the posterior tube current provided similar dose reduction as ODM for breast and eye lens, increased dose to the spine, with noise effects ranging from 2% noise reduction to 16% noise increase. At noise equal to SmartmA, ODM increased the estimated effective dose by 4% and 8% for chest and head scans, respectively. Increasing the posterior tube current further increased the effective dose by 15% (chest) and 18% (head) relative to SmartmA. Conclusions: ODM reduced dose in all experimental and simulation studies over a range of phantoms, while increasing noise. The results suggest a net dose/noise benefit for breast and eye lens for all studied phantoms, negligible lung dose effects for two phantoms, increased lung dose and/or noise for eight phantoms, and increased dose and/or noise for brain and spine for all studied phantoms compared to the reference protocol.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

Investigation of practical initial attenuation image estimates in TOF-MLAA reconstruction for PET/MR

Purpose: Time-of-flight joint attenuation and activity positron emission tomography reconstruction requires additional calibration (scale factors) or constraints during or post-reconstruction to produce a quantitative μ-map. In this work, the impact of various initializations of the joint reconstruction was investigated, and the initial average mu-value (IAM) method was introduced such that the forward-projection of the initial μ-map is already very close to that of the reference μ-map, thus reducing/minimizing the offset (scale factor) during the early iterations of the joint reconstruction. Consequently, the accuracy and efficiency of unconstrained joint reconstruction such as time-of-flight maximum likelihood estimation of attenuation and activity (TOF-MLAA) can be improved by the proposed IAM method. Methods: 2D simulations of brain and chest were used to evaluate TOF-MLAA with various initial estimates which include the object filled with water uniformly (conventional initial estimate), bone uniformly, the average μ-value uniformly (IAM magnitude initialization method), and the perfect spatial μ-distribution but with a wrong magnitude (initialization in terms of distribution). 3D GATE simulation was also performed for the chest phantom under a typical clinical scanning condition, and the simulated data were reconstructed with a fully corrected list-mode TOF-MLAA algorithm with various initial estimates. The accuracy of the average μ-values within the brain, chest, and abdomen regions obtained from the MR derived μ-maps was also evaluated using computed tomography μ-maps as the gold-standard. Results: The estimated μ-map with the initialization in terms of magnitude (i.e., average μ-value) was observed to reach the reference more quickly and naturally as compared to all other cases. Both 2D and 3D GATE simulations produced similar results, and it was observed that the proposed IAM approach can produce quantitative μ-map/emission when the corrections for physical effects such as scatter and randoms were included. The average μ-value obtained from MR derived μ-map was accurate within 5% with corrections for bone, fat, and uniform lungs. Conclusions: The proposed IAM-TOF-MLAA can produce quantitative μ-map without any calibration provided that there are sufficient counts in the measured data. For low count data, noise reduction and additional regularization/rescaling techniques need to be applied and investigated. The average μ-value within the object is prior information which can be extracted from MR and patient database, and it is feasible to obtain accurate average μ-value using MR derived μ-map with corrections as demonstrated in this work.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2014

Harmonic analysis for the characterization and correction of geometric distortion in MRI

Purpose: Magnetic resonance imaging (MRI) is gaining widespread use in radiation therapy planning, patient setup verification, and real-time guidance of radiation delivery. Successful implementation of these technologies relies on the development of simple and efficient methods to characterize and monitor the geometric distortions arising due to system imperfections and gradient nonlinearities. To this end, the authors present the theory and validation of a novel harmonic approach to the quantification of system-related distortions in MRI. Methods: The theory of spatial encoding in MRI is applied to demonstrate that the 3D distortion vector field (DVF) is given by the solution of a second-order boundary value problem (BVP). This BVP is comprised of Laplace’s equation and a limited measurement of the distortion on the boundary of a specified region of interest (ROI). An analytical series expansion solving this BVP within a spherical ROI is obtained, and a statistical uncertainty analysis is performed to determine how random errors in the boundary measurements propagate to the ROI interior. This series expansion is then evaluated to obtain volumetric DVF mappings that are compared to reference data obtained on a 3 T full-body scanner. This validation is performed within two spheres of 20 cm diameter (one centered at the scanner origin and the other offset +3 cm along each of the transverse directions). Initially, a high-order mapping requiring measurements at 5810 boundary points is used. Then, after exploring the impact of the boundary sampling density and the effect of series truncation, a reduced-order mapping requiring measurements at 302 boundary points is evaluated. Results: The volumetric DVF mappings obtained from the harmonic analysis are in good agreement with the reference data. Following distortion correction using the high-order mapping, the authors estimate a reduction in the mean distortion magnitude from 0.86 to 0.42 mm and from 0.93 to 0.39 mm within the central and offset ROIs, respectively. In addition, the fraction of points with a distortion magnitude greater than 1 mm is reduced from 35.6% to 2.8% and from 40.4% to 1.5%, respectively. Similarly, following correction using the reduced-order mapping, the mean distortion magnitude reduces to 0.45–0.42 mm within the central and offset ROIs, and the fraction of points with a distortion magnitude greater than 1 mm is reduced to 2.8% and 1.5%, respectively. Conclusions: A novel harmonic approach to the characterization of system-related distortions in MRI is presented. This method permits a complete and accurate mapping of the DVF within a specified ROI using a limited measurement of the distortion on the ROI boundary. This technique eliminates the requirement to exhaustively sample the DVF at a dense 3D array of points, thereby permitting the design of simple, inexpensive phantoms that may incorporate additional modules for auxiliary quality assurance objectives.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2023

Fast and accurate reduced-order modeling of a MOOSE-based additive manufacturing model with operator learning

One predominant challenge in additive manufacturing (AM) is to achieve specific material properties by manipulating manufacturing process parameters during the runtime. Such manipulation tends to increase the computational load imposed on existing simulation tools employed in AM. The goal of the present work is to construct a fast and accurate reduced-order model (ROM) for an AM model developed within the Multiphysics Object-Oriented Simulation Environment (MOOSE) framework, ultimately reducing the time/cost of AM control and optimization processes. Our adoption of the operator learning (OL) approach enabled us to learn a family of differential equations produced by altering process variables in the laser’s Gaussian point heat source. More specifically, we used the Fourier neural operator (FNO) and deep operator network (DeepONet) to develop ROMs for time-dependent responses. Furthermore, we benchmarked the performance of these OL methods against a conventional deep neural network (DNN)-based ROM. Ultimately, we found that OL methods offer comparable performance and, in terms of accuracy and generalizability, even outperform DNN at predicting scalar model responses. The DNN-based ROM afforded the fastest training time. Furthermore, all the ROMs were faster than the original MOOSE model yet still provided accurate predictions. FNO had a smaller mean prediction error than DeepONet, with a larger variance for time-dependent responses. Unlike DNN, both FNO and DeepONet were able to simulate time series data without the need for dimensionality reduction techniques. Finally, the present work can help facilitate the AM optimization process by enabling faster execution of simulation tools while still preserving evaluation accuracy.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Mapping structural heterogeneity at the nanoscale with scanning nano-structure electron microscopy (SNEM)

Here, in this work, we explore the use of scanning electron diffraction (also known as 4D-STEM) coupled with electron atomic pair distribution function analysis (ePDF) to understand the local order (structure and chemistry) as a function of position in a complex multicomponent system, a hot rolled, Ni-encapsulated, Zr 65 Cu 17.5 Ni 10 Al 7.5 bulk metallic glass (BMG), with a spatial resolution of 3 nm. We show that it is possible to gain insight into the chemistry and chemical clustering/ordering tendency in different regions of the sample, including in the vicinity of nano-scale crystallites that are identified from virtual dark field images and in heavily deformed regions at the edge of the BMG. In addition to simpler analysis, unsupervised machine learning was used to extract partial PDFs from the material, modeled as a quasi-binary alloy, and map them in space. These maps allowed key insights not only into the local average composition, as validated by EELS, but also a unique insight into chemical short-range ordering tendencies in different regions of the sample during formation. The experiments are straightforward and rapid and, unlike spectroscopic measurements, don’t require energy filters on the instrument. We spatially map different quantities of interest (QoI’s), defined as scalars that can be computed directly from positions and widths of ePDF peaks or parameters refined from fits to the patterns. We developed a flexible and rapid data reduction and analysis software framework that allows experimenters to rapidly explore images of the sample on the basis of different QoI’s. The power and flexibility of this approach are explored and described in detail. Because of the fact that we are getting spatially resolved images of the nanoscale structure obtained from ePDFs we call this approach scanning nano-structure electron microscopy (SNEM), and we believe that it will be powerful and useful extension of current 4D-STEM methods.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Data-driven models of nonautonomous systems

Nonautonomous dynamical systems are characterized by time-dependent inputs, which complicates the discovery of predictive models describing the spatiotemporal evolution of the state variables of quantities of interest from their temporal snapshots. When dynamic mode decomposition (DMD) is used to infer a linear model, this difficulty manifests itself in the need to approximate the time-dependent Koopman operators. Our approach is to approximate the original nonautonomous system with a modified system derived via a local parameterization of the time-dependent inputs. The modified system comprises a sequence of local parametric systems, which are subsequently approximated by a parametric surrogate model using the DRIPS (dimension reduction and interpolation in parameter space) framework. The offline step of DRIPS relies on DMD to build a linear surrogate model, endowed with reduced-order bases for the observables mapped from training data. The online step interpolates on suitable manifolds to construct a sequence of iterative parametric surrogate models; the target/test parameter points on these manifolds are specified by a local parameterization of the test time-dependent inputs. Here, we use numerical experimentation to demonstrate the robustness of our method and compare its performance with that of deep neural networks.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2014

Improving Thermal Model Prediction Through Statistical Analysis of Irradiation and Post-Irradiation Data from AGR Experiments

As part of the High Temperature Reactors (HTR) R&D program, a series of irradiation tests, designated as Advanced Gas-cooled Reactor (AGR), have been defined to support development and qualification of fuel design, fabrication process, and fuel performance under normal operation and accident conditions. The AGR tests employ fuel compacts placed in a graphite cylinder shrouded by a steel capsule and instrumented with thermocouples (TC) embedded in graphite blocks enabling temperature control. While not possible to obtain by direct measurements in the tests, crucial fuel conditions (e.g., temperature, neutron fast fluence, and burnup) are calculated using core physics and thermal modeling codes. This paper is focused on AGR test fuel temperature predicted by the ABAQUS code's finite element-based thermal models. The work follows up on a previous study, in which several statistical analysis methods were adapted, implemented in the NGNP Data Management and Analysis System (NDMAS), and applied for qualification of AGR-1 thermocouple data. Abnormal trends in measured data revealed by the statistical analysis are traced to either measuring instrument deterioration or physical mechanisms in capsules that may have shifted the system thermal response. The main thrust of this work is to exploit the variety of data obtained in irradiation and post-irradiation examination (PIE) for assessment of modeling assumptions. As an example, the uneven reduction of the control gas gap in Capsule 5 found in the capsule metrology measurements in PIE helps identify mechanisms other than TC drift causing the decrease in TC readings. This suggests a more physics-based modification of the thermal model that leads to a better fit with experimental data, thus reducing model uncertainty and increasing confidence in the calculated fuel temperatures of the AGR-1 test.

99 GENERAL AND MISCELLANEOUS↗