Reconstruction of plasma equilibrium plays an important role in the analysis and simulation of plasma experiments. The kinetic equilibrium reconstruction with pressure and edge current constraints has been employed on...Reconstruction of plasma equilibrium plays an important role in the analysis and simulation of plasma experiments. The kinetic equilibrium reconstruction with pressure and edge current constraints has been employed on EAST tokamak. However, the internal safety factor(q) profile is not accurate. This paper proposes a new way of incorporating q profile constraints into kinetic equilibrium reconstruction. The q profile is yielded from the Polarimeter Interferometer(POINT)reconstruction. Virtual probes containing information on q profile constraints are added to inputs of the kinetic equilibrium reconstruction program to obtain the final equilibrium. The new equilibrium produces a more accurate internal q profile. This improved method would help analyze EAST experiments.展开更多
Plasma equilibrium reconstruction provides essential information for tokamak operation and physical analysis.An extensive and reliable set of magnetic diagnostics is required to obtain accurate plasma equilibrium.This...Plasma equilibrium reconstruction provides essential information for tokamak operation and physical analysis.An extensive and reliable set of magnetic diagnostics is required to obtain accurate plasma equilibrium.This study designs and optimizes the magnetic diagnostics layout for the reconstruction of the equilibrium of the plasma according to the scientific objectives,engineering design parameters,and limitations of the Chinese Fusion Engineering Test Reactor(CFETR).Based on the CFETR discharge simulation,magnetic measurement data are employed to reconstruct consistent plasma equilibrium parameters,and magnetic diagnostics'number and position are optimized by truncated Singular value decomposition,verifying the redundancy reliability of the magnetic diagnostics layout design.This provides a design solution for the layout of the magnetic diagnostics system required to control the plasma equilibrium of CFETR,and the developed design and optimization method can provide effective support to design magnetic diagnostics systems for future magnetic confinement fusion devices.展开更多
Reconstruction of experimental axisymmetric equilibria is an important part of tokamak data analysis. Fourier expansion is applied to reconstruct the vessel current distribution in EFIT code. Benchmarking and testing ...Reconstruction of experimental axisymmetric equilibria is an important part of tokamak data analysis. Fourier expansion is applied to reconstruct the vessel current distribution in EFIT code. Benchmarking and testing calculations are performed to evaluate and validate this algorithm. Two cases for circular and non-circular plasma discharges are presented. Fourier expansion used to fit the eddy current is a robust method and the real time EFIT can be introduced to the plasma control system in the coming campaign.展开更多
To achieve real-time control of tokamak plasmas, the equilibrium reconstruction has to be completed sufficiently quickly. For the case of an EAST tokamak experiment, real-time equilibrium reconstruction is generally r...To achieve real-time control of tokamak plasmas, the equilibrium reconstruction has to be completed sufficiently quickly. For the case of an EAST tokamak experiment, real-time equilibrium reconstruction is generally required to provide results within 1ms. A graphic processing unit(GPU) parallel Grad–Shafranov(G-S) solver is developed in P-EFIT code,which is built with the CUDA? architecture to take advantage of massively parallel GPU cores and significantly accelerate the computation. Optimization and implementation of numerical algorithms for a block tri-diagonal linear system are presented. The solver can complete a calculation within 16 μs with 65×65 grid size and 27 μs with 129×129 grid size, and this solver supports that P-EFIT can fulfill the time feasibility for real-time plasma control with both grid sizes.展开更多
The EFIT program is integrated with the high resolution laser polarimeter interferometer system(POLARIS), the soft X-ray imaging diagnostic system(SXR) and the electron cyclotron emission radiometer(ECE) in the ...The EFIT program is integrated with the high resolution laser polarimeter interferometer system(POLARIS), the soft X-ray imaging diagnostic system(SXR) and the electron cyclotron emission radiometer(ECE) in the J-TEXT tokamak. Then some internal information about Faraday angle and the position of safety factor q=1 can be obtained as a constraint to EFIT. The modified EFIT code is used to calculate the internal parameters such as flux function, safety factor q, pressure and current density.展开更多
The equilibrium reconstruction is important to study the tokamak plasma physical processes.To analyze the contribution of fast ions to the equilibrium,the kinetic equilibria at two time-slices in a typical H-mode disc...The equilibrium reconstruction is important to study the tokamak plasma physical processes.To analyze the contribution of fast ions to the equilibrium,the kinetic equilibria at two time-slices in a typical H-mode discharge with different auxiliary heatings are reconstructed by using magnetic diagnostics,kinetic diagnostics and TRANSP code.It is found that the fast-ion pressure might be up to one-third of the plasma pressure and the contribution is mainly in the core plasma due to the neutral beam injection power is primarily deposited in the core region.The fast-ion current contributes mainly in the core region while contributes little to the pedestal current.A steep pressure gradient in the pedestal is observed which gives rise to a strong edge current.It is proved that the fast ion effects cannot be ignored and should be considered in the future study of EAST.展开更多
The equilibrium fitting code (EFIT) and its application in the J-TEXT tokamak are integrated by the Matlab language. The function of analysis and visualization to the results is added. In addition, the experiment da...The equilibrium fitting code (EFIT) and its application in the J-TEXT tokamak are integrated by the Matlab language. The function of analysis and visualization to the results is added. In addition, the experiment data measured by soft X-ray (SXR) are used to calculate plasma equilibrium as a constraint condition. The improved EFIT code is used for J-TEXT discharge and the profiles of plasma parameters such as flux function, safety factor q, pressure and current density are obtained from the reconstructed configurations.展开更多
Real-time equilibrium reconstruction is crucially important for plasma shape control in the process of tokamak plasma discharge.However,as the reconstruction algorithm is computationally intensive,it is very difficult...Real-time equilibrium reconstruction is crucially important for plasma shape control in the process of tokamak plasma discharge.However,as the reconstruction algorithm is computationally intensive,it is very difficult to improve its accuracy and reduce the computation time,and some optimizations need to be done.This article describes the three most important aspects of this optimization:(1) compiler optimization;(2) some optimization for middle-scale matrix multiplication on the graphic processing unit and an algorithm which can solve the block tri-diagonal linear system efficiently in parallel;(3) a new algorithm to locate the X&O point on the central processing unit.A static test proves the correctness and a dynamic test proves the feasibility of using the new code for real-time reconstruction with 129?×?129 grids;it can complete one iteration around 575 μs for each equilibrium reconstruction.The plasma displacements from real-time equilibrium reconstruction are compared with the experimental measurements,and the calculated results are consistent with the measured ones,which can be used as a reference for the real-time control of HL-2 A discharge.展开更多
In the reversed field pinch(RFP),plasmas exhibit various self-organized states.Among these,the three-dimensional(3D)helical state known as the“quasi-single-helical”(QSH)state enhances RFP confinement.However,accurat...In the reversed field pinch(RFP),plasmas exhibit various self-organized states.Among these,the three-dimensional(3D)helical state known as the“quasi-single-helical”(QSH)state enhances RFP confinement.However,accurately describing the equilibrium is challenging due to the presence of 3D structures,magnetic islands,and chaotic regions.It is difficult to obtain a balance between the available diagnostic and the real equilibrium structure.To address this issue,we introduce KTX3DFit,a new 3D equilibrium reconstruction code specifically designed for the Keda Torus eXperiment(KTX)RFP.KTX3DFit utilizes the stepped-pressure equilibrium code(SPEC)to compute 3D equilibria and uses polarimetric interferometer signals from experiments.KTX3DFit is able to reconstruct equilibria in various states,including axisymmetric,doubleaxis helical(DAx),and single-helical-axis(SHAx)states.Notably,this study marks the first integration of the SPEC code with internal magnetic field data for equilibrium reconstruction and could be used for other 3D configurations.展开更多
The EFIT (Equilibrium Fitting) code is modified for the equilibrium configuration reconstruction in HL-2A. Signals from Langmuir probe (LP) at the divertor target plates are employed in the reconstruction of diver...The EFIT (Equilibrium Fitting) code is modified for the equilibrium configuration reconstruction in HL-2A. Signals from Langmuir probe (LP) at the divertor target plates are employed in the reconstruction of divertor configurations. The results show that discharge #2895 starts with a limiter configuration and develops gradually into a divertor configuration after t = 230 ms. This transition process is clearly demonstrated by the LP signals for the reconstruction. The profiles of plasma parameters such as safety factor q, pressure and current density as well as the evolution of major shape parameters of plasma, such as the boundary magnetic fluxes, the positions of both x-point and magnetic axis, are calculated from the reconstructed configurations. The possibility to apply the method to the swing of strike point on the target plate is discussed.展开更多
An online plasma shape reconstruction, based on the offiine version of the EFIT code and MPI library, can be carried out between two adjacent shots in EAST. It combines online data acquisition, parallel calculation, a...An online plasma shape reconstruction, based on the offiine version of the EFIT code and MPI library, can be carried out between two adjacent shots in EAST. It combines online data acquisition, parallel calculation, and data storage together. The program on the master node of the cluster detects the termination of the discharge promptly, reads diagnostic data from the EAST mdsplus server on the completion of data storing, and writes the results onto the EFIT mdsplus server after the calculation is finished. These processes run automatically on a nine-nodes IBM blade center. The total time elapsed is about 1 second to several minutes, depending on the duration of the shot. With the results stored in the mdsplus server, it is convenient for operators and physicists to analyze the behavior of plasma using visualization tools.展开更多
In recent EAST experiments,current profile broadening characterized by reduced internal inductance has been achieved by utilizing radio-frequency current drives(RFCD).In contrast to previous density scan experiments,w...In recent EAST experiments,current profile broadening characterized by reduced internal inductance has been achieved by utilizing radio-frequency current drives(RFCD).In contrast to previous density scan experiments,which showed an outward shift of the current density profile of lower hybrid current drive(LHCD)in higher plasma density,the core electron temperature(T_(e)(0))is found to affect the LHCD current profile as well.According to equilibrium reconstruction,a significant increase in on-axis safety factor(q0)from 2.05 to 3.41 is observed by careful arrangement of RFCD.Simulations using ray-tracing code GENRAY and Fokker–Planck code CQL3D have been performed to thoroughly analyze the LHCD current profile,revealing the sensitivity of the LHCD current profile to T_(e)(0).The LHCD current density tends to accumulate in the plasma core with higher current drive efficiency benefiting from higher T_(e)(0).With a lower T_(e)(0),the LHCD current profile broadens due to off-axis deposition of power density.The sensitivity of the power deposition and current profile of LHCD to T_(e)(0)provides a promising way to effectively optimize current profile via control of the core electron temperature.展开更多
Plasmas with vertically elongated cross-sections tend to be unstable to an axis-symmetric instability. This paper studies the magnetohydrodynamic equilibria in elongated plasmas after failure of vertical feedback cont...Plasmas with vertically elongated cross-sections tend to be unstable to an axis-symmetric instability. This paper studies the magnetohydrodynamic equilibria in elongated plasmas after failure of vertical feedback control by using magnetic data for EAST device. Vertical forces on the vessel due to the induced polodial and toroidal currents are evaluated. The maximum force of the Fzpol in vertical displacement events for EAST designed parameters is given.展开更多
基金supported by National Key R&D Program of China(Nos.2019YFE03040004 and 2017YFE0300404)supported by Comprehensive Research Facility for Fusion Technology Program of China(No.2018000052-73-01-001228)。
文摘Reconstruction of plasma equilibrium plays an important role in the analysis and simulation of plasma experiments. The kinetic equilibrium reconstruction with pressure and edge current constraints has been employed on EAST tokamak. However, the internal safety factor(q) profile is not accurate. This paper proposes a new way of incorporating q profile constraints into kinetic equilibrium reconstruction. The q profile is yielded from the Polarimeter Interferometer(POINT)reconstruction. Virtual probes containing information on q profile constraints are added to inputs of the kinetic equilibrium reconstruction program to obtain the final equilibrium. The new equilibrium produces a more accurate internal q profile. This improved method would help analyze EAST experiments.
基金Project supported by the National MCF Energy Research and Development Program of China (Grant Nos.2022YFE03010002,2018YFE0302100,and 2018YFE0301105)the National Natural Science Foundation of China (Grant Nos.11875291,11805236,11905256,and 12075285)。
文摘Plasma equilibrium reconstruction provides essential information for tokamak operation and physical analysis.An extensive and reliable set of magnetic diagnostics is required to obtain accurate plasma equilibrium.This study designs and optimizes the magnetic diagnostics layout for the reconstruction of the equilibrium of the plasma according to the scientific objectives,engineering design parameters,and limitations of the Chinese Fusion Engineering Test Reactor(CFETR).Based on the CFETR discharge simulation,magnetic measurement data are employed to reconstruct consistent plasma equilibrium parameters,and magnetic diagnostics'number and position are optimized by truncated Singular value decomposition,verifying the redundancy reliability of the magnetic diagnostics layout design.This provides a design solution for the layout of the magnetic diagnostics system required to control the plasma equilibrium of CFETR,and the developed design and optimization method can provide effective support to design magnetic diagnostics systems for future magnetic confinement fusion devices.
基金supported by National Natural Science Foundation of China (No. 10725523)
文摘Reconstruction of experimental axisymmetric equilibria is an important part of tokamak data analysis. Fourier expansion is applied to reconstruct the vessel current distribution in EFIT code. Benchmarking and testing calculations are performed to evaluate and validate this algorithm. Two cases for circular and non-circular plasma discharges are presented. Fourier expansion used to fit the eddy current is a robust method and the real time EFIT can be introduced to the plasma control system in the coming campaign.
基金supported by the National Magnetic Confinement Fusion Research Program of China(Grant No.2014GB103000)the National Natural Science Foundation of China(Grant No.11575245)the National Natural Science Foundation of China for Youth(Grant No.11205191)
文摘To achieve real-time control of tokamak plasmas, the equilibrium reconstruction has to be completed sufficiently quickly. For the case of an EAST tokamak experiment, real-time equilibrium reconstruction is generally required to provide results within 1ms. A graphic processing unit(GPU) parallel Grad–Shafranov(G-S) solver is developed in P-EFIT code,which is built with the CUDA? architecture to take advantage of massively parallel GPU cores and significantly accelerate the computation. Optimization and implementation of numerical algorithms for a block tri-diagonal linear system are presented. The solver can complete a calculation within 16 μs with 65×65 grid size and 27 μs with 129×129 grid size, and this solver supports that P-EFIT can fulfill the time feasibility for real-time plasma control with both grid sizes.
文摘The EFIT program is integrated with the high resolution laser polarimeter interferometer system(POLARIS), the soft X-ray imaging diagnostic system(SXR) and the electron cyclotron emission radiometer(ECE) in the J-TEXT tokamak. Then some internal information about Faraday angle and the position of safety factor q=1 can be obtained as a constraint to EFIT. The modified EFIT code is used to calculate the internal parameters such as flux function, safety factor q, pressure and current density.
基金supported by National Key R&D Program of China under Grant No.2017YFE0300400National Natural Science Foundation of China under Grant Nos.11475220,11405218,11575248+1 种基金the National Magnetic Confinement Fusion Science Program of China under Contracts No.2014GB106001sponsored in part by Youth Innovation Promotion Association Chinese Academy of Sciences (Grant No.2016384)
文摘The equilibrium reconstruction is important to study the tokamak plasma physical processes.To analyze the contribution of fast ions to the equilibrium,the kinetic equilibria at two time-slices in a typical H-mode discharge with different auxiliary heatings are reconstructed by using magnetic diagnostics,kinetic diagnostics and TRANSP code.It is found that the fast-ion pressure might be up to one-third of the plasma pressure and the contribution is mainly in the core plasma due to the neutral beam injection power is primarily deposited in the core region.The fast-ion current contributes mainly in the core region while contributes little to the pedestal current.A steep pressure gradient in the pedestal is observed which gives rise to a strong edge current.It is proved that the fast ion effects cannot be ignored and should be considered in the future study of EAST.
文摘The equilibrium fitting code (EFIT) and its application in the J-TEXT tokamak are integrated by the Matlab language. The function of analysis and visualization to the results is added. In addition, the experiment data measured by soft X-ray (SXR) are used to calculate plasma equilibrium as a constraint condition. The improved EFIT code is used for J-TEXT discharge and the profiles of plasma parameters such as flux function, safety factor q, pressure and current density are obtained from the reconstructed configurations.
基金supported in part by National Natural Science Foundation of China under Grant Nos.11505050 and 11775071the ITER-CN under Grant No.2014GB107003
文摘Real-time equilibrium reconstruction is crucially important for plasma shape control in the process of tokamak plasma discharge.However,as the reconstruction algorithm is computationally intensive,it is very difficult to improve its accuracy and reduce the computation time,and some optimizations need to be done.This article describes the three most important aspects of this optimization:(1) compiler optimization;(2) some optimization for middle-scale matrix multiplication on the graphic processing unit and an algorithm which can solve the block tri-diagonal linear system efficiently in parallel;(3) a new algorithm to locate the X&O point on the central processing unit.A static test proves the correctness and a dynamic test proves the feasibility of using the new code for real-time reconstruction with 129?×?129 grids;it can complete one iteration around 575 μs for each equilibrium reconstruction.The plasma displacements from real-time equilibrium reconstruction are compared with the experimental measurements,and the calculated results are consistent with the measured ones,which can be used as a reference for the real-time control of HL-2 A discharge.
基金supported by National Natural Science Foundation of China(Nos.12175227 and 12375226)the National Magnetic Confinement Fusion Program of China(No.2022YFE03100004)+1 种基金the Fundamental Research Funds for the Central Universities(No.USTC 20210079)the Collaborative Innovation Program of Hefei Science Center,CAS(No.2022HSC-CIP022)。
文摘In the reversed field pinch(RFP),plasmas exhibit various self-organized states.Among these,the three-dimensional(3D)helical state known as the“quasi-single-helical”(QSH)state enhances RFP confinement.However,accurately describing the equilibrium is challenging due to the presence of 3D structures,magnetic islands,and chaotic regions.It is difficult to obtain a balance between the available diagnostic and the real equilibrium structure.To address this issue,we introduce KTX3DFit,a new 3D equilibrium reconstruction code specifically designed for the Keda Torus eXperiment(KTX)RFP.KTX3DFit utilizes the stepped-pressure equilibrium code(SPEC)to compute 3D equilibria and uses polarimetric interferometer signals from experiments.KTX3DFit is able to reconstruct equilibria in various states,including axisymmetric,doubleaxis helical(DAx),and single-helical-axis(SHAx)states.Notably,this study marks the first integration of the SPEC code with internal magnetic field data for equilibrium reconstruction and could be used for other 3D configurations.
基金supported by National Natural Science Foundation of China(No.10635010)Specialized Research Fund for the Doctoral Program of Higher Education of China(No.20090041110026)
文摘The EFIT (Equilibrium Fitting) code is modified for the equilibrium configuration reconstruction in HL-2A. Signals from Langmuir probe (LP) at the divertor target plates are employed in the reconstruction of divertor configurations. The results show that discharge #2895 starts with a limiter configuration and develops gradually into a divertor configuration after t = 230 ms. This transition process is clearly demonstrated by the LP signals for the reconstruction. The profiles of plasma parameters such as safety factor q, pressure and current density as well as the evolution of major shape parameters of plasma, such as the boundary magnetic fluxes, the positions of both x-point and magnetic axis, are calculated from the reconstructed configurations. The possibility to apply the method to the swing of strike point on the target plate is discussed.
基金supported by National Natural Science Foundation of China (No.10835009)the Knowledge Innovation Program of The Chinese Academy of Sciences (No.KJCX3.SYW.N4)
文摘An online plasma shape reconstruction, based on the offiine version of the EFIT code and MPI library, can be carried out between two adjacent shots in EAST. It combines online data acquisition, parallel calculation, and data storage together. The program on the master node of the cluster detects the termination of the discharge promptly, reads diagnostic data from the EAST mdsplus server on the completion of data storing, and writes the results onto the EFIT mdsplus server after the calculation is finished. These processes run automatically on a nine-nodes IBM blade center. The total time elapsed is about 1 second to several minutes, depending on the duration of the shot. With the results stored in the mdsplus server, it is convenient for operators and physicists to analyze the behavior of plasma using visualization tools.
基金supported by the National MCF Energy R&D Program of China (No. 2019YFE0304000)National Natural Science Foundation of China (Nos. 12005262 and 11975274)+3 种基金the Anhui Provincial Natural Science Foundation (No. 2108085J06)the Users with Excellence Program of Hefei Science Center CAS (Nos. 2021HSC-UE018 and 2020HSC-UE011)the External Cooperation Program of Chinese Academy of Sciences (No. 116134KYSB20180035)the Science Foundation of Institute of Plasma Physics, Chinese Academy of Sciences (No. DSJJ-2021-04)
文摘In recent EAST experiments,current profile broadening characterized by reduced internal inductance has been achieved by utilizing radio-frequency current drives(RFCD).In contrast to previous density scan experiments,which showed an outward shift of the current density profile of lower hybrid current drive(LHCD)in higher plasma density,the core electron temperature(T_(e)(0))is found to affect the LHCD current profile as well.According to equilibrium reconstruction,a significant increase in on-axis safety factor(q0)from 2.05 to 3.41 is observed by careful arrangement of RFCD.Simulations using ray-tracing code GENRAY and Fokker–Planck code CQL3D have been performed to thoroughly analyze the LHCD current profile,revealing the sensitivity of the LHCD current profile to T_(e)(0).The LHCD current density tends to accumulate in the plasma core with higher current drive efficiency benefiting from higher T_(e)(0).With a lower T_(e)(0),the LHCD current profile broadens due to off-axis deposition of power density.The sensitivity of the power deposition and current profile of LHCD to T_(e)(0)provides a promising way to effectively optimize current profile via control of the core electron temperature.
基金supported by the National Natural Science Foundation of China (Grant No 10725523)the Ministry of Science and Technology of China (Grant No 2007DFA11290)
文摘Plasmas with vertically elongated cross-sections tend to be unstable to an axis-symmetric instability. This paper studies the magnetohydrodynamic equilibria in elongated plasmas after failure of vertical feedback control by using magnetic data for EAST device. Vertical forces on the vessel due to the induced polodial and toroidal currents are evaluated. The maximum force of the Fzpol in vertical displacement events for EAST designed parameters is given.