Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “sampling statistics”

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

Analysis of Harrell Monosodium Titanate Lot #46000824120

Monosodium titanate (MST) for use in the Actinide Removal Process (ARP) must be qualified and verified in advance. A single qualification sample for each batch of material is sent to SRNL for analysis, as well as a statistical sampling of verification samples. The original Harrell Industries Lot #46000824120 qualification and 16 verification samples received in September 2012 failed to meet the specification for weight percent solids. All of the pails sampled and tested contained less than 15 wt % MST solids. The lot was returned to the vendor, and in February 2014 a new qualification sample and set of 14 verification samples were received from this lot. The new lot met each of the selected specification requirements that were tested and, consequently, the material is acceptable for use in the ARP process.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Analysis of Harrell Monosodium Titanate Lot #46000908120

Monosodium titanate (MST) for use in the Actinide Removal Process (ARP) must be qualified and verified in advance. A single qualification sample for each batch of material is sent to SRNL for analysis, as well as a statistical sampling of verification samples. The original Harrell Industries Lot #46000908120 qualification and 16 verification samples received in October 2012 failed to meet the specification for weight percent solids. All of the pails sampled and tested contained less than 15 wt % MST solids. The lot was returned to the vendor, and in February 2014 a new qualification sample and set of 16 verification samples were received from this lot. The new lot met each of the selected specification requirements that were tested and, consequently, the material is acceptable for use in the ARP process.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Parametric description of intermittent probability distribution functions in solar wind and magnetohydrodynamic turbulence

ABSTRACT In this work, we find empirical evidence that the scale-dependent statistical properties of solar wind and magnetohydrodynamic (MHD) turbulence can be described in terms of a family of parametric probability distribution functions (PDFs) known as Normal Inverse Gaussian (NIG). Understanding these PDFs is one of the most important goals in turbulence theory, as they are inherently connected to the intermittent properties of solar wind turbulence. We investigate the properties of PDFs of Elsasser increments based on a large statistical sample from solar wind observations and high-resolution numerical simulations of MHD turbulence. In order to measure the PDFs and their corresponding properties, three experiments are presented: fast and slow solar wind for experimental data and a simulation of reduced MHD (RMHD) turbulence. Conditional statistics on a 23-yr-long sample of WIND data near 1 au and high-resolution pseudo-spectral simulation of steadily driven RMHD turbulence on a $2048^3$ mesh are used to construct scale-dependent PDFs. The empirical PDFs are fitted to NIG distributions, which depend on four free parameters. Our analysis shows that NIG distributions accurately capture the evolution of the PDFs, with scale-dependent parameters, from large scales characterized by a Gaussian distribution, turning to exponential tails within the inertial range and stretched exponentials at dissipative scales. We also show that empirically-measured NIG parameters exhibit well-defined scaling properties that are similar across the three empirical data sets, which may be indicative of universal behaviour.

Astronomy & Astrophysics↗

Evolution of Galaxies and the Star Formation Rate in the Infrared

A central goal of extragalactic observational astronomy is to understand how normal galaxies evolve with redshift, and particularly when galaxies formed their stars. While optical and rest-frame UV (ultraviolet) observations have begun to address these issues, the interpretation of such data is particularly challenging because of the sensitivity to dust obscureness (at optical and UV wavelengths). The absorbed light is re-radiated at IR (infrared) wavelengths, hence the optimal indicators of the star formation rate (SFR) is at a rest-frame of (lambda) (is approximately equal to) 60 microns. The SIRTF (space infrared telescope facility) mission will revolutionize the study of the global evolution of the SFR by providing mass-selected, complete samples of galaxies and fares estimators of the SFR. This research program is to study the SFR using statistical samples of galaxies in the local universe, at intermediate redshifts, and set the stage for continuing studies up to z = 5.

Pahre, Michael, A.↗

New giant planet beyond the snow line for an extended MOA exoplanet microlens sample

Characterizing a planet detected by microlensing is hard if the planetary signal is weak or the lens-source relative trajectory is far from caustics. However, statistical analyses of planet demography must include those planets to accurately determine occurrence rates. As part of a systematic modelling effort in the context of a >10-yr retrospective analysis of MOA’s survey observations to build an extended MOA statistical sample, we analyse the light curve of the planetary microlensing event MOA-2014-BLG-472. This event provides weak constraints on the physical parameters of the lens, as a result of a planetary anomaly occurring at low magnification in the light curve. We use a Bayesian analysis to estimate the properties of the planet, based on a refined Galactic model and the assumption that all Milky Way’s stars have an equal planet-hosting probability. We find that a lens consisting of a 1.9(+2.2,−1.2)M(J) giant planet orbiting a 0.31(+0.36,−0.19)Mꙩ host at a projected separation of 0.75±0.24au is consistent with the observations and is most likely, based on the Galactic priors. The lens most probably lies in the Galactic bulge, at 7.2(+0.6,−1.7)kpc from Earth. The accurate measurement of the measured planet-to-host star mass ratio will be included in the next statistical analysis of cold planet demography detected by microlensing.

