Next Article in Journal
Challenges of Ablatively Cooled Hybrid Rockets for Satellites or Upper Stages
Previous Article in Journal
A Review on the Current Status of Icing Physics and Mitigation in Aviation
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Ion Source—Thermal and Thermomechanical Simulation

by
Victoria V. Svotina
*,
Andrey I. Mogulkin
and
Alexandra Y. Kupreeva
Research Institute of Applied Mechanics and Electrodynamics of Moscow Aviation Institute, 125080 Moscow, Russia
*
Author to whom correspondence should be addressed.
Aerospace 2021, 8(7), 189; https://doi.org/10.3390/aerospace8070189
Submission received: 20 May 2021 / Revised: 24 June 2021 / Accepted: 7 July 2021 / Published: 14 July 2021

Abstract

:
The main purpose of this work is to conduct ground development testing of the ion source intended for use the space debris contactless transportation system. In order to substantiate the operating capability of the developed ion source, its thermal and thermomechanical simulation was carried out. The ion source thermal model should verify the ion source operating capability under thermal loading conditions, and demonstrate the conditions for ion source interfacing with the systems of the service spacecraft with the ion source installed as a payload. The mechanical and mathematical simulation for deformation of the ion source ion-extraction system profiled electrodes under thermal loading in conjunction with the prediction of the strained state based on the numerical simulation of the ion source ion-extraction system units, making it possible to ensure the stability of the ion source performance. Good agreement between the thermal and thermo-mechanical ion source simulation results and experimental data has been demonstrated. It is shown that the developed ion source will be functional in outer space and can be used as an element of the space debris contactless transportation system into graveyard orbits.

1. Introduction

Non-functioning spacecraft (SC), or their fragments are understood to be space debris (SD). The problem of near-Earth space pollution by SD has acquired an exceptional urgency at the present time as it significantly complicates and encumbers further space-related activities [1]. As the simulation of the processes of near-Earth space operations and pollution shows, an annual active removal of 5–20 large-sized objects from near-Earth space can result in reducing the growth of space debris population caused by collisions and fragmentation of available, both functioning and non-functioning, space objects [2]. By now, a large number of different methods to remove SD from the vicinity of the geostationary orbit and from low Earth orbits have been proposed [1,3,4,5,6,7,8,9,10,11,12,13,14,15,16]. However, the optimum way to clear the near-Earth space has not yet been found while technological aspects to solve the problem of near-Earth space pollution are under development and demonstration testing.
The service spacecraft (SSC) that approach the SD to be removed and that tug it into a graveyard orbit by means of mechanical contact or by contactless impact can be used for SD removal. A large-sized SD can have great residual angular momenta, or can be tumbling, making it difficult to remove it by contact methods, especially in cases with large solar panels and/or uncontrolled high angular velocities. These circumstances increase the risks of destruction for SSC structural elements or SD to be removed, which can lead to an increase in the number of SD fragments with a simultaneous reduction in their size, and, consequently, to growing difficulties in their detection and control. In this regard, contactless removal methods are of great interest.
Transportation of space objects by an ion beam was proposed in order to perform transportation in the “Earth-Moon” system [17], and to remove space debris [18,19,20]. This SD removal method was named the Ion Beam Shepherd. The use of anion source (IS) for transportation of dead SC or SD appears to be the most promising due to the SD material ablation observed during SD removal by IS beam.
It is necessary to show the IS sustainable operating capability in outer space under the conditions of thermal and thermomechanical loading to assess the possibility and efficiency of contactless SD transportation method.

2. The Status of the Problem

The IS thermal model development, taking into account the requirements of reliable functioning and integration with the SSC systems, is considered to be one of the most challenging research and engineering issues in the development of a system for contactless transportation in space [21,22,23,24,25,26,27,28,29].
The complexity of the problem is in the need to develop reliable computational models for IS thermo-physical processes that would take processes of different natures into account, including plasma processes, driven by ionization of the plasma-forming gas in a self-sustained inductive discharge, taking into account thermal dissipation characteristics during interaction of plasma with the SSC structural elements, as well as the heterogeneity of the materials used in the IS structure.
Structurally, the IS ion-extraction system (IES) is the most complex unit. It is necessary to consider the electrode performance features while developing the IES design [30,31,32]. This is driven by the complexity in providing a gap between the IES electrodes. IES electrodes are most often made in the form of thin-walled slightly sloping rotary shells with close-packed perforation. Peculiarities of IES electrode operation include uneven heating along the radius and thickness, which creates prerequisites for the reduction of output IES performance and of those of IS as a whole, and also enhances the likelihood of electrical breakdown in the inter-electrode gap. As a result, internal temperature stresses arise even in freely fixed electrodes, which inevitably leads to deformation of the electrodes.
One should expect the occurrence of additional deflections of the IES electrodes taking the following into account: a difference in electrode thicknesses, usage of various construction materials, and operation under conditions of temperature gradients. There is a need to solve the problem of ensuring stable operation of the IES electrodes and their stability under nonlinear bending. The problem is further complicated by the tendency for the electrode thickness to be decreased down to 0.3 mm.
In order to solve the problematic issues indicated, the task was set to develop an IES thermo-mechanical model to account for the thinness of walls and non-uniform thickness of electrodes, the effect of perforation, uneven heating along the radius of the electrodes; as well as the task to develop an IS thermal model taking into account the usage of various materials and operation in heterogeneous thermal load conditions. During ground experimental development, it is necessary to demonstrate operating capability of the IS produced on the basis of simulation results, and to show consistency between simulation results and experimental data.
In order to implement the contactless transportation method, force coupling between the SSC and SD is ensured using a quasi-neutral ion beam injected by the IS from the SSC. The “SSC-SD” virtual assembly orbital motion is controlled with the help of a compensating propulsion system. The diagram of propelling force action arising during transportation of SD by an ion beam is shown in Figure 1. The compensating propulsion system comprises two thrusters mounted on a two-axis gimbal unit. This configuration allows us to change thrust of the compensating propulsion system (FEP) and its projections on mutually perpendicular axes by the thruster gimballing.
The ion beam acts on the SD by a momentum transfer to it. In this case, the SSC must escort the SD at some distance (d). To that end, the SSC and the SD average accelerations must be equal, which can be achieved under the following condition: F E P = ( 1 + k M S C M S D ) F I S , where: MSC is the SSC mass; MSD is the SD mass; k = Fd/FIS ≤ 1 is a coefficient equal to the ratio of momentum Fd transferred by the IS ion beam to the SD, to the momentum generated by the IS. The contactless method efficiency is largely determined by the k coefficient value. Many publications have been devoted to consideration of the quasi-neutral ion beam behavior, for example, [33].
The current surface that limits 95% of the ion current is generally considered to be the beam boundary. Dimensions of the near zone are 2–3 diameters of the ion beam at the IS exit plane. Ion beam characteristics at a distance equal to 2–3 diameters of the ion beam from the IS exit plane are defined by the conditions for the ion beam formation by the IS electrostatic system. Distribution of the ion current density along the radius and angle in the IS exit plane is an initial condition to analyze the ion beam behavior in the near zone. Similar characteristics at the outer boundary of the near zone are the initial parameters for considering the problem regarding the ion beam evolution during its further propagation in space. Further behavior of the ion beam is dictated by the presence of the so-called ambipolar electric field in the beam. This field emerges due to the presence of electrons in the beam that are introduced into the beam by the neutralizer (in flight conditions), or come from the vacuum chamber walls (during ground development testing) to offset the volume charge of ions in the beam. In this case, the condition of equality for electron and ion densities is locally satisfied in the volume of the beam: ne = ni, and the condition of equality for the ion and electron currents flowing out through the cross section at the near zone outer boundary is metal so: Ii = Ie.

3. Balance of Particles and Power in the Discharge

