[Table of Contents]

Plasma and Fusion Research

Volume 5, S1028 (2010)

Regular Articles

Practical Calculation of Nuclear Fusion Power for a Toroidal Plasma Device with Magnetic Confinement
St. Petersburg Polytechnic University, 195251, Russia
(Received 10 January 2009 / Accepted 19 August 2009 / Published 26 March 2010)


An algorithm has been developed and realized as a FORTRAN code to calculate the volume integral power of a magnetic confinement nuclear fusion reactor and the local fusion rate function. It takes the nuclei energy distributions, fusion cross-sections, and magnetic surface geometry as input data. Two fast, simple analytic models of magnetic flux surfaces were used, and the corresponding Jacobian determinants were found. Gaussian kernel empirical probability density estimation has been proposed to reconstruct the ion energy probability density function from experimentally obtained random samples of escaping neutral atom energies. The influence of ion heating and fast ion confinement on high-energy distribution tails, and thus on the neutron yield and fusion power, can be calculated. The code has been applied to obtain radial profiles of the nuclear fusion reaction rate and volume integral power for both Maxwellian and suprathermal D and T particle distributions. A fast neutral particle diagnostic database may serve as a basis for an experimentally confirmed calculation technique for reactor power and ignition criterion.


fusion reactor power, nuclear fusion rate, ion distribution function, high-energy particle, non-Maxwellian rate coefficient, empirical probability density

DOI: 10.1585/pfr.5.S1028


  • [1] G.H. Miley, H. Towner and N. Ivich, Rept. COO-2218-17, University of Illinois (1974).
  • [2] B.H. Duane, Rept. BNWL-1685, Brookhaven National Laboratory (1972).
  • [3] H.-S. Bosch and G.M. Hale, Nucl. Fusion 32, 611 (1992).
  • [4] L.E. Zakharov, Preprint IAE-4114/6, Moscow (1985).
  • [5] R.L. Miller, M.S. Chu, J.M. Greene, Y.R. Lin-Liu and R.E. Waltz, Phys. Plasmas 5, 973 (1998).
  • [6] J.G. Cordey and M.J. Houghton, Nucl. Fusion 13, 215 (1973).
  • [7] P.R. Goncharov, T. Ozaki et al., Rev. Sci. Instrum. 79, 10F311 (2008).
  • [8] P.R. Goncharov, T. Ozaki et al., Rev. Sci. Instrum. 79, 10F312 (2008).
  • [9] P.R. Goncharov, T. Ozaki et al., Plasma Fusion Res. 3, S1083 (2008).
  • [10] M. Matsumoto and T. Nishimura, ACM Trans. Model. Comput. Simul. 8, 3 (1998).
  • [11] S. Murakami et al., A Demonstration of Magnetic Field Optimization in LHD, IAEA-CN-94/ EX/C5-3, 19th IAEA Fusion Energy Conference, Lyon, France, October, 2002.
  • [12] C. Gormezano, Plasma Phys. Control. Fusion 41, B367 (1999).
  • [13] R.J. Hawryluk et al., Phys. Plasmas 5, 1577 (1998).
  • [14] M. Matsukawa et al., Engineering Feature in the Design of JT-60SA, FT/P7-5, 21st IAEA Fusion Energy Conference, Chengdu, China, October 2006.
  • [15] M. Shimada et al., eds., Nucl. Fusion 47, S1 (2007).

This paper may be cited as follows:

Pavel R. GONCHAROV, Plasma Fusion Res. 5, S1028 (2010).