Search NASA⌕ Search

SEARCH · Search NASA

Results for “Kernel 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 415 records · Page 23
All sourcesNASA NTRSDOE OSTI
DOE OSTI2012

Solution to PDEs using radial basis function finite-differences (RBF-FD) on multiple GPUs

This paper presents parallelization strategies for the radial basis function-finite difference (RBF-FD) method. As a generalized finite differencing scheme, the RBF-FD method functions without the need for underlying meshes to structure nodes. It offers high-order accuracy approximation and scales as O(N) per time step, with N being with the total number of nodes. To our knowledge, this is the first implementation of the RBF-FD method to leverage GPU accelerators for the solution of PDEs. Additionally, this implementation is the first to span both multiple CPUs and multiple GPUs. OpenCL kernels target the GPUs and inter-processor communication and synchronization is managed by the Message Passing Interface (MPI). We verify our implementation of the RBF-FD method with two hyperbolic PDEs on the sphere, and demonstrate up to 9x speedup on a commodity GPU with unoptimized kernel implementations. On a high performance cluster, the method achieves up to 7x speedup for the maximum problem size of 27,556 nodes.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Methods in PES-Learn: Direct-Fit Machine Learning of Born–Oppenheimer Potential Energy Surfaces

The release of PES-L EARN version 1.0 as an open-source software package for the automatic construction of machine learning models of semi-global molecular potential energy surfaces (PESs) is presented. Improvements to PES-L EARN ’s interoperability are stressed with new Python API that simplifies workflows for PES construction via interaction with QCSchema input and output infrastructure. In addition, a new machine learning method is introduced to PES-L EARN : kernel ridge regression (KRR). The capabilities of KRR are emphasized with examination of select semi-global PESs. All machine learning methods available in PES-L EARN are benchmarked with benzene and ethanol datasets from the rMD17 database to illustrate PES-L EARN ’s performance ability. Fitting performance and timings are assessed for both systems. Finally, the ability to predict gradients with neural network models is presented and benchmarked with ethanol and benzene. PES-L EARN is an active project and welcomes community suggestions and contributions.

kernel ridge regression↗
DOE OSTI2018

Stepwise approach to SNP-set analysis illustrated with the Metabochip and colorectal cancer in Japanese Americans of the Multiethnic Cohort

Background: Common variants have explained less than the amount of heritability expected for complex diseases, which has led to interest in less-common variants and more powerful approaches to the analysis of whole-genome scans. Because of low frequency (low statistical power), less-common variants are best analyzed using SNP-set methods such as gene-set or pathway-based analyses. However, there is as yet no clear consensus regarding how to focus in on potential risk variants following set-based analyses. We used a stepwise, telescoping approach to analyze commonand rare-variant data from the Illumina Metabochip array to assess genomic association with colorectal cancer (CRC) in the Japanese sub-population of the Multiethnic Cohort (676 cases, 7180 controls). We started with pathway analysis of SNPs that are in genes and pathways having known mechanistic roles in colorectal cancer, then focused on genes within the pathways that evidenced association with CRC, and finally assessed individual SNPs within the genes that evidenced association. Pathway SNPs downloaded from the dbSNP database were cross-matched with Metabochip SNPs and analyzed using the logistic kernel machine regression approach (logistic SNP-set kernel-machine association test, or sequence kernel association test; SKAT) and related methods. Results: The TGF-β and WNT pathways were associated with all CRC, and the WNT pathway was associated with colon cancer. Individual genes demonstrating the strongest associations were TGFBR2 in the TGF-β pathway and SMAD7 (which is involved in both the TGF-β and WNT pathways). As partial validation of our approach, a known CRC risk variant in SMAD7 (in both the TGF-β and WNT pathways: rs11874392) was associated with CRC risk in our data. We also detected two novel candidate CRC risk variants (rs13075948 and rs17025857) in TGFBR2, a gene known to be associated with CRC risk. Conclusions: A stepwise, telescoping approach identified some potentially novel risk variants associated with colorectal cancer, so it may be a useful method for following up on results of set-based SNP analyses. Further work is required to assess the statistical characteristics of the approach, and additional applications should aid in better clarifying its utility.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2015

Relativistic four-component calculations of indirect nuclear spin-spin couplings with efficient evaluation of the exchange-correlation response kernel

