Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “accelerator modeling”

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 145 records · Page 8

Accelerating Multiscale Materials Modeling with Machine Learning

The focus of this project is to accelerate and transform the workflow of multiscale materials modeling by developing an integrated toolchain seamlessly combining DFT, SNAP, LAMMPS, (shown in Figure 1-1) and a machine-learning (ML) model that will more efficiently extract information from a smaller set of first-principles calculations. Our ML model enables us to accelerate first-principles data generation by interpolating existing high fidelity data, and extend the simulation scale by extrapolating high fidelity data (10 2 atoms) to the mesoscale (10 4 atoms). It encodes the underlying physics of atomic interactions on the microscopic scale by adapting a variety of ML techniques such as deep neural networks (DNNs), and graph neural networks (GNNs). We developed a new surrogate model for density functional theory using deep neural networks. The developed ML surrogate is demonstrated in a workflow to generate accurate band energies, total energies, and density of the 298K and 933K Aluminum systems. Furthermore, the models can be used to predict the quantities of interest for systems with more number of atoms than the training data set. We have demonstrated that the ML model can be used to compute the quantities of interest for systems with 100,000 Al atoms. When compared with 2000 Al system the new surrogate model is as accurate as DFT, but three orders of magnitude faster. We also explored optimal experimental design techniques to choose the training data and novel Graph Neural Networks to train on smaller data sets. These are promising methods that need to be explored in the future.

36 MATERIALS SCIENCE↗

A GPU‐Accelerated Generative Adversarial Model for Causal Inference

We develop a GPU-accelerated machine learning generative adversarial model designed to facilitate causal inferences from observational data. Our model's theoretical framework is conceptualized in a manner that is amenable to being operable and scalable for high-performance computing platforms. We leverage GPU acceleration to develop a parallel evolutionary algorithm to achieve large-scale parallel computation of the model within a now widely accessible computing platform. This capability both enhances computational speedup and efficiency and also extends the use of the model to a broader range of substantive research domains while maintaining the underlying theoretical properties of the model.

GPU↗

Results and lessons learned from accelerating radio frequency modeling using machine learning [slides]

The “advanced tokamak” reactor concept is a leading candidate for a steady state fusion pilot plant. An advanced tokamak (AT) sustains a majority of the required plasma current with effects resulting from maintenance of the peaked pressure at the device center. This current is augmented by auxiliary current drive sources. These auxiliary actuators may consist of neutral particle beams and/or radio frequency (RF) systems such as lower hybrid current drive (LHCD) and high harmonic fast wave (HHFW) current drive using radio and microwaves from antennas. The primary focus of this work is to develop models of RF current profile control suitable for use in integrated modeling frameworks and for real-time control in experiments. Direct physics models of RF current drive can be computationally intensive. In order to achieve predictive times appropriate for the thousands of calls needed in real-time control of experiments and for use in integrated models, we will apply modern machine learning (ML) techniques to accelerate these models and interpolate their results. To generate the fast and accurate models for use in control level algorithms and integrated modeling we need to replace present models with high dimensional interpolation of their results. We will perform additional simulations across a broader parameter range for EAST and other tokamaks in different physics regimes (Alcator C-Mod, DIII-D, WEST, CFETR, ARC, ITER) and combine them into a larger database for training and testing of the ML models. Further testing of the control level models with experimental current profile data from EAST and C-Mod tokamaks will provide additional confirmation of the control level model before integration in a tokamak control system or integrated modeling suite. ML will be used to optimize the selection of training data consisting of RF current driven at different values of density profile, temperature profile, plasma current, and wavenumber. ML will also be used to facilitate classification of current drive from these input data. The output of this effort will be a validated classifier capable of determining the current drive profiles for HHFW CD and LHCD on a mille-second timescale. This will provide a breakthrough capability enabling real-time control of RF driven current profiles in experiments including ITER ICRF and use integrated modeling frameworks requiring thousands of current profile calculations in discharge simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Two-way coupled particle-turbulence interaction: Effect of numerics and resolution on fluid and particle statistics

