Sticking and erosion at carbon-containing plasma-facing materials in fusion reactors

Controlled nuclear fusion is one of the most promising sources of energy for the future. Before this goal can be achieved, one must be able to control the enormous energy densities which are present in the core plasma in a fusion reactor. In order to be able to predict the evolution and thereby the lifetime of different plasma facing materials under reactor-relevant conditions, the interaction of atoms and molecules with plasma first wall surfaces have to be studied in detail. In this thesis, the fundamental sticking and erosion processes of carbon-based materials, the nature of hydrocarbon species released from plasma-facing surfaces, and the evolution of the components under cumulative bombardment by atoms and molecules have been investigated by means of molecular dynamics simulations using both analytic potentials and a semi-empirical tight-binding method. The sticking cross-section of CH3 radicals at unsaturated carbon sites at diamond (111) surfaces is observed to decrease with increasing angle of incidence, a dependence which can be described by a simple geometrical model. The simulations furthermore show the sticking cross-section of CH3 radicals to be strongly dependent on the local neighborhood of the unsaturated carbon site. The erosion of amorphous hydrogenated carbon surfaces by helium, neon, and argon ions in combination with hydrogen at energies ranging from 2 to 10 eV is studied using both non-cumulative and cumulative bombardment simulations. The results show no significant differences between sputtering yields obtained from bombardment simulations with different noble gas ions. The final simulation cells from the 5 and 10 eV ion bombardment simulations, however, show marked differences in surface morphology. In further simulations the behavior of amorphous hydrogenated carbon surfaces under bombardment with D^+, D^+2, and D^+3 ions in the energy range from 2 to 30 eV has been investigated. The total chemical sputtering yields indicate that molecular projectiles lead to larger sputtering yields than atomic projectiles. Finally, the effect of hydrogen ion bombardment of both crystalline and amorphous tungsten carbide surfaces is studied. Prolonged bombardment is found to lead to the formation of an amorphous tungsten carbide layer, regardless of the initial structure of the sample. In agreement with experiment, preferential sputtering of carbon is observed in both the cumulative and non-cumulative simulations

Computer simulation of H and He effects in fusion reactor materials

Fusion power is an appealing source of clean and abundant energy. The radiation resistance of reactor materials is one of the greatest obstacles on the path towards commercial fusion power. These materials are subject to a harsh radiation environment, and cannot fail mechanically or contaminate the fusion plasma. Moreover, for a power plant to be economically viable, the reactor materials must withstand long operation times, with little maintenance. The fusion reactor materials will contain hydrogen and helium, due to deposition from the plasma and nuclear reactions because of energetic neutron irradiation. The first wall divertor materials, carbon and tungsten in existing and planned test reactors, will be subject to intense bombardment of low energy deuterium and helium, which erodes and modifies the surface. All reactor materials, including the structural steel, will suffer irradiation of high energy neutrons, causing displacement cascade damage. Molecular dynamics simulation is a valuable tool for studying irradiation phenomena, such as surface bombardment and the onset of primary damage due to displacement cascades. The governing mechanisms are on the atomic level, and hence not easily studied experimentally. In order to model materials, interatomic potentials are needed to describe the interaction between the atoms. In this thesis, new interatomic potentials were developed for the tungsten-carbon-hydrogen system and for iron-helium and chromium-helium. Thus, the study of previously inaccessible systems was made possible, in particular the effect of H and He on radiation damage. The potentials were based on experimental and ab initio data from the literature, as well as density-functional theory calculations performed in this work. As a model for ferritic steel, iron-chromium with 10% Cr was studied. The difference between Fe and FeCr was shown to be negligible for threshold displacement energies. The properties of small He and He-vacancy clusters in Fe and FeCr were also investigated. The clusters were found to be more mobile and dissociate more rapidly than previously assumed, and the effect of Cr was small. The primary damage formed by displacement cascades was found to be heavily influenced by the presence of He, both in FeCr and W. Many important issues with fusion reactor materials remain poorly understood, and will require a huge effort by the international community. The development of potential models for new materials and the simulations performed in this thesis reveal many interesting features, but also serve as a platform for further studies.

