Method, device, terminal and storage medium for identifying electron cyclotron wave trajectory

Abstract

Embodiments of the present invention propose a method, device, terminal, and storage medium for identifying electron cyclotron wave trajectories. The method comprises: obtaining normalized equilibrium poloidal flux under multiple different preset plasma configurations; determining the total magnetic field based on the normalized equilibrium poloidal flux; and integrating the wave trajectory equation in cylindrical coordinates over the entire interval using the fourth-order Runge-Kutta method with an indefinite step size based on the total magnetic field to sequentially determine information about each point in the electron cyclotron wave propagation trajectory. This solution determines the total magnetic field based on the magnetic flux, and further determines the electron cyclotron wave propagation trajectory based on the total magnetic field, providing a novel approach for simulating heating and driving of ECWs in plasmas.

Description

Technical Field

[0001] The present invention relates to the technical field of high energy physics, and in particular to a method, device, terminal and storage medium for identifying electron cyclotron wave trajectories. Background Art

[0002] Electron cyclotron waves (ECWs) are a type of radio frequency wave. ECW injection is crucial for initiating discharges and maintaining steady-state operation in magnetic confinement fusion plasma devices (such as the HL-2A / EAST tokamak). Research on electron cyclotron wave physics has led to research on wave antennas and driving currents. The calculation of electron cyclotron wave current drive relies on the identification of wave trajectories. Methods for identifying wave trajectories in plasmas have been developed. These methods can be divided into three categories based on their underlying principles:

[0003] One method is to identify wave paths based on ray tracing. This method uses geometric optics approximation and ordinary differential methods to solve the wave path equations to obtain wave path information in thermal plasma. There are a large number of literature and patents on this method, such as Alberti S et al., Full absorption of third harmonic ECH in TCV tokamak plasmas in the presence of second harmonic ECCD, Nuclear fusion, 2002, 42(1):42. Lin Liu YR et al., Electron cyclotron current drive efficiency in general tokamak geometry [J]. Phys. of Plasmas, 2003, 10(10), 4064-4071.

[0004] Another method is to identify wave paths based on the Gaussian beam method. Specifically, when interference or diffraction is important, the Gaussian beam method can address the shortcomings of ray tracing simulation programs. For example, when the beam is emitted with strong focus, the geometric optics of the ray tracing code may cause all rays to intersect at a point in space in an unphysical way, while the Gaussian beam code will find a non-zero beam width under the non-zero diffraction limit. There are also a large number of literatures on this type of method, such as Poli E et al., TORBEAM 2.0, a paraxial beam tracing code for electron-cyclotron beams infusion plasmas for extended physics applications. 2018 Computer Physics Communications 225, 36-46. GVPereverzev et al., Beam tracing in inhomogeneous anisotropic plasmas, 1998, Physics of Plasmas 5, 3529.

[0005] The third method is to identify the wave path based on the quasi-optical method. The quasi-optical approximation method tracks the trajectory to simulate the Gaussian beam, but unlike the ray tracing method, it retains the interaction between the light rays to take into account the interference and diffraction effects. This method reduces the limitations of the above-mentioned Gaussian beam method. Even when the absorption and refraction of the entire beam are uneven, Gaussian distribution fitting is required. There are also relevant literature on this research, such as D. Farina et al., A Quasi-Optical Beam-Tracing Code for Electron Cyclotron Absorption and Current Drive: GRAY, Fusion Sci. Technol. 52 (2007) 154.

[0006] Although there are a large number of literatures and patents on the propagation trajectory of electron cyclotron waves in plasma, the existing methods do not consider the propagation trajectory and absorption of electron cyclotron waves under various configurations, resulting in limitations in practical application processes. Summary of the Invention

[0007] In view of this, the present invention proposes a method, device, terminal and storage medium for identifying electron cyclotron wave trajectories to solve the problems in the prior art.

[0008] Specifically, the present invention proposes the following specific embodiments:

[0009] An embodiment of the present invention provides a method for identifying electron cyclotron wave trajectories, comprising:

[0010] Obtain the normalized equilibrium poloidal magnetic flux under various preset plasma configurations;

[0011] determining a total magnetic field based on the normalized equilibrium poloidal flux;

[0012] Based on the total magnetic field, the wave track equation in cylindrical coordinates is fully integrated using the fourth-order Runge-Kutta method with an indefinite step size to sequentially determine information of each point in the propagation trajectory of the electron cyclotron wave.

