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 1,549 records · Page 86
All sourcesNASA NTRSDOE OSTI
DOE OSTI2017

Sensitivity Analysis of Neutron Multiplicity Counting Statistics Using First-Order Perturbation Theory and Application to a Subcritical Plutonium Metal Benchmark

It is frequently important to estimate the uncertainty and sensitivity of measured and computed detector responses in subcritical experiments and simulations. These uncertainties arise from the physical construction of the experiment, uncertainties in the transport parameters, and counting uncertainties. Perturbation theory enables sensitivity analysis (SA) and uncertainty quantification on integral quantities like detector responses. Here the aim of our work is to apply SA to the statistics of subcritical neutron multiplicity counting distributions. Current SA methods have only been applied to mean detector responses and the keff eigenvalue. For multiplicity counting experiments, knowledge of the higher-order counting moments and their uncertainties are essential for a complete SA. We apply perturbation theory to compute the sensitivity of neutron multiplicity counting moments to arbitrarily high order. Each moment is determined by solving an adjoint transport equation with a source term that is a function of the adjoint solutions for lower-order moments. This enables moments of arbitrarily high order to be sequentially determined, and it shows that each moment is sensitive to the uncertainties of all lower-order moments. To close our SA of the moments, we derive forward transport equations that are functions of the forward flux and lower-order moment adjoint fluxes. We verify our calculations for the first three moments by comparison with multiplicity counting measurements of a subcritical plutonium metal sphere. For the first three moments, the most influential parameters are ranked, and the validity of first-order perturbation theory is demonstrated by examining the series truncation error. This enables a detailed SA of subcritical multiplicity counting measurements of fissionable material based on transport theory.

42 ENGINEERING↗
NASA NTRS1992

Statistics of the cosmic Mach number from numerical simulations of a cold dark matter universe

Results are presented of an analysis of the cosmic Mach number, M, the ratio of the streaming velocity, v, to the random velocity dispersion, sigma, of galaxies in a given patch of the universe, which was performed on the basis of hydrodynamical simulations of the cold dark matter scenario. Galaxy formation is modeled by application of detailed physical processes rather than by the ad hoc assumption of 'bias' between dark matter and galaxy fluctuations. The correlation between M and sigma is found to be very weak for both components. No evidence is found for a physical 'velocity bias' in the quantities which appear in the definition of M. Standard cold-dark-matter-dominated universes are in conflict, at a statistically significant level, with the available observation, in that they predict a Mach number considerably lower than is observed.

Suto, Yasushi↗
NASA NTRS1976

Statistical simulations of clusters of galaxies

By comparing observed and simulated rich clusters of galaxies, it is shown that the observed clusters actually possess physical cores. The accuracy with which core radii can be determined is found. It is also shown that the observations of density profiles of galaxies in the clusters give no significant evidence for a dynamical reason as the cause of the anomalously close resemblance found previously between such density profiles and the isothermal gas distribution.

Avni, Y.↗
NASA NTRS1984

Evaluation of the radar response to land surfaces and volumes: Examination of theoretical models, target statistics, and applications

Four areas of L-band radar remote sensing of terrain were examined: (1) the behavior of the radar backscatter coefficient of distributed surface and volumes as a function of the targets' dielectric and geometric parameters and as a fnction of their physical parameters; (2) the correspondence of the angular behavior of the relative backscatter coefficient as extracted from SIR-B digital imagery and truck mounted L-band scatterometer measurements for about 100 fields; (3) the statistical behavior of SIR-B image density for targets that appear homogeneous on Thematic Mapper (TM) optical imagery and/or color IR photography; and (4) the applicability of SIR-B imagery both alone and in conjunction with TM imagery for the classification and monitoring of land cover and renewable resources.

Ulaby, F. T.↗
NASA NTRS2000

Comments on "The Sensitivity Study of Radiative-Convective Equilibrium in the Tropics with a Convective Resolving Model"

