等离子体断路开关的一维磁流体力学数值模拟

 采用磁流体力学的方法对等离子体端路开关进行了一维数值模拟，通过对给定的进入等离子体的电流波形模型的计算，给出了轴向上等离子体的运动及其不同时刻的密度分布，分析了磁场在等离子体内部的异常渗透过程，指出考虑Hall效应是产生这一异常现象的最为主要的因素。     

Measurements of Flow Velocity in the Plasma Generator PSI‐2

The plasma flow velocity in the Plasma Generator PSI‐2 has been investigated by using of Mach probe. PSI‐2 is a stationary high‐current arc discharge in which the quasi‐neutral plasma expands along the magnetic field lines. The low‐temperature (T_e < 20 eV), medium density (n_e ∼ 10¹⁸— 10¹⁹ m^—3 ) plasma in the discharge is similar to the plasma in the divertor region of tokamaks. From the ratio of ion saturation currents collected from opposite sides of the probe the flow velocities (Mach numbers) in argon and hydrogen discharges are obtained.     

Anomalous Diffusion in a Linear Plasma Generator

Plasma production and particle transport of the magnetized plasma in the linear device PSI‐2 are investigated by analyzing the radial density and electron temperature pro.les obtained from Langmuir probe measurements. Additional information on the atomic (H) and proton temperatures is obtained from high resolution Doppler spectroscopy. The density profiles are found to be hollow, but do extend radially far beyond the visible rim with an exponential decay length of 3 cm. They can be explained by parallel and perpendicular diffusion in combination with a cylindrical source extending along the magnetic field all over the device. A large perpendicular diffusion coefficient in the range of 5 m²/s is inferred from the experimental results. (© 2008 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Role of poloidal E × B drift in divertor heat transport in DIII-D

Simulations for DIII-D high confinement mode plasmas with the multifluid code UEDGE show a strong role of poloidal E × B drifts on divertor heat transport, challenging the paradigm of conduction-limited scrape-off layer (SOL) transport. While simulations with reduced drift magnitude are well aligned with the assumption that electron heat conduction dominates the SOL heat transport, simulations with drifts predict that the poloidal convective E × B heat transport dominates over electron heat conduction in both attached and detached conditions. As poloidal E × B flow propagates across magnetic field lines, poloidal transport with shallow magnetic pitch angles can reach values that are of the same order as would be provided by sonic flows parallel to the field lines. These flows can lead to strong convection-dominated divertor heat transport, increasing the poloidal volume of radiative power front, consistent with previous measurements at DIII-D. Due to these convective flows, the Lengyel integral approach, assuming zero convective fraction, is expected to provide a pessimistic estimate for the radiative capability of impurities in the divertor. For the DIII-D simulations shown here, the Lengyel integral approach underestimates the radiated power by a factor of 6, indicating that, for reliable DIII-D divertor power exhaust predictions, full two-dimensional (2D) calculations, including drifts, would be necessary.     

The Analysis of Probe I‐V Characteristics in a Magnetized Low‐Temperature Plasma

The validity of probe theories based on the classical electron transport to probes in low‐temperature magnetized plasma of the toroidal device “Blaamann” has been demonstrated. The analysis was carried out for the conditions when global transport of charged particles in the device is anomalous, namely for magnetic field up to 0.3 T and pressure range 0.1—1 Pa. It was shown also that for the magnetic field larger than 0.1 T probes longer than 15 mm provide electron saturation current practically independent of probe potential, hence more accurate measurements of the plasma parameters. The experiments have revealed that application of long probes oriented parallel to the magnetic field may cause an anomaly of the I‐V characteristics in the sense that a local increase of the electron current appears near the plasma potential.     

Investigation of Hydrodynamic Properties of Hot Dense Plasma

Hydrodynamic properties of hot dense plasma bunch have been studied in terms of their application to the problem of inertial confinement fusion. The two-temperature hydrodynamic equations taking into account the electron heat conduction and kinetics of the charge composition are presented. The following model problems are considered: expansion of a plasma bunch in a vacuum, compression of the plasma by shock waves, interaction of the plasma bunch with a laser radiation and a beam of heavy ions. The some problems arising in the problem of inertial confinement fusion are analyzed on the basis of the developed numerical methods.     

Rotating relativistic electron beam-plasma interaction and formation of a field-reversed configuration

Experimental results on interaction of a rotating relativistic electron beam with plasma and neutral gas are presented. The rotating relativistic electron beam has been propagated up to a distance of 150 cm in a plasma. The response of the plasma to the rotating electron beam is found to be of magnetic diffusion type over a plasma density range 10¹¹–10¹³ cm⁻³. Excitation of the axial and azimuthal return currents by the rotating beam and subsequent trapping of the azimuthal return current layer by the magnetic mirror field are observed. A field-reversed configuration has been formed by the rotating relativistic electron beam when injected into neutral hydrogen gas. We have observed field reversal up to three times the initial field in an axial length of 100 cm.     

Thermal dynamics in the flux-coordinate independent turbulence code GRILLIX

GRILLIX employs the flux-coordinate independent approach (FCI), which allows us to study boundary plasma turbulence in realistic diverted configurations. Recently, the physical model in GRILLIX has been extended to a global drift-reduced Braginskii model, without any separation between background and fluctuations. It includes electromagnetic and thermal dynamics with hot ions, relaxation of the Boussinesq approximation and non-linear parametric dependencies. This contribution presents solutions to associated issues, that is, the ion diamagnetic polarization and the stiff parallel heat conduction. Simulations based on parameters characteristic for the Alcator C-Mod tokamak were carried out. In circular geometry, the self-consistent electric field contains zonal flows and geodesic acoustic modes in the confined region. In the scrape-off layer, the electron parallel heat conduction and its boundary condition determine the temperature and electric field, leading to sheared flows at the last closed flux surface.     

托卡马克等离子体不同平衡位形下的电子回旋波电流驱动

在给定的等离子体总电流和中心电流密度条件下,数值求解平衡方程,求出不同拉长比和三角形变因子的托卡马克等离子体温度、密度、磁场分布,然后通过求解波迹方程和Fokker-Planck方程,分别计算这些位形下的电子回旋波波迹和电流驱动.结果表明：电子回旋波X模从顶部发射时,随着拉长比的增大,波迹会向弱场侧偏移.电子回旋波X模从弱场侧发射时,电子回旋波在等离子体中传播沉积的功率份额随着拉长比的增大而增加,驱动电流位置随着三角形变因子的增大向等离子体中心移动.驱动电流位置随环向和极向发射角的减小向中心移动,对应的电流密度峰值也变大.     

A Self-Similar Model for Conduction in the Plasma Erosion Opening Switch

The conduction phase of the plasma erosion opening switch (PEOS) is characterized by combining a 1-D fluid model for plasma hydrodynamics, Maxwell's equations, and a 2-D electron-orbit analysis. A self-similar approximation for the plasma and field variables permits analytic expressions for their space and time variations to be derived. It is shown that a combination of axial MHD compression and magnetic insulation of high-energy electrons emitted from the switch cathode can control the character of switch conduction. The analysis highlights the need to include additional phenomena for accurate fluid modeling of PEOS conduction.     

Hot Ion Plasma Heating Experiments in SUMMA

Initial empirical results are presented for the hot-ion plasma heating experiments conducted in the new SUMMA (Superconducting Magnetic Mirror Apparatus) at NASA Lewis Research Center. A discharge was formed by applying a radially inward DC electric field near the mirror throats. Data were obtained at midplane magnetic flux densities from 1.0 to 3.5 tesla. Charge-exchange neutral particle energy analyzer data were reduced to ion temperatures using a plasma model that included a Maxwellian energy distribution super-imposed on an azimuthal drift, finite ion orbits, and radial variations in density and electric field. Using this plasma model, the highest ion temperatures computed were 5 keV, 1.2 keV, and 1 keV for He+, H2+, and H+, respectively. These were obtained at a mid-plane magnetic flux density of 1.6 T. Ion temperature was found to scale roughly as (P/B)n, where P/B is the ratio of power input to magnetic flux density and n is about 1 for hydrogen and 2 for helium. Optical spectroscopy line-broadening measurements yielded ion temperatures about 15 percent higher than the charge-exchange neutral particle analyzer results for hydrogen and about 50 percent higher for helium. Spectroscopically obtained electron temperatures ranged from 3 to 30 eV.     

Velocity distribution of plasma electrons in the negative H2- and He-glow with superimposed longitudinal magnetic field

The negative glow plasma has been found nearly field free in axial direction. Therefore plasma electrons in the stationary glow can thermalize down to the temperature of the neutral gas, whenever their diffusion—and recombination—lifetime is high enough. Applying Boltzmann's equation to this problem, the conditions of thermalization of plasma electrons are derived as a function of the outer parameters of the plasma: vessel diameter 2R, neutral gas pressurep and longitudinal magnetic fieldB. — If plasma electrons have a too short diffusion—and recombination—lifetime to be in thermal equilibrium with the neutral gas, the electron energy increases. For this case the distribution function of plasma electrons is derived using Boltzmann's equation. Approximating the calculated energy distribution by a Maxwellian distribution function, the electron temperature in the glow is obtained as a function of the parameters:R, p, B. OurT _e -measurements carried out in the H₂- and He-glows of different tube diameters, neutral gas pressures and magnetic fields agree closely with the theoretical results. TheT _e -measurements have been performed with Langmuir probes and by the method of reversal of the radial ambipolar electric field in a longitudinal magnetic field.     

Experimental observation and numerical analysis of arc plasmas diffused by magnetism

At atmospheric pressure, an intensified charge-coupled device (ICCD) camera with a narrow-band filter is used to capture the unsaturated images of a magnetically rotating arc. Comparison of its configurations with different arc current and external axial magnetic field (AMF) strength shows that the strong electromagnetic force may impel the arc to diffuse. Under the fully diffuse mode, the plasma is distributed throughout the electrode gap and no anode attachment can be seen in the cross-section of the torch. The fully diffuse plasma runs more steadily and its intensity distribution is more uniform, while its voltage fluctuation is reduced significantly. Using a commercial CFD (computational fluid dynamics) code FLUENT, the fluid flow and heat transfer of the fully diffuse plasma in an assumed magnetron torch have been simulated for qualitatively discussing the AMF effects. Numerical results show that the AMF significantly impels the plasma to retract axially and expand radially. As a result, the plasma intensity distribution on the cross section of the torch gets to be more uniform.     

Fluid Simulation of the Conduction Phase of the Plasma Erosion Opening Switch

The conduction phase of the plasma erosion opening switch (PEOS) is studied using a 1?-D electromagnetic two-fluid code. The focus of this work is on understanding how two effects, a current-limiting model of electron emission, and the magnetic insulation of electrons at the cathode, determine current conduction in the plasma. Simulations are performed in the parameter regimes of the Gamble I, POP, and PBFA II pulsed power generators, and previous low-density, short-rise time simulations of the PEOS. Fluid code results are compared to a 1-D analytic theory and to the Gamble I and POP experiments. Good agreement between theory and simulation, but mixed agreement between simulation and experiment is found. Experimental B-field measurements on POP show weaker j × B compression than the simulation. Current penetration and plasma current channels qualitatively similar to experimental observation are found in the Gamble I regime. However, magnetic insulation of electrons emitted from the cathode bunches the electron flow into narrower current channels than observed experimentally. In several cases, the presence of an electron-scattering or energy-loss mechanism near the cathode must be invoked to overcome magnetic insulation and widen the current channels.     

Implicit Collisional Three-Fluid Simulation of the Plasma Erosion Opening Switch

The plasma erosion opening switch (PEOS) has been studied with the aid of the ANTHEM implicit simulation code. This switch consists of fill plasma injected into a transmission line. The plasma is ultimately removed by self-electrical forces, permitting energy delivery to a load. Here, ANTHEM treats the ions and electrons of the fill plasma and the electrons emitted from the transmission-line cathode as three distinct Eulerian fluids-with electron inertia retained. This permits analysis of charge separation effects, and avoids the singularities that plague conventional MHD codes at low density. E and B fields are computed by the implicit moment method, allowing for time steps well in excess of the electron plasma period ?t >> ?p-1, and cells much wider than a Debye length, ?x >> ?D. Switch dynamics are modeled as a function of the driving electrical pulse characteristics, the fill plasma parameters, and the emission properties of the transmission line walls-for both collisionless and anomalously collisional electrons. Our low-fill-density (ne ? 4 × 1012 electrons/cm3) collisionless calculations are in accord with earlier particle code results. Our high-density computations (ne ? 2 × 1013 electrons/cm3) show the opening of the switch proceeding through both ion erosion and magnetic pressure effects. The addition of anomalous electron collisions is found to diffuse the driving B field into the fill plasma, producing broad current channels and reduced magnetic pressure effects, in some agreement with NRL experimental measurements.     

Most electron heat transport is not anomalous; it is a paleoclassical process in toroidal plasmas

It is hypothesized that radial electron heat transport in magnetically confined toroidal plasmas results from paleoclassical Coulomb collision processes (parallel electron heat conduction and magnetic field diffusion). In such plasmas the electron temperature is equilibrated along magnetic field lines a long length L (> poloidal periodicity length piR0q), which is the minimum of the electron collision length and an effective field line length. Thus, diffusing field lines induce a radial electron heat diffusivity M identical with L/(piR0q) approximately 10>1 times the magnetic field diffusivity eta/mu0 approximately nue(c/omegap)2.     

Electron diamagnetic effect on axial force in an expanding plasma: experiments and theory

The axial force imparted from a magnetically expanding current-free plasma is directly measured for three different experimental configurations and compared with a two-dimensional fluid theory. The force component solely resulting from the expanding field is directly measured and identified as an axial force produced by the azimuthal current due to an electron diamagnetic drift and the radial component of the magnetic field. The experimentally measured forces are well described by the theory.     

Dielectric Characteristics for Radio Frequency Waves in a Laboratory Dipole Plasma

Transverse and parallel dielectric permittivity elements have been derived for radio frequency waves in a laboratory dipole magnetic field plasma. Vlasov equation is resolved for both the trapped and untrapped particles as a boundary value problem to define their separate contributions to the dielectric tensor components. To estimate the wave power absorbed in the plasma volume the perturbed electric field and current density components are decomposed in a Fourier series over the poloidal angle. In this case, the dielectric characteristics can be analyzed independently of the solution of the Maxwell's equations. As usual, imaginary part of the parallel permittivity elements is necessary to estimate the electron Landau damping of radio frequency waves, whereas imaginary part of the transverse permittivity elements is important to estimate the wave dissipation by the cyclotron resonances. Computations of the imaginary part of the parallel permittivity elements are carried out in a wide range of the wave frequencies. (© 2006 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Vacuum Arc Plasma Centrfuge for Element and Isotope Separation

Vacuum arcs have been studied extensively in the past several decades with applications primarily in the areas of switching, vacuum remelting, and vapor deposition. Application of the vacuum arc for element and isotope separation has been studied recently and is reviewed in this paper. An arc was produced in a 30-cm-diameter 4-m-long cylindrical chamber with coaxially mounted electromagnets providing a 2.6-m-long constant axial magnetic field of up to 6 kG. The vacuum discharge between a solid cathode and a mesh anode was triggered electrically. A pulse-forming network (PFN) of 70-m? impedance provided nearly constant-current discharge pulses of several kiloamps and 6-12-ms duration. The magnetized plasma column, flowing axially from the anode with a typical velocity of 106 cm/s, rotated nearly as a solid body. This rotation was due to the E × B drift, produced by the axial magnetic field and the radial electric field across the column. A typical rotation frequency was 105 rad/s. The centrifugal effect due to the rotation caused a radial redistribution of ions within the plasma column, thereby producing elemental and isotope enrichment. The separation was observed to increase exponentially with the square of the radius. Enrichments of up to 300 percent were measured in a Cu-Zn plasma. The radial plasma density profile was found to be roughly Gaussian, with central electron densities of about 1013 cm-3. The radial potential profile across the column was measured and found to be parabolic with radius.     

脉冲驱动磁化等离子体内爆升温点火的数值模拟北大核心CSCD

通过建立物理模型,用一维三温拉氏磁流体力学程序分析了由强电流(MA)脉冲驱动的金属套筒内爆压缩磁化等离子体的升温点火及能量增益过程。分析了脉冲驱动的金属套筒内爆、不同驱动源对金属套筒内爆升温的影响、Z箍缩过程中内嵌磁场和预加热温度对磁化等离子体升温的影响,以及点火需要的初态参数和点火后的能量输出。此外,对该过程中磁场增加α粒子能量沉积、降低电子离子热传导能量损失的物理机制做了介绍和分析。磁流体数值模拟结果显示:当初始的内嵌磁场和燃料的预加热温度分别取5T和250eV时,即可获得超过4keV的升温,初始参数包括内嵌磁场、预加热温度、燃料密度、套筒尺度、驱动脉冲幅值、加载时间等。在一定的条件下,点火成功,可产生kT量级的强磁场,并获得百kJ/mm量级的能量输出。