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,225 records · Page 68

Muon Neutrino Disappearance in NOvA with a Deep Convolutional Neural Network Classifier

The NuMI Off-axis Neutrino Appearance Experiment (NOvA) is designed to study neutrino oscillation in the NuMI (Neutrinos at the Main Injector) beam. NOvA observes neutrino oscillation using two detectors separated by a baseline of 810 km; a 14 kt Far Detector in Ash River, MN and a functionally identical 0.3 kt Near Detector at Fermilab. The experiment aims to provide new measurements of Δm 2 and θ23 and has potential to determine the neutrino mass hierarchy as well as observe CP violation in the neutrino sector. Essential to these analyses is the classification of neutrino interaction events in NOvA detectors. Raw detector output from NOvA is interpretable as a pair of images which provide orthogonal views of particle interactions. A recent advance in the field of computer vision is the advent of convolutional neural networks, which have delivered top results in the latest image recognition contests. This work presents an approach novel to particle physics analysis in which a convolutional neural network is used for classification of particle interactions. The approach has been demonstrated to improve the signal efficiency and purity of the event selection, and thus physics sensitivity. Early NOvA data has been analyzed (2.74×10 20 POT, 14 kt equivalent) to provide new best- fit measurements of sin 2 (θ23) = 0.43 (with a statistically-degenerate compliment near 0.60) and |Δm2 | = 2.48 × 10 -3 eV 2 .

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Search for physics beyond the standard model in events with overlapping photons and jets

Results are reported from a search for new particles that decay into a photon and two gluons, in events with jets. Novel jet substructure techniques are developed that allow photons to be identified in an environment densely populated with hadrons. The analyzed proton-proton collision data were collected by the CMS experiment at the LHC, in 2016 at $ \sqrt{s}$ = 13 TeV, and correspond to an integrated luminosity of 35.9 fb -1 . The spectra of total transverse hadronic energy of candidate events are examined for deviations from the standard model predictions. No statistically significant excess is observed over the expected background. The first cross section limits on new physics processes resulting in such events are set. The results are interpreted as upper limits on the rate of gluino pair production, utilizing a simplified stealth supersymmetry model. The excluded gluino masses extend up to 1.7 TeV, for a neutralino mass of 200 GeV and exceed previous mass constraints set by analyses targeting events with isolated photons.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Effect of CLIQ on Training of HL-LHC Quadrupole Magnets

The high-luminosity LHC upgrade requires stronger than LHC low-beta quadrupole magnets to reach the luminosity goals of the project. The project is well advanced and HL-LHC quadrupole magnets are currently being commissioned in US Labs (MQXFA magnets) and CERN (MQXFB magnets). Those are the first Nb 3 Sn magnets to be used in any large particle accelerator. At development stages, many Nb 3 Sn accelerator sub-scale models showed relatively slow training and MQXFA magnets were projected to have low tens of quenches before reaching operational field. Recently it was shown that dedicated capacitor-based de-vices can affect Nb 3 Sn magnet training, and it was suggested that CLIQ, a capacitor-based device intended for quench protection, can do too. Here, the present paper investigates effects on training likely induced by CLIQ, using the base fact that only half the coils in a quadrupole experience upward current modulation at quench because of capacitor discharge. The study encompasses all MQXFA production magnets trained at BNL to date. No other high-statistics data from identical magnets (series) with CLIQ protection exist so far. Implications and opportunities stemming from data analysis are discussed and conclusions drawn.

43 PARTICLE ACCELERATORS↗

Statistics of link blockage due to cloud cover for free-space optical communications using NCDC surface weather observation data

Cloud opacity is one of the main atmospheric physical phenomena that can jeopardize the successful completion of an optical link between a spacecraft and a ground station. Hence, the site location chosen for a telescope used for optical communications must rely on knowledge of weather and cloud cover statistics for the geographical area where the telescope itself is located.

optical communication cloud cover↗

Impact of Fractionation and Dose in a Multivariate Model for Radiation-Induced Chest Wall Pain

