Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Gaussian process”

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

Uniform Beam Simulation Technique for Beam Scans and Machine Learning Studies at Fermilab

Fermilab's neutrino facilities, including NuMI and the upcoming LBNF, use proton beams to produce positively and negatively charged pions and kaons. Detailed simulations are necessary to study particle interactions and beam propagation. To efficiently analyze beam scan effects, we propose a technique to generate multiple simulation samples with high statistics. These samples can be used to develop beamline simulation based machine learning applications. In this technique, we generate a uniformly distributed single simulation data sample. We calculate Gaussian weights for each beam configurations and apply them to post-processing measurements. In this poster, we demonstrate the proposed simulation technique. This technique reduces simulation time and computing resources significantly.

Wickremasinghe, Athula↗

Diffusion of energetic electrons in turbulent plasmas of the solar wind

A method of calculation of the diffusion coefficients D(v) of particles in velocity space, based on the statistical analysis of the motion of a great number of test electrons, is proposed. In the case of Langmuir turbulence developing in plasmas with fluctuating density inhomogeneities such as the solar wind, simulations provide coefficients D(v) which mainly depend on the Langmuir wave spectra and agree well with the analytical predictions D{sub th}(v) of the quasilinear theory of weak turbulence. Nevertheless, some noticeable differences exist with this theory: in the range of phase velocity of the short waves where the main part of the wave energy is concentrated, D(v) is noticeably smaller than D{sub th}(v), due to the scattering, the reflection, and the focusing effects encountered by the Langmuir waves when they interact with the plasma density inhomogeneities. Moreover, the probability of large velocity jumps in the particles’ trajectories essentially exceeds the probability of a Gaussian distribution. These large jumps, which are connected with the waves’ transformation processes, modify the nature of the particle diffusion, which is no more classical. These higher order effects cause the discrepancies observed with the quasilinear theory, which does not take them into account in its perturbative approach. The solar wind plasmas, which present fluctuating density inhomogeneities of noticeable average levels, are a very good laboratory to study such diffusion processes, which can eventually influence significantly on the development of essential physical phenomena, as electromagnetic radio emissions by type III solar radio bursts, for example.

79 ASTRONOMY AND ASTROPHYSICS↗

Roundoff noise and scaling in the digital implementation of control compensators

Researchers in digital signal processing have examined at length the effects of finite wordlength in the design of digital filters. The issues that have been considered apply to any digital system. In particular, the design of digital control systems must consider these issues. The paper will use, adapt, and extend the ideas developed in digital signal processing to the issue of roundoff noise in digital linear-quadratic-Gaussian (LQG) compensators. The roundoff noise effects for a particular LQG example and several different implementation structures will be then examined.

Moroney, P.↗

Constraints on variations in inflaton decay rate from modulated preheating

Modulated (p)reheating is thought to be an alternative mechanism for producing super-horizon curvature perturbations in CMB. But large non-gaussianity and iso-curvature perturbations produced by this mechanism rule out its acceptability as the sole process responsible for generating CMB perturbations. We explore the situation where CMB perturbations are mostly generated by usual quantum fluctuations of inflaton during inflation, but a modulated coupling constant between inflaton and a secondary scalar affects the preheating process and produces some extra curvature perturbations. If the modulating scalar field is considered to be a dark matter candidate, coupling constant between the fields has to be unnaturally fine tuned in order to keep the local-form non-gaussianity and the amplitude of iso-curvature perturbations within observational limit; otherwise parameters of the models have to be tightly constrained. Those constraints imply that the curvature perturbations generated by modulated preheating should be less than 15% of the total observed CMB perturbations. On the other hand if the modulating scalar field is not a dark matter candidate, parameters of the models could not be constrained, but the constraints on the maximum amount of the curvature perturbations coming from modulated preheating remain valid.

79 ASTRONOMY AND ASTROPHYSICS↗

Technical Note: Range verification system using edge detection method for a scintillator and a CCD camera system

