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,009 records · Page 56
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

AuriDESI: mock catalogues for the DESI Milky Way Survey

The Dark Energy Spectroscopic Instrument Milky Way Survey (DESI MWS) will explore the assembly history of the Milky Way by characterizing remnants of ancient dwarf galaxy accretion events and improving constraints on the distribution of dark matter in the outer halo. We present mock catalogues that reproduce the selection criteria of MWS and the format of the final MWS data set. These catalogues can be used to test methods for quantifying the properties of stellar halo substructure and reconstructing the Milky Way’s accretion history with the MWS data, including the effects of halo-to-halo variance. The mock catalogues are based on a phase-space kernel expansion technique applied to star particles in the Auriga suite of six high-resolution lambda-cold dark matter magnetohydrodynamic zoom-in simulations. They include photometric properties (and associated errors) used in DESI target selection and the outputs of the MWS spectral analysis pipeline (radial velocity, metallicity, surface gravity, and temperature). They also include information from the underlying simulation, such as the total gravitational potential and information on the progenitors of accreted halo stars. We discuss how the subset of halo stars observable by MWS in these simulations corresponds to their true content and properties. These mock Milky Ways have rich accretion histories, resulting in a large number of substructures that span the whole stellar halo out to large distances and have substantial overlap in the space of orbital energy and angular momentum.

dynamics↗
DOE OSTI2024

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↗
NASA NTRS1972

Development and verification of design methods for ducts in a space nuclear shield

A practical method for computing the effectiveness of a space nuclear shield perforated by small tubing and cavities is reported. Performed calculations use solutions for a two dimensional transport code and evaluate perturbations of that solution using last flight estimates and other kernel integration techniques. In general, perturbations are viewed as a change in source strength of scattered radiation and a change in attenuation properties of the region.

Cerbone, R. J.↗
NASA NTRS1986

Aerodynamics via acoustics - Application of acoustic formulas for aerodynamic calculations

Prediction of aerodynamic loads on bodies in arbitrary motion is considered from an acoustic point of view, i.e., in a frame of reference fixed in the undisturbed medium. An inhomogeneous wave equation which governs the disturbance pressure is constructed and solved formally using generalized function theory. When the observer is located on the moving body surface there results a singular linear integral equation for surface pressure. Two different methods for obtaining such equations are discussed. Both steady and unsteady aerodynamic calculations are considered. Two examples are presented, the more important being an application to propeller aerodynamics. Of particular interest for numerical applications is the analytical behavior of the kernel functions in the various integral equations.

Farassat, F.↗
NASA NTRS1986

Aerodynamics Via Acoustics: Application of Acoustic Formulas for Aerodynamic Calculations

Prediction of aerodynamic loads on bodies in arbitrary motion is considered from an acoustic point of view, i.e., in a frame of reference fixed in the undisturbed medium. An inhomogeneous wave equation which governs the disturbance pressure is constructed and solved formally using generalized function theory. When the observer is located on the moving body surface there results a singular linear integral equation for surface pressure. Two different methods for obtaining such equations are discussed. Both steady and unsteady aerodynamic calculations are considered. Two examples are presented, the more important being an application to propeller aerodynamics. Of particular interest for numerical applications is the analytical behavior of the kernel functions in the various integral equations.

Farassat, F.↗
DOE OSTI2022

Just-in-Time Compilation and Link-Time Optimization for OpenMP Target Offloading

Following the mass adoption of external accelerators for high performance computing, the overall performance of many applications has become increasingly dependent on relatively small accelerated kernels. As static analysis is fundamentally limited by dynamic values and external definitions, standard ahead-of-time compilation is not always sufficient to achieve the best performance. Furthermore, many users looking to port an existing application to run on an external accelerator will not want to fundamentally restructure their programs. These and other problems can be addressed through both link-time optimization (LTO) and just-in-time (JIT) compilation, but until now had sparse and inconsistent support from the compiler. In this work, we present a new compilation method that enables device-side LTO as well as a transparent JIT compilation tool-chain for OpenMP target offloading. Our contributions include an entirely new device linking and embedding scheme to enable LTO as well as a novel JIT engine to efficiently optimize OpenMP offloading regions at run-time. We also introduce a persistent caching system to improve end-to-end runtime using the JIT engine and minimize kernel launching overheads. We measure the performance of our LTO and JIT implementation via several real-world scientific applications. With our optimizations we observe significant improvements through LTO on large applications as well as significant end-to-end execution time improvement using JIT.

