| CPC G06F 30/25 (2020.01) [G06F 17/13 (2013.01); G06F 30/28 (2020.01); G06F 2113/08 (2020.01)] | 1 Claim |
|
1. A simulation method for electron temperature evolution caused by an EAST tokamak radiofrequency wave, which realizes self-consistently nonlinear evolution of an electron temperature profile after injection of a radiofrequency wave in a real three-dimensional magnetic field configuration of an EAST tokamak, to obtain the electron temperature profile at any moment and provide space distribution situations of the electron temperature; comprising the following steps:
step 1: according to geometrical configuration of an EAST tokamak device, meshing a core high-temperature plasma area in a discharge experiment, and storing physical quantity generated in the process of numerical simulation through nodes obtained by meshing;
step 2: according to a radiofrequency system adopted by the EAST tokamak device, describing a heating source Sec of a radiofrequency wave by adopting a gaussian function, specifically:
 where, R is a horizontal coordinate; z is a vertical coordinate; ξ is circumferential coordinate; S0 is an amplitude of heating source calculated according to transmitting power of the radiofrequency system; R0, Z0 and ξ0 are respectively positions in a horizontal direction, a vertical direction and a circumferential direction of a radiofrequency wave action area, and ΔRd, ΔZd and Δξd are respectively action widths in the horizontal direction, the vertical direction and the circumferential direction;
initializing parameters related to the heating source according to physical size of a radiofrequency system antenna and the transmitting power of the radiofrequency system, comprising a space distribution function, a strength amplitude and a heating area of the heating source; initializing transport coefficients in an electron heat transport equation and spitzer resistivity coefficients in a spitzer resistivity equation in a tokamak configuration according to experimental discharge parameters; and at the same time, obtaining the initial disturbed electron temperature or δTe(0);
Step 3: measuring the initial magnetic field configuration, i.e. distribution information on a magnetic flux function in three-dimensional space by adopting a magnetic probe and a magnetic flux ring on the EAST tokamak device, and converting the measured magnetic flux function into a grid divided in step 1 through a cubic spline interpolation method, thereby obtaining the initial magnetic flux function ψ(0) and storing into a mesh node;
step 4: calculating the evolution with time of the disturbed electron temperature along the current magnetic field configuration according to the electron heat transport equation in the tokamak configuration, and obtaining the disturbed electron temperature δTe(1) at the next moment;
wherein the solved heat transport equation is:
 where, Te is electron temperature and Te=Te0+δTe, Te0 and δTe are respectively equilibrium electron temperature and disturbed electron temperature, t is time, ν is magnetohydrodynamic velocity, κ// and κ⊥ are respectively transport coefficients of the electrons in parallel and vertical magnetic field line directions, the subscripts of // and ⊥ indicates parallel direction and vertical direction respectively, and sec indicates the heating source;
step 5: according to the current electron temperature obtained by calculating the electron heat transport equation in step 4, calculating current plasma resistivity η(1) according to a spitzer resistivity equation;
wherein the solved spitzer resistivity equation is:
η=κSTe3/2
where, η is plasma resistivity, and κS is a spitzer resistivity coefficient;
step 6: substituting the plasma resistivity obtained in step 5 into a magnetic flux evolution equation, calculating the magnetic flux at the next moment and obtaining the magnetic flux ψ(1) after nonlinear evolution;
wherein the solved magnetic flux evolution equation is:
 where, ψ is magnetic flux, ϕ is potential, η is plasma resistivity, j is total current of the plasmas, and jbs is bootstrap current;
step 7: outputting the three-dimensional space distribution information on the disturbed electron temperature δTe into a computer to store; and
step 8: calculating the current magnetic field configuration according to the magnetic flux obtained in step 6, and repeating steps 4-8, thereby obtaining the three-dimensional space distribution δTe(n) of the disturbed electron temperature at any moment.
|