Purpose: Three-dimensional irradiation with a scanned carbon-ion beam has been performed from 2011 at the authors’ facility. The authors have developed the rotating-gantry equipped with the scanning irradiation system. The number of combinations of beam properties to measure for the commissioning is more than 7200, i.e., 201 energy steps, 3 intensities, and 12 gantry angles. To compress the commissioning time, quick and simple range verification system is required. In this work, the authors develop a quick range verification system using scintillator and charge-coupled device (CCD) camera and estimate the accuracy of the range verification. Methods: A cylindrical plastic scintillator block and a CCD camera were installed on the black box. The optical spatial resolution of the system is 0.2 mm/pixel. The camera control system was connected and communicates with the measurement system that is part of the scanning system. The range was determined by image processing. Reference range for each energy beam was determined by a difference of Gaussian (DOG) method and the 80% of distal dose of the depth-dose distribution that were measured by a large parallel-plate ionization chamber. The authors compared a threshold method and a DOG method. Results: The authors found that the edge detection method (i.e., the DOG method) is best for the range detection. The accuracy of range detection using this system is within 0.2 mm, and the reproducibility of the same energy measurement is within 0.1 mm without setup error. Conclusions: The results of this study demonstrate that the authors’ range check system is capable of quick and easy range verification with sufficient accuracy.

60 APPLIED LIFE SCIENCES↗

Probing non-Gaussian stochastic gravitational wave backgrounds with LISA

The stochastic gravitational wave background (SGWB) contains a wealth of information on astrophysical and cosmological processes. A major challenge of upcoming years will be to extract the information contained in this background and to disentangle the contributions of different sources. In this paper we provide the formalism to extract, from the correlation of three signals in the Laser Interferometer Space Antenna (LISA), information about the tensor three-point function, which characterizes the non-Gaussian properties of the SGWB. This observable can be crucial to discriminate whether a SGWB has a primordial or astrophysical origin. Compared to the two-point function, the SGWB three-point function has a richer dependence on the gravitational wave momenta and chiralities. It can be used therefore as a powerful discriminator between different models. For the first time we provide the response functions of LISA to a general SGWB three-point function. As examples, we study in full detail the cases of an equilateral and squeezed SGWB bispectra, and provide the explicit form of the response functions, ready to be convoluted with any theoretical prediction of the bispectrum to obtain the observable signal. We further derive the optimal estimator to compute the signal-to-noise ratio. Our formalism covers general shapes of non-Gaussianity, and can be extended straightaway to other detector geometries. Finally, we provide a short overview of models of the early universe that can give rise to a non-Gaussian SGWB.

79 ASTRONOMY AND ASTROPHYSICS↗

A Closed-Loop Hardware Simulation of Decentralized Satellite Formation Control

In recent years, there has been significant interest in the use of formation flying spacecraft for a variety of earth and space science missions. Formation flying may provide smaller and cheaper satellites that, working together, have more capability than larger and more expensive satellites. Several decentralized architectures have been proposed for autonomous establishment and maintenance of satellite formations. In such architectures, each satellite cooperatively maintains the shape of the formation without a central supervisor, and processing only local measurement information. The Global Positioning System (GPS) sensors are ideally suited to provide such local position and velocity measurements to the individual satellites. An investigation of the feasibility of a decentralized approach to satellite formation flying was originally presented by Carpenter. He extended a decentralized linear-quadratic-Gaussian (LQG) framework proposed by Speyer in a fashion similar to an extended Kalman filter (EKE) which processed GPS position fix solutions. The new decentralized LQG architecture was demonstrated in a numerical simulation for a realistic scenario that is similar to missions that have been proposed by NASA and the U.S. Air Force. Another decentralized architecture was proposed by Park et al. using carrier differential-phase GPS (CDGPS). Recently, Busse et al demonstrated the decentralized CDGPS architecture in a hardware-in-the-loop simulation on the Formation Flying TestBed (FFTB) at Goddard Space Flight Center (GSFC), which features two Spirent Cox 16 channel GPS signal generator. Although representing a step forward by utilizing GPS signal simulators for a spacecraft formation flying simulation, only an open-loop performance, in which no maneuvers were executed based on the real-time state estimates, was considered. In this research, hardware experimentation has been extended to include closed-loop integrated guidance and navigation of multiple spacecraft formations using GPS receivers and real-time vehicle telemetry. A hardware closed-loop simulation has been performed using the decentralized LQG architecture proposed by Carpenter in the GPS test facility at the Center for Space Research (CSR). This is the first presentation using this type of hardware for demonstration of closed-loop spacecraft formation flying.

Ebimuma, Takuji↗

Using a coupled dispersion model to estimate depletion of a tritium oxide plume by a forest

