Optimized tomography methods for plasma emissivity reconstruction at the ASDEX Upgrade tokamak.

The soft X-ray (SXR) emission provides valuable insight into processes happening inside of high-temperature plasmas. A standard method for deriving the local emissivity profiles of the plasma from the line-of-sight integrals measured by pinhole cameras is the tomographic inversion. Such an inversion is challenging due to its ill-conditioned nature and because the reconstructed profiles depend not only on the quality of the measurements but also on the inversion algorithm used. This paper provides a detailed description of several tomography algorithms, which solve the inversion problem of Tikhonov regularization with linear computational complexity in the number of basis functions. The feasibility of combining these methods with the minimum Fisher information regularization is demonstrated, and various statistical methods for the optimal choice of the regularization parameter are investigated with emphasis on their reliability and robustness. Finally, the accuracy and the capability of the methods are demonstrated by reconstructions of experimental SXR profiles, featuring poloidal asymmetric impurity distributions as measured at the ASDEX Upgrade tokamak.

[1]  A. Tikhonov,et al.  Use of the regularization method in non-linear problems , 1965 .

[2]  H. Akaike.  A new look at the statistical model identification , 1974 .

[3]  G. Wahba.  Practical Approximate Solutions to Linear Operator Equations When the Data are Noisy , 1977 .

[4]  S. Jardin,et al.  An iterative metric method for solving the inverse tokamak equilibrium problem , 1979 .

[5]  C. Loan.  Computing the CS and the generalized singular value decompositions , 1985 .

[6]  D. Girard.  Optimal Regularized Reconstruction in Computerized Tomography , 1987 .

[7]  T. Kitagawa.  A deterministic approach to optimal regularization —The finite dimensional case— , 1987 .

[8]  M. Bertero,et al.  Linear inverse problems with discrete data: II. Stability and regularisation , 1988 .

[9]  R. Granetz,et al.  X-ray tomography on JET , 1988 .

[10]  T. Tsukishima,et al.  Phillips–Tikhonov regularization of plasma image reconstruction with the generalized cross validation , 1989 .

[11]  R. Korde,et al.  The effect of neutron irradiation on silicon photodiodes , 1989 .

[12]  Clifford M. Hurvich,et al.  Regression and time series model selection in small samples , 1989 .

[13]  Per Christian Hansen,et al.  Analysis of Discrete Ill-Posed Problems by Means of the L-Curve , 1992, SIAM Rev..

[14]  P. Hansen.  Numerical tools for analysis and solution of Fredholm integral equations of the first kind , 1992 .

[15]  H. Weisen,et al.  Relative Calibration of Photodiodes in the Soft-X-Ray Spectral Range , 1995 .

[16]  P. J. Paris,et al.  X-ray tomography on the TCV tokamak , 1996 .

[17]  W. Linden,et al.  Maximum entropy based reconstruction of soft X-ray emissivity profiles in W7-AS , 1996 .

[18]  Yohsuke Hosoda,et al.  Linear algebraic algorithms for high speed and stable reconstruction of plasma image , 1997 .

[19]  J. Wesson.  Poloidal distribution of impurities in a rotating tokamak plasma , 1997 .

[20]  G. Fehmers,et al.  An algorithme for quadratic optimization with one quadratic constraint and bounds on the variables , 1998 .

[21]  F. Scholze,et al.  Mean energy required to produce an electron-hole pair in silicon for photons of energies between 50 and 1500 eV , 1998 .

[22]  Clifford M. Hurvich,et al.  Smoothing parameter selection in nonparametric regression using an improved Akaike information criterion , 1998 .

[23]  L. C. Ingesson,et al.  Soft X ray tomography during ELMs and impurity injection in JET , 1998 .

[24]  C. Görner.  Tomographische Untersuchung von globalen Alfven-Eigenmoden am Stellarator Wendelstein 7-AS , 1998 .

