Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “statistical models”

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,045 records · Page 58

Comparison of a 2D Photochemical Model to Data Using Statistical Trend Analysis

We have analyzed our 23-year merged ozone data set for variability and trends with a statistical time-series model. To assist in that analysis, we have analyzed the Goddard ED photochemical model for the same time period with the same time-series model. Multiple runs of the photochemical model allow us to separate the effects of various terms on ozone, such as solar cycle and volcanic eruptions. We use this to help us separate those signals from each other in the data. We also used a smoothed version of the photochemical model's prediction of global ozone change in place of a simple linear trend. We find a number of interesting results. This particular photochemical model is more sensitive to chlorine perturbations than the atmosphere appears to be. It is less sensitive to solar cycle. It predicts an effect from the Pinatubo eruption that is nearly symmetric in the two hemispheres, but the data appears to have not responded to Pinatubo in the southern mid-latitudes. These results and their uncertainties will be discussed.

Stolarski, Richard S.↗

Comparison of a 2D Photochemical Model to Data Using Statistical Trend Analysis

We have analyzed our 23-year merged ozone data set for variability and trends with a statistical time-series model. To assist in that analysis, we have analyzed the Goddard 2D photochemical model for the same time period with the same time-series model. Multiple runs of the photochemical model allow us to separate the effects of various terms on ozone, such as solar cycle and volcanic eruptions. We use this to help us separate those signals from each other in the data. We also used a smoothed version of the photochemical model's prediction of global ozone change in place of a simple linear trend. We find a number of interesting results. This particular photochemical model is more sensitive to chlorine perturbations than the atmosphere appears to be. It is less sensitive to solar cycle. It predicts an effect from the Pinatubo eruption that is nearly symmetric in the two hemispheres, but the data appears to have not responded to Pinatubo in the southern mid-latitudes. These results and their uncertainties will be discussed.

Bhartia, P. K.↗

Statistics of internal stress fluctuations in dislocated crystals and relevance to density-based dislocation dynamics models

A statistical analysis of internal stress fluctuations, defined as the difference between the local mean stress and stress on dislocations, is presented for deforming crystals with 3D discrete dislocation systems. Dislocation realizations are generated using dislocation dynamics simulations and the associated stress field is computed as a superposition of a regularized stress field of dislocation lines within the domain of the solution and a complementary stress field computed via a finite-element boundary value problem. The internal stress fluctuations of interest are defined by an ensemble of the difference between the stress on dislocation lines and the local mean field stress in the crystal. The latter is established in a piecewise fashion over small voxels in the crystal thus allowing the difference between the local average stress and stress on segments to be easily estimated. The results show that the Schmid stress (resolved shear stress) and Escaig stress fluctuations on various slip systems sampled over a random set of points follow a Cauchy (Lorentz) distribution at all strain levels, with the amplitude and width of the distribution being dependent on the strain. Finally, the implications of the Schmid and Escaig internal stress fluctuations are discussed from the points of view of dislocation cross-slip and the dislocation motion in continuum dislocation dynamics.

36 MATERIALS SCIENCE↗

How to optimally constrain galaxy assembly bias: supplement projected correlation functions with count-in-cells statistics

Most models for the statistical connection between galaxies and their haloes ignore the possibility that galaxy properties may be correlated with halo properties other than halo mass, a phenomenon known as galaxy assembly bias. And yet, it is known that such correlations can lead to systematic errors in the interpretation of survey data that are analysed using traditional halo occupation models. At present, the degree to which galaxy assembly bias may be present in the real Universe, and the best strategies for constraining it remain uncertain. We study the ability of several observables to constrain galaxy assembly bias from redshift survey data using the decorated halo occupation distribution (dHOD), an empirical model of the galaxy-halo connection that incorporates assembly bias. We cover an expansive set of observables, including the projected two-point correlation function w p ( r p ), the galaxy-galaxy lensing signal ΔΣ( r p ), the void probability function VPF( r ), the distributions of counts-in-cylinders P(NCIC), and counts-in-annuli P(NCIA), and the distribution of the ratio of counts in cylinders of different sizes P ( N 2 /N 5 ). We find that despite the frequent use of the combination w p ( r p ) + ΔΣ( r p ) in interpreting galaxy data, the count statistics, P ( N CIC ) and P ( N CIA ), are generally more efficient in constraining galaxy assembly bias when combined with w p ( r p ). Constraints based upon w p ( r p ) and ΔΣ(rp) share common degeneracy directions in the parameter space, while combinations of r p w p ( r p ) with the count statistics are more complementary. Therefore, we strongly suggest that count statistics should be used to complement the canonical observables in future studies of the galaxy-halo connection.

