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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 649 records · Page 36

All sources NASA NTRS DOE OSTI

DOE OSTI2020

Program Verification for Extreme-Scale Applications (Final Scientific/Technical Report)

The project seeks to develop tools and techniques to help software developers verify the correctness and accuracy of their code. The target domain is scientific software of the kind widely used and developed in the Department of Energy research community, with a particular focus on "extreme scale" programs - those that are expected to involve possibly millions of parallel threads of execution. The report covers the University of Delaware contribution to the collaborative project. The project had a number of successful outcomes, especially regarding the development and extension of the CIVL software verification framework. CIVL is a verification tool for C or Fortran programs that use MPI, OpenMP, CUDA, and/or Pthreads for parallelization. CIVL went through vast improvements and extensions, and was successfully applied to a number of challenging codes. It found a subtle bug in the Devito PDE framework. It was able to verify the functional correctness of a conjugate gradient solver using a novel probabilistic technique with vanishingly small chance of error. CIVL was used very successfully in verification competitions, and the CIVL solutions to the competition challenges were presented and published. CIVL was also used successfully by other researchers on a computational chemistry kernel.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1994

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

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

Hu, Fang Q.↗

DOE OSTI2024

Efficient Parallelization of Irregular Applications on GPU Architectures

With the enlarging computation capacity of general Graphics Processing Units (GPUs), leveraging GPUs to accelerate parallel applications has become a critical topic in academia and industry. However, a wide range of irregular applications with the computation-/memory-intensive nature cannot easily achieve high GPU utilization. The challenges mainly involve the following aspects: first, data dependence leads to coarse-grained kernel and inefficient parallelism; second, heavy GPU memory usage may cause frequent memory evictions and extra overhead of I/O; third, specific computation patterns produce memory redundancies; last, workload balance and data reusability conjunctly benefit the overall performance, but there may exist a dynamic trade-off between them. Targeting these challenges, this dissertation proposes multiple optimizations to accelerate two real-world applications: many-body correlation functions to simulate nuclear physics in a large-scale scientific system; the other is the eALS-based matrix factorization recommendation system. To accelerate the calculations of many-body correlation functions, this dissertation presents three frameworks in GPU memory management and multi-GPU scheduling. Firstly, an optimized systematic GPU memory management framework, MemHC, utilizes a series of new memory reduction designs in GPU memory allocation, CPU/GPU communications, and GPU memory oversubscription. Secondly, an enhanced multi-GPU scheduling framework, MICCO, particularly by taking both data dimension (e.g., data reuse and data eviction) and computation dimension into account. MICCO designs a heuristic scheduling algorithm and a machine learning-based regression model to generate the optimal settings of a proposed new concept to manage the trade-off. Thirdly, a locality-aware multi-GPU scheduling framework. This scheduler leverages pipeline batch generation with a looking-ahead strategy by building local dependency graphs for memory transfer reduction and better data reuse, achieving up to 79.92% memory cost reduction and 1.67x speedup. To parallelize the eALS-based recommendation system, this dissertation proposes an efficient CPU/GPU heterogeneous recommendation system, HEALS. HEALS employs newly designed architecture-adaptive data formats to achieve load balance and good data locality on CPU and GPU. To mitigate the data dependence, HEALS presents a CPU/GPU collaboration model for both task parallelism and data parallelism with multiple kernel computation optimizations. In summary, this dissertation efficiently accelerates two typical irregular applications on GPUs by building four frameworks, including CPU/GPU collaboration, GPU memory management, and multi-GPU scheduling.

Wang, Qihan↗

DOE OSTI2020

Density-Matrix Based Extended Lagrangian Born–Oppenheimer Molecular Dynamics