In this work, we report on the development and implementation of a new scheme for efficient calculation of indirect nuclear spin-spin couplings in the framework of four-component matrix Dirac-Kohn-Sham approach termed matrix Dirac-Kohn-Sham restricted magnetic balance resolution of identity for J and K, which takes advantage of the previous restricted magnetic balance formalism and the density fitting approach for the rapid evaluation of density functional theory exchange-correlation response kernels. The new approach is aimed to speedup the bottleneck in the solution of the coupled perturbed equations: evaluation of the matrix elements of the kernel of the exchange-correlation potential. The performance of the new scheme has been tested on a representative set of indirect nuclear spin-spin couplings. The obtained results have been compared with the corresponding results of the reference method with traditional evaluation of the exchange-correlation kernel, i.e., without employing the fitted electron densities. Overall good agreement between both methods was observed, though the new approach tends to give values by about 4%-5% higher than the reference method. On the average, the solution of the coupled perturbed equations with the new scheme is about 8.5 times faster compared to the reference method.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Examining the microtexture evolution in a hole-edge punched into 780 MPa grade hot-rolled steel

The deformation behavior in the hole-edge of 780 MPa grade hot-rolled steel during the punching process was investigated via microstructure characterization and computational simulation. Microstructure characterization was conducted to observe the edges of punched holes through the thickness direction, and electron back-scattered diffraction (EBSD) was used to analyze the heterogeneity of the deformation. Finite element analysis (FEA) that could account for a ductile fracture criterion was conducted to simulate the deformation and fracture behaviors of 780 MPa grade hot-rolled steel during the punching process. Calculation of rotation rate fields at the edges of the punched holes during the punching process revealed that metastable orientations in Euler space were confined to specific orientation groups. Rotation-rate fields effectively explained the stability of the initial texture components in the hole-edge region during the punching process. A visco-plastic self-consistent (VPSC) polycrystal model was used to calculate the microtexture evolution in the hole-edge region during the punching process. FEA revealed that the heterogeneous effective strain was closely related to the heterogeneity of the Kernel average misorientation (KAM) distribution in the hole-edge region. A simulation of the deformation microtexture evolution in the hole-edge region using a VPSC model was in good agreement with the experimental results. - Highlights: •We analyzed the microstructure in a hole-edge punched in HR 780HB steel. •Rotation rate fields revealed the stability of the initial texture components. •Heterogeneous effective stain was closely related to the KAM distribution. •VPSC model successfully simulated the deformation microtexture evolution.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Propagation Pattern for Moment Representation of the Lattice Boltzmann Method

A propagation pattern for the moment representation of the regularized lattice Boltzmann method (LBM) in three dimensions is presented. Using effectively lossless compression, the simulation state is stored as a set of moments of the lattice Boltzmann distribution function, instead of the distribution function itself. An efficient cache-aware propagation pattern for this moment representation has the effect of substantially reducing both the storage and memory bandwidth required for LBM simulations. This article extends recent work with the moment representation by expanding the performance analysis on central processing unit (CPU) architectures, considering how boundary conditions are implemented, and demonstrating the effectiveness of the moment representation on a graphics processing unit (GPU) architecture.

42 ENGINEERING↗
NASA NTRS1994

A spectral boundary integral equation method for the 2-D Helmholtz equation

In this paper, we present a new numerical formulation of solving the boundary integral equations reformulated from the Helmholtz equation. The boundaries of the problems are assumed to be smooth closed contours. The solution on the boundary is treated as a periodic function, which is in turn approximated by a truncated Fourier series. A Fourier collocation method is followed in which the boundary integral equation is transformed into a system of algebraic equations. It is shown that in order to achieve spectral accuracy for the numerical formulation, the nonsmoothness of the integral kernels, associated with the Helmholtz equation, must be carefully removed. The emphasis of the paper is on investigating the essential elements of removing the nonsmoothness of the integral kernels in the spectral implementation. The present method is robust for a general boundary contour. Aspects of efficient implementation of the method using FFT are also discussed. A numerical example of wave scattering is given in which the exponential accuracy of the present numerical method is demonstrated.

Hu, Fang Q.↗
DOE OSTI2016

Approximate but accurate quantum dynamics from the Mori formalism: I. Nonequilibrium dynamics

