CN112966453B - Simulation method for EAST tokamak radio frequency waveguide induced electronic temperature evolution - Google Patents

Simulation method for EAST tokamak radio frequency waveguide induced electronic temperature evolution Download PDF

Info

Publication number
CN112966453B
CN112966453B CN202110273883.3A CN202110273883A CN112966453B CN 112966453 B CN112966453 B CN 112966453B CN 202110273883 A CN202110273883 A CN 202110273883A CN 112966453 B CN112966453 B CN 112966453B
Authority
CN
China
Prior art keywords
evolution
radio frequency
electron temperature
electron
temperature
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.)
Active
Application number
CN202110273883.3A
Other languages
Chinese (zh)
Other versions
CN112966453A (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.)
Dalian University of Technology
Original Assignee
Dalian University of Technology
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 Dalian University of Technology filed Critical Dalian University of Technology
Priority to CN202110273883.3A priority Critical patent/CN112966453B/en
Publication of CN112966453A publication Critical patent/CN112966453A/en
Application granted granted Critical
Publication of CN112966453B publication Critical patent/CN112966453B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • G06F30/28Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/11Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
    • G06F17/13Differential equations
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/10Geometric CAD
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00Details relating to CAD techniques
    • G06F2111/10Numerical modelling
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2113/00Details relating to the application field
    • G06F2113/08Fluids

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Mathematical Physics (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Computational Mathematics (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Algebra (AREA)
  • Data Mining & Analysis (AREA)
  • Fluid Mechanics (AREA)
  • Computing Systems (AREA)
  • Operations Research (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • Plasma Technology (AREA)

Abstract

The invention discloses a simulation method for EAST tokamak radio frequency waveguide induced electron temperature evolution. When numerical simulation radio frequency wave is injected into an EAST tokamak device to cause an electronic temperature evolution process, the form of a heating source is determined according to hardware parameters of a radio frequency system, the evolution of the disturbance electronic temperature under a given magnetic field configuration is calculated according to an electronic transport equation under the tokamak configuration, a plasma resistor is corrected according to the evolved electronic temperature, the evolution of the magnetic field configuration is calculated by adopting the corrected resistor, the evolution of the disturbance electronic temperature is continuously calculated according to the magnetic field configuration, and repeated calculation is carried out to achieve the simulation of long-time evolution. The invention realizes the evolution of the electron temperature caused by simulating the radio frequency wave under the EAST tokamak real three-dimensional magnetic field configuration, obtains the distribution condition of the disturbance electron temperature at any moment in the three-dimensional space, makes up the defect of low precision of experimental measurement time, and has high calculation efficiency and strong numerical stability.

Description

Simulation method for EAST tokamak radio frequency waveguide induced electronic temperature evolution
Technical Field
The invention relates to numerical simulation of tokamak device discharge in the field of magnetic confinement controllable nuclear fusion, in particular to a simulation method for EAST tokamak radio frequency waveguide induced electron temperature evolution.
Background
The energy problem is a troublesome problem which is commonly faced by all people at present. At present, the energy used by people is mainly fossil fuel energy. Fossil fuels have a limited energy reserve and produce a large amount of harmful gases during combustion, thereby causing serious environmental pollution problems. To solve the energy problem fundamentally, the current method accepted by the scientific community is controllable nuclear fusion. In order to realize the controllable nuclear fusion, many proposals have been made, among which the experimental device which is most promising to be realized first is the tokamak device in the magnetic confinement method. Scientists in China are actively investing in the research of tokamak. At present, EAST super loop (EAST) of Anhui fertilizer in China is a full-superconducting Tokamak device, and all operation parameters are in the front of the world. To achieve steady state operation of EAST tokamak devices, the magnetic fluid instability in the tokamak configuration must be controlled. Wherein the core plasma can be heated by means of radio frequency wave injection, and the resistivity of the plasma is reduced. This suppresses resistive magnetic fluid instability.
In the EAST Tokamak discharge experiment, when radio frequency waves are injected, the distribution of the temperature of disturbance electrons in the space needs to be detected through soft X-ray diagnostic equipment, so that the instability of the magnetic fluid is better controlled. However, the electron temperature distribution measured by this method has limited space-time accuracy, and particularly, the time accuracy is not high enough, and it is difficult to satisfy the requirement of real-time measurement. Therefore, the three-dimensional spatial distribution of the temperature of the disturbance electrons at any time needs to be calculated by means of auxiliary processing of numerical simulation. Meanwhile, the magnetic field configuration in the experiment evolves along with time, and the evolution of the electronic temperature profile is closely related to the magnetic field configuration. Therefore, a simulation method capable of simulating the self-consistent nonlinear evolution of the temperature of the disturbance electrons caused by the injection of the radio frequency wave along the true three-dimensional magnetic field configuration along with time is needed. The method provided by the invention can exactly meet the requirement, can accurately describe the three-dimensional space distribution condition of the electronic temperature profile at any moment in higher time precision, has high calculation efficiency and strong numerical stability, and is an accurate and efficient numerical simulation method.
Disclosure of Invention
The invention aims to make up the defect of low time precision when the three-dimensional space distribution of the disturbed electron temperature is caused by measuring radio frequency waves in experiments, realize the nonlinear evolution of the simulated electron temperature under the EAST tokamak real three-dimensional magnetic field configuration, and obtain the distribution condition of the disturbed electron temperature in the three-dimensional space at any time with high time precision.
The technical scheme of the invention is as follows:
a simulation method for EAST tokamak radio frequency waveguide induced electron temperature evolution realizes self-consistent nonlinear evolution of an electron temperature profile under EAST real magnetic field configuration, can describe the space distribution condition of disturbance electron temperature by adopting a mode with higher time precision, has high calculation efficiency and strong numerical stability, and specifically comprises the following steps:
step 1: according to the geometric configuration of the EAST tokamak device, grid division is carried out on a core high-temperature plasma region in a discharge experiment, and physical quantities generated in the numerical simulation process are stored through nodes obtained through grid division;
step 2: according to the RF system used by the EAST Tokamak device, a Gaussian function is used to describe the RF wave heating source S ec The concrete form is as follows:
Figure BDA0002975707620000021
wherein R is a horizontal coordinate, Z is a vertical coordinate,
Figure BDA0002975707620000022
as a circumferential coordinate, S 0 The amplitude of the heating source can be calculated according to the transmitting power of a radio frequency system, R 0 、Z 0 And
Figure BDA0002975707620000023
respectively the horizontal direction, the vertical direction and the circumferential direction of the RF wave action region, Delta R d 、ΔZ d And
Figure BDA0002975707620000024
respectively horizontal direction, vertical direction and circumferential action width.
Initializing relevant parameters of a heating source according to experimental parameters such as the physical size of an antenna of the radio frequency system, the transmitting power of the radio frequency system and the like, wherein the relevant parameters comprise a spatial distribution function, an intensity amplitude value, a heating area and the like of the heating source. And initializing basic parameters such as a transport coefficient in an electron heat transport equation under the Tokamak configuration and a Sibirez resistance coefficient in a Sibirez resistance equation according to the experimental discharge parameters. Simultaneously obtaining the disturbance electron temperature delta T at the initial moment e (0)
And 3, step 3: the initial magnetic field configuration (namely the distribution information of the magnetic flux function in the three-dimensional space) is measured by adopting diagnostic equipment such as a magnetic detection ring, a magnetic flux ring and the like on an EAST Tokamak device, and a cubic spline interpolation method is usedThe measured magnetic flux function is converted to the grid divided in step 1, so as to obtain the magnetic flux function psi at the initial moment (0) And stored in the grid nodes.
And 4, step 4: calculating the evolution of the temperature of the disturbance electron along the magnetic field configuration at the current moment along with the time according to the electron heat transport equation under the Tokamak configuration, and calculating to obtain the temperature delta T of the disturbance electron at the next moment e (1)
The heat transport equation solved here is:
Figure BDA0002975707620000031
wherein, T e Is electron temperature and T e =T e0 +δT e ,T e0 And δ T e Equilibrium electron temperature and disturbance electron temperature, respectively, t is time, v is magnetofluid velocity, k // And kappa The transport coefficients of electrons along the parallel and vertical magnetic lines of force respectively,// and ^ subscript denote the parallel and vertical directions respectively, S ec Representing a heating source.
And 5: according to the electron temperature of the current moment obtained by calculating the electron heat transport equation in the step 4, calculating the plasma resistance eta of the current moment according to the sibez resistance equation (1)
The stezizer resistance equation solved here is:
Figure BDA0002975707620000032
where eta is the plasma resistivity, kappa S Is the sbez resistivity.
And 6: substituting the plasma resistance obtained in the step 5 into the magnetic fluid equation, numerically calculating the magnetic flux at the next moment and obtaining the magnetic flux psi after nonlinear evolution (1)
The evolution equation of the magnetic flux solved here is:
Figure BDA0002975707620000033
where ψ is the magnetic flux, φ is the electric potential, η is the plasma resistivity, j is the total plasma current, j is the plasma current bs Is the bootstrap current.
And 7: will disturb the electron temperature deltaT e Outputting the three-dimensional space distribution information to a file for storage;
and 8: calculating the magnetic field configuration at the current moment according to the magnetic flux obtained in the step 6, and repeating the steps 4-8 to obtain the three-dimensional space distribution delta T of the disturbance electron temperature at any moment e (n)
Has the advantages that: the method realizes the nonlinear evolution of the electron temperature profile after the injection of the numerical simulation radio frequency wave under the EAST tokamak real three-dimensional magnetic field configuration, can obtain the distribution condition of the disturbance electron temperature at any moment in a three-dimensional space, makes up the defect of low experimental measurement time precision, has high calculation efficiency and strong numerical stability, and is a stable and efficient numerical simulation method.
Drawings
FIG. 1 is a schematic cross-sectional view of an EAST tokamak apparatus to which the present invention is applicable.
FIG. 2 is a cross-sectional view of the initial magnetic field configuration used in the numerical simulation of the present invention.
FIG. 3(a) is a cross-sectional view of the simulation result of the temperature evolution of the perturbed electrons at a time 2ms after the injection of the RF wave in the EAST tokamak apparatus to which the present invention is applied.
FIG. 3(b) is a cross-sectional view of the simulation result of the temperature evolution of the perturbed electrons at 4ms after the injection of the RF wave in the EAST tokamak apparatus.
FIG. 4 is a main flow chart for numerically simulating the temperature evolution of the perturbed electrons according to the present invention.
Detailed Description
The following further describes a specific embodiment of the present invention with reference to the drawings and technical solutions.
The EAST tokamak device has a cross-sectional configuration as shown in fig. 1, where the core location is the region containing the plasma during the experimental discharge. EAST tokamak is equipped with a soft X-ray system to measure the spatial distribution of perturbed electron temperature. The core plasma region is first gridded according to the EAST geometry shown in fig. 1. In EAST tokamak discharge experiments, suppression of magnetic islands is one of the main effects of rf wave heating, so the initial magnetic field configuration usually includes a magnetic island structure, as shown in fig. 2, where magnetic lines of force are coupled to form a magnetic island chain structure including three magnetic islands in one period of polar direction. The rf wave injection site is generally centered on the magnetic island. The electron temperature of disturbance caused by radio frequency wave injection can be subjected to self-consistent nonlinear evolution along a specific magnetic field configuration, and finally, steady-state distribution is achieved in space. This embodiment describes the heating source using a gaussian function. EAST experiments usually use electron cyclotron emission systems to generate rf waves, and gaussian functions can better describe the heating source. Then, the heating source is initialized and the initial magnetic field information and the initial electron temperature spatial distribution information are stored in the grid nodes. And calculating the nonlinear evolution of the temperature of the disturbance electrons along a specific magnetic line according to an electron heat transport equation in the Tokamak configuration. And calculating the resistivity of the plasma at the moment according to the electron temperature at the current moment through a Spanish resistance equation. Substituting the calculated plasma resistivity into a magnetic flux evolution equation to calculate the correction of the magnetic field at the last moment by the injection of the radio frequency wave, and obtaining the corrected magnetic field. And finally, taking the corrected magnetic field as the magnetic field configuration at the current moment, calculating the nonlinear evolution of the temperature of the disturbed electrons, and repeating the previous steps for multiple times to obtain the three-dimensional stable distribution of the electron temperature at any moment with high time precision, as shown in fig. 3, wherein a cross section diagram of the spatial distribution of the disturbed electron temperature at two moments is shown.
The specific implementation steps are as follows:
step 1: according to the geometric configuration of the EAST tokamak device, grid division is carried out on a core high-temperature plasma region in a discharge experiment, and physical quantities generated in the numerical simulation process are stored through nodes obtained through grid division;
step 2: mining according to EAST tokamak apparatusRF system for use therein, where the RF wave heating source S is described by a Gaussian function ec The concrete form is as follows:
Figure BDA0002975707620000051
wherein R is a horizontal coordinate, Z is a vertical coordinate,
Figure BDA0002975707620000052
as a circumferential coordinate, S 0 The amplitude of the heating source can be calculated according to the transmitting power of a radio frequency system, R 0 、Z 0 And
Figure BDA0002975707620000053
the horizontal direction, the vertical direction and the circumferential direction of the radio frequency wave action region, delta R d 、ΔZ d And
Figure BDA0002975707620000054
respectively horizontal direction, vertical direction and circumferential action width.
Initializing relevant parameters of a heating source according to experimental parameters such as the physical size of an antenna of the radio frequency system, the transmitting power of the radio frequency system and the like, wherein the relevant parameters comprise a spatial distribution function, an intensity amplitude value, a heating area and the like of the heating source. And initializing basic parameters such as a transport coefficient in an electron heat transport equation under the Tokamak configuration and a Sibirez resistance coefficient in a Sibirez resistance equation according to the experimental discharge parameters. The specific parameters used in this example are as follows:
S 0 =1.5MW、R 0 =1.9m、Z 0 =0.3m、
Figure BDA0002975707620000055
ΔR d =0.1m、ΔZ d 0.1m and
Figure BDA0002975707620000056
in this embodiment, the RF system is turned on at the initial time becauseThe temperature delta T of the disturbance electrons caused by the radio frequency system at the initial moment e (0) Is zero;
and step 3: measuring initial magnetic field configuration (namely distribution information of the magnetic flux function in a three-dimensional space) by adopting diagnostic equipment such as a magnetic probe ring and a magnetic flux ring on an EAST tokamak device, and converting the measured magnetic flux function to the grid divided in the step 1 by a cubic spline interpolation method, thereby obtaining the magnetic flux function psi at the initial moment (0) And stored in grid nodes, the magnetic flux function psi (0) The contour plot is shown in figure 2.
And 4, step 4: calculating the evolution of the temperature of the disturbance electron along the magnetic field configuration at the current moment along with the time according to the electron heat transport equation under the Tokamak configuration, and calculating to obtain the temperature delta T of the disturbance electron at the next moment e (1)
The heat transport equation solved here is:
Figure BDA0002975707620000061
wherein, T e Is electron temperature and T e =T e0 +δT e ,T e0 And δ T e Equilibrium electron temperature and disturbance electron temperature, respectively, t is time, v is magnetofluid velocity, k // And kappa The transport coefficients of electrons along the parallel and vertical magnetic lines of force respectively,// and ^ subscript denote the parallel and vertical directions respectively, S ec Representing a heating source.
The specific parameters used in this example are as follows:
v~400m/s,κ // =1.6×10 6 m 2 s and kappa =1.6×10 -3 m 2 /s,S ec Can be calculated by the data in step 2.
And 5: according to the electron temperature of the current moment obtained by calculating the electron heat transport equation in the step 4, calculating the plasma resistance eta of the current moment according to the sibez resistance equation (1)
The stezizer resistance equation solved here is:
Figure BDA0002975707620000062
where eta is the plasma resistivity, kappa S When the value is 1.0, the sbuzer resistivity is obtained.
Step 6: substituting the plasma resistance obtained in the step 5 into the magnetic fluid equation, numerically calculating the magnetic flux at the next moment and obtaining the magnetic flux psi after nonlinear evolution (1)
The evolution equation of the magnetic flux solved here is:
Figure BDA0002975707620000063
where ψ is the magnetic flux, φ is the potential, the potential is measured by Langmuir probe, η is the plasma resistivity, j is the total plasma current, j is the plasma resistivity bs For bootstrap currents, the current profile can be obtained by inversion of the equilibrium profile.
And 7: will disturb the electron temperature deltaT e Outputting the three-dimensional space distribution information to a file for storage;
and 8: calculating the magnetic field configuration at the current moment according to the magnetic flux obtained in the step 6, and repeating the steps 4-8 to obtain the three-dimensional space distribution delta T of the disturbance electron temperature at any moment e (n)
In this example, dt is 1 × 10 -9 s time steps for simulation (i.e. 1X 10 simulation time passes for each round of 4-8 steps) -9 s), passing through 2X 10 6 The round simulation can obtain the disturbance electron temperature space distribution at the time t-2 ms, as shown in fig. 3(a), and then 2 × 10 6 The round simulation can obtain the situation of the disturbance electron temperature space distribution at the moment t-4 ms, as shown in fig. 3 (b).
The above description is further detailed in connection with the preferred embodiments of the present invention, and it is not intended to limit the practice of the invention to these descriptions. It will be apparent to those skilled in the art that various modifications, additions, substitutions, and the like can be made without departing from the spirit of the invention.

Claims (1)

1. A simulation method for EAST tokamak radio frequency waveguide-induced electron temperature evolution realizes self-consistent nonlinear evolution of an electron temperature profile after radio frequency wave injection under a real EAST tokamak three-dimensional magnetic field configuration, obtains an electron temperature profile at any moment and provides an electron temperature spatial distribution condition; the method is characterized by comprising the following steps:
step 1: according to the geometric configuration of the EAST tokamak device, grid division is carried out on a core high-temperature plasma region in a discharge experiment, and physical quantities generated in the numerical simulation process are stored through nodes obtained through grid division;
step 2: according to the radio frequency system adopted by the EAST tokamak device, a Gaussian function is adopted to describe the radio frequency wave heating source S ec The concrete form is as follows:
Figure FDA0002975707610000011
wherein, R is a horizontal coordinate; z is a vertical coordinate;
Figure FDA0002975707610000012
is a circumferential coordinate; s 0 Calculating the amplitude of the heating source according to the transmitting power of the radio frequency system; r 0 、Z 0 And
Figure FDA0002975707610000013
respectively the horizontal direction, the vertical direction and the circumferential direction of the RF wave action region, Delta R d 、ΔZ d And
Figure FDA0002975707610000014
the action widths in the horizontal direction, the vertical direction and the annular direction are respectively;
initializing relevant parameters of a heating source according to the physical size of an antenna of the radio frequency system and the transmitting power of the radio frequency system, wherein the relevant parameters comprise a spatial distribution function, an intensity amplitude and a heating area of the heating source; initializing a transport coefficient in an electronic heat transport equation under a Tokamak potential form and a Sibirez resistance coefficient in a Sibirez resistance equation according to experimental discharge parameters; simultaneously obtaining the disturbance electron temperature delta T at the initial moment e (0)
And step 3: measuring the initial magnetic field configuration, namely the distribution information of the magnetic flux function in the three-dimensional space by adopting a magnetic detection ring and a magnetic flux ring on an EAST tokamak device, and converting the measured magnetic flux function to the grids divided in the step 1 by a cubic spline interpolation method, thereby obtaining the magnetic flux function psi at the initial moment (0) And stored in the grid nodes;
and 4, step 4: calculating the evolution of the temperature of the disturbance electron along the magnetic field configuration at the current moment along with the time according to the electron heat transport equation under the Tokamak configuration, and calculating to obtain the temperature delta T of the disturbance electron at the next moment e (1)
The heat transport equation solved is:
Figure FDA0002975707610000021
wherein, T e Is electron temperature and T e =T e0 +δT e ,T e0 And δ T e Equilibrium electron temperature and disturbance electron temperature, respectively, t is time, v is magnetofluid velocity, k // And kappa The transport coefficients of electrons along the parallel and vertical magnetic lines of force respectively,// and ^ subscript denote the parallel and vertical directions respectively, S ec Represents a heating source;
and 5: according to the electron temperature of the current moment obtained by calculating the electron heat transport equation in the step 4, calculating the plasma resistance eta of the current moment according to the sibez resistance equation (1)
The stezizer resistance equation solved is:
Figure FDA0002975707610000022
where eta is the plasma resistivity, kappa S Has a Sibirez resistivity;
step 6: substituting the plasma resistivity obtained in the step 5 into a magnetic flux evolution equation, calculating the magnetic flux at the next moment and obtaining the magnetic flux psi after nonlinear evolution (1)
The magnetic flux evolution equation solved is:
Figure FDA0002975707610000023
where ψ is the magnetic flux, φ is the electric potential, η is the plasma resistivity, j is the total plasma current, j is the plasma current bs Is a bootstrap current;
and 7: will disturb the electron temperature deltaT e Outputting the three-dimensional space distribution information to a computer for storage;
and 8: calculating the magnetic field configuration at the current moment according to the magnetic flux obtained in the step 6, and then repeating the steps 4-8 to obtain the disturbance electron temperature three-dimensional space distribution delta T at any moment e (n)
CN202110273883.3A 2021-03-15 2021-03-15 Simulation method for EAST tokamak radio frequency waveguide induced electronic temperature evolution Active CN112966453B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110273883.3A CN112966453B (en) 2021-03-15 2021-03-15 Simulation method for EAST tokamak radio frequency waveguide induced electronic temperature evolution

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110273883.3A CN112966453B (en) 2021-03-15 2021-03-15 Simulation method for EAST tokamak radio frequency waveguide induced electronic temperature evolution

Publications (2)

Publication Number Publication Date
CN112966453A CN112966453A (en) 2021-06-15
CN112966453B true CN112966453B (en) 2022-09-20

Family

ID=76278976

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110273883.3A Active CN112966453B (en) 2021-03-15 2021-03-15 Simulation method for EAST tokamak radio frequency waveguide induced electronic temperature evolution

Country Status (1)

Country Link
CN (1) CN112966453B (en)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114896854B (en) * 2022-05-12 2024-06-14 大连理工大学 Simulation method for plasma radiation evolution in tokamak
CN116167247B (en) * 2023-04-21 2023-07-11 中国科学院合肥物质科学研究院 GS equation numerical calculation method based on Fengha Gnoff method

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108733946A (en) * 2018-05-24 2018-11-02 大连理工大学 A kind of analogy method to develop for electron cyclotron Current profile in tokamak
CN111259599A (en) * 2020-01-07 2020-06-09 大连理工大学 Test particle simulation method for analyzing constraint performance of high-energy particles in tokamak

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108733946A (en) * 2018-05-24 2018-11-02 大连理工大学 A kind of analogy method to develop for electron cyclotron Current profile in tokamak
CN111259599A (en) * 2020-01-07 2020-06-09 大连理工大学 Test particle simulation method for analyzing constraint performance of high-energy particles in tokamak

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
基于Corsica的EAST等离子体平衡计算;程曹军等;《合肥工业大学学报(自然科学版)》;20110428(第04期);全文 *

Also Published As

Publication number Publication date
CN112966453A (en) 2021-06-15

Similar Documents

Publication Publication Date Title
WO2022193041A1 (en) Method for simulating electron temperature evolution caused by east tokamak radio frequency wave
CN112966453B (en) Simulation method for EAST tokamak radio frequency waveguide induced electronic temperature evolution
CN108733946B (en) Simulation method for electron cyclotron current profile evolution in tokamak
CN107657137B (en) Fractional order electromagnetic anomalous diffusion three-dimensional simulation method for rational function approximation
CN110276109B (en) Simulation method for electromagnetic property of plasma sheath of hypersonic aircraft
Speckhofer et al. Experimental and theoretical investigation of high-pressure arcs. II. The magnetically deflected arc (three-dimensional modeling)
Neiman et al. On the question of taking into account the main dimensions when selecting type of electromagnet according to value of constructive factor
Liu et al. Two‐step method for identifying photovoltaic grid‐connected inverter controller parameters based on the adaptive differential evolution algorithm
CN114626268B (en) High-precision time domain calculation method for strong electromagnetic pulse propagation process
Gurusinghe et al. Efficient algorithms for real‐time monitoring of transmission line parameters and their performance with practical synchrophasors
Xia et al. The influence of magnetic circuit and operating parameters on the plasma property of 2 cm ECRIT ion source
Dan et al. Parameters estimation of horizontal multilayer soils using a heuristic algorithm
Tinguely et al. High-resolution disruption halo current measurements using Langmuir probes in Alcator C-Mod
Sanpei et al. A data-assimilation based method for equilibrium reconstruction of magnetic fusion plasma and its application to reversed field pinch
CN117013561A (en) New energy system weak point positioning method and device based on port time domain characteristics
WO2024027015A1 (en) Grounding grid topological structure measurement method under extremely cold conditions
Yang et al. High‐Temperature Superconducting Cable Fault Location Method Based on Improved Time‐Frequency Domain Reflection Method and EEMD Noise Reduction
Eichhorn et al. On quench propagation, quench detection and second sound in SRF cavities
Xu et al. Optimisation of phase factor in phase retrieval for reflector antennas with active surface
CN106485071A (en) A kind of quick far to field transformation method of many layer packet structures
Huang et al. Implementation of a Kalman filter-based eddy current estimator for the P-EFIT magnetic equilibrium reconstruction code
Frerichs et al. Field line reconstruction for edge transport modeling in non-axisymmetric tokamaks configurations
Liu et al. Modelling and application of a new method to measure the non-thermal electron current in the edge of magnetically confined plasma
Bettini et al. Dynamic identification of plasma magnetic contour in fusion machines
CN110579790B (en) Method for evaluating electrostatic effect of self-powered neutron detector

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