128 Matching Results

Search Results

Advanced search parameters have been applied.

Three dimensional simulations of space charge dominated heavy ion beams with applications to inertial fusion energy

Description: Heavy ion fusion requires injection, transport and acceleration of high current beams. Detailed simulation of such beams requires fully self-consistent space charge fields and three dimensions. WARP3D, developed for this purpose, is a particle-in-cell plasma simulation code optimized to work within the framework of an accelerator`s lattice of accelerating, focusing, and bending elements. The code has been used to study several test problems and for simulations and design of experiments. Two applications are drift compression experiments on the MBE-4 facility at LBL and design of the electrostatic quadrupole injector for the proposed ILSE facility. With aggressive drift compression on MBE-4, anomalous emittance growth was observed. Simulations carried out to examine possible causes showed that essentially all the emittance growth is result of external forces on the beam and not of internal beam space-charge fields. Dominant external forces are the dodecapole component of focusing fields, the image forces on the surrounding pipe and conductors, and the octopole fields that result from the structure of the quadrupole focusing elements. Goal of the design of the electrostatic quadrupole injector is to produce a beam of as low emittance as possible. The simulations show that the dominant effects that increase the emittance are the nonlinear octopole fields and the energy effect (fields in the axial direction that are off-axis). Injectors were designed that minimized the beam envelope in order to reduce the effect of the nonlinear fields. Alterations to the quadrupole structure that reduce the nonlinear fields further were examined. Comparisons were done with a scaled experiment resulted in very good agreement.
Date: November 1, 1994
Creator: Grote, D.P.
Partner: UNT Libraries Government Documents Department

Novel Methods in the Particle-In-Cell Accelerator Code-Framework Warp

Description: The Particle-In-Cell (PIC) Code-Framework Warp is being developed by the Heavy Ion Fusion Science Virtual National Laboratory (HIFS-VNL) to guide the development of accelerators that can deliver beams suitable for high energy density experiments and implosion of inertial fusion capsules. It is also applied in various areas outside the Heavy Ion Fusion program to the study and design of existing and next-generation high-energy accelerators, including the study of electron cloud effects and laser wakefield acceleration for example. This paper presents an overview of Warp’s capabilities, summarizing recent original numerical methods that were developed by the HIFS-VNL (including Particle-In-Cell with Adaptive Mesh Refinement, a large-timestep “drift-Lorentz” mover for arbitrarily magnetized species, a relativistic Lorentz invariant leapfrog particle pusher, simulations in Lorentz boosted frames, an electromagnetic solver with tunable numerical dispersion and efficient stride20 based digital filtering), with great emphasis on the description of the mesh refinement capability. Selected examples of applications of the methods to the abovementioned fields are given.
Date: September 1, 2011
Creator: Vay, J.-L.; Grote, D. P.; Cohen, R. H.; Friedman, A.; Grote, D. P.; Cohen, R. H. et al.
Partner: UNT Libraries Government Documents Department

Large-timestep mover for particle simulations of arbitrarilymagnetized species

Description: For self-consistent ion-beam simulations including electron motion, it is desirable to be able to follow electron dynamics accurately without being constrained by the electron cyclotron timescale. To this end, we have developed a particle-advance that interpolates between full particle dynamics and drift motion. By making a proper choice of interpolation parameter, simulation particles experience physically correct parallel dynamics, drift motion, and gyroradius when the timestep is large compared to the cyclotron period, though the effective gyro frequency is artificially low; in the opposite timestep limit, the method approaches a conventional Boris particle push. By combining this scheme with a Poisson solver that includes an interpolated form of the polarization drift in the dielectric response, the movers utility can be extended to higher-density problems where the plasma frequency of the species being advanced exceeds its cyclotron frequency. We describe a series of tests of the mover and its application to simulation of electron clouds in heavy-ion accelerators.
Date: March 26, 2007
Creator: Cohen, R.H.; Friedman, A.; Grote, D.P. & Vay, J-L.
Partner: UNT Libraries Government Documents Department

An Implicit "Drift-Lorentz" Particle Mover for Plasma and Beam Simulations

Description: In order to efficiently perform particle simulations in systems with widely varying magnetization, we have developed a 'drift-Lorentz mover', which interpolates between full particle dynamics and drift kinetics in such a way as to preserve a physically correct gyroradius and particle drifts for both large and small ratios of the timestep to the cyclotron period. In order to extend applicability of the mover to systems with plasma frequency exceeding the cyclotron frequency - such as one may have with fully neutralized drift compression of a heavy ion beam - we have developed an implicit version of the mover. A first step in this direction, in which the polarization charge was added to the field solver, was described previously. Here we describe a fully implicit algorithm (which is analogous to the direct-implicit method for conventional particle-in-cell simulation), a stability analysis of it, its implementation in the WARP code, and several tests of the resultant code. The fully implemented version is electrostatic; we are beginning development of an electromagnetic version, and describe also the status of that effort.
Date: March 7, 2008
Creator: Cohen, R H; Friedman, A; Grote, D P & Vay, J
Partner: UNT Libraries Government Documents Department

Compact High-Current Heavy-Ion Injector

Description: To provide a compact high-brightness heavy-ion beam source for Heavy Ion Fusion (HIF), we have been experimenting with merging multi-beamlets in an injector which uses an RF plasma source. An array of converging beamlets was use to produce a beam with the envelope radius, convergence, and ellipticity matched to an electrostatic quadrupole (ESQ) channel. Experimental results were in good quantitative agreement with simulation and have demonstrated the feasibility of this concept. The size of a driver-scale injector system using this approach will be several times smaller than one designed using traditional single large-aperture beams. The success of this experiment has possible significant economical and technical impacts on the architecture of HIF drivers.
Date: April 13, 2006
Creator: Westenskow, G A; Grote, D P; Kwan, J W & Bieniosek, F
Partner: UNT Libraries Government Documents Department

DOD-SBIR Structured Multi-Resolution PIC Code for Electromagnetic Plasma Simulations, Final Report

Description: A novel electromagnetic solver with mesh refinement capability was implemented in Warp. The solver allows for calculations in 2-1/2 and 3 dimensions, includes the standard Yee stencil, and the Cole-Karkkainen stencil for lower numerical dispersion along the principal axes. Warp implementation of the Cole-Karkkainen stencil includes an extension to perfectly matched layers (PML) for absorption of waves, and is preserving the conservation property of charge conserving current deposition schemes, like the Buneman-Villanesor and Esirkepov methods. Warp's mesh refinement framework (originally developed for electrostatic calculations) was augmented to allow for electromagnetic capability, following the methodology presented in [1] extended to an arbitrary number of refinement levels. Other developments include a generalized particle injection method, internal conductors using stair-cased approximation, and subcycling of particle pushing. The solver runs in parallel using MPI message passing, with a choice at runtime of 1D, 2D and 3D domain decomposition, and is shown to scale linearly on a test problem up-to 32,768 CPUs. The novel solver was tested on the modeling of filamentation instability, fast ignition, ion beam induced plasma wake, and laser plasma acceleration.
Date: April 22, 2010
Creator: Vay, J L; Grote, D P & Friedman, A
Partner: UNT Libraries Government Documents Department

Large-Timestep Mover for Particle Simulations of Arbitrarily Magnetized Species

Description: For self-consistent ion-beam simulations including electron motion, it is desirable to be able to follow electron dynamics accurately without being constrained by the electron cyclotron timescale. To this end, we have developed a particle-advance that interpolates between full particle dynamics and drift motion. By making a proper choice of interpolation parameter, simulation particles experience physically correct parallel dynamics, drift motion, and gyroradius when the timestep is large compared to the cyclotron period, though the effective gyro frequency is artificially low; in the opposite timestep limit, the method approaches a conventional Boris particle push. By combining this scheme with a Poisson solver that includes an interpolated form of the polarization drift in the dielectric response, the movers utility can be extended to higher-density problems where the plasma frequency of the species being advanced exceeds its cyclotron frequency. We describe a series of tests of the mover and its application to simulation of electron clouds in heavy-ion accelerators.
Date: June 16, 2006
Creator: Cohen, R. H.; Friedman, A.; Grote, D. P. & Vay, J.
Partner: UNT Libraries Government Documents Department

Report Initial Work on Developing Plasma Modeling Capability in WARP for NDCX Experiments

Description: This milestone has been accomplished. The Heavy Ion Fusion Science Virtual National Laboratory (HIFS-VNL) has developed and implemented an initial beam-in-plasma implicit modeling capability in Warp; has carried out tests validating the behavior of the models employed; has compared the results of electrostatic and electromagnetic models when applied to beam expansion in an NDCX-I relevant regime; has compared Warp and LSP results on a problem relevant to NDCX-I; has modeled wave excitation by a rigid beam propagating through plasma; and has implemented and begun testing a more advanced implicit method that correctly captures electron drift motion even when timesteps too large to resolve the electron gyro-period are employed. The HIFS-VNL is well on its way toward having a state-of-the-art source-to-target simulation capability that will enable more effective support of ongoing experiments in the NDCX series and allow more confident planning for future ones.
Date: December 14, 2007
Creator: Friedman, A; Cohen, R H; Grote, D P & Vay, J
Partner: UNT Libraries Government Documents Department

Synthesizing a four-dimensional beam particle distribution frommultiple two-dimensional views

Description: The transverse dynamics of a nearly-monoenergetic particle beam are described by the evolution of the 4D distribution f(x,y,x',y'), where x and y are the transverse spatial coordinates and x' {triple_bond} p{sub x}/p{sub z} and y' {triple_bond} p{sub y}/p{sub z} are the corresponding momentum components divided by the longitudinal momentum component. In present-day experimental practice, such beams are often diagnosed by passing them through an axially-separated pair of slits parallel to the y axis. This selects for x and x' and integrates over y and y'. A sequence of pulses (with the slits at various x positions) yields a 2D projection of the beam phase space, f(x,x'). Another scanner might yield f(y,y') or, using crossed slits, f(x,y). The challenge is that a small set of such 2D scans does not uniquely specify f(x,y,x',y'); correlations in planes other than those measured are unknown. We have developed Monte-Carlo methods and formulated physically-motivated constraints to synthesize a ''reasonable'' set of particles having 2D projectional densities consistent with the experimental data. Such a set may be used to initialize simulations of the downstream beam. The methods and their performance on model problems are described.
Date: February 20, 2002
Creator: Friedman, A.; Grote, D. P.; Celata, C. M. & Staples, J. W.
Partner: UNT Libraries Government Documents Department

Progress in the study of mesh refinement for particle-in-cell plasma simulations and its application to heavy ion fusion

Description: The numerical simulation of the driving beams in a heavy ion fusion power plant is a challenging task, and, despite rapid progress in computer power, one must consider the use of the most advanced numerical techniques. One of the difficulties of these simulations resides in the disparity of scales in time and in space which must be resolved. When these disparities are in distinctive zones of the simulation region, a method which has proven to be effective in other areas (e.g. fluid dynamics simulations) is the Adaptive-Mesh-Refinement (AMR) technique. We follow in this article the progress accomplished in the last few months in the merging of the AMR technique with Particle-In-Cell (PIC) method. This includes a detailed modeling of the Lampel-Tiefenback solution for the one-dimensional diode using novel techniques to suppress undesirable numerical oscillations and an AMR patch to follow the head of the particle distribution. We also report new results concerning the modeling of ion sources using the axisymmetric WARPRZ-AMR prototype showing the utility of an AMR patch resolving the emitter vicinity and the beam edge.
Date: September 15, 2002
Creator: Vay, J.-L.; Friedman, A. & Grote, D.P.
Partner: UNT Libraries Government Documents Department

Effects of longitudinal space charge in beams for heavy-ion fusion

Description: A new analytic model is presented that accurately estimates the radially averaged axial component of the space-charge field of an axisymmetric heavy-ion beam in a cylindrical beam pipe. The model recovers details of the field near the beam ends that are overlooked by simpler models, and the results compare well to exact solutions of Poisson`s equation. Field values are shown for several simple beam profiles and are compared with values obtained from simpler models. The model has been implemented in the fluid/envelope code CIRCE and used to study longitudinal confinement in beams with a variety of axial profiles. The effects of errors in the longitudinal-control fields are presented.
Date: December 27, 1995
Creator: Sharp, W.M.; Friedman, A. & Grote, D.P.
Partner: UNT Libraries Government Documents Department