We present a formalism that explicitly unifies the commonly used Nakajima-Zwanzig approach for reduced density matrix dynamics with the more versatile Mori theory in the context of nonequilibrium dynamics. Employing a Dyson-type expansion to circumvent the difficulty of projected dynamics, we obtain a self-consistent equation for the memory kernel which requires only knowledge of normally evolved auxiliary kernels. To illustrate the properties of the current approach, we focus on the spin-boson model and limit our attention to the use of a simple and inexpensive quasi-classical dynamics, given by the Ehrenfest method, for the calculation of the auxiliary kernels. For the first time, we provide a detailed analysis of the dependence of the properties of the memory kernels obtained via different projection operators, namely, the thermal (Redfield-type) and population based (NIBA-type) projection operators. We further elucidate the conditions that lead to short-lived memory kernels and the regions of parameter space to which this program is best suited. Via a thorough analysis of the different closures available for the auxiliary kernels and the convergence properties of the self-consistently extracted memory kernel, we identify the mechanisms whereby the current approach leads to a significant improvement over the direct usage of standard semi- and quasi-classical dynamics.

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

Sparse Matrix-Matrix Multiplication on Multilevel Memory Architectures: Algorithms and Experiments

Architectures with multiple classes of memory media are becoming a common part of mainstream supercomputer deployments. So called multi-level memories offer differing characteristics for each memory component including variation in bandwidth, latency and capacity. This paper investigates the performance of sparse matrix multiplication kernels on two leading highperformance computing architectures — Intel's Knights Landing processor and NVIDIA's Pascal GPU. We describe a data placement method and a chunking-based algorithm for our kernels that exploits the existence of the multiple memory spaces in each hardware platform. We evaluate the performance of these methods w.r.t. standard algorithms using the auto-caching mechanisms Our results show that standard algorithms that exploit cache reuse performed as well as multi-memory-aware algorithms for architectures such as Ki\iLs where the memory subsystems have similar latencies. However, for architectures such as GPUS where memory subsystems differ significantly in both bandwidth and latency, multi-memory-aware methods are crucial for good performance. In addition, our new approaches permit the user to run problems that require larger capacities than the fastest memory of each compute node without depending on the software-managed cache mechanisms.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

Scalable Extended Dynamic Mode Decomposition Using Random Kernel Approximation

The Koopman operator is a linear, infinite-dimensional operator that governs the dynamics of system observables; extended dynamic mode decomposition (EDMD) is a data-driven method for approximating the Koopman operator using functions (features) of the system state snapshots. This paper investigates an approach to EDMD in which the features used provide random approximations to a particular kernel function. The objective of this is computational economy for large data sets: EDMD is generally ill-suited for problems with large state dimension, and its dual kernel formulation (KDMD) is well-suited for such problems only if the number of data snapshots is relatively small. We discuss two specific methods for generating features: random Fourier features and the Nÿstrom method. The first method is a data-independent method for translation-invariant kernels only and involves random sampling in feature space; the second method is a data-dependent empirical method that may be used for any kernel and involves random sampling of data. We first discuss how these ideas may be applied in an EDMD context, as well as a means for adaptively adding random Fourier features. As a result, we demonstrate these methods on two example problems and conclude with an analysis of the relative benefits and drawbacks of each method.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Method of Green’s Potentials for Elliptic PDEs in Domains with Random Apertures

Problems with topological uncertainties appear in many fields ranging from nano-device engineering to the design of bridges. In many of such problems, a part of the domains boundaries is subjected to random perturbations making inefficient conventional schemes that rely on discretization of the whole domain. Here, we study elliptic PDEs in domains with boundaries comprised of a deterministic part and random apertures, and apply the method of modified potentials with Green’s kernels defined on the deterministic part of the domain. This approach allows to reduce the dimension of the original differential problem by reformulating it as a boundary integral equation posed on the random apertures only. The multilevel Monte Carlo method is then applied to this modified integral equation and its optimal ϵ -2 asymptotical complexity is shown. Finally, we provide the qualitative analysis of the proposed technique and support it with numerical results.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Delegated Regressor, A Robust Approach for Automated Anomaly Detection in the Soil Radon Time Series Data

