Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “applied 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 901 records · Page 50

Trajectory Design for a Single-String Impactor Concept

This paper introduces a trajectory design for a secondary spacecraft concept to augment science return in interplanetary missions. The concept consist of a single-string probe with a kinetic impactor on board that generates an artificial plume to perform in-situ sampling. The trajectory design was applied to a particular case study that samples ejecta particles from the Jovian moon Europa. Results were validated using statistical analysis. Details regarding the navigation, targeting and disposal challenges related to this concept are presented herein.

Monte Carlo↗

Thermodynamic Spectrum of Solar Flares Based on SDO/EVE Observations: Techniques and First Results

The Solar Dynamics Observatory (SDO)/EUV Variability Experiment (EVE) provides rich information on the thermodynamic processes of solar activities, particularly on solar flares. Here, we develop a method to construct thermodynamic spectrum (TDS) charts based on the EVE spectral lines. This tool could potentially be useful for extreme ultraviolet (EUV) astronomy to learn about the eruptive activities on distant astronomical objects. Through several cases, we illustrate what we can learn from the TDS charts. Furthermore, we apply the TDS method to 74 flares equal to or greater than the M5.0 class, and reach the following statistical results. First, EUV peaks are always behind the soft X-ray (SXR) peaks and stronger flares tend to have faster cooling rates. There is a power-law correlation between the peak delay times and the cooling rates, suggesting a coherent cooling process of flares from SXR to EUV emissions. Second, there are two distinct temperature drift patterns, called Type I and Type II. For Type I flares, the enhanced emission drifts from high to low temperature like a quadrilateral, whereas for Type II flares the drift pattern looks like a triangle. Statistical analysis suggests that Type II flares are more impulsive than Type I flares. Third, for late-phase flares, the peak intensity ratio of the late phase to the main phase is roughly correlated with the flare class, and the flares with a strong late phase are all confined. We believe that the re-deposition of the energy carried by a flux rope, which unsuccessfully erupts out, into thermal emissions is responsible for the strong late phase found in a confined flare. Furthermore, we show the signatures of the flare thermodynamic process in the chromosphere and transition region in the TDS charts. These results provide new clues to advance our understanding of the thermodynamic processes of solar flares and associated solar eruptions, e.g., coronal mass ejections.

Wang, Yuming↗

Examining Dense Data Usage near the Regions with Severe Storms in All-Sky Microwave Radiance Data Assimilation and Impacts on GEOS Hurricane Analyses

Many numerical weather prediction (NWP) centers assimilate radiances affected by clouds and precipitation from microwave sensors, with the expectation that these data can provide critical constraints on meteorological parameters in dynamically sensitive regions to make significant impacts on forecast accuracy for precipitation. The Global Modeling and Assimilation Office (GMAO) at NASA Goddard Space Flight Center assimilates all-sky microwave radiance data from various microwave sensors such as all-sky GPM Microwave Imager (GMI) radiance in the Goddard Earth Observing System (GEOS) atmospheric data assimilation system (ADAS), which includes the GEOS atmospheric model, the Gridpoint Statistical Interpolation (GSI) atmospheric analysis system, and the Goddard Aerosol Assimilation System (GAAS). So far, most of NWP centers apply same large data thinning distances, that are used in clear-sky radiance data to avoid correlated observation errors, to all-sky microwave radiance data. For example, NASA GMAO is applying 145 km thinning distances for most of satellite radiance data including microwave radiance data in which all-sky approach is implemented. Even with these coarse observation data usage in all-sky assimilation approach, noticeable positive impacts from all-sky microwave data on hurricane track forecasts were identified in GEOS-5 system. The motivation of this study is based on the dynamic thinning distance method developed in our all-sky framework to use of denser data in cloudy and precipitating regions due to relatively small spatial correlations of observation errors. To investigate the benefits of all-sky microwave radiance on hurricane forecasts, several hurricane cases selected between 2016-2017 are examined. The dynamic thinning distance method is utilized in our all-sky approach to understand the sources and mechanisms to explain the benefits of all-sky microwave radiance data from various microwave radiance sensors like Advanced Microwave Sounder Unit (AMSU-A), Microwave Humidity Sounder (MHS), and GMI on GEOS-5 analyses and forecasts of various hurricanes.

Kim, Min-Jeong↗

Multiscale Studies of the Three-Dimensional Dayside X-Line