Tian, Shilei↗
DOE OSTI2016

SELF-GRAVITATIONAL FORCE CALCULATION OF INFINITESIMALLY THIN GASEOUS DISKS ON NESTED GRIDS

We extend the work of Yen et al. and develop second-order formulae to accommodate a nested grid discretization for the direct self-gravitational force calculation for infinitesimally thin gaseous disks. This approach uses a two-dimensional kernel that is derived for infinitesimally thin disks and is free of artificial boundary conditions. The self-gravitational force calculation is presented in generalized convolution forms for a nested grid configuration. A numerical technique derived from a fast Fourier transform is employed to reduce the computational complexity to be nearly linear. By comparing with analytic potential–density pairs associated with the generalized Maclaurin disks, the extended approach is verified to be of second-order accuracy when using numerical simulations. The proposed method is accurate, computationally fast, and has the potential to be applied to studies of planetary migration and the gaseous morphology of disk galaxies.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Data-driven Minimum Entropy Control for Stochastic Nonlinear Systems using the Cumulant-Generating Function

Here, we present a novel minimum entropy control algorithm for a class of stochastic nonlinear systems subjected to non-Gaussian noises. The entropy control can be considered as an optimization problem for the system randomness attenuation, but the mean value has to be considered separately. To overcome this disadvantage, a new representation of the system stochastic properties was given using the cumulant-generating function based on the moment-generating function, in which the mean value and the entropy was reflected by the shape of the cumulant-generating function. Based on the samples of the system output and control input, a time-variant linear model was identified, and the minimum entropy optimization was transformed to system stabilization. Then, an optimal control strategy was developed to achieve the randomness attenuation, and the boundedness of the controlled system output was analyzed. The effectiveness of the presented control algorithm was demonstrated by a numerical example. In this paper, a data-driven minimum entropy design is presented without pre-knowledge of the system model; entropy optimization is achieved by the system stabilization approach in which the stochastic distribution control and minimum entropy are unified using the same identified structure; and a potential framework is obtained since all the existing system stabilization methods can be adopted to achieve the minimum entropy objective.

42 ENGINEERING↗
DOE OSTI2022

Physics Discovery in Nanoplasmonic Systems via Autonomous Experiments in Scanning Transmission Electron Microscopy

Abstract Physics‐driven discovery in an autonomous experiment has emerged as a dream application of machine learning in physical sciences. Here, this work develops and experimentally implements a deep kernel learning (DKL) workflow combining the correlative prediction of the target functional response and its uncertainty from the structure, and physics‐based selection of acquisition function, which autonomously guides the navigation of the image space. Compared to classical Bayesian optimization (BO) methods, this approach allows to capture the complex spatial features present in the images of realistic materials, and dynamically learn structure–property relationships. In combination with the flexible scalarizer function that allows to ascribe the degree of physical interest to predicted spectra, this enables physical discovery in automated experiment. Here, this approach is illustrated for nanoplasmonic studies of nanoparticles and experimentally implemented in a truly autonomous fashion for bulk‐ and edge plasmon discovery in MnPS 3 , a lesser‐known beam‐sensitive layered 2D material. This approach is universal, can be directly used as‐is with any specimen, and is expected to be applicable to any probe‐based microscopic techniques including other STEM modalities, scanning probe microscopies, chemical, and optical imaging.

42 ENGINEERING↗
DOE OSTI2021

Determining the leading-order contact term in neutrinoless double β decay

