In this work, a jet of cold plasma is generated in a device supplied in helium and powered with a high-voltage nanopulse power supply, hence generating guided streamers. We focus on the interaction between these guided streamers and two targets placed in a series: a metal mesh target (MM) at floating potential followed by a metal plate target (MP) grounded by a 1500 Ω resistor. We demonstrate that such an experimental setup allows to shift from a physics of streamer repeatability to a physics of streamer self-organization, i.e., from the repetition of guided streamers that exhibit fixed spatiotemporal constants to the emergence of self-organized guided streamers, each of which is generated on the rising edge of a high-voltage pulse. Up to five positive guided streamers can be self-organized one after the other, all distinct in space and time. While self-organization occurs in the capillary and up to the MM target, we also demonstrate the existence of transient emissive phenomena in the inter-target region, especially a filamentary discharge whose generation is directly correlated with complexity order Ω. The mechanisms of the self-organized guided streamers are deciphered by correlating their optical and electrical properties measured by fast ICCD camera and current-voltage probes, respectively. For the sake of clarity, special attention is paid to the case where three self-organized guided streamers (α, \b{eta} and γ) propagate at vα = 75.7 km.s-1, v\b{eta} = 66.5 km.s-1 and vγ = 58.2 km.s-1), before being accelerated in the vicinity of the MM target.
Fabio Veronese, Piero Agostinetti, Christian Hopf
et al.
The Beam Driven Plasma Neutralizer (BDPN) has been proposed as a more efficient alternative to the gas neutralizer for negative-ion based Neutral Beam Injection (NNBI). In this paper we model the performance of an entire NNBI beamline with a BDPN. We simultaneously consider all the relevant physics and engineering aspects, the most important being the plasma density and degree of ionization inside the BDPN as a function of its geometry and feed gas flow, the geometrical transmission of the beamline, the dependence of the neutral gas distribution in the beamline on the geometry of the beamline components and gas flows, and the species evolution of the extracted D$^-$ beam through this neutral and charged particle distribution. Furthermore, we calculate the heat loads expected on the BDPN parts and on the NBI components located downstream of it and study the effect of the magnetic cusp field across the BDPN entrance on beamline transmission. While our results constitute an optimization only under the applied boundary conditions, we find that the beamline with a BDPN increases the system's wall plug efficiency by about 13% to 0.34 from the 0.30 estimated for a gas neutralizer.
Particle-in-cell codes usually represent large groups of particles as a single macroparticle. These codes are computationally efficient but lose information about the internal structure of the macroparticle. To improve the accuracy of these codes, this work presents a method in which, as well as tracking the macroparticle, the moments of the macroparticle are also tracked. Although the equations needed to track these moments are known, the coordinate transformations for moments where the space and time coordinates are mixed cannot be calculated using the standard method for representing moments. These coordinate transformations are important in astrophysical plasma, where there is no preferred coordinate system. This work uses the language of Schwartz distributions to calculate the coordinate transformations of moments. Both the moment tracking and coordinate transformation equations are tested by modelling the motion of uncharged particles in a circular orbit around a black hole in both Schwarzschild and Kruskal-Szekeres coordinates. Numerical testing shows that the error in tracking moments is small, and scales quadratically. This error can be improved by including higher order moments. By choosing an appropriate method for using these moments to deposit the charge back onto the grid, a full particle-in-cell code can be developed.
Magnetohydrodynamics (MHD) is a subject concerned with the dynamics of electrically conducting fluids (plasma) and can be applied in electric power generation. As a unique technology for producing direct-current electricity without moving parts, it can be utilized within a high-temperature topping power cycle to be combined with a traditional bottoming power cycle, forming a combined-cycle MHD system. This study presents governing equations for the electric field and current density field within a moving plasma subject to an applied magnetic field. The modeling equations are described at four descending levels of complexity. Starting with the first level of modeling, which is the most general case, where no assumptions are made regarding the electric field, plasma velocity field, applied magnetic field, or electrode geometry. In the second level of modeling, the magnetic field is treated as one-dimensional. In the third level of modeling, a specific Faraday-type magnetohydrodynamics plasma generator channel is considered, having two continuous electrodes acting as parallel constant-voltage terminals. In the fourth (and simplest) level of modeling, an additional approximation is made by setting the Hall parameter to zero and replacing all vector fields with scalar quantities. For that simplest model, a representative set of operation conditions (electric conductivity 20 S/m, temperature 2800 K, supersonic plasma gas speed 2000 m/s with Mach 2.134, and magnetic flux density 5 T) shows that the optimum idealized electric power that can be extracted from a unit volume of plasma is estimated as 500 MW/m3. This is a much larger volumetric power density than typical values encountered in reciprocating piston-type engines (0.2 MW/m3) or rotary gas turbine engines (0.5 MW/m3). Such an extremely high power density enables very compact power generation units.
We apply inverse design methods to produce two-dimensional triangular-lattice plasma metamaterial (PMM) devices which are then constructed and demonstrated experimentally. Finite difference frequency domain simulations are used along with forward-mode automatic differentiation to optimize the plasma densities of each of the plasma elements in the PMM to perform beam steering and demultiplexing under transverse magnetic polarization. The optimal device parameters are then used to assign plasma density values to elements that make up an experimental version of the device. Device performance is evaluated against both the simulated results and human-designed alternatives, showing the benefits and disadvantages of in-silico inverse design and paving the way for future fully in-situ optimization.
This research finds an equation in a continuous domain and a discrete equation governing the system of cold bosonic atoms (CBA) in a zigzag optical lattice using a continuum approximation. Many solutions to the equation were obtained using two distinct methods: the three-wave approach (multi-wave interaction, rational solutions, and rational solution interaction) and the extended sub-equation method. These analytical approaches are more effective, consistent, and comprehensive mathematical tools for obtaining various exact closed-form solutions for a wide range of fractional space-time nonlinear evolution equations encountered in optical physics, condensed matter physics, and plasma physics. The solutions generated are in the form of hyperbolic and trigonometric solutions, and other-form solutions are obtained. Three-dimensional graphics and contour plots are often used to depict the graphical representations of the combined soliton solutions. These findings will aid our understanding of the dynamics of the zigzag optical grids and many other structures formed by colder bosonic atoms. The applied approaches are more simple, efficient, and straightforward to obtain the closed-form solutions for various nonlinear evolution equations in the fields of nonlinear sciences and physical engineering.
Growing attention has been paid to nonthermal plasma treatment technology and its effects on the degradation of organic matter, especially for antibiotics. However, the majority of the conducted research has focused on the experimental results. Rare attempts were made to analyze the reaction mechanism at the microscopic level. In this paper, molecular dynamics simulation and reactive forcefields were used to investigate the reaction mechanism of different plasma particle interactions with azithromycin molecules. The simulation results indicated that the degradation of azithromycin was caused by the destruction of C-H and C-C bonds, followed by the formation of C=C and C=O bonds when reacted with the active particles. It was also found that the ability of degrading azithromycin varied among the different types of active particles. The oxygen atoms had the strongest ability to decompose the azithromycin molecule, with 38.61% of the C-H bonds broken as compared with other oxygenated species. The findings from this computational simulation could provide theoretical support and guidance for subsequent practical experiments.
The ITER Diagnostic Residual Gas Analyzer (DRGA) will measure the distribution of gas species, i.e., deuterium (D), tritium (T), and impurities, in the divertor exhaust stream and in the plasma periphery, with time resolution relevant to fusion plasma–wall particle dynamics. The uniqueness of the DRGA, over previous implementations of plasma dynamics residual gas analysis, is an integrated approach, combining mass and low-temperature plasma-activated optical spectroscopy, in a differentially pumped analysis station. A further unique feature of the ITER divertor-specific DRGA is an ~8-m separation of the analysis station from the sampled pumping duct, while still maintaining a ~1-s response time for hydrogen isotopic concentrations. ITER DRGA final design activities are strongly benefiting from testing of prototypical DRGA components and methods on present fusion devices, most currently on JET and W7-X. DRGA systems are implemented on both these devices and include sensors (and pumping methods) that are directly relevant to the ITER DRGA design. The recent JET-DTE2 campaign has provided the first experience on operating the combined ITER DRGA sensors with D-T plasmas. While enhancing system design for ITER, this experience on operating devices has also revealed additional engineering challenges, which further guide the continuing final design project. Meanwhile, the recent determination that the ITER DRGA, with slight optimization, will resolve the helium isotopes well enough to support an ITER pre-DT, He-3-based heating scheme, has greatly increased ITER Research Program interest in the DRGA and its implementation well ahead of the DT phase.
Magnetohydrodynamic (MHD) simulations of the solar corona have become more popular with the increased availability of computational power. Modern computational plasma codes, relying upon computational fluid dynamics (CFD) methods, allow the coronal features to be resolved using solar surface magnetograms as inputs. These computations are carried out in a full three-dimensional domain and, thus, selection of the correct mesh configuration is essential to save computational resources and enable/speed up convergence. In addition, it has been observed that for MHD simulations close to the hydrostatic equilibrium, spurious numerical artefacts might appear in the solution following the mesh structure, which makes the selection of the grid also a concern for accuracy. The purpose of this paper is to discuss and trade off two main mesh topologies when applied to global solar corona simulations using the unstructured ideal MHD solver from the COOLFluiD platform. The first topology is based on the geodesic polyhedron and the second on $UV$ mapping. Focus is placed on aspects such as mesh adaptability, resolution distribution, resulting spurious numerical fluxes and convergence performance. For this purpose, first a rotating dipole case is investigated, followed by two simulations using real magnetograms from the solar minima (1995) and solar maxima (1999). It is concluded that the most appropriate mesh topology for the simulation depends on several factors, such as the accuracy requirements, the presence of features near the polar regions and/or strong features in the flow field in general. If convergence is of concern and the simulation contains strong dynamics, then grids which are based on the geodesic polyhedron are recommended compared with more conventionally used $UV$-mapped meshes.
The phenomenon of cathode spot (CS) retrograde motion has been attracting extensive research in the field of vacuum arc. Although a number of theoretical models have been proposed to explain this phenomenon, up to now, the first-principle understanding is still missing. In this work, based on the generalized Ohm’s law, it is proposed that there is a transport current in the explosive emission center when a transverse magnetic field (TMF) is applied. The investigation of the influence of the transport current on the strength of the electric field in the cathode sheath suggests that the transport current affects the ignition of new CS in the retrograde direction. In the framework of the ecton model, combined with the dynamics of explosive emission plasma, the relationship between the strength of applied TMF and the ignition probability of new CS in different directions is established, and the retrograde motion characteristics of CS are simulated using a statistical model. The simulation results show that the directionality of retrograde motion of CS is enhanced with an increase in TMF, and the retrograde motion speed of CS increases with the increase of the applied TMF. When the strength of TMF is small, the motion speed of a single CS has a linear relationship with the strength of TMF, but with the increase of the strength of TMF, the motion speed of a single CS is gradually saturated, which is consistent with experimental facts. The feasibility of the model is verified by comparing the simulation results with relevant experimental results.
The magnetosphere–ionosphere dynamics comprises processes both directly related to solar wind variability and of purely internal origin. The latter represent a huge drawback for correctly forecasting the magnetosphere–ionosphere dynamics during geomagnetic storms and substorms. Here, we use wavelet analysis to further characterize the storm–substorm relationship through the use of the AL and SYM-H geomagnetic indices. We focus our analysis on one of the strongest geomagnetic storms of solar cycle 23 that occurred on 20 November 2003. Our findings suggest that, during disturbed periods, a significant amount of information comes from the interactions between geomagnetic storms and magnetospheric substorms. Thus, predicting the intensity and the duration of a geomagnetic storm requires information coming not only from the solar wind variability but also from the nonlinear variability of the magnetosphere–ionosphere system occurring on short time scales. Our results are also discussed in the framework of Space Weather, suggesting an extended use of non-traditional dynamical systems approaches (such as those based on extreme value statistics and tipping point analysis) to deal with emergent behaviours coming from different sources during geomagnetic storms and magnetospheric substorms.
Today, partially magnetized low-temperature plasmas (LTP) in an ExB configuration, where the applied magnetic field is perpendicular to the self-consistent electric field, have important industrial applications. Hall thrusters, a type of electrostatic plasma propulsion, are one of the main LTP technologies whose advancement is hindered by the not-fully-understood underlying physics of operation, particularly, with respect to the plasma instabilities and the associated electron cross-field transport. The development of Hall thrusters with unconventional magnetic field topologies has imposed further questions regarding the instabilities' characteristics and the electrons' dynamics in these modern cross-field configurations. Accordingly, we present in this effort a series of studies on the influence of four factors on the plasma processes in the radial-azimuthal coordinates of a Hall thruster, namely, the magnetic field gradient, Secondary Electron Emission, electron-neutral collisions, and plasma number density. The studies are carried out using the reduced-order particle-in-cell (PIC) code developed by the authors. The setup of the radial-azimuthal simulations largely follows a well-defined benchmark case from the literature in which the magnetic field is oriented along the radius and a constant axial electric field is applied perpendicular to the simulation plane. The salient finding from our investigations is that, in the studied cases corresponding to elevated plasma densities, an inverse energy cascade leads to the formation of a long-wavelength, high-frequency azimuthal mode. Moreover, in the presence of strong magnetic field gradients, this mode is fully developed and induces a significant electron cross-field transport as well as a notable heating of the ion population.
Relevant uncertainties on theoretical atomic data are vital to determine the accuracy of plasma diagnostics in a number of areas including in particular the astrophysical study. We present a new calculation of the uncertainties on the present theoretical ion-impact charge exchange atomic data and X-ray spectra based on a set of comparisons with the existing laboratory data obtained in historical merged-beam, cold-target recoil-ion momentum spectroscopy, and electron beam ion traps experiments. The average systematic uncertainties are found to be 35-88% on the total cross sections, and 57-75% on the characteristic line ratios. The model deviation increases as the collision energy decreases. The errors on total cross sections further induce a significant uncertainty to the calculation of ionization balance for low temperature collisional plasmas. Substantial improvements of the atomic database and dedicated laboratory measurements are needed to get the current models ready for the X-ray spectra from the next X-ray spectroscopic mission.
Superconducting transformer, which used for power supply of the large current superconductor, is faced with the conductor sample current attenuation problem since the existence of joint resistance. A solution is proposed to solve this attenuation problem to make the conductor sample current controllable. The adaptive Proportion Integration Differentiation (PID) control strategy based on learning rates (η) Radial Basis Function (RBF) neural network (NN) is an improvement over the conventional PID algorithm. The algorithm is applied to obtain the optimal parameters for the PID controller to the working variations arising from transformer nonlinear dynamics. The controller system is carried out on the superconducting conductor test platform of the Institute of Plasma Physics, Chinese Academy of Sciences (ASIPP). The results showed that the overall deviation of the system was within 0.2%. The performance of the proposed approach is validated by superconductor experiments under practical conditions.
A complete 28-GHz/50-kW gyrotron system was developed by the Institute of Applied Electronics, China Academy of Engineering Physics, Mianyang, China, as a plasma heating source for the EXL50 fusion device, a compact spherical tokamak constructed by the ENN Group in 2019. The gyrotron employs a triode magnetron injection gun, a TE02-mode cylindrical cavity, a built-in quasi-optical mode converter, a single-stage depressed collector, and a single-disk boron nitride window. The gyrotron system features a match optical unit, a superconducting magnet, a calorimetric dummy load, a high-voltage power supply system, and other auxiliary components, was delivered to the EXL50 site. In the acceptance test, five successive shots of each power over 50 kW were generated for 30 s at 28 GHz. An average power of 51 kW, a maximum power of 55 kW, and 46% overall efficiency were obtained. The gyrotron system was assembled in the EXL50 and applied to a plasma discharge experiment. The gyrotron is the first practical domestic gyrotron produced for fusion applications in China. The design and test results of this instrument are presented in this article.
Miguel Jiménez-Redondo, Audrey Chatain, Olivier Guaitella
et al.
In this work, we present the results of simulations carried out for N2-H2 capacitively coupled radio-frequency discharges, running at low pressure (0.3-0.9 mbar), low power (5-20 W), and for amounts of H2 up to 5 pct. Simulations are performed using a hybrid code that couples a two-dimensional time-dependent fluid module, describing the dynamics of the charged particles in the discharge, to a zero-dimensional kinetic module, that solves the Boltzmann equation and describes the production and destruction of neutral species. The model accounts for the production of several vibrationally and electronic excited states, and contains a detailed surface chemistry that includes recombination processes and the production of NHx molecules. The results obtained highlight the relevance of the interactions between plasma and surface, given the role of the secondary electron emission in the electrical parameters of the discharge and the critical importance of the surface production of ammonia to the neutral and ionic chemistry of the discharge.
Abstract Electroporation processes affect the permeability of cell membranes, which can be utilized for the delivery of plasma species in cancer therapy. By means of computational dynamics, many aspects of membrane electroporation have been unveiled at the atomic level for lipid membranes. Herein, a molecular dynamics simulation study was performed on native and oxidized membrane systems with transversal electric fields. The simulation result shows that the applied electric field mainly affects the membrane properties so that electroporation takes place and these pores are lined by hydrophilic headgroups of the lipid components. The calculated hydrophobic thickness, lateral diffusion and pair correlation revealed the role of 5α-CH in creation of water-pore in an oxidized membrane. Additionally, the permeability of reactive oxygen species was examined through these electroporated systems. The permeability study suggested that water pores in the membrane facilitate the penetration of these species across the membrane to the interior of the cell. These findings may have significance in experimental applications in vivo as once the reactive oxygen species reaches the interior of the cell, they may cause oxidative stress and induce apoptosis. Communicated by Ramaswamy H. Sarma
Protoplanetary discs are made of gas and dust orbiting a young star. They are also the birth place of planetary systems, which motivates a large amount of observational and theoretical research. In these lecture notes, I present a review of the magnetic mechanisms applied to the outer regions ($R\gtrsim 1\ \mathrm {AU}$) of these discs, which are the planet-formation regions. In contrast to usual astrophysical plasmas, the gas in these regions is noticeably cold ($T < 300\ \mathrm {K}$) and dense, which implies a very low ionisation fraction close to the disc midplane. In these notes, I deliberately ignore the innermost $(R\sim 0.1\ \mathrm {AU})$ region, which is influenced by the star–disc interaction and various radiative effects. I start by presenting a short overview of the observational evidence for the dynamics of these objects. I then introduce the methods and approximations used to model these plasmas, including non-ideal magnetohydrodynamics, and the uncertainties associated with this approach. In this framework, I explain how the global dynamics of these discs is modelled, and I present a stability analysis of this plasma in the local approximation, introducing the non-ideal magneto-rotational instability. Following this mostly analytical part, I discuss numerical models that have been used to describe the saturation mechanisms of this instability, and the formation of large-scale structures by various saturation mechanisms. Finally, I show that local numerical models are insufficient because magnetised winds are also emitted from the surface of these objects. After a short introduction on wind physics, I present global models of protoplanetary discs, including both a large-scale wind and the non-ideal dynamics of the disc.
These lecture notes are based on a tutorial given in 2017 at a plasma physics winter school in Les Houches. Their aim is to provide a self-contained graduate-student level introduction to the theory and modelling of the dynamo effect in turbulent fluids and plasmas, blended with a review of current research in the field. The primary focus is on the physical and mathematical concepts underlying different (turbulent) branches of dynamo theory, with some astrophysical, geophysical and experimental context disseminated throughout the document. The text begins with an introduction to the rationale, observational and historical roots of the subject, and to the basic concepts of magnetohydrodynamics relevant to dynamo theory. The next two sections discuss the fundamental phenomenological and mathematical aspects of (linear and nonlinear) small- and large-scale MHD dynamos. These sections are complemented by an overview of a selection of current active research topics in the field, including the numerical modelling of the geo- and solar dynamos, shear dynamos driven by turbulence with zero net helicity, and MHD-instability-driven dynamos such as the magnetorotational dynamo. The difficult problem of a unified, self-consistent statistical treatment of small and large-scale dynamos at large magnetic Reynolds numbers is also discussed throughout the text. Finally, an excursion is made into the relatively new but increasingly popular realm of magnetic-field generation in weakly-collisional plasmas. A short discussion of the outlook and challenges for the future of the field concludes the presentation.