Search NASA⌕ Search

SEARCH · Search NASA

Results for “collocation techniques”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

Distributed caching for processing raw arrays

As applications continue to generate multi-dimensional data at exponentially increasing rates, fast analytics to extract meaningful results is becoming extremely important. The database community has developed array databases that alleviate this problem through a series of techniques. In-situ mechanisms provide direct access to raw data in the original format---without loading and partitioning. Parallel processing scales to the largest datasets. In-memory caching reduces latency when the same data are accessed across a workload of queries. However, we are not aware of any work on distributed caching of multi-dimensional raw arrays. In this paper, we introduce a distributed framework for cost-based caching of multi-dimensional arrays in native format. Given a set of files that contain portions of an array and an online query workload, the framework computes an effective caching plan in two stages. First, the plan identifies the cells to be cached locally from each of the input files by continuously refining an evolving R-tree index. In the second stage, an optimal assignment of cells to nodes that collocates dependent cells in order to minimize the overall data transfer is determined. We design cache eviction and placement heuristic algorithms that consider the historical query workload. A thorough experimental evaluation over two real datasets in three file formats confirms the superiority -- by as much as two orders of magnitude -- of the proposed framework over existing techniques in terms of cache overhead and workload execution time.

Zhao, Weijie↗
DOE OSTI2023

Use of Sobol’ Variance-Based Global Sensitivity Analysis and Multidimensional Legendre Polynomial Fitting for Reduced Order Modeling

Sandia National Laboratories (SNL) has developed a novel reduced order modeling approach. Prioritization of inputs is accomplished using Sobo' indices obtained through a more efficient variance-based global sensitivity analysis. To determine the Sobo' functions, simulated input values are aligned to collocation points to permit the use of Gauss-Lobatto integration, thereby reducing the number of simulation trials needed by more than an order of magnitude compared to standard Monte Carlo approaches. Furthermore, by leveraging the orthogonality of Legendre polynomials in conjunction with those same simulations at the collocation nodes, an efficient fitting method is developed to represent the Sobo' functions from which a reduced order model (ROM) is constructed. The developed method is both more efficient computationally, and the resulting ROM is more accurate. The efficacy of this technique is demonstrated on a nonlinear polynomial test function as well as the nonlinear Ishigami and Sobo' g functions.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Using optimal estimation to retrieve winds from velocity-azimuth display (VAD) scans by a Doppler lidar

Abstract. Low-powered commercially available coherent Doppler lidar (CDL) wind profilers provide continuous measurement of vertical profiles of wind in the lower troposphere, usually close to or up to the top of the planetary boundary layer. The vertical extent of these wind profiles is limited by the availability of scatterers and thus varies substantially throughout the day and from one day to the next. This makes it challenging to develop continuous products that rely on CDL-observed wind profiles. In order to overcome this problem, we have developed a new method for wind profile retrievals from CDL that combines the traditional velocity-azimuth display (VAD) technique with optimal estimation (OE) to provide continuous wind profiles up to 3 km. The new method exploits the level-to-level covariance present in the wind profile to fill in the gaps where the signal-to-noise ratio of the CDL return is too low to provide reliable results using the traditional VAD method. Another advantage of the new method is that it provides the full error covariance matrix of the solution and profiles of information content, which more easily facilitates the assimilation of the observed wind profiles into numerical weather prediction models. This method was tested using yearlong CDL measurements at the Atmospheric Radiation Measurement (ARM) Southern Great Plains (SGP) Central Facility in 2019. Comparison with the ARM operational CDL wind profile product and collocated radiosonde wind measurements shows excellent agreement (R2>0.99) with no degradation in results where the traditional VAD provided a valid solution. In the region where traditional VAD does not provide results, the OE wind speed and wind vector have uncertainties of 3.44 and 4.33 m s−1, respectively. As a result, the new method provides additional information over the standard technique and increases the effective range of existing CDL systems without the need for additional hardware.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Pseudospectral convex optimization for on-ramp merging control of connected vehicles