Extended Lagrangian Born–Oppenheimer molecular dynamics [ Phys. Rev. Lett. 2008, 100, 123004] is presented for Hartree–Fock theory, where the extended electronic degrees of freedom are represented by a density matrix, including fractional occupation numbers at elevated electronic temperatures. In contrast to regular direct Born–Oppenheimer molecular dynamics simulations, no iterative self-consistent field optimization is required prior to the force evaluations. To sample regions of the potential energy landscape where the gap is small or vanishing, which leads to particular convergence problems in regular direct Born–Oppenheimer molecular dynamics simulations, an adaptive integration scheme for the extended electronic degrees of freedom is presented. The integration scheme is based on a tunable, low-rank approximation of a fourth-order kernel, K, that determines the metric tensor, T ≡ K T K, used in the extended harmonic oscillator of the Lagrangian that generates the dynamics of the electronic degrees of freedom. Here, the formulation and algorithms provide a general guide to implement extended Lagrangian Born–Oppenheimer molecular dynamics for quantum chemistry, density functional theory, and semiempirical methods using a density matrix formalism.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

Filtering, averaging, and scale dependency in homogeneous variable density turbulence

We investigate relationships between statistics obtained from filtering and from ensemble or Reynolds-averaging turbulence flow fields as a function of length scale. Generalized central moments in the filtering approach are expressed as inner products of generalized fluctuating quantities, q ' ( ξ , x ) = q ( ξ ) - q ¯ ( x ) , representing fluctuations of a field q ( ξ ) , at any point ξ, with respect to its filtered value at x. For positive-definite filter kernels, these expressions provide a scale-resolving framework, with statistics and realizability conditions at any length scale. In the small-scale limit, scale-resolving statistics become zero. In the large-scale limit, scale-resolving statistics and realizability conditions are the same as in the Reynolds-averaged description. Using direct numerical simulations (DNS) of homogeneous variable density turbulence, we diagnose Reynolds stresses, T i j , resolved kinetic energy, kr, turbulent mass-flux velocity, a i , and density-specific volume covariance, b, defined in the scale-resolving framework. These variables, and terms in their governing equations, vary smoothly between zero and their Reynolds-averaged definitions at the small and large scale limits, respectively. At intermediate scales, the governing equations exhibit interactions between terms that are not active in the Reynolds-averaged limit. For example, in the Reynolds-averaged limit, b follows a decaying process driven by a destruction term; at intermediate length scales, it is a balance between production, redistribution, destruction, and transport, where b grows as the density spectrum develops and then decays when mixing becomes strong enough. This work supports the notion of a generalized, length-scale adaptive model that converges to DNS at high resolutions and to Reynolds-averaged statistics at coarse resolutions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2020

Custom NEON AOP reflectance mosaics and maps of shade masks, canopy water content

This mosaic of visible to shortwave infrared (VSWIR) data was derived from the assignable asset NEON AOP radiance data that was collected by LBNL’s Watershed Function SFA during the summer of 2018 (DOI: 10.15485/1617204). This atmospheric correction was completed to take into account site-specific terrain variability in the 334 km2 survey area centered around Crested Butte, CO. The atmospheric correction was completed using ACORN atmospheric correction software executed on 200 x 200 pixel kernels rather than line by line in order to capture local flight altitude conditions. Manual cloud delineation removed any small shaded areas that occurred within the data collection areas. Mosaics were developed based on preference for days that were in close timing proximity to the coincident ground campaign and with mosaicing criteria that minimized the angle between the sun and the sensor to retrieve the most consistent reflectance (min_phase_refl_tiled). Finally, we generated shade masks based on the geometry between the sun angle at time of flight, the ground surface, and the sensor. Here we provide the orthomosaiced estimated reflectance data (min_phase_refl_tiled.tif), canopy water content (wtrl), estimated atmospheric water vapor (wtrv), the observational data (obs) for this particular mosaic, the estimated visibility (vis), shade masks derived from the digital surface elevation model (shade) and the digital terrain model (shade_tch) added to the canopy height model (DOI: 10.15485/1617203), and a wavelength metadata file. Each image file is provided as a GeoTiff with internal tiling and LZW compression. These data can also be found on Google Earth Engine for extraction, download, and analysis of smaller extents (Google Earth Engine data is linked in external link table).