The IS distinguishing characteristic is that the only thermal power source is the electromagnetic field initiating an inductive discharge. The main processes that stipulate the IS thermophysical state are the following:
-
Absorption of HF power by plasma due to plasma electrical conductivity;
-
Plasma formation driven by ionization and heating of electrons and leading to an increase in the ionized propellant gas internal energy;
-
Carry-over of a portion of plasma power by ions and electrons onto the walls, including ionization energy released during recombination of ions and electrons of the propellant gas, which is also converted into heat.
The IS operation principle is that high-frequency (HF) currents in the inductor generate a HF magnetic field with induction B H F = 0.5 B H F 0 sin ω t in the IS gas discharge chamber (GDC) volume. In turn, the magnetic field induces a HF electric field with intensity E I N D = 0.5 B H F 0 r ω cos ω t where ω is HF generator frequency.
Electrons oscillating with the electric field frequency are accelerated under the action of the electric field. In the case of inelastic collisions of electrons with propellant gas atoms or ions, their excitation and/or ionization occurs. Simultaneously, as a result of ionization, free electrons are generated, while some electrons escape the discharge and deposit on the surfaces limiting the discharge (the GDC walls and/or IES electrodes). The RF discharge in the IS is inductive and self-sustained. It does not require a permanently operating cathode as a source of electrons inside the IS GDC. After feeding the propellant gas and supplying power from the RF generator, the initial pulse injection of electrons emitted by an external source of electrons (for example, a neutralizer) can be used to ignite the discharge.
Intensities of electric and magnetic fields are relatively low. Consequently, due to the significant mass of propellant gas ions, it can be assumed that the electric and magnetic fields do not act on the propellant gas ions. Ions and atoms have a low temperature governed by the GDC wall temperature (on the order of 250–300 °C).
In accordance with the law of J.C. Maxwell and L.E. Boltzmann, the main mechanism to establish an equilibrium electron energy distribution are electron–electron collisions (thermalization). In this case, cold electrons formed as a result of an inelastic collision acquire the temperature of plasma electrons heated by an HF electric field. According to [34], electrons in plasma have the Maxwellian distribution with temperature Te and density exceeding ne = 1017 m−3. Te = 3–4 eV in the inductive discharge volume, and Te = 5–7 eV in the near-wall discharge area. The electron temperature weakly depends on the value of HF power supplied to the discharge [35,36,37,38].
A relatively low temperature of electrons, on one hand, leads to an increase in specific energy consumption to ionize the propellant gas within the volume, and on the other hand, it causes a low plasma potential relative to the GDC walls, and, consequently, a low energy of ions depositing onto the GDC inner surface and the surface of the IES screen grid facing the plasma. Low ion energy does not result in erosion of these surfaces.
The IS plasma is weakly ionized and the ionization degree does not exceed 5%. The fraction of doubly charged and multiply charged ions does not exceed 1%. Their presence is normally neglected when performing thermal simulation. The plasma-forming gas pressure in the GDC is about 10−1 Pa. The representative electron density values are about ne ~ 4 × 1017 m−3. The mean free path of electrons exceeds the characteristic GDC volume dimensions.
It has been experimentally proven that the HF discharge plasma in RIT (radio-frequency ion thruster) is rarefied and non-equilibrium. The photons emitted by the excited particles reach the GDC walls and are absorbed by it without interacting with the propellant gas atoms. The consequence of the non-equilibrium state is that the electron temperature Te is much higher than the temperature of atoms Ta and ions Ti.
In order to describe the plasma state correctly, it is necessary to know spatial distributions for the densities of atoms and electrons (densities of electrons and ions are equal under quasi-neutral conditions), as well as for the electron temperature. Simulation of processes in the GDC is based on the establishment of balance of particles and power in the plasma volume limited by the GDC walls and by the screen electrode surface [39].
In a self-sustained HF discharge, the balance of electrons is defined by the rate of their formation as a result of ionization (the ionization rate) and the rate of their escape (deposition) onto the GDC walls and into the holes in the IES electrodes, since the probability for recombination of ions and electrons in the GDC volume is close to zero. The deposition rate depends on the equilibrium plasma potential relative to the GDC walls, which is established automatically.
The balance model is developed similarly to the one described in [40] with an additional simplification with respect to neglecting the processes of stepwise ionization through metastable states, since the probability of such ionization is low.
The electron energy distribution function in all models is assumed to be the Maxwellian one. Obviously, only the electrons from the Maxwellian distribution tail area have an energy exceeding the atom ionization energy for almost any propellant that can be used in an IS. On the other hand, atom excitation energies are lower than the ionization energy, which leads to the fact that there can be several excitation events per an ionization event [41].
The excitation energy is primarily spent for radiation. The effective energy spent for the production of one ion in the discharge is CID = Vi + CR, where: CR are total radiation losses in all the events of excitation-quenching pertaining to one event of the propellant gas atom ionization; Vi is the first potential of the propellant gas atom ionization.
A relative amount of HF power consumption for the excitation of propellant gas atoms is defined by the ratio of cross sections of ionization and excitation of atomic levels of the propellant gas atoms, which is calculated by integrating over the Maxwellian distribution of the energy transferred by electrons to atoms in inelastic collisions, taking the cross section of atomic level excitation into account.
Specific energy losses for the formation of one propellant gas (Xenon) ion CR in the plasma discharge as a function of electron temperature Te are given in [42]. With Te growth, the specific radiation losses and, accordingly, the cost of an ion in the discharge decrease. Since the electron temperature in an inductive high-frequency discharge is 4–6 eV, the radiation losses amount to 4–6 eV per one ionization event. An increase in specific radiation losses with Te decrease is explained by an increase in the number of excitation events per one ionization event. The formation of a doubly charged xenon ion with the Maxwellian distribution of electrons with a temperature of 5 eV is unlikely. In an HF induction discharge, the fraction of doubly charged ions does not exceed 1% [43], which allows us to discard doubly charged ions in the discharge power balance.
Variations of the main plasma parameters in the discharge along the GDC radius is presented in [43]. The electron density decreases from the GDC axis towards its walls. The electron temperature is maximum on the GDC walls. The ion current density changes insignificantly, which makes it possible to use approximate relations for a homogeneous plasma with electron density and electron temperature being uniformly distributed over the volume. For electron density and electron temperature, their mean values that define integral plasma parameters derived from them (density, ion beam current, near-wall plasma potential drop near the GDC surface and emission electrode) are assumed. The total ion production in a plasma discharge is:
I u = n e n a σ a i V e V d = 0.43 n e 2 e T e M S G D C
where Ve is the electron velocity; σ a i V e is the averaging of the ionization cross section over electron velocities Ve; Vd is the GDC plasma volume; M is the propellant gas atom (ion) mass; SGDC is the internal area of GDC surface surrounding the discharge volume.
The electron balance equation is as follows:
n e n a σ a i V e V d = 0.25 8 e T e π m e exp ( Δ U T e ) S G D C
Here: ∆U is the near-wall plasma potential drop.
Equations (1) and (2) describe the balance between the particles depositing onto the walls (the left-hand sides) per unit time, and the particles being formed (the right-hand sides) per unit time in the discharge volume.
The total power can be determined by summing over the volume and surface. The power balance equation has the following form:
N R F = n e n a σ a i V e C i d V d = 0.25 8 e T e π m e exp ( Δ U T e ) 5 2 T e S G D C + 0.43 n e 2 e T e M Δ U S G D C
where NRF is the HF power absorbed by plasma.

4. Development of the Ion Source Thermal Model

