Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Functional Data Analysis”

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

RMCProfile7 : reverse Monte Carlo for multiphase systems

This work introduces a completely rewritten version of the programRMCProfile(version 7), big-box, reverse Monte Carlo modelling software for analysis of total scattering data. The major new feature ofRMCProfile7is the ability to refine multiple phases simultaneously, which is relevant for many current research areas such as energy materials, catalysis and engineering. Other new features include improved support for molecular potentials and rigid-body refinements, as well as multiple different data sets. An empirical resolution correction and calculation of the pair distribution function as a back-Fourier transform are now also available.RMCProfile7is freely available for download at https://rmcprofile.ornl.gov/.

Chemistry↗

Data and scripts associated with the manuscript "Using Mutual Information for Global Sensitivity Analysis on Watershed Modeling"

This package contains the data and scripts used in "Using Mutual Information for Global Sensitivity Analysis on Watershed Modeling" (Jiang et al., 2022). The ARW_SWAT.zip file contains the SWAT simulation and the sensitivity analysis result in the American River Watershed. The Portage_SWAT.zip file contains the SWAT simulation and the sensitivity analysis result in the Portage River Watershed. The src folder contains the source code of sensitivity analysis and plotting functions. The notebooks folder contains the Jupyter notebooks for performing sensitivity analysis and other post analyses.Global Sensitivity Analysis (GSA) often is applied to assess the sensitivity of model outputs to their inputs using ensemble simulations. However, increasing model complexity and the associated computational cost have limited the use of most GSA approaches for process-based watershed models. We propose to use Mutual Information (MI) as a computationally efficient GSA method for watershed modeling. Such MI computed from several hundred realizations usually can capture nonlinear relationships between inputs and outputs of interest. We perform MI-based watershed sensitivity analyses in studies of the Portage River Watershed in Ohio and the American River Watershed in Washington. In these studies, MI is used to evaluate the sensitivity of river discharges simulated by the Soil and Water Assessment Tool to no less than 20 parameters for each watershed. Our MI-based sensitivity analyses achieved convergence with about 300~500 realizations, a small fraction of the ensemble size (i.e., several thousand) required by the Sobol method. Nevertheless, MI yields similar sensitivity ranking compared to the Sobol method, especially for sensitive parameters. Our study thus sheds new light on the use of MI as an affordable GSA method for computationally intensive models such as the hyper-resolution, watershed hydrobiogeochemical models.

54 ENVIRONMENTAL SCIENCES↗

SRNL Analysis of ICCWR LCM and WAMS data for Corrosion and Cracking (FY2020 Progress Report)

The development of algorithms for machine learning and data analysis for the 3013 Surveillance Program is a collaborative effort by the Savannah River National Laboratory (SRNL) and the University of South Carolina (USC). For corrosion detection, Laser Confocal Microscope (LCM) or Wide Area 3D Measurement System (WAMS) data is extracted from large binary files, with software written to convert the data to physical attributes (e.g. height, color and grayscale values; all as functions of a location in a plane projection). It is the objective of this project to produce a user-friendly interface that incorporated all operations needed to perform surface examination. For this reason, a Matlab-based Graphical User Interface (GUI) was created to integrate data input with software developed for processing and evaluation. In summary, the GUI permits selective downloading of binary data, interrogation of attributes, data labeling, flagging of significant features, execution of Machine Learning (ML) algorithms, output of parameters for trained ML algorithms, reports of ML model accuracy with respect to labeled data, and generation of graphical representations of various analyses.

3013 corrosion↗

FY21 Progress Report: SRNL Analysis of ICCWR LCM and WAMS data for Corrosion and Cracking