We review recent experience from the Cluster, Double Star, and THEMIS missions for lessons that apply to the upcoming Magnetospheric Multiscale Mission (MMS) being developed for launch in 2014. On global scales, simulation and statistical studies lead to mean congurations of dayside reconnection, implying specific relative alignments of the inflow magnetic fields and X-line, with implications for MMS operations designed to maximize the number of close encounters with the diffusion region. At intermediate MHD-to-ion scales, reconstruction of features created by one or two X-lines have developed to the point where data from a cluster of spacecraft can determine their temporal trends and the approximate three-dimensional X-line structure. Recent petascale particle-in-cell (PIC) simulations of reconnection encompass three spatial dimensions with excellent resolution, and make striking predictions of electron scale physics that creates complex interacting flux ropes under component reconnection. High time resolution measurements from MMS will determine the detailed electron scale kinetics embedded within the global and MHDion scale contexts. These developments will lead to the refinement of our three-dimensional multiscale picture of reconnection, yielding improved understanding of the global, MHD, and local physics controlling the onset or quenching, variability, and mean rate of reconnection. This in turn will enable improved predictability of the structural features created by transient reconnection, and their space weather consequences.

Reconnection↗

A statistical method for time-continuous assimilation of remote sounding temperatures

A time-continuous statistical method is presented for the four-dimensional assimilation of remote sounding temperatures based on radiometric measurements from polar-orbiting satellites. This method is applied to DST-6 data from the NOAA-4 and Nimbus-6 satellites. Experiments are reported in which the state of the atmosphere throughout the test period was determined using a varying amount of satellite data and in which different methods were used for their assimilation. Data from the NOAA-4 satellite only, from Nimbus-6 only, and from both satellites together were used; the methods tested include different variations of the statistical method as well as more traditional methods. The conclusions are that: (1) satellite-derived temperature data can have a modest, but statistically significant positive impact on numerical weather prediction in the two-to-three day range; (2) this impact is highly sensitive to the quantity of data available; and (3) the assimilation method plays a major role in the magnitude of the impact for the same data.

Ghil, M.↗

An Evaluation of a Neutron Time Correlated Interrogation Method for Measurement of Fissile Content in Research Reactor Spent Fuel Assemblies

There are many research reactors worldwide that have been, or are being converted from the use of high enrichment uranium (HEU) to low enrichment uranium (LEU, < 20% enriched). The verification of fissile content and initial enrichment of the spent fuel is needed for the effective safeguards of the fuel. The advanced experimental fuel counter (AEFC) was developed for the measurement of spent fuel rods and assemblies from research reactors for safeguards verification. This measurement system contains components for active neutron interrogation, passive neutron totals counting, neutron coincidence counting, and gross gamma-ray counting. This report presents the first application of the time correlated interrogation technique for the measurement of the 235 U content in research reactor spent fuel assemblies. The technique, called time correlated induced fission (TCIF), uses a 252 Cf neutron source to irradiate the fuel assembly, and the subsequent induced fission events in the fissile material are measured by coincidence counting. The doubles rates are enhanced by having the neutron trigger events from both the 252 Cf source and the induced fission neutrons in the same time gate in the coincidence analysis. The average neutrons per fission of the 252 Cf source is 3.76 and the induced fission neutrons for 235 U is 2.44, so the number of neutrons that are produced is higher than for random neutron interrogation. This high effective neutron number increases the multiplicity counting rates and reduces the statistical error. The background coincidence counts from the 252 Cf are reduced by the water in the sample cavity and the polyethylene surrounding the 3 He detector tubes. This method of active neutron interrogation has been applied to the measurement of spent research reactor (IRT) fuel assemblies. The advanced experimental fuel counter (AEFC) was used to compare the TCIF method with the typically used AmLi neutron interrogation source that emits neutrons that are random in time. The statistical uncertainty for the use of the random neutron source (AmLi) and the time correlated source ( 252 Cf) for spent fuel interrogations was evaluated.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Hierarchical Inference with Bayesian Neural Networks: An Application to Strong Gravitational Lensing

In the past few years, approximate Bayesian Neural Networks (BNNs) have demonstrated the ability to produce statistically consistent posteriors on a wide range of inference problems at unprecedented speed and scale. However, any disconnect between training sets and the distribution of real-world objects can introduce bias when BNNs are applied to data. This is a common challenge in astrophysics and cosmology, where the unknown distribution of objects in our universe is often the science goal. In this work, we incorporate BNNs with flexible posterior parameterizations into a hierarchical inference framework that allows for the reconstruction of population hyperparameters and removes the bias introduced by the training distribution. We focus on the challenge of producing posterior PDFs for strong gravitational lens mass model parameters given Hubble Space Telescope–quality single-filter, lens-subtracted, synthetic imaging data. We show that the posterior PDFs are sufficiently accurate (statistically consistent with the truth) across a wide variety of power-law elliptical lens mass distributions. We then apply our approach to test data sets whose lens parameters are drawn from distributions that are drastically different from the training set. We show that our hierarchical inference framework mitigates the bias introduced by an unrepresentative training set's interim prior. Simultaneously, we can precisely reconstruct the population hyperparameters governing our test distributions. Our full pipeline, from training to hierarchical inference on thousands of lenses, can be run in a day. The framework presented here will allow us to efficiently exploit the full constraining power of future ground- and space-based surveys (https://github.com/swagnercarena/ovejero).