The ion beam current (I0) is formed by ions entering the holes of the IES screen electrode (SE). Specific power consumption (Ci) to produce ion current (“the ion cost”) is [44,45]: Ci = NRF/I0, where: NRF is the RF generator power.
The ion cost is: Ci = CidI0/Iu.
The major portion of the RF power input (80–90%) is converted into heat when plasma components—ions, electrons, and photons—deposit on the internal surfaces of the GDC and the surface of the IES screen electrode facing the GDC. The value of the specific power flux converted into heat is defined by the flow density for ions and electrons depositing on the GDC walls.
In the conditions with unchanging plasma parameters, the number of particles formed in plasma per unit time should be equal to the number of particles leaving plasma per unit time. If dielectric or metal walls under a floating potential are present, the conditions for particles to escape towards the wall must satisfy the condition of density equality for the electron and ion currents:
ji = je
Since the electron mobility substantially exceeds the ion mobility, a transition (near-wall) layer, in which the quasi-neutrality condition is violated, should be established near the GDC surface. Between the GDC wall and undisturbed plasma, a potential barrier for electrons (∆U) is established, thanks to which the electrons are retained in plasma. The ∆U value is driven by the condition of the current continuity (a balance between the generation of charged particles in plasma and their deposition on the GDC walls).
It follows from Equation (3) taking Equations (1) and (2) into account that:
e n e 8 e T e π m e exp ( Δ U T e ) = 0.43 e n e 2 e T e M
whence it follows that the near-wall potential drop in the vicinity of the GDC dielectric surface is:
Δ U G D C = T e ln ( 1 0.84 M π m e )
The potential drop at the IES SE surface is defined by the following equation:
Δ U S E = T e ln ( 1 σ S E 0.84 M π m e )
where σSE is transparency of the screen electrode (SE) (the ratio of the electrode total area to the area of its perforation).
For Xenon, the value of the near-wall potential drop in the vicinity of the GDC dielectric surface is: ∆UGDC = 5.8 Te, and the potential drop at the SE surface is: ∆USE = 5.1 Te.
The IES electrodes are exposed to heat fluxes from ions, electrons, and photons. During ion source operation, the charge-exchange ions formed in the volume of the outward ion flow, both in the inter-electrode gap and in the beam neutralization zone located beyond the output cross-section of the accelerating electrode (AE), deposit on the AE surface. The IES electrodes are exposed to thermal radiation coming from the internal surface of the GDC, heated up to temperatures of about 300 °C.
The resulting charge-exchange ions have thermal velocity and about-zero energy. They are accelerated in the electrostatic field, deposit on the negatively charged AE (usually minus 200–500 V), which leads to formation of a potential barrier for electrons from the neutralization zone. The charge-exchange ions accelerating in the AE electrostatic field bombard its surface with an energy that is a function of the potential at the point of their formation. If the acquired energy of the charge-exchange ion exceeds the sputtering threshold of the AE material (30–40 eV), its surface is sputtered. It is this mechanism of AE surface erosion that dictates, to a large extent, the IS lifetime [32,46,47]. Heat fluxes brought by the charge-exchange ions to the AE are determined by the charge-exchange ion current and by the average energy transferred by them to the AE surface atoms. Experimentally, these currents are determined as the current in the AE power-supply circuit (10–11% of the ion beam current).
As already noted, the electromagnetic field is the only supplier of thermal power to the IS. Electromagnetic energy is absorbed by plasma and carried over by ions and electrons to the GDC wall, being converted into heat, while a part of it is absorbed by the IS structural elements and dissipates as heat. Since the IS operates in vacuum, there is no convective heat transfer. In these conditions, it is sufficient to factor in conductive and radiative heat exchange in order to describe the IS thermal state (the source of heat dissipation is a gas-discharge plasma). The conductive heat exchange inside and between the IS structural elements is described by the Fourier law:
( λ T ) = Q i
where is the differential operator; λ is the thermal conductivity coefficient; Qi is a function of volumetric heat sources (heat dissipation in the IS caused by the induced eddy currents);and T is the temperature.
In the presented thermal model, for any structural element the thermal conductivity coefficient is assumed to be constant and independent of temperature; then the equation is transformed into the P.–S. Laplace equation:
Δ T = 0
The radiative heat transfers between the surfaces of the IS structure units, between the external surfaces of the IS structure units and the environment is described by the law of J. Stefan, L.E. Boltzmann:
Qrad = ε(Grad + GplσT4),
where ε is the emissivity of the material (it is equal to the absorption capacity over the entire wave length spectrum according to the gray body assumption [48]); Grad is the parameter that determines the input heat flux of the mutual radiative heat exchange between the units of the IS structure; Gpl is the parameter that determines the power flux from the discharge plasma; and σ is the Stefan–Boltzmann constant.
Generalization of thermal processes occurring in the IS is presented in [40].
Equations (4) and (6) fully describe the IS thermal state. Due to the overwhelming use of thin-walled elements in the IS structure, it is not necessary to take temperature distribution over the thickness of structural elements into account, i.e., the internal heat sources can be factored in as an integral part of Gpl, the more so since such sources are present in a thin skin layer during the inductive heating.
The main problem is the need to take contact resistances between different structural elements into account, which is determined as applicable, depending on the configuration under consideration.
The thermal state calculation procedure includes the following:
-
Calculate discharge plasma parameters in the GDC;
-
Determine heat fluxes from plasma over the surface of the IS structure;
-
Solve the problem of heat transfer from the IS structural elements using Equation (5) by numerical methods with account of the boundary conditions (environmental conditions), heat fluxes from the structural elements in view of contact resistances and input heat fluxes from plasma.
The IS design is transformed into a simplified model, which is covered by a computational grid to subsequently solve the problem by the method of discrete elements, in each of which the temperature, thermal conductivity, and degree of blackness (emissivity factor) are assumed to be the same. Accordingly, the accuracy of the performed calculation will be determined by the number of discrete, which in its turn depends on the computing system capacity available to the computing engineer. In our case of the IS thermo-physical state analysis, the structure was divided into more than 100,000 discrete elements (Figure 2), for each of which the following was assigned: spatial coordinates, the degree of blackness (emissivity factor), the thermal conductivity coefficient.
Radiation fluxes entering each element are presented in [40] and are determined by the following equations:
i = 1 M ( δ j i ε i F j i 1 ε i ε i ) 1 A i Q i R A D = i = 1 N ( δ j i F j i ) σ T i 4
where ε i is the degree of blackness (emissivity factor) for the surfaces of the elements; δ j i is the L. Kronecker symbol—the indicator of equality of the elements: a function of two integer variables, which is equal to one if the elements are equal, and equal to zero in any other case; A i are the areas of the elements; F j i is the matrix of the element mutual radiation coefficients calculated by:
F j i = 1 A i A i A j cos θ i cos θ j π r 2 d ( A j ) d ( A i )
where θ i are the angles between the normal to element i and the axis connecting the centers of elements i and j.
The heat balance equation for each partition element has the form:
Q i = ( Q j i R A D + Q j i C o n d ) + W i p l a s m a + W i n d
where Q i is the thermal power flux coming out of the i-th element of the computational grid; Q j i R A D is the radiative heat flux coming to the i-th element from the j-th element derived from Q r a d by (7); W i p l a s m a is the aggregate thermal power flux reaching leaving the surface of the element (in our case, these are the heat fluxes onto the GDC wall:
W G D C = j B S G D C ( Δ U G D C + V i + 2.5 T e ) ; onto the screen electrode surface facing the GDC: W S E = j B S K ( Δ U S E + V i ϕ e ) + 0.92 j B ( Δ U S E + V i ϕ e ) + I e S E ( 2.5 T e + ϕ e ) ; onto the AE surface facing the SE: W i A E = I i A E ( α E C P ϕ e A E + 0.5 V i ) I —in some cases, they turn to zero); W i n d is the thermal power flux from the inductor, which is dissipated on the active resistance of the inductor and is defined by eddy currents in the inductor (in some cases, they turn to zero); S G D C is the internal area of GDC surface; I i A E = 0.01 I 0 is the charge-exchange ion current; E C P is the average energy of the incident charge-exchange ions; α = 0.9 is the ion energy accommodation coefficient; ϕ e A E is the electronic work function of the AE material; S K is the area of the SE annular section; ϕ e is the electronic work function of the SE conductive material; Q j i C o n d are the thermal conductive fluxes between the neighboring elements i and j, which are defined in the following way:
-
If discrete elements belong to one and the same structural element—by using the finite-difference expression of the Fourier thermal conductivity law, where, instead of the temperature gradient, the temperature difference of adjacent elements attributed to the effective distance between the elements is used [40]:
Q j i C o n d = λ ¯ Δ T j i Δ x j i
-
If discrete elements belong to different structural elements—this heat flux is software-defined, and the effective value of the thermal conductivity coefficient is used as λ ¯ that is derived from λ i and λ j of the contacting surfaces [49] in the case of an ideal thermal contact; in the case of a real thermal contact, the heat flux is derived from the relation:
Q j i C o n d = 1 R C Δ T j i
where R c is the numerical value of the contact resistance.
The resulting equation defining the HF power balance has the following form:
W H F = I 0 ( Δ U S E + V i ) + I 0 γ σ S G D C S 0 ( Δ U G D C + V i + 2.5 T e ) + I 0 1 γ σ γ σ ( Δ U S E + V i ) + 2.5 I 0 γ σ T e + 0.2 W H F
Model calculations for the developed IS design have been conducted. Physical parameters of the materials, of which the IS structural elements are made, and initial temperatures of such elements (room temperature) were assumed as initial data. The results of the IS analysis are shown in the diagrams below (Figure 3, Figure 4 and Figure 5).
In situ measurements of IS temperature patterns using a thermal imager were performed to verify accuracy of the developed thermal model. As an example, the photo shows a IS GDC thermogram (Figure 6). It can be seen that the GDC is almost homogeneously heated. The temperature falls near the gas inlet. It is difficult to define temperature profile for the AE perforated part (Figure 7) due to the insufficient resolution of the thermal imaging camera. However, the conclusion can be made that the data obtained through the analytical model are correct.
The performed temperature measurements demonstrate consistency between the predicted and experimental temperature values for the IS surfaces with a high engineering accuracy, which makes it possible to use the proposed thermal model developed for predicting temperature distributions while making calculations for predicting heating of the IS structural elements.