2018 NEON and 2025 CHESS Campaigns↗

DOE OSTI2019

The Borg Cube Simulation: Cosmological Hydrodynamics with CRK-SPH

A challenging requirement posed by next-generation observations is a firm theoretical grasp of the impact of baryons on structure formation. Cosmological hydrodynamic simulations modeling gas physics are vital in this regard. A high degree of modeling flexibility exists in this space, making it important to explore a range of methods in order to gauge the accuracy of simulation predictions. We present results from the first cosmological simulation using Conservative Reproducing Kernel Smoothed Particle Hydrodynamics (CRK-SPH). We employ two simulations: one evolved purely under gravity, and the other with nonradiative hydrodynamics. Each contains 2 × 2304 3 cold dark matter plus baryon particles in an 800 h -1 Mpc box. We compare statistics to previous nonradiative simulations including power spectra, mass functions, baryon fractions, and concentration. We find self-similar radial profiles of gas temperature, entropy, and pressure and show that a simple analytic model recovers these results to better than 40% over two orders of magnitude in mass. We quantify the level of nonthermal pressure support in halos and demonstrate that hydrostatic mass estimates are biased low by 24% (10%) for halos of mass 10 15 (10 13 ) h -1 M⊙. Here, we compute angular power spectra for the thermal and kinematic Sunyaev–Zel’dovich effects and find good agreement with the low-ℓ Planck measurements. Finally, artificial scattering between particles of unequal mass is shown to have a large impact on the gravity-only run, and we highlight the importance of better understanding this issue in hydrodynamic applications. This is the first in a simulation campaign using CRK-SPH, with future work including subresolution gas treatments.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRSDate not supplied

Assessing Potential Geophysical and Environmental Impacts from Frequent Rocket Launch Missions at Kennedy Space Center

Kennedy Space Center (KSC) in Florida has been utilized for space missions over many years with a gradually increasing number of rocket launches. There have been multiple studies where air-coupled acoustic waves and infrasound originating from launched rockets were used for operational purposes, such as locating booster trajectories as a function of changing atmospheric conditions. However, no study has utilized the acoustoelastic waves as a signal source for subsurface seismic investigations. We conducted a dispersion analysis using the seismic energy recorded from the Artemis I rocket launch at KSC in November 2022, and from these results we generated depth-sensitivity kernels at different wave frequencies. The kernels were compared with sedimentary core data to verify boundaries of carbonate layers above the Floridan Aquifer System. Accumulative information of bedrock-sediment boundary across the sedimentary platform could especially provide geo-structural evidence that manifests the configuration of coastal features. Continuous dispersion analysis of seismic recordings from consecutive rocket launches also has potential to identify non-stationary environmental effects, such as reorientation of sedimentary structures and fluctuation of the saltwater/groundwater lens from gravitational tides, which may affect erosional susceptibility of coastal features.

Han Byul Woo↗

NASA NTRSDate not supplied

Assessing Potential Geophysical and Environmental Impacts from Frequent Rocket Launch Missions at Kennedy Space Center

Rocket Launch Seismicity↗

DOE OSTI2016

The Suite for Embedded Applications and Kernels

Many applications of high performance embedded computing are limited by performance or power bottlenecks. We have designed SEAK, a new benchmark suite, (a) to capture these bottlenecks in a way that encourages creative solutions to these bottlenecks? and (b) to facilitate rigorous, objective, end-user evaluation for their solutions. To avoid biasing solutions toward existing algorithms, SEAK benchmarks use a mission-centric (abstracted from a particular algorithm) and goal-oriented (functional) specification. To encourage solutions that are any combination of software or hardware, we use an end-user blackbox evaluation that can capture tradeoffs between performance, power, accuracy, size, and weight. The tradeoffs are especially informative for procurement decisions. We call our benchmarks future proof because each mission-centric interface and evaluation remains useful despite shifting algorithmic preferences. It is challenging to create both concise and precise goal-oriented specifications for mission-centric problems. This paper describes the SEAK benchmark suite and presents an evaluation of sample solutions that highlights power and performance tradeoffs.

