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,513 records · Page 84
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Parallel algorithms for hyperdynamics and local hyperdynamics

Hyperdynamics (HD) is a method for accelerating the timescale of standard molecular dynamics (MD). It can be used for simulations of systems with an energy potential landscape that is a collection of basins, separated by barriers, where transitions between basins are infrequent. HD enables the system to escape from a basin more quickly while enabling a statistically accurate renormalization of the simulation time, thus effectively boosting the timescale of the simulation. In [Kim, Perez, Voter, J Chem Phys, 139:144110, 2013)1, a local version of HD was formulated, which exploits the intrinsic locality characteristic typical of most systems to mitigate the poor scaling properties of standard HD as the system size is increased. In this paper, we discuss how both HD and local HD can be formulated to run efficiently in parallel. We have implemented these ideas in the LAMMPS MD code, which means HD can be used with any interatomic potential LAMMPS supports. Together, these parallel methods allow simulations of any size to achieve the time acceleration offered by HD (which can be orders of magnitude), at a cost 3-5x that of standard MD. As examples, we performed two simulations of a million-atom system to model the diffusion and clustering of Pt adatoms on a large patch of Pt(100) surface for 80 and 160 μs.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

EXPERIMENTAL STUDIES OF NONLINEAR INTEGRABLE OPTICS

State-of-the-art accelerators at energy and intensity frontiers require increasingly bright and powerful particle beams. In conventional linear lattices, intense beams suffer from collective instabilities, resulting in beam losses and maximum beam intensity limits. This thesis presents experimental studies of a novel lattice design concept, the nonlinear integrable optics (NIO), aimed at enhancing beam stability limits with little to no beam dynamics degradation. Single-particle beam dynamics measurements of two NIO devices, the quasi-integrable octupole system and the fully integrable Danilov-Nagaitsev system, were carried out at the purpose-built Integrable Optics Test Accelerator (IOTA) at Fermilab. Their simulation, hardware design, and commissioning process are presented. Extensive model and analysis algorithm development and benchmarking is described. Electron beam data from two scientific runs is analyzed, yielding frequency and phase space dynamics consisten t with m odels. These results demonstrate viability and advantages of the NIO design, providing the groundwork for proton studies in the strong space-charge regime and future integrable accelerators.

Kuklev, Nikita↗
NASA NTRS1992

Precipitation of fast ion beams from the plasma sheet boundary layer

This paper presents a model of precipitated fluxes from the PSBL and CPS. Simulation results and data from Aureol-3 spacecraft indicate the presence of velocity dispersed precipitated ion structures (VDIS) at the poleward edge of the auroral oval. These structures are associated with fast ion beams in the PSBL region of the earth's magnetotail, confirming previous experimental results. The simulations also reveal possible substructuring of the VDIS. The bulk of the PSBL population which is not precipitated is very effectively thermalized and quasi-isotropized after multiple interactions with the magnetotail current layer. After each reflection cycle some part of the distribution is precipitated and forms multiple 'echoes' of VDIS. The CPS distributions occurring as a result of scattering, convection, multiple reflections and Fermi acceleration appear isotropic in the simulation model. This paper portrays the important role of the VDIS auroral region medium for complicated and energetically significant processes occurring in different regions of the distant magnetotail.

Ashour-Abdalla, M.↗
NASA NTRS2005

Dusty-Plasma Particle Accelerator