Tritium processing facilities may release tritium oxide (HTO) to the atmosphere which poses potential health risks to exposed co-located workers and to offsite individuals. Most radiological consequence analyses determine HTO dose by applying Gaussian plume models to simulate the transport of HTO. Within these models, deposition velocity is used to assess the sum of all deposition processes acting on the plume. While this may account for vegetative and soil uptake or respiration processes, it may currently lack inclusion of the complex interactions within heterogeneous forested environments. In this complex morphology, dispersion patterns are significantly altered by changing flow regimes above and below the forest canopy and by the transfer of plume material across the canopy boundary. To determine the effects of a heterogeneous forest canopy on an airborne HTO plume, a Gaussian plume model coupled with an advection-diffusion plume model was applied to estimate transport in the free atmosphere above the forest and within the forest canopy and understory. During 2012, wind speed and wind direction measurements taken at 5 heights, ranging from 2-m to 28-m, on an instrumented meteorological tower located in a loblolly pine forest at the Department of Energy (DOE) Savannah River Site (SRS), near Aiken, SC. From these measurements, model predictions were made over a full spectrum of meteorological conditions. Deposition and resuspension velocities were calculated based on the model-predicted flux of plume material across the top of the forest canopy. Additionally, net deposition velocity of the plume material was calculated as the difference between the deposition and resuspension velocities. The 1st and 5th percentile net deposition velocities were estimated to be 0.7 cm s-1 and 1.2 cm s-1, respectively.

63 RADIATION, THERMAL, AND OTHER ENVIRON. POLLUTAN↗

Evolution of intermittent filaments in the scrape-off layer of NSTX

