Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “statistical model”

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,027 records · Page 57

Improving ChemCam LIBS long-distance elemental compositions using empirical abundance trends

The ChemCam instrument on the Curiosity rover provides chemical compositions of Martian rocks and soils using remote laser-induced breakdown spectroscopy (LIBS). The elemental calibration is stable as a function of distance for Ti, Fe, Mg, and Ca. The calibration shows small, systematically increasing abundance trends as a function of distance for Al, Na, K, and to some extent, Si. The distance effect is known to be due to a dependence with distance on the relative strengths of atomic transition lines. Emission lines representing transitions from relatively low energy levels remain intense at longer distances while emission lines representing transitions from higher energy levels decrease in intensity more rapidly as a function of distance. The multivariate algorithms used to determine elemental compositions rely on a large number of emission lines in many cases, so rather than trying to correct all emission lines, a study was made of the predicted compositions as a function of distance, in order to determine an empirical correction. Abundance trends can be well approximated by a linear trend with distance within the ranges of abundances and distances observed up to ~6 m. Data from 11 distinct geological members and data groups of the Murray formation in Gale crater, Mars, were used to form the model, selecting the members and data groups yielding the best statistics. The model was tested using data from several targets observed from two different distances, and using data from the Kimberley formation, the composition of which is significantly different from the Murray formation, showing that the model works on other compositions beyond those used to build the model. For long-distance observations up to ~6 m, corrections can be made back to an equivalent composition at the median distance of ChemCam observations (2.6 m). Finally, the model has been validated up to 6.2 m, although ChemCam is able to observe bedrock targets to >7 m, and iron meteorites to distances of >9 m.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Computationally efficient Bayesian estimation of graphical networks for omics data

Graphical networks are useful, widely-used modeling approaches to represent complex biological processes with biological measurements generated by platforms such as mass spectrometry. Bayesian analyses of graphical networks for omics data have several advantages over their frequentist counterparts, such as the inclusion of prior knowledge in the estimation of models. However, Bayesian approaches to date have only been feasible for data with a couple hundred biomolecules due to prohibitive computational time, but omics data often contains tens of thousands of biomolecules. Here, we present and illustrate a more computationally efficient approach named BPlane (Bayesian PseudoLikelihood-based Algorithm for Network Estimation) to extend Bayesian modeling capabilities for larger-sized datasets, such as most untargeted proteomics data. Via simulation, we demonstrate that BPlane produces substantial computational savings over a current state-of-the-art Bayesian algorithm while maintaining competitive edge detection accuracy. On a SARS-CoV2 proteomics data with 7000 proteins, the competing algorithm takes three times as long to complete the first iteration as BPlane takes to converge after over 100 iterations.

EM algorithm↗

Gas Hydrate Stability at Low Temperatures and High Pressures with Applications to Mars and Europa

Gas hydrates are implicated in the geochemical evolution of both Mars and Europa [1- 3]. Most models developed for gas hydrate chemistry are based on the statistical thermodynamic model of van der Waals and Platteeuw [4] with subsequent modifications [5-8]. None of these models are, however, state-of-the-art with respect to gas hydrate/electrolyte interactions, which is particularly important for planetary applications where solution chemistry may be very different from terrestrial seawater. The objectives of this work were to add gas (carbon dioxide and methane) hydrate chemistries into an electrolyte model parameterized for low temperatures and high pressures (the FREZCHEM model) and use the model to examine controls on gas hydrate chemistries for Mars and Europa.

Marion, G. M.↗

Evaluating the Contribution of NASA Remotely-Sensed Data Sets on a Convection-Allowing Forecast Model

The Short-term Prediction Research and Transition (SPoRT) Center is a collaborative partnership between NASA and operational forecasting partners, including a number of National Weather Service forecast offices. SPoRT provides real-time NASA products and capabilities to help its partners address specific operational forecast challenges. One challenge that forecasters face is using guidance from local and regional deterministic numerical models configured at convection-allowing resolution to help assess a variety of mesoscale/convective-scale phenomena such as sea-breezes, local wind circulations, and mesoscale convective weather potential on a given day. While guidance from convection-allowing models has proven valuable in many circumstances, the potential exists for model improvements by incorporating more representative land-water surface datasets, and by assimilating retrieved temperature and moisture profiles from hyper-spectral sounders. In order to help increase the accuracy of deterministic convection-allowing models, SPoRT produces real-time, 4-km CONUS forecasts using a configuration of the Weather Research and Forecasting (WRF) model (hereafter SPoRT-WRF) that includes unique NASA products and capabilities including 4-km resolution soil initialization data from the Land Information System (LIS), 2-km resolution SPoRT SST composites over oceans and large water bodies, high-resolution real-time Green Vegetation Fraction (GVF) composites derived from the Moderate-resolution Imaging Spectroradiometer (MODIS) instrument, and retrieved temperature and moisture profiles from the Atmospheric Infrared Sounder (AIRS) and Infrared Atmospheric Sounding Interferometer (IASI). NCAR's Model Evaluation Tools (MET) verification package is used to generate statistics of model performance compared to in situ observations and rainfall analyses for three months during the summer of 2012 (June-August). Detailed analyses of specific severe weather outbreaks during the summer will be presented to assess the potential added-value of the SPoRT datasets and data assimilation methodology compared to a WRF configuration without the unique datasets and data assimilation.

