Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “variance 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 343 records · Page 19

In Situ Inference for Earth System Predictability

Focal Area: Focal Area 3: Insight gleaned from complex simulated data using AI, big data analytics, and other advanced methods, including explainable AI and physics- or knowledge-guided AI. Science Challenge: An understanding of future evolution in precipitation extremes is critical to numerous DOE mission questions. Extreme events are by nature short time-scale events that are difficult to diagnose in available model data. Accurate modeling of extreme events necessarily requires high spatial resolution at the storm scale locally. However, the environment in which storms grow is dependent on global, remote, processes. These complex spatiotemporal relationships are impossible to diagnose at resolutions required to accurately model storms responsible for extreme precipitation. At exascale, climate simulations will produce results at fine enough resolution to investigate these relationships. However, the resulting data from these simulations will be far too large to save for post-simulation analysis. We advocate for fitting statistical models inside the simulations as they run, a context known as in situ, which will facilitate scientific investigations using the full fine-scale data stream. Figure 1 shows an example of the type of model we could consider, a Bayesian hierarchical spatial regression model. Precipitation extremes at each grid cell are modeled using extreme value distributions. Since extremes are rare, fitting models to individual grid cells can result in high variance and poor estimates. Instead, the model can be made more robust by smoothing the parameters of the extreme value model across space. Additionally, the parameters themselves can be functionally linked to other variables elsewhere in the simulation. Thus, we can use the fine-scale data to build more robust models for extremes that link extreme behavior to other climate patterns.

54 ENVIRONMENTAL SCIENCES↗

Improving the Nomenclature Around Uncertainty [Slides]

This presentation finds that defining the marginal probability density function (PDF) for nuclear data is important. Additionally, the vocabulary of “means and covariances” and new GNDS 2.0 formats are limited to Gaussian (normal) representations— always incorrect—but clearly of practical significance when uncertainties are large (>40%). Finally, the Triage Solution: declare our current data as containing best estimate (mode) plus variance for a truncated normal or lognormal.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Advances in the Assessment of Wind Turbine Operating Extreme Loads via More Efficient Calculation Approaches

A new adaptive stratified importance sampling (ASIS) method is proposed as an alternative approach for the calculation of the 50 year extreme load under operational conditions, as in design load case 1.1 of the the International Electrotechnical Commission design standard. ASIS combines elements of the binning and extrapolation technique, currently described by the standard, and of the importance sampling (IS) method to estimate load probability of exceedances (POEs). Whereas a Monte Carlo (MC) approach would lead to the sought level of POE with a daunting number of simulations, IS-based techniques are promising as they target the sampling of the input parameters on the parts of the distributions that are most responsible for the extreme loads, thus reducing the number of runs required. We compared the various methods on select load channels as output from FAST, an aero-hydro-servo-elastic tool for the design and analysis of wind turbines developed by the National Renewable Energy Laboratory (NREL). Our newly devised method, although still in its infancy in terms of tuning of the subparameters, is comparable to the others in terms of load estimation and its variance versus computational cost, and offers great promise going forward due to the incorporation of adaptivity into the already powerful importance sampling concept.

17 WIND ENERGY↗

The geoid spectrum from altimetry

Satellite altimetry information from the world's major oceans was analyzed to arrive at a geoid power spectrum. Using the equivalent of about 7 revolutions of data (mostly from GEOS-3) the power spectrum of the sea surface generally follows the expected values from Kaula's rule applied to the geoid. Analysis of overlapping altimetry arcs (and oceanographic data) shows that the surface spectrum is dominated by the geoid to about 500 cycles (40 km half wavelength) but that sea state departures are significant starting at about 250 cycles (80 km). Estimates of geopotential variances from a derived (smooth) geoid spectrum show significantly less power than Kaula's rule to about 60 cycles, but somewhat more from there to about 400 cycles. At less than 40 km half wavelength, the total power in the marine geoid may be negligible.

Wagner, C. A.↗

Optimization of fringe-type laser anemometers for turbine engine component testing