Filamentary structures naturally arise from background turbulence in the scrape-off layer (SOL) of plasmas, leading to significant particle and heat transport that can degrade overall plasma confinement. This enhanced transport can contribute to unacceptably high heat loads on plasma-facing components. As such, understanding the physics of SOL plasma filaments is critical for predicting and mitigating their effects in future fusion devices. On the National Spherical Torus Experiment (NSTX), plasma filaments—commonly referred to as “blobs”—were investigated using the gas-puff imaging (GPI) diagnostic in the edge and SOL regions. The analysis involved identifying, segmenting, and tracking the characteristic contours of the blobs in each frame of the GPI video sequences. Their evolution was characterized through shape descriptors, velocity, and angular velocity derived from their contour coordinates. The results indicate that as the blob area increases, their shapes become more concave and less circular, suggesting reduced structural stability in larger blobs. This result aligns with previous theoretical results where it was shown that larger blobs are more susceptible to instabilities [Krasheninnikov et al., J. Plasma Phys. 74, 679–717 (2008) and D'Ippolito et al., Phys. Plasmas 18, 060501 (2011)]. A positive correlation was observed between radial velocity and radial position, suggesting radially outward acceleration of the filaments, potentially driven by decreasing viscous drag toward the far SOL. Interestingly, blobs in background SOL turbulence exhibited minimal spinning in contrast to filaments originating from edge localized modes, which show substantial rotation during their paths [Lampert et al., Phys. Plasmas 29, 102502 (2022)]. Statistical analysis of the solidity and total curvature shape descriptors, along with their temporal evolution, revealed relatively broad, near-Gaussian distributions. This suggests that blob morphology is strongly influenced by stochastic turbulent processes in the surrounding plasma environment. Blob parameters were also compared with bulk plasma and radial profile measurements. Notable trends were found between blob rotation and poloidal velocity with collisionality and line-integrated density. These findings contribute to a deeper understanding of blob dynamics and provide valuable insights for refining SOL turbulence models.

Covariance and correlation↗

Tracking atomic structure evolution during directed electron beam induced Si-atom motion in graphene via deep machine learning

Using electron beam manipulation, we enable deterministic motion of individual Si atoms in graphene along predefined trajectories. Structural evolution during the dopant motion was explored, providing information on changes of the Si atom neighborhood during atomic motion and providing statistical information of possible defect configurations. The combination of a Gaussian mixture model and principal component analysis applied to the deep learning-processed experimental data allowed disentangling of the atomic distortions for two different graphene sublattices. This approach demonstrates the potential of e-beam manipulation to create defect libraries of multiple realizations of the same defect and explore the potential of symmetry breaking physics. The rapid image analytics enabled via a deep learning network further empowers instrumentation for e-beam controlled atom-by-atom fabrication. Here, the analysis described in the paper can be reproduced via an interactive Jupyter notebook at https://git.io/JJ3Bx.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Need for Caution in Interpreting Extreme Weather Statistics

Given the reality of anthropogenic global warming, it is tempting to seek an anthropogenic component in any recent change in the statistics of extreme weather. This paper cautions that such efforts may, however, lead to wrong conclusions if the distinctively skewed and heavy-tailed aspects of the probability distributions of daily weather anomalies are ignored or misrepresented. Departures of several standard deviations from the mean, although rare, are far more common in such a distinctively non-Gaussian world than they are in a Gaussian world. This further complicates the problem of detecting changes in tail probabilities from historical records of limited length and accuracy. A possible solution is to exploit the fact that the salient non-Gaussian features of the observed distributions are captured by so-called stochastically generated skewed (SGS) distributions that include Gaussian distributions as special cases. SGS distributions are associated with damped linear Markov processes perturbed by asymmetric stochastic noise and as such represent the simplest physically based prototypes of the observed distributions. The tails of SGS distributions can also be directly linked to generalized extreme value (GEV) and generalized Pareto (GP) distributions. The Markov process model can be used to provide rigorous confidence intervals and to investigate temporal persistence statistics. The procedure is illustrated for assessing changes in the observed distributions of daily wintertime indices of large-scale atmospheric variability in the North Atlantic and North Pacific sectors over the period 1872–2011. No significant changes in these indices are found from the first to the second half of the period.

54 ENVIRONMENTAL SCIENCES↗

Kaon Femtoscopy with Lévy-Stable Sources from $\sqrt{s_{NN}}$ = 200 GeV Au+Au Collisions at RHIC

Femtoscopy has the capacity to probe the space-time geometry of the particle-emitting source in heavy-ion collisions. In particular, femtoscopy of like-sign kaon pairs may shed light on the origin of non-Gaussianity of the spatial emission probability density. The momentum correlations between like-sign kaon pairs are measured in data recorded by the STAR experiment, from $\sqrt{s_{NN}}$ = 200 GeV Au+Au collisions at RHIC, BNL. Preliminary results hint at the possible existence of non-Gaussian, Lévy-stable sources, with the likely presence of an anomalous diffusion process in the signal for the identically charged kaon pairs so produced. More statistically significant studies at lower centre-of-mass energies may contribute to the search for the critical end point of QCD.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Control law design to meet constraints using SYNPAC-synthesis package for active controls

Major features of SYNPAC (Synthesis Package for Active Controls) are described. SYNPAC employs constrained optimization techniques which allow explicit inclusion of design criteria (constraints) in the control law design process. Interrelationships are indicated between this constrained optimization approach, classical and linear quadratic Gaussian design techniques. Results are presented that were obtained by applying SYNPAC to the design of a combined stability augmentation/gust load alleviation control law for the DAST ARW-2.

Adams, W. M., Jr.↗

Recent Results on "Approximations to Optimal Alarm Systems for Anomaly Detection"

An optimal alarm system and its approximations may use Kalman filtering for univariate linear dynamic systems driven by Gaussian noise to provide a layer of predictive capability. Predicted Kalman filter future process values and a fixed critical threshold can be used to construct a candidate level-crossing event over a predetermined prediction window. An optimal alarm system can be designed to elicit the fewest false alarms for a fixed detection probability in this particular scenario.

Martin, Rodney Alexander↗

Adaptive Key Component Control of Nonlinear Evolving Flexible Structures

An optimal alarm system and its approximations may use Kalman filtering for univariate linear dynamic systems driven by Gaussian noise to provide a layer of predictive capability. Predicted Kalman filter future process values and a fixed critical threshold can be used to construct a candidate level-crossing event over a predetermined prediction window. An optimal alarm system can be designed to elicit the fewest false alarms for a fixed detection probability in this particular scenario.

Balas, Mark↗

Initial conditions for slow-roll inflation in a random Gaussian landscape

In the landscape perspective, our Universe begins with a quantum tunneling from an eternally-inflating parent vacuum, followed by a period of slow-roll inflation. We investigate the tunneling process and calculate the probability distribution for the initial conditions and for the number of e-folds of slow-roll inflation, modeling the landscape by a small-field one-dimensional random Gaussian potential. We find that such a landscape is fully consistent with observations, but the probability for future detection of spatial curvature is rather low, P ∼ 10{sup −3}.

79 ASTRONOMY AND ASTROPHYSICS↗

Large Deviations for Level Sets of a Branching Brownian Motion and Gaussian Free Fields

We study deviation probabilities for the number of high positioned particles in branching Brownian motion and confirm a conjecture of Derrida and Shi. We also solve the corresponding problem for the two-dimensional discrete Gaussian free field. Our method relies on an elementary inequality for inhomogeneous Galton–Watson processes.

97 MATHEMATICS AND COMPUTING↗