• <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)士”
    99国产精品免费福利视频| 欧美黄色片欧美黄色片| 一边摸一边抽搐一进一出视频| 久久热在线av| tocl精华| 人妻一区二区av| 精品一区二区三区视频在线观看免费 | 亚洲国产欧美一区二区综合| tube8黄色片| 在线观看66精品国产| 一级毛片精品| 久久久久久亚洲精品国产蜜桃av| 大型av网站在线播放| 午夜福利免费观看在线| 美女午夜性视频免费| 亚洲中文字幕日韩| 丰满迷人的少妇在线观看| 91国产中文字幕| 人人妻人人澡人人爽人人夜夜| 午夜福利乱码中文字幕| 国产亚洲一区二区精品| 777米奇影视久久| 99riav亚洲国产免费| 国产精品一区二区在线观看99| 动漫黄色视频在线观看| 黄色毛片三级朝国网站| 首页视频小说图片口味搜索| 性色av乱码一区二区三区2| av天堂在线播放| 国产免费男女视频| 国产单亲对白刺激| a级毛片在线看网站| 午夜免费观看网址| 丝袜在线中文字幕| 老司机亚洲免费影院| 久久九九热精品免费| 国产精品.久久久| 国产一区有黄有色的免费视频| 在线观看66精品国产| 看片在线看免费视频| av一本久久久久| 在线观看免费日韩欧美大片| 久久国产精品大桥未久av| 日韩 欧美 亚洲 中文字幕| 欧美日韩亚洲高清精品| 老司机深夜福利视频在线观看| 国产精品免费大片| 欧美乱妇无乱码| 叶爱在线成人免费视频播放| 国产欧美日韩一区二区三| 亚洲 国产 在线| 激情视频va一区二区三区| 欧美乱色亚洲激情| 国产精品电影一区二区三区 | 少妇的丰满在线观看| 国产成人影院久久av| 亚洲九九香蕉| 成年人免费黄色播放视频| 亚洲精品久久成人aⅴ小说| 国产精品影院久久| 国产无遮挡羞羞视频在线观看| 看黄色毛片网站| 欧美精品av麻豆av| 可以免费在线观看a视频的电影网站| 婷婷精品国产亚洲av在线 | av视频免费观看在线观看| 日韩 欧美 亚洲 中文字幕| 一a级毛片在线观看| 91麻豆精品激情在线观看国产 | 老熟女久久久| 亚洲美女黄片视频| 欧美在线黄色| 国产蜜桃级精品一区二区三区 | 法律面前人人平等表现在哪些方面| 日韩欧美在线二视频 | 在线观看免费视频网站a站| 80岁老熟妇乱子伦牲交| 天堂动漫精品| 国产精品九九99| 99国产综合亚洲精品| 午夜亚洲福利在线播放| 人成视频在线观看免费观看| 99riav亚洲国产免费| 国产欧美日韩综合在线一区二区| 午夜精品久久久久久毛片777| 久久香蕉激情| 欧美丝袜亚洲另类 | 亚洲中文av在线| 精品电影一区二区在线| 亚洲成人免费电影在线观看| 9191精品国产免费久久| 黄片播放在线免费| 欧美日韩亚洲国产一区二区在线观看 | 国产亚洲精品久久久久久毛片 | 亚洲全国av大片| 757午夜福利合集在线观看| 99国产精品99久久久久| av一本久久久久| 国产单亲对白刺激| 国产国语露脸激情在线看| 丝袜人妻中文字幕| 亚洲aⅴ乱码一区二区在线播放 | 美女高潮喷水抽搐中文字幕| 亚洲少妇的诱惑av| 国产精品综合久久久久久久免费 | 香蕉丝袜av| 村上凉子中文字幕在线| 欧美不卡视频在线免费观看 | 亚洲第一欧美日韩一区二区三区| 老鸭窝网址在线观看| 一个人免费在线观看的高清视频| 久久久久久久午夜电影 | 欧美中文综合在线视频| 亚洲熟女精品中文字幕| 亚洲色图 男人天堂 中文字幕| 午夜福利在线观看吧| 国产精品av久久久久免费| 99国产精品免费福利视频| 不卡一级毛片| 精品国产一区二区三区久久久樱花| 国产蜜桃级精品一区二区三区 | 久久久久久免费高清国产稀缺| 精品国产一区二区久久| 免费高清在线观看日韩| 中文字幕人妻丝袜一区二区| 国产又爽黄色视频| 国产在线一区二区三区精| 丰满饥渴人妻一区二区三| 久久国产精品影院| 久久国产亚洲av麻豆专区| √禁漫天堂资源中文www| x7x7x7水蜜桃| 黄频高清免费视频| 国产一区二区三区综合在线观看| 日本欧美视频一区| 久久精品国产a三级三级三级| 身体一侧抽搐| 欧美激情久久久久久爽电影 | 中文字幕人妻丝袜制服| 久久久久精品人妻al黑| 99热只有精品国产| 80岁老熟妇乱子伦牲交| 亚洲欧美激情在线| 9191精品国产免费久久| 黄色毛片三级朝国网站| 国产片内射在线| 人人妻人人澡人人看| 中国美女看黄片| 亚洲aⅴ乱码一区二区在线播放 | 国产真人三级小视频在线观看| 在线观看www视频免费| 热re99久久国产66热| 大码成人一级视频| 午夜精品在线福利| 久久久久久久精品吃奶| 亚洲九九香蕉| 国产精品九九99| 99精品欧美一区二区三区四区| 日韩欧美一区二区三区在线观看 | 久热爱精品视频在线9| 日韩免费高清中文字幕av| 一级片'在线观看视频| 欧美日韩一级在线毛片| 亚洲熟女精品中文字幕| 深夜精品福利| 国产xxxxx性猛交| 在线观看66精品国产| 国产精品免费视频内射| 国产精品免费视频内射| 精品国产乱子伦一区二区三区| 国产高清国产精品国产三级| 国产成人啪精品午夜网站| 国产亚洲欧美98| 国产国语露脸激情在线看| 69精品国产乱码久久久| 欧美亚洲 丝袜 人妻 在线| 视频在线观看一区二区三区| 亚洲一码二码三码区别大吗| 99国产综合亚洲精品| 国产精品一区二区在线观看99| 国产三级黄色录像| 日韩三级视频一区二区三区| 久久香蕉国产精品| 99re6热这里在线精品视频| 99国产精品一区二区三区| 成人18禁在线播放| 国产成人啪精品午夜网站| 少妇猛男粗大的猛烈进出视频| 国产人伦9x9x在线观看| 国产一区二区激情短视频| 免费高清在线观看日韩| 男人的好看免费观看在线视频 | 身体一侧抽搐| 国产一区有黄有色的免费视频| 国产主播在线观看一区二区| 天天影视国产精品| 亚洲伊人色综图| 久久久精品区二区三区| 九色亚洲精品在线播放| 亚洲va日本ⅴa欧美va伊人久久| xxx96com| 中文字幕最新亚洲高清| 欧美黄色淫秽网站| 在线观看日韩欧美| 99riav亚洲国产免费| 日本a在线网址| 日韩欧美免费精品| 黄色视频,在线免费观看| 丁香欧美五月| 国产亚洲精品久久久久久毛片 | 亚洲伊人色综图| 亚洲欧美激情在线| 国产精品国产高清国产av | 久久久精品免费免费高清| 精品一区二区三卡| 成人国产一区最新在线观看| 精品一区二区三区视频在线观看免费 | 一进一出好大好爽视频| 色94色欧美一区二区| 99riav亚洲国产免费| 人妻久久中文字幕网| 夜夜躁狠狠躁天天躁| 757午夜福利合集在线观看| 好男人电影高清在线观看| 宅男免费午夜| 国产成人精品无人区| 色婷婷av一区二区三区视频| 高潮久久久久久久久久久不卡| 亚洲色图av天堂| av线在线观看网站| 免费少妇av软件| 亚洲av日韩精品久久久久久密| 国产黄色免费在线视频| 精品欧美一区二区三区在线| 色播在线永久视频| 热re99久久精品国产66热6| av欧美777| 丝袜美腿诱惑在线| 一区在线观看完整版| 欧美中文综合在线视频| av电影中文网址| 亚洲专区中文字幕在线| a级片在线免费高清观看视频| 一夜夜www| 老司机午夜十八禁免费视频| 丰满人妻熟妇乱又伦精品不卡| 天天躁日日躁夜夜躁夜夜| www.自偷自拍.com| 啦啦啦 在线观看视频| 男人的好看免费观看在线视频 | 亚洲第一av免费看| 久久精品人人爽人人爽视色| 国产色视频综合| 日韩欧美在线二视频 | 亚洲欧美一区二区三区黑人| x7x7x7水蜜桃| 亚洲精品久久午夜乱码| 亚洲三区欧美一区| 国产淫语在线视频| 叶爱在线成人免费视频播放| 亚洲一卡2卡3卡4卡5卡精品中文| 国产高清国产精品国产三级| 色婷婷久久久亚洲欧美| 黑人巨大精品欧美一区二区蜜桃| 久久精品国产清高在天天线| 性色av乱码一区二区三区2| 久久久久久亚洲精品国产蜜桃av| 亚洲国产欧美网| 亚洲专区中文字幕在线| 日韩欧美三级三区| 亚洲欧美一区二区三区黑人| 女警被强在线播放| 怎么达到女性高潮| 欧美黄色片欧美黄色片| 老司机亚洲免费影院| 免费日韩欧美在线观看| 女人爽到高潮嗷嗷叫在线视频| 在线国产一区二区在线| 午夜日韩欧美国产| 黄色毛片三级朝国网站| 一区二区三区国产精品乱码| 国产成人免费无遮挡视频| 91老司机精品| 夜夜躁狠狠躁天天躁| 精品高清国产在线一区| 丝袜美足系列| 一本综合久久免费| 欧美日韩瑟瑟在线播放| 亚洲欧美精品综合一区二区三区| 欧美成人午夜精品| 免费少妇av软件| 国产精品av久久久久免费| 青草久久国产| 怎么达到女性高潮| 999久久久精品免费观看国产| 在线免费观看的www视频| 十八禁人妻一区二区| 中文字幕最新亚洲高清| 精品国产亚洲在线| 十八禁网站免费在线| 69精品国产乱码久久久| 成人精品一区二区免费| 午夜免费成人在线视频| 极品教师在线免费播放| 欧美黄色淫秽网站| 欧美激情久久久久久爽电影 | 亚洲片人在线观看| 久久人妻熟女aⅴ| 久久久久精品国产欧美久久久| 亚洲精品成人av观看孕妇| 亚洲美女黄片视频| 69av精品久久久久久| 亚洲第一欧美日韩一区二区三区| www日本在线高清视频| 一级a爱片免费观看的视频| 狠狠狠狠99中文字幕| av网站在线播放免费| 又黄又爽又免费观看的视频| 亚洲片人在线观看| 亚洲五月婷婷丁香| 一级毛片女人18水好多| 三上悠亚av全集在线观看| 老汉色av国产亚洲站长工具| 国产国语露脸激情在线看| 亚洲专区字幕在线| 精品一品国产午夜福利视频| 男人的好看免费观看在线视频 | 精品人妻在线不人妻| 69av精品久久久久久| 久久国产精品人妻蜜桃| 少妇粗大呻吟视频| 久久久久精品人妻al黑| 美女 人体艺术 gogo| 精品久久久精品久久久| 久久国产精品男人的天堂亚洲| 久久天躁狠狠躁夜夜2o2o| 国产片内射在线| 每晚都被弄得嗷嗷叫到高潮| 人人妻人人澡人人爽人人夜夜| 免费不卡黄色视频| 欧美日韩瑟瑟在线播放| 一a级毛片在线观看| 久久性视频一级片| 亚洲精品av麻豆狂野| 午夜免费成人在线视频| 日本a在线网址| 成人影院久久| 欧美日韩亚洲高清精品| 亚洲九九香蕉| 法律面前人人平等表现在哪些方面| 精品国产亚洲在线| 成人亚洲精品一区在线观看| 99久久综合精品五月天人人| av免费在线观看网站| 亚洲第一青青草原| 操出白浆在线播放| 亚洲专区字幕在线| 日韩精品免费视频一区二区三区| 亚洲av欧美aⅴ国产| 亚洲欧洲精品一区二区精品久久久| 欧美亚洲日本最大视频资源| 老司机影院毛片| 精品一区二区三区四区五区乱码| 99国产精品一区二区三区| 精品国产美女av久久久久小说| 午夜精品在线福利| 国产成人免费无遮挡视频| 亚洲国产精品sss在线观看 | 免费黄频网站在线观看国产| 性色av乱码一区二区三区2| 夫妻午夜视频| 精品一区二区三区av网在线观看| 精品国产亚洲在线| 久久久国产精品麻豆| 校园春色视频在线观看| 天堂俺去俺来也www色官网| 日韩三级视频一区二区三区| 十分钟在线观看高清视频www| 国产精品一区二区精品视频观看| 91精品三级在线观看| x7x7x7水蜜桃| 午夜激情av网站| 深夜精品福利| 黄色片一级片一级黄色片| 中文字幕精品免费在线观看视频| 搡老岳熟女国产| 母亲3免费完整高清在线观看| 看免费av毛片| av天堂在线播放| 精品一区二区三区av网在线观看| 国产高清视频在线播放一区| 女人久久www免费人成看片| 国产男女内射视频| ponron亚洲| 人成视频在线观看免费观看| xxx96com| 亚洲精品国产区一区二| 黄色视频,在线免费观看| 日韩欧美在线二视频 | 亚洲精品中文字幕一二三四区| 免费久久久久久久精品成人欧美视频| 亚洲中文字幕日韩| 亚洲人成77777在线视频| 欧美精品av麻豆av| 欧美精品啪啪一区二区三区| 黄色女人牲交| 国产国语露脸激情在线看| www日本在线高清视频| 成年动漫av网址| 成人免费观看视频高清| 久久草成人影院| 丝袜美足系列| 日韩欧美国产一区二区入口| 国产精品久久久久久人妻精品电影| 夜夜躁狠狠躁天天躁| 国产日韩欧美亚洲二区| 国产精品久久久久久人妻精品电影| 看免费av毛片| 精品一区二区三卡| 午夜精品久久久久久毛片777| 成年版毛片免费区| 啦啦啦 在线观看视频| 老司机福利观看| 12—13女人毛片做爰片一| 亚洲aⅴ乱码一区二区在线播放 | 午夜福利一区二区在线看| 欧美乱妇无乱码| 女人爽到高潮嗷嗷叫在线视频| 国产欧美日韩一区二区三| 久久久精品国产亚洲av高清涩受| 成年版毛片免费区| 老汉色av国产亚洲站长工具| 少妇的丰满在线观看| 久热爱精品视频在线9| 亚洲av电影在线进入| 久久午夜综合久久蜜桃| 十八禁高潮呻吟视频| 十八禁人妻一区二区| 女人爽到高潮嗷嗷叫在线视频| 91九色精品人成在线观看| 中文亚洲av片在线观看爽 | 成人国语在线视频| 亚洲成人免费电影在线观看| 久久国产精品大桥未久av| 9色porny在线观看| 无遮挡黄片免费观看| 久久午夜亚洲精品久久| 好看av亚洲va欧美ⅴa在| 国产一卡二卡三卡精品| 视频区图区小说| 久久国产精品影院| 国产av精品麻豆| 久久人人爽av亚洲精品天堂| 日韩成人在线观看一区二区三区| 在线免费观看的www视频| 欧美日韩亚洲综合一区二区三区_| 色尼玛亚洲综合影院| 国产亚洲精品第一综合不卡| 后天国语完整版免费观看| 精品福利观看| 午夜福利在线免费观看网站| 国产日韩一区二区三区精品不卡| 一级a爱片免费观看的视频| 日韩有码中文字幕| 亚洲精品成人av观看孕妇| 99国产精品99久久久久| 丁香欧美五月| 无人区码免费观看不卡| 一区二区三区激情视频| 国产男女超爽视频在线观看| 脱女人内裤的视频| 少妇被粗大的猛进出69影院| av片东京热男人的天堂| 国产无遮挡羞羞视频在线观看| 久久久久久久久久久久大奶| 国产亚洲精品一区二区www | 黑人猛操日本美女一级片| 国产日韩欧美亚洲二区| 桃红色精品国产亚洲av| 女性被躁到高潮视频| 国产免费现黄频在线看| 国产在线精品亚洲第一网站| av在线播放免费不卡| 精品无人区乱码1区二区| 超碰97精品在线观看| 久久精品91无色码中文字幕| 无限看片的www在线观看| videosex国产| 黄色怎么调成土黄色| 热99久久久久精品小说推荐| 久久久久国产一级毛片高清牌| 女人被躁到高潮嗷嗷叫费观| 久久天躁狠狠躁夜夜2o2o| 不卡av一区二区三区| 大陆偷拍与自拍| 人妻丰满熟妇av一区二区三区 | 99国产综合亚洲精品| 人妻丰满熟妇av一区二区三区 | 午夜福利在线观看吧| 亚洲一区二区三区不卡视频| 久久久水蜜桃国产精品网| 欧洲精品卡2卡3卡4卡5卡区| 无遮挡黄片免费观看| 精品无人区乱码1区二区| 成人av一区二区三区在线看| 这个男人来自地球电影免费观看| 欧美日韩国产mv在线观看视频| 99在线人妻在线中文字幕 | 国产精品久久久人人做人人爽| 成熟少妇高潮喷水视频| 午夜免费成人在线视频| 中文字幕人妻熟女乱码| 身体一侧抽搐| 大香蕉久久成人网| 亚洲欧美日韩另类电影网站| 高清黄色对白视频在线免费看| 午夜免费成人在线视频| 国产色视频综合| 一级黄色大片毛片| 交换朋友夫妻互换小说| www日本在线高清视频| 免费观看人在逋| 热re99久久国产66热| 亚洲精品粉嫩美女一区| 国产精品一区二区免费欧美| 宅男免费午夜| 高清在线国产一区| 91成人精品电影| 亚洲欧美激情综合另类| 免费在线观看影片大全网站| 欧美在线一区亚洲| 精品一区二区三区av网在线观看| 黄色 视频免费看| 欧美黄色片欧美黄色片| 69精品国产乱码久久久| 国产成人av激情在线播放| 午夜精品国产一区二区电影| 亚洲成av片中文字幕在线观看| 午夜精品久久久久久毛片777| 高清在线国产一区| videos熟女内射| 成人黄色视频免费在线看| 一级黄色大片毛片| 久久久久久久精品吃奶| 国产精品国产高清国产av | 窝窝影院91人妻| 啦啦啦 在线观看视频| 高潮久久久久久久久久久不卡| 免费在线观看完整版高清| 最新美女视频免费是黄的| 国产av一区二区精品久久| 亚洲精品国产一区二区精华液| 亚洲少妇的诱惑av| 免费在线观看亚洲国产| 国产精品影院久久| 国产成人一区二区三区免费视频网站| 18禁观看日本| 久久久精品国产亚洲av高清涩受| 亚洲专区中文字幕在线| 亚洲五月婷婷丁香| 在线观看www视频免费| 91国产中文字幕| 妹子高潮喷水视频| 别揉我奶头~嗯~啊~动态视频| 一进一出抽搐gif免费好疼 | 亚洲色图 男人天堂 中文字幕| 动漫黄色视频在线观看| 91国产中文字幕| 免费在线观看影片大全网站| 中文字幕人妻熟女乱码| 精品一区二区三区视频在线观看免费 | 高潮久久久久久久久久久不卡| www.999成人在线观看| 午夜免费鲁丝| 亚洲成av片中文字幕在线观看| 国产在线精品亚洲第一网站| 欧美黄色片欧美黄色片| 淫妇啪啪啪对白视频| 黄色视频不卡| 啦啦啦在线免费观看视频4| 视频区欧美日本亚洲| 人妻丰满熟妇av一区二区三区 | 久久天躁狠狠躁夜夜2o2o| 在线观看66精品国产| 久久久久国内视频| 女同久久另类99精品国产91| 欧美日韩成人在线一区二区| 不卡一级毛片| 69精品国产乱码久久久| 韩国精品一区二区三区| 女警被强在线播放| 十八禁高潮呻吟视频| 中文字幕色久视频| 亚洲精品久久成人aⅴ小说| www.精华液| av在线播放免费不卡| 国产片内射在线| 国产精品二区激情视频| xxxhd国产人妻xxx| 午夜福利免费观看在线| 午夜91福利影院| 欧美日韩亚洲综合一区二区三区_| 又黄又粗又硬又大视频| 高清黄色对白视频在线免费看| 欧美日韩av久久| 中亚洲国语对白在线视频| 在线观看www视频免费| 91国产中文字幕| 久久久久久久国产电影| 国产麻豆69| 免费黄频网站在线观看国产| 精品久久蜜臀av无| av欧美777| 欧美中文综合在线视频| 国产精品国产av在线观看|