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

    Void Fraction Distributions in Cold-gassed and Hot-sparged Three Phase Stirred Tanks with Multi-impeller*

    2009-05-15 01:40:16CHENLei陳雷BAOYuyun包雨云andGAOZhengming高正明
    關(guān)鍵詞:陳雷

    CHEN Lei (陳雷), BAO Yuyun (包雨云) and GAO Zhengming (高正明)

    ?

    Void Fraction Distributions in Cold-gassed and Hot-sparged Three Phase Stirred Tanks with Multi-impeller*

    CHEN Lei (陳雷), BAO Yuyun (包雨云) and GAO Zhengming (高正明)**

    School of Chemical Engineering, Beijing University of Chemical Technology, Beijing 100029, China

    Vertical distributions of void fraction in gas-liquid and gas-liquid-solid stirred tanks have been measured in a fully baffled dished base vessel of 0.48 m diameter, using a conductivity probe. The impeller configuration (a hollow half elliptical blade dispersing turbine below two up-pumping wide blade hydrofoils, identified as HEDT+2WHU) recommended in previous work has been used in this work. The operating temperatures were 24°C and 81°C, identified as cold and hot respectively. The effects of superficial gas velocity, agitator speed and the corresponding power input on the local void fraction in two-phase systems are investigated and discussed. Results show that the increasing of agitator speed or gas flow rate leads to an increase in local void fraction at the majority of measurement points in both cold and hot systems. However, the uniformity of gas dispersion does not always increase as the raising of agitator speed and power input. In either cold or hot sparged conditions, the two- and three-phase systems have similar vertical profiles for void fraction, with maxima in similar locations; however, the void fractions are significantly lower in hot sparging than with cold. In cold operation the presence of particles leads to a lower void fraction at most points, although the local void fractions increase a little with the addition of solid particles at high temperature, in good agreement with the global gas holdup results, and the possible reasons are discussed in this paper. This work can give a better understanding of the differences between cold-gassed and hot-sparged three phase stirred tanks.

    void fraction, conductivity probe, hot sparged reactor, dispersed gas distribution, multi-impeller

    1 INTRODUCTION

    Gas-liquid-solid stirred reactors are widely used in the chemical, mineral, and biochemical industries and wastewater treatment. Gas dispersion is not uniform among different regions of a stirred tank, even when stirred by a multi-impeller agitator. Bubble coalescence and breakup rates are believed to be different in various regions of the tank, resulting in different bubble sizes and gas-liquid mass transfer coefficients. However, a relatively uniform distribution of gas is desirable in many stirred reactors—such as fermenters in which aerobic bacteria normally require a certain range of dissolved oxygen concentration which is in turn closely related to the local void fraction. Complete analysis and optimized design cannot be realized if based only on the global information on performance properties. The distributions of void fraction throughout the tank would help to understand these differences in gas dispersion and mass transfer. On the other hand, development in computer technique and computational fluid dynamics (CFD) software encourages more and more researchers to focus on numerical simulation of gas dispersion in stirred tanks [1-5]. The validation of CFD investigations requires detailed information, such as the distribution of voids and solid particles, together with bubble sizes, all of which demand more investigation than simply the measurement of global properties.

    During the last several years there appear a few investigations of gas dispersion in industrial gas-liquid reactors [6-10]. The selection of measurement technique is an important part in these studies. The measurement techniques of local gas-phase characteristics in multi-phase systems can be classified as non-invasive and invasive techniques [11]. Among these techniques, the conductivity probe is a traditional and convenient one [8, 12]. Conductivity probes are invasive techniques, but they will not disturb the flow field seriously if the volume of the probe is kept small compared to the scale of the whole flow field. The conductivity probe is suitable for both opaque and dense dispersions, although the minimum detectable bubble size is limited by the size of the probe tip. Gao. [13] measured the vertical void fraction distribution in cold gassed, hot sparged and boiling gas-liquid systems with multiple impeller agitators, using the conductivity probe method. Two distinct peaks are observed in the void fraction profiles, and the local void fractions are always less in the hot sparged case than that in cold sparging. However, the effects of agitation speed and gas flow rate on void fraction distribution, which are of great importance in industrial design, are not included in their investigations.

    Three-phase systems are used in many industrial processes, including hydrogenation, oxidation, fermentation, evaporative crystallization and froth flotation. Extensive studies [14-17] of three phase systems have been reported during the last 30 years. All these investigations were carried out in ambient (.. “cold”) conditions where the vapor pressure of the liquid phase can be assumed to have little effect on the process. For many exothermic (.. “hot”) processes,.. oxidation, hydrogenation and polymerization, the vapor pressure should be taken into account. In our previous work [18-20], radical differences were reported between cold-gassed and hot-sparged stirred tanks. The absolute value of sparged relative power demand is generally higher in hot conditions than that in cold operation. As solid concentration is increased, there is little change in the total gas holdup at higher temperature, whereas the total gas holdup at ambient temperature decreases significantly. At a given total gas flow rate and power input, the total gas holdup in a hot sparged system is lower than that in cold conditions, but the difference is reduced from 50% to about 25% with increasing solid concentration. However, only these global measurements are not sufficient. The distribution of local void fraction in a three phase stirred tank is still an open question for both cold and hot-sparged systems.

    This paper reports the effects of agitation speed and gas flow rate on local void fractions in multi-impeller gas-liquid system at both low and high temperatures. The relationship between the uniformity of local void fraction and the operating parameters, including the agitation speed, gas rates and temperature, are also discussed on basis of the voidage distribution. The effects of solid particles on the distribution of local void fraction in both cold and hot systems are also investigated.

    2 EXPERIMENTAL

    2.1 Equipment

    Figure 1 Schematic of the experimental setup

    Figure 2 Impellers

    Four cylindrical heaters (each rated at 3 kW on full power) were mounted vertically in the base of the tank and located midway between two baffles. The vapor generated in the vessel was condensed in the condenser and returned to the stirred tank, maintaining the liquid inventory in the tank.

    Air and deionized water were used for all the experiments. In order to ensure sufficient conductivity in the deionized water, a little phosphoric acid was added. The solid particles were glass beads of 90 μm mean diameter and density 2500 kg·m-3at a volumetric concentration of 9%. Air was purified by passing through three-stage filters and then sparged into the tank. The total gas rates, including the vaporized water, ranged from 10 to 25 m3·h-1(with corresponding superficial velocities were 0.0156 to 0.039 m·s-1based on the tank cross section). Agitator speeds ranged from 4 to 8 s-1. The liquid bulk temperatures in cold and hot conditions were 24oC and 81oC respectively, maintained within ±0.2oC by a thermal resistance thermostat.

    2.2 Measurement technique

    In this work, the local void fraction was measured by using a needle point conductivity probe technique co-developed with Institute of Process Engineering of the Chinese Academy of Sciences. The hard stainless steel needle of 0.2 mm diameter, coated with insulating varnish, was used to serve as the sensor probe, as shown in Fig. 3. The varnish is hydrophobic, giving the probe a very short response time. A very short length was left uncoated at each needle tip for exposure to the experimental fluids as an electrode of the probe. The probe was supported by a stainless steel tube which serves as another electrode. In order to prevent the probe becoming polarized, an AC voltage was used.

    Figure 3 Construction of the electric conductivity probe

    The principle of this technique is based on the difference in electrical conductivity between liquid and gas phases. During the measurement, the bubbles are pierced by the probe placed at a given position in gas-liquid flow. When the probe is in the gas phase, the measurement circuit is open, and a high voltage output signal is produced. When the probe is in the liquid phase, the measurement circuit is closed, giving a low voltage output. After amplifier, rectifier and smoothing circuits, square-wave signals are generated, as shown in Fig. 4.

    Figure 4 Processed signals from the electric conductivity probe

    The gas holdup is defined as the volume fraction of gas phase. If it is assumed that the local flow is homogeneous, the probe tip has an equal probability of piercing any point on the projected frontal area of an approaching bubble and the bubble-probe intercept varies from zero to the bubble diameter, the local void fraction corresponds to the probability that the probe tip is in the gas phase. The local void fraction was estimated using the time averaged quantity,, given by

    wheregiis theth bubble passage time,bis the number of bubbles in the sample, andis the total sampling period.

    In this work a sampling frequency of 10 kHz was used to collect the data from the sensor probe. In every group used in the measurement more than 500 bubbles were detected to obtain statistically representative data.

    2.3 Saturation percentage of air in the hot sparged system

    Throughout this paper the total volumetric gas flow rate in a hot-sparged system has been calculated by Eq. (3), which has been used and discussed in our previous papers [19, 20]:

    where*is the various hot sparged “equilibrium temperature” ranging in this paper from 24 to 81oC,0is the temperature of the gas at inlet in centigrade,KandK0are “equilibrium temperatures” and inlet gas temperature in Kelvin, respectively.gis the inlet gas flow rate,0is the ambient pressure, andvis the saturated vapour pressure at*calculated by the Antoine equation. For the control temperatures of 24oC and 81oC in this work, the vapour pressure is 3 kPa and 50 kPa respectively.g+vis the total flow rate of the saturated gas.

    3 RESULTS AND DISCUSSION

    3.1 Effect of agitation speed on void fraction in cold and hot systems

    Figure 5 Measurement points arrangement

    /s-1: □?4; ○?5; △?6; ▽?7; ◇?8

    The comparison of data obtained in cold and hot systems at the same operating conditions [Fig. 6 (a). 6 (c) and Fig. 6 (b). 6 (d)] suggests that the gas void fraction at ambient temperature is generally higher than that at the same location but at higher temperature. This result is in good agreement with Gao. [13] who provided possible reasons to explain the difference of gas holdup between cold and hot systems. It should be noted that, as seen in Fig. 6, when the agitation speed is less than 5 s-1, the local gas void fraction near the bottom of the tank in a hot system is higher than that at ambient temperature.

    This may be explained by a difference in the size of cavities behind the turbine blades in hot and cold systems. When a lot of gas is introduced into the tank with agitation at a relatively low power input, the cavities are larger in cold conditions than in hot liquid [20]. The pumping capacity of the impellers is therefore lower in cold system at low agitation speed, so less gas is circulated to the bottom of the tank. The large cavities in the cold system are gradually dispersed at higher agitator speeds. The pumping capacity of the impellers is then enhanced and the local gas holdup near the bottom increases.

    /:■?0.85;□ 0.7

    The uniformity of gas dispersion in different regions of the stirred vessel cannot be evaluated from the overall gas holdup, though being of great importance in gas-liquid dispersion. In this work, the distributions of local void fraction in the bulk of the stirred tank were measured, providing the possibility of relating the uniformity of gas dispersion to different operating conditions. The dimensionless variation coefficient,, was used, as given in Eq. (4), to evaluate the uniformity of the gas dispersion.

    Figure 8 shows that given a constant superficial gas velocity, increasing the agitation speed at first leads to a rapid decrease in the variation coefficients of the void fraction distributions, but later the coefficients remain relatively constant and stable in both cold and hot systems. The higher value ofin a hot sparged system implies that the gas void fractions are less uniform than those in cold conditions at the same agitation speed or power input, as is seen in Fig. 8 (b). These stable variation coefficients mean that increasing the agitation speed or power input is not always helpful for the gas dispersion, especially with the gas superficial velocity of 0.031 m·s-1when the power input is already at a pretty high level such as 1.7 or 1.8 W·kg-1. Of course, this critical power input is closely related to the gas rate, and will decrease when less gas is introduced into the vessel: this will be discussed in the next section.

    3.2 Effect of gas rate on void fraction in both cold and hot systems

    Increasing the gas rate will increase the overall gas holdup. This is true of the local void fraction at most measurement points in the stirred tank as shown in Figs. 9 (a), 9 (b) and 9 (c), 9 (d) for both cold and hot conditions. However, the local void fraction near the bottom of the tank changes in a different way when gas rates are increased. It seems that the local void fraction there does not increase as much as in other regions of the tank, and even decreases a little in a cold system. On one hand, buoyancy will help bubbles to rise faster when more gas is introduced from the sparger, especially when the bottom impeller does not disperse the introduced gas well, which is exactly the case here.

    ■?cold systems;□?hot systems

    s/m·s-1:□?0.0156;○?0.023;△?0.031;▽?0.039

    ■?cold systems;□?hot systems

    Figure 10 also shows that the variation coefficientapproaches the stable value only at the lowest gas rate, which implies that this agitator speed (or power input) is not high enough to distribute higher gas rates. With a constant, relatively low, agitator speed, the cavities behind the turbine blades grow with an increased superficial gas velocity, weakening the pumping capacity of the impellers and also reducing the power input. When the superficial gas velocity is large enough, the gas dispersion gradually changes as the loaded impeller floods. In addition, gas distribution is even worse in a hot system than in cold, as reflected in the higher variation coefficientat similar gas rates.

    Looking at Fig. 8 (a), it can be seen that with a superficial gas velocity of 0.031 m·s-1, when the agitation speed is increased from 6 to 7 s-1the variation coefficientdecreases from 0.347 to 0.288 very close to the constant value. This can be defined as the critical agitation speed with reference to reach a limiting value of. When the superficial velocity is increased to 0.039 m·s-1, the critical speed increases to 8 s-1where the variation coefficientis 0.290. The corresponding power inputs are also shown in parentheses in Fig. 10. When the superficial velocity is as low as 0.0156 m·s-1and the power input is 1.27 W·kg-1, the value ofis about 0.3. In order to reach a similar gas dispersion state at the superficial velocities of 0.031 and 0.039 m·s-1the power levels have to be increased to 1.75 and 2.53 W·kg-1, respectively. This information will be helpful in determining the power demand needed for industrial applications.

    3.3 Effect of particles on local void fraction distribution in both cold and hot systems

    Vertical distributions of void fraction in cold-gassed and hot-sparged two phase systems are shown in Fig. 11. As mentioned before the profiles are very similar, with two distinct peaks evident above the top and bottom impeller planes. Given the same total gas flow rates the local void fractions are always less in the hot-sparged system than those in cold conditions. This is in good agreement with the results of Gao. [13] who included some discussion about the difference between hot and cold systems, though that paper did not consider the effect of solids.

    ■?cold systems;□?hot systems

    v: ■?9%;□?0%

    Our previous works [19, 20] have reported that the effects of solids on the overall gas holdup are quite different in cold and hot systems. The results showed that in a cold system the averaged gas holdup decreased significantly with an increase in solid concentration, but, in hot conditions, it was almost independent of the solid concentration, or even increased a little. Figs. 12 and 13 show the effect of solids on vertical void fraction distributions in cold and hot systems. It can be seen that particles at a volumetric concentration of 9% have little effect on the vertical void fraction distribution. Both two- and three-phase systems have similar void fraction profiles, with maxima at similar locations. Fig. 12 demonstrates that the presence of particles leads to a tendency of void fractions to be reduced at most measurement points in the cold system, which is similar to the results of our previous paper on overall gas holdup [19]. The probable reason for this was attributed to a change in bubble size between two and three phase conditions, a point which needs further study through the direct measurement of bubble size. The quantitative comparison for both local and overall gas hold-ups will be discussed later.

    v: ■?9%;□?0%

    Table 1 Comparison between overall gas holdup andvertical line average gas holdup

    The local void fractions are always less in hot conditions than when cold, as is the rule in two phase systems. However, the temperature effect on void fractions in three phase systems is less than that in two phase ones, which can be seen by comparing the data in Table 1. When the temperature is increased from 24 to 81°C, the mean void fraction in a three phase system decreases from 9.7% to 7.1%, about 36%, but it changes by more than 60%, from 10.6% to 6.6% in the absence of solids.

    4 CONCLUSIONS

    The effects of agitation speed, superficial gas velocity and particles on the local void fraction in cold- gassed (24oC) and hot-sparged (81oC) gas-liquid dispersion system have been investigated in an agitated vessel of 0.48 m diameter agitated by a multi-impeller combination. A point conductivity probe was used for all measurements.

    (1) With an increase in agitation speed or power input, the void fraction increases at most points in both cold and hot systems. With a fixed superficial gas velocity, an increased agitator speed leads to at first a more uniform gas void fraction distribution in the vessel, followed by a relatively stable level at relatively high power input levels when an increase in agitation speed or power input does not always help in distributing the gas. The gas void fractions in hot-sparged systems are less uniform than those in cold-gassed systems at similar agitation speeds or power input.

    (2) An increase in gas rate increases the local void fraction at most measurement points in both cold and hot systems. However, at a constant agitation speed, a higher sparge rate leads to a less uniform distribution and a greater power input is needed to distribute the gas more uniformly.

    (3) The presence of solid particles at a volumetric concentration of 9% has little effect on the vertical void fraction distributions. Both two- and three-phase systems have similar void fraction profiles, with maxima at similar locations. However, the presence of particles leads to tendency to decrease the void fraction at most locations in the cold system, against a trend of increase in hot conditions. The temperature effect on the void fractions in three-phase systems is not as great as that in two-phase ones.

    ACKNOWLEDGEMENTS

    The authors sincerely acknowledge the helpfuldiscussion with Prof. John M. Smith [Fluids and SystemsResearch Centre, School of Engineering (J2), Universityof Surrey, Guildford, GU2 7XH, UK].

    NOMENCLATURE

    variation coefficient

    diameter of impeller, m

    height of liquid in tank (general), m

    agitation speed, s-1

    bnumber of bubbles in the sample

    ginlet gas flow rate, m3·s-1

    g+vtotal gas flow rate including both air and vapour, m3·s-1

    radius of the tank, m

    distance between the measurement point and the center of the tank, m

    diameter of the tank, m

    0ambient temperature of the air,°C

    *equilibrium temperature in hot sparged system,°C

    total sampling period, s

    githeth bubble passage time, s

    ssuperficial gas velocity, m·s-1

    axial height of measurement points, m

    local void fraction

    standard deviation

    1 Wang, W.J., Mao, Z.S., “Numerical simulation of gas-liquid flow in a stirred tank with a Rushton impeller”,...., 10, 385-395 (2002).

    2 Lane, G.L., Schwarz, M.P., Evans, G.M., “Numerical modeling of gas-liquid flow in stirred tanks”,..., 60, 2203-2214 (2005).

    3 Khopkar, A.R., Kasat, G.R., Pandit, A.B., Ranade, V.V., “CFD simulation of mixing in tall gas-liquid stirred vessel: Role of local flow patterns”,..., 61, 2921-2929 (2006).

    4 Min, J., Bao, Y., Chen, L., Gao, Z., Smith, J.M., “Numerical simulation of gas dispersion in an aerated stirred reactor with multiple impellers”,...., 47, 7112-7117 (2008).

    5 Murthy, B.N., Kasundra, R.B., Joshi, J.B., “Hollow self-inducing impellers for gas-liquid-solid dispersion: Experimental and computational study”,..., 141, 332-345 (2008).

    6 Lu, W., Ju, S., “Local gas holdup, mean liquid velocity and turbulence in an aerated stirred tank using hot-film anemometry”,..., 35, 9-17 (1987).

    7 Gao, Z., Wang, Y., Shi, L., Fu, J., “Theoretical and experimental studies on bubble diameter and gas holdup in aerated stirred tanks”,...., 4 (4), 283-289 (1996).

    8 Barigou, M., Greaves, M., “Gas holdup and interfacial area distributions in a mechanically agitated gas-liquid contactor”,., 74 (A), 397-405 (1996).

    9 Laakkonen, M., Honkanen, M., Saarenrinne, P., Aittamaa, J., “Local bubble size distributions, gas-liquid interfacial areas and gas holdups in a stirred vessel with particle image velocimetry”,..., 109, 37-47 (2005).

    10 Gao, Z., Smith, J.M., Müller-Steinhagen, H., “Gas dispersion in sparged and boiling reactors”,...., 79, 973-978 (2001).

    11 Boyer, C., Duquenne, A.M., Wild, G., “Measuring techniques in gas-liquid and gas-liquid-solid reactors”,..., 57, 3185-3215 (2002).

    12 Sun, K., Zhang, M., Chen, X., “Local measurement of gas-liquid bubbly flow with a double-sensor probe”,...., 8 (1), 33-40 (2000).

    13 Gao, Z., Smith, J.M., Muller-Steinhagen, H., “Void fraction distribution in sparged and boiling reactors with modern impeller configuration”,..., 40, 489-497 (2001).

    14 Nienow, A.W., Bujalski, W., “Recent studies on agitated three-phase (gas-solid-liquid) systems in the turbulent regime”,., 80 (A), 832-838 (2002).

    15 Wen, J., Zhou, H., Chen, Y., “Local gas phase flow characteristics of a gas-liquid-solid three-phase reversed flow jet loop reactor”,...., 10 (1), 119-122 (2002).

    16 Bao, Y., Hao, Z., Gao, Z., Shi, L., Smith, J.M., “Suspension of buoyant particles in a three phase stirred tank”,..., 60, 2283-2292 (2005).

    17 Bao, Y., Hao, Z., Gao, Z., Shi, L., Smith, J.M., Thorpe, R.B., “Gas dispersion and solid suspension in a three-phase stirred tank with multiple impellers”,..., 193, 801-825 (2006).

    18 Gao, Z., Shi, L., “Effect of temperature on gas hold-up in aerated stirred tanks”,...., 11 (2), 204-207 (2003).

    19 Bao, Y.Y., Chen, L., Gao, Z.M., Zhang, X., Smith, J.M., Kirkby, N.F., “Temperature effects on gas dispersion and solid suspension in a three-phase stirred reactor”,...., 47, 4270-4277 (2008).

    20 Bao, Y.Y., Zhang, X., Gao, Z.M., Chen, L., Smith, J.M., Kirkby, N.F., “Gas dispersion and solid suspension in a hot sparged multi-impeller stirred tank”,...., 47, 2049-2055 (2008).

    2009-06-19,

    2009-10-16.

    the National Natural Science Foundation of China (20576009, 20821004) and the National Basic Research Program of China (2007CB714300).

    ** To whom correspondence should be addressed. E-mail: gaozm@mail.buct.edu.cn

    猜你喜歡
    陳雷
    Numerical simulation study of ionization characteristics of argon dielectric barrier discharge
    Detecting the possibility of a type of photon number splitting attack in decoy-state quantum key distribution
    間斷吸唾技術(shù)對(duì)根管治療患者舒適度的影響
    陳雷膠漆
    命中注定
    水利部部長(zhǎng)陳雷在我省調(diào)研水利工作
    山西水利(2016年5期)2017-01-20 08:51:22
    時(shí)尚“吃播”陳雷:邊吃邊秀邊掙錢
    陳雷主持召開國(guó)家防總會(huì)商會(huì)安排臺(tái)風(fēng)“天鵝”和新一輪強(qiáng)降雨防范工作
    陳雷講“三嚴(yán)三實(shí)”專題黨課
    陳雷主持召開全國(guó)水庫(kù)安全度汛視頻會(huì)議
    人妻久久中文字幕网| 美女cb高潮喷水在线观看 | 97碰自拍视频| 亚洲片人在线观看| 免费观看精品视频网站| 国产亚洲欧美98| 好男人在线观看高清免费视频| 露出奶头的视频| 免费看美女性在线毛片视频| 精品国产乱码久久久久久男人| a在线观看视频网站| 国产成人影院久久av| 久99久视频精品免费| 久久久国产成人精品二区| 欧美中文综合在线视频| 美女扒开内裤让男人捅视频| 丰满人妻熟妇乱又伦精品不卡| 男女那种视频在线观看| 日本一本二区三区精品| 色综合婷婷激情| 亚洲成av人片在线播放无| 亚洲欧美日韩东京热| 白带黄色成豆腐渣| 老司机深夜福利视频在线观看| 日韩欧美在线乱码| 欧美极品一区二区三区四区| 久久久成人免费电影| 国产高清视频在线观看网站| 19禁男女啪啪无遮挡网站| 中文字幕熟女人妻在线| 亚洲七黄色美女视频| 国产精品久久视频播放| 免费观看的影片在线观看| www.www免费av| 1024手机看黄色片| 一本一本综合久久| 狂野欧美白嫩少妇大欣赏| 夜夜夜夜夜久久久久| 熟女少妇亚洲综合色aaa.| 亚洲欧洲精品一区二区精品久久久| 久久国产精品人妻蜜桃| 国产探花在线观看一区二区| 国产探花在线观看一区二区| 欧美日本视频| 欧美乱色亚洲激情| 色老头精品视频在线观看| 九九在线视频观看精品| 黑人巨大精品欧美一区二区mp4| av福利片在线观看| 国产精品久久久久久久电影 | 三级男女做爰猛烈吃奶摸视频| 久久国产精品影院| 一区二区三区高清视频在线| 此物有八面人人有两片| 中文字幕人成人乱码亚洲影| 天天躁狠狠躁夜夜躁狠狠躁| 五月伊人婷婷丁香| 久久久成人免费电影| 亚洲av免费在线观看| 欧美日韩中文字幕国产精品一区二区三区| 亚洲av电影在线进入| 制服丝袜大香蕉在线| 免费看美女性在线毛片视频| 日韩av在线大香蕉| 俄罗斯特黄特色一大片| 免费大片18禁| 国产精品久久视频播放| 日本 欧美在线| 这个男人来自地球电影免费观看| 嫩草影院精品99| 亚洲av第一区精品v没综合| 成年女人永久免费观看视频| 久久久国产精品麻豆| 最近最新中文字幕大全免费视频| a在线观看视频网站| 熟女人妻精品中文字幕| 非洲黑人性xxxx精品又粗又长| 少妇的丰满在线观看| 国内久久婷婷六月综合欲色啪| 一级毛片高清免费大全| 亚洲欧美日韩东京热| 亚洲中文字幕日韩| 18禁观看日本| 欧美又色又爽又黄视频| 婷婷六月久久综合丁香| 深夜精品福利| 国产精品女同一区二区软件 | 又爽又黄无遮挡网站| 高清在线国产一区| 观看免费一级毛片| 中文资源天堂在线| 亚洲国产精品合色在线| 午夜视频精品福利| 天天躁狠狠躁夜夜躁狠狠躁| 国产一区二区激情短视频| 叶爱在线成人免费视频播放| 久久性视频一级片| 国产日本99.免费观看| 亚洲人成网站高清观看| 搡老熟女国产l中国老女人| 精品福利观看| 1000部很黄的大片| 一区福利在线观看| 欧美一区二区国产精品久久精品| 岛国在线免费视频观看| 在线观看日韩欧美| 成人三级做爰电影| 九色国产91popny在线| 国产成人精品久久二区二区免费| 又黄又爽又免费观看的视频| 动漫黄色视频在线观看| 亚洲真实伦在线观看| 在线观看一区二区三区| 亚洲国产欧洲综合997久久,| 日韩欧美免费精品| 国产亚洲精品久久久com| 国产淫片久久久久久久久 | 999精品在线视频| 亚洲真实伦在线观看| 少妇丰满av| 欧美黑人欧美精品刺激| 男女床上黄色一级片免费看| 成在线人永久免费视频| 国产成人影院久久av| 99久久精品国产亚洲精品| 伦理电影免费视频| 在线国产一区二区在线| 一本久久中文字幕| 99久久无色码亚洲精品果冻| bbb黄色大片| 校园春色视频在线观看| 人妻丰满熟妇av一区二区三区| 日本熟妇午夜| 日本黄大片高清| 九九久久精品国产亚洲av麻豆 | 熟女人妻精品中文字幕| svipshipincom国产片| 国产亚洲精品av在线| 国产精品1区2区在线观看.| 精品电影一区二区在线| 亚洲自拍偷在线| 中文字幕人妻丝袜一区二区| 51午夜福利影视在线观看| 国产真实乱freesex| 午夜成年电影在线免费观看| 日日摸夜夜添夜夜添小说| 国产精品国产高清国产av| 一进一出抽搐gif免费好疼| 中文字幕人成人乱码亚洲影| 哪里可以看免费的av片| 狠狠狠狠99中文字幕| 欧美日韩福利视频一区二区| av片东京热男人的天堂| 精品久久久久久成人av| 日日夜夜操网爽| 亚洲欧美精品综合一区二区三区| 99国产极品粉嫩在线观看| 丰满人妻一区二区三区视频av | 国产不卡一卡二| 成人特级av手机在线观看| 18美女黄网站色大片免费观看| 久99久视频精品免费| 精品久久久久久久久久久久久| 一区二区三区国产精品乱码| 曰老女人黄片| 国产三级中文精品| 久久精品夜夜夜夜夜久久蜜豆| 中文字幕熟女人妻在线| 这个男人来自地球电影免费观看| 免费看日本二区| 中文字幕久久专区| 麻豆成人av在线观看| 一级a爱片免费观看的视频| 国产精品九九99| 免费看十八禁软件| 亚洲欧美精品综合久久99| 欧美高清成人免费视频www| 亚洲 欧美 日韩 在线 免费| 精品久久久久久,| 久久久久久大精品| 成人av在线播放网站| 两个人看的免费小视频| 全区人妻精品视频| 丰满人妻一区二区三区视频av | 久久久水蜜桃国产精品网| 后天国语完整版免费观看| 少妇的丰满在线观看| 国产在线精品亚洲第一网站| 国产免费男女视频| 97碰自拍视频| 久久精品国产清高在天天线| 欧美激情在线99| 欧美日韩福利视频一区二区| 久久久精品大字幕| 欧美日韩综合久久久久久 | 午夜视频精品福利| 久久亚洲精品不卡| 91av网站免费观看| tocl精华| 欧洲精品卡2卡3卡4卡5卡区| 亚洲最大成人中文| 久久久国产欧美日韩av| 中文字幕熟女人妻在线| 国产成人aa在线观看| 国产成人av激情在线播放| 麻豆成人av在线观看| 一级a爱片免费观看的视频| 国内毛片毛片毛片毛片毛片| 亚洲成人久久爱视频| 一级毛片女人18水好多| 夜夜看夜夜爽夜夜摸| 国产精品,欧美在线| 免费人成视频x8x8入口观看| 国产亚洲欧美98| 国产熟女xx| 国产精品影院久久| 女生性感内裤真人,穿戴方法视频| 欧美成狂野欧美在线观看| 精品久久久久久成人av| 757午夜福利合集在线观看| 啦啦啦免费观看视频1| 黄色视频,在线免费观看| 99精品久久久久人妻精品| 男女那种视频在线观看| 久久香蕉精品热| 婷婷六月久久综合丁香| 日韩av在线大香蕉| 亚洲男人的天堂狠狠| 国产精品1区2区在线观看.| 婷婷六月久久综合丁香| 99久久无色码亚洲精品果冻| 亚洲精品456在线播放app | 夜夜看夜夜爽夜夜摸| 久久精品国产清高在天天线| 国产成人影院久久av| 亚洲精品乱码久久久v下载方式 | 亚洲欧美日韩东京热| 亚洲自拍偷在线| 亚洲性夜色夜夜综合| 国产激情欧美一区二区| 波多野结衣巨乳人妻| 久久精品人妻少妇| 亚洲国产欧美人成| 免费观看的影片在线观看| 制服人妻中文乱码| 一二三四在线观看免费中文在| 啦啦啦免费观看视频1| 亚洲av第一区精品v没综合| 成年免费大片在线观看| 欧美乱码精品一区二区三区| 日韩精品青青久久久久久| 亚洲国产欧美人成| bbb黄色大片| 免费一级毛片在线播放高清视频| 男人舔女人的私密视频| 免费看美女性在线毛片视频| 91麻豆av在线| 小说图片视频综合网站| 免费av毛片视频| 久久九九热精品免费| 亚洲在线自拍视频| 身体一侧抽搐| 97超视频在线观看视频| www.www免费av| av黄色大香蕉| 99视频精品全部免费 在线 | 美女扒开内裤让男人捅视频| 小说图片视频综合网站| 精品国产美女av久久久久小说| 母亲3免费完整高清在线观看| 中文在线观看免费www的网站| 一区二区三区激情视频| 久久天堂一区二区三区四区| 久久久久国内视频| 制服人妻中文乱码| 国产精品98久久久久久宅男小说| 精华霜和精华液先用哪个| cao死你这个sao货| 精品久久久久久成人av| 国产精品一区二区三区四区免费观看 | 男女午夜视频在线观看| 久久香蕉国产精品| 久久午夜亚洲精品久久| 国产伦精品一区二区三区四那| 欧美一区二区精品小视频在线| 综合色av麻豆| 免费在线观看成人毛片| 熟妇人妻久久中文字幕3abv| 日韩 欧美 亚洲 中文字幕| 免费人成视频x8x8入口观看| www.熟女人妻精品国产| 小说图片视频综合网站| 91麻豆精品激情在线观看国产| 亚洲成av人片免费观看| 免费在线观看日本一区| 我的老师免费观看完整版| 床上黄色一级片| 午夜激情福利司机影院| 亚洲午夜理论影院| 美女午夜性视频免费| 又大又爽又粗| 亚洲欧洲精品一区二区精品久久久| 熟女电影av网| 日韩三级视频一区二区三区| 免费在线观看成人毛片| 国产成人精品久久二区二区免费| 午夜免费成人在线视频| 黄频高清免费视频| 亚洲avbb在线观看| 动漫黄色视频在线观看| 国产午夜精品久久久久久| aaaaa片日本免费| 757午夜福利合集在线观看| 国内揄拍国产精品人妻在线| 高清毛片免费观看视频网站| 夜夜躁狠狠躁天天躁| 国产精品综合久久久久久久免费| 19禁男女啪啪无遮挡网站| 精品无人区乱码1区二区| 日韩国内少妇激情av| 国产一区二区在线观看日韩 | 久久久久久久久久黄片| 国产亚洲精品久久久com| 一级毛片高清免费大全| 无遮挡黄片免费观看| 久久亚洲真实| 一级毛片女人18水好多| 亚洲成人精品中文字幕电影| 中亚洲国语对白在线视频| 脱女人内裤的视频| 97碰自拍视频| 少妇裸体淫交视频免费看高清| 精品午夜福利视频在线观看一区| h日本视频在线播放| 波多野结衣高清作品| 国产美女午夜福利| 欧美绝顶高潮抽搐喷水| 在线观看一区二区三区| www.www免费av| 一区二区三区国产精品乱码| 亚洲自拍偷在线| 日本 欧美在线| 免费观看精品视频网站| 999久久久国产精品视频| 偷拍熟女少妇极品色| 国产成人啪精品午夜网站| 成年女人看的毛片在线观看| 欧美日韩中文字幕国产精品一区二区三区| 亚洲成a人片在线一区二区| 国产毛片a区久久久久| 一二三四社区在线视频社区8| x7x7x7水蜜桃| 国产精品99久久久久久久久| avwww免费| 久久久色成人| 亚洲国产欧美网| 天天躁狠狠躁夜夜躁狠狠躁| 嫁个100分男人电影在线观看| 美女被艹到高潮喷水动态| 国产三级中文精品| 国产成人精品无人区| 99久久成人亚洲精品观看| 国产三级中文精品| 中亚洲国语对白在线视频| 一级毛片精品| 成人三级黄色视频| 日韩高清综合在线| 日韩欧美在线乱码| 久久亚洲精品不卡| 99国产精品一区二区蜜桃av| 国产亚洲欧美98| 一级黄色大片毛片| 女人高潮潮喷娇喘18禁视频| 极品教师在线免费播放| 制服人妻中文乱码| 日本五十路高清| 午夜福利在线观看吧| 看片在线看免费视频| 国产精品香港三级国产av潘金莲| 国产欧美日韩一区二区精品| 91在线观看av| 国产不卡一卡二| 日本一二三区视频观看| 亚洲av片天天在线观看| 最新中文字幕久久久久 | 99热精品在线国产| 男女下面进入的视频免费午夜| 黄色片一级片一级黄色片| 757午夜福利合集在线观看| 亚洲精品美女久久久久99蜜臀| 欧美黄色淫秽网站| 亚洲精品在线美女| 黄频高清免费视频| 亚洲av成人精品一区久久| 亚洲自拍偷在线| 中国美女看黄片| 老司机深夜福利视频在线观看| 老汉色av国产亚洲站长工具| 又爽又黄无遮挡网站| 黑人巨大精品欧美一区二区mp4| 色综合站精品国产| 国产精品久久电影中文字幕| 日本黄色视频三级网站网址| 成年女人永久免费观看视频| 久久久久国产一级毛片高清牌| 精品一区二区三区视频在线观看免费| 久久久国产精品麻豆| 欧美另类亚洲清纯唯美| 成人av一区二区三区在线看| а√天堂www在线а√下载| 国内精品久久久久精免费| 国产精品香港三级国产av潘金莲| 亚洲国产色片| www日本在线高清视频| www.999成人在线观看| 成年女人毛片免费观看观看9| 精品不卡国产一区二区三区| 欧美一区二区国产精品久久精品| 欧美日韩亚洲国产一区二区在线观看| 亚洲一区二区三区色噜噜| 性色av乱码一区二区三区2| 国产高清三级在线| 少妇人妻一区二区三区视频| 国产精华一区二区三区| 婷婷亚洲欧美| 亚洲精品色激情综合| 久久久久性生活片| 日本成人三级电影网站| 亚洲av成人av| 色综合欧美亚洲国产小说| 成人av一区二区三区在线看| 亚洲aⅴ乱码一区二区在线播放| 国产1区2区3区精品| 国产高清激情床上av| 日韩有码中文字幕| 日本成人三级电影网站| 观看美女的网站| 啦啦啦观看免费观看视频高清| 热99在线观看视频| 久久久久性生活片| 两个人的视频大全免费| 啪啪无遮挡十八禁网站| 搡老熟女国产l中国老女人| 国产成年人精品一区二区| 婷婷丁香在线五月| 又粗又爽又猛毛片免费看| 日本一二三区视频观看| 精品福利观看| 亚洲欧美一区二区三区黑人| 男女那种视频在线观看| 国产精品香港三级国产av潘金莲| 无人区码免费观看不卡| 国产精品1区2区在线观看.| 国产精品av久久久久免费| 亚洲九九香蕉| 成人午夜高清在线视频| 少妇的逼水好多| www.999成人在线观看| www.熟女人妻精品国产| 日韩欧美精品v在线| 亚洲人成伊人成综合网2020| 精品国产超薄肉色丝袜足j| 母亲3免费完整高清在线观看| 搡老熟女国产l中国老女人| 哪里可以看免费的av片| 亚洲天堂国产精品一区在线| www日本黄色视频网| 久久久久性生活片| 久久精品人妻少妇| 久久午夜亚洲精品久久| 变态另类成人亚洲欧美熟女| 亚洲黑人精品在线| 熟妇人妻久久中文字幕3abv| 亚洲无线在线观看| 夜夜躁狠狠躁天天躁| 午夜免费观看网址| 法律面前人人平等表现在哪些方面| 天天躁狠狠躁夜夜躁狠狠躁| 免费在线观看成人毛片| 嫩草影视91久久| 麻豆成人午夜福利视频| 免费av不卡在线播放| 亚洲七黄色美女视频| 在线观看舔阴道视频| 国产一区二区激情短视频| 一级毛片女人18水好多| 欧美黄色片欧美黄色片| 伦理电影免费视频| av国产免费在线观看| 老司机午夜福利在线观看视频| 男人舔奶头视频| 亚洲欧美日韩高清专用| 国产亚洲精品久久久久久毛片| 精品久久久久久久末码| 99国产精品一区二区三区| 亚洲九九香蕉| 一本综合久久免费| 亚洲av电影在线进入| 日韩欧美在线乱码| 国产精品九九99| 啦啦啦韩国在线观看视频| 岛国视频午夜一区免费看| 日韩欧美国产在线观看| 国产免费av片在线观看野外av| www.精华液| 一区福利在线观看| 亚洲av电影不卡..在线观看| 国产成人系列免费观看| 欧美黄色淫秽网站| 亚洲乱码一区二区免费版| 97超级碰碰碰精品色视频在线观看| 国产精品 欧美亚洲| 美女大奶头视频| 国产亚洲欧美98| 国产三级中文精品| 亚洲人成电影免费在线| 动漫黄色视频在线观看| a在线观看视频网站| 在线观看免费视频日本深夜| 久久久久国内视频| 一二三四在线观看免费中文在| 最近最新中文字幕大全免费视频| 亚洲自偷自拍图片 自拍| 一个人免费在线观看电影 | 97人妻精品一区二区三区麻豆| e午夜精品久久久久久久| 婷婷精品国产亚洲av在线| 99热6这里只有精品| 18禁观看日本| bbb黄色大片| 久久久久亚洲av毛片大全| 久久久久性生活片| 亚洲自拍偷在线| 亚洲精品456在线播放app | 岛国视频午夜一区免费看| 日韩欧美国产在线观看| 国产蜜桃级精品一区二区三区| 亚洲国产精品久久男人天堂| 亚洲片人在线观看| 国产精品 欧美亚洲| 久久亚洲精品不卡| 久久伊人香网站| 婷婷六月久久综合丁香| 精品99又大又爽又粗少妇毛片 | 哪里可以看免费的av片| 午夜精品一区二区三区免费看| 香蕉av资源在线| 午夜久久久久精精品| 色吧在线观看| 夜夜躁狠狠躁天天躁| 变态另类丝袜制服| 99久久精品一区二区三区| 久久香蕉国产精品| 在线国产一区二区在线| 欧美性猛交╳xxx乱大交人| 一级毛片高清免费大全| 亚洲av日韩精品久久久久久密| 国产亚洲精品综合一区在线观看| 最近视频中文字幕2019在线8| 亚洲国产欧美网| 丰满人妻一区二区三区视频av | 亚洲成人久久爱视频| 好男人电影高清在线观看| 99久久国产精品久久久| 亚洲精品一区av在线观看| 国内精品久久久久久久电影| 午夜免费成人在线视频| 女人被狂操c到高潮| 亚洲自拍偷在线| 亚洲avbb在线观看| 一进一出好大好爽视频| 制服人妻中文乱码| 岛国视频午夜一区免费看| 成人一区二区视频在线观看| 国产精品九九99| 91九色精品人成在线观看| 欧美日韩乱码在线| 亚洲va日本ⅴa欧美va伊人久久| 69av精品久久久久久| 精品不卡国产一区二区三区| 国产精华一区二区三区| 每晚都被弄得嗷嗷叫到高潮| 99热6这里只有精品| 丁香欧美五月| 校园春色视频在线观看| 男女做爰动态图高潮gif福利片| 一本精品99久久精品77| 日韩三级视频一区二区三区| 国产亚洲精品综合一区在线观看| 午夜亚洲福利在线播放| 成人国产综合亚洲| 日本一二三区视频观看| 国产午夜福利久久久久久| 在线免费观看不下载黄p国产 | 19禁男女啪啪无遮挡网站| 亚洲专区字幕在线| av视频在线观看入口| 亚洲成av人片在线播放无| 免费大片18禁| 国内精品一区二区在线观看| 亚洲国产欧美一区二区综合| 国产一区二区三区视频了| 美女 人体艺术 gogo| av国产免费在线观看| 午夜久久久久精精品| 亚洲熟妇中文字幕五十中出| 日韩欧美在线乱码| 欧美激情在线99| 99视频精品全部免费 在线 | 亚洲 欧美 日韩 在线 免费| 成人av在线播放网站| 少妇人妻一区二区三区视频| 精品99又大又爽又粗少妇毛片 | 一级毛片精品| 色综合亚洲欧美另类图片|