Euler-Lagrange point-particle simulation has emerged as a premier methodology for studying dispersed particle-laden flows. This method's popularity stems from its ability to resolve fine-scale fluid structures while also tracking individual particles at reduced cost using an appropriate particle acceleration model. However, the point-particle model has known convergence issues in that refinement of the fluid grid can lead to changes in the predicted statistics. The reasons for non-convergence are twofold: the point-particle two-way coupling force in the Navier-Stokes equations requires a numerical regularization and, without careful implementation, yields a singular force on the fluid with grid refinement. The second factor that yields grid-dependent statistics is that the point-particle force model in general depends on the undisturbed fluid velocity. When the undisturbed fluid velocity is not robustly modeled in a grid-insensitive way, the calculated force for both particles and fluid will be grid-dependent, contaminating their respective statistics. While the first issue regarding regularizing the point-particle source term has received attention in the literature, the consequences of robustly modeling the undisturbed velocity in the context of grid refinement of turbulence has received little attention. Here, we consider decaying homogeneous isotropic turbulence laden with particles at different Stokes numbers. For a given Stokes number, we systematically refine the grid and demonstrate that explicitly modeling the undisturbed fluid velocity yields relative grid insensitivity for the energy of the particle and fluid phases, as well as acceleration of the particles. We also demonstrate that an appropriately defined dissipation rate is also grid-insensitive when an undisturbed fluid velocity correction is used. In contrast, when the undisturbed fluid velocity is modeled using the conventional approach of interpolating the local fluid velocity to the particle location, we show this procedure yields divergent statistics with grid refinement. In particular, we show that higher-order interpolation of the fluid velocity in two-way coupled problems is worse than lower-order interpolation, in the absence of a correction procedure to estimate the undisturbed fluid velocity. We also examine velocity derivative statistics of the fluid phase and demonstrate that these statistics are not in general convergent even when the undisturbed fluid velocity is explicitly modeled. Collectively, the observations in this paper are used to present a philosophy on the types of questions which are answerable with point-particle methods.

42 ENGINEERING↗

Testing Solar Flare Models with BATSE

We propose to use high-sensitivity Burst and Transient Source Experiment (BATSE) hard X-ray observations to test the thick-target and electric field acceleration models of solar flares. We will compare the predictions made by these models with hard X-ray spectral observations obtained with BATSE and simultaneous soft X-ray Ca XIX emission observed with the Yohkoh Bragg Crystal Spectrometer (BCS). The increased sensitivities of the BATSE and BCS (relative to previous detectors) permits a renewed study of the relationship between heating and dynamical motions during the crucial rise phase of flares. With these observations, we will: (1) investigate the ability of the thick-target model to explain the temporal evolution of hard X-ray emission relative to the soft X-ray blueshift during the earliest stages of the impulsive phase; and (2) search for evidence of electric-field acceleration as implied by temporal correlations between hard X-ray spectral breaks and the Ca XIX blueshift. The proposed study will utilize hard X-ray lightcurve and spectral measurements in the 10-100 keV energy range obtained with the BATSE Large Area Detectors (LAD). The DISCLA and CONT data will be the primary data products used in this analysis.

Zarro, Dominic M.↗

Physics through the 1990s: Gravitation, cosmology and cosmic-ray physics

The volume contains recommendations for space-and ground-based programs in gravitational physics, cosmology, and cosmic-ray physics. The section on gravitation examines current and planned experimental tests of general relativity; the theory behind, and search for, gravitational waves, including sensitive laser-interferometric tests and other observations; and advances in gravitation theory (for example, incorporating quantum effects). The section on cosmology deals with the big-bang model, the standard model from elementary-particle theory, the inflationary model of the Universe. Computational needs are presented for both gravitation and cosmology. Finally, cosmic-ray physics theory (nucleosynthesis, acceleration models, high-energy physics) and experiment (ground and spaceborne detectors) are discussed.

Source record↗

Modeling of Particle Acceleration at Multiple Shocks via Diffusive Shock Acceleration: Preliminary Results

Successful forecasting of energetic particle events in space weather models require algorithms for correctly predicting the spectrum of ions accelerated from a background population of charged particles. We present preliminary results from a model that diffusively accelerates particles at multiple shocks. Our basic approach is related to box models in which a distribution of particles is diffusively accelerated inside the box while simultaneously experiencing decompression through adiabatic expansion and losses from the convection and diffusion of particles outside the box. We adiabatically decompress the accelerated particle distribution between each shock by either the method explored in Melrose and Pope (1993) and Pope and Melrose (1994) or by the approach set forth in Zank et al. (2000) where we solve the transport equation by a method analogous to operator splitting. The second method incorporates the additional loss terms of convection and diffusion and allows for the use of a variable time between shocks. We use a maximum injection energy (E(sub max)) appropriate for quasi-parallel and quasi-perpendicular shocks and provide a preliminary application of the diffusive acceleration of particles by multiple shocks with frequencies appropriate for solar maximum (i.e., a non-Markovian process).

