Contributions to Plasma Physics [SJR: 0.717] [H-I: 36] [3 followers] Follow Hybrid journal (It can contain Open Access articles) ISSN (Print) 0863-1042 - ISSN (Online) 1521-3986 Published by John Wiley and Sons [1611 journals] |
- Study on a New Ignition Method of the Spark Gap Based on Plasma Ejection
in Air- Authors: D. Huang; L.J. Yang, P. Huo, J.B. Ma, S. Liu, W. Wang, W.D. Ding
Abstract: A new ignition method of the spark gap based on plasma ejection is proposed in this paper, as the conventional trigatron spark gap performs poorly under the low working coefficient (the ratio of the charging voltage to the self‐breakdown voltage) in air. The plasma is generated by the capillary discharge, which has high pressure, high temperature and high velocity. The capillary discharge device is placed inside the low voltage electrode. As long as the triggering signal is sent to the device, a column of the plasma flow is ejected in axial direction and develops rapidly towards the high voltage electrode. Subsequently, the gap is broken down and a high resistive channel is formed, where the thermal ionization takes place and the arc across the whole gap is generated and develops into a well conductive channel. The process of the thermal ionization of the high resistive channel varies with the change of the spark gap distance. The breakdown delay and the delay jitter of the spark gap increase with the spark gap distance, as both parameters are mainly determined by the developing process of the plasma ejection. The characteristics of the plasma flow determine the possibility of the breakdown of the spark gap under the low working coefficient. The ignition method based on capillary plasma ejection has been proved by the preliminary experiments, which indicate that under the gap length of 8 cm and the working coefficient of less than 3%, the effective ignition is still achievable.
PubDate: 2016-09-12T01:50:11.389532-05:
DOI: 10.1002/ctpp.201600037
- Authors: D. Huang; L.J. Yang, P. Huo, J.B. Ma, S. Liu, W. Wang, W.D. Ding
- An Innovative Simple Method for Study of the Characteristics of the
Trigatron Plasma Switch- Authors: M. Aram; N. Morshedian, S. Behrouzinia, M. Namnabat
Abstract: A simple Semi‐Analytical method used to fit the experimentally recorded current of the closed switch discharge circuit into the free running under damped LC oscillator model and the arc plasma conductivity, electron density of the plasma and efficiency of the spark gap switch energy transfer have been driven from the model by some simple calculations. The charging voltage, switch pressure and the spark gap has been changed between 6 to 15 kV, 1 to 1.5 bar and 1.4 or 2.8 mm, respectively. The obtained values for plasma conductivity and electron density of plasma are (5‐25)(Ω mm)–1, (0.5–3.5)·1024m–3, respectively. The efficiency of the switch is plotted Vs. Pd and E/P which in both cases the peak value has been about 80 present.
PubDate: 2016-09-12T01:40:17.938978-05:
DOI: 10.1002/ctpp.201600023
- Authors: M. Aram; N. Morshedian, S. Behrouzinia, M. Namnabat
- Issue Information: Contrib. Plasma Phys. 6–8/2016
- PubDate: 2016-08-23T04:34:17.069572-05:
DOI: 10.1002/ctpp.201690009
- PubDate: 2016-08-23T04:34:17.069572-05:
- Cover Picture: Contrib. Plasma Phys. 6–8/2016
- Abstract: Influence of the shear‐layer width on the nonlinear oscillations, spatiotemporal plot of the instantaneous convective flux Qconv with the shear layer‐width dE = 0.02. For details see paper of M. Leconte et al.
PubDate: 2016-08-23T04:34:12.514418-05:
DOI: 10.1002/ctpp.201690008
- Abstract: Influence of the shear‐layer width on the nonlinear oscillations, spatiotemporal plot of the instantaneous convective flux Qconv with the shear layer‐width dE = 0.02. For details see paper of M. Leconte et al.
- Nonlinear Dynamics of Circularly Polarized Laser Pulse Propagating in a
Magnetized Plasma with q ‐Nonextensive Velocity Distributions- Authors: R. Etemadpour; N. Sepehri Javan, D. Dorranian
Abstract: The nonlinear dynamics of a circularly polarized laser pulse propagating in magnetized plasma contains hot nonextensive q ‐distributed electrons and ions is studied theoretically. A nonlinear equation which describes the dynamics of the slowly varying amplitude electromagnetic wave is obtained using the relativistic two‐fluids model. Some nonlinear phenomena include modulational instability, self‐focusing, soliton formation, and longitudinal and transversal evolutions of laser pulse in nonextensive plasma medium are investigated. Results show that the nonextensivity of particles can substantially change the nonlinearity of medium. The external magnetic field enhances the modulation instability growth rate of right‐hand polarization wave but for the left‐hand polarization the growth rate decreases. The spot size of the laser pulse is strongly affected by the plasma nonextensivity. (© 2016 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)
PubDate: 2016-08-09T02:40:04.699318-05:
DOI: 10.1002/ctpp.201600010
- Authors: R. Etemadpour; N. Sepehri Javan, D. Dorranian
- A Simplified MHD Model of Capillary Z‐Pinch Compared with
Experiments- Authors: A. A. Shapolov; M. Kiss, S. V. Kukhlevsky
Abstract: The most accurate models of the capillary Z‐pinches used for excitation of soft X‐ray lasers and photolithography XUV sources currently are based on the magnetohydrodynamics theory (MHD). The output of MHD‐based models greatly depends on details in the mathematical description, such as initial and boundary conditions, approximations of plasma parameters, etc. Small experimental groups who develop soft X‐ray/XUV sources often use the simplest Z‐pinch models for analysis of their experimental results, despite of these models are inconsistent with the MHD equations. In the present study, keeping only the essential terms in the MHD equations, we obtained a simplified MHD model of cylindrically symmetric capillary Z‐pinch. The model gives accurate results compared to experiments with argon plasmas, and provides simple analysis of temporal evolution of main plasma parameters. The results clarify the influence of viscosity, heat flux and approximations of plasma conductivity on the dynamics of capillary Z‐pinch plasmas. The model can be useful for researchers, especially experimentalists, who develop the soft X‐ray/XUV sources. (© 2015 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)
PubDate: 2016-07-11T01:30:08.623346-05:
DOI: 10.1002/ctpp.201600025
- Authors: A. A. Shapolov; M. Kiss, S. V. Kukhlevsky
- Parametric Study of Nonplanar Dust Acoustic Solitons in Two‐Dimensional
Superthermal Dusty Plasmas- Authors: S. Khan; Ata-ur-Rahman, S. Ali, M. Khalid, N. A. D. Khattak, A. Zeb
Abstract: The nonlinear properties of two dimensional low‐frequency electrostatic excitations of charged dust particles (or defects) are studied in a collisionless, unmagnetized dusty plasma. A fully ionized three‐component model plasma consisting of kappa distributed electrons, Maxwellian ions, and negatively charged massive dust grains is considered. In this regard, the well known reductive perturbation technique is used to the hydrodynamical equations and the Poisson equation, obtaining the cylindrical Kadomtsev–Petviashvili (CKP) equation. A parametric investigation indicates that the structural characteristics of these nonlinear excitations (width, amplitude) are significantly affected by the plasma nonthermality as well as by the relevant plasma parameters, such as dust concentration and dust temperature. (© 2015 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)
PubDate: 2016-07-11T01:30:07.783243-05:
DOI: 10.1002/ctpp.201600016
- Authors: S. Khan; Ata-ur-Rahman, S. Ali, M. Khalid, N. A. D. Khattak, A. Zeb
- Electrostatic Fluid Accelerator under High‐Speed Free Air Stream
- Authors: T. Y. Wen; K. Lindgren, A. V. Mamishev
Abstract: The electrostatic fluid accelerator (EFA) generates ionic wind with a simple structure that barely obstructs the free air stream or produce excessive noise. This paper presents the velocity characteristics of an EFA under a high speed free air stream to simulate an EFA‐powered propulsor. The results show that when the EFA generates identical velocity to the free air stream, the EFA contributes 25% of the resultant velocity. When the EFA is replaced by a rotary fan that generates identical velocity to the free air stream, the fan contributes only 13.4% of the resultant velocity. (© 2015 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)
PubDate: 2016-07-11T01:30:06.895214-05:
DOI: 10.1002/ctpp.201600013
- Authors: T. Y. Wen; K. Lindgren, A. V. Mamishev
- Discharge Characteristics of Microwave Plasma in Ethanol Solution
- Authors: X. Zhao; B. Sun, X. Zhu, Z. Yan, H. Liu, Y. Liu
Abstract: An apparatus was designed for generating plasma in ethanol solution. The plasma was generated on the top of the electrode by applying microwave radiation of 2.45 GHz. The results showed that ignition power decreased with increasing temperature of ethanol solution. However, ignition power increased with increasing pressure and point electrode radius of curvature. Plasma and bubbles were generated periodically in the same manner. The electron temperature of the plasma increased with increasing power, while it decreased with increasing pressure. (© 2015 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)
PubDate: 2016-07-11T01:30:05.901933-05:
DOI: 10.1002/ctpp.201600011
- Authors: X. Zhao; B. Sun, X. Zhu, Z. Yan, H. Liu, Y. Liu
- Long‐Time Evolution of a Low Pressure Laboratory Plasma after
Application of Transient high Voltage Positive Pulses- Authors: M. Choudhary; S. Kar, S. Mukherjee
Abstract: The long‐time evolution of weakly‐collisional plasma with application of high voltage positive pulses to an electrode immersed in plasma, with pulse widths less than as well as more than ion plasma periods, is studied. The plasma is produced by electron impact ionization of argon or helium gas, where electrons are coming out from dc biased hot thoriated tungsten filaments. It is observed that during the temporal evolution of argon plasma, a beam component exists along with temporal bulk electrons giving rise to a double hump profile of transient Electron Distribution Function (EDF). However, in the case of temporal evolution of helium plasma, only a bulk electron population is present. The obtained results are explained by understanding the role played by thermionically emitted electrons during the plasma evolution, the role of the difference of ionization rates of helium and argon, and the higher temporal plasma potential. (© 2016 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)
PubDate: 2016-06-24T05:10:07.936983-05:
DOI: 10.1002/ctpp.201500146
- Authors: M. Choudhary; S. Kar, S. Mukherjee
- Transport and Energy Deposition of hot Electrons through the Shock
Ignition pre Compressed Target- Authors: S. Rezaei; A. H. Farahbod, M. J. Jafari
Abstract: Shock ignition as an alternative scheme of the laser fusion has the potential of achieving efficient implosion. However, hot electrons produced in result of ignitor‐corona interaction may penetrate deep into the fuel making the compression less effective. Transport and energy deposition of hot electron beam into the dense pre compressed of HiPER target by means of Monte Carlo approach are discussed considering the influence of real density and electron beam characteristics. The target parameters before igniting the hot spot have been extracted from a fluid code and used as the initial profile for Monte Carlo simulations. In comparison with simplified step like density profile, electrons penetrate slightly deeper in the case of real shaped density profile. In addition, deposition zone of a broad spectrum electron beam is wider while, monoenergetic electrons depose their energy locally resulting more maximum energy deposition value. (© 2016 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)
PubDate: 2016-06-24T04:30:14.591074-05:
DOI: 10.1002/ctpp.201600026
- Authors: S. Rezaei; A. H. Farahbod, M. J. Jafari
- Excitation of Highly Charged Ions in Plasma by Ultrashort Electromagnetic
Pulses- Authors: V. A. Astapenko; S. V. Sakhno, V. S. Lisitsa
Abstract: Excitation of highly charged ions in plasmas under the action of ultrashort electromagnetic pulses is investigated theoretically in the frame of perturbation theory. The study is based on analytical expression for probability of bound‐bound transition which is derived for Gaussian pulse and Doppler broadening of spectral line shape of the transition. Electronic transitions from the ground state of hydrogen‐like ions are considered with account for fine splitting of upper energy levels. The main attention is paid to the dependence of the excitation probability on pulse duration for various ion charges and carrier frequencies of the pulse. The results obtained are of interest for plasma diagnostics based on ion excitation from ground state by ultrashort pulses. (© 2016 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)
PubDate: 2016-06-24T04:30:13.824457-05:
DOI: 10.1002/ctpp.201600018
- Authors: V. A. Astapenko; S. V. Sakhno, V. S. Lisitsa
- Second Harmonic Generation of Self‐Focused Cosh‐Gaussian Laser Beam in
Thermal Quantum Plasma by Excitation of an Electron Plasma Wave- Authors: N. Gupta; A. Singh
Abstract: This paper presents a scheme for second harmonic generation (SHG) of an intense Cosh‐Gaussian (ChG) laser beam in thermal quantum plasmas. Moment theory approach in W.K.B approximation has been adopted in deriving the differential equation governing the propagation characteristics of the laser beam with distance of propagation. The effect of relativistic increase in electron mass on propagation dynamics of laser beam has been incorporated. Due to relativistic nonlinearity in the dielectric properties of the plasma, the laser beam gets self‐focused and produces density gradients in the transverse direction. The generated density gradients excite electron plasma wave (EPW) at pump frequency that interacts with the incident laser beam to produce its second harmonics. Numerical simulations have been carried out to investigate the effects of laser parameters on selffocusing of the laser beam and hence on the conversion efficiency of its second harmonics. Simulation results predict that within a specific range of decentered parameter the ChG laser beams show smaller divergence as they propagate and, thus, lead to enhanced conversion efficiency of second harmonics. (© 2016 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)
PubDate: 2016-06-24T04:30:12.988933-05:
DOI: 10.1002/ctpp.201600008
- Authors: N. Gupta; A. Singh
- Control of High Power Pulsed Magnetron Discharge by Monitoring the Current
Voltage Characteristics- Authors: A. von Keudell; A. Hecimovic, C. Maszl
Abstract: Discharge current voltage (IV) curves are directly measured at the target of a high impulse power magnetron sputtering (HiPIMS) plasma for the target materials aluminium, chromium, titanium and copper. These discharge IV curves have been correlated with ICCD camera images of the plasma torus. A clear connection between the change in the discharge IV curve slopes at specific currents and the appearance of localized ionization zones, so‐called spokes, in a HiPIMS plasma is identified. These spokes appear above typical target current densities of 2 A/cm2. The slope of the discharge IV curves, at current densities when spokes are formed, depends on the mass of the target atoms with a higher plasma conductivity for higher mass target materials. This is explained by the momentum transfer from the sputter wind to the argon background gas, which leads to higher plasma densities for heavier target materials. The change in the VI curve slope can be used to identify the spokes regime for HiPIMS plasmas, as being mandatory for deposition of good quality materials by HiPIMS. Consequently, the discharge IV curve slope monitoring can be regarded an essential control approach of any industrial HiPIMS process, where discharge IV curves are much easier accessible compared to more complex diagnostics such as time and space resolved ICCD camera measurements. (© 2016 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)
PubDate: 2016-06-24T04:30:12.078361-05:
DOI: 10.1002/ctpp.201600004
- Authors: A. von Keudell; A. Hecimovic, C. Maszl
- Usage of Ray Tracing Transfer Matrix to Mitigate the Stray Light for ITER
Spectroscopy- Authors: S. Kajita; E. Veshchev, R. Barnsley, M. Walsh
Abstract: Stray light formed by the reflection of photons on inner wall from a bright divertor region can be a serious issue in spectroscopic measurement systems in ITER. In this study, we propose a method to mitigate the influence of stray light using a ray tracing analysis. Usually, a ray tracing simulation requires a time consuming runs. We constructed transfer matrices based on the ray tracing simulation results and used them to demonstrate the influence of stray light. It is shown that the transfer matrix can be used to reconstruct the emission profile by considering the influence of the stray light without any additional ray tracing runs. Mitigation of the stray light in ITER divertor impurity monitor was demonstrated, and a method of prediction of the stray light level for the scrape off layer spectroscopy from divertor region was proposed. (© 2016 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)
PubDate: 2016-05-04T06:10:09.086432-05:
DOI: 10.1002/ctpp.201500124
- Authors: S. Kajita; E. Veshchev, R. Barnsley, M. Walsh
- Numerical Modeling on Heat Transport Across Stochastic Magnetic Field
- Authors: Q. Wang; X. Zha, H. Lu, X. Wang, B. Wu, S. Zhu
Abstract: The heat diffusion across the stochastic magnetic field is studied numerically. The stochastic field is generated by the overlap of two magnetic islands. The parameter w/wc, is found tobe an important parameter in charactering the transport, where w is the magnetic island width, and wc is the critical island width for flattening the electron temperature across an island. For w/wc < 1, the enhanced radial heat diffusivity χr is proportional to the parallel heat diffusivity χ∥, and the heat transport is dominated by the additive effect of individual islands. For w/wc > 3, χr is also proportional to χ∥ and the additional degradation of the energy confinement due to stochastic magnetic field becomes apparent. (© 2016 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)
PubDate: 2016-05-04T06:10:08.318582-05:
DOI: 10.1002/ctpp.201500087
- Authors: Q. Wang; X. Zha, H. Lu, X. Wang, B. Wu, S. Zhu
- Simulation and Experimental Study for Degradation of Organic Dyes Using
Dual pin‐to‐plate Corona Discharge Plasma reactors for Industrial
Wastewater Treatment- Authors: A. EL-Tayeb; A.H. EL-Shazly, M.F. Elkady, A. Abdel-Rahman
Abstract: This work researches the possibility of increasing the dye removal efficiency from wastewater using nonthermal plasma. A study for the optimal air gap distance between dual pin and surface of Acid Blue 25 dye solution and thickness of ground plate is carried out using 3D‐EM simulator to find maximum electric field intensity at the tip of both pins. The consequences display that the best gap for corona discharge is approximately 5 mm using 15 kV source. In addition, the optimum plate thickness is 0.1 mm. These distance and thickness were mentioned are constant during the study of other factors. Dual pin‐to‐plate high‐voltage corona discharge plasma system is presented to investigation experimentally the gap distance, thickness of ground plate, initial dye concentration, pH solution and conductivity on the amount of Acid Blue 25 dye color removal efficiency from wastewater. There is a large consensus among the simulation and experimental work in the air gap and thickness of ground plate. Where the decolorization for air gap 5 mm is 95.74 at time 35 min compared with 91% and 17% for 1 mm and 20 mm gap distance respectively. Also, the discharge energy at each air gap are calculated. Measurement results for the impact of thickness of an Aluminum ground plate on color removal competence showed color removal efficiencies of 86.3%, 90.78% and 98.06%, after treatment time 15 min for thicknesses of 2, 0.5 and 0.1 mm respectively. The decolorization behavior utilizing dual pin‐to‐plate corona discharge plasma system display 82% pigment evacuation proficiency inside 11min. The complete decolorization was accomplished within 28min for distinctive examined introductory color focuses 5 ppm up to 100 ppm. Likewise, the impacts of conductivity by utilizing diverse salts as AlCl3, CaCl2, KCl and NaCl and with distinctive focuses have been explored. The rising of the solution conductivity leads to the reduction of decolorization efficiency. The decolorization efficiency and discharge energy are calculated at different concentration molarity for AlCl3, CaCl2, KCl and NaCl. It was observed that the presence of salts at the same concentration level substantially decreased the rate and the extent of decolorization. The results indicate that the optimum pH for the decolorization of Acid Blue 25 dye is in the range between 3 and 6. Furthermore the conductivity and discharge energy were measurement at each value of pH. Energy yield for decolorization and Electrical Energy per Order (EE/O) under different initial pH value were calculated. A kinetic model is used to define the performance of corona discharge system under different value of pH. The model of pseudo ‐zero, pseudo‐first order, and pseudo‐second order reactions kinetic are utilized to investigate the decolorization of Acid Blue 25 dye. The rate of degradation of Acid Blue 25 dye follows the pseudo‐first order kinetics in the dye concentration. Energy consumption requirements for decolorization was considered. The outcomes will be useful for designing the plasma treatment systems suitable for industrial wastewaters. (© 2016 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)
PubDate: 2016-05-04T06:10:07.212529-05:
DOI: 10.1002/ctpp.201500080
- Authors: A. EL-Tayeb; A.H. EL-Shazly, M.F. Elkady, A. Abdel-Rahman
- Statistical Model of Dielectronic Recombination of Heavy Ions in Plasmas
- Authors: D. S. Leontyev; V. S. Lisitsa
Abstract: Calculation of the total dielectronic recombination (DR) rates was done in the frame of a statistical model of atoms. The model is based on the idea of collective excitations of atomic electrons with the local plasma frequency, which depends on atomic electrons density distribution. The electron density is described in a frame of the Thomas‐Fermi model of atoms. Simple scaling laws for temperature Te and nuclear charge Z dependences follow from the statistical model of DR. Results of the statistical model were compared with other numerical data following detailed level‐by‐level computations for different multielectron ions. The specific attention is paid to Ni‐like ion sequences of different chemical elements in order to check the Z ‐dependence of DR rates. A comparison with numerical data of Flexible Atomic Code (FAC) is presented for tungsten ions. The reasonable correspondence between the statistical model and the detailed numerical data is demonstrated. The application of the statistical model provides very simple and fast calculations of the DR rates useful in modern plasma modelling. (© 2016 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)
PubDate: 2016-05-04T06:10:06.34635-05:0
DOI: 10.1002/ctpp.201500075
- Authors: D. S. Leontyev; V. S. Lisitsa
- The Organosilicon thin Film Deposited Using an Atmospheric Pressure
Dual‐Frequency 50 kHz/33 MHz Frequency MicroPlasma Jet- Authors: Q. H. Yuan; X. M. Wang, G. Q. Yin, J. Li, C. Z. Dong
Abstract: In this paper, an atmospheric pressure dual‐frequency (50 kHz/33 MHz) micro‐plasma jet was used to deposit organosilicon film. The discharge generated in atmospheric environment. Plasma composition was characterized by optical emission spectroscopy. With introduction of tetraethyl orthosilicate, we observed various spectra, for example Si(251.6 nm), OH(308.9 nm), C(247.8 nm), O(777.5 nm). Abundant reactive radical species which are benefit to film deposition were generated in plasma. The deposited film was characterized by scanning electron microscopy, X‐ray photoelectron spectroscopy and Fourier transform infrared spectroscopy. The film is mostly composed of Si and O. The film has Si‐O‐Si backbone with a small number of organic component (‐CHx). (© 2016 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)
PubDate: 2016-05-04T06:10:05.50579-05:0
DOI: 10.1002/ctpp.201500066
- Authors: Q. H. Yuan; X. M. Wang, G. Q. Yin, J. Li, C. Z. Dong
- Factorization of 3‐Point Static Structure Functions in 3D Yukawa
Liquids- Authors: P. Magyar; P. Hartmann, G.J. Kalman, K. I. Golden, Z. Donkó
Abstract: In many‐body systems the convolution approximation states that the 3‐point static structure function, S(3)(k1, k2), can approximately be “factorized” in terms of the 2‐point counterpart, S(2)(k1). We investigate the validity of this approximation in 3‐dimensional strongly‐coupled Yukawa liquids: the factorization is tested for specific arrangements of the wave vectors k1 and k2, with molecular dynamics simulations. With the increase of the coupling parameter we find a breakdown of factorization, of which a notable example is the appearance of negative values of S(3)(k1, k2), whereas the approximate factorized form is restricted to positive values. These negative values – based on the quadratic Fluctuation‐Dissipation Theorem – imply that the quadratic part of the density response of the system changes sign with wave number. Our simulations that incorporate an external potential energy perturbation clearly confirm this behavior. (© 2016 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)
PubDate: 2016-05-04T06:10:04.411432-05:
DOI: 10.1002/ctpp.201500062
- Authors: P. Magyar; P. Hartmann, G.J. Kalman, K. I. Golden, Z. Donkó
- Contents: Contrib. Plasma Phys. 6–8/2016
- First page: 487
PubDate: 2016-08-23T04:34:15.263582-05:
DOI: 10.1002/ctpp.201610028
- First page: 487
- Preface: Contrib. Plasma Phys. 6–8/2016
- First page: 493
PubDate: 2016-08-23T04:34:11.900191-05:
DOI: 10.1002/ctpp.201610029
- First page: 493
- Recent Developments in Plasma Edge Theory
- Authors: W. M. Stacey
First page: 495
Abstract: Recent developments in electromagnetic particle pinch, ion orbit loss, intrinsic rotation, rotation theory and radial electric field theory in the tokamak plasma edge are described. (© 2016 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)
PubDate: 2016-07-12T05:50:08.893602-05:
DOI: 10.1002/ctpp.201610060
- Authors: W. M. Stacey
- Connecting Collisionless Landau Fluid Closures to Collisional Plasma
Physics Models- Authors: I. Joseph; A. M. Dimits
First page: 504
Abstract: In this work, collisionless Landau fluid models are combined with collisional fluid models using a Padé approximation that is accurate in both collisionality limits. The collisionless models capture Landau damping through a nonlocal integro‐differential sub‐diffusion operator with ballistic characteristics. Collisional extensions of Landau fluid models are derived by analyzing the higher order moment equations which combine a Landau closure with collisional friction forces. The model derived here evolves fluid moments for density, parallel velocity, and anisotropic pressure and includes the frictional heat flux, the parallel thermal force and anisotropic electrical conductivity. Since anisotropies must vanish in the collisional limit, a simple closure can be derived if the friction force neglects pitch‐angle scattering for the closure moments themselves. The resulting plasma physics model is potentially quite useful for applications in magnetic fusion and astrophysics. (© 2016 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)
PubDate: 2016-07-12T05:40:14.223165-05:
DOI: 10.1002/ctpp.201610043
- Authors: I. Joseph; A. M. Dimits
- Turbulence Dynamics with the Coupling of Density Gradient and Parallel
Velocity Gradient in the Edge Plasmas- Authors: Y. Kosuga; S.-I. Itoh, K. Itoh
First page: 511
Abstract: Theory of parallel shear flow driven instability (PSFI) and its impact on turbulence dynamics and transport are presented. The mode is linearly unstable when the parallel flow shear exceeds a critical value. The quasilinear particle flux contains both outward and inward components. Nonlinear dynamics is formulated in terms of hydrodynamic helicity balance. The result implies that once excited, PSFI with helicity may spread from the excited region to stable regions. Implication for turbulence in scrape off layer plasmas is discussed. (© 2016 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)
PubDate: 2016-06-24T04:40:11.161937-05:
DOI: 10.1002/ctpp.201610047
- Authors: Y. Kosuga; S.-I. Itoh, K. Itoh
- Basic Consideration of Monte Carlo Algorithm to Solve Fluid Equations for
SOL/divertor Plasmas- Authors: R. Tatsumi; Y. Homma, S. Yamoto, A. Hatayama
First page: 516
Abstract: Monte Carlo method is thought to be effective to solve fluid plasma equations for SOL/divertor plasmas, especially for three dimensional simulation. In the Monte Carlo algorithm based on a Lagrangian scheme, how to treat the Monte Carlo test particles at the calculation boundaries is not always trivial. In this paper, 1D diffusion equation with source terms has been solved with several different treatments of the boundaries in relatively a simple model. Comparison between the results and analytic solutions show that careful treatment of the boundary seems to be needed. (© 2016 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)
PubDate: 2016-07-12T05:50:08.031111-05:
DOI: 10.1002/ctpp.201610066
- Authors: R. Tatsumi; Y. Homma, S. Yamoto, A. Hatayama
- On the Origin of Steep and Localized Radial Electric Field in the
Transport Barrier at Plasma Edge- Authors: K. Itoh; S.-I. Itoh, T. Kobayashi, K. Kamiya, T. Ido, T. Tokuzawa, S. Inagaki, N. Kasuya
First page: 522
Abstract: Theories to understand the steep and localized radial electric field in the edge of toroidal plasma, which appears in conjunction with H‐mode, is revisited based on the electric field bifurcation model. Key elements in the models of the L‐H transition (including the toroidal effects on the dielectric constant and the effects of the curvature of radial electric field on turbulence suppression) are assessed. Results are applied to tokamak and helical plasmas, for which data with high‐resolution have been obtained recently. The status of quantitative tests on various mechanisms through comparison with experimental observations is also addressed. (© 2016 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)
PubDate: 2016-06-24T04:40:08.346523-05:
DOI: 10.1002/ctpp.201610042
- Authors: K. Itoh; S.-I. Itoh, T. Kobayashi, K. Kamiya, T. Ido, T. Tokuzawa, S. Inagaki, N. Kasuya
- Understanding of Hysteresis Behaviors at the L‐H‐L Transitions in
Tokamak Plasma Based on Bifurcation Concept- Authors: B. Chatthong; T. Onjun
First page: 528
Abstract: The hysteresis behaviour at the L‐H‐L transitions in tokamak plasma is investigated based on bifurcation concept. The formation of an edge transport barrier (ETB) is modeled via thermal and particle transport equations with the flow shear suppression effect on anomalous transport included. The anomalous transport is modeled based on critical gradients threshold and the flow shear is calculated from the force balance equation, couples the two transport equations leading to a non‐linear behaviour. Analytical investigation reveals that the fluxes versus gradients space exhibits bifurcation behaviour with s ‐curve soft bifurcation type. Apparently, the backward H‐L transition occurs at lower values than that of the forward L‐H transition, illustrating hysteresis behaviour. The hysteresis properties, i.e. locations of threshold fluxes, gradients and their ratios are analyzed as a function of neoclassical and anomalous transport values and critical gradients. It is found that the minimum heat flux for maintaining H ‐mode depends on several plasma parameters including the strength of anomalous transport and neoclassical transport. In particular, the hysteresis depth becomes larger when neoclassical transport decreases or anomalous transport increases. (© 2016 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)
PubDate: 2016-06-24T04:30:15.950077-05:
DOI: 10.1002/ctpp.201610034
- Authors: B. Chatthong; T. Onjun
- Gyrokinetic Theory and Dynamics of the Tokamak Edge
- Authors: B. Scott
First page: 534
Abstract: The validity of modern gyrokinetic field theory is assessed for the tokamak edge. The basic structure of the Lagrangian and resulting equations and their conservation laws is reviewed. The conventional microturbulence ordering for expansion is small potential/arbitrary wavelength. The equilibrium ordering for expansion is long wavelength/arbitrary amplitude. The long‐wavelength form of the conventional Lagrangian is derived in detail. The two Lagrangians are shown to match at long wavelength if the E × B Mach number is small enough for its corrections to the gyroaveraging to be neglected. Therefore, the conventional derivation and its Lagrangian can be used at all wavelengths if these conditions are satisfied. Additionally, dynamical compressibility of the magnetic field can be neglected if the plasma beta is small. This allows general use of a shear‐Alfvén Lagrangian for edge turbulence and self consistent equilibrium‐scale phenomena for flows, currents, and heat fluxes for conventional tokamaks without further modification by higher‐order terms. Corrections in polarisation and toroidal angular momentum transport due to these higher‐order terms for global edge turbulence computations are shown to be small. (© 2016 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)
PubDate: 2016-07-18T01:50:09.339327-05:
DOI: 10.1002/ctpp.201610057
- Authors: B. Scott
- Rubrik
- Authors: N. Miyato
First page: 543
Abstract: An energy‐conserving reduced kinetic model is formulated which is valid for perturbations with large amplitude and long wavelength beyond the standard gyrokinetic ordering. The standard gyrokinetic model is elegantly formulated by the two‐step phase space transformation consisting of the guiding‐center and gyro‐center transformations, but the amplitude of perturbations is assumed to be small in return for use of the two‐step transformation. Besides, the wavelength of perturbations is assumed to be short in deriving the gyrokinetic Poisson equation. Both limitations are relaxed by using the modified guiding‐center 1‐form for flowing plasmas [Miyato et al., J. Phys. Soc. Jpn. 78, 104501 (2009)]. (© 2016 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)
PubDate: 2016-06-24T04:40:13.573355-05:
DOI: 10.1002/ctpp.201610051
- Authors: N. Miyato
- Gyrokinetic Simulations of the Tokamak Plasma Edge in Circular Limiter
Configuration- Authors: T. Korpilo; T. P. Kiviniemi, S. Leerink, P. Niskala, R. Rochford
First page: 549
Abstract: The full‐f gyrokinetic code ELMFIRE is used to simulate the impact of turbulent tokamak plasma transport on the edge plasma flow and scrape‐off‐layer width. The simulation is performed in the circular limiter configuration and extends from the magnetic axis to the material surface. The results show that the sheath potential and parallel Mach number are in agreement with theoretical predictions, while the E × B dynamics are strongly affected by the sheath boundary. The radial fall‐off of density and temperature profiles show non‐exponential behaviour in the scrape‐off‐layer. Numerical issues appearing in the ELMFIRE simulation of edge plasma transport are discussed. (© 2016 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)
PubDate: 2016-07-12T05:40:13.449384-05:
DOI: 10.1002/ctpp.201610046
- Authors: T. Korpilo; T. P. Kiviniemi, S. Leerink, P. Niskala, R. Rochford
- Interchange Turbulence Model for the Edge Plasma in SOLEDGE2D‐EIRENE
- Authors: H. Bufferand; G. Ciraolo, Ph. Ghendrih, Y. Marandet, J. Bucalossi, C. Colin, N. Fedorczak, D. Galassi, J. Gunn, R. Leybros, E. Serre, P. Tamain
First page: 555
Abstract: Cross‐field transport in edge tokamak plasmas is known to be dominated by turbulent transport. A dedicated effort has been made to simulate this turbulent transport from first principle models but the numerical cost to run these simulations on the ITER scale remains prohibitive. Edge plasma transport study relies mostly nowadays on so‐called transport codes where the turbulent transport is taken into account using effective ad‐hoc diffusion coeffecients. In this contribution, we propose to introduce a transport equation for the turbulence intensity in SOLEDGE2D‐EIRENE to describe the interchange turbulence properties. Going beyond the empirical diffusive model, this system automatically generates profiles for the turbulent transport and hence reduces the number of degrees of freedom for edge plasma transport codes. We draw inspiration from the k‐epsilon model widely used in the neutral fluid community. (© 2016 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)
PubDate: 2016-07-12T05:40:17.322606-05:
DOI: 10.1002/ctpp.201610033
- Authors: H. Bufferand; G. Ciraolo, Ph. Ghendrih, Y. Marandet, J. Bucalossi, C. Colin, N. Fedorczak, D. Galassi, J. Gunn, R. Leybros, E. Serre, P. Tamain
- On the Effect of Electron Temperature Fluctuations on Turbulent Heat
Transport in the Edge Plasma of Tokamaks- Authors: C. Baudoin; P. Tamain, G. Ciraolo, R. Futtersack, A. Gallo, P. Ghendrih, Y. Marandet, N. Nace, C. Norscini
First page: 563
Abstract: In this paper we study the impact of electron temperature fluctuations in a two‐dimensional turbulent model. This modification adds a second linear instability, known as sheath‐driven conducting‐wall instability, with respect to the previous isothermal model only driven by the interchange instability. Non‐linear simulations, backed up by the linear analysis, show that the additional mechanism can change drastically the dynamics of turbulence (scales, density‐potential correlation, and statistical momentum). Moreover, its importance relatively to the interchange instability should be more significant in the private flux region than in the main scrape of layer. Its effect on heat transport is also investigated for different regimes of parameters, results show that both instabilities are at play in the heat transport. Finally, the sheath neagtive resistance instability could be responsible for the existence of corrugated heat flux profiles in the scrape‐off layer leading to a multiple decay length. (© 2016 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)
PubDate: 2016-07-12T05:40:15.832864-05:
DOI: 10.1002/ctpp.201610030
- Authors: C. Baudoin; P. Tamain, G. Ciraolo, R. Futtersack, A. Gallo, P. Ghendrih, Y. Marandet, N. Nace, C. Norscini
- Interplay between Plasma Turbulence and Particle Injection in 3D Global
Simulations- Authors: P. Tamain; H. Bufferand, L. Carbajal, Y. Marandet, C. Baudoin, G. Ciraolo, C. Colin, R. Futtersack, D. Galassi, P. Ghendrih, N. Nace, F. Schwander, E. Serre
First page: 569
Abstract: The impact of a 3D localized particle source on the edge plasma in 3D global turbulence simulations is investigated using the TOKAM3X fluid code. Results apply to advanced fueling methods such as Supersonic Molecular Beam Injection (SMBI) or pellets injection. The fueling source is imposed as a volumetric particle source in the simulations so that the physics leading to the ionization of particles and its localization are not taken into account. As already observed in experiments, the localized particle source strongly perturbs both turbulence and the large scale organization of the edge plasma. The localized increase of the pressure generated by the source drives sonic parallel flows in the plasma, leading to a poloidal redistribution of the particles on the time scale of the source duration. However, the particle deposition also drives localized transverse pressure gradients which impacts the stability of the plasma with respect to interchange processes. The resulting radial transport occurs on a sufficiently fast time scale to compete with the parallel redistribution of particles, leading to immediate radial losses of a significant proportion of the injected particles. Low Field Side (LFS) and High Field Side (HFS) injections exhibit different dynamics due to their interaction with curvature. In particular, HFS particle deposition drives an inward flux leading to differences in the particle deposition efficiency (higher for HFS than LFS). These results demonstrate the importance of taking into account plasma transport in a self‐consistent manner when investigating fueling methods. (© 2015 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)
PubDate: 2016-07-04T04:40:05.717577-05:
DOI: 10.1002/ctpp.201610063
- Authors: P. Tamain; H. Bufferand, L. Carbajal, Y. Marandet, C. Baudoin, G. Ciraolo, C. Colin, R. Futtersack, D. Galassi, P. Ghendrih, N. Nace, F. Schwander, E. Serre
- First Principle Modelling of Interplay between Langmuir Probes and Plasma
Turbulence- Authors: R. Futtersack; C. Colin, P. Tamain, G. Ciraolo, Ph. Ghendrih, Y. Marandet, F. Schwander, E. Serre
First page: 575
Abstract: The interplay between Langmuir probes (LP) and Scrape‐Off‐Layer plasma turbulence is numerically investigated with the TOKAM2D and TOKAM3X fluid codes. The LP is modelled by biasing a part of the target plates surface; we then study its impact on the turbulent transport 1) in presence of electron temperature fluctuations and 2) with a complete description of the parallel dynamics. We find that a biased probe can disturb local plasma parameters as well as turbulent transport in its vicinity, by polarizing the connected flux tube and thus driving a strong ExB vortex. Moreover, electron temperature fluctuations are found to account significantly those of floating potential, but with a limited impact on flux measurements depending on the probe's exact geometry. The 3D study of the problem shows the attenuation, but the persistency, of these perturbations induced by the presence of the LP. (© 2016 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)
PubDate: 2016-07-12T05:40:16.529766-05:
DOI: 10.1002/ctpp.201610038
- Authors: R. Futtersack; C. Colin, P. Tamain, G. Ciraolo, Ph. Ghendrih, Y. Marandet, F. Schwander, E. Serre
- Turbulence Interaction with Driven Transport Barriers in the
Scrape‐Off Layer of Tokamaks- Authors: N. Nace; M. Perin, P. Tamain, G. Ciraolo, C. Baudoin, R. Futtersack, Ph. Ghendrih, C. Norscini
First page: 581
Abstract: Two mechanisms able to generate a transport barrier by the stabilization of turbulence are analysed: the radial shear of the E × B velocity and the radial gradient of this shear. These two mechanisms are artificially forced in dedicated 2D non‐linear fluid simulations to study their impact on turbulence in the Scrape‐Off Layer. Generated barriers are characterized by their turbulence stopping capability and a key parameter, the barrier efficiency, is pointed out. The first two moments of this criterion can be studied to quantify the global efficiency, shape and width of the barrier as well as its intermittency. It is found that both mechanisms can drive intermittent or stable barriers but the radial shear mechanism generates stronger relaxation events. Small scale structures are damped by the E × B velocity shear while the radial shear gradient acts mainly on large scales. (© 2016 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)
PubDate: 2016-07-12T05:50:10.433275-05:
DOI: 10.1002/ctpp.201610052
- Authors: N. Nace; M. Perin, P. Tamain, G. Ciraolo, C. Baudoin, R. Futtersack, Ph. Ghendrih, C. Norscini
- Modeling of ITER Edge Plasma in the Presence of Resonant Magnetic
Perturbations- Authors: V. Rozhansky; E. Kaveeva, I. Veselova, S. Voskoboynikov, D. Coster
First page: 587
Abstract: The modeling of the ITER edge is performed with the use of the code B2SOLPS5.2 in the presence of the electron conductivity caused by RMPs as well as for the reference case with the same input parameters but without RMPs. The radial electric field close to the neoclassical one is obtained without RMPs. Even the modest level of RMPs changes the direction of the electric field and causes the toroidal spin‐up of the edge plasma. At the same time the pump‐out effect is small. (© 2016 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)
PubDate: 2016-07-12T05:50:09.652085-05:
DOI: 10.1002/ctpp.201610056
- Authors: V. Rozhansky; E. Kaveeva, I. Veselova, S. Voskoboynikov, D. Coster
- Development of a Drift‐Kinetic Simulation Code for Estimating
Collisional Transport Affected by RMPs and Radial Electric Field- Authors: R. Kanno; M. Nunami, S. Satake, S. Matsuoka, H. Takamaru
First page: 592
Abstract: A drift‐kinetic 8f simulation code is developed for estimating collisional transport in quasi‐steady state of toroidal plasma affected by resonant magnetic perturbations and radial electric field. In this paper, validity of the code is confirmed through several test calculations. It is found that radial electron flux is reduced by positive radial‐electric field, although radial diffusion of electron is strongly affected by chaotic field‐lines under an assumption of zero electric field. (© 2016 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)
PubDate: 2016-06-24T04:40:09.451311-05:
DOI: 10.1002/ctpp.201610044
- Authors: R. Kanno; M. Nunami, S. Satake, S. Matsuoka, H. Takamaru
- Modeling of Linear Divertor Plasma Simulator Experiments with
Three‐dimensional Target Structure by Using EMC3‐EIRENE Code- Authors: T. Kuwabara; H. Tanaka, G. Kawamura, N. Ohno, H. Nishikata, M. Kobayashi, Y. Feng
First page: 598
Abstract: We have adapted the EMC3‐EIRENE code for modeling of a linear divertor plasma simulator in order to demonstrate plasma‐wall interactions with three‐dimensional (3D) effects. 3D distributions of hydrogen plasma and neutrals can be successfully calculated for four different types of target plates: a V‐shaped target, inclined targets with open and closed structures, and a planer target. Hydrogen atoms and molecules are accumulated more effectively in the V‐shaped target plate, leading to a higher electron density with lower electron temperature than the planar target plate. (© 2016 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)
PubDate: 2016-06-24T04:40:11.925191-05:
DOI: 10.1002/ctpp.201610049
- Authors: T. Kuwabara; H. Tanaka, G. Kawamura, N. Ohno, H. Nishikata, M. Kobayashi, Y. Feng
- Effect of Statistical Noise on Simulation Results with a Plasma Fluid Code
Coupled to a Monte Carlo Kinetic Neutral Code- Authors: Y. Marandet; H. Bufferand, G. Ciraolo, P. Genesio, P. Meliga, J. Rosato, E. Serre, P. Tamain
First page: 604
Abstract: Power exhaust is one of the major challenges that future devices such as ITER and DEMO will face. Because of the lack of identified scaling parameters, predictions for divertor plasma conditions in these devices have to rely on detailed modelling. Most plasma edge simulations carried out so far rely on transport codes, which most of the times consist of a fluid code for the plasma coupled to a kinetic Monte Carlo (MC) code for neutral particles. One of the main difficulties in interpreting code results is the statistical noise from the MC procedure, which makes it difficult to define a convergence criterion for the simulations. In this work, we elaborate on similarities between noisy transport code simulations and turbulence simulations, and argue that the time averaged solution is a well defined stationary solution for the system. We illustrate these ideas with a simple slab test case with fluid neutrals, to which we add synthetic noise. In this case, the effects of noise are found to be significant only at high noise levels and for large enough correlations times.
PubDate: 2016-07-13T10:20:08.255218-05:
DOI: 10.1002/ctpp.201610009
- Authors: Y. Marandet; H. Bufferand, G. Ciraolo, P. Genesio, P. Meliga, J. Rosato, E. Serre, P. Tamain
- Fluid Neutral Model for Use in Hybrid Neutral Simulations of a Detached
Case- Authors: N. Horsten; W. Dekeyser, G. Samaey, P. Börner, M. Baelmans
First page: 610
Abstract: Two fluid neutral models are compared with a finite volume solution of the kinetic neutral equation for a simplified 1D detached case. The neutral equations are coupled to the plasma equations. The first fluid neutral model assumes that neutral and ion temperature are equal while the second adds a separate neutral energy equation. It is illustrated that both the AMJUEL‐HYDHEL databases and TRIM code can be incorporated. The first model is already accurate giving a maximum relative error of 10% for the target fluxes. Adding an energy equation further decreases this maximum relative error to 5%.
PubDate: 2016-07-13T10:20:03.950494-05:
DOI: 10.1002/ctpp.201610005
- Authors: N. Horsten; W. Dekeyser, G. Samaey, P. Börner, M. Baelmans
- Accuracy and Convergence of Coupled Finite‐Volume / Monte‐Carlo Codes
for Plasma Edge Simulations- Authors: K. Ghoos; W. Dekeyser, G. Samaey, P. Börner, D. Reiter, M. Baelmans
First page: 616
Abstract: Present computational techniques for coupled finite‐volume/Monte‐Carlo codes for plasma edge modeling under ITER or DEMO conditions face serious challenges with respect to computational time and accuracy. In this paper, scaling laws for different error contributions are assessed and practical procedures for error estimation are proposed. First results on a 1D and a 2D test case are discussed.
PubDate: 2016-07-13T10:20:07.570666-05:
DOI: 10.1002/ctpp.201610010
- Authors: K. Ghoos; W. Dekeyser, G. Samaey, P. Börner, D. Reiter, M. Baelmans
- Modelling of Impurity Transport and Plasma–Wall Interaction in Fusion
- Authors: A. Kirschner; D. Tskhakaya, G. Kawamura, D. Borodin, S. Brezinsek, R. Ding, Ch. Linsmeier, J. Romazanov
First page: 622
Abstract: The 3D ERO code, which simulates plasma–wall interaction and impurity transport in magnetically confined fusion‐relevant devices is described. As application, prompt deposition of eroded tungsten has been simulated at surfaces with shallow magnetic field of 3 T. Dedicated PIC simulations have been performed to calculate the characteristics of the sheath in front of plasma–exposed surfaces to use as input for these ERO simulations. Prompt deposition of tungsten reaches 100% at the highest electron temperature and density. In comparison to more simplified assumptions for the sheath the amount of prompt deposition is in general smaller if the PIC–calculated sheath is used. Due to friction with the background plasma the impact energy of deposited tungsten can be significantly larger than the energy gained in the sheath potential.
PubDate: 2016-07-13T10:20:04.739718-05:
DOI: 10.1002/ctpp.201610014
- Authors: A. Kirschner; D. Tskhakaya, G. Kawamura, D. Borodin, S. Brezinsek, R. Ding, Ch. Linsmeier, J. Romazanov
- EMC3‐EIRENE Simulation of Impurity Transport in Comparison with EUV
Emission Measurements in the Stochastic Layer of LHD: Effects of Force
Balance and Transport Coefficients- Authors: S. Dai; M. Kobayashi, G. Kawamura, S. Morita, T. Oishi, H. M. Zhang, X.L. Huang, Y. Feng, D. Z. Wang, Y. Suzuki, the LHD experiment group
First page: 628
Abstract: The transport properties and line emissions of the intrinsic carbon in the stochastic layer of the Large Helical Device have been investigated with the three‐dimensional edge transport code EMC3‐EIRENE. The simulations of impurity transport and emissivity have been performed to study the dedicated experiment in which the carbon emission distributions are measured by a space‐resolved EUV spectrometer system. A discrepancy of the CIV impurity emission between the measurement and simulation is obtained, which is studied with the variation of the ion thermal force, friction force and the perpendicular diffusivity in the impurity transport model. An enhanced ion thermal force or a reduced friction force in the modelling can increase the CIV impurity emission at the inboard X‐point region. Furthermore, the impact of the perpendicular diffusivity Dimp is studied which shows that the CIV impurity emission pattern is very sensitive to Dimp. It is found that the simulation results with the increased Dimp tend to be closer to the experimental observation. (© 2016 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)
PubDate: 2016-06-24T04:30:10.830718-05:
DOI: 10.1002/ctpp.201610036
- Authors: S. Dai; M. Kobayashi, G. Kawamura, S. Morita, T. Oishi, H. M. Zhang, X.L. Huang, Y. Feng, D. Z. Wang, Y. Suzuki, the LHD experiment group
- Extended Numerical Modeling of Impurity Neoclassical Transport in Tokamak
Edge Plasmas- Authors: H. Inoue; Y. Homma, S. Yamoto, A. Hatayama
First page: 634
Abstract: Understanding of impurity transport in tokamaks is an important issue in order to reduce the impurity contamination in fusion core plasmas. Recently, a new kinetic numerical scheme of impurity classical/neoclassical transport has been developed. This numerical scheme makes it possible to include classical self‐diffusion (CL SD), classical inward pinch (CL IWP), and classical temperature screening effect (CL TSE) of impurity ions. However, impurity neoclassical transport has been modeled only in the case where background plasmas are in the Pfirsch‐Schluter (PS) regime. The purpose of this study is to extend our previous model to wider range of collisionality regimes, i.e., not only the PS regime, but also the plateau regime. As in the previous study, a kinetic model with Binary Collision Monte‐Carlo Model (BMC) has been adopted. We focus on the modeling of the neoclassical self‐diffusion (NC SD) and the neoclassical inward pinch (NC IWP). In order to simulate the neoclassical transport with the BCM, velocity distribution of background plasma ions has been modeled as a deformed Maxwell distribution which includes plasma density gradient. Some test simulations have been done. As for NC SD of impurity ions, our scheme reproduces the dependence on the collisionality parameter in wide range of collisionality regime. As for NC IWP, in cases where test impurity ions and background ions are in the PS and plateau regimes, parameter dependences have been reproduced. (© 2016 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)
PubDate: 2016-07-12T05:40:14.876106-05:
DOI: 10.1002/ctpp.201610041
- Authors: H. Inoue; Y. Homma, S. Yamoto, A. Hatayama
- An Analytical Expression for the Electric Field and Particle Tracing in
Modelling of Be Erosion Experiments at the JET ITER‐like Wall- Authors: I. Borodkina; D. Borodin, A. Kirschner, I.V. Tsvetkov, V.A. Kurnaev, M. Komm, R. Dejarnac, JET Contributors
First page: 640
Abstract: A new analytical approximation for the electric potential profile in the presence of an oblique magnetic field and the analytical solution for the particle motion just before the impact with a plasma‐facing surface are presented. These approximations are in good agreement with fluid solutions and the corresponding PIC simulations. These expressions were applied to provide effective physical erosion yields for Be, which have in a second step been used in ERO code simulations of spectroscopy at Be limiters of the JET ITER‐like wall. These new analytical expressions lead to an increase of the effective physical sputtering yields of Be by deuteron impact up to 30% in comparison with earlier pure numerical simulations. (© 2016 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)
PubDate: 2016-06-24T04:30:15.26933-05:0
DOI: 10.1002/ctpp.201610032
- Authors: I. Borodkina; D. Borodin, A. Kirschner, I.V. Tsvetkov, V.A. Kurnaev, M. Komm, R. Dejarnac, JET Contributors
- Effects of Classical and Neo‐classical Cross‐field Transport of
Tungsten Impurity in Realistic Tokamak Geometry- Authors: S. Yamoto; Y. Homma, H. Inoue, Y. Sawada, K. Hoshino, A. Hatayama, X. Bonnin, D. Coster, R. Schneider
First page: 646
Abstract: The initial simulation study of the neoclassical perpendicular self‐diffusion transport in the SOL/Divertor regions for a realistic tokamak geometry with the IMPGYRO code has been performed in this paper. One of the most unique features of the IMPGYRO code is calculating exact Larmor orbit of the test particle instead of assuming guiding center approximation. Therefore, effects of the magnetic drifts in realistic tokamaks are naturally taken into account in the IMPGYRO code. This feature makes it possible to calculate neoclassical transport processes, which possibly become large in the SOL/divertor plasma. Indeed, neoclassical self‐diffusion process, the resultant effect of the combination of magnetic drift and Coulomb collisions with background ions, has already been included in the IMPGYRO model. In the present paper, prior to implementing the detailed model of neoclassical transport process into IMPGYRO, we have investigated the effect of neoclassical selfdiffusion in a realistic tokamak geometry with lower single null X‐point. We also use a model with guiding center approximation in order to compare with the IMPGYRO full orbit model. The preliminary calculation results of each model have shown differences in the perpendicular average velocity of impurity ions at the top region of the SOL. The mechanism which leads to the difference has been discussed. (© 2015 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)
PubDate: 2016-07-04T04:40:09.631922-05:
DOI: 10.1002/ctpp.201610068
- Authors: S. Yamoto; Y. Homma, H. Inoue, Y. Sawada, K. Hoshino, A. Hatayama, X. Bonnin, D. Coster, R. Schneider
- Simulation Analysis of Carbon Deposition Profile in the Closed Helical
Divertor Configuration in the Large Helical Device- Authors: M. Shoji; G. Kawamura, H. Tanaka, I. Watanabe, M. Kobayashi, M. Tokitani, S. Masuzaki, LHD Experiment Group
First page: 651
Abstract: Recent long pulse plasma discharges in the Large Helical Device have been interrupted by radiation collapse induced by large amounts of dusts released from a closed divertor region. Traces of exfoliation of carbon‐rich mixed material deposition layers were found in the divertor region after the experimental campaigns, indicating that the thickly accumulated deposited layers were exfoliated and broken into dusts. Simulation of the density profiles of carbon deposition in the divertor region by EMC3‐EIRENE reasonably explains the observation of the exfoliation of the deposition layers. It shows that the positions of the exfoliated layers correspond to the area where the density of the carbon deposition is high. It also proves that change of the configuration of target plates in the closed divertor is effective for controlling the accumulation of the carbon deposition, and carbon atoms released from the new target plates does not influence the main plasma confinement region. (© 2015 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)
PubDate: 2016-07-04T04:40:06.817162-05:
DOI: 10.1002/ctpp.201610059
- Authors: M. Shoji; G. Kawamura, H. Tanaka, I. Watanabe, M. Kobayashi, M. Tokitani, S. Masuzaki, LHD Experiment Group
- Photon Trapping Effects in DEMO Divertor Plasma
- Authors: K. Hoshino; K. Sawada, R. Idei, S. Tokunaga, N. Asakura, K. Shimizu, N. Ohno
First page: 657
Abstract: In the DEMO divertor, the neutral density becomes high to produce the full detachment and therefore the photon trapping can become important. In this paper, effects of the photon trapping on the DEMO divertor plasma has been studied. The pre‐evaluation of the photon trapping effects on the fixed background plasma profile was carried out by using an iterative self‐consistent collisional radiative model. The recombining plasma near the inner target and the private region changed to the ionizing plasma by the photon‐excitation. Based on the preevaluation result, the database of the effective ionization rate coefficient including the photon transport inside a 2 mm sphere. Advantage of the 2 mm sphere approximation is that the extra calculation cost is not necessary. Iterative calculation of the SONIC including the photon trapping effects was carried out. While the electron density increased and the neutral density decreased in the wide region, the electron density decreases close to the inner strike point. This may be due to decrease in the ionization rate by decrease in the neutral density. (© 2016 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)
PubDate: 2016-07-12T05:40:18.011422-05:
DOI: 10.1002/ctpp.201610039
- Authors: K. Hoshino; K. Sawada, R. Idei, S. Tokunaga, N. Asakura, K. Shimizu, N. Ohno
- Hybrid Formulation of Radiation Transport in Optically Thick Divertor
Plasmas- Authors: J. Rosato; Y. Marandet, H. Bufferand, D. Reiter, R. Stamm
First page: 663
Abstract: Kinetic Monte Carlo simulations of coupled atom‐radiation transport in optically thick divertor plasmas can be computationally very demanding, in particular in ITER relevant conditions or even larger devices, e.g. for power plant divertor studies. At high (∼ 1015 cm–3) atomic densities, it can be shown that sufficiently large divertors behave in certain areas like a black body near the first resonance line of hydrogen (Lyman α). This suggests that, at least in part, the use of continuum model (radiation hydrodynamics) can be sufficiently accurate, while being less time consuming. In this work, we report on the development of a hybrid model devoted to switch automatically between a kinetic and a continuum description according to the plasma conditions. Calculations of the photo‐excitation rate in a homogeneous slab are performed as an illustration. The outlined hybrid concept might be also applicable to neutral atom transport, due to mathematical analogy of transport equations for neutrals and radiation.
PubDate: 2016-07-13T10:20:06.908829-05:
DOI: 10.1002/ctpp.201610011
- Authors: J. Rosato; Y. Marandet, H. Bufferand, D. Reiter, R. Stamm
- A Model of Self‐similar Radiative Transfer in Resonance Lines for
Testing the Edge Plasma Codes- Authors: P. A. Sdvizhenskii; S. I. Krasheninnikov, A. B. Kukushkin
First page: 669
Abstract: Possibilities of obtaining the analytical solutions of the one‐dimensional equation of radiative transfer in resonance atomic/ionic lines in the Biberman‐Holstein model for testing the edge plasma codes are investigated. It is shown that for some types of similarity of spatial profiles of three characteristics, namely, background plasma density, line shape width and non‐radiation source of atomic excitation, the profile of excited atoms density may be described analytically in terms of the similarity of the above‐mentioned profiles. (© 2015 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)
PubDate: 2016-07-04T04:40:07.910788-05:
DOI: 10.1002/ctpp.201610058
- Authors: P. A. Sdvizhenskii; S. I. Krasheninnikov, A. B. Kukushkin
- Effects of Multi‐Species Ions on Sheath and Presheath in a Magnetic
Field Decreasing toward a Wall- Authors: A. Fukano; S. Nishioka, I. Goto, A. Hatayama
First page: 675
Abstract: In the divertor region, the plasma contains deuterium (D) ion and tritium (T) ion in an open magnetic field. The effect of D ion and T ion on the electric potential near the wall in plasma with magnetic field decreasing toward a wall is investigated analytically. The distribution of the electric potential is obtained by plasma‐sheath equation, where D ion and T ion are considered. The potential distribution depends on the parameters such as the profile of the magnetic field, the temperatures of D ion and T ion, and the amount rate of T ion to D ion. The particle distribution also is shown. (© 2016 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)
PubDate: 2016-06-24T04:30:09.897332-05:
DOI: 10.1002/ctpp.201610037
- Authors: A. Fukano; S. Nishioka, I. Goto, A. Hatayama
- Effect of PSI on Momentum Input to Plasma‐Facing Material Surfaces
- Authors: S. Takamura; Y. Uesugi, T. Kuwabara
First page: 681
Abstract: A comprehensive formulation for momentum input to plasma‐facing material surfaces is newly tried, including a momentum increase due to ion reflection on surfaces with a new definition of momentum reflection coefficient, a momentum due to reaction of evaporated atoms, an increase due to physical sputtering, and the effect of plasma flow. Experimentally observed momentum input is evaluated with a simple formula as well as a mechanical cantilever model. In addition, the measurement with dial tension gauge in off‐plasma condition gives a total force, the value of which is compared with the above two estimations. (© 2015 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)
PubDate: 2016-07-04T04:40:06.451093-05:
DOI: 10.1002/ctpp.201610061
- Authors: S. Takamura; Y. Uesugi, T. Kuwabara
- Electrostatic Characteristic of Spherical Dust on PFM in Sheath Field
- Authors: Y. Tomita; G. J. Niu, G. N. Luo
First page: 687
Abstract: The electrostatic characteristics of the spherical dust particle on the PFM were investigated analytically by introducing the bipolar coordinate in the plasma sheath field. The smaller dust compared to the Debye length λD is applicable to this model. It is clarified that the effect of finite size of the dust becomes remarkable at the larger dust than 0.2, Rd/λD > 0.2, where the change of the e‐folding length amounts to 10% of λD. The surface charge density on the dust is disturbed on the surface facing to the plasma due to the larger dust. (© 2016 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)
PubDate: 2016-07-12T05:50:06.531158-05:
DOI: 10.1002/ctpp.201610065
- Authors: Y. Tomita; G. J. Niu, G. N. Luo
- Visualization of Dust Particle Data with Plasma Simulation Results Using
Virtual‐Reality System- Authors: H. Ohtani; M. Shoji, N. Ohno, Y. Suzuki, S. Ishiguro, A. Kageyama, Y. Tamura
First page: 692
Abstract: We develop the environment in which the time‐sequence experimental data of three‐dimensional dust particle positions observed by the stereoscopic fast framing cameras is visualized in one virtual‐reality space with the results of the equilibrium plasma simulation in the vessel device data integrally using a virtual‐reality system. Since their three‐dimensional relative positions can be displayed, it is possible to analyze the three‐dimensional relationship between the dust trajectories, the magnetic‐field structure, and the devices in the vessel. It is confirmed that the observed dust particles are distributed in the periphery region. Most dusts move along the magnetic‐field line, but some particles move radially across the lines with sharply curved trajectories. The transport direction varies from place to place. This system promoted the viewers' understanding. (© 2015 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)
PubDate: 2016-07-04T04:40:08.19699-05:0
DOI: 10.1002/ctpp.201610054
- Authors: H. Ohtani; M. Shoji, N. Ohno, Y. Suzuki, S. Ishiguro, A. Kageyama, Y. Tamura
- Kinetic Modelling of the Plasma Recombination
- Authors: D. Tskhakaya
First page: 698
Abstract: Implementation of self‐consistent model of plasma recombination into the BIT1 PIC code and the simulation of detached SOL plasma are described. Our simulations indicate that in a strongly recombining plasma edge the sheath properties do not change qualitatively. The most affected parameter is the sheath heat transmission coefficient, which can increase by order of magnitude.
PubDate: 2016-07-18T07:30:03.784297-05:
DOI: 10.1002/ctpp.201611004
- Authors: D. Tskhakaya
- Simulations of Tungsten Re‐deposition Using a Particle‐In‐Cell Code
with Non‐uniform Super Particle Sizes- Authors: K. Ibano; S. Togo, T. L. Lang, Y. Ogawa, H.T. Lee, Y. Ueda, T. Takizuka
First page: 705
Abstract: Tungsten (W), a promising candidate as divertor plasma facing material in magnetic fusion devices, is anticipated to promptly redeposit when sputtered or evaporated from surface due to its small ionization energy and long gyro radius. Using an artificial factor for the reaction cross sections, effects of ionization lengths to the re‐deposition rate was studied by a newly developed particle‐in‐cell code. Treating numbers of particles in a super particle, electric charge, and mass as particle variables in the code, a special scheme for ionization and recombination was developed and used for the calculation. Simulations on W test particles with imaginary properties (neglecting the electric force) revealed the effects of ejection angles. Simulations with secondary electrons from surface showed that the sheath potential is weakened and the re‐deposition rate becomes small. It was found that the multi‐ionization as well as the ionization mean‐free‐path influences the re‐deposition rate in both simulations. (© 2016 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)
PubDate: 2016-06-24T04:30:06.540127-05:
DOI: 10.1002/ctpp.201610040
- Authors: K. Ibano; S. Togo, T. L. Lang, Y. Ogawa, H.T. Lee, Y. Ueda, T. Takizuka
- The Role of “Momentum Removal” in Divertor Detachment
- Authors: A. S. Kukushkin; H. D. Pacher
First page: 711
Abstract: The role of “momentum removal” (the drag force on the plasma ion flow) in divertor detachment is considered and analysed in detail. This analysis of the 2D modelling results shows that the drag force cannot reduce the power and particle flux to the target directly. However, it is essential for creating the conditions for efficient radiation and volumetric plasma recombination, which in turn do the job. (© 2015 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)
PubDate: 2016-07-04T04:30:04.850304-05:
DOI: 10.1002/ctpp.201610048
- Authors: A. S. Kukushkin; H. D. Pacher
- Detailed Analysis of Plasma Resistivity in Detached Recombining Plasmas
- Authors: H. Nishikata; Y. Hayashi, N. Ohno, S. Kajita, T. Kuwabara
First page: 717
Abstract: We have investigated plasma resistivity in a detached recombining plasma (DRP) by using a modified double probe theory, in which the influence of the plasma resistivity on the current‐voltage characteristics of the double probe, is taken into account. The plasma resistivity in DRP, obtained by comparing between the modified double probe theory and experimental results in the linear plasma device NAGDIS‐II, is much larger than both Spitzer resistivity and plasma resistivity caused by electron‐neutral collision including Rydberg atom's effect. (© 2015 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)
PubDate: 2016-07-04T04:40:09.352014-05:
DOI: 10.1002/ctpp.201610053
- Authors: H. Nishikata; Y. Hayashi, N. Ohno, S. Kajita, T. Kuwabara
- Statistical Analysis of Particle Flux Flowing into the End‐Target in
between Attached and Detached States in the Linear Divertor Plasma
Simulator NAGDIS‐II- Authors: H. Tanaka; N. Ohno, T. Onda, K. Takeyama, S. Kajita, T. Kuwabara, Y. Tsuji
First page: 723
Abstract: We have investigated the particle flux flowing into the axisymmetric end‐target in the transient state from attached to detached divertor conditions in the linear plasma device NAGDIS‐II. In the transient state, a dramatic decrease of the mean particle flux and a large‐amplitude fluctuation with negative and positive spikes were observed. We have analyzed the fluctuation with a newly suggested analysis technique: pre‐multiplied cubic spectrum with the wavelet transform. Analysis result indicates that these spikes consist of a few kilohertz components. (© 2015 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)
PubDate: 2016-07-04T04:40:04.91549-05:0
DOI: 10.1002/ctpp.201610064
- Authors: H. Tanaka; N. Ohno, T. Onda, K. Takeyama, S. Kajita, T. Kuwabara, Y. Tsuji
- Simulation Study of Detached Plasmas by Using One‐Dimensional
SOL‐Divertor Fluid Code with Virtual Divertor Model- Authors: S. Togo; T. Takizuka, M. Nakamura, K. Hoshino, K. Ibano, T. L. Lang, Y. Ogawa
First page: 729
Abstract: The detached plasmas due to the volume recombination are studied by using one‐dimensional (1D) scrape‐off‐layer‐divertor (SOL‐DIV) plasma fluid code with virtual divertor (VD) model. By introducing the anisotropic ion temperature, the parallel momentum transport equation becomes the first‐order differential and the Mach number at the sheath entrance is determined self‐consistently by the upstream condition. The total particle flux at the divertor plates and the flux amplification factors are shown as functions of the plasma density at the stagnation point and the dependence of these parameters on the heat flux from the core plasma, radial width of the flux tube in the divetor region and the strength of the impurity radiation is investigated. (© 2015 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)
PubDate: 2016-07-04T04:40:04.205197-05:
DOI: 10.1002/ctpp.201610067
- Authors: S. Togo; T. Takizuka, M. Nakamura, K. Hoshino, K. Ibano, T. L. Lang, Y. Ogawa
- Ginzburg‐Landau Model in a Finite Shear‐Layer and Onset of Transport
Barrier Nonlinear Oscillations: A Paradigm for TypeIII ELMs- Authors: M. Leconte; Y.M. Jeon, G. S. Yun
First page: 736
Abstract: We study a Reaction‐Diffusion model describing the nonlinear oscillations of a transport barrier in a finite shear‐layer (width dE ≪ a), where a is the plasma minor radius, based on a 1D reduced model derived to explain nonlinear barrier oscillations in 3D turbulence simulations [P. Beyer, S. Benkadda, G. Fuhr‐Chaudier et al., Phys. Rev. Lett. 94, 105001 (2005)]. We show that this single nonlinear equation encompasses most of the physics of these barrier relaxations. The nonlinear oscillations have common characteristics with type‐III edge localized modes (ELMs), such as a repetition frequency which decreases with increasing power. In addition to the flow shear, the shear‐layer width is also shown to control the nature of the oscillations. (© 2016 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)
PubDate: 2016-06-24T04:40:12.800868-05:
DOI: 10.1002/ctpp.201610050
- Authors: M. Leconte; Y.M. Jeon, G. S. Yun
- Modeling of Small‐sized ELMs in Detached Divertor Plasmas
- Authors: A. Yu. Pigarov; S.I. Krasheninnikov, E. M. Hollmann, T. D. Rognlien
First page: 742
Abstract: Self‐consistent modeling of edge plasma transport and deuterium wall inventory including multiple transient events was performed with the multi‐physics 2‐D transport code UEDGE‐MB‐W. In agreement with experimental data trends on DIII‐D, the modeling results show that relatively small‐sized and frequent type‐I Edge Localized Mode (ELM) events, which are typical for high‐power H‐mode discharges with strong deuterium gas‐puff fueling on this tokamak, are not “burning through” the formed detached plasma in the inner and, with further increase in fueling, in the outer divertor. In the latter case, the divertors are filled by sub‐eV, high‐density, strongly‐recombining and highly impurity contaminated plasma. Time‐dependent experimental data supporting the view that volumetric plasma recombination is enhanced during small‐sized ELMs penetrating into the detached plasma are discussed and the confirmatory results from UEDGE‐MB‐W modeling are presented. (© 2016 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)
PubDate: 2016-06-24T04:40:18.841279-05:
DOI: 10.1002/ctpp.201610055
- Authors: A. Yu. Pigarov; S.I. Krasheninnikov, E. M. Hollmann, T. D. Rognlien
- Core‐SOL Modelling of Neon Seeded JET Discharges with the
ITER‐like Wall- Authors: G. Telesca; I. Ivanova-Stanik, R. Zagörski, S. Brezinsek, A. Czarnecka, P. Drewelow, C. Giroud, A. Huber, S. Wiesen, JET EFDA contributors
First page: 748
Abstract: Five ELMy H‐mode Ne seeded JET pulses have been simulated with the self‐consistent core‐SOL model COREDIV. In this five pulse series only the Ne seeding rate was changed shot by shot, allowing a thorough study of the effect of Ne seeding on the total radiated power and of its distribution between core and SOL tobe made. The increase in the simulations of the Ne seeding rate level above that achieved in experiments shows saturation of the total radiated power at a relatively low radiated‐heating power ratio (frad = 0.60) and a further increase of the ratio of SOL to core radiation, in agreement with the reduction of W release at high Ne seeding level. In spite of the uncertainties caused by the simplified SOL model of COREDIV (neutral model, absence of ELMs and slab model for the SOL), the increase of the perpendicular transport in the SOL with increasing Ne seeding rate, which allows to reproduce numerically the experimental distribution core‐SOL of the radiated power, appears to be of general applicability.
PubDate: 2016-07-13T10:20:05.443647-05:
DOI: 10.1002/ctpp.201610013
- Authors: G. Telesca; I. Ivanova-Stanik, R. Zagörski, S. Brezinsek, A. Czarnecka, P. Drewelow, C. Giroud, A. Huber, S. Wiesen, JET EFDA contributors
- Effect of PFC Recycling Conditions on JET Pedestal Density
- Authors: S. Wiesen; S. Brezinsek, D. Harting, T. Dittmar, E. de la Luna, D. Matveev, K. Schmid, JET contributors
First page: 754
Abstract: There is experimental evidence that the pedestal dynamics in type‐I ELMy H‐mode discharges is significantly affected by a change in the recycling conditions at the tungsten plasma‐facing components (W‐PFCs) after an ELM event. The integrated code JINTRAC has been employed to assess the impact of recycling conditions during type‐I ELMs in JET ITER‐like wall H‐mode discharges. By employing a heuristic approach, a model to mimic the physical processes leading to formation and release (i.e. outgassing) of finite near‐surface fuel reservoirs in W‐PFCs has been implemented into the EDGE2D‐EIRENE plasma‐wall interaction code being part of JINTRAC. As main result it is shown, that a delay in the density pedestal build‐up after an ELM event can be provoked by reduced recycling induced by depleted W‐PFC particle near‐surface reservoirs. However the pedestal temperature evolution is barely affected by the change in recycling parameters suggesting that the presented model is incomplete.
PubDate: 2016-07-13T10:20:11.290013-05:
DOI: 10.1002/ctpp.201610015
- Authors: S. Wiesen; S. Brezinsek, D. Harting, T. Dittmar, E. de la Luna, D. Matveev, K. Schmid, JET contributors
- COREDIV and SOLPS Numerical Simulations of the Nitrogen Seeded JET ILW
L‐mode Discharges- Authors: I. Ivanova-Stanik; L. Aho-Mantila, M. Wischmeier, R. Zagörski, JET contributors
First page: 760
Abstract: In this paper we present the comparison of simulations with the numerical codes COREDIV and SOLPS5.0 for JET L‐mode discharges with ITER like wall (ILW). The simulations have been performed for L‐mode shots with and without nitrogen seeding (#82291 ‐ 9) which are characterised by relatively low auxiliary heating power (PNBI = 1.1 MW) and low electron density (ne = 2.35 × 1019 m–3). Comparisons are made to the experimental measurements (e.g. radiation levels, plasma profiles) and the differences between the results from the two codes (e.g. temperature and density profiles at the outer divertor plate) are shown and discussed.
PubDate: 2016-07-13T10:20:10.591547-05:
DOI: 10.1002/ctpp.201610006
- Authors: I. Ivanova-Stanik; L. Aho-Mantila, M. Wischmeier, R. Zagörski, JET contributors
- Modelling of the JET DT Experiments in Carbon and ITER‐like Wall
Configurations- Authors: R. Zagörski; W. Stȩpniewski, I. Ivanova-Stanik, S. Brezinsek, JET contributors
First page: 766
Abstract: In this paper numerical simulations with the self‐consistent COREDIV code of the planned JET DT experiments have been performed. First, record shot from the 1997 experiments was simulated and good agreement with experimental data has been found. Direct extrapolation of the carbon wall results to the new ILW configuration (discharge parameters as for the shot #42746) shows very good core plasma performance with even higher fusion power but with too large power to the divertor. However, with the neon seeding the heat load and plate temperatures can be efficiently reduced keeping good the plasma performance. Investigations have been done also for the planned DT operation scenario based on a conventional ELMy H‐mode at high plasma current and magnetic field. Simulations for the reference ELMy H‐mode shot #87412 show good agreement with the experimental data but the direct extrapolation of the DD results to deuterium‐tritium operation shows relatively poor performance in terms of the achieved fusion power. The situation improves, if the highest heating power is assumed (41 MW) and fusion powers in the excess of 12 MW can be achieved. All the high performance shots require the heat load control by neon seeding which shows rather beneficial effect on the plasma performance allowing for relatively wide operational window in terms of the amount of the allowed neon influx.
PubDate: 2016-07-13T10:20:09.861657-05:
DOI: 10.1002/ctpp.201610007
- Authors: R. Zagörski; W. Stȩpniewski, I. Ivanova-Stanik, S. Brezinsek, JET contributors
- Impurity Seeding in ASDEX Upgrade Tokamak Modeled by COREDIV Code
- Authors: K. Gała̧zka; I. Ivanova-Stanik, M. Bernert, A. Czarnecka, A. Kallenbach, R. Zagörski, ASDEX Upgrade Team
First page: 772
Abstract: The self‐consistent COREDIV code is used to simulate discharges in a tokamak plasma, especially the influence of impurities during nitrogen and argon seeding on the key plasma parameters. The calculations are performed with and without taking into account the W prompt redeposition in the divertor area and are compared to the experimental results acquired on ASDEX Upgrade tokamak (shots #29254 and #29257).For both impurities the modeling shows a better agreement with the experiment in the case without prompt redeposition. It is attributed to higher average tungsten concentration, which on the other hand seriously exceeds the experimental value. By turning the prompt redeposition process on, the W concentration is lowered, what, in turn, results in underestimation of the radiative power losses. By analyzing the influence of the transport coefficients on the radiative power loss and average W concentration it is concluded that the way to compromise the opposing tendencies is to include the edge‐localized mode flushing mechanism into the code, which dominates the experimental particle and energy balance. Also performing the calculations with both anomalous and neoclassical diffusion transport mechanisms included is suggested.
PubDate: 2016-07-13T10:20:09.077102-05:
DOI: 10.1002/ctpp.201610008
- Authors: K. Gała̧zka; I. Ivanova-Stanik, M. Bernert, A. Czarnecka, A. Kallenbach, R. Zagörski, ASDEX Upgrade Team
- Simulation of Radiative Divertor Plasmas by Ar Seeding with the Full
W‐Wall in JT‐60SA- Authors: H. Kawashima; K. Shimizu, K. Hoshino, T. Nakano, N. Asakura
First page: 778
Abstract: Radiative divertor plasmas for JT‐60SA with a full tungsten (W) wall, which is one of options in future, have been simulated with a SOL/divertor integrated code, SONIC. A conventional modified‐coronal radiation (MCR) model with a finite confinement time is used for both Ar and W for the purpose of wide‐range parameter surveys for the divertor plasma to obtain the required conditions (qt ≤ 10 MW/m2, nSepe–mid = 3∼8×1019 m–3,Prad< ∼ 30 MW), saving the calculation time. At low W density ratio (nW /ni = 1×10–5), due to low radiative power from W ions, Ar density ratio (nAr/ni ≥ 1.0× 10–3) and a strong gas puff (Γp ≥ 3.0×1022 s–1) are inevitable to suppress the divertor heat flux down to 10 MW/m2. Increasing nW/ni to 1×10–3 in the divertor region, the divertor heat load becomes low and the operative regions are expanded. While, the W production shall be suppressed since the W radiation is increased with replacement of Ar radiation and the particle recycling decreased. A Monte‐Carlo module (IMPMC) implemented in SONIC for Ar seeding reveals that the spatial distribution of Ar ions is predominantly determined by shell structures of the Ar ions. The consistency between IMPMC and MCR calculations is demonstrated for the averaged n Ar/n i ratio, the electron density and temperature profiles on the divertor target and typical parameter such as the divertor heat load. It shows that the detailed analysis with IMPMC model can be speedily obtained, using a steady state solution obtained by MCR model as an initial state. (© 2016 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)
PubDate: 2016-06-24T04:40:10.371915-05:
DOI: 10.1002/ctpp.201610045
- Authors: H. Kawashima; K. Shimizu, K. Hoshino, T. Nakano, N. Asakura
- Numerical Simulation Study of Plasma Flow in the GAMMA10/PDX
End‐cell Using a Fluid Code- Authors: H. Takeda; Y. Nakashima, A. Hatayama, M. S. Islam, K. Ichimura, M. M. Islam, K. Shimuzu, K. Fukui, M. Ohuchi, M. Sakamoto, T. Imai
First page: 784
Abstract: In order to understand the effect of impurity injection into divertor plasma, a multi‐fluid code was developed and applied to the region from west plug/barrier‐cell to end cell in GAMMA 10/PDX. In this research, Ar is applied as impurity neutrals. To investigate the effect of impurity injection on plasma parameters, interactions between Ar impurity and hydrogen bulk plasma were included in this model. Under the condition of injecting neutral hydrogen and impurity Ar, the charge exchange loss between hydrogen ion and neutrals reaches to about 95% and the radiation loss for electron is enhanced to about 45%. Therefore, hydrogen ion and electron temperatures decrease from about 55 eV to 4 eV and 15 eV to 4 eV, respectively. The heat load on the target plate is also reduced by nearly 80% compared with the condition of low injection rate of neutral hydrogen and impurity. It is also observed that the particle flux increases with neutral density and then shows a tendency to saturate in the higher neutral density (> 2.1 × 1018 m–3). (© 2016 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)
PubDate: 2016-07-12T05:50:07.277041-05:
DOI: 10.1002/ctpp.201610062
- Authors: H. Takeda; Y. Nakashima, A. Hatayama, M. S. Islam, K. Ichimura, M. M. Islam, K. Shimuzu, K. Fukui, M. Ohuchi, M. Sakamoto, T. Imai
- Reduced Physics Models in SOLPS for Reactor Scoping Studies
- Authors: D. P. Coster
First page: 790
Abstract: Heat exhaust is a challenge for ITER and becomes even more of an issue for devices beyond ITER. The main reason for this is that the power produced in the core scales as R3 while relying on standard exhaust physics results in the heat exhaust scaling as R1 (R is the major radius). ITER has used SOLPS (B2‐EIRENE) to design the ITER divertor, as well as to provide a database that supports the calculations of the ITER operational parameter space. The typical run time for such SOLPS runs is of the order 3 months (for D+C+He using EIRENE to treat the neutrals kinetically with an extensive choice of atomic and molecular physics). Future devices will be expected to radiate much of the power before it crosses the separatrix, and this requires treating extrinsic impurities such as Ne, Ar, Kr and Xe — the large number of charge states puts additional pressure on SOLPS, further slowing down the code.For design work of future machines, fast models have been implemented in system codes but these are usually unavoidably restricted in the included physics. As a bridge between system studies and detailed SOLPS runs, SOLPS offers a number of possibilities to speed up the code considerably at the cost of reducing the fidelity of the physics. By employing a fluid neutral model, aggressive bundling of the charge state of impurities, and reducing the size of the grids used, the run time for one second of physics time (which is often enough for the divertor to come to a steady state) can be reduced to approximately one day. This work looks at the impact of these trade‐offs in the physics by comparing key parameters for different simulation assumptions. (© 2016 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)
PubDate: 2016-06-24T04:30:16.689067-05:
DOI: 10.1002/ctpp.201610035
- Authors: D. P. Coster
- Magnetic Field Models and their Application in Optimal Magnetic Divertor
Design- Authors: M. Blommaert; M. Baelmans, H. Heumann, Y. Marandet, H. Bufferand, N. R. Gauger, D. Reiter
First page: 796
Abstract: In recent automated design studies, optimal design methods were introduced to successfully reduce the often excessive heat loads that threaten the divertor target surface. To this end, divertor coils were controlled to improve the magnetic configuration. The divertor performance was then evaluated using a plasma edge transport code and a “vacuum approach” for magnetic field perturbations. Recent integration of a free boundary equilibrium (FBE) solver allows to assess the validity of the vacuum approach. It is found that the absence of plasma response currents significantly limits the accuracy of the vacuum approach. Therefore, the optimal magnetic divertor design procedure is extended to incorporate full FBE solutions. The novel procedure is applied to obtain first results for the new WEST (Tungsten Environment in Steady‐state Tokamak) divertor currently under construction in the Tore Supra tokamak at CEA (Commissariat à l'Energie Atomique, France). The sensitivities and the related divertor optimization paths are strongly affected by the extension of the magnetic model. (© 2016 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)
PubDate: 2016-07-12T05:40:10.405111-05:
DOI: 10.1002/ctpp.201610031
- Authors: M. Blommaert; M. Baelmans, H. Heumann, Y. Marandet, H. Bufferand, N. R. Gauger, D. Reiter
- TECXY Study of a Liquid Lithium Divertor for DEMO
- Authors: G. Pełka; P. Chmielewski, R. Zagórski, V. Pericoli-Ridolfini, B. Viola
First page: 802
Abstract: Divertor targets made out of a capillary porous system (CPS) filled with liquid lithium, have been proposed as an alternative to standard, solid state plates. In the current work we simulate the DEMO edge plasma in either a standard single‐null or snowflake divertor configuration. Our tool is the 2D code TECXY.Lithium ablated from the target plate surface and released into the plasma is shown here to partially screen the incoming heat flux. Lithium's moderate SOL radiation levels suggest additional seeding to be beneficial. Very high heat fluxes to the divertor need to be avoided, as intensive lithium evaporation might unacceptably pollute the plasma.
PubDate: 2016-07-13T10:20:06.111229-05:
DOI: 10.1002/ctpp.201610012
- Authors: G. Pełka; P. Chmielewski, R. Zagórski, V. Pericoli-Ridolfini, B. Viola