CN113539525A - Method for identifying spatial position and structure of tearing die magnetic island in tokamak - Google Patents

Method for identifying spatial position and structure of tearing die magnetic island in tokamak Download PDF

Info

Publication number
CN113539525A
CN113539525A CN202110670960.9A CN202110670960A CN113539525A CN 113539525 A CN113539525 A CN 113539525A CN 202110670960 A CN202110670960 A CN 202110670960A CN 113539525 A CN113539525 A CN 113539525A
Authority
CN
China
Prior art keywords
tearing
electron temperature
channel
magnetic
delta
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.)
Granted
Application number
CN202110670960.9A
Other languages
Chinese (zh)
Other versions
CN113539525B (en
Inventor
蒋敏
石中兵
陈伟
施培万
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Southwestern Institute of Physics
Original Assignee
Southwestern Institute of Physics
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Southwestern Institute of Physics filed Critical Southwestern Institute of Physics
Priority to CN202110670960.9A priority Critical patent/CN113539525B/en
Publication of CN113539525A publication Critical patent/CN113539525A/en
Application granted granted Critical
Publication of CN113539525B publication Critical patent/CN113539525B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

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 spatial position and the structure of a tearing mode magnetic island in a Tokamak, which is used for obtaining the polar section temperature of a Tokamak plasma, obtaining the characteristic frequency of a tearing mode corresponding to the maximum power, determining the relative electron temperature disturbance after smoothing, determining the physical surface position of a tearing mode magnetic field, judging the X point and the O point of the tearing mode magnetic island, and obtaining the spatial position of the tearing mode magnetic island, the rotation direction of the tearing mode and the polar modulus. Based on direct experimental measurement data without setting any hypothesis and approximate conditions, the accuracy is improved, the method has the advantages of simplicity, rapidness, intuition and high efficiency, and can provide an important basis for accurately controlling the tearing die in the Tokamak plasma.

Description

Method for identifying spatial position and structure of tearing die magnetic island in tokamak
Technical Field
The invention belongs to a magnetic confinement nuclear fusion technology, and particularly relates to a method for identifying the spatial position and structure of a tearing mode magnetic island in a tokamak.
Background
Magnetic confinement nuclear fusion, as a clean energy source with sufficient fuel reserves, is considered to be one of the important ways that can ultimately solve the human energy crisis. Various instabilities are common in magnetically confined nuclear fusion tokamak plasmas, which can severely degrade plasma confinement performance. The tearing mode is one of the most threatening instabilities, and can destroy a well-constrained nested magnetic surface in the tokamak, and a magnetic 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 and radially transported from the inner side of the magnetic island to the outer side of the magnetic island along the direction of the magnetic line. When the magnetic islands are large enough, these charged particles can directly impact the outer wall of the device causing the fusion reaction to terminate and damage the device wall material, causing significant economic loss. Theoretical and experimental research results on tear pattern control show that effective control of the tear pattern can only be achieved if the position of the magnetic islands (particularly the position of the O-points of the magnetic islands) and the spatial structure are accurately determined.
Disclosure of Invention
The invention aims to provide a method for identifying the spatial position and the structure of a tearing die magnetic island in a tokamak, which can simply, quickly, accurately and efficiently obtain the spatial position and the two-dimensional structure of the magnetic island.
The technical scheme of the invention is as follows:
a method for identifying the spatial position and structure of a tearing die magnetic island in a Tokamak comprises the following steps:
step 1, obtaining polar section temperature T of Tokamak plasmae(Rj,Zk,tp)；
Wherein, TeIs the electron temperature; j is the channel number in the horizontal direction, RjIs the horizontal coordinate of the jth channel, j ═ 1,2,3 …, M; k is a vertical channel serial number, j is 1,2,3 …, N; zkIs the vertical coordinate of the k channel, tpTime, 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 mold corresponding to the maximum powerTM(ii) a Obtaining a relative electron temperature disturbance delta (R)j，Zk，tp)；
Step 3, determining the smoothed relative electronic temperature disturbance delta 3 (R)ii，Zjj，tp)
Step 4, smoothing the relative electronic temperature disturbance delta3(Rii，Zjj，tp) In [ R, Z ]]Generating four color images under a coordinate system, determining the rational surface position of a tearing die magnetic field where a tearing die is positioned in each image, and judging an X point and an O point of a tearing die magnetic island;
step 5, determining the spatial 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 measured temperature T of each measuring channel in the microwave electron cyclotron radiation imaging system is measurede(i,tp) Polar cross-sectional temperature T mapped to Tokamak plasmae(Rj,Zk,tp) (ii) a i is the channel number of the microwave electron cyclotron radiation imaging system, i is 1,2, … L, L is M N, and L is the total number of channels.
In the step 2, the temperature T of the antipodal sectione(Rj,Zk,tp) Fourier transform is carried out to obtain the characteristic frequency f of the tearing mode corresponding to the maximum powerTM。
In the step 2, each group Te(Rj,Zk,tp) The data is subjected to relative electron temperature disturbance analysis to obtain relative electron temperature disturbance delta (R)j，Zk，tp) The calculation formula is as follows:
Figure BDA0003119206760000021
wherein R isjIs the horizontal coordinate of the j-th channel, ZkIs the vertical coordinate of the k channel, δ (R)j，Zk，tp) Is located at RjPosition, ZkThe relative electron temperature perturbation at location and time t,<Te(Rj，Zk，tp)>is tp-1/fTMTo tp+1/fTMThe relative electron temperature over the time period perturbs the average.
The step 3 specifically comprises
Step 3.1) disturbance of the relative electron temperature delta (R) obtained in step 2j，Zk，tp) Performing band-pass filtering analysis to obtain band-pass filtered relative electron temperature disturbance delta 1 (R)j，Zk，tp)；
Step 3.2) the filtered relative electron temperature disturbance delta 1 (R) with band-passj，Zk，tp) Two-dimensional interpolation processing is carried out to obtain the relative electron temperature disturbance delta 2 (R) after the difference valueii，Zjj，tp)；
Wherein, ii is 1,2,3 … MM-1, jj is 1,2,3 … NN-1, MM and NN are the total number of horizontal channels and vertical channels after difference respectively;
step 3.3) disturbance delta 2 (R) of the relative electron temperature after the differenceii，Zjj，tp) Two-dimensional space smoothing is carried out to obtain the smoothed relative electron temperature disturbance delta 3 (R)ii，Zjj，tp) The calculation formula is as follows:
Figure BDA0003119206760000031
Riiis the horizontal position of the ii-th passage, Rii-1Is the horizontal position of the ii-1 th channel, Rii+1Is the horizontal position of the ii +1 th channel, ZjjIs the vertical position of the jj channel, Zjj-1Is the vertical position of the jj-1 th channel, Zjj+1Is the vertical position of the jj +1 th channel.
In the step 3.1), the filtering range is [ fTM-Δf，fTM+Δf]Δ f is fTM/2。
In the step 3.2), the spatial coordinate range of the interpolation is [ R ]min，Rmax]And [ Zmin，Zmax]The space coordinate interval after interpolation is 0.5 cm, Rmin、Zmin、RmaxAnd ZmaxThe minimum and maximum coordinates in the horizontal and vertical directions are measured for electron cyclotron radiation imaging, respectively.
The range of MM is 16-60, and the range of NN is 40-120.
In the step 4, four times t are selected in one period of the relative electronic temperature disturbance1,t2,t3And t4Perturbing the smoothed relative electron temperature obtained in step 3 by δ 3 (R)ii，Zjj，tp) In [ R, Z ]]Four color images are generated under a coordinate system, different colors represent different smoothed relative electron temperature disturbance values, space points of which the values are reversed 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 the tearing die magnetic island, namely the position of a magnetic field physical surface where the tearing die is located.
The reversal 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 is disturbed by delta 3 (R)ii，Zjj，tp) The maximum value appears at the inner side of a certain point or points on the position curve of the rational surface, and the minimum value appears at the outer side, and the point or the points are X points of the tearing mode magnetic island; if the smoothed relative electron temperature is disturbed by delta 3 (R)ii，Zjj，tp) And if a minimum value appears at the inner side and a maximum value appears at the outer side of a certain point or certain points on the rational surface position curve, the point or the points are O points of the tearing mode magnetic island.
In the step 5, the maximum value or the minimum value of the smoothed relative electron temperature disturbance at the same radial position in the step 4 is determined as the rotation direction of the tearing die in the direction of polar rotation with time.
The polar modulus m of the tearing die in the step 6 is calculated by the formula
Figure BDA0003119206760000041
Wherein, thetaMaximum valueAnd thetaMinimum valueFor smoothed relative electron temperature disturbance delta 3 (R)ii，Zjj，tp) The maximum and minimum values of (a) and the polar phase angle at which the minimum is located.
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 O/X point of the magnetic island, the width of the magnetic island, the polar modulus and other information. 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 hypothesis and approximate conditions, and improves the accuracy. Meanwhile, the method can be directly compared with the analytic result of the heat transport equation based on the magnetic island, so that the parameter setting of the analytic equation is optimized. The method is based on the microwave electron cyclotron radiation imaging measurement result, has the advantages of simplicity, rapidness, intuition and high efficiency, and can provide an important basis for accurately controlling the tearing mode in the Tokamak plasma.
Drawings
FIG. 1 is a schematic diagram of a polar cross section of a Tokamak and coordinate definition, and a blue dot matrix is a spatial measurement point of a microwave electron cyclotron radiation imaging system;
FIG. 2a is a schematic view of a tear die magnetic island projected in a poloidal cross-section;
FIG. 2b is a graph illustrating the evolution of relative electron temperature perturbations at different locations of the magnetic island for the measurement points inside and outside the magnetic island of FIG. 2 a;
FIG. 3 shows a signal obtained by band-pass filtering the relative electron temperature perturbation of a certain channel in microwave electron cyclotron radiation imaging, t1,t2,t3And t4Four moments in a cycle;
FIG. 4 is a diagram of the result of identifying information such as the position of the magnetic island, the X/O point, the width of the magnetic island, and the rotation direction on the microwave ECR imaging diagram.
Detailed Description
The invention is further illustrated by the accompanying drawings and the detailed description.
Tokamak is a device for confining high-temperature plasma in a ring-shaped container by utilizing a strong magnetic field to generate thermonuclear fusion reaction, the cross section of the device is shown as figure 1, R and Z represent the horizontal direction and the vertical direction respectively, R, theta and phi represent the radial radius and the polar angle respectivelyAnd the circumferential angle, O is the center of the Tokamak ring, R0Is the plasma center and a is the plasma radius.
The point in fig. 1 is indicated as the measurement position of a microwave electron cyclotron radiation imaging system, which is a microwave imaging system, is used for measuring the evolution of the electron temperature (or temperature fluctuation) along with the time in the fusion reaction process of a magnetic confinement nuclear fusion tokamak device, is a very important diagnostic tool on tokamak, and belongs to commonly used known equipment.
Fig. 2a is a schematic view of a projection of a tearing-mode magnetic island with a polar modulus of 2 in a tokamak plasma in a polar section, a dotted line represents an embedded magnetic line before the tearing-mode magnetic island is formed, a solid line represents a magnetic line after the magnetic island is formed, the center of the magnetic island is a magnetic island point O, the magnetic line is reconnected at a magnetic island point X, and the maximum distance from the inner side to the outer side along the radial direction in a region where the magnetic line is reconnected is the width w of the magnetic island.
According to the principle that the electron temperature of the tokamak plasma monotonically decreases from inside to outside and the electron temperature on the same magnetic surface is equal, the schematic curve of the evolution of the electron temperature disturbance (the average value of the electron temperature minus the electron temperature) at positions (1 → 8) passing through different magnetic island positions at the inner side and the outer side of the magnetic island can be deduced as shown in fig. 2 b. It can be seen that the internal and external electron temperatures of the magnetic island are in opposite phases. At point O of the magnetic island, the inside electron temperature perturbation is negative and the outside electron temperature perturbation is positive, at point X of the magnetic island, the inside electron temperature perturbation is positive and the outside electron temperature perturbation is negative. The electron temperature perturbation is negative maximum at the inner boundary of the magnetic island and positive maximum at the outer boundary of the magnetic island as viewed along a radial direction passing 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 other processing to generate an image disturbed by the relative electron temperature, information such as the space position, the two-dimensional structure, the width of the tearing-mode magnetic island, the polar modulus and the like is clearly identified from the image, and partial results are shown in fig. 4.
Step 1: temperature measurement of each measurement channel in microwave electron cyclotron radiation imaging systemDegree Te(i,tp) Polar cross-sectional temperature T mapped to Tokamak plasmae(Rj,Zk,tp)
TeTaking the electron temperature, i is the channel number of the microwave electron cyclotron radiation imaging system, and i is 1,2, … L; j is the channel number in the horizontal direction, RjIs the horizontal coordinate of the jth channel, j ═ 1,2,3 …, M; k is a vertical channel serial number, j is 1,2,3 …, N; zkIs the vertical coordinate of the k channel, tpTime, p ═ 1,2,3 …, C; wherein, L is M N, L is total channel number, M is horizontal channel number, N is vertical channel number, C is time point number.
Step 2: obtaining the characteristic frequency f of the tearing mode corresponding to the maximum powerTMObtaining a relative electron temperature disturbance delta (R)j，Zk，tp)；
For the polar section temperature T obtained in the step 1e(Rj,Zk,tp) Fourier transform is carried out to obtain the characteristic frequency f of the tearing mode corresponding to the maximum powerTM. For each group Te(Rj,Zk,tp) And (3) carrying out relative electron temperature disturbance analysis on the data, wherein the calculation formula is as follows:
Figure BDA0003119206760000071
wherein R isjIs the horizontal position of the jth channel, ZkIs the vertical position of the k-th channel, δ (R)j，Zk，tp) Is located at RjPosition, ZkThe relative electron temperature perturbation at location and time t,<Te(Rj，Zk，tp)>is tp-1/fTMTo tp+1/fTMThe relative electron temperature over the time period perturbs the average.
And step 3: determining a smoothed relative electron temperature disturbance δ 3 (R)ii，Zjj，tp)
3.1) first, the relative electron temperature perturbation obtained in step 2 is measuredMoving delta (R)j，Zk，tp) Performing band-pass filter analysis with a filter range of fTM-Δf，fTM+Δf]Δ f is f TM2, obtaining the delta 1 (R) after band-pass filteringj，Zk，tp)。
FIG. 3 shows the evolution over time of a certain channel after the filtering of the relative electron temperature disturbance, t1,t2,t3And t4Four moments in the same period.
3.2) secondly, on δ 1 (R)j，Zk，tp) Relative electron temperature disturbance delta 2 (R) after two-dimensional interpolation processing is carried out to obtain difference valueii，Zjj，tp) Wherein ii is 1,2,3 … MM-1, jj is 1,2,3 … NN-1, and MM and NN are the total number of horizontal channels and vertical channels after difference, respectively.
The interpolated spatial coordinate range is [ R ]min，Rmax]And [ Zmin，Zmax]The interpolated spatial coordinate interval is typically taken to be 0.5 cm, Rmin、Zmin、RmaxAnd ZmaxMinimum and maximum coordinates in the horizontal and vertical directions are measured for the electron cyclotron radiation imaging, respectively, and ii and jj are channel numbers in the horizontal and vertical directions after interpolation.
3.3) finally, on δ 2 (R)ii，Zjj，tp) Two-dimensional space smoothing is carried out to obtain the smoothed relative electron temperature disturbance delta 3 (R)ii，Zjj，tp) The calculation formula is as follows:
Figure BDA0003119206760000081
Riiis the horizontal position of the ii-th passage, Rii-1Is the horizontal position of the ii-1 th channel, Rii+1Is the horizontal position of the ii +1 th channel, ZjjIs the vertical position of the jj channel, Zjj-1Is the vertical position of the jj-1 th channel, Zjj+1Is the vertical position of the jj +1 th channel.
And 4, step 4: smoothing the relative electronsTemperature disturbance delta 3 (R)ii，Zjj，tp) In [ R, Z ]]Generating four color images under a coordinate system, determining the rational surface position of a tearing die magnetic field where the tearing die is positioned in each image, and judging the X point and the O point of the tearing die magnetic island
Four typical times are selected within one period of the relative electron temperature perturbation, e.g. t-t in fig. 31,t2,t3And t4Perturbing the smoothed relative electron temperature obtained in step 3 by δ 3 (R)ii，Zjj，tp) In [ R, Z ]]Four color images were generated under the coordinate system, with different colors representing different δ 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 δ 3 is reversed (positive → negative or negative → positive) along the radial direction are found on each image, and the connected curves are the tracks of the rotation of the center of the tearing-die magnetic island, namely the position of the magnetic field facet where the tearing die is located, as shown by the dotted lines in fig. 4.
If the relative electron temperature disturbance δ 3 has a maximum value inside and a minimum value outside a certain (or some) point on the curve, this point is the X point of the tear mode magnetic island, which is marked by the white letter "X" in the first sub-graph of fig. 4. Conversely, if the relative electron temperature disturbance δ 3 has a minimum value inside and a maximum value outside a certain (or some) point on the curve, the point is the O point of the tear-mode magnetic island, i.e. the center of the magnetic island, and is marked by the white letter "O" in the third sub-graph of fig. 4.
And 5: determining spatial location of tear die magnetic islands
And (4) connecting the space position corresponding to the minimum value of the relative electron temperature disturbance delta 3 at the inner side of the O point of the tearing-die magnetic island determined in the step (4) and the X point of the adjacent magnetic island 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 and the adjacent X point into a curve, wherein the curve is the outer boundary line of the tearing die magnetic island.
Minimum value of relative electron temperature disturbance delta 3 at inner side and outer side of O point of tearing mode magnetic islandThe distance between the point and the maximum value point is the width w of the magnetic island, as shown in the third sub-graph t of FIG. 43Shown by the double arrow in the moment.
In this step, the center and boundaries of the magnetic island are determined, i.e., the spatial location of the magnetic island is determined. Compared with the traditional method based on magnetic measurement inversion or one-dimensional electron cyclotron radiometer, the method for determining the spatial position of the magnetic island is only based on the two-dimensional electron temperature data measured by experiments, is simpler and more intuitive to implement, has higher spatial precision, does not need to set any hypothesis and approximate conditions, and improves the accuracy.
Step 6: taking the maximum value (or minimum value) of the relative electron temperature disturbance delta 3 at the same radial position in the step 4 as t1→t4The direction of rotation in the polar direction is determined as the direction of rotation of the tearing die, as shown in the first drawing t of FIG. 41White one-way arrows in time.
And 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 mold, and calculating according to the formula
Figure BDA0003119206760000091
Wherein, thetaMaximum valueAnd thetaMinimum valueFor smoothed relative electron temperature disturbance delta 3 (R)ii，Zjj，tp) The maximum and minimum values of (a) and the polar phase angle at which the minimum is located.
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 analytic result of the heat transport equation based on the magnetic island, so that the parameter setting of the analytic equation is optimized.
In addition, the space position and the structure of the magnetic island obtained in the steps 4, 5 and 7 provide important basis for accurately controlling the tearing mode magnetic island 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 (14)

