Two-dimensional simulation of inductive plasma sources with self-consistent power deposition

A time averaged two-dimensional fluid model including an electromagnetic module with self-consistent power deposition was developed to simulate the transport of a low pressure radio frequency inductively coupled plasma source. Comparsions with experiment and previous simulation results show, that the fluid model is feasible in a certain range of gas pressure. In addition, the effects of gas pressure and power input have been discussed.

Parameter-estimation of a spherical top inductively-coupled plasma source with global-model

Using spatially averaged global model, we succeed in obtaining some plasma parameters for a low pressure inductively coupled plasma source of our laboratory. As far as the global balance is concerned, the models can give reasonable results of the parameters, such as the global electron temperature and the ion impacting energy, etc. It is found that the ion flow is hardly affected by the neutral gas pressure. Finally, the magnetic effects are calculated by means of the method. The magnetic field can play an important role to increase plasma density and ion current.

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.

A note on the heat-transfer to nonspherical particles from plasma

Heat transfer from plasma to a nonspherical partical in the free-molecular regime is studied in the present paper under thin plasma sheath condition. Analytical expressions for the floating potential charge and heat fluxes of an ellipsoid particle of revolution are derived and curves are given for key parameters for arbitrary plasma flow direction. On the basis of these results, an equivalent sphere with the same surface area as the nonspherical particle is suggested to be used for calculating the total heat flux of nonspherical particle in engineering application with acceptable accuracy. Furthermore, the effects of particle rotation, which occurs in most aerosol systems, on the heat transfer are also discussed.

The positive floating potential of high-temperature particles in plasma

It has been predicted that the floating potential of particles in plasma may become positive when the particle surface temperature is high enough, but, to our knowledge, no positive floating potential has been obtained yet. In the present paper the floating potential theory of high-temperature particles in plasma is developed to cover the positive potential range for the first time, and a general approximate analytical formula for the positive floating potential with a thin plasma sheath and subsonic plasma flow is derived from the new model recently proposed by the authors. The results show that when the floating potential is positive, the net flux of charge incident on the particle approaches a constant similar to the 'electron saturation' phenomena in the case of the electric probes.

Numerical-simulation of ac plasma-arc thermodynamics

A mathematical model and approximate analysis for the energy distribution of an ac plasma arc with a moving boundary is developed. A simplified electrical conductivity function is assumed so that the dynamic behavior of the arc may be determined, independent of the gas type. The model leads to a reduced set of non-linear partial differential equations which governs the quasi-steady ac arc. This system is solved numerically and it is found that convection plays an important role, not only in the temperature distribution, but also in arc disruptions. Moreover, disruptions are found to be influenced by convection only for a limited frequency range. The results of the present studies are applicable to the frequency range of 10-10(2) Hz which includes most industry ac arc frequencies. (C) 1994 Academic Press, Inc.

A new model for the floating potential of fine particles in plasma

In the plasma processing of ultrafine particles of material, the heat transfer and force are considerably affected by particle charging. In this communication a new model, including thermal electron emission and incorporating the effect of electric field near the particle surface, is developed for metallic spherical particles under the condition of a thin plasma sheath. Based on this model, the particle floating potential, and thus the heat transfer and force, can be detemined more accurately and more realistically than previously.

Fusion energy-production from a deuterium-tritium plasma in the jet tokamak

Describes a series of experiments in the Joint European Torus (JET), culminating in the first tokamak discharges in deuterium-tritium fuelled mixture. The experiments were undertaken within limits imposed by restrictions on vessel activation and tritium usage. The objectives were: (i) to produce more than one megawatt of fusion power in a controlled way; (ii) to validate transport codes and provide a basis for accurately predicting the performance of deuterium-tritium plasmas from measurements made in deuterium plasmas; (iii) to determine tritium retention in the torus systems and to establish the effectiveness of discharge cleaning techniques for tritium removal; (iv) to demonstrate the technology related to tritium usage; and (v) to establish safe procedures for handling tritium in compliance with the regulatory requirements. A single-null X-point magnetic configuration, diverted onto the upper carbon target, with reversed toroidal magnetic field was chosen. Deuterium plasmas were heated by high power, long duration deuterium neutral beams from fourteen sources and fuelled also by up to two neutral beam sources injecting tritium. The results from three of these high performance hot ion H-mode discharges are described: a high performance pure deuterium discharge; a deuterium-tritium discharge with a 1% mixture of tritium fed to one neutral beam source; and a deuterium-tritium discharge with 100% tritium fed to two neutral beam sources. The TRANSP code was used to check the internal consistency of the measured data and to determine the origin of the measured neutron fluxes. In the best deuterium-tritium discharge, the tritium concentration was about 11% at the time of peak performance, when the total neutron emission rate was 6.0 × 10¹⁷ neutrons/s. The integrated total neutron yield over the high power phase, which lasted about 2 s, was 7.2 × 10¹⁷ neutrons, with an accuracy of ±7%. The actual fusion amplification factor, Q_DT was about 0.15

Magnetohydrodynamics equilibrium of a self-confined elliptic plasma ball

A variational principle is applied to the problem of magnetohydrodynamics (MHD) equilibrium of a self-contained elliptical plasma ball, such as elliptical ball lightning. The principle is appropriate for an approximate solution of partial differential equations with arbitrary boundary shape. The method reduces the partial differential equation to a series of ordinary differential equations and is especially valuable for treating boundaries with nonlinear deformations. The calculations conclude that the pressure distribution and the poloidal current are more uniform in an oblate self-confined plasma ball than that of an elongated plasma ball. The ellipticity of the plasma ball is obviously restricted by its internal pressure, magnetic field, and ambient pressure. Qualitative evidence is presented for the absence of sighting of elongated ball lightning.

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.

The induction plasma chemical reactor: Part II. Kinetic model

A kinetic model has been developed for the prediction of the concentration gelds in an rf plasma reactor. A sample calculation for a SiCl4/H2 system is then performed. The model considers the mixing processes along with the kinetics of seven reactions involving the decomposition of these reactants. The results obtained are compared to those assuming chemical equilibrium. The predictions indicate that an equilibrium assumption will result in lower predicted temperature fields in the reactor. Furthermore, for the chemical system considered here, while differences exist between the concentration fields obtained by the two models, the differences are not substantial.

Magnetohydrodynamic equilibrium of plasma ball lightning

Thin foil observations using transmission electron microscopy reveal that the density of dislocations within the band is extremely high and the tangled arrangement of dislocations tends to align along the length of the shear band. The grains in the band were also elongated along the shear band and clearly exhibited a crystallographic nature.

Adiabatic compression of plasma with purely toroidal rotation

A variational method is developed to find approximate solutions to the generalized Grad-Shafranov equations for an adiabatic compression of the plasma with toroidal rotation, via the expansion in Fourier series in poloidal angle of the flux surface coordinates. The numerical results, which are carried out by the present method and by the usual two-dimensional method for a static equilibrium state, agree well.

Variational method for adiabatic compression of plasma with poloidal and toroidal rotation

A variational method is developed for adiabatic compression of plasma with both poloidal and toroidal rotation.