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

    Simulation of liquid cone formation on the tip apex of indium field emission electric propulsion thrusters

    2024-04-24 06:17:24YimingSUN孫逸鳴HanwenDENG鄧涵文XinyuLIU劉欣宇andXiaomingKANG康小明
    Plasma Science and Technology 2024年4期

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

    1 School of Mechanical Engineering,Shanghai Jiao Tong University,Shanghai 200240,People’s Republic of China

    2 Innovation Academy for Microsatellites,Chinese Academy of Sciences,Shanghai 201210,People’s Republic of China

    Abstract Field emission electric propulsion (FEEP) thrusters possess excellent characteristics,such as high specific impulse,low power requirements,compact size and precise pointing capabilities,making them ideal propulsion devices for micro-nano satellites.However,the detection of certain aspects,such as the evolution process of the liquid cone and the physical quantities at the cone apex,proves challenging due to the minute size of the needle tip and the vacuum environment in which they operate.Consequently,this paper introduces a computational fluid dynamics (CFD) model to gain insight into the formation process of the liquid cone on the tip apex of indium FEEP.The CFD model is based on electrohydrodynamic (EHD) equations and the volume of fluid (VOF) method.The entire cone formation process can be divided into three stages,and the time-dependent characteristics of the physical quantities at the cone apex are investigated.The influences of film thickness,apex radius size and applied voltage are compared.The results indicate a gradual increase in the values of electrostatic stress and surface tension stress at the cone apex over an initial period,followed by a rapid escalation within a short duration.Apex configurations featuring a small radius,thick film and high voltage exhibit a propensity for liquid cone formation,and the cone growth time decreases as the film thickness increases.Moreover,some unstable behavior is observed during the cone formation process.

    Keywords: FEEP,needle emitter,liquid cone formation,CFD simulation

    1.Introduction

    As a novel electric propulsion technology,field emission electric propulsion (FEEP) offers precise thrust ranging from severalμN to hundreds ofμN while maintaining low thrust noise [1].FEEP possesses several advantages,including high specific impulse,compact form factor and long service life,making it highly suitable for drag compensation,attitude control and orbit maintenance for nano-satellites [2].According to the wetting properties of the propellant,FEEP can be categorized into needle,capillary,and porous types [3–6].The solid needle-type FEEP features a micro-sized tungsten tip,which exhibits high electrical impedance,operates at a low firing voltage,and facilitates thrust clustering to generate substantial propulsion [7].Indium serves as a recommended propellant due to its melting point (156 °C),enabling it to solidify and withstand vibrations during launch.

    The development of indium-needle FEEP is based on the technology of the liquid metal ion source (LMIS) [8].LMIS comprises a heater assembly,an extractor,and a needle emitter that is wetted by liquid metal,as illustrated in figure 1.When the heater assembly heats liquid metal propellant to its melting point and a high electric potential is applied within the space between the emitter and extractor,the liquid bumps on the tip apex due to the combined effects of electrostatic force and surface tension stress.Subsequently,it undergoes growth,forming a conical cone,a so-called Taylor cone,with a cusp protrusion at the apex,until the electric field at the cusp apex reaches a critical level,with ion emission occurring [9].Simultaneously,the presence of fluid negative pressure at the cone apex propels the liquid metal to flow along microgrooves on the surface towards the apex of the cone,which serves to replenish the already ejected ions [10,11].

    The formation and evolution of the liquid cone on the tip apex constitute the fundamental physical processes underlying FEEP.These processes involve complex interactions among electrostatic forces,surface tension and fluid flow.Taylor was the first to establish a theoretical analysis of liquid cone formation,considering the static equilibrium between electrostatic forces and surface tension pressure and ignoring liquid flow and ion emission [12].Under these assumptions,the half-cone angle was derived as 49.3° when a critical voltage was applied between the emitter and extractor.Kingham and Swanson adopted a model of a jet-like protrusion on the end of a Taylor cone and found that the cusp apex radii were typically in the order of a few tens of angstroms with protrusion lengths ranging from 400 ? to 500 ? [13].Early models of cone-jet formation focused on the quasi-equilibrium of the liquid cone surface.In later studies,Burton and Taborek simulated the dynamic evolution process of charged inviscid droplets and observed that the droplets become unstable with small perturbations,ultimately forming a lemon shape with a sharp tip [14].Suvorov presented an electrohydrodynamic model that described the dynamic behavior of the Taylor cone using the Navier–Stokes equations and considered free liquid boundaries subject to Maxwell field stress,surface tension stress and viscous effects [15].However,the calculations in this model assumed a flat needle tip,whereas the fabricated conical needle tip typically has an approximately hemispherical apex.Several scholars have indirectly or directly observed the formation of the liquid cone.Wagneret alobserved the rapidly frozen liquid cone on the tip apex while operating at approximately 20μA [16].In addition,dynamic instabilities were also found on the needle shank behind the Taylor cone by a 200 keV transmission electron microscope,resulting in microdroplet emission.Praprotniket alobserved the in situ shape of the liquid cone in indium LMIS with ion emission using a 1 MeV electron transmission microscope [17].It was noted that the jet length increased and the cone angle decreased with increasing emission current.However,due to the interference of electric fields and the limitation of exposure speed,precise observation of the transient response during operation using an electron microscope is challenging.Recently,some CFD codes have been employed to simulate the dynamic behavior of cone formation of aqueous solution or ionic liquid.Senet alsimulated the cone-jet electrospray on carbon nanofiber emitters using a CFD code(FLOW-3D) and demonstrated the feasibility of these simulations [18].Liuet alsimulated the meniscus formation of an ionic liquid in the electrospray process and found that the alternation of polarity can delay the cone growth time [19].Heet alsimulated the droplet emission process on capillary standing waves under the dual action of electric field and vibration by a CFD software FLUENT [20,21].In addition to CFD methods,molecular dynamics has been utilized to investigate the formation of the liquid cone by computing the interaction between atom and field.Zhanget alstudied Taylor cone formation and the electrospray jet process of capillary emitters at an atomistic scale [22].However,due to computational costs,the model was considerably simplified relative to the real system using scaling laws [23].

    Despite significant effort dedicated to investigating the liquid cone phenomenon,little attention has been given to the cone formation process and the evolution of physical quantities at the cone apex over time.This lack of attention can be attributed to the challenges associated with observation and measurement.Since the needle emitter cannot be operated and directly observed in the atmosphere due to the liquid metal’s susceptibility to oxidation,this creates the issues of focus and capture in the vacuum chamber.The apex radius of the needle is typically several micrometers,and the distance of the counter electrodes is in the order of tens of microns.The jet diameter of the cusp on the cone apex in the range of sub-microns is too small to be observed through optical cameras.In addition,measuring physical quantities,such as the electric field on the fluid surface and pressure beneath the cone apex,is nearly impossible.Therefore,numerical simulation provides a more feasible approach for studying the formation of the liquid cone in FEEP.In this regard,we have developed a computational fluid dynamics (CFD) model of FEEP utilizing a solid needle emitter.This model aims to enhance the understanding of the fundamental characteristics and formation process of the liquid cone at the tip apex.

    The structure of the paper is organized as follows.Section 2 provides an overview of the electrohydrodynamic theory,which is employed to describe the behavior of the liquid during the cone formation process.Section 3 presents the analysis and establishment of the CFD model.In section 4,the evolution of the liquid cone on the tip apex is investigated,and the time-dependent characteristics of physical quantities at the cone apex are revealed.In addition,the influence of film thickness,apex radius and applied voltage on cone growth time and free surface shape is examined and compared.This section also addresses the discussion of instabilities observed during the formation process of the Taylor cone.Finally,in section 5,a summary of the work is provided.

    2.Theory

    The CFD simulation is based on the Taylor–Melcher dielectric fluid model and the EHD equations [24,25].To precisely capture the free surface between the air and liquid phases,the volume of fluid (VOF) method [26] is applied.This is facilitated through the utilization of FLOW-3D,a dedicated CFD software.FLOW-3D offers various physics modules that are internally integrated,allowing users to select and couple multiple modules for calculations.In this simulation,the modules of electro-mechanics,surface tension and viscosity are specifically employed to analyze the liquid cone formation process on the tip apex of the FEEP thruster.

    The equation of fluid continuity and the Navier–Stokes equation were used to describe the dynamic behavior of the liquid cone formation.Assuming that the liquid is incompressible (of constant densityρ) and has the constant viscosityμ.The gravity term can be omitted,considering the size of the liquid cone (several micrometers).It can be described as,

    whereuis the fluid velocity,Pis the hydrodynamic pressure andFeis the body force.The convective effects induced by temperature primarily stem from viscosity and conductivity.Given the constant operational temperature of the thruster,the size of the liquid cone (several micrometers),and the transient formation process (several microseconds),it is reasonable to neglect variations in viscosity and conductivity,treating them as constants.Therefore,electromechanical forces are the main driving force inFe,and the electrostatic force is expressed as,

    whereρeis the net free charge density,Eis the electric field strength,ε0is the vacuum permittivity andεlis the fluid dielectric constant.The first term on the right is Coulomb force due to the free charge on the liquid surface,and the second term is polarization force arising from polarization effects.

    The governing equation for the electric potential and electric field in the fluid is given by the Laplace equation:

    The law of charge conservation on the vacuum-liquid interface is as follows:

    wherenis the unit normal vector andσis the electric conductivity of the liquid.At the interface,there is a jump in the electric field in the normal direction due to surface charge density,while the electric field in the tangential direction is continuous:

    whereφlis the electric potential in the liquid,φvis the electric potential in the vacuum,εvis the relative permittivity of the vacuum,qis the electric charge quantity andtis the unit tangential vector.The mechanical equilibrium condition at the normal direction of the interface is,

    whereTeis the electrostatic stress tensor,Tfis the fluid stress tensor andγis the surface tension coefficient.In the VOF method,the motion of the interface is governed by the kinematic equation:

    The interface is tracked by the volume fraction functionF:

    where 1 represents the cell occupied by the liquid,and 0 is the vacuum.Equations are solved in FLOW-3D by choosing appropriate boundaries and initial conditions.

    3.Modelling and validation

    3.1. Modelling of the liquid cone formation process

    In the FEEP system,the needle emitter is typically fabricated by electrochemical etching.The apex radius of the tip ranges from several micrometers,as shown in figure 2(b),while the liquid film thickness is in the range of tens to hundreds of nanometers [10].However,the dimensions of the needle tip length and extractor diameter are in the order of several millimeters,as depicted in figure 2(a).This significant difference in dimensions makes it challenging to simulate the entire flow on the needle tip surface and apex.

    In this study,our focus is on the process by which the liquid film forms a liquid cone on the needle tip.Consequently,the needle tip geometry can be approximated by a solid hemisphere uniformly covered by a liquid film,as illustrated in figure 3(a).The physical parameters of liquid indium [27–29] including fluid densityρ,viscosityμ,electrical conductivityσand surface tension coefficientγare listed in table 1.Figure 3(b) depicts the solid disk extractor placed 100μm above the emitter.A positive voltageΦ=V0is applied to the tip hemisphere,while the extractor is set to zero voltageΦ=0,creating a high electric field between the two electrodes.The no-slip wall condition of fluidu=0 is imposed on the tip surface.As the tip is approached,the area near the apex has no grooves [30],as shown in figure 2(b).

    Table 1.Physical properties of the liquid indium.

    Thus,the tip apex can be considered smooth in this model.The geometry in this model exhibits cylindrical symmetry around thez-axis,allowing the use of cylindrical coordinates and mesh for computational simplification.Mesh convergence study was conducted in figure 4,in order to ensure that the results are mesh independent.The results indicate that a mesh size less than or equal to 0.25μm has sufficient refinement to obtain a better resolution.Considering the balance between resolution and computational efficiency,the fluid domain employs a mesh size of 0.25μm,while the spatial domain is assigned a mesh size of 0.5μm.The simulation domain encompasses a radial dimension of 6μm and an axial dimension of 120μm.Boundaries 1 and 2 are set as the outflow boundary,and boundary 3 represents the boundary wall.

    3.2. Validation of the established model

    To verify the accuracy of the established model,an experimental approach was conducted by freezing the liquid cone on the tip apex and observing it using an electron microscope.In the experiment,a tungsten tip with an apex radius of 5μm was wetted with liquid indium,resulting in a film thickness of 0.4μm at the tip apex estimated from the SEM image.The distance between the counter electrode and the emitter was adjusted to 100μm.The liquid cone was frozen while operating at emission current of 18μA,a voltage of 7.5 kV and a pressure of 5×10-4Pa.The simulation used the same parameters as those in the experimental setup.

    Figure 5(a) displays the scanning electron microscope(SEM) image of the frozen Taylor cone on the tip apex obtained from the experiment.Figure 5(b) illustrates the simulation results of the liquid cone shape.The cone shape obtained from the simulation exhibits a better consistency than the experimental results.The measured cone angle of the frozen cone is 98.5°,with an error of 6.7% when compared to the simulated value of 91.9°.Furthermore,the model undergoes validation by simulating the onset voltage,and the results are compared with the experimental data presented in figure 6.To ensure consistency,identical wetting processes and parameters are employed,maintaining a uniform film thickness of 0.4μm on the needle tip,with a fixed distance of 200μm between the counter electrode and the emitter.The results show a slight discrepancy between simulation and experimentation,with the simulation values averaging approximately 4% higher than the experimental results.These results indicate a reasonable agreement between the simulation and experimental observations,validating the accuracy of the established model.

    3.3. Parameters set in the simulation

    The tip apex radius can be fabricated at a desired value through electrochemical etching and a roughening process.However,it is common for the apex radius to be processed below 10μm due to considerations of low onset voltage and high mass efficiency [31,32].For the purpose of comparison,a typical tip radius of 5μm was chosen,while 3μm and 7μm were used as alternative values.The film thickness represents the degree of wetting,with a thicker film indicating a better wetting degree.A minimum film thickness of 0.3μm is required [33].Consequently,a film thickness ranging from 0.3 to 0.7μm was employed to examine the impact of different film thicknesses.In the simulation,the applied voltage is set near the onset voltage corresponding to the apex radius.Furthermore,different voltages are applied to the same apex radius in order to investigate the influence of voltage.

    4.Results and discussion

    4.1. Evolution of the liquid cone

    The parameters utilized in the current model are as follows:an apex radius of 5μm,a film thickness of 0.5μm,an applied voltage of 6210 V.The counter electrode is positioned at a distance of 100μm from the tip apex.The formation process of the liquid cone can be divided into three stages,as depicted in figure 7.Stage 1: the liquid film is a hemisphere shape at the initial stage.Over time (0.6–1.2μs),the liquid gradually converges from both sides of the tip towards the apex under the influence of the electric field.At 2.09μs,a thicker fluid accumulates at the apex,resulting in a bump shape.It is observed that the formed conical liquid cone has an approximate half angle of 42° at 2.4μs,slightly smaller than the theoretical value of the Taylor cone (49.3°).This deviation is attributed to the fact that the system is not in static equilibrium during the movement of the liquid.In this first stage,the initial hemisphere liquid film redistributes due to the electrically driven process.Stage 2: from 2.4 to 2.71μs,the height of the liquid continues to increase,and the cone apex becomes sharper.By 2.71μs,the Taylor cone is essentially formed.During this stage,the apex radius of the cone gradually decreases,approaching the characteristics of a Taylor cone.Stage 3: from 2.71 to 2.74μs,the cone cusp is further protruded due to the high electric field acting on the small surface of the cusp apex.At 2.74μs,the cone jet is formed with a half angle of 38.7°.Subsequently,microdroplets are emitted in the cone jet area.In this stage,the cone jet is fully developed,and microdroplet emission occurs.Figure 8 illustrates the electric field intensity along thez-axis (symmetry axis) during the evolution of the liquid cone,and the center of the hemisphere is thez-zero point.Initially,the electric field intensity is low,but as the liquid cone gradually forms,the electric field at the cone apex increases,reaching values ranging from 108to 109V·m-1.The electric field at the cone apex reaches 1.3×109V·m-1at 2.7μs,achieving the ionized field threshold.It is important to note that FLOW-3D is a simulation software primarily designed for fluid dynamics.As such,it can effectively simulate the entire formation process of the liquid cone.However,it should be highlighted that the subsequent ion emission mechanism is field emission.Consequently,FLOW-3D can only accurately simulate the formation process of the cone jet,and after that it will transit into the droplet emission stage instead of the actual ion emission state.

    4.2. Time-dependent characteristic of physical quantities at the cone apex

    Figure 9(a) illustrates the time dependence of electrostatic stress and surface tension stress at the cone apex.These physical quantities exhibit a distinct characteristic of gradually increasing values over an extended period,followed by a rapid escalation within a short duration.The evolution of stress growth can also be divided into three stages corresponding to liquid cone growth.In the initial stage,stress increases slowly over time,with the shape of the liquid cone also changing slowly during this time.Advancing to the second stage,due to the liquid cone shape basic formation,stress increases at a quicker rate.Finally,in the third stage,because of the cone cusp protrusion and high electric field action,stress experiences a rapid and run-away escalation.It can be observed at 2.74μs that the electrostatic stress surpasses and is much larger than the surface tension stress.The Bernoulli equation at the cone apex can be expressed as,

    wherePis the pressure,ρis the density,uis the fluid velocity,Teis the electrostatic stress,Tfis the surface tension stress and δ is the strain due to the fluid stretching.From the equation and figure 9(a),it can be analyzed that the unequal values of electrostatic and surface tension stress include contributions from kinetic energy density and negative pressure beneath the cusp.The negative pressure plays a role in driving the flow of liquid on the surface of the needle.Figure 10 presents the pressure distribution within the liquid cone during its evolution.Notably,there is a region of negative pressure beneath the cone apex.Figure 9(b) is the shape of the liquid level obtained by approximating the fluid mesh points.The horizontal axis is in the radial direction and the vertical axis is in thezdirection (symmetry axis).By comparing the surface shapes throughout the evolution process in figure 9(b),it becomes apparent that the onset of the run-away behavior coincides with the formation of a quasi-conical surface shape.As the conical liquid cone develops,the cone apex area and radius shrink with increasing cone height,leading to the rapid growth of electrostatic stress and surface tension stress.

    Figure 1.Schematic of the needle-type FEEP setup.

    Figure 2.Comparing the macro size of the FEEP thruster with the micro size of the tip apex.(a) Schematic setup of the FEEP thruster and (b) SEM image of the needle tip and tip apex.

    Figure 3.Modeling and meshing of the emission process.(a) Model of the needle tip wetted by a film of liquid metal and (b) meshes and boundary conditions.

    Figure 4.Maximum electric field strength versus mesh size.

    Figure 5.Comparison of the cone shape in a frozen Taylor cone with the simulation result.(a) Frozen Taylor cone of liquid metal indium at the tungsten tip and (b) the simulation result.

    Figure 6.Onset voltage versus tip apex radius.

    Figure 7.The formation process of the Taylor cone at the apex radius of 5 μm,applied voltage of 6210 V and film thickness of 0.5 μm.(a) First stage,(b) second stage and (c) third stage.

    Figure 8.Evolution of the electric field strength along the z-axis in the liquid cone formation.

    Figure 9.Time-dependent behavior at the cone apex.(a) Electrostatic stress and the surface tension stress at the cone apex and (b) variation of the free surface shape.

    Figure 10.The pressure distribution inside the liquid cone during evolution at 2.56 μs.

    In figure 11(a),the cone-apex velocity along thez-axis exhibits behavior similar to the time-dependent stress acting on the liquid surface.During the stage of liquid accumulation,the velocity increases gradually over an extended period and then rapidly escalates.Prior to the occurrence of emission,the velocity of the cone apex reaches a value of 33 m·s-1.In addition,the position of the cone apex also follows the same pattern of motion.Although the ion emission process cannot be represented with FLOW-3D,the charged microdroplets generated during the actual highcurrent situation can be reflected and investigated in the simulation.Figure 11(b) illustrates the velocity distribution of the charged microdroplets along thez-axis.It is observed that the velocity of the microdroplets jetted from the cusp exhibits a linear distribution and the maximum velocity is 35100 m·s-1.

    Figure 11.Motion of liquid surface and microdroplet emission at the cone apex.(a) Variation of the fluid position velocity of the cone apex and (b) velocity distribution of the jetted microdroplet.

    4.3. Influence of the film thickness,apex radius and applied voltage

    Figure 12 presents the influence of film thickness on the cone growth time,considering an apex radius of 5μm and an applied voltage of 6210 V.The results indicate that the cone growth time decreases as the film thickness increases.In other words,a thinner film requires more time to form the Taylor cone.This behavior can be attributed to the difficulty of liquid flow towards the apex when a thin film is present,due to the viscous drag near the solid boundary.At thinner film thicknesses,the viscous effect becomes dominant,resulting in a longer time for the liquid cone to form.As the film thickness increases,the influence of viscosity on the formation time gradually weakens.

    Figure 12.Formation time of the Taylor cone as a function of film thickness.

    Figure 13.Formations of the liquid cone with different apex radii of 3 μm,5 μm and 7 μm.

    Figure 14.Formation of the Taylor cone at different applied voltages.

    Figure 15.Unstable formations of the liquid cone.(a) Microdroplets occurring at the bottom of the Taylor cone with the film thickness of 0.8 μm and voltage of 7200 V and (b) fragmentation of the liquid cone in the formation process with the film thickness of 0.5 μm and voltage of 7980 V.

    In figure 13,the formation of the liquid cone at different apex radii is depicted,considering a film thickness of 0.5μm and an applied voltage of 5820 V.It can be seen that the tip apex with a radius of 5μm exactly forms the Taylor cone,as the applied voltage approaches the onset voltage.Furthermore,the tip apex with a smaller radius of 3μm is more prone to forming the Taylor cone.However,the tip apex with a larger radius fails to generate a liquid cone.This discrepancy can be explained by considering the interplay between the liquid film and the walls.In the case of a thick film (0.5μm) relative to the apex radius of 3μm,the influence of viscosity can be neglected.Moreover,the sharp apex at the initial stage results in a higher electric field at the liquid surface and promotes the redistribution process driven by the electrostatic force.Consequently,the combination of a blunt apex and a thin film leads to the unsuccessful formation of the liquid cone at an apex radius of 7μm.

    In figure 14,the influence of different voltages on the cone formation time is investigated,considering an apex radius of 5μm and a film thickness of 0.5μm.The results show that the cone growth time decreases with the applied voltage.Because of the higher voltage,the electrostatic stress at the liquid surface impels the faster formation of the liquid cone.In addition,at a low voltage of 5820 V,the surface of the liquid cone exhibits a convex shape,while the jet of the liquid cone becomes slender,and the height rises with the increase in voltage.This phenomenon can be explained by the self-shaping characteristic of LMIS,where changes in the field strength not only extend the cone-jet length but also result in a decrease in the half angle,and a concave shape of the cone surface.The simulation results presented in figure 14 are in good agreement with the experiments conducted by Benassayaget al[34],who observed the dynamic formation of the liquid cone on the tip of Ga-LMIS using anin situ3 MeV electron transmission microscope.Furthermore,Praprotniket al[17] also observed the dynamic behavior of the indium liquid cone and the cone jet process using anin situmethod.

    4.4. Instabilities during the formation of the Taylor cone

    It is found in the simulations that Taylor cones do not form stably under all conditions.Two typical cases analyzed below demonstrate the potential instabilities that can occur during the formation of the Taylor cone,highlighting the importance of parameters,such as film thickness,applied voltage and wetting properties,in achieving stable cone formation and emission processes.

    The first case of the instabilities of the Taylor cone in the simulation is investigated with a film thickness of 0.8μm and an applied voltage of 7200 V,as shown in figure 15(a).When the liquid film is thick and the applied voltage (emission current) is high,the liquid surface is unstable and turbulent flow occurs in the cone formation process.As a result,microdroplets,known as Faraday droplets,with diameters of a few tenths of a micron are emitted from the shank of the liquid cone [16,35].The instability is likely triggered by surface waves that form during the liquid cone formation process.It is found that this situation occurs in a specific combination of liquid film thickness and high-voltage parameters.

    In the second case,depicted in figure 15(b),a film thickness of 0.5μm and an applied voltage of 7980 V are used.When the liquid film is relatively thin and the voltage is high(or the emission current is high),the speed at which the bump of the liquid at the apex moves is greater than the speed at which the surrounding fluid flows towards the apex,resulting in cone jet fragmentation and the formation of microdroplets.This results in the liquid being unable to replenish the fluid needed for cone formation,leading to the collapse of the cone and the emission of Rayleigh droplets[36,37].The second instability case has smaller scale of emitted droplets and higher frequency compared with the first case.This situation can be observed in needles with a poor wetting film when there is oxide contamination on the liquid metal surface or a discontinuous crack liquid film due to thermo stress during cyclic heating of the thruster,causing high flow impedance and forming unstable emission.

    5.Conclusion

    The CFD simulations presented in this work are demonstrated to be capable of predicting the physical quantities in the liquid cone,and provide valuable reference in FEEP design and performance evaluation.The following conclusions can be drawn:

    (1) Employing EHD equations and the VOF method to describe fluid behavior and track the free surface,the model accurately captures the formation of the liquid cone on the micro tip apex.The formation process can be divided into three stages: liquid hump at tip apex,liquid cone basic formation,cone apex protrusion and cone jet formation.

    (2) The simulation results reveal interesting behavior of various physical quantities at the cone apex such as electrostatic stress,surface tension stress and velocity.These quantities exhibit a run-away behavior,slowly increasing over a long time and then rapidly becoming very large within a short time.This behavior is closely related to the formation of the quasi-cone structure.

    (3) The film thickness is identified as a significant factor influencing the formation of the liquid cone on the tip apex.The cone growth time decreases as the film thickness increases due to the viscous drag acting close to the solid boundary.However,when the liquid film becomes sufficiently thick,the influence of viscosity becomes negligible.

    (4) The simulation also indicates that the formation of the liquid cone is more likely with a smaller apex radius under the same film thickness and applied voltage.The higher the applied voltage,the shorter the cone growth time.Furthermore,increasing the applied voltage leads to an increase in the protrusion length of the cusp and a concave shape of the free surface,resulting in a decrease in the cone angle.

    (5) In simulation,microdroplets are sometimes emitted from the shank of the liquid cone.In addition,it is noted that the discontinuity of the liquid metal film caused by poor wetting or cyclic heating can lead to the failed formation and collapse of the Taylor cone.

    Overall,the CFD simulation presented in this work provides valuable insight into the behavior of the liquid cone formation and its dependence on various parameters,contributing to the understanding and optimization of FEEP systems.

    Acknowledgments

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

    av国产精品久久久久影院| 久久人人爽人人爽人人片va| 日日啪夜夜撸| 国产永久视频网站| 高清午夜精品一区二区三区| 永久免费av网站大全| 自拍偷自拍亚洲精品老妇| 亚洲av中文av极速乱| 国产乱人视频| 日韩av不卡免费在线播放| 亚洲美女视频黄频| 黄色怎么调成土黄色| 亚洲成人久久爱视频| 国产高清三级在线| 精品国产乱码久久久久久小说| 久久精品国产自在天天线| 久久精品综合一区二区三区| 老司机影院成人| 直男gayav资源| 99re6热这里在线精品视频| 三级国产精品欧美在线观看| 久久久国产一区二区| 能在线免费看毛片的网站| 婷婷色av中文字幕| 欧美老熟妇乱子伦牲交| 亚洲真实伦在线观看| 色哟哟·www| 成人无遮挡网站| 嘟嘟电影网在线观看| 日韩在线高清观看一区二区三区| 亚洲精品视频女| 国产男人的电影天堂91| 国精品久久久久久国模美| 欧美xxxx性猛交bbbb| 久久久久久久亚洲中文字幕| 午夜福利网站1000一区二区三区| 又粗又硬又长又爽又黄的视频| 欧美成人a在线观看| 麻豆精品久久久久久蜜桃| 80岁老熟妇乱子伦牲交| 久久精品人妻少妇| 日本猛色少妇xxxxx猛交久久| 最后的刺客免费高清国语| 18禁动态无遮挡网站| 欧美高清成人免费视频www| 国产免费一级a男人的天堂| 久久久久久久久大av| 亚洲欧美日韩卡通动漫| 成人亚洲精品av一区二区| 人人妻人人爽人人添夜夜欢视频 | 精品久久久精品久久久| 国产欧美日韩一区二区三区在线 | 久久97久久精品| 青春草国产在线视频| 国产伦精品一区二区三区四那| 大话2 男鬼变身卡| 国产男女内射视频| 免费观看的影片在线观看| 欧美成人午夜免费资源| 亚洲欧美日韩东京热| 亚洲精品久久久久久婷婷小说| 久久久国产一区二区| 久久国内精品自在自线图片| 极品少妇高潮喷水抽搐| 你懂的网址亚洲精品在线观看| 国产精品秋霞免费鲁丝片| 久久99热这里只频精品6学生| 免费黄色在线免费观看| 熟女av电影| 男人舔奶头视频| 亚洲av成人精品一二三区| 卡戴珊不雅视频在线播放| 日本黄色片子视频| 在线精品无人区一区二区三 | a级毛片免费高清观看在线播放| 99久国产av精品国产电影| 高清在线视频一区二区三区| 亚洲电影在线观看av| 少妇人妻久久综合中文| 少妇裸体淫交视频免费看高清| 日韩一区二区三区影片| 人人妻人人爽人人添夜夜欢视频 | 成人午夜精彩视频在线观看| 日产精品乱码卡一卡2卡三| 成人特级av手机在线观看| 麻豆成人午夜福利视频| 波多野结衣巨乳人妻| 三级经典国产精品| 最近最新中文字幕免费大全7| 久久午夜福利片| 欧美日韩亚洲高清精品| 精品国产一区二区三区久久久樱花 | 欧美成人精品欧美一级黄| 一级av片app| 国产一区亚洲一区在线观看| 亚洲av.av天堂| 国产亚洲最大av| 一区二区三区四区激情视频| 一本一本综合久久| 丝瓜视频免费看黄片| 高清在线视频一区二区三区| 日韩一区二区视频免费看| 亚洲精品一二三| 亚洲在久久综合| 在线天堂最新版资源| 亚洲精品日韩在线中文字幕| 亚洲成人中文字幕在线播放| 各种免费的搞黄视频| 黄色怎么调成土黄色| 热re99久久精品国产66热6| 久久精品国产亚洲av涩爱| av播播在线观看一区| 18禁在线播放成人免费| 91精品国产九色| 国产一区二区三区综合在线观看 | 欧美日韩精品成人综合77777| 国产精品女同一区二区软件| 日本黄色片子视频| 国产一区二区三区综合在线观看 | 777米奇影视久久| 一区二区三区乱码不卡18| 中文字幕制服av| 伊人久久精品亚洲午夜| 视频区图区小说| 欧美zozozo另类| 中文字幕免费在线视频6| 久久鲁丝午夜福利片| 3wmmmm亚洲av在线观看| 国内精品宾馆在线| 人人妻人人看人人澡| 国产精品偷伦视频观看了| 18禁在线播放成人免费| 久久久久久久久久久丰满| 老司机影院毛片| 成人国产av品久久久| 少妇的逼水好多| 亚洲精品乱码久久久久久按摩| 久久久精品欧美日韩精品| 超碰97精品在线观看| 99re6热这里在线精品视频| 精品久久久久久久久亚洲| 亚洲色图综合在线观看| 老女人水多毛片| 免费观看在线日韩| 激情五月婷婷亚洲| 国产91av在线免费观看| av.在线天堂| 大片电影免费在线观看免费| 国产69精品久久久久777片| 麻豆成人av视频| 最新中文字幕久久久久| av国产久精品久网站免费入址| 新久久久久国产一级毛片| 亚洲电影在线观看av| 777米奇影视久久| 国产成人aa在线观看| 人妻一区二区av| 国产一区二区在线观看日韩| 舔av片在线| 久久久久久久精品精品| 看十八女毛片水多多多| 国产男女超爽视频在线观看| 成人亚洲欧美一区二区av| 天天一区二区日本电影三级| 久久女婷五月综合色啪小说 | 午夜免费男女啪啪视频观看| 又黄又爽又刺激的免费视频.| 精品一区二区三卡| 成年女人在线观看亚洲视频 | 亚洲精品国产成人久久av| 一区二区三区免费毛片| 两个人的视频大全免费| 国产精品99久久久久久久久| 国产成人a区在线观看| 欧美3d第一页| 国产中年淑女户外野战色| 狂野欧美激情性bbbbbb| 一级av片app| 亚洲丝袜综合中文字幕| h日本视频在线播放| 日韩av不卡免费在线播放| 国产精品一及| 国产一区二区亚洲精品在线观看| 亚洲天堂av无毛| 1000部很黄的大片| 成年版毛片免费区| 大陆偷拍与自拍| 人妻制服诱惑在线中文字幕| 亚洲经典国产精华液单| 纵有疾风起免费观看全集完整版| 黄片wwwwww| av免费在线看不卡| 国产精品国产三级国产av玫瑰| 少妇被粗大猛烈的视频| 中文字幕av成人在线电影| 日韩人妻高清精品专区| 我的老师免费观看完整版| 大话2 男鬼变身卡| 我的老师免费观看完整版| 在线观看一区二区三区| 综合色av麻豆| 国产精品熟女久久久久浪| 狠狠精品人妻久久久久久综合| 午夜免费男女啪啪视频观看| 三级男女做爰猛烈吃奶摸视频| 欧美xxⅹ黑人| 成人亚洲欧美一区二区av| 国产老妇女一区| 久久久久久久久久人人人人人人| 国内精品美女久久久久久| 精品久久久久久久久亚洲| 国产成人精品婷婷| 又粗又硬又长又爽又黄的视频| 亚洲综合色惰| 美女主播在线视频| 成人黄色视频免费在线看| 亚洲国产av新网站| 国产高清国产精品国产三级 | 国产成人aa在线观看| 国产毛片在线视频| 视频区图区小说| 国内精品宾馆在线| 国产综合懂色| 欧美激情国产日韩精品一区| 特级一级黄色大片| 麻豆成人av视频| 中文资源天堂在线| 99re6热这里在线精品视频| 亚洲av二区三区四区| 一区二区三区精品91| 在线观看三级黄色| 噜噜噜噜噜久久久久久91| 欧美少妇被猛烈插入视频| 最新中文字幕久久久久| 日韩中字成人| 菩萨蛮人人尽说江南好唐韦庄| 男女边吃奶边做爰视频| 亚洲国产精品成人久久小说| 亚洲丝袜综合中文字幕| 99久久九九国产精品国产免费| 人妻制服诱惑在线中文字幕| 亚洲欧美成人精品一区二区| 少妇的逼水好多| 久久久久久久久久久免费av| 国产有黄有色有爽视频| 亚洲精品亚洲一区二区| 亚洲欧洲国产日韩| 中文欧美无线码| 国产精品成人在线| 97在线视频观看| 一区二区三区乱码不卡18| 亚洲欧美中文字幕日韩二区| 国产精品三级大全| 欧美性感艳星| 一级毛片久久久久久久久女| 美女cb高潮喷水在线观看| 国产精品国产三级国产专区5o| 亚洲三级黄色毛片| .国产精品久久| 亚洲av国产av综合av卡| 成年女人看的毛片在线观看| 国产乱人视频| 日韩av不卡免费在线播放| 亚洲美女视频黄频| 听说在线观看完整版免费高清| 黄色视频在线播放观看不卡| 不卡视频在线观看欧美| 精品人妻视频免费看| 国产精品国产三级国产av玫瑰| 麻豆精品久久久久久蜜桃| 久久久久国产精品人妻一区二区| 国产精品久久久久久久久免| 欧美老熟妇乱子伦牲交| 你懂的网址亚洲精品在线观看| 黄色日韩在线| 久久韩国三级中文字幕| 99久国产av精品国产电影| 美女视频免费永久观看网站| 午夜精品一区二区三区免费看| 久久热精品热| 午夜福利视频1000在线观看| 成人一区二区视频在线观看| 夫妻性生交免费视频一级片| 男的添女的下面高潮视频| 婷婷色麻豆天堂久久| 亚洲成人精品中文字幕电影| 精品国产三级普通话版| 亚洲国产欧美人成| 69人妻影院| 免费看av在线观看网站| 国产久久久一区二区三区| 亚洲三级黄色毛片| 亚洲欧美成人精品一区二区| 国产免费福利视频在线观看| 哪个播放器可以免费观看大片| 成人特级av手机在线观看| 观看美女的网站| 狂野欧美白嫩少妇大欣赏| 久久久久久久大尺度免费视频| 18禁动态无遮挡网站| 久久久久久久亚洲中文字幕| 十八禁网站网址无遮挡 | 少妇 在线观看| 国产高清不卡午夜福利| 亚洲色图综合在线观看| 伦精品一区二区三区| 2018国产大陆天天弄谢| 国产 一区 欧美 日韩| 99热国产这里只有精品6| 嫩草影院新地址| 香蕉精品网在线| 麻豆乱淫一区二区| 日日啪夜夜爽| 国产日韩欧美在线精品| 亚洲性久久影院| 久久精品人妻少妇| 18禁在线播放成人免费| 一级爰片在线观看| 永久免费av网站大全| 日日撸夜夜添| 国产男人的电影天堂91| 国产男人的电影天堂91| 亚洲欧洲国产日韩| 久久人人爽人人爽人人片va| 亚洲精品久久午夜乱码| 午夜免费观看性视频| av免费在线看不卡| 别揉我奶头 嗯啊视频| 亚洲欧美日韩无卡精品| a级毛色黄片| 成人毛片a级毛片在线播放| 嫩草影院入口| 天堂网av新在线| 精品人妻偷拍中文字幕| freevideosex欧美| 日本wwww免费看| 国产一区亚洲一区在线观看| 亚洲aⅴ乱码一区二区在线播放| 成人鲁丝片一二三区免费| 欧美激情久久久久久爽电影| 深爱激情五月婷婷| 别揉我奶头 嗯啊视频| 国产有黄有色有爽视频| 国内少妇人妻偷人精品xxx网站| 久久99精品国语久久久| 国产成人免费观看mmmm| 亚洲国产精品成人综合色| 黄色视频在线播放观看不卡| 午夜视频国产福利| 日本午夜av视频| 久久精品夜色国产| 婷婷色麻豆天堂久久| av国产精品久久久久影院| 亚洲不卡免费看| 欧美老熟妇乱子伦牲交| 最近中文字幕2019免费版| 视频区图区小说| 国产精品精品国产色婷婷| 国产老妇伦熟女老妇高清| 日日啪夜夜爽| 免费观看在线日韩| 国产亚洲91精品色在线| 久久综合国产亚洲精品| av女优亚洲男人天堂| 国产v大片淫在线免费观看| 大又大粗又爽又黄少妇毛片口| 大香蕉久久网| 日韩不卡一区二区三区视频在线| 亚洲欧美日韩另类电影网站 | 欧美区成人在线视频| 日韩不卡一区二区三区视频在线| av又黄又爽大尺度在线免费看| 校园人妻丝袜中文字幕| 国产成人aa在线观看| 精品99又大又爽又粗少妇毛片| 天天一区二区日本电影三级| 2021天堂中文幕一二区在线观| 欧美一区二区亚洲| 日日啪夜夜爽| 三级男女做爰猛烈吃奶摸视频| 日韩一区二区视频免费看| 亚洲天堂av无毛| 熟女人妻精品中文字幕| 国产探花极品一区二区| 各种免费的搞黄视频| 又爽又黄无遮挡网站| 午夜福利在线观看免费完整高清在| 成人漫画全彩无遮挡| 97在线人人人人妻| 日本与韩国留学比较| 少妇猛男粗大的猛烈进出视频 | 国产精品福利在线免费观看| a级毛色黄片| 久久久亚洲精品成人影院| 久久久久网色| 男的添女的下面高潮视频| 在线免费十八禁| 在线观看av片永久免费下载| 亚洲欧美一区二区三区国产| 2021天堂中文幕一二区在线观| 久久久久性生活片| 亚洲真实伦在线观看| a级毛色黄片| 只有这里有精品99| 美女视频免费永久观看网站| 老司机影院毛片| 五月玫瑰六月丁香| 人妻系列 视频| 一级黄片播放器| 黄色视频在线播放观看不卡| 日日摸夜夜添夜夜添av毛片| 看非洲黑人一级黄片| 女人十人毛片免费观看3o分钟| 国产亚洲av片在线观看秒播厂| 久久热精品热| 激情五月婷婷亚洲| 亚洲成人久久爱视频| 欧美97在线视频| 在线a可以看的网站| 视频区图区小说| 在现免费观看毛片| av线在线观看网站| 亚洲精品成人久久久久久| 免费看日本二区| 3wmmmm亚洲av在线观看| 亚洲不卡免费看| 亚洲国产精品专区欧美| 亚洲成人久久爱视频| 国产淫语在线视频| 2022亚洲国产成人精品| 国产亚洲午夜精品一区二区久久 | 欧美潮喷喷水| 日本黄大片高清| 亚洲av.av天堂| 国产精品久久久久久av不卡| 免费黄频网站在线观看国产| 国产女主播在线喷水免费视频网站| 亚洲精品国产成人久久av| 国产色婷婷99| av一本久久久久| 欧美老熟妇乱子伦牲交| 最近中文字幕高清免费大全6| 中文字幕人妻熟人妻熟丝袜美| 纵有疾风起免费观看全集完整版| 日韩强制内射视频| 在线免费十八禁| 国产淫语在线视频| 男女边摸边吃奶| 91久久精品国产一区二区三区| 国产v大片淫在线免费观看| tube8黄色片| 亚洲精品乱码久久久v下载方式| 亚洲经典国产精华液单| 亚洲国产av新网站| 老司机影院成人| 看黄色毛片网站| 免费播放大片免费观看视频在线观看| 国产淫片久久久久久久久| 伊人久久精品亚洲午夜| 国产精品人妻久久久影院| 中文字幕亚洲精品专区| 免费观看在线日韩| 我要看日韩黄色一级片| 在线亚洲精品国产二区图片欧美 | 亚洲,欧美,日韩| 一区二区av电影网| 国产精品一区www在线观看| 大香蕉久久网| 欧美性猛交╳xxx乱大交人| 成人一区二区视频在线观看| 日本午夜av视频| 18禁在线无遮挡免费观看视频| 亚洲国产成人一精品久久久| 男的添女的下面高潮视频| 国产精品嫩草影院av在线观看| 欧美日韩综合久久久久久| 校园人妻丝袜中文字幕| 成人欧美大片| 国产成人精品婷婷| 亚洲人成网站高清观看| 精品国产三级普通话版| 日韩一区二区三区影片| 亚洲欧美一区二区三区黑人 | 国精品久久久久久国模美| av专区在线播放| 七月丁香在线播放| 成年免费大片在线观看| 99视频精品全部免费 在线| 夫妻午夜视频| 一级av片app| 国产乱人视频| 亚洲国产最新在线播放| 欧美激情久久久久久爽电影| 国产亚洲最大av| 热re99久久精品国产66热6| 欧美三级亚洲精品| 狂野欧美激情性xxxx在线观看| 精品亚洲乱码少妇综合久久| 久久99精品国语久久久| 男插女下体视频免费在线播放| 国产毛片在线视频| 久久久国产一区二区| 亚洲图色成人| 日韩伦理黄色片| 亚洲av成人精品一区久久| 国产亚洲最大av| 国产精品.久久久| 精品少妇久久久久久888优播| 欧美日本视频| 国产亚洲91精品色在线| 女的被弄到高潮叫床怎么办| 两个人的视频大全免费| 欧美成人a在线观看| 熟女av电影| 七月丁香在线播放| 亚洲欧美日韩另类电影网站 | 成年女人在线观看亚洲视频 | 中文字幕久久专区| 男女下面进入的视频免费午夜| 久久久久国产精品人妻一区二区| 免费看光身美女| 菩萨蛮人人尽说江南好唐韦庄| 麻豆精品久久久久久蜜桃| 黄色欧美视频在线观看| 啦啦啦在线观看免费高清www| 亚洲av成人精品一区久久| 国产视频首页在线观看| 久久午夜福利片| av天堂中文字幕网| 久久国内精品自在自线图片| 欧美xxⅹ黑人| 精品久久久精品久久久| 免费观看av网站的网址| 男人和女人高潮做爰伦理| 日韩欧美 国产精品| 日韩 亚洲 欧美在线| 亚洲精品日韩在线中文字幕| 韩国高清视频一区二区三区| 久久久久精品久久久久真实原创| 国产91av在线免费观看| 国产伦精品一区二区三区视频9| 伊人久久国产一区二区| 久久久精品94久久精品| 久久久久久久亚洲中文字幕| 日本猛色少妇xxxxx猛交久久| 亚洲在久久综合| 69av精品久久久久久| 身体一侧抽搐| 97精品久久久久久久久久精品| 欧美bdsm另类| 日日啪夜夜爽| 国产精品爽爽va在线观看网站| 久久久欧美国产精品| 国产精品女同一区二区软件| 老司机影院毛片| 日本欧美国产在线视频| 黄片无遮挡物在线观看| 久久久久久久国产电影| 国产成人精品婷婷| 亚洲无线观看免费| 久久6这里有精品| 人妻夜夜爽99麻豆av| 蜜桃亚洲精品一区二区三区| 五月玫瑰六月丁香| 人妻一区二区av| 少妇猛男粗大的猛烈进出视频 | 大香蕉97超碰在线| 国产高潮美女av| 国产精品久久久久久久久免| 亚洲av福利一区| 不卡视频在线观看欧美| 亚洲最大成人手机在线| av播播在线观看一区| 国产高清三级在线| 久久人人爽av亚洲精品天堂 | 日韩不卡一区二区三区视频在线| 男女无遮挡免费网站观看| 日本爱情动作片www.在线观看| 亚洲国产高清在线一区二区三| 欧美丝袜亚洲另类| 一区二区av电影网| 日韩欧美一区视频在线观看 | 91久久精品国产一区二区三区| 日日摸夜夜添夜夜添av毛片| 一区二区三区乱码不卡18| 久久久午夜欧美精品| 亚洲av在线观看美女高潮| 嫩草影院精品99| 国产精品一区二区三区四区免费观看| 在线 av 中文字幕| 99热国产这里只有精品6| 午夜免费观看性视频| 亚洲婷婷狠狠爱综合网| 插逼视频在线观看| 国产亚洲一区二区精品| av卡一久久| 国产永久视频网站| 亚洲精品乱码久久久v下载方式| 五月玫瑰六月丁香| 精品酒店卫生间| 国产美女午夜福利| 少妇人妻久久综合中文| kizo精华| 最近手机中文字幕大全| 国产黄色免费在线视频| 国产 一区 欧美 日韩| 精品少妇久久久久久888优播| 人妻夜夜爽99麻豆av| 亚洲精品乱码久久久v下载方式| 狂野欧美激情性xxxx在线观看| 女人久久www免费人成看片| 最近的中文字幕免费完整| av免费观看日本| 国产淫语在线视频| 国产av不卡久久|