Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Statistical fluctuations”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Optimisation of the Search for CP-symmetry Violation at the Deep Underground Neutrino Experiment

The Deep Underground Neutrino Experiment (DUNE) is a next-generation long baseline experiment, which will be situated in South Dakota. Its detectors will utilise liquid-argon time projection chamber technology, which is able to capture neutrino interactions with an incredible spatial and calorimetric resolution. With what will become the world’s most intense neutrino beam, a highly capable near detector, and four (10kt fducial mass) far detector modules, DUNE will be able to achieve an ambitious physics programme. Most notably, DUNE will determine whether charge-parity symmetry is broken in neutrino oscillations - a finding that would have significant implications for the understanding of the matter-antimatter asymmetry in our Universe.This thesis presents the optimisation of a CP-violation analysis at DUNE using thePandora pattern-recognition software. The analysis assumes a 3.5 year exposure (1.36 × 1023 protons on target) to a neutrino and an antineutrino beam (7 year total). Only the predicted data of the far detector modules is used; near detector samples are not included. The initial sensitivity to CP-violation is found to be 3.8σ+0.9σ−1.1σ in an estimate that includes oscillation parameter uncertainties, systematic uncertainties and statistical fluctuations, and assumes a normal-ordering of the neutrino mass hierarchy. The performance of the Pandora event reconstruction is linked to that of the analysis, which is found to be limited by the reconstruction of the initial track-like region of electrons and photons. The ShowerRefinement algorithm is developed in response to this and its implementation into the analysis workflow results in an improved sensitivity to CP-violation of 4.6σ+0.9σ−1.0σ. With a perfected neutrino interaction vertex placement, this is further increased to 5.1σ+1.0σ−1.1σ.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

Long-term stability-degradation analysis of DUNE SiPMs in Liquid Nitrogen

We developed a cryogenic photon readout system to monitor arrays of DUNE SiPMs operating at various over-voltages in liquid nitrogen (LN2) for over three months. Photoelectron signals were read out simultaneously via weak capacitive coupling to the micro-FEMB, a 32-channel charge sensitive readout board designed for 77 K to 300 K operation in liquid argon time projection chambers (LAr-TPCs). A waveform snippet acquisition scheme was implemented to capture waveform signals whenever they exceeded a predetermined trigger level; consequently, empty waveforms were discarded. Selected SiPM parameters were monitored to detect any deviations beyond statistical fluctuations. While some parameters exhibited ~1σ variations over the 3-month test period, there was no evidence of drift being enhanced when operating at 4 V or 5 V compared to a 3 V over-voltage. However, a discernible PDE drop was observed across all channels, which warrants further investigation. In addition, we demonstrated that our system can simultaneously read out single photons from a selected group of SiPMs at room temperature.

Tsang, Thomas [Brookhaven]↗
DOE OSTI2024

Avalanche statistics of fluctuation-induced fluxes from the SLPM and the W7-AS stellarator

Measurements of fluctuating floating potentials and ion saturation currents at different radial locations in the Santander Linear Plasma Machine (Castellanos et al 2005 Plasma Phys. Control. Fusion47 2067) and at the edge of the W7-AS stellarator by means of radially movable Langmuir probes allow to infer the corresponding fluctuation-induced radial flux temporal series. Avalanche-like transport events are identified in the time series and statistically characterized in terms of avalanche size/duration/quiet-time distributions and size-duration scaling relations. Transport is diffusive in the inner and intermediate radial region of the SLPM r < r tr ≈ 2.6 cm, undergoing a transition at r tr , becoming non-diffusive in the outermost region of the device, r > r tr . Here, the results obtained at the edge of the W7-AS stellarator are similar to those found in SLPM for r > r tr , i.e. consistent with what would be expected for scale-free, self-similar plasma transport dynamics near a critical state.

Avalanches↗
DOE OSTI2022

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

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

36 MATERIALS SCIENCE↗
DOE OSTI2021

Outlier classification using autoencoders: application for fluctuation driven flows in fusion plasmas

Understanding the statistics of fluctuation driven flows in the boundary layer of magnetically confined plasmas is desired to accurately model the lifetime of the vacuum vessel components. Mirror Langmuir probes (MLPs) are a novel diagnostic that uniquely allow us to sample the plasma parameters on a time scale shorter than the characteristic time scale of their fluctuations. Sudden large-amplitude fluctuations in the plasma degrade the precision and accuracy of the plasma parameters reported by MLPs for cases in which the probe bias range is of insufficient amplitude. While some data samples can readily be classified as valid and invalid, we find that such a classification may be ambiguous for up to 40% of data sampled for the plasma parameters and bias voltages considered in this study. In this contribution, we employ an autoencoder (AE) to learn a low-dimensional representation of valid data samples. By definition, the coordinates in this space are the features that mostly characterize valid data. Ambiguous data samples are classified in this space using standard classifiers for vectorial data. In this way, we avoid defining complicated threshold rules to identify outliers, which require strong assumptions and introduce biases in the analysis. By removing the outliers that are identified in the latent low-dimensional space of the AE, we find that the average conductive and convective radial heat fluxes are between approximately 5% and 15% lower as when removing outliers identified by threshold values. For contributions to the radial heat flux due to triple correlations, the difference is up to 40%.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Statistical analysis of the fluctuations of an initial-state model with independently distributed hot spots