It can be a daunting task for human drivers to merge into highways because of the intricate vehicle negotiations and potential risk within limited time and space. Connected vehicle (CV) technologies could be a solution to this problem and offer many benefits to the road safety, traffic mobility, and energy efficiency. However, real-time optimal control of CVs is still an open challenge, due to the nonlinear vehicle dynamics, non-convex fuel consumption model, and highly dynamic uncertain inter-vehicle interactions. To tackle these issues, a novel real-time optimal control approach that balances the computational efficiency and solution optimality is proposed for the purpose of onboard application. To this end, the pseudospectral collocation method is integrated with a sequential convex programming approach to develop two new optimization algorithms, which are implemented within a model predictive control (MPC) framework to allow for real-time generation of optimal merging speed profiles. One algorithm leverages the line search technique to improve convergence, and the other benefits from the trust region method for better computational efficiency. The optimality and convergence process of both proposed algorithms are investigated by comparing their solutions with a popular non-linear solver. Furthermore, simulation results show that the proposed methods outperform the benchmark in terms of computational cost, fuel consumption, and traffic efficiency. In particular, the proposed fuel-economy merging rule can save 57.1% fuel consumption on average on four different traffic volumes. Meanwhile, the proposed optimal control algorithms can reduce 2.2% travel time on average comparing to the “first-in-first-out” merging rule.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2017

Mixed Presence Collaboration using Scalable Visualizations in Heterogeneous Display Spaces

Mixed presence collaboration involves remote collaboration between multiple collocated groups. This paper presents the design and results of a user study that focused on mixed presence collaboration using large-scale tiled display walls. The research was conducted in order to compare data synchronization schemes for multi-user visualization applications. Our study compared three techniques for sharing data between display spaces with varying constraints and affordances. The results provide empirical evidence that using data sharing techniques with continuous synchronization between the sites lead to improved collaboration for a search and analysis task between remotely located groups. We have also identified aspects of synchronized sessions that result in increased remote collaborator awareness and parallel task coordination. It is believed that this research will lead to better utilization of large-scale tiled display walls for distributed group work.

data-conferencing↗
DOE OSTI2016

The Nocturnal Avian Migration Experiment Final Campaign Report

Remote sensing techniques are playing a greater role in ornithology, and radar has proven a valuable tool for high resolution, long-term observations of airborne animals. The major disadvantage in radar remote sensing is the current inability to gain taxonomic information from these measurements. One solution is the incorporation of collocated acoustic monitoring that can provide recordings of species-specific nocturnal flight calls of migrating birds in flight. In addition, by taking multichannel recordings of these calls, the position of the calling bird can be calculated and linked to collocated radar measurements.

47 OTHER INSTRUMENTATION↗
DOE OSTI2023

Computation of high-frequency magnetoelastic waves in layered materials

Here, the direct calculation of magnetoelastic wave dispersion in layered media is presented using an efficient, accurate computational technique. The governing, coupled equations for elasticity and magnetism, the Navier and Landau-Lifshitz equations, respectively, are linearized to form a quadratic eigenvalue problem that determines a complex web of wave-number–frequency dispersion branches and their corresponding mode profiles. Numerical discretization of the eigenvalue problem via a spectral collocation method (SCM) is employed to determine the complete dispersion maps for both a single, finite-thickness magnetic layer and a finite magnetic-nonmagnetic double-layer. The SCM, previously used to study elastic waves in nonmagnetic media, is fast, accurate, and adaptable to a variety of sample configurations and geometries. Emphasis is placed on the extremely high-frequency regimes being accessed in ultrafast magnetism experiments. The dispersion maps and modes provide insight into how energy propagates through the coupled system, including how energy can be transferred between elastic- and magnetic-dominated waves as well as between different layers. The numerical computations for a single layer are further understood by a simplified analytical calculation in the high-frequency, exchange-dominated regime where the resonance condition required for energy exchange (an anticrossing) between quasi-elastic and quasi-magnetic dispersion branches is determined. Nonresonant interactions are shown to be well approximated by the dispersion of uncoupled elastic and magnetic waves. The methods and results provide fundamental theoretical tools to model and understand current and future magnetic devices powering spintronic innovation.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

S-OPT: A Points Selection Algorithm for Hyper-Reduction in Reduced Order Models

While projection-based reduced order models can reduce the dimension of full order solutions, the resulting reduced models may still contain terms that scale with the full order dimension. Hyper-reduction techniques are sampling-based methods that further reduce this computational complexity by approximating such terms with a much smaller dimension. The goal of this work is to introduce the points selection algorithm developed by Shin and Xiu as a hyper-reduction method. The selection algorithm was originally proposed as a stochastic collocation method for uncertainty quantification. Since the algorithm aims at maximizing a quantity $\mathcal{S}$ that measures both the column orthogonality and the determinant, we refer to the algorithm as S-OPT. Numerical examples are provided to demonstrate the performance of S-OPT and to compare its performance with a gappy proper orthogonal decomposition (POD) algorithm. Here, we found that using the S-OPT algorithm is shown to predict the full order solutions with higher accuracy than gappy POD especially when the number of sampling points is small, although we note that S-OPT shows slow asymptotic convergence with respect to the number of samples for some applications, e.g., Lagrangian hydrodynamics.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Dual-frequency (Ka-band and G-band) radar estimates of liquid water content profiles in shallow clouds

