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,279 records · Page 71

Solar total irradiance and sunspot area in 1981

Precise observations of the total solar irradiance in 1980 from the solar-maximum mission, showed a strong correlation with the projected total sunspot area. This correlation explained about half of the variance (r of about 0.75) in the total-irradiance data, leaving the other half for other solar phenomena and errors in the data. The analysis of the intervening three years of spin-mode data, which have reduced coverage, has begun. It is found that the correlation persists at about the same qualitative level, but with increased scatter that can be attributed to the smaller amount of data. The flatness of the distribution of areas of sunspot groups makes it possible to estimate PSI approximately from only the large groups.

Hudson, H. S.↗

Evaluating importance maps for TRIPOLI-4{sup R} using deterministic or on-line methods - 25446

Variance reduction is a key ingredient for solving radiation-protection problems with Monte-Carlo particle-transport codes. Many variance-reduction methods require the definition of an importance map and exhibit optimal performance if the importance map is given by the adjoint flux, the solution of the adjoint Boltzmann equation. This paper presents the implementation of the Consistent Adjoint-Driven Importance Sampling (CADIS) methodology via a coupling between the TRIPOLI-4 R Monte-Carlo particle-transport code and the IDT deterministic solver, Additionally, we describe the implementation of a new TRIPOLI-4 R score that makes it possible to estimate the adjoint flux during a direct Monte-Carlo calculation. These new features are expected to simplify the solution of difficult shielding problems. (authors)

61 RADIATION PROTECTION AND DOSIMETRY↗

Data driven discovery and quantification of hyperspectral leaf reflectance phenotypes across a maize diversity panel

Abstract Estimates of plant traits derived from hyperspectral reflectance data have the potential to efficiently substitute for traits, which are time or labor intensive to manually score. Typical workflows for estimating plant traits from hyperspectral reflectance data employ supervised classification models that can require substantial ground truth datasets for training. We explore the potential of an unsupervised approach, autoencoders, to extract meaningful traits from plant hyperspectral reflectance data using measurements of the reflectance of 2151 individual wavelengths of light from the leaves of maize ( Zea mays ) plants harvested from 1658 field plots in a replicated field trial. A subset of autoencoder‐derived variables exhibited significant repeatability, indicating that a substantial proportion of the total variance in these variables was explained by difference between maize genotypes, while other autoencoder variables appear to capture variation resulting from changes in leaf reflectance between different batches of data collection. Several of the repeatable latent variables were significantly correlated with other traits scored from the same maize field experiment, including one autoencoder‐derived latent variable (LV8) that predicted plant chlorophyll content modestly better than a supervised model trained on the same data. In at least one case, genome‐wide association study hits for variation in autoencoder‐derived variables were proximal to genes with known or plausible links to leaf phenotypes expected to alter hyperspectral reflectance. In aggregate, these results suggest that an unsupervised, autoencoder‐based approach can identify meaningful and genetically controlled variation in high‐dimensional, high‐throughput phenotyping data and link identified variables back to known plant traits of interest.

Tross, Michael C.↗

Characterizing the structure of diffuse emission in Hi-GAL maps

We present a study of the structure of the Galactic interstellar medium (ISM) through the Δ-variance technique, related to the power spectrum and the fractal properties of infrared/submillimeter maps. Through this method, it is possible to provide quantitative parameters, which are useful for characterizing different morphological and physical conditions, and better constraining the theoretical models. In this respect, the Herschel Infrared Galactic Plane Survey, carried out at five photometric bands from 70 to 500 μm, constitutes a unique database for applying statistical tools to a variety of regions across the Milky Way. In this paper, we derive a robust estimate of the power-law portion of the power spectrum of four contiguous 2° × 2° Hi-GAL tiles located in the third Galactic quadrant (217° ≲ ℓ ≲ 225°, –2° ≲ b ≲ 0°). The low level of confusion along the line of sight, testified by CO observations, makes this region an ideal case. We find very different values for the power spectrum slope from tile to tile but also from wavelength to wavelength (2 ≲ β ≲ 3), with similarities between fields attributable to components located at the same distance. Thanks to comparisons with models of turbulence, an explanation of the determined slopes in terms of the fractal geometry is also provided, and possible relations with the underlying physics are investigated. In particular, an anti-correlation between ISM fractal dimension and star formation efficiency is found for the two main distance components observed in these fields. A possible link between the fractal properties of the diffuse emission and the resulting clump mass function is discussed.

