Search NASA⌕ Search

SEARCH · Search NASA

Results for “Function Estimation”

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 91 records · Page 5

nbodykit: An Open-source, Massively Parallel Toolkit for Large-scale Structure

We present nbodykit, an open-source, massively parallel Python toolkit for analyzing large-scale structure (LSS) data. Using Python bindings of the Message Passing Interface, we provide parallel implementations of many commonly used algorithms in LSS. nbodykit is both an interactive and scalable piece of scientific software, performing well in a supercomputing environment while still taking advantage of the interactive tools provided by the Python ecosystem. Existing functionality includes estimators of the power spectrum, two- and three-point correlation functions, a friends-of-friends grouping algorithm, mock catalog creation via the halo occupation distribution technique, and approximate N-body simulations via the FastPM scheme. The package also provides a set of distributed data containers, insulated from the algorithms themselves, that enables nbodykit to provide a unified treatment of both simulation and observational data sets. nbodykit can be easily deployed in a high-performance computing environment, overcoming some of the traditional difficulties of using Python on supercomputers. We provide performance benchmarks illustrating the scalability of the software. The modular, component-based approach of nbodykit allows researchers to easily build complex applications using its tools. The package is extensively documented at http://nbodykit.readthedocs.io, which also includes an interactive set of example recipes for new users to explore. As open-source software, we hope nbodykit provides a common framework for the community to use and develop in confronting the analysis challenges of future LSS surveys.

79 ASTRONOMY AND ASTROPHYSICS↗

Coding the Computing Continuum: Fluid Function Execution in Heterogeneous Computing Environments

Advances in network technologies have greatly decreased barriers to accessing physically distributed computers. This newfound accessibility coincides with increasing hardware specialization, creating exciting new opportunities to dispatch workloads to the best resource for a specific purpose, rather than those that are closest or most easily accessible. We present Delta, a service designed to intelligently schedule function-based workloads across a distributed set of heterogeneous computing resources. Delta implements an extensible architecture in which different predictors and scheduling algorithms can be integrated to provide dynamically evolving estimates of function execution times on different resources-estimates that can be used to determine the most appropriate location for execution. We describe predictors for function runtime, data transfer time, and cold-start resource provisioning and configuration delay; dynamic learning methods that update predictor models over time; and scheduling strategies that take into account both function and endpoint information. We show that these methods can halve workload makespan when compared with a strategy that selects the fastest resource, and decrease makespan by a factor of five when compared to a round robin strategy, when deployed on a heterogeneous testbed with resources ranging from a Raspberry Pi to a GPU node in an academic cloud.

Computing continuum↗

Cluster Sampling Filters for Non-Gaussian Data Assimilation

This paper presents a fully non-Gaussian filter for sequential data assimilation. The filter is named the “cluster sampling filter”, and works by directly sampling the posterior distribution following a Markov Chain Monte-Carlo (MCMC) approach, while the prior distribution is approximated using a Gaussian Mixture Model (GMM). Specifically, a clustering step is introduced after the forecast phase of the filter, and the prior density function is estimated by fitting a GMM to the prior ensemble. Using the data likelihood function, the posterior density is then formulated as a mixture density, and is sampled following an MCMC approach. Four versions of the proposed filter, namely C ℓ MCMC , C ℓ HMC , MC- C ℓ HMC , and MC- C ℓ HMC are presented. C ℓ MCMC uses a Gaussian proposal density to sample the posterior, and C ℓ HMC is an extension to the Hamiltonian Monte-Carlo (HMC) sampling filter. MC- C ℓ MCMC and MC- C ℓ HMC are multi-chain versions of the cluster sampling filters C ℓ MCMC and C ℓ HMC respectively. The multi-chain versions are proposed to guarantee that samples are taken from the vicinities of all probability modes of the formulated posterior. The new methodologies are tested using a simple one-dimensional example, and a quasi-geostrophic (QG) model with double-gyre wind forcing and bi-harmonic friction. Here, numerical results demonstrate the usefulness of using GMMs to relax the Gaussian prior assumption especially in the HMC filtering paradigm.

97 MATHEMATICS AND COMPUTING↗

Functional diversity and redundancy across fish gut, sediment and water bacterial communities