Clément Ranc↗

Optimization of Microphone Locations for Acoustic Liner Impedance Eduction

Two impedance eduction methods are explored for use with data acquired in the NASA Langley Grazing Flow Impedance Tube. The first is an indirect method based on the convected Helmholtz equation, and the second is a direct method based on the Kumaresan and Tufts algorithm. Synthesized no-flow data, with random jitter to represent measurement error, are used to evaluate a number of possible microphone locations. Statistical approaches are used to evaluate the suitability of each set of microphone locations. Given the computational resources required, small sample statistics are employed for the indirect method. Since the direct method is much less computationally intensive, a Monte Carlo approach is employed to gather its statistics. A comparison of results achieved with full and reduced sets of microphone locations is used to determine which sets of microphone locations are acceptable. For the indirect method, each array that includes microphones in all three regions (upstream and downstream hard wall sections, and liner test section) provides acceptable results, even when as few as eight microphones are employed. The best arrays employ microphones well away from the leading and trailing edges of the liner. The direct method is constrained to use microphones opposite the liner. Although a number of arrays are acceptable, the optimum set employs 14 microphones positioned well away from the leading and trailing edges of the liner. The selected sets of microphone locations are also evaluated with data measured for ceramic tubular and perforate-over-honeycomb liners at three flow conditions (Mach 0.0, 0.3, and 0.5). They compare favorably with results attained using all 53 microphone locations. Although different optimum microphone locations are selected for the two impedance eduction methods, there is significant overlap. Thus, the union of these two microphone arrays is preferred, as it supports usage of both methods. This array contains 3 microphones in the upstream hard wall section, 14 microphones opposite the liner, and 3 microphones in the downstream hard wall section.

Jones, M. G.↗

Orbital Debris: Past, Present, and Future

In the early days of spaceflight, the gBig Sky h theory was the near universally accepted paradigm for dealing with collisions of orbiting objects. This theory was also used during the early years of the aviation industry. Just as it did in aviation, the gBig Sky h theory breaks down as more and more objects accumulate in the environment. Fortunately, by the late 1970 fs some visionaries in NASA and the US Department of Defense (DoD) realized that trends in the orbital environment would inevitably lead to increased risks to operational spacecraft from collisions with other orbiting objects. The NASA Orbital Debris Program was established at and has been conducted at Johnson Space Center since 1979. At the start of 1979, fewer than 5000 objects were being tracked by the US Space Surveillance Network and very few attempts had been made to sample the environment for smaller sizes. Today, the number of tracked objects has quadrupled. Ground ]based and in situ measurements have statistically sampled the LEO environment over most sizes and mitigation guidelines and requirements are common among most space faring nations. NASA has been a leader, not only in defining the debris environment, but in promoting awareness of the issues in the US and internationally, and in providing leadership in developing policies to address the issue. This paper will discuss in broad terms the evolution of the NASA debris program from its beginnings to its present broad range of debris related research. The paper will discuss in some detail current research topics and will attempt to predict future research trends.

Stansbery, Gene↗

Statistical interpretation of the Amerithrax “morph” assay results

In the Amerithrax investigation PCR–based “morph assays” were used to link the anthrax letters with the RMR–1029 flask at USAMRIID. Quantitative data reported for several of these assays are not consistent with Poisson sampling statistics, but instead exhibit “Taylor's Law” behavior where the variance greatly exceeds the mean. A plausible statistical model for this behavior can explain the large number of observed negative and “inconclusive” findings, and implies a high likelihood that a repository sample could contain a “morph” mutant at concentrations well above the nominal detection limit but nonetheless give a negative or inconclusive test result. A Bayesian framework relates the assay results to the probability that a sample actually contains all four morph mutants, even though it tested negative for at least one. The analysis implies that the observed false negative rate actually does not significantly weaken the conclusion that the morph assays correctly excluded all but the stocks derived from RMR–1029 as possible sources of the letter powders, at least when the test results were unambiguous. In conclusion, these findings expand upon and resolve some of the issues cited in recent reviews, and indicate the importance of developing a rigorous statistical framework for interpreting “morph” assay data.

