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

    色av中文字幕| 国产精品 欧美亚洲| 男男h啪啪无遮挡| 久久国产精品影院| 18禁观看日本| 国产成人精品在线电影| a级毛片在线看网站| 最新在线观看一区二区三区| 国产成人精品久久二区二区免费| ponron亚洲| 深夜精品福利| 国产精品 国内视频| 亚洲成国产人片在线观看| 99精品在免费线老司机午夜| 欧美+亚洲+日韩+国产| 亚洲,欧美精品.| 两个人视频免费观看高清| 操美女的视频在线观看| 欧美乱色亚洲激情| 午夜福利欧美成人| 久久天躁狠狠躁夜夜2o2o| 精品午夜福利视频在线观看一区| 久久精品亚洲熟妇少妇任你| 亚洲片人在线观看| a在线观看视频网站| 色尼玛亚洲综合影院| 欧美国产精品va在线观看不卡| 亚洲精品一区av在线观看| av有码第一页| 很黄的视频免费| 男男h啪啪无遮挡| 午夜精品国产一区二区电影| 亚洲无线在线观看| 亚洲一区高清亚洲精品| 少妇熟女aⅴ在线视频| 村上凉子中文字幕在线| 搡老妇女老女人老熟妇| 成熟少妇高潮喷水视频| 91老司机精品| 精品少妇一区二区三区视频日本电影| 欧美在线黄色| 88av欧美| 人人妻人人澡人人看| 国产精品野战在线观看| 最好的美女福利视频网| 高清黄色对白视频在线免费看| 性少妇av在线| 女性被躁到高潮视频| 97人妻精品一区二区三区麻豆 | 国产成人精品无人区| 久久久久久亚洲精品国产蜜桃av| 超碰成人久久| 亚洲国产毛片av蜜桃av| 变态另类丝袜制服| 国产精华一区二区三区| 日日摸夜夜添夜夜添小说| 成人av一区二区三区在线看| av在线天堂中文字幕| 免费不卡黄色视频| 少妇熟女aⅴ在线视频| 亚洲午夜理论影院| 亚洲 国产 在线| 久久久久久大精品| 亚洲欧美精品综合久久99| 日韩国内少妇激情av| 日韩成人在线观看一区二区三区| 久久人妻熟女aⅴ| 老司机午夜福利在线观看视频| 老司机深夜福利视频在线观看| 欧美日本亚洲视频在线播放| 日本 欧美在线| 国产精品永久免费网站| 国产午夜福利久久久久久| 亚洲精品美女久久久久99蜜臀| 两性夫妻黄色片| 脱女人内裤的视频| 午夜影院日韩av| 黄频高清免费视频| 国产精品 欧美亚洲| 岛国在线观看网站| 国产免费av片在线观看野外av| 日本 欧美在线| 国产精品一区二区在线不卡| 精品电影一区二区在线| 日本a在线网址| 欧美另类亚洲清纯唯美| 色哟哟哟哟哟哟| 日韩三级视频一区二区三区| 性欧美人与动物交配| 久久精品亚洲熟妇少妇任你| 在线永久观看黄色视频| 99香蕉大伊视频| 久久 成人 亚洲| 久久久久亚洲av毛片大全| 亚洲九九香蕉| 丰满人妻熟妇乱又伦精品不卡| 一边摸一边做爽爽视频免费| 美女高潮到喷水免费观看| 99久久国产精品久久久| 成年版毛片免费区| 亚洲精品美女久久久久99蜜臀| 在线观看免费午夜福利视频| www日本在线高清视频| avwww免费| 亚洲精品久久国产高清桃花| 精品国产乱子伦一区二区三区| 18禁观看日本| 亚洲av熟女| svipshipincom国产片| 黑人巨大精品欧美一区二区mp4| 欧美激情 高清一区二区三区| 亚洲色图 男人天堂 中文字幕| 免费观看精品视频网站| 国产男靠女视频免费网站| 国产亚洲精品久久久久5区| 精品日产1卡2卡| 成年版毛片免费区| 男女床上黄色一级片免费看| 一进一出抽搐动态| 成年版毛片免费区| 国产一区二区在线av高清观看| 亚洲人成伊人成综合网2020| 亚洲av五月六月丁香网| 亚洲视频免费观看视频| 美女 人体艺术 gogo| 女人爽到高潮嗷嗷叫在线视频| 日本欧美视频一区| 久久香蕉激情| 亚洲va日本ⅴa欧美va伊人久久| 久久亚洲精品不卡| 精品熟女少妇八av免费久了| 搞女人的毛片| 精品人妻在线不人妻| 十八禁人妻一区二区| 看免费av毛片| 99re在线观看精品视频| 免费看十八禁软件| 人妻丰满熟妇av一区二区三区| 午夜福利高清视频| av超薄肉色丝袜交足视频| xxx96com| 亚洲美女黄片视频| 母亲3免费完整高清在线观看| 制服人妻中文乱码| 乱人伦中国视频| 久久精品国产亚洲av高清一级| www.熟女人妻精品国产| 午夜视频精品福利| 日日干狠狠操夜夜爽| 国产av一区在线观看免费| 黄色 视频免费看| 亚洲人成77777在线视频| 成熟少妇高潮喷水视频| 日本免费一区二区三区高清不卡 | 成人国产综合亚洲| 99在线人妻在线中文字幕| 国产成人一区二区三区免费视频网站| 久久精品亚洲熟妇少妇任你| 天堂√8在线中文| 精品一区二区三区四区五区乱码| 一边摸一边做爽爽视频免费| av在线天堂中文字幕| 午夜福利成人在线免费观看| 午夜两性在线视频| 国产一区二区三区在线臀色熟女| 色尼玛亚洲综合影院| 成人av一区二区三区在线看| 一个人观看的视频www高清免费观看 | 在线播放国产精品三级| 亚洲 欧美 日韩 在线 免费| 搡老妇女老女人老熟妇| 啦啦啦观看免费观看视频高清 | 人妻久久中文字幕网| 日韩中文字幕欧美一区二区| 亚洲七黄色美女视频| 国产高清videossex| 老熟妇仑乱视频hdxx| 亚洲精品中文字幕一二三四区| 欧美大码av| 自线自在国产av| 久久久久久久精品吃奶| 国产激情欧美一区二区| 两个人免费观看高清视频| 男人操女人黄网站| 国产精品98久久久久久宅男小说| 国产一区二区在线av高清观看| 亚洲第一欧美日韩一区二区三区| 日韩免费av在线播放| 亚洲七黄色美女视频| 亚洲在线自拍视频| 日韩欧美一区视频在线观看| 操美女的视频在线观看| 久久婷婷成人综合色麻豆| 免费观看人在逋| 最新美女视频免费是黄的| 97碰自拍视频| 久久性视频一级片| 老熟妇乱子伦视频在线观看| 亚洲av电影在线进入| 多毛熟女@视频| 一进一出抽搐gif免费好疼| 黄频高清免费视频| 国产1区2区3区精品| 国产真人三级小视频在线观看| 后天国语完整版免费观看| 欧美性长视频在线观看| 日日摸夜夜添夜夜添小说| 亚洲全国av大片| 啪啪无遮挡十八禁网站| 午夜福利影视在线免费观看| 长腿黑丝高跟| 午夜成年电影在线免费观看| 午夜影院日韩av| 禁无遮挡网站| 国语自产精品视频在线第100页| 国产精品九九99| 亚洲男人的天堂狠狠| 亚洲国产高清在线一区二区三 | 老鸭窝网址在线观看| 久久国产亚洲av麻豆专区| 男人舔女人下体高潮全视频| 免费不卡黄色视频| 午夜精品国产一区二区电影| 啪啪无遮挡十八禁网站| 一本大道久久a久久精品| 88av欧美| 正在播放国产对白刺激| 国产真人三级小视频在线观看| 欧美成狂野欧美在线观看| 久久精品亚洲熟妇少妇任你| 精品久久久久久久毛片微露脸| 国产精品一区二区在线不卡| av欧美777| x7x7x7水蜜桃| 中文亚洲av片在线观看爽| 香蕉丝袜av| 50天的宝宝边吃奶边哭怎么回事| 国产免费av片在线观看野外av| aaaaa片日本免费| 亚洲欧美激情在线| 亚洲国产精品成人综合色| 日日摸夜夜添夜夜添小说| 精品人妻在线不人妻| 日本vs欧美在线观看视频| 动漫黄色视频在线观看| 啪啪无遮挡十八禁网站| 天天一区二区日本电影三级 | 亚洲自偷自拍图片 自拍| 亚洲av成人av| 亚洲天堂国产精品一区在线| 男女下面进入的视频免费午夜 | 法律面前人人平等表现在哪些方面| 纯流量卡能插随身wifi吗| 久久精品91无色码中文字幕| 在线观看舔阴道视频| 99热只有精品国产| 亚洲欧美精品综合久久99| 亚洲 国产 在线| 一本久久中文字幕| 欧美性长视频在线观看| 精品国产乱码久久久久久男人| 国产成人av激情在线播放| 99re在线观看精品视频| www国产在线视频色| 日本免费一区二区三区高清不卡 | 亚洲伊人色综图| 午夜福利视频1000在线观看 | 香蕉国产在线看| 亚洲午夜精品一区,二区,三区| 女性生殖器流出的白浆| 欧美一级a爱片免费观看看 | 波多野结衣av一区二区av| 欧美色视频一区免费| 九色亚洲精品在线播放| 高清黄色对白视频在线免费看| 丝袜美足系列| 国产精品九九99| 久久影院123| 成人三级黄色视频| 性少妇av在线| 国产精品久久久久久亚洲av鲁大| 99在线视频只有这里精品首页| av免费在线观看网站| 欧美日韩乱码在线| 欧美乱码精品一区二区三区| 欧美日韩中文字幕国产精品一区二区三区 | 久久国产亚洲av麻豆专区| 国产精品免费一区二区三区在线| 日韩欧美在线二视频| 亚洲自拍偷在线| 可以在线观看的亚洲视频| 夜夜看夜夜爽夜夜摸| 久久亚洲真实| 又黄又粗又硬又大视频| 国产av精品麻豆| 久久 成人 亚洲| 欧美黑人欧美精品刺激| 久久人人爽av亚洲精品天堂| 黄片小视频在线播放| 一卡2卡三卡四卡精品乱码亚洲| av天堂在线播放| 在线观看免费午夜福利视频| 午夜视频精品福利| 欧美日韩一级在线毛片| 麻豆久久精品国产亚洲av| 久久精品91无色码中文字幕| 欧美日韩瑟瑟在线播放| 国产成人精品无人区| 国产伦一二天堂av在线观看| 久久久精品欧美日韩精品| www.熟女人妻精品国产| 色综合亚洲欧美另类图片| 成人亚洲精品av一区二区| 亚洲自拍偷在线| 女警被强在线播放| 身体一侧抽搐| 后天国语完整版免费观看| 亚洲欧美激情在线| a级毛片在线看网站| 岛国在线观看网站| 69av精品久久久久久| 免费在线观看视频国产中文字幕亚洲| 免费高清在线观看日韩| 欧美乱码精品一区二区三区| 久久人妻福利社区极品人妻图片| 搡老熟女国产l中国老女人| 女性生殖器流出的白浆| 亚洲中文日韩欧美视频| 亚洲av电影在线进入| 久久 成人 亚洲| 欧美日本中文国产一区发布| 国产成人av激情在线播放| 欧美老熟妇乱子伦牲交| 黄片播放在线免费| www.999成人在线观看| 亚洲精品国产色婷婷电影| 巨乳人妻的诱惑在线观看| 国产精品亚洲av一区麻豆| 精品人妻1区二区| 国产国语露脸激情在线看| 999久久久精品免费观看国产| 久久久久九九精品影院| 久久精品91蜜桃| 99riav亚洲国产免费| 侵犯人妻中文字幕一二三四区| 高清黄色对白视频在线免费看| 搞女人的毛片| 男女下面插进去视频免费观看| 亚洲自偷自拍图片 自拍| 久久人妻福利社区极品人妻图片| 熟妇人妻久久中文字幕3abv| 欧美成人性av电影在线观看| 9热在线视频观看99| cao死你这个sao货| 一级a爱视频在线免费观看| 一级片免费观看大全| 又黄又爽又免费观看的视频| 老鸭窝网址在线观看| 国产乱人伦免费视频| 亚洲aⅴ乱码一区二区在线播放 | 中文字幕精品免费在线观看视频| 日本a在线网址| 亚洲va日本ⅴa欧美va伊人久久| 黄频高清免费视频| 手机成人av网站| 男女下面插进去视频免费观看| 色婷婷久久久亚洲欧美| 久久久久久久精品吃奶| 侵犯人妻中文字幕一二三四区| 一级,二级,三级黄色视频| 国产精品99久久99久久久不卡| 亚洲黑人精品在线| 激情在线观看视频在线高清| 亚洲五月天丁香| 亚洲精品一卡2卡三卡4卡5卡| av网站免费在线观看视频| 曰老女人黄片| 欧美成人一区二区免费高清观看 | 欧美午夜高清在线| 成人18禁在线播放| 色综合亚洲欧美另类图片| 级片在线观看| 亚洲第一青青草原| 欧美一级a爱片免费观看看 | 99精品欧美一区二区三区四区| xxx96com| 国产三级黄色录像| 香蕉久久夜色| 久久精品国产亚洲av高清一级| 亚洲国产精品合色在线| 麻豆av在线久日| 精品少妇一区二区三区视频日本电影| 欧美乱妇无乱码| 国产一区二区三区视频了| 真人一进一出gif抽搐免费| 99久久久亚洲精品蜜臀av| 免费在线观看黄色视频的| 97碰自拍视频| 狂野欧美激情性xxxx| 一级a爱片免费观看的视频| 国产亚洲欧美精品永久| 色综合亚洲欧美另类图片| 国产日韩一区二区三区精品不卡| 国产一区二区三区综合在线观看| 亚洲av片天天在线观看| 18禁美女被吸乳视频| 欧美国产精品va在线观看不卡| 欧美在线一区亚洲| 亚洲中文字幕日韩| 99香蕉大伊视频| 国产精品久久久久久亚洲av鲁大| 精品久久蜜臀av无| av视频免费观看在线观看| 制服人妻中文乱码| 亚洲精品国产精品久久久不卡| 一区二区三区国产精品乱码| 欧美日韩乱码在线| 国产精品1区2区在线观看.| 一边摸一边抽搐一进一小说| 国产一级毛片七仙女欲春2 | 女人爽到高潮嗷嗷叫在线视频| 性色av乱码一区二区三区2| 亚洲午夜精品一区,二区,三区| 1024视频免费在线观看| 欧美亚洲日本最大视频资源| 久久久精品国产亚洲av高清涩受| 97人妻精品一区二区三区麻豆 | 日日夜夜操网爽| 国产精品爽爽va在线观看网站 | 国产亚洲av高清不卡| 久久午夜亚洲精品久久| 人人妻,人人澡人人爽秒播| 少妇的丰满在线观看| 国产高清激情床上av| 国产精品一区二区在线不卡| 女生性感内裤真人,穿戴方法视频| 亚洲午夜理论影院| 69av精品久久久久久| 黑人巨大精品欧美一区二区mp4| 欧美成狂野欧美在线观看| 18禁观看日本| 老司机午夜福利在线观看视频| 国产主播在线观看一区二区| 黄片播放在线免费| 嫩草影院精品99| 女性被躁到高潮视频| 在线观看午夜福利视频| 纯流量卡能插随身wifi吗| 国产区一区二久久| 日韩欧美一区视频在线观看| 啦啦啦观看免费观看视频高清 | 国产精品免费视频内射| tocl精华| 日本五十路高清| 欧美成人免费av一区二区三区| 国产高清激情床上av| 亚洲在线自拍视频| 黄色视频,在线免费观看| 美女国产高潮福利片在线看| 国产精品乱码一区二三区的特点 | 一进一出抽搐gif免费好疼| av片东京热男人的天堂| 日韩欧美一区视频在线观看| 脱女人内裤的视频| 午夜亚洲福利在线播放| 一边摸一边做爽爽视频免费| 99久久国产精品久久久| 久久狼人影院| 亚洲欧美日韩高清在线视频| 成在线人永久免费视频| 亚洲欧美精品综合一区二区三区| 久久久国产成人免费| 国产男靠女视频免费网站| 精品久久久久久成人av| 天堂√8在线中文| 亚洲一区二区三区不卡视频| 久久精品91无色码中文字幕| 日本免费一区二区三区高清不卡 | 日韩欧美三级三区| 亚洲av片天天在线观看| 电影成人av| 12—13女人毛片做爰片一| 免费不卡黄色视频| 色综合亚洲欧美另类图片| 国产激情欧美一区二区| 亚洲七黄色美女视频| 黄网站色视频无遮挡免费观看| 热99re8久久精品国产| 精品国产一区二区三区四区第35| 国产成人av激情在线播放| 大陆偷拍与自拍| 一边摸一边抽搐一进一出视频| 久久九九热精品免费| 日韩精品青青久久久久久| 亚洲精品国产一区二区精华液| 国产成年人精品一区二区| 99久久99久久久精品蜜桃| 久久久久久久久久久久大奶| 久久久久久人人人人人| www国产在线视频色| 欧美老熟妇乱子伦牲交| 亚洲国产日韩欧美精品在线观看 | 老司机深夜福利视频在线观看| 欧美午夜高清在线| 91大片在线观看| 精品国产美女av久久久久小说| 黑人欧美特级aaaaaa片| 欧美黄色片欧美黄色片| 成人18禁高潮啪啪吃奶动态图| 久热这里只有精品99| 久久草成人影院| 亚洲专区中文字幕在线| 男人的好看免费观看在线视频 | 50天的宝宝边吃奶边哭怎么回事| 国产欧美日韩精品亚洲av| 一级毛片精品| 一个人免费在线观看的高清视频| 中文字幕精品免费在线观看视频| av视频在线观看入口| 国产精品免费视频内射| 国产亚洲av嫩草精品影院| 久热爱精品视频在线9| 黄片大片在线免费观看| 久久婷婷成人综合色麻豆| 日韩有码中文字幕| 亚洲伊人色综图| 欧美成人免费av一区二区三区| 色综合亚洲欧美另类图片| 成人18禁高潮啪啪吃奶动态图| 一二三四社区在线视频社区8| 黄色毛片三级朝国网站| 国产亚洲欧美在线一区二区| 在线国产一区二区在线| 淫秽高清视频在线观看| 欧美乱码精品一区二区三区| 日韩三级视频一区二区三区| 脱女人内裤的视频| 国产一区二区三区综合在线观看| 欧美久久黑人一区二区| 国产91精品成人一区二区三区| 精品午夜福利视频在线观看一区| 亚洲成人精品中文字幕电影| 精品久久蜜臀av无| av欧美777| 欧美国产日韩亚洲一区| 黄色女人牲交| 黄色丝袜av网址大全| 久久中文看片网| 9色porny在线观看| 丰满的人妻完整版| 丰满人妻熟妇乱又伦精品不卡| 亚洲片人在线观看| 日本在线视频免费播放| 亚洲精品久久成人aⅴ小说| 亚洲av电影不卡..在线观看| 国产在线观看jvid| 亚洲中文字幕日韩| 好男人在线观看高清免费视频 | 999久久久国产精品视频| 久99久视频精品免费| 日本五十路高清| 91麻豆av在线| 日韩 欧美 亚洲 中文字幕| 亚洲欧洲精品一区二区精品久久久| 88av欧美| 亚洲午夜理论影院| 人成视频在线观看免费观看| 日韩中文字幕欧美一区二区| 一区二区三区精品91| 久久久久久免费高清国产稀缺| 一级,二级,三级黄色视频| 国产精品98久久久久久宅男小说| 亚洲aⅴ乱码一区二区在线播放 | 色精品久久人妻99蜜桃| 麻豆成人av在线观看| 日韩大码丰满熟妇| 又紧又爽又黄一区二区| 久久久久九九精品影院| 亚洲精品国产区一区二| 禁无遮挡网站| 亚洲国产欧美日韩在线播放| 久9热在线精品视频| 女人被躁到高潮嗷嗷叫费观| 热99re8久久精品国产| 精品人妻1区二区| 欧美丝袜亚洲另类 | 国产成+人综合+亚洲专区| 久久人妻熟女aⅴ| 超碰成人久久| 啪啪无遮挡十八禁网站| 在线观看66精品国产| 91国产中文字幕| 国产欧美日韩一区二区三区在线| bbb黄色大片| 97人妻精品一区二区三区麻豆 | www.www免费av| 日本黄色视频三级网站网址| 99在线视频只有这里精品首页| 纯流量卡能插随身wifi吗| 非洲黑人性xxxx精品又粗又长| 91字幕亚洲| 亚洲成a人片在线一区二区| 美女大奶头视频| 美女扒开内裤让男人捅视频| 亚洲精品久久成人aⅴ小说| 亚洲国产精品久久男人天堂| 精品国产美女av久久久久小说| 看免费av毛片| 黄频高清免费视频| 啦啦啦观看免费观看视频高清 | 日日爽夜夜爽网站| 国内毛片毛片毛片毛片毛片| 757午夜福利合集在线观看| 中国美女看黄片| av免费在线观看网站|