Here, this work explores the functional diversity and redundancy in a bacterial metacommunity constituted of three habitats (sediment, water column and fish gut) in a coastal lagoon under anthropogenic pressure. Comprehensive functional gene arrays covering a wide range of ecological processes and stress resistance genes to estimate the functional potential of bacterial communities were used. Then, diversity partitioning was used to characterize functional diversity and redundancy within (α), between (β) and across (γ) habitats. It was showed that all local communities exhibit a highly diversified potential for the realization of key ecological processes and resistance to various environmental conditions, supporting the growing evidence that macro-organisms microbiomes harbour a high functional potential and are integral components of functional gene dynamics in aquatic bacterial metacommunities. Several levels of functional redundancy at different scales of the bacterial metacommunity were observed (within local communities, within habitats and at the metacommunity level). The findings indicated a high potential for the realization of spatial ecological insurance within this ecosystem, that is, the functional compensation among microorganisms for the realization and maintenance of key ecological processes, within and across habitats. Finally, the role of macro-organisms as dispersal vectors of microbes and their potential influence on marine metacommunity dynamics were discussed.

59 BASIC BIOLOGICAL SCIENCES↗

climextremes v0.2.1

Software package for R and Python for analysis of extremes in climate data. Functions for fitting GEV and POT (via point process fitting) models for extremes in climate data, providing return values, return probabilities, and return periods for stationary and nonstationary models. Also provides differences in return values and differences in log return probabilities for contrasts of covariate values. Functions for estimating risk ratios for event attribution analyses, including uncertainty. Under the hood, many of the functions use functions from 'extRemes', including for fitting the statistical models.

Krishnan, Harinarayan↗

Estimation of conditional cumulative incidence functions under generalized semiparametric regression models with missing covariates, with application to analysis of biomarker correlates in vaccine trials

Herein, this article presents generalized semiparametric regression models for conditional cumulative incidence functions with competing risks data when covariates are missing by sampling design or happenstance. A doubly robust augmented inverse probability weighted (AIPW) complete-case approach to estimation and inference is investigated. This approach modifies IPW complete-case estimating equations by exploiting the key features in the relationship between the missing covariates and the phase-one data to improve efficiency. An iterative numerical procedure is derived to solve the nonlinear estimating equations. The asymptotic properties of the proposed estimators are established. A simulation study examining the finite-sample performances of the proposed estimators shows that the AIPW estimators are more efficient than the IPW estimators. The developed method is applied to the RV144 HIV-1 vaccine efficacy trial to investigate vaccine-induced IgG binding antibodies to HIV-1 as correlates of acquisition of HIV-1 infection while taking account of whether the HIV-1 sequences are near or far from the HIV-1 sequences represented in the vaccine construct.

97 MATHEMATICS AND COMPUTING↗

CLIP Representation Enhanced Predicate Estimation

This software provides functions to perform visual predicate estimation by leveraging the pre-trained vision-language model CLIP. In addition to training the prediction estimator with a labeled dataset, it also performs performance evaluation and visualization of results.

Jayaraman Thiagarajan, Jayaraman↗

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↗

QED corrections to the thermal neutrino interaction rate

Abstract Motivated by precision computations of neutrino decoupling at MeV temperatures, we show how QED corrections to the thermal neutrino interaction rate can be related to the electron-positron spectral function as well as an effective$$\overline{\nu }\nu \gamma $$vertex. The spectral function is needed both in a timelike and in a spacelike domain, and for both of its physical polarization states (transverse and longitudinal with respect to spatial momentum). Incorporating an NLO evaluation of this spectral function, an estimate of the$$\overline{\nu }\nu \gamma $$vertex, and HTL resummation of scatterings mediated by soft Bose-enhancedt-channel photons, we compute the interaction rate as a function of the neutrino momentum and flavour. Effects on the – (0. . .2)% level are found, noticeably smaller than a previous estimate of a related quantity.

Physics↗

Wafer-scale waveguide sidewall roughness scattering loss characterization by image processing