59 BASIC BIOLOGICAL SCIENCES↗

Implementation and Testing of Turbulence Models for the F18-HARV Simulation

This report presents three methods of implementing the Dryden power spectral density model for atmospheric turbulence. Included are the equations which define the three methods and computer source code written in Advanced Continuous Simulation Language to implement the equations. Time-history plots and sample statistics of simulated turbulence results from executing the code in a test program are also presented. Power spectral densities were computed for sample sequences of turbulence and are plotted for comparison with the Dryden spectra. The three model implementations were installed in a nonlinear six-degree-of-freedom simulation of the High Alpha Research Vehicle airplane. Aircraft simulation responses to turbulence generated with the three implementations are presented as plots.

Yeager, Jessie C.↗

CRN5EXP: Expert system for statistical quality control

The purpose of the Expert System CRN5EXP is to assist in checking the quality of the coils at two very important mills: Hot Rolling and Cold Rolling in a steel plant. The system interprets the statistical quality control charts, diagnoses and predicts the quality of the steel. Measurements of process control variables are recorded in a database and sample statistics such as the mean and the range are computed and plotted on a control chart. The chart is analyzed through patterns using the C Language Integrated Production System (CLIPS) and a forward chaining technique to reach a conclusion about the causes of defects and to take management measures for the improvement of the quality control techniques. The Expert System combines the certainty factors associated with the process control variables to predict the quality of the steel. The paper presents the approach to extract data from the database, the reason to combine certainty factors, the architecture and the use of the Expert System. However, the interpretation of control charts patterns requires the human expert's knowledge and lends to Expert Systems rules.

Hentea, Mariana↗

Neutron/Gamma-ray discrimination through measures of fit

Statistical tests and their underlying measures of fit can be utilized to separate neutron/gamma-ray pulses in a mixed radiation field. In this article, first the application of a sample statistical test is explained. Fit measurement-based methods require true pulse shapes to be used as reference for discrimination. This requirement makes practical implementation of these methods difficult; typically another discrimination approach should be employed to capture samples of neutrons and gamma-rays before running the fit-based technique. In this article, we also propose a technique to eliminate this requirement. These approaches are applied to several sets of mixed neutron and gamma-ray pulses obtained through different digitizers using stilbene scintillator in order to analyze them and measure their discrimination quality. (authors)

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Analysis of the Diurnal Cycle of Precipitation and its Relation to Cloud Radiative Forcing using TRMM Products

By incorporating the Tropical Rain Measuring Mission (TRMM) satellite orbital information into the geodesic version of the Colorado State University General Circulation Model (CSU GCM), we are able to fly a satellite in the GCM, and sample the simulated atmosphere in the same way as the TRMM sensors sample the real atmosphere. The TRMM sampling statistics of precipitation and radiative fluxes at annual, intraseasonal, monthly-mean and composited diurnal time scales are evaluated by comparing the satellite-sampled against fully-sampled simulated atmospheres. This information provides a valuable guidance for efficient usage of TRMM data and future satellite mission planning.

Randall, David A.↗

Measurements of three-flavor neutrino oscillations from a PISCES two-detector fit to the NOvA Experiment data

NOvA is a long-baseline neutrino oscillation experiment with two functionally identical detectors: a Near Detector (ND) at Fermilab, placed 1 km from the neutrino source, and a Far Detector (FD) located 810 km away from the ND in Minnesota. NOvA s primary physics goals are to measure the neutrino oscillation parameters $\theta_{23}$ and $\Delta m^2_{32}$ with high precision, determine the neutrino mass hierarchy, and constrain the value of $\delta_{CP}$, primarily via the study of muon neutrino to electron neutrino oscillation. Extracting values for oscillation parameters from fits to data usually relies on treating systematic uncertainties as nuisance parameters, a strategy that suffers from poor scalability as the number of uncertainties becomes larger. This work introduces PISCES (Parameter Inference with Systematic Covariance and Exact Statistics), a novel method that circumvents this scalability problem by encoding systematic uncertainties into a covariance matrix. PISCES utilizes a nested minimization in which optimal systematic pulls are first computed using the covariance matrix in an inner minimization step, then the oscillation parameters are profiled over in the outer minimization. PISCES also uses a Poisson Likelihood term, making it ideal for the inclusion of low-statistic samples in the fits. PISCES is a flexible framework that also supports complex fits, such as a joint Near and Far detector fit. In the standard NOvA analysis, oscillation parameters are extracted using an extrapolation technique in which the ND data indirectly constrain the FD prediction via a ratio method. PISCES, on the other hand, enables a simultaneous ND+FD fit, allowing the high-statistics ND data to directly constrain systematic uncertainties across all samples. This thesis presents the full PISCES joint ND+FD fit for the NOvA three-flavor analysis, details its implementation, and evaluates its performance through extensive robustness tests and fake data studies. It also provides a comparison between the PISCES joint ND+FD results and the standard NOvA extrapolation method using the full NOvA 10-year data set. The results demonstrate that PISCES can successfully fit NOvA data while incorporating the constraints from the ND detectors consistently, using physically motivated systematic uncertainties to account for data/MC discrepancies.

