圆截面轴对称等离子体环平衡的分析解和极限比压〈β_r〉

本文根据Mercier的基本方程用反纵横比ε的暴级数展开方法求出了胖Tokamak的平衡分析解。为了考查ε的影响,本文计算到了ε~3的量级,并在没有负压的限制下求得了相应的极限比压〈β_r〉。求得的平衡解亦可作为不稳定分析之用。

Jet:progress in performance and understanding

In 1990 JET operated with a number of technical improvements which led to advances in performance and permitted the carrying out of experiments specifically aimed at improving physics understanding of selected topics relevant to the "NEXT STEP". The new facilities include beryllium antenna screens, a prototype lower hybrid current drive system, and modification of the NI system to enable the injection of He-3 and He-4. Continued investigation of the hot-ion H-mode produced a value of n(D)(0)tau-E(T)(i)(0) = 9 x 10(20)m-3s keV, which is near conditions required for Q(DT) = 1, while a new peaked density profile H-mode was developed with only slightly lower performance. Progress towards steady state operation has been made by achieving ELMy H-modes under certain operating conditions, while maintaining good tau-E values. Experimental simulation of He ash transport indicates effective removal of alpha-particles from the plasma core for both L and H mode plasmas. Detailed analyses of particle and energy transport have helped establish a firmer link between particle and energy transport, and have suggested a connection between reduced energy transport and reversed shear. Numerical and analytic studies of divertor physics carried out for the pumped divertor phase of JET have helped clarify the key parameters governing impurity retention, and an intensive model validation effort has begun. Experimental simulation of alpha-particle effects with beta-fast up to 8% have shown that the slowing down processes are classical, and have given no evidence of deleterious collective effects.

Robust Sliding Mode Control for Tokamaks

Nuclear fusion has arisen as an alternative energy to avoid carbon dioxide emissions, being the tokamak a promising nuclear fusion reactor that uses a magnetic field to confine plasma in the shape of a torus. However, different kinds of magnetohydrodynamic instabilities may affect tokamak plasma equilibrium, causing severe reduction of particle confinement and leading to plasma disruptions. In this sense, numerous efforts and resources have been devoted to seeking solutions for the different plasma control problems so as to avoid energy confinement time decrements in these devices. In particular, since the growth rate of the vertical instability increases with the internal inductance, lowering the internal inductance is a fundamental issue to address for the elongated plasmas employed within the advanced tokamaks currently under development. In this sense, this paper introduces a lumped parameter numerical model of the tokamak in order to design a novel robust sliding mode controller for the internal inductance using the transformer primary coil as actuator.

Quantifying chaos: practical estimation of the correlation dimension

Be it a physical object or a mathematical model, a nonlinear dynamical system can display complicated aperiodic behavior, or "chaos." In many cases, this chaos is associated with motion on a strange attractor in the system's phase space. And the dimension of the strange attractor indicates the effective number of degrees of freedom in the dynamical system.

In this thesis, we investigate numerical issues involved with estimating the dimension of a strange attractor from a finite time series of measurements on the dynamical system.

Of the various definitions of dimension, we argue that the correlation dimension is the most efficiently calculable and we remark further that it is the most commonly calculated. We are concerned with the practical problems that arise in attempting to compute the correlation dimension. We deal with geometrical effects (due to the inexact self-similarity of the attractor), dynamical effects (due to the nonindependence of points generated by the dynamical system that defines the attractor), and statistical effects (due to the finite number of points that sample the attractor). We propose a modification of the standard algorithm, which eliminates a specific effect due to autocorrelation, and a new implementation of the correlation algorithm, which is computationally efficient.

Finally, we apply the algorithm to chaotic data from the Caltech tokamak and the Texas tokamak (TEXT); we conclude that plasma turbulence is not a low- dimensional phenomenon.

Topology of toroidal helical fields in non-circular cross-sectional tokamaks

