Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “accelerator simulation”

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,729 records · Page 96
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

Ion acceleration in laser generated megatesla magnetic vortex

We report that Magnetic Vortex Acceleration (MVA) from near critical density targets is one of the promising schemes of laser-driven ion acceleration. 3D particle-in-cell simulations are used to explore a more extensive laser-target parameter space than previously reported in the literature as well as to study the laser pulse coupling to the target, the structure of the fields, and the properties of the accelerated ion beam in the MVA scheme. The efficiency of acceleration depends on the coupling of the laser energy to the self-generated channel in the target. The accelerated proton beams demonstrate a high level of collimation with achromatic angular divergence, and carry a significant amount of charge. For petawatt-class lasers, this acceleration regime provides a favorable scaling of the maximum ion energy with the laser power for the optimized interaction parameters. Lastly, the megatesla-level magnetic fields generated by the laser-driven coaxial plasma structure in the target are a prerequisite for accelerating protons to the energy of several hundred mega-electron-volts.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2017

Observation of Betatron X-Ray Radiation in a Self-Modulated Laser Wakefield Accelerator Driven with Picosecond Laser Pulses

We investigate a new regime for betatron x-ray emission that utilizes kilojoule-class picosecond lasers to drive wakes in plasmas. When such laser pulses with intensities of ~ 5 × 1 0 18 W / cm 2 are focused into plasmas with electron densities of ~ 1 × 1 0 19 cm - 3 , they undergo self-modulation and channeling, which accelerates electrons up to 200 MeV energies and causes those electrons to emit x rays. The measured x-ray spectra are fit with a synchrotron spectrum with a critical energy of 10–20 keV, and 2D particle-in-cell simulations were used to model the acceleration and radiation of the electrons in our experimental conditions

43 PARTICLE ACCELERATORS↗
DOE OSTI2017

Irradiation of materials with short, intense ion pulses at NDCX-II

Here, we present an overview of the performance of the Neutralized Drift Compression Experiment-II (NDCX-II) accelerator at Berkeley Lab, and report on recent target experiments on beam-driven melting and transmission ion energy loss measurements with nanosecond and millimeter-scale ion beam pulses and thin tin foils. Bunches with around 10 11 ions, 1 mm radius, and 2–30 ns full width at half maximum duration have been created with corresponding fluences in the range of 0.1–0.7 J/cm 2 . To achieve these short pulse durations and mm-scale focal spot radii, the 1.1 MeV [megaelectronvolt (10 6 eV)] He + ion beam is neutralized in a drift compression section, which removes the space charge defocusing effect during final compression and focusing. The beam space charge and drift compression techniques resemble necessary beam conditions and manipulations in heavy ion inertial fusion accelerators. In conclusion, quantitative comparison of detailed particle-in-cell simulations with the experiment plays an important role in optimizing accelerator performance

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

Laser-driven collisionless shock accelerated ion beams (Final Report)

The major goal of the project was to study laser-driven collisionless shock ion acceleration.This acceleration mechanism has been demonstrated using CO 2 lasers, but the current state of the art CO 2 lasers cannot access peak normalized vector potentials of much greater than $a_0$ ~ 2. Therefore the favorable scalings of ion energies with $a_0$ observed using numerical modeling cannot be tested experimentally with these systems. Therefore this project is to develop an experimental platform to investigate laser-driven shock ion acceleration using 1.053 µm lasers, which can access extremely relativistic conditions, $a_0$ > 20. A near-critical cluster jet target was developed for studying this regime and the target was characterized under a number of conditions. This included a cooled cluster regime where a knife edge was used to achieve front-side scale lengths of < 60 µm. The experimental Tcubed facility was used to perform experiments and successfully measured accelerated protons in a narrow divergence beam. Particle-in-cell simulations have been performed to show shocks capable of ion acceleration can be produced. In contrast to the "ideal'' target scalelength profile described by Fiuza et al., where the electrons are strongly heated in a localized region near peak density, we show electron heating within the underdense region is rapid enough to enable collisionless shock formation. This is potentially good news for future experiments since it relaxes the strict, and difficult to assemble target requirements.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

