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

    Burst pulses for positive corona discharges in atmospheric air:the collective movement of charged species

    2021-06-21 02:00:36YongkangPENG彭永康XiaoyueCHEN陳小月LeiLAN藍(lán)磊HaoyuZHAN詹皓予YongcongLIU劉永聰WanglingHE何旺齡andXishanWEN文習(xí)山
    Plasma Science and Technology 2021年6期
    關(guān)鍵詞:永康

    Yongkang PENG(彭永康),Xiaoyue CHEN(陳小月),Lei LAN(藍(lán)磊),Haoyu ZHAN(詹皓予),Yongcong LIU(劉永聰),Wangling HE(何旺齡) and Xishan WEN(文習(xí)山)

    1 School of Electrical Engineering and Automation,Wuhan University,Wuhan 430074,People’s Republic of China

    2 School of Electrical and Electronic Engineering,North China Electric Power University,Baoding 071000,People’s Republic of China

    Abstract Positive corona burst pulses are an unstable pulse mode.They appear in a small range of the onset stage,and their current pulses result from the collective movement of charged species.This paper focused on the connections between these pulses and the collective movement of charged species.The movement of species is divided into four parts with respect to time:the(1)initial growth of species,(2)formation and development of the streamer region and negative ion sheath,(3)dead time(the time interval between the pulses),and(4)rapid re-growth of species.The movement of the species in the four parts and the correspondence with the current pulse were analyzed.The numerical results indicated the following:the rapid rising of the species matched the rising edge of the pulses,the streamer region,and negative ion sheath appeared in the falling edge of the primary pulse,and the rapid re-growth of species matched the re-ignition of the pulses.The results were in qualitative agreement with deductions and experimental observations in the literature.

    Keywords:positive corona discharge,burst pulses,pre-onset pulses,collective movement of charged species

    1.Introduction

    Positive corona discharge is widely known to involve four stages[1–3],where these stages are marked by burst pulses[4],streamer pulses[5],glow coronas[6],and breakdown pulses[7].The burst pulses appear in a small range of the onset stage,where the pulse mode is a large primary pulse followed by several secondary pulses.According to the venerable body of research on the connections between the burst pulses and the collective movement of charged species,the current pulses result from the collective movement of charged species.

    Trichel[8],who first detailed the burst pulses of positive coronas in 1939,proposed a mechanism involving electron avalanche progress in the initial pulses and the role of negative ions in continuous glow[8].Morrow[9]stated that‘Trichel deduced most of the details of the mechanism proposed for positive glow corona’[10].According to Amin’s observation[4]regarding burst pulses,a pre-onset positive corona has three critical potentials for different pulse states,and the new burst pulse starts when the space charge has crossed 80%of the gap.Miyoshi,Hosokawa[3]and Trinh[2]also observed the typical pulse train.Miyoshi,Hosokawa discussed the relationships between the burst pulse and streamer pulse,and they proposed that the successor electron avalanches supply the initiating electrons to ignite the secondary pulse(figure 2 of[3]).Trinh synopsized the entire discharge progress of positive coronas,presenting the typical burst pulses in figure 7 of[2].

    In recent decades,simulation methods[11,12]have been used to analyze corona discharge together with experimental ones.The research of Morrow and his coworkers is canonical in the related field.They reported a series of theories[11,13–15],methods[16,17],and classical documentation[18,19]with a fluid model for positive and negative coronas in air and some other gases.Two previous studies[15,20]formulated a theory for positive coronas,which describes the development of the first primary pulse and the first secondary pulse when a positive voltage is abruptly applied to a point in the air at atmospheric pressure.The initial changes in the pulse,distribution of charged species,and electric field were studied.Notably,the propagation of positive ions and electrons in the streamer was elucidated.

    Liu et al[21–24]advanced Morrow’s research by simulating the transformation process.Their findings were consistent with those of Morrow in terms of positive corona simulation.Many other researchers analyzed the pulses of positive corona or burst pulses[25–28].He et al[25]studied burst pulses in the positive half of alternating current corona.Li et al[26]studied the detailed characteristics of positive corona.Liu et al simulated pulses of positive corona in[27].He et al also proposed a theoretical model based on the movement of the space charged clouds and explained some physical mechanism of the pulses of positive corona with the model[28].More recently,Chen et al[29]studied the mode transition of negative coronas from Trichel pulses to pulseless glow.He emphasized the great importance of gas heating.The gas dynamics and a detachment reaction are taken into account in the model.

    Notably,in[15],the charged species were plotted on a single graph(figure 13 in[14]),which indicated an obvious region where positive ions have the same order of magnitude as electrons;this region corresponded to the so-called streamer region.A later study also indicated the same phenomenon(see figure 6.10(a)in[30]).

    Building on these previous studies,the present study aimed to investigate(i)the collective movement of the charged species and(ii)the charged species’connections with the pulse current in the stage marked by a positive corona burst pulse.

    2.Model description

    A one dimensional model with a DC 30 kV supplied voltage used in the present paper had been calculated in our former work[31].An in-depth study on the basis of its simulation results was done in the present paper.The focus in this paper is the coupled relationship between the collective movement of the charged species and the electric field.Some important details for the model are recapitulate here.

    The continuity equations for electrons,positive ions,negative ions,are as follows:

    wherene,np,nn,μe,μp,μn,De,Dp,Dn,are the number densities(in m?3),mobilities(in m2(V?1s?1)),and diffusion coefficients(in m2s?1)of electrons,positive ions and negative ions,respectively.Sph,S0are the photoionization source term and initial source term.Ionization,attachment of electrons to neutral molecules,recombination of electrons with positive ions and recombination of positive and negative ions,are with the reaction coefficientsα(m-1),η(m-1)andβepandβnprespectively.

    To make the electric source term more richly descriptive,the role of metastable molecules was also considered in our model.The consideration of such a role was used by Lowke[22]to account for the effect of electron detachment from negative ions

    αmis the creation coefficient for metastable molecules impact with electrons,kdis the detachment coefficient for detaching electrons from negative ions,andkqis quenching rate constant for metastable molecules.

    The electric fieldand electric potentialφcan be calculated using Poisson’s equation

    The parameters used in the simulation are shown in table 1.

    Table 1.The transport parameters for the charged species.

    The photoionization source termSphsolved by the Helmholtz equations[32]:

    The values ofAjandλjare shown in table 2.Po2is the partial pressure of O2in Pa andis the emission of ionizing UV photons.

    Table 2.Constant coefficients for equation(6).

    The initial number densities of charged particles were given artificially by the Gaussian distribution in equation(8)

    where Nmax=1016m?3,r0=5×10?4m,and s0=1×10?4m.Zero initial density was assumed for negative ions and metastable molecules,as shown in figure 1.

    Figure 1.Initial distribution of species and the electric field in the discharge gap(t=0).

    The electrode system comprised a solid cylinder with a radius of 0.5 mm and a coaxial steel mesh cage with a radius of 200 mm.A 30 kV DC voltage was applied to the anode(the inner electrode)when the cathode(the outer electrode)was grounded.

    E-fieldis the legend of electric field.We put it with the densities of charged particles in the same figure to show the corresponding relationship between charged particle density and electric field.

    The mesh of the model was divided into four parts,as shown in table 3.

    Table 3.Mesh of the model.

    Figure 2 illustrates the flow of the computing program used for the simulation.The finite element method was used to solve these equations.Specifically,when the computing program is initiated,each variable is given an initial value.Subsequently,the electric field is calculated by solving Poisson’s equation,andSphis obtained by solving the Helmholtz equations.Thereafter,the densities of electrons,positive ions,negative ions,and metastable molecules—denoted asne,Np,Nn,andnm,respectively—are obtained by solving the continuity equations.The requisite reaction coefficients take those values noted in[31].Error checking was conducted to verify the convergence characteristics.To investigate the successive burst pulse trains,tmaxwas set to 200 μs.In addition,the time step,Δt,depends on the absolute tolerance;Δt had the orders of 10?9–10?12s in a single pulse stage and 10?8s in the interval stage.

    Figure 2.Flowchart of the computing program.

    3.Simulation results and discussion

    The obtained positive corona burst pulses at 30 kV are shown in figure 3.The primary pulse and secondary pulses had an alternating pattern;this finding is in qualitative agreement with the experimental observations in[2,3]and numerical results of[20].

    The collective movement of charged species was analyzed.This collective movement was divided into four parts with respect to time:the(1)initial growth of the species;(2)formation and development of the streamer region and negative ion sheath;(3)dead time(the time interval between the pulses);and(4)rapid re-growth of the species.These parts match the different stages of the pulses(figure 4).

    Figure 3.Positive corona burst pulse.

    Figure 4.Time divisions of burst pulse.

    3.1.Initial growth of species

    From t=0,the charged species began to grow all over the discharged channel.The positive ions and electrons around the anode both rapidly increased from 1016m?3order to 1019m?3order.The densities of the positive ions,negative ions,and electrons then reached their peak values of 1.54×1019m?3,4.458×1017m?3,and 1.48×1019m?3at 0.445 ns,0.725 ns,and 0.915 ns,respectively.Their distributions and corresponding positions are presented in figure 5.

    Figure 5.Initial growth and peak value of the charged species near the anode.

    Figures 1 and 5(a)illustrate the rising edge of the first pulse from ta=0 to tb=0.445 ns.In this period,there were more positive ions than electrons near the anode.Because the electrons were absorbed by the anode,positive ions attained their peak first,resulting in the trend for the electric field near the anode falling to become a valley.The current increased due to the rapid increase of species and the decrease of the electric field.Subsequently,the low electric field decreased and choked the collisional ionization at the corresponding location,which decreased the current and consequently formed the peak of the current pulse.

    3.2.Formation and development of the streamer region and negative ion sheath

    3.2.1.Formation and development of the streamer region.Although the ionization at the remaining part near the anode was still strong,the electrons and negative ions were able to accumulate to their peaks where the electric field was low,as depicted in figures 5(b)and(c).

    From t=0.725 ns,the density of electrons near the anode reached the same order of magnitude as positive ions did,thus forming a streamer region[6,15,33].As illustrated in figure 6,the streamer region was gradually developing toward the direction of the cathode,and the densities of electrons and positive ions in the streamer region decreased gradually.In the streamer region,the density of positive ions was slightly greater than that of electrons(enlarged part of figure 6).This resulted in the electric field at the corresponding location decreasing to approximately 103V m?1.In other words,the electric field in the streamer region was low—a similar outcome as that from extending the radius of the anode.

    Figure 6.Densities of particles and electric field in the discharge gap at t=20 ns.

    In the present study,we propose a criterion to describe the little difference between positive ions and electrons in the streamer region.The length of the streamer region was included in the analysis if the corresponding density difference between positive ions and electrons was <10%of the density of positive ions in the streamer region(equation(9))

    Figure 7 marks the boundary of the streamer region and its direction of movement with a black dotted line.The r coordinate and the corresponding electric field are labeled in parentheses.The streamer region was 0.011 mm at t=0.725 ns,and it was developing toward the cathode until 60 ns.Subsequently,the direction of movement of the streamer region turned back toward the anode.At approximately t=8 μs,the boundary of the streamer region arrived at the anode,where the streamer region completely disappeared.

    Figure 7.Movement of the streamer region.

    The length of the streamer region versus time is plotted in figure 8.The trend had a steep rising edge and a gentle falling edge.

    Figure 8.Length of the streamer region.

    Compared with the current pulse,the extending stage of the streamer region matched the main falling edge of the current from tb=0.725 ns to t=60 ns.During this period,the current fell to 10%of the peak value.The shrinking stage matched the remaining part of the current falling edge.At t=8 μs,the remaining current of the pulse was almost 0;this indicated that the current on the falling edge of the pulse was mainly formed by the migration of streamers.

    3.2.2.Formation and development of the negative ion sheath.When the moving direction of the streamer region turned back,the number of negative ions began to increase at the end of the streamer region.As depicted in the enlarged drawing of figure 9(a),the negative ions exhibited an obvious peak at the cross of the electrons and negative ions.As shown in figures 9(a)and(b),this small peak grew gradually and reached the same order of magnitude as the positive ions did at the end of the streamer region.The negative ions formed a stable sheath near the streamer region,which was similar to the so-called glow layer discussed by much literature[2,3,9,34].This meant that the negative ions also played an important role in the primary burst pulse.After the formation of the negative ion sheath,the negative ions migrated to and disappeared at the anode after the streamer region.It was shown in figures 9(c)and(d).

    Figure 9.Densities of charged particles and electric field in the discharge gap.

    The migrated track of the negative ion sheath is presented by the shadow between the two curves in figure 10.We define that the negative ion sheath is where the difference between positive ions and negative ions was less than 10% of the density of positive ions.The distance between the left and right edges is the length of the negative ion sheath.When the streamer region began to shrink back to the anode,the length of the negative ion sheath was extending,and moving toward the anode.

    Figure 10.Position of the negative ion sheath.

    3.3.Dead time

    As shown in figure 4,the time from tcto tdis the dead time between the first two pulses,and the time from tfto tgis the dead time between the second and third pulses.During the first dead time,the current was maintained at a low level of<5 μA(15% of the secondary pulse).The dead time is characterized by the migration of positive ions,recovery of the electric field,and tentative growth of electrons.

    The migration of positive ions could be divided into two stages.The first stage did not belong to the dead time,and it occurred after the formation of the streamer region(at t=0.1–10 μs).As indicated by the red dotted line of figure 9,the high-density region of positive ions gradually extended from the end of the streamer region to the cathode in this first stage.The second stage constituted the main part of the dead time,as indicated in figures 9(d)and 11.The peak density of positive ions near the anode migrated to the cathode,which led to the electric field recovering gradually;this migration determines the dead time length.

    With the recovery of the electric field,electrons began to develop tentatively.Initially,density of electrons rises up at random positions in the ionization region of the discharging channel owing to the photoionization and the initial avalanche.The electrons then rapidly developed and formed more electron avalanches,which adhered to the anode and absorbed by the anode.The density of electrons rises up and declines in the channel as the avalanches formed and subsided.The accompanying positive ions stayed and moved to the cathode under the force of the electric field.

    As noted in figure 11,several positive ion waves of different sizes were scattered in the migration area,these waves were the residual effect of the electron avalanches at the same location dozens of ns or several hundred ns before.The scales of these tentative electron avalanches increased as this process repeated itself.

    Figure 11.Tentative electron avalanches in corona dead times before current pulses of different versions.

    During this period,despite the low current,the movement of positive ions and the recovery of the electric field were creating the conditions for the ignition of the next pulse.The tentative electron avalanches were probing if the conditions for ionization were prepared.

    3.4.Rapid re-growth of the species

    When the electric field and species conditions were prepared,the pulses began to emerge.The rapid re-growth of the charged species had two variants.The first had no obvious streamer region and was created by one of the tentative electron avalanches.Different from the conventional rise of tentative avalanches,this mutant rise drives the simultaneous growth of both positive and negative ions.A fragment of its life circle development is displayed in figure 12.

    Figure 12.Development of life cycle of a small-particle pulse.

    The avalanche emerged at td=16.10 μs,as shown in figure 12(a).The negative ions emerged at t=16.16 μs.The species reached their peak values at te=16.2 μs.Subsequently,the densities of electrons and negative ions began to decrease rapidly,and the positive ions slowly drifted toward the cathode.The electrons density reduced to a low level which cannot be seen in the figure 12(g)at tf=16.32 μs,and the negative ions became unobservable after t=19.40 μs.The life of the first secondary pulse ended at t=3.30 μs;the pulse grew in the first 0.1 μs and decayed thereafter.

    The rising edge of the first secondary pulse was from tdto te(figure 4).Similar to the first primary pulse,the rising edge of the first secondary steamer pulse matched the rise of the species.Conversely,owing to the absence of a streamer region,the falling edge of the secondary pulse matched the falling of electrons from teto tf.After the electrons of the first secondary pulse disappeared,the pulse current entered the second dead time at tfto tg.

    Similar to the first primary pulse,the second variant of the re-growth of the species was able to trigger a streamer.It emerged at tg=23 μs(figure 4).Its streamer region at the peak time is illustrated in figure 13.From a comparison of figures 13 with 11(b),it is evident that the densities of electrons and positive ions rapidly increased.At t=23.14 μs,both densities reached a 1019m?3order of magnitude at the same time,and the density of negative ions reached 1017m?3;this behavior of the charged species was identical to the behaviors described in section 3.1.

    Figure 13.Rapid growth of the second primary streamer region at t=23.14 μs.

    Owing to the artificial initial distribution of the species,the first primary pulse skipped over a part of the accumulated progress of the species,doing so while this process had been completed in the second primary pulse from tg=23 μs to t=23.12 μs;this phenomenon was studied in our former research[31].Although the current rise was slight during this period,species growth had already commenced.Such growth is necessary for the last jump of the rising edge of the primary pulse current.

    4.Summary and conclusion

    The collective migration of electrons,positive ions,and negative ions in positive corona discharges in the air was analyzed in this paper.Our findings yielded the following.

    (1)The rapid growth of the species for the first pulse is from tato tb.This growth matches the rising edge of the pulse current.When the positive ions peak at tb,the current peaks at the same time.

    (2)The streamer region formed rapidly at t=0.725 ns in the present study.The region extends to the direction of the cathode from 0.725 ns to 60 ns and shrinks back to the anode from 60 ns to 8 μs.The extending stage of the streamer matches the first 90%of the current fade of the pulse current falling edge in 59.275 ns,and the shrinking stage matches the remaining 10% current drops at approximately 7.4 μs.

    (3)The negative ion sheath migrates to the anode and shrinks to the anode on the heel of the streamer region;it appears in the second stage of the pulse current falling edge.

    (4)In the present model,the current in dead time was <5 μA.The main species activities in dead time are the migration of positive ions and the tentative growth of electrons.The tentative growth of electrons is probed,and thus we can inspect whether the electric field has recovered to a level sufficient for ionization.

    (5)The rapid re-growth of the charged species has two variants.The first has no streamer region,grows in the first 0.1 μs,decays in the following 3.2 μs,and matches the first secondary pulse.The second behaves similarly to the first primary pulse,and it replaces the artificial initial distribution of the species with a longer rising edge.

    Acknowledgments

    This work was supported by National Natural Science Foundation of China(No.51907145).

    猜你喜歡
    永康
    吳永康雕塑作品選登
    名作欣賞(2023年17期)2023-06-28 00:48:08
    論姚永康瓷塑《世紀(jì)娃》的藝脈文心
    中國典型環(huán)境地區(qū)用戶模型搭建及應(yīng)用
    北京汽車(2021年2期)2021-05-07 04:00:20
    試論中國陶藝“博古藏今”之圖式——以姚永康捏塑《世紀(jì)娃》為例證
    陸永康男兒膝下是理想的天空
    中國工人(2019年9期)2019-09-10 07:22:44
    永康日報(bào)
    傳媒評論(2019年3期)2019-06-18 10:59:14
    “店小二”王永康,這個(gè)月可真忙!
    窺視花樣年華,引得春光乍泄
    攝影之友(2018年3期)2018-04-12 01:53:12
    51歲兒子懷抱82歲母親江邊嘮嗑
    晚晴(2016年10期)2016-11-16 22:13:59
    流變與傳承
    ——秦漢時(shí)期“伏日”考論
    古代文明(2013年4期)2013-07-18 12:08:45
    高清不卡的av网站| 成人特级av手机在线观看| 亚洲成人手机| 麻豆乱淫一区二区| 国产成人午夜福利电影在线观看| 男的添女的下面高潮视频| 欧美bdsm另类| 亚洲精品国产色婷婷电影| 欧美精品国产亚洲| av国产精品久久久久影院| videossex国产| 妹子高潮喷水视频| 91狼人影院| freevideosex欧美| 偷拍熟女少妇极品色| 狂野欧美激情性bbbbbb| 亚洲精华国产精华液的使用体验| 婷婷色麻豆天堂久久| 国产亚洲欧美精品永久| 亚洲四区av| 新久久久久国产一级毛片| 三级经典国产精品| 亚洲欧美成人综合另类久久久| 久久人人爽av亚洲精品天堂 | 成人亚洲精品一区在线观看 | 人妻少妇偷人精品九色| 午夜视频国产福利| 亚洲真实伦在线观看| 激情五月婷婷亚洲| 国产视频内射| 亚洲精品一二三| 午夜福利高清视频| 国国产精品蜜臀av免费| 免费在线观看成人毛片| 欧美激情国产日韩精品一区| 少妇丰满av| 免费观看性生交大片5| 青春草国产在线视频| av视频免费观看在线观看| 1000部很黄的大片| 日韩av免费高清视频| 亚洲av免费高清在线观看| 乱码一卡2卡4卡精品| 久久国产精品大桥未久av | 久久国产亚洲av麻豆专区| 香蕉精品网在线| 99re6热这里在线精品视频| 国产精品福利在线免费观看| 一二三四中文在线观看免费高清| 欧美成人午夜免费资源| 熟女电影av网| 国产国拍精品亚洲av在线观看| 2021少妇久久久久久久久久久| 中文字幕人妻熟人妻熟丝袜美| 老司机影院毛片| 九九久久精品国产亚洲av麻豆| 少妇裸体淫交视频免费看高清| 能在线免费看毛片的网站| 最近中文字幕高清免费大全6| av在线蜜桃| 国产白丝娇喘喷水9色精品| 精品国产露脸久久av麻豆| 国产有黄有色有爽视频| 亚洲精品自拍成人| 精品一区二区免费观看| 97在线视频观看| 六月丁香七月| 乱系列少妇在线播放| 亚洲av综合色区一区| 亚洲国产毛片av蜜桃av| 久久久久久久久久成人| 寂寞人妻少妇视频99o| 国产精品无大码| 欧美精品国产亚洲| 亚洲成人手机| 欧美丝袜亚洲另类| 夜夜爽夜夜爽视频| 蜜臀久久99精品久久宅男| 18禁在线播放成人免费| 亚洲人成网站在线播| 亚洲中文av在线| 亚洲欧洲国产日韩| 国内精品宾馆在线| 免费看不卡的av| 免费人妻精品一区二区三区视频| 99久久综合免费| 久久精品国产鲁丝片午夜精品| av天堂中文字幕网| 亚洲一级一片aⅴ在线观看| 国内精品宾馆在线| 99久久精品一区二区三区| 亚洲一区二区三区欧美精品| 国产高清有码在线观看视频| 新久久久久国产一级毛片| 久久人妻熟女aⅴ| 一级片'在线观看视频| 女性被躁到高潮视频| 成人一区二区视频在线观看| 国产精品福利在线免费观看| videossex国产| 亚洲av成人精品一区久久| 1000部很黄的大片| 精品一区在线观看国产| 欧美日韩亚洲高清精品| 亚洲电影在线观看av| 自拍欧美九色日韩亚洲蝌蚪91 | 欧美极品一区二区三区四区| 亚洲第一区二区三区不卡| 欧美精品国产亚洲| 97精品久久久久久久久久精品| 中文资源天堂在线| 久久久成人免费电影| 国内精品宾馆在线| 另类亚洲欧美激情| 亚洲欧洲国产日韩| 免费观看的影片在线观看| 精品视频人人做人人爽| 日韩中字成人| 最近手机中文字幕大全| 日韩成人伦理影院| 99精国产麻豆久久婷婷| 日韩三级伦理在线观看| 人妻少妇偷人精品九色| 大香蕉久久网| 夫妻午夜视频| 国产69精品久久久久777片| www.av在线官网国产| 菩萨蛮人人尽说江南好唐韦庄| 久久国内精品自在自线图片| 国产精品不卡视频一区二区| 又大又黄又爽视频免费| 国产人妻一区二区三区在| 蜜桃在线观看..| 少妇人妻一区二区三区视频| 好男人视频免费观看在线| 久久国产精品男人的天堂亚洲 | 国产爽快片一区二区三区| 黄色配什么色好看| 国产亚洲91精品色在线| 免费看日本二区| 波野结衣二区三区在线| 国产男人的电影天堂91| 久热久热在线精品观看| 国产高潮美女av| 中国国产av一级| 久久久久久久大尺度免费视频| 欧美日韩综合久久久久久| 中文精品一卡2卡3卡4更新| 国产成人aa在线观看| 在线精品无人区一区二区三 | 亚洲av中文字字幕乱码综合| 久久99热这里只有精品18| 成年人午夜在线观看视频| 国产精品国产av在线观看| 少妇人妻一区二区三区视频| 国产亚洲精品久久久com| 午夜免费观看性视频| 国产一区二区三区av在线| 亚洲成人手机| 国产伦精品一区二区三区四那| 午夜福利在线在线| 插阴视频在线观看视频| 人人妻人人爽人人添夜夜欢视频 | 夫妻午夜视频| 直男gayav资源| 交换朋友夫妻互换小说| 亚洲高清免费不卡视频| 欧美变态另类bdsm刘玥| 少妇人妻久久综合中文| 国产乱来视频区| 成人免费观看视频高清| 日韩强制内射视频| 国产精品麻豆人妻色哟哟久久| 久久久亚洲精品成人影院| 午夜免费鲁丝| 久久久精品免费免费高清| 亚洲精品日韩av片在线观看| 欧美精品亚洲一区二区| 人体艺术视频欧美日本| 在线观看美女被高潮喷水网站| 老熟女久久久| 欧美日韩在线观看h| 久久午夜福利片| 日韩一区二区视频免费看| 一区在线观看完整版| 亚洲综合精品二区| 高清视频免费观看一区二区| 高清av免费在线| 国产精品久久久久久精品古装| 九九在线视频观看精品| 日韩成人av中文字幕在线观看| 国产乱人偷精品视频| 在线观看av片永久免费下载| 日韩免费高清中文字幕av| 国产成人午夜福利电影在线观看| 免费不卡的大黄色大毛片视频在线观看| 51国产日韩欧美| 亚洲经典国产精华液单| 男女边吃奶边做爰视频| 国产 一区 欧美 日韩| 只有这里有精品99| 久久久色成人| 高清视频免费观看一区二区| 有码 亚洲区| 国产爽快片一区二区三区| 国产爱豆传媒在线观看| 精品久久久精品久久久| 中国美白少妇内射xxxbb| av免费在线看不卡| 欧美成人精品欧美一级黄| 国产免费福利视频在线观看| 亚洲aⅴ乱码一区二区在线播放| 日韩电影二区| 精品一区二区免费观看| 欧美 日韩 精品 国产| 制服丝袜香蕉在线| av卡一久久| 人妻 亚洲 视频| 欧美bdsm另类| 欧美3d第一页| 国产精品一区二区性色av| 久久99热这里只频精品6学生| 亚洲国产毛片av蜜桃av| 九九久久精品国产亚洲av麻豆| 99热6这里只有精品| 亚洲av二区三区四区| 亚洲精品国产色婷婷电影| 国产日韩欧美亚洲二区| 精品国产乱码久久久久久小说| 国产片特级美女逼逼视频| 日韩国内少妇激情av| 777米奇影视久久| 久久久久久久国产电影| 国产精品福利在线免费观看| 欧美xxxx性猛交bbbb| 国产男女超爽视频在线观看| 啦啦啦中文免费视频观看日本| 久久精品熟女亚洲av麻豆精品| 久久精品久久久久久噜噜老黄| 国产精品99久久99久久久不卡 | 看十八女毛片水多多多| 偷拍熟女少妇极品色| 国产色婷婷99| 亚洲av.av天堂| 国产无遮挡羞羞视频在线观看| 精品熟女少妇av免费看| 免费观看性生交大片5| av专区在线播放| 亚洲伊人久久精品综合| 91精品国产九色| 亚洲av欧美aⅴ国产| 久久97久久精品| 久久久精品免费免费高清| 亚洲av.av天堂| 久热久热在线精品观看| av网站免费在线观看视频| 天天躁日日操中文字幕| 亚洲国产日韩一区二区| 综合色丁香网| 亚洲欧美成人精品一区二区| 久久久精品94久久精品| 国产色婷婷99| 日韩中文字幕视频在线看片 | 十分钟在线观看高清视频www | 国产免费又黄又爽又色| 精品久久久久久久久亚洲| 欧美精品亚洲一区二区| 精品熟女少妇av免费看| 国语对白做爰xxxⅹ性视频网站| 精品少妇黑人巨大在线播放| 偷拍熟女少妇极品色| 久久婷婷青草| 亚洲自偷自拍三级| 极品少妇高潮喷水抽搐| 亚洲美女视频黄频| 国产无遮挡羞羞视频在线观看| 中文天堂在线官网| 麻豆成人午夜福利视频| 成年美女黄网站色视频大全免费 | 我要看日韩黄色一级片| 午夜免费观看性视频| 亚洲一级一片aⅴ在线观看| 国产欧美日韩一区二区三区在线 | 男女免费视频国产| 国产精品久久久久成人av| 成人免费观看视频高清| 一二三四中文在线观看免费高清| 免费大片黄手机在线观看| 丰满乱子伦码专区| 国产欧美日韩精品一区二区| 大片电影免费在线观看免费| 国产乱人偷精品视频| av国产免费在线观看| 午夜激情福利司机影院| av网站免费在线观看视频| 春色校园在线视频观看| 欧美日韩一区二区视频在线观看视频在线| 日韩欧美精品免费久久| 久久久久久久久久人人人人人人| 啦啦啦中文免费视频观看日本| 日韩在线高清观看一区二区三区| 伦理电影免费视频| 久久精品熟女亚洲av麻豆精品| 亚洲一级一片aⅴ在线观看| 欧美高清性xxxxhd video| 只有这里有精品99| 天堂俺去俺来也www色官网| 91久久精品电影网| 黄色视频在线播放观看不卡| 久久久久久九九精品二区国产| 亚洲综合精品二区| 夜夜看夜夜爽夜夜摸| 王馨瑶露胸无遮挡在线观看| 国产黄片美女视频| 18禁在线无遮挡免费观看视频| 亚洲久久久国产精品| 91久久精品国产一区二区三区| 精品视频人人做人人爽| 色哟哟·www| 国产一区有黄有色的免费视频| 亚洲美女视频黄频| 亚洲国产精品国产精品| 亚洲怡红院男人天堂| 欧美极品一区二区三区四区| 婷婷色综合大香蕉| 精品99又大又爽又粗少妇毛片| 少妇人妻一区二区三区视频| 亚洲精品456在线播放app| 国产亚洲5aaaaa淫片| 狂野欧美激情性bbbbbb| 这个男人来自地球电影免费观看 | 毛片女人毛片| 在线看a的网站| 美女cb高潮喷水在线观看| 欧美日韩综合久久久久久| 美女福利国产在线 | 深爱激情五月婷婷| 亚洲无线观看免费| 视频中文字幕在线观看| 99国产精品免费福利视频| 美女xxoo啪啪120秒动态图| av专区在线播放| 一级片'在线观看视频| 成人美女网站在线观看视频| 午夜激情久久久久久久| 秋霞伦理黄片| 97在线视频观看| 精品一品国产午夜福利视频| 多毛熟女@视频| 日韩一区二区视频免费看| 国产成人精品婷婷| 人妻 亚洲 视频| 国国产精品蜜臀av免费| 亚洲欧美成人精品一区二区| 青青草视频在线视频观看| 亚洲精品一二三| 搡老乐熟女国产| 国产又色又爽无遮挡免| 青春草亚洲视频在线观看| 小蜜桃在线观看免费完整版高清| 亚洲精品国产色婷婷电影| 另类亚洲欧美激情| 久久人人爽av亚洲精品天堂 | 亚洲综合精品二区| 日本黄大片高清| 在线 av 中文字幕| 日韩av在线免费看完整版不卡| 丝袜脚勾引网站| 日本免费在线观看一区| 亚洲四区av| 身体一侧抽搐| 国产美女午夜福利| 蜜臀久久99精品久久宅男| 99热这里只有是精品在线观看| 亚洲av国产av综合av卡| 日本黄色片子视频| 精品国产乱码久久久久久小说| 久久久久久久国产电影| 亚洲综合色惰| 一级爰片在线观看| 国国产精品蜜臀av免费| 尾随美女入室| 精品少妇久久久久久888优播| 国产一区二区三区综合在线观看 | 亚洲第一区二区三区不卡| 99久国产av精品国产电影| 精品熟女少妇av免费看| 黄色视频在线播放观看不卡| 免费人妻精品一区二区三区视频| 欧美日韩亚洲高清精品| 色婷婷av一区二区三区视频| 国产色爽女视频免费观看| 国产精品一及| 日韩电影二区| 我要看日韩黄色一级片| 久久鲁丝午夜福利片| 久久6这里有精品| 性高湖久久久久久久久免费观看| 亚洲国产最新在线播放| 日本黄大片高清| 人妻少妇偷人精品九色| 日本与韩国留学比较| 热re99久久精品国产66热6| 国产真实伦视频高清在线观看| 国产精品久久久久久精品古装| 日韩中字成人| 国产深夜福利视频在线观看| 久久久成人免费电影| 内地一区二区视频在线| 一区二区三区精品91| 亚洲高清免费不卡视频| 18禁在线无遮挡免费观看视频| 国产高清三级在线| 国产精品成人在线| 欧美成人一区二区免费高清观看| 国产精品不卡视频一区二区| 一区二区av电影网| 国产91av在线免费观看| 中国三级夫妇交换| 精品一品国产午夜福利视频| 欧美性感艳星| 亚洲经典国产精华液单| 国产精品熟女久久久久浪| 亚洲av欧美aⅴ国产| 韩国av在线不卡| 亚洲精品日韩av片在线观看| 高清在线视频一区二区三区| 大片电影免费在线观看免费| 日日摸夜夜添夜夜添av毛片| 国产高清有码在线观看视频| 插逼视频在线观看| 国产美女午夜福利| 免费观看a级毛片全部| 国产91av在线免费观看| 国产欧美另类精品又又久久亚洲欧美| 久久久精品免费免费高清| 深爱激情五月婷婷| .国产精品久久| 日本欧美国产在线视频| 蜜臀久久99精品久久宅男| 如何舔出高潮| 草草在线视频免费看| 涩涩av久久男人的天堂| 有码 亚洲区| 人人妻人人爽人人添夜夜欢视频 | av国产精品久久久久影院| 国产男女超爽视频在线观看| 亚洲精品日本国产第一区| 日日啪夜夜撸| 纯流量卡能插随身wifi吗| 七月丁香在线播放| 赤兔流量卡办理| 精品少妇久久久久久888优播| 成人亚洲精品一区在线观看 | 欧美最新免费一区二区三区| 中文字幕精品免费在线观看视频 | 亚洲美女搞黄在线观看| 精品午夜福利在线看| 国产亚洲午夜精品一区二区久久| 国产成人午夜福利电影在线观看| 色婷婷久久久亚洲欧美| av在线蜜桃| 99热这里只有是精品在线观看| 精品国产乱码久久久久久小说| 老女人水多毛片| 国产国拍精品亚洲av在线观看| 老女人水多毛片| www.色视频.com| 一个人免费看片子| 国产成人精品婷婷| 男人添女人高潮全过程视频| 久久人人爽人人爽人人片va| 日韩在线高清观看一区二区三区| 在线观看免费高清a一片| 制服丝袜香蕉在线| 美女高潮的动态| 视频区图区小说| av在线app专区| 麻豆国产97在线/欧美| 久久av网站| 男男h啪啪无遮挡| 国产精品一区二区性色av| 麻豆精品久久久久久蜜桃| 少妇被粗大猛烈的视频| 新久久久久国产一级毛片| 多毛熟女@视频| 黄色配什么色好看| 在线 av 中文字幕| 一区在线观看完整版| 最近最新中文字幕大全电影3| 国产亚洲午夜精品一区二区久久| 国产欧美亚洲国产| 国产片特级美女逼逼视频| 欧美精品一区二区免费开放| 在线观看人妻少妇| 久久ye,这里只有精品| 亚洲色图av天堂| 麻豆成人午夜福利视频| 美女主播在线视频| 精品一区二区三卡| 亚洲精品久久午夜乱码| 国产有黄有色有爽视频| 建设人人有责人人尽责人人享有的 | 亚洲精品,欧美精品| 久久精品国产a三级三级三级| 色哟哟·www| 免费大片黄手机在线观看| 国产精品熟女久久久久浪| 黄片无遮挡物在线观看| av在线老鸭窝| 免费大片18禁| 国产精品久久久久久久电影| 亚洲精品第二区| 久久精品久久久久久噜噜老黄| 欧美97在线视频| 日本av手机在线免费观看| 夜夜看夜夜爽夜夜摸| 国产亚洲一区二区精品| 日日撸夜夜添| kizo精华| 国产在视频线精品| 黄片无遮挡物在线观看| 国产精品麻豆人妻色哟哟久久| 偷拍熟女少妇极品色| 在线天堂最新版资源| 国产精品一二三区在线看| 亚洲av.av天堂| 亚洲欧美成人精品一区二区| 精品久久久精品久久久| 国产黄频视频在线观看| 国产高清三级在线| 天美传媒精品一区二区| av在线蜜桃| 亚洲国产成人一精品久久久| 少妇的逼好多水| 亚洲中文av在线| 国产成人aa在线观看| 少妇人妻 视频| 久久久精品免费免费高清| 国产精品久久久久久av不卡| 日韩人妻高清精品专区| 丝瓜视频免费看黄片| 激情五月婷婷亚洲| 亚洲成人av在线免费| 少妇人妻精品综合一区二区| www.色视频.com| 日韩成人伦理影院| 2021少妇久久久久久久久久久| 亚洲熟女精品中文字幕| 国产日韩欧美在线精品| 久久6这里有精品| 欧美精品人与动牲交sv欧美| 亚洲不卡免费看| 噜噜噜噜噜久久久久久91| 欧美日韩综合久久久久久| 韩国高清视频一区二区三区| 色吧在线观看| 免费看光身美女| 国产人妻一区二区三区在| 亚洲美女视频黄频| 熟女av电影| 夜夜骑夜夜射夜夜干| 国产真实伦视频高清在线观看| 精品亚洲乱码少妇综合久久| 最近最新中文字幕免费大全7| 青春草视频在线免费观看| 亚洲天堂av无毛| 亚洲色图综合在线观看| 人妻制服诱惑在线中文字幕| av在线播放精品| 五月开心婷婷网| 激情 狠狠 欧美| 国产探花极品一区二区| 亚洲aⅴ乱码一区二区在线播放| 欧美最新免费一区二区三区| 春色校园在线视频观看| 久久久久久九九精品二区国产| 国产视频内射| 国内揄拍国产精品人妻在线| av在线老鸭窝| 国国产精品蜜臀av免费| 亚洲无线观看免费| 久久韩国三级中文字幕| av.在线天堂| 亚洲欧美日韩另类电影网站 | 2018国产大陆天天弄谢| 国产av一区二区精品久久 | 中文字幕制服av| 色吧在线观看| 午夜视频国产福利| 日本爱情动作片www.在线观看| 人人妻人人澡人人爽人人夜夜| 在现免费观看毛片| 亚洲精品视频女| 国产精品伦人一区二区| www.色视频.com| 美女中出高潮动态图| 成人18禁高潮啪啪吃奶动态图 | 精品一区二区三卡| 亚洲婷婷狠狠爱综合网| 亚洲精品,欧美精品| 美女中出高潮动态图| 一区二区三区免费毛片| 欧美日韩视频精品一区| 亚洲av中文字字幕乱码综合| 夜夜骑夜夜射夜夜干| 99热国产这里只有精品6| 亚洲精品一二三| 成年人午夜在线观看视频| 少妇人妻久久综合中文| 日本黄色日本黄色录像| 日韩 亚洲 欧美在线| 一本—道久久a久久精品蜜桃钙片| 黄色一级大片看看| 国产精品一二三区在线看|