Computer simulation of multi-elemental fusion reactor materials

Thermonuclear fusion is a sustainable energy solution, in which energy is produced using similar processes as in the sun. In this technology hydrogen isotopes are fused to gain energy and consequently to produce electricity. In a fusion reactor hydrogen isotopes are confined by magnetic fields as ionized gas, the plasma. Since the core plasma is millions of degrees hot, there are special needs for the plasma-facing materials. Moreover, in the plasma the fusion of hydrogen isotopes leads to the production of high energetic neutrons which sets demanding abilities for the structural materials of the reactor. This thesis investigates the irradiation response of materials to be used in future fusion reactors. Interactions of the plasma with the reactor wall leads to the removal of surface atoms, migration of them, and formation of co-deposited layers such as tungsten carbide. Sputtering of tungsten carbide and deuterium trapping in tungsten carbide was investigated in this thesis. As the second topic the primary interaction of the neutrons in the structural material steel was examined. As model materials for steel iron chromium and iron nickel were used. This study was performed theoretically by the means of computer simulations on the atomic level. In contrast to previous studies in the field, in which simulations were limited to pure elements, in this work more complex materials were used, i.e. they were multi-elemental including two or more atom species. The results of this thesis are in the microscale. One of the results is a catalogue of atom species, which were removed from tungsten carbide by the plasma. Another result is e.g. the atomic distributions of defects in iron chromium caused by the energetic neutrons. These microscopic results are used in data bases for multiscale modelling of fusion reactor materials, which has the aim to explain the macroscopic degradation in the materials. This thesis is therefore a relevant contribution to investigate the connection of microscopic and macroscopic radiation effects, which is one objective in fusion reactor materials research.

Simulation and experimental study on compositional evolution of Li-Co in LiCoO(2) thin films during sputter deposition

The compositional evolution in sputter deposited LiCoO(2) thin films is influenced by process parameters involved during deposition. The electrochemical performance of these films strongly depends on their microstructure, preferential orientation and stoichiometry. The transport process of sputtered Li and Co atoms from the LiCoO(2) target to the substrate, through Ar plasma in a planar magnetron configuration, was investigated based on the Monte Carlo technique. The effect of sputtering gas pressure and the substrate-target distance (d(st)) on Li/Co ratio, as well as, energy and angular distribution of sputtered atoms on the substrate were examined. Stable Li/Co ratios have been obtained at 5 Pa pressure and d(st) in the range 5 11 cm. The kinetic energy and incident angular distribution of Li and Co atoms reaching the substrate have been found to be dependent on sputtering pressure. Simulations were extended to predict compositional variations in films prepared at various process conditions. These results were compared with the composition of films determined experimentally using x-ray photoelectron spectroscopy (XPS). Li/Co ratio calculated using XPS was in moderate agreement with that of the simulated value. The measured film thickness followed the same trend as predicted by simulation. These studies are shown to be useful in understanding the complexities in multicomponent sputtering. (C) 2011 American Institute of Physics. doi:10.1063/1.3597829]

Simulation and experimental study on compositional evolution of Li-Co in LiCoO2 thin films during sputter deposition

The compositional evolution in sputter deposited LiCoO2 thin films is influenced by process parameters involved during deposition. The electrochemical performance of these films strongly depends on their microstructure, preferential orientation and stoichiometry. The transport process of sputtered Li and Co atoms from the LiCoO2 target to the substrate, through Ar plasma in a planar magnetron configuration, was investigated based on the Monte Carlo technique. The effect of sputtering gas pressure and the substrate-target distance (dst) on Li/Co ratio, as well as, energy and angular distribution of sputtered atoms on the substrate were examined. Stable Li/Co ratios have been obtained at 5 Pa pressure and dst in the range 5−11 cm. The kinetic energy and incident angular distribution of Li and Co atoms reaching the substrate have been found to be dependent on sputtering pressure. Simulations were extended to predict compositional variations in films prepared at various process conditions. These results were compared with the composition of films determined experimentally using x-ray photoelectron spectroscopy (XPS). Li/Co ratio calculated using XPS was in moderate agreement with that of the simulated value. The measured film thickness followed the same trend as predicted by simulation. These studies are shown to be useful in understanding the complexities in multicomponent sputtering.

