GE FULL PROBLEM (with diff. num. particles)


The following runs use the original XPDC1 cross-sectional data. All the runs were conducted on the Cray T3E MPP machine in batch mode using the parallel version of XPDC1. We used 16 processors for the 400 K particle runs. All the other runs used 64 processors. The simulation starts off with initial Bessel function electron and ion density profiles with axial density = 1.0e11/cc. The number of particles in the simulation is varied by changing the parameter nc2p (ratio of number of real particles vs. computer particles). For 400 K particles, nc2p = 5.0e7. The following is a link to the input file for 400 K particles: Input File.
The simulations were run for 3.6 microseconds.




The following is the result of doubling the time step in the 400 K particle run. This simulation was run for 2.2 microsecs. Doubling the time step does not reduce the heating at the axis. Instead the axial peak growth rate increases.