79 ASTRONOMY AND ASTROPHYSICS↗

Aircraft NO/x/ emissions and stratospheric ozone reductions - Another look

New estimates for stratospheric ozone perturbations attributable to supersonic transport (SST) emissions are presented. First, a review is given of recent data pointing to lower OH concentrations below 30 km, as compared to the values predicted by photochemical models. The evidence for lower OH comes from a wide range of laboratory and atmospheric studies. The sensitivity of theoretical estimates of ozone change to OH abundances, and the coupling mechanisms between the O(x)-NO(x)-HO(x)-Cl(x) families which are responsible for the sensitivity, are discussed. Updated calculations for SST-induced ozone alterations are compared with older predictions. For example, assuming continuous aircraft injection of NO2 at 20 km at a rate of 1 x 10 to the 9th kg per year (globally), a 4% ozone decrease, is now calculated where earlier a 3% ozone increase was found. This large variance from previous forecasts suggests that new assessments of certain other polluting agents, particularly nitrogen fertilizers, are needed.

Turco, R. P.↗

Sensitivity Analysis for the Component Design App: Analysis of Success Assured Data

A new tool has been developed to perform variance-based global sensitivity analysis (VBGSA) on data from a set-based concurrent engineering software called Success Assured (SA). The tool is part of a digital component design app, which is currently in production as an Accelerated Digital Engineering Pathfinder at Sandia National Laboratories. When working with complex digital models, it is important to understand relationships between inputs and outputs, i.e., how “sensitive” model outputs are to changes in model inputs. After extensive research and trials of various sensitivity analysis methods, it was determined that estimation of Sobol’ indices for VBGSA with Monte Carlo simulation, paired with simple surrogate models, produces the best results for SA data. This tool increases understanding of SA models and streamlines the creation of SA datasets. This report details the methodology and implementation of this sensitivity analysis tool so others can understand it and implement it.

97 MATHEMATICS AND COMPUTING↗

Measurement of heat and moisture fluxes at the top of the rain forest during ABLE

Observations are presented of turbulent heat, moisture, and momentum transport made at two levels, approximately 5 and 10 m above the Amazon rain forest canopy. Data acquired at 10 Hz included variances and some mixed third moments of vertical velocity, temperature, and humidity. Two features of the data appear to question the displacement height hypothesis: (1) The characteristic dissipation length scale in the near-canopy layer varied between 20 m in stable conditions to approximately 150 m during afternoon convective conditions, generally larger scales than would be expected; and (2) No appreciable difference in dissipation scales was seen at the two observed levels. Observed peaks in vertical velocity-temperature cospectra lead to similar length scale estimates for dominant eddies. Heat budgets on selected days show that frequent periods with negative heat flux concurrent with continuing positive moisture flux occur in early afternoon, and this is believed to indicate the patchy nature of canopy-atmosphere coupling. Vertical velocity skewness was observed to be negative on three successive days and exhibited a sharp positive gradient.

Fitzjarrald, David R.↗

Soil moisture estimation using GOES-VISSR infrared data - A case study with a simple statistical method

Five days of clear sky observations of Kansas and Nebraska are used to examine the statistical relationship between soil moisture and infrared surface temperature observations taken from a geosynchronous satellite. Linear regression is used to relate soil moisture to surface temperature and other variables that represent wind speed, vegetation cover, and low-level temperature advection. Results show good agreement between estimated and observed soil moisture features on each of the 5 days. The average coefficient of determination for five pseudoindependent tests in which the test day is held out of the regression is 0.71. It is shown that a depletion coefficient of 0.92, when used to compute antecedent precipitation index (API), produces the best correlation between API and soil moisture as inferred from GOES thermal infrared data. By averaging daily predicted values over the 5-day rain-free case study period, 92 percent of the variance of the morning surface temperature change is explained by a simple multiple linear regression with all independent variables, or, alternatively, 85 percent of the observed variance in API is explained. It is concluded that this approach can distinguish at least four classes of soil wetness, but the necessity for measurement of surface advection may limit its usefulness in remote areas.

Wetzel, Peter J.↗

Systematic methods for knowledge acquisition and expert system development

