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

    Dynamic Behavior and Heat Transfer Characteristics of Non-spherical Ice Crystals in High-Temperature Air Flow

    2023-05-21 02:37:20,,,,,,,2*

    ,,,,,,,2*

    1.School of Energy and Power Engineering,Xi’an Jiaotong University,Xi’an 710049,P.R.China;

    2.Anti/De-Icing Key Laboratory,China Aerodynamics Research and Development Center(CARDC),Mianyang 621000,P.R.China

    Abstract: This study investigates the dynamic behavior and heat transfer of ice crystals in high-temperature air flow after being ingested by an aero-engine.A numerical model under the Lagrange framework is built to simulate ice crystal motion,heat and mass transfer.Different ice crystal melting models and ice crystals with various shapes and sizes are also compared.The results show that when using the “naked” ice particle model,it takes longer to completely melt compared with the water cover model.It takes 0.04 s for a 20 μm particle to melt,but the melting rate of a “naked” ice particle is faster before 0.023 s.Under the same condition,the liquid water content of a spherical ice crystal is high,and that of a 20 μm spherical one is 49.05% at the outlet.The liquid water contents of ellipsoidal and hexagonal plate ice crystals are close,which are roughly 40% at the outlet.The smaller the size,the earlier the crystals start to melt,and the higher the liquid water content of the crystals.The liquid water content is 60.4% of a 20 μm ice crystal and only 15.5% of a 40 μm one at the outlet.

    Key words:ice crystals;aero-engine;Lagrange method;heat and mass transfer;non-spherical

    0 Introduction

    Aircraft icing problems were recognized as early as the beginning of modern flight.The first reported aircraft icing can be traced back to 1920s[1].Researchers first discovered and studied the phenomena of supercooled water droplets icing.Related prediction and protection have also been widely applied in practice[2].However,since 1990s,there has emerged an increasing number of accidents caused by engine thrust loss at the altitude of 6.7 km even though supercooled water droplets rarely exist there.According to the accident investigation and analysis conducted by the relevant authorities,the main cause of accidents is the adhesion and freezing of ice crystals after being ingested by the engine[3].

    In recent years,intensive investigations have been conducted to resolve the problem of ice crystal icing inside engines via theoretical and numerical modelling,and experiment[4-5].For the ice crystal melting and freezing process,Villedieu et al.[6]and Trontin et al.[7]proposed the ice crystal freezing model based on the supercooled water droplet freezing model,which divided the ice crystal melting phase transition into three processes,namely sublimation of supercooled ice crystal,ice crystal melting and evaporation of droplets after melting completed.Lou et al.[8]proposed three models of ice crystal melting,which are respectively the “naked”ice model,the water film cover model and the porous model.They compared the first two models and showed that the water film cover model had a faster melting rate.Norde[9]and Nilamdeen et al.[10]used the Euler method to simulate ice crystal freezing in airfoils and engines,considering the complex processes of ice crystal adhesion,fragmentation and erosion of ice crystals.Veres et al.[11-13]applied an engine icing computational tool to simulate the flow through a turbofan engine and assess the risk of ice accretion.They found that the ice accumulation was affected by the local wet bulb temperature,ice crystal speed and ice crystal liquid water content.Feulner et al.[14-15]used the Monte Carlo particle trajectory simulation method and studied the motion and crushing process of ice crystal particles ingested by the engine,and found that the ice crystal fragmentation mainly occurred on the low-pressure compressor blades.The internal flow field,ice crystal inertia and centrifugal effect of the compressor had great effects on the radial distribution of ice crystals.Nilamdeen et al.[16-17]extended the original supercooled water icing software and developed a mixed-phase icing calculation program for air,water droplets and ice crystals,and conducted the three-dimensional modeling on the ice crystal freezing in an engine to obtain the evolution process of ice shape on the surface of stator blades.

    However,relevant research in China is still in its infancy stage.Especially,the motion,heat transfer and icing characteristics of ice crystal ingested in the turbofan engine are rarely studied.Tan[18]used the Euler method to simulate the freezing of airfoil ice crystals,used the Messinger model to obtain the ice shape,and found that freezing is easier to occur in a low-pressure environment.Bu et al.[19]and Guo et al.[20]modified the Messinger model to calculate the icing of NACA0012 airfoil under mixed-phase conditions.The results showed that liquid water content,airflow temperature and ice crystal adhesion effect have great influences on the amount and profile of icing.Jiang et al.[21]used the Lagrange method to study the motion,heat transfer and adhesion of ice crystals to the engine.They obtained ice profile of the surface area of stator.Guo et al.[22]simulated the movement and heat transfer of ice crystals in a typical large icy wind tunnel configuration based on Euler method,and clarified the thermal/force balance characteristics of medium ice crystals in large icy wind tunnels.Ma et al.[23]developed an airfoil icing model on the NNWICE platform,and analyzed the influence of ice crystal shape and flow parameters on the phase transition and adhesion of moving ice crystals.

    Most of the reported studies simplify ice crystals as spherical particles.Ice crystals in clouds can undergo freezing and condensation,leading to complex shapes and size distributions.Therefore,in simulations of ice crystal freezing,it is crucial to take into account the impact of irregular particles which can affect the behavior of ice crystals.Naturally-formed ice crystals also have different melting phase transition processes under high-speed air flows,so it is necessary to compare and analyze different melting models.In this study,by using the Lagrange method,the dynamic behavior and heat and mass transfer of ice crystals after being ingested into an engine are investigated numerically.The heat transfer characteristics of ice crystals inside an aircraft engine are obtained.In the model,different melting models of crystals are compared,and the effects of shape irregularity of ice crystals and compressor blade thickness are examined.This study consolidates the foundation for further numerical simulation on ice crystal icing phenomena in aircraft engines.

    1 Numerical Modeling

    According to the two-phase flow theory,the icing condition of aircraft ice crystals is in the range of unidirectional coupled dilute two-phase flow,and the following assumptions are made.

    (1)Only consider the effect of air flow on ice crystals,and neglect the effect of ice crystals to air flow.

    (2)Ice crystals neither condense nor break in motion.

    (3)Neglect internal temperature variation of an ice crystal,and the viscosity and density remain unchanged during the movement.

    1.1 Motion and heat transfer model of ice crystals

    In the engine flow channel,the ice crystal particles are affected by a variety of forces,including viscous forceFD,gravity and buoyancyFG,Magnus forceFM,Basset forceFB,etc.According to Newton’s second law

    whereMiceis the mass of the ice crystal and duice/dtthe acceleration of the ice crystal.Miceis a function of the radius of ice crystals shown as

    whereρiceis the ice density andRicethe radius of ice crystals.

    Since ice density is much greater than that of air,and the size of an ice crystal is small,the viscous force is two to three orders of magnitude greater than gravity,only the viscous force is considered,and other forces such as gravity with smaller magnitude are negligible.Combined with the drag coefficient,the equation of motion can be expressed as

    whereρairis the air density,CDthe drag coefficient of ice crystals,uairthe airflow velocity,anduicethe ice crystal velocity.

    Lawson et al.[24]measured and generalized the size and shape of ice crystal particles in cirrus clouds.The classification of the ice crystal shape is shown in Fig.1(a).They are categorized into the following mathematical models:Spherical,ellipsoidal,columns,hexagonal plate ice crystals,as shown in Fig.1(b).

    Fig.1 Classification of the shape of ice crystal particles in clouds and simplifying mathematical models[24]

    The sphericity is defined byΦ=Ap/A,whereApis the surface area of the sphere with the same volume as the ice crystal,andAthe actual surface area of the ice crystal.Φis used to measure the degree to which the non-spherical ice crystal approximates the spherical shape,and then further correct the drag coefficient and heat and mass transfer coefficient of the movement.By comparing theCDcorrection methods of various non-spherical particles proposed in the studies[9,25-26],the empirical correlation of spherical,ellipsoidal and Columns is

    whereA(Φ)=exp(2.328 8-6.458 1+2.448 6Φ2);B(Φ)=0.096 4+0.556 5Φ;C(Φ)=exp(4.905-13.894 4Φ+18.422 2Φ2-10.259 9Φ3);D(Φ)=exp(1.468 1+12.258 4Φ-20.732 2Φ2+15.885 5Φ3).

    Calculation empirical correlation ofCDwith a hexagonal plate ice crystal

    whereReis the relative Reynolds number.

    The ice crystals that are sucked into the engine experience a temperature difference with the surrounding hot air.As a result,the outer surface of the ice crystal melts into a water film,which either falls off or covers the surface of the ice crystal and continues to evaporate.Due to this melting process,the diameter of the ice crystal particles decreases.This phase change melting process can be simplified into the following two idealized models.

    The first model is the “naked” ice particle model,where the ice crystal particles are always in contact with air without any external water layer.In this model,once the outer surface melts,the resulting water film is blown away from the surface of the ice particles.The second model is the water cover model,where the water film resulting from melting covers the outer surface of the ice crystal particles,as shown in Fig.2.

    Fig.2 Schematic diagram of two melting models

    The initial stage of both models involves the sublimation process of supercooled ice crystals,with a heat transfer relationship as

    wherehis the convective heat transfer coefficient on the surface of the ice crystal,Cpithe specific heat capacity of ice crystals,hmthe mass transfer coefficient on the surface of ice crystals,ρgathe density of water vapor in the air,ρgithe saturation vapor density on the surface of the ice crystal,Tairthe airflow temperature,Ticethe ice temperature,andLsgthe sublimation heat of ice.

    For the “naked” ice model,after the ice crystals begin to melt,all the heat absorbed from the air is used for the melting process,and the energy balance equation is

    whereLslis the melting heat of ice crystals.

    For the water cover model,after the ice crystals begin to melt,the heat absorbed from the air is used for the melting of ice crystals and the evaporation of the water film.The energy conservation relationship therefore is

    whereLlgis the heat of vaporization of water,Lslthe heat of melting of ice crystals,andMgthe mass of water film evaporated into vapor.

    The mass of the ice-water mixture is the summation of those of ice and water,which is also equivalent to the mass of the ice-water mixture at the initial moment of melting minus the mass of evaporated water vapor

    whereMmixis the mass of the ice-water mixture,Mwaterthe mass of water in the ice-water mixture,andMsmthe mass of the ice-water mixture at the initial moment of melting.

    The equivalent densityρmixof ice-water mixtures can be determined by the mass of ice and water as

    Assuming that the water film produced after the melting of ice crystals evenly covers the surface of the particles,the equivalent radiusRmixof the ice-water mixture can be calculated from the equivalent density

    Similarly,the heat transfer of non-spherical particles is corrected,and the number of non-spherical particlesNuis calculated as[26]

    wherePris the Prandtl number.

    The mass transfer coefficient of ice crystal or ice-water mixture is obtained by heat and mass transfer analogy as

    whereCpais the specific heat capacity of air,andLethe Lewis number.

    The density of ice crystals is 919.1 kg/m3in this study,and the thermophysical parameters of water and air including density,viscosity,specific heat capacity and thermal conductivity are calculated according to temperature interpolation.

    1.2 Engine flow channel model

    This study uses the typical low-pressure stage compressor model of Ref.[27]for calculation.Fig.3 shows a five-stage compressor planar blade grid,with each stage composed of rotors and stators.

    Fig.3 Schematic diagram of ice crystals passing through the interior of a low-pressure stage compressor[27]

    The NACA airfoil family thickness distribution is adopted[28].In order to facilitate the calculation,a polynomial expression is selected as

    wherec=5% is the relative thickness,andxandyare dimensionless horizontal and ordinate coordinates of the blade profile line relative to the chord length,respectively.

    Ref.[22]is used to perform triangulation on the inlet and outlet of the rotor.This provides the axial and circumferential velocities of airflow within the flow channel.Interpolation and calculation are then used to determine the temperature and pressure of the airflow at all stages.The calculation follows a discrete time and progressive space method,as depicted in Fig.4.

    Fig.4 Calculation flow chart

    2 Result Analysis

    2.1 Comparison of two melting models

    The “naked” ice particle model and the water-covered model are both simplified.It is necessary to compare the two models as they both exist in actual processes.

    In Fig.5,a comparison is made between the melting processes of the “naked” ice particle model and the water-covered model under the same high-temperature air flow environment and with the same initial conditions of ice crystals.The specific running parameters are:constant airflow environment withUair=50 m/s,Tair=50 ℃,h=10%;Rice=20 μm,Tice=-20 ℃ for ice crystal initial conditions.The melting trends of the two models are similar,but the “naked” ice particle model takes longer to completely melt under the same working conditions.Specifically,the 20 μm particle takes 0.038 s to completely melt in the “naked” ice particle model.This is because,in the “naked” ice particle model,the radius of the ice crystal decreases over time,resulting in a reduction of the surrounding air convection heat transfer area.This,in turn,it leads to lower convective heat transfer during the later stages of melting compared to earlier stages.In contrast,when covered with a water film,the outside of the ice particle is wrapped with a water film,which keeps the heat exchange area and heat exchange rate essentially constant.Therefore,the melting time is shorter in the water-covered model.

    Fig.5 Comparison of melting processes in the “naked” ice particle model and the water cover model

    Furthermore,it can be observed that,during the initial melting stages before 0.023 s,the “naked” ice particle model actually melts slightly faster than the water-covered model.This can be attributed to the fact that,in the water-covered model,the presence of the water film means that some of the convective heat transfer must be used to evaporate the water film.This results in less heat being available for melting the ice crystals.

    Therefore,in the early stages of melting,the presence of the water film slows down the melting rate due to the energy required for evaporation.However,in the middle and late stages of melting,the water film increases the convective heat transfer area such that it is much larger than that of the anhydrous film in the “naked” ice particle model.This leads to a more significant heat exchange,thereby accelerating the melting speed.Overall,the presence of the water film covering accelerates the melting of ice crystals.

    2.2 Effect of initial conditions of ice crystals on motion heat transfer

    Fig.6 shows the movement of ice crystals(equivalent volume diameter is 40 μm)of different shapes inside the compressor after being sucked in from the same position.The columns ice crystals collide in turn with the third and the fourth stage stator,with a liquid water content of 15.3% when colliding with the fourth stage stator.The hexagonal plate ice crystal collides with the third stage rotor and third stage stator with the liquid water content of 12.3%.

    Fig.6 Trajectory of ice crystals with different shapes after inhalation

    In Fig.7,the change in liquid water content along the axial position is shown after ice crystals of different shapes,but at the same initial position,are sucked into the compressor.The percentage content of liquid water mass of different shapes of ice crystals along the axial direction gradually increases.The liquid water content of spherical ice crystals is the largest at the same position,while the liquid water content of ellipsoidal and hexagonal ice plate crystals is similar.Specifically,at the outlet of the low-pressure compressor,the liquid water content of spherical,ellipsoidal,column-shaped,and hexagonal ice plate crystals is 49.05%,44.32%,29.74%,and 37.13%,respectively.

    Fig.7 Changes in liquid water content after inhalation of ice crystals with different shapes

    Fig.8 shows the speed variation of ice crystals along the axial positions after different shapes of ice crystals are sucked in from the same position.The gray dotted line is the axial velocity of air,and collisions between the first stage of rotor and stator are absent before 400 mm.The ice crystals continue to accelerate in the air flow,and their final speed approaches that of air.The subsequent velocity change follows the same trend as air.The occurrence of collision can also be found in the distribution of axial velocity,such as ellipsoidal ice crystals(red curve)violently colliding at the head of the third stage rotor.This results in greatly reduced axial velocity and an almost total loss of axial momentum due to a direct impact at the front end of the fifth stage rotor.

    Fig.8 Changes in axial velocity of ice crystals with different shapes after suction

    It is evident that there is a noticeable difference between spherical and non-spherical ice crystals in terms of their motion during the phase transition.Combining Figs.7 and 8,it can be found that when there is no collision effect,the speed of columnar particles increases the fastest,and that of spherical particles increases slower,because non-spherical particles have greater viscous force and greater acceleration in the air flow.Since spherical particles move more slowly in the air flow,the same axial displacement takes longer,which increases the heat and mass transferred between the hot air and the ice crystals,resulting in the maximum water content of spherical liquids and the lowest water content of columnar liquids.Due to the consistency of motion and melting laws,it can also be concluded that nonspherical characteristic indirectly affects the melting behavior through the particle trajectory and velocity in this study.

    Fig.9 shows the variations in heat transfer coefficient for three ice crystal shapes in the engine runner.The heat transfer coefficient of the three ice crystal shapes are basically the same following the axial position trend with oscillatory fluctuations.The heat transfer coefficient of spherical ice crystals is large,and that of column ice crystals is small,attributed by the large difference between spherical ice crystals and air velocity.This difference not only affects the movement of ice crystals but also influences heat exchange.At an axial position of approximately 0.8 m,the ellipsoidal ice crystal had a violent collision,resulting in a significant difference between its axial velocity and air velocity,and a sudden increase in the heat transfer coefficient.

    Fig.9 Changes in heat transfer coefficient after inhalation of ice crystals with different shapes

    Fig.10 shows trajectories of ice crystals as they move within the compressor at various radii after being ingested.The larger the radius of the ice crystal,the stronger the inertia which results in upward drift of the trajectory.Ice crystals with a radius of 10 μm(yellow trajectory)collide with the third rotor and fourth stage stator where the liquid water content reaches 86.96%.As a result,the possibility of adhesion of 10 μm ice crystals on the fourth stage stator is greater.Ice crystals with a radius of 20 μm(purple trajectory)collide with the third stage stator,the fourth stage rotor and the fifth stage rotor.The liquid water content at the collision point is 45.53%,so the ice crystal with a radius of 20 μm is more likely to adhere to the third stator.Ice crystals with a radius of 100 μm(blue trace)collide with the first stage rotor,the first stator,and the third stage rotor,with the liquid water content still less than 10% at the outlet of the compressor.Therefore,ice crystals with a radius of 100 μm may not adhere nor freeze inside the low-pressure stage compressor.

    Fig.10 Trajectory of ice crystals with different radii after inhalation

    Fig.11 shows the axial temperature variation inside the compressor after inhalation of ice crystals with different radii and the same initial temperature.Ice crystals are warmed up from supercooled to 0 ℃along the axis before starting to melt.The smaller the ice crystals,the faster the temperature increases to reach 0 ℃.Ice crystals with a large size,such as those with a radius of 100 μm(purple curve),travel farther inside the compressor and enter the fifth stage of rotor and stator before they start melting.Due to their low liquid water content,the likelihood of adhesion and freezing is also low.

    Fig.11 Temperature changes after inhalation of ice crystals with different radii

    Figs.12 and 13 respectively show the melting rate and percentage variation of liquid water mass in the compressor after inhalation of ice crystals with different radii.The melting rates of ice crystals with different radii remain stable after melting starts.The liquid water content also increases slowly after melting starts.However,ice crystals with different radii begin to melt at different locations.The smaller the size,the earlier the crystals begin to melt.The melting rate for a larger size is faster,which is due to the larger surface area for heat transfer.

    Fig.12 Melting rate of ice crystals with different radii changes after inhalation

    The 80 μm particles exhibit a unique behavior in the front of the third stage stator(at approximately 0.6 m),which is due to their violent collision with the front of the blade.This collision results in a great decrease in the axial velocity of ice crystals,causing its speed to be greatly different from air flow speed,thereby enhancing heat exchange.Furthermore,due to the decreased velocity,the heat transfer time is longer under the same displacement,resulting in a sudden increase in the temperature,melting speed and liquid water content of the particle.Fig.13 shows that the liquid water content of 20,40,60,80,and 100 μm at the outlet of the low-pressure compressor is 60.4%,15.5%,6.3%,6.3%,and 1.9%,respectively.The liquid water content of the 60 μm ice crystal is less than 10% as it reaches the outlet.Therefore,ice crystals with size greater than 60 μm are difficult to adhere and freeze in the flow channel of the low-pressure compressor.

    Fig.13 Percentage content of liquid water mass after inhalation of ice crystals with different radii changes

    3 Conclusions

    The study utilized a numerical model under the Lagrange framework to analyze the dynamic behavior and heat and mass transfer of ice crystals in the flow channel of an aeroengine after ingestion.The main conclusions are as follows:

    (1)Two melting models of ice crystal particles were compared and analyzed.The results showed that the “naked” ice particle model took longer to fully melt compared to the water-covered model.The 20 μm particle took 0.04 s to melt,while the“naked” ice particle melted faster before 0.023 s.

    (2)The liquid water content of spherical ice crystals was higher compared to other shapes under the same conditions.Specifically,at the outlet,the liquid water content of 20 μm spherical ice crystals was calculated to be 60.7%.On the other hand,the liquid water content of ellipsoidal and hexagonal plate ice crystals was close,with both shapes having a liquid water content around 40% at the outlet.

    (3)The initial diameter affected motion and heat transfer of ice crystals in the compressor.The smaller the size,the earlier it started to melt,and the higher the liquid water content.The 20 μm crystal’s liquid water content at the outlet was as high as 60.7%,while the 40 μm one was only 17.8% under the same condition.The liquid water content at the outlet of the 60 μm ice crystal was close to 10%,enabling difficulty for the crystal to adhere to the ice.

    婷婷六月久久综合丁香| 免费黄网站久久成人精品 | 男女视频在线观看网站免费| 18美女黄网站色大片免费观看| 少妇高潮的动态图| 亚洲精品影视一区二区三区av| 亚洲成av人片免费观看| 91av网一区二区| 他把我摸到了高潮在线观看| 亚洲精品乱码久久久v下载方式| 精品久久久久久久久久免费视频| 露出奶头的视频| 欧美日本视频| bbb黄色大片| 成人鲁丝片一二三区免费| 亚洲av日韩精品久久久久久密| 亚洲国产精品成人综合色| 我要搜黄色片| 亚洲精品456在线播放app | 91麻豆精品激情在线观看国产| a级一级毛片免费在线观看| 天堂av国产一区二区熟女人妻| 欧美日韩综合久久久久久 | 久久久久久久久久成人| 麻豆久久精品国产亚洲av| 国产视频内射| 九色国产91popny在线| 国内毛片毛片毛片毛片毛片| 久久九九热精品免费| 精品人妻1区二区| 久久久成人免费电影| 亚洲天堂国产精品一区在线| 在线播放无遮挡| 国产爱豆传媒在线观看| 人妻久久中文字幕网| 亚洲在线自拍视频| 91久久精品电影网| 欧美高清性xxxxhd video| 久久中文看片网| 国产精品自产拍在线观看55亚洲| 12—13女人毛片做爰片一| 日韩成人在线观看一区二区三区| 91在线精品国自产拍蜜月| 99热6这里只有精品| 精品久久久久久久人妻蜜臀av| 在线a可以看的网站| 热99在线观看视频| 日韩高清综合在线| 久久伊人香网站| 亚洲人成电影免费在线| 一二三四社区在线视频社区8| 国产毛片a区久久久久| 在线观看舔阴道视频| 激情在线观看视频在线高清| 免费无遮挡裸体视频| 永久网站在线| 成人毛片a级毛片在线播放| 此物有八面人人有两片| 久久精品国产自在天天线| 国产成年人精品一区二区| 欧美性感艳星| 一区二区三区免费毛片| 18禁黄网站禁片午夜丰满| 亚洲人成电影免费在线| 国产精品美女特级片免费视频播放器| 久久久久久国产a免费观看| 99久久99久久久精品蜜桃| 久久久久久久久中文| 久久天躁狠狠躁夜夜2o2o| 天堂动漫精品| 国产精品日韩av在线免费观看| 欧美黑人巨大hd| 午夜精品久久久久久毛片777| 如何舔出高潮| 成人美女网站在线观看视频| 国产亚洲精品综合一区在线观看| 亚洲人成网站在线播放欧美日韩| 久久精品人妻少妇| 久久亚洲精品不卡| 夜夜爽天天搞| 又黄又爽又刺激的免费视频.| 日本黄色视频三级网站网址| 亚洲精品色激情综合| 美女cb高潮喷水在线观看| 淫妇啪啪啪对白视频| 日韩欧美 国产精品| 国产精品三级大全| 国产单亲对白刺激| 一个人看视频在线观看www免费| 亚洲美女黄片视频| 日韩高清综合在线| 99精品在免费线老司机午夜| 亚洲美女黄片视频| 国产成人影院久久av| 又黄又爽又刺激的免费视频.| 欧美区成人在线视频| 久久婷婷人人爽人人干人人爱| 国产精品久久视频播放| 啦啦啦观看免费观看视频高清| 12—13女人毛片做爰片一| 欧美一区二区国产精品久久精品| 老司机福利观看| 亚洲色图av天堂| 久久99热这里只有精品18| 午夜老司机福利剧场| 最近最新中文字幕大全电影3| 亚洲电影在线观看av| 男女床上黄色一级片免费看| 亚洲欧美激情综合另类| 亚洲精品一卡2卡三卡4卡5卡| 国产私拍福利视频在线观看| 亚洲人与动物交配视频| 91狼人影院| 日韩中字成人| 久9热在线精品视频| 色综合婷婷激情| 美女大奶头视频| 69av精品久久久久久| 亚洲美女搞黄在线观看 | 高潮久久久久久久久久久不卡| 1000部很黄的大片| 51午夜福利影视在线观看| 国产精品永久免费网站| 亚洲精品亚洲一区二区| 日本a在线网址| 亚洲美女搞黄在线观看 | 国产精品乱码一区二三区的特点| 亚洲激情在线av| 男插女下体视频免费在线播放| 国产探花在线观看一区二区| 欧美+亚洲+日韩+国产| 超碰av人人做人人爽久久| 欧美成人免费av一区二区三区| 久久国产精品影院| 国产一区二区在线av高清观看| 99热6这里只有精品| 亚洲最大成人中文| 天堂√8在线中文| 热99re8久久精品国产| 91狼人影院| 国内揄拍国产精品人妻在线| 18禁黄网站禁片午夜丰满| 久久久久久久久久成人| 最新中文字幕久久久久| 99riav亚洲国产免费| 啦啦啦观看免费观看视频高清| 久久欧美精品欧美久久欧美| 夜夜爽天天搞| 午夜免费男女啪啪视频观看 | 国产爱豆传媒在线观看| 啪啪无遮挡十八禁网站| 观看免费一级毛片| 久久国产乱子免费精品| 久久午夜亚洲精品久久| 久久久久久久久久成人| 长腿黑丝高跟| 欧美极品一区二区三区四区| 日本免费一区二区三区高清不卡| 国模一区二区三区四区视频| 看免费av毛片| avwww免费| 国产亚洲av嫩草精品影院| 99久久精品一区二区三区| 国产野战对白在线观看| 亚洲五月婷婷丁香| 免费看光身美女| 亚洲综合色惰| 亚洲一区高清亚洲精品| 永久网站在线| 亚洲国产精品成人综合色| 精品久久久久久久久久免费视频| 日韩欧美精品v在线| 久久久久国内视频| 好看av亚洲va欧美ⅴa在| 久久精品国产亚洲av香蕉五月| 欧美色视频一区免费| 欧美日韩国产亚洲二区| 国产伦精品一区二区三区视频9| 一本一本综合久久| 内地一区二区视频在线| 我的老师免费观看完整版| 久久久精品欧美日韩精品| 一卡2卡三卡四卡精品乱码亚洲| 午夜两性在线视频| 一个人免费在线观看电影| 蜜桃久久精品国产亚洲av| 免费在线观看影片大全网站| 美女xxoo啪啪120秒动态图 | 亚洲熟妇熟女久久| 亚洲美女搞黄在线观看 | 国产亚洲精品久久久久久毛片| 国产三级黄色录像| 国产av一区在线观看免费| 99国产精品一区二区三区| 成人毛片a级毛片在线播放| 国产午夜精品久久久久久一区二区三区 | .国产精品久久| 熟女电影av网| 精品国产三级普通话版| 日本a在线网址| 国产精品亚洲av一区麻豆| 国产成+人综合+亚洲专区| 国产精品久久久久久亚洲av鲁大| 90打野战视频偷拍视频| 久99久视频精品免费| 国产大屁股一区二区在线视频| 亚洲欧美日韩东京热| 亚洲欧美日韩无卡精品| 天美传媒精品一区二区| 久久草成人影院| 国产乱人伦免费视频| 国产精品三级大全| 在线播放国产精品三级| 精华霜和精华液先用哪个| 男人的好看免费观看在线视频| 欧美高清成人免费视频www| 日韩欧美一区二区三区在线观看| 欧美黄色淫秽网站| 91午夜精品亚洲一区二区三区 | 99热这里只有精品一区| 国产男靠女视频免费网站| a在线观看视频网站| 中文字幕高清在线视频| 久久精品国产亚洲av香蕉五月| 午夜激情欧美在线| 日本黄色视频三级网站网址| 久久亚洲真实| 99国产精品一区二区三区| 天堂√8在线中文| 亚洲avbb在线观看| 人人妻人人澡欧美一区二区| 亚洲欧美日韩高清专用| 国产精品1区2区在线观看.| 亚洲最大成人av| 丰满人妻熟妇乱又伦精品不卡| 99国产极品粉嫩在线观看| 18美女黄网站色大片免费观看| 欧美日本视频| 国产一区二区亚洲精品在线观看| 国产aⅴ精品一区二区三区波| 午夜福利18| 久久性视频一级片| 国产精品99久久久久久久久| 有码 亚洲区| 国产精品永久免费网站| 无遮挡黄片免费观看| 久久性视频一级片| 国产老妇女一区| 俄罗斯特黄特色一大片| 99精品在免费线老司机午夜| 亚洲av免费在线观看| 精品日产1卡2卡| 好男人电影高清在线观看| 国模一区二区三区四区视频| 老熟妇乱子伦视频在线观看| 精品免费久久久久久久清纯| 中文资源天堂在线| 丰满乱子伦码专区| 美女被艹到高潮喷水动态| 久久这里只有精品中国| 久久草成人影院| 好男人电影高清在线观看| 久久久久国内视频| .国产精品久久| 免费大片18禁| 久久久久久久久中文| 97超视频在线观看视频| 直男gayav资源| 亚洲欧美日韩无卡精品| 久久99热6这里只有精品| av在线蜜桃| 亚洲国产精品成人综合色| 日韩欧美在线乱码| 日韩欧美三级三区| 亚洲天堂国产精品一区在线| 久久精品国产亚洲av香蕉五月| 色视频www国产| 国产一级毛片七仙女欲春2| 亚洲黑人精品在线| 黄色配什么色好看| 男女之事视频高清在线观看| 国产伦一二天堂av在线观看| 色综合站精品国产| 美女被艹到高潮喷水动态| 久久精品国产清高在天天线| 国产av一区在线观看免费| 97人妻精品一区二区三区麻豆| 一级作爱视频免费观看| 成人特级黄色片久久久久久久| 91久久精品国产一区二区成人| 搡老熟女国产l中国老女人| av欧美777| 午夜老司机福利剧场| 黄色丝袜av网址大全| 免费看光身美女| 亚洲精品在线观看二区| 国产免费av片在线观看野外av| 91在线精品国自产拍蜜月| 国产av一区在线观看免费| 国内精品美女久久久久久| 欧美性猛交黑人性爽| 国产精品亚洲av一区麻豆| 亚洲精品一卡2卡三卡4卡5卡| 亚洲熟妇中文字幕五十中出| 国内精品久久久久久久电影| 亚洲色图av天堂| 久久精品影院6| 国产日本99.免费观看| 欧美日韩瑟瑟在线播放| 香蕉av资源在线| 国产麻豆成人av免费视频| 伊人久久精品亚洲午夜| 精品午夜福利在线看| 国产成人啪精品午夜网站| 久久国产乱子免费精品| 亚洲第一欧美日韩一区二区三区| 91av网一区二区| 欧美色欧美亚洲另类二区| 成人毛片a级毛片在线播放| 国产亚洲av嫩草精品影院| 变态另类成人亚洲欧美熟女| 久久午夜福利片| 内射极品少妇av片p| 亚洲无线在线观看| 国产在线男女| 国产亚洲欧美在线一区二区| 欧美潮喷喷水| 琪琪午夜伦伦电影理论片6080| 亚洲人成网站高清观看| 88av欧美| 国产 一区 欧美 日韩| 男女下面进入的视频免费午夜| 首页视频小说图片口味搜索| 最近最新中文字幕大全电影3| 国产野战对白在线观看| 性欧美人与动物交配| 成人国产一区最新在线观看| 非洲黑人性xxxx精品又粗又长| 午夜福利高清视频| 午夜激情欧美在线| av天堂在线播放| 午夜激情福利司机影院| 亚洲真实伦在线观看| 中文字幕精品亚洲无线码一区| 国产精品久久久久久久电影| 九九热线精品视视频播放| 天堂影院成人在线观看| 国产免费av片在线观看野外av| 蜜桃久久精品国产亚洲av| 美女免费视频网站| 亚洲avbb在线观看| 久久久成人免费电影| 日本在线视频免费播放| 高清日韩中文字幕在线| 国产综合懂色| 99热只有精品国产| 国产成人a区在线观看| 丰满人妻一区二区三区视频av| 日韩欧美国产在线观看| 国产成年人精品一区二区| 国产色爽女视频免费观看| 免费电影在线观看免费观看| 久久人妻av系列| 一进一出好大好爽视频| 国产综合懂色| 亚洲av成人av| 亚洲中文日韩欧美视频| 国内毛片毛片毛片毛片毛片| 免费看光身美女| 日本a在线网址| 亚洲欧美日韩东京热| 69人妻影院| 99国产精品一区二区三区| 日本a在线网址| 人妻制服诱惑在线中文字幕| 日韩av在线大香蕉| 婷婷精品国产亚洲av| 丝袜美腿在线中文| 黄色日韩在线| 99riav亚洲国产免费| 熟女电影av网| 最近在线观看免费完整版| 悠悠久久av| 麻豆国产av国片精品| 一进一出抽搐gif免费好疼| 成人av在线播放网站| 12—13女人毛片做爰片一| 精品人妻熟女av久视频| 波多野结衣高清无吗| 午夜激情欧美在线| 婷婷精品国产亚洲av在线| 两人在一起打扑克的视频| 无遮挡黄片免费观看| 国产精品伦人一区二区| 欧美bdsm另类| 中文在线观看免费www的网站| 色综合亚洲欧美另类图片| 黄色丝袜av网址大全| 美女黄网站色视频| 中文字幕久久专区| 亚洲最大成人av| 国产成+人综合+亚洲专区| 嫩草影院新地址| 国产精品一区二区三区四区免费观看 | 亚洲av成人不卡在线观看播放网| 在线看三级毛片| 亚洲成人精品中文字幕电影| 免费黄网站久久成人精品 | 亚洲中文日韩欧美视频| 制服丝袜大香蕉在线| 久久久久国产精品人妻aⅴ院| av视频在线观看入口| 国产精品98久久久久久宅男小说| 高清毛片免费观看视频网站| 色哟哟哟哟哟哟| 熟女人妻精品中文字幕| av中文乱码字幕在线| 美女高潮的动态| 婷婷色综合大香蕉| 丰满人妻熟妇乱又伦精品不卡| 亚洲,欧美精品.| 99热精品在线国产| 18禁在线播放成人免费| 国产白丝娇喘喷水9色精品| 亚洲精华国产精华精| 又紧又爽又黄一区二区| av女优亚洲男人天堂| 精品人妻偷拍中文字幕| 两人在一起打扑克的视频| 一级a爱片免费观看的视频| 欧洲精品卡2卡3卡4卡5卡区| 国产精品免费一区二区三区在线| 国产精品野战在线观看| 国产黄色小视频在线观看| 又黄又爽又免费观看的视频| 亚洲欧美清纯卡通| 老熟妇仑乱视频hdxx| 国产日本99.免费观看| 美女免费视频网站| 中文在线观看免费www的网站| 一二三四社区在线视频社区8| 中文字幕高清在线视频| 69人妻影院| 亚洲专区国产一区二区| 一本一本综合久久| 亚洲精品一区av在线观看| 成熟少妇高潮喷水视频| 久久久久久久久大av| 怎么达到女性高潮| 亚洲狠狠婷婷综合久久图片| 伦理电影大哥的女人| 一边摸一边抽搐一进一小说| 九九在线视频观看精品| 精品欧美国产一区二区三| 国产精品综合久久久久久久免费| 97超级碰碰碰精品色视频在线观看| 国产真实伦视频高清在线观看 | 亚洲专区国产一区二区| 别揉我奶头~嗯~啊~动态视频| 久久精品人妻少妇| 如何舔出高潮| 美女高潮喷水抽搐中文字幕| netflix在线观看网站| 一夜夜www| 婷婷精品国产亚洲av| 久久午夜亚洲精品久久| 日本与韩国留学比较| 欧美极品一区二区三区四区| 亚洲18禁久久av| 成人精品一区二区免费| 久久久久久久精品吃奶| 精品一区二区三区人妻视频| 日本撒尿小便嘘嘘汇集6| 3wmmmm亚洲av在线观看| 国产成人a区在线观看| 噜噜噜噜噜久久久久久91| 日本熟妇午夜| 亚洲最大成人手机在线| 人人妻人人澡欧美一区二区| 在线播放国产精品三级| 久久久久久久午夜电影| 激情在线观看视频在线高清| 1024手机看黄色片| 日韩成人在线观看一区二区三区| 观看免费一级毛片| 99国产综合亚洲精品| 激情在线观看视频在线高清| 最近在线观看免费完整版| 午夜影院日韩av| 国内久久婷婷六月综合欲色啪| 亚洲18禁久久av| 搡老熟女国产l中国老女人| 亚洲精华国产精华精| 成人毛片a级毛片在线播放| 国产精品伦人一区二区| 亚洲av免费在线观看| netflix在线观看网站| 老鸭窝网址在线观看| 91久久精品电影网| 成熟少妇高潮喷水视频| 国产av不卡久久| 日韩国内少妇激情av| 美女cb高潮喷水在线观看| 亚洲精品在线观看二区| 国产精品永久免费网站| 九九热线精品视视频播放| 99久久九九国产精品国产免费| 亚洲欧美日韩高清专用| 18禁黄网站禁片免费观看直播| 亚洲天堂国产精品一区在线| 我要看日韩黄色一级片| .国产精品久久| 最后的刺客免费高清国语| 午夜福利在线观看免费完整高清在 | 久久久久久大精品| 精品久久久久久久久久免费视频| 五月玫瑰六月丁香| 99热精品在线国产| 精品乱码久久久久久99久播| 十八禁国产超污无遮挡网站| 亚洲性夜色夜夜综合| 麻豆成人av在线观看| 日韩 亚洲 欧美在线| 日本 av在线| 免费人成在线观看视频色| 两个人视频免费观看高清| av在线蜜桃| 精品午夜福利视频在线观看一区| 深夜精品福利| 国产免费男女视频| 国产亚洲av嫩草精品影院| 国产午夜精品久久久久久一区二区三区 | 国产av不卡久久| 天天一区二区日本电影三级| 最近在线观看免费完整版| 欧美在线黄色| 18禁裸乳无遮挡免费网站照片| 国内精品美女久久久久久| 男人狂女人下面高潮的视频| 十八禁网站免费在线| 成人高潮视频无遮挡免费网站| 亚洲欧美日韩高清专用| 午夜福利欧美成人| a级一级毛片免费在线观看| 欧美黑人巨大hd| 日韩精品青青久久久久久| 亚洲久久久久久中文字幕| 变态另类丝袜制服| h日本视频在线播放| 在线免费观看不下载黄p国产 | 俺也久久电影网| 午夜影院日韩av| 久久久精品大字幕| 成人午夜高清在线视频| 国产乱人视频| 亚洲第一电影网av| 亚洲最大成人av| 亚洲成人免费电影在线观看| xxxwww97欧美| 美女黄网站色视频| 日本熟妇午夜| 国产成人aa在线观看| 久久99热这里只有精品18| netflix在线观看网站| 亚洲欧美清纯卡通| 精华霜和精华液先用哪个| 欧美潮喷喷水| 一区二区三区免费毛片| 亚洲乱码一区二区免费版| 成人美女网站在线观看视频| 亚洲aⅴ乱码一区二区在线播放| 中文字幕av成人在线电影| 精品久久久久久久久久久久久| 18禁在线播放成人免费| 国内精品美女久久久久久| 亚洲av电影不卡..在线观看| 90打野战视频偷拍视频| 一区二区三区激情视频| 如何舔出高潮| a级毛片免费高清观看在线播放| 亚洲专区国产一区二区| 最新在线观看一区二区三区| 亚洲av成人av| 亚洲精品一卡2卡三卡4卡5卡| 两人在一起打扑克的视频| 天堂动漫精品| 亚洲av二区三区四区| 国产亚洲欧美98| 天堂√8在线中文| 亚洲精品一卡2卡三卡4卡5卡| 两人在一起打扑克的视频| 国产欧美日韩精品亚洲av| 国产精品98久久久久久宅男小说| 久久99热6这里只有精品| 国产欧美日韩精品亚洲av| 亚洲经典国产精华液单 | 国产欧美日韩一区二区三| 久久久久免费精品人妻一区二区| 国产欧美日韩一区二区精品| 久久婷婷人人爽人人干人人爱| 国产伦在线观看视频一区| 国产精品嫩草影院av在线观看 | 久久人人爽人人爽人人片va | 天美传媒精品一区二区| 久久热精品热| 日韩欧美国产在线观看| 岛国在线免费视频观看| 久久精品人妻少妇| av在线老鸭窝| 久久亚洲真实| 无遮挡黄片免费观看| 看片在线看免费视频| 一边摸一边抽搐一进一小说| 精华霜和精华液先用哪个| 日韩欧美在线乱码|