We propose a new method based on the idea of delegating regressors for predicting the soil radon gas concentration (SRGC) and anomalies in radon or any other time series data. The proposed method is compared to different traditional boosting e.g., Extreme Gradient Boosting (EGB) and simple regression methods e.g., support vector regressors with linear kernel and radial kernel in terms of accurate predictions. R language has been used for the statistical analysis of radon time series (RTS) data. The results obtained show that the proposed methodology predicts SRGC more accurately when compared to different traditional boosting and regression methods. The best correlation is found between the actual and predicted radon concentration for window size of 2 i.e., two days before and after the start of seismic activities. RTS data was collected from 05 February 2017 to 16 February 2018, including 7 seismic events recorded during the study period. Findings of study show that the proposed methodology predicts the SRGC with more precision, for all the window sizes, by overlapping predicted with the actual radon time series concentrations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Contributions From Cloud Morphological Changes to the Interannual Shortwave Cloud Feedback Based on MODIS and ISCCP Satellite Observations

The surface temperature-mediated change in cloud properties, referred to as the cloud feedback, continues to dominate the uncertainty in climate projections. A larger number of contemporary global climate models (GCMs) project a higher degree of warming than the previous generation of GCMs. This greater projected warming has been attributed to a less negative cloud feedback in the Southern Ocean. Here, we apply a novel “double decomposition method” that employs the “cloud radiative kernel” and “cloud regime” concepts, to two data sets of satellite observations to decompose the interannual cloud feedback into contributions arising from changes within and shifts between cloud morphologies. Our results show that contributions from the latter to the cloud feedback are large for certain regimes. We then focus on interpreting how both changes within and between cloud morphologies impact the shortwave cloud optical depth feedback over the Southern Ocean in light of additional observations. Results from the former cloud morphological changes reveal the importance of the wind response to warming increases low- and mid-level cloud optical thickness in the same region. Results from the latter cloud morphological changes reveal that a general shift from thick storm-track clouds to thinner oceanic low-level clouds contributes to a positive feedback over the Southern Ocean that is offset by shifts from thinner broken clouds to thicker mid- and low-level clouds. Our novel analysis can be applied to evaluate GCMs and potentially diagnose shortcomings pertaining to their physical parameterizations of particular cloud morphologies.

58 GEOSCIENCES↗
DOE OSTI2015

Ion aggregation in high salt solutions. III. Computational vibrational spectroscopy of HDO in aqueous salt solutions

The vibrational frequency, frequency fluctuation dynamics, and transition dipole moment of the O—D stretch mode of HDO molecule in aqueous solutions are strongly dependent on its local electrostatic environment and hydrogen-bond network structure. Therefore, the time-resolved vibrational spectroscopy the O—D stretch mode has been particularly used to investigate specific ion effects on water structure. Despite prolonged efforts to understand the interplay of O—D vibrational dynamics with local water hydrogen-bond network and ion aggregate structures in high salt solutions, still there exists a gap between theory and experiment due to a lack of quantitative model for accurately describing O—D stretch frequency in high salt solutions. To fill this gap, we have performed numerical simulations of Raman scattering and IR absorption spectra of the O—D stretch mode of HDO in highly concentrated NaCl and KSCN solutions and compared them with experimental results. Carrying out extensive quantum chemistry calculations on not only water clusters but also ion-water clusters, we first developed a distributed vibrational solvatochromic charge model for the O—D stretch mode in aqueous salt solutions. Furthermore, the non-Condon effect on the vibrational transition dipole moment of the O—D stretch mode was fully taken into consideration with the charge response kernel that is non-local polarizability density. From the fluctuating O—D stretch mode frequencies and transition dipole vectors obtained from the molecular dynamics simulations, the O—D stretch Raman scattering and IR absorption spectra of HDO in salt solutions could be calculated. The polarization effect on the transition dipole vector of the O—D stretch mode is shown to be important and the asymmetric line shapes of the O—D stretch Raman scattering and IR absorption spectra of HDO especially in highly concentrated NaCl and KSCN solutions are in quantitative agreement with experimental results. We anticipate that this computational approach will be of critical use in interpreting linear and nonlinear vibrational spectroscopies of HDO molecule that is considered as an excellent local probe for monitoring local electrostatic and hydrogen-bonding environment in not just salt but also other confined and crowded solutions.

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

Data-Driven Nonlinear Subspace Modeling for Prediction and Control of Molten Iron Quality Indices in Blast Furnace Ironmaking

