Full Paper
Contents  Return  Previous Page  Next Page  Index

Return To: Session TH/P3 - Edge Transport & Stability, Waves
Prev Page: (TH/P3-20) Variations of the Magnetic Moment of Fast
Next Page: (TH/P3-22) Selfconsistent RF Driven and Bootstrap Currents

(TH/P3-21) Tera-Scale Computation of Wave-Plasma Interactions in Multidimensional Fusion Plasmas

D.B. Batchelor1), L.A. Berry1), M.D. Carter1), C.K. Phillips2), R.J. Dumont2), P.T. Bonoli3), J.C. Wright3), D.N. Smithe4), R.W. Harvey5), D.A. D’Ippolito6), J.R. Myra6), E. D’Azevedo7)
1) Fusion Energy Division, Oak Ridge National Laboratory, Oak Ridge, Tennessee, USA
2) Princeton Plasma Physics Laboratory, Princeton, New Jersey, USA
3) Massachusetts Institute of Technology, Cambridge, Massachusetts, USA
4) Mission Research Corporation, Newington, Virginia, USA
5) CompX, Del Mar, California, USA
6) Lodestart Research Corporation, Boulder, Colorado, USA
7) Computer Science and Mathematics Division, Oak Ridge National Laboratory, Oak Ridge, Tennessee, USA

Abstract.  With support from the Scientific Discovery Through Advanced Computation (SciDAC) program, we have established a multi-institutional partnership between plasma physicists and computational scientists, with the overarching goal to obtain quantitatively accurate predictive understanding of electromagnetic wave processes important for heating, current drive, stability and transport applications in fusion-relevant plasmas. Activities during the first year of the project have focussed on massive parallelization and acceleration of computer-intensive full-wave RF field solver codes, extension of all-orders methods to two and three-dimensional plasmas, benchmarking and code comparison, and application to wave propagation problems in Alcator C-Mod, NSTX and LHD. The TORIC code has been modified to use an out-of-core linear solver thereby reducing memory requirements by a factor of 25 and processing time by a factor of 3. Calculation of D(3 He) mode conversion in Alcator C-Mod have been carried out with up to 256 poloidal modes, whereas it was previously infeasible to use as many as 100 modes. The IBW wave fields appear to be converged with this number of modes and now exhibit the proper localization to the high field side. Also the partitioning of p

Read the full paper in PDF format.

IAEA 2003