The profile of the liquid water content (LWC) in clouds provides fundamental information for understanding the internal structure of clouds, their radiative effects, propensity to precipitate, and degree of entrainment and mixing with the surrounding environment. In principle, differential absorption techniques based on coincident dual-frequency radar reflectivity observations have the potential to provide the LWC profile. Previous differential frequency radar reflectivity (DFR) efforts were challenged by the fact that the measurable differential attenuation for small quantities of LWC is usually comparable to the system measurement error. This typically renders the retrieval impractical, as the uncertainty can become many times greater than the retrieved value itself. Theoretically, this drawback can be mitigated following two interconnected approaches: (1) increasing the frequency separation between the dual-frequency radar system to measure greater differential attenuation and (2) increasing the radar operating frequency to reduce the instrument measurement random error. Our recently developed 239 GHz radar was deployed during the Eastern Pacific Cloud Aerosol Precipitation Experiment (EPCAPE) along with a variety of collocated remote sensing and in situ instruments. We have combined Ka-band (35 GHz) and G-band (239 GHz) observations to retrieve the LWC from more than 100 vertical profiles of shallow clouds with typical amounts of LWC smaller than 1 g m -3 . We theoretically and experimentally demonstrate that the Ka-band and G-band pair of frequencies offers at least a 65 % relative improvement in the LWC retrieval sensitivity compared to previous works reported in the literature using lower-frequency radars. This new technique provides a missing capability to determine the LWC in the challenging low liquid water path (LWP) range (< 200 g m -2 ) and suggests a way forward to characterize microphysical and dynamical processes more precisely in shallow clouds.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2013

MULTILEVEL ACCELERATION OF STOCHASTIC COLLOCATION METHODS FOR PDE WITH RANDOM INPUT DATA

Stochastic Collocation (SC) methods for stochastic partial differential equa- tions (SPDEs) suffer from the curse of dimensionality, whereby increases in the stochastic dimension cause an explosion of computational effort. To combat these challenges, multilevel approximation methods seek to decrease computational complexity by balancing spatial and stochastic discretization errors. As a form of variance reduction, multilevel techniques have been successfully applied to Monte Carlo (MC) methods, but may be extended to accelerate other methods for SPDEs in which the stochastic and spatial degrees of freedom are de- coupled. This article presents general convergence and computational complexity analysis of a multilevel method for SPDEs, demonstrating its advantages with regard to standard, single level approximation. The numerical results will highlight conditions under which multilevel sparse grid SC is preferable to the more traditional MC and SC approaches.

Webster, Clayton G↗
DOE OSTI2020

A Polarimetric Radar Analysis of Ice Microphysical Processes in Melting Layers of Winter Storms Using S-Band Quasi-Vertical Profiles

Quasi-vertical profiles (QVPs) obtained from a database of U.S. WSR-88D data are used to document polarimetric characteristics of the melting layer (ML) in cold-season storms with high vertical resolution and accuracy. A polarimetric technique to define the top and bottom of the ML is first introduced. Using the QVPs, statistical relationships are developed to gain insight into the evolution of microphysical processes above, within, and below the ML, leading to a statistical polarimetric model of the ML that reveals characteristics that reflectivity data alone are not able to provide, particularly in regions of weak reflectivity factor at horizontal polarization Z H . QVP ML statistics are examined for two regimes in the ML data: Z H ≥ 20 dB Z and Z H < 20 dB Z . Regions of Z H ≥ 20 dB Z indicate locations of MLs collocated with enhanced differential reflectivity Z DR and reduced copolar correlation coefficient ρ hv , while for Z H < 20 dB Z a well-defined ML is difficult to discern using Z H alone. Evidence of large Z DR up to 4 dB, backscatter differential phase δ up to 8°, and low ρ hv down to 0.80 associated with lower Z H (from -10 to 20 dB Z ) in the ML is observed when pristine, nonaggregated ice falls through it. Positive correlation is documented between maximum specific differential phase K DP and maximum Z H in the ML; these are the first QVP observations of K DP in MLs documented at S band. Negative correlation occurs between minimum ρ hv in the ML and ML depth and between minimum ρ hv in the ML and the corresponding enhancement of Z H (Δ Z H = Z H max - Z H rain ).

