Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 469 records · Page 26

A structure function representation theorem with applications to frequency stability estimation

Random processes with stationary nth differences serve as models for oscillator phase noise. A theorem which obtains the structure function (covariance of the nth differences) of such a process in terms of the differences of a single function of one time variable is proven. In turn, this function can easily be obtained from the spectral density of the process. The theorem is used for computing the variance of two estimators of frequency stability.

Greenhall, C. A.↗

Probabilistic assessment of cloud fraction using Bayesian blending of independent datasets: Feasibility study of a new method

We describe and evaluate a novel method to blend two observed cloud fraction (CF) datasets through Bayesian posterior estimation. The research reported here is a feasibility study designed to explore the method. In this proof‐of‐concept study, we illustrate the approach using specific observational datasets from the U. S. Department of Energy Atmospheric Radiation Measurement Program's Southern Great Plains site in the central United States, but the method is quite general and is readily applicable to other datasets. The total sky image (TSI) camera observations are used to determine the prior distribution. A regression model and the active remote sensing of clouds (ARSCL) radar/lidar observations are used to determine the likelihood function. The posterior estimate is a probability density function (pdf) of the CF whose mean is taken to be the optimal blend of the two observations. The data at hourly, daily, 5‐day, monthly, and annual time scales are considered. Some physical and probabilistic properties of the CFs are explored from radar/lidar, camera, and satellite observations and from simulations using the Community Atmosphere Model (CAM5). Our results imply that (a) the Beta distribution is a reasonable model for CF for both short‐ and long‐time means, the 5‐day data are skewed right, and the annual data are almost normally distributed, and (b) the Bayesian method developed successfully yields a pdf of CF, rather than a deterministic CF value, and it is feasible to blend the TSI and ARSCL data with a capability for bias correction.

58 GEOSCIENCES↗

Public charging infrastructure for plug-in electric vehicles: What is it worth?

A lack of charging infrastructure is an important barrier to the growth of the plug-in electric vehicle (PEV) market. Public charging infrastructure has tangible and intangible value, such as reducing range anxiety or building confidence in the future of the PEV market. Quantifying the value of public charging infrastructure can inform analysis of investment decisions and can help predict the impact of charging infrastructure on future PEV sales. Estimates of willingness to pay (WTP) based on stated preference surveys are limited by consumers' lack of familiarity with PEVs. As an alternative, we focus on quantifying the tangible value of public PEV chargers in terms of their ability to displace gasoline use for PHEVs and to enable additional electric (e-) vehicle miles for BEVs, thereby mitigating the limitations of shorter range and longer recharging time. Simulation studies provide data that can be used to quantify e-miles enabled by public chargers and the value of additional e-miles can be inferred from econometric estimates of WTP for increased vehicle range. Functions are synthesized that estimate the WTP for public charging infrastructure by plug-in hybrid and battery electric vehicles, conditional on vehicle range, annual vehicle travel, pre-existing charging infrastructure, energy prices, vehicle efficiency, and household income. A case study based on California's public charging network in 2017 indicates that, to the purchaser of a new BEV with a 100-mile range and home recharging, existing public fast chargers are worth about $1500 for intraregional travel, and fast chargers along intercity routes are valued at over $6500.

33 ADVANCED PROPULSION SYSTEMS↗

THE STRUCTURE OF THE MERGING RCS 231953+00 SUPERCLUSTER AT z {approx} 0.9

The RCS 2319+00 supercluster is a massive supercluster at z = 0.9 comprising three optically selected, spectroscopically confirmed clusters separated by <3 Mpc on the plane of the sky. This supercluster is one of a few known examples of the progenitors of present-day massive clusters (10{sup 15} M{sub Sun} by z {approx} 0.5). We present an extensive spectroscopic campaign carried out on the supercluster field resulting, in conjunction with previously published data, in 1961 high-confidence galaxy redshifts. We find 302 structure members spanning three distinct redshift walls separated from one another by {approx}65 Mpc ({Delta} z = 0.03). The component clusters have spectroscopic redshifts of 0.901, 0.905, and 0.905. The velocity dispersions are consistent with those predicted from X-ray data, giving estimated cluster masses of {approx}10{sup 14.5}-10{sup 14.9} M{sub Sun }. The Dressler-Shectman test finds evidence of substructure in the supercluster field and a friends-of-friends analysis identified five groups in the supercluster, including a filamentary structure stretching between two cluster cores previously identified in the infrared by Coppin et al. The galaxy colors further show this filamentary structure to be a unique region of activity within the supercluster, comprised mainly of blue galaxies compared to the {approx}43%-77% red-sequence galaxies present in the other groups and cluster cores. Richness estimates from stacked luminosity function fits result in average group mass estimates consistent with {approx}10{sup 13} M{sub Sun} halos. Currently, 22% of our confirmed members reside in {approx}> 10{sup 13} M{sub Sun} groups/clusters destined to merge onto the most massive cluster, in agreement with the massive halo galaxy fractions important in cluster galaxy pre-processing in N-body simulation merger tree studies.

