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 1,369 records · Page 76

General approaches for shear-correcting coordinate transformations in Bragg coherent diffraction imaging. Part II

X-ray Bragg coherent diffraction imaging (BCDI) has been demonstrated as a powerful 3D microscopy approach for the investigation of sub-micrometre-scale crystalline particles. The approach is based on the measurement of a series of coherent Bragg diffraction intensity patterns that are numerically inverted to retrieve an image of the spatial distribution of the relative phase and amplitude of the Bragg structure factor of the diffracting sample. This 3D information, which is collected through an angular rotation of the sample, is necessarily obtained in a non-orthogonal frame in Fourier space that must be eventually reconciled. To deal with this, the approach currently favored by practitioners (detailed in Part I) is to perform the entire inversion in conjugate non-orthogonal real- and Fourier-space frames, and to transform the 3D sample image into an orthogonal frame as a post-processing step for result analysis. In this article, which is a direct follow-up of Part I, two different transformation strategies are demonstrated, which enable the entire inversion procedure of the measured data set to be performed in an orthogonal frame. The new approaches described here build mathematical and numerical frameworks that apply to the cases of evenly and non-evenly sampled data along the direction of sample rotation ( i.e. the rocking curve). The value of these methods is that they rely on the experimental geometry, and they incorporate significantly more information about that geometry into the design of the phase-retrieval Fourier transformation than the strategy presented in Part I. Finally, two important outcomes are (1) that the resulting sample image is correctly interpreted in a shear-free frame and (2) physically realistic constraints of BCDI phase retrieval that are difficult to implement with current methods are easily incorporated. Computing scripts are also given to aid readers in the implementation of the proposed formalisms.

36 MATERIALS SCIENCE↗

An iterative method to deblend AGN-Host contributions for Integral Field spectroscopic observations

ABSTRACT We present a new iterative deblending method to separate the host galaxy (HG) and their Active Galactic Nuclei (AGNs) emission with the use of Integral Field spectroscopic (IFS) data. The method decomposes the resolved HG emission from the unresolved AGN emission by modelling the two-dimensional surface brightness (SB) profile of the point-spread function (PSF) and the two-dimensional SB HG continuum simultaneously per each monochromatic slide. Our method does not require any prior information about the observed SB profile or a detailed fitting of the PSF, making it ideal for the automatic analysis of large galaxy samples. In this work, we test the quality of our method, its advantages, and its disadvantages. We test our method by using a set of IFS mock data cubes to quantify the reliability of our deblending process and further compare our method with the qdblend3d analysis tool. Furthermore, we applied our method to three data cubes selected from the MaNGA survey according to the dominance of either its HG or its AGN. We show that our deblending method is capable of disengaging the bright, non-resolved AGN emission from the HG continuum and its narrow emission lines. However, the decoupling depends on how well the IFS spatially resolves the PSF, and on the relative flux intensity of the HG-AGN. Therefore, the method is ideal for disentangling the bright-flux contribution from AGN-dominated spectra.

Ibarra-Medel, H. (ORCID:0000000297906313)↗

Distributed Optimization for Nonrigid Nano-Tomography

Resolution level and reconstruction quality in nano-computed tomography (nano-CT) are in part limited by the stability of microscopes, because the magnitude of mechanical vibrations during scanning becomes comparable to the imaging resolution, and the ability of the samples to resist radiation induced deformations during data acquisition. In such cases, there is no incentive in recovering the sample state at different time steps like in time-resolved reconstruction methods, but instead the goal is to retrieve a single reconstruction at the highest possible spatial resolution and without any imaging artifacts. Here we propose a distributed optimization solver for tomographic imaging of samples at the nanoscale. Our approach solves the tomography problem jointly with projection data alignment, nonrigid sample deformation correction, and regularization. Projection data consistency is regulated by dense optical flow estimated by Farneback's algorithm, leading to sharp sample reconstructions with less artifacts. Synthetic data tests show robustness of the method to Poisson and low-frequency background noise. We accelerated the solver on multi-GPU systems and validated the method on three nano-imaging experimental data sets.

97 MATHEMATICS AND COMPUTING↗