We present a method to determine the leading-order (LO) contact term contributing to the nn → ppe - e - amplitude through the exchange of light Majorana neutrinos. Our approach is based on the representation of the amplitude as the momentum integral of a known kernel (proportional to the neutrino propagator) times the generalized forward Compton scattering amplitude n ( p 1 ) n ( p 2 ) W + ( k ) → \( p\left({p}_1^{\prime}\right)p\left({p}_2^{\prime}\right){W}^{-}(k) \) , in analogy to the Cottingham formula for the electromagnetic contribution to hadron masses. We construct model-independent representations of the integrand in the low- and high-momentum regions, through chiral EFT and the operator product expansion, respectively. We then construct a model for the full amplitude by interpolating between these two regions, using appropriate nucleon factors for the weak currents and information on nucleon-nucleon ( NN ) scattering in the 1 S 0 channel away from threshold. By matching the amplitude obtained in this way to the LO chiral EFT amplitude we obtain the relevant LO contact term and discuss various sources of uncertainty. We validate the approach by computing the analog I = 2 NN contact term and by reproducing, within uncertainties, the charge-independence-breaking contribution to the 1 S 0 NN scattering lengths. While our analysis is performed in the \( \overline{\mathrm{MS}} \) scheme, we express our final result in terms of the scheme-independent renormalized amplitude \( {\mathcal{A}}_{\nu}\left(\left|\mathbf{p}\right|,\left|\mathbf{p}^{\prime}\right|\right) \) at a set of kinematic points near threshold. We illustrate for two cutoff schemes how, using our synthetic data for \( {\mathcal{A}}_{\nu } \) , one can determine the contact-term contribution in any regularization scheme, in particular the ones employed in nuclear-structure calculations for isotopes of experimental interest.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

A compute-bound formulation of Galerkin model reduction for linear time-invariant dynamical systems

This work aims to advance computational methods for projection-based reduced-order models (ROMs) of linear time-invariant (LTI) dynamical systems. For such systems, current practice relies on ROM formulations expressing the state as a rank-1 tensor (i.e., a vector), leading to computational kernels that are memory bandwidth bound and, therefore, ill-suited for scalable performance on modern architectures. This weakness can be particularly limiting when tackling many-query studies, where one needs to run a large number of simulations. This work introduces a reformulation, called rank-2 Galerkin, of the Galerkin ROM for LTI dynamical systems which converts the nature of the ROM problem from memory bandwidth to compute bound. We present the details of the formulation and its implementation, and demonstrate its utility through numerical experiments using, as a test case, the simulation of elastic seismic shear waves in an axisymmetric domain. We quantify and analyze performance and scaling results for varying numbers of threads and problem sizes. In conclusion, we present an end-to-end demonstration of using the rank-2 Galerkin ROM for a Monte Carlo sampling study. We show that the rank-2 Galerkin ROM is one order of magnitude more efficient than the rank-1 Galerkin ROM (the current practice) and about 970 times more efficient than the full-order model, while maintaining accuracy in both the mean and statistics of the field.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2003

Presumed PDF Modeling of Early Flame Propagation in Moderate to Intense Turbulence Environments

The present paper describes the results obtained from a one-dimensional time dependent numerical technique that simulates early flame propagation in a moderate to intense turbulent environment. Attention is focused on the development of a spark-ignited, premixed, lean methane/air mixture with the unsteady spherical flame propagating in homogeneous and isotropic turbulence. A Monte-Carlo particle tracking method, based upon the method of fractional steps, is utilized to simulate the phenomena represented by a probability density function (PDF) transport equation. Gaussian distributions of fluctuating velocity and fuel concentration are prescribed. Attention is focused on three primary parameters that influence the initial flame kernel growth: the detailed ignition system characteristics, the mixture composition, and the nature of the flow field. The computational results of moderate and intense isotropic turbulence suggests that flames within the distributed reaction zone are not as vulnerable, as traditionally believed, to the adverse effects of increased turbulence intensity. It is also shown that the magnitude of the flame front thickness significantly impacts the turbulent consumption flame speed. Flame conditions studied have fuel equivalence ratio s in the range phi = 0.6 to 0.9 at standard temperature and pressure.

Carmen, Christina↗
DOE OSTI2011

Neutron Cross Section Covariances for Structural Materials and Fission Products

We describe neutron cross section covariances for 78 structural materials and fission products produced for the new US evaluated nuclear reaction library ENDF/B-VII.1. Neutron incident energies cover full range from 10{sup -5} eV to 20 MeV and covariances are primarily provided for capture, elastic and inelastic scattering as well as (n,2n). The list of materials follows priorities defined by the Advanced Fuel Cycle Initiative, the major application being data adjustment for advanced fast reactor systems. Thus, in addition to 28 structural materials and 49 fission products, the list includes also {sup 23}Na which is important fast reactor coolant. Due to extensive amount of materials, we adopted a variety of methodologies depending on the priority of a specific material. In the resolved resonance region we primarily used resonance parameter uncertainties given in Atlas of Neutron Resonances and either applied the kernel approximation to propagate these uncertainties into cross section uncertainties or resorted to simplified estimates based on integral quantities. For several priority materials we adopted MF32 covariances produced by SAMMY at ORNL, modified by us by adding MF33 covariances to account for systematic uncertainties. In the fast neutron region we resorted to three methods. The most sophisticated was EMPIRE-KALMAN method which combines experimental data from EXFOR library with nuclear reaction modeling and least-squares fitting. The two other methods used simplified estimates, either based on the propagation of nuclear reaction model parameter uncertainties or on a dispersion analysis of central cross section values in recent evaluated data files. All covariances were subject to quality assurance procedures adopted recently by CSEWG. In addition, tools were developed to allow inspection of processed covariances and computed integral quantities, and for comparing these values to data from the Atlas and the astrophysics database KADoNiS.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2015