In general, there are two broad scientific objectives when using cloud resolving models (CRMs or cloud ensemble models-CEMs) to study tropical convection. The first one is to use them as a physics resolving models to understand the dynamic and microphysical processes associated with the tropical water and energy cycles and their role in the climate system. The second approach is to use the CRMs to improve the representation of moist processes and their interaction with radiation in large-scale models. In order to improve the credibility of the CRMs and achieve the above goals, CRMs using identical initial conditions and large-scale influences need to produce very similar results. Two CRMs produced different statistical equilibrium (SE) states even though both used the same initial thermodynamic and wind conditions. Sensitivity tests to identify the major physical processes that determine the SE states for the different CRM simulations were performed. Their results indicated that atmospheric horizontal wind is treated quite differently in these two CRMs. The model that had stronger surface winds and consequently larger latent and sensible heat fluxes from the ocean produced a warmer and more humid modeled thermodynamic SE state. In addition, the domain mean thermodynamic state is more unstable for those experiments that produced a warmer and more humid SE state. Their simulated wet (warm and humid) SE states are thermally more stable in the lower troposphere (from the surface to 4-5 km in altitude). The large-scale horizontal advective effects on temperature and water vapor mixing ratio are needed when using CRMs to perform long-term integrations to study convective feedback under specified large-scale environments. In addition, it is suggested that the dry and cold SE state simulated was caused by enhanced precipitation but not enough surface evaporation. We find some problems with the interpretation of these three phenomena.

Tao, W.-K.↗
DOE OSTI2018

Statistical relationship between metabolic decomposition and chemical uptake predicts bioconcentration factor data for diverse chemical exposures

Background: A challenge of in vitro to in vivo extrapolation (IVIVE) is to predict the physical state of organisms exposed to chemicals in the environment from in vitro exposure assay data. Although toxicokinetic modeling approaches promise to bridge in vitro screening data with i Danio rerio ) based upon particle swarm optimizations (PSO) of the chemical uptake and degradation rates that predict bioconcentration factors (BCF) for a broad range of chemicals. PSO reveals a relationship between chemical uptake and decomposition parameter values that predicts chemical-specific BCF values with moderate statistical agreement to a limited yet diverse chemical dataset, and all without a need to retrain the model to new data. Conclusions: The presented model requires only the octanol-water partitioning ratio to predict BCFs to a fidelity consistent with existing QSAR models. This success begs re-evaluation of the modeling assumptions; specifically, it suggests that chemical uptake into arterial blood may be limited by transport across gill membranes (diffusion) rather than by counter-current flow between gill lamellae (convection). Therefore, more detailed molecular modeling of aquatic respiration may further improve predictive accuracy of the rTK approach.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Unified bulk semiclassical theory for intrinsic thermal transport and magnetization currents

In this work, we reveal the unexpected role of the material inhomogeneity in unifying the formulation of intrinsic thermal and thermoelectric transport as well as magnetization currents. The smooth inhomogeneity leads to the position dependent local band dispersion and phase-space Berry curvature, enabling a general and rapid access to transport and magnetization currents displaying the momentum-space Berry curvature physics. Our theory does not invoke the boundary current, the thermodynamic approach to magnetization or any mechanical counterpart of statistical forces. By introducing a fictitious inhomogeneity, it applies to homogeneous samples as well, promoting the inhomogeneity to be a basic trick in semiclassical transport theories. Such a trick works regardless of the driving force of transport, e.g., temperature gradient, in contrast to the trick of fictitious gravitational field in quantum transport theories. We thus include more general mechanical driving forces and establish the Mott relation between the resulting transport thermal and electric currents, whereas this relation for these two currents was previously only known when an electric field is the driving force.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Using data-science approaches to unravel insights for enhanced transport of lithium ions in single-ion conducting polymer electrolyte

