Plasma Fluid Model Research Articles

Multi-stage electroporation dynamics induced by ionization wave–substrate interaction in cold plasma jet

The reactive species-independent nature of cold plasma's electric field is pivotal in biomedical applications. This work proposes to connect the plasma fluid model and the asymptotic Smoluchowski model for electroporation, providing a unified framework to investigate the evolution of the electric field in the biological substrate and the multi-stage electroporation response of the human cell. Two common substrates with distinct dielectric properties, namely, the cultivation medium and epidermis, are selected to report three stages of ionization wave (IW)–substrate interaction. The three-stage streamer discharge dynamics (restrike, axial-radial transition, and radial expansion of IW) induce three-stage cell electroporation dynamics (slow charging, fast charging, and electroporation), though the two processes are asynchronous. Specifically, the inner membrane covered the cell nucleus with ultra-short charging time that undergoes only the first two discharge stages in both substrates. Whether the cell membrane is exposed to the third stage of discharge depends on the permittivity of the substrate. The asynchrony can be attributed to the difference in the charging time of the cell membranes and substrates affected by the substrate permittivity. The presented model can provide quantitative insights into the cell electroporation induced by the IW–substrate interaction and theoretical guidance for plasma biomedical applications.

Numerical study on positive streamer in parallel-rod dielectric barrier discharge in atmospheric air

In this work, a parallel-rod dielectric barrier discharge (DBD) operating in atmospheric air is investigated through the two-dimensional plasma fluid model. The effects of applied voltage (Vp), secondary electron emission coefficient (γ), and photoionization are examined. Photoionization can significantly influence streamer dynamics by accelerating and broadening both volumetric and surface streamers and enhance the impact of the applied voltage. Without photoionization, the propagation distance of the surface streamer along the curved dielectric surface is limited to 0.1–0.2 mm under applied voltages of 8–8.5 kV. In contrast, with photoionization, this distance can extend to 0.3–0.6 mm. Achieving the same distance requires much higher voltages (10–11 kV) if without photoionization. The “double-layer” structure of the surface streamer is investigated, revealing that γ predominantly affects the surface branch with little impact on the volumetric branch. The critical charge density for streamer onset is found to be about 1018 m−3, and the volume-to-surface streamer transition is attributed to the lateral electric field provided by the space charges. This work provides insights into the regulation strategies and underlying mechanisms of streamer dynamics in parallel-rod DBDs in atmospheric air.

Numerical Simulation of Corona Discharge Plasma Affecting the Surface Behavior of Polymer Insulators

Corona discharge is a significant problem in the operation of high-voltage transmission and distribution systems, particularly for polymer insulators. Numerical simulation has become an effective tool for investigating the underlying physical mechanisms and optimizing the design of insulators. In this paper, we present a two-dimensional numerical simulation study on corona discharge plasma affecting the surface behavior of polymer insulators. The simulation was performed with the Comsol Multiphysic software and is based on the finite element method and the fluid plasma model, which considers ionization, recombination, and the transport of plasma species. The numerical results are analyzed to study the spatial and temporal characteristics of the corona discharge and its effect on the surface behavior of polymer insulators. The results show that the electric field is affected not only by the volume charge density but also by the surface charge density, which in turn depends on the densities of the charge carriers migrating on the insulator surface. However, the electric field drops drastically when one or two grading rings are installed. But one grading ring is not enough to limit the discharge.

Asymptotic-Preserving Discretization of Three-Dimensional Plasma Fluid Models

Asymptotic-Preserving Discretization of Three-Dimensional Plasma Fluid Models

Computational study of a helium-propellant microwave electrothermal thruster