Context-aware learning of hierarchies of low-fidelity models for multi-fidelity uncertainty quantification

Multi-fidelity Monte Carlo methods leverage low-fidelity and surrogate models for variance reduction to make tractable uncertainty quantification even when numerically simulating the physical systems of interest with high-fidelity models is computationally expensive. This work proposes a context-aware multi-fidelity Monte Carlo method that optimally balances the costs of training low-fidelity models with the costs of Monte Carlo sampling. It generalizes the previously developed context-aware bi-fidelity Monte Carlo method to hierarchies of multiple models and to more general types of low-fidelity models. When training low-fidelity models, the proposed approach takes into account the context in which the learned low-fidelity models will be used, namely for variance reduction in Monte Carlo estimation, which allows it to find optimal trade-offs between training and sampling to minimize upper bounds of the mean-squared errors of the estimators for given computational budgets. This is in stark contrast to traditional surrogate modeling and model reduction techniques that construct low-fidelity models with the primary goal of approximating well the high-fidelity model outputs and typically ignore the context in which the learned models will be used in upstream tasks. Further, the proposed context-aware multi-fidelity Monte Carlo method applies to hierarchies of a wide range of types of low-fidelity models such as sparse-grid and deep-network models. Numerical experiments with the gyrokinetic simulation code Gene show speedups of up to two orders of magnitude compared to standard estimators when quantifying uncertainties in small-scale fluctuations in confined plasma in fusion reactors. This corresponds to a runtime reduction from 72 days to four hours on one node of the Lonestar6 supercomputer at the Texas Advanced Computing Center.

42 ENGINEERING↗

Forward modeling fluctuations in the DESI LRGs target sample using image simulations

We use the forward modeling pipeline, Obiwan, to study the imaging systematics of the Luminous Red Galaxies (LRGs) targeted by the Dark Energy Spectroscopic Instrument (DESI). Imaging systematics refers to the false fluctuation of galaxy densities due to varying observing conditions and astrophysical foregrounds corresponding to the imaging surveys from which DESI LRG target galaxies are selected. We update the Obiwan pipeline, which we previously developed to simulate the optical images used to target DESI data, to further simulate WISE images in the infrared. This addition allows simulating the DESI LRGs sample, which utilizes WISE data in the target selection. Deep DESI imaging data combined with a method to account for biases in their shapes is used to define a truth sample of potential LRG targets. We inject these data evenly throughout the DESI Legacy Imaging Survey footprint at declinations between -30 and 32.375 degrees. We simulate a total of 15 million galaxies to obtain a simulated LRG sample (Obiwan LRGs) that predicts the variations in target density due to imaging properties. We find that the simulations predict the trends with depth observed in the data, including how they depend on the intrinsic brightness of the galaxies. We observe that faint LRGs are the main contributing source of the imaging systematics trend induced by depth. We also find significant trends in the data against Galactic extinction that are not predicted by Obiwan. These trends depend strongly on the particular map of Galactic extinction chosen to test against, implying systematic contamination in the Galactic extinction maps is a likely root cause (e.g., Cosmic-Infrared Background, dust temperature correction). We additionally observe a morphological change of the DESI LRGs population evidenced by a correlation between OII emission line average intensity and the size of the z-band PSF. This effect most likely results from uncertainties in background subtraction. The detailed findings we present should be used to guide any observational systematics mitigation treatment for the clustering of the DESI LRGs sample.

79 ASTRONOMY AND ASTROPHYSICS↗

Maskless Fourier transform holography

Fourier transform holography is a lensless imaging technique that retrieves an object's exit-wave function with high fidelity. It has been used to study nanoscale phenomena and spatio-temporal dynamics in solids, with sensitivity to the phase component of electronic and magnetic textures. However, the method requires an invasive and labor-intensive nanopatterning of a holography mask directly onto the sample, which can alter the sample properties, forces a fixed field-of-view, and leads to a low signal-to-noise ratio at high resolution. In this work, we propose using wavefront-shaping diffractive optics to create a structured probe with full control of its phase at the sample plane, circumventing the need for a mask. We demonstrate in silico that the method can image nanostructures and magnetic textures and validate our approach with a visible light-based experiment. The method enables investigation of a plethora of phenomena at the nanoscale including magnetic and electronic phase coexistence in solids, with further uses in soft and biological matter research.