The fringe type laser anemometer is analyzed using the Cramer-Rao bound for the variance of the estimate of the Doppler frequency as a figure of merit. Mie scattering theory is used to calculate the Doppler signal wherein both the amplitude and phase of the scattered light are taken into account. The noise from wall scatter is calculated using the wall bidirectional reflectivity and the irradiance of the incident beams. A procedure is described to determine the optimum aperture mask for the probe volume located a given distance from a wall. The expected performance of counter type processors is also discussed in relation to the Cramer-Rao bound. Numerical examples are presented for a coaxial backscatter anemometer.

Seasholtz, R. G.↗

Remote sensing of floe size distribution and surface topography

Floe size can be measured by several properties p- for instance, area or mean caliper diameter. Two definitions of floe size distribution seem particularly useful. F(p), the fraction of area covered by floes no smaller than p; and N(p), the number of floes per unit area no smaller than p. Several summertime distributions measured are a graph, their slopes range from -1.7 to -2.5. The variance of an estimate is also calculated.

Rothrock, D. A.↗

Optimization of fringe-type laser anemometers for turbine engine component testing

The fringe type laser anemometer is analyzed using the Cramer-Rao bound for the variance of the estimate of the Doppler frequency as a figure of merit. Mie scattering theory is used to calculate the Doppler signal wherein both the amplitude and phase of the scattered light are taken into account. The noise from wall scatter is calculated using the wall bidirectional reflectivity and the irradiance of the incident beams. A procedure is described to determine the optimum aperture mask for the probe volume located a given distance from a wall. The expected performance of counter type processors is also discussed in relation to the Cramer-Rao bound. Numerical examples are presented for a coaxial backscatter anemometer. Previously announced in STAR as N84-25019

Seasholtz, R. G.↗

Rendezvous BET Program

Computes relative positions of two vehicles in concentric orbits. LRBET3 program best-estimate-of-trajectory (BET) calculation for postflight trajectory analysis of Shuttle orbital rendezvous maneuvers. LRBET3 produces estimated measurements for reconstructing relative positions of two vehicles. Kalman filter and smoothing filter applied to relative measurement input data to estimate state vector, reduce noise, and produce BET output. BET calculation minimizes variances of all trajectory estimation errors. LRBET3 written in FORTRAN IV for batch execution.

Lear, W. M.↗

Real-time flutter identification with close mode resolution

Real-time flutter prediction including close modes can be effectively estimated from turbulence or on-board excitation with an Extended Kalman Filter (EKF) approach. A physically based model form enables prediction of the damping rate as well as damping, giving a time to instability estimate with its variance. The approach is recursive and can operate asynchronously to drop data outliers and hence is quite robust. Its speed is reasonable for on-line application but can also be used effectively as an off-line analysis tool for application to any modal testing situation.

Roy, R. H.↗

Observations of interstellar absorption lines towards hot white dwarfs

A method of detecting and determining parameters for narrow absorption lines in IUE high-resolution spectra of hot white dwarfs was developed. It provides a statistical estimate of the variance of the parameters measured. The statistical accuracy of narrow line measurements with IUE spectra is discussed. The method was employed to determine the column densities of low-ionization species (including N I, Si II, and C II) in the line of sight to seven hot white dwarfs. High-ionization species (N V, Si IV, and C IV) are detected towards four of the stars. The temperature of the low-ionization gas towards 4 of the stars is constrained to be less than 50,000 K. The temperature of the high-ionization gas to two of the stars is constrained to be under 100,000 K, which is well below the equilibrium ionization temperature of a collisionally ionized gas. This suggests that this gas is either in a photoionized region (most likely near the star) or is produced in the photosphere of the star.

Jelinsky, Patrick↗

Third generation earth radiation budget measurements; ERBE in the context of earlier systems

The Earth Radiation Budget Experiment (ERBE) observations are just becoming available for scientific use. These represent the third generation of measurements with steadily improving accuracy and resolution. Beginning in the 1960's observations by spherical detectors established the mean albedo of the Earth near 30 percent in substantial variance from presatellite estimates. The Nimbus 6 and 7 wide field of view ERB measurements represent a long-term climatology of measurements at 1000 km resolution. The ERBE measurements introduce higher accuracy and higher space and time resolution result. Comparisons will be presented of several April ERB measurements to illustrate what this improvement in resolution and accuracy can yield. Simultaneous ERBE and Nimbus 7 measurements for April 1985 show nearly identical results on the large scale. Comparison of measurements of direct solar energy from ERBE, Solar Max Mission and Nimbus 7 suggest a 'solar constant' value of 1368 w/sq m for the 1979 to 1986 period. The long-term record of earth radiation budget (Aprils from 1976 to 1985) over large regions is shown to have interannual variation of plus or minus 20 to 30 w/sq. m. The new ERBE data will allow this climate record measurement to continue.

