Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “inference accelerators”

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 361 records · Page 20

MOOSE ProbML: Parallelizable Probabilistic Machine Learning and Uncertainty Quantification Capabilities

The Multiphysics Object Oriented Simulation Environment (MOOSE) is a widely used open- source finite element software for performing multiphysics multiscale simulations in a massively parallel fashion. Recently, the computational team at Idaho National Laboratory (INL) has implemented Probabilistic Machine Learning (ProbML) capabilities in MOOSE—in a parallelized fashion—and enable active learning with large-scale computational models for tasks such as surrogate model development, scale bridging, forward/inverse uncertainty quantification (UQ), Bayesian optimization, etc. This presentation summarizes these developments in MOOSE along with demonstrations on several real applications relevant to nuclear energy. At the fundamental level, samplers like Monte Carlo/Latin Hypercube, variance reduction, parallelized Markov Chain Monte Carlo (MCMC) support uncertainty propagation in both forward and inverse settings. These samplers can be integrated with the Gaussian processes (GP) suite in MOOSE, which offer several variants like scalar GPs, multi-output GPs, and deep GPs, to enable active learning. These GPs can be tuned using gradient-based optimization methods like Adam and its variants or gradient-free methods like the elliptical slice sampler (a variant of MCMC adept under Gaussian settings) for more complex covariance kernels or likelihoods whose gradient computations can be cumbersome. A variety of batch acquisition functions permit parallelized evaluation of the computational model and support different learning objectives with high efficiency like Bayesian inference, global surrogate development, optimization, etc. Furthermore, libtorch integration supports training, evaluation, and re-training of neural networks and other complex machine learning models in active learning settings. The impacts of these developments are shown on several real applications: (1) nuclear fuel inverse UQ and model inadequacy assessment using the Kennedy O’Hagan framework; (2) uncertainty aware surrogate modeling for additive manufacturing to predict field quantities; (3) nuclear reactor rare events analysis; and (4) complex fluid flow prediction using a global surrogate with quantified prediction uncertainty. Finally, the outlook of MOOSE ProbML is discussed for both outer-loop and inner-loop computations in the broad view to accelerate fuels and materials qualification, address gaps in knowledge and data, and assess new reactor/fuel systems.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Evaluation of data driven low-rank matrix factorization for accelerated solutions of the Vlasov equation

Low-rank methods have shown success in accelerating simulations of a collisionless plasma described by the Vlasov equation, but still rely on computationally costly linear algebra every time step. We propose a data-driven factorization method using artificial neural networks, specifically with convolutional layer architecture, that trains on existing simulation data. At inference time, the model outputs a low-rank decomposition of the distribution field of the charged particles, and we demonstrate that this step is faster than the standard linear algebra technique. Numerical experiments show that the method achieves comparable reconstruction accuracy for interpolation tasks, generalizing to unseen test data in a manner beyond just memorizing training data; patterns in factorization also inherently followed the same numerical trend as those within algebraic methods (e.g., truncated singular-value decomposition). However, when training on the first 70% of a time-series data and testing on the remaining 30%, the method fails to meaningfully extrapolate. Despite this limiting result, the technique may have benefits for simulations in a statistical steady-state or otherwise showing temporal stability. These results suggest that while the model offers a computationally efficient alternative for datasets with temporal stability, its current formulation is best suited for interpolation rather than for predicting future states in time-evolving systems. This study thus lays the groundwork for further refinement of neural network-based approaches to low-rank matrix factorization in high-dimensional plasma simulations.

97 MATHEMATICS AND COMPUTING↗

Study plasma interactions in the auroral ionosphere

Analyzed data from rocket flight, 29.007UE is presented. In a discrete electron arc the measured upward moving electrons are well accounted for by secondaries produced in collisional scattering of the measured downcoming electrons. No collective mechanisms need to invoke. The low energy downcoming electrons are accounted for by thermal plasma accelerated through a potential drop of a few kV that specularly reflects upward-moving lower energy electrons. No low altitude collective effects need to invoke in the arc. Simultaneous measurements of electric field by double probes on 29.007 and the Chatanika Radar allow one to infer that there are upward drifting ions above the discrete electron arc, and there is a westward neutral wind in the discrete arc. Two rocket payloads were built to investigate plasma effects in the pulsating aurora.

Anderson, H. R.↗

Spin-Up Instability of a Levitated Molten Drop in MHD-Flow Transition to Turbulence

