Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “PLASMA CONFINEMENT”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 235 records · Page 13

A note on the application of the Prigogine theorem to rotation of tokamak-plasmas in absence of external torques

Rotation of tokamak-plasmas, not at the mechanical equilibrium, is investigated using the Prigogine thermodynamic theorem. This theorem establishes that, for systems confined in rectangular boxes, the global motion of the system with barycentric velocity does not contribute to dissipation. This result, suitably applied to toroidally confined plasmas, suggests that the global barycentric rotations of the plasma, in the toroidal and poloidal directions, are pure reversible processes. In case of negligible viscosity and by supposing the validity of the balance equation for the internal forces, we show that the plasma, even not in the mechanical equilibrium, may freely rotate in the toroidal direction with an angular frequency, which may be higher than the neoclassical estimation. In addition, its toroidal rotation may cause the plasma to rotate globally in the poloidal direction at a speed faster than the expression found by the neoclassical theory. The eventual configuration is attained when the toroidal and poloidal angular frequencies reaches the values that minimize dissipation. The physical interpretation able to explain the reason why some layers of plasma may freely rotate in one direction while, at the same time, others may freely rotate in the opposite direction, is also provided. Invariance properties, herein studied, suggest that the dynamic phase equation might be of the second order in time. We then conclude that a deep and exhaustive study of the invariance properties of the dynamical and thermodynamic equations is the most correct and appropriate way for understanding the triggering mechanism leading to intrinsic plasma-rotation in toroidal magnetic configurations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Development of a radio frequency ion source with multi-helicon plasma injectors for neutral beam injection system of Versatile Experiment Spherical Torus

Despite of high plasma density, helicon plasma has not yet been applied to a large area ion source such as a driver for neutral beam injection (NBI) system due to intrinsically poor plasma uniformity in the discharge region. In this study, a radio-frequency (RF) ion source with multi-helicon plasma injectors for high plasma density with good uniformity has been designed and constructed for the NBI system of Versatile Experiment Spherical Torus at Seoul National University. The ion source consists of a rectangular plasma expansion chamber (120 × 120 × 120 mm{sup 3}), four helicon plasma injectors with annular permanent magnets and RF power system. Main feature of the source is downstream plasma confinement in the cusp magnetic field configuration which is generated by arranging polarities of permanent magnets in the helicon plasma injectors. In this paper, detailed design of the multi-helicon plasma injector and plasma characteristics of the ion source are presented.

43 PARTICLE ACCELERATORS↗

Confinement improvement by fluctuating input power

An interesting effect of the fluctuating input power on plasma confinement is reported. Specifically, the addition of the sinusoidal perturbation, Asin(wt), to constant power q{sub 0} is shown to promote the confinement, leading to the L-H transition at a lower value of q{sub 0}, as compared to the case of constant q{sub 0} without the sinusoidal perturbation. In general, higher amplitude (A) and lower frequency (w) are found to be more favorable for the L-H transition while an interesting linear relation between A and w leading to the L-H transition is established for different values of q{sub 0}.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A flux-balanced fluid model for collisional plasma edge turbulence: Model derivation and basic physical features

We propose a new reduced fluid model for the study of the drift wave–zonal flow dynamics in magnetically confined plasmas. Our model can be viewed as an extension of the classic Hasegawa-Wakatani (HW) model and is based on an improved treatment of the electron dynamics parallel to the field lines, to guarantee a balanced electron flux on the magnetic surfaces. Our flux-balanced HW (bHW) model contains the same drift-wave instability as previous HW models, but unlike these models, it converges exactly to the modified Hasegawa-Mima model in the collisionless limit. We rely on direct numerical simulations to illustrate some of the key features of the bHW model, such as the enhanced variability in the turbulent fluctuations and the existence of stronger and more turbulent zonal jets than the jets observed in other HW models, especially for high plasma resistivity. Our simulations also highlight the crucial role of the feedback of the third-order statistical moments in achieving a statistical equilibrium with strong zonal structures. Finally, we investigate the changes in the observed dynamics when more general dissipation effects are included and, in particular, when we include the reduced model for ion Landau damping originally proposed by Wakatani and Hasegawa.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Nonlinear control of the minimum safety factor in tokamaks by optimal allocation of spatially moving electron cyclotron current drive