A dusty-plasma apparatus is being investigated as means of accelerating nanometer- and micrometer-sized particles. Applications for the dusty-plasma particle accelerators fall into two classes: Simulation of a variety of rapidly moving dust particles and micrometeoroids in outer-space environments that include micrometeoroid streams, comet tails, planetary rings, and nebulae and Deposition or implantation of nanoparticles on substrates for diverse industrial purposes that could include hardening, increasing thermal insulation, altering optical properties, and/or increasing permittivities of substrate materials. Relative to prior apparatuses used for similar applications, dusty-plasma particle accelerators offer such potential advantages as smaller size, lower cost, less complexity, and increased particle flux densities. A dusty-plasma particle accelerator exploits the fact that an isolated particle immersed in plasma acquires a net electric charge that depends on the relative mobilities of electrons and ions. Typically, a particle that is immersed in a low-temperature, partially ionized gas, wherein the average kinetic energy of electrons exceeds that of ions, causes the particle to become negatively charged. The particle can then be accelerated by applying an appropriate electric field. A dusty-plasma particle accelerator (see figure) includes a plasma source such as a radio-frequency induction discharge apparatus containing (1) a shallow cup with a biasable electrode to hold the particles to be accelerated and (2) a holder for the substrate on which the particles are to impinge. Depending on the specific design, a pair of electrostatic-acceleration grids between the substrate and discharge plasma can be used to both collimate and further accelerate particles exiting the particle holder. Once exposed to the discharge plasma, the particles in the cup quickly acquire a negative charge. Application of a negative voltage pulse to the biasable electrode results in the initiation of a low-current, high-voltage cathode spot. Plasma pressure associated with the cathode spot as well as the large voltage drop at the cathode spot accelerates the charged particles toward the substrate. The ultimate kinetic energy attained by particles exiting the particle holder depends in part on the magnitude of the cathode spot sheath potential difference, which is proportional to the magnitude of the voltage pulse, and the on the electric charge on the dust. The magnitude of the voltage pulse can be controlled directly, whereas the particle s electric charge can be controlled indirectly by controlling the operating parameters of the plasma apparatus.

Foster, John E.↗
DOE OSTI2014

Ion probe beam experiments and kinetic modeling in a dense plasma focus Z-pinch

The Z-pinch phase of a dense plasma focus (DPF) emits multiple-MeV ions in a ∼cm length. The mechanisms through which these physically simple devices generate such high energy beams in a relatively short distance are not fully understood. We are exploring the origins of these large gradients using measurements of an ion probe beam injected into a DPF during the pinch phase and the first kinetic simulations of a DPF Z-pinch. To probe the accelerating fields in our table top experiment, we inject a 4 MeV deuteron beam along the z-axis and then sample the beam energy distribution after it passes through the pinch region. Using this technique, we have directly measured for the first time the acceleration of an injected ion beam. Our particle-in-cell simulations have been benchmarked on both a kJ-scale DPF and a MJ-scale DPF. They have reproduced experimentally measured neutron yields as well as ion beams and EM oscillations which fluid simulations do not exhibit. Direct comparisons between the experiment and simulations enhance our understanding of these plasmas and provide predictive design capability for accelerator and neutron source applications.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2002

Hubble Space Telescope Pointing Performance Due to Micro-Dynamic Disturbances from the NICMOS Cryogenic Cooler

The Vibration Admittance Test (VET) was performed to measure the emitted disturbances of the Near Infrared Camera and Multi-Object Spectrometer (NICMOS) Cryogenic Cooler (NCC) in preparation for NCC installation onboard the Hubble Space Telescope (HST) during Servicing Mission 3B (SM3B). Details of the VET ground-test are described, including facility characteristics, sensor complement and configuration, NCC suspension, and background noise measurements. Kinematic equations used to compute NCC mass center displacements and accelerations from raw measurements are presented, and dynamic equations of motion for the NCC VET system are developed and verified using modal test data. A MIMO linear frequency-domain analysis method is used to compute NCC-induced loads and HST boresight jitter from VET measurements. These results are verified by a nonlinear time-domain analysis approach using a high-fidelity structural dynamics and pointing control simulation for HST. NCC emitted acceleration levels not exceeding 35 micro-g rms were measured in the VET and analysis methods herein predict 3.1 milli-areseconds rms jitter for HST on-orbit. Because the NCC is predicted to become the predominant disturbance source for HST, VET results indicate that HST will continue to meet the 7 milli-arcsecond pointing stability mission requirement in the post-SM3B era.

Clapp, Brian R.↗
DOE OSTI2012

Plasma Power Station with Quasi Spherical Direct Drive Capsule for Fusion Yield and Inverse Diode for Driver-Target Coupling