1. A method for identifying the spatial position and structure of a tearing die magnetic island in a Tokamak is characterized by comprising the following steps:
step 1, obtaining polar section temperature T of Tokamak plasmae(Rj,Zk,tp)；
Wherein, TeIs the electron temperature; j is the channel number in the horizontal direction, RjIs the horizontal coordinate of the jth channel, j ═ 1,2,3 …, M; k is a vertical channel serial number, j is 1,2,3 …, N; zkIs the vertical coordinate of the k channel, tpTime, 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 mold corresponding to the maximum powerTM(ii) a Obtaining a relative electron temperature disturbance delta (R)j，Zk，tp)；
Step 3, determining the smoothed relative electronic temperature disturbance delta 3 (R)ii，Zjj，tp)
Step 4, smoothing the relative electronic temperature disturbance delta 3 (R)ii，Zjj，tp) In [ R, Z ]]Generating four color images under a coordinate system, determining the rational surface position of a tearing die magnetic field where a tearing die is positioned in each image, and judging an X point and an O point of a tearing die magnetic island;
step 5, determining the spatial 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.
2. The method of claim 1, wherein the method further comprises the steps of: in the step 1, the measured temperature T of each measuring channel in the microwave electron cyclotron radiation imaging system is measurede(i,tp) Polar cross-sectional temperature T mapped to Tokamak plasmae(Rj,Zk,tp) (ii) a i is the channel number of the microwave electron cyclotron radiation imaging system, i is 1,2, … L, L is M N, and L is the total number of channels.
3. The method of claim 1, wherein the method further comprises the steps of: in the step 2, the temperature T of the antipodal sectione(Rj,Zk,tp) Fourier transform is carried out to obtain the characteristic frequency f of the tearing mode corresponding to the maximum powerTM。
4. The method of claim 1, wherein the method further comprises the steps of: in the step 2, each group Te(Rj,Zk,tp) The data is subjected to relative electron temperature disturbance analysis to obtain relative electron temperature disturbance delta (R)j，Zk，tp) The calculation formula is as follows:
Figure FDA0003119206750000021
wherein R isjIs the horizontal coordinate of the j-th channel, ZkIs the vertical coordinate of the k channel, δ (R)j，Zk，tp) Is located at RjPosition, ZkThe relative electron temperature perturbation at location and time t,<Te(Rj，Zk，tp)>is tp-1/fTMTo tp+1/fTMThe relative electron temperature over the time period perturbs the average.
5. The method of claim 1, wherein the method further comprises the steps of: the step 3 specifically comprises
Step 3.1) disturbance of the relative electron temperature delta (R) obtained in step 2j，Zk，tp) Performing band-pass filtering analysis to obtain band-pass filtered relative electronsTemperature disturbance delta 1 (R)j，Zk，tp)；
Step 3.2) the filtered relative electron temperature disturbance delta 1 (R) with band-passj，Zk，tp) Two-dimensional interpolation processing is carried out to obtain the relative electron temperature disturbance delta 2 (R) after the difference valueii，Zjj，tp)；
Wherein, ii is 1,2,3 … MM-1, jj is 1,2,3 … NN-1, MM and NN are the total number of horizontal channels and vertical channels after difference respectively;
step 3.3) disturbance delta 2 (R) of the relative electron temperature after the differenceii，Zjj，tp) Two-dimensional space smoothing is carried out to obtain the smoothed relative electron temperature disturbance delta 3 (R)ii，Zjj，tp) The calculation formula is as follows:
Figure FDA0003119206750000031
Riiis the horizontal position of the ii-th passage, Rii-1Is the horizontal position of the ii-1 th channel, Rii+1Is the horizontal position of the ii +1 th channel, ZjjIs the vertical position of the jj channel, Zjj-1Is the vertical position of the jj-1 th channel, Zjj+1Is the vertical position of the jj +1 th channel.
6. The method of claim 5, wherein the method further comprises the step of identifying the spatial location and structure of the magnetic islands of tearing mode in the tokamak: in the step 3.1), the filtering range is [ fTM-Δf，fTM+Δf]Δ f is fTM/2。
7. The method of claim 5, wherein the method further comprises the step of identifying the spatial location and structure of the magnetic islands of tearing mode in the tokamak: in the step 3.2), the spatial coordinate range of the interpolation is [ R ]min，Rmax]And [ Zmin，Zmax]The space coordinate interval after interpolation is 0.5 cm, Rmin、Zmin、RmaxAnd ZmaxThe minimum and maximum coordinates in the horizontal and vertical directions are measured for electron cyclotron radiation imaging, respectively.
8. The method of claim 5, wherein the method further comprises the step of identifying the spatial location and structure of the magnetic islands of tearing mode in the tokamak: the range of MM is 16-60, and the range of NN is 40-120.
9. The method of claim 1, wherein the method further comprises the steps of: in the step 4, four times t are selected in one period of the relative electronic temperature disturbance1,t2,t3And t4Perturbing the smoothed relative electron temperature obtained in step 3 by δ 3 (R)ii，Zjj，tp) In [ R, Z ]]Four color images are generated under a coordinate system, different colors represent different smoothed relative electron temperature disturbance values, space points of which the values are reversed 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 the tearing die magnetic island, namely the position of a magnetic field physical surface where the tearing die is located.
10. The method of claim 9, wherein the method further comprises the step of identifying the spatial location and structure of the magnetic islands of tearing mode in the tokamak: the reversal in the radial direction means that the value changes from positive to negative or from negative to positive.
11. The method of claim 1, wherein the method further comprises the steps of: in the step 4, if the smoothed relative electron temperature is disturbed by delta 3 (R)ii，Zjj，tp) The maximum value appears at the inner side of a certain point or points on the position curve of the rational surface, and the minimum value appears at the outer side, and the point or the points are X points of the tearing mode magnetic island; if the smoothed relative electron temperature is disturbed by delta 3 (R)ii，Zjj，tp) The minimum value appears at the inner side of a certain point or certain points on the position curve of the rational surface and appears at the outer sideAt the maximum, this point or these points are the O points of the tearing mode magnetic islands.
12. The method of claim 1, wherein the method further comprises the steps of: in the step 5, the maximum value or the minimum value of the smoothed relative electron temperature disturbance at the same radial position in the step 4 is determined as the rotation direction of the tearing die in the direction of polar rotation with time.
13. The method of claim 1, wherein the method further comprises the steps of: the polar modulus m of the tearing die in the step 6 is calculated by the formula
Figure FDA0003119206750000041
Wherein, thetaMaximum valueAnd thetaMinimum valueFor smoothed relative electron temperature disturbance delta 3 (R)ii，Zjj，tp) The maximum and minimum values of (a) and the polar phase angle at which the minimum is located.
14. The method of claim 1, wherein the method further comprises the steps of: 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.
CN202110670960.9A 2021-06-17 2021-06-17 Method for identifying space position and structure of tearing die magnetic island in tokamak Active CN113539525B (en)

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 true CN113539525A (en) 2021-10-22
CN113539525B 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)