The microwave electrothermal thruster (MET) utilizes wave-exited microdischarges to heat gas flows, enhancing the specific impulse of the thruster. Our computational study investigates a 17.5 GHz helium-propellant MET, employing a two-dimensional, axisymmetric fluid model of plasma coupled with electromagnetic wave and gas flows. The discharges operate in the glow regime, remaining weakly ionized, and in thermal non-equilibrium. The plasma densities reach approximately 1020m−3, and the gas temperature is around 2000 K. Even a slight off-resonant frequency operation results in a significantly lower plasma density and gas temperature. Gas heating, primarily driven by electromagnetic Joule heating, plays a critical role in influencing the overall discharge behavior. The measured peak thrust and specific impulse are 8.24 mN and 292 s, respectively, at a mass flow rate of 3.2 mg/s with 30 W of power. Compared to a cold gas thruster, there is a significant increase in the specific impulse by a factor of approximately 1.7. The enhanced performance trades off with propulsive efficiency, which decreases by a factor of 1.5 from the peak 65% at 10 W. This is due to higher energy losses to cavity walls from heat conduction with increased power. These findings underscore the critical balance between the input power and mass flow rate to improve the MET performance, highlighting the importance of power management to maximize thrust and efficiency. Furthermore, the predicted thrust and specific impulse agree well with experimental values for nominally similar MET thruster studies in the literature.

Self-pulsing of dielectric barrier discharges at low driving frequencies

This paper investigates the self-pulsing of dielectric barrier discharges (DBDs) at low driving frequencies. In particular, (a) the dependence of current on the product pd of gas pressure p and the gas gap length d, (b) the effects of lossy dielectrics (in resistive discharges) and large dielectric permittivity (in ferroelectrics) on current dynamics, (c) the transition from Townsend to a dynamic capacitively coupled plasma (CCP) discharge with changing pd values, and (d) the transition from Townsend to a high-frequency CCP regime with increasing the driving frequency. A one-dimensional fluid model of argon plasma is coupled to an equivalent RC circuit for lossy dielectrics. Our results show multiple current pulses per AC period in Townsend and CCP discharge modes which are explained by uncoupled electron–ion transport in the absence of quasineutrality and surface charge deposition at dielectric interfaces. The number of current pulses decreases with an increasing applied frequency when the Townsend discharge transforms into the CCP discharge. The resistive barrier discharge with lossy dielectrics exhibits Townsend and glow modes for the same pd value (7.6 Torr cm) for higher and lower resistances, respectively. Finally, we show that ferroelectric materials can amplify discharge current in DBDs. Similarities between current pulsing in DBD, Trichel pulses in corona discharges, and subnormal oscillations in DC discharges are discussed.

Self Consistent Modeling of Relativistic Runaway Electron Beams Giving Rise to Terrestrial Gamma‐Rays Flashes

AbstractTerrestrial Gamma Ray Flashes (TGFs) are short bursts of gamma rays occurring during thunderstorms. They are believed to be produced by Relativistic Runaway Electron Avalanches (RREAs). In this paper, we present a new numerical model based on the Particle‐In‐Cell (PIC) method to simulate the interactions between the electromagnetic fields and the electron avalanche self‐consistently. The code uses a cylindrical Yee lattice to numerically solve the electromagnetic fields, a Monte Carlo approach to simulate collisions with air molecules, and a plasma fluid model to calculate the effects of low‐energy electrons and ions. The model is first tested through the reproduction of dispersion relations in a hot plasma. RREAs propagating under a homogeneous background electric field are then simulated. Owing to the self‐consistent nature of this description, we report here new physical properties such as saturation processes in the electron density and in the number of high‐energy electrons, detailed dynamical screening of the electric field in the ion trail of the avalanche, and the associated electric currents. We find that the saturation of RREAs is obtained when the numbers of high‐energy electrons and photons is consistent with those believed to be representative of TGF sources.

Correction to: Coupled nonlinear Schrödinger (CNLS) equations for two interacting electrostatic wavepackets in a non-Maxwellian fluid plasma model

Correction to: Coupled nonlinear Schrödinger (CNLS) equations for two interacting electrostatic wavepackets in a non-Maxwellian fluid plasma model

Deicing and status characteristics of dual-side pulsed surface dielectric barrier discharge

The deicing process and its status characteristics of dual-side pulsed surface dielectric barrier discharge (SDBD) are studied via electro-optical diagnostics, thermal properties, and numerical simulation. Experimental results show that the dual-side pulsed SDBD can remove the glaze ice compared to the traditional pulsed SDBD under the applied pulse voltage of 8 kV and a pulse frequency of 1 kHz. The maximal temperature of dual-side pulsed SDBD reaches 39.5 °C under the discharge time of 800 s, while the maximal temperature of traditional pulsed SDBD is still below ice point about −7.8 °C. Surface temperatures of dual-side pulsed SDBD demonstrate that the SDBD with a gap of 1 mm possesses prospects in deicing. The maximal surface temperature reaches 37.1 °C under the pulse of 8 kV after the discharge time of 90 s. Focusing on the thermal effect, a two-dimensional plasma fluid model is implemented, and the results also indicate that the dual-side pulsed SDBD with a gap of 1 mm produces a highest heat density among the three different configurations. Comparing the spatial-temporal evolutions of plasma on both dielectric sides, primary positive streamer has a longer propagation length of 8.6 mm than the secondary negative streamer, the primary negative streamer, and the secondary positive streamer, which induces a long heat covered area. Four stages of deicing process are analyzed through a series of electrical parameters under different covered ice conditions.

Dynamic evolutions of Bohm sheaths and pre-sheaths

The time evolutions of the Bohm sheath and the related pre-sheath are analyzed as an initial value problem. The standard classical fluid model of a collisionless plasma is used with cold ions and Boltzmann distributed electrons. Numerical solutions of the basic equations show that a stationary plasma sheath itself is established within a few ion plasma periods. It is demonstrated analytically that for infinitely extended homogeneous plasmas in one spatial dimension, the only possible pre-sheath solution is dynamic, giving a steady expansion: no stationary solution exists for this one-dimensional case. The potential and density drops along the dynamic pre-sheath differ notably from the values found for stationary models suggested in the literature. Collisions give only formal changes to the collisionless results, and their substance remains the same in one spatial dimension. Cylindrical and spherical geometries, on the other hand, allow for physically acceptable, fully nonlinear, stationary solutions with analytical results given. These are supported by numerical solutions.

Physics-informed neural networks based on source term decoupled and its application in discharge plasma simulation

In recent years, the artificial intelligence computing paradigm represented by physics-informed neural networks (PINNs) has received great attention in the field of plasma numerical simulation. However, the plasma chemical system considered in related research is relatively simplified, and the research on solving the more complex multi-particle low-temperature fluid model based on PINNs is still blank. In more complex chemical systems, the coupling relationship between particle densities and between particle densities and mean electron energy become more intricate. Therefore, the applicability of PINNs in dealing with sophisticated reaction systems needs further exploring and improving. In this work, we propose a general PINN framework (source term decoupled PINNs, Std-PINNs) for solving multi-particle low-temperature plasma fluid model. By introducing equivalent positive ions and replacing each particle transport equation with the current continuity equation as a physical constraint, Std-PINN splits the entire solution process into the training processes of two neural networks, realizing the decoupling of the source term of the heavy particle transport equation from the electron density and mean electron energy, which greatly reduces the complexity of neural network training. In this work, the application of Std-PINNs to solving multi-particle low-temperature plasma fluid models is demonstrated through two classic discharge cases with different complexity of reaction systems (low-pressure argon glow discharge and atmospheric-pressure helium glow discharge) and the performance of Std-PINN is compared with that of conventional PINN and finite element method (FEM). The results show that the training results output from the traditional PINN are completely incorrect due to the strong coupling correlation of each physical variable through the source terms of each particle transport equation, while the <i>L</i><sub>2</sub> relative error between Std-PINN and FEM results can reach up to ~10<sup>–2</sup> , thus verifying the feasibility of Std-PINN in simulating multi-particle plasma fluid model. Std-PINN expands the application of deep learning method to modeling complex physical systems and provides new ideas for conducting low-temperature plasma simulations. In addition, this study provides novel insights into the field of artificial intelligence scientific computing: the mathematical form that describes the state of a physical system is not unique. By introducing equivalent physical variables, equations suitable for neural network solutions can be derived and combined with observable data to simplify problems.

