Search NASA⌕ Search

SEARCH · Search NASA

Results for “Variance to Mean”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 91 records · Page 5

The neutron number probability distribution in coupled lumped assemblies

Here, the validity of the gamma distribution in describing the neutron number probability distribution function for both isolated and coupled multiplying assemblies when constrained to reproduce the true mean and variance is investigated in lumped geometry by numerical comparison with kinetic Monte Carlo simulations. The mean and variance are obtained from numerical solution of moment equations constructed from the relevant forward Master equation with assembly coupling coefficients obtained from a view factor method. Numerical results for a two-group, two coupled assemblies model, with static and dynamic reactivity insertion, show that except for subcritical assemblies, the gamma distribution well-approximates the number distribution. Differences in the fast and thermal neutron population shapes are explained in terms of effective source strengths due to downscatter and coupling.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

A new multi-model absolute difference-based sensitivity (MMADS) analysis method to screen non-influential processes under process model and parametric uncertainty

Process-based models have been widely used for hydrologic modeling, and it is a common practice to use sensitivity analysis methods for excluding non-influential hydrologic processes from further investigation and/or model improvement. This study develops a new method called multi-model absolute difference-based sensitivity (MMADS) analysis method to screen non-influential system processes and parameters. MMADS is conceptually similar to the Morris method for addressing parametric uncertainty, but has a unique feature to address both process model uncertainty (i.e., a process may be represented by multiple process models) and process model parameter uncertainty (i.e., parameters associated with a process model are random). MMADS first evaluates absolute differences of a quantity of interest (i.e., a system model output) by varying process models and/or process model parameter values, and then calculates the mean and variance of the differences for investigating process influence. The mean measures overall influence of the process on the quantity of interest, and the variance estimates influence of nonlinear effects of the process and/or its interactions with other processes. MMADS is an extension of the Morris method from a parameter space to a joint parameter-model space for explicitly addressing both process model uncertainty and model parameter uncertainty. The performance of MMADS is evaluated by using two numerical experiments. One experiment is based on Sobol’s G*-function with ten product elements, and has analytical solutions of the MMADS mean and variance of absolute differences. The other experiment is for groundwater flow modeling which considers three processes (i.e., recharge, geology, and snowmelt) that interact with each other. Finally, results indicate that MMADS is computationally efficient and can identify non-influential processes of complex hydrological systems.

54 ENVIRONMENTAL SCIENCES↗

Stochastic Kinetics of Nanocatalytic Systems

Catalytic reaction events occurring on the surface of a nanoparticle constitute a complex stochastic process. Although advances in modern single-molecule experiments enable direct measurements of individual catalytic turnover events occurring on a segment of a single nanoparticle, we do not yet know how to measure the number of catalytic sites in each segment or how the catalytic turnover counting statistics and the catalytic turnover time distribution are related to the microscopic dynamics of catalytic reactions. Here, we address these issues by presenting a stochastic kinetics for nanoparticle catalytic systems. We propose a new experimental measure of the number of catalytic sites in terms of the mean and variance of the catalytic event count. By considering three types of nanocatalytic systems, we investigate how the mean, the variance, and the distribution of the catalytic turnover time depend on the catalytic reaction dynamics, the heterogeneity of catalytic activity, and communication among catalytic sites. This work enables accurate quantitative analyses of single-molecule experiments for nanocatalytic systems and enzymes with multiple catalytic sites.

Statistical Physics↗

Data-driven predictive probability density function control of fiber length stochastic distribution shaping in refining process

