CN115512791A - Particle evolution simulation method based on multi-component plasma radiation effect in tokamak - Google Patents
Particle evolution simulation method based on multi-component plasma radiation effect in tokamak Download PDFInfo
- Publication number
- CN115512791A CN115512791A CN202211133612.9A CN202211133612A CN115512791A CN 115512791 A CN115512791 A CN 115512791A CN 202211133612 A CN202211133612 A CN 202211133612A CN 115512791 A CN115512791 A CN 115512791A
- Authority
- CN
- China
- Prior art keywords
- density distribution
- radiation
- plasma
- evolution
- impurity
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 44
- 239000002245 particle Substances 0.000 title claims abstract description 18
- 238000004088 simulation Methods 0.000 title claims abstract description 16
- 230000000191 radiation effect Effects 0.000 title claims abstract description 8
- 238000009826 distribution Methods 0.000 claims abstract description 151
- 230000005855 radiation Effects 0.000 claims abstract description 107
- 239000012535 impurity Substances 0.000 claims abstract description 71
- 150000002500 ions Chemical class 0.000 claims abstract description 50
- 238000004364 calculation method Methods 0.000 claims abstract description 31
- 238000012937 correction Methods 0.000 claims abstract description 14
- 239000011553 magnetic fluid Substances 0.000 claims abstract description 4
- 230000004907 flux Effects 0.000 claims description 30
- 238000009792 diffusion process Methods 0.000 claims description 12
- 238000001228 spectrum Methods 0.000 claims description 11
- 230000005467 cyclotron radiation Effects 0.000 claims description 7
- 238000002474 experimental method Methods 0.000 claims description 7
- 230000008569 process Effects 0.000 claims description 5
- 230000008878 coupling Effects 0.000 claims description 4
- 238000010168 coupling process Methods 0.000 claims description 4
- 238000005859 coupling reaction Methods 0.000 claims description 4
- 238000001816 cooling Methods 0.000 claims description 3
- 230000004927 fusion Effects 0.000 abstract description 7
- 230000000694 effects Effects 0.000 abstract description 2
- 210000002381 plasma Anatomy 0.000 description 60
- 241000282414 Homo sapiens Species 0.000 description 2
- 230000003993 interaction Effects 0.000 description 2
- 230000007246 mechanism Effects 0.000 description 2
- 241000894007 species Species 0.000 description 2
- 230000005461 Bremsstrahlung Effects 0.000 description 1
- 230000002159 abnormal effect Effects 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 238000007792 addition Methods 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000006866 deterioration Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000004083 survival effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
- G16C60/00—Computational materials science, i.e. ICT specially adapted for investigating the physical or chemical properties of materials or phenomena associated with their design, synthesis, processing, characterisation or utilisation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
- G06T17/20—Finite element generation, e.g. wire-frame surface description, tesselation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/06—Power analysis or power optimisation
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E30/00—Energy generation of nuclear origin
- Y02E30/10—Nuclear fusion reactors
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Computing Systems (AREA)
- Geometry (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Computer Hardware Design (AREA)
- Life Sciences & Earth Sciences (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Computer Graphics (AREA)
- Software Systems (AREA)
- Plasma Technology (AREA)
Abstract
The invention discloses a particle evolution simulation method based on multi-component plasma radiation effect in tokamak, and belongs to the technical field of magnetic confinement controlled nuclear fusion. The evolution of the impurity density distribution, the electron density distribution and the ion density distribution is firstly calculated, then the continuous radiation and the radiation of the radiation are calculated, and then various plasma radiations are accurately added in the magnetic island. And then calculating the magnetic field configuration under the plasma radiation correction, and continuously calculating the total plasma radiation and the evolution of related parameters under the advanced magnetic field configuration, so that the calculation is repeated, and the simulation effect of long-time evolution is achieved. The invention not only realizes the calculation of real-time dynamic distribution of impurities, electrons and ions, but also can calculate the evolution of various plasma radiations under the true three-dimensional magnetic field configuration so as to obtain the plasma related parameter profile at any moment; and simultaneously, the space distribution condition of plasma radiation and the instability condition of the magnetic fluid under the advanced magnetic field configuration are more accurately described.
Description
Technical Field
The invention belongs to the technical field of magnetic confinement controlled nuclear fusion, and particularly relates to a particle density distribution evolution method for multi-component plasma radiation effect in tokamak.
Background
In recent years, productivity and production technology have been effectively improved due to the rapid increase of the world population, and thus human consumption of energy has been rapidly increased, resulting in the increasing exhaustion of non-renewable fossil energy. Obviously, this has led to the gradual deterioration of the environment on which human beings rely for survival in recent years, which has resulted in the rapid aggravation of the greenhouse effect, and the creation of abnormal weather and severe natural phenomena. Therefore, the search and use of clean energy is a continuous pursuit of scientists. At present, nuclear energy is relatively efficient and pollution-free energy, and the nuclear energy strategy proposed in China is as follows: the thermal neutron reactor is built in the near term, the fast neutron reactor is used in the middle term, and the fusion reactor is built in the long term. The fusion reactor refers to magnetic confinement controlled nuclear fusion. To realize the magnetic confinement controlled nuclear fusion, the steady state operation mechanism is completely mastered, and meanwhile, the instability caused by the interaction between internal plasmas, the interaction between the plasmas and the wall of the tokamak and other complex problems which can promote energy loss to cause the discharge to be extinguished are fully understood and solved. Therefore, it is important to study the tokamak energy loss in the most advanced magnetic field configurations.
Advanced magnetic field configurations can more easily confine the plasma, thereby achieving stable long pulse discharge and ultimately its commercial value. In tokamak plasma, plasma radiation, which is one of the principal energy loss mechanisms of tokamak, refers to electromagnetic waves emitted from the plasma in each frequency band. The spectrum of the radiation mainly comprises a continuous spectrum and a line spectrum, and the corresponding radiation also comprisesIt is called continuous radiation and line radiation. The main representatives of continuous radiation are bremsstrahlung radiation, which mainly includes electron cyclotron radiation and characteristic radiation. The characteristic radiation is mainly from impurity ion radiation. In general, the radiation spectrum covers all bands from infrared to X-rays. In future large magnetically confined nuclear fusion devices, high z eff Is very important and has a strong influence on the magneto-hydrodynamic instability in the plasma. Therefore, it is important to study plasma radiation and clearly grasp characteristics of the plasma radiation to realize a long pulse discharge. Therefore, the invention provides a simulation method which can be used for coupling a plurality of plasma radiation effects together under the condition of the configuration of the advanced magnetic field, simultaneously considers the influence of the distribution evolution of impurity particles, electrons and ions on the total radiation of the plasma, thereby carrying out numerical simulation research on the aspect of magnetohydrodynamic instability in Tokamak discharge.
Disclosure of Invention
In order to fill the technical blank of numerical simulation of plasma radiation in tokamak, the invention provides a simulation method for particle density distribution evolution of plasma radiation in tokamak under an advanced magnetic field configuration, namely, particles used in plasma radiation are evolved along with time. The method can be well combined with experiments, real-time dynamic distribution of impurities, electrons and ions is obtained through calculation of impurity density distribution evolution, electron density distribution evolution and ion density distribution evolution, plasma radiation can be predicted and simulated more accurately, three-dimensional simulation of nonlinear evolution of plasma radiation is realized, further influence of the plasma radiation on Tokamak discharge under an advanced magnetic field configuration can be studied systematically, and meanwhile reference can be provided for experimental results of the plasma radiation.
The technical scheme adopted by the invention is as follows:
a particle evolution simulation method based on multi-component plasma radiation effects in tokamak realizes the evolution of electron, ion and impurity density distribution, further can obtain the nonlinear evolution of various plasma radiations under a real three-dimensional magnetic field configuration, and can also obtain the evolution condition of the plasma radiation at any time and the evolution section of related parameters of the plasma, thereby more accurately describing the spatial distribution condition of the tokamak plasma under the combined action of the various plasma radiations under an advanced magnetic field configuration. The method specifically comprises the following steps:
step 1: the plasma area in the discharge experiment of the tokamak device is divided into grids, and the total plasma radiation value, the magnetic flux function value and the like obtained in the plasma radiation evolution process can be stored by the divided grid nodes.
Step 2: respectively calculating initial impurity density distribution according to highest ionization state, far infrared interferometer and Thomson scattering in experiment
Initial electron density distribution
And initial electron temperature distribution
And step 3: obtaining the configuration of the initial magnetic field by adopting equipment such as a magnetic flux loop and the like in Tokamak discharge, and calculating the initial magnetic flux psi by a numerical simulation method and the like (0) And stored in the grid nodes.
And 4, step 4: the initial magnetic flux psi (0) Initial impurity density distribution
Initial electron density distribution
And initial ion density distribution
Respectively carrying out calculation in a magnetofluid equation, an impurity density distribution evolution equation, an electron density distribution evolution equation and an ion density distribution evolution equation to respectively obtain the magnetic flux psi at the next moment (1) Impurity density distribution
Electron density distribution
And ion density distribution
The method comprises the following specific steps:
the evolution equation of the impurity density distribution is:
the evolution equation of the electron density distribution is:
the evolution equation for the ion density distribution is:
wherein n is z Is the particle density distribution of the impurity, n e Denotes the electron density distribution, n i Represents an ion density distribution; phi is the potential and t is the time; d z|| And D z⊥ Impurity diffusion coefficients in the parallel direction and the perpendicular direction of the impurity, respectively; d e|| And D e⊥ Impurity diffusion coefficients in the parallel direction and the perpendicular direction of electrons, respectively; d i|| And D i⊥ Impurity diffusion coefficients in the parallel direction and the perpendicular direction of the ions, respectively;
and
gradient operators in the parallel direction and the vertical direction respectively; s radz Is a source term of impurities, S rade Being a source item of electrons, S radi Is the source term of the ions.
The specific calculation method of the evolution equation of the impurity density distribution, the electron density distribution and the ion density distribution is as follows:
step 4.1: the impurity density distribution, electron density distribution and ion density distribution are expressed by a spectrum method:
wherein, (m, n) is the module of the circumferential direction and the polar direction; r 0 Is the large radius of tokamak; theta is a polar angle; z represents the column direction; f is the particle distribution, including impurity density distribution, electron density distribution, and ion density distribution.
Step 4.2: the time advance calculation of various particle distributions was performed using a two-step prediction-correction method. The calculation format of the two-step prediction-correction method is as follows:
and (3) prediction:
and (3) correction:
wherein H represents a magnetic flux and a magnetic currentParameters in the volume equation; v is the diffusion coefficient; t is time;
represents half a time step; dt represents a time step; the subscript rhs represents the right-hand term of the magnetofluid equation,
representing the gradient in the direction of the perpendicular magnetic field.
The evolution equations of impurity density distribution, electron density distribution and ion density distribution can be solved through the steps 4.1 and 4.2. The reason for adopting the spectrum method and the prediction-correction method is that the spectrum method and the prediction-correction method can more accurately, quickly and stably calculate the impurity density distribution, the electron density distribution and the ion density distribution at each moment, further couple the evolution equation of the impurity density distribution, the evolution equation of the electron density distribution and the evolution equation of the ion density distribution with the magnetofluid equation, and can calculate the distribution of impurities, electrons and ions at any moment in a self-consistent manner.
And 5: calculating the impurity density distribution according to the step 4, and when the impurity density distribution of the Tokamak core reaches a set threshold value, further opening a calculation module of plasma radiation, thereby calculating the evolution of the magnetic field configuration along with the time and obtaining the total radiation of the plasma at the moment
The method comprises the following specific steps:
step 5.1: according to electron temperature distribution T e The corresponding radiation cooling rate L (T) is calculated according to the value of (A) and the generated impurity species e ) The numerical value of (c).
And step 5.2: calculating the total radiation of the plasma, including continuous radiation and line radiation.
According to the formula of continuous radiation:
and electron cyclotron radiation formula:
P c ＝5.4×10 -25 n e B 2 T e erg s -1 cm -3
and the calculation formula of the characteristic line radiation:
P z ＝n e n z ∑L(T e )erg s -1 cm -3
calculating the magnitude of the three kinds of radiation, and then accumulating the three kinds of radiation to obtain the total radiation of the plasma
Wherein, P b Representing the magnitude of the continuous radiation, P c Representing the magnitude of electron cyclotron radiation, P z Represents the magnitude of the line radiation; z is a radical of eff Is the effective charge distribution.
And 6: coupling the calculated plasma total radiation into a magnetofluid equation, calculating the evolution of magnetic flux at a unit time step length, and obtaining the magnetic flux added with the evolution of the plasma total radiation
And 7: the total radiation of the plasma obtained by the calculation in the step 5 is used
And outputting the three-dimensional space distribution information.
And 8: according to the magnetic flux which is obtained after the calculation in the step 6 and takes the total radiation of the plasma into consideration
Further calculating the current magnetic field configuration, and further continuously repeating the steps 5-8 to obtain the total plasma radiation at any moment
And magnetic flux after total radiation of plasma
The invention has the beneficial effects that: the method can calculate the electron density distribution evolution, the ion density distribution evolution and the impurity density distribution evolution, can represent the evolution process of each parameter in the tokamak plasma in real time, and can calculate the radiation of various plasmas, thereby being capable of carrying out numerical simulation and analysis on the instability of the magnetofluid in the discharge of the tokamak device under the advanced magnetic field configuration, providing reference for the experimental result of the radiation of the plasmas, having high calculation efficiency and strong numerical stability, and being a stable and efficient numerical simulation method.
Drawings
Fig. 1 is a schematic diagram of a three-dimensional real magnetic field configuration of a tokamak experimental apparatus suitable for the present invention, including total radiation of magnetic islands and plasma.
FIG. 2 is a calculation result of the spatial distribution of the total plasma radiation over time, which is calculated by the method after the evolution of impurities, electrons and ions under the configuration of the advanced magnetic field; wherein (a) is 55 tau A Spatial distribution of radiation in time, (b) is 60. Tau A Spatial distribution of radiation in time, (c) is 65 τ A Spatial distribution of the radiation in time, (d) is 75. Tau A Spatial distribution of radiation. Wherein tau is A Is the alfen time.
FIG. 3 is a flow chart of the present invention for simulating and calculating various plasma radiation effects under advanced magnetic field configurations.
Detailed Description
In order to make the objects, technical solutions and advantages of the present invention more apparent, embodiments of the present invention will be described in detail below with reference to the accompanying drawings.
The tokamak device is a toroidal device shaped like a tire, schematically shown in fig. 1, in which magnetic surfaces are nested one upon another, and magnetic lines of force are wound in circumferential and polar directions, which may generate magnetic islands in case of instability. Fig. 2 shows that the spatial scale of plasma radiation is continuously evolved and is increased along with the increase of the width of the magnetic island, and simultaneously, the impurity density distribution, the electron density distribution and the ion density distribution are evolved along with time, so that the method not only can calculate the evolution process of the impurity density distribution, the electron density distribution and the ion density distribution in tokamak along with time, but also can completely add various plasma radiations into the magnetic island, so that the plasma radiations are self-consistent and nonlinearly evolved along with the evolution of the width of the magnetic island.
The specific implementation steps are as follows:
step 1: the plasma area in the discharge experiment of the tokamak device is divided into grids, and the total plasma radiation value, the magnetic flux function value and the like obtained in the plasma radiation evolution process can be stored by the divided grid nodes.
Step 2: respectively calculating initial impurity density distribution according to highest ionization state, far infrared interferometer and Thomson scattering in experiment
Initial electron density distribution
And initial electron temperature distribution
And step 3: the configuration of an initial magnetic field is obtained by adopting a magnetic flux loop in Tokamak discharge, and the initial magnetic flux psi is obtained by calculating through a numerical simulation method (0) And stored in the grid nodes.
And 4, step 4: will initiate a magnetic flux v (0) Initial impurity density distribution
Initial electron density distribution
And initial ion density distribution
Respectively carrying out calculation in a magnetofluid equation, an impurity density distribution evolution equation, an electron density distribution evolution equation and an ion density distribution evolution equation to obtain the magnetic flux psi at the next moment (1) Impurity density distribution
Electron density distribution
And ion density distribution
The specific calculation method of the evolution equation of the impurity density distribution, the electron density distribution and the ion density distribution is as follows:
the evolution equation for the impurity density distribution is:
the evolution equation of the electron density distribution is:
the evolution equation for the ion density distribution is:
wherein n is z Is the particle density distribution of the impurity, n e Denotes the electron density distribution, n i Represents an ion density distribution; phi is the potential and t is the time; d z|| And D z⊥ Impurity diffusion coefficients in the parallel direction and the perpendicular direction of the impurity, respectively; d e|| And D e⊥ Impurity diffusion coefficients in the parallel direction and the perpendicular direction of electrons, respectively; d i|| And D i⊥ Impurity diffusion coefficients in the parallel direction and the perpendicular direction of the ions, respectively;
and
gradient operators in the parallel direction and the perpendicular direction respectively; s radz Is a source term of impurities, is a source term of electrons, S radi Is the source term of the ions.
The specific calculation method of the evolution equation of the impurity density distribution, the electron density distribution and the ion density distribution is as follows:
step 4.1: the impurity density distribution, electron density distribution and ion density distribution are expressed by a spectrum method:
wherein, (m, n) is the module of the circumferential direction and the polar direction; r is 0 Is the large radius of tokamak; theta is a polar angle; z represents the column direction; f is the particle distribution, including impurity density distribution, electron density distribution, and ion density distribution.
And 4.2: the time advance calculation of various particle distributions was performed using a two-step prediction-correction method. The calculation format of the two-step prediction-correction method is as follows:
and (3) prediction:
and (3) correction:
wherein H represents the magnetic flux and parameters in the magnetic fluid equation; v is the diffusion coefficient; t is time;
represents half a time step; dt represents a time step; the subscript rhs represents the right-hand term of the magnetofluid equation,
representing the gradient in the direction of the perpendicular magnetic field.
The evolution equations of impurity density distribution, electron density distribution and ion density distribution can be solved through the steps 4.1 and 4.2. The reason for adopting the spectrum method and the prediction-correction method is that the spectrum method and the prediction-correction method can more accurately, quickly and stably calculate the impurity density distribution, the electron density distribution and the ion density distribution at each moment, further couple the evolution equation of the impurity density distribution, the evolution equation of the electron density distribution and the evolution equation of the ion density distribution with the magnetofluid equation, and can calculate the distribution of impurities, electrons and ions at any moment in a self-consistent manner.
And 5: the impurity density distribution is calculated according to step 4, when the impurity density distribution of the tokamak core reaches a set threshold value (which needs to be determined according to the specific study, case, for example, n) z ＝1×10 -5 n e ) And then opening a calculation module of plasma radiation, thereby calculating the evolution of the magnetic field configuration along with the time and obtaining the total plasma radiation at the moment
The method comprises the following specific steps:
step 5.1: according to electron temperature distribution T e The corresponding radiation cooling rate L (T) is calculated according to the value of (A) and the generated impurity species e ) The numerical value of (c).
Step 5.2: calculating the total radiation of the plasma, including continuous radiation and line radiation.
According to the formula of continuous radiation:
and electron cyclotron radiation formula:
P c ＝5.4×10 -25 n e B 2 T e erg s -1 cm -3
and the calculation formula of the characteristic line radiation:
P z ＝n e n z ∑L(T e )erg s -1 cm -3
calculating the magnitude of the three kinds of radiation, and then accumulating the three kinds of radiation to obtain the total radiation of the plasma
Wherein, P b Representing the magnitude of the continuous radiation, P c Representing the magnitude of electron cyclotron radiation, P z Representing the magnitude of the radiation of the line, z eff Is the effective charge distribution.
Step 6: coupling the calculated plasma total radiation into a magnetofluid equation, calculating the evolution of magnetic flux under a unit time step length, and obtaining the magnetic flux added with the plasma total radiation evolution
And 7: the total radiation of the plasma obtained by the calculation in the step 5 is used
And outputting the three-dimensional space distribution information.
And step 8: according to the magnetic flux which is obtained after the calculation in the step 6 and takes the total radiation of the plasma into consideration
Further calculating the current magnetic field configuration, and further continuously repeating the steps 5-8 to obtain the total plasma radiation at any moment
And magnetic flux after total radiation of plasma
The above is a detailed description of an example of the present invention for numerical simulation calculation of multiple plasma radiations in an advanced magnetic field configuration tokamak, and the specific implementation of the present invention is not considered to be limited 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 (3)
1. A particle evolution simulation method based on multi-component plasma radiation effect in tokamak is characterized by comprising the following steps:
step 1: carrying out grid division on a plasma region in a discharge experiment of the tokamak device, and storing data of a total plasma radiation value and a magnetic flux function value obtained in the plasma radiation evolution process by the divided grid nodes;
and 2, step: respectively calculating initial impurity density distribution according to highest ionization state, far infrared interferometer and Thomson scattering in experiment
Initial electron density distribution
And initial electron temperature distribution
And step 3: obtaining the configuration of an initial magnetic field by adopting a magnetic flux loop in Tokamak discharge, and calculating to obtain an initial magnetic flux psi (0) And stored in the grid nodes;
and 4, step 4: the initial magnetic flux psi (0) Initial impurity density distribution
Initial electron density distribution
And a given initial ion density distribution
Respectively carrying into a magnetic fluid equation, an impurity density distribution evolution equation, an electron density distribution evolution equation and an ion density distribution evolution equation for calculation to respectively obtain the magnetic flux psi at the next moment (1) Impurity density distribution
Electron density distribution
And ion density distribution
The method comprises the following specific steps:
the evolution equation of the impurity density distribution is:
the evolution equation for the electron density distribution is:
the evolution equation for the ion density distribution is:
wherein n is z Is the particle density distribution of the impurity, n e Representing electron densityCloth, n i Represents an ion density distribution; phi is the potential and t is the time; d z|| And D z⊥ Impurity diffusion coefficients in the parallel direction and the perpendicular direction of the impurity, respectively; d e|| And D e⊥ Impurity diffusion coefficients in the parallel direction and the perpendicular direction of electrons, respectively; d i|| And D i⊥ Impurity diffusion coefficients in the parallel direction and the perpendicular direction of the ions, respectively;
and
gradient operators in the parallel direction and the vertical direction respectively; s radz Is a source term of impurities, S rade Being a source item of electrons, S radi A source term for an ion;
and 5: calculating the impurity density distribution according to the step 4, and when the impurity density distribution of the Tokamak core reaches a set threshold value, further opening a calculation module of plasma radiation, thereby calculating the evolution of the magnetic field configuration along with the time and obtaining the total radiation of the plasma at the moment
Step 6: coupling the calculated plasma total radiation into a magnetofluid equation, calculating the evolution of magnetic flux at a unit time step length, and obtaining the magnetic flux added with the evolution of the plasma total radiation
And 7: the total radiation of the plasma obtained by the calculation in the step 5 is used
Outputting the three-dimensional space distribution information;
and 8: according to the magnetic flux which is obtained after the calculation in the step 6 and takes the total radiation of the plasma into consideration
Further calculating the current magnetic field configuration, and further continuously repeating the steps 5-8 to obtain the total plasma radiation at any moment
And magnetic flux after total radiation of plasma
2. The method as claimed in claim 1, wherein in step 4, the specific calculation method of the evolution equation of the impurity density distribution, the electron density distribution and the ion density distribution is as follows:
(1) The impurity density distribution, electron density distribution and ion density distribution are expressed by a spectrum method:
wherein, (m, n) is the module of the circumferential direction and the polar direction; r 0 Is the large radius of tokamak; theta is a polar angle; z represents the column direction; f is particle distribution, including impurity density distribution, electron density distribution and ion density distribution;
(2) Adopting a two-step prediction-correction method to carry out time advance calculation of various particle distributions; the calculation format of the two-step prediction-correction method is as follows:
and (3) prediction:
and (3) correction:
wherein H represents the magnetic flux and parameters in the magnetic fluid equation; v is the diffusion coefficient; t is time;
represents half a time step; dt represents a time step; the subscript rhs represents the right-hand term of the magnetofluid equation,
representing the gradient in the direction of the perpendicular magnetic field.
3. The method according to claim 1 or 2, wherein the specific steps of the step 5 are as follows:
step 5.1: according to electron temperature distribution T e The corresponding radiation cooling rate L (T) is calculated according to the value of (A) and the generated impurity species e ) The value of (d);
step 5.2: calculating the total radiation of the plasma, including continuous radiation and line radiation;
according to the formula of continuous radiation:
and electron cyclotron radiation formula:
P c ＝5.4×10 -25 n e B 2 T e erg s -1 cm -3
and the calculation formula of the characteristic line radiation:
wherein, P b Representing the magnitude of the continuous radiation, P c Representing the magnitude of electron cyclotron radiation, P z Represents the magnitude of the line radiation; z is a radical of eff Is available electricityDistributing the load;
calculating the magnitude of three kinds of radiation, and then accumulating the three kinds of radiation to obtain the total radiation of the plasma
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202211133612.9A CN115512791A (en) | 2022-09-19 | 2022-09-19 | Particle evolution simulation method based on multi-component plasma radiation effect in tokamak |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202211133612.9A CN115512791A (en) | 2022-09-19 | 2022-09-19 | Particle evolution simulation method based on multi-component plasma radiation effect in tokamak |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN115512791A true CN115512791A (en) | 2022-12-23 |
Family
ID=84504239
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202211133612.9A Pending CN115512791A (en) | 2022-09-19 | 2022-09-19 | Particle evolution simulation method based on multi-component plasma radiation effect in tokamak |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN115512791A (en) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116679863A (en) * | 2023-07-28 | 2023-09-01 | 中国科学院合肥物质科学研究院 | Data processing system suitable for EAST radiation divertor inflation system |
| CN117875218A (en) * | 2024-03-11 | 2024-04-12 | 西南交通大学 | Method and system for inhibiting three-dimensional balance magnetic island in quasi-ring symmetrical star simulator |
-
2022
- 2022-09-19 CN CN202211133612.9A patent/CN115512791A/en active Pending
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116679863A (en) * | 2023-07-28 | 2023-09-01 | 中国科学院合肥物质科学研究院 | Data processing system suitable for EAST radiation divertor inflation system |
| CN116679863B (en) * | 2023-07-28 | 2023-10-24 | 中国科学院合肥物质科学研究院 | Data processing system suitable for EAST radiation divertor inflation system |
| CN117875218A (en) * | 2024-03-11 | 2024-04-12 | 西南交通大学 | Method and system for inhibiting three-dimensional balance magnetic island in quasi-ring symmetrical star simulator |
| CN117875218B (en) * | 2024-03-11 | 2024-05-24 | 西南交通大学 | Method and system for inhibiting three-dimensional balance magnetic island in quasi-ring symmetrical star simulator |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN115512791A (en) | Particle evolution simulation method based on multi-component plasma radiation effect in tokamak | |
| Peysson et al. | Numerical Simulations of the Radio-Frequency–driven Toroidal Current in Tokamaks | |
| Van Devender | Inertial confinement fusion with light-on beams | |
| Gota et al. | Overview of C-2W: high temperature, steady-state beam-driven field-reversed configuration plasmas | |
| Kawamura et al. | First EMC3‐EIRENE Simulations with Divertor Legs of LHD in Realistic Device Geometry | |
| Nakashima et al. | Recent results of divertor simulation experiments using D-module in the GAMMA 10/PDX tandem mirror | |
| Hatayama et al. | Present status of numerical modeling of hydrogen negative ion source plasmas and its comparison with experiments: Japanese activities and their collaboration with experimental groups | |
| Wu et al. | Transport simulation of EAST long-pulse H-mode discharge with integrated modeling | |
| Vainionpaa et al. | High yield neutron generators using the DD reaction | |
| Qian et al. | Nonlinear δf simulation studies of intense ion beam propagation through an alternating-gradient quadrupole focusing field | |
| Kikuchi et al. | Beam instability induced by space charge oscillation during final beam bunching<? format?> for heavy ion inertial fusion | |
| CN114896854B (en) | Simulation method for plasma radiation evolution in tokamak | |
| Kulyk et al. | Modelling of radio-frequency wall conditioning in short pulses in a stellarator | |
| Abe et al. | Analysis of H− extraction in the Linac4 negative ion source by 2.5 D particle simulation | |
| Vshivkova et al. | Numerical modeling of plasma phenomena using the PIC-method | |
| Liu et al. | Seeding the m= 0 instability in dense plasma focus Z-pinches with a hollow anode | |
| Liu et al. | Toroidal modeling of runaway avalanche in DIII-D discharges | |
| Pointon et al. | 3-D PIC simulations of electron flow in the vacuum power flow section of the Z accelerator | |
| Vakili et al. | Conceptual design of 30 kW-NBI injector using a multi-cusp ion source for heating of D-shaped Damavand tokamak plasma | |
| Nishioka et al. | Integrated modeling of the beam formation and extraction in the Linac4 hydrogen negative ion source | |
| Sorokina et al. | Testing of the SOLPS-ITER code at Globus-M2 spherical tokamak with detached divertor | |
| Kobayashi et al. | Experimental study of high density plasma operation in Heliotron J | |
| Johnson | Ionization Cooling | |
| Boozer | Simulation Center for Runaway Electron Avoidance and Mitigation. Final report | |
| Douai et al. | Modelling of Ion Cyclotron Wall Conditioning plasmas |
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 | ||
| CB03 | Change of inventor or designer information | ||
| CB03 | Change of inventor or designer information |
Inventor after: Wang Jincheng Inventor after: Xu Gui Inventor after: Jiang Shuai Inventor after: Wang Zhengxiong Inventor after: Liu Tong Inventor after: Tang Weikang Inventor after: Wei Lai Inventor before: Wang Zhengxiong Inventor before: Jiang Shuai Inventor before: Liu Tong Inventor before: Tang Weikang Inventor before: Wei Lai |