Processing math: 0%
Advanced Search+
Jiangbo ZHANG, Hongxu GAO, Fei XIAO, Wei LIU, Taixin LIANG, Zhongliang MA, Jian WU. Effect of shock wave formation on propellant ignition in capillary discharge[J]. Plasma Science and Technology, 2022, 24(6): 065504. DOI: 10.1088/2058-6272/ac5e71
Citation: Jiangbo ZHANG, Hongxu GAO, Fei XIAO, Wei LIU, Taixin LIANG, Zhongliang MA, Jian WU. Effect of shock wave formation on propellant ignition in capillary discharge[J]. Plasma Science and Technology, 2022, 24(6): 065504. DOI: 10.1088/2058-6272/ac5e71

Effect of shock wave formation on propellant ignition in capillary discharge

More Information
  • Corresponding author:

    Hongxu GAO, E-mail: gordon888@163.com

    Zhongliang MA, E-mail: mazhongliang@nuc.edu.cn

  • Received Date: January 04, 2022
  • Revised Date: February 12, 2022
  • Accepted Date: March 15, 2022
  • Available Online: December 12, 2023
  • Published Date: May 18, 2022
  • In order to study the effect of shock wave formation on propellant ignition in capillary discharge, the shock wave formation process was analyzed using experimental and theoretical methods; the plasma jet temperature was measured, and closed bomb and 30 mm gun experiments were carried out. The results show that the first shock wave has a smaller value and larger range of influence, while the second shock wave has a larger value and smaller range of influence. A plasma jet can generate a shock wave at the nozzle according to the calculated plasma pressure and velocity, which is well confirmed by experiments and calculations. The plasma jet temperature is high during the formation of a shock wave and then decreases sharply. Plasma ignition can increase the burning rate of a propellant by about 30% by increasing the burning surface area of the propellant. Compared to conventional ignition, the average maximum chamber pressure and average muzzle velocity of plasma ignition are increased by 9.1 MPa and 29.3 m·s-1 (~3%), respectively, in a 30 mm gun. Plasma ignition has strong ignition ability and short ignition delay time due to the generation of a shock wave. By increasing the burning rate of the propellant, the muzzle velocity can be greatly improved when the maximum chamber pressure increases a little. The characteristics of the shock wave can be applied in the application of the capillary discharge plasma. For example, it can be applied in fusion, launching and combustion.

  • Because ion-acoustic (IA) waves propagate through positively charged ions, they frequently interact with their electromagnetic fields. Basically, these are longitudinal perturbations of plasma density leading to the propagation of electrons and dynamical ions in a phase space [1]. The coupling of ions and electrons is provided by an electric field. The theoretical prediction of IA waves was carried out in 1929 by Tonks and Langmuir [2] whereas their experimental verification was done in 1933 by Rewans [3]. Subsequently, other scientists [46] performed many experiments and confirmed the existence of these waves. Their discovery has motivated the researchers towards investigation of their propagation and excitation in different plasma environments, e.g., see references [712]. Ion-acoustic waves exhibit many applications in various regimes of the plasmas. Many authors [1315] showed the significant damping of IA waves in collisionless plasmas for same temperatures of both ions and electrons. Later, Bhadra and Varma [16] examined the properties of IA waves by including the collisional effects of ion-ion interactions. In further studies, Buti [17] also incorporated the effects of Coulomb collisions on the propagation of IA waves and found that such collisions increase damping. The prediction of the existence of super nonlinear waves was initially made in a two-component electron-ion plasma featuring Maxwellian-distributed electrons [18]. Recently, the nonlinear characteristics of ion acoustic waves that arise from the interaction between the ionospheric plasma of Mars and the solar wind are investigated [19].

    Ocean circulation, quantum vortices in liquid nitrogen, spiral galaxies in the Milky Way, solar winds, typhoon vortices, optical and electromagnetic fields, and other examples of vortex-like formations can all be found in nature. The vortex beams are characterized by the Laguerre-Gaussian (LG) modes through Hilbert factor exp(i), where “l” represents the azimuthal mode index and carries finite amount of orbital angular momentum (OAM) [20]. The vortex beam is a classical tool for investigating the characteristics of optical vortices because it can be easily generated in the laboratory [21]. The presence of angular momentum and co-rotating solar wind associated with the plasma species has been found in the study [22] of comet Halley. It is interesting to note that when rotating solar wind interacts with a cometary plasma, it induces the generation of angular momentum in the plasma particles. This interaction gives rise to helical field structures, which are commonly observed in various astrophysical plasmas such as interplanetary clouds, coronal plasmas, solar jets, and other similar phenomena.

    The exchange of energy and momentum occurs when electromagnetic waves (laser) interact with the matter [23]. It is now confirmed that both spin and orbital parts contribute in the total angular momentum of electromagnetic radiations. The OAM is an extrinsic characteristic of electromagnetic (EM) radiations, exhibiting the helical or non-uniform wavefronts. In contrast, spin angular momentum (SAM) is associated with planar wavefronts and is considered an inherent property of EM radiations. Beth [24] was the first to perform an experiment of transferring angular momentum to a half wave plate. Allen et al [25] theoretically proposed the vortex beams with finite OAM states and identified the propagation of optical vortices in paraxial beams. Yao and Padgett [26] investigated the mechanism of generating non-planar structures and corresponding OAM states. Mendonça [27] identified the excitation of twisted modes due to finite OAM states. Later, various environments of space and astrophysical plasmas are identified where distinct and finite OAM states might contribute to plasma excitation. The examples of these environments include the co-rotating solar winds and spinning regions of the Jovian planets. The presence of helical phase structures in space plasmas results in notable effects of OAM on the propagation of helicon waves [28]. Recent studies indicated an increasing interest on the inclusion of OAM states in the study of plasma modes [2933]. These studies include the investigation of the characteristics of the nonlinear twisted waves [30], plasmons with twisted effects [31], magnetic field generation by LG light beams [32], phonons with finite OAM states [33], excitation of twisted waves [34], twisted DA waves with self-gravitational effects [35], instability of ion-acoustic twisted waves in permeating and superthermal plasmas [36, 37].

    In present work, a kinetic theory of the plasma physics is employed for the study of unstable ion-acoustic twisted waves in a plasma with (a) single and (b) double species of electrons. A linearized set of Vlasov and Poisson equations is used for the investigation of twisted waves. Analytical expressions for plasma response function are derived for the examination of twisted IA waves. The instability conditions for two distinct cases are also obtained. We particularly concentrate on the influence of finite OAM states on the instability conditions and growth rate of waves in a plasma which is in thermal equilibrium. Such thermal equilibrium may be attained after many successive collisions of the plasma species. The impact of an extra specie of electrons on the growth rate is highlighted.

    We assume an unmagnetized and collisionless plasma composed of drifting electrons and non-drifting positive ions. Any type of perturbation in a plasma from equilibrium corresponds to perturb the distribution function and plasma densities, fα = fα0 + fα1 and nα = nα0 + nα1 with α = e, i refers to electrons and ions. The digits 0 and 1 in the subscripts represent the quantities in equilibrium (unperturbed) and non-equilibrium (perturbed) states, respectively. The physical model employed in the present study is based on the collisionless Vlasov-Poisson equations in which the particle distribution function fα(r, v, t) evolves according to the Vlasov equation,

    \left(\frac{\partial}{\partial t}+V_{\alpha}\cdot\nabla\right)f_{\alpha1}\left(\boldsymbol{r},\boldsymbol{\ v},\ t\right)=\frac{q_{\alpha}}{m_{\alpha}}\boldsymbol{E}_1\cdot\nabla_{v\alpha}f_{\alpha0}\left(\boldsymbol{v}\right), (1)

    where qα (mα) is the charge (mass) of α species, respectively. After many successive interparticle collisions, a thermal equilibrium is attained in a plasma and particle distribution function in such a situation will always exhibit the Maxwellian velocity distribution as

    f_{\alpha0}\left(\boldsymbol{v}\right)=n_{\alpha}\left(\frac{m_{\alpha}}{2\text{π}k_{\mathrm{B}}T_{\alpha}}\right)^{\frac{1}{2}}\mathrm{exp}\left[-\frac{m_{\alpha}\left\{v_r^2+v_{\varphi}^2+\left(v_z-v_{\alpha0}\right)^2\right\}}{2k_{\mathrm{B}}T_{\alpha}}\right], (2)

    where nα, vα and Tα are the density, velocity and temperature of the α-species, respectively. The Vlasov equation can be more conveniently solved by expressing the electric field in the form of electrostatic potential as, {E}_{1}=-\nabla {\phi }_{1} and is examined by the Poisson’s equation

    \nabla^2\phi_1\left(\boldsymbol{r},\ t\right)=\mathrm{ }-\sum_{\alpha}^{ }4\text{π}q_{\alpha}\int_{ }^{ }f_{\alpha1}\left(\boldsymbol{r},\boldsymbol{\ v},\ t\right)\mathrm{d}\boldsymbol{v}_{\alpha}, (3)

    with fα1(r, v, t)\, \ll \,[fα0(v)] being the perturbed distribution function (DF) [Maxwellian equilibrium DF] and \phi_1\left(r,\ t\right) refers to the perturbed shielding potential. Throughout the present analysis, a field free plasma is assumed accounting for zero external electric (E0) and magnetic (B0) fields and the plasma equilibrium state is characterized by charge neutrality condition \sum n_\alpha q_\alpha =0.

    To analyze the effect of finite OAM states [27] on the electrostatic IA waves, the electric field can be expressed as E_1=-\mathrm{i}k_{\mathrm{e}\mathrm{f}\mathrm{f}}\phi_1\left(r,\ t\right) in cylindrical coordinates system due to helical trajectories via new effective wave vector {k}_{\mathrm{e}\mathrm{f}\mathrm{f}}= \dfrac{-\mathrm{i}}{{\psi }_{pl}}{\partial }_{r}{\psi }_{pl}\hat{r}+\dfrac{1}{r}\hat{\varphi }+\left(k-\dfrac{\mathrm{i}}{{\psi }_{pl}}{\partial }_{ {\textit{z}} }{\psi }_{pl}\right)\hat{ {\textit{z}} } with \hat{r}, \hat{\varphi } , \hat{z} being the unit vectors. The perturbed electrostatic potential {\phi }_{1} is now proportional to \mathrm{exp}\left[\mathrm{i}\left(l\varphi +k {\textit{z}} -\omega t\right)\right] with an extra effect of azimuthal angle “\varphi ” due to phase structure which is always absent in plane wave solutions.

    The perturbed distribution function [27] can be written in terms of generalized LG beam solution in this form

    f_{\alpha1}\left(\boldsymbol{r},\boldsymbol{\ v},\ t\right)=f_{pl}\left(\boldsymbol{v}\right)\psi_{pl}\left(r,\ z\right)\mathrm{exp}\left[\mathrm{i}\left(l\varphi+kz-\omega t\right)\right], (4)

    where the LG mode functions and normalization constant are defined as {\psi }_{pl}\left(r,z\right)=C_{pl}L_p^{|l |} (\zeta )\zeta^{|l |} \;{\mathrm{exp}}[-\zeta/2] and Cpl = (1/2)π−1/2[(l + p)!/p!]1/2, respectively. Moreover, fpl(v) is mode amplitude whereas L_p^{|l|}(\zeta)=\left(\dfrac{\zeta^{|-l|}{\mathrm{exp}}(\zeta)} {p!} \right)\dfrac{{\mathrm{d}}p}{{\mathrm{d}}\zeta ^p}[\zeta^{p+|l|}\;{\mathrm{exp}}(\zeta)] indicates the Laguerre polynomials in terms of Rodrigues formula with \zeta = [r/w(z)]2 and w(z) the beam waist. Likewise, the twisted electrostatic potential {\phi }_{1} can be expressed by taking into account the helical phase structures, as

    \phi_1\left(\boldsymbol{r},\ t\right)=\phi_{pl}\psi_{pl}\left(r,\text{ z}\right)\mathrm{exp}\left[\mathrm{i}\left(l\varphi+kz-\omega t\right)\right], (5)

    here, p (= 0, 1, 2, ...) is the radial mode index and l (= 0, ±1, ±2, ...) represents the azimuthal mode index or quantum number which is proportional to the OAM. It can be seen that equations (4) and (5) may reduce to the Gaussian plane beam solution by assuming p = l = 0, which is the lowest order mode of the LG beams.

    To proceed further, we adopt the standard procedure [38, 39] for obtaining the generalized dielectric response function with contribution of finite OAM in an electrostatic collisionless electron-ion plasma, as

    \begin{split} D\left(k,\ lq_{\varphi},\ \omega\right)= & 1+\frac{2\omega_{\mathrm{p}\mathrm{e}}^2}{k^2v_{\mathrm{t}\mathrm{e}}^2}\left[2+\mu_{z\mathrm{e}}Z\left(\mu_{z\mathrm{e}}\right)+\mu_{\varphi\mathrm{e}}Z\left(\mu_{\varphi\mathrm{e}}\right)\right] \\ & +\frac{2\omega_{\mathrm{p}\mathrm{i}}^2}{k^2v_{\mathrm{t}\mathrm{i}}^2}\left[2+\mu_{z\mathrm{i}}Z\left(\mu_{z\mathrm{i}}\right)+\mu_{\varphi\mathrm{i}}Z\left(\mu_{\varphi\mathrm{i}}\right)\right], \end{split} (6)

    where, the axial and azimuthal dispersion functions are respectively, defined as Z\left({\mu }_{z\alpha }\right)={{\text{π}}}^{-\frac{1}{2}}{\Large \int }_{-\mathrm{\infty }}^{+\mathrm{\infty }}\dfrac{\mathrm{d}q{\mathrm{e}}^{-{q}^{2}}}{q-{\mu }_{z\alpha }} and Z\left({\mu }_{\varphi \alpha }\right)={{\text{π}}}^{-\frac{1}{2}}{\Large \int }_{-\mathrm{\infty }}^{+\mathrm{\infty }}\dfrac{\mathrm{d}q{\mathrm{e}}^{-{q}^{2}}}{q-{\mu }_{\varphi \alpha }} with arguments {\mu }_{z\mathrm{e}} ≡ ( \omega kve0)/kvte taking into account drifting electrons and {\mu }_{z\mathrm{i}}=\omega /k{v}_{\mathrm{t}\mathrm{i}} for non-streaming positive ions while {\mu }_{\varphi \alpha } =\omega /l{q}_{\varphi }{v}_{\mathrm{t}\alpha } represents the azimuthal argument. The third term in both square brackets comes from the contribution of OAM states associated with nonplanar wavefronts and by ignoring it, we can obtain the expression for planar waves. For twisted IA waves, the resonance condition is modified as \omega =k{v}_{z}-l{q}_{\varphi }{v}_{\alpha }. It is interesting that, in the limit l→0, the dielectric constant (6) is reduced to the planar case which shows its generality.

    The main aim of the present study is to investigate the twisted IA modes and the effects of finite OAM states on instability growth rate. For this, we employ constraints on the plasma waves of the form vti \ll vph \left(= \dfrac{\omega }{k}\right) \ll vte, with vte(vti) is the thermal speed of electrons(ions). In present case, the axial and azimuthal arguments for plasma species obey the following conditions; {\mu }_{z\mathrm{e}}\ll 1 , {\mu }_{z\mathrm{i}}\gg 1 and {\mu }_{\varphi \alpha }\gg1. For further proceedings of the mathematical analysis, asymptotic expansions are used for electrons and ions. As a result, the dielectric constant D\left(k,\ lq_{\varphi},\ \omega\right) takes the form

    D\left(k,\ lq_{\varphi},\ \omega\right)=1+\sum_{\alpha=\mathrm{e},\mathrm{\ i}}^{ }\chi_{\alpha}\left(k,\ lq_{\varphi},\ \omega\right), (7)

    where sucesseptibilities of the electrons \chi_{\mathrm{e}}\left(k,\ lq_{\varphi},\ \omega\right) and ions \chi_{\mathrm{i}}\left(k,\ lq_{\varphi},\ \omega\right) are respectively defined as

    \chi_{\mathrm{e}}\left(k,\ lq_{\varphi},\ \omega\right)=\frac{1}{k^2\lambda_{\mathrm{D}\mathrm{e}}^2}+\frac{\mathrm{i}\sqrt{\text{π}}}{k^2\lambda_{\mathrm{D}\mathrm{e}}^2}\left\{\frac{\omega-kv_{\mathrm{e}0}}{v_{\mathrm{t}\mathrm{e}}k}+\frac{\beta\omega}{v_{\mathrm{t}\mathrm{e}}k}\mathrm{exp}\left(-\frac{\beta^2\omega^2}{v_{\mathrm{t}\mathrm{e}}^2k^2}\right)\right\}, (8)

    and

    \begin{split} {\chi }_{\mathrm{i}}\left(k, l{q}_{\varphi }, \omega \right)=&\frac{\mathrm{i}\sqrt{{\text{π} }}}{{k}^{2}{\lambda }_{\mathrm{D}\mathrm{i}}^{2}}\left\{\frac{\omega }{{v}_{\mathrm{t}\mathrm{i}}k}\mathrm{exp}\left(-\frac{{\omega }^{2}}{{v}_{\mathrm{t}\mathrm{i}}^{2}{k}^{2}}\right)+\frac{\beta \omega }{{v}_{\mathrm{t}\mathrm{i}}k}\mathrm{exp}\left(-\frac{{\beta }^{2}{\omega }^{2}}{{v}_{\mathrm{t}\mathrm{i}}^{2}{k}^{2}}\right)\right\} \\&- \frac{1}{{k}^{2}{\lambda }_{\mathrm{D}\mathrm{i}}^{2}}\left\{\frac{{k}^{2}{v}_{\mathrm{t}\mathrm{i}}^{2}}{2{\omega }^{2}}\left(1+\frac{3{k}^{2}{v}_{\mathrm{t}\mathrm{i}}^{2}}{{\omega }^{2}}\right)\right.\\&\left.+\frac{{k}^{2}{v}_{\mathrm{t}\mathrm{i}}^{2}}{2{\omega }^{2}{\beta }^{2}}\left(1+\frac{3{k}^{2}{v}_{\mathrm{t}\mathrm{i}}^{2}}{{\omega }^{2}{\beta }^{2}}\right)\right\}. \end{split} (9)

    Here, \lambda_{\mathrm{D}\mathrm{e}}\left[=\left(\dfrac{\varepsilon_0T_{\mathrm{e}}}{n_{\mathrm{e}0}e^2}\right)^{\Large\frac{1}{2}}\right] and \lambda_{\mathrm{D}\mathrm{i}}\left[=\left(\dfrac{\varepsilon_0T_{\mathrm{i}}}{n_{\mathrm{i}0}e^2}\right)^{\Large\frac{1}{2}}\right] are the Debye lengths of electrons and ions, respectively. The twist parameter β\left (=\dfrac{k}{l{q}_{\varphi }}\right ) is the ratio of longitudinal to azimuthal wave number. The real part, i.e. Re[ D\left(k,\ lq_{\varphi},\ \omega\right) ] = 0 of the dielectric constant leads to the wave frequency, as

    \frac{\omega_{\mathrm{i}\mathrm{a}}}{k}=\frac{C_{\mathrm{i}\mathrm{a}}}{\sqrt{1+k^2\lambda_{\mathrm{D}\mathrm{e}}^2}}\left\{\frac{1+\beta^2}{\beta^2}+\frac{3\beta^{-2}}{\text{Δ}}\left(\frac{1+\beta^4}{1+\beta^2}\right)\right\}^{\Large\frac{1}{2}}, (10)

    where, {C}_{\mathrm{i}\mathrm{a}} = {{\omega }_{\mathrm{p}\mathrm{i}}\lambda }_{\mathrm{D}\mathrm{e}} \sqrt{\dfrac{{T}_{\mathrm{e}}}{{m}_{\mathrm{i}}}} is the IA speed and {\text{Δ}} =\dfrac{{T}_{\mathrm{e}}}{{T}_{\mathrm{i}}} is the electron-to-ion temperature ratio. The dispersion relation (10) is significantly modified with twist parameter β. The approximation β→∞, leads equation (10) in its plane wave solution.

    In plasmas, the constituent particles may have their speeds either slower or faster than the wave speed. In case of larger number of faster particles, the resonant particles loose net amount of energy by transferring it to the wave which exhibits growth rate of instability. Otherwise, resonant particles gain energy from the wave resulting in the damping of wave. The instability growth rate of IA waves carrying OAM states is obtained, as

    \frac{\gamma_{\mathrm{i}\mathrm{a}}}{\omega_{\mathrm{i}\mathrm{a}}}=\mathrm{ }\mathrm{ }\mathrm{ }\mathrm{ }-\sqrt{\frac{\rm{\pi}}{8}}\sqrt{\frac{1+\beta^2}{\beta^2\left(1+k^2\lambda_{\mathrm{D}\mathrm{e}}^2\right)^3}}\left[\sqrt{\frac{m_{\mathrm{e}}}{m_{\mathrm{i}}}}\left(1-\frac{v_{\mathrm{e}0}}{v_{\mathrm{p}\mathrm{h}}}+\Gamma \right)\right], (11)

    where Γ is defined as

    \Gamma=\beta\sqrt{\frac{m_{\mathrm{e}}}{m_{\mathrm{i}}}}\mathrm{exp}\left(-\frac{\beta^2\omega_{\mathrm{i}\mathrm{a}}^2}{v_{\mathrm{t}\mathrm{e}}^2k^2}\right)+ {\text{Δ}} ^{\Large\frac{3}{2}}\left\{\mathrm{exp}\left(-\frac{\omega_{\mathrm{i}\mathrm{a}}^2}{v_{\mathrm{t}\mathrm{i}}^2k^2}\right) +\beta\mathrm{exp}\left(-\frac{\beta^2\omega_{\mathrm{i}\mathrm{a}}^2}{v_{\mathrm{t}\mathrm{i}}^2k^2}\right)\right\}. (12)

    Equation (11) shows that the growth rate is modified significantly by the streaming speed ve0 and parameter β. The characteristic behavior of the wave, i.e. either it is growing or is damped depends on the sign of (1 − ve0/vph) in equation (11). Hence, the cutoff value for ve0, above which instability arises is calculated as

    \begin{split} {v}_{\mathrm{e}0} > & \left[\beta \mathrm{exp}\left(-\frac{{\beta }^{2}{\omega }_{\mathrm{i}\mathrm{a}}^{2}}{{v}_{\mathrm{t}\mathrm{e}}^{2}{k}^{2}}\right)+\sqrt{{ {\text{Δ}} }^{3}\frac{{m}_{\mathrm{i}}}{{m}_{\mathrm{e}}}}\right.\\& \left. \times \left\{\mathrm{exp}\left(-\frac{{\omega }_{\mathrm{i}\mathrm{a}}^{2}}{{v}_{\mathrm{t}\mathrm{i}}^{2}{k}^{2}}\right)+\beta \mathrm{exp}\left(-\frac{{\beta }^{2}{\omega }_{\mathrm{i}\mathrm{a}}^{2}}{{v}_{\mathrm{t}\mathrm{i}}^{2}{k}^{2}}\right)\right\}+1 \right] . \end{split} (13)

    Equation (13) is the threshold limit for instability of IA waves with helical structures above which instability steps in invoking \gamma_{\mathrm{ia}} >0 .

    Here, we examine the twisted IA waves in a magnetic field free plasma, whose constituent particles are the ions, cold electrons and hot electrons. The phase speed of the wave is restricted as vti \ll ω/k < vtc, vth. The dispersion functions Z( {\mu }_{z\alpha } ) and Z({\mu }_{\varphi \alpha }) in equation (6) can be asymptotically expanded according to the constraints {\mu }_{z\mathrm{i}}\gg 1, ( {\mu }_{z\mathrm{c}} , {\mu }_{z\mathrm{h}} ) \ll 1 and {\mu }_{\varphi \alpha }\gg1 for ions, cold and hot electrons, respectively. Now, following the same procedure as in previous case, we eventually arrive at

    D\left(k,\ lq_{\varphi},\ \omega\right)=1+\sum_{\alpha=\mathrm{e},\mathrm{\ i}}^{ }\chi_{\alpha}+\chi_{\mathrm{i}}, (14)

    where

    \begin{split}{\chi }_{\alpha }=&\frac{1}{{k}^{2}{\lambda }_{\mathrm{D}\mathrm{e}\mathrm{f}}^{2}}+\mathrm{i}\sqrt{{\text{π}}}\frac{1}{{k}^{2}{\lambda }_{\mathrm{D}\mathrm{h}}^{2}}\left\{\frac{\omega -k{v}_{0}}{{v}_{\mathrm{t}\mathrm{h}}k}+\frac{\beta \omega }{k{v}_{\mathrm{t}\mathrm{h}}}\mathrm{exp}\left(-\frac{{\beta }^{2}{\omega }^{2}}{{v}_{\mathrm{t}\mathrm{h}}^{2}{k}^{2}}\right)\right\} \\&+\mathrm{i}\sqrt{\mathrm{\pi }}\frac{1}{{k}^{2}{\lambda }_{\mathrm{D}\mathrm{e}}^{2}}\left\{\frac{\omega }{{v}_{\mathrm{t}\mathrm{c}}k}+\frac{\beta \omega }{k{v}_{\mathrm{t}\mathrm{c}}}\mathrm{exp}\left(-\frac{{\beta }^{2}{\omega }^{2}}{{v}_{\mathrm{t}\mathrm{c}}^{2}{k}^{2}}\right)\right\}, \end{split} (15)

    and {\chi }_{\mathrm{i}} has the same value as given in equation (9). Here, {\lambda }_{\mathrm{D}\mathrm{e}\mathrm{f}}={\left[{\lambda }_{\mathrm{D}\mathrm{h}}^{-2}+{\lambda }_{\mathrm{D}\mathrm{c}}^{-2}\right]}^{-\frac{1}{2}} is the effective Debye length with \lambda_{\mathrm{D}\mathrm{h}}=\sqrt{\dfrac{\varepsilon_0T_{\mathrm{h}}}{n_{\mathrm{h}0}e^2}} and \mathrm{\lambda}_{\mathrm{D}\mathrm{c}}=\sqrt{\dfrac{\varepsilon_0T_{\mathrm{c}}}{n_{\mathrm{c}0}e^2}}. In present model, only hot electron streaming speed (vh0 = v0) is taken into account while drift speed of cold electrons is zero, i.e., vc0 = 0. Equation (14) represents the generalized form of dielectric constant of twisted IA waves in two-temperature electron plasma exhibiting finite OAM states. After the decomposition of equation (14) into real and imaginary parts and assuming D_{\mathrm{r}}\left(k,\ lq_{\varphi},\ \omega\right) = 0, we come up with the following dispersion relation

    \begin{split} \frac{{\omega }_{\mathrm{i}\mathrm{a}}}{k}=&{C'_{\mathrm{ia}}}{\Bigg\{ \frac{1}{1+{k}^{2}{\lambda }_{\mathrm{D}\mathrm{e}\mathrm{f}}^{2}}\frac{1+{\beta }^{2}}{{\beta }^{2}}+3{\beta }^{-2}} \\& \times \Bigg(\frac{f}{1+f}\frac {T_{\mathrm{i}}}{T_{\mathrm{h}}}+\frac{1}{1+f}\frac{T_{\mathrm{i}}}{T_{\mathrm{c}}}\Bigg)\Bigg(\frac{1+\beta^4}{1+\beta^2} \Bigg)\Bigg\}^{\Large \frac{1}{2}}.\end{split} (16)

    Note that the dispersive characteristics of twisted IA waves are now significantly affected by the modified IA speed in a two-temperature electron plasma. The modified IA speed is expressed as

    C'_{\mathrm{i}\mathrm{a}}=\omega_{\mathrm{p}\mathrm{i}}\lambda_{\mathrm{D}\mathrm{e}\mathrm{f}}\equiv\left(\frac{T_{\mathrm{c}}}{m_{\mathrm{i}}}\frac{1+f}{1+\frac{f}{\sigma}}\right)^{\Large\frac{1}{2}}, (17)

    where σ(=Th/Tc) and f(=nh0/nc0) are the hot-to-cold electron temperature and density ratios, respectively. The assumption nh0 = 0 simply implies that f = 0 while for Th = 0, one may obtain λDh = 0 and λDef = λDc. Thus, replacing the subscript c with e, equation (16) simply reduces to the case for IA waves in a plasma with single species of electron.

    The growth rate with two-temperature electrons retaining hot electron streaming effects can be written, as

    \begin{split}\frac{{\gamma }_{\mathrm{i}\mathrm{a}}}{{\omega }_{\mathrm{i}\mathrm{a}}}=&\sqrt{\frac{{\text{π}}}{8}}{\left(\frac{1+{\beta }^{2}}{{\beta }^{2}}\right)}^{\frac{1}{2}}{\left(\frac{{T}_{\mathrm{c}}}{{T}_{\mathrm{i}}}\right)}^{\frac{3}{2}}\frac{{\left(1+f\right)}^{\frac{3}{2}})}{{\left(1+\frac{f}{\sigma }\right)}^{\frac{3}{2}}{\left(1+{k}^{2}{\lambda }_{\mathrm{D}\mathrm{e}\mathrm{f}}^{2}\right)}^{\frac{3}{2}}}\\& \times \left[\frac{f}{1+f}\sqrt{\frac{{m}_{\mathrm{e}}}{{m}_{\mathrm{i}}}}{\left(\frac{{T}_{\mathrm{i}}}{{T}_{\mathrm{h}}}\right)}^{\frac{3}{2}}\left(1-\frac{{v}_{\mathrm{e}0}}{{v}_{\mathrm{p}\mathrm{h}}}+{\Gamma'}\right)\right], \end{split} (18)

    with

    \begin{split} \Gamma'=&\beta \frac{f}{1+f}\sqrt{\frac{{m}_{\mathrm{e}}}{{m}_{\mathrm{i}}}}{\left(\frac{{T}_{\mathrm{i}}}{{T}_{\mathrm{h}}}\right)}^{\frac{3}{2}}\mathrm{exp}\left(-\frac{{{\beta }^{2}\omega }_{\mathrm{i}\mathrm{a}}^{2}}{{v}_{\mathrm{t}\mathrm{h}}^{2}{k}^{2}}\right)+\frac{1}{1+f}\sqrt{\frac{{m}_{\mathrm{e}}}{{m}_{\mathrm{i}}}}{\left(\frac{{T}_{\mathrm{i}}}{{T}_{\mathrm{c}}}\right)}^{\frac{3}{2}} \\&\times \left\{1+\beta \mathrm{exp}\left(-\frac{{{\beta }^{2}\omega }_{\mathrm{i}\mathrm{a}}^{2}}{{v}_{\mathrm{t}\mathrm{c}}^{2}{k}^{2}}\right)\right\}+\Bigg\{\mathrm{exp}\left(-\frac{{\omega }_{\mathrm{i}\mathrm{a}}^{2}}{{v}_{\mathrm{t}\mathrm{i}}^{2}{k}^{2}}\right)\\&+\beta \mathrm{exp}\left(-\frac{{\beta }^{2}{\omega }_{\mathrm{i}\mathrm{a}}^{2}}{{v}_{\mathrm{t}\mathrm{i}}^{2}{k}^{2}}\right)\Bigg\} . \end{split} (19)

    The growth rate is now significantly modified with β, the ratio of hot-to-cold electrons density f and temperature ratio σ. In a limiting case, when nh0 = v0 = 0 and subscript c is simply replaced by e, the expression [40] for decay rate of twisted IA waves is recovered. Equation (18) satisfies the following instability condition

    {v}_{0} > {v}_{\mathrm{p}\mathrm{h}}\left[1+\frac{\Gamma'}{\large \frac{f}{1+f}\sqrt{ \frac{{m}_{\mathrm{e}}}{{m}_{\mathrm{i}}}}{\left(\frac{{T}_{\mathrm{i}}}{{T}_{\mathrm{h}}}\right)}^{\frac{3}{2}}}\right]. (20)

    which is the onset of the instability. The threshold limit for instability of twisted IA waves is modified with nh0 and Th. The ignorance of the contribution of OAM states assumes l → 0, which corresponds to β → ∞. In this situation, the factors \mathrm{exp}\left(-\dfrac{{{\beta }^{2}\omega }_{\mathrm{i}\mathrm{a}}^{2}}{{v}_{\mathrm{t}\mathrm{h}}^{2}{k}^{2}}\right) , \mathrm{exp}\left(-\dfrac{{{\beta }^{2}\omega }_{\mathrm{i}\mathrm{a}}^{2}}{{v}_{\mathrm{t}\mathrm{c}}^{2}{k}^{2}}\right) and \mathrm{exp}\left(-\dfrac{{{\beta }^{2}\omega }_{\mathrm{i}\mathrm{a}}^{2}}{{v}_{\mathrm{t}\mathrm{i}}^{2}{k}^{2}}\right) that appear as a multiple of twist parameter in equation (19) become zero more quickly as compared to β → ∞. This leads to vanish the products in equation (19) and at condition (l → 0), the analytical result is valid for planar IA waves.

    In this section, the numerical calculations in the kinetic framework are used to examine the results of equations (10), (11), (16) and (18). With non-zero streaming speed of electrons ve0 and finite OAM states, the wave frequency, its growth rate and cutoff limit for instability are drastically altered. Special attention and care are needed for the examination of the influence of these parameters on wave dynamics. It is interesting to mention that, the smaller values of azimuthal wave number l{q}_{\varphi }, i.e., l{q}_{\varphi } \ll k leads to the larger values of β, i.e. β \gg 1. In the limit of larger β values (i.e. β → ∞), the solution for straight propagating IA waves is retrieved. For the validation of obtained results, plasma parameters corresponding to the physical situations [41, 42] are used.

    The pictorial illustration of the wave frequency of IA waves exhibiting finite OAM via twist parameter [while ignoring spin angular momentum] is shown in figure 1, in which the impact of β on wave frequency is depicted. It is interesting to mention that parameter β (= k/l{q}_{\varphi }) describes the helical phase structures and is ultimately controlled by azimuthal number l (= 0, 1, 2, …).

    Figure  1.  Real frequency of twisted IA waves “ωia” versus normalized wave number “k” plots at different values of twist parameter “β”, i.e. β = 0.7 (black dashed curve), β = 1.0 (red solid curve), β = 1.5 (blue dotted curve) and β = 50 (green dashed curve) at ∆ = 100.

    The important point is that the handedness depends on the sign of l (i.e. either positive or negative) and is normally taken by convention. For larger l, i.e., l → ±∞, the twist parameter β → 0 and consequently at β = 0, a phase singularity arises. An increasing amplitude of the wave frequency is noticed in higher momentum states as smaller β values indicate the enhanced OAM effects. Therefore, it is easier to conclude that the real frequency of twisted IA waves is more as compared to the planar waves. Furthermore, for smaller wavelengths, i.e. at longer wave numbers, all curves admit large number of waves through a certain point in unit time.

    Figure 2 presents the contour plots of growth rate against wave number “k” and twist parameter “β” showing the impact of finite OAM. It has already been described in above discussion that for larger values of azimuthal wave number, β reduces which shows an enhanced twisted effect. Figure 2 shows the threshold at β = 10, above which IA wave is Landau damped while below this value it suffers the instability. Furthermore, the amplitude of the unstable wave increases dramatically at smaller values of β, i.e. β → 0. For a smaller value of β (higher value of azimuthal number “l”), a large number of resonant particles contribute their energy to the wave. It is also well-known that larger values of “β” correspond to the planar waves.

    Figure  2.  Contour plots of instability growth rate against wave number and azimuthal quantum number (l) via twist parameter β for ve0/Cia = 10.8 and ∆ = 95.

    Therefore, it can be expected that twisted waves exhibit more instability than planar waves. In our present problem, the source of free energy is the streaming speed of inertial electrons. In the absence of streaming speed, the transmission of energy from resonant species to the wave becomes almost zero and twisted ion-acoustic waves show damping for the whole range of β. The contour plots also highlighted another interesting fact that for shorter wavelengths (longer wave numbers), the damping rate vanishes and only growth rate of instability appears. This implies that the threshold limit, i.e. equation (13) for unstable waves is modified significantly with the twist parameter and wave number k.

    Figure 3 depicts the growth rate “γia” plotted against the twist parameter “β” and wave number k illustrating the changes in instability for three distinct values of the streaming speed “ve0”. The amplitude of the unstable waves first increases with wave number, reaches at the maximum value and then starts to decrease [all curves show almost same trend]. It can be observed that twisted IA waves gain more energy from the resonant plasma species for larger values of electron’s streaming speed.

    Figure  3.  Three dimensional plots of instability growth rate against twist parameter and normalized wave number at ve0 = 20.8Cia (upper curve with green grid lines), ve0 = 15.8Cia (middle curve with blue grid lines) and ve0 = 10.8Cia (bottom curve with red grid lines) at ∆ = 95.

    In plasmas, kinetic effects associated with the streaming particles, such as Landau damping and wave-particle resonances, become more significant at higher streaming speeds. These effects can modify the dispersion properties of the acoustic waves, potentially leading to increased instability. This observation confirms that the free energy source available for the excitation of twisted ion-acoustic waves is streaming speed. Figure 3 also reveals that the maxima of instability growth rate is shifted towards right with streaming speed (see 3D curve with green grid lines).

    Figure 4 indicates the growth rate of twisted IA waves against normalized wave number at distinct values of streaming speed for the case where k < l{q}_{\varphi }, so that β = k/l{q}_{\varphi } = 0.1. Here, we have noticed that twisted IA waves are Landau damped in longer wavelength limit while they suffer growth rate of instability for shorter wavelengths.

    Figure  4.  Growth rate of instability “γia” versus wave number “k” plots for distinct values of streaming speed ve0, i.e. ve0/Cia = 7.0 (black dashed curve), ve0/Cia = 7.5 (red dotted curve), ve0/Cia = 8.0 (green solid curve) and ve0/Cia = 8.5 (blue dashed curve) at fixed β = 0.1 and ∆ = 95.

    Thus, a cutoff exists at certain value of wave number above which waves become unstable instead of exhibiting damping. It can be seen that the threshold limit for unstable waves reduces by increasing the streaming speed. This shows that waves are more unstable when larger number of electrons with high streaming speeds are available. It is important to note that the increase in streaming speed of particles can influence the instability of acoustic waves through various mechanisms, including Doppler shifts, resonant interactions, specific instability mechanisms and kinetic effects. The specific nature and magnitude of the instability depend on the characteristics of the system and the interplay between the streaming particles and the acoustic waves.

    A comparison of the plots of growth rate of twisted IA waves for (a) single species of electrons and (b) two species of electrons is presented in figure 5. The significant impact of an extra species of electrons on the dynamics of twisted IA waves is observed. An enhancement in the magnitude of unstable wave is identified in the presence of two species electrons. It demonstrates that twisted IA waves are more unstable in a two-species electron plasma than in a single-species electron plasma.

    Figure  5.  Contours of instability growth rate “γia” against wave number and azimuthal quantum number (l) via twist parameter β for (a) single species of electron and (b) two species of electrons.

    To conclude, we have investigated the instability of IA waves with finite OAM states in a plasma which is collisionless and magnetic-field free. The solutions of the model equations of kinetic theory are obtained with Laguerre-Gaussian (LG) profiles of perturbed quantities. The frequency and growth rate of IA plasma waves with non-zero OAM states were obtained and numerically analyzed. It is highlighted that the wave characteristics and instability of low frequency IA twisted wave are significantly modified with finite OAM states via twist parameter β(= k/l{q}_{\varphi }). The analysis illustrates a clear threshold at β = 10, signifying that when the value of β surpasses this threshold, the IA wave experiences Landau damping. Conversely, when β is below this critical value, the IA wave exhibits instability, characterized by significant growth in its amplitude.

    Moreover, the enhanced instability of twisted IA waves is attributed to the increase in electron drift speed, which serves as a recognized source of instability. The results indicate that when the electron drift speed ve0 is at 20.8 times the ion-acoustic phase velocity Cia, the twisted ion-acoustic waves exhibit a higher degree of instability. Furthermore, a comparison of instability growth rate of twisted IA waves for (a) single species of electrons and (b) two species of electrons is presented to showing the significant modification in the magnitude of growth rate in the presence of an extra electron’s species. The presence of angular momentum and co-rotating solar wind associated with plasma species, as identified in the study of comet Halley [22], suggests that there are dynamic plasma phenomena occurring in the vicinity of the comet. These findings open the possibility for further investigations to explore and expand our understanding of the plasma environment surrounding comet Halley. Our results may have relevancy in that environment. Furthermore, the study also has its relevance to the plasma waves propagating in the spinning regions of space plasmas such as those found in the Earth’s magnetosphere, planetary magnetospheres and interplanetary medium.

  • [1]
    DeSilva A W and Stamper J A 1967 Phys. Rev. Lett. 19 1027 doi: 10.1103/PhysRevLett.19.1027
    [2]
    Ghofraniha N et al 2012 Phys. Rev. Lett. 109 243902 doi: 10.1103/PhysRevLett.109.243902
    [3]
    Tsang A C H and Kanso E 2016 Phys. Rev. Lett. 116 048101 doi: 10.1103/PhysRevLett.116.048101
    [4]
    Vlases G C 1964 Phys. Rev. Lett. 12 43 doi: 10.1103/PhysRevLett.12.43
    [5]
    Ehrlich Y et al 1996 Phys. Rev. Lett. 77 4186 doi: 10.1103/PhysRevLett.77.4186
    [6]
    Chu H Y and Lee H K 2011 Phys. Rev. Lett. 107 225001 doi: 10.1103/PhysRevLett.107.225001
    [7]
    Luther B M et al 2004 Phys. Rev. Lett. 92 235002 doi: 10.1103/PhysRevLett.92.235002
    [8]
    Sobral H et al 2000 Appl. Phys. Lett. 77 3158 doi: 10.1063/1.1324986
    [9]
    Dunbar T J and Cetinkaya C 2007 Appl. Phys. Lett. 91 051912 doi: 10.1063/1.2754359
    [10]
    Yi R X et al 2016 J. Anal. At. Spectrom. 31 961 doi: 10.1039/C5JA00500K
    [11]
    Kim J U and Suk H 2002 Appl. Phys. Lett. 80 368 doi: 10.1063/1.1435807
    [12]
    Sueda T, Katsuki S and Akiyama H 1996 Appl. Phys. Lett. 68 1766 doi: 10.1063/1.116660
    [13]
    Hickstein D D et al 2014 Phys. Rev. Lett. 112 115004 doi: 10.1103/PhysRevLett.112.115004
    [14]
    Krall N A and Book D 1969 Phys. Rev. Lett. 23 574 doi: 10.1103/PhysRevLett.23.574
    [15]
    Saeki K and Ikezi H 1972 Phys. Rev. Lett. 29 253 doi: 10.1103/PhysRevLett.29.253
    [16]
    Ehrlich Y et al 1994 Appl. Phys. Lett. 64 3542 doi: 10.1063/1.111970
    [17]
    Kukhlevsky S V et al 1999 Appl. Phys. Lett. 74 2779 doi: 10.1063/1.124011
    [18]
    Mishin G I 1997 Appl. Phys. Lett. 71 49 doi: 10.1063/1.119484
    [19]
    Winfrey A L et al 2012 IEEE Trans. Plasma Sci. 40 843 doi: 10.1109/TPS.2011.2179985
    [20]
    Winfrey A L, Mittal S and Abd Al-Halim M A 2015 IEEE Trans. Plasma Sci. 43 2104 doi: 10.1109/TPS.2015.2417878
    [21]
    Winfrey A L et al 2015 IEEE Trans. Plasma Sci. 43 2195 doi: 10.1109/TPS.2015.2417879
    [22]
    Alimi R and Berdichevsky V 2008 Prop. Explos. Pyrotech. 33 118 doi: 10.1002/prep.200800215
    [23]
    Alimi R and Berdichevsky V 2007 Prop. Explos. Pyrotech. 32 514 doi: 10.1002/prep.200700047
    [24]
    Hang Y H et al 2016 J. Phys. D: Appl. Phys. 49 075201 doi: 10.1088/0022-3727/49/7/075201
    [25]
    Kappen K and Beyer R 2003 Prop. Explos. Pyrotech. 28 32 doi: 10.1002/prep.200390004
    [26]
    Xiao Z G, Ying S J and Xu F M 2015 Prop. Explos. Pyrotech. 40 484 doi: 10.1002/prep.201400215
    [27]
    Wang Q et al 2019 J. Phys. D: Appl. Phys. 52 334002 doi: 10.1088/1361-6463/ab2437
    [28]
    Zhang J B et al 2018 J. Phys. D: Appl. Phys. 51 235204 doi: 10.1088/1361-6463/aac0b0
    [29]
    Hankins O E et al 1993 IEEE Trans. Magn. 29 1158 doi: 10.1109/20.195745
    [30]
    Hankins O E, Bourham M and Mann D 1997 IEEE Trans. Magn. 33 295 doi: 10.1109/20.559972
  • Related Articles

    [1]Dongyu WANG (王栋煜), Ming ZHANG (张明), Shaoxiang MA (马少翔), Shu YANG (杨舒), Kexun YU (于克训), Yuan PAN (潘垣). Design and control of the accelerator grid power supply-conversion system applied to CFETR N-NBI prototype[J]. Plasma Science and Technology, 2020, 22(8): 85601-085601. DOI: 10.1088/2058-6272/ab8d9c
    [2]ZENG Wubing(曾武兵), DING Yonghua(丁永华), YI Bin(易斌), XU Hangyu(许航宇), RAO Bo(饶波), ZHANG Ming(张明), LIU Minghai(刘明海). New Current Control Method of DC Power Supply for Magnetic Perturbation Coils on J-TEXT[J]. Plasma Science and Technology, 2014, 16(11): 1074-1078. DOI: 10.1088/1009-0630/16/11/14
    [3]HUANG Haihong(黄海宏), YIN Ming(殷明), WANG Haixin(王海欣). Design of Controller for New EAST Fast Control Power Supply[J]. Plasma Science and Technology, 2014, 16(11): 1068-1073. DOI: 10.1088/1009-0630/16/11/13
    [4]WU Jinglin(吴景林), LONG Jiaojiao(龙佼佼), LIU Xiaoning(刘小宁). Study on a Highly Stabilized Power Supply for Hybrid-Magnet Superconducting Outsert[J]. Plasma Science and Technology, 2014, 16(9): 890-896. DOI: 10.1088/1009-0630/16/9/15
    [5]HUANG Haihong(黄海宏), YAN Teng(晏腾), WANG Haixin(王海欣). Application of a Current and Voltage Mixed Control Mode for the New Fast Control Power Supply at EAST[J]. Plasma Science and Technology, 2014, 16(4): 420-423. DOI: 10.1088/1009-0630/16/4/22
    [6]LIU Hui (刘辉), TANG Ke (唐柯), GAO Ge (高格), FU Peng (傅鹏), et al.. Study of the EAST Fast Control Power Supply Based on Carrier Phase-Shift PWM[J]. Plasma Science and Technology, 2013, 15(9): 950-954. DOI: 10.1088/1009-0630/15/9/22
    [7]K. SHIMADA, T. TERAKADO, K. YAMAUCHI, M. MATSUKAWA, O. BAULAIGUE, R. COLETTI, A. COLETTI, L. NOVELLO. Minimization of Reactive Power Fluctuation in JT-60SA Magnet Power Supply[J]. Plasma Science and Technology, 2013, 15(2): 184-187. DOI: 10.1088/1009-0630/15/2/22
    [8]ZHU Zhe (朱哲), ZHU Yinfeng (朱银锋), HUANG Ronglin (黄荣林), FU Peng (傅鹏), DING Yixiao(丁逸骁). Study on the Current-sharing Control System of the TF Power Supply for a Superconducting Tokamak[J]. Plasma Science and Technology, 2012, 14(10): 941-946. DOI: 10.1088/1009-0630/14/10/16
    [9]CHEN Wenguang (陈文光), RAO Jun (饶军), LI Bo (李波), LEI Guangjiu (雷光玖), CAO Jianyong (曹建勇), WANG Mingwei (王明伟), KANG Zihua (康自华), FENG Kun (冯鲲), HL-A NBI Group. Technical Design of Arc-Discharge and Deceleration Power Supply for MW Level NBI System on HL-2A Tokamak[J]. Plasma Science and Technology, 2012, 14(10): 936-940. DOI: 10.1088/1009-0630/14/10/15
    [10]XU Weidong (徐伟东), XUAN Weimin (宣伟民), YAO Lieying (姚列英), WANG Yingqiao (王英翘). Development of 8 MW Power Supply Based on Pulse Step Modulation Technique for Auxiliary Heating System on HL-2A[J]. Plasma Science and Technology, 2012, 14(3): 263-268. DOI: 10.1088/1009-0630/14/3/14
  • Cited by

    Periodical cited type(3)

    1. Yang, Z., Huang, C., Luo, Y. et al. Regulation of the pantograph-catenary arc by an alternating transverse magnetic field. Plasma Science and Technology, 2025, 27(3): 035402. DOI:10.1088/2058-6272/adab4a
    2. Rath, S., Kar, S. Microwave atmospheric pressure plasma jet: A review. Contributions to Plasma Physics, 2025, 65(2): e202400036. DOI:10.1002/ctpp.202400036
    3. Suresh, I., Srikar, P.S.N.S.R., Gangwar, R.K. Integration of ML methods with CR model-based optical diagnostic for the estimation of electron temperature in Ga laser produced plasma. Physics of Plasmas, 2024, 31(11): 113501. DOI:10.1063/5.0223030
    1. Yang, Z., Huang, C., Luo, Y. et al. Regulation of the pantograph-catenary arc by an alternating transverse magnetic field. Plasma Science and Technology, 2025, 27(3): 035402. DOI:10.1088/2058-6272/adab4a
    2. Rath, S., Kar, S. Microwave atmospheric pressure plasma jet: A review. Contributions to Plasma Physics, 2025, 65(2): e202400036. DOI:10.1002/ctpp.202400036
    3. Suresh, I., Srikar, P.S.N.S.R., Gangwar, R.K. Integration of ML methods with CR model-based optical diagnostic for the estimation of electron temperature in Ga laser produced plasma. Physics of Plasmas, 2024, 31(11): 113501. DOI:10.1063/5.0223030

    Other cited types(0)

Catalog

    Figures(15)  /  Tables(2)

    Article views (72) PDF downloads (37) Cited by(3)

    /

    DownLoad:  Full-Size Img  PowerPoint
    Return
    Return