[0013] In a specific embodiment, the normalized balanced poloidal flux is expressed by the following formula:

[0014]

[0015] Where ψ(r, z) is the normalized equilibrium poloidal flux, r is the radial coordinate, z is the vertical coordinate, ψ0 is the magnetic axis position, R0 is the lateral radius of the tokamak device, R x , E and τ are the position of the intersection of magnetic flux lines, the elongation ratio and the three angles respectively.

[0016] In a specific embodiment, determining the total magnetic field based on the normalized balanced poloidal magnetic flux includes:

[0017] Derivative the normalized balanced polar magnetic flux in the vertical direction to obtain a magnetic field in the radial direction;

[0018] Derivative the normalized balanced polar magnetic flux in the radial direction to obtain a magnetic field in the vertical direction;

[0019] The total magnetic field is determined based on the magnetic field in the radial direction, the magnetic field in the vertical direction, and the toroidal magnetic field.

[0020] In a specific embodiment, the total magnetic field is expressed by the following formula:

[0021]

[0022]

[0023]

[0024]

[0025]

[0026] Where B is the total magnetic field, B r is the magnetic field in the r direction, B zis the magnetic field in the Z direction, B φ is the toroidal magnetic field, q0 is the central safety factor, B0 is the central magnetic field, R0 is the lateral radius of the tokamak device, and r is the radial coordinate.

[0027] In a specific embodiment, the plasma configurations include: tokamak, spherical tokamak and FRC configurations; wherein different plasma configurations correspond to different magnetic field parameters; the magnetic field parameters include: the position of the intersection of magnetic flux lines, elongation ratio, three angles, center safety factor, center magnetic field and lateral radius of the tokamak device.

[0028] In a specific embodiment, the wave trajectory equation includes:

[0029]

[0030]

[0031]

[0032] Where r is the radial coordinate, φ is the annular angle, z is the vertical coordinate, n = kc / ω = (n r ,m=rn φ , n z ), k is the wave vector, c is the speed of light, ω is the angular frequency of the wave; D0 is the dispersion function calculated using the cold plasma approximation;

[0033]

[0034] D0=D r +iD i ;D r is the real part of plasma dispersion; Di is the imaginary part of plasma dispersion;

[0035] D0=AN 4 -BN 2 +C=0;

[0036] A=Ssin 2 θ+Pcos 2 θ;

[0037] B=RLsin 2 θ+PS(1+cos 2 θ);

[0038] C = PRL;

[0039]

[0040]

[0041]

[0042] is the cyclotron frequency of the particle, B is the total magnetic field, ω s is the plasma frequency, q s is the charge of the particle, m s is the mass of the particle, n s0 is the center density of the particle, ∈0 is the conductivity, θ is the angle between the magnetic field and the wave vector; r, φ, z, k r , n φ , k z The information that makes up the point.

[0043] In a specific embodiment, the method further includes: determining the absorption power of the electron cyclotron wave at each of the points based on the imaginary part of the dispersion function in the wave trajectory equation;

[0044] The absorbed power is calculated based on the following formula:

[0045]

[0046]

[0047] Where, P(l) is the absorbed power, k r is the radial wave vector; D r is the real part of plasma dispersion; D i is the imaginary part of plasma dispersion; P(0) is the total absorbed power.

[0048] The embodiment of the present invention further provides a device for identifying electron cyclotron wave trajectories in different plasma configurations, comprising:

[0049] The magnetic flux module is used to obtain the normalized equilibrium poloidal magnetic flux under various preset plasma configurations;

[0050] a magnetic field module for determining a total magnetic field based on the normalized balanced poloidal flux;

[0051] The trajectory point module is used to perform full interval integration of the wave track equation in cylindrical coordinates based on the total magnetic field and using the fourth-order Runge-Kutta method with an indefinite step size, so as to sequentially determine the information of each point in the propagation trajectory of the electron cyclotron wave.

[0052] An embodiment of the present invention further provides a terminal, comprising: a processor and a memory, wherein the memory stores a computer program, and the processor implements the above method when executing the computer program.

[0053] An embodiment of the present invention further provides a storage medium, wherein the storage medium stores a computer program, and the computer program implements the above method when executed.

[0054] Thus, embodiments of the present invention propose a method, apparatus, terminal, and storage medium for identifying electron cyclotron wave trajectories. The method includes: obtaining normalized equilibrium poloidal magnetic flux under a variety of different preset plasma configurations; determining the total magnetic field based on the normalized equilibrium poloidal magnetic flux; and integrating the wave trajectory equation in cylindrical coordinates over the entire interval using the fourth-order Runge-Kutta method with an indefinite step size based on the total magnetic field to sequentially determine information about each point in the propagation trajectory of the electron cyclotron wave. This solution determines the total magnetic field based on the magnetic flux, and further determines the propagation trajectory of the electron cyclotron wave based on the total magnetic field, providing a new approach for simulating the heating and driving of ECWs in plasmas. BRIEF DESCRIPTION OF THE DRAWINGS

[0055] In order to more clearly illustrate the technical solutions of the embodiments of the present invention, the following briefly introduces the drawings required for use in the embodiments. It should be understood that the following drawings only illustrate certain embodiments of the present invention and therefore should not be regarded as limiting the scope. For ordinary technicians in this field, other relevant drawings can be obtained based on these drawings without making any creative efforts.

[0056] Figure 1 A schematic diagram of a process framework of a method for identifying electron cyclotron wave trajectories in different plasma configurations proposed in an embodiment of the present invention;

[0057] Figure 2A and Figure 2B and Figure 2C A schematic diagram comparing the method proposed in an embodiment of the present invention and the calculation of wave trajectories under a tokamak configuration using GENRAY;

[0058] Figure 3A and Figure 3B Schematic diagram of the propagation trajectory of the 60 GHz electron cyclotron wave in FRC plasma obtained based on this scheme;

[0059] Figure 4A 、 Figure 4B 、 Figure 4C 、 Figure 4D Schematic diagram of the propagation trajectory of 28GHz electron cyclotron wave in the spherical tokamak obtained based on this scheme;

[0060] Figure 5 A schematic diagram of the framework structure of a device for identifying electron cyclotron wave trajectories in different plasma configurations proposed in an embodiment of the present invention. DETAILED DESCRIPTION

[0061] Hereinafter, various embodiments of the present disclosure will be described more fully. The present disclosure may have various embodiments, and adjustments and changes may be made therein. However, it should be understood that there is no intention to limit the various embodiments of the present disclosure to the specific embodiments disclosed herein, but rather that the present disclosure should be construed to encompass all adjustments, equivalents, and / or alternatives falling within the spirit and scope of the various embodiments of the present disclosure.

[0062] The terms used in the various embodiments of the present disclosure are only used to describe the purpose of specific embodiments and are not intended to limit the various embodiments of the present disclosure. As used herein, the singular form is intended to also include the plural form, unless the context clearly indicates otherwise. Unless otherwise specified, all terms used herein (including technical terms and scientific terms) have the same meaning as those generally understood by those skilled in the art of the various embodiments of the present disclosure. The terms (such as those defined in generally used dictionaries) will be interpreted as having the same meaning as the contextual meaning in the relevant technical field and will not be interpreted as having an idealized meaning or an overly formal meaning, unless clearly specified in the various embodiments of the present disclosure.

[0063] In practical applications, the calculation of electron cyclotron wave trajectories in plasmas with different configurations is often required. For example, electron cyclotron waves can serve as discharge and current drivers in tokamaks, current drivers in spheroidal maks, and even as effective heating methods in field reverse pinch devices. Therefore, considering the propagation of electron cyclotron waves in various configurations is crucial. This proposal presents a method for identifying electron cyclotron wave trajectories in different plasma configurations. This method not only identifies the trajectory of plasma electron cyclotron waves but also calculates parameters such as the power absorption of the ECW in different configurations.

[0064] Specifically, embodiment 1 of the present invention discloses a method for identifying electron cyclotron wave trajectories in different plasma configurations, such as Figure 1 As shown, the following steps are included:

[0065] Step 101: obtaining normalized equilibrium poloidal magnetic fluxes under a plurality of different preset plasma configurations;

[0066] Specifically, we first construct an analytical Solovev equilibrium, which includes multiple plasma configurations, such as tokamaks, spherical tokamaks, and FRC configurations, in the same model. This is also the solution of the Grad-Shafranov MHD equilibrium equation, from which we obtain the normalized equilibrium poloidal flux. Specifically, the normalized equilibrium poloidal flux is expressed by the following formula:

[0067]

