The present invention relates to numerical simulation of a tokamak device in the field of magnetically controlled fusion, and particularly relates to a simulation method for three-dimensional full-space plasma response in EAST tokamak experimental discharge process.
Controlled fusion is an ultimate permanent solution to solve the energy problem. The most promising solution of multiple controlled fusion solutions is the magnetic confinement method. As the most promising device in magnetically controlled fusion, tokamak is favored by researchers, EAST tokamak in Hefei, Anhui province is a full superconducting tokamak device operated in China currently. All parameters of the EAST device are at the forefront of the world. In the tokamak discharge experimental process, to detect a plasma state, the plasma response needs to be excited by an external perturbation current field, and the plasma state is judged by analyzing the plasma response.
At present, plasma magnetic response signals are generally measured by a low field side magnetic probe array on the EAST tokamak device. However, only the outermost magnetic perturbation signal can be obtained by an experimental measurement means, and magnetic signals of plasma core cannot be obtained by the experimental means currently. To obtain the three-dimensional full-space magnetic response distribution of the plasma core and obtain the plasma state information more intuitively, a simulation method used with experiments urgently needed. The method proposed in the present invention can just satisfy this need. In the method, a current signal applied in an external coil on an experiment is directly used as an input parameter for numerical simulation. Firstly, linear equations are solved in the frequency domain space, and then the information in the frequency domain space is converted into a time domain space. The method has high computational efficiency, small numerical error, good convergence and strong numerical stability, and is an accurate and efficient numerical simulation method.
The invention purpose of the present invention is to realize the numerical simulation of three-dimensional full-space plasma magnetic response signals in EAST tokamak discharge experiments, to make up for the deficiency that a magnetic probe array can only measure the magnetic response signals at limited positions.
The technical solution of the present invention is as follows:
A simulation method for three-dimensional full-space plasma response in EAST Tokamak is provided. A disturbance current signal applied in an external coil of an EAST tokamak device is directly used as input information to solve single frequency three-dimensional full-space plasma magnetic response signals, and finally, the magnetic signals in a frequency domain space are converted into a time domain space and superimposed together to ultimately obtain the distribution of the magnetic response signals in a three-dimensional full space at any time. The method has high calculation efficiency and accurate numerical results, and specifically comprises the following steps:
Step 1: according to the cross-section shape of EAST tokamak, meshing a core high-temperature plasma area in an experiment, and storing physical quantity involved in a computational process through nodes obtained by meshing.
Step 2: designing a waveform of an external perturbation current field δJ according to experimental needs, wherein δJ is generally formed by superimposing a group of perturbation current with different frequencies and amplitudes, i.e., δJ=δJ1+δJ2+. . . +δJn.
Step 3: measuring initial equilibrium information (intensity of pressure P0, magnetic field B0 and current J0) by using a diagnostic device on an EAST device and storing in the mesh nodes.
Step 4: storing an external perturbation current signal δJ1 applied in step 2 into a mesh node in a position corresponding to the external coil.
Step 5: solving linearized magnetohydrodynamic equations comprising external driving terms and obtaining a magnetic response signal δB1.
The equations solved here are:
iωρν1=−∇P1+(j1+δJ1)×B0+J0×δB1
iωδB1=∇×(ν1×B0)
iωP1=−ν1·∇P0−ΓP0∇·ν1
j1+δJ1=∇×δB1
where i is an imaginary number symbol; Γ is an adiabatic coefficient; ω is the frequency of an external perturbation current field δJ1; ρ is plasma density; P1 is intensity of pressure of perturbation; j1 is perturbation current density; ν1 is perturbation velocity; and δB1 is the plasma magnetic response signals to be solved.
Step 6: replacing the external perturbation current signal δJ1 in step 4 with δJ2 . . . δJn and repeating step 4 and step 5 for several times until all single frequency magnetic response signals δB1 . . . δBn are obtained.
Step 7: converting the magnetic response signals δB1 . . . δBn, in the frequency domain space into the time domain space, and superimposing the signals together to obtain the three-dimensional full-space plasma magnetic response signal δBt=δB1t+δB2t+. . . +δBnt at any time.
The beneficial effects of the present invention: the present invention realizes the simulation of the three-dimensional full-space plasma magnetic response signals in discharge experiments of the EAST tokamak device, and makes up for the deficiency that the outermost magnetic response signal can only be obtained by an experimental measurement means. At the same time, a single frequency linear model is used for numerical simulation. Finally, the response signals are converted into the time domain space, so that the numerical processing accuracy is higher, and the magnetic response signals at any time can be obtained. The method has high calculation efficiency and good convergence, and is an efficient and stable numerical simulation method.
Specific embodiments of the present invention are further described below in combination with the drawings and the technical solution.
An EAST tokamak device has a cross-sectional shape as shown in
Specific implementation steps are as follows:
Step 1: according to the cross-section shape of EAST tokamak, meshing a core high-temperature plasma area in an experiment, and storing physical quantity involved in a computational process through nodes obtained by meshing.
Step 2: designing a waveform of an external perturbation current field δJ according to experimental needs, wherein δJ is generally formed by superimposing a group of perturbation current with different frequencies and amplitudes, i.e., δJ=δJ1+δJ2+. . . +δJn.
Step 3: measuring initial equilibrium information (intensity of pressure P0, magnetic field B0 and current J0) by using a diagnostic device on an EAST device and storing in the mesh nodes.
Step 4: storing an external perturbation current signal δJ1 applied in step 2 into a mesh node in a position corresponding to the external coil.
Step 5: solving linearized magnetohydrodynamic equations comprising external driving terms and obtaining a magnetic response signal δB1.
The equations solved here are:
iωρν1=−∇P1+(j1+δJ1)×B0+J0×δB1
iωδB1=∇×(ν1×B0)
iωP1=−ν1·∇P0−ΓP0∇·ν1
j1+δJ1=∇×δB1
where i is an imaginary number symbol; Γ is an adiabatic coefficient; ω is the frequency of an external perturbation current field δJ1; ρ is plasma density; P1 is intensity of pressure of perturbation; j1 is perturbation current density; ν1 is perturbation velocity; and δB1 is the plasma magnetic response signals to be solved.
Step 6: replacing the external perturbation current signal δJ1 in step 4 with δJ2 . . . δJn and repeating step 4 and step 5 for several times until all single frequency magnetic response signals δB1 . . . δBn are obtained.
Step 7: converting the magnetic response signals δB1 . . . δBn in the frequency domain space into the time domain space, and superimposing the signals together to obtain the three-dimensional full-space plasma magnetic response signal δBt=δB1t+δB2t+. . . +δBnt at any time.
The above contents are further detailed descriptions of the present invention in combination with preferred technical solutions. The specific implementation of the present invention shall not be considered to be only limited to these descriptions. For those ordinary skilled in the art to which the present invention belongs, several simple deductions and replacements may be made without departing from the conception of the present invention, all of which shall be considered to belong to the protection scope of the present invention.
Number | Date | Country | Kind |
---|---|---|---|
202110032259.4 | Jan 2021 | CN | national |
Filing Document | Filing Date | Country | Kind |
---|---|---|---|
PCT/CN2022/070014 | 1/4/2022 | WO |