The ordinary differential magnetic field line equations are solved numerically; the tokamak magnetic structure is studied on Hefei Tokamak-7 Upgrade (HT-7U) when the equilibrium field with a monotonic q-profile is perturbed by a helical magnetic field. We find that a single mode (m, n) helical perturbation can cause the formation of islands on rational surfaces with q = m/n and q = (m +/- 1, +/- 2, +/- 3,...)/n due to the toroidicity and plasma shape (i.e. elongation and triangularity), while there are many undestroyed magnetic surfaces called Kolmogorov-Arnold-Moser (KAM) barriers on irrational surfaces. The islands on the same rational surface do not have the same size. When the ratio between the perturbing magnetic field B-r(r) and the toroidal magnetic field amplitude B(phi)0 is large enough, the magnetic island chains on different rational surfaces will overlap and chaotic orbits appear in the overlapping area, and the magnetic field becomes stochastic. It is remarkable that the stochastic layer appears first in the plasma edge region.

Low Effort L-i Nuclear Fusion Plasma Control Using Model Predictive Control Laws

One of the main problems of fusion energy is to achieve longer pulse duration by avoiding the premature reaction decay due to plasma instabilities. The control of the plasma inductance arises as an essential tool for the successful operation of tokamak fusion reactors in order to overcome stability issues as well as the new challenges specific to advanced scenarios operation. In this sense, given that advanced tokamaks will suffer from limited power available from noninductive current drive actuators, the transformer primary coil could assist in reducing the power requirements of the noninductive current drive sources needed for current profile control. Therefore, tokamak operation may benefit from advanced control laws beyond the traditionally used PID schemes by reducing instabilities while guaranteeing the tokamak integrity. In this paper, a novel model predictive control (MPC) scheme has been developed and successfully employed to optimize both current and internal inductance of the plasma, which influences the L-H transition timing, the density peaking, and pedestal pressure. Results show that the internal inductance and current profiles can be adequately controlled while maintaining the minimal control action required in tokamak operation.

Extreme ultraviolet emission lines of Ca XV in solar and laboratory spectra

New R-matrix calculations of electron impact excitation rates in Ca XV are used to derive theoretical electron density diagnostic emission line intensity ratios involving 2s(2)2p(2)- 2s2p(3) transitions, specifically R-1 = I(208.70 Angstrom)/I(200.98 Angstrom), R-2 = I(181.91 Angstrom)/I(200.98 Angstrom), and R-3 = I(215.38 Angstrom)/I(200.98 Angstrom), for a range of electron temperatures (T-e = 10(6.4)-10(6.8) K) and densities (Ne = 10(9)-10(13) cm(-3)) appropriate to solar coronal plasmas. Electron densities deduced from the observed values of R-1, R-2, and R-3 for several solar flares, measured from spectra obtained with the Naval Research Laboratory's S082A spectrograph on board Skylab, are found to be consistent. In addition, the derived electron densities are in excellent agreement with those determined from line ratios in Ca XVI, which is formed at a similar electron temperature to Ca XV. These results provide some experimental verification for the accuracy of the line ratio calculations, and hence the atomic data on which they are based. A set of eight theoretical Ca XV line ratios involving 2s(2)2p(2)-2s2p(3) transitions in the wavelength range similar to140-216 Angstrom are also found to be in good agreement with those measured from spectra of the TEXT tokamak plasma, for which the electron temperature and density have been independently determined. This provides additional support for the accuracy of the theoretical line ratios and atomic data.

Extreme ultraviolet transitions of Fe XXI in solar, stellar and laboratory spectra

