Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Stochastic”

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

Implicit learning of convective organization explains precipitation stochasticity

Accurate prediction of precipitation intensity is crucial for both human and natural systems, especially in a warming climate more prone to extreme precipitation. Yet, climate models fail to accurately predict precipitation intensity, particularly extremes. One missing piece of information in traditional climate model parameterizations is subgrid-scale cloud structure and organization, which affects precipitation intensity and stochasticity at coarse resolution. Here, using global storm-resolving simulations and machine learning, we show that, by implicitly learning subgrid organization, we can accurately predict precipitation variability and stochasticity with a low-dimensional set of latent variables. Using a neural network to parameterize coarse-grained precipitation, we find that the overall behavior of precipitation is reasonably predictable using large-scale quantities only; however, the neural network cannot predict the variability of precipitation (R 2 ~ 0.45) and underestimates precipitation extremes. The performance is significantly improved when the network is informed by our organization metric, correctly predicting precipitation extremes and spatial variability (R 2 ~ 0.9). The organization metric is implicitly learned by training the algorithm on a high-resolution precipitable water field, encoding the degree of subgrid organization. The organization metric shows large hysteresis, emphasizing the role of memory created by subgrid-scale structures. We demonstrate that this organization metric can be predicted as a simple memory process from information available at the previous time steps. These findings stress the role of organization and memory in accurate prediction of precipitation intensity and extremes and the necessity of parameterizing subgrid-scale convective organization in climate models to better project future changes of water cycle and extremes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2012

Stochastic orbit loss of neutral beam ions from NSTX due to toroidal Alfvén eigenmode avalanches

Short toroidal Alfvén eigenmode (TAE) avalanche bursts in the National Spherical Torus Experiment (NSTX) cause a drop in the neutron rate and could also cause a loss of neutral beam ions at or near the full injection energy over an extended range of pitch angles. The simultaneous loss of wide ranges of pitch angle suggests stochastic transport of the beam ions takes place. When beam ion orbits are followed with a guiding centre code that incorporates the plasma's magnetic equilibrium plus the measured modes, the predicted ranges of lost pitch angle are like those seen in the experiment, with distinct populations of trapped and passing orbits lost. These correspond to domains where the stochasticity extends in the orbit phase space from the region of beam ion deposition to the loss boundary and the trajectories along which modes may transport particles extend from the deposition volume to the loss boundary.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Stochastic and a posteriori optimization to mitigate coil manufacturing errors in stellarator design

Abstract It was recently shown in Wechsung et al (2022 Proc. Natl Acad. Sci. USA 119 e2202084119) that there exist electromagnetic coils that generate magnetic fields, which are excellent approximations to quasi-symmetric fields and have very good particle confinement properties. Using a Gaussian process-based model for coil perturbations, we investigate the impact of manufacturing errors on the performance of these coils. We show that even fairly small errors result in noticeable performance degradation. While stochastic optimization yields minor improvements, it is not possible to mitigate these errors significantly. As an alternative to stochastic optimization, we then formulate a new optimization problem for computing optimal adjustments of the coil positions and currents without changing the shapes of the coil. These a-posteriori adjustments are able to reduce the impact of coil errors by an order of magnitude, providing a new perspective for dealing with manufacturing tolerances in stellarator design.

Physics↗
DOE OSTI2019

Stochastic particle production in a de Sitter background

Here, we explore non-adiabatic particle production in a de Sitter universe for a scalar spectator field, by allowing the effective mass m 2 (t) of this field and the cosmic time interval between non-adiabatic events to vary stochastically. Two main scenarios are considered depending on the (non-stochastic) mass M of the spectator field: the conformal case with M 2 =2H 2 , and the case of a massless field. We make use of the transfer matrix formalism to parametrize the evolution of the system in terms of the "occupation number", and two phases associated with the transfer matrix; these are used to construct the evolution of the spectator field. Assuming short-time interactions approximated by Dirac-delta functions, we numerically track the change of these parameters and the field in all regimes: sub- and super-horizon with weak and strong scattering. In all cases a log-normally distributed field amplitude is observed, and the logarithm of the field amplitude approximately satisfies the properties of a Wiener process outside the horizon. We derive a Fokker-Planck equation for the evolution of the transfer matrix parameters, which allows us to calculate analytically non-trivial distributions and moments in the weak-scattering limit.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Training models using forces computed by stochastic electronic structure methods