Solid polymer electrolytes have yet to achieve the an ionic conductivity > 1 mS/cm at room temperature for realistic applications. This target implies the need to reduce the effective energy barriers of ion transport in polymer electrolytes to around 20 kJ/mol. In this work, we combine information extracted from existing experimental results with theoretical calculations to provide insights into ion transport in single-ion conductors (SICs) with a focus on lithium ion SICs. Through the analysis of temperature-dependent ionic conductivity data obtained from the literature, we evaluate different methods of extracting energy barriers for lithium transport. The traditional Arrhenius fit to the temperature-dependent ionic conductivity data indicates that the Meyer-Neldel rule holds for SICs. However, the values of the fitting parameters remain unphysical. Our modified approach based on recent work (Macromolecules, 56, 15, 6051(2023)), which incorporates a fixed pre-exponential factor, reveals that the energy barriers exhibit temperature dependence over a wide range of temperatures. Using this approach, we identify a series of anions leading to the energy barriers less than 30 kJ/mol, which include trifluoromethane sulfonimide (TFSI), fluoromethane sulfonimide (FSI), and boron-based organic anions. In our efforts to design the next generation of anions, which can exhibit the energy barriers less than 20 kJ/mol, we focused on boron-containing SICs, and performed density functional theory (DFT) based calculations to connect the chemical structures via the binding energy of cation (lithium)-anion pairs with the experimentally derived effective energy barriers for ion transport. Not only have we identified a correlation between the binding energy and the energy barriers, but we also propose a strategy to design new boron-based anions by using the correlation. This combined approach involving experiments and theoretical calculations is capable of facilitating the identification of promising new anions, which can exhibit ionic conductivity $> 1$ mS/cm near room temperature, thereby expediting the development of novel superionic single-ion conducting polymer electrolytes. The published datasets include all the temperature-dependent ionic conductivity collected from the literature with literature DOIs, DFT calculated binding energies, and python scripts to analyze data, construct statistical models, and generate plots.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Bayesian Estimation of Oscillator Parameters: Toward Anomaly Detection and Cyber-Physical System Security

Cyber-physical system security presents unique challenges to conventional measurement science and technology. Anomaly detection in software-assisted physical systems, such as those employed in additive manufacturing or in DNA synthesis, is often hampered by the limited available parameter space of the underlying mechanism that is transducing the anomaly. As a result, the formulation of anomaly detection for such systems often leads to inverse or ill-posed problems, requiring statistical treatments. Here, we present Bayesian inference of unknown parameters associated with a generic actuator considered as a representative vital element of a cyber-physical system. Via a series of experimental input-output measurements, a transfer function for the actuator is obtained numerically, which serves as our model for the proposed method. Linear, nonlinear, and delayed dynamics may be assumed for the actuator response. By devising a code-based malicious signal, we study the efficacy of Bayesian inference for its potential to produce a detection, including uncertainty quantification, with a remarkably small number of input data points. Our approach should be adaptable to a variety of real-time cyber-physical anomaly detection scenarios.

47 OTHER INSTRUMENTATION↗
DOE OSTI2018

Exact Maximum-Entropy Estimation with Feynman Diagrams

A longstanding open problem in statistics is finding an explicit expression for the probability measure which maximizes entropy with respect to given constraints. In this paper a solution to this problem is found, using perturbative Feynman calculus. The explicit expression is given as a sum over weighted trees.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Estimating Large Global Significances with a New Monte Carlo Extrapolation Method

In particle physics, it is needed to evaluate the possibility that excesses of events in mass spectra are due to statistical fluctuations as quantified by the standards of local and global significances. Without prior knowledge of a particle’s mass, it is especially critical to estimate its global significance. The usual approach is to count the number of times a significance limit is exceeded in a collection of simulated Monte Carlo (MC) “toy experiments”. To demonstrate this conventional method for global significance, we performed simulation studies according to a recent Compact Muon Solenoid (CMS) result to show its effectiveness. However, this counting method is not practical for computing large global significances. To address this problem, we developed a new “extrapolation” method to evaluate the global significance. We compared the global significance estimated by our new method with that of the conventional approach, and verified its feasibility and effectiveness. This method is also applicable for cases where only small toy MC samples are available. In this approach, the significance is calculated based on p-values, assuming symmetrical Gaussian distributions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2016