Parker, L. Neergaard↗

Acceleration processes in the quasi-steady magnetoplasmadynamic discharge

The flow field characteristics within the discharge chamber and exhaust of a quasi-steady magnetoplasmadynamic (MPD) arcjet were examined to clarify the nature of the plasma acceleration process. The observation of discharge characteristics unperturbed by insulator ablation and terminal voltage fluctuations, first requires the satisfaction of three criteria: the use of refractory insulator materials; a mass injection geometry tailored to provide propellant to both electrode regions of the discharge; and a cathode of sufficient surface area to permit nominal MPD arcjet operation for given combinations of arc current and total mass flow. The axial velocity profile and electromagnetic discharge structure were measured for an arcjet configuration which functions nominally at 15.3 kA and 6 g/sec argon mass flow. An empirical two-flow plasma acceleration model is advanced which delineates inner and outer flow regions and accounts for the observed velocity profile and calculated thrust of the accelerator.

Boyle, M. J.↗

Correlation between the He/H ratios in upstream particle events and in the solar wind

The ratio of helium to proton differential intensities at a fixed energy per charge in the diffuse ion population during long-lived upstream particle events is compared with the simultaneously measured solar wind He/H ratio. It is found that the He/H intensity ratio at about 30 keV in upstream events is highly correlated with and on the average directly proportional to the He/H ratio in the solar wind, and larger than the He/H ratio in the solar wind by an average factor of 1.6. It is concluded that the solar wind is the parent population for the diffuse ions. The He abundance in the diffuse component is on average observed to be larger than in the solar wind. The observations are consistent with the predictions of Fermi acceleration models for a seed particle population drawn from the solar wind with little compositional bias and moderately enhanced in He during the acceleration process.

Ipavich, F. M.↗

Comparison of Computed and Measured Performance of a Pulsed Inductive Thruster Operating on Argon Propellant

Pulsed inductive plasma accelerators are electrodeless space propulsion devices where a capacitor is charged to an initial voltage and then discharged through a coil as a high-current pulse that inductively couples energy into the propellant. The field produced by this pulse ionizes the propellant, producing a plasma near the face of the coil. Once a plasma is formed if can be accelerated and expelled at a high exhaust velocity by the Lorentz force arising from the interaction of an induced plasma current and the magnetic field. A recent review of the developmental history of planar-geometry pulsed inductive thrusters, where the coil take the shape of a flat spiral, can be found in Ref. [1]. Two concepts that have employed this geometry are the Pulsed Inductive Thruster (PIT)[2, 3] and the Faraday Accelerator with Radio-frequency Assisted Discharge (FARAD)[4]. There exists a 1-D pulsed inductive acceleration model that employs a set of circuit equations coupled to a one-dimensional momentum equation. The model was originally developed and used by Lovberg and Dailey[2, 3] and has since been nondimensionalized and used by Polzin et al.[5, 6] to define a set of scaling parameters and gain general insight into their effect on thruster performance. The circuit presented in Fig. 1 provides a description of the electrical coupling between the current flowing in the thruster I1 and the plasma current I2. Recently, the model was upgraded to include an equation governing the deposition of energy into various modes present in a pulsed inductive thruster system (acceleration, magnetic flux generation, resistive heating, etc.)[7]. An MHD description of the plasma energy density evolution was tailored to the thruster geometry by assuming only one-dimensional motion and averaging the plasma properties over the spatial dimensions of the current sheet to obtain an equation for the time-evolution of the total energy. The equation set governing the dynamics of the coupled electrodynamic-current sheet system is composed of first-order, coupled ordinary differential equations that can be easily solved numerically without having to resort to much more complex 2-D finite element plasma simulations.

Polzin, Kurt A.↗

REACTER: A Versatile Tool for Large-Scale Reactive Molecular Dynamics

