Search NASA⌕ Search

SEARCH · Search NASA

Results for “Kernel function”

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 289 records · Page 16
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

Global adjoint tomography: first-generation model

We present the first-generation global tomographic model constructed based on adjoint tomography, an iterative full-waveform inversion technique. Synthetic seismograms were calculated using GPU-accelerated spectral-element simulations of global seismic wave propagation, accommodating effects due to 3-D anelastic crust & mantle structure, topography & bathymetry, the ocean load, ellipticity, rotation, and self-gravitation. Fréchet derivatives were calculated in 3-D anelastic models based on an adjoint-state method. The simulations were performed on the Cray XK7 named ‘Titan’, a computer with 18 688 GPU accelerators housed at Oak Ridge National Laboratory. The transversely isotropic global model is the result of 15 tomographic iterations, which systematically reduced differences between observed and simulated three-component seismograms. Our starting model combined 3-D mantle model S362ANI with 3-D crustal model Crust2.0. We simultaneously inverted for structure in the crust and mantle, thereby eliminating the need for widely used ‘crustal corrections’. We used data from 253 earthquakes in the magnitude range 5.8 ≤ M w ≤ 7.0. We started inversions by combining ~30 s body-wave data with ~60 s surface-wave data. The shortest period of the surface waves was gradually decreased, and in the last three iterations we combined ~17 s body waves with ~45 s surface waves. We started using 180 min long seismograms after the 12th iteration and assimilated minor- and major-arc body and surface waves. The 15th iteration model features enhancements of well-known slabs, an enhanced image of the Samoa/Tahiti plume, as well as various other plumes and hotspots, such as Caroline, Galapagos, Yellowstone and Erebus. Furthermore, we see clear improvements in slab resolution along the Hellenic and Japan Arcs, as well as subduction along the East of Scotia Plate, which does not exist in the starting model. Point-spread function tests demonstrate that we are approaching the resolution of continental-scale studies in some areas, for example, underneath Yellowstone. Here, this is a consequence of our multiscale smoothing strategy in which we define our smoothing operator as a function of the approximate Hessian kernel, thereby smoothing gradients less wherever we have good ray coverage, such as underneath North America.

58 GEOSCIENCES↗
DOE OSTI2022

MemHC: An Optimized GPU Memory Management Framework for Accelerating Many-body Correlation

The many-body correlation function is a fundamental computation kernel in modern physics computing applications, e.g., Hadron Contractions in Lattice quantum chromodynamics (QCD). This kernel is both computation and memory intensive, involving a series of tensor contractions, and thus usually runs on accelerators like GPUs. Existing optimizations on many-body correlation mainly focus on individual tensor contractions (e.g., cuBLAS libraries and others). In contrast, this work discovers a new optimization dimension for many-body correlation by exploring the optimization opportunities among tensor contractions. More specifically, it targets general GPU architectures (both NVIDIA and AMD) and optimizes many-body correlation’s memory management by exploiting a set of memory allocation and communication redundancy elimination opportunities: first, GPU memory allocation redundancy: the intermediate output frequently occurs as input in the subsequent calculations; second, CPU-GPU communication redundancy: although all tensors are allocated on both CPU and GPU, many of them are used (and reused) on the GPU side only, and thus, many CPU/GPU communications (like that in existing Unified Memory designs) are unnecessary; third, GPU oversubscription: limited GPU memory size causes oversubscription issues, and existing memory management usually results in near-reuse data eviction, thus incurring extra CPU/GPU memory communications.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

Compadre Toolkit

The Compadre Toolkit solves a minimization problem, that once solved allows a user to reconstruct a function from sample data collected from a cloud of data sites. The solution to this minimization can also be used to assemble a linear system that can be solved numerically using some other software. The minimization problem the toolkit solves is to find the optimal finite polynomial in a user specified finite dimensional space with regards to minimizing an objective function which measures a linear sampling functional applied to the polynomial, compared to the same sampling functional applied to the source function, weighted by a positive kernel that has compact support.

SciDAC↗
DOE OSTI2022

LaplaceInterpolation.jl: A Julia package for fast interpolation on a grid