Keskinbora, Kahraman↗

Enhanced Ionization Charge Identification in the Short-Baseline Neutrino Program Neutrino Detectors with Deep Neural Networks

We present a deep neural net-based region of interest detection method (DNN ROI) for signal processing in the liquid argon time projection chambers of the Short-Baseline Neutrino (SBN) Program, SBND and ICARUS. DNN ROI addresses limitations of the traditional wire-by-wire thresholding algorithm by leveraging the full two-dimensional detector readout and cross-plane matching information. To account for detector performance variations, we explore training with augmented samples. We find that DNN ROI outperforms the traditional method in both low-level ROI identification performance and high-level reconstruction metrics for high-energy cosmic and accelerator neutrino interaction products, while also being more robust against detector variations, with or without sample augmentation.

Abratenko, P. [Tufts U.]↗

Three Photon Excited Image Scanning Microscopy for in-Depth Super-Resolution Studies of Biological Samples

Multiphoton laser scanning microscopy is a powerful tool for deep imaging of thick biological samples. Image scanning microscopy (ISM) has demonstrated significant improvements in the signal-to-noise ratio in confocal laser scanning microscopy, while at the same time improving upon the effectively attainable resolution. Two-photon excitation (2PE), combined with ISM, has been shown to allow for deep tissue imaging with enhanced resolution compared to 2PE microscopy. Three-photon excitation (3PE) has enabled record imaging depth and contrast for multiphoton imaging, due to the superior suppression of out-of-focus signal generation. In this paper, we demonstrate super-resolution 3PE ISM. This is achieved using a single-photon avalanche detector array, and 1040-nm pulses for 3PE of blue fluorescence. This method enables subdiffraction limited resolution imaging of biological samples stained with blue fluorescent markers, such as mouse myocardial and spinal cord tissues stained with 4 ′ , 6 -diamidino-2-phenylindole. Deconvolution improves the resolving power further and allows for imaging with better than λ / 8 resolution with respect to the 3PE wavelength λ . With the ISM pixel reassignment procedure, we demonstrate a resolution enhancement of ∼ 1.6 laterally, compared to the resolution attained using a photomultiplier tube in a non-descanned detection arrangement, and a factor of ∼ 1.8 enhancement in axial resolution. The experimentally measured three-dimensional point spread function volume is shrunk ∼ 4.4 -fold, which is close to the theoretically expected enhancement. Published by the American Physical Society 2024

47 OTHER INSTRUMENTATION↗

Quasi-In-Situ Analysis of Electrode Top Atomic Layers via High-Sensitivity Low-Energy Ion Scattering and Potential-Controlled Sample Transfer

Electrocatalytic reactions involve interfacial interactions between the surfaces of electrodes and reactive species at an electrolyte interface. There are presently no universal or unambiguous methods to directly assay the active top atomic layer composition that influences the reactivity of these electrodes under relevant operating conditions. Low-energy ion scattering (LEIS) spectroscopy is a surface characterization technique that yields compositional analysis of the outermost atomic layer of a material, but it must be performed in ultrahigh vacuum (UHV). Application of LEIS measurements to electrochemical materials that are removed from ambient liquid-phase environments thus leaves an open question as to whether the surface that is transferred to UHV is truly the surface that manifested during the electrochemical reaction. Toward the goal of preserving the active surface state, we developed a sample transfer workflow for LEIS enabling air-free removal and drying of an electrode from an electrochemical cell while maintaining control of the potential using an auxiliary electrode. The potential-controlled emersion method was demonstrated to give distinct potential-dependent surface compositions for a Cu−Pd alloy relative to removal after uncontrolled return to open-circuit potential. A Cu-enriched surface was found at anodic potential and a Pd-enriched surface at cathodic potential, suggesting that the approach can be used to retain representative atomic configurations during transfer. Since adsorbates will often persist from the reaction environment, conventional sample pretreatment methods for removal, including atomic O and atomic H exposure, were also contrasted. Both methods were found to differ with results from incidental low-dose depth profiling by the LEIS primary ion source, which removes adventitious species and surface atoms during the course of repeated measurements. These depth profiles were found to be sensitive to sample history and thus qualitatively informative, despite the possible changes induced by ion damage. The results exhibit (i) the need for complete control over the polarization state of the sample at all times (no excursions to open circuit during transfer) and (ii) the utility of low-dose depth profiling to capture changes in the near-surface composition.

