Alex's Application


Heavy Ion Fusion beam simulations

We are working to develop heavy-ion particle accelerators as drivers for Inertial Fusion Energy power production. Plasma Particle-In-Cell (PIC) methods have proven very effective for studies of sections of such accelerators, but integrated full-machine simulations have been beyond our reach. One of our goals is to transform an existing accelerator code, WARP3d, into a qualitatively improved tool, which will be efficient on new-generation computers with non-uniform memory access. The need for a self-consistent field makes this challenging. One or more elliptic or hyperbolic PDE?s must be solved on a 3-D computational grid, using source terms which must be accumulated from the particles. This is nontrivial in a distributed-computing environment. Load-balancing may require adaptive domain decomposition. Another goal is to simulate the beam's passage from the driver through the fusion chamber and onto the fusion target; this involves interaction with background plasma, and a fully electromagnetic (PIC) description is required. Finally, we would like to develop continuum (phase-space fluid) methods for solving the Vlasov equation deterministically on a 4-D or 5-D grid.

A solid knowledge of available tools, especially ACTS tools, will be of great benefit to this research effort. We are especially interested in remote steering and visualization (in particular, CUMULVS) using either AVS or OPENDX. Also, as Python users (our principal accelerator code, WARP, uses a Python interface) we are interested in SILOON (Scripting Interface Languages for Object-Oriented Numerics) which supports Python and for which FORTRAN support may be forthcoming.