METHOD FOR MEASURING PLASMA ION NONEXTENSIVE PARAMETER

Abstract

The present invention relates to a method for measuring the ion nonextensive parameter of plasma includes the following steps: describe the plasma with nonextensive statistical mechanics, obtain the equation describing the relationship between the geodesic acoustic mode frequency and the ion acoustic speed of plasma; collect the measurement data of the geodesic acoustic mode frequencies and plasma temperature in the device where the plasma is to be measured; the obtained equation describing the relationship between the geodesic acoustic mode frequency and the ion acoustic speed of plasma is used to linearly fit the collected measured data of the geodesic acoustic mode frequency and the plasma temperature in the device where the plasma is to be measured to obtain the slope value; based on the derived equation and the obtained slope values, and combining with the safety factor of the device where the plasma is to be measured, the ion nonextensive parameter is solved numerically. The present invention fills the gap where the electron nonextensive parameter can be measured with the nonextensive single electric probe, but the corresponding ion nonextensive parameter cannot be diagnosed yet in the field of nonextensive parameters diagnosis.

Claims

1. A method for measuring an ion nonextensive parameter of a plasma, comprises: step 1, describing the plasma with nonextensive statistical mechanics to obtain an equation describing the relationship between the geodesic acoustic mode frequency and the ion acoustic speed of plasma; step 2, collecting the measurement data of the geodesic acoustic mode frequencies and plasma temperature in the device where the plasma is to be measured; step 3, using the equation describing the relationship between the geodesic acoustic mode frequency and the ion acoustic speed of plasma derived in step 1 to obtain the slope value by linearly fitting the measurement data of the geodesic acoustic mode frequency and plasma temperature collected in step 2 in the device where the plasma to be measured; step 4, according to the formula obtained in step 1 and the slope value obtained in step 3, combined with the safety factor of the device where the plasma is to be measured to obtain the ion nonextensive parameter by numerical solution.

2. The method according to claim 1, wherein a formula of describing the relationship between the geodesic acoustic mode frequency and the ion nonextensive parameter of plasma in step 1 is: $f_{\mathrm{GAM}}=\sqrt{S\left(q_{F_i},q\right)}\frac{c_s}{2\pi R_0}=\sqrt{\frac{7}{4\left(3q_{F_i}-1\right)}\left(1+\sqrt{1+\frac{4}{q}\frac{\frac{23}{2\left(5q_{F_i}-3\right)\left(3q_{F_i}-1\right)}}{{\left[-\frac{7}{2\left(3q_{F_i}-1\right)}\right]}^2}}\right)}\frac{v_{\mathrm{ti}}}{2\pi R_0}$ where v.sub.ti is ion effective thermal speed, c.sub.s is ion acoustic speed, q.sub.F.sub.i is ion nonextensive parameter, R.sub.0 is the major radius of tokamak plasma, q is the safety factor of the device.

3. The method according to claim 1, wherein the measuring the ion nonextensive parameter of plasma, comprises: step 5, drawing a SSE−q.sub.F.sub.i curve using a sum of squares due to error (SSE) as a goodness of fit measure; step 6, obtaining a value of the ion nonextensive parameter corresponding to a minimum value of the SSE according to the SSE−q.sub.F.sub.i, curve; step 7, drawing a R.sup.2−q.sub.F.sub.i, curve using a nonlinear fitting coefficient of determination R.sup.2 as a goodness of fit measure; step 8, obtaining a value of the ion nonextensive parameter corresponding to a maximum value of the R.sup.2 according to the R.sup.2-qq.sub.F.sub.i □ curve; step 9, comparing the values of the ion nonextensive parameter respectively corresponding to the minimum value of the SSE obtained in the step 6 and the maximum value of the R.sup.2 obtained in the step 8; step 10, determining whether the values of the ion nonextensive parameter compared in the step 9 are consistent; step 11, confirming any one of the values of the ion nonextensive parameter as an optimal value of the ion nonextensive parameter when the values of the ion nonextensive parameter are determined to be consistent; step 12, substituting the optimal value of the ion nonextensive parameter obtained in the step 11 into the data set in the step 3 and step 4 to obtain a target data set corresponding to the optimal value of the ion nonextensive parameter; step 13, outputting a measurement result report.

