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

    WAVE DISSIPATING PERFORMANCE OF AIR BUBBLE BREAKWATERS WITH DIFFERENT LAYOUTS*

    2010-05-06 08:22:18ZHANGChengxing

    ZHANG Cheng-xing

    State Key Laboratory of Coastal and Offshore Engineering, Dalian University of Technology, Dalian 116024, China

    College of Urban Planning and Environmental Science, Xuchang University, Xuchang 461000, China, E-mail: zhangchengxing19@163.com

    WANG Yong-xue, WANG Guo-yu

    State Key Laboratory of Coastal and Offshore Engineering, Dalian University of Technology, Dalian 116024, China

    YU Long-mei

    Dalian Port Construction Supervision and Consultation Co. Ltd., Dalian 116015, China

    WAVE DISSIPATING PERFORMANCE OF AIR BUBBLE BREAKWATERS WITH DIFFERENT LAYOUTS*

    ZHANG Cheng-xing

    State Key Laboratory of Coastal and Offshore Engineering, Dalian University of Technology, Dalian 116024, China

    College of Urban Planning and Environmental Science, Xuchang University, Xuchang 461000, China, E-mail: zhangchengxing19@163.com

    WANG Yong-xue, WANG Guo-yu

    State Key Laboratory of Coastal and Offshore Engineering, Dalian University of Technology, Dalian 116024, China

    YU Long-mei

    Dalian Port Construction Supervision and Consultation Co. Ltd., Dalian 116015, China

    (Received December 28, 2009, Revised March 30, 2010)

    The wave dissipating performance of air bubble breakwaters with different layouts is studied by experimental and numerical methods in this article. Based on the assumpation that the mixture of air and water is regarded as a variable density fluid, the mathematical model of the air bubble breakwater is built. The numerical simulation results are compared with the experimental data, which shows that the mathematical model is reasonable for the transmission coefficient Ctm. The influencing factors are studied experimentally and numerically, including the incident wave heightHi, the incidentt wave period T , the air amount Qm, the submerged pipe depth D and the single or double air discharging pipe structure. Some valuable conclusions are obtained for further research of the mechanism and practical applications of air bubble breakwaters.

    Air bubble breakwater, wave dissiptating performance, experiments, variable density fluid

    1. Introduction

    The air bubble breakwater is composed of an air compressor and pipes with orifices. The pipes are often placed at the sea bottom, or at a certain depth in a deep water area. The air compressor is usually installed on the shore or on a boat[1]. The air bubble breakwater is used for protecting a harbor, the entrance of a port, a part of sea water against wave induced wrecks. The air bubble breakwater is one of mobile breakwaters, with several specific features such as mobility, temporality, and low cost. Besides, it may be pointed out that, with the use of the air bubble breakwater, from an environment point of view, one need not interrupt the exchange of water in harbors[2].

    The first air bubble breakwater appeared in 1907, which was used by Brasher of New York to protect civil engineering works. The second was built in 1915 at El Segundo, California, under Brasher’s patent. With some more projects of that nature, the results, however, were not promising and the method fell into disuse. Fundamental tests were conducted by the Admiralty in 1924, by Russian scientists in 1935 and by Professor Thysse of Delft in 1936. During the Second World War, the research and analysis on air bubble and water jet breakwaters were carried out by White and Taylor[3].

    In the early 1950’s, a wide publicity was given to an air bubble breakwater designed by Laurie to protect a train ferry dock entrance at Dover. Before that, Carr and Schiff studied the air bubbles breakwater in a model tank at the California Institute of Technology. With the use of the Dover installation, Evans conducted a series of model tests in 1954 in conjunction with the theoretical work by Taylor. Bulson reviewed the analytical and experimental studies carried out by himself and others, and in was concluded that the surface currents produced by the air bubble motion were the main mechanism of wave dissipation in the system of the air bubble breakwater. The empirical formulae were presented for the surface velocity and the thickness of the horizontal current produced by an air bubble curtain; and for the amount of free air required to suppress waves with certain length and height. Zhang et al.[4]studied the horizontal current generated by air bubble curtain in still water with a numerical simulation technique, which gave a further explanation of the wave dissipating mechanism of the air bubble breakwater. Furthermore, Zhang[5]carried out a preliminary study of the performance of the air bubble breakwater, and the results demonstrated the effects of air amount as well as incident wave period on the air bubble breakwater system.

    Despite of those experimental and theoretical researches done in past years, due to the complexity of the interactions between air bubbles and water, the air bubble breakwater remains a topic of studies and the practical application data are very much in demand. This article presents our experimental and numerical studies of the wave dissipating performance of the air bubble breakwaters with different layouts.

    Fig.1 The schematic diagram of the tests for the air bubble breakwater

    2. Experiment

    2.1 Equipment

    The length scales 1:10 and 1:15 are adopted in our experiments in the large wave tank of the State Key Laboratory of Coastal and Offshore Engineering, Dalian University of Technology. The large wave tank is 69.0 m long, 2.0 m wide and 1.8 m deep. The pipe to discharge air with orifices is made of plexiglass, with the length of 2.0 m, the orifice diameter of 0.0008 m and the distance between two orifices of 0.01 m. The compressed air is supplied by the air compressor up to the maximum amount of 3.0 m3/min. Figure 1 showes the schematic diagram of the tests for the air bubble breakwater.

    Four wave gauges of capacitance type are installed as shown in Fig.1. The transmission coefficient Ctmof the air bubble breakwater is calculated by the Goda method[6]. In the tests, the same case is tested for three times, to obtain the average value.

    Table 1 The parameters of tests with different length scales

    Fig.2 The apparatus of the air bubble breakwater system

    The prototype parameters include the incident wave height H=3.55 m, the incident wave period T=4s , 5 s and 6 s, and the water depth d=12 m. With respect to two length scales, the incident waves are governed by the law of gravity similarity, the corresponding parameters in model tests are shown in Table 1. The subscript of m and p means the model and prototype, respectively, in the article.

    During the tests, the air is pushed by the air compressor to the pipe, and then the air bubbles are produced continuously from the orifices and an air bubble curtain will appear across the tank width. The incident wave passes the air bubble curtain, and the wave dissipation occurs due to the wave deformation and the wave break. In the tests, the flow meter is installed at different locations along the pipe for the measurement of air amount Qm. Figure 2 shows the apparatus of the air bubble breakwater.

    2.2 Results and discussions

    The incident wave height Hiand the transmitive wave height Htcan be obtained from the test data by the Goda’s method. The transmission coefficient Ctmis expressed as Ctm=Ht/Hi. The impacts of the incident wave period T and the air amount Qmon the wave dissipating performance of the air bubble breakwater are compared and analyzed.

    Figure 3 shows the relations of the transmission coefficient Ctmversus the air amount Qmwith different length scales and different incident wave periods T. It is found that the wave dissipating performance of the air bubble breakwater is closely related with the air amout Qm, and it can be improved by the increase of air amount Qm.

    From Fig.3, it is seen that the incident wave period T is the key factor influencing the wave dissipating performance of the air bubble breakwater. Take Fig.3(a) for example, the incident waves with three different wave periods are dissipated with various air amounts Qm. In the case of the incident wave period T =1.26s, the transmission coefficient Ctmchanges from 0.87 to 0.48 when the air amount increases from 5 m3/h·m to 20 m3/h·m. However, the

    transmission coefficient Ctmchanges but little for the cases of the incident wave period T =1.58s and T =1.91s. It is demonstrated that the wave dissipating performance of the air bubble breakwater is affected evidently for short period waves with a given air amount Qm, and the wave dissipating performance of the air bubble breakwater is good for short period waves.

    Fig.3 Test data of Ctmagainst different air amountsQm

    3. Numerical simulation

    3.1 Mathematical formulation

    The mixture of air and water is considered as a fluid of variable density. The density and the dynamic viscosity coefficient of the fluid are defined as

    where a0,ρ0, μ0are the volume fraction, density and the dynamic viscosity coefficient of water, and a1,ρ1, μ1are the corresponding parameters for air. The volume fractions of both phases aq(q =0,1) can be solved from Eqs.(3) and (4)

    The continutiy equation is expressed as Eq.(5), and the momentum equation is the Reynolds averaging equation in the form[7-10], as shown in Eqs.(6) and (7), in which the density and the dynamic viscosity coefficient of the fluid is defined by Eqs.(1) and (2). The pressure pattern and the velocity pattern take the average values.

    where μtis the turbulent viscosity coefficient, μ =ρC k2/ε, C is a constant, C=0.09, the

    tμμμ turbulent kinetic energy k and the turbulent diffusionε are governed by Eqs.(8) - (9)[11,12]. As mentioned previously, it is assumed that the mixture of air and water is regarded as a variable density single phase liquid, so the standard k?ε equation can be used in this condition.

    where Gkis the turbulent kinetic energy produced

    Fig.4 A schematic diagram of the numerical simulation model of the air bubble breakwater

    3.3 Model validation

    In order to verify the rationality of the numerical model of the air bubble breakwater presented in the article, numerical computations are conducted and the schematic diagram of the numerical simulation model is shown in Fig.4. The length scales 1:10 and 1:15 are selected, the wave parameters and the air amounts Qmare the same as the experimental cases, as shown in Table 1. transmission coefficient Ctmof the test results would be smaller than those in other cases.

    It is seen from Figs.5 and 6 that the air amount Qmplays an important role in the wave dissipating performance of the air bubble breakwater, at the same time, the transmission coefficient Ctmis reduced with the increase of air amout Qm. It is also found that the incident wave period T affects the performance of wave dissipating, and the air bubble breakwater would have good wave dissipating performance when the incident wave has a short period T. All these are verified by experimental results.

    3.4 Results and discussions

    The factors that influence the wave dissipating performance of the air bubble breakwater are studied by numerical simulations, including the incident wave height H, the pipe submerged depth D and the spacing between two air discharging pipes dsfor the air bubble breakwater with double air discharging pipes. Based on the results obtained from numerical simulations, various influencing factors are analyzed to determine their impacts on the wave dissipating performance of the air bubble breakwater. A schematic diagram of the numerical wave tank is shown in Fig.4. The wave tank is 330 m long and 15 m deep, and the grids are refined in the region from X =110 m to X =130 m, with the air bubble curtain located at X =120 m. The transmission coefficientCtpis obtained at the monitor #3 as the numerical results.

    3.4.1 Incident wave height

    The cases considered in the numerical simulations include various incident wave heights H, incident wave period T and air amount Qpas shown in Table 2, and the results are illustrated in Fig.7.

    Figure 7(a) shows the variation of the transmission coefficient Ctpwith incident wave heights Hiand the air amount Qpfor T=4s. It is found that the transmission coefficientCtpdecreases with the increase of air amount Qp. But the reduction of transmission coefficient Ctpvaries with different incident wave heights Hi, the smaller the incident wave height Hi, the smaller the transmission coefficientCtpbecomes. The same results can be seen from the Figs.7(b) and 7(c) for T=5s and T=6s.

    Fig.5 The comparison of Ctmbetween experimental experimental and numerical simulations(Length scale, 1:10)

    Fig.6 The comparison of Ctmbetween and numerical simulations (Length scale, 1:15)

    Figures 5 and 6 show the comparisons of transmission coefficient Ctmbetween experimental and numerical results. It can be seen from Figs.5 and 6 that the maximum deviation of transmission coefficient Ctmbetween numerical results and experimental results is 9% except the case of the length scale 1:10, the incident wave periodT =1.26s and the air amount Q=20 m3/h· m . In the

    m experiments, the incident wave will be broken with the horizontal flow produced by a large air amount Qm, but the wave broken phenomenon is not considered in the numerical simulation. So the

    Table 2 Cases for numerical simulations with different incident wave heights Hi

    Fig.7 Transmission coefficient Ctpversus air amounts Qpfor different incident wave heights Hi

    It is demonstrated from Fig.7 that the transmission coefficient Ctpis reduced by a large amount with a small incident wave height Hifor the same air amountQp. which means that the wave dissipating performance of the air bubble breakwater is good for the small incident wave height Hiwith a given air amount Qp, the reason of which may be illustrated from the view of wave energy.

    When the water depth d and the incident wave period T are constants, the greater the incident wave height Hi, the greater the wave energy will be. The wave dissipating performance of the air bubble breakwater is not good when the incident wave height is large at a certain air amountQp.

    3.4.2 Submerged pipe depth

    The different submerged pipe depths D are considered with different air amounts Qpand different incident wave periods. The submerged pipe depth is defined as the vertical distance from the location of air discharging pipe in the water to the water surface. Three submerged pipe depths are considered in the article. The cases with different parameters in numerical simulations are shown in Table 3.

    Figure 8 shows the variation of the transmission coefficient Ctpwith the increase of air amount Qpfor different submerged pipe depths D and the incident wave periods T=4 s, 5 s and 6 s, respectively.

    Figure 8(a) shows the variation of the transmission coefficient Ctpversus the air amount Qpfor various submerged pipe depths D. It is found that the transmission coefficient Ctpdecreases with the increase of air amount Qpfor different submerged pipe depths D. The greater the submerged pipe depth D , the smaller the transmission coefficient Ctpwill be. The sameresults can be seen from Figs.8(b) and 8(c). transmission coefficient Ctpdecreases with the increase of the submerged pipe depth D for the same air amount Qp, which means that the submerged pipe depth D is one of influencing factors related to wave dissipation in the air bubble breakwater system. It is also seen that when the incident wave period T is short, the effect of the submerged pipe depth D is significant. Furthermore, the air amount Qpplays an important role in the wave dissipating performance of the air bubble breakwater, which is consistent with the experiment results.

    Table 3 The numerical simulation cases with different submerged depths of air discharging pipe D

    Fig.8 Transmission coefficient Ctpversus air amounts Qpfor different submerged pipe depths D

    Fig.9 The profile of Vmversus water depth ( X =134 m, t=67 s , Q=500 m3/h· m )

    Figure 9 shows the variation of the horizontal flow velocity Vmgenerated by the air bubble curtain versus the water depth d at the locationX =134 m in the calm water, when the water depth d =12 m,air amount Q=500 m3/h· m , the time

    p

    t=67 s and the submerged pipe depthD=6 m and 12 m, respectively. It can be found that the flow velocity Vmproduced by the pipe with the submergence D=12 m is greater than that with the submergence D=6 m under the same air amountQpin the calm water. According to the wave dissipating mechanism of the air bubble breakwater, as the effects of the horizontal flow velocity Vmon the incident waves, the greater the horizontal flow velocity Vm, the better the wave dissipating performance will be. It is seen that the wave dissipation of the air bubble breakwater is affected by the submerged pipe depthD.

    3.4.3 Double air discharging pipes

    The air bubble breakwater with double air discharging pipes is designed and its schematic diagram is shown in Fig.10 for studying the wave dissipating performance of the air bubble breakwater. Two different spacings between air discharging pipes ds=6 m and 10 m are considered in the article. The locations of air discharging pipes are atX =117 m and X=123m for ds=6 m and atX =115 m and X =125 m for ds=10 m , respectively.

    Fig.10 Schematic diagram of the air bubbles breakwater with double air discharging pipes

    The wave dissipating performance of the air bubble breakwater with double air discharging pipes are studied, considering different spacings between double air discharging pipes dsand different air amounts Qp, which are taken as the total air amounts of both pipes for the case of double air discharging pipes. The incident wave period T=5s, the incident wave height H=3.55 m and the water depth d=12 mare adopted in the numerical simulations. The total air amountQpis selected from 100 m3/h·m to 580 m3/h·m with increment of 80 m3/h·m.

    Figure 11 shows the comparison of transmission coefficient Ctpbetween single and double air discharging pipe structures with various total air amounts Qp, where (a) for ds=6 m and (b) for ds=10 m . It is found from Fig.11(a) that the transmission coefficientCtpof the single and double air discharging pipe structures sees no significant difference for various total air amounts Qpfor ds=6 m , which means that the wave dissipating effects of the air bubbles breakwaters of both structures are similar. It is also found from Fig.11(b) that the transmission coefficient Ctpof the single air discharging pipe structure is smaller compared with that of the double air discharging pipe structure for ds=10 m for various total air amountsQp, which means that the spacing between double air discharging pipes affects the wave dissipating performance to a certain extent.

    Fig.11 The comparison of transmission coefficient Ctpbetween single and double air discharging pipe structures with various total air amounts Qp

    Figure 12 shows the pattern of velocity vector generated by the double air discharging pipes in the calm water at the time t=67 s with the air amount Q=500 m3/h· m , the spacing d=6 m and 10 m,

    p s and the water depth d=12 m. Figure 13 shows the distribution of the horizontal flow velocity Vmalong the water depth at the position X =134 m with the the same condition referred above. It is found from Fig.12(a) that the pattern of the circulation flow produced by the air bubble’s motion is similar to that by the single air discharging pipe for the spacing ds=6 m , so the wave dissipating capability isnearly the same between two kinds of structures. But for the spacing ds= 10 m, on the other hand, as shown in Fig.12(b), there are two minor vortexes between twoair discharging pipes. The existence of the two minor vortexes affects the distribution of the part of the total energy generated by the air bubble’s motion, which would reduce the horizontal flow velocity Vm. Figure 13 also shows obviously that the horizontal flow velocityVmdecreases with the increase of the spacing ds. It is seen that the wave dissipating performance of the air bubble breakwater with double air discharging pipes may not be improved as compared with the air bubble breakwater with a single air discharging pipe, especially, when the spacing between double air discharging pipes dsis large.

    Fig.12 The pattern of velocity vector with different spacings between double air discharging pipes dsin the calm water (t =67 s , Q=500 m3/h· m )p

    Fig.13 The profile of Vmversus water depth ( X =134 m, t=67 s ,Q=500 m3/h· m )p

    4. Conclusions

    In the article the wave dissipating performance of the air bubble breakwater and the influencing factors are studied by experiments and numerical computations. In the numerical model, the mixture of water and air is regarded as a variable density liquid, the Reynolds averaging equation and the standard k?ε equations are adopted as the governing equations, the method of VOF is used to track the two-phase interface; an additional mass source is added in the continuity equation and an additional momentum source is added in the momentum equations. The following conclusions are drawn.

    (1) By the comparison of experimental results and numerical simulation results, it is demonstrated that the mathematical model in this article is reasonable.

    (2) The wave dissipating performance of the air bubble breakwater depends closely on the air amount Q. With the increase of the air amount Q, the wave dissipating performance of the air bubble breakwater is improved.

    (3) The wave dissipating performance of the air bubble breakwater is affected significantly by the incident wave period T, especially, for the incident wave with short period T.

    (4) With the increase of the incident wave height Hi, the wave dissipating performance of the air bubble breakwater becomes worse.

    (5) The increase of the submerged pipe depth D would improve the wave dissipating performance of the air bubble breakwater.

    (6) With the same total air amount Q, the wave dissipating performance of the air bubble breakwater with double air discharging pipes may not be improved in comparison with the air bubble breakwater with a single air discharging pipe, especially, when the spacing between the double air discharging pipesdsis large. It may be a good choice to use the single air discharging pipe structure instead of the double air discharging pipe structure.

    [1] WANG Guo-yu. Investigation on the structure type and performance of the special breakwaters[D]. Ph. D. Thesis, Dalian: Dalian University of Technology, 2005(in Chinese).

    [2] TONG Chao-feng, YANG Yi-xin. Study on the decaying wave feature of floating breakwater[J]. Port and Waterway Engineering, 2002, 343(8): 32-35(in Chinese).

    [3] WANG Guo-yu, WANG Yong-xue and LI Guang-wei. Experimental study on the performance of air bubbles breakwater[J]. Shipbuilding of China, 2004, 45(z1): 103-109(in Chinese).

    [4] ZHANG Cheng-xing, WANG Yong-xue and WANG Guo-yu et al. Numerical simulation study on the horizontal current generated by air bubbles curtain in still water[J]. Chinese Journal of Hydrodynamics, 2010, 25(1): 59-66(in Chinese).

    [5] ZHANG Cheng-xing, WANG Guo-yu and WANG Yong-xue. Study on the peirformances of air curtain breakwater by numerical simulation study[J]. Chinese Journal of Hydrodynamics, 2009, 24(5): 543-549(in Chinese).

    [6] YU Yu-xiu. Random wave theories and engineering applications[M]. Dalian: Dalian University of Technology Press, 2000(in Chinese).

    [7] HAN Zhan-zhong, WANG Jing and LAN Xiao-ping. Fluids engineering simulation examples and applications of Fluent[M]. Beijing: Beijing Institute of Technology Press, 2004(in Chinese).

    [8] CHANG K. A., HSU T. J. and LIU P. L. F. Vortex generation and evolution in water waves propagating over a submerged rectangular obstacle: Part II: Cnoidal waves[J]. Coastal Engineering, 2005, 52(3): 257-283.

    [9] CHEN G., KHARIF C. and ZALESKI S. et al. Two-dimensional Navier-Stokes simulation of breaking waves[J]. Physics of Fluids, 1999, 11(1): 121-133.

    [10] LUBIN P., VINCENT S. and ABADIE S. et al. Three-dimensional large eddy simulation of air entrainment under plunging breaking waves[J]. Coastal Engineering, 2006, 53(8): 631-655.

    [11] LIU Zhen, HYUN Beom-soo and KIM Moo-rong et al. Experimental and numerical study for hydrodynamic characteristics of an oscillating hydrofoil[J]. Journal of Hydrodynamics, 2008, 20(3): 280-287.

    [12] WANG Rui-jin, ZHANG Kai and WANG Gang. The technical bases and application examples of Fluent[M]. Beijing: Tsinghua University Press, 2007(in Chinese).

    [13] LI Ling. The numerical simulation of interaction of water waves and floating structures in a viscous fluid[D]. Master Thesis, Shanghai: Shanghai Jiao Tong University, 2007(in Chinese).

    [14] LI Sheng-zhong. Study on 2-D numerical wave tank based on the software Fluent[D]. Master Thesis, Harbin: Harbin Institute of Technology, 2006(in Chinese).

    [15] LI Ling, LIN Zhao-wei and YOU Yun-xiang et al. The numerical wave flume of the viscous fluid based on the momentum source method[J]. Journal of Hydrodynamics, Ser. A, 2007, 22(1): 76-82(in Chinese).

    [16] WANG Yan, YOU Yun-xiang. Numerical simulation of interaction of viscous wave fields with a semi-submersible platform[J]. Chinese Journal of Hydrodynamics, 2009, 24(6): 793-799(in Chinese).

    [17] ZHANG Hong-sheng, YU Xiao-wei and YANG Jian-min et al. Verification and application of the improved numerical model for nonlinear wave propagation[J]. Chinese Journal of Hydrodynamics, 2009, 24(3): 364-373(in Chinese).

    [18] DONG Zhi, ZHAN Jie-min. Comparison of existing methods for wave generating and absorbing in VOF-based numerical tank[J]. Journal of Hydrodynamics, Ser. A, 2009, 24(1): 15-21(in Chinese).

    [19] LU Chang-na, WANG Ru-yun and CHEN Ping-ping. A keeping volume fraction method for moving interfaces reconstruction of VOF using triangle meshes[J]. Chinese Journal of Hydrodynamics, 2008, 23(3): 255-260(in Chinese).

    [20] DONG Zhi, ZHAN Jie-min. Numerical modeling of wave evolution and runup in shallow water[J]. Journal of Hydrodynamics, 2009, 21(6): 731-738.

    10.1016/S1001-6058(09)60102-5

    * Project supported by the National Natural Science Foundation of China (Grant No. 50809015).

    Biography: ZHANG Cheng-xing (1977-), Male, Ph. D., Lecturer

    99在线视频只有这里精品首页| 内地一区二区视频在线| 国产伦理片在线播放av一区 | 3wmmmm亚洲av在线观看| av免费在线看不卡| а√天堂www在线а√下载| 久久午夜亚洲精品久久| 精品人妻熟女av久视频| 亚洲天堂国产精品一区在线| 熟妇人妻久久中文字幕3abv| 99国产极品粉嫩在线观看| 日日撸夜夜添| 免费av毛片视频| 亚洲精品国产成人久久av| 丝袜美腿在线中文| 国产一级毛片在线| 午夜免费男女啪啪视频观看| 国内少妇人妻偷人精品xxx网站| 亚洲成人中文字幕在线播放| 欧美性猛交╳xxx乱大交人| 国产精品久久视频播放| 欧美成人免费av一区二区三区| 99热这里只有是精品50| 久久久久久久久久久免费av| av国产免费在线观看| 变态另类丝袜制服| 亚洲四区av| 免费黄网站久久成人精品| 国产成人午夜福利电影在线观看| 一区二区三区高清视频在线| av免费在线看不卡| 欧美成人一区二区免费高清观看| 久久精品国产鲁丝片午夜精品| 国产毛片a区久久久久| 亚洲国产精品sss在线观看| 91久久精品国产一区二区三区| 国产真实乱freesex| 成人综合一区亚洲| 此物有八面人人有两片| 观看美女的网站| 国产精品一二三区在线看| 亚洲成av人片在线播放无| 男女啪啪激烈高潮av片| 美女xxoo啪啪120秒动态图| 美女 人体艺术 gogo| 国产欧美日韩精品一区二区| 天天躁夜夜躁狠狠久久av| 丰满人妻一区二区三区视频av| 午夜亚洲福利在线播放| 人妻系列 视频| 99久久成人亚洲精品观看| 国产伦在线观看视频一区| 老熟妇乱子伦视频在线观看| 悠悠久久av| 久久99热这里只有精品18| 美女xxoo啪啪120秒动态图| 国产精品一区www在线观看| 精品久久久久久久久久久久久| 内射极品少妇av片p| 久久热精品热| 亚洲国产欧美人成| 精品一区二区三区视频在线| 国产在视频线在精品| 日日摸夜夜添夜夜添av毛片| 91aial.com中文字幕在线观看| 久久人人爽人人爽人人片va| 一级黄片播放器| 青春草视频在线免费观看| 欧美激情久久久久久爽电影| 1000部很黄的大片| 国产男人的电影天堂91| 亚洲精品国产av成人精品| 亚洲av电影不卡..在线观看| 可以在线观看毛片的网站| 成年版毛片免费区| 在线免费观看的www视频| 一级毛片我不卡| 啦啦啦啦在线视频资源| 精品久久久久久成人av| 哪里可以看免费的av片| 99热只有精品国产| 欧美一区二区亚洲| 成人无遮挡网站| 国产美女午夜福利| 国产黄色视频一区二区在线观看 | 欧美性猛交黑人性爽| a级毛片a级免费在线| 一区二区三区四区激情视频 | 极品教师在线视频| 老司机福利观看| 精品久久久久久久久久久久久| 女人被狂操c到高潮| 成熟少妇高潮喷水视频| 亚洲图色成人| 国产色婷婷99| 亚洲av中文字字幕乱码综合| 99久久精品一区二区三区| 在线播放无遮挡| 久久久久久久久久久免费av| 成年女人看的毛片在线观看| 欧美+亚洲+日韩+国产| 亚洲精品亚洲一区二区| 日韩精品有码人妻一区| 国产成人福利小说| 亚洲真实伦在线观看| 亚洲精品日韩av片在线观看| 亚洲欧美精品专区久久| 看非洲黑人一级黄片| 精品国产三级普通话版| 91麻豆精品激情在线观看国产| 欧美精品国产亚洲| АⅤ资源中文在线天堂| 91久久精品国产一区二区三区| 校园人妻丝袜中文字幕| 日韩强制内射视频| 日日干狠狠操夜夜爽| 久久国内精品自在自线图片| 午夜免费激情av| 色尼玛亚洲综合影院| 少妇熟女aⅴ在线视频| 国产亚洲5aaaaa淫片| 国国产精品蜜臀av免费| 亚洲最大成人手机在线| 日韩视频在线欧美| 亚洲国产高清在线一区二区三| 神马国产精品三级电影在线观看| 国产视频首页在线观看| 色视频www国产| 国产午夜精品一二区理论片| av在线观看视频网站免费| or卡值多少钱| 国产黄色视频一区二区在线观看 | 人人妻人人看人人澡| 悠悠久久av| 长腿黑丝高跟| 丝袜美腿在线中文| 美女高潮的动态| 成人永久免费在线观看视频| 村上凉子中文字幕在线| 99久久无色码亚洲精品果冻| 亚洲国产精品久久男人天堂| h日本视频在线播放| 在线观看一区二区三区| 欧美成人a在线观看| 午夜福利成人在线免费观看| 天天一区二区日本电影三级| 精品一区二区免费观看| 91在线精品国自产拍蜜月| 青春草国产在线视频 | 99久久人妻综合| www.av在线官网国产| 黄片wwwwww| 久久午夜福利片| 一边亲一边摸免费视频| 91久久精品国产一区二区成人| 级片在线观看| 中文欧美无线码| 日韩av在线大香蕉| 欧美一级a爱片免费观看看| 日韩人妻高清精品专区| 国产成人a∨麻豆精品| 男女那种视频在线观看| 大型黄色视频在线免费观看| 性插视频无遮挡在线免费观看| 久久久久久久亚洲中文字幕| 精华霜和精华液先用哪个| 最后的刺客免费高清国语| 一级毛片aaaaaa免费看小| 亚洲国产精品sss在线观看| 日韩高清综合在线| 色5月婷婷丁香| 丰满的人妻完整版| 亚洲一区高清亚洲精品| 成年女人看的毛片在线观看| 日日摸夜夜添夜夜爱| 久久久久久久亚洲中文字幕| 少妇被粗大猛烈的视频| 国产美女午夜福利| 精品不卡国产一区二区三区| 国产一区二区亚洲精品在线观看| 在线播放国产精品三级| 国产精品福利在线免费观看| 毛片女人毛片| 国产白丝娇喘喷水9色精品| 少妇的逼好多水| 性欧美人与动物交配| 全区人妻精品视频| 久久综合国产亚洲精品| 国产日韩欧美在线精品| 亚洲美女视频黄频| 国产免费一级a男人的天堂| or卡值多少钱| 国产日韩欧美在线精品| 日韩欧美精品v在线| 日本一二三区视频观看| 久久久久国产网址| 亚洲精品国产成人久久av| 蜜桃久久精品国产亚洲av| 悠悠久久av| av在线播放精品| 亚洲在线观看片| 欧美日韩在线观看h| 亚洲精品久久国产高清桃花| 搞女人的毛片| 久久久久久大精品| 日本一本二区三区精品| 国产 一区精品| 最近手机中文字幕大全| www.av在线官网国产| 精品久久久久久久久av| 天堂影院成人在线观看| av国产免费在线观看| 99在线视频只有这里精品首页| 少妇熟女aⅴ在线视频| 99久久中文字幕三级久久日本| 午夜激情欧美在线| 欧美三级亚洲精品| 色播亚洲综合网| 欧美潮喷喷水| 国产成人91sexporn| 精品无人区乱码1区二区| 在线观看66精品国产| 亚洲精品国产成人久久av| 国产乱人偷精品视频| 国产极品天堂在线| 99在线视频只有这里精品首页| 久久精品国产亚洲av涩爱 | 久久午夜亚洲精品久久| 国产视频内射| 麻豆一二三区av精品| 91久久精品电影网| 亚洲国产欧美在线一区| 一边摸一边抽搐一进一小说| 蜜桃久久精品国产亚洲av| av福利片在线观看| 久久这里只有精品中国| 一区二区三区免费毛片| 国产乱人偷精品视频| 岛国毛片在线播放| 久久精品国产亚洲av香蕉五月| 嫩草影院入口| 99久久成人亚洲精品观看| 午夜福利在线观看免费完整高清在 | 午夜福利高清视频| 国产视频首页在线观看| 免费看日本二区| 亚洲经典国产精华液单| 在线观看美女被高潮喷水网站| 亚洲av成人精品一区久久| 日韩中字成人| 欧美xxxx性猛交bbbb| 在线a可以看的网站| 嫩草影院精品99| 真实男女啪啪啪动态图| 欧美日韩精品成人综合77777| 国产精品一区www在线观看| 变态另类成人亚洲欧美熟女| 一区二区三区四区激情视频 | 岛国在线免费视频观看| 亚洲,欧美,日韩| 免费不卡的大黄色大毛片视频在线观看 | 99精品在免费线老司机午夜| 色哟哟哟哟哟哟| 黑人高潮一二区| 国产精品一区www在线观看| 人人妻人人看人人澡| 欧美日本亚洲视频在线播放| 一个人免费在线观看电影| 欧美另类亚洲清纯唯美| 美女脱内裤让男人舔精品视频 | 乱系列少妇在线播放| 成人午夜精彩视频在线观看| 国产三级中文精品| 免费看美女性在线毛片视频| 中文字幕免费在线视频6| 小说图片视频综合网站| 中文精品一卡2卡3卡4更新| 欧美色欧美亚洲另类二区| 女同久久另类99精品国产91| 精品久久国产蜜桃| 成人二区视频| 国产精品一区二区在线观看99 | 三级经典国产精品| 国产精品爽爽va在线观看网站| 日本一本二区三区精品| 亚洲一级一片aⅴ在线观看| 精品日产1卡2卡| 高清午夜精品一区二区三区 | 最近手机中文字幕大全| 尤物成人国产欧美一区二区三区| 日韩成人伦理影院| 免费看a级黄色片| 国产成人精品久久久久久| 久久亚洲国产成人精品v| 精品不卡国产一区二区三区| 美女脱内裤让男人舔精品视频 | 久久久色成人| 欧美激情在线99| 久久人妻av系列| 久久久a久久爽久久v久久| 国产色爽女视频免费观看| 人人妻人人看人人澡| 免费观看精品视频网站| 美女 人体艺术 gogo| 久久午夜福利片| 国产亚洲av嫩草精品影院| 99久久无色码亚洲精品果冻| 少妇裸体淫交视频免费看高清| 在现免费观看毛片| 麻豆av噜噜一区二区三区| 老女人水多毛片| 免费av观看视频| 精品一区二区三区人妻视频| 久久国产乱子免费精品| 日日撸夜夜添| 美女脱内裤让男人舔精品视频 | 寂寞人妻少妇视频99o| 国产女主播在线喷水免费视频网站 | 国产在线精品亚洲第一网站| 内射极品少妇av片p| 夜夜看夜夜爽夜夜摸| 欧美日韩精品成人综合77777| 有码 亚洲区| 99riav亚洲国产免费| 精品日产1卡2卡| 亚洲人成网站在线观看播放| 一个人看的www免费观看视频| 亚洲真实伦在线观看| 99热只有精品国产| 久久综合国产亚洲精品| 欧美xxxx性猛交bbbb| 观看免费一级毛片| 人妻制服诱惑在线中文字幕| 欧美另类亚洲清纯唯美| 尾随美女入室| 欧美高清成人免费视频www| 国产探花在线观看一区二区| 婷婷精品国产亚洲av| 精华霜和精华液先用哪个| 久久精品国产自在天天线| 亚洲国产欧美在线一区| 亚洲自偷自拍三级| 国产视频首页在线观看| 成人av在线播放网站| 亚洲不卡免费看| 欧美xxxx性猛交bbbb| 国产精品蜜桃在线观看 | 高清毛片免费观看视频网站| 18禁在线无遮挡免费观看视频| 美女内射精品一级片tv| 国产精品久久久久久av不卡| 欧美xxxx性猛交bbbb| 观看免费一级毛片| 久久久久久久久久久免费av| 又粗又爽又猛毛片免费看| 3wmmmm亚洲av在线观看| 爱豆传媒免费全集在线观看| 国产精品女同一区二区软件| 国产片特级美女逼逼视频| 免费黄网站久久成人精品| 白带黄色成豆腐渣| 欧美日本亚洲视频在线播放| 高清日韩中文字幕在线| 99九九线精品视频在线观看视频| 亚洲四区av| 十八禁国产超污无遮挡网站| 亚洲,欧美,日韩| 在线天堂最新版资源| 亚洲欧美日韩无卡精品| 黄色视频,在线免费观看| 秋霞在线观看毛片| 久久久欧美国产精品| 免费观看在线日韩| 赤兔流量卡办理| 国产极品精品免费视频能看的| 亚洲一级一片aⅴ在线观看| 18禁裸乳无遮挡免费网站照片| 级片在线观看| 一级黄色大片毛片| 欧美又色又爽又黄视频| 九九久久精品国产亚洲av麻豆| 99九九线精品视频在线观看视频| 亚洲人与动物交配视频| 麻豆一二三区av精品| 国产免费一级a男人的天堂| АⅤ资源中文在线天堂| 99久国产av精品国产电影| 久久人妻av系列| 永久网站在线| 精品午夜福利在线看| 好男人视频免费观看在线| 91精品国产九色| 国产探花极品一区二区| 97人妻精品一区二区三区麻豆| 国产视频首页在线观看| 三级毛片av免费| 亚洲欧美成人精品一区二区| 男的添女的下面高潮视频| 亚洲内射少妇av| 最近2019中文字幕mv第一页| 国内精品久久久久精免费| 1000部很黄的大片| 亚洲av成人av| 久久久国产成人免费| 在线播放无遮挡| 国产爱豆传媒在线观看| 欧美精品国产亚洲| 99在线视频只有这里精品首页| 精品久久久噜噜| 波多野结衣高清无吗| 亚洲欧美日韩无卡精品| av又黄又爽大尺度在线免费看 | 我的老师免费观看完整版| 精品久久久久久久久av| 色视频www国产| 亚洲成av人片在线播放无| 国产真实伦视频高清在线观看| 欧美一区二区亚洲| 欧美另类亚洲清纯唯美| 在线观看免费视频日本深夜| 国产伦精品一区二区三区视频9| 一本久久中文字幕| 亚洲av不卡在线观看| kizo精华| 国产成人精品久久久久久| 丰满人妻一区二区三区视频av| 亚洲四区av| 国产午夜精品论理片| 少妇丰满av| 中国美女看黄片| 精品人妻一区二区三区麻豆| av在线蜜桃| 又粗又硬又长又爽又黄的视频 | a级毛色黄片| 不卡视频在线观看欧美| 国产精品,欧美在线| 最近中文字幕高清免费大全6| 国产亚洲欧美98| 国产成人a区在线观看| 夜夜看夜夜爽夜夜摸| 成人特级av手机在线观看| 1024手机看黄色片| 一本精品99久久精品77| 成人特级黄色片久久久久久久| 久久精品综合一区二区三区| 婷婷色av中文字幕| 在线观看免费视频日本深夜| 午夜a级毛片| 午夜福利视频1000在线观看| 成人一区二区视频在线观看| 中文字幕av成人在线电影| 日韩成人av中文字幕在线观看| 如何舔出高潮| 亚洲一区二区三区色噜噜| 少妇的逼好多水| 婷婷色av中文字幕| 人体艺术视频欧美日本| 日韩在线高清观看一区二区三区| 免费看av在线观看网站| www.av在线官网国产| 国产一区亚洲一区在线观看| 岛国在线免费视频观看| 极品教师在线视频| 不卡一级毛片| 嫩草影院新地址| 精品99又大又爽又粗少妇毛片| 国产精品野战在线观看| 91久久精品电影网| 国产精品久久久久久精品电影| 欧美一区二区亚洲| 国产 一区 欧美 日韩| 国产精品福利在线免费观看| 欧美xxxx黑人xx丫x性爽| 看黄色毛片网站| avwww免费| 久久久久网色| 国内精品一区二区在线观看| 国产成年人精品一区二区| 久99久视频精品免费| 精品少妇黑人巨大在线播放 | 国产精品久久久久久亚洲av鲁大| 一边摸一边抽搐一进一小说| 国产精品不卡视频一区二区| 黄片无遮挡物在线观看| 久久九九热精品免费| 天天一区二区日本电影三级| 综合色av麻豆| 国产成人精品一,二区 | 99久国产av精品国产电影| 91狼人影院| 色吧在线观看| 亚洲欧美精品综合久久99| 日本黄大片高清| 岛国在线免费视频观看| 青春草视频在线免费观看| 12—13女人毛片做爰片一| 国产精品久久久久久av不卡| 欧美高清成人免费视频www| 国产一区二区在线观看日韩| 能在线免费观看的黄片| 一级黄色大片毛片| 我要搜黄色片| 成人一区二区视频在线观看| 成人国产麻豆网| 麻豆精品久久久久久蜜桃| 国产精品久久久久久久电影| 搡老妇女老女人老熟妇| 成人欧美大片| 日韩中字成人| 麻豆国产av国片精品| 白带黄色成豆腐渣| 综合色av麻豆| av在线老鸭窝| 亚洲成人精品中文字幕电影| 午夜福利视频1000在线观看| 国产探花极品一区二区| 有码 亚洲区| 精品久久久久久久久久久久久| 久久99精品国语久久久| 日本欧美国产在线视频| 免费看光身美女| 一级毛片电影观看 | 精品久久久久久久久av| 国产av一区在线观看免费| 国产午夜精品久久久久久一区二区三区| 亚洲欧美清纯卡通| 久久久久国产网址| 国产一区二区亚洲精品在线观看| 国产在视频线在精品| 精品人妻视频免费看| 亚洲国产精品成人久久小说 | 免费大片18禁| 国内精品久久久久精免费| 菩萨蛮人人尽说江南好唐韦庄 | 3wmmmm亚洲av在线观看| 亚洲成人久久爱视频| 插逼视频在线观看| 日日摸夜夜添夜夜爱| 欧美变态另类bdsm刘玥| av免费在线看不卡| 国产69精品久久久久777片| 两个人的视频大全免费| 精品久久久久久久末码| 国产精品一区二区三区四区免费观看| 成人漫画全彩无遮挡| 成人一区二区视频在线观看| 国产精品,欧美在线| 亚洲美女搞黄在线观看| 欧美日韩精品成人综合77777| 麻豆乱淫一区二区| 日本-黄色视频高清免费观看| 国产精品久久视频播放| 欧美激情久久久久久爽电影| 我要搜黄色片| 亚洲人成网站在线观看播放| 99久久中文字幕三级久久日本| 内地一区二区视频在线| 亚洲欧美日韩无卡精品| 亚洲成人久久爱视频| 全区人妻精品视频| 一个人观看的视频www高清免费观看| 日韩在线高清观看一区二区三区| 成年女人永久免费观看视频| 精品久久久久久久久久久久久| 国产精品av视频在线免费观看| 国产精品一二三区在线看| 又粗又爽又猛毛片免费看| 好男人在线观看高清免费视频| 全区人妻精品视频| videossex国产| 久久久久久久久久久免费av| 欧美性猛交╳xxx乱大交人| 成人一区二区视频在线观看| 男女下面进入的视频免费午夜| 一本一本综合久久| 亚洲在久久综合| 深夜精品福利| 99热6这里只有精品| 黄色日韩在线| 久久草成人影院| 亚洲人成网站在线播放欧美日韩| 国产精品三级大全| 蜜桃久久精品国产亚洲av| 亚洲一区二区三区色噜噜| 最近中文字幕高清免费大全6| 国产综合懂色| 精品人妻视频免费看| 亚洲欧美日韩高清专用| 国产成人freesex在线| 久久久久久久久久黄片| 性欧美人与动物交配| av.在线天堂| 青青草视频在线视频观看| 99热这里只有是精品50| 亚洲精品日韩av片在线观看| 亚洲欧美中文字幕日韩二区| 国产熟女欧美一区二区| 一边摸一边抽搐一进一小说| 亚洲成人av在线免费| 高清毛片免费观看视频网站| 成人特级黄色片久久久久久久| 欧美成人a在线观看| 夜夜看夜夜爽夜夜摸| 男人和女人高潮做爰伦理| 99久久久亚洲精品蜜臀av| 成人性生交大片免费视频hd| 麻豆久久精品国产亚洲av| 国产片特级美女逼逼视频| 日韩欧美 国产精品| 简卡轻食公司| 国产三级中文精品| 青春草国产在线视频 | 我的女老师完整版在线观看| av在线老鸭窝| 天堂√8在线中文|