Recent R-matrix calculations of electron impact excitation rates for transitions among the 2s(2)2p(2), 2s2p(3) and 2p(4) levels of Fe XXI are used to derive theoretical electron density (N-e) sensitive emission-line ratios involving 2S2(2)p(2)-2s2p(3) transitions in the similar to 98-146 Angstrom wavelength range. A comparison of these with observations from the PLT tokamak plasma, for which the electron density has been independently determined, reveals generally very good agreement between theory and experiment, and in some instances removes discrepancies found previously. The observed Fe XXI ratios for a solar flare, obtained with the OSO-5 satellite, imply electron densities which are consistent, with discrepancies that do not exceed 0.2 dex. In addition, the derived values of N-e are similar to those estimated for the high-temperature regions of other solar flares. The good agreement between theory and observation, in particular for the tokamak spectra, provides experimental support for the accuracy of the present line-ratio calculations, and hence for the atomic data on which they are based.

Improved data for solar flare X-ray spectral analysis

The inclusion of collisional rates for He-like Fe and Ca ions is discussed with reference to the analysis of solar flare Fe XXV and Ca XIX line emission, particularly from the Yohkoh Bragg Crystal Spectrometer (BCS). The new data are a slight improvement on calculations presently used in the BCS analysis software in that the discrepancy in the Fe XXV y and z line intensities (observed larger than predicted) is reduced. Values of electron temperature from satellite-to-resonance line ratios are slightly reduced (by up to 1 MK) for a given observed ratio. The new atomic data will be incorporated in the Yohkoh BCS databases. The data should also be of interest for the analysis of high-resolution, non-solar spectra expected from the Constellation-X and Astro-E space missions. A comparison is made of a tokamak S XV spectrum with a synthetic spectrum using atomic data in the existing software and the agreement is found to be good, so validating these data for particularly high-n satellite wavelengths close to the S XV resonance line. An error in a data file used for analyzing BCS Fe XXVI spectra is corrected, so permitting analysis of these spectra.

Spectroscopy of astrophysical plasmas in the laboratory

Most of the matter in the universe is in the few form of a plasma. Over the past years physicists have produced laboratory plasmas that can mimic those observed in astrophysics. The best known is probably the tokamak, which has similar physical conditions and plasma processes to those found in collisionally dominated solar and stellar transition regions and coronae. Spectroscopy of such laboratory plasmas, in, particular at, ultraviolet and X-ray wavelengths, has greatly aided our understanding of their astrophysical counterparts. More recently, experiments have been performed on the Z Machine at the Sandia National Laboratory in the USA with the aim of creating, for the first time, steady-state photoionization-dominated plasmas that recreate the conditions found in some accretion-powered X-ray sources, such as X-ray binaries. In the future, experiments are envisaged with laser-produced plasmas at AWE Aldermaston that may be able to mimic the steady-state conditions found in high-energy accretion-powered sources, including the central regions of active galaxies.

Progress at JET in integrating ITER-relevant core and edge plasmas within the constraints of an ITER-like wall

This paper reports the progress made at JET-ILW on integrating the requirements of the reference ITER baseline scenario with normalized confinement factor of 1, at a normalized pressure of 1.8 together with partially detached divertor whilst maintaining these conditions over many energy confinement times. The 2.5 MA high triangularity ELMy H-modes are studied with two different divertor configurations with D-gas injection and nitrogen seeding. The power load reduction with N seeding is reported. The relationship between an increase in energy confinement and pedestal pressure with triangularity is investigated. The operational space of both plasma configurations is studied together with the ELM energy losses and stability of the pedestal of unseeded and seeded plasmas. The achievement of stationary plasma conditions over many energy confinement times is also reported.

Electron-impact excitation and ionization of W 3+ for the determination of tungsten influx in a fusion plasma

