Search

SIDON (SImulator of raDiO-frequency Networks) is an in-house developed Radio-Frequency (RF) network solver that has been implemented to cross-validate the design of WEST ICRF launchers and simulate their impedance matching algorithm while considering all mutual couplings and asymmetries. In this paper, the authors illustrate the theory of SIDON as well as results of its calculations. The authors have built time-varying plasma scenarios (a sequence of launchers front-faces L-mode and H-mode Z-matrices), where at each time step (1 millisecond here), SIDON solves the RF network. At the same time, when activated, the impedance matching algorithm controls the matching elements (vacuum capacitors) and thus their corresponding S-matrices. Typically a 1-second pulse requires around 10 seconds of computational time on a desktop computer. These tasks can be hardly handled by commercial RF software. This innovative work allows identifying strategies for the launchers future operation while insuring the limitations on the currents, voltages and electric fields, matching and Load-Resilience, as well as the required straps voltage amplitude/phase balance. In this paper, a particular attention is paid to the simulation of the launchers behavior when arcs appear at several locations of their circuits using SIDON calculator. This latter work shall confirm or identify strategies for the arc detection using various RF electrical signals. One shall note that the use of such solvers in not limited to ICRF launchers simulations but can be employed, in principle, to any linear or linearized RF problem.

The paper reports results of MCNP-5 calculations to assess the effect on the Tritium Breeding Ratio (TBR) when integrating a distributed Ion Cyclotron Range of Frequencies (ICRF) antenna in the blanket of DEMO fusion power reactor. The calculations consider different parameters such as the ICRF covering ratio and the type of breeding blanket including the Helium Cooled Pebble Bed (HCPB) and the Helium Cooled Lithium Lead (HCLL) concepts. For an antenna with a full toroidal circumference of 360°, located poloidally at 40° with a poloidal extension of 1 m, the reduction of the TBR is −0.349% for the HCPB blanket and −0.532% for the HCLL blanket. The distributed ICRF antenna is thus shown to have only a marginal effect on the TBR of the DEMO reactor.

RF sheaths rectification can be the reason for operational limits for Ion Cyclotron Range of Frequencies (ICRF) heating systems via impurity production or excessive heat loads. To simulate this process in realistic geometry, the Self-consistent Sheaths and Waves for Ion Cyclotron Heating (SSWICH) code is a minimal set of coupled equations that computes self-consistently wave propagation and DC plasma biasing. The present version of its wave propagation module only deals with the Slow Wave assumed to be the source of RF sheath oscillations. However the ICRF power coupling to the plasma is due to the fast wave (FW). This paper proposes to replace this one wave equation module by a full wave module in either 2D or 3D as a first step towards integrated modelling of RF sheaths and wave coupling. Since the FW is propagative in the main plasma, Perfectly Matched Layers (PMLs) adapted for plasmas were implemented at the inner side of the simulation domain to absorb outgoing waves and tested numerically with tilted BD in Cartesian geometry, by either rotating the cold magnetized plasma dielectric tensors in 2D or rotating the coordinate vector basis in 3D. The PML was further formulated in cylindrical coordinates to account for for the toroidal curvature of the plasma. Toroidal curvature itself does not seem to change much the coupling. A detailed 3D geometrical description of Tore Supra and ASDEX Upgrade (AUG) antennas was included in the coupling code. The full antenna structure was introduced, since its toroidal symmetry with respect to the septum plane is broken (FS bars, toroidal phasing, non-symmetrical structure). Reliable convergence has been obtained with the density profile up to the leading edge of antenna limiters. Parallel electric field maps have been obtained as an input for the present version of SSWICH.

The paper reports on latest developments for the ITER Ion Cyclotron Heating and Current Drive system: imminent acceptance tests of a prototype power supply at full power; successful factory acceptance of candidate RF amplifier tubes which will be tested on dedicated facilities; further design integration and experimental validation of transmission line components under 6MW hour-long pulses. The antenna Faraday shield thermal design has been validated above requirements by cyclic high heat flux tests. R&D on ceramic brazing is under way for the RF vacuum windows. The antenna port plug RF design is stable but major evolution of the mechanical design is in preparation to achieve compliance with the load specification, warrant manufacturability and incorporate late interface change requests. The antenna power coupling capability predictions have been strengthened by showing that, if the plasma scrape-off layer turns out to be steep and the edge density low, the reference burning plasma can realistically be displaced to improve the coupling.

In tokamak fusion plasmas, coherent fluctuations in the form of blobs or filaments are routinely observed in the scrape-off layer. In this paper we develop an analytical formalism for the scattering of radio frequency waves by filaments which are cylindrical with their major axis aligned along the toroidal magnetic field lines. Since the magnitude of the ratio of the density inside the filaments to the background density is generally of order 1, the geometric optics approximation cannot be used to describe the scattering. A full-wave model is formulated which assumes that the plasma is cold and that the plasma in the cylindrical filament has uniform density. The background plasma, in which the filament is present, is also assumed to be cold and uniform. The theoretical framework applies to the scattering of any plasma wave.

