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 577 records · Page 32
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

ConKer: An algorithm for evaluating correlations of arbitrary order

Context. High order correlations in the cosmic matter density have become increasingly valuable in cosmological analyses. However, computing these correlation functions is computationally expensive. Aims. We aim to circumvent these challenges by developing a new algorithm called ConKer for estimating correlation functions. Methods. This algorithm performs convolutions of matter distributions with spherical kernels using FFT. Since matter distributions and kernels are defined on a grid, it results in some loss of accuracy in the distance and angle definitions. We study the algorithm setting at which these limitations become critical and suggest ways to minimize them. Results. ConKer is applied to the CMASS sample of the SDSS DR12 galaxy survey and corresponding mock catalogs, and is used to compute the correlation functions up to correlation order n = 5. We compare the n = 2 and n = 3 cases to traditional algorithms to verify the accuracy of the new algorithm. We perform a timing study of the algorithm and find that three of the four distinct processes within the algorithm are nearly independent of the catalog size N , while one subdominant component scales as O ( N ). The dominant portion of the calculation has complexity of O ( N c 4/3 log N c ), where N c is the of cells in a three-dimensional grid corresponding to the matter density. Conclusions. We find ConKer to be a fast and accurate method of probing high order correlations in the cosmic matter density, then discuss its application to upcoming surveys of large-scale structure.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

Analytic Evolution of Singular Distribution Amplitudes in QCD

We describe a method of analytic evolution of distribution amplitudes (DA) that have singularities, such as non-zero values at the end-points of the support region, jumps at some points inside the support region and cusps. We illustrate the method by applying it to the evolution of a flat (constant) DA, anti-symmetric at DA and then use it for evolution of the two-photon generalized distribution amplitude. Our approach has advantages over the standard method of expansion in Gegenbauer polynomials, which requires infinite number of terms in order to accurately reproduce functions in the vicinity of singular points, and over a straightforward iteration of an initial distribution with evolution kernel. The latter produces logarithmically divergent terms at each iteration, while in our method the logarithmic singularities are summed from the start, which immediately produces a continuous curve, with only one or two iterations needed afterwards in order to get rather precise results.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2018

Unsupervised learning about 4D features of microparticle motion

Material clusters of different sizes are known to exist in high-temperature plasmas due to plasma-wall interactions. The facts that these clusters, ranging from sub-microns to above mm in size, can move from one location to another quickly and that there are a lot of them make high-speed imaging and tracking one of the best, effective, and sometimes only diagnostic. An unsupervised machine learning technique based on deconvolutional neural networks is developed to analyze two-camera videos of high-temperature microparticles generated from exploding wires. The neural network utilizes a locally competitive algorithm to infer representations and optimize a dictionary composed of kernels, or basis vectors, for image analysis. Our primary goal is to use this method for feature recognition and prediction of the time-dependent three-dimensional (or “4D”) microparticle motion. Features equivalent to local velocity vectors have been identified as the dictionary kernels or “building blocks” of the scene. The dictionary elements from the left and right camera views are found to be strongly correlated and satisfy the projection geometrical constraints. The results show that unsupervised machine learning techniques are promising approaches to process large sets of images for high-temperature plasmas and other scientific experiments. In conclusion, machine learning techniques can be useful to handle the large amount of data and therefore aid the understanding of plasma-wall interaction.

47 OTHER INSTRUMENTATION↗
DOE OSTI2021

A variational multiscale immersed meshfree method for heterogeneous materials

Abstract We introduce an immersed meshfree formulation for modeling heterogeneous materials with flexible non-body-fitted discretizations, approximations, and quadrature rules. The interfacial compatibility condition is imposed by a volumetric constraint, which avoids a tedious contour integral for complex material geometry. The proposed immersed approach is formulated under a variational multiscale based formulation, termed the variational multiscale immersed method (VMIM). Under this framework, the solution approximation on either the foreground or the background can be decoupled into coarse-scale and fine-scale in the variational equations, where the fine-scale approximation represents a correction to the residual of the coarse-scale equations. The resulting fine-scale solution leads to a residual-based stabilization in the VMIM discrete equations. The employment of reproducing kernel (RK) approximation for the coarse- and fine-scale variables allows arbitrary order of continuity in the approximation, which is particularly advantageous for modeling heterogeneous materials. The effectiveness of VMIM is demonstrated with several numerical examples, showing accuracy, stability, and discretization efficiency of the proposed method.