Vonderhaar, Thomas H.↗

Toward a More Robust Pruning Procedure for MLP Networks

Choosing a proper neural network architecture is a problem of great practical importance. Smaller models mean not only simpler designs but also lower variance for parameter estimation and network prediction. The widespread utilization of neural networks in modeling highlights an issue in human factors. The procedure of building neural models should find an appropriate level of model complexity in a more or less automatic fashion to make it less prone to human subjectivity. In this paper we present a Singular Value Decomposition based node elimination technique and enhanced implementation of the Optimal Brain Surgeon algorithm. Combining both methods creates a powerful pruning engine that can be used for tuning feedforward connectionist models. The performance of the proposed method is demonstrated by adjusting the structure of a multi-input multi-output model used to calibrate a six-component wind tunnel strain gage.

Stepniewski, Slawomir W.↗

Information Content of the Near-Field I: Two-Dimensional Samples

Limits on the effective resolution of many optical near-field experiments are investigated. The results are applicable to variants of total-internal-reflection microscopy (TIRM), photon-scanning-tunneling microscopy (PSTM), and near-field-scanning-optical microscopy (NSOM) in which the sample is weakly scattering and the direction of illumination may be controlled. Analytical expressions for the variance of the estimate of the complex susceptibility of an unknown two-dimensional object as a function of spatial frequency are obtained for Gaussian and Poisson noise models, and a model-independent measure is examined. The results are used to explore the transition from near-zone to far-zone detection. It is demonstrated that the information content of the measurements made at a distance of even one wavelength away from the sample is already not much different from the information content of the far field. Copyright 2004 Optical Society of America

Frazin, Richard A.↗

Algorithm Performance Dataset from NASA Open-Source Software

NASA Langley Research Center has recently developed and released the open-source software Multi Model Monte Carlo with Python (MXMCPy- LAR-19756-1) as a general capability for computing the statistics of outputs from an expensive, high-fidelity model by leveraging faster, low-fidelity models for speedup. Given a fixed computational budget and a collection of models with varying cost/accuracy, multi model Monte Carlo (MC) seeks a sample allocation strategy across the models that results in an estimator with optimal variance reduction. MXMCPy is a versatile tool that enables convenient access to many existing multi-model MC approaches (over a dozen algorithms available) within one modular and extensible package [1]. With MXMCPy, users can easily compare existing methods to determine the best choice for their particular problem,while developers have a basis for implementing and sharing new variance reduction approaches. However,there is currently very little understanding about which algorithm will perform best for a given problem (defined by the correlation between and relative cost of the available models) without a brute force search.

Geoffrey F Bomarito↗

Assessing Risk Due to Small Sample Size in Probability of Detection Analysis Using Tolerance Intervals

Small sample size (e.g.6-30) poses risk in results of probability of detection (POD) analysis using tolerance intervals. This method is also called as the limited sample or LS POD. The analysis is performed either during NDE procedure qualification or for assessment of reliability of an NDE procedure. The risk is primarily due to sampling error. Smaller samples are not likely to be random to the population or representative of the population. The small samples are likely to be biased. Biased samples have smaller standard deviation compared to the population. POD analysis with small biased sample can lead to overestimation of POD. Many sampling schemes are available in statistics to mitigate sampling risk. Primary objective of POD analysis is to determine a decision threshold from signal response measurements of a sample such that it is less than or equal to population decision threshold for 90% POD. Sampling error implies that this NDE reliability condition is violated. One of sampling types is called a representative sample. Representative samples reduce variance in POD estimates but also reduce magnitude of the error. Sampling sensitivity analysis for some sampling types is performed here using repetitive random sampling or Monte Carlo method. Six sampling types are considered for comparison. Some of the sampling types are similar to drawing a representative sample. LS POD model assumes random sampling. Therefore, random sampling is used as a basis for comparison with each sampling type. The sampling types used in the analysis are, A. Nominal and worst-case sampling, B. Worst-case sampling, C. Nominal case sampling, D. Random sampling, E. Random target, and sub-target sampling. F. Nominal target and sub-target sampling. Results of Monte Carlo simulation indicate that type F sampling can mitigate sampling risk and is also more practical to implement. Type A sampling may also mitigate the sampling risk, but it may be less practical to implement.

