Aluminum is the primary structural material in nuclear engineering,and its cross section induced by 14-MeV neutrons is of great significance.To address the issue of insufficient accuracy for the^(27)Al(n,2n)^(26)Al re...Aluminum is the primary structural material in nuclear engineering,and its cross section induced by 14-MeV neutrons is of great significance.To address the issue of insufficient accuracy for the^(27)Al(n,2n)^(26)Al reaction cross section,the activation method and accelerator mass spectrometry(AMS)technique were used to determine the^(27)Al(n,2n)^(26)Al cross section,which could be used as a D-T plasma ion temperature monitor in fusion reactors.At the China Academy of Engineering Physics,neutron activation was performed using a K-400 neutron generator produced by the T(d,n)4He reaction.The^(26)Al∕^(27)Al isotope ratios were measured using the newly installed GYIG 1 MV AMS at the Institute of Geochemistry,Chinese Academy of Sciences.The neutron flux was monitored by measuring the activity of 92mNb produced by the 93Nb(n,2n)92mNb reaction.The measured results were compared with available data in the experimental nuclear reaction database,and the measured values showed a reasonable degree of consistency with partially available literature data.The newly acquired cross-sectional data at 12 neutron energy points through systematic measurements clarified the divergence,which has two different growth trends from the existing experimental values.The obtained results are also compared with the corresponding evaluated database,and the newly calculated excitation functions with TALYS−1.95 and EMPIRE−3.2 codes,the agreement with CENDL−3.2,TENDL-2021 and EMPIRE−3.2 results are generally acceptable.A substantial improvement in the knowledge of the^(27)Al(n,2n)^(26)Al reaction excitation function was obtained in the present work,which will lay the foundation for the diagnosis of the fusion ion temperature,testing of the nuclear physics model,evaluation of nuclear data,etc.展开更多
Many magnetohydrodynamic stability analyses require generation of a set of equilibria with a fixed safety factor q-profile while varying other plasma parameters.A neural network(NN)-based approach is investigated that...Many magnetohydrodynamic stability analyses require generation of a set of equilibria with a fixed safety factor q-profile while varying other plasma parameters.A neural network(NN)-based approach is investigated that facilitates such a process.Both multilayer perceptron(MLP)-based NN and convolutional neural network(CNN)models are trained to map the q-profile to the plasma current density J-profile,and vice versa,while satisfying the Grad–Shafranov radial force balance constraint.When the initial target models are trained,using a database of semianalytically constructed numerical equilibria,an initial CNN with one convolutional layer is found to perform better than an initial MLP model.In particular,a trained initial CNN model can also predict the q-or J-profile for experimental tokamak equilibria.The performance of both initial target models is further improved by fine-tuning the training database,i.e.by adding realistic experimental equilibria with Gaussian noise.The fine-tuned target models,referred to as fine-tuned MLP and fine-tuned CNN,well reproduce the target q-or J-profile across multiple tokamak devices.As an important application,these NN-based equilibrium profile convertors can be utilized to provide a good initial guess for iterative equilibrium solvers,where the desired input quantity is the safety factor instead of the plasma current density.展开更多
Measurements of the total radiated power and its spatial distribution are crucial for fusion research.On the experimental advanced superconducting tokamak(EAST),both the metal foil resistive bolometer and the absolute...Measurements of the total radiated power and its spatial distribution are crucial for fusion research.On the experimental advanced superconducting tokamak(EAST),both the metal foil resistive bolometer and the absolute extreme ultraviolet(AXUV)photodiodes have been used to quantify the radiated power.This article introduces the latest improvement of the bolometer diagnostic system on EAST.It also details the successful design and installation of new divertor AXUV cameras,which are dedicated to the investigation of divertor physics.The shielding components of the bolometer detector have been refined,and the article provides a detailed exposition of the double shielding structures that have been verified as effective in microwave shielding.Additionally,the changes in the radiated power distribution in the divertor region during the plasma detachment process are measured using the divertor AXUV camera.Finally,the radiated power measured by the AXUV detector and metal foil resistive bolometer are compared,and different detector performances are presented.展开更多
In a tokamak fusion reactor operated at steady state,the equilibrium magnetic field is likely to have reversed shear in the core region,as the noninductive bootstrap current profile generally peaks off-axis.The revers...In a tokamak fusion reactor operated at steady state,the equilibrium magnetic field is likely to have reversed shear in the core region,as the noninductive bootstrap current profile generally peaks off-axis.The reversed shear Alfvén eigenmode(RSAE)as a unique branch of the shear Alfvén wave in this equilibrium,can exist with a broad spectrum in wavenumber and frequency,and be resonantly driven unstable by energetic particles(EP).After briefly discussing the RSAE linear properties in burning plasma condition,we review several key topics of the nonlinear dynamics for the RSAE through both wave-EP resonance and wave-wave coupling channels,and illustrate their potentially important role in reactor-scale fusion plasmas.By means of simplified hybrid MHD-kinetic simulations,the RSAEs are shown to have typically broad phase space resonance structure with both circulating and trapped EP,as results of weak/vanishing magnetic shear and relatively low frequency.Through the route of wave-EP nonlinearity,the dominant saturation mechanism is mainly due to the transported resonant EP radially decoupling with the localized RSAE mode structure,and the resultant EP transport generally has a convective feature.The saturated RSAEs also undergo various nonlinear couplings with other collective oscillations.Two typical routes as parametric decay and modulational instability are studied using nonlinear gyrokinetic theory,and applied to the scenario of spontaneous excitation by a finite amplitude pump RSAE.Multiple RSAEs could naturally couple and induce the spectral energy cascade into a low frequency Alfvénic mode,which may effectively transfer the EP energy to fuel ions via collisionless Landau damping.Moreover,zero frequency zonal field structure could be spontaneously excited by modulation of the pump RSAE envelope,and may also lead to saturation of the pump RSAE by both scattering into stable domain and local distortion of the continuum structure.展开更多
The hybrid scenario,which has good confinement and moderate MHD instabilities,is a proposed operation scenario for international thermonuclear experimental reactor(ITER).In this work,the effect of plasma rotation on t...The hybrid scenario,which has good confinement and moderate MHD instabilities,is a proposed operation scenario for international thermonuclear experimental reactor(ITER).In this work,the effect of plasma rotation on the HL-3 hybrid scenario is analyzed with the integrated modeling framework OMFIT.The results show that toroidal rotation has no obvious effect on confinement with a high line averaged density of n_(bar)~(7)×10^(19)m^(-3).In this case,the ion temperature only changes from 4.7 keV to 4.4 keV with the rotation decreasing from 10^(5) rad/s to 10^(3) rad/s,which means that the turbulent heat transport is not dominant.While in the scenarios characterized by lower densities,such as n_(bar)~4×10^(19)m^(-3),turbulent transport becomes dominant in determining heat transport.The ion temperature rises from 3.8 keV to 6.1 keV in the core as the rotation velocity increases from 10^(3) rad/s to 10^(5) rad/s.Despite the ion temperature rising,the rotation velocity does not obviously affect electron temperature or density.Additionally,it is noteworthy that the variation in rotation velocity does not significantly affect the global confinement of plasma in scenarios with low density or with high density.展开更多
An accurate plasma current profile has irreplaceable value for the steady-state operation of the plasma.In this study,plasma current tomography based on Bayesian inference is applied to an HL-2A device and used to rec...An accurate plasma current profile has irreplaceable value for the steady-state operation of the plasma.In this study,plasma current tomography based on Bayesian inference is applied to an HL-2A device and used to reconstruct the plasma current profile.Two different Bayesian probability priors are tried,namely the Conditional Auto Regressive(CAR)prior and the Advanced Squared Exponential(ASE)kernel prior.Compared to the CAR prior,the ASE kernel prior adopts nonstationary hyperparameters and introduces the current profile of the reference discharge into the hyperparameters,which can make the shape of the current profile more flexible in space.The results indicate that the ASE prior couples more information,reduces the probability of unreasonable solutions,and achieves higher reconstruction accuracy.展开更多
The linear and nonlinear simulations are carried out using the gyrokinetic code NLT for the electrostatic instabilities in the core region of a deuterium plasma based on the International Thermonuclear Experimental Re...The linear and nonlinear simulations are carried out using the gyrokinetic code NLT for the electrostatic instabilities in the core region of a deuterium plasma based on the International Thermonuclear Experimental Reactor(ITER)baseline scenario.The kinetic electron effects on the linear frequency and nonlinear transport are studied by adopting the adiabatic electron model and the fully drift-kinetic electron model in the NLT code,respectively.The linear simulations focus on the dependence of linear frequency on the plasma parameters,such as the ion and electron temperature gradientsκ_(Ti,e)≡R=L_(Ti,e),the density gradientκ_(n)≡R/L_(n)and the ion-electron temperature ratioτ=T_(e)=T_(i).Here,is the major radius,and T_(e)and T_(i)denote the electron and ion temperatures,respectively.L_(A)=-(δ_(r)lnA)^(-1)is the gradient scale length,with denoting the density,the ion and electron temperatures,respectively.In the kinetic electron model,the ion temperature gradient(ITG)instability and the trapped electron mode(TEM)dominate in the small and large k_(θ)region,respectively,wherek_(θ)is the poloidal wavenumber.The TEMdominant region becomes wider by increasing(decreasing)κ_(T_(e))(κ_(T_(i)))or by decreasingκ_(n).For the nominal parameters of the ITER baseline scenario,the maximum growth rate of dominant ITG instability in the kinetic electron model is about three times larger than that in the adiabatic electron model.The normalized linear frequency depends on the value ofτ,rather than the value of T_(e)or T_(i),in both the adiabatic and kinetic electron models.The nonlinear simulation results show that the ion heat diffusivity in the kinetic electron model is quite a lot larger than that in the adiabatic electron model,the radial structure is finer and the time oscillation is more rapid.In addition,the magnitude of the fluctuated potential at the saturated stage peaks in the ITGdominated region,and contributions from the TEM(dominating in the higher k_(θ)region)to the nonlinear transport can be neglected.In the adiabatic electron model,the zonal radial electric field is found to be mainly driven by the turbulent energy flux,and the contribution of turbulent poloidal Reynolds stress is quite small due to the toroidal shielding effect.However,in the kinetic electron model,the turbulent energy flux is not strong enough to drive the zonal radial electric field in the nonlinear saturated stage.The kinetic electron effects on the mechanism of the turbulence-driven zonal radial electric field should be further investigated.展开更多
First mirror(FM)cleaning,using radio frequency(RF)plasma,has been proposed to recover FM reflectivity in nuclear fusion reactors such as the International Thermonuclear Experimental Reactor(ITER).To investigate the in...First mirror(FM)cleaning,using radio frequency(RF)plasma,has been proposed to recover FM reflectivity in nuclear fusion reactors such as the International Thermonuclear Experimental Reactor(ITER).To investigate the influence of simultaneous cleaning of two mirrors on mirror cleaning efficiency and uniformity,experiments involving single-mirror cleaning and dual-mirror cleaning were conducted using RF capacitively coupled plasma in the laboratory.For the test and simultaneous cleaning of two mirrors,the FM and second mirror(SM),both measuring 110 mm×80 mm,were placed inside the first mirror unit(FMU).They were composed of 16 mirror samples,each with a dimension of 27.5 mm×20 mm.These mirror samples consist of a titanium-zirconium-molybdenum alloy substrate,a 500 nm Mo intermediate layer and a 30 nm Al_(2)O_(3) surface coating as a proxy for Be impurities.The cleaning of a single first mirror(SFM)and the simultaneous cleaning of the FM and SM(DFM and DSM)lasted for 9 h using Ar plasma at a pressure of 1 Pa.The total reflectivity of mirror samples on the DSM did not fully recover and varied with location,with a self-bias of−140 V.With a self-bias of−300 V,the total reflectivity of mirror samples on the SFM and DFM was fully recovered.The energy dispersive spectrometer results demonstrated that the Al_(2)O_(3) coating had been completely removed from these mirror samples.However,the mass loss of each mirror sample on the SFM and DFM before and after cleaning varied depending on its location,with higher mass loss observed for mirror samples located in the corners and lower loss for those in the center.Compared with SM cleaning,the simultaneous cleaning of two mirrors reduced the difference between the highest and lowest mass loss.Furthermore,this mass loss for the mirror samples of the DFM facing the DSM was increased.This indicated that mirror samples cleaned face to face in the FMU simultaneously could influence each other,highlighting the need for special attention in future studies.展开更多
Tokamak plasmas with elongated cross sections are susceptible to vertical displacement events(VDEs),which can damage the first wall via heat flux or electromagnetic(EM)forces.We present a 3D nonlinear reduced magnetoh...Tokamak plasmas with elongated cross sections are susceptible to vertical displacement events(VDEs),which can damage the first wall via heat flux or electromagnetic(EM)forces.We present a 3D nonlinear reduced magnetohydrodynamic(MHD)simulation of CFETR plasma during a cold VDE following the thermal quench,focusing on the relationship among the EM force,plasma displacement,and the n=1 mode.The dominant mode,identified as m/n=2/1,becomes destabilized when most of the current is contracted within the q=2 surface.The displacement of the plasma current centroid is less than that of the magnetic axis due to the presence of SOL current in the open field line region.Hence,the symmetric component of the induced vacuum vessel current is significantly mitigated.The direction of the sideways force keeps a constant phase approximately compared to the asymmetric component of the vacuum vessel current and the SOL current,which in turn keep in-phase with the dominant 2/1 mode.Their amplitudes are also closely associated with the growth of the dominant mode.These findings provide insights into potential methods for controlling the phase and amplitude of sideways forces during VDEs in the future.展开更多
Field reversed configuration(FRC)is widely considered as an ideal target plasma for magnetoinertial fusion.However,its confinement and stability,both proportional to the radius,will deteriorate inevitably during radia...Field reversed configuration(FRC)is widely considered as an ideal target plasma for magnetoinertial fusion.However,its confinement and stability,both proportional to the radius,will deteriorate inevitably during radial compression.Hence,we propose a new fusion approach based on axial compression of a large-sized FRC.The axial compression can be made by plasma jets or plasmoids converging onto the axial ends of the FRC.The parameter space that can reach the ignition condition while preserving the FRC's overall quality is studied using a numerical model based on different FRC confinement scalings.It is found that ignition is possible for a large FRC that can be achieved with the current FRC formation techniques if compression ratio is greater than 50.A more realistic compression is to combine axial with moderate radial compression,which is also presented and calculated in this work.展开更多
A three-fluid equilibrium plasma with bulk plasma and energetic electrons has been observed on the Xuanlong-50(EXL-50) spherical torus, where the energetic electrons play a crucial role in sustaining the plasma curren...A three-fluid equilibrium plasma with bulk plasma and energetic electrons has been observed on the Xuanlong-50(EXL-50) spherical torus, where the energetic electrons play a crucial role in sustaining the plasma current and pressure. In this study, the equilibrium of a multi-fluid plasma was investigated by analyzing the relationship between the external vertical magnetic field(B_(V)),plasma current(I_(p)), the poloidal ratio(β_(p)) and the Shafranov formula. Remarkably, our research demonstrates some validity of the Shafranov formula in the presence of multi-fluid plasma in EXL-50 spherical torus. This finding holds significant importance for future reactors as it allows for differentiation between alpha particles and background plasma. The study of multi-fluid plasma provides a significant reference value for the equilibrium reconstruction of burning plasma involving alpha particles.展开更多
离子回旋射频(Ion Cyclotron Range of Frequencies,ICRF)波加热是托卡马克装置上至关重要的辅助加热方式之一。托卡马克装置中国环流三号(HL-3,原名HL-2M)拟安装加热功率为6 MW的ICRF加热系统。本工作利用TRANSP程序,模拟并研究了ICRF...离子回旋射频(Ion Cyclotron Range of Frequencies,ICRF)波加热是托卡马克装置上至关重要的辅助加热方式之一。托卡马克装置中国环流三号(HL-3,原名HL-2M)拟安装加热功率为6 MW的ICRF加热系统。本工作利用TRANSP程序,模拟并研究了ICRF加热的频率和功率对聚变中子产额以及快离子分布的影响。研究结果表明:ICRF的频率和功率对中子产额有显著影响,固定ICRF频率时,中子产额与加热功率成正比关系,而在固定ICRF加热功率的情况下,中子产额的增加幅度显著依赖ICRF的频率,在研究参数范围内,30 MHz的ICRF对中子产额的增加具有最显著的增强作用。快离子分布的模拟结果显示,在考虑ICRF加热后,中性束和ICRF的协同加热机制能够将快离子加热至最高1 MeV,有效地提高了中子产额。此外,基于中子相机诊断的概念对中子信号进行了仿真。结果表明,中子相机能够有效地测量到由ICRF加热导致的中子产额高低和分布剖面的变化,这为将来优化中子相机诊断系统设计和测量中子空间分布提供了一定的参考。展开更多
基金the Open Project of Guangxi Key Laboratory of Nuclear Physics and Nuclear Technology(NLK 2022-04)the Central Government Guidance Funds for Local Scientific and Technological Development,China(No.Guike,ZY22096024)+1 种基金the National Natural Science Foundation of China(12065003)Guangxi Key R&D Project(2023AB07029).
文摘Aluminum is the primary structural material in nuclear engineering,and its cross section induced by 14-MeV neutrons is of great significance.To address the issue of insufficient accuracy for the^(27)Al(n,2n)^(26)Al reaction cross section,the activation method and accelerator mass spectrometry(AMS)technique were used to determine the^(27)Al(n,2n)^(26)Al cross section,which could be used as a D-T plasma ion temperature monitor in fusion reactors.At the China Academy of Engineering Physics,neutron activation was performed using a K-400 neutron generator produced by the T(d,n)4He reaction.The^(26)Al∕^(27)Al isotope ratios were measured using the newly installed GYIG 1 MV AMS at the Institute of Geochemistry,Chinese Academy of Sciences.The neutron flux was monitored by measuring the activity of 92mNb produced by the 93Nb(n,2n)92mNb reaction.The measured results were compared with available data in the experimental nuclear reaction database,and the measured values showed a reasonable degree of consistency with partially available literature data.The newly acquired cross-sectional data at 12 neutron energy points through systematic measurements clarified the divergence,which has two different growth trends from the existing experimental values.The obtained results are also compared with the corresponding evaluated database,and the newly calculated excitation functions with TALYS−1.95 and EMPIRE−3.2 codes,the agreement with CENDL−3.2,TENDL-2021 and EMPIRE−3.2 results are generally acceptable.A substantial improvement in the knowledge of the^(27)Al(n,2n)^(26)Al reaction excitation function was obtained in the present work,which will lay the foundation for the diagnosis of the fusion ion temperature,testing of the nuclear physics model,evaluation of nuclear data,etc.
基金supported by National Natural Science Foundation of China (Nos. 12205033, 12105317, 11905022 and 11975062)Dalian Youth Science and Technology Project (No. 2022RQ039)+1 种基金the Fundamental Research Funds for the Central Universities (No. 3132023192)the Young Scientists Fund of the Natural Science Foundation of Sichuan Province (No. 2023NSFSC1291)
文摘Many magnetohydrodynamic stability analyses require generation of a set of equilibria with a fixed safety factor q-profile while varying other plasma parameters.A neural network(NN)-based approach is investigated that facilitates such a process.Both multilayer perceptron(MLP)-based NN and convolutional neural network(CNN)models are trained to map the q-profile to the plasma current density J-profile,and vice versa,while satisfying the Grad–Shafranov radial force balance constraint.When the initial target models are trained,using a database of semianalytically constructed numerical equilibria,an initial CNN with one convolutional layer is found to perform better than an initial MLP model.In particular,a trained initial CNN model can also predict the q-or J-profile for experimental tokamak equilibria.The performance of both initial target models is further improved by fine-tuning the training database,i.e.by adding realistic experimental equilibria with Gaussian noise.The fine-tuned target models,referred to as fine-tuned MLP and fine-tuned CNN,well reproduce the target q-or J-profile across multiple tokamak devices.As an important application,these NN-based equilibrium profile convertors can be utilized to provide a good initial guess for iterative equilibrium solvers,where the desired input quantity is the safety factor instead of the plasma current density.
基金supported by the National MCF Energy R&D Program of China(Nos.2022YFE03040003 and 2019YFE03010003).
文摘Measurements of the total radiated power and its spatial distribution are crucial for fusion research.On the experimental advanced superconducting tokamak(EAST),both the metal foil resistive bolometer and the absolute extreme ultraviolet(AXUV)photodiodes have been used to quantify the radiated power.This article introduces the latest improvement of the bolometer diagnostic system on EAST.It also details the successful design and installation of new divertor AXUV cameras,which are dedicated to the investigation of divertor physics.The shielding components of the bolometer detector have been refined,and the article provides a detailed exposition of the double shielding structures that have been verified as effective in microwave shielding.Additionally,the changes in the radiated power distribution in the divertor region during the plasma detachment process are measured using the divertor AXUV camera.Finally,the radiated power measured by the AXUV detector and metal foil resistive bolometer are compared,and different detector performances are presented.
基金supported by National Natural Science Foundation of China (Nos. 12205251, 12275236 and 12261131622)Italian Ministry for Foreign Affairs and International Cooperation Project (No. CN23GR02)+2 种基金the National Key Research and Development Program of China (Nos. 2019YFE03020003 and 2017YFE0301900)Users of Excellence program of Hefei Science Center CAS (No. 2021HSC-UE016)funded by the European Union via the Euratom Research and Training Programme (No. 101052200–EUROfusion)
文摘In a tokamak fusion reactor operated at steady state,the equilibrium magnetic field is likely to have reversed shear in the core region,as the noninductive bootstrap current profile generally peaks off-axis.The reversed shear Alfvén eigenmode(RSAE)as a unique branch of the shear Alfvén wave in this equilibrium,can exist with a broad spectrum in wavenumber and frequency,and be resonantly driven unstable by energetic particles(EP).After briefly discussing the RSAE linear properties in burning plasma condition,we review several key topics of the nonlinear dynamics for the RSAE through both wave-EP resonance and wave-wave coupling channels,and illustrate their potentially important role in reactor-scale fusion plasmas.By means of simplified hybrid MHD-kinetic simulations,the RSAEs are shown to have typically broad phase space resonance structure with both circulating and trapped EP,as results of weak/vanishing magnetic shear and relatively low frequency.Through the route of wave-EP nonlinearity,the dominant saturation mechanism is mainly due to the transported resonant EP radially decoupling with the localized RSAE mode structure,and the resultant EP transport generally has a convective feature.The saturated RSAEs also undergo various nonlinear couplings with other collective oscillations.Two typical routes as parametric decay and modulational instability are studied using nonlinear gyrokinetic theory,and applied to the scenario of spontaneous excitation by a finite amplitude pump RSAE.Multiple RSAEs could naturally couple and induce the spectral energy cascade into a low frequency Alfvénic mode,which may effectively transfer the EP energy to fuel ions via collisionless Landau damping.Moreover,zero frequency zonal field structure could be spontaneously excited by modulation of the pump RSAE envelope,and may also lead to saturation of the pump RSAE by both scattering into stable domain and local distortion of the continuum structure.
基金Project supported by the National Magnetic Confinement Fusion Program of China (Grants Nos.2019YFE03040002 and 2018YFE0301101)the Talent Project of China National Nuclear Corporation,China (Grant No.2022JZYF-01)。
文摘The hybrid scenario,which has good confinement and moderate MHD instabilities,is a proposed operation scenario for international thermonuclear experimental reactor(ITER).In this work,the effect of plasma rotation on the HL-3 hybrid scenario is analyzed with the integrated modeling framework OMFIT.The results show that toroidal rotation has no obvious effect on confinement with a high line averaged density of n_(bar)~(7)×10^(19)m^(-3).In this case,the ion temperature only changes from 4.7 keV to 4.4 keV with the rotation decreasing from 10^(5) rad/s to 10^(3) rad/s,which means that the turbulent heat transport is not dominant.While in the scenarios characterized by lower densities,such as n_(bar)~4×10^(19)m^(-3),turbulent transport becomes dominant in determining heat transport.The ion temperature rises from 3.8 keV to 6.1 keV in the core as the rotation velocity increases from 10^(3) rad/s to 10^(5) rad/s.Despite the ion temperature rising,the rotation velocity does not obviously affect electron temperature or density.Additionally,it is noteworthy that the variation in rotation velocity does not significantly affect the global confinement of plasma in scenarios with low density or with high density.
基金supported by the National MCF Energy R&D Program of China (Nos. 2018 YFE0301105, 2022YFE03010002 and 2018YFE0302100)the National Key R&D Program of China (Nos. 2022YFE03070004 and 2022YFE03070000)National Natural Science Foundation of China (Nos. 12205195, 12075155 and 11975277)
文摘An accurate plasma current profile has irreplaceable value for the steady-state operation of the plasma.In this study,plasma current tomography based on Bayesian inference is applied to an HL-2A device and used to reconstruct the plasma current profile.Two different Bayesian probability priors are tried,namely the Conditional Auto Regressive(CAR)prior and the Advanced Squared Exponential(ASE)kernel prior.Compared to the CAR prior,the ASE kernel prior adopts nonstationary hyperparameters and introduces the current profile of the reference discharge into the hyperparameters,which can make the shape of the current profile more flexible in space.The results indicate that the ASE prior couples more information,reduces the probability of unreasonable solutions,and achieves higher reconstruction accuracy.
基金supported by the National MCF Energy R&D Program of China(No.2019YFE03060000)National Natural Science Foundation of China(Nos.12005063,12375215 and 12175034)the Collaborative Innovation Program of Hefei Science Center,CAS(No.2022HSC-CIP008).
文摘The linear and nonlinear simulations are carried out using the gyrokinetic code NLT for the electrostatic instabilities in the core region of a deuterium plasma based on the International Thermonuclear Experimental Reactor(ITER)baseline scenario.The kinetic electron effects on the linear frequency and nonlinear transport are studied by adopting the adiabatic electron model and the fully drift-kinetic electron model in the NLT code,respectively.The linear simulations focus on the dependence of linear frequency on the plasma parameters,such as the ion and electron temperature gradientsκ_(Ti,e)≡R=L_(Ti,e),the density gradientκ_(n)≡R/L_(n)and the ion-electron temperature ratioτ=T_(e)=T_(i).Here,is the major radius,and T_(e)and T_(i)denote the electron and ion temperatures,respectively.L_(A)=-(δ_(r)lnA)^(-1)is the gradient scale length,with denoting the density,the ion and electron temperatures,respectively.In the kinetic electron model,the ion temperature gradient(ITG)instability and the trapped electron mode(TEM)dominate in the small and large k_(θ)region,respectively,wherek_(θ)is the poloidal wavenumber.The TEMdominant region becomes wider by increasing(decreasing)κ_(T_(e))(κ_(T_(i)))or by decreasingκ_(n).For the nominal parameters of the ITER baseline scenario,the maximum growth rate of dominant ITG instability in the kinetic electron model is about three times larger than that in the adiabatic electron model.The normalized linear frequency depends on the value ofτ,rather than the value of T_(e)or T_(i),in both the adiabatic and kinetic electron models.The nonlinear simulation results show that the ion heat diffusivity in the kinetic electron model is quite a lot larger than that in the adiabatic electron model,the radial structure is finer and the time oscillation is more rapid.In addition,the magnitude of the fluctuated potential at the saturated stage peaks in the ITGdominated region,and contributions from the TEM(dominating in the higher k_(θ)region)to the nonlinear transport can be neglected.In the adiabatic electron model,the zonal radial electric field is found to be mainly driven by the turbulent energy flux,and the contribution of turbulent poloidal Reynolds stress is quite small due to the toroidal shielding effect.However,in the kinetic electron model,the turbulent energy flux is not strong enough to drive the zonal radial electric field in the nonlinear saturated stage.The kinetic electron effects on the mechanism of the turbulence-driven zonal radial electric field should be further investigated.
基金supported by the National Key R&D Project of China(No.2022YFE03030000)National Natural Science Foundation of China(Nos.11975269,12275306 and 12075279)+3 种基金the Youth Innovation Promotion Association of the Chinese Academy of Sciences(No.2022452)the Anhui Provincial Natural Science Foundation(No.2208085J40)the CASHIPS Director’s Fund(Nos.YZJJQY202302 and BJPY2023B03)the Comprehensive Research Facility for Fusion Technology Program of China(No.2018-000052-73-01-001228).
文摘First mirror(FM)cleaning,using radio frequency(RF)plasma,has been proposed to recover FM reflectivity in nuclear fusion reactors such as the International Thermonuclear Experimental Reactor(ITER).To investigate the influence of simultaneous cleaning of two mirrors on mirror cleaning efficiency and uniformity,experiments involving single-mirror cleaning and dual-mirror cleaning were conducted using RF capacitively coupled plasma in the laboratory.For the test and simultaneous cleaning of two mirrors,the FM and second mirror(SM),both measuring 110 mm×80 mm,were placed inside the first mirror unit(FMU).They were composed of 16 mirror samples,each with a dimension of 27.5 mm×20 mm.These mirror samples consist of a titanium-zirconium-molybdenum alloy substrate,a 500 nm Mo intermediate layer and a 30 nm Al_(2)O_(3) surface coating as a proxy for Be impurities.The cleaning of a single first mirror(SFM)and the simultaneous cleaning of the FM and SM(DFM and DSM)lasted for 9 h using Ar plasma at a pressure of 1 Pa.The total reflectivity of mirror samples on the DSM did not fully recover and varied with location,with a self-bias of−140 V.With a self-bias of−300 V,the total reflectivity of mirror samples on the SFM and DFM was fully recovered.The energy dispersive spectrometer results demonstrated that the Al_(2)O_(3) coating had been completely removed from these mirror samples.However,the mass loss of each mirror sample on the SFM and DFM before and after cleaning varied depending on its location,with higher mass loss observed for mirror samples located in the corners and lower loss for those in the center.Compared with SM cleaning,the simultaneous cleaning of two mirrors reduced the difference between the highest and lowest mass loss.Furthermore,this mass loss for the mirror samples of the DFM facing the DSM was increased.This indicated that mirror samples cleaned face to face in the FMU simultaneously could influence each other,highlighting the need for special attention in future studies.
基金supported by the National MCF Energy R&D Program of China(Grant Nos.2019YFE03010001 and 2018YFE0311300).
文摘Tokamak plasmas with elongated cross sections are susceptible to vertical displacement events(VDEs),which can damage the first wall via heat flux or electromagnetic(EM)forces.We present a 3D nonlinear reduced magnetohydrodynamic(MHD)simulation of CFETR plasma during a cold VDE following the thermal quench,focusing on the relationship among the EM force,plasma displacement,and the n=1 mode.The dominant mode,identified as m/n=2/1,becomes destabilized when most of the current is contracted within the q=2 surface.The displacement of the plasma current centroid is less than that of the magnetic axis due to the presence of SOL current in the open field line region.Hence,the symmetric component of the induced vacuum vessel current is significantly mitigated.The direction of the sideways force keeps a constant phase approximately compared to the asymmetric component of the vacuum vessel current and the SOL current,which in turn keep in-phase with the dominant 2/1 mode.Their amplitudes are also closely associated with the growth of the dominant mode.These findings provide insights into potential methods for controlling the phase and amplitude of sideways forces during VDEs in the future.
基金supported by National Natural Science Foundation of China(No.12175226)。
文摘Field reversed configuration(FRC)is widely considered as an ideal target plasma for magnetoinertial fusion.However,its confinement and stability,both proportional to the radius,will deteriorate inevitably during radial compression.Hence,we propose a new fusion approach based on axial compression of a large-sized FRC.The axial compression can be made by plasma jets or plasmoids converging onto the axial ends of the FRC.The parameter space that can reach the ignition condition while preserving the FRC's overall quality is studied using a numerical model based on different FRC confinement scalings.It is found that ignition is possible for a large FRC that can be achieved with the current FRC formation techniques if compression ratio is greater than 50.A more realistic compression is to combine axial with moderate radial compression,which is also presented and calculated in this work.
文摘A three-fluid equilibrium plasma with bulk plasma and energetic electrons has been observed on the Xuanlong-50(EXL-50) spherical torus, where the energetic electrons play a crucial role in sustaining the plasma current and pressure. In this study, the equilibrium of a multi-fluid plasma was investigated by analyzing the relationship between the external vertical magnetic field(B_(V)),plasma current(I_(p)), the poloidal ratio(β_(p)) and the Shafranov formula. Remarkably, our research demonstrates some validity of the Shafranov formula in the presence of multi-fluid plasma in EXL-50 spherical torus. This finding holds significant importance for future reactors as it allows for differentiation between alpha particles and background plasma. The study of multi-fluid plasma provides a significant reference value for the equilibrium reconstruction of burning plasma involving alpha particles.
文摘离子回旋射频(Ion Cyclotron Range of Frequencies,ICRF)波加热是托卡马克装置上至关重要的辅助加热方式之一。托卡马克装置中国环流三号(HL-3,原名HL-2M)拟安装加热功率为6 MW的ICRF加热系统。本工作利用TRANSP程序,模拟并研究了ICRF加热的频率和功率对聚变中子产额以及快离子分布的影响。研究结果表明:ICRF的频率和功率对中子产额有显著影响,固定ICRF频率时,中子产额与加热功率成正比关系,而在固定ICRF加热功率的情况下,中子产额的增加幅度显著依赖ICRF的频率,在研究参数范围内,30 MHz的ICRF对中子产额的增加具有最显著的增强作用。快离子分布的模拟结果显示,在考虑ICRF加热后,中性束和ICRF的协同加热机制能够将快离子加热至最高1 MeV,有效地提高了中子产额。此外,基于中子相机诊断的概念对中子信号进行了仿真。结果表明,中子相机能够有效地测量到由ICRF加热导致的中子产额高低和分布剖面的变化,这为将来优化中子相机诊断系统设计和测量中子空间分布提供了一定的参考。