Purpose: To determine the role of patient/tumor characteristics, radiation dose, and fractionation using the linear-quadratic (LQ) model to predict stereotactic body radiation therapy–induced grade ≥2 chest wall pain (CWP2) in a larger series and develop clinically useful constraints for patients treated with different fraction numbers. Methods and Materials: A total of 316 lung tumors in 295 patients were treated with stereotactic body radiation therapy in 3 to 5 fractions to 39 to 60 Gy. Absolute dose–absolute volume chest wall (CW) histograms were acquired. The raw dose-volume histograms (α/β = ∞ Gy) were converted via the LQ model to equivalent doses in 2-Gy fractions (normalized total dose, NTD) with α/β from 0 to 25 Gy in 0.1-Gy steps. The Cox proportional hazards (CPH) model was used in univariate and multivariate models to identify and assess CWP2 exposed to a given physical and NTD. Results: The median follow-up was 15.4 months, and the median time to development of CWP2 was 7.4 months. On a univariate CPH model, prescription dose, prescription dose per fraction, number of fractions, D83cc, distance of tumor to CW, and body mass index were all statistically significant for the development of CWP2. Linear-quadratic correction improved the CPH model significance over the physical dose. The best-fit α/β was 2.1 Gy, and the physical dose (α/β = ∞ Gy) was outside the upper 95% confidence limit. With α/β = 2.1 Gy, V{sub NTD99Gy} was most significant, with median V{sub NTD99Gy} = 31.5 cm{sup 3} (hazard ratio 3.87, P<.001). Conclusion: There were several predictive factors for the development of CWP2. The LQ-adjusted doses using the best-fit α/β = 2.1 Gy is a better predictor of CWP2 than the physical dose. To aid dosimetrists, we have calculated the physical dose equivalent corresponding to V{sub NTD99Gy} = 31.5 cm{sup 3} for the 3- to 5-fraction groups.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

How data science methods can improve the quality and efficiency of ICF and HEDP research

Data Science methods (many that are Bayesian based) are widely used in the physical sciences to estimate model parameters from experimental data, synthesize heterogeneous data, calibrate models, design experiments, and determine statistical significance of data. These methods provide a wealth of advantages over traditional analysis techniques because: 1) uncertainties are rigorously defined and propagated naturally through complex systems including covariance, 2) prior information is captured within the analysis framework (including rad-MHD and rad-hydro simulations), 3) competing models can be selected and/or ruled out using quantitative criteria, and 4) complex, heterogeneous data can be incorporated simultaneously. While these methods have been widely adopted as the gold standard in fields such as particle physics, astronomy, and biology, they have been slow to catch on in Inertial Confinement Fusion (ICF) and High Energy Density Physics (HEDP) research. Recently, several teams at LLNL, SNL, LANL, and the LLE have been exploring the use of these tools in their research and have found success. Here we propose that a concerted effort to consolidate these independent research efforts by developing and deploying common tools for use across the complex can revolutionize the way we approach data analysis, assimilation of theory and experiment, and decision making. The Bayesian formalism provides a means to accomplish this, but we are lacking certain infrastructure to make it happen on a large scale. Furthermore, once adopted, these techniques can be used to develop standards by which discoveries can be judged, similar to the so-called 5σ rule in high energy particle physics. Such standards may be used in the future to address the issue of unknown reproducibility in ICF and HED experiments caused by low shot rate and high cost per experiment. Our goals as a group are to advance the state of the art in HED measurement science by enabling: 1) better inferences from data with well-defined uncertainties, 2) better use of the data we have and continue to collect, 3) intelligent synthesis of data, 4) evaluation of the statistical significance of our data, and 5) informed decision making regarding the design of new experiments and instruments.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Energy dependence of light hypernuclei production in heavy-ion collisions from a coalescence and statistical-thermal model perspective

A comparison of light hypernuclei production, from Ultra-relativistic Quantum Molecular Dynamics (UrQMD)+coalescence and the thermal model, in heavy ion collisions over a wide range of beam energies and system sizes is presented. We find that both approaches provide generally similar results, with differences in specific details. Especially the ratios of hypertriton to Λ are affected by both the source radius Δr of the coalescence procedure as well as canonical effects. On the other hand, the double ratio S 3 is almost independent of canonical effects, which is in contrast to coalescence. Thus, both the beam energy dependence and centrality dependence of S 3 can be used to constrain the hypertriton source radius. To do so the currently available data are not yet sufficient. Elliptic flow is shown to be unaffected by the source size of the nuclei and an almost perfect mass scaling of the elliptic flow is observed. Furthermore, our predictions further suggest that the existence of the H-dibaryon(ΛΛ) seems ruled out by ALICE data.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Strongly coupled fluid-particle flows in vertical channels. I. Reynolds-averaged two-phase turbulence statistics