[0068] Where ψ(r,z) is the normalized equilibrium poloidal flux, r is the radial coordinate, z is the vertical coordinate, ψ0 is the magnetic axis position, R0 is the lateral radius of the tokamak device, R x , E and τ are the position of the intersection of magnetic flux lines, the elongation ratio and the three angles respectively.

[0069] Step 102: determining a total magnetic field based on the normalized balanced poloidal magnetic flux;

[0070] In a specific embodiment, the determining of the total magnetic field based on the normalized balanced poloidal flux in step 102 includes:

[0071] Derivative the normalized balanced polar magnetic flux in the vertical direction to obtain a magnetic field in the radial direction;

[0072] Derivative the normalized balanced polar magnetic flux in the radial direction to obtain a magnetic field in the vertical direction;

[0073] The total magnetic field is determined based on the magnetic field in the radial direction, the magnetic field in the vertical direction, and the toroidal magnetic field.

[0074] Specifically, the total magnetic field is expressed by the following formula:

[0075]

[0076]

[0077]

[0078]

[0079]

[0080] Where B is the total magnetic field, B r is the magnetic field in the r direction, B z is the magnetic field in the z direction, B φ is the toroidal magnetic field, q0 is the central safety factor, B0 is the central magnetic field, R0 is the lateral radius of the tokamak device, and r is the radial coordinate.

[0081] Specifically, at the intersection of magnetic flux lines, B z (R x ,Z x )=0, and we have:

[0082]

[0083] The expression of ψ0 is: q0 is the central safety factor and B0 is the central magnetic field.

[0084] The toroidal magnetic field is:

[0085]

[0086] In obtaining B r ,B z , and B φ Then, the total magnetic field can be obtained

[0087] The plasma configurations include: tokamak, spherical tokamak and FRC configurations; wherein different plasma configurations correspond to different magnetic field parameters; the magnetic field parameters include: the position of the intersection point of magnetic flux lines, elongation ratio, three angles, center safety factor, center magnetic field and the lateral radius of the tokamak device.

[0088] Here, taking the FRC configuration as an example, in order to construct the FRC configuration, we can set τ = 0, R x =0 and B φ = 0, which means Hill-vortex equilibrium is generated. In summary, for tokamaks, spherical tokamaks and FRC configurations, by giving a specific R x , E, B0, R0, and τ (typical values of these parameters for the three configurations are shown in Table 1). Substituting these quantities into Equation (1) in the two-dimensional (r, z) plane yields the plasma poloidal flux for each configuration. By derivatizing Equations (2)-(3), we obtain the magnetic field components and, ultimately, the total magnetic field. In subsequent calculations, we only need to select the total magnetic field information for one configuration as the initial condition for the calculation.

[0089] Table 1

[0090] Parameters / Configuration Tokamak Spherical Tokamak FRC <![CDATA[R x ]]> 0.75 0.17 0.0 E 1.8 1.5 2.0 <![CDATA[B0(T)]]> 1.78 0.32 0.5 τ 0.5 0.8 0.0 <![CDATA[R0]]> 1.87 0.64 0.5 <![CDATA[q0]]> 1.50 1.58 /

[0091] Step 103: Based on the total magnetic field, the wave trajectory equation in cylindrical coordinates is integrated over the entire interval using a fourth-order Runge-Kutta method with an indefinite step size to sequentially determine information of each point in the propagation trajectory of the electron cyclotron wave;

[0092] The wave trajectory equation includes:

[0093]

[0094]

[0095]

[0096] Where r is the radial coordinate, φ is the annular angle, z is the vertical coordinate, n = kc / ω = (n r ,m=rn φ , nz ), k is the wave vector, c is the speed of light, ω is the angular frequency of the wave; D0 is the dispersion function calculated using the cold plasma approximation; n r is n in the r direction; n z is n, n in the r direction φ is the circumferential n.

[0097]

[0098] D0=D r +iD i ;D r is the real part of plasma dispersion; D i is the imaginary part of plasma dispersion;

[0099] D0=AN 4 -BN 2 +C=0;

[0100] A=Ssin 2 θ+Pcos 2 θ;

[0101] B=RLsin 2 θ+PS(1+cos 2 θ);

[0102] C = PRL;

[0103]

[0104]

[0105]

[0106] is the cyclotron frequency of the particle, B is the total magnetic field, ω s is the plasma frequency, q s is the charge of the particle, m s is the mass of the particle, n s0 is the center density of the particle, ∈0 is the conductivity, θ is the angle between the magnetic field and the wave vector; r, φ, z, k r , n φ , k z The information that makes up the point.

