Modelling impurity control at JET

Simonini, R.; Taroni, A.; Keilhacker, M.; Radford, Gail; Spence, James; Vlases, G.; Watkins, M.L.; Weber, S.

doi:10.1016/s0022-3115(06)80062-0

Cited by 35 publications

(37 citation statements)

References 2 publications

Supporting

Mentioning

Contrasting

Order By: Relevance

“…Presently, the favoured solutions to these problems are based on the extrapolation of the radiative divertor regimes observed in experiments, where volumetric losses reduce the power and ion¯ux to the target achieving the detached divertor regime. Two-dimensional plasma¯uid codes coupled to Monte-Carlo or diusive codes for neutral species such as B2-EIRENE [2,3], EDGE2D-NIMBUS [4,5] and UEDGE [6,7] have been developed to perform realistic calculations of the divertor and scrape-o layer (SOL) plasma parameters. These codes are routinely used to optimise the design of divertors for existing and next step devices.…”

Section: Introductionmentioning

confidence: 99%

Comparison of B2-EIRENE calculations with multi-machine experimental measurements

Loarte

Kukushkin²,

Pacher

et al. 1999

Journal of Nuclear Materials

View full text Add to dashboard Cite

The B2-EIRENE code has been compared with plasma edge measurements from Alcator C-MOD, DIII-D, JET, JT-60U in L-mode and H-mode regimes. Similarities in L-mode discharges between JET and JT-60U are well reproduced by the code as well as the sudden approach to detachment seen in JT-60U when a 1% chemical sputtering yield for carbon is assumed. Comparisons between similar L-mode regimes in JET, JT-60U and Alcator C-MOD reveal the larger importance of anomalous energy convection in Alcator C-MOD with regard to larger tokamaks. H-mode discharges of DIII-D and JET can be satisfactorily modelled with similar values of the transport coecients (approximately an order magnitude smaller than in L-mode) when the power¯ow into the scrape-o layer (SOL) is dominated by the ion channel, as previous EDGE2D simulations for JET suggested. Ó

show abstract

Section: Introductionmentioning

confidence: 99%

Comparison of B2-EIRENE calculations with multi-machine experimental measurements

Loarte

Kukushkin²,

Pacher

et al. 1999

Journal of Nuclear Materials

View full text Add to dashboard Cite

show abstract

“…In most 2D multi-fluid divertor codes [2][3][4], impurity transport is solved simultaneously as a fluid species. On the other hand, some elaborate Monte Carlo codes have been developed for the interpretation of impurity behavior in the divertor, e.g., DIVIMP [17] and IMPMC [6].…”

Section: Impmc Codementioning

confidence: 99%

“…This code includes the modeling of the dissociation process of methane and the dynamics of dissociation products (CD + 3 , CD 3 , CD + 2 , etc). For simplicity, these complicated dissociation processes of CD 4 are not taken into account in the present calculations. Instead, we treat chemically sputtered CD 4 as a carbon C with a low energy (∼1 eV).…”

Section: Analysis Of X-point Marfe In the Jt-60umentioning

confidence: 99%

See 1 more Smart Citation

Development of Integrated SOL/Divertor Code and Simulation Study in JAEA

Kawashima

Shimizu

Takizuka

et al. 2006

Plasma and Fusion Research

View full text Add to dashboard Cite

An integrated SOL/divertor code is being developed by the JAEA (Japan Atomic Energy Agency) for interpretation and prediction studies of the behavior of plasmas, neutrals, and impurities in the SOL/divertor region. A code system consists of the 2D fluid code for plasma (SOLDOR), the neutral Monte-Carlo code (NEUT2D), the impurity Monte-Carlo code (IMPMC), and the particle simulation code (PARASOL). The physical processes of neutrals and impurities are studied using the Monte Carlo (MC) code to accomplish highly accurate simulations. The so-called divertor code, SOLDOR/NEUT2D, has the following features: 1) a high-resolution oscillation-free scheme for solving fluid equations, 2) neutral transport calculation under the condition of fine meshes, 3) successful reduction of MC noise, and 4) optimization of the massive parallel computer. As a result, our code can obtain a steady state solution within 3 ∼ 4 hours even in the first run of a series of simulations, allowing the performance of an effective parameter survey. The simulation reproduces the X-point MARFE (multifaceted asymmetric radiation from edge) in the JT-60U. It is found that the chemically sputtered carbon at the dome causes radiation peaking near the X-point. The performance of divertor pumping in the JT-60U is evaluated based on particle balances. In regard to the divertor design of the next tokamak of JT-60U, the simulation indicates the dependencies of pumping efficiency on the divertor geometry and operational conditions. The efficiency is determined by the balance between the incident and back-flow fluxes into and from the exhaust chamber.

show abstract

“…A majority, such as B2 [11], UEDGE [12], EDGE2D [13], and SOLDOR [14] are fluid codes, which reduce the dimensionality by assuming a Maxwellian velocity distribution. However, the fluid model is not self-consistent; it requires closure parameters and boundary conditions inside the sheath region to be given by external models.…”

Section: Introductionmentioning

confidence: 99%

Kinetic Particle Simulation Study of Parallel Heat Transport in Scrape-off Layer Plasmas over a Wide Range of Collisionalities

Froese

Takizuka

Yagi

2010

Plasma and Fusion Research

View full text Add to dashboard Cite

Fluid models are not generally applicable to fusion edge plasmas without external provision of kinetic factors: closure parameters and boundary conditions inside the sheath region. We explain the PARASOL-1D simulation, a particle-in-cell code with a binary collision Monte-Carlo model, and use it to determine four kinetic factors commonly needed in fluid codes. These are the electron and ion heat flux limiting factors, α e and α i , the ion adiabatic index, γ A , and the electron and ion temperature anisotropy, T || /T ⊥ . We survey these factors over a wide range of collisionalities and find that, as predicted, the conductive heat flux is accurately described by the Spitzer-Härm expression in the collisional limit and asymptotes to a constant value in the collisionless limit. However, unique behavior occurs in the weakly collisional regime when the ratio of the mean free path to connection length is 0.1 < λ mfp /L || < 10, when the SOL is between the conduction-and sheath-limited regimes. We find that α e can peak, becoming larger than the collisionless limit, γ A is less than unity, and only the ions are anisotropic. The effects of electron energy radiation and Langevin heating are explored. Finally, the strong deviations of the energy distribution function from Maxwellian in the weakly collisional and collisionless regimes are explained.

show abstract

Modelling impurity control at JET

Cited by 35 publications

References 2 publications

Comparison of B2-EIRENE calculations with multi-machine experimental measurements

Comparison of B2-EIRENE calculations with multi-machine experimental measurements

Development of Integrated SOL/Divertor Code and Simulation Study in JAEA

Kinetic Particle Simulation Study of Parallel Heat Transport in Scrape-off Layer Plasmas over a Wide Range of Collisionalities

Contact Info

Product

Resources

About