Abstract Quantum Monte Carlo (QMC) can play a very important role in generating accurate data needed for constructing potential energy surfaces. We argue that QMC has advantages in terms of a smaller systematic bias and an ability to cover phase space more completely. The stochastic noise can ease the training of the machine learning model. We discuss how stochastic errors affect the generation of effective models by analyzing the errors within a linear least squares procedure, finding that there is an advantage to having many relatively imprecise data points for constructing models. We then analyze the effect of noise on a model of many-body silicon finding that noise in some situations improves the resulting model. We then study the effect of QMC noise on two machine learning models of dense hydrogen used in a recent study of its phase diagram. The noise enables us to estimate the errors in the model. We conclude with a discussion of future research problems.

Ceperley, David M. (ORCID:0000000150826271)↗
DOE OSTI2020

Online accelerator optimization with a machine learning-based stochastic algorithm

Abstract Online optimization is critical for realizing the design performance of accelerators. Highly efficient stochastic optimization algorithms are needed for many online accelerator optimization problems in order to find the global optimum in the non-linear, coupled parameter space. In this study, we propose to use the multi-generation Gaussian process optimizer for online accelerator optimization and demonstrate that the algorithm is significantly more efficient than other stochastic algorithms that are commonly used in the accelerator community.

Zhang, Zhe (ORCID:0000000281430381)↗
DOE OSTI2024

Differentiable stochastic halo occupation distribution

ABSTRACT In this work, we demonstrate how differentiable stochastic sampling techniques developed in the context of deep reinforcement learning can be used to perform efficient parameter inference over stochastic, simulation-based, forward models. As a particular example, we focus on the problem of estimating parameters of halo occupation distribution (HOD) models that are used to connect galaxies with their dark matter haloes. Using a combination of continuous relaxation and gradient re-parametrization techniques, we can obtain well-defined gradients with respect to HOD parameters through discrete galaxy catalogue realizations. Having access to these gradients allows us to leverage efficient sampling schemes, such as Hamiltonian Monte Carlo, and greatly speed up parameter inference. We demonstrate our technique on a mock galaxy catalogue generated from the Bolshoi simulation using a standard HOD model and find near-identical posteriors as standard Markov chain Monte Carlo techniques with an increase of ∼8× in convergence efficiency. Our differentiable HOD model also has broad applications in full forward model approaches to cosmic structure and cosmological analysis.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

Connecting massive galaxies to dark matter haloes in BOSS – I. Is galaxy colour a stochastic process in high-mass haloes?

In this work, we use subhalo abundance matching (SHAM) to model the stellar mass function (SMF) and clustering of the Baryon Oscillation Spectroscopic Survey (BOSS) 'CMASS' sample at z ~ 0.5. We introduce a novel method which accounts for the stellar mass incompleteness of CMASS as a function of redshift, and produce CMASS mock catalogues which include selection effects, reproduce the overall SMF, the projected two-point correlation function w p , the CMASS dn/dz, and are made publicly available. We study the effects of assembly bias above collapse mass in the context of 'age matching' and show that these effects are markedly different compared to the ones explored by Hearin et al. at lower stellar masses. We construct two models, one in which galaxy colour is stochastic ('AbM' model) as well as a model which contains assembly bias effects ('AgM' model). By confronting the redshift dependent clustering of CMASS with the predictions from our model, we argue that that galaxy colours are not a stochastic process in high-mass haloes. Lastly, our results suggest that the colours of galaxies in high-mass haloes are determined by other halo properties besides halo peak velocity and that assembly bias effects play an important role in determining the clustering properties of this sample.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Stochastic Waveform Estimation at the Fundamental Quantum Limit

Although measuring the deterministic waveform of a weak classical force is a well-studied problem, estimating a random waveform, such as the spectral density of a stochastic signal field, is much less well understood despite it being a widespread task at the frontier of experimental physics. State-of-the-art precision sensors of random forces must account for the underlying quantum nature of the measurement but the optimal quantum protocol for interrogating such linear sensors is not known. We derive the fundamental precision limit: the extended-channel quantum Cramér-Rao bound. In the experimentally relevant regime in which losses dominate, we prove that non-Gaussian-state preparation and measurement are required to achieve this fundamental limit and we determine numerically the optimal non-Gaussian protocol. We discuss how this scheme could accelerate searches for signatures of quantum gravity, stochastic gravitational waves, and axionic dark matter.

Axions↗
DOE OSTI2024

Stochastic Maxwell-Bloch equations for modeling amplified spontaneous emission

An approach for studying atom-radiation interaction has been developed, associating quantum operators with stochastic variables governed by discrete Heisenberg equations. This framework models general multilevel atomic systems using the paraxial approximation in both single- and multipass configurations. Simulation results for x-ray laser oscillators with parameters of interest are presented. Finally, the stochastic model is efficient to solve, as the required computational resources scale linearly with the number of emitters, and it may be applicable to various other quantum systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Numerical modeling of a proof-of-principle experiment on optical stochastic cooling at an electron storage ring

Cooling of beams circulating in storage rings is critical for many applications including particle colliders and synchrotron light sources. A method enabling unprecedented beam-cooling rates, optical stochastic cooling (OSC), was recently demonstrated in the Integrable Optics Test Accelerator (IOTA) electron storage ring at Fermilab [J. Jarvis , ]. This paper describes the numerical implementation of the OSC process in the particle-tracking program and discusses the validation of the developed model with available experimental data. The model is also employed to highlight some features associated with different modes of operation of OSC. The developed simulation tool should be valuable in guiding future configurations of optical stochastic cooling and, more broadly, modeling self-field-based beam manipulations. Published by the American Physical Society 2024

43 PARTICLE ACCELERATORS↗
DOE OSTI2023

Topology, criticality, and dynamically generated qubits in a stochastic measurement-only Kitaev model

Here we consider a paradigmatic solvable model of topological order in two dimensions, Kitaev's honeycomb Hamiltonian, and turn it into a measurement-only dynamics consisting of stochastic measurements of two-qubit bond operators. We find an entanglement phase diagram that resembles that of the Hamiltonian problem in some ways, while being qualitatively different in others. When one type of bond is dominantly measured, we find area-law entangled phases that protect two topological qubits (on a torus) for a time exponential in system size. This generalizes the recently proposed idea of Floquet codes, where logical qubits are dynamically generated by a time-periodic measurement schedule, to a stochastic setting. When all types of bonds are measured with comparable frequency, we find a critical phase with a logarithmic violation of the area law, which sharply distinguishes it from its Hamiltonian counterpart. The critical phase has the same set of topological qubits, as diagnosed by the tripartite mutual information, but protects them only for a time polynomial in system size. Furthermore, we observe an unusual behavior for the dynamical purification of mixed states, characterized at late times by the dynamical exponent z = $\frac{1}{2}$ a superballistic dynamics made possible by measurements.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Two-nucleon 𝑆-wave interactions at the SU(3) flavor-symmetric point with 𝑚 𝑢⁢𝑑 ≈ 𝑚$^{phys}_𝑠$: A first lattice QCD calculation with the stochastic Laplacian Heaviside method

We report on the first application of the stochastic Laplacian Heaviside method for computing multiparticle interactions with lattice QCD to the two-nucleon system. Like the Laplacian Heaviside method, this method allows for the construction of interpolating operators which can be used to construct a set of positive-definite two-nucleon correlation functions, unlike nearly all other applications of lattice QCD to two nucleons in the literature. It also allows for a variational analysis in which optimal linear combinations of the interpolating operators are formed that couple predominantly to the eigenstates of the system. Utilizing such methods has become of paramount importance to help resolve the discrepancy in the literature on whether two nucleons in either isospin channel form a bound state at pion masses heavier than physical, with the discrepancy persisting even in the SU(3)-flavor-symmetric point with all quark masses near the physical strange quark mass. This is the first in a series of papers aimed at resolving this discrepancy. In the present work, we employ the stochastic Laplacian Heaviside method without a hexaquark operator in the basis at a lattice spacing of 𝑎 ≈0.086 fm, lattice volume of 𝐿 = 48⁢𝑎 ≈ 4.1 fm and pion mass 𝑚 𝜋 ≈ 714 MeV. With this setup, the observed spectrum of two-nucleon energy levels strongly disfavors the presence of a bound state in either the deuteron or dineutron channel.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Stochastic hydrodynamics and hydro-kinetics: Similarities and differences

The hydro-kinetic formalism has been used as a complementary approach to solving the Stochastic Differential Equations (SDE) corresponding to noisy hydrodynamics. The hydro-kinetic formalism consists of a deterministic set of relaxation type equations that tracks the evolution of 2-point correlation functions of stochastic hydrodynamic quantities. Hence they are comparatively easier to solve than the SDEs, which are computationally intensive and need to deal with arbitrarily large gradients. This work compares the two approaches for the propagation and diffusion of conserved charge fluctuations in the Bjorken hydrodynamic model. For white noise, the two approaches agree. For colored Catteneo noise, which is causal, the two approaches diverge. This is because white noise only induces two-point correlations, while Catteneo noise also induces higher-order correlations. Furthermore, this difference is quantified from the effects of causal evolution and influence from higher-order correlations induced by the Catteneo noise.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Multimessenger cosmology: Correlating cosmic microwave background and stochastic gravitational wave background measurements