3D particle simulations of space-charge-dominated beams in HIF accelerator experiments

Description: The development of a high current, heavy-ion beam for inertial confinement fusion requires a detailed understanding of the behavior of the beam, including effects of the large self-fields. This necessity makes particle-in-cell (PIC) simulation the appropriate tool, and for this reason, the three-dimensional PIC/accelerator code WARP3d is being developed. WARP3d has been used extensively to study the creation and propagation of ion beams both to support experiments and for the understanding of basic beam physics. An overview of the structure of the code is presented along with a discussion of features that make the code an effective tool in the understanding of space-charge dominated beam behavior. A number of applications where WARP3d has played an important role is discussed, emphasizing the need of three-dimensional, first principles simulations. Results and comparisons with experiment are presented.
Date: May 1, 1997
Creator: Grote, D.P.; Friedman, A.; Lund, S.M. & Haber, I.
Partner: UNT Libraries Government Documents Department

Manipulation of high-current pulses for heavy-ion fusion

Description: For efficient induction-driven heavy-ion fusion, the current profile along a pulse must be modified in a non-selfsimilar manner between the accelerator and the target. In the accelerator, the pulse should have a duration of at least 50 ns in order to make efficient use of the induction cores, and the current should by nearly uniform along the pulse to minimize the aperture. In contrast, the optimal current profile on target consists of a main pulse of about 10 ns preceded by a longer low-current `foot.` This pulse-shape manipulation must be carried out at the final pulse energy (5-10 GeV for 200 amu ions) in the presence of a large nonlinear longitudinal space-charge field. A straightforward method is presented here for doing the required pulse shaping. Induction-ceU voltages are generated using idealized beam profiles both in the accelerator and on target, and they are verified and checked for error sensitivity using the fluid/envelope code CIRCE.
Date: October 28, 1996
Creator: Sharp, W.M.; Callahan, D.A.; Griedman, A. & Grote, D.P.
Partner: UNT Libraries Government Documents Department