Accurately describing reactive events over long length and time scales remains a grand challenge of computational materials science. REACTER is a general protocol for modeling chemical reactions using classical force fields, and is implemented in the popular molecular dynamics software LAMMPS. REACTER has a growing user base and has been used as a model-building tool for a variety of materials, including thermoplastics, thermosets, glassy materials and composites. The method has also been applied to accelerated modeling of reversible chemical reactions, such as the formation of electrochemical components for batteries. Recently, the REACTER protocol has received some major upgrades to enhance its ability to predict when reactions occur and to make it easier to use. Force field parameters can now be automatically assigned to newly created bonds, angles and other interactions. Advanced reaction constraints have also been added, including an Arrhenius constraint to enforce an effective activation energy, a root-mean-square-deviation option for complex geometrical constraints, and a constraint based on the total potential energy of the atoms involved in a reactive site. This potential energy constraint allows for the accurate reproduction of DFT-based tight-binding (DFTB3) predicted bond dissociation curves for mechanically induced bond breaking.

polymer simulations, molecular dynamics↗

Relativistic Electron Beams Accelerated By an Interplanetary Shock

Context: Collisionless shock waves have long been considered to be among the most prolific particle accelerators in the universe. Shocks alter the plasma they propagate through, and often exhibit complex evolution across multiple scales. Interplanetary (IP) traveling shocks have been recorded in situ for over half a century and act as a natural laboratory for experimentally verifying various aspects of large-scale collisionless shocks. A fundamentally interesting problem in both heliophysics and astrophysics is the acceleration of electrons to relativistic energies (> 300 keV) by traveling shocks. Aims: The reason for an incomplete understanding of electron acceleration at IP shocks is due to scale-related challenges and a lack of instrumental capabilities. This Letter presents the first observations of field-aligned beams of relativistic electrons upstream of an IP shock, observed thanks to the instrumental capabilities of Solar Orbiter. This study presents the characteristics of the electron beams close to the source and contributes to the understanding of their acceleration mechanism. Methods: On 25 July 2022, Solar Orbiter encountered an IP shock at 0.98 AU. The shock was associated with an energetic storm particle event, which also featured upstream field-aligned relativistic electron beams observed 14 min prior to the actual shock crossing. The distance of the beam’s origin was investigated using a velocity dispersion analysis (VDA). Peak-intensity energy spectra were analyzed and compared with those obtained from a semi-analytical fast-Fermi acceleration model. Results: By leveraging Solar Orbiter’s high temporal resolution Energetic Particle Detector (EPD), we successfully showcase an IP shock’s ability to accelerate relativistic electron beams. Our proposed acceleration mechanism offers an explanation for the observed electron beam and its characteristics, while we also explore the potential contributions of more complex mechanisms.

acceleration of particles↗

Addressing GPU memory limitations for Graph Neural Networks in High-Energy Physics applications

Introduction Reconstructing low-level particle tracks in neutrino physics can address some of the most fundamental questions about the universe. However, processing petabytes of raw data using deep learning techniques poses a challenging problem in the field of High Energy Physics (HEP). In the Exa.TrkX Project, an illustrative HEP application, preprocessed simulation data is fed into a state-of-art Graph Neural Network (GNN) model, accelerated by GPUs. However, limited GPU memory often leads to Out-of-Memory (OOM) exceptions during training, due to the large size of models and datasets. This problem is exacerbated when deploying models on High-Performance Computing (HPC) systems designed for large-scale applications. Methods We observe a high workload imbalance issue during GNN model training caused by the irregular sizes of input graph samples in HEP datasets, contributing to OOM exceptions. We aim to scale GNNs on HPC systems, by prioritizing workload balance in graph inputs while maintaining model accuracy. Our paper introduces diverse balancing strategies aimed at decreasing the maximum GPU memory footprint and avoiding the OOM exception, across various datasets. Results Our experiments showcase memory reduction of up to 32.14% compared to the baseline. We also demonstrate the proposed strategies can avoid OOM in application. Additionally, we create a distributed multi-GPU implementation using these samplers to demonstrate the scalability of these techniques on the HEP dataset. Discussion By assessing the performance of these strategies as data loading samplers across multiple datasets, we can gauge their effectiveness in both single-GPU and distributed environments. Our experiments, conducted on datasets of varying sizes and across multiple GPUs, broaden the applicability of our work to various GNN applications that handle input datasets with irregular graph sizes.

Lee, Claire Songhyun↗

Voyager energetic particle observations at interplanetary shocks and upstream of planetary bow shocks - 1977-1990