Characterizing the physical properties of the stochastic gravitational wave background (SGWB) is a key step towards identifying the nature of its possible origin. We focus our analysis on SGWB anisotropies. The existence of a nontrivial primordial scalar-tensor-tensor (STT) correlation in the squeezed configuration may be inferred from the effect that a long wavelength scalar mode has on the gravitational wave power spectrum: an anisotropic contribution. Crucially, such a contribution is correlated with temperature anisotropies in the cosmic microwave background (CMB). This paper shows that, for inflationary models that generate suitably large STT non-Gaussianities, cross correlating the CMB with the stochastic background of gravitational waves is a very effective probe of early universe physics. The resulting signal can be a smoking gun for primordial SGWB anisotropies.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Gravitational waves from stochastic scalar fluctuations

Here, we present a novel mechanism for gravitational wave generation in the early Universe. Light spectator scalar fields during inflation can acquire a blue-tilted power spectrum due to stochastic effects. We show that this effect can lead to large curvature perturbations at small scales (induced by the spectator field fluctuations) while maintaining the observed, slightly red-tilted curvature perturbations at large cosmological scales (induced by the inflaton fluctuations). Along with other observational signatures, such as enhanced dark matter substructure, large curvature perturbations can induce a stochastic gravitational wave background (SGWB). The predicted strength of SGWB in our scenario, Ω GW ⁢ℎ 2 ≃10 -20 -10 -15 , can be observed with future detectors, operating between 10 -5 Hz and 10 Hz. We note that, in order to accommodate the newly reported NANOGrav observation, one could consider the same class of spectator models. At the same time, one would need to go beyond the simple benchmark considered here and consider a regime in which a misalignment contribution is also important.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Modeling electrokinetic flows with the discrete ion stochastic continuum overdamped solvent algorithm

In this article we develop an algorithm for the efficient simulation of electrolytes in the presence of physical boundaries. In previous work the discrete ion stochastic continuum overdamped solvent (DISCOS) algorithm was derived for triply periodic domains, and was validated through ion-ion pair correlation functions and Debye-Hückel-Onsager theory for conductivity, including the Wien effect for strong electric fields. In extending this approach to include an accurate treatment of physical boundaries we must address several important issues. First, the modifications to the spreading and interpolation operators necessary to incorporate interactions of the ions with the boundary are described. Next we discuss the modifications to the electrostatic solver to handle the influence of charges near either a fixed potential or dielectric boundary. An additional short-ranged potential is also introduced to represent interaction of the ions with a solid wall. Finally, the dry diffusion term is modified to account for the reduced mobility of ions near a boundary, which introduces an additional stochastic drift correction. Herein, several validation tests are presented confirming the correct equilibrium distribution of ions in a channel. Additionally, the methodology is demonstrated using electro-osmosis and induced-charge electro-osmosis, with comparison made to theory and other numerical methods. Notably, the DISCOS approach achieves greater accuracy than a continuum electrostatic simulation method. We also examine the effect of under-resolving hydrodynamic effects using a “dry diffusion” approach, and find that considerable computational speedup can be achieved with a negligible impact on accuracy.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Stochastic density functional theory combined with Langevin dynamics for warm dense matter

Here, this study overviews and extends a recently developed stochastic finite-temperature Kohn-Sham density functional theory to study warm dense matter using Langevin dynamics, specifically under periodic boundary conditions. The method's algorithmic complexity exhibits nearly linear scaling with system size and is inversely proportional to the temperature. Additionally, a linear-scaling stochastic approach is introduced to assess the Kubo-Greenwood conductivity, demonstrating exceptional stability for dc conductivity. Utilizing the developed tools, we investigate the equation of state, radial distribution, and electronic conductivity of hydrogen at a temperature of 30 000 K. As for the radial distribution functions, we reveal a transition of hydrogen from gaslike to liquidlike behavior as its density exceeds 4 g/cm 3 . As for the electronic conductivity as a function of the density, we identified a remarkable isosbestic point at frequencies around 7 eV, which may be an additional signature of a gas-liquid transition in hydrogen at 30 000 K.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