Photonic integrated circuits (PICs) are vital for developing affordable, high-performance optoelectronic devices that can be manufactured at an industrial scale, driving innovation and efficiency in various applications. Optical loss of modes in thin film waveguides and devices is a critical measure of their performance. Thin film growth, lithography, masking, and etching processes are imperfect processes that introduce significant sidewall and top-surface roughness and cause dominating optical losses in waveguides and photonic structures. This roughness, as perturbations couple light from guided to far-field radiation modes, leads to scattering losses that can be estimated from theoretical models. Typically, with UV-based lithography, sidewall roughness is significantly larger than wafer-top surface roughness. Atomic force microscopy (AFM) imaging measurement gives a 3D and high-resolution roughness profile, but the measurement is inconvenient, costly, and unscalable for large-scale PICs and at wafer-scale. Here, we evaluate the sidewall roughness profile based on 2D high-resolution scanning electron microscope (SEM) imaging. We characterized the loss on two homemade nitride and oxide films on 3-inch silicon wafers with 12 waveguide devices on each and correlated the scattering loss estimated from a 2D image-based sidewall profile and theoretical Payne model. The lowest loss of guided fundamental transverse electric (TE 0 ) mode is found at 0.075 dB/cm at 633 nm across 24 devices, a record at visible wavelength. Our work shows 100% success (edge continuity span exceeding 95% of image width/height) in edge detection in image processing of all images to estimate autocorrelation function and optical mode loss. These demonstrations offer valuable insights into waveguide sidewall roughness and a comparison of experimental and 2D SEM image processing based loss estimations with applications in loss characterization at wafer-scale PICs.

Khurana, Mohit (ORCID:0000000207898672)↗

Photovoltaic (PV) System Levelized Cost of Energy (LCOE) Evaluation with Grid Support Function Valuation and Service Lifetime Estimation

Photovoltaic (PV) systems play a critical role in renewable energy resource grid integration, and levelized cost of energy (LCOE) is commonly used to evaluate PV system feasibility in modern power grids. In this work, a revised PV system LCOE calculation model is derived to quantify the potential of LCOE reduction. Particularly, the grid support functions are valuated to offset the investment and operation costs of PV systems, which thereby reduces the LCOE. Meanwhile, PV system service lifetime is also estimated with the derived PV inverter reliability model, considering the critical components (i.e., semiconductor devices and capacitors). The case studies with field datasets are conducted to validate the effectiveness of the developed LCOE calculation model.

Zhao, Shijia↗

Propagation of priors for more accurate and efficient spectroscopic functional fits and their application to ferroelectric hysteresis

Multi-dimensional spectral-imaging is a mainstay of the scanning probe and electron microscopies, micro-Raman, and various forms of chemical imaging. In many cases, individual spectra can be fit to a specific functional form, with the model parameter maps, providing direct insight into material properties. Since spectra are often acquired across a spatial grid of points, spatially adjacent spectra are likely to be similar to one another; yet, this fact is almost never used when considering parameter estimation for functional fits. On datasets tried here, we show that by utilizing proximal information, whether it be in the spatial or spectral domains, it is possible to improve the reliability and increase the speed of such functional fits by ~2-3x, as compared to random priors. We explore and compare three distinct new methods: (1) spatially averaging neighborhood spectra, and propagating priors based on functional fits to the averaged case, (2) hierarchical clustering-based methods where spectra are grouped hierarchically based on response, with the priors propagated progressively down the hierarchy, and (3) regular clustering without hierarchical methods with priors propagated from fits to cluster means. Our results highlight that utilizing spatial and spectral neighborhood information is often critical for accurate parameter estimation in noisy environments, which we show for ferroelectric hysteresis loops acquired on a prototypical PbTiO3 thin film with piezoresponse spectroscopy. This method is general and applicable to any spatially measured spectra where functional forms are available. Examples include exploring the superconducting gap with tunneling spectroscopy, using the Dynes formula, or current-voltage curve fits in conductive atomic force microscopy mapping. Here we explore the problem for ferroelectric hysteresis, which, given its large parameter space, constitutes a more difficult task than, for example, fitting current-voltage curves with a Schottky emission formula.

42 ENGINEERING↗

A simplified framework for prediction of sensor network coverage in real-time structural health monitoring of plate-like structures