galaxies: evolution↗

Modeling regional precipitation over the Indus River basin of Pakistan using statistical downscaling

Complex processes govern spatiotemporal distribution of precipitation within the high-mountainous headwater regions (commonly known as the upper Indus basin (UIB)), of the Indus River basin of Pakistan. Reliable precipitation simulations particularly over the UIB present a major scientific challenge due to regional complexity and inadequate observational coverage. Here, we present a statistical downscaling approach to model observed precipitation of the entire Indus basin, with a focus on UIB within available data constraints. Taking advantage of recent high altitude (HA) observatories, we perform precipitation regionalization using K-means cluster analysis to demonstrate effectiveness of low-altitude stations to provide useful precipitation inferences over more uncertain and hydrologically important HA of the UIB. We further employ generalized linear models (GLM) with gamma and Tweedie distributions to identify major dynamic and thermodynamic drivers from a reanalysis dataset within a robust cross-validation framework that explain observed spatiotemporal precipitation patterns across the Indus basin. Final statistical models demonstrate higher predictability to resolve precipitation variability over wetter southern Himalayans and different lower Indus regions, by mainly using different dynamic predictors. The modeling framework also shows an adequate performance over more complex and uncertain trans-Himalayans and the northwestern regions of the UIB, particularly during the seasons dominated by the westerly circulations. However, the cryosphere-dominated trans-Himalayan regions, which largely govern the basin hydrology, require relatively complex models that contain dynamic and thermodynamic circulations. Furthermore, we also analyzed relevant atmospheric circulations during precipitation anomalies over the UIB, to evaluate physical consistency of the statistical models, as an additional measure of reliability. Overall, our results suggest that such circulation-based statistical downscaling has the potential to improve our understanding towards distinct features of the regional-scale precipitation across the upper and lower Indus basin. Additionally, such understanding should help to assess the response of this complex, data-scarce, and climate-sensitive river basin amid future climatic changes, to serve communal and scientific interests.

54 ENVIRONMENTAL SCIENCES↗

Structural model optimization using statistical evaluation

The results of research in applying statistical methods to the problem of structural dynamic system identification are presented. The study is in three parts: a review of previous approaches by other researchers, a development of various linear estimators which might find application, and the design and development of a computer program which uses a Bayesian estimator. The method is tried on two models and is successful where the predicted stiffness matrix is a proper model, e.g., a bending beam is represented by a bending model. Difficulties are encountered when the model concept varies. There is also evidence that nonlinearity must be handled properly to speed the convergence.

Collins, J. D.↗

Modeling Soot Oxidation and Gasification with Bayesian Statistics

This paper presents a statistical method for model calibration using data collected from literature. The method is used to calibrate parameters for global models of soot consumption in combustion systems. This consumption is broken into two different submodels: first for oxidation where soot particles are attacked by certain oxidizing agents; second for gasification where soot particles are attacked by H 2 O or CO 2 molecules. Rate data were collected from 19 studies in the literature and evaluated using Bayesian statistics to calibrate the model parameters. Bayesian statistics are valued in their ability to quantify uncertainty in modeling. The calibrated consumption model with quantified uncertainty is presented here along with a discussion of associated implications. The oxidation results are found to be consistent with previous studies. Significant variation is found in the CO 2 gasification rates.

01 COAL, LIGNITE, AND PEAT↗

Graph link prediction in computer networks using Poisson matrix factorisation

Graph link prediction is an important task in cybersecurity: relationships between entities within a computer network, such as users interacting with computers or system libraries and the corresponding processes that use them, can provide key insights into adversary behaviour. Poisson matrix factorisation (PMF) is a popular model for link prediction in large networks, particularly useful for its scalability. In this article PMF is extended to include scenarios that are commonly encountered in cybersecurity applications. Specifically, an extension is proposed to explicitly handle binary adjacency matrices and include known categorical covariates associated with the graph nodes. A seasonal PMF model is also presented to handle seasonal networks. To allow the methods to scale to large graphs, variational methods are discussed for performing fast inference. The results show an improved performance over the standard PMF model and other statistical network models.

97 MATHEMATICS AND COMPUTING↗

Statistical and Prediction modeling of the Ka Band Using Experimental Results from ACTS Propagation Terminals at 20.185 and 27.505 GHZ