4. The method according to claim 1, wherein the measuring the ion nonextensive parameter of plasma, comprises: in step 2, the device where the plasma to be measured is T-10 tokamak.

Description

BRIEF DESCRIPTION OF DRAWINGS

[0026] FIG. 1 is a schematic diagram of the frequency measurement method HIBP and nonextensive electric probe used in the T-10 tokamak device of the present invention;

[0027] FIG. 2 is the experimental data analysis of the geodesic acoustic model on the T-10 tokamak device of the present invention;

[0028] FIG. 3 is a measurement diagram of ion nonextensive parameter of plasma generated by 36815 shot on T-10 tokamak device in the present invention;

[0029] FIG. 4 is the analysis graphs of statistic SSE for the optimal ion nonextensive parameters in the present invention;

[0030] FIG. 5 is the analysis graphs of statistic R.sup.2 for the optimal ion nonextensive parameters in the present invention.

[0031] In the picture, 1 is ion gun, 2 is defection plate, 3 is primary beam, 4 is secondary beam, 5 is primary beam detector, 6 is electrostatic analyzer.

DETAILED DESCRIPTION OF EMBODIMENTS

[0032] The invention is described in detail below in conjunction with the attached drawings.

[0033] See FIGS. 1-5.

[0034] The invention provides a method for measuring ion nonextensive parameter of plasma. In order to obtain the geodesic acoustic mode theory under the nonextensive statistical framework consistent with the experiment, we extend the geodesic acoustic mode theory under the Boltzmann-Gibbs statistical framework to the theory under the nonextensive statistical framework. The obtained geodesic acoustic mode dispersion relationship under the nonextensive statistical framework is as follows:

[00002]$\begin{array}{cc}{f}_{\mathrm{GAM}}=\frac{{\omega }_{\mathrm{GAM}}}{2\pi }=\sqrt{S\left(q_{F_i},q\right)}\frac{c_s}{2\pi R_0}& \left(1\right)\end{array}$$\mathrm{where}$$\begin{array}{cc}S\left(q_{F_i},q\right)=\frac{7}{4\left(3q_{F_i}-1\right)}\left(1+\sqrt{1+\frac{4}{q^2}\frac{\frac{23}{2\left(5q_{F_i}-3\right)\left(3q_{F_i}-1\right)}}{{\left(-\frac{7}{2\left(3q_{F_i}-1\right)}\right)}^2}}\right),\left(q_{F_i}>\frac{1}{3}\right)& \left(2\right)\end{array}$

[0035] The above formula is clearly illustrated as follows: the geodesic acoustic mode frequency is proportional to the ion acoustic speed, which is supported by the fluid and kinetic theory and experimental data, while what is different from the theory under the Boltzmann-Gibbs (extensive) statistical framework is that the proportional coefficient is not only a function of the safety factor, but also a function of nonextensive parameters: the proportional coefficient decreases with the increase of the safety factor, and also decreases with the increase of the ion nonextensive parameter; substitute Eq. (2) into Eq. (1), and c.sub.s approximately equal to the v.sub.ti, the following formula is obtained:

[00003]$f_{\mathrm{GAM}}=\sqrt{S\left(q_{F_i},q\right)}\frac{c_s}{2\pi R_0}=\sqrt{\frac{7}{4\left(3q_{F_i}-1\right)}\left(1+\sqrt{1+\frac{4}{q^2}\frac{\frac{23}{2\left(5q_{F_i}-3\right)\left(3q_{F_i}-1\right)}}{{\left[-\frac{7}{2\left(3q_{F_i}-1\right)}\right]}^2}}\right)}\frac{v_{\mathrm{ti}}}{2\pi R_0}$

[0036] This indicates that the geodesic acoustic mode f.sub.GAM−c.sub.s/2πR.sub.0 curve (the theoretical cornerstone of ion nonextensive parameter diagnosis) has a kind of complicated dependence on nonextensive parameters, which is different from the traditional (excluding nonextensive parameters) geodesic acoustic mode theory; in addition, we found that at the extensive limit (q.sub.F.sub.i=1), the above results all return to the traditional theory based on Boltzmann-Gibbs statistical framework, which supports the correctness and universality of the nonextensive theory (namely a larger scope of application).

[0037] The above analysis has shown that nonextensive parameters have an influence on the geodesic acoustic mode f.sub.GAM−c.sub.s/2πR.sub.0 curve; based on this theory, next, we will explain how to measure ion nonextensive parameters that cannot be measured even with a nonextensive single electric probe.

