CN113539525B - Method for identifying space position and structure of tearing die magnetic island in tokamak - Google Patents
Method for identifying space position and structure of tearing die magnetic island in tokamak Download PDFInfo
- Publication number
- CN113539525B CN113539525B CN202110670960.9A CN202110670960A CN113539525B CN 113539525 B CN113539525 B CN 113539525B CN 202110670960 A CN202110670960 A CN 202110670960A CN 113539525 B CN113539525 B CN 113539525B
- Authority
- CN
- China
- Prior art keywords
- electron temperature
- magnetic island
- tearing
- tokamak
- relative electron
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Classifications
-
- G—PHYSICS
- G21—NUCLEAR PHYSICS; NUCLEAR ENGINEERING
- G21B—FUSION REACTORS
- G21B1/00—Thermonuclear fusion reactors
- G21B1/05—Thermonuclear fusion reactors with magnetic or electric plasma confinement
- G21B1/057—Tokamaks
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E30/00—Energy generation of nuclear origin
- Y02E30/10—Nuclear fusion reactors
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Plasma & Fusion (AREA)
- General Engineering & Computer Science (AREA)
- High Energy & Nuclear Physics (AREA)
- Analysing Materials By The Use Of Radiation (AREA)
- Radiation Pyrometers (AREA)
Abstract
The invention belongs to a magnetic confinement nuclear fusion technology, and particularly relates to a method for identifying the space position and structure of a tearing die magnetic island in a Tokamak, wherein the polar cross section temperature of Tokamak plasma is obtained, the frequency corresponding to the maximum power is the characteristic frequency of the tearing die, the smooth relative electron temperature disturbance is determined, the rational surface position of a tearing die magnetic field is determined, the X point and the O point of the tearing die magnetic island are determined, and the space position of the tearing die magnetic island, the rotation direction of the tearing die and the polar modulus are obtained. Based on direct experimental measurement data without setting any assumption and approximate conditions, the accuracy is improved, the method has the advantages of simplicity, rapidness, intuitiveness and high efficiency, and an important basis can be provided for accurately controlling the tearing die in the Tokamak plasma.
Description
Technical Field
The invention belongs to the magnetic confinement nuclear fusion technology, and particularly relates to a method for identifying the space position and structure of a tearing die magnetic island in a tokamak.
Background
Magneto-constrained nuclear fusion is considered as a clean energy source with sufficient fuel reserves, and is one of important ways in which human energy crisis can be finally solved. There are often various instabilities in magnetically confined nuclear fusion tokamak plasmas that can severely degrade plasma confinement performance. The tearing die is one of the most threatening instabilities, and can destroy the nested magnetic surface well-constrained in the tokamak, and the magnetic force line reconnection phenomenon occurs to form a magnetic island structure (the center of the magnetic island is defined as an O point, and the reconnection point is defined as an X point), so that charged particles are rapidly transported radially from the inner side of the magnetic island to the outer side of the magnetic island along the direction of the magnetic force line. When the magnetic islands are large enough, these charged particles can directly strike the outer walls of the device causing fusion reactions to terminate and destroy the device wall material, resulting in significant economic loss. Theoretical and experimental research results on the control of the tearing die show that only the position (particularly the position of the O point of the magnetic island) and the space structure of the magnetic island are accurately determined, so that the tearing die can be effectively controlled.
Disclosure of Invention
The invention aims to provide a method for identifying the space position and the structure of a tearing die magnetic island in a tokamak, which can simply, quickly, accurately and efficiently obtain the space position and the two-dimensional structure of the magnetic island.
The technical scheme of the invention is as follows:
a method for identifying the space position and structure of a tearing die magnetic island in a tokamak comprises the following steps:
step 1, obtaining the polar cross section temperature T of the Tokamak plasma e (R j ,Z k ,t p )；
Wherein T is e Is electron temperature; j is the channel number in the horizontal direction, R j J=1, 2,3 …, M for the horizontal coordinate of the j-th channel; k is the vertical channel number, j=1, 2,3 …, N; z is Z k Is the vertical direction coordinate of the kth channel, t p For time, p=1, 2,3 …, C; m is the number of horizontal channels, N is the number of vertical channels, and C is the number of time points;
step 2, obtaining the characteristic frequency f of the tearing mode corresponding to the maximum power TM The method comprises the steps of carrying out a first treatment on the surface of the Obtaining a relative electron temperature disturbance delta (R) j ，Z k ，t p )；
Step 3, determining the smoothed relative electron temperature disturbance delta 3 (R ii ，Z jj ，t p )
Step 4, smoothing the relative electron temperature disturbance delta 3 (R ii ，Z jj ，t p ) At [ R, Z]Four color images are generated under a coordinate system, the position of a physical surface of a tearing die magnetic field where a tearing die is positioned is determined in each color image, and X points and O points of a tearing die magnetic island are determined;
step 5, determining the space position of the tearing die magnetic island, namely an inner boundary line and an outer boundary line;
step 6, determining the rotation direction of the tearing die;
and 7, determining the polar modulus m of the tearing die.
In the step 1, the measurement temperature T of each measurement channel in the microwave electron cyclotron radiation imaging system is measured e (i,t p ) Polar cross-section temperature T mapped to Tokamak plasma e (R j ,Z k ,t p ) The method comprises the steps of carrying out a first treatment on the surface of the i is the channel number of the microwave electron cyclotron radiation imaging system, i=1, 2, … L, l=m×n, and L is the total channel number.
In the step 2Temperature T of opposite polar section e (R j ,Z k ,t p ) Fourier transform is carried out to obtain the characteristic frequency f of which the frequency corresponding to the maximum power is the tearing mode TM 。
In the step 2, each group T e (R j ,Z k ,t p ) The data is subjected to relative electron temperature disturbance analysis to obtain relative electron temperature disturbance delta (R j ，Z k ，t p ) The calculation formula is as follows:
wherein R is j Is the horizontal coordinate of the jth channel, Z k Is the vertical coordinate of the kth channel, δ (R j ，Z k ，t p ) Is positioned at R j Position, Z k The relative electron temperature at the location and time t is disturbed,<T e (R j ，Z k ，t p )>at t p -1/f TM To t p +1/f TM The relative electron temperature over a period of time perturbs the average.
In the step 3, the method specifically comprises the following steps of
Step 3.1) perturbation of the relative electron temperature obtained in step 2 by delta (R) j ，Z k ，t p ) And performing band-pass filtering analysis to obtain relative electron temperature disturbance delta 1 (R) j ，Z k ，t p )；
Step 3.2) perturbation of the band-pass filtered relative electron temperature by δ1 (R) j ，Z k ，t p ) Performing two-dimensional interpolation to obtain relative electron temperature disturbance delta 2 (R ii ，Z jj ，t p )；
Wherein ii=1, 2,3 … MM-1, jj=1, 2,3 … NN-1, MM and NN are the total number of horizontal and vertical channels, respectively, after difference;
step 3.3) perturbation of the relative electron temperature after the difference by δ2 (R) ii ，Z jj ，t p ) Performing two-dimensional space smoothing to obtain smoothed relative electron temperature disturbance delta3(R ii ，Z jj ，t p ) The calculation formula is as follows:
R ii is the horizontal position of the ii th channel, R ii-1 Is the horizontal position of the ii-1 th channel, R ii+1 Is the horizontal position of the ii+1 channel, Z jj For the vertical position of the jj-th channel, Z jj-1 Is the vertical position of the jj-1 th channel, Z jj+1 Is the vertical position of the jj+1 th channel.
In the step 3.1), the filtering range is [ f ] TM -Δf，f TM +Δf]Δf is f TM /2。
In the step 3.2), the spatial coordinate range of the interpolation is [ R ] min ，R max ]And [ Z ] min ，Z max ]The space coordinate interval after interpolation is 0.5 cm, R min 、Z min 、R max And Z max The minimum and maximum coordinates in the horizontal and vertical directions are measured for electron cyclotron radiation imaging, respectively.
The MM has a value range of 16-60 and NN has a value range of 40-120.
In the step 4, four moments t are selected in one period of the relative electron temperature disturbance 1 ,t 2 ,t 3 And t 4 The smoothed relative electron temperature disturbance δ3 (R) ii ，Z jj ，t p ) At [ R, Z]Four color images are generated under a coordinate system, different colors represent different smoothed relative electronic temperature disturbance values, space points with the values inverted along the radial direction are found on each image, and a curve formed by connecting the space points is a track of the rotation of the center of a magnetic island of the tearing die, namely the position of a magnetic field rational surface where the tearing die is located.
Inversion in the radial direction means that the value changes from positive to negative, or from negative to positive.
In the step 4, if the smoothed relative electron temperature disturbance δ3 (R ii ，Z jj ，t p ) At the position ofA point or points on the rational surface position curve have maxima on the inner side and minima on the outer side, and are X points of the tearing mode magnetic island; if the smoothed relative electron temperature disturbance delta 3 (R ii ，Z jj ，t p ) A minimum value occurs inside and a maximum value occurs outside at a point or points on the rational face location curve, which is the O point of the tearing mode magnetic island.
In step 5, the maximum or minimum value of the smoothed relative electron temperature disturbance at the same radial position in step 4 is determined as the rotation direction of the tearing die in the direction of the polar rotation along time.
The polar modulus m of the tearing die in the step 6 is calculated as the following formula
Wherein θ Maximum value And theta Minimum value of Is the smoothed relative electron temperature disturbance delta 3 (R ii ，Z jj ，t p ) The maximum and minimum values of (a) are at the polar phase angles.
The value range of M is 8-20, the value range of N is 20-30, and the value range of C is 2000-200000.
The invention has the following remarkable effects: the method can identify the information such as the O/X point of the magnetic island, the width of the magnetic island, the polar modulus and the like. Compared with the traditional method based on magnetic measurement inversion or one-dimensional electron cyclotron radiometer, the method is based on direct experimental measurement data without setting any assumption and approximation conditions, and improves accuracy. Meanwhile, the method can be directly compared with the analysis result of the heat transport equation based on the magnetic island, so that the parameter setting of the analysis equation is optimized. The method provided by the invention is based on the microwave electron cyclotron radiation imaging measurement result, has the advantages of simplicity, rapidness, intuitiveness and high efficiency, and can provide an important basis for accurately controlling the tearing die in the Tokamak plasma.
Drawings
FIG. 1 is a schematic view of a Tokamak polar cross section and coordinate definition, wherein a blue dot matrix is a spatial measurement point of a microwave electron cyclotron radiation imaging system;
FIG. 2a is a schematic illustration of a tear-away molded magnetic island projected in polar cross section;
FIG. 2b is a schematic graph showing the evolution of the relative electron temperature disturbance at different locations of the magnetic island at the measurement points inside and outside the magnetic island in FIG. 2 a;
FIG. 3 shows a signal obtained by bandpass filtering the relative electron temperature disturbance of a certain channel in microwave electron cyclotron radiation imaging, t 1 ,t 2 ,t 3 And t 4 Four times in a period;
fig. 4 is a graph showing the result of identifying information such as the position of the magnetic islands, the X/O points, the width of the magnetic islands, and the rotation direction on the microwave electron cyclotron radiation imaging map.
Detailed Description
The invention is further illustrated by the following figures and detailed description.
Tokamak is a device for restraining high-temperature plasma in a ring-shaped container by using strong magnetic field to generate thermonuclear fusion reaction, the cross section of the device is shown in figure 1, R and Z respectively represent horizontal direction and vertical direction, R, theta and phi respectively represent radial radius, polar angle and circumferential angle, O is the center of the Tokamak ring, R 0 Is the center of the plasma, and a is the radius of the plasma.
The points in fig. 1 are shown as measurement positions of a microwave electron cyclotron radiation imaging system, which is a microwave imaging system for measuring the evolution of electron temperature (or temperature fluctuation) with time in the fusion reaction process of a magnetic confinement nuclear fusion tokamak device, and the system is a very important diagnostic tool on tokamak, and belongs to common known equipment.
Fig. 2a is a schematic view of a projection of a tearing mode magnetic island with a polar modulus of 2 in tokamak plasma on a polar cross section, a broken line represents nested magnetic lines of force before the tearing mode magnetic island is formed, a solid line is a magnetic line of force after the magnetic island is formed, the center of the magnetic island is a magnetic island O point, the magnetic lines of force are in reconnection at a magnetic island X point, and the maximum distance from the inner side to the outer side along the radial direction in a region where the reconnection of the magnetic lines of force occurs is the width w of the magnetic island.
According to the principle that the electron temperature of Tokamak plasma monotonically decreases from inside to outside and the electron temperature on the same magnetic surface is equal, the evolution schematic curve of the electron temperature disturbance (the average value of the electron temperature subtracted by the electron temperature) at the inner side and the outer side of the magnetic island at the position (1-8) passing through different magnetic islands can be deduced, as shown in fig. 2 b. It can be seen that the electron temperature phase is opposite inside and outside the magnetic island. At the point of the magnetic island O, the inside electron temperature disturbance is negative, the outside electron temperature disturbance is positive, and at the point of the magnetic island X, the inside electron temperature disturbance is positive, the outside electron temperature disturbance is negative. The electron temperature disturbance is negative maximum at the inner boundary of the magnetic island and positive maximum at the outer boundary of the magnetic island, as viewed in the radial direction through the O-point of the magnetic island.
According to the principle, the electron temperature measured by the microwave electron cyclotron radiation imaging system is subjected to two-dimensional space mapping, band-pass filtering, interpolation, space smoothing and the like to generate an image of relative electron temperature disturbance, and the information of the space position, the two-dimensional structure, the magnetic island width, the polar modulus and the like of the tearing mode magnetic island is clearly identified from the image, and part of the result is shown in figure 4.
Step 1: measurement temperature T of each measurement channel in microwave electron cyclotron radiation imaging system e (i,t p ) Polar cross-section temperature T mapped to Tokamak plasma e (R j ,Z k ,t p )
T e For electron temperature, i is the channel number of the microwave electron cyclotron radiation imaging system, i=1, 2, … L; j is the channel number in the horizontal direction, R j J=1, 2,3 …, M for the horizontal coordinate of the j-th channel; k is the vertical channel number, j=1, 2,3 …, N; z is Z k Is the vertical direction coordinate of the kth channel, t p For time, p=1, 2,3 …, C; where l=m×n, L is the total channel number, M is the horizontal channel number, N is the vertical channel number, and C is the time point number.
Step 2: obtaining the characteristic frequency f of the tearing mode corresponding to the maximum power TM Obtaining relative electron temperature disturbance delta (R j ，Z k ，t p )；
For the polar cross section temperature T obtained in step 1 e (R j ,Z k ,t p ) Fourier transform is carried out to obtain the characteristic frequency f of which the frequency corresponding to the maximum power is the tearing mode TM . For each group T e (R j ,Z k ,t p ) The data is subjected to relative electron temperature disturbance analysis, and the calculation formula is as follows:
wherein R is j Is the horizontal position of the jth channel, Z k For the vertical position of the kth channel, δ (R j ，Z k ，t p ) Is positioned at R j Position, Z k The relative electron temperature at the location and time t is disturbed,<T e (R j ，Z k ，t p )>at t p -1/f TM To t p +1/f TM The relative electron temperature over a period of time perturbs the average.
Step 3: determining the smoothed relative electron temperature disturbance delta 3 (R ii ，Z jj ，t p )
3.1 First, the relative electron temperature disturbance δ (R) obtained in step 2 j ，Z k ，t p ) Band-pass filtering analysis is carried out, and the filtering range is [ f ] TM -Δf，f TM +Δf]Δf is f TM And/2 to obtain δ1 (R j ，Z k ，t p )。
FIG. 3 shows the evolution of a channel relative to the electron temperature disturbance with time, t 1 ,t 2 ,t 3 And t 4 Four times in the same period.
3.2 Next, for δ1 (R) j ，Z k ，t p ) The two-dimensional interpolation is carried out to obtain the relative electron temperature disturbance delta 2 (R ii ，Z jj ，t p ) Wherein ii=1, 2,3 … MM-1, jj=1, 2,3 … NN-1, MM and NN are the total number of horizontal and vertical channels, respectively, after the difference.
The interpolated spatial coordinate range is [ R min ，R max ]And [ Z ] min ，Z max ]The interpolated space coordinate interval is typically taken to be 0.5 cm, R min 、Z min 、R max And Z max The minimum and maximum coordinates in the horizontal and vertical directions are measured for electron cyclotron radiation imaging, respectively, and ii and jj are the channel numbers in the horizontal and vertical directions after interpolation.
3.3 Finally, to delta 2 (R) ii ，Z jj ，t p ) Performing two-dimensional space smoothing to obtain smoothed relative electron temperature disturbance delta 3 (R ii ，Z jj ，t p ) The calculation formula is as follows:
R ii is the horizontal position of the ii th channel, R ii-1 Is the horizontal position of the ii-1 th channel, R ii+1 Is the horizontal position of the ii+1 channel, Z jj For the vertical position of the jj-th channel, Z jj-1 Is the vertical position of the jj-1 th channel, Z jj+1 Is the vertical position of the jj+1 th channel.
Step 4: the smoothed relative electron temperature disturbance delta 3 (R ii ，Z jj ，t p ) At [ R, Z]Four color images are generated under a coordinate system, the position of the rational surface of the tearing die magnetic field where the tearing die is positioned is determined in each color image, and the X point and the O point of the tearing die magnetic island are determined
Four typical moments are selected within one period of the relative electron temperature perturbation, as t=t in fig. 3 1 ,t 2 ,t 3 And t 4 The smoothed relative electron temperature disturbance δ3 (R) ii ，Z jj ，t p ) At [ R, Z]Four color images are generated in the coordinate system, with different colors representing different delta 3 values, as shown in fig. 4.
According to the principle in fig. 2, spatial points where the sign of the relative electron temperature disturbance delta 3 is inverted (positive- & gt negative or negative- & gt positive) along the radial direction are found on each image, and the curve formed by connecting the spatial points is the track of the center rotation of the magnetic island of the tearing die, namely the position of the magnetic field rational surface where the tearing die is located, as shown by the dotted line in fig. 4.
If the relative electron temperature disturbance delta 3 occurs with a maximum on the inside and a minimum on the outside of a certain (or some) point on the curve, that point is the X point of the tearing mode magnetic island, marked by the white letter "X" in the first plot of fig. 4. Conversely, when the relative electron temperature disturbance δ3 has a minimum value inside and a maximum value outside at a certain (or some) point(s) on the curve, the point is the O point of the tearing mode magnetic island, i.e. the center of the magnetic island, and is marked by the white letter "O" in the third subgraph in fig. 4.
Step 5: determining the spatial position of a tear-mode magnetic island
And (3) connecting the space position corresponding to the minimum value of the relative electron temperature disturbance delta 3 inside the tearing die magnetic island O point determined in the step (4) with the adjacent magnetic island X point determined in the step (4) into a curve, wherein the curve is the inner boundary line of the tearing die magnetic island.
And connecting the space position corresponding to the maximum value of the relative electron temperature disturbance delta 3 outside the O point of the tearing die magnetic island with the adjacent X point to form a curve, wherein the curve is the outer boundary line of the tearing die magnetic island.
The distance between the minimum value point and the maximum value point of the relative electron temperature disturbance delta 3 at the inner side and the outer side of the O point of the tearing die magnetic island is the width w of the magnetic island, as shown in a third subgraph t of fig. 4 3 Shown by the double-headed arrow in time.
In this step the center and boundaries of the magnetic islands, i.e. the spatial positions of the magnetic islands, are determined. Compared with the traditional method based on magnetic measurement inversion or one-dimensional electron cyclotron radiometer, the method for determining the space position of the magnetic island is only based on the two-dimensional electron temperature data of experimental measurement, is simpler and more visual to implement, has higher space precision, does not need to set any assumption and approximate conditions, and improves the accuracy.
Step 6: the maximum value (or minimum value) of the relative electron temperature disturbance delta 3 at the same radial position in the step 4 is followed by t 1 →t 4 The direction of rotation in the polar direction is determined as the direction of rotation of the tearing die, as shown in the first sub-graph t of FIG. 4 1 White unidirectional arrow in time.
Step 7: determining the polar phase difference between the maximum value and the minimum value of the smoothed relative electron temperature disturbance as the polar modulus m of the tearing mode, wherein the calculation formula is as follows
Wherein θ Maximum value And theta Minimum value of Is the smoothed relative electron temperature disturbance delta 3 (R ii ，Z jj ，t p ) The maximum and minimum values of (a) are at the polar phase angles.
The polar modulus, i.e. the structure of the magnetic islands, is determined. The magnetic island structure determined by the method can be directly compared with the analysis result of the heat transport equation based on the magnetic island, so that the parameter setting of the analysis equation is optimized.
In addition, based on the space positions and structures of the magnetic islands obtained in the steps 4, 5 and 7, an important basis is provided for accurately controlling the tearing mode magnetic islands in the Tokamak plasma.
Wherein, the value range of M is 8-20, the value range of N is 20-30, the value range of C is 2000-200000, the value range of MM is 16-60, and the value range of NN is 40-120.
Claims (10)
1. The method for identifying the space position and the structure of the magnetic island of the tearing die in the tokamak is characterized by comprising the following steps:
step 1, obtaining the polar cross section temperature T of the Tokamak plasma e (R j ，Z k ，t p )；
Wherein T is e Is electron temperature; j is the channel number in the horizontal direction, R j J=1, 2,3 …, M for the horizontal coordinate of the j-th channel; k is the vertical channel number, j=1, 2,3 …, N; z is Z k Is the vertical direction coordinate of the kth channel, t p For time, p=1, 2,3 …, C; m is the number of horizontal channels, N is the number of vertical channels, and C is the number of time points;
step 2, obtaining the characteristic frequency f of the tearing mode corresponding to the maximum power TM The method comprises the steps of carrying out a first treatment on the surface of the Obtaining the obtainedObtaining relative electron temperature disturbance delta (R) j ，Z k ，t p )；
Step 3, determining the smoothed relative electron temperature disturbance delta 3 (R ii ，Z jj ，t p )
Step 4, smoothing the relative electron temperature disturbance delta 3 (R ii ，Z jj ，t p ) At [ R, Z]Four color images are generated under a coordinate system, the position of a physical surface of a tearing die magnetic field where a tearing die is positioned is determined in each color image, and X points and O points of a tearing die magnetic island are determined;
step 5, determining the space position of the tearing die magnetic island, namely an inner boundary line and an outer boundary line;
connecting the space position corresponding to the minimum value of the relative electron temperature disturbance delta 3 inside the tearing die magnetic island O point determined in the step 4 with the adjacent magnetic island X point determined in the step 4 into a curve, wherein the curve is the inner boundary line of the tearing die magnetic island;
connecting a space position corresponding to the maximum value of the relative electron temperature disturbance delta 3 outside the O point of the tearing die magnetic island with an adjacent X point to form a curve, wherein the curve is an outer boundary line of the tearing die magnetic island;
the distance between the minimum value point and the maximum value point of the relative electron temperature disturbance delta 3 at the inner side and the outer side of the tearing die magnetic island O point is the width w of the magnetic island;
step 6, determining the rotation direction of the tearing die;
determining the maximum value or minimum value of the smoothed relative electron temperature disturbance at the same radial position in the step 4 as the rotation direction of the tearing die in the direction of the polar rotation along with time;
step 7, determining the polar modulus m of the tearing die;
the polar modulus m of the tearing die in the step 7 is calculated as the following formula
Wherein θ Maximum value And theta Minimum value of Is the smoothed relative electron temperature disturbance delta 3 (R ii ，Z jj ，t p ) A kind of electronic deviceThe maximum and minimum are located at the polar phase angle;
in the step 4, four moments t are selected in one period of the relative electron temperature disturbance 1 ，t 2 ，t 3 And t 4 The smoothed relative electron temperature disturbance δ3 (R) ii ，Z jj ，t p ) At [ R, Z]Generating four color images under a coordinate system, wherein different colors represent different smoothed relative electronic temperature disturbance values, and finding out space points of the values, which are inverted along the radial direction, on each image, wherein a curve formed by connecting the space points is a track of the rotation of the center of a magnetic island of a tearing die, namely the position of a magnetic field rational surface where the tearing die is located;
in the step 4, if the smoothed relative electron temperature disturbance δ3 (R ii ，Z jj ，t p ) A maximum value appears inside a certain point or a certain point on the rational surface position curve and a minimum value appears outside the rational surface position curve, and the point or the points are X points of the tearing mode magnetic island; if the smoothed relative electron temperature disturbance delta 3 (R ii ，Z jj ，t p ) A minimum value occurs inside and a maximum value occurs outside at a point or points on the rational face location curve, which is the O point of the tearing mode magnetic island.
2. A method of identifying the spatial location and structure of a tear film magnetic island in a tokamak as claimed in claim 1, wherein: in the step 1, the measurement temperature T of each measurement channel in the microwave electron cyclotron radiation imaging system is measured e (i，t p ) Polar cross-section temperature T mapped to Tokamak plasma e (R j ，Z k ，t p ) The method comprises the steps of carrying out a first treatment on the surface of the i is the channel number of the microwave electron cyclotron radiation imaging system, i=1, 2,..l, l=m×n, L is the total channel number.
3. A method of identifying the spatial location and structure of a tear film magnetic island in a tokamak as claimed in claim 1, wherein: in the step 2, the temperature T of the opposite polar section e (R j ，Z k ，t p ) Making Fourier transformThe frequency corresponding to the maximum power is obtained as the characteristic frequency f of the tearing mode by transformation TM 。
4. A method of identifying the spatial location and structure of a tear film magnetic island in a tokamak as claimed in claim 1, wherein: in the step 2, each group T e (R j ，Z k ，t p ) The data is subjected to relative electron temperature disturbance analysis to obtain relative electron temperature disturbance delta (R j ，Z k ，t p ) The calculation formula is as follows:
wherein R is j Is the horizontal coordinate of the jth channel, Z k Is the vertical coordinate of the kth channel, δ (R j ，Z k ，t p ) Is positioned at R j Position, Z k The relative electron temperature at the location and time t is disturbed,<T e (R j ，Z k ，t p )>at t p -1/f TM To t p +1/f TM The relative electron temperature over a period of time perturbs the average.
5. A method of identifying the spatial location and structure of a tear film magnetic island in a tokamak as claimed in claim 1, wherein: in the step 3, the method specifically comprises the following steps of
Step 3.1) perturbation of the relative electron temperature obtained in step 2 by delta (R) j ，Z k ，t p ) And performing band-pass filtering analysis to obtain relative electron temperature disturbance delta 1 (R) j ，Z k ，t p )；
Step 3.2) perturbation of the band-pass filtered relative electron temperature by δ1 (R) j ，Z k ，t p ) Performing two-dimensional interpolation to obtain relative electron temperature disturbance delta 2 (R ii ，Z jj ，t p )；
Wherein ii=1, 2,3 … MM-1, jj=1, 2,3 … NN-1, MM and NN are the total number of horizontal and vertical channels, respectively, after difference;
step 3.3) perturbation of the relative electron temperature after the difference by δ2 (R) ii ，Z jj ，t p ) Performing two-dimensional space smoothing to obtain smoothed relative electron temperature disturbance delta 3 (R ii ，Z jj ，t p ) The calculation formula is as follows:
R ii is the horizontal position of the ii th channel, R ii-1 Is the horizontal position of the ii-1 th channel, R ii+1 Is the horizontal position of the ii+1 channel, Z jj For the vertical position of the jj-th channel, Z jj-1 Is the vertical position of the jj-1 th channel, Z jj+1 Is the vertical position of the jj+1 th channel.
6. The method for identifying the spatial location and structure of a tear film magnetic island in a tokamak according to claim 5, wherein: in the step 3.1), the filtering range is [ f ] TM -Δf，f TM +Δf]Delta f is f TM /2。
7. The method for identifying the spatial location and structure of a tear film magnetic island in a tokamak according to claim 5, wherein: in the step 3.2), the spatial coordinate range of the interpolation is [ R ] min ，R max ]And [ Z ] min ，Z max ]The space coordinate interval after interpolation is 0.5 cm, R min 、Z min 、R max And Z max The minimum and maximum coordinates in the horizontal and vertical directions are measured for electron cyclotron radiation imaging, respectively.
8. The method for identifying the spatial location and structure of a tear film magnetic island in a tokamak according to claim 5, wherein: the MM has a value range of 16-60 and NN has a value range of 40-120.
9. A method of identifying the spatial location and structure of a tear film magnetic island in a tokamak as claimed in claim 1, wherein: inversion in the radial direction means that the value changes from positive to negative, or from negative to positive.
10. A method of identifying the spatial location and structure of a tear film magnetic island in a tokamak as claimed in claim 1, wherein: the value range of M is 8-20, the value range of N is 20-30, and the value range of C is 2000-200000.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110670960.9A CN113539525B (en) | 2021-06-17 | 2021-06-17 | Method for identifying space position and structure of tearing die magnetic island in tokamak |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110670960.9A CN113539525B (en) | 2021-06-17 | 2021-06-17 | Method for identifying space position and structure of tearing die magnetic island in tokamak |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113539525A CN113539525A (en) | 2021-10-22 |
| CN113539525B true CN113539525B (en) | 2023-08-15 |
Family
ID=78096157
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110670960.9A Active CN113539525B (en) | 2021-06-17 | 2021-06-17 | Method for identifying space position and structure of tearing die magnetic island in tokamak |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113539525B (en) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114896854A (en) * | 2022-05-12 | 2022-08-12 | 大连理工大学 | Simulation method for plasma radiation evolution in tokamak |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4767590A (en) * | 1986-04-25 | 1988-08-30 | The United States Of America As Represented By The United States Department Of Energy | Anomalous - viscosity current drive |
| CN104678330A (en) * | 2013-11-26 | 2015-06-03 | 核工业西南物理研究院 | High precision magnetic field tilt angle measuring system applied to tokamak device |
| CN105277906A (en) * | 2014-06-05 | 2016-01-27 | 核工业西南物理研究院 | High-precision magnetic field inclination angle measuring system applied to magnetic confinement fusion device |
| CN108802650A (en) * | 2018-07-17 | 2018-11-13 | 中国科学院合肥物质科学研究院 | The real-time detection system and method for tearing mode magnetic island position in a kind of plasma |
| CN111308942A (en) * | 2020-03-31 | 2020-06-19 | 中国科学技术大学 | Signal acquisition and processing system and method for active control of field-reversed pinch nuclear fusion device |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US11087891B2 (en) * | 2017-12-21 | 2021-08-10 | General Atomics | Methods and apparatus for mitigating plasma disruption in fusion devices |
-
2021
- 2021-06-17 CN CN202110670960.9A patent/CN113539525B/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4767590A (en) * | 1986-04-25 | 1988-08-30 | The United States Of America As Represented By The United States Department Of Energy | Anomalous - viscosity current drive |
| CN104678330A (en) * | 2013-11-26 | 2015-06-03 | 核工业西南物理研究院 | High precision magnetic field tilt angle measuring system applied to tokamak device |
| CN105277906A (en) * | 2014-06-05 | 2016-01-27 | 核工业西南物理研究院 | High-precision magnetic field inclination angle measuring system applied to magnetic confinement fusion device |
| CN108802650A (en) * | 2018-07-17 | 2018-11-13 | 中国科学院合肥物质科学研究院 | The real-time detection system and method for tearing mode magnetic island position in a kind of plasma |
| CN111308942A (en) * | 2020-03-31 | 2020-06-19 | 中国科学技术大学 | Signal acquisition and processing system and method for active control of field-reversed pinch nuclear fusion device |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113539525A (en) | 2021-10-22 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN113539525B (en) | Method for identifying space position and structure of tearing die magnetic island in tokamak | |
| Klinger et al. | Performance and properties of the first plasmas of Wendelstein 7-X | |
| Fuchert et al. | Blob properties in L-and H-mode from gas-puff imaging in ASDEX upgrade | |
| Jankowski et al. | Approximate analytical solution for induction heating of solid cylinders | |
| Peng et al. | Algorithms for on-line monitoring of micro spheres in an optical tweezers-based assembly cell | |
| Cho et al. | Characterization of SiO2 etching profiles in pulse-modulated capacitively coupled plasmas | |
| Narayanan et al. | Support vector machine and convolutional neural network based approaches for defect detection in fused filament fabrication | |
| Thomas | Driven dust acoustic waves with thermal effects: Comparison of experiment to fluid theory | |
| Jia et al. | Control of three dimensional particle flux to divertor using rotating RMP in the EAST tokamak | |
| CN101074880B (en) | Method for scanning entry safety landing area in moon detector suspension stage | |
| Zuin et al. | Experimental observation of microtearing modes in a toroidal fusion plasma | |
| Krat et al. | MEPHIST-0 tokamak toroidal magnetic field system | |
| Jinbang et al. | Observation of blobs using a gas puff imaging diagnostic on the HL-2A tokamak | |
| Moreno et al. | Real-time investigation of recovery, recrystallization and austenite transformation during annealing of a cold-rolled steel using high energy X-ray diffraction (HEXRD) | |
| Muller et al. | Evolution strategies for film cooling optimization | |
| Maciusowicz et al. | Identification of Grain Oriented SiFe Steels Based on Imaging the Instantaneous Dynamics of Magnetic Barkhausen Noise Using Short-Time Fourier Transform and Deep Convolutional Neural Network | |
| CN204302229U (en) | A kind of interior crossing type giant magnetoresistance array probe of nonferromagnetic light-wall pipe defects detection | |
| Liu et al. | Radially inward particle transport driven by low-frequency instability in cylindrical magnetized plasma | |
| Morimatsu et al. | Development of a scanning nanopipette probe microscope for fine processing using atmospheric pressure plasma jet | |
| Ivo et al. | Study and classification of the Crystallographic Orientation Distribution Function of a non-grain oriented electrical steel using computer vision system | |
| Garland et al. | The structure and poloidal dynamics of blob filaments in TJ-K | |
| Finken et al. | Overview of experiments with the Dynamic Ergodic Divertor on TEXTOR | |
| Shah et al. | Turbulent transport in the core of a trailing half-delta-wing vortex | |
| Zhang et al. | Technical Parameter Analysis of High-Frequency Induction Heating Applied to Steel Plate Bending. | |
| CN104457653A (en) | Displacement detection method for valve element of linear electromagnetic valve |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |