Simulation method for three-dimensional full-space plasma response in east tokamak

A simulation method for three-dimensional full-space plasma response in EAST Tokamak, wherein 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 simulation method 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 mesh nodes obtained by meshing;
• step 2:
• designing a waveform of an external perturbation current field oJ according to experimental needs, wherein oJ is formed by superimposing a group of perturbation current with different frequencies and amplitudes, i.e., oJ=oJ1+oJ2+ . . . +oJn;
• 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 oJ1 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 oB1;
• the solved equations are:
• wpv1=−∇P1+(j1+oJ1)×B0+J0×oB1
• woB1=∇×(v1×B0)
• wP1=−v1·∇P0−I P0∇·v1
• where is an imaginary number symbol;
• I′ is an adiabatic coefficient;
• w is the frequency of an external perturbation current field oJ1;
• p is plasma density;
• P1 is intensity of pressure of perturbation;
• j1 is perturbation current density;
• v1 is perturbation velocity;
• and oB1 is the plasma magnetic response signals to be solved;
• step 6:
• replacing the external perturbation current signal oJ1 in step 4 with oJ2 . . . oJn and repeating step 4 and step 5 for several times until all single frequency magnetic response signals oB1 . . . oBn are obtained;
• and
• step 7:
• converting the magnetic response signals oB1 . . . oBn 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 oBt=oB1t+oB2t+ . . . +oBnt at any time.