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

    Experimental study on plasma actuation characteristics of nanosecond pulsed dielectric barrier discharge

    2022-02-15 11:08:40HaoZHENG鄭浩HuaLIANG梁華JieCHEN陳杰HaohuaZONG宗豪華XiangzheMENG孟祥喆LikeXIE謝理科andYinghongLI李應紅
    Plasma Science and Technology 2022年1期
    關鍵詞:陳杰理科

    Hao ZHENG (鄭浩), Hua LIANG (梁華), Jie CHEN (陳杰),Haohua ZONG (宗豪華), Xiangzhe MENG (孟祥喆), Like XIE (謝理科) and Yinghong LI (李應紅)

    1 Science and Technology on Plasma Dynamics Laboratory, Air Force Engineering University, Xi’an 710038, People’s Republic of China

    2 Chinese Flight Test Establishment, Xi’an 710089, People’s Republic of China

    3 School of Mechanical Engineering,Xi’an Jiaotong University,Xi’an 710049,People’s Republic of China

    Abstract Combining high-speed schlieren technology and infrared imaging technology, related research has been carried out on the influence of parameters such as actuation voltage, repetition frequency, and electrode size of an actuator on the discharge characteristics, induced flow field characteristics,and thermal characteristics of nanosecond pulsed dielectric barrier discharge.The results show that increasing the value of the actuation voltage can significantly increase the actuation intensity, and the plasma discharge area is significantly extended.Increasing the repetition frequency can increase the number of discharges per unit time.Both will cause more energy input and induce more changes in the flow field.The effect of temperature rise is more significant.The width of the covered electrode will affect the potential distribution during the discharge process, which in turn will affect the extension process of the plasma discharge filament.Under the same actuation intensity,the wider the covered electrode,the larger range the induced flow field and temperature rise is.Preliminary experimental analyses of high-frequency actuation characteristics,temperature field characteristics,flow field characteristics and actuation parameter settings provide support for the parameter selection and partial mechanism analysis of plasma anti-icing.

    Keywords: NS-DBD, actuation voltage, repetition frequency, induced flow field,temperature rise

    1.Introduction

    Nowadays, the dielectric barrier discharge (DBD) based on sinusoidal alternating current high frequency and high voltage actuation (AC-DBD) and dielectric barrier discharge based on nanosecond pulsed voltage actuation (NS-DBD) are the most studied in the plasma application field[1–4].These two discharge methods are widely used in flow field control.The first to attract researchers is AC-DBD,which can induce a wall jet with a speed ranging from 0 to 11 m s?1,and it can control the flow separation of the flow field with an incoming flow velocity up to 100 m s?1[5–7].Later,studies have shown that NS-DBD has stronger flow control capabilities, so in recent years it has also attracted more and more attention from researchers.For NS-DBD, the flow control mechanism is mainly the rapid heat release effect caused by the rapid heating of local air [8–10].

    During the study of the actuated discharge characteristics of NS-DBD plasma, the influence of the actuation parameters, the configuration of the actuator,the type of the insulating dielectric layer on the energy consumption, the voltage and current waveform, and the discharge modes are mostly discussed.Jiang Hui et al[11–14]studied the influence of actuation voltage peak, actuation frequency, electrode width, and electrode spacing on the discharge characteristics of NS-DBD.The results showed that increasing the actuation voltage could make the discharge filament more uniform, and meanwhile the plasma filament length would also increase.Zhou Yang et al [15]compared the discharge characteristics of nanosecond pulsed plasma actuation and microsecond pulsed plasma actuation at different actuation voltages.The results showed that with the actuation voltage amplitude increasing, the current peak value and energy consumption of the surface DBD at the two actuation conditions increase.Ndong et al [16] conducted experimental studies on the voltage and current waveforms and the changes of single pulsed energy at different electrical parameters.Compared with the negative polarity pulsed voltage, the positive polarity pulsed voltage actuation can bring greater energy deposition.Takashima et al [17] used high-speed schlieren technology to observe the flow field characteristics of the NS-DBD plasma actuation from the lateral and longitudinal sides of the actuator respectively and observed the shock wave structure.Zhao et al[18,19]used a dynamic pressure sensor with a higher response frequency and a two-dimensional phase-locked PIV to measure the flow field characteristics of NS-DBD actuation.The NSDBD plasma actuation can induce an instantaneous flow of the flow field, and a continuous vortex with a speed less than 0.3 m s?1will be generated near the wall area.The heat generated by NS-DBD plasma actuation can be divided into two parts:the rapid heating of the air near the wall and the heating of the dielectric layer material.At the same time,there is a process of convective heat transfer between the dielectric layer material and the heated air on the upper layer[20].Unfer et al[21]simulated and calculated that nanosecond pulsed plasma actuation would instantaneously (at the time scale of a nanosecond) heat the air near the wall to about 1000 K.The simulation results of Zhu Yifei et al [22, 23] showed that the near-wall air would be rapidly heated to 1100 K and then rapidly dropped to 310 K at a speed of k = ?225/t Ks?1, and the heating area was approximately equal to the extension area of plasma discharge filaments.The simulation results of Zhu Yifei et al showed that the nearwall air would be rapidly heated to 1100 K and then rapidly dropped to 310 K at a speed of k = ?225/t Ks?1and the heating area was approximately equal to the extension area of plasma discharge filaments.Popov et al [24] proposed a rapid heating model of breakdown discharge in the nitrogen–oxygen mixture environment.Under different reduced field intensities,the physical and chemical reactions that play a major role in gas heating are the dissociation reaction of oxygen caused by electron collision and the dissociation reaction of nitrogen molecules under the action of electron collision.

    Researchers have studied the effects of peak actuation voltage, pulsed actuation frequency, and electrode size of the actuator on the discharge characteristics and induced flow field characteristics of NS-DBD.However, there are few researches on the infrared radiant thermal characteristics of NS-DBD.With the development of NS-DBD plasma actuation,it can not only be used for flow control,but also has application prospects in the field of aircraft anti-icing[25,26].Therefore, the law of its thermal effect and induced flow field effect with actuation parameters is worthy of in-depth research.

    Figure 1.Surface dielectric barrier discharge plasma actuator.

    In this work, the effects of the pulsed voltage, actuation frequency and the width of the covered electrode on the discharge characteristics, flow field characteristics and infrared radiation thermal characteristics of NS-DBD have been studied,and the influence of various parameters on the actuation characteristics of NS-DBD plasma has been obtained.

    2.Experimental setup

    2.1.Plasma actuator and power supply

    The configuration of the plasma actuator is shown in figure 1,including an upper electrode,a covering electrode,an insulating medium (Kapton), and a base plate.Copper foil with a thickness of 0.027 mm is used as the electrode material of the actuator.The upper electrode is exposed in the air and the width is 5 mm.The width of the covering electrode w is selected as three sizes:5 mm,10 mm,and 20 mm respectively.The lengths of the upper electrode and the covering electrode are both 60 mm.The dielectric layer is made of polyimide tape with a thickness of 0.18 mm and a dielectric constant of 3.4.The substrate plate is an aluminum nitride ceramic sheet with an area of 100 mm × 100 mm and a thickness of 1 mm.

    The nanosecond pulsed power supply (Parametric Highvoltage Pulsed Power Supply, Xi’an Smart Maple Electronic Technology Co.,Ltd,Xi’an,China)is used in the experiment.The output voltage waveform is shown in figure 2.The peak output voltage range of the power supply is continuously adjustable from 0 to 20 kV.The pulsed repetition frequency is continuously adjustable from 1 Hz to 20 kHz.the pulsed width tpis adjustable from 0 ns to 1 ms, and the rising edge time trand falling edge time tdare both adjustable.

    2.2.Discharge characteristic test system

    Figure 2.Schematic diagram of the output voltage waveform of the nanosecond pulsed power supply.

    Figure 3.Layout of NS-DBD electrical characteristics test system.

    The NS-DBD plasma actuation discharge characteristic test system mainly includes a voltage probe, a current probe, an oscilloscope, and a CCD camera, as shown in figure 3.The voltage and current applied to the upper electrode and the covering electrode of the actuator were measured by a voltage probe (Tektronix, P6015A) and a current probe (Tektronix,TCP0030A), respectively.The accuracy of the current probe is 1 mA.A digital oscilloscope(Tektronix,DPO4104)is used to display and store voltage–current signals at a sampling frequency of 5 GHz.The discharge image was taken by a Nikon D7000 digital camera with an exposure time of 0.04 s.The single pulsed energy of NS-DBD was computed by the integration of product the voltage and current acquired through oscilloscope over the single pulsed time: =W

    2.3.Flow field characteristic test system

    The high-speed schlieren technology is used to observe the induced flow field of NS-DBD plasma actuation in still air.As shown in figure 4,the high-speed schlieren test system mainly includes a xenon lamp, lens, knife edge, and high-speed camera.The xenon lamp is used as the light source to produce a parallel light path between the two concave lenses.The NSDBD plasma actuator is placed in the parallel light path to observe the NS-DBD plasma actuation induced flow field.A high-speed camera (Phantom—V2512) was used to capture the schlieren image induced by NS-DBD actuation.The camera frame rate is 40 000 fps.During the experiment, the resolution is 1024 × 1024 pixel2, and the camera exposure time is 0.000 002 s.The field of view for schlieren is 25 mm high and 30 mm long, and the actuator is perpendicular to it and parallel to the optical path.Through the DG535 timedelay pulsed generator,the high-voltage pulsed power supply,and the high-speed camera are triggered synchronously to shoot.

    Figure 4.High-speed schlieren experimental system diagram.

    Figure 5.Schematic diagram of positive and negative peak current and single pulsed energy change with voltage.

    2.4.Infrared radiant heat characteristic test system

    The infrared thermal imager used in the experiment is FLIR SC7300M, and the temperature measurement range is from?20 °C to 2000 °C.The infrared thermal imager is close to the surface of the actuator, so the attenuation of radiated power caused by the absorption and scattering of various components in the air can be ignored.The emissivity of polyimide tape was set as 0.83.The sampling frequency of the infrared thermal imager in the experiment is 25 Hz.

    3.Experimental results and analyses

    3.1.Influence of pulsed voltage peak

    The pulsed actuation frequency of 2 kHz, rising edge time of 150 ns,a pulsed width of 100 ns,falling edge time of 150 ns,and covering electrode width of 20 mm were selected to remain unchanged, and the output voltage peak value was changed.

    Figure 6.Discharge images of different voltage peaks (exposure time 0.04 s).

    Figure 7.Variation process of the flow field induced by nanosecond pulsed actuation with different pulsed voltage peaks: (a) 8 kV, (b)10 kV, (c) 12 kV, (d) 14 kV.

    Figure 5 shows the variation of single pulsed energy and the peak value of positive and negative current with the increase of the pulsed voltage peak value.It can be seen that with the increase of the peak value of the pulsed voltage, the peak value of the positive and negative pulsed current and the single pulsed energy both increase nonlinearly, which is consistent with the research results in the literature [27].

    Figure 6 shows the discharge images at different voltage peaks.It can be seen that the voltage peak has a greater impact on the discharge area and discharge pattern.When the peak voltage is 8 kV, the discharge is mainly concentrated in the range less than 2 mm from the edge of the electrode, the discharge area is small, and the brightness is darker from the discharge image.With the increase of the voltage peak value,the discharge filaments continue to grow, the discharge area gradually increases, and the discharge becomes more intense and the brightness is greater.Therefore, a larger range of plasma discharge can be obtained by increasing the peak value of the pulsed voltage.

    Figure 8.Surface temperature distribution of the actuator after continuous discharge of 120 s(t = 120 s)with different pulsed voltage peaks.

    Figure 7 shows the change process of the flow field induced by nanosecond pulsed actuation with different pulsed voltage peaks.It can be seen that the greater the peak value of the pulsed voltage, the greater the actuation intensity and the greater the disturbance to the still air.In addition, it can be concluded that by observing the change in the length of the tail as the peak value of the pulsed voltage increases,the tail of the induced shock wave structure becomes longer.Combined with the discharge image, it can be obtained that the length of the tail of the shock wave structure is related to the length of the discharge filament, which is consistent with the previous conclusions.

    Figure 9.The surface temperature changes under different pulsed voltage peak values along the X-axis at Y = 4 mm direction(t = 120 s).

    Figure 10.Changes of surface temperature along the Y-axis at X = 30 mm at different pulsed voltage peaks (t = 120 s).

    Figure 8 shows the distribution of the surface temperature of the actuator with different pulsed voltage peaks at the time of 120 s.It can be seen that as the peak value of the pulsed voltage increases, the surface temperature of the dielectric layer will increase significantly,and the high-temperature area will also increase significantly.

    Figure 9 shows the changes of the surface temperature along the X-axis at Y = 4 mm under different pulsed voltage peak conditions (t = 120 s).When the pulsed voltage peak value increases from 8 to 10 kV, the surface temperature increases by about 5 °C.When the peak value of the pulsed voltage increases from 10 to 12 kV,the temperature increases by nearly 30 °C.

    Figure 11.Variation of surface temperature with time at point A(X = 30 mm, Y = 4 mm) at different pulsed voltage peaks.

    The results in figures 10 and 11 also show that the surface temperature significantly increases with the increase of the pulsed voltage peak value.It can be seen from the graph of surface temperature changing with time in figure 11 that the higher the peak value of the pulsed voltage,the faster the surface temperature rises.Besides, it should be pointed out that in figure 11, when the pulsed voltage peak value is 14 kV, the sudden drop in temperature at the time of 72 s is caused by the breakdown of the actuator.The insulating medium (Kapton) near the upper electrode reaches the ignition point, causing the insulating material to ignite and breakdown.A plasma arc was generated at the location where it was pierced, and a strong current was formed, causing the circuit to short-circuit.

    Figure 12.Schematic diagram of positive and negative peak current and single pulsed energy change with frequency.

    Figure 13.Discharge images of different pulsed frequencies (exposure time 0.04 s).

    Figure 14.Variation process of the flow field induced by nanosecond pulsed actuation with different pulsed actuation frequencies: (a)500 Hz, (b) 1000 Hz, (c) 2000 Hz, (d) 4000 Hz.

    Increasing the pulsed voltage peak increases the actuation intensity of the single pulsed energy,so the input energy is also very high.The intensity of the electric field increases as the voltage increases, so that the gas far away from the upper electrode can also be broken down.Thus,the discharge filament becomes longer with the increase of voltage.The higher the energy is, the faster the surrounding flow field flows, the faster the temperature rises and the wider the actuation range is.

    3.2.Influence of pulsed actuation frequency

    The pulsed voltage peak value is 12 kV,the rising edge time is 150 ns, the pulsed width is 100 ns, the falling edge time is 150 ns,and the covering electrode width is 20 mm.Further study was carried out by changing the pulsed actuation frequency.

    Figure 12 shows the variation of the single pulsed energy and the peak value of positive and negative current with the increase of pulsed actuation frequency.It can be seen that the change of the pulsed actuation frequency has almost no effect on the peak value of the positive and negative current and the single pulsed energy,indicating that the change of the pulsed actuation frequency only changes the number of discharges per unit time, but has no effect on the single discharge process.This conclusion is consistent with the conclusion in the literature[25],but is not consistent with the conclusion in the literature[17].The difference in experimental results is due to the different internal structures and design principles of the high-voltage pulsed power supply.

    Figure 15.Surface temperature distribution of the actuator after continuous discharge of 120 s (t = 120 s) with different pulsed actuation frequencies.

    Figure 16.The surface temperature changes under different pulsed actuation frequencies along the X-axis at Y = 4 mm direction(t = 120 s).

    Figure 17.Changes of surface temperature along the Y-axis at X = 30 mm at different pulsed actuation frequencies (t = 120 s).

    Figure 18.Variation of surface temperature with time at point A(X = 30 mm, Y = 4 mm) at different pulsed actuation frequencies.

    Figure 19.Variation of single pulsed energy with rising edge time,falling edge time,and pulsed width with different covered electrode widths.

    Figure 13 shows the discharge images of different pulsed actuation frequencies.By comparison,as the pulsed frequency increases, the filament discharge brightness gradually increases, but the filament length and discharge area hardly change.The reason is that the actuation frequency does not affect the single discharge process.With the increase of frequency, the discharge times per unit time increase,which is reflected in the discharge images of multiple pulsed integrals,so the brightness increases.

    Figure 20.Discharge images with different covered electrode widths (exposure time 0.04 s).

    Figure 14 shows the change process of the flow field induced by nanosecond pulsed actuation with different pulsed actuation frequencies.It can be seen from the figure that the single shock wave structure almost does not change with the increase of pulsed actuation frequency, which is consistent with the previous conclusion that changing the pulsed actuation frequency has almost no effect on a single discharge.Besides, as the pulsed actuation frequency increases,the more discharge times per unit time, the faster the energy accumulation,and the faster the induced flow field tends to a steady state.

    Figure 15 shows the surface temperature distribution cloud diagram under different pulsed actuation frequencies after the discharge time is 120 s.Figures 16 and 17 show the temperature distribution along the X-axis and Y-axis, respectively.It can be seen that increasing the pulsed actuation frequency will keep the heated area unchanged, but the surface temperature of the dielectric layer will increase accordingly.Because the pulsed actuation frequency is increased,the number of discharges per unit time increases.However,the single pulsed energy remains unchanged, so the energy loaded per unit time increases and the surface temperature of the dielectric layer increases.According to the results given in figure 18, as the pulsed actuation frequency increases, the temperature rises faster.

    The pulsed voltage peak remains unchanged to change the actuation frequency, but the single pulsed energy of the discharge does not change.The width of the discharge filaments is unchanged, but the brightness is brighter.This is because the increase of frequency makes the discharge filament that can be recorded increase,and the brightness is more obvious.Increasing the frequency increases the input power and the energy injected into the flow field, strengthening its effect on the surrounding flow field and increasing the thermal effect accordingly.

    Figure 21.Variation process of the flow field induced by nanosecond pulsed actuation with different covered electrode widths: (a) 5 mm,(b) 10 mm, (c) 20 mm.

    3.3.Influence of the covered electrode width

    The peak value of the pulsed voltage is 12 kV, the pulsed excitation frequency is 2 kHz, and the width of the covered electrode is changed.

    Figure 19 shows the variation of the single pulsed energy with the change of the rising edge time,falling edge time,and pulsed width with different coverage electrode widths.It can be seen that the width of the covered electrode has a greater influence on the single pulsed energy,and the wider the width of the covered electrode,the greater the single pulsed energy.Under different coverage electrode widths,the changing trend of single pulsed energy with a rising edge, falling edge, and pulsed width is not obvious.However, the smaller the electrode width, the smaller the influence of the electrical parameters on the single pulsed energy change.

    Figure 22.Surface temperature distribution of the actuator after continuous discharge of 120 s (t = 120 s) with different covered electrode widths.

    Figure 20 shows the discharge images of different electrode widths when the pulsed voltage peak is 12 kV, the pulsed frequency is 2 kHz, the rising edge time is 150 ns,the pulsed width is 500 ns,and the falling edge time is 150 ns.It can be seen that the electrode width also has a greater influence on the discharge area and discharge pattern.When the electrode width is 5 mm,the discharge area is smaller and the brightness is relatively dark.As the electrode width increases,the discharge filaments gradually grow, the discharge area gradually increases, and the discharge becomes more intense and brighter.It is shown that the width of the covered electrode is a key parameter that restricts the size of the discharge area,and the length of the discharge filament will not exceed the width of the lower electrode.In order to obtain a larger plasma discharge area,the width of the covered electrode can be appropriately increased.

    Figure 21 shows the change process of the flow field induced by nanosecond pulsed actuation under different covered electrode widths.The results show that the greater the electrode width,the greater the discharge intensity,the greater the power consumption, and the stronger the disturbance of the plasma actuation to the flow field.

    The pulsed voltage peak value is 12 kV, the pulsed excitation frequency is 2 kHz, the rising edge time is 150 ns,the pulsed width is 100 ns,and the falling edge time is 150 ns.The width of the covered electrode is changed.

    Figure 23.The surface temperature changes under different covered electrode widths along the X-axis at Y = 4 mm direction(t = 120 s).

    Figure 24.Changes of surface temperature along the Y-axis at X = 30 mm at different covered electrode widths (t = 120 s).

    Figure 22 shows the surface temperature distribution of different covered electrode widths when the discharge time is 120 s.It can be seen that when the actuation parameters remain unchanged, increasing the electrode width can increase the surface temperature of the dielectric layer.The results in figures 23–25 also confirm this conclusion.The width of the covered electrode limits the extension of the discharge filament so that the length of the discharge filament does not exceed the width of the covered electrode.Therefore,under the actuation parameters in this experiment, increasing the width of the covered electrode will increase the injection energy,and the energy converted into heating the wall surface will also increase,resulting in a higher surface temperature of the dielectric layer.

    When the voltage and frequency are constant, as the covered electrode width is wider, the discharge filament is longer,and the single pulsed energy is higher.The rising and falling edges have little effect on the single pulsed energy,and the single pulsed energy is the highest when the pulsed width is 500 ns.The longer the covered electrode width, the higher the single pulsed energy, the greater the range and intensity of the flow field and temperature field after discharge.

    Figure 25.Variation of surface temperature with time at point A(X = 30 mm, Y = 4 mm) at different covered electrode widths.

    4.Conclusion

    NS-DBD plasma actuation is an active flow control method commonly used to control flow separation.Many characteristics and control mechanisms of NS-DBD plasma actuation,including discharge characteristics and induced flow field characteristics, have been deeply studied.For the high-frequency NS-DBD, it is found that its thermal effect is more obvious, and it has a good effect in the field of plasma antiicing.This paper tests the plasma discharge characteristics of NS-DBD under high frequency conditions, which provides a reference for anti-icing research.

    According to the influence of pulsed voltage peak,pulsed actuation frequency, and covered electrode width on NSDBD discharge characteristics, induced flow field characteristics, and infrared radiation thermal characteristics, the rules of surface DBD on nanosecond pulsed are analyzed.The conclusion is as follows: there exists NS-DBD in the process of pulsed voltage rising and falling, and the air is heated rapidly and the shock wave and jet with lower velocity are induced.The direction of the induced jet is greatly affected by various parameters, and the nanosecond pulsed actuation is not only simple rapid heating but also has the effect of volume force.The covered electrode width will affect the potential distribution during the discharge process, and then affect the extension process of the plasma discharge filament,and the discharge area will be restricted by the electrode width.Increasing the peak value of the pulsed voltage can significantly increase the actuation intensity, and the plasma discharge area is significantly increased.Increasing the pulsed actuation frequency can increase the number of discharges per unit time.Both will cause more energy input, and the temperature rising effect will be more obvious.The surface temperature rising of the dielectric layer is closely related to the amount of input energy.The more the input energy, the more obvious the surface temperature rising effect of the dielectric layer.The surface temperature distribution is directly affected by different discharge patterns, and the temperature rising is more obvious in the area with strong discharge and more concentrated filaments.

    Acknowledgments

    This paper was supported by the National Key R&D Program of China (No.2019YFA0405300) and National Natural Science Foundation of China (Nos.51907205 and 12002363).

    猜你喜歡
    陳杰理科
    Dynamics of magnetic microbubble transport in blood vessels
    陳杰
    牡丹(2023年7期)2023-04-13 13:44:38
    封面攝影
    牡丹(2022年9期)2022-05-17 17:32:03
    和理科男談戀愛也太“有趣”啦
    意林(2021年21期)2021-11-26 20:27:37
    陳杰
    《高師理科學刊》征稿簡則
    文科不懂理科的傷悲
    放風箏
    2017年天津卷理科第19題的多種解法
    紀實
    河南電力(2017年7期)2017-11-30 07:09:16
    男人舔女人下体高潮全视频| 日韩欧美三级三区| av视频在线观看入口| 别揉我奶头~嗯~啊~动态视频| 欧洲精品卡2卡3卡4卡5卡区| 亚洲专区字幕在线| 无遮挡黄片免费观看| 99精品欧美一区二区三区四区| cao死你这个sao货| 午夜福利成人在线免费观看| 国产精品一区二区三区四区免费观看 | 亚洲性夜色夜夜综合| 欧美成狂野欧美在线观看| 欧洲精品卡2卡3卡4卡5卡区| 久热爱精品视频在线9| 巨乳人妻的诱惑在线观看| 国产高清视频在线观看网站| 精品欧美国产一区二区三| 日韩欧美在线二视频| 欧美 亚洲 国产 日韩一| 午夜福利在线在线| 免费看美女性在线毛片视频| xxxwww97欧美| 国产伦一二天堂av在线观看| 日韩中字成人| 最近中文字幕高清免费大全6| 久久精品国产亚洲网站| 国产成人一区二区在线| 哪个播放器可以免费观看大片| 国产视频首页在线观看| 国产精品一二三区在线看| 赤兔流量卡办理| 成人毛片a级毛片在线播放| 99热只有精品国产| 午夜爱爱视频在线播放| 91精品一卡2卡3卡4卡| 亚洲精品自拍成人| 在线观看美女被高潮喷水网站| 桃色一区二区三区在线观看| 不卡视频在线观看欧美| 91aial.com中文字幕在线观看| 99九九线精品视频在线观看视频| 91精品国产九色| 久久久久久久久久成人| 亚洲精品成人久久久久久| 日本爱情动作片www.在线观看| 欧美3d第一页| 三级经典国产精品| 菩萨蛮人人尽说江南好唐韦庄 | 日本欧美国产在线视频| 看黄色毛片网站| 91午夜精品亚洲一区二区三区| 国产一区二区在线av高清观看| 好男人在线观看高清免费视频| 久久人妻av系列| 亚洲人成网站在线播放欧美日韩| 精品午夜福利在线看| 日本欧美国产在线视频| 97超碰精品成人国产| 此物有八面人人有两片| 69人妻影院| 特大巨黑吊av在线直播| 美女脱内裤让男人舔精品视频 | 青春草国产在线视频 | 麻豆国产97在线/欧美| 亚洲国产精品成人久久小说 | 亚洲七黄色美女视频| 人妻久久中文字幕网| 九草在线视频观看| 色吧在线观看| 国产精品美女特级片免费视频播放器| 欧美激情国产日韩精品一区| 美女国产视频在线观看| 久99久视频精品免费| 午夜视频国产福利| 国产真实伦视频高清在线观看| 中国国产av一级| 亚洲一区高清亚洲精品| 久久久久久久久久久免费av| 亚洲精品色激情综合| 日本欧美国产在线视频| 91精品一卡2卡3卡4卡| 91aial.com中文字幕在线观看| 69人妻影院| 国产精品久久电影中文字幕| 欧美区成人在线视频| 久久久久久伊人网av| 99riav亚洲国产免费| 精品久久久久久久久久免费视频| 精品人妻偷拍中文字幕| 观看美女的网站| 国产v大片淫在线免费观看| 大香蕉久久网| 国产高清视频在线观看网站| 亚洲人与动物交配视频| 亚洲性久久影院| 亚洲成人中文字幕在线播放| 美女cb高潮喷水在线观看| 国产成人精品一,二区 | 午夜久久久久精精品| 日本-黄色视频高清免费观看| 国产成人精品婷婷| 女人被狂操c到高潮| 人妻制服诱惑在线中文字幕| 亚洲一级一片aⅴ在线观看| 好男人在线观看高清免费视频| 天堂av国产一区二区熟女人妻| 亚洲精品日韩在线中文字幕 | 九九在线视频观看精品| 91久久精品国产一区二区三区| 美女高潮的动态| 国内精品一区二区在线观看| av免费在线看不卡| 少妇熟女欧美另类| 三级经典国产精品| 亚洲国产欧美人成| 校园春色视频在线观看| 三级国产精品欧美在线观看| 观看美女的网站| 高清毛片免费看| 亚洲第一区二区三区不卡| 精品久久国产蜜桃| 欧美极品一区二区三区四区| 人妻夜夜爽99麻豆av| 淫秽高清视频在线观看| 欧洲精品卡2卡3卡4卡5卡区| 亚洲中文字幕一区二区三区有码在线看| 岛国在线免费视频观看| 久久草成人影院| 国语自产精品视频在线第100页| 久久久国产成人免费| 午夜久久久久精精品| 边亲边吃奶的免费视频| av视频在线观看入口| 又爽又黄a免费视频| 三级国产精品欧美在线观看| 久久久精品欧美日韩精品| 极品教师在线视频| 女同久久另类99精品国产91| 日韩强制内射视频| av黄色大香蕉| 久久99蜜桃精品久久| 男女视频在线观看网站免费| a级一级毛片免费在线观看| 亚洲自拍偷在线| 91久久精品国产一区二区成人| 国产亚洲av片在线观看秒播厂 | 久久久久久久久中文| 亚洲国产精品sss在线观看| 好男人视频免费观看在线| 国产成人a区在线观看| 日韩av不卡免费在线播放| 日本爱情动作片www.在线观看| 国产在视频线在精品| 国产精品综合久久久久久久免费| 欧美一区二区亚洲| 亚洲av成人精品一区久久| 国产乱人视频| 亚洲婷婷狠狠爱综合网| 国产成人午夜福利电影在线观看| 久久久久久久久大av| 简卡轻食公司| 九草在线视频观看| 成人特级av手机在线观看| 欧美成人免费av一区二区三区| 美女脱内裤让男人舔精品视频 | 黄色欧美视频在线观看| 看非洲黑人一级黄片| 午夜福利视频1000在线观看| 久久久久免费精品人妻一区二区| 五月玫瑰六月丁香| 中出人妻视频一区二区| 亚洲国产高清在线一区二区三| 日韩强制内射视频| 亚洲精品乱码久久久v下载方式| 人妻夜夜爽99麻豆av| 国产黄色小视频在线观看| 日韩欧美国产在线观看| 亚洲欧美日韩无卡精品| 国产片特级美女逼逼视频| 亚洲无线观看免费| 国产精品,欧美在线| 最近最新中文字幕大全电影3| 国产精品综合久久久久久久免费| 欧美xxxx黑人xx丫x性爽| 国产精品野战在线观看| 久久久久久伊人网av| 狂野欧美白嫩少妇大欣赏| 黄色一级大片看看| 麻豆国产97在线/欧美| 成熟少妇高潮喷水视频| 2022亚洲国产成人精品| 亚洲av熟女| 成人综合一区亚洲| 国产伦精品一区二区三区视频9| 午夜爱爱视频在线播放| 欧美区成人在线视频| 舔av片在线| 精品久久久久久久久久久久久| 日韩一本色道免费dvd| 国产高清有码在线观看视频| 中文字幕av在线有码专区| 国产高清激情床上av| 国产黄色视频一区二区在线观看 | 久久婷婷人人爽人人干人人爱| 国产老妇伦熟女老妇高清| 日韩,欧美,国产一区二区三区 | 人妻少妇偷人精品九色| 色综合色国产| 超碰av人人做人人爽久久| 国产大屁股一区二区在线视频| 国产av在哪里看| 亚洲精品久久国产高清桃花| 日本熟妇午夜| 91久久精品国产一区二区三区| 国产av一区在线观看免费| 亚洲精品亚洲一区二区| 一本一本综合久久| 成人亚洲欧美一区二区av| 中国国产av一级| 最近最新中文字幕大全电影3| 黄片无遮挡物在线观看| 一进一出抽搐动态| 国产精品1区2区在线观看.| a级毛片a级免费在线| 在线播放国产精品三级| 国产午夜精品久久久久久一区二区三区| 国产成人91sexporn| 少妇人妻精品综合一区二区 | 一区二区三区四区激情视频 | 狠狠狠狠99中文字幕| 日韩大尺度精品在线看网址| 男女做爰动态图高潮gif福利片| 午夜激情福利司机影院| 国产老妇女一区| 一级毛片久久久久久久久女| 天堂av国产一区二区熟女人妻| 国产蜜桃级精品一区二区三区| 午夜激情福利司机影院| 国产高潮美女av| 性插视频无遮挡在线免费观看| 免费av观看视频| 国产成人91sexporn| 啦啦啦啦在线视频资源| 亚洲精品亚洲一区二区| 国产真实乱freesex| h日本视频在线播放| 中文字幕熟女人妻在线| 成人欧美大片| 精品午夜福利在线看| 蜜桃亚洲精品一区二区三区| 人人妻人人澡欧美一区二区| 三级毛片av免费| 插逼视频在线观看| 校园人妻丝袜中文字幕| 亚洲成av人片在线播放无| 在线免费观看的www视频| 晚上一个人看的免费电影| 亚洲精品乱码久久久久久按摩| 一边摸一边抽搐一进一小说| 国产黄片美女视频| 亚洲av熟女| 国产成人午夜福利电影在线观看| 国产精品三级大全| 久久久精品大字幕| 国产私拍福利视频在线观看| 国内久久婷婷六月综合欲色啪| 自拍偷自拍亚洲精品老妇| 亚洲欧美日韩高清专用| 波野结衣二区三区在线| 男女视频在线观看网站免费| 大香蕉久久网| 日日干狠狠操夜夜爽| 欧美变态另类bdsm刘玥| 成年女人永久免费观看视频| 国产综合懂色| 你懂的网址亚洲精品在线观看 | 中出人妻视频一区二区| 只有这里有精品99| 亚洲成人久久爱视频| 美女脱内裤让男人舔精品视频 | 午夜视频国产福利| 五月玫瑰六月丁香| 国产高清三级在线| 在线免费观看的www视频| 国产91av在线免费观看| 欧美日韩精品成人综合77777| 99热6这里只有精品| 亚洲一区二区三区色噜噜| 精品久久久久久久末码| 白带黄色成豆腐渣| 岛国在线免费视频观看| 欧美高清成人免费视频www| 99热只有精品国产| 别揉我奶头 嗯啊视频| 国产真实伦视频高清在线观看| 欧美xxxx性猛交bbbb| 日本-黄色视频高清免费观看| 天天躁日日操中文字幕| 免费无遮挡裸体视频| 我要搜黄色片| 搞女人的毛片| 欧美日韩乱码在线| 免费av观看视频| 色哟哟·www| 国产色爽女视频免费观看| avwww免费| 干丝袜人妻中文字幕| 亚洲激情五月婷婷啪啪| 国产成人aa在线观看| 不卡一级毛片| 精品欧美国产一区二区三| 欧美区成人在线视频| 免费搜索国产男女视频| 伊人久久精品亚洲午夜| 两个人视频免费观看高清| 中国美白少妇内射xxxbb| 久久久精品94久久精品| 人人妻人人看人人澡| 中文精品一卡2卡3卡4更新| 一本精品99久久精品77| 色吧在线观看| 亚洲在线自拍视频| 亚州av有码| 在线播放国产精品三级| 91麻豆精品激情在线观看国产| 亚洲不卡免费看| 亚洲精品成人久久久久久| 干丝袜人妻中文字幕| 人人妻人人澡欧美一区二区| 欧美3d第一页| 五月伊人婷婷丁香| 久久草成人影院| 中文欧美无线码| 99热全是精品| 男的添女的下面高潮视频| 一夜夜www| 国产精品野战在线观看| 丰满乱子伦码专区| 中文在线观看免费www的网站| 又粗又爽又猛毛片免费看| 国产老妇伦熟女老妇高清| 亚洲中文字幕一区二区三区有码在线看| 观看美女的网站| 中文资源天堂在线| 在线观看美女被高潮喷水网站| 国产伦理片在线播放av一区 | 可以在线观看毛片的网站| 亚洲18禁久久av| 男女那种视频在线观看| 亚洲国产欧洲综合997久久,| 美女被艹到高潮喷水动态| 床上黄色一级片| 免费人成在线观看视频色| 久久精品影院6| 日韩大尺度精品在线看网址| 亚洲欧洲日产国产| 国产精品不卡视频一区二区| 日韩成人av中文字幕在线观看| 尾随美女入室| 嫩草影院新地址| 欧美激情国产日韩精品一区| 99国产精品一区二区蜜桃av| 亚洲av二区三区四区| 亚洲精品影视一区二区三区av| 99久久九九国产精品国产免费| 一区二区三区四区激情视频 | 老司机影院成人| 日本-黄色视频高清免费观看| 国产亚洲av嫩草精品影院| 全区人妻精品视频| 国产淫片久久久久久久久| 一进一出抽搐gif免费好疼| av女优亚洲男人天堂| 成人毛片60女人毛片免费| 老司机影院成人| 国产女主播在线喷水免费视频网站 | 色综合亚洲欧美另类图片| 在线播放国产精品三级| 国产精品国产三级国产av玫瑰| 中文在线观看免费www的网站| 国产精品免费一区二区三区在线| 欧美一级a爱片免费观看看| 亚洲熟妇中文字幕五十中出| 午夜视频国产福利| 少妇裸体淫交视频免费看高清| 两性午夜刺激爽爽歪歪视频在线观看| 亚洲aⅴ乱码一区二区在线播放| 国产精品国产高清国产av| 性插视频无遮挡在线免费观看| 夫妻性生交免费视频一级片| 国产亚洲av嫩草精品影院| 三级国产精品欧美在线观看| 热99re8久久精品国产| 欧美又色又爽又黄视频| 韩国av在线不卡| 国产午夜精品久久久久久一区二区三区| 国产精品美女特级片免费视频播放器| 蜜臀久久99精品久久宅男| 色视频www国产| 国产真实伦视频高清在线观看| 欧美成人免费av一区二区三区| 成人国产麻豆网| 天堂影院成人在线观看| 少妇的逼好多水| 国产激情偷乱视频一区二区| 中国国产av一级| 不卡视频在线观看欧美| 国模一区二区三区四区视频| 我的老师免费观看完整版| 一进一出抽搐动态| 深夜a级毛片| 国产91av在线免费观看| 亚洲一级一片aⅴ在线观看| 国产色婷婷99| 乱人视频在线观看| 黄色一级大片看看| 禁无遮挡网站| 亚洲激情五月婷婷啪啪| 91精品国产九色| 男女做爰动态图高潮gif福利片| 乱码一卡2卡4卡精品| 精品午夜福利在线看| 国产伦精品一区二区三区四那| 亚洲最大成人av| 欧美人与善性xxx| 色播亚洲综合网| 女同久久另类99精品国产91| 成人av在线播放网站| 国产精品久久久久久亚洲av鲁大| 久久草成人影院| 男女做爰动态图高潮gif福利片| 国产 一区 欧美 日韩| 九色成人免费人妻av| 中文字幕av在线有码专区| 少妇裸体淫交视频免费看高清| 婷婷亚洲欧美| 久久国内精品自在自线图片| 国产三级中文精品| 精品不卡国产一区二区三区| 性色avwww在线观看| 成人毛片60女人毛片免费| 久久人人爽人人片av| 亚洲欧美日韩无卡精品| 高清在线视频一区二区三区 | 嫩草影院新地址| 波野结衣二区三区在线| 国产真实乱freesex| 免费看a级黄色片| 一卡2卡三卡四卡精品乱码亚洲| 国产成人91sexporn| 免费观看精品视频网站| 一级毛片久久久久久久久女| 在线观看66精品国产| 国产老妇伦熟女老妇高清| av天堂在线播放| 极品教师在线视频| 精品久久久久久成人av| 亚洲av免费在线观看| 小蜜桃在线观看免费完整版高清| 日日干狠狠操夜夜爽| 九九在线视频观看精品| 亚洲人与动物交配视频| 国产一区二区在线观看日韩| 亚洲在线观看片| 亚洲无线在线观看| 国产精品.久久久| 毛片女人毛片| 免费看a级黄色片| 国产精品伦人一区二区| 久久精品国产清高在天天线| 欧美变态另类bdsm刘玥| 卡戴珊不雅视频在线播放| 永久网站在线| 美女高潮的动态| 日韩制服骚丝袜av| 99热精品在线国产| 观看美女的网站| 国产老妇女一区| 亚洲人与动物交配视频| 日韩精品有码人妻一区| 日韩大尺度精品在线看网址| 日本免费a在线| 九草在线视频观看| 精品欧美国产一区二区三| 2021天堂中文幕一二区在线观| 亚洲图色成人| 久久久久久久午夜电影| 国产精品一区二区三区四区免费观看| 久久午夜亚洲精品久久| 国产三级在线视频| 一本精品99久久精品77| or卡值多少钱| 一级黄色大片毛片| 久久人人爽人人片av| 亚洲欧洲国产日韩| 国产欧美日韩精品一区二区| 日日干狠狠操夜夜爽| 天堂中文最新版在线下载 | 欧美精品一区二区大全| 99久国产av精品| 亚洲自偷自拍三级| 内射极品少妇av片p| 一本久久精品| 久久99精品国语久久久| 一本精品99久久精品77| 国产精品精品国产色婷婷| 网址你懂的国产日韩在线| 国产淫片久久久久久久久| 日本色播在线视频| 国内精品久久久久精免费| 日本免费a在线| 欧美日韩一区二区视频在线观看视频在线 | 美女被艹到高潮喷水动态| 91在线精品国自产拍蜜月| 国产成人精品久久久久久| 色视频www国产| 色吧在线观看| 成人永久免费在线观看视频| 久久久久九九精品影院| 免费看a级黄色片| 中文字幕人妻熟人妻熟丝袜美| 成人特级av手机在线观看| h日本视频在线播放| 亚洲人与动物交配视频| 九九久久精品国产亚洲av麻豆| 亚洲丝袜综合中文字幕| 日日摸夜夜添夜夜爱| 亚洲成人久久性| 精品日产1卡2卡| 又粗又硬又长又爽又黄的视频 | 我要看日韩黄色一级片| 99视频精品全部免费 在线| 秋霞在线观看毛片| av免费在线看不卡| 麻豆成人av视频| 美女cb高潮喷水在线观看| 中文欧美无线码| 亚洲人成网站高清观看| 国产黄片视频在线免费观看| 三级国产精品欧美在线观看| 午夜福利成人在线免费观看| 色哟哟·www| 国语自产精品视频在线第100页| 国产麻豆成人av免费视频| 九九久久精品国产亚洲av麻豆| 国产成人午夜福利电影在线观看| 能在线免费看毛片的网站| 国产亚洲5aaaaa淫片| 最近最新中文字幕大全电影3| 超碰av人人做人人爽久久| 国产一级毛片七仙女欲春2| 如何舔出高潮| 成年女人永久免费观看视频| 在线观看66精品国产| 国产视频内射| 国产精品国产三级国产av玫瑰| 特级一级黄色大片| 亚洲国产精品成人久久小说 | 美女被艹到高潮喷水动态| 欧美成人a在线观看| 99国产极品粉嫩在线观看| 国产又黄又爽又无遮挡在线| 天美传媒精品一区二区| 一区二区三区高清视频在线| 九九爱精品视频在线观看| 国产视频首页在线观看| 夜夜看夜夜爽夜夜摸| 嘟嘟电影网在线观看| 国产一区二区在线av高清观看| 婷婷精品国产亚洲av| 午夜激情福利司机影院| 美女黄网站色视频| 久久久久国产网址| 三级毛片av免费| 国产精品免费一区二区三区在线| 两性午夜刺激爽爽歪歪视频在线观看| 成人av在线播放网站| 精品一区二区免费观看| 国产成人一区二区在线| 国产淫片久久久久久久久| 国产伦精品一区二区三区视频9| 欧美日本视频| 久久国内精品自在自线图片| 日韩中字成人| 18+在线观看网站| 给我免费播放毛片高清在线观看| 久久精品国产自在天天线| 小蜜桃在线观看免费完整版高清| 一本久久精品| 3wmmmm亚洲av在线观看| 在线国产一区二区在线| 亚洲最大成人手机在线| 成人国产麻豆网| 亚洲国产日韩欧美精品在线观看| 一本一本综合久久| 国产伦精品一区二区三区视频9| 久久久久免费精品人妻一区二区| 国产精品嫩草影院av在线观看| avwww免费| 日日摸夜夜添夜夜添av毛片| 日韩亚洲欧美综合| 国产高清激情床上av| 中文资源天堂在线| 国产淫片久久久久久久久| 国产精品久久久久久精品电影| 两性午夜刺激爽爽歪歪视频在线观看| 熟女人妻精品中文字幕| 我的女老师完整版在线观看| 亚洲欧美日韩东京热| 美女 人体艺术 gogo| 日韩三级伦理在线观看| 国产黄片美女视频| 婷婷色av中文字幕| 亚洲国产欧美人成|