Plasma Model Research Articles

Propagation of circularly polarized electromagnetic wave in magnetized spin plasma

The spin effects on the propagation characteristic of circularly polarized electromagnetic (EM) wave in high density strongly magnetized plasma are discussed based on the the classical hydrodynamical model of relativistic spin plasma. The dielectric coefficients for right-hand circularly polarized (RCP) and left-hand circularly polarized (LCP) waves are obtained. Results show that the spin effects can affect the propagation characteristic of circularly polarized EM wave dramatically. Provided the spin effect is strong enough, LCP waves can also propagate in the magnetized over-dense plasma, while RCP waves may not. The strength of spin effects can be enhanced by increasing the plasma density or/and EM wave intensity.

Behavior of plasma parameters under mirror and cusp magnetic fields in DC glow discharge – a numerical study

Abstract The plasma parameters in DC glow discharge can easily be controlled by externally applied magnetic field without disturbing internal / built-in process parameters. Though several works have been carried out to study the influence of the magnetic field at different configurations on the plasma parameters, a complete understanding on the behaviour of the plasma discharge under the mirror and cusp fields could not be achieved. Further studies on the same are needed to improve the efficiency of the DC glow discharge system for existing applications and to find the new applications. In this work, a 2D axis-symmetric non-equilibrium plasma model with drift-diffusion approach is developed to study the characteristics of the plasma in DC glow discharge through various plasma parameters under mirror and cusp fields. The effects of current and position of magnetic coils on electric potential, ionization rate, electron temperature and electron number density are predicted and discussed. The distribution of electron number density at various coil currents and positions under both mirror and cusp fields are presented and the operating conditions favorable for applications in surface modification are suggested.

Studies of electron-ion resonant recombination of li-like Si11+ions

Abstract Theoretical investigations of L-shell △n = 0, 1, and 2 (2s → 2p, 3 l ′ , and 4 l ′ ) dielectronic recombination (DR), as well as K-shell △n = 1, 2, and 3 (1s → 2 l ′ , 3 l ′ , and 4 l ′ ) DR of Li-like Si 11+ ions are performed, using the relativistic distorted-wave approach. Detailed resonance energies, resonance strengths, and recombination rate coefficients are calculated including complex configuration mixing, Breit interaction, and QED corrections. Good agreement is achieved between the theoretical rate coefficients and the experimental results. The plasma rate coefficients are also calculated for electron temperatures ranging from 10−1 to 106 eV, and an analytical formula is presented to fit the total plasma rate coefficients for convenient modeling of astrophysical and fusion plasmas.

Cross-calibration and first vertical ECE measurement of electron energy distribution in the TCV tokamak

Abstract This paper describes the first vertical electron cyclotron emission measurement of non-thermal electron distributions in the Tokamak à Configuration Variable. These measurements were conducted in runaway electron scenarios and in the presence of electron cyclotron current drive. Measured intensities of linearly polarized X- and O-mode radiation from fast electrons allow the analysis of the energy distribution. The measurements were made possible through the creation of an operational regime for the diagnostic that is free of thermal background radiation, in relaxed electron density operations. This operational regime notably enables the cross-calibration of the diagnostic system, relying on thermal plasma measurements and modeling with the ray-tracing code SPECE.

Kinetic model of a three-component plasma

In this work, a system of kinetic equations is constructed to study the processes in a three-component plasma. It is an analogue of the Krook model, which is widely used in the dy-namics of rarefied gases. The model is supposed to be used to study the processes in the channels of rocket electric thrusters.

Development and Diagnostic Study of the RF nitrogen Atom Source

With the development of thin film technology, particularly the increasing demand for research on group III nitride materials, the need for N ions has significantly increased. A wider and more efficient RF nitrogen atom source specifically designed for nitride film growth has been developed. To evaluate the nitrogen atom density from the RF nitrogen source, a global model of RF nitrogen gas discharge plasma was created. This model was calibrated by measuring the emission spectrum of RF nitrogen gas discharge, and the relationships between nitrogen atom density and actual working pressure, power, and discharge chamber size were determined. It was found that the nitrogen atom density increases with increasing power and initially increases with pressure before decreasing. Finally, the nitrogen atom density was determined, achieving an 18% nitrogen atom dissociation rate and a nitrogen atom density of 5×1012/cm³ at an RF power of 400 W. This density is significantly higher than that of nitrogen ions, indicating that this is a nitrogen atom source with high purity.

A Broadband X-Ray Investigation of Fast-spinning Intermediate Polar CTCV J2056–3014

