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

    Investigation on Liquid Breakup Regimes and Atomization Performance of Slinger

    2022-05-20 06:52:26WenlongTan

    Wenlong Tan

    (Chinese Flight Test Establishment, Xi’an 710089, China)

    Abstract: A theoretical and experimental investigation on the liquid breakup regimes of slinger is presented in this paper. The theoretical results show that in the round orifice of the slinger, if the contact angle between the liquid and the wall is less than 90°, the liquid surface tension will inhibit the liquid lateral flowing under the Coriolis force. The formula used to predict the liquid regime leaving the orifice in the past was proved to be inaccurate, which is explained by two assumptions presented in this paper. The breakup regimes of the liquid film are related to the Weber number, and the characteristic length of the Weber number is equivalent to the thickness of the liquid film. The liquid film performed subcritical breakup when the Weber number was small and performed supercritical breakup when the Weber number was large. But when the Weber number was in the middle state, the breakup regimes could not be determined. The liquid breakup regimes affected the primary breakup of the liquid, but a lot of tiny droplets were produced in the secondary breakup, which led to poor correlation between the spray Sauter Mean Diameter (SMD) and the Weber number corresponding to various liquid breakup regimes. On the whole, the spray SMD performs a power relationship with the Weber number corresponding to the liquid film regime, and the correlation coefficient is close to 0.9.

    Keywords: rotary atomization; slinger; liquid breakup regime; spray SMD

    0 Introduction

    Slinger is a mechanical atomizing device and mainly used in the combustion chamber of small and microturbine engines in virtue of simple structure and lightweight. The liquid leaves the slinger orifice in various regimes and the liquid regimes will affect the atomization performance of the slinger.

    Dahm et al.[1-3]investigated the liquid breakup regimes of the slinger. The liquid leaves the orifice of the slinger in two modes: film mode and stream mode. If the liquid leaves the slinger along the circumference of the orifice and forms a film, it is termed film mode; if the liquid accumulates on the side of the orifice opposite to the slinger rotation direction when leaving, it is termed stream mode. Likewise, the breakup regime of the film mode is divided into two kinds: subcritical breakup and supercritical breakup. Subcritical breakup means the liquid leaves the orifice in the film regime and then aggregates into a liquid column under liquid surface tension. Supercritical breakup means the liquid film breaks into fragments directly by aerodynamic force when leaving the orifice. Ro is the ratio of centrifugal force to Coriolis force and is used to estimate the liquid regime when it leaves the orifice. The larger the Ro, the more likely it is to form stream mode. The breakup regime of film mode is correlated to the Weber number corresponding to the liquid film. The larger the Weber number, the more likely the liquid film performs supercritical breakup. The spray Sauter Mean Diameter (SMD) is correlated to the Weber number with a characteristic length determined by the liquid breakup regimes. The characteristic length of the film-subcritical regime is the diameter of the liquid column aggregated by the liquid film, while that of film-supercritical regime is the thickness of the film and that of the stream mode is the orifice diameter.

    Besides, most of the existing research is focused on the liquid spray characteristics of the slinger.Morishita[4]studied the atomization characteristics of three slingers with different structures. He found the rotary speed of the slinger is the most important factor that affects the atomization, and presented the relationship formula between spray SMD and rotary speed. Choi et al.[5-9]investigated the spray spatial distribution and how the crossflow affected the atomization when the slinger was in a crossflow whose velocity direction is parallel to the axis of the slinger. Sescu et al.[10-13]studied the liquid regimes in the orifice of slinger by numerical simulation, and used Phase Doppler Particle Analyzer (PDPA) to measure the spray SMD at different rotary speeds. Jones[14]designed a new type of slinger that uses a swirl nozzle to replace the traditional orifice and found the atomization performance became worse. Kottala[15]performed a numerical simulation on the slinger, which is mainly focused on the primary and secondary breakup of the liquid near the orifice. Ye et al.[16]studied the liquid regimes of in and out the slinger by three-dimensional numerical simulation, and found the total liquid flowrate of the slinger varies significantly with the time and the liquid flowrates in each orifice are different in a short time. Chakraborty and Sahu[17]studied the liquid regimes in the orifice of slinger when the rotary speed is high. The results show that the Coriolis force cannot be ignored in most cases and the liquid will accumulate on the side opposite to the rotary direction when leaving the orifice. Rezayat et al.[18]studied the liquid breakup regimes of the slinger, and found when the rotary speed is up to 4000 r/min, the liquid performs subcritical breakup; when the rotary speed is more than 5000 r/min, a phenomenon of liquid jet bifurcation is observed; when the rotary speed is up to 10000 r/min, the liquid performs supercritical breakup.

    Among all the research about the slinger atomization of past scholars, Werner J. A. Dahm’s research is very valuable and has been cited in many papers. Though Dahm et al.[1-3]carried out a detailed and deep study on liquid breakup regimes and atomization of a slinger, the effect of the liquid surface tension on the liquid regimes is not involved, and the conclusions presented in the papers have not been verified by reliable experimental data. This paper will perform theoretical analysis and experimental research on the slinger to verify, correct, and extend Dahm’s conclusions.

    1 Effect of Structure Tension

    If the shape of the slinger orifice is not round, such as square, slot, etc., the liquid surface tension will drive the liquid to accumulate on the corners of the orifice for their small curvature radius[1-3]. However, if the orifice is round whose curvature radius is the same everywhere, how the liquid surface tension affects the liquid regime needs to be studied.

    1.1 Liquid Regime in the Orifice (Excluding Coriolis Force)

    In Fig.1,R1denotes the inner diameter of the slinger andR0denotes the outer diameter. The orifice of the slinger is round, andd,ldenote its diameter and length respectively, wherel=R0-R1.ωdenotes the rotary speed of the slinger. According to the research of Dahm et al.[3], the liquid will form a film in a very short distance after entering the orifice, as shown in Fig.2, in whichx,ydenote the coordinate axes;Rdenotes the distance from the calculation point to the slinger axis;h(R) is the thickness of the liquid film that has a distance ofRaway from the axis, which in most cases is very small; a denotes the centrifugal acceleration, wherea=ω2R.

    Fig.1 Schematic diagram of the slinger[3]

    Fig.2 Schematic diagram of the liquid on the inner wall of the orifice (axis in the picture denotes the slinger axis)[3]

    The Reynolds number of the liquid filmis Re=(ρLuaveh)/μL, whereρLdenotes liquid density,uavedenotes the average velocity of the liquid in the film in thex-axis direction, andμLdenotes liquid dynamic viscosity. Forhis very small and Re is usually less than 1000, the liquid can be regarded as a laminar flow. The velocity of the liquid iny-axis direction can be neglected when the liquid film becomes stabilized. Dahm et al.[3]found

    (1)

    whereudenotes the velocity of the liquid in thex-direction,ωdenotes the rotary speed of the slinger.

    The overall liquid massflowrate of the slinger isM, the number of the orifices on the slinger isN, and the mass flowrate corresponding to a single orifice ism=M/N.The thickness of the film is calculated as

    (2)

    TheRof the orifice can be regarded as constant forR0?lusually, andR=(R0+R1)/2, thush(R) is constant in radial direction as well.

    Substituting Eq.(1) into Eq.(2), the following equation is obtained:

    (3)

    The average velocity of the liquid in thex-direction is calculated by Eq.(3):

    (4)

    Assumecis the ratio of film thickness to orifice diameter, then

    (5)

    Based on the practical application of the slinger in turbine engines[19], the value of the parameters in the previous formula is shown as follows:

    In Table 1,σdenotes liquid surface tension,ρadenotes air density. The maximum rotary speed of the slingerωmaxis inversely proportional to the slinger radiusR, whereR=5 cm,ωmax=40000 r/min; orR=1 cm,ωmax=100000 r/min.

    Table 1 Parameters value

    The fuel is kerosene when the slinger practically works, but in the experiment, it is dangerous to atomize kerosene, thus water is usually used instead. The parameters of the two liquids are shown in Table 2.

    Table 2 Liquids parameters

    With reference to Table 1, Table 2, and the actual working conditions of the slinger, the thickness of the liquid film in the orifice is usually much smaller than the diameter of the orifice, and the range ofccan be estimated as follows:

    0.003≤c≤0.1

    (6)

    1.2 Liquid Regime in the Orifice (Including Coriolis Force)

    According to the research of Dahm et al.[1-3], after entering the orifice, the liquid will form a regime as shown in Fig.3 in a short radial distance. Then the liquid will flow to the side of the orifice reverse to the slinger rotary direction due to Coriolis force and finally form a liquid regime as shown in Fig.4. When the liquid regime changes from Fig.3 to Fig.4, the change of the surface energy isWσas the surface area of the liquid changes; the change of the wall adhesion work isWsvas the area of the wall covered by the liquid changes.

    Fig.3 Schematic diagram of liquid lateral movement in the orifice (ω denotes rotation direction)

    The liquid area in Fig.4 is a sector with a central angleθ, and is equivalent to the liquid area in Fig.3.

    θ-sinθ=2π(2c-c2)

    (7)

    Fig.4 Schematic diagram of liquid distribution in steady state under Coriolis force (ω denotes rotation direction)

    According to Eq.(6), the high-order terms in Eq.(7) can be ignored:

    θ-sinθ≈4πc

    (8)

    Since Eq.(7) cannot be solved directly, the equation is firstly solved numerically, and then an approximate alternative solution of the equation is found by curve fitting. In the range of Eq.(6), 1000 points is equally divided, and Eq.(7) is numerically solved at each point.

    Fig.5 is the fitting curve of Eq.(7), and the fitting formula is

    θ=4.608c0.35

    (9)

    The correlation coefficient is 0.9998, thus Eq.(9) can be used to replace Eq.(7).

    According to Eq.(6) and (7),

    0.613≤θ≤2.071

    (10)

    Fig.5 Fitting curve of Eq.(7)

    1.3 Influence Law of Surface Tension

    From Fig.3 to Fig.4, the change of surface energy is

    (11)

    As shown in Fig.6, the contact angle of liquid and wall isα, where 0<α<π.When the wall covered by liquid changes to be covered by gas, it is necessary to overcome the wall adhesion. Per unit area wall covered by liquid changes to be covered by saturated vapor, and the adhesion work is[20]

    Wsv'=σ(1+cosα)

    (12)

    Fig.6 Diagram of the contact angle between liquid and wall

    As the air near the wall in the orifice can be considered as saturated vapor of the liquid,with reference to Fig.3-4 and Eq.(12), the following equation is obtained:

    Wsv=σr(1+cosα)(2π-θ)

    (13)

    whereWσis positive work,Wsvis negative work, and the total work of the surface tensionσis

    (14)

    Differentiating Eq.(14) byσ, the following equation is obtained:

    (15)

    Perform Taylor series expansionon sinθ, then

    (16)

    In Eq.(16), the larger theθ, the larger theo(θ5).Combine them with Eq.(10), then

    (17)

    Therefore, the 5th-order Taylor series can be used to replace sinθ.Substituting Eq.(16) into Eq.(8), the following equation is obtained:

    (18)

    Substituting sin(θ/2) into Eq.(15) with a 5th-order Taylor series and referring to Eq.(18), Eq.(19) is obtained:

    (19)

    By reference to Eq.(6),(8),(9), and (10) and numerical calculation, the following inequality is obtained:

    (20)

    Whenα≤90°, which is equivalent to cosα≥0, it is obvious ?W/?σ<0.In this case, as the surface tension increases, the liquid in the orifice becomes more difficult to accumulate on one side of the orifice.

    Whenα>90°, which is equivalent to cosα<0, if ?W/?σ>0,

    (21)

    cosαwill reach the maximum ascreaches its maximum. Combine it with Eq.(6) and (10), the following is obtained:

    cosα<-0.06

    (22)

    Whenα>93.5°, ?W/?σ>0.In this case, as the surface tension increases, it becomes easier for the liquid in the orifice to accumulate on one side of the orifice.

    In summary, there is a critical angleα0, where

    (23)

    Ifα<α0, the surface tension of the liquid inhibits the liquid from accumulating on one side of the orifice, and the smaller theα, the greater the effect of surface tension. Ifα=α0, liquid surface tension exerts no effect on liquid aggregation. Ifα>α0, the surface tension will promote the liquid to accumulate on one side of the orifice, and the larger theα, the greater the effect of surface tension. With reference to Eqs.(2),(9), and (23), it is found that the value ofα0is related toc, and the larger thec, the larger theα0, which means increases in liquid flowrate or decreases in diameter of the orifice or rotary speed will extend the range of the contact angle in which the surface tension will inhibit the liquid lateral movement in the orifice.

    2 Experimental Investigation

    2.1 Test Apparatus

    As shown in Fig.7, the test slinger is driven by a motor mounted on a bracket and controlled by a frequency regulator. The rotary speed of the motor was in the range of 0-24000 r/min, and the minimum adjustment unit was 6 r/min. A transparent shield surrounded the test slinger to collect and return the liquid into the liquid tank. The liquid in the tank was driven by a pump, and sequentially passed through filter, pump, valve, mass flowmeter, and inlet pipe, then hit on the wall of the test slinger, which was thrown out and started a new cycle. A fine-turning valve was used to adjust the liquid flowrate precisely. The measuring range of the mass flowmeter was 0-40 g/s and the minimum measurement unit corresponded to 0.1 g/s. A high-speed camera was used to capture the liquid regime and the maximum capture frequency was 2500 fps. The spray SMD was measured by a PDPA. The measurement error of PDPA was related to the types of measured particles when liquid spray was being measured, and the relative error of particle size measurement was 0.5%[21], which is too little to be shown in the chart of this paper.

    Fig.7 Diagram of test apparatus

    The structure of the test slinger is shown in Fig.1. There was only one orifice on the slinger to prevent mutual interference between the liquids from plurality orifices. The parameters of the test slinger are shown in Table 3.

    Table 3 Parameter of test slinger( mm)

    2.2 Liquid Regime in Orifice

    Fig.8 is the pictures of the liquid regimes in the orifice. Besides film mode (II) and stream mode (III) found by Dahm et al.[3], the liquid may fill the orifice if the rotary speed is low or the liquid flowrate is large, which is termed full mode (I).

    Fig.8 Liquid regimes in the orifice (rotary direction is anticlockwise)

    Fig.9 is the breakup mode of liquid film. As shown in Fig.8 and Fig.9, different liquid regimes can be distinguished clearly except a few middle states. In the experiment, multiple pictures were taken for the liquid regime at every state in order to eliminate manual errors.

    Film-subcritical Film-supercritical

    Rois the ratio of centrifugal force to Coriolis force on the liquid in the orifice:

    Ro=Rω/uave

    (24)

    δis the liquid lateral displacement driven by Coriolis force when the liquid is entering or leaving the orifice. Dahm believeδ/dcan be used to estimate whether the liquid regime is film mode or stream mode.

    δ/d≈0.5(l/d)(l/R)Ro

    (25)

    The lager theδ/d, the more likely the liquid performs stream mode. Conversely, the smaller theδ/d, the more likely it performs film mode[3].

    Fig.10 is the relationship between liquid regimes andδ/d.Theoretically,δ≤d, ifδ/d>1, the liquid regime must be stream mode. Whenδ/d<0.24, the three liquid regimes all appear. In full mode, the orifice is filled with the liquid, indicating it is reasonable thatδ/dis small. But neither the appearance of the stream mode whenδ/dis small nor the existence of the film mode whenδ/dis close to 1 can be explained by the theory of Dahm et al.[3].

    Fig.10 Relationship between liquid regimes and δ/d (On the ordinate, 1 denotes full mode, 2 denotes film mode, 3 denotes stream mode)

    The effects of liquid viscosity and surface tension were not considered in Eq.(25). During the transition from Fig.3 to Fig.4, the liquid viscosity prevents the liquid from flowing to one side of the orifice to form stream mode. The liquid surface tension will inhibit the liquid lateral flow if the contact angleα<90°.In most cases, the contact angle between the liquid and the slinger is smaller than 90°, therefore under the combined effect of liquid viscosity and surface tension, the liquid can keep the film mode even whenδ/dis close to 1.

    Fig.11 is the relationship betweenmandδ/dof stream mode and film mode. When the liquid performs stream mode, the larger them, the smaller theδ/dwhen the stream mode first appears. Once the stream mode appears, increases in rotary speed will not change the liquid regime, and ifm>12 g/s, increases inδ/dwill not change the liquid regime either. In the same orifice,δ/dcorresponding to stream mode is larger than that of film mode, and once the liquid regime changes from film mode to stream mode, the liquid regime will keep stream mode with increases inδ/d.Whenm<6 g/s, film mode and stream mode appear alternatively with increases inδ/d, which means theδ/dcorresponding to film mode and stream mode are different in different orifices, andδ/dcannot be used to estimate the liquid regime in various orifices.

    The film mode appears onlywhenm≤12 g/s. Ifm>12 g/s, film mode disappears and only full mode or stream mode occurs. Whenδ/dis small, the liquid regime performs full mode; asδ/dincreases, the liquid regime switches from full mode to stream mode directly without going through film mode, even whenδ/dis very small.

    Besides, it is found the liquid regime will change from film mode to stream mode suddenly under certain conditions, as shown in Fig.12. The rotary speed only increases by 60 r/min, and the liquid regime switches from film mode to stream mode without any intermediate state.

    Fig.11 Liquid regimes verses m and δ/d

    For the phenomena above cannot be explained by the theories of Dahm et al.[2-3], two hypotheses are presented to provide a reasonable explanation.

    2.2.1DrainPoolHypothesis

    When the liquid flowrate entering the orifice is greater than the maximum flowrate at which the orifice can pass through, a liquid film will be formed in the slinger and cover the entrance of the orifice, as shown in picture (a) of Fig.13. The process of the liquid passing through the orifice can be regarded as a drained pool.

    (a) 3300 r/min

    (b) 3360 r/min

    The liquid in the orifice will be accelerated by centrifugal force. According to the conservation of mass, an air recirculation zone may form at the exit of the orifice.Besides centrifugal force, the liquid in the orifice is also affected by the Coriolis force. Fig.14 is the liquid radial speed distribution at the exit of the orifice calculated by numerical simulation. A cylinder zone with a diameter of 1 mm and a height of 3 mm is created to simulate the fluid zone of the orifice, the boundary condition of the inlet is mass flow inlet and the mass flowrate is 4 g/s. The boundary condition of outlet is pressure outlet and the relative pressure is 0. The fluid zone rotates around the axis of the slinger, and the rotary speed equals the rotary speed of slinger. The turbulent mode isk-ε, the mesh number is 0.1 million, and the calculation state is steady.

    According to the results, it is found that the liquid part with low velocity is concentrated on the same side as the rotation direction, where the air is most likely to return and form an air recirculation zone.At the same time, a vortex appears on the liquid film surface inner the slinger, as shown in picture (b) of Fig.13. As the rotary speed increases, the vortex and the air recirculation zone become larger and larger, and are eventually connected, as shown in picture (c)-(d) of Fig.13.

    The liquid regimes corresponding to each figure of Fig.13 is as follows: full mode corresponds to picture (a), stream mode corresponds to others, and as the rotary speed increases, the liquid accumulation of stream mode is increasing. The above information provides a reasonable explanation for the phenomenon that the liquid regime changes from full mode to stream mode directly, and even whenδ/dis very small, stream mode can also be formed.

    (a) Full covered by liquid (b) Air recirculation zone occurs

    (c) Joint of in and out air zone (d) Full development of air zone

    (a) 1500 r/min (b) 3000 r/min (c) 5000 r/min

    2.2.2TangentialRelativeSlipHypothesis

    The rotary speed of the liquid in the cavity of the slinger is not exactly equivalent to the slinger rotary speed. There is a tangential relative slip between the liquid and the wall of the slinger. Assumekis the tangential slip ratio, where

    k=(Vθ-vθ)/Vθ=(ωR-vθ)/ωR

    (26)

    whereVθdenotes the tangential velocity of the slinger andvθdenotes the tangential velocity of the liquid. The research of Ahmed suggests that the tangential relative slip ratio between the liquid and the rotary atomizer is about 18%-40%[22].

    When the tangential relative slip velocity between the liquid and the wall in the slinger is low, the liquid flows closely on the wall of the orifice after entering the orifice under the combined effect of Coanda effect, surface tension, viscous and centrifugal force, and forms the liquid regime as shown in Fig.3 within a short radial distance. However, when the tangential relative slip velocity reaches a certain level, the liquid will enter the orifice in the form of a parabola, as shown in Fig.15. Then the liquid hits on the side of the orifice opposite to the slinger rotation direction, therefore the liquid regime outside the orifice will change from film mode to stream mode suddenly.

    Fig.15 Schematic diagram of liquid regime entering the orifice (the thick line indicates the liquid, the thin line indicates the wall, ω denotes rotation direction)

    2.3 Subcritical and Supercritical Breakup

    Dahm et al.[1-3]considered the breakup regimes of film mode are related to Weber number corresponding to the film, and the decision parameter isT, where

    (27)

    whereucdenotes the tangential velocity at the edge of the slinger. The smaller theT, the more likely the film performs subcritical breakup; the larger theT, the more likely the liquid performs supercritical breakup.

    2.4 Spray SMD

    Ohtdenotes Ohnesorge number, andWetdenotes Weber number, where

    (28)

    (29)

    Fig.16 Relationship between breakup regime and Weber number (On the ordinate, 1 denotes subcritical, 2 denotes supercritical)

    Dahm et al. analyzed the experimental data from the perspective of Morishita’s research[3-4]. The results show that the correlation coefficient of SMD/d~Wedor SMD/d’~Wed’is 0.81, and one of SMD/h~Wehis only 0.46. This indicates the liquid performs film-subcritical breakup or stream mode in Morishita’s experiment, and film-supercritical breakup hardly appears. But the theoretical results of the research of Dahm et al.[1-3]indicate the liquid in Morishita’s experiment mainly performs film-supercritical breakup. He considered the data of spray SMD was inaccurate because the samples of the spray were collected on a silicon oil film on a glass plate and then were photographed to measure the size of the droplets[4]. The method cannot get the spray SMD accurately. However, according to the research in Section 2.2 of this paper, the liquid in Morishita’s experiment performed stream mode in most cases, and the conclusion of Dahm’s research on the liquid regimes about Morishita’s experiment is right.

    Fig.17 is the correlation between spray SMD and Weber numbers corresponding to various liquid breakup regimes based on the experimental data in this paper. The correlation coefficient of film-subcritical is 0.081, film-supercritical is 0.434, stream mode is 0.323, all of which show a bad correspondence between spray SMD and Weber number. The liquid breakup regime only affects the primary breakup of the liquid, and in the process of the secondary breakup, a lot of tiny droplets are produced when the large droplets are split into small droplets by aerodynamics force, and little correlation exists between these tiny droplets and the liquid breakup regimes. These tiny droplets make the correlation between the spray SMD and liquid breakup regime worse, therefore the correlation coefficients are poor in Fig.17.

    (a)Film with subcritical breakup

    (b)Film with supercritical breakup

    (c)Stream mode

    The measurement method of Morishita ignored the tiny droplets[4]. The spray SMD measured by PDPA in this paper is smaller than Morishita’s under the same operating condition. The size of bigger droplets in the spray is affected by the primary breakup, that is affected by liquid breakup regimes. Morishita summed up the spray SMD by measuring the size of bigger droplets, therefore the spray SMD coincidentally shows good correlation to the Weber number corresponding to stream mode.

    Though the liquid breakup regime only affects the primary breakup, the spray SMD is affected by Weber number. Fig.18 shows the fitting curve between spray SMD andWeh, where

    Fig.18 Correlation between spray SMD and Weh of slinger

    The correlation coefficient is 0.897. Compared withWedandWed’, the correlation coefficient of Eq.(30) is the largest, which indicates the characteristic length of the liquid corresponds to the thickness of the film calculated by Eq.(2) in most cases.

    The liquid breakup regime in the experiment of this paper is mainly stream mode. Fig.19 shows two common breakup regimes of stream mode. In Fig.19(a), although the liquid has accumulated on one side of the orifice, it re-opens into a liquid film under the action of air after leaving the orifice and then performs primary breakup. In Fig.19(b), the tangential speed of the liquid is very high, and the liquid stream will be torn into lots of fragments violently by aerodynamic force after leaving the orifice. The characteristic lengths of both breakup regimes are close toh, therefore the correlation between spray SMD andWehis the best.

    (a) Re-open to film (b) Violently breakup

    3 Conclusions

    A theoretical analysis was performed on slinger, which shows the liquid surface tension and the contact angleαbetween liquid and slinger wall play an important role in liquid regimes. If the orifice is round whose curvature radius is the same everywhere, there is a critical angleα0>90°; ifα<α0, the liquid surface tension will inhibit the liquid from flowing to one side of the orifice driven by Coriolis force; ifα>α0, the liquid surface tension will promote the liquid to flow to the side.

    Three liquid regimes appeared in the orifice: full mode, film mode, and stream mode, andδ/dcould not be used to estimate the liquid regime accurately. The reason may be that when liquid flowrate is large enough, the process of the liquid passing through the orifice is just like a drained pool, the air outside the orifice returns from the same side of the orifice as the rotation direction, and the liquid regime changes from full mode to stream mode directly. If the flowrate is low, a tangential relative slip exists between the liquid and the wall in the slinger. When the slip velocity is large enough, the liquid may enter the orifice in a parabola form, and the liquid regime outside the orifice changes from film mode to stream mode within a little increase in rotary speed.

    The liquid breakup regimes could affect the primary breakup of the liquid, but a lot of tiny droplets were produced in the secondary breakup with little correlation to the liquid breakup regime. Therefore, the correlation between SMD/tandWetwas poor.

    On the whole, the spray SMD has a power relationship withWeh, and the correlation coefficient is close to 0.9.

    女的被弄到高潮叫床怎么办| 丰满乱子伦码专区| 久久人妻av系列| 久久九九热精品免费| 黄色视频,在线免费观看| 久久热精品热| 精品少妇黑人巨大在线播放 | 一级黄色大片毛片| 97超碰精品成人国产| 蜜桃亚洲精品一区二区三区| 久久6这里有精品| 日韩一本色道免费dvd| 99热精品在线国产| 亚洲欧美中文字幕日韩二区| 91久久精品国产一区二区三区| 国产精品乱码一区二三区的特点| 亚洲天堂国产精品一区在线| 免费观看的影片在线观看| or卡值多少钱| 在线a可以看的网站| 国内精品久久久久精免费| 美女xxoo啪啪120秒动态图| 美女免费视频网站| 网址你懂的国产日韩在线| 老师上课跳d突然被开到最大视频| 国产欧美日韩一区二区精品| 一级毛片久久久久久久久女| 久久久精品94久久精品| 亚洲18禁久久av| 久久久久久久久久黄片| 久久久久久久久大av| 国产蜜桃级精品一区二区三区| 日本五十路高清| 久久精品国产亚洲av涩爱 | 国产黄色视频一区二区在线观看 | 插逼视频在线观看| 国产综合懂色| 亚洲国产日韩欧美精品在线观看| 国产探花极品一区二区| 看黄色毛片网站| 美女高潮的动态| 性插视频无遮挡在线免费观看| 五月玫瑰六月丁香| 2021天堂中文幕一二区在线观| 变态另类成人亚洲欧美熟女| 精品久久久久久久久久久久久| 内地一区二区视频在线| 插阴视频在线观看视频| 丰满乱子伦码专区| 国产伦精品一区二区三区视频9| 国产精品综合久久久久久久免费| 看十八女毛片水多多多| 国产亚洲精品综合一区在线观看| 欧洲精品卡2卡3卡4卡5卡区| 久久精品夜色国产| 最后的刺客免费高清国语| 一级毛片电影观看 | 91久久精品国产一区二区三区| 菩萨蛮人人尽说江南好唐韦庄 | 国产成年人精品一区二区| 欧美日韩在线观看h| av在线老鸭窝| 亚洲精品影视一区二区三区av| 国产国拍精品亚洲av在线观看| 蜜臀久久99精品久久宅男| 特级一级黄色大片| 日韩大尺度精品在线看网址| 免费黄网站久久成人精品| 日本黄色片子视频| 直男gayav资源| 能在线免费观看的黄片| 床上黄色一级片| 久久久久久久午夜电影| 日韩制服骚丝袜av| 深夜a级毛片| 麻豆久久精品国产亚洲av| 女同久久另类99精品国产91| 日韩,欧美,国产一区二区三区 | 熟女电影av网| 寂寞人妻少妇视频99o| 亚洲最大成人中文| 国产精品99久久久久久久久| 伦精品一区二区三区| 亚洲最大成人中文| 又爽又黄a免费视频| 久久精品国产亚洲av天美| 在线免费观看不下载黄p国产| 日日摸夜夜添夜夜爱| 国产成人福利小说| 色综合色国产| 熟妇人妻久久中文字幕3abv| 亚洲久久久久久中文字幕| 成年av动漫网址| 少妇裸体淫交视频免费看高清| 亚洲性夜色夜夜综合| 欧美一区二区亚洲| 色吧在线观看| 三级国产精品欧美在线观看| 伦理电影大哥的女人| 亚洲成人av在线免费| a级一级毛片免费在线观看| 麻豆一二三区av精品| 99久久成人亚洲精品观看| 午夜久久久久精精品| 亚洲一区高清亚洲精品| 晚上一个人看的免费电影| 久99久视频精品免费| 校园人妻丝袜中文字幕| 草草在线视频免费看| 寂寞人妻少妇视频99o| 又黄又爽又免费观看的视频| 欧美极品一区二区三区四区| 国产成人aa在线观看| 非洲黑人性xxxx精品又粗又长| 国产黄色视频一区二区在线观看 | 欧美人与善性xxx| 久久久久久久久久久丰满| 亚洲美女搞黄在线观看 | 久久欧美精品欧美久久欧美| 永久网站在线| 一进一出抽搐gif免费好疼| 男人舔奶头视频| 午夜福利18| 在线观看一区二区三区| 久久99热这里只有精品18| 国产黄a三级三级三级人| 中文字幕熟女人妻在线| 岛国在线免费视频观看| 免费搜索国产男女视频| 午夜精品在线福利| 精品无人区乱码1区二区| 在线播放国产精品三级| 日日摸夜夜添夜夜添小说| 精品久久久久久久久久免费视频| 久久久久久久久中文| 欧洲精品卡2卡3卡4卡5卡区| 老司机影院成人| 精品欧美国产一区二区三| 麻豆成人午夜福利视频| 国产美女午夜福利| 乱码一卡2卡4卡精品| 午夜激情欧美在线| 女人被狂操c到高潮| 日本一二三区视频观看| 亚洲自拍偷在线| 日本五十路高清| 国内精品美女久久久久久| 亚洲成人av在线免费| 熟女电影av网| 亚洲第一区二区三区不卡| 国产乱人偷精品视频| 久久精品国产亚洲网站| 51国产日韩欧美| 亚洲成人久久爱视频| 国产精品久久久久久久久免| 亚洲欧美精品综合久久99| 欧美中文日本在线观看视频| 亚洲成人久久爱视频| 精品一区二区三区视频在线| 国产成人影院久久av| 午夜亚洲福利在线播放| 久久天躁狠狠躁夜夜2o2o| 久久精品国产99精品国产亚洲性色| 免费人成在线观看视频色| 久久精品91蜜桃| 午夜福利18| 日韩欧美一区二区三区在线观看| 精品国产三级普通话版| 床上黄色一级片| 久久人人爽人人爽人人片va| 长腿黑丝高跟| 免费看美女性在线毛片视频| 俺也久久电影网| 老司机影院成人| 午夜激情欧美在线| 欧美极品一区二区三区四区| 成人综合一区亚洲| 中国美白少妇内射xxxbb| 大香蕉久久网| 大香蕉久久网| 国内精品美女久久久久久| 日本五十路高清| 亚洲精品日韩av片在线观看| 亚洲18禁久久av| 97热精品久久久久久| 久久人人爽人人片av| 成年女人毛片免费观看观看9| 欧美成人精品欧美一级黄| 国产69精品久久久久777片| 亚洲图色成人| 免费人成在线观看视频色| 搡老熟女国产l中国老女人| 小说图片视频综合网站| 97碰自拍视频| 久久精品久久久久久噜噜老黄 | 夜夜爽天天搞| 在线观看美女被高潮喷水网站| 亚洲最大成人av| 成人漫画全彩无遮挡| av视频在线观看入口| 日日摸夜夜添夜夜爱| 午夜福利视频1000在线观看| videossex国产| av女优亚洲男人天堂| 久99久视频精品免费| 亚洲精品乱码久久久v下载方式| 美女被艹到高潮喷水动态| 成人毛片a级毛片在线播放| 长腿黑丝高跟| av免费在线看不卡| 午夜亚洲福利在线播放| 亚洲精品乱码久久久v下载方式| 人妻夜夜爽99麻豆av| av国产免费在线观看| 亚洲aⅴ乱码一区二区在线播放| 午夜激情欧美在线| 国产色婷婷99| av天堂在线播放| 一级毛片久久久久久久久女| 少妇熟女欧美另类| 男人的好看免费观看在线视频| 又爽又黄无遮挡网站| 看黄色毛片网站| 嫩草影视91久久| 亚洲最大成人av| 啦啦啦啦在线视频资源| 一本久久中文字幕| 一卡2卡三卡四卡精品乱码亚洲| 97超碰精品成人国产| 亚洲精品亚洲一区二区| 男女下面进入的视频免费午夜| 亚洲精品国产成人久久av| 啦啦啦啦在线视频资源| 欧美+日韩+精品| 日本黄色视频三级网站网址| 亚洲18禁久久av| 国产欧美日韩一区二区精品| 99热精品在线国产| 免费看a级黄色片| 超碰av人人做人人爽久久| 又爽又黄无遮挡网站| 搞女人的毛片| 午夜a级毛片| 欧美xxxx黑人xx丫x性爽| 国内久久婷婷六月综合欲色啪| 全区人妻精品视频| 免费高清视频大片| 精品一区二区三区视频在线| 99久久精品一区二区三区| 又粗又爽又猛毛片免费看| 九九热线精品视视频播放| 在线看三级毛片| 两性午夜刺激爽爽歪歪视频在线观看| 深夜a级毛片| 久久国内精品自在自线图片| 国产探花极品一区二区| 日本精品一区二区三区蜜桃| 男人和女人高潮做爰伦理| 久久精品国产亚洲av天美| 免费黄网站久久成人精品| 波多野结衣巨乳人妻| 久久久久国产网址| 韩国av在线不卡| 国产精品美女特级片免费视频播放器| 日日干狠狠操夜夜爽| 丰满人妻一区二区三区视频av| 一级毛片电影观看 | 一级毛片aaaaaa免费看小| 亚洲精品日韩在线中文字幕 | 亚洲精品影视一区二区三区av| 国产又黄又爽又无遮挡在线| 12—13女人毛片做爰片一| 日韩强制内射视频| 久久精品国产99精品国产亚洲性色| 99久久九九国产精品国产免费| 91麻豆精品激情在线观看国产| 在线免费十八禁| 日日摸夜夜添夜夜添av毛片| 国产三级在线视频| 亚洲五月天丁香| 直男gayav资源| 国产视频一区二区在线看| 免费不卡的大黄色大毛片视频在线观看 | 美女免费视频网站| 久久精品夜色国产| 人妻夜夜爽99麻豆av| 在线观看午夜福利视频| 日韩国内少妇激情av| 一级毛片电影观看 | 午夜精品在线福利| 熟女电影av网| 99riav亚洲国产免费| 人妻制服诱惑在线中文字幕| 偷拍熟女少妇极品色| 成人欧美大片| 少妇的逼好多水| 国产aⅴ精品一区二区三区波| 久久精品久久久久久噜噜老黄 | 欧美另类亚洲清纯唯美| 精品日产1卡2卡| 亚洲国产精品国产精品| 一区福利在线观看| 晚上一个人看的免费电影| 日产精品乱码卡一卡2卡三| 99在线视频只有这里精品首页| 欧美区成人在线视频| 岛国在线免费视频观看| 麻豆一二三区av精品| 十八禁网站免费在线| 久久草成人影院| 亚洲精品亚洲一区二区| 亚洲专区国产一区二区| 精品乱码久久久久久99久播| 韩国av在线不卡| 国产精品国产高清国产av| 国产精品av视频在线免费观看| 偷拍熟女少妇极品色| 国产乱人偷精品视频| 精品人妻偷拍中文字幕| 久久综合国产亚洲精品| 色尼玛亚洲综合影院| av福利片在线观看| 91av网一区二区| 丰满人妻一区二区三区视频av| 久久亚洲国产成人精品v| 中文字幕av在线有码专区| av免费在线看不卡| 老司机福利观看| 亚洲精品亚洲一区二区| 男插女下体视频免费在线播放| 国内精品一区二区在线观看| 亚洲国产欧美人成| 精品久久久久久久久久久久久| 久久国产乱子免费精品| 可以在线观看的亚洲视频| 在线观看免费视频日本深夜| 日本免费a在线| 青春草视频在线免费观看| 91久久精品电影网| 色吧在线观看| 一本久久中文字幕| 中文字幕av在线有码专区| 嫩草影视91久久| 国产精品一区二区三区四区久久| 亚洲aⅴ乱码一区二区在线播放| 色综合色国产| 97碰自拍视频| 露出奶头的视频| 欧洲精品卡2卡3卡4卡5卡区| 一个人看视频在线观看www免费| 五月伊人婷婷丁香| 日韩高清综合在线| 真人做人爱边吃奶动态| www日本黄色视频网| 久久久色成人| 久久久欧美国产精品| 97碰自拍视频| 在线观看美女被高潮喷水网站| 亚洲人成网站在线播| 欧美区成人在线视频| 日本黄大片高清| 少妇猛男粗大的猛烈进出视频 | 一本精品99久久精品77| 成人三级黄色视频| 一进一出抽搐动态| 欧美日韩一区二区视频在线观看视频在线 | 亚洲美女搞黄在线观看 | 97热精品久久久久久| 可以在线观看毛片的网站| 一区二区三区高清视频在线| 午夜福利18| 久久6这里有精品| 亚洲欧美精品综合久久99| 老熟妇乱子伦视频在线观看| 国产成年人精品一区二区| а√天堂www在线а√下载| 最新在线观看一区二区三区| 国产伦精品一区二区三区四那| 欧美三级亚洲精品| 久久这里只有精品中国| 久久久久久久午夜电影| 国内精品宾馆在线| 人妻丰满熟妇av一区二区三区| 高清午夜精品一区二区三区 | 少妇被粗大猛烈的视频| 一级毛片久久久久久久久女| 欧美成人精品欧美一级黄| 国产精品国产高清国产av| 高清日韩中文字幕在线| 成人av在线播放网站| 中文资源天堂在线| 九九久久精品国产亚洲av麻豆| 亚洲人与动物交配视频| 成人永久免费在线观看视频| 精品乱码久久久久久99久播| 99热6这里只有精品| 久久鲁丝午夜福利片| 免费高清视频大片| 欧美xxxx性猛交bbbb| 黄色一级大片看看| 国产国拍精品亚洲av在线观看| 成人鲁丝片一二三区免费| 成年女人看的毛片在线观看| 日韩成人av中文字幕在线观看 | 婷婷色综合大香蕉| 给我免费播放毛片高清在线观看| 淫妇啪啪啪对白视频| 自拍偷自拍亚洲精品老妇| 久久久久九九精品影院| 国产欧美日韩精品一区二区| 久久国产乱子免费精品| 亚洲欧美日韩高清专用| av专区在线播放| 午夜a级毛片| 日韩 亚洲 欧美在线| 成人特级黄色片久久久久久久| 一级毛片久久久久久久久女| 日本一本二区三区精品| 亚洲欧美日韩无卡精品| 亚洲欧美日韩高清专用| 亚洲av二区三区四区| 免费看美女性在线毛片视频| 国产成人a区在线观看| 精品熟女少妇av免费看| 亚洲最大成人av| 别揉我奶头 嗯啊视频| 亚洲精品一区av在线观看| 亚洲婷婷狠狠爱综合网| 国内精品久久久久精免费| 亚洲欧美日韩无卡精品| 日本一二三区视频观看| 国产毛片a区久久久久| a级毛片免费高清观看在线播放| 久久精品国产亚洲av天美| 亚洲精品日韩在线中文字幕 | 九九热线精品视视频播放| 成熟少妇高潮喷水视频| 免费av不卡在线播放| 国产亚洲91精品色在线| 欧美另类亚洲清纯唯美| 3wmmmm亚洲av在线观看| 日韩亚洲欧美综合| 夜夜看夜夜爽夜夜摸| 欧美xxxx黑人xx丫x性爽| 欧美激情久久久久久爽电影| 亚洲av成人av| 欧美高清性xxxxhd video| 少妇人妻一区二区三区视频| 在线免费观看不下载黄p国产| 精品不卡国产一区二区三区| 少妇裸体淫交视频免费看高清| 国产在线精品亚洲第一网站| 亚洲精品成人久久久久久| 午夜爱爱视频在线播放| 禁无遮挡网站| 狂野欧美白嫩少妇大欣赏| 国产精品乱码一区二三区的特点| 久久精品国产清高在天天线| 在线看三级毛片| 在线天堂最新版资源| 亚洲国产欧美人成| 久久草成人影院| 成人av一区二区三区在线看| av在线天堂中文字幕| 国产精品永久免费网站| 嫩草影视91久久| av视频在线观看入口| 成人综合一区亚洲| 2021天堂中文幕一二区在线观| 国产成人a∨麻豆精品| 欧美最黄视频在线播放免费| 亚洲精品在线观看二区| 日韩欧美精品免费久久| 国产在线男女| 搡女人真爽免费视频火全软件 | 校园春色视频在线观看| 久久精品91蜜桃| 欧美zozozo另类| 欧美潮喷喷水| 久久亚洲精品不卡| 国产黄色视频一区二区在线观看 | 男女视频在线观看网站免费| 联通29元200g的流量卡| 亚洲成人中文字幕在线播放| 欧美色视频一区免费| 看免费成人av毛片| 长腿黑丝高跟| 中文字幕免费在线视频6| 婷婷精品国产亚洲av在线| 在线观看午夜福利视频| 尾随美女入室| 国产一区亚洲一区在线观看| 人妻久久中文字幕网| 成人av一区二区三区在线看| 成人午夜高清在线视频| 久久久欧美国产精品| 美女黄网站色视频| 男人狂女人下面高潮的视频| 久久久久久九九精品二区国产| 春色校园在线视频观看| 欧美+日韩+精品| 床上黄色一级片| 亚洲国产精品国产精品| 看非洲黑人一级黄片| 国产黄色视频一区二区在线观看 | 18禁在线无遮挡免费观看视频 | 国产欧美日韩精品亚洲av| 免费黄网站久久成人精品| 欧美+亚洲+日韩+国产| 99riav亚洲国产免费| 精品少妇黑人巨大在线播放 | 亚洲三级黄色毛片| 日本 av在线| 一级黄片播放器| 精品一区二区三区视频在线观看免费| 性插视频无遮挡在线免费观看| 国产精品人妻久久久久久| 99热只有精品国产| 免费看日本二区| 嫩草影视91久久| 日本一二三区视频观看| 大又大粗又爽又黄少妇毛片口| 国产精品久久久久久亚洲av鲁大| 最近手机中文字幕大全| 国产毛片a区久久久久| 乱人视频在线观看| a级一级毛片免费在线观看| 色吧在线观看| 精品一区二区三区av网在线观看| 亚洲国产精品sss在线观看| 最新在线观看一区二区三区| av在线蜜桃| 国产成人影院久久av| 欧美区成人在线视频| 狠狠狠狠99中文字幕| 在线播放无遮挡| 国产一区二区在线av高清观看| 啦啦啦韩国在线观看视频| 亚洲av成人av| 麻豆一二三区av精品| 有码 亚洲区| 久久久精品大字幕| 香蕉av资源在线| 中文字幕免费在线视频6| 人人妻,人人澡人人爽秒播| 麻豆成人午夜福利视频| 日本黄色视频三级网站网址| 亚洲精品影视一区二区三区av| 国产毛片a区久久久久| 搡老妇女老女人老熟妇| 97超碰精品成人国产| 成年女人永久免费观看视频| 午夜精品一区二区三区免费看| www.色视频.com| 老司机福利观看| 亚洲国产精品sss在线观看| 性欧美人与动物交配| 最好的美女福利视频网| 一级黄片播放器| 日韩av不卡免费在线播放| 成人性生交大片免费视频hd| 不卡一级毛片| 国产亚洲av嫩草精品影院| 国产亚洲精品久久久久久毛片| 国产精品久久久久久av不卡| 国产免费男女视频| 男女啪啪激烈高潮av片| 亚洲欧美日韩卡通动漫| 看免费成人av毛片| 免费人成视频x8x8入口观看| 两性午夜刺激爽爽歪歪视频在线观看| 欧美性感艳星| 在线国产一区二区在线| 中国国产av一级| 熟妇人妻久久中文字幕3abv| 国产精品三级大全| 国产伦精品一区二区三区四那| 九色成人免费人妻av| 波多野结衣巨乳人妻| 欧美一区二区国产精品久久精品| 中文字幕av在线有码专区| 露出奶头的视频| 97超碰精品成人国产| 国产精品女同一区二区软件| 国产精品99久久久久久久久| 亚洲七黄色美女视频| 亚洲国产精品合色在线| 久久6这里有精品| 欧美最新免费一区二区三区| 最后的刺客免费高清国语| 日本黄色视频三级网站网址| 国内精品宾馆在线| 国产乱人视频| 99热全是精品| 亚洲乱码一区二区免费版| 国产69精品久久久久777片| 精品人妻视频免费看| 色av中文字幕| 99久久精品热视频| 久久久a久久爽久久v久久| 亚洲国产欧洲综合997久久,| 欧美潮喷喷水| 日韩强制内射视频| 欧美人与善性xxx| 丝袜美腿在线中文| 成人亚洲欧美一区二区av| 日本欧美国产在线视频| 国产女主播在线喷水免费视频网站 | 国内揄拍国产精品人妻在线| 直男gayav资源| 亚洲av免费在线观看| 亚洲三级黄色毛片| 九九热线精品视视频播放| 欧美日韩国产亚洲二区| 亚洲aⅴ乱码一区二区在线播放|