Development of SVD Algorithm for Turbulence Tomography.pdf

1、To realize a nuclear fusion reactor，magnetic confinement of plasma has been investigated for a long time to produce high performance (or high density and temperature) plasma. It has been found that plasma turbulence degr

2、ades the property of plasma confinement. Therefore，plasma turbulence has been one of the important problems to be studied in the field of magnetically confinement fusion. In the present research of plasma turbulence，it i

3、s necessary to observe the turbulence not only in local region of plasma but also over the entire plasma because the plasma turbulence consists of the fluctuations of various scales from ion Lamor radius to the full size

4、 of plasma cross-section. Tomography is one of possible diagnostics to observe the entire plasma and to realizelocal and fine measurements of turbulence. We have constructed a prototype of such tomography system and inst

5、alled on a linear cylindrical device，named Plasma Assenbly for Nonlinear Turbulence Analysis (PANTA) in Kyushu University. We are now developing a tomography algorithm using Singular Value Decomposition (SVD) in addition

6、 to the already developed ones，Maximum Likelihood-Expectation Method (ML-EM)，Algebraic Reconstruction Technique (ART)，F(xiàn)ourier-Bessel series and Cormack expansion. The paper reports simulation and experimental results of

7、the tomography using SVD algorithm.
　　Chapter 2 describes the mathematical base of SVD method together with ML-EM method that is used for our prototype system as a standard. In tomography，the relation between the line

8、-integrated values giand local emissionεi on i-th grid is expressed asg=hε,where h is a matrix. The tomography reconstruction can be performed if the inverse relation is obtained. The SVD is used to obtain the relation b

9、y factorizing the matrix h. More precisely saying，the emission can be obtained with a least square method to minimize the modifiedχ2-function Φ，with assuming Φ=∥g-hε∥2/M+?∥Cε∥2,whereγand∥∥represent the regularization pa

10、rameter and the norm，respectively，and C is selected to be the Laplacian matrix here. The solution is expressed as ε=(hTh+MγCTC)-1hTg，where M is the number of data. In this process the SVD is used and the regularization p

11、arameterγneeds to be determined using Generalized Cross Validation (GCV).
　　Chapter 3 presents the comparison results between newly applied tomography algorithm，SVD and already developed one，ML-EM.In order to evaluate

12、 new algorithm，SVD method，the SVD method is applied on sinograms (or a set of line-integrate data) calculated from given emission images. The reconstruction is performed under the;nbsp;condition that the detector arrays

13、are located at 16 different azimuthal angle positions with 101 channels for each array，respectively，thus，the number of observation channels is assumed to totally 1616. In comparison between SVD and MLEM，both algorithms h

14、ave sufficient ability to reconstruct overall structure and sufficient spatial resolution to resolve a fine structure represented by a sharp peak. More detailed comparison is made in evaluating the residual inference err

15、or.
　　Chapter 4 presents the results of the application of SVD method on the real experimental data. A prototype system of turbulence tomography is explained with the linear cylindrical device on which the prototype i

16、s installed. The cylindrical plasma has the diameter of approximately 100 mm and the length of 4 m. In the tomography system，the plasma emission is measured with four light-guide arrays installed at azimuthal angle，0°，45

17、°，90°，and 135°. Each light-guide is equipped with four sets of detector arrays，each of which has 33ch viewpoints lined up in every 5 mm，and covers the entire plasma. In the experiment the target is argon plasma，and the e

18、missions of ArI and ArII are measured simultaneously. The tomography reconstruction of SVD method is successfully performed with optimizing the regularization parameterthe GCV.
　　Finally，summary is given in chapter 5.

19、 Tomography algorithm based on SVD method is being developed，and its application to an assumed emission data well reproduces the original image. Moreover，the tomographic reconstruction on the experimental data is perform

20、ed in PANTA and the results are successfully obtained with optimizing the regularization parameterto minimize GCV. In comparison with MLEM method iteration process，the SVD method should be a faster algorithm in tomograph

21、y reconstruction of temporal series of sonogram data. This is because the SVD just repeats the calculation using the same matrix obtained in advance，while the MLEM needs to apply iterative calculations to each slice of t