With the increase in demand for wireless communication services, most of the operating frequency bands have become very congested. The increase of wireless costumers is only fractional contribution to this phenomenon. The demand for more services such as video streams and internet explorer which require a lot of band width has been a more significant contributor to the congestion in a communication system. One way to increase the amount of information or data per unit of time transmitted with in a wireless communication system is to use a higher radio frequency. However in spite the advantage available in the using higher frequency bands such as, the Ka-band, higher frequencies also implies short wavelengths. And shorter wavelengths are more susceptible to rain attenuation. Until the Advanced Communication Technology Satellite (ACTS) was launched, the Ka- band frequency was virtually unused - the majority of communication satellites operated in lower frequency bands called the C- and Ku- bands. Ka-band is desirable because its higher frequency allows wide bandwidth applications, smaller spacecraft and ground terminal components, and stronger signal strength. Since the Ka-band is a high frequency band, the millimeter wavelengths of the signals are easily degraded by rain. This problem known as rain fade or rain attenuation The Advanced Communication Technology Satellite (ACTS) propagation experiment has collected 5 years of Radio Frequency (RF) attenuation data from December 1993 to November 1997. The objective of my summer work is to help develop the statistics and prediction techniques that will help to better characterize the Ka Frequency band. The statistical analysis consists of seasonal and cumulative five-year attenuation statistics for the 20.2 and 27.5 GHz. The cumulative five-year results give the link outage that occurs for a given link margin. The experiment has seven ground station terminals that can be attributed to a unique rain zone climate. The locations are White Sands, NM, Tampa, Fly Clarksburg, MD, Norman, OK, Ft. Collins, COY Vancouver, BC, and Fairbanks, AK. The analysis will help us to develop and define specific parameters that will help system engineers develop the appropriated instrumentation and structure for a Ka-band wireless communication systems and networks.

Ogunwuyi, Oluwatosin O.↗

A statistical rain attenuation prediction model with application to the advanced communication technology satellite project. 1: Theoretical development and application to yearly predictions for selected cities in the United States

A rain attenuation prediction model is described for use in calculating satellite communication link availability for any specific location in the world that is characterized by an extended record of rainfall. Such a formalism is necessary for the accurate assessment of such availability predictions in the case of the small user-terminal concept of the Advanced Communication Technology Satellite (ACTS) Project. The model employs the theory of extreme value statistics to generate the necessary statistical rainrate parameters from rain data in the form compiled by the National Weather Service. These location dependent rain statistics are then applied to a rain attenuation model to obtain a yearly prediction of the occurrence of attenuation on any satellite link at that location. The predictions of this model are compared to those of the Crane Two-Component Rain Model and some empirical data and found to be very good. The model is then used to calculate rain attenuation statistics at 59 locations in the United States (including Alaska and Hawaii) for the 20 GHz downlinks and 30 GHz uplinks of the proposed ACTS system. The flexibility of this modeling formalism is such that it allows a complete and unified treatment of the temporal aspects of rain attenuation that leads to the design of an optimum stochastic power control algorithm, the purpose of which is to efficiently counter such rain fades on a satellite link.

Manning, Robert M.↗

Data-driven prediction of scaling and ignition of inertial confinement fusion experiments

Recent advances in inertial confinement fusion (ICF) at the National Ignition Facility (NIF), including ignition and energy gain, are enabled by a close coupling between experiments and high-fidelity simulations. Neither simulations nor experiments can fully constrain the behavior of ICF implosions on their own, meaning pre- and postshot simulation studies must incorporate experimental data to be reliable. Linking past data with simulations to make predictions for upcoming designs and quantifying the uncertainty in those predictions has been an ongoing challenge in ICF research. We have developed a data-driven approach to prediction and uncertainty quantification that combines large ensembles of simulations with Bayesian inference and deep learning. The approach builds a predictive model for the statistical distribution of key performance parameters, which is jointly informed by past experiments and physics simulations. The prediction distribution captures the impact of experimental uncertainty, expert priors, design changes, and shot-to-shot variations. We have used this new capability to predict a 10× increase in ignition probability between Hybrid-E shots driven with 2.05 MJ compared to 1.9 MJ, and validated our predictions against subsequent experiments. We describe our new Bayesian postshot and prediction capabilities, discuss their application to NIF ignition and validate the results, and finally investigate the impact of data sparsity on our prediction results.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Southern Ocean Uptake in the MPAS-Ocean Model

The research added to the understanding of Southern Ocean Uptake in the following ways. 1) No single process in the model could be “tuned” to improve Southern Ocean uptake of heat, CO 2 and CFCs. 2) More comprehensive studies revealed that multiple aspects of the model are important; in particular the general circulation control of deep (>200m) stratification (Small et al., 2019) and salinity contributions to the shallow (< 200 m) stratification (DuVivier et al, 2018). These two aspects are not strongly connected. The first can be very much improved by increasing horizontal resolution, primarily because of more transport of high salinity water from the Indian Ocean into the Southern Ocean south of Africa, so we now know that “parameterizing” this salinity source will be one necessary step to improving the Southern Ocean uptake at low resolution. 3) However, the second was found to be more than just a resolution issue, which may explain why many studies have found that higher resolution improves only some aspects of model solutions, but degrades others. Even in cases where the surface wind is too strong, the shallow stratification is sufficient to inhibit vertical mixing even though the salinity stabilization is weak compared to observations. The conclusion is that the process responsible for mixing through the stable salinity barrier above 200m is missing from the model. Therefore, a second necessary step is to find the process then incorporate it into the model. The project focused on ocean surface wave driven Langmuir turbulence as the source of additional mixing by studying Large Eddy Simulations (LES) of Southern Ocean boundary layers as they eroded a realistic salinity barrier. The wind, buoyancy and wave driven boundary layer turbulence was found to be well described by Monin-Obukhov similarity by developing similarity functions of the wave driving, and combining these with the established wind and buoyancy functions (Large et al., 2018a). These functions and a novel analysis of the momentum flux and vertical shear fields from the LES revealed two distinct aspects where the boundary layer turbulence was much more organized than expected, and hence much more amenable to implementation in the model (Large et al., 2018b). First, there is a robust non-local across-shear transport of momentum throughout the ocean boundary layer; a hitherto missing process. Furthermore, this transport is related to the across-shear component of the wind stress in the same way that non-local buoyancy transport is related to the surface buoyancy flux, both in the ocean and atmosphere. Even the coefficients are similar in the convective limit. Second, the turbulent fluxes of momentum and buoyancy are distributed through the ocean boundary layer of the LES in a very orderly fashion given by a single shape function that appears to be applicable to both momentum and buoyancy, with or without surface wave driving. Thus, a very technically effective and economic path to incorporating surface wave driving in the vertical mixing scheme of the model is now clear. The high degree of order found for ocean boundary layer turbulence, including Langmuir turbulence, is more amenable to empirical models than to statistical turbulence models, and the computational cost is very much less. Fortuitously, the empirical K-Profile Parameterization is already implemented in the model and the wave-driving similarity function, the non-local momentum transport and the vertical shape function are straightforward enhancements.

58 GEOSCIENCES↗

High-Resolution Regional Wave Hindcast for the U.S. Alaska Coast

Regional long-term wave hindcast with state-of-art third-generation models is essential for characterization of regional wave resources. This report summarizes modeling efforts for the simulation of the wave climate along southern Alaska, which used an unstructured, nested-grid modeling approach that incorporates a global-regional nested grid using WAVEWATCH III ® and the high-resolution SWAN (Simulating Waves Nearshore) model. Wave resource and wave bulk parameters were simulated for a 32-year period from 1979 to 2010 and were subsequently validated with wave buoy data in the model domain. Error statistics of model skills were calculated at all buoy stations. Overall, model results, both resource and bulk parameters, match observations well. The model hindcast was able to reproduce the seasonal variation of the sea state with large waves that occur in the winter and early spring months when wind forcing is strong and the calm sea state during the summer when wind forcing is weak. The hindcast also considers the effect of sea ice on wave growth, which affects the southern Bering Sea. The hindcast results, including resource and bulk parameters, were archived at every model grid point for the entire model domain at three-hour intervals. In addition, spectral results were saved at National Data Buoy Center buoy stations and selected virtual buoy locations at hourly intervals. The nested-grid modeling framework employed in this study provides a powerful and efficient modeling approach for accurately simulating wave climate at regional and long-term temporal scales with sufficiently fine resolutions in the nearshore region.

54 ENVIRONMENTAL SCIENCES↗

Properties of hadron matter. II - Dense baryon matter and neutron stars.

In this article we have provided certain details of a nuclear-matter computation, based on the Brueckner-Bethe-Goldstone theory of nuclear reaction, which leads to an equation of state for matter in the density region of 10 to 500 trillion g/cu cm. We also explore the possibilities that at very high baryon densities or for very short baryon separations, the net baryon-baryon interaction may be negligible so that the results of dynamical models, like the statistical bootstrap model and the dual-resonance model, may be applicable to the study of dense baryon matter. Several plausible equations of state are constructed, and their effect on the limiting mass of the neutron star is examined.

Leung, Y. C.↗

Model optimization using statistical estimation

Program revises initial or prior estimate of stiffness and mass parameters to parameters yielding frequency and mode characteristics in agreement with test data. Variances are also calculated and consequently define uncertainties of final estimates.

Collins, J. D.↗