• <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
    精品国产乱子伦一区二区三区| 欧洲精品卡2卡3卡4卡5卡区| av电影中文网址| 麻豆成人av在线观看| 99国产极品粉嫩在线观看| 婷婷六月久久综合丁香| 丝袜美腿诱惑在线| 男女之事视频高清在线观看| 黄色丝袜av网址大全| 午夜福利成人在线免费观看| 色av中文字幕| 国产人伦9x9x在线观看| 国产成人欧美在线观看| e午夜精品久久久久久久| 国产97色在线日韩免费| www.www免费av| 男人操女人黄网站| 1024手机看黄色片| 国产亚洲欧美98| 欧美另类亚洲清纯唯美| 级片在线观看| 免费在线观看日本一区| 欧美 亚洲 国产 日韩一| 欧美不卡视频在线免费观看 | avwww免费| 免费看美女性在线毛片视频| 性欧美人与动物交配| 99久久99久久久精品蜜桃| 亚洲专区中文字幕在线| 亚洲av成人av| 国产极品粉嫩免费观看在线| 久久久久久久久中文| 一进一出抽搐gif免费好疼| 99热只有精品国产| 国产熟女午夜一区二区三区| 国产精品电影一区二区三区| www.自偷自拍.com| 91av网站免费观看| 精品欧美一区二区三区在线| 欧美日本视频| 日韩 欧美 亚洲 中文字幕| 香蕉丝袜av| 很黄的视频免费| 日日爽夜夜爽网站| av天堂在线播放| 亚洲av日韩精品久久久久久密| 麻豆成人午夜福利视频| 日韩大尺度精品在线看网址| 91字幕亚洲| 亚洲精品美女久久av网站| 成人手机av| 国产高清videossex| 欧美精品亚洲一区二区| 国产精品免费一区二区三区在线| 后天国语完整版免费观看| 九色国产91popny在线| 中文亚洲av片在线观看爽| 少妇被粗大的猛进出69影院| 禁无遮挡网站| 久久热在线av| 男人舔女人的私密视频| 久久中文字幕一级| 欧美黄色淫秽网站| 白带黄色成豆腐渣| 黄片小视频在线播放| 在线观看午夜福利视频| 中文字幕最新亚洲高清| 12—13女人毛片做爰片一| 国产精品一区二区三区四区久久 | 女性被躁到高潮视频| 欧美av亚洲av综合av国产av| 日本 欧美在线| 日本三级黄在线观看| 亚洲男人的天堂狠狠| 国产熟女xx| 久久人人精品亚洲av| 欧美日韩亚洲综合一区二区三区_| x7x7x7水蜜桃| 757午夜福利合集在线观看| 久久狼人影院| 18禁裸乳无遮挡免费网站照片 | 国产亚洲精品第一综合不卡| 亚洲七黄色美女视频| 免费看日本二区| 日韩国内少妇激情av| 波多野结衣av一区二区av| 性色av乱码一区二区三区2| 精品欧美一区二区三区在线| 777久久人妻少妇嫩草av网站| 一区福利在线观看| 一边摸一边抽搐一进一小说| 一级毛片女人18水好多| www.自偷自拍.com| 国产精品久久久久久人妻精品电影| 亚洲av电影不卡..在线观看| 亚洲 欧美 日韩 在线 免费| 精品久久久久久久毛片微露脸| 免费在线观看成人毛片| 女性被躁到高潮视频| 婷婷六月久久综合丁香| 欧美日韩福利视频一区二区| 国产伦在线观看视频一区| 亚洲成人精品中文字幕电影| 国产精品久久视频播放| 国产激情欧美一区二区| 欧美+亚洲+日韩+国产| 久久久久久久久久黄片| 三级毛片av免费| 1024手机看黄色片| 久久午夜亚洲精品久久| 欧美成狂野欧美在线观看| 老司机午夜十八禁免费视频| 国产97色在线日韩免费| 日本一本二区三区精品| 香蕉久久夜色| 亚洲国产欧美日韩在线播放| 天天躁狠狠躁夜夜躁狠狠躁| 国产成年人精品一区二区| 国内毛片毛片毛片毛片毛片| 午夜免费激情av| www国产在线视频色| 免费观看精品视频网站| 精品久久久久久久久久免费视频| 亚洲欧美日韩高清在线视频| 中文字幕久久专区| 岛国在线观看网站| 欧美黑人欧美精品刺激| 国产97色在线日韩免费| 夜夜看夜夜爽夜夜摸| 99在线人妻在线中文字幕| 亚洲成人久久爱视频| 成人永久免费在线观看视频| 久久精品国产清高在天天线| 免费在线观看亚洲国产| 黄色成人免费大全| 色综合站精品国产| 亚洲av电影不卡..在线观看| 日韩欧美 国产精品| 一个人免费在线观看的高清视频| 丰满人妻熟妇乱又伦精品不卡| 中文亚洲av片在线观看爽| 亚洲精品粉嫩美女一区| 精品午夜福利视频在线观看一区| 亚洲电影在线观看av| 欧美激情极品国产一区二区三区| 校园春色视频在线观看| 色哟哟哟哟哟哟| 久久久国产精品麻豆| 首页视频小说图片口味搜索| 夜夜看夜夜爽夜夜摸| 琪琪午夜伦伦电影理论片6080| 久久草成人影院| 国产精品免费视频内射| 午夜福利免费观看在线| 在线播放国产精品三级| 亚洲真实伦在线观看| 母亲3免费完整高清在线观看| 制服人妻中文乱码| 村上凉子中文字幕在线| 最近最新免费中文字幕在线| 一级毛片高清免费大全| 香蕉av资源在线| 午夜精品久久久久久毛片777| 搡老熟女国产l中国老女人| 精品国产一区二区三区四区第35| 亚洲一卡2卡3卡4卡5卡精品中文| 美国免费a级毛片| 国产又爽黄色视频| 两性夫妻黄色片| 亚洲成国产人片在线观看| 欧美 亚洲 国产 日韩一| 亚洲成人国产一区在线观看| 国产精品亚洲av一区麻豆| 激情在线观看视频在线高清| www日本黄色视频网| 亚洲第一青青草原| 国产精品九九99| 欧美乱妇无乱码| 色综合欧美亚洲国产小说| 成人18禁高潮啪啪吃奶动态图| 日本熟妇午夜| 女人高潮潮喷娇喘18禁视频| 亚洲熟妇熟女久久| 青草久久国产| 日本三级黄在线观看| 一本大道久久a久久精品| 欧美+亚洲+日韩+国产| 国产精品二区激情视频| 亚洲一码二码三码区别大吗| 亚洲欧美激情综合另类| 久久 成人 亚洲| 亚洲一区二区三区色噜噜| 黄色女人牲交| 老鸭窝网址在线观看| 国产97色在线日韩免费| 91在线观看av| 人人妻人人澡欧美一区二区| 91成年电影在线观看| 真人一进一出gif抽搐免费| 亚洲熟妇熟女久久| 男人舔奶头视频| 日韩欧美一区二区三区在线观看| 国内毛片毛片毛片毛片毛片| 俺也久久电影网| 色av中文字幕| 女生性感内裤真人,穿戴方法视频| 琪琪午夜伦伦电影理论片6080| 免费人成视频x8x8入口观看| xxxwww97欧美| 真人做人爱边吃奶动态| 999精品在线视频| 宅男免费午夜| 免费搜索国产男女视频| 在线av久久热| 国产亚洲精品一区二区www| 日本五十路高清| 国产单亲对白刺激| 中文字幕精品免费在线观看视频| 欧美又色又爽又黄视频| 视频区欧美日本亚洲| 久久性视频一级片| 正在播放国产对白刺激| 久久精品91蜜桃| 国产91精品成人一区二区三区| 亚洲熟女毛片儿| 国产精品亚洲一级av第二区| 亚洲激情在线av| 久久人妻福利社区极品人妻图片| 在线观看日韩欧美| 国产在线精品亚洲第一网站| 国产精品98久久久久久宅男小说| 99re在线观看精品视频| 最近最新中文字幕大全免费视频| 日本撒尿小便嘘嘘汇集6| 搡老熟女国产l中国老女人| 亚洲真实伦在线观看| 午夜福利免费观看在线| 亚洲人成伊人成综合网2020| 香蕉久久夜色| 国产三级在线视频| 久久人妻福利社区极品人妻图片| 中出人妻视频一区二区| 亚洲最大成人中文| 国产av一区二区精品久久| 久久精品成人免费网站| 国产亚洲精品一区二区www| 日本精品一区二区三区蜜桃| 国产人伦9x9x在线观看| 巨乳人妻的诱惑在线观看| 中文字幕人妻熟女乱码| 亚洲专区字幕在线| 午夜福利18| 91国产中文字幕| 制服人妻中文乱码| 女人高潮潮喷娇喘18禁视频| 国产麻豆成人av免费视频| 欧美成人午夜精品| 免费搜索国产男女视频| 欧美国产精品va在线观看不卡| 国产欧美日韩精品亚洲av| 亚洲男人的天堂狠狠| av天堂在线播放| 后天国语完整版免费观看| 欧美精品亚洲一区二区| 大型av网站在线播放| 久久国产亚洲av麻豆专区| 宅男免费午夜| 国产片内射在线| 国产成人系列免费观看| 最新在线观看一区二区三区| 婷婷精品国产亚洲av在线| 精品欧美国产一区二区三| 激情在线观看视频在线高清| 日日干狠狠操夜夜爽| 国产激情偷乱视频一区二区| 中文字幕人妻丝袜一区二区| 日韩欧美国产在线观看| 国产免费av片在线观看野外av| 最新在线观看一区二区三区| 婷婷精品国产亚洲av在线| 精品欧美国产一区二区三| 亚洲专区国产一区二区| 69av精品久久久久久| 人人妻人人看人人澡| 免费在线观看完整版高清| 精品午夜福利视频在线观看一区| 亚洲国产精品sss在线观看| 三级毛片av免费| 国产一区二区三区视频了| 久久天躁狠狠躁夜夜2o2o| 无限看片的www在线观看| videosex国产| 国产激情欧美一区二区| 久久久水蜜桃国产精品网| 在线播放国产精品三级| 少妇被粗大的猛进出69影院| 国内久久婷婷六月综合欲色啪| 村上凉子中文字幕在线| 国产区一区二久久| 女人高潮潮喷娇喘18禁视频| aaaaa片日本免费| 久久久久久人人人人人| 国产真实乱freesex| 熟女电影av网| 嫁个100分男人电影在线观看| 亚洲色图av天堂| 99国产综合亚洲精品| 窝窝影院91人妻| 国产精品二区激情视频| 亚洲欧美日韩高清在线视频| 亚洲国产欧洲综合997久久, | 午夜免费成人在线视频| 国产成人av教育| 999久久久精品免费观看国产| 视频在线观看一区二区三区| 麻豆成人av在线观看| 色播亚洲综合网| 超碰成人久久| 国产亚洲精品久久久久5区| 色播在线永久视频| 天天添夜夜摸| 好看av亚洲va欧美ⅴa在| 国产黄a三级三级三级人| 人人澡人人妻人| 国产私拍福利视频在线观看| 午夜久久久在线观看| 久久国产精品男人的天堂亚洲| 国产日本99.免费观看| 国产高清视频在线播放一区| 亚洲成人久久爱视频| 丝袜在线中文字幕| 欧美日韩福利视频一区二区| 大型黄色视频在线免费观看| 美女 人体艺术 gogo| 免费在线观看影片大全网站| 亚洲人成电影免费在线| ponron亚洲| 精品欧美国产一区二区三| 亚洲电影在线观看av| 国产成人影院久久av| 少妇 在线观看| 亚洲一码二码三码区别大吗| 白带黄色成豆腐渣| 1024视频免费在线观看| 午夜两性在线视频| 久久性视频一级片| 999久久久精品免费观看国产| a在线观看视频网站| 日本成人三级电影网站| 亚洲精品久久国产高清桃花| 亚洲五月色婷婷综合| 久久天躁狠狠躁夜夜2o2o| 午夜激情av网站| 两个人视频免费观看高清| 欧美乱色亚洲激情| av视频在线观看入口| 亚洲va日本ⅴa欧美va伊人久久| 少妇裸体淫交视频免费看高清 | 一二三四在线观看免费中文在| 国产黄色小视频在线观看| 欧美乱色亚洲激情| 啦啦啦韩国在线观看视频| 宅男免费午夜| 精品国产乱子伦一区二区三区| 天堂动漫精品| 欧美最黄视频在线播放免费| 欧美日韩中文字幕国产精品一区二区三区| 麻豆一二三区av精品| 亚洲精品国产精品久久久不卡| 看免费av毛片| 美女高潮到喷水免费观看| 亚洲精品av麻豆狂野| 哪里可以看免费的av片| 国产成人精品久久二区二区免费| or卡值多少钱| 好看av亚洲va欧美ⅴa在| 欧美日韩亚洲综合一区二区三区_| 国产亚洲精品一区二区www| 丝袜美腿诱惑在线| 欧美激情久久久久久爽电影| 久久人妻福利社区极品人妻图片| 亚洲一区二区三区不卡视频| 中文字幕另类日韩欧美亚洲嫩草| 日韩精品中文字幕看吧| 人成视频在线观看免费观看| 在线观看免费视频日本深夜| 好男人在线观看高清免费视频 | 成年女人毛片免费观看观看9| 丁香欧美五月| 校园春色视频在线观看| 欧美午夜高清在线| 91老司机精品| 国产黄a三级三级三级人| 看片在线看免费视频| 美女高潮喷水抽搐中文字幕| 人人澡人人妻人| 老汉色av国产亚洲站长工具| 国产成人精品久久二区二区免费| 日本免费a在线| 欧美性长视频在线观看| 最近最新免费中文字幕在线| 久久久久久免费高清国产稀缺| 中国美女看黄片| 一边摸一边做爽爽视频免费| 久久久国产成人免费| 18禁观看日本| x7x7x7水蜜桃| 国产一区在线观看成人免费| 黄色丝袜av网址大全| 此物有八面人人有两片| 成人国语在线视频| 久久人人精品亚洲av| 国产v大片淫在线免费观看| 日韩精品免费视频一区二区三区| www.熟女人妻精品国产| 欧美午夜高清在线| 国产激情欧美一区二区| 叶爱在线成人免费视频播放| 老熟妇乱子伦视频在线观看| 中国美女看黄片| 在线观看www视频免费| 国产真实乱freesex| 久久欧美精品欧美久久欧美| 首页视频小说图片口味搜索| 婷婷六月久久综合丁香| 精品国产乱子伦一区二区三区| 美女高潮到喷水免费观看| 黑人操中国人逼视频| 国产伦人伦偷精品视频| 侵犯人妻中文字幕一二三四区| 国产精品 国内视频| 一区二区三区激情视频| 热99re8久久精品国产| 99国产极品粉嫩在线观看| 成年女人毛片免费观看观看9| 国产主播在线观看一区二区| 色尼玛亚洲综合影院| 久久精品国产亚洲av香蕉五月| 国产高清videossex| www.999成人在线观看| 亚洲中文日韩欧美视频| 国产精品精品国产色婷婷| 两性午夜刺激爽爽歪歪视频在线观看 | 欧美性长视频在线观看| 久久久久久久精品吃奶| 国产亚洲精品久久久久久毛片| 久久午夜亚洲精品久久| 亚洲精品一区av在线观看| 人妻久久中文字幕网| 91麻豆精品激情在线观看国产| 亚洲在线自拍视频| 精品久久久久久久毛片微露脸| 在线十欧美十亚洲十日本专区| 母亲3免费完整高清在线观看| 欧美三级亚洲精品| 黄色片一级片一级黄色片| 午夜影院日韩av| 亚洲男人的天堂狠狠| 亚洲av美国av| 热re99久久国产66热| 亚洲欧美一区二区三区黑人| 亚洲熟女毛片儿| 一进一出抽搐动态| 中文字幕人成人乱码亚洲影| 91大片在线观看| 国产精品日韩av在线免费观看| 国产熟女xx| 母亲3免费完整高清在线观看| 妹子高潮喷水视频| 国产亚洲欧美98| 一区二区日韩欧美中文字幕| 黄色a级毛片大全视频| 91成年电影在线观看| 久久久久国产精品人妻aⅴ院| 日韩大码丰满熟妇| 亚洲专区中文字幕在线| 精品日产1卡2卡| 香蕉av资源在线| 免费一级毛片在线播放高清视频| 国产亚洲精品久久久久5区| 麻豆成人av在线观看| av福利片在线| 成人三级黄色视频| 中文在线观看免费www的网站 | 亚洲人成电影免费在线| 亚洲国产毛片av蜜桃av| 两个人看的免费小视频| 久久久久久大精品| 欧美精品亚洲一区二区| 亚洲成人久久爱视频| 久久国产精品男人的天堂亚洲| 九色国产91popny在线| 露出奶头的视频| 老熟妇乱子伦视频在线观看| 午夜精品久久久久久毛片777| 日韩成人在线观看一区二区三区| 国产高清videossex| 国产亚洲欧美在线一区二区| 又黄又爽又免费观看的视频| 国产av又大| 久久久国产欧美日韩av| 九色国产91popny在线| or卡值多少钱| 又紧又爽又黄一区二区| 欧美人与性动交α欧美精品济南到| 波多野结衣av一区二区av| 国产aⅴ精品一区二区三区波| av有码第一页| 在线国产一区二区在线| 麻豆av在线久日| 嫩草影视91久久| 法律面前人人平等表现在哪些方面| 亚洲熟妇熟女久久| 日韩精品中文字幕看吧| 99精品久久久久人妻精品| 18禁黄网站禁片免费观看直播| 国产精品永久免费网站| 丰满的人妻完整版| 成人18禁高潮啪啪吃奶动态图| av中文乱码字幕在线| 男人舔女人的私密视频| 亚洲五月天丁香| 精品福利观看| 国产亚洲精品久久久久5区| 日日摸夜夜添夜夜添小说| 欧美日本亚洲视频在线播放| 精品熟女少妇八av免费久了| 在线观看免费视频日本深夜| 亚洲aⅴ乱码一区二区在线播放 | 精品国内亚洲2022精品成人| 欧美激情极品国产一区二区三区| 日日爽夜夜爽网站| 国产一区二区三区在线臀色熟女| 国产成人精品久久二区二区免费| 成在线人永久免费视频| 超碰成人久久| 母亲3免费完整高清在线观看| 一二三四在线观看免费中文在| 国产成人影院久久av| 欧美日韩黄片免| 人成视频在线观看免费观看| 欧美精品亚洲一区二区| www.www免费av| 久久精品91蜜桃| 在线观看日韩欧美| 久久久久亚洲av毛片大全| 国产亚洲欧美98| 中文字幕人成人乱码亚洲影| 午夜成年电影在线免费观看| 午夜福利免费观看在线| 久久香蕉精品热| 日韩一卡2卡3卡4卡2021年| 亚洲成人国产一区在线观看| 中文字幕精品免费在线观看视频| 真人一进一出gif抽搐免费| 国产一区在线观看成人免费| 成人欧美大片| 久久中文字幕一级| 757午夜福利合集在线观看| 久久久国产成人免费| 日韩精品免费视频一区二区三区| 亚洲av中文字字幕乱码综合 | 日韩免费av在线播放| 国产精品 欧美亚洲| 日韩视频一区二区在线观看| 女同久久另类99精品国产91| 在线av久久热| 免费高清视频大片| 亚洲色图av天堂| 亚洲人成伊人成综合网2020| 国产精品,欧美在线| 欧美乱色亚洲激情| 别揉我奶头~嗯~啊~动态视频| 欧美黑人精品巨大| 91在线观看av| 制服诱惑二区| 亚洲片人在线观看| 丁香六月欧美| 男女视频在线观看网站免费 | 老汉色av国产亚洲站长工具| 99国产精品一区二区三区| 麻豆国产av国片精品| 人人妻人人澡欧美一区二区| 国产在线精品亚洲第一网站| 国产高清有码在线观看视频 | 天天躁狠狠躁夜夜躁狠狠躁| 亚洲 欧美一区二区三区| 欧美日韩亚洲综合一区二区三区_| 国产精品国产高清国产av| 黄片大片在线免费观看| 亚洲精品国产区一区二| 欧美成狂野欧美在线观看| 麻豆成人av在线观看| 久久亚洲真实| 一级a爱视频在线免费观看| 99riav亚洲国产免费| 国产高清视频在线播放一区| 十分钟在线观看高清视频www| 日韩欧美一区视频在线观看| 男女那种视频在线观看| 亚洲熟妇熟女久久| 精品国内亚洲2022精品成人| 99国产综合亚洲精品| 午夜激情av网站| 国产色视频综合| ponron亚洲| 少妇的丰满在线观看| 亚洲片人在线观看| 亚洲欧美日韩高清在线视频| 精品久久久久久久毛片微露脸| 亚洲一区中文字幕在线| 午夜免费激情av| 日韩精品青青久久久久久| 成熟少妇高潮喷水视频| 岛国视频午夜一区免费看|