Nine cooperating rule-based systems, collectively called AUTOCREW which were designed to automate functions and decisions associated with a combat aircraft's subsystems, are discussed. The organization of tasks within each system is described; performance metrics were developed to evaluate the workload of each rule base and to assess the cooperation between the rule bases. Simulation and comparative workload results for two mission scenarios are given. The scenarios are inbound surface-to-air-missile attack on the aircraft and pilot incapacitation. The methodology used to develop the AUTOCREW knowledge bases is summarized. Issues involved in designing the navigation sensor selection expert in AUTOCREW's NAVIGATOR knowledge base are discussed in detail. The performance of seven navigation systems aiding a medium-accuracy INS was investigated using Kalman filter covariance analyses. A navigation sensor management (NSM) expert system was formulated from covariance simulation data using the analysis of variance (ANOVA) method and the ID3 algorithm. ANOVA results show that statistically different position accuracies are obtained when different navaids are used, the number of navaids aiding the INS is varied, the aircraft's trajectory is varied, and the performance history is varied. The ID3 algorithm determines the NSM expert's classification rules in the form of decision trees. The performance of these decision trees was assessed on two arbitrary trajectories, and the results demonstrate that the NSM expert adapts to new situations and provides reasonable estimates of the expected hybrid performance.

Belkin, Brenda L.↗

Early prediction of antigenic transitions for influenza A/H3N2

Influenza A/H3N2 is a rapidly evolving virus which experiences major antigenic transitions every two to eight years. Anticipating the timing and outcome of transitions is critical to developing effective seasonal influenza vaccines. Using a published phylodynamic model of influenza transmission, we identified indicators of future evolutionary success for an emerging antigenic cluster and quantified fundamental trade-offs in our ability to make such predictions. The eventual fate of a new cluster depends on its initial epidemiological growth rate––which is a function of mutational load and population susceptibility to the cluster––along with the variance in growth rate across co-circulating viruses. Logistic regression can predict whether a cluster at 5% relative frequency will eventually succeed with ~80% sensitivity, providing up to eight months advance warning. As a cluster expands, the predictions improve while the lead-time for vaccine development and other interventions decreases. However, attempts to make comparable predictions from 12 years of empirical influenza surveillance data, which are far sparser and more coarse-grained, achieve only 56% sensitivity. By expanding influenza surveillance to obtain more granular estimates of the frequencies of and population-wide susceptibility to emerging viruses, we can better anticipate major antigenic transitions. This provides added incentives for accelerating the vaccine production cycle to reduce the lead time required for strain selection.

59 BASIC BIOLOGICAL SCIENCES↗

Behavior and mechanisms of Doppler wind lidar error in varying stability regimes

Abstract. Wind lidars are widespread and important tools in atmospheric observations. An intrinsic part of lidar measurement error is due to atmospheric variability in the remote-sensing scan volume. This study describes and quantifies the distribution of measurement error due to turbulence in varying atmospheric stability. While the lidar error model is general, we demonstrate the approach using large ensembles of virtual WindCube V2 lidar performing a profiling Doppler-beam-swinging scan in quasi-stationary large-eddy simulations (LESs) of convective and stable boundary layers. Error trends vary with the stability regime, time averaging of results, and observation height. A systematic analysis of the observation error explains dominant mechanisms and supports the findings of the empirical results. Treating the error under a random variable framework allows for informed predictions about the effect of different configurations or conditions on lidar performance. Convective conditions are most prone to large errors (up to 1.5 m s−1 in 1 Hz wind speed in strong convection), driven by the large vertical velocity variances in convective conditions and the high elevation angle of the scanning beams (62∘). Range-gate weighting induces a negative bias into the horizontal wind speeds near the surface shear layer (−0.2 m s−1 in the stable test case). Errors in the horizontal wind speed and direction computed from the wind components are sensitive to the background wind speed but have negligible dependence on the relative orientation of the instrument. Especially during low winds and in the presence of large errors in the horizontal velocity estimates, the reported wind speed is subject to a systematic positive bias (up to 0.4 m s−1 in 1 Hz measurements in strong convection). Vector time-averaged measurements can improve the behavior of the error distributions (reducing the 10 min wind speed error standard deviation to <0.3 m s−1 and the bias to <0.1 m s−1 in strong convection) with a predictable effectiveness related to the number of decorrelated samples in the time window. Hybrid schemes weighting the 10 min scalar- and vector-averaged lidar measurements are shown to be effective at reducing the wind speed biases compared to cup measurements in most of the simulated conditions, with time averages longer than 10 min recommended for best use in some unstable conditions. The approach in decomposing the error mechanisms with the help of the LES flow field could be extended to more complex measurement scenarios and scans.

17 WIND ENERGY↗

Uncertainty Quantification and Global Sensitivity Analysis of Subsurface Flow Parameters to Gravimetric Variations During Pumping Tests in Unconfined Aquifers