We implement a linear-time algorithm for interpolation on a regular multidimensional grid in the Julia language. The algorithm is an approximate Laplace interpolation (Press, 1992) when no parameters are given; and when parameters m ∈ Z and ϵ > 0 are set, the interpolant approximates a Matérn kernel, of which radial basis functions and polyharmonic splines are a special case. We implement, in addition, Neumann, Dirichlet (trivial), and average boundary conditions with potentially different aspect ratios in the different dimensions. The interpolant functions in arbitrary dimensions.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Advanced Finite-Volume Numerics and Source Term Assumptions for Kernel and G-Equation Modelling of Propane/Air Flames

Here G-Equation models represent propagating flame fronts with an implicit two-dimensional surface representation (level-set). Level-set methods are fast, as transport source terms for the implicit surface can be solved with finite-volume operators on the finite-volume domain, without having to build the actual surface. However, they include approximations whose practical effects are not properly understood. In this study, we improved the numerics of the FRESCO CFD code’s G-Equation solver and developed a new method to simulate kernel growth using signed distance functions and the analytical sphere-mesh overlap. We analyzed their role for simulating propane/air flames, using three well-established constant-volume configurations: a one-dimensional, freely propagating laminar flame; a disc-shaped, constant-volume swirl combustor; and torch-jet flame development through an orifice from a two-chamber device. We tested the explicit (sub-cycled) vs. implicit formulation for the standard transport operators (advection, diffusion, compressibility). In addition to the accurate flame swept-volume method for chemistry and species source term, we developed a more accurate estimator for the burnt/unburnt split cell composition. Then, we developed a signed-distance-function (SDF) based method which provides a more stable reinitialization of the level-set field at every time-step. We found that simplifying assumptions common to several G-Equation implementations, for straightforward terms such as compressibility and advection, lead to large errors in predicting the propagation of even laminar flames, with deviations up to ~300% in simulated vs. formulated flame speed. Conversely, the enhanced numerics enabled through the SDF field reinitialization and improved chemistry source term improve simulation stability and smooth flame propagation even with significantly larger solver time-steps.

42 ENGINEERING↗
NASA NTRS1992

Calculation of stress intensity factors in an isotropic multicracked plate: Part 2: Symbolic/numeric implementation

