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

    Effects of Expansion Chamber on the Flow Pattern and Performance of Cyclones

    2017-08-07 21:32:09HanXiaoSunGuogangSunZhanpengHanXiaopeng
    中國煉油與石油化工 2017年2期

    Han Xiao; Sun Guogang,2; Sun Zhanpeng; Han Xiaopeng

    (1. State Key Laboratory of Heavy Oil Processing, China University of Petroleum, Beijing 102249; 2. Beijing Key Laboratory of Process Fluid Filtration and Separation, Beijing 102249)

    Effects of Expansion Chamber on the Flow Pattern and Performance of Cyclones

    Han Xiao1; Sun Guogang1,2; Sun Zhanpeng1; Han Xiaopeng1

    (1. State Key Laboratory of Heavy Oil Processing, China University of Petroleum, Beijing 102249; 2. Beijing Key Laboratory of Process Fluid Filtration and Separation, Beijing 102249)

    To improve the performance of the Stairmand cyclone, the effects of an expansion chamber on the fow feld, the pressure drop and the separation efficiency were investigated numerically and experimentally. The experimental results showed that compared with the Stairmand cyclone, the cyclone with an upper expansion chamber worked better at low inlet velocity (less than 14 m/s in this study), while the cyclone with a lower expansion chamber achieved higher effciency at a relatively high inlet velocity (14-20 m/s). The presence of an expansion chamber can generally result in a slight decrease in the cyclone pressure drop. The simulated results, which were used to further analyze the reason behind the experimental phenomena, suggested that the expansion chamber had insignifcant effects on the tangential velocity profles in the cylindrical part of cyclones. While in the cone part, the expansion chamber and the dipleg, the tangential velocity slightly decreased. Nevertheless, the expansion chamber decreased the possibility of the vortex end to sweep the wall and then reduce the particle re-entrainment. Therefore, installing the expansion chamber at a proper position could improve the separation performance of Stairmand cyclones. Both the experimental and simulated results represent a potential improvement of the Stairmand cyclone performance.

    expansion chamber; numerical simulation; fow feld; separation effciency

    1 Introduction

    Cyclones are typical dust collection devices that utilize the centrifugal force to capture particles. Because of their simplicity in construction, low cost in operating and maintenance, as well as adaptability to harsh circumstance, cyclones are widely used in petrochemical industry[1].

    The cyclone performance could be optimized by adjusting different geometrical parameters including the cyclone length, the inlet, and the vortex fnder geometries. Some researchers have paid their attention to effects of the lower part of cyclones. Obermair, et al.[2]tested five cyclones with different dust outlet geometries, and proved that changing the dust outlet geometry could influence the separation effciency signifcantly. Besides, by using a downcomer tube between the cyclone cone and the dustbin, a high separation efficiency could be obtained with an acceptable pressure drop. Also, Obermair, et al.[3]reaffirmed that the best separation efficiency could be achieved by a cyclone with a downcomer tube and bin, due to the “additional” separation process taking place within it. Elsayed and Lacor[4]examined the effect of the dustbin on the flow pattern and the performance of cyclones, and reported that the maximum tangential velocity was almost unchanged while the axial velocity could exhibit an inverted V or inverted W profile. Gao Cuizhi[5]made the comparison between four cases, viz.: the cyclone without dustbin, the cyclone with half a dustbin, with one dustbin, and with double dustbins. The results showed that at high velocity the cyclone with half a dustbin could obtain the best separation performance.

    The particle re-entrainment from the dustbin is an important factor infuencing the separation effciency[6-7]. Mothes[8]found that higher particle concentration was near the cyclone bottom rather than at the inlet, which indicated a requirement for control of particle reentrainment from the dustbin. Researchers have optimizedcyclones for installing the apex cone and down-comer tubes for the purpose of reducing particle re-entrainment. Yoshida[7,9]investigated the effect of apex cone shape and height on the separation performance. Sakura, et al.[10]investigated the performance of cyclones modified with down-comer tubes and proved the effectiveness of downcomer tubes in reducing the particle re-entrainment. Qian Fuping, et al.[6]examined the vertical tube length and found that the tangential velocity, axial velocity and the turbulent kinetic energy in the dustbin reduced greatly when the vertical tube was connected to the dust outlet, which could avoid the re-entrainment effectively.

    In order to decrease the particle re-entrainment, the expansion chamber has been introduced in this paper to study the effect of the expansion chamber on the flow pattern, velocity profles, and separation performance and to compare different axial positions of expansion chamber at different inlet velocities. This approach serves two purposes: (1) It gives information for designers about the effect of the expansion chamber on the fow pattern and separation performance of the Stairmand high-effciency cyclones. (2) It determines the proper working range of the Stairmand high-effciency cyclones with and without an expansion chamber.

    2 Materials and Methods

    2.1 Experimental setup

    Three tested cyclones are shown in Figure 1. The S cyclone is a standard Stairmand cyclone with a dipleg underneath, while both the SA cyclone and the SB cyclone have an expansion chamber respectively located above and below the connecting line (CL) of the Stairmand cyclone between the cone and the dipleg. In the SA cyclone, the expansion chamber is located above CL while the SB cyclone is located below CL. The Stairmand cyclone is used as the benchmark. Table 1 gives the dimensions of cyclones. The test facility is illustrated in Figure 2. The experiments were carried out under negative pressure by a suction fan. Standard pitot tubes were used to measure the inlet velocity and the pressure drop between the inlet and the outlet was tested by a U-tube. The separation effciency could be achieved by weighing the powder before and after the tests. The solid particles used were talcum powder, with a physical density of 2 700 kg/m3and a volume median diameter of 13.59 μm. The feed material concentration was fxed at 20 g/m3.

    Figure 1 Schematic diagram of three tested cyclones

    Table 1 Geometrical dimensions of the three cyclones

    Figure 2 Schematic diagram of experiment setup

    2.2 CFD simulation

    Computational Fluid Dynamics (CFD) simulations were performed to investigate the effects of the expansion chamber on fow feld using the commercial CFD code Fluent 14.0. The Gambit 2.4.6 was used as the preprocessing tool of Fluent for establishing the solid model and mesh generation. Pure-phase simulation has been proved to be an effective way to analyze the flow field, from the results of which the velocity profiles and the static pressure distribution could be observed and this simulation method has been employed in many studies to support the experimental results[6,11-13]. Therefore, the pure-phase simulation was employed in this paper.

    For the three-dimensional turbulent fow in cyclones, the key to the success of CFD lies in the accurate description of the turbulent behavior of the flow. The Reynolds stress turbulence model (RSM) has been proved to be an accurate prediction on the tangential velocity, the axial velocity and the pressure drop in the cyclone[14]. It is capable of predicting the combined vortex in accordance with the experimental data[15-16]. In this study, the Reynolds stress turbulence model (RSM) was adopted to simulate the fow feld in cyclones.

    The boundary condition at inlet section was the velocityinlet with an inlet velocity of 18 m/s. An outflow boundary condition was used at the outlet. The noslip boundary condition was used at other boundaries. The PRESTO scheme has been used for the pressure interpolation, and the SIMPLE algorithm has been used for pressure velocity coupling, while the QUICK scheme has been used to deal with the momentum equations, the turbulent kinetic energy, the turbulent dissipation rate, and the Reynolds stresses.

    Figure 3 shows the surface meshes of the three cyclones. In our preliminary computations, we have run the same simulation on three progressively finer grids with 129 840, 182 730 and 269 565 cells used for SA cyclone, respectively. The computational results on the three grids are shown in Table 2. Since the maximum difference between the results is less than 5%, so it could be argued that the simulations are grid independent.

    Furthermore, to verify the precision of simulation computation, comparison with the tangential velocity distribution measured by Hu Liyuan[17]with LDV is made at the sectionZ=340 mm. As Figure 4 indicates a good agreement between the experimental and the numerical results, it is concluded that the numerical simulation could be used with acceptable accuracy.

    Figure 3 Surface meshes for three cyclones

    Table 2 Results of grid independent study for SA cyclone

    Figure 4 Comparison of the tangential velocity distribution between experimental and simulated results

    3 Results and Discussion

    3.1 Experimental results

    As shown in Figure 5 (a), the expansion chamber improved the separation efficiency by a maximum of 1.37%, and decreased the pressure drop by a maximum of 16%. For ease of clear expression, we separate the efficiency curve into three regions. Separation effciencies generally increase to a maximum frst and then decrease with the increase of inlet velocity. It is of interest to note that the value of optimal inlet velocity for different cyclones differs greatly. For Region 1, SA cyclone has a higher effciency. With further increasing the inlet velocity, the effciency of the SA cyclone begins to decrease, while the effciency of other two types of cyclones continues to rise. When the inlet velocity is between 15- 20 m/s (Region 2), the SB cyclone behaves better in terms of the particle separation. But once the inlet velocity exceeds 20 m/s (Region 3), the efficiency of SB cyclone decreases dramatically and the S cyclone obtains a highest effciency. The above results indicate that we should move upward the expansion chamber in the region of low inlet velocity (<15 m/s). On the other hand, when the inlet velocity varies from 14 m/s to 20 m/s, it is suitable to move downward the expansion chamber. But if the cyclone works at a high inlet velocity (>20 m/s), the expansion chamber causes negative effect on efficiency and we would better not install the expansion chamber under such condition.

    Figure 5 Effect of an expansion chamber on separation efficiency (a) and pressure drop (b)

    Figure 5 (b) shows the effect of an expansion chamber on the pressure drops. When the inlet velocity is less than 18 m/s, the variation of pressure drop caused by expansion chamber is almost negligible. However, when the inlet velocity exceeds 20 m/s, moving up the expansion chamber can decrease the pressure drop; otherwise, the pressure drop would increase.

    3.2 CFD results

    3.2.1 The effect of expansion chamber on the velocity profiles

    The tangential and axial velocities make a big difference on the separation efficiency. Especially, the tangential velocity is the dominant component of the gas flow in cyclones, which results in the centrifugal force for particle separation[18]. To make it convenient to signal different sections, Figure 6 is the sketch on different sections of the cyclone. Figure 7 and Figure 8 show the tangential and axial velocity profles in different sections at an inlet velocity of 18 m/s, and it is found that there is no big difference among tangential velocities in the cylinder and the conic space, and the highest tangential velocity value is almost the same. This conclusion is consistent with the literature report of Elsayed and Lacor[19].

    Figure 6 Different sections of cyclone for capturing the velocity field

    Figure 7 Radial profiles for tangential velocity in different sections operating at an inlet velocity of 18 m/s

    Figure 8 Radial profiles for axial velocities in different sections operating at an inlet velocity of 18 m/s

    The tangential velocities of SA cyclone are inferior to those of the other two cyclones entering the cone, illustrating that the expansion chamber can decrease the tangential velocity subtly and lower the tangential velocity, resulting in smaller centrifugal force making fewer particles captured by the wall, which in turn would lead to lower separation efficiency. This explains well why SA cyclone has inferior separation effciency at an inlet velocity of 18 m/s. For the case of SB cyclone, owing to the low velocity, the tangential velocity decreases dramatically after entering the expansion chamber. Meanwhile, it is obvious to find that SB cyclone features a better symmetry for tangential velocity than the other two cyclones atZ=-1 400 mm. This situation can to some extent help reduce the dust-entrainment and enhance the separation effciency.

    As for the axial velocity, it is found that the axial velocity varies between each other. The axial velocity of S cyclone appears to be the shape of an inverted V profile, and the SA and SB cyclones have the shape of inverted W. The expansion chamber has provided SA cyclone with a higher downward velocity at outer vortex, thus making downward flow as a main factor causing severe dust re-entrainment. So the SA cyclone shows an inferior performance at high inlet velocity. In a similar way, the main rotational flow extends downwards with a growing inlet velocity, while the same axial profles turns up for the SB cyclone, which can also explain why SB cyclone is inferior to S cyclone at higher inlet velocity.

    As the section extends downwards, it is clear that SB cyclone has the lowest axial velocity along the central axis, which reveals that the central upward fow has the longest residence time for the secondary separation. Besides, for S cyclone and SA cyclone, the axial velocities turn up to change from double-peak to the single-peak and then weaken in symmetry, leading to instability of the fow feld. As for the SB cyclone, the axial velocity profiles maintain the appearance of double-peak, which is quite good in symmetry until entering the dipleg. Therefore, the SB cyclone is superior to the other two types of cyclones at an inlet velocity of 18 m/s.

    3.2.2 The effect of expansion chamber on the static pressure

    Figure 9 Contours of static pressure atVin=18 m/s andVin=25 m/s

    Figure 9 shows the contour of static pressure at an inlet velocity(Vin) of 18 m/s and 25 m/s, respectively. It can be seen from Figure 9 that the vortex end of S cyclone falls between the other two cyclones at an inlet velocity of 18 m/s. In comparison with the S cyclone, the vortex achieves its longest length in the SB cyclone. On onehand, the SB cyclone makes its vortex end located in the dipleg, which reduces the dust re-entrainment. On the other hand, the SB cyclone obtains longer residence time for particles, resulting in a better separation performance. This phenomenon matches the experimental results well. In a similar way, when the inlet velocity increases to a relatively high value, the vortex end goes through the expansion chamber in SB cyclone and the expansion chamber weakens the swirl of vortex thus shortening the natural vortex length, resulting in shorter natural vortex length and low separation effciency. This fact has been proved by simulated results depicted in Figure 9. When the inlet velocity reaches 25 m/s, the vortex achieves its longest length in S cyclone and a shortest length in SB cyclone. Therefore, the S cyclone achieves a highest separation efficiency, which matches the experimental results well.

    Figure 10 shows the variations of static pressure at different inlet velocity ranging from 10 m/s to 25 m/s. It is suggested that with an increasing inlet velocity, the vortex end is approaching the bottom of dipleg, resulting in the extension of separation space. However, there appears to be swing of vortex end under such condition, resulting in instability of the fow feld and decrease of the separation effciency.

    Figure 10 Contours of static pressure for SA cyclone operating at different inlet velocities

    With the increase in inlet velocity, the vortex end of SA cyclone extends out of the expansion chamber, which would result in more possibilities of contact with dipleg wall and dust re-entrainment, which could decrease the separation efficiency. Figure 10 illustrates that a fierce swing of the vortex end happens through the whole dipleg at an inlet velocity of 25 m/s, which would bring about much dust-entrainments. This is why the SA cyclone fails to enhance efficiency in spite of a comparatively larger vortex length as compared with S cyclone at an inlet velocity of 25 m/s. Therefore, for the case of SA cyclone, an increased inlet velocity is not ideal for the separation effciency.

    Besides, it is hard for SA cyclone to form turbulent fow based on a weaker cone contraction effect. Also, an upper expansion chamber would effectively extend the vortex at low inlet velocities. Therefore, the SA cyclone has a superior performance at low inlet velocities. An excessively large volume of downward gas fow at high inlet velocities would bring on severe dust re-entrainment and decreased effciency.

    3.2.3 The effect of expansion chamber on the axial velocity vectors

    Figure 11 Distribution of axial velocity vector for three cyclones atVin=10 m/s

    Besides the tangential velocity, the dust-entrainment is also a main factor affecting the separation efficiency. It can be seen from Figure 11 that the dust-entrainment in different degree appears in three cyclones at the inlet velocity of 10 m/s. Despite the longer vortex length, the separation efficiencies of S cyclone and SB cyclone arereduced because of severe dust-entrainment appearing as eddies in the axial velocity vector diagram (as seen from the red circles), which would bring the collected dust off the wall to enter the exit tube. As a consequence, SA cyclone achieves the best separation performance with less dust-entrainment at an inlet velocity of 10 m/s, which agrees well with the experimental results.

    3.2.4 The effect of expansion chamber on the central static pressure

    Figure 12 Position versus static pressure at different inlet velocities

    To further analyze the reason why the SB cyclone has an advantage over SA cyclone at the inlet velocity of 18 m/s but fails to assume a leading position at an inlet velocity exceeding 18 m/s, it is important to state that the main separation domain covers from the vortex fnder to the end of vortex, so the vortex length is critical for the separation effciency. At the end of vortex, there is a sharp change of static pressure. As shown in Figure 12, the vortex end of SA cyclone is far away from the other two cyclones at the inlet velocity of 18 m/s, thus obtaining the lowest effciency. While, it is hard to distinguish the break points of static pressure among S cyclone and SB cyclone at the inlet velocity of 18 m/s, but the SB cyclone demonstrates a relatively stable fow feld after the break point with a minimum static pressure variation. Upon taking the vortex length and dust re-entrainment into consideration, the SB cyclone obtains the best separation performance at the inlet velocity of 18 m/s, which matches the experimental results well. The tangential velocity, the vortex length and the dust re-entrainment are main factors for weighing the efficiency. Therefore, the tangential velocity, the vortex length and the degree of dust-entrainment are vital factors adopted to estimate the separation effciency. When the inlet velocity reaches 25 m/s, it is clear that the S cyclone obtains the longest vortex length, and besides this situation there appears to be a second jump in static pressure after the break point for SA cyclone. It brings the collected particles off the dipleg to re-enter the inner fow, thus decreasing the separation effciency.

    4 Conclusions

    Three different types of cyclones have been investigated using CFD and experimental method. The following conclusions have been drawn up to give more information about the effects of the expansion chamber:

    (1) The tangential velocities obtained in different cyclones are very similar in the cylindrical part, but some differences appear in the conical part, the expansion chamber, and the dipleg. The expansion chamber could decrease the tangential velocity and also reduce the possibilities of the vortex end to sweep the internal wall, leading to enhancement of the separation efficiency. Therefore, the tangential velocity is not the main factor deciding the separation efficiency, which is assessed based on the gas re-entrainment.

    (2) The axial velocity patterns appear to be similar in the cylindrical part for the three cyclones. For areas below the separation zone, the cyclone with the lower expansion chamber has the lowest axial velocity along the central axis, which has revealed that the central upward fow has the longest residence time for the secondary separation.

    (3) Compared with the Stairmand cyclone, the cyclone with an upper expansion chamber works better at low inlet velocity (less than 14 m/s in this study), while the cyclone with a lower expansion chamber obtains higher effciency when the inlet velocity varies from 14 m/s to 20 m/s. Thus, with regard to the design of a Stairmandcyclone, we believe that it would be possible to design a cyclone with improved separation performance by considering the expansion chamber position with a proper inlet gas velocity.

    Acknowledgement: The authors gratefully acknowledge the financial assistance from the National Natural Science Foundation (Grant No.21276274) and the National Key Project of Basic Research of the Ministry for Science and Technology of China (Grant No. 2014CB744304).

    [1] Ji Z L, Wu X L, Shi M X. Experimental research on the natural turning length in the cyclone[J]. Acta Petrolei Sinica (Petroleum Processing Edition), 1993(4): 86-91 (in Chinese)

    [2] Obermair S, Woisetschl?ger J, Staudinger G. Investigation of the flow pattern in different dust outlet geometries of a gas cyclone by laser Doppler anemometry[J]. Powder Technology, 2003, 138(2): 239-251

    [3] Obermair S, Staudinger G. The dust outlet of a gas cyclone and its effects on separation efficiency[J]. Chemical Engineering & Technology, 2001, 24(12): 1259-1263

    [4] Elsayed K, Lacor C. The effect of the dust outlet geometry on the performance and hydrodynamics of gas cyclones[J]. Computers & Fluids, 2012, 68: 134–147

    [5] Gao C Z, Flow Characteristics of Inner Vortex and Its Effect on Separation Performance of Cyclones[D]. China University of Petroleum (Beijing), 2011

    [6] Qian F P, Zhang J G, Zhang M Y. Effects of the prolonged vertical tube on the separation performance of a cyclone.[J]. Journal of Hazardous Materials, 2006, 136(3): 822-829

    [7] Yoshida H, Yang K S, Fukui K, et al. Effect of apex cone height on particle classifcation performance of a cyclone separator [J]. Advanced Powder Technology, 2003, 14(3): 263-278

    [8] Helmut M, Friedrich L. Bewegung und Abscheidung der Partikeln im Zyklon[J]. Chemie Ingenieur Technik, 1984, 56(9): 714-715

    [9] Yoshida H, Nishimura Y, Fukui K, et al. Effect of apex cone shape on fne particle classifcation of gas-cyclone[J]. Powder Technology, 2010, 204(1): 54–62

    [10] Ganegama B S, Leung A Y. Improvements of the cyclone separator performance by down-comer tubes [J]. Journal of Hazardous Materials, 2016, 311: 100-114

    [11] Fu X, Sun G G, Liu S X et.al. Numerical simulation of flow field in annular space of single-and double-inlet cyclones[J]. Petroleum Refnery Engineering, 2010, 40(8): 26-30 (in Chinese)

    [12] Kaya F, Karagoz I. Numerical investigation of performance characteristics of a cyclone prolonged with a dipleg[J]. Chemical Engineering Journal, 2009, 151(1–3): 39-45

    [13] Fu X, Sun G G, Liu J, et al. Effect of induced velocity on separation effciency and pressure drop of cyclones caused by vortex in vortex-tube inlet pipe[J]. Journal of Chemical Industry and Engineering, 2011, 62(7): 1927-1932

    [14] Hoffmann A C, Stein L E. Gas Cyclones and Swirl Tubes: Principles, Design and Operation[M]. Springer, 2004: 1-34

    [15] Hoekstra A J. Gas Flow Field and Collection Effciency of Cyclone Separators [D]. 2000

    [16] Kaya F, Karagoz I, Kaya F, et al. Performance analysis of numerical schemes in highly swirling turbulent flows in cyclones[J]. Current Science, 2008, 94(10): 1273-1278

    [17] Hu L Y, Shi M X. Study on flow structure of strong turbulent flow field in cyclone separator. Journal of Chemical Industry and Engineering, 2003(54): 549-556

    [18] Xiang R B, Lee K W. Numerical study of flow field in cyclones of different height[J]. Chemical Engineering & Processing Process Intensifcation, 2005, 44(8): 877-883

    [19] Elsayed K, Lacor C. The Effect of Cyclone Dustbin on the Flow Pattern and Performance[C]// Tenth International Congress of Fluid Dynamics, 2010.

    date: 2017-01-18; Accepted date: 2017-03-14.

    Prof. Sun Guogang, Telephone: +86-10- 89734820; E-mail: ggsunbj@163.com.

    两性午夜刺激爽爽歪歪视频在线观看| 韩国av一区二区三区四区| 国内精品久久久久精免费| 一个人看视频在线观看www免费 | 欧美黑人巨大hd| 亚洲人与动物交配视频| 不卡一级毛片| 国产一区二区在线观看日韩 | 两个人看的免费小视频| 激情在线观看视频在线高清| 一区福利在线观看| 国产精品久久视频播放| 91麻豆av在线| 99久久99久久久精品蜜桃| 亚洲成人久久性| 香蕉av资源在线| 欧美zozozo另类| 国产69精品久久久久777片| 亚洲av中文字字幕乱码综合| 神马国产精品三级电影在线观看| 久久久国产精品麻豆| 欧美成狂野欧美在线观看| 特大巨黑吊av在线直播| 激情在线观看视频在线高清| 精品国产美女av久久久久小说| 超碰av人人做人人爽久久 | 成人特级黄色片久久久久久久| 搡女人真爽免费视频火全软件 | 天堂网av新在线| 又紧又爽又黄一区二区| 岛国视频午夜一区免费看| 99国产精品一区二区蜜桃av| 性欧美人与动物交配| 9191精品国产免费久久| 精品久久久久久,| 国产一区二区三区在线臀色熟女| 日本撒尿小便嘘嘘汇集6| 成人一区二区视频在线观看| 亚洲专区中文字幕在线| 99国产综合亚洲精品| 国产美女午夜福利| 好看av亚洲va欧美ⅴa在| 波多野结衣高清无吗| 国产一区二区亚洲精品在线观看| 国产麻豆成人av免费视频| 怎么达到女性高潮| 日韩欧美 国产精品| 天堂动漫精品| tocl精华| 人人妻人人看人人澡| 日韩免费av在线播放| 波多野结衣高清作品| 桃红色精品国产亚洲av| 免费人成在线观看视频色| 国产免费一级a男人的天堂| 岛国视频午夜一区免费看| 日韩高清综合在线| 脱女人内裤的视频| 午夜福利在线在线| 中文字幕熟女人妻在线| 老司机午夜福利在线观看视频| 成人永久免费在线观看视频| 午夜a级毛片| 色播亚洲综合网| 内地一区二区视频在线| 一个人看的www免费观看视频| 精品国产三级普通话版| 亚洲精华国产精华精| 亚洲国产精品久久男人天堂| 九九在线视频观看精品| 99在线人妻在线中文字幕| 久久午夜亚洲精品久久| 一区二区三区国产精品乱码| 97人妻精品一区二区三区麻豆| a级毛片a级免费在线| 日韩高清综合在线| 成年人黄色毛片网站| 美女高潮喷水抽搐中文字幕| 十八禁网站免费在线| 一个人免费在线观看的高清视频| 最近在线观看免费完整版| 亚洲国产中文字幕在线视频| 国产精品乱码一区二三区的特点| 男人舔奶头视频| 久久久久久久久大av| 偷拍熟女少妇极品色| 少妇人妻精品综合一区二区 | 一进一出好大好爽视频| 听说在线观看完整版免费高清| 国产真实乱freesex| 国产精品久久久久久久久免 | 日韩av在线大香蕉| 欧美一区二区国产精品久久精品| av天堂中文字幕网| 国产精品一区二区三区四区免费观看 | 精品无人区乱码1区二区| 一本精品99久久精品77| 欧美日韩瑟瑟在线播放| 日本五十路高清| 欧美一区二区国产精品久久精品| 麻豆国产av国片精品| 国产aⅴ精品一区二区三区波| 亚洲av二区三区四区| 黄色日韩在线| 国产一区在线观看成人免费| 欧美性猛交╳xxx乱大交人| 熟女人妻精品中文字幕| xxx96com| 亚洲av五月六月丁香网| 久久国产乱子伦精品免费另类| 久久中文看片网| 精品乱码久久久久久99久播| 亚洲av一区综合| 嫁个100分男人电影在线观看| 内地一区二区视频在线| 国产午夜精品论理片| 三级国产精品欧美在线观看| 亚洲美女黄片视频| 在线观看美女被高潮喷水网站 | 国产又黄又爽又无遮挡在线| www.熟女人妻精品国产| 激情在线观看视频在线高清| 男人舔奶头视频| 在线观看一区二区三区| 麻豆国产97在线/欧美| 一个人观看的视频www高清免费观看| 国产欧美日韩一区二区精品| 亚洲精品乱码久久久v下载方式 | 午夜福利免费观看在线| 一区二区三区免费毛片| 亚洲人与动物交配视频| 国产精品电影一区二区三区| 女同久久另类99精品国产91| 日韩欧美精品v在线| 国产av不卡久久| 一本久久中文字幕| 日本在线视频免费播放| 最近最新中文字幕大全电影3| 色av中文字幕| 丁香六月欧美| 欧美乱色亚洲激情| 熟女电影av网| 日韩 欧美 亚洲 中文字幕| 男人舔奶头视频| 久久精品国产亚洲av香蕉五月| 99国产综合亚洲精品| 久久性视频一级片| 亚洲男人的天堂狠狠| 波野结衣二区三区在线 | 麻豆国产av国片精品| 天堂网av新在线| 中文字幕av成人在线电影| 亚洲性夜色夜夜综合| 美女 人体艺术 gogo| 国产精品99久久久久久久久| 国产一区二区在线观看日韩 | a级一级毛片免费在线观看| 嫩草影视91久久| 乱人视频在线观看| 欧美日韩亚洲国产一区二区在线观看| 国产精品1区2区在线观看.| 国产精品一及| 欧美乱妇无乱码| av专区在线播放| 免费看a级黄色片| av在线天堂中文字幕| 国产97色在线日韩免费| 国产亚洲精品综合一区在线观看| 又粗又爽又猛毛片免费看| 99在线人妻在线中文字幕| 久9热在线精品视频| 啦啦啦韩国在线观看视频| 每晚都被弄得嗷嗷叫到高潮| 亚洲,欧美精品.| 亚洲精华国产精华精| 国产高清videossex| 亚洲 欧美 日韩 在线 免费| 热99在线观看视频| 老鸭窝网址在线观看| 欧美另类亚洲清纯唯美| 天堂av国产一区二区熟女人妻| 亚洲成人免费电影在线观看| 久久精品夜夜夜夜夜久久蜜豆| 亚洲熟妇中文字幕五十中出| 国产精品久久视频播放| 性欧美人与动物交配| 午夜激情福利司机影院| 成年女人毛片免费观看观看9| 人妻丰满熟妇av一区二区三区| 女同久久另类99精品国产91| 美女大奶头视频| 九九久久精品国产亚洲av麻豆| 免费观看人在逋| 精品欧美国产一区二区三| 人妻久久中文字幕网| 国产激情欧美一区二区| av天堂中文字幕网| 亚洲国产精品sss在线观看| 亚洲无线在线观看| 一本久久中文字幕| 99久国产av精品| 欧美乱码精品一区二区三区| 女同久久另类99精品国产91| 久久精品夜夜夜夜夜久久蜜豆| 18禁在线播放成人免费| 亚洲第一欧美日韩一区二区三区| 小说图片视频综合网站| 一二三四社区在线视频社区8| 男人和女人高潮做爰伦理| 中文字幕高清在线视频| 久久欧美精品欧美久久欧美| 日韩欧美精品v在线| 亚洲成av人片在线播放无| 老汉色∧v一级毛片| 亚洲欧美日韩无卡精品| 国产精品av视频在线免费观看| 免费无遮挡裸体视频| 激情在线观看视频在线高清| 国产av在哪里看| 男女床上黄色一级片免费看| 日韩高清综合在线| 最新美女视频免费是黄的| 精品国产亚洲在线| 欧美一区二区国产精品久久精品| 欧美中文日本在线观看视频| 亚洲欧美日韩东京热| 日韩 欧美 亚洲 中文字幕| 免费观看人在逋| 蜜桃亚洲精品一区二区三区| 亚洲天堂国产精品一区在线| 成年女人看的毛片在线观看| АⅤ资源中文在线天堂| 日本黄色片子视频| 婷婷亚洲欧美| 男女床上黄色一级片免费看| 日本撒尿小便嘘嘘汇集6| 午夜日韩欧美国产| 欧美+亚洲+日韩+国产| 亚洲国产欧美人成| 国产午夜精品久久久久久一区二区三区 | www.熟女人妻精品国产| 亚洲人成网站在线播放欧美日韩| 美女黄网站色视频| 可以在线观看的亚洲视频| 国产伦一二天堂av在线观看| 丰满的人妻完整版| 久久精品91无色码中文字幕| 高清在线国产一区| 香蕉久久夜色| 啦啦啦免费观看视频1| 久久精品91无色码中文字幕| 天堂网av新在线| 午夜福利在线在线| 岛国在线免费视频观看| 国内久久婷婷六月综合欲色啪| 老熟妇乱子伦视频在线观看| www日本黄色视频网| 18禁美女被吸乳视频| 久久精品国产99精品国产亚洲性色| 小说图片视频综合网站| 深夜精品福利| 欧美日韩福利视频一区二区| 色综合婷婷激情| 国产亚洲精品久久久久久毛片| 五月伊人婷婷丁香| 制服丝袜大香蕉在线| 午夜精品久久久久久毛片777| 最近最新中文字幕大全电影3| 国语自产精品视频在线第100页| 女人高潮潮喷娇喘18禁视频| 久久久久久国产a免费观看| 国产中年淑女户外野战色| 午夜视频国产福利| 最好的美女福利视频网| 欧美一区二区国产精品久久精品| 最近在线观看免费完整版| 日日夜夜操网爽| 亚洲真实伦在线观看| 久久精品综合一区二区三区| 操出白浆在线播放| 99久久精品一区二区三区| 亚洲专区中文字幕在线| 俄罗斯特黄特色一大片| 99久久综合精品五月天人人| 99视频精品全部免费 在线| 国产成人福利小说| 亚洲欧美日韩高清在线视频| 亚洲精品色激情综合| 国产精品久久电影中文字幕| 亚洲自拍偷在线| 亚洲成人久久爱视频| 亚洲国产高清在线一区二区三| 日韩av在线大香蕉| 少妇高潮的动态图| 亚洲七黄色美女视频| 看黄色毛片网站| 亚洲熟妇中文字幕五十中出| 天堂网av新在线| 一进一出好大好爽视频| 亚洲第一电影网av| 搡女人真爽免费视频火全软件 | 国产午夜福利久久久久久| 亚洲黑人精品在线| 国产探花极品一区二区| 国产精品女同一区二区软件 | 久久久精品大字幕| 国产精品,欧美在线| 欧美极品一区二区三区四区| 久久精品国产自在天天线| 首页视频小说图片口味搜索| av中文乱码字幕在线| 在线观看日韩欧美| 日韩免费av在线播放| 真人一进一出gif抽搐免费| 在线观看免费视频日本深夜| 性色avwww在线观看| 欧美精品啪啪一区二区三区| 欧美成人性av电影在线观看| 在线十欧美十亚洲十日本专区| 国产精品女同一区二区软件 | 脱女人内裤的视频| 九色成人免费人妻av| 国产精品一区二区免费欧美| 99精品欧美一区二区三区四区| 精品日产1卡2卡| 日本与韩国留学比较| 午夜激情欧美在线| 午夜免费成人在线视频| 国产视频内射| 午夜影院日韩av| 免费在线观看成人毛片| 久久欧美精品欧美久久欧美| 在线看三级毛片| 伊人久久大香线蕉亚洲五| 麻豆成人午夜福利视频| 日韩免费av在线播放| av女优亚洲男人天堂| 亚洲一区二区三区色噜噜| 日韩高清综合在线| 日本黄色视频三级网站网址| 午夜福利视频1000在线观看| 久久精品91蜜桃| 久久精品91无色码中文字幕| 99视频精品全部免费 在线| 久久精品91无色码中文字幕| 国产精品亚洲一级av第二区| 一本久久中文字幕| 国产一区二区三区视频了| 免费观看的影片在线观看| 免费av不卡在线播放| 欧美zozozo另类| 色老头精品视频在线观看| 69人妻影院| 久久精品亚洲精品国产色婷小说| av在线天堂中文字幕| 亚洲午夜理论影院| 少妇裸体淫交视频免费看高清| 一本综合久久免费| 午夜亚洲福利在线播放| 国产精品一区二区三区四区免费观看 | xxx96com| АⅤ资源中文在线天堂| 午夜视频国产福利| 欧美性感艳星| 亚洲最大成人中文| 日本 欧美在线| 欧美最新免费一区二区三区 | 午夜激情福利司机影院| 日本成人三级电影网站| 九九热线精品视视频播放| 日韩人妻高清精品专区| 麻豆国产97在线/欧美| 国产麻豆成人av免费视频| 日韩国内少妇激情av| 国产精品三级大全| 亚洲精品乱码久久久v下载方式 | 美女黄网站色视频| 少妇的丰满在线观看| 成年女人看的毛片在线观看| 亚洲欧美一区二区三区黑人| 啦啦啦观看免费观看视频高清| 午夜福利18| 夜夜爽天天搞| 欧美日韩瑟瑟在线播放| 国产一区二区三区在线臀色熟女| 国产精品久久电影中文字幕| 12—13女人毛片做爰片一| a在线观看视频网站| 国产成+人综合+亚洲专区| 久久久久久久久中文| 亚洲成人精品中文字幕电影| 麻豆国产av国片精品| 麻豆成人av在线观看| 露出奶头的视频| 国产成人啪精品午夜网站| 一进一出抽搐gif免费好疼| 亚洲av日韩精品久久久久久密| 窝窝影院91人妻| 欧美又色又爽又黄视频| 久久精品91无色码中文字幕| 看黄色毛片网站| 中文字幕人成人乱码亚洲影| 亚洲一区高清亚洲精品| 国内精品美女久久久久久| 性色av乱码一区二区三区2| 黑人欧美特级aaaaaa片| 国产三级中文精品| 欧美日韩一级在线毛片| 成人国产综合亚洲| 18禁在线播放成人免费| 欧美最黄视频在线播放免费| 精品国产美女av久久久久小说| 国产成人欧美在线观看| 动漫黄色视频在线观看| 日本 欧美在线| 久久久久国产精品人妻aⅴ院| 国产黄a三级三级三级人| 亚洲天堂国产精品一区在线| avwww免费| 非洲黑人性xxxx精品又粗又长| 亚洲 欧美 日韩 在线 免费| 国产一区二区亚洲精品在线观看| 啦啦啦韩国在线观看视频| 日韩中文字幕欧美一区二区| 免费人成视频x8x8入口观看| 少妇的逼水好多| 国产av不卡久久| 高潮久久久久久久久久久不卡| 亚洲欧美日韩无卡精品| 亚洲精品456在线播放app | 可以在线观看的亚洲视频| 99视频精品全部免费 在线| 免费高清视频大片| 国产成人av激情在线播放| 又粗又爽又猛毛片免费看| 熟女人妻精品中文字幕| 亚洲成人久久性| 精品人妻一区二区三区麻豆 | 亚洲五月婷婷丁香| 三级毛片av免费| 久久久久国产精品人妻aⅴ院| 99精品欧美一区二区三区四区| 法律面前人人平等表现在哪些方面| av天堂在线播放| 好男人在线观看高清免费视频| 亚洲欧美日韩高清在线视频| 亚洲真实伦在线观看| 最近在线观看免费完整版| 亚洲精品一卡2卡三卡4卡5卡| or卡值多少钱| 国内精品久久久久精免费| 2021天堂中文幕一二区在线观| 中国美女看黄片| 日韩国内少妇激情av| 亚洲内射少妇av| 亚洲人成网站在线播| 国产高潮美女av| 成人无遮挡网站| 伊人久久精品亚洲午夜| 国产一区二区在线av高清观看| www.熟女人妻精品国产| 五月伊人婷婷丁香| 久久久久性生活片| 国产一级毛片七仙女欲春2| 久99久视频精品免费| 精品人妻1区二区| 人妻久久中文字幕网| 美女黄网站色视频| 欧美高清成人免费视频www| 少妇的逼水好多| АⅤ资源中文在线天堂| 国产精品亚洲美女久久久| 最近最新免费中文字幕在线| 每晚都被弄得嗷嗷叫到高潮| 免费看日本二区| 一进一出好大好爽视频| 久久久久久久久中文| 女人十人毛片免费观看3o分钟| 亚洲成人免费电影在线观看| 老司机午夜十八禁免费视频| 九九久久精品国产亚洲av麻豆| 脱女人内裤的视频| 性色av乱码一区二区三区2| 国产视频一区二区在线看| bbb黄色大片| 欧美bdsm另类| 真实男女啪啪啪动态图| 黄色视频,在线免费观看| 婷婷精品国产亚洲av| 久9热在线精品视频| 亚洲欧美日韩无卡精品| 两个人看的免费小视频| 18禁在线播放成人免费| 女同久久另类99精品国产91| 一级黄片播放器| 搡女人真爽免费视频火全软件 | 男女做爰动态图高潮gif福利片| 一级毛片女人18水好多| 日本a在线网址| 午夜两性在线视频| 国内精品久久久久久久电影| av欧美777| 国产激情欧美一区二区| 丝袜美腿在线中文| 人妻夜夜爽99麻豆av| 又爽又黄无遮挡网站| 免费电影在线观看免费观看| 日韩欧美在线乱码| 亚洲精品久久国产高清桃花| 99国产极品粉嫩在线观看| 欧美国产日韩亚洲一区| 国产三级黄色录像| 成人国产综合亚洲| 熟女人妻精品中文字幕| 亚洲美女视频黄频| 国产激情欧美一区二区| 丝袜美腿在线中文| 日本免费一区二区三区高清不卡| www国产在线视频色| 欧美黄色片欧美黄色片| 又黄又粗又硬又大视频| 久久九九热精品免费| 久久国产乱子伦精品免费另类| 观看免费一级毛片| 在线观看美女被高潮喷水网站 | 婷婷精品国产亚洲av在线| 欧美中文综合在线视频| 九九热线精品视视频播放| 国产精品98久久久久久宅男小说| 三级毛片av免费| 99久久精品热视频| netflix在线观看网站| 亚洲精华国产精华精| 男人舔奶头视频| 尤物成人国产欧美一区二区三区| a级毛片a级免费在线| 欧美最新免费一区二区三区 | svipshipincom国产片| 亚洲av日韩精品久久久久久密| 国产黄片美女视频| x7x7x7水蜜桃| 在线免费观看的www视频| 国产成+人综合+亚洲专区| 青草久久国产| 欧美av亚洲av综合av国产av| 久久久久久久亚洲中文字幕 | 色综合欧美亚洲国产小说| 国产精品一区二区三区四区久久| 欧美日韩瑟瑟在线播放| 欧美日韩一级在线毛片| 婷婷六月久久综合丁香| 欧美日本视频| 欧美乱妇无乱码| 欧美在线黄色| 国产激情欧美一区二区| 亚洲五月天丁香| 此物有八面人人有两片| 在线播放国产精品三级| 国内毛片毛片毛片毛片毛片| 国内精品美女久久久久久| 成熟少妇高潮喷水视频| 亚洲av日韩精品久久久久久密| 在线a可以看的网站| 老司机午夜十八禁免费视频| 99热只有精品国产| 在线十欧美十亚洲十日本专区| 99精品久久久久人妻精品| 日本 av在线| 91在线观看av| 亚洲av免费高清在线观看| 国产一区二区在线观看日韩 | 岛国视频午夜一区免费看| 久久久久免费精品人妻一区二区| 男人的好看免费观看在线视频| 91av网一区二区| 国产激情欧美一区二区| 中国美女看黄片| 无遮挡黄片免费观看| 欧美成人免费av一区二区三区| 有码 亚洲区| 三级男女做爰猛烈吃奶摸视频| 国产69精品久久久久777片| 十八禁网站免费在线| 精品福利观看| 免费在线观看亚洲国产| 中文字幕人妻丝袜一区二区| 白带黄色成豆腐渣| 99久久精品国产亚洲精品| 中文字幕人妻丝袜一区二区| 白带黄色成豆腐渣| 97超级碰碰碰精品色视频在线观看| 欧美黄色片欧美黄色片| 婷婷亚洲欧美| 99久久精品国产亚洲精品| 亚洲av成人av| 免费av不卡在线播放| 国产精品亚洲av一区麻豆| 国产欧美日韩精品亚洲av| 亚洲精品一卡2卡三卡4卡5卡| 在线观看免费午夜福利视频| 91麻豆精品激情在线观看国产| 国产精品乱码一区二三区的特点| 一进一出抽搐动态| 97超视频在线观看视频| 在线十欧美十亚洲十日本专区| 在线免费观看不下载黄p国产 | 一区福利在线观看| 亚洲无线观看免费| 精品熟女少妇八av免费久了| 91麻豆精品激情在线观看国产| 成人高潮视频无遮挡免费网站| 18禁黄网站禁片午夜丰满| 欧美黄色片欧美黄色片| 亚洲色图av天堂| 欧美3d第一页|