The minimum value of the safety factor profile is related to the magnetohydrodynamic (MHD) stability of the plasma confined in a tokamak. Therefore, active control of the minimum safety factor may mitigate MHD instabilities that can degrade or even terminate plasma confinement. Typically, in most tokamak scenarios, the minimum safety factor evolves spatially with time, i.e., the location at which the safety factor achieves the minimum value changes with time. In addition to the inherent nonlinearities in the minimum safety factor evolution, its spatial variation makes the control design challenging. In particular, complexity in control design may arise from the need for time-dependent nonlinear models that account for spatial variation of the minimum safety factor. Furthermore, the minimum safety factor may drift to locations where the actuator authority is low. The problem of minimum safety factor control with target location tracking and moving electron cyclotron current drive (ECCD) is addressed in this work. A nonlinear time-dependent model that incorporates the spatial variation of the minimum safety factor is presented. A nonlinear controller based on optimal feedback linearization is developed to track a target minimum safety factor. The proposed controller treats the ECCD position as a controllable variable. In other words, the controller prescribes the ECCD position (in addition to the non-inductive powers) in real time based on an optimal criterion that is defined a priori. This work also presents the steps necessary to integrate the minimum safety factor controller with a total energy controller to achieve multiple control objectives simultaneously. In conclusion, the proposed integrated control algorithm is tested using nonlinear simulations in the Control Oriented Transport SIMulator (COTSIM) for a DIII-D tokamak scenario.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

MPEX AI Digital Twins

All magnetically confined plasma fusion power plant concepts (Tokamak, Spherical Tokamak, Stellarator, Mirror, ...) must exhaust the heat and plasma from the core confinement region to the material walls. The primary channel for this exhaust is through a plasma divertor which directs plasma along open magnetic field lines to a material target. The Material Plasma Exposure eXperiment (MPEX) illustrated in Figure 1, is a high-power, steady-state linear plasma device designed to produce the plasma material interaction (PMI) conditions of the divertor of future magnetic confinement fusion power plants: energy flux 20MW/m 2 , ion fluence 1031/m 2 , pulse duration 106 sec. These goals of plasma exposure in MPEX are well beyond those achieved in magnetic fusion experimental devices. Successfully achieving these high power steady state conditions for long pulses requires operational control of the heating and particle sources and the plasma flux to the walls and target. The MPEX AI Hot Spot Controller, proposed in this project, will help achieve the operational milestones of MPEX. The MPEX device will begin commissioning at the end of FY26. A smaller proto-MPEX was operated for 14,666 plasma discharges and will resume operation in September of 2025 as proto-MPEX-lite, with reduced capability, to test a new window for the Helicon plasma source. The proto-MPEX data has undergone surrogate modeling with machine learning methods (R. Archibald, 2022 IEEE International Conference on Big Data). This proto-MPEX data will be used to begin development of the AI digital twins described in this white paper. The scientific mission of MPEX is to qualify materials of different composition for use in the high energy and plasma flux conditions of a fusion power plant. The materials exposed in MPEX will in some cases be exposed to high neutron fluxes at other ORNL facilities to measure the changes to their PMI properties. The targets exposed in MPEX will be transported under vacuum to a Surface Analysis Station (SAS). The SAS will be equipped with the following diagnostics: Focused Ion Beam (FIB) for trench milling, 100-400 angstrom resolution scanning electron microscope (SEM), surface mapping x-ray spectrometer, high resolution camera, and a future upgrade to a laser induced breakdown spectroscopy quadruple mass spectrometer (LIBS-QMS). The MPEX experiments will generate diverse pre- and post-exposure measurement data of detailed material properties down to the crystal grain level in 3D for post-exposure assessment of PMI damage (e.g. cracking, melting, erosion and redeposition of the material). Physics models for the PMI, and how the material composition and manufacturing impact its performance under high energy plasma exposure, need to be validated with MPEX data to guide the selection of new candidate materials. Our vision for the MPEX AI Digital Twins project is to supply experimental and physics model simulation data to train Artificial Intelligence (AI) models for data processing, analysis, operational control, PMI and materials simulation to maximize the scientific output of the MPEX device. Ultimately, an AI digital twin of MPEX material assessment metrics for tested and synthetic material types with simulated PMI will be trained by the AI Modeling Teams on the experimental and physics simulation data submitted to the American Science Cloud by this project. A purely empirical search for the best material is inefficient given the finite number of samples that can be tested on MPEX. In order to expand the material properties database for training the MPEX Material Assessment AI Digital Twin, and to gain physics understanding of the PMI processes, physics models of the material properties and PMI processes are required. The physics simulations provide detailed simulation data, like impact angles for plasma ions, sputtering yields, transport of the ionized sputtered target material in the plasma, and redeposition locations. This simulation data expands the measurement data for deeper physics understanding. The experimental data is essential to validate the PMI and material structure simulation models. The validated models can then be used to generate new simulation data of MPEX material assessments for synthetic material compositions that have not been exposed in MPEX. These predictive simulations, plus the whole experimental dataset, will be used to train the MPEX Material Assessment AI Digital Twin allowing a rapid generative AI search for new materials with reduced PMI damage by interpolating the domain of the training set. These new optimum materials can be simulated with the physics codes and/or tested in MPEX. The ability of AI neural networks to interpolate multi-dimensional parameter spaces and generate virtual data is exploited for a more efficient search for optimum materials. The advent of the Transformational AI Models Consortium (TAIMC) is an opportunity to engage with state of the art private and public AI developers to achieve the goals of the AI digital twins and AI accelerated physics models proposed in this project. Our partners at ORNL from the Advance Scientific Computing Research (ASCR) organization will collaborate in accelerating the integrated plasma material interaction simulation framework. This simulation framework will provide a platform for generating simulation data across a range of physical fidelities, including hybrid methods that produce multi-fidelity results. This data will be leveraged for AI model development, both for generation of surrogates and the automation of simulation campaigns. A part of the research below will include collaborative efforts with the TAIMC to (i) adapt data storage approaches to ensure AI-readiness, (ii) provide a protypical exemplar to inform and exercise constructed workflows, and (iii) generate and share data, using the TAIMC unified AI data standard, for foundational models that will be trained from multiple sources across the DOE complex. We will also collaborate with the TAIMC, as well as the planned AI modeling teams, to develop approaches for reducing the cost of data generation. These include tailored multi-fidelity approaches as well as fine-tuning strategies to augment general, large-scale foundational models.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Analysis of the linear and nonlinear stability of Alfven eigenmodes and fish-bones in JET DT discharges: mode identification and shear flows generation

The plasma in future nuclear fusion reactors will be heated by neutral beam injectors (NBIs) and high frequency electromagnetic waves as well as fusion born alpha particles. Energetic particles (EPs), with energies up to two orders of magnitude larger than the thermal plasma, can trigger EP driven modes and induce harmful EP losses, reducing the plasma heating efficiency and the economical viability of the reactor. The present study is dedicated to analyze the Alfven Eigenmode (AE) activity in JET D–T discharges, the closest experiment to reactor-like operation performed until now. There, EP driven modes are induced by the combined effect of tangential NBIs and ion cyclotron resonance heating (ICRH) driven EP. Linear and nonlinear simulations are performed with the gyro-fluid FAR3d code to analyze the AE activity observed in the discharge 99896. The linear simulations reproduce the unstable n = 3 to 5 toroidal AEs (TAE) at the inner plasma region observed in the experiment, triggered by highly energetic passing deuterium populations injected by the tangential NBIs, further accelerated by the effect of the ICRH up to 1 MeV. In addition, fish-bones triggered by energetic trapped hydrogen induced by the ICRH are also reproduced. On the other hand, the alpha particles density is too small to destabilize AEs in the experiment. Nonetheless, increasing artificially the alpha density by one order of magnitude, an n = 1 beta induced AE can be destabilized in the inner plasma region. Nonlinear simulations indicate the generation of zonal structures during the AE/fish-bone saturation phase. TAE and fish-bones causes a rather weak increase of the passing D and trapped H EP (around 2%), respectively. Shear flows and zonal currents are generated during the saturation of TAE and fish-bones. Nonlinear simulations performed for D–T and pure deuterium thermal plasma indicate AE/fish-bone activity is weaker and shear flows are less intense in the pure deuterium case, trends consistent with the experimental observations that also indicates a deterioration of the thermal plasma confinement. Therefore, both numerical studies and experimental evidence indicate the generation of shear flows by AE/fish-bones could be connected with an improvement of the thermal plasma confinement.

