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 1,099 records · Page 61

Viking Lander image analysis of Martian atmospheric dust

We have reanalyzed three sets of Viking Lander 1 and 2 (VL1 and VL2) images of the Martian atmosphere to better evaluate the radiative properties of the atmospheric dust particles. The properties of interest are the first two moments of the size distribution, the single-scattering albedo, the dust single-scattering phase function, and the imaginary index of refraction. These properties provide a good definition of the influence that the atmospheric dust has on heating of the atmosphere. Our analysis represents a significant improvement over past analyses (Pollack et al. 1977, 1979) by deriving more accurate brightness closer to the sun, by carrying out more precise analyses of the data to acquire the quantities of interest, and by using a better representation of scattering by nonspherical particles. The improvements allow us to better define the diffraction peak and hence the size distribution of the particles. For a lognormal particle size distribution, the first two moments of the size distribution, weighted by the geometric cross section, are found. The geometric cross-section weighted mean radius r(sub eff) is found to be 1.85 +/- 0.3 micrometers at VL2 during northern summer when dust loading was low and 1.52 +/- 0.3 micrometers at VL1 during the first dust storm. In both cases the best cross-section weighted mean variance nu(sub eff) of the size distribution is equal to 0.5 +/- 0.2 micrometers. The changes in size distribution, and thus radiative properties, do not represent a substantial change in solar energy deposition in the atmosphere over the Pollak et al. (1977, 1979) estimates.

Pollack, James B.↗

A Feasibility Study for Simultaneous Measurements of Water Vapor and Precipitation Parameters using a Three-frequency Radar

The radar return powers from a three-frequency radar, with center frequency at 22.235 GHz and upper and lower frequencies chosen with equal water vapor absorption coefficients, can be used to estimate water vapor density and parameters of the precipitation. A linear combination of differential measurements between the center and lower frequencies on one hand and the upper and lower frequencies on the other provide an estimate of differential water vapor absorption. The coupling between the precipitation and water vapor estimates is generally weak but increases with bandwidth and the amount of non-Rayleigh scattering of the hydrometeors. The coupling leads to biases in the estimates of water vapor absorption that are related primarily to the phase state and the median mass diameter of the hydrometeors. For a down-looking radar, path-averaged estimates of water vapor absorption are possible under rain-free as well as raining conditions by using the surface returns at the three frequencies. Simulations of the water vapor attenuation retrieval show that the largest source of error typically arises from the variance in the measured radar return powers. Although the error can be mitigated by a combination of a high pulse repetition frequency, pulse compression, and averaging in range and time, the radar receiver must be stable over the averaging period. For fractional bandwidths of 20% or less, the potential exists for simultaneous measurements at the three frequencies with a single antenna and transceiver, thereby significantly reducing the cost and mass of the system.

Meneghini, R.↗

Resource Management in the Microgravity Science Division

In the Microgravity Science Division, the primary responsibilities of the Business Management Office are resource management and data collection. Resource management involves working with a budget to do a number of specific projects, while data collection involves collecting information such as the status of projects and workforce hours. This summer in the Business Management Office I assisted Margie Allen with resource planning and the implementation of specific microgravity projects. One of the main duties of a Project Control Specialists, such as my mentor, is to monitor and analyze project manager s financial plans. Project managers work from the bottom up to determine how much money their project will cost. They then set up a twelve month operating plan which shows when money will be spent. I assisted my mentor in checking for variances in her data against those of the project managers. In order to successfully check for those variances, we had to understand: where the project is including plans vs. actual performance, why it is in its present condition, and what the future impact will be based on known budgetary parameters. Our objective was to make sure that the plan, or estimated resources input, are a valid reflection of the actual cost. To help with my understanding of the process, over the course of my tenure I had to obtain skills in Microsoft Excel and Microsoft Access.

Casselle, Justine↗

Uncertainty analysis of terrestrial net primary productivity and net biome productivity in China during 1901-2005