Tungsten will be employed as a plasma facing material in the ITER fusion reactor under construction in Cadarache, France; therefore, there is a significant need for accurate electron-impact excitation and ionization data for the ions of tungsten. We report on the results of extensive calculations of ionization and excitation for W 3+ that are intended to provide the atomic data needed for the determination of impurity influx diagnostics of tungsten in several existing tokamak reactors. The electron-impact excitation rate coefficients for this study were determined using the relativistic R -matrix method. The contribution to direct electron-impact ionization was determined using the distorted-wave approximation, the accuracy of which was verified by an R -matrix with pseudo states calculation. Contributions to total ionization from excitation autoionization were also generated from the relativistic R -matrix method. These results were then employed to calculate values of ionization per emitted photon, or SXB ratios, for four carefully selected spectral lines; these data will allow the determination of impurity influx from tungsten facing surfaces. For the range of densities of importance in the edge region of a tokamak reactor, these SXB ratios are found to be nearly independent of electron density but vary significantly with electron temperature.

The effects of radiative cascades on the x-ray diagnostic lines of Fe16+

We present complete collisional-radiative modelling results for the soft x-ray emission lines of Fe16+ in the 15 Å–17 Å range. These lines have been the subject of much controversy in the astrophysical and laboratory plasma community. Radiative transition rates are generated from fully relativistic atomic structure calculations. Electron-impact excitation cross sections are determined using a fully relativistic R-matrix method employing 139 coupled atomic levels through n = 5. We find that, in all cases, using a simple ratio of the collisional rate coefficient times a radiative branching factor is not sufficient to model the widely used diagnostic line ratios. One has to include the effects of collisional-radiative cascades in a population model to achieve accurate line ratios. Our line ratio results agree well with several previous calculations and reasonably well with tokamak experimental measurements, assuming a Maxwellian electron-energy distribution. Our modelling results for four EBIT line ratios, assuming a narrow Gaussian electron-energy distribution, are in generally poor agreement with all four NIST measurements but are in better agreement with the two LLNL measurements. These results suggest the need for an investigation of the theoretical polarization calculations that are required to interpret the EBIT line ratio measurements.

Application of Imaging to Study the Evolution and Dynamics of Laser Produced Plasma from Solid and Thin Film Li Targets

The present thesis report the results obtained from the studies carried out on the laser blow off plasma (LBO) from LiF-C (Lithium Fluoride with Carbon) thin film target, which is of particular importance in Tokamak plasma diagnostics. Keeping in view of its significance, plasma generated by the irradiation of thin film target by nanosecond laser pulses from an Nd:YAG laser over the thin film target has been characterized by fast photography using intensified CCD. In comparison to other diagnostic techniques, imaging studies provide better understanding of plasma geometry (size, shape, divergence etc) and structural formations inside the plume during different stages of expansion.

Estudi i aplicació de llibreries de física per a motors de videojocs

En el món dels videojocs el realisme és un punt molt important a tenir en compte ja que dónamés sensació a l’usuari d’estar immers en el videojoc. Això passa en part per aconseguir realisme en la dinàmica dels objectes i fer que aquests segueixin les lleis de la física de Newton. Per això s’han desenvolupat diverses llibreries que s’anomenen “motors de física” (physics engines), que empren variables com la massa, la velocitat, la fricció i la resistència del vent. Els objectius d’aquest projecte seran l’estudi de diferents llibreries físiques existents, la seva comparació i com s’integren en els motors de jocs. A més a més , la generació de contingut amb comportament que respongui a les funcions definides a aquestes llibreries no és trivial i per aquest motiu també es desenvoluparà una aplicació per generar murs de forma semiautomàtica que respongui a impactes. Per assolir aquests objectius caldrà: d’ una banda, comparar els cossos rígids, unions i funcionament en general de diferents llibreries físiques: Newton Game Dynamics, NVIDIA PhysX Technology, Open Dynamics Engine, Bullet Physics Library, Tokamak Physics Engine i Havok i d’ altra banda, implementar una aplicació que donant-li una imatge en planta d’una paret o conjunt de parets en format vectorial i les mides d’un maó, generi murs que puguin reaccionar de forma adequada quan rebin l’impacte d’una massa determinada. L’aplicació s’implementarà en C++ i amb l’entorn de desenvolupament Microsoft Visual Studio 2005. La visualització serà amb OpenGL