79 ASTRONOMY AND ASTROPHYSICS↗

Uncertainty Quantification for Point Spread Functions in X-Ray Radiography

X-ray radiographers have tools for estimating the point spread function (PSF) which determines the extent of the blur in images. Knowledge of uncertainties in their estimations can be extremely valuable, and existing methods do not provide this information. We have developed a new technique using a Bayesian framework that will provide useful uncertainties to the estimated PSF.

point spread function↗

Mars - VLA observations of the northern hemisphere and the north polar region at wavelengths of 2 and 6 cm

Calculations based on 2- and 6-cm observations of Mars with the A configuration of the VLA have yielded a whole-disk effective dielectric constant of 2.34 + or - 0.05, implying a subsurface density of 1.24 + or - 0.11 g/cu cm at 2 cm, as well as 1.45 + or - 0.10 g/cu cm effective density and 2.70 + or - 0.10 dielectric constant at 6 cm. These parameters have also been estimated as a function of latitude over the 15 deg S - 60 deg N range; subsurface radio absorption length was estimated to be about 15 wavelengths at most of these latitudes. Most of the subsurface density calculations yielded results in the 1-2-g/cu cm range, implying that the subsurface is not very different from the surface observed by Viking and Mariner spacecraft; the decrease in correlation with depth is in keeping with slow variation of the subsurface in the near-subsurface region.

Rudy, Donald J.↗

A Cellular Automaton Simulation for Predicting Phase Evolution in Solid-State Reactions

New computational tools for solid-state synthesis recipe design are needed in order to accelerate the experimental realization of novel functional materials proposed by high-throughput materials discovery workflows. This work contributes a cellular automaton simulation framework for predicting the time-dependent evolution of intermediate and product phases during solid-state reactions as a function of precursor choice and amount, reaction atmosphere, and heating profile. The simulation captures the effects of reactant particle spatial distribution, particle melting, and reaction atmosphere. Reaction rates based on rudimentary kinetics are estimated using density functional theory data from the Materials Project and machine learning estimators for the melting point and the vibrational entropy component of the Gibbs free energy. The resulting simulation framework allows for the prediction of the likely outcome of a reaction recipe before any experiments are performed. We analyze five experimental solid-state recipes for BaTiO 3 , CaZrN 2 , and YMnO 3 found in the literature to illustrate the performance of the model in capturing reaction selectivity and reaction pathways as a function of temperature and precursor choice. This simulation framework offers an easier way to optimize existing recipes, aid in the identification of intermediates, and design effective recipes for yet unrealized inorganic solids in silico .

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

GMLS-Nets: A Framework for Learning from Unstructured Data

Data fields sampled on irregularly spaced points arise in many applications in the sciences and engineering. For regular grids, Convolutional Neural Networks (CNNs) have been successfully used to gaining benefits from weight sharing and invariances. We generalize CNNs by introducing methods for data on unstructured point clouds based on Generalized Moving Least Squares (GMLS). GMLS is a nonparametric technique for estimating linear bounded functionals from scattered data, and has recently been used in the literature for solving partial differential equations. By parameterizing the GMLS estimator, we obtain learning methods for operators with unstructured stencils. In GMLS-Nets the necessary calculations are local, readily parallelizable, and the estimator is supported by a rigorous approximation theory. We show how the framework may be used for unstructured physical data sets to perform functional regression to identify associated differential operators and to regress quantities of interest. The results suggest the architectures to be an attractive foundation for data-driven model development in scientific machine learning applications.

42 ENGINEERING↗

GMLS-Nets: Scientific Machine Learning Methods for Unstructured Data

Data fields sampled on irregularly spaced points arise in many applications in the sciences and engineering. For regular grids, Convolutional Neural Networks (CNNs) have been successfully used to gaining benefits from weight sharing and invariances. We generalize CNNs by introducing methods for data on unstructured point clouds based on Generalized Moving Least Squares (GMLS). GMLS is a nonparametric technique for estimating linear bounded functionals from scattered data, and has recently been used in the literature for solving partial differential equations. By parameterizing the GMLS estimator, we obtain learning methods for operators with unstructured stencils. In GMLS-Nets the necessary calculations are local, readily parallelizable, and the estimator is supported by a rigorous approximation theory. We show how the framework may be used for unstructured physical data sets to perform functional regression to identify associated differential operators and to regress quantities of interest. The results suggest the architectures to be an attractive foundation for data-driven model development in scientific machine learning applications