Here, despite the importance of net primary productivity (NPP) and net biome productivity (NBP), estimates of NPP and NBP for China are highly uncertain. To investigate the main sources of uncertainty, we synthesized model estimates of NPP and NBP for China from published literature and the Multi-scale Synthesis and Terrestrial Model Intercomparison Project (MsTMIP). The literature-based results showed that total NPP and NBP in China were 3.35 ± 1.25 and 0.14 ± 0.094 Pg C yr –1 , respectively. Classification and regression tree analysis based on literature data showed that model type was the primary source of the uncertainty, explaining 36% and 64% of the variance in NPP and NBP, respectively. Spatiotemporal scales, land cover conditions, inclusion of the N cycle, and effects of N addition also contributed to the overall uncertainty. Results based on the MsTMIP data suggested that model structures were overwhelmingly important (>90%) for the overall uncertainty compared to simulations with different combinations of time-varying global change factors. The interannual pattern of NPP was similar among diverse studies and increased by 0.012 Pg C yr –1 during 1981–2000. In addition, high uncertainty in China's NPP occurred in areas with high productivity, whereas NBP showed the opposite pattern. Our results suggest that to significantly reduce uncertainty in estimated NPP and NBP, model structures should be substantially tested on the basis of empirical results. To this end, coordinated distributed experiments with multiple global change factors might be a practical approach that can validate specific structures of different models.

54 ENVIRONMENTAL SCIENCES↗

Validation of SAS4A/SASSYS-1 for predicting steady-state single-phase natural circulation

The validation of system analysis codes for nuclear reactor systems is required for the development and application of these computational tools. Designed as a comprehensive system analysis code for advanced nuclear reactors, SAS4A/SASSYS-1 requires validation of its physics model for capturing single-phase natural circulation behavior. To support the validation of SAS4A/SASSYS-1, high-precision experiments are performed capturing steady-state single-phase natural circulation on a scaled facility with comprehensive instrumentation. In this work, dedicated tests are performed quantifying the critical modeling parameters, and a single-phase natural circulation benchmark dataset is obtained with well-documented uncertainty and comprehensive facility description. The validation is then performed against the dataset examining the capability of SAS4A/SASSYS-1 in simulating steady-state single-phase natural circulation. The experimental facility is modeled in the candidate code. Solution verification is performed using Richardson-extrapolation-based estimators which quantify and restrict numerical errors from discretization. Input uncertainty provided by the benchmark dataset is forward propagated through the candidate code, quantifying the output uncertainty in a Monte Carlo approach. The composition of the output uncertainty is also quantified through a variance-based sensitivity analysis. With the uncertainty quantified for each individual condition, a detailed comparison between the simulation results and experimental data is performed covering the whole dataset. The results show consistent agreement for all primary parameters. The current validation activity provides a valuable benchmark dataset for the validation of system analysis codes in capturing single-phase natural circulation and demonstrates satisfactory prediction capability of SAS4A/SASSYS1 for steady-state single-phase natural circulation.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Optimal multifrequency weighting for CMB lensing

We report that extragalactic foregrounds in cosmic microwave background (CMB) temperature maps lead to significant biases in CMB lensing reconstruction if not properly accounted for. Combinations of multifrequency data have been used to minimize the overall map variance (internal linear combination, or ILC), or specifically null a given foreground, but these are not tailored to CMB lensing. In this paper, we derive an optimal multifrequency combination to jointly minimize CMB lensing noise and bias. We focus on the standard lensing quadratic estimator, as well as the "shear-only"and source-hardened estimators, whose responses to foregrounds differ. We show that an optimal multifrequency combination is a compromise between the ILC and joint deprojection, which nulls the thermal Sunyaev-Zel'dovich (tSZ) and cosmic infrared background (CIB) contributions. In particular, for a Simons Observatory-like experiment with ℓ max,T =3000, we find that profile hardening alone (with the standard ILC) reduces the bias to the lensing power amplitude by 40%, at a 20% cost in noise, while the bias to the cross-correlation with a LSST-like sample is reduced by nearly an order of magnitude at a 10% noise cost, relative to the standard quadratic estimator. With a small amount of joint deprojection the bias to the profile hardened estimator can be further reduced to less than half the statistical uncertainty on the respective amplitudes, at a 20% and 5% noise cost for the auto- and cross-correlation respectively, relative to the profile hardened estimator with the standard ILC weights. Finally, we explore possible improvements with more aggressive masking and varying ℓ max,T .