We study the contribution of typically uncertain subsurface flow parameters to gravity changes that can be recorded during pumping tests in unconfined aquifers. We do so in the framework of a Global Sensitivity Analysis and quantify the effects of uncertainty of such parameters on the first four statistical moments of the probability distribution of gravimetric variations induced by the operation of the well. System parameters are grouped into two main categories, respectively, governing groundwater flow in the unsaturated and saturated portions of the domain. We ground our work on the three-dimensional analytical model proposed by Mishra and Neuman (2011), which fully takes into account the richness of the physical process taking place across the unsaturated and saturated zones and storage effects in a finite radius pumping well. The relative influence of model parameter uncertainties on drawdown, moisture content, and gravity changes are quantified through (a) the Sobol' indices, derived from a classical decomposition of variance and (b) recently developed indices quantifying the relative contribution of each uncertain model parameter to the (ensemble) mean, skewness, and kurtosis of the model output. Our results document (i) the importance of the effects of the parameters governing the unsaturated flow dynamics on the mean and variance of local drawdown and gravity changes; (ii) the marked sensitivity (as expressed in terms of the statistical moments analyzed) of gravity changes to the employed water retention curve model parameter, specific yield, and storage, and (iii) the influential role of hydraulic conductivity of the unsaturated and saturated zones to the skewness and kurtosis of gravimetric variation distributions. The observed temporal dynamics of the strength of the relative contribution of system parameters to gravimetric variations suggest that gravity data have a clear potential to provide useful information for estimating the key hydraulic parameters of the system.

54 ENVIRONMENTAL SCIENCES↗

Meta-analysis of the Alpha/Beta Ratio for Prostate Cancer in the Presence of an Overall Time Factor: Bad News, Good News, or No News?

Purpose: To present a novel method for meta-analysis of the fractionation sensitivity of tumors as applied to prostate cancer in the presence of an overall time factor. Methods and Materials: A systematic search for radiation dose-fractionation trials in prostate cancer was performed using PubMed and by manual search. Published trials comparing standard fractionated external beam radiation therapy with alternative fractionation were eligible. For each trial the {alpha}/{beta} ratio and its 95% confidence interval (CI) were extracted, and the data were synthesized with each study weighted by the inverse variance. An overall time factor was included in the analysis, and its influence on {alpha}/{beta} was investigated. Results: Five studies involving 1965 patients were included in the meta-analysis of {alpha}/{beta}. The synthesized {alpha}/{beta} assuming no effect of overall treatment time was -0.07 Gy (95% CI -0.73-0.59), which was increased to 0.47 Gy (95% CI -0.55-1.50) if a single highly weighted study was excluded. In a separate analysis, 2 studies based on 10,808 patients in total allowed extraction of a synthesized estimate of a time factor of 0.31 Gy/d (95% CI 0.20-0.42). The time factor increased the {alpha}/{beta} estimate to 0.58 Gy (95% CI -0.53-1.69)/1.93 Gy (95% CI -0.27-4.14) with/without the heavily weighted study. An analysis of the uncertainty of the {alpha}/{beta} estimate showed a loss of information when the hypofractionated arm was underdosed compared with the normo-fractionated arm. Conclusions: The current external beam fractionation studies are consistent with a very low {alpha}/{beta} ratio for prostate cancer, although the CIs include {alpha}/{beta} ratios up to 4.14 Gy in the presence of a time factor. Details of the dose fractionation in the 2 trial arms have critical influence on the information that can be extracted from a study. Studies with unfortunate designs will supply little or no information about {alpha}/{beta} regardless of the number of subjects enrolled.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Smoothing Lexis diagrams using kernel functions: A contemporary approach

Lexis diagrams are rectangular arrays of event rates indexed by age and period. Analysis of Lexis diagrams is a cornerstone of cancer surveillance research. Typically, population-based descriptive studies analyze multiple Lexis diagrams defined by sex, tumor characteristics, race/ethnicity, geographic region, etc. Inevitably the amount of information per Lexis diminishes with increasing stratification. Several methods have been proposed to smooth observed Lexis diagrams up front to clarify salient patterns and improve summary estimates of averages, gradients, and trends. In this article, we develop a novel bivariate kernel-based smoother that incorporates two key innovations. First, for any given kernel, we calculate its singular values decomposition, and select an optimal truncation point—the number of leading singular vectors to retain—based on the bias-corrected Akaike information criterion. Second, we model-average over a panel of candidate kernels with diverse shapes and bandwidths. The truncated model averaging approach is fast, automatic, has excellent performance, and provides a variance-covariance matrix that takes model selection into account. We present an in-depth case study (invasive estrogen receptor-negative breast cancer incidence among non-Hispanic white women in the United States) and simulate operating characteristics for 20 representative cancers. The truncated model averaging approach consistently outperforms any fixed kernel. Our results support the routine use of the truncated model averaging approach in descriptive studies of cancer.

