• <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.

    婷婷色麻豆天堂久久| 热99国产精品久久久久久7| 成年免费大片在线观看| kizo精华| 日产精品乱码卡一卡2卡三| 一本色道久久久久久精品综合| 欧美精品国产亚洲| 啦啦啦在线观看免费高清www| av一本久久久久| 国产成人a∨麻豆精品| 超碰97精品在线观看| 少妇的逼水好多| 色网站视频免费| 日韩伦理黄色片| 欧美日韩视频高清一区二区三区二| 欧美精品亚洲一区二区| 99久久精品一区二区三区| 亚洲av二区三区四区| 91aial.com中文字幕在线观看| 免费av中文字幕在线| 六月丁香七月| 亚洲精品一区蜜桃| 久久女婷五月综合色啪小说| 日本-黄色视频高清免费观看| 久久99精品国语久久久| 免费av中文字幕在线| 国产片特级美女逼逼视频| 国产大屁股一区二区在线视频| 国产白丝娇喘喷水9色精品| 日本与韩国留学比较| 精品久久久久久电影网| 亚洲va在线va天堂va国产| 亚洲美女搞黄在线观看| 久久久久久久久久人人人人人人| 亚洲av男天堂| 国内精品宾馆在线| 成年女人在线观看亚洲视频| 在线观看av片永久免费下载| 亚洲精品国产成人久久av| 中文乱码字字幕精品一区二区三区| 国产亚洲av片在线观看秒播厂| 亚洲色图av天堂| 最黄视频免费看| 一级二级三级毛片免费看| 最近的中文字幕免费完整| 国产免费又黄又爽又色| 国产av码专区亚洲av| 日韩制服骚丝袜av| 99久久精品国产国产毛片| 国产乱人偷精品视频| 国产成人精品福利久久| 九九久久精品国产亚洲av麻豆| 亚洲aⅴ乱码一区二区在线播放| 久热久热在线精品观看| 亚洲久久久国产精品| 在线观看免费日韩欧美大片 | 最近最新中文字幕免费大全7| 毛片一级片免费看久久久久| 久久久a久久爽久久v久久| 九九在线视频观看精品| 777米奇影视久久| 日本一二三区视频观看| 国产成人a区在线观看| 国产熟女欧美一区二区| 日韩中文字幕视频在线看片 | 亚洲激情五月婷婷啪啪| 偷拍熟女少妇极品色| a级毛色黄片| 欧美高清成人免费视频www| 男女啪啪激烈高潮av片| 天堂8中文在线网| 国产亚洲一区二区精品| 久久久欧美国产精品| 男人舔奶头视频| 成人特级av手机在线观看| 最近的中文字幕免费完整| 80岁老熟妇乱子伦牲交| 久久99蜜桃精品久久| 伦理电影免费视频| 97热精品久久久久久| 国产人妻一区二区三区在| 少妇人妻精品综合一区二区| kizo精华| 国产精品成人在线| 久久6这里有精品| 欧美 日韩 精品 国产| 久久久久国产网址| 高清日韩中文字幕在线| 美女cb高潮喷水在线观看| 国产 一区 欧美 日韩| 夜夜看夜夜爽夜夜摸| 18禁裸乳无遮挡动漫免费视频| 久久国内精品自在自线图片| 久久久精品免费免费高清| 日韩一本色道免费dvd| 国国产精品蜜臀av免费| 午夜激情福利司机影院| 亚洲欧美日韩东京热| 美女主播在线视频| 天天躁日日操中文字幕| 亚洲av在线观看美女高潮| 亚洲成人av在线免费| 日本wwww免费看| 国模一区二区三区四区视频| 国产欧美日韩精品一区二区| 国产黄频视频在线观看| 啦啦啦中文免费视频观看日本| 久久精品久久久久久久性| 午夜福利网站1000一区二区三区| 中文精品一卡2卡3卡4更新| 麻豆成人午夜福利视频| 人人妻人人看人人澡| 自拍偷自拍亚洲精品老妇| 观看美女的网站| 欧美 日韩 精品 国产| 国产精品一区二区三区四区免费观看| 各种免费的搞黄视频| 一本久久精品| 午夜免费观看性视频| 午夜日本视频在线| 啦啦啦视频在线资源免费观看| 久久精品国产亚洲av天美| 黄色一级大片看看| 精品午夜福利在线看| 日韩亚洲欧美综合| 日韩成人伦理影院| 3wmmmm亚洲av在线观看| 久久国内精品自在自线图片| 熟女av电影| 亚洲熟女精品中文字幕| 国产精品久久久久久av不卡| 中文字幕久久专区| 哪个播放器可以免费观看大片| 日日摸夜夜添夜夜爱| 亚洲精品乱码久久久v下载方式| 欧美日韩综合久久久久久| 日韩三级伦理在线观看| 亚洲人与动物交配视频| 久久av网站| 美女主播在线视频| 五月天丁香电影| 少妇的逼水好多| 精品人妻视频免费看| 日日摸夜夜添夜夜爱| 老女人水多毛片| 男女免费视频国产| 我要看日韩黄色一级片| 精品亚洲乱码少妇综合久久| 久久精品人妻少妇| 国产av国产精品国产| 亚洲国产精品国产精品| 亚洲精品日韩在线中文字幕| 男男h啪啪无遮挡| 丝袜喷水一区| 黄片无遮挡物在线观看| 亚洲va在线va天堂va国产| 秋霞伦理黄片| 久久精品久久久久久噜噜老黄| 丰满少妇做爰视频| 中国美白少妇内射xxxbb| 日韩伦理黄色片| 久久国内精品自在自线图片| 国产精品不卡视频一区二区| 视频区图区小说| 99久久精品国产国产毛片| 国产成人精品一,二区| 久久人妻熟女aⅴ| 日韩伦理黄色片| 亚洲av电影在线观看一区二区三区| 国产黄色视频一区二区在线观看| 国产色爽女视频免费观看| 看十八女毛片水多多多| 国产一区二区三区av在线| 大码成人一级视频| 在线观看国产h片| 日本一二三区视频观看| 久久人人爽人人片av| 久久久久久久久久久丰满| 国产精品国产av在线观看| 亚洲最大成人中文| 精品人妻偷拍中文字幕| 久久久久精品久久久久真实原创| 国产精品免费大片| 国产 精品1| 亚洲久久久国产精品| 亚洲综合色惰| 免费看av在线观看网站| 亚洲精品aⅴ在线观看| 亚洲国产精品一区三区| 97在线人人人人妻| 久久99精品国语久久久| 国产老妇伦熟女老妇高清| 国产成人一区二区在线| 国产亚洲91精品色在线| 99热这里只有精品一区| 亚洲精品日韩在线中文字幕| 99久久精品国产国产毛片| 欧美日本视频| 菩萨蛮人人尽说江南好唐韦庄| 麻豆成人午夜福利视频| 久久国产精品大桥未久av | 精品久久国产蜜桃| 九九爱精品视频在线观看| 日本黄色日本黄色录像| 少妇人妻 视频| 亚洲国产色片| 老司机影院毛片| 免费av中文字幕在线| 欧美高清性xxxxhd video| 男的添女的下面高潮视频| 久久青草综合色| 精品久久久噜噜| 欧美日韩视频精品一区| 下体分泌物呈黄色| 欧美日韩综合久久久久久| 噜噜噜噜噜久久久久久91| 国产乱人偷精品视频| 午夜日本视频在线| 一区二区三区乱码不卡18| 午夜免费男女啪啪视频观看| 色婷婷av一区二区三区视频| 国产一级毛片在线| 国产爱豆传媒在线观看| 夜夜骑夜夜射夜夜干| 亚洲欧美成人精品一区二区| 五月玫瑰六月丁香| 熟妇人妻不卡中文字幕| 亚洲激情五月婷婷啪啪| 国产精品久久久久成人av| 亚洲久久久国产精品| 亚洲av福利一区| 美女视频免费永久观看网站| 午夜日本视频在线| 国产片特级美女逼逼视频| 国产在线免费精品| 成人特级av手机在线观看| 久久精品国产鲁丝片午夜精品| 国产精品国产三级专区第一集| 精品少妇久久久久久888优播| 最近最新中文字幕大全电影3| 亚洲欧美日韩无卡精品| 一本一本综合久久| 国产av精品麻豆| 免费大片18禁| 一级爰片在线观看| 精品国产三级普通话版| 国产成人aa在线观看| 中国美白少妇内射xxxbb| 成人黄色视频免费在线看| 亚洲精品久久午夜乱码| 少妇人妻精品综合一区二区| 国产亚洲av片在线观看秒播厂| www.av在线官网国产| 亚洲av国产av综合av卡| 久久99热这里只有精品18| 国产亚洲午夜精品一区二区久久| 久久久国产一区二区| 国产免费一区二区三区四区乱码| 国产无遮挡羞羞视频在线观看| 欧美97在线视频| 亚洲精品久久久久久婷婷小说| 18禁在线播放成人免费| 亚洲av综合色区一区| 成年美女黄网站色视频大全免费 | 亚洲精品国产成人久久av| videos熟女内射| 少妇人妻一区二区三区视频| 国产人妻一区二区三区在| 欧美最新免费一区二区三区| 嫩草影院新地址| 日韩欧美一区视频在线观看 | 精品久久久精品久久久| 女人久久www免费人成看片| 亚洲精品国产色婷婷电影| 国产精品麻豆人妻色哟哟久久| 少妇丰满av| 人妻少妇偷人精品九色| 男女免费视频国产| 好男人视频免费观看在线| 尾随美女入室| 国产亚洲最大av| 在线观看人妻少妇| 在线观看一区二区三区| 美女高潮的动态| 在线看a的网站| 亚洲国产精品一区三区| 国产精品国产av在线观看| 一本—道久久a久久精品蜜桃钙片| 欧美日韩亚洲高清精品| 深夜a级毛片| 国产精品一区www在线观看| .国产精品久久| 久久综合国产亚洲精品| kizo精华| 久久ye,这里只有精品| 99热这里只有是精品50| 少妇人妻 视频| 欧美3d第一页| av播播在线观看一区| 国产 精品1| 国产精品久久久久久久电影| 一区二区三区免费毛片| 国产伦理片在线播放av一区| 一边亲一边摸免费视频| 中文字幕久久专区| 自拍欧美九色日韩亚洲蝌蚪91 | 99re6热这里在线精品视频| 国产国拍精品亚洲av在线观看| 精品亚洲成a人片在线观看 | 人妻 亚洲 视频| 亚州av有码| 久久久久久久国产电影| 国产一区二区三区av在线| 秋霞在线观看毛片| 亚洲精品视频女| 中文乱码字字幕精品一区二区三区| 美女脱内裤让男人舔精品视频| 啦啦啦视频在线资源免费观看| 伦精品一区二区三区| 黑人猛操日本美女一级片| 又黄又爽又刺激的免费视频.| 九九久久精品国产亚洲av麻豆| 纯流量卡能插随身wifi吗| 日本与韩国留学比较| 少妇 在线观看| 国产高潮美女av| 日本欧美国产在线视频| 亚洲第一区二区三区不卡| videossex国产| 久久久久久久久久久免费av| 六月丁香七月| 成人亚洲欧美一区二区av| 国产免费视频播放在线视频| 久久久久久久久久久免费av| 99久久综合免费| 岛国毛片在线播放| 亚洲aⅴ乱码一区二区在线播放| 精品少妇黑人巨大在线播放| 91在线精品国自产拍蜜月| 国产伦在线观看视频一区| 国产精品免费大片| 一级毛片aaaaaa免费看小| 自拍欧美九色日韩亚洲蝌蚪91 | 国产免费一级a男人的天堂| 插阴视频在线观看视频| 黑丝袜美女国产一区| 国产成人a∨麻豆精品| 岛国毛片在线播放| 天天躁日日操中文字幕| 欧美日韩综合久久久久久| 美女内射精品一级片tv| 欧美区成人在线视频| 亚洲婷婷狠狠爱综合网| 成人亚洲欧美一区二区av| 精品久久久久久久末码| 全区人妻精品视频| 日本欧美国产在线视频| 中文天堂在线官网| 国产成人aa在线观看| 国产永久视频网站| 午夜激情福利司机影院| 国产精品成人在线| 欧美成人精品欧美一级黄| 国产午夜精品久久久久久一区二区三区| 秋霞在线观看毛片| 中文欧美无线码| 欧美3d第一页| 欧美bdsm另类| 尤物成人国产欧美一区二区三区| 亚州av有码| 国产成人精品一,二区| 国产精品人妻久久久影院| 免费黄频网站在线观看国产| 精品午夜福利在线看| 免费少妇av软件| 国产精品人妻久久久影院| 国产日韩欧美在线精品| 777米奇影视久久| 国产成人a区在线观看| 国产精品人妻久久久影院| 免费少妇av软件| 又粗又硬又长又爽又黄的视频| 精品亚洲成国产av| 一本色道久久久久久精品综合| 免费黄频网站在线观看国产| 各种免费的搞黄视频| 99久久人妻综合| 美女国产视频在线观看| av线在线观看网站| 亚洲怡红院男人天堂| 成年人午夜在线观看视频| 在线看a的网站| 在线免费十八禁| 国产免费一级a男人的天堂| 午夜福利视频精品| 国产精品.久久久| 91精品伊人久久大香线蕉| 久久av网站| 久久久久久久精品精品| 精品酒店卫生间| 欧美日韩视频高清一区二区三区二| 内射极品少妇av片p| 国国产精品蜜臀av免费| 97超碰精品成人国产| 一级毛片久久久久久久久女| 18+在线观看网站| 精品久久国产蜜桃| 精品一品国产午夜福利视频| 久久久国产一区二区| 午夜视频国产福利| 欧美一区二区亚洲| 免费看光身美女| 国产精品爽爽va在线观看网站| 亚洲国产精品999| 国产男人的电影天堂91| 精品视频人人做人人爽| 久久影院123| 在线精品无人区一区二区三 | 中文字幕久久专区| 亚洲精品国产色婷婷电影| 秋霞伦理黄片| 一级毛片久久久久久久久女| 赤兔流量卡办理| 深爱激情五月婷婷| 五月开心婷婷网| 嫩草影院新地址| 亚洲经典国产精华液单| 一二三四中文在线观看免费高清| 全区人妻精品视频| 亚洲国产色片| 男的添女的下面高潮视频| 久久精品国产亚洲网站| 久久热精品热| 91精品国产九色| 99久久精品热视频| 久久久久久久大尺度免费视频| 99久久精品热视频| 久久精品国产自在天天线| 国产在线男女| 少妇裸体淫交视频免费看高清| 黄色配什么色好看| 亚洲综合精品二区| 成人毛片a级毛片在线播放| 亚洲精品第二区| 国产精品久久久久久久电影| 少妇猛男粗大的猛烈进出视频| 99re6热这里在线精品视频| 成人18禁高潮啪啪吃奶动态图 | 日本vs欧美在线观看视频 | 男人和女人高潮做爰伦理| 黄色日韩在线| 亚洲精品国产色婷婷电影| 久久婷婷青草| 亚洲成人中文字幕在线播放| 少妇熟女欧美另类| 国产高清不卡午夜福利| 2018国产大陆天天弄谢| 干丝袜人妻中文字幕| 免费观看av网站的网址| 啦啦啦视频在线资源免费观看| 国产精品不卡视频一区二区| 一级a做视频免费观看| 内地一区二区视频在线| 久久99蜜桃精品久久| 国产精品熟女久久久久浪| 美女脱内裤让男人舔精品视频| 最近中文字幕高清免费大全6| videos熟女内射| 深夜a级毛片| 舔av片在线| 在线观看一区二区三区激情| 夜夜骑夜夜射夜夜干| 高清视频免费观看一区二区| 国产精品爽爽va在线观看网站| 少妇精品久久久久久久| 亚洲国产精品专区欧美| 美女高潮的动态| 自拍偷自拍亚洲精品老妇| 18禁裸乳无遮挡免费网站照片| 女性生殖器流出的白浆| 国产老妇伦熟女老妇高清| 边亲边吃奶的免费视频| 国产精品伦人一区二区| 亚洲欧洲国产日韩| 一级毛片aaaaaa免费看小| 色吧在线观看| 熟女av电影| 亚洲高清免费不卡视频| 国产有黄有色有爽视频| 在线免费观看不下载黄p国产| 久久久午夜欧美精品| 国产91av在线免费观看| 777米奇影视久久| 日韩一区二区视频免费看| 欧美一区二区亚洲| 看免费成人av毛片| 亚洲四区av| 日本与韩国留学比较| 国产精品久久久久久久久免| 亚洲av男天堂| 联通29元200g的流量卡| 亚洲精品久久午夜乱码| av.在线天堂| 在线 av 中文字幕| 日本av免费视频播放| 午夜精品国产一区二区电影| 高清av免费在线| 伊人久久精品亚洲午夜| 精品久久久久久久末码| 亚洲精品视频女| 亚洲成人一二三区av| 日韩大片免费观看网站| 免费看不卡的av| 视频中文字幕在线观看| 天堂中文最新版在线下载| 日韩 亚洲 欧美在线| 日韩伦理黄色片| 亚洲三级黄色毛片| 狂野欧美激情性xxxx在线观看| 黄色配什么色好看| 一本久久精品| av播播在线观看一区| 少妇 在线观看| 国产真实伦视频高清在线观看| 观看免费一级毛片| 欧美国产精品一级二级三级 | 22中文网久久字幕| 国产 精品1| 久久99精品国语久久久| 国产大屁股一区二区在线视频| 777米奇影视久久| 不卡视频在线观看欧美| 欧美日韩国产mv在线观看视频 | 久久久久国产精品人妻一区二区| 久久久久精品久久久久真实原创| 精品亚洲成国产av| 欧美性感艳星| 三级经典国产精品| 最近的中文字幕免费完整| 国产精品偷伦视频观看了| 国产亚洲av片在线观看秒播厂| 国产国拍精品亚洲av在线观看| 免费大片黄手机在线观看| 欧美成人一区二区免费高清观看| 久久久午夜欧美精品| 女的被弄到高潮叫床怎么办| 自拍欧美九色日韩亚洲蝌蚪91 | 91狼人影院| 一区二区三区乱码不卡18| 身体一侧抽搐| 婷婷色av中文字幕| av又黄又爽大尺度在线免费看| 亚洲av中文字字幕乱码综合| 大片免费播放器 马上看| 高清黄色对白视频在线免费看 | 人人妻人人看人人澡| 免费久久久久久久精品成人欧美视频 | 久久久久久久久久成人| 亚洲内射少妇av| 老师上课跳d突然被开到最大视频| 日韩不卡一区二区三区视频在线| 欧美3d第一页| 亚洲精品,欧美精品| 国产成人免费无遮挡视频| 一区二区三区乱码不卡18| 寂寞人妻少妇视频99o| 午夜福利影视在线免费观看| 国产一区二区三区综合在线观看 | 大码成人一级视频| 久久国产精品男人的天堂亚洲 | 日韩中字成人| 国产精品蜜桃在线观看| 久久99精品国语久久久| 日韩 亚洲 欧美在线| 国产亚洲5aaaaa淫片| 水蜜桃什么品种好| 日韩一本色道免费dvd| 久久99热6这里只有精品| 男的添女的下面高潮视频| 免费观看的影片在线观看| 亚洲第一区二区三区不卡| 亚洲国产精品成人久久小说| 国产av码专区亚洲av| 王馨瑶露胸无遮挡在线观看| 九九在线视频观看精品| 王馨瑶露胸无遮挡在线观看| 不卡视频在线观看欧美| 久久韩国三级中文字幕| 看免费成人av毛片| av国产久精品久网站免费入址| 看十八女毛片水多多多| 男男h啪啪无遮挡| 久久人妻熟女aⅴ| 日本欧美国产在线视频| 国产精品一区二区性色av| 夜夜骑夜夜射夜夜干| av.在线天堂| 麻豆精品久久久久久蜜桃| 大码成人一级视频| 91狼人影院| 亚洲欧美中文字幕日韩二区| 精华霜和精华液先用哪个| 午夜免费观看性视频| av在线播放精品| 久久热精品热| 97精品久久久久久久久久精品| 在线看a的网站| 人妻少妇偷人精品九色| 国产女主播在线喷水免费视频网站| 亚洲精品成人av观看孕妇| 欧美精品人与动牲交sv欧美| 久久午夜福利片| 免费观看在线日韩| 精品久久久久久久末码| 国产探花极品一区二区|