Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “sampling 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 739 records · Page 41

Waveform resampling with LMN method

In this article, resampling is a common technique applied in digital signal processing. Based on the Fast Fourier Transformation (FFT), we apply an optimization called here the LMN method to achieve fast and robust re-sampling. In addition to performance comparisons with some other popular methods, we illustrate the effectiveness of this LMN method in a particle physics experiment: re-sampling of waveforms from Liquid Argon Time Projection Chambers.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Excess free volume and structural properties of inert gas condensation synthesized nanoparticles based CuZr nanoglasses

Abstract Nanoglass (NG) as a new structure-tunable material has been investigated using both experiments and computational modeling. Experimentally, inert gas condensation (IGC) is commonly employed to prepare metallic glass (MG) nanoparticles that are consolidated using cold compression to generate an NG. In computational modeling, various methods have been used to generate NGs. However, due to the high computational cost involved, heretofore modeling investigations have not followed the experimental synthesis route. In this work, we use molecular dynamics simulations to generate an NG model by consolidating IGC-prepared Cu 64 Zr 36 nanoparticles following a workflow similar to that of experiments. The resulting structure is compared with those of NGs produced following two alternative procedures previously used: direct generation employing Voronoi tessellation and consolidation of spherical nanoparticles carved from an MG sample. We focus on the characterization of the excess free volume and the Voronoi polyhedral statistics in order to identify and quantify contrasting features of the glass-glass interfaces in the three NG samples prepared using distinct methods. Results indicate that glass-glass interfaces in IGC-based NGs are thicker and display higher structural contrast with their parent MG structure. Nanoparticle-based methods display excess free volume exceeding 4%, in agreement with experiments. IGC-prepared nanoparticles, which display Cu segregation to their surfaces, generate the highest glass-glass interface excess free volume levels and the largest relative interface volume with excess free volume higher than 3%. Voronoi polyhedral analysis indicates a sharp drop in the full icosahedral motif fraction in the glass-glass interfaces in nanoparticle-based NG as compared to their parent MG.

36 MATERIALS SCIENCE↗

In-Field Alpha Spectrometer Development FY2022 Mid-Year Report

In response to needs identified by the International Atomic Energy Agency (IAEA), Idaho National Laboratory (INL) has developed an In-Field Alpha Spectrometer (IFAS) to allow IAEA safeguards inspectors to collect samples of uranium hexafluoride (UF 6 ) at processing facilities, to perform field assessments to verify uranium enrichment. For sample collection, the IFAS method uses Single-Use Destructive Assay (SUDA) samplers, developed at Pacific Northwest National Laboratory (PNNL), which contain thin zeolite coatings that trap UF6 gas and convert it to the safer, more stable form uranyl fluoride (as a dihydrate, UO 2 F 2 ·2H 2 O). For alpha spectrometry, the IFAS instrument employs a large area silicon semiconductor transducer to detect and record alpha particle energy-deposition events. Over the course of this project INL, PNNL, Oak Ridge National Laboratory worked to optimize the SUDA sample design for alpha spectrometry, to optimize processes for manufacturing, loading, and shipping SUDA samples, and to allow the measurement of loaded SUDA samples. In fiscal year (FY) 2022 INL has bene investigating the reproducibility and measurement precision of IFAS measurements and the associated determination of uranium enrichment.

46 - INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AN↗

Bilinear dynamic mode decomposition for quantum control

Abstract Data-driven methods for establishing quantum optimal control (QOC) using time-dependent control pulses tailored to specific quantum dynamical systems and desired control objectives are critical for many emerging quantum technologies. We develop a data-driven regression procedure, bilinear dynamic mode decomposition (biDMD), that leverages time-series measurements to establish quantum system identification for QOC. The biDMD optimization framework is a physics-informed regression that makes use of the known underlying Hamiltonian structure. Further, the biDMD can be modified to model both fast and slow sampling of control signals, the latter by way of stroboscopic sampling strategies. The biDMD method provides a flexible, interpretable, and adaptive regression framework for real-time, online implementation in quantum systems. Further, the method has strong theoretical connections to Koopman theory, which approximates nonlinear dynamics with linear operators. In comparison with many machine learning paradigms minimal data is needed to construct a biDMD model, and the model is easily updated as new data is collected. We demonstrate the efficacy and performance of the approach on a number of representative quantum systems, showing that it also matches experimental results.