Abstract We report on XMM-Newton, NuSTAR, and NICER X-ray observations of CTCV J2056–3014, a cataclysmic variable (CV) with one of the fastest-spinning white dwarfs (WDs) at P = 29.6 s. While previously classified as an intermediate polar, CJ2056 also exhibits the properties of WZ Sge–type CVs, such as dwarf novae and superoutbursts. With XMM-Newton and NICER, we detected the spin period up to ∼2 keV with 7σ significance. We constrained its derivative to | P ̇ | < 1.8 × 10 − 12 s s−1 after correcting for binary orbital motion. The pulse profile is characterized by a single broad peak with ∼25% modulation. NuSTAR detected a fourfold increase in unabsorbed X-ray flux coincident with an optical flare, in 2022 November. The XMM-Newton and NICER X-ray spectra at 0.310 keV are best characterized by an absorbed, optically thin three-temperature thermal plasma model (kT = 0.3, 1.0, and 4.9 keV), while the NuSTAR spectra at 3–30 keV are best fit by a single-temperature thermal plasma model (kT = 8.4 keV), both with Fe abundance Z Fe/Z ⊙ = 0.3. CJ2056 exhibits similarities to other fast-spinning CVs, such as low plasma temperatures and no significant X-ray absorption at low energies. As the WD’s magnetic field strength is unknown, we applied both nonmagnetic and magnetic CV spectral models (MKCFLOW and MCVSPEC) to determine the WD mass. The derived WD mass range (M = 0.7–1.0 M ⊙) is above the centrifugal breakup mass limit of 0.56 M ⊙ and consistent with the mean WD mass of local CVs (M ≈ 0.8–0.9 M ⊙).

Triple-oxygen isotopes of stony micrometeorites by secondary ion mass spectrometry (SIMS): Olivine, basaltic glass and iron oxide matrix effects for sensitive high-mass resolution ion microprobe-stable isotope (SHRIMP-SI).

Micrometeorites are extraterrestrial particles smaller than ~2mm in diameter, most of which melted during atmospheric entry and crystallised or quenched to form 'cosmic spherules'. Their parentage among meteorite groups can be inferred from triple-oxygen isotope compositions, for example, by secondary ion mass spectrometry (SIMS). This method uses sample efficiently, preserving spherules for other investigations. While SIMS precisions are improving steadily, application requires assumptions about instrumental mass fractionation, which is controlled by sample chemistry and mineralogy (matrix effects). We have developed a generic SIMS method using sensitive high-mass resolution ion micro probe-stable isotope (SHRIMP-SI) that can be applied to finely crystalline igneous textures as in cosmic spherules. We correct for oxygen isotope matrix effects using the bulk chemistry of samples obtained by laser ablation inductively coupled plasma mass spectrometry (LA-ICP-MS) and model bulk chemical compositions as three-component mixtures of olivine, basaltic glass and Fe-oxide (magnetite), finding a unique matrix correction for each target. Our first results for cosmic spherules from East Antarctica compare favourably with established micrometeorite groups defined by precise and accurate but consumptive bulk oxygen isotope methods. The Fe-oxide content of each spherule is the main control on magnitude of oxygen isotope ratio bias, with effects on δ18O up to ~6‰. Our main peak in compositions closely coincides with so-called 'Group 1' objects identified by consumptive methods. The magnitude of SIMS matrix effects we find is similar to the previous intraspherule variations, which are now the limiting factor in understanding their compositions. The matrix effect for each spherule should be assessed quantitatively and individually, especially addressing Fe-oxide content. We expect micrometeorite triple-oxygen isotope compositions obtained by SIMS to converge on the main clusters (Groups 1 to 4) after correction firstly for magnetite content and secondarily for other phases (e.g., basaltic glass) in each target.

Excited ion-scale turbulence by a magnetic island in fusion plasmas

The characteristics of ion-scale turbulence in the presence of a magnetic island are numerically investigated using a gyrokinetic model in fusion plasma. We observe that in the absence of the usual ion temperature gradient (ITG) drive gradient, a magnetic island and its flatten effect could drive ITG instability. The magnetic island (MI) not only drives high-n modes of ITG instability but also induces low-n modes of vortex flow. Moreover, as the magnetic island width increases, the width of the vortex flow also increases. This implies that wider islands may more easily induce vortex flows. The study further indicates that the saturated amplitude and transport level of MI-induced ITG turbulence vary with different magnetic island widths. In general, larger magnetic islands enhance both particle and heat transport. When the magnetic island width reaches to 21ρi, the turbulence-driven transport becomes the same level with the cases that ITG is driven by pressure gradients. Our findings indicate the presence of intricate nonlinear effects in the modulation of plasma turbulence by MIs. These effects are of significant importance for comprehending the phenomenon of nonlinear coupling in forthcoming tokamaks such as ITER.