When an alternating magnetic field interacts with induced eddy currents in a conducting body, there will be a repulsive force between the body and the driving coil system generating the field. This repulsive force is the basis of electromagnetic levitation, which allows containerless processing of different materials. The eddy currents in the conducting body also generate Joule heating. Axial rotation of electromagnetically levitated objects is a common observation in levitation systems and often an undesirable side effect of such experiments on 1-g and -g. There have been recent efforts to use magnetic damping and suppress this tendency of body rotation. The first report of rotation in EML drops was attributed to a slight asymmetry of the shape and location of the levitation coils could change the axis and speed of rotation. Other theories of sample rotation include a frequency difference in the traveling electromagnetic waves and a phase difference in two different applied fields of the same frequency. All of these different mechanisms share the following characteristics: the torque is small, constant for constant field strength, and very weakly dependent on the sample's temperature and phase (solid or liquid). During experiments on the MSL-1 (First Microgravity Science Laboratory) mission of the Space Shuttle (STS-83 and STS-94, April and July 1997), a droplet of palladium-silicon alloy was electromagnetically levitated for viscosity measurements. For the non-deforming droplet, the resultant MHD flow inside the drop is inferred from motion of impurities on the surface. These observations indicate formation of a pair of co-rotating toroidal flow structures inside the spheroidal levitated drop that undergo secondary flow instabilities. As rise in the fluid temperature rises, the viscosity falls and the internal flow accelerates and becomes oscillatory; and beyond a point in the experiments, the surface impurities exhibit non-coherent chaotic motion signifying emergence of turbulence inside the drop. In this work, a background of these set of observations will be given followed by a presentation of our results on the digital particle tracking analysis that has been performed on a number of available videos. The analysis indicates that the levitated drop attains a constant rotational speed during the melting phase and formation of the co-rotating axi-symmetric laminar toroidal structures. However, the rate of axial rotation increases dramatically during the deformation of the toroidal structures anti their breakup into chaotic entities. This new data suggests an interaction between the flow inside the levitated molten drop and the driving coils in the experiments. Possible mechanisms for this interaction will be reviewed. The data will also be used to make an assessment of existing theories on droplet rotation.

Abedian, B.↗

Comparative Study of Three Reconnection X Line Models at the Earth's Dayside Magnetopause Using in Situ Observations

This work examines the large-scale aspects of magnetic field reconnection at the Earth's dayside magnetopause. We use two sets of reconnection events, which are identified mostly by the in situ detection of accelerated and Alfvenic plasma flows. We intercompare three analytical models that predict the reconnection X line location and orientation, namely, the Trattner et al. (2007) and Swisdak and Drake (2007) models and also a modified version of the component merging model. In the first set of reconnection observations, we show three fortuitous, quasi-simultaneous dayside magnetopause crossing events where two widely separated spacecraft detect reconnection signatures, and the X line location and orientation can be inferred from the observations. We compare X line model predictions to those inferred from observations. These three reconnection events indicate the presence of an extended (greater than 7 Earth radii in length), component-type reconnection X line on Earth's dayside magnetopause connecting and structuring the reconnection signatures at locations far apart. In the second set of reconnection events, we analyze the X line models' performance in predicting the observed reconnection outflow direction, i.e., its north-south and/or east-west senses, in a total of 75 single, rather than multiple and quasi-simultaneous, magnetopause crossing events, where reconnection-associated plasma flows were clearly present. We found that the Swisdak and Drake's (2007) X line model performs slightly better, albeit not statistically significant, when predicting both accelerated plasma flow north-south and east-west components in 73% and 53% of the cases, respectively, as compared to the Trattner et al. (2007) model (70% north-south and 42% east-west) and the modified component merging model (66% north-south and 50% east-west).

Souza, V. M.↗

DeepLensSBI: Deep inference of simulated strong lenses in ground-based surveys

This code is used to train and test machine learning models and generate results and plots presented in 2501.08524 [astro-ph.IM]. The code is written in python. The goal of this work is to train ML models trained on simulated images of strong gravitational lenses. The trained model can then quickly infer properties of the lensed objects with uncertainty quantification.

Poh, Jason [Univ. of Chicago, IL (United States)] ↗

A Strong Gravitational Lens Is Worth a Thousand Dark Matter Halos: Inference on Small-scale Structure Using Sequential Methods