Rajaoalisoa, Miriama [Cincinnati U.]↗

Statistics of the geomagnetic secular variation for the past 5Ma

A new statistical model is proposed for the geomagnetic secular variation over the past 5Ma. Unlike previous models, the model makes use of statistical characteristics of the present day geomagnetic field. The spatial power spectrum of the non-dipole field is consistent with a white source near the core-mantle boundary with Gaussian distribution. After a suitable scaling, the spherical harmonic coefficients may be regarded as statistical samples from a single giant Gaussian process; this is the model of the non-dipole field. The model can be combined with an arbitrary statistical description of the dipole and probability density functions and cumulative distribution functions can be computed for declination and inclination that would be observed at any site on Earth's surface. Global paleomagnetic data spanning the past 5Ma are used to constrain the statistics of the dipole part of the field. A simple model is found to be consistent with the available data. An advantage of specifying the model in terms of the spherical harmonic coefficients is that it is a complete statistical description of the geomagnetic field, enabling us to test specific properties for a general description. Both intensity and directional data distributions may be tested to see if they satisfy the expected model distributions.

Constable, C. G.↗

Statistics of the geomagnetic secular variation for the past 5 m.y

A new statistical model is proposed for the geomagnetic secular variation over the past 5Ma. Unlike previous models, the model makes use of statistical characteristics of the present day geomagnetic field. The spatial power spectrum of the non-dipole field is consistent with a white source near the core-mantle boundary with Gaussian distribution. After a suitable scaling, the spherical harmonic coefficients may be regarded as statistical samples from a single giant Gaussian process; this is the model of the non-dipole field. The model can be combined with an arbitrary statistical description of the dipole and probability density functions and cumulative distribution functions can be computed for declination and inclination that would be observed at any site on Earth's surface. Global paleomagnetic data spanning the past 5Ma are used to constrain the statistics of the dipole part of the field. A simple model is found to be consistent with the available data. An advantage of specifying the model in terms of the spherical harmonic coefficients is that it is a complete statistical description of the geomagnetic field, enabling us to test specific properties for a general description. Both intensity and directional data distributions may be tested to see if they satisfy the expected model distributions.

Constable, C. G.↗

Flying the TRMM Satellite in a GCM

By incorporating the Tropical Rainfall Measurement Mission (TRMM) satellite orbital information into the Colorado State University General Circulation Model (CSU GCM), we are able to 'fly' a satellite in the GCM, and sample the simulated atmosphere in the same way as the TRMM sensors sample the real atmosphere. The TRMM-sampled statistics for precipitation and radiative fluxes at annual, intraseasonal, monthly-mean and seasonal-mean diurnal time scales are evaluated by comparing the satellite-sampled against fully-sampled simulated atmospheres. The sampling rates of the TRMM sensors are significantly affected by the sensors' swath widths. The TRMM Microwave Imager (TMI) and the Visible Infrared Scanner (VIRS) sample each 2.25 x 2.25 degree grid box in the Tropics and subtropics about once per day, but at a different local time every day, while the Precipitation Radar (PR) and the Clouds and the Earth's Radiant Energy System (CERES) sensor visit each grid box about once every three days and twice per day, respectively. Besides inadequate samplings resulting from sensors' swath widths, there is a large, systematic diurnal undersampling associated with TRMM's orbital geometry for grid boxes away from the Equator. When only one month of TRMM data are used, this diurnal undersampling can lead to more daytime samples relative to nighttime samples in one hemisphere, and more nighttime samples relative to daytime samples in the other hemisphere. The resulting sampling biases (3-6 W m(exp-2)) are very pronounced in outgoing longwave radiation (OLR) over the subtropical land masses. The sampling errors in OLR monthly- and seasonal-means are less than 8 W m(exp-2) (5%) for each 2.25 x 2.25 degree grid box. The OLR monthly- and seasonal-means are not sensitive to diurnal undersamplings associated with the TRMM orbits and sensors' swath widths. However, this is not the case for total precipitation. Diurnal undersampling could produce errors as large as 20% in the Tropics and 40% in the subtropics, for the zonally averaged monthly-mean rain rates. The TRMM orbits sample each 2.25 x 2.25 degree grid box in the Tropics and subtropics 1-6 times for each hour of the day within a single season. The seasonal-mean diurnal cycles of precipitation and OLR are not well sampled for any one grid box. By either accumulating the satellite data for a long enough period, or averaging the data over a large area with a relatively uniform diurnal signal, the diurnal cycles of precipitation and OLR can be satisfactorily sampled. The effects TRMM sampling errors on the inferred tropical-mean hydrologic cycle and radiative fluxes are also evaluated. There are strong spurious oscillations associated with TRMM's orbital geometry, with periods of 23 days and 3-4 months, in tropical-mean daily and monthly precipitation. While the relative fluctuations of the sampled-OLR are negligible, the relative fluctuations of the sampled precipitation have magnitudes similar to those of the observed climate variability. Caution must therefore be used when applying TRMM observations of tropical-mean precipitation to interpret climate variations at intraseasonal and interannual scales.

