BRIEF RESEARCH REPORT article
Sec. Low-Temperature Plasma Physics
Volume 10 - 2022 | https://doi.org/10.3389/fphy.2022.951960
On the onset of breathing mode in Hall thrusters and the role of electron mobility fluctuations
- 1Department of Civil and Industrial Engineering, University of Pisa, Pisa, Italy
- 2SITAEL S.p.A., Pisa, Italy
- 3Department of Physics, University of Pisa, Pisa, Italy
Breathing mode is an ionization instability which is observed ubiquitously in the operation of Hall thrusters. It is recognized as a relatively low frequency (10–30 kHz) longitudinal oscillation of the discharge current and the plasma parameters. Although breathing instability is widely studied in the literature, the conditions for its origin are not fully understood. In this work we investigate the mechanisms responsible for the origin of the breathing mode in Hall thrusters by using a numerical model, allowing us to highlight the importance of electron mobility fluctuations for the onset and self-sustenance of the instability. Our one-dimensional, fully fluid model of the thruster channel is calibrated against the measured discharge current signal for a 5 kW-class Hall thruster operating in a condition where breathing mode is fully developed. The corresponding steady, unstable configuration (base state) is numerically computed by applying the Selective Frequency Damping (SFD) method. Then, a series of numerical tests is performed to show the existence of a feedback loop involving fluctuations around the base state of the neutral density, electron mobility, and electric field. We show that oscillations of the electron mobility are mainly caused by variations of the neutral density and are in phase with them; this, in turn, induces oscillations of the electric field, which are in phase opposition. The electric field acts simultaneously on the electron temperature and on the ion dynamics, promoting the depletion and replenishment of neutrals in the chamber.
Hall thrusters are plasma-based electric propulsion devices for space applications. The cross-field configuration of these devices allows for an efficient (and synergic) ionization and acceleration of the propellant, without space-charge limiting effects that are typical of a configuration based on electrostatic acceleration. However, the strong coupling between the processes of ionization and acceleration is associated with high-amplitude, low frequency (5–30 kHz) discharge current oscillations, the so-called breathing mode. This mode has detrimental effects on the thruster performance and can cause significant coupling issues with the driving electronics. Additional details on the design and operation of Hall thrusters can be found, for instance, in the work of .
Several theoretical and experimental research efforts have been dedicated to the investigation of breathing mode. From the experimental point of view, breathing mode is directly diagnosed by sampling at high frequency the thruster discharge current. Measurements performed with intrusive [2–4] and non-intrusive [5–8] diagnostics have shown that breathing oscillations manifest as longitudinal coherent waves propagating in the plasma domain without significant frequency dispersion. From the modelling point of view, investigations were performed by means of zero-dimensional models [9–11] as well as with multidimensional simulations . Zero-dimensional simulations have demonstrated that the origin of the breathing oscillations is associated with the periodic ionization and replenishment cycle of the propellant . However, the proposed mechanisms could not reproduce self-sustained oscillations . One-dimensional fluid models, on the other hand, have shown to be able to simulate the breathing mode. The authors of  employed a one dimensional formulation and selectively fixed in time different contributors to the discharge dynamics. Their results, for an uncalibrated case, suggest the importance of the electron temperature oscillations and of the non-linearity in the electron power absorption on the onset of breathing mode. As described in , these models can be specifically calibrated using the discharge current signal, in which the breathing mode oscillations are the most distinctive feature, indicative of the thruster operative condition. In the present work, the calibrated model of  is used to further investigate the mechanisms responsible for the onset of such self-sustained breathing oscillations in Hall thrusters.
To describe the thruster discharge we consider an unsteady 1D fluid model . The domain extends along the thruster symmetry axis z, from z = 0 (anode) to z = zf = 2 Lc, where Lc is the length of the thruster channel. Sensitivity of the results to different choices of the domain length has been verified and it is negligible for the results presented here. The xenon plasma is assumed as quasi-neutral (ions and electrons have equal density n) and involves three species, neutrals, singly-charged ions and electrons. It is assumed that neutrals have a constant axial velocity, ions are unmagnetized and cold, and electrons inertia is neglected, thus relying on a drift-diffusion approximation. Finally, the external magnetic field is assumed to be purely radial, so that the model may not be suited for investigating modern Hall thrusters employing magnetic shielding. The plasma dynamics is described by the following set of equations:
where n, u and T stand for density, axial velocity and temperature, and kB is the Boltzmann constant. Subscripts n, i and e indicate neutrals, ions and electrons, respectively. The discharge current density J is computed by integrating Eq. 4 between the anode and the end of the domain. The resulting integral equation contains the applied potential difference ΔV, which is specified in the operating condition. The ionization-rate coefficient kI and the collisional energy loss coefficient K are computed as a function of the electron internal energy using the LXCat database  and the Bolsig+  solver. Bolsig+ is used to parse the cross-section data to compute Maxwellian rates. The terms
where μ0 = e/(νe me) is the unmagnetized mobility and Ω = ωe/νe is the Hall parameter, ωe = (eBr)/me being the local cyclotron frequency of electrons and Br being the external magnetic field intensity. The term νe represents the momentum transfer collision frequency and is given by the sum of three contributions:
where νew is the electron-wall collision frequency, proportional to a wall interaction coefficient and the Bohm velocity, and νc takes into account collisions between electrons and neutral atoms (elastic, excitation and ionization collisions) and it is given by the product of the neutral density and the relevant reaction rate kx computed using LXCat and Bolsig+. Coulomb electron-ion (e-i) collision frequency has not been accounted for after having verified that it is always negligible in comparison to νe. More specifically, for the investigated case the e-i collision frequency is one order of magnitude lower than νc in the ionization and acceleration regions, where the core of the breathing instability is supposed to be localized. The comparison against the experiments of the calibrated model further support this and the other modeling assumptions. However, we underline that in different thrusters the proposed modeling assumptions should be critically reconsidered. Finally, the anomalous collision frequency νa is modeled as
where β(z), which is time independent, is the same function of the axial position already used in . Several attempts to provide a closure function for the anomalous collision frequency, i.e., to identify a relation between β and the local plasma fields, have been performed (see, for instance ), but to date no self-consistent closure function exists. Even if in specific cases different spatial profiles have been adopted (see [18, 19]), the particular choice followed in our model, as shown in , leads to reasonable agreement against the experiments for the considered thruster.
For the objectives of the present work it is important to underline that, in our model, the mobility μ depends only on two plasma quantities, the neutrals density and the electron temperature; this aspect is made explicit, when needed for clarity, by writing μ as μ(nn, Te). Finally, we assume that the electron population follows an isotropic Maxwellian distribution with temperature Te. Thus, the electron pressure pe is related to Te as follows:
As concerns the boundary conditions, at the inner boundary (z = 0) the neutral density is fixed based on the injected mass flow rate and accounting for the ion recombination at the anode surface. Additionally, ions must enter the anode sheath edge with a velocity greater than or equal to the local sound speed (for the electron-repelling sheath to form), and the electron energy flux is imposed according to the analysis of a sheath in front of a biased electrode. At the outer boundary (z = zf), instead, we fix Te = 2eV. For additional details on the adopted boundary conditions and on the finite-volume method used to discretize in space the PDEs of the model refer to .
The objective of the present analysis is to investigate the mechanisms at the basis of the self-sustenance of the breathing mode using the model described above. To this purpose a particular configuration is selected in which breathing is evident and well observable. The considered thruster is a 5 kW-class Hall thruster, SITAEL’s HT5k, operating with a conventional magnetic field topology (almost radial magnetic lines at the channel exit) at nearly 2.6 kW, with a discharge voltage of 300 V and a mass flow rate of 8 mg/s. This particular configuration has been investigated experimentally and detailed measurements of the plasma properties are available . Some details on the magnetic topology can be found in .
In order to increase the realism of the numerical model, a calibration has been carried out against the discharge current signal available from the experiments, and the resulting estimated properties of the plasma have been compared against the available measurements . More precisely, the neutral velocity, wall interaction coefficient, and anomalous collision frequency profile were progressively varied until the simulated discharge current signal matched as close as possible the measured one in terms of dominant frequency, mean value, and amplitude of the oscillations. The model resulting from the calibration and its accuracy in representing the plasma properties are discussed in details in . We underline that different calibration strategies may be applied, for instance by taking into account local measurements of plasma properties. The one adopted here has the advantage of relying only on a single measurable quantity which is usually available in existing experimental databases. Moreover, although the kind of calibration may affect the resulting model, the core mechanism leading to breathing mode should not be significantly affected because it is by far the dominant dynamical process at play in the selected calibration database. For this reason, in the present work we investigate breathing mode by using the calibrated model of .
Since the breathing mode has the characteristics of a self-sustained instability, the first step is to find the steady (unstable) configuration on top of which it develops. This particular configuration, denoted in the following as base configuration, is computed by searching for a steady solution of the calibrated model. This is achieved numerically by applying the Selective Frequency Damping (SFD) method, which is commonly used to find steady solutions in high-dimensional dynamical systems, as for instance in the field of computational fluid dynamics [21, 22]. This technique relies on the application of a proportional feedback control which vanishes on the steady state and makes it stable, so that it can be computed by classical time-advancing methods. In the present case the controller is applied only to the neutral continuity equation, given the importance of neutral density oscillations in the breathing mode dynamics. Equation 1 becomes
where χ ≃ 107 Hz,
and Δ ≃ 10−4 s. Equation 12 updates the target solution
FIGURE 1. Discharge current signal of the calibrated model with (red) and without (blue) the application of the SFD technique.
We now investigate the self-sustenance process driving the instability by means of ad-hoc numerical simulations. We start by calculating the base configuration; then, starting from the base configuration, we have verified that the instability spontaneously emerges, leading to an exponential growth of its amplitude. The instability then saturates in the non linear phase finally settling on a periodic limit cycle (Test T1, Figure 2A). This evidence supports the idea that the breathing mode is a linear Hopf-like instability of the base configuration.
FIGURE 2. Simulated discharge current signal for tests (A) T1, (B) T2, (C) T3 (D) T4 starting from a perturbation of the base configuration and tests (E) T2 and (F) T3 starting from the saturated breathing instability; results from tests T6 and T7 are reported in (G) and (H), respectively.
We now focus on the role of mobility on the self-sustenance mechanism of the breathing mode. The idea that mobility can play a dominant role in the instability originates from the observation that breathing mode can be suppressed, both in experiments [8, 23] and in the present model (not shown here), by increasing the intensity of the magnetic field. Note that in the model, Br enters only via the mobility μ.
Considering now the calibrated case, a second artificial test (T2) has been carried out by starting from the base configuration with the spatial distribution of μ artificially fixed to the value computed on the base configuration, i.e., μ = μbc = μ(nn,bc, Te,bc) (where the subscript bc indicates quantities taken in the base configuration). In this way we avoid that μ varies due to plasma fluctuations
Considering that nn,bc and Te,bc are fixed quantities, the mobility
A fourth test T4 is carried out by considering only variations of mobility induced by fluctuations of neutral density,
In order to further investigate this point, a simplified model for μ is derived by fixing Te = Te,bc and linearizing its dependence on
where km is the sum of the reaction rates for momentum transfer collisions between electrons and neutrals.
A fifth test, T5, is carried out using
FIGURE 3. Discharge current signal for test T5 (red line); green and blue lines are the envelopes of the discharge current signals obtained by imposing
With the support of the tests described above, it is now possible to propose a physical interpretation concerning the feedback mechanism that leads to the breathing instability in the investigated case. To this purpose we look at variations around the base configuration and we focus on the region where ionization is concentrated. From Eq. 4 we note first that the electric field
When the neutral density decreases, the same happens for the electron mobility and, through the first term on the right side of Eq. 14, this provides a positive contribution to the electric field. As a matter of fact, during the simulated oscillations, the electric field reaches its peak when the neutral density is at its minimum, and the two quantities oscillate in phase opposition, as depicted in Figure 4A. This phenomenon promotes two parallel effects:
(E1) As concerns Te, when the neutral density decreases and the electric field increases, the work performed by the electric field on the electron fluid also increases, inducing an increase of electron temperature and, in turn, of ionization rate (see Figure 4B). This creates a positive feedback as it contributes to further decrease the neutral density.
(E2) As concerns the ion dynamics, when replenishment of neutrals begins, the electric field is at its maximum and ions are accelerated towards the exit. As a consequence, during replenishment plasma density decreases and this leads to a decrease of the ionization rate, thus promoting the channel replenishment and creating again a positive feedback on the neutral density.
FIGURE 4. Oscillations in time of neutral density (blue) and electric field (red) around their base state values at point z = Lc/2 (A) at incipient instability and (B) of neutral density and temperature (black) on the limit cycle.
The electron temperature fluctuations play an intricate role in the system dynamics, with possible feedback on all species. The performed tests have indicated that their contribution to the variations of mobility are not relevant for the development nor for the sustenance of breathing mode oscillations. Nevertheless, to assess the importance of temperature fluctuations
In this work the mechanisms leading to the onset of breathing mode in Hall thrusters are further investigated numerically and theoretically. A one-dimensional fluid model of the thruster channel is used, which is calibrated and validated against the experimental data available for SITAEL’s HT5K thruster. A crucial element for the proposed analysis is the identification of the base configuration, i.e., the unstable steady configuration of the thruster which undergoes breathing instability, so that breathing can be investigated both during its initial onset and in its final saturated state. Results demonstrate that fluctuations of the electron mobility around its distribution in the base configuration are necessary for the onset of the breathing mode. In particular, a linear dependence between oscillations of neutral density and electron mobility is sufficient for triggering the instability in the proposed model, inducing fluctuations of the electric field in phase opposition. This establishes a positive feedback on the ionization rate via two different mechanisms, one involving the electron energy and one, which is dominant in the case analyzed, involving the ion dynamics. These observations provide a new perspective in the interpretation of the ionization instability at the basis of breathing mode.
Data availability statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
Conceptualization, LL, SC, VG, and TA; methodology, LL, SC, and VG; software, VG, MMS, LL, and SC; analysis, LL, VG, SC, TA, and FC; resources, SC, TA, and FC; writing—original draft preparation, SC, LL, VG, and TA; writing—review and editing, LL, SC, VG, FC, and TA. All authors have read and agreed to the published version of the manuscript.
The work described in this paper has been funded by the European Union under the H2020 Programme ASPIRE-GA 101004366.
Conflict of interest
Authors VG and TA were employed by the company SITAEL S.p.A.
The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
1. Boeuf J-P. Tutorial: Physics and modeling of Hall thrusters. J Appl Phys (2017) 121:011101. doi:10.1063/1.4972269
2. Giannetti V, Saravia MM, Andreussi T. Measurement of the breathing mode oscillations in Hall thruster plasmas with a fast-diving triple Langmuir probe. Phys Plasmas (2020) 27. doi:10.1063/5.0022928
3. Lobbia R, Gallimore A. A method of measuring transient plume properties. In: proceedings of the 44th AIAA/ASME/SAE/ASEE Joint Propulsion Conference & Exhibit (2008). doi:10.2514/6.2008-4650
4. Lobbia RB, Gallimore AD. High-speed dual Langmuir probe. Rev Scientific Instr (2010) 81:073503. doi:10.1063/1.3455201
5. Chaplin VH, Lobbia RB, Lopez Ortega A, Mikellides IG, Hofer RR, Polk JE, et al. Time-resolved ion velocity measurements in a high-power Hall thruster using laser-induced fluorescence with transfer function averaging. Appl Phys Lett (2020) 116:234107. doi:10.1063/5.0007161
6. Dale ET, Jorns BA. Non-invasive time-resolved measurements of anomalous collision frequency in a Hall thruster. Phys Plasmas (2019) 26:013516. doi:10.1063/1.5077008
7. Fabris AL, Young CV, Cappelli MA. Time-resolved laser-induced fluorescence measurement of ion and neutral dynamics in a Hall thruster during ionization oscillations. J Appl Phys (2015) 118:233301. doi:10.1063/1.4937272
8. Sekerak MJ, Gallimore AD, Brown DL, Hofer RR, Polk JE. Mode transitions in Hall-effect thrusters induced by variable magnetic field strength. J Propulsion Power (2016) 32:903–17. doi:10.2514/1.B35709
9. Dale ET, Jorns BA. Two-zone Hall thruster breathing mode mechanism, Part I: Theory. In: 36th International Electric Propulsion Conference (2019). p. 1.
10. Fife JM, Martinez-Sanchez M, Szabo J. A numerical study of low-frequency discharge oscillations in Hall thrusters. In: proceedings of the 33rd Joint Propulsion Conference and Exhibit. American Institute of Aeronautics and Astronautics Inc, AIAA (1997). doi:10.2514/6.1997-3052
11. Hara K, Sekerak MJ, Boyd ID, Gallimore AD. Perturbation analysis of ionization oscillations in Hall effect thrusters. Phys Plasmas (2014) 21:122103. doi:10.1063/1.4903843
12. Hara K. An overview of discharge plasma modeling for Hall effect thrusters. Plasma Sourc Sci Technol (2019) 28:044001. doi:10.1088/1361-6595/ab0f70
13. Lafleur T, Chabert P, Bourdon A. The origin of the breathing mode in Hall thrusters and its stabilization. J Appl Phys (2021) 130. doi:10.1063/5.0057095
14. Giannetti V, Saravia MM, Leporini L, Camarri S, Andreussi T. Numerical and experimental investigation of longitudinal oscillations in Hall thrusters. Aerospace (2021) 8:148. doi:10.3390/aerospace8060148
15. [Dataset] Biagi S. Biagi database (2021). Available from: www.lxcat.net (accessed March 9, 2021).
16. Hagelaar GJM, Pitchford LC. Solving the Boltzmann equation to obtain electron transport coefficients and rate coefficients for fluid models. Plasma Sourc Sci Technol (2005) 14:722–33. doi:10.1088/0963-0252/14/4/011
17. Lafleur T, Baalrud SD, Chabert P. Theory for the anomalous electron transport in Hall effect thrusters. II. kinetic model. Phys Plasmas (2016) 23. doi:10.1063/1.4948496
18. Hara K, Mikellides IG. Characterization of low frequency ionization oscillations in Hall thrusters using a one-dimensional fluid model. In: 2018 Joint Propulsion Conference. American Institute of Aeronautics and Astronautics (2018). doi:10.2514/6.2018-4904
19. Lopez Ortega A, Mikellides IG, Sekerak MJ, Jorns BA. Plasma simulations in 2-d (r-z) geometry for the assessment of pole erosion in a magnetically shielded Hall thruster. J Appl Phys (2019) 125:033302. doi:10.1063/1.5077097
20. Andreussi T, Giannetti V, Leporini A, Saravia MM, Andrenucci M. Influence of the magnetic field configuration on the plasma flow in Hall thrusters. Plasma Phys Control Fusion (2018) 60:014015. doi:10.1088/1361-6587/aa8c4d
21. Jordi B, Cotter C, Sherwin S. Encapsulated formulation of the selective frequency damping method. Phys Fluids (2014) 26:034101. doi:10.1063/1.4867482
22. Åkervik E, Brandt L, Henningson D, Hæpffner J, Marxen O, Schlatter P, et al. Steady solutions of the Navier-Stokes equations by selective frequency damping. Phys Fluids (2006) 18:068102. doi:10.1063/1.2211705
23. Giannetti V, Piragino A, Paissoni CA, Ferrato E, Estublier D, Andreussi T, et al. Experimental scaling laws for the discharge oscillations and performance of Hall thrusters. J Appl Phys (2022) 131. doi:10.1063/5.0070945
24. Lopez Ortega A, Mikellides IG, Chaplin V. Numerical simulations for the assessment of erosion in the 12.5-kW Hall effect rocket with magnetic shielding (HERMeS). In: proceedings of the 35th International Electric Propulsion Conference. IEPC-2017-154 (2017).
Keywords: plasma physics, low-temperature plasmas, electric propulsion, Hall thrusters, breathing mode, ionization instability, electron mobility
Citation: Leporini L, Giannetti V, Saravia MM, Califano F, Camarri S and Andreussi T (2022) On the onset of breathing mode in Hall thrusters and the role of electron mobility fluctuations. Front. Phys. 10:951960. doi: 10.3389/fphy.2022.951960
Received: 24 May 2022; Accepted: 12 July 2022;
Published: 09 August 2022.
Edited by:Francesco Taccogna, National Research Council (CNR), Italy
Reviewed by:Mirko Magarotto, University of Padua, Italy
Stefano Boccelli, University of Ottawa, Canada
Alejandro Lopez Ortega, NASA Jet Propulsion Laboratory (JPL), United States
Copyright © 2022 Leporini , Giannetti , Saravia , Califano , Camarri and Andreussi . This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: S. Camarri, email@example.com