The physics of wave heating and current drive processes in the bulk hot plasma are generally well identified. However, details of the wave-plasma interaction in the cold plasma edge are still not fully understood. In this paper a spectroscopic technique allowing for measurements of the electric field driving wave-plasma interactions in the edge region for LH current drive and ICRH systems will be discussed. Experimental results obtained in Tore Supra near the LH C3 antenna and in the magnetized capacitively coupled RF sheath from an experiment designed to mimic the ICRF near-field are presented. In conclusion, future plans to implement a laser-based spectroscopic technique to acquire the high resolution measurements needed for model validation is discussed.

With the return to full metal walls - a necessary step towards viable fusion machines - and due to the high power densities of current-day ICRH (Ion Cyclotron Resonance Heating) or RF (radio frequency) antennas, there is ample renewed interest in exploring the reasons for wave-induced sputtering and formation of hot spots. Moreover, there is experimental evidence on various machines that RF waves influence the density profile close to the wave launchers so that waves indirectly influence their own coupling efficiency. The present study presents a return to first principles and describes the wave-particle interaction using a 2-time scale model involving the equation of motion, the continuity equation and the wave equation on each of the time scales. Through the changing density pattern, the fast time scale dynamics is affected by the slow time scale events. In turn, the slow time scale density and flows are modified by the presence of the RF waves through quasilinear terms. Although finite zero order flows are identified, the usual cold plasma dielectric tensor - ignoring such flows - is adopted as a first approximation to describe the wave response to the RF driver. The resulting set of equations is composed of linear and nonlinear equations and is tackled in 1D in the present paper. Whereas the former can be solved using standard numerical techniques, the latter require special handling. At the price of multiple iterations, a simple ’derivative switch-on’ procedure allows to reformulate the nonlinear problem as a sequence of linear problems. Analytical expressions allow a first crude assessment - revealing that the ponderomotive potential plays a role similar to that of the electrostatic potential arising from charge separation - but numerical implementation is required to get a feeling of the full dynamics. A few tentative examples are provided to illustrate the phenomena involved.

ITER is an international project to design and build an experimental fusion reactor based on the “tokamak” concept. ITER relies upon localized electron cyclotron current drive (ECCD) at the rational safety factor q=2 to suppress or stabilize the expected poloidal mode m=2, toroidal mode n=1 neoclassical tearing mode (NTM) islands. Such islands if unmitigated degrade energy confinement, lock to the resistive wall (stop rotating), cause loss of “H-mode” and induce disruption. The International Tokamak Physics Activity (ITPA) on MHD, Disruptions and Magnetic Control joint experiment group MDC-8 on Current Drive Prevention/Stabilization of Neoclassical Tearing Modes started in 2005, after which assessments were made for the requirements for ECCD needed in ITER, particularly that of rf power and alignment on q=2 [1]. Narrow well-aligned rf current parallel to and of order of one percent of the total plasma current is needed to replace the “missing” current in the island O-points and heal or preempt (avoid destabilization by applying ECCD on q=2 in absence of the mode) the island [2-4]. This paper updates the advances in ECCD stabilization on NTMs learned in DIII-D experiments and modeling during the last 5 to 10 years as applies to stabilization by localized ECCD of tearing modes in ITER. This includes the ECCD (inside the q=1 radius) stabilization of the NTM “seeding” instability known as sawteeth (m/n=1/1) [5]. Recent measurements in DIII-D show that the ITER-similar current profile is classically unstable, curvature stabilization must not be neglected, and the small island width stabilization effect from helical ion polarization currents is stronger than was previously thought [6]. The consequences of updated assumptions in ITER modeling of the minimum well-aligned ECCD power needed are all-in-all favorable (and well-within the ITER 24 gyrotron capability) when all effects are included. However, a “wild card” may be broadening of the localized ECCD by the presence of the island; various theories predict broadening could occur and there is experimental evidence for broadening in DIII-D. Wider than now expected ECCD in ITER would make alignment easier to do but weaken the stabilization and thus require more rf power. In addition to updated modeling for ITER, advances in the ITER-relevant DIII-D ECCD gyrotron launch mirror control system hardware and real-time plasma control system have been made [7] and there are plans for application in DIII-D ITER demonstration discharges.

Several experiments on different machines and in different fast wave (FW) heating regimes, such as hydrogen minority heating and high harmonic fast waves, have found strong interactions between radio-frequency (RF) waves and the scrape-off layer (SOL) region. This paper examines the propagation and the power loss in the SOL by using the full wave code AORSA, in which the edge plasma beyond the last closed flux surface (LCFS) is included in the solution domain and a collisional damping parameter is used as a proxy to represent the real, and most likely nonlinear, damping processes. 3D AORSA results for the National Spherical Torus eXperiment (NSTX), where a full antenna spectrum is reconstructed, are shown, confirming the same behavior found for a single toroidal mode results in Bertelli et al, Nucl. Fusion, 54 083004, 2014, namely, a strong transition to higher SOL power losses (driven by the RF field) when the FW cut-off is moved away from in front of the antenna by increasing the edge density. Additiona...