• <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ì)議
    日本黄色日本黄色录像| 大又大粗又爽又黄少妇毛片口| 婷婷成人精品国产| 国产一区二区三区av在线| 在线观看免费日韩欧美大片 | 午夜影院在线不卡| 最近中文字幕2019免费版| 国产一区二区在线观看日韩| a级毛片免费高清观看在线播放| 久久综合国产亚洲精品| 久久青草综合色| 18禁在线播放成人免费| videossex国产| 午夜激情av网站| 亚洲五月色婷婷综合| 婷婷色综合www| 最后的刺客免费高清国语| 欧美xxxx性猛交bbbb| 在线免费观看不下载黄p国产| 亚洲精品日韩av片在线观看| 我的老师免费观看完整版| 精品国产一区二区三区久久久樱花| 精品国产露脸久久av麻豆| 免费观看在线日韩| 午夜福利网站1000一区二区三区| 18在线观看网站| 亚州av有码| 亚洲高清免费不卡视频| 国产亚洲精品第一综合不卡 | 一级,二级,三级黄色视频| 99久久中文字幕三级久久日本| 99久久综合免费| 18禁动态无遮挡网站| 看免费成人av毛片| 亚洲国产av影院在线观看| 卡戴珊不雅视频在线播放| 国产一级毛片在线| 国产精品.久久久| 一区在线观看完整版| 欧美性感艳星| 亚洲精品456在线播放app| 久久亚洲国产成人精品v| 两个人的视频大全免费| 久久人人爽人人爽人人片va| 亚洲欧美一区二区三区国产| 久久久久久人妻| 日本av手机在线免费观看| 亚洲精品乱码久久久久久按摩| 日韩人妻高清精品专区| 人成视频在线观看免费观看| 天堂中文最新版在线下载| 久久久久久久久久久免费av| 国产熟女欧美一区二区| 色网站视频免费| 插逼视频在线观看| 国产精品一区二区在线观看99| 黄色一级大片看看| 五月开心婷婷网| 成人免费观看视频高清| 国产免费一级a男人的天堂| 免费观看a级毛片全部| 午夜激情久久久久久久| 狂野欧美激情性xxxx在线观看| 啦啦啦视频在线资源免费观看| 18禁在线无遮挡免费观看视频| 2021少妇久久久久久久久久久| 国产精品熟女久久久久浪| 丰满饥渴人妻一区二区三| av线在线观看网站| 精品人妻在线不人妻| 女人久久www免费人成看片| 美女大奶头黄色视频| 黑人高潮一二区| 99久久精品国产国产毛片| av在线老鸭窝| 国产极品天堂在线| av国产精品久久久久影院| 美女视频免费永久观看网站| 亚洲经典国产精华液单| 日韩av免费高清视频| 欧美精品亚洲一区二区| 国产精品99久久久久久久久| 亚洲国产av新网站| 亚洲四区av| 在线天堂最新版资源| 少妇人妻久久综合中文| 亚洲高清免费不卡视频| 久久国产精品男人的天堂亚洲 | av播播在线观看一区| 午夜免费鲁丝| 美女视频免费永久观看网站| 精品国产一区二区三区久久久樱花| 18禁裸乳无遮挡动漫免费视频| 免费观看的影片在线观看| 涩涩av久久男人的天堂| 日韩av在线免费看完整版不卡| 国产片内射在线| 日韩一区二区三区影片| 人人妻人人澡人人爽人人夜夜| 国产精品秋霞免费鲁丝片| 毛片一级片免费看久久久久| 久久国产精品男人的天堂亚洲 | 有码 亚洲区| 熟女人妻精品中文字幕| 欧美激情国产日韩精品一区| 日本vs欧美在线观看视频| 国产在视频线精品| 日韩一区二区视频免费看| 男女啪啪激烈高潮av片| 卡戴珊不雅视频在线播放| 涩涩av久久男人的天堂| 18在线观看网站| xxxhd国产人妻xxx| 免费黄网站久久成人精品| 一区二区日韩欧美中文字幕 | 亚洲五月色婷婷综合| 久久精品久久久久久久性| 国产深夜福利视频在线观看| 国产日韩欧美亚洲二区| av在线老鸭窝| kizo精华| 肉色欧美久久久久久久蜜桃| 制服诱惑二区| 国产亚洲欧美精品永久| 亚洲婷婷狠狠爱综合网| 精品人妻一区二区三区麻豆| 日本黄大片高清| 国产视频首页在线观看| 九九爱精品视频在线观看| 99国产综合亚洲精品| 日韩三级伦理在线观看| 亚洲人与动物交配视频| 性高湖久久久久久久久免费观看| 日韩三级伦理在线观看| 香蕉精品网在线| 国产成人精品在线电影| 自拍欧美九色日韩亚洲蝌蚪91| 最近的中文字幕免费完整| 人人妻人人爽人人添夜夜欢视频| 成人国语在线视频| 亚洲精品456在线播放app| 成人亚洲欧美一区二区av| 午夜激情久久久久久久| 久久久久久久国产电影| 久久久久久久久久久久大奶| 免费观看无遮挡的男女| 成人午夜精彩视频在线观看| 亚洲美女视频黄频| 最近的中文字幕免费完整| 不卡视频在线观看欧美| 大片电影免费在线观看免费| 国产爽快片一区二区三区| 国产有黄有色有爽视频| 久久午夜综合久久蜜桃| 91成人精品电影| 男的添女的下面高潮视频| 青青草视频在线视频观看| 亚洲人成网站在线播| 大话2 男鬼变身卡| 啦啦啦中文免费视频观看日本| 亚洲国产精品国产精品| 美女大奶头黄色视频| 精品少妇久久久久久888优播| 中文字幕人妻熟人妻熟丝袜美| 天天操日日干夜夜撸| 丰满迷人的少妇在线观看| 国产爽快片一区二区三区| 肉色欧美久久久久久久蜜桃| 亚洲精品久久午夜乱码| av在线app专区| 日韩一本色道免费dvd| 一区二区三区四区激情视频| 在线观看免费视频网站a站| 国产精品无大码| 18禁在线无遮挡免费观看视频| 精品视频人人做人人爽| 国产伦精品一区二区三区视频9| 亚洲国产精品成人久久小说| 看十八女毛片水多多多| 国产视频首页在线观看| videos熟女内射| 91午夜精品亚洲一区二区三区| 久久午夜综合久久蜜桃| 天堂俺去俺来也www色官网| 国产精品麻豆人妻色哟哟久久| av视频免费观看在线观看| 久久久久久久久久人人人人人人| 交换朋友夫妻互换小说| 日韩中字成人| 极品人妻少妇av视频| 欧美另类一区| 亚洲婷婷狠狠爱综合网| 欧美日韩av久久| 自线自在国产av| 精品一区二区三卡| 国产高清不卡午夜福利| 国产精品成人在线| av福利片在线| 人人妻人人添人人爽欧美一区卜| 亚洲av综合色区一区| 人妻一区二区av| 成人手机av| 国产视频首页在线观看| 久久久久精品性色| 日韩视频在线欧美| 国产精品久久久久久久电影| 亚洲精品日韩av片在线观看| 一级二级三级毛片免费看| 下体分泌物呈黄色| 亚洲人成网站在线播| 精品人妻熟女毛片av久久网站| 亚洲成色77777| 日韩人妻高清精品专区| 国产精品人妻久久久影院| av电影中文网址| 精品亚洲成a人片在线观看| 九九在线视频观看精品| 有码 亚洲区| 最近手机中文字幕大全| 欧美激情 高清一区二区三区| 国产成人免费观看mmmm| 嘟嘟电影网在线观看| 你懂的网址亚洲精品在线观看| 在线观看国产h片| 伊人久久国产一区二区| 国产 一区精品| av不卡在线播放| 欧美+日韩+精品| 欧美精品国产亚洲| 极品少妇高潮喷水抽搐| 亚洲天堂av无毛| 日韩亚洲欧美综合| 国产精品.久久久| 国产片内射在线| 黑人巨大精品欧美一区二区蜜桃 | 中国三级夫妇交换| 国产精品一二三区在线看| 国产av码专区亚洲av| 久久ye,这里只有精品| 亚洲av在线观看美女高潮| 午夜日本视频在线| 最后的刺客免费高清国语| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 伊人久久精品亚洲午夜| 在线观看国产h片| 国产一区有黄有色的免费视频| 国产熟女午夜一区二区三区 | 国产精品无大码| 色婷婷久久久亚洲欧美| 成人亚洲欧美一区二区av| 免费大片18禁| 国产不卡av网站在线观看| 日日撸夜夜添| 乱人伦中国视频| 国产精品国产三级专区第一集| 亚洲四区av| 女性生殖器流出的白浆| 啦啦啦中文免费视频观看日本| 久久这里有精品视频免费| 亚洲美女视频黄频| 国产欧美日韩一区二区三区在线 | 久久久久久人妻| 国产极品天堂在线| 老司机影院成人| 波野结衣二区三区在线| 亚洲,一卡二卡三卡| 一区二区日韩欧美中文字幕 | 在线天堂最新版资源| 免费观看a级毛片全部| 综合色丁香网| 国产精品久久久久久精品古装| 少妇被粗大的猛进出69影院 | 少妇精品久久久久久久| 亚洲在久久综合| 街头女战士在线观看网站| 国产男女内射视频| 下体分泌物呈黄色| 老熟女久久久| 久久女婷五月综合色啪小说| 成人国产麻豆网| 精品卡一卡二卡四卡免费| 国产欧美日韩一区二区三区在线 | 国产片内射在线| 亚洲精品456在线播放app| 99国产综合亚洲精品| 97在线视频观看| 韩国高清视频一区二区三区| 久久久欧美国产精品| 免费高清在线观看视频在线观看| 制服人妻中文乱码| 老熟女久久久| 在线天堂最新版资源| 亚洲综合色网址| 国产老妇伦熟女老妇高清| 中文字幕久久专区| 婷婷色综合大香蕉| 国产在视频线精品| 一区二区日韩欧美中文字幕 | 亚洲精华国产精华液的使用体验| 在线观看三级黄色| 免费黄频网站在线观看国产| 最近最新中文字幕免费大全7| 亚洲av综合色区一区| 全区人妻精品视频| 久久久国产一区二区| 高清欧美精品videossex| 国产亚洲av片在线观看秒播厂| 永久网站在线| 日日撸夜夜添| 国产成人aa在线观看| 成年人免费黄色播放视频| 国产一区二区在线观看av| 曰老女人黄片| 男人添女人高潮全过程视频| 亚洲精品aⅴ在线观看| 男女国产视频网站| 两个人免费观看高清视频| 在线观看三级黄色| 成人综合一区亚洲| 国产国拍精品亚洲av在线观看| 精品人妻一区二区三区麻豆| 亚洲怡红院男人天堂| 国产精品国产三级国产专区5o| 亚洲中文av在线| 日日撸夜夜添| 2021少妇久久久久久久久久久| av视频免费观看在线观看| 日韩av在线免费看完整版不卡| 国产黄片视频在线免费观看| 女的被弄到高潮叫床怎么办| 精品少妇黑人巨大在线播放| 精品国产一区二区久久| 亚洲精品国产av成人精品| 久久久精品免费免费高清| 黄片无遮挡物在线观看| 精品久久久精品久久久| 18+在线观看网站| 国产一区二区在线观看av| 三上悠亚av全集在线观看| 亚洲精品日韩av片在线观看| 亚洲国产成人一精品久久久| 国产精品一国产av| 亚洲精品中文字幕在线视频| 亚洲色图 男人天堂 中文字幕 | 热re99久久精品国产66热6| 欧美亚洲 丝袜 人妻 在线| 欧美精品亚洲一区二区| freevideosex欧美| 热re99久久精品国产66热6| 高清黄色对白视频在线免费看| 色视频在线一区二区三区| 亚洲精品456在线播放app| 欧美三级亚洲精品| 久久综合国产亚洲精品| 亚洲精品456在线播放app| 国产男女超爽视频在线观看| 精品国产一区二区久久| 精品国产国语对白av| 乱码一卡2卡4卡精品| 亚洲精品乱码久久久久久按摩| 99九九在线精品视频| 伦理电影大哥的女人| 99久久人妻综合| 久久精品国产自在天天线| 亚洲av综合色区一区| 日本与韩国留学比较| 麻豆乱淫一区二区| freevideosex欧美| 亚洲精品中文字幕在线视频| 精品国产乱码久久久久久小说| 久久久久久久大尺度免费视频| av免费观看日本| 成人二区视频| 人人妻人人澡人人看| 亚洲国产成人一精品久久久| 亚洲国产av新网站| 亚洲国产毛片av蜜桃av| 十分钟在线观看高清视频www| 国产日韩欧美亚洲二区| 亚洲国产欧美在线一区| 一区二区日韩欧美中文字幕 | 久久99热这里只频精品6学生| 国产一区二区三区综合在线观看 | 内地一区二区视频在线| 久久久久久久久大av| 国产色爽女视频免费观看| 十分钟在线观看高清视频www| 熟妇人妻不卡中文字幕| 菩萨蛮人人尽说江南好唐韦庄| 亚洲国产精品国产精品| 日本wwww免费看| 国产乱来视频区| 在线观看一区二区三区激情| 天天影视国产精品| 欧美精品人与动牲交sv欧美| 99九九线精品视频在线观看视频| a 毛片基地| 成人午夜精彩视频在线观看| 黑人欧美特级aaaaaa片| 91精品三级在线观看| 久久97久久精品| 国产欧美亚洲国产| 18禁在线无遮挡免费观看视频| 国产av精品麻豆| 大香蕉97超碰在线| 老司机影院毛片| 亚洲情色 制服丝袜| 婷婷色av中文字幕| 97精品久久久久久久久久精品| 18+在线观看网站| 日韩制服骚丝袜av| 国产成人aa在线观看| 大陆偷拍与自拍| 欧美性感艳星| 久久国产精品男人的天堂亚洲 | 中文字幕人妻丝袜制服| av在线播放精品| 亚洲久久久国产精品| 久久久久久久精品精品| 亚洲人成77777在线视频| 日韩欧美精品免费久久| 寂寞人妻少妇视频99o| a级毛片黄视频| 成人毛片a级毛片在线播放| 精品人妻熟女av久视频| 大香蕉久久成人网| 99热这里只有精品一区| 日日撸夜夜添| 日韩欧美一区视频在线观看| 一级二级三级毛片免费看| 成年人免费黄色播放视频| 美女内射精品一级片tv| 两个人的视频大全免费| 少妇熟女欧美另类| 永久免费av网站大全| 青春草亚洲视频在线观看| 国产精品免费大片| 国产老妇伦熟女老妇高清| 免费看光身美女| xxx大片免费视频| 啦啦啦啦在线视频资源| 一级毛片黄色毛片免费观看视频| 欧美国产精品一级二级三级| 久热久热在线精品观看| 99视频精品全部免费 在线| 中文字幕最新亚洲高清| 久久精品熟女亚洲av麻豆精品| 中文欧美无线码| 夜夜骑夜夜射夜夜干| 国产伦理片在线播放av一区| 少妇人妻久久综合中文| 久久精品久久久久久久性| 波野结衣二区三区在线| 成人18禁高潮啪啪吃奶动态图 | 精品久久久噜噜| 亚洲av男天堂| 自线自在国产av| 3wmmmm亚洲av在线观看| 美女视频免费永久观看网站| 97超视频在线观看视频| 一级黄片播放器| 免费观看在线日韩| 99九九在线精品视频| 亚洲国产日韩一区二区| 亚洲欧洲精品一区二区精品久久久 | 午夜老司机福利剧场| 高清黄色对白视频在线免费看| 99热6这里只有精品| 精品久久久久久久久亚洲| 精品亚洲乱码少妇综合久久| 精品久久国产蜜桃| 看免费成人av毛片| 亚洲av免费高清在线观看| 男女啪啪激烈高潮av片| 中文字幕免费在线视频6| 我的女老师完整版在线观看| 18+在线观看网站| 一二三四中文在线观看免费高清| 亚洲av日韩在线播放| 99热这里只有是精品在线观看| 国产国语露脸激情在线看| 久久久a久久爽久久v久久| 国产精品三级大全| 精品少妇黑人巨大在线播放| 夜夜骑夜夜射夜夜干| 亚洲一区二区三区欧美精品| 内地一区二区视频在线| 精品久久久久久久久av| 国产欧美亚洲国产| 女性被躁到高潮视频| 街头女战士在线观看网站| 97在线人人人人妻| 国产视频首页在线观看| .国产精品久久| 水蜜桃什么品种好| 日本午夜av视频| 乱码一卡2卡4卡精品| a级毛片免费高清观看在线播放| 久久久久国产网址| 国产老妇伦熟女老妇高清| 91aial.com中文字幕在线观看| 亚洲图色成人| 一级毛片电影观看| 成人亚洲精品一区在线观看| 国产片特级美女逼逼视频| 男人操女人黄网站| 天天影视国产精品| 校园人妻丝袜中文字幕| 免费观看在线日韩| 超碰97精品在线观看| 十分钟在线观看高清视频www| 国产成人午夜福利电影在线观看| 丝袜美足系列| 午夜日本视频在线| 日日啪夜夜爽| 搡老乐熟女国产| 国产在线免费精品| 成年人午夜在线观看视频| 亚洲天堂av无毛| 亚洲第一av免费看| 亚洲精品久久午夜乱码| 三级国产精品欧美在线观看| 91精品国产九色| 中文精品一卡2卡3卡4更新| 国产毛片在线视频| 午夜91福利影院| 亚洲一区二区三区欧美精品| 精品亚洲乱码少妇综合久久| 一级毛片 在线播放| 一级爰片在线观看| 国国产精品蜜臀av免费| 伦理电影免费视频| 日韩精品免费视频一区二区三区 | 大陆偷拍与自拍| 色婷婷久久久亚洲欧美| 欧美日韩国产mv在线观看视频| 亚洲人成77777在线视频| 女性被躁到高潮视频| 亚洲精品视频女| 婷婷成人精品国产| 日韩不卡一区二区三区视频在线| 十八禁高潮呻吟视频| 国产精品人妻久久久影院| 亚洲第一av免费看| 在线观看三级黄色| 两个人免费观看高清视频| 亚洲欧美成人精品一区二区| 一边亲一边摸免费视频| 久久狼人影院| 丰满乱子伦码专区| 在线播放无遮挡| 国国产精品蜜臀av免费| 考比视频在线观看| 欧美少妇被猛烈插入视频| 日韩一区二区视频免费看| 国产亚洲欧美精品永久| 又黄又爽又刺激的免费视频.| 少妇人妻 视频| 飞空精品影院首页| 九九久久精品国产亚洲av麻豆| 美女内射精品一级片tv| 久久久久久久久久久丰满| 色婷婷av一区二区三区视频| 久久久久久久亚洲中文字幕| 丰满乱子伦码专区| 亚洲精品美女久久av网站| av电影中文网址| 边亲边吃奶的免费视频| 在线免费观看不下载黄p国产| 国产黄频视频在线观看| 国产爽快片一区二区三区| 日韩大片免费观看网站| 午夜视频国产福利| 天堂8中文在线网| 成人免费观看视频高清| 毛片一级片免费看久久久久| 我的老师免费观看完整版| 国产一级毛片在线| 如日韩欧美国产精品一区二区三区 | 日本av手机在线免费观看| 午夜福利,免费看| 亚洲国产最新在线播放| 在线天堂最新版资源| 中文乱码字字幕精品一区二区三区| 日本-黄色视频高清免费观看| 卡戴珊不雅视频在线播放| 亚洲精品aⅴ在线观看| 97精品久久久久久久久久精品| 亚洲色图 男人天堂 中文字幕 | 日韩 亚洲 欧美在线| 国产欧美日韩综合在线一区二区| 久久久久久久久久久丰满| 欧美国产精品一级二级三级| 亚洲精品美女久久av网站| 国模一区二区三区四区视频| 久久99一区二区三区| 亚洲欧美中文字幕日韩二区| 欧美精品一区二区免费开放| 91成人精品电影| 多毛熟女@视频| 国产探花极品一区二区| av国产久精品久网站免费入址| 国产成人精品久久久久久| 日本黄色片子视频| 亚洲国产av影院在线观看| 日本欧美视频一区| 男女边吃奶边做爰视频| 亚洲精品一二三| 免费观看无遮挡的男女| 国产精品一区二区在线观看99| 免费观看的影片在线观看| 22中文网久久字幕| 国产视频内射| 最近最新中文字幕免费大全7| 国产精品国产av在线观看|