97 MATHEMATICS AND COMPUTING↗

Strategic Placement and Sizing of Distributed Generation for Resilience Enhancement of Distribution Grids With Microgrid Formation

The rise in frequency and severity of extreme weather events highlights the need for resilient power distribution networks. Microgrids can help improve the resilience of distribution grids by providing continuous power supply using local distribution generation (DG) when the distribution grid fails. In this paper, we propose an approach for optimal placement and sizing of DG to form multiple microgrids throughout the distribution network by restoration actions such as switching operations in case of distribution grid outages caused by extreme weather events. Considering the randomness of damaged distribution lines, the DG placement and sizing problem is formulated as a two-stage stochastic mixed-integer program, with the first stage determining the placement and size of DG, and the second stage focusing on minimizing the amount of load shedding through network restoration and microgrid formations for each scenario. Due to the large number of scenarios, the sample average approximation (SAA) method is employed to solve the problem. The results of case studies on a modified IEEE 33 bus distribution grid demonstrate the effectiveness of the proposed DG placement and sizing strategy in improving the resilience of distribution grids by allowing the formation of multiple microgrids. In addition, the robustness and accuracy of the SAA method are validated through various case studies.

Distributed generation planning↗

Improved Infrared Optical Constants from Pressed Pellets: II. Ellipsometric n and k Values for Ammonium Sulfate with Variability Analysis

Infrared spectral imaging is facilitated via the comparison of spectra recorded in situ to a databank of actual or synthetic infrared reflectance spectra. It has recently been shown that reference spectra corresponding to the many different morphological forms of the same chemical can be generated synthetically using the imaginary, k, and real, n, components of the complex refractive index, n~ = n + ik. One method to obtain the n and k vectors is infrared ellipsometry, which measures the changes in amplitude, tan ?, and phase, ?, of polarized light reflected from the sample both as a function of wavenumber and angle of incidence. The method requires specularly reflected light so best results are usually obtained with polished planar samples of large surface area. Due to the difficulties of obtaining such samples, however, we investigate the possibility of pressing powders of neat materials and obtaining the corresponding optical constants from the pellets. In this paper, variability in the sample pellet and preparation method is investigated, as is variability in the fitting procedure for the derived optical constants. The n/k vectors are derived from the measured ellipsometric parameters, tan ? and ?, as they are fit by an oscillator model which yield n(?) and k(?) vectors as a function of wavenumber, ?. Construction of the oscillator model is not automatic and depends on significant input from the analyst as well as the sample’s physical characteristics. For the pellet pressing, the experimental variability was found to be minimized for size-selected powdered samples as gauged by the minimal variance in ? and ? for three different pellets; similarly, the analytical precision for multiple measurements of the same pellet was also quite good, suggesting that a pressed pellet is a viable sample preparation method. The greatest variability came in the analytic fitting procedure with differences up to 20% arising from four different fits to the same data set. The final ellipsometric n/k data are compared to literature values obtained from crystalline (NH4)2SO4 samples as well as single-angle reflectance measurements that also used pressed pellets. Comparison with the previous literature values shows generally good agreement although higher k-values are observed for the independent sets of data derived from pressed pellets. These data are suggested as an improved set of optical constants for (NH4)2SO4.

optical constants, ellipsometry, ammonium sulfate,↗

Computational catalyst discovery: Active classification through myopic multiscale sampling