Zavodsky, Bradley T.↗

Dynamics Explorer guest investigator

A data base of satellite particle, electric field, image, and plasma data was used to determine correlations between the fields and the particle auroral boundaries. A data base of 8 days of excellent coverage from all instruments was completed. The geomagnetic conditions associated with each of the selected data periods, the number of UV image passes per study day that were obtained, and the total number of UV images for each day are given in tabular form. For each of the days listed in Table 1, both Vector Electric Field Instrument (VEFI) electric potential data and LAPI integrated particle energy fluxes were obtained. On the average, between 8 and 11 passes of useful data per day were obtained. These data are displayed in a format such that either the statistical electric field model potential or the statistical precipitation energy flux could be superimposed. The Heppner and Maynard (1987) and Hardy et al. (1987) models were used for the electric potential and precipitation, respectively. In addition, the auroral image intensity along the Dynamics Explorer-2 satellite pass could be computed and plotted along with the LAPI precipitation data and Hardy et al. (1987) values.

Sojka, Jan J.↗

Incorporating Yearly Derived Winter Wheat Maps Into Winter Wheat Yield Forecasting Model

Wheat is one of the most important cereal crops in the world. Timely and accurate forecast of wheat yield and production at global scale is vital in implementing food security policy. Becker-Reshef et al. (2010) developed a generalized empirical model for forecasting winter wheat production using remote sensing data and official statistics. This model was implemented using static wheat maps. In this paper, we analyze the impact of incorporating yearly wheat masks into the forecasting model. We propose a new approach of producing in season winter wheat maps exploiting satellite data and official statistics on crop area only. Validation on independent data showed that the proposed approach reached 6% to 23% of omission error and 10% to 16% of commission error when mapping winter wheat 2-3 months before harvest. In general, we found a limited impact of using yearly winter wheat masks over a static mask for the study regions.

MODIS↗

Experimental validation of a high fidelity Monte Carlo neutron transport model of the MIT graphite exponential pile

High-fidelity modeling and simulation were performed for the MIT graphite exponential pile (MGEP) using Monte Carlo neutron transport codes OpenMC and MCNP, and the results were validated by experimental data. The MGEP is being used as the test bed for the design of an autonomous control system for the pile's neutron flux distribution. The main contribution of this work is to generate the training data sets of neutron flux distributions with different locations of control rods that perturb the neutron flux profiles. First, code -to-code cross verification between OpenMC and MCNP was performed to ensure consistency of the numerical modeling within statistical uncertainties. To validate the accuracy of this high-fidelity model, a series of neutron flux measurements were conducted using a Helium-3 (He-3) neutron detector on a mobile platform that is placed inside the pile. Second, the neutron flux profiles were measured in four vertical layers of interest, and compared to the corresponding simulation results. The comparison results shows that the root mean square error is less than 2.5% in the two upper layers, and less than 4.5% in all four measured layers. Here the results validated the accuracy of the modeling and simulation. Finally, the relative change of the neutron flux profiles from moving control rods was analyzed, which identified the layer that has the best sensitivity regarding the control rods movements. Thus, this work identified and provided training data sets of both simulated and experimental neutron flux profiles in the most sensitive layer, paving the path forward to the real-time experimental demonstration of the autonomous control system.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Ripening of Rh Nanoparticle Catalysts in Reverse Water–Gas Shift via a Data-Driven Model Combining Physics, Theory, and Experiment