Alloys↗

Characterization of Preliminary Powder Production and Consolidation for Advanced LEU Fuel Concepts

A set of novel uranium-based alloys has been proposed for use in advanced LEU fuel systems that utilize powder metallurgical methods for fabrication of the fuel forms. This study seeks to evaluate the effectiveness of current fabrication methods to produce powder of each of the candidate alloys and consolidate composite systems. Powder samples were produced via atomization and analyzed for particle size, shape, and microstructure utilizing sieving methods and scanning electron microscopy (SEM). Preliminary pellet consolidation tests were conducted with the atomized powder using a press and sinter method. The resulting pellets were analyzed via SEM to evaluate the morphology of the composite system. These preliminary tests demonstrate the ability to produce powder but did not yield the desired quality of powder and impacted the ability to fabricate a fully dense pellet. A refined process is under development to produce a higher yield of small (<250 µm) spherical powder.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Thermal Analysis of High Entropy Rare Earth Oxides

Phase transformations in multicomponent rare earth sesquioxides were studied by splat quenching from the melt, high temperature differential thermal analysis and synchrotron X-ray diffraction on laser-heated samples. Three compositions were prepared by the solution combustion method: (La,Sm,Dy,Er,RE) 2 O 3 , where all oxides are in equimolar ratios and RE is Nd or Gd or Y. After annealing at 800 °C, all powders contained mainly a phase of C-type bixbyite structure. After laser melting, all samples were quenched in a single-phase monoclinic B-type structure. Thermal analysis indicated three reversible phase transitions in the range 1900–2400 °C, assigned as transformations into A, H, and X rare earth sesquioxides structure types. Unit cell volumes and volume changes on C-B, B-A, and H-X transformations were measured by X-ray diffraction and consistent with the trend in pure rare earth sesquioxides. The formation of single-phase solid solutions was predicted by Calphad calculations. The melting point was determined for the (La,Sm,Dy,Er,Nd) 2 O 3 sample as 2456 ± 12 °C, which is higher than for any of constituent oxides. An increase in melting temperature is probably related to nonideal mixing in the solid and/or the melt and prompts future investigation of the liquidus surface in Sm 2 O 3 -Dy 2 O 3 , Sm 2 O 3 -Er 2 O 3 , and Dy 2 O 3 -Er 2 O 3 systems.

36 MATERIALS SCIENCE↗

Neural Posterior Estimation for Cataloging Astronomical Images with Spatially Varying Backgrounds and Point Spread Functions

Neural posterior estimation (NPE), a type of amortized variational inference, is a computationally efficient means of constructing probabilistic catalogs of light sources from astronomical images. To date, NPE has not been used to perform inference in models with spatially varying covariates. However, ground-based astronomical images exhibit spatially varying sky backgrounds and point spread functions (PSFs), and accounting for this variation is essential for constructing accurate catalogs of imaged light sources. In this work, we introduce a novel NPE-based cataloging method that trains an inference network with semisynthetic astronomical images generated using PSFs and backgrounds sampled from the Sloan Digital Sky Survey. In experiments with semisynthetic images, we evaluate the method on key cataloging tasks: light source detection, star/galaxy separation, and flux measurement. A “generalist” inference network—trained with diverse PSFs and backgrounds—performs as well as a “specialist” network even when both are evaluated on the specialist’s particular PSF/background combination. This result suggests that a single NPE network can generalize across spatial variations, eliminating the need for retraining on each observational condition.

astronomy image processing↗

An ICP-OES method for the precise and accurate quantification of rare earth elements in natural water: A comparative study from mine waste sites in New Mexico, USA

