Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Kernel”

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 OSTI2018

Multithreaded sparse matrix-matrix multiplication for many-core and GPU architectures

Sparse matrix-matrix multiplication is a key kernel that has applications in several domains such as scientific computing and graph analysis. Several algorithms have been studied in the past for this foundational kernel. In this work, we develop parallel algorithms for sparse matrix-matrix multiplication with a focus on performance portability across different high performance computing architectures. The performance of these algorithms depend on the data structures used in them. We compare different types of accumulators in these algorithms and demonstrate the performance difference between these data structures. Furthermore, we develop a meta-algorithm, kkSpGEMM, to choose the right algorithm and data structure based on the characteristics of the problem. We show performance comparisons on three architectures and demonstrate the need for the community to develop two phase sparse matrix-matrix multiplication implementations for efficient reuse of the data structures involved.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Achieving performance portability in Gaussian basis set density functional theory on accelerator based architectures in NWChemEx

The numerical integration of the exchange–correlation (XC) potential is one of the primary computational bottlenecks in Gaussian basis set Kohn–Sham density functional theory (KS-DFT). To achieve optimal performance and accuracy, care must be taken in this numerical integration to preserve local sparsity as to allow for near linear weak scaling with system size. This leads to an integration scheme with several performance critical kernels which must be hand optimized for each architecture of interest. As the set of available accelerator hardware goes more diverse, a key challenge for developers of KS-DFT software is to maintain performance portability across a wide range of computational architectures. In this article, we examine a modular software design pattern which decouples the implementation details of performance critical kernels from the expression of high-level algorithmic workflows in a device-agnostic language such as C++; thus allowing for developers to target existing and emerging accelerator hardware within a single code base. We consider the efficacy of such a design pattern in the numerical integration of the XC potential by demonstrating its ability to achieve performance portability across a set of accelerator architectures which are representative of those on current and future U.S. Department of Energy Leadership Computing Facilities.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

High performance sparse multifrontal solvers on modern GPUs

Here, we have ported the numerical factorization and triangular solve phases of the sparse direct solver STRUMPACK to GPU. STRUMPACK implements sparse LU factorization using the multifrontal algorithm, which performs most of its operations in dense linear algebra operations on so-called frontal matrices of various sizes. Our GPU implementation off-loads these dense linear algebra operations, as well as the sparse scatter–gather operations between frontal matrices. For the larger frontal matrices, our GPU implementation relies on vendor libraries such as cuBLAS and cuSOLVER for NVIDIA GPUs and rocBLAS and rocSOLVER for AMD GPUs. For the smaller frontal matrices we developed custom CUDA and HIP kernels to reduce kernel launch overhead. Overall, high performance is achieved by identifying submatrix factorizations corresponding to sub-trees of the multifrontal assembly tree which fit entirely in GPU memory. The multi-GPU setting uses SLATE (Software for Linear Algebra Targeting Exascale) as a modern GPU-aware replacement for ScaLAPACK. On 4 nodes of SUMMIT the code runs ~10X faster when using all 24 V100 GPUs compared to when it only uses the 168 POWER9 cores. On 8 SUMMIT nodes, using 48 V100 GPUs, the sparse solver reaches over 50TFlop/s. Compared to SuperLU, on a single V100, for a set of 17 matrices our implementation is faster for all but one matrix, and is on average 5X (median 4X) faster

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

High volume packing fraction TRISO-based fuel in light water reactors