Methods used in WARP3d, a three-dimensional PIC/accelerator code

Description: WARP-3d(1,2), a three-dimensional PIC/accelerator code, has been developed over several years and has played a major role in the design and analysis of space-charge dominated beam experiments being carried out by the heavy-ion fusion programs at LLNL and LBNL. Major features of the code will be reviewed, including: residence corrections which allow large timesteps to be taken, electrostatic field solution with subgrid scale resolution of internal conductor boundaries, and a beat beam algorithm. Emphasis will be placed on new features and capabilities of the code, which include: a port to parallel processing environments, space-charge limited injection, and the linking of runs covering different sections of an accelerator. Representative applications in which the new features and capabilities are used will be presented along with the important results.
Date: February 28, 1997
Creator: Grote, D.P.; Friedman, A. & Haber, I.
Partner: UNT Libraries Government Documents Department

New Development in WARP: Progress Toward End-to-End Simulation

Description: The development of a high current, heavy-ion beam driver for inertial confinement fusion requires a detailed understanding of the behavior of the beam, including effects of the strong self-fields. The necessity of including the self-fields of the beam makes particle-in-cell (PIC) simulation techniques ideal, and for this reason, the multi-dimensional PIC/accelerator code WARP has been developed. WARP has been used extensively to study the creation and propagation of ion beams both in experiments and for the understanding of basic beam physics. An overview of the structure of the code will be presented along with a discussion of features that make the code an effective tool in the understanding of space-charge dominated beam behavior. Much development has been done on WARP increasing its flexibility and generality. Major additions include a generalized field description, an efficient steady-state modelling technique, a transverse slice model with a bending algorithm, further improvement of the parallel processing version, and capabilities for linking to chamber transport codes. With these additions, the capability of modeling a large scale accelerator from end-to-end comes closer to reality.
Date: November 7, 1997
Creator: Grote, D.P.; Friedman, A.; Haber, I.; Fawley, W. & Luc Vay, J.
Partner: UNT Libraries Government Documents Department

