The use of an all-metal first wall in future magnetic fusion reactors equipped with a divertor may impose severe limitations on the capabilities of optical diagnostics in the main chamber because of ...the divertor stray light (DSL) produced by reflections of the intense light emitted in the divertor. Here, we introduce a synthetic H-alpha diagnostics to estimate the errors of solutions of the inverse problems aimed at recovering the neutral hydrogen parameters (density and isotope ratio) in the scrape-off layer (SOL) with allowance for (a) strong DSL on the observation chords in the main chamber, (b) substantial deviation of the neutral atom velocity distribution function from a Maxwellian in the SOL, and (c) the data from the direct observation of the divertor. The results of recovering the relative contributions of all three sources to the signal along an observation chord in the main chamber (namely, from the high-field-side and low-field-side SOL sections of the observation chord, and the DSL), together with the isotope ratios in the SOL, are presented for the flattop stage of Q = 10 inductive operation of ITER.
In ITER, as in any tokamak, the first wall and divertor plasma-facing components (PFC) must provide adequate protection of in-vessel structures, sufficient heat exhaust capability and be compatible ...with the requirements of plasma purity. These functions take on new significance in ITER, which will combine long pulse, high power operation with severe restrictions on permitted core impurity concentrations and which, in addition, will produce transient energy loads on a scale unattainable in today’s devices. The current ITER PFC design has now reached a rather mature stage following the 2007 ITER Design Review. This paper presents the key elements of the design, reviews the physics drivers, essentially thermal load specifications, which have defined the concept and discusses a selection of material and design issues.
The new SOLPS-ITER code package Wiesen, S.; Reiter, D.; Kotov, V. ...
Journal of nuclear materials,
08/2015, Letnik:
463
Journal Article
Recenzirano
The SOLPS-ITER package represents a renewed coupling of the up-to-date parallelized EIRENE Monte-Carlo code and the B2.5 edge plasma solver. Each part of the fluid-kinetic code system has been ...dismantled and analysed. For SOLPS-ITER, full backwards compatibility with previous SOLPS4.3 ITER reference simulations is critical. We present the results of the ongoing comprehensive benchmark activity and discuss how SOLPS-ITER can be tuned to reproduce ITER simulations performed with SOLPS4.3. In particular in SOLPS-ITER a quantitative “metric” to measure the level of code convergence has been implemented, with corresponding code diagnostics carried over from SOLPS4.3, allowing code run times to be accurately compared. Good agreement is found between the two codes in the tests so far. Remaining discrepancies in the solution in front of the divertor targets are likely related to different formulations of the equations and the different discretization schemes used by the codes.
► This is a review of the development of edge plasma modeling at ITER and of its interaction with the evolving divertor design. ► The SOLPS (B2-Eirene) code has been developed for the evaluation and ...the design of the ITER divertor for the last 15 years. ► Divertor modeling has developed into essential design tool synthesizing information from theoretical analysis, experimental studies, and engineering intuition. ► A number of examples given in the paper illustrate this process.
The paper presents a review of the development of edge plasma modeling at ITER and of its interaction with the evolving divertor design. The SOLPS (B2-Eirene) code has been developed for, and applied to, the evaluation and the design of the ITER divertor for the last 15 years. With respect to the physics and engineering design, divertor modeling had started as an evaluation tool and has developed into essential design tool synthesizing information from theoretical analysis, experimental studies, and engineering intuition. Examples given in the paper illustrate this process.
In the present work, a new method for obtaining silicon carbide of the cubic polytype 3C-SiC with silicon vacancies in a stable state is proposed theoretically and implemented experimentally. The ...idea of the method is that the silicon vacancies are first created by high-temperature annealing in a silicon substrate Si(111) doped with boron B, and only then is this silicon converted into 3C-SiC(111), due to a chemical reaction with carbon monoxide CO. A part of the silicon vacancies that have bypassed “chemical selection” during this transformation get into the SiC. As the process of SiC synthesis proceeds at temperatures of ~1350 °C, thermal fluctuations in the SiC force the carbon atom C adjacent to the vacancy to jump to its place. In this case, an almost flat cluster of four C atoms and an additional void right under it are formed. This stable state of the vacancy, by analogy with NV centers in diamond, is designated as a C4V center. The C4V centers in the grown 3C-SiC were detected experimentally by Raman spectroscopy and spectroscopic ellipsometry. Calculations performed by methods of density-functional theory have revealed that the C4V centers have a magnetic moment equal to the Bohr magneton μB and lead to spin polarization in the SiC if the concentration of C4V centers is sufficiently high.
In this work, silicon carbide layers containing silicon vacancies are grown by the Method of Coordinated Substitution of Atoms (MCSA). The main idea of this fundamentally new method is that silicon ...vacancies are first created in silicon, which is much simpler, and only then is silicon converted into silicon carbide by chemical reaction with carbon monoxide. The dielectric function of silicon carbide containing silicon vacancies, grown on both n- and p-type silicon substrates, is measured for the first time. The density functional method in the spin-polarized approximation is used to calculate the dielectric function of silicon carbide containing silicon vacancies. It is shown that the influence of the magnetic moment of vacancies on the dielectric function is decisive. Qualitative correspondence of the computational model to the obtained experimental data is demonstrated. It is discovered that silicon vacancies make silicon carbide much less transparent. It is shown that the imaginary part of the dielectric function is described as a sum of oscillatory peaks in the form of the Gaussian functions. Vacancies lead, as a rule, to one or two additional peaks. According to the amplitude and position of the additional peaks, it is possible to qualitatively estimate the concentration of vacancies and their charge.
Thin films of single-crystal silicon carbide of cubic polytype with a thickness of 40-100 nm, which were grown from the silicon substrate material by the method of coordinated substitution of atoms ...by a chemical reaction of silicon with carbon monoxide CO gas, have been studied by spectral ellipsometry in the photon energy range of 0.5-9.3 eV. It has been found that a thin intermediate layer with the dielectric constant corresponding to a semimetal is formed at the 3C-SiC(111)/Si(111) interface. The properties of this interface corresponding to the minimum energy have been calculated using quantum chemistry methods. It has turned out that silicon atoms from the substrate are attracted to the interface located on the side of the silicon carbide (SiC) film. The symmetry group of the entire system corresponds to
. The calculations have shown that Si atoms in silicon carbide at the interface, which are the most distant from the Si atoms of the substrate and do not form a chemical bond with them (there are only 12% of them), provide a sharp peak in the density of electronic states near the Fermi energy. As a result, the interface acquires semimetal properties that fully correspond to the ellipsometry data.
Plasma wall interaction (PWI) is important for the material choice in ITER and for the plasma scenarios compatible with material constraints. In this paper, different aspects of the PWI are assessed ...in their importance for the initial wall materials choice: CFC for the strike point tiles, W in the divertor and baffle and Be on the first wall. Further material options are addressed for comparison, such as W divertor/Be first wall and all-W or all-C. One main parameter in this evaluation is the particle flux to the main vessel wall. One detailed plasma scenario exists for a
Q
=
10 ITER discharge G. Federici et al., J. Nucl. Mater. 290–293 (2001) 260 which was taken as the basis of further erosion and tritium retention evaluations. As the assessment of steady state wall fluxes from a scaling of present fusion devices indicates that global wall fluxes may be a factor of 4
±
3 higher, this margin has been adopted as uncertainty of the scaling. With these wall and divertor fluxes, important PWI processes such as erosion and tritium accumulation have been evaluated: It was found that the steady state erosion is no problem for the lifetime of plasma-facing divertor components. Be wall erosion may pose a problem in case of a concentration of the wall fluxes to small wall areas. ELM erosion may drastically limit the PFC lifetime if ELMs are not mitigated to energies below 0.5
MJ. Dust generation is still a process which requires more attention. Conversion from gross or net erosion to dust and the assessment of dust on hot surfaces need to be investigated. For low-
Z materials the build-up of the tritium inventory is dominated by co-deposition with eroded wall atoms. For W, where erosion and tritium co-deposition are small, the implantation, diffusion and bulk trapping constitute the dominant retention processes. First extrapolations with models based on laboratory data show small contributions to the inventory. For later ITER phases and the extrapolation to DEMO additional tritium trapping sites due to neutron-irradiation damage need to be taken into account. Finally, the expected values for erosion and tritium retention are compared to the ITER administrative limits for the lifetime, dust and tritium inventory.
A possibility of up–down asymmetry of the profiles of the edge plasma parameters in a symmetric double-null divertor configuration in the case of impurity seeding is found by computational modeling. ...The physical mechanism of such a symmetry violation, related to radiation-condensation instability, is proposed. This effect can reduce the would-be advantage of the double-null configuration in the power load spreading over a larger area, which is vital for fusion reactor projects.
The effect of thermodynamic nonequilibrium of hydrogen isotope recycling, which consists in significant deviation of the velocity distribution function (VDF) of neutral hydrogen isotope atoms from ...the Maxwellian VDF, on the passive signal of the charge-exchange recombination spectroscopy (CXRS diagnostics) of edge plasma in tokamaks is analyzed using the ITER tokamak reactor as an example. The calculations of the VDF for atomic and molecular deuterium are carried out using the three-dimensional Monte Carlo code EIRENE as a part of the two-dimensional transport code SOLPS and the semi-analytical Ballistic Model BM1D2V, which is one-dimensional in coordinate and two-dimensional in projections of atomic velocity. The source of level population of the Be IV beryllium ion, the line radiation of which is used in the CXRS diagnostics, is calculated for the above two methods of calculating the VDF and the average atomic kinetic energy, as well as in the case of neglecting the thermodynamic nonequilibrium of recycling, when the temperature of atoms is equal to the temperature of the plasma ions. The calculations of the cross section for the charge-exchange reaction of beryllium ions Be V on deuterium atoms using the ARSENY code and other codes and the erosion kinetics of the beryllium first wall using the three-dimensional code ERO2.0 are used. The calculations of the level populations of the Be IV beryllium ion using the
nl
-KinRyd code, the luminosity of the ions, and the passive charge-exchange radiation intensity on the observation chord are carried out for the three indicated cases of the VDF of deuterium atoms. A significant dependence of the results of the conducted predictive modeling on the thermodynamic nonequilibrium of recycling is shown.