Search NASA⌕ Search

SEARCH · Search NASA

Results for “Statistical Physics”

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 775 records · Page 43

Modeling Flow and Transport in Fractured Rock using Machine Learning [Slides]

Fractured systems in the subsurface play a role in many natural and engineered applications such as geologic carbon sequestration, hydraulic fracturing and underground nuclear test detection. Structural information (fracture size, orientation, etc.) plays a key role in governing the dominant physics for these systems but can only be known statistically. Traditional modeling approaches either ignore or idealize structural information at these larger scales because we lack a computational framework that utilizes it in its entirety. The work presented here integrates computational physics, machine learning and graph theory to make a paradigm shift from computationally intensive high-fidelity models to coarse-scale graphs without loss of critical structural information. We exploit the underlying discrete structure of fracture networks in systems considering flow through fractures and fracture propagation. We demonstrate that compact graph representations require significantly fewer degrees of freedom (dof) to capture micro-fracture information and further accelerate these models with Machine Learning.

58 GEOSCIENCES↗

Physics through the 1990s: Elementary-particle physics

The volume begins with a non-mathematical discussion of the motivation behind, and basic ideas of, elementary-particle physics theory and experiment. The progress over the past two decades with the quark model and unification of the electromagnetic and weak interactions is reviewed. Existing theoretical problems in the field, such as the origin of mass and the unification of the fundamental forces, are detailed, along with experimental programs to test the new theories. Accelerators, instrumentation, and detectors are described for both current and future facilities. Interactions with other areas of both theoretical and applied physics are presented. The sociology of the field is examined regarding the education of graduate students, the organization necessary in large-scale experiments, and the decision-making process involved in high-cost experiments. Finally, conclusions and recommendations for maintaining US excellence in theory and experiment are given. Appendices list both current and planned accelerators, and present statistical data on the US elementary-particle physics program. A glossary is included.

Source record↗

Effects of roughness on density-weighted particle statistics in turbulent channel flows

The distribution of inertial particles in turbulent flows is strongly influenced by the characteristics of the coherent turbulent structures which develop in the carrier flow field. In wall-bounded flows, these turbulent structures, which control the turbulent regeneration cycles, are strongly affected by the roughness of the wall, nevertheless its effects on the particle transport in two-phase turbulent flows has been still poorly investigated. The issue is discussed here by addressing DNS combined with LPT to obtain statistics of velocity and preferential accumulation of a dilute dispersion of heavy particles in a turbulent channel flow, bounded by irregular two-dimensional rough surfaces, in the one-way coupling regime.

42 ENGINEERING↗

The Sloan Digital Sky Survey Reverberation Mapping Project: An investigation of biases in C iv emission line properties

We investigate the dependence on data quality of quasar properties measured from the C iv emission line region at high redshifts. Our measurements come from 32 epochs of Sloan Digital Sky Survey Reverberation Mapping Project spectroscopic observations of 482 z>1.46 quasars. We compare the differences between measurements made from the single-epoch (SE) and coadded spectra, focusing on the C iv λ1549 emission line because of its importance for studies of high-redshift quasar demographics and physical properties, including black hole masses. In addition to statistical errors increasing (by factors of ∼2–4), we find increasing systematic offsets with decreasing signal-to-noise ratio (S/N). The systematic difference (measurement uncertainty) in our lowest-S/N (<5) subsample between the SE and coadded spectrum (i) C iv equivalent width is 17 Å (31 Å), (ii) centroid wavelength is <1 Å (2 Å), and fractional velocity widths, ΔV/V, characterized by (iii) the line dispersion, σ {sub l}, is 0.104 (0.12), and (iv) the mean absolute deviation (MAD) is 0.072 (0.11). These remain smaller than the 1σ measurement uncertainties for all subsamples considered. The MAD is found to be the most robust line-width characterization. Offsets in the C iv FWHM velocity width and the C iv profile characterized by FWHM/σ {sub l} are only smaller than the statistical uncertainties when S/N > 10, although offsets in lower-S/N spectra exceed the statistical uncertainties by only a factor of ∼1.5 and may depend on the type of functional fit to the line. Characterizing the C iv line profile by the kurtosis is the least robust property investigated, as the median systematic coadded–SE measurement differences are larger than the statistical uncertainties for all S/N subsamples.