[0107] Furthermore, the method also includes: determining the absorption power of the electron cyclotron wave at each of the points based on the imaginary part of the dispersion function in the wave trajectory equation.

[0108] The absorbed power is calculated based on the following formula:

[0109]

[0110]

[0111] Where, P(l) is the absorbed power, k r is the radial wave vector; D r is the real part of plasma dispersion; D i is the imaginary part of plasma dispersion; P(0) is the total absorbed power.

[0112] Through the above method, each point in the propagation trajectory of the electron cyclotron wave can be determined, and a complete trajectory can be summarized based on each point. In addition, based on the absorption power of each point obtained by this scheme, the propagation and power absorption of the wave trajectory in the entire plasma area can be obtained.

[0113] To verify the practical effects of the present invention, we used electron cyclotron injection in a tokamak configuration for testing. This configuration is known (and can be obtained using the internationally recognized GENRAY program. GENRAY is limited in that it can only calculate wave trajectories in tokamak configurations). Therefore, we can compare the two results to verify the accuracy and effectiveness of the present invention. First, with initial conditions of r = 2.28, φ = 1.21, and z = -0.0304, the propagation trajectory of a 110 GHz electron cyclotron wave in the EAST tokamak plasma is as follows: Figure 2A and Figure 2B. The continuous lines in the three figures represent the results calculated using GENRAY, while the dashed lines represent the results calculated using this solution. It can be seen that the results calculated using this solution are consistent with those given by GENRAY in terms of both wave trajectory and absorption efficiency, demonstrating the reliability of this solution.

[0114] Then, based on this scheme, the identification of electron cyclotron waves is extended to FRC configurations. Figure 3A Figure 3B shows the propagation trajectory of the 60 GHz electron cyclotron wave in the FRC plasma. It can be seen that the wave trajectory is an approximately straight line propagation, which is due to the small circumferential emission angle of this solution.

[0115] In addition, based on this scheme, this method is used to identify the wave trajectory under the spherical tokamak configuration. Figure 4A 4B, 4C and 4D are the propagation of 28GHz electron cyclotron wave in spherical tokamak. Figure 4A is the wave trace under the small cross section, Figure 4B is the two-dimensional density distribution of plasma, Figure 4C is the wave trace under large cross section, Figure 4DThe temporal variation of the wave number components is shown in Figure 1. Due to the large parallel refractive index, the wave trajectory undergoes a significant deflection and returns. This demonstrates that the proposed method can successfully identify electron cyclotron wave trajectories in tokamaks, spherical tokamaks, and FRC configurations. Compared to traditional methods, the present invention offers advantages such as applicability to any configuration and the ability to determine absorption rates.

[0116] In order to further illustrate this solution, the second embodiment of the present invention further discloses a device for identifying electron cyclotron wave trajectories, such as Figure 5 Shown, including:

[0117] The magnetic flux module 201 is used to obtain the normalized equilibrium poloidal magnetic flux under a variety of different preset plasma configurations;

[0118] A magnetic field module 202 is configured to determine a total magnetic field based on the normalized balanced poloidal flux;

[0119] The trajectory point module 203 is used to perform full interval integration of the wave track equation in cylindrical coordinates based on the total magnetic field and using the fourth-order Runge-Kutta method with an indefinite step size to sequentially determine the information of each point in the propagation trajectory of the electron cyclotron wave.

[0120] In a specific embodiment, the normalized balanced poloidal flux is expressed by the following formula:

[0121]

[0122] Where ψ(r, z) is the normalized equilibrium poloidal flux, r is the radial coordinate, z is the vertical coordinate, ψ0 is the magnetic axis position, R0 is the lateral radius of the tokamak device, R x , E and τ are the position of the intersection of magnetic flux lines, the elongation ratio and the three angles respectively.

[0123] In a specific embodiment, the magnetic field module 202 includes:

[0124] a radial module, configured to derive the normalized balanced polar magnetic flux in the vertical direction to obtain a magnetic field in the radial direction;

[0125] A vertical module, configured to derive the normalized balanced polar magnetic flux in a radial direction to obtain a magnetic field in a vertical direction;

[0126] The determination module is configured to determine the total magnetic field based on the magnetic field in the radial direction, the magnetic field in the vertical direction, and the toroidal magnetic field.

[0127] In a specific embodiment, the total magnetic field is expressed by the following formula:

[0128]

[0129]