We report the recent boom in computational chemistry has enabled several projects aimed at discovering useful materials or catalysts. We acknowledge and address two recurring issues in the field of computational catalyst discovery. First, calculating macro-scale catalyst properties is not straightforward when using ensembles of atomic-scale calculations [e.g., density functional theory (DFT)]. We attempt to address this issue by creating a multi-scale model that estimates bulk catalyst activity using adsorption energy predictions from both DFT and machine learning models. The second issue is that many catalyst discovery efforts seek to optimize catalyst properties, but optimization is an inherently exploitative objective that is in tension with the explorative nature of early-stage discovery projects. In other words, why invest so much time finding a “best” catalyst when it is likely to fail for some other, unforeseen problem? We address this issue by relaxing the catalyst discovery goal into a classification problem: “What is the set of catalysts that is worth testing experimentally?” Here, we present a catalyst discovery method called myopic multiscale sampling, which combines multiscale modeling with automated selection of DFT calculations. It is an active classification strategy that seeks to classify catalysts as “worth investigating” or “not worth investigating” experimentally. Our results show an ~7–16 times speedup in catalyst classification relative to random sampling. These results were based on offline simulations of our algorithm on two different datasets: a larger, synthesized dataset and a smaller, real dataset.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Challenges in the Use of AI-Driven Non-Destructive Spectroscopic Tools for Rapid Food Analysis

Routine, remote, and process analysis for foodstuffs is gaining attention and can provide more confidence for the food supply chain. A new generation of rapid methods is emerging both in the literature and in industry based on spectroscopy coupled with AI-driven modelling methods. Current published studies using these advanced methods are plagued by weaknesses, including sample size, abuse of advanced modelling techniques, and the process of validation for both the acquisition method and modelling. This paper aims to give a comprehensive overview of the analytical challenges faced in research and industrial settings where screening analysis is performed while providing practical solutions in the form of guidelines for a range of scenarios. After extended literature analysis, we conclude that there is no easy way to enhance the accuracy of the methods by using state-of-the-art modelling methods and the key remains that capturing good quality raw data from authentic samples in sufficient volume is very important along with robust validation. A comprehensive methodology involving suitable analytical techniques and interpretive modelling methods needs to be considered under a tailored experimental design whenever conducting rapid food analysis.

59 BASIC BIOLOGICAL SCIENCES↗

Accelerating Markov Chain Monte Carlo sampling with diffusion models

Global fits of physics models require efficient methods for exploring high-dimensional and/or multimodal posterior functions. We introduce a novel method for accelerating Markov Chain Monte Carlo (MCMC) sampling by pairing a Metropolis-Hastings algorithm with a diffusion model that can draw global samples with the aim of approximating the posterior. We briefly review diffusion models in the context of image synthesis before providing a streamlined diffusion model tailored towards low-dimensional data arrays. We then present our adapted Metropolis-Hastings algorithm which combines local proposals with global proposals taken from a diffusion model that is regularly trained on the samples produced during the MCMC run. Our approach leads to a significant reduction in the number of likelihood evaluations required to obtain an accurate representation of the Bayesian posterior across several analytic functions, as well as for a physical example based on a global fit of parton distribution functions. Our method is extensible to other MCMC techniques, and we briefly compare our method to similar approaches based on normalising flows. A code implementation can be found at https://github.com/NickHunt-Smith/MCMC-diffusion.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Divide and conquer: using RhizoVision Explorer to aggregate data from multiple root scans using image concatenation and statistical methods

Roots are important in agricultural and natural systems for determining plant productivity and soil carbon inputs. Sometimes, the amount of roots in a sample is too much to fit into a single scanned image, so the sample is divided among several scans, and there is no standard method to aggregate the data. Here, we describe and validate two methods for standardizing measurements across multiple scans: image concatenation and statistical aggregation. We developed a Python script that identifies which images belong to the same sample and returns a single, larger concatenated image. These concatenated images and the original images were processed with RhizoVision Explorer, a free and open-source software. An R script was developed, which identifies rows of data belonging to the same sample and applies correct statistical methods to return a single data row for each sample. These two methods were compared using example images from switchgrass, poplar, and various tree and ericaceous shrub species from a northern peatland and the Arctic. Most root measurements were nearly identical between the two methods except median diameter, which cannot be accurately computed by statistical aggregation. We believe the availability of these methods will be useful to the root biology community.

59 BASIC BIOLOGICAL SCIENCES↗

Characterization of Acoustic Emissions From Analogue Rocks Using Sparse Regression‐DMDc