We report that for a decade, fully ceramic microencapsulated (FCM) fuel, containing tri-structural isotropic (TRISO) fuel particles in a silicon carbide (SiC) matrix, has been investigated as an accident-tolerant fuel for light water reactors (LWRs). Other examples exist of TRISO-based concepts for LWR fuels with different matrix materials. Previous studies assumed TRISO particle volume packing of approximately 0.44 in SiC (or another) matrix, the highest realistic packing fractions possible with conventional manufacturing. Recent advances in advanced manufacturing have yielded the development and demonstration of a fuel form that consists of conventionally manufactured TRISO particles in a 3D-printed SiC matrix with significantly higher possible TRISO packing fractions (0.5–0.7). This increased uranium loading enhances the viability of using TRISO-based particle fuel forms in LWRs. The viability of high-packing-fraction TRISO-based particle fuel forms in LWRs is assessed from the perspective of fuel cycle length, achievable fuel burnup, reactivity coefficients, and fuel cycle performance. Higher-packing-fraction TRISO-based fuel enables either longer cycle lengths (by ~25% at a packing fraction of 0.55 relative to 0.44) at a constant enrichment or decreased enrichments (by ~25% at a packing fraction of 0.55 relative to 0.44) at a constant cycle length. Studies of different fuel kernel types (uranium nitride, uranium oxycarbide, and uranium carbide) yield similar results, although the cycle length of uranium oxycarbide is shorter than for uranium nitride or uranium carbide (due to the lower density of uranium oxide). This work also characterized the production of 14 C resulting from neutron absorption in 14 N during operation for uranium mononitride fuel kernels; the ratio of 14 C/N was 1–2 at. % at discharge. For the fuel cycle evaluation, the activity of spent nuclear fuel and high-level waste at 100 and 100,000 years was lower for high-packing-fraction fuels than for conventional LWR fuel. Environmental impact metrics were similar overall, but higher on the front end of the fuel cycle and lower on the back end of the fuel cycle. Reactivity coefficients of higher-packing-fraction TRISO-based fuel were reasonable compared with those of conventional fuels.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2018

CFD-DEM solution verification: Fixed-bed studies

Fixed (frozen) particle configurations are used to study CFD-DEM convergence behavior in the limit of infinite grid resolution. A filtering scheme based on solving a diffusion equation, i.e., a Gaussian kernel, is used to interpolate particle information to the continuum grid, free of cell size limitations. A regression method is used for Richardson extrapolation of the discrete solutions to the grid-free solution. The regression method is shown to predict accurate extrapolation curves in the presence of the observed oscillatory convergence behavior. The extrapolation curve is used to complete the solution verification process by quantifying the numerical uncertainty present in the calculations. Several numerical and model parameters are studied including the spatial discretization method, the kernel width and the mean drag law. Here, a manageable level of discretization error is found with relative discretization errors generally <4%.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Data-Efficient Strategies for Probabilistic Voltage Envelopes under Network Contingencies

This work presents an efficient data-driven method to construct probabilistic voltage envelopes (PVE) using power flow learning in grids with network contingencies. First, a network-aware Gaussian process (GP) termed Vertex-Degree Kernel (VDK-GP), developed in prior work, is used to estimate voltage–power functions for a few network configurations. The paper introduces a novel multi-task vertex degree kernel (MT-VDK) that amalgamates the learned VDK-GPs to determine power flows for unseen networks, with a significant reduction in the computational complexity and hyperparameter requirements compared to alternate approaches. Simulations on the IEEE 30-Bus network demonstrate the retention and transfer of power flow knowledge in both N-1 and N-2 contingency scenarios. The MT-VDK-GP approach achieves over 50 % reduction in mean prediction error for novel N-1 contingency network configurations in low training data regimes (50–250 samples) over VDK-GP. Additionally, MT-VDK-GP outperforms a hyper-parameter based transfer learning approach in over 75 % of N-2 contingency network structures, even without historical N-2 outage data. Furthermore, the proposed method demonstrates the ability to achieve PVEs using sixteen times fewer power flow solutions compared to Monte-Carlo sampling-based methods.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2023

Seeing the whole picture: Methods for getting the most from micro X-ray computed tomography of TRISO nuclear fuel particles

Tristructural isotropic (TRISO) coated fuel particles are a nuclear fuel form under extensive study for use in advanced nuclear reactor concepts. TRISO fuels are subjected to high temperature neutron irradiations and then examined to assess their performance by determining fission product retention and studying morphological changes. Micro X-ray computed tomography is one method of nondestructively studying the effects of TRISO performance. This work addresses the need for image processing to remove X-ray tomographic reconstruction artifacts that prevent the study of TRISO features, as the TRISO particles’ high Z kernel can introduce metal artifacts that degrade the image quality in the surrounding low Z coating layers. These metal artifacts were reduced by imaging the TRISO particles with both high- and low-energy X-rays and applying a mask to the radiographs obtained with low-energy X-rays to digitally remove the dense fuel kernel region. These masked radiographs were then used to produce a tomographic reconstruction which was combined with the tomographic reconstruction of the high-energy data. This enabled the relatively-low-density TRISO buffer layer to be examined in more detail, providing information on irradiation induced dimensional changes of the coatings. This methodology, which helps see the full picture of a TRISO particle, is not limited to nuclear fuels but can be applied to systems that contain highly attenuating material surrounded by less dense materials.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2015