Abstract We determine the uncorrelated modes that characterize the fluctuations in a semi-realistic model for the initial state of high-energy nuclear collisions, consisting of hot spots whose positions are distributed independently. Varying the number of hot spots, their size, and the weights with which they contribute to the initial state, we find that the parameter that has the largest influence on the relative importance of the fluctuation modes is the source size, with more extended hot spots leading to a more marked predominance of the principal modes.

Borghini, Nicolas↗
DOE OSTI2023

Modified Fermi’s golden rule rate expressions

Fermi’s golden rule (FGR) serves as the basis for many expressions of spectroscopic observables and quantum transition rates. The utility of FGR has been demonstrated through decades of experimental confirmation. However, there still remain important cases where the evaluation of a FGR rate is ambiguous or ill-defined. Examples are cases where the rate has divergent terms due to the sparsity in the density of final states or time dependent fluctuations of system Hamiltonians. Strictly speaking, assumptions of FGR are no longer valid for such cases. However, it is still possible to define modified FGR rate expressions that are useful as effective rates. The resulting modified FGR rate expressions resolve a long standing ambiguity often encountered in using FGR and offer more reliable ways to model general rate processes. Finally, simple model calculations illustrate the utility and implications of new rate expressions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Forces from Stochastic Density Functional Theory under Nonorthogonal Atom-Centered Basis Sets

We develop a formalism for calculating forces on the nuclei within the linear-scaling stochastic density functional theory (sDFT) in a nonorthogonal atom-centered basis set representation (Fabian et al. Wiley Interdiscip. Rev.: Comput. Mol. Sci. 2019, 9, e1412, 10.1002/wcms.1412) and apply it to the Tryptophan Zipper 2 (Trp-zip2) peptide solvated in water. We use an embedded-fragment approach to reduce the statistical errors (fluctuation and systematic bias), where the entire peptide is the main fragment and the remaining 425 water molecules are grouped into small fragments. We analyze the magnitude of the statistical errors in the forces and find that the systematic bias is of the order of 0.065 eV/Å (~1.2 × 10 -3 E h / a 0 ) when 120 stochastic orbitals are used, independently of system size. This magnitude of bias is sufficiently small to ensure that the bond lengths estimated by stochastic DFT (within a Langevin molecular dynamics simulation) will deviate by less than 1% from those predicted by a deterministic calculation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

First-Principles Statistical Mechanics Study of Magnetic Fluctuations and Order–Disorder in the Spinel LiNi 0.5 Mn 1.5 O 4 Cathode

While significant magnetic interactions exist in lithium transition metal oxides, commonly used as Li-ion cathodes, the interplay between magnetic couplings, disorder, and redox processes remains poorly understood. In this work, we focus on the high-voltage spinel LiNi 0.5 Mn 1.5 O 4 (LNMO) cathode as a model system on which to apply a computational framework that uses first principles-based statistical mechanics methods to predict the finite temperature magnetic properties of materials and provide insights into the complex interplay between magnetic and chemical degrees of freedom. Density functional theory calculations on multiple distinct Ni–Mn orderings within the LNMO system, including the ordered ground-state structure (space group P4332), reveal a preference for a ferrimagnetic arrangement of the Ni and Mn sublattices due to strong antiferromagnetic superexchange interactions between neighboring Mn 4+ and Ni 2+ ions and ferromagnetic Mn–Mn and Ni–Ni couplings, as revealed by magnetic cluster expansions. These results are consistent with qualitative predictions using the Goodenough-Kanamori-Anderson rules. Simulations of the finite temperature magnetic properties of LNMO are conducted using Metropolis Monte Carlo. We find that a “semiclassical” Monte Carlo sampling method based on the Heisenberg Hamiltonian accurately predicts experimental magnetic transition temperatures observed in magnetometry measurements. This study highlights the importance of a robust computational toolkit that accurately captures the complex chemomagnetic interactions and predicts finite temperature magnetic behavior to help analyze experimental magnetic and magnetic resonance spectroscopy data acquired ex situ and operando.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Measuring the Speed of Sound Using Cumulants of Baryon Number