Inductively coupled plasma techniques such as ICP-OES and ICP-MS are routinely used to determine the concentrations of rare earth elements (REE) in water samples. However, their performance for the determination of REE concentration in mine drainage waters from epithermal vein and porphyry copper mining districts has not been evaluated extensively. In this work, we develop an REE analysis method on an Agilent 5900 ICP-OES instrument and assess the accuracy and precision for the quantification of REE in the natural waters collected from mine adits and an acid seep of mine sites in the Steeple Rock and Hillsboro mining districts, New Mexico, USA. The total REE concentrations in the water samples were measured using the methods we developed for both ICP-OES and ICP-MS. The power of the new ICP-OES method lies in routine analysis of μg/L level concentrations normally analyzed using ICP-MS, including a U.S. Geological Survey standard reference sample, laboratory blank samples spiked with a National Institute of Standards and Technology traceable standard, and surface water samples from mine waste sites. This ICP-OES method achieves low quantification limits ranging from 0.2 to 5 μg/L and excellent analytical accuracy and precision for REE analysis. The precision of light (La-Gd) and heavy (Tb-Lu) REE analysis using this method are better than 5% at average concentrations above 5 ± 4 μg/L and 3 ± 2 μg/L, respectively, and 3% at average concentrations above 10 ± 9 μg/L and 5 ± 4 μg/L, respectively. This method also shows excellent sensitivity and reproducibility for our laboratory and field samples.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A machine learning decision criterion for reducing scan time for hyperspectral neutron computed tomography systems

We present the first machine learning-based autonomous hyperspectral neutron computed tomography experiment performed at the Spallation Neutron Source. Hyperspectral neutron computed tomography allows the characterization of samples by enabling the reconstruction of crystallographic information and elemental/isotopic composition of objects relevant to materials science. High quality reconstructions using traditional algorithms such as the filtered back projection require a high signal-to-noise ratio across a wide wavelength range combined with a large number of projections. This results in scan times of several days to acquire hundreds of hyperspectral projections, during which end users have minimal feedback. To address these challenges, a golden ratio scanning protocol combined with model-based image reconstruction algorithms have been proposed. This novel approach enables high quality real-time reconstructions from streaming experimental data, thus providing feedback to users, while requiring fewer yet a fixed number of projections compared to the filtered back projection method. In this paper, we propose a novel machine learning criterion that can terminate a streaming neutron tomography scan once sufficient information is obtained based on the current set of measurements. Our decision criterion uses a quality score which combines a reference-free image quality metric computed using a pre-trained deep neural network with a metric that measures differences between consecutive reconstructions. The results show that our method can reduce the measurement time by approximately a factor of five compared to a baseline method based on filtered back projection for the samples we studied while automatically terminating the scans.

97 MATHEMATICS AND COMPUTING↗

Evaluation of Radiography for TRISO Buffer Layer Density Measurement

Tristructural isotropic (TRISO) fuel particles consist of a central uranium-bearing kernel and a series of coating layers designed to retain fission products and to ensure fuel performance. Several parameters such as thickness and density must be measured for these coating layers to show that they conform with fuel specifications. Current methods for measuring the density of pyrolytic carbon and silicon carbide layers (liquid gradient density column) and the buffer layer (mercury porosimetry) generate Resource Conservation and Recovery Act (RCRA) radiological-mixed waste. In addition, measurement of buffer and inner pyrolytic carbon layer densities require hot sampling or interrupted coating runs and the mercury porosimetry method used for buffer density measurement only measures the mean buffer density, not the interparticle distribution. A new approach has been evaluated to measure the density of coating layers in TRISO particles based on the dependence of x-ray attenuation in radiographs on material density. This method does not generate RCRA mixed waste, measures density on a particle-by-particle basis, and in principle is capable of measuring the density of all coating layers in a single process. Initial results using thinned TRISO particle sections to evaluate radiography measurement of density as a quality control characterization method are reported herein. In this work, the primary focus is on measurement of the density of the buffer layer; however, with appropriate calibration the method should be applicable to other coating layers. Improvements to the initial method and a full demonstration of the method on the remaining coating layers may be pursued as a future effort.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Facile One-Pot Nanoproteomics for Label-Free Proteome Profiling of 50–1000 Mammalian Cells