MHD arcjet onset boundary of high specific impulse Plasma Thrusters

Magnetoplasmadynamic thrusters are known to enter a strongly unstable regime, calledas onset in the literature, under high specific impulse operation. This paper probes the early signs of onset in relatively moderate specific impulse operation by a single fluid plasma thruster simulation. The procedure involves solving the combined Maxwell’s-Navier-Stokes equation, with an onset criterion of radial current reaching close to zero values near the electrodes. Thruster parameters are varied starting from voltage potential, plasma temperature and cathodic radius. Onset curves are plotted which can provide important engine-specific information in order to understand the onset performance of the plasma thruster.

Atmospheric pressure plasma chemical vapor deposition reactor for 100 mm wafers, optimized for minimum contamination at low gas flow rates

Gas discharge plasmas used for thinfilm deposition by plasma-enhanced chemical vapor deposition (PECVD) must be devoid of contaminants, like dust or active species which disturb the intended chemical reaction. In atmospheric pressure plasma systems employing an inert gas, the main source of such contamination is the residual air inside the system. To enable the construction of an atmospheric pressure plasma (APP) system with minimal contamination, we have carried out fluid dynamic simulation of the APP chamber into which an inert gas is injected at different mass flow rates. On the basis of the simulation results, we have designed and built a simple, scaled APP system, which is capable of holding a 100 mm substrate wafer, so that the presence of air (contamination) in the APP chamber is minimized with as low a flow rate of argon as possible. This is examined systematically by examining optical emission from the plasma as a function of inert gas flow rate. It is found that optical emission from the plasma shows the presence of atmospheric air, if the inlet argon flow rate is lowered below 300 sccm. That there is minimal contamination of the APP reactor built here, was verified by conducting an atmospheric pressure PECVD process under acetylene flow, combined with argon flow at 100 sccm and 500 sccm. The deposition of a polymer coating is confirmed by infrared spectroscopy. X-ray photoelectron spectroscopy shows that the polymer coating contains only 5% of oxygen, which is comparable to the oxygen content in polymer deposits obtained in low-pressure PECVD systems. (C) 2015 AIP Publishing LLC.

Rossby Vortex Simulation on a Paraboloidal Coordinate System Using the Lattice Boltzmann Method

In this paper, we apply our compressible lattice Boltzmann model to a rotating parabolic coordinate system to simulate Rossby vortices emerging in a layer of shallow water flowing zonally in a rotating paraboloidal vessel. By introducing a scaling factor, nonuniform curvilinear mesh can be mapped to a flat uniform mesh and then normal lattice Boltzmann method works. Since the mass per unit area on the two-dimensional (2D) surface varies with the thickness of the water layer, the 2D flow seems to be "compressible" and our compressible model is applied. Simulation solutions meet with the experimental observations qualitatively. Based on this research, quantitative solutions and many natural phenomena simulations in planetary atmospheres, oceans, and magnetized plasma, such as the famous Jovian Giant Red Spot, the Galactic Spiral-vortex, the Gulf Stream, and the Kuroshio Current, etc,, can be expected.

Study of Gliding Arc Discharge Plasma

本文研究了滑动弧放电过程中电参数的变化，并对滑动弧等离子体中的非平衡度和各参数之间关系进行了讨论。应用了双通道电弧模型 ，对电弧在气流作用下的运动规律进行了数值模拟。模拟的结果有助于分析滑动弧非平衡等离子体的产生机理。 The elelctric parameters change during discharge is studied and the relationship between non-equilibrium degree and parameters is discussed for gliding arc discharge.Using two-channel model, the rules of arc moving due to effect of the airflow is simulated.The numerical simulation results can help analyzing the generation mechanism of gliding arc non-equilibrium plasma.

Modelling of Laminar Plasma Jet Impinging on a Flatplate with Approximate Box Relaxation Method