Ajay M Koshti↗

SPT-3G D1: Foreground-Robust Lensing Templates for Primordial Gravitational Wave Searches

Gravitational lensing of the cosmic microwave background (CMB) generates B-mode polarization that acts as a source of contamination to searches for B modes generated by primordial gravitational waves (PGWs). The strongest constraint on PGW B modes is already significantly limited by lensing B modes, as shown in the most recent BICEP result. In this work, we present CMB lensing B-mode templates constructed using SPT-3G and Planck data, which characterize the lensing B modes and can be used to improve PGW B-mode searches. We use SPT-3G data from the 2019 and 2020 observing seasons for the E modes and the CMB-reconstructed lensing potential, and a cosmic infrared background (CIB) map from Planck as an external lensing tracer. To test for extragalactic foreground biases in the lensing template, we consider CMB lensing reconstruction variants with different levels of foreground immunity: the standard and profile-hardened global minimum variance (GMV) quadratic estimators, and a polarization-only quadratic estimator. We validate the template construction using Gaussian simulations and Agora simulations with realistic non-Gaussian foregrounds. From simulations, we find that foreground-induced biases are strongly suppressed for the template constructed with the profile-hardened GMV + CIB tracer, with residual bias below 10% of the statistical uncertainty on the template power spectrum. Data difference tests on this template similarly show no evidence for significant foreground contamination. This foreground-immune lensing template achieves delensed residual BB power of $A_{\rm lens}^{\rm res} \simeq 0.48$ averaged over $20 \leq \ell \leq 200$, the highest delensing efficiency lensing template to date. These results demonstrate and validate a method to construct foreground-robust lensing templates which will be used in upcoming delensed PGW B-mode analyses of BICEP data.

Nakato, Y. [Stanford U.]↗

Convergence Study in Global Sensitivity Analysis

Monte Carlo (MC) sampling is a common method used to randomly sample a range of scenarios. The associated error follows a predictable rate of convergence of $1/\sqrt{N}$, such that quadrupling the sample size halves the error. This method is often employed in performing global sensitivity analysis which computes sensitivity indices, measuring fractional contributions of uncertain model inputs to the total output variance. In this study, several models are used to observe the rate of decay in the MC error in the estimation of the conditional variance, the total variance in the output, and the global sensitivity indices. The purpose is to examine the rate of convergence of the error in existing specialized, albeit MC-based, sampling methods for estimation of the sensitivity indices. It was found that the conditional variances and sensitivity indices all follow the $1/\sqrt{N}$ convergence rate. Future work will test the convergence of observables from more complex models such as ignition time in combustion.

97 MATHEMATICS AND COMPUTING↗

Multigrid deflation for Lattice QCD

Computing the trace of the inverse of large matrices is typically addressed through statistical methods. Deflating out the lowest eigenvectors or singular vectors of the matrix reduces the variance of the trace estimator. This work summarizes our efforts to reduce the computational cost of computing the deflation space while achieving the desired variance reduction for Lattice QCD applications. Previous efforts computed the lower part of the singular spectrum of the Dirac operator by using an eigensolver preconditioned with a multigrid linear system solver. Despite the improvement in performance in those applications, as the problem size grows the runtime and storage demands of this approach will eventually dominate the stochastic estimation part of the computation. In this work, we propose to compute the deflation space in one of the following two ways. First, by using an inexact eigensolver on the Hermitian, but maximally indefinite, operator. Second, by exploiting the fact that the multigrid prolongator for this operator is rich in components toward the lower part of the singular spectrum. We show experimentally that the inexact eigensolver can approximate the lower part of the spectrum even for ill-conditioned operators. Also, the deflation based on the multigrid prolongator is more efficient to compute and apply, and, despite its limited ability to approximate the fine level spectrum, it obtains similar variance reduction on the trace estimator as deflating with approximate eigenvectors from the fine level operator.

97 MATHEMATICS AND COMPUTING↗