The Meier-Moir economic model for Pulsed Power Driven Inertial Fusion Energy shows at least two approaches for fusion energy at 7 to 8 cents/kw-hr: One with large yield at 0.1 Hz and presented by M. E. Cuneo at ICENES 2011 and one with smaller yield at 3 Hz presented in this paper. Both use very efficient and low cost Linear Transformer Drivers (LTDs) for the pulsed power. Here, we report the system configuration and end-to-end simulation for the latter option, which is called the Plasma Power Station (PPS), and report the first results on the two, least mature, enabling technologies: a magnetically driven Quasi Spherical Direct Drive (QSDD) capsule for the fusion yield and an Inverse Diode for coupling the driver to the target. In addition, we describe the issues and propose to address the issues with a prototype of the PPS on the Saturn accelerator and with experiments on a short pulse modification of the Z accelerator test the validity of simulations showing megajoule thermonuclear yield with DT on a modified Z.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2024

RADIATION DOSE SIMULATIONS FOR JEFFERSON LAB’S PERMANENT MAGNET RESILIENCY LDRD STUDY

In late 2023, Thomas Jefferson National Accelerator Facil- ity (Jefferson Lab) funded a Laboratory Directed Research and Development (LDRD) grant dedicated to investigating the impact of radiation on permanent magnet materials. This research initiative is specifically geared towards assessing materials slated for use in the CEBAF energy upgrade. The experimental approach involves strategically placing perma- nent magnet samples throughout the accelerator, exposing them to varying radiation doses. The simulation code BD- SIM is used to first validate the data and then to simulate the effects on future higher energy passes to study the degra- dation effects on the permanent magnets. In this paper we present the progress of that work.

Deitrick, K.↗
DOE OSTI2024

Accelerating transients with NekRS: GPU overlapping domain implementation and multi-rate timestepping

The simulation of nuclear transients using Computational Fluid Dynamics (CFD) presents significant computational challenges due to the inherent complexity and the wide separation in temporal scales between various flow physical phenomena. These disparities lead to high computational costs, often making the simulation of transients impractical without advanced techniques. Consequently, multiple research initiatives are being pursued by the NEAMS thermal-hydraulic area, some driven by academic institutions and some by national laboratories. Overall, they are exploring novel methods to make transient simulations more feasible and efficient. This report delves into recent advancements within the CFD code NekRS, specifically those achieved in Fiscal Year 2024 under the CONNECT effort, aimed at improving the performance and feasibility of transient simulations. The first major advancement involves the porting of NekRS to Aurora, one of the Department of Energy’s (DOE) most powerful supercomputers. Additionally, the report discusses the implementation of an overlapping domain capability within NekRS. This novel GPU-accelerated capability allows different spatial regions of the domain to be solved independently, enhancing the code’s efficiency, particularly when running large-scale simulations in complex domains. The scalability of this approach is demonstrated, highlighting its potential to transform how transients are approached in CFD simulations. Lastly, the report focuses on how this overlapping domain capability specifically accelerates transient simulations through multi-rate timestepping. By decoupling different regions and facilitating faster computations, this method offers a promising pathway to making nuclear transient simulations more computationally feasible, addressing one of the critical bottlenecks in the field. Together, these advancements represent a significant leap forward in transient simulation technology, bringing closer the possibility of handling highly complex nuclear scenarios with greater efficiency and accuracy.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2016

Fermilab Booster Transition Crossing Simulations and Beam Studies

The Fermilab Booster accelerates beam from 400 MeV to 8 GeV at 15 Hz. In the PIP (Proton Improvement Plan) era, it is required that Booster deliver 4.2 x $10^{12}$ protons per pulse to extraction. One of the obstacles for providing quality beam to the users is the longitudinal quadrupole oscillation that the beam suffers from right after transition. Although this oscillation is well taken care of with quadrupole dampers, it is important to understand the source of these oscillations in light of the PIP II requirements that require 6.5 x $10^{12}$ protons per pulse at extraction. This paper explores the results from machine studies, computer simulations and solutions to prevent the quadrupole oscillations after transition.

Bhat, C. M.↗
NASA NTRS2009

Validation of Finite Element Crash Test Dummy Models for Predicting Orion Crew Member Injuries During a Simulated Vehicle Landing