AE↗

Fast camera studies at an electron cyclotron resonance table plasma generator

A simple table-size ECR plasma generator operates in the ATOMKI without axial magnetic trap and without any particle extraction tool. Radial plasma confinement is ensured by a NdFeB hexapole. The table-top ECR is a simplified version of the 14 GHz ATOMKI-ECRIS. Plasma diagnostics experiments are planned to be performed at this device before installing the measurement setting at the “big” ECRIS. Recently, the plasma generator has been operated in pulsed RF mode in order to investigate the time evolution of the ECR plasma in two different ways. (1) The visible light radiation emitted by the plasma was investigated by the frames of a fast camera images with 1 ms temporal resolution. Since the visible light photographs are in strong correlation with the two-dimensional spatial distribution of the cold electron components of the plasma it can be important to understand better the transient processes just after the breakdown and just after the glow. (2) The time-resolved ion current on a specially shaped electrode was measured simultaneously in order to compare it with the visible light photographs. The response of the plasma was detected by changing some external setting parameters (gas pressure and microwave power) and was described in this paper.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Transport studies in high-performance field reversed configuration plasmas

A significant improvement of field reversed configuration (FRC) lifetime and plasma confinement times in the C-2 plasma, called High Performance FRC regime, has been observed with neutral beam injection (NBI), improved edge stability, and better wall conditioning [Binderbauer et al., Phys. Plasmas 22, 056110 (2015)]. A Quasi-1D (Q1D) fluid transport code has been developed and employed to carry out transport analysis of such C-2 plasma conditions. The Q1D code is coupled to a Monte-Carlo code to incorporate the effect of fast ions, due to NBI, on the background FRC plasma. Numerically, the Q1D transport behavior with enhanced transport coefficients (but with otherwise classical parametric dependencies) such as 5 times classical resistive diffusion, classical thermal ion conductivity, 20 times classical electron thermal conductivity, and classical fast ion behavior fit with the experimentally measured time evolution of the excluded flux radius, line-integrated density, and electron/ion temperature. The numerical study shows near sustainment of poloidal flux for nearly 1 ms in the presence of NBI.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Verification of electromagnetic simulation capabilities in global gyrokinetic particle-in-cell code GTS

Recently, the numerical scheme presented by Mishchenko et al. enabled explicit gyrokinetic simulations of low-frequency electromagnetic instabilities in tokamaks at experimentally relevant values of plasma β⁠. This scheme resolved the long-standing cancellation problem that previously hindered gyrokinetic particle-in-cell code simulations of magnetohydrodynamic phenomena with inherently small parallel electric fields. Moreover, the scheme did not employ approximations that eliminate critical tearing-type instabilities. Here, we report on the implementation of this numerical scheme in the global gyrokinetic particle-in-cell code GTS. This implementation allows for a more complete and accurate picture of interaction between small scale turbulence and MHD modes in tokamaks. Additionally, we present a comprehensive set of verification simulations of numerous electromagnetic instabilities relevant to present-day tokamaks. These simulations encompass the kinetic ballooning mode, the internal kink mode, the tearing mode, the micro-tearing mode, and the toroidal Alfven eigenmode destabilized by energetic ions, which are all instrumental in understanding tokamak physics. We will also showcase the preliminary nonlinear simulations of kinetic ballooning instabilities and (2,1) island formation due to tearing mode instability. These simulations validate the accuracy of the scheme implementation and pave the way for studying how these instabilities affect plasma confinement and performance.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Spontaneous focusing of plasma flow in a weak perpendicular magnetic field

Structure formation of high-beta plasma flow in a perpendicular magnetic field is investigated in the ion kinetic regime by a fully kinetic particle-in-cell simulation. We demonstrate that directional plasma flow is spontaneously focused to form a sharp density structure. The primary focusing process comes from field-aligned electron inflow associated with the whistler mode and plasma confinement due to a self-generated magnetic field. The resulting concave magnetic field lines modulate ion gyration to cause a secondary focusing process with significant plasma concentration. Required conditions for these processes are determined by a dimensionless parameter α ≡ β{sub i0}(Δ{sub W}/ρ{sub i0}), where β{sub i0}, Δ{sub W}, and ρ{sub i0} denote the plasma kinetic beta, window size, and ion gyration radius, respectively. The focusing process is apparent for small α, whereas diamagnetic expansion is dominant for large α. This condition describes a transition between diamagnetic cavity formation and the focusing process.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Multimachine analysis of intermittent fluctuations in the scrape-off layer of magnetically confined fusion plasmas