5. Thermomechanical Model Development

Heating of the IS structural elements causes stress and strain (deformation) in these elements. Stability of the IES geometric parameters plays an especially important role in the IS operating capability. The IES unit is the most critical and complex IS element from the structural and technological stand point. More than 95% of power supplied to the IS is used in the IES [50]. To provide for the IES unit reliable operation, it is necessary to ensure the inter-electrode gap stability in various IS operation modes and the minimum possible deformation of IES electrodes. Operation of electrodes as part of the IES features significantly non-uniform heating and deformation that can cause deterioration in the source output characteristics (the ion current density and thrust), and also lead to a high-voltage breakdown [46,51,52].
One of the proven technical solutions to ensure stability of IES unit operation is the usage of profiled electrodes with an initially designed deflection (bend). Non-uniform heating leads to additional deflections of electrodes. Due to the above issues, it is important to perform, at the design engineering stage, a preliminary stress–strain electrode behavior analysis, as well as an analysis of their linear and angular displacements during the thermal deformation process.
The approach to calculate IES electrodes is based on application of the electrode continuous model representing it as a structurally orthotropic round plate, or a shallow spherical panel (Figure 8) [53]. The electrode element located at distance r from the electrode center, with initial thickness h, and final thickness h + dh, with depth dr and width dφ is assumed as a calculation unit. This element is decomposed into sublayers with thickness dz, in each of which the radial σr and tangential σφ stress components are given.
The structural orthotropy coefficients for such model are the cylindrical section and the meridian section filling factors, kr and kφ, respectively, in a thin-walled structure with close-packed perforation. The use of these dimensionless factors made it possible to utilize the thermo-elasticity equations [53] in combination with solutions developed to calculate perforated shells loaded with pressure, membrane forces, or bending moments [54]. The solutions were additionally adjusted by introducing additional non-uniform heating in the radial direction with simultaneous application of active loads, which allowed us to calculate nonlinear deformation for both flat and spherical electrodes [55]. Based on such a thermomechanical model and on mathematical methods of the boundary integral and integro-differential equations [56], a method to calculate thermal stability, post-critical (post-buckling) behavior and nonlinear bending for the flat IES electrodes has been developed [57].
The basic differential equations for axisymmetric thermo-elastic deformation and nonlinear bending of structurally orthotropic annular plates and of shallow spherical segments with radially variable stiffness, being under the non-uniform heating conditions along the radius and thickness, were obtained with the following assumptions:
-
The elasticity modulus and thermal linear expansion coefficient depend on temperature in the radial direction only;
-
The Kirchhoff—Love’s hypotheses are valid [58] since we consider thin electrodes;
-
A rectilinear shell element that is normal to the median surface remains to be rectilinear, normal to the deformed surface and maintains its length after deformation (the straight normal hypothesis);
-
Normal stresses in area elements parallel to the median shell surface are low compared to others and can be neglected;
-
Tangent stresses are static and drive equilibrium conditions;
-
The segments under consideration have an axisymmetric initial deflection.
With such an approach, calculation errors do not exceed δ/R, where δ is the shell thickness and R is the shell radius.
Stringent requirements are imposed on IES electrodes with respect to minimization of variations in the IES inter-electrode gaps during IS operation. The requirement relating to the permissible additional electrode deflections ∆w(0) (Figure 9) caused by the temperature fields during IS operation, is expressed as a fraction of the specified inter-electrode gap in the ES.
The electrodes are primarily impacted by the following:
-
Temperature fields in the IES, varying in the radial direction with a temperature differential of tens of degrees between the electrode center and the outer envelope;
-
Initially designed deflections of spherical electrodes w0(0) represented as spherical segments that can reach a height of 10–30 mm and more;
-
Electrode material linear thermal expansion coefficients;
-
Electrode structural orthotropy coefficients kr and kφ depending on the perforation geometry (they are constant over thickness since the electrodes have through holes).
With such initial data, thermo-elasticity equations for a structurally orthotropic modeling round-shape plate or a spherical segment that are made of an isotropic material, become as follows [50]:
ε ϕ z ( r , z ) = 1 E ( r ) [ σ ϕ ( r , z ) μ k r ( r ) σ r ( r , z ) ] + α ( r ) T ( r , z )
ε r z ( r , z ) = 1 E ( r ) [ σ r ( r , z ) μ k ϕ ( r ) σ ϕ ( r , z ) ] + α ( r ) T ( r , z )
where: ε r z ( r , z ) and ε ϕ z ( r , z ) are the relative deformations in radial and circumferential directions for the plate points located at distance z from the median surface; σ r ( r , z ) and σ ϕ ( r , z ) are the normal stresses in the plate acting along the edges of the elementary layer filled with material; E(r) is the elasticity modulus of the electrode material; α(r) is the electrode material thermal linear expansion coefficient; μ is the S.D. Poisson’s coefficient; T(r, z) is the change in the electrode temperature during its heating relative to the temperature of its assembly as a part of the IES.
The problem concerning an equilibrium of a structurally orthotropic shallow shell spatial element is reduced to the problem concerning an equilibrium of a median surface element, to which the following forces and moments are applied: internal forces (the transverse force Qr in the shallow spherical segment cylindrical section, the normal force Nr in the cylindrical section of the shallow spherical segment, the normal force Nφ in the meridian section, and the tangential stress τ r z in the cylindrical section of the shallow spherical segment), and the moments (the radial and circumferential bending moments Mr, Mφ); the positive directions of such forces and moments are shown in Figure 10.
The element equilibrium equations given the initial electrode deflection w 0 are written as follows [50]:
d d r ( r N r ) N ϕ = 0
d d r { r [ Q r N r d ( w + w 0 ) d r ] } = 0
d d r ( r M r ) M ϕ + r Q r = 0
The membrane forces and bending moments are defined by the following relations:
N r = k r [ ( ε r + μ k ϕ ε ϕ ) D N ( 1 + μ k ϕ ) N T ]
N ϕ = k ϕ [ ( ε ϕ + μ k r ε r ) D N ( 1 + μ k r ) N T ]
M r = k r [ ( d 2 w d r 2 + μ k ϕ r d w d r ) D M + ( 1 + μ k ϕ ) M T ]
M ϕ = k ϕ [ ( 1 r d w d r + μ k r d 2 w d r 2 ) D M + ( 1 + μ k r ) M T ]
D N ( r ) = E ( r ) h ( r ) 1 μ 2 k r ( r ) k ϕ ( r )
D M ( r ) = E ( r ) h 3 ( r ) 12 [ 1 μ 2 k r ( r ) k ϕ ( r ) ]
N T ( r ) = E ( r ) α ( r ) 1 μ 2 k r ( r ) k ϕ ( r ) h / 2 h / 2 T ( r , z ) d z
M T ( r ) = E ( r ) α ( r ) 1 μ 2 k r ( r ) k ϕ ( r ) h / 2 h / 2 T ( r , z ) z d z
where DN is the cylindrical tensile stiffness of the electrode structural-orthotropic model; DM is the cylindrical bending stiffness of the electrode structural-orthotropic model; NT is the membrane force caused by thermal action on the electrode; MT is the bending moment caused by thermal action on the electrode.
The resolving integral equations describing nonlinear bending for electrodes with radially variable stiffness and with non-uniform heating along the radius and thickness, are presented in [50], whereas the temperature pattern of electrodes is represented as T ( ρ ) = T ( 1 ) + Δ T r f T 0 ( ρ ) , where: T(1) is the temperature on the outer electrode envelope; f T 0 ( ρ ) is the law of temperature drop changes along the electrode radius.
The peculiarities of IES electrode loading and of their fixing are specified in boundary conditions, for example, for the case of loose fixing we have: N r ( ρ a ) = N a , N r ( 1 ) = N b , M r ( ρ a ) = M a , M r ( 1 ) = M b , where N a and N b is the intensity of radial compressive forces that are uniformly distributed on the inner and outer envelopes; ρ a = a b , a is the inner envelope radius; b is the outer envelope radius; M a and M b are the intensities of the radial bending moments on the inner and outer envelopes, respectively.
Calculations were carried out for a carbon–carbon composite high-density material based on the “IPRESSKON” nonwoven frame. The IS IES electrodes additional deflection values under various temperature loading for the case of loose IES electrode fixing are shown in the graphs (Figure 11 and Figure 12) and in Table 1.
It is obvious that the used material ensures sufficiently low values of additional electrode deflections that occurring the IES structure even with significant temperature gradients along the electrode radii, which are within 1.5 mm for a temperature gradient of 80 °C. According to the above calculations for temperature fields that will occur in the IS structure during its operation in the space flight conditions, the temperature gradient will not exceed 10 °C for the AE and 81 °C for the SE, which corresponds to the values of additional deflections at the level of 0.005 mm for the AE, and at the level of 1.5 mm for the SE.