[0038] As shown in FIG. 2, in order to measure ion nonextensive parameter, we can first measure the geodesic acoustic mode frequency and electron temperature of the plasma to be measured in a particular tokamak device by the present method, the present embodiment adopts the frequency measurement method HIBP and nonextensive electric probe already used by the T-10 tokamak device as shown in FIG. 1. In this invention, 4 experimental data points were obtained for the plasma generated by 36815 discharges of the T-10 tokamak device.

[0039] As shown in FIG. 3, Then the least square method is used to fit this set of experimental data. Since both the theory and the experiment have proved the relationship between f.sub.GAM and c.sub.s/2πR.sub.0 is the direct proportional, we make a linear fitting of the direct proportional function without intercept. It turns out that the optimal slope is √{square root over (S)}=1.003, as shown in Table 1 below.

TABLE-US-00001 TABLE 1 Parameters related to plasma generated by 36815 shot on T-10 Tokamak device. √{square root over (S)} SSE R.sup.2 q q.sub.Fi 1.003 0.284 0.992 3.3 1.565

[0040] Then, according to equation (2), the relationship between the safety factor q and the ion nonextensive parameter q.sub.F.sub.i is known. Since q=3.3 is the safety factor of 36815 shot on T-10 device, the corresponding ion nonextensive parameter q.sub.F.sub.i=1.565 can be solved. The four curves are the q.sub.F.sub.i−q figure given by Eq. (2) combining with the slop √{square root over (S)}=1.003 obtained from figure f.sub.GAM−c.sub.s/2πR.sub.0. As you can see from the diagram, there are four groups of solutions, and here, from a physical point of view, we only consider the group of solutions represented by the red real curve. q.sub.F.sub.i=1.565 is the ion nonextensive parameter of plasma generated by 36815 shot (safety factor q=3.3) on T-10 device.

[0041] As shown in FIG. 4, In order to illustrate that q.sub.F.sub.i=1.565 is the optimal ion nonextensive parameter, statistics SSE and R.sup.2 are analyzed. We found that when the ion nonextensive parameter q.sub.F takes different values, SSE also takes different values and gets the minimum value when q.sub.F.sub.i=1.565, indicating that q.sub.F.sub.i=1.565 is the optimal ion nonextensive parameter. In order to confirm the reliability of the results measured with the statistic SSE.

[0042] As shown in FIG. 5, we also analyzed another independent indicator R.sup.2, and made an R.sup.2−q.sub.F.sub.i graph, and also found that when q.sub.F.sub.i=1.565, R.sup.2 achieves the maximum value, which confirms q.sub.F.sub.i=1.565 is the optimal ion nonextensive parameter, that is, the measurement result of the ion nonextensive parameter of the plasma generated by 36815 shot on T-10 tokamak device is q.sub.F.sub.i=1.565.

[0043] The results show the effectiveness of a method for measuring ion nonextensive parameters in a tokamak device. Recent studies have shown that replacing Boltzmann-Gibbs statistical mechanics with nonextensive statistical mechanics has a strong advantage in describing plasmas. The diagnosis error without using nonextensive statistical mechanics to describe the plasma may be as high as 83.91%. We established the theory of nonextensive geodesic acoustic mode by introducing nonextensive statistical mechanics to take into account the system nonextensivity which has been proven by a large number of facts. This theory not only can obtain the related results of the traditional geodesic acoustic mode at the extensive limit, which proves the correctness of the nonextensive geodesic acoustic mode theory, but also can measure the ion nonextensive parameter (q.sub.F.sub.i=1.565, see FIG. 3) that cannot be measured even with a nonextensive single electric probe by combining the existing diagnosis methods of geodesic acoustic mode frequency and plasma electron temperature (such as nonextensive electric probe).

[0044] The present invention fills the gap where the electron nonextensive parameter can be measured with the nonextensive single electric probe, but the corresponding ion nonextensive parameter cannot be diagnosed yet in the field of nonextensive parameters diagnosis.

[0045] The above is only an embodiment of the invention, does not therefore limit the scope of the patent of the invention, all equivalent transformations made by the contents of the description and drawings of the invention or direct or indirect application in related technical fields, are included in the patent protection scope of the invention.