[0130]

[0131]

[0132]

[0133] Where B is the total magnetic field, B r is the magnetic field in the r direction, B z is the magnetic field in the Z direction, B φ is the toroidal magnetic field, q0 is the central safety factor, B0 is the central magnetic field, R0 is the lateral radius of the tokamak device, and r is the radial coordinate.

[0134] In a specific embodiment, the plasma configurations include: tokamak, spherical tokamak and FRC configurations; wherein different plasma configurations correspond to different magnetic field parameters; the magnetic field parameters include: the position of the intersection of magnetic flux lines, elongation ratio, three angles, center safety factor, center magnetic field and lateral radius of the tokamak device.

[0135] In a specific embodiment, the wave trajectory equation includes:

[0136]

[0137]

[0138]

[0139] Where r is the radial coordinate, φ is the annular angle, z is the vertical coordinate, n = kc / ω = (n r ,m=rn φ , n z ), k is the wave vector, c is the speed of light, ω is the angular frequency of the wave; D0 is the dispersion function calculated using the cold plasma approximation;

[0140]

[0141] D0=D r +iD i ;D r is the real part of plasma dispersion; D i is the imaginary part of plasma dispersion;

[0142] D0=AN 4 -BN 2 +C=0;

[0143] A=Ssin2 θ+Pcos 2 θ;

[0144] B=RLsin 2 θ+PS(1+cos 2 θ);

[0145] C = PRL;

[0146]

[0147]

[0148]

[0149] is the cyclotron frequency of the particle, B is the total magnetic field, ω s is the plasma frequency, q s is the charge of the particle, m s is the mass of the particle, n s0 is the center density of the particle, ∈0 is the conductivity, θ is the angle between the magnetic field and the wave vector; r, φ, z, k r , n φ , k z The information that makes up the point.

[0150] In a specific embodiment, the method further includes: an absorption power module 204 for determining the absorption power of the electron cyclotron wave at each of the points based on the imaginary part of the dispersion function in the wave trajectory equation;

[0151] The absorbed power is calculated based on the following formula:

[0152]

[0153]

[0154] Where, P(l) is the absorbed power, k r is the radial wave vector; D r is the real part of plasma dispersion; D i is the imaginary part of plasma dispersion; P(0) is the total absorbed power.

[0155] Embodiment 3 of the present invention further discloses a terminal, including: a processor and a memory, wherein the memory stores a computer program, and when the processor executes the computer program, the method described in embodiment 1 is implemented.

[0156] Embodiment 4 of the present invention further discloses a storage medium, in which a computer program is stored. When the computer program is executed, the method described in embodiment 1 is implemented.

[0157] Embodiments of the present invention propose a method, device, terminal, and storage medium for identifying electron cyclotron wave trajectories. The method includes: obtaining normalized equilibrium poloidal magnetic flux under a variety of different preset plasma configurations; determining the total magnetic field based on the normalized equilibrium poloidal magnetic flux; and integrating the wave trajectory equation in cylindrical coordinates over the entire interval using the fourth-order Runge-Kutta method with an indefinite step size based on the total magnetic field to sequentially determine information about each point in the propagation trajectory of the electron cyclotron wave. This solution determines the total magnetic field based on the magnetic flux, and further determines the propagation trajectory of the electron cyclotron wave based on the total magnetic field, providing a new approach for simulating the heating and driving of ECWs in plasmas.

[0158] Those skilled in the art will understand that the accompanying drawings are merely schematic diagrams of a preferred implementation scenario, and the modules or processes in the accompanying drawings are not necessarily required to implement the present invention.

[0159] Those skilled in the art will appreciate that the modules in the devices in the implementation scenario can be distributed in the devices of the implementation scenario according to the implementation scenario description, or can be modified accordingly and located in one or more devices different from the implementation scenario. The modules in the above implementation scenario can be combined into one module or further split into multiple submodules.

[0160] The above serial numbers of the present invention are for description only and do not represent the advantages or disadvantages of the implementation scenarios.

[0161] The above disclosures are only several specific implementation scenarios of the present invention. However, the present invention is not limited thereto. Any changes that can be conceived by those skilled in the art should fall within the scope of protection of the present invention.

Claims