Experimental characterization and 3D simulations of turbulent flames assisted by nanosecond plasma discharges

This paper presents quantitative experimental data generated for the validation of plasma-assisted combustion (PAC) simulations. These data are then used to validate the phenomenological model of Castela et al. They are also useful to test other PAC models. In the experiment presented here, Nanosecond Repetitively Pulsed (NRP) discharges are applied to a lean-premixed turbulent methane-air flame initially near the lean blow-off limit. The discharges significantly enhance the combustion and stabilize the flame after a few pulses. Electrical and optical diagnostics are employed to extensively quantify the transient and steady state of the plasma-stabilization process. The flame shape is characterized by OH* chemiluminescence imaging. In the discharge region, OH density profiles are obtained by 1D laser-induced fluorescence, and the gas temperature is measured by optical emission spectroscopy measurements. The local gas temperature increases by 1250 K, and the OH number density rises sevenfold when NRP discharges are applied. These results evidence the cumulative thermal and chemical effects of NRP discharges, which are especially challenging to replicate numerically. A Large Eddy Simulation (LES) of the experiment is performed. Combustion chemistry is modeled by an analytically reduced mechanism, while the plasma discharge is described by the low-CPU cost phenomenological model of Castela et al., which aims to capture the main thermal and chemical effects induced by the discharges. The model of Castela et al. is validated in the burnt gases by the remarkable agreement between the simulations and the experiments regarding the flame shape, the local gas temperature, and the OH number density. More generally, this work demonstrates the relevance of simplified plasma models in LES solvers to simulate complex plasma-assisted burners.

Numerical Experiment on the Influence of Granulation-induced Waves on Solar Chromosphere Heating and Plasma Outflows in a Magnetic Arcade

This paper offers a fresh perspective on solar chromosphere heating and plasma outflows, focusing on the contribution of waves generated by solar granulation. Utilizing a 2.5D numerical experiment for the partially ionized lower solar atmosphere, we investigate the dissipation of these waves and their impact on plasma outflows and chromospheric heating via ion-neutral collisions. Employing the JOint ANalytical and Numerical Approach code, we adopt two-fluid model equations, examining partially ionized hydrogen plasma dynamics, including protons+electrons and neutrals, treated as two separate fluids that are coupled through ion-neutral collisions. Our investigation focuses on a quiet solar chromosphere region characterized by gravitational stratification and magnetic confinement by an initially set single magnetic arcade. The primary source of the waves is the solar convection beneath the photosphere. Our results demonstrate that ion-neutral collisions result in the dissipation of such waves, releasing thermal energy that heats the chromosphere plasma. Notably, this is accompanied by upward-directed plasma flows. Finally, we conclude that wave dissipation due to ion-neutral collisions in the two-fluid plasma model induces chromosphere heating and plasma outflows.

Numerical and experimental study of TIG welding arc in high frequency longitudinal magnetic field

The phenomenon of arc compression in TIG welding with high-frequency longitudinal magnetic field (HF-LMF) has been observed. However, its mechanism remains unclear, and there is a lack of numerical investigation. This work presents a two-dimension axisymmetrical model of arc plasma in high-frequency longitudinal magnetic field based on the magnetohydrodynamics (MHD). The distributions of temperature, velocity and pressure of arc plasma in HF-LMF are investigated at different applied LMF frequencies ranging from 0 to 5000 Hz. In order to make the simulation results more convincing, the induced current and the electric force were considered. The results show that the HF-LMF could generate a circumferential electric field, causing particles rotating. However, since HF-LMF is variable, the direction of the electric field also changes, which causes the arc to rotate back and forth. When the frequency of HF-LMF reaches 2000 Hz, the negative pressure zone above the center of the anode disappeared and the arc was compressed. Thus, as the frequency of LMF increases, the distribution of current density, anodic heat flux, and arc pressure changes from a bimodal distribution to a unimodal distribution, and the peak value increases. The theoretical predictions were in good agreement with the experimental results.

Formation of Alfvénic supernonlinear waves in a plasma containing double spectral distributed electrons

The characteristics of nonlinear and supernonlinear Alfvén waves propagating in a multicomponent plasma composed of a double spectral electron distribution and positive and negative ions were investigated. The Sagdeev technique was employed, and an energy equation was derived. Our findings show that the proposed system reveals the existence of a double-layer solution, periodic, supersoliton, and superperiodic waves. The phase portrait and potential analysis related to these waves were investigated to study the main features of existing waves. It was also found that decreasing the electron temperature helps the superperiodic structure to be excited in our plasma model. Our results help interpret the nonlinear and supernonlinear features of the recorded Alfvén waves propagating in the ionosphere D-region.