In this article, the coverage area prediction of piezoelectric sensor network for detecting a specific type of under-surface crack in plate-like structures is addressed. In particular, this article proposes a simplified framework to estimate the coverage of any given sensor network arrangement when a critical defect is known. Based on numerical results from finite element methods (FEM), a simplified framework to estimate coverage area of any given network arrangement is developed. Using such a simplified framework, one can avoid time-consuming procedure of evaluating numerous FEM models in estimating sensor network coverage. Back-scatter fields of partial cracks are estimated using a proposed function, whose parameters are estimated from the results of a limited number of FEM simulations. The proposed function efficiently predicts the back-scattered field of any combination of transmitters and receivers for a given crack geometry. Superposition is used to estimate the coverage area of an arbitrary piezoelectric (e.g., PZT) sensor network. It is shown that the coverage area of a sensor network depends on both sensor network geometry and defect properties (e.g., crack inclination) and it is not necessarily a linear function of the number of sensors. Furthermore, it is shown that the network arrangement has an important effect on the geometry of the coverage area. Experimental results of a network of 14 PZTs in two clusters confirm the accuracy of the method.

Engineering↗

Implications of theoretical analysis to explore the functional core dimension in one-rod core microstructured optical fibers

A new development in optical fiber technology that has attracted interest for multitude of novel and potential applications is an all-in-silica microstructured optical fibre (MOF) with one-rod core. We aim to estimate the functional core dimension for one-rod core triangular MOF (T-MOF) from its geometrical characteristics by using an alternative analytical field model. Moreover, to scrutinize the validity of the choice of effective core dimension (by adapting V-parameter, i.e., the normalized frequency) for T-MOF, we have explored the rudimentary propagation characteristics of the T-MOF in single mode guidance realm. We made comparisons with those which are based on the experimental and numerical studies. Relative errors are also reported.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

encore : an O ( N g2) estimator for galaxy N -point correlation functions

ABSTRACT We present a new algorithm for efficiently computing the N-point correlation functions (NPCFs) of a 3D density field for arbitrary N. This can be applied both to a discrete spectroscopic galaxy survey and a continuous field. By expanding the statistics in a separable basis of isotropic functions built from spherical harmonics, the NPCFs can be estimated by counting pairs of particles in space, leading to an algorithm with complexity $\mathcal {O}(N_\mathrm{g}^2)$ for Ng particles, or $\mathcal {O}(N_\mathrm{FFT}\log N_\mathrm{FFT})$ when using a Fast Fourier Transform with NFFT grid-points. In practice, the rate-limiting step for N > 3 will often be the summation of the histogrammed spherical harmonic coefficients, particularly if the number of radial and angular bins is large. In this case, the algorithm scales linearly with Ng. The approach is implemented in the encore code, which can compute the 3PCF, 4PCF, 5PCF, and 6PCF of a BOSS-like galaxy survey in ${\sim}100$ CPU-hours, including the corrections necessary for non-uniform survey geometries. We discuss the implementation in depth, along with its GPU acceleration, and provide practical demonstration on realistic galaxy catalogues. Our approach can be straightforwardly applied to current and future data sets to unlock the potential of constraining cosmology from the higher point functions.

79 ASTRONOMY AND ASTROPHYSICS↗

Subsurface Hydrogen Assessment, Storage, and Technology Acceleration (SHASTA) - Hydrogen Estimator for Logistical Planning (HELP) User’s Manual

This user’s manual describes version 1.0 of the SHASTA-HELP (Subsurface Hydrogen Assessment, Storage, and Technology Acceleration (SHASTA) - Hydrogen Estimator for Logistical Planning (HELP)) tool and provides instructions for use. The purpose of the SHASTA-HELP tool is to integrate the suite of tools developed by the SHASTA research team in a web-based framework for ease of access and more advanced, integrated analysis. This tool contains functionality for estimating the storage potential of pure and blended natural gas-hydrogen mixtures in various subsurface formations and is intended to be used for pre-characterization or site screening purposes. Detailed storage analysis should be conducted by full reservoir simulation models.

08 HYDROGEN↗

Secondary battery management system

A method and system for managing a battery system. The method including receiving at least one measured characteristic of the battery over a pre-defined time horizon from the at least one sensor, receiving at least one estimated characteristic of the battery from a electrochemical-based battery model based on differential algebraic equations, determining a cost function of a Moving Horizon Estimation based on the at least one measured characteristic and the at least one estimated characteristic, updating the electrochemical-based battery model based on the cost function, estimating at least one state of the at least one battery cell by applying the electrochemical-based battery model, and regulating at least one of charging or discharging of the battery based on the estimation of the at least one state of the at least one battery cell.

MirTabatabaei, Anahita↗