Strong gravitational lenses are a singular probe of the Universe’s small-scale structure—they are sensitive to the gravitational effects of low-mass (<10 10 M ⊙ ) halos even without a luminous counterpart. Recent strong-lensing analyses of dark matter structure rely on simulation-based inference (SBI). Modern SBI methods, which leverage neural networks as density estimators, have shown promise in extracting the halo-population signal. However, it is unclear whether the constraints from these models are limited by the methodology or the data. In this study, we introduce an accelerator-optimized simulation pipeline that can generate lens images with realistic subhalo populations in milliseconds. Leveraging this simulator, we identify the main limitation of our fiducial SBI analysis: training set size. We then adopt a sequential neural posterior estimation (SNPE) approach, allowing us to refine the training distribution to align with the observed data. Using only one-fifth as many mock Hubble Space Telescope images, SNPE matches the constraints on the low-mass halo population produced by our best nonsequential model. Our experiments suggest that an over 3 order-of-magnitude increase in training set size and GPU hours would be required to achieve an equivalent result without sequential methods. While the full potential of the existing lens sample remains to be explored, the notable improvement in constraining power enabled by our sequential approach highlights that current constraints are limited primarily by methodology and not the data itself. Moreover, our results emphasize the need to treat training set generation and model optimization as interconnected stages of any cosmological analysis using SBI.

79 ASTRONOMY AND ASTROPHYSICS↗

The plumes of Titan

The plasma density enhancements observed during the Voyager 1 Titan encounter are interpreted in terms of a single plume of plasma wrapped around Saturn by corotation. Consideration is given to the radial motions of the plume driven by fluctuations in solar wind pressure and the dispersal of the plasma by the centrifugal interchange instability, by heating, and by azimuthal acceleration. It is shown that Saturn cannot readily impose corotation directly on the plume and, incidentally, that the total dissipation associated with the Titan-magnetosphere interaction is insufficient to supply the power to produce the observed Titan ultraviolet airglow. A pickup velocity of 8 km/s is inferred on the basis of the observed velocity and a standing wave model. It is found that the Voyager plasma observations are consistent with the predictions of the model.

Eviatar, A.↗

Precipitating Electron Energy Flux and Characteristic Energies in Jupiter's Main Auroral Region as Measured by Juno/JEDI