Far scrape-off layer fluctuations, also known as blobs/filaments, are compared for the first time across six tokamaks: Alcator C-Mod, DIII-D, TCV, KSTAR, MAST, and MAST-U, under the context of stochastic modeling in Ohmic and low confinement mode scenarios. In this model, which well describes these fluctuations, measurements from Langmuir probes are modeled as a superposition of uncorrelated, two-sided exponential temporal pulses. The key parameters of this model are the mean waiting times between consecutive pulses, the pulse duration times and the mean of the pulse amplitudes. The study examines the influence of line-averaged density, plasma current, magnetic field and aspect ratio. The mean amplitudes and pulse waiting times increase with the line-averaged density across devices, while the pulse duration times remain constant. Whereas the variation in waiting times is gradual across the density scan, the mean amplitudes increase abruptly as the density profile broadens and flattens for Alcator C-Mod. The blob waiting times decrease with increasing plasma current and the toroidal magnetic field, while they increase with aspect ratio for conventional tokamaks. The duration times follow the same trends as the waiting times in these parameters but vary more slowly, resulting in the trends in the waiting times dominating the intermittency level of the fluctuations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Mode transitions and spoke structures in E × B Penning discharge

Two-dimensional particle-in-cell simulations in the (radial-azimuthal) plane perpendicular to the axial direction of a cylindrical $E$ x $B$ Penning discharge are presented. The low-pressure discharge is self-consistently supported by plasma ionization from the electron beam injected axially, along the direction of the external magnetic field. It is shown that with the increasing strength of the external magnetic field, the discharge undergoes a sequence of transitions between several azimuthal modes. Azimuthal m > 1 spiral arm structures are excited at low magnetic field values as plasma confinement improves and the radial density profile becomes peaked. With a larger field, spiral arms with m > 1 are replaced by the m = 1 spoke mode, most clearly seen in plasma density. A transition from spiral arms to the spoke regime occurs when the plasma potential in the center changes from weakly positive (or zero) to negative. Further increase in the magnetic field results in a well-developed m = 1 spoke mode with additional small-scale higher-frequency m > 1 structures inside and around the spoke. It is shown that while ionization and collisions affect some characteristics of the observed fluctuations, the basic features of the spoke and m > 1 spiral structure remained similar without ionization. The role of energy conservation in small-scale high-frequency modes and spoke dynamics is discussed. It is demonstrated that in regimes with the m = 1 spoke mode, additional m = 4 harmonics of the ion and electron fluxes to the wall appear due to the square boundary. The frequency of the m = 1 mode is weakly affected by the geometry of the boundary.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Experiments on the transportation of a magnetized plasma stream in the GOL-3 facility

The program of the deep upgrade of the GOL-3 multiple-mirror trap is presented. The upgrade is aimed at creating a new GOL-NB open trap located at the GOL-3 site and intended to directly demonstrate the efficiency of using multiple-mirror magnetic cells to improve longitudinal plasma confinement in a gasdynamic open trap. The GOL-NB device will consist of a new central trap, adjoint cells with a multiple-mirror magnetic field, and end tanks (magnetic flux expanders). Plasma in the central trap will be heated by neutral beam injection with a power of up to 1.5 MW and duration of 1 ms. At present, physical experiments directed at developing plasma technologies that are novel for this facility are being carried out using the 6-m-long autonomous part of the GOL-3 solenoid. The aim of this work was to develop a method for filling the central trap with a low-temperature start plasma. Transportation of a plasma stream from an arc source over a distance of 3 m in a uniform magnetic field with an induction of 0.5–4.5 T is demonstrated. In these experiments, the axial plasma density was (1–4) × 10{sup 20} m{sup –3} and the mirror ratio varied from 5 to 60. In general, the experiments confirmed the correctness of the adopted decisions for the start plasma source of the GOL-NB device.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Design and implementation of a portable diagnostic system for Thomson scattering and optical emission spectroscopy measurements