Lin, Xin↗

Use of TRMM Data to Test an Improved Parameterization of Stratiform Precipitation

By incorporating the TRMM satellite orbital information into the most recent geodesic version of the Colorado State University General Circulation Model (CSU GCM), we are able to fly a satellite in the GCM, and sample the simulated atmosphere in the same way as the TRMM sensors sample the real atmosphere. The TRMM sampling statistics of precipitation and radiative fluxes at annual, intra-seasonal, monthly-mean and composited diurnal time scales are evaluated by comparing the satellite-sampled against fully-sampled simulated atmospheres. This information provides a valuable guidance for efficient usage of TRMM data and future satellite mission planning. We have also evaluated the effects of TRMM sampling errors on the inferred tropical-mean hydrologic cycle and radiative fluxes. We have found that there are strong spurious oscillations associated with the TRMM orbital geometry, with periods of 23 days and 3-4 months, in tropical-mean daily and monthly precipitation. Caution must therefore be used when applying TRMM observations of tropical-mean precipitation to interpret climate variations at intraseasonal and interannual scales.

Randall, David A.↗

Robust parallel iterative solvers for linear and least-squares problems (Final Technical Report)

The primary goal of this project is to study and develop robust iterative methods for solving linear systems of equations and least squares systems. The focus of the Minnesota team is on algorithms development, robustness issues, and on tests and validation of the methods on realistic problems. 1. The project begun with an investigation on how to practically update a preconditioner obtained from an ILU-type factorization, when the coefficient matrix changes. 2. We investigated strategies to improve robustness in parallel preconditioners in a specific case of a PDE with discontinuous coefficients. 3. We explored ways to adapt standard preconditioners for solving linear systems arising from the Helmholtz equation. These are often difficult linear systems to solve by iterative methods. 4. We have also worked on purely theoretical issues related to the analysis of Krylov subspace methods for linear systems. 5. We developed an effective strategy for performing ILU factorizations for the case when the matrix is highly indefinite. The strategy uses shifting in some optimal way. The method was extended to the solution of Helmholtz equations by using complex shifts, yielding very good results in many cases. 6. We addressed the difficult problem of preconditioning sparse systems of equations on GPUs. 7. A by-product of the above work is a software package consisting of an iterative solver library for GPUs based on CUDA. This was made publicly available. It was the first such library that offers complete iterative solvers for GPUs. 8. We considered another form of ILU which blends coarsening techniques from Multigrid with algebraic multilevel methods. 9. We have released a new version on our parallel solver - called pARMS [new version is version 3]. As part of this we have tested the code in complex settings - including the solution of Maxwell and Helmholtz equations and for a problem of crystal growth.10. As an application of polynomial preconditioning we considered the problem of evaluating f(A)v which arises in statistical sampling. 11. As an application to the methods we developed, we tackled the problem of computing the diagonal of the inverse of a matrix. This arises in statistical applications as well as in many applications in physics. We explored probing methods as well as domain-decomposition type methods. 12. A collaboration with researchers from Toulouse, France, considered the important problem of computing the Schur complement in a domain-decomposition approach. 13. We explored new ways of preconditioning linear systems, based on low-rank approximations.

97 MATHEMATICS AND COMPUTING↗