Author: Samulyak, V.
Paper Title Page
MOPBA06 Algorithms and Self-consistent Simulations of Beam-induced Plasma in Muon Cooling Devices 186
 
  • V. Samulyak
    BNL, Upton, Long Island, New York, USA
  • M. Chung, A.V. Tollestrup, K. Yonehara
    Fermilab, Batavia, USA
  • R.D. Ryne
    LBNL, Berkeley, California, USA
  • K. Yu
    SBU, Stony Brook, USA
 
  Funding: Research is partially supported by the DOE MAP program
Interaction of muon beams with plasma generated in muon cooling absorbers is an important issue affecting the efficiency of muon cooling. We have developed numerical algorithms and parallel software for self-consistent simulation of the plasma production and its interaction with particle beams and external electromagnetic fields. Simulations support the FNAL experimental program on dense hydrogen gas filled RF cavities proposed for muon beam phase space cooling and acceleration. The core code uses the particle-in-cell (PIC) method for the Maxwell equations coupled to the dynamics of particles. Electromagnetic PIC methods are combined with probabilistic treatment of atomic physics processes responsible for the plasma production. The PIC code supports the dynamics of multiple particle species undergoing rapid acceleration / deceleration (variable relativistic factor) and uses accurate charge and current conservation methods and symplectic discretization schemes. It is fully parallel and runs on multicore supercomputers. Benchmarks and simulations of experiments on gas-filled RF cavities will be discussed.
 
 
TUPBA09 Simulation of High Power Mercury Jet Targets for Neutrino Factory and Muon Collider 541
 
  • V. Samulyak, H.G. Kirk
    BNL, Upton, Long Island, New York, USA
  • H.C. Chen
    SBU, Stony Brook, USA
  • K.T. McDonald
    PU, Princeton, New Jersey, USA
 
  Funding: Research has been supported by the DOE MAP project.
Hydrodynamic behavior of high power targets for future particle accelerators and, in particular, for the proposed Neutrino Factory and Muon Collider has been investigated via numerical simulations. The target will contain a series of mercury jet pulses of about 1 cm in diameter, interacting with strong proton pulses in 15 – 20 T magnetic fields. Simulations used the Lagrangian particle / smooth particle hydrodynamics code designed to accurately resolve free surface 3D hydrodynamic flows. Simulation results have been compared with existing experimental data and previous simulations performed with the Front tracking code FronTier. Both codes use realistic equations of state for mercury that support cavitation under critical tension. They are parallelized and optimized for the use on large distributed memory supercomputers. Simulations include the power range of neutrino factory (4MW, 8 GeV proton beam containing 150 bunches per second), the muon collider (the same power proton beam but delivered in 15 bunches per second), and proton beams exceeding the muon collider power range.