Pulp is the most important raw material in paper in paper industries, whose Fiber Length Stochastic Distribution (FLSD) shaping directly determines the energy consumption and paper quality of the subsequent papermaking processes. However, the mean and variance are insufficient to describe the output FLSD shaping, which display non-Gaussian distribution properties. Therefore, the traditional control, method based on the mean and variance of fiber length is difficult to control the FLSD shaping effectively. In this paper, a novel data-driven predictive probability density function (PDF) control method is proposed for output FLSD shaping in refining process. Primarily, in order to improve the approximation accuracy of the PDF employing RBF neural network, the parameters (center value and width) of each RBF basis functions are tuned via utilizing iterative learning control (ILC) law, and the corresponding estimations of weights law can be obtained. Secondly, considering that the conventional linear model of weights vector has drawbacks of low accuracy and weak generalization ability, random vector functional link networks (RVFLNs) based data driven nonlinear modelling method is employed to characterize the prediction model between the input variables and weights vector. Finally, in order to reduce the randomness of the output FLSD, a minimum entropy control method under mean constraint is employed to design the predictive PDF controller for the output FLSD shaping with the help of the established nonlinear stochastic distribution model. Both simulations and experiments are given to illustrate the effectiveness and practicability of the proposed method.

Probability debsity function, Fiber length distrib↗

Inverse sequential procedures for the monitoring of time series

Climate changes traditionally have been detected from long series of observations and long after they happened. The 'inverse sequential' monitoring procedure is designed to detect changes as soon as they occur. Frequency distribution parameters are estimated both from the most recent existing set of observations and from the same set augmented by 1,2,...j new observations. Individual-value probability products ('likelihoods') are then calculated which yield probabilities for erroneously accepting the existing parameter(s) as valid for the augmented data set and vice versa. A parameter change is signaled when these probabilities (or a more convenient and robust compound 'no change' probability) show a progressive decrease. New parameters are then estimated from the new observations alone to restart the procedure. The detailed algebra is developed and tested for Gaussian means and variances, Poisson and chi-square means, and linear or exponential trends; a comprehensive and interactive Fortran program is provided in the appendix.

Radok, Uwe↗

Global Gravity Wave Variances from Aura MLS: Characteristics and Interpretation

The gravity wave (GW)-resolving capabilities of 118-GHz saturated thermal radiances acquired throughout the stratosphere by the Microwave Limb Sounder (MLS) on the Aura satellite are investigated and initial results presented. Because the saturated (optically thick) radiances resolve GW perturbations from a given altitude at different horizontal locations, variances are evaluated at 12 pressure altitudes between 21 and 51 km using the 40 saturated radiances found at the bottom of each limb scan. Forward modeling simulations show that these variances are controlled mostly by GWs with vertical wavelengths z 5 km and horizontal along-track wavelengths of y 100-200 km. The tilted cigar-shaped three-dimensional weighting functions yield highly selective responses to GWs of high intrinsic frequency that propagate toward the instrument. The latter property is used to infer the net meridional component of GW propagation by differencing the variances acquired from ascending (A) and descending (D) orbits. Because of improved vertical resolution and sensitivity, Aura MLS GW variances are 5?8 times larger than those from the Upper Atmosphere Research Satellite (UARS) MLS. Like UARS MLS variances, monthly-mean Aura MLS variances in January and July 2005 are enhanced when local background wind speeds are large, due largely to GW visibility effects. Zonal asymmetries in variance maps reveal enhanced GW activity at high latitudes due to forcing by flow over major mountain ranges and at tropical and subtropical latitudes due to enhanced deep convective generation as inferred from contemporaneous MLS cloud-ice data. At 21-28-km altitude (heights not measured by the UARS MLS), GW variance in the tropics is systematically enhanced and shows clear variations with the phase of the quasi-biennial oscillation, in general agreement with GW temperature variances derived from radiosonde, rocketsonde, and limb-scan vertical profiles.

MLS (MICROWAVE LIMB SOUNDER)↗

Convergence of oscillator spectral estimators for counted-frequency measurements.

A common intermediary connecting frequency-noise calibration or testing of an oscillator to useful applications is the spectral density of the frequency-deviating process. In attempting to turn test data into predicts of performance characteristics, one is naturally led to estimation of statistical values by sample-mean and sample-variance techniques. However, sample means and sample variances themselves are statistical quantities that do not necessarily converge (in the mean-square sense) to actual ensemble-average means and variances, except perhaps for excessively large sample sizes. This is especially true for the flicker noise component of oscillators. This article shows, for the various types of noises found in oscillators, how sample averages converge (or do not converge) to their statistical counterparts. The convergence rate is shown to be the same for all oscillators of a given spectral type.