A new classification scheme of plastic wastes based upon recycling labels

Highlights: • PET, HPDE or PP types of plastics are considered. • An automated classification of plastic bottles based on the feature extraction and classification methods is performed. • The decision mechanism consists of PCA, Kernel PCA, FLDA, SVD and Laplacian Eigenmaps methods. • SVM is selected to achieve the classification task and majority voting technique is used. - Abstract: Since recycling of materials is widely assumed to be environmentally and economically beneficial, reliable sorting and processing of waste packaging materials such as plastics is very important for recycling with high efficiency. An automated system that can quickly categorize these materials is certainly needed for obtaining maximum classification while maintaining high throughput. In this paper, first of all, the photographs of the plastic bottles have been taken and several preprocessing steps were carried out. The first preprocessing step is to extract the plastic area of a bottle from the background. Then, the morphological image operations are implemented. These operations are edge detection, noise removal, hole removing, image enhancement, and image segmentation. These morphological operations can be generally defined in terms of the combinations of erosion and dilation. The effect of bottle color as well as label are eliminated using these operations. Secondly, the pixel-wise intensity values of the plastic bottle images have been used together with the most popular subspace and statistical feature extraction methods to construct the feature vectors in this study. Only three types of plastics are considered due to higher existence ratio of them than the other plastic types in the world. The decision mechanism consists of five different feature extraction methods including as Principal Component Analysis (PCA), Kernel PCA (KPCA), Fisher’s Linear Discriminant Analysis (FLDA), Singular Value Decomposition (SVD) and Laplacian Eigenmaps (LEMAP) and uses a simple experimental setup with a camera and homogenous backlighting. Due to the giving global solution for a classification problem, Support Vector Machine (SVM) is selected to achieve the classification task and majority voting technique is used as the decision mechanism. This technique equally weights each classification result and assigns the given plastic object to the class that the most classification results agree on. The proposed classification scheme provides high accuracy rate, and also it is able to run in real-time applications. It can automatically classify the plastic bottle types with approximately 90% recognition accuracy. Besides this, the proposed methodology yields approximately 96% classification rate for the separation of PET or non-PET plastic types. It also gives 92% accuracy for the categorization of non-PET plastic types into HPDE or PP.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2018

Differential formulation of the viscous history force on a particle for efficient and accurate computation

It is well known that the computation of the Basset-like history force is very demanding in terms of CPU and memory requirements, since it requires the evaluation of a history integral. We use the recent rational theory of Beylkin & Monzón ( Appl. Comput. Harmon. Anal. , vol. 19, 2005, pp. 17–48) to approximate the history kernel in the form of exponential sums to reformulate the viscous history force in a differential form. This theory allows us to approximate the history kernel in terms of exponential sums to any desired order of accuracy. This removes the need for long-time storage of the acceleration histories of the particle and the fluid. The proposed differential form approximation is applied to compute the history force on a spherical particle in a synthetic turbulent flow and a wall-bounded turbulent channel flow. Particles of various diameters are considered, and results obtained using the present technique are in reasonable agreement with those achieved using the full history integral.

Mechanics↗
DOE OSTI2024

Non-locality of mean scalar transport in two-dimensional Rayleigh–Taylor instability using the macroscopic forcing method

The importance of non-locality of mean scalar transport in two-dimensional Rayleigh–Taylor Instability (RTI) is investigated. The macroscopic forcing method is utilized to measure spatio-temporal moments of the eddy diffusivity kernel representing passive scalar transport in the ensemble averaged fields. Presented in this work are several studies assessing the importance of the higher-order moments of the eddy diffusivity, which contain information about non-locality, in models for RTI. First, it is demonstrated through a comparison of leading-order models that a purely local eddy diffusivity is insufficient to capture the mean field evolution of the mass fraction in RTI. Therefore, higher-order moments of the eddy diffusivity operator are not negligible. Models are then constructed by utilizing the measured higher-order moments. It is demonstrated that an explicit operator based on the Kramers–Moyal expansion of the eddy diffusivity kernel is insufficient. An implicit operator construction that matches the measured moments is shown to offer improvements relative to the local model in a converging fashion.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

Acquisition of STEM Images by Adaptive Compressive Sensing