[25]  Chao Yang,et al.  ARPACK users' guide - solution of large-scale eigenvalue problems with implicitly restarted Arnoldi methods , 1998, Software, environments, tools.

[26]  P. McCarthy,et al.  The CLISTE Interpretive Equilibrium Code , 1999 .

[27]  R. Neu,et al.  Z dependence of the core impurity transport in ASDEX Upgrade H mode discharges , 1999 .

[28]  G. Janeschitz.  Plasma–wall interaction issues in ITER , 2001 .

[29]  J. Lister,et al.  Investigation of the consistency of magnetic and soft X-ray plasma position measurements on TCV by means of a rapid tomographic inversion algorithm , 2003 .

[30]  Per Christian Hansen,et al.  REGULARIZATION TOOLS: A Matlab package for analysis and solution of discrete ill-posed problems , 1994, Numerical Algorithms.

[31]  R. Neu,et al.  Tungsten: an option for divertor and main chamber plasma facing components in future fusion devices , 2005 .

[32]  K. Tobita,et al.  Case study on tritium inventory in the fusion DEMO plant at JAERI , 2006 .

[33]  Gene H. Golub,et al.  Generalized cross-validation as a method for choosing a good ridge parameter , 1979, Milestones in Matrix Computation.

[34]  L. C. Ingesson,et al.  Chapter 7: Tomography Diagnostics: Bolometry and Soft-X-Ray Detection , 2008 .

[35]  YANQING CHEN,et al.  Algorithm 8 xx : CHOLMOD , supernodal sparse Cholesky factorization and update / downdate ∗ , 2006 .

[36]  J. Stober,et al.  Ten years of W programme in ASDEX Upgrade—challenges and conclusions , 2009 .

[37]  M. Maraschek,et al.  Hotlink based Soft X-ray Diagnostic on ASDEX Upgrade , 2010 .

[38]  R. Neu,et al.  Local effects of ECRH on argon transport in L-mode discharges at ASDEX Upgrade , 2011 .

[39]  R. Neu,et al.  Calculation and experimental test of the cooling factor of tungsten , 2010 .

[40]  A. Hubbard,et al.  Poloidal variation of high-Z impurity density due to hydrogen minority ion cyclotron resonance heating on Alcator C-Mod , 2012 .

[41]  D. Moreau,et al.  Soft x-ray tomography for real-time applications: present status at Tore Supra and possible future developments. , 2012, The Review of scientific instruments.

[42]  A. Murari,et al.  Modern numerical methods for plasma tomography optimisation , 2012 .

[43]  D. Mazon,et al.  Absolute spectral characterization of silicon barrier diode: Application to soft X-ray fusion diagnostics at Tore Supra , 2013 .

[44]  Soft X-ray imaging techniques on Tore Supra: Present status and possible future developments , 2013 .

[45]  D. Clayton,et al.  Fast Nickel and Iron Density Estimation Using Soft X-Ray Measurements in Tore Supra: Preliminary Study , 2013 .

[46]  J. Rice,et al.  Non-neoclassical up/down asymmetry of impurity emission on Alcator C-Mod , 2013 .

[47]  Investigation of 3D tungsten distributions in (1,1) kink modes induced by toroidal plasma rotation , 2015 .

[48]  A. Gude.  Identification of MHD Instabilities in Experiments , 2015 .

[49]  Judith C. Hill,et al.  Regularization of soft-X-ray imaging in the DIII-D tokamak , 2015, J. Comput. Phys..

[50]  V. Kiptily,et al.  Phillips-Tikhonov regularization with a priori information for neutron emission tomographic reconstruction on Joint European Torus. , 2015, The Review of scientific instruments.

[51]  Martin Imrisek,et al.  Minimum Fisher Tikhonov Regularization Adapted to Real-Time Tomography , 2016 .

[52]  A. Gude,et al.  Non-monotonic growth rates of sawtooth precursors evidenced with a new method on ASDEX Upgrade , 2016 .