Tausworthe, R. C.↗

Redshift data and statistical inference

Frequency histograms and the 'power spectrum analysis' (PSA) method, the latter developed by Yu & Peebles (1969), have been widely employed as techniques for establishing the existence of periodicities. We provide a formal analysis of these two classes of methods, including controlled numerical experiments, to better understand their proper use and application. In particular, we note that typical published applications of frequency histograms commonly employ far greater numbers of class intervals or bins than is advisable by statistical theory sometimes giving rise to the appearance of spurious patterns. The PSA method generates a sequence of random numbers from observational data which, it is claimed, is exponentially distributed with unit mean and variance, essentially independent of the distribution of the original data. We show that the derived random processes is nonstationary and produces a small but systematic bias in the usual estimate of the mean and variance. Although the derived variable may be reasonably described by an exponential distribution, the tail of the distribution is far removed from that of an exponential, thereby rendering statistical inference and confidence testing based on the tail of the distribution completely unreliable. Finally, we examine a number of astronomical examples wherein these methods have been used giving rise to widespread acceptance of statistically unconfirmed conclusions.

Newman, William I.↗

Nonlinear probabilistic finite element models of laminated composite shells

A probabilistic finite element analysis procedure for laminated composite shells has been developed. A total Lagrangian finite element formulation, employing a degenerated 3-D laminated composite shell with the full Green-Lagrange strains and first-order shear deformable kinematics, forms the modeling foundation. The first-order second-moment technique for probabilistic finite element analysis of random fields is employed and results are presented in the form of mean and variance of the structural response. The effects of material nonlinearity are included through the use of a rate-independent anisotropic plasticity formulation with the macroscopic point of view. Both ply-level and micromechanics-level random variables can be selected, the latter by means of the Aboudi micromechanics model. A number of sample problems are solved to verify the accuracy of the procedures developed and to quantify the variability of certain material type/structure combinations. Experimental data is compared in many cases, and the Monte Carlo simulation method is used to check the probabilistic results. In general, the procedure is quite effective in modeling the mean and variance response of the linear and nonlinear behavior of laminated composite shells.

Engelstad, S. P.↗

Dynamic and Regression Modeling of Ocean Variability in the Tide-Gauge Record at Seasonal and Longer Periods

Comparison of monthly mean tide-gauge time series to corresponding model time series based on a static inverted barometer (IB) for pressure-driven fluctuations and a ocean general circulation model (OM) reveals that the combined model successfully reproduces seasonal and interannual changes in relative sea level at many stations. Removal of the OM and IB from the tide-gauge record produces residual time series with a mean global variance reduction of 53%. The OM is mis-scaled for certain regions, and 68% of the residual time series contain a significant seasonal variability after removal of the OM and IB from the tide-gauge data. Including OM admittance parameters and seasonal coefficients in a regression model for each station, with IB also removed, produces residual time series with mean global variance reduction of 71%. Examination of the regional improvement in variance caused by scaling the OM, including seasonal terms, or both, indicates weakness in the model at predicting sea-level variation for constricted ocean regions. The model is particularly effective at reproducing sea-level variation for stations in North America, Europe, and Japan. The RMS residual for many stations in these areas is 25-35 mm. The production of "cleaner" tide-gauge time series, with oceanographic variability removed, is important for future analysis of nonsecular and regionally differing sea-level variations. Understanding the ocean model's strengths and weaknesses will allow for future improvements of the model.

Hill, Emma M.↗

A Lumped Stochastic Model of Coupled Neutronic Assemblies