Compressive Sensing (CS) allows a signal to be sparsely measured first and accurately recovered later in software [1]. In scanning transmission electron microscopy (STEM), it is possible to compress an image spatially by reducing the number of measured pixels, which decreases electron dose and increases sensing speed [2,3,4]. The two requirements for CS to work are: (1) sparsity of basis coefficients and (2) incoherence of the sensing system and the representation system. However, when pixels are missing from the image, it is difficult to have an incoherent sensing matrix. Nevertheless, dictionary learning techniques such as Beta-Process Factor Analysis (BPFA) [5] are able to simultaneously discover a basis and the sparse coefficients in the case of missing pixels. On top of CS, we would like to apply active learning [6,7] to further reduce the proportion of pixels being measured, while maintaining image reconstruction quality. Suppose we initially sample 10% of random pixels. We wish to select the next 1% of pixels that are most useful in recovering the image. Now, we have 11% of pixels, and we want to decide the next 1% of “most informative” pixels. Active learning methods are online and sequential in nature. Our goal is to adaptively discover the best sensing mask during acquisition using feedback about the structures in the image. In the end, we hope to recover a high quality reconstruction with a dose reduction relative to the non-adaptive (random) sensing scheme. In doing this, we try three metrics applied to the partial reconstructions for selecting the new set of pixels: (1) variance, (2) Kullback-Leibler (KL) divergence using a Radial Basis Function (RBF) kernel, and (3) entropy. Figs. 1 and 2 display the comparison of Peak Signal-to-Noise (PSNR) using these three different active learning methods at different percentages of sampled pixels. At 20% level, all the three active learning methods underperform the original CS without active learning. However, they all beat the original CS as more of the “most informative” pixels are sampled. One can also argue that CS equipped with active learning requires less sampled pixels to achieve the same value of PSNR than CS with pixels randomly sampled, since all the three PSNR curves with active learning grow at a faster pace than that without active learning. For this particular STEM image, by observing the reconstructed images and the sensing masks, we find that while the method based on RBF kernel acquires samples more uniformly, the one on entropy samples more areas of significant change, thus less uniformly. The KL-divergence method performs the best in terms of reconstruction error (PSNR) for this example [8].

Xie, Weiyi↗
DOE OSTI2019

Ab Initio Computation of Rotationally-Averaged Pump–Probe X-ray and Electron Diffraction Signals

We develop a new algorithm for the computation of the rotationally averaged elastic molecular diffraction signal for the cases of perpendicular or parallel pump–probe geometries. The algorithm first collocates the charge density from an arbitrary ab initio wave function onto a Becke quadrature grid [A. Becke, J. Chem. Phys. 1988, 88, 2457], providing a high-fidelity multiresolution representation of the charge density. A double sum is then performed over the Becke grid points, and the interaction between points computed using the scattering kernels of Williamson and Zewail [J. C. Williamson and A. H. Zewail, J. Phys. Chem. 1994, 98, 2766]. These kernels analytically average over the molecular orientations with the cos 2 γ selection factor appropriate for one-photon dipole absorption in a perpendicular pump–probe geometry. We demonstrate that the method is converged with small grids containing <500 points/atom. We implement the algorithm on a GPU for increased efficiency and emonstrate the algorithm for molecules with up to a few dozen atoms. We explore the accuracy of the independent atom model (IAM) by comparison with our new and more accurate method. We also investigate the possibility of detecting signatures of electronic transitions in polyatomic pump–probe diffraction experiments.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Spatial Signatures of Electron Correlation in Least-Squares Tensor Hypercontraction

Least Squares Tensor Hypercontraction (LS-THC) has received some attention in recent years as an approach to reduce the significant computational costs of wavefunc- tion based methods in quantum chemistry. However, previous work has demonstrated that the LS-THC factorization performs disproportionately worse in the description of wavefunction components (e.g. cluster amplitudes T 2 ) than Hamiltonian compo- nents (e.g. electron repulsion integrals (pq|rs)). This work develops novel theoretical methods to study the source of these errors in the context of the real-space T 2 kernel, and reports, for the first time, the existence of a “correlation feature” in the errors of the LS-THC representation of the “exchange-like” correlation energy EX and T 2 that is remarkably consistent across ten molecular species, three correlated wavefunctions, and four basis sets. This correlation feature portends the existence of a “pair-point kernel” missing in the usual LS-THC representation of the wavefunction, which critically depends upon pairs of grid points situated close to atoms and with inter-pair distances between one and two Bohr radii. These findings point the way for future LS-THC developments to address these shortcomings.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Source of Bright Near-Infrared Luminescence in Gold Nanoclusters

