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

    Effect of tongue clearance on hydraulic performance of double support vortex pump

    2019-10-31 07:08:40XingyngLiWnqingChenShengdunZho
    Defence Technology 2019年4期

    Xing-yng Li ,Wn-qing Chen ,Sheng-dun Zho

    a Pump Equipment Engineering Research Center of Shaanxi Province,Xi'an Aeronautical University,Xi'an,710077,China

    b School of Mechanical Engineering,Xi'an Jiaotong University,Xi'an,710049,China

    Keywords:Vortex pump Tongue clearance Non-uniform flow passage Numerical simulation

    A B S T R A C T In order to study influence of tongue clearance on the hydraulic performance of double support vortex pump,three dimensional flow model with symmetric variable size tongue clearance was established by changing section parameter of tongue clearance,and applied to numerically simulate the steady inner flow in a vortex pump based on fixed working conditions and impeller through ANSYS FLUENT.The numerical results show the influence of tongue channel with two sections(non-uniform and uniform sections)on flow field characteristics in vortex pump.Firstly,the variation of tongue channel section changes the vortex structure distribution in circumferential flow channel.Specifically the spiral forward vortex structure in the circumferential direction shifts to large radius side with increasing tongue width.Secondly,the circumferential velocity gradient and axial pressure gradient both diminish with the increase of the tongue section,and the inlet/outlet pressure difference and velocity distribution also reduce with increasing tongue channel section.Finally,for vortex pump with non-uniform section of tongue channel,the head diminishes and the efficiency remains constant approximately with the increase of the tongue width,while the head and the efficiency both diminish for uniform section of tongue channel.

    1. Introduction

    Due to many advantages of simple structure,low power consumption,small size and high dependability,dual-support vortex pumps have been widely used in various occasions such as civil aviation,ships,military armored vehicles,and chemical industries[1],in order to replace other types of pump which face many restrictions in the context of use[2—4].Since the development of the vortex pump theory,due to the constraints of low efficiency in the use of traditional vortex pumps, it has become an important research on how to improve the coupling mechanism between vortex pump rotors and flow channels,in order to improve the output characteristics of vortex pumps under the same rotor conditions,and reduce the internal flow loss for improving working efficiency.In the open literatures,scholars at home and abroad reveal the coupling mechanism of impeller and fluid inside vortex pump, and analyze the energy transfer law between rotating impeller and fluid,based on the study of the flow theoretical model of vortex pump[3,4].At the same time,a large number of scholars[4—7]analyze the flow field structure of different vortex pumps in detail by the numerical simulation method,and obtain the influence rules of design parameters of vortex pump on its performance.Furthermore,the general effects of the structure parameter,such as number of blades and the size of flow channel,on the output performance of vortex pump are analyzed[8—10]by means of the experimental method,which provides an important reference for optimal design of vortex pump. In addition, many domestic research institutions such as Jiangsu university and Lanzhou university of technology have carried out a lot of works on impeller and runner structures [11—16], and Northwest A&F university[17—19]makes remarkable achievement in many aspects such as Modeling oscillation modal interaction in a hydroelectric generating system and Hamiltonian modeling of multi-hydro-turbine governing systems with sharing common penstock.The theoretical analysis and verification conclusions in their published literature provide important references for the design of vortex pumps.

    For the research on flow field mechanism and structure optimization of vortex pump,published literatures show that a lot of works are mainly focused on analyzing mechanism about the influence of impeller structure on internal flow,while few researches have been made on flow field and output performance for the mature application of structure expansion for vortex pump product.Based on the previous research methods and the double support vortex pump with typical closed rectangular symmetric flow passage,this paper studies the effect of two kinds of continuously variable section of tongue flow passage,namely the uniform and non-uniform sections,on the hydraulic performance of the vortex pump.In this paper,a three-dimensional flow field calculation model is designed through Pro-E software,and the steady flow is simulated by ANSYS-FLUENT based on non-structured tetrahedral grid.Based on the structure of the impeller and the inlet/outlet,as well as the design parameters,this paper gives the characteristics of the output performance of the vortex pump about the tongue channel with the uniform and non-uniform cross-sections,and shows the influence of the flow passage section characteristic parameters on the inner flow field in the vortex pump.And then by comparing with the experimental results of the typical dualsupport vortex pump design model,this paper analyzes the influence mechanism of two kinds of sectional size parameters on the head and efficiency of vortex pump,which provides important reference for the optimal design of vortex pump.

    2. Double support vortex pump model

    The vortex pump studied in this paper is a typical closed double support structure,whose structure is shown in Fig.1(a).The vortex pump is composed of a pump body,a pump cover and an impeller.The inlet and outlet of the vortex pump are located at the top of the pump cover.The Fig.1(b)is the internal flow domain of the vortex pump,in which D2is impeller diameter.

    The hydraulic structure about the cross-section of the flow passage of the vortex pump is shown in Fig.2(a),and the parameters in Fig.2(a)are given in Table 1 and Table 2.Fig.2(b)is the axial crosssection diagram of the vortex pump,which shows the location of the inlet/outlet,the distribution pattern of the blades and the location of the tongue.The rotations of the blades,which are distributed averagely along the radial direction,make the fluid move forward along the circumferential direction, when the fluid passes the spinning wheel of the impeller after entering the passage of the vortex pump.Meanwhile,the rotating kinetic energy is transformed into the fluid pressure energy,and finally the fluid discharges from the outlet of the vortex pump.The research in this paper is mainly the effect of the tongue clearance on the performance of the vortex pump,including the characteristic parameter of the tongue channel with variable and uniform cross-sections which are shown in Fig.2(c)and(d),respectively.Table 3 shows the working condition about the cross-sectional size of the tongue channel.In Fig.2(c),through keeping the height of the tongue channel G fixed and changing the width E1~E3,it is studied that the performance and inner flow field of the vortex pump are affected by the tongue clearance with the channel section increasing under the condition of the variable cross-sectional area A of the tongue channel.In Fig.2(d),through keeping the cross-sectional area of the tongue channel fixed and changing the width E1~E3 as well as the height G1~G3,it is studied that the performance and inner flow field of the vortex pump are affected by the tongue clearance under the condition of the fixed cross-sectional area A of the tongue channel.

    3. Numerical method

    3.1. Computational mesh

    Due to the complex structure of the dual support vortex pump fluid domain,the non-structured tetrahedral grid is obtained by ICEM software.The grid distributions of the whole fluid domain and its flow channel section are shown in Fig.3.Local mesh refinement is carried out for the region near blade,interface and channel shape with large gradient.Meanwhile,the domain nearby the wall at the tongue gives a small volume division to prevent the excessive size of the mesh in the variable-section region from affecting the overall grid quality and calculation accuracy.In order to improve the accuracy of the numerical simulation,the grid independence of the model is analyzed by calculating the percentage of the head of the vortex pump based on five different grids.The result is shown in Fig.4.It is indicated clearly that,with the number of grid gradually increasing,the difference in the percentage of the head gradually decreases.Considering the deviation of the calculation results and computation efficiency,the number of computational grid is set as 5.58×106,and the maximum deviation of the head is approximately 2.87%.

    3.2. Governing equations and boundary conditions

    Fig.1.Dual support vortex pump with closed impeller.

    Fig.2.Parameters for the cross-section of the flow channel and the variation character of parameter for the cross-section of the tongue channel of the typical double support vortex pump.

    Table 1 Design parameters of vortex pump.

    Table 2 Structure parameters of vortex pump.

    Based on the three-dimensional steady Reynolds-average N-S equations and RNG k-epsilon turbulence model,the flow in the double support vortex pump with rectangular channel was simulated with the segregated implicit method and the SIMPLEC algorithm of the pressure-velocity coupling. The Pressure term is computed by the two-order central difference scheme,and the velocity term,turbulent kinetic energy term and turbulence viscosity coefficient are computed by the two-order upwind scheme.For the numerical iterative computation,the convergence precision is 10-5[12,13].

    About the numerical simulation for flow in the vortex pump based on the ANSYS-FLUENT,the boundary conditions are as follows:(1)The inlet is the boundary mass flow rate,and the flow rate was 0.6 m3/h according to the design parameters in Tables 1 and 2.(2)The water was used and its density is 998kg/m3.(3)The outlet is the condition natural outflow,and the flow rate is 1.(4)The boundary condition is no slip at the wall of the passage and the tongue region,and the standard wall function was used with the relative velocity ω being 0 rad/s.(5)The secondary boundary condition based on pressure was used for three-dimensional simulation,and the undeveloped turbulent flow was corrected by the wall function method.(6)The flow between the impeller and volute was disposed by the dynamic-static double coordinate systems.The region of impeller channel was based on the rotating coordinate system,and the angular velocity of rotation at the design condition ω is set to 345.4 rad/s according to the design speed.The angular velocity of rotation at impeller wallω is 0 rad/s.

    3.3. Key performance parameters

    The head of the vortex pump[12],with the unit being M,is calculated through the following expression:

    In the Eq.(1),Po,Pi,νo2and νi2are pump exit static pressure,import static pressure,export speed,and import speed,respectively.The height difference between imports and exports can be considered as zero. And for working medium water, ρ=998kg/m3,g=9.98m/s2.

    The expression about the efficiency of the vortex pump η[12]is

    In the Eq.(2),Q and M are the mass flow and the torque with theunit being m3/h and Nm,respectively,ω is the rotation angular velocity,and the expression is ω=(n×2π)/60.

    Table 3 Working conditions for numerical simulation on vortex pump.

    Fig.3.Gird distribution for double support vortex pump.

    4. Numerical results of the flow field in vortex pump with tongue channel of non-uniform cross section

    4.1. Streamline structure in the channel of vortex pump

    Based on the working conditions for carrying out the numerical simulation of the flow field in vortex pump with variable cross-section tongue channel,which are shown in Table 3,the velocity streamline characteristics in the vortex pump are analyzed by increasing the cross-section width of the tongue channel(E)under the condition of constant height of the tongue channel(G).Fig.5 shows the velocity streamline of the working conditions case 12 and case 19,which are corresponding to the small and large cross-section of the tongue channel respectively.

    Fig.5(a)and(b)are the velocity streamlines starting at the inlet with different sectional area of tongue channel respectively.The fluid enters the circumferential passage of the vortex pump from the inlet along the radial direction,a spiral forward vortex along flow direction is formed and its size and strength increase gradually.The reason is that with the impeller rotating,the fluid at the root of the blade enters the clearance for the periodicity of blade which rotates with the impeller,and the fluid flows along radial direction with rotating centrifugal force,so that the spiral forward vortex is formed along the direction of rotation of the blade.The rotational kinetic energy of the impeller is converted into the pressure energy in the process of the fluid squeezed out periodically from the groove to the channel,and the pressure energy is increased with the influence of impeller when the fluid flows along the circumferential direction.Fig.5(a)shows that the streamline in the inlet region is affected by cross-section of the tongue channel,and the effect of the tongue clearance to backflow is small which results from the leaking fluid in the outlet region flowing into the inlet region.Therefore,the rotational kinetic energy of the impeller is converted to the fluid,so that the head increases.While for the large tongue channel,because the fluid is leaked to the inlet region from the tongue channel,the high pressure backflow has a strong import on the flow in the inlet region and impeller,as shown in Fig.5(b).It makes the fluid not effectively flowing into the impeller to achieve energy conversion,which affects the energy output efficiency of the vortex pump.Fig.6(c)and(d)are the velocity streamline distributions along the horizontal center section(Y=0)for tongue channel of different sectional area respectively.It can be seen that a high speed flow zone cannot be formed at the root of the blade for the tongue channel with larger area.A rotating vortex moving forward helically is not formed.The energy output action is weakened from the impeller to the fluid.And a strong fluid leakage occurs in the tongue clearance,a longitudinal vortex cannot be formed in the blade clearance for the backflow so as not to achieve efficient energy transfer process.While for the tongue channel with small area,there is obvious vortex structure characteristic between the blades and the fluid.It makes the rated capacity cavity to fully reflect the energy transfer of the fluid driven by the rotating blade,so that the fluid enters the flow passage and continuously outputs the pressure energy along the flow direction,thus improving the output head of the vortex pump.In Fig.5,the distribution law of velocity streamlines in different sections shows that:the distribution of tongue clearance has a great influence on the flow field structure inside the vortex pump.By setting a reasonable tongue clearance,the energy transfer between the impeller and the fluid can be effectively improved.

    Fig.5.Velocity streamlines of vortex pump affected by tongue channel with variable cross-section.

    4.2. Pressure and velocity field of the vortex pump

    Fig.6 shows the pressure distribution at the axial cross-section of the vortex pump with tongue channel of two different crosssectional areas which correspond to case 12 and case 19.It can be seen that as the pressure gradually increases along the circumferential direction from the inlet to the outlet,the increase amplitude and law of increase about the pressure are similar basically at the inlet and outlet.Due to the influence of the cross-sectional width of tongue channel,a large area of low pressure is formed at the root of blades which are nearby the inlet,because the fluid leaks to the inlet along the sides of the tongue which is caused by the high pressure at the outlet,and it is so difficult for the fluid at the inlet to enter the impeller as to not form the effective energy exchange,and reduce the energy transfer efficiency. Meanwhile, the tongue channel with large cross-sectional area has disadvantages of not improving the inlet/outlet pressure differential and the head of the cortex pump.For the tongue channel with small cross-sectional area,the inlet/outlet pressure differential is larger,and there is an obvious pressure differential on both sides of the tongue region and increasing pressure gradient along the circumferential direction.These show that the small tongue clearance is beneficial to the energy transfer,and the pressure and efficiency increase normally.

    Fig.6.Pressure distribution of cross-section of vortex pump with non-uniform cross-section tongue channel(z=-2mm).

    Fig.7(a)and(b)are the velocity vector distributions of the variable cross-sectional tongue channel vortex pump,and two cross-sectional locations y=0 and x=0 are shown in Fig.5.By comparing the left figures in Fig.7(a)and(b),a symmetrical longitudinal vortex structure is formed on the central section of the blade clearance,the fluid enters the blade clearance from the outside of the impeller along the root of the blade,and flows toward the outside along the radial direction,then through both sides of the passage,finally goes back to the clearance at the root of the blade again.Therefore,a periodic spiral vortex is formed.For the tongue channel with small cross-sectional area,the cross-sectional velocity gradient is larger,and larger energy transfer is obtained when the fluid flows along the impeller.While for the tongue channel with large cross-sectional area,the radial velocity gradient is smaller,and the energy of the fluid cannot be increased because of the rotation of the impeller.By comparing the right figures in Fig.7(a)and(b),there is high speed leaked fluid at the cross-section of tongue channel with large area,which is bad to improvement of output pressure and energy output.It can be seen from velocity distribution at the axial section which shown in Fig.8(c)and(d),radial sectional velocity gradient is small for the tongue channel with large cross-section area,it is difficult to form a transverse vortex structure in the blade clearance,this is consistent with the distribution laws of the pressure field and velocity streamline,which is not conducive to the improvement of the output head of the vortex pump.

    5. Numerical results of the flow field in vortex pump with tongue channel of uniform cross-section

    5.1. Streamline structure in the channel of vortex pump

    Based on the working conditions for carrying out the numerical simulation of the flow field in vortex pump with uniform crosssection tongue channel,which are shown in Table 3,the velocity streamline characteristic in the vortex pump are analyzed by adjusting the cross-section width(E)and height(G)of the tongue channel under the condition of constant cross-sectional area of the tongue channel.

    Fig.8(a)and(b)show the velocity streamline of the working conditions case12 and case19,which are corresponding to the narrow and wide cross-section of the tongue channel respectively.By comparing Fig.8(a)with Fig.8(b),it can be seen that the streamline distribution standing for higher flow speed is formed on the outside channel of the impeller for the condition of the tongue channel with narrow cross-section,the fluid could not enter the impeller channel from both sides of impeller so as not to form flow vortex,but there is a large amount of fluid entering the both sides of the impeller for the tongue channel with wide cross-section,so that there is a longitudinal vortex structure distribution in impeller channel.Vortex structure for the condition of the wide crosssection has a relatively small impact on the inflow,while the outside channel of impeller for the condition of narrow crosssection has strong impact on the inflow,which weakens the energy exchange process between fluid and impeller.Therefore,it is effective to avoid the situation that the leak flow is formed at the outside diameter clearance of impeller because of the high pressure at the outlet by increasing the width of tongue channel,and then the influence of the outlet on the flow at the inlet region is reduced,so that the energy increasing process of the vortex pump along the circumferential direction is more efficient and continuous,finally the output efficiency and the head are improved.

    5.2. Pressure and velocity field of the vortex pump

    Fig.9(a)and(b)shows the radial cross-section pressure distributions of the vortex pump corresponding to the narrow and wide cross-sectional tongue channels respectively.The tongue width contributes to the improvement of the circumferential pressure gradient.The reduction of the tongue width leads to the outside clearance of the impeller increasing and the radial outside flow having strong leak,so as to affect the flow at the inlet and the improvement of the head and efficiency.

    Fig.10(a)and(b)are the velocity vector characteristics of vortex pumps with two kinds of uniform cross-section tongue channel.The high speed leak occurs in the outside channel of the impeller for the condition of the narrow cross-section,and it is difficult for the fluid in inlet to enter the impeller for energy exchange.While for the condition of the wide cross-section,the longitudinal vortex structure is formed on both sides of the impeller and moves ahead along the circumferential direction in tongue region.This reduces the influence of the high pressure at the outlet on the flowfield in the inlet,and improves the head and efficiency of the vortex pump.According to the circular cross-sectional velocity distribution shown in Fig.10(c)and(d),the velocity gradients on the radial cross-section are similar for the narrow and wide cross-sections.However,the velocity gradient on the circular cross-section of the tongue region is larger,and the transverse vortex structure which moves along the flow direction is formed between the blades.This improves the energy transfer and the head of the vortex pump.

    Fig.7.Velocity distribution of cross-section of vortex pump with non-uniform cross-section tongue channel(y=0 mm、x=0 mm).

    Fig.8.Velocity streamlines of vortex pump with uniform cross-sectional tongue channel.

    Fig.9.Pressure distribution of cross-section of vortex pump with uniform cross-section tongue channel(z=-2mm).

    Fig.10.Velocity distribution of cross-section of vortex pump with uniform cross-section tongue channel(y=0 mm、x=0 mm).

    6. Effect of the cross-section of the tongue on the performance of vortex pump

    Fig.11 shows the effect of cross-sectional width(E)on the head of vortex pump with the tongue channel of the two kinds of crosssection.It can be seen from the Fig.11,as the width increases,the head of vortex pump reduces.And the laws of the influence of the cross-section on the output head are coincident basically.Fig.12 shows the effect of cross-sectional width on the output efficiency of vortex pump with the tongue channel of the two kinds of crosssection.It can be seen from the Fig.12,with the increase of the cross-sectional width of the tongue,the efficiency remains basically unchanged for the condition of variable cross-section,while reduces for the condition of uniform cross-section.

    7. Conclusion

    In this paper,the effects of tongue channel with different crosssection on the internal flowfield and performance of vortex pump were studied under the conditions of the fixed flow channel and impeller,and the conclusions were given as follows.

    (1)The effects of the cross-section characteristic of the tongue channel on the performance of the vortex pump were studied.The variation of cross-section of the tongue channel could affect the vortex structure distribution of the radial inlet flow of the vortex pump in the flow channel,so that the energy exchange efficiency between the impeller and the fluid and the pressure distribution of the vortex pump are changed.

    (2)With the increase of the width of the tongue,the inlet flow of the vortex pump will be more difficult to enter the channel of the impeller.And the spiral forward vortex structure of the vortex pump will shift along large radius.Furthermore,the output efficiency and the head of the vortex pump are both reduced.

    (3)Increasing the cross-section of the tongue channel will reduce the velocity gradient along the circumferential direction and axial pressure gradient.And the radial vortex structure shifts towards the side of the large radius.This is not conductive to the improvement of the head and the efficiency.

    (4)When the width of the tongue channel increases,the head of the vortex pump reduces and the efficiency is not changed basically for the condition of the variable cross-section,while the both reduce for the condition of the uniform crosssection.

    Fig.11.Effect of the width of tongue passage on the head

    Fig.12.Effect of the width of tongue passage on the efficiency.

    According to the results of numerical simulation,this paper only studies the influence of single features including width and height of the rectangular tongue channel of vortex pump.The output performance and internal flow field of vortex pump are obtained,but the influence of many coupling factors in the vortex pump design cannot be taken into account.The subsequent researches include the coupling analysis of design and structure and iterative simulation of structure multi-parameters. At the same time,through the numerical simulation for the double support vortex pump,the fast model of its specific structure can be established and the flow field characteristics can be analyzed.In the study of numerical simulation for the double support vortex pump,the reliability and accuracy of the numerical simulation can be constantly improved by modifying its computing method and boundary conditions based on the existing test data.Then the basic experimental research of the similar model for the flow field structure of the vortex pump and the design of the visualization research method of the flow field are carried out.Finally the law about the influence of the internal flow field on the output performance of the vortex pump is verified by the tests,which provides important data support for fully mastering the optimization design method of the vortex pump.

    Acknowledgment

    We would like to thank the National Natural Science Foundation of China(No.51575420);the Local Service Special Project of the Education Department of Shaanxi Province,China(No.14 JF011);the Special Fund Project of Major Science and Technology Innovation of Shaanxi Province,China(No.2017ZKC01-24,No.2018TP-17);Xi'an Science and Technology project,China(No.201805032YD10CG16(5))for the support.

    日韩欧美在线乱码| 国产精品三级大全| 日日撸夜夜添| 国产精品女同一区二区软件 | 精品久久久久久久久久免费视频| 国产一区二区三区av在线 | 欧美国产日韩亚洲一区| 中文字幕久久专区| 一区二区三区激情视频| 男人舔女人下体高潮全视频| 欧美日韩中文字幕国产精品一区二区三区| 亚洲va在线va天堂va国产| 免费在线观看成人毛片| 久久久精品大字幕| 亚洲精品一卡2卡三卡4卡5卡| 午夜福利高清视频| 免费高清视频大片| 日韩欧美精品免费久久| 久久精品国产亚洲av天美| 真人做人爱边吃奶动态| 国产女主播在线喷水免费视频网站 | 中文字幕高清在线视频| 91久久精品电影网| 亚洲熟妇熟女久久| 1024手机看黄色片| 日韩大尺度精品在线看网址| 国产一级毛片七仙女欲春2| 婷婷色综合大香蕉| 男女边吃奶边做爰视频| 在线观看午夜福利视频| 热99re8久久精品国产| 亚洲人成伊人成综合网2020| 成熟少妇高潮喷水视频| 又黄又爽又刺激的免费视频.| .国产精品久久| 精品人妻视频免费看| a级毛片a级免费在线| 欧美日韩精品成人综合77777| 国产私拍福利视频在线观看| 最后的刺客免费高清国语| 欧美成人免费av一区二区三区| 18+在线观看网站| 免费在线观看影片大全网站| 少妇的逼好多水| 午夜免费激情av| 国产精品久久电影中文字幕| 国产精品98久久久久久宅男小说| 99久久精品热视频| 俺也久久电影网| 免费黄网站久久成人精品| 欧美激情国产日韩精品一区| www日本黄色视频网| 久久久久久大精品| 日本欧美国产在线视频| 亚洲三级黄色毛片| 男女啪啪激烈高潮av片| 欧美日韩中文字幕国产精品一区二区三区| 少妇人妻精品综合一区二区 | 色精品久久人妻99蜜桃| 精品久久久久久久久久久久久| 国产探花极品一区二区| 国产精品无大码| 亚洲精品国产成人久久av| 久久久国产成人免费| 九色国产91popny在线| 欧美不卡视频在线免费观看| 一进一出好大好爽视频| 欧美激情在线99| 女的被弄到高潮叫床怎么办 | 成人无遮挡网站| 91精品国产九色| 亚洲精华国产精华液的使用体验 | 看十八女毛片水多多多| 热99在线观看视频| 久久久久国内视频| 国产 一区 欧美 日韩| 亚洲经典国产精华液单| 国产成年人精品一区二区| 欧美黑人巨大hd| 久久精品国产99精品国产亚洲性色| 成人精品一区二区免费| 无遮挡黄片免费观看| 我要看日韩黄色一级片| 啦啦啦啦在线视频资源| 国产免费男女视频| 男女下面进入的视频免费午夜| 两个人视频免费观看高清| 性欧美人与动物交配| 中国美白少妇内射xxxbb| 99国产极品粉嫩在线观看| 久久久久久大精品| 一级av片app| 在线播放无遮挡| 久久精品人妻少妇| 亚洲无线在线观看| aaaaa片日本免费| 国产在线男女| 国产精品福利在线免费观看| 观看美女的网站| xxxwww97欧美| 亚洲,欧美,日韩| 国产精品久久久久久亚洲av鲁大| 国产精品女同一区二区软件 | 免费av毛片视频| 一区二区三区免费毛片| 欧洲精品卡2卡3卡4卡5卡区| 国产免费一级a男人的天堂| 亚洲自偷自拍三级| 无人区码免费观看不卡| 成人高潮视频无遮挡免费网站| 日韩欧美精品v在线| 我要看日韩黄色一级片| 桃色一区二区三区在线观看| 成人高潮视频无遮挡免费网站| 国产三级在线视频| 国内毛片毛片毛片毛片毛片| 美女高潮的动态| 色5月婷婷丁香| 久久热精品热| 免费看美女性在线毛片视频| 在线免费观看不下载黄p国产 | 一区二区三区四区激情视频 | 亚洲欧美清纯卡通| 国产在线精品亚洲第一网站| 日韩中文字幕欧美一区二区| 婷婷六月久久综合丁香| 俺也久久电影网| av在线老鸭窝| ponron亚洲| 变态另类丝袜制服| 乱系列少妇在线播放| 日日夜夜操网爽| 国产女主播在线喷水免费视频网站 | 成人国产麻豆网| 日韩av在线大香蕉| 色尼玛亚洲综合影院| 在线观看66精品国产| 他把我摸到了高潮在线观看| 春色校园在线视频观看| 日韩大尺度精品在线看网址| 亚洲精品一区av在线观看| 国产一区二区三区在线臀色熟女| 亚洲不卡免费看| av.在线天堂| 国产成人福利小说| 午夜福利成人在线免费观看| 12—13女人毛片做爰片一| a级毛片a级免费在线| 99久国产av精品| 日本成人三级电影网站| 露出奶头的视频| 一级av片app| 国产高清不卡午夜福利| 露出奶头的视频| 精品日产1卡2卡| 日韩中文字幕欧美一区二区| 我要搜黄色片| 最近最新中文字幕大全电影3| 国产高潮美女av| 免费大片18禁| 成人av一区二区三区在线看| 99国产精品一区二区蜜桃av| 人妻丰满熟妇av一区二区三区| 久久久午夜欧美精品| 国产高清有码在线观看视频| 国产精品永久免费网站| 一区二区三区高清视频在线| 成人午夜高清在线视频| 亚洲三级黄色毛片| 成人三级黄色视频| 日韩亚洲欧美综合| 人妻制服诱惑在线中文字幕| 天堂√8在线中文| xxxwww97欧美| 成年女人看的毛片在线观看| 日韩强制内射视频| 欧美另类亚洲清纯唯美| 国产精品久久久久久亚洲av鲁大| 久久久久久久精品吃奶| 最近最新中文字幕大全电影3| 国产免费男女视频| 18+在线观看网站| avwww免费| 亚洲av一区综合| 午夜爱爱视频在线播放| 99久久九九国产精品国产免费| av在线观看视频网站免费| 午夜福利在线观看免费完整高清在 | 22中文网久久字幕| 国产精品99久久久久久久久| 国产在线男女| 亚洲第一电影网av| 午夜爱爱视频在线播放| 久久精品影院6| 婷婷六月久久综合丁香| av在线观看视频网站免费| 国产亚洲精品久久久久久毛片| 国产精品人妻久久久久久| 日本在线视频免费播放| 能在线免费观看的黄片| 免费人成视频x8x8入口观看| 91在线精品国自产拍蜜月| 九九久久精品国产亚洲av麻豆| 国产高清视频在线观看网站| 欧美性感艳星| 啪啪无遮挡十八禁网站| 午夜爱爱视频在线播放| 午夜福利成人在线免费观看| 啦啦啦观看免费观看视频高清| 中文字幕高清在线视频| 日韩精品中文字幕看吧| 在线观看免费视频日本深夜| 99久久久亚洲精品蜜臀av| 日韩人妻高清精品专区| 亚洲自拍偷在线| 国产精品一及| 午夜a级毛片| 99在线人妻在线中文字幕| 亚洲最大成人av| 中亚洲国语对白在线视频| 最近中文字幕高清免费大全6 | 国产成人一区二区在线| 免费看日本二区| 精品一区二区三区视频在线| 99久国产av精品| netflix在线观看网站| 国产乱人视频| 九九热线精品视视频播放| 欧美高清性xxxxhd video| 亚洲av五月六月丁香网| 男人舔女人下体高潮全视频| 高清毛片免费观看视频网站| 亚洲在线自拍视频| 国产亚洲精品久久久久久毛片| 俄罗斯特黄特色一大片| 97超级碰碰碰精品色视频在线观看| 精品久久久久久久人妻蜜臀av| 中文字幕av成人在线电影| 国产一区二区三区在线臀色熟女| 亚洲久久久久久中文字幕| 亚洲色图av天堂| 午夜福利欧美成人| av女优亚洲男人天堂| 天堂av国产一区二区熟女人妻| 精品久久久久久久久亚洲 | 亚洲三级黄色毛片| 3wmmmm亚洲av在线观看| 禁无遮挡网站| 尤物成人国产欧美一区二区三区| 国产精品福利在线免费观看| 国产成人影院久久av| 日韩欧美在线二视频| 欧美+亚洲+日韩+国产| 噜噜噜噜噜久久久久久91| 夜夜夜夜夜久久久久| 欧美日韩综合久久久久久 | 亚洲五月天丁香| 又爽又黄a免费视频| 在线天堂最新版资源| 国产av不卡久久| bbb黄色大片| 99精品在免费线老司机午夜| 我的女老师完整版在线观看| 人妻少妇偷人精品九色| 欧美日韩国产亚洲二区| 国产不卡一卡二| bbb黄色大片| 亚洲午夜理论影院| 精品国产三级普通话版| 在现免费观看毛片| 91在线观看av| 床上黄色一级片| 美女 人体艺术 gogo| 欧美高清成人免费视频www| 无遮挡黄片免费观看| 久久久久久久精品吃奶| 又黄又爽又刺激的免费视频.| 久久精品国产自在天天线| 一级二级三级毛片免费看| 嘟嘟电影网在线观看| 九草在线视频观看| 国产欧美另类精品又又久久亚洲欧美| 最近的中文字幕免费完整| 日韩一区二区视频免费看| 精品国产三级普通话版| 日韩在线高清观看一区二区三区| 一本—道久久a久久精品蜜桃钙片| 老司机影院成人| 日本欧美视频一区| 在线看a的网站| 中国国产av一级| 性色av一级| 精品酒店卫生间| 色吧在线观看| 狠狠精品人妻久久久久久综合| 久久国产精品大桥未久av | 午夜激情久久久久久久| 国产中年淑女户外野战色| 色5月婷婷丁香| 国产真实伦视频高清在线观看| 久久久久精品性色| 蜜桃在线观看..| 国产亚洲5aaaaa淫片| 成人漫画全彩无遮挡| 成年人午夜在线观看视频| 成年女人在线观看亚洲视频| 日本午夜av视频| 欧美成人一区二区免费高清观看| 亚洲第一av免费看| 男女啪啪激烈高潮av片| 国产老妇伦熟女老妇高清| 蜜桃在线观看..| 看非洲黑人一级黄片| 丰满少妇做爰视频| 欧美激情国产日韩精品一区| 色婷婷av一区二区三区视频| 亚洲精品第二区| 日本黄色日本黄色录像| av免费在线看不卡| 熟妇人妻不卡中文字幕| 成人毛片a级毛片在线播放| 免费av中文字幕在线| 男人狂女人下面高潮的视频| 久久久久久久久久人人人人人人| 久久青草综合色| 国产在线一区二区三区精| 久久人妻熟女aⅴ| 少妇的逼好多水| 99热全是精品| 成人亚洲精品一区在线观看 | 一区二区三区乱码不卡18| 国产精品麻豆人妻色哟哟久久| 成人毛片a级毛片在线播放| 精品少妇久久久久久888优播| 熟女av电影| 亚洲欧洲国产日韩| 国产淫片久久久久久久久| 亚洲国产精品999| 3wmmmm亚洲av在线观看| 国产一级毛片在线| 夜夜看夜夜爽夜夜摸| 精品久久久久久久久av| 精品99又大又爽又粗少妇毛片| 国产午夜精品久久久久久一区二区三区| av卡一久久| 国产淫语在线视频| 日本午夜av视频| 一本久久精品| 亚洲欧美一区二区三区黑人 | 国产一区二区在线观看日韩| 大香蕉97超碰在线| 欧美日韩一区二区视频在线观看视频在线| 一级二级三级毛片免费看| 极品少妇高潮喷水抽搐| 天天躁夜夜躁狠狠久久av| 啦啦啦中文免费视频观看日本| 新久久久久国产一级毛片| 亚洲av欧美aⅴ国产| 成人影院久久| 九九爱精品视频在线观看| 男人舔奶头视频| 久久6这里有精品| 成年女人在线观看亚洲视频| 国产精品人妻久久久影院| 天堂俺去俺来也www色官网| 国产高清不卡午夜福利| 成人18禁高潮啪啪吃奶动态图 | 成人毛片60女人毛片免费| 美女cb高潮喷水在线观看| 熟女av电影| 亚洲伊人久久精品综合| 少妇人妻精品综合一区二区| 亚洲精品视频女| 久久 成人 亚洲| 精品视频人人做人人爽| 五月伊人婷婷丁香| 国产深夜福利视频在线观看| 亚洲精品日韩在线中文字幕| 国产在视频线精品| 亚洲成人av在线免费| 久久久久人妻精品一区果冻| 老师上课跳d突然被开到最大视频| 午夜福利网站1000一区二区三区| 狂野欧美白嫩少妇大欣赏| av福利片在线观看| 欧美国产精品一级二级三级 | 亚洲三级黄色毛片| 婷婷色综合大香蕉| 国产 精品1| 亚洲av欧美aⅴ国产| 免费观看无遮挡的男女| 中国美白少妇内射xxxbb| 久久久久人妻精品一区果冻| 黄色欧美视频在线观看| 日韩制服骚丝袜av| 国产中年淑女户外野战色| 欧美区成人在线视频| 激情五月婷婷亚洲| 丝袜喷水一区| 777米奇影视久久| 国产欧美日韩一区二区三区在线 | 国产精品嫩草影院av在线观看| 在线天堂最新版资源| 美女xxoo啪啪120秒动态图| 新久久久久国产一级毛片| 欧美日韩在线观看h| 超碰97精品在线观看| 永久免费av网站大全| 免费观看性生交大片5| 联通29元200g的流量卡| 狂野欧美激情性xxxx在线观看| 老女人水多毛片| 国产精品国产三级专区第一集| 一区二区三区乱码不卡18| 人人妻人人看人人澡| 一级二级三级毛片免费看| 国产一区有黄有色的免费视频| 日本av免费视频播放| 亚洲av不卡在线观看| 亚洲精品日韩在线中文字幕| 国产69精品久久久久777片| 国产精品一区二区在线观看99| 一级毛片我不卡| 女性生殖器流出的白浆| av免费在线看不卡| 新久久久久国产一级毛片| 亚洲精品,欧美精品| 在线免费十八禁| 国产探花极品一区二区| 丰满乱子伦码专区| 亚洲人成网站高清观看| 只有这里有精品99| 在线亚洲精品国产二区图片欧美 | 国产伦精品一区二区三区四那| 亚洲av.av天堂| 极品教师在线视频| 日韩制服骚丝袜av| 久热久热在线精品观看| 亚洲av中文av极速乱| 亚洲av综合色区一区| 91久久精品电影网| 国内精品宾馆在线| 亚洲中文av在线| 欧美性感艳星| 国产欧美日韩精品一区二区| 亚洲国产欧美在线一区| 蜜桃久久精品国产亚洲av| 日日啪夜夜撸| 国产成人aa在线观看| 日韩国内少妇激情av| 美女内射精品一级片tv| 不卡视频在线观看欧美| 97在线视频观看| 亚洲欧美日韩另类电影网站 | 日日啪夜夜撸| 蜜桃亚洲精品一区二区三区| 中文字幕av成人在线电影| 国产男人的电影天堂91| 国产无遮挡羞羞视频在线观看| 国产乱人视频| 免费观看的影片在线观看| 亚洲性久久影院| 欧美成人a在线观看| 97超碰精品成人国产| 免费久久久久久久精品成人欧美视频 | 性高湖久久久久久久久免费观看| 免费大片18禁| 高清毛片免费看| 日韩大片免费观看网站| 欧美极品一区二区三区四区| av播播在线观看一区| 国产欧美日韩一区二区三区在线 | 亚洲美女搞黄在线观看| 汤姆久久久久久久影院中文字幕| 2021少妇久久久久久久久久久| 国内少妇人妻偷人精品xxx网站| 国产免费一区二区三区四区乱码| 国产精品.久久久| 美女xxoo啪啪120秒动态图| 三级经典国产精品| 综合色丁香网| 人人妻人人爽人人添夜夜欢视频 | 中国三级夫妇交换| 99热这里只有精品一区| 日产精品乱码卡一卡2卡三| 国产精品99久久99久久久不卡 | 日韩大片免费观看网站| 蜜桃亚洲精品一区二区三区| 国产乱人视频| 水蜜桃什么品种好| 91久久精品电影网| 久久久久精品性色| 精品人妻偷拍中文字幕| 国产老妇伦熟女老妇高清| 国产精品成人在线| 亚洲精品成人av观看孕妇| 欧美三级亚洲精品| 亚洲,一卡二卡三卡| 久久人妻熟女aⅴ| 美女福利国产在线 | 性高湖久久久久久久久免费观看| 一级二级三级毛片免费看| .国产精品久久| 午夜免费鲁丝| 日韩成人伦理影院| 午夜老司机福利剧场| 九草在线视频观看| 日本av手机在线免费观看| kizo精华| 天天躁日日操中文字幕| 成人高潮视频无遮挡免费网站| 舔av片在线| 国产男女超爽视频在线观看| 一区二区三区免费毛片| 国产精品.久久久| av女优亚洲男人天堂| 成年美女黄网站色视频大全免费 | 精品酒店卫生间| 亚洲精品国产色婷婷电影| 三级经典国产精品| 91精品国产国语对白视频| 久久av网站| 午夜福利视频精品| 亚洲经典国产精华液单| 18禁在线无遮挡免费观看视频| 22中文网久久字幕| 狂野欧美白嫩少妇大欣赏| 日韩av不卡免费在线播放| 99热这里只有精品一区| 搡女人真爽免费视频火全软件| 蜜桃久久精品国产亚洲av| 欧美亚洲 丝袜 人妻 在线| 国产乱人偷精品视频| 一区二区av电影网| 久久久久久人妻| 夫妻午夜视频| 日本欧美视频一区| av又黄又爽大尺度在线免费看| 另类亚洲欧美激情| 亚洲内射少妇av| 美女中出高潮动态图| 精品一品国产午夜福利视频| 99热网站在线观看| 亚洲美女搞黄在线观看| 精品久久久久久久久av| 亚洲丝袜综合中文字幕| 日韩av免费高清视频| 国产一区二区三区综合在线观看 | av不卡在线播放| 国产精品人妻久久久久久| 毛片一级片免费看久久久久| 人妻制服诱惑在线中文字幕| 国产黄片视频在线免费观看| 欧美xxⅹ黑人| 欧美成人一区二区免费高清观看| 国产男女超爽视频在线观看| 美女内射精品一级片tv| 久久久午夜欧美精品| 一级黄片播放器| 这个男人来自地球电影免费观看 | 在线观看一区二区三区激情| 最近最新中文字幕免费大全7| 国产无遮挡羞羞视频在线观看| 亚州av有码| 美女中出高潮动态图| 国产国拍精品亚洲av在线观看| 日本黄色日本黄色录像| 亚洲欧美精品自产自拍| 久久久色成人| 精品人妻偷拍中文字幕| 国产精品国产三级国产专区5o| 亚洲av男天堂| 一个人免费看片子| 久久韩国三级中文字幕| 亚洲精品亚洲一区二区| 婷婷色av中文字幕| 黄片无遮挡物在线观看| 成人综合一区亚洲| 欧美另类一区| 丝袜喷水一区| 最近中文字幕高清免费大全6| 看免费成人av毛片| 啦啦啦视频在线资源免费观看| 纵有疾风起免费观看全集完整版| 国产精品一区www在线观看| 人妻夜夜爽99麻豆av| 免费黄网站久久成人精品| 亚洲av成人精品一二三区| 老熟女久久久| 美女xxoo啪啪120秒动态图| 日韩欧美一区视频在线观看 | 美女国产视频在线观看| 尤物成人国产欧美一区二区三区| 国产乱人视频| 伊人久久国产一区二区| 国产欧美另类精品又又久久亚洲欧美| 久久国产精品男人的天堂亚洲 | 男女国产视频网站| 最近中文字幕2019免费版| 高清欧美精品videossex| 女人久久www免费人成看片| 丝袜脚勾引网站| 国产爱豆传媒在线观看| 免费黄色在线免费观看| 中文天堂在线官网| 午夜激情福利司机影院| 亚洲av中文av极速乱| 99热6这里只有精品| 国产欧美另类精品又又久久亚洲欧美| 精品一品国产午夜福利视频| 黑人高潮一二区| 国产又色又爽无遮挡免| 国产成人a∨麻豆精品| 国产伦精品一区二区三区四那| 日韩成人av中文字幕在线观看|