MO-FG-204-09: High Spatial Resolution and Artifact-Free CT Bone Imaging at Off-Centered Positions: An Application of Model-Based Iterative Reconstruction

Purpose: Although the anatomy of interest should be positioned as close as possible to the isocenter of CT scanners, off-centering may be inevitable during certain exams in clinical practice such as lumbar spine and elbow imaging. Off-centering degrades image sharpness, generates streak artifacts, and sometimes creates blooming artifacts due to truncation. The purpose of this work was to investigate whether the use of model-based image reconstruction (MBIR) can alleviate the negative impacts of off-centering to achieve high quality CT bone imaging. Methods: Both an anthropomorphic phantom and an ex vivo swine elbow sample were scanned at centered and off-centered positions using clinical CT bone scan protocols. The magnitude of off-centering was determined from localizer radiographs. Both FBP and MBIR reconstructions were performed. For FBP, both standard and Bone Plus kernels commonly used in bone imaging were used. Objective assessment of image sharpness, noise standard deviation, and noise nonuniformity were performed. Additionally, we retrospectively analyzed human subject data acquired under off-centered conditions as a validation study. Results: In FBP images of the phantom, off-centering by 10 cm led to a 14% increase in noise (p<1e-3) and a 68% increase in noise nonuniformity (p<0.02). A visible drop in bone sharpness was observed. In contrast, no significant difference in the noise magnitude or the noise nonuniformity between the centered and off-centered MBIR images was found. The image sharpness of off-centered MBIR images outperformed that of FBP images reconstructed with the Bone Plus kernel. In images of the swine elbow off-centered by 20 cm, not only was the noise and spatial resolution performance improved by MBIR, truncation artifacts were also elliminated. The human subject study generated similar results, in which the visibility of the off-centered lumbar spine was significantly improved. Conclusion: High quality CT bone imaging at off-centered positions can be achieved using MBIR. This work was partially supported by an NIH grant R01CA169331 and GE Healthcare. K. Li, D. Gomez-Cardona: Nothing to disclose. G.-H. Chen: Research funded, GE Healthcare; Research funded, Siemens AX. A. Budde, J. Hsieh: Employee, GE Healthcare.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

Nuclear Forensics for Irradiated Materials using Kernel Regression Models