Blast furnace (BF) in ironmaking is a nonlinear dynamic process with complicated physical-chemical reactions, where multi-phase and multi-field coupling and large time delay occur during its operation. In BF operation, the molten iron temperature (MIT) as well as Si, P and S contents of molten iron are the most essential molten iron quality (MIQ) indices, whose measurement, modeling and control have always been important issues in metallurgic engineering and automation field. This paper develops a novel data-driven nonlinear state space modeling for the prediction and control of multivariate MIQ indices by integrating hybrid modeling and control techniques. First, to improve modeling efficiency, a data-driven hybrid method combining canonical correlation analysis and correlation analysis is proposed to identify the most influential controllable variables as the modeling inputs from multitudinous factors would affect the MIQ indices. Then, a Hammerstein model for the prediction of MIQ indices is established using the LS-SVM based nonlinear subspace identification method. Such a model is further simplified by using piecewise cubic Hermite interpolating polynomial method to fit the complex nonlinear kernel function. Compared to the original Hammerstein model, this simplified model can not only significantly reduce the computational complexity, but also has almost the same reliability and accuracy for a stable prediction of MIQ indices. Last, in order to verify the practicability of the developed model, it is applied in designing a genetic algorithm based nonlinear predictive controller for multivariate MIQ indices by directly taking the established model as a predictor. Industrial experiments show the advantages and effectiveness of the proposed approach.

Blast furnace (BF)↗
DOE OSTI2024

Radiological Source Term Estimation and Isotopic Identification with Parallel Log Domain Particle Filters

This paper presents a parallel log-domain particle filtering algorithm combined with gamma spectrum unfolding to perform localization, identification, and evaluation of multiple point sources of various isotopes in an environment with attenuating obstacles. The method uses sets of precomputed attenuation kernels that map the attenuation characteristics of the environment. These kernels are specific to the energy level of a photopeak of interest. The spectral measurements are deconvolved into count measurements of each photopeak. These count measurements are fed into a set of parallel particle filters using attenuation kernels computed for that photopeak’s energy level. The individual regularized particle filters perform all likelihood calculations in the logarithmic domain to mitigate the effects of particle degeneracy. The output of each particle filter is combined to estimate which isotopes are present as well as their positions and strengths. The performance of the algorithm is characterized in a lab-scale environment using a mobile robot equipped with a gamma ray spectrometer in the presence of up to three different radioactive isotopes simultaneously. The sources were localized to within 10 cm, and their strengths were estimated within 10% of their true values. Furthermore, the isotopes were all correctly identified, and no spurious sources were reported.

42 ENGINEERING↗
DOE OSTI2016

SU-F-J-08: Quantitative SPECT Imaging of Ra-223 in a Phantom

Purpose: Ra-223 therapy of prostate cancer bone metastases is being used to treat patients routinely. However, the absorbed dose distribution at the macroscopic and microscopic scales remains elusive, due to the inability to image the small activities injected. Accurate activity quantification through imaging is essential to calculate the absorbed dose in organs and sub-units in radiopharmaceutical therapy, enabling personalized absorbed dose-based treatment planning methodologies and more effective and optimal treatments. Methods: A 22 cm diameter by 20 cm long cylindrical phantom, containing a 3.52 cm diameter sphere, was used. A total of 2.01 MBq of Ra-223 was placed in the phantom with 177.6 kBq in the sphere. Images were acquired on a dual-head Siemens Symbia T16 gamma camera using three 20% full-width energy windows and centered at 84, 154, and 269 keV (120 projections, 360° rotation, 45 s per view). We have implemented reconstruction of Ra-223 SPECT projections using OS-EM (up to 20 iterations of 10 subsets) with compensation for attenuation using CT-based attenuation maps, collimator-detector response (CDR) (including septal penetration, scatter and Pb x-ray modeling), and scatter in the patient using the effective source scatter estimation (ESSE) method. The CDR functions and scatter kernels required for ESSE were computed using the SIMIND MC simulation code. All Ra-223 photon emissions as well as gamma rays from the daughters Rn-219 and Bi-211 were modeled. Results: The sensitivity of the camera in the three combined windows was 107.3 cps/MBq. The visual quality of the SPECT images was reasonably good and the activity in the sphere was 27% smaller than the true activity. This underestimation is likely due to partial volume effect. Conclusion: Absolute quantitative Ra-223 SPECT imaging is achievable with careful attention to compensate for image degrading factors and system calibration.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2021

BASIC HARTREE-FOCK PROXY APPLICTION

This proxy application simulates the compute load and data-movement of the kernel of the Hartree-Fock method in quantum chemistry. The proxy features a simplified algorithm for computing electron repulsion integrals that is easily offloaded to GPUs.

FLETCHER, GRAHAMD↗