The development of algorithms for machine learning and data analysis for the 3013 Surveillance Program is a collaborative effort by the Savannah River National Laboratory (SRNL) and the University of South Carolina (USC). For corrosion detection, Laser Confocal Microscope (LCM) or Wide Area 3D Measurement System (WAMS) data is extracted from large binary files, with software written to convert the data to physical attributes (e.g., height, color and grayscale values; all as functions of a location in a plane projection). A user-friendly Matlab Graphical User Interface (GUI) that reads data from either LCM or WAMS files was developed to integrate input data with software developed for processing and evaluation. The GUI can selectively download binary data, interrogate data attributes, label data, flag significant features, execute Machine Learning (ML) algorithms, output parameters for trained ML algorithms, report ML model accuracy with respect to labeled data, and generate graphical representations for various analyses. Features can be called out by user-specified thresholds, manual labeling or machine learning algorithms when they have been completed. The ability to rapidly label data is important because of the volume of data required for training machine learning algorithms. The GUI has the flexibility to allow addition of improved ML algorithms, methods for data visualization, and statistical computations. Statistical analyses via the GUI include areas of pits within a defined range of pit depths, correlations between Red-Green-Blue (RGB) or grayscale intensity and relative surface height, covariances between values associated with features, and feature histograms. The development of supervised machine learning algorithms, however, has been hindered by a lack of training data. The machine learning algorithms for crack identification are being refined but require improvements to the true positive rate for crack detection. This shortcoming is an artifact of the limited training data currently available, perhaps more so than the structure of the neural networks. At present, the best results are had from a consensus over an ensemble of randomly generated Deep Neural Network (DNN) or Convolutional Neural Network (CNN) algorithms. Although the consensus accuracy method has yielded optimum true positive and true negative rates in excess of 80%, additional validation testing is necessary. In addition to the suite of LCM data that was initially used, and which represents the majority of the work presented in this report, WAMS image data was also reviewed at a preliminary level. The review included a comparison between image resolution and dynamic range for each method. WAMS (ZON file) image data was found to have a pixel pitch of 3.69μm compared to 1 μm for the LCM (vk4 file) data, which implies a lower resolution for the WAMS images. Conversely, the ratio of dynamic range of the WAMS data to the LCM data was approximately 41:20 for height data, suggesting that information from WAMS should more accurately determine the depth of pits. At present, the significance of the greater dynamic range of the WAMS data relative to the LCM data has not yet been evaluated.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Examining Constituent Redistribution in U-19Pu-10Zr Fuel as it Evolves with Local Burnup

While constituent redistribution is a known irradiation behavior in U-Pu-Zr fuel, new data have shown it is more complex than our current understanding and predictive capabilities. The size and composition of redistributed rings evolve as a function of pin composition, burnup, geometry, and irradiation temperature. In this work, we extract microstructural information from optical microscopy conducted on U-19Pu-10Zr pins (irradiated between 1.9 at. % and 11.6 at. % peak burnup). Both manual image analysis techniques and machine learning-assisted segmentation are used to quantify the thicknesses of the cladding, fuel-cladding interaction layers, and rings of fuel constituent redistribution in addition to pore distribution. These microstructural features and individual redistributed regions affect local thermomechanical properties, and identifying the relationship between burnup and constituent redistribution will improve accurate prediction of advanced reactor fuel performance.

Constituent Redistribution↗

Comparison of Rotor Wake Identification and Characterization Methods for the Analysis of Wake Dynamics and Evolution

Optimal wind power plant design requires understanding of wind turbine wake physics and validation of engineering wake models under wake-controlled operating conditions. In this work, we have developed and investigated several different wake identification and characterization methods for analyzing wake evolution and dynamics. The accuracy and robustness of these methods, based on Gaussian function fitting and adaptive contour identification, have been assessed by application to a large-eddy simulation data set. A new contour-based method based on downstream momentum deficit has been considered. Uncertainties arising from wake-identification errors result in characterizations of the wake expansion, recovery, and meandering motion that differ by 19% of the rotor area, 4% of the freestream, and 15% rotor diameter, respectively.

17 WIND ENERGY↗

Bayesian calibration and uncertainty quantification of a rate-dependent cohesive zone model for polymer interfaces

