Engineering Papers⌕ Search

SEARCH · Engineering Papers

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,117 records · Page 62

End-Use Load Profiles for the U.S. Building Stock

The End-Use Load Profiles for the U.S. Building Stock project uses New ResStock and ComStock models to statistically represent the energy use of U.S. buildings. The project's hybrid approach combines best-available ground truth data, such as submetering studies and statistical disaggregation of whole-building interval meter data, with the reach, cost-effectiveness, and granularity of physics-based and data-driven building stock modeling to deliver a nationally-comprehensive data set at a fraction of the historical cost.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Simulation of Particle-Material Interactions

This paper gives an overview of the particle transport theory essentials, the basics of particle-material interaction simulation, physical quantities needed to simulate particle transport and interactions in materials, Monte Carlo simulation flow, response of additive detectors, statistical weights and other techniques to minimize statistical errors. Effects in materials under irradiation, materials response related to component lifetime and performance are considered with a focus on high-energy and high-power accelerator applications. Implementation of simulation of particle-material interactions in the modern Monte Carlo codes along with the code s main features and results of recent benchmarking are described.

43 PARTICLE ACCELERATORS↗

Analysis of the multicomponent relativistic Boltzmann equation for electron scattering in big bang nucleosynthesis

Big-bang nucleosynthesis (BBN) is valuable as a means to constrain the physics of the early universe and it is the only probe of the radiation-dominated epoch. A fundamental assumption in BBN is that the nuclear velocity distributions obey Maxwell-Boltzmann (MB) statistics as they do in stars. Specifically, the BBN epoch is characterized by a dilute baryon plasma for which the velocity distribution of nuclei is mainly determined by the dominant Coulomb elastic scattering with mildly relativistic electrons. One must therefore deduce the momentum distribution for reacting nuclei from the multicomponent relativistic Boltzmann equation. However, the full multicomponent relativistic Boltzmann equation has only recently been analyzed and its solution has only been worked out in special cases. Moreover, a variety of schemes have been proposed that introduce nonthermal components into the BBN environment which can alter the thermal distribution of reacting nuclei. Here, we construct the relativistic Boltzmann equation in the context of BBN. We also derive a Langevin model and perform relativistic Monte-Carlo simulations which clarify the baryon distribution during BBN and can be used to analyze any relaxation from a nonthermal injection. Further, we show by these analyses that the thermalization process leads to a nuclear distribution function that remains very close to MB statistics even during the most relativistic environment relevant to BBN. Hence, the predictions of standard BBN remain unchanged.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Second quantization in bit-string physics

Using a new fundamental theory based on bit-strings, a finite and discrete version of the solutions of the free one particle Dirac equation as segmented trajectories with steps of length h/mc along the forward and backward light cones executed at velocity +/- c are derived. Interpreting the statistical fluctuations which cause the bends in these segmented trajectories as emission and absorption of radiation, these solutions are analogous to a fermion propagator in a second quantized theory. This allows us to interpret the mass parameter in the step length as the physical mass of the free particle. The radiation in interaction with it has the usual harmonic oscillator structure of a second quantized theory. How these free particle masses can be generated gravitationally using the combinatorial hierarchy sequence (3,10,137,2(sup 127) + 136), and some of the predictive consequences are sketched.

Noyes, H. Pierre↗

Heterogeneous physicochemistry of the winter polar stratosphere

Present chemical theories of the Antarctic ozone hole assume that heterogeneous reactions involving polar stratospheric clouds (PSCs) are the precursor of springtime ozone depletions. However, none of the theories quantify the rates of proposed heterogeneous processed, and none utilize the extensive data base on PSC's. Thus, all of the theories must be considered incomplete until the heterogeneous mechanisms are properly defined. A unified treatment developed of the cloud related processes, both physical and chemical, and the importance of these processes using observation data is calibrated. The rates are compared competitive heterogeneous processes to place reasonable limits on critical mechanisms such as the denitrification and dechlorination of the polar winter stratosphere. Among the subjects addressed here are the physical/chemical properties of PSC's including their relevant microphysical, optical and compositional characteristics, mass transfer rates of gaseous constituents to cloud particles, adsorption, accommodation and sticking coefficients on cloud particles, time constants for condensation, absorption and other microphysical processes, effects of solubility and vapor pressure on cloud composition, the statistics of cloud processing of chemically active condensible species, rate limiting steps in heterogeneous chemical reactions, and the nonlinear dependence of ozone loss on physical and chemical parameters.