The objective of the present work is to apply the Harris ansatz to coupled assemblies in a lumped setting, i.e., ignoring the neutron phase dependence, which hitherto has not been considered. A Master equation for the single-time joint distribution of neutron number in each assembly, accounting for transfer of neutrons from one assembly to another via leakage transition probabilities and intercept or view factors is derived from which equations for the mean and variance in each assembly are developed. Because of the neutronic coupling, it is necessary to additionally compute the correlation function between pairs of assemblies which gives an expanded set of moment equations to solve. The reactivity of each assembly is allowed to vary with time, although a formal mechanism that is responsible for reactivity coupling between assemblies is not considered or proposed. Following Harris, it is then hypothesized that the number distribution in each assembly is a gamma distribution with assembly-specific mean and variance computed exactly from the moment equations. Illustrative numerical results are presented for the moments of a system of four coupled assemblies under various reactivity and neutronic coupling scenarios. The gamma distribution is then used to compute the probability that the neutron population exceeds a threshold value in any assembly. The report concludes with a discussion of recommended further work on this problem.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Multilevel Monte Carlo Predictions of First Passage Times in Three-Dimensional Discrete Fracture Networks: A Graph-Based Approach

In this paper, we present a method combining multilevel Monte Carlo (MLMC) and a graph-based primary subnetwork identification algorithm to provide estimates of the mean and variance of the distribution of first passage times in fracture media at significantly lower computational cost than standard Monte Carlo (MC) methods. Simulations of solute transport are performed using a discrete fracture network (DFN), and instead of using various grid resolutions for levels in the MLMC, which is standard practice in MLMC, we identify a hierarchy of subnetworks in the DFN based on the shortest topological paths through the network using a graph-based method. While the mean of these ensembles is of critical importance, the variance is also essential in fractured media where uncertainty is an overarching theme, and understanding variability across an ensemble is a requirement for safety assessments. The method provides good estimates of the mean and variance at two orders of magnitude lower computational cost than MC.

54 ENVIRONMENTAL SCIENCES↗

Statistical moments of VIIRS night-time lights

Annual VIIRS global night-time lights traditionally report mean radiance levels after filtering to remove moonlit and cloudy data followed by outlier removal. The mean is the first statistical moment, but there are three other statistical moments (variance, skew and kurtosis) that are currently not produced globally. In this study, we generated VIIRS day/night band multiyear and annual statistical moments for a widely dispersed set of test areas. The moments were calculated from 15 arc second nightly temporal profiles spanning 2012-2020, filtered to exclude cloudy and sunlit data, with radiance adjustments to reduce view angle and lunar illuminance effects. The moment data were examined in two ways: (1) Geospatial grids-which reveal zonation and temporal changes present in urban areas and (2) Scattergrams of moment pairs. The variance versus mean scattergram exhibits several distinct data clusters used to define five zones: core lighting, dark-erratics, mid-erratics, bright-erratics and bright and seteady. The results indicate that the information content of VIIRS night-time lights can be augmented if all four statistical moments are calculated.

47 OTHER INSTRUMENTATION↗

Future of the MUSiC Experiment Data

The Measurements of Uranium Subcritical and Critical (MUSiC) experiment was a highly enriched uranium (HEU) experiment performed at the National Criticality Experiments Research Center (NCERC) executed between December 2020 and April 2021. The experiment intended to measure criticality and reactor kinetics parameters in a bare HEU system. The experiment concurrently measured radiation signatures from the system while utilizing different neutron source types. This was an attempt to benchmark both detectors and analysis techniques against one another for identical measurements. The experiment consisted of the Rocky Flats HEU hemi-shells constructed into ten configurations spanning between deeply subcritical (about 14 kgs) to supercritical (about 60 kgs). Two of the ten configurations were supercritical and the other eight were subcritical. The two supercritical configurations, often referred to as the critical configurations, were documented into an International Criticality Safety Benchmark Evaluation Project (ICSBEP) benchmark evaluation and submitted to the technical review group (TRG). The subcritical configurations of the MUSiC experiment were examined using three different detection systems. The systems include: the NoMAD He-3 neutron detector which consists of 15 He-3 tubes surrounded by a polyethylene matrix, a liquid scintillator system called the Rossi-α Measurement Rapid Organic Discriminating Detector (RAM-RODD), and a trans-stilbene organic scintillator array (OSCAR) provided by the University of Michigan. The NoMAD and RAM-RODD data are planned to be evaluated in two separate ICSBEP evaluations utilizing neutron noise methods such as Feynman Variance-to-Mean, Rossi-α and the pulsed neutron source method. Each of the subcritical configurations were examined using three different source types including: a Cf-252 source in the center, with only the intrinsic neutron source in the HEU, and with an external D-T neutron generator.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Accurate Simulation of Both Sensitivity and Variability for Amazonian Photosynthesis: Is It Too Much to Ask?