Analytical derivations of stress intensity factors (SIF's) of a multicracked plate can be complex and tedious. Recent advances, however, in intelligent application of symbolic computation can overcome these difficulties and provide the means to rigorously and efficiently analyze this class of problems. Here, the symbolic algorithm required to implement the methodology described in Part 1 is presented. The special problem-oriented symbolic functions to derive the fundamental kernels are described, and the associated automatically generated FORTRAN subroutines are given. As a result, a symbolic/FORTRAN package named SYMFRAC, capable of providing accurate SIF's at each crack tip, was developed and validated. Simple illustrative examples using SYMFRAC show the potential of the present approach for predicting the macrocrack propagation path due to existing microcracks in the vicinity of a macrocrack tip, when the influence of the microcrack's location, orientation, size, and interaction are taken into account.

Arnold, S. M.↗
NASA NTRS1997

A General Reversible Hereditary Constitutive Model: Theoretical Developments - Part 1

Using an internal-variable formalism as a starting point, we describe the viscoelastic extension of a previously-developed viscoplasticity formulation of the complete potential structure type. It is mainly motivated by experimental evidence for the presence of rate/time effects in the so-called quasilinear, reversible, material response range. Several possible generalizations are described, in the general format of hereditary-integral representations for non-equilibrium, stress-type, state variables, both for isotropic as well as anisotropic materials. In particular, thorough discussions are given on the important issues of thermodynamic admissibility requirements for such general descriptions, resulting in a set of explicit mathematical constraints on the associated kernel (relaxation and creep compliance) functions. In addition, a number of explicit, integrated forms are derived, under stress and strain control to facilitate the parametric and qualitative response characteristic studies reported here, as well as to help identify critical factors in the actual experimental characterizations from test data that will be reported in Part II.

Saleeb, A. F.↗
NASA NTRS2016

Using Ground Targets to Validate S-NPP VIIRS Day-Night Band Calibration

In this study, the observations from S-NPP VIIRS Day-Night band (DNB) and Moderate resolution bands (M bands) of Libya 4 and Dome C over the first four years of the mission are used to assess the DNB low gain calibration stability. The Sensor Data Records produced by NASA Land Product Evaluation and Algorithm Testing Element (PEATE) are acquired from nearly nadir overpasses for Libya 4 desert and Dome C snow surfaces. A kernel-driven bidirectional reflectance distribution function (BRDF) correction model is used for both Libya 4 and Dome C sites to correct the surface BRDF influence. At both sites, the simulated top-of-atmosphere (TOA) DNB reflectances based on SCIAMACHY spectral data are compared with Land PEATE TOA reflectances based on modulated Relative Spectral Response (RSR). In the Libya 4 site, the results indicate a decrease of 1.03% in Land PEATE TOA reflectance and a decrease of 1.01% in SCIAMACHY derived TOA reflectance over the period from April 2012 to January 2016. In the Dome C site, the decreases are 0.29% and 0.14%, respectively. The consistency between SCIAMACHY and Land PEATE data trends is good. The small difference between SCIAMACHY and Land PEATE derived TOA reflectances could be caused by changes in the surface targets, atmosphere status, and on-orbit calibration. The reflectances and radiances of Land PEATE DNB are also compared with matching M bands and the integral M bands based on M4, M5, and M7. The fitting trends of the DNB to integral M bands ratios indicate a 0.75% decrease at the Libya 4 site and a 1.89% decrease at the Dome C site. Part of the difference is due to an insufficient number of sampled bands available within the DNB wavelength range. The above results indicate that the Land PEATE VIIRS DNB product is accurate and stable. The methods used in this study can be used on other satellite instruments to provide quantitative assessments for calibration stability.

Chen, Xuexia↗
DOE OSTI2015

TMD factorization and evolution at large $b_T$

In using transverse-momentum-dependent (TMD) parton densities and fragmentation functions, important non-perturbative information is at large transverse position $b_T$. This concerns both the TMD functions and their evolution. Fits to high energy data tend to predict too rapid evolution when extrapolated to low energies where larger values of $b_T$ dominate. I summarize a new analysis of the issues. It results in a proposal for much weaker $b_T$ dependence at large $b_T$ for the evolution kernel, while preserving the accuracy of the existing fits. The results are particularly important for using transverse-spin-dependent functions like the Sivers function.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Exact Gaussian processes for massive datasets via non-stationary sparsity-discovering kernels

Abstract A Gaussian Process (GP) is a prominent mathematical framework for stochastic function approximation in science and engineering applications. Its success is largely attributed to the GP’s analytical tractability, robustness, and natural inclusion of uncertainty quantification. Unfortunately, the use of exact GPs is prohibitively expensive for large datasets due to their unfavorable numerical complexity of $$O(N^3)$$ O ( N 3 ) in computation and $$O(N^2)$$ O ( N 2 ) in storage. All existing methods addressing this issue utilize some form of approximation—usually considering subsets of the full dataset or finding representative pseudo-points that render the covariance matrix well-structured and sparse. These approximate methods can lead to inaccuracies in function approximations and often limit the user’s flexibility in designing expressive kernels. Instead of inducing sparsity via data-point geometry and structure, we propose to take advantage of naturally-occurring sparsity by allowing the kernel to discover—instead of induce—sparse structure. The premise of this paper is that the data sets and physical processes modeled by GPs often exhibit natural or implicit sparsities, but commonly-used kernels do not allow us to exploit such sparsity. The core concept of exact, and at the same time sparse GPs relies on kernel definitions that provide enough flexibility to learn and encode not only non-zero but also zero covariances. This principle of ultra-flexible, compactly-supported, and non-stationary kernels, combined with HPC and constrained optimization, lets us scale exact GPs well beyond 5 million data points.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2016

SU-F-SPS-06: Implementation of a Back-Projection Algorithm for 2D in Vivo Dosimetry with An EPID System

Purpose: To implement a back-projection algorithm for 2D dose reconstructions for in vivo dosimetry in radiation therapy using an Electronic Portal Imaging Device (EPID) based on amorphous silicon. Methods: An EPID system was used to calculate dose-response function, pixel sensitivity map, exponential scatter kernels and beam hardenig correction for the back-projection algorithm. All measurements were done with a 6 MV beam. A 2D dose reconstruction for an irradiated water phantom (30×30×30 cm{sup 3}) was done to verify the algorithm implementation. Gamma index evaluation between the 2D reconstructed dose and the calculated with a treatment planning system (TPS) was done. Results: A linear fit was found for the dose-response function. The pixel sensitivity map has a radial symmetry and was calculated with a profile of the pixel sensitivity variation. The parameters for the scatter kernels were determined only for a 6 MV beam. The primary dose was estimated applying the scatter kernel within EPID and scatter kernel within the patient. The beam hardening coefficient is σBH= 3.788×10{sup −4} cm{sup 2} and the effective linear attenuation coefficient is µAC= 0.06084 cm{sup −1}. The 95% of points evaluated had γ values not longer than the unity, with gamma criteria of ΔD = 3% and Δd = 3 mm, and within the 50% isodose surface. Conclusion: The use of EPID systems proved to be a fast tool for in vivo dosimetry, but the implementation is more complex that the elaborated for pre-treatment dose verification, therefore, a simplest method must be investigated. The accuracy of this method should be improved modifying the algorithm in order to compare lower isodose curves.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2025

Transverse-momentum-dependent pion structures from lattice QCD: Collins-Soper kernel, soft factor, TMDWF, and TMDPDF

We present the first lattice quantum chromodynamics (QCD) calculation of the pion valence-quark transverse-momentum-dependent parton distribution function (TMDPDF) within the framework of large-momentum effective theory (LaMET). Using correlators fixed in the Coulomb gauge (CG), we computed the quasi-TMD beam function for a pion with a mass of 300 MeV, a fine lattice spacing of 𝑎 =0.06 fm, and multiple large momenta up to 3 GeV. The intrinsic soft functions in the CG approach are extracted from form factors with large momentum transfer, and as a byproduct, we also obtain the corresponding Collins-Soper (CS) kernel. Our determinations of both the soft function and the CS kernel agree with perturbation theory at small transverse separations (𝑏 ⊥ ) between the quarks. At larger 𝑏 ⊥ , the CS kernel remains consistent with recent results obtained using both CG and gauge-invariant TMD correlators in the literature. By combining next-to-leading logarithmic factorization of the quasi-TMD beam function and the soft function, we obtain an 𝑥-dependent pion valence-quark TMDPDF for transverse separations 𝑏 ⊥ ≳1 fm. Interestingly, we find that the 𝑏 ⊥ dependence of the phenomenological parametrizations of TMDPDF for moderate values of 𝑥 are in reasonable agreement with our QCD determinations. In addition, we present results for the transverse-momentum-dependent wave function for a heavier pion with 670 MeV mass.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2015

SU-F-BRD-06: Robust Dose Calculation in Intensity Modulated Proton Therapy

Purpose: Commissioning data for intensity modulated proton therapy (IMPT) must be post-processed by fits to ad-hoc functions to derive the dose calculation kernel parameters in a treatment planning system (TPS). Whether from experimental measurement or Monte Carlo simulation, the limited and noisy nature of such data makes this task very challenging. We present a method to improve the modeling of the lateral dose distribution of clinical energy proton beams in water to commission an in-house IMPT dose calculation engine. Methods: A linear sum of three Gaussian distribution functions was fitted to the lateral dose data in logarithmic scale. Starting values of fitting solutions were determined from the Generalized Highland Approximation. We exhaustively optimized the combinations of data weights with upper bounds of the fitting solutions to minimize confidence intervals of the fitting solutions while maintaining the coefficient of determination (R{sup 2}). Results: Across all energies, average confidence bounds improved 72.88% [Max: 88.28%, Min: 55.05%] for small angle coulomb scattering, 114.25% [409.13%, 66.72%,] for nuclear scattering, and 68.66% [141.09%, 33.27%] for large angle coulomb scattering, while the coefficients of determination of the fits (R{sup 2}) remained comparable. On average R {sup 2} only changed 0.18% and were very close to 1 (approx. 0.999). Wilcoxon signed rank tests comparing unweighted/unbounded fits with weighted/bounded fits averaged 0.0146 (Max: 0.177, Min: 7.05×10−{sup 7}) for small angle Coulomb, 0.0903 (0.945, 7.05×10−{sup 7}) for nuclear, and 0.254 (0.871, 1.86×10−{sup 6}) for large angle Coulomb scattering. This allows rejection of the null hypothesis for small angle Coulomb scattering at the 0.015 level and nuclear interaction at the 0.1 level. Conclusion: Optimal weights assigned to IMPT lateral dose data minimized fitting to stochastic noise in the tail region. Optimizing the upper bounds of fitting parameters improved the robustness of fitting solutions. Both led to a more physically meaningful dose calculation kernel for the TPS. NIH/NCI K25CA168984, Eagles Cancer Research Career Development, The Lawrence W. and Marilyn W. Matteson Fund for Cancer Research, Mayo ASU Seed Grant, and The Kemper Marley Foundation.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2020

Hole–hole Tamm–Dancoff-approximated density functional theory: A highly efficient electronic structure method incorporating dynamic and static correlation

The study of photochemical reaction dynamics requires accurate as well as computationally efficient electronic structure methods for the ground and excited states. While time-dependent density functional theory (TDDFT) is not able to capture static correlation, complete active space self-consistent field methods neglect much of the dynamic correlation. Hence, inexpensive methods that encompass both static and dynamic electron correlation effects are of high interest. Here, we revisit hole–hole Tamm–Dancoff approximated (hh-TDA) density functional theory for this purpose. The hh-TDA method is the hole–hole counterpart to the more established particle–particle TDA (pp-TDA) method, both of which are derived from the particle–particle random phase approximation (pp-RPA). In hh-TDA, the N-electron electronic states are obtained through double annihilations starting from a doubly anionic (N+2 electron) reference state. In this way, hh-TDA treats ground and excited states on equal footing, thus allowing for conical intersections to be correctly described. Furthermore, the treatment of dynamic correlation is introduced through the use of commonly employed density functional approximations to the exchange-correlation potential. Additionally, we show that hh-TDA is a promising candidate to efficiently treat the photochemistry of organic and biochemical systems that involve several low-lying excited states—particularly those with both low-lying ππ* and nπ* states where inclusion of dynamic correlation is essential to describe the relative energetics. In contrast to the existing literature on pp-TDA and pp-RPA, we employ a functional-dependent choice for the response kernel in pp- and hh-TDA, which closely resembles the response kernels occurring in linear response and collinear spin-flip TDDFT.

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

Helicity evolution at small x: the single-logarithmic contribution

We calculate single-logarithmic corrections to the small-x flavor-singlet helicity evolution equations derived recently [1–3] in the double-logarithmic approximation. The new single-logarithmic part of the evolution kernel sums up powers of α s ln(1/x), which are an important correction to the dominant powers of α s ln 2 (1/x) summed up by the double-logarithmic kernel from [1–3] at small values of Bjorken x and with α s the strong coupling constant. The single-logarithmic terms arise separately from either the longitudinal or transverse momentum integrals. Consequently, the evolution equations we derive employing the light-cone perturbation theory simultaneously include the small-x evolution kernel and the leading-order polarized DGLAP splitting functions. We further enhance the equations by calculating the running coupling corrections to the kernel.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1987

An orthogonal oriented quadrature hexagonal image pyramid

An image pyramid has been developed with basis functions that are orthogonal, self-similar, and localized in space, spatial frequency, orientation, and phase. The pyramid operates on a hexagonal sample lattice. The set of seven basis functions consist of three even high-pass kernels, three odd high-pass kernels, and one low-pass kernel. The three even kernels are identified when rotated by 60 or 120 deg, and likewise for the odd. The seven basis functions occupy a point and a hexagon of six nearest neighbors on a hexagonal sample lattice. At the lowest level of the pyramid, the input lattice is the image sample lattice. At each higher level, the input lattice is provided by the low-pass coefficients computed at the previous level. At each level, the output is subsampled in such a way as to yield a new hexagonal lattice with a spacing sq rt 7 larger than the previous level, so that the number of coefficients is reduced by a factor of 7 at each level. The relationship between this image code and the processing architecture of the primate visual cortex is discussed.

Watson, Andrew B.↗
DOE OSTI2022

Advanced stationary and nonstationary kernel designs for domain-aware Gaussian processes

Gaussian process regression is a widely-applied method for function approximation and uncertainty quantification. The technique has gained popularity recently in the machine learning community due to its robustness and interpretability. The mathematical methods we discuss in this paper are an extension of the Gaussian-process framework. We are proposing advanced kernel designs that only allow for functions with certain desirable characteristics to be elements of the reproducing kernel Hilbert space (RKHS) that underlies all kernel methods and serves as the sample space for Gaussian process regression. These desirable characteristics reflect the underlying physics; two obvious examples are symmetry and periodicity constraints. In addition, non-stationary kernel designs can be defined in the same framework to yield flexible multi-task Gaussian processes. We will show the impact of advanced kernel designs on Gaussian processes using several synthetic and two scientific data sets. The results of our research show that including domain knowledge, communicated through advanced kernel designs, has a significant impact on the accuracy and relevance of the function approximation.

97 MATHEMATICS AND COMPUTING↗