Degradation via sintering is an ongoing challenge that impedes the broad commercial success of supported metallic nanoparticle catalysts. To mitigate degradation via informed catalyst design and process operations, here we aim to disambiguate the underlying mechanisms of sintering by combining theory and experiment in a quantitative framework. While mechanistic sintering models exist, they only model a single sintering pathway, even though multiple sintering mechanisms can occur simultaneously or dominate at different stages of the process. Data-driven machine learning models have emerged as a means to represent complex processes through data regression. However, machine learning models have very large data needs and lack mechanistic insights due to their black-box encoding. To develop an interpretive model of catalyst degradation via sintering, we constructed a hybrid model combining mechanistic “physics-based” models and data-driven methods to obtain both reliable predictions and mechanistic insights regarding experimentally observed sintering phenomena. Focusing on nanoparticle sintering in the Rh–TiO 2 catalyst for the reverse water–gas shift (RWGS) reaction, the hybrid model couples a mechanistic term for Ostwald ripening with energy values calculated via density functional theory (DFT) with a parametric, data-driven discrepancy function term for unmodeled mechanisms. The hybrid model is trained using Bayesian inference with data collected from small-angle X-ray scattering (SAXS) in situ experiments wherein average nanoparticle diameter versus time was measured at three relevant operating temperatures. The calibrated hybrid model results show that an Ostwald ripening-only model parameterized with fixed DFT energies does not fully capture the time and temperature dependence of the SAXS-observed sintering kinetics, and that an additional functional contribution, or DFT energy calibration, is required to reconcile simulation and experiment. Analysis of the hybrid-model error confirms that the hybrid model outperforms both the purely mechanistic and purely data-driven alternatives in terms of expected predictive accuracy for time-evolving average particle sizes. Furthermore, the results support the hypothesis that the Ostwald ripening mechanism is less important for explaining the sintering phenomena as operating temperature increases under an assumed fixed DFT parameterization. This could be explained in one of two ways: either latent, unmodeled sintering mechanisms dominate at higher temperatures, or the DFT uncertainty increases with temperature. The proposed modeling approach directly links theory to experiments and simulations via a statistical hybrid modeling framework and can be extended to other catalytic systems to improve predictive models and mechanistic understanding.

Bayesian hybrid modeling↗

Does the emission measure decrease during the start of a soft X-ray flare

It is suggested that the apparent emission dip found in soft X-ray flare observations, interpreted as the emission from a single temperature plasma, is an artifact. According to this hypothesis, the flare phenomenon represents a combination of flare plus a stable active region plasma. A small amount of hot flare plasma will dominate the total emission, and the calculated isothermal plasma parameters will trend towards those of the hot flare plasma, so that a decrease in emission measure will be inferred. It is claimed that observed data fit the model of a stable plasma plus a flare plasma better than it does the model of a single plasma, in accordance with a chi square statistic. The model of stable plasma plus flare plasma is also supported by the observation that the plasma parameters in the postflare condition are indistinguishable from those during preflare.

Wolfson, C. J.↗

Interpreting forest biome productivity and cover utilizing nested scales of image resolution and biogeographical analysis

The objective was to relate spectral imagery of varying resolution with ground-based data on forest productivity and cover, and to create models to predict regional estimates of forest productivity and cover with a quantifiable degree of accuracy. A three stage approach was outlined. In the first stage, a model was developed relating forest cover or productivity to TM surface reflectance values (TM/FOREST models). The TM/FOREST models were more accurate when biogeographic information regarding the landscape was either used to stratigy the landscape into more homogeneous units or incorporated directly into the TM/FOREST model. In the second stage, AVHRR/FOREST models that predicted forest cover and productivity on the basis of AVHRR band values were developed. The AVHRR/FOREST models had statistical properties similar to or better than those of the TM/FOREST models. In the third stage, the regional predictions were compared with the independent U.S. Forest Service (USFS) data. To do this regional forest cover and forest productivity maps were created using AVHRR scenes and the AVHRR/FOREST models. From the maps the county values of forest productivity and cover were calculated. It is apparent that the landscape has a strong influence on the success of the approach. An approach of using nested scales of imagery in conjunction with ground-based data can be successful in generating regional estimates of variables that are functionally related to some variable a sensor can detect.

Iverson, Louis R.↗

Estimation of the Ocean Skin Temperature using the NASA GEOS Atmospheric Data Assimilation System

