• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    The far-field plasma characteristics of LHT40 low-power Hall thruster for commercial aerospace applications

    2021-10-31 08:15:00XinweiCHEN陳新偉ZuoGU顧左JunGAO高俊
    Plasma Science and Technology 2021年10期
    關(guān)鍵詞:劉超

    Xinwei CHEN (陳新偉), Zuo GU (顧左), Jun GAO (高俊),

    Shangmin WANG (王尚民), Tao CHEN (陳燾), Ning GUO (郭寧),Sanxiang YANG (楊三祥) and Chao LIU (劉超)

    Science and Technology on Vacuum Technology and Physics Laboratory, Lanzhou Institute of Physics,Lanzhou 730000, People’s Republic of China

    Abstract To achieve a better insight into the far-field plasma spatial distribution and evolution characteristics of the 300 W class low-power Hall thruster (LHT) for commercial aerospace applications, a dedicated and integrated plasma diagnostic system composed of seventeen Faraday cups (FC)and two triple Langmuir probes (TLP)is established to investigate the timeaveraged in situ spatial distribution characteristics of far-field ions and electrons.The ion current density (ICD), plasma potential, plasma density, and electron temperature at 1000 mm downstream of 300 W class LHT for commercial aerospace applications in the azimuthal angle range of ?90° to 90° were investigated under the conditions of different anode mass flow rates and discharge voltages.The results demonstrated that ICD, beam divergence angle, and mass utilization efficiency increased with increasing anode mass rate.The double-wings phenomenon was observed in the spatial distribution of ICD at large angles from the thruster axis, which is attributed to charge exchange collisions at increasing vacuum backpressure.The plasma electron temperature, electron density, and plasma potential parameters derived from the TLP decreased rapidly in the angle range from 0° to 30° and did not exhibit significant variations above 30°,which was also in good agreement with the results of the measured divergence angle of the FC.The discrepancy of average ion speed was calculated.The maximum error is better than 31.5%which checks the consistency between the TLP’s results and that of FC to some extent.

    Keywords: low power Hall thruster, plasma characteristics, commercial aerospace, plasma diagnostic system

    1.Introduction

    It is generally known that electric propulsion has the advantages of high specific impulse compared with chemical propulsion, which can significantly save the mass of propellant and reduce the launch cost [1].Among several electric propulsion systems, Hall thrusters (HT) have many advantages such as simple structure, low thrust, high specific impulse,long life, high integration, and the high thrust-to-power ratio[2,3],which have realized a wide range of space applications and are mainly used for various orbital missions such as attitude control, orbital hoisting, and deep space exploration.In recent years, compared to large communication satellite platforms that have achieved in-orbit applications using a variety of HT in the 1–5 kW power range [4, 5], the microsatellites [6] have drawn significant attention from satellite manufacturers because of their wide range of applications and high economic profitability, and have made rapid developments.With speeding up the emergence of low earth orbit(LEO) Internet constellations and the demand for orbital high-precision maneuvers, the well-known commercial space companies represented by OneWeb [7] and SpaceX [8] have vigorously developed LEO Internet constellations with lowpower Hall thruster (LHT) as propulsion systems for highprecision station keeping and orbit insertion.Attracted by the enormous demand for propulsion systems in commercial aerospace, the 300 W class LHT with an outer diameter of 40 mm (LHT40) was specifically developed by the Lanzhou Institute of Physics (LIP).

    Nevertheless, LHT [9, 10], whose primary performance parameters, such as specific impulse, efficiency, and thrust,depend significantly on the structure of the beam plasma,has large surface-to-volume ratios resulting in poorer performance of LHT than medium and high-power HT.Therefore, the correlation between the plume characteristics and the dimensions of the LHT is highly essential to improve the scaling relationship of the LHT and is of great value for the successful integration of thrusters and spacecraft on-orbit and plume shielding.The plume generated by the electric thruster causes interaction between the spacecraft and the surrounding environment, reducing the performance of the spacecraft subsystems and on-board sensors.The plume containing energetic particles, electrons, neutral atoms, etc will have a variety of adverse effects on the spacecraft: (1) space charge effects result in a charged spacecraft surface; (2) various charged particles change the electrical and thermal properties of the spacecraft surface;(3)plume deposition will contribute to the performance degradation of solar cell panel, thermal control surfaces,and optical windows,etc.Many studies have been conducted in identifying the performance and efficiency of LHT [11–14] with Faraday probe [15–17] and cylindrical Langmuir probe [18–20] which are commonly used to perform medium and high-power HT beam plasma studies.But,a single Faraday probe [15–17] combined with a moving mechanism was employed in previous studies, which cannot determine the simultaneous spatial distribution characteristics of beam plasma.In addition, there are rarely studies on the far-field plasma spatial distribution and evolution characteristics of LHT.The plume plasma characteristics of LHT are an external manifestation of LHT.The plume characteristics of LHT are of great interest to study because of their large surface area to volume ratio.Reconstructing the spatial plasma parameter distribution is of great value for understanding the operating mechanism of LHT.

    In this work, a dedicated and integrated plasma diagnostic system (PDS) composed of seventeen Faraday cups(FC) and two triple Langmuir probes (TLP) is established to investigate the time-averagedin situspatial distribution characteristics of far-field ions and electrons of the 300 W class LHT for commercial aerospace applications.The ion current density (ICD), plasma potential (Vsp), plasma density(Ne), and electron temperature (Te) at 1000 mm downstream of LHT in the azimuthal angle range of ?90° to 90° were investigated under the conditions of different anode mass flow rates and discharge voltages.

    2.Experimental set-up

    2.1.LHT40 HT

    The LHT40-V1, composed of an anode, discharge chamber,magnetic circuit, and a hollow cathode, etc, is the primary version of the 300 W class LHT developed independently by the LIP, which is a type of stationary state plasma thruster.The picture of the LHT40-V1 HT is shown in figure 1.

    The discharge chamber consists of a BN-SiO2annular discharge channel with an external diameter of 40 mm and a width of 10 mm.The propellant gas is uniformly injected inside the channel through the gas distributor and the anode,which consists of a non-magnetic stainless steel ring placed at the bottom of the discharge channel and adjacent to the outer ceramic wall.The magnetic field is generated by electromagnetic coil components(four outer coils, one internal coil,one internal trim coil, and magnetic screen, etc).The magnetic field lines of LHT40-V1 are concave and symmetric about the discharge chamber centerline, which is similar to the NASA-173Mv1 HT[21].Furthermore, the peak of radial magnetic field magnitude of the channel middle diameter(Br)located at the exit of the thruster discharge channel is approximately 140 G and the distribution of magnetic field magnitude at its channel middle diameter shows a positive gradient, as shown in figure 2.The cathode is a 3 A emitting current LaB6hollow hot cathode, located outside the four excitation coils.The typical operating parameters of LHT40-V1 HT are shown in table 1.

    2.2.Plasma diagnostic system (PDS)

    Table 1.Typical operation parameters of LHT40-V1 Hall thruster.

    The dedicated and integrated PDS [22] consists of seventeen FC,two TLP,a rotary rake,an integrated control system,and an advanced post-processing system used to perform investigation of ICD,plasma potential,plasma density,and electron temperature.The rotary rake controlled by the motor is rotated around thez-axis defined by the thruster exit plane(the thruster coordinates are defined as shown in figure 3),and its rotation radius is 1000 mm.The rotation angle range(θ)of the rake is ±90° with rotation precision ±0.1°.The FC is distributed at 0°, ±10°, ±15°, ±20°, ±25°, ±35°, ±45°,±60°, ±75° positions of the rake, as shown by the yellow arrows in figure 3(b) and two TLP are symmetrically distributed at the 0° axis of the rake on both sides, which are at an angle of 10 with theY-axis of the thruster exit plane.The experiment was conducted in the vacuum system (chamber TS-6C) of LIP, which consists of a vacuum chamber with a diameter of 4.2 m and a length of 10 m and equipped with six external cryogenic pumps(1×104l s?1for nitrogen)and 20 internal cryogenic pumps (1.49×104l s?1for xenon gas).The no-load ultimate vacuum pressure of the vacuum system is better than 5.0×10?5Pa, and the operating pressure is better than 1×10?3Pa at a xenon mass flow rate of 8 sccm.One KF-40 ion gauge , which is parallel to the thruster exit plane, is mounted on the vacuum wall.The schematic of the vacuum facility and measuring apparatus is shown in figure 3(a).The relative positions of the PDS and thruster in the vacuum chamber are shown in figure 3(b).

    2.2.1.Faraday cups.The FC is composed of a collector,gathering the ions coming from the plasma source,a polarized shield, acting as a filter for the electrons, and an ammeter measuring the current from the probe collector to the vacuum chamber ground,as shown in figure 4.The main output of an array of FC is the ICD distribution(ICAD).The FC collector and external shield (EXT SHLD) are set to ?40 V.The collection area of a FC is 100 mm2.The FC can be utilized to obtain the total beam ion current(Ii,tot).Whenr>10d(rand

    drepresent the distance between the entrance hole of the probe and the exit plane of the thruster,and outer diameter of the thruster exit cross-section), the plasma plume can be considered to originate from the center of the thruster exit plane (a point source) [23].The total ion current [24, 25] is therefore equal to the integral over the ICDjihemisphere,i.e.

    whereR, θ, andjidenote the distance of the probe collection surface from the center of the thruster exit cross-section, the angle between the central axis of the collection surface of the probe, and the ICD is collected by the probe, respectively.Ii,totwas calculated by the central FC (at 0°) during the rake rotation in the range of ±90° in this paper.

    In this case, the plume divergence half-angle (PDHA,θdiv) is defined as the plume boundary containing 90% of the total ion current (Ii,90%) [25, 26], as shown in equation (2)

    2.2.2.Triple Langmuir probe.The TLP used here, derived from Chen and Sekiguchi’s design [27], consists of three cylindrical tungsten wires (P1, P2, P3) with a diameter of 4 mm and a length of 57 mm.The distance between each two of the three tips is 8 mm.The TLP placed in the plasma bulk consists of three voltage-biased electrodes and allows the evaluation of plasma characteristic parameters such as plasma potential (Vsp), plasma density (Ne), and electron temperature(Te).The picture and principle scheme of TLP is shown in figure 5.The TLP, composed of two bias electrodes (P1, P3)to which a fixed voltage of 32 V is applied and an additional third electrode (as a floating electrode), has the advantage of instantaneous availability of plasma parameters without a complex scanning circuit comparable to that of single Langmuir probe.The TLP is widely used in plasma diagnostics for pulsed plasma thrusters, HT, ion thrusters,and other thrusters.Plasma potential (Vsp), plasma density(Ne), and electron temperature (Te) are derived by a specially developed post-processing program.

    The electron temperature (Te) and electron number density (Ne) can be represented as equations (3) and (4)[22], respectively

    HereV21andIrepresent the voltage between probe electrodes P1 and P2 and the current between electrodes P1 and P3,respectively.Thek,e,S, andmirepresent the Boltzmann constant,the single electron charge,the ion collection area of the probe, and the mass of the xenon atom, respectively.

    The plasma potential (Vsp) can be written as the below[22]:

    whereV21represents the potential difference between probe electrodes P1 and P2,?V1refis the acquired voltage between probe electrode P1 and the external reference, respectively.

    3.Results and discussion

    The estimation of the plasma parameters within the far-field plume of the 300 W LHT (LHT40-V1) is displayed for diverse working conditions.The ICD (Ii), beam divergence angle (θ), mass utilization efficiency (ηm), plasma potential(Vsp),plasma density(Ne),and electron temperature(Te)were investigated under various discharge voltages (100–350 V)and anode mass rates(5–8 sccm).The xenon propellant with a purity of 99.9995% was delivered to the anode and cathode through a dedicated mass flow controller with an accuracy of 1% of full scale (10 sccm).For each operating condition, the discharge current oscillation was recorded.Comparing the discharge current and discharge current oscillations with different probe positions in the plume, it was observed that there was no influence of the discharge behavior regardless of the probe positions.Before the beam plasma diagnosis, the thruster was carried out for a continuous discharge period up to 2 h.The thruster was essentially considered to reach thermal stability,due to the discharge current remaining constant.During the plume test, the entire period of each operating condition (θ from ?90° to +90°) was about 5 min.

    3.1.FC measurement results

    3.1.1.Influence of the anode mass flow rate.The influence of the anode mass flow rate on the plasma parameters of the far-field plume was investigated at a fixed anode voltage of 300 V (Ud=300 V), a cathode mass flow rate of 1 sccm(mc.=1 sccm),and an excited current of 3.2 A(Iec=3.2 A).The range of the rotary rake is ±90° and the angular resolution of the data collected during the rake rotation is 3°.The ICD distribution of the thruster is obtained by the discrete ICDs of the probes at different locations using a linear differential smoothing algorithm.Figure 6 shows the spatial distribution profile of the ICD with anode flow rate(ma.)from 5 to 8 sccm(interval 1 sccm)discharge voltage for a cathode mass flow rate of 1 sccm and a discharge voltage of 300 V.In figure 6, the probe position represents the specific physical distribution and the angle represents the rotation angle in respect to theZ-axis (as shown in figure 3).

    Figure 1.The picture of the LHT40-V1 Hall thruster.

    Figure 2.The distribution of magnetic field strength at the channel middle diameter of the LHT40-V1 Hall thruster (L represents the channel length,and z represents the distance from the anode surface,z/L is a dimensionless number).

    Figure 3.Schematic of the vacuum facility and measuring apparatus(a), the relative positions of the PDS and thruster in the vacuum chamber (b).

    Figure 4.The picture of FC.

    Figure 5.The picture (a) and principle scheme (b) of TLP.

    Figure 6.The distributions of ICD with different anode mass flow rates under cathode mass flow rate of 1 sccm and discharge voltage of 300 V.

    Figure 7.The distributions of ICD with different background pressures.

    Figure 8.The profiles of beam divergence angle (a) and mass utilization efficiency (b) with anode mass flow rates under the cathode mass flow rate of 1 sccm and discharge voltage of 300 V.

    Figure 9.The distribution of ion current density with discharge voltage under cathode mass flow rate of 1 sccm and anode mass flow rate of 7 sccm.

    Figure 10.The profile of beam divergence angle (a), mass utilization efficiency (b), and current utilization efficiency (c) with different discharge voltages under cathode mass flow rate of 1 sccm and anode mass flow rate of 7 sccm.

    Figure 11.The profiles of electron temperature (a), plasma density (b), and plasma potential (c) with different anode mass flow rates at different angles.

    Figure 12.The profiles of electron temperature (a),plasma density(b),and plasma potential(c)at different discharge voltages with TLP 0°position.

    Figure 13.The curves of average ions velocity with the different discharge voltages.

    As can be seen from figure 6, the ICD is symmetrically distributed with the center of the thruster.The maximum value of ICD is located at the thruster axis.The ICD decreases sharply off-axis (e.g.the ICD is 0.78 mA cm?2when FC is located atφθ= =0,0 in figure 6(b).However, the ICD decreases to 0.28 mA cm?2when the angle increases to 42°).As the anode mass flow rate increases from 5 to 8 sccm, the ICD increases from 0.54 to 1.26 mA cm?2, an approximate percentage increase of 133.3%.The main reason for the analysis may be that with the increasing density of neutral atoms in the HT discharge channel, the probability of collision of electrons and atoms increases.To some extent,the more frequent the probability of collisions between electrons and atoms, the higher the number of ions in the plume.The ICD increases with increasing anode mass flow rate under other constant conditions.

    Besides,the wings were observed at approximately±60°for each operating condition and this structure has been identified as a result of the backpressure effect [23] and the shape of this ion current is similar to other medium and high power HT plume profiles [28–30].The distribution of ICD with different background pressure is shown in figure 7.As the anode mass flow rate increases, the back pressure gradually increases and the ‘wings’ phenomenon becomes more and more obvious.It was also noticed that the ICD distribution profile widens with increasing anode mass flow rate.This phenomenon has been observed in plume measurements of the multiple HT, e.g.SPT-100 [28],PPS1350-G [29], and BHT-600 [30].The problem [23] was attributed to the charge exchange (CEX) collisions as a consequence of the increase of the vacuum chamber background pressure as the anode mass flow rate increases for a certain pumping speed.The CEX collision is the result of a fast-moving xenon ion () colliding slow-moving neutral atoms (xeslow), giving a fast atom (xefast) and a slow ion ().In the HT plume, the large curvature of the magnetic field lines[21]leads to thegradual increase in the off-axis region of the radial electric field gradient, which further contributes tobeing accelerated to large angular positions.So, background pressure is an essential consideration in measuring ICD.As the anode mass flow rate increases,the probability of CEX occurring is greater.Therefore, the total ion current is overvalued as the background pressure rises,leading to an overvaluation of the thruster performance.

    The mass utilization efficiency (ηm) is an important parameter to evaluate the design of HT and represents the percentage of neutral atoms converted to ions, which is expressed in equation (6) [23, 25].

    wherem.i,Iiandm.adenote the mass flow rate of ions in the beam,ion current(considering only single charged ions),and atom mass, respectively.

    As can be seen from figure 6,the ICAD is approximately symmetrical with the center of the thruster.TheIi,totwas calculated by the central FC (φ= 0) during the rake rotation in the range of ±90° in this paper.Besides, the total ion current (Ii,tot-2) was also derived from the ICD measured by the remaining FC.The difference betweenIi,totandIi,tot-2was less than 1%.Given that the ‘wings’ are located at a large angle (characterized by a minor ICD), the method of direct integration in the process of calculating the total current is employed.

    The profiles of beam divergence angle and mass utilization efficiency with anode mass flow rates under the cathode mass flow rate of 1 sccm and discharge voltage of 300 V are shown in figure 8.As can be seen from figure 8,the anode mass flow rate increases from 5 to 8 sccm, the mass utilization efficiency increases from 77% to 93%, which is consistent with the ICAD at different anode mass flow rates in figure 6.The beam divergence angle increases from 21.3° to 23.8°, which is associated with the increase in the backpressure of the vacuum chamber due to the increase in the anode mass flow rate further causing an increase in the CEX collisions, improving the beam loss efficiency (θcos2).

    3.1.2.Influence of the discharge voltage.The influence of the discharge voltage on the plasma parameters of the thruster far-field plume was investigated at the fixed anode mass flow rate of 7 sccm (ma.=7 sccm), cathode mass flow rate of 1 sccm(mc.=1 sccm),and an excited current of 3.2 A(Iec=3.2 A).The contours of the spatial distribution of the ICD with the discharge voltage for discharge voltages from 150 to 300 V (50 V interval) are shown in figure 9.

    As expected, it can be noticed from figure 9 that the ICAD of the thruster shows axisymmetric characteristics under each operating condition,with the maximum ICD at the thruster axis and decreasing off-axis in ascending order.Interestingly, as the anode voltage increases from 150 to 300 V,the ion current at the axis keeps increasing.However,when the anode voltage continues to increase to 350 V, the shape and magnitude of the ion current do not continue to increase, but a maximum occurs.

    One of the main reasons analyzed is that the high discharge voltage operation favors the generation of multicharged ions, where the number of charged ions at large angles of the axis increases.It was also observed in the studies conducted by Reid [31], Dannenmaye [32].Another reason could be that the electric field distribution inside the thruster discharge channel is affected by the increase in discharge voltage,which causes the ionization zone inside the thruster to move toward the anode position and the acceleration region to be enlarged [33].A small fraction of the ions at low potentials are accelerated to large angular positions.

    The profiles of beam divergence angle and mass utilization at different discharge voltages under cathode mass flow rate are shown in figure 10.As can be seen from figure 10,discharge voltage increases from 150 to 350 V,the beam divergence angle decreases from 27.36° to 22.8°, mass utilization efficiency decreases from 91%to 82%and current utilization efficiency (ηc=Ii,tol/Id) decreases from 76% to 68%.It could be attributed to increasing the discharge voltage, the ionization zone inside the discharge channel shifting upstream and the acceleration region widening leading to an increase in the proportion of multi-charged ions in the beam.Therefore, the increase in electron temperature leads to an increase in electron leakage to the anode, a decrease in the ion current ratio, and a decrease in mass utilization efficiency.The current utilization efficiency decreases with increasing discharge current, which supports this deduction.

    3.2.TLP measurement results

    Given that the FC measurements show that the thruster plume plasma is symmetrically distributed to the central axis.During the experiments, the TLP measurements of the plasma parameters were scanned in the range of 0°–90°.Since the angle between the TLP is 10°,to satisfy the probe scan range of 0°–90°,the rake rotation is rotated 0°,5°,20°,25°,40°,45°,60°,65°, 80° respectively.

    3.2.1.Influence of the anode mass flow rate.The effects of the anode mass flow rate (from 5 to 8 sccm) on the plasma potential (Vsp), plasma density (Ne) and electron temperature(Te)of the thruster far-field plume were investigated at a fixed anode voltage of 300 V (Ud=300 V), a cathode mass flow rate of 1 sccm(mc.=1 sccm),and an excited current of 3.2 A(Iec=3.2 A).The experiment results of the TLP are shown in figure 11.

    It is evident from figure 11 that the plasma(Vsp),plasma density(Ne),and electron temperature(Te)show a decreasing trend with the increasing angle at the same mass flow rate.Besides, the three parameters have the largest values at the thruster axis position(0°),sharply decrease with an increasing angle when changing from 0° to 30°, and the values of each parameter decreasing from 30° to 90° gradually become smaller.At different anode flow rates, the electron temperature, electron density, and plasma potential decrease sequentially in the range of 0°–30° as the anode flow rate decreases, and the differences between the three values are not significant from 30° to 70°.This is mainly due to the increase of the beam plasma parameter as the anode flow rate increases,the neutral atoms in the discharge channel increase,and the ionization rate increases (details of the measurement results are presented in section 3.1.1).From the FC measurements, it is evident that the beam divergence angle at different anode flow rates is mainly concentrated in the range of 0°–30°,which means that the beam plasma is mainly located within 30°.Therefore, the plasma electron temperature,electron density,and space potential decrease rapidly in the range of angles from 0°to 30°,and there is no significant change in parameters beyond 30° which also verifies the consistency of the FC and TLP probe measurements.

    3.2.2.Influence of the discharge voltage.The profiles of plasma electron temperature (Te), electron density (Ne), and plasma potential (Vsp) at different discharge voltages with TLP 0°position are shown in figure 12.It is observed that the plasma electron temperature (Te), electron density (Ne), and space potential increase with increasing discharge voltage at the same mass flow rate.Besides, the plasma electron temperature, electron density, and space potential change faster with different anode voltages as the anode mass flow rate increases.The probable reason for the analysis is that as the discharge voltage increases,the higher the electric field of the thruster discharge channel, the higher the ion kinetic energy of the plasma.The electron temperature at higher voltages may be due to the higher kinetic energy of the ejected ions which neutralize electrons emitted from the cathode to maintain quasi-neutrality.

    ICD (θj( )), plasma potential (Vsp), and plasma density(Ne)were obtained by FC and TLP,respectively.The average ion speed(AIV)was deduced from the measured data of TLP and FC.The AIV (Vprobe) (considering only single charged ions) can be represented as equation (7)

    From the principle of HT operation,it is well known that the ions are accelerated through the electric field.The maximum AIV(Vacc)by the electric field of ions(considering only single charged ions) can be represented as equation (8)

    To check the consistency between the TLP results and FC results,the AIV(Vacc,Vprobe)obtained by the two methods were calculated and compared.Figure 13 shows the profiles of the AIV with the different discharge voltages at the fixed anode mass flow rate of 7 sccm, cathode mass flow rate of 1 sccm, and an excited current of 3.2 A.

    As can be seen from figure 12, the range of AIV calculated by different methods for different discharge voltages is from 8 to 22.3 km s?1.Besides,the AIV increases gradually with the increase of discharge voltage.The discrepancy of AIV calculated by the two methods is less than 10%with the discharge voltage from 150 to 350 V.The maximum error is better than 31.5%for a discharge voltage of 100 V, which checks the consistency between the TLP’s results and that of FC to some extent.The main reason for this error may be the low ionization efficiency of the thruster and the ineffective acceleration of the ions in the plume at this discharge voltage, both of which combined lead to a large error in the plasma parameters calculated by the probe.

    4.Conclusions

    The dedicated and integrated PDS composed of seventeen FC and two TLP is established to investigate simultaneously and instantaneously the far-field time-averagedin situplasma characteristics of 300 W class LHT for commercial aerospace applications.The ICD(Ii),beam divergence angle(θ),plasma potential (Vsp), and plasma density (Ne), as well as electron temperature (Te) in the range of ?90° to 90° at an angle of 1000 mm downstream of the thruster exit, were investigated under the conditions of anode mass flow rate(5–8 sccm)and discharge voltage (100–350 V).The ICD, beam divergence angle, and mass utilization efficiency increase with anode mass flow rate.The wings were also observed, which is related to the increase of the anode mass flow rate and the increase of the background pressure leading to the increase of the CEX collisions.The beam divergence angle and mass utilization efficiency decrease with discharge voltage, which may be attributed to the fact the high discharge voltage operation favors the generation of multi-charged ions.The plasma electron temperature, electron density, and space potential decrease rapidly in the range of angles from 0° to 30°, and there is no significant change in parameters beyond 30°,which also verifies the thruster beam divergence angle is less than 30° under various operating conditions and the consistency of the FC and TLP probe measurements.The discrepancy of AIV was calculated by the two methods.The maximum error is better than 31.5%for a discharge voltage of 100 V, which checks the consistency between the TLP’s results and that of FC to some extent.The experimental data provide insight into the physics of the plasma of the crossfield discharge of the LHT, provide experimental data boundary input to the thruster plume model.At the same time,it is also of great importance for the optimal design and space flight application of the thruster.

    Acknowledgments

    The authors would like to acknowledge National Natural Science Foundation of China (Nos.12005087 and 61901204), the Science and Technology Plan of Gansu Province (No.20JR10RA478), the Military Test Instruments Program (No.2006ZCTF0054) and the Key Laboratory Funds for Science and Technology on Vacuum Technology and Physics Laboratory (No.HTKJ2019KL510003).

    猜你喜歡
    劉超
    老吳的拉面館
    Theoretical study of novel B-C-O compounds with non-diamond isoelectronic
    Removal of GaN film over AlGaN with inductively coupled BCl3/Ar atomic layer etch
    征程萬(wàn)里,初心如一
    雷鋒(2021年12期)2021-04-12 00:57:22
    《站立式起跑》教學(xué)設(shè)計(jì)
    小數(shù)與小數(shù)點(diǎn)
    小數(shù)與小數(shù)點(diǎn)
    Design and experiment of the centrifugal pump impellers with twisted inlet vice blades *
    A simple method for estimating bed shear stress in smooth and vegetated compound channels*
    A randomized, controlled trial of Shutangluo fang in treating painful diabetic peripheral neuropathy
    久久久水蜜桃国产精品网| 亚洲免费av在线视频| 国产精品偷伦视频观看了| 国产成人欧美在线观看 | 国产真人三级小视频在线观看| 久久天堂一区二区三区四区| 满18在线观看网站| 久久草成人影院| 真人做人爱边吃奶动态| 在线观看66精品国产| 波多野结衣av一区二区av| 久久中文看片网| 最近最新中文字幕大全免费视频| 精品久久蜜臀av无| 无限看片的www在线观看| 欧美日韩亚洲高清精品| 一二三四社区在线视频社区8| 人成视频在线观看免费观看| 成人18禁在线播放| 国产亚洲av高清不卡| 午夜福利乱码中文字幕| 午夜福利在线观看吧| 18禁裸乳无遮挡动漫免费视频| 男男h啪啪无遮挡| 国产高清视频在线播放一区| 亚洲中文av在线| 中文欧美无线码| 成人影院久久| 天堂俺去俺来也www色官网| 久久亚洲精品不卡| 大型av网站在线播放| 久久久国产成人精品二区 | 老熟妇仑乱视频hdxx| 亚洲av美国av| 亚洲一卡2卡3卡4卡5卡精品中文| 成人精品一区二区免费| 日韩免费av在线播放| 亚洲精品美女久久av网站| 如日韩欧美国产精品一区二区三区| 最近最新中文字幕大全电影3 | 日本a在线网址| 精品无人区乱码1区二区| 啦啦啦视频在线资源免费观看| 999精品在线视频| 亚洲,欧美精品.| 亚洲精品自拍成人| 91老司机精品| 女人被躁到高潮嗷嗷叫费观| 91成年电影在线观看| 黄色片一级片一级黄色片| 窝窝影院91人妻| 欧美日韩黄片免| 免费在线观看完整版高清| 这个男人来自地球电影免费观看| 国产人伦9x9x在线观看| 老汉色∧v一级毛片| 99精品欧美一区二区三区四区| 十八禁人妻一区二区| 国产高清videossex| 亚洲专区国产一区二区| 十八禁网站免费在线| 国产av又大| 黄色a级毛片大全视频| 999久久久精品免费观看国产| 欧美激情 高清一区二区三区| 高清毛片免费观看视频网站 | 亚洲 欧美一区二区三区| 超碰97精品在线观看| 精品人妻熟女毛片av久久网站| 免费观看a级毛片全部| 久久青草综合色| 激情视频va一区二区三区| 日韩欧美三级三区| 99久久国产精品久久久| 美国免费a级毛片| videos熟女内射| 少妇裸体淫交视频免费看高清 | 国产精品一区二区免费欧美| 好男人电影高清在线观看| 色94色欧美一区二区| 国产精品秋霞免费鲁丝片| 亚洲av片天天在线观看| 99re在线观看精品视频| 久久天躁狠狠躁夜夜2o2o| 脱女人内裤的视频| 热re99久久精品国产66热6| 老汉色av国产亚洲站长工具| 亚洲成a人片在线一区二区| 咕卡用的链子| 在线永久观看黄色视频| 99riav亚洲国产免费| 亚洲成a人片在线一区二区| 亚洲av成人av| 欧美精品亚洲一区二区| 亚洲欧美日韩另类电影网站| 成年版毛片免费区| 精品久久久精品久久久| 成人永久免费在线观看视频| 一边摸一边抽搐一进一出视频| 亚洲精品中文字幕在线视频| 久久亚洲精品不卡| 俄罗斯特黄特色一大片| 成人手机av| 色综合欧美亚洲国产小说| 丝袜美足系列| 成人特级黄色片久久久久久久| 久久国产精品男人的天堂亚洲| 国产午夜精品久久久久久| 电影成人av| 国产精品亚洲av一区麻豆| 法律面前人人平等表现在哪些方面| 黄片大片在线免费观看| 亚洲色图 男人天堂 中文字幕| 男女之事视频高清在线观看| 麻豆av在线久日| 久久久久久亚洲精品国产蜜桃av| 久久ye,这里只有精品| 咕卡用的链子| 欧美亚洲日本最大视频资源| 欧美在线黄色| 色94色欧美一区二区| 91国产中文字幕| √禁漫天堂资源中文www| 亚洲av成人一区二区三| 国产99白浆流出| 99re在线观看精品视频| 精品国产一区二区三区久久久樱花| 这个男人来自地球电影免费观看| 老司机福利观看| 一本综合久久免费| 999久久久国产精品视频| 高清视频免费观看一区二区| 91精品三级在线观看| 女人爽到高潮嗷嗷叫在线视频| 午夜视频精品福利| 天堂动漫精品| 久久午夜亚洲精品久久| 一级作爱视频免费观看| 成年女人毛片免费观看观看9 | 免费在线观看视频国产中文字幕亚洲| 婷婷成人精品国产| 久久久国产欧美日韩av| 国产又色又爽无遮挡免费看| 久久人妻av系列| 久久久久久久午夜电影 | 19禁男女啪啪无遮挡网站| 一a级毛片在线观看| 亚洲国产看品久久| 亚洲第一青青草原| 天天操日日干夜夜撸| 十分钟在线观看高清视频www| 精品国产乱子伦一区二区三区| 欧美精品亚洲一区二区| 国产精品.久久久| 久99久视频精品免费| 91大片在线观看| 黄色怎么调成土黄色| 美女国产高潮福利片在线看| 老熟妇仑乱视频hdxx| 国产在线一区二区三区精| 99精品欧美一区二区三区四区| 一级,二级,三级黄色视频| 91麻豆av在线| 亚洲va日本ⅴa欧美va伊人久久| 久久精品aⅴ一区二区三区四区| 午夜精品在线福利| 丁香六月欧美| 亚洲一区中文字幕在线| 女同久久另类99精品国产91| 欧美日韩av久久| 国产亚洲精品久久久久久毛片 | 欧美色视频一区免费| 老司机福利观看| 亚洲自偷自拍图片 自拍| 韩国精品一区二区三区| 久热这里只有精品99| 亚洲九九香蕉| 少妇被粗大的猛进出69影院| 国产成人欧美| 91成年电影在线观看| 亚洲成人手机| 亚洲精品在线观看二区| 深夜精品福利| 人成视频在线观看免费观看| 国产在线精品亚洲第一网站| 欧美亚洲日本最大视频资源| av有码第一页| 久久99一区二区三区| 国产片内射在线| 亚洲欧美日韩另类电影网站| 老汉色av国产亚洲站长工具| 国产成人一区二区三区免费视频网站| 午夜日韩欧美国产| 日韩成人在线观看一区二区三区| 交换朋友夫妻互换小说| 黑人猛操日本美女一级片| 色老头精品视频在线观看| 久久久久久久国产电影| 人妻 亚洲 视频| 丝袜美腿诱惑在线| 久久午夜亚洲精品久久| 亚洲国产中文字幕在线视频| 久久精品国产综合久久久| 超碰97精品在线观看| 日韩欧美一区视频在线观看| 黑丝袜美女国产一区| 69精品国产乱码久久久| 国产高清国产精品国产三级| 不卡av一区二区三区| 久久精品国产a三级三级三级| 亚洲男人天堂网一区| 99国产精品一区二区蜜桃av | 久久精品亚洲精品国产色婷小说| 中文字幕制服av| 亚洲男人天堂网一区| 亚洲精华国产精华精| 国产97色在线日韩免费| 日日夜夜操网爽| 动漫黄色视频在线观看| 一区福利在线观看| 欧美激情高清一区二区三区| 99热国产这里只有精品6| 久久影院123| 又紧又爽又黄一区二区| 国产淫语在线视频| 久久中文字幕一级| 很黄的视频免费| 啦啦啦在线免费观看视频4| 亚洲五月色婷婷综合| av在线播放免费不卡| 久久久久久久精品吃奶| 中文字幕最新亚洲高清| 女同久久另类99精品国产91| 老汉色av国产亚洲站长工具| 国产亚洲欧美98| 韩国av一区二区三区四区| 国产深夜福利视频在线观看| 91精品三级在线观看| 中文欧美无线码| 国产日韩一区二区三区精品不卡| 美女国产高潮福利片在线看| 丁香欧美五月| 国内毛片毛片毛片毛片毛片| bbb黄色大片| 人妻丰满熟妇av一区二区三区 | 国产一卡二卡三卡精品| 国产精品久久久久久人妻精品电影| 国产97色在线日韩免费| 搡老熟女国产l中国老女人| 男人舔女人的私密视频| 天堂动漫精品| 亚洲色图 男人天堂 中文字幕| 国产日韩一区二区三区精品不卡| 亚洲专区国产一区二区| 日韩欧美一区视频在线观看| 久久久久久久精品吃奶| 欧美一级毛片孕妇| 欧美精品啪啪一区二区三区| 波多野结衣av一区二区av| 久久人妻福利社区极品人妻图片| 久久久精品国产亚洲av高清涩受| 亚洲中文日韩欧美视频| 丝瓜视频免费看黄片| 99国产精品免费福利视频| 大香蕉久久成人网| xxx96com| 麻豆乱淫一区二区| 国产一区二区三区视频了| 国产欧美日韩一区二区三| 日本精品一区二区三区蜜桃| 天天操日日干夜夜撸| 夫妻午夜视频| 少妇 在线观看| 两个人免费观看高清视频| 久久久久久久久久久久大奶| 国产成人啪精品午夜网站| 岛国在线观看网站| 日韩精品免费视频一区二区三区| av网站免费在线观看视频| 午夜久久久在线观看| 亚洲一区中文字幕在线| 亚洲片人在线观看| av网站免费在线观看视频| 国产精品久久电影中文字幕 | 91成人精品电影| 男女之事视频高清在线观看| 91麻豆av在线| 777久久人妻少妇嫩草av网站| 精品视频人人做人人爽| 国产精品 国内视频| 亚洲精品国产色婷婷电影| 女同久久另类99精品国产91| 亚洲,欧美精品.| 国产欧美日韩一区二区精品| 精品国产乱子伦一区二区三区| 欧美 日韩 精品 国产| av有码第一页| 操美女的视频在线观看| 不卡av一区二区三区| 一区福利在线观看| 1024视频免费在线观看| 国产高清videossex| 一边摸一边做爽爽视频免费| 亚洲va日本ⅴa欧美va伊人久久| 97人妻天天添夜夜摸| 免费少妇av软件| 欧美精品人与动牲交sv欧美| 法律面前人人平等表现在哪些方面| 久久久久久人人人人人| 国产一区二区激情短视频| 青草久久国产| 黑人欧美特级aaaaaa片| 大型黄色视频在线免费观看| 国产有黄有色有爽视频| 欧美日韩视频精品一区| 精品欧美一区二区三区在线| 精品少妇久久久久久888优播| 欧美性长视频在线观看| 每晚都被弄得嗷嗷叫到高潮| 国产精品久久久av美女十八| 欧美黑人欧美精品刺激| 十八禁网站免费在线| 69av精品久久久久久| 视频区图区小说| 免费av中文字幕在线| 精品卡一卡二卡四卡免费| 午夜两性在线视频| 国产精品久久久av美女十八| 成年人免费黄色播放视频| 热99re8久久精品国产| 亚洲午夜精品一区,二区,三区| 精品久久久久久,| 午夜福利免费观看在线| 久99久视频精品免费| 国产99久久九九免费精品| 九色亚洲精品在线播放| 亚洲国产看品久久| 男女免费视频国产| 中文字幕色久视频| 欧美人与性动交α欧美软件| 91成人精品电影| 欧美日韩亚洲综合一区二区三区_| 精品无人区乱码1区二区| 一级毛片精品| 亚洲国产欧美日韩在线播放| 亚洲久久久国产精品| 国产色视频综合| 精品一区二区三区av网在线观看| 久久香蕉精品热| 宅男免费午夜| 又大又爽又粗| 亚洲欧洲精品一区二区精品久久久| 下体分泌物呈黄色| 精品一品国产午夜福利视频| 免费人成视频x8x8入口观看| 国产精品免费大片| 欧美不卡视频在线免费观看 | xxx96com| 妹子高潮喷水视频| 天堂动漫精品| 在线观看日韩欧美| 99国产精品99久久久久| 99久久人妻综合| 精品国产一区二区三区久久久樱花| 精品无人区乱码1区二区| 久久精品国产综合久久久| 男女午夜视频在线观看| 国产亚洲精品久久久久5区| 首页视频小说图片口味搜索| 欧美 亚洲 国产 日韩一| 亚洲国产精品sss在线观看 | 久久久国产欧美日韩av| 国产亚洲精品久久久久5区| 777米奇影视久久| 欧美日韩精品网址| 国产精品国产av在线观看| 桃红色精品国产亚洲av| 亚洲情色 制服丝袜| 悠悠久久av| 一个人免费在线观看的高清视频| 夜夜夜夜夜久久久久| 免费在线观看视频国产中文字幕亚洲| 国产真人三级小视频在线观看| 在线看a的网站| 91麻豆av在线| 手机成人av网站| 国产av一区二区精品久久| 在线观看免费日韩欧美大片| 免费在线观看日本一区| 老熟妇仑乱视频hdxx| 国产男女超爽视频在线观看| 成人国语在线视频| 欧美乱色亚洲激情| 男女午夜视频在线观看| 18禁国产床啪视频网站| 在线观看免费日韩欧美大片| 亚洲一区高清亚洲精品| 成人精品一区二区免费| 男女午夜视频在线观看| 男女下面插进去视频免费观看| 国内久久婷婷六月综合欲色啪| videos熟女内射| 久99久视频精品免费| 欧美激情久久久久久爽电影 | 最近最新中文字幕大全电影3 | 精品福利观看| 欧美午夜高清在线| av福利片在线| 中文字幕精品免费在线观看视频| 久热这里只有精品99| 国产av又大| 女人爽到高潮嗷嗷叫在线视频| 99re在线观看精品视频| 桃红色精品国产亚洲av| 久久国产亚洲av麻豆专区| 国产成人系列免费观看| 一级,二级,三级黄色视频| 一级毛片女人18水好多| 精品亚洲成国产av| 麻豆av在线久日| 看片在线看免费视频| 巨乳人妻的诱惑在线观看| 18在线观看网站| 美女福利国产在线| 国产不卡一卡二| 成人黄色视频免费在线看| 99久久99久久久精品蜜桃| 一区福利在线观看| 两性夫妻黄色片| 人人妻人人澡人人爽人人夜夜| 免费在线观看黄色视频的| 不卡av一区二区三区| 欧美黄色片欧美黄色片| 亚洲人成电影免费在线| 又黄又粗又硬又大视频| 国产91精品成人一区二区三区| 国产国语露脸激情在线看| 美女国产高潮福利片在线看| 久久久久精品人妻al黑| 欧美久久黑人一区二区| 国产欧美日韩一区二区三| 成人18禁在线播放| 一区二区日韩欧美中文字幕| av片东京热男人的天堂| 9热在线视频观看99| 中文字幕人妻丝袜一区二区| 老司机午夜十八禁免费视频| 亚洲自偷自拍图片 自拍| av视频免费观看在线观看| 99精品久久久久人妻精品| 欧美人与性动交α欧美软件| 国产欧美日韩一区二区三| 久99久视频精品免费| 欧美精品啪啪一区二区三区| 精品乱码久久久久久99久播| 精品久久久久久电影网| 不卡av一区二区三区| 亚洲精品国产一区二区精华液| 国产人伦9x9x在线观看| 国产亚洲一区二区精品| 久久久久久久午夜电影 | 18禁美女被吸乳视频| 黑人猛操日本美女一级片| 国产无遮挡羞羞视频在线观看| 亚洲av日韩精品久久久久久密| 久久精品国产99精品国产亚洲性色 | 免费观看a级毛片全部| 黄网站色视频无遮挡免费观看| 国产黄色免费在线视频| 国产乱人伦免费视频| 51午夜福利影视在线观看| 精品国产国语对白av| av有码第一页| 午夜久久久在线观看| 国产精品久久久av美女十八| 美女 人体艺术 gogo| 成人av一区二区三区在线看| 高清毛片免费观看视频网站 | 91av网站免费观看| svipshipincom国产片| 久久久久久久国产电影| 亚洲精品中文字幕一二三四区| 久久ye,这里只有精品| 中文字幕制服av| 在线看a的网站| 大码成人一级视频| 日韩中文字幕欧美一区二区| 欧美黄色片欧美黄色片| 亚洲色图 男人天堂 中文字幕| 看黄色毛片网站| 无遮挡黄片免费观看| 又大又爽又粗| 丝袜美足系列| 国产蜜桃级精品一区二区三区 | 手机成人av网站| 欧美激情极品国产一区二区三区| 久久久国产成人精品二区 | 色综合婷婷激情| avwww免费| 一区在线观看完整版| 在线十欧美十亚洲十日本专区| 99精品欧美一区二区三区四区| 熟女少妇亚洲综合色aaa.| 国产在线一区二区三区精| 午夜精品在线福利| 国产精品99久久99久久久不卡| videosex国产| 国产精品一区二区在线不卡| 国产又爽黄色视频| 免费在线观看亚洲国产| 精品国产超薄肉色丝袜足j| 国产精品久久久久久人妻精品电影| 精品卡一卡二卡四卡免费| 99国产精品免费福利视频| 国产精品欧美亚洲77777| 精品久久久久久,| 久久久精品国产亚洲av高清涩受| 精品国产超薄肉色丝袜足j| 老鸭窝网址在线观看| 亚洲自偷自拍图片 自拍| 色94色欧美一区二区| 淫妇啪啪啪对白视频| 久久久精品免费免费高清| 可以免费在线观看a视频的电影网站| 天天躁夜夜躁狠狠躁躁| 久久九九热精品免费| 国产精品一区二区精品视频观看| 中文字幕人妻丝袜制服| 亚洲av片天天在线观看| 国产欧美亚洲国产| 91大片在线观看| 久久精品国产亚洲av香蕉五月 | 国产亚洲精品第一综合不卡| 中文字幕人妻熟女乱码| 亚洲欧美激情在线| 在线看a的网站| 亚洲国产欧美网| 在线免费观看的www视频| 成年女人毛片免费观看观看9 | 搡老岳熟女国产| 人人妻人人澡人人看| 麻豆乱淫一区二区| 狠狠狠狠99中文字幕| 在线国产一区二区在线| 欧美激情 高清一区二区三区| 婷婷丁香在线五月| 韩国精品一区二区三区| 国产色视频综合| 国产精品1区2区在线观看. | 人成视频在线观看免费观看| 黑人操中国人逼视频| 捣出白浆h1v1| 男女高潮啪啪啪动态图| 免费少妇av软件| 中文字幕av电影在线播放| 黄片小视频在线播放| 少妇的丰满在线观看| 嫩草影视91久久| 麻豆国产av国片精品| 在线观看免费视频日本深夜| 高清黄色对白视频在线免费看| 亚洲午夜理论影院| 欧美激情极品国产一区二区三区| 最近最新中文字幕大全电影3 | 成人国语在线视频| 久久人人97超碰香蕉20202| 免费观看a级毛片全部| 久久久精品区二区三区| 电影成人av| 在线观看免费视频日本深夜| 亚洲国产欧美一区二区综合| 精品福利观看| 大型av网站在线播放| 曰老女人黄片| 欧美精品人与动牲交sv欧美| 成年人午夜在线观看视频| 国产成人精品无人区| 免费在线观看视频国产中文字幕亚洲| 大陆偷拍与自拍| 91九色精品人成在线观看| 国产精品一区二区在线观看99| 免费观看精品视频网站| 亚洲视频免费观看视频| 久久久久国产一级毛片高清牌| 乱人伦中国视频| 在线观看免费午夜福利视频| 成人特级黄色片久久久久久久| 亚洲精品自拍成人| 超碰成人久久| 久久中文字幕一级| 色播在线永久视频| 国产免费现黄频在线看| 黄色毛片三级朝国网站| 天堂√8在线中文| 国产成人av激情在线播放| 亚洲,欧美精品.| 男女高潮啪啪啪动态图| 午夜福利乱码中文字幕| 国产亚洲欧美精品永久| 亚洲美女黄片视频| 在线av久久热| videos熟女内射| e午夜精品久久久久久久| 日韩欧美国产一区二区入口| 国产欧美日韩综合在线一区二区| 亚洲全国av大片| 操美女的视频在线观看| 波多野结衣av一区二区av| 99久久综合精品五月天人人| 一a级毛片在线观看| 免费久久久久久久精品成人欧美视频| 久久精品国产亚洲av香蕉五月 | 久久久国产欧美日韩av| 欧美日韩亚洲国产一区二区在线观看 | 纯流量卡能插随身wifi吗|