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 1,243 records · Page 69
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

Petascale Orbital-Free Density Functional Theory Enabled by Small-Box Algorithms

Orbital-free density functional theory (OFDFT) is a quantum-mechanics-based method that utilizes electron density as its sole variable. The main computational cost in OFDFT is the ubiquitous use of the fast Fourier transform (FFT), which is mainly adopted to evaluate the kinetic energy density functional (KEDF) and electron–electron Coulomb interaction terms. We design and implement a small-box FFT (SBFFT) algorithm to overcome the parallelization limitations of conventional FFT algorithms. We also propose real-space truncation of the nonlocal Wang–Teter KEDF kernel. The scalability of the SBFFT is demonstrated by efficiently simulating one full optimization step (electron density, energies, forces, and stresses) of 1,024,000 lithium (Li) atoms on up to 65,536 cores. We perform other tests using Li as a test material, including calculations of physical properties of different phases of bulk Li, geometry optimizations of nanocrystalline Li, and molecular dynamics simulations of liquid Li. All of the tests yield excellent agreement with the original OFDFT results, suggesting that the OFDFT-SBFFT algorithm opens the door to efficient first-principles simulations of materials containing millions of atoms.

Chemistry↗
DOE OSTI2025

SYCL for Performance Portability: Application Experience with Coupled Cluster Formalism in Quantum Chemistry on Exascale Systems

The exascale computing has brought unprecedented heterogeneity in node architectures, with systems such as Frontier and Aurora featuring diverse GPU accelerators, network connectivity among others. Ensuring performance portability across these platforms is a key challenge. To address this, we employ the SYCL programming model to develop portable, high-performance quantum chemistry workloads. As a representative application, we focus on the non-iterative Triples component of the coupled-cluster CCSD(T) method, a key driver in quantum chemistry. In this work, we report on our experience deploying SYCL-based implementations using both DPC++ and AdaptiveCPP across two flagship exascale platforms: OLCF Frontier with AMD MI250X GPUs and ALCF Aurora with Intel GPUs. Our results demonstrate that SYCL enables efficient, single-source implementations that scale to thousands of nodes, delivering performance on par with vendor-optimized HIP solutions. We highlight key insights into runtime behavior, kernel portability, and scaling characteristics, showing that SYCL offers a viable path for performance-portable computing.