36 MATERIALS SCIENCE↗
DOE OSTI2013

Determination and verification of a 2D pencil-beam kernel for a radiosurgery system with cones

The quality and correctness of dosimetric data of small fields in stereotactic radiosurgery (SRS) depends significantly on the election of the detector employed in the measurements. This work provides an independent method of verification of these data through the determination of a polyenergetic 2-dimensional pencil-beam kernel for a BrainLAB SRS system with cones, employing the deconvolution/convolution of a reference experimental off-axis ratio (OAR) profile (cone diameter c{sub 0} = 35 mm). The kernel in real space k{sub c{sub 0}}(r,z{sub 0}) is convolved with the ideal fluence Φ for the cones 7.5 to 35 mm in diameter to obtain the OAR profiles, and the total scatter factors, S{sub t}, which are compared with experimental values of the same quantities. The experimental OARs and S{sub t} factors are measured in water with a PTW 60003 diamond detector. Additionally, the reference OAR is corrected for beam divergence and spectral fluence fluctuations defining a function of boundary correction factors (BF). The BF and Φ functions are transformed to the conjugate space with the zeroth-order Hankel transform, appropriated to the radial symmetry of the cones. Therefore, the kernel in real space k{sub c{sub 0}}(r,z{sub 0}) is the inverse Hankel transform of the ratio of the Hankel transforms of BF and Φ. Finally, an uncertainty analysis according to the Guide to the Expression of Uncertainty in Measurement is carried out for 3 different values of k{sub c{sub 0}}(r,z{sub 0}). Calculated and measured OARs agree within the dose/distance-to-agreement criteria of 2%/0.12 mm; while, S{sub t} factors agree within 2%. This procedure supplies an independent method to validate the dosimetric data necessary to feed treatment planning systems for SRS with cones.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2015

SU-E-T-745: The Use of Electron Dense Materials for Quality Assurance Testing in a Magnetic Field

Purpose: Due to the effect of the magnetic field on the dose deposition (skewing the beam), machine quality assurance (QA) tests for the MR-linac (MRI combined with a linear accelerator) need to be redesigned. In this work we focus on the redesign of QA tests that address geometrical accuracy of the system. Methods: Using electron dense materials (e.g. copper in our experiment) the dose kernel is minimized and thereby the effect of the magnetic field on the dose distribution. This approach is supported by Monte-Carlo simulations and can be used in practice with film measurements. Two examples of QA tests are presented: beam profile and star-shot measurements. Results: The novel method was verified by performing both measurements on a conventional linac and the MR-linac with a film that was sandwiched between copper layers. Measurements were compared with a reference setup which was similar to setup used in clinical practice. On a conventional linac the experimental outcome showed good agreement between the reference and the new setup for both QA tests. The results from the MR-linac showed that the symmetry of the beam profile was restored in presence of the copper layers in the setup and that the isocenter size can be determined accurately with the introduced star-shot setup (see supporting material). Conclusion: The use of electron dense materials for QA tests was shown to be a simple and effective method to remove the effects on the dose distribution enabling assessment of geometrical accuracy of a MR-linac system. The use of high dense materials is not limited to the presented QA tests only, but has a broad applicability for beam specific QA tests in presence of a magnetic field.

43 PARTICLE ACCELERATORS↗
DOE OSTI2018

Nonlinear Integral Equations with Potential-Type Kernels on a Segment

We study various classes of nonlinear equations containing operators of potential type (Riesz potential). By the method of monotone operators in the Lebesgue spaces of real-valued functions L{sub p}(a, b) we prove global theorems on the existence, uniqueness, estimates, and methods of construction of their solutions. We present applications that illustrate the results obtained.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2016

Application of mixed kernels function (MKF) based support vector regression model (SVR) for CO 2 – Reservoir oil minimum miscibility pressure prediction