79 ASTRONOMY AND ASTROPHYSICS↗

Likelihood-Based Climate Model Evaluation

Climate models are deterministic, mathematical descriptions of the physics of climate. Confidence in predictions of future climate is increased if the physics are verifiably correct. A necessary, (but not sufficient) condition is that past and present climate be simulated well. Quantify the likelihood that a (summary statistic computed from a) set of observations arises from a physical system with the characteristics captured by a model generated time series. Given a prior on models, we can go further: posterior distribution of model given observations.

climate models↗

Applications of the BIOPHYS Algorithm for Physically-Based Retrieval of Biophysical, Structural and Forest Disturbance Information

Canopy reflectance model inversion using look-up table approaches provides powerful and flexible options for deriving improved forest biophysical structural information (BSI) compared with traditional statistical empirical methods. The BIOPHYS algorithm is an improved, physically-based inversion approach for deriving BSI for independent use and validation and for monitoring, inventory and quantifying forest disturbance as well as input to ecosystem, climate and carbon models. Based on the multiple-forward mode (MFM) inversion approach, BIOPHYS results were summarized from different studies (Minnesota/NASA COVER; Virginia/LEDAPS; Saskatchewan/BOREAS), sensors (airborne MMR; Landsat; MODIS) and models (GeoSail; GOMS). Applications output included forest density, height, crown dimension, branch and green leaf area, canopy cover, disturbance estimates based on multi-temporal chronosequences, and structural change following recovery from forest fires over the last century. Good correspondences with validation field data were obtained. Integrated analyses of multiple solar and view angle imagery further improved retrievals compared with single pass data. Quantifying ecosystem dynamics such as the area and percent of forest disturbance, early regrowth and succession provide essential inputs to process-driven models of carbon flux. BIOPHYS is well suited for large-area, multi-temporal applications involving multiple image sets and mosaics for assessing vegetation disturbance and quantifying biophysical structural dynamics and change. It is also suitable for integration with forest inventory, monitoring, updating, and other programs.

Peddle, Derek R.↗

The joint space-time statistics of macroweather precipitation, space-time statistical factorization and macroweather models

Over the range of time scales from about 10 days to 30–100 years, in addition to the familiar weather and climate regimes, there is an intermediate “macroweather” regime characterized by negative temporal fluctuation exponents: implying that fluctuations tend to cancel each other out so that averages tend to converge. We show theoretically and numerically that macroweather precipitation can be modeled by a stochastic weather-climate model (the Climate Extended Fractionally Integrated Flux, model, CEFIF) first proposed for macroweather temperatures and we show numerically that a four parameter space-time CEFIF model can approximately reproduce eight or so empirical space-time exponents. In spite of this success, CEFIF is theoretically and numerically difficult to manage. We therefore propose a simplified stochastic model in which the temporal behavior is modeled as a fractional Gaussian noise but the spatial behaviour as a multifractal (climate) cascade: a spatial extension of the recently introduced ScaLIng Macroweather Model, SLIMM. Both the CEFIF and this spatial SLIMM model have a property often implicitly assumed by climatologists that climate statistics can be “homogenized” by normalizing them with the standard deviation of the anomalies. Physically, it means that the spatial macroweather variability corresponds to different climate zones that multiplicatively modulate the local, temporal statistics. This simplified macroweather model provides a framework for macroweather forecasting that exploits the system's long range memory and spatial correlations; for it, the forecasting problem has been solved. We test this factorization property and the model with the help of three centennial, global scale precipitation products that we analyze jointly in space and in time.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Improving the Fidelity of Lattice QCD Calculations

