Cyclotron beam dynamic simulations in MATLAB
MATLAB is useful for beam dynamic simulations in cyclotrons. Programming in an easy-to-use environment permits creation of models in a short space of time. Advanced graphical tools of MATLAB give good visualization features to created models. The beam dynamic modeling results with an example of two different cyclotron designs are presented. Programming with MATLAB opens wide possibilities of the development of the complex program, able to perform complete block of calculations for the design of the accelerators
Parallel beam dynamics simulation of linear accelerators
Qiang, Ji; Ryne, Robert D.
2002-01-01
In this paper we describe parallel particle-in-cell methods for the large scale simulation of beam dynamics in linear accelerators. These techniques have been implemented in the IMPACT (Integrated Map and Particle Accelerator Tracking) code. IMPACT is being used to study the behavior of intense charged particle beams and as a tool for the design of next-generation linear accelerators. As examples, we present applications of the code to the study of emittance exchange in high intensity b...
Kikuchi Takashi
2013-11-01
Full Text Available In a final stage of an accelerator system for heavy ion inertial fusion (HIF, pulse shaping and beam current increase by bunch compression are required for effective pellet implosion. A compact simulator with an electron beam was constructed to understand the beam dynamics. In this study, we investigate theoretically and numerically the beam dynamics for the extreme bunch compression in the final stage of HIF accelerator complex. The theoretical and numerical results implied that the compact experimental device simulates the beam dynamics around the stagnation point for initial low temperature condition.
Periodic Poisson model for beam dynamics simulation
Dohlus, M.; Henning, Ch.
2016-03-01
A method is described to solve the Poisson problem for a three dimensional source distribution that is periodic into one direction. Perpendicular to the direction of periodicity a free space (or open) boundary condition is realized. In beam physics, this approach allows us to calculate the space charge field of a continualized charged particle distribution with periodic pattern. The method is based on a particle-mesh approach with equidistant grid and fast convolution with a Green's function. The periodic approach uses only one period of the source distribution, but a periodic extension of the Green's function. The approach is numerically efficient and allows the investigation of periodic- and pseudoperiodic structures with period lengths that are small compared to the source dimensions, for instance of laser modulated beams or of the evolution of micro bunch structures. Applications for laser modulated beams are given.
Beam dynamics simulations using a parallel version of PARMILA
The computer code PARMILA has been the primary tool for the design of proton and ion linacs in the United States for nearly three decades. Previously it was sufficient to perform simulations with of order 10000 particles, but recently the need to perform high resolution halo studies for next-generation, high intensity linacs has made it necessary to perform simulations with of order 100 million particles. With the advent of massively parallel computers such simulations are now within reach. Parallel computers already make it possible, for example, to perform beam dynamics calculations with tens of millions of particles, requiring over 10 GByte of core memory, in just a few hours. Also, parallel computers are becoming easier to use thanks to the availability of mature, Fortran-like languages such as Connection Machine Fortran and High Performance Fortran. We will describe our experience developing a parallel version of PARMILA and the performance of the new code. (author)
Beam dynamics simulations using a parallel version of PARMILA
The computer code PARMILA has been the primary tool for the design of proton and ion linacs in the United States for nearly three decades. Previously it was sufficient to perform simulations with of order 10000 particles, but recently the need to perform high resolution halo studies for next-generation, high intensity linacs has made it necessary to perform simulations with of order 100 million particles. With the advent of massively parallel computers such simulations are now within reach. Parallel computers already make it possible, for example, to perform beam dynamics calculations with tens of millions of particles, requiring over 10 GByte of core memory, in just a few hours. Also, parallel computers are becoming easier to use thanks to the availability of mature, Fortran-like languages such as Connection Machine Fortran and High Performance Fortran. We will describe our experience developing a parallel version of PARMILA and the performance of the new code
Beam dynamics simulations of post low energy beam transport section in RAON heavy ion accelerator
Jin, Hyunchang, E-mail: hcjin@ibs.re.kr; Jang, Ji-Ho; Jang, Hyojae; Hong, In-Seok [Institute for Basic Science, Yuseong-daero 1689-gil, Yuseong-gu, Daejeon (Korea, Republic of)
2016-02-15
RAON (Rare isotope Accelerator Of Newness) heavy ion accelerator of the rare isotope science project in Daejeon, Korea, has been designed to accelerate multiple-charge-state beams to be used for various science programs. In the RAON accelerator, the rare isotope beams which are generated by an isotope separation on-line system with a wide range of nuclei and charges will be transported through the post Low Energy Beam Transport (LEBT) section to the Radio Frequency Quadrupole (RFQ). In order to transport many kinds of rare isotope beams stably to the RFQ, the post LEBT should be devised to satisfy the requirement of the RFQ at the end of post LEBT, simultaneously with the twiss parameters small. We will present the recent lattice design of the post LEBT in the RAON accelerator and the results of the beam dynamics simulations from it. In addition, the error analysis and correction in the post LEBT will be also described.
Beam dynamics simulations and measurements at the Project X Test Facility
Project X, under study at Fermilab, is a multitask high-power superconducting RF proton beam facility, aiming to provide high intensity protons for rare processes experiments and nuclear physics at low energy, and simultaneously for the production of neutrinos, as well as muon beams in the long term. A beam test facility - former known as High Intensity Neutrino Source (HINS) - is under commissioning for testing critical components of the project, e.g. dynamics and diagnostics at low beam energies, broadband beam chopping, RF power generation and distribution. In this paper we describe the layout of the test facility and present beam dynamics simulations and measurements.
Beam dynamics simulations and measurements at the Project X Test Facility
Gianfelice-Wendt, E.; Scarpine, V.E.; Webber, R.C.; /Fermilab
2011-03-01
Project X, under study at Fermilab, is a multitask high-power superconducting RF proton beam facility, aiming to provide high intensity protons for rare processes experiments and nuclear physics at low energy, and simultaneously for the production of neutrinos, as well as muon beams in the long term. A beam test facility - former known as High Intensity Neutrino Source (HINS) - is under commissioning for testing critical components of the project, e.g. dynamics and diagnostics at low beam energies, broadband beam chopping, RF power generation and distribution. In this paper we describe the layout of the test facility and present beam dynamics simulations and measurements.
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.
Beam dynamics simulations for linacs driving short-wavelength FELs
The fast code HOMDYN has been recently developed, in the framework of the TTF (Tesla test facility) collaboration, in order to study the beam dynamics of linacs delivering high brightness beams as those needed for short wavelength Fel experiments. These linacs are typically driven by radio-frequency photo-injectors, where correlated time dependent space charge effects are of great relevance: these effects cannot be studied by standard beam optics codes (TRACE3D, etc.) and they have been modeled so far by means of multi-particle (Pic or quasistatic) codes requiring heavy cpu time and memory allocations. HOMDYN is able to describe the beam generation at the photo-cathode and the emittance compensation process in the injector even running on a laptop with very modest running rimes (less than a minute). In this paper it is showed how this capability of the code is exploited so to model a whole linac up to the point where the space charge dominated regime is of relevance (200 MeV)
Many-beam dynamical simulation of electron backscatter diffraction patterns.
Winkelmann, Aimo; Trager-Cowan, Carol; Sweeney, Francis; Day, Austin P; Parbrook, Peter
2007-01-01
We present an approach for the simulation of complete electron backscatter diffraction (EBSD) patterns where the relative intensity distributions in the patterns are accurately reproduced. The Bloch wave theory is applied to describe the electron diffraction process. For the simulation of experimental patterns with a large field of view, a large number of reflecting planes has to be taken into account. This is made possible by the Bethe perturbation of weak reflections. Very good agreement is obtained for simulated and experimental patterns of gallium nitride GaN{0001} at 20kV electron energy. Experimental features like zone-axis fine structure and higher-order Laue zone rings are accurately reproduced. We discuss the influence of the diffraction of the incident beam in our experiment. PMID:17126489
Simulations of longitudinal beam dynamics of space-charge dominated beams for heavy ion fusion
The longitudinal instability has potentially disastrous effects on the ion beams used for heavy ion driven inertial confinement fusion. This instability is a open-quotes resistive wallclose quotes instability with the impedance coining from the induction modules in the accelerator used as a driver. This instability can greatly amplify perturbations launched from the beam head and can prevent focusing of the beam onto the small spot necessary for fusion. This instability has been studied using the WARPrz particle-in-cell code. WARPrz is a 2 1/2 dimensional electrostatic axisymmetric code. This code includes a model for the impedance of the induction modules. Simulations with resistances similar to that expected in a driver show moderate amounts of growth from the instability as a perturbation travels from beam head to tail as predicted by cold beam fluid theory. The perturbation reflects off the beam tail and decays as it travels toward the beam head. Nonlinear effects cause the perturbation to steepen during reflection. Including the capacitive component of the, module impedance. has a partially stabilizing effect on the longitudinal instability. This reduction in the growth rate is seen in both cold beam fluid theory and in simulations with WARPrz. Instability growth rates for warm beams measured from WARPrz are lower than cold beam fluid theory predicts. Longitudinal thermal spread cannot account for this decrease in the growth rate. A mechanism for coupling the transverse thermal spread to decay of the longitudinal waves is presented. The longitudinal instability is no longer a threat to the heavy ion fusion program. The simulations in this thesis have shown that the growth rate for this instability will not be as large as earlier calculations predicted
Contact dynamics of elasto-plastic thin beams simulated via absolute nodal coordinate formulation
Wang, Qing-Tao; Tian, Qiang; Hu, Hai-Yan
2015-12-01
Under the frame of multibody dynamics, the contact dynamics of elasto-plastic spatial thin beams is numerically studied by using the spatial thin beam elements of absolute nodal coordinate formulation (ANCF). The internal force of the elasto-plastic spatial thin beam element is derived under the assumption that the plastic strain of the beam element depends only on its longitudinal deformation. A new body-fixed local coordinate system is introduced into the spatial thin beam element of ANCF for efficient contact detection in the contact dynamics simulation. The linear isotropic hardening constitutive law is used to describe the elasto-plastic deformation of beam material, and the classical return mapping algorithm is adopted to evaluate the plastic strains. A multi-zone contact approach of thin beams previously proposed by the authors is also introduced to detect the multiple contact zones of beams accurately, and the penalty method is used to compute the normal contact force of thin beams in contact. Four numerical examples are given to demonstrate the applicability and effectiveness of the proposed elasto-plastic spatial thin beam element of ANCF for flexible multibody system dynamics.
Start-to-end simulations for beam dynamics in the RISP heavy-ion accelerator
Kim, Eun-San, E-mail: eskim1@knu.ac.kr [Department of Physics, Kyungpook National University, Deagu (Korea, Republic of); Bahng, JungBae [Department of Physics, Kyungpook National University, Deagu (Korea, Republic of); Hwang, Ji-Gwang [KIRAMS, Seoul (Korea, Republic of); Choi, Bong-Hyuk; Kim, Hye-Jin; Jeon, Dong-O [Institute for Basic Science, Daejeon (Korea, Republic of)
2015-09-11
RAON has been designed as a rare isotope accelerator facility for the Rare Isotope Science Project (RISP). The main accelerator for the in-flight system accelerates uranium and proton beams to 200 MeV/u and 660 MeV, respectively, with a beam power of 400 kW. The front-end system consists of two 28 GHz electron cyclotron resonance ion sources (10 keV/u), a low-energy beam transport (LEBT) line with two 90° bends, a multi-harmonic buncher with three different rf frequencies, a radio-frequency quadrupole (RFQ), and a medium-energy beam transport line (MEBT) with three rebunchers and eight quadrupoles. A driver linac system consisting of Linac-1 and Linac-2 has been designed to optimize the beam and accelerator parameters so as to meet the required design goals. A charge stripper section is located between Linac-1 and Linac-2. To optimize these designs, we performed start-to-end simulations with the beams from the LEBT to Linac-2 using 1 million macroparticles. We present the resulting beam dynamics to evaluate the performance of the accelerator. Our simulation results predict that the transmission rate of the uranium beam is 85.8% from the LEBT to Linac-2. The designed facility is expected to achieve the required beam loss condition of less than 1 W/m. The RAON driver linac lattice design was developed and an overview of the beam dynamics is presented.
SciDAC Advances in Beam Dynamics Simulation: From Light Sources to Colliders
In this paper, we report on progress that has been made in beam dynamics simulation, from light sources to colliders, during the first year of the SciDAC-2 accelerator project 'Community Petascale Project for Accelerator Science and Simulation (ComPASS).' Several parallel computational tools for beam dynamics simulation are described. Also presented are number of applications in current and future accelerator facilities (e.g., LCLS, RHIC, Tevatron, LHC, and ELIC). Particle accelerators are some of most important tools of scientific discovery. They are widely used in high-energy physics, nuclear physics, and other basic and applied sciences to study the interaction of elementary particles, to probe the internal structure of matter, and to generate high-brightness radiation for research in materials science, chemistry, biology, and other fields. Modern accelerators are complex and expensive devices that may be several kilometers long and may consist of thousands of beamline elements. An accelerator may transport trillions of charged particles that interact electromagnetically among themselves, that interact with fields produced by the accelerator components, and that interact with beam-induced fields. Large-scale beam dynamics simulations on massively parallel computers can help provide understanding of these complex physical phenomena, help minimize design cost, and help optimize machine operation. In this paper, we report on beam dynamics simulations in a variety of accelerators ranging from next generation light sources to high-energy ring colliders that have been studied during the first year of the SciDAC-2 accelerator project.
XIAO Chen; HE Yuan; YUAN You-Jin; YAO Qing-Gao; WANG Zhi-Jun; CHANG Wei; LIU Yong; XIA Jia-Wen
2011-01-01
A new SSC-linac system (injector into separated sector cyclotron) is being designed in the HIRFL (heavy ion research facility of Lanzhou). As part of SSC-Linac, the LEBT (low energy beam transport) consists of seven solenoids, four quadrupoles, a bending magnet and an extra multi-harmonic buncher. The total length of this segment is about 7 meters. The beam dynamics in this LEBT has been studied using three-dimensional PIC (particle-in-cell) code BEAMPATH. The simulation results show that the continuous beam from the ion source is first well analyzed by a charge-to-mass selection system, and the beam of the selected charge-to-mass ratio is then efficiently pre-bunched by a multi-harmonic buncher and optimally matched into the RFQ (radio frequency quadrupole) for further acceleration. The principles and effects of the solenoid collimation channel are discussed, and it could limit the beam emittance by changing the aperture size.
Beam dynamics simulation studies of a 100 MeV RF electron linac
We have studied the beam dynamics of a 100 MeV rf electron linear accelerator for neutron generation purposes at APPD, BARC, Mumbai. The linac is designed to provide optimal beams to generate 1012-1014 n/sec for measurement of neutron cross-section of (n, gamma), (n, xn) and (n, f) reactions. An optium choice of beam parameters is performed to reduce the transverse emittance in order to limit the radioactivity and the cost of the linac itself and increase the transmission efficiency. The tracking simulations in the linac include the self field effects in the injector section. The bunch compressor is designed to reduce the bunch length with a careful choice of beam parameters and accelerator parameters such that the transverse emittance is not significantly increased. The linac consists of electron gun followed by an accelerating section. The accelerating section consists of 1.3 GHz TESLA cavities of length 1 m. Beam focusing is ensured by solenoids, quadrupoles. The bunch compressor is a 4-dipole chicane magnet to compress the bunch. This beam dynamics study with ASTRA and ELLEGANT explores the optimized beam parameters and predict the beam quality for the linac operation. (author)
A particle-in-cell mode beam dynamics simulation of medium energy beam transport for the SSC-Linac
A new linear accelerator system, called the SSC-Linac injector, is being designed at HIRFL (the heavy ion research facility of Lanzhou). As part of the SSC-Linac, the medium energy beam transport (MEBT) consists of seven magnetic quadrupoles, a re-buncher and a diagnose box. The total length of this segment is about 1.75 m. The beam dynamics simulation in MEBT has been studied using the TRACK 3D particle- in-cell code, and the simulation result shows that the beam accelerated from the radio frequency quadrupole (RFQ) matches well with the acceptance of the following drift tube linac (DTL) in both the transverse and longitudinal phase spaces, and that most of the particles can be captured by the final sector focusing cyclotron for further acceleration. The longitudinal emittance of the RFQ and the longitudinal acceptance of the DTL was calculated in detail, and a multi-particle beam dynamics simulation from the ion source to the end of the DTL was done to verify the original design. (authors)
Design and beam dynamics simulations of an S-band photocathode rf gun
Kumar, Arvind; Pant, K. K.; Krishnagopal, S.
2002-10-01
We are building an S-band photocathode rf gun as an injector to a 30MeV electron linac for FEL applications. Here we discuss details of design simulations performed using superfish and gdfidl and compare with results of cold tests performed on prototype cells of the photocathode rf gun. We also discuss beam dynamics simulations performed using parmela and report results from simulations to achieve a normalized transverse rms emittance of about 1π mm mrad for a 10ps pulse with 1nC charge in the presence of a solenoid magnetic field used for emittance compensation.
Comparison of Parmela and MAFIA Simulations of Beam Dynamics in High Current Photoinjector
Kurennoy, Sergey S
2004-01-01
A high-current RF photoinjector producing low-emittance electron beam is an important technology for high-power CW FEL. LANL-AES team designed a 2.5-cell, pi-mode, 700-MHz normal-conducting RF photoinjector with magnetic emittance compensation. With the electric field gradients of 7, 7, and 5 MV/m in the three subsequent cells, the photoinjector will produce a 2.5-MeV electron beam with 3-nC charge per bunch and the transverse rms emittance 7 mm-mrad. Beam dynamics in the photoinjector has been modeled in details. In addition to the usual approach, with fields calculated by Superfish-Poisson and beam simulations performed by Parmela, we also used MAFIA group of codes, both to calculate cavity fields and to model beam dynamics with its particle-in-cell module TS. The second way naturally includes wake-field effects into consideration. The simulation results and comparison between two approaches will be presented.
Object-Oriented Parallel Particle-in-Cell Code for Beam Dynamics Simulation in Linear Accelerators
Qiang, J.; Ryne, R.D.; Habib, S.; Decky, V.
1999-11-13
In this paper, we present an object-oriented three-dimensional parallel particle-in-cell code for beam dynamics simulation in linear accelerators. A two-dimensional parallel domain decomposition approach is employed within a message passing programming paradigm along with a dynamic load balancing. Implementing object-oriented software design provides the code with better maintainability, reusability, and extensibility compared with conventional structure based code. This also helps to encapsulate the details of communications syntax. Performance tests on SGI/Cray T3E-900 and SGI Origin 2000 machines show good scalability of the object-oriented code. Some important features of this code also include employing symplectic integration with linear maps of external focusing elements and using z as the independent variable, typical in accelerators. A successful application was done to simulate beam transport through three superconducting sections in the APT linac design.
Efficiency optimization of a fast Poisson solver in beam dynamics simulation
Zheng, Dawei; Pöplau, Gisela; van Rienen, Ursula
2016-01-01
Calculating the solution of Poisson's equation relating to space charge force is still the major time consumption in beam dynamics simulations and calls for further improvement. In this paper, we summarize a classical fast Poisson solver in beam dynamics simulations: the integrated Green's function method. We introduce three optimization steps of the classical Poisson solver routine: using the reduced integrated Green's function instead of the integrated Green's function; using the discrete cosine transform instead of discrete Fourier transform for the Green's function; using a novel fast convolution routine instead of an explicitly zero-padded convolution. The new Poisson solver routine preserves the advantages of fast computation and high accuracy. This provides a fast routine for high performance calculation of the space charge effect in accelerators.
IMPACT simulation and the SNS linac beam
Multi-particle tracking simulations for the SNS linac beam dynamics studies are performed with the IMPACT code. Beam measurement results are compared with the computer simulations, including beam longitudinal halo and beam losses in the superconducting linac, transverse beam Courant-Snyder parameters and the longitudinal beam emittance in the linac. In most cases, the simulations show good agreement with the measured results
Gladkikh, P I; Karnaukhov, I M
2002-01-01
The feasibility of the development of intense X-ray sources based on Compton scattering in laser-electron storage rings is discussed. The results of the electron beam dynamics simulation involving Compton and intrabeam scattering are presented.
2D beam dynamics simulation in linear mode LPWA channel with pre-modulation stage
Laser plasma wakefield acceleration (LPWA) is one of most popular novel methods of acceleration. The acceleration process differs significantly for linear LPWA mode and bubble (non-linear) modes. The LPWA has two serous disadvantages as very high energy spread and low part of electrons trapped into acceleration. The energy spectrum better than 10% does not observed anyone in simulations or experiments without of especial plasma density distribution. Such simulations and first experiments was done for bubble mode with different injection methods as varying of plasma density into bunching sub-stage, ponderomotive injection, etc. But linear mode LPWA is also very interesting to design a compact hundreds-MeV accelerator. 2D beam dynamics in linear mode LPWA is discussed in this report. The waveguide and klystron type beam pre-modulation schemes are studied. The simulation shows that the klystron type pre-modulation can to gives the energy spectrum better than 1.5% for 200...300 MeV beam and to achieve the capturing coefficient 70...80%
Synergia an advanced object-oriented framework for beam dynamics simulation
Dechow, Douglas R; Spentzouris, Panagiotis; Stoltz, Peter
2005-01-01
Synergia is a 3-D, parallel, particle-in-cell beam dynamics simulation toolkit. At heart of the software development effort is the integration of two extant object-oriented accelerator modeling frameworksImpact written in Fortran 90 and mxyptlk written in C++so that they may be steered by a third, a more flexible human interface framework, written in Python. Recent efforts are focused on the refactoring of the Impact-Fortran 90 codes in order to expose more loosely-coupled interfaces to the Python interface framework.
Beam dynamics simulation of the EURISOL radioactive beam superconducting LINAC option
The EURISOL project is aimed at designing a preliminary setup of the next-generation European Rare Isotope Separation On-Line facility. For most of the scientific applications a post-accelerator is needed in order to accelerate the rare isotopes to a wide range of energies. The EURISOL main dictated requirements for the post accelerator are: a: acceleration of all possible radioactive beams. b: to all possible final energies up to 100 MeV/u, c: finely tunable final energy. d: very good beam-quality up to at least 10 MeV/u. e: high transmission efficiency aiming to 100%. f: mass resolution capable of isobar separation and g: affordable construction and operation cost. For the design of the post-accelerator a test case of a radioactive beam of 132Sn25+, generated by a charge breeder, was dictated. In this work we propose a design of an independently phased superconducting cavities linear accelerator for this EURISOL post-accelerator. Recent developments in the field allow high transmission efficiency after stripping due to multi-charge beam transport [2] and high acceleration gradient. Both achievements enable the design of a competitive construction and low operation cost post-accelerator. The rare isotopes are produced in a conventional ISOL target by fragmentation, spallation, or fission techniques. Upon extraction, a desired isotope, as 132Sn, is selected by a high-resolution mass separator and introduced into a charge breeder. The ions are then accelerated by the post-accelerator
A Fast Parallel Poisson Solver on Irregular Domains Applied to Beam Dynamic Simulations
Adelmann, A; Ineichen, Y
2009-01-01
We discuss the scalable parallel solution of the Poisson equation within a Particle-In-Cell (PIC) code for the simulation of electron beams in particle accelerators of irregular shape. The problem is discretized by Finite Differences. Depending on the treatment of the Dirichlet boundary the resulting system of equations is symmetric or `mildly' nonsymmetric positive definite. In all cases, the system is solved by the preconditioned conjugate gradient algorithm with smoothed aggregation (SA) based algebraic multigrid (AMG) preconditioning. We investigate variants of the implementation of SA-AMG that lead to considerable improvements in the execution times. We demonstrate good scalability of the solver on distributed memory parallel processor with up to 2048 processors. We also compare our SAAMG-PCG solver with an FFT-based solver that is more commonly used for applications in beam dynamics.
J. Rodnizki, D. Berkovits, K. Lavie, I. Mardor, A. Shor and Y. Yanay (Soreq NRC, Yavne), K. Dunkel, C. Piel (ACCEL, Bergisch Gladbach), A. Facco (INFN/LNL, Legnaro, Padova), V. Zviagintsev (TRIUMF, Vancouver)
AbstractBeam dynamics simulations of SARAF (Soreq Applied Research Accelerator Facility) superconducting RF linear accelerator have been performed in order to establish the accelerator design. The multi-particle simulation includes 3D realistic electromagnetic field distributions, space charge forces and fabrication, misalignment and operation errors. A 4 mA proton or deuteron beam is accelerated up to 40 MeV with a moderated rms emittance growth and a high real-estate gradient of 2 MeV/m. An envelope of 40,000 macro-particles is kept under a radius of 1.1 cm, well below the beam pipe bore radius. The accelerator design of SARAF is proposed as an injector for the EURISOL driver accelerator. The Accel 176 MHZ β0=0.09 and β0=0.15 HWR lattice was extended to 90 MeV based on the LNL 352 MHZ β0=0.31 HWR. The matching between both lattices ensures smooth transition and the possibility to extend the accelerator to the required EURISOL ion energy.
The target of the present paper is the study of chirality effects in molecular dynamics from both a theoretical and an experimental point of view under the hypothesis of a molecular dynamics mechanism as the origin of chiral discrimination. This is a fundamental problem per se, and of possible relevance for the problem of the intriguing homochirality in Nature, so far lacking satisfactory explanations. We outline the steps that have been taken so far toward this direction, motivated by various experimental studies of supersonic molecular beams carried out in this laboratory, such as the detection of aligned oxygen in gaseous streams and further evidence on nitrogen, benzene and various hydrocarbons, showing the insurgence of molecular orientation in the dynamics of molecules in flows and in molecular collisions. Chiral effects are theoretically demonstrated to show up in the differential scattering of oriented molecules, also when impinging on surfaces. Focus on possible mechanisms for chiral bio-stereochemistry of oriented reactants may be of pre-biotical interest, for example when flowing in atmospheres of rotating bodies, specifically the planet Earth, as well as in vortex motions of celestial objects. Molecular dynamics simulations and experimental verifications of the hypothesis are reviewed and objectives of future research activity proposed.
Aquilanti, Vincenzo; Grossi, Gaia; Lombardi, Andrea; Maciel, Glauciete S; Palazzetti, Federico [Dipartimento di Chimica, Universita di Perugia, Via Elce di Sotto 8, 06123 Perugia (Italy)], E-mail: abulafia@dyn.unipg.it
2008-10-15
The target of the present paper is the study of chirality effects in molecular dynamics from both a theoretical and an experimental point of view under the hypothesis of a molecular dynamics mechanism as the origin of chiral discrimination. This is a fundamental problem per se, and of possible relevance for the problem of the intriguing homochirality in Nature, so far lacking satisfactory explanations. We outline the steps that have been taken so far toward this direction, motivated by various experimental studies of supersonic molecular beams carried out in this laboratory, such as the detection of aligned oxygen in gaseous streams and further evidence on nitrogen, benzene and various hydrocarbons, showing the insurgence of molecular orientation in the dynamics of molecules in flows and in molecular collisions. Chiral effects are theoretically demonstrated to show up in the differential scattering of oriented molecules, also when impinging on surfaces. Focus on possible mechanisms for chiral bio-stereochemistry of oriented reactants may be of pre-biotical interest, for example when flowing in atmospheres of rotating bodies, specifically the planet Earth, as well as in vortex motions of celestial objects. Molecular dynamics simulations and experimental verifications of the hypothesis are reviewed and objectives of future research activity proposed.
Three-dimensional (3D) ordering of a charged-particle beams circulating in a storage ring is systematically studied with a molecular-dynamics simulation code. An ion beam can exhibit a 3D ordered configuration at ultralow temperature as a result of powerful 3D laser cooling. Various unique characteristics of the ordered beams, different from those of crystalline beams, are revealed in detail, such as the single-particle motion in the transverse and longitudinal directions, and the dependence of the tune depression and the Coulomb coupling constant on the operating points
Yuri, Yosuke, E-mail: yuri.yosuke@jaea.go.jp [Takasaki Advanced Radiation Research Institute, Japan Atomic Energy Agency, 1233 Watanuki-machi Takasaki, Gunma 370-1292 Japan (Japan)
2015-06-29
Three-dimensional (3D) ordering of a charged-particle beams circulating in a storage ring is systematically studied with a molecular-dynamics simulation code. An ion beam can exhibit a 3D ordered configuration at ultralow temperature as a result of powerful 3D laser cooling. Various unique characteristics of the ordered beams, different from those of crystalline beams, are revealed in detail, such as the single-particle motion in the transverse and longitudinal directions, and the dependence of the tune depression and the Coulomb coupling constant on the operating points.
Computer simulation of low-energy high-current electron beam dynamics in a long plasma-filled diode
Results of computer simulation of low-energy high-current electron beam dynamics in a low-impedance system consisting of a diode with a long plasma anode, just siding with an explosive emission cathode and an auxiliary thermionic cathode are presented. Plasma anode plays simultaneously a role of the transport channel providing charge neutralization of high-current beam and is created by means of the residual gas ionisation by low-current, low-voltage electron beam emitted from the auxiliary cathode in an external longitudinal magnetic field. The main peculiarities of the beam-plasma system are discussed: 1) the formation of the beam of currents exceeding the limiting Alven's ones; 2) the formation of paramagnetic states of the beam under condition of beam charge density close to the plasma density. These peculiarities complicate beam-plasma interaction significantly due to sharp nonuniform distribution of the beam current density, significant transverse motion of the beam electrons and redistribution of ion plasma density under the influence of high-current electron beam fields. Computer simulation was performed using electromagnetic PIC code KARAT
Development of RFQ particle dynamics simulation tools and validation with beam tests
Two different strategies of designing RFQs have been introduced. The analytic description of the electric fields inside the quadrupole channel has been derived and the two term simplification was shown as well as the limitation of these approaches. The main work of this thesis was the implementation and analysis of a multigrid Poisson solver to describe the potential and electric field of RFQs which are needed to simulate the particle dynamics accurately. The main two ingredients of a multigrid Poisson solver are the ability of a Gauss-Seidel iteration method to smooth the error of an approximation within a few iteration steps and the coarse grid principle. The smoothing corresponds to a damping of the high frequency components of the error. After the smoothing, the error term can well be approximated on a coarser grid in which the low frequency components of the error on the fine grid are converted to high frequency errors on the coarse grid which can be damped further with the same Gauss-Seidel method. After implementation, the multigrid Poisson solver was analyzed using two different type of test problems: with and without a charge density. As a charge density, a homogeneously charged ball and cylinder were used to represent the bunched and unbunched beam and placed inside a quadruple channel. The solver showed a good performance. Next, the performance of the solver to calculate the external potentials (and fields) of RFQs was analyzed. Closing the analysis of the external field, the transmission and fraction of accelerated particles of the set of 12 RFQs for the two different methods were shown. In the last chapter of this thesis some experimental work on the MAFF (Munich Accelerator for Fission Fragments) IH-RFQ is described. The MAFF RFQ was designed to accelerate very neutron-rich fission fragments for various experiments. The machine was assembled in Frankfurt and a beam test stand was built. As a part of this thesis the shunt impedance of the structure was
Development of RFQ particle dynamics simulation tools and validation with beam tests
Maus, Johannes M.
2010-07-01
Two different strategies of designing RFQs have been introduced. The analytic description of the electric fields inside the quadrupole channel has been derived and the two term simplification was shown as well as the limitation of these approaches. The main work of this thesis was the implementation and analysis of a multigrid Poisson solver to describe the potential and electric field of RFQs which are needed to simulate the particle dynamics accurately. The main two ingredients of a multigrid Poisson solver are the ability of a Gauss-Seidel iteration method to smooth the error of an approximation within a few iteration steps and the coarse grid principle. The smoothing corresponds to a damping of the high frequency components of the error. After the smoothing, the error term can well be approximated on a coarser grid in which the low frequency components of the error on the fine grid are converted to high frequency errors on the coarse grid which can be damped further with the same Gauss-Seidel method. After implementation, the multigrid Poisson solver was analyzed using two different type of test problems: with and without a charge density. As a charge density, a homogeneously charged ball and cylinder were used to represent the bunched and unbunched beam and placed inside a quadruple channel. The solver showed a good performance. Next, the performance of the solver to calculate the external potentials (and fields) of RFQs was analyzed. Closing the analysis of the external field, the transmission and fraction of accelerated particles of the set of 12 RFQs for the two different methods were shown. In the last chapter of this thesis some experimental work on the MAFF (Munich Accelerator for Fission Fragments) IH-RFQ is described. The MAFF RFQ was designed to accelerate very neutron-rich fission fragments for various experiments. The machine was assembled in Frankfurt and a beam test stand was built. As a part of this thesis the shunt impedance of the structure was
Developing models for simulation of pinched-beam dynamics in heavy ion fusion. Revision 1
For heavy-ion fusion energy applications, Mark and Yu have derived hydrodynamic models for numerical simulation of energetic pinched-beams including self-pinches and external-current pinches. These pinched-beams are applicable to beam propagation in fusion chambers and to the US High Temperature Experiment. The closure of the Mark-Yu model is obtained with adiabatic assumptions mathematically analogous to those of Chew, Goldberger, and Low for MHD. Features of this hydrodynamic beam model are compared with a kinetic treatment
Numerical simulations of driving beam dynamics in the plasma wakefield accelerator
Novel plasma based acceleration devices have become the subject of active research because of their ability to support acceleration gradients in excess of 10 GeV/m. The plasma wakefield accelerator (PWFA) is one such device which consists of an intense electron beam (the primary beam) whose purpose is to excite a plasma wave which, in turn, accelerates a trailing electron bunch (the secondary beam). Two issues of current interest in the PWFA are (1) the equilibrium and stability of the driving beam and (2) the effect of the wakefield on the quality of the trailing electron bunch. In the UCLA experiment, a question of particular interest is the equilibrium state of the driving electron beam. Two intriguing suggestions have been made. The first is that in the limit that the beam density greatly exceeds the plasma density, the plasma electrons will be completely expelled from the axis. The second is that, in parameter regimes of interest, the driving beam will experience a severe radial pinching force. In order to investigate these assertions, the authors first consider the envelope equation for an electron beam propagating in a plasma with nb ≥ np. They then compare numerical solutions of this equation to results obtained via two-dimensional axisymmetric (r,z) particle simulation using the GRIEZR particle simulation code
Simulation of Electron Beam Dynamics in a Nonmagnetized High-Current Vacuum Diode
Anishchenko, Sergey
2016-01-01
The electron beam dynamics in a nonmagnetized high-current vacuum diode is analyzed for different cathode-anode gap geometries. The conditions enabling to achieve the minimal {initial} momentum spread in the electron beam are found out. A drastic rise of current density in a vacuum diode with a ring-type cathode is described. The effect is shown to be caused by electrostatic repulsion.
Gao Xiguang; Song Yingdong; Sun Zhigang; Hu Xuteng
2010-01-01
A multiscale method for simulating the dynamic response of ceramic matrix composite (CMC) with matrix cracks is developed.At the global level,the finite element method is employed to simulate the dynamic response ofa CMC beam.While at the local level,the multiscale mechanical method is used to estimate the stress/strain response of the material.A distributed computing system is developed to speed up the simulation.The simulation of dynamic response of a Nicalon/CAS-Ⅱ beam being subjected to harmonic loading is performed as a numerical example.The results show that both the stress/strain responses under tension and compressive loading are nonlinear.These conditions result in a different response compared with that of elastic beam,such as:1) the displacement response is not symmetric about the axis of time;2) in the condition of small external load,the response at first order natural frequency is limited within a finite range;3) decreasing the matrix crack space will increase the displacement response of the beam.
Macroparticle simulation plays an important role in modern accelerator design and operation. Most linear rf accelerators have been designed based on macroparticle simulations using longitudinal position as the independent variable. In this paper, we have done a systematic comparison between using longitudinal position as the independent variable and using time as the independent variable in macroparticle simulations. We have found that, for an rms-matched beam, the maximum relative moment difference for second, fourth moments and beam maximum amplitudes between these two types of simulations is 0.25 percent in a 10 m reference transport system with physical parameters similar to the Spallation Neutron Source linac design. The maximum z-to-t transform error in the space-charge force calculation of the position dependent simulation is about 0.1 percent in such a system. This might cause a several percent error in a complete simulation of a linac with a length of hundreds of meters. Furthermore, the error may be several times larger in simulations of mismatched beams. However, if such errors are acceptable to the linac designer, then one is justified in using position dependent macroparticle simulations in this type of linac design application
Beam dynamics simulation of injector for high power CW electron linac in PNC
The injector consists of a 200 kV DC gun, a RF chopper, a chopper slit, a prebuncher and a buncher. Solenoid coils covered from the exit of gun to accelerating tube 1 except between the RF chopper and chopper slit. Beam trajectories are simulated by PARMELA in order to design the injector. In this report, two simulation results are shown. One is for a beam trajectory from gun to solenoid coils. There is thick concrete wall between gun to RF chopper. Low energy electrons are transported through long solenoid coil area. The other is for a chopper part. The novel chopper system is designed to reduce the emittance growth. (author)
Non-Linear Time-Domain Simulations of the RF Station-Beam Dynamics Interaction for the LHC
Non-linear time-domain simulations have been developed for the Positron-Electron Project (PEP-II) and the Large Hadron Collider (LHC). These simulations capture the dynamic behavior of the RF station-beam interaction and are structured to reproduce the technical characteristics of the system (noise contributions, non-linear elements, and more). As such, they provide useful results and insight for the development and design of future LLRF feedback systems. They are also a valuable tool for the study of diverse longitudinal beam dynamics effects such as coupled-bunch impedance driven instabilities and single bunch longitudinal emittance growth. Results from these studies and related measurements from PEP-II and LHC have been presented in multiple places. This report presents an example of the time-domain simulation implementation for the LHC.
Study on Size-Dependent Young’s Modulus of a Silicon Nano beam by Molecular Dynamics Simulation
Young’s modulus of a silicon nano beam with a rectangular cross-section is studied by molecular dynamics method. Dynamic simulations are performed for doubly clamped silicon nano beams with lengths ranging from 4.888 to 12.491 nm and cross-sections ranging from 1.22 nm ×1.22 nm to 3.39 nm × 3.39 nm. The results show that Young’s moduli of such small silicon nano beams are much higher than the value of Young’s modulus for bulk silicon. Moreover, the resonant frequency and Young’s modulus of the Si nano beam are strongly dependent not only on the size of the nano beam but also on surface effects. Young’s modulus increases significantly with the decreasing of the thickness of the silicon nano beam. This result qualitatively agrees with one of the conclusions based on a semi continuum model, in which the surface relaxation and the surface tension were taken into consideration. The impacts of the surface reconstruction with (2 ×1) dimmers on the resonant frequency and Young’s modulus are studied in this paper too. It is shown that the surface reconstruction makes the silicon nano beam stiffer than the one without the surface reconstruction, resulting in a higher resonant frequency and a larger Young’s modulus
Beam dynamics simulations for Gaussian and flat-top laser pulses at PITZ
The Photo Injector Test Facility at DESY, Zeuthen site (PITZ), has been built in order to develop and optimize electron sources for Free Electron Lasers (FELs) like FLASH and the European XFEL. The electron beam is generated by photoemission initiated with a laser pulse having a flat-top temporal profile. Compared to the Gaussian one, such a flat-top profile yields smaller transverse projected emittance. In order to estimate the difference between the two cases, systematic simulations for 500 pC bunch charge are presented. Dependences of electron beam properties, like beam momentum, transverse beam size, phase space, emittance, on various machine parameters, e.g. gun phase, solenoid current, are shown as well.
Kim, Hojin; Strachan, Alejandro [School of Materials Engineering and Birck Nanotechnology Center, Purdue University, West Lafayette, Indiana 47907 (United States)
2015-11-28
We use large-scale molecular dynamics (MD) to characterize fluid damping between a substrate and an approaching beam. We focus on the near contact regime where squeeze film (where fluid gap is comparable to the mean free path of the gas molecules) and many-body effects in the fluid become dominant. The MD simulations provide explicit description of many-body and non-equilibrium processes in the fluid as well as the surface topography. We study how surface roughness and beam width increases the damping coefficient due to their effect on fluid mobility. We find that the explicit simulations are in good agreement with prior direct simulation Monte Carlo results except at near-contact conditions where many-body effects in the compressed fluid lead the increased damping and weaker dependence on beam width. We also show that velocity distributions near the beam edges and for short gaps deviate from the Boltzmann distribution indicating a degree of local non-equilibrium. These results will be useful to parameterize compact models used for microsystem device-level simulations and provide insight into mesoscale simulations of near-contact damping.
无
2007-01-01
Based on the electron dynamic diffraction, phase shift of the exit wave function vs misorientation of the incident electron beam from the exact zone axis has been calculated for the [001] oriented copper. The result shows that the peak of phase shift is the maximum at the atom position as the electron beam along the exact [001] zone axis, and the peak value of phase shift decreases as increases of the misorientation. At small misorientation, i.e. less than 5 degree, change of the phase shift is minimal. The peak value of phase shift decreases significantly when the incident beam deviates form the zone axis over 10 degree and the exit wave has a planar configuration as the misoriention angle arrives ～17 degree. The effect of this phase shift characteristics on the information extracted from the hologram has also been considered.
Simulative research on the anode plasma dynamics in the high-power electron beam diode
Anode plasma generated by electron beams could limit the electrical pulse-length, modify the impedance and stability of diode, and affect the generator to diode power coupling. In this paper, a particle-in-cell code is used to study the dynamics of anode plasma in the high-power electron beam diode. The effect of gas type, dynamic characteristic of ions on the diode operation with bipolar flow model are presented. With anode plasma appearing, the amplitude of diode current is increased due to charge neutralizations of electron flow. The lever of neutralization can be expressed using saturation factor. At same pressure of the anode gas layer, the saturation factor of CO2 is bigger than the H2O vapor, namely, the generation rate of C+ ions is larger than the H+ ions at the same pressure. The transition time of ions in the anode-cathode gap could be used to estimate the time of diode current maximum
Ekdahl, Carl August Jr. [Los Alamos National Lab. (LANL), Los Alamos, NM (United States)
2014-10-14
Beam dynamics issues are assessed for a new linear induction electron accelerator being designed for flash radiography of large explosively driven hydrodynamic experiments. Special attention is paid to equilibrium beam transport, possible emittance growth, and beam stability. It is concluded that a radiographic quality beam will be produced possible if engineering standards and construction details are equivalent to those on the present radiography accelerators at Los Alamos.
A Fast Parallel Poisson Solver on Irregular Domains Applied to Beam Dynamic Simulations
Adelmann, A.; P. Arbenz; Ineichen, Y.
2009-01-01
We discuss the scalable parallel solution of the Poisson equation within a Particle-In-Cell (PIC) code for the simulation of electron beams in particle accelerators of irregular shape. The problem is discretized by Finite Differences. Depending on the treatment of the Dirichlet boundary the resulting system of equations is symmetric or `mildly' nonsymmetric positive definite. In all cases, the system is solved by the preconditioned conjugate gradient algorithm with smoothed aggregation (SA) b...
Beam dynamics simulation of HEBT for the SSC-linac injector
LI Xiao-Ni; YUAN You-Jin; XIAO Chen; HE Yuan; WANG Zhi-Jun; SHENG Li-Na
2012-01-01
The SSC-linac (a new injector for the Separated Sector Cyclotron) is being designed in the HIRFL (Heavy Ion Research Facility in Lanzhou) system to accelerate 238U34+ from 3.72 keV/u to 1.008 MeV/u.As a part of the SSC-linac injector,the HEBT (high energy beam transport) has been designed by using the TRACE-3D code and simulated by the 3D PIC (particle-in-cell) Track code.The total length of the HEBT is about 12 meters and a beam line of about 6 meters are shared with the exiting beam line of the HIRFL system.The simulation results show that the particles can be delivered efficiently in the HEBT and the particles at the exit of the HEBT well match the acceptance of the SSC for further acceleration.The dispersion is eliminated absolutely in the HEBT.The space-charge effect calculated by the Track code is inconspicuous.According to the simulation,more than 60 percent of the particles from the ion source can be transported into the acceptance of the SSC.
Numerical Simulations of Nonlinear Dynamics of Electron Cyclotron Maser with a Straight Beam
KONG Ling-Bao; HOU Zhi-Ling
2011-01-01
An electron cyclotron maser based on anomalous Doppler effect (ADECM) with an initially axial beam velocity is considered,and the nonlinear equation of beam-wave interaction is presented.With the numerical methods,the nonlinear dynamics of the ADECM is investigated.It is shown that the saturated interaction efficiency of the ADECM approaches 90％ and the interaction length for the saturated efficiency spans about 5-20cm.The results may be of importance for designing a compact device in applications in microwave generations or microwave heating of ceramic laminates.In the late 1950s,the theoretical studies on the instability of electron cyclotron maser based on normal Doppler effect (NDECM) were performed almost simultaneously by Gaponov,[1] Twiss,[2] and Schneider.[3] Their discoveries have resulted in the most successful fast-wave devices such as the gyrotron and variants.[4,5] The possible applications of microwaves span a wide range of technologies such as in thermonuclear fusion energy,charged particle accelerations,radar systems,and processing of advanced ceramics.[6-16]%An electron cyclotron maser based on anomalous Doppler effect (ADECM) with an initially axial beam velocity is considered, and the nonlinear equation of beam-wave interaction is presented. With the numerical methods, the nonlinear dynamics of the ADECM is investigated. It is shown that the saturated interaction efficiency of the ADECM approaches 90% and the interaction length for the saturated efficiency spans about 5-20 cm. The results may be of importance for designing a compact device in applications in microwave generations or microwave heating of ceramic laminates.
A computer program LEADS, written in about 6600 statements of MS FORTRAN 5.1 language, is suitable to run in IBM PC and other compatible computers. Program LEADS can make simulation for beam optical systems consisting of three-tube einzel lenses, three-aperture einzel lenses, two-tube lenses, uniform field DC accelerating tubes, magnetic and electrostatic quadrupoles, bending magnets, single-gap rf resonators, two-gap rf resonators (QWR) and three-gap rf resonators (SLR). Multi-particle tracking and matrix multiplication are used to calculate the beam transport. Monte Carlo techniques are adopted to generate randomly the initial particle coordinates in the phase spaces, and Powell nonlinear optimization routines are incorporated in the codes to search the given optical conditions. The calculated results can be displayed graphically on the computer monitors
Beam dynamics, efficiency and power of the SLAC lasertron: simulation results
Results are described for the computer simulation of the SLAC proof of principle lasertron device with a conventional single gap output cavity, using the 2D relativistic field and particle code called MASK. The rf to beam power efficiency is calculated for different power levels, dc voltages and optical pulse lengths. The calculated efficiency at the initial operating point of 50 MW beam power, 400 kV, and with 60 picosecond optical pulse duration, is 66%. The maximum rf power at 400 kV is about 50 MW. At 600 kV the maximum power increases to about 110 MW, but the efficiency at low power is not much changed from what it was at 400 kV. The simulation calculation does not take into account loss of rf power due to backscattered electrons nor the full effects of the impedance of the accelerating gap. A calculation of the efficiency of the lasertron with a double output cavity has been carried out, and generally yields efficiencies about 10 percentage points higher than the single cavity simulation
Tecker, F.
2016-01-01
The course gives a summary of longitudinal beam dynamics for both linear and circular accelerators. After discussing different types of acceleration methods and synchronism conditions, it focuses on the particle motion in synchrotrons.
Tecker, F
2014-01-01
The course gives a summary of longitudinal beam dynamics for both linear and circular accelerators. After discussing different types of acceleration methods and synchronism conditions, it focuses on the particle motion in synchrotrons.
Simulative research on the anode plasma dynamics in the high-power electron beam diode
Cai, Dan; Liu, Lie; Ju, Jin-Chuan; Zhang, Tian-Yang; Zhao, Xue-Long; Zhou, Hong-Yu [College of Optoelectric Science and Engineering, National University of Defense Technology, Hunan 410073 (China)
2015-07-15
Anode plasma generated by electron beams could limit the electrical pulse-length, modify the impedance and stability of diode, and affect the generator to diode power coupling. In this paper, a particle-in-cell code is used to study the dynamics of anode plasma in the high-power electron beam diode. The effect of gas type, dynamic characteristic of ions on the diode operation with bipolar flow model are presented. With anode plasma appearing, the amplitude of diode current is increased due to charge neutralizations of electron flow. The lever of neutralization can be expressed using saturation factor. At same pressure of the anode gas layer, the saturation factor of CO{sub 2} is bigger than the H{sub 2}O vapor, namely, the generation rate of C{sup +} ions is larger than the H{sup +} ions at the same pressure. The transition time of ions in the anode-cathode gap could be used to estimate the time of diode current maximum.
Beam dynamic simulation and optimization of the CLIC positron source and the capture linac
Bayar, C.; Doebert, S.; Ciftci, A. K.
2016-03-01
The CLIC Positron Source is based on the hybrid target composed of a crystal and an amorphous target. Simulations have been performed from the exit of the amorphous target to the end of pre-injector linac which captures and accelerates the positrons to an energy of 200 MeV. Simulations are performed by the particle tracking code PARMELA. The magnetic field of the AMD is represented in PARMELA by simple coils. Two modes are applied in this study. The first one is accelerating mode based on acceleration after the AMD. The second one is decelerating mode based on deceleration in the first accelerating structure. It is shown that the decelerating mode gives a higher yield for the e+ beam in the end of the Pre-Injector Linac.
Computer simulations of the chaotic dynamics of the Pierce beam endash plasma system
The cold fluid theory of the Pierce beam endash plasma system is modified by the incorporation of warm plasma effects. The controlling parameter α in the cold theory, where α=Lωp/V0, L=diode width, ωp=plasma frequency, and V0=beam velocity at injection, is replaced in the warm theory by an effective value of α involving the thermal velocity. The theory is verified by means of a fluid simulation code; the phase states for a cold plasma, including the chaotic state, are recovered for a warm plasma, but with a shift in values of the bifurcation parameter. Furthermore, in order to include plasma kinetic effects, an extensive electrostatic particle simulation code is developed to model the Pierce system. Among the new physical effects arising in this particle model are the local and global thermalization of electrons by electrostatic waves, and blocking oscillations due to particle reflection and trapping. As the parameter α is decreased, the electric field at the injection point typically changes state as follows: blocking oscillation implies small fluctuations implies quasisteady oscillation implies prechaotic oscillation implies chaos implies blocking oscillation implies dc electric field. The mechanics of chaotic oscillations in the Pierce system are examined with particular regard to kinetic effects. copyright 1996 American Institute of Physics
Assessment of dynamic mechanical behaviour of reinforced concrete beams using a blast simulator
Peroni, Marco; Solomos, George; Caverzan, Alessio; Larcher, Martin; Valsamos, Georgios
2015-09-01
Critical infrastructures may become the target of terrorist bombing attacks or may have to withstand explosive loads due to accidents. The impulsive load connected to explosions is delivered to the structure in a few milliseconds forcing it to respond or fail in a peculiar mode. With reference to the above scientific framework this work presents an innovative apparatus designed and developed at the European Laboratory for Structural Assessment to reproduce a blast pressure history without using explosives. This apparatus is practically a hybrid nitrogen-spring-driven actuator that accelerates masses of up to 100 kg to a maximum velocity of about 25 m/s that impact against the tested structure. The pressure-load history applied to the structure is modulated and reshaped using appropriate layers of elastic soft materials (such as polymeric foams) placed between the specimen and the impacting masses. Specific instrumentation has extensively been utilised to investigate the blast simulator performance and to precisely measure the pressure loads applied to the specimen. A series of tests on real scale reinforced concrete beams/columns (250 × 250 × 2200 mm) has been performed to efficiently assess the performance and potentiality of the new blast simulator. Results are under evaluation. In addition to the experimental work, a series of numerical simulations by means of the explicit FEM code EUROPLEXUS have been carried out to support and improve the equipment design.
Impurity Dynamics under Neutral Beam Injection at TJ-II (simulation)
In this study the simulations of plasma transport under NBI for TJ-II, previously performed, are extended. Since than a considerable number of important modifications have been introduced in the model: change of magnetic configuration, use of experimental initial profiles, expansion of the Data base from NBI calculations and, mainly, a detailed handling of impurities with inclusion of sputtering effects. Moreover there is now a particular emphasis on the analysis of the conditions for discharge collapse and on the possible effects of single beam injection. This analysis of impurity behaviour with sputtering shows that in the expected usual cases there is no radioactive collapse and that if the recycling coefficients remain lower the unity it is always possible to find a strategy for external gas puffing leading to a stationary state, with densities below the limit and efficient NBI absorption (>50%). The radioactive collapse can appear either at high densities (central value higher than 1.4x10''20 m''3), excessive influx of impurities (i. e. with sputtering rates higher than twice the expected values) o for insufficient injected beam power (less than 45 kW). The present study analyses only the 1004464 configuration of TJ-II, but future works will start a systematic scan of configuration using this same model. (Author) 12 Refs
Dynamic acoustic tractor beams
Mitri, F. G.
2015-03-01
Pulling a sphere and vibrating it around an equilibrium position by amplitude-modulation in the near-field of a single finite circular piston transducer is theoretically demonstrated. Conditions are found where a fluid hexane sphere (with arbitrary radius) chosen as an example, centered on the axis of progressive propagating waves and submerged in non-viscous water, experiences an attractive (steady) force pulling it towards the transducer, as well as an oscillatory force forcing it to vibrate back-and-forth. Numerical predictions for the dynamic force illustrate the theory and suggest an innovative method in designing dynamic acoustical tractor beams.
Dynamic acoustic tractor beams
Mitri, F. G., E-mail: F.G.Mitri@ieee.org [Chevron, Area 52 Technology – ETC, Santa Fe, New Mexico 87508 (United States)
2015-03-07
Pulling a sphere and vibrating it around an equilibrium position by amplitude-modulation in the near-field of a single finite circular piston transducer is theoretically demonstrated. Conditions are found where a fluid hexane sphere (with arbitrary radius) chosen as an example, centered on the axis of progressive propagating waves and submerged in non-viscous water, experiences an attractive (steady) force pulling it towards the transducer, as well as an oscillatory force forcing it to vibrate back-and-forth. Numerical predictions for the dynamic force illustrate the theory and suggest an innovative method in designing dynamic acoustical tractor beams.
Ben-Zvi I.; Kuczewski A.; Altinbas, Z.; Beavis, D.; Belomestnykh,; Dai, J. et al
2012-07-01
The Collider-Accelerator Department at Brookhaven National Laboratory is building a high-brightness 500 mA capable Energy Recovery Linac (ERL) as one of its main R&D thrusts towards eRHIC, the polarized electron - hadron collider as an upgrade of the operating RHIC facility. The ERL is in final assembly stages, with injection commisioning starting in October 2012. The objective of this ERL is to serve as a platform for R&D into high current ERL, in particular issues of halo generation and control, Higher-Order Mode (HOM) issues, coherent emissions for the beam and high-brightness, high-power beam generation and preservation. The R&D ERL features a superconducting laser-photocathode RF gun with a high quantum efficiency photoccathode served with a load-lock cathode delivery system, a highly damped 5-cell accelerating cavity, a highly flexible single-pass loop and a comprehensive system of beam instrumentation. In this ICFA Beam Dynamics Newsletter article we will describe the ERL in a degree of detail that is not usually found in regular publications. We will discuss the various systems of the ERL, following the electrons from the photocathode to the beam dump, cover the control system, machine protection etc and summarize with the status of the ERL systems.
Beam Dynamics and Beam Losses - Circular Machines
Kain, V
2016-01-01
A basic introduction to transverse and longitudinal beam dynamics as well as the most relevant beam loss mechanisms in circular machines will be presented in this lecture. This lecture is intended for physicists and engineers with little or no knowledge of this subject.
Lalayan, M V; Sobenin, N P; Shvedunov, V I; Zavadtsev, A A; Dohlus, M
2001-01-01
Some features of accelerating section field computation presented by the development of power and high order mode couplers for TESLA linear collider are considered. The devices mentioned produce electromagnetic field asymmetry in the beam area, thus causing transverse kick. For this kick and its influence on beam under acceleration parameters estimation the dynamics modelling calculations were done. 3D-simulation code MAFIA was used for field computation. These data were further used in beam dynamics calculations by means of TRMTrace code. Standing wave mode was simulated while considering HOM couplers, and travelling wave in case of power couplers. Transverse kicks and focussing forces are calculated for one HOM coupler design and two coaxial FM couplers.
After a bibliographic research on field emission, photoemission and photo-field emission, the principle of the field equations (Poisson and Maxwell's) resolution by the finite element method is developed. The PRIAM program is shown to be efficient (adaptive mesh and refinement in the selected area). Several possibilities exist to reduce the effect of space charge such as the decrease of the laser pulse duration, the increase of the electric field and the application of a magnetic field. Calculations of the transverse emittance for a metallic plan photocathode have been made at different moments of the emission: transverse emittance is small at the beginning and at the end of the emission. It passes by a maximum which can be the origin of the electronic beam explosion for strong field. If a small emittance is wanted, one must illuminate the photocathode by a short pulsed laser
Single-particle beam dynamics in Boomerang
We describe simulations of the beam dynamics in the storage ring (Boomerang), a 3-GeV third-generation light source being designed for the Australian Synchrotron Project[1]. The simulations were performed with the code Goemon[2]. They form the basis for design specifications for storage ring components (apertures, alignment tolerances, magnet quality, etc.), and for determining performance characteristics such as coupling and beam lifetime
Introduction to Transverse Beam Dynamics
Holzer, B J
2014-01-01
In this chapter we give an introduction to the transverse dynamics of the particles in a synchrotron or storage ring. The emphasis is more on qualitative understanding rather than on mathematical correctness, and a number of simulations are used to demonstrate the physical behaviour of the particles. Starting from the basic principles of how to design the geometry of the ring, we review the transverse motion of the particles, motivate the equation of motion, and show the solutions for typical storage ring elements. Following the usual treatment in the literature, we present a second way to describe the particle beam, using the concept of the emittance of the particle ensemble and the beta function, which reflects the overall focusing properties of the ring. The adiabatic shrinking due to Liouville's theorem is discussed as well as dispersive effects in the most simple case.
Simulating Transient Effects of Pulsed Beams on Beam Intercepting Devices
Richter, Herta; Noah Messomo, Etam
2011-01-01
The development in the physics community towards higher beam power through the possibilities of particle accelerators lead to challenges for the developers of elements which are exposed to effect of particle beams (beam intercepting devices = BIDs). For the design of BIDs, the increasing heat load onto these devices due to energetic and focused beams and - in most cases - their highly pulsed nature has to be taken into account. The physics requirements are sometimes opposed to the current state of the art. As one possibility of many in combining the different aspects for these ambitious demands, two highly developed computer programs, namely FLUKA and ANSYS AUTODYN, were joined for this dissertation. The former is a widely enhanced Monte-Carlo-code which specializes on the interaction of particles with static matter, while the latter is a versatile explicit code for the simulation of highly dynamic processes. Both computer programs were developed intensively over years and are still continuously enhanced in o...
Benedetto, E; Borburgh, J; Carli, C; Martini, M; Forte, V
2014-01-01
The CERN PS Booster will be upgraded with an H- injection system. The chicanemagnets for the injection bump ramp-down in 5 ms and generate eddy currents in the inconel vacuum chamber which perturb the homogeneity of the magnetic field. The multipolar field components are extracted from 3D OPERA simulations and are included in the lattice model. The -beating correction is computed all along the ramp and complete tracking simulations including space-charge are performed to evaluate the impact of these perturbations and correction on beam dynamics.
Simulation based analysis of laser beam brazing
Dobler, Michael; Wiethop, Philipp; Schmid, Daniel; Schmidt, Michael
2016-03-01
Laser beam brazing is a well-established joining technology in car body manufacturing with main applications in the joining of divided tailgates and the joining of roof and side panels. A key advantage of laser brazed joints is the seam's visual quality which satisfies highest requirements. However, the laser beam brazing process is very complex and process dynamics are only partially understood. In order to gain deeper knowledge of the laser beam brazing process, to determine optimal process parameters and to test process variants, a transient three-dimensional simulation model of laser beam brazing is developed. This model takes into account energy input, heat transfer as well as fluid and wetting dynamics that lead to the formation of the brazing seam. A validation of the simulation model is performed by metallographic analysis and thermocouple measurements for different parameter sets of the brazing process. These results show that the multi-physical simulation model not only can be used to gain insight into the laser brazing process but also offers the possibility of process optimization in industrial applications. The model's capabilities in determining optimal process parameters are exemplarily shown for the laser power. Small deviations in the energy input can affect the brazing results significantly. Therefore, the simulation model is used to analyze the effect of the lateral laser beam position on the energy input and the resulting brazing seam.
Beam simulation studies of ECR beam extraction and low energy beam transport for FRIB
Ren, Haitao, E-mail: ren@frib.msu.edu; Pozdeyev, Eduard; Lund, Steven M.; Machicoane, Guillaume; Wu, Xiaoyu; Morgan, Glenn [Facility for Rare Isotope Beams, Michigan State University, East Lansing, Michigan 48824 (United States)
2016-02-15
To meet the beam power requirements of 400 kW at the fragmentation target for facility for Rare Isotope Beams (FRIB), simultaneous acceleration of two-charge states should be used for heavier ions. These intense multi-charged ion beams will be produced by a 28 GHz electron cyclotron resonance (ECR) ion source at a high voltage of 35 kV. After extraction, the ion beam will be pre-accelerated to 12 keV/u with a 50 kV platform, transported down to an achromatic charge state selection (CSS) system followed by a vertical transport line, and then injected into a radio frequency quadrupole accelerator. The TRACK code developed at ANL is used to perform the simulations of the ECR beam extraction and low energy beam transport for FRIB. In this study, we include the magnetic field of ECR ion source into simulations. Different initial beam conditions as well as different space charge neutralization levels are tested for the ECR beamline. The beam loss in CSS system and the corresponding protective measures are discussed. The detailed results about the beam dynamic simulation and beam loss in CSS system will be presented in this paper.
Beam simulation studies of ECR beam extraction and low energy beam transport for FRIB
Ren, Haitao; Pozdeyev, Eduard; Lund, Steven M.; Machicoane, Guillaume; Wu, Xiaoyu; Morgan, Glenn
2016-02-01
To meet the beam power requirements of 400 kW at the fragmentation target for facility for Rare Isotope Beams (FRIB), simultaneous acceleration of two-charge states should be used for heavier ions. These intense multi-charged ion beams will be produced by a 28 GHz electron cyclotron resonance (ECR) ion source at a high voltage of 35 kV. After extraction, the ion beam will be pre-accelerated to 12 keV/u with a 50 kV platform, transported down to an achromatic charge state selection (CSS) system followed by a vertical transport line, and then injected into a radio frequency quadrupole accelerator. The TRACK code developed at ANL is used to perform the simulations of the ECR beam extraction and low energy beam transport for FRIB. In this study, we include the magnetic field of ECR ion source into simulations. Different initial beam conditions as well as different space charge neutralization levels are tested for the ECR beamline. The beam loss in CSS system and the corresponding protective measures are discussed. The detailed results about the beam dynamic simulation and beam loss in CSS system will be presented in this paper.
Beam simulation studies of ECR beam extraction and low energy beam transport for FRIB
To meet the beam power requirements of 400 kW at the fragmentation target for facility for Rare Isotope Beams (FRIB), simultaneous acceleration of two-charge states should be used for heavier ions. These intense multi-charged ion beams will be produced by a 28 GHz electron cyclotron resonance (ECR) ion source at a high voltage of 35 kV. After extraction, the ion beam will be pre-accelerated to 12 keV/u with a 50 kV platform, transported down to an achromatic charge state selection (CSS) system followed by a vertical transport line, and then injected into a radio frequency quadrupole accelerator. The TRACK code developed at ANL is used to perform the simulations of the ECR beam extraction and low energy beam transport for FRIB. In this study, we include the magnetic field of ECR ion source into simulations. Different initial beam conditions as well as different space charge neutralization levels are tested for the ECR beamline. The beam loss in CSS system and the corresponding protective measures are discussed. The detailed results about the beam dynamic simulation and beam loss in CSS system will be presented in this paper
Electron beam linac simulations, 1985
The results of particle simulation studies of electron beam problems for the year 1985 are summarized. Although the main subjects this year were mostly the same as in previous years, namely foilless diodes, accelerating gaps, and beam extraction, a major new element this year was the addition of IFR (ion-focused-regime) channels to the various components of a linac beam line. Of particular interest may be the new results obtained for IFR autoacceleration gaps (relevant to MIMI), and the studies of extraction from a guide magnetic field onto an IFR propagation channel (relevant to RADLAC II). 49 figs., 20 refs
RIA Beam Dynamics Comparing TRACK to IMPACT
Mustapha, Brahim; Ostroumov, Peter; Qiang, Ji; Ryne, Robert D
2005-01-01
In order to benchmark the newly developed beam dynamics code TRACK we have performed comparisons with well established existing codes. During code development, codes like TRANSPORT, COSY, GIOS and RAYTRACE were used to check TRACK's implementation of the different beam line elements. To benchmark the end-to-end simulation of the RIA driver linac, the simulation of the low-energy part (from the ion source to the entrance of the SC linac) was compared with PARMTEQ and found to agree well. For the simulation of the SC linac the code IMPACT is used. Prior to these simulations, the code IMPACT had to be updated to meet the special requirements of the RIA driver linac. Features such as multiple charge state acceleration, stripper simulation and beam collimation were added to the code. IMPACT was also modified to support new types of rf cavities and to include fringe fields for all the elements. This paper will present a comparison of the beam dynamics simulation in the RIA driver linac between the codes TRACK and I...
A contemporary guide to beam dynamics
A methodological discussion is given for single particle beam dynamics in circular machines. The discussions are introductory, but (or, even therefore) we avoid to rely on too much simplified concepts. We treat things from a very general and fundamental point of view, because this is the easiest and rightest way to teach how to simulate particle motion and how to analyze its results. We give some principles of particle tracking free from theoretical prejudices. We also introduce some transparent methods to deduce the necessary information from the tracking: many of the traditional beam-dynamics concepts can be abstracted from them as approximate quantities which are valid in certain limiting cases
A contemporary guide to beam dynamics
A methodological discussion is given for single particle beam dynamics in circular machines. The discussions are introductory, but (or, even therefore) we avoid to rely on too much simplified concepts. We treat things from a very general and fundamental point of view, because this is the easiest and rightest way to teach how to simulate particle motion and how to analyze its results. We give some principles of particle tracking free from theoretical prejudices. We also introduce some transparent methods to deduce the necessary information from the tracking: many of the traditional beam-dynamics concepts can be abstracted from them as approximate quantities which are valid in certain limiting cases. (author)
We study the effects of small-scale parameter on the buckling loads and strains of nanobeams, based on nonlocal Timoshenko beam model. However, the lack of higher order boundary conditions leads to inconsistencies in critical buckling loads. In this paper, we apply a novel approach based on nonlocal Timoshenko kinematics, strain gradient approach and variational methods for deriving all classical and higher-order boundary conditions as well as governing equations. Therefore, closed-form and exact critical buckling loads of nanobeams with various end conditions are investigated. Moreover, the dependence of buckling loads on the small-scale parameter as well as shear deformation coefficient is studied using these new boundary conditions. Then, numerical results from this new beam model are presented for carbon nanotubes (CNTs). They illustrate a more accurate buckling response as compared to the previous works. Furthermore, the critical strains are compared with results obtained from molecular dynamic simulations as well as Sanders shell theory and are found to be in good agreement. Results show that unlike the other beam theories, this model can capture correctly the small-scale effects on buckling strains of short CNTs for the shell-type buckling. Moreover, the value of nonlocal constant is calculated for CNTs using molecular dynamic simulation results. (author)
Beam Simulations for IRE and Driver-Status and Strategy
The methods and codes employed in the U.S. Heavy Ion Fusion program to simulate the beams in an Integrated Research Experiments (IRE) facility and a fusion driver are presented in overview. A new family of models incorporating accelerating module impedance, multi-beam, and self-magnetic effects is described, and initial WARP3d particle simulations of beams using these models are presented. Finally, plans for streamlining the machine-design simulation sequence, and for simulating beam dynamics from the source to the target in a consistent and comprehensive manner, are described
Beam and spin dynamics of hadron beams in intermediate-energy ring accelerators
In this thesis beam and spin dynamics of ring accelerators are described. After a general theoretical treatment methods for the beam optimization and polarization conservation are discussed. Then experiments on spin manipulation at the COSY facility are considered. Finally the beam simulation and accelerator lay-out for the HESR with regards to the FAIR experiment are described. (HSI)
Dynamic Bowtie for Fan-beam CT
Liu, Fenglin; Cong, Wenxiang; Hsieh, Scott; Pelc, Norbert
2013-01-01
A bowtie is a filter used to shape an x-ray beam and equalize its flux reaching different detector channels. For development of spectral CT with energy-discriminative photon-counting (EDPC) detectors, here we propose and evaluate a dynamic bowtie for performance optimization based on a patient model or a scout scan. Our dynamic bowtie modifies an x-ray beam intensity profile by mechanical rotation and adaptive adjustment of the x-ray source flux. First, a mathematical model for dynamic bowtie filtering is established for an elliptical section in fan-beam geometry, and the contour of the optimal bowtie is derived. Then, numerical simulation is performed to compare the performance of the dynamic bowtie in the cases of an ideal phantom and a realistic cross-section relative to the counterparts without any bowtie and with a fixed bowtie respectively. Our dynamic bowtie can equalize the expected numbers of photons in the case of an ideal phantom. In practical cases, our dynamic bowtie can effectively reduce the dy...
Beam dynamics in CIME for third harmonic
This report presents the results from simulations for beam dynamics in CIME third harmonics. Details are given regarding the procedures to reach the adaptation at the inflector exit. The aim of these simulation is to determine, for any given ion, the beam correlations at the inflector exit as well as the current values in the isochronous coils for all the field levels. Although not all the steps of the simulation are thoroughly displayed, the report gathers all the the elements necessary for CIME control. Information useful for controlling the Very Low Energy line, the main field and the isochronous coils are also presented. The report has the following content: I. Introduction. II. The field maps and the used codes. A. The maps of CIME magnetic fields; B. The 3D map of CIME electric potentials; C. The maps of 3D electric potentials in the CIME central region; D. Code LIONS and sorting codes. III. Central region. A. An outlook. B.Central rays; IV. Determination of beam correlations. A. Analytical calculation of adaptation conditions; B. Calculation of adaptation conditions based on particle distributions; C. Creating the beam matrices. D. Calculation method for inverse return correlations. V. Results of simulations. VI. Interpolation of isochronous coils at the referential frequency. VII. The interpolation code PARAM. VIII. Conclusions. The paper is supplemented by 4 appendices. The harmonics 2, 4 and 5 are currently under way and the results will be reported in a future paper
Simulating transient effects of pulsed beams on beam intercepting devices
The development in the physics community towards higher beam power through the possibilities of particle accelerators lead to challenges for the developers of elements which are exposed to effect of particle beams (beam intercepting devices =BIDs). For the design of BIDs, the increasing heat load onto these devices due to energetic and focused beams and - in most cases - their highly pulsed nature has to be taken into account. The physics requirements are sometimes opposed to the current state of the art. As one possibility of many in combining the different aspects for these ambitious demands, two highly developed computer programs, namely FLUKA and ANSYS AUTODYN, were joined for this dissertation. The former is a widely enhanced Monte-Carlo-code which specializes on the interaction of particles with static matter, while the latter is a versatile explicit code for the simulation of highly dynamic processes. Both computer programs were developed intensively over years and are still continuously enhanced in order to achieve their best potential. As a consequence of their separate development histories, their combination requires a large amount of work - at the physics limits of their application as well as at the frontier of computing technology. The current work did not touch all different points needed for a full integration, but it is a first step towards their coupling within a feasible time frame. For the simulation of metallic targets irradiated with highly energetic uranium ions different material models have been combined and one parameter describing the damage of the material was varied. In the case of two copper targets, this procedure led to a qualitative agreement between simulations and experimental results. (author)
Salau Tajudeen A.O.
2014-01-01
Full Text Available This study reported a simulation approach to the understanding of the interactions between a buried pipe and the soil system by computing the bending stress variation of harmonically-excited buried pipes. The established principles of linear dynamics theory and simple beam theory were utilised in the analysis of the problem of buried pipe bending stress accumulation and its dynamics. With regards to the parameters that influence the bending stress variations, the most important are the isolation factor, uniform external load, and the corresponding limiting conditions. The simulated mathematical expressions, containing static and dynamic parameters of the buried pipe and earth, were coded in Fortran programming language and applied in the simulation experiment. The results obtained showed that harmonically-excited buried thick-walled pipe became stable and effective when the ratio of the natural frequency of vibration to the forced frequency is greater than 2.0, whenever the damped factor is used as the control parameter for the maximum bending stress. The mirror image of the stress variation produces variation in the location of the maximum bending stress in quantitative terms. The acceptable pipe materials for the simulated cases must have yield strength in bending greater than or equal to 13.95 MPa. The results obtained in this work fill a gap in the literature and will be useful to pipeline engineers and designers, as well as to environmental scientists in initialising and controlling environmental issues and policy formulation concerning the influence of buried pipe on the soil and water in the environment.
Beam dynamics studies for photocathode RF gun
Photocathode RF guns are very popular choice as injector for low emittance beams especially to light sources world wide. In demand for these gun is increasing steadily and efforts are on to make 2.6 cell RF Gun as SAMEER as proto type for future use at various laboratories. The base design of this 2.6 cell RF Gun is ready and fabrication is planned in near future. In this paper, we present beam dynamic study results of the gun and methodology to arrive at the operating point. Simulation results for Gaussian with nano-second pulse length will be discussed in detail and proposal for generation of few MeV beam will be presented. (author)
Irradiation of nanotube structures with electron and ion beams has been used to produce functionalized nanotubes and fundamentally new structures, including junctions. Here, we build on previous studies to investigate the low-energy electron and ion (Ar and CF3) beam irradiation of triple-walled carbon nanotubes that consist entirely of either chiral or armchair tubes. Effective incident energies of 50 eV/ion and 50 keV/electron are considered. The approach is classical molecular dynamic simulations using reactive empirical bond-order potentials and the primary knock-on atom approach to model the effects of electron irradiation. The results indicate that these various irradiation processes produce local damage to the nanotubes that includes crosslinking, that the degree of damage depends to some degree on the chirality of the nanotubes, and that the radial distribution of crosslinks depends significantly on the irradiating particle. Importantly, the effect of these crosslinks, and their radial distribution along the circumference of the nanotube, on the tendency of multiwalled nanotubes to fail by the sword-in-sheath mechanism is examined
Beam dynamics issues for linear colliders
Ruth, R.D.
1987-09-01
In this paper we discuss various beam dynamics issues for linear colliders. The emphasis is to explore beam dynamics effects which lead to an effective dilution of the emittance of the beam and thus to a loss of luminosity. These considerations lead to various tolerances which are evaluated for a particular parameter set.
Dynamics of beam halo in mismatched beams
High-power proton linacs for nuclear materials transmutation and production, and new accelerator-driven neutron spallation sources must be designed to control beam-halo formation, which leads to beam loss. The study of particle-core models is leading to a better understanding of the causes and characteristics of beam halo produced by space-charge forces in rms mismatched beams. Detailed studies of the models have resulted in predictions of the dependence of the maximum amplitude of halo particles on a mismatch parameter and on the space-charge tune-depression ratio. Scaling formulas have been derived which will provide guidance for choosing the aperture radius to contain the halo without loss. (author)
Beam stability and nonlinear dynamics. Summary report
A open-quotes Beam Stability and Nonlinear Dynamicsclose quotes Symposium was held October 3-5, 1996 at the Institute for Theoretical Physics (ITP) in Santa Barbara. This was one of the 3 symposia hosted by the ITP and supported by its sponsor, the National Science Foundation, as part of our open-quotes New Ideas for Particle Acceleratorsclose quotes program. The symposia was organized and chaired by Dr. Zohreh Parsa of ITP/ Brookhaven National Laboratory. The purpose of this symposium was to deal with some of the fundamental theoretical problems of accelerator physics by bringing together leaders from accelerator physics communities, mathematics, and other fields of physics. The focus was on nonlinear dynamics and beam stability. The symposium began with some defining talks on relevant mathematical topics such as single-particle Hamiltonian dynamics, chaos, and new ideas in symplectic integrators. The physics topics included single-particle and many-particle dynamics. These topics concern circular accelerators in which particles circulate for a very large number of turns as well as linear accelerators where space charge and wakefields induced in accelerating cavities play a strong role. A major question is to determine the best model for numerical simulations in order to accurately reproduce behavior of beams in real accelerators and to predict long-term or long distance stability. Comparison with experiment is recognized as an important tool in improving models
Controllable circular Airy beams via dynamic linear potential
Zhong, Hua; Belić, Milivoj R; Li, Changbiao; Wen, Feng; Zhang, Zhaoyang; Zhang, Yanpeng
2016-01-01
We investigate controllable spatial modulation of circular autofocusing Airy beams, under action of different dynamic linear potentials, both theoretically and numerically. We introduce a novel treatment method in which the circular Airy beam is represented as a superposition of narrow azimuthally-modulated one-dimensional Airy beams that can be analytically treated. The dynamic linear potentials are appropriately designed, so that the autofocusing effect can either be weakened or even eliminated when the linear potential exerts a "pulling" effect on the beam, or if the linear potential exerts a "pushing" effect, the autofocusing effect can be greatly strengthened. Numerical simulations agree with the theoretical results very well.
Spacecraft Dynamic Characteristics While Deploying Flexible Beams
程绪铎; 李俊峰; 樊勇; 王照林
2002-01-01
The attitude dynamic equations of a spacecraft while deploying two flexible beams and the beam equations were developed from momentum theory. The dynamic equations were solved numerically using the Runge-Kutta method to calculate the vibration amplitudes of the flexible beams and the attitude angular velocity. The results show that the vibration amplitudes increase as the beam length increases or as the initial attitude angular velocity increases. The results also show that the vibration amplitudes decrease as the deployment velocity increases.
Large deformation dynamic bending of composite beams
Derian, Edward J.
1985-01-01
The large deformation response of composite beams subjected to a dynamic axial load was studied. The beams were loaded with a moderate amount of eccentricity to promote bending. The study was primarily experimental but some finite element results were obtained. Both the deformation and the failure of the beams were of interest. The static response of the beams was also studied in order to determine the difference between the static and dynamic failure. Twelve different la...
Simulations of beam-beam and beam-wire interactions in RHIC
Kim, Hyung J.; Sen, Tanaji; /Fermilab; Abreu, Natalia P.; Fischer, Wolfram; /Brookhaven
2009-02-01
The beam-beam interaction is one of the dominant sources of emittance growth and luminosity lifetime deterioration. A current carrying wire has been proposed to compensate long-range beam-beam effects in the LHC and strong localized long-range beam-beam effects are experimentally investigated in the RHIC collider. Tune shift, beam transfer function, and beam loss rate are measured in dedicated experiments. In this paper, they report on simulations to study the effect of beam-wire interactions based on diffusive apertures, beam loss rates, and beam transfer function using a parallelized weak-strong beam simulation code (BBSIMC). The simulation results are compared with measurements performed in RHIC during 2007 and 2008.
Müller, Kei W.; Meier, Christoph; Wall, Wolfgang A.
2015-12-01
Networks of crosslinked biopolymer filaments such as the cytoskeleton are the subject of intense research. Oftentimes, mechanics on the scale of single monomers (∼ 5 nm) govern the mechanics of the entire network (∼ 10 μm). Until now, one either resolved the small scales and lost the big (network) picture or focused on mechanics above the single-filament scale and neglected the molecular architecture. Therefore, the study of network mechanics influenced by the entire spectrum of relevant length scales has been infeasible so far. We propose a method that reconciles both small and large length scales without the otherwise inevitable loss in either numerical efficiency or geometrical (molecular) detail. Both explicitly modeled species, filaments and their crosslinkers, are discretized with geometrically exact beam finite elements of Simo-Reissner type. Through specific coupling conditions between the elements of the two species, mechanical joints can be established anywhere along a beam's centerline, enabling arbitrary densities of chemical binding sites. These binding sites can be oriented to model the monomeric architecture of polymers. First, we carefully discuss the method and then demonstrate its capabilities by means of a series of numerical examples.
Molecular dynamics simulations
Tarmyshov, Konstantin B.
2007-01-01
Molecular simulations can provide a detailed picture of a desired chemical, physical, or biological process. It has been developed over last 50 years and is being used now to solve a large variety of problems in many different fields. In particular, quantum calculations are very helpful to study small systems at a high resolution where electronic structure of compounds is accounted for. Molecular dynamics simulations, in turn, are employed to study development of a certain molecular ensemble ...
Statics and rotational dynamics of composite beams
Ghorashi, Mehrdaad
2016-01-01
This book presents a comprehensive study of the nonlinear statics and dynamics of composite beams and consists of solutions with and without active elements embedded in the beams. The static solution provides the initial conditions for the dynamic analysis. The dynamic problems considered include the analyses of clamped (hingeless) and articulated (hinged) accelerating rotating beams. Two independent numerical solutions for the steady state and the transient responses are presented. The author illustrates that the transient solution of the nonlinear formulation of accelerating rotating beam converges to the steady state solution obtained by the shooting method. Other key areas considered include calculation of the effect of perturbing the steady state solution, coupled nonlinear flap-lag dynamics of a rotating articulated beam with hinge offset and aerodynamic damping, and static and dynamic responses of nonlinear composite beams with embedded anisotropic piezo-composite actuators. The book is intended as a t...
Interactive Dynamic-System Simulation
Korn, Granino A
2010-01-01
Showing you how to use personal computers for modeling and simulation, Interactive Dynamic-System Simulation, Second Edition provides a practical tutorial on interactive dynamic-system modeling and simulation. It discusses how to effectively simulate dynamical systems, such as aerospace vehicles, power plants, chemical processes, control systems, and physiological systems. Written by a pioneer in simulation, the book introduces dynamic-system models and explains how software for solving differential equations works. After demonstrating real simulation programs with simple examples, the author
Beam dynamics in magnetic quadrupole triplets
The Frankfurt Neutron source at the Stern-Gerlach-Zentrum (FRANZ) will produce high intensity neutron pulses in the energy range of 1 to 500 keV at a very short repetition rate. The neutrons are gained from 7Li(p,n)7Be reactions induced by 2 MeV protons and will be used to examine the nucleosynthesis during the s-process as it occurs in stars, cross sections of neutron capture reactions as well as the behaviour of non-neutral plasmas. In the linear accelerator section, consisting of a 4-rod-radio-frequency-quadrupole and a H-type drift tube linac, the proton pulses are accelerated to 2.03 MeV. Inside the drift tube cavity a magnetic quadrupole triplet will be integrated, in order to compensate transversally defocussing effects and therefore avoid losses. Behind the linear accelerator section the proton beam ist rebunched in a 5-cell CH-rebuncher which is framed by two more quadrupole triplets. To investigate the beam dynamics inside the magnetic quadrupole triplets, various magnetostatic and particle tracking codes like CST Studio and LORASR were used to simulate the beam transport properties of the magnets and compare the individual magnetic field distributions with the ones measured at the magnet laboratory at GSI. In doing so, important aspects to be considered are the longitudinal and transversal fringe fields and the saturation effects which all possibly cause emittance growth and geometrical aberrations.
Beam dynamics in magnetic quadrupole triplets
Claessens, Christine; Heilmann, Manuel; Meusel, Oliver; Podlech, Holger; Ratzinger, Ulrich; Wiesner, Christoph [IAP, Frankfurt University, Frankfurt am Main (Germany)
2013-07-01
The Frankfurt Neutron source at the Stern-Gerlach-Zentrum (FRANZ) will produce high intensity neutron pulses in the energy range of 1 to 500 keV at a very short repetition rate. The neutrons are gained from {sup 7}Li(p,n){sup 7}Be reactions induced by 2 MeV protons and will be used to examine the nucleosynthesis during the s-process as it occurs in stars, cross sections of neutron capture reactions as well as the behaviour of non-neutral plasmas. In the linear accelerator section, consisting of a 4-rod-radio-frequency-quadrupole and a H-type drift tube linac, the proton pulses are accelerated to 2.03 MeV. Inside the drift tube cavity a magnetic quadrupole triplet will be integrated, in order to compensate transversally defocussing effects and therefore avoid losses. Behind the linear accelerator section the proton beam ist rebunched in a 5-cell CH-rebuncher which is framed by two more quadrupole triplets. To investigate the beam dynamics inside the magnetic quadrupole triplets, various magnetostatic and particle tracking codes like CST Studio and LORASR were used to simulate the beam transport properties of the magnets and compare the individual magnetic field distributions with the ones measured at the magnet laboratory at GSI. In doing so, important aspects to be considered are the longitudinal and transversal fringe fields and the saturation effects which all possibly cause emittance growth and geometrical aberrations.
S C L Srivastava; S V L S Rao; P Singh
2007-10-01
A code for 2D space-charge dominated beam dynamics study in beam transport lines is developed. The code is used for particle-in-cell (PIC) simulation of -uniform beam in a channel containing solenoids and drift space. It can also simulate a transport line where quadrupoles are used for focusing the beam. Numerical techniques as well as the results of beam dynamics studies are presented in the paper.
Computer simulation of surface modification with ion beams
Niobium surface modification dynamics treated by cluster ion irradiation was studied based on atomistic and mesoscopic simulation methods and the results were compared to experiments. A surface smoothening method was proposed consisting of a treatment of the Nb cavity surfaces by accelerated gas (argon) cluster ion beams (GCIB) that is capable of reducing the surface roughness up to the theoretical limit
Beam dynamics design of the Compact Linear Collider Drive Beam injector
In the Compact Linear Collider (CLIC) the RF power for the acceleration of the Main Beam is extracted from a high-current Drive Beam that runs parallel to the main linac. The longitudinal and transverse beam dynamics of the Drive Beam injector has been studied in detail and optimized. The injector consists of a thermionic gun followed by a bunching system, some accelerating structures, and a magnetic chicane. The bunching system contains three sub-harmonic bunchers, a prebuncher, and a traveling wave buncher all embedded in a solenoidal magnetic field. The main characteristic of the Drive Beam injector is the phase coding process done by the sub-harmonic bunching system operating at half the acceleration frequency. This process is essential for the frequency multiplication of the Drive Beam. During the phase coding process the unwanted satellite bunches are produced that adversely affects the machine power efficiency. The main challenge is to reduce the population of particles in the satellite bunches in the presence of strong space-charge forces due to the high beam current. The simulation of the beam dynamics has been carried out with PARMELA with the goal of optimizing the injector performance compared to the existing model studied for the Conceptual Design Report (CDR). The emphasis of the optimization was on decreasing the satellite population, the beam loss in the magnetic chicane and limiting the beam emittance growth in transverse plane
Beam dynamics design of the Compact Linear Collider Drive Beam injector
Hajari, Sh. Sanaye, E-mail: ssanayeh@cern.ch [Institute for Research in Fundamental Sciences (IPM), School of Particles and Accelerators, P.O. Box 19395-5531, Tehran (Iran, Islamic Republic of); European Organization for Nuclear Research (CERN), BE Department, CH-1211 Geneva 23 (Switzerland); Shaker, H. [Institute for Research in Fundamental Sciences (IPM), School of Particles and Accelerators, P.O. Box 19395-5531, Tehran (Iran, Islamic Republic of); European Organization for Nuclear Research (CERN), BE Department, CH-1211 Geneva 23 (Switzerland); Doebert, S. [European Organization for Nuclear Research (CERN), BE Department, CH-1211 Geneva 23 (Switzerland)
2015-11-01
In the Compact Linear Collider (CLIC) the RF power for the acceleration of the Main Beam is extracted from a high-current Drive Beam that runs parallel to the main linac. The longitudinal and transverse beam dynamics of the Drive Beam injector has been studied in detail and optimized. The injector consists of a thermionic gun followed by a bunching system, some accelerating structures, and a magnetic chicane. The bunching system contains three sub-harmonic bunchers, a prebuncher, and a traveling wave buncher all embedded in a solenoidal magnetic field. The main characteristic of the Drive Beam injector is the phase coding process done by the sub-harmonic bunching system operating at half the acceleration frequency. This process is essential for the frequency multiplication of the Drive Beam. During the phase coding process the unwanted satellite bunches are produced that adversely affects the machine power efficiency. The main challenge is to reduce the population of particles in the satellite bunches in the presence of strong space-charge forces due to the high beam current. The simulation of the beam dynamics has been carried out with PARMELA with the goal of optimizing the injector performance compared to the existing model studied for the Conceptual Design Report (CDR). The emphasis of the optimization was on decreasing the satellite population, the beam loss in the magnetic chicane and limiting the beam emittance growth in transverse plane.
Screw theoretic view on dynamics of spatially compliant beam
Xi-lun DING; J.M.SELIG
2010-01-01
Beams with spatial compliance can be deformed as bending in a plane,twisting,and extending.In terms of the screw theory on rigid body motions,the concept of"deflection screw"is introduced,a spatial compliant beam theory via the deflection screw is proposed,and the spatial compliance of such a beam system is presented and analysed based on the material theory and fundamental kinematic assumptions.To study the dynamics of the spatially compliant beam,the potential energy and the kinetic energy of the beam are discussed by using the screw theory to obtain the Lagrangian.The Rayleigh-Ritz method is used to compute the vibrational frequencies based on discussions of boundary conditions and shape functions.The eigenfrequencies of the beam with spatial compliance are compared with those of individual deformation cases,pure bending,extension,or torsion.Finally,dynamics of a robot with two spatial compliant links and perpendicular joints is studied using the spatial compliant beam theory.Coupling between the joint rigid body motions and the deformations of spatial compliant links can easily be found in dynamic simulation.The study shows the effectiveness of using the screw theory to deal with the problems of dynamic modeling and analysis of mechanisms with spatially compliant links.
W-Band Sheet Beam Klystron Simulation
Colby, E.R.; Caryotakis, G.; Fowkes, W.R.; /SLAC; Smithe, D.N.; /Mission Res., Newington
2005-09-12
With the development of ever higher energy particle accelerators comes the need for compactness and high gradient, which in turn require very high frequency high power rf sources. Recent development work in W-band accelerating techniques has spurred the development of a high-power W-band source. Axisymmetric sources suffer from fundamental power output limitations (P{sub sat} {approx} {lambda}{sup 2}) brought on by the conflicting requirements of small beam sizes and high beam current. The sheet beam klystron allows for an increase in beam current without substantial increase in the beam current density, allowing for reduced cathode current densities and focusing field strengths. Initial simulations of a 20:1 aspect ratio sheet beam/cavity interaction using the 3 dimensional particle-in-cell code Magic3D have demonstrated a 35% beam-power to RF power extraction efficiency. Calculational work and numerical simulations leading to a prototype W-band sheet beam klystron will be presented, together with preliminary cold test structure studies of a proposed RF cavity geometry.
Beam Simulation Studies of the LEBT for RIA Driver Linac
The low energy beam transport (LEBT) system in the front-end of the Rare Isotope Accelerator (RIA) uses a 70 kV platform to pre-accelerate the ion beam from a 30 kV Electron Cyclotron Resonance (ECR) ion source, followed by an achromatic charge selection system. The selected beam is then pre-bunched and matched into the entrance of a Radio Frequency Quadrupole (RFQ) with a multi-harmonic buncher. To meet the beam power requirements for heavy ions, high current (several mA), multi-species beams will be extracted from the ECR. Therefore, it is crucial to control space charge effects in order to obtain the low emittance beam required for RIA. The PARMELA code is used to perform the LEBT simulations for the multi-species beams with 3D space charge calculations. The results of the beam dynamics simulations are presented, and the key issues of emittance growth in the LEBT and its possible compensation are discussed
Quantum fluctuations in beam dynamics.
Kim, K.-J.
1998-06-04
Quantum effects could become important for particle and photon beams used in high-luminosity and high brightness applications in the current and next generation accelerators and radiation sources. This paper is a review of some of these effects.
Quantum fluctuations in beam dynamics
Quantum effects could become important for particle and photon beams used in high-luminosity and high brightness applications in the current and next generation accelerators and radiation sources. This paper is a review of some of these effects
Dynamics of a high-current relativistic electron beam
The dynamics of a high-current relativistic electron beam is studied experimentally and by numerical simulation. The beam is formed in a magnetically insulated diode with a transverse-blade explosive-emission cathode. It is found experimentally that the radius of a 500-keV beam with a current of 2 kA and duration of 500 ns decreases with time during the beam current pulse. The same effect was observed in numerical simulations. This effect is explained by a change in the shape of the cathode plasma during the current pulse, which, according to calculations, leads to a change in the beam parameters, such as the electron pitch angle and the spread over the longitudinal electron momentum. These parameters are hard to measure experimentally; however, the time evolution of the radial profile of the beam current density, which can be measured reliably, coincides with the simulation results. This allows one to expect that the behavior of the other beam parameters also agrees with numerical simulations
Beam dynamics in high energy particle accelerators
Wolski, Andrzej
2014-01-01
Particle accelerators are essential tools for scientific research in fields as diverse as high energy physics, materials science and structural biology. They are also widely used in industry and medicine. Producing the optimum design and achieving the best performance for an accelerator depends on a detailed understanding of many (often complex and sometimes subtle) effects that determine the properties and behavior of the particle beam. Beam Dynamics in High Energy Particle Accelerators provides an introduction to the concepts underlying accelerator beam line design and analysis, taking an approach that emphasizes the elegance of the subject and leads into the development of a range of powerful techniques for understanding and modeling charged particle beams.
Simulation studies for ion beam extraction systems
Abdelrahman, M.M.; Zakhary, S.G. [Atomic Energy, Cairo (Egypt). Nuclear Research Center. Accelerators and Ion Sources Dept.], e-mail: moustafa82003@yahoo.com
2009-06-15
The characteristics of the ion beam extracted from an ion sources were investigated using computer code SIMION 3 D Version 7.0. It has been used to evaluate the extraction system in order to produce an ion beam with high current and low beam emittance. The results show that the shape of the extraction electrode plays an important role in ion beam formation. Comparison has been made between two extraction systems, Pierce extraction electrode and spherical extraction electrode. The results show that the spherical extraction system yields ion extraction beam with lower emittance and radius than that the Pierce system. The simulation can provide the basis for optimizing the extraction system and the acceleration gap for ion source. (author)
The δf algorithm for beam dynamics
An algorithm is developed to study particle dynamics of beams including collective interaction with high accuracy and low noise. Particle dynamics with collective interactions is treated through particle simulation, where the main or average distribution f0 and the deviation away from it δf are separately followed. The main distribution f0 is handled by an analytic equilibrium solution and the perturbation away from it δf is followed by the method of characteristics. We call this the δf algorithm. We specifically model a synchrotron collider which includes the collision section where collective effects of collisions are simulated by this δf algorithm and the rest of the collider where single particle dynamics are treated by simple harmonic transport. The most important target of this simulation is to understand and predict the long-time behavior of the beam luminosity and lifetime. The δf method allows the study the effect of small perturbations over long timescales on beam lifetime by eliminating the numerical noise problem inherent in Particle-in-Cell techniques. In the δf code using the reference parameters of the SSC (Superconducting Super Collider), beam blow-up near resonances and oscillations in the tune shift, Δν, far from resonances are observed. In studying long timescale particle diffusion in the phase space of the beams away from resonances, the δf code performance is compared with a tracking code which does not incorporate collective interaction
Molecular dynamics simulations
The molecular dynamics computer simulation discovery of the slow decay of the velocity autocorrelation function in fluids is briefly reviewed in order to contrast that long time tail with those observed for the stress autocorrelation function in fluids and the velocity autocorrelation function in the Lorentz gas. For a non-localized particle in the Lorentz gas it is made plausible that even if it behaved quantum mechanically its long time tail would be the same as the classical one. The generalization of Fick's law for diffusion for the Lorentz gas, necessary to avoid divergences due to the slow decay of correlations, is presented. For fluids, that generalization has not yet been established, but the region of validity of generalized hydrodynamics is discussed. 20 refs., 5 figs
Photon collider beam simulation with CAIN
Aleksander Filip Żarnecki
2007-11-01
The CAIN simulation program was used to study the outgoing beam profile for the photon collider at ILC. The main aim of the analysis was to verify the feasibility of the photon linear collider running with 20 mrad electron beam crossing angle. The main problem is the distorted electron beam, which has to be removed from the interaction region. It is shown that with a new design of the final dipole, it should be possible to avoid large energy losses at the face of the magnet.
Simulation of beam-induced plasma for the mitigation of beam-beam effects
Ma, J.; Wang, G.; Samulyak, R.; Yu, K.; Litvinenko, V.
2015-05-03
One of the main challenges in the increase of luminosity of circular colliders is the control of the beam-beam effect. In the process of exploring beam-beam mitigation methods using plasma, we evaluated the possibility of plasma generation via ionization of neutral gas by proton beams, and performed highly resolved simulations of the beam-plasma interaction using SPACE, a 3D electromagnetic particle-in-cell code. The process of plasma generation is modelled using experimentally measured cross-section coefficients and a plasma recombination model that takes into account the presence of neutral gas and beam-induced electromagnetic fields. Numerically simulated plasma oscillations are consistent with theoretical analysis. In the beam-plasma interaction process, high-density neutral gas reduces the mean free path of plasma electrons and their acceleration. A numerical model for the drift speed as a limit of plasma electron velocity was developed. Simulations demonstrate a significant reduction of the beam electric field in the presence of plasma. Preliminary simulations using fully-ionized plasma have also been performed and compared with the case of beam-induced plasma.
Simulating Electron Cloud Effects in Heavy-Ion Beams
Stray electrons can be introduced in heavy ion fusion accelerators as a result of ionization of ambient gas or gas released from walls due to halo-ion impact, or as a result of secondary-electron emission. We summarize here results from several studies of electron-cloud accumulation and effects: (1) Calculation of the electron cloud produced by electron desorption from computed beam ion loss; the importance of ion scattering is shown; (2) Simulation of the effect of specified electron cloud distributions on ion beam dynamics. We find electron cloud variations that are resonant with the breathing mode of the beam have the biggest impact on the beam (larger than other resonant and random variations), and that the ion beam is surprisingly robust, with an electron density several percent of the beam density required to produce significant beam degradation in a 200-quadrupole system. We identify a possible instability associated with desorption and resonance with the breathing mode. (3) Preliminary investigations of a long-timestep algorithm for electron dynamics in arbitrary magnetic fields
Radio Frequency Station - Beam Dynamics Interaction in Circular Accelerators
The longitudinal beam dynamics in circular accelerators is mainly defined by the interaction of the beam current with the accelerating Radio Frequency (RF) stations. For stable operation, Low Level RF (LLRF) feedback systems are employed to reduce coherent instabilities and regulate the accelerating voltage. The LLRF system design has implications for the dynamics and stability of the closed-loop RF systems as well as for the particle beam, and is very sensitive to the operating range of accelerator currents and energies. Stability of the RF loop and the beam are necessary conditions for reliable machine operation. This dissertation describes theoretical formalisms and models that determine the longitudinal beam dynamics based on the LLRF implementation, time domain simulations that capture the dynamic behavior of the RF station-beam interaction, and measurements from the Positron-Electron Project (PEP-II) and the Large Hadron Collider (LHC) that validate the models and simulations. These models and simulations are structured to capture the technical characteristics of the system (noise contributions, non-linear elements, and more). As such, they provide useful results and insight for the development and design of future LLRF feedback systems. They also provide the opportunity to study diverse longitudinal beam dynamics effects such as coupled-bunch impedance driven instabilities and single bunch longitudinal emittance growth. Coupled-bunch instabilities and RF station power were the performance limiting effects for PEP-II. The sensitivity of the instabilities to individual LLRF parameters, the effectiveness of alternative operational algorithms, and the possible tradeoffs between RF loop and beam stability were studied. New algorithms were implemented, with significant performance improvement leading to a world record current during the last PEP-II run of 3212 mA for the Low Energy Ring. Longitudinal beam emittance growth due to RF noise is a major concern for LHC
Radio Frequency Station - Beam Dynamics Interaction in Circular Accelerators
Mastoridis, Themistoklis [Stanford Univ., CA (United States)
2010-08-01
The longitudinal beam dynamics in circular accelerators is mainly defined by the interaction of the beam current with the accelerating Radio Frequency (RF) stations. For stable operation, Low Level RF (LLRF) feedback systems are employed to reduce coherent instabilities and regulate the accelerating voltage. The LLRF system design has implications for the dynamics and stability of the closed-loop RF systems as well as for the particle beam, and is very sensitive to the operating range of accelerator currents and energies. Stability of the RF loop and the beam are necessary conditions for reliable machine operation. This dissertation describes theoretical formalisms and models that determine the longitudinal beam dynamics based on the LLRF implementation, time domain simulations that capture the dynamic behavior of the RF station-beam interaction, and measurements from the Positron-Electron Project (PEP-II) and the Large Hadron Collider (LHC) that validate the models and simulations. These models and simulations are structured to capture the technical characteristics of the system (noise contributions, non-linear elements, and more). As such, they provide useful results and insight for the development and design of future LLRF feedback systems. They also provide the opportunity to study diverse longitudinal beam dynamics effects such as coupled-bunch impedance driven instabilities and single bunch longitudinal emittance growth. Coupled-bunch instabilities and RF station power were the performance limiting effects for PEP-II. The sensitivity of the instabilities to individual LLRF parameters, the effectiveness of alternative operational algorithms, and the possible tradeoffs between RF loop and beam stability were studied. New algorithms were implemented, with significant performance improvement leading to a world record current during the last PEP-II run of 3212 mA for the Low Energy Ring. Longitudinal beam emittance growth due to RF noise is a major concern for LHC
Radio Frequency Station - Beam Dynamics Interaction in Circular Accelerators
Mastoridis, Themistoklis; /Stanford U., Elect. Eng. Dept. /SLAC
2011-03-01
The longitudinal beam dynamics in circular accelerators is mainly defined by the interaction of the beam current with the accelerating Radio Frequency (RF) stations. For stable operation, Low Level RF (LLRF) feedback systems are employed to reduce coherent instabilities and regulate the accelerating voltage. The LLRF system design has implications for the dynamics and stability of the closed-loop RF systems as well as for the particle beam, and is very sensitive to the operating range of accelerator currents and energies. Stability of the RF loop and the beam are necessary conditions for reliable machine operation. This dissertation describes theoretical formalisms and models that determine the longitudinal beam dynamics based on the LLRF implementation, time domain simulations that capture the dynamic behavior of the RF station-beam interaction, and measurements from the Positron-Electron Project (PEP-II) and the Large Hadron Collider (LHC) that validate the models and simulations. These models and simulations are structured to capture the technical characteristics of the system (noise contributions, non-linear elements, and more). As such, they provide useful results and insight for the development and design of future LLRF feedback systems. They also provide the opportunity to study diverse longitudinal beam dynamics effects such as coupled-bunch impedance driven instabilities and single bunch longitudinal emittance growth. Coupled-bunch instabilities and RF station power were the performance limiting effects for PEP-II. The sensitivity of the instabilities to individual LLRF parameters, the effectiveness of alternative operational algorithms, and the possible tradeoffs between RF loop and beam stability were studied. New algorithms were implemented, with significant performance improvement leading to a world record current during the last PEP-II run of 3212 mA for the Low Energy Ring. Longitudinal beam emittance growth due to RF noise is a major concern for LHC
Nonlinear beam dynamics experimental program at SPEAR
Since nonlinear effects can impose strict performance limitations on modern colliders and storage rings, future performance improvements depend on further understanding of nonlinear beam dynamics. Experimental studies of nonlinear beam motion in three-dimensional space have begun in SPEAR using turn-by-turn transverse and longitudinal phase-space monitors. This paper presents preliminary results from an on-going experiment in SPEAR
Longitudinal beam dynamics for heavy ion fusion
The longitudinal wall impedance instability is of great interest for a heavy ion fusion (HIF) driver because complete stabilization of this mode via momentum spread is impractical due to requirements of focusing the beam onto the inertial confinement fusion target. This instability is being studied with the WARPrz particle-in-cell code. The impedance of the induction linac modules is modeled as a wall impedance corresponding to a continuum of resistors and capacitors in parallel. We discuss simulations of the this instability, including reflections of perturbations off the beam end and the effects of finite temperature, and simulations of errors in intermittently-applied axial confining fields as a seed for this instability. We also present very long simulations in which we look for beam equilibria
Recent progress in the simulation of heavy ion beams
Production of electric power by using a beam of heavy ions to ignite an inertially-confined fusion target requires the focusing of high-power beams onto a small spot several meters distant from the final lens system. Beams with the necessary intensity generally behave like warm nonneutral bounded plasmas where beam kinetic temperatures are sufficiently high that a cold-plasma description can be inadequate for describing the collective space-charge modes. In view of the complexity of the self-consistent nonlinear dynamics, analytic study has largely been limited to the singular Kapchinskij - Vladimirskij (K - V) distribution. Numerical simulations, primarily using the WARP [D. P. Grote, A. Friedman, I. Haber, W. Fawley, and J. L. Vay, Nucl. Instrum. Methods Phys. Res. A 415, 428 (1998)] particle-in-cell (PIC)/accelerator code, have been employed to identify the degree to which the analytic results, especially the predictions of unstable modes, are applicable to realistic beam distributions. During extensive benchmarking of the code against experiments at Lawrence Berkeley National Laboratory, Lawrence Livermore National Laboratory, and the University of Maryland, a particular feature of the beam which has been seen in both experiment and simulation is the launching, in the source region, of collective warm-plasma oscillations similar to those predicted on the basis of the K - V analysis. copyright 1999 American Institute of Physics
The Beam Break-Up Numerical Simulator
Travish, G.A.
1989-11-01
Beam Break-Up (BBU) is a severe constraint in accelerator design, limiting beam current and quality. The control of BBU has become the focus of much research in the design of the next generation collider, recirculating and linear induction accelerators and advanced accelerators. Determining the effect on BBU of modifications to cavities, the focusing elements or the beam is frequently beyond the ability of current analytic models. A computer code was written to address this problem. The Beam Break-Up Numerical Simulator (BBUNS) was designed to numerically solve for beam break-up (BBU) due to an arbitrary transverse wakefield. BBUNS was developed to be as user friendly as possible on the Cray computer series. The user is able to control all aspects of input and output by using a single command file. In addition, the wakefield is specified by the user and read in as a table. The program can model energy variations along and within the beam, focusing magnetic field profiles can be specified, and the graphical output can be tailored. In this note we discuss BBUNS, its structure and application. Included are detailed instructions, examples and a sample session of BBUNS. This program is available for distribution. 50 refs., 18 figs., 5 tabs.
The Beam Break-Up Numerical Simulator
Beam Break-Up (BBU) is a severe constraint in accelerator design, limiting beam current and quality. The control of BBU has become the focus of much research in the design of the next generation collider, recirculating and linear induction accelerators and advanced accelerators. Determining the effect on BBU of modifications to cavities, the focusing elements or the beam is frequently beyond the ability of current analytic models. A computer code was written to address this problem. The Beam Break-Up Numerical Simulator (BBUNS) was designed to numerically solve for beam break-up (BBU) due to an arbitrary transverse wakefield. BBUNS was developed to be as user friendly as possible on the Cray computer series. The user is able to control all aspects of input and output by using a single command file. In addition, the wakefield is specified by the user and read in as a table. The program can model energy variations along and within the beam, focusing magnetic field profiles can be specified, and the graphical output can be tailored. In this note we discuss BBUNS, its structure and application. Included are detailed instructions, examples and a sample session of BBUNS. This program is available for distribution. 50 refs., 18 figs., 5 tabs
Study of Effect of Ion Source Energy Spread on RFQ Beam Dynamics at REX-ISOLDE
Fraser, M A
2013-01-01
With an upgrade to the Electron Beam Ion Source (EBIS) at REX under consideration a study was launched in order to understand the effect of an increased energy spread from the ion source on the beam dynamics of the RFQ. Due to the increased electron beam potential needed to achieve the upgrade’s charge breeding specification it is expected that the energy spread of the beam will increase from today’s estimated value of approximately +-0.1%. It is shown through beam dynamics simulations that the energy spread can be increased to +-1% without significant degradation of the beam quality output by the RFQ.
Jacobi equations and particle accelerator beam dynamics
Torrome, Ricardo Gallego
2012-01-01
A geometric formulation of the linear beam dynamics in accelerator physics is presented. In particular, it is proved that the linear transverse and longitudinal dynamics can be interpret geometrically as an approximation to the Jacobi equation of an affine averaged Lorentz connection. We introduce a specific notion reference trajectory as integral curves of the main velocity vector field. A perturbation caused by the statistical nature of the bunch of particles is considered.
Muawia A. Magzoub
2014-07-01
Full Text Available This study presents the design of a Fuzzy Static (FS and a Fuzzy Dynamic (FD Sliding-Mode Controllers (SMC for both basic and complete ball on beam system. At first, the FSSMC was designed for the simplified and the complete models. Then, the FDSMC was designed on the simplified and the comprehensive models of the system in which the ball is placed on a beam as well. In addition, the lyapunov stability and linearization were used to check the stability of the system. There is an in-built issue of chattering with (FSSMC. However, (FDSMC counter it well. Also, FDSMC is effective with respect to matched disturbance rejection. It has been found out from this research study that the designs of the models which utilize a FDSMC with a comprehensive model of the system were more efficient than the designs that utilize the basic system’s prototype. Lastly, a comprehensive comparative analysis is provided and MATLAB/SIMULINK outcomes confirm the dominance of FDSMC.
Beam Delivery Simulation - Recent Developments and Optimization
AUTHOR|(INSPIRE)INSPIRE-00232566; Boogert, Stewart Takashi; Garcia-Morales, H; Gibson, Stephen; Kwee-Hinzmann, Regina; Nevay, Laurence James; Deacon, Lawrence Charles
2015-01-01
Beam Delivery Simulation (BDSIM) is a particle tracking code that simulates the passage of particles through both the magnetic accelerator lattice as well as their interaction with the material of the accelerator itself. The Geant4 toolkit is used to give a full range of physics processes needed to simulate both the interaction of primary particles and the production and subsequent propagation of secondaries. BDSIM has already been used to simulate linear accelerators such as the International Linear Collider (ILC) and the Compact Linear Collider (CLIC), but it has recently been adapted to simulate circular accelerators as well, producing loss maps for the Large Hadron Collider (LHC). In this paper the most recent developments, which extend BDSIM’s functionality as well as improve its efficiency are presented. Improvement and refactorisation of the tracking algorithms are presented alongside improved automatic geometry construction for increased particle tracking speed.
Guasp, J.; Fuentes, C.; Liniers, M.
2001-07-01
In this study the simulations of plasma transport under NBI for TJ-II, previously performed, are extended. Since than a considerable number of important modifications have been introduced in the model: change of magnetic configuration, use of experimental initial profiles, expansion of the Data base from NBI calculations and, mainly, a detailed handling of impurities with inclusion of sputtering effects. Moreover there is now a particular emphasis on the analysis of the conditions for discharge collapse and on the possible effects of single beam injection. This analysis of impurity behaviour with sputtering shows that in the expected usual cases there is no radioactive collapse and that if the recycling coefficients remain lower the unity it is always possible to find a strategy for external gas puffing leading to a stationary state, with densities below the limit and efficient NBI absorption (>50%). The radioactive collapse can appear either at high densities (central value higher than 1.4x10''20 m''3), excessive influx of impurities (i. e. with sputtering rates higher than twice the expected values) o for insufficient injected beam power (less than 45 kW). The present study analyses only the 100{sub 4}4{sub 6}4 configuration of TJ-II, but future works will start a systematic scan of configuration using this same model. (Author) 12 Refs.
Notes on beam dynamics in linear accelerators
Gluckstern, R.L.
1980-09-01
A collection of notes, on various aspects of beam dynamics in linear accelerators, which were produced by the author during five years (1975 to 1980) of consultation for the LASL Accelerator Technology (AT) Division and Medium-Energy Physics (MP) Division is presented.
Beam dynamic issues in TESLA damping ring
In this paper we study general requirements on impedances of the linear collider TESLA damping ring design. Quantitative consideration is performed for 17-km long ''dog-bone'' ring. Beam dynamics in alternative options of 6.3 and 2.3-km long damping rings is briefly discussed. 5 refs., 2 tabs
Beam dynamics studies in a tesla positron pre-accelerator
Moiseev, V A; Flöttmann, K
2001-01-01
The TESLA linear collider is based on superconducting accelerating cavities.Behind the positron production target normal conducting cavities have to be used in order to cope with high particle losses and with focusing solenoid surrounding the cavities.The main purpose of this pre-accelerator is to provide maximum capture efficiency for the useful part of the totally acceptable positron beam with technically reasonable parameters of the linac.The coupled optimization of the capture optics behind the target and pre-accelerator rf-operation has been carried out.The beam dynamics simulation results as well as the pre-accelerator peculiarities are presented.
Radiative transfer simulations of magnetar flare beaming
van Putten, T; Baring, M G; Wijers, R A M J
2016-01-01
Magnetar giant flares show oscillatory modulations in the tails of their light curves, which can only be explained via some form of beaming. The fireball model for magnetar bursts has been used successfully to fit the phase-averaged light curves of the tails of giant flares, but so far no attempts have been made to fit the pulsations. We present a relatively simple numerical model to simulate beaming of magnetar flare emission. In our simulations, radiation escapes from the base of a fireball trapped in a dipolar magnetic field, and is scattered through the optically thick magnetosphere of the magnetar until it escapes. Beaming is provided by the presence of a relativistic outflow, as well as by the geometry of the system. We find that a simple picture for the relativistic outflow is enough to create the pulse fraction and sharp peaks observed in pulse profiles of magnetar flares, while without a relativistic outflow the beaming is insufficient to explain giant flare rotational modulations.
Beam dynamics of alternating-phase-focused linac
Iwata, Y; Kapin, V
2004-01-01
A simple method to find an array of synchronous phases for alternating-phase-focused (APF) linacs is presented. The phase array is described with a smooth function having free parameters. With a set of the parameters, a simulation on the beam dynamics was made and distributions of the six-dimensional phase spaces were calculated for each set of the parameters. The parameters were varied, and numbers of the simulations have been performed. An optimum set of the parameters were determined so that the simulations of the beam dynamics yield large acceptances and small emittances of the extracted beams. Since the APF linac can provide both axial and radial stability of beams just with the rf acceleration-field, no additional focusing element inside of drift tubes are necessary. Comparing with conventional linacs having focusing elements, it has advantage in construction and operation costs as well as its acceleration rate. Therefore, the APF linacs would be suited for an injector of medical synchrotrons. A practic...
Introduction to Longitudinal Beam Dynamics
Holzer, B J
2014-01-01
This chapter gives an overview of the longitudinal dynamics of the particles in an accelerator and, closely related to that, the issue of synchronization between the particles and the accelerating field. Beginning with the trivial case of electrostatic accelerators, the synchronization condition is explained for a number of driven accelerators like Alvarez linacs, cyclotrons and finally synchrotrons and storage rings, where it plays a crucial role. In the case of the latter, the principle of phase focusing is motivated qualitatively as well as on a mathematically more correct level and the problem of operation below and above the transition energy is discussed. Throughout, the main emphasis is more on physical understanding rather than on a mathematically rigorous treatment.
An interactive beam position monitor system simulator
A system simulator has been implemented to aid the development of the RHIC position monitor system. Based on the LabVIEW software package by National Instruments, this simulator allows engineers and technicians to interactively explore the parameter space of a system during the design phase. Adjustable parameters are divided into three categories: beam, pickup, and electronics. The simulator uses these parameters in simple formulas to produce results in both time-domain and frequencydomain. During the prototyping phase, these simulated results can be compared to test data acquired with the same software package. The RHIC position monitor system is presented as an example, but the software is applicable to several other systems as well
Modelling and simulation of beam formation in electron guns
This paper describes a new PC version of the software package GUN-EBT for computer simulation of beam formation in rotationally symmetric electron guns with thermionic cathodes. It is based on a self-consistent physical model which takes into account the beam space charge and the initial velocity effects. The theoretical framework used for both the formulation of the model and for the interpretation of the results of numerical experiments is the formalism of the charged particle dynamics in phase space. This enables not only a trajectory analysis (ray tracing) but also a phase-space analysis of beams to be performed. The package can be used as an effective tool for computer aided design and optimization of electron guns in various electron-optical systems. The operation of the package is illustrated with a typical example. (orig.)
Stochastic beam dynamics in storage rings
In this thesis several approaches to stochastic dynamics in storage rings are investigated. In the first part the theory of stochastic differential equations and Fokker-Planck equations is used to describe the processes which have been assumed to be Markov processes. The mathematical theory of Markov processes is well known. Nevertheless, analytical solutions can be found only in special cases and numerical algorithms are required. Several numerical integration schemes for stochastic differential equations will therefore be tested in analytical solvable examples and then applied to examples from accelerator physics. In particular the stochastically perturbed synchrotron motion is treated. For the special case of a double rf system several perturbation theoretical methods for deriving the Fokker-Planck equation in the action variable are used and compared with numerical results. The second part is concerned with the dynamics of electron storage rings. Due to the synchrotron radiation the electron motion is influenced by damping and exciting forces. An algorithm for the computation of the density function in the phase space of such a dissipative stochastically excited system is introduced. The density function contains all information of a process, e.g. it determines the beam dimensions and the lifetime of a stored electron beam. The new algorithm consists in calculating a time propagator for the density function. By means of this propagator the time evolution of the density is modelled very computing time efficient. The method is applied to simple models of the beam-beam interaction (one-dimensional, round beams) and the results of the density calculations are compared with results obtained from multiparticle tracking. Furthermore some modifications of the algorithm are introduced to improve its efficiency concerning computing time and storage requirements. Finally, extensions to two-dimensional beam-beam models are described. (orig.)
Parallel beam dynamics calculations on high performance computers
Faced with a backlog of nuclear waste and weapons plutonium, as well as an ever-increasing public concern about safety and environmental issues associated with conventional nuclear reactors, many countries are studying new, accelerator-driven technologies that hold the promise of providing safe and effective solutions to these problems. Proposed projects include accelerator transmutation of waste (ATW), accelerator-based conversion of plutonium (ABC), accelerator-driven energy production (ADEP), and accelerator production of tritium (APT). Also, next-generation spallation neutron sources based on similar technology will play a major role in materials science and biological science research. The design of accelerators for these projects will require a major advance in numerical modeling capability. For example, beam dynamics simulations with approximately 100 million particles will be needed to ensure that extremely stringent beam loss requirements (less than a nanoampere per meter) can be met. Compared with typical present-day modeling using 10,000-100,000 particles, this represents an increase of 3-4 orders of magnitude. High performance computing (HPC) platforms make it possible to perform such large scale simulations, which require 10's of GBytes of memory. They also make it possible to perform smaller simulations in a matter of hours that would require months to run on a single processor workstation. This paper will describe how HPC platforms can be used to perform the numerically intensive beam dynamics simulations required for development of these new accelerator-driven technologies
Dynamic two-dimensional beam-pattern steering technique
Zhou, Shaomin; Yeh, Pochi; Liu, Hua-Kuang
1993-01-01
A dynamic two-dimensional laser-beam-pattern steering technique using photorefractive holograms in conjunction with electrically addressed spatial light modulators is proposed and investigated. The experimental results demonstrate the dynamic steering of random combinations of basis beam patterns. The proposed method has the advantages of random beam-pattern combination, good beam intensity uniformity, and higher diffraction efficiency compared with conventional methods.
Electron-beam dynamics for an advanced flash-radiography accelerator
Ekdahl, Carl August Jr. [Los Alamos National Laboratory
2015-06-22
Beam dynamics issues were assessed for a new linear induction electron accelerator. Special attention was paid to equilibrium beam transport, possible emittance growth, and beam stability. Especially problematic would be high-frequency beam instabilities that could blur individual radiographic source spots, low-frequency beam motion that could cause pulse-to-pulse spot displacement, and emittance growth that could enlarge the source spots. Beam physics issues were examined through theoretical analysis and computer simulations, including particle-in cell (PIC) codes. Beam instabilities investigated included beam breakup (BBU), image displacement, diocotron, parametric envelope, ion hose, and the resistive wall instability. Beam corkscrew motion and emittance growth from beam mismatch were also studied. It was concluded that a beam with radiographic quality equivalent to the present accelerators at Los Alamos will result if the same engineering standards and construction details are upheld.
Molecular dynamics simulation of diffusivity
Juanfang LIU; Danling ZENG; Qin LI; Hong GAO
2008-01-01
Equilibrium molecular dynamics simulation was performed on water to calculate its diffusivity by adopting different potential models. The results show that the potential models have great influence on the simulated results. In addition, the diffusivities obtained by the SPCE model conform well to the experimental values.
Investigations of the Dynamics of Space Charged Dominated Beams
York, Richard C.
2002-08-01
We propose to perform investigations of the dynamics of space charge dominated beams. These investigations will support present activities such as the electron ring project at the University of Maryland as well as provide an improved basis for future accelerator designs. Computer simulations will provide the primary research element with improved code development being an integral part of the activities during the first period. We believe that one of the code development projects provides a unique strategy for the inclusion of longitudinal dynamics, and that this concept should provide a computationally rapid research tool.
Beam Delivery Simulation: BDSIM - Development & Optimization
Nevay, Laurence James; Garcia-Morales, H; Gibson, S M; Kwee-Hinzmann, R; Snuverink, J; Deacon, L C
2014-01-01
Beam Delivery Simulation (BDSIM) is a Geant4 and C++ based particle tracking code that seamlessly tracks particles through accelerators and detectors, including the full range of particle interaction physics processes from Geant4. BDSIM has been successfully used to model beam loss and background conditions for many current and future linear accelerators such as the Accelerator Test Facility 2 (ATF2) and the International Linear Collider (ILC). Current developments extend its application for use with storage rings, in particular for the Large Hadron Collider (LHC) and the High Luminosity upgrade project (HL-LHC). This paper presents the latest results from using BDSIM to model the LHC as well as the developments underway to improve performance.
Beam dynamics calculations for fault-tolerance
The European Transmutation Demonstration requires a high-power proton accelerator operating in CW mode. This accelerator is also expected to have a very limited number of unexpected beam interruptions per year. To reach such an ambitious goal, it is clear that reliability-oriented design practices need to be followed from the early stage of components design and fault-tolerance capabilities have to be introduced to the maximum extent. The goal of this document is precisely to investigate in more details the fault-tolerance capability of the XT-ADS linac. From previous analysis, it appears that if nothing is done, a cavity's failure leads in nearly all the cases to a complete beam loss, due to the non-relativistic varying velocity of the particles. To avoid such a total beam loss, it is clear that some kind of retuning has to be performed to compensate the lack of acceleration due to the faulty cavity. We have to identify and develop fast failure recovery scenarios to ensure that such retuning can be performed in less than 1 second. 2 ways are investigated. The first way is to stop the beam to achieve the retuning (Scenario 1). The other way is to try to perform the retuning without stopping the beam (Scenario 2). The present analysis demonstrates on the beam dynamics point of view that a fast retuning procedure can be envisaged without stopping the beam (Scenario 2). Nevertheless, this Scenario 2 implies stringent specifications, especially on: - the fault detection time, that has to be extremely short (order of magnitude: 100 μs) and - the margins required on the accelerating field and RF power point of view, that are higher than in Scenario 1
Longitudinal beam dynamics for heavy ion fusion using WARPrz
WARPrz is a 2.5 dimensional, cylindrically symmetric, electrostatic, particle-in-cell code. It is part of the WARP family of codes which has been developed to study heavy ion fusion driver issues. WARPrz is being used to study the longitudinal dynamics of heavy ion beams including a longitudinal instability that is driven by the impedance of the LINAC accelerating modules. This instability is of concern because it can enhance longitudinal momentum spread; chromatic abhoration in the lens system restricts the amount of momentum spread allowed in the beam in the final focusing system. The impedance of the modules is modeled by a continuum of resistors and capacitors in parallel in WARPrz. We discuss simulations of this instability including the effect of finite temperature and reflection of perturbations off the beam ends. We also discuss intermittency of axial confining fields (''ears'' fields) as a seed for this instability
Transverse beam dynamics studies of a heavy ion induction linac
The multiple beam induction linac experiment (MBE-4) was built to study the accelerator physics of the low energy, electrostatically focused end of a driver for heavy ion inertial confinement fusion. In this machine four beams of Cs+ ions are accelerated through 24 common induction gaps while being focused in separate AG focusing channels. Each channel consists of a syncopated FODO lattice of 30 periods. The authors report results of the most recent studies of the transverse beam dynamics of a single drifting (180 keV) beam in this machine. The dependence of the emittance on the zero-current phase advance shows systematic variations which may be understood in the light of previous theoretical work on this topic. This result, unique to the beam parameters of a linac for heavy ion fusion, is discussed in the context of its implications for a driver design. In addition they discuss recent measurements of the motion of the beam centroid through the linac. These measurements, coupled with simulations, have proven to be a powerful tool in determining the presence of misalignment errors in the lattice of the accelerator
De Martino, Salvatore; De Siena, Silvio; Illuminati, Fabrizio
1998-01-01
A recent proposal (see quant-ph/9803068) to simulate semiclassical corrections to classical dynamics by suitable classical stochastic fluctuations is applied to the specific instance of charged beam dynamics in particle accelerators. The resulting picture is that the collective beam dynamics, at the leading semiclassical order in Planck constant can be described by a particular diffusion process, the Nelson process, which is time-reversal invariant. Its diffusion coefficient $\\sqrt{N}\\lambda_...
SIMULATION OF INTENSE BEAMS FOR HEAVY ION FUSION
Friedman, A
2004-06-10
Computer simulations of intense ion beams play a key role in the Heavy Ion Fusion research program. Along with analytic theory, they are used to develop future experiments, guide ongoing experiments, and aid in the analysis and interpretation of experimental results. They also afford access to regimes not yet accessible in the experimental program. The U.S. Heavy Ion Fusion Virtual National Laboratory and its collaborators have developed state-of-the art computational tools, related both to codes used for stationary plasmas and to codes used for traditional accelerator applications, but necessarily differing from each in important respects. These tools model beams in varying levels of detail and at widely varying computational cost. They include moment models (envelope equations and fluid descriptions), particle-in-cell methods (electrostatic and electromagnetic), nonlinear-perturbative descriptions (''{delta}f''), and continuum Vlasov methods. Increasingly, it is becoming clear that it is necessary to simulate not just the beams themselves, but also the environment in which they exist, be it an intentionally-created plasma or an unwanted cloud of electrons and gas. In this paper, examples of the application of simulation tools to intense ion beam physics are presented, including support of present-day experiments, fundamental beam physics studies, and the development of future experiments. Throughout, new computational models are described and their utility explained. These include Mesh Refinement (and its dynamic variant, Adaptive Mesh Refinement); improved electron cloud and gas models, and an electron advance scheme that allows use of larger time steps; and moving-mesh and adaptive-mesh Vlasov methods.
Issues and opportunities: beam simulations for heavy ion fusion
UCRL- JC- 134975 PREPRINT code offering 3- D, axisymmetric, and ''transverse slice'' (steady flow) geometries, with a hierarchy of models for the ''lattice'' of focusing, bending, and accelerating elements. Interactive and script- driven code steering is afforded through an interpreter interface. The code runs with good parallel scaling on the T3E. Detailed simulations of machine segments and of complete small experiments, as well as simplified full- system runs, have been carried out, partially benchmarking the code. A magnetoinductive model, with module impedance and multi- beam effects, is under study. experiments, including an injector scalable to multi- beam arrays, a high- current beam transport and acceleration experiment, and a scaled final- focusing experiment. These ''phase I'' projects are laying the groundwork for the next major step in HIF development, the Integrated Research Experiment (IRE). Simulations aimed directly at the IRE must enable us to: design a facility with maximum power on target at minimal cost; set requirements for hardware tolerances, beam steering, etc.; and evaluate proposed chamber propagation modes. Finally, simulations must enable us to study all issues which arise in the context of a fusion driver, and must facilitate the assessment of driver options. In all of this, maximum advantage must be taken of emerging terascale computer architectures, requiring an aggressive code development effort. An organizing principle should be pursuit of the goal of integrated and detailed source- to- target simulation. methods for analysis of the beam dynamics in the various machine concepts, using moment- based methods for purposes of design, waveform synthesis, steering algorithm synthesis, etc. Three classes of discrete- particle models should be coupled: (1) electrostatic/ magnetoinductive PIC simulations should track the beams from the source through the final- focusing optics, passing details of the time- dependent distribution function to
Issues and opportunities: beam simulations for heavy ion fusion
Friedman, A
1999-07-15
UCRL- JC- 134975 PREPRINT code offering 3- D, axisymmetric, and ''transverse slice'' (steady flow) geometries, with a hierarchy of models for the ''lattice'' of focusing, bending, and accelerating elements. Interactive and script- driven code steering is afforded through an interpreter interface. The code runs with good parallel scaling on the T3E. Detailed simulations of machine segments and of complete small experiments, as well as simplified full- system runs, have been carried out, partially benchmarking the code. A magnetoinductive model, with module impedance and multi- beam effects, is under study. experiments, including an injector scalable to multi- beam arrays, a high- current beam transport and acceleration experiment, and a scaled final- focusing experiment. These ''phase I'' projects are laying the groundwork for the next major step in HIF development, the Integrated Research Experiment (IRE). Simulations aimed directly at the IRE must enable us to: design a facility with maximum power on target at minimal cost; set requirements for hardware tolerances, beam steering, etc.; and evaluate proposed chamber propagation modes. Finally, simulations must enable us to study all issues which arise in the context of a fusion driver, and must facilitate the assessment of driver options. In all of this, maximum advantage must be taken of emerging terascale computer architectures, requiring an aggressive code development effort. An organizing principle should be pursuit of the goal of integrated and detailed source- to- target simulation. methods for analysis of the beam dynamics in the various machine concepts, using moment- based methods for purposes of design, waveform synthesis, steering algorithm synthesis, etc. Three classes of discrete- particle models should be coupled: (1) electrostatic/ magnetoinductive PIC simulations should track the beams from the source through the final- focusing optics
Automated analysis for detecting beams in laser wakefield simulations
Ushizima, Daniela M.; Rubel, Oliver; Prabhat, Mr.; Weber, Gunther H.; Bethel, E. Wes; Aragon, Cecilia R.; Geddes, Cameron G.R.; Cormier-Michel, Estelle; Hamann, Bernd; Messmer, Peter; Hagen, Hans
2008-07-03
Laser wakefield particle accelerators have shown the potential to generate electric fields thousands of times higher than those of conventional accelerators. The resulting extremely short particle acceleration distance could yield a potential new compact source of energetic electrons and radiation, with wide applications from medicine to physics. Physicists investigate laser-plasma internal dynamics by running particle-in-cell simulations; however, this generates a large dataset that requires time-consuming, manual inspection by experts in order to detect key features such as beam formation. This paper describes a framework to automate the data analysis and classification of simulation data. First, we propose a new method to identify locations with high density of particles in the space-time domain, based on maximum extremum point detection on the particle distribution. We analyze high density electron regions using a lifetime diagram by organizing and pruning the maximum extrema as nodes in a minimum spanning tree. Second, we partition the multivariate data using fuzzy clustering to detect time steps in a experiment that may contain a high quality electron beam. Finally, we combine results from fuzzy clustering and bunch lifetime analysis to estimate spatially confined beams. We demonstrate our algorithms successfully on four different simulation datasets.
Automated analysis for detecting beams in laser wakefield simulations
Laser wakefield particle accelerators have shown the potential to generate electric fields thousands of times higher than those of conventional accelerators. The resulting extremely short particle acceleration distance could yield a potential new compact source of energetic electrons and radiation, with wide applications from medicine to physics. Physicists investigate laser-plasma internal dynamics by running particle-in-cell simulations; however, this generates a large dataset that requires time-consuming, manual inspection by experts in order to detect key features such as beam formation. This paper describes a framework to automate the data analysis and classification of simulation data. First, we propose a new method to identify locations with high density of particles in the space-time domain, based on maximum extremum point detection on the particle distribution. We analyze high density electron regions using a lifetime diagram by organizing and pruning the maximum extrema as nodes in a minimum spanning tree. Second, we partition the multivariate data using fuzzy clustering to detect time steps in a experiment that may contain a high quality electron beam. Finally, we combine results from fuzzy clustering and bunch lifetime analysis to estimate spatially confined beams. We demonstrate our algorithms successfully on four different simulation datasets
Study of longitudinal dynamics in space-charge dominated beams
Tian, Kai
Modern accelerator applications, such as heavy ion fusion drivers, pulsed neutron sources, electron injectors for high-energy linear colliders, and X-ray Free Electron Lasers, demand beams with high intensity, low emittance and small energy spread. At low (non-relativistic) energies, the "electrostatic", collective interactions from space-charge forces existing in such intense beams play the dominant role; we characterize these beams as space-charge dominated beams. This dissertation presents numerous new findings on the longitudinal dynamics of a space-charge dominated beam, particularly on the propagation of density perturbations. In order to fully understand the complex physics of longitudinal space-charge waves, we combine the results of theory, computer simulation, and experiment. In the Long Solenoid Experimental system (LSE), with numerous diagnostic tools and techniques, we have, for the first time, experimentally measured the detailed energy profiles of longitudinal space-charge waves at different locations, both near the beam source and at the end of the transport system. Along with the current profiles, we have a complete set of experimental data for the propagation of space-charge waves. We compare these measured results to a 1-D theory and find better agreement for beams with perturbations in the linear regime, where the perturbation strength is less than 10%, than those with nonlinear perturbations. Using fast imaging techniques that we newly developed, we have, for the first time, obtained the progressive time-resolved images of longitudinal slices of a space-charge dominated beam. These images not only provide us time-resolved transverse density distribution of the beam, but also enable us to take time-resolved transverse phase space measurement using computerized tomography. By combining this information with the longitudinal energy measurement, we have, for the first time, experimentally constructed the full 6-D phase space. Part of the results
Beam dynamics issues in the FCC
AUTHOR|(CDS)2067437; Benedikt, Michael; Besana, Maria Ilaria; Bruce, Roderik; Bruning, Oliver; Buffat, Xavier; Burkart, Florian; Burkhardt, Helmut; Calatroni, Sergio; Cerutti, Francesco; Fartoukh, Stephane; Fiascaris, Maria; Garion, Cedric; Goddard, Brennan; Hofle, Wolfgang; Holzer, Bernhard; Jowett, John; Kersevan, Roberto; Martin, Roman; Mether, Lotta Maria; Milanese, Attilio; Pieloni, Tatiana; Redaelli, Stefano; Rumolo, Giovanni; Salvant, Benoit; Schaumann, Michaela; Schulte, Daniel; Chapochnikova, Elena; Stoel, Linda; Tambasco, Claudia; Tomas Garcia, Rogelio; Tommasini, Davide; Zimmermann, Frank; Guillermo Canton, Gerardo; Kornilov, Vladimir; Boine-Frankenheim, Oliver; Niedermayer, Uwe; Mitsuhashi, Toshiyuki; Ohmi, Kazuhito; Chance, Antoine; Dalena, Barbara; Payet, Jacques; Bambade, Philip; Faus-Golfe, Angeles; Molson, James; Biarrotte, Jean-Luc; Lachaize, Antoine; Fox, John D; Stupakov, Gennady; Abelleira, Jose; Cruz Alaniz, Emilia; Seryi, Andrei; Appleby, Robert Barrie; Boscolo, Manuela; Collamati, Francesco; Drago, Alessandro; Barranco Garcia, Javier; Khan, Shaukat; Riemann, Bernhard
2016-01-01
The international Future Circular Collider (FCC) study is designing hadron, lepton and lepton-hadron colliders based on a new 100 km tunnel in the Geneva region. The main focus and ultimate goal of the study are high luminosity proton-proton collisions at a centre-of-mass energy of 100 TeV, using 16 T Nb3Sn dipole magnets. Specific FCC beam dynamics issues are related to the large circumference, the high brightness—made available by radiation damping —, the small geometric emittance, unprecedented collision energy and luminosity, the huge amount of energy stored in the beam, large synchrotron radiation power, plus the injection scenarios. In addition to the FCC-hh proper, also a High-Energy LHC (HE-LHC) is being explored, using the FCC-hh magnet technology in the existing LHC tunnel, which can yield a centre-of-mass energy around 25 TeV.
Modeling and simulation of LHC beam-based collimator setup
Valentino, G; Assmann, R W; Burkart, F; Redaelli, S; Rossi, A; Lari, L
2012-01-01
In the 2011 Large Hadron Collider run, collimators were aligned for proton and heavy ion beams using a semiautomatic setup algorithm. The algorithm provided a reduction in the beam time required for setup, an elimination of beam dumps during setup and better reproducibility with respect to manual alignment. A collimator setup simulator was developed based on a Gaussian model of the beam distribution as well as a parametric model of the beam losses. A time-varying beam loss signal can be simulated for a given collimator movement into the beam. The simulation results and comparison to measurement data obtained during collimator setups and dedicated fills for beam halo scraping are presented. The simulator will then be used to develop a fully automatic collimator alignment algorithm.
Photorefractive dynamic holography using self-pumped phase conjugate beam
Arun Anand; C S Narayanamurthy
2006-03-01
Dynamic holography in photorefractive materials using self-pumped phase conjugate beam of the object beam itself as the other writing beam is proposed. Our detailed theoretical analysis shows four-fold increase in the diffraction efficiency of dynamic holograms if recorded using this geometry even in photorefractive crystal like BTO (having low optical activity) without applying external field. Detailed theoretical analysis is given.
Summary of session 3 on synchrotron radiation and beam dynamics
Shiltsev, V.; /Fermilab; Metral, E.; /CERN
2010-12-01
We summarize presentations, discussions and general conclusions of the Workshop session on 'Beam Dynamics Issues'. Major subjects include effects due to synchrotron radiation (SR), cryogenic loads, electron cloud, impedances, intra-beam scattering (IBS) and beam-beam interactions.
Simulations with dynamical HISQ quarks
Bazavov, A; DeTar, C; Freeman, W; Gottlieb, Steven; Heller, U M; Hetrick, J E; Laiho, J; Levkova, L; Oktay, M; Osborn, J; Sugar, R L; Toussaint, D; Van de Water, R S
2010-01-01
We report on the status of a program of generating and using configurations with four flavors of dynamical quarks, using the HISQ action. We study the lattice spacing dependence of physical quantities in these simulations, using runs at several lattice spacings, but with the light quark mass held fixed at two tenths of the strange quark mass. We find that the lattice artifacts in the HISQ simulations are much smaller than those in the asqtad simulations at the same lattice spacings and quark masses. We also discuss methods for setting the scale, or assigning a lattice spacing to ensembles run at unphysical parameters.
Beam stability & nonlinear dynamics. Formal report
Parsa, Z. [ed.
1996-12-31
his Report includes copies of transparencies and notes from the presentations made at the Symposium on Beam Stability and Nonlinear Dynamics, December 3-5, 1996 at the Institute for Theoretical Physics, University of California, Santa Barbara California, that was made available by the authors. Editing, reduction and changes to the authors contributions were made only to fulfill the printing and publication requirements. We would like to take this opportunity and thank the speakers for their informative presentations and for providing copies of their transparencies and notes for inclusion in this Report.
Beam stability ampersand nonlinear dynamics. Formal report
This report includes copies of transparencies and notes from the presentations made at the Symposium on Beam Stability and Nonlinear Dynamics, December 3-5, 1996 at the Institute for Theoretical Physics, University of California, Santa Barbara California, that was made available by the authors. Editing, reduction and changes to the authors contributions were made only to fulfill the printing and publication requirements. We would like to take this opportunity and thank the speakers for their informative presentations and for providing copies of their transparencies and notes for inclusion in this Report
Molecular beam studies of reaction dynamics
Lee, Y.T.
1987-03-01
Purpose of this research project is two-fold: (1) to elucidate detailed dynamics of simple elementary reactions which are theoretically important and to unravel the mechanism of complex chemical reactions or photo chemical processes which play an important role in many macroscopic processes and (2) to determine the energetics of polyatomic free radicals using microscopic experimental methods. Most of the information is derived from measurement of the product fragment translational energy and angular distributions using unique molecular beam apparati designed for these purposes.
Simulation of Electron-Beam Generating Plasma at Atmospheric Pressure
OUYANG Liang; LI Hong; LI Benben; ZHOU Junqing; YAN Hong; SU Tie; WANG Huihui; LIUWandong
2007-01-01
As electron-beam generating plasma is widely applied,the software tool EGS4(Electron-Gamma Shower) was used to simulate the transmission and energy deposition of electron-beam in air.The simulation results indicated that the range of the electron-beam was inversely proportional to the gas pressure in a wide range of gas pressure,and the electron-beam of 200 keV could generate a plasma with a density 1011 cm-3 in air of latm.In addition,the energy distribution of the beam-electron and plasma density profile produced by the beam were achieved.
The dynamics of radiation damage by focused ion beams in the ion beam synthesis
The following topics were covered: direct ion implantation, focused ion beams, cobalt silicides (CoSi2), RBS, ion beam synthesis, CoSi2 ion beam synthesis by focused ion beams in Si(111), germanium FIB implantation in Si(111), radiation damage at FIB implantation, models and simulation. (WL)
Parallel beam dynamics calculations on high performance computers
Faced with a backlog of nuclear waste and weapons plutonium, as well as an ever-increasing public concern about safety and environmental issues associated with conventional nuclear reactors, many countries are studying new, accelerator-driven technologies that hold the promise of providing safe and effective solutions to these problems. Proposed projects include accelerator transmutation of waste (ATW), accelerator-based conversion of plutonium (ABC), accelerator-driven energy production (ADEP), and accelerator production of tritium (APT). Also, next-generation spallation neutron sources based on similar technology will play a major role in materials science and biological science research. The design of accelerators for these projects will require a major advance in numerical modeling capability. For example, beam dynamics simulations with approximately 100 million particles will be needed to ensure that extremely stringent beam loss requirements (less than a nanoampere per meter) can be met. Compared with typical present-day modeling using 10,000 endash 100,000 particles, this represents an increase of 3 endash 4 orders of magnitude. High performance computing (HPC) platforms make it possible to perform such large scale simulations, which require 10 close-quote s of GBytes of memory. They also make it possible to perform smaller simulations in a matter of hours that would require months to run on a single processor workstation. This paper will describe how HPC platforms can be used to perform the numerically intensive beam dynamics simulations required for development of these new accelerator-driven technologies. copyright 1997 American Institute of Physics
On the dynamics of viscous masonry beams
Lucchesi, M.; Pintucchi, B.; Šilhavý, M.; Zani, N.
2015-05-01
In this paper, we consider the longitudinal and transversal vibrations of the masonry beams and arches. The basic motivation is the seismic vulnerability analysis of masonry structures that can be modeled as monodimensional elements. The Euler-Bernoulli hypothesis is employed for the system of forces in the beam. The axial force and the bending moment are assumed to consist of the elastic and viscous parts. The elastic part is described by the no-tension material, i.e., the material with no resistance to tension and which accounts for the cases of limitless, as well as bounded compressive strength. The adaptation of this material to beams has been developed in Orlandi (Analisi non lineare di strutture ad arco in muratura. Thesis, 1999) and Zani (Eur J Mech A/Solids 23:467-484, 2004). The viscous part amounts to the Kelvin-Voigt damping depending linearly on the time derivatives of the linearized strain and curvature. The dynamical equations are formulated, and a mathematical analysis of them is presented. Specifically, following Gajewski et al. (Nichtlineare Operatorgleichungen und Operatordifferentialgleichungen. Akademie-Verlag, Berlin, 1974), the theorems of existence, uniqueness and regularity of the solution of the dynamical equations are recapitulated and specialized for our purposes, to support the numerical analysis applied previously in Lucchesi and Pintucchi (Eur J Mech A/Solids 26:88-105, 2007 ). As usual, for that the Galerkin method has been used. As an illustration, two numerical examples (slender masonry tower and masonry arch) are presented in this paper with the applied forces corresponding to the acceleration in the earthquake in Emilia Romagna in May 29, 2012.
Investigation of the transverse beam dynamics in the thermal wave model with a functional method
Jang, Ji-Ho; Cho, Yong-sub; Kwon, Hyeok-jung
2006-01-01
We investigated the transverse beam dynamics in a thermal wave model by using a functional method. It can describe the beam optical elements separately with a kernel for a component. The method can be applied to general quadrupole magnets beyond a thin lens approximation as well as drift spaces. We found that the model can successfully describe the PARMILA simulation result through an FODO lattice structure for the Gaussian input beam without space charge effects.
Simulation of ion beam injection and extraction in an EBIS
An example simulation of Au+ charge breeding using FAR-TECH’s integrated EBIS (electron beam ion source) modeling toolset is presented with the emphasis on ion beam injection and extraction. The trajectories of injected ions are calculated with PBGUNS (particle beam gun simulation) self-consistently by including the space charges from both ions and electrons. The ion beam, starting with initial conditions within the 100% acceptance of the electron beam, is then tracked by EBIS-PIC (particle-in-cell EBIS simulation code). In the trap, the evolution of the ion charge state distribution is estimated by charge state estimator. The extraction of charge bred ions is simulated with PBGUNS. The simulations of the ion injections show significant ion space charge effects on beam capture efficiency and the ionization efficiency
Simulation of ion beam injection and extraction in an EBIS
Zhao, L., E-mail: zhao@far-tech.com; Kim, J. S. [FAR-TECH, Inc., San Diego, California 92121 (United States)
2016-02-15
An example simulation of Au+ charge breeding using FAR-TECH’s integrated EBIS (electron beam ion source) modeling toolset is presented with the emphasis on ion beam injection and extraction. The trajectories of injected ions are calculated with PBGUNS (particle beam gun simulation) self-consistently by including the space charges from both ions and electrons. The ion beam, starting with initial conditions within the 100% acceptance of the electron beam, is then tracked by EBIS-PIC (particle-in-cell EBIS simulation code). In the trap, the evolution of the ion charge state distribution is estimated by charge state estimator. The extraction of charge bred ions is simulated with PBGUNS. The simulations of the ion injections show significant ion space charge effects on beam capture efficiency and the ionization efficiency.
Simulation of ion beam injection and extraction in an EBIS
Zhao, L.; Kim, J. S.
2016-02-01
An example simulation of Au+ charge breeding using FAR-TECH's integrated EBIS (electron beam ion source) modeling toolset is presented with the emphasis on ion beam injection and extraction. The trajectories of injected ions are calculated with PBGUNS (particle beam gun simulation) self-consistently by including the space charges from both ions and electrons. The ion beam, starting with initial conditions within the 100% acceptance of the electron beam, is then tracked by EBIS-PIC (particle-in-cell EBIS simulation code). In the trap, the evolution of the ion charge state distribution is estimated by charge state estimator. The extraction of charge bred ions is simulated with PBGUNS. The simulations of the ion injections show significant ion space charge effects on beam capture efficiency and the ionization efficiency.
Simulation of intense beam bunching using 3D PIC method
Most of the ion sources produce continuous beam of charged particles. In a cyclotron using such an external ion source, only a small fraction of the injected continuous beam is accepted in the central region for further acceleration. By transforming the continuous beam into a suitably bunched beam using a buncher prior to injection, the amount of accepted particles in the central region of cyclotron can be increased. To compress the continuous beam longitudinally one needs to impose a velocity modulation at the buncher gap which results in density modulation as the beam advances. In the case of low beam current the velocity modulation of the beam has very little effect on the transverse envelope of the beam. However, in the case of high intensity beams, the space-charge force introduces much collective behaviour and increase of current in the specified bunch width affects the transverse dynamics
Study of beam dynamics at cooler synchrotron TARN-II
Several kinds of beam diagnostic instruments, have been developed at cooler-synchrotron TARN-II. These are intended to study beam dynamics at low beam current of several microamperes and then have high sensitivity of good S/N ratio. In addition, the acceleration system, especially low level RF system, has been improved to attain the maximum beam energy. With the successful performance of these instrumentations, the study of beam dynamics are presently being carried out. For example, the synchrotron acceleration of the light ions was achieved up to 220 MeV/u without any beam loss. (author)
Simulation of wire-compensation of long range beam beam interaction in high energy accelerators
Full text: We present weak-strong simulation results for the effect of long-range beam-beam (LRBB) interaction in LHC as well as for proposed wire compensation schemes or wire experiments, respectively. In particular, we discuss details of the simulation model, instability indicators, the effectiveness of compensation, the difference between nominal and PACMAN bunches for the LHC, beam experiments, and wire tolerances. The simulations are performed with the new code BBTrack. (author)
Simulation and Experimental Studies on Composite Beams
M. Abhinay
2014-09-01
Full Text Available A composite beam a one dimensional structure or a rod all of them are sectional dimensions in which width and height are much smaller in comparison to the structure. In structural applications longer beams are more frequently used. In this work a composite beam is manufactured with glass and epoxy combination. And stress analysis is carried out using derived analytical expressions. This research work carried out will enable to determine the beam strength due to bending loads. The importance of fiber reinforcement in the manufacturing of the beam is studied in terms of bending strength of the beam. Mat lab codes are generated to implement analytical equations of the composite beam. The analytical results are validated by performing experiments on composite beams. In this investigation, two different composite beams have been tested and compared the experimental results with the analytical results.
Numerical simulation to study the transient self focusing of laser beam in plasma
In this paper, we present the numerical simulation for the coupled system of equations governing the dynamics of laser and Ion Acoustic Wave (IAW) in a collisionless plasma, when the coupling between the waves is through ponderomotive non-linearity. The nonlinear evolution of the laser beam is studied when the pump laser is perturbed by a periodic perturbation. By changing the perturbation wave number, we have studied its effect on the nonlinear evolution pattern of laser beam. In order to have a physical insight into the nonlinear dynamics of laser beam evolution in time and space, we have studied the laser and IAW spectra containing spatial harmonics. The magnitude of these harmonics changes with time and leads to time dependent localization of laser beam in spatial domain. The nonlinear dynamics of this localization is investigated in detail by using simulation and a semi-analytical model
Vehicle dynamics modeling and simulation
Schramm, Dieter; Bardini, Roberto
2014-01-01
The authors examine in detail the fundamentals and mathematical descriptions of the dynamics of automobiles. In this context different levels of complexity will be presented, starting with basic single-track models up to complex three-dimensional multi-body models. A particular focus is on the process of establishing mathematical models on the basis of real cars and the validation of simulation results. The methods presented are explained in detail by means of selected application scenarios.
Implementation of geometrically exact beam element for nonlinear dynamics modeling
Wang, Jielong, E-mail: jelon.wang@gmail.com [Commercial Aircraft Corporation of China, Ltd., Beijing Aeronautical Science & Technology Research Institute (China)
2015-12-15
The paper develops a new type of geometrically exact beam element featuring large displacements and rotations together with small warping. The dimension reduction approach based on variational asymptotic method has been explored, and a linear two-dimensional finite element procedure has been implemented to predict the cross-sectional stiffness and recover the cross-sectional strain fields of the beam. The total and incremental variables mixed formula of governing equations of motion is presented, in which the Wiener–Milenković parameters are selected to vectorize the finite rotation. The dynamic problem of geometrically exact beam has been solved by the implicit Radau IIA algorithms, the time histories of large translations and rotations with small three-dimensional warping have been integrated. Numerical simulations have been performed and the results have been compared to those of commercial software LS-DYNA. It can be concluded that the current modeling approach features high accuracy and that the new geometrically exact beam with warping is robust enough to predict large deformations with small strain.
Implementation of geometrically exact beam element for nonlinear dynamics modeling
The paper develops a new type of geometrically exact beam element featuring large displacements and rotations together with small warping. The dimension reduction approach based on variational asymptotic method has been explored, and a linear two-dimensional finite element procedure has been implemented to predict the cross-sectional stiffness and recover the cross-sectional strain fields of the beam. The total and incremental variables mixed formula of governing equations of motion is presented, in which the Wiener–Milenković parameters are selected to vectorize the finite rotation. The dynamic problem of geometrically exact beam has been solved by the implicit Radau IIA algorithms, the time histories of large translations and rotations with small three-dimensional warping have been integrated. Numerical simulations have been performed and the results have been compared to those of commercial software LS-DYNA. It can be concluded that the current modeling approach features high accuracy and that the new geometrically exact beam with warping is robust enough to predict large deformations with small strain
Simulation and Experimental Studies on Composite Beams
M. Abhinay; P.Sampath Rao
2014-01-01
A composite beam a one dimensional structure or a rod all of them are sectional dimensions in which width and height are much smaller in comparison to the structure. In structural applications longer beams are more frequently used. In this work a composite beam is manufactured with glass and epoxy combination. And stress analysis is carried out using derived analytical expressions. This research work carried out will enable to determine the beam strength due to bending loads. ...
De Martino, S; Illuminati, F; Martino, Salvatore De; Siena, Silvio De; Illuminati, Fabrizio
1999-01-01
A recent proposal (see quant-ph/9803068) to simulate semiclassical corrections to classical dynamics by suitable classical stochastic fluctuations is applied to the specific instance of charged beam dynamics in particle accelerators. The resulting picture is that the collective beam dynamics, at the leading semiclassical order in Planck constant can be described by a particular diffusion process, the Nelson process, which is time-reversal invariant. Its diffusion coefficient $\\sqrt{N}\\lambda_{c}$ represents a semiclassical unit of emittance (here $N$ is the number of particles in the beam, and $\\lambda_{c}$ is the Compton wavelength). The stochastic dynamics of the Nelson type can be easily recast in the form of a Schroedinger equation, with the semiclassical unit of emittance replacing Planck constant. Therefore we provide a physical foundation to the several quantum-like models of beam dynamics proposed in recent years. We also briefly touch upon applications of the Nelson and Schroedinger formalisms to inc...
Molecular dynamics simulation of benzene
Trumpakaj, Zygmunt; Linde, Bogumił B. J.
2016-03-01
Intermolecular potentials and a few models of intermolecular interaction in liquid benzene are tested by Molecular Dynamics (MD) simulations. The repulsive part of the Lennard-Jones 12-6 (LJ 12-6) potential is too hard, which yields incorrect results. The exp-6 potential with a too hard repulsive term is also often used. Therefore, we took an expa-6 potential with a small Gaussian correction plus electrostatic interactions. This allows to modify the curvature of the potential. The MD simulations are carried out in the temperature range 280-352 K under normal pressure and at experimental density. The Rayleigh scattering of depolarized light is used for comparison. The results of MD simulations are comparable with the experimental values.
Dynamic simulations of tissue welding
Maitland, D.J.; Eder, D.C.; London, R.A.; Glinsky, M.E. [and others
1996-02-01
The exposure of human skin to near-infrared radiation is numerically simulated using coupled laser, thermal transport and mass transport numerical models. The computer model LATIS is applied in both one-dimensional and two-dimensional geometries. Zones within the skin model are comprised of a topical solder, epidermis, dermis, and fatty tissue. Each skin zone is assigned initial optical, thermal and water density properties consistent with values listed in the literature. The optical properties of each zone (i.e. scattering, absorption and anisotropy coefficients) are modeled as a kinetic function of the temperature. Finally, the water content in each zone is computed from water diffusion where water losses are accounted for by evaporative losses at the air-solder interface. The simulation results show that the inclusion of water transport and evaporative losses in the model are necessary to match experimental observations. Dynamic temperature and damage distributions are presented for the skin simulations.
SPP Beamline Design and Beam Dynamics
Turemen, G; Yasatekin, B; Yildiz, V; Celik, M; Alacakir, A; Unel, G; Mete, O
2014-01-01
The Radio Frequency Quadrupole of SANAEM Project Prometheus will be a demonstration and educational machine which will accelerate protons from 20 keV to 1.5 MeV. The project is funded by Turkish Atomic Energy Authority and it will be located at Saraykoy Nuclear Research and Training Center in Ankara. The SPP beamline consists of a multi-cusp H+ ion source, a Low Energy Beam Transport line and a four-vane RFQ operating at 352.2 MHz. The design studies for the multi-cusp ion source (RF or DC) were performed with IBSimu and SIMION software packages. The source has already been produced and currently undergoes extensive testing. There is also a preliminary design for the solenoid based LEBT, POISSON and PATH were used in parallel for the preliminary design. Two solenoid magnets are produced following this design. The RFQ design was made using LIDOS.RFQ.Designer and it was crosschecked with a home-grown software package, DEMIRCI. The initial beam dynamics studies have been performed with both LIDOS and TOUTATIS. T...
Synchrotron radiation damping, intrabeam scattering and beam-beam simulations for HE-LHC
Valishev, A
2011-01-01
The proposed High-Energy LHC project presents an unusual combination of strong synchrotron radiation (SR) damping and intrabeam scattering (IBS), which is not seen in present-day hadron colliders. The subject of investigation reported in this paper was the simulation of beam-beam effect for the HE-LHC parameters. Parameters of SR and IBS are calculated, and the luminosity evolution is simulated in the absence of beambeam interaction. Then, a weak-strong numerical simulation is used to predict the effect of beam-beam interaction on particle losses and emittance evolution.
Dynamic bowtie filter for cone-beam/multi-slice CT.
Fenglin Liu
Full Text Available A pre-patient attenuator ("bowtie filter" or "bowtie" is used to modulate an incoming x-ray beam as a function of the angle of the x-ray with respect to a patient to balance the photon flux on a detector array. While the current dynamic bowtie design is focused on fan-beam geometry, in this study we propose a methodology for dynamic bowtie design in multi-slice/cone-beam geometry. The proposed 3D dynamic bowtie is an extension of the 2D prior art. The 3D bowtie consists of a highly attenuating bowtie (HB filled in with heavy liquid and a weakly attenuating bowtie (WB immersed in the liquid of the HB. The HB targets a balanced flux distribution on a detector array when no object is in the field of view (FOV. The WB compensates for an object in the FOV, and hence is a scaled-down version of the object. The WB is rotated and translated in synchrony with the source rotation and patient translation so that the overall flux balance is maintained on the detector array. First, the mathematical models of different scanning modes are established for an elliptical water phantom. Then, a numerical simulation study is performed to compare the performance of the scanning modes in the cases of the water phantom and a patient cross-section without any bowtie and with a dynamic bowtie. The dynamic bowtie can equalize the numbers of detected photons in the case of the water phantom. In practical cases, the dynamic bowtie can effectively reduce the dynamic range of detected signals inside the FOV. Furthermore, the WB can be individualized using a 3D printing technique as the gold standard. We have extended the dynamic bowtie concept from 2D to 3D by using highly attenuating liquid and moving a scale-reduced negative copy of an object being scanned. Our methodology can be applied to reduce radiation dose and facilitate photon-counting detection.
Molecular beam studies of reaction dynamics
Lee, Y.T. [Lawrence Berkeley Laboratory, CA (United States)
1993-12-01
The major thrust of this research project is to elucidate detailed dynamics of simple elementary reactions that are theoretically important and to unravel the mechanism of complex chemical reactions or photochemical processes that play important roles in many macroscopic processes. Molecular beams of reactants are used to study individual reactive encounters between molecules or to monitor photodissociation events in a collision-free environment. Most of the information is derived from measurement of the product fragment energy, angular, and state distributions. Recent activities are centered on the mechanisms of elementary chemical reactions involving oxygen atoms with unsaturated hydrocarbons, the dynamics of endothermic substitution reactions, the dependence of the chemical reactivity of electronically excited atoms on the alignment of excited orbitals, the primary photochemical processes of polyatomic molecules, intramolecular energy transfer of chemically activated and locally excited molecules, the energetics of free radicals that are important to combustion processes, the infrared-absorption spectra of carbonium ions and hydrated hydronium ions, and bond-selective photodissociation through electric excitation.
Beam-beam simulations with non-Gaussian distributions for SLC and SLC-2000
Due to various upstream beam manipulations, the longitudinal bunch shape at the interactions point of the Stanford Linear Collider (SLC) is highly non-Gaussian. In this paper, we report beam-beam simulations with realistic longitudinal bunch shapes for the present SLC parameters and for the SLC-2000 luminosity upgrade. The simulation results allow us to estimate the luminosity enhancement due to the pinch effect and to find optimum parameter settings for the bunch compressor and the linac
Simulation of Head-on Beam-Beam Limitations in Future High Energy Colliders
Buffat, Xavier; Florio, Adrien; Pieloni, Tatiana; Tambasco, Claudia
2016-01-01
The Future Circular Hadron Collider (FCC-hh) project calls for studies in a new regime of beam-beam interactions. While the emittance damping due to synchrotron radiation is still slower than in past or existing lepton colliders, it is significantly larger than in other hadron colliders. The slow reduction of the emittance is profitable for higher luminosity in term of transverse beam size at the interaction points and also to mitigate long-range beam-beam effects, potentially allowing for a reduction of the crossing angle between the beams during the operation. In such conditions, the strength of head-on beam-beam interactions increases, potentially limiting the beam brightness. 4D weak-strong and strong-strong simulations are performed in order to assess these limitations.
Longitudinal Beam Dynamics Studies at CTF3 And Pulse Compressor Controlling
Shaker, S H
2009-01-01
The aim of the CLIC Test Facility CTF3, built at CERN by an international collaboration, is to address the main feasibility issues of the CLIC electron-positron linear collider technology by 2010. One key-issue studied at CTF3 is the generation of the very high current drive beam, used in CLIC as the RF power source. It is particularly important to simulate and control the drive beam longitudinal dynamics in the drive beam generation complex, since it directly affects the efficiency and stability of the RF power production process. In this thesis how to use pulse compressors to achieve a high RF power in the drive beam accelerator is discussed. We also describe the ongoing effort in modelling the longitudinal evolution of the CTF3 drive beam and compare the simulations with experimental results. Our study is based on the single bunch simulation.
Direct simulation of ion beam induced stressing and amorphization of silicon
Beardmore, Keith M.; Gronbech-Jensen, Niels
1999-01-01
Using molecular dynamics (MD) simulation, we investigate the mechanical response of silicon to high dose ion-irradiation. We employ a realistic and efficient model to directly simulate ion beam induced amorphization. Structural properties of the amorphized sample are compared with experimental data and results of other simulation studies. We find the behavior of the irradiated material is related to the rate at which it can relax. Depending upon the ability to deform, we observe either the ge...
Beam simulation tools for GEANT4 (and neutrino source applications)
Geant4 is a tool kit developed by a collaboration of physicists and computer professionals in the High Energy Physics field for simulation of the passage of particles through matter. The motivation for the development of the Beam Tools is to extend the Geant4 applications to accelerator physics. Although there are many computer programs for beam physics simulations, Geant4 is ideal to model a beam going through material or a system with a beam line integrated to a complex detector. There are many examples in the current international High Energy Physics programs, such as studies related to a future Neutrino Factory, a Linear Collider, and a very Large Hadron Collider
Ewen, Hamish Maclean; Riccardo, Bartolini
The non-linear beam dynamics of a circular accelerator, such as the Large Hadron Collider, can have a significant impact on its operation. In order to avoid limitations on the performance reach of the accelerator, and ensure machine protection, it is vital that the beam dynamics are well understood and controlled. This thesis presents the results of studies of non-linear beam dynamics undertaken on the Large Hadron Collider at CERN, during the 2010 to 2013 period. It sets out to quantify the understanding of the non-linear beam dynamics through the comparison of beam-based measurements to simulation, and where able and appropriate seeks to explain deviations of measurement from the model, and define corrections for relevant aspects of the dynamics. The analyses presented in this thesis represent considerable advances in the understanding of the LHC beam dynamics which should allow for an improved operation of the machine in the coming years.
Simulations for CLIC Drive Beam Linac
Aksoy, Avni
2012-01-01
The Drive Beam Linac of the Compact Linear Collider (CLIC) has to accelerate an electron beam with 4.2 A up to 2.4 GeV in almost fully-loaded structures. The pulse contains about 70000 bunches, one in every second rf bucket, and has a length of 140 $\\mu$s. The beam stability along the beamline is of concern for such a high current and pulse length. We present different options for the lattice of the linac based on FODO, triplet and doublet cells and compare the transverse instability for each lattice including the effects of beam jitter, alignment and beam-based correction. Additionally longitudinal stability is discussed for different bunch compressors using FODO type of lattice.
Computer simulation of scalar vortex beams LG0L in time-varying random inhomogeneous media
Sennikov, Victor A.; Konyaev, Petr A.; Lukin, Vladimir P.
2015-11-01
Computer modeling of vortex beams LG0L propagation in random inhomogeneous media using the dynamic time varying algorithm of media evolution is presented. The temporal power spectrum of vortex beam LG0L has been investigated by numerical simulation of propagation through atmospheric turbulence. The split-step Fourier method was used for solving scalar wave equation for a randomly inhomogeneous medium with a power-law Kolmogorov spectrum.
Spin dynamics of electron beams in circular accelerators
Experiments using high energy beams of spin polarized, charged particles still prove to be very helpful in disclosing a deeper understanding of the fundamental structure of matter. An important aspect is to control the beam properties, such as brilliance, intensity, energy, and degree of spin polarization. In this context, the present studies show various numerical calculations of the spin dynamics of high energy electron beams in circular accelerators. Special attention has to be paid to the emission of synchrotron radiation, that occurs when deflecting charged particles on circular orbits. In the presence of the fluctuation of the kinetic energy due to the photon emission, each electron spin moves non-deterministically. This stochastic effect commonly slows down the speed of all numeric estimations. However, the shown simulations cover - using appropriate approximations - trackings for the motion of thousands of electron spins for up to thousands of turns. Those calculations are validated and complemented by empirical investigations at the electron stretcher facility ELSA of the University of Bonn. They can largely be extended to other boundary conditions and thus, can be consulted for new accelerator layouts.
Simulations of the LEDA LEBT H+ beam
The computer codes TRACE and SCHAR model the Low-Energy Demonstration Accelerator (LEDA) Low-Energy Beam Transport (LEBT) for 75-keV, 110-mA, dc H+ beams. Solenoid-lens location studies verify that the proposed LEBT design gives a near-optimum match to the LEDA RFQ. The desired RFQ transmission (≥ 90%) and output emittance (≤ 0.22 π mm mrad, transverse) are obtained when PARMTEQM transports the file for the SCHAR-generated optimum beam through the RFQ
Large Dynamic Range Beam Profile Measurements with Low Current Electron Beams
Large dynamic range [Peak/Noise > 105] beam profile measurements are routinely performed in the Hall-B beamline at Jefferson Lab. These measurements are made with a 1 to 10nA electron beam current with energies between 1 to 6 GeV. The electron beam scatters off of a thin [25 mu-m] W or Fe wire and the scattered particle/shower is detected via scintillation or Cerenkov light several meters downstream of the wire. This light is converted to an electrical pulse via photomultiplier tubes [PMT]. The PMT readout and wire motion are controlled and synchronized by VME electronics. This report describes results on increasing the dynamic range by using multiple wires of varying diameters. Profile measurements with this large dynamic range can be of use for machines with very large beam currents (ERL) where any FR-actional beam loss represents a significant amount of beam power [1,2
Chaotic dynamics of flexible beams with piezoelectric and temperature phenomena
Krysko, V.A. [Department of Mathematics and Modeling, Saratov State Technical University, Politehnicheskaya 77, 410054 Saratov (Russian Federation); Awrejcewicz, J., E-mail: awrejcew@p.lodz.pl [Lodz University of Technology, Department of Automation and Biomechanics, 1/15 Stefanowski St., 90-924 Lodz (Poland); Warsaw University of Technology, Department of Vehicles, 84 Narbutta St., 02-524 Warsaw (Poland); Kutepov, I.E.; Zagniboroda, N.A.; Papkova, I.V.; Serebryakov, A.V. [Department of Mathematics and Modeling, Saratov State Technical University, Politehnicheskaya 77, 410054 Saratov (Russian Federation); Krysko, A.V. [Department of Applied Mathematics and Systems Analysis, Saratov State Technical University, Politehnicheskaya 77, 410054 Saratov (Russian Federation)
2013-11-01
The Euler–Bernoulli kinematic model as well as the von Kármán geometric non-linearity are used to derive the PDEs governing flexible beam vibrations. The beam is embedded into a 2D temperature field, and its surface is subjected to action of the electric potential. We report how an increase of the exciting load amplitude yields the beam turbulent behavior, and how the temperature changes a scenario from a regular/laminar to spatio-temporal/turbulent dynamics. Both classical Fourier analysis and Morlet wavelets are used to monitor a strong influence of temperature on regular and chaotic beam dynamics.
Humanoid Robot Simulator: A Realistic Dynamics Approach
Lima, José; Gonçalves, José; Costa, Paulo; Moreira, António
2008-01-01
This paper describes a humanoid robot simulator with realistic dynamics. As simulation is a powerful tool for speeding up the control software development, the suggested accurate simulator allows to accomplish this goal. The simulator, based on the Open Dynamics Engine and GLScene graphics library, provides instant visual feedback and allows the user to test any control strategy without damaging the real robot in the early stages of the development. The proposed simulator also captures some c...
Design and simulation of a superconducting gun cavity with 3½ cells have been studied in order to give the first push to the electron beam for the linear accelerating system at The Institute of Accelerator Technologies at Ankara University. Electrons are accelerated through the gun cavity with the help of the Radiofrequency power suppliers from cryogenic systems. Accelerating gradient should be as high as possible to accelerate electron beam inside the cavity. In this study, electron beam reaches to 9.17 MeV energy at the end of the gun cavity with the accelerating gradient; Ec=19.21 MV/m. 1.3 GHz gun cavity consists of three TESLA-like shaped cells while the special designed gun-cell includes a cathode plug. Optimized important beam parameters inside the gun cavity, average beam current 3 mA, transverse emittance 2.5 mm mrad, repetition rate 30 MHz and other parameters are obtained for the SASE-FEL System. The Superfish/Poisson program is used to design each cell of the superconducting cavity. Superconducting gun cavity and Radiofrequency properties are studied by utilizing 2D Superfish/Poisson, 3D Computer Simulation Technology Microwave Studio, and 3D Computer Simulation Technology Particle Studio. Superfish/Poisson is also used to optimize the geometry of the cavity cells to get the highest accelerating gradient. The behavior of the particles along the beamline is included in this study. ASTRA Code is used to track the particles
Yildiz, H. Duran; Cakir, R.; Porsuk, D.
2015-06-01
Design and simulation of a superconducting gun cavity with 3½ cells have been studied in order to give the first push to the electron beam for the linear accelerating system at The Institute of Accelerator Technologies at Ankara University. Electrons are accelerated through the gun cavity with the help of the Radiofrequency power suppliers from cryogenic systems. Accelerating gradient should be as high as possible to accelerate electron beam inside the cavity. In this study, electron beam reaches to 9.17 MeV energy at the end of the gun cavity with the accelerating gradient; Ec=19.21 MV/m. 1.3 GHz gun cavity consists of three TESLA-like shaped cells while the special designed gun-cell includes a cathode plug. Optimized important beam parameters inside the gun cavity, average beam current 3 mA, transverse emittance 2.5 mm mrad, repetition rate 30 MHz and other parameters are obtained for the SASE-FEL System. The Superfish/Poisson program is used to design each cell of the superconducting cavity. Superconducting gun cavity and Radiofrequency properties are studied by utilizing 2D Superfish/Poisson, 3D Computer Simulation Technology Microwave Studio, and 3D Computer Simulation Technology Particle Studio. Superfish/Poisson is also used to optimize the geometry of the cavity cells to get the highest accelerating gradient. The behavior of the particles along the beamline is included in this study. ASTRA Code is used to track the particles.
Yildiz, H. Duran, E-mail: hdyildiz@ankara.edu.tr [Institute of Accelerator Technologies, Ankara University, Ankara (Turkey); Cakir, R. [Nanotechnology Engineering Department, Recep Tayyip Erdogan University, Rize (Turkey); Porsuk, D. [Physics Department, Dumlupinar University, Kutahya (Turkey)
2015-06-11
Design and simulation of a superconducting gun cavity with 3½ cells have been studied in order to give the first push to the electron beam for the linear accelerating system at The Institute of Accelerator Technologies at Ankara University. Electrons are accelerated through the gun cavity with the help of the Radiofrequency power suppliers from cryogenic systems. Accelerating gradient should be as high as possible to accelerate electron beam inside the cavity. In this study, electron beam reaches to 9.17 MeV energy at the end of the gun cavity with the accelerating gradient; E{sub c}=19.21 MV/m. 1.3 GHz gun cavity consists of three TESLA-like shaped cells while the special designed gun-cell includes a cathode plug. Optimized important beam parameters inside the gun cavity, average beam current 3 mA, transverse emittance 2.5 mm mrad, repetition rate 30 MHz and other parameters are obtained for the SASE-FEL System. The Superfish/Poisson program is used to design each cell of the superconducting cavity. Superconducting gun cavity and Radiofrequency properties are studied by utilizing 2D Superfish/Poisson, 3D Computer Simulation Technology Microwave Studio, and 3D Computer Simulation Technology Particle Studio. Superfish/Poisson is also used to optimize the geometry of the cavity cells to get the highest accelerating gradient. The behavior of the particles along the beamline is included in this study. ASTRA Code is used to track the particles.
Simulations of proton beam depth-dose distributions
Proton beams are successfully used in radiotherapy. A correct modification of beam parameters allows one to spare normal surrounding tissues from radiation action. Our work is focused on passive beam-shaping techniques, which are used to modify the proton beam properties. The beam passes through the scattering system, which consists of scattering materials, energy degraders, drift spaces and collimators. In order to model the proton beam transport through the scattering system, the new Monte Carlo (MC) computer code Track has been developed. The code Track can predict output proton beam parameters modulated by various system adjustments and helps to optimize them. It calculates a beam profile, creates beam emittance diagram at a specified position of the system and predicts proton beam depth-dose distribution in a water phantom. In addition it calculates beam losses on individual components. We present a physical model of the beam transport calculations and algorithm implemented in a code Track. We compared the Track code calculations of depth-dose distributions in water phantom with experimental data and with a set of MC calculations in the FLUKA code. The accuracy of simulation results and calculation time in Track code are observed
Egorov, E. N.; Koronovskii, A. A.; Kurkin, S. A.; Hramov, A. E.
2013-11-01
Results of numerical simulations and analysis of the formation and nonlinear dynamics of the squeezed state of a helical electron beam in a vircator with a magnetron injection gun as an electron source and with additional electron deceleration are presented. The ranges of control parameters where the squeezed state can form in such a system are revealed, and specific features of the system dynamics are analyzed. It is shown that the formation of a squeezed state of a nonrelativistic helical electron beam in a system with electron deceleration is accompanied by low-frequency longitudinal dynamics of the space charge.
Egorov, E. N., E-mail: evgeniy.n.egorov@gmail.com; Koronovskii, A. A.; Kurkin, S. A.; Hramov, A. E. [Chernyshevsky Saratov State University (Russian Federation)
2013-11-15
Results of numerical simulations and analysis of the formation and nonlinear dynamics of the squeezed state of a helical electron beam in a vircator with a magnetron injection gun as an electron source and with additional electron deceleration are presented. The ranges of control parameters where the squeezed state can form in such a system are revealed, and specific features of the system dynamics are analyzed. It is shown that the formation of a squeezed state of a nonrelativistic helical electron beam in a system with electron deceleration is accompanied by low-frequency longitudinal dynamics of the space charge.
Longitudinal motion in the SNS external beam simulation experiment
Using the envelope equation of Neuffer for longitudinal bunch motion, basic parameters are obtained for the SNS/HIF external beam simulation experiment. Solutions of the equation are calculated for zero space charge and zero external focussing and parameters given for bunch compression in the SNS synchrotron and in the external beam line. (U.K.)
Simulations of LEIR Injection Line Beam Position Monitors
Maltseva, Mariya
2016-01-01
In this paper sensitivity characteristics of a beam position monitor are described. Characteristics are obtained during the simulations in CST Studio, the results are compared with the calculated values. The results for a low-beta beam and with a wire are compared.
Beam Dynamics Based Design of Solenoid Channel for TAC Proton Linac
Kisoglu, H F
2014-01-01
Today a linear particle accelerator (linac), in which electric and magnetic fields are of vital importance, is one of the popular energy generation sources like Accelerator Driven System (ADS). A multipurpose, including primarily ADS, proton linac with energy of ~2 GeV is planned to constitute within the Turkish Accelerator Center (TAC) project collaborated by more than 10 Turkish universities. A Low Energy Beam Transport (LEBT) channel with two solenoids is a subcomponent of this linac. It transports the proton beam ejected by an ion source, and matches it with the Radio Frequency Quadrupole (RFQ) that is an important part of the linac. The LEBT channel would be consisted of two focusing solenoids and some diagnostic elements such as faraday cup, BC transformers, etc. This paper includes a beam dynamical design and optimization study of LEBT channel for TAC proton linac done by using a beam dynamics simulation code PATH MANAGER and comparing of the simulation results with the theoretical expectations.
Simulation of Fault Zone Dynamics
Mora, P.; Abe, S.; Place, D.
2002-12-01
Particle models such as the discrete element model for granular assemblies and the lattice solid model provide a means to study the dynamics of fault zones. The lattice solid model was developed with the aim of progressively building up the capacity to simulate all relevent physical processes in fault zones. The present implementation of the model is able to simulate the dynamics of a granular lattice consisting of bonded or unbonded circular (2D) or spherical (3D) particles. Thermal effects (frictional hear generation, thermal expansion, heat flow) and pore fluid effects (heat induced pore pressure gradients and the consequent Darcian flow and impact on effective friction) can be modelled. Past work involving both circular particles and non-circular grains constructed as groups of bonded particles have demonstrated that grain shape has a fundamental impact on zero-th order behaviour. When circular particles are used, rolling is the most efficient means to accomodate slip of a simulated fault gouge layer leading to unrealistically low friction, typically around 0.2. This is consistent with laboratory results by Mair and Marone which have demonstrated that gouge consisting entirely of spherical beads shows a lower coefficient of friction than gouge containing irregular shaped particles. Recent work comparing quasi-2D laboratory results using pasta (Marone) with 2D numerical results (Morgan) have confirmed that numerical and laboratory results with circular ``particles'' are in agreement. When irregular grains are modelled at the lowest scale, the friction of simulated gouge layers matches with laboratory observations of rock friction (μ ~ 0.6) and is insentitive to the value used for interparticle friction (Mora et al, 2000). This indicates a self-regulation mechanism is occurring in which the group behaviour of the gouge layer remains constant at around 0.6 by balancing the amount of slip and rolling of grains within the gouge layer. A limitation of these studies
APPLICATION OF UAL TO HIGH INTENSITY BEAM DYNAMICS STUDIES IN THE SNS ACCUMULATOR RING
The SNS Ring off-line parallel simulation environment based on the Unified Accelerator Libraries (UAL) has been implemented and used for extensive full-scale beam dynamics studies arising in high-intensity rings. The paper describes the structure of this environment and its application to the development and analysis of the SNS accumulator ring beam loss model including a complex combination of several physical effects
Atomic oxygen beam source for erosion simulation
Cuthbertson, J. W.; Langer, W. D.; Motley, R. W.
1990-01-01
A device for production of low-energy (5-10 eV) neutral atomic beams for surface modification studies, which recreates the flux of atomic oxygen in LEO, is described. The beam is produced by acceleration of plasma ions onto a negatively biased plate of high-Z metal; the ions are neutralized and reflected by the surface, retaining a large fraction of their incident kinetic energy, forming a beam of atoms. The device is based on a magnetically confined (3-4 kG) coaxial plasma source and the atom energy can be varied by adjusting the bias voltage. The source provides a neutral flux of roughly 5 x 10 to the 16th/sq cm/s at a distance of 10 cm and a fluence of roughly 10 to the 21st/sq cm in five hours. The source has been characterized with plasma diagnostics and by measuring the energy of an atomic argon beam using a mass spectrometer. Samples of carbon film, carbon-based paint, Kapton, Mylar, and Teflon exposed to atomic O beams show erosion quite similar to those observed in orbit on the Space Shuttle.
Beam dynamics study of RFQ for CADS with a 3D space-charge-effect
The ADS (accelerator driven subcritical system) project was proposed by the Chinese Academy of Sciences. The initial proton beams delivered from an electron cyclotron resonance ion source can be effectively accelerated by 162.5 MHz 4.2 m long room temperature radio-frequency-quadrupoles (RFQ) operating in CW mode. To test the feasibility of this physical design, a new Fortran code for RFQ beam dynamics study, which is space charge dominated, was developed. This program is based on Particle-In-Cell (PIC) technique in the time domain. Using the RFQ structure designed for the CADS project, the beam dynamics behavior is performed. The well-known simulation code TRACK is used for benchmarks. The results given by these two codes show good agreements. Numerical techniques as well as the results of beam dynamics studies are presented in this paper. (authors)
DTL cavity design and beam dynamics for a TAC linear proton accelerator
Caliskan, A.; Yılmaz, M.
2012-02-01
A 30 mA drift tube linac (DTL) accelerator has been designed using SUPERFISH code in the energy range of 3-55 MeV in the framework of the Turkish Accelerator Center (TAC) project. Optimization criteria in cavity design are effective shunt impedance (ZTT), transit-time factor and electrical breakdown limit. In geometrical optimization we have aimed to increase the energy gain in each RF gap of the DTL cells by maximizing the effective shunt impedance (ZTT) and the transit-time factor. Beam dynamics studies of the DTL accelerator have been performed using beam dynamics simulation codes of PATH and PARMILA. The results of both codes have been compared. In the beam dynamical studies, the rms values of beam emittance have been taken into account and a low emittance growth in both x and y directions has been attempted.
DTL cavity design and beam dynamics for a TAC linear proton accelerator
A 30 mA drift tube linac (DTL) accelerator has been designed using SUPERFISH code in the energy range of 3-55 MeV in the framework of the Turkish Accelerator Center (TAC) project. Optimization criteria in cavity design are effective shunt impedance (ZTT), transit-time factor and electrical breakdown limit. In geometrical optimization we have aimed to increase the energy gain in each RF gap of the DTL cells by maximizing the effective shunt impedance (ZTT) and the transit-time factor. Beam dynamics studies of the DTL accelerator have been performed using beam dynamics simulation codes of PATH and PARMILA. The results of both codes have been compared. In the beam dynamical studies, the rms values of beam emittance have been taken into account and a low emittance growth in both x and y directions has been attempted. (authors)
Beam dynamics study of RFQ for CADS with a 3D space-charge-effect
Li, Chao; Zhang, Zhi-Lei; Qi, Xin; Xu, Xian-Bo; He, Yuan; Yang, Lei
2014-03-01
The ADS (accelerator driven subcritical system) project was proposed by the Chinese Academy of Sciences. The initial proton beams delivered from an electron cyclotron resonance ion source can be effectively accelerated by 162.5 MHz 4.2 m long room temperature radio-frequency-quadrupoles (RFQ) operating in CW mode. To test the feasibility of this physical design, a new Fortran code for RFQ beam dynamics study, which is space charge dominated, was developed. This program is based on Particle-In-Cell (PIC) technique in the time domain. Using the RFQ structure designed for the CADS project, the beam dynamics behavior is performed. The well-known simulation code TRACK is used for benchmarks. The results given by these two codes show good agreements. Numerical techniques as well as the results of beam dynamics studies are presented in this paper.
Beam-Beam Simulation of Crab Cavity White Noise for LHC Upgrade
Qiang, J; Pieloni, Tatiana; Ohmi, Kazuhito
2015-01-01
High luminosity LHC upgrade will improve the luminosity of the current LHC operation by an order of magnitude. Crab cavity as a critical component for compensating luminosity loss from large crossing angle collision and also providing luminosity leveling for the LHC upgrade is being actively pursued. In this paper, we will report on the study of potential effects of the crab cavity white noise errors on the beam luminosity lifetime based on strong-strong beam-beam simulations.
Dynamic cone beam CT angiography of carotid and cerebral arteries using canine model
Cai Weixing; Zhao Binghui; Conover, David; Liu Jiangkun; Ning Ruola [Department of Imaging Sciences, University of Rochester, 601 Elmwood Avenue, Rochester, New York 14642 (United States); Department of Radiology, Shanghai 6th People' s Hospital, 600 Yishan Road, Xuhui, Shanghai (China); Koning Corporation, Lennox Tech Enterprise Center, 150 Lucius Gordon Drive Suite 112, West Henrietta, New York 14586 (United States); Department of Imaging Sciences, University of Rochester, 601 Elmwood Avenue, Rochester, New York 14642 (United States); Department of Imaging Sciences, University of Rochester, 601 Elmwood Avenue, Rochester, New York 14642 (United States) and Koning Corporation, Lennox Tech Enterprise Center, 150 Lucius Gordon Drive Suite 112, West Henrietta, New York 14586 (United States)
2012-01-15
Purpose: This research is designed to develop and evaluate a flat-panel detector-based dynamic cone beam CT system for dynamic angiography imaging, which is able to provide both dynamic functional information and dynamic anatomic information from one multirevolution cone beam CT scan. Methods: A dynamic cone beam CT scan acquired projections over four revolutions within a time window of 40 s after contrast agent injection through a femoral vein to cover the entire wash-in and wash-out phases. A dynamic cone beam CT reconstruction algorithm was utilized and a novel recovery method was developed to correct the time-enhancement curve of contrast flow. From the same data set, both projection-based subtraction and reconstruction-based subtraction approaches were utilized and compared to remove the background tissues and visualize the 3D vascular structure to provide the dynamic anatomic information. Results: Through computer simulations, the new recovery algorithm for dynamic time-enhancement curves was optimized and showed excellent accuracy to recover the actual contrast flow. Canine model experiments also indicated that the recovered time-enhancement curves from dynamic cone beam CT imaging agreed well with that of an IV-digital subtraction angiography (DSA) study. The dynamic vascular structures reconstructed using both projection-based subtraction and reconstruction-based subtraction were almost identical as the differences between them were comparable to the background noise level. At the enhancement peak, all the major carotid and cerebral arteries and the Circle of Willis could be clearly observed. Conclusions: The proposed dynamic cone beam CT approach can accurately recover the actual contrast flow, and dynamic anatomic imaging can be obtained with high isotropic 3D resolution. This approach is promising for diagnosis and treatment planning of vascular diseases and strokes.
Dynamic cone beam CT angiography of carotid and cerebral arteries using canine model
Purpose: This research is designed to develop and evaluate a flat-panel detector-based dynamic cone beam CT system for dynamic angiography imaging, which is able to provide both dynamic functional information and dynamic anatomic information from one multirevolution cone beam CT scan. Methods: A dynamic cone beam CT scan acquired projections over four revolutions within a time window of 40 s after contrast agent injection through a femoral vein to cover the entire wash-in and wash-out phases. A dynamic cone beam CT reconstruction algorithm was utilized and a novel recovery method was developed to correct the time-enhancement curve of contrast flow. From the same data set, both projection-based subtraction and reconstruction-based subtraction approaches were utilized and compared to remove the background tissues and visualize the 3D vascular structure to provide the dynamic anatomic information. Results: Through computer simulations, the new recovery algorithm for dynamic time-enhancement curves was optimized and showed excellent accuracy to recover the actual contrast flow. Canine model experiments also indicated that the recovered time-enhancement curves from dynamic cone beam CT imaging agreed well with that of an IV-digital subtraction angiography (DSA) study. The dynamic vascular structures reconstructed using both projection-based subtraction and reconstruction-based subtraction were almost identical as the differences between them were comparable to the background noise level. At the enhancement peak, all the major carotid and cerebral arteries and the Circle of Willis could be clearly observed. Conclusions: The proposed dynamic cone beam CT approach can accurately recover the actual contrast flow, and dynamic anatomic imaging can be obtained with high isotropic 3D resolution. This approach is promising for diagnosis and treatment planning of vascular diseases and strokes.
Nonlinear δf Simulations of Collective Effects in Intense Charged Particle Beams
A nonlinear delta(f) particle simulation method based on the Vlasov-Maxwell equations has been recently developed to study collective processes in high-intensity beams, where space-charge and magnetic self-field effects play a critical role in determining the nonlinear beam dynamics. Implemented in the Beam Equilibrium, Stability and Transport (BEST) code [H. Qin, R.C. Davidson, and W.W. Lee, Physical Review -- Special Topics on Accelerator and Beams 3 (2000) 084401; 3 (2000) 109901.], the nonlinear delta(f) method provides a low-noise and self-consistent tool for simulating collective interactions and nonlinear dynamics of high-intensity beams in modern and next-generation accelerators and storage rings, such as the Spallation Neutron Source and heavy ion fusion drivers. A wide range of linear eigenmodes of high-intensity charged-particle beams can be systematically studied using the BEST code. Simulation results for the electron-proton two-stream instability in the Proton Storage Ring experiment [R. Macek, et al., in Proc. of the Particle Accelerator Conference, Chicago, 2001 (IEEE, Piscataway, NJ, 2001), Vol. 1, p. 688.] at the Los Alamos National Laboratory agree well with experimental observations. Large-scale parallel simulations have also been carried out for the ion-electron two-stream instability in the very-high-intensity heavy ion beams envisioned for heavy ion fusion applications. In both cases, the simulation results indicate that the dominant two-stream instability has a dipole-mode (hose-like) structure and can be stabilized by a modest axial momentum spread of the beam particles
Beam dynamics with new booster dipoles
New bending magnets are being designed for the booster synchrotron at RRCAT, Indore with the same effective length and field which will be installed in the existing ring with the same configuration of drifts and quadrupole magnets. Presently sector type dipoles are in use. It is easier to fabricate parallel edge (rectangular type) dipoles but the beam optics gets modified due to edges which provide additional focusing. The effect on tune point can be corrected using two quadrupole families. Studies indicate that the beam emittance is lower in the optics with rectangular type dipoles but the beam injection and extraction are more difficult. In this paper, the beam optics, beam emittance, injection and extraction with two configurations of the dipole magnets are compared. (author)
Nonlinear delta f Simulations of Collective Effects in Intense Charged Particle Beams
Hong Qi
2003-01-01
A nonlinear delta(f) particle simulation method based on the Vlasov-Maxwell equations has been recently developed to study collective processes in high-intensity beams, where space-charge and magnetic self-field effects play a critical role in determining the nonlinear beam dynamics. Implemented in the Beam Equilibrium, Stability and Transport (BEST) code [H. Qin, R.C. Davidson, and W.W. Lee, Physical Review -- Special Topics on Accelerator and Beams 3 (2000) 084401; 3 (2000) 109901.], the nonlinear delta(f) method provides a low-noise and self-consistent tool for simulating collective interactions and nonlinear dynamics of high-intensity beams in modern and next-generation accelerators and storage rings, such as the Spallation Neutron Source and heavy ion fusion drivers. A wide range of linear eigenmodes of high-intensity charged-particle beams can be systematically studied using the BEST code. Simulation results for the electron-proton two-stream instability in the Proton Storage Ring experiment [R. Macek, ...
Enabling Strain Hardening Simulations with Dislocation Dynamics
Arsenlis, A; Cai, W
2006-12-20
Numerical algorithms for discrete dislocation dynamics simulations are investigated for the purpose of enabling strain hardening simulations of single crystals on massively parallel computers. The algorithms investigated include the /(N) calculation of forces, the equations of motion, time integration, adaptive mesh refinement, the treatment of dislocation core reactions, and the dynamic distribution of work on parallel computers. A simulation integrating all of these algorithmic elements using the Parallel Dislocation Simulator (ParaDiS) code is performed to understand their behavior in concert, and evaluate the overall numerical performance of dislocation dynamics simulations and their ability to accumulate percents of plastic strain.
Mastorides, T; Rivetta, C.; Fox, J.D.; Winkle, D.Van; /SLAC; Tytelman, D.; /Dimtel, Redwood City
2008-07-07
A dynamics simulation model is used to estimate limits of performance of the Positron-Electron Project (PEP-II). The simulation captures the dynamics and technical limitations of the Low Level Radio Frequency (LLRF) system, the high-power RF components and the low-order mode coupled bunch longitudinal beam dynamics. Simulation results showing the effect of non-linearities on the LLRF loops, and studies of the effectiveness of technical component upgrades are reported, as well as a comparison of these results with PEP-II measurements. These studies have led to the estimation of limits and determining factors in the maximum stored current that the Low Energy Ring/High Energy Ring (LER/HER) can achieve, based on system stability for different RF station configurations and upgrades. In particular, the feasibility of the PEP-II plans to achieve the final goal in luminosity, which required an increase of the beam currents to 4A for LER and 2.2A for HER, is studied. These currents are challenging in part because they would push the longitudinal low-order beam mode stability to the limit, and the klystron forward power past a level of satisfactory margin. An acceptable margin is defined in this paper, which in turn determines the corresponding klystron forward power limitation.
Simulation of ion beam scattering in a gas stripper
Maxeiner, Sascha; Suter, Martin; Christl, Marcus; Synal, Hans-Arno
2015-10-01
Ion beam scattering in the gas stripper of an accelerator mass spectrometer (AMS) enlarges the beam phase space and broadens its energy distribution. As the size of the injected beam depends on the acceleration voltage through phase space compression, the stripper becomes a limiting factor of the overall system transmission especially for low energy AMS system in the sub MV region. The spatial beam broadening and collisions with the accelerator tube walls are a possible source for machine background and energy loss fluctuations influence the mass resolution and thus isotope separation. To investigate the physical processes responsible for these effects, a computer simulation approach was chosen. Monte Carlo simulation methods are applied to simulate elastic two body scattering processes in screened Coulomb potentials in a (gas) stripper and formulas are derived to correctly determine random collision parameters and free path lengths for arbitrary (and non-homogeneous) gas densities. A simple parametric form for the underlying scattering cross sections is discussed which features important scaling behaviors. An implementation of the simulation was able to correctly model the data gained with the TANDY AMS system at ETH Zurich. The experiment covered transmission measurements of uranium ions in helium and beam profile measurements after the ion beam passed through the He-stripper. Beam profiles measured up to very high stripper densities could be understood in full system simulations including the relevant ion optics. The presented model therefore simulates the fundamental physics of the interaction between an ion beam and a gas stripper reliably. It provides a powerful and flexible tool for optimizing existing AMS stripper geometries and for designing new, state of the art low energy AMS systems.
By performing full Particle-In-Cell simulations, we examined the transient response of electrons released for the charge neutralization of a local ion beam emitted from an ion engine which is one of the electric propulsion systems. In the vicinity of the engine, the mixing process of electrons in the ion beam region is not so obvious because of large difference of dynamics between electrons and ions. A heavy ion beam emitted from a spacecraft propagates away from the engine and forms a positive potential region with respect to the background. Meanwhile electrons emitted for a neutralizer located near the ion engine are electrically attracted or accelerated to the core of the ion beam. Some electrons with the energy lower than the ion beam potential are trapped in the beam region and move along with the ion beam propagation with a multi-streaming structure in the beam potential region. Since the locations of the neutralizer and the ion beam exit are different, the above-mentioned bouncing motion of electrons is also observed in the direction of the beam diameter
Strength and behavior in shear of reinforced concrete deep beams under dynamic loading conditions
Highlights: ► Effects of wider range of loading rates on dynamic shear behavior of RC deep beams. ► Experimental investigation of RC deep beam with and without shear reinforcements. ► Verification of experimental results with truss model and FE simulation results. ► Empirical equations are proposed to predict the dynamic increase factor of maximum resistance. -- Abstract: Research on reinforced concrete (RC) deep beams has seen considerable headway over the past three decades; however, information on the dynamic shear strength and behavior of RC deep beams under varying rates of loads remains limited. This paper describes the experimental results of 24 RC deep beams with and without shear reinforcements under varying rates of concentrated loading. Results obtained serve as useful data on shear resistance, failure patterns and strain rates corresponding to varying loading rates. An analytical truss model approach proves its efficacy in predicting the dynamic shear resistance under varying loading rates. Furthermore, three-dimensional nonlinear finite element (FE) model is described and the simulation results are verified with the experimental results. A parametric study is then conducted to investigate the influence of longitudinal reinforcement ratio, transverse reinforcement ratio and shear span to effective depth ratio on shear behavior. Subsequently, two empirical equations were proposed by integrating the various parameters to assess the dynamic increase factor (DIF) of maximum resistance under varying rates of concentrated loading
A Particle In Cell code development for high current ion beam transport and plasma simulations
Joshi, N
2016-01-01
A simulation package employing a Particle in Cell (PIC) method is developed to study the high current beam transport and the dynamics of plasmas. This package includes subroutines those are suited for various planned projects at University of Frankfurt. In the framework of the storage ring project (F8SR) the code was written to describe the beam optics in toroidal magnetic fields. It is used to design an injection system for a ring with closed magnetic field lines. The generalized numerical model, in Cartesian coordinates is used to describe the intense ion beam transport through the chopper system in the low energy beam section of the FRANZ project. Especially for the chopper system, the Poisson equation is implemented with irregular geometries. The Particle In Cell model is further upgraded with a Monte Carlo Collision subroutine for simulation of plasma in the volume type ion source.
Numerical simulation program of multicomponent ion beam transport from ECR ion source
MA Lei; SONG Ming-Tao; CAO Yun; ZHAO Hong-Wei; ZHANG Zi-Min; LI Xue-Qian; LI Jia-Cai
2004-01-01
In order to research multi-component ion beam transport process and improve transport efficiency, a special simulating program for ECR beam is becoming more and more necessary. We have developed a program written by Visual Basic to be dedicated to numerical simulation of the highly charged ion beam and to optimization of beam dynamics in transport line. In the program the exchange of electrons between highly charged ions and low chargedions or neutral atoms (residual gas in transport line) is taken into account, adopting classical molecular over-barrier model and Monte Carlo method, so the code can easily give the change of charge state distribution along the transmission line. The main advantage of the code is the ability to simultaneously simulate a large quantity of ions with different masses and charge states, and particularly, to simulate the loss of highly charged ions and the increase of low charged ions due to electron exchange in the whole transport process. Some simulations have been done to study the transmission line of LECR3[1] which is an ECR ion source for highly charged ion beam at IMP. Compared with experimental results, the simulations are considered to be successful.
Three-dimensional particle simulation of heavy-ion fusion beams
The beams in a heavy-ion-beam-driven inertial fusion (HIF) accelerator are collisionless, nonneutral plasmas, confined by applied magnetic and electric fields. These space-charge-dominated beams must be focused onto small (few mm) spots at the fusion target, and so preservation of a small emittance is crucial. The nonlinear beam self-fields can lead to emittance growth, and so a self-consistent field description is needed. To this end, a multidimensional particle simulation code, WARP [Friedman et al., Part. Accel. 37-38, 131 (1992)], has been developed and is being used to study the transport of HIF beams. The code's three-dimensional (3-D) package combines features of an accelerator code and a particle-in-cell plasma simulation. Novel techniques allow it to follow beams through many accelerator elements over long distances and around bends. This paper first outlines the algorithms employed in WARP. A number of applications and corresponding results are then presented. These applications include studies of: beam drift-compression in a misaligned lattice of quadrupole focusing magnets; beam equilibria, and the approach to equilibrium; and the MBE-4 experiment [AIP Conference Proceedings 152 (AIP, New York, 1986), p. 145] recently concluded at Lawrence Berkeley Laboratory (LBL). Finally, 3-D simulations of bent-beam dynamics relevant to the planned Induction Linac Systems Experiments (ILSE) [Fessenden, Nucl. Instrum. Methods Plasma Res. A 278, 13 (1989)] at LBL are described. Axially cold beams are observed to exhibit little or no root-mean-square emittance growth at midpulse in transiting a (sharp) bend. Axially hot beams, in contrast, do exhibit some emittance growth
Simulation and Measurements of Beam Losses on LHC Collimators During Beam Abort Failures
Lari, L; Bruce, R; Goddard, B; Redaelli, S; Salvachua, B; Valentino, G; Faus-Golfe, A
2013-01-01
One of the main purposes of tracking simulations for collimation studies is to produce loss maps along the LHC ring, in order to identify the level of local beam losses during nominal and abnormal operation scenarios. The SixTrack program is the standard tracking tool used at CERN to perform these studies. Recently, it was expanded in order to evaluate the proton load on different collimators in case of fast beam failures. Simulations are compared with beam measurements at 4 TeV. Combined failures are assumed which provide worst-case scenarios of the load on tungsten tertiary collimators.
Online external beam radiation treatment simulator
Radiation therapy is an effective and widely accepted form of treatment for many types of cancer that requires extensive computerized planning. Unfortunately, current treatment planning systems have limited or no visual aid that combines patient volumetric models extracted from patient-specific CT data with the treatment device geometry in a 3D interactive simulation. We illustrate the potential of 3D simulation in radiation therapy with a web-based interactive system that combines novel standards and technologies. We discuss related research efforts in this area and present in detail several components of the simulator. An objective assessment of the accuracy of the simulator and a usability study prove the potential of such a system for simulation and training. (orig.)
Numerical Simulation Multicomponent Ion Beam Transport form ECR Ion Source
MaLei; SongMingtao; ZhangZimin; CaoYun
2003-01-01
In order to simulate the transport of multi-components ion beam extracted from an ECR ion source, we have developed a multi-charged ion beam transport program named MCIBS 1.0. The program is dedicated to numerical simulation of the behavior of highly-charged ion beam and optimization of beam optics in transport lines and is realized on a PC with Windows user interface of Microsoft Visual Basic. Among all the ions with different charge states in the beam, the exchanges of electrons between highly charged ions and low charged ions or neutral,atoms of residual gas are taken into account by using classical Molecular Over-barrier Model and Monte Carlo method. An advanced Windows graphical interface makes it; comfortable and friendly for the user to operate in an interactive mode. The present program is used for the numerical calculation and optimization of beam optics in a transport line consisting of various magnetic elements, such as dipole magnet, quadrupole and so on. It is possible to simultaneously simulate 200,000 particles, in a transport line of 340 m at most, and show every particle orbit. Beam cross section graphics and emittance phase pictures can be also shown at any position in the transport line.
Nonlinear hybrid simulation of internal kink with beam ion effects in DIII-D
Shen, Wei; Sheng, Zheng-Mao [Department of Physics, Institute for Fusion Theory and Simulation, Zhejiang University, Hangzhou 310027 (China); Fu, G. Y.; Tobias, Benjamin [Princeton Plasma Physics Laboratory, Princeton, New Jersey 08543 (United States); Zeeland, Michael Van [General Atomics, San Diego, California 92186-5608 (United States); Wang, Feng [School of Physics and Optoelectronic Engineering, Dalian University of Technology, Dalian 116024 (China)
2015-04-15
In DIII-D sawteething plasmas, long-lived (1,1) kink modes are often observed between sawtooth crashes. The saturated kink modes have two distinct frequencies. The mode with higher frequency transits to a fishbone-like mode with sufficient on-axis neutral beam power. In this work, hybrid simulations with the global kinetic-magnetohydrodynamic (MHD) hybrid code M3D-K have been carried out to investigate the linear stability and nonlinear dynamics of the n = 1 mode with effects of energetic beam ions for a typical DIII-D discharge where both saturated kink mode and fishbone were observed. Linear simulation results show that the n = 1 internal kink mode is unstable in MHD limit. However, with kinetic effects of beam ions, a fishbone-like mode is excited with mode frequency about a few kHz depending on beam pressure profile. The mode frequency is higher at higher beam power and/or narrower radial profile consistent with the experimental observation. Nonlinear simulations have been performed to investigate mode saturation as well as energetic particle transport. The nonlinear MHD simulations show that the unstable kink mode becomes a saturated kink mode after a sawtooth crash. With beam ion effects, the fishbone-like mode can also transit to a saturated kink mode with a small but finite mode frequency. These results are consistent with the experimental observation of saturated kink mode between sawtooth crashes.
Molecular dynamics simulation by atomic mass weighting
Mao, Boryeu; Friedman, Alan R.
1990-01-01
A molecular dynamics-based simulation method in which atomic masses are weighted is described. Results from this method showed that the capability for conformation search in molecular dynamics simulation of a short peptide (FMRF-amide) is significantly increased by mass weighting.
Mathematical simulation of precision scanning by a H- beam
A method and a computer program are developed for computing the dynamics of precision beams. The numerical method used for integrating differential equations is conservative. Within a paraxial approximation, the method ensures strict conservation of the phase volume at every step of calculations. Geometrical aberrations up to the third order are included and the space charge is taken into consideration within the paraxial approximation. The method is valid for relativistic beams. Controlling fields can have arbitrary longitudinal distributions. An accompanying frame of reference with a moving base is used where a trajectory of one of the particles is chosen to be a curvilinear axis. If a beam axis does not coincide with that of the device through which the beam travels, or if the angle of the beam axis deviation from that of the device (e.g. a scanner) considerably exceeds the angular divergence, within the given accuracy the method allows to be restricted with calculations of aberrations of lower order than in the case of a fixed frame of reference with a rectilinear axis. The method ensures simultaneous calculations of dynamics of any number of particles. For the case when a set consists of a large number of particles, a gain in the calculation time is obtained in comparison with times for calculation of the dynamics of all the individual particles. Dynamics of a broad H- beam having a small angular divergence at the inlet is studied in the scanner field. In order to keep the angular divergence at a small level at the outlet either distortion of the angular divergence of the beam should not be allowed or appearing angular distortions should be incessantly suppressed. Distributions of fields acting upon the beam while scanning are found when a small angular divergence at the outlet is provided
Multi-Bunch Beam Dynamics Studies for the European XFEL
Baboi, N
2004-01-01
In the X-ray free electron laser planned to be built at DESY (TESLA XFEL) the acceleration of the electron bunches will be made with 9-cell superconducting cavities. These cavities have been initially developed within the TESLA linear collider study. The impact of the higher order modes (HOM) has been shown to be within the acceptable beam dynamics limits for the collider. For the XFEL the dynamics is relaxed from point of view of multi-bunch effects (e.g. shorter length, higher emittance). However the lower energy and different time structure of the beam make the study of the HOM effects in the XFEL linac necessary. Multi-bunch beam dynamics studies are ongoing. The results of the HOM measurements at the TESLA Test Facility are used. Several options for the beam structure, as necessary for various applications, are studied. The results will be discussed.
Temporal nonlinear beam dynamics in infiltrated photonic crystal fibers
Bennet, Francis; Rosberg, Christian Romer; Neshev, Dragomir N.;
-sensing as well as active devices for all-optical switching at low (mW) laser powers. Commercially available PCFs infiltrated with liquids also provide a versatile and compact tool for exploration of the fundamentals of nonlinear beam propagation in periodic photonic structures. To explore the full...... scientific and technological potential of liquid-infiltrated PCFs it is important to understand the temporal dynamics of nonlinear beam propagation in such structures. In this work we consider thermally induced spatial nonlinear effects in infiltrated photonic crystal fibers. We experimentally study the...... temporal dynamics of nonlinear beam reshaping occurring on a short time scale before the establishment of a steady state regime. In experiment, a 532nm laser beam can be injected into a single hole of an infiltrated PCF cladding structure, and the temporal dynamics of the nonlinear response is measured by...
Transverse and longitudinal beam dynamics studies at the Fermilab photoinjector
The Fermilab photoinjector produces electron bunches of 1-12 nC charge with an energy of 16-18 MeV. Detailed measurements and optimization of the transverse emittance have been carried out for a number of beam line optics conditions, and at a number of beam line locations. The length of the bunches has also been measured, first for an uncompressed beam (as a function of the charge) and then for a compressed beam of 8 nC charge (as a function of the 9-cell cavity phase). These measurements are presented and compared with the simulation codes HOMDYN and ASTRA
Numerical simulation and analytical study of glulam timber beams
Themistoklis Tsalkatidis
2014-03-01
Full Text Available Glulam beams or glued-laminated beams consist of sawn lumber laminations (timber bonded with an adhesive material. This paper, through the mathematical description of the contact conditions that apply at the interfaces of glulam beams and the development of two three-dimensional finite element models by the use of the ANSYS software package, studies the flexural properties of unreinforced (UGB and reinforced (RGB glulam beams. The first computational model presents an unreinforced glulam beam that has been produced by three wood laminations of dimensions 6 by 3.6 by 176 cm. The latter one describes a reinforced glulam beam, which has been produced by gluing a 0.15 cm thick steel plate at the bottom edge of the previously described beam. The computational analysis indicates that the two glulam beams have significantly different bearing capacities under the same load and support conditions. The failure mode of the UGB is brittle whereas the one of the RGB is ductile. The numerical results of both models are in close agreement with experimental ones from the international literature. Keywords: Glulam Timber Beams, Numerical Simulation, Contact.
Beam Dynamics Challenges for FCC-ee
AUTHOR|(SzGeCERN)442987; Benedikt, Michael; Oide, Katsunobu; Bogomyagkov, Anton; Levichev, Evgeny; Migliorati, Mauro; Wienands, Uli
2015-01-01
The goals of FCC-ee include reaching luminosities of up to a few 1036 cm-2s-1 per interaction point at the Z pole or some 1034 cm-2s-1 at the ZH production peak, and pushing the beam energy up to ≥175 GeV, in a ring of 100 km circumference, with a total synchrotron-radiation power not exceeding 100 MW. A parameter baseline as well as high-luminosity crab-waist options were described in [1] and [2], respectively. The extremely high luminosity and resulting short beam lifetime (due to radiative Bhabha scattering) are sustained by top-up injection. The FCC-ee design status and typical beam parameters for different modes of operation are reported in [3]. One distinct feature of the FCC-ee design is its conception as a double ring, with separate beam pipes for the two counter-rotating (electron and positron) beams, resembling, in this aspect, the high-luminosity B factories PEP-II, KEKB and SuperKEKB as well as the LHC. The two separate rings do not only permit operation with a large number of bunches, up to a f...
Computer simulations of a low energy proton beam tomograph
This work presents the recent development of a low energy proton beam tomograph. The proton tomograph prototype (involving UTFPR, UERJ, UFRJ and IEN/CNEN) has been installed and tested at the cyclotron CV-28 of IEN/CNEN. New computer simulations were performed in order to optimize the performance of the scattered proton beam and its aluminum collimator energy losses. The computer code simulates the tomographic measurements with two aluminum collimators (variable aperture from 0.2 mm to 0.4 mm in diameter and variable thickness from 4 mm to 8 mm), a water phantom and a Si(Li) detector. The analysis of the exit beam energy spectra in comparison with a perfectly collimated proton beam made it possible to achieve the best quality of reconstructed tomographic images of water phantom. (author)
Computer simulations of a low energy proton beam tomograph
Milhoretto, E.; Schelin, H.R.; Setti, J.A.P.; Denyak, V.; Paschuk, S.A.; Basilio, A.C.; Rocha, R.; Ribeiro Junior, S. [Universidade Tecnologica Federal do Parana (UTFPR), Curitiba, PR (Brazil). Curso de Pos-Graduacao em Engenharia Eletrica e Informatica Industrial (CPGEI)]. E-mails: sergei@utfpr.edu.br; edneymilhoretto@yahoo.com; schelin@cpgei.cefetpr.br; Evseev, I.; Yevseyeva, O. [Universidade Estadual do Rio de Janeiro (UERJ), Nova Friburgo, RJ (Brazil)]. E-mail: evseev@iprj.uerj.br; Lopes, R.T. [Universidade Federal do Rio de Janeiro (UFRJ), Rio de Janeiro, RJ (Brazil). Coordenacao dos Programas de Pos-graducao em Engenharia (COPPE). Lab. de Instrumentacao Nuclear]. E-mail: ricardo@lin.ufrj.br; Vinagre Filho, U.M. [Instituto de Energia Nuclear (IEN/CNEN-RJ), Rio de Janeiro, RJ (Brazil)
2007-07-01
This work presents the recent development of a low energy proton beam tomograph. The proton tomograph prototype (involving UTFPR, UERJ, UFRJ and IEN/CNEN) has been installed and tested at the cyclotron CV-28 of IEN/CNEN. New computer simulations were performed in order to optimize the performance of the scattered proton beam and its aluminum collimator energy losses. The computer code simulates the tomographic measurements with two aluminum collimators (variable aperture from 0.2 mm to 0.4 mm in diameter and variable thickness from 4 mm to 8 mm), a water phantom and a Si(Li) detector. The analysis of the exit beam energy spectra in comparison with a perfectly collimated proton beam made it possible to achieve the best quality of reconstructed tomographic images of water phantom. (author)
Measurements and simulations of focused beam for orthovoltage therapy
Purpose: Megavoltage photon beams are typically used for therapy because of their skin-sparing effect. However, a focused low-energy x-ray beam would also be skin sparing, and would have a higher dose concentration at the focal spot. Such a beam can be produced with polycapillary optics. MCNP5 was used to model dose profiles for a scanned focused beam, using measured beam parameters. The potential of low energy focused x-ray beams for radiation therapy was assessed. Methods: A polycapillary optic was used to focus the x-ray beam from a tungsten source. The optic was characterized and measurements were performed at 50 kV. PMMA blocks of varying thicknesses were placed between optic and the focal spot to observe any variation in the focusing of the beam after passing through the tissue-equivalent material. The measured energy spectrum was used to model the focused beam in MCNP5. A source card (SDEF) in MCNP5 was used to simulate the converging x-ray beam. Dose calculations were performed inside a breast tissue phantom. Results: The measured focal spot size for the polycapillary optic was 0.2 mm with a depth of field of 5 mm. The measured focal spot remained unchanged through 40 mm of phantom thickness. The calculated depth dose curve inside the breast tissue showed a dose peak several centimeters below the skin with a sharp dose fall off around the focus. The percent dose falls below 10% within 5 mm of the focus. It was shown that rotating the optic during scanning would preserve the skin-sparing effect of the focused beam. Conclusions: Low energy focused x-ray beams could be used to irradiate tumors inside soft tissue within 5 cm of the surface
Beam equipment electromagnetic interaction in accelerators: simulation and experimental benchmarking
Passarelli, Andrea; Vaccaro, Vittorio Giorgio; Massa, Rita; Masullo, Maria Rosaria
One of the most significant technological problems to achieve the nominal performances in the Large Hadron Collider (LHC) concerns the system of collimation of particle beams. The use of collimators crystals, exploiting the channeling effect on extracted beam, has been experimentally demonstrated. The first part of this thesis is about the optimization of UA9 goniometer at CERN, this device used for beam collimation will replace a part of the vacuum chamber. The optimization process, however, requires the calculation of the coupling impedance between the circulating beam and this structure in order to define the threshold of admissible intensity to do not trigger instability processes. Simulations have been performed with electromagnetic codes to evaluate the coupling impedance and to assess the beam-structure interaction. The results clearly showed that the most concerned resonance frequencies are due solely to the open cavity to the compartment of the motors and position sensors considering the crystal in o...
Absolute nodal coordinate plane beam formulation for multibody systems dynamics
A new plane beam dynamic formulation for constrained multibody system dynamics is developed. Flexible multibody system dynamics includes rigid body dynamics and superimposed vibratory motions. The complexity of mechanical system dynamics originates from rotational kinematics, but the natural coordinate formulation does not use rotational coordinates, so that simple dynamic formulation is possible. These methods use only translational coordinates and simple algebraic constraints. A new formulation for plane flexible multibody systems are developed utilizing the curvature of a beam and point masses. Using absolute nodal coordinates, a constant mass matrix is obtained and the elastic force becomes a nonlinear function of the nodal coordinates. In this formulation, no infinitesimal or finite rotation assumptions are used and no assumption on the magnitude of the element rotations is made. The distributed body mass and applied forces are lumped to the point masses. Closed loop mechanical systems consisting of elastic beams can be modeled without constraints since the loop closure constraints can be substituted as beam longitudinal elasticity. A curved beam is modeled automatically. Several numerical examples are presented to show the effectiveness of this method.
Absolute nodal coordinate plane beam formulation for multibody systems dynamics
Souh, Byungyil, E-mail: bysouh@dyu.ac.kr [Dongyang University (Korea, Republic of)
2013-06-15
A new plane beam dynamic formulation for constrained multibody system dynamics is developed. Flexible multibody system dynamics includes rigid body dynamics and superimposed vibratory motions. The complexity of mechanical system dynamics originates from rotational kinematics, but the natural coordinate formulation does not use rotational coordinates, so that simple dynamic formulation is possible. These methods use only translational coordinates and simple algebraic constraints. A new formulation for plane flexible multibody systems are developed utilizing the curvature of a beam and point masses. Using absolute nodal coordinates, a constant mass matrix is obtained and the elastic force becomes a nonlinear function of the nodal coordinates. In this formulation, no infinitesimal or finite rotation assumptions are used and no assumption on the magnitude of the element rotations is made. The distributed body mass and applied forces are lumped to the point masses. Closed loop mechanical systems consisting of elastic beams can be modeled without constraints since the loop closure constraints can be substituted as beam longitudinal elasticity. A curved beam is modeled automatically. Several numerical examples are presented to show the effectiveness of this method.
Study of transient self-consistent beam dynamics in RF linacs using a particle tracing code
Mytrochenko, V. V.; Opanasenko, A.
2006-03-01
The paper describes a simulation technique for study of unsteady self-consistent dynamics of charged particles in RF linacs that consist of cavities and travelling wave sections. The approach proposed is based on unsteady theories of the excitation of cavities and waveguides by a beam of charged particles and RF feeders. The theory of waveguide excitation is extended to the case of spatially inhomogeneous travelling-wave structures. The SUPERFISH code is used to evaluate the characteristics of the axisymmetric travelling-wave sections. The PARMELA code is applied for simulation of particle motion and to obtain data required for solving the equations for excitation of the RF structures by the beam.