Carbon dioxide (CO 2 ) injection into oil reservoirs is considered a mature enhanced oil recovery (EOR) technique for conventional reservoirs. The local displacement efficiency of the CO 2 -EOR process is highly dependent on the minimum miscibility pressure (MMP), estimating this parameter is critical to design of the CO 2 injection process. Traditional empirical methods to test the CO 2 -oil MMP are time consuming and expensive; derived correlations are fast but not accurate. Therefore, an efficient and reliable method to determine MMP is beneficial. Here in this study, a mixed kernels function (MKF) based support vector regression (SVR) model was developed and used to predict the MMP for both pure and impure CO 2 injection cases. Four parameters were chosen as input parameters: (1) reservoir temperature; (2) average critical temperature; (3) molecular weight of pentane plus (C 5+ ) fraction of crude oil, and; (4) the ratio of volatile components to intermediate components in crude oil. MMP was selected as the desired output parameter to train and test this newly developed model. The performance of basic kernels function based SVR model is compared with that of this newly developed MKF-SVR model. The well-trained MFK-SVR was compared with three well-established published correlations, demonstrated the highest correlation coefficient (R of 0.9381), lowest root mean square error (RMSE of 1.9151), smallest average absolute error (AAE of 1.1406) and maximum absolute error (MAE of 4.6291). We believe that the proposed MFK-SVM model is a more reliable and stable regression model to predict MMP. In addition, a sensitivity analysis was conducted to evaluate the physical correctness. It indicates that the predicted results from the newly developed model are in excellent agreement with previous empirical work.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2022

Regularized inversion of aerosol hygroscopic growth factor probability density function: application to humidity-controlled fast integrated mobility spectrometer measurements

Abstract. Aerosol hygroscopic growth plays an important role in atmospheric particle chemistry and the effects of aerosol on radiation and hence climate. The hygroscopic growth is often characterized by a growth factor probability density function (GF-PDF), where the growth factor is defined as the ratio of the particle size at a specified relative humidity to its dry size. Parametric, least-squares methods are the most widely used algorithms for inverting the GF-PDF from measurements of the humidified tandem differential mobility analyzer (HTDMA) and have been recently applied to the GF-PDF inversion from measurements of the humidity-controlled fast integrated mobility spectrometer (HFIMS). However, these least-squares methods suffer from noise amplification due to the lack of regularization in solving the ill-posed problem, resulting in significant fluctuations in the retrieved GF-PDF and even occasional failures of convergence. In this study, we introduce nonparametric, regularized methods to invert the aerosol GF-PDF and apply them to HFIMS measurements. Based on the HFIMS kernel function, the forward convolution is transformed into a matrix-based form, which facilitates the application of the nonparametric inversion methods with regularizations, including Tikhonov regularization and Twomey's iterative regularization. Inversions of the GF-PDF using the nonparameteric methods with regularization are demonstrated using HFIMS measurements simulated from representative GF-PDFs of ambient aerosols. The characteristics of reconstructed GF-PDFs resulting from different inversion methods, including previously developed least-squares methods, are quantitatively compared. The result shows that Twomey's method generally outperforms other inversion methods. The capabilities of Twomey's method in reconstructing the pre-defined GF-PDFs and recovering the mode parameters are validated.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2017

Applications of Anomaly Detection and Precursor Identification in Airspace Operations

As we continue to advance the U.S. National Airspace into the next generation of air traffic, we face challenges in both increase in complexity, as well as, a significant growth in traffic volume. Addressing these challenges, while maintaining the same level of safety is an important application of data mining. Because of these significant shifts in airspace design and usage there is a need to identify current and emergent safety risks along with their potential precursors. In recent years NASA has made advancements in developing scalable methods to address this effort in the Big Data paradigm. Multiple kernel anomaly detection approaches have been employed on both surveillance radar data and flight operational quality assurance data to identify operationally significant safety risks. Additionally, events have been explored with a recently developed precursor identification tool to discover states that reveal an increased probability of a safety event. These tools can be used to discover emerging safety risks that may not be currently monitored, which allows for mitigation tactics to be employed and ultimately make the overall airspace safer. This talk will discuss an overview of these methods and a discussion of the findings.

anomaly detection↗
DOE OSTI2012

Nonparametric Method for Genomics-Based Prediction of Performance of Quantitative Traits Involving Epistasis in Plant Breeding