Simulations of strongly coupled (i.e., high-mass-loading) fluid-particle flows in vertical channels are performed with the purpose of understanding the fundamental physics of wall-bounded multiphase turbulence. The exact Reynolds-averaged (RA) equations for high-mass-loading suspensions are presented, and the unclosed terms that are retained in the context of fully developed channel flow are evaluated in an Eulerian–Lagrangian (EL) framework for the first time. A key distinction between the RA formulation presented in the current work and previous derivations of multiphase turbulence models is the partitioning of the particle velocity fluctuations into spatially correlated and uncorrelated components, used to define the components of the particle-phase turbulent kinetic energy (TKE) and granular temperature, respectively. The adaptive spatial filtering technique developed in our previous work for homogeneous flows [J. Capecelatro, O. Desjardins, and R. O. Fox, “Numerical study of collisional particle dynamics in cluster-induced turbulence,” J. Fluid Mech. 747, R2 (2014)] is shown to accurately partition the particle velocity fluctuations at all distances from the wall. Strong segregation in the components of granular energy is observed, with the largest values of particle-phase TKE associated with clusters falling near the channel wall, while maximum granular temperature is observed at the center of the channel. The anisotropy of the Reynolds stresses both near the wall and far away is found to be a crucial component for understanding the distribution of the particle-phase volume fraction. In Part II of this paper, results from the EL simulations are used to validate a multiphase Reynolds-stress turbulence model that correctly predicts the wall-normal distribution of the two-phase turbulence statistics.

42 ENGINEERING↗

An Illumination- and Temperature-Dependent Analytical Model for Copper Indium Gallium Diselenide (CIGS) Solar Cells

In this study, we present a physics-based analytical model for copper indium gallium diselenide (CIGS) solar cells that describes the illumination- and temperature-dependent current-voltage (I-V) characteristics and accounts for the statistical shunt variation of each cell. The model is derived by solving the drift-diffusion transport equation so that its parameters are physical and, therefore, can be obtained from independent characterization experiments. The model is validated against CIGS I-V characteristics as a function of temperature and illumination intensity. This physics-based model can be integrated into a large-scale simulation framework to optimize the performance of solar modules, as well as predict the long-term output yields of photovoltaic farms under different environmental conditions.

14 SOLAR ENERGY↗

A method for obtaining a statistically stationary turbulent free shear flow

The long-term goal of the current research is the study of Large-Eddy Simulation (LES) as a tool for aeroacoustics. New algorithms and developments in computer hardware are making possible a new generation of tools for aeroacoustic predictions, which rely on the physics of the flow rather than empirical knowledge. LES, in conjunction with an acoustic analogy, holds the promise of predicting the statistics of noise radiated to the far-field of a turbulent flow. LES's predictive ability will be tested through extensive comparison of acoustic predictions based on a Direct Numerical Simulation (DNS) and LES of the same flow, as well as a priori testing of DNS results. The method presented here is aimed at allowing simulation of a turbulent flow field that is both simple and amenable to acoustic predictions. A free shear flow is homogeneous in both the streamwise and spanwise directions and which is statistically stationary will be simulated using equations based on the Navier-Stokes equations with a small number of added terms. Studying a free shear flow eliminates the need to consider flow-surface interactions as an acoustic source. The homogeneous directions and the flow's statistically stationary nature greatly simplify the application of an acoustic analogy.

Timson, Stephen F.↗

Statistical Analysis of the First Passage Path Ensemble of Jump Processes