Coupled nonlinear Schrödinger (CNLS) equations for two interacting electrostatic wavepackets in a non-Maxwellian fluid plasma model

The nonlinear dynamics of two co-propagating electrostatic wavepackets in a one-dimensional non-magnetized plasma fluid model is considered, from first principles. The coupled waves are characterized by different (carrier) wavenumbers and amplitudes. A plasma consisting of non-thermalized (κ\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\kappa $$\\end{document}-distributed) electrons evolving against a cold (stationary) ion background is considered. The original model is reduced, by means of a multiple-scale perturbation method, to a pair of coupled nonlinear Schrödinger (CNLS) equations for the dynamics of the wavepacket envelopes. For arbitrary wavenumbers, the resulting CNLS equations exhibit no known symmetry and thus intrinsically differ from the Manakov system, in general. Exact analytical expressions have been derived for the dispersion, self-modulation (nonlinearity) and cross-modulation (coupling) coefficients involved in the CNLS equations, as functions of the wavenumbers (k1\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$k_1$$\\end{document}, k2\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$k_2$$\\end{document}) and of the spectral index κ\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\kappa $$\\end{document} characterizing the electron profile. An analytical investigation has thus been carried out of the modulational instability (MI) properties of this pair of wavepackets, focusing on the role of the intrinsic (variable) parameters. Modulational instability is shown to occur in most parts of the parameter space. The instability window(s) and the corresponding growth rate are calculated numerically in a number of case studies. Two-wave interaction favors MI by extending its range of occurrence and by enhancing its growth rate. Growth rate patterns obtained for different κ\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\kappa $$\\end{document} index (values) suggest that deviation from thermal (Maxwellian) equilibrium, for low κ\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\kappa $$\\end{document} values, leads to enhance MI of the interacting wave pair. To the best of our knowledge, the dynamics of two co-propagating wavepackets in a plasma described by a fluid model with κ\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\kappa $$\\end{document}-distributed electrons is investigated thoroughly with respect to their MI properties as a function of κ\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\kappa $$\\end{document} for the first time, in the framework of an asymmetric CNLS system whose coefficients present no obvious symmetries for arbitrary k1\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$k_1$$\\end{document} and k2\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$k_2$$\\end{document}. Although we have focused on electrostatic wavepacket propagation in non-thermal (non-Maxwellian) plasma, the results of this study are generic and may be used as basis to model energy localization in nonlinear optics, in hydrodynamics or in dispersive media with Kerr-type nonlinearities where MI is relevant.

Gas heating by inductively coupled low-pressure chlorine process plasmas

The mechanism of gas heating in low-pressure inductively coupled chlorine plasma is analyzed using a self-consistent two-dimensional axisymmetric fluid plasma model that is coupled with the compressible Navier–Stokes equations. For gas pressures of 10 and 20 mTorr and the discharge power in the range 0.1–1.3 kW, the main reactions contributing to gas heating were the ion–ion recombination reactions and the quenching of electronically excited chlorine atoms. At the same time, the energy released by the electron impact dissociation reaction of molecular chlorine is negligible due to its high degree of dissociation within the plasma bulk. The comparison between the results of our simulations and the fitting equation proposed in the literature show qualitative agreement, although there is significant quantitative discrepancy.

Hermes-3: Multi-component plasma simulations with BOUT++