Abstract Estimates of Amazon rainforest gross primary productivity (GPP) differ by a factor of 2 across a suite of three statistical and 18 process models. This wide spread contributes uncertainty to predictions of future climate. We compare the mean and variance of GPP from these models to that of GPP at six eddy covariance (EC) towers. Only one model's mean GPP across all sites falls within a 99% confidence interval for EC GPP, and only one model matches EC variance. The strength of model response to climate drivers is related to model ability to match the seasonal pattern of the EC GPP. Models with stronger seasonal swings in GPP have stronger responses to rain, light, and temperature than does EC GPP. The model to data comparison illustrates a trade‐off inherent to deterministic models between accurate simulation of a mean (average) and accurate responsiveness to drivers. The trade‐off exists because all deterministic models simplify processes and lack at least some consequential driver or interaction. If a model's sensitivities to included drivers and their interactions are accurate, then deterministically predicted outcomes have less variability than is realistic. If a GPP model has stronger responses to climate drivers than found in data, model predictions may match the observed variance and seasonal pattern but are likely to overpredict GPP response to climate change. High or realistic variability of model estimates relative to reference data indicate that the model is hypersensitive to one or more drivers.

54 ENVIRONMENTAL SCIENCES↗

Analysis of the MUSIC 3 He Multiplicity Data

A measurement campaign called the Measurement of Uranium Subcritical and Critical (MUSiC) was performed on a range of configurations of highly-enriched uranium (HEU) from December 2020 through April of 2021. While part of the focus was to measure reactor kinetics parameters on delayed supercritical systems, an additional focus was performing neutron noise measurements on subcritical configurations from deeply subcritical to nearly delayed critical. Multiple detector systems were used to perform these measurements, such as a 3 He multiplicity detector called the Neutron Multiplicity Array Detector (NoMAD) and a liquid scintillator system called the Rossi-α Measurement Rapid Organic Discriminating Detector (RAM-RODD). Also included were a scintillator system from the University of Michigan and a set of small 3He tubes that have previously been used to measure Rossi-α values on near-critical systems. The focus of this paper will be a comparison of prospective analysis methods for the NoMAD measurements. Previous subcritical measurements at the National Criticality Experiments Research Center (NCERC) submitted to the International Criticality Safety Benchmark Evaluation Project (ICSBEP) used the Hage-Cifarelli formalism of the Feynman Variance-to-Mean method. This relies on the time correlations of neutron detections to infer the spontaneous fission rate and neutron multiplication of a system through binning the time tagged detections and analyzing resulting histograms of the numbers of counts. However, there are other neutron noise methods that rely on similar processes, such as the Hansen-Dowdy formalism which uses a slightly different methodology to extract multiplication from the neutron multiplicity counting moments. Comparisons can be made between these experimental results and those obtained through simulations to validate or identify deficiencies in analysis, detection methods, or the underlying nuclear data. Different time gating strategies and their effects on count rate uncertainties are also investigated.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Estimation of parameters in linear structural relationships: Sensitivity to the choice of the ratio of error variances

Maximum likelihood estimation of parameters in linear structural relationships under normality assumptions requires knowledge of one or more of the model parameters if no replication is available. The most common assumption added to the model definition is that the ratio of the error variances of the response and predictor variates is known. The use of asymptotic formulae for variances and mean squared errors as a function of sample size and the assumed value for the error variance ratio is investigated.

Lakshminarayanan, M. Y.↗

Estimation of parameters in linear structural relationships Sensitivity to the choice of the ratio of error variances

Maximum likelihood estimation of parameters in linear structural relationships under normality assumptions requires knowledge of one or more of the model parameters if no replication is available. The most common assumption added to the model definition is that the ratio of the error variances of the response and predictor variates is known. This paper investigates the use of asymptotic formulae for variances and mean squared errors as a function of sample size and the assumed value for the error variance ratio.

Lakshminarayanan, M. Y.↗