79 ASTRONOMY AND ASTROPHYSICS↗

ADVANTG An Automated Variance Reduction Parameter Generator

The AutomateD VAriaNce reducTion Generator (ADVANTG) software automates the generation of variance reduction (VR) parameters for continuous-energy Monte Carlo simulations of fixed-source neutron, photon, and coupled neutron-photon transport problems using MCNP5 (X-5 Monte Carlo Team 2003). ADVANTG generates space- and energy dependent mesh-based weight-window bounds and biased source distributions from three dimensional (3-D) discrete ordinates (SN) calculations that are performed by the Denovo package (Evans et al. 2010). The deterministic calculations can be performed using multiple cores and/or processors (e.g., on multi-core desktop systems and clusters). The final variance reduction parameters are output in a format that can be used with unmodified versions of MCNP. The primary objective of the development of ADVANTG has been to reduce both the user effort and the computational time required to obtain accurate and precise tally estimates across a broad range of challenging transport application areas.

ADVANTG↗

Mean sea surface height of the Antarctic Circumpolar Current from Geosat data: Method and application

The mean sea surface height across the Antarctic Circumpolar Current (ACC) has been reconstructed from height variability measured by the Geosat altimeter without assuming prior knowledge of the geoid. For this study, an automated technique has been developed to estimate mean sea surface height for each satellite ground track using a meandering Gaussian jet model, and errors have been estimated using Monte Carlo simulation. The results are objectively mapped to produce a picture of the mean Subantarctic and Polar Fronts, which together comprise the major components of the Antarctic Circumpolar Current. The meandering jet model explains between 40% and 70% of the height variance along the jet axes. The results show that the fronts are substantially steered by topography and that the jets have an average Gaussian width of about 44 km in the meridional direction and meander about 75 km to either side of their mean locations. The average height difference across the Subantarctic Front (SAF) is 0.7 m and across the Polar Front (PF) 0.6 m. The mean widths of the fronts are correlated with the size of the baroclinic Rossby radius.

Gille, Sarah T.↗

The Influence of Observation Errors on Analysis Error and Forecast Skill Investigated with an Observing System Simulation Experiment

The Global Modeling and Assimilation Office (GMAO) observing system simulation experiment (OSSE) framework is used to explore the response of analysis error and forecast skill to observation quality. In an OSSE, synthetic observations may be created that have much smaller error than real observations, and precisely quantified error may be applied to these synthetic observations. Three experiments are performed in which synthetic observations with magnitudes of applied observation error that vary from zero to twice the estimated realistic error are ingested into the Goddard Earth Observing System Model (GEOS-5) with Gridpoint Statistical Interpolation (GSI) data assimilation for a one-month period representing July. The analysis increment and observation innovation are strongly impacted by observation error, with much larger variances for increased observation error. The analysis quality is degraded by increased observation error, but the change in root-mean-square error of the analysis state is small relative to the total analysis error. Surprisingly, in the 120 hour forecast increased observation error only yields a slight decline in forecast skill in the extratropics, and no discernable degradation of forecast skill in the tropics.

Observation Error↗

Estimation of distributions via multilevel Monte Carlo with stratified sampling

We design and implement a novel algorithm for computing a multilevel Monte Carlo (MLMC) estimator of the joint cumulative distribution function (CDF) of a vector-valued quantity of interest in problems with random input parameters and initial conditions. Our approach combines MLMC with stratified sampling of the input sample space by replacing standard Monte Carlo at each level with stratified Monte Carlo initialized with proportionally allocated samples. We show that the resulting stratified MLMC (sMLMC) algorithm is more efficient than its standard MLMC counterpart due to the additional variance reduction provided by the stratification of the random parameter's domain, especially at the coarsest levels. Additional computational cost savings are obtained by smoothing the indicator function with a Gaussian kernel, which proves to be an efficient and robust alternative to recently developed polynomial-based techniques.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Interannual Sea Level Variations in the Tropical Indian Ocean from Geosat and Shallow Water Simulations

