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

    99久久综合免费| 十八禁网站网址无遮挡| 伊人亚洲综合成人网| 一区二区三区精品91| 亚洲av福利一区| 国产成人免费无遮挡视频| 久久精品久久久久久久性| 男女无遮挡免费网站观看| 黑丝袜美女国产一区| 久久久久精品久久久久真实原创| 亚洲 欧美一区二区三区| 国产精品国产三级国产专区5o| 国产精品香港三级国产av潘金莲 | 尾随美女入室| 日韩av在线免费看完整版不卡| 七月丁香在线播放| 精品少妇内射三级| 国产成人欧美| 美女高潮到喷水免费观看| 制服诱惑二区| 黑人欧美特级aaaaaa片| 国产精品久久久久久久久免| 国产在线免费精品| 国产麻豆69| 国产成人免费无遮挡视频| 欧美成人精品欧美一级黄| 亚洲,欧美精品.| 国产福利在线免费观看视频| 夫妻午夜视频| 日韩中文字幕欧美一区二区 | 国产一区有黄有色的免费视频| 国产高清不卡午夜福利| 少妇的丰满在线观看| 人妻人人澡人人爽人人| 精品亚洲成a人片在线观看| 国产成人精品在线电影| 久久久精品国产亚洲av高清涩受| 国语对白做爰xxxⅹ性视频网站| 狠狠精品人妻久久久久久综合| 国产探花极品一区二区| 亚洲av成人精品一二三区| 少妇猛男粗大的猛烈进出视频| 欧美 日韩 精品 国产| 在线精品无人区一区二区三| 国产极品粉嫩免费观看在线| 91在线精品国自产拍蜜月| 成人二区视频| 国产一区二区三区综合在线观看| 最新中文字幕久久久久| 午夜老司机福利剧场| 新久久久久国产一级毛片| 夫妻性生交免费视频一级片| 免费不卡的大黄色大毛片视频在线观看| 一级毛片我不卡| 新久久久久国产一级毛片| 99国产综合亚洲精品| 日本vs欧美在线观看视频| 精品国产国语对白av| 久久精品国产鲁丝片午夜精品| 色94色欧美一区二区| 少妇被粗大猛烈的视频| 一级毛片我不卡| 老女人水多毛片| av.在线天堂| 少妇的丰满在线观看| 亚洲四区av| 亚洲欧美色中文字幕在线| 肉色欧美久久久久久久蜜桃| www.精华液| 下体分泌物呈黄色| 欧美黄色片欧美黄色片| 青草久久国产| 久久久国产一区二区| 国产精品嫩草影院av在线观看| 久久久久国产精品人妻一区二区| 18+在线观看网站| 婷婷色麻豆天堂久久| 成人影院久久| 中文字幕制服av| 五月开心婷婷网| 国产熟女午夜一区二区三区| 欧美日韩综合久久久久久| 少妇人妻 视频| 嫩草影院入口| 久久国产精品男人的天堂亚洲| 在线天堂中文资源库| 亚洲精品aⅴ在线观看| 哪个播放器可以免费观看大片| 一区在线观看完整版| 最新中文字幕久久久久| 免费不卡的大黄色大毛片视频在线观看| 欧美国产精品一级二级三级| 日本av手机在线免费观看| 精品久久久精品久久久| 性色av一级| 日本91视频免费播放| 国产又色又爽无遮挡免| 五月开心婷婷网| 亚洲av日韩在线播放| 国产国语露脸激情在线看| 男女高潮啪啪啪动态图| 久久久久久久国产电影| 精品国产一区二区三区久久久樱花| 女人精品久久久久毛片| 久久久久国产精品人妻一区二区| 伦理电影大哥的女人| 亚洲国产av新网站| 国产成人aa在线观看| 国产精品不卡视频一区二区| 日本wwww免费看| 日韩免费高清中文字幕av| 亚洲第一av免费看| 欧美精品一区二区免费开放| 毛片一级片免费看久久久久| 国产淫语在线视频| 久久久亚洲精品成人影院| 色哟哟·www| 男男h啪啪无遮挡| 亚洲国产欧美日韩在线播放| 国产成人午夜福利电影在线观看| 永久网站在线| 极品人妻少妇av视频| 中文字幕av电影在线播放| 亚洲伊人色综图| 狠狠婷婷综合久久久久久88av| 啦啦啦在线观看免费高清www| 一二三四在线观看免费中文在| 久久亚洲国产成人精品v| 波多野结衣一区麻豆| 国产免费视频播放在线视频| 国产亚洲av片在线观看秒播厂| 满18在线观看网站| 丰满迷人的少妇在线观看| 丝袜在线中文字幕| 亚洲第一区二区三区不卡| √禁漫天堂资源中文www| 久久精品久久久久久久性| 女性生殖器流出的白浆| 高清av免费在线| 麻豆精品久久久久久蜜桃| 国产在线视频一区二区| 女人精品久久久久毛片| 激情视频va一区二区三区| 亚洲精品美女久久av网站| 久久久精品国产亚洲av高清涩受| 一级a爱视频在线免费观看| 亚洲,欧美精品.| 视频区图区小说| 丰满少妇做爰视频| 日韩电影二区| 边亲边吃奶的免费视频| 女人精品久久久久毛片| 少妇猛男粗大的猛烈进出视频| 国产精品一二三区在线看| 久久久久久久亚洲中文字幕| 午夜福利,免费看| 69精品国产乱码久久久| 99久国产av精品国产电影| 热99久久久久精品小说推荐| 一区二区av电影网| 国产成人aa在线观看| 十八禁网站网址无遮挡| 成年女人毛片免费观看观看9 | √禁漫天堂资源中文www| 999精品在线视频| 成人亚洲精品一区在线观看| 免费久久久久久久精品成人欧美视频| 飞空精品影院首页| 国产成人精品久久二区二区91 | 好男人视频免费观看在线| 午夜影院在线不卡| 波多野结衣一区麻豆| 又黄又粗又硬又大视频| 久久精品国产亚洲av涩爱| 美女国产视频在线观看| 亚洲一区中文字幕在线| 欧美激情极品国产一区二区三区| 国产男女内射视频| 欧美成人午夜免费资源| 超色免费av| 少妇猛男粗大的猛烈进出视频| 精品一区二区三卡| 观看av在线不卡| 免费在线观看视频国产中文字幕亚洲 | 欧美精品一区二区免费开放| 日韩免费高清中文字幕av| 国产成人午夜福利电影在线观看| 一本—道久久a久久精品蜜桃钙片| 最近的中文字幕免费完整| 亚洲中文av在线| 午夜激情久久久久久久| 99久久人妻综合| 国产免费福利视频在线观看| 卡戴珊不雅视频在线播放| 男人添女人高潮全过程视频| 99热国产这里只有精品6| 亚洲一区中文字幕在线| 成人亚洲欧美一区二区av| 亚洲精品国产色婷婷电影| 在线亚洲精品国产二区图片欧美| 男女啪啪激烈高潮av片| 久久国内精品自在自线图片| 久久久久久久久久久免费av| 久久人人爽av亚洲精品天堂| 亚洲三级黄色毛片| 久久久久久伊人网av| 国产精品免费视频内射| 天堂中文最新版在线下载| 国产精品一区二区在线不卡| 青春草视频在线免费观看| 日本午夜av视频| 在现免费观看毛片| 色网站视频免费| 中国国产av一级| 免费黄频网站在线观看国产| 免费人妻精品一区二区三区视频| 不卡av一区二区三区| 精品人妻熟女毛片av久久网站| 国产一区亚洲一区在线观看| 国产在线一区二区三区精| 久久 成人 亚洲| 久久狼人影院| 亚洲美女视频黄频| 精品国产超薄肉色丝袜足j| 一级毛片黄色毛片免费观看视频| 一本色道久久久久久精品综合| 国产在视频线精品| 丝袜在线中文字幕| 国产毛片在线视频| 大话2 男鬼变身卡| 各种免费的搞黄视频| 国产精品久久久久成人av| 精品亚洲成a人片在线观看| 久久久久久久久久久免费av| 免费播放大片免费观看视频在线观看| 人人妻人人添人人爽欧美一区卜| 日韩中文字幕欧美一区二区 | 国产在视频线精品| 91成人精品电影| 99国产综合亚洲精品| 99久久中文字幕三级久久日本| videossex国产| 国产精品久久久av美女十八| 国产熟女午夜一区二区三区| 日韩中文字幕欧美一区二区 | 90打野战视频偷拍视频| 青青草视频在线视频观看| 精品国产一区二区三区四区第35| 精品一区二区免费观看| 女人被躁到高潮嗷嗷叫费观| 在线观看免费日韩欧美大片| 最近的中文字幕免费完整| 叶爱在线成人免费视频播放| 少妇猛男粗大的猛烈进出视频| 18+在线观看网站| 免费不卡的大黄色大毛片视频在线观看| 丝袜在线中文字幕| 国产精品成人在线| 曰老女人黄片| 日本-黄色视频高清免费观看| 免费观看性生交大片5| 久久这里有精品视频免费| 男女啪啪激烈高潮av片| 毛片一级片免费看久久久久| 看免费成人av毛片| 免费不卡的大黄色大毛片视频在线观看| 亚洲一区中文字幕在线| 成年美女黄网站色视频大全免费| 丝袜美足系列| 国产一区二区激情短视频 | 亚洲av电影在线进入| 9色porny在线观看| 久久这里有精品视频免费| 99热网站在线观看| 女人被躁到高潮嗷嗷叫费观| 国产 一区精品| 欧美精品亚洲一区二区| 欧美xxⅹ黑人| av卡一久久| 免费在线观看完整版高清| 午夜激情久久久久久久| 国产成人精品一,二区| 国产麻豆69| 黄色 视频免费看| 午夜日本视频在线| 777久久人妻少妇嫩草av网站| 精品国产一区二区三区久久久樱花| 日本爱情动作片www.在线观看| 少妇的逼水好多| 亚洲av福利一区| 97在线人人人人妻| 精品一区二区免费观看| 久久精品国产亚洲av天美| 香蕉国产在线看| 欧美精品亚洲一区二区| av网站免费在线观看视频| 性色avwww在线观看| 免费大片黄手机在线观看| 91久久精品国产一区二区三区| 丝袜喷水一区| 色吧在线观看| videos熟女内射| 最黄视频免费看| 99热全是精品| freevideosex欧美| 欧美成人午夜精品| 日产精品乱码卡一卡2卡三| 下体分泌物呈黄色| 伊人亚洲综合成人网| 99久国产av精品国产电影| 亚洲人成77777在线视频| 人人妻人人添人人爽欧美一区卜| 在线精品无人区一区二区三| 午夜免费鲁丝| 日本黄色日本黄色录像| 青青草视频在线视频观看| 赤兔流量卡办理| 大码成人一级视频| 国产色婷婷99| 一级片'在线观看视频| 丁香六月天网| 日日撸夜夜添| 免费黄频网站在线观看国产| 一本—道久久a久久精品蜜桃钙片| 黑人巨大精品欧美一区二区蜜桃| 看非洲黑人一级黄片| 超碰97精品在线观看| freevideosex欧美| 国产日韩欧美视频二区| 国产有黄有色有爽视频| 一本久久精品| 美女脱内裤让男人舔精品视频| 国产成人免费无遮挡视频| 捣出白浆h1v1| 青青草视频在线视频观看| 日韩在线高清观看一区二区三区| 丰满迷人的少妇在线观看| 热99久久久久精品小说推荐| 制服人妻中文乱码| 欧美黄色片欧美黄色片| 日韩成人av中文字幕在线观看| 尾随美女入室| 日本av手机在线免费观看| 国产成人一区二区在线| 日本爱情动作片www.在线观看| 欧美+日韩+精品| 久久午夜综合久久蜜桃| 伦精品一区二区三区| 久久热在线av| 91精品三级在线观看| 久久综合国产亚洲精品| 捣出白浆h1v1| 国产精品久久久久久av不卡| 老司机影院毛片| 大片免费播放器 马上看| 看十八女毛片水多多多| 毛片一级片免费看久久久久| 人体艺术视频欧美日本| 99久国产av精品国产电影| 免费女性裸体啪啪无遮挡网站| 国产在视频线精品| 最近最新中文字幕大全免费视频 | 久久午夜福利片| 欧美最新免费一区二区三区| 一本久久精品| 高清黄色对白视频在线免费看| 我要看黄色一级片免费的| 日本免费在线观看一区| 热99国产精品久久久久久7| 黄频高清免费视频| 亚洲一区中文字幕在线| 亚洲一区二区三区欧美精品| 黄网站色视频无遮挡免费观看| 欧美亚洲 丝袜 人妻 在线| 97精品久久久久久久久久精品| 亚洲激情五月婷婷啪啪| 激情视频va一区二区三区| 黄片无遮挡物在线观看| 韩国av在线不卡| 欧美日韩成人在线一区二区| 国产一级毛片在线| 亚洲欧洲国产日韩| 国产亚洲av片在线观看秒播厂| 十分钟在线观看高清视频www| 超碰成人久久| 一级片免费观看大全| 日韩 亚洲 欧美在线| 最近最新中文字幕大全免费视频 | 18在线观看网站| 午夜日韩欧美国产| 美女视频免费永久观看网站| 免费女性裸体啪啪无遮挡网站| 天美传媒精品一区二区| 欧美在线黄色| 久久精品久久久久久噜噜老黄| 最近的中文字幕免费完整| 熟女电影av网| 精品少妇久久久久久888优播| 亚洲精品久久成人aⅴ小说| 日韩精品有码人妻一区| 少妇人妻精品综合一区二区| 熟女电影av网| 一区二区三区四区激情视频| 久久毛片免费看一区二区三区| 一本色道久久久久久精品综合| 18禁动态无遮挡网站| 亚洲精品日韩在线中文字幕| 婷婷色综合www| 女性生殖器流出的白浆| 成人18禁高潮啪啪吃奶动态图| 一区福利在线观看| 久久久久人妻精品一区果冻| 精品少妇黑人巨大在线播放| 嫩草影院入口| 亚洲国产精品一区二区三区在线| 一级黄片播放器| 91国产中文字幕| 麻豆精品久久久久久蜜桃| 国产 一区精品| 久久久久久人妻| 亚洲国产av影院在线观看| 国产野战对白在线观看| 色哟哟·www| 久久亚洲国产成人精品v| 校园人妻丝袜中文字幕| 久久久久久久久久人人人人人人| a级毛片黄视频| 久久免费观看电影| 美女视频免费永久观看网站| 色婷婷av一区二区三区视频| 国产色婷婷99| 亚洲国产精品一区三区| videosex国产| 老司机亚洲免费影院| 男的添女的下面高潮视频| 国产乱来视频区| 哪个播放器可以免费观看大片| 丝瓜视频免费看黄片| 久久久久视频综合| 免费看不卡的av| 国产成人av激情在线播放| www.av在线官网国产| 王馨瑶露胸无遮挡在线观看| 久久精品国产亚洲av高清一级| 97精品久久久久久久久久精品| 国产精品无大码| 日韩制服丝袜自拍偷拍| 精品午夜福利在线看| 99热国产这里只有精品6| 亚洲av在线观看美女高潮| 各种免费的搞黄视频| 日韩一区二区视频免费看| 国产精品一区二区在线不卡| 伦理电影大哥的女人| 又黄又粗又硬又大视频| 丝袜人妻中文字幕| 熟妇人妻不卡中文字幕| 欧美精品国产亚洲| 久久av网站| 少妇猛男粗大的猛烈进出视频| av.在线天堂| 一级a爱视频在线免费观看| 欧美日韩视频精品一区| av在线app专区| 国产激情久久老熟女| 人妻一区二区av| 丝袜美足系列| 男的添女的下面高潮视频| 妹子高潮喷水视频| 精品一区二区免费观看| 两个人免费观看高清视频| 人人澡人人妻人| 免费高清在线观看日韩| 日韩 亚洲 欧美在线| 一区在线观看完整版| 丁香六月天网| 丰满乱子伦码专区| 女人久久www免费人成看片| 不卡av一区二区三区| 亚洲精品久久久久久婷婷小说| 午夜精品国产一区二区电影| 亚洲综合色网址| 午夜福利,免费看| 久久午夜福利片| 男的添女的下面高潮视频| 国产一级毛片在线| 日韩中文字幕视频在线看片| 免费在线观看完整版高清| 精品卡一卡二卡四卡免费| 满18在线观看网站| 国产成人欧美| 大码成人一级视频| 亚洲久久久国产精品| 国语对白做爰xxxⅹ性视频网站| 久久免费观看电影| 精品99又大又爽又粗少妇毛片| 亚洲成国产人片在线观看| 热99久久久久精品小说推荐| 热re99久久国产66热| 欧美精品一区二区大全| 欧美日韩亚洲高清精品| 国产男人的电影天堂91| 亚洲久久久国产精品| 午夜福利影视在线免费观看| 欧美亚洲 丝袜 人妻 在线| 久久精品国产a三级三级三级| 国产在线视频一区二区| 一本大道久久a久久精品| 天天躁夜夜躁狠狠躁躁| 久久久国产一区二区| 精品人妻熟女毛片av久久网站| 一本大道久久a久久精品| 欧美日韩成人在线一区二区| 97在线人人人人妻| 制服人妻中文乱码| 熟女电影av网| 精品国产超薄肉色丝袜足j| av免费在线看不卡| 黄片小视频在线播放| 王馨瑶露胸无遮挡在线观看| 国产黄色免费在线视频| 十八禁高潮呻吟视频| 久久ye,这里只有精品| 午夜福利在线免费观看网站| 亚洲精品第二区| 日韩不卡一区二区三区视频在线| 精品人妻在线不人妻| 在线观看免费高清a一片| 人成视频在线观看免费观看| 狠狠精品人妻久久久久久综合| 亚洲精品国产一区二区精华液| 免费大片黄手机在线观看| 男人爽女人下面视频在线观看| 制服诱惑二区| 日韩视频在线欧美| 日日爽夜夜爽网站| 日韩人妻精品一区2区三区| 自线自在国产av| 欧美中文综合在线视频| 国产在线免费精品| 最近2019中文字幕mv第一页| 成人国产麻豆网| 国产av一区二区精品久久| 大片免费播放器 马上看| 亚洲av中文av极速乱| 精品亚洲乱码少妇综合久久| 久久久久人妻精品一区果冻| 中文字幕另类日韩欧美亚洲嫩草| 国产精品久久久久久精品古装| 久久久久国产一级毛片高清牌| 十八禁网站网址无遮挡| a级毛片黄视频| 又大又黄又爽视频免费| 久久这里有精品视频免费| 精品少妇内射三级| 国产极品天堂在线| 亚洲国产毛片av蜜桃av| 日本av手机在线免费观看| 男女下面插进去视频免费观看| 久久久久久人人人人人| 欧美另类一区| 叶爱在线成人免费视频播放| 婷婷色av中文字幕| 国产精品免费大片| 亚洲人成77777在线视频| 麻豆精品久久久久久蜜桃| 爱豆传媒免费全集在线观看| 亚洲国产日韩一区二区| 日韩 亚洲 欧美在线| 少妇人妻久久综合中文| 尾随美女入室| 五月天丁香电影| 日韩不卡一区二区三区视频在线| 国产不卡av网站在线观看| 婷婷色综合www| 五月伊人婷婷丁香| 午夜福利一区二区在线看| 亚洲欧美一区二区三区久久| 国产精品无大码| 日韩精品免费视频一区二区三区| 亚洲成人手机| 久久99蜜桃精品久久| 亚洲三级黄色毛片| 色播在线永久视频| 中文欧美无线码| 精品视频人人做人人爽| 少妇猛男粗大的猛烈进出视频| freevideosex欧美| 国产精品久久久av美女十八| 国产1区2区3区精品| 婷婷色麻豆天堂久久| 久久久久久久亚洲中文字幕| 啦啦啦视频在线资源免费观看| 一级片免费观看大全| 丁香六月天网| 一区二区三区四区激情视频| 一级,二级,三级黄色视频| 免费黄频网站在线观看国产| 国精品久久久久久国模美| 久久免费观看电影| 午夜福利视频在线观看免费| 91在线精品国自产拍蜜月| 1024视频免费在线观看| 亚洲国产毛片av蜜桃av| 国产1区2区3区精品| 免费人妻精品一区二区三区视频| 天天操日日干夜夜撸| 国产乱人偷精品视频| 美女xxoo啪啪120秒动态图| 亚洲精品aⅴ在线观看| 97在线人人人人妻| 啦啦啦中文免费视频观看日本| 亚洲av中文av极速乱| 久久国产精品男人的天堂亚洲| 少妇猛男粗大的猛烈进出视频|