A series of crash tests were conducted with dummies during simulated Orion crew module landings at the Wright-Patterson Air Force Base. These tests consisted of several crew configurations with and without astronaut suits. Some test results were collected and are presented. In addition, finite element models of the tests were developed and are presented. The finite element models were validated using the experimental data, and the test responses were compared with the computed results. Occupant crash data, such as forces, moments, and accelerations, were collected from the simulations and compared with injury criteria to assess occupant survivability and injury. Some of the injury criteria published in the literature is summarized for completeness. These criteria were used to determine potential injury during crew impact events.

Tabiei, Al↗
DOE OSTI2015

Technical Note: Motion-perturbation method applied to dosimetry of dynamic MLC target tracking—A proof-of-concept

Purpose: Previous studies show that dose to a moving target can be estimated using 4D measurement-guided dose reconstruction based on a process called virtual motion simulation, or VMS. A potential extension of VMS is to estimate dose during dynamic multileaf collimator (MLC)-tracking treatments. The authors introduce a modified VMS method and quantify its performance as proof-of-concept for tracking applications. Methods: Direct measurements with a moving biplanar diode array were used to verify accuracy of the VMS dose estimates. A tracking environment for variably sized circular MLC apertures was simulated by sending preprogrammed control points to the MLC while simultaneously moving the accelerator treatment table. Sensitivity of the method to simulated tracking latency (0–700 ms) was also studied. Potential applicability of VMS to fast changing beam apertures was evaluated by modeling, based on the demonstrated dependence of the cumulative dose on the temporal dose gradient. Results: When physical and virtual latencies were matched, the agreement rates (2% global/2 mm gamma) between the VMS and the biplanar dosimeter were above 96%. When compared to their own reference dose (0 induced latency), the agreement rates for VMS and biplanar array track closely up to 200 ms of induced latency with 10% low-dose cutoff threshold and 300 ms with 50% cutoff. Time-resolved measurements suggest that even in the modulated beams, the error in the cumulative dose introduced by the 200 ms VMS time resolution is not likely to exceed 0.5%. Conclusions: Based on current results and prior benchmarks of VMS accuracy, the authors postulate that this approach should be applicable to any MLC-tracking treatments where leaf speeds do not exceed those of the current Varian accelerators.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2012

FEA Simulations of Magnets with Grain Oriented Steel

One of the potential successors of the Large Hadron Collider is a Muon Collider. Muons are short-lived particles, which therefore require fast acceleration. One potential avenue is a very fast cycling cyclotron, where the bending is sup- plied by a combination of fixed-field superconducting magnets and fast ramping normal conducting iron-cored coils. Due to the high ramping rate (around 1 kHz) eddy current and hysteresis losses are a concern. One way to overcome these is by using grain-oriented soft-iron, which promises superior magnetic properties in the direction of the grains. This note summarizes efforts to include the anisotropic material properties of grain-oriented steel in finite element analysis to predict the behaviour of the dipole magnets for this accelerator. It was found that including anisotropic material properties has a detrimental effect on model convergence. During this study it was not possible to include grain oriented steel with an accuracy necessary to study the field quality of a dipole magnet.

43 PARTICLE ACCELERATORS↗
DOE OSTI2012

Dosimetric Impact of Breathing Motion in Lung Stereotactic Body Radiotherapy Treatment Using Image-Modulated Radiotherapy and Volumetric Modulated Arc Therapy

Purpose: The objective of this study was to investigate the influence of tumor motion on dose delivery in stereotactic body radiotherapy (SBRT) for lung cancer, using fixed field intensity- modulated radiotherapy (IMRT) and volumetric modulated arc therapy (VMAT). Methods and Materials: For each of 10 patients with stage I/II non-small-cell pulmonary tumors, a respiration-correlated four-dimensional computed tomography (4DCT) study was carried out. The internal target volume was delineated on the maximum intensity projection CT, which was reconstructed from the 4DCT dataset. A 5-mm margin was used for generation of the planning target volume. VMAT and five-field IMRT plans were generated using Pinnacle{sup 3} SmartArc and direct machine parameter optimization, respectively. All plans were generated for an Elekta Synergy linear accelerator using 6-MV photons. Simulation was performed to study the interplay between multileaf collimator (MLC) sequences and target movement during the delivery of VMAT and IMRT. For each plan, 4D dose was calculated using deformable image registration of the 4DCT images. Target volume coverage and doses to critical structures calculated using 4D methodology were compared with those calculated using 3D methodology. Results: For all patients included in this study, the interplay effect was found to present limited impact (less than 1% of prescription) on the target dose distribution, especially for SBRT, in which fewer fractions (three fractions) are delivered. Dose to the gross tumor volume (GTV) was, on average, slightly decreased (1% of prescription) in the 4D calculation compared with the 3D calculation. The motion impact on target dose homogeneity was patient-dependent and relatively small. Conclusions: Both VMAT and IMRT plans experienced negligible interplay effects between MLC sequence and tumor motion. For the most part, the 3D doses to the GTV and critical structures provided good approximations of the 4D dose calculations.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2018

