Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userDevelopment of CFETR scenarios with self-consistent core-pedestal coupled simulations
59
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
This paper develops two non-inductive steady state scenarios for larger size configuration of China Fusion Engineering Test Reactor (CFETR) with integrated modeling simulations. A self-consistent core-pedestal coupled workflow for CFETR is developed under integrated modeling framework OMFIT, which allows more accurate evaluation of CFETR performance. The workflow integrates equilibrium code EFIT, transport codes ONETWO and TGYRO, and pedestal code EPED. A fully non-inductive baseline phase I scenario is developed with the workflow, which satisfies the minimum goal of Fusion Nuclear Science Facility. Compared with previous work, which proves the larger size and higher toroidal field CFETR configuration than has the advantages of reducing heating and current drive requirements, lowering divertor and wall power loads, allowing higher bootstrap current fraction and better confinement. A fully non-inductive high-performance phase II scenario is developed, which explores the alpha-particle dominated self-heating regime. Phase II scenario achieves the target of fusion power Pfus>1GW and fusion gain Qfus>20, and it largely reduces auxiliary heating and current drive power. Moreover, the large neutron production of phase II increases the energy generation power and tritium breeding rate.
rate research
Read More
We present the first successful simulation of a induced disruption in ASDEX Upgrade from massive material injection (MMI) up to established runaway electron (RE) beam, thus covering pre-thermal quench, thermal quench and current quench (CQ) of the discharge. For future high-current fusion devices such as ITER, the successful suppression of REs through MMI is of critical importance to ensure the structural integrity of the vessel. To computationally study the interplay between MMI, background plasma response, and RE generation, a toolkit based on the 1.5D transport code coupling ASTRA-STRAHL is developed. Electron runaway is described by state-of-the-art reduced kinetic models in the presence of partially ionized impurities. Applied to argon MMI in ASDEX Upgrade discharge #33108, key plasma parameters measured experimentally, such as temporal evolution of the line averaged electron density, plasma current decay rate and post-CQ RE current, are well reproduced by the simulation presented. Impurity ions are transported into the central plasma by the combined effect of neoclassical processes and additional effects prescribed inside the $q = 2$ rational surface to explain experimental time scales. Thus, a thermal collapse is induced through strong impurity radiation, giving rise to a substantial RE population as observed experimentally.
The Particle-in-Cell (PIC) method was used to study two different ion thruster concepts - Stationary Plasma Thrusters (SPT) and High Efficiency Multistage Plasma Thrusters (HEMP-T), in particular the plasma properties in the discharge chamber due to the different magnetic field configurations. Special attention was paid to the simulation of plasma particle fluxes on the thrusters channel surfaces. In both cases, PIC proved itself as a powerful tool, delivering important insight into the basic physics of the different thruster concepts. The simulations demonstrated that the new HEMP thruster concept allows for a high thermal efficiency due to both minimal energy dissipation and high acceleration efficiency. In the HEMP thruster the plasma contact to the wall is limited only to very small areas of the magnetic field cusps, which results in much smaller ion energy flux to the thruster channel surface as compared to SPT. The erosion yields for dielectric discharge channel walls of SPT and HEMP thrusters were calculated with the binary collision code SDTrimSP. For SPT, an erosion rate on the level of 1 mm of sputtered material per hour was observed. For HEMP, thruster simulations have shown that there is no erosion inside the dielectric discharge channel.
A new fluid model for runaway electron simulation based on fluid description is introduced and implemented in the magnetohydrodynamics code M3D-C1, which includes self-consistent interactions between plasma and runaway electrons. The model utilizes the method of characteristics to solve the continuity equation for the runaway electron density with large convection speed, and uses a modified Boris algorithm for pseudo particle pushing. The model was employed to simulate magnetohydrodynamics instabilities happening in a runaway electron final loss event in the DIII-D tokamak. Nonlinear simulation reveals that a large fraction of runaway electrons get lost to the wall when kink instabilities are excited and form stochastic field lines in the outer region of the plasma. Plasma current converts from runaway electron current to Ohmic current, and get pinched at the magnetic axis. Given the good agreement with experiment, the simulation model provides a reliable tool to study macroscopic plasma instabilities in existence of runaway electron current, and can be used to support future studies of runaway electron mitigation strategies in ITER.
This paper discusses the capture of an ion beam in a magnetized plasma of an Electron Cyclotron Resonance Ion Source based Charge Breeder, as modelled by numerical simulations. As a relevant step forward with respect to previous works, here the capture is modeled by considering a plasma structure determined in a self-consisent way. The plasmoid-halo structure of the ECR plasma - that is consisting of a dense core (the plasmoid) surrounded by a rarefied halo - is further confirmed by the self-consistent calculations, having also some fine structures affected by the electromagnetic field distribution and by the magnetostatic field profile. The capture of Rb1+ ions has been investigated in details, vs. various plasma parameters, and then compared to experimental results.
The understanding of the halo current properties during disruptions is key to design and operate large scale tokamaks in view of the large thermal and electromagnetic loads that they entail. For the first time, we present a fully self-consistent model for halo current simulations including neutral particles and sheath boundary conditions. The model is used to simulate Vertical Displacement Events (VDEs) occurring in the COMPASS tokamak. Recent COMPASS experiments have shown that the parallel halo current density at the plasma-wall interface is limited by the ion saturation current during VDE-induced disruptions. We show that usual MHD boundary conditions can lead to the violation of this physical limit and we implement this current density limitation through a boundary condition for the electrostatic potential. Sheath boundary conditions for the density, the heat flux, the parallel velocity and a realistic parameter choice (e.g. Spitzer $eta$ and Spitzer-Harm $chi_parallel$ values) extend present VDE simulations beyond the state of the art. Experimental measurements of the current density, temperature and heat flux profiles at the COMPASS divertor are compared with the results obtained from axisymmetric simulations. Since the ion saturation current density ($J_{sat}$) is shown to be essential to determine the halo current profile, parametric scans are performed to study its dependence on different quantities such as the plasma resistivity and the particle and heat diffusion coefficients. In this respect, the plasma resistivity in the halo region broadens significantly the $J_{sat}$ profile, increasing the halo width at a similar total halo current.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