The transition mechanism of jump processes between two different subsets in state space reveals important dynamical information of the processes and therefore has attracted considerable attention in the past years. In this paper, we study the first passage path ensemble of both discrete-time and continuous-time jump processes on a finite state space. The main approach is to divide each first passage path into nonreactive and reactive segments and to study them separately. The analysis can be applied to jump processes which are non-ergodic, as well as continuous-time jump processes where the waiting time distributions are non-exponential. In the particular case that the jump processes are both Markovian and ergodic, our analysis elucidates the relations between the study of the first passage paths and the study of the transition paths in transition path theory. We provide algorithms to numerically compute statistics of the first passage path ensemble. The computational complexity of these algorithms scales with the complexity of solving a linear system, for which efficient methods are available. Several examples demonstrate the wide applicability of the derived results across research areas.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Optimizing time integration for accurate recovery of shockwave interface location in radiography

We present simulations and experiments of time integrated radiographic imaging of a moving 1D shock wave front and a quantitative method for determining the statistical error in locating the shock front as a function of integration time and noise in the radiograph. We discuss the trade-off between increasing motion blur, which leads to decreased shock front location certainty, and increasing signal-to-noise, which leads to improved image quality with increasing integration time. We find an optimum integration time between a short integration time, where noise limits the error, and a long integration time, where motion blurring limits the error. This methodology can be used to tune experimental configurations to obtain the highest quality radiograph for a given experimental configuration.

Bremsstrahlung↗

Chromospheric Evaporation and Decimetric Radio Emission in Solar Flares

We have discovered decimetric signatures of the chromospheric evaporation process. Evidence for the radio detection of chromospheric evaporation is based on the radio-inferred values of (1) the electron density, (2) the propagation speed, and (3) the timing, which are found to be in good agreement with statistical values inferred from the blueshifted Ca xix soft X-ray line. The physical basis of our model is that free-free absorption of plasma emission is strongly modified by the steep density gradient and the large temperature increase in the upflowing flare plasma. The steplike density increase at the chromospheric evaporation front causes a local discontinuity in the plasma frequency, manifested as almost infinite drift rate in decimetric type III bursts. The large temperature increase of the upflowing plasma considerably reduces the local free-free opacity (due to the T-(exp -3/2) dependence) and thus enhances the brightness of radio bursts emitted at the local plasma frequency near the chromospheric evaporation front, while a high-frequency cutoff is expected in the high-density regions behind the front, which can be used to infer the velocity of the upflowing plasma. From model calculations we find strong evidence that decimetric bursts with a slowly drifting high-frequency cutoff are produced by fundamental plasma emission, contrary to the widespread belief that decimetric bursts are preferentially emitted at the harmonic plasma level. We analyzed 21 flare episodes from 1991-1993 for which broadband (100-3000 MHz) radio dynamic spectra from Phoenix, hard X-ray data from BATSE/CGRO, and soft X-ray data from GOES were available. We detected slowly drifting high-frequency cutoffs between 1.1 and 3.0 GHz, with drift rates of -41 +/- 32 MHz/s, extending over time intervals of 24 +/- 23 s. Developing a density model for type III-emitting flare loops based on the statistically observed drift rate of type III bursts by Alvarez & Haddock, we infer velocities of up to 360 km/s for the upflowing plasma, with an average of v(sub CE) = 236 +/- 130 km /s for episodes with 5-15 s duration. The mean electron density of the upflowing plasma is n(sub e) = 5.2(+/-3.1) x 10(exp 10) /cu cm when it is first detected in radio, at coronal altitudes of h(sub 0) = 9.2 +/- 2.3 Mm.

Aschwanden, Markus J.↗

Molecular-based analysis of nanoparticle solvation: Classical density functional approach

Proper statistical mechanics understanding of nanoparticle solvation processes requires an accurate description of the molecular structure of the solvent. Achieving this goal with standard molecular dynamics (MD) simulation methods is challenging due to large length scales. An alternative approach to this problem can be formulated using classical density functional theory (cDFT), where a full configurational description of the positions of all the atoms is replaced by collective atomic site densities in the molecule. Using an example of the negatively charged silica-like system in an aqueous polar environment represented by a two-site water model, we demonstrate here that cDFT can reproduce MD data at a fraction of the computational cost. An important implication of this result is the ability to understand how the solvent molecular features may affect the system’s properties at the macroscopic scale. A concrete example highlighted in this work is the analysis of nanoparticle interactions with sizes of up to 100 nm in diameter.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Preface to special topic: The High Repetition Rate Frontier in High-Energy-Density Physics