Meteorology & Atmospheric Sciences↗
DOE OSTI2016

Simplex-stochastic collocation method with improved scalability

The Simplex-Stochastic Collocation (SSC) method is a robust tool used to propagate uncertain input distributions through a computer code. However, it becomes prohibitively expensive for problems with dimensions higher than 5. The main purpose of this paper is to identify bottlenecks, and to improve upon this bad scalability. In order to do so, we propose an alternative interpolation stencil technique based upon the Set-Covering problem, and we integrate the SSC method in the High-Dimensional Model-Reduction framework. In addition, we address the issue of ill-conditioned sample matrices, and we present an analytical map to facilitate uniformly-distributed simplex sampling.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

Parasol: An Architecture for Cross-Cloud Federated Graph Querying

Large scale data fusion of multiple datasets can often provide in- sights that examining datasets individually cannot. However, when these datasets reside in different data centers and cannot be collocated due to technical, administrative, or policy barriers, a unique set of problems arise that hamper querying and data fusion. To ad- dress these problems, a system and architecture named Parasol is presented that enables federated queries over graph databases residing in multiple clouds. Parasol’s design is flexible and requires only minimal assumptions for participant clouds. Query optimization techniques are also described that are compatible with Parasol’s lightweight architecture. Experiments on a prototype implementation of Parasol indicate its suitability for cross-cloud federated graph queries.

Graph database, federated database, data fusion, s↗
DOE OSTI2012

Derivation of new 3D discrete ordinate equations

The Sn equations have been the workhorse of deterministic radiation transport calculations for many years. Here we derive two new angular discretizations of the 3D transport equation. The first set of equations, derived using Lagrange interpolation and collocation, retains the classical Sn structure, with the main difference being how the scattering source is calculated. Because of the formal similarity with the classical S n equations, it should be possible to modify existing computer codes to take advantage of the new formulation. In addition, the new S n-like equations correctly capture delta function scattering. The second set of equations, derived using a Galerkin technique, does not retain the classical Sn structure because the streaming term is not diagonal. However, these equations can be cast into a form similar to existing methods developed to reduce ray effects. Numerical investigation of both sets of equations is under way. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2023

Additive Polynomial Time Integrators, Part I: Framework and Fully Implicit-Explicit Collocation Methods

In this paper we generalize the polynomial time integration framework to additively partitioned initial value problems. The framework we present is general and enables the construction of many new families of additive integrators with arbitrary order-of-accuracy and varying degree of implicitness. In this first work, we focus on a new class of implicit-explicit polynomial block methods that are based on fully implicit Runge–Kutta methods with Radau nodes and that possess high stage order. Here, we show that the new fully implicit-explicit (FIMEX) integrators have improved stability compared to existing IMEX Runge–Kutta methods, while also being more computationally efficient due to recent developments in preconditioning techniques for solving the associated systems of nonlinear equations. For PDEs on periodic domains where the implicit component is trivial to invert, we will show how parallelization of the right-hand side evaluations can be exploited to obtain significant speedup compared to existing serial IMEX Runge–Kutta methods. For parallel (in space) finite element discretizations, the new methods can achieve orders of magnitude better accuracy than existing IMEX Runge–Kutta methods and/or achieve a given accuracy several times times faster in terms of computational runtime.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2012

An Adaptive Sparse-Grid High-Order Stochastic Collocation Method for Bayesian Inference in Groundwater Reactive Transport Modeling