Updated constraints on velocity and momentum-dependent asymmetric dark matter

We present updated constraints on dark matter models with momentum-dependent or velocity-dependent interactions with nuclei, based on direct detection and solar physics. We improve our previous treatment of energy transport in the solar interior by dark matter scattering, leading to significant changes in fits to many observables. Based on solar physics alone, DM with a spin-independent q{sup 4} coupling provides the best fit to data, and a statistically satisfactory solution to the solar abundance problem. Once direct detection limits are accounted for however, the best solution is spin-dependent v{sup 2} scattering with a reference cross-section of 10{sup −35} cm{sup 2} (at a reference velocity of v{sub 0}=220 km s{sup −1}), and a dark matter mass of about 5 GeV.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Organic Electrochemical Transistor Channel Materials: Copolymerization Versus Physical Mixing of Glycolated and Alkoxylated Polymers

Organic electrochemical transistors (OECTs) feature a polymer channel capable of conducting both ions and electronic charges. The choice of the channel material is critical for OECT performance. Many efforts have focused on improving performance via the chemical tunability of conjugated polymers – through backbone, side chain, and molar mass engineering – leading to useful design principles for accumulation-mode OECT materials. However, tuning the chemical structure of conjugated polymers often requires time-consuming optimization of the synthesis route. Meanwhile, variations in molar mass, dispersity, structural defects, and metal content present challenges when attempting to analyze the detailed effects of structural modifications, as multiple performance-determining factors are often (unintentionally) changed at the same time. Therefore, this study explores blended channel materials obtained by physically mixing glycolated and alkoxylated polymers in different ratios, and compares their OECT performance with the corresponding statistical copolymers. It is shown that mixing two well-performing materials creates blends that enable rational tuning of the transistor properties without compromising on performance. Thus, channels based on blends of alkoxylated and glycolated polymers hold promise for OECT technology with tailored response, as only two materials are needed to achieve any desired side chain ratio, simplifying the optimization of OECT characteristics.

copolymerization↗
DOE OSTI2019

Modelling the Fluctuations of Reactive Shock Waves in Heterogeneous Solid Explosives as Stochastic Processes

Although current phenomenological burn models are useful for describing the average or bulk reactive flow behaviour of heterogeneous explosives, one fundamental weakness inherent to these models is the loss of detailed microstructural information at the scale of the calculation. In order to include the effects of the microstructure, and in particular the underlying material heterogeneities that influence the build-up to detonation, a new paradigm is put forth for modelling sub-grid, reaction-induced fluctuations (i.e. “hot spots”) at the continuum level. This modelling approach assumes that the reaction rate is stochastic, rather than deterministic, and it uses Langevin-type equations with a mathematical framework built upon Itô calculus and Lambourn's CIM model for shocked heterogeneous explosives. This approach follows directly from our previous letter, Ref. [1], and is inspired by the probability density function (pdf) methods used in turbulent reactive flows. Here, the stochastic burn model is derived, implemented, and exercised far beyond what has been shown in previous work. New hydrocode simulation results demonstrate the role of stochastic fluctuations during shock initiation; these fluctuations are approximated by collections of discrete particles, that evolve with drift (i.e. deterministic) and diffusion (i.e. stochastic) coefficients. Furthermore, the particle values are propagated and averaged to calculate the heat release, yield strength, and material impedance in each computational cell. Hydrocode simulation results further show how the fluctuating hot spot energy may or may not be transmitted to the wave front, and result in a detonation wave-like structure. The fundamental stochastic nature of this model permits simulations to have varying outcomes with the same initial conditions; this allows for go/no-go estimation (e.g., marginal or failed detonations), which might possibly be calibrated using the statistical distributions from real materials. Mesoscale calculations of shocked heterogeneous explosives also show that these fluctuations are physically justified (i.e., Ref. [2]), and it is hypothesized that the pdf functions provide a link between the meso(grain) and continuum scales for practical engineering calculations. Hence, our approach is a paradigm shift for building efficient, reduced order continuum burn models, that represent the sub-grid stochastic behaviour of shocked heterogeneous solid explosives.