High-repetition-rate (HRR) experiments can collect large datasets with high temporal, spatial, and/or parametric resolution or large numbers of repeat measurements for statistics. HRR experiments also enable new experimental designs, including active feedback control loops and novel diagnostics, that can improve the reproducibility as well as the quantity of measurements. Together, these attributes make HRR experiments ideal for performing high-quality repeatable science. Until recently, these techniques have not been applied to high-energy-density–physics (HEDP) experiments, which are typically restricted to repetition rates of a few per day. However, recent advancements in lasers, pulsed-power drivers, target fabrication, and diagnostics are starting to change this fact, opening an exciting new frontier of HRR HEDP experiments. A mini-conference on this subject at the 2021 meeting of the American Physical Society Division of Plasma Physics brought together members of this growing community. As a result, the “High Repetition Rate Frontier in High-Energy-Density Physics” special topic in Physics of Plasmas highlights current progress in this exciting area.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Statistical Study of Primary Knock-on Atoms under Neutron Irradiation in RMC Code

Materials in high flux neutron condition will suffer from an evolution of damage which ultimately restricts the lifetime of the fuel clad and reactor pressure vessel. Atoms in irradiated material transform to primary knock-on atoms (PKAs) by absorbing the kinetic energy transferred from neutrons, and the energy released from nuclear reactions. The secondary collision cascades which are induced by PKAs create point defects within the material, with residual vacancies and interstitials eventually causing swelling and embrittlement that degrades mechanical behavior. For a long time, Displacements per Atom (DPA) index has been identified as a measure of radiation damage in material. Nowadays, several advanced methods, such as Molecular Dynamics (MD) and Kinetic Monte Carlo (KMC) simulations, have been developed for more accurately simulating the evolution of the damage. However, all of these techniques need more comprehensive information of PKA, such as the type, energy, direction and spatial distribution. The limitations of DPA emerge and a method should be provided to simulate the generative process of PKAs. Recently, a deterministic analysis code SPECTRA-PKA has been developed by Gilbert to calculate the energy spectra of almost all types PKAs by solving the group-to-group cross-section matrices equation which includes neutron scattering, light charged particles, and recoils of the residual nuclei. But this deterministic analysis code could not give the direction, spatial and time distribution of PKAs. Based on the Geant4 program, Cui Zhenguo has simulated the irradiation by neutrons with 1 MeV energy in zirconium, while only the two adjacent elastic collisions between injected neutron and target atoms are taken into consideration. In this paper, Monte Carlo method to calculate PKAs information during neutron transport in RMC code is introduced. RMC is a continuous-energy Reactor Monte Carlo neutron and photon transport code being developed by Department of Engineering Physics at Tsinghua University, Beijing. PKAs in zirconium are calculated in fixed neutron source mode. Statistical information of PKAs in PWR fuel cell and SFR fuel cell is shown and compared. The results show that Monte Carlo method can provide more comprehensive information of PKAs. (authors)

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Spectral broadening from turbulence in multiscale lower hybrid current drive simulations

Abstract The scattering of lower hybrid (LH) waves due to scrape-off layer filaments is investigated. It is revealed that scattering can account for the LH spectral gap without any ad hoc modification to the wave-spectrum. This is shown using a multiscale simulation approach which allows, for the first time, the inclusion of full-wave scattering physics in ray-tracing/Fokker–Planck calculations. In this approach, full-wave scattering probabilities are calculated for a wave interacting with a statistical ensemble of filaments. These probabilities are coupled to ray-tracing equations using radiative transfer theory. This allows the modeling of scattering along the entire ray-trajectory, which can be important in the multi-pass regime. Simulations are conducted for LH current drive (LHCD) in Alcator C-Mod, resulting in excellent agreement with experimental current and hard x-ray profiles. A region in filament parameter space is identified in which the impact of scattering on LHCD is saturated. Such a state coincides with experimental LHCD measurements, suggesting saturation indeed occurs in C-Mod, and therefore the exact statistical properties of the filaments are not important.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