Although Bayesian analysis has become vital to the quantification of prediction uncertainty in groundwater modeling, its application has been hindered due to the computational cost associated with numerous model executions needed for exploring the posterior probability density function (PPDF) of model parameters. This is particularly the case when the PPDF is estimated using Markov Chain Monte Carlo (MCMC) sampling. In this study, we develop a new approach that improves computational efficiency of Bayesian inference by constructing a surrogate system based on an adaptive sparse-grid high-order stochastic collocation (aSG-hSC) method. Unlike previous works using first-order hierarchical basis, we utilize a compactly supported higher-order hierar- chical basis to construct the surrogate system, resulting in a significant reduction in the number of computational simulations required. In addition, we use hierarchical surplus as an error indi- cator to determine adaptive sparse grids. This allows local refinement in the uncertain domain and/or anisotropic detection with respect to the random model parameters, which further improves computational efficiency. Finally, we incorporate a global optimization technique and propose an iterative algorithm for building the surrogate system for the PPDF with multiple significant modes. Once the surrogate system is determined, the PPDF can be evaluated by sampling the surrogate system directly with very little computational cost. The developed method is evaluated first using a simple analytical density function with multiple modes and then using two synthetic groundwater reactive transport models. The groundwater models represent different levels of complexity; the first example involves coupled linear reactions and the second example simulates nonlinear ura- nium surface complexation. The results show that the aSG-hSC is an effective and efficient tool for Bayesian inference in groundwater modeling in comparison with conventional MCMC sim- ulations. The computational efficiency is expected to be more beneficial to more computational expensive groundwater problems.

Zhang, Guannan↗
DOE OSTI2015

How Well Can We Estimate Areal-Averaged Spectral Surface Albedo from Ground-Based Transmission in an Atlantic Coastal Area?

Areal-averaged albedos are particularly difficult to measure in coastal regions, because the surface is not homogenous, consisting of a sharp demarcation between land and water. With this difficulty in mind, we evaluate a simple retrieval of areal-averaged surface albedo using ground-based measurements of atmospheric transmission alone under fully overcast conditions. To illustrate the performance of our retrieval, we find the areal-averaged albedo using measurements from the Multi-Filter Rotating Shadowband Radiometer (MFRSR) at five wavelengths (415, 500, 615, 673, and 870 nm). These MFRSR data are collected at a coastal site in Graciosa Island, Azores supported by the U.S. Department of Energy’s (DOE’s) Atmospheric Radiation Measurement (ARM) Program. The areal-averaged albedos obtained from the MFRSR are compared with collocated and coincident Moderate Resolution Imaging Spectroradiometer (MODIS) white-sky albedo at four nominal wavelengths (470, 560, 670 and 860 nm). These comparisons are made during a 19-month period (June 2009 - December 2010). We also calculate composite-based spectral values of surface albedo by a weighted-average approach using estimated fractions of major surface types observed in an area surrounding this coastal site. Taken as a whole, these three methods of finding albedo show spectral and temporal similarities, and suggest that our simple, transmission-based technique holds promise, but with estimated errors of about ±0.03. Additional work is needed to reduce this uncertainty in areas with inhomogeneous surfaces.

Kassianov, Evgueni I.↗
DOE OSTI2019

Vectoring into Potential Blind Geothermal Systems in the Granite Springs Valley Area, Western Nevada: Application of the Play Fairway Analysis at Multiple Scales

The Great Basin region of the western USA is capable of generating much greater amounts of geothermal energy than currently produced. Most geothermal resources in this region are blind, and thus the favorable characteristics for geothermal activity must be synthesized and methodologies developed to discover new commercial-grade systems. The geothermal play fairway concept involves integration of multiple parameters indicative of geothermal activity as a means of identifying promising areas for new development. In the Nevada play fairway project, nine geologic, geochemical, and geophysical parameters were initially synthesized to produce a new geothermal potential map of 96,000 km2. Granite Springs Valley in western Nevada is a particularly promising site selected for detailed study. It contains several favorable structural settings, including terminations of major Quaternary normal faults and fault intersections. Geologic, geophysical, and geochemical techniques were employed to define the most likely sites for high permeability and select drilling targets for temperature-gradient holes. Local and intermediate permeability models were revised to reflect results of detailed analyses and generate new detailed play fairway maps of the area. The most promising site lies in the northeastern part of the basin directly east of Adobe Flat, where the horse-tailing termination of a major normal fault (as revealed by new gravity data), opaline sinter deposits, warm shallow wells, and a low-resistivity anomaly are collocated. Geothermometry suggests a blind system with temperatures as high as ~130°C in this area. Four new temperature-gradient holes document temperatures of ~80°C at ~150 m depth. Lessons learned in the detailed studies of this project include: 1) initially identified sites commonly include multiple favorable settings at a finer scale; 2) promising sites in Cenozoic basins cannot be well defined without detailed geophysical surveys; and 3) play fairway analysis is critical at multiple scales, providing a means to select regional prospects as well as vectoring into drilling targets at individual sites.

15 GEOTHERMAL ENERGY↗