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

    Comparative analysis of the arc characteristics inside the convergingdiverging and cylindrical plasma torches

    2020-04-24 05:50:10JianghongSUN孫江宏SurongSUN孫素蓉HaixingWANG王海興ChongNIU牛沖andTaoZHU朱濤
    Plasma Science and Technology 2020年3期

    Jianghong SUN (孫江宏),Surong SUN (孫素蓉),Haixing WANG (王海興),Chong NIU (牛沖) and Tao ZHU (朱濤)

    School of Astronautics,Beijing University of Aeronautics and Astronautics,Beijing 100191,People’s Republic of China

    Abstract

    Keywords:plasma torch,thermodynamic non-equilibrium,arc characteristics

    1.Introduction

    As an important high temperature heat source,non-transferred arc plasma torch has a wide range of applications in material processing,spraying,biomass gasification and waste disposal.In these applications,the torch performance,such as high enthalpy,good stability,and long service life,especially electrodes life,are of primary concerns for optimal design of plasma torch in making these devices commercially viable[1,2].In spite of significant progresses over the past decades,the physical processes inside plasma torch are still not fully understood.This lack of understanding is a consequence of the complexity prevailing in plasma torch caused by the strongly interactions of electrical,magnetic,thermal,and fluid dynamic effects [3–6].

    Traditional direct current (DC) plasma torch has a cylindrical structure as shown in figure 1(a).Over the past several decades,intensive studies have been conducted to investigate the arc behavior in such kinds of torches [7–13].The main advantage of cylindrical torch is that its structure is relatively simple while it also suffers two problems.Firstly,because the essential nature of arc always tends to find the shortest path between the cathode and the anode[14–16],this usually results in a short arc and low arc voltage.Secondly,the gas flow parallels to the anode of cylindrical torch.The arc anode attachment is directly exposed to the action of gas flow and is very susceptible to axial gas dynamic forces.The location of the arc anode attachment is determined by a balance of gas dynamic,electric,and magnetic forces as well as thermal effects [17,18].Under the influence of the superimposed flow the arc attachment at the anode may move upstream or downstream,which leads to that the torch is often operated in a ‘fluctuating mode’ [19–21].The arc root movement inside the torch often relates to the voltage fluctuations and then directly translates into power fluctuations which then affect heating of injected powders in material processing.

    Recently,an improved structure of torch with converging-diverging anode is developed which attempts to overcome these disadvantages[22,23].The basic structure of this DC non-transferred arc plasma torch is schematically shown in figure 1(b) which consists of converging section,constrictor and diverging section.It is expected that the converging-diverging torch offers several advantages over cylindrical plasma torch.Firstly,the arc length would be increased due to the presence of the constrictor,and this leads to the increase of arc voltage and thus the input power and enthalpy of plasma jets also increase.Secondly,due to the higher arc voltage,amplitude of the voltage fluctuations becomes small relative to the total arc voltage.Thirdly,it has been found from experiment that in the case of a specific combination of operating parameters,the arc passes through the constrictor and attaches to the expanded section of the anode in a diffusion mode with low current density [24].Although intensive experimental and analytical studies have been carried out for this converging-diverging torch[22,23],the evolution and development of arcs in this kind of torch is still unclear.Since the experimental diagnosis inside the plasma torch is quite difficult,it is necessary to obtain the characteristics of plasma flow in the converging-diverging torch by numerical simulation analysis,especially investigating the effect of the constrictor on the arc characteristics.

    The purpose of this paper is to obtain a deep insight into the effect of torch structure on arc characteristics through comparative analysis of the two kinds of plasma torch.In order to avoid specifying the lowest conductivity for conducting current on the anode wall of the plasma torch under local thermal equilibrium (LTE) assumption,a two-temperature thermodynamic non-equilibrium method is used in this numerical simulation.Section 2 presents the physical and mathematic model based on two-temperature chemical equilibrium assumption,the computational domain and boundary conditions.Numerical results of the modeling of the arc behavior inside two kinds of structure plasma torches,including comparisons with experimental measurements,are shown in section 3.Conclusions are drawn in section 4.

    2.Modeling approach

    2.1.Computation domain

    The numerical simulation of this paper is based on two structures of plasma torches.The first is a conventional cylindrical torch which was numerically simulated in [3].In this study,it is used as the baseline torch.As in [3],the simplified structure of cylindrical torch used in the simulation is shown in figure 2(a).It is noted that only the half of the torch geometry is presented since axial symmetry is assumed.The torch has a length of 32 mm (GF) and a radius of 6.35 mm (EF).The cathode is assumed as a cylindrical rod with a radius of 3 mm.The cathode current density is assumed to be uniformly distributed on cathode spot (CD)with a radius of 0.5 mm.For the rest of cathode plane(BC),it is assumed that there is no current flow through this area.The second is a converging-diverging torch which was reported in[22,23].A schematic sketch of the torch is shown in figure 2(b),the length of the torch is 62.82 mm.The inlet JA is 6.65 mm in radius with a length of 2 mm(JI).The converging section(IH)of the torch is 10 mm in length with a half angle of 28°.The diverging section of the anode(FG)is conical with a half-angle of 7°,and the nozzle-exit radius(EF)is 7.1 mm.The constrictor (HG) is 3.22 mm in length and 1.4 mm in radius.The cathode(ABC)is a cylinder rod with a radius of 3.54 mm,with a flat tip (CD) of a radius of 0.5 mm.

    2.2.Basic assumptions

    The following assumptions are used in this simulation:(1)The plasma is assumed to be axially symmetric,laminar,and steady.(2) The plasma is considered in chemical equilibrium and thermodynamic non-equilibrium,i.e.the electron temperature Teand heavy particle temperature Thcan be different and are separately solved by each energy equation.Maxwellian distributions are assumed for both electrons and heavy particles.(3) Thermodynamic and transport data are determined by Teand Th,which can be calculated by the Chapman–Enskog theory and detailed procedure is described in [25,26].(4) The plasma is assumed to be optically thin,and radiative losses are treated by the net emission coefficient method.(5) Plasma sheath is not included.

    Figure 2.Computational domain of torches with (a) the cylindrical anode and (b) the converging-diverging anode (unit:mm).

    Table 1.Boundary conditions for two-temperature model.

    The laminar flow assumption adopted in this simulation is due to the low Reynolds number inside this convergingdiverging torch.According to the geometrical and operating parameters of the plasma torch,the estimated Reynolds number is on the order of 1000.The estimated Reynolds number for cylindrical torch under the flow conditions used in this study is around 200,which is on the same order with the estimation of Westhoff[3].For this condition,the assumption of laminar is reasonable for cylindrical torch.In previous experimental study of plasma torch [22,27,28]the uniform and diffused arc attachments at the anode have been observed with similar operating conditions used in this study.This infers that the arc can be considered as axially symmetric for this case.Although there are small fluctuations in the voltage under the present operational conditions,the measured voltage and electro-thermal efficiency are time averaged and the same input current and gas flow rate are used for simulations.Hence the steady state is assumed in the present model.

    It is noted that two-temperature thermodynamic assumption is used in this study.It is because that in the LTE model,the electron temperature is equal to the heavy particle temperature,which is low near the electrodes due to the strong cooling they experience.This causes the low value of the equilibrium electrical conductivity near the electrodes which is not enough to support electrical current flow to the anode.In order to solve this problem,numerical simulation methods based on LTE assumptions usually use various artificial methods to improve conductivity [3,7,8].One of main advantages of two-temperature model is that it does not require any special treatment of the electrical conductivity in cold boundary layer,because the electron temperature will typically remain high in the arc attachment regions.

    2.3.Boundary conditions

    Figure 3.The predicted temperature distribution (lines) with reported results from [3],I=250 A,Q=9.83 slm.

    The boundary conditions for different variables are shown in table 1.The inlet axial velocity is assumed to be constant which is determined by gas flow Q and inlet area A and the radial velocity at inlet is zero.The inlet temperatures for electrons and heavy particles are equal to 1000 K.Zero electron temperature gradient condition is defined in other boundaries.On the axis of symmetry DE,all variables have zero gradient.At the torch outlet,the pressure is set to be 1 atm.For cathode,a linear rising of heavy particles temperature from 1000 K to 3500 K is set along ABC,while at the cathode tip(CD)temperature for the heavy particle is set to be 3500 K.For the sake of comparison,the same constant current density of 3×107A m-2is given at flat tips of cathode for the two torches,which was also used in[3].Zero current density is set on the rest part of the cathode.For the anode boundary,the electric potential is set to be zero.These boundary conditions are typical for plasma torch simulations and the corresponding discussion was presented in detail in[29].For the sake of comparison in the numerical simulation,we use the same gas flow,operating current and cathode tip current density for two kinds of torches.

    2.4.Governing equations

    The governing equations used in this simulation consist of mass continuity,momentum,electron and heavy particle energy conversation equations and Maxwell equations.The electromagnetic field is closely coupled by the Lorentz force and the Ohmic heating with the momentum equation and the energy equation.Mass continuity:

    where ρ represents the mass density and u is the flow velocity Momentum conservation:

    Energy conservation for electrons:

    where necorresponds to the electron number density.kBis the Boltzmann constant.represents a particle charged ζ times.represents the ionization energy of this particle,κethe electron conductivity,κrthe reaction thermal conductivity,the enthalpy work of electrons,Qjthe Joule heating,Eechthe electron-heavy-particle energy exchange term,Qrthe volumetric net radiation losses.

    Energy conservation for heavy particles:

    where nhis the heavy particle density,κhthe heavy particle conductivity.

    The magnetic vector potential A is obtained from the Poisson equation:

    whereμ0denotes the magnetic permeability,j the current density.The magnetic field B which causes the Lorentz force FLcan be derived fromB=?×A.

    Current continuity equation:

    Equation (6) together with Ohm’s law (j=σE) and electric potential equation (E=-?V) are solved,where E is the electric field,σ the electrical conductivity,V the electric potential.

    3.Results

    3.1.Code validation

    The numerical code used in this paper is first used for numerical simulation of the baseline torch to verify the reliability of the code.The cylindrical torch shown in figure 2(a),was previously simulated to investigate the heat transfer,and fluid flow phenomena inside plasma torches by Westhoff’s [3]based on local thermodynamic equilibrium assumption.In order to treat the problem of current conduction in the low temperature region of the anode boundary layer,they adopted an artificially given dependence function of electrical conductivity and temperature below 9000 K.For temperatures above 9000 K the conductivity is equal to the calculated(LTE) values.The arc current and gas flow rate in [3]were respectively set to be 250 A and 9.83 slm.

    In order to test our modeling approach,the same geometry,operating parameters of the plasma torch and also based on LTE assumption with[3]are used in this simulation.The comparison between our calculation results of temperature distribution inside plasma torch and the correspondingresults of [3]is shown in figure 3 and comparison of some important typical parameters is given in table 2.As can be seen from figure 3 and table 2,the numerical predictions are in good agreement with Westhoff’s work [3].

    In order to further verify the numerical simulation approach used in this paper,the code is extended to simulate the converging-diverging plasma torch based on two-temperature model.Comparison of the predicted axial temperature distribution with the experimental measurement results is shown in figure 4.In the experiment,the converging-diverging torch was operated at arc current of 90 A and argon flow rate of 9.8 slm.Spectral measurement of atomic excitation temperature was conducted with a similar method used in[30].Several 0.9 mm diameter holes were drilled into the diverging section of the plasma torch to provide side-on optical access to the internal flow.Axial distribution of atomic excitation temperature of argon was determined from emission line ratios.As can be seen from the figure 4,the numerical simulation results are consistent with the experimental measurements.It is noted that both the electron and heavy particle temperatures obtained from the two-temperature numerical simulation are plotted in the figure 4,and they are very close to each other,indicating that the plasma is close to the thermodynamic equilibrium along the torch axis.

    The measured axial temperatures of mixture of Ar–N2with the converging-diverging torch are presented in figure 5.Similar to the experimental measurement of excitation temperature of pure argon,the axial variation of excitation temperature of argon in Ar–N2mixture was also obtained by relative line intensity ratio method.As can be seen from figure 5 that the addition of nitrogen to the mixture leads to a significant increase in the excitation temperature along the axis of torch.It is mainly due to that adding nitrogen to the mixed gas can significantly increase the specific enthalpy of the gas [31,32].In order to further verify the code,we also performed a numerical simulation using an Ar–N2mixture as the working gas.The calculated axial distribution of excitation temperature of Ar–N2mixture is also presented in figure 5.Considering the uncertainties (~25% [33]) of relative line intensity method used in the experiment,the measured and predicted temperature can be considered to be reasonably agreed.

    Figure 4.Comparison of the predicted axial temperatures with experimental measurement in the converging-diverging torch(argon,90 A,9.8 slm).

    Figure 5.Comparison of the predicted axial temperatures with experimental measurement in the converging-diverging torch (Ar–N2 (65.2% N2),100 A,16.1 slm).

    Figures 6 and 7 show the arc attachment conditions on anode surface inside the converging-diverging torch which were taken with digital camera with exposure time of 250 μs and further ICCD pictures also show the same arc root attachment tendency with the exposure time of 5 μs.A specially designed copper mirror with a boron nitride film coated on its surface central-region is employed.The copper mirror was put at the downstream of the torch exit to transmit the upward-view image of the arc plasma torch into the intensified charge coupled device(ICCD)camera[22,24,25].It can be seen from figures 6 and 7 that the arcs are attached to the anode surface in a diffused and uniform mode in the converging-diverging torch for pure argon and even for the mixture of argon and nitrogen.The experimentally observed arc attachment modes demonstrate the reasonability of axisymmetric assumption in the simulation.

    Figure 6.The image of the diffused arc attachment for converging-diverging torch with gas flow rates of(a)7,(b)11.2,(c) 14 slm of pure argon,arc current of 80 A.

    Figure 7.The image of the diffused arc attachment for converging-diverging torch with the mixture of 2.8 slm argon and 12 slm nitrogen,arc current of (a) 90,(b) 100,(c) 110 A.

    Figure 8.The electron and heavy particle temperature distributions of (a) the cylindrical torch and (b) converging-diverging torch.

    Figure 9.Calculated velocity and pressure distributions of the cylindrical torch (a),(b) and converging-diverging torch (c),(d).

    3.2.Temperature and flow fields

    In order to facilitate comparison,the same operating parameters,arc current of 90 A and gas flow rates of 9.8 slm of pure Ar,are chosen for conventional cylindrical torch and converging-diverging torch in the following simulations.The calculated electron temperature and heavy particle temperature inside the cylindrical torch and converging-diverging torch are given in figure 8.It can be seen that the temperature distributions share two features in common.Firstly,in the central region of the torch,the electron temperature is close to the heavy particle temperature,which infers that the plasma is close to thermodynamic equilibrium state.However,in arc fringes,the electron temperature is higher than the heavy particle temperature,and this indicates that the plasma deviates from the thermodynamic equilibrium state.Secondly,the highest temperatures in the two kinds of torches are almost equal and appear downstream of the cathode tip,indicating that the maximum temperature is primarily determined by the current density distribution of the cathode tip.An important difference between the two temperature distributions is that the length of the high temperature zone in the convergingdiverging torch is much longer than that of cylindrical torch.It is further noted that the length of the converging-diverging torch is approximately 63 mm,which is much larger than 32 mm of the cylindrical torch.It demonstrates that a longer high temperature arc column can be formed inside the converging-diverging torch.The numerical results show that the predicted arc voltage of the converging-diverging torch is 24 V,which is much higher than the 7 V of the cylindrical torch.In the case of the same arc current and gas flow rate,the input power of the converging-diverging torch is much higher than that of the cylindrical torch.

    The velocity and pressure distributions in the cylindrical torch and the converging-diverging torch are respectively given in figures 9(a)–(d).It can be seen from the figure 9 that the maximum gas velocity in the converging-diverging torch is 1400 m s-1,which is much higher than 40 m s-1of the cylindrical torch.The main reason for this difference in gas velocities is caused by the pressure difference between the inlet and outlet of the two structural torches.For the converging-diverging torch,the decrease of gas flow cross section area at converging section and the choking effects caused by arc heating inside the constrictor are the main reasons of the increase in the pressure difference between the inlet and outlet of the torch [34].It can be seen that the pressure difference of converging-diverging torch is on the order of 22 000 Pa,while the pressure difference of the cylindrical torch is only around 100 Pa.In the convergingdiverging torch,the higher gas velocity acts to stretch the arc,increase the arc voltage and the power input.

    Figure 10.Axial temperature (a) and velocity (b) distributions of two different torches.

    Figure 10 shows the velocity distribution and temperature distribution along the axis of the two kinds of torches.It should be noted that in order to facilitate comparison,the axial distances are normalized to their respective length of torch.It is seen that the axial plasma temperature in the convergingdiverging torch decays much slower than that of the cylindrical torch even though the actual length of the converging-diverging torch is much longer.As can be seen from figure 10(b),the positions at which the highest speeds appear inside the two kinds of torches are different.For cylindrical torch,the maximum of gas velocity occurs near the downstream of the cathode tip,which is primarily driven by j×B pumping effect(j is the current density and B the magnetic induction due to the arc current) in the vicinity of the cathode.However,the maximal velocity inside converging-diverging torch appears around the exit of constrictor.It is because that the gas flow in the constrictor is still subsonic,arc heating continues to accelerate the gas,reaching a maximum at the beginning of the diverging section of the torch.The predicted highest Mach number of approximately 0.56 also appears near the exit of the constrictor infers that the flow velocity in diverging-converging torch is lower than the speed of sound.Since the plasma torch is operated at a back pressure of 1 atm,the gas cannot be fully expanded in the diverging section of the torch unlike the arcjet thruster operated under vacuum conditions[35–41],so both the gas velocity and the Mach number in the diverging section of the torch decrease as the axial distance increases.

    3.3.Arc characteristics

    For the converging-diverging torch,the evolution and development of the arc between the cathode tip and the constrictor exit are critical,so comparison and analysis of arc characteristics of the two kinds of torches are focused in this region.Figure 11 presents the radial distribution of electron,heavy particle temperature and degree of thermodynamic non-equilibrium(Te/Th)at the axial distance 0.5,2 and 6 mm downstream of cathode tip,which corresponds to axial position of the near cathode,the inlet and the outlet of constrictor for the converging-diverging torch.It is noted here that in order to facilitate comparison the radial distances are also normalized according to the respective arc channel radius.

    As can be seen from the figures 11(a) and (c),for the cylindrical torch,the thermodynamic equilibrium area of the arc occupies a small proportion of the channel of torch,which is about 0.3 of the dimensionless radius.In most of the radial region of the torch channel,the plasma deviates significantly from the thermodynamic equilibrium state.As the axial distance increases from 0.2 to 6 mm,the variation of arc column diameter is small since the channel diameter does not change.

    For the converging-diverging torch,the evolution characteristics of the arc before and after the constrictor are different from that of cylindrical torch.It can be seen from the figures 11(b)and(d),as the axial distance increases from 0.2 to 6 mm,the heavy particles and electron temperatures increase across the entire section of channel.It is shown by the dot–dash lines in figures 11(a) and (b),if we define a heavy particles temperature of 6000 K as the boundary of the arc column [34],for the cylindrical torch,the dimensionless radius of the arc column increases from 0.4 to 0.55.However,for the converging-diverging torch the dimensionless radius of the arc column increases significantly,which increases from 0.4 to 0.75.On the other hand,the dimensionless radial thermodynamic equilibrium region of the arc within the converging-diverging torch is also significantly increased relative to that of the cylindrical torch.Similar phenomenon has also been found in the numerical simulation of an arcjet thruster with a converging-diverging nozzle [34]that the radius of the arc core in the constrictor also increases with increasing distance.It can be seen from figure 2 that the radius of the constrictor is only 1.4 mm,which is much smaller than 6.35 mm,the channel radius of the cylindrical torch.The gas is compressed in the constrictor due to the limitation of cold wall,which enhances the mixing of the low temperature gas around the arc column and the high temperature core region.This increases the average temperature of the gas across the constrictor section and improves the degree of thermodynamic equilibrium of the plasma.

    Figure 11.Normalized radial distribution of Te,Th(a),(b)and the non-equilibrium temperature ratio(c),(d)at different distances away from the cathode tip for the cylindrical (a),(c) and converging-diverging torch (b),(d).

    3.4.Analysis of radial momentum balance

    The evolution and development of plasma arc in the torch are determined by the interaction of different forces.Therefore,analyzing the radial distribution of different forces in the developing stage of the arc can deepen our understanding of the influence of torch structure on the arc characteristics[42].Figure 12 presents the normalized radial profiles of different forces acting on the arcs for two kinds of torches.It can be seen from figures 12(a)–(c),for the cylindrical torch,the forces that dominate the radial evolution of the arc are radial gas dynamic forces and Lorentz forces which are in the opposite direction.With the increase of the axial distance from 0.2 to 6 mm,the radial pressure gradient and arc current density decrease with the arc development.This results in the decrease of peak values of radial gas dynamics force and Lorentz force.

    The radial evolution of forces in the converging-diverging torch is significantly different from that in the cylindrical torch.At the downstream of cathode tip,i.e.axial distance of 0.5 mm,the radial gas dynamic force is mainly balanced with the Lorentz force when the dimensionless radius is small,i.e.corresponding to the arc column region.While for the region with dimensionless radius larger than 0.3,the radial gas dynamic force becomes negative which also acts as to compress gas.This is due to the fact that in the converging section of the torch,most of the gas has to pass through the low temperature gas region in arc fringes,thereby increasing the gas pressure gradient in this region and creating a large radial force pointing to arc center.At the inlet of constrictor,i.e.axial distance of 2 mm,the radial gas dynamic and Lorentz force forces are still balanced in the arc column region.In the arc fringes,the radius of the converging section decreases to the radius of constrictor and this leads to the significant increase of the radial gas dynamic force that points to arc center.This radial gas dynamic forces directing toward the arc center in the converging section and the constrictor serve to enhance the mixing of the gas in the high temperature zone of the arc column with the gas in the cold boundary layer.At the axial distance of 6 mm,i.e.the exit of constrictor and also the upstream of diverging section of the torch,the dominant force is the gas dynamic force that causes the arc to expand radially,and the Lorentz force is almost negligible due to the decrease in current density.This facilitates the arc attachment to the anode in this area.

    The experimental observation of arc root attachment for this converging-diverging torch was reported in [22].The pictures recorded by ICCD camera show that in most situations the mixture of nitrogen–argon arc attaches at the downstream of constrictor,i.e.upstream of diverging section of torch.The experimentally observed arc attachment area is qualitatively consistent with this numerical simulation results.It should be mentioned,however,that it is very difficult to accurately predict the location and area of the arc attachment at anode.Arc attachment not only depends on the evolution of radial forces but also on the balance of axial gas dynamic forces and Lorentz forces at the arc attachment region.In arc attachment region,the plasma deviates from both the thermodynamic and chemical equilibrium simultaneously,making the prediction of current density distribution very difficult.There is a need to develop a reasonable model to describe chemical non-equilibrium and species transport process that can accurately capture arc fringe characteristics.

    4.Conclusions

    Arc characteristics inside the conventional cylindrical torch and converging-diverging torch are investigated in this study.In order to avoid the use of artificially given conductivity near the surface of the anode,a two-temperature model is used in the numerical simulation.The reliability of the code is verified by comparing the predicted results with the numerical simulation results of the baseline torch and with the measured axial distribution of temperature in the converging-diverging torch.

    It is found that that due to the presence of the constrictor,the arc characteristics in the converging-diverging torch are quite different from those in the cylindrical torch.In the converging-diverging torch,the pressure difference between the inlet and outlet is much larger than that of cylindrical torch due to the reduction of the flow cross-sectional area and the thermal choking effects of arc heating inside constrictor.Therefore a much higher plasma flow velocity and relatively longer high temperature region can be generated and maintained inside the converging-diverging torch.Correspondingly,the higher gas velocity acts to stretch the arc,increase the arc voltage and the power input.

    The other important role of the constrictor is to enhance the mixing of high temperature gas in the arc column and the low temperature gas in the boundary layer,and to increase the average temperature of the channel section.The modeling results show that in the constrictor,the radius of the arc column increases with the increase of the axial distance,and the thermodynamic equilibrium degree is also improved.This corresponds to the decrease of thickness of the cold boundary layer in the constrictor,which is beneficial to the arc attaching to the anode in the diverging section.

    The momentum balance analysis before and after the constrictor allows us to have a deeper understanding of the evolution process of the arc.For converging section of torch,from the downstream of the cathode tip to the inlet of the constrictor,there exists a strong radial gas dynamic force which points to arc center.At the exit of constrictor and also the diverging section of the torch,the dominant force is the aerodynamic force that causes the arc to expand radially and this is also very beneficial for the arc attaching to this area with diffusive mode.

    Acknowledgments

    This work was supported by National Natural Science Foundation of China (Nos.11575273,11735004,11575019)and the National Postdoctoral Program for Innovative Talents(BX20180029).The authors especially thank Pan Wenxia for providing the experimental data and photos.

    日韩av不卡免费在线播放| av在线app专区| 亚洲不卡免费看| a级片在线免费高清观看视频| 亚洲精品aⅴ在线观看| 国产伦精品一区二区三区视频9| 男女国产视频网站| 热99国产精品久久久久久7| 久久久久久久国产电影| 伊人久久精品亚洲午夜| 日本av免费视频播放| 亚洲精品国产色婷婷电影| 国产精品嫩草影院av在线观看| 国产日韩欧美亚洲二区| .国产精品久久| 最新中文字幕久久久久| 免费观看无遮挡的男女| 亚洲av不卡在线观看| 麻豆成人av视频| 免费高清在线观看视频在线观看| 国产精品女同一区二区软件| 亚洲av综合色区一区| 免费看不卡的av| 美女大奶头黄色视频| 亚洲一级一片aⅴ在线观看| 91久久精品国产一区二区成人| 国产乱来视频区| 午夜视频国产福利| 交换朋友夫妻互换小说| 高清av免费在线| 老司机亚洲免费影院| 欧美精品高潮呻吟av久久| 九九在线视频观看精品| 女人精品久久久久毛片| 午夜福利视频精品| 亚洲不卡免费看| 色94色欧美一区二区| 国产一区二区在线观看日韩| 91精品伊人久久大香线蕉| 欧美 日韩 精品 国产| 久久精品国产亚洲av涩爱| 日本与韩国留学比较| 久久影院123| 亚洲精品成人av观看孕妇| 两个人免费观看高清视频| 天天操日日干夜夜撸| 国产成人精品婷婷| 黄色视频在线播放观看不卡| 久久久欧美国产精品| 日韩av在线免费看完整版不卡| 午夜久久久在线观看| 老司机影院成人| 国产免费又黄又爽又色| 久久99一区二区三区| 亚洲欧美清纯卡通| 精品亚洲成a人片在线观看| 亚洲成人av在线免费| 九九久久精品国产亚洲av麻豆| av卡一久久| 高清毛片免费看| 久久久久久久久久久免费av| 日韩不卡一区二区三区视频在线| 亚洲欧美清纯卡通| 美女中出高潮动态图| 在线看a的网站| 高清av免费在线| 欧美精品亚洲一区二区| 亚洲精华国产精华液的使用体验| 如日韩欧美国产精品一区二区三区 | 国产精品免费大片| 最近的中文字幕免费完整| av视频免费观看在线观看| 日韩一本色道免费dvd| 丰满乱子伦码专区| 麻豆精品久久久久久蜜桃| 尾随美女入室| 久久精品久久久久久噜噜老黄| 91久久精品国产一区二区成人| 亚洲人成网站在线播| 国产免费现黄频在线看| 一区二区av电影网| av电影中文网址| 亚洲国产毛片av蜜桃av| 中文天堂在线官网| 中文字幕最新亚洲高清| 另类精品久久| 亚洲人成77777在线视频| 日日撸夜夜添| 大香蕉97超碰在线| 精品人妻在线不人妻| 汤姆久久久久久久影院中文字幕| 最近中文字幕高清免费大全6| 丝袜在线中文字幕| 人人妻人人澡人人爽人人夜夜| 亚洲伊人久久精品综合| 日本欧美国产在线视频| 午夜视频国产福利| 免费观看a级毛片全部| 亚洲欧美日韩卡通动漫| 亚洲少妇的诱惑av| 我的老师免费观看完整版| 夫妻午夜视频| 欧美 日韩 精品 国产| 久久精品人人爽人人爽视色| 自拍欧美九色日韩亚洲蝌蚪91| 欧美日韩av久久| 亚洲精华国产精华液的使用体验| 人人妻人人澡人人看| 最后的刺客免费高清国语| 一个人免费看片子| a级毛片在线看网站| 久久午夜综合久久蜜桃| 三上悠亚av全集在线观看| 国产精品一区二区在线观看99| av免费在线看不卡| 亚洲欧美一区二区三区国产| 黄色配什么色好看| 精品一区二区免费观看| 亚洲国产色片| 欧美日本中文国产一区发布| 久久久久久伊人网av| 国产精品嫩草影院av在线观看| 80岁老熟妇乱子伦牲交| 亚洲熟女精品中文字幕| 午夜福利影视在线免费观看| 国产不卡av网站在线观看| 国产黄频视频在线观看| 99久久综合免费| 少妇精品久久久久久久| 婷婷色麻豆天堂久久| 综合色丁香网| 国产色爽女视频免费观看| 久久久久久久久久久久大奶| 啦啦啦视频在线资源免费观看| 女的被弄到高潮叫床怎么办| 国产精品麻豆人妻色哟哟久久| 人妻一区二区av| 高清黄色对白视频在线免费看| 我的女老师完整版在线观看| 亚洲综合精品二区| 久久久精品区二区三区| 免费看av在线观看网站| av线在线观看网站| 99国产精品免费福利视频| 国产极品天堂在线| 中国三级夫妇交换| 久久影院123| 一区二区三区免费毛片| 欧美亚洲日本最大视频资源| 女的被弄到高潮叫床怎么办| 久久久a久久爽久久v久久| 国产午夜精品一二区理论片| 精品酒店卫生间| 亚洲欧美精品自产自拍| 啦啦啦在线观看免费高清www| 最近中文字幕高清免费大全6| 一级黄片播放器| 晚上一个人看的免费电影| 国产精品蜜桃在线观看| 国产av一区二区精品久久| 在线天堂最新版资源| 亚洲激情五月婷婷啪啪| 一边亲一边摸免费视频| 日韩欧美精品免费久久| 麻豆乱淫一区二区| 日韩电影二区| 夫妻性生交免费视频一级片| 成人手机av| 亚洲综合色网址| 一级毛片我不卡| 国产视频内射| 国产av一区二区精品久久| 久久久久国产网址| 国产精品成人在线| 爱豆传媒免费全集在线观看| 黄片播放在线免费| 亚洲精品成人av观看孕妇| 大香蕉久久网| 久久久精品94久久精品| 18禁观看日本| 18在线观看网站| 亚洲成人手机| 午夜免费男女啪啪视频观看| 最近最新中文字幕免费大全7| 日韩 亚洲 欧美在线| 日韩,欧美,国产一区二区三区| 桃花免费在线播放| 最近中文字幕2019免费版| 日韩中文字幕视频在线看片| 特大巨黑吊av在线直播| av网站免费在线观看视频| xxxhd国产人妻xxx| 欧美一级a爱片免费观看看| 国产一区二区三区av在线| 插阴视频在线观看视频| 国产成人一区二区在线| 日日撸夜夜添| 中文字幕亚洲精品专区| 伊人亚洲综合成人网| 久久国产亚洲av麻豆专区| 美女中出高潮动态图| 国产熟女午夜一区二区三区 | 寂寞人妻少妇视频99o| 美女xxoo啪啪120秒动态图| 国产亚洲精品第一综合不卡 | 欧美日韩视频高清一区二区三区二| 日日啪夜夜爽| 日本午夜av视频| 热99国产精品久久久久久7| 亚洲在久久综合| 一级毛片 在线播放| 日韩一区二区三区影片| 日本-黄色视频高清免费观看| 国产精品一区二区在线不卡| 欧美97在线视频| 一区在线观看完整版| 我要看黄色一级片免费的| 伦精品一区二区三区| 一本—道久久a久久精品蜜桃钙片| 亚洲欧美日韩另类电影网站| 最近中文字幕2019免费版| 国产日韩一区二区三区精品不卡 | 久久精品久久久久久噜噜老黄| 人妻少妇偷人精品九色| 免费av不卡在线播放| 国产片内射在线| av在线app专区| 国产精品成人在线| 中文乱码字字幕精品一区二区三区| 十八禁网站网址无遮挡| 欧美另类一区| 天天影视国产精品| 久久久久精品久久久久真实原创| 啦啦啦啦在线视频资源| 亚洲精品aⅴ在线观看| 国产成人a∨麻豆精品| 亚洲人成网站在线观看播放| 又粗又硬又长又爽又黄的视频| 18禁动态无遮挡网站| a级毛片在线看网站| 午夜91福利影院| 国产精品人妻久久久久久| 色婷婷av一区二区三区视频| 中文字幕亚洲精品专区| 日产精品乱码卡一卡2卡三| 日韩免费高清中文字幕av| 人妻一区二区av| 午夜激情久久久久久久| 十八禁高潮呻吟视频| 大香蕉久久成人网| a 毛片基地| 人人妻人人澡人人爽人人夜夜| 99九九线精品视频在线观看视频| 亚洲精品久久成人aⅴ小说 | 国产在线一区二区三区精| 亚洲精品,欧美精品| 老司机影院成人| 成人免费观看视频高清| 国产成人av激情在线播放 | 国产一区二区在线观看av| 天堂俺去俺来也www色官网| 永久网站在线| 久久午夜福利片| 欧美激情 高清一区二区三区| 在线观看三级黄色| 十八禁高潮呻吟视频| 日本av免费视频播放| 简卡轻食公司| 高清黄色对白视频在线免费看| av卡一久久| 高清视频免费观看一区二区| 制服人妻中文乱码| 一区二区日韩欧美中文字幕 | 伦理电影免费视频| 日韩精品免费视频一区二区三区 | 97在线视频观看| 日日啪夜夜爽| 亚洲国产av新网站| 18禁观看日本| 午夜av观看不卡| 亚洲欧美一区二区三区国产| 91久久精品电影网| 一级毛片我不卡| 亚洲精品一二三| 精品国产乱码久久久久久小说| av国产精品久久久久影院| 久久久久久久国产电影| 亚洲av二区三区四区| 我的女老师完整版在线观看| 99久久精品国产国产毛片| 国产日韩欧美亚洲二区| 岛国毛片在线播放| 夫妻午夜视频| 秋霞伦理黄片| 一边亲一边摸免费视频| 久久国产精品男人的天堂亚洲 | 午夜福利视频精品| 久久久久国产网址| 黄色怎么调成土黄色| 亚洲人成网站在线观看播放| 久久久精品94久久精品| 国产淫语在线视频| 亚洲av福利一区| 亚洲av在线观看美女高潮| av免费观看日本| 丰满饥渴人妻一区二区三| 在线看a的网站| av电影中文网址| 色94色欧美一区二区| 日韩成人av中文字幕在线观看| av视频免费观看在线观看| 午夜激情久久久久久久| 成人综合一区亚洲| 国产一区亚洲一区在线观看| av在线老鸭窝| 亚州av有码| 亚洲国产精品国产精品| 国产国语露脸激情在线看| 亚洲精品456在线播放app| 2022亚洲国产成人精品| 97精品久久久久久久久久精品| 国产免费福利视频在线观看| 国语对白做爰xxxⅹ性视频网站| 五月开心婷婷网| 黄色怎么调成土黄色| 秋霞伦理黄片| 男女国产视频网站| 免费看不卡的av| 黑人猛操日本美女一级片| 精品视频人人做人人爽| 日本猛色少妇xxxxx猛交久久| 亚洲欧美色中文字幕在线| 色婷婷av一区二区三区视频| 国产亚洲精品第一综合不卡 | 91精品伊人久久大香线蕉| 日韩,欧美,国产一区二区三区| 欧美性感艳星| 国产精品99久久久久久久久| 日本av手机在线免费观看| 亚洲精品一二三| 亚洲精品国产色婷婷电影| 日韩一区二区视频免费看| 自拍欧美九色日韩亚洲蝌蚪91| 国产高清三级在线| 中文精品一卡2卡3卡4更新| 午夜福利视频在线观看免费| 大码成人一级视频| 精品一区二区三卡| 国产国拍精品亚洲av在线观看| 亚洲av成人精品一区久久| 十分钟在线观看高清视频www| 久久精品国产鲁丝片午夜精品| 久久精品国产亚洲av涩爱| 一级毛片电影观看| 80岁老熟妇乱子伦牲交| 夫妻午夜视频| 91精品一卡2卡3卡4卡| 亚洲精品国产色婷婷电影| 午夜免费观看性视频| 国产精品免费大片| 一区二区日韩欧美中文字幕 | 国产视频首页在线观看| 午夜免费男女啪啪视频观看| 狂野欧美激情性xxxx在线观看| 久久毛片免费看一区二区三区| 欧美另类一区| 国产毛片在线视频| a级毛片免费高清观看在线播放| 亚洲三级黄色毛片| 高清av免费在线| 亚洲国产毛片av蜜桃av| 爱豆传媒免费全集在线观看| 又黄又爽又刺激的免费视频.| 成人国语在线视频| 美女中出高潮动态图| 日韩av免费高清视频| 成人综合一区亚洲| 男女国产视频网站| 欧美激情国产日韩精品一区| 日日爽夜夜爽网站| 欧美日韩综合久久久久久| av国产久精品久网站免费入址| 色网站视频免费| 插阴视频在线观看视频| 成人毛片60女人毛片免费| 曰老女人黄片| 18在线观看网站| 国产精品一二三区在线看| 欧美日韩成人在线一区二区| 久热久热在线精品观看| 亚洲av成人精品一二三区| 久久精品国产鲁丝片午夜精品| 久久久久久久久大av| 乱码一卡2卡4卡精品| 我的女老师完整版在线观看| 美女中出高潮动态图| 国国产精品蜜臀av免费| 久久亚洲国产成人精品v| 成人漫画全彩无遮挡| 亚洲在久久综合| 亚洲精品av麻豆狂野| 女性生殖器流出的白浆| 午夜福利影视在线免费观看| 成人毛片a级毛片在线播放| 狂野欧美激情性bbbbbb| 日本av免费视频播放| 晚上一个人看的免费电影| 成人国产麻豆网| 亚洲国产精品国产精品| 自拍欧美九色日韩亚洲蝌蚪91| h视频一区二区三区| 亚洲欧洲精品一区二区精品久久久 | 国产免费福利视频在线观看| 狂野欧美白嫩少妇大欣赏| 超色免费av| 成人二区视频| 日韩精品有码人妻一区| 午夜福利视频在线观看免费| 国产精品欧美亚洲77777| kizo精华| 啦啦啦中文免费视频观看日本| 亚洲中文av在线| 成人无遮挡网站| 天堂中文最新版在线下载| 午夜激情av网站| 国产黄色免费在线视频| 亚洲av不卡在线观看| www.色视频.com| 天天影视国产精品| 久久久久人妻精品一区果冻| 看免费成人av毛片| 婷婷色综合www| 边亲边吃奶的免费视频| 欧美3d第一页| 亚洲伊人久久精品综合| 亚洲丝袜综合中文字幕| 9色porny在线观看| 波野结衣二区三区在线| 曰老女人黄片| 国产精品一国产av| 久久精品国产亚洲av天美| 各种免费的搞黄视频| 精品人妻在线不人妻| 一级二级三级毛片免费看| 国产伦理片在线播放av一区| 久久 成人 亚洲| 高清午夜精品一区二区三区| 熟女人妻精品中文字幕| 哪个播放器可以免费观看大片| 内地一区二区视频在线| 精品久久蜜臀av无| 熟女av电影| 免费观看性生交大片5| 亚洲av二区三区四区| 亚洲av免费高清在线观看| 三级国产精品欧美在线观看| 人妻系列 视频| 久久精品国产自在天天线| 亚洲国产精品国产精品| 一级a做视频免费观看| 亚洲性久久影院| 91久久精品国产一区二区成人| 日本免费在线观看一区| 麻豆乱淫一区二区| 草草在线视频免费看| 免费少妇av软件| 中国国产av一级| 少妇高潮的动态图| 国产精品嫩草影院av在线观看| 18在线观看网站| 18禁观看日本| 国产亚洲欧美精品永久| 人人妻人人爽人人添夜夜欢视频| 大片免费播放器 马上看| 乱人伦中国视频| 十分钟在线观看高清视频www| 国产淫语在线视频| 午夜免费鲁丝| 欧美日韩一区二区视频在线观看视频在线| 亚洲人与动物交配视频| 久久精品国产自在天天线| 亚洲欧美色中文字幕在线| 成人影院久久| 国产午夜精品一二区理论片| 一级毛片电影观看| 成人午夜精彩视频在线观看| 亚洲av成人精品一二三区| 亚洲人成网站在线播| 国产精品成人在线| 亚洲五月色婷婷综合| 国产精品偷伦视频观看了| 国产精品蜜桃在线观看| 欧美日韩视频精品一区| 狂野欧美激情性bbbbbb| 在线看a的网站| 免费日韩欧美在线观看| 免费高清在线观看日韩| 国产成人免费无遮挡视频| 精品少妇黑人巨大在线播放| 精品久久国产蜜桃| 欧美成人精品欧美一级黄| 亚洲一区二区三区欧美精品| 大香蕉97超碰在线| 黄色一级大片看看| 91久久精品国产一区二区三区| 精品一品国产午夜福利视频| 一本色道久久久久久精品综合| 蜜桃国产av成人99| 国产精品国产三级专区第一集| 狂野欧美激情性bbbbbb| 中文天堂在线官网| 女性被躁到高潮视频| 成年美女黄网站色视频大全免费 | 国产白丝娇喘喷水9色精品| 一区二区三区乱码不卡18| 伊人亚洲综合成人网| 亚洲av不卡在线观看| 国产男人的电影天堂91| 久久鲁丝午夜福利片| 国产精品嫩草影院av在线观看| 国产精品一国产av| 免费少妇av软件| 黄色怎么调成土黄色| 欧美xxⅹ黑人| 亚洲国产欧美日韩在线播放| 在线观看三级黄色| 人妻制服诱惑在线中文字幕| 久久久久网色| 观看av在线不卡| 中文字幕久久专区| 亚洲色图综合在线观看| 欧美精品一区二区大全| 亚洲人成网站在线播| 亚洲欧洲日产国产| 女性生殖器流出的白浆| 51国产日韩欧美| 日韩不卡一区二区三区视频在线| 午夜视频国产福利| 午夜91福利影院| 欧美亚洲 丝袜 人妻 在线| 国产日韩欧美亚洲二区| 9色porny在线观看| 亚洲欧美精品自产自拍| 女的被弄到高潮叫床怎么办| 午夜久久久在线观看| 国产白丝娇喘喷水9色精品| 亚洲人与动物交配视频| 精品一区二区免费观看| 欧美激情 高清一区二区三区| 午夜福利影视在线免费观看| 国模一区二区三区四区视频| 国产精品 国内视频| av线在线观看网站| 91久久精品国产一区二区三区| 国产高清不卡午夜福利| 国产精品一二三区在线看| 97超视频在线观看视频| 欧美日本中文国产一区发布| 黄色视频在线播放观看不卡| 国产伦精品一区二区三区视频9| 久久久欧美国产精品| 插阴视频在线观看视频| 午夜av观看不卡| 美女中出高潮动态图| 国产成人免费无遮挡视频| 男女边摸边吃奶| 极品少妇高潮喷水抽搐| 国产成人av激情在线播放 | av国产久精品久网站免费入址| 日产精品乱码卡一卡2卡三| 丝袜喷水一区| 日本爱情动作片www.在线观看| 天天操日日干夜夜撸| 久久精品国产亚洲av天美| 亚洲内射少妇av| 人人妻人人澡人人爽人人夜夜| 国产精品国产三级国产专区5o| 中国三级夫妇交换| 人妻制服诱惑在线中文字幕| 亚洲欧美一区二区三区国产| 桃花免费在线播放| 高清视频免费观看一区二区| 亚洲国产色片| 少妇人妻久久综合中文| 麻豆成人av视频| 亚洲欧美中文字幕日韩二区| 亚洲欧美成人综合另类久久久| 日韩在线高清观看一区二区三区| 欧美 亚洲 国产 日韩一| 看非洲黑人一级黄片| 热99久久久久精品小说推荐| 中文字幕精品免费在线观看视频 | 91精品国产国语对白视频| 这个男人来自地球电影免费观看 | 在线观看美女被高潮喷水网站| 亚洲国产精品一区二区三区在线| 蜜桃久久精品国产亚洲av| 亚洲精品国产色婷婷电影| 狠狠婷婷综合久久久久久88av| 如何舔出高潮| 激情五月婷婷亚洲| 精品熟女少妇av免费看| 亚洲av不卡在线观看| 狂野欧美激情性xxxx在线观看| 制服丝袜香蕉在线| 亚洲人成网站在线播| 国产日韩一区二区三区精品不卡 | 精品久久久精品久久久| 亚洲美女黄色视频免费看| 国产淫语在线视频| 日韩欧美精品免费久久| 日本午夜av视频| 18+在线观看网站| av不卡在线播放| 少妇人妻 视频|