Steering algorithms for a small recirculating heavy-ion accelerator

Description: Beam-steering algorithms are proposed for a small recirculating induction accelerator being built at the Lawrence Livermore National Laboratory. The principal problem is that the transverse position and velocity of the beam must be inferred from capacitive position monitors, and this determination is complicated by the limited probe resolution and by the lattice errors within steering modules. The fluid/envelope code CIRCE is used to evaluate these algorithms.
Date: November 7, 1997
Creator: Sharp, W.M.; Grote, D.P. & Hemandez, G.W.
Partner: UNT Libraries Government Documents Department

Longitudinal dynamics and stability in beams for heavy-ion fusion

Description: Successful transport of induction-driven beams for heavy-ion fusion requires careful control of the longitudinal space charge. The usual control technique is the periodic application of time-varying longitudinal electric fields, called `ears`, that on the average, balance the space-charge field. this technique is illustrated using a fluid/envelope code CIRCE, and the sensitivity of the method to errors in these ear fields is illustrated. The possibility that periodic ear fields also excite the longitudinal instability is examined.
Date: January 5, 1996
Creator: Sharp, W.M.; Callahan, D.A. & Grote, D.P.
Partner: UNT Libraries Government Documents Department

Three-dimensional simulations of high-current beams in induction accelerators with WARP3d

