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

    欧美xxⅹ黑人| 免费播放大片免费观看视频在线观看| 18禁裸乳无遮挡免费网站照片| 一级毛片电影观看| 国产一区二区三区av在线| 国产av码专区亚洲av| 亚洲精华国产精华液的使用体验| 纯流量卡能插随身wifi吗| 亚洲人与动物交配视频| 中文字幕精品免费在线观看视频 | 少妇的逼水好多| 免费黄网站久久成人精品| 最近手机中文字幕大全| 国产成人a∨麻豆精品| 偷拍熟女少妇极品色| 伦理电影免费视频| av免费观看日本| 欧美xxⅹ黑人| 精品国产露脸久久av麻豆| 黄色配什么色好看| 欧美日本视频| 国产亚洲一区二区精品| 一区二区三区乱码不卡18| 激情五月婷婷亚洲| 99久久精品热视频| 欧美亚洲 丝袜 人妻 在线| 91久久精品电影网| 国产熟女欧美一区二区| 久久精品国产亚洲av涩爱| 在线观看人妻少妇| 日本欧美国产在线视频| 91在线精品国自产拍蜜月| 国产高清国产精品国产三级 | 最近2019中文字幕mv第一页| 国产在视频线精品| videossex国产| av不卡在线播放| videossex国产| 国产免费一区二区三区四区乱码| 久久影院123| 九色成人免费人妻av| 汤姆久久久久久久影院中文字幕| 久久久午夜欧美精品| 永久免费av网站大全| 久久国产精品大桥未久av | 麻豆成人av视频| 亚洲av在线观看美女高潮| 国产一级毛片在线| 91精品伊人久久大香线蕉| 亚洲av福利一区| 日韩精品有码人妻一区| 午夜日本视频在线| 成人亚洲精品一区在线观看 | 国产 一区 欧美 日韩| 九九爱精品视频在线观看| 超碰97精品在线观看| 国内精品宾馆在线| 91久久精品国产一区二区成人| 日韩成人av中文字幕在线观看| 有码 亚洲区| 最近最新中文字幕免费大全7| 欧美bdsm另类| 午夜日本视频在线| 日韩欧美 国产精品| 两个人的视频大全免费| 日本免费在线观看一区| 三级经典国产精品| 亚洲三级黄色毛片| 国产亚洲av片在线观看秒播厂| 久久精品国产a三级三级三级| 日本-黄色视频高清免费观看| 99精国产麻豆久久婷婷| 极品少妇高潮喷水抽搐| 中文字幕免费在线视频6| 最近2019中文字幕mv第一页| 狂野欧美白嫩少妇大欣赏| 久久久色成人| 亚洲美女搞黄在线观看| 精品少妇久久久久久888优播| 亚洲国产精品一区三区| 超碰av人人做人人爽久久| 国产高清有码在线观看视频| 美女福利国产在线 | 国产老妇伦熟女老妇高清| 美女内射精品一级片tv| 水蜜桃什么品种好| 大片电影免费在线观看免费| 青春草国产在线视频| 91精品伊人久久大香线蕉| av国产久精品久网站免费入址| 毛片一级片免费看久久久久| 少妇猛男粗大的猛烈进出视频| 国产成人免费无遮挡视频| 嫩草影院新地址| 久久99热这里只频精品6学生| 18禁裸乳无遮挡免费网站照片| 国产精品熟女久久久久浪| 成人无遮挡网站| 插阴视频在线观看视频| 在线免费十八禁| 国产视频内射| 国产精品久久久久久精品电影小说 | 亚洲av成人精品一二三区| 精品国产三级普通话版| 各种免费的搞黄视频| 丰满乱子伦码专区| 在线观看美女被高潮喷水网站| 国国产精品蜜臀av免费| 国产视频首页在线观看| 精品久久国产蜜桃| 九色成人免费人妻av| 日韩欧美一区视频在线观看 | 18禁在线无遮挡免费观看视频| 天堂8中文在线网| 国产精品免费大片| 欧美精品国产亚洲| 十分钟在线观看高清视频www | 一本色道久久久久久精品综合| www.色视频.com| 九九久久精品国产亚洲av麻豆| 1000部很黄的大片| 最近中文字幕2019免费版| 嫩草影院新地址| 亚洲精品色激情综合| 婷婷色综合大香蕉| 亚洲真实伦在线观看| 成人黄色视频免费在线看| 欧美日韩在线观看h| 免费高清在线观看视频在线观看| 中文天堂在线官网| 国产在线免费精品| 麻豆成人av视频| 午夜视频国产福利| 国产欧美日韩一区二区三区在线 | 色视频www国产| 六月丁香七月| 久久久久久久精品精品| 超碰av人人做人人爽久久| 女人十人毛片免费观看3o分钟| 一级a做视频免费观看| 丰满人妻一区二区三区视频av| 菩萨蛮人人尽说江南好唐韦庄| 十分钟在线观看高清视频www | 亚洲国产成人一精品久久久| 亚洲国产精品成人久久小说| 熟女人妻精品中文字幕| 青春草亚洲视频在线观看| 日韩不卡一区二区三区视频在线| 少妇人妻一区二区三区视频| 在线观看av片永久免费下载| 国产精品福利在线免费观看| 精品酒店卫生间| 久久av网站| 国产探花极品一区二区| 网址你懂的国产日韩在线| 国产成人精品一,二区| 晚上一个人看的免费电影| 一个人免费看片子| 少妇精品久久久久久久| 色视频www国产| 久热久热在线精品观看| 久久久久久久大尺度免费视频| 97在线人人人人妻| 亚洲第一av免费看| 久久婷婷青草| 日韩欧美一区视频在线观看 | 久久99热6这里只有精品| 欧美日韩国产mv在线观看视频 | 99热这里只有精品一区| 亚洲av在线观看美女高潮| 一个人免费看片子| av.在线天堂| 国产欧美亚洲国产| 亚洲久久久国产精品| 尾随美女入室| 伦精品一区二区三区| 亚洲色图av天堂| 身体一侧抽搐| 国产精品国产三级专区第一集| 看非洲黑人一级黄片| 久久人人爽人人爽人人片va| 黑人高潮一二区| 丰满乱子伦码专区| 久久精品久久精品一区二区三区| 六月丁香七月| 91精品国产国语对白视频| 卡戴珊不雅视频在线播放| 毛片一级片免费看久久久久| 久久精品久久精品一区二区三区| www.色视频.com| 国产精品免费大片| 国产高清国产精品国产三级 | 99热全是精品| 美女视频免费永久观看网站| 97在线人人人人妻| 国产精品一二三区在线看| 久久精品久久久久久久性| 丝瓜视频免费看黄片| 国产色婷婷99| 亚洲成色77777| 日韩av不卡免费在线播放| av国产免费在线观看| 国产成人91sexporn| 麻豆国产97在线/欧美| 五月玫瑰六月丁香| 有码 亚洲区| 青春草亚洲视频在线观看| 免费观看无遮挡的男女| 亚洲国产最新在线播放| 91在线精品国自产拍蜜月| 久久6这里有精品| 精品久久久久久久末码| 97热精品久久久久久| 卡戴珊不雅视频在线播放| 卡戴珊不雅视频在线播放| 18禁裸乳无遮挡免费网站照片| av女优亚洲男人天堂| 亚洲国产精品成人久久小说| 国产中年淑女户外野战色| 熟妇人妻不卡中文字幕| 成人综合一区亚洲| 日韩国内少妇激情av| 久久人妻熟女aⅴ| 在线观看免费日韩欧美大片 | 国产一区二区在线观看日韩| 看非洲黑人一级黄片| 黑人高潮一二区| 中文精品一卡2卡3卡4更新| 久久人人爽人人爽人人片va| 国产精品女同一区二区软件| 欧美成人a在线观看| 国产精品熟女久久久久浪| 久久久久性生活片| 中文字幕亚洲精品专区| 国产一区有黄有色的免费视频| 91精品一卡2卡3卡4卡| 黑丝袜美女国产一区| 交换朋友夫妻互换小说| 亚洲成人中文字幕在线播放| 全区人妻精品视频| 高清黄色对白视频在线免费看 | 久久精品熟女亚洲av麻豆精品| 国产伦理片在线播放av一区| 久久女婷五月综合色啪小说| 午夜日本视频在线| 亚洲精品乱码久久久久久按摩| 在线观看一区二区三区| 国产国拍精品亚洲av在线观看| 久久 成人 亚洲| 99热6这里只有精品| 亚洲国产精品国产精品| 大香蕉久久网| 欧美精品一区二区免费开放| 丰满人妻一区二区三区视频av| 久久影院123| 超碰av人人做人人爽久久| 午夜老司机福利剧场| 51国产日韩欧美| 一级爰片在线观看| 成人午夜精彩视频在线观看| 国产午夜精品一二区理论片| 少妇高潮的动态图| 制服丝袜香蕉在线| 最新中文字幕久久久久| 日韩强制内射视频| 在线观看免费高清a一片| 日韩 亚洲 欧美在线| .国产精品久久| 中文天堂在线官网| 一边亲一边摸免费视频| 一级毛片我不卡| 国产一级毛片在线| 色综合色国产| 久久久久精品久久久久真实原创| 色吧在线观看| 国产精品爽爽va在线观看网站| 国模一区二区三区四区视频| 国产精品伦人一区二区| 老师上课跳d突然被开到最大视频| 国产毛片在线视频| 午夜福利在线观看免费完整高清在| 成人国产麻豆网| 色哟哟·www| 18禁裸乳无遮挡免费网站照片| 欧美3d第一页| 亚洲av中文字字幕乱码综合| 成人亚洲精品一区在线观看 | 欧美精品人与动牲交sv欧美| 久久99热这里只有精品18| 91在线精品国自产拍蜜月| 韩国av在线不卡| 老女人水多毛片| 一个人看的www免费观看视频| 国产精品一区二区性色av| 精品国产三级普通话版| 亚洲内射少妇av| 毛片一级片免费看久久久久| 国产成人aa在线观看| av在线app专区| 日本黄色日本黄色录像| 国产精品国产三级专区第一集| 在线观看一区二区三区| 亚洲精品国产色婷婷电影| 国产精品久久久久成人av| 成人二区视频| 免费少妇av软件| 免费看光身美女| 免费高清在线观看视频在线观看| 爱豆传媒免费全集在线观看| 欧美zozozo另类| 少妇熟女欧美另类| 日韩视频在线欧美| 亚洲久久久国产精品| 亚洲欧美一区二区三区黑人 | 中文字幕久久专区| 免费av不卡在线播放| 国产精品无大码| 男女边吃奶边做爰视频| 美女cb高潮喷水在线观看| 狠狠精品人妻久久久久久综合| 中国美白少妇内射xxxbb| 午夜福利高清视频| 国产精品国产三级国产av玫瑰| 免费观看无遮挡的男女| 嘟嘟电影网在线观看| 秋霞伦理黄片| 两个人的视频大全免费| 高清日韩中文字幕在线| 毛片一级片免费看久久久久| 亚洲av综合色区一区| 亚洲av中文字字幕乱码综合| 亚洲第一区二区三区不卡| 日本色播在线视频| 国产精品免费大片| 五月天丁香电影| 久久久久久久亚洲中文字幕| 成人亚洲精品一区在线观看 | 久久热精品热| 亚洲精品色激情综合| 一级毛片我不卡| 免费黄网站久久成人精品| 日日摸夜夜添夜夜添av毛片| 欧美精品人与动牲交sv欧美| 精品久久久久久久久亚洲| 国产精品女同一区二区软件| 国产人妻一区二区三区在| 夜夜爽夜夜爽视频| 九九在线视频观看精品| 精品一区二区免费观看| 国产成人a∨麻豆精品| 亚洲不卡免费看| 伊人久久精品亚洲午夜| 黄片无遮挡物在线观看| 免费久久久久久久精品成人欧美视频 | 王馨瑶露胸无遮挡在线观看| 久久精品夜色国产| 成人一区二区视频在线观看| 高清黄色对白视频在线免费看 | 精华霜和精华液先用哪个| 青青草视频在线视频观看| 日本黄色片子视频| 大香蕉久久网| 99热国产这里只有精品6| 亚洲精品国产av成人精品| 大又大粗又爽又黄少妇毛片口| 你懂的网址亚洲精品在线观看| 日韩精品有码人妻一区| 亚洲精品国产色婷婷电影| 国产精品国产三级国产专区5o| 久久亚洲国产成人精品v| 只有这里有精品99| 舔av片在线| 国产视频首页在线观看| 亚洲欧美中文字幕日韩二区| 国产精品不卡视频一区二区| 亚洲精品自拍成人| 激情 狠狠 欧美| 日韩一区二区三区影片| 80岁老熟妇乱子伦牲交| 久久久久久久久久人人人人人人| 欧美日韩精品成人综合77777| 午夜福利视频精品| 欧美极品一区二区三区四区| 亚洲欧美日韩东京热| 少妇人妻一区二区三区视频| 国产免费视频播放在线视频| 91精品国产九色| 久久久久国产网址| av免费观看日本| 国产成人免费无遮挡视频| 国产精品久久久久久久久免| 不卡视频在线观看欧美| 国产 一区精品| 午夜福利视频精品| 97在线人人人人妻| 全区人妻精品视频| 久久人人爽人人爽人人片va| 精品人妻熟女av久视频| 秋霞伦理黄片| 蜜臀久久99精品久久宅男| 亚洲精品久久午夜乱码| 国产v大片淫在线免费观看| 国产成人a区在线观看| 日本黄大片高清| 有码 亚洲区| 国产黄色视频一区二区在线观看| 国产精品国产av在线观看| 亚洲美女黄色视频免费看| 亚洲婷婷狠狠爱综合网| 国产亚洲av片在线观看秒播厂| 纵有疾风起免费观看全集完整版| 一级a做视频免费观看| 韩国高清视频一区二区三区| 看十八女毛片水多多多| 亚州av有码| 欧美精品亚洲一区二区| 97超碰精品成人国产| 国产精品成人在线| 啦啦啦中文免费视频观看日本| 蜜桃亚洲精品一区二区三区| 国产欧美亚洲国产| 精品一区二区免费观看| 一本一本综合久久| 久久婷婷青草| 国产精品国产三级专区第一集| 欧美一区二区亚洲| 国产精品蜜桃在线观看| 亚洲成人av在线免费| 久久99热这里只有精品18| 久久久久性生活片| 伦理电影大哥的女人| 十分钟在线观看高清视频www | 91aial.com中文字幕在线观看| 国产av国产精品国产| 成人综合一区亚洲| 一区二区三区精品91| 欧美xxxx黑人xx丫x性爽| 精品一品国产午夜福利视频| 精品人妻视频免费看| 久久精品久久精品一区二区三区| 老师上课跳d突然被开到最大视频| 大话2 男鬼变身卡| 亚洲欧洲日产国产| 亚洲精品一区蜜桃| 大话2 男鬼变身卡| 99九九线精品视频在线观看视频| 久久久久久人妻| 人妻夜夜爽99麻豆av| 在线观看美女被高潮喷水网站| 久久精品国产自在天天线| 成人黄色视频免费在线看| 少妇猛男粗大的猛烈进出视频| 尤物成人国产欧美一区二区三区| 久久久久久久久久人人人人人人| 九色成人免费人妻av| 久久久久久久久久久丰满| 边亲边吃奶的免费视频| 色视频在线一区二区三区| 一个人免费看片子| 免费高清在线观看视频在线观看| 国产黄片视频在线免费观看| tube8黄色片| 久久99精品国语久久久| 欧美xxxx黑人xx丫x性爽| 国产深夜福利视频在线观看| 精品人妻偷拍中文字幕| 成人无遮挡网站| 国产在线男女| 一区二区av电影网| 卡戴珊不雅视频在线播放| 国产色爽女视频免费观看| 男人舔奶头视频| av国产久精品久网站免费入址| 中文字幕久久专区| 国产一区二区三区av在线| 一区二区三区乱码不卡18| 大又大粗又爽又黄少妇毛片口| av免费观看日本| 精品久久久久久久末码| 国产欧美日韩精品一区二区| 一级毛片 在线播放| 国产美女午夜福利| 黄色一级大片看看| 午夜激情久久久久久久| 亚洲av综合色区一区| 人妻夜夜爽99麻豆av| 久久精品久久精品一区二区三区| 久久久a久久爽久久v久久| 欧美日韩在线观看h| 国精品久久久久久国模美| 久久这里有精品视频免费| 国产午夜精品久久久久久一区二区三区| 欧美极品一区二区三区四区| 国产av码专区亚洲av| 少妇裸体淫交视频免费看高清| 国产午夜精品一二区理论片| 日韩不卡一区二区三区视频在线| 亚洲欧美一区二区三区国产| 少妇高潮的动态图| 国产一区亚洲一区在线观看| 免费看日本二区| 男女免费视频国产| 久久青草综合色| 中文字幕亚洲精品专区| 在线观看av片永久免费下载| 91在线精品国自产拍蜜月| 中文字幕久久专区| 国产真实伦视频高清在线观看| 一区二区三区乱码不卡18| 久久av网站| 久久久国产一区二区| 熟妇人妻不卡中文字幕| 边亲边吃奶的免费视频| 免费观看无遮挡的男女| 日韩欧美一区视频在线观看 | 亚洲第一区二区三区不卡| 久久影院123| 免费黄频网站在线观看国产| 街头女战士在线观看网站| 高清欧美精品videossex| 国产av码专区亚洲av| 亚洲精品久久午夜乱码| 成人美女网站在线观看视频| 精品久久久久久久久亚洲| 狂野欧美白嫩少妇大欣赏| 久久久久久久久久久免费av| 免费av中文字幕在线| 丝瓜视频免费看黄片| 国产精品偷伦视频观看了| 一区二区三区免费毛片| 午夜视频国产福利| 久久久国产一区二区| 老熟女久久久| 国产精品福利在线免费观看| 国产片特级美女逼逼视频| 五月开心婷婷网| 狂野欧美激情性xxxx在线观看| 观看免费一级毛片| 91aial.com中文字幕在线观看| 精品一区二区三区视频在线| 免费av不卡在线播放| 亚洲精品久久午夜乱码| 亚洲av成人精品一区久久| 日本与韩国留学比较| 亚洲国产色片| 国产精品一区二区在线观看99| 亚洲欧洲日产国产| 免费黄网站久久成人精品| 国产精品av视频在线免费观看| 在线免费观看不下载黄p国产| 肉色欧美久久久久久久蜜桃| 久久 成人 亚洲| 国产伦在线观看视频一区| 成人特级av手机在线观看| 久久久a久久爽久久v久久| 亚洲欧洲国产日韩| 亚洲不卡免费看| 哪个播放器可以免费观看大片| 最近的中文字幕免费完整| 天堂中文最新版在线下载| 王馨瑶露胸无遮挡在线观看| 内射极品少妇av片p| 国产亚洲午夜精品一区二区久久| 一级毛片aaaaaa免费看小| 日日撸夜夜添| 亚洲精品国产av蜜桃| 成人毛片a级毛片在线播放| 国语对白做爰xxxⅹ性视频网站| 大香蕉久久网| 久久97久久精品| av福利片在线观看| 亚洲国产最新在线播放| 国产熟女欧美一区二区| 99久国产av精品国产电影| 两个人的视频大全免费| 国产精品久久久久久精品电影小说 | 亚洲欧美精品自产自拍| 下体分泌物呈黄色| 蜜臀久久99精品久久宅男| 97精品久久久久久久久久精品| 国产精品国产三级专区第一集| 欧美3d第一页| 少妇被粗大猛烈的视频| 老熟女久久久| xxx大片免费视频| 成人午夜精彩视频在线观看| 亚洲av电影在线观看一区二区三区| 亚洲久久久国产精品| 久久韩国三级中文字幕| 国产av码专区亚洲av| 亚洲精品日本国产第一区| 亚洲精品亚洲一区二区| 超碰97精品在线观看| 成人国产麻豆网| 老女人水多毛片| 大码成人一级视频| 中文字幕人妻熟人妻熟丝袜美| 成人国产麻豆网| 亚洲国产成人一精品久久久| av播播在线观看一区| 亚洲内射少妇av| 成人综合一区亚洲| 亚洲欧洲国产日韩| 国产91av在线免费观看| 精品久久久久久久久亚洲| 超碰97精品在线观看| 午夜日本视频在线| 老司机影院成人| 午夜福利在线在线| 人妻少妇偷人精品九色| 午夜福利高清视频| 精品人妻熟女av久视频| 午夜福利高清视频| 啦啦啦啦在线视频资源| 亚洲欧美精品专区久久|