Source record↗

DOE OSTI2024

Vacancy-Induced Tunable Kondo Effect in Twisted Bilayer Graphene

In single sheets of graphene, vacancy-induced states have been shown to host an effective spin-1/2 hole that can be Kondo screened at low temperatures. Here, we show how these vacancy-induced impurity states survive in twisted bilayer graphene (TBG), which thus provides a tunable system to probe the critical destruction of the Kondo effect in pseudogap hosts. Ab initio calculations and atomic-scale modeling are used to determine the nature of the vacancy states in the vicinity of the magic angle in TBG, demonstrating that the vacancy can be treated as a quantum impurity. Utilizing this insight, we construct an Anderson impurity model with a TBG host that we solve using the numerical renormalization group combined with the kernel polynomial method. We determine the phase diagram of the model and show how there is a strict dichotomy between vacancies in the AA/BB versus AB/BA tunneling regions. In AB/BA vacancies, the Kondo temperature at the magic angle develops a broad distribution with a tail to vanishing temperatures due to multifractal wave functions at the magic angle. Finally, we argue that scanning tunneling microscopy in the vicinity of the vacancy can act as a probe of both the critical single-particle states and the underlying many-body ground state in magic-angle TBG.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2016

TU-H-BRC-05: Stereotactic Radiosurgery Optimized with Orthovoltage Beams

Purpose: To achieve improved stereotactic radiosurgery (SRS) dose distributions using orthovoltage energy fluence modulation with inverse planning optimization techniques. Methods: A pencil beam model was used to calculate dose distributions from the institution’s orthovoltage unit at 250 kVp. Kernels for the model were derived using Monte Carlo methods as well as measurements with radiochromic film. The orthovoltage photon spectra, modulated by varying thicknesses of attenuating material, were approximated using open-source software. A genetic algorithm search heuristic routine was used to optimize added tungsten filtration thicknesses to approach rectangular function dose distributions at depth. Optimizations were performed for depths of 2.5, 5.0, and 7.5 cm, with cone sizes of 8, 10, and 12 mm. Results: Circularly-symmetric tungsten filters were designed based on the results of the optimization, to modulate the orthovoltage beam across the aperture of an SRS cone collimator. For each depth and cone size combination examined, the beam flatness and 80–20% and 90–10% penumbrae were calculated for both standard, open cone-collimated beams as well as for the optimized, filtered beams. For all configurations tested, the modulated beams were able to achieve improved penumbra widths and flatness statistics at depth, with flatness improving between 33 and 52%, and penumbrae improving between 18 and 25% for the modulated beams compared to the unmodulated beams. Conclusion: A methodology has been described that may be used to optimize the spatial distribution of added filtration material in an orthovoltage SRS beam to result in dose distributions at depth with improved flatness and penumbrae compared to standard open cones. This work provides the mathematical foundation for a novel, orthovoltage energy fluence-modulated SRS system.

60 APPLIED LIFE SCIENCES↗

NASA NTRS1987

Approximating the linear quadratic optimal control law for hereditary systems with delays in the control

The fundamental control synthesis issue of establishing a priori convergence rates of approximation schemes for feedback controllers for a class of distributed parameter systems is addressed within the context of hereditary systems. Specifically, a factorization approach is presented for deriving approximations to the optimal feedback gains for the linear regulator-quadratic cost problem associated with time-varying functional differential equations with control delays. The approach is based on a discretization of the state penalty which leads to a simple structure for the feedback control law. General properties of the Volterra factors of Hilbert-Schmidt operators are then used to obtain convergence results for the controls, trajectories and feedback kernels. Two algorithms are derived from the basic approximation scheme, including a fast algorithm, in the time-invariant case. A numerical example is also considered.

Milman, Mark H.↗

NASA NTRS1988