Turco, R. P.↗

Equations of State in Computational Physics: What, Why, and How

Summary: 1. Equations of state furnish thermodynamic properties of a material. 2. Usually involve P, p (or V), and T, or P, r, and E 3. Necessary ingredient to solve compressible flow problems. 4. Varieties: i. Analytical expressions (ideal gas, Grüneisen) – fast, convenient for verification and test problems ii. Graphs and charts – good for hand calculations iii. Tables and databases – generally most accurate, require interpolation i. Steam tables ii. NIST database iii. SESAME EOS Library 5. Constructed by using semi-empirical models from Statistical Mechanics and fitting to available data

36 MATERIALS SCIENCE↗

Data Selection Improvement For MicroBooNE

Data selection is an extremely important part of data analysis for any experiment. Finding a physics result is often the result of sifting through a massive amount of data, keeping data that we believe to be signal, and throwing out data we do not. This process is called data selection. Creating a selection algorithm is an intensive process that must balance keeping enough data to have statistics and maximizing the signal purity of that data. We also need to choose the right reconstruction method, a tool to take raw data from the detector and convert it into physics results. In this study, we used three different reconstruction tools, Pandora, WireCell, and LANTERN, for the MicroBooNE experiment in conjunction to improve the selection algorithm for analysis. For the case of this study, we look into the charged current N proton 0 pions (CCNp0$\pi$) interaction channel. This is the dominant channel for the Short Baseline Neutrino (SBN) program and is expected to be a large contributor to the Deep Underground Neutrino Experiment (DUNE). We first investigated each of the three tools to find out more about their strengths and weaknesses as reconstructions, and compared them to the truth information directly from the MicroBooNE simulation pipeline. We then put together a direct comparison of the three methods to find which method or combination of methods would return the best result for us. While the study is ongoing, we have learned a lot about data selection for the experiment and the differences between the reconstruction tools.

Dillon, Brayden [Fermilab]↗

Methodology for physics-informed generation of synthetic neutron time-of-flight measurement data

Accurate neutron cross section data are a vital input to the simulation of nuclear systems for a wide range of applications from energy production to national security. The evaluation of experimental data is a key step in producing accurate cross sections. There is a widely recognized lack of reproducibility in the evaluation process due to its artisanal nature and therefore there is a call for improvement within the nuclear data community. This can be realized by automating/standardizing viable parts of the process, namely, parameter estimation by fitting theoretical models to experimental data. This automation effort could greatly benefit from a synthetic data resource. This work leverages problem-specific physics, Monte Carlo sampling, and a general methodology for data synthesis to generate unlimited, labelled experimental cross-section data that is statistically indistinguishable to the observed data. Heuristic and, where applicable, rigorous statistical comparisons to observed data support this claim. The demonstration is based on/limited to transmission measurements at Rensselaer Polytechnic Institute (RPI) and energy-differential cross sections in the resolved resonance region (RRR). An open-source software is published alongside this article that executes the complete methodology to produce high-utility synthetic datasets. The goal of this work is to provide an approach and corresponding tool that will allow the evaluation community to begin exploring more data-driven, ML-based solutions to long-standing challenges in the field.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Comparison of CERES-MODIS Stratus Cloud Properties with Ground-Based Measurements at the DOE ARM Southern Great Plains Site