In this work we present a rate-dependent cohesive zone model for the fracture of polymeric interfaces and performs a Bayesian calibration, an uncertainty quantification, and a sensitivity analysis for the model. The proposed cohesive zone model accounts for both reversible elastic and irreversible rate-dependent separation sliding deformation at the interface. The viscous dissipation due to the irreversible opening at the interface is modeled using elastic-viscoplastic kinematics that incorporates the effects of strain rate. Inverse calibration of parameters for such complex models through trial and error is challenging due to the large number of parameters of the model. Moreover, the calibrated parameter values are often non-unique and uncertain when the available experimental data is limited. To tackle this challenge, we employ a Bayesian calibration approach to identify parameters from experimental data, the resulting parameters significantly enhance the accuracy of the model. To quantify the uncertainty associated with the inverse parameter estimation, a modular Bayesian approach is employed to calibrate the unknown model parameters, accounting for the parameter uncertainty of the cohesive zone model. The advantages of the Bayesian calibration over a deterministic parameter fit are demonstrated. Further, to quantify the model uncertainties, such as incorrect assumptions or missing physics, a discrepancy function is introduced, which significantly improves the model’s prediction. Finally, the total uncertainty of the model is quantified in a predictive setting. A sensitivity analysis is performed to assess how changes in the input variables of the model affect the peak load, facilitating the identification of a concise set of highly influential parameters. The present approach can be used for calibration and uncertainty quantification for other complex computational mechanics models. It should also facilitate the designing of interface materials under uncertainty.

42 ENGINEERING↗

Dark energy survey year 3 results: weak lensing shape catalogue

ABSTRACT We present and characterize the galaxy shape catalogue from the first 3 yr of Dark Energy Survey (DES) observations, over an effective area of 4143 deg2 of the southern sky. We describe our data analysis process and our self-calibrating shear measurement pipeline metacalibration, which builds and improves upon the pipeline used in the DES Year 1 analysis in several aspects. The DES Year 3 weak-lensing shape catalogue consists of 100 204 026 galaxies, measured in the riz bands, resulting in a weighted source number density of neff = 5.59 gal arcmin−2 and corresponding shape noise σe = 0.261. We perform a battery of internal null tests on the catalogue, including tests on systematics related to the point spread function (PSF) modelling, spurious catalogue B-mode signals, catalogue contamination, and galaxy properties.

79 ASTRONOMY AND ASTROPHYSICS↗

Enabling Modular Autonomous Feedback‐Loops in Materials Science through Hierarchical Experimental Laboratory Automation and Orchestration

Abstract Materials acceleration platforms (MAPs) operate on the paradigm of integrating combinatorial synthesis, high‐throughput characterization, automatic analysis, and machine learning. Within a MAP, one or multiple autonomous feedback loops may aim to optimize materials for certain functional properties or to generate new insights. The scope of a given experiment campaign is defined by the range of experiment and analysis actions that are integrated into the experiment framework. Herein, the authors present a method for integrating many actions within a hierarchical experimental laboratory automation and orchestration (HELAO) framework. They demonstrate the capability of orchestrating distributed research instruments that can incorporate data from experiments, simulations, and databases. HELAO interfaces laboratory hardware and software distributed across several computers and operating systems for executing experiments, data analysis, provenance tracking, and autonomous planning. Parallelization is an effective approach for accelerating knowledge generation provided that multiple instruments can be effectively coordinated, which the authors demonstrate with parallel electrochemistry experiments orchestrated by HELAO. Efficient implementation of autonomous research strategies requires device sharing, asynchronous multithreading, and full integration of data management in experimental orchestration, which to the best of the authors’ knowledge, is demonstrated for the first time herein.

36 MATERIALS SCIENCE↗

Measurements of inclusive and differential cross sections for the Higgs boson production and decay to four-leptons in proton-proton collisions at $ \sqrt{s} $ = 13 TeV