Gold nanoclusters with near-infrared (NIR) photoluminescence (PL) have great potential as sensing and imaging materials in biomedical and bioimaging applications. In this work, Au 21 (S-Adm) 15 and Au 38 S 2 (S-Adm) 20 are used to unravel the underlying mechanisms for the improved quantum yields (QY), large Stokes shifts and long PL lifetimes in gold nanoclusters. Both nanoclusters show decent PL QY. In particular, the Au 38 S 2 (S-Adm) 20 nanocluster shows a bright NIR PL at 900 nm with QY up to 15% in normal solvents (such as toluene) at ambient conditions. The relatively lower QY for Au 21 (S-Adm) 15 (4%) compared to Au 38 S 2 (S-Adm) 20 is attributed to the lowest-lying excited state being symmetry-disallowed, as evidenced by the pressure-dependent anti-spectral shift of the absorption spectra compared to PL. Yet, Au 21 (S-Adm) 15 maintains some emissive properties due to a nearby symmetry-allowed excited state. Furthermore, our results show that suppression of non-radiative decay due to the surface “lock rings” which encircle the Au kernel and the surface “lock atoms” which bridge the fundamental Au-kernel units (e.g., tetrahedra, icosahedra, etc.) is the key to obtain high QYs in gold nanoclusters. Here, the complicated excited-state processes and the small absorption coefficient of the band-edge transition lead to the large Stokes shifts and the long PL lifetimes that are widely observed in gold nanoclusters.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2018

Local radiative feedbacks over the Arctic based on observed short-term climate variations

We compare various radiative feedbacks in response to short-term climate variations estimated from reanalysis, satellite, and global climate model datasets using the combined Kernel-Gregory approach over the Arctic (60-90°N). Except for the shortwave cloud feedback, all feedbacks are positive over the Arctic. Lapse rate feedback is comparable with surface albedo feedback. Relative to the tropics (30°S-30°N), the temperature feedback (lapse rate feedback plus Planck feedback deviation from its global mean) is the largest contributor to Arctic amplification among all feedbacks, followed by surface albedo feedback. The choice of meteorological fields has a small impact on the cloud feedback estimate, but different top-of-the-atmosphere clear-sky fluxes datasets can result in a wide spread of estimated cloud feedbacks. The net cloud feedback over the Arctic is positive due to a larger positive in longwave than the negative shortwave feedback, but has the largest uncertainty. Relative to the tropics, the longwave cloud feedback largely contributes to the Arctic amplification.

Zhang, Rudong↗
DOE OSTI2020

Assessing Global and Local Radiative Feedbacks Based on AGCM Simulations for 1980–2014/2017

In order to avoid contamination of unobservable and uncertain effective radiative forcing (ERF) on the diagnosis of radiative feedbacks based on short-term climate variability, we verify the Kernel-Gregory feedback calculation method using atmospheric model experiments with prescribed ERF. We show that both clear-sky radiative fluxes and all-sky radiative feedbacks have a closure between model simulations and kernel derivations. A near-zero global mean net cloud feedback is found in the simulation with prescribed ERF, which results in a more negative global net climate feedback, -2 W m -2 K -1 . Consistent with AMIP6 ensemble mean results, the lapse rate feedback is the largest contributor among all feedbacks to temperature amplification over the three poles (Arctic, Antarctic and Tibetan Plateau), followed by surface albedo feedback and Planck feedback deviation from its global mean. Except for the higher surface albedo feedback in the Antarctic, other feedbacks are almost same between the Arctic and Antarctic.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Bayesian force fields from active learning for simulation of inter-dimensional transformation of stanene

Abstract We present a way to dramatically accelerate Gaussian process models for interatomic force fields based on many-body kernels by mapping both forces and uncertainties onto functions of low-dimensional features. This allows for automated active learning of models combining near-quantum accuracy, built-in uncertainty, and constant cost of evaluation that is comparable to classical analytical models, capable of simulating millions of atoms. Using this approach, we perform large-scale molecular dynamics simulations of the stability of the stanene monolayer. We discover an unusual phase transformation mechanism of 2D stanene, where ripples lead to nucleation of bilayer defects, densification into a disordered multilayer structure, followed by formation of bulk liquid at high temperature or nucleation and growth of the 3D bcc crystal at low temperature. The presented method opens possibilities for rapid development of fast accurate uncertainty-aware models for simulating long-time large-scale dynamics of complex materials.

Chemistry↗