Overcast stratus cloud properties derived for the Clouds and the Earth's Radiant Energy system (CERES) Project using Terra and Aqua Moderate Resolution Imaging Spectroradiometer (MODIS) data are compared with observations taken at the Atmospheric Radiation Measurement (ARM) Southern Great Plains site from March 2000 through December 2004. Retrievals from ARM surface-based data were averaged over a 1-hour interval centered at the time of each satellite overpass, and the CERES-MODIS cloud properties were averaged within a 30-km x 30 km box centered on the ARM SGP site. Two datasets were analyzed: all of the data (ALL) which include multilayered, single-layered, and slightly broken stratus decks and a subset, single-layered unbroken decks (SL). The CERES-MODIS effective cloud heights were determined from effective cloud temperature using a lapse rate method with the surface temperature specified as the 24-h mean surface air temperature. For SL stratus, they are, on average, within the ARM radar-lidar estimated cloud boundaries and are 0.534 +/- 0.542 km and 0.108 +/- 0.480 km lower than the cloud physical tops and centers, respectively, and are comparable for day and night observations. The mean differences and standard deviations are slightly larger for ALL data, but not statistically different to those of SL data. The MODIS-derived effective cloud temperatures are 2.7 +/- 2.4 K less than the surface-observed SL cloud center temperatures with very high correlations (0.86-0.97). Variations in the height differences are mainly caused by uncertainties in the surface air temperatures, lapse rates, and cloud-top height variability. The biases are mainly the result of the differences between effective and physical cloud top, which are governed by cloud liquid water content and viewing zenith angle, and the selected lapse rate, -7.1 K km(exp -1). Based on a total of 43 samples, the means and standard deviations of the differences between the daytime Terra and surface retrievals of effective radius r(sub e), optical depth, and liquid water path for SL stratu are 0.1 +/- 1.9 micrometers (1.2 +/- 23.5%), -1.3 +/- 9.5 (-3.6 +/-26.2%), and 0.6 +/- 49.9 gm (exp -2) (0.3 +/- 27%), respectively, while the corresponding correlation coefficients are 0.44, 0.87, and 0.89. For Aqua, they are 0.2 +/- 1.9 micrometers (2.5 +/- 23.4%), 2.5 +/- 7.8 (7.8 +/- 24.3%), and 28.1 +/- 52.7 gm (exp -2) (17.2 +/- 32.2%), as well as 0.35, 0.96, and 0.93 from a total of 21 cases. The results for ALL cases are comparable. Although a bias in R(sub e) was expected because the satellite retrieval of effective radius only represents the top of the cloud, the surface-based radar retrievals revealed that the vertical profile of r(sub e) is highly variable with smaller droplets occurring at cloud top in some cases. The larger bias in optical depth and liquid water path for Aqua is due, at least partially, to differences in the Terra and Aqua MODIS visible channel calibrations. methods for improving the cloud-top height and microphysical property retrievals are suggested.

Dong, Xiquan↗

Evaluation of an efficient statistical transport model for turbulent droplet dispersion in dilute combusting sprays

Evaluation of an efficient statistical transport model for turbulent droplet dispersion is made for a dilute spray of methanol droplets injected into a turbulent, axisymmnetric methane-fuelled diffusion flame burning in stagnant air. In the dispersion model, a computational parcel representing a group of physical particles (droplets) is considered to have a normal probability density function (pdf) in space. The mean is determined by Lagrangian tracking through a sequence of stochastically generated turbulent eddies and the variance is evaluated from a statistical formulation based on the linearized particle equations of motion. The basic validity of this model is established through comparison with available experimental data and with theoretical predictions using a conventional stochastic direct modeling approach. The conclusion of the evaluation is that the proposed dispersion model compares favorably with experimental data and provides a valid technique for simulating turbulent combusting sprays with significant computational savings over conventional methods.

Litchford, Ron J.↗

Overview of the Tolerance Limit Calculations with Application to TSURFER

