Journal of Plasma Physics最新文献

The particle-in-cell (PIC) method is successfully used to study magnetized plasmas. However, this requires large computational costs and limits simulations to short physical run times and often to set-ups of less than three spatial dimensions. Traditionally, this is circumvented either via hybrid-PIC methods (adopting massless electrons) or via magneto-hydrodynamic-PIC methods (modelling the background plasma as a single charge-neutral magneto-hydrodynamical fluid). Because both methods preclude modelling important plasma-kinetic effects, we introduce a new fluid-PIC code that couples a fully explicit and charge-conserving multi-fluid solver to the PIC code SHARP through a current-coupling scheme and solve the full set of Maxwell's equations. This avoids simplifications typically adopted for Ohm's law and enables us to fully resolve the electron temporal and spatial scales while retaining the versatility of initializing any number of ion, electron or neutral species with arbitrary velocity distributions. The fluid solver includes closures emulating Landau damping so that we can account for this important kinetic process in our fluid species. Our fluid-PIC code is second-order accurate in space and time. The code is successfully validated against several test problems, including the stability and accuracy of shocks and the dispersion relation and damping rates of waves in unmagnetized and magnetized plasmas. It also matches growth rates and saturation levels of the gyro-scale and intermediate-scale instabilities driven by drifting charged particles in magnetized thermal background plasmas in comparison with linear theory and PIC simulations. This new fluid-SHARP code is specially designed for studying high-energy cosmic rays interacting with thermal plasmas over macroscopic time scales.

This paper explores the feasibility of a break-even-class mirror referred to as BEAM (break-even axisymmetric mirror): a neutral-beam-heated simple mirror capable of thermonuclear-grade parameters and $Qsim 1$ conditions. Compared with earlier mirror experiments in the 1980s, BEAM would have: higher-energy neutral beams, a larger and denser plasma at higher magnetic field, both an edge and a core and capabilities to address both magnetohydrodynamic and kinetic stability of the simple mirror in higher-temperature plasmas. Axisymmetry and high-field magnets make this possible at a modest scale enabling a short development time and lower capital cost. Such a $Qsim 1$ configuration will be useful as a fusion technology development platform, in which tritium handling, materials and blankets can be tested in a real fusion environment, and as a base for development of higher-$Q$ mirrors.

Electron-cyclotron resonance heating (ECRH) is the main heating mechanism in the Wendelstein 7-X (W7-X) stellarator. Although second-harmonic ECRH (X2) has been used routinely for plasma startup, startup at third harmonic (X3) is known to be much more difficult. In this work, we investigate the energy gain of particles during nonlinear wave–particle interaction for conditions relevant to second- and third-harmonic startups in W7-X. We take into account both the beam and the ambient magnetic field inhomogeneities. The latter is shown to significantly increase the mean energy gain resulting from a single wave–particle resonant interaction. In W7-X-like conditions, the improvement in maximum gained energy is up to 4 times the analogous uniform magnetic field case. However, this improvement is not enough to ensure X3 startup. The optimal magnetic field inhomogeneity length scale for average energy gain and start up in W7-X-like conditions is found to be in the range of $1$ to $3 {rm km}^{-1}$. A possibility of using multiple beams with neighbouring resonances is also considered. A considerable enhancement of the energy gain is demonstrated.

Beam–target reactions in fusion plasmas play an important role in both magnetic confinement fusion and inertial confinement fusion in the condition of low-density plasmas with high-velocity interactions. The traditional method for calculating beam–target reaction rate neglects the transport process of incident particles in inhomogeneous plasmas, leading to errors providing that the temperature and density in the transport path of incident particles vary obviously. An improved method considering the transport process is proposed in this paper to eliminate the deficiencies. Then the method is employed in high-speed plasma collision studies. When the initial plasma density and temperature are set to $0.5,{rm g},{rm cm}^{-3}$ and 100 eV, it is found that the beam–target reaction rate calculated by the traditional method is almost identical to that by our method if the collision velocity is less than 600 km s$^{-1}$. However, the traditional method is not suitable for study as the collision velocity gets higher, inducing obvious differences, which can reach 70 % at 1000 km s$^{-1}$. The improved method will make large corrections to evaluate the importance of the non-negligible beam–target reaction for inertial confinement fusion schemes with large implosion velocity such as double-cone ignition and impact ignition, in which the high-speed plasmas collide with each other to realize plasma ignition.

Our earth is immersed in the near-earth space plasma environment, which plays a vital role in protecting our planet against the solar-wind impact and influencing space activities. It is significant to investigate the physical processes dominating the environment, for deepening our scientific understanding of it and improving the ability to forecast the space weather. As a crucial part of the National Major Scientific and Technological Infrastructure–Space Environment Simulation Research Infrastructure (SESRI) in Harbin, the Space Plasma Environment Research Facility (SPERF) builds a system to replicate the near-earth space plasma environment in the laboratory. The system aims to simulate the three-dimensional (3-D) structure and processes of the terrestrial magnetosphere for the first time in the world, providing a unique platform to reveal the physics of the 3-D asymmetric magnetic reconnection relevant to the earth's magnetopause, wave–particle interaction in the earth's radiation belt, particles’ dynamics during the geomagnetic storm, etc. The paper will present the engineering design and construction of the near-earth space plasma simulation system of the SPERF, with a focus on the critical technologies that have been resolved to achieve the scientific goals. Meanwhile, the possible physical issues that can be studied based on the apparatus are sketched briefly. The earth-based system is of great value in understanding the space plasma environment and supporting space exploration.

The ability to quickly refresh gas-jet targets without cycling the vacuum chamber makes them a promising candidate for laser-accelerated ion experiments at high repetition rate. Here we present results from the first high repetition rate ion acceleration experiment on the VEGA-3 PW-class laser at CLPU. A near-critical density gas-jet target was produced by forcing a 1000 bar H$_2$ and He gas mix through bespoke supersonic shock nozzles. Proton energies up to 2 MeV were measured in the laser forward direction and 2.2 MeV transversally. He$^{2+}$ ions up to 5.8 MeV were also measured in the transverse direction. To help maintain a consistent gas density profile over many shots, nozzles were designed to produce a high-density shock at distances larger than 1 mm from the nozzle exit. We outline a procedure for optimizing the laser–gas interaction by translating the nozzle along the laser axis and using different nozzle materials. Several tens of laser interactions were performed with the same nozzle which demonstrates the potential usefulness of gas-jet targets as high repetition rate particle source.

The Hamiltonian dynamics of a single particle in a rotating plasma column, interacting with an magnetic multipole is perturbatively solved for up to second order, using the method of Lie transformations. First, the exact Hamiltonian is expressed in terms of canonical action-angle variables, and then an approximate integrable Hamiltonian is introduced, using another set of actions and angles, which describe the centre of oscillation for the particle. The perturbation introduces an effective ponderomotive potential, which to leading order is positive. At the second order, the pseudopotential consists of a sum of terms of the Miller form, and can have either sign. Additionally, at second order, the ponderomotive interaction introduces a modification to the particle effective mass, when considering the motion along the column axis. It is found that particles can be axially confined by the ponderomotive potentials, but acquire radial excursions which scale as the confining potential. The radial excursions of the particle along its trajectory are investigated, and a condition for the minimal rotation frequency for which the particle remains radially confined is derived. Last, we comment on the changes to the aforementioned solution to the pseudopotentials and particle trajectory in the case of resonant motion, that is, a motion which has the same periodicity as the perturbation.

Angular momentum coupling between a rotating magnetized plasma and torsional Alfvén waves carrying orbital angular momentum (OAM) is examined. It is demonstrated not only that rotation is the source of Fresnel–Faraday rotation – or orbital Faraday rotation effects – for OAM-carrying Alfvén waves, but also that angular momentum from an OAM-carrying Alfvén wave can be transferred to a rotating plasma through the inverse process. For the direct process, the transverse structure angular rotation frequency is derived by considering the dispersion relation for modes with opposite OAM content. For the inverse process, the torque exerted on the plasma is derived as a function of wave and plasma parameters.

Fishbone bursts have been observed to strongly correlate to internal transport barrier (ITB) formation in a number of tokamak devices. A simple model incorporating the fishbone dynamics and ion pressure gradient evolution is proposed in order to investigate the key physics parameters assisting the triggering of ITB. The time evolution of fishbone is described by the well-known predator–prey model. For each burst cycle, the energetic particles (EPs) resonantly interact with fishbone and are radially expelled from inner region leading to a radial current. A compensating bulk plasma return current and, hence, poloidal flow can be induced if the fishbone cycle frequency is greater than the poloidal flow damping rate. When the shear of the poloidal flow exceeds a critical value, the turbulent fluctuations are suppressed and the bulk ion pressure gradient transits to the high-confinement state. It is shown that this process is only sensitive to the deposition rate of the trapped EPs within the $q=1$ surface, but not sensitive to other parameters. A quantitative formula for the shearing rate of poloidal flow induced by fishbone bursts is derived and verified numerically.