Recent advances in sample preparation enable label-free MS-based proteome profiling of small numbers of mammalian cells. However, specific devices are often required to downscale sample processing volume from the standard 50-200 µL to sub-µL for effective nanoproteomics, which greatly impedes the implementation of current nanoproteomics methods by broad proteomics research community. Here we report a facile one-pot nanoproteomics method termed SOPs-MS (Surfactant-assisted One-Pot sample processing at the standard volume coupled with MS) for convenient proteome profiling of 50-1000 mammalian cells. Building upon our recent development of SOP-MS for label-free single-cell proteomics at low µL volume (Commun Bio 2021, 4, 265), we have systematically evaluated its processing volume at 10-200 µL using 100 human cells for robust reproducible nanoproteomic analysis. The processing volume of 50 µL which is in the range of volume for standard proteomics sample preparation, has been selected for easy sample handling with benchtop micropipette. Using the commonly accessible LC-MS platform, SOPs-MS allows for reliable label-free quantification of ~1200-2700 protein groups from 50-1000 MCF10A cells. When applied to small subpopulations of mouse colon crypt cells, SOPs-MS can reveal distinct protein signatures between any two subpopulation cells with identification of ~1500-2500 protein groups for each subpopulation. SOPs-MS may pave the way for routine deep proteome profiling of small numbers of cells as well as low-input samples.

59 BASIC BIOLOGICAL SCIENCES↗

Implementation of a practical Markov chain Monte Carlo sampling algorithm in PyBioNetFit

Abstract Summary Bayesian inference in biological modeling commonly relies on Markov chain Monte Carlo (MCMC) sampling of a multidimensional and non-Gaussian posterior distribution that is not analytically tractable. Here, we present the implementation of a practical MCMC method in the open-source software package PyBioNetFit (PyBNF), which is designed to support parameterization of mathematical models for biological systems. The new MCMC method, am, incorporates an adaptive move proposal distribution. For warm starts, sampling can be initiated at a specified location in parameter space and with a multivariate Gaussian proposal distribution defined initially by a specified covariance matrix. Multiple chains can be generated in parallel using a computer cluster. We demonstrate that am can be used to successfully solve real-world Bayesian inference problems, including forecasting of new Coronavirus Disease 2019 case detection with Bayesian quantification of forecast uncertainty. Availability and implementation PyBNF version 1.1.9, the first stable release with am, is available at PyPI and can be installed using the pip package-management system on platforms that have a working installation of Python 3. PyBNF relies on libRoadRunner and BioNetGen for simulations (e.g. numerical integration of ordinary differential equations defined in SBML or BNGL files) and Dask.Distributed for task scheduling on Linux computer clusters. The Python source code can be freely downloaded/cloned from GitHub and used and modified under terms of the BSD-3 license (https://github.com/lanl/pybnf). Online documentation covering installation/usage is available (https://pybnf.readthedocs.io/en/latest/). A tutorial video is available on YouTube (https://www.youtube.com/watch?v=2aRqpqFOiS4&t=63s). Supplementary information Supplementary data are available at Bioinformatics online.

59 BASIC BIOLOGICAL SCIENCES↗

Haar-Like Wavelets on Hierarchical Trees

Here, discrete wavelet methods, originally formulated in the setting of regularly sampled signals, can be adapted to data defined on a point cloud if some multiresolution structure is imposed on the cloud. A wide variety of hierarchical clustering algorithms can be used for this purpose, and the multiresolution structure obtained can be encoded by a hierarchical tree of subsets of the cloud. Prior work introduced the use of Haar-like bases defined with respect to such trees for approximation and learning tasks on unstructured data. This paper builds on that work in two directions. First, we present an algorithm for constructing Haar-like bases on general discrete hierarchical trees. Second, with an eye towards data compression, we present thresholding techniques for data defined on a point cloud with error controlled in the $L$ $\infty$ norm and in a Hölder-type norm. In a concluding trio of numerical examples, we apply our methods to compress a point cloud dataset, study the tightness of the $L$ $\infty$ error bound, and use thresholding to identify MNIST classifiers with good generalizability.

97 MATHEMATICS AND COMPUTING↗