This report documents the status of the development of a sea surface temperature (SST) analysis for the Goddard Earth Observing System (GEOS) Version-5 atmospheric data assimilation system (ADAS). Its implementation is part of the steps being taken toward the development of an integrated earth system analysis. Currently, GEOS-ADAS SST is a bulk ocean temperature (from ocean boundary conditions), and is almost identical to the skin sea surface temperature. Here we describe changes to the atmosphere-ocean interface layer of the GEOS-atmospheric general circulation model (AGCM) to include near surface diurnal warming and cool-skin effects. We also added SST relevant Advanced Very High Resolution Radiometer (AVHRR) observations to the GEOS-ADAS observing system. We provide a detailed description of our analysis of these observations, along with the modifications to the interface between the GEOS atmospheric general circulation model, gridpoint statistical interpolation-based atmospheric analysis and the community radiative transfer model. Our experiments (with and without these changes) show improved assimilation of satellite radiance observations. We obtained a closer fit to withheld, in-situ buoys measuring near-surface SST. Evaluation of forecast skill scores corroborate improvements seen in the observation fits. Along with a discussion of our results, we also include directions for future work.

AGCM↗

Fade Dynamics and Its Evolution: The Other Part of the Acts Rain Prediction Model

The inception of the Advanced Communication Technology Satellite (ACTS) Project has required as similarly advanced statistical mathematical modeling formalism to describe the behavior of the 30/20 GHz links emanating to and from the earth terminals through the deleterious effects of the earth's atmosphere. The resulting ACTS Rain Attenuation Prediction Model has been thoroughly described in [Manning, 1990]. In the present paper, the basic rudiments of this model will be reviewed; Section 1 covers the static or time-independent portion of the model and Section 2 covers the dynamic of time-dependent portion. The results of Section 2 are then applied to a new approximate solution of the famous problem of the time duration tau of a fade of random process below some threshold. This is known as the fade duration. The new approximate solution was published in Russian [Denisenko] and, unfortunately, was never published into English. Hence, this work is restated following [Denisenko] in Section 3 which is immediately applied to the random rain fade process. The results for all five ACTS propagation sites as well as Tampa, FL are then given.

Manning, Robert M.↗

Fade Dynamics and its Evolution: The Other Part of the ACTS Rain Prediction Model

The inception of the Advanced Communication Technology Satellite (ACTS) Project has required a similarly advanced statistical mathematical modeling formalism to describe the behavior of the 30/20 GHz links emanating to and from the earth terminals through the deleterious effects of the earth's atmosphere. The resulting ACTS Rain Attenuation Prediction Model has been thoroughly described in (Manning). In the present paper, the basic rudiments of this model will be reviewed; Section 1 covers the static or time-independent portion of the model and Section 2 covers the dynamic or time-dependent portion. The results of Section 2 are then applied to a new approximate solution of the famous problem of the time duration tau of a fade of a random process below some threshold. This is known as the fade duration. The new approximate solution was published in Russian (Denisenko) and, unfortunately, was never published into English. Hence, this work is restated following (Denisenko) in Section 3 which is immediately applied to the random rain fade process. The results for all five ACTS propagation sites as well as Tampa, FL are then given.

Manning, Robert M.↗

Multi-Frequency Analysis of Simulated versus Observed Variability in Tropospheric Temperature

Studies seeking to identify a human-caused global warming signal generally rely on climate model estimates of the “noise” of intrinsic natural variability. Assessing the reliability of these noise estimates is of critical importance. We evaluate here the statistical significance of differences between climate model and observational natural variability spectra for global-mean mid- to upper-tropospheric temperature (TMT). We use TMT information from satellites and large multimodel ensembles of forced and unforced simulations. Our main goal is to explore the sensitivity of model-versus-data spectral comparisons to a wide range of subjective decisions. These include the choice of satellite and climate model TMT datasets, the method for separating signal and noise, the frequency range considered, and the statistical model used to represent observed natural variability. Of particular interest is the amplitude of the interdecadal noise against which an anthropogenic tropospheric warming signal must be detected. We find that on time scales of 5–20 years, observed TMT variability is (on average) overestimated by the last two generations of climate models participating in the Coupled Model Intercomparison Project. This result is relatively insensitive to different plausible analyst choices, enhancing confidence in previous claims of detectable anthropogenic warming of the troposphere and indicating that these claims may be conservative. A further key finding is that two commonly used statistical models of short-term and long-term memory have deficiencies in their ability to capture the complex shape of observed TMT spectra.

54 ENVIRONMENTAL SCIENCES↗

Assessing concurrent effects of climate change on hydropower supply, electricity demand, and greenhouse gas emissions in the Upper Yangtze River Basin of China