A diagnostic system, which has a design goal of high-portability, has been designed at Oak Ridge National Laboratory (ORNL). This project aims at providing measurements of key plasma parameters (ne, Te, ni, Ti) for fusion-relevant devices, utilizing Thomson scattering (TS) and optical emission spectroscopy (OES). The innovative design employs mostly commercial off-the-shelf instrumentation and a traveling team of researchers to conduct measurements at various magnetic-confinement plasma devices. The TS diagnostic uses a Quantel Q-smart 1500 Nd:YAG laser with a 2ω harmonic generator to produce up to 850 mJ of 532 nm laser pulses at 10 Hz. Collection optics placed at the detection port consists of an 11 × 3 optical fiber bundle, where the TS diagnostic uses an 11 × 1 subset array of the fibers, the OES diagnostic uses another 11 fibers, and the remaining fibers are available to the host institution. The detection system is comprised of two separate IsoPlane-320 spectrometers with triple-grating turrets of various line spacing and two PI-MAX 4 intensified CCD detectors, used simultaneously to measure a broad range of ion, impurity, and electron parameters. The self-contained diagnostic package also includes a data processing and storage system. The design and initial implementation of the TS-OES diagnostic system are described. The experiments from the proof-of-principle operation of the portable package on a high density (∼2.5 × 1022 m−3) and low-temperature (∼5 eV) electrothermal arc source at ORNL are also discussed.

Kafle, N. (ORCID:000000028144787X)↗

DESIGN AND TESTING OF DIMES CARBON ABLATION RODS IN THE DIII-D TOKAMAK

We present the design of ATJ graphite rods developed for ablation experiments under high heat flux (up to 50 MW/m2) in the lower divertor of the DIII-D tokamak, a magnetic plasma confinement device. This work is motivated by the need to test ablation models relevant to carbon-based thermal shields used in high-speed spacecraft atmospheric entries, where the heat fluxes encountered can be comparable to those achieved in the DIII-D divertor plasma. Several different designs for the flowfacing side of the rod are analyzed, including "sharp nose," "blunt," and "concave". The last shape is studied for its potential to lower heat fluxes at the rod surface by increased radiation from trapped neutrals and reduced parallel plasma pressure. We also analyze the possibility of applying a thin (approximately 30 microns) layer of silicon carbide (SiC) to the exposed part of several carbon ablation rods to benchmark its erosion calculations and lifetime predictions. Such calculations are of interest as SiC represents a promising material for both thermal protection systems (TPS) and a fusion plasma-facing material (PFM). Preliminary results from the DIII-D rod ablation experiments are also discussed.

Orlov, Dmitri M.↗

Plasma instability in the afterglow of electron cyclotron resonance discharge sustained in a mirror trap

The work presented in this article is devoted to time-resolved diagnostics of non-linear effects observed during the afterglow plasma decay of a 14 GHz electron cyclotron resonance ion source operated in pulsed mode. Plasma instabilities that cause perturbations of the extracted ion current during the decay were observed and studied. It is shown that these perturbations are associated with precipitation of high energy electrons along the magnetic field lines and strong bursts of bremsstrahlung emission. The effect of ion source settings on the onset of the observed instabilities was investigated. Based on the experimental data and estimated plasma properties, it is assumed that the instabilities are of cyclotron type. The conclusion is supported by a comparison to other types of plasma devices which exhibit similar characteristics but which operate in a different plasma confinement regime.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Study of lower hybrid wave propagation in ionized gas by Hamiltonian theory

In order to find an approximate solution to the Vlasov-Maxwell equation system describing the lower hybrid wave propagation in magnetic confined plasmas, the use of the WKB method leads to the ray tracing equations. The Hamiltonian character of the ray tracing equations is investigated analytically and numerically in order to deduce the physical properties of the wave propagating without absorption in the confined plasma. The consequences of the Hamiltonian character of the equations on the travelling wave, in particular, on the evolution of the parallel wavenumber along the propagation path have been accounted and the chaotic diffusion of the timeaveraged parallel wave-number towards higher values has been evaluated. Numerical analysis by means of a Runge-Kutta based algorithm implemented in a ray tracing code supplies the analytical considerations. A numerical tool based on the symplectic integration of the ray trajectories has been developed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