79 ASTRONOMY AND ASTROPHYSICS↗

Effective techniques for the identification and accommodation of disturbances

The successful control of dynamic systems such as space stations, or launch vehicles, requires a controller design methodology that acknowledges and addresses the disruptive effects caused by external and internal disturbances that inevitably act on such systems. These disturbances, technically defined as uncontrollable inputs, typically vary with time in an uncertain manner and usually cannot be directly measured in real time. A relatively new non-statistical technique for modeling, and (on-line) identification, of those complex uncertain disturbances that are not as erratic and capricious as random noise is described. This technique applies to multi-input cases and to many of the practical disturbances associated with the control of space stations, or launch vehicles. Then, a collection of smart controller design techniques that allow controlled dynamic systems, with possible multi-input controls, to accommodate (cope with) such disturbances with extraordinary effectiveness are associated. These new smart controllers are designed by non-statistical techniques and typically turn out to be unconventional forms of dynamic linear controllers (compensators) with constant coefficients. The simplicity and reliability of linear, constant coefficient controllers is well-known in the aerospace field.

Johnson, C. D.↗

Trajectory Ensemble Methods Provide Single-Molecule Statistics for Quantum Dynamical Systems

The emergence of experiments capable of probing quantum dynamics at the single-molecule level requires the development of new theoretical tools capable of simulating and analyzing these dynamics beyond an ensemble-averaged description. In this article, we present an efficient method for sampling and simulating the dynamics of the individual quantum systems that make up an ensemble and apply it to study the nonequilibrium dynamics of the ubiquitous spin-boson model. We generate an ensemble of single-system trajectories, and we analyze this trajectory ensemble using tools from classical statistical mechanics. Our results demonstrate that the dynamics of quantum coherence is highly heterogeneous at the single-system level due to variations in the initial bath configuration, which significantly affects the transient exchange of coherence between the system and its bath. Here, we observe that single systems tend to retain coherence over time scales longer than that of the ensemble. We also compute a novel thermodynamic entanglement entropy that quantifies a thermodynamic driving force favoring system–bath entanglement.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

FAST.Farm Development and Validation of Structural Load Prediction Against Large Eddy Simulations

FAST.Farm is a mid-fidelity engineering tool developed by the National Renewable Energy Laboratory targeted at accurately and efficiently predicting wind turbine power production and structural loading in wind farm settings, including wake interactions between turbines. FAST.Farm is based on several principles of the dynamic wake meandering (DWM) model, but also addresses limitations of previous DWM implementations. Previous FAST.Farm studies have shown the similarities and differences between FAST.Farm and large eddy simulations for rigid turbine cases. The objective of this work is to quantify the ability of FAST.Farm to accurately predict turbine structural response in a small wind farm. This is done by comparing FAST.Farm structural response to SOWFA-OpenFAST results for three laterally-aligned turbines. The purpose of this study is to characterize the similarities and differences between FAST.Farm and a higher fidelity model for predicting turbine structural response in a wind farm for differing atmospheric inflows. Strong statistical agreement was found between FAST.Farm and SOWFA-OpenFAST structural response for the non-waked upstream turbine, and good agreement was found for the downstream turbines for most structural quantities. Higher differences were seen for downstream turbines with low ambient turbulence intensity or yawed turbines, suggesting areas for FAST.Farm wake dynamics modeling improvements. For all cases and turbines, small statistical differences were seen between blade deflections and bending moments, with larger differences for tower-top and tower-base bending moments. Overall, the results establish confidence for applying FAST.Farm to wind farm power and loads analyses and identify areas where further model validations and model improvements should be targeted.

