• <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蜜桃| 久久久久久久久中文| 欧美三级亚洲精品| 精品国产三级普通话版| 国精品久久久久久国模美| 国产精品国产三级国产专区5o| 亚洲在线自拍视频| 你懂的网址亚洲精品在线观看| 26uuu在线亚洲综合色| 国产视频首页在线观看| 免费在线观看成人毛片| 波多野结衣巨乳人妻| 大陆偷拍与自拍| 日韩欧美精品v在线| 精品一区在线观看国产| 人妻系列 视频| 可以在线观看毛片的网站| 免费观看在线日韩| 国产 亚洲一区二区三区 | 国产精品女同一区二区软件| 秋霞在线观看毛片| 尾随美女入室| av免费观看日本| av一本久久久久| 国产精品久久久久久精品电影小说 | 欧美三级亚洲精品| 中文精品一卡2卡3卡4更新| 欧美xxxx黑人xx丫x性爽| 亚洲国产成人一精品久久久| 国产黄片美女视频| 高清日韩中文字幕在线| 日本欧美国产在线视频| 亚洲丝袜综合中文字幕| 午夜日本视频在线| 青春草亚洲视频在线观看| 水蜜桃什么品种好| av在线亚洲专区| 亚洲av国产av综合av卡| 嫩草影院精品99| 最近2019中文字幕mv第一页| 国产在线一区二区三区精| 高清视频免费观看一区二区 | 如何舔出高潮| 成年版毛片免费区| 男女下面进入的视频免费午夜| av在线天堂中文字幕| 久久久色成人| 午夜激情欧美在线| 亚洲无线观看免费| 精品人妻视频免费看| av在线天堂中文字幕| 建设人人有责人人尽责人人享有的 | 亚洲精品国产成人久久av| 精品人妻熟女av久视频| 人妻一区二区av| 亚洲精品日韩在线中文字幕| av一本久久久久| 一级黄片播放器| 欧美高清成人免费视频www| 久久久亚洲精品成人影院| 亚洲av成人精品一二三区| 国产亚洲午夜精品一区二区久久 | 国产精品嫩草影院av在线观看| 男女那种视频在线观看| 午夜福利在线观看免费完整高清在| 99久久中文字幕三级久久日本| 国语对白做爰xxxⅹ性视频网站| 亚洲天堂国产精品一区在线| 国产女主播在线喷水免费视频网站 | 女的被弄到高潮叫床怎么办| .国产精品久久| 亚洲国产av新网站| 男人爽女人下面视频在线观看| 欧美97在线视频| 午夜福利网站1000一区二区三区| 听说在线观看完整版免费高清| 麻豆av噜噜一区二区三区| 91狼人影院| 街头女战士在线观看网站| 深爱激情五月婷婷| 极品教师在线视频| 久久久久久久久久人人人人人人| 2018国产大陆天天弄谢| 亚洲av免费高清在线观看| 夜夜爽夜夜爽视频| 久久99蜜桃精品久久| 嫩草影院新地址| 18禁动态无遮挡网站| 国产成人精品福利久久| 亚洲熟妇中文字幕五十中出| 三级国产精品片| 国产成人一区二区在线| 亚洲伊人久久精品综合| 成人av在线播放网站| 看黄色毛片网站| 汤姆久久久久久久影院中文字幕 | 秋霞伦理黄片| 男女边摸边吃奶| 听说在线观看完整版免费高清| 干丝袜人妻中文字幕| 日韩,欧美,国产一区二区三区| 亚洲精品日韩av片在线观看| 午夜老司机福利剧场| 亚洲在线观看片| 91久久精品电影网| 亚洲精品久久午夜乱码| 国产色爽女视频免费观看| 两个人的视频大全免费| 我要看日韩黄色一级片| 99热这里只有精品一区| 欧美日韩综合久久久久久| 亚洲欧美精品专区久久| 国产男人的电影天堂91| 天天躁日日操中文字幕| 秋霞在线观看毛片| 中文字幕av在线有码专区| 国产有黄有色有爽视频| 国产 一区精品| 亚洲18禁久久av| 91精品伊人久久大香线蕉| 国产精品一二三区在线看| 欧美一级a爱片免费观看看| 久久久久久久久久久免费av| 国产av在哪里看| 一边亲一边摸免费视频| 嫩草影院入口| 色视频www国产| 免费人成在线观看视频色| 色网站视频免费| 91av网一区二区| 有码 亚洲区| 免费观看的影片在线观看| 精品久久久久久成人av| 国产伦一二天堂av在线观看| 欧美bdsm另类| 精品人妻熟女av久视频| 十八禁网站网址无遮挡 | 亚洲av成人精品一区久久| 天天躁夜夜躁狠狠久久av| 青青草视频在线视频观看| 欧美成人精品欧美一级黄| 欧美最新免费一区二区三区| 国产精品一区二区三区四区久久| 久久人人爽人人片av| 黄色配什么色好看| 亚洲综合色惰| 国产精品嫩草影院av在线观看| kizo精华| 久久精品久久久久久噜噜老黄| 一区二区三区乱码不卡18| 可以在线观看毛片的网站| 亚洲自拍偷在线| 大片免费播放器 马上看| 国产精品麻豆人妻色哟哟久久 | 亚洲国产成人一精品久久久| 夜夜看夜夜爽夜夜摸| 男女国产视频网站| 国产精品久久久久久久久免| 在线a可以看的网站| 国产亚洲91精品色在线| 久久99热6这里只有精品| 网址你懂的国产日韩在线| 亚洲18禁久久av| 久久99热6这里只有精品| 人人妻人人澡人人爽人人夜夜 | 国产乱人偷精品视频| 一区二区三区高清视频在线| 偷拍熟女少妇极品色| 中文乱码字字幕精品一区二区三区 | 亚洲伊人久久精品综合| 国产精品一区www在线观看| 麻豆成人av视频| 久久99热这里只有精品18| 国产精品一区www在线观看| 内射极品少妇av片p| 七月丁香在线播放| 国产黄色视频一区二区在线观看| 国产成人freesex在线| 最近最新中文字幕免费大全7| 亚洲欧美一区二区三区国产| 久久鲁丝午夜福利片| 国产一区亚洲一区在线观看| 建设人人有责人人尽责人人享有的 | 国产精品久久久久久av不卡| 成年女人在线观看亚洲视频 | 两个人视频免费观看高清| 观看免费一级毛片| 国产极品天堂在线| 精品久久久久久久末码| 成人午夜精彩视频在线观看| 26uuu在线亚洲综合色| 好男人在线观看高清免费视频| 尾随美女入室| 高清日韩中文字幕在线| 日韩一区二区三区影片| 国产白丝娇喘喷水9色精品| 亚洲av免费在线观看| 99热这里只有是精品在线观看| 久久久久久久久久人人人人人人| 国产精品日韩av在线免费观看| 99久国产av精品| 国产一级毛片七仙女欲春2| 午夜亚洲福利在线播放| 国产一区有黄有色的免费视频 | 国产美女午夜福利| 黄色一级大片看看| 少妇的逼好多水| 午夜福利成人在线免费观看| 1000部很黄的大片| 日本av手机在线免费观看| 日本一本二区三区精品| 国产精品无大码| 99久国产av精品国产电影| 欧美人与善性xxx| 亚洲精品色激情综合| 超碰97精品在线观看| 亚洲第一区二区三区不卡| 午夜亚洲福利在线播放| 一级爰片在线观看| 免费观看a级毛片全部| 日韩亚洲欧美综合| 日本熟妇午夜| 2018国产大陆天天弄谢| 精品99又大又爽又粗少妇毛片| av网站免费在线观看视频 | 91狼人影院| 亚洲av免费在线观看| 国产精品伦人一区二区| 1000部很黄的大片| 色播亚洲综合网| 欧美潮喷喷水| 九九久久精品国产亚洲av麻豆| 亚洲精品成人av观看孕妇| 在线a可以看的网站| 99九九线精品视频在线观看视频| 十八禁国产超污无遮挡网站| 亚洲欧美精品专区久久| 亚洲最大成人手机在线| 高清视频免费观看一区二区 | 欧美激情国产日韩精品一区| 亚洲自偷自拍三级| 欧美 日韩 精品 国产| 成人毛片a级毛片在线播放| 一二三四中文在线观看免费高清| 成人午夜精彩视频在线观看| 久久久久久久久大av| 啦啦啦啦在线视频资源| 两个人的视频大全免费| 男女视频在线观看网站免费| 欧美变态另类bdsm刘玥| 一级毛片 在线播放| 成人二区视频| 一区二区三区高清视频在线| 好男人视频免费观看在线| 日韩成人伦理影院| 国产精品一及| 国产69精品久久久久777片| 寂寞人妻少妇视频99o| 欧美zozozo另类| 欧美高清成人免费视频www| 亚洲精品色激情综合| 中文字幕免费在线视频6| 午夜福利在线在线| 亚洲四区av| 久久久久久久久久久免费av| 美女xxoo啪啪120秒动态图| 99热这里只有是精品50| 亚洲精品成人久久久久久| 国产黄色视频一区二区在线观看| 九九在线视频观看精品| 91在线精品国自产拍蜜月| 成人午夜精彩视频在线观看| 日韩欧美三级三区| 欧美极品一区二区三区四区| 亚洲av中文av极速乱| 午夜福利视频1000在线观看| 国产白丝娇喘喷水9色精品| 中国国产av一级| 国产成人精品一,二区| 久久这里只有精品中国| 色视频www国产| ponron亚洲| 婷婷色综合www| 黄色一级大片看看| 成人漫画全彩无遮挡| 亚洲三级黄色毛片| 菩萨蛮人人尽说江南好唐韦庄| 欧美3d第一页| 日韩人妻高清精品专区| 精品久久久精品久久久| 国产毛片a区久久久久| 我要看日韩黄色一级片| 国产淫语在线视频| 高清日韩中文字幕在线| 欧美极品一区二区三区四区| 亚洲乱码一区二区免费版| 久久久久网色| 国产有黄有色有爽视频| 久久韩国三级中文字幕| 人妻一区二区av| 色哟哟·www| 亚洲av成人精品一区久久| 99久久精品热视频| 亚洲av电影不卡..在线观看| 国产黄频视频在线观看| 水蜜桃什么品种好| 大话2 男鬼变身卡| 一级毛片我不卡| 白带黄色成豆腐渣| 综合色丁香网| 免费av观看视频| 精品久久久久久久久av| 男女那种视频在线观看| 国产精品久久久久久精品电影| 亚洲av电影不卡..在线观看| 国产精品无大码| 久久人人爽人人片av| 国产精品综合久久久久久久免费| 99久国产av精品国产电影| 欧美日韩综合久久久久久| 亚洲国产精品成人综合色| 熟妇人妻久久中文字幕3abv| 听说在线观看完整版免费高清| 成人午夜精彩视频在线观看| 2018国产大陆天天弄谢| 国产高清三级在线| 国产精品一及| 亚洲国产精品国产精品| 热99在线观看视频| 亚洲欧美精品专区久久| 乱系列少妇在线播放| 欧美三级亚洲精品| 丰满少妇做爰视频| 人妻制服诱惑在线中文字幕| 日本免费a在线| 成人综合一区亚洲| 深夜a级毛片| 日韩 亚洲 欧美在线| 亚洲第一区二区三区不卡| 乱人视频在线观看| 一级爰片在线观看| 爱豆传媒免费全集在线观看| 国产在视频线在精品| 婷婷六月久久综合丁香| 女的被弄到高潮叫床怎么办| 久久国产乱子免费精品| 国产精品一二三区在线看| 2022亚洲国产成人精品| 国产成人aa在线观看| 男的添女的下面高潮视频| 婷婷色综合www| 国产av码专区亚洲av| 能在线免费看毛片的网站| 亚洲av二区三区四区| 色播亚洲综合网| 国产成人精品福利久久| 国精品久久久久久国模美| 亚洲最大成人中文| 校园人妻丝袜中文字幕| 国产探花极品一区二区| 偷拍熟女少妇极品色| av在线亚洲专区| 免费看美女性在线毛片视频| 看黄色毛片网站| 最新中文字幕久久久久| 精品一区二区三区视频在线| 精品少妇黑人巨大在线播放| 亚洲欧美成人综合另类久久久| 成人性生交大片免费视频hd| 免费看光身美女| 久久久精品94久久精品| 久久久久久久大尺度免费视频| 亚洲电影在线观看av| 日本av手机在线免费观看| 美女被艹到高潮喷水动态| eeuss影院久久| 国产精品一区www在线观看| 国产探花极品一区二区| 国产男人的电影天堂91| 亚洲av不卡在线观看| ponron亚洲| 大片免费播放器 马上看| 亚洲欧美成人精品一区二区| 九九久久精品国产亚洲av麻豆| 少妇熟女欧美另类| 高清午夜精品一区二区三区| 少妇裸体淫交视频免费看高清| 亚洲在线观看片| 3wmmmm亚洲av在线观看| 最近2019中文字幕mv第一页| 欧美最新免费一区二区三区| 乱人视频在线观看| 欧美一级a爱片免费观看看| 亚洲av成人精品一二三区| 91精品国产九色| 神马国产精品三级电影在线观看| 男人和女人高潮做爰伦理| 永久免费av网站大全| 人妻少妇偷人精品九色| 天天躁夜夜躁狠狠久久av| 三级毛片av免费| 免费看a级黄色片| 成人二区视频| 男女下面进入的视频免费午夜| 日本三级黄在线观看| www.色视频.com| 十八禁国产超污无遮挡网站| 久久久国产一区二区| 午夜视频国产福利| 一级二级三级毛片免费看| 国产精品99久久久久久久久| 欧美激情在线99| 国产精品1区2区在线观看.| 最后的刺客免费高清国语| 国产午夜福利久久久久久| 中文乱码字字幕精品一区二区三区 | 久久久精品94久久精品| 国产日韩欧美在线精品| 91av网一区二区| 精品人妻熟女av久视频| 少妇高潮的动态图| 舔av片在线| 国产色婷婷99| 久久99热这里只有精品18| 亚洲精品中文字幕在线视频 | 男的添女的下面高潮视频| 国语对白做爰xxxⅹ性视频网站| 黄色配什么色好看| 欧美性猛交╳xxx乱大交人| 大又大粗又爽又黄少妇毛片口| 99热网站在线观看| 亚洲精品中文字幕在线视频 | 亚洲乱码一区二区免费版| 日韩 亚洲 欧美在线| 国产熟女欧美一区二区| 卡戴珊不雅视频在线播放| 国产精品不卡视频一区二区| xxx大片免费视频| 国产综合懂色| 久久久久精品久久久久真实原创| 干丝袜人妻中文字幕| 神马国产精品三级电影在线观看| 日韩亚洲欧美综合| 国产精品一区二区在线观看99 | 女人被狂操c到高潮| 亚洲精品乱码久久久久久按摩| 一本一本综合久久| 成人亚洲精品一区在线观看 | 国产成人午夜福利电影在线观看| 777米奇影视久久| 天堂影院成人在线观看| eeuss影院久久| 国产单亲对白刺激| 只有这里有精品99| 亚洲人与动物交配视频| 成人毛片60女人毛片免费| 亚洲精品国产av成人精品| 波野结衣二区三区在线| 日日摸夜夜添夜夜爱| 看免费成人av毛片| 18禁在线播放成人免费| 精品酒店卫生间| 99热这里只有精品一区| 永久网站在线| 国产精品国产三级国产专区5o| 免费观看a级毛片全部| 又爽又黄a免费视频| 2018国产大陆天天弄谢| 亚洲在久久综合| av国产久精品久网站免费入址| 身体一侧抽搐| 国产亚洲一区二区精品| 日韩制服骚丝袜av| 国产伦精品一区二区三区四那| 综合色丁香网| 国产av码专区亚洲av| 日韩 亚洲 欧美在线| 亚洲成人中文字幕在线播放| 男女那种视频在线观看| 乱人视频在线观看| 老司机影院毛片| 国产人妻一区二区三区在| 精品人妻视频免费看| 老师上课跳d突然被开到最大视频| 亚洲欧洲日产国产| 国产探花极品一区二区| 高清欧美精品videossex| 神马国产精品三级电影在线观看| 蜜桃久久精品国产亚洲av| 国产 亚洲一区二区三区 | 99久久人妻综合| 成人国产麻豆网| 久久精品国产鲁丝片午夜精品| 精品久久久久久电影网| 2018国产大陆天天弄谢| 国产乱人偷精品视频| 少妇熟女aⅴ在线视频| 久久人人爽人人爽人人片va| 亚洲欧美一区二区三区国产| av在线蜜桃| 欧美xxxx性猛交bbbb| 亚洲精品日韩在线中文字幕| 超碰97精品在线观看| 欧美日韩精品成人综合77777| 亚洲精品中文字幕在线视频 | 日韩视频在线欧美| 淫秽高清视频在线观看| 色5月婷婷丁香| 国产成人a∨麻豆精品| 你懂的网址亚洲精品在线观看| 成年人午夜在线观看视频 | 久久久精品94久久精品| 日韩成人伦理影院| 最近中文字幕高清免费大全6| 熟妇人妻久久中文字幕3abv| 69人妻影院| 色综合色国产| 午夜老司机福利剧场| 日本黄大片高清| 99热全是精品| 久久精品国产自在天天线| 精品人妻一区二区三区麻豆| 国产精品人妻久久久影院| 波多野结衣巨乳人妻| 综合色丁香网| 三级男女做爰猛烈吃奶摸视频| 久久久色成人| 婷婷色麻豆天堂久久| 嘟嘟电影网在线观看| 国产伦精品一区二区三区视频9| 国产成人a区在线观看| 国产成人福利小说| 久久亚洲国产成人精品v| 久久精品久久精品一区二区三区| 网址你懂的国产日韩在线| 国产在视频线精品| 纵有疾风起免费观看全集完整版 | 中文资源天堂在线| 嫩草影院新地址| 亚洲成人av在线免费| 在现免费观看毛片| 国产黄色免费在线视频| 99久国产av精品| 国产 亚洲一区二区三区 | 国产日韩欧美在线精品| 国产亚洲91精品色在线| 联通29元200g的流量卡| 久久久久久国产a免费观看| 91狼人影院| 午夜福利视频精品| 国产精品久久久久久久电影| av女优亚洲男人天堂| 中文字幕av成人在线电影| 久久久久性生活片| 欧美日韩在线观看h| 亚洲不卡免费看| 亚洲18禁久久av| 欧美日本视频| 国产免费视频播放在线视频 | 男女视频在线观看网站免费| 日本免费在线观看一区| av女优亚洲男人天堂| 蜜臀久久99精品久久宅男| 国产91av在线免费观看| 99久久精品热视频| 99久久精品一区二区三区| 欧美性感艳星| 一个人看的www免费观看视频| 午夜精品在线福利| 午夜日本视频在线| 免费看日本二区| av免费在线看不卡| 成人综合一区亚洲| 99九九线精品视频在线观看视频| 白带黄色成豆腐渣| 国内精品宾馆在线| 2022亚洲国产成人精品| 成年av动漫网址| 亚洲精品色激情综合| 人人妻人人看人人澡| 色播亚洲综合网| 我的老师免费观看完整版| 久久99热6这里只有精品| 亚洲怡红院男人天堂| 欧美一级a爱片免费观看看| 草草在线视频免费看| av黄色大香蕉| 麻豆国产97在线/欧美| 久久久欧美国产精品| 69人妻影院| 亚洲欧美精品自产自拍| 亚洲人与动物交配视频| 久久精品久久久久久噜噜老黄| 欧美一级a爱片免费观看看| 777米奇影视久久| 中文在线观看免费www的网站| 国产伦精品一区二区三区视频9| 夜夜爽夜夜爽视频| 成人特级av手机在线观看| 欧美人与善性xxx| 天美传媒精品一区二区| 欧美一级a爱片免费观看看| 精品一区在线观看国产| 一个人看的www免费观看视频| 能在线免费观看的黄片| 亚洲自拍偷在线| 亚洲精品456在线播放app| 日韩欧美精品v在线| 亚洲av电影不卡..在线观看| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 |