Sea level variations of the Indian Ocean north of 20 deg S are analyzed from Geosat satellite altimeter data over April 1985-September 1989. These variations are compared and interpreted with numerical simulations derived from a reduced gravity model forced by FSU observed winds over the same period. After decomposition into complex empirical orthogonal functions, the low-frequency anomalies are described by the first two modes for observations as well as for simulations. The sums of the two modes contain 34% and 40% of the observed and simulated variances, respectively. Averaged over the basin, the observed and simulated sea level changes are correlated by 0.92 over 1985-1988. The strongest change happens during the El Ninio 1986-1987: between winter 1986 and summer 1987 the basin-averaged sea level rises by approx. 1 cm. These low-frequency variations can partly be explained by changes in the Sverdrup circulation. The southern tropical Indian Ocean between 1O deg and 20 deg S is the domain where those changes are strongest: the averaged sea level rises by approx. 4.5 cm between winter 1986 and winter 1987. There, the signal propagates southwestward across the basin at a speed similar to free Rossby waves. Sensitivity of observed anomalies is examined over 1987-1988, with different orbit ephemeris, tropospheric corrections, and error reduction processes. The uncertainty of the basin-averaged sea level estimates is mostly due to the way the orbit error is reduced and reaches approx. 1 cm. Nonetheless, spatial correlation is good between the various observations and better than between observations and simulations. Sensitivity of simulated anomalies to the wind uncertainty, examined with Former Soviet Union (FSU) and European Center for Medium-Range Weather Forecasting (ECMWF) forcings over 1985-1988, shows that the variance of the simulations driven by ECMWF is 52% smaller, as FSU winds are stronger than ECMWF. Results show that the wind strength also affects the dynamic response of the ocean: anomalies propagate westward across the basin more than twice as fast with FSU than with ECMWF. It is found that the discrepancy is larger between ECMWF and FSU simulations than between observations and FSU simulations.

Perigaud, Claire↗

On estimating the effects of clock instability with flicker noise characteristics

A scheme for flicker noise generation is given. The second approach is that of successive segmentation: A clock fluctuation is represented by 2N piecewise linear segments and then converted into a summation of N+1 triangular pulse train functions. The statistics of the clock instability are then formulated in terms of two sample variances at N+1 specified averaging times. The summation converges very rapidly that a value of N 6 is seldom necessary. An application to radio interferometric geodesy shows excellent agreement between the two approaches. Limitations to and the relative merits of the two approaches are discussed.

Wu, S. C.↗

Biomass Retrieval from L-Band Polarimetric UAVSAR Backscatter and PRISM Stereo Imagery