Bagusetty, Abhishek [Argonne National Laboratory (↗
DOE OSTI2017

Quasi parton distributions and the gradient flow

We propose a new approach to determining quasi parton distribution functions (PDFs) from lattice quantum chromodynamics. By incorporating the gradient flow, this method guarantees that the lattice quasi PDFs are finite in the continuum limit and evades the thorny, and as yet unresolved, issue of the renormalization of quasi PDFs on the lattice. In the limit that the flow time is much smaller than the length scale set by the nucleon momentum, the moments of the smeared quasi PDF are proportional to those of the lightfront PDF. Finally, we use this relation to derive evolution equations for the matching kernel that relates the smeared quasi PDF and the light-front PDF.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2012

NEUTRONICS STUDIES OF URANIUM-BASED FULLY CERAMIC MICRO-ENCAPSULATED FUEL FOR PWRs

This study evaluates the core neutronics and fuel cycle characteristics that result from employing uranium-based fully ceramic micro-encapsulated (FCM) fuel in a pressurized water reactor (PWR). Specific PWR bundle designs with FCM fuel have been developed, which by virtue of their TRISO particle based elements, are expected to safely reach higher fuel burnups while also increasing the tolerance to fuel failures. The SCALE 6.1 code package, developed and maintained at ORNL, was the primary software employed to model these designs. Analysis was performed using the SCALE double-heterogeneous (DH) fuel modeling capabilities. For cases evaluated with the NESTLE full-core three-dimensional nodal simulator, because the feature to perform DH lattice physics branches with the SCALE/TRITON sequence is not yet available, the Reactivity-Equivalent Physical Transformation (RPT) method was used as workaround to support the full core analyses. As part of the fuel assembly design evaluations, fresh feed lattices were modeled to analyze the within-assembly pin power peaking. Also, a color-set array of assemblies was constructed to evaluate power peaking and power sharing between a once-burned and a fresh feed assembly. In addition, a parametric study was performed by varying the various TRISO particle design features; such as kernel diameter, coating layer thicknesses, and packing fractions. Also, other features such as the selection of matrix material (SiC, Zirconium) and fuel rod dimensions were perturbed. After evaluating different uranium-based fuels, the higher physical density of uranium mononitride (UN) proved to be favorable, as the parametric studies showed that the FCM particle fuel design will need roughly 12% additional fissile material in comparison to that of a standard UO2 rod in order to match the lifetime of an 18-month PWR cycle. Neutronically, the FCM fuel designs evaluated maintain acceptable design features in the areas of fuel lifetime, temperature coefficients of reactivity, as well as pin cell and assembly peaking factors. Key Words: FCM, TRISO, Uranium Mononitride, PWR

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2013

Detailed Simulations of Shock-Bifurcation and Ignition of an Argon-diluted Hydrogen/Oxygen Mixture in a Shock Tube

Detailed simulations of the bifurcation and ignition of an Argon-diluted Hydrogen/Oxygen mixture in the two-stage weak ignition regime are performed. An adaptive mesh-refinement (AMR) technique is employed to resolve all relevant physical scales that are associated with the viscous boundary-layer, the reaction front, and the shock-wave. A high-order hybrid WENO/central-differencing method is used as spatial discretization scheme, and a detailed chemical mechanism is employed to describe the combustion of the H2/O2 mixture. The operating conditions considered in this study are p = 5 bar and T = 1100 K, and fall in the third explosion limit. The computations show that the mixing of the thermally stratified fluid, carrying different momentum and enthalpy, introduces inhomogeneities in the core-region behind the reflected shock. These inhomogeneities act as localized ignition kernels. During the induction period, these kernels slowly expand and eventually transition to a detonation wave that rapidly consumes the unburned mixture.In competition with this detonation wave are the presence of secondary ignition kernels that appear in the unreacted core-region between reflected shock and detonation wave.

Ihme, Matthias↗
NASA NTRS2022

Snow Microphysical Retrieval from the NASA D3R Radar During ICE-POP 2018

A method is developed to use both polarimetric and dual-frequency radar measurements to retrieve microphysical properties of falling snow. It is applied to the Ku- and Ka-band measurements of the NASA dual-polarization, dual-frequency Doppler radar (D3R) obtained during the International Collaborative Experiments for PyeongChang 2018 Olympic and Paralympic winter games (ICE-POP 2018) field campaign and incorporates the Atmospheric Radiative Transfer Simulator (ARTS) microwave singles cattering property database for oriented particles. The retrieval uses optimal estimation to solve for several parameters that describe the particle size distribution (PSD), relative contribution of pristine, aggregate, and rimed ice species, and the orientation distribution along an entire radial simultaneously. Examination of Jacobian matrices and averaging kernels shows that the dual-wavelength ratio (DWR) measurements provide information regarding the characteristic particle size, and to a lesser extent, the rime fraction and shape parameter of the size distribution, whereas the polarimetric measurements provide information regarding the mass fraction of pristine particles and their characteristic size and orientation distribution. Thus, by combining the dual-frequency and polarimetric measurements, some ambiguities can be resolved that should allow a better determination of the PSD and bulk microphysical properties (e.g., snowfall rate) than can be retrieved from single-frequency polarimetric measurements or dual-frequency, single-polarization measurements. The D3R ICE-POP retrievals were validated using Precipitation Imaging Package (PIP) and Pluvio weighing gauge measurements taken nearby at the May Hills ground site. The PIP measures the snow PSD directly, and its measurements can be used to derived the snowfall rate (volumetric and water equivalent), mean volume-weighted particle size, and effective density, as well as particle aspect ratio and orientation. Four retrieval experiments were performed to evaluate the utility of different measurement combinations: Kuonly, DWR-only, Ku-pol, and All-obs. In terms of correlation, the volumetric snowfall rate (r D 0:95) and snow water equivalent rate (r D 0:92) were best retrieved by the Ku-pol method, while the DWR-only method had the lowest magnitude bias for these parameters (-31% and -8 %, respectively). The methods that incorporated DWR also had the best correlation to particle size (r D 0:74 and r D 0:71 for DWR-only and All-obs, respectively), although none of the methods retrieved density particularly well (r D 0:43 for Allobs). The ability of the measurements to retrieve mean aspect ratio was also inconclusive, although the polarimetric methods (Ku-pol and All-obs) had reduced biases and mean absolute error (MAE) relative to the Ku-only and DWR-only methods. The significant biases in particle size and snowfall rate appeared to be related to biases in the measured DWR, emphasizing the need for accurate DWR measurements and frequent calibration in future D3R deployments.

S Joseph Munchak↗
DOE OSTI2018

Thermionic Fuel Element with Cermet Fuel

A method of calculating and the computational results obtained for the time variation of the radial deformation of the fuel-element cladding of a single-element EGC based on cermet with a polycrystalline molybdenum matrix in which uranium dioxide particles are distributed are presented. It is shown that the maximum admissible deformation of cladding comprised of 1-mm thick, hardened, single-crystalline, tungsten-tantalum alloy with a 10 mm in diameter central channel in the kernel is reached in 10 years.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2022

SPEL: Software tool for Porting E3SM Land Model with OpenACC in a Function Unit Test Framework

Most high-end computers adopt hybrid architecture, porting a large-scale scientific code onto accelerators is necessary. The paper presents a generic method for porting large-scale scientific code onto accelerators using compiler directives within a modularized function unit test platform. We have implemented the method and designed a software tool (SPEL) to port the E3SM Land Model (ELM) onto the GPUs in the Summit computer. SPEL automatically generates GPU-ready test modules for all ELM functions, such as CanopyFlux, SoilTemperature, and EcosystemDynamics. SPEL breaks the ELM into a collection of standalone unit test programs for easy code verification and further performance improvement. We further optimize several ELM test modules with advanced techniques, including memory reduction, reconstructed parallel loops, and asynchronous GPU kernel launch. We hope our study will inspire new toolkit developments that expedite large-scale scientific code porting with compiler directives.

Schwartz, Peter↗
DOE OSTI2024

Radiation Damage Mitigation in FeCrAl Alloy at Sub-Recrystallization Temperatures

Traditional defect recovery methods rely on high-temperature annealing, often exceeding 750 °C for FeCrAl. In this study, we introduce electron wind force (EWF)-assisted annealing as an alternative approach to mitigate irradiation-induced defects at significantly lower temperatures. FeCrAl samples irradiated with 5 MeV Zr 2+ ions at a dose of 10 14 cm −2 were annealed using EWF at 250 °C for 60 s. We demonstrate a remarkable transformation in the irradiated microstructure, where significant increases in kernel average misorientation (KAM) and low-angle grain boundaries (LAGBs) typically indicate heightened defect density; the use of EWF annealing reversed these effects. X-ray diffraction (XRD) confirmed these findings, showing substantial reductions in full width at half maximum (FWHM) values and a realignment of peak positions toward their original states, indicative of stress and defect recovery. To compare the effectiveness of EWF, we also conducted traditional thermal annealing at 250 °C for 7 h, which proved less effective in defect recovery as evidenced by less pronounced improvements in XRD FWHM values.

FeCrAl alloys↗
NASA NTRS1990

Using a multifrontal sparse solver in a high performance, finite element code

We consider the performance of the finite element method on a vector supercomputer. The computationally intensive parts of the finite element method are typically the individual element forms and the solution of the global stiffness matrix both of which are vectorized in high performance codes. To further increase throughput, new algorithms are needed. We compare a multifrontal sparse solver to a traditional skyline solver in a finite element code on a vector supercomputer. The multifrontal solver uses the Multiple-Minimum Degree reordering heuristic to reduce the number of operations required to factor a sparse matrix and full matrix computational kernels (e.g., BLAS3) to enhance vector performance. The net result in an order-of-magnitude reduction in run time for a finite element application on one processor of a Cray X-MP.

King, Scott D.↗
NASA NTRS2007

NASA Langley's Approach to the Sandia's Structural Dynamics Challenge Problem

The objective of this challenge is to develop a data-based probabilistic model of uncertainty to predict the behavior of subsystems (payloads) by themselves and while coupled to a primary (target) system. Although this type of analysis is routinely performed and representative of issues faced in real-world system design and integration, there are still several key technical challenges that must be addressed when analyzing uncertain interconnected systems. For example, one key technical challenge is related to the fact that there is limited data on target configurations. Moreover, it is typical to have multiple data sets from experiments conducted at the subsystem level, but often samples sizes are not sufficient to compute high confidence statistics. In this challenge problem additional constraints are placed as ground rules for the participants. One such rule is that mathematical models of the subsystem are limited to linear approximations of the nonlinear physics of the problem at hand. Also, participants are constrained to use these models and the multiple data sets to make predictions about the target system response under completely different input conditions. Our approach involved initially the screening of several different methods. Three of the ones considered are presented herein. The first one is based on the transformation of the modal data to an orthogonal space where the mean and covariance of the data are matched by the model. The other two approaches worked solutions in physical space where the uncertain parameter set is made of masses, stiffnesses and damping coefficients; one matches confidence intervals of low order moments of the statistics via optimization while the second one uses a Kernel density estimation approach. The paper will touch on all the approaches, lessons learned, validation 1 metrics and their comparison, data quantity restriction, and assumptions/limitations of each approach. Keywords: Probabilistic modeling, model validation, uncertainty quantification, kernel density

Horta, Lucas G.↗
DOE OSTI2015

Many-body Green’s function theory for electron-phonon interactions: The Kadanoff-Baym approach to spectral properties of the Holstein dimer

We present a Kadanoff-Baym formalism to study time-dependent phenomena for systems of interacting electrons and phonons in the framework of many-body perturbation theory. The formalism takes correctly into account effects of the initial preparation of an equilibrium state and allows for an explicit time-dependence of both the electronic and phononic degrees of freedom. The method is applied to investigate the charge neutral and non-neutral excitation spectra of a homogeneous, two-site, two-electron Holstein model. This is an extension of a previous study of the ground state properties in the Hartree (H), partially self-consistent Born (Gd) and fully self-consistent Born (GD) approximations published in Säkkinen et al. [J. Chem. Phys. 143, 234101 (2015)]. Here, the homogeneous ground state solution is shown to become unstable for a sufficiently strong interaction while a symmetry-broken ground state solution is shown to be stable in the Hartree approximation. Signatures of this instability are observed for the partially self-consistent Born approximation but are not found for the fully self-consistent Born approximation. By understanding the stability properties, we are able to study the linear response regime by calculating the density-density response function by time-propagation. This amounts to a solution of the Bethe-Salpeter equation with a sophisticated kernel. The results indicate that none of the approximations is able to describe the response function during or beyond the bipolaronic crossover for the parameters investigated. Overall, we provide an extensive discussion on when the approximations are valid and how they fail to describe the studied exact properties of the chosen model system.

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

PyFly: A fast, portable aerodynamics simulator

Here, we present a fast, user-friendly implementation of a potential flow solver based on the unsteady vortex lattice method (UVLM), namely PyFly. UVLM computes the aerodynamic loads applied on lifting surfaces while capturing the unsteady effects such as the added mass forces, the growth of bound circulation, and the wake while assuming that the flow separation location is known a priori. This method is based on discretizing the body surface into a lattice of vortex rings and relies on the Biot–Savart law to construct the velocity field at every point in the simulated domain. We introduce the pointwise approximation approach to simulate the interactions of the far-field vortices to overcome the computational burden associated with the classical implementation of UVLM. The computational framework uses the Python programming language to provide an easy to handle user interface while the computational kernels are written in Fortran. The mixed language approach enables high performance regarding solution time and great flexibility concerning easiness of code adaptation to different system configurations and applications. The computational tool predicts the unsteady aerodynamic behavior of multiple moving bodies (e.g., flapping wings, rotating blades, suspension bridges) subject to incoming air. The aerodynamic simulator can also deal with enclosure effects, multi-body interactions, and B-spline representation of body shapes. Finally, we simulate different aerodynamic problems to illustrate the usefulness and effectiveness of PyFly.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

On finite-dimensional smoothed-particle Hamiltonian reductions of the Vlasov equation

The inclusion of spatial smoothing in finite-dimensional particle-based Hamiltonian reductions of the Vlasov equation and related models is considered. Here, this work investigates the underlying Hamiltonian structure of such smoothed particle-based methods for Hamiltonian systems and the small-scale regularization such methods implicitly make in approximating the continuum theory. In the context of the Vlasov–Poisson equation and other mean-field Lie–Poisson systems, of which Vlasov–Poisson is a special case, smoothing amounts to a convolutive regularization of the Hamiltonian. This regularization may be interpreted as a change of the inner product structure used to identify the dual space in the Lie–Poisson Hamiltonian formulation. In particular, the shape function used for spatial smoothing may be identified as the kernel function of a reproducing kernel Hilbert space whose inner product is used to define the Lie–Poisson Hamiltonian structure. It is likewise possible to introduce smoothing in the Vlasov–Maxwell system, but in this case the Poisson bracket must be modified rather than the Hamiltonian. The smoothing applied to the Vlasov–Maxwell system is incorporated by inserting smoothing in the map from canonical to kinematic coordinates. In the filtered system, the Lorentz force law and the current, the two terms coupling the Vlasov equation with Maxwell’s equations, are spatially smoothed.

Hamiltonian mechanics↗
DOE OSTI2013

UNVEILING THE NATURE OF THE UNIDENTIFIED GAMMA-RAY SOURCES. V. ANALYSIS OF THE RADIO CANDIDATES WITH THE KERNEL DENSITY ESTIMATION

Nearly one-third of the γ-ray sources detected by Fermi are still unidentified, despite significant recent progress in this area. However, all of the γ-ray extragalactic sources associated in the second Fermi-LAT catalog have a radio counterpart. Motivated by this observational evidence, we investigate all the radio sources of the major radio surveys that lie within the positional uncertainty region of the unidentified γ-ray sources (UGSs) at a 95% level of confidence. First, we search for their infrared counterparts in the all-sky survey performed by the Wide-field Infrared Survey Explorer (WISE) and then we analyze their IR colors in comparison with those of the known γ-ray blazars. We propose a new approach, on the basis of a two-dimensional kernel density estimation technique in the single [3.4] – [4.6] – [12] μm WISE color-color plot, replacing the constraint imposed in our previous investigations on the detection at 22 μm of each potential IR counterpart of the UGSs with associated radio emission. The main goal of this analysis is to find distant γ-ray blazar candidates that, being too faint at 22 μm, are not detected by WISE and thus are not selected by our purely IR-based methods. We find 55 UGSs that likely correspond to radio sources with blazar-like IR signatures. An additional 11 UGSs that have blazar-like IR colors have been found within the sample of sources found with deep recent Australia Telescope Compact Array observations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2010

Quantification of arterial plaque and lumen density with MDCT

Purpose: This study aimed to derive a mathematical correction function in order to normalize the CT number measurements for small volume arterial plaque and small vessel mimicking objects, imaged with multidetector CT (MDCT). Methods: A commercially available calcium plaque phantom (QRM GmbH, Moehrendorf, Germany) and a custom built cardiovascular phantom were scanned with 320 and 64 MDCT scanners. The calcium hydroxyapatite plaque phantom contained objects 0.5-5.0 mm in diameter with known CT attenuation nominal values ranging 50-800 HU. The cardiovascular phantom contained vessel mimicking objects 1.0-5.0 mm in diameter with different contrast media. Both phantoms were scanned using clinical protocols for CT angiography and images were reconstructed with different filter kernels. The measured CT number (HU) and diameter of each object were analyzed on three clinical postprocessing workstations. From the resultant data, a mathematical formula was derived based on absorption function exp(-{mu}{sup *}d) to demonstrate the relation between measured CT numbers and object diameters. Results: The percentage reduction in measured CT number (HU) for the group of selected filter kernels, apparent during CT angiography, is dependent only on the object size (plaque or vessel diameter). The derived formula of the form 1-c{sup *}exp(-a{sup *}d{sup b}) showed reduction in CT number for objects between 0.5 and 5 mm in diameter, with asymptote reaching background noise for small objects with diameters nearing the CT in-plane resolution (0.35 mm). No reduction was observed for the objects with diameters equal or larger than 5 mm. Conclusions: A clear mathematical relationship exists between object diameter and reduction in measured CT number in HU. This function is independent of exposure parameters and inherent attenuation properties of the objects studied. Future developments include the incorporation of this mathematical model function into quantification software in order to automatically generate a true assessment of measured CT number (HU) corresponding to plaque physical density {rho} (g/cm{sup 3}). This is a significant development for the accurate, noninvasive classification of noncalcified arterial plaque.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2025

Agentic AI vs ML-Based Autotuning: A Comparative Study for Loop Reordering Optimization

High Performance Computing (HPC) applications rely heavily on code optimizations to achieve good performance on modern CPU and GPU architectures. Traditional Machine Learning auto-tuning approaches have demonstrated success in exploring high-dimensional spaces, but they often require expensive compile-run evaluations and lack adaptability for large HPC applications. The recent advances in Large Language Models (LLMs) and Agentic AI systems raise intriguing questions about the potential of these approaches to address specific optimization methodologies. This work aims to answer an essential question for the HPC community: “How Agentic AI Systems Compare to Traditional ML Autotuning Techniques?” To address this question, we present a comparative analysis between a traditional ML-based optimization approach and an Agentic AI system, evaluating their respective capabilities and limitations for loop-level optimization. In addition, we introduced a new Agentic AI system named LoopGen-AI using three different Large Language Models: GPT-4.1, Claude 4.0, and Gemini 2.5. A key finding is that LoopGen-AI achieves competitive per-formance with only a few program runs, the reasoning logs from the agents revealed that their decisions rely heavily on the combination of semantic understanding of the target kernel with dynamic feedback from the environment, highlighting a promising new dimension in performance tuning. In contrast, ML-based autotuners focus on statistical exploration, and require orders of magnitude more runs to reach peak performance. Additionally, our analysis shows that prompt engineering, particularly using Persona + Context Manager patterns, significantly impacts the effectiveness of Agentic AI. Our results indicate that while Agentic AI systems are not yet a complete replacement for ML-based autotuners, it can effectively complement traditional methods.

Rosas, Miguel Romero↗
DOE OSTI2023

Real-space inversion and super-resolution of ultrafast scattering

Ultrafast scattering using x-rays or electrons is an emerging method to obtain structure dynamics at the atomic lengthscales and timescales. However, directly resolving in real-space atomic motions is inherently limited by the finite detector range and the probe energy. So, as a result, the time-resolved signal interpretation is mostly done in reciprocal space and relies on modeling and simulations of specific structures and processes. Here, we introduce a model-free approach to directly resolve scattering signals in real space, surpassing the diffraction limit, using scattering kernels and signal priors that naturally arise from the measurement constraints. We demonstrate the approach on simulated and experimental data, recover multiple atomic motions at sub-angstrom resolutions, and discuss the recovery accuracy and resolution limits versus signal fidelity. The approach offers a robust path to obtain high-resolution real-space information of atomic-scale structure dynamics using current time-resolved x-ray or electron scattering sources.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