The Voyager 1 and 2 vehicles include instrumentation that makes comprehensive electron and ion measurements in several energy channels with good energy, temporal, and compositional resolution. Data gathered from 1977 to 1988, including observations downstream and upstream of four planetary bow shocks (earth, Saturn, Uranus, Jupiter) and numerous interplanetary shocks to about 30 AU, are analyzed in the context of the Fermi and shock drift acceleration models. Overall results indicate that electrons and ions observed upstream of planetary bow shocks have their source inside the parent magnetosphere, with first order Fermi acceleration playing a secondary role at best.

Krimigis, S. M.↗

Launch Alaska Transportation and Energy Accelerator (LATEA)

The Launch Alaska Transportation and Energy Accelerator (LATEA), funded through the U.S. Department of Energy Office of Technology Commercialization’s Energy Program for Innovation Clusters (EPIC),advanced deployment of innovative and efficient transportation and energy technology in Alaska from October 2021 through June 2025. The project was designed to leverage Launch Alaska’s accelerator model to identify, recruit, and support transportation technology companies with novel solutions to market needs while building the stakeholder networks, demonstration opportunities, and institutional capacity necessary to accelerate commercialization in one of the most challenging operating environments in the United States.

08 HYDROGEN↗

Hard Gamma Ray Emission from the Starburst Galaxy NGC 253

We have completed the study to search for hard gamma ray emission from the starburst galaxy NGC 253. Since supernovae are thought to provide the hard gamma ray emission from the Milky Way, starburst galaxies, with their extraordinarily high supernova rates, are prime targets to search for hard gamma ray emission. We conducted a careful search for hard gamma ray emission from NGC 253 using the archival data from the EGRET experiment aboard the CGRO. Because this starburst galaxy happens to lie near the South Galactic Pole, the Galactic gamma ray background is minimal. We found no significant hard gamma ray signal toward NGC 253, although a marginal signal of about 1.5 sigma was found. Because of the low Galactic background, we obtained a very sensitive upper limit to the emission of greater than 100 MeV gamma-rays of 8 x 10(exp -8) photons/sq cm s. Since we expected to detect hard gamma ray emission, we investigated the theory of gamma ray production in a dense molecular medium. We used a leaky-box model to simulate diffusive transport in a starburst region. Since starburst galaxies have high infrared radiation fields, we included the effects of self-Compton scattering, which are usually ignored. By modelling the expected gamma-ray and synchrotron spectra from NGC 253, we find that roughly 5 - 15% of the energy from supernovae is transferred to cosmic rays in the starburst. This result is consistent with supernova acceleration models, and is somewhat larger than the value derived for the Galaxy (3 - 10%). Our calculations match the EGRET and radio data very well with a supernova rate of 0.08/ yr, a magnetic field B approx. greater than 5 x 10(exp -5) G, a density n approx. less than 100/sq cm, a photon density U(sub ph) approx. 200 eV/sq cm, and an escape time scale tau(sub 0) approx. less than 10 Myr. The models also suggest that NGC 253 should be detectable with only a factor of 2 - 3 improvement in sensitivity. Our results are consistent with the standard picture of gamma-ray acceleration by supernovae.

Jackson, James M.↗

On solar cosmic rays

The report summarizes the major results obtained over the last four years in the following areas: (1) transport in interplanetary space which includes solar particle observations combined with theoretical efforts to develop models which simulate the dynamics of the event; (2) transport in the solar corona and release of particles in the interplanetary medium; and (3) characteristics of source spectra and acceleration models. Modelling the propagation of solar cosmic rays through interplanetary space and establishing the magnitude and the radial and energy dependence of the diffusion coefficient K are considered. An important result due to statistical observation is that the coronal propagation must be independent of or weakly dependent on both rigidity and energy. Increasing evidence is found for prolonged injection at the sun.

Van Hollebeke, M. A. I.↗

Resistive Wall Simulations for the DARHT Multi-Pulse Test Line

The multi-pulse test line (MPTL) presently under T development at the Dual Axis Radiography Hydrodynamic Test (DARHT) facility may have several meters of drifting beam transport without external magnetic focusing field. The absence of external magnetic focusing is the worst case for growth of the resistive-wall instability. This instability is usually considered to be strictly a long-pulse (low-frequency) problem, but has been theoretically shown to couple the multiple pulses in a long pulse train, such as are expected to be the subject of experimentation on the MPTL. In this note, we explore the resistive wall instability in the MPTL parameter range using our LAMDA beam dynamics code. (LAMDA stands for Linear Accelerator Model for DARHT)

43 PARTICLE ACCELERATORS↗