The forest above-ground biomass (AGB) and spatial distribution of vegetation elements have profound effects on the productivity and biodiversity of terrestrial ecosystems. In this paper, we evaluated biomass estimation from L-band Synthetic Aperture Radar (SAR) data acquired by National Aeronautics and Space Administration (NASA) Uninhabited Aerial Vehicle SAR (UAVSAR) and the improvement of accuracy by adding canopy height information derived from stereo imagery acquired by Japan Aerospace Exploration Agency (JAXA) Panchromatic Remote Sensing Instrument for Stereo Mapping (PRISM) on-board the Advanced Land Observing Satellite (ALOS). Various models for prediction of forest biomass from UAVSAR data were investigated at pixel sizes of 1/4 ha (50 m x 50 m) and 1 ha. The variance inflation factor (VIF) was calculated for each of the explanatory variables in multivariable regression models to assess the multi-collinearity between explanatory variables. In addition, the t-and p-values were used to interpret the significance of the coefficients of each explanatory variables. The R(exp. 2), Root Mean Square Error (RMSE), bias and Akaike information criterion (AIC), and leave-one-out cross-validation (LOOCV) and bootstrapping were used to validate models. At 1/4-ha scale, the R(exp. 2) and RMSE of biomass estimation from a model using a single track of polarimetric UAVSAR data were 0.59 and 52.08 Mg/ha. With canopy height from PRISM as additional independent variable, R(exp. 2) increased to 0.76 and RMSE decreased to 39.74 Mg/ha (28.24%). At 1-ha scale, the RMSE of biomass estimation based on UAVSAR data of a single track was 39.42 Mg/ha with a R(exp. 2) of 0.77. With the canopy height from PRISM, R(exp. 2) increased to 0.86 and RMSE decreased to 29.47 Mg/ha (20.18%). The models using UAVSAR data alone underestimated biomass at levels above approximately 150 Mg/ha showing the saturation phenomenon. Adding canopy height from PRISM stereo imagery significantly improved the biomass estimation and elevated the saturation level in estimating biomass. Combined use of UAVSAR data acquired from opposite directions (odd and even tracks) slightly improved the biomass estimation.Combined use of UAVSAR data acquired from opposite directions (odd and even tracks) slightly improved the biomass estimation at 1/4-ha scale, R(exp. 2) increased from 0.59 to 0.66 and RMSE reduced from 52.08 to 48.57 Mg/ha. Averaging multiple acquisitions of UAVSAR data from the same look azimuth direction did not improve biomass estimation. A biomass map derived from NASA's LVIS (Laser Vegetation Imaging System) wave-form data was used as a reference for evaluation of the biomass maps from these models. The study has also shown that the errors decreased when deciduous, evergreen, and mixed forests were modeled separately but the improvement was not significant

Zhang, Zhiyu↗

On the scaling behavior of hardness with ligament diameter of nanoporous-Au: Constrained motion of dislocations along the ligaments

The scaling behavior of hardness with ligament diameter and vacancy defect concentration in nanoporous Au (np-Au) has been investigated using a combination of Vickers Hardness, Scanning electron microscopy, and positron lifetime measurements. It is shown that for np-Au, the hardness scales with the ligament diameter with an exponent of −0.3, that is, at variance with the conventional Hall-Petch exponent of −0.5 for bulk systems, as seen in the controlled experiments on cold worked Au with varying grain size. The hardness of np-Au correlates with the vacancy concentration C{sub V} within the ligaments, as estimated from positron lifetime experiments, and scales as C{sub V}{sup 1/2}, pointing to the interaction of dislocations with vacancies. The distinctive Hall-Petch exponent of −0.3 seen for np-Au, with ligament diameters in the range of 5–150 nm, is rationalized by invoking the constrained motion of dislocations along the ligaments.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

Choosing a reliability inspection plan for interval censored data

Reliability test plans are important for producing precise and accurate assessment of reliability characteristics. This paper explores different strategies for choosing between possible inspection plans for interval censored data given a fixed testing timeframe and budget. A new general cost structure is proposed for guiding precise quantification of total cost in inspection test plan. Multiple summaries of reliability are considered and compared as the criteria for choosing the best plans using an easily adapted method. Different cost structures and representative true underlying reliability curves demonstrate how to assess different strategies given the logistical constraints and nature of the problem. Results show several general patterns exist across a wide variety of scenarios. Given the fixed total cost, plans that inspect more units with less frequency based on equally spaced time points are favored due to the ease of implementation and consistent good performance across a large number of case study scenarios. Plans with inspection times chosen based on equally spaced probabilities offer improved reliability estimates for the shape of the distribution, mean lifetime, and failure time for a small fraction of population only for applications with high infant mortality rates. The paper uses a Monte Carlo simulation based approach in addition to the common evaluation based on the asymptotic variance and offers comparison and recommendation for different applications with different objectives. Additionally, the paper outlines a variety of different reliability metrics to use as criteria for optimization, presents a general method for evaluating different alternatives, as well as provides case study results for different common scenarios.

97 MATHEMATICS AND COMPUTING↗

Allometric and Mobile Terrestrial LiDAR Modeling of Aboveground Woody Biomass of Populus in Coppice Production

