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

    Investigation of the performance of CF3I/c-C4F8/N2 and CF3I/c-C4F8/CO2 gas mixtures from electron transport parameters

    2020-06-14 08:45:28RuishuangZHONG鐘蕊霜SuZHAO趙謖DengmingXIAO肖登明HuiWANG王輝XiuchenJIANG江秀臣ZhongminYU余鐘民andYunkunDENG鄧云坤
    Plasma Science and Technology 2020年5期
    關(guān)鍵詞:王輝

    Ruishuang ZHONG(鐘蕊霜),Su ZHAO(趙謖),Dengming XIAO(肖登明),4,Hui WANG (王輝),4,Xiuchen JIANG (江秀臣),Zhongmin YU (余鐘民) and Yunkun DENG (鄧云坤)

    1 Department of Electrical Engineering,Shanghai Jiao Tong University,Shanghai 200030,People’s Republic of China

    2 State Grid Shanghai Electric Power Company,Shanghai 200122,People’s Republic of China

    3 Yunnan Power Grid Corporation,Kunming 650217,People’s Republic of China

    4 Authors to whom any correspondence should be addressed.

    Abstract

    Keywords:CF3I/c-C4F8/N2 gas mixtures,CF3I/c-C4F8/CO2 gas mixtures,Boltzmann method,electron transport parameters,synergistic effect

    1.Introduction

    Sulfur hexafluoride (SF6) is a colorless,odorless,non-toxic,non-flammable gas with stable chemical properties and excellent insulation properties.It is widely used in power systems[1].However,SF6is also a strong greenhouse gas;its global warming potential (GWP) is approximately 23 900.To maintain sustainable development of the environment and resources,SF6was classified as a greenhouse gas under the Kyoto Protocol in 1997,which required a gradual reduction of its use [2].In recent years,researchers have found a number of alternatives to SF6.

    Figure 1.The molecular structure of gases:(a) CF3I;(b) c-C4F8.

    Trifluoroiodomethane (CF3I) has been proposed as a potential replacement for SF6owing to its excellent physicochemical properties and outstanding dielectric strength.The insulation strength of CF3I is approximately 1.2 times that of SF6in a slightly non-uniform electrical field.It has good oil solubility and material compatibility.Its GWP is approximately 5 and its ozone destruction potential (ODP) is approximately 0.008,which is negligible [3,4].Its molecular structure is shown in figure 1(a).The carbon–iodine chemical bond in the molecular structure of CF3I is relatively active,and easily broken under solar radiation,leading to the decomposition of CF3I.This property limits its diffusion to the stratosphere and reduces its greenhouse effect.However,the liquefaction temperature of CF3I is approximately -22.5°C under normal temperature and pressure.A buffer gas is needed to reduce its liquefaction temperature,and then CF3I gas mixtures can be suitable for applications in electrical equipment [5].

    Perfluorocyclobutane (c-C4F8) is another kind of potential electronegative replacement gas for SF6.Its electric insulation performance is better than that of SF6.It has a GWP of 8 700,which is only 36% of that of SF6,and its ODP is zero [6].The molecular structure of c-C4F8is shown in figure 1(b).Owing to its symmetrical structure,and having two more fluorine atoms than SF6,its electronegativity is stronger than that of SF6and it is chemically stable.However,similar to CF3I,c-C4F8also has a high liquefaction temperature.The materials must be mixed with buffer gases to reduce the liquefaction temperature [7,8].

    In research into electronegative gas mixtures,N2and CO2are the most commonly used buffer gases.They have stable physico-chemical properties and their liquefaction temperatures are -195.8 °C and -78.5 °C [9].They can greatly improve the liquefaction temperature performance of electronegative gases.

    In previous research,De Urquijo et al used a pulsed Townsend (PT) method to measure the microscopic parameters of CF3I and CF3I-N2[10].Katagiri et al found that the dielectric strength of CF3I-CO2(30%–70%)is approximately 0.75 to 0.8 times of that of SF6[5].Li et al analyzed the insulation characteristics of CF3I mixtures with CF4,CO2,N2,O2,and air,and found that CF3I/N2and CF3I/air have better insulation performance than the other three mixtures [2].Zhao et al studied CF3I gas mixtures and found that the synergistic effect coefficient of CF3I/CO2gas mixtures was higher than that of CF3I/N2[11].For c-C4F8,Xiao et al analyzed the insulation characteristics of c-C4F8gas mixtures by Monte Carlo and Boltzmann methods [7,12].Man et al studied the insulating properties of c-C4F8/N2and c-C4F8/CO2under an alternating current and found that 30 % c-C4F8/70% N2and 30 % c-C4F8/70% CO2mixed gases had 75% and 79% respectively of the capacity of pure SF6at 0.1 MPa [13].For ternary gas mixtures,Xu et al studied the CF3I binary and ternary gas mixtures with SF6,CO2,N2,and CF4,and found synergetic effects in CF3I with N2and CO2,whereas the CF3I/SF6system had a slight negative synergetic effect [14].Cheng et al studied the breakdown characteristics of binary and ternary gas mixtures of CF3I/SF6/N2and CF3I/SF6/CO2under DC voltage[15].At present,there has been little research on ternary mixed gases of CF3I and c-C4F8.Therefore,studies of the electron transport parameters of ternary gas mixtures of CF3I and c-C4F8might provide basic data to support the application of these gases as alternatives to SF6.

    In this paper,based on a Boltzmann simulation method,CF3I/c-C4F8/N2and CF3I/c-C4F8/CO2gas mixtures are investigated in terms of the electron energy distribution function(EEDF),reduced effective ionization coefficient,critical reduced electrical field strength,electric field sensitivity coefficient,electron drift velocity,and synergistic effect.

    2.Theoretic framework

    2.1.Boltzmann equation method

    The process of solving the Boltzmann equation is to build a self-consistent model reflecting the interaction between an electric field and particles.This complex model takes into account both elastic and inelastic collision processes,and the effects of external electric fields on charged particles[16].By solving the Boltzmann equation,the EEDFs of gas mixtures at different volume ratios are obtained,and the electron group parameters and transport parameters in the discharge process of the gas mixtures are calculated.The critical reduced field strength of the gas mixtures is determined.Next,we analyze the insulation properties of the gas mixtures.

    For a particle that exists in a six-dimensional phase space(r,v),we define the distribution functionf(r,v,t).Its physical meaning is:the number of particles in the six-dimensional volume elementr vd d3 3 around the timetand the phase space (r,v).If the particle is an electron in a weakly ionized gas,its distribution can be described by the Boltzmann equation as follows:

    where r is the electron space position,v is the electron velocity,eis the electron charge (approximately ×1.6 10-19C),mis the electron mass,E is the strength of the electric field in space,andandC f[ ]are collision terms,defining the rate of change of the function of the electron during various collisions [17].

    Figure 2.The electron impact cross-sections of gases:(a) CF3I;(b) c-C4F8.

    To obtain the electron transport parameters,it is necessary to solve the Boltzmann equations.In this paper,the Boltzmann equations were solved with the use of BOLSIG + software by a two-term approximation method.This method has high accuracy and a rapid operation speed,which makes it the most common method for solving Boltzmann equations.The two-term approximation method requires the conditions that the electric field is uniformly distributed in space and the collision probability is equal.Electrons in velocity space are symmetric with respect to the direction of the electric field,and the distribution in the position space changes only along the direction of the electric field[18,19].In a spherical coordinate system in velocity space,equation (1) is written as follows:

    wherevis the magnitude of velocity,θ is the angle between the velocity and the electric field,andzis the position of the electric field.

    We can expand the energy distribution functionfas a Legendre polynomial aboutθcos .In a two-term approximation method,only two terms are expanded.

    wheref0represents the isotropic part andf1represents the anisotropic part [20].

    2.2.Collision cross-sections

    To calculate the Boltzmann equation and solve the electron transport parameters,the gas collision cross-section data are required,including the attachment cross-sectionQa,the ionization cross-sectionQi,the elastic cross-sectionQel,and the electronic excitation cross-sectionQex.Most of the collision cross-section data(CF3I,SF6,CO2,N2)in this paper are taken from the LXcat website [21],which provides a comprehensive and reliable database for such data.The crosssection data of CF3I are taken from the TRINITI database,as shown in figure 2(a).For c-C4F8,the attachment cross-sectionQaand ionization cross-sectionQiare taken from [22],whereas the elastic cross-sectionQeland electronic excitation cross-sectionQexare taken from [7],as shown in figure 2(b).

    According to the selected cross-section,and based on Boltzmann calculations,the reduced effective ionization coefficients(α - η)/N of pure CF3I and c-C4F8are shown in figures 3(a) and (b),respectively.Data obtained from experiments and other calculation methods [10,23–27]are listed in table 1 for comparison.According to table 1,for CF3I,the critical reduced electrical field strength (E/N)crof pure CF3I in our work is 432.7 Townsends (Td),approximately 1.22 times that of SF6,and similar to the results from the PT method,time-of-flight(TOF)method,and steady-state Townsend (SST) method.For c-C4F8,the (E/N)crof pure c-C4F8of is 406.2 Td,approximately 1.14 times that of SF6,and similar to the results from Monte Carlo simulation,and SST and PT methods.The calculated results in this paper are consistent with the experimental results.The microscopic discharge parameters in this paper are accurately solved with the use of the two-term Boltzmann equation under the current calculation conditions,and the results can be used for analysis of gas discharge processes with high reliability.

    2.3.Synergistic effect

    The electrical characteristics of insulating gas mixtures are often better than those of a single gas,owing to synergies between the gases.To compare and analyze the synergies of ternary gas mixtures more intuitively,a synergies factorSis introduced in this paper,defined as follows [11]:

    whereUmixis the breakdown voltage of gas mixtures,ξiis the proportion of the gas in the gas mixtures,andUiis the breakdown voltage of a single gas in the gas mixtures.

    Figure 3.Comparison of reduced effective ionization coefficients of pure gases:(a) CF3I;(b) c-C4F8.

    Table 1.Comparison data of (E/N)cr of pure gases.

    If >S0,the actual breakdown voltage of the gas mixtures is greater than the theoretical breakdown voltage,indicating that the gas mixtures have a synergistic effect or even a positive synergistic effect.A largerSvalue indicates a more pronounced synergistic effect.If=S0,there is no synergistic effect and the breakdown voltage of the gas mixtures varies linearly.If <S0,the gas mixtures have a negative synergistic effect(a Penning effect),which is relevant in the field of plasma applications.

    3.Results and discussion

    In this paper,the interactions between two strongly electronegative gases in ternary gas mixtures are studied,and the effect of the ratio of the strongly electronegative gas on the buffer gas is discussed.Hence,in the following the ratio of buffer gas is fixed.According to the requirements of the liquefaction temperature and insulation strength,the proportion of buffer gas is determined to be 70%,which generally meets the requirements for applications in power systems.A parameter=kiχ( c‐ C4F8) /(χ( c‐ C4F8) +χ(C F3I) +χ(b uffer gas))(i=1,buffer gas is N2;=i2,buffer gas is CO2)is used to represent the ratio of c-C4F8.The temperature is set at 300 K.

    3.1.EEDF of gas mixtures

    The EEDFs of CF3I/c-C4F8/70% N2and CF3I/c-C4F8/70%CO2gas mixtures at E/N=300 Td are shown in figure 4.Only slight differences between the EEDFs were observed,and when the electron energy was in the range of 0–1 eV,the number of low-energy electrons increased as the c-C4F8content increased and the CF3I content decreased.We attribute this phenomenon to c-C4F8having a stronger adsorption capacity and to many electrons being concentrated in the low energy region owing to absorption.Conversely,in the highenergy region,the curves of 0%and 30%are higher than any other curves.This result might be attributed to synergistic effects in the gas mixtures.Moreover,the curve of SF6is lower than other curves in the low-energy region but higher than other curves in the high-energy region,which indicates that the CF3I/c-C4F8gas mixtures are more electronegative than the SF6gas mixtures.

    The EEDF values of CF3I/c-C4F8/70% N2and CF3I/c-C4F8/70% CO2having the same c-C4F8contents(0%,10%,20%,30%) are compared with 30% SF6/70% N2and 30% SF6/70% CO2,as shown in figure 5.The curves of CO2gas mixtures are lower than the curves of N2gas mixtures of the same proportion in the low-energy region.The reason for this is that the N2molecule has a large vibrational excitation cross-section in the low-energy region,and frequent excitation collisions decrease the kinetic energy of free electrons.Thus,CO2is more effective than N2in reducing electron energy in the high-energy region.

    3.2.Electron transport parameters of gas mixtures

    The reduced effective ionization coefficient (α - η)/N,critical reduced electrical field strength (E/N)cr,electric field sensitivity coefficient c,and electron drift velocity Vecan be determined from the EEDFs of gas mixtures.

    Figure 4.The EEDF of CF3I/c-C4F8/70% N2 and CF3I/c-C4F8/70% CO2 gas mixtures at 300 Td.(a) CF3I/c-C4F8/70% N2,(b) CF3I/c-C4F8/70% CO2.

    Figure 5.The comparison of EEDF of CF3I/c-C4F8/70% N2 and CF3I/c-C4F8/70% CO2 gas mixtures at 300 Td.

    Figure 6.Relationship between effective ionization coefficient and E/N in CF3I/c-C4F8/70% N2 and CF3I/c-C4F8/70% CO2 gas mixtures at different ratio k.(a) CF3I/c-C4F8/70% N2;(b) CF3I/c-C4F8/70% CO2.

    The relationship between the effective ionization coefficient and E/N in CF3I/c-C4F8/70% N2and CF3I/c-C4F8/70% CO2gas mixtures at different k ratios is shown in figure 6.With the increase of E/N,free electrons in the gas mixtures can obtain more energy under acceleration by an electric field.As a result,the probability of collision ionization increases and the effective ionization coefficient of the gas mixtures also increases.Research on binary gas mixtures[2]has shown that as the proportion of the strong electronegative gas in the gas mixtures increases (i.e.,giving stronger electronegativity of the gas mixtures),(α - η)/N decreases under the conditions of the same E/N.We apply this conclusion to analyze the data in the figures,and 30%SF6/70% N2gas mixtures show stronger electronegativity than CF3I/c-C4F8/70% N2gas mixtures,whereas CF3I/c-C4F8/70% CO2gas mixtures show stronger electronegativity than 30%SF6/70%N2when E/N is approximately 300–600 Td.

    Figure 7.The relationship between (E/N)cr and the c-C4F8 content.

    When the reduced effective ionization coefficient(α - η)/N=0,the reduced electrical field strength E/N is the critical reduced electrical field strength (E/N)cr.Figure 7 shows the relationship between (E/N)crand the content of c-C4F8.The(E/N)crof CF3I/c-C4F8/70%N2gas mixtures is higher than that of CF3I/c-C4F8/70%CO2,and the difference between these two gas mixtures gradually increases as the proportion of c-C4F8gas increases.Moreover,the (E/N)crof CF3I/c-C4F8/70% N2gas mixtures varies little as the proportion of c-C4F8changes,and is approximately 0.88 times that of 30%bSF6/70%bN2.When the c-C4F8content was 17%,the insulation strength of CF3I/c-C4F8/70%N2was the highest.Furthermore,CF3I/c-C4F8/70% CO2had excellent insulation performance.The (E/N)crof CF3I/c-C4F8/70%CO2was higher than that of 30% SF6/70% SF6for c-C4F8contents of less than 17%.

    Table 2 lists values of the electric field sensitivity coefficient c of the CF3I/c-C4F8/70% N2and CF3I/c-C4F8/70%CO2gas mixtures.The electric field sensitivity of CF3I/c-C4F8/70% N2gas mixtures is lower than that of CF3I/c-C4F8/70% CO2gas mixtures,indicating that CF3I/c-C4F8/70% N2gas mixtures have excellent performance in terms of electric field sensitivity.Moreover,a higher content of c-C4F8promoted a lower sensitivity of the gas mixtures to an electric field.Compared with SF6gas mixtures,for CF3I/c-C4F8/70% N2gas mixtures at a c-C4F8content greater than 10%,the electric field sensitivity of the CF3I/c-C4F8/70% N2gas mixtures was higher than that of 30%SF6/70% N2gas mixtures (c=20.1 kV-1).For CF3I/c-C4F8/70% CO2gas mixtures,the electric field sensitivity of the CF3I/c-C4F8/70% CO2gas mixtures was higher than that of 30% SF6/70% CO2gas mixtures(c=28.6 kV-1) at all ratios.

    The electron drift velocities Vein gas mixtures,as a function of E/N at different c-C4F8gas mixture ratios k,are given in figure 8.As for CF3I/c-C4F8/70% N2,the electron drift velocity of the gas mixtures of all ratios increased almost linearly as E/N increased.Conversely,as the ratio of c-C4F8increased the electron drift velocity rapidly decreased,particularly in areas where E/N was high.The reason for this is that when E/N is high,after intense collisions between electrons and molecules,the electron energy decreases and enters a region with a large adsorption area by CF3I and c-C4F8.This enhances the adsorption capacity of gas mixtures and weakens the drift motion of electrons in the electric field.Moreover,in the range of 300–900 Td,the electron drift velocity of 30% SF6/70% N2overlaps with that of 5%c-C4F8/25%CF3I/70%N2gas mixtures.For CF3I/c-C4F8/70%CO2gas mixtures,the electron drift velocity is discussed based on the three regimes.In the range of 100–250 Td,when the ratio of c-C4F8increases,the electron drift velocity slowly increases.A reversal occurs at approximately 250 Td.In the range of 300–900 Td,the electron drift velocity rapidly decreases as the ratio of c-C4F8increases.The curve of electron drift velocity of 30%SF6/70%CO2is between the electron drift velocity curves of 5% CF3I/25% c-C4F8/70% CO2and 10% CF3I/20%c-C4F8/70% CO2.

    3.3.Synergistic effect of gas mixtures

    The synergistic effect,as described by the factor S,increased in the ternary gas mixtures (E/N)cras the content of c-C4F8was increased.We studied the synergistic effects of ternary gas mixtures;hence,formula (4) is rewritten as:

    whereξ1,ξ2,ξ3are the proportions of the three gases in the gas mixtures,andξ1+ξ2+ξ3=1,[(E/N)cr]1,[(E/N)cr]2,[(E/N)cr]3are the critical reduced electrical field strengths of the three single gases under the same conditions.(ξ1[ (E/N)cr]1+ξ2[(E/N)cr]2+ξ3[(E/N)cr]3) represents the theoretical critical reduced electrical field strength of the gas mixtures,calculated according to the weighting of the mixing proportion.The above equation calculates the synergy factor S of ternary gas mixtures.

    Figure 9 shows the relationship between the synergy factor S and the c-C4F8content.The synergistic effect of CF3I/c-C4F8/70% CO2was much greater than that of CF3I/c-C4F8/70%N2when the c-C4F8content was less than 25%.Moreover,as the c-C4F8content increases,the gap of the synergy factors of the two gas mixtures decreased.Furthermore,the maximum synergy factor value of the CF3I/c-C4F8/70%N2gas mixtures occurred when the c-C4F8content was approximately 20%.The synergistic effect for CF3I/c-C4F8/70% CO2was most obvious when the c-C4F8content was approximately 10%.

    Table 2.Values of electric field sensitivity coefficient c of gas mixtures.

    Figure 8.Electron drift velocities Ve in gas mixtures as a function of E/N at different c-C4F8 gas mixture ratio k.(a) CF3I/c-C4F8/70%N2.(b) CF3I/c-C4F8/70% CO2.

    3.4.Comparison of electron transport parameters

    Table 3 includes a comparison of parameters of the gas mixtures.Taken together,for CF3I/c-C4F8/70% N2gas mixtures,when the content of c-C4F8is the 15%–20%range,the gas mixtures show the best insulation performance.However,for CF3I/c-C4F8/70% CO2,the gas mixtures have excellent performance when the c-C4F8content is 10%–15%.

    Figure 9.The relationship between synergistic factor S and the c-C4F8 content.

    4.Conclusions

    In this paper,the insulating properties and synergistic effects of ternary gas mixtures CF3I/c-C4F8/N2and CF3I/c-C4F8/CO2were analyzed.The gas mixtures of c-C4F8and CF3I with ratios of 0%,5%,10%,15%,20%,25%,and 30%are calculated with the use of two-term approximation Boltzmann equations.

    The main conclusions are as follows.

    (1) The critical electric field strength(E/N)crof CF3I/c-C4F8/70% N2is 251.8 Td when the proportion of c-C4F8is 17%,which is the highest among the systems studied here.This value is approximately 0.89 times that of 30%SF6/70% N2.When the content of CF3I was higher than 17%,the CF3I/c-C4F8/70% CO2critical electric field strength was higher than that of 30% SF6/70% CO2.

    (2) For electric field sensitivity,an increase of the content of c-C4F8decreased the sensitivity of the gas mixtures to the electric field.Moreover,the electric field sensitivities of the gas mixtures with N2buffer gas were lower than those with CO2.In terms of electron drift velocity,when the content of c-C4F8was greater than 5%,the electron velocity of CF3I/c-C4F8/70% N2was lower than that of 30%SF6/70% N2.For CF3I/c-C4F8/70% CO2,when the content of c-C4F8was greater than 10%,the electron drift velocity of CF3I/c-C4F8/70%CO2was lower than that of 30% SF6/70% CO2.

    Table 3.Comparison of parameters of the gas mixtures.

    (3) The synergistic effect of CF3I/c-C4F8/70% CO2was greater than that of CF3I/c-C4F8/70% N2when the CF3I content was greater than 5%.The synergistic effects for CF3I/c-C4F8/70% N2were most obvious when the c-C4F8content was approximately 20%,and for CF3I/c-C4F8/70% CO2when the c-C4F8content was approximately 10%.

    In summary,in our studies of electron transport parameters and synergistic effects,CF3I/c-C4F8/70% N2shows the best insulation performance when the c-C4F8content is in the 15%–20% range.For CF3I/c-C4F8/70% CO2,when the c-C4F8content is 10%–15%,the gas mixtures have excellent performance.These gas systems might be potential alternatives to SF6as new insulating media.

    Acknowledgments

    This work is supported by National Natural Science Foundation of China (No.51337006).

    ORCID iDs

    猜你喜歡
    王輝
    本期面孔
    遼河(2022年12期)2023-01-29 13:24:58
    Spin transport properties for B-doped zigzag silicene nanoribbons with different edge hydrogenations
    劇作家王輝
    火花(2021年10期)2021-11-04 09:23:52
    Characterization of size effect of natural convection in melting process of phase change material in square cavity?
    Chinese Acupuncture
    竹之韻
    A novel male-associated marker for sex identifi cation in Ussuri catfishPseudobagrus ussuriensis*
    The Thought on PPP in China
    弟子規(guī)
    山西老西兒
    国产精品久久久久久久电影 | 亚洲va日本ⅴa欧美va伊人久久| 日本黄大片高清| 最新中文字幕久久久久| 欧美性感艳星| 亚洲精品456在线播放app | 成年女人毛片免费观看观看9| 制服人妻中文乱码| 又黄又粗又硬又大视频| 日韩精品青青久久久久久| 久久精品国产清高在天天线| 免费搜索国产男女视频| 日本在线视频免费播放| 免费无遮挡裸体视频| 久久国产乱子伦精品免费另类| 90打野战视频偷拍视频| xxxwww97欧美| 熟妇人妻久久中文字幕3abv| 色在线成人网| 亚洲国产精品成人综合色| 又爽又黄无遮挡网站| 亚洲成av人片在线播放无| 大型黄色视频在线免费观看| 淫秽高清视频在线观看| 桃色一区二区三区在线观看| 中文字幕高清在线视频| 日韩欧美国产一区二区入口| 99国产综合亚洲精品| 热99re8久久精品国产| 国产探花在线观看一区二区| 久久精品人妻少妇| 舔av片在线| 三级毛片av免费| 窝窝影院91人妻| 身体一侧抽搐| 麻豆国产97在线/欧美| 亚洲最大成人手机在线| 国产色爽女视频免费观看| 三级男女做爰猛烈吃奶摸视频| 男人舔女人下体高潮全视频| 亚洲欧美一区二区三区黑人| 法律面前人人平等表现在哪些方面| 国产探花极品一区二区| 久久这里只有精品中国| 亚洲国产精品成人综合色| 欧美日韩亚洲国产一区二区在线观看| 欧美最新免费一区二区三区 | 十八禁网站免费在线| 男人和女人高潮做爰伦理| 日本一本二区三区精品| 精品一区二区三区人妻视频| 老汉色av国产亚洲站长工具| x7x7x7水蜜桃| 日本熟妇午夜| 久久精品国产清高在天天线| 久久九九热精品免费| 亚洲人成网站高清观看| 午夜福利高清视频| 舔av片在线| 国产视频内射| 欧美乱色亚洲激情| 亚洲成人中文字幕在线播放| 亚洲国产精品成人综合色| 亚洲久久久久久中文字幕| 日韩精品中文字幕看吧| 久久这里只有精品中国| 十八禁人妻一区二区| 国产精品 欧美亚洲| 成人国产综合亚洲| 久久精品夜夜夜夜夜久久蜜豆| 我要搜黄色片| 国产av一区在线观看免费| 亚洲专区国产一区二区| 国产私拍福利视频在线观看| 99在线视频只有这里精品首页| 精品国内亚洲2022精品成人| 国产三级在线视频| 亚洲七黄色美女视频| 午夜福利成人在线免费观看| 在线观看66精品国产| 免费av毛片视频| 18禁裸乳无遮挡免费网站照片| 国产 一区 欧美 日韩| 亚洲精品在线美女| 淫秽高清视频在线观看| 99热这里只有精品一区| 国产黄片美女视频| 亚洲av成人av| 亚洲一区高清亚洲精品| 久久久久国产精品人妻aⅴ院| 少妇丰满av| 母亲3免费完整高清在线观看| 色尼玛亚洲综合影院| 女同久久另类99精品国产91| 97超视频在线观看视频| 亚洲欧美激情综合另类| 网址你懂的国产日韩在线| 香蕉av资源在线| 欧美黑人巨大hd| 中文字幕人妻熟人妻熟丝袜美 | 狂野欧美白嫩少妇大欣赏| 国产欧美日韩精品一区二区| av国产免费在线观看| 久久亚洲真实| 99精品在免费线老司机午夜| 亚洲 国产 在线| or卡值多少钱| 色尼玛亚洲综合影院| 一区二区三区激情视频| 九色成人免费人妻av| 三级毛片av免费| 国产伦一二天堂av在线观看| 他把我摸到了高潮在线观看| 一进一出抽搐动态| 伊人久久精品亚洲午夜| 美女黄网站色视频| 午夜激情福利司机影院| 成年版毛片免费区| АⅤ资源中文在线天堂| 在线看三级毛片| 国产伦一二天堂av在线观看| 美女黄网站色视频| 午夜激情欧美在线| 欧美精品啪啪一区二区三区| 天堂动漫精品| 亚洲一区高清亚洲精品| 国语自产精品视频在线第100页| 国产欧美日韩一区二区三| 日本一二三区视频观看| 一区二区三区免费毛片| 高清在线国产一区| 亚洲 国产 在线| 欧美日韩乱码在线| 国产高清videossex| 亚洲一区二区三区色噜噜| 国产一区二区亚洲精品在线观看| 校园春色视频在线观看| 欧美一区二区精品小视频在线| 乱人视频在线观看| 中文字幕人妻熟人妻熟丝袜美 | 99久国产av精品| 国产精品一区二区三区四区免费观看 | 操出白浆在线播放| 亚洲色图av天堂| 美女黄网站色视频| 少妇的丰满在线观看| 网址你懂的国产日韩在线| 欧美黄色片欧美黄色片| 少妇的丰满在线观看| 国产老妇女一区| 国产高清视频在线观看网站| 成年女人毛片免费观看观看9| 欧美日韩综合久久久久久 | 亚洲aⅴ乱码一区二区在线播放| 18禁黄网站禁片午夜丰满| 亚洲一区二区三区不卡视频| 国产免费av片在线观看野外av| 亚洲欧美日韩高清专用| 欧美精品啪啪一区二区三区| 亚洲精品色激情综合| 国产成人av教育| 两性午夜刺激爽爽歪歪视频在线观看| 免费高清视频大片| 国产精品久久电影中文字幕| 精品一区二区三区人妻视频| 99久久无色码亚洲精品果冻| 在线国产一区二区在线| 国产成人系列免费观看| 亚洲av免费高清在线观看| 久久精品综合一区二区三区| 三级毛片av免费| 中文字幕av在线有码专区| 大型黄色视频在线免费观看| 91麻豆精品激情在线观看国产| 日本a在线网址| xxxwww97欧美| 国产亚洲精品av在线| 在线观看av片永久免费下载| 亚洲国产欧洲综合997久久,| 天天躁日日操中文字幕| 亚洲欧美激情综合另类| 真实男女啪啪啪动态图| 日韩精品中文字幕看吧| 欧美bdsm另类| av视频在线观看入口| 亚洲男人的天堂狠狠| 老熟妇乱子伦视频在线观看| 国产一区二区在线观看日韩 | 别揉我奶头~嗯~啊~动态视频| 国产精品野战在线观看| 大型黄色视频在线免费观看| 国产一区二区三区在线臀色熟女| 成年免费大片在线观看| 一级毛片高清免费大全| 身体一侧抽搐| 美女高潮喷水抽搐中文字幕| 听说在线观看完整版免费高清| 精品欧美国产一区二区三| 蜜桃久久精品国产亚洲av| 又黄又爽又免费观看的视频| 国产午夜精品久久久久久一区二区三区 | 国产高潮美女av| 亚洲国产欧洲综合997久久,| 床上黄色一级片| 两个人视频免费观看高清| 91麻豆av在线| 特级一级黄色大片| 神马国产精品三级电影在线观看| 精品久久久久久久毛片微露脸| 婷婷丁香在线五月| 精品人妻偷拍中文字幕| 一本久久中文字幕| 少妇高潮的动态图| 18禁在线播放成人免费| 免费人成在线观看视频色| 91av网一区二区| 禁无遮挡网站| 99国产精品一区二区三区| 五月玫瑰六月丁香| 日韩欧美在线二视频| 成人国产一区最新在线观看| 日本一二三区视频观看| 又黄又爽又免费观看的视频| 俄罗斯特黄特色一大片| av专区在线播放| 一a级毛片在线观看| 亚洲成人久久爱视频| 激情在线观看视频在线高清| 欧美性感艳星| 国产av麻豆久久久久久久| 丁香欧美五月| 美女 人体艺术 gogo| 三级男女做爰猛烈吃奶摸视频| 久久久久久久午夜电影| 国产高清三级在线| 欧美在线黄色| 国产视频内射| 欧美一区二区精品小视频在线| 亚洲精华国产精华精| 久久精品国产综合久久久| 久久中文看片网| 国产97色在线日韩免费| 亚洲 国产 在线| 18禁裸乳无遮挡免费网站照片| 国产精品99久久久久久久久| 久久精品综合一区二区三区| www.熟女人妻精品国产| 国产精品香港三级国产av潘金莲| 成年免费大片在线观看| 欧美一区二区精品小视频在线| 中亚洲国语对白在线视频| 午夜日韩欧美国产| АⅤ资源中文在线天堂| 久久国产精品人妻蜜桃| 色综合亚洲欧美另类图片| 首页视频小说图片口味搜索| 夜夜看夜夜爽夜夜摸| 黄片小视频在线播放| 丝袜美腿在线中文| 日本黄色片子视频| 99国产精品一区二区三区| 欧美bdsm另类| 最后的刺客免费高清国语| 久久久久久久亚洲中文字幕 | 日本五十路高清| 制服丝袜大香蕉在线| 欧美另类亚洲清纯唯美| 免费av不卡在线播放| 久久精品国产清高在天天线| 给我免费播放毛片高清在线观看| 99久久精品一区二区三区| 神马国产精品三级电影在线观看| 乱人视频在线观看| 狂野欧美激情性xxxx| 老熟妇乱子伦视频在线观看| 国产91精品成人一区二区三区| 精品一区二区三区视频在线 | 九九热线精品视视频播放| 欧美黄色淫秽网站| 国产乱人伦免费视频| 久久久久久久久大av| 禁无遮挡网站| a级一级毛片免费在线观看| 日韩欧美国产一区二区入口| 叶爱在线成人免费视频播放| 亚洲国产精品久久男人天堂| 亚洲欧美日韩无卡精品| 国产日本99.免费观看| 麻豆久久精品国产亚洲av| 嫩草影院入口| 亚洲国产精品合色在线| 观看美女的网站| 美女被艹到高潮喷水动态| 黄片小视频在线播放| 亚洲国产中文字幕在线视频| 亚洲精品成人久久久久久| 非洲黑人性xxxx精品又粗又长| 999久久久精品免费观看国产| 老鸭窝网址在线观看| 欧美乱色亚洲激情| 黄色视频,在线免费观看| 丁香六月欧美| 少妇的逼好多水| 中文字幕人妻熟人妻熟丝袜美 | 亚洲欧美日韩高清在线视频| 99热只有精品国产| 国产高清三级在线| 亚洲av美国av| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 国产私拍福利视频在线观看| 国产高清videossex| 男女下面进入的视频免费午夜| 性欧美人与动物交配| 国产精品影院久久| 老司机午夜十八禁免费视频| 日韩精品青青久久久久久| 久9热在线精品视频| 女同久久另类99精品国产91| www.www免费av| 一个人看的www免费观看视频| 免费观看人在逋| 少妇的逼水好多| 久久精品影院6| 国产 一区 欧美 日韩| 男女做爰动态图高潮gif福利片| 精品午夜福利视频在线观看一区| 3wmmmm亚洲av在线观看| 特大巨黑吊av在线直播| 欧美bdsm另类| 亚洲国产中文字幕在线视频| 久久精品国产亚洲av香蕉五月| 成人亚洲精品av一区二区| 亚洲自拍偷在线| aaaaa片日本免费| 日韩人妻高清精品专区| 熟女人妻精品中文字幕| 好男人电影高清在线观看| 亚洲欧美日韩高清专用| 色哟哟哟哟哟哟| 久久九九热精品免费| 女同久久另类99精品国产91| 桃红色精品国产亚洲av| 亚洲 欧美 日韩 在线 免费| 午夜福利视频1000在线观看| 动漫黄色视频在线观看| 一本一本综合久久| 欧美一区二区精品小视频在线| 好男人在线观看高清免费视频| 熟女电影av网| 给我免费播放毛片高清在线观看| 亚洲成人久久性| 午夜精品在线福利| 国产高潮美女av| 99在线视频只有这里精品首页| 亚洲中文字幕一区二区三区有码在线看| 久久久久九九精品影院| 亚洲国产高清在线一区二区三| 十八禁人妻一区二区| 免费av不卡在线播放| 天天躁日日操中文字幕| 99视频精品全部免费 在线| 亚洲精品久久国产高清桃花| 亚洲人成网站在线播放欧美日韩| 午夜福利视频1000在线观看| 身体一侧抽搐| 一进一出抽搐动态| 国产精品久久久久久亚洲av鲁大| 亚洲成人久久爱视频| 国内少妇人妻偷人精品xxx网站| 国产淫片久久久久久久久 | 免费在线观看亚洲国产| 2021天堂中文幕一二区在线观| 在线免费观看不下载黄p国产 | 在线观看免费视频日本深夜| 亚洲五月婷婷丁香| 不卡一级毛片| 国产精品美女特级片免费视频播放器| 97碰自拍视频| 美女高潮的动态| 精品久久久久久,| 三级毛片av免费| 国产精品影院久久| 亚洲精品456在线播放app | 亚洲 国产 在线| 美女 人体艺术 gogo| 亚洲精品在线观看二区| 性欧美人与动物交配| 丁香欧美五月| 国内久久婷婷六月综合欲色啪| 精品国产美女av久久久久小说| 日本在线视频免费播放| 搡女人真爽免费视频火全软件 | 国产不卡一卡二| 亚洲男人的天堂狠狠| 黄色女人牲交| 精品一区二区三区av网在线观看| 九色成人免费人妻av| 精品久久久久久久人妻蜜臀av| 亚洲无线在线观看| 无限看片的www在线观看| 91在线精品国自产拍蜜月 | 麻豆国产av国片精品| 人妻夜夜爽99麻豆av| 麻豆成人午夜福利视频| 午夜福利在线在线| 亚洲精品456在线播放app | 村上凉子中文字幕在线| 老司机福利观看| 国产精品av视频在线免费观看| 国产在视频线在精品| 色吧在线观看| 丁香六月欧美| 在线观看美女被高潮喷水网站 | 国产一区二区三区视频了| 久9热在线精品视频| 又粗又爽又猛毛片免费看| 亚洲成a人片在线一区二区| 在线视频色国产色| 国产久久久一区二区三区| 色吧在线观看| a级一级毛片免费在线观看| 天天添夜夜摸| 少妇高潮的动态图| 免费观看的影片在线观看| 身体一侧抽搐| 99久久精品国产亚洲精品| 亚洲成人精品中文字幕电影| 成熟少妇高潮喷水视频| 国产精品一及| 窝窝影院91人妻| 丰满乱子伦码专区| 久久精品国产综合久久久| 18禁美女被吸乳视频| 国产精品99久久99久久久不卡| 亚洲五月婷婷丁香| 欧美+亚洲+日韩+国产| 国产精品亚洲av一区麻豆| 国产精品香港三级国产av潘金莲| 黑人欧美特级aaaaaa片| 色综合站精品国产| 欧美+日韩+精品| 最近最新中文字幕大全免费视频| 一级黄色大片毛片| 亚洲人成电影免费在线| 成人三级黄色视频| 熟女少妇亚洲综合色aaa.| 久久中文看片网| 成人欧美大片| 日韩欧美免费精品| 亚洲av二区三区四区| 日韩欧美精品v在线| 国产成人欧美在线观看| 免费av毛片视频| 精品国内亚洲2022精品成人| 超碰av人人做人人爽久久 | 欧美性猛交╳xxx乱大交人| 在线十欧美十亚洲十日本专区| 亚洲精品影视一区二区三区av| 国产精品98久久久久久宅男小说| 欧美色视频一区免费| 亚洲精品亚洲一区二区| 精品日产1卡2卡| 久久精品国产清高在天天线| 午夜视频国产福利| 欧美成人性av电影在线观看| 久久精品91蜜桃| 狂野欧美白嫩少妇大欣赏| 97人妻精品一区二区三区麻豆| 法律面前人人平等表现在哪些方面| 欧美黄色片欧美黄色片| 国产亚洲精品av在线| 色播亚洲综合网| 动漫黄色视频在线观看| av片东京热男人的天堂| 搡女人真爽免费视频火全软件 | 高清日韩中文字幕在线| 国产av麻豆久久久久久久| 免费观看人在逋| 变态另类成人亚洲欧美熟女| 国产精品美女特级片免费视频播放器| 亚洲精品成人久久久久久| 成人18禁在线播放| 琪琪午夜伦伦电影理论片6080| 国产精品亚洲美女久久久| 亚洲国产精品成人综合色| 精品一区二区三区视频在线观看免费| 99热6这里只有精品| 国产精品香港三级国产av潘金莲| 日韩大尺度精品在线看网址| 叶爱在线成人免费视频播放| 小说图片视频综合网站| 中文字幕精品亚洲无线码一区| 中文字幕人成人乱码亚洲影| xxx96com| 欧美性感艳星| 亚洲人成网站在线播放欧美日韩| x7x7x7水蜜桃| 丰满的人妻完整版| 麻豆一二三区av精品| 国产精品精品国产色婷婷| 网址你懂的国产日韩在线| 精品一区二区三区视频在线 | 真实男女啪啪啪动态图| 听说在线观看完整版免费高清| 国产精品1区2区在线观看.| 1024手机看黄色片| 亚洲成人中文字幕在线播放| 夜夜躁狠狠躁天天躁| 亚洲av一区综合| 又紧又爽又黄一区二区| x7x7x7水蜜桃| 国产精品久久久久久精品电影| 国产激情偷乱视频一区二区| 亚洲欧美日韩卡通动漫| 日韩欧美 国产精品| 欧美一级毛片孕妇| 天堂√8在线中文| 国产探花极品一区二区| 国产精品嫩草影院av在线观看 | 国产免费男女视频| 久久久久九九精品影院| 午夜影院日韩av| 国产 一区 欧美 日韩| 99国产综合亚洲精品| 小说图片视频综合网站| 久久久久性生活片| 午夜免费观看网址| 日韩欧美在线乱码| 国产精品久久视频播放| 久久欧美精品欧美久久欧美| 精品久久久久久久久久久久久| 国产精品爽爽va在线观看网站| 成年版毛片免费区| 无人区码免费观看不卡| 国产精品久久久久久精品电影| 午夜免费男女啪啪视频观看 | ponron亚洲| 91麻豆av在线| 有码 亚洲区| 老司机在亚洲福利影院| 51国产日韩欧美| 最新中文字幕久久久久| 日日夜夜操网爽| 免费人成视频x8x8入口观看| 人妻夜夜爽99麻豆av| 亚洲av不卡在线观看| 亚洲熟妇熟女久久| 老司机在亚洲福利影院| 国产中年淑女户外野战色| 国产欧美日韩一区二区三| 欧美日本视频| 日本一二三区视频观看| 97人妻精品一区二区三区麻豆| 怎么达到女性高潮| 亚洲国产精品成人综合色| 大型黄色视频在线免费观看| 国产av一区在线观看免费| 亚洲av日韩精品久久久久久密| av福利片在线观看| e午夜精品久久久久久久| 久久久久久久久久黄片| 熟妇人妻久久中文字幕3abv| 亚洲av免费高清在线观看| 啦啦啦免费观看视频1| 国产欧美日韩精品一区二区| 日本 欧美在线| 丁香欧美五月| 免费一级毛片在线播放高清视频| 男女视频在线观看网站免费| 97碰自拍视频| 久久国产精品人妻蜜桃| 狠狠狠狠99中文字幕| 每晚都被弄得嗷嗷叫到高潮| 成人亚洲精品av一区二区| 欧美一区二区国产精品久久精品| 国产黄片美女视频| xxxwww97欧美| www.999成人在线观看| 欧美激情久久久久久爽电影| 操出白浆在线播放| 欧美午夜高清在线| 搞女人的毛片| 熟女人妻精品中文字幕| 久久精品国产清高在天天线| 搡老熟女国产l中国老女人| 性色avwww在线观看| 丁香欧美五月| 看片在线看免费视频| 十八禁人妻一区二区| 亚洲成人中文字幕在线播放| 天堂av国产一区二区熟女人妻| 日韩欧美精品免费久久 | 狠狠狠狠99中文字幕| 日本免费一区二区三区高清不卡| 国产精品永久免费网站| 制服丝袜大香蕉在线| 久久久久久九九精品二区国产| 丰满人妻熟妇乱又伦精品不卡| 国内毛片毛片毛片毛片毛片| 我的老师免费观看完整版| 内地一区二区视频在线| 欧美性感艳星| 91久久精品电影网| 一卡2卡三卡四卡精品乱码亚洲| 成人特级av手机在线观看| 一区二区三区激情视频| 欧美极品一区二区三区四区| 国产精品爽爽va在线观看网站| 精品一区二区三区视频在线观看免费| 熟女少妇亚洲综合色aaa.| 听说在线观看完整版免费高清| 国产一级毛片七仙女欲春2| 国产成人啪精品午夜网站| 国产精品女同一区二区软件 | 国产亚洲精品久久久久久毛片|