Cited By (1)

* Cited by examiner, † Cited by third party
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 (6)

* Cited by examiner, † Cited by third party
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
US20190198182A1 (en) * 2017-12-21 2019-06-27 General Atomics Methods and apparatus for mitigating plasma disruption in fusion devices
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

Patent Citations (6)

* Cited by examiner, † Cited by third party
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
US20190198182A1 (en) * 2017-12-21 2019-06-27 General Atomics Methods and apparatus for mitigating plasma disruption in fusion devices
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

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114896854A (en) * 2022-05-12 2022-08-12 大连理工大学 Simulation method for plasma radiation evolution in tokamak

Also Published As

Publication number Publication date
CN113539525B (en) 2023-08-15

Similar Documents

Publication Publication Date Title
Jakubowski et al. Change of the Magnetic-Field Topology by an Ergodic Divertor and the Effect<? format?> on the Plasma Structure and Transport
Klinger et al. Performance and properties of the first plasmas of Wendelstein 7-X
Xu et al. Frequency-Resolved Nonlinear Turbulent Energy Transfer into Zonal Flows<? format?> in Strongly Heated L-Mode Plasmas in the HL-2A Tokamak
Chandra et al. Dynamics and symmetries of flow reversals in turbulent convection
CN113539525A (en) Method for identifying spatial position and structure of tearing die magnetic island in tokamak
Narayanan et al. Support vector machine and convolutional neural network based approaches for defect detection in fused filament fabrication
Aga et al. Experimental flow structure investigation of compound angled film cooling
Sechrest et al. Flow and shear behavior in the edge and scrape-off layer of L-mode plasmas in National Spherical Torus Experiment
Spolaore et al. H-mode achievement and edge features in RFX-mod tokamak operation
Zuin et al. Experimental observation of microtearing modes in a toroidal fusion plasma
CN105823797A (en) Inductive thermal image nondestructive testing device based on commensal type magnet yoke coil
Lee et al. Topology optimization for design of segmented permanent magnet arrays with ferromagnetic materials
Dong et al. A study on plasma edge boundary in ergodic layer of LHD based on radial profile measurement of impurity line emissions
Xia et al. The influence of magnetic circuit and operating parameters on the plasma property of 2 cm ECRIT ion source
Shao et al. Velocimetry of edge turbulence during the dithering L–H transition with dynamic programming based time-delay estimation technique in the EAST superconducting tokamak
Chen et al. Finite element analysis of axial flux permanent magnetic hysteresis dampers based on vector Jiles-Atherton model
Krat et al. MEPHIST-0 tokamak toroidal magnetic field system
Takahashi et al. Scrape-Off-Layer Current Model for Filament Structure Observed<? format?> during Edge-Localized Modes in the DIII-D Tokamak
Feng Up-scaling the island divertor along the W7-stellarator line
Muller et al. Evolution strategies for film cooling optimization
CN204302229U (en) A kind of interior crossing type giant magnetoresistance array probe of nonferromagnetic light-wall pipe defects detection
CN103699752A (en) Coupling method for processing moving boundary problems in electromagnetic field based on edge element method
Finken et al. Overview of experiments with the Dynamic Ergodic Divertor on TEXTOR
Smets et al. On the relationship between quadrupolar magnetic field and collisionless reconnection
Iida et al. Analysis of aerodynamic sound source with measurement of static-pressure fluctuation

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