Hydropower importantly provides flexible low-carbon electricity, however, climate change will affect the hydropower system through altering hydrologic regimes while also affecting electricity demands for heating and cooling that hydropower resources serve. This study assesses the effect of climate change on hydropower and electricity demand in the Upper Yangtze River Basin (UYRB) in China on the regional net electric load and greenhouse gas (GHG) emissions. This is accomplished by using climate projections from five global climate models (GCMs) to simultaneously force (1) a physically-based hydrological model and a statistically-based hydropower model to estimate the future generating capacity of 21 large hydropower plants in the UYRB and (2) an empirical electricity demand model accounting for socioeconomic and climatic factors. Under climate change, the projected hydropower generation in the UYRB tends to increase in the 21st century but is far less than the increase in electricity demand, increasing the gap between demand and supply. Future increases in overall electricity demand are driven by GDP growth, but climate change will alter the distribution of the seasonal electricity demand. Climate warming decreases electricity demand for heating in winter and increases electricity demand for cooling in summer, but ultimately increases demand. Meanwhile, there is an increasing mismatch between electricity demand and hydropower supply associated with inter- and intra-annual variations, owing to the temporal climate change and increase in compound climate extremes (droughts and heatwaves). Finally, meeting the gap between supply and demand due to climate change is estimated to contribute 79.0–184.6 and 50.6–316.2 MMT CO 2e /yr of additional GHG emissions by the mid and end of 21st century, respectively.

13 HYDRO ENERGY↗

Beam energy dependence of net- Λ fluctuations measured by the STAR experiment at the BNL Relativistic Heavy Ion Collider

The measurements of particle multiplicity distributions have generated considerable interest in understanding the fluctuations of conserved quantum numbers in the quantum chromodynamics (QCD) hadronization regime, in particular near a possible critical point and near the chemical freeze-out. Net-protons and net-kaons have been used as proxies for the net-baryon number and net-strangeness, respectively. We report the measurement of efficiency- and centrality-bin width-corrected cumulant ratios (C 2 /C 1 , C 3 /C 2 ) of net- Λ distributions, in the context of both strangeness and baryon number conservation, as a function of collision energy, centrality, and rapidity. The results are for Au + Au collisions at five beam energies ( √s NN = 19.6 , 27, 39, 62.4, and 200 GeV) recorded with the Solenoidal Tracker at RHIC (STAR). We compare our results to the Poisson and negative binomial (NBD) expectations, as well as to ultrarelativistic quantum molecular dynamics (UrQMD) and hadron resonance gas (HRG) model predictions. Both NBD and Poisson baselines agree with data within the statistical and systematic uncertainties. UrQMD describes the measured net- Λ C 1 and C 3 at 200 GeV reasonably well but deviates from C 2 , and the deviation increases as a function of collision energy. The ratios of the measured cumulants show no features of critical fluctuations. The chemical freeze-out temperatures extracted from a recent HRG calculation, which was successfully used to describe the net-proton, net-kaon, and net-charge data, indicate Λ freeze-out conditions similar to those of kaons. However, large deviations are found when comparing with temperatures obtained from net-proton fluctuations. The net- Λ cumulants show a weak but finite dependence on the rapidity coverage in the acceptance of the detector, which can be attributed to quantum number conservation.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Coupled fission fragment angular momenta

Here, nuclear fission produces fragments whose spins are coupled to the relative angular motion via angular momentum conservation. It is shown how ensembles of such spins can readily be obtained by either direct microcanonical sampling or by sampling of the associated normal modes of rotation. The resulting distribution of the spin-spin opening angle is illustrated in various three- and two-dimensional scenarios and it is demonstrated how recent mutually conflicting model calculations can be well reproduced with different assumptions about the scission geometry.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Empirical Acoustic Source Model for Chemical Explosions in Air

Chemical explosions generate pressure disturbances in air that radiate as nonlinear shock waves near the source and transition into acoustic waves with distance. Because low-frequency acoustic waves generally travel large distances without significant loss of energy, they are often used for explosion monitoring and yield estimation. However, quantitative relationships between acoustic energy and explosion yields are required for accurate yield estimation. Here, we develop an empirical acoustic source model for chemical explosions from experimental data. The empirical model returns the acoustic pressure waveform for the detonation of 1 kg of trinitrotoluene, which is conventionally used to represent the explosive release of 4.184 MJ of explosion energy. The full-waveform model can be used to predict acoustic signals for an arbitrary yield of a high-explosive detonation based on the standard scaling law and to estimate acoustic energies in a specific frequency range. We evaluate the accuracy of the acoustic source model independently by estimating the yield of other explosive events that are not included in the model development. Finally, statistical characteristics of the model and their implications for the uncertainty quantification of estimated yields are discussed.

58 GEOSCIENCES↗