Here, we show that the values of the first three cumulants of the baryon number distribution can be used to calculate the isothermal speed of sound and its logarithmic derivative with respect to the baryon-number density. We discuss applications of this result to heavy-ion collision experiments and address possible challenges, including effects due to baryon-number conservation, differences between proton and baryon cumulants, and the influence of finite-number statistics on fluctuation observables in both ex periment and hadronic transport simulations. In particular, we investigate the relation between quantities calculated in infinite, continuous matter, and observables obtained in simulations using a finite number of particles.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Thermal-fluctuation effects on small-scale statistics in turbulent gas flow

Kolmogorov's theory of turbulence assumes that the small-scale turbulent structures in the energy cascade are universal and are determined by the energy dissipation rate and the kinematic viscosity alone. However, thermal fluctuations, absent from the continuum description, terminate the energy cascade near the Kolmogorov length scale. Here, we propose a simple superposition model to account for the effects of thermal fluctuations on small-scale turbulence statistics. Furthermore, for compressible Taylor–Green vortex flow, we demonstrate that the superposition model in conjunction with data from direct numerical simulation of the Navier–Stokes equations yields spectra and structure functions that agree with the corresponding quantities computed from the direct simulation Monte Carlo method of molecular gas dynamics, verifying the importance of thermal fluctuations in the dissipation range.

42 ENGINEERING↗
DOE OSTI2024

Hydrodynamic fluctuations near a Hopf bifurcation: Stochastic onset of vortex shedding behind a circular cylinder

Here, we investigate hydrodynamic fluctuations in the flow past a circular cylinder near the critical Reynolds number Re c for the onset of vortex shedding. Starting from the fluctuating Navier-Stokes equations, we perform a perturbation expansion around Re c to derive analytical expressions for the statistics of the fluctuating lift force. Molecular-level simulations using the direct simulation Monte Carlo method support the theoretical predictions of the lift power spectrum and amplitude distribution. Notably, we have been able to collect sufficient statistics at distances Re ⁡/ Re c – 1 = O ⁡(10 –3 ) from the instability that confirm the appearance of non-Gaussian fluctuations, and we observe that they are associated with intermittent vortex shedding. These results emphasize how unavoidable thermal-noise-induced fluctuations become dramatically amplified in the vicinity of oscillatory flow instabilities and that their onset is fundamentally stochastic.

42 ENGINEERING↗
DOE OSTI2022

Surface fluctuating hydrodynamics methods for the drift-diffusion dynamics of particles and microstructures within curved fluid interfaces

Here we introduce fluctuating hydrodynamics approaches on surfaces for capturing the drift-diffusion dynamics of particles and microstructures immersed within curved fluid interfaces of spherical shape. We take into account the interfacial hydrodynamic coupling, traction coupling with the surrounding bulk fluid, and thermal fluctuations. For fluid-structure interactions, we introduce Immersed Boundary Methods (IBM) and related Stochastic Eulerian-Lagrangian Methods (SELM) for curved surfaces. We use these approaches to investigate the statistics of surface fluctuating hydrodynamics and microstructures. For velocity autocorrelations, we find characteristic power-law scalings $τ^{-1}, τ^{-2}$, and plateaus can emerge. This depends on the physical regime associated with the geometry, surface viscosity, and bulk viscosity. This differs from the characteristic $τ^{-3/2}$ scaling for bulk three dimensional fluids. We develop theory explaining these observed power-laws associated with time-scales for dissipation within the fluid interface and coupling to the surrounding fluid. We then use our introduced methods to investigate a few example systems and roles of hydrodynamic coupling and thermal fluctuations including for the kinetics of passive particles and active microswimmers in curved fluid interfaces.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

A Goldstone boson equivalence for inflation

The effective field theory of single-field inflation characterizes the inflationary epoch in terms of a pattern of symmetry breaking. An operator acquires a time-dependent vacuum expectation value, defining a preferred spatial slicing. In the absence of dynamical gravity, the fluctuations around the time-dependent background are described by the Goldstone boson associated with this symmetry breaking process. With gravity, the Goldstone is eaten by the metric, becoming the scalar metric fluctuation. In this paper, we will show that in general single-field inflation, the statistics of scalar metric fluctuations are given by the statistics of this Goldstone boson decoupled from gravity up to corrections that are controlled as an expansion in slow-roll parameters. This even holds in the presence of additional parameters, like the speed of sound, that naively enhance the impact of the gravitational terms. In the process, we derive expressions for leading and sub-leading gravitational corrections to all-orders in the Goldstone boson.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Searching for the QCD Critical Point Along the Pseudo-critical/Freeze-out Line Using Padé-resummed Taylor Expansions of Cumulants of Conserved Charge Fluctuations