Poplars ( Populus spp.) and their hybrids are increasingly being grown in coppice production to generate bioenergy feedstocks at frequent intervals. Allometric equations are re-quired to predict aboveground biomass (AGB) of coppiced individuals with minimal field measurements. Likewise, remote sensing tools like LiDAR (light detection and ranging) can be used if models are available to predict AGB from point cloud data. Therefore, this study sought to develop equations to predict dry woody AGB from field measurements and LiDAR data from coppiced poplar field trials containing eastern cottonwood ( P. del-toides ) and hybrid poplar taxa. We found that taxa-specific allometric models containing the summed basal area of the three largest stems in the coppice provided the best predictive model, with stem height and stem count failing to provide additional explanatory power. The best predictive LiDAR-based model was independent of taxa but had slightly lower adjusted R 2 and higher RMSE than the allometric model. It contained four parameters including crown volume, leaf area index, variance of height returns, and the top point density (i.e., density metric 9 or the proportion of points in the highest point interval when the point cloud is evenly divided into ten vertical intervals). In total, these models can be used to quickly and efficiently estimate dry woody AGB of Populus coppice systems for bioenergy feedstock production.

AGB↗

Efficient Reliability Analysis using Generalized Multifidelity Modeling and Explainable Active Learning

To assess the reliability of critical technologies like nuclear plants and infrastructure systems and improve the robustness of design, engineers have to quantify the uncertainties surrounding the system behavior accurately. However, the complexity of the problem can make standard reliability analysis algorithms prohibitively expensive, primarily due to the high computational cost of estimating the system response at each iteration. This cost can be greatly reduced by using multi-fidelity modeling and machine learning to build a surrogate model to replace the expensive response function. We propose a general and robust method for building surrogates from multiple Low Fidelity (LF) models coupled with machine learning to retain accuracy. Our framework first constructs “Corrected Low Fidelity models” (CLFs) by coupling a High Fidelity (HF) model inferred Gaussian Process correction term with each of the LF models. It then uses the correction terms to assign model probabilities to each of these CLFs in an explainable way before using them to assemble the final surrogate. No assumptions are made about the type of the LF models or their correlation with the HF model. The proposed surrogate modeling framework is used within the subset simulation algorithm (a variance-reduced MCMC-based reliability analysis algorithm) for enhanced efficiency. Additionally, an active learning step is added to the algorithm to adaptively decide when the surrogate is not sufficiently accurate, at which point the HF model is called and used to refine the surrogate. Through a frame buckling example, our method is shown to be highly efficient at reducing the expensive HF model calls while accurately estimating the failure probability.

97 MATHEMATICS AND COMPUTING↗

Global variation of transpiration and soil evaporation and the role of their major climate drivers

Although global variation in actual evapotranspiration has been widely investigated, it remains unclear how its two major components, transpiration and soil evaporation, are driven by climate drivers across global land surface. This paper uses a well-validated, process-based model that estimates transpiration and soil evaporation, and for the first time investigates and quantifies how the main global drivers, associated to vegetation process and the water and energy cycle, drive the spatiotemporal variation of the two components. The results show that transpiration and soil evaporation dominate the variance of actual evapotranspiration in wet and dry regions, respectively. Dry southern hemisphere from 13°S to 27°S is highlighted since it contributes to 21% global soil evaporation variance, with only 11% global land area. In wet regions, particularly in the humid tropics, there are strong correlations between transpiration, actual evapotranspiration, and potential evapotranspiration, with precipitation playing a relatively minor role, and available radiative energy is the major contributor to the interannual variability in transpiration and actual evapotranspiration in Amazonia. Conversely in dry regions, there are strong correlations between soil evaporation, actual evapotranspiration, and precipitation. Lastly, our findings highlight that ecohydrological links are highly related to climate regimes, and the small region such as Australia has important contribution to interannual variation in global soil evaporation and evapotranspiration, and anthropogenic activities strongly influence the variances in irrigation regions

54 ENVIRONMENTAL SCIENCES↗