Genomic selection (GS) procedures have proven useful in estimating breeding value and predicting phenotype with genome-wide molecular marker information. However, issues of high dimensionality, multicollinearity, and the inability to deal effectively with epistasis can jeopardize accuracy and predictive ability. We, therefore, propose a new nonparametric method, pRKHS, which combines the features of supervised principal component analysis (SPCA) and reproducing kernel Hilbert spaces (RKHS) regression, with versions for traits with no/low epistasis, pRKHS-NE, to high epistasis, pRKHS-E. Instead of assigning a specific relationship to represent the underlying epistasis, the method maps genotype to phenotype in a nonparametric way, thus requiring fewer genetic assumptions. SPCA decreases the number of markers needed for prediction by filtering out low-signal markers with the optimal marker set determined by cross-validation. Principal components are computed from reduced marker matrix (called supervised principal components, SPC) and included in the smoothing spline ANOVA model as independent variables to fit the data. The new method was evaluated in comparison with current popular methods for practicing GS, specifically RR-BLUP, BayesA, BayesB, as well as a newer method by Crossa et al., RKHS-M, using both simulated and real data. Results demonstrate that pRKHS generally delivers greater predictive ability, particularly when epistasis impacts trait expression. Beyond prediction, the new method also facilitates inferences about the extent to which epistasis influences trait expression.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2014

Software optimization for electrical conductivity imaging in polycrystalline diamond cutters

We previously reported on an electrical conductivity imaging instrument developed for measurements on polycrystalline diamond cutters. These cylindrical cutters for oil and gas drilling feature a thick polycrystalline diamond layer on a tungsten carbide substrate. The instrument uses electrical impedance tomography to profile the conductivity in the diamond table. Conductivity images must be acquired quickly, on the order of 5 sec per cutter, to be useful in the manufacturing process. This paper reports on successful efforts to optimize the conductivity reconstruction routine, porting major portions of it to NVIDIA GPUs, including a custom CUDA kernel for Jacobian computation.

36 MATERIALS SCIENCE↗
DOE OSTI2019

FVCOM-plume – A Three-dimensional Lagrangian Outfall Plume Dilution and Transport Model for Dynamic Tidal Environments: Model Development

Although numerous steady-state effluent plume dilution models are in use, their application in tidal environments remains a challenge. Three-dimensional dynamic circulation models are also inadequate, often due to the lack of required resolution and simplifying assumptions. To overcome these limitations, we present FVCOM-plume—an outfall plume dilution and transport model that operates within the Finite Volume Coastal Ocean Model (FVCOM) framework. It provides simultaneous inclusion of near-field dilution and far-field plume transport processes. The near-field is based on UM3 model using Lagrangian Control Volume approach to compute buoyant plume trajectory and dilution from multiport diffusers. The far-field uses neutrally buoyant particles with point masses and the random walk method to solve unsteady advection-diffusion processes. A density kernel approach is used to compute concentrations at point locations and analyze far-field plume characteristics. The results demonstrate the ability of FVCOM-plume to simultaneously capture near-field and far-field effluent plume dynamics in tidal environments.

Plume Dilution Models↗
DOE OSTI2021

Improving the astrometric solution of the Hyper Suprime-Cam with anisotropic Gaussian processes

Context. We study astrometric residuals from a simultaneous fit of Hyper Suprime-Cam images. Aims. We aim to characterize these residuals and study the extent to which they are dominated by atmospheric contributions for bright sources. Methods. We used Gaussian process interpolation with a correlation function (kernel) measured from the data to smooth and correct the observed astrometric residual field. Results. We find that a Gaussian process interpolation with a von Kármán kernel allows us to reduce the covariances of astrometric residuals for nearby sources by about one order of magnitude, from 30 mas 2 to 3 mas 2 at angular scales of ~1 arcmin. This also allows us to halve the rms residuals. Those reductions using Gaussian process interpolation are similar to recent result published with the Dark Energy Survey dataset. We are then able to detect the small static astrometric residuals due to the Hyper Suprime-Cam sensors effects. We discuss how the Gaussian process interpolation of astrometric residuals impacts galaxy shape measurements, particularly in the context of cosmic shear analyses at the Rubin Observatory Legacy Survey of Space and Time.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

Incorporation of memory effects in coarse-grained modeling via the Mori-Zwanzig formalism