49 EE - Wind and Water Power Program - Wind (EE-4W↗

VADER: A Tool for Criticality Safety Validation

The purpose of criticality safety is to prevent any inadvertent criticality from occurring during the handling or storage of fissile material. Calculations are frequently used to demonstrate that a sufficient subcritical margin exists. Validation is a key aspect of the evaluation process, establishing the suitability, accuracy, and associated uncertainty of the computational method and data to be used for the intended application. The validation process is performed by comparing the results of critical experiments with the calculated results from models of the experiments using the computational method to be validated. Laboratory critical experiments are controlled systems that achieve a k eff of approximately 1 in order to investigate the parameters at which such a critical condition is achieved. The validation parameters that are traditionally applied to safety analysis calculations are the bias and the bias uncertainty . The bias is the deviation of the average k eff of the validation suite from unity. The bias uncertainty accounts for the statistical uncertainty in the bias based on the standard deviation, sample size, and distribution of k eff values of the validation suite. The values of bias and bias uncertainty ensure that the systems predicted to be subcritical by the computational method will indeed be subcritical. The bias and bias uncertainty are often combined to determine an upper subcritical limit (USL) or computational margin that can then be applied to safety analysis calculations. Many methods have been developed by different organizations to calculate the bias and bias uncertainty for various types of criticality analyses. Each of these methods typically requires that the validity of various underpinning statistical assumptions be confirmed to demonstrate that the method is appropriate for the analysis of a given validation suite. An example of the validation decision making flow is shown in Fig.1. As shown in Fig. 1, the analyst performing the validation fits a trend line to the data and performs a test to determine if the trend was a statistically better representation of the data than if it were treated as an uncorrelated sample. If the trend line is a better representation of the data, then the analyst uses any one of a number of trending techniques to determine the bias and bias uncertainty. If a trend is not an appropriate representation of the data, then the analyst proceeds to perform a normality assessment for the data. If the normal assumption can be shown to be acceptable, then the analyst calculates the bias and bias uncertainty with the parametric technique. If the assumption of normality cannot be justified, then the nonparametric technique is used. Once the decision flow has been followed and the appropriate technique has been selected, the bias and bias uncertainty is typically combined with an administrative margin to determine a USL below which calculated values of k eff for safety analysis models can be considered subcritical. The calculations used in each decision are often performed with spreadsheets or with small programs available at various sites performing criticality analyses. Expertise in understanding and interpreting the results must be maintained to perform these calculations. This can often be an error-prone process. Oak Ridge National Laboratory (ORNL) is currently developing the Validation and Data Evaluation Resource (VADER) to simplify and automate the criticality safety validation process and to provide a software quality assurance pedigree to the calculational methods used. This paper discusses the use of the Fulcrum user interface with VADER, the anticipated initial capabilities of VADER to perform validation analyses, and the output from the code.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Statistical Analysis of the Limiting Dynamics of Two dimensional Boussinesq Turbulent Systems

In the current study, we investigate rotational and stratified turbulent systems within the framework of the Boussinesq equations. The effects of rotation and stratification give rise to a variety of dynamical regimes in atmospheric and oceanic turbulence. To develop a better understanding of the fundamental flow properties through a more straightforward mathematical framework, we commence our research with a simplified two-dimensional Boussinesq model that incorporates rotation and stratification effects. We explore two distinct limiting dynamics, considering pure rotation and pure stratification separately. We examine these intriguing dynamics through numerical investigations, applying the exact solution theory of Boussinesq equations and examining small-scale perturbations to the exact solution under rotation and stratification. We conduct an analysis of statistical quantities associated with the turbulent state variables. The results indicate that, under weak rotation at the equilibrium statistical state, the flow field exhibits vortex flow. However, strong rotation transforms the flow into vertical shear flow. We observe both downscale and upscale energy transfers with a decay rate of k −3 in the absence of external forcing and dissipation. Furthermore, we delve into the study of internal gravity wave mode interactions in the simplified two-dimensional Boussinesq system, which is a crucial aspect of geophysical turbulence due to gravity. Our findings reveal that the time series of the mode coefficients exhibit wave-like interactions between wave modes.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Interpolating Fields of Carbon Monoxide Data Using a Hybrid Statistical-Physical Model

Atmospheric Carbon Monoxide (CO) is a pollutant gas of which the US congress has mandated regular monitoring, and satellite sensors can be used to retrieve regional concentrations of CO over several vertical layers. However, CO at cloudy locations cannot be observed and have to be estimated from the observed data set, resulting in an interpolation problem. The current state-of-the-art solution is to combine prior information, computed by a deterministic physical model, with observations. However, the deterministic model may introduce uncertainties that do not derive from the data. While sharing certain features with the physical model, this paper presents a Bayesian hierarchical model for interpolating CO on a 3-dimensional spatial grid, across time. To our knowledge such a model has not been considered before. The model is applied to a hypothetical air-quality monitoring scenario, and is compared to existing interpolation methods. The results provide motivation for the use of the statistical model for regional to local applications.

Arellano, A. A.↗

Modern analogs for understanding pollen-vegetation dynamics in a Mediterranean mosaic landscape (Balearic Islands, Western Mediterranean)

This paper presents the study of modern pollen analogs from the Balearic Islands. While similar studies have been largely applied to mainland areas, research focused on modern vegetation dynamics on Mediterranean islands remains very rare. In this research, we combine vegetation surveys, pollen analysis and multivariate statistics to understand landscape composition. The main objectives of are: (1) to examine pollen-vegetation relationships in relation to environmental and land-use variables; (2) to understand modern pollen representation in a mosaic landscape structure; and (3) to propose pollen indicators that characterize the primary vegetation types from the Balearic Islands to better interpret past pollen records in Mediterranean island environments. Pollen results and Redundancy Analysis (RDA) distinguish three major groups: (a) Holm oak and box formations; (b) maquis and garrigues; and (c) anthropogenic and open habitats. Landscape form, mean decadal rainfall, mean decadal temperature, fire activity, trampling, slope percentage, wet/flooded soil, saline soil, distance to agropastoral cells, gHM index, domestic herbivory presence, agropastoral use, and soil type are the major variables explaining modern pollen assemblage variation in our research. Poaceae undiff., Plantago sp., Apiaceae undiff., Cerealia-t, and Cichorieae are highly correlated to human activities but should be interpreted cautiously when occurring in low values. Quercus ilex-t, Hypericum, and Buxus are correlated to humid locations while Pistacia, Pinus, Juniperus-t, and Olea to high mean decadal temperatures. Our study indicates how pollen analysis and multivariate analysis are powerful tools for characterizing the mosaic landscape, with special focus on the main vegetation types of the Balearic Islands.

Geology↗

Statistics of primordial density perturbations from discrete seed masses

The statistics of density perturbations for general distributions of seed masses with arbitrary matter accretion is examined. Formal expressions for the power spectrum, the N-point correlation functions, and the density distribution function are derived. These results are applied to the case of uncorrelated seed masses, and power spectra are derived for accretion of both hot and cold dark matter plus baryons. The reduced moments (cumulants) of the density distribution are computed and used to obtain a series expansion for the density distribution function. Analytic results are obtained for the density distribution function in the case of a distribution of seed masses with a spherical top-hat accretion pattern. More generally, the formalism makes it possible to give a complete characterization of the statistical properties of any random field generated from a discrete linear superposition of kernels. In particular, the results can be applied to density fields derived by smoothing a discrete set of points with a window function.

Scherrer, Robert J.↗

A New Normalized Difference Cloud Retrieval Technique Applied to Landsat Radiances Over the Oklahoma ARM Site

We suggest a new approach to cloud retrieval, using a normalized difference of nadir reflectivities (NDNR) constructed from a non-absorbing and absorbing (with respect to liquid water) wavelength. Using Monte Carlo simulations we show that this quantity has the potential of removing first order scattering effects caused by cloud side illumination and shadowing at oblique Sun angles. Application of the technique to TM (Thematic Mapper) radiance observations from Landsat-5 over the Southern Great Plains site of the ARM (Atmospheric Radiation Measurement) program gives very similar regional statistics and histograms, but significant differences at the pixel level. NDNR can be also combined with the inverse NIPA (Nonlocal Independent Pixel Approximation) of Marshak (1998) which is applied for the first time on overcast Landsat scene subscenes. We demonstrate the sensitivity of the NIPA-retrieved cloud fields on the parameters of the method and discuss practical issues related to the optimal choice of these parameters.

Orepoulos, Lazaros↗

Automated SEM Modal Analysis Applied to the Diogenites

Analysis of volume proportions of minerals, or modal analysis, is routinely accomplished by point counting on an optical microscope, but the process, particularly on brecciated samples such as the diogenite meteorites, is tedious and prone to error by misidentification of very small fragments, which may make up a significant volume of the sample. Precise volume percentage data can be gathered on a scanning electron microscope (SEM) utilizing digital imaging and an energy dispersive spectrometer (EDS). This form of automated phase analysis reduces error, and at the same time provides more information than could be gathered using simple point counting alone, such as particle morphology statistics and chemical analyses. We have previously studied major, minor, and trace-element chemistry of orthopyroxene from a suite of diogenites. This abstract describes the method applied to determine the modes on this same suite of meteorites and the results of that research. The modal abundances thus determined add additional information on the petrogenesis of the diogenites. In addition, low-abundance phases such as spinels were located for further analysis by this method.

Bowman, L. E.↗