60 APPLIED LIFE SCIENCES↗

Parameter estimation applied to Nimbus 6 wide-angle longwave radiation measurements

A parameter estimation technique was used to analyze the August 1975 Nimbus 6 Earth radiation budget data to demonstrate the concept of deconvolution. The longwave radiation field at the top of the atmosphere is defined from satellite data by a fifth degree and fifth order spherical harmonic representation. The variations of the major features of the radiation field are defined by analyzing the data separately for each two-day duty cycle. A table of coefficient values for each spherical harmonic representation is given along with global mean, gradients, degree variances, and contour plots. In addition, the entire data set is analyzed to define the monthly average radiation field.

Green, R. N.↗

Software Estimates Costs of Testing Rocket Engines

Simulation-Based Cost Model (SiCM) is a computer program that simulates pertinent aspects of the testing of rocket engines for the purpose of estimating the costs of such testing during time intervals specified by its users. A user enters input data for control of simulations; information on the nature of, and activity in, a given testing project; and information on resources. Simulation objects are created on the basis of this input. Costs of the engineering-design, construction, and testing phases of a given project are estimated from numbers and labor rates of engineers and technicians employed in each phase, the duration of each phase; costs of materials used in each phase; and, for the testing phase, the rate of maintenance of the testing facility. The three main outputs of SiCM are (1) a curve, updated at each iteration of the simulation, that shows overall expenditures vs. time during the interval specified by the user; (2) a histogram of the total costs from all iterations of the simulation; and (3) table displaying means and variances of cumulative costs for each phase from all iterations. Other outputs include spending curves for each phase.

Source record↗

Software Estimates Costs of Testing Rocket Engines

Simulation-Based Cost Model (SiCM), a discrete event simulation developed in Extend , simulates pertinent aspects of the testing of rocket propulsion test articles for the purpose of estimating the costs of such testing during time intervals specified by its users. A user enters input data for control of simulations; information on the nature of, and activity in, a given testing project; and information on resources. Simulation objects are created on the basis of this input. Costs of the engineering-design, construction, and testing phases of a given project are estimated from numbers and labor rates of engineers and technicians employed in each phase, the duration of each phase; costs of materials used in each phase; and, for the testing phase, the rate of maintenance of the testing facility. The three main outputs of SiCM are (1) a curve, updated at each iteration of the simulation, that shows overall expenditures vs. time during the interval specified by the user; (2) a histogram of the total costs from all iterations of the simulation; and (3) table displaying means and variances of cumulative costs for each phase from all iterations. Other outputs include spending curves for each phase.

Smith, C. L.↗

Dark Energy Survey Year 3 results: cosmology with moments of weak lensing mass maps – validation on simulations

ABSTRACT We present a simulated cosmology analysis using the second and third moments of the weak lensing mass (convergence) maps. The second moment, or variances, of the convergence as a function of smoothing scale contains information similar to standard shear two-point statistics. The third moment, or the skewness, contains additional non-Gaussian information. The analysis is geared towards the third year (Y3) data from the Dark Energy Survey (DES), but the methodology can be applied to other weak lensing data sets. We present the formalism for obtaining the convergence maps from the measured shear and for obtaining the second and third moments of these maps given partial sky coverage. We estimate the covariance matrix from a large suite of numerical simulations. We test our pipeline through a simulated likelihood analyses varying 5 cosmological parameters and 10 nuisance parameters and identify the scales where systematic or modelling uncertainties are not expected to affect the cosmological analysis. Our simulated likelihood analysis shows that the combination of second and third moments provides a 1.5 per cent constraint on S8 ≡ σ8(Ωm/0.3)0.5 for DES Year 3 data. This is 20 per cent better than an analysis using a simulated DES Y3 shear two-point statistics, owing to the non-Gaussian information captured by the inclusion of higher order statistics. This paper validates our methodology for constraining cosmology with DES Year 3 data, which will be presented in a subsequent paper.

79 ASTRONOMY AND ASTROPHYSICS↗