Approximate Box Relaxation method was used t'o simulate a plasma jet flow impinging on a flatplate at atmospheric pressure, to achieve a better understanding of the characteristics of plasma jet in materials surface treating. The flow fields under different conditions were simulated and analyzed. The distributions of temperature, velocity and pressure were obtained by modelling. Computed results indicate that this numerical method is suitable for simulation of the flow characteristics of plasma jet: and is helpful for understanding of the mechanism of the plasma-material processing.

Study on the ions' behavior in an electron cyclotron resonance plasma

The energy, velocity, angle distribution of ions in magnetoactive electron cyclotron resonance plasma have been studied with a two-dimension hybrid mode. The dependence of these distribution functions versus position and pressure are discussed. Our simulation results are in good agreement with many experimental measurements. (C) 1997 American Institute of Physics.

Behavior of ar plasma formed in a high-density plasma source - an ecr reactor

In order to develop the ultra-large scale integration(ULSI), low pressure and high density plasma apparatus are required for etching and deposit of thin films. To understand critical parameters such as the pressure, temperature, electrostatic potential and energy distribution of ions impacting on the wafer, it is necessary to understand how these parameters are influenced by the power input and neutral gas pressure. In the present work, a 2-D hybrid electron fluid-particle ion model has been developed to simulate one of the high density plasma sources-an Electron Cyclotron Resonance (ECR) plasma system with various pressures and power inputs in a non-uniform magnetic field. By means of numerical simulation, the energy distributions of argon ion impacting on the wafer are obtained and the plasma density, electron temperature and plasma electrostatic potential are plotted in 3-D. It is concluded that the plasma density depends mainly on both the power input and neutral gas pressure. However, the plasma potential and electron temperature can hardly be affected by the power input, they seem to be primarily dependent on the neutral gas pressure. The comparison shows that the simulation results are qualitatively in good agreement with the experiment measurements.

The induction plasma chemical reactor: Part I. Equilibrium model

A mathematical model is presented for the numerical simulation of the flow, temperature, and concentration fields in an rf plasma chemical reactor. The simulation is performed assuming chemical equilibrium. The extent of validity of this assumption is discussed. The system considered is the reaction of SiCl4 and NH3 for the production of Si3N4.

Plasma-Arc Technology For The Thermal Treatment Of Chemical Wastes

Plasma-arc technology was developed to dispose of chemical wastes from a chemical plant by the Institute of Mechanics, Chinese Academy of Sciences (CAS-IMECH). A pilot plant system with this technology was constructed to destroy two types of chemical wastes. The system included shredding, mixing, and feeding subsystems, a plasma-arc reactor of 150 kW, an off-gas burning subsystem, and a scrubbing subsystem. The additives (CaO, SiO2, and Fe) were added into the reactor to form vitrified slag and capture the hazardous elements. The molten slag was quickly quenched to form an amorphous glassy structure. A direct current (DC) experimental facility of 30kW with plasma-arc technology was also set up to study the pyrolysis process in the laboratory, and the experimental results showed the cooling speed is the most important factor for good vitrified structure of the slag. According to previous tests, the destruction and removal efficiency (DRE) for these chemical wastes was more than 99.999%, and the polychlorinated biphenyls (PCBs) concentration in the solid residues was in the range of 1.28 to 12.9mg/kg, which is far below the Chinese national emission limit for the hazardous wastes. A simplified electromagneto model for numerical simulation was developed to predict the temperature and velocity fields. This model can make satisfactory maximum temperature and velocity distributions in the arc region, as well as the results by the magneto hydrodynamic approach.

Plasma channeling by multiple short-pulse lasers

Channeling by a train of laser pulses into homogeneous and inhomogeneous plasmas is studied using particle-in-cell simulation. When the pulse duration and the interval between the successive pulses are appropriate, the laser pulse train can channel into the plasma deeper than a single long-pulse laser of similar peak intensity and total energy. The increased penetration distance can be attributed to the repeated actions of the ponderomotive force, the continuous between-pulse channel lengthening by the inertially evacuating ions, and the suppression of laser-driven plasma instabilities by the intermittent laser-energy cut-offs.