97 MATHEMATICS AND COMPUTING↗

GMLS-Nets: A Framework for Learning from Unstructured Data

Data fields sampled on irregularly spaced points arise in many applications in the sciences and engineering. For regular grids, Convolutional Neural Networks (CNNs) have been successfully used to gaining benefits from weight sharing and invariances. We generalize CNNs by introducing methods for data on unstructured point clouds based on Generalized Moving Least Squares (GMLS). GMLS is a non-parametric technique for estimating linear bounded functionals from scattered data, and has recently been used in the literature for solving partial differential equations. By parameterizing the GMLS estimator, we obtain learning methods for operators with unstructured stencils. In GMLS-Nets the necessary calculations are local, readily parallelizable, and the estimator is supported by a rigorous approximation theory. We show how the framework may be used for unstructured physical data sets to perform functional regression to identify associated differential operators and to regress quantities of interest. The results suggest the architectures to be an attractive foundation for data-driven model development in scientific machine learning applications

97 MATHEMATICS AND COMPUTING↗

Statistical estimates of the binary properties of rotational variables

ABSTRACT We present a model to estimate the average primary masses, companion mass ranges, the inclination limit for recognizing a rotational variable, and the primary mass spreads for populations of binary stars. The model fits a population’s binary mass function distribution and allows for a probability that some mass functions are incorrectly estimated. Using tests with synthetic data, we assess the model’s sensitivity to each parameter, finding that we are most sensitive to the average primary mass and the minimum companion mass, with less sensitivity to the inclination limit and little to no sensitivity to the primary mass spread. We apply the model to five populations of binary spotted rotational variables identified in ASAS-SN, computing their binary mass functions using RV data from APOGEE. Their average primary mass estimates are consistent with our expectations based on their CMD locations ($\sim 0.75 \, {\rm M}_{\odot }$ for lower main sequence primaries and $\sim 0.9$–$1.2 \, {\rm M}_{\odot }$ for RS CVn and sub-subgiants). Their companion mass range estimates allow companion masses down to $M_2/M_1\simeq 0.1$, although the main sequence population may have a higher minimum mass fraction ($\sim 0.4$). We see weak evidence of an inclination limit $\gtrsim 50^{\circ }$ for the main sequence and sub-subgiant groups and no evidence of an inclination limit in the other groups. No groups show strong evidence for a preferred primary mass spread. We conclude by demonstrating that the approach will provide significantly better estimates of the primary mass and the minimum mass ratio and reasonable sensitivity to the inclination limit with 10 times as many systems.

Phillips, Anya (ORCID:000900051914974X)↗

Nucleon momentum fraction, helicity and transversity from 2+1-flavor lattice QCD

A detailed analysis of the systematic uncertainties in the calculation of the isovector momentum fraction, ( x ) u - d , helicity moment, ( x ) Δ u - Δ d , and the transversity moment, ( x ) δu - δd , of the nucleon is presented using high-statistics data on seven ensembles of gauge configurations generated by the JLab/W&M/LANL/MIT collaborations using 2 + 1-flavors of dynamical Wilson-clover quarks. The much higher statistics have facilitated better control over all systematics compared to previous lattice calculations. The least understood systematic — excited-state contamination — is quantified by studying the variation of the results as a function of different estimates of the mass gap of the first excited state, obtained from two- and three-point correlation functions, and as a function of the pion mass M π . The final results are obtained using a simultaneous fit in the lattice spacing a , pion mass M π and the finite volume parameter M π L keeping leading order corrections. The data show no significant dependence on the lattice spacing and some evidence for finite-volume corrections. Our final results, in the $ \overline{\mathrm{MS}} $ scheme at 2 GeV, are ( x ) u - d = 0 . 155(17)(20), ( x ) Δ u - Δ d = 0 . 183(14)(20) and ( x ) δu - δd = 0 . 220(18)(20), where the first error is the overall analysis uncertainty assuming excited-state contributions have been removed, and the second is an additional systematic uncertainty due to possible residual excited-state contributions. These results are consistent with phenomenological global fit values.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

The Containment Assurance Risk Framework of the Mars Sample Return Program