Abstract Moisture loss in rock is known to generate acoustic emissions (AE). Phenomena that result in AE during drying are related to the movement of fluids through the pores and induced‐cracks that arise from differential mineral shrinkage, especially in clay‐bearing rock. AE from the movement of fluids occurs from the reconfiguration of fluid interfaces during drying, while AE from mineral shrinkage involves the debonding within or between minerals. Here, analogue rock samples were used to examine the differences in the AE signatures when one or both AE source‐types are present. An unsupervised sparse regression model, Dynamic Mode Decomposition with control, that extends Dynamic Mode Decomposition is used to characterize the AE signals recorded during the drying of porous analogue rock samples fabricated with ordinary Portland cement, with and without clay. This method can effectively and accurately reconstruct acoustic signals emitted from samples that only experience moisture loss without cracking. However, the method struggles to reconstruct signals from samples with intricate crack networks that formed during drying because AE generating mechanisms can emit contemporaneously, and the resulting waves propagate through drying‐induced cracks that can lead to multiple internal reflections. Thus, the differential reconstruction accuracy of time series generated by different underlying physical processes provides a robust filter for reducing large data catalogs. In general, both dynamics and sparse initiating events are learned directly from data and this method exposes a data hierarchy based on the complexity of the intrinsic dynamics.

58 GEOSCIENCES↗

BayFlux: A Bayesian method to quantify metabolic Fluxes and their uncertainty at the genome scale

Metabolic fluxes, the number of metabolites traversing each biochemical reaction in a cell per unit time, are crucial for assessing and understanding cell function. 13 C Metabolic Flux Analysis ( 13 C MFA) is considered to be the gold standard for measuring metabolic fluxes. 13 C MFA typically works by leveraging extracellular exchange fluxes as well as data from 13 C labeling experiments to calculate the flux profile which best fit the data for a small, central carbon, metabolic model. However, the nonlinear nature of the 13 C MFA fitting procedure means that several flux profiles fit the experimental data within the experimental error, and traditional optimization methods offer only a partial or skewed picture, especially in “non-gaussian” situations where multiple very distinct flux regions fit the data equally well. Here, we present a method for flux space sampling through Bayesian inference (BayFlux), that identifies the full distribution of fluxes compatible with experimental data for a comprehensive genome-scale model. This Bayesian approach allows us to accurately quantify uncertainty in calculated fluxes. We also find that, surprisingly, the genome-scale model of metabolism produces narrower flux distributions (reduced uncertainty) than the small core metabolic models traditionally used in 13 C MFA. The different results for some reactions when using genome-scale models vs core metabolic models advise caution in assuming strong inferences from 13 C MFA since the results may depend significantly on the completeness of the model used. Based on BayFlux, we developed and evaluated novel methods (P- 13 C MOMA and P- 13 C ROOM) to predict the biological results of a gene knockout, that improve on the traditional MOMA and ROOM methods by quantifying prediction uncertainty.

59 BASIC BIOLOGICAL SCIENCES↗

CADIS and FW-CADIS Variance Reduction in Gamma Transport for Predicting Prompt Forensics Signatures

The goal of prompt nuclear forensics is to determine the characteristics of a nuclear detonation based on the signatures available almost immediately after the explosion. An important characteristic is the reaction time history (RTH), a measure of the device’s rate of neutron multiplication. The RTH can be estimated by observation of the gamma radiation emitted from the detonation, which can be detected directly or observed indirectly as Teller light. Gamma transport simulations used to predict these radiation fields are often modeled stochastically using the Monte Carlo N-Particle (MCNP) code, which can be a computationally demanding task due to the number of particle histories needed to achieve statistical convergence. In an attempt to improve the efficiency of these calculations, we evaluate two variance reduction techniques: Consistent Adjoint-Driven Importance Sampling (CADIS) and Forward-Weighted Consistent Adjoint-Driven Importance Sampling (FW-CADIS). These methods use a deterministically calculated adjoint flux to create weight windows and source biasing that guide MCNP sampling. We study the utility of CADIS and FW-CADIS for their use in MCNP gamma transport for nuclear forensics prediction simulations. Furthermore, the results demonstrate that both CADIS and FW-CADIS improve the accuracy for forensics-focused simulations, with CADIS being most beneficial in direct detection and FW-CADIS being ideal for computing a global Teller light source.

CADIS↗

Multiclass Classification Using Bayesian Multivariate Adaptive Regression Splines