A new open source tool for fluid simulation of multi-component plasmas is presented, based on a flexible software design that is applicable to scientific simulations in a wide range of fields. Hermes-3 is built on plasma simulation framework BOUT++, consolidating earlier SD1D and Hermes models into a single code that can be configured at run-time to solve plasma models in 1D, 2D or 3D, either for transport (steady-state) or turbulent (time-evolving) problems, with an arbitrary number of ion and neutral species. We describe the improved numerical algorithms and software design that have been implemented in Hermes-3.To demonstrate the capabilities of this tool, applications relevant to the boundary of tokamak plasmas are presented: 1D simulations of diveror plasmas evolving equations for all charge states of neon and deuterium; 2D transport simulations of tokamak equilibria in single-null X-point geometry with plasma ion and neutral atom species; and simulations of the time-dependent propagation of plasma filaments (blobs).Hermes-3 is publicly available on Github under the GPL-3 open source license. The repository includes documentation and a suite of unit, integrated and convergence tests. Program summaryProgram Title: Hermes-3CPC Library link to program files:https://doi.org/10.17632/rtkvc3r743.1Developer's repository link:https://github.com/bendudson/hermes-3Licensing provisions: GPLv3Programming language: C++14Nature of problem: Simulation of dynamics and steady-state solutions of multiple ion and neutral species in magnetised plasmas using drift-reduced fluid plasma models in 1 to 3 spatial dimensions. Focused on but not restricted to the simulation of turbulence in the boundary of tokamaks.Solution method: Hermes-3 implements a system of flexible software components, that are configured at run-time and coupled by passing a state object between them (the Encapsulate State and Command patterns). Hermes-3 builds on BOUT++ [1], employing the Method of Lines with implicit and explicit time-integration methods in curvilinear coordinates on block-structured meshes, making use of libraries including SUNDIALS, PETSc and Hypre. Hermes-3 implements drift-reduced plasma fluid equations using conservative finite difference methods, and atomic processes that couple plasma and neutral fluids.

Combining the hybrid mimetic mixed method with the Scharfetter-Gummel scheme for magnetised transport in plasmas

In this paper, we propose a numerical scheme for fluid models of magnetised plasmas. One important feature of the numerical scheme is that it should be able to handle the anisotropy induced by the magnetic field. In order to do so, we propose the use of the hybrid mimetic mixed (HMM) scheme for diffusion. This is combined with a hybridised variant of the Scharfetter-Gummel (SG) scheme for advection. The proposed hybrid scheme can be implemented very efficiently via static condensation. Numerical tests are then performed to show the applicability of the combined HMM-SG scheme, even for highly anisotropic magnetic fields.

Numerical simulation of magnetization induction coupled discharge plasma discharge process

Abstract First, this paper analyzes the plasma discharge and fluid model, and constructs the plasma discharge model by drift-diffusion approximation control equation, heavy particle component control equation, electric field distribution and volume force calculation, and plasma chemical kinetic model. Next, the coupling mechanism of inductively coupled RF plasma and its discharge characteristics are analyzed. Finally, the magnetized inductively coupled plasma discharge is simulated numerically. The results demonstrate that the current flowing on an inductor coil develops quicker at 0.045T and then calms down with an increase in the supplied constant dynamic magnetic field power, but the coil voltage exhibits the reverse effect.

Design and Study of a Large-Scale Microwave Plasma Torch with Four Ports

The microwave plasma torch (MPT) has gained popularity in industrial applications due to its high energy density, ionization levels, and high temperature. However, the non-uniform and unstable plasma generated by microwave plasma sources has limited the production of large-scale MPTs. This paper proposes a novel MPT device utilizing a four-port microwave source (2.45 GHz, 4 kW) to address these issues. The improved plasma uniformity and stability are achieved through the new structure, and the microwave efficiency is enhanced by introducing the focusing dielectric in the coupled cavity. Using a 3D electromagnetic field model, microwave plasma model and fluid model, the paper optimizes the geometry and inlet mode of the MPT device. Experimental results show that the novel MPT device can generate a plasma torch with a maximum height of 545 mm, a working range of 10–95 L/min, and a microwave efficiency up to 86%. The proposed device not only competently meets industrial requirements, but also provides design ideas and methods for future MPT devices.

Generalized Hamiltonian drift-fluid and gyrofluid reductions

We provide a procedure for deriving Hamiltonian reduced fluid models for plasmas, starting from a Hamiltonian gyrokinetic system in the approximation. The procedure generalizes, to a considerable extent, previous results. In particular, the evolution of moments with respect to the magnetic moment coordinate is also taken into account, together with background density and magnetic inhomogeneities. In the limit of vanishing Finite Larmor Radius (FLR) effects, an infinite family of reduced electron drift-fluid equations is derived, evolving all the electron moments , with and , where N and M are arbitrary non-negative integers counting the maximum order of the moments taken with respect to the parallel velocity and to the magnetic moment coordinates, respectively. An analogous result is found in the gyrofluid case and applied to the ion species. The gyrofluid result holds for , finite FLR effects and a low ratio β e between electron internal pressure and magnetic guide field pressure. In both the drift and gyrofluid case, the key for the identification of the Hamiltonian structure resides in changes of variables based on orthogonal matrices that diagonalize the Jacobi matrices associated with Hermite and Laguerre polynomials. In terms of the transformed variables, the drift and gyrofluid equations are cast in a simple form, which reduces to advection equations for Lagrangian invariants in the two-dimensional case with homogeneous background. Because the procedure requires to evolve, for a particle species s, all and only the moments , with and , not every choice for the set of moments is admissible, for fixed N and M. This might also explain the scarcity of Hamiltonian reduced fluid models obtained so far which account for anisotropic temperature fluctuations.

Stabilized bi-cubic Hermite Bézier finite element method with application to gas-plasma interactions occurring during massive material injection in tokamaks

Development of a numerical tool based upon the high-order, high-resolution Galerkin finite element method (FEM) often encounters two challenges: First, the Galerkin FEMs give central approximations to the differential operators and their use in the simulation of the convection-dominated flows may lead to the dispersion errors yielding entirely wrong numerical solutions. Secondly, high-order, high-resolution numerical methods are known to produce high wave-number oscillations in the vicinity of shocks/discontinuities in the numerical solution adversely affecting the stability of the method. We present the stabilized finite element method for plasma fluid models to address the two challenges. The numerical stabilization is based on two strategies: Variational Multiscale (VMS) and the shock-capturing approach. The former strategy takes into account (the approximation of) the effect of the unresolved scales onto resolved scales to introduce upwinding in the Galerkin FEM. The latter adaptively adds the artificial viscosity only in the vicinity of shocks. These numerical stabilization strategies are applied to stabilize the bi-cubic Hermite Bézier FEM in the computational framework of the nonlinear magnetohydrodynamics (MHD) code JOREK. The application of the stabilized FEM to the challenging simulation of Shattered Pellet Injection (SPI) in JET-like plasma is presented. It is shown that the developed numerical stabilization model improves the stability of the underlying numerical algorithm and the computational cost required to reveal the complex physics is reduced. The physical and numerical models presented can be used to perform expensive simulations of the plasma applications in large computational domains such as JET, and ITER.

Cold plasma modes in the chiral Maxwell-Carroll-Field-Jackiw electrodynamics

In this work, we study the propagation and absorption of plasma waves in the chiral Maxwell-Carroll-Field-Jackiw (MCFJ) electrodynamics. The Maxwell equations are rewritten for a cold, uniform, and collisionless fluid plasma model, allowing us to determine the new refractive indices and propagating modes. The cases of propagation parallel and orthogonal to the magnetic field are examined considering a purely timelike CFJ background that plays the role of the magnetic conductivity chiral parameter. The collective electromagnetic modes are associated with four distinct refractive indices associated with right-circularly polarized and left-circularly polarized waves. For each index, the propagation and absorption zones are illustrated for some specific parameter values. In low-frequency regime, we have obtained modified helicons with right- and left-circularly polarizations. The optical behavior is investigated by means of the rotatory power (RP) and dichroism coefficient. The existence of a negative refraction zone enhances the rotatory power. It is also observed RP sign reversal, a feature of rotating plasmas.