Implementation of a Particle Accelerator Beam Dynamics Code on Multi-Node GPUs

Particle accelerators play an important role in a wide range of scientific discoveries and industrial applications. The self-consistent multi-particle simulation based on the particle-in-cell (PIC) method has been used to study charged particle beam dynamics inside those accelerators. However, the PIC simulation is time-consuming and needs to use modern parallel computers for high-resolution applications. In this paper, we implemented a parallel beam dynamics PIC code on multi-node hybrid architecture computers with multiple Graphics Processing Units (GPUs). We used two methods to parallelize the PIC code on multiple GPUs and observed that the replication method is a better choice for moderate problem size and current computer hardware while the domain decomposition method might be a better choice for large problem size and more advanced computer hardware that allows direct communications among multiple GPUs. Using the multi-node hybrid architectures at Oak Ridge Leadership Computing Facility (OLCF), the optimized GPU PIC code achieves a reasonable parallel performance and scales up to 64 GPUs with 16 million particles.

43 PARTICLE ACCELERATORS↗
DOE OSTI2014

A Spin-Light Polarimeter for Multi-GeV Longitudinally Polarized Electron Beams

The physics program at the upgraded Jefferson Lab (JLab) and the physics program envisioned for the proposed electron-ion collider (EIC) include large efforts to search for interactions beyond the Standard Model (SM) using parity violation in electroweak interactions. These experiments need precision electron polarimetry with an uncertainty of <; 0.5 %. The spin dependent Synchrotron radiation, called "spin-light," can be used to monitor the electron beam polarization. Here, we develop a conceptual design for a "spin-light" polarimeter that can be used at a high intensity, multi-GeV electron accelerator. We have also built a Geant4 based simulation for a prototype device and report some of the results from these simulations

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Validation of edge turbulence codes against the TCV-X21 diverted L-mode reference case

Self-consistent full-size turbulent-transport simulations of the divertor and scrape-off-layer (SOL) of existing tokamaks have recently become feasible. This enables the direct comparison of turbulence simulations against experimental measurements. In this work, we perform a series of diverted ohmic L-mode discharges on the tokamak à configuration variable (TCV) tokamak, building a first-of-a-kind dataset for the validation of edge turbulence models. This dataset, referred to as TCV-X21, contains measurements from five diagnostic systems from the outboard midplane (OMP) to the divertor targets—giving a total of 45 one- and two-dimensional comparison observables in two toroidal magnetic field directions. The experimental dataset is used to validate three flux-driven 3D fluid-turbulence models—GBS, GRILLIX and TOKAM3X. With each model, we perform simulations of the TCV-X21 scenario, individually tuning the particle and power source rates to achieve a reasonable match of the upstream separatrix value of density and electron temperature. We find that the simulations match the experimental profiles for most observables at the OMP—both in terms of profile shape and absolute magnitude—while a comparatively poorer agreement is found towards the divertor targets. The match between simulation and experiment is seen to be sensitive to the value of the resistivity, the heat conductivities, the power injection rate and the choice of sheath boundary conditions. Additionally, despite targeting a sheath-limited regime, the discrepancy between simulations and experiment also suggests that the neutral dynamics should be included. The results of this validation show that turbulence models are able to perform simulations of existing devices and achieve reasonable agreement with experimental measurements. Where disagreement is found, the validation helps to identify how the models can be improved. By publicly releasing the experimental dataset and validation analysis, this work should help to guide and accelerate the development of predictive turbulence simulations of the edge and SOL.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2011

Efficient generation of fast ions from surface modulated nanostructure targets irradiated by high intensity short-pulse lasers