The relationship between electron energy flux and the characteristic energy of electron distributions in the main auroral loss cone bridges the gap between predictions made by theory and measurements just recently available from Juno. For decades such relationships have been inferred from remote sensing observations of the Jovian aurora, primarily from the Hubble Space Telescope, and also more recently from Hisaki. However, to infer these quantities, remote sensing techniques had to assume properties of the Jovian atmospheric structure - leading to uncertainties in their profile. Juno's arrival and subsequent auroral passes have allowed us to obtain these relationships unambiguously for the first time, when the spacecraft passes through the auroral acceleration region. Using Juno /Jupiter Energetic particle Detector Instrument (JEDI), an energetic particle instrument, we present these relationships for the 30-kiloelectronvolts to 1-megaelectronvolts electron population. Observations presented here show that the electron energy flux in the loss cone is a nonlinear function of the characteristic or mean electron energy and supports both the predictions from Knight (1973, https://doi.org/10.1016/0032-0633(73)90093-7) and magnetohydrodynamic turbulence acceleration theories (e.g., Saur et al., 2003, https://doi.org/10.1029/2002GL015761). Finally, we compare the in situ analyses of Juno with remote Hisaki observations and use them to help constrain Jupiter's atmospheric profile. We find a possible solution that provides the best agreement between these data sets is an atmospheric profile that more efficiently transports the hydrocarbons to higher altitudes. If this is correct, it supports the previously published idea (e.g., Parkinson et al., 2006, https://doi.org/10.1029/2005JE002539) that precipitating electrons increase the hydrocarbon eddy diffusion coefficients in the auroral regions.

Auroral↗

Electron Scattering by High-Frequency Whistler Waves at Earth's Bow Shock

Electrons are accelerated to non-thermal energies at shocks in space and astrophysical environments. While different mechanisms of electron acceleration have been proposed, it remains unclear how non-thermal electrons are produced out of the thermal plasma pool. Here, we report in situ evidence of pitch-angle scattering of non-thermal electrons by whistler waves at Earths bow shock. On 2015 November 4, the Magnetospheric Multiscale (MMS) mission crossed the bow shock with an Alfvn Mach number is approximately 11 and a shock angle of approximately 84deg. In the ramp and overshoot regions, MMS revealed bursty enhancements of non-thermal (0.52 keV) electron flux, correlated with high-frequency (0.2 - 0.4 Omega(sub ce), where Omega(sub ce) is the cyclotron frequency) parallel-propagating whistler waves. The electron velocity distribution (measured at 30 ms cadence) showed an enhanced gradient of phase-space density at and around the region where the electron velocity component parallel to the magnetic field matched the resonant energy inferred from the wave frequency range. The flux of 0.5 keV electrons (measured at 1ms cadence) showed fluctuations with the same frequency. These features indicate that non-thermal electrons were pitch-angle scattered by cyclotron resonance with the high-frequency whistler waves. However, the precise role of the pitch-angle scattering by the higher-frequency whistler waves and possible nonlinear effects in the electron acceleration process remains unclear.

Oka, M.↗

A Nonparametric Method for the Inference of Halo Occupation Distributions

The galaxy–halo connection traces processes by which galaxies form and evolve. The halo occupation distribution (HOD) describes the relationship between galaxies and their host dark matter haloes. Measurements of the galaxy two-point correlation function (2PCF) allow us to extract information about the HODs of observed galaxy samples. Several parametric HOD models have been proposed in the literature, but the choice of parameterization restricts the space of possible HODs. To resolve this issue, we introduce a nonparametric HOD fitting method in which we train an emulator to learn the mappings among the galaxy 2PCF, physical properties used to select galaxy samples, and the HOD, all obtained from simulated past light cones constructed with the Santa Cruz semianalytic model. Implementing this emulator within a likelihood analysis framework, we derive constraints on the HOD of a galaxy sample when provided with a measurement of its 2PCF. Using the emulator to accelerate likelihood evaluations, we test the nonparametric HOD approach on a set of 2PCFs for mock galaxy samples drawn from the TNG100-1 simulation and selected above threshold values of stellar mass and star formation rate. Our framework is able to recover TNG100-1 HODs within 0.2 dex. We use the TNG100-1 mocks to tune the reported uncertainties to estimate those expected in the analysis of observations. Comparing to parametric HOD modelling routines applied to the same mock galaxy samples, our approach consistently infers the HOD with comparable or greater precision and accuracy.

Kennedy, Jacob [Rutgers Univ., Piscataway, NJ (Uni↗

Effect of bars on evolution of SDSS spiral galaxies

ABSTRACT We explore the significance of bars in triggering central star formation (SF) and active galactic nucleus (AGN) activity for spiral galaxy evolution using a volume-limited sample with 0.020 < z < 0.055, Mr < −19.5, and σ > 70 km s−1 selected from Sloan Digital Sky Survey Data Release 7. On a central SF rate–σ plane, we measure the fraction of galaxies with strong bars in our sample and also the AGN fractions for barred and non-barred galaxies, respectively. The comparison between the bar and AGN fractions reveals a causal connection between the two phenomena of SF quenching and AGN activity. A massive black hole and abundant gas fuels are sufficient conditions to trigger AGNs. We infer that the AGNs triggered by satisfying the two conditions drive the strong AGN feedback, suddenly suppressing the central SF and leaving the SF sequence. We find that in galaxies where either of the two conditions is not sufficient, bars are a great help for the AGN triggering, accelerating the entire process of evolution, which is particularly evident in pseudo-bulge galaxies. All of our findings are obtained only when plotted in terms of their central velocity dispersion and central SFR (not galactic scale SFR), indicating that the AGN-driven SF quenching is confined in the central kpc region.

Kim, Minbae↗

Materials Learning Algorithms (MALA): Scalable machine learning for electronic structure calculations in large-scale atomistic simulations

We present the Materials Learning Algorithms (MALA) package, a scalable machine learning framework designed to accelerate density functional theory (DFT) calculations suitable for large-scale atomistic simulations. Using local descriptors of the atomic environment, MALA models efficiently predict key electronic observables, including local density of states, electronic density, density of states, and total energy. The package integrates data sampling, model training and scalable inference into a unified library, while ensuring compatibility with standard DFT and molecular dynamics codes. We demonstrate MALA's capabilities with examples including boron clusters, aluminum across its solid-liquid phase boundary, and predicting the electronic structure of a stacking fault in a large beryllium slab. Scaling analyses reveal MALA's computational efficiency and identify bottlenecks for future optimization. With its ability to model electronic structures at scales far beyond standard DFT, MALA is well suited for modeling complex material systems, making it a versatile tool for advanced materials research.

Density functional theory↗

Measurement of U 238 ( n , n ' γ ) cross section data and their impact on reaction models

A better knowledge of (n,xn) reaction cross sections is important for both reaction modeling and energy applications. This article focuses on inelastic scattering of neutrons off 238 U for which improvements are needed to better constrain evaluations and solve inconsistencies in nuclear power reactor calculations. A new precise measurement of (n,xnγ) reaction cross sections on 238 U has been performed at the GELINA (Geel Electron LINear Accelerator) neutron facility operated by EC-JRC-Geel (Belgium) with the GRAPhEME (GeRmanium array for Actinides PrEcise MEasurements) setup. The prompt γ-ray spectroscopy method coupled to time-of-flight measurements is used to extract (n,xnγ) cross section values which can be further combined to infer the total neutron inelastic scattering cross section. Cross section data for 18 γ transitions (five never measured before) are presented and compared to the data in the literature. Here, emphasis is especially given to the uncertainty determination to produce partial cross section data as accurate as possible. Due to intrinsic limitations of the experimental method, the use of additional nuclear structure information coupled with theoretical modeling is required to determine the total (n, n') cross section over the whole neutron energy range. We have investigated modeling aspects of the 238 U(n,n'γ) cross sections related to the description of compound nucleus and preequilibirum mechanisms as well as the discrete part of nuclear structure. Through comparison between experimental and calculated (n,n'γ) cross sections, we pinpoint inaccuracies in the description of specific reaction mechanisms and challenge recently implemented models. This helps improving the whole modeling of the (n, n') reaction.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Ionospheric Flow and Escape of Ions from Titan and Venus

Titan and Venus are unmagnetized bodies that interact directly with the high speed plasmas flowing around them. The similarities of these interactions are used to reinforce the interpretations of measurements made at each body from different measurement sites. In particular, observations of plasma properties at Titan and Venus from Voyager I and Pioneer Venus, respectively, when considered together, tend to reinforce the concept that ions of ionospheric origin escape down the ionotails of each body. The plasma measurements at Titan were made in its ionotail, well above its ionosphere. They revealed plasma flowing from Titan and escaping down its ionotail. On the other hand, the measurements at Venus were made in its ionosphere, where ionospheric ions were inferred to be flowing upward toward Venus' ionotail. When these processes are applied to Titan's ionosphere, upward flow toward the ionotail is found to be possible, consistent with the plasma observed escaping further down the ionotail. Applying similar reasoning to Venus, the upward ionospheric flow is expected to accelerate and escape down its ionotail. The latter result is reinforced by the recent detection, from SOHO, of cold ions in the distant wake (at 1 AU), which were interpreted to originate in the ionosphere of Venus.

Hartle, R. E.↗

Hydra: computer vision for data quality monitoring

Hydra, initially developed for Hall-D in 2019, is a system that utilizes computer vision to perform near real time data quality monitoring. Since then, it has been deployed across all experimental halls at Jefferson Lab, with the CLAS12 collaboration in Hall-B being the first outside of GlueX to fully utilize Hydra. The system comprises back end processes that manage the models, their inferences, and the data flow. Finally, the front-end components, accessible via web pages, allow detector experts and shift crews to view and interact with the system.

47 OTHER INSTRUMENTATION↗

Consequences and inferences from tidal interactions in the solar system

A general theory of the rotation of solid bodies in the solar system leads to the following conclusions: No special initial conditions on the primordial rotation states of most tidally evolved bodies are necessary to account for their present rotation. The possible exception is the rotation of Venus, which must have been initially retrograde unless a liquid core-solid mantle interaction or an accelerating atmospheric tide could dominate the ordinary gravitational tides. Tides eventually drive all bodies in precessing orbits to 'Cassini states', where the spin vector, normal to the orbit and normal to the invariable plane remain coplanar as the first two precess about the latter. The constant obliquity in this configuration, together with the lowest order gravitational harmonics and the amplitude of the physical libration allows the determination of the central condensation and the extent of a liquid-core for Mercury. Criteria for the existence of observable amplitudes of the free wobble, free precession and libration of the moon are established which do not exclude current, non-zero values.

Peale, S. J.↗

Observation of Undulator Radiation Generated by a Single Electron Circulating in a Storage Ring and Possible Applications

An experimental study into the undulator radiation, generated by a single electron was carried out at the Integrable Optics Test Accelerator (IOTA) storage ring at Fermilab. The individual photons were detected by a Single Photon Avalanche Diode (SPAD) at an average rate of 1 detection per 300 revolutions in the ring. The detection events were continuously recorded by a picosecond event timer for as long as 1 minute at a time. The collected data were used to test if there is any deviation from the classically predicted Poissonian photostatistics. It was motivated by the observation * of sub-Poissonian statistics in a similar experiment. The observation * could be an instrumentation effect related to low detection efficiency and long detector dead time. In our experiment, the detector (SPAD) has a much higher efficiency (65%) and a much lower dead time. In addition, we show that the collected data (recorded detection times) can be used to study the synchrotron motion of a single electron and infer some parameters of the ring. For example, by comparing the results of simulation and measurement for the synchrotron motion we were able to estimate the magnitude of the RF phase jitter.

43 PARTICLE ACCELERATORS↗