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

    Simulation of buoyancy-induced turbulent flow from a hot horizontal jet*

    2014-06-01 12:30:00ElAMINSUNShuyuSALAMAmgad

    El-AMIN M. F., SUN Shuyu, SALAM Amgad

    Computational Transport Phenomena Laboratory (CTPL), Division of Physical Sciences and Engineering, King Abdullah University of Science and Technology, Thuwal 23955-6900, Kingdom of Saudi Arabia,

    E-mail: mohamed.elamin@kaust.edu.sa

    Simulation of buoyancy-induced turbulent flow from a hot horizontal jet*

    El-AMIN M. F., SUN Shuyu, SALAM Amgad

    Computational Transport Phenomena Laboratory (CTPL), Division of Physical Sciences and Engineering, King Abdullah University of Science and Technology, Thuwal 23955-6900, Kingdom of Saudi Arabia,

    E-mail: mohamed.elamin@kaust.edu.sa

    (Received September 14, 2013, Revised November 11, 2013)

    Experimental visualizations and numerical simulations of a horizontal hot water jet entering cold water into a rectangular storage tank are described. Three different temperature differences and their corresponding Reynolds numbers are considered. Both experimental visualization and numerical computations are carried out for the same flow and thermal conditions. The realizablekε- model is used for modeling the turbulent flow while the buoyancy is modeled using the Boussinesq approximation. Polynomial approximations of the water properties are used to compare with the Boussinesq approximation. Numerical solutions are obtained for unsteady flow while pressure, velocity, temperature and turbulence distributions inside the water tank as well as the Froude number are analyzed. The experimental visualizations are performed at intervals of five seconds for all different cases. The simulated results are compared with the visualized results, and both of them show the stratification phenomena and buoyancy force effects due to temperature difference and density variation. After certain times, depending on the case condition, the flow tends to reach a steady state.

    turbulent flow, realizablekε- turbulence model, heat transfer, jet, heat storage

    Introduction

    Thermal stratification in the storage tanks has a strong influence on the thermal performance of solar heating systems. The water entering the solar water storage as a jet causes mixing and destroys the stratifications. CFD-software has been used widely to simulate the components of solar heating systems. Knudsen et al.[1]studied the effect of thermal stratification phenomena in the storage tank by means of experiments and numerical computation. They concluded that the effect of thermal stratification in storage tanks is very important for the thermal performance. Shah and Furbo[2]presented a numerical and experimental analysis of water jets entering a solar storage tank. Three inlet designs with different inlet flow rates were simulated out to illustrate the varying behavior of the thermal conditions in a solar store. Their results showed how the inlet design influences the flow patterns in the tank and how the energy flux in a hot water tank is reduced with a poor inlet design. Velocity and temperature fields around a cold water inlet device of small solar domestic hot water tanks were investigated by Jordan and Furbo[3]using the CFD tool Fluent. The simulation results were compared with temperature measurements inside a commercial storage tank. Knudsen et al.[4]analyzed the flow structure and heat transfer in a vertical mantle heat exchanger. The flow structure and velocities in the inner tank and in the mantle were measured using a particle image velocimetry (PIV) system. In turbulent flow, the Reynoldsaveraged Navier-Stokes (RANS) tewchnique is usually adopted in order to make the system amenable to solution. The problem of using RANS approach, however, is that till now, there is no unifying set of equations to model all kinds of turbulent flows and heat transfer scenarios. Therefore, it is important to choose the model which suites the case under investigation and even to calibrate its coefficients in order to fit experimental results. The flow structure of a water jet entering into a rectangular storage tank was investigated by experimental visualizations as well as numerical CFD calculations for the problem have been introduced by El-Amin et al.[5]. El-Amin et al.[6]have pre-sented the 2-D upward, axisymmetric turbulent confined jet and developed several models to describe flow patterns using the realizablekε- turbulence model. Both realizable and RNGkε- turbulence models were later calibrated[7]. Several experimental works were conducted to highlight the interesting patterns and the governing parameters pertinent to this kind of flows[8,9]. For example, Arakeri et al.[10]introduced the problem of bifurcation in a buoyant horizontal laminar jet. O’Hern et al.[11]performed experimental work on a turbulent buoyant helium plume. El-Amin and Kanayama[12,13]studied buoyant jet resulting from hydrogen leakage. They developed the similarity formulation and solutions of the centerline quantities such as velocity and concentration. Moreover, El-Amin[14]performed a numerical investigation of a vertical axisymmetric non-Boussinesq buoyant jet resulting from hydrogen leakage in air as an example of injecting a low-density gas into high-density ambient. On the other hand, the mechanics of buoyant jet flows issuing with a general three-dimensional geometry into an unbounded ambient environment with uniform density or stable density stratification and under stagnant or steady sheared current conditions is investigated by Jirka[15,16]extended this work to also encounter plane buoyant jet dynamics resulting from the interaction of multiple buoyant jet effluxes spaced along a diffuser line.

    This paper introduces an analysis for three temperature differences with the three corresponding Reynolds numbers of horizontal hot water jets entering a rectangular solar storage tank filled with cold water. In order to investigate the stratification phenomena inside the tank, an experimental visualization was performed for the different cases in sequential time steps, followed by numerical investigations under the same conditions with wide-ranging of analyses for fields of pressure, velocity, temperature and turbulence inside the water store.

    Fig.1. Sketch of the water tank

    1. Experimental flow visualizations

    Visualization has been used to observe the flow patterns and thermal effects in a storage tank with dimensions in meterX×Y×Z=0.486 m×0.275 m× 0.29 m. The upper boundary is open with no cover for the tank, thus, the outflow of the tank is on the top. It is worth mentioning that, initially the tank was filled by water and this is useful for describing the upper boundary condition as pressure outlet with avoiding open surface boundary condition. So, this simplification will help in the CFD simulation such that the upper boundary will be considered as outflow boundary condition. In Fig.1 a systematic diagram for the problem is drawn. Red dye is applied to visualize the charging process of the tank. In this experiment the density of the water coloring is slightly different from the water density. So, the amount of change of density will be negligibly small. The colored hot water enters the tank through a horizontal inlet nozzle that has a diameter of 0.007 m. The sketch in Fig.1 indicates that the entering pipe has a converging nozzle, i.e., a nozzle diameter less than the pipe itself. The desired Reynolds number (Re) is determined by flow rate and temperature of the nozzle. In this study we need different flow rates high or low. So, the choice of this design of the nozzle can provide a quite large flow rate that may be not obtained by the normal pipe.

    Table 1 Overview of the experiments

    Fig.2 Visualization of flow structures for ΔT=0 K (1st row), ΔT=20 K (2nd row) and ΔT=40 K (3rd row)

    The parameters used for the three cases of studies are listed in Table 1, where Re is the Reynolds num-ber,Tcoldthe temperature of the water in the tank at the beginning of the experiment,Tinthe temperature of the entering water, ΔT=Tin-Tcoldthe temperature difference between the hot water and the cold water andtthe total time of the visualization. For the three cases, the flow rate is 0.7 l/min and the inlet velocity is 0.3 m/s.

    Visualizations are performed by using a digital video camera every five seconds. Figure 2 shows an example of visualization for the flow patterns during a draw-off test with an inlet flow rate of 0.7 l/min with three different inlet temperatures for the three cases. The buoyancy effects due to natural convection are clearly illustrated in Fig.2, such that for the isothermal case (ΔT=0 K), there is no buoyancy effect because there is no heat transfer, while for the second and third cases it can be observed that the red colored water (hot water) moves to the top of the tank depending on the temperature difference. Also, one can note that the buoyancy force of Case 3 ΔT=40 K is greater than that of Case 2 ΔT=20 K.

    2. Mathematical formulations

    The problem of a 3-D vertical water jet entering water in a rectangular storage tank is considered. The RANS with the realizablekε- model are used for the modeling of the turbulent flow. The realizablekε- model included an eddy-viscosity formula and a model equation for dissipationεbased on the dynamic equation of the mean-square vorticity fluctuation. So, the governing equations of mass, momentum, energy and turbulence take the form:

    Continuity equation

    Turbulence kinetic energy ()kequation

    whereSis the modulus of the tensor of the mean strain rate which is defined by

    is the generation of turbulence due to buoyancy.Prtis the turbulent Prandtl number for energy andβis the thermal expansion coefficient. The model constants of thek-εmodel are established to ensure that the model performs well for certain canonical flows.C1ε=1.44,C2=1.9,σk=1.0,Prt=0.85 andσε=1.2.

    The turbulence intensity can be defined as the ratio of the magnitude of the root mean square (rms) turbulent fluctuations to the reference velocity

    wherekis the turbulence kinetic energy andVis the reference velocity. The reference value specified should be the mean velocity magnitude for the flow. The velocity magnitude can be defined as the magnitude of the rms velocity components. Note that turbulence intensity can be defined in different ways. It is clear that the unit quantity ofVis velocity whileIis a dimensionless quantity, becausekhas a velocity square unit.

    Fig.3 Model outline of the symmetrical half tank

    3. Numerical investigations

    The current problem is to visualize and to simulate the stratification of the hot water inside a water store. Symmetry is assumed in the vertical central plane of the tank, so for this model only half of the tank is modeled. Figure 3 shows the model outline for the symmetrical half tank with the dimensionsX×Y×Z=0.486 m×0.1375 m×0.29 m. Furthermore, the heat transfer through the walls of tank is not taken into consideration, i.e., the tank walls are adiabatic. These simplifications are considered to reduce the total number of grid elements in the numerical mesh. The meshes are built up of hexahedral map cells, and the number of grid elements used in the model is 50,975 cells for the “fine grid” or “grid-1” which was used in all calculations. Also, two other grid cases are used to compare the results with the fine grid, one of them is a coarse grid which has 15 310 cells while the other is fine also (named “grid-2”) and has 46 715 cells. Figures 4(a) and 4(b) illustrate the fine grid faces of theXZ- plane and theXY- plane, respectively. For the coarse grid, faces of theXZ- plane and theXY-plane are shown in Figs.5(a) and 5(b), respectively. The CFD code Fluent with the grid generation tool Gambit is used to model the flow in the tank by solving the momentum, turbulence and energy equations. The realizablek-εmodel was used for modeling of the turbulent flow and the buoyancy is modeled using the Boussinesq approximation.

    Fig.4 Fine grid

    Because the water inside the tank was at rest before starting the water injection, the initial values of all dependent variables are taken to be zero except the initial temperature which is taken to be 293.15oK.

    Fig.5 Coarse grid

    Numerical solutions are obtained for unsteady flow while pressure, velocity, temperature and turbulence distributions inside the water tank are analyzed. In order to achieve convergence, the under-relaxation parameters are used on pressure, velocities, energy, turbulent viscosity, turbulence kinetic energy and turbulent dissipation rate.

    Body force weighted discretization was used for pressure and the velocity-pressure coupling was treated using the PISO algorithm with the Skewness-Neighbor coupling correction parameters. A secondorder upwind scheme was used in the equations of momentum, energy, turbulence kinetic energy and turbulence dissipation rate. A 3-D segregated implicit solver with the implicit second-order scheme is used for unsteady formulations.

    Fig.6 Residuals of the dependent variables for the first 1 000 iterations at =t1.0 s for Case 1

    A time step of 0.1 s was used during the calculations. The heat losses from the tank have not been modeled. The number of iterations used in this problem is 100 iterations for each time step. That was enough to achieve the convergence for all physical quantities for the three different cases. Figure 6 shows the residuals of all dependent variables for the first 1 000 iterations (t=1s) for Case 1.

    In this study, we will define two lines, and one of them is parallel to thex-axis with the dimensions (x=0-0.36,y=0.0001,z=0.1), this line lies in the core of the jet and will denoted byx. The other is parallel to thez-axis with the dimensions (x=0.1,y=0.0001,z=0-0.29) and is close to the wall opposite the jet entrance, which is also an important area.

    Fig.7 Temperature (a) and turbulence intensity (b) againstzwith different grid size at =t105 s, of Case 3

    Figure 7 shows a comparison for the temperature and turbulence intensity, plotted as a function ofz, between the two different fine grid sizes “grid-1” and“grid-2” as defined above. These figures indicate that the solutions are not depending on the mesh for these cases. Also, Figs.8(a) and 8(b) illustrate a comparison between the fine grid size case and the coarse grid case for the static pressure and the temperature, respectively. From these figures it can be seen that the solutions for the two cases have the same behavior with only small quantitative differences.

    Fig.8 Static pressure (a) temperature profiles (b) as a function of x, with different grid size at =t75 s, of Case 3

    Fig.9 Static pressure distribution as a function of x with various values of time for Case 1

    4. Results and discussion

    Static pressure along of the defined straight linexwith various values of the time for Case 1, ΔT=0 K, is shown in Fig.9. This figure indicates that the static pressure decreases as time increases alongx. Also, it can be seen that the static pressure in the jet area is greater than it in close to the wall.

    The area between the inlet and the opposite wall has a small static pressure. It is interesting to note from Fig.9 that with advancing time steps the effect of time on pressure tends to vanish. This means that unsteady state tends to approach steady state.

    Fig.10 Turbulence intensity profiles as a function ofz, with various values of time for Case 1

    Fig.11 Velocity magnitude profiles as a function ofzwith various values of time for Case 2

    Figure 10 shows the turbulence intensity profiles as a function ofzwith various values of time for Case 1. It is observed that the highest turbulence intensity lies in the jet area. The top and the down parts of the tank have small turbulence intensity. Time effect on the turbulence intensity is non-uniformly depending on the turbulence state inside the tank. The distribution of the velocity magnitude as a function ofzwith various values of time for Case 2 is shown in Fig.11. It is noteworthy that the middle part (jet area) of the tank has the highest velocity magnitude while it is smaller in the top part and tends to vanish in the bottom part of the tank.

    Fig.12 Temperature profiles as a function of z with various values of time for Case 2

    Figure 12 shows the temperature distribution as a function ofzwith various values of time for Case 2. It is clear that the top part of the tank has the highest temperature and the opposite is true for the bottom part where the temperature seems to be the same as the initial temperature of the tank (293 K). From the same figure it can be seen that the temperature increases as time increases.

    Fig.13 Static pressure profiles as a function of x with temperature differences =TΔ0 K,20 K,40 K at t=85 s

    Fig.14 Velocity Magnitude profiles as a function of (a) x, (b) z with different temperature differences at t=85 s

    A comparison between the static pressure profiles for the three cases of temperature differences, ΔT=0 K, 20 K, 40 K as a function ofxare shown in Fig.13 att=85 s. Figure 13 shows that the static pressure profiles along the linexdecrease as the temperature difference increases, and that the highest pressure for each case lies in the inlet area.

    In Figs.14(a) and 14(b) the velocity magnitude profiles are plotted as a function ofxandz, respectively plotted for different temperature differences ΔT=0 K, 20 K, 40 K att=85 s. Figure 14(a) indicates that the velocity magnitude alongxclose to the wall for Case 1 is higher than for the second and third cases. This because the linexlies at the center of the jet for Case 1, while for the second and third cases the position of the jet area lies abovexaccording to buoyancy effect. From Fig.14(b), it is noteworthy that the position of the highest velocity magnitude depends on the position of the jet area that is shifted by buoyancy force in terms of temperature difference.

    The temperature profiles as a function ofzare shown in Fig.15 for three temperature differences ΔT=0 K, 20 K, 40 K att=85 s. From Fig.15 it can be seen that the temperature is increased by raising the temperature difference especially in the top of the tank, while there is no difference at the tank bottom.

    Fig.16 Turbulence intensity profiles as a function of z for different temperature differences at t=85 s

    Figure 16 shows the turbulence intensity profiles as a function ofzfor different temperature differences ΔT=0 K, 20 K, 40 K att=85 s. From Fig.16 one can see that the position of the highest turbulence intensity depends on the temperature difference. The highest turbulence intensity lies in the jet area, thus forCase 1, ΔT=0 K, in the position around the inlet (because in this case there is no buoyancy force). The position of the highest turbulence intensity of Case 2 ΔT=20 K lies above the position of Case 1, and the position of the highest turbulence intensity of Case 3 ΔT=40 K lies above the position of Case 2.

    Fig.17 Velocity vectors in parallel vertical XY-plane x=0.1, x=0.2, x=0.3 and x=0.4 for Case 3, at 85 s

    Figure 17 shows the velocity vectors on the parallel verticalXYplanesx=0.1,x=0.2,x=0.3 andx=0.4 for Case 3 at 85 s. It is obvious that, the jet area has the highest velocity. After impacting the opposite wall the water returns above and below the jet.

    Fig.18 Internal Froude number as a function of z for Case 2

    P hysical pro perties o f water such a s den sity, dynamicviscosity,thermalconductivityandspecific heat are expressed as functions of temperature using the interpolation method in the temperature range 293.15 K-313.15 K for Case 2:

    Fig.19 Comparison between Boussinesq and polynomial approximation at t=85 s for Case 2

    The purpose of this section is to compare between this approximation and the Boussinesq approximation which is used above as shown in Fig.19. The results indicate that the distribution of velocity and temperature have the same behavior with small quantitative differences for the two approximations, while the behavior of pressure distribution as a function ofzis different (see Fig.19(a)). This is because in the Boussinesq approximation density varies only in the buoyancy force term in momentum equation and is considered constant in the other terms, while in the polynomial approximation the density varies always as a function of the temperature. Since the pressure term is divided by density, it should be affected by variation of density, which leads to a more realistic pressure distribution.

    Fig.20 Contours of the temperature, for times 10 s, 1 min and 2 min, for Case 3

    A good prediction of the stratification phenomena is shown in Fig.20. Contours of the temperature, at times 10 s, 1 min, 2 min of injection for Case 3 are plotted. It is observed from this figure that the layers of temperature appear early after a few seconds and become more clearly after a while. Moreover, this figure indicates that the realizablekε- model can model this problem very well.

    Figure 21 shows pictures of the experimental visualization, and the contours of the velocity magnitude through a vertical section (the symmetry plane), the vertical section =x0.1, the horizontal sections (=z0.06, =z0.1 and =z0.16 for Case 1, =z0.1, =z0.16 and =z0.27 for Case 3), at 65 s for Case 1 and at 75 s for Case 3. These figures provide a complete understanding of jet flow under the current conditions. It is noteworthy that the non-buoyant jet is a pure jet, or in other words it is momentum-dominated jet. However, the buoyancy effect takes place when we have a temperature difference and the jet takes a plume (or forced plume) shape in which both momentum and buoyancy have effect.

    Fig.21 Experimental visualization, and contours of the velocity magnitude through vertical and horizontal sections

    5. Conclusions

    CFD treatments and experimental visualizations for the problem of a hot water jet entering horizontally into a water store containing cold water have been performed to illustrate the varying behavior of the thermal conditions in a solar store. Three temperature differences with the corresponding Reynolds numbers are considered. From this work one can conclude that:

    (1) When these problems are treated as unsteadystate problems, after certain time the unsteady-state develops into the steady state.

    (2) Jet entering the tank at the same temperature results in a regular spread of the jet inside the tank, while entering of a hot jet into the tank make the spread of the jet inside the tank depend on buoyancy forces due to the temperature difference.

    (3) Thekε- turbulence model is suitable to decribe turbulent flow for this kind of problems, the Boussinesq approximation was used to model the buoyancy effect.

    (4) Thermal stratification is verified through the temperature distribution.

    (5) A polynomial approximation for the waterproperties as functions of temperature was compared to the results with the Boussinesq approximation.

    [1] KNUDSEN S., FURBO S. and SHAH L. J. Design of the inlet to the mantle in a vertical mantle storage tank[C]. Proceedings ISES 2001 Solar World Congress. Adelaide, Australia, 2001.

    [2] SHAH L. J., FURBO S. Entrance effects in solar storage tanks[J]. Solar Energy, 2003, 75(4): 337-348.

    [3] JORDAN U., FURBO S. Investigations of the flow into a storage tank by means of advanced experimental and theoretical methods[C]. ISES Solar World Congress. Goteborg, Sweden, 2003.

    [4] KNUDSEN S., MORRISON G. L. and BEHNIA M. et al. Analysis of the flow structure and heat transfer in a vertical mantle heat exchanger[C]. ISES Solar World Congress. Goteborg, Sweden, 2003.

    [5] El-AMIN M. F., HEIDEMANN W. and MüLLERSTEINHAGEN H. Turbulent jet flow into a water store[C]. Proceedings Heat Transfer in Components and Systems for Sustainable Energy Technologies. Grenoble, France, 2005, 345-350.

    [6] El-AMIN M. F., SUN S. and HEIDEMANN W. et al. Analysis of a turbulent buoyant confined jet modeled using realizable model[J]. Heat Mass Transfer, 2010, 46(8): 943-960.

    [7] PANTHALOOKARAN V., El-AMIN M. F. and HEIDEMANN W. et al. Calibrated models for simulation of stratified hot water heat stores[J]. International Journal of Energy Research, 2008, 32(7): 661-676.

    [8] FUKUSHIMA C., AANEN L. and WESTERWEEL J. Investigation of the mixing process in an axisymmetric turbulent jet using PIV and LIF[C]. Proceedings 10th International Symposium Applications of Laser Techniques Fluid Mechanics, Lisbon, Portugal, 2000.

    [9] AGRAWAL A., PRASAD A. K. Integral solution for the mean flow profiles of turbulent jets, plumes, and wakes[J]. Journal of Fluids Engineering, 2003, 125(5): 813-822.

    [10] ARAKERI J. H., DAS D. and SRINIVASAN J. Bifurcation in a buoyant horizontal laminar jet[J]. Journal of Fluid Mechanics, 2000, 412: 61-73.

    [11] O’HERN T. J., WECKMAN E. J.and GERHART A. L. et al. Experimental study of a turbulent buoyant helium plume[J]. Journal of Fluid Mechanics, 2005, 544: 143-171.

    [12] El-AMIN M. F., KANAYAMA H. Integral solutions for selected turbulent quantities of small-scale hydrogen leakage: A non-buoyant jet or momentum-dominated buoyant jet regime[J]. International Journal of Hydrogen Energy, 2009, 34(3): 1607-1612.

    [13] El-AMIN M. F., KANAYAMA H. Similarity consideration of the buoyant jet resulting from hydrogen leakage[J]. International Journal of Hydrogen Energy, 2009, 34(14): 5803-5809.

    [14] El-AMIN M. F. Non-Boussinesq turbulent buoyant jet resulting from hydrogen leakage in air[J]. Internatinal Journal of Hydrogen Energy, 2009, 34(3): 7873-7882

    [15] JIRKA G. H. Integral model for turbulent buoyant jets in unbounded stratified flows, Part 1: Single round jet[J]. Environmental Fluid Mechanics, 2004, 4(1): 1-56.

    [16] JIRKA G. H. Integral model for turbulent buoyant jets in unbounded stratified flows, Part 2: Plane jet dynamics resulting from multiport diffuser jets[J]. Environmentai Fluid Mechanics, 2006, 6(1): 43-100.

    10.1016/S1001-6058(14)60012-3

    * Biography: EL-AMIN M. F. (1971-), Male, Ph. D., Research Scientist

    男人舔女人的私密视频| 国产主播在线观看一区二区| 在线观看66精品国产| 久久国产亚洲av麻豆专区| 后天国语完整版免费观看| 欧美黑人欧美精品刺激| 欧美精品av麻豆av| 国产精品美女特级片免费视频播放器 | 成人国语在线视频| 欧美日韩成人在线一区二区| 精品视频人人做人人爽| 久久久精品区二区三区| 精品高清国产在线一区| 18禁裸乳无遮挡动漫免费视频| 电影成人av| 久久中文字幕人妻熟女| 国产精品久久久人人做人人爽| 国产亚洲一区二区精品| 少妇猛男粗大的猛烈进出视频| 夫妻午夜视频| 99热国产这里只有精品6| 久久天躁狠狠躁夜夜2o2o| 五月开心婷婷网| 在线观看免费视频日本深夜| 成熟少妇高潮喷水视频| 美女高潮到喷水免费观看| 久久狼人影院| 国产日韩欧美亚洲二区| 亚洲午夜精品一区,二区,三区| 啦啦啦在线免费观看视频4| 欧美黄色淫秽网站| 亚洲人成电影免费在线| av网站免费在线观看视频| 久久精品人人爽人人爽视色| avwww免费| 极品少妇高潮喷水抽搐| 激情视频va一区二区三区| 国产精品国产高清国产av | 欧美午夜高清在线| 日韩制服丝袜自拍偷拍| 老司机亚洲免费影院| 少妇猛男粗大的猛烈进出视频| 精品一区二区三区av网在线观看| 成人18禁在线播放| 久久久久视频综合| 中国美女看黄片| 免费观看人在逋| 亚洲熟妇熟女久久| 日本wwww免费看| 亚洲一码二码三码区别大吗| 女人久久www免费人成看片| 一进一出好大好爽视频| 无遮挡黄片免费观看| 老熟妇仑乱视频hdxx| 久久久精品区二区三区| www日本在线高清视频| 精品第一国产精品| 亚洲综合色网址| 99久久国产精品久久久| 亚洲av日韩在线播放| av有码第一页| 欧美黄色片欧美黄色片| 国产主播在线观看一区二区| 久久久精品免费免费高清| 亚洲第一av免费看| 99国产极品粉嫩在线观看| 在线观看www视频免费| 免费一级毛片在线播放高清视频 | 黄网站色视频无遮挡免费观看| av不卡在线播放| 免费一级毛片在线播放高清视频 | 久久久久久久国产电影| 成年人午夜在线观看视频| 亚洲中文日韩欧美视频| 久久青草综合色| 午夜成年电影在线免费观看| 欧美丝袜亚洲另类 | 水蜜桃什么品种好| 丝袜美足系列| 正在播放国产对白刺激| 两人在一起打扑克的视频| 99精品久久久久人妻精品| 夜夜爽天天搞| 18禁裸乳无遮挡动漫免费视频| 99riav亚洲国产免费| 国产97色在线日韩免费| 在线观看免费日韩欧美大片| 国产一卡二卡三卡精品| 不卡av一区二区三区| 啦啦啦在线免费观看视频4| 国产精品永久免费网站| 久久九九热精品免费| 天天操日日干夜夜撸| 黑丝袜美女国产一区| 精品人妻熟女毛片av久久网站| 国产高清激情床上av| av天堂在线播放| 丰满饥渴人妻一区二区三| 老司机亚洲免费影院| 国产成人精品久久二区二区免费| 男女高潮啪啪啪动态图| 成年人黄色毛片网站| 亚洲成av片中文字幕在线观看| 国产无遮挡羞羞视频在线观看| 男男h啪啪无遮挡| 精品久久久久久久久久免费视频 | 久久精品国产综合久久久| 国产麻豆69| 麻豆av在线久日| 免费人成视频x8x8入口观看| 国产免费现黄频在线看| 另类亚洲欧美激情| 在线视频色国产色| 免费在线观看视频国产中文字幕亚洲| 欧美中文综合在线视频| 亚洲一区高清亚洲精品| 久久久精品国产亚洲av高清涩受| 亚洲人成电影观看| 久久天躁狠狠躁夜夜2o2o| 大陆偷拍与自拍| videosex国产| 亚洲成人国产一区在线观看| 两个人免费观看高清视频| 亚洲精品在线观看二区| 午夜精品国产一区二区电影| 丰满饥渴人妻一区二区三| 亚洲欧美激情在线| 午夜激情av网站| 欧美性长视频在线观看| 黑人操中国人逼视频| 亚洲av成人一区二区三| 窝窝影院91人妻| 国产成人精品久久二区二区免费| 99精品在免费线老司机午夜| 我的亚洲天堂| 免费观看精品视频网站| 黄片播放在线免费| 91老司机精品| 久久热在线av| 精品国产乱子伦一区二区三区| 久久国产精品人妻蜜桃| 最新美女视频免费是黄的| 制服人妻中文乱码| 亚洲在线自拍视频| 夜夜躁狠狠躁天天躁| 看免费av毛片| 99国产精品99久久久久| 色综合婷婷激情| 色播在线永久视频| 亚洲专区字幕在线| 亚洲av成人一区二区三| 王馨瑶露胸无遮挡在线观看| 一本综合久久免费| 超碰成人久久| 热99国产精品久久久久久7| 欧美激情久久久久久爽电影 | 一边摸一边抽搐一进一小说 | 中文字幕最新亚洲高清| 久久久国产成人精品二区 | 大香蕉久久成人网| 日本a在线网址| 免费观看人在逋| 免费观看人在逋| 少妇裸体淫交视频免费看高清 | 久久国产亚洲av麻豆专区| 亚洲在线自拍视频| 国产精品永久免费网站| 国产精品永久免费网站| 欧洲精品卡2卡3卡4卡5卡区| 熟女少妇亚洲综合色aaa.| 99久久国产精品久久久| 国产97色在线日韩免费| 国产aⅴ精品一区二区三区波| 狂野欧美激情性xxxx| bbb黄色大片| 男人的好看免费观看在线视频 | 亚洲色图综合在线观看| 久久人人爽av亚洲精品天堂| 国产激情欧美一区二区| 亚洲色图综合在线观看| 免费一级毛片在线播放高清视频 | www.自偷自拍.com| 国产深夜福利视频在线观看| 99国产极品粉嫩在线观看| 新久久久久国产一级毛片| 男女午夜视频在线观看| 精品亚洲成国产av| 99re在线观看精品视频| 国产一卡二卡三卡精品| 精品人妻熟女毛片av久久网站| 国产亚洲欧美在线一区二区| 人人妻人人爽人人添夜夜欢视频| 男女床上黄色一级片免费看| www.999成人在线观看| 免费高清在线观看日韩| videosex国产| 国产精品自产拍在线观看55亚洲 | 老司机福利观看| 亚洲avbb在线观看| 黄片播放在线免费| 亚洲中文日韩欧美视频| 757午夜福利合集在线观看| 变态另类成人亚洲欧美熟女 | 国产有黄有色有爽视频| 男女下面插进去视频免费观看| 亚洲国产欧美一区二区综合| a级毛片在线看网站| 青草久久国产| 精品无人区乱码1区二区| 亚洲国产中文字幕在线视频| 精品国产乱子伦一区二区三区| 色综合婷婷激情| 9热在线视频观看99| 久久香蕉精品热| 亚洲专区中文字幕在线| 丁香欧美五月| 热99久久久久精品小说推荐| 人人妻人人爽人人添夜夜欢视频| 亚洲熟妇熟女久久| 日韩欧美免费精品| 99riav亚洲国产免费| 国产麻豆69| 多毛熟女@视频| 99国产精品99久久久久| 亚洲色图av天堂| 欧美 日韩 精品 国产| 午夜福利乱码中文字幕| 精品国产美女av久久久久小说| 老汉色av国产亚洲站长工具| 黄色视频不卡| 一夜夜www| 欧美亚洲 丝袜 人妻 在线| 在线观看日韩欧美| 精品久久蜜臀av无| 久久影院123| 一级毛片精品| 久久精品亚洲熟妇少妇任你| av免费在线观看网站| 动漫黄色视频在线观看| 99久久人妻综合| 一区二区日韩欧美中文字幕| 亚洲成人国产一区在线观看| 视频区欧美日本亚洲| 亚洲av第一区精品v没综合| 亚洲五月婷婷丁香| 女人精品久久久久毛片| 在线播放国产精品三级| 91老司机精品| 亚洲自偷自拍图片 自拍| 久热爱精品视频在线9| 久久久久久久国产电影| 欧美精品高潮呻吟av久久| 午夜免费观看网址| 久久国产精品大桥未久av| 一级片免费观看大全| 亚洲免费av在线视频| 国产高清国产精品国产三级| 国产成+人综合+亚洲专区| 亚洲国产欧美日韩在线播放| 日韩 欧美 亚洲 中文字幕| 女人久久www免费人成看片| 亚洲一区二区三区欧美精品| 女人被狂操c到高潮| 青草久久国产| 后天国语完整版免费观看| 大片电影免费在线观看免费| 美女福利国产在线| 热re99久久精品国产66热6| av欧美777| 极品教师在线免费播放| 亚洲成国产人片在线观看| 一级片'在线观看视频| svipshipincom国产片| 午夜福利影视在线免费观看| 亚洲五月色婷婷综合| 日本欧美视频一区| 午夜福利在线免费观看网站| 欧美日韩视频精品一区| 国产成人啪精品午夜网站| 在线观看舔阴道视频| 悠悠久久av| a级片在线免费高清观看视频| 日韩熟女老妇一区二区性免费视频| 91精品三级在线观看| 午夜免费观看网址| 激情在线观看视频在线高清 | 这个男人来自地球电影免费观看| 热re99久久精品国产66热6| 色婷婷久久久亚洲欧美| 免费看a级黄色片| 久久香蕉国产精品| 欧美日韩成人在线一区二区| 大码成人一级视频| 精品少妇久久久久久888优播| 一级a爱视频在线免费观看| 亚洲国产毛片av蜜桃av| 日韩熟女老妇一区二区性免费视频| 国产精品久久久久久人妻精品电影| 夜夜夜夜夜久久久久| xxxhd国产人妻xxx| 国产亚洲av高清不卡| av在线播放免费不卡| 午夜福利在线免费观看网站| 制服诱惑二区| 欧美另类亚洲清纯唯美| 国产在视频线精品| 三上悠亚av全集在线观看| 国产精品免费一区二区三区在线 | 别揉我奶头~嗯~啊~动态视频| 丰满人妻熟妇乱又伦精品不卡| 又黄又粗又硬又大视频| 999久久久国产精品视频| 男女床上黄色一级片免费看| 亚洲一区二区三区欧美精品| 欧美成人免费av一区二区三区 | 丰满饥渴人妻一区二区三| avwww免费| 亚洲欧美一区二区三区黑人| 日本vs欧美在线观看视频| svipshipincom国产片| 黑人操中国人逼视频| 亚洲国产欧美日韩在线播放| 少妇被粗大的猛进出69影院| 交换朋友夫妻互换小说| 最新在线观看一区二区三区| 少妇猛男粗大的猛烈进出视频| 亚洲三区欧美一区| 精品国产一区二区三区四区第35| 午夜福利免费观看在线| 精品国产一区二区久久| 精品一区二区三区av网在线观看| 一二三四在线观看免费中文在| 国产淫语在线视频| 精品国产一区二区三区久久久樱花| 一区在线观看完整版| 国产激情久久老熟女| 欧美乱色亚洲激情| 欧美成人午夜精品| 国产淫语在线视频| 久久中文字幕人妻熟女| 日本撒尿小便嘘嘘汇集6| 丝袜美足系列| 欧美+亚洲+日韩+国产| 美女高潮到喷水免费观看| 国产精华一区二区三区| av天堂久久9| 精品少妇一区二区三区视频日本电影| 国产一区在线观看成人免费| 99国产综合亚洲精品| 精品久久久久久电影网| 亚洲中文日韩欧美视频| 男女午夜视频在线观看| 丝袜美腿诱惑在线| 久久精品国产99精品国产亚洲性色 | 精品久久蜜臀av无| 欧美 亚洲 国产 日韩一| 色在线成人网| 91av网站免费观看| 婷婷丁香在线五月| 国产精品乱码一区二三区的特点 | 香蕉国产在线看| 欧美大码av| 亚洲熟女毛片儿| 久久久久精品人妻al黑| 欧美日韩福利视频一区二区| 国精品久久久久久国模美| 三上悠亚av全集在线观看| 亚洲av日韩精品久久久久久密| 别揉我奶头~嗯~啊~动态视频| 日日爽夜夜爽网站| 国产精品欧美亚洲77777| 欧美不卡视频在线免费观看 | 日本vs欧美在线观看视频| 咕卡用的链子| 亚洲第一av免费看| 亚洲国产精品一区二区三区在线| 国产男女内射视频| 身体一侧抽搐| 黄片小视频在线播放| 麻豆国产av国片精品| 精品一品国产午夜福利视频| 黄色视频不卡| 国产精品成人在线| 啦啦啦在线免费观看视频4| 韩国精品一区二区三区| 99国产综合亚洲精品| 亚洲第一欧美日韩一区二区三区| 老司机亚洲免费影院| 他把我摸到了高潮在线观看| 午夜精品在线福利| 精品高清国产在线一区| 久久久国产欧美日韩av| 久久久久久久久久久久大奶| 搡老乐熟女国产| 久久久久久免费高清国产稀缺| 在线免费观看的www视频| 99热国产这里只有精品6| 91在线观看av| 久久狼人影院| 嫁个100分男人电影在线观看| 精品一区二区三区av网在线观看| 一进一出抽搐gif免费好疼 | 精品久久久久久,| av网站免费在线观看视频| 黑人欧美特级aaaaaa片| www.自偷自拍.com| 狠狠狠狠99中文字幕| 最新美女视频免费是黄的| 女人被躁到高潮嗷嗷叫费观| 露出奶头的视频| 精品午夜福利视频在线观看一区| 久久国产精品影院| 啦啦啦 在线观看视频| 亚洲精品乱久久久久久| 欧美激情久久久久久爽电影 | 久久ye,这里只有精品| 国产精品综合久久久久久久免费 | 99久久国产精品久久久| 99re在线观看精品视频| 黄色视频,在线免费观看| 99久久综合精品五月天人人| 一级a爱视频在线免费观看| 天堂√8在线中文| 久久久久久亚洲精品国产蜜桃av| 老司机在亚洲福利影院| 高清视频免费观看一区二区| 无人区码免费观看不卡| 亚洲国产精品合色在线| 人妻一区二区av| 欧美中文综合在线视频| 国产99白浆流出| videosex国产| 多毛熟女@视频| 女人被躁到高潮嗷嗷叫费观| 国产精品自产拍在线观看55亚洲 | 99国产极品粉嫩在线观看| 国产精品久久久久久人妻精品电影| 久久人人爽av亚洲精品天堂| 自线自在国产av| av网站免费在线观看视频| 看免费av毛片| 国产精品久久久av美女十八| 高清毛片免费观看视频网站 | 国产国语露脸激情在线看| 不卡av一区二区三区| 在线永久观看黄色视频| bbb黄色大片| av免费在线观看网站| 国产免费现黄频在线看| 欧美中文综合在线视频| 精品电影一区二区在线| av超薄肉色丝袜交足视频| 热re99久久国产66热| 另类亚洲欧美激情| 操美女的视频在线观看| 中文亚洲av片在线观看爽 | 中文字幕av电影在线播放| av网站在线播放免费| 国产乱人伦免费视频| 午夜福利免费观看在线| 国产精品永久免费网站| 老熟女久久久| 欧美亚洲 丝袜 人妻 在线| 香蕉久久夜色| 午夜福利,免费看| 亚洲第一欧美日韩一区二区三区| 黄色 视频免费看| 亚洲黑人精品在线| 国产激情欧美一区二区| 少妇被粗大的猛进出69影院| 操美女的视频在线观看| 欧美精品亚洲一区二区| 国产av精品麻豆| 亚洲精品国产色婷婷电影| 成人av一区二区三区在线看| 90打野战视频偷拍视频| 女人被狂操c到高潮| 黄色毛片三级朝国网站| 麻豆乱淫一区二区| 亚洲精品久久成人aⅴ小说| 国产成人精品久久二区二区91| 精品一区二区三区四区五区乱码| 日日夜夜操网爽| av在线播放免费不卡| 欧美精品一区二区免费开放| 国产免费男女视频| 美女高潮喷水抽搐中文字幕| 欧美成人免费av一区二区三区 | 两个人免费观看高清视频| 欧美日韩一级在线毛片| 水蜜桃什么品种好| 成在线人永久免费视频| 国产人伦9x9x在线观看| 国产一区有黄有色的免费视频| 久久久国产成人免费| 国产精品98久久久久久宅男小说| 国产精品自产拍在线观看55亚洲 | 久久亚洲精品不卡| 国产麻豆69| 久久久久久免费高清国产稀缺| 最近最新中文字幕大全免费视频| 一级毛片精品| 18禁国产床啪视频网站| 亚洲精品乱久久久久久| 欧美亚洲日本最大视频资源| 色精品久久人妻99蜜桃| 国产精品久久久久成人av| 国产欧美亚洲国产| 国产片内射在线| 色在线成人网| 国产欧美日韩综合在线一区二区| 国产精品1区2区在线观看. | 麻豆国产av国片精品| 手机成人av网站| videosex国产| 91老司机精品| 国产精品久久久人人做人人爽| 国产成人一区二区三区免费视频网站| 亚洲五月色婷婷综合| 日本a在线网址| 亚洲av成人一区二区三| 亚洲av美国av| svipshipincom国产片| 国产精品久久久av美女十八| 精品无人区乱码1区二区| 黄色 视频免费看| 免费看十八禁软件| 一二三四在线观看免费中文在| 午夜亚洲福利在线播放| 色综合婷婷激情| 91精品三级在线观看| 国产人伦9x9x在线观看| 亚洲第一av免费看| 性少妇av在线| 国产高清视频在线播放一区| 电影成人av| 9热在线视频观看99| 欧美日韩亚洲综合一区二区三区_| 亚洲色图综合在线观看| 人人妻人人爽人人添夜夜欢视频| 成人18禁在线播放| 亚洲国产看品久久| 最新在线观看一区二区三区| 18禁观看日本| 国产熟女午夜一区二区三区| 精品少妇一区二区三区视频日本电影| 一夜夜www| 欧美人与性动交α欧美软件| 一个人免费在线观看的高清视频| 淫妇啪啪啪对白视频| 99精品欧美一区二区三区四区| 丰满人妻熟妇乱又伦精品不卡| 亚洲成人免费av在线播放| 免费在线观看亚洲国产| 18禁观看日本| av福利片在线| 精品国产亚洲在线| 国产色视频综合| 亚洲中文字幕日韩| 午夜福利在线观看吧| 99国产精品一区二区三区| 国产一区二区激情短视频| 精品国内亚洲2022精品成人 | 亚洲av日韩精品久久久久久密| 国产成人av激情在线播放| avwww免费| 国产激情久久老熟女| 国产真人三级小视频在线观看| 免费在线观看亚洲国产| 午夜视频精品福利| 亚洲人成电影免费在线| 国产精品综合久久久久久久免费 | 午夜免费鲁丝| 18禁国产床啪视频网站| 久久人人爽av亚洲精品天堂| 国产91精品成人一区二区三区| 丰满人妻熟妇乱又伦精品不卡| 久久精品国产亚洲av香蕉五月 | 久久国产精品男人的天堂亚洲| 少妇的丰满在线观看| 亚洲欧美色中文字幕在线| 91麻豆av在线| 天天躁狠狠躁夜夜躁狠狠躁| 狂野欧美激情性xxxx| 国产真人三级小视频在线观看| 人妻久久中文字幕网| 国产不卡一卡二| 多毛熟女@视频| 欧美午夜高清在线| 国产在线精品亚洲第一网站| 深夜精品福利| 淫妇啪啪啪对白视频| 人人妻人人爽人人添夜夜欢视频| 我的亚洲天堂| 国产欧美亚洲国产| 午夜日韩欧美国产| 国产精品一区二区在线观看99| 韩国av一区二区三区四区| 香蕉丝袜av| 欧美激情久久久久久爽电影 | 欧美 亚洲 国产 日韩一| 狠狠狠狠99中文字幕| 97人妻天天添夜夜摸| 1024视频免费在线观看| 日韩中文字幕欧美一区二区| 欧美乱妇无乱码| 欧美日韩亚洲高清精品| 精品无人区乱码1区二区| 午夜影院日韩av| 亚洲综合色网址| 久久精品aⅴ一区二区三区四区| 国产高清videossex| 天堂动漫精品| 看黄色毛片网站| 亚洲九九香蕉|