The ultimate purpose of nuclear forensics is to locate the origin of nuclear and other radioactive materials out of regulatory control used in nuclear security events such as classical forensics. For this purpose, some specific characteristics of interdicted materials (referred to as 'signatures') should be defined to compare the known sample database. In some cases, involved with irradiated material, the signatures could be specific amount of fission product inventory depending on reactor operational histories such as reactor type, initial enrichment, burn up rate of nuclear fuel and cooling period after irradiation in the reactor. In order to make inference model for irradiated material for nuclear forensics, we suggested 3 ways as follows; (i) analyzing ORIGEN produced nuclear substances for the fission product inventory, (ii) classifying the range of reactor type, (iii) predicting the operational history from the fission product inventory. In terms of nuclear substance analysis, the construction of spent fuel database has been internationally carried out since ITWG (Nuclear Smuggling International Technical Working Group) was established in 1995. However, enough sample data cannot be obtained in Korea to determine reactor type and operational histories because analyzing nuclear substances have been limited due to the international regulation. Thus, this study does not mention how to obtain mass or ratios of nuclides for the spent fuel. We just have generated simulation sample data as a fission product inventory using ORIGEN to construct spent fuel database. When unidentified radioactive materials are discovered, it is assumed that mass or ratios of the nuclides are analyzed by certain methods which are available in Korea as well as other countries. Then, a reactor type is determined, and finally, regression model can locate the origin or source of the radioactive materials by inferring operational histories using the fission product inventory. There are some studies on the classification of reactor type using statistical methods. As a part of predicting the operational history, this paper focuses on the development of spent fuel inference algorithm for nuclear forensics to estimate reactor operational histories that reflect characteristics of a reactor from the fission product inventory. The spent fuel inference algorithm for nuclear forensics has been developed on the basis of Kernel Regression (KR) which can conditionally generate the estimated output data by giving weight depending on geometric distance between sample data and input data. Also the KR which comprises Auto-Associative KR and Inferential KR have some benefits that can estimate only using data without assumption for specific physical model and re-learning process. Auto-associative KR (AAKR) was used to sort out outliers of input data and Inferential KR (IKR) was carried out to predict the operational histories for spent fuel. All of regression algorithms were performed in R programming tool. From the validation of both AAKR and IKR, it is expected that kernel regression has enough accuracy to be applied to nuclear forensics for spent fuel. Once the nuclear substance is measured to fill the fission product inventory with a certain reactor type, this kernel regression algorithm would predict the operational histories of spent fuel with high accuracy. Finally, this information is expected to be useful to find the origin of irradiated materials. (authors)

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Non-uniform active learning for Gaussian process models with applications to trajectory informed aerodynamic databases

The ability to non-uniformly weight the input space is desirable for many applications, and has been explored for space-filling approaches. Increased interests in linking models, such as in a digital twinning framework, increases the need for sampling emulators where they are most likely to be evaluated. In particular, here we apply non-uniform sampling methods for the construction of aerodynamic databases. This paper combines non-uniform weighting with active learning for Gaussian Processes (GPs) to develop a closed-form solution to a non-uniform active learning criterion. We accomplish this by utilizing a kernel density estimator as the weight function. We demonstrate the need and efficacy of this approach with an atmospheric entry example that accounts for both model uncertainty as well as the practical state space of the vehicle, as determined by forward modeling within the active learning loop.

42 ENGINEERING↗
DOE OSTI2017

Optimal smoothing length scale for actuator line models of wind turbine blades based on Gaussian body force distribution: Wind energy, actuator line model

The actuator line model (ALM) is a commonly used method to represent lifting surfaces such as wind turbine blades within large-eddy simulations (LES). In the ALM, the lift and drag forces are replaced by an imposed body force that is typically smoothed over several grid points using a Gaussian kernel with some prescribed smoothing width e. To date, the choice of e has most often been based on numerical considerations related to the grid spacing used in LES. However, especially for finely resolved LES with grid spacings on the order of or smaller than the chord length of the blade, the best choice of e is not known. In this work, a theoretical approach is followed to determine the most suitable value of e, based on an analytical solution to the linearized inviscid flow response to a Gaussian force. We find that the optimal smoothing width eopt is on the order of 14%-25% of the chord length of the blade, and the center of force is located at about 13%-26% downstream of the leading edge of the blade for the cases considered. These optimal values do not depend on angle of attack and depend only weakly on the type of lifting surface. It is then shown that an even more realistic velocity field can be induced by a 2-D elliptical Gaussian lift-force kernel. Some results are also provided regarding drag force representation.

17 WIND ENERGY↗
DOE OSTI2012

Experiences with High-Level Programming Directives for Porting Applications to GPUs. In: Facing the Multicore--Challenge II, Lecture Notes in Computer Science.

HPC systems now exploit GPUs within their compute nodes to accelerate program performance. As a result, high-end application development has become extremely complex at the node level. In addition to restructuring the node code to exploit the cores and specialized devices, the programmer may need to choose a programming model such as OpenMP or CPU threads in conjunction with an accelerator programming model to share and manage the different node resources. This comes at a time when programmer productivity and the ability to produce portable code has been recognized as a major concern. In order to offset the high development cost of creating CUDA or OpenCL kernels, directives have been proposed for programming accelerator devices, but their implications are not well known. In this paper, we evaluate the state of the art accelerator directives to program several applications kernels, explore transformations to achieve good performance, and examine the expressivity and performance penalty of using high-level directives versus CUDA. We also compare our results to OpenMP implementations to understand the benefits of running the kernels in the accelerator versus CPU cores.

CUDA Implementation↗