To establish confidence in the results of computerized physics models, a key regulatory requirement is to develop a scientifically defendable process. The methods employed for confidence, characterization, and consolidation, or C3, are statistically involved and are often accessible only to avid statisticians. This manuscript serves as a pedagogical presentation of the C3 process to all stakeholders—including researchers, industrial practitioners, and regulators—to impart an intuitive understanding of the key concepts and mathematical methods entailed by C3. The primary focus is on calculation of tolerance limits, which is the overall goal of the C3 process. Tolerance limits encode the confidence in the calculation results as communicated to the regulator. Understanding the C3 process is especially critical today, as the nuclear industry is considering more innovative ways to assess new technologies, including new reactor and fuel concepts, via an integrated approach that optimally combines modeling and simulation and minimal targeted validation experiments. This manuscript employs intuitive, analytical, numerical, and visual representations to explain how tolerance limits may be calculated for a wide range of configurations, and it also describes how their values may be interpreted. Various verification tests have been developed to test the calculated tolerance limits and to help delineate their values. The manuscript demonstrates the calculation of tolerance limits for TSURFER, a computer code developed by the Oak Ridge National Laboratory for criticality safety applications. The goal is to evaluate the tolerance limit for TSURFER-determined criticality biases to support the determination of upper, subcritical limits for regulatory purposes.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

New physics in single resonant top quarks

Abstract Searches for new physics in the top quark sector are of great theoretical interest, yet some powerful avenues for discovery remain unexplored. We characterize the expected statistical power of the LHC dataset to constrain the single production of heavy top partnersTdecaying to a top quark and a photon or a top quark and a gluon. We describe an effective interaction which could generate such production, though the limits apply to a range of theoretical models. We find sensitivity to cross sections in the 10 2 − 10 5 fb range, forTmasses between 300 and 1000 GeV, depending on decay mode.

Physics↗

Machine Learning to Select Experiments Driven by Fundamental Science and Applications for Targeted Nuclear Data Improvement

This work describes a blueprint for a process that accelerates progress in science by quantitatively answering the following question: What is the optimal combination of fundamental-science and application-driven experiments to maximally reduce pertinent data uncertainties? Answering this question entails solving a high-dimensional and complex optimization problem that is best solved with advanced statistic techniques often classified as machine learning. We apply this process within the framework of nuclear data with the aim to select an experiment combination that will reduce uncertainties in 239 Pu nuclear data for neutron energies between 1 and 600 keV. In this field, fundamental-physics driven data, called differential, look at one nuclear physics observable at a time. They are contrasted to application-driven, integral, data where one or few resulting values inform a broad set of nuclear data across several nuclides and energies. The candidates for integral experiments are criticality measurements that were refined by a genetic algorithm to be maximally sensitive to 239 Pu fission cross sections in the desired energy range. Twenty-three candidate differential experiments were investigated and span multiple nuclear physics observables (e.g., total, capture cross sections) for isotopes appearing in the integral experiments. The optimal combination among these candidate experiments was investigated via generalized least squares fitting, augmented with Gaussian processes to ameliorate statistical irregularities in data, and the D-optimality criterion. The latter evaluates for each pair of candidates the joint reduction in uncertainties of all 12200 nuclear data appearing in the integral experiments compared to the knowledge we have from 168 past experiments, theory, and nuclear data. We chose as differential measurements those that investigate 63 Cu and 239 Pu total cross sections, based on D-optimality rank and feasibility constraints. Two integral (criticality) experiments were selected: An experiment with Al 2 ⁢O 3 and graphite interleaved with Pu and a thick Cu reflector explores 1–30 keV, while we target the 30–600 keV range with an experiment that swaps boron in place of graphite with a different geometry.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Bayesian estimation of the $S$ factor and thermonuclear reaction rate for 16 O(p, γ) 17 F