Measurements of the inclusive and differential fiducial cross sections for the Higgs boson production in the H → ZZ → 4ℓ (ℓ = e, μ) decay channel are presented. The results are obtained from the analysis of proton-proton collision data recorded by the CMS experiment at the CERN LHC at a center-of-mass energy of 13 TeV, corresponding to an integrated luminosity of 138 fb$^{−1}$. The measured inclusive fiducial cross section is 2.73 ± 0.26 fb, in agreement with the standard model expectation of 2.86 ± 0.1 fb. Differential cross sections are measured as a function of several kinematic observables sensitive to the Higgs boson production and decay to four leptons. A set of double-differential measurements is also performed, yielding a comprehensive characterization of the four leptons final state. Constraints on the Higgs boson trilinear coupling and on the bottom and charm quark coupling modifiers are derived from its transverse momentum distribution. All results are consistent with theoretical predictions from the standard model.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

PythonFOAM: In-situ data analyses with OpenFOAM and Python

Here, we outline the development of a general-purpose Python-based data analysis tool for OpenFOAM. Our implementation relies on the construction of OpenFOAM applications that have bindings to data analysis libraries in Python. Double precision data in OpenFOAM is cast to a NumPy array using the NumPy C-API and Python modules may then be used for arbitrary data analysis and manipulation on flow-field information. We highlight how the proposed wrapper may be used for an in-situ online singular value decomposition (SVD) implemented in Python and accessed from the OpenFOAM solver PimpleFOAM. Here, 'in-situ' refers to a programming paradigm that allows for a concurrent computation of the data analysis on the same computational resources utilized for the partial differential equation solver. In addition, to demonstrate parallel deployments, we deploy a distributed SVD, which collects snapshot data across the ranks of a distributed simulation to compute the global left singular vectors. Crucially, both OpenFOAM and Python share the same message passing interface (MPI) communicator for this deployment which allows Python objects and functions to exchange NumPy arrays across ranks. Subsequently, we provide scaling assessments of this distributed SVD on multiple nodes of Intel Broadwell and KNL architectures for canonical test cases such as the large eddy simulations of a backward facing step and a channel flow at friction Reynolds number of 395. Finally, we demonstrate the deployment of a deep neural network for compressing the flow-field information using an autoencoder to demonstrate an ability to use state-of-the-art machine learning tools in the Python ecosystem.

97 MATHEMATICS AND COMPUTING↗

Impedance-Based Detection of NO 2 Using Ni-MOF-74: Influence of Competitive Gas Adsorption

Chemically robust, low-power sensors are needed for the direct electrical detection of toxic gases. Metal–organic frameworks (MOFs) offer exceptional chemical and structural tunability to meet this challenge, though further understanding is needed regarding how coadsorbed gases influence or interfere with the electrical response. To probe the influence of competitive gases on trace NO 2 detection in a simulated flue gas stream, a combined structure–property study integrating synchrotron powder diffraction and pair distribution function analyses was undertaken, to elucidate how structural changes associated with gas binding inside Ni-MOF-74 pores correlate with the electrical response from Ni-MOF-74-based sensors. Data were evaluated for 16 gas combinations of N 2 , NO 2 , SO 2 , CO 2 , and H 2 O at 50 °C. Fourier difference maps from a rigid-body Rietveld analysis showed that additional electron density localized around the Ni-MOF-74 lattice correlated with large decreases in Ni-MOF-74 film resistance of up to a factor of 6 × 10 3 , observed only when NO 2 was present. These changes in resistance were significantly amplified by the presence of competing gases, except for CO 2 . Without NO 2 , H 2 O rapidly (<120 s) produced small (1–3×) decreases in resistance, though this effect could be differentiated from the slower adsorption of NO 2 by the evaluation of the MOF’s capacitance. Furthermore, samples exposed to H 2 O displayed a significant shift in lattice parameters toward a larger lattice and more diffuse charge density in the MOF pore. Evaluating the Ni-MOF-74 impedance in real time, NO 2 adsorption was associated with two electrically distinct processes, the faster of which was inhibited by competitive adsorption of CO 2 . Together, this work points to the unique interaction of NO 2 and other specific gases (e.g., H 2 O, SO 2 ) with the MOF’s surface, leading to orders of magnitude decrease in MOF resistance and enhanced NO2 detection. Finally, understanding and leveraging these coadsorbed gases will further improve the gas detection properties of MOF materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Machine learned interatomic potential for dispersion strengthened plasma facing components