The Mori-Zwanzig formalism for coarse-graining a complex dynamical system typically introduces memory effects. The Markovian assumption of delta-correlated fluctuating forces is often employed to simplify the formulation of coarse-grained (CG) models and numerical implementations. However, when the time scales of a system are not clearly separated, the memory effects become strong and the Markovian assumption becomes inaccurate. To this end, we incorporate memory effects into CG modeling by preserving non-Markovian interactions between CG variables, and the memory kernel is evaluated directly from microscopic dynamics. For a specific example, molecular dynamics (MD) simulations of star polymer melts are performed while the corresponding CG system is defined by grouping many bonded atoms into single clusters. Then, the effective interactions between CG clusters as well as the memory kernel are obtained from the MD simulations. The constructed CG force field with a memory kernel leads to a non-Markovian dissipative particle dynamics (NM-DPD). Quantitative comparisons between the CG models with Markovian and non-Markovian approximations indicate that including the memory effects using NM-DPD yields similar results as the Markovian-based DPD if the system has clear time scale separation. However, for systems with small separation of time scales, NM-DPD can reproduce correct short-time properties that are related to how the system responds to high-frequency disturbances, which cannot be captured by the Markovian-based DPD model.

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

Insights into the spurious long-range nature of local rs-dependent non-local exchange-correlation kernels

A systematic route to go beyond the exact exchange plus random phase approximation (RPA) is to include a physical exchange-correlation kernel in the adiabatic-connection fluctuation-dissipation theorem. Previously, [D. Lu, J. Chem. Phys. 140, 18A520 (2014)], we found that non-local kernels with a screening length depending on the local Wigner-Seitz radius, r s (r), suffer an error associated with a spurious long-range repulsion in van der Waals bounded systems, which deteriorates the binding energy curve as compared to RPA. Here, we analyze the source of the error and propose to replace r s (r) by a global, average r s in the kernel. Exemplary studies with the Corradini, del Sole, Onida, and Palummo kernel show that while this change does not affect the already outstanding performance in crystalline solids, using an average r s significantly reduces the spurious long-range tail in the exchange-correlation kernel in van der Waals bounded systems. Finally, when this method is combined with further corrections using local dielectric response theory, the binding energy of the Kr dimer is improved three times as compared to RPA.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Risk Reduction of CO 2 Storage with Stochastic Simulations

The purpose of this project was to create an efficient and flexible tool based on principal component analysis (PCA) for a generation of differentiable realizations of porosity and permeability fields. Of particular interest were data which had significant connectivity between patches of the same rock type, which is called “binary”. The efficiency requirement needed a so called kernel PCA, and the binary images were found with the method called optimization based PCA (OPCA) for which interpretation is provided. The ability to honor given data at particular locations is also incorporated. The tool worked very well and efficiently for two dimensional (2-D) simulations. Furthermore, the tool was applied to a three-dimensional (3-D) dataset and determined that conditioning can be used to maintain the connectivity between vertical layers.

47 OTHER INSTRUMENTATION↗
NASA NTRS1995

Algorithms for parallel flow solvers on message passing architectures

The purpose of this project has been to identify and test suitable technologies for implementation of fluid flow solvers -- possibly coupled with structures and heat equation solvers -- on MIMD parallel computers. In the course of this investigation much attention has been paid to efficient domain decomposition strategies for ADI-type algorithms. Multi-partitioning derives its efficiency from the assignment of several blocks of grid points to each processor in the parallel computer. A coarse-grain parallelism is obtained, and a near-perfect load balance results. In uni-partitioning every processor receives responsibility for exactly one block of grid points instead of several. This necessitates fine-grain pipelined program execution in order to obtain a reasonable load balance. Although fine-grain parallelism is less desirable on many systems, especially high-latency networks of workstations, uni-partition methods are still in wide use in production codes for flow problems. Consequently, it remains important to achieve good efficiency with this technique that has essentially been superseded by multi-partitioning for parallel ADI-type algorithms. Another reason for the concentration on improving the performance of pipeline methods is their applicability in other types of flow solver kernels with stronger implied data dependence. Analytical expressions can be derived for the size of the dynamic load imbalance incurred in traditional pipelines. From these it can be determined what is the optimal first-processor retardation that leads to the shortest total completion time for the pipeline process. Theoretical predictions of pipeline performance with and without optimization match experimental observations on the iPSC/860 very well. Analysis of pipeline performance also highlights the effect of uncareful grid partitioning in flow solvers that employ pipeline algorithms. If grid blocks at boundaries are not at least as large in the wall-normal direction as those immediately adjacent to them, then the first processor in the pipeline will receive a computational load that is less than that of subsequent processors, magnifying the pipeline slowdown effect. Extra compensation is needed for grid boundary effects, even if all grid blocks are equally sized.

Vanderwijngaart, Rob F.↗