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.展开更多
The turbulence characteristics of plasmas with internal transport barriers in the HL-2A tokamak are analyzed by means of linear gyrokinetic simulations. It is found that turbulence is dominated by the ion temperature ...The turbulence characteristics of plasmas with internal transport barriers in the HL-2A tokamak are analyzed by means of linear gyrokinetic simulations. It is found that turbulence is dominated by the ion temperature gradient(ITG)mode together with large-scale modes characterized by high-frequency electromagnetic fluctuation, which are destabilized by the steep ion temperature gradient in the weak magnetic shear regime. Comparison with solutions of analytical dispersion relations shows that their linear features match well with the beta-induced Alfvén eigenmode branch of the shear Alfvénic spectrum. It is further clarified that the large population of fast ions in these plasmas plays a stabilization role through the dilution mechanism in high-n ITG mode regimes.展开更多
Global linear gyrokinetic simulations using realistic DIII-D tokamak geometry and plasma profiles find co-existence of unstable reversed shear Alfvén eigenmodes(RSAE)with low toroidal mode number n and electromag...Global linear gyrokinetic simulations using realistic DIII-D tokamak geometry and plasma profiles find co-existence of unstable reversed shear Alfvén eigenmodes(RSAE)with low toroidal mode number n and electromagnetic ion temperature gradient(ITG)instabilities with higher toroidal mode number n.For intermediate n?=?[10,12],RSAE and ITG co-exist and overlap weakly in the radial domain with similar growth rates but different real frequencies.Both RSAE and ITG growth rates decrease less than 5%when compressible magnetic perturbations are neglected in the simulations.The ITG growth rates increase less than 7%when fast ions are not included in the simulations.Finally,the effects of trapped electrons on the RSAE are negligible.展开更多
Field-aligned coordinates have been implemented in the gyrokinetic semi-Lagrangian code NLT, Ye et al （2016 J. Comput. Phys. 316 180）, to improve the computational efficiency for the numerical simulations of tokamak...Field-aligned coordinates have been implemented in the gyrokinetic semi-Lagrangian code NLT, Ye et al （2016 J. Comput. Phys. 316 180）, to improve the computational efficiency for the numerical simulations of tokamak turbulence and transport. 4D B-spline interpolation in field- aligned coordinates is applied to solve the gyrokinetic Vlasov equation. A fast iterative algorithm is proposed for efficiently solving the quasi-neutrality equation. A pseudo transform method is used for the numerical integration of the gyro-average operator for perturbations with a high toroidal mode number. The new method is shown to result in an improved code performance for reaching a given accuracy. Some numerical tests are presented to illustrate the new methods.展开更多
We present a new model for simulating the electromagnetic fluctuations with frequencies much lower than the ion cyclotron frequency in plasmas confined in general magnetic configurations.This novel model(termed as GK-...We present a new model for simulating the electromagnetic fluctuations with frequencies much lower than the ion cyclotron frequency in plasmas confined in general magnetic configurations.This novel model(termed as GK-E&B)employs nonlinear gyrokinetic equations formulated in terms of electromagnetic fields along with momentum balance equations for solving fields.It,thus,not only includes kinetic effects,such as wave-particle interaction and microscopic(ion Larmor radius scale)physics;but also is computationally more efficient than the conventional formulation described in terms of potentials.As a benchmark,we perform linear as well as nonlinear simulations of the kinetic Alfvén wave;demonstrating physics in agreement with the analytical theories.展开更多
A method for gyrokinetic simulation of low frequency(lower than the cyclotron frequency)magnetic compressional modes in general geometry is presented.The gyrokinetic-Maxwell system of equations is expressed fully in t...A method for gyrokinetic simulation of low frequency(lower than the cyclotron frequency)magnetic compressional modes in general geometry is presented.The gyrokinetic-Maxwell system of equations is expressed fully in terms of the compressional component of the magnetic perturbation,δBk,with finite Larmor radius effects.This introduces a"gyro-surface"averaging ofδBk in the gyrocenter equations of motion,and similarly in the perpendicular Ampere’s law,which takes the form of the perpendicular force balance equation.The resulting system can be numerically implemented by representing the gyro-surface averaging by a discrete sum in the configuration space.For the typical wavelength of interest(on the order of the gyroradius),the gyro-surface averaging can be reduced to averaging along an effective gyro-orbit.The phase space integration in the force balance equation can be approximated by summing over carefully chosen samples in the magnetic moment coordinate,allowing for an efficient numerical implementation.展开更多
The parameter dependence of transition between electrostatic instabilities is studied using gyrokinetic simulation based on a real discharge of steady-state scenario in the Experimental Advanced Superconducting Tokama...The parameter dependence of transition between electrostatic instabilities is studied using gyrokinetic simulation based on a real discharge of steady-state scenario in the Experimental Advanced Superconducting Tokamak.The scan of radial locations shows that trapped electron mode(TEM)dominates around the core while the ion temperature gradient mode(ITG)simultaneously dominates outside.The maximum growth rate of TEM appears aroundρ=0.24,where the maximum electron temperature gradient R/LTelocates,ρis the normalized poloidal flux.Effects of the parameters on the transition between TEM and ITG instability are studied atρ=0.24.It is found that TEM dominates in the scanning with individually changing R/LTe from 2.50 to 25.02 or the density gradient R/L_(n)from 1.38 to 13.76.Meanwhile,the electron-ion temperature ratio T_(e)/T_(i)is found to destabilize TEM,the effect of Teis more sensitive than that of T_(i).The dominant instability diagrams in the(R/L_(Te),R/L_(Ti))plane at different T_(e)/T_(i)and R/Lnare numerically obtained,which clearly show the parameter range of the dominant TEM or dominant ITG instability region.It is found that the dominant TEM region becomes narrower in the plane by decreasing R/L_(n)when T_(e)/T_(i)>0.5.展开更多
The theoretical and numerical studies on kinetic micro-instabilities,including ion temperature gradient（ITG） driven modes,trapped electron modes（TEMs） in the presence of impurity ions as well as impurity modes（IM...The theoretical and numerical studies on kinetic micro-instabilities,including ion temperature gradient（ITG） driven modes,trapped electron modes（TEMs） in the presence of impurity ions as well as impurity modes（IMs）,induced by impurity density gradient alone,in toroidal magnetized plasmas,such as tokamak and reversed-field pinch（RFP） are reviewed briefly.The basic theory for IMs,the electrostatic instabilities in tokamak and RFP plasmas are discussed.The observations of hybrid and coexistence of the instabilities are categorized systematically.The effects of impurity ions on electromagnetic instabilities such as ITG modes,the kinetic ballooning modes（KBMs） and kinetic shear Alfvén modes induced by impurity ions in tokamak plasmas of finite β（=plasma pressure/magnetic pressure） are analyzed.The interesting topics for future investigation are suggested.展开更多
A 1-Dlinear gyrokinetic code called AWECS is developed to study the kinetic excitation of Alfvenic instabilities in a high-β tokamak plasma,with β being the ratio of thermal to magnetic pressure.It is designed to de...A 1-Dlinear gyrokinetic code called AWECS is developed to study the kinetic excitation of Alfvenic instabilities in a high-β tokamak plasma,with β being the ratio of thermal to magnetic pressure.It is designed to describe physics associated with a broad range of frequencies and wavelengths.For example,AWECS is capable of simulating kinetic ballooning modes,Alfvenic ion-temperature-gradient-driven modes,as well as Alfven instabilities due to energetic particles.In addition,AWECS may be used to study drift-Alfven instabilities in the low-β regime.Here,the layout of the code and the numerical methods used are described.AWECS is benchmarked against other codes and a convergence study is carried out.展开更多
The full torus electromagnetic gyrokinetic particle simulations using the hybrid model with kinetic electrons in the presence of magnetic shear is presented.The fluid-kinetic electron hybrid model employed in this pap...The full torus electromagnetic gyrokinetic particle simulations using the hybrid model with kinetic electrons in the presence of magnetic shear is presented.The fluid-kinetic electron hybrid model employed in this paper improves numerical properties by removing the tearing mode,meanwhile,preserves both linear and nonlinear wave-particle resonances of electrons with Alfven wave and ion acoustic wave.展开更多
The massively parallel,nonlinear,three-dimensional(3D),toroidal,electrostatic,gyrokinetic,particle-in-cell(PIC),Cartesian geometry UCANcode,with particle ions and adiabatic electrons,has been successfully exercised to...The massively parallel,nonlinear,three-dimensional(3D),toroidal,electrostatic,gyrokinetic,particle-in-cell(PIC),Cartesian geometry UCANcode,with particle ions and adiabatic electrons,has been successfully exercised to identify non-diffusive transport characteristics in present day tokamak discharges.The limitation in applying UCAN to larger scale discharges is the 1D domain decomposition in the toroidal(or z-)direction for massively parallel implementation using MPI which has restricted the calculations to a few hundred ion Larmor radii or gyroradii per plasma minor radius.To exceed these sizes,we have implemented 2D domain decomposition in UCANwith the addition of the y-direction to the processor mix.This has been facilitated by use of relevant components in the P2LIB library of field and particle management routines developed for UCLA’s UPIC Framework of conventional PIC codes.The gyroaveraging specific to gyrokinetic codes is simplified by the use of replicated arrays for efficient charge accumulation and force deposition.The 2D domain-decomposed UCAN2 code reproduces the original 1D domain nonlinear results within round-off.Benchmarks of UCAN2 on the Cray XC30 Edison at NERSC demonstrate ideal scaling when problem size is increased along with processor number up to the largest power of 2 available,namely 131,072 processors.These particle weak scaling benchmarks also indicate that the 1 nanosecond per particle per time step and 1 TFlops barriers are easily broken by UCAN2 with 1 billion particles or more and 2000 or more processors.展开更多
A new predictive computer simulation tool targeting the development of the H-mode pedestal at the plasma edge in tokamaks and the triggering and dynamics of edge localizedmodes(ELMs)is presented in this report.This to...A new predictive computer simulation tool targeting the development of the H-mode pedestal at the plasma edge in tokamaks and the triggering and dynamics of edge localizedmodes(ELMs)is presented in this report.This tool brings together,in a coordinated and effective manner,several first-principles physics simulation codes,stability analysis packages,and data processing and visualization tools.A Kepler workflow is used in order to carry out an edge plasma simulation that loosely couples the kinetic code,XGC0,with an ideal MHD linear stability analysis code,ELITE,and an extended MHD initial value code such as M3D or NIMROD.XGC0 includes the neoclassical ion-electron-neutral dynamics needed to simulate pedestal growth near the separatrix.The Kepler workflow processes the XGC0 simulation results into simple images that can be selected and displayed via the Dashboard,a monitoring tool implemented in AJAX allowing the scientist to track computational resources,examine running and archived jobs,and view key physics data,all within a standard Web browser.The XGC0 simulation is monitored for the conditions needed to trigger an ELM crash by periodically assessing the edge plasma pressure and current density profiles using the ELITE code.If an ELM crash is triggered,the Kepler workflow launches the M3D code on a moderate-size Opteron cluster to simulate the nonlinear ELM crash and to compute the relaxation of plasma profiles after the crash.This process is monitored through periodic outputs of plasma fluid quantities that are automatically visualized with AVS/Express and may be displayed on the Dashboard.Finally,the Kepler workflow archives all data outputs and processed images using HPSS,as well as provenance information about the software and hardware used to create the simulation.The complete process of preparing,executing and monitoring a coupled-code simulation of the edge pressure pedestal buildup and the ELM cycle using the Kepler scientific workflow system is described in this paper.展开更多
基金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 Research and Development Program of China (Grant No. 2017YFE0301201)partially by the National Natural Science Foundation of China (Grant Nos. U1967206 and 11775069)。
文摘The turbulence characteristics of plasmas with internal transport barriers in the HL-2A tokamak are analyzed by means of linear gyrokinetic simulations. It is found that turbulence is dominated by the ion temperature gradient(ITG)mode together with large-scale modes characterized by high-frequency electromagnetic fluctuation, which are destabilized by the steep ion temperature gradient in the weak magnetic shear regime. Comparison with solutions of analytical dispersion relations shows that their linear features match well with the beta-induced Alfvén eigenmode branch of the shear Alfvénic spectrum. It is further clarified that the large population of fast ions in these plasmas plays a stabilization role through the dilution mechanism in high-n ITG mode regimes.
基金supported by the China National Magnetic Confinement Fusion Science Program(Grant No.2018YFE0304100)the US Department of Energy,Office of Science,Office of Advanced Scientific Computing Research and Office of Fusion Energy Sciences,and the Scientific Discovery through Advanced Computing(Sci DAC)program under Award No.DE-SC0018270(Sci DAC ISEP Center)+2 种基金the China Scholarship Council(Grant No.201806010067)used resources of the Oak Ridge Leadership Computing Facility at the Oak Ridge National Laboratory(DOE Contract No.DEAC05-00OR22725)the National Energy Research Scientific Computing Center(DOE Contract No.DE-AC02-05CH11231)
文摘Global linear gyrokinetic simulations using realistic DIII-D tokamak geometry and plasma profiles find co-existence of unstable reversed shear Alfvén eigenmodes(RSAE)with low toroidal mode number n and electromagnetic ion temperature gradient(ITG)instabilities with higher toroidal mode number n.For intermediate n?=?[10,12],RSAE and ITG co-exist and overlap weakly in the radial domain with similar growth rates but different real frequencies.Both RSAE and ITG growth rates decrease less than 5%when compressible magnetic perturbations are neglected in the simulations.The ITG growth rates increase less than 7%when fast ions are not included in the simulations.Finally,the effects of trapped electrons on the RSAE are negligible.
基金supported by National Natural Science Foundation of China under Grant Nos.11505240,11375196 and 11405174the National ITER program of China under Contract No.2014GB113000
文摘Field-aligned coordinates have been implemented in the gyrokinetic semi-Lagrangian code NLT, Ye et al （2016 J. Comput. Phys. 316 180）, to improve the computational efficiency for the numerical simulations of tokamak turbulence and transport. 4D B-spline interpolation in field- aligned coordinates is applied to solve the gyrokinetic Vlasov equation. A fast iterative algorithm is proposed for efficiently solving the quasi-neutrality equation. A pseudo transform method is used for the numerical integration of the gyro-average operator for perturbations with a high toroidal mode number. The new method is shown to result in an improved code performance for reaching a given accuracy. Some numerical tests are presented to illustrate the new methods.
基金supported by the National Science Foundation of China(Grant Nos.11235009,and 11905097)funding from the Euratom Research And Training Programme 2014-2018 and 2019-2020(Grant Agreement No.633053,Project No.WP19-ER/ENEA-05)。
文摘We present a new model for simulating the electromagnetic fluctuations with frequencies much lower than the ion cyclotron frequency in plasmas confined in general magnetic configurations.This novel model(termed as GK-E&B)employs nonlinear gyrokinetic equations formulated in terms of electromagnetic fields along with momentum balance equations for solving fields.It,thus,not only includes kinetic effects,such as wave-particle interaction and microscopic(ion Larmor radius scale)physics;but also is computationally more efficient than the conventional formulation described in terms of potentials.As a benchmark,we perform linear as well as nonlinear simulations of the kinetic Alfvén wave;demonstrating physics in agreement with the analytical theories.
基金supported by National Science Foundation and U.S.Department of Energy grants.
文摘A method for gyrokinetic simulation of low frequency(lower than the cyclotron frequency)magnetic compressional modes in general geometry is presented.The gyrokinetic-Maxwell system of equations is expressed fully in terms of the compressional component of the magnetic perturbation,δBk,with finite Larmor radius effects.This introduces a"gyro-surface"averaging ofδBk in the gyrocenter equations of motion,and similarly in the perpendicular Ampere’s law,which takes the form of the perpendicular force balance equation.The resulting system can be numerically implemented by representing the gyro-surface averaging by a discrete sum in the configuration space.For the typical wavelength of interest(on the order of the gyroradius),the gyro-surface averaging can be reduced to averaging along an effective gyro-orbit.The phase space integration in the force balance equation can be approximated by summing over carefully chosen samples in the magnetic moment coordinate,allowing for an efficient numerical implementation.
基金supported by the National MCF Energy R&D Program of China(Nos.2019YFE03060000,2019YFE03050000 and 2019YFE03020004)National Natural Science Foundation of China(Nos.12005063 and 11875131)+1 种基金Users with Excellence Program of Hefei Science Center CAS(Nos.2020HSC-UE011 and 2021HSC-UE015)Anhui Provincial Natural Science Foundation(No.2008085Jo4)。
文摘The parameter dependence of transition between electrostatic instabilities is studied using gyrokinetic simulation based on a real discharge of steady-state scenario in the Experimental Advanced Superconducting Tokamak.The scan of radial locations shows that trapped electron mode(TEM)dominates around the core while the ion temperature gradient mode(ITG)simultaneously dominates outside.The maximum growth rate of TEM appears aroundρ=0.24,where the maximum electron temperature gradient R/LTelocates,ρis the normalized poloidal flux.Effects of the parameters on the transition between TEM and ITG instability are studied atρ=0.24.It is found that TEM dominates in the scanning with individually changing R/LTe from 2.50 to 25.02 or the density gradient R/L_(n)from 1.38 to 13.76.Meanwhile,the electron-ion temperature ratio T_(e)/T_(i)is found to destabilize TEM,the effect of Teis more sensitive than that of T_(i).The dominant instability diagrams in the(R/L_(Te),R/L_(Ti))plane at different T_(e)/T_(i)and R/Lnare numerically obtained,which clearly show the parameter range of the dominant TEM or dominant ITG instability region.It is found that the dominant TEM region becomes narrower in the plane by decreasing R/L_(n)when T_(e)/T_(i)>0.5.
基金supported by National Natural Science Foundation of China(Nos.11475057 and 11575158)the National Key R&D Program of China under Grant No.2017YFE0300405
文摘The theoretical and numerical studies on kinetic micro-instabilities,including ion temperature gradient（ITG） driven modes,trapped electron modes（TEMs） in the presence of impurity ions as well as impurity modes（IMs）,induced by impurity density gradient alone,in toroidal magnetized plasmas,such as tokamak and reversed-field pinch（RFP） are reviewed briefly.The basic theory for IMs,the electrostatic instabilities in tokamak and RFP plasmas are discussed.The observations of hybrid and coexistence of the instabilities are categorized systematically.The effects of impurity ions on electromagnetic instabilities such as ITG modes,the kinetic ballooning modes（KBMs） and kinetic shear Alfvén modes induced by impurity ions in tokamak plasmas of finite β（=plasma pressure/magnetic pressure） are analyzed.The interesting topics for future investigation are suggested.
基金This research is supported by U.S.DoE Grant DE-AC02-CH0-3073,NSF Grant ATM-0335279,and in part by SciDAC GSEP.
文摘A 1-Dlinear gyrokinetic code called AWECS is developed to study the kinetic excitation of Alfvenic instabilities in a high-β tokamak plasma,with β being the ratio of thermal to magnetic pressure.It is designed to describe physics associated with a broad range of frequencies and wavelengths.For example,AWECS is capable of simulating kinetic ballooning modes,Alfvenic ion-temperature-gradient-driven modes,as well as Alfven instabilities due to energetic particles.In addition,AWECS may be used to study drift-Alfven instabilities in the low-β regime.Here,the layout of the code and the numerical methods used are described.AWECS is benchmarked against other codes and a convergence study is carried out.
基金This work is supported by Department of Energy(DOE)Cooperative Agreement No.DE-FC02-04ER54796Grant No.DE-FC02-06ER54860in part by SciDAC GPS,GSEP,and CPES centers。
文摘The full torus electromagnetic gyrokinetic particle simulations using the hybrid model with kinetic electrons in the presence of magnetic shear is presented.The fluid-kinetic electron hybrid model employed in this paper improves numerical properties by removing the tearing mode,meanwhile,preserves both linear and nonlinear wave-particle resonances of electrons with Alfven wave and ion acoustic wave.
基金This work was supported in part in the USA by Grant No.DE-FG02-04ER54741 to the University of Alaska,Fairbanks,AK,from the Office of Fusion Energy Sciences,Office of Science,United States Department of EnergyIt was also supported in part at Univer-sidad Carlos III,Madrid,Spain,by Spanish National Project No.ENE2009-12213-C03-03+1 种基金This research used resources of the National Energy Research Scientific Computing Center(NERSC)which is supported by the Office of Science of the U.S.Department of Energy under Contract No.DE-AC02-05CH11231.
文摘The massively parallel,nonlinear,three-dimensional(3D),toroidal,electrostatic,gyrokinetic,particle-in-cell(PIC),Cartesian geometry UCANcode,with particle ions and adiabatic electrons,has been successfully exercised to identify non-diffusive transport characteristics in present day tokamak discharges.The limitation in applying UCAN to larger scale discharges is the 1D domain decomposition in the toroidal(or z-)direction for massively parallel implementation using MPI which has restricted the calculations to a few hundred ion Larmor radii or gyroradii per plasma minor radius.To exceed these sizes,we have implemented 2D domain decomposition in UCANwith the addition of the y-direction to the processor mix.This has been facilitated by use of relevant components in the P2LIB library of field and particle management routines developed for UCLA’s UPIC Framework of conventional PIC codes.The gyroaveraging specific to gyrokinetic codes is simplified by the use of replicated arrays for efficient charge accumulation and force deposition.The 2D domain-decomposed UCAN2 code reproduces the original 1D domain nonlinear results within round-off.Benchmarks of UCAN2 on the Cray XC30 Edison at NERSC demonstrate ideal scaling when problem size is increased along with processor number up to the largest power of 2 available,namely 131,072 processors.These particle weak scaling benchmarks also indicate that the 1 nanosecond per particle per time step and 1 TFlops barriers are easily broken by UCAN2 with 1 billion particles or more and 2000 or more processors.
基金This work is part of the ongoing research activities within the SciDAC Fusion Simulation Prototype(FSP)Center for Plasma Edge Simulationwhich is supported by the Office of Fusion Energy Sciences and the Office of Advanced Scientific Computing Research within the U.S.Department of Energy.
文摘A new predictive computer simulation tool targeting the development of the H-mode pedestal at the plasma edge in tokamaks and the triggering and dynamics of edge localizedmodes(ELMs)is presented in this report.This tool brings together,in a coordinated and effective manner,several first-principles physics simulation codes,stability analysis packages,and data processing and visualization tools.A Kepler workflow is used in order to carry out an edge plasma simulation that loosely couples the kinetic code,XGC0,with an ideal MHD linear stability analysis code,ELITE,and an extended MHD initial value code such as M3D or NIMROD.XGC0 includes the neoclassical ion-electron-neutral dynamics needed to simulate pedestal growth near the separatrix.The Kepler workflow processes the XGC0 simulation results into simple images that can be selected and displayed via the Dashboard,a monitoring tool implemented in AJAX allowing the scientist to track computational resources,examine running and archived jobs,and view key physics data,all within a standard Web browser.The XGC0 simulation is monitored for the conditions needed to trigger an ELM crash by periodically assessing the edge plasma pressure and current density profiles using the ELITE code.If an ELM crash is triggered,the Kepler workflow launches the M3D code on a moderate-size Opteron cluster to simulate the nonlinear ELM crash and to compute the relaxation of plasma profiles after the crash.This process is monitored through periodic outputs of plasma fluid quantities that are automatically visualized with AVS/Express and may be displayed on the Dashboard.Finally,the Kepler workflow archives all data outputs and processed images using HPSS,as well as provenance information about the software and hardware used to create the simulation.The complete process of preparing,executing and monitoring a coupled-code simulation of the edge pressure pedestal buildup and the ELM cycle using the Kepler scientific workflow system is described in this paper.