Using high-statistics datasets generated in (2+1)-flavor QCD calculations at finite temperature we construct estimators for the radius of convergence from an eighth order series expansion of the pressure as well as the number density. We show that the estimator for pressure and number density will be identical in the asymptotic limit. In the vicinity of the pseudo-critical temperature, T pc ≃ 156.5 MeV, we find the estimator of the radius of convergence to be µ B /T ≳ 3 for strangeness-neutral matter. We also present results for the pole structure of the Pad´e approximants for the pressure at non-zero values of the baryon chemical potential and show that the pole structure of the [4,4] Pad´e is consistent with not having a critical point at temperatures larger than 135 MeV and a baryon chemical potential smaller than µ B /T ~ 2.5.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Differentiating multi-MeV, multi-ion spectra with CR-39 solid-state nuclear track detectors

The development of high intensity petawatt lasers has created new possibilities for ion acceleration and nuclear fusion using solid targets. In such laser-matter interaction, multiple ion species are accelerated with broad spectra up to hundreds of MeV. To measure ion yields and for species identification, CR-39 solid-state nuclear track detectors are frequently used. However, these detectors are limited in their applicability for multi-ion spectra differentiation as standard image recognition algorithms can lead to a misinterpretation of data, there is no unique relation between track diameter and particle energy, and there are overlapping pit diameter relationships for multiple particle species. In this report, we address these issues by first developing an algorithm to overcome user bias during image processing. Second, we use calibration of the detector response for protons, carbon and helium ions (alpha particles) from 0.1 to above 10 MeV and measurements of statistical energy loss fluctuations in a forward-fitting procedure utilizing multiple, differently filtered CR-39, altogether enabling high-sensitivity, multi-species particle spectroscopy. To validate this capability, we show that inferred CR-39 spectra match Thomson parabola ion spectrometer data from the same experiment. Filtered CR-39 spectrometers were used to detect, within a background of ~ 2 × 10 11 sr -1 J -1 protons and carbons, (1.3 ± 0.7) × 10 8 sr -1 J -1 alpha particles from laser-driven proton-boron fusion reactions.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2024

Phase-space methods for neutrino oscillations: Extension to multibeams

The phase-space approach (PSA), which was originally introduced in Lacroix [] to describe neutrino flavor oscillations for interacting neutrinos emitted from stellar objects is extended to describe arbitrary numbers of neutrino beams. The PSA is based on mapping the quantum fluctuations into a statistical treatment by sampling initial conditions followed by independent mean-field evolution. A new method is proposed to perform this sampling that allows treating an arbitrary number of neutrinos in each neutrino beams. We validate the technique successfully and confirm its predictive power on several examples where a reference exact calculation is possible. We show that it can describe many-body effects, such as entanglement and dissipation induced by the interaction between neutrinos. Due to the complexity of the problem, exact solutions can only be calculated for rather limited cases, with a limited number of beams and/or neutrinos in each beam. The PSA approach considerably reduces the numerical cost and provides an efficient technique to accurately simulate arbitrary numbers of beams. Examples of PSA results are given here, including up to 200 beams with time-independent or time-dependent Hamiltonians. We anticipate that this approach will be useful to bridge exact microscopic techniques with more traditional transport theories used in neutrino oscillations. It will also provide important reference calculations for future quantum computer applications where other techniques are not applicable to classical computers. Published by the American Physical Society 2024

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Characterizing and Mitigating Intraday Variability: Reconstructing Source Structure in Accreting Black Holes with mm-VLBI

The extraordinary physical resolution afforded by the Event Horizon Telescope has opened a window onto the astrophysical phenomena unfolding on horizon scales in two known black holes, M87* and Sgr A*. However, with this leap in resolution has come a new set of practical complications. Sgr A* exhibits intraday variability that violates the assumptions underlying Earth aperture synthesis, limiting traditional image reconstruction methods to short timescales and data sets with very sparse (u, v) coverage. We present a new set of tools to detect and mitigate this variability. We develop a data-driven, model-agnostic procedure to detect and characterize the spatial structure of intraday variability. This method is calibrated against a large set of mock data sets, producing an empirical estimator of the spatial power spectrum of the brightness fluctuations. We present a novel Bayesian noise modeling algorithm that simultaneously reconstructs an average image and statistical measure of the fluctuations about it using a parameterized form for the excess variance in the complex visibilities not otherwise explained by the statistical errors. These methods are validated using a variety of simulated data, including general relativistic magnetohydrodynamic simulations appropriate for Sgr A* and M87*. We find that the reconstructed source structure and variability are robust to changes in the underlying image model. We apply these methods to the 2017 EHT observations of M87*, finding evidence for variability across the EHT observing campaign. The variability mitigation strategies presented are widely applicable to very long baseline interferometry observations of variable sources generally, for which they provide a data-informed averaging procedure and natural characterization of inter-epoch image consistency.

79 ASTRONOMY AND ASTROPHYSICS↗