The 16 O(p, γ) 17 F reaction is the slowest hydrogen-burning process in the CNO mass region. Its thermonuclear rate sensitively impacts predictions of oxygen isotopic ratios in a number of astrophysical sites, including AGB stars. The reaction has been measured several times at low bombarding energies using a variety of techniques. The most recent evaluated experimental rates have a reported uncertainty of about 7.5% below 1 GK. However, the previous rate estimate represents a best guess only and was not based on rigorous statistical methods. We apply a Bayesian model to fit all reliable 16 O(p, γ) 17 F cross section data, and take into account independent contributions of statistical and systematic uncertainties. The nuclear reaction model employed is a single-particle potential model involving a Woods-Saxon potential for generating the radial bound state wave function. The model has three physical parameters, the radius and diffuseness of the Woods-Saxon potential, and the asymptotic normalization coefficients (ANCs) of the final bound state in 17 F. Here, we find that performing the Bayesian S -factor fit using ANCs as scaling parameters has a distinct advantage over adopting spectroscopic factors instead. Based on these results, we present the first statistically rigorous estimation of experimental 16 O(p, γ) 17 F reaction rates, with uncertainties (±4.2%) of about half the previously reported values.

6 ≤ A ≤ 19↗

Comments on foliated gauge theories and dualities in 3+1d

We investigate the properties of foliated gauge fields and construct several foliated field theories in 3+1d that describe foliated fracton orders both with and without matter, including the recent hybrid fracton models. These field theories describe Abelian or non-Abelian gauge theories coupled to foliated gauge fields, and they fall into two classes of models that we call the electric models and the magnetic models. We show that these two classes of foliated field theories enjoy a duality. We also construct a model (using foliated gauge fields and an exactly solvable lattice Hamiltonian model) for a subsystem-symmetry protected topological (SSPT) phase, which is analogous to a one-form symmetry protected topological phase, with the subsystem symmetry acting on codimension-two subregions. We construct the corresponding gauged SSPT phase as a foliated two-form gauge theory. Some instances of the gauged SSPT phase are a variant of the X-cube model with the same ground state degeneracy and the same fusion, but different particle statistics.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Simulations and models for the Richtmyer– Meshkov instability with broadband perturbations

The Richtmyer–Meshkov instability (RMI) is shock driven and affects many phenomena from inertial fusion to supernova explosions. The behavior of single-modes in the RMI has been studied extensively but less is known with the broadband perturbations that occur in applications. Here, we describe extensive numerical simulations and modeling of the RMI with broadband perturbations with an initial power spectrum of the form P∝ k m , where k is the wavenumber. The hydrodynamic condition is the same as in the $\theta$-Group Collaboration [Thornber et al., “Late-time growth rate, mixing, and anisotropy in the multimode narrowband Richtmyer–Meshkov instability: The $\theta$-group collaboration,” Phys. Fluids 29, 105107 (2017)] with a Mach ~1.86 shock and gamma-law = 5/3 fluids with Atwood number A = 0.5. The bubble amplitude h B is found to grow in two stages. Initially, h B undergoes a phase-inversion and grows linearly in time (t) at a rate consistent with a linear theory. Asymptotically, h B grows as a power law $~t^{\theta}$ when $\langle{k}\rangle{|h_B|}$ > O(1), where $\langle{k}\rangle$ is the average wavenumber for the initial spectrum. The RMI behavior in simulations and modeling agree over a wide range of exponent m, spectral width, initial amplitude, and time. The agreement is quantified objectively using statistical analysis.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Electric charge and strangeness-dependent directed flow splitting of produced quarks in Au+Au collisions

We report directed flow (v 1 ) of multistrange baryons (Ξ and Ω) and improved v1 data for $K¯$, $\bar{p}$, $\bar{Λ}$ and φ in Au+Au collisions at $\sqrt{s_{NN}} =27$ and 200 GeV from the STAR at the Relativistic Heavy Ion Collider (RHIC). We focus on particles whose constituent quarks are not transported from beam rapidity rather produced in the collisions. In midcentral collisions, we observe a coalescence sum rule for hadron combinations with identical quark content and a difference (“splitting”) in the slope of v 1 vs. rapidity for combinations having nonidentical quark content. The splitting strength appears to increase with the electric charge difference and strangeness content difference of the constituent quarks in the combinations, consistent with an electromagnetic effect. The peripheral collision statistics are insufficient to draw firm conclusions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