Approximating the linear quadratic optimal control law for hereditary systems with delays in the control

The fundamental control synthesis issue of establishing a priori convergence rates of approximation schemes for feedback controllers for a class of distributed parameter systems is addressed within the context of hereditary schemes. Specifically, a factorization approach is presented for deriving approximations to the optimal feedback gains for the linear regulator-quadratic cost problem associated with time-varying functional differential equations with control delays. The approach is based on a discretization of the state penalty which leads to a simple structure for the feedback control law. General properties of the Volterra factors of Hilbert-Schmidt operators are then used to obtain convergence results for the controls, trajectories and feedback kernels. Two algorithms are derived from the basic approximation scheme, including a fast algorithm, in the time-invariant case. A numerical example is also considered.

Milman, Mark H.↗

NASA NTRS1993

A dual reciprocal boundary element formulation for viscous flows

The advantages inherent in the boundary element method (BEM) for potential flows are exploited to solve viscous flow problems. The trick is the introduction of a so-called dual reciprocal technique in which the convective terms are represented by a global function whose unknown coefficients are determined by collocation. The approach, which is necessarily iterative, converts the governing partial differential equations into integral equations via the distribution of fictitious sources or dipoles of unknown strength on the boundary. These integral equations consist of two parts. The first is a boundary integral term, whose kernel is the unknown strength of the fictitious sources and the fundamental solution of a convection-free flow problem. The second part is a domain integral term whose kernel is the convective portion of the governing PDEs. The domain integration can be transformed to the boundary by using the dual reciprocal (DR) concept. The resulting formulation is a pure boundary integral computational process.

Lafe, Olu↗

NASA NTRS1974

Nonparametric maximum likelihood estimation of probability densities by penalty function methods

When it is known a priori exactly to which finite dimensional manifold the probability density function gives rise to a set of samples, the parametric maximum likelihood estimation procedure leads to poor estimates and is unstable; while the nonparametric maximum likelihood procedure is undefined. A very general theory of maximum penalized likelihood estimation which should avoid many of these difficulties is presented. It is demonstrated that each reproducing kernel Hilbert space leads, in a very natural way, to a maximum penalized likelihood estimator and that a well-known class of reproducing kernel Hilbert spaces gives polynomial splines as the nonparametric maximum penalized likelihood estimates.

Demontricher, G. F.↗

NASA NTRS1980

VLA observations of impulsive solar flares at 4.9 GHz

A series of snapshot synthesis maps constructed as a function of time provided a complete mapping of source structure in one dimension and absolute location in two dimensions of five impulsive solar microwave bursts observed during July 15-18, 1978. It was found that all five flares were dominated by a single 12-18 arcsec source located between the H-alpha kernels, near the magnetic neutral line. Spectral data from Sagamore Hill (1.4-35 GHz) show that the spectral turnover at low frequencies results from the absorption by the ambient active region; these results, together with the observed magnetic configuration, suggest a model which could explain the evolution in source structure at 4.9 GHz in terms of magnetic reconnection which starts below the optical depth of the ambient active region equal to 1.

Marsh, K. A.↗

NASA NTRS1995

Error analysis for the ground-based microwave ozone measurements during STOIC

We present a formal error analysis and characterization of the microwave measurements made during the Stratospheric Ozone Intercomparison Campaign (STOIC). The most important error sources are found to be determination of the tropospheric opacity, the pressure-broadening coefficient of the observed line, and systematic variations in instrument response as a function of frequency ('baseline'). Net precision is 4-6% between 55 and 0.2 mbar, while accuracy is 6-10%. Resolution is 8-10 km below 3 mbar and increases to 17km at 0.2 mbar. We show the 'blind' microwave measurements from STOIC and make limited comparisons to other measurements. We use the averaging kernels of the microwave measurement to eliminate resolution and a priori effects from a comparison to SAGE 2. The STOIC results and comparisons are broadly consistent with the formal analysis.

Connor, Brian J.↗