6. The Study of Electrode Thermal Deformation

The ion-extraction system is the most critical and complicated part of an ion thruster/source in terms of its manufacturing and design. This is stipulated by the fact that more than 95% of supplied power is consumed by the ES, as well as by the challenging task to maintain a gap between the screen and acceleration electrodes. The electrodes operate under uneven thermal loading. Heating of electrodes leads to their deformation causing deterioration in IES performance and probable occurrence of a breakdown in the inter-electrode gap. Therefore, the gap should vary within a rigorously limited range.
The decision was made to assess the level of electrode deformation during IS operation. Similarly to the described above, this issue can be studied in two ways:
-
By using a well-known powerful resource, such as, for instance, the ANSYS computing system;
-
By using the developed simplified algorithm to evaluate the magnitude of additional deflection;
-
By determining the manufactured electrode deformation experimentally.
That is, it is necessary to determine additional deflection of the electrode, which occurs under uneven thermal heating, as a function of its initial shape and kind of material, and to define correlating between the magnitudes of electrode deformation obtained by the developed algorithm and by the numerical simulation in the ANSYS engineering simulation software with account of the data obtained through the electrode thermal deformation experimental study.
The IES perforated electrodes were selected for the sake of this experimental study. A numerical simulation of electrode deflection was conducted by a simplified algorithm with the use of the following basic design data: for SE − T = 330°; λ3kp = 35.5 (Figure 13), for AE − T = 272°; λ3kp = 36.4 (Figure 14). The graphs illustrate a possible choice of the initial electrode deflection required to ensure additional deflection (the additional deflection as a function of the initially designed one).
When performing numerical simulations in the ANSYS software environment, the number of holes was reduced while maintaining structural orthotropy to calculate loads and deformations (Figure 15 and Figure 16). Figure 17 shows the numerical simulation of the electrode thermal fields.
Figure 18 shows deflection as a function of the temperature differential along the electrode radius for a numerical simulation by a simplified algorithm and in the ANSYS software environment.
A high degree of consistency between the results of a numerical simulation and those obtained by simplified algorithm (over 90%) has been demonstrated with a significant reduction to the calculation time.
However, in order to assess experimental data on deformation of the manufactured AE as a part of the high-frequency IS structure, thermal fields of the IS structure as a whole were numerically simulated.
Figure 19 shows the IS thermal fields with the 350 W power supplied to the inductor.
Two pins were installed and their displacements were assessed in order to assess thermal deformation of the electrodes. The temperature of the electrodes was measured by thermocouples in the center and at the periphery. The AE temperature was additionally measured with a thermal imager. The following values for AE deformation were obtained with the power of 350 W supplied to the inductor:
-
According to experimental data, AE deformation for a spherical segment was 0.15 mm;
-
The deformation numerical simulation with loose AE fixing for a spherical segment without taking the cone height into account resulted in 0.13 mm;
-
The total deformation of AE with loose fixing was 0.16 mm;
-
The total deformation of AE with elastic fixing was 0.11 mm.
A high degree of consistency between the numerical modeling results both when using a simplified algorithm and the ANSYS software environment, and the experimental results for the AE deformation has been demonstrated.

7. Results and Discussion

A thermal model has been developed that links specific heat fluxes on the surfaces of GDC, SE and AE to the IS integral parameters that can be experimentally measured with high accuracy. It is shown that the ion beam current (current in the SE circuit) should be taken as the thermal state governing parameter, through the value of which the plasma parameters and specific power fluxes to the GDC walls and to the SE are defined. It is possible to define numerical values of the electron temperature more accurately, in order to derive the power balance on the GDC walls and on the SE surface, and the plasma potential, on the one hand, from the overall heat dissipation in the GDC volume compared with the incident RF power of the RF generator taking the efficiency of power input into plasma in to account, and, on the other hand, from the experiment concerning the beam ion energy spread. The model allows for its gradual complication by increasing the IS level of detail, which brings the simplified IS representation closer to the real one. The expected heating of IS surfaces can reach 300–400 °C, which is critical in terms of performance of structural materials. For this reason, it is impossible to use approximation of the smallness of the working surface temperature changes that allows the use of linear approximation and accepts analytical solutions to the thermal problem. The thermal model factors in the heat exchange between the structural units, either by radiation, or by thermal conductance. In the latter case, the model allows to account for thermal contact resistances in the interface surfaces.
The IS integral output characteristics (Table 2) providing the minimum ion beam divergence angle were determined by tests.

8. Conclusions

Thus, it is shown that the IS developed and manufactured on the basis of the conducted modeling results operates steadily. Its output characteristics correspond to the calculated values with high accuracy, and, therefore, the physical and mathematical IS models developed in the course of the project implementation can be extensively used going forward in the development of similar technical devices with a different standard size.
Based on numerical modeling results, and upon obtaining a high level of agreement with experimental data, we can conclude the following:
-
Correlation between the results obtained through numerical modeling using the ANSYS software system and by the developed algorithm was defined;
-
The additional deflection value calculated by the proposed algorithm agrees with the results obtained through the numerical simulation in ANSYS and with the experimental data within 10%;
-
The conducted study does not allow us to precisely assess the accuracy of the numerical modeling results, but a fairly good agreement between the calculation results obtained by two independent methods and the experimental data provides more confidence in the reliability of the developed model.

Author Contributions

Conceptualization, methodology, validation, writing—review and editing, supervision, project administration, funding acquisition V.V.S., software, A.I.M., A.Y.K.; formal analysis, investigation, resources, data curation, writing—original draft preparation, visualization, V.V.S., A.I.M., A.Y.K. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the Government of Russian Federation (Grant number 075-15-2019-1894).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Acknowledgments

This research work was conducted in the Laboratory Space Flight Mechanics of the Moscow Aviation Institute.

Conflicts of Interest

The authors declare no conflict of interests.

Abbreviations

AEaccelerating electrode
GDCgas discharge chamber
HFhigh-frequency
IESion-extraction system
ISion source
RITradio-frequency ion thruster
SCspacecraft
SDspace debris
SSCservice spacecraft
SEscreen electrode

Nomenclature

Indexes
AEindex of accelerating electrode parameter
EPindex of compensating propulsion system parameter
GDCindex of gas discharge chamber parameter
ISindex of ion source parameter
SCindex of spacecraft parameter
SDindex of space debris parameter
SEindex of screen electrode parameter
Parameters
BHFinduction of high-frequency magnetic field
CIspecific power consumption
CIDeffective energy spent for the production of one ion in the discharge
CRtotal radiation loss
ddistance between service spacecraft and space debris
eelectron charge
ECPaverage energy of the incident charge-exchange ions
EINDintensity of high-frequency electric field
Fjimatrix of the element mutual radiation coefficient
Fdforce acting on space debris
FEPforce acting on the service spacecraft from compensating propulsion system
Gradparameter that determines the input heat flux of the mutual radiative heat exchange between the units of the ion source structure
Gplparameter that determines the power flux from the discharge plasma
Ieelectron current
Iiion current
IiAEcharge-exchange ion current
I0ion beam current
kratio between force acting on space debris and force acting on the service spacecraft from compensating propulsion system
Mpropellant gas atom (ion) mass
MSCmass of service spacecraft
MSDmass of space debris
neelectron density
niion density
NRFhigh-frequency generator power
Qifunction of volumetric heat source
Qradradiative heat
RCnumerical value of the contact resistance
SGDCinternal area of gas discharge chamber surface
SKarea of the screen electrode annular section
Ttemperature
Taatom temperature
Teelectron temperature
Tiion temperature
Veelectron velocity
Vdgas discharge chamber plasma volume
Vifirst potential of the propellant gas atomionization
WEEthermal power flux onto screen electrode surface facing the gas discharge chamber
WiAEthermal power flux onto accelerating electrode surface facing the screen electrode
WGDCthermal power flux onto the gas discharge chamber wall
Windthermal power flux from the inductor
differential operator
αion energy accommodation coefficient
δL. Kronecker symbol
εemissivity of the material
φeelectronic work function of the screen electrode conductive material
φeAEelectronic work function of the accelerating electrode material
σStefan-Boltzmann constant
σSEtransparency of the screen electrode
λthermal conductivity coefficient
θIangle between the normal to element I and the axis connecting the centers of elements i and j
ωhigh-frequency generator frequency
Unear-wall plasma potential drop
UGDCplasma potential drop near gas discharge chamber surface
USEplasma potential drop near screen electrode surface