The Mars Sample Return campaign aims at bringing rock and atmospheric samples from Mars to Earth through a series of robotic missions. These missions would collect the samples being cached and deposited on Martian soil by the Perseverance rover, place them in a container, and launch them into Martian orbit for subsequent capture by an orbiter that would bring them back. Given there exists a non-zero probability that the samples contain biological material, precautions are being taken to design systems that would break the chain of contact between Mars and Earth. These include techniques such as sterilization of Martian particles, redundant containment vessels, and a robust reentry capsule capable of accurate landings without a parachute. Requirements exist that the probability of containment not assured of Martian-contaminated material into Earth’s biosphere be less than one in a million. To demonstrate compliance with this strict requirement, a statistical framework was developed to assess the likelihood of containment loss during each sample return phase and make a statement about the total combined mission probability of containment not assured. The work presented here describes this framework, which considers failure modes or fault conditions that can initiate failure sequences ultimately leading to containment not assured. Reliability estimates are generated from databases, design heritage, component specifications, or expert opinion in the form of probability density functions or point estimates and provided as inputs to the mathematical models that simulate the different failure sequences. The probabilistic outputs are then combined following the logic of several fault trees to compute the ultimate probability of containment not assured. Given the multidisciplinary nature of the problem and the different types of mathematical models used, the statistical tools needed for analysis are required to be computationally efficient. While standard Monte Carlo approaches are used for fast models, a multi-fidelity approach to rare event probabilities is proposed for expensive models. In this paradigm, inexpensive low-fidelity models are developed for computational acceleration purposes while the expensive high-fidelity model is kept in the loop to retain accuracy in the results. This work presents an example of end-to-end application of this framework highlighting the computational benefits of a multi-fidelity approach.

Giuseppe Cataldo↗

Completely monotone regression estimates of software failure rates

A method for estimating the present failure rate of a program is presented. A crude nonparameter estimate of the failure rate function is obtained from past failure times. This estimate is then smoothed by fitting a completely monotonic function, which is the solution of a quadratic programming problem. The value of the smoothed function at present time is used as the estimate of present failure rate. Results of a Monte Carlo study of performance are given.

Miller, D. R.↗

A structure function representation theorem with applications to frequency stability estimation

Random processes with stationary nth differences serve as models for oscillator phase noise. The theorem proved here allows one to obtain the structure function (covariances of the nth differences) of such a process in terms of the differences of a single function of one time variable. In turn, this function can easily be obtained from the spectral density of the process. The theorem is used for computing the variance of two estimators of frequency stability.

Greenhall, C. A.↗

Approximation of Optimal Infinite Dimensional Compensators for Flexible Structures

The infinite dimensional compensator for a large class of flexible structures, modeled as distributed systems are discussed, as well as an approximation scheme for designing finite dimensional compensators to approximate the infinite dimensional compensator. The approximation scheme is applied to develop a compensator for a space antenna model based on wrap-rib antennas being built currently. While the present model has been simplified, it retains the salient features of rigid body modes and several distributed components of different characteristics. The control and estimator gains are represented by functional gains, which provide graphical representations of the control and estimator laws. These functional gains also indicate the convergence of the finite dimensional compensators and show which modes the optimal compensator ignores.

Gibson, J. S.↗

An Interactive MATLAB Program for Fitting Transfer Functions to Frequency Responses

A computer program called FRFit (Frequency Response Fitting) is described for fitting single-input single-output transfer function models to empirical frequency response data. The program is interactive in that the user specifies ``elementary factors'' (gain, delay, pure differentiators and integrators, and first- and second-order zeros and poles) by entering numerical values or moving sliders in a graphical user interface. A nonlinear optimization can then be performed to obtain maximum likelihood estimates of transfer function parameters and uncertainties to provide feedback on the modeling and refine estimates. Several examples are discussed, including the identification of aircraft pitch dynamics from simulation data and data reported in the literature, approximating Theodorsen's function of unsteady aerodynamics, and obtaining a reduced-order model of a computational fluid dynamics code describing the unsteady aerodynamics around an aeroelastic wing. The program has some usefulness as a teaching aid, and can be applied to model structure determination, reduced-order modeling, preliminary analysis, and simple system identification problems. The program was written in MATLAB and is planned for public release through the NASA Software Catalog.

System identification↗

Capturing Real-Time Power System Dynamics: Opportunities and Challenges

The power grid evolves towards a new mix of generation and consumption that introduces new dynamic and stochastic behaviors. These emerging grid behaviors would invalidate the steady-state assumption in today’s state estimation – an essential function for real-time power grid operation. This paper examines this steady-state assumption and identifies the need for estimating dynamic states. Supporting technologies are presented as well as a proposed formulation for estimating dynamic states. Metrics for evaluating methods for solving the dynamic state estimation problem are proposed, with example results to illustrate the use of these metrics. The overall objective of this paper is to provide a basis that more research on this topic can follow.

Huang, Zhenyu↗