I begin with an overview of the various statistical and systematic uncertainties that enter into lattice calculations in nuclear physics, and the hierarchy of computations needed to investigate them. I emphasize in particular the need to obtain as precise a statistical error as possible on a given ensemble of calculations, and the various variance reduction and algorithmic methods that have been proposed recently with the aim of minimizing the statistical uncertainties. I conclude with a description of the fitting methods used in lattice QCD calculations.

Richards, David↗

Engineering Super–Poissonian Photon Statistics of Spatial Light Modes

The nature of light sources is defined by the statistical fluctuations of the electromagnetic field. As such, the photon statistics of light sources are typically associated with distinct emitters. Here, the possibility of producing light beams with various photon statistics through the spatial modulation of coherent light is demonstrated. This is achieved by the sequential encoding of controllable Kolmogorov phase screens in a digital micromirror device. Interestingly, the flexibility of this scheme allows for the shaping of spatial light modes with engineered photon statistics at different spatial positions. The performance of this scheme is assessed through the photon-number-resolving characterization of different families of spatial light modes with engineered photon statistics. Furthermore, it is believed that the possibility of controlling the photon fluctuations of the light field at arbitrary spatial locations has important implications for quantum spectroscopy, sensing, and imaging.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Recurrence theorems: A unified account

I discuss classical and quantum recurrence theorems in a unified manner, treating both as generalisations of the fact that a system with a finite state space only has so many places to go. Along the way, I prove versions of the recurrence theorem applicable to dynamics on linear and metric spaces and make some comments about applications of the classical recurrence theorem in the foundations of statistical mechanics.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Parameters Impacting Columnated Granular Soil Pneumatic Seal Performance

The ability of a column of loose granular soil to form a pneumatic seal was investigated by varying the diameter of the soil column, the effective column height, and the level of compaction in the soil. The soil column diameter was tested at three levels using pipes with inner diameters measuring 5.08, 10.16, and 15.24 centimeters (2, 4, and 6 inches). Soil was filled in each pipe to form 15, 30, and 45 centimeter (6, 12, and 18 inch) tall soil columns. Each diameter/height configuration was also tested at three levels of soil compaction, compared by calculating the bulk density of the soil with mass and volume measurements. GRC-1a simulant was used, with approximate low/medium/high bulk densities of 1.6, 1.75, and 1.9 g/cc achieved with a combination of vibration and tamping. The order of tests for a given column diameter was randomized and repeated three times. With the top of the soil open to atmosphere at room temperature, compressed air was injected through a small diffuser at the column base with several small downward-facing holes. The number and size of these holes was scaled such that a constant total inlet orifice area to column cross sectional area ratio was maintained for each column diameter. Inlet air pressure was slowly increased via a precision regulator to preserve quasi-static equilibrium in the soil column to minimize the impact of dynamics. Air pressure was increased all the way through the static and bubbling regimes until slugging or turbulent behavior was observed in the soil to ensure that the entire static regime had been captured during data collection. A three-factor, three-level analysis of variance (ANOVA) statistical analysis was performed on the resulting data to determine the extent to which each physical parameter impacted the soil column seal performance. It was concluded that there is statistically significant evidence that column height, and the interaction between column height and diameter impact soil seal performance. In all other cases there was insufficient data to identify a statistically significant causal relationship. Additionally, plots were generated comparing experimental data to the predictive formula developed by Ogino et al. for fluidized beds in 1993. Because this model was developed for industrial spouted fluidized beds, the accuracy of its output prior to fluidization in the static seal ‘edge case’ is unknown, especially considering in this application the working gas was diffused across the column base rather than being injected through a spout. Further, the fidelity of the model had not yet been tested with lunar soil simulants. Plotting the Ogino et al. model alongside test data allows for a more intuitive sense of the impact of test parameters on soil seal performance, as well as providing a quick means to further tune this predictive model for more accurate use with static seals across granular lunar soil simulants. The model provided by Ogino et al. was further tuned using test data to determine the degree to which a spouted bed model could be applied to a slightly modified set of testing conditions: lunar soil simulants and a more diffuse, homogeneous application of pneumatic pressure. A Matlab script was created to test different values for the leading coefficient and exponents in Ogino’s formula. The script swept preset ranges, then iterated with higher resolutions across narrower ranges to converge on the optimal value for each parameter. The resulting adjusted model was compared to the original, as well as test data with noticeable improvements across the entire test domain.

Jack Stewart↗

Bivariate Kumaraswamy distribution with an application on earthquake data

Bivariate Kumaraswamy (BK) distribution whose marginals are Kumaraswamy distributions has been recently introduced. However, its statistical properties are not studied in detail. In this study, statistical properties of the BK distribution are investigated. We suggest that the BK could provide suitable description for the earthquakes characteristics of Turkey. We support this argument using earthquakesoccurred in Turkey between 1900 and 2009. We also find that the BK distribution simulates earthquakes well.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

On thermalization of electron-positron-photon plasma

Recently a progress has been made in understanding thermalization mechanism of relativistic plasma starting from a non-equilibrium state. Relativistic Boltzmann equations were solved numerically for homogeneous isotropic plasma with collision integrals for two- and three-particle interactions calculated from the first principles by means of QED matrix elements. All particles were assumed to fulfill Boltzmann statistics. In this work we follow plasma thermalization by accounting for Bose enhancement and Pauli blocking in particle interactions. Our results show that particle in equilibrium reach Bose-Einstein distribution for photons, and Fermi-Dirac one for electrons, respectively.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Is solar neutrino capture rate correlated with sunspot number?

The statistical significance of the apparent correlation between sunspots and the observed neutrino rate is quantified. It is shown that the correlation depends almost entirely upon four low neutrino capture rates near the beginning of 1980. A calculation based on standard electroweak theory and neutrino production processes demonstrates that a correlation, if real, would be extremely puzzling on energetic grounds alone. It is concluded that measurements with the Cl-37 detector during the next sunspot cycle will be needed to show that there is a physical correlation, since the existing data are not statistically significant at a definitive level.

Bahcall, J. N.↗

Modeling single-molecule stretching experiments using statistical thermodynamics

Single-molecule stretching experiments are widely utilized within the fields of physics and chemistry to characterize the mechanics of individual bonds or molecules, as well as chemical reactions. Analytic relations describing these experiments are valuable, and these relations can be obtained through the statistical thermodynamics of idealized model systems representing the experiments. Since the specific thermodynamic ensembles manifested by the experiments affect the outcome, primarily for small molecules, the stretching device must be included in the idealized model system. Though the model for the stretched molecule might be exactly solvable, including the device in the model often prevents analytic solutions. In the limit of large or small device stiffness, the isometric or isotensional ensembles can provide effective approximations, but the device effects are missing. Here a dual set of asymptotically correct statistical thermodynamic theories are applied to develop accurate approximations for the full model system that includes both the molecule and the device. In conclusion, the asymptotic theories are first demonstrated to be accurate using the freely jointed chain model and then using molecular dynamics calculations of a single polyethylene chain.

74 ATOMIC AND MOLECULAR PHYSICS↗

To have value, comparisons of high-throughput phenotyping methods need statistical tests of bias and variance

The gap between genomics and phenomics is narrowing. The rate at which it is narrowing, however, is being slowed by improper statistical comparison of methods. Quantification using Pearson’s correlation coefficient ( r ) is commonly used to assess method quality, but it is an often misleading statistic for this purpose as it is unable to provide information about the relative quality of two methods. Using r can both erroneously discount methods that are inherently more precise and validate methods that are less accurate. These errors occur because of logical flaws inherent in the use of r when comparing methods, not as a problem of limited sample size or the unavoidable possibility of a type I error. A popular alternative to using r is to measure the limits of agreement (LOA). However both r and LOA fail to identify which instrument is more or less variable than the other and can lead to incorrect conclusions about method quality. An alternative approach, comparing variances of methods, requires repeated measurements of the same subject, but avoids incorrect conclusions. Variance comparison is arguably the most important component of method validation and, thus, when repeated measurements are possible, variance comparison provides considerable value to these studies. Statistical tests to compare variances presented here are well established, easy to interpret and ubiquitously available. The widespread use of r has potentially led to numerous incorrect conclusions about method quality, hampering development, and the approach described here would be useful to advance high throughput phenotyping methods but can also extend into any branch of science. The adoption of the statistical techniques outlined in this paper will help speed the adoption of new high throughput phenotyping techniques by indicating when one should reject a new method, outright replace an old method or conditionally use a new method.