42 ENGINEERING↗
DOE OSTI2019

Model reduction for fractured porous media: a machine learning approach for identifying main flow pathways

Discrete fracture networks (DFN) are often used to model flow and transport in fractured porous media. The accurate resolution of flow and transport behavior on a large DFN involving thousands of fractures is computationally expensive. This makes uncertainty quantification studies of quantities of interest such as travel time through the network computationally intractable, since hundreds to thousands of runs of the DFN model are required to get bounds on the uncertainty of the predictions. Prior works on the subject demonstrated that the complexity of a DFN could be reduced by considering a sub-network of it (often termed a “backbone” sub-network), one whose flow and transport properties were then shown to be similar to that of the full network. The technique is tantamount to partitioning the complete set of fractures of a network into two disjoint sets, one of which is the backbone sub-network while the other its complement. It is in this context that we present a system-reduction technique for DFNs using supervised machine learning via a Random Forest Classifier that selects a backbone sub-network from the full set of fractures. The in-sample errors (in terms of precision and recall scores) of the trained classifier are found to be very accurate indicators of the out-of-sample errors, thus exhibiting that the classifier generalizes well to test data. Moreover, this system-reduction technique yields sub-networks as small as 12% of the full DFN that still recover transport characteristics of the full network such as the peak dosage and tailing behavior for late times. Most importantly, the sub-networks remain connected, and their size can be controlled by a single dimensionless parameter. Furthermore, measures of KL-divergence and KS-statistic for the breakthrough curves of the sub-networks with respect to the full network show physically realistic trends in that the measures decrease monotonically as the size of the sub-networks increase. In conclusion, the computational efficiency gained by this technique depends on the size of the sub-network, but large reductions in computational time can be expected for small sub-networks, yielding as much as 90% computational savings for sub-networks that are as small as 10-12% of the full network.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2015

THE YOUNG SOLAR ANALOGS PROJECT. I. SPECTROSCOPIC AND PHOTOMETRIC METHODS AND MULTI-YEAR TIMESCALE SPECTROSCOPIC RESULTS

This is the first in a series of papers presenting methods and results from the Young Solar Analogs Project, which began in 2007. This project monitors both spectroscopically and photometrically a set of 31 young (300–1500 Myr) solar-type stars with the goal of gaining insight into the space environment of the Earth during the period when life first appeared. From our spectroscopic observations we derive the Mount Wilson S chromospheric activity index (S{sub MW}), and describe the method we use to transform our instrumental indices to S{sub MW} without the need for a color term. We introduce three photospheric indices based on strong absorption features in the blue-violet spectrum—the G-band, the Ca i resonance line, and the Hydrogen-γ line—with the expectation that these indices might prove to be useful in detecting variations in the surface temperatures of active solar-type stars. We also describe our photometric program, and in particular our “Superstar technique” for differential photometry which, instead of relying on a handful of comparison stars, uses the photon flux in the entire star field in the CCD image to derive the program star magnitude. This enables photometric errors on the order of 0.005–0.007 magnitude. We present time series plots of our spectroscopic data for all four indices, and carry out extensive statistical tests on those time series demonstrating the reality of variations on timescales of years in all four indices. We also statistically test for and discover correlations and anti-correlations between the four indices. We discuss the physical basis of those correlations. As it turns out, the “photospheric” indices appear to be most strongly affected by emission in the Paschen continuum. We thus anticipate that these indices may prove to be useful proxies for monitoring emission in the ultraviolet Balmer continuum. Future papers in this series will discuss variability of the program stars on medium (days–months) and short (minutes to hours) timescales.

79 ASTRONOMY AND ASTROPHYSICS↗