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

    Formation of air-entraining vortices at horizontal intakes without approach flow induced circulation*

    2016-10-14 09:18:36MustafaGOGUSMeteKOKENAliBAYKARA

    Mustafa GOGUS, Mete KOKEN, Ali BAYKARA

    ?

    Formation of air-entraining vortices at horizontal intakes without approach flow induced circulation*

    Mustafa GOGUS, Mete KOKEN, Ali BAYKARA

    Civil Engineering Department, Middle East Technical University, Ankara,Turkey, E-mail: mgogus@metu.edu.tr(Received December 10, 2014, Revised October 5, 2015)

    The aim of this experimental study is to investigate the effects of hydraulic parameters on the formation of air-entraining vortices at horizontal intake structures without approach flow induced circulation. Six intake pipes of different diameters were tested in the study. The intake pipe to be tested was horizontally mounted to the front side of a large reservoir and then for a wide range of discharges experiments were conducted and critical submergences were detected with adjustable approach channel sidewalls. Empiri- cal equations were derived for the dimensionless critical submergence as a function of the relevant dimensionless parameters. Availa- ble data is also checked for the possible scale effect. Then, these obtained equations were compared with the similar ones in the literature which showed a quite good agreement.

    air-entraining vortex, horizontal intakes, model scale effect, sidewall clearance, vortex

    Introduction

    Vortex formation is a phenomenon occurs on the free surface in front of the intakes, it has great impo- rtance in the design criteria of the intakes. The loca- tion and direction of the intake should be so arranged that the water level should be well above the intake to prevent vortex formation even if the reservoir is at dead or minimum storage level. Alden Research La- boratory, ARL, defines six types of vortex formation[1]. The most critical one among these six vortex types is Type 6 where vortex formation causes air ingestion continuously. The vertical distance between the intake and the free surface is called submergence and the submergence where air-entraining vortex just occurs is defined as the “critical submergence”. Depending on the direction of the intake, the critical submergence can be measured from the center or summit point of the intake. For an effective intake, the submergence should be large enough to prevent occurrence of air- entraining vortices.

    The main problems of air-entraining vortex for- mations are increase of head loss, reduction of intake discharge, reduction of efficiency of hydraulic machi- nes due to the lower discharge, operational problems at hydraulic machines, vibrations and cavitation[2]. Air-entraining can be removed by increasing the dista- nce of streamlines between the intake structure and free water surface, improving approach flow conditio- ns and applying some anti-vortex devices[1]. Borghei and Kabiri-Samani[3]and Kabiri-Samani and Borghei[4]investigated the effect of anti-vortex plates on forma- tion of air entraining vortex and critical submergence. Sarkarhed et al.[5]investigated the effect of intake headwall and trash rack on the type and strength of the intake vortices.

    In literature, many researchers have worked to figure out and solve formation of vortices at intake st- ructures by comprehensive analytical and numerical methods. Many models and existing installations were examined by studies conducted at different intake ar- rangements, boundary and approach flow conditions. Researchers like Gordon[6], Knauss[1], Jiming et al.[7], Li et al.[8], Ahmad et al.[9]and Gürbüzdal[10]proposed empirical equations to find critical submergence. Ta?tan and Yildirim[11,12]and Yildirim et al.[13]studied flow conditions and formation of air entraining vorti- ces at vertical intakes. In addition, some of them men- tioned limits of the hydraulic parameters above which vortex formation is independent of these parameters.

    1. Theoretical analysis

    Since vortex formation is the consequence of the complex interaction between the geometry of the in- take zone and approach channel, the flow velocity and the liquid properties, parameters influencing vortex formation at a horizontal intake shown in Fig.1 may be expressed by

    Fig.1 Basic illustration of the experimental intake

    Based on dimensional analysis, critical submer- gence ratio,, was found to be the function of the following dimensionless parameters

    In the experimental setup, the bottom clearance,, was so arranged that it is equal to zero all the time. Therefore,can be dropped from Eq.(2). In addi- tion, experiments were conducted at symmetrical flow conditions and consequently,andare equal. That is why these two parameters were combined and wri- tten asin the reduced form of Eq.(2). Circula- tion depends on the characteristics of the approach flow, the geometry of the intake chamber and the dis- charge. In the present study, since there is no imposed circulation and the flow and geometrical parameters of the intake are already in Eq.(1), the circulation para- meter,, could be removed from Eq.(2). In literatu- re, Weber number is accepted as it is effective in weak dimple like vortices and some of the researchers like Anwar et al.[14]stated that the effect of surface tension and viscous forces can be removed from the air-entrai- ning vortex equation ifandare greater than 104and 3×104, respectively. However, in this study, both effects of Weber and Reynolds number were investigated because these limiting values cover a wide range and they are valid only for their own studies. Froude number is the primary parameter affe- cting vortex formation because air-entraining vortex formation is a free surface phenomenon which is affe- cted by gravity. Also, model studies are generally undertaken based on Froude similitude. By conside- ring aforementioned explanations, the final form of the critical submergence ratio,, can be written in the following form

    2. Experimental equipment and procedure

    An experimental concrete rectangular reservoir was built to analyze the effects of hydraulic paramete- rs on the formation of air-entraining vortices at hori- zontal intakes[15]. The reservoir, shown in Fig.2, is 3.1 m in length and width and 2.2 m in depth and con- sisting of a dead-volume section, an active reservoir section and an intake section. A diaphragm slab was built to create a dead-volume section and active rese- rvoir section by dividing the reservoir into two parts in the horizontal plane and a 0.9 m space was left be- tween the slab and the rear wall of the reservoir so that water could reach the active reservoir part once the dead volume is filled. Screens were placed over the end of the slab and fastened to the reservoir side walls to provide energy dissipation of water and uniform motion of water without any waves and circulation in the reservoir. A part of front of the dead-end above the diaphragm slab, 2.2 m in length, 0.3 m in width and 1.2 m in depth, was extended from the reservoir front plane and was made of plexiglass panels to obtain a good observation behind the intake entrance. Two ple- xiglass side walls were placed orthogonally between the extended dead-end and screens and their location could be altered so that the effect of the side wall clea- rance on the formation of air-entraining vortices could be observed. Plexiglass intake pipes of diameters 0.3 m, 0.25 m, 0.194 m, 0.144 m, 0.1 m and 0.05 m were installed to the extended dead-end for each set of experiments in this order. After mounting the pipe of the desired diameter on this hole, a plexiglass plate was placed to the bottom of the reservoir so that the bottom clearance of the pipe would be zero. The in- take pipes were connected to a 10.0 kW centrifugal pump which conveys water to a steel pipe of 0.165 m diameter. Water fall freely from the steel pipe to an open channel which ends with a rectangular sharp-cre- sted weir. The calibration of the sharp-crested weir was done by using an acoustic discharge measurement device before the experiments were started. The requi- red flow is provided from a large elevated-constant head tank via an inflow pipe of 0.300 m diameter. In addition, a drainage pipe, to equalize the inflow to outflow, is connected to the front dead-end section.

    Fig.2 Plan (top) view and side view (bottom) of the physical model

    After installation of each intake pipe, the plexi- glass side walls were adjusted at required distances into the reservoir and the plexiglass plane at the bo- ttom was levelled accordingly, the experimental setup became ready for the experiment. In the experiments six different pipe diameters were used and for each pipe, six different side wall clearances (0.4 m, 0.6 m, 0.8 m, 1.0 m, 1.2 m and 1.4 m) were tested for a series of discharges. In each experiment, the discha- rge and water surface level were measured when an air-entraining vortex initiated.

    To start the experiment, the inflow pipe was ope- ned and the reservoir was filled to a level which is above the critical submergence. Then, by turning on the pump and opening the valve on the pumping line, water was conveyed from the reservoir to the open channel. The openings of the valve on the pumping line, the inflow pipe valve and the drainage pipe valve were so arranged that the inflow was almost the same as the outflow. To make the observations under steady conditions, the water level was kept constant in the re- servoir for a considerable time, 5 min-10 min. After observing that there was no vortex formation in front of the inlet, by opening the drainage pipe valve at a very small amount, the water level was decreased gra- dually and then kept constant again for a considerable time for the second observation session of the free sur- face. This process was repeated several times until ob- serving the formation of swirls, surface depressions and finally air-entraining vortices on the free surface (This process took 30 min to 2 h for different experi- ments). At this point, the water level and the location of the formation of the air-entraining vortex were noted for that discharge. By changing the opening of the valve on the pumping line and the valve of the in- flow pipe, various values of discharges and the dimen- sionless numbers were obtained for each side wall opening and pipe diameter. This procedure was sta- rted with the maximum discharge provided by the pump with fully open valve on the discharge pipe and continued with reduced discharges by closing the valve at small rates until getting the minimum discha- rge which resulted in an air-entraining vortex. How- ever, for the plexiglass intake pipes with diameters of 0.10 m and 0.05 m, the tests were conducted at lower discharges in order to avoid cavitation.

    3. Results and discussions

    The summary of the present study in terms of the related important dimensionless parameters affecting the critical submergence are presented in Table 1[15,16]. The variation ofwith the related dimensionless parameters,,andare given in Figs.3-5 as a functionfor each pipe tested. From these figures the following assessments can be made.

    Table 1 Summary of the important parameters related to the critical submergence

    Fig.3 Variations of with for various values for pipes

    Fig.4 Variations of for various values for pipes

    Fig.5 Variations of for various values for pipes

    For the remaining values oftested the va- riation ofwithas a function of,andis not very significant and becomes almost negligible as the related parameter,,or, increases. The data offor these intermediate values ofcoincide with each other at the ma- ximum values of,andtested. These ma- ximum values of,,and corresponding li- miting values ofabove whichis in-dependent of, as a function of intake diameter ofwere determined from the data of Figs.3-5 and listed in Table 2. From Table 2, Fig.6 which shows the limit values ofwas prepared and plotted. Ac- cording to Fig.6,is independent of the side wall clearance as long asvalues are larger than those stated on the figure for a given intake pipe dia- meter, and is only function of,andifis larger than those corresponding values given on the figure (the limit values of other parameters,and, are not given on the figure). The area below the curve given in Fig.6 belongs to a zone whereis a function ofand, while in the area above the curve whereis larger than that stated on the curve,is independent ofand fun- ction offor an intake of known. This figure also reveals another important information for the in- take pipes of larger diameters not tested in this study. From the general trend of the curve, it can be conclu- ded that as the intake pipe diameter gets larger than 0.3 m, the rate of change ofwithgets smaller and maybe approaches to zero at much higher values of. At the same time, at largervalues, the limit values forand, above whichwill be independent of, get smaller than those of the intake pipe of. Table 3 showsvalues of the extreme cases which yield the maximum and minimumdata, as well as intermediatedata as a function ofandfor the intake pipes tested.

    Table 2 Limit values of above which is in- dependent of as a function of , and for the intake pipe of diameter of

    Table 2 Limit values of above which is in- dependent of as a function of , and for the intake pipe of diameter of

    /m 0.05012.0010.68373 00038 400 0.10010.006.45635 00056 000 0.1446.943.20545 00028 500 0.1945.161.53409 00011 900 0.2504.000.73284 0004 460 0.3003.330.50258 0003 000

    Fig.6 Limit values of above which is indepe- ndent of and function of , and

    Table 3 values of the intake pipes of resulted in maximum, minimum and intermediate

    Table 3 values of the intake pipes of resulted in maximum, minimum and intermediate

    /m/m 0.0500.1000.1440.1940.2500.300 0.4842.782.061.61.33 0.61264.173.092.42.00 0.81685.564.123.22.67 1.0-106.945.164.03.33 1.2-128.336.194.84.00 1.4--9.727.135.64.67 values result in intermediate values result in minimum values result in maximum

    3.1

    Based on the information of Table 3, these three groups are treated separately and the empirical equa- tions for dimensionless critical submergence were provided. In order to determine the equation of the best fit curve to the data ofas a function of,,andconsidering Eq.(3), various types of empirical equations were tested by using regression analysis. Finally, the type of equation given in Eq.(4) was found to be the one having the highest correlation coefficient.

    The derived equations for the maximum, mini- mum and intermediatevalues are as follows

    Equation (5) was derived from thedata of smallvalues tested in the study and results in largevalues which are not desired in practical applications. Equation (7) gives intermediatevalues and corresponds to largevalues. In this zone,values become almost independent ofafter certain limit values ofas a function of intake diameter used. Therefore, this equation can be preferred in the design of intakes. As for Eq.(6), it yields the smallestvalues but theirva- lues are not as large as those of the intermediate zone and even at large,is function of.

    In order to see the effect of,andon the value of, these terms were excluded in a row from Eq.(7) and the regression analysis was app- lied to the available intermediate data and the follo- wing equations were derived.

    From the correlation coefficients and error per- centages of the equations given above, it can be con- cluded that as the related dimensionless terms,and, are dropped from the initial general equation, the error in the calculation ofwill in- crease slightly.

    Finally, if only the limiting data values given on the curve of Fig.6 are used,values with thestated on the curve with the correspondingandfor an intake of known diameter whereis independent of, in the regression analysis, the following equations are derived:

    Figures 7 and 8 show the variation of the measu- red and calculatedvalues which are obtained from Eqs.(11) and (12), respectively. As the paramete- rsandare excluded from the general equation of, Eq.(11), the correlation coefficient of the new equation, Eq.(12), decreases slightly.

    Fig.7 Comparison of measured and calculated values from the application of Eq.(11) to the available interme- diate data

    Fig.8 Comparison of measured and calculated values from the application of Eq.(12) to the available interme- diate data

    3.2

    In hydraulic model studies the aim is to provide a complete similarity between the prototype and the model. In this case one can get correct prototype data directly from the measured or calculated model data. The complete similarity between the model and proto- type is satisfied if all the corresponding dimensionless parameters of the governing equation of the phenome- non are the same. In this study the governing equation is Eq.(3).

    In practice it is not possible to satisfy this condi- tion because it results in model length scale of unity. Therefore in the study of small scale hydraulic models only the equality of the most important dimensionless flow property is considered. This predominant para- meter is Froude number in vortex models and it is as- sumed that the effect of the other parameters conside- ring Eq.(3), the Reynolds and Weber numbers, on the dependent dimensionless parameter,, is negli- gible. Due to this assumption it is obvious that the complete similarity between the two systems cannot be satisfied. Therefore, the measured values offrom the model cannot be directly converted into the prototype value without considering the scale effect of the model. The scale effect of a model decreases as the model length ratio decreases.

    In order to find out the scale effects on the values of, the available experimental results were ana- lyzed and the pipe and flow parameters for the same Froude numbers were determined. Among these data, the ones having almost the samevalues were selected. Each pair of data for a known Froude num- ber belongs to two pipes of different diameters. The first pipe which is larger in diameter than the second pipe is considered as the prototype of the second one. In the observation of the scale effect on the dimen- sionless critical submergence, the following paramete- rs were used: model length scale, the ratio of the dimensionless critical submergence values for the model and the prototype,, the ratios of the Reynolds numbers,, and the ratios of the Weber numbers,.

    Fig.9 Variation of with Froude number

    Fig.10 Variation of with

    Fig.11 Variation of with

    Fig.12 Variation of with

    From the following argument another estimation for the calculation ofcan be made. If it is assumed that Eq.(11) is applicable to prototypes as well as models, one can write the following relation- ship for.

    If Eq.(14) is plotted in Fig.10 for the zone ofwhere the test data are available, it is seen that the curve of Eq.(14) lies between the envelope curves. For a rough calculation of, Eq.(14) may be used within the zone oftested in this study. Additional experiments must be conducted at smallervalues to show that Eq.(14) is applicable or not.

    In a similar way to Fig.10, Figs.11 and 12 show the variation ofwithandwith, respectively. The general trends of the data values of these figures also imply that asandincrease,increase and approach to 1.0 forand.

    Finally, from the analysis of all these figures, it can be concluded that as the model scale,, gets smaller, the scale effect on the value ofbe- comes more significant.

    3.3Validation of the proposed Emprical equations with the related relationships presented in Gürbüzdal’s study

    Gürbüzdal[10]conducted similar experiments to those presented in this study with four pipes of diffe- rent diameters. The experiments were conducted in a similar setup with constantfor a wide range of discharges.values used in the experi- ments and the maximum Froude numbers tested with each intake pipe used are shown in Fig.13 along with the corresponding values of this study which were presented in Fig.6. From Fig.13, it is clearly seen that the data of Gürbüzdal[10]are below the curve of the present study which was named as “intermediate zone” whereis function of,,and. In order to show how well the empirical equation, Eq.(7), proposed in this study fits the data of Gürbüzdal[10], Fig.14 is plotted. In this figure, the cal- culatedvalues were determined using the re- quired values of,,andfrom the data of Gürbüzdal[10]and plotted with respect to those observed by Gürbüzdal[10]as measured. Figure 14 reveals that there is a good agreement between the proposed equation, Eq.(7), and Gürbüzdal’s[10]data which lie withinerror lines.

    Fig.13 For present and Gürbüzdal[10] study, limit values of above which is independent of and function of , and

    Fig.14 Comparison of measured and calculated (from Eq.(7)) intermediate values of Gürbüzdal[10]

    3.4Gordon’s, Reddy and Pickford’s and Gulliver and Rindels’ studies

    Gulliver and Rindel[18]compiled all the available data, together with Gordon[6]’s data, on existing insta- llations and model studies of proposed installations and presented it in a figure, Fig.15, to obtain the most accurate values of submergence as a function of. Two envelope curves which result from Gordon’s[6]criteria are also given in this figure. From this figure, it is seen that neither Gordon’s[6]nor Reddy and Pickford’s[17]design criteria are sufficient to avoid vortex problems. However, there is a region where free-surface vortices are less likely to occur, that is the region segmented by a dimensionless submergence,, greater than 0.7 and an intake Froude number,, less than 0.5. However, if a given intake installation has extremely poor approach flow condi- tions, vortices are still possible in this “safe” region. In order to compare the validity of Eq.(12) with the available data ofin literature, Eq.(12) was plotted into this and presented as Fig.15. The curve of Eq.(12) lies below all the other curves given in the figure. The agreement between Eq.(12) and the rela- tions of Gordon[6]is quite good for the Froude number up to about 0.6. For larger Froude numbers, Eq.(12) estimates much lowervalues than the others. Large Froude numbers were achieved in the present study by using mostly pipes of small diameters. Whereas, most of the data presented in Fig.15 were provided from the existing installations. Therefore, the data of the small diameter pipes used in the present study represent only the data of small scale models, from which one cannot directly convert the model data to the prototype values. The scale effects due to neglecting the similarity of,and probably the wall clearance, are the reasons for not having a good agreement between the model data and prototype data.

    Fig.15 Dimensionless plot of the data obtained from the existi- ng intakes, field installations and model studies[18] with the relationships proposed by Gordon[6], Reddy and Ve Pickford[17], Eq.(12) and corrected Eq.(12) for 1/10

    As a rough estimation, if it is assumed that the experimental setups tested in this study were models of the prototypes, with a length scale of, from Eq.(14), the value ofis calculated as 0.5, from which. This means that, if the validity of Eq.(14) is assumed for a model of, the critical submergence of the proto- type may be selected as twice the corresponding model value. In the light of this informationvalues of Eq.(12) can be multiplied by 2 to determine thevalues to be expected in the prototype which are shown in Fig.15 as “corrected Eq.(12)” for this specific example. As it is seen from the figure, the curve of the “corrected Eq.(12)” lies above most of the data having vortex problems.

    4. Conclusions

    In the present study the effects of the Froude, Reynolds, Weber numbers and side wall clearance on the formation of air-entraining vortices at horizontal intakes were investigated experimentally. Empirical equations for dimensionless critical submergence,, were derived from the experimental data as a fun- ction of the relevant dimensionless parameters and the validity of them was checked with the other data in the literature. Considering the model-prototype con- cept, the available data were analyzed and the scale effect on the values of dimensionless critical submer- gence was investigated. From this experimental study, the following conclusions can be drawn:

    (4) As a function of intake pipe diameter, there are limit values for,,andabove whichis independent ofand function of,and.

    (5) For all three cases described above, yielding maximum, minimum and intermediatevalues as a function of, and for only limit values of,and, several empirical equations were derived forwith quite high correlation coeffi- cients.

    (6) Due to the significant scale effects,values to be obtained from model studies should be multiplied by a correction coefficient before calcula- ting the correspondingvalues of the prototype directly from the length scale of the model. As the model length scale,, decreases, the value of the co- rrection coefficient increases.

    Acknowledgements

    This work was supported by the TUBITAK (The Scientific and Technological Research Council of Turkey) (Grant No. 110M676), which is gratefully ac- knowledged here.

    References

    [1] KNAUSS J. Swirling flow problems at intakes (Hy- draulic structures design manual)[M]. Rotterdam, The Netherlands: A. A. Balkema, 1987.

    [2] CHANG Yu-chi, HU Chin-ning and TU Jing-chin et al. Experimental investigation and numerical prediction of cavitation incurred on propeller surfaces[J]. Journal of Hydrodynamics, 2010, 22(5Suppl.): 764-769.

    [3] BORGHEI S. M., KABIRI-SAMANI A. R. Effects of anti-vortex plates on critical submergence at a vertical in- take[J]. Scientia Iranica, 2010, 17(2): 89-95.

    [4] KABIRI-SAMANI A. R., BORGHEI S. M. Effects of anti vortex plates on air entrainment by free vortex[J]. Scientia Iranica, 2013, 20(2): 251-258

    [5] SARKARDEH H., ZARRATI A. R. and ROSHAN R. Ef- fect of intake head wall and trash rack on vortices[J]. Journal of Hydraulic Research, 2010, 48(1): 108-112.

    [6] GORDON J. L. Vortices at intakes[J]. Water Power, 1970, 22(4): 137-138.

    [7] JIMING M., YUANBO L. and and JITANG H. Minimum submergence before double-entrance pressure intakes[J]. Journal of Hydraulic Engineering, ASCE, 2000, 126(8): 628-631.

    [8] LI Hai-feng, CHEN Hong-xun and MA Zheng et al. Expe- rimental and numerical investigation of free surface vor- tex[J]. Journal of Hydrodynamics, 2008, 20(4): 485-491.

    [9] AHMAD Z., RAO K. V. and and MITTAL M. K. Critical submergence for horizontal intakes in open channel flows[J]. Journal of Dam Engineering, 2008, 19(3): 169- 184.

    [10] GüRBüZDAL F. Scale effects on the formation of vor- tices at intake structures[D]. Master Thesis,, Ankara, Turkey: Middle East Technical University, 2009.

    [11] TA?TAN K., YILDIRIM N. Effects of dimensionless pa- rameters on an air-entrainig vortices[J]. Journal of Hy- draulic Research, 2010, 48(1): 57-64.

    [12] TA?TAN K., YILDIRIM N. Effects of Froude, Reynolds, and Weber numbers on an air-entrainig vortex[J]. Journal of Hydraulic Research, 2014, 52(3): 421-425.

    [13] YILDIRIM N., EYüPO?LU A. S. and TA?TAN K. Criti- cal submergence for dual rectangular intakes[J]. Journal of Energy Engineering, 2012, 138(4): 237-245.

    [14] ANWAR H. O., WELLER J. A. and AMPHLETT M. B. Similarity of free-vortex at horizontal intake[J]. Journal of Hydraulic Research, 1978, 16(2): 95-105.

    [15] BAYKARA A. Effect of hydraulic parameters on the formation of vortices at intake structures[D]. Master Thesis, Ankara, Turkey: Middle East Technical University (METU), 2013.

    [16] GOGUS M., KOKEN M. Formation of vortices at water intake structures and determination of devices to eliminate them[R]. Turkish Science Foundation Report No:110M676, 2013.

    [17] REDDY Y. R. VE PICKFORD J. A. Vortices at intakes in conventional sumps[J]. Water Power, 1972, 24(3): 108- 109.

    [18] GULLIVER J. S., RINDELS A. J. An experimental study of critical submergence to avoid free-surface vortices at vertical intakes[R]. St. Anthony Falls Laboratory, Univer- sity of Minnesota, Project No. 224, 1983.


    * Biography: Mustafa GOGUS (1952-), Male, Ph. D.,Professor

    Corresponging author: Mete KOKEN, E-mail: mkoken@metu.edu.tr

    10.1016/S1001-6058(16)60612-1 2016,28(1):102-113

    久久久国产成人免费| 欧美一区二区精品小视频在线| 人人妻人人澡人人爽人人夜夜 | 国产精品嫩草影院av在线观看| 亚洲成a人片在线一区二区| 国产日韩欧美在线精品| 成人特级黄色片久久久久久久| 美女大奶头视频| 国国产精品蜜臀av免费| 精品欧美国产一区二区三| 蜜桃久久精品国产亚洲av| 中国美白少妇内射xxxbb| 九九在线视频观看精品| 国产高清激情床上av| 国产精品三级大全| 熟女人妻精品中文字幕| 十八禁国产超污无遮挡网站| 插逼视频在线观看| 久久久久久伊人网av| 成人午夜高清在线视频| 国产精品久久久久久精品电影小说 | 深夜a级毛片| 最好的美女福利视频网| 激情 狠狠 欧美| 精品人妻熟女av久视频| 久久这里有精品视频免费| 在线观看66精品国产| 日韩欧美精品v在线| 一边亲一边摸免费视频| 毛片女人毛片| 久久人妻av系列| 国产毛片a区久久久久| 2021天堂中文幕一二区在线观| 国产国拍精品亚洲av在线观看| 日本爱情动作片www.在线观看| 51国产日韩欧美| 真实男女啪啪啪动态图| 欧美日韩综合久久久久久| 亚洲精品日韩av片在线观看| 青春草亚洲视频在线观看| 亚洲国产色片| 欧美潮喷喷水| 又粗又爽又猛毛片免费看| 亚洲熟妇中文字幕五十中出| 久久久国产成人免费| 男女那种视频在线观看| 寂寞人妻少妇视频99o| 日日撸夜夜添| 老司机影院成人| 国产 一区 欧美 日韩| 看十八女毛片水多多多| 久久亚洲精品不卡| 欧美在线一区亚洲| 色尼玛亚洲综合影院| 亚洲在久久综合| 久久热精品热| 国产av不卡久久| 男女视频在线观看网站免费| 中文字幕免费在线视频6| 波多野结衣巨乳人妻| 身体一侧抽搐| 九色成人免费人妻av| 午夜精品在线福利| 国产欧美日韩精品一区二区| 久久久午夜欧美精品| 欧美一区二区国产精品久久精品| 日韩一区二区三区影片| 日日摸夜夜添夜夜添av毛片| 美女大奶头视频| 少妇裸体淫交视频免费看高清| 人妻系列 视频| 成人午夜精彩视频在线观看| 男女啪啪激烈高潮av片| 欧美性猛交╳xxx乱大交人| 能在线免费观看的黄片| 波多野结衣巨乳人妻| 午夜福利高清视频| 美女xxoo啪啪120秒动态图| 午夜激情欧美在线| 日韩欧美在线乱码| 久久人人精品亚洲av| 嫩草影院精品99| 国产真实乱freesex| 亚洲在线观看片| 精品熟女少妇av免费看| 国产一区二区亚洲精品在线观看| 亚洲国产精品成人久久小说 | 美女xxoo啪啪120秒动态图| 久久久久性生活片| 91狼人影院| 亚洲成人精品中文字幕电影| 校园人妻丝袜中文字幕| 日本在线视频免费播放| av.在线天堂| 亚洲一级一片aⅴ在线观看| 性色avwww在线观看| 亚洲天堂国产精品一区在线| 久久精品91蜜桃| 国产一区亚洲一区在线观看| 国产精品一区www在线观看| 可以在线观看毛片的网站| 两个人的视频大全免费| 国产精品爽爽va在线观看网站| 国产麻豆成人av免费视频| 日韩大尺度精品在线看网址| 91午夜精品亚洲一区二区三区| 欧美一区二区国产精品久久精品| 久久久国产成人精品二区| 欧美+日韩+精品| 日韩欧美精品v在线| 最近中文字幕高清免费大全6| 国产一级毛片在线| av福利片在线观看| 亚洲av免费在线观看| 一区二区三区免费毛片| 免费观看精品视频网站| 中文字幕av成人在线电影| 综合色av麻豆| 国产精品伦人一区二区| 欧美区成人在线视频| 美女国产视频在线观看| 能在线免费观看的黄片| 亚洲最大成人中文| 丰满人妻一区二区三区视频av| 亚洲精品国产av成人精品| 亚洲不卡免费看| 精品不卡国产一区二区三区| 国产成人福利小说| 一级av片app| 岛国毛片在线播放| 日韩三级伦理在线观看| 久久人人爽人人片av| 人人妻人人看人人澡| 免费看a级黄色片| 日本免费a在线| 人妻久久中文字幕网| 成年版毛片免费区| АⅤ资源中文在线天堂| 蜜桃亚洲精品一区二区三区| 精品不卡国产一区二区三区| 久久久久久国产a免费观看| 韩国av在线不卡| 亚洲人成网站在线观看播放| 精华霜和精华液先用哪个| 丝袜美腿在线中文| 日韩大尺度精品在线看网址| 久久国产乱子免费精品| 春色校园在线视频观看| 久久久久免费精品人妻一区二区| 99热网站在线观看| 少妇人妻精品综合一区二区 | 中文精品一卡2卡3卡4更新| 国产真实乱freesex| 欧美日本视频| 在线观看一区二区三区| a级毛片a级免费在线| 男人狂女人下面高潮的视频| 夜夜看夜夜爽夜夜摸| 有码 亚洲区| 国产视频内射| 国产精品.久久久| 亚洲av.av天堂| 精品一区二区三区视频在线| 91狼人影院| 精品国产三级普通话版| 日韩欧美精品免费久久| a级一级毛片免费在线观看| 神马国产精品三级电影在线观看| 久久久精品欧美日韩精品| 特大巨黑吊av在线直播| 国产一级毛片七仙女欲春2| 国产一区二区激情短视频| 99热这里只有是精品在线观看| 少妇熟女aⅴ在线视频| 最新中文字幕久久久久| 精品欧美国产一区二区三| 国产精品人妻久久久久久| 美女大奶头视频| 午夜福利视频1000在线观看| 99久国产av精品| 国产成人freesex在线| 在线观看免费视频日本深夜| 久久99精品国语久久久| 国内揄拍国产精品人妻在线| 亚洲精品国产成人久久av| 国产黄色视频一区二区在线观看 | 看黄色毛片网站| 两性午夜刺激爽爽歪歪视频在线观看| 三级国产精品欧美在线观看| 美女脱内裤让男人舔精品视频 | 精品久久国产蜜桃| 国产一区二区在线av高清观看| 69人妻影院| 国产精品伦人一区二区| 久久精品国产亚洲av涩爱 | 日本欧美国产在线视频| 特大巨黑吊av在线直播| 中国美白少妇内射xxxbb| 国产女主播在线喷水免费视频网站 | 成人毛片60女人毛片免费| 国产三级中文精品| 自拍偷自拍亚洲精品老妇| 欧美精品国产亚洲| 性色avwww在线观看| 亚洲久久久久久中文字幕| 少妇熟女aⅴ在线视频| 欧美日本亚洲视频在线播放| 国产高清激情床上av| avwww免费| 久久久精品94久久精品| 欧美在线一区亚洲| 国产精品一二三区在线看| 午夜视频国产福利| 你懂的网址亚洲精品在线观看 | 久久精品夜夜夜夜夜久久蜜豆| 综合色av麻豆| 午夜激情欧美在线| 久久久久久久亚洲中文字幕| 中文资源天堂在线| av天堂中文字幕网| 欧美一区二区精品小视频在线| 亚洲精品日韩av片在线观看| 国产一区亚洲一区在线观看| 免费搜索国产男女视频| 国产淫片久久久久久久久| 99久国产av精品| a级一级毛片免费在线观看| 亚洲成人久久性| 国产三级中文精品| 黑人高潮一二区| 免费观看的影片在线观看| 国产综合懂色| 亚洲人成网站在线播放欧美日韩| 久久中文看片网| 精品久久久久久久末码| 国产成人午夜福利电影在线观看| 久久中文看片网| 秋霞在线观看毛片| 亚洲第一区二区三区不卡| 久久综合国产亚洲精品| 一级av片app| 国产精品乱码一区二三区的特点| 我的女老师完整版在线观看| 国产片特级美女逼逼视频| 免费av观看视频| 色综合亚洲欧美另类图片| 免费在线观看成人毛片| 午夜福利在线观看免费完整高清在 | 国产色爽女视频免费观看| 2021天堂中文幕一二区在线观| 一本久久精品| 亚洲色图av天堂| 成年女人永久免费观看视频| 久久久a久久爽久久v久久| 日韩欧美精品v在线| 一本一本综合久久| 人人妻人人澡人人爽人人夜夜 | 特大巨黑吊av在线直播| 成人亚洲欧美一区二区av| 九九热线精品视视频播放| 久久精品国产鲁丝片午夜精品| 波多野结衣高清无吗| 成人无遮挡网站| 国产精品人妻久久久久久| 国产精品一及| 婷婷色av中文字幕| 亚洲成人久久性| 欧美成人a在线观看| 在线播放无遮挡| 亚洲自偷自拍三级| 亚洲美女搞黄在线观看| 国产 一区精品| 给我免费播放毛片高清在线观看| 免费无遮挡裸体视频| 淫秽高清视频在线观看| 一级毛片aaaaaa免费看小| 99久久精品热视频| 啦啦啦啦在线视频资源| 国产精品久久久久久亚洲av鲁大| 1024手机看黄色片| 婷婷精品国产亚洲av| 美女xxoo啪啪120秒动态图| 国产色婷婷99| 美女黄网站色视频| 伊人久久精品亚洲午夜| 高清日韩中文字幕在线| 成人二区视频| 看免费成人av毛片| 久久精品国产亚洲av香蕉五月| 亚洲高清免费不卡视频| 亚洲图色成人| 老熟妇乱子伦视频在线观看| 久久精品国产清高在天天线| 毛片女人毛片| 色尼玛亚洲综合影院| 国产av一区在线观看免费| 婷婷精品国产亚洲av| 美女黄网站色视频| 男女视频在线观看网站免费| 麻豆乱淫一区二区| 日韩一区二区视频免费看| 精品人妻偷拍中文字幕| 我要看日韩黄色一级片| 99热这里只有是精品在线观看| 欧美潮喷喷水| 国产高潮美女av| 亚洲欧洲日产国产| 久久精品国产亚洲av香蕉五月| 亚洲精品国产成人久久av| 天堂网av新在线| 亚洲欧美日韩高清在线视频| 大型黄色视频在线免费观看| 午夜老司机福利剧场| 精品国内亚洲2022精品成人| 亚洲综合色惰| 成人二区视频| 欧美+亚洲+日韩+国产| 国产亚洲av嫩草精品影院| 欧美变态另类bdsm刘玥| 精华霜和精华液先用哪个| 久久久色成人| 在线a可以看的网站| 欧美日韩乱码在线| av在线亚洲专区| 成年女人永久免费观看视频| 亚洲在线观看片| 在线观看一区二区三区| 国产单亲对白刺激| 亚洲真实伦在线观看| 午夜福利在线观看吧| 少妇高潮的动态图| 我要搜黄色片| 色5月婷婷丁香| 久久精品影院6| 99热这里只有精品一区| 中文欧美无线码| 最近的中文字幕免费完整| 大香蕉久久网| 天堂√8在线中文| 国产91av在线免费观看| 国产极品天堂在线| 热99re8久久精品国产| 欧美激情国产日韩精品一区| 嫩草影院精品99| 国产精品野战在线观看| 国产黄a三级三级三级人| 又黄又爽又刺激的免费视频.| 精品国产三级普通话版| 欧美另类亚洲清纯唯美| 男女下面进入的视频免费午夜| 亚洲成a人片在线一区二区| 国产爱豆传媒在线观看| 久久久久网色| 国产极品精品免费视频能看的| 国模一区二区三区四区视频| 日本成人三级电影网站| 成人无遮挡网站| 网址你懂的国产日韩在线| 嫩草影院新地址| 午夜免费男女啪啪视频观看| 久久精品综合一区二区三区| 日韩一区二区视频免费看| 国产精品国产三级国产av玫瑰| 一级av片app| 日本与韩国留学比较| 51国产日韩欧美| 亚洲第一区二区三区不卡| 噜噜噜噜噜久久久久久91| 精品人妻一区二区三区麻豆| 哪个播放器可以免费观看大片| 一级黄片播放器| 亚洲中文字幕一区二区三区有码在线看| 超碰av人人做人人爽久久| 成人av在线播放网站| 日韩国内少妇激情av| 一级二级三级毛片免费看| 国内精品久久久久精免费| 国产午夜精品久久久久久一区二区三区| 能在线免费观看的黄片| 欧美最黄视频在线播放免费| 国产黄a三级三级三级人| 欧美成人a在线观看| 亚洲成av人片在线播放无| 亚洲国产欧美人成| 级片在线观看| 久久99热这里只有精品18| 一本一本综合久久| 一进一出抽搐gif免费好疼| 两个人视频免费观看高清| 成人永久免费在线观看视频| 蜜臀久久99精品久久宅男| 午夜激情福利司机影院| 国产探花极品一区二区| 亚洲av熟女| 国产精品伦人一区二区| 久久精品国产自在天天线| 人妻久久中文字幕网| 成人毛片a级毛片在线播放| 男插女下体视频免费在线播放| 观看免费一级毛片| 久久久成人免费电影| 99热全是精品| av在线观看视频网站免费| 18禁在线播放成人免费| 国产精品乱码一区二三区的特点| 日韩欧美一区二区三区在线观看| 午夜激情福利司机影院| 啦啦啦韩国在线观看视频| 亚州av有码| 亚洲av不卡在线观看| 国产精品无大码| 国产片特级美女逼逼视频| 精品一区二区三区人妻视频| 国产精品久久久久久av不卡| 久久久久久久久中文| 国产午夜精品久久久久久一区二区三区| 久99久视频精品免费| 在线观看av片永久免费下载| 成人无遮挡网站| 国产极品精品免费视频能看的| 亚洲欧美日韩高清专用| 熟女电影av网| 精品一区二区三区视频在线| 亚洲电影在线观看av| 亚洲av免费高清在线观看| 免费av观看视频| 美女xxoo啪啪120秒动态图| 久久久久久久久中文| 99在线人妻在线中文字幕| 联通29元200g的流量卡| 少妇被粗大猛烈的视频| 最近手机中文字幕大全| 99久久成人亚洲精品观看| 欧美高清成人免费视频www| 久久精品国产亚洲av涩爱 | 亚洲国产精品成人综合色| 国产精品久久久久久精品电影小说 | 91久久精品国产一区二区三区| 亚洲国产精品成人久久小说 | 国产精品国产三级国产av玫瑰| 午夜福利视频1000在线观看| 亚洲成av人片在线播放无| 日韩欧美精品v在线| 哪里可以看免费的av片| 亚洲国产精品成人久久小说 | 精品人妻视频免费看| 亚洲精品久久国产高清桃花| 中文资源天堂在线| 午夜福利在线在线| 亚洲精品国产av成人精品| 国产探花极品一区二区| 午夜福利在线观看吧| 国产成人一区二区在线| 卡戴珊不雅视频在线播放| 免费av不卡在线播放| 特大巨黑吊av在线直播| 国产爱豆传媒在线观看| 亚洲国产色片| 最好的美女福利视频网| 婷婷色av中文字幕| 亚洲乱码一区二区免费版| 午夜精品国产一区二区电影 | 韩国av在线不卡| 国产精品久久久久久久电影| 久久久久久大精品| 十八禁国产超污无遮挡网站| 国产三级中文精品| 99riav亚洲国产免费| 国产色爽女视频免费观看| 永久网站在线| 长腿黑丝高跟| 国产淫片久久久久久久久| 亚洲,欧美,日韩| 97超视频在线观看视频| 久久九九热精品免费| 亚洲av中文字字幕乱码综合| 日韩欧美一区二区三区在线观看| 简卡轻食公司| 亚洲经典国产精华液单| 久久精品国产清高在天天线| 啦啦啦啦在线视频资源| 白带黄色成豆腐渣| 国产国拍精品亚洲av在线观看| 亚洲人成网站在线播| 亚洲欧美精品自产自拍| 日本爱情动作片www.在线观看| 国产高清三级在线| 精品一区二区免费观看| 人妻少妇偷人精品九色| 波野结衣二区三区在线| av又黄又爽大尺度在线免费看 | 亚洲av不卡在线观看| 看免费成人av毛片| 色综合站精品国产| 久久久久久久久久久丰满| 欧美最新免费一区二区三区| 国产男人的电影天堂91| 如何舔出高潮| 人人妻人人澡欧美一区二区| 美女cb高潮喷水在线观看| 欧美人与善性xxx| 黑人高潮一二区| 亚洲图色成人| 在线a可以看的网站| 一级毛片久久久久久久久女| 99精品在免费线老司机午夜| 亚洲精品亚洲一区二区| 久久九九热精品免费| 国产精品.久久久| 一级毛片久久久久久久久女| 日韩一区二区视频免费看| 精品午夜福利在线看| 国产v大片淫在线免费观看| 亚洲自拍偷在线| 97热精品久久久久久| 欧美高清成人免费视频www| 91麻豆精品激情在线观看国产| 夜夜夜夜夜久久久久| 久久久久免费精品人妻一区二区| 精华霜和精华液先用哪个| 99热这里只有是精品在线观看| 国产白丝娇喘喷水9色精品| 欧美性感艳星| 中文字幕人妻熟人妻熟丝袜美| 久久精品国产自在天天线| 精品欧美国产一区二区三| 麻豆成人av视频| 九色成人免费人妻av| 日本撒尿小便嘘嘘汇集6| 黄色一级大片看看| 国产一区二区在线av高清观看| 一本精品99久久精品77| av天堂在线播放| 麻豆成人午夜福利视频| 97超视频在线观看视频| 欧美3d第一页| 日本黄大片高清| 国产麻豆成人av免费视频| 久久久久久九九精品二区国产| 国产亚洲欧美98| 亚洲精品乱码久久久久久按摩| 老司机福利观看| eeuss影院久久| 级片在线观看| 欧美一区二区亚洲| 国产伦一二天堂av在线观看| 午夜福利在线观看吧| 91久久精品电影网| 色播亚洲综合网| 一区二区三区高清视频在线| 久久久国产成人免费| 91aial.com中文字幕在线观看| 美女高潮的动态| 国产高清视频在线观看网站| 午夜精品一区二区三区免费看| 免费观看人在逋| 欧美三级亚洲精品| 欧美色视频一区免费| 国产一区二区三区av在线 | 亚洲精品自拍成人| 亚洲久久久久久中文字幕| 26uuu在线亚洲综合色| 中出人妻视频一区二区| 悠悠久久av| 可以在线观看的亚洲视频| 国产在线男女| 日韩强制内射视频| 国内精品久久久久精免费| 亚洲av成人av| 午夜福利成人在线免费观看| 天天躁夜夜躁狠狠久久av| av视频在线观看入口| 国产精品乱码一区二三区的特点| 色哟哟·www| 波多野结衣巨乳人妻| 国产精品一及| 亚洲人成网站高清观看| 精品久久久久久久久亚洲| 亚洲久久久久久中文字幕| av免费观看日本| 亚洲图色成人| 看非洲黑人一级黄片| 国产又黄又爽又无遮挡在线| 日韩一区二区三区影片| 精品久久久噜噜| 综合色av麻豆| 国产午夜精品一二区理论片| 嫩草影院新地址| 免费看光身美女| av视频在线观看入口| 在线播放无遮挡| 听说在线观看完整版免费高清| 免费人成在线观看视频色| 丰满乱子伦码专区| 久久鲁丝午夜福利片| 我要看日韩黄色一级片| 看黄色毛片网站| 亚洲国产欧美在线一区| 国产激情偷乱视频一区二区| 国产高清视频在线观看网站| 看片在线看免费视频| 精品一区二区免费观看| 午夜精品一区二区三区免费看| 我的女老师完整版在线观看| av在线老鸭窝| 久久久成人免费电影| 青春草视频在线免费观看| av又黄又爽大尺度在线免费看 | 国产私拍福利视频在线观看| 国产精品麻豆人妻色哟哟久久 | 国产一区二区在线观看日韩| 日本撒尿小便嘘嘘汇集6| 久久久成人免费电影| 成人美女网站在线观看视频| 国产成人91sexporn|