59 BASIC BIOLOGICAL SCIENCES↗

Proton imaging of stochastic magnetic fields

Recent laser-plasma experiments report the existence of dynamically significant magnetic fields, whose statistical characterization is essential for a complete understanding of the physical processes these experiments are attempting to investigate. In this paper, we show how a proton-imaging diagnostic can be used to determine a range of relevant magnetic-field statistics, including the magnetic-energy spectrum. To achieve this goal, we explore the properties of an analytic relation between a stochastic magnetic field and the image-flux distribution created upon imaging that field. This ‘Kugland image-flux relation’ was previously derived under simplifying assumptions typically valid in actual proton-imaging set-ups. We conclude that, as with regular electromagnetic fields, features of the beam’s final image-flux distribution often display a universal character determined by a single, field-scale dependent parameter – the contrast parameter $\unicode[STIX]{x1D707}\equiv d_{s}/{\mathcal{M}}l_{B}$ – which quantifies the relative size of the correlation length $l_{B}$ of the stochastic field, proton displacements $d_{s}$ due to magnetic deflections and the image magnification ${\mathcal{M}}$ . For stochastic magnetic fields, we establish the existence of four contrast regimes, under which proton-flux images relate to their parent fields in a qualitatively distinct manner. These are linear, nonlinear injective, caustic and diffusive. The diffusive regime is newly identified and characterized. The nonlinear injective regime is distinguished from the caustic regime in manifesting nonlinear behaviour, but as in the linear regime, the path-integrated magnetic field experienced by the beam can be extracted uniquely. Thus, in the linear and nonlinear injective regimes we show that the magnetic-energy spectrum can be obtained under a further statistical assumption of isotropy. This is not the case in the caustic or diffusive regimes. We discuss complications to the contrast-regime characterization arising for inhomogeneous, multi-scale stochastic fields, which can encompass many contrast regimes, as well as limitations currently placed by experimental capabilities on one’s ability to extract magnetic-field statistics. Furthermore, the results presented in this paper are of consequence in providing a comprehensive description of proton images of stochastic magnetic fields, with applications for improved analysis of proton-flux images.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A variance deconvolution estimator for efficient uncertainty quantification in Monte Carlo radiation transport applications

Monte Carlo simulations are at the heart of many high-fidelity simulations and analyses for radiation transport systems. As is the case with any complex computational model, it is important to propagate sources of input uncertainty and characterize how they affect model output. Unfortunately, uncertainty quantification (UQ) is made difficult by the stochastic variability that Monte Carlo transport solvers introduce. The standard method to avoid corrupting the UQ statistics with the transport solver noise is to increase the number of particle histories, resulting in very high computational costs. In this contribution, we propose and analyze a sampling estimator based on the law of total variance to compute UQ variance even in the presence of residual noise from Monte Carlo transport calculations. Here, we rigorously derive the statistical properties of the new variance estimator, compare its performance to that of the standard method, and demonstrate its use on neutral particle transport model problems involving both attenuation and scattering physics. We illustrate, both analytically and numerically, the estimator’s statistical performance as a function of available computational budget and the distribution of that budget between UQ samples and particle histories. We show analytically and corroborate numerically that the new estimator is unbiased, unlike the standard approach, and is more accurate and precise than the standard estimator for the same computational budget.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