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

    Energy properties and spatial plume profile of ionic liquid ion sources based on an array of porous metal strips

    2022-01-10 14:51:38XinyuLIU劉欣宇XiaomingKANG康小明HanwenDENG鄧涵文andYimingSUN孫逸鳴
    Plasma Science and Technology 2021年12期

    Xinyu LIU (劉欣宇), Xiaoming KANG (康小明),2, Hanwen DENG (鄧涵文) and Yiming SUN (孫逸鳴)

    1State Key Laboratory of Mechanical System and Vibration,School of Mechanical Engineering,Shanghai Jiao Tong University, Shanghai 200240, People’s Republic of China

    2 Shanghai Frontier Science Center for Gravitational Wave Detection, School of Mechanical Engineering,Shanghai Jiao Tong University, Shanghai 200240, People’s Republic of China

    Abstract An ionic liquid ion source (ILIS) is a kind of high brightness ion source capable of providing high-speed positive or negative ion beams.This paper presents a miniaturized ILIS based on an array of porous metal strips.The porous emitter array,integrated with seven 10 mm long strips,is fabricated using wire electrical discharge machining(WEDM)combined with electrochemical etching.The assembled ILIS is 30 mm×30 mm×17.5 mm in size and weighs less than 25 g.A series of experiments,including an I–V characteristic test,a retarding potential analyzer(RPA)test, and a spatial plume distribution test, have been conducted in vacuo to characterize the performance of the ILIS.Results show that the emitted current is up to about 800 μA and ion transparency is as high as 94%.Besides, RPA curves reveal that the total fragmentation rate of the emitted particles accounts for 48.8%in positive mode and 59.8%in negative mode.Further,with the increase in applied acceleration voltage, the voltage loss rises while the energy efficiency decreases.It is also found that the plume perpendicular to the strips has a higher divergence than the one parallel to the strips.A numerical simulation by COMSOL reveals that the electric field distribution between the two electrodes results in such a spatial plume profile.

    Keywords:electrospray,ionic liquid ion source,retarding potential analyzer,energy distribution,spatial plume profile

    1.Introduction

    An ionic liquid ion source (ILIS) is a promising ion source based on the electrospray phenomenon.Conductive fluids subjected to strong electric fields deform from a rounded shape into a conical structure and emit charged droplets,ions,or a mixture of both [1], known as electrospray.Traditional electrospray sources emit combinations of ions and charged droplets from organic solvents [2–4], such as doped glycerol and formamide.However, ILISs produce purely positive or negative ion beams from ionic liquids (ILs).ILs are molten salts at room temperature.Unlike those solvents, ILs are composed of cations and anions, exhibiting relatively high electrical conductivities,low volatility as well as high thermal stability[5].ILs were first applied in electrospray by Romero Sanzet alin 2003 [6], and it was discovered that pure ion emission could be achieved with EMI-BF4at low flow rates.Apart from EMI-BF4, EMI-Im and EMI-(HF)2.3F have been reported to reach the purely ionic regime (PIR) [7, 8].

    ILISs can be applied to a wide range of areas, from micro-fabrication to space propulsion[9–12].There are many benefits of using ILISs as electric thrusters.Firstly, the emitted ions in ILISs are accelerated to a high velocity of the order of magnitude of 105m s?1, contributing to a high specific impulse.Besides, no gas discharge is required to produce charged particles compared with plasma sources[13–15].The liquid propellant can be passively fed by capillary force.The removal of additional pumps, valves and pipelines makes the whole system light and compact.Furthermore, compared with liquid metal ion sources (LMISs),ILISs need no heaters and operate at a lower voltage[16,17],resulting in a considerable reduction in power consumption.

    An ILIS device basically consists of an emitter, an extractor, a storage tank, and housing.A sufficiently high voltage difference is applied between the emitter and the downstream extractor so that a strong electric field is formed in the vicinity of the emitter tip.Ions are extracted from the electrified menisci formed on the emitter tip and accelerated by the electrostatic field.Ion emission from the menisci of ionic liquid formedin vacuois taken to obey the kinetic law[18]

    wherejeis the current emitted per unit area of the meniscus surface,σis the density of surface charge,kBis the Boltzmann’s constant,Tis the liquid temperature,his the Planck’s constant,DGis the activation energy barrier for solvated species,is the normal component of the electric field acting at the meniscus exterior on the vacuum side, andG (Env) is the reduction ofDGdue to the presence of the electric field.G( Env) is assumed to take the formin whichε0is the permittivity of the vacuum andqis the charge of the emitted ions.Significant fluxes can only begin whenG (Env) approachesDG[19].The solvation energy for ionic liquid is about 1–2 eV, thusEnvis estimated to be 0.1–1 V nm?1.

    To amplify the electric field intensity,emitters are usually terminated by sharp geometries with an apex curvature radius of several to tens of microns.Further, to scale up the total emitted current, emitters are fabricated into linear or planar arrays [20–22].There are generally three types of emitters:capillary emitters, externally wetted emitters, and porous emitters, as shown in figure 1.Capillaries have been used as emitters since the 1960s.However,micro-fabricated capillary arrays would be unable to produce stable emissions from all capillaries since the hydraulic impedance is not high enough.Though increasing the capillary aspect ratio seems to be an effective way to solve this problem, the method suffers from larger operating voltages and capillary clogging.Another drawback of capillary emitters is the addition of a liquid supply system, increasing the volume and weight of ILIS.As to externally wetted emitters,clogging can be avoided using such an open architecture.Another advantage of externally wetted emitters is that ionic liquid can be passively fed without any pumps or valves.But challenges lie in the improvement of the wettability of emitters and the sustainability of liquid supply.

    Among the three types of emitter configurations, porous emitters are adopted by an increasing number of ILIS applications.In a porous ILIS, inner pore structures provide multiple paths to arrive at the tips,and liquid can be passively fed via capillary force.Porous metals are important materials for ILIS emitters.Leggeet alfirstly proposed linear emitter arrays fabricated from porous tungsten by electrochemical etching [22].Porous emitters were shown to generate a wide range of emitted currents up to levels generally greater than those in externally wetted emitters.Further, a planar array of 480 individual emitter tips was achieved on a 1 cm2porous nickel substrate by Courtneyet alto improve the integration of the emitters[21].In the subsequent studies,the substrate of the emitters was expanded to porous glass and carbon xerogels [23, 24], and laser ablation was used to fabricate the emitter array.Except for sharp tips, emitters are also fabricated into an array of strips terminated with triangular prisms.Courtneyet almanufactured a porous glass ILIS with strip emitters using conventional mill[20].Chenet aldeveloped a porous ceramic ILIS with strip emitters by sintering [25].However, porous metal strip emitters are difficult to machine by electrochemical micromachining(suffering from geometry nonuniformity)or laser ablation(problems in pore blocking).Thus, few investigations of ILISs with an array of porous metal strips have been reported in the literature before.

    This paper presents a miniaturized ILIS based on an array of porous stainless steel strips.Section 2 gives the description of the processes of the emitter preparation and ILIS fabrication.To characterize the developed ILIS,anI–Vcharacteristic test, a retarding potential analyzer (RPA) test, and a spatial plume distribution test are conducted.Section 3 describes the experimental apparatus of the tests, and the results are revealed and discussed in section 4.

    2.Emitter preparation and ILIS fabrication

    The schematic diagram of the designed ILIS is depicted in figure 2.The ILIS includes a porous emitter array, an extractor, a storage tank, a distal electrode and housing.The emitter is a seven-strip array integrated on a porous stainless steel substrate.The trigonous strip is 1 mm tall and the stepped substrate has a total height of 3 mm.The adjacent strips are spaced by 2 mm, and each of them is 10 mm long.Ionic liquid EMI-BF4(1-ethyl-3-methylimidazolium tetrafluoroborate, and the molecular formula is C6H11BF4N2) is used in this paper.EMI-BF4is passively fed from the porous reservoir to the emitter by capillary force.The high voltage power supply represents either a positive or a negative high voltage source.When operating the ILIS with a positive bias,positive particles are emitted.On the contrary,when operated in negative polarity, negative particles are emitted.

    The emitter array is prepared using wire electrical discharge machining (WEDM) combined with electrochemical etching, as proposed in [26].The fabricated emitter array is observed by KEYENCE VHX-5000 digital microscope with a large depth of field and accurate measurement capabilities.The apex curvature radii of the strips range from 6 μm to 10 μm,as shown in figure 3,and the surface roughness along the apex of the emitter strip is around 10 μm.

    The extractor is laser cut from a 120 μm thick stainless steel sheet.Seven linear extractor holes are 0.6 mm in width with a pitch distance of 2 mm.The storage tank is integrated with housing which is computerized numerical control(CNC)machined from a polyetheretherketone (PEEK) bulk.After being wettedin vacuo,the emitter array is assembled with the extractor and the housing through screw fastening.The emitter and the extractor hole are aligned with the aid of the microscope, as shown in figure 4(a).Actually, the geometry between the emitter and the extractor is difficult to observe using a normal optics microscope.At hundreds of times magnification, the microscope can only focus on one plane,while the areas on other planes are blurry.Since the emitter apex and the extractor are not in the same plane, it is very hard to capture the geometry between the two electrodes.To solve the problem, the super deep scene digital microscope(KEYENCE VHX-5000) is used, which can capture the geometries at hundreds of separate focusing planes by automatically lifting the lens,and then reconstruct a 3D geometry by its internal image processing technology.In the 3D reconstruction process, the entity margins can be well captured, while the void margins are distorted to some extent.Thus,the scan geometry looks like a‘W’shape in figure 4(a).However, what we are concerned about, such as the distance between the two electrodes and the alignment of the emitter,can be successfully achieved by the 3D geometry.The points 1,2 and 3 are measurement points, taken from the two edges of the extractor and the emitter apex, respectively.The horizontal distances between the emitter apex and the two edges of the extractor are 295.82 μm and 303.50 μm, showing a good symmetry.The distance between the emitter apex and the upper surface of the extractor is 262.26 μm,meaning that the distance from the emitter apex to the lower surface of the extractor is around 140 μm.The assembled thruster measures 30 mm×30 mm×17.5 mm with a wet weight of less than 25 g, as shown in figure 4(b).

    3.Experimental apparatus

    The experiments are conductedin vacuoin the State Key Laboratory of Mechanical System and Vibration in Shanghai Jiao Tong University.The experimental facilities and instruments are shown in figure 5, including a vacuum chamber,two high voltage power suppliers, a digital multimeter, a signal generator as well as a high voltage amplifier.The vacuum chamber is 0.6 m in diameter and 1 m in length with two mechanical pumps, a roots pump, and two turbo molecular pumps, which can achieve a base pressure of 5×10?4Pa.The positive high voltage power supplier (Spellman SL30P30) and the negative one (Boher LAS-20kV-200N) provide acceleration voltages for the ILIS.The digital multimeter(Keithley DMM 7510) is used to measure and record the voltages and currents.Besides,the signal generator(Tektronix AFG 3022B), together with the high voltage amplifier (Matsusada AMPS-20B20),can provide sweeping voltage for the RPA test.

    3.1.I–V characteristic test

    Either the positive or negative high voltage supply is connected to the distal electrode,and the extractor is grounded.In the positive mode, the voltage is regulated from 0 to+3000 V, and in the negative mode, the voltage is regulated from 0 to ?3060 V.Both the emitted current, which is the current drawn by the emitter from the high voltage supply,and the intercepted current, which is defined as the current measured from the extractor to ground, are obtained by measuring the voltage drop across a 1 kΩ resistor.

    3.2.Retarding potential analysis test

    The RPA test is performed to measure the energy distribution of the emitted particles.The actual acceleration voltage is usually smaller than the acceleration voltage applied between the distal electrode and the extractor from the high-voltage power supply.Thus, the RPA test can obtain the nonkinetic voltage losses.Further,ion fragmentation,which is the breakup of the emitted ion clusters, occurs during the ionic liquid electrospray process, and RPA test can measure the rates of fragmentation of ion clusters by analyzing the RPA curves.

    As suggested in [27], a four-grid design is adopted in the RPA.The RPA device consists of a four-layer grid assembly,a metal shielding, and a collector, as depicted in figure 6.The first layer of the grid is grounded,and both the second and the fourth layers of the grids are connected to a voltage of ?40 V to suppress the secondary electron emission.The third grid is the retarding grid,connected to the sweeping voltage generated from a signal generator and a high voltage amplifier.The RPA device is installed 95 mm away from the source.The center of the ILIS and the center of the RPA aperture are mounted coaxially.The flying particles enter the RPA device through the front aperture, pass the four layers of grids, and finally arrive at the collector.The collected currents are measured and recorded by a Keithley DMM 7510 digital multimeter.

    3.3.Spatial plume distribution measurement

    The spatial plume distribution is tested by measuring the collected current when rotating the ILISs either parallel to the strips or perpendicular to the strips,as shown in figure 7.The ILIS is installed on a rotating platform, and the rotating shaft traverses the emission plane.The rotary motion of the platform is driven by a turbine worm mechanism connected to a stepper motor, and the rotational speed can be regulated by the frequency of the pulse signal input.The rotational angular speed is set to six degrees per second in the test.

    The RPA collector is used to collect the emitted current,mounted 155 mm downstream from the source.Unlike the RPA test,all four layers of grids,including the retarding grid,are grounded.A similar method using the RPA device to collect the ion plume with the sweeping grid grounded was implemented in [23].The ILIS center, aperture center, and collector center are installed coaxially.By rotating the ILIS,the plume with a certain divergence angle can selectively pass through the RPA aperture.Thus, the spatial plume distribution can then be measured from the collected current.

    4.Results and discussion

    4.1.I–V characteristics

    The variation of both emitted current and the intercepted current with the applied acceleration voltage is shown in figure 8.The startup voltageVstart,which is defined as the minimum voltage that starts an electrospray emission, is 1280 V.Apart from the experimental value, the startup voltage can also be estimated by the following equation [28]

    whereγis the surface tension coefficient,Rtis the apex curvature radius of the emitter,rpis the mean pore radius of the emitter,anddis the distance between two electrodes.Given that the apex curvature radii and pore sizes fluctuate in a range, an accurate estimation of the startup voltage is not expected.Here,we take some typical values in the range to calculate the startup voltage.Takingγ=45.2 mN m?1,d=140 μm,rp=3 μm,andRt=9 μm,the calculated startup voltage is 1217 V,with an error of 4.9% compared with the experimental value.

    Both the emitted current and the intercepted current grow with the increase of the applied acceleration voltage.In the positive polarity, the curve is similar to quadratic from the startup voltage to 2500 V while it is similar to linear form 2500 V to 3000 V.The emitted current reaches 804 μA at+3000 V, and ?790 μA at ?3060 V, and the intercepted current is within 6% of the emitted current, meaning that the ion transparency is up to 94%.

    4.2.Fragmentation of the ILIS

    The normalized RPA curves for both positive and negative modes are plotted in figure 9.As the sweeping voltage of the retarding grid increases from zero to the acceleration voltage,the normalized collected current first holds steady, then experiences three stages of decline, and finally decreases to zero in both polarities.The overall RPA curves are generally consistent with some other RPA studies of ILISs[25,29,30].According to the working principles of the ILISs,the emitted charged particles gain kinetic energy through electrostatic acceleration between emitter and extractor.If all the particles have the same kinetic energy, the collected current will maintain a steady value and dramatically decrease at a specific voltage.Therefore,the RPA curves indicate that the plume is not monoenergetic.

    The divergence of the energy distribution may be derived from ion extraction or ion fragmentation.Our previous time of flight test showed that the plume is mainly composed of monomers [31] (single ions, EMI+or-BF4) and dimers(single ions combined with a molecule, (EMI-BF4)EMI+or(EMI-BF4)-BF4), very little, if any, trimer (single ions combined with two molecules, (EMI-BF4)2EMI+or (EMI-BF4)2BF-4) or droplets.The relative molecular mass, molecular mass, charge, and charge to mass ratio of the ion species are listed in table 1.The inexistence of droplets contributes to concentrated extraction voltages.The energy of ions to be extracted from the ionic liquid is small (1–2 eV), which means there is little possibility of such significant energy dispersion (hundreds of eV) caused by ion extraction.Thus,only ion fragmentation effects need further consideration.

    Table 1.Relative molecular mass, molecular mass, charge, and charge-to-mass ratio of the ion species.

    If the fragmentation occurs in the field-free region, as shown in figure 10, which means the acceleration process is completed, according to the principle of conservation of energy

    wherem0is the mass of a molecule,m1is the mass of a monomer ion,v1is the velocity after electrostatic acceleration,q1is the charge of a monomer ion, andVais the applied acceleration voltage.

    Supposing that the velocity of the neutral molecule and the monomer ion remains unchanged after fragmentation, the accelerated ion travels at an invariable speed until it is stopped by the retarding grid

    whereVRPAis the potential of the retarding grid.The ratio ofVRPAandVacan be obtained from equations (3) and (4)

    The ratio is a constant only depending on the masses of the monomer ion and the molecule.For EMI-BF4,VRPA/Vais 0.36 in positive mode and 0.30 in negative mode.The segment indicating fragmentation in the field-free region is marked in red in the RPA curve in figure 10.

    If the dissociation of ion clusters takes place in the acceleration region, as shown in figure 11, both the molecule and the monomer ion are accelerated until they get to the departure position

    wherev1,1is the velocity of the molecule and the monomer ion, andV1is the electric potential of the departure position.

    After departure, the ion sustains acceleration up to leaving the acceleration region

    wherev1is the velocity of the monomer ion after ending the acceleration process.The accelerated ion travels at an invariable velocity until it is stopped by the retarding grid

    The ratio ofVRPAandVacan be derived from equations (6)–(8)

    The value ofV1ranges from 0 toVa,thusVRPA/Vais within the limits ofIn figure 10, the segment representing fragmentation in the acceleration region is marked in green in the RPA curve.

    If no fragmentation occurs, the theoretical value ofVRPA/Vais 1.The blue drop lines of the RPA curve in figures 10 and 11 denote the fraction of the ions without fragmentation.

    Through the above analysis, the three segments of decline in the RPA curve imply different types of fragmentation.Since all the ion species are singly charged, the fragmentation rate at a certain applied acceleration voltage can be valued as the current fraction.From figure 9, there are no strong signs that ion fragmentation is affected by the acceleration voltage in the same polarity.The average fragmentation rates are calculated and presented in figure 12.The total fragmentation rate of the emitted particles accounts for 48.8%in positive mode and 59.8% in negative mode.Coles and Miller proposed and validated that ion clusters composed of complex ions tend to be more stable [29, 32].The better stability of (EMI-BF4)EMI+than (EMI-BF4)-BF4is probably due to a more complex structure of(EMI-BF4)EMI+despite different polarities.

    4.3.Energy distribution and voltage losses of the ILIS

    Energy spectra of the ILIS can be derived from the derivative of the RPA curve.Figure 13 shows a series of normalized energy spectra of the ILIS at different acceleration voltages.It is found that there are two prominent peaks in each separate curve.This is due to the two sharp drop lines in the RPA curve, caused by fractions fragmented in the field-free region and fractions without fragmentation, respectively.After the calculation of the derivative, the signal intensity is therefore relatively high in the vicinity of these two positions, indicating that the plume contains two concentrated energy distributions.

    Figure 1.Three types of electrospray emitters: (a) capillary, (b) externally wetted, and (c) porous.

    Figure 2.Schematic diagram of the designed ILIS.

    Figure 3.Observation of the emitter under the digital microscope: (a) overall geometry and (b) measurement of the radii of apex curvature(rainbow bar denotes the relative height of the strip, μm).

    Figure 4.Assembly of the ionic liquid electrospray thruster:(a)alignment under microscope(rainbow bar denotes the relative height of the scanned geometry) and (b) the assembled thruster.

    Figure 5.Experimental facilities and instruments used in the tests.

    Figure 6.(a) Schematic diagram of the RPA setup and (b) RPA collector.

    Figure 7.Schematic diagram of the spatial plume distribution setup:(a)rotating parallel to the emission strips and(b)rotating perpendicular to the emission strips.

    Figure 8.I–V characteristic curve of the ILIS.

    Figure 9.Normalized RPA curves: (a) positive mode and (b) negative mode.

    Figure 10.The segment marked in red represents the fractions fragmented in the field-free region.

    Figure 11.The segment marked in green represents the fractions fragmented in the acceleration region.

    The voltage of the second peak,which indicates particles without fragmentation, is slightly smaller than the corresponding applied acceleration voltage.Since the unbroken particles obtained kinetic energy only by electrostatic field acceleration, the voltage loss of the ILIS can be obtained by analyzing the difference between the voltage of the second peak and the applied acceleration voltage.Further, energy efficiency, defined as the ratio of the actual acceleration voltage and the applied acceleration voltage, can also be calculated.Figure 14 illustrates voltage losses and energy efficiencies of the ILIS in positive and negative modes.Results indicate that the voltage loss is increased with the applied acceleration voltage.Correspondingly, the energy efficiency versus the applied acceleration voltage is decreased, though the decrease is relatively slight due to the simultaneous increase of the applied acceleration voltage.A similar trend was mentioned in references [30], though the material of the emitter is different.The voltage loss is probably due to the resistance of the propellant liquid or circuit.

    Figure 12.Fragmentation rates in (a) positive mode and (b) negative mode.

    Figure 13.Energy spectra of the ILIS in (a) positive mode and(b) negative mode.

    Figure 14.Voltage losses and energy efficiencies of the ILIS in positive and negative modes.

    Figure 15.Spatial plume profiles:(a)perpendicular to the strips in positive mode,(b)parallel to the strips in positive mode,(c)perpendicular to the strips in negative mode, and (d) parallel to the strips in positive modes.

    Figure 16.3D COMSOL modelling: (a) model establishment, (b) boundary conditions and (c) meshes.

    Figure 17.Simulated particle trajectories at t=60 ns.(a) Va=3000 V, perpendicular to the strips, (b) Va=3000 V, parallel to the strips,(c) Va=?3000 V, perpendicular to the strips, and (d) Va=?3000 V, parallel to the strips (right rainbow bar denotes the velocity of the particles, m s?1).

    4.4.Spatial plume profile of the ILIS

    The spatial plume profiles both perpendicular and parallel to the strips of the ILIS are shown in figure 15.The collected current has been normalized to the height of the peak at an applied acceleration voltage of 3000 V or ?3000 V.The main source of error is approximating the finite dimensions of the ILIS as a point source,which is estimated to be less than six degrees.

    It is found that the plume cross profile is approximately symmetric in the direction of either perpendicular or parallel to the strips, and better symmetry is found at higher acceleration voltages.This phenomenon is likely caused by the uneven distribution of the emission sites at lower voltages.At lower applied voltages,some parts of the emitters can form emission sites while others cannot [33], and this is likely due to the variation among the emission strips.This variability can be attributed to differences in apex geometry(mainly radii),emitter-to-extractor distance,and material properties (such as pore size) along the emitter apex.These differences lead to different startup voltages for emission site formation at the different parts of the strips.Thus, at lower voltages, only a part of the emission sites can produce emission and the distribution of the emission sites is uneven.

    Besides, another interesting finding is that the plume perpendicular to the strips shows a higher divergence in both positive and negative modes compared with the plume parallel to the strips.Specifically, taking the angle containing 95%of the overall current as the plume angle,the half plume angle perpendicular to the strips is 41.2° at +3000 V and 40.1°at ?3000 V,while the plume distribution parallel to the strips is more concentrated, 19.4° at +3000 V and 19.7° at?3000 V.To get an in-depth understanding of the plume divergence difference in two cross sections, numerical simulation is carried out via the software package COMSOL Multiphysics version 5.5.

    In the simulation,a 3D model is established,as shown in figure 16(a).The sizes, materials, and circumstances are set the same as the experiments to reflect the actual conditions of the experiments.The combination of the extractor and the emitter array is placed in a sphere, which is defined as the calculation zone.The extractor is grounded,while the potential of the emitter array is set to either a positive or negative acceleration voltage.Particles (monomers and dimers) are released at an initial velocity of zero along the seven strips of the upper surface on the emitter array apex.The particle release position is marked in green in figure 16(b).The sequence type of meshing is set as a physics-controlled mesh,and the element size is selected as an extra-fine level so that the mesh of the emitter apex is fine enough to reflect the real geometry, as depicted in figure 16(c).

    Figure 17 shows the simulated trajectories of monomers and dimers at bothVa=3000 V andVa=?3000 V.Results further confirm that the plume divergence perpendicular to the strips is greater than parallel to the strips.From the COMSOL output file, the positions of the particles can be obtained.Taking the angle containing 95%of the total emitted particles as the plume angle, the simulated half plume angle perpendicular to the strips is 40.1° at +3000 V (an error of 2.7% with the measured value) and 39.7° at ?3000 V (an error of 1.0% with the measured value), while this value parallel to the strips is 17.6° at +3000 V (an error of 9.3%with the measured value) and 18.1° at ?3000 V (an error of 8.1% with the measured value).The simulated results are in good agreement with the experimental values.

    Figure 18.Arrow line graph of the electric field between two electrodes (the arrow length is proportional to the magnitude of the electric field intensity).

    The accelerating particles are driven by the electric force provided by the electric field.Investigating the electric field distribution between two electrodes as shown in figure 18,the main electric field lines match well with particle trajectories perpendicular to the strips.Thus, it can be inferred that the spatial plume profile is influenced by the electric field distribution between the emitter and the extractor to a great extent.

    5.Conclusions

    This paper presents a miniaturized ILIS based on an array of porous metal strips fabricated using WEDM combined with electrochemical etching.A series of experiments have been performed to gain an insight into some significant characteristics,such as theI–Vcharacteristic,energy properties,and the spatial plume profile of the ILIS.Results show that the ILIS is capable of producing ionic beams of up to 0.8 mA.Besides,the fractions of particles fragmented in the field-free region, fragmented in the acceleration region as well as without fragmentation, have been obtained by analyzing the RPA curves.The fragmentation rate of the emitted particles accounts for 48.8%in positive mode and 59.8%in negative mode.Further,voltage losses and energy efficiencies of the ILIS have been calculated, and results reveal that with the increase in applied acceleration voltage, the voltage loss rises, and the energy efficiency decreases.As to the spatial plume distribution, it is found that the plume perpendicular to the strips has a higher divergence than that parallel to the strips.The numerical simulation reveals that the electric field distribution of the ILIS may well lead to such a spatial plume profile.

    Acknowledgments

    This work is supported by National Natural Science Foundation of China (No.52075334).

    99久久中文字幕三级久久日本| 99久久精品热视频| 99热精品在线国产| 婷婷亚洲欧美| 人人妻人人澡欧美一区二区| 久久亚洲精品不卡| 久久这里只有精品中国| 校园人妻丝袜中文字幕| 亚洲四区av| 国产精品蜜桃在线观看 | 一本精品99久久精品77| 欧美人与善性xxx| 插逼视频在线观看| 3wmmmm亚洲av在线观看| 久久久精品欧美日韩精品| 丰满乱子伦码专区| 热99在线观看视频| 性插视频无遮挡在线免费观看| 简卡轻食公司| 国产成人a区在线观看| 白带黄色成豆腐渣| 亚洲精品自拍成人| 亚洲欧洲国产日韩| 午夜亚洲福利在线播放| 久久午夜亚洲精品久久| 啦啦啦观看免费观看视频高清| 两性午夜刺激爽爽歪歪视频在线观看| 日韩人妻高清精品专区| 国产蜜桃级精品一区二区三区| 国产国拍精品亚洲av在线观看| 少妇的逼好多水| 18禁裸乳无遮挡免费网站照片| 青春草国产在线视频 | 三级国产精品欧美在线观看| 国产精品福利在线免费观看| 国产精品三级大全| 日本爱情动作片www.在线观看| 一进一出抽搐动态| 国产探花在线观看一区二区| 美女脱内裤让男人舔精品视频 | 麻豆久久精品国产亚洲av| 日韩欧美一区二区三区在线观看| 免费搜索国产男女视频| 亚洲精品456在线播放app| 色5月婷婷丁香| 在线天堂最新版资源| 亚洲精品自拍成人| 少妇裸体淫交视频免费看高清| 99久久精品国产国产毛片| 婷婷色av中文字幕| 午夜免费男女啪啪视频观看| 高清午夜精品一区二区三区 | 欧美另类亚洲清纯唯美| 精品少妇黑人巨大在线播放 | 身体一侧抽搐| 亚洲综合色惰| 最后的刺客免费高清国语| 亚洲成a人片在线一区二区| 国产在视频线在精品| 国产伦理片在线播放av一区 | 欧美日韩乱码在线| 午夜激情福利司机影院| 99精品在免费线老司机午夜| 亚洲在久久综合| 久久热精品热| 国产亚洲91精品色在线| av卡一久久| av女优亚洲男人天堂| 一级毛片aaaaaa免费看小| 日本成人三级电影网站| www日本黄色视频网| 国产黄色小视频在线观看| 久久韩国三级中文字幕| 成人av在线播放网站| 久久人人爽人人片av| 99久久成人亚洲精品观看| 三级男女做爰猛烈吃奶摸视频| 好男人在线观看高清免费视频| 亚洲欧美日韩高清在线视频| 在线天堂最新版资源| 日韩中字成人| 欧美成人一区二区免费高清观看| 中文字幕免费在线视频6| 亚洲精华国产精华液的使用体验 | 成人无遮挡网站| 美女高潮的动态| 久久久久久久亚洲中文字幕| 91午夜精品亚洲一区二区三区| 99久久中文字幕三级久久日本| 一区二区三区四区激情视频 | 免费观看a级毛片全部| 九色成人免费人妻av| 夜夜爽天天搞| 国内久久婷婷六月综合欲色啪| 91精品一卡2卡3卡4卡| 国产极品精品免费视频能看的| 国产一区二区三区av在线 | 国产伦精品一区二区三区四那| 亚洲欧美日韩东京热| 亚洲激情五月婷婷啪啪| 有码 亚洲区| 亚洲欧美日韩卡通动漫| 国产单亲对白刺激| 联通29元200g的流量卡| 午夜福利在线观看免费完整高清在 | 女同久久另类99精品国产91| www.色视频.com| 最近中文字幕高清免费大全6| 午夜精品国产一区二区电影 | 精品国内亚洲2022精品成人| 久久久久久九九精品二区国产| 爱豆传媒免费全集在线观看| 69人妻影院| 亚洲不卡免费看| 老熟妇乱子伦视频在线观看| 国内揄拍国产精品人妻在线| 天堂√8在线中文| 美女cb高潮喷水在线观看| 欧美成人精品欧美一级黄| 国产精品电影一区二区三区| 亚洲成人久久性| 99久久精品热视频| 欧美xxxx黑人xx丫x性爽| 国产精品伦人一区二区| 欧美日韩国产亚洲二区| av在线观看视频网站免费| 中文精品一卡2卡3卡4更新| 岛国在线免费视频观看| 亚洲av电影不卡..在线观看| 免费看光身美女| 久久99精品国语久久久| 麻豆国产97在线/欧美| 18禁裸乳无遮挡免费网站照片| 久久精品国产亚洲av天美| 国产精品乱码一区二三区的特点| 菩萨蛮人人尽说江南好唐韦庄 | 久久中文看片网| 亚洲av熟女| 丝袜美腿在线中文| 欧美成人免费av一区二区三区| 人人妻人人看人人澡| 麻豆久久精品国产亚洲av| 在线免费观看的www视频| 嫩草影院精品99| 精品无人区乱码1区二区| 成人美女网站在线观看视频| .国产精品久久| 亚洲精品成人久久久久久| АⅤ资源中文在线天堂| 久久久久网色| 国产精品电影一区二区三区| 国产极品天堂在线| 有码 亚洲区| 亚洲国产色片| 看片在线看免费视频| 九九在线视频观看精品| 亚洲av免费在线观看| 亚洲高清免费不卡视频| 亚洲天堂国产精品一区在线| 亚洲人成网站在线播放欧美日韩| 全区人妻精品视频| 2022亚洲国产成人精品| 久久综合国产亚洲精品| 麻豆精品久久久久久蜜桃| 少妇的逼水好多| 两个人的视频大全免费| 久久精品影院6| 久久九九热精品免费| 老师上课跳d突然被开到最大视频| 国产精品.久久久| 少妇人妻一区二区三区视频| 久久鲁丝午夜福利片| 91aial.com中文字幕在线观看| 中文字幕熟女人妻在线| 亚洲人与动物交配视频| 亚洲人成网站在线观看播放| 日本色播在线视频| 99久久无色码亚洲精品果冻| 国产色爽女视频免费观看| 中文亚洲av片在线观看爽| 免费av毛片视频| 毛片一级片免费看久久久久| 级片在线观看| 三级男女做爰猛烈吃奶摸视频| 热99在线观看视频| 日韩高清综合在线| 亚洲在线观看片| 一级黄色大片毛片| 久久中文看片网| 一本一本综合久久| 久久久欧美国产精品| 亚洲精品日韩在线中文字幕 | 国产在线男女| 国内精品一区二区在线观看| 国产不卡一卡二| 成人漫画全彩无遮挡| 非洲黑人性xxxx精品又粗又长| 久久久久久久久中文| 身体一侧抽搐| 国产成人a∨麻豆精品| 久久久久久久亚洲中文字幕| eeuss影院久久| 91在线精品国自产拍蜜月| 精品人妻一区二区三区麻豆| 插逼视频在线观看| 又粗又硬又长又爽又黄的视频 | 亚洲成av人片在线播放无| 2022亚洲国产成人精品| 日本av手机在线免费观看| 嘟嘟电影网在线观看| 内射极品少妇av片p| 日本色播在线视频| 一区二区三区四区激情视频 | 亚洲精品自拍成人| 国产精品精品国产色婷婷| 成年免费大片在线观看| 岛国毛片在线播放| 亚洲国产精品sss在线观看| 我要搜黄色片| 尤物成人国产欧美一区二区三区| 国产色婷婷99| 国产精品久久久久久亚洲av鲁大| 噜噜噜噜噜久久久久久91| 久久人人爽人人片av| 老司机福利观看| 国产黄片美女视频| 久久婷婷人人爽人人干人人爱| 高清毛片免费观看视频网站| 欧美潮喷喷水| 99国产极品粉嫩在线观看| 免费看美女性在线毛片视频| 亚洲国产精品国产精品| 悠悠久久av| 国产成人福利小说| 99热这里只有是精品在线观看| 国产伦在线观看视频一区| 少妇人妻一区二区三区视频| 免费搜索国产男女视频| 亚洲在线自拍视频| 亚洲精品亚洲一区二区| 国产黄色小视频在线观看| 亚洲国产精品久久男人天堂| 天堂中文最新版在线下载 | 免费不卡的大黄色大毛片视频在线观看 | 村上凉子中文字幕在线| 日日摸夜夜添夜夜爱| 一本精品99久久精品77| 国产av在哪里看| 亚洲av熟女| 一级黄片播放器| 久久久久久久久久久免费av| 亚洲欧洲日产国产| 久久精品国产99精品国产亚洲性色| 国产探花极品一区二区| 在线免费观看的www视频| 国产美女午夜福利| 少妇熟女aⅴ在线视频| avwww免费| 春色校园在线视频观看| 亚洲欧美成人精品一区二区| 国产精品三级大全| 国产精品一二三区在线看| 国产精品久久久久久久电影| 亚洲最大成人av| 桃色一区二区三区在线观看| 午夜福利高清视频| 亚洲aⅴ乱码一区二区在线播放| 日韩欧美国产在线观看| 久久99精品国语久久久| 一个人观看的视频www高清免费观看| 中文字幕人妻熟人妻熟丝袜美| 人妻夜夜爽99麻豆av| 一级二级三级毛片免费看| 人妻久久中文字幕网| 你懂的网址亚洲精品在线观看 | 综合色丁香网| 亚洲三级黄色毛片| 精品久久久久久久久亚洲| 在线观看午夜福利视频| 中文字幕av在线有码专区| 精品人妻视频免费看| 亚洲五月天丁香| 亚洲欧美日韩高清在线视频| 久久国产乱子免费精品| 噜噜噜噜噜久久久久久91| 人妻久久中文字幕网| 天堂中文最新版在线下载 | 日韩av在线大香蕉| 亚洲三级黄色毛片| 精品国产三级普通话版| 秋霞在线观看毛片| 九九在线视频观看精品| 免费观看人在逋| 尤物成人国产欧美一区二区三区| 亚洲av中文av极速乱| 男女视频在线观看网站免费| 黄色欧美视频在线观看| 久久精品夜夜夜夜夜久久蜜豆| 毛片女人毛片| 两性午夜刺激爽爽歪歪视频在线观看| 日韩欧美三级三区| 日本欧美国产在线视频| 欧美+日韩+精品| 久久久久久伊人网av| 99热6这里只有精品| 色综合色国产| 亚洲欧洲日产国产| 欧美高清性xxxxhd video| 身体一侧抽搐| 亚洲欧美日韩卡通动漫| 亚洲丝袜综合中文字幕| www日本黄色视频网| 日韩亚洲欧美综合| 看免费成人av毛片| 亚洲无线在线观看| 真实男女啪啪啪动态图| 久久国产乱子免费精品| 精品人妻熟女av久视频| 99视频精品全部免费 在线| 色综合站精品国产| 在线播放国产精品三级| 亚洲av中文av极速乱| 精品久久久久久久末码| 久久久欧美国产精品| 热99在线观看视频| 午夜激情欧美在线| 人体艺术视频欧美日本| 综合色丁香网| 99久久精品一区二区三区| 精品久久久久久久久亚洲| av女优亚洲男人天堂| 两性午夜刺激爽爽歪歪视频在线观看| 校园人妻丝袜中文字幕| 成人综合一区亚洲| 亚洲精品成人久久久久久| 给我免费播放毛片高清在线观看| 伦精品一区二区三区| 黑人高潮一二区| 蜜桃亚洲精品一区二区三区| 又粗又爽又猛毛片免费看| 我要搜黄色片| 天堂av国产一区二区熟女人妻| 男人舔女人下体高潮全视频| 午夜精品在线福利| 国产白丝娇喘喷水9色精品| 国产成人福利小说| 欧美不卡视频在线免费观看| 一个人观看的视频www高清免费观看| 久久鲁丝午夜福利片| 久久精品综合一区二区三区| 看片在线看免费视频| 国产一区二区在线av高清观看| 亚洲欧美日韩高清在线视频| 亚洲国产欧美在线一区| 午夜精品一区二区三区免费看| 日韩欧美三级三区| 亚洲三级黄色毛片| 色综合站精品国产| 日日摸夜夜添夜夜添av毛片| 一级毛片久久久久久久久女| 日日摸夜夜添夜夜添av毛片| 国产一区二区在线av高清观看| 色视频www国产| 国产精品综合久久久久久久免费| 久久午夜福利片| 国产伦在线观看视频一区| 国产男人的电影天堂91| 欧美成人免费av一区二区三区| 中文字幕av成人在线电影| 男人的好看免费观看在线视频| 亚洲,欧美,日韩| 亚洲综合色惰| 如何舔出高潮| 日产精品乱码卡一卡2卡三| 亚洲精品国产av成人精品| 亚洲成av人片在线播放无| av又黄又爽大尺度在线免费看 | 日本免费a在线| 亚洲在线自拍视频| av免费观看日本| 国内精品宾馆在线| 免费看光身美女| 中文字幕熟女人妻在线| 国产欧美日韩精品一区二区| 深夜精品福利| 别揉我奶头 嗯啊视频| 国产午夜精品久久久久久一区二区三区| 亚洲四区av| 国产黄片美女视频| 国产一区二区在线av高清观看| 秋霞在线观看毛片| 亚洲成人久久性| 日韩欧美精品免费久久| 久久鲁丝午夜福利片| 综合色av麻豆| 日韩大尺度精品在线看网址| 久久久久久伊人网av| 晚上一个人看的免费电影| 国产免费一级a男人的天堂| 天天躁日日操中文字幕| 国内精品一区二区在线观看| 免费看美女性在线毛片视频| 中文字幕久久专区| 午夜a级毛片| 热99re8久久精品国产| 一个人看的www免费观看视频| 听说在线观看完整版免费高清| 天堂网av新在线| 欧美bdsm另类| 国产视频内射| 天天躁夜夜躁狠狠久久av| 亚洲va在线va天堂va国产| 看免费成人av毛片| 精品免费久久久久久久清纯| 精品久久久久久久人妻蜜臀av| 欧美日韩在线观看h| 人妻少妇偷人精品九色| 夜夜看夜夜爽夜夜摸| 国产成人a∨麻豆精品| 成人午夜精彩视频在线观看| 久久综合国产亚洲精品| 欧美日本视频| 波多野结衣高清作品| 亚洲第一区二区三区不卡| 99国产精品一区二区蜜桃av| 日韩在线高清观看一区二区三区| 国产精品一区二区三区四区久久| 欧洲精品卡2卡3卡4卡5卡区| 亚州av有码| 国产成人精品一,二区 | 直男gayav资源| av免费观看日本| 夜夜夜夜夜久久久久| 久久久久久九九精品二区国产| 日韩欧美精品免费久久| 麻豆成人av视频| 一级毛片久久久久久久久女| 18+在线观看网站| 麻豆国产97在线/欧美| 成人美女网站在线观看视频| 亚州av有码| 蜜桃久久精品国产亚洲av| 国产精品久久久久久亚洲av鲁大| 国产黄色视频一区二区在线观看 | 欧美区成人在线视频| 日日摸夜夜添夜夜爱| 一个人看视频在线观看www免费| 最近视频中文字幕2019在线8| 99精品在免费线老司机午夜| 国产一区亚洲一区在线观看| 中文在线观看免费www的网站| 欧美日韩在线观看h| 国产精品一区www在线观看| 日本黄色视频三级网站网址| 美女国产视频在线观看| 欧美成人免费av一区二区三区| 日韩精品青青久久久久久| 国产精品久久久久久精品电影小说 | 亚洲欧美精品自产自拍| 毛片女人毛片| 成人特级黄色片久久久久久久| 国产精品久久久久久av不卡| 啦啦啦观看免费观看视频高清| 亚洲av免费在线观看| 99久久成人亚洲精品观看| 草草在线视频免费看| 久久99热6这里只有精品| 中文字幕av成人在线电影| 日韩亚洲欧美综合| 亚洲欧美精品自产自拍| 久久久色成人| 国产午夜福利久久久久久| av卡一久久| 亚洲国产精品合色在线| 国产精品女同一区二区软件| 亚洲av中文字字幕乱码综合| 天堂√8在线中文| 成人综合一区亚洲| 亚洲国产欧洲综合997久久,| 中文在线观看免费www的网站| 国产精品久久久久久精品电影| 性欧美人与动物交配| av在线观看视频网站免费| 三级毛片av免费| 日本-黄色视频高清免费观看| 亚洲高清免费不卡视频| 免费无遮挡裸体视频| 日本av手机在线免费观看| 国产av一区在线观看免费| 中文欧美无线码| 免费看日本二区| 免费观看a级毛片全部| 久久精品久久久久久久性| 亚洲激情五月婷婷啪啪| 国产精品电影一区二区三区| 久久国产乱子免费精品| 热99re8久久精品国产| 国产黄片美女视频| 久99久视频精品免费| 亚洲精品成人久久久久久| av在线老鸭窝| 精品少妇黑人巨大在线播放 | 欧美丝袜亚洲另类| 看非洲黑人一级黄片| 亚洲欧美日韩高清在线视频| 国产精品女同一区二区软件| 人妻制服诱惑在线中文字幕| 国产老妇伦熟女老妇高清| 亚洲自拍偷在线| 久久人妻av系列| 我要看日韩黄色一级片| 人体艺术视频欧美日本| a级毛色黄片| 亚洲av.av天堂| 欧美变态另类bdsm刘玥| 亚洲婷婷狠狠爱综合网| 久久欧美精品欧美久久欧美| 欧洲精品卡2卡3卡4卡5卡区| 变态另类丝袜制服| 国产精品.久久久| 国产精品一区www在线观看| 26uuu在线亚洲综合色| 免费观看精品视频网站| 99国产精品一区二区蜜桃av| 成人毛片a级毛片在线播放| 91精品国产九色| 国产男人的电影天堂91| 不卡视频在线观看欧美| 麻豆国产av国片精品| 精品人妻视频免费看| 国产日韩欧美在线精品| 日日干狠狠操夜夜爽| 永久网站在线| 国产视频内射| 国产精品免费一区二区三区在线| 久久亚洲国产成人精品v| 高清毛片免费观看视频网站| 黄片无遮挡物在线观看| 精品久久久噜噜| 日本免费a在线| 日本黄大片高清| 女的被弄到高潮叫床怎么办| 国产成人精品婷婷| 久久久久国产网址| 简卡轻食公司| 九九久久精品国产亚洲av麻豆| 久久久欧美国产精品| 午夜老司机福利剧场| 国产精品福利在线免费观看| 亚洲在线自拍视频| 亚洲av熟女| 国产乱人视频| 亚洲第一电影网av| 国产一区二区三区在线臀色熟女| 网址你懂的国产日韩在线| 熟女人妻精品中文字幕| 久久久成人免费电影| 国产精品人妻久久久久久| 午夜福利在线观看吧| 亚洲av中文av极速乱| 国产亚洲av片在线观看秒播厂 | 村上凉子中文字幕在线| 国产成人a区在线观看| 日日干狠狠操夜夜爽| 22中文网久久字幕| 日本在线视频免费播放| 18禁在线播放成人免费| 人妻夜夜爽99麻豆av| 国产精品久久电影中文字幕| av又黄又爽大尺度在线免费看 | 色哟哟哟哟哟哟| 如何舔出高潮| 欧美一区二区国产精品久久精品| 国产精品1区2区在线观看.| 欧美色欧美亚洲另类二区| а√天堂www在线а√下载| 国产老妇女一区| 久久精品国产99精品国产亚洲性色| 99在线人妻在线中文字幕| 国产精品一区二区在线观看99 | 欧美性猛交╳xxx乱大交人| 国产午夜福利久久久久久| 午夜精品国产一区二区电影 | 美女被艹到高潮喷水动态| 国产精品久久视频播放| 国产精品无大码| 一级毛片aaaaaa免费看小| videossex国产| 国产 一区精品| 欧美最黄视频在线播放免费| 免费搜索国产男女视频| av在线观看视频网站免费| 成人综合一区亚洲| 精品久久久久久久末码| 欧美日韩在线观看h| 日韩一区二区视频免费看| 国产精品久久视频播放| 国产真实伦视频高清在线观看| 日韩大尺度精品在线看网址| 嫩草影院入口| 国产真实伦视频高清在线观看| 国产精品野战在线观看| 中国美女看黄片| 国产一区二区三区av在线 | 91久久精品国产一区二区成人| 成人毛片a级毛片在线播放| 中文字幕人妻熟人妻熟丝袜美| 国产高清有码在线观看视频| 在线观看av片永久免费下载| 国产色爽女视频免费观看| 欧美在线一区亚洲| 日日啪夜夜撸| 精品国内亚洲2022精品成人| 日韩高清综合在线| 国产不卡一卡二|