Beam-Dynamics Analysis of Long-Range Wakefield Effects on the SCRF Cavities at the Fast Facility

Long-range wakefields in superconducting RF (SCRF) cavities create complicated effects on beam dynamics in SCRF-based FEL beamlines. The driving bunch excites effectively an infinite number of structure modes (including HOMs) which oscillate within the SCRF cavity. Couplers with loads are used to damp the HOMs. However, these HOMs can persist for long periods of time in superconducting structures, which leads to long-range wakefields. Clear understanding of the long-range wakefield effects is a critical element for risk mitigation of future SCRF accelerators such as XFEL at DESY, LCLS-II XFEL, and MaRIE XFEL. We are currently developing numerical tools for simulating long-range wakefields in SCRF accelerators and plan to experimentally verify the tools by measuring these wakefields at the Fermilab Accelerator Science and Technology (FAST) facility. This paper previews the experimental conditions at the FAST 50 MeV beamline based on the simulation results.

43 PARTICLE ACCELERATORS↗
DOE OSTI2012

LINAC BEAM DYNAMICS SIMULATIONS WITH PY-ORBIT

Linac dynamics simulation capabilities of the PyORBIT code are discussed. PyORBIT is an open source code and a further development of the original ORBIT code that was created and used for design, studies, and commissioning of the SNS ring. The PyORBIT code, like the original one, has a two-layer structure. C++ is used to perform time-consuming computations, and the program flow is controlled from a Python language shell. The flexible structure makes it possible to use PyORBIT also for linac dynamics simulations. A benchmark of PyORBIT with Parmila and the XAL Online model is presented.

43 PARTICLE ACCELERATORS↗
DOE OSTI2021

Modeling laser-driven ion acceleration with deep learning

Developments in machine learning promise to ameliorate some of the challenges of modeling complex physical systems through neural-network-based surrogate models. High-intensity, short-pulse lasers can be used to accelerate ions to mega-electronvolt energies, but to model such interactions requires computationally expensive techniques such as particle-in-cell simulations. Multilayer neural networks allow one to take a relatively sparse ensemble of simulations and generate a surrogate model that can be used to rapidly search the parameter space of interest. In this study, we created an ensemble of over 1,000 simulations modeling laser-driven ion acceleration and developed a surrogate to study the resulting parameter space. A neural-network-based approach allows for rapid feature discovery not possible for traditional parameter scans given the computational cost. A notable observation made during this study was the dependence of ion energy on the pre-plasma gradient length scale. While this methodology harbors great promise for ion acceleration, it has ready application to all topics in which large-scale parameter scans are restricted by significant computational cost or relatively large, but sparse, domains.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

lasy v0.1

Lasy is a Python library that facilitates numerical simulations of the interaction of complex laser pulses with plasmas. More specifically, lasy offers many ways to define complex laser pulses (e.g. from commonly-known analytical formulas, from experimental measurements, etc.) and then computes and saves the corresponding laser field in a standardized file. This file can then be read by external plasma simulation codes, that will then simulate the interaction between the complex laser pulse and a plasma. Lasy is useful for instance for numerical simulations of plasma-based particle accelerators, plasma-based laser amplification, etc. Lasy avoids having to duplicate the code that defines complex laser pulses in each plasma simulation codes, and instead provides a trusted, centralized implementation.

Lehe, Remi↗