References

  1. Kaplan, M.H. Space Debris Realities and Removal. SOSTC Improving Space Operations Workshop Spacecraft Collision Avoidance and Co–Location; The Johns Hopkins University, Applied Physics Laboratory: Baltimore, MD, USA, 2010. [Google Scholar]
  2. Liou, J.C.; Johnson, N.L.; Hill, N.M. Controlling the Growth of Future LEO Debris Population with Active Debris Removal. Acta Astronaut. 2010, 66, 648–653. [Google Scholar] [CrossRef]
  3. Debus, T.J.; Dougherty, S.P. Overview and Performance of the Front-end Robotics Enabling Near-Term Demonstration (FREND) Robotic Arm. In Proceedings of the AIAA Infotech@Aerospace Conference, Seattle, WA, USA, 6–9 April 2009. [Google Scholar]
  4. Kaiser, C.; Sjoberg, F.; Delcura, J.M.; Eilertsen, B. SMART-OLEV—An Orbital Life Extension Vehicle for Servicing Commercial Spacecrafts in GEO. Acta Astronaut. 2008, 63, 400–410. [Google Scholar] [CrossRef]
  5. Melamed, N.; Griffice, C.P.; Chobotov, V. Survey of GEO Debris Removal Concepts. In Proceedings of the 59th International Astronautical Congress, Glasgow, UK, 29 September–3 October 2008. [Google Scholar]
  6. Bosse, A.B.; Barnds, W.J.; Brown, M.A.; Creamer, N.G.; Feerst, A.; Hensgaw, C.G.; Hope, A.S.; Kelm, B.E.; Klein, P.A.; Pipitone, F.; et al. SUMO: Spacecraft for the Universal Modification of Orbits. In Proceedings of the International Symposium on Defense and Security, Orlando, FL, USA, 12–16 April 2004. [Google Scholar]
  7. Norman, T.L.; Gaskin, D. Conceptual Design of an Orbital Debris Defense System; NASA-CR-197211; NASA: Morgantown, WV, USA, 1994.
  8. Murdoch, N.; Izzo, D.; Bombardelli, C.; Carnelli, I.; Hilgers, A.; Rodgers, D. Electrostatic Tractor for Near Earth Object Deflection. In Proceedings of the 59th International Astronautical Congress, Glasgow, UK, 29 September–3 October 2008. [Google Scholar]
  9. Bondarenko, S.G.; Lyagushin, S.F.; Shifrin, G.A. Prospects of Using Lasers and Military Space Technology for Space Debris Removal. In Proceedings of the 2nd European Conference on Space Debris, Darmstadt, Germany, 17–19 March 1997. [Google Scholar]
  10. Nishida, S.-I.; Kawamoto, S.; Okawa, Y.; Terui, F. Space debris removal system using a small satellite. Acta Astronaut. 2009, 65, 95–102. [Google Scholar] [CrossRef]
  11. Kawamoto, S.; Nishida, S.; Kibe, S. Research on a Space Debris Removal System; NAL Research Progress: Tokyo, Japan, 2003. [Google Scholar]
  12. Starke, J.; Bischof, B.; Foth, W.-O.; Gunther, J.J. Roger a Potential Orbital Space Debris Removal System. In Proceedings of the 38th COSPAR Scientific Assembly, Bremen, Germany, 15–18 July 2010. [Google Scholar]
  13. Absalamov, S.K.; Andreev, V.B.; Colbert, T.; Day, M.; Egorov, V.V.; Gnizdor, R.U.; Kaufman, H.; Kim, V.; Koryakin, A.I.; Kozubsky, K.N.; et al. Measurements of Plasma Parameters in the SPT-100 Plume and its effect on Spacecraft Components. In Proceedings of the 28 Joint Propulsion Conference and Exhibit, Nashville, TN, USA, 6–8 July 1992. [Google Scholar]
  14. Raikunov, G. (Ed.) Prevention of Space Debris Formation/Under the Science; Ed. Doct. Tech. Sciences; FIZMATLIT: Moscow, Russia, 2014; (In Russian). ISBN 9785922115049. [Google Scholar]
  15. Cherkasova, M.; Mogulkin, A.; Nigmatzyanov, V.; Peysakhovich, O.; Popov, G. Computational and experimental investigation of a weakly diverging ion beam source. In IOP Conference Series: Materials Science and Engineering, Proceedings of the 18th International Conference “Aviation and Cosmonautics” (AviaSpace-2019), Moscow, Russia, 18–22 November 2019; IOP: Bristol, UK, 2020; Volume 868, p. 012019. [Google Scholar]
  16. Obukhov, V.A.; Mogulkin, A.I.; Peysakhovich, O.D.; Nigmatzyanov, V.V.; Svotina, V.V. Modeling and experimental study of an ion source with a weakly diverging ion beam. In IOP Conference Series: Materials Science and Engineering, Proceedings of the II International Conference «MIP: Engineering-2020: Modernization, Innovations, Progress: Advanced Technologies in Material Science, Mechanical and Automation Engineering», Krasnoyarsk, Russia, 16–18 April 2020; IOP: Bristol, UK, 2020; Volume 862, p. 062082. [Google Scholar]
  17. Brown, I.; Lane, J.; Youngquist, R. A Lunar-Based Spacecraft Propulsion Concept—The Ion Beam Sail. Acta Astronaut. 2007, 60, 834–845. [Google Scholar] [CrossRef] [Green Version]
  18. Kitamura, S. Large Space Debris Reorbiter Using Ion Beam Irradiation. In Proceedings of the 61th International Astronautical Congress, Prague, Czech Republic, 27 September–1 October 2010. [Google Scholar]
  19. Bombardelli, C.; Peláez, J. Ion Beam Shepherd for Contactless Space Debris Removal. J. Guid. Control Dyn. 2011, 34, 916–920. [Google Scholar] [CrossRef] [Green Version]
  20. Bombardelli, C.; Pelaez, J. Universidad Politecnica de Madrid. Sistema de Modificación de la Posicióny Actitud de Cuerpos en Órbita por Medio de Satélites Guía. ES Patent 2365394B2, 3 October 2011. [Google Scholar]
  21. Edwards, C.H.; Wallace, N.C. The T5 Ion Propulsion Assembly for Drag Compensation on GOCE. In Proceedings of the 2nd International GOCE User Workshop GOCE, The Geoid and Oceanography, Frascati, Italy, 8–10 March 2004. [Google Scholar]
  22. Löb, H.W. Ein elektrostatisches Raketentriebwerk mit Hochfrequenzionenquelle. Astronaut. Acta VIII 1962, 1, 49. [Google Scholar]
  23. Loeb, H.W. Principle of Radio-Frequency Ion Thrusters RIT. In Workshop RIT-22 Demonstrator Test at Astrium at Giessen University; Giessen University: Giessen, Germany, 2010. [Google Scholar]
  24. Groh, K.H.; Loeb, H.W. State-of-the-art of radio-frequency ion thrusters. J. Propuls. Power 1991, 7, 573–579. [Google Scholar] [CrossRef]
  25. Killinger, R.; Kukies, R.; Surauer, M.; Tomasetto, A.; van Holtz, L. ARTEMIS Orbit Raising In-flight Experience with Ion Propulsion. Acta Astronaut. 2003, 53, 607–621. [Google Scholar] [CrossRef]
  26. Killinger, R.; Leiter, H.; Kukies, R. RITA Ion Propulsion Systems for Commercial and Scientific Applications. In Proceedings of the 43rd AIAA/ASME/SAE/ASEE Joint PropulsionConference & Exhibit, Cincinnati, OH, USA, 8–11 July 2007. [Google Scholar]
  27. Muller, H.; Kukies, R.; Bassner, H. EMC Tests on the RITA Ion Propulsion Assembly for the ARTEMIS Satellite. In Proceedings of the 28th AIAA, SAE, ASME, ASEE Joint Propulsion Conference and Exhibit, Nashville, NT, USA, 6–8 July 1992. [Google Scholar]
  28. Angelo, N.G.; Clark, S.; Wallace, N.; Collingwood, C.; Guarducci, F. Qualification of the T6 Ion Thruster for the BepiColombo Mission to the Planet Mercury. In Proceedings of the 32nd International Electric Propulsion Conference, Wiesbaden, Germany, 11–15 September 2011. [Google Scholar]
  29. Loeb, H.W.; Feili, D.; Popov, G.A.; Obukhov, V.A.; Balashov, V.V.; Mogulkin, A.I.; Murashko, V.M.; Nesterenko, A.N.; Khartov, S. Design of High-Power High-Specific Impulse RF-Ion Thruster. In Proceedings of the International Electric Propulsion Conference, Wiesbaden, Germany, 11–15 September 2011. [Google Scholar]
  30. Riaby, V.A.; Masherov, P.E.; Obukhov, V.A.; Savinov, V.P. Effectiveness of Radiofrequency Inductively Coupled Plasma Sources for Space Propulsion. High Volt. Eng. 2013, 39, 2077–2088. [Google Scholar]
  31. Abgaryan, V.K.; Akhmetzhanov, R.V.; Loeb, H.W.; Obukhov, V.A.; Cherkasova, M.V. High-Frequency Ion Sources of Inert and Chemically Active Gases. J. Surf. Investig. 2012, 6, 693–698. [Google Scholar] [CrossRef]
  32. Abgaryan, V.K.; Akhmetzhanov, R.V.; Loeb, H.W.; Obukhov, V.A.; Cherkasova, M.V. Numerical Simulation of a High-Perveance Ion-Extraction System with a Plasma Emitter. J. Surf. Investig. 2013, 7, 1092–1099. [Google Scholar] [CrossRef]
  33. Bombardelli, C.; Merino-Martinez, M.; Galilea, E.A.; Pelaez, J.; Urrutxua, H.; Herrera-Montojo, J.; Iturri-Torrea, A. Ariadna Call for Ideas: Active Removal of Space Debris. Ion Beam Shepherd for Contactless Debris Removal; Advanced Concept Team; Ariadna ID: 10-4611; ESTEC-ESA: Madrid, Spain, 2011. [Google Scholar]
  34. Godyak, V.A.; Piejak, R.B.; Alexandrovich, B.M. Experimental Setup and Electrical Characteristics of an Inductively Coupled Plasma. J. Appl. Phys. 1999, 85, 703–712. [Google Scholar] [CrossRef]
  35. Kozhevnikov, V.V.; Khartov, S.A. Investigation of Local Plasma Parameters in the Discharge Chamber of a Low-Power High-Frequency Ion Thruster. Bull. Russ. Acad. Sci. Energy 2017, 3, 13–20. (In Russian) [Google Scholar]
  36. Kozhevnikov, V.V.; Khartov, S.A. Investigation of Local Plasma Parameters in the Discharge Chamber of a Low-Power High-Frequency Ion Engine with Multi-Electrode Probes. Bull. Russ. Acad. Sci. Energy 2016, 2, 26–33. (In Russian) [Google Scholar]
  37. Ryabyi, V.A.; Masherov, P.E. Integral and Local Diagnostics of a Model of an Energy-Efficient RF Source of an Ion Beam. Bull. Russ. Acad. Sci. Energy 2016, 2, 46–57. (In Russian) [Google Scholar]
  38. Riaby, V.A.; Obukhov, V.A.; Masherov, P.E. On the Objectivity of Plasma Diagnostics Using Langmuir Probes. High Volt. Eng. 2012, 38, 790–793. [Google Scholar]
  39. Obukhov, V.A.; Grigoryan, V.G.; Latyshev, L.A. Sources of Heavy Ions. In Plasma Accelerators and Ion Injectors; Science: Moscow, Russia, 1984. (In Russian) [Google Scholar]
  40. Kruglov, K.I. Modeling of Thermophysical Processes in a High-Frequency Ion Thruster. Ph.D. Thesis, Moscow Aviation Institute, Moscow, Russia, 2017. (In Russian). [Google Scholar]
  41. Abgaryan, V.K.; Abgaryan, M.V.; Nadiradze, A.B.; Nigmatzyanov, V.V.; Semenov, A.A. Radiative Heat Flux from Excited Atoms in a Radiofrequency Ion Thruster. Tech. Phys. Lett. 2021, 47, 201–204. [Google Scholar] [CrossRef]
  42. Bugrov, G.E.; Kralkina, E.A. Calculation of Specific Energy Consumption for Radiation in a Low-Pressure Xenon Discharge. In Proceedings of the VI All-Union Conference on Plasma Accelerators and Ion Injectors, Dnepropetrovsk, Russia, 26–28 September 1986. (In Russian). [Google Scholar]
  43. Lotz, B. Plasma Physical and Material Physical Aspects of the Application of Atmospheric Gases as a Propellant for Ion-Thruster of the RIT-Type. Ph.D. Thesis, Justus-Liebig Universität, Bunkissen, Germany, 2013. [Google Scholar]
  44. Abgaryan, V.K.; Kruglov, K.I.; Mogulkin, A.I.; Obukhov, V.A.; Riaby, V.A. On the Cost of Ion Production in Radio-Frequency Ion Thrusters and Technological Ion Sources. J. Surf. Investig. 2017, 11, 1239–1245. [Google Scholar] [CrossRef]
  45. Abgaryan, V.K.; Riaby, V.A.; Yamashev, G.G. Power Loss to Generating Ion Current in Radio-Frequency Ion Thrusters and Technological Sources of Ions. J. Surf. Investig. 2017, 11, 1008–1013. [Google Scholar] [CrossRef]
  46. Antipov, E.A.; Balashov, V.V.; Weber, A.V.; Kuftyrev, R.Y.; Mogulkin, A.I.; Nigmatzyanov, V.V.; Pankov, A.I.; Popov, G.A.; Sitnikov, S.A.; Khartov, S.A. The Choice of Structural Materials for High-Frequency Ion Thruster. Electron. J. Proc. MAI 2013, 65. (In Russian) [Google Scholar]
  47. Abgaryan, V.K.; Akhmetzhanov, R.V.; Loeb, H.W.; Obukhov, V.A.; Cherkasova, M.V. Numerical Modeling of the Primary Ion Beam and Secondary Ion Flux in the Ion-Optical System of the Ion Thruster. Electron. J. Proc. MAI 2013, 71. (In Russian) [Google Scholar]
  48. Avduevsky, V.S.; Galitseysky, B.M.; Glebovetal, G.A. Basics of Heat Transfer in Aviation and Rocket and Space Technologies; Mashinostroenie: Moscow, Russia, 1992. (In Russian) [Google Scholar]
  49. Gartner, W.; Lotz, B.; Meyer, B.K. 3D Thermal Simulation of a μN-RIT. In Proceedings of the 33d International Electric Propulsion Conference, The George Washington University, Washington, DC, USA, 6–10 October 2013. [Google Scholar]
  50. Mogulkin, A.I. Mechanical–Mathematical Model of Deformations of Profiled Electrodes of Ion Thrusters. Ph.D. Thesis, Moscow Aviation Institute, Moscow, Russia, 2015. [Google Scholar]
  51. Balashov, V.V.; Mogulkin, A.I.; Nigmatzyanov, V.V.; Obukhov, V.A.; Smirnova, M.E.; Khartov, S.A. Ion Source. Russia Patent RU116273U1, 2 December 2011. (In Russian). [Google Scholar]
  52. Balashov, V.V.; Popov, G.A.; Antipov, E.A.; Ionov, A.V.; Mogulkin, A.I. Method of Manufacturing Electrodes of the Ion-Optical System. Russia Patent RU2543063C1, 9 October 2013. (In Russian). [Google Scholar]
  53. Birger, I.A. Round Plates and Shells of Rotation; Oborongiz: Moscow, Russia, 1961. (In Russian) [Google Scholar]
  54. Grigolyuk, E.I.; Filshtinsky, L.A. Perforated Plates and Shells; Nauka: Moscow, Russia, 1970. (In Russian) [Google Scholar]
  55. Gurov, A.F.; Fedorov, V.A. Investigation of Stability and Bending of Electrodes of Ion-Optical Systems. In Proceedings of the 3rd Interuniversity Conference on the Research and Design of Medium-Power Electric Propulsion Systems, Kharkov, Ukraine, 15–20 June 1975. (In Russian). [Google Scholar]
  56. Birger, I.A. Some Mathematical Methods for Solving Engineering Problems; Oborongiz: Moscow, Russia, 1956. (In Russian) [Google Scholar]
  57. Fedorov, V.A. Investigation of Thermal Stability, Supercritical Behavior and Nonlinear Bending of Densely Perforated Annular Plates of Variable Stiffness. Izv. SSSR Rigid Body Mech. 1976, 5, 197. (In Russian) [Google Scholar]
  58. Besukhov, N.N.; Goldenblat, I.I.; Badjanov, V.L. Calculations for Strength, Stability and Fluctuations in High Temperature Conditions; Mechanical Engineering: Moscow, Russia, 1965. (In Russian) [Google Scholar]