Role of nonthermal solar wind protons in the excitation of electromagnetic proton-cyclotron instability: a kinetic theory based exact numerical investigation

Free transverse kinetic energy i.e. perpendicular temperature anisotropy of protons excite the electromagnetic ion/proton cyclotron instability which is pertained to waves associated with prevalent electromagnetic ion/proton cyclotron emissions in various natural regions of plasmas. The transverse dielectric response function of left hand circularly polarized electromagnetic proton cyclotron (EPC) instability is calculated for two models of nonthermal Cairns distributed plasmas. These models are distinguished according to the effective thermal velocities of protons. For the energetic nonthermal protons populations, nonthermality dependent effective temperature model is proposed which significantly contributes in the excitation of aforementioned plasma mode and cause an appreciable enhancement in the instability growth rate. Exact numerical solution of dispersion relation yields oscillatory real frequency and growth rate of instability. A comparative analysis is also carried out to examine the instability behavior in distinct nonthermal and thermal plasma models. Contemporary numerical investigations are highly beneficial to understand the intricate dynamics of space plasmas.

Impact of ionization and transport on pedestal density structure in DIII-D and Alcator C-Mod

Abstract This paper investigates the role of ionization on the pedestal structure using both measurements and modeling for H-mode plasmas on DIII-D and Alcator C-Mod to enhance our ability to predict pedestal behavior in future pilot plants. The impact of the neutral penetration depth on the pedestal density is investigated using dimensionally matching hydrogen and deuterium DIII-D H-mode discharges at low and high electron density. The DIII-D Lyman-α diagnostic measurements show that hydrogen neutrals penetrate deeper inside the plasma on both the high field and low field side, while the pedestal electron density structure is similar for both isotopes. However, as the opaqueness increases we observe that the pedestal density gradient becomes stiff, similar to prior observations on DIII-D and C-Mod (Mordijck 2020 Nuclear Fusion 60 082006). In addition, these results also confirm prior measured and modeled poloidal asymmetries in neutral densities, indicating that to make transport predictions, 2D neutral modeling is necessary. The first direct validation of SOLPS-ITER for the measured brightness, emissivity and neutral densities for three different confinement regimes on C-Mod is introduced. The SOLPS-ITER model shows good agreement, within the constrains of the model for all regimes. In addition, a comparison of SOLPS-ITER modeling for DIII-D and C-Mod shows that as opaqueness increases, the role of divertor fueling and thus poloidal asymmetries in the neutral density profiles decreases. Based on these experimental and modeling results we estimate the size of a potential particle pinch using typical values for the diffusion coefficient for both DIII-D and C-Mod H-mode discharges.

Simulations of the stationary Q = 10 and the exit phase from the flat-top of an ITER 15MA baseline scenario: predictive JINTRAC simulation with a consistent treatment of D and T in the whole plasma

Abstract Designing a robust termination scenario for a burning ITER plasma is a challenge that requires extensive core plasma and divertor modelling. The presented work consists of coupled core/edge/SOL/divertor simulations, performed with the JINTRAC code, to study the Q = 10 flat-top phase and exit phase of the ITER 15 MA/5.3 T DT scenario. The modelling utilizes the recently implemented option to treat deuterium and tritium separately in the SOL/divertor, enabling a consistent treatment of deuterium and tritium in the whole plasma volume, which is a unique capability of JINTRAC. In addition, these are the first JINTRAC simulations of this scenario that use a first-principles transport model to self-consistently model the ECRH power deposition and to include tungsten while keeping track of tungsten sputtering and accumulation. The flat-top simulations demonstrate the possibility of sustaining a steady state fusion Q of 10 using pure deuterium gas puffs together with DT mixed pellets, which is an option to make a more effective use of tritium. Simulations of the exit phase are set up sequentially, with each phase providing initial conditions for the next, starting with a density decay at full current and auxiliary power, and demonstrate the possibility of reducing the density robustly within a few seconds. Following the density decay, a subsequent auxiliary power ramp-down in H-mode is performed with a late H–L transition at low auxiliary power, which may provide an option for the optimization of the plasma termination. The final ramp-down phase consists of a current ramp-down in L-mode to 3.75 MA.

Target-site cefiderocol pharmacokinetics in soft tissues of healthy volunteers.

Cefiderocol may potentially be used to treat skin and soft tissue infections (SSTIs). However, the pharmacokinetics of cefiderocol in human soft tissues have not yet been determined. The objective of the present PK study was to investigate whether target-site concentrations of cefiderocol are sufficiently high for the treatment of SSTIs. In this pharmacokinetic study, a single intravenous dose of 2 g cefiderocol was administered to eight healthy male volunteers. Drug concentrations were determined in plasma, muscle and subcutis over 8 h. Free plasma concentrations were calculated using the plasma protein binding determined with ultrafiltration. Free tissue concentrations were obtained using microdialysis. Penetration ratios were calculated as AUC0-8h_free_tissue/AUC0-8h_free_plasma. A population pharmacokinetic model was developed, and the probability of target attainment (PTA) was determined using Monte Carlo simulations. Cefiderocol showed good tissue penetration, with mean penetration ratios ± standard deviation of 0.99 ± 0.33 and 0.92 ± 0.30 for subcutis and muscle, respectively. Cefiderocol pharmacokinetics in plasma were best described with a two-compartment model, and tissue concentrations were described by scaling the tissue concentrations to concentrations in the peripheral compartment of the plasma model. For a thrice-daily regimen with 2 g doses intravenously infused over 3 h, PTA was ≥90% for MIC values up to 4 mg/L, both based on free plasma and soft tissue pharmacokinetics. This study indicates that a dose of 2 g cefiderocol achieves concentrations in plasma considered sufficient for treating relevant bacterial species. Assuming a comparable PK/PD target for soft tissues, sufficiently high concentrations would also be achieved in soft tissues.

Apparent effective ionization coefficient in N2 and O2 gas mixtures determined with two separate methods

N2:O2 gas mixtures are important for applications of atmospheric pressure plasmas such as ozone production, air purification from VOCs and NOx and surface treatments. Fundamental parameters such as the effective ionization coefficient are inputs for theoretical plasma models for applications and must thus be accurately known. This work determined the apparent effective ionization coefficient in N2:O2 mixtures in a broad reduced electric field strength E/N range of 150–1200 Td with two separate methods and compared with BOLSIG+ calculations of reduced effective ionization coefficient. Additionally, the equilibrium distance required to establish a steady-state electron energy distribution was estimated from spatial profiles of optical emission.

A model for fast electron-driven high-density plasma in the double-cone ignition scheme

Abstract A model for fast electron-driven high-density plasma is proposed to describe the effect of injected fast electrons on the temperature and inner pressure of the plasma in the fast heating process of the double-cone ignition (DCI) scheme. Due to the collision of the two low-density plasmas, the density and volume of the high-density plasma vary. Therefore, the ignition temperature and energy requirement of the high-density plasma vary at different moments, and the required energy for hot electrons to heat the plasma also changes. In practical experiments, the energy input of hot electrons needs to be considered. To reduce the energy input of hot electrons, the optimal moment and the shortest time for injecting hot electrons with minimum energy are analyzed. In this paper, it is proposed to inject hot electrons for a short time to heat the high-density plasma to a relatively high temperature. Then, the alpha particles with the high heating rate and PdV work heat the plasma to the ignition temperature, further reducing the energy required to inject hot electrons. The study of the injection time of fast electrons can reduce the energy requirement of fast electrons for the high-density plasma and increase the probability of successful ignition of the high-density plasma.

Full-scale modeling and experimental study of a gas neutron tube with a Penning ion source

This paper studies full-scale gas neutron tube modeling, including the following processes: gas discharge combustion in a Penning ion source, particle motion in an ion optical system, and modeling the in-target processes, such as sputtering, diffusion, thermal desorption of hydrogen isotopes, and nuclear reactions. Plasma modeling in quadrupole electric and axial magnetic fields was based on the electrostatic particle-in-cell method with molecular kinetic processes. The TechX Vsim software package was used. The neutron tube element sputtering by ions was simulated using SRIM/TRIM software based on Monte–Carlo methods. The OpenFOAM, as an open integrated platform for numerical simulation in continuum mechanics, was used to calculate the hydrogen thermal desorption activated by ion irradiation. The time-dependent neutron yield modeling was performed using the Geant4 software based on Monte–Carlo methods with CHIPS-TPT VNIIA-developed library. In addition, an experimental study of a gas neutron tube with a Penning ion source was conducted here as well. Details are given on the experiment and measurement technique used in this study. The operating characteristics for the gas neutron tube, including amplitude-time characteristics of current flashes (discharge and extraction currents), were determined. The neutron flux dependencies on the discharge current at various accelerating voltages were also obtained. Finally, a comparison between the experimental and calculated results is presented.