Tungsten (W) is a material of choice for the divertor material due to its high melting temperature, thermal conductivity, and sputtering threshold. However, W has a very high brittle-to-ductile transition temperature, and at fusion reactor temperatures (≥1000 K), it may undergo recrystallization and grain growth. Dispersion-strengthening W with zirconium carbide (ZrC) can improve ductility and limit grain growth, but much of the effects of the dispersoids on microstructural evolution and thermomechanical properties at high temperatures are still unknown. We present a machine learned Spectral Neighbor Analysis Potential for W–ZrC that can now be used to study these materials. In order to construct a potential suitable for large-scale atomistic simulations at fusion reactor temperatures, it is necessary to train on ab initio data generated for a diverse set of structures, chemical environments, and temperatures. Further accuracy and stability tests of the potential were achieved using objective functions for both material properties and high temperature stability. Additionally, validation of lattice parameters, surface energies, bulk moduli, and thermal expansion is confirmed on the optimized potential. Tensile tests of W/ZrC bicrystals show that although the W(110)–ZrC(111) C-terminated bicrystal has the highest ultimate tensile strength (UTS) at room temperature, observed strength decreases with increasing temperature. At 2500 K, the terminating C layer diffuses into the W, resulting in a weaker W–Zr interface. Meanwhile, the W(110)–ZrC(111) Zr-terminated bicrystal has the highest UTS at 2500 K.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Study of e + e − → π + π − π 0 at s from 2.00 to 3.08 GeV at BESIII

With the data samples taken at center-of-mass energies from 2.00 to 3.08 GeV with the BESIII detector at the BEPCII collider, a partial wave analysis on the e + e − → π + π − π 0 process is performed. The Born cross sections for e + e − → π + π − π 0 and its intermediate processes e + e − → ρ π and ρ ( 1450 ) π are measured as functions of s . The results for e + e − → π + π − π 0 are consistent with previous results measured with the initial state radiation method within one standard deviation, and improve the uncertainty by a factor of ten. By fitting the line shapes of the Born cross sections for the e + e − → ρ π and e + e − → ρ ( 1450 ) π , a structure with mass M = 2119 ± 11 ± 15 MeV / c 2 and width Γ = 69 ± 30 ± 5 MeV is observed with a significance of 5.9 σ , where the first uncertainties are statistical and the second ones are systematic. This structure can be interpreted as an excited ω state. Published by the American Physical Society 2024

Astronomy & Astrophysics↗

VISILIENCE: An Interactive Visualization Framework for Resilience Analysis using Control-Flow Graph

Soft errors have become one of the major concerns for the error resilience of HPC applications, as those errors can cause HPC applications to generate serious outcomes such as Silent Data Corruptions (SDCs). A large body of approaches has been proposed to analyze the resilience of HPC applications. However, existing studies rarely address the challenges of the analysis result perception. Specifically, resilience analysis techniques often produce a massive volume of unstructured data, making it difficult for programmers to conduct the resilience analysis due to non-intuitive raw data. Furthermore, different analysis models produce diverse results with multiple levels of details, which may create hurdles to compare and explore the resilience of HPC program execution. To this end, we present VISILIENCE, an interactive VISual resILIENCE analysis framework to allow programmers to facilitate the resilience analysis of HPC applications. In particular, VISILIENCE leverages an effective visualization approach Control Flow Graph (CFG) to present a function execution. In addition, three widely-used models for resilience analysis (i.e., Y-Branch, IPAS, and TRIDENT) are seamlessly embedded into the framework for resilience analysis and result comparison. Multiple case studies have been conducted to demonstrate the effectiveness of our proposed framework VISILIENCE.

Jiang, Hailong↗

Using Simulation to Interpret LArS's Results

Scintillation light analysis in liquid argon based neutrino detectors is restrained in capability due to uncertainty in fundamental constants critical to the analysis process. One such property is the Rayleigh scattering length of liquid argon. In the fall of 2023, the TallBo cryostat, located in the Noble Liquid Testing Facility (NLTF) at Fermilab, was used to study the scattering length of liquid argon in the Liquid Argon Scattering (LArS) experiment. Due to systematic errors unknown during measurement analysis, LArS s measurements were quite uncertain. By simulating the LArS experiment, we found that the downturn in detector count rate as a function of liquid argon height at low heights was caused by a misplaced silicon photo multiplier (SiPM). With concentrated effort, we may be able to successfully correct this effect by understanding the relationship between the specified Rayleigh scattering length and the measured attenuation length. With this information and further progression in analysis, we may be able to extract corrected measurements from the LArS data and attain a tangible experimental measurement of the scattering length of liquid argon.

Breaux, Auto D. [Tulane U.]↗

WM26 Paper Multi-Robot Collaboration for Hazardous Environments

Hazardous nuclear and industrial facilities are rarely designed for robots. Work in these domains demand precise manipulation and robust mobility in cluttered, constrained spaces where off-the-shelf platforms struggle and “one-size-fits-all” machines become costly and complex. Idaho National Laboratory (INL) is developing an autonomous, multi-robot inspection system that coordinates task-specific platforms rather than relying on a single omni-tool robot. An electric truck serves as a power and compute hub for a custom manipulator co-developed with Florida International University (FIU), a commercial mini crawler, a pan–tilt–zoom camera, and a Nexxis Argus LiDAR mapping system. Working in concert, these robots generate spatial, radiation, and temperature maps of the pit environments at the Hanford Waste Tank Farms. These systems will capture visual records and environmental telemetry to allow for analysis post inspection. The system architecture uses Robot Operating System 2 (ROS 2) for publish/subscribe integration, NVIDIA Isaac Sim and Unity for simulation and visualization, and algorithms such as NVBlox to fuse data into unified 3D overlays. This robot-agnostic approach reduces operator burden by enabling autonomy across heterogeneous platforms and lets each robot be used where it is strongest. Having autonomous functions means operators don’t have to fully control multiple different components. The ease of use could allow for more widespread adoption of advanced robotics at waste management sites that see continued use. By coordinating simpler, purpose-built mechanisms, the approach lowers design and manufacturing complexity, reduces capital risk in contaminated settings, and improves controllability for complex inspection and manipulation tasks. We present the architecture, early results, and lessons learned from building and deploying this coordinated multi-robot system, with the goal of accelerating safe, cost-effective adoption of advanced robotics at waste-management sites.

42 - ENGINEERING↗

Uncertainty about the Uncertainty [Slides]

A golden standard in science is to repeat an experiment a statistically significant number of times, recording data using the same set of detectors and the same data analysis methodology. In such case experimental error includes both the range of true values generated by repetitions of the experiment, and measurement uncertainty caused by the detector. They are independent. It is a huge and too frequently used simplification, to assume that one can measure multiple repetitions of an identical experiment, resulting in identical true experimental value. Repetitions, as similar is it is experimentally achievable result in a range of the true values rather than in a single value. When modern, very sensitive and well calibrated measurement systems are used, this range is not negligible, and sometimes dominates, in comparison to the measurement uncertainty. Range of true values depends on the physics of the experiment, while measurement uncertainty depends on the measurement method (properties of the detector not of the experiment). When data from one–of–a kind experiment are analyzed, only the measurement uncertainty is reported. It gives no information about the range, in which the true values of experiment would spread if the experiment was repeated. A frequently used approximation, that if a physical quantity is measured as a function of time, only measurement of this quantity, produces uncertainty is also in some real experiments fare to strong. Example: In reaction history time measurement uncertainty propagated to alpha dominated under certain conditions over the flux measurement uncertainty propagated to alpha. Reliability of a data point is in general independent from its measurement uncertainty. However, in practice reliable measurement methods frequently have high measurement uncertainty, while low reliability methods are applied to limit measurement uncertainty. Comparison of reliable data with high measurement uncertainty to not so reliable data measured with low uncertainty is discussed – in different scenarios different data analysis methods are applicable.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