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

    Orthogonal analysis of the influencing factors of gas-solid two-phase jet particles

    2020-11-02 13:34:56SONGXinhuaLIUZhenfengLIXiaojieYANHonghao
    實驗流體力學 2020年5期
    關(guān)鍵詞:數(shù)量儀器顆粒

    SONG Xinhua, LIU Zhenfeng, LI Xiaojie, YAN Honghao

    (State Key Laboratory of structural analysis for industrial equipment, Dalian University of Technology, Dalian Liaoning 116024, China)

    Abstract:Since the gas-solid two-phase particle injection technology is widely used, it is of great engineering significance to study its influencing factors.Firstly, to inverstigate two factors, pre-ssure and mass, experiments are arranged by the orthogonal experimental design principle.Then, based on the CFD-DEM (Computational Fluid Dynamics / Discrete Element Method) model, the FLUENT software is used for numerical simulation, and the results are analyzed by range analysis and variance analysis.Finally, a set of jet experiments is designed to obtain the particle injection trajectory, and the photographs are processed by color histogram.The following conclusions are drawn: the influence of pressure is greater than that of mass by range analysis and variance analysis, and the mean value is the smallest when the pressure is 0.3 MPa and the mass is 2 g.When the pressure is 0.3 MPa and the mass is 2 g, the injection effect is optimized, and the shooting effect is the best, which is consistent with the results of numerical simulation.

    Keywords:gas-solid two-phase;orthogonal analysis;CFD-DEM;numerical simulation;color histogram

    0 Introduction

    Gas-solid two-phase particle injection is widely used in agriculture, industry and engineering.Dintwa et al.[1〗described a simulation model of a spinning disc fertilizer spreader, developed for incorporation into an automatic control system, and obtained the movement track and distribution characteristics of fertilizer throwing.Villette et al.[2]defined a mathe-matical model of fertilizer particle motion on a spinning disc, and the study results in an analytical solution for a concave disc equipped with pitched straight vanes.Lei et al.[3]designed a screw combined centralized fertilizer-feeding device for granular fertilizer and determined the key structural parameters of inclined spiral hole.Najjar et al.[4]developed a flexible simulation framework, in which multiphase flow computations are performed that include three-way coupling among phases (mixture-droplet-smoke), conservative coupling approach, and full heat release for the burning mechanisms.Liu and Xing[5]used Euler-Lagrangian method to simulate the gas-solid two-phase flow field in the thrust vector nozzle of solid rocket motor, and studied the effect of solid particles on the thrust vector performance of the nozzle.Xie et al.[6]obtained the structural design formula of pneumatic conveying gas-solid ejector by using the hydrodynamics and empirical formulas of gas-solid two-phase flow, and studied the influence of structure and parameters on the performance of gas-solid ejector.In the primary air passage of a new type of once-through burner of a coal-fired boiler, Huang used a three-dimensional particle dynamic analyzer (3D-PDA) to measure the gas-particle two-phase flow.The gas-solid two-phase velocity distri-bution, turbulent kinetic energy distribution, particle diameter and particle concentration distribution were obtained in the primary air passage[7].Li et al.[8]characterized the effect of oxidation temperatures on the microstructure and micro hardness of Ni/Al2O3composite coatings, as-cold-sprayed and vacuum heat-treated.Hua et al.used discrete phase model and spray rate model to simulate the release and flow process of fire extinguishing agent particles, and analyzed the movement characteristics of ultra-fine particle fire extinguishing agent with particle size of 1~12 μm under obstacles.And the reference valueR(t) and evaluationAwere introduced to quantitatively evaluate the movement characteristics of 1~12 μm ultrafine particle fire extinguishers[9].Li et al.[10]proposed a hybrid (deterministic/stochastic) fundamental model for the major physico-chemical processes involved in an industrial HVOF thermal spray process (Diamond Jet hybrid gun, Sulzer Metco, Westbury, NY, USA).Zhang et al.simu-lated the gas-solid ejector based on the particle trajectory model.The ejection performance was analyzed by using the number of particles ejected from the ejector, the number of remaining particles and the trajectory of particles.It is found that the inlet diameter of the contraction nozzle in the gas-solid ejector has a certain influence on the injection performance[11].Wang et al.[12]used FLUENT software to simulate the hydrodynamic characteristics of gas-solid two-phase flow in a boiling bed of titanium slag spouted bed, and analyzed the flow characteristics and zoning characteristics of gas-solid two-phase flow under different blowing speeds.Wang et al.used a three-component particle-dynamics anemometer to measure the characteristics of two-phase gas-particle flows in a primary air nozzle using a gas/particle two-phase test facility.The jet trajectory of primary air velocities, Reynolds stress and particle concentration profiles were obtained.Based on this, the advantages and disadvantages of increasing the offset angle were analyzed to optimize performance[13].Wang et al.[14]had presented a measurement method to obtain the solid phase particle size and the concentration of gas-solid two-phase flow, based on optical principle and image processing, and an experimental system was built.

    With the development of gas-solid two-phase particle injection technology, it is of great value to study the factors affecting particle injection.Many scholars have studied particle injection factors such as nozzle exit velocity, nozzle distance, nozzle size, injection angle[15-17].However, the factors of the injection pressure and particle quality are not studied.Therefore, based on previous research, this paper combines numerical simulation with orthogonal design theory, and analyzes the influence of spray pressure and particle quality on gas-solid two-phase particle injection.Moreover, it is proved by experiments that this method provides a theoretical reference for the practical application of injection technology.

    1 Fundamental theory

    1.1 Orthogonal experimental design principle

    Orthogonal experimental design is a multi-factor and multi-level design method.According to the orthogonality, some representative points are selected from the comprehensive experiment to carry out the experiment.These points have the advantages of uniform dispersion, neatness and comparability[18-20].In this paper, there are two factors, the injection pressure and the particle mass, and each factor has four levels.The specified parameters of orthogonal table used is listed in Table 1.

    Table 1 Specified parameters表1 具體參數(shù)的選取

    1.2 CFD-DEM model

    The CFD-DEM model uses the averaged Navier-Stokes equation to describe the motion of the gas phase, and uses Newton’s Second Law to describe the motion of the particle phase[21].The governing equation for the CFD-DEM model is[22]:

    The equation of motion of the particles:

    (1)

    Where,mais the particle quality,rais the particle position,Vais the particle volume,pis the pressure,βis the interphase drag coefficient,εgis the gas phase concentration,ugis the gas phase macro speed,vsis the speed of a single particle,gis the acceleration of gravity,Fcont,ais the contact force from the particles collision.Fvdw,ais van der Waals force,Tais the moment,Iais the moment of inertia,Ωais the rotational acceleration, andωais the rotational speed.

    The collision force between two particles is given by,

    Fab,n=-knδnnab-ηnvab,n

    δn=(Ra+Rb)-|rb-ra|

    (2)

    Where,Fab,nis the normal contact force when the particles collide,knis the direction stiffness coefficient,δnis the shape variable when the particles collide,nabis the normal unit vector,ηnis the normal damping coefficient,vab,nis the normal relative speed,RaandRbrepresent the particle radii of the particle a and b, respectively.raandrbare the positions of the particle a and b, respectively.

    The mass conservation equation and the momentum conservation equation for the gas phase are as follows,

    -εgp-Sp+·(εgτg)

    (3)

    Where,ρgis the gas phase density andτgis the gas phase stress.

    The gas-solid phase drag force is given by,

    (4)

    VCis the volume of the calculation grid, and the phase drag coefficientβis given by,

    (5)

    The ideal gas state equation is as follow,

    (6)

    Where,dpis the particle diameter,Mgis the molar mass of the gas,Ris the molar gas constant, andTgis the gas phase temperature.

    The gas phase stress constitutive relationship is as follow,

    (7)

    Where,μgis the shear viscosity of the gas phase, andIrepresents an identity matrix.

    2 Numerical simulation

    2.1 Fluid calculation model

    In this study, the average diameter of granules used in numerical simulation and experiment is 10.67 μm, spherical starch-like granules.When the pressure is 0.2~0.5 MPa, according to the Bernoulli equation[23], the energy conservation equation of one-dimensional defined flow under isentropic conditions is[24-25]:

    (8)

    The isentropic equation is:

    p=Aργ

    (9)

    According to Eq.(8) and Eq.(9), the jet velocity is:

    (10)

    For air,γ=1.4[24],ρ0=1.185 kg/m3[26],p0=0.1 MPa.Whenp=0.2 MPa, the substitution formula (10) calculatesv0.2=359.7 m/s.Similarly, we can getv0.3=466.7 m/s,v0.4=535.7 m/s,v0.5=587.3 m/s.In addition, the diameter of the air compressor outlet piped1=0.008 m, and the diameter of the nozzled2=0.05 m.v′0.2≈v0.2×(d1/d2)2≈9.2 m/s, and the effective injection velocitiesv′ are 9.2, 11.9, 13.7 and 15.0 m/s, respectively.Since the flow rate is less than 0.6 to 0.8 times of the speed of sound, the air can be regarded as an incompressible fluid during the simulation[27], and the reference atmospheric pressure is one atmosphere.The FLUENT is used to calculate the transient jet flow field under each pressure, and the CFD-DEM coupling method is used for the simulation calculation.

    In actual engineering, the air viscosity at room temperatureυ=15.06×10-6m2/s[28], and the equivalent diameterD=0.05 m.The Reynolds number is calculated:Re=vD/υ>4000, so the realizablek-εturbulence model equation is used[29].The kinetic energykand the turbulent dissipation transport equation are shown in Eq.(11):

    (11)

    2.2 Mesh and boundary conditions of model

    2.2.1Meshgeneration

    The trajectory of the particle jet into the box is studied.In order to facilitate the calculation, the model is calculated and analyzed by a two- dimensional model.The schematic diagram of the model is shown in Fig.1(a).The geometric model has two parts, the upper part is a rectangle of 1.0 m, and the lower part of the nozzle has a diameter of 0.05 m and a length of 0.1 m.

    Fig.1 Schematic diagram of the model

    In view of the uncomplicated characteristics of the particle injection model, a common grid structure is divided by the commonly used quadrilateral mesh, the mesh is divided by the workbench mesh, and the FLUENT software is imported for simulation calculation.The grid case is as follows: cell is 1268, face is 2611, and node is 1344.Fig.1(b) is a schematic diagram of the mesh division.

    2.2.2Boundaryconditions

    Each boundary condition is set according to the conditions of the sprayed particles.The initialization method selects “Standard Initialization” and selects “All-Zones”.The inlet boundary of the particle is set to “Velocity-Inlet”, and the boundary of the other face is set to “Wall”.The solid wall is no slip condition, the top wall is set to “Trap”, and the other walls are set to “Reflect”.

    2.3 Analysis and discussion of simulation results

    Particles were put into a rectangular frame with a diameter of 0.05 m and a length of 0.1 m.The simulation experiments were carried out according to the parameters in Table 1.The calculation steps were 500 and the particles were tracked.The simulation calculates the percentage of particle loss in each experiment, that is to say, the ratio of the number of tracked particles to the total number of tracked particles is the amount of particle loss in the simulation process.The data of the 16 computational simulations are shown in Table 2, and the percentage of loss is shown in Fig.2.

    Table 2 Tracking the number of particles trapping and the total number of particles (500 steps)表2 追蹤顆粒逃逸的數(shù)量和追蹤總顆粒數(shù)量(500步)

    Fig.2 Particle loss in orthogonal experiments

    By performing the range analysis of the 16 experiments, the mean value of pressure and mass was the smallest when the pressure was 0.3 MPa and the mass was 2 g.The range of the pressure and the mass are 20.1% and 12.4%, respectively, indicating that the magnitude of the pressure has the greatest impact.Then the variance analysis was performed on the two factors of the pressure and the mass, as shown in Table 3.TheFratios of the pressure and the mass obtained from the variance analysis table were 1.431 and 0.569, respectively, indicating that the influence of pressure was greater.

    Table 3 Variance analysis of pressure and quality impact 表3 壓強和質(zhì)量影響方差分析表

    3 Experimental verification

    3.1 Experimental device diagram

    This experiment mainly observes the flow trajectory of particles under different injection pressures (The dimensions of main instruments in the experiment are shown in Fig.1).The experimental diagram is shown in Fig.3, where A is a compressed air mercury, B is a gas storage chamber, C is a pulse valve, and D is an imaging system, E is where the particles are stored.Using the pulse valve C to eject particles at E.And the instruments used in this experiment are shown in Table 4.

    Table 4 Experimental measurement instrument list表4 實驗測量儀器清單

    Fig.3 Schematic diagram of injection experiment

    3.2 Analysis of experimental results

    Through the above numerical simulation calcu-lation and experimental orthogonal analysis, the results are obtained.When the pressure is 0.3 MPa and the mass is 2 g, the loss optimization is achieved.The experimental parameters were taken as the above results, and the spray experiment was carried out with the designed particle ejection device, and then the image was taken by the imaging system.

    Finally, the image color histogram processing was performed on the captured image[30-31], and the ejection trajectory diagram is shown in the left of Fig.4.Select the captured particle ejection images to perform image color histogram processing.Considering that the RGB color space does not conform to the human perception of color, the HSV space is visually sensed to quantify the HSV space and then its histogram is calculated to reduce the amount of calculation.The color space is converted from RGB to HSV by the rgb2hsv function in MATLAB.The HSV grayscale image obtained by processing the image is shown in the right of Fig.4.

    From Fig.2, we can see that Exp.1, Exp.6 and Exp.8 are the three groups with the lowest particle loss, which are 12.46%, 9.76% and 8.88% respectively.However, the optimal combination obtained by variance analysis is Exp.6.Then the injection experiment was taken on the spot when the pressure was 0.3 MPa and the mass was 2 g.At this time, the whole process of particle injection from the beginning to the top of the box could be captured, so the shooting effect was considered to be the best.The results are in agreement with those obtained by orthogonal experimental simulation, as shown in Fig.4.

    Fig.4 Numerical simulations of the sprayed particles and the trajectory of the shot (P=0.3 MPa, m=2 g)

    4 Conclusion

    The two factors, the pressure and the mass, which affect the gas-solid two-phase jet particles were studied and analyzed.The 16 simulation experiments were arranged based on the principle of orthogonal experimental design.Then, based on CFD-DEM model, FLUENT was used for numerical simulation, and the simulation results were analyzed by orthogonal test.Finally, a set of gas-solid injection device was designed to carry out the injection experiment, and the color histogram of the image was processed.The following conclusions are drawn:

    The range of the pressure and the mass is 20.1% and 12.4%, respectively.TheFratios of pressure and mass are 1.431 and 0.569, respectively.Both methods show that the effect of pressure is greater than that of mass, and the mean value is the smallest when the pressure is 0.3 MPa and the mass is 2 g.

    The image of the particle injection trajectory captured by the injection experiment was processed by color histogram and compared with the trajectory obtained by numerical simulation.The results show that when the pressure is 0.3 MPa and the mass is 2 g, the injection effect is optimized, and the shooting effect is the best.

    Acknowledgements:This project was financially supported by the National Natural Science Foundation of China (11672068, 11672067).

    猜你喜歡
    數(shù)量儀器顆粒
    《現(xiàn)代儀器與醫(yī)療》2022年征訂回執(zhí)
    《現(xiàn)代儀器與醫(yī)療》2022年征訂回執(zhí)
    Efficacy and safety of Mianyi granules (免疫Ⅱ顆粒) for reversal of immune nonresponse following antiretroviral therapy of human immunodeficiency virus-1:a randomized,double-blind,multi-center,placebo-controlled trial
    要讓顆粒都歸倉
    心聲歌刊(2019年1期)2019-05-09 03:21:32
    統(tǒng)一數(shù)量再比較
    我國古代的天文儀器
    疏風定喘顆粒輔料的篩選
    中成藥(2017年4期)2017-05-17 06:09:29
    頭發(fā)的數(shù)量
    我國博物館數(shù)量達4510家
    連花清瘟顆粒治療喉瘖30例
    丰满饥渴人妻一区二区三| 久久ye,这里只有精品| 老汉色∧v一级毛片| 免费在线观看黄色视频的| tocl精华| 成人国产av品久久久| 美女中出高潮动态图| 国产淫语在线视频| 成人免费观看视频高清| 欧美精品人与动牲交sv欧美| 在线观看免费午夜福利视频| 中国国产av一级| 热re99久久国产66热| 手机成人av网站| 18禁观看日本| 男女床上黄色一级片免费看| 欧美亚洲 丝袜 人妻 在线| 婷婷丁香在线五月| 中文字幕精品免费在线观看视频| av天堂在线播放| 亚洲成av片中文字幕在线观看| 亚洲av欧美aⅴ国产| 夜夜夜夜夜久久久久| 丰满迷人的少妇在线观看| 我要看黄色一级片免费的| 久久久久久久大尺度免费视频| 另类亚洲欧美激情| 视频在线观看一区二区三区| 久久久精品94久久精品| 国产精品一区二区精品视频观看| 在线观看一区二区三区激情| 国产男女超爽视频在线观看| 老司机午夜十八禁免费视频| 欧美国产精品一级二级三级| 亚洲av成人一区二区三| 久久久久久久国产电影| 岛国毛片在线播放| 免费在线观看黄色视频的| 欧美久久黑人一区二区| 另类精品久久| 两性午夜刺激爽爽歪歪视频在线观看 | 日韩 亚洲 欧美在线| 欧美日韩福利视频一区二区| 亚洲,欧美精品.| 国产高清视频在线播放一区 | 久久久欧美国产精品| 亚洲精品国产区一区二| 99国产精品99久久久久| 爱豆传媒免费全集在线观看| 嫁个100分男人电影在线观看| 国产精品九九99| 婷婷成人精品国产| 一区二区日韩欧美中文字幕| 中国美女看黄片| 中文精品一卡2卡3卡4更新| 欧美乱码精品一区二区三区| av在线老鸭窝| 三上悠亚av全集在线观看| 亚洲av男天堂| 菩萨蛮人人尽说江南好唐韦庄| 欧美人与性动交α欧美软件| 一边摸一边做爽爽视频免费| 男女床上黄色一级片免费看| 午夜两性在线视频| 午夜精品久久久久久毛片777| 51午夜福利影视在线观看| 国产亚洲欧美在线一区二区| 美女视频免费永久观看网站| 亚洲欧美激情在线| 我要看黄色一级片免费的| 婷婷丁香在线五月| 午夜91福利影院| 亚洲精华国产精华精| 成人亚洲精品一区在线观看| 亚洲国产精品999| 自线自在国产av| 亚洲精品国产区一区二| 人妻人人澡人人爽人人| 久久久久久久精品精品| 久久久久久久精品精品| 我的亚洲天堂| 搡老乐熟女国产| 国产精品亚洲av一区麻豆| 老司机靠b影院| 久久国产精品影院| 国产成+人综合+亚洲专区| 午夜日韩欧美国产| 久久99热这里只频精品6学生| 亚洲美女黄色视频免费看| 亚洲av欧美aⅴ国产| 一区二区日韩欧美中文字幕| 性色av一级| 五月天丁香电影| 久久精品国产综合久久久| 看免费av毛片| 高潮久久久久久久久久久不卡| 亚洲成国产人片在线观看| 中文欧美无线码| 亚洲av日韩精品久久久久久密| 国产精品久久久人人做人人爽| 国产精品久久久久久精品古装| 手机成人av网站| 午夜91福利影院| 色视频在线一区二区三区| 免费av中文字幕在线| 两人在一起打扑克的视频| 免费人妻精品一区二区三区视频| 国产精品.久久久| 日本欧美视频一区| 国产精品1区2区在线观看. | 一本大道久久a久久精品| 亚洲精品成人av观看孕妇| 久久亚洲精品不卡| 人成视频在线观看免费观看| 久久这里只有精品19| av国产精品久久久久影院| 韩国高清视频一区二区三区| 国产又色又爽无遮挡免| 在线看a的网站| 一区在线观看完整版| 国产成人免费观看mmmm| 国产成人欧美在线观看 | 亚洲国产日韩一区二区| 国产av精品麻豆| 十八禁人妻一区二区| 久久这里只有精品19| 多毛熟女@视频| 在线观看舔阴道视频| 亚洲欧美精品综合一区二区三区| 999久久久精品免费观看国产| 精品人妻一区二区三区麻豆| 五月开心婷婷网| 欧美在线一区亚洲| 免费观看a级毛片全部| 久久久国产一区二区| 国产99久久九九免费精品| 亚洲欧美日韩另类电影网站| 久久狼人影院| 国产成人精品无人区| 中文字幕色久视频| 久久精品亚洲熟妇少妇任你| 每晚都被弄得嗷嗷叫到高潮| 久久毛片免费看一区二区三区| 俄罗斯特黄特色一大片| 后天国语完整版免费观看| 正在播放国产对白刺激| 9热在线视频观看99| 在线观看www视频免费| 91麻豆精品激情在线观看国产 | 法律面前人人平等表现在哪些方面 | 老司机在亚洲福利影院| 亚洲成人免费av在线播放| 黄色视频不卡| 亚洲一区二区三区欧美精品| 国产亚洲欧美精品永久| 精品第一国产精品| 免费日韩欧美在线观看| 999精品在线视频| 1024香蕉在线观看| 亚洲性夜色夜夜综合| 欧美日韩一级在线毛片| 亚洲第一欧美日韩一区二区三区 | 老司机午夜福利在线观看视频 | 久久精品久久久久久噜噜老黄| 伦理电影免费视频| 咕卡用的链子| 女人被躁到高潮嗷嗷叫费观| 搡老岳熟女国产| 久久女婷五月综合色啪小说| tocl精华| svipshipincom国产片| 日韩免费高清中文字幕av| 国产日韩欧美视频二区| 精品人妻1区二区| 高清视频免费观看一区二区| 国产精品 国内视频| 久久久久视频综合| 久久 成人 亚洲| 色老头精品视频在线观看| 亚洲 国产 在线| 国产在线观看jvid| 一个人免费在线观看的高清视频 | 俄罗斯特黄特色一大片| 久久天躁狠狠躁夜夜2o2o| 精品欧美一区二区三区在线| 91精品国产国语对白视频| 国产亚洲av高清不卡| 天天影视国产精品| 国产在线免费精品| 美女大奶头黄色视频| 久久久久国产一级毛片高清牌| 在线观看舔阴道视频| 免费在线观看影片大全网站| 人人妻人人爽人人添夜夜欢视频| 欧美xxⅹ黑人| 亚洲国产av新网站| 国产区一区二久久| 91大片在线观看| av不卡在线播放| 嫩草影视91久久| 日韩精品免费视频一区二区三区| 91成年电影在线观看| 日韩有码中文字幕| 亚洲情色 制服丝袜| 九色亚洲精品在线播放| 一个人免费看片子| 大陆偷拍与自拍| 啪啪无遮挡十八禁网站| 久久久精品国产亚洲av高清涩受| av天堂久久9| 亚洲精品国产区一区二| 欧美日本中文国产一区发布| 欧美精品一区二区免费开放| 高清在线国产一区| a 毛片基地| 亚洲成人国产一区在线观看| 男女免费视频国产| 高清av免费在线| 亚洲欧美精品综合一区二区三区| 一二三四社区在线视频社区8| 午夜精品久久久久久毛片777| www.精华液| 水蜜桃什么品种好| 国产欧美日韩一区二区三 | 国产成人啪精品午夜网站| 真人做人爱边吃奶动态| 亚洲少妇的诱惑av| 亚洲av电影在线进入| 欧美亚洲 丝袜 人妻 在线| 91大片在线观看| 亚洲国产欧美日韩在线播放| 日本黄色日本黄色录像| 中文欧美无线码| 啦啦啦视频在线资源免费观看| 日韩 亚洲 欧美在线| 一个人免费看片子| 成年人免费黄色播放视频| 免费看十八禁软件| 国产精品亚洲av一区麻豆| 99久久99久久久精品蜜桃| 老司机靠b影院| 精品人妻熟女毛片av久久网站| 亚洲黑人精品在线| 自拍欧美九色日韩亚洲蝌蚪91| 久久久久精品国产欧美久久久 | 久久久久久久大尺度免费视频| 欧美黑人精品巨大| 黑丝袜美女国产一区| 亚洲美女黄色视频免费看| 首页视频小说图片口味搜索| 91精品国产国语对白视频| 另类精品久久| 美女午夜性视频免费| 午夜福利在线免费观看网站| 亚洲人成电影观看| 午夜成年电影在线免费观看| 亚洲av成人一区二区三| www.自偷自拍.com| 最近中文字幕2019免费版| 各种免费的搞黄视频| 日韩 亚洲 欧美在线| 亚洲第一青青草原| 99精品久久久久人妻精品| 日日摸夜夜添夜夜添小说| 日韩一区二区三区影片| 无遮挡黄片免费观看| 90打野战视频偷拍视频| 亚洲人成电影观看| 首页视频小说图片口味搜索| 男男h啪啪无遮挡| 秋霞在线观看毛片| 考比视频在线观看| 国产一区二区三区在线臀色熟女 | 国产区一区二久久| 超碰成人久久| 国产亚洲av片在线观看秒播厂| 国产一区二区三区综合在线观看| 久久久久久久久久久久大奶| 国产成人欧美在线观看 | 老汉色av国产亚洲站长工具| 亚洲七黄色美女视频| 桃花免费在线播放| 人妻人人澡人人爽人人| 婷婷色av中文字幕| 人妻一区二区av| 母亲3免费完整高清在线观看| 老司机午夜十八禁免费视频| 91精品伊人久久大香线蕉| 蜜桃国产av成人99| 侵犯人妻中文字幕一二三四区| 久9热在线精品视频| 桃花免费在线播放| 丝袜脚勾引网站| a级片在线免费高清观看视频| 极品人妻少妇av视频| 97精品久久久久久久久久精品| 91精品三级在线观看| 日韩大片免费观看网站| 99久久99久久久精品蜜桃| 在线亚洲精品国产二区图片欧美| 啦啦啦免费观看视频1| 欧美日韩亚洲综合一区二区三区_| 99香蕉大伊视频| 国产精品一区二区免费欧美 | 国产av一区二区精品久久| cao死你这个sao货| 老司机午夜福利在线观看视频 | 日本wwww免费看| 91老司机精品| 一区二区三区激情视频| avwww免费| 满18在线观看网站| 色婷婷av一区二区三区视频| 真人做人爱边吃奶动态| 亚洲精品第二区| 国产精品国产av在线观看| 免费高清在线观看视频在线观看| 99久久精品国产亚洲精品| 大片免费播放器 马上看| 欧美 日韩 精品 国产| a级片在线免费高清观看视频| 天天操日日干夜夜撸| 爱豆传媒免费全集在线观看| 精品国产国语对白av| 一个人免费看片子| 国产精品自产拍在线观看55亚洲 | 午夜免费成人在线视频| 国产免费现黄频在线看| 后天国语完整版免费观看| 久久久久久久精品精品| 婷婷成人精品国产| 国产在线免费精品| 最黄视频免费看| 国产日韩一区二区三区精品不卡| 黄色视频在线播放观看不卡| 国产不卡av网站在线观看| 欧美精品人与动牲交sv欧美| 欧美日韩亚洲国产一区二区在线观看 | 高清欧美精品videossex| 精品亚洲成国产av| 天天操日日干夜夜撸| 狂野欧美激情性xxxx| 在线 av 中文字幕| 国产精品久久久av美女十八| 黄频高清免费视频| av网站在线播放免费| 亚洲精品第二区| 一本一本久久a久久精品综合妖精| 国产av又大| 啦啦啦视频在线资源免费观看| 欧美+亚洲+日韩+国产| 成人免费观看视频高清| 菩萨蛮人人尽说江南好唐韦庄| 久久国产亚洲av麻豆专区| 十八禁高潮呻吟视频| 在线看a的网站| 成人影院久久| 国产成人欧美| 一个人免费在线观看的高清视频 | 亚洲av成人一区二区三| 午夜免费成人在线视频| 91av网站免费观看| 久久久精品94久久精品| 国产无遮挡羞羞视频在线观看| 精品少妇黑人巨大在线播放| 最新在线观看一区二区三区| 久久久久国产精品人妻一区二区| 国产一区二区激情短视频 | 欧美精品高潮呻吟av久久| 男女免费视频国产| 丰满人妻熟妇乱又伦精品不卡| 亚洲精品中文字幕在线视频| 免费女性裸体啪啪无遮挡网站| 天堂8中文在线网| 日本欧美视频一区| 久久香蕉激情| 我要看黄色一级片免费的| 1024香蕉在线观看| 国产欧美日韩综合在线一区二区| 国产成+人综合+亚洲专区| 久久精品国产亚洲av香蕉五月 | 亚洲熟女精品中文字幕| 亚洲成人手机| 日本欧美视频一区| 国产黄频视频在线观看| 99久久99久久久精品蜜桃| 久久久精品免费免费高清| 人人妻人人澡人人看| av免费在线观看网站| 久久亚洲国产成人精品v| 精品久久久久久电影网| 在线 av 中文字幕| 日韩熟女老妇一区二区性免费视频| 国产男女内射视频| 亚洲欧美精品综合一区二区三区| 久久亚洲国产成人精品v| 精品熟女少妇八av免费久了| 咕卡用的链子| 宅男免费午夜| 精品国产乱码久久久久久小说| 国产在视频线精品| 1024香蕉在线观看| 日本撒尿小便嘘嘘汇集6| 国产精品1区2区在线观看. | 国产成人欧美在线观看 | 亚洲专区国产一区二区| 男女床上黄色一级片免费看| 亚洲精品乱久久久久久| 午夜两性在线视频| 精品久久久久久久毛片微露脸 | 老熟妇乱子伦视频在线观看 | 中文字幕人妻熟女乱码| 侵犯人妻中文字幕一二三四区| 精品免费久久久久久久清纯 | 久久狼人影院| 久久精品亚洲av国产电影网| 热re99久久国产66热| 精品视频人人做人人爽| 欧美日韩视频精品一区| 精品久久久精品久久久| 午夜福利一区二区在线看| 桃红色精品国产亚洲av| 精品卡一卡二卡四卡免费| 亚洲第一青青草原| 久久久久视频综合| 另类精品久久| 午夜福利在线免费观看网站| 国产激情久久老熟女| 欧美日韩黄片免| 大型av网站在线播放| 一级毛片女人18水好多| 丝袜喷水一区| 18在线观看网站| 黑人欧美特级aaaaaa片| 免费少妇av软件| 美女主播在线视频| 精品国产国语对白av| 亚洲精品久久成人aⅴ小说| 日韩制服骚丝袜av| 亚洲精品第二区| 精品国产一区二区三区久久久樱花| 成人18禁高潮啪啪吃奶动态图| 美女午夜性视频免费| 国产97色在线日韩免费| 久久国产精品影院| 亚洲九九香蕉| 久久久精品94久久精品| 一本大道久久a久久精品| 日本欧美视频一区| 久久久久国内视频| 男女高潮啪啪啪动态图| 欧美日韩亚洲高清精品| 五月开心婷婷网| 精品欧美一区二区三区在线| 国产精品一二三区在线看| 狠狠狠狠99中文字幕| av天堂久久9| 午夜福利视频在线观看免费| 免费高清在线观看视频在线观看| 国产区一区二久久| 欧美少妇被猛烈插入视频| 韩国精品一区二区三区| 亚洲精品自拍成人| 91九色精品人成在线观看| 十八禁网站免费在线| svipshipincom国产片| 精品亚洲乱码少妇综合久久| 香蕉国产在线看| 欧美精品亚洲一区二区| 国产av精品麻豆| 日本av手机在线免费观看| h视频一区二区三区| 性少妇av在线| 少妇粗大呻吟视频| 亚洲av国产av综合av卡| 超碰成人久久| 亚洲av日韩精品久久久久久密| 久久影院123| 精品一区二区三卡| 99热网站在线观看| 欧美日韩亚洲高清精品| 中国美女看黄片| 两性夫妻黄色片| 久久精品亚洲熟妇少妇任你| 视频区欧美日本亚洲| 91大片在线观看| 考比视频在线观看| 可以免费在线观看a视频的电影网站| 国产高清视频在线播放一区 | 久久国产亚洲av麻豆专区| 亚洲精品一二三| 欧美一级毛片孕妇| 精品亚洲乱码少妇综合久久| 日韩有码中文字幕| 黄网站色视频无遮挡免费观看| 69av精品久久久久久 | 超碰成人久久| 午夜福利在线免费观看网站| 精品国内亚洲2022精品成人 | 亚洲欧美成人综合另类久久久| 黑人巨大精品欧美一区二区mp4| 欧美在线一区亚洲| 91九色精品人成在线观看| www.av在线官网国产| 亚洲av电影在线观看一区二区三区| 久久精品aⅴ一区二区三区四区| 精品欧美一区二区三区在线| 女人高潮潮喷娇喘18禁视频| 菩萨蛮人人尽说江南好唐韦庄| 亚洲精品国产av蜜桃| 亚洲美女黄色视频免费看| 黑人巨大精品欧美一区二区mp4| 天天操日日干夜夜撸| 岛国在线观看网站| 国产免费福利视频在线观看| netflix在线观看网站| 日韩三级视频一区二区三区| 我的亚洲天堂| 亚洲,欧美精品.| av天堂久久9| 亚洲第一青青草原| 国产人伦9x9x在线观看| 国产主播在线观看一区二区| 中文精品一卡2卡3卡4更新| 日韩电影二区| 色综合欧美亚洲国产小说| 国产一区有黄有色的免费视频| 免费不卡黄色视频| 久久久久久久国产电影| 操出白浆在线播放| 欧美亚洲日本最大视频资源| 男人操女人黄网站| 少妇的丰满在线观看| 新久久久久国产一级毛片| 国产成人av教育| 国产日韩欧美视频二区| 视频在线观看一区二区三区| av免费在线观看网站| 欧美精品人与动牲交sv欧美| 欧美日韩av久久| 亚洲av欧美aⅴ国产| 老熟妇乱子伦视频在线观看 | 亚洲精品国产一区二区精华液| 国产精品久久久av美女十八| 亚洲成国产人片在线观看| 精品人妻一区二区三区麻豆| 亚洲国产av新网站| 国产日韩欧美亚洲二区| 18禁裸乳无遮挡动漫免费视频| 久久久国产欧美日韩av| 大码成人一级视频| 日本五十路高清| 一区在线观看完整版| 丝瓜视频免费看黄片| 日韩一卡2卡3卡4卡2021年| www.自偷自拍.com| 极品人妻少妇av视频| 欧美+亚洲+日韩+国产| 欧美日韩福利视频一区二区| 国产高清视频在线播放一区 | 老熟妇仑乱视频hdxx| 男人操女人黄网站| 好男人电影高清在线观看| 日韩制服骚丝袜av| 国产在线观看jvid| 亚洲综合色网址| 成年动漫av网址| 免费久久久久久久精品成人欧美视频| 久久精品熟女亚洲av麻豆精品| 欧美另类亚洲清纯唯美| 亚洲精品久久成人aⅴ小说| 国产日韩一区二区三区精品不卡| 亚洲中文av在线| 我要看黄色一级片免费的| 亚洲精华国产精华精| 精品国产超薄肉色丝袜足j| 精品国产乱码久久久久久小说| 十八禁高潮呻吟视频| 日韩大片免费观看网站| 青春草视频在线免费观看| 国产男人的电影天堂91| 青春草视频在线免费观看| 2018国产大陆天天弄谢| 国产xxxxx性猛交| 久久久国产一区二区| 欧美在线黄色| 午夜激情久久久久久久| tube8黄色片| 自拍欧美九色日韩亚洲蝌蚪91| 精品一区二区三卡| 精品欧美一区二区三区在线| 久久久久久久国产电影| 99热国产这里只有精品6| 女人高潮潮喷娇喘18禁视频| 可以免费在线观看a视频的电影网站| 日韩 欧美 亚洲 中文字幕| 免费av中文字幕在线| 久久av网站| 看免费av毛片| 国产一区二区激情短视频 | 国产又色又爽无遮挡免| 一级毛片电影观看| 久久天堂一区二区三区四区| 成人影院久久| 捣出白浆h1v1| 91精品国产国语对白视频| 免费高清在线观看日韩| 99国产精品一区二区三区| 午夜福利在线观看吧| 久久久久国产精品人妻一区二区| 亚洲精品久久久久久婷婷小说| 建设人人有责人人尽责人人享有的| 亚洲人成电影观看| 丝袜在线中文字幕| 国产麻豆69|