We present a new Bayesian model for the problem of multiclass classification. In this model, the probabilities of class membership of a given observation are determined by the mean of a latent Gaussian distribution. The mean functions of this latent distribution consist of combinations of highly flexible basis functions of the inputs: multivariate adaptive regression splines (MARS), first developed for multiple regression. We use reversible jump Markov chain Monte Carlo to make inference on the classification model, including the number of basis functions. We compare the probabilistic classification performance of our proposed approach to existing methods on simulated and benchmark data, and compare uncertainty estimates on simulated data. Our proposed method compares favorably with existing Bayesian and frequentist multiclass classification methods in out-of-sample probabilistic classification, and uncertainty estimation of these probabilistic classifications. We examine the fit of the proposed method to a data set of hurricane storm surge levels near Delaware Bay, US, and conclude that sea level rise is a key contributor to damage delivered by storm surge.

97 MATHEMATICS AND COMPUTING↗

Coarse-Grained Density Functional Theory Predictions via Deep Kernel Learning

Scalable electronic predictions are critical for soft materials design. Recently, the Electronic Coarse-Graining (ECG) method was introduced to renormalize all-atom quantum chemical (QC) predictions to coarse-grained (CG) resolutions using deep neural networks (DNNs). While DNNs can learn complex representations that prove challenging for kernel-based methods, they are susceptible to overfitting and the overconfidence of uncertainty estimations. Here, we develop ECG within a GPU-accelerated Deep Kernel Learning (DKL) framework to enable CG QC predictions using range-separated hybrid density functional theory (DFT), obtaining a 107 speedup relative to naive all-atom QC. By treating the predicted electronic properties as random Gaussian Processes, DKL incorporates CG mapping degeneracy by learning the distribution of electronic energies as a function of CG configuration. DKL-ECG accurately reproduces molecular orbital energies from range-separated DFT while facilitating efficient training via active learning using the uncertainties provided by DKL. Further, we show that while active learning algorithms enable efficient sampling of a more diverse configurational space relative to random sampling, all explored query methods exhibit comparable performance for the examined system. We attribute this result to the significant overlap of the feature space and output property distributions across multiple temperatures.

97 MATHEMATICS AND COMPUTING↗

A Rapid Sequential Separation Method for Determination of Actinides and TENORM in Fracking Wastes - 20223

With recent advances in unconventional drilling technology in the US and around the world, the risks of environmental contamination and exposure due to technologically enhanced naturally occurring radioactive materials (TENORM) from fracking wastes have also surfaced. Although the Permian Basin (southeastern New Mexico, western Texas) has long been a modest producer of oil, the advent of fracking technology a decade ago turned it into hot property for producers. A 2018 assessment of undiscovered, technically recoverable continuous oil and gas resources by the U.S. Geological Survey estimates an average 46.3 billion barrels of oil and 281 trillion cubic feet of gas. This is likely to increase the use of unconventional drilling and TENORM generation. TENORM concentrations in oil and gas exploration and production waste can be several hundred to several thousand times higher than background TENORM concentrations. Both the thorium and uranium decay chains contribute to airborne radionuclides arising from the radon gas escaping the ground and subsequently decaying as airborne particulates. Additionally, radon decay products {sup 210}Pb and {sup 210}Po can build up in scale on the internal surfaces of oil and gas handling pipes and in sludge in refineries, becoming potential inhalation and ingestion hazards for workers. Furthermore, southeastern NM is also home to world's only licensed and operating transuranic nuclear waste repository, the DoE's Waste Isolation Pilot Plant (WIPP). Plutonium isotopes ({sup 239+240}Pu) and {sup 241}Am, are expected to account for more than 99% of the total radioactivity scheduled for disposal in the WIPP repository. Thus, an improved understanding of the environmental fate and transport of TENORM, liberated by unconventional drilling, is essential to assess how best to protect individuals and the environment. In this context, accurate measurement of TENORM and actinides in environmental and biological samples is essential. In this presentation, a new sequential method for the separation and pre-concentration of actinides (Pu, Am, Np) and TENORM (Po, U and Th) derived from oil and gas exploration is proposed. The TEVA method involves a rapid co-precipitation step to remove matrix interferences from the samples, followed by plutonium oxidation state adjustment to Pu (IV) and an incubation period of ∼ 1 hour at 50-60 deg. C to allow the resultant Po (II) to oxidize into Po(IV). The polonium, neptunium, thorium and plutonium are then separated on a TEVA column, while americium and uranium are separated on a TRU column. After separation, the alpha counting source was prepared by micro-precipitation with copper sulfide (CuS) for polonium and neodymium fluoride (NdF{sub 3}) for actinides. The efficiency and reliability of the procedures were tested by analyzing filter, drinking water and frack sand samples. The method is simple, robust and can be performed quickly with excellent removal of interferences, high chemical recovery and very good alpha peak resolution. (authors)