Description: For many issues relevant to acceleration and propagation of heavy-ion beams for inertial confinement fusion, understanding the behavior of the beam requires the self-consistent inclusion of the self-fields of the beams in multiple dimensions. For these reasons, the three-dimensional simulation code WARP3d A.Friedman was developed. The code combines the particle-in-cell plasma simulation technique with a realistic description of the elements which make up an accelerator. In this paper, the general structure of the code is reviewed and details of two ongoing applications are presented along with a discussion of simulation techniques used. The most important results of this work are presented.
Date: September 1, 1995
Creator: Grote, D.P.; Friedman, A. & Haber, I.
Partner: UNT Libraries Government Documents Department

Three dimensional simulations of a small induction recirculator accelerator

Description: A recirculating induction accelerator has potential cost advantages for a heavy-ion fusion driver. In order to explore the physics and technological issues, a small prototype recirculator is being built. The three dimensional particle-in-cell code, WARP3d, is being used in the design and analysis of the experiments. WARP3d is used to examine the behavior of the beam in the electric dipoles and in the non-linear fields associated with the accelerator lattice elements. The dipoles have focusing and fringe fields which can adversely affect the beam quality. Both single particle and full beam dynamics are examined in the dipoles using realistic geometries. Dipole plate designs which minimize the adverse effects are described. The non-linear fields associated with the permanent magnetic quadrupoles have been included in the simulations. They were found to have little effect on the quality of the beam.
Date: April 24, 1995
Creator: Grote, D.P.; Friedman, A. & Haber, I.
Partner: UNT Libraries Government Documents Department

Beam dynamics studies of the Heavy Ion Fusion Accelerator injector

Description: A driver-scale injector for the Heavy Ion Fusion Accelerator project has been built at LBL. This machine has exceeded the design goals of high voltage (> 2 MV), high current (> 0.8 A of K{sup +}) and low normalized emittance (< 1 {pi} mm-mr). The injector consists of a 750 keV diode pre-injector followed by an electrostatic quadrupole accelerator (ESQ) which provides strong (alternating gradient) focusing for the space-charge dominated beam and simultaneously accelerates the ions to 2 MeV. The fully 3-D PIC code WARP together with EGUN and POISSON were used to design the machine and analyze measurements of voltage, current and phase space distributions. A comparison between beam dynamics characteristics as measured for the injector and corresponding computer calculations will be presented.
Date: April 1, 1995
Creator: Henestroza, E.; Yu, S.S.; Eylon, S. & Grote, D.P.
Partner: UNT Libraries Government Documents Department

Effects of space charge in beams for heavy ion fusion

Description: A new analytic model is presented that accurately estimates the radially averaged axial component of the space-charge field of an axisymmetric heavy-ion beam in a cylindrical beam pipe. The model recovers details of the field near the beam ends that are overlooked by simpler models, and the results compare well to exact solutions of Poisson`s equation. Field values are shown for several simple beam profiles and are compared with values obtained from simpler models. The model has been implemented in the fluid/envelope code CIRCE and used to study longitudinal confinement in beams with a variety of axial profiles. The effects of errors in the longitudinal-control fields are presented.
Date: September 1, 1995
Creator: Sharp, W.M.; Friedman, A. & Grote, D.P.
Partner: UNT Libraries Government Documents Department

Acceleration schedules for a recirculating heavy-ion accelerator

Description: The recent development of miniature inductive adders has made it feasible to design programmable, high-repetition-rate pulsers with a substantially higher voltage than is possible using a conventional field-effect transistor architecture. Prototype pulsers using the new technology are being developed as part of a series of experiments at LLNL to test the concept of a recirculating induction accelerator. Preliminary numerical work is reported here to determine what effects the higher-voltage pulsers would have on the beam quality of the LLNL small recirculator.
Date: June 1, 1999
Creator: Grote, D. P. & Sharp, W.
Partner: UNT Libraries Government Documents Department