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

    Effect of Different Types of Structural Configuration on Air Distribution in a Compact Purification Device

    2020-02-01 09:05:00YANGXuebinLIANGLe梁樂(lè)LIUHao劉浩LIULin劉琳LIANGZhen梁珍WANGFang王方
    關(guān)鍵詞:劉琳王方劉浩

    YANGXuebin,LIANGLe(梁樂(lè)),LIUHao(劉浩),LIULin(劉琳),LIANGZhen(梁珍),WANGFang(王方)

    1 College of Environmental Science and Engineering, Donghua University, Shanghai 201620, China 2 College of Urban Construction and Safety Engineering, Shanghai Institute of Technology, Shanghai 201418, China 3 School of Energy and Environment, Zhongyuan University of Technology, Zhengzhou 450007, China

    Abstract: In some old industrial plants, in order to meet the increasingly strict requirements of pollutant emission limits, it is necessary to install the compact filtration and/or purification devices in a given narrow machine room. Different types of structural configuration might influence air distribution inside these devices. The unreasonable air distribution might lead each part of filtration or purification media to operating at largely different air flow rates. Based on a computational fluid dynamics (CFD) model, this study explores the influence of different outlet positions and different upper heights on the flow field inside chamber. The porous medium model is employed to simulate the air flow in porous media. The changing structural configurations include three positioning cases of the outlet opening and eight height cases of the upper chamber. The root mean square is defined as the non-uniformity coefficient to evaluate the uniformity of air flow distribution. The results show that the farther distance between inlet and outlet openings will bring more uniform air distribution, and the increasing height of upper chamber totally trends to exhibit more uniform air distribution.

    Key words: structural configuration; air distribution; uniform air distribution; compact purification device; coefficient of velocity inhomogeneity

    Introduction

    In order to deal with the particulate or gaseous pollutants produced in painting workshops, some filters or porous adsorbents are usually used to intercept or adsorb pollutants to meet the national or local emission standards. For the treatment process of gaseous pollutants, once the porous media reaches its saturation, the incoming hot air desorbs the pollutants adsorbed in media. The desorbed pollutants with high concentration is treated by thermal combustion or catalytic oxidation, and the relatively clean gas must meet the concentration emission limits specified in the standards before being discharged to the external atmospheric environment.

    To effectively control the pollutants, some materials[1-2]or technologies[3]have been developed and applied to improve purification efficiency. In the treatment process, many materials possessed some disadvantages. It has been reported that the desorption temperature of zeolite molecular sieve[4]can increase to more than 180 ℃. The higher desorption temperature brings higher desorption efficiency and better operating safety. However, the price of zeolite is about five times higher than that of activated carbon. Higher initial investment cost requires more reliable operation of equipment.

    Many researchers have investigated the uniformity of flow distribution in industrial equipment, involving the inlet structural configuration, various geometric parameters and so on.

    Reviewed from the published literatures, some changes have been conducted in inlet structural configurations to improve flow characteristics.

    The perforated plates emerging from a wide-angle diffuser have been used to improve the air flow in electrostatic precipitator systems. Also, the inlet velocity profiles could influence the air velocity distribution in the electrostatic precipitator. The flow pressure losses in pipes depend on the change of the attack angle[5]. Different configurations of an asymmetric wide-angle diffuser could change the velocity distribution downstream[6]. Bayazitetal.[7]studied the influence of plate geometry on air flow uniform distribution in wind tunnels. Choietal.[8]investigated the flow uniform distribution in a a two- dimensional channel based on non-uniform porosity distribution.

    The box complex method was used to explore the optimal positioning of guide plates in flow manifolds[9]. The shape optimization of a tee guide vane was conducted to improve the air flow in ventilation ducts[10].

    Huangetal.[11]investigated the effects of geometric structures of air deflectors on thermo-flow performances in air-cooled condensers. The flow-guide device has been used to adjust the pressure distribution of duct branches in the exhaust ventilation system[12]. Zhaoetal.[13]explored the effects of air deflectors on air flow mechanism in dry cooling towers.

    A lot of literatures have investigated various geometric parameters to improve flow characteristics. Czetanyetal.[14]investigated the uniform air distribution on the outlets for rectangular ducts at a constant height and variable widths. For the steam ejector for recovery of waste heat, Wuetal.[15]analyzed the optimal structure parameters based on computational fluid dynamics (CFD) model and orthogonal tests. For the grooved channel with curved flow deflectors, Lorenzini-Gutierrezetal.[16]revealed the effect of different geometric parameters (deflector radius, deflector position and channel height) on the thermal and hydraulic behaviors. Nakhchi[17]researched the effects of geometric and flow parameters on pressure drop and thermal performance in sinusoidal wavy channels. Tahanietal.[18]investigated the geometrical parameters of chord, twist distributions and airfoils at different sections to capture the maximum amount of available power. Amanowicz[19]investigated the flow characteristics and total pressure losses to the change of multi-pipe exchanger geometry. For the heat exchanger tubes fitted with perforated conical rings, Nakhchi and Esfahani[20]explored the influence of different geometrical parameters of perforated conical rings on flow structure and heat transfer. For the S-type Pitot measuring the stack gas velocity, Nguyenetal.[21]investigated the effect of various geometries on velocity changes in actual smokestacks. For the subcooled flow boiling in confined microchannels, Hongetal.[22]analyzed the effect of the geometric parameters on bubble behaviors and heat transfer performance. For the throat-tapped flow nozzle, Furuichi and Terao[23]revealed the influence of geometric parameters on the discharge coefficient.

    In this study, the air distribution inside a chamber is investigated by changing the structural configuration. Section 1 overviews the current devices to achieve uniform flow in the published literatures. Section 2 develops a CFD model and verifies the model computational accuracy. Section 3 establishes the simulation model of a purification device. Section 4 investigates the influence of different outlet positions and different upper heights on air distribution.

    1 Numerical Calculation and Verification

    CFD employs numerical methods to directly solve the main flow governing equations in order to discover the phenomena of air distribution. Compared with analytical and experimental fluid dynamics, CFD can provide detailed, visualized and comprehensive information. The ANSYS Fluent 19.2 software is employed as the numerical calculation tool.

    The standardk-εmodel is effectively used to simulate different types of flows such as jets and mixtures. In view of the mixed pollutant gas simulated in this study, the standardk-εmodel is used without considering the heat exchange between gas and device walls.

    Free outflow is selected as the boundary conditions of air inlet and outlet openings. The standard non-slip wall boundary conditions are selected for wall conditions. All surfaces except air inlet and outlet openings are set as wall conditions. Large meshes are used in the numerical model. The parameter settings of the control equation are shown in Table 1.

    Table 1 Parameter settings of control equation in numerical simulation model

    1.1 Verification of numerical calculation

    The numerical calculation should be verified to identify and quantify the errors in solutions, and to assess how accurately the computational results meet with the experimental data[24]. The best way to check whether parameter settings of the control equation, boundary conditions and modeling process are correct or not, is to use the field measured data. In most cases, however, the field measurement of full-scale equipment is time-consuming and labor-intensive. In this study, the validated data are from Ref.[25]. The scale-reduction wind tunnel experiment is with a scaled ratio of 1∶30. A numerical model is developed under the same conditions. As shown in Fig. 1, the difference between simulated and experimental velocities at each measuring point is less than 0.5 m/s, and the average difference is 0.19 m/s. The correlation coefficientR2is 0.982. Some non-ideal factors in the simulating process, as well as the great difference between the finite element model and the actual flow scenario, might lead to a certain degree of deviation between the simulated and measured values.

    Fig. 1 Comparison between simulated and measured value

    1.2 Analysis of computational errors

    The errors between numerical calculation and experimental values are from the approximate error of the physical model, the truncation error of difference equations, the discrete error, the iteration error and the rounding error of the solving region.

    Generally, the discrete error decreases with the decrease of mesh size, but the number of discrete points increases with the decrease of the mesh size, and the rounding error increases with the increasing mesh size. More meshes are not always better, and too large or too small might lead to great calculating errors. The mesh number within a certain range will bring close calculated results to the experimental values.

    This requires that the mesh should be roughly divided and solved according to the parameters of the existing model. The simulated results are compared with the experimental values. Then the mesh is partially encrypted or reduced according to the comparison. If the difference between two last calculation results is small, the computational results are considered acceptable.

    2 CFD Simulation Model

    In section 2, the parameter settings of the control equation, boundary conditions and modeling process in the numerical simulation model, have been checked by the experimental data from Ref.[25]. This section establishes the simulated model of a purification device.

    2.1 Geometric structure

    The geometric structure and spatial dimensions of an adsorbent device is shown in Fig. 2.His the distance from the adsorbent surface to the chamber top (mm), andLis the length of the adsorbent chamber (mm). The inlet opening is located atH=1 500 mm andL=50 mm on the top of the device, and the air outlet is located atH=0 andL=950 mm. The sizes of both inlet and outlet openings are 600 mm× 600 mm. The size of the upper chamber is 1 600 mm× 900 mm × 450 mm, the size of porous media is 1 600 mm × 900 mm × 600 mm, and the size of lower chamber is 1 600 mm × 900 mm × 450 mm.

    1-deflector; 2-inlet opening; 3-upper chamber; 4-purification material; 5-lower chamber; 6-outlet opening.Fig. 2 Geometric structure and spatial dimensions of purification device

    2.2 Simulation model

    Reynolds number is the criterion to judge whether the fluid flow belongs to laminar or turbulent flow. The calculated Reynolds number is 39 534 which is far greater than the critical value 2 320[26], and the air flow is turbulent. The air flow rate through the purification device is 5 000 m3/h, and thus the face velocity at the inlet opening is 3.86 m/s. The air temperature is nearly 20 ℃, and the working condition is much close to the normal atmospheric pressure. The air density is 1.293 kg/m3, and the aerodynamic viscosity coefficient is 1.79×10-5Pa·s.

    Porous medium model (PMM) and porous jump model (PJM) are two theoretical models to deal with filtration, heat transfer and air flow related to porous media. PMM can simulate the internal flow in porous material and can reflect their main airflow characteristics. The calculating results of resistance coefficient are shown in Table 2.

    Table 2 Calculating results of resistance coefficient

    The air flow can be regarded as an incompressible and non-high-speed fluid, and thus the SIMPLE pressure-velocity coupling algorithm is employed in the solver algorithm. The hybrid meshes are used in the mesh generation, and the second order upwind discrete format is adopted for each item.

    2.3 Mesh independence verification

    The finite element model is used in the numerical computation. The mesh quality directly affects the computing accuracy. The lower mesh density will lead to larger errors and even lower computational accuracy. Too higher mesh density will improve the computational accuracy, but will lead to too much calculation.

    The model using structured meshes has the advantages of fast calculation speed, good convergence and relatively higher computational accuracy. For the model with a complex geometric structure, unfortunately, it is difficult to divide all parts into structured meshes. On the contrary, the hybrid meshes combining unstructured meshes with structured meshes are used in some areas. In order to ensure that the simulated results are independent of the number of meshes, it is necessary to test the mesh independence in order to minimize the error between numerical calculation and experimental measurements.

    Given the boundary conditions and solution parameters (the inlet velocity is 3.85 m/s, the number of deflectors is 7, and the deflector width is 50 mm), the model calculations are compared by three mesh cases including the mesh numbers of 1.4 million, 3 million and 6.4 million. The residual of each variable will converge at about 2 000 iterations. The calculated results show that the average velocity of the three mesh models is 4.31 m/s at the highest position planeH=1.3 m. Their maximum difference of the air velocity is 0.09 m/s and the error is 4.28%. The calculated results are basically independent of the number of meshes and thus 1.4 million meshes are employed in the subsequent numerical model.

    3 Case Study

    Three cross sections are atW=0.45 m,L=0.35 m andH=1.06 m in the adsorbent device, respectively. The cross sectionW=0.45 m is intercepted at the plane in the center of the device, the cross sectionL=0.35 m is at the plane in the center of the inlet opening, and the cross sectionH=1.06 m is located at the inlet airflow edge of porous media. The flow characteristic at sectionH=1.06 m represents the air distribution entering the porous media.

    3.1 Different position of outlet openings

    This section studies the influence of the outlet position on the internal flow field of the purification device. Given that the air volume is 5 000 m3/h, the average velocity at the inlet opening is 3.86 m/s, the deflector number is 11, the width of deflectors is 100 mm, and the upper chamber height is 0.45 m.

    Three positioning cases of the outlet opening include: (1) right below the inlet andL=0.35 m; (2) center of bottom of the chamber andL=0.80 m; and (3) farthest below the inlet andL=1.25 m.

    3.1.1Evaluationofairdistribution

    Figure 3 illustrates the velocity distribution ofH-direction on the planeW=0.45 m at three different positions of the outlet opening, which includes right below, central bottom and farthest below. It can be seen from Fig. 3 (a) that the velocity on the surface of porous media changes obviously, and the velocity of the outlet position is the same as that of the inlet position. The air velocities in theL-axis region from 0 to 0.4 m are higher than those in other regions. This means that more polluting gas flows through the left side of the purification medium. After long-time operation, the left region will reach the service life earlier than other regions.

    As depicted in Fig. 3(b), the position of the outlet is in the middle of bottom. The inconsistency of the airflow velocity on the upper surface of porous media is improved but there still exists small parts where the downward velocities are higher than those in other regions.

    Fig. 3 Distribution of H-direction air velocity vector on the plane W=0.45 m at different central positions of outlet opening: (a) right below the inlet and L=0.35 m; (b) central bottom of the chamber and L=0.80 m; (c) farthest below the inlet and L=1.25 m

    As shown in Fig. 3(c), the downward velocity near the surface of porous media is basically about 1 m/s. The air distribution is so even that each part of adsorbent is used equally and no parts are used excessively. In addition, the air velocity in the range of -1.5 m/s to -1.0 m/s is obviously larger than those at the outlet located on the left side and the middle of the bottom. The diffusion distance is much longer, and the air flow is more uniform.

    Figure 4 shows the downward velocity distribution on the planeH=1.06 m at three different positions of the outlet opening, which includes right below, central bottom and farthest below. On the planeH=1.06 m, the air velocity in the range of -1.5 m/s to -1.0 m/s is gradually enlarged if the position of the outlet opening is getting farther from the inlet. In some regions near the chamber wall, the air velocities range from -0.5 m/s to 0. The main reason is that these regions exist some eddy currents which are formed by the deflected airflow when the air rushes into the chamber and collides into porous media.

    Fig. 4 Distribution of H-direction air velocity vector on the plane H=1.06 m at different central positions of outlet opening: (a) right below the inlet; (b) central bottom of the chamber; (c) farthest below the inlet

    Generally speaking, in the upper chamber, the velocity at the surface of porous media tends to distribute more and more evenly if the distance between inlet and outlet openings is getting farther and farther. In the upper region far from the upper inlet, there are no obvious changes in the vortices with upward velocity. The closer the area is to the side of chamber, the faster upward speed is. The maximum velocity is more than 1.0 m/s. For the velocity ranging from -2.0 m/s to -1.5 m/s, theL-axis region of the right bottom outlet ranges from 0 to 1.1 m, while that of the farthest bottom is from 0 to 1.2 m.

    In the lower chamber, the downward velocities are greater than 1.0 m/s at the region close to the outlet, and the maximum distance away from the edge of the outlet inLdirection is less than 0.2 m. Although the position of the outlet opening is various, the velocity distribution near the outlet is not changed obviously. This indicates that different outlet positions have little influence on air distribution in the lower chamber.

    3.1.2Evaluationofairdistribution

    The root mean square (RMS), which is more sensitive to the non-uniform response of the airflow field, is defined as the non-uniformity coefficient to evaluate the air distribution on the inlet plane of purification media. The RMS value is calculated by the standard deviation of individual air velocity to the average section velocity at the entrance section of purification media. Statistically, the RMS value represents the dispersion degree of the single airflow velocity from the average velocity. Greater RMS value means larger difference among the velocities and less uniformity. The RMS values less than 0.25 can be deemed as qualified, and those greater than 0.25 are unqualified.

    The non-uniformity coefficient is formulated as

    (1)

    Table 3 shows the average surface velocity, the global velocity interval and the non-uniformity coefficient of the downward velocity on the planeH=1.06 m at various outlet locations. For the outlet openings located at the position (L=0.35, 0.80, 1.25 m), the RMS values ofH-direction velocities are 0.232, 0.217, and 0.180, respectively. The best airflow uniformity is the position where the inlet and outlet openings are arranged on the opposite side of the chamber, and the worst one is the position where the inlet and outlet openings are arranged in the same side.

    Table 3 Comparison of H-direction air velocity on the plane H=1.06 m with outlet openings at different positions of outlet openings

    3.2 Different heights in upper chamber

    Air deflectors can effectively improve the uniformity of air distribution in the purification device, but due to the limitation of chamber size, the flow distance in the chamber seems too short to limit the full diffusion of air flow. The whole air diffusion region can be enlarged by changing the height and the size of the upper chamber.

    Given that the air flow rate through the purification device is 5 000 m3/h, thus the average velocity at inlet opening is 3.86 m/s. The outlet opening is on the bottom of the chamber and farthest below the inlet. And 11 deflectors with a width of 100 mm are set up on the inlet opening. As listed in Table 4, the height of the upper chamber is defined as eight cases from 350 mm to 700 mm with an interval of 50 mm.

    Table 4 Different height cases of upper chamber

    3.2.1Distributionofairvelocity

    The ideal condition for the most uniform airflow distribution is that the number of deflectors or perforated plates should be as large as possible and each deflector is long enough to reach the entrance of filtering and purifying media materials. That is, there is enough chamber space between deflectors and media materials to fully diffuse the air, and thus the pressure and velocity distribution on the inlet surface of media materials can be as uniform as possible. The chamber space is enough to avoid the influence of the inlet flow.

    Figure 5 indicates the distribution of theH-direction air velocity vector on the planeL=0.35 m at eight types of upper chamber heights. It is assumed that the velocity along the positiveH-axis direction is positive and the velocity along the negativeH-axis direction is negative.

    Totally, higher chamber height leads to longer air flow path through deflectors towards the media materials. In the region (L=0-1.3 m), the air velocity is from -1.0 m/s to 1.6 m/s. This region tends to expand much more with the increasing chamber height. In the region (L=1.3-1.6 m), the air velocity at the inlet plane of filtration material is from -1.0 m/s to 0. In the upper part of the chamber, there exists a vortex region with upward airflow velocity, and the maximum value is even up to 2.0 m/s. This vortex region with upward flow velocity will increase with the increasing chamber height.

    Fig. 5 Distribution of H-direction air velocity vector on the plane L=0.35 m at different upper chamber heights: (a) 350 mm; (b) 400 mm; (c) 450 mm; (d) 500 mm; (e) 550 mm; (f) 600 mm; (g) 650 mm; (h) 700 mm

    Figure 6 illustrates the distribution of theH-direction air velocity on the planeH=1.06 m under different upper chamber heights. The negativeH-direction velocity in the range of -2.5 m/s to -1.3 m/s mainly concentrates in the region (L=0-0.6 m) where this velocity gradually enlarges with the increasing height of the upper chamber. The velocity in the region (L=0.6-1.4 m) is mainly in the range of -1.3 m/s to -0.4 m/s. In the region from 0 to 0.8 m/s near the chamber wall, the positiveH-direction velocity decreases with the increase of the chamber height. Generally speaking, the velocity distribution trends to be more uniform with the increasing height of the upper chamber.

    Fig. 6 Distribution of H-direction air velocity vector on the plane H=1.06 m at different upper chamber heights: (a) 350 mm; (b) 400 mm; (c) 450 mm; (d) 500 mm; (e) 550 mm; (f) 600 mm; (g) 650 mm; (h) 700 mm

    The distribution of theH-direction air velocity vector on the planeL=0.35 m is shown in Fig. 7. In the upper chamber where the pollutant gas enters, the region of negativeH-direction velocity is mainly converged the regionW=0.2-0.7 m. In the vicinity ofH=1.06 m near the entrance surface of media materials, the velocity of the region (W=0.2-0.7 m) is from -1.5 m/s to -1.0 m/s. In the region near the chamber wall, there will be positiveH-direction velocity due to the existence of vortices. Totally, with the increasing heights of the upper chamber, the area of negativeH-direction velocity from -1.5 m/s to -1.0 m/s decreases, but the area of velocity from -1.0 m/s to -0.5 m/s increases. Similarly, in the region near the chamber wall, the vortices region increases and the maximum of the positiveH-direction velocity decreases.

    Fig. 7 Distribution of H-direction air velocity vector on the plane L=0.35 m at different upper chamber heights: (a) 350 mm; (b) 400 mm; (c) 450 mm; (d) 500 mm; (e) 550 mm; (f) 600 mm; (g) 650 mm; (h) 700 mm

    3.2.2Evaluationofairdistribution

    Assuming that the direction of the air velocity is the same as that ofH-axis, the negative velocity means their direction along the negativeH-axis. Table 5 shows the average surface velocity, velocity interval and non-uniformity coefficient at different upper chamber heights. The average surface velocity is basically maintained at -0.94 m/s, because their total airflow rates are identical and given as 5 000 m3/h.

    When the height of the upper box decreases gradually, the global velocity range shows a trend of increasing at first and then decreasing. When the height of the upper chamber increases gradually, the global velocity interval decreases firstly, then increases and finally decreases. At the same time, when the height of the upper chamber decreases, the non-uniformity coefficient decreases slightly. When the chamber height increases to 550 mm, non-uniformity coefficient reaches 0.225, which exhibits the most non-uniformity of the flow field inside the equipment. When the height of the box exceeds 550 mm, the non-uniformity coefficient begins to decrease gradually. When the height reaches 700 mm, the RMS value is 0.132, and 41.33% less than that at 550 mm.

    Table 5 Comparison of H-direction air velocity on the plane H=1.06 m at different heights

    4 Conclusions

    This study investigates the air distribution under different structural configurations, such as changing the position of the outlet opening and the height in the upper chamber. The non-uniformity coefficient is used to evaluate the air distribution.

    (1) Farther distance between inlet and outlet openings will lead to more uniform air distribution. The case where inlet and outlet openings are arranged on the opposite side of chamber exhibits the best airflow uniformity.

    (2) Increasing the height of the upper chamber properly can improve the uniformity of air flow field. With the increasing height of the upper chamber, the non-uniformity coefficient of velocity distribution totally trends to increase firstly and then decrease.

    (3) The results can contribute to the selection of the compact structural configuration of purification devices. Farther distance between inlet and outlet openings, as well as the chamber height of 350 mm, are recommended in the studied cases.

    猜你喜歡
    劉琳王方劉浩
    才子大比拼
    情景劇小品
    High-fidelity resonant tunneling passage in three-waveguide system
    Theoretical analysis and numerical simulation of acoustic waves in gas hydrate-bearing sediments?
    THE EXTENSION OPERATORS ON Bn+1 AND BOUNDED COMPLETE REINHARDT DOMAINS*
    光影空間
    赤水源(2018年2期)2018-04-21 06:42:42
    Developing Learner Autonomy and Second Language Context
    劉浩藝術(shù)作品欣賞
    消除“鈍”感肌就是這樣滑!
    Coco薇(2015年3期)2015-12-24 03:06:17
    講臺(tái)上的“鋼鐵戰(zhàn)士”
    国产精品久久久久成人av| 人体艺术视频欧美日本| 日本av免费视频播放| 国国产精品蜜臀av免费| a 毛片基地| 中文字幕人妻熟人妻熟丝袜美| 国产精品不卡视频一区二区| 91精品一卡2卡3卡4卡| 人妻一区二区av| 亚洲欧美一区二区三区黑人 | 熟女人妻精品中文字幕| 美女中出高潮动态图| 你懂的网址亚洲精品在线观看| av福利片在线| 亚洲婷婷狠狠爱综合网| av黄色大香蕉| 国产有黄有色有爽视频| 亚洲精品色激情综合| 精品卡一卡二卡四卡免费| 亚洲精品中文字幕在线视频 | 777米奇影视久久| 91精品国产国语对白视频| 色视频www国产| 国产精品国产av在线观看| 精品熟女少妇av免费看| 日本猛色少妇xxxxx猛交久久| 中文字幕精品免费在线观看视频 | 免费看av在线观看网站| 久久久久国产网址| 婷婷色综合大香蕉| 亚洲成人手机| 久久韩国三级中文字幕| 欧美成人精品欧美一级黄| 亚洲无线观看免费| 日本午夜av视频| 噜噜噜噜噜久久久久久91| 亚洲欧美一区二区三区国产| 中文资源天堂在线| 婷婷色综合大香蕉| videossex国产| 亚洲精品乱码久久久v下载方式| 国产欧美亚洲国产| 日本wwww免费看| 久久综合国产亚洲精品| 亚洲精品中文字幕在线视频 | 亚洲内射少妇av| 男人和女人高潮做爰伦理| 丁香六月天网| 伊人亚洲综合成人网| 国产男人的电影天堂91| 精品一区二区三区视频在线| 人人妻人人看人人澡| 久久久久久久久久久久大奶| 青春草亚洲视频在线观看| 午夜激情福利司机影院| 熟女av电影| 在线观看免费日韩欧美大片 | 国产69精品久久久久777片| 亚洲精品一二三| 午夜日本视频在线| 久久精品久久精品一区二区三区| av线在线观看网站| 制服丝袜香蕉在线| 亚州av有码| 色婷婷久久久亚洲欧美| 久久久欧美国产精品| 久久久亚洲精品成人影院| 精品久久国产蜜桃| 韩国高清视频一区二区三区| 国产在线免费精品| 亚洲欧洲国产日韩| 日韩欧美精品免费久久| 亚洲丝袜综合中文字幕| 边亲边吃奶的免费视频| 99九九线精品视频在线观看视频| 91精品国产九色| 日本黄色日本黄色录像| 极品人妻少妇av视频| 亚洲婷婷狠狠爱综合网| 中文字幕久久专区| 只有这里有精品99| 免费观看无遮挡的男女| 又黄又爽又刺激的免费视频.| 肉色欧美久久久久久久蜜桃| 亚洲精品国产成人久久av| 久久人人爽人人片av| 国产又色又爽无遮挡免| 精华霜和精华液先用哪个| 一区二区三区四区激情视频| 九九爱精品视频在线观看| 色网站视频免费| 各种免费的搞黄视频| 老司机亚洲免费影院| 中文字幕久久专区| 亚洲美女黄色视频免费看| 视频中文字幕在线观看| 亚洲精品国产av蜜桃| 久久久久久久久久久免费av| 大片电影免费在线观看免费| 午夜激情福利司机影院| 性高湖久久久久久久久免费观看| 大片电影免费在线观看免费| 另类精品久久| 亚洲成人av在线免费| 最近最新中文字幕免费大全7| 国产亚洲精品久久久com| 精品久久久噜噜| 22中文网久久字幕| 日本av手机在线免费观看| 三上悠亚av全集在线观看 | 国产 精品1| 日韩在线高清观看一区二区三区| 人人妻人人澡人人看| 好男人视频免费观看在线| 老女人水多毛片| 亚洲国产精品一区二区三区在线| 国产精品女同一区二区软件| 亚洲精品一区蜜桃| 三级国产精品欧美在线观看| 噜噜噜噜噜久久久久久91| 亚洲精品456在线播放app| 日韩熟女老妇一区二区性免费视频| 亚洲av男天堂| 亚洲人与动物交配视频| 欧美日本中文国产一区发布| 亚洲av福利一区| 51国产日韩欧美| 国产精品99久久久久久久久| 女人精品久久久久毛片| 日韩 亚洲 欧美在线| 久久av网站| 99热国产这里只有精品6| 欧美日韩国产mv在线观看视频| 欧美一级a爱片免费观看看| 日韩人妻高清精品专区| 中文字幕av电影在线播放| 久久精品国产a三级三级三级| 国产有黄有色有爽视频| 9色porny在线观看| 久久精品国产鲁丝片午夜精品| 夜夜看夜夜爽夜夜摸| 午夜激情久久久久久久| 在线观看www视频免费| 国产毛片在线视频| 久久久久久伊人网av| 亚洲婷婷狠狠爱综合网| 久久狼人影院| 26uuu在线亚洲综合色| 人妻人人澡人人爽人人| 中文资源天堂在线| 久久久久久久久大av| 中文字幕精品免费在线观看视频 | 精品人妻一区二区三区麻豆| 男女边摸边吃奶| 免费观看性生交大片5| 国产日韩欧美亚洲二区| 亚洲精品,欧美精品| av卡一久久| 久久狼人影院| 少妇人妻 视频| 一本一本综合久久| 蜜桃在线观看..| 99热国产这里只有精品6| 亚洲av免费高清在线观看| av视频免费观看在线观看| 如日韩欧美国产精品一区二区三区 | 国产精品不卡视频一区二区| 极品人妻少妇av视频| 韩国高清视频一区二区三区| 欧美最新免费一区二区三区| 噜噜噜噜噜久久久久久91| 蜜桃久久精品国产亚洲av| av福利片在线观看| 偷拍熟女少妇极品色| 精品亚洲乱码少妇综合久久| av在线观看视频网站免费| 男男h啪啪无遮挡| 精品视频人人做人人爽| 欧美日韩在线观看h| 99久久综合免费| 亚洲丝袜综合中文字幕| 日本vs欧美在线观看视频 | 欧美日韩一区二区视频在线观看视频在线| 久久久精品94久久精品| videos熟女内射| 亚洲精品乱码久久久久久按摩| 国产高清有码在线观看视频| 亚洲三级黄色毛片| 成人国产麻豆网| 五月天丁香电影| 欧美激情国产日韩精品一区| 免费久久久久久久精品成人欧美视频 | 国产 一区精品| 能在线免费看毛片的网站| 少妇被粗大的猛进出69影院 | 亚洲中文av在线| 2018国产大陆天天弄谢| 久久精品国产亚洲av天美| 日韩欧美 国产精品| 偷拍熟女少妇极品色| 国产乱人偷精品视频| 交换朋友夫妻互换小说| 久久 成人 亚洲| 亚洲三级黄色毛片| 国产一区二区在线观看av| 久久久久国产精品人妻一区二区| 全区人妻精品视频| 人人妻人人添人人爽欧美一区卜| 两个人免费观看高清视频 | 国产精品麻豆人妻色哟哟久久| 国产精品国产三级专区第一集| 99久久综合免费| 久久毛片免费看一区二区三区| 成人国产麻豆网| 亚洲人与动物交配视频| 午夜福利在线观看免费完整高清在| 亚洲情色 制服丝袜| 涩涩av久久男人的天堂| 久久国产乱子免费精品| 国产精品熟女久久久久浪| 高清视频免费观看一区二区| 亚洲欧洲日产国产| 久久久a久久爽久久v久久| 免费高清在线观看视频在线观看| 国产精品熟女久久久久浪| 国产一区二区三区综合在线观看 | 国产极品粉嫩免费观看在线 | 国产精品久久久久久久久免| 久久热精品热| 美女cb高潮喷水在线观看| 久久精品国产亚洲av涩爱| 亚州av有码| 国产精品久久久久久久久免| 亚洲真实伦在线观看| 国产淫语在线视频| 国产真实伦视频高清在线观看| 精品久久久精品久久久| 亚洲国产精品999| 你懂的网址亚洲精品在线观看| 欧美国产精品一级二级三级 | 一边亲一边摸免费视频| 久久久国产欧美日韩av| 欧美日韩在线观看h| 久久久久网色| 日韩熟女老妇一区二区性免费视频| 国产精品欧美亚洲77777| 黑人巨大精品欧美一区二区蜜桃 | 永久网站在线| 好男人视频免费观看在线| 尾随美女入室| 久久久精品94久久精品| 青春草视频在线免费观看| 国产欧美日韩综合在线一区二区 | 日本-黄色视频高清免费观看| 久久久久久久精品精品| 亚洲精品一二三| 成人国产av品久久久| 51国产日韩欧美| 大陆偷拍与自拍| 偷拍熟女少妇极品色| 国产免费又黄又爽又色| 久久久久久久久大av| 亚洲av综合色区一区| 欧美精品人与动牲交sv欧美| 性色avwww在线观看| 高清视频免费观看一区二区| 啦啦啦在线观看免费高清www| 成人二区视频| 成人毛片60女人毛片免费| 少妇猛男粗大的猛烈进出视频| 欧美一级a爱片免费观看看| 又大又黄又爽视频免费| 91精品国产九色| 亚洲精品视频女| 99久久精品国产国产毛片| 伦理电影大哥的女人| 成人漫画全彩无遮挡| 中文字幕久久专区| av又黄又爽大尺度在线免费看| 午夜精品国产一区二区电影| 哪个播放器可以免费观看大片| 人妻少妇偷人精品九色| 精品国产一区二区三区久久久樱花| 黄色怎么调成土黄色| 久热这里只有精品99| 七月丁香在线播放| 97在线人人人人妻| 亚洲国产av新网站| 狠狠精品人妻久久久久久综合| 免费人妻精品一区二区三区视频| 国产黄色视频一区二区在线观看| 男女免费视频国产| 国产精品熟女久久久久浪| 黄色日韩在线| 91精品国产国语对白视频| 麻豆成人av视频| 男女边吃奶边做爰视频| 久久精品熟女亚洲av麻豆精品| 少妇人妻久久综合中文| 亚洲精品456在线播放app| 午夜日本视频在线| 午夜激情福利司机影院| 男女边摸边吃奶| 男人和女人高潮做爰伦理| 老司机影院毛片| 日日爽夜夜爽网站| 日本黄色日本黄色录像| 免费观看性生交大片5| a级片在线免费高清观看视频| 久久人人爽人人爽人人片va| 成人综合一区亚洲| 国产成人精品婷婷| 日本欧美视频一区| 国产一区亚洲一区在线观看| 国产亚洲av片在线观看秒播厂| 精品久久久久久久久亚洲| 国产乱来视频区| 中文天堂在线官网| 国内精品宾馆在线| 亚洲欧美精品自产自拍| 久久国内精品自在自线图片| 国产熟女欧美一区二区| 亚洲va在线va天堂va国产| 亚洲经典国产精华液单| 在线观看免费高清a一片| 色5月婷婷丁香| 一本—道久久a久久精品蜜桃钙片| 在现免费观看毛片| 国产精品人妻久久久久久| 哪个播放器可以免费观看大片| av卡一久久| 乱人伦中国视频| 熟女电影av网| 成人特级av手机在线观看| 人妻系列 视频| 大香蕉久久网| 亚洲内射少妇av| 特大巨黑吊av在线直播| 免费观看av网站的网址| 久久午夜综合久久蜜桃| 国产日韩欧美亚洲二区| 91午夜精品亚洲一区二区三区| 免费观看性生交大片5| tube8黄色片| 国产精品福利在线免费观看| 欧美日韩国产mv在线观看视频| h视频一区二区三区| 日韩精品有码人妻一区| 亚洲国产色片| 国产一区二区三区综合在线观看 | 婷婷色麻豆天堂久久| 亚洲人与动物交配视频| freevideosex欧美| 色视频www国产| 一区二区三区精品91| 少妇人妻精品综合一区二区| 亚洲av中文av极速乱| 欧美xxxx性猛交bbbb| 交换朋友夫妻互换小说| 少妇 在线观看| 天天躁夜夜躁狠狠久久av| 男人爽女人下面视频在线观看| 天堂中文最新版在线下载| 高清欧美精品videossex| 亚洲成人av在线免费| videossex国产| av天堂中文字幕网| 亚洲综合精品二区| 天堂中文最新版在线下载| 高清欧美精品videossex| av.在线天堂| av播播在线观看一区| 99热这里只有是精品在线观看| 汤姆久久久久久久影院中文字幕| 另类精品久久| 3wmmmm亚洲av在线观看| 乱人伦中国视频| 日本vs欧美在线观看视频 | 丰满人妻一区二区三区视频av| 亚洲av在线观看美女高潮| 精品久久久久久电影网| 国产日韩欧美亚洲二区| 国产男人的电影天堂91| 国产在线男女| 人人澡人人妻人| 色5月婷婷丁香| 视频中文字幕在线观看| 少妇高潮的动态图| 精品国产国语对白av| 欧美区成人在线视频| 亚洲欧美清纯卡通| 人妻夜夜爽99麻豆av| 女性被躁到高潮视频| 亚洲精品aⅴ在线观看| 国产美女午夜福利| 只有这里有精品99| 午夜精品国产一区二区电影| 日日撸夜夜添| 成人美女网站在线观看视频| 国产亚洲午夜精品一区二区久久| 日韩在线高清观看一区二区三区| 精品少妇久久久久久888优播| 人妻制服诱惑在线中文字幕| 嫩草影院入口| 永久网站在线| 免费黄色在线免费观看| 国产精品偷伦视频观看了| 男人狂女人下面高潮的视频| 亚洲欧美成人综合另类久久久| 99久久人妻综合| 色吧在线观看| 久久这里有精品视频免费| 男女国产视频网站| 亚洲在久久综合| 久久 成人 亚洲| 黑丝袜美女国产一区| 国产亚洲精品久久久com| 永久免费av网站大全| 日韩 亚洲 欧美在线| 免费高清在线观看视频在线观看| 国产精品蜜桃在线观看| 七月丁香在线播放| 精品亚洲成a人片在线观看| 欧美精品一区二区免费开放| 欧美日韩av久久| 永久网站在线| 久久人人爽人人片av| 欧美日韩在线观看h| 中文字幕免费在线视频6| 草草在线视频免费看| av在线app专区| 自拍欧美九色日韩亚洲蝌蚪91 | 国产av国产精品国产| 伦理电影大哥的女人| 一级毛片我不卡| 三级国产精品欧美在线观看| 国产精品国产三级国产专区5o| av播播在线观看一区| 嘟嘟电影网在线观看| 九九久久精品国产亚洲av麻豆| 中文在线观看免费www的网站| 亚洲美女视频黄频| 欧美97在线视频| 亚洲欧洲日产国产| 亚洲无线观看免费| 国产成人免费无遮挡视频| 九九在线视频观看精品| av女优亚洲男人天堂| 国产精品国产三级专区第一集| 国产精品久久久久久久久免| 日韩中字成人| 狠狠精品人妻久久久久久综合| 99热网站在线观看| 日本欧美视频一区| 久久影院123| 亚洲精品一二三| 免费观看无遮挡的男女| 熟女电影av网| 热99国产精品久久久久久7| www.av在线官网国产| 午夜福利视频精品| 亚洲第一区二区三区不卡| 水蜜桃什么品种好| 午夜av观看不卡| 伊人久久国产一区二区| 精品午夜福利在线看| 欧美97在线视频| 国产亚洲最大av| 伦精品一区二区三区| 99久久综合免费| 精品国产露脸久久av麻豆| 水蜜桃什么品种好| 日本色播在线视频| videos熟女内射| xxx大片免费视频| 美女中出高潮动态图| 亚洲情色 制服丝袜| www.av在线官网国产| 免费观看的影片在线观看| 国产成人精品无人区| 国产 精品1| 成人黄色视频免费在线看| 亚洲av在线观看美女高潮| 2018国产大陆天天弄谢| 中国美白少妇内射xxxbb| 免费观看性生交大片5| 欧美最新免费一区二区三区| 啦啦啦视频在线资源免费观看| 国产毛片在线视频| 午夜91福利影院| 欧美xxⅹ黑人| 国产片特级美女逼逼视频| 亚洲电影在线观看av| 老司机影院成人| 国内精品宾馆在线| 久久久午夜欧美精品| 久久久久久人妻| 99热网站在线观看| 国产在线男女| 久久综合国产亚洲精品| av有码第一页| 男人爽女人下面视频在线观看| 久久久午夜欧美精品| 蜜桃久久精品国产亚洲av| 日韩中文字幕视频在线看片| 国产精品不卡视频一区二区| 成人无遮挡网站| 国产永久视频网站| 中文字幕免费在线视频6| 久久久久久伊人网av| 另类亚洲欧美激情| 午夜福利视频精品| 久久热精品热| 欧美 日韩 精品 国产| 五月玫瑰六月丁香| 色5月婷婷丁香| 久久韩国三级中文字幕| 亚洲欧美精品自产自拍| 在线观看www视频免费| 女人精品久久久久毛片| 国产精品人妻久久久久久| 亚洲婷婷狠狠爱综合网| 天天操日日干夜夜撸| 欧美成人精品欧美一级黄| 久久久久久久久大av| 蜜桃在线观看..| 一二三四中文在线观看免费高清| 成人毛片60女人毛片免费| 韩国高清视频一区二区三区| 丰满人妻一区二区三区视频av| 国产av一区二区精品久久| 毛片一级片免费看久久久久| 三级国产精品欧美在线观看| 欧美日韩亚洲高清精品| 亚洲国产欧美在线一区| 国产成人一区二区在线| 性色avwww在线观看| 婷婷色麻豆天堂久久| 乱人伦中国视频| 夫妻午夜视频| 亚洲av二区三区四区| 乱系列少妇在线播放| 成人国产麻豆网| 成人美女网站在线观看视频| 久久青草综合色| av.在线天堂| 男人狂女人下面高潮的视频| 日韩中文字幕视频在线看片| 国产一区二区三区综合在线观看 | 午夜免费鲁丝| 九九在线视频观看精品| 在线观看www视频免费| 亚洲国产精品专区欧美| 亚洲国产精品999| 王馨瑶露胸无遮挡在线观看| 欧美亚洲 丝袜 人妻 在线| 精品亚洲乱码少妇综合久久| 亚洲av电影在线观看一区二区三区| 国产熟女欧美一区二区| 一级黄片播放器| 日韩不卡一区二区三区视频在线| 国产在线男女| 一区二区三区乱码不卡18| 日产精品乱码卡一卡2卡三| 涩涩av久久男人的天堂| 少妇裸体淫交视频免费看高清| 免费观看av网站的网址| 嘟嘟电影网在线观看| 啦啦啦在线观看免费高清www| 大香蕉久久网| 亚洲精品成人av观看孕妇| 精品国产一区二区久久| 熟女av电影| 99久国产av精品国产电影| 亚洲天堂av无毛| 成人综合一区亚洲| 日日摸夜夜添夜夜爱| 男人舔奶头视频| 国产日韩欧美亚洲二区| 日韩成人av中文字幕在线观看| 制服丝袜香蕉在线| 亚洲人成网站在线播| 欧美国产精品一级二级三级 | 乱系列少妇在线播放| 少妇人妻 视频| 18+在线观看网站| 久久国内精品自在自线图片| 一区在线观看完整版| 久热这里只有精品99| 亚洲婷婷狠狠爱综合网| 黄片无遮挡物在线观看| 久久久久久久亚洲中文字幕| 啦啦啦视频在线资源免费观看| a级毛片在线看网站| 精品亚洲成a人片在线观看| 亚洲情色 制服丝袜| 国产女主播在线喷水免费视频网站| 欧美成人午夜免费资源| 成人影院久久| 亚洲av日韩在线播放| 欧美日韩亚洲高清精品| 国产亚洲最大av| 全区人妻精品视频| 中文字幕亚洲精品专区| 成人影院久久| 免费观看av网站的网址| 99热这里只有精品一区| 噜噜噜噜噜久久久久久91| 日本-黄色视频高清免费观看| 女的被弄到高潮叫床怎么办| 97在线视频观看| 国产精品久久久久久久电影| 亚洲成人一二三区av| 国国产精品蜜臀av免费| 插逼视频在线观看|