It's shown that the imposition of sub-laser wavelength relief structures on the surface of mass-limited-targets results into several folds higher short-pulse laser absorption, and consequently the efficient generation of fast ions. The optimum relief parameters for enhanced short-pulse laser absorption and higher ion acceleration are estimated numerically by particle-in-cell simulations and then corroborated by analytical scalings. The stability of the pre-imposed surface modulation during the laser pulse foil interaction is also examined.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2010

Gravity Modeling Effects on Surface-Interacting Vehicles in Supersonic Flight

A vehicle simulation is "surface-interacting" if the state of the vehicle (position, velocity, and acceleration) relative to the surface is important. Surface-interacting simulations per-form ascent, entry, descent, landing, surface travel, or atmospheric flight. The dynamics of surface-interacting simulations are influenced by the modeling of gravity. Gravity is the sum of gravitation and the centrifugal acceleration due to the world s rotation. Both components are functions of position relative to the world s center and that position for a given set of geodetic coordinates (latitude, longitude, and altitude) depends on the world model (world shape and dynamics). Thus, gravity fidelity depends on the fidelities of the gravitation model and the world model and on the interaction of these two models. A surface-interacting simulation cannot treat gravitation separately from the world model. This paper examines the actual performance of different pairs of world and gravitation models (or direct gravity models) on the travel of a supersonic aircraft in level flight under various start-ing conditions.

Madden, Michael M.↗
NASA NTRS2009

Further Investigations of Gravity Modeling on Surface-Interacting Vehicle Simulations

A vehicle simulation is "surface-interacting" if the state of the vehicle (position, velocity, and acceleration) relative to the surface is important. Surface-interacting simulations perform ascent, entry, descent, landing, surface travel, or atmospheric flight. The dynamics of surface-interacting simulations are influenced by the modeling of gravity. Gravity is the sum of gravitation and the centrifugal acceleration due to the world s rotation. Both components are functions of position relative to the world s center and that position for a given set of geodetic coordinates (latitude, longitude, and altitude) depends on the world model (world shape and dynamics). Thus, gravity fidelity depends on the fidelities of the gravitation model and the world model and on the interaction of the gravitation and world model. A surface-interacting simulation cannot treat the gravitation separately from the world model. This paper examines the actual performance of different pairs of world and gravitation models (or direct gravity models) on the travel of a subsonic civil transport in level flight under various starting conditions.

Madden, Michael M.↗
DOE OSTI2025

Massively parallel phase-field simulations targeting exascale

The interface thickness in the phase-field (PF) method limits its simulation scales. Consequently, large-scale PF simulations become prohibitively expensive for resolving the extremely fine microstructures that typically form during rapid solidification processing. This challenge is significant in predicting microstructure evolution in metal additive manufacturing and has been identified by the United States Department of Energy’s Exascale Computing Project. Here, to address this, we develop a multi-GPU and MPI-based massively parallel simulation code, utilizing state-of-the-art algorithms, software, and libraries, for large-scale three-dimensional (3D) PF simulations. We report the first GPU-parallel PF simulations on Frontier (currently the second TOP500 exascale cluster) and Summit machines, taking dendritic growth as an example problem. We evaluate the parallel performance of our implementation using scaling studies with more than 24 000 GPUs (among the largest known computations to date) and the acceleration performance using large-scale simulations of dendritic growth in 3D. Finally, massively parallel GPUs in these supercomputers enabled the first coupled multiscale simulations of laser melting and subsequent dendritic solidification on the scale of a full melt-pool, demonstrating the feasibility of performing PF simulations with a point total over 2 billion grid points within an acceptable time.

Exascale↗
DOE OSTI2023

FASDA: An FPGA-Aided, Scalable and Distributed Accelerator for Range-Limited Molecular Dynamics

Conducting long-term simulations of small molecules using Molecular Dynamics (MD) is crucial in drug design. However, traditional methods to accelerate the process including ASICs or GPUs, have limitations. ASICs are difficult to access, while publicly accessible GPU packages for MD have issues with more GPUs when processing small molecules, causing decreased performance. Conversely, FPGAs are renowned for their scalability and have been commonly used as accelerators for various applications, making them suitable for the task. Moreover, with the availability of FPGAs in the cloud, pharmaceutical developers can now access FPGA devices. The only missing piece is a design deployed for MD acceleration. Therefore, we present FASDA, the first FPGA-based MD accelerator available for community development. FASDA is designed to evaluate range-limited MD, which is the most resource-intensive and computationally demanding component in MD. It outperforms the state of-the-art GPU solution by 4.67x, significantly reducing the time for MD in drug design.

Wu, Chunshu↗
DOE OSTI2016

Revisiting the thermal effect on shock wave propagation in weakly ionized plasmas

Many researchers have investigated shock propagation in weakly ionized plasmas and observed the following anomalous effects: shock acceleration, shock recovery, shock weakening, shock spreading, and splitting. It was generally accepted that the thermal effect can explain most of the experimental results. However, little attention was paid to the shock recovery. In this paper, the shock wave propagation in weakly ionized plasmas is studied by fluid simulation. It is found that the shock acceleration, weakening, and splitting appear after it enters the plasma (thermal) region. The shock splits into two parts right after it leaves the thermal region. The distance between the splitted shocks keeps decreasing until they recover to one. This paper can explain a whole set of features of the shock wave propagation in weakly ionized plasmas. It is also found that both the shock curvature and the splitting present the same photoacoustic deflection (PAD) signals, so they cannot be distinguished by the PAD experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

Induced core formation time in subcritical magnetic clouds by large-scale trans-Alfvénic flows

We clarify the mechanism of accelerated core formation by large-scale nonlinear flows in subcritical magnetic clouds by finding a semi-analytical formula for the core formation time and describing the physical processes that lead to them. Recent numerical simulations show that nonlinear flows induce rapid ambipolar diffusion that leads to localized supercritical regions that can collapse. Here, we employ non-ideal magnetohydrodynamic simulations including ambipolar diffusion for gravitationally stratified sheets threaded by vertical magnetic fields. One of the horizontal dimensions is eliminated, resulting in a simpler two-dimensional simulation that can clarify the basic process of accelerated core formation. A parameter study of simulations shows that the core formation time is inversely proportional to the square of the flow speed when the flow speed is greater than the Alfvén speed. We find a semi-analytical formula that explains this numerical result. The formula also predicts that the core formation time is about three times shorter than that with no turbulence, when the turbulent speed is comparable to the Alfvén speed.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2012

Effect of target composition on proton acceleration in ultraintense laser-thin foil interaction

The interactions of ultraintense circularly polarized laser pulses with a mixed solid target and a double-layer target are studied by two-dimensional particle-in-cell simulations. Different carbon and proton compositions in the targets are used in the simulations. It is shown that the proton acceleration mechanisms in both targets are very sensitive to the ion density ratios between protons and carbon ions. For a mixed solid target, a relatively low proton density gives rise to monoenergetic peaks in the proton energy spectrum while a high proton density leads to a large cut-off energy and wide energy spread. With the increase of the ratio, the so-called directed-Coulomb-explosion becomes dominated over the radiation pressure. Surprisingly, for a double-layer target with a front proton layer and an ultrathin rear carbon layer, a highly monoenergetic proton beam with a peak energy of 1.7 GeV/u, an energy spread of {approx}4%, and a divergency angle of 2 Degree-Sign can be obtained, which might have diverse applications in medical therepy and proton imaging in future.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Verification and Validation of Spectral Element Code for Supercritical CO2 Flow in Vertical Heated Tubes

The investigation of heat transfer in supercritical CO2 (sCO2) has garnered considerable attention in recent decades, given sCO2's potential as a promising working fluid for advanced power conversion cycles. Despite previous research efforts, there are still gaps in our understanding of sCO2 heat transfer, particularly in conditions associated with heat transfer deterioration. To delve into sCO2 heat transfer more comprehensively, we propose employing the high-fidelity computational fluid dynamics code NekRS to simulate sCO2 flow using the large eddy simulation technique. Through graphics processing unit acceleration, NekRS achieves a higher computational speed than traditional CPU-based systems. However, before using NekRS in practical applications involving sCO2, it is imperative to perform verification and validation. Here, this paper presents our efforts to verify and validate the NekRS code's capability for simulating sCO2 using heated vertical tubes, where heat transfer deterioration usually happens. To accommodate the unique properties of sCO2, we have modified the NekRS code by integrating third-party property modules, such as REFPROP and PROPATH. Our simulations are compared with experimental and numerical data from the literature, instilling confidence in leveraging NekRS for future engineering applications. Our simulations also reveal that the accuracy of the property module significantly impacts the results, with REFPROP outperforming PROPATH for sCO2 properties. Additionally, we observed that, depending on the flow direction, buoyancy can either enhance or suppress turbulence in sCO2 flow. In upward flow, under certain conditions, the suppressed turbulence leads to heat transfer deterioration, resulting in elevated wall temperatures.

NekRS↗
DOE OSTI2024

Celeritas R&D Report: Accelerating Geant4

Celeritas is a new Monte Carlo (MC) detector simulation code designed for computationally intensive applications on high-performance heterogeneous architectures. In the past two years Celeritas has advanced from prototyping a Graphics Processing Unit (GPU)-based single physics model in infinite medium to implementing a full set of electromagnetic (EM) physics processes in complex geometries. The current release of Celeritas, version 0.4, has incorporated full device-based navigation, an event loop in the presence of magnetic fields, and detector hit scoring. New functionality incorporates a scheduler to offload electromagnetic physics to the GPU within a Geant4-driven simulation, enabling straightforward integration of Celeritas into the high energy physics (HEP) experimental frameworks CMSSW and ATLAS FullSimLight. On the Perlmutter supercomputer, Celeritas performs EM physics between 3× and 18× faster using the machine’s Nvidia GPUs compared to using only CPUs, corresponding to an electrical power efficiency up to a factor of 5. When running a multithreaded Geant4 ATLAS test beam application with full hadronic physics, using Celeritas to accelerate the EM physics results in an overall simulation speedup of 1.7–2.2× on GPU and 1.2× on CPU. In a CMS test application using tt¯ events and the prototype Run 4 configuration, compared to Geant4 CPU, Celeritas with a Nvidia A100 improves overall throughput up to a factor of 2.7× but cannot be efficiently shared with more than 8 cores.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2011

Sensitivity Analysis of the Off-Normal Conditions of the SPIDER Accelerator

In the context of the development of the 1 MV neutral beam injector for the ITER tokamak, the study on beam formation and acceleration has considerable importance. This effort includes the ion source and accelerator SPIDER (Source for Production of Ions of Deuterium Extracted from an Rf plasma) ion source, planned to be built in Padova, and designed to extract and accelerate a 355 A/m{sup 2} current of H{sup -}(or 285 A/m{sup 2} D{sup -}) up to 100 kV. Exhaustive simulations were already carried out during the accelerator optimization leading to the present design. However, as it is expected that the accelerator shall operate also in case of pre-programmed or undesired off-normal conditions, the investigation of a large set of off-normal scenarios is necessary. These analyses will also be useful for the evaluation of the real performances of the machine, and should help in interpreting experimental results, or in identifying dangerous operating conditions.The present contribution offers an overview of the results obtained during the investigation of these off-normal conditions, by means of different modeling tools and codes. The results, showed a good flexibility of the device in different operating conditions. Where the consequences of the abnormalities appeared to be problematic further analysis were addressed.

43 PARTICLE ACCELERATORS↗