1. A method for identifying electron cyclotron wave trajectories, characterized in that: include: Obtain the normalized equilibrium poloidal magnetic flux under various preset plasma configurations; determining a total magnetic field based on the normalized equilibrium poloidal flux; Based on the total magnetic field, a full-interval integration of the wave trajectory equation in cylindrical coordinates is performed using a fourth-order Runge-Kutta method with an indefinite step size to sequentially determine information of each point in the propagation trajectory of the electron cyclotron wave; The normalized equilibrium poloidal flux is expressed by the following formula: Where ψ(r,z) is the normalized equilibrium poloidal flux, r is the radial coordinate, z is the vertical coordinate, ψ0 is the magnetic axis position, R0 is the lateral radius of the tokamak device, R x , E and τ are the position of the intersection of magnetic flux lines, the elongation ratio and the three angles respectively.

2. The method according to claim 1, wherein Determining the total magnetic field based on the normalized balanced poloidal magnetic flux comprises: Derivative the normalized balanced polar magnetic flux in the vertical direction to obtain a magnetic field in the radial direction; Derivative the normalized balanced polar magnetic flux in the radial direction to obtain a magnetic field in the vertical direction; The total magnetic field is determined based on the magnetic field in the radial direction, the magnetic field in the vertical direction, and the toroidal magnetic field.

3. The method according to claim 2, wherein The total magnetic field is expressed by the following formula: Where B is the total magnetic field, B r is the magnetic field in the r direction, B z is the magnetic field in the z direction, B φ is the toroidal magnetic field, q0 is the central safety factor, B0 is the central magnetic field, R0 is the lateral radius of the tokamak device, and r is the radial coordinate.

4. The method according to claim 1, wherein The plasma configurations include: tokamak, spherical tokamak and FRC configurations; wherein different plasma configurations correspond to different magnetic field parameters; the magnetic field parameters include: the position of the intersection point of magnetic flux lines, elongation ratio, three angles, center safety factor, center magnetic field and the lateral radius of the tokamak device.

5. The method according to claim 1, wherein The wave trajectory equation includes: Where r is the radial coordinate, φ is the annular angle, z is the vertical coordinate, n = kc / ω = (n r ,m=rn φ ,n z ), k is the wave vector, c is the speed of light, ω is the angular frequency of the wave; D0 is the dispersion function calculated using the cold plasma approximation; D0=D r +iD i ;D r is the real part of plasma dispersion; D i is the imaginary part of plasma dispersion; D0=An 4 -Bn 2 +C=0; A=Ssin 2 θ+Pcos 2 I; B=RLsin 2 θ+PS(1+cos 2 i); C = PRL; is the cyclotron frequency of the particle, B is the total magnetic field, ω s is the plasma frequency, q s is the charge of the particle, m s is the mass of the particle, n s0 is the central density of the particle, ∈0 is the conductivity, θ is the angle between the magnetic field and the wave vector; r, φ, z, k r ,n φ ,k z The information that constitutes the point; k r is the radial wave vector.

6. The method according to claim 1 or 5, wherein: Also includes: determining the absorption power of the electron cyclotron wave at each of the points based on the imaginary part of the dispersion function in the wave trajectory equation; The absorbed power is calculated based on the following formula: Where, P(l) is the absorbed power, k r is the radial wave vector; D r is the real part of plasma dispersion; D i is the imaginary part of plasma dispersion; P(0) is the total absorbed power.

7. A device for identifying electron cyclotron wave trajectories, characterized in that: include: The magnetic flux module is used to obtain the normalized equilibrium poloidal magnetic flux under various preset plasma configurations; a magnetic field module for determining a total magnetic field based on the normalized balanced poloidal flux; a trajectory point module, configured to perform full interval integration of the wave trajectory equation in cylindrical coordinates based on the total magnetic field and using a fourth-order Runge-Kutta method with an indefinite step size, so as to sequentially determine information of each point in the propagation trajectory of the electron cyclotron wave; The normalized equilibrium poloidal flux is expressed by the following formula: Where ψ(r,z) is the normalized equilibrium poloidal flux, r is the radial coordinate, z is the vertical coordinate, ψ0 is the magnetic axis position, R0 is the lateral radius of the tokamak device, R x , E and τ are the position of the intersection of magnetic flux lines, the elongation ratio and the three angles respectively.

8. A terminal, characterized in that: include: A processor and a memory, wherein the memory stores a computer program, and when the processor executes the computer program, the method according to any one of claims 1 to 6 is implemented.

9. A storage medium, characterized in that: The storage medium stores a computer program, which implements the method according to any one of claims 1 to 6 when executed.

Images