07 ISOTOPE AND RADIATION SOURCES↗

A multi-scale temperature-based strategy to map hydrologic exchange flows in highly dynamic systems

Mapping and quantifying hydrologic exchange flows (HEFs) is critical to environmental monitoring and remediation at contaminated sites; however, these objectives are challenging in highly dynamic systems, e.g., along dam-regulated rivers, where HEFs vary rapidly. Direct seepage measurements are labor-intensive and difficult to automate, whereas indirect (e.g., thermal) and remote sensing methods have potential to allow continuous monitoring with limited field effort. We present a preliminary assessment of a multi-scale temperature-based strategy for monitoring HEFs along the Hanford Reach of the Columbia River, in eastern WA, United States. Five thermal methods were assessed. First, a vertical temperature profile (VTP) was installed into the streambed. The VTP data were analyzed using a data assimilation algorithm designed for automated real-time estimation in dynamic systems. Second, a thermal infrared (TIR) camera was used in roving surveys to identify seeps. Third, a TIR camera was stationed at the VTP site to collect images at 1-h intervals. Together, the two TIR datasets provided a basis to assess the potential for drone-based TIR. Fourth, temperature was measured at the sediment/water interface to assess fiber-optic distributed temperature sensing. Fifth, imagery from the ECOSTRESS satellite mission was acquired to assess the potential of spaceborne thermal monitoring. Based on our preliminary assessment, VTP, TIR, and bed temperature measurements provide complementary spatial coverage, temporal sampling, and resolution; these methods have potential for long-term, automated monitoring of HEFs. The publicly available spaceborne imagery, however, proved inadequate because of insufficient spatial resolution and data gaps resulting from cloud cover and revisit frequency.

54 ENVIRONMENTAL SCIENCES↗

Weak Gravitational Lensing around Low Surface Brightness Galaxies in the DES Year 3 Data

We present galaxy-galaxy lensing measurements using a sample of low surface brightness galaxies (LSBGs) drawn from the Dark Energy Survey Year 3 (Y3) data as lenses. LSBGs are diffuse galaxies with a surface brightness dimmer than the ambient night sky. These dark-matter-dominated objects are intriguing due to potentially unusual formation channels that lead to their diffuse stellar component. Given the faintness of LSBGs, using standard observational techniques to characterize their total masses proves challenging. Weak gravitational lensing, which is less sensitive to the stellar component of galaxies, could be a promising avenue to estimate the masses of LSBGs. Our LSBG sample consists of 23,790 galaxies separated into red and blue color types at g - i ≥ 0.60 and g - i < 0.60 , respectively. Combined with the DES Y3 shear catalog, we measure the tangential shear around these LSBGs and find signal-to-noise ratios of 6.67 for the red sample, 2.17 for the blue sample, and 5.30 for the full sample. We use the clustering redshifts method to obtain redshift distributions for the red and blue LSBG samples. Assuming all red LSBGs are satellites, we fit a simple model to the measurements and estimate the host halo mass of these LSBGs to be log(M host /M ⊙ ) = $12.98^{+0.10}_{-0.11}$. We place a 95% upper bound on the subhalo mass at log(M sub /M ⊙ ) < 11.51. By contrast, we assume the blue LSBGs are centrals, and place a 95% upper bound on the halo mass at log(M host /M ⊙ ) < 11.84. We find that the stellar-to-halo mass ratio of the LSBG samples is consistent with that of the general galaxy population. This work illustrates the viability of using weak gravitational lensing to constrain the halo masses of LSBGs.

79 ASTRONOMY AND ASTROPHYSICS↗