Figure 1. Forces arising in the “SSC-SD” virtual assembly (FIS is the force acting on the SSC (3) from the IS; FEP is the force acting on the SSC from the compensating propulsion system (1,2); Fd is the force acting on space debris).
Figure 1. Forces arising in the “SSC-SD” virtual assembly (FIS is the force acting on the SSC (3) from the IS; FEP is the force acting on the SSC from the compensating propulsion system (1,2); Fd is the force acting on space debris).
Aerospace 08 00189 g001
Figure 2. IS discrete elements.
Figure 2. IS discrete elements.
Aerospace 08 00189 g002
Figure 3. Predicted temperature patterns in the IS from the Slit-Type IES end.
Figure 3. Predicted temperature patterns in the IS from the Slit-Type IES end.
Aerospace 08 00189 g003
Figure 4. Predicted temperature patterns in the IS from the GDC end (sectional view).
Figure 4. Predicted temperature patterns in the IS from the GDC end (sectional view).
Aerospace 08 00189 g004
Figure 5. Predicted temperature patterns in the IS for the screen electrode.
Figure 5. Predicted temperature patterns in the IS for the screen electrode.
Aerospace 08 00189 g005
Figure 6. IS GDC thermogram.
Figure 6. IS GDC thermogram.
Aerospace 08 00189 g006
Figure 7. IS acceleration electrode thermogram.
Figure 7. IS acceleration electrode thermogram.
Aerospace 08 00189 g007
Figure 8. Plate element.
Figure 8. Plate element.
Aerospace 08 00189 g008
Figure 9. Deflection diagram.
Figure 9. Deflection diagram.
Aerospace 08 00189 g009
Figure 10. Internal forces and moments.
Figure 10. Internal forces and moments.
Aerospace 08 00189 g010
Figure 11. IS screen electrode additional deflection (∆w(0)) as a function of initially designed deflection (w0(0)) at different levels of temperature drop along SE radius.
Figure 11. IS screen electrode additional deflection (∆w(0)) as a function of initially designed deflection (w0(0)) at different levels of temperature drop along SE radius.
Aerospace 08 00189 g011
Figure 12. IS acceleration electrode additional deflection (∆w(0)) as a function of initially designed deflection (w0(0)) at different levels of temperature drop along AE radius.
Figure 12. IS acceleration electrode additional deflection (∆w(0)) as a function of initially designed deflection (w0(0)) at different levels of temperature drop along AE radius.
Aerospace 08 00189 g012
Figure 13. Additional deflection (∆w(0)) for the SE made of titanium alloy with 160 mm diameter as a function of its initially designed deflection (w0(0)).
Figure 13. Additional deflection (∆w(0)) for the SE made of titanium alloy with 160 mm diameter as a function of its initially designed deflection (w0(0)).
Aerospace 08 00189 g013
Figure 14. Additional deflection (∆w(0)) for the AE made of titanium alloy with 160 mm diameter as a function of its initially designed deflection (w0(0)).
Figure 14. Additional deflection (∆w(0)) for the AE made of titanium alloy with 160 mm diameter as a function of its initially designed deflection (w0(0)).
Aerospace 08 00189 g014
Figure 15. Numerical simulations of electrode geometry in the ANSYS software.
Figure 15. Numerical simulations of electrode geometry in the ANSYS software.
Aerospace 08 00189 g015
Figure 16. Element-by-element partition of electrodes in the ANSYS software—mesh model development.
Figure 16. Element-by-element partition of electrodes in the ANSYS software—mesh model development.
Aerospace 08 00189 g016
Figure 17. Thermal patterns of electrodes in the ANSYS software.
Figure 17. Thermal patterns of electrodes in the ANSYS software.
Aerospace 08 00189 g017
Figure 18. Electrode deflection w(0) as a function of temperature differential along the electrode radius calculated by simplified algorithm and in the ANSYS software.
Figure 18. Electrode deflection w(0) as a function of temperature differential along the electrode radius calculated by simplified algorithm and in the ANSYS software.
Aerospace 08 00189 g018
Figure 19. IS thermal fields with 350 W power supplied to the inductor.
Figure 19. IS thermal fields with 350 W power supplied to the inductor.
Aerospace 08 00189 g019
Table 1. Electrode design deflection simulation results as a function of initially designed deflection and temperature gradient along electrode radius.
Table 1. Electrode design deflection simulation results as a function of initially designed deflection and temperature gradient along electrode radius.
ΔTr, °C3040506080
w0(0), mmΔw(0), mm
Screen electrode
50.00710.00940.01180.01410.0188
7.50.00470.00630.00790.00940.0126
9.50.00370.004980.00620.00750.0099
12.50.00280.00380.00470.00570.0076
150.00240.00320.00390.00470.0063
200.00180.00240.0030.00360.0047
250.00140.00190.00240.00280.0038
300.00120.00160.0020.00240.0032
Acceleration electrode
50.00570.00760.00950.01140.0152
7.50.00430.00570.00710.00850.0114
9.50.00350.00470.00580.00700.0093
12.50.00270.00360.00460.00550.0073
150.00230.00310.00380.00460.0061
200.00170.00230.00290.00350.0047
250.00140.00190.00230.00280.0038
300.00120.00160.00200.00240.0031
Table 2. Experimental results for the IS with the slit-type perforation.
Table 2. Experimental results for the IS with the slit-type perforation.
Mass Flow Rate, mg/sPRFG, WUSE, VISE, AUAE, VIAE, AIsp, sIp, AP, mNW, keV
0.77122025000.1902400.00820830.1816.02.50
0.97631024000.4002200.00133960.433.12.40
0.97633024000.4102200.00134800.4133.92.40
0.48810030000.1003000.00818960.099.283.00
0.78115030000.1503000.00917780.1413.93.00
0.78116030000.1703000.00920150.1615.83.00
0.78118030000.1903000.00922520.1817.63.00
0.91722030000.2403000.01024210.2322.33.00
0.94723530000.2603000.01025420.2524.13.00
0.94723530000.2803000.01027370.2725.973.00
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Svotina, V.V.; Mogulkin, A.I.; Kupreeva, A.Y. Ion Source—Thermal and Thermomechanical Simulation. Aerospace 2021, 8, 189. https://0-doi-org.brum.beds.ac.uk/10.3390/aerospace8070189

AMA Style

Svotina VV, Mogulkin AI, Kupreeva AY. Ion Source—Thermal and Thermomechanical Simulation. Aerospace. 2021; 8(7):189. https://0-doi-org.brum.beds.ac.uk/10.3390/aerospace8070189

Chicago/Turabian Style

Svotina, Victoria V., Andrey I. Mogulkin, and Alexandra Y. Kupreeva. 2021. "Ion Source—Thermal and Thermomechanical Simulation" Aerospace 8, no. 7: 189. https://0-doi-org.brum.beds.ac.uk/10.3390/aerospace8070189

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop