• <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片| 国产有黄有色有爽视频| 边亲边吃奶的免费视频| 亚洲一区二区三区欧美精品| 久久久亚洲精品成人影院| 最新的欧美精品一区二区| 婷婷色麻豆天堂久久| 午夜福利网站1000一区二区三区| 交换朋友夫妻互换小说| 男女边摸边吃奶| 啦啦啦在线观看免费高清www| 美女内射精品一级片tv| 精品亚洲乱码少妇综合久久| 99国产综合亚洲精品| 久久精品国产自在天天线| 久久久久人妻精品一区果冻| 日韩一区二区三区影片| 国产精品国产av在线观看| 满18在线观看网站| 亚洲成国产人片在线观看| 久久午夜福利片| 精品99又大又爽又粗少妇毛片| 伦精品一区二区三区| 日本免费在线观看一区| 欧美精品国产亚洲| 美女国产高潮福利片在线看| 日韩欧美一区视频在线观看| 欧美精品一区二区大全| 一本色道久久久久久精品综合| 大片免费播放器 马上看| 天天躁夜夜躁狠狠躁躁| 波野结衣二区三区在线| 日韩大片免费观看网站| 亚洲四区av| 91成人精品电影| 久久女婷五月综合色啪小说| 精品国产一区二区久久| 在线观看免费高清a一片| 女的被弄到高潮叫床怎么办| 欧美最新免费一区二区三区| 成人国产麻豆网| 高清不卡的av网站| 久久精品国产a三级三级三级| 大香蕉久久网| 久久久久久伊人网av| 亚洲精品久久午夜乱码| 老司机影院成人| 欧美丝袜亚洲另类| 美女国产高潮福利片在线看| 女性被躁到高潮视频| 九色亚洲精品在线播放| 免费看av在线观看网站| 国产亚洲av片在线观看秒播厂| av视频免费观看在线观看| 女人久久www免费人成看片| 午夜福利网站1000一区二区三区| 国产麻豆69| 精品少妇久久久久久888优播| 一级毛片我不卡| 男女边摸边吃奶| 久久国产精品大桥未久av| 男女下面插进去视频免费观看 | 成人国产av品久久久| 男人爽女人下面视频在线观看| 国产精品久久久av美女十八| 日韩av在线免费看完整版不卡| 亚洲欧洲日产国产| 在线观看www视频免费| 久久久久久久久久成人| 一区二区av电影网| 中文字幕精品免费在线观看视频 | 一边亲一边摸免费视频| 欧美精品av麻豆av| 天美传媒精品一区二区| 三级国产精品片| 欧美老熟妇乱子伦牲交| 观看av在线不卡| 亚洲,一卡二卡三卡| 中文字幕制服av| 日韩熟女老妇一区二区性免费视频| 视频区图区小说| 亚洲婷婷狠狠爱综合网| 亚洲精品一区蜜桃| av.在线天堂| 天天影视国产精品| 亚洲国产精品999| 成人无遮挡网站| 久热久热在线精品观看| 人妻一区二区av| 亚洲性久久影院| 国产男女超爽视频在线观看| 午夜日本视频在线| 99久久中文字幕三级久久日本| 亚洲av电影在线观看一区二区三区| 又黄又粗又硬又大视频| 国产在线一区二区三区精| 久久久国产欧美日韩av| 亚洲精品一二三| 国产高清不卡午夜福利| 欧美日韩综合久久久久久| 亚洲性久久影院| 国产免费一级a男人的天堂| 精品人妻偷拍中文字幕| 91午夜精品亚洲一区二区三区| 制服诱惑二区| 天天躁夜夜躁狠狠久久av| 黄片无遮挡物在线观看| xxx大片免费视频| 亚洲成人手机| 欧美3d第一页| 亚洲精品久久成人aⅴ小说| 国产不卡av网站在线观看| 人人妻人人澡人人爽人人夜夜| 最新的欧美精品一区二区| 久久久久视频综合| 国产精品女同一区二区软件| 9热在线视频观看99| 韩国精品一区二区三区 | av国产精品久久久久影院| 内地一区二区视频在线| 亚洲av欧美aⅴ国产| 久久ye,这里只有精品| 国产成人欧美| 欧美最新免费一区二区三区| 一级毛片电影观看| 国产亚洲一区二区精品| 国产av码专区亚洲av| 中文精品一卡2卡3卡4更新| 亚洲精品日本国产第一区| 午夜免费观看性视频| 2022亚洲国产成人精品| 美女福利国产在线| 日本欧美视频一区| 午夜av观看不卡| 久久国产亚洲av麻豆专区| 伊人久久国产一区二区| 午夜日本视频在线| 日韩成人av中文字幕在线观看| 久久久久久人妻| 熟女电影av网| 久久综合国产亚洲精品| 99久久人妻综合| 大话2 男鬼变身卡| 成人18禁高潮啪啪吃奶动态图| 国产欧美日韩一区二区三区在线| 男女下面插进去视频免费观看 | 人人妻人人澡人人爽人人夜夜| 少妇的丰满在线观看| 最近2019中文字幕mv第一页| 久久精品人人爽人人爽视色| 欧美另类一区| 亚洲图色成人| 1024视频免费在线观看| 欧美性感艳星| 欧美成人午夜免费资源| 亚洲av国产av综合av卡| 午夜福利影视在线免费观看| 欧美激情 高清一区二区三区| 欧美人与性动交α欧美精品济南到 | 2022亚洲国产成人精品| av电影中文网址| 考比视频在线观看| 日韩制服骚丝袜av| 日韩一区二区三区影片| 性高湖久久久久久久久免费观看| 国产片内射在线| 18禁动态无遮挡网站| 欧美日韩一区二区视频在线观看视频在线| 91久久精品国产一区二区三区| 女的被弄到高潮叫床怎么办| 一级片免费观看大全| 黄片播放在线免费| 亚洲第一av免费看| 精品国产乱码久久久久久小说| 高清欧美精品videossex| 亚洲人成网站在线观看播放| 国产成人91sexporn| 春色校园在线视频观看| 国产在线视频一区二区| 久久这里只有精品19| 91精品伊人久久大香线蕉| 欧美日韩av久久| 纯流量卡能插随身wifi吗| 巨乳人妻的诱惑在线观看| 99精国产麻豆久久婷婷| 搡女人真爽免费视频火全软件| 91在线精品国自产拍蜜月| 久久精品久久久久久久性| 熟女av电影| 亚洲精品自拍成人| 国产精品免费大片| 免费日韩欧美在线观看| 午夜福利视频在线观看免费| 黑丝袜美女国产一区| 精品一区二区三区四区五区乱码 | 99久久综合免费| 国产深夜福利视频在线观看| 国产成人精品在线电影| 精品第一国产精品| 青春草视频在线免费观看| 香蕉国产在线看| 国产精品久久久久成人av| 日韩成人伦理影院| 免费不卡的大黄色大毛片视频在线观看| 免费黄色在线免费观看| 成年美女黄网站色视频大全免费| 在线观看免费高清a一片| 久久婷婷青草| 日韩不卡一区二区三区视频在线| 韩国高清视频一区二区三区| 人成视频在线观看免费观看| 王馨瑶露胸无遮挡在线观看| 久久久久久久久久成人| 亚洲精品自拍成人| 亚洲四区av| 久久久久精品人妻al黑| www.熟女人妻精品国产 | 九色成人免费人妻av| 欧美亚洲日本最大视频资源| 亚洲精品乱码久久久久久按摩| 美女福利国产在线| 欧美人与善性xxx| 亚洲精品色激情综合| 超碰97精品在线观看| 免费观看av网站的网址| 91精品伊人久久大香线蕉| 两个人免费观看高清视频| 精品一区二区三区四区五区乱码 | 亚洲国产精品一区二区三区在线| 少妇猛男粗大的猛烈进出视频| 久久精品熟女亚洲av麻豆精品| 日本午夜av视频| 国产精品久久久久久精品电影小说| 亚洲精品美女久久av网站| 国产欧美另类精品又又久久亚洲欧美| 女的被弄到高潮叫床怎么办| 母亲3免费完整高清在线观看 | 国产成人a∨麻豆精品| 欧美 日韩 精品 国产| 大陆偷拍与自拍| 久久久精品94久久精品| 日韩视频在线欧美| 美女国产视频在线观看| 男人舔女人的私密视频| 大香蕉久久成人网| 熟女人妻精品中文字幕| 日本免费在线观看一区| 亚洲色图综合在线观看| 亚洲,欧美精品.| 久久久久久伊人网av| 免费看av在线观看网站| 2022亚洲国产成人精品| 丝袜脚勾引网站| 久久国产精品大桥未久av| 久久影院123| 老女人水多毛片| 久久99蜜桃精品久久| 免费日韩欧美在线观看| 岛国毛片在线播放| 久久综合国产亚洲精品| 国产成人午夜福利电影在线观看| 2022亚洲国产成人精品| 女性被躁到高潮视频| 欧美激情极品国产一区二区三区 | 国产一区二区激情短视频 | 欧美老熟妇乱子伦牲交| 汤姆久久久久久久影院中文字幕| 啦啦啦中文免费视频观看日本| 精品亚洲成国产av| 美女大奶头黄色视频| 欧美另类一区| 51国产日韩欧美| 久久精品熟女亚洲av麻豆精品| 巨乳人妻的诱惑在线观看| 亚洲av国产av综合av卡| 美女内射精品一级片tv| tube8黄色片| 欧美精品一区二区大全| 久久人人爽av亚洲精品天堂| 久热久热在线精品观看| 人体艺术视频欧美日本| 国产视频首页在线观看| 97在线人人人人妻| 日韩av在线免费看完整版不卡| 亚洲欧美日韩卡通动漫| 精品人妻一区二区三区麻豆| 天天躁夜夜躁狠狠久久av| 久久99热这里只频精品6学生| 新久久久久国产一级毛片| 90打野战视频偷拍视频| 五月玫瑰六月丁香| 热99国产精品久久久久久7| 涩涩av久久男人的天堂| 在线天堂最新版资源| 女人精品久久久久毛片| 亚洲精品国产av成人精品| 欧美精品高潮呻吟av久久| 狠狠婷婷综合久久久久久88av| 啦啦啦在线观看免费高清www| 97人妻天天添夜夜摸| 人妻一区二区av| 免费在线观看黄色视频的| 咕卡用的链子| 美女中出高潮动态图| 亚洲精品自拍成人| 观看美女的网站| 午夜精品国产一区二区电影| 夫妻性生交免费视频一级片| 五月玫瑰六月丁香| 在线免费观看不下载黄p国产| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 国产精品久久久av美女十八| 国产成人精品无人区| 狂野欧美激情性xxxx在线观看| 国产乱来视频区| 久久狼人影院| 永久网站在线| 亚洲国产av影院在线观看| 午夜福利影视在线免费观看| 青春草视频在线免费观看| 国产黄色免费在线视频| 色视频在线一区二区三区| 国产高清三级在线| 亚洲久久久国产精品| 色婷婷av一区二区三区视频| 曰老女人黄片| 精品人妻熟女毛片av久久网站| 婷婷色综合大香蕉| 在现免费观看毛片| 中国美白少妇内射xxxbb| 成人国产麻豆网| 国产精品一区www在线观看| 伦精品一区二区三区| 十分钟在线观看高清视频www| 国产精品国产三级国产专区5o| 黑丝袜美女国产一区| 男男h啪啪无遮挡| 在线精品无人区一区二区三| 天堂8中文在线网| 久久久精品免费免费高清| 五月玫瑰六月丁香| 夜夜爽夜夜爽视频| 国产激情久久老熟女| 国产色爽女视频免费观看| 极品人妻少妇av视频| 中文字幕亚洲精品专区| 免费播放大片免费观看视频在线观看| 日韩制服丝袜自拍偷拍| www.色视频.com| 国产爽快片一区二区三区| 国产精品一国产av| 久热久热在线精品观看| 又大又黄又爽视频免费| 香蕉精品网在线| 制服诱惑二区| 中文字幕免费在线视频6| 国产日韩一区二区三区精品不卡| 在线观看免费视频网站a站| 亚洲综合色网址| 国产免费又黄又爽又色| 九色成人免费人妻av| 欧美97在线视频| 亚洲av.av天堂| 久久精品久久久久久久性| 午夜免费男女啪啪视频观看| 久久久国产精品麻豆| 国产国拍精品亚洲av在线观看| 51国产日韩欧美| 香蕉丝袜av| 青青草视频在线视频观看| 国产高清三级在线| 亚洲美女黄色视频免费看| 美国免费a级毛片| 看免费成人av毛片| 一级a做视频免费观看| 在线观看一区二区三区激情| 日韩精品有码人妻一区| 亚洲精品成人av观看孕妇| 中文字幕另类日韩欧美亚洲嫩草| 美女xxoo啪啪120秒动态图| 亚洲国产精品国产精品| 天天操日日干夜夜撸| 综合色丁香网| tube8黄色片| 国产一区二区在线观看日韩| 亚洲欧美中文字幕日韩二区| 欧美日本中文国产一区发布| 国产一区亚洲一区在线观看| 18禁动态无遮挡网站| 91午夜精品亚洲一区二区三区| 18禁在线无遮挡免费观看视频| 亚洲精华国产精华液的使用体验| 美女中出高潮动态图| 中文字幕人妻熟女乱码| av女优亚洲男人天堂| 少妇的逼水好多| 各种免费的搞黄视频| 亚洲成人一二三区av| 欧美日韩一区二区视频在线观看视频在线| 久久97久久精品| 日韩精品免费视频一区二区三区 | 欧美少妇被猛烈插入视频| videos熟女内射| 两性夫妻黄色片 | 91精品三级在线观看| 五月玫瑰六月丁香| 久久精品国产a三级三级三级| 少妇被粗大猛烈的视频| 韩国高清视频一区二区三区| 交换朋友夫妻互换小说| 久久国产亚洲av麻豆专区| 国产一区有黄有色的免费视频| 亚洲精品,欧美精品| 9热在线视频观看99| 久久久久久久久久人人人人人人| 搡老乐熟女国产| 成人综合一区亚洲| av福利片在线| 国产精品一二三区在线看| 国产免费福利视频在线观看| 大片免费播放器 马上看| 人人妻人人爽人人添夜夜欢视频| 伦理电影免费视频| 久热这里只有精品99| 免费不卡的大黄色大毛片视频在线观看| 亚洲精品国产av蜜桃| 国产一区二区三区av在线| 国产乱来视频区| 一级片免费观看大全| 婷婷色综合www| 亚洲成色77777| 色视频在线一区二区三区| 中国美白少妇内射xxxbb| 亚洲综合色网址| 久久久久久人妻| 日本与韩国留学比较| 亚洲欧美中文字幕日韩二区| 成年美女黄网站色视频大全免费| 99久久综合免费| 精品少妇黑人巨大在线播放| 国产精品欧美亚洲77777| 亚洲国产毛片av蜜桃av| 国产欧美日韩综合在线一区二区| 少妇精品久久久久久久| 欧美日韩亚洲高清精品| 国产成人精品久久久久久| 成人漫画全彩无遮挡| a级片在线免费高清观看视频| av又黄又爽大尺度在线免费看| 国产深夜福利视频在线观看| 欧美精品一区二区免费开放| 国产精品熟女久久久久浪| 国产高清不卡午夜福利| 国产日韩一区二区三区精品不卡| 成年动漫av网址| 亚洲欧洲国产日韩| 欧美 日韩 精品 国产| 亚洲精品视频女| 视频中文字幕在线观看| 亚洲图色成人| 一级毛片 在线播放| 丝瓜视频免费看黄片| 欧美xxⅹ黑人| 毛片一级片免费看久久久久| 亚洲丝袜综合中文字幕| 少妇的丰满在线观看| 999精品在线视频| 男的添女的下面高潮视频| 99久久精品国产国产毛片| 宅男免费午夜| 午夜免费鲁丝| 欧美97在线视频| 国产精品国产三级国产av玫瑰| 亚洲伊人色综图| a级毛片在线看网站| 一级,二级,三级黄色视频| 最后的刺客免费高清国语| 精品少妇黑人巨大在线播放| 午夜久久久在线观看| 成年女人在线观看亚洲视频| 亚洲国产av新网站| 国产黄色视频一区二区在线观看| 天堂中文最新版在线下载| 国产成人a∨麻豆精品| 岛国毛片在线播放| 亚洲精品久久午夜乱码| 午夜福利视频在线观看免费| 中文字幕人妻熟女乱码| 人妻少妇偷人精品九色| 久久国产精品大桥未久av| 啦啦啦视频在线资源免费观看| 欧美激情国产日韩精品一区| 中文字幕最新亚洲高清| 亚洲精品久久成人aⅴ小说| 有码 亚洲区| 成人综合一区亚洲| 十分钟在线观看高清视频www| 国产熟女午夜一区二区三区| 侵犯人妻中文字幕一二三四区| 国产永久视频网站| www日本在线高清视频| a级片在线免费高清观看视频| 国产在线视频一区二区| 欧美少妇被猛烈插入视频| 欧美变态另类bdsm刘玥| 男女高潮啪啪啪动态图| 成年动漫av网址| 国产亚洲午夜精品一区二区久久| 又黄又爽又刺激的免费视频.| 国产成人91sexporn| 亚洲美女搞黄在线观看| 日韩中字成人| 精品人妻一区二区三区麻豆| 搡女人真爽免费视频火全软件| 韩国精品一区二区三区 | 久久青草综合色| 不卡视频在线观看欧美| 久久热在线av| 伦精品一区二区三区| 午夜免费鲁丝| 久久久久精品人妻al黑| 久久国产亚洲av麻豆专区| 人妻一区二区av| 色94色欧美一区二区| 丝袜美足系列| 国产精品久久久久久久电影| 午夜福利影视在线免费观看| 精品国产露脸久久av麻豆| 亚洲国产最新在线播放| 午夜福利在线观看免费完整高清在| 亚洲精品日本国产第一区| a级毛片在线看网站| 狂野欧美激情性bbbbbb| 成人毛片60女人毛片免费| 九色亚洲精品在线播放| 香蕉丝袜av| 久久精品国产综合久久久 | 草草在线视频免费看| 亚洲国产成人一精品久久久| 美女主播在线视频| 久久精品国产鲁丝片午夜精品| 免费看光身美女| 一级a做视频免费观看| 一级黄片播放器| 国产熟女午夜一区二区三区| 久久精品国产鲁丝片午夜精品| 亚洲精品乱久久久久久| 午夜影院在线不卡| 国产亚洲欧美精品永久| 亚洲精品一二三| 在线观看国产h片| 亚洲久久久国产精品| 日日爽夜夜爽网站| 国产av精品麻豆| 五月天丁香电影| 母亲3免费完整高清在线观看 | 亚洲国产精品一区三区| 黄色一级大片看看| 日韩欧美精品免费久久| 男的添女的下面高潮视频| 免费看av在线观看网站| 五月伊人婷婷丁香| 午夜福利视频在线观看免费| 国产激情久久老熟女| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 免费观看a级毛片全部| 亚洲少妇的诱惑av| 亚洲精品色激情综合| 免费看不卡的av| 91精品伊人久久大香线蕉| 亚洲成人av在线免费| 国产老妇伦熟女老妇高清| 婷婷色麻豆天堂久久| 如何舔出高潮| 菩萨蛮人人尽说江南好唐韦庄| 国产色婷婷99| 精品一区在线观看国产| 99国产精品免费福利视频| 在线观看免费高清a一片| 免费观看av网站的网址| 国产黄色视频一区二区在线观看| freevideosex欧美| 在线亚洲精品国产二区图片欧美| 中国三级夫妇交换| 少妇被粗大的猛进出69影院 | 激情视频va一区二区三区| 亚洲av成人精品一二三区| 午夜福利视频在线观看免费| 夫妻午夜视频| 深夜精品福利| 91精品国产国语对白视频| 国产有黄有色有爽视频| 欧美日韩av久久| 亚洲色图 男人天堂 中文字幕 | 亚洲精品国产色婷婷电影| 91成人精品电影| 日韩三级伦理在线观看| 熟女av电影| 91精品国产国语对白视频| 国产男女超爽视频在线观看| 一本大道久久a久久精品| 精品国产一区二区三区久久久樱花| 一级毛片我不卡| a级毛片在线看网站| 亚洲av免费高清在线观看| 中文天堂在线官网| 国产精品国产三级国产专区5o| 69精品国产乱码久久久| 男女午夜视频在线观看 | 女人被躁到高潮嗷嗷叫费观| 亚洲国产av新网站| 成人综合一区亚洲| 久久ye,这里只有精品| 伦精品一区二区三区|