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

    A dual-route optical emission spectroscopy diagnostic with wide spectral range and high wavelength resolution on HL-2A tokamak

    2024-04-06 07:15:56YueCHEN陳越JikunGAO高繼昆TingLONG龍婷LinNIE聶林JinmingGAO高金明YaoMA馬堯YuanHUANG黃淵WenjingTIAN田文靜YanminLIU劉延民XiaodongZHU朱曉東GeZHUANG莊革WulyuZHONG鐘武律andMinXU許敏
    Plasma Science and Technology 2024年3期

    Yue CHEN (陳越) ,Jikun GAO (高繼昆) ,Ting LONG (龍婷) ,Lin NIE (聶林) ,Jinming GAO (高金明) ,Yao MA (馬堯) ,Yuan HUANG (黃淵),Wenjing TIAN (田文靜),3,Yanmin LIU (劉延民),Xiaodong ZHU (朱曉東),*,Ge ZHUANG (莊革),Wulyu ZHONG (鐘武律) and Min XU (許敏)

    1 School of Physical Sciences,School of Nuclear Sciences and Technology,University of Science and Technology of China,Hefei 230026,People’s Republic of China

    2 Southwestern Institute of Physics,Chengdu 610041,People’s Republic of China

    3 Tsinghua University,Beijing 100084,People’s Republic of China

    Abstract A dual-route optical emission spectroscopy (D-OES) diagnostic is newly developed to monitor the optical emission from the X-point plasma region on the HL-2A tokamak.This diagnostic is composed of an imaging system,a beam-splitting system for dual-route measurements,fiber bundles,a spectrometer system,and a control and acquisition system.One route is used to obtain wide-spectral-range spectra,and the other route is used to acquire high-wavelengthresolution line shapes.The spectral resolution of the wide-range spectrometers is 0.8 nm with a coverage of 800 nm (@200-1000 nm).The spectral resolution of the high-resolution spectrometer is 0.01 nm with a coverage of 6 nm (@200-660 nm).The spatial resolution of each route of D-OES is about 4 cm with 11 channels.The temporal resolution is 16 ms at maximum in the single-channel mode.Wide-range spectra (containing Balmer series and a Fulcher band)and highly resolved Hα line shapes are obtained by D-OES in the hydrogen glow discharge in the lab.D-OES measurements are carried out in the high-density deuterium experiments of HL-2A.The electron density ne and deuterium temperature TD in the X-point multifaceted asymmetric radiation from the edge (MARFE) region are derived simultaneously by fitting the measured Dα shape.The density ne is observed to increase from~ 8.7×1018 m-3 to~7.8×1019 m-3,and the temperature TD drops from~14.4 eV to~2.3 eV after the onset of MARFE in the discharge #38260.

    Keywords: optical emission spectroscopy,Balmer series,tokamak

    1.Introduction

    The atomic and molecular kinetics and radiation behavior at the boundary of tokamaks strongly affect the plasma-wall interaction [1],divertor detachment and impurity flux into the main plasma [2],which are critical issues in magnetic fusion research.Optical emission spectroscopy (OES)provides a useful diagnostic tool for understanding atomic and molecular processes via the measurements of the optical spectra from the emission of atoms,molecules and ions [3-8].Comprehensive information can be derived from these spectra.The species can be distinguished from the characteristic lines,and their concentration can be obtained from the radiation intensity.Meanwhile,plasma parameters are related to the shapes of spectral lines [9,10].A spectral line with an infinitely sharp line shape cannot be observed due to the natural broadening from Heisenberg’s uncertainty principle.In addition,Doppler broadening from random thermal motion and Stark broadening from the effects of collision and the electric field are important and often not negligible[3,11-13].By fitting the spectral line shape,the plasma temperature from the Doppler width and the electron density from the Stark width can be derived.

    Spectroscopic investigations of plasma in the divertor and edge regions have been carried out on many tokamaks[14-20].Maoet almonitored the line radiation of W I to investigate the tungsten flux in the divertor of EAST [21].By measuring the intensity ratio of C III and C IV,Koubitiet alobtained the electron temperature profile and showed the existence of a recombining zone around the X-point multifaceted asymmetric radiation from the edge (MARFE)in JT-60U [22].Further,the high-resolution spectrometer of Alcator C-Mod was used to obtain the electron density by measuring the line shapes of the deuterium Balmer series[2].Xuet alobtained the electron density in the divertor region by measuring the line width of Dδon EAST [23].The OES diagnostic at the boundary of tokamaks requires the measurement capability of both wide spectral range and high wavelength resolution,because of the existence of different radiating species and various plasma parameters.The timeresolved OES integrated with wide-spectral-range and highwavelength-resolution measurements for the same detecting area is not routinely applied,since the higher resolving power of dispersing elements for spectroscopic techniques usually allows a lower wavelength coverage.

    In this work,we develop a dual-route optical emission spectroscopy (D-OES) diagnostic on the HL-2A tokamak to monitor the optical emission from the X-point region.The dual routes,one for wide-range spectra and the other for highresolution line shapes,are synthesized via a dedicated beamsplitting system.The remainder of this paper is organized as follows.Section 2 introduces the overall design and test results of D-OES.Section 3 reports the derivation of density and temperature from D-OES measurements on the HL-2A.Section 4 presents the conclusion.

    2.Design and test of D-OES diagnostic

    2.1. Layout of D-OES on HL-2A

    A dual-route optical emission spectroscopy (D-OES) diagnostic with wide spectral range and high spectral resolution is newly developed on the HL-2A,which is used to detect the light emitted from the X-point region.The HL-2A is a medium-size tokamak device with a major radius ofR0=1.65 m and a minor radius ofa=0.4 m [24].The DOES is installed at the SSW (i.e.between the south and southwest) port of the HL-2A,which is 30 cm above the midplane.The 11 lines of sight of D-OES in the cross-section of the HL-2A are shown in figure 1.

    Figure 1.The layout of D-OES and some other diagnostics mapped on the cross-section of HL-2A tokamak.FIR: far-infrared laser interferometer;Bol: bolometer.

    2.2. Design and construction of D-OES

    A schematic diagram of the D-OES diagnostic is shown in figure 2 .The D-OES diagnostic consists of an imaging system,a beam-splitting system,fiber bundles,a spectrometer system,and a control and acquisition system.The imaging system contains a shutter,a protecting plate,a reflecting mirror,a vacuum window and imaging lenses.The shutter is in front of a protecting plate.The plate is made of fused silica with a 50 mm diameter,which is installed inside the vacuum chamber of HL-2A.It can avoid contamination by plasma emission.The 75 mm mirror next to the protecting plate is used to change the direction of the optical axis to be coaxial with the subsequent optical components.The window on the vacuum chamber is made of fused silica with a 50 mm diameter.The imaging lenses coaxial with the vacuum window have a focal length of 50 mm.When the shutter is opened,the plasma emission from the object plane near the X-point passes through the protecting plate,gets reflected by the mirror,and then goes through the vacuum window.Later,the emission is collected by the imaging lenses and focused on the image planes.The field of view of the imaging system is 16.8 degrees.The magnification of the lenses is about 40.The diameter of the object plane is about 50 cm,and the diameter of the image plane is 12.5 mm.

    Figure 2.Overall sketch map of D-OES diagnostic.The insert is the enlarged view of fiber patch panel.

    To achieve dual-route measurements,a beam-splitting system after the imaging lenses is used to divide the plasma emission into two beams.The splitting system consists of a non-polarizing 50:50 beam-splitter cube,two filter components and two fiber patch panels,as shown in figure 2.The splitter cube is used to obtain two emission beams.The filter component contains neutral density filters (NDFs) and narrow bandpass filters (NBFs).The NDFs are used to attenuate the emission intensity to avoid overexposure,and the NBFs are used to select the spectral range for measurement.The two fiber patch panels,which are made of aluminum,are installed and aligned at two image planes of the lenses,respectively.Each beam divided by the splitter cube then passes through the filter component.The two filtered beams are transmitted onto the two fiber patch panels and enter the fiber bundles for the wide-range spectrometers and highresolution spectrometer,respectively.There are 5×11 fiber ports on each patch panel for spatial resolved measurements,as shown in the insert of figure 2.The spatial resolution is about 4 cm.The observation spots can be flexibly changed by plugging the fibers into different ports on the patch panel.The 1×11 fiber ports in the middle are used in the experiments.

    The spectrometer system is composed of 11 wide-range spectrometers and one high-resolution spectrometer.The wide-range spectrometers (Thorlabs CCS200) have gratings with 600 grooves/mm,covering a wavelength range of 200-1000 nm.The resolution of the corresponding CCD (Charge Coupled Device) is 1×3648 with an 8μm pitch.The linear CCD can only work in single-channel mode.The highest temporal resolution is 5 ms and the spectral resolution is about 0.8 nm with a fixed entrance slit of 20 μm.The highresolution spectrometer (SOL MSDD1004i) is a double dispersion monochromator-spectrograph with an effective focal length of 1 m.The 2400 grooves/mm grating is used to resolve the emission with a wavelength range of 200-660 nm.The spectral resolution is 0.01 nm and the coverage is 6 nm with an entrance slit of 30 μm.The resolution of the corresponding CCD camera (SOL HS101H) is 506×2048 with a pixel size of 12×12 μm2.In the single-channel mode,this CCD can achieve a 16 ms temporal resolution at maximum.The temporal resolution is about 56 ms in the multi-channel mode.For the synchronous measurements of D-OES,the temporal resolutions of the wide-range spectrometer and high-resolution spectrometer are usually set to be the same.Thus,the temporal resolution of the D-OES is limited by the high-resolution spectrometer.For the whole spectrometer system,a temporal resolution of 16 ms can be achieved in the single-channel mode.

    The control and acquisition system are realized by a local/remote PC.Control software is developed based on the software development kits provided by the manufacturers.Using this software,the control of spectrometers and data acquisition are automatically achieved at the same time.The total D-OES diagnostic can be synchronized with the external trigger signal.

    2.3. Test results of D-OES in the lab

    The D-OES performance is tested with a DC hydrogen glow discharge [25].As shown in the insert of figure 3(a),DC hydrogen plasma is generated between two plate electrodes at 1.3×104Pa with a discharge current of 13 A.The lines of sight of D-OES are set to pass through the middle of the two electrodes.An NBF with a center wavelength of 660 nm and full width at half maximum of 10 nm is used for the highresolution route.The measured optical spectrum in the range of 200-1000 nm is shown in figure 3(a),where the hydrogen Balmer series (Hα656.3 nm,Hβ486.1 nm,Hγ434.0 nm) and hydrogen Fulcher alpha band (590-640 nm)are presented.At the same time,the line shape of Hαis shown in figure 3(b) and the total line width is about 52 pm.

    Figure 3.The spectrum obtained by D-OES in DC hydrogen glow discharge experiment.(a) Spectrum of DC hydrogen glow discharge obtained by wide-range spectrometer,(b) Hα line shape of DC hydrogen glow discharge obtained by high-resolution spectrometer.

    The instrument width of the high-resolution spectrometer is estimated in the lab [26].Since the gas temperature and electron density of the DC hydrogen plasma above are high,the corresponding Doppler width and Stark width are much larger than the instrument width in the total line shape in figure 3(b).Instead,a deuterium electron cyclotron resonance (ECR) plasma at low pressure is used to directly measure the instrument width.The microwave power and plasma pressure are 200 W and 0.1 Pa,respectively.Due to the low electron density and deuterium temperature (near room temperature) of this ECR plasma,the Stark width and Doppler width are below 1 pm.These are much lower than the instrument width (typically dozens of pm).Thus,it is reasonable to regard this measured Dαline width as the instrument width at 656.1 nm.Basically,the instrument broadening is caused by the non-zero slit width and aberrations of the spectrometer [27].The fitted line profile is shown in figure 4 with an entrance slit of 30 μm.The total width is composed of a Gaussian component (9.6±0.3 pm)and a Lorentz component (1.7±0.4 pm) .Therefore,the instrument width is estimated to be about 10 pm with a Gaussian line shape in experiments.

    Figure 4.The measured spectrum and fitted line shape for deuterium ECR plasma.

    3.D-OES measurements on HL-2A

    The D-OES diagnostic is used to monitor the plasma emission in the X-point region of the HL-2A tokamak with divertor configuration.An NBF with a center wavelength of 660 nm and full width at half maximum of 10 nm is used for the high-resolution route to obtain the deuterium Balmer line shape.The Dβand Dγin HL-2A plasmas are extremely weak compared with Dα.Therefore,the high-resolution spectrometer of D-OES is used only to measure the DαBalmer line. Dαconsists of seven fine components distributed within 20 pm.However,their structures are difficult to distinguish because the broadening of their spectral lines (Doppler,Stark,Zeeman and instrumental effects)overlaps with each other.Meanwhile,the measured spectral width on the HL-2A tokamak,which is shown below,exceeds 100 pm.It is much larger than 20 pm.Therefore,the fine structure of Dαis not considered in the OES measurements.

    The main parameters in deuterium Ohmic discharge#38260 are shown in figure 5.The plasma currentIP,the central line-averaged density by the FIR system,the neutral gas pressure measured by a fast ionization gauge,and radiation from bolometer (Bol) measurements are shown in figures 5(a)-(d),respectively.A MARFE characterized by the localized enhanced radiation near the X-point appears at 1246 ms.

    Figure 5.(a) Plasma current,(b) line-averaged density,(c) neutral gas pressure,and (d) bolometer radiation,for discharge #38260.

    Figures 6(a) and (b) display the visible light images.It is found that the light emission near the X-point,which is marked with the dashed box,is greatly enhanced when MARFE is generated.For the D-OES,high-resolution spectrometer measurements are conducted,while the measurements of wide-range spectrometers are not available due to their low sensitivity.An upgrade of the wide-range spectrometers is planned for future work.The Dαline shape measured by the high-resolution spectrometer of D-OES at 1040-1240 ms without MARFE is shown in figure 6(c).The Dαline shape at 1300-1500 ms with MARFE is shown in figure 6(d).

    The Dαline shape is mainly affected by the instrument broadening,Doppler broadening,Stark broadening and Zeeman broadening mechanisms [28].The instrument width ΔλInstrumentis about 10 pm with a Gaussian line shape,as shown in section 2.3.The Doppler line shape originates from the thermal motion of D atoms,and presents a Gaussian line shape.The Doppler width ΔλDis expressed as

    whereMis the relative atomic mass of the D atom,and λ is the central wavelength in nm. ΔλDis only related to the temperatureTDof D atoms.The Zeeman broadening with a Lorentz shape originates from the energy level splitting by external magnetic field with a width of ΔλZas

    whereBis the magnetic field.Ambient charged particles influence the emission process of radiators.This results in a broadening and shift of spectral lines with a Lorentz shape,which is termed Stark broadening [29,30].The Stark width ΔλSis mainly related to the local electron densityne,

    Figure 6.Visible light images and Dα line shape fitting for the cases without (top row) and with (bottom row) MARFE.

    The measured Dαline shape is the convolution of the Lorentz line shape and Gaussian line shape with different functions [28].A nonlinear least-square fitting routine is used to obtain the fitted Lorentz and Gaussian line shape separately.The Doppler width is obtained by extracting the instrument broadening effect from the total broadening ofthe Gaussian line shape,ΔλD=The Stark width is obtained by extracting the Zeeman broadening effect from the total broadening of the Lorentz line shape,ΔλS=ΔλLorentz-ΔλZ.The temperature and density are then estimated from the Doppler width and Stark width,respectively [10,11,30].

    The local electron density and deuterium temperature near the X-point are investigated by analyzing the Dαline shape before and after the MARFE onset.At 1040-1240 ms without MARFE,the fitted Stark width is ΔλS=8 pm with a fitting uncertainty of ±3 pm,and the calculated electron density is 8.7×1018m-3with an uncertainty of±4.1×1018m-3.The fitted Doppler width is ΔλD=126 pm with a fitting uncertainty of ±21 pm,and the calculated deuterium temperature is 14.4 eV with an uncertainty of±4.6 eV.These are shown in figure 6(c).At 1300-1500 ms with MARFE,the Dαintensity is enhanced by about three times,as shown in figure 6(d).The fitted Stark width is ΔλS=33 pm with a fitting uncertainty of ±6 pm,and the calculated electron density is 7.8×1019m-3with an uncertainty of ±2.0×1019m-3.The electron density rises after MARFE appears,and is even higher than the line-averaged density.This feature of high local electron density in the MARFE region was also reported in Alcator C-Mod and EAST [31].Meanwhile,the fitted Doppler width is ΔλD=51 pm with an uncertainty of ±8 pm,and the calculated deuterium temperature drops to 2.3 eV with an uncertainty of ±0.7 eV.The localized high density,low temperature and enhanced radiation of MARFE are in agreement with previous observations on tokamaks [2,23,31,32].

    4.Conclusion

    We have developed a D-OES diagnostic on the HL-2A tokamak.One route is used to measure wide-range spectra,and the other route is used to acquire high-resolution line shapes.The D-OES diagnostic consists of a dedicated imaging system,a beam-splitting system,fiber bundles,a spectrometer system,and a control and acquisition PC.The spectral coverage of the wide-range spectrometers is 200-1000 nm with a resolution of 0.8 nm,while the spectral resolution/coverage of the high-resolution spectrometer is 0.01 nm/6 nm.The D-OES can measure the optical spectra near the Xpoint region of the HL-2A with both the wide spectral range and the high-wavelength resolution.The spatial resolution of each route of D-OES is about 4 cm with 11 channels.The temporal resolution is 16 ms at maximum in the single-channel mode.The D-OES diagnostic is tested using a DC hydrogen glow discharge in the lab.The high-resolution Hαline shape and the wide-range spectra including the hydrogen Balmer series and hydrogen Fulcher alpha band are observed.By conducting D-OES measurements in the highdensity operations of the HL-2A tokamak,the electron density and deuterium temperature in the X-point MARFE region are derived at the same time from the fitted Dαshape,and instrument broadening,Doppler broadening,Stark broadening and Zeeman broadening are taken into account.After the onset of MARFE in discharge #38260,the electron density is observed to increase from~ 8.7×1018m-3to~7.8×1019m-3,and the deuterium temperature drops from~14.4 eV to~2.3 eV.

    Acknowledgments

    The authors thank the HL-2A team for their support with the experiments and discussion.This work was supported by the National MCF Energy R&D Program of China (Nos.2018YFE0301102,2022YFE03100004 and 2018YFE 0303102),National Natural Science Foundation of China(Nos.12375210 and 12305238),and the Sichuan Natural Science Foundation (Nos.2022NSFSC1791,2022JDRC0014 and 2022TFQCCXTD).

    国产精品99久久99久久久不卡| 午夜亚洲福利在线播放| 99在线视频只有这里精品首页| 精品欧美一区二区三区在线| 国产精品野战在线观看| 全区人妻精品视频| 国内揄拍国产精品人妻在线| 狠狠狠狠99中文字幕| 99久久精品国产亚洲精品| 国产精品电影一区二区三区| 日韩免费av在线播放| 欧美黑人精品巨大| 欧美黑人巨大hd| 热99re8久久精品国产| 国产成人av激情在线播放| 小说图片视频综合网站| 99re在线观看精品视频| 国产亚洲精品久久久久5区| 热99re8久久精品国产| 成人一区二区视频在线观看| 999精品在线视频| 免费在线观看完整版高清| 亚洲aⅴ乱码一区二区在线播放 | 色噜噜av男人的天堂激情| 免费看美女性在线毛片视频| 久久草成人影院| 日韩欧美在线乱码| 午夜影院日韩av| 亚洲无线在线观看| 非洲黑人性xxxx精品又粗又长| 中文字幕久久专区| 亚洲 欧美 日韩 在线 免费| 变态另类丝袜制服| 亚洲欧美日韩高清专用| 国产av麻豆久久久久久久| 999久久久精品免费观看国产| 色综合亚洲欧美另类图片| 卡戴珊不雅视频在线播放| 一级av片app| 亚洲成av人片在线播放无| 国内少妇人妻偷人精品xxx网站| 国产一区亚洲一区在线观看| 国产老妇伦熟女老妇高清| 丝袜喷水一区| 婷婷精品国产亚洲av| 男女下面进入的视频免费午夜| 精品免费久久久久久久清纯| 国产精品久久电影中文字幕| 在线天堂最新版资源| 免费av毛片视频| 午夜福利在线观看吧| 高清毛片免费看| 一个人看的www免费观看视频| 久久精品国产自在天天线| 天美传媒精品一区二区| 久久精品国产99精品国产亚洲性色| 日韩成人伦理影院| 日韩精品有码人妻一区| 一级毛片我不卡| 99久久成人亚洲精品观看| 欧美日韩综合久久久久久| 九九久久精品国产亚洲av麻豆| 一区二区三区四区激情视频 | 国产精品人妻久久久影院| 久久久久九九精品影院| 国产 一区 欧美 日韩| 日本成人三级电影网站| 亚洲综合色惰| 1024手机看黄色片| 亚洲av中文字字幕乱码综合| 国产 一区精品| 日韩成人av中文字幕在线观看| 久久久国产成人精品二区| 麻豆成人av视频| 人人妻人人澡欧美一区二区| 久久6这里有精品| 女的被弄到高潮叫床怎么办| 国产精品一区二区在线观看99 | av福利片在线观看| av又黄又爽大尺度在线免费看 | 99热这里只有是精品在线观看| 在线播放无遮挡| 中文欧美无线码| 超碰av人人做人人爽久久| 蜜桃亚洲精品一区二区三区| 在线观看66精品国产| 免费电影在线观看免费观看| 亚洲四区av| 成人综合一区亚洲| 欧美一区二区精品小视频在线| 国产一区二区激情短视频| 午夜激情福利司机影院| 日韩,欧美,国产一区二区三区 | 欧美成人a在线观看| 亚洲精品乱码久久久v下载方式| 老熟妇乱子伦视频在线观看| 亚洲av免费在线观看| 麻豆成人午夜福利视频| 日韩国内少妇激情av| 日韩中字成人| 欧美xxxx黑人xx丫x性爽| 国产色爽女视频免费观看| 国内精品美女久久久久久| 少妇的逼好多水| 亚洲成a人片在线一区二区| 欧美人与善性xxx| 26uuu在线亚洲综合色| 国产精品久久久久久av不卡| 日韩视频在线欧美| 少妇人妻一区二区三区视频| 禁无遮挡网站| kizo精华| 欧美bdsm另类| 精品一区二区三区视频在线| 97超碰精品成人国产| 91av网一区二区| 色播亚洲综合网| 老师上课跳d突然被开到最大视频| 寂寞人妻少妇视频99o| 亚洲内射少妇av| 精品无人区乱码1区二区| 亚洲在久久综合| 亚洲成人中文字幕在线播放| 99久国产av精品国产电影| 久久这里有精品视频免费| 男人舔女人下体高潮全视频| 日韩高清综合在线| 亚洲欧美成人精品一区二区| 亚洲美女视频黄频| 国产亚洲精品久久久com| 免费人成在线观看视频色| av天堂在线播放| 亚洲av一区综合| 99热这里只有精品一区| 两个人的视频大全免费| 国产精品爽爽va在线观看网站| 蜜桃亚洲精品一区二区三区| 三级国产精品欧美在线观看| 免费观看人在逋| 内地一区二区视频在线| 久久久久久九九精品二区国产| 久久久精品欧美日韩精品| 中文字幕制服av| 插逼视频在线观看| 91久久精品国产一区二区三区| 女人被狂操c到高潮| 国产精品国产高清国产av| 高清日韩中文字幕在线| 欧美3d第一页| 我的老师免费观看完整版| 久久久精品欧美日韩精品| 中文字幕制服av| 亚洲国产高清在线一区二区三| 亚洲欧洲日产国产| 亚洲精品456在线播放app| 波野结衣二区三区在线| 成人特级av手机在线观看| 性欧美人与动物交配| 人人妻人人澡人人爽人人夜夜 | 99久久久亚洲精品蜜臀av| 综合色av麻豆| 亚洲精品乱码久久久v下载方式| 永久网站在线| 69人妻影院| 乱码一卡2卡4卡精品| 成人毛片a级毛片在线播放| 国产亚洲91精品色在线| 久久午夜福利片| 欧美性猛交╳xxx乱大交人| 国产精品人妻久久久影院| 午夜福利视频1000在线观看| 色尼玛亚洲综合影院| 18+在线观看网站| 人人妻人人澡欧美一区二区| 成人亚洲欧美一区二区av| av天堂中文字幕网| 亚洲欧美成人综合另类久久久 | 91av网一区二区| videossex国产| 看片在线看免费视频| 国产精品无大码| 国模一区二区三区四区视频| 91精品一卡2卡3卡4卡| 男人舔女人下体高潮全视频| 国产麻豆成人av免费视频| www.av在线官网国产| 免费不卡的大黄色大毛片视频在线观看 | 日韩中字成人| 好男人在线观看高清免费视频| 日日摸夜夜添夜夜添av毛片| 国产一区亚洲一区在线观看| 国产精品久久电影中文字幕| 男人舔奶头视频| 国产精品久久久久久av不卡| 久久久久久久久久成人| 日韩欧美三级三区| 午夜精品国产一区二区电影 | 亚洲精品影视一区二区三区av| 亚洲丝袜综合中文字幕| 色噜噜av男人的天堂激情| 久久久国产成人免费| 日韩,欧美,国产一区二区三区 | 老司机影院成人| 国产伦在线观看视频一区| 久久精品国产清高在天天线| 嘟嘟电影网在线观看| 日本黄色视频三级网站网址| 精品免费久久久久久久清纯| 国产精品野战在线观看| 丰满人妻一区二区三区视频av| 午夜激情福利司机影院| av在线蜜桃| 一夜夜www| 久久久久久久亚洲中文字幕| 男人舔女人下体高潮全视频| 成年免费大片在线观看| 欧美一区二区精品小视频在线| 亚洲国产精品成人综合色| 国产成人精品一,二区 | 青春草视频在线免费观看| 日本色播在线视频| 99热全是精品| 国产女主播在线喷水免费视频网站 | 美女大奶头视频| 99久国产av精品国产电影| 两性午夜刺激爽爽歪歪视频在线观看| 哪里可以看免费的av片| 级片在线观看| 精华霜和精华液先用哪个| 在现免费观看毛片| 在线a可以看的网站| or卡值多少钱| 男女那种视频在线观看| 久久99热这里只有精品18| 成人性生交大片免费视频hd| 国产老妇女一区| 99久国产av精品国产电影| 国产大屁股一区二区在线视频| 亚洲精品久久久久久婷婷小说 | 麻豆乱淫一区二区| 欧美精品一区二区大全| 免费观看的影片在线观看| 成人永久免费在线观看视频| av又黄又爽大尺度在线免费看 | 日韩大尺度精品在线看网址| 国产精品,欧美在线| 美女高潮的动态| 永久网站在线| 欧美在线一区亚洲| 国产成人一区二区在线| 天堂√8在线中文| 国产v大片淫在线免费观看| 成人毛片a级毛片在线播放| 六月丁香七月| 日韩成人伦理影院| 一级毛片电影观看 | 免费一级毛片在线播放高清视频| 国产一区亚洲一区在线观看| 日韩欧美国产在线观看| 亚洲欧美精品专区久久| 少妇熟女欧美另类| 国产精品综合久久久久久久免费| 男的添女的下面高潮视频| 变态另类丝袜制服| 最近最新中文字幕大全电影3| av在线观看视频网站免费| 亚洲欧美清纯卡通| 婷婷精品国产亚洲av| 免费大片18禁| 国产激情偷乱视频一区二区| 亚洲自偷自拍三级| 又爽又黄a免费视频| 欧美日本亚洲视频在线播放| 麻豆国产av国片精品| 精品少妇黑人巨大在线播放 | 国产精品久久久久久精品电影小说 | 亚洲国产高清在线一区二区三| av在线播放精品| 小蜜桃在线观看免费完整版高清| 91麻豆精品激情在线观看国产| 晚上一个人看的免费电影| 22中文网久久字幕| 在线免费十八禁| av福利片在线观看| 日日摸夜夜添夜夜爱| 日韩成人伦理影院| 激情 狠狠 欧美| 国产成人精品一,二区 | 色5月婷婷丁香| 变态另类丝袜制服| 联通29元200g的流量卡| 亚洲av成人av| 国产色婷婷99| 久久久久久久久久黄片| 一本精品99久久精品77| 日日撸夜夜添| 国产乱人偷精品视频| 国产 一区精品| 非洲黑人性xxxx精品又粗又长| 永久网站在线| 男人和女人高潮做爰伦理| 夫妻性生交免费视频一级片| 亚洲国产精品成人综合色| 色5月婷婷丁香| 色综合站精品国产| av专区在线播放| 99久久中文字幕三级久久日本| 校园春色视频在线观看| 26uuu在线亚洲综合色| 色视频www国产| 日韩av不卡免费在线播放| 在现免费观看毛片| 色哟哟哟哟哟哟| 午夜免费激情av| 欧美不卡视频在线免费观看| 日日干狠狠操夜夜爽| 91狼人影院| 成人特级av手机在线观看| 亚洲三级黄色毛片| 亚洲国产欧洲综合997久久,| 日韩视频在线欧美| 哪里可以看免费的av片| 啦啦啦韩国在线观看视频| 色吧在线观看| 中文在线观看免费www的网站| 少妇被粗大猛烈的视频| 晚上一个人看的免费电影| 国产三级在线视频| 女人被狂操c到高潮| 日产精品乱码卡一卡2卡三| 最近2019中文字幕mv第一页| 22中文网久久字幕| 午夜激情欧美在线| 亚洲精品久久久久久婷婷小说 | 亚洲久久久久久中文字幕| 精品人妻视频免费看| 成人性生交大片免费视频hd| 婷婷亚洲欧美| 国产三级中文精品| 国产精品1区2区在线观看.| 丰满的人妻完整版| 好男人视频免费观看在线| 非洲黑人性xxxx精品又粗又长| a级一级毛片免费在线观看| 少妇人妻一区二区三区视频| 在线观看午夜福利视频| 欧美日本视频| 日韩欧美一区二区三区在线观看| 能在线免费观看的黄片| 精品少妇黑人巨大在线播放 | 美女大奶头视频| 午夜福利视频1000在线观看| 老师上课跳d突然被开到最大视频| 此物有八面人人有两片| 性色avwww在线观看| 在线播放无遮挡| 99热只有精品国产| 国产一级毛片在线| 美女脱内裤让男人舔精品视频 | 色播亚洲综合网| 国产 一区 欧美 日韩| 久久久久久九九精品二区国产| 黄色视频,在线免费观看| 观看免费一级毛片| 日本与韩国留学比较| 亚洲欧美清纯卡通| 国产成人一区二区在线| 国产精品免费一区二区三区在线| 日本黄色片子视频| 熟女电影av网| 噜噜噜噜噜久久久久久91| 直男gayav资源| 国产av在哪里看| 男插女下体视频免费在线播放| 国产伦在线观看视频一区| 深夜精品福利| 午夜精品国产一区二区电影 | 九九久久精品国产亚洲av麻豆| 亚洲人成网站在线播| 亚洲欧洲国产日韩| 老女人水多毛片| 亚州av有码| 欧美一级a爱片免费观看看| 欧美精品国产亚洲| 国产在线男女| 成人漫画全彩无遮挡| 亚洲,欧美,日韩| 最近中文字幕高清免费大全6| 成人高潮视频无遮挡免费网站| 亚洲av熟女| 看十八女毛片水多多多| 亚洲精品国产av成人精品| 亚洲第一区二区三区不卡| 好男人在线观看高清免费视频| 99久久久亚洲精品蜜臀av| 乱系列少妇在线播放| 日韩av不卡免费在线播放| 少妇熟女欧美另类| 亚洲成人中文字幕在线播放| 国内久久婷婷六月综合欲色啪| www.色视频.com| 国产色爽女视频免费观看| 中文精品一卡2卡3卡4更新| 国产成人精品一,二区 | 免费无遮挡裸体视频| 亚洲婷婷狠狠爱综合网| 国产成人精品久久久久久| 国产精品久久久久久精品电影小说 | 亚洲精品自拍成人| 成人二区视频| 久久精品国产亚洲av香蕉五月| 中文在线观看免费www的网站| 99热只有精品国产| 成人午夜精彩视频在线观看| 欧美日韩乱码在线| 国产精品1区2区在线观看.| 国产成人影院久久av| 国内精品美女久久久久久| 乱系列少妇在线播放| 久久久午夜欧美精品| 国产亚洲欧美98| 欧美日韩精品成人综合77777| 亚洲不卡免费看| 免费av毛片视频| 国产成人a区在线观看| 床上黄色一级片| 午夜免费男女啪啪视频观看| 搡老妇女老女人老熟妇| 日本一本二区三区精品| 久久久精品94久久精品| 一个人看视频在线观看www免费| 免费人成视频x8x8入口观看| 久久精品91蜜桃| 乱系列少妇在线播放| 亚洲精品亚洲一区二区| 中出人妻视频一区二区| 久久精品国产清高在天天线| 日韩一区二区视频免费看| 最近视频中文字幕2019在线8| 麻豆成人av视频| 免费搜索国产男女视频| 身体一侧抽搐| 色视频www国产| 色噜噜av男人的天堂激情| 日韩欧美国产在线观看| 99久久中文字幕三级久久日本| 观看免费一级毛片| 久久久久久久午夜电影| 特级一级黄色大片| 在线观看午夜福利视频| 欧美一区二区精品小视频在线| 国产精品一二三区在线看| 久久精品国产自在天天线| 国产大屁股一区二区在线视频| 国产成人福利小说| 成人永久免费在线观看视频| 欧美xxxx黑人xx丫x性爽| 人妻夜夜爽99麻豆av| 麻豆乱淫一区二区| 国产高清激情床上av| 亚洲在线自拍视频| 黄色视频,在线免费观看| 99久久九九国产精品国产免费| 最近手机中文字幕大全| 亚洲av中文av极速乱| 国产精品蜜桃在线观看 | 亚洲精品456在线播放app| 日韩制服骚丝袜av| 久久这里有精品视频免费| 少妇裸体淫交视频免费看高清| 免费av不卡在线播放| 精品一区二区三区视频在线| 国产真实伦视频高清在线观看| 国内久久婷婷六月综合欲色啪| 长腿黑丝高跟| 麻豆成人av视频| 欧美xxxx黑人xx丫x性爽| 欧美精品一区二区大全| 欧美成人精品欧美一级黄| 久久中文看片网| 日韩大尺度精品在线看网址| 欧美日韩乱码在线| 全区人妻精品视频| 麻豆乱淫一区二区| 99久久精品热视频| 99久久精品一区二区三区| 国产黄色小视频在线观看| 九九在线视频观看精品| 韩国av在线不卡| 亚洲美女视频黄频| 免费一级毛片在线播放高清视频| 成人午夜高清在线视频| 色哟哟哟哟哟哟| 亚洲一区高清亚洲精品| 久久中文看片网| 18+在线观看网站| 少妇猛男粗大的猛烈进出视频 | 久久精品国产99精品国产亚洲性色| 丰满的人妻完整版| 人人妻人人看人人澡| 久久韩国三级中文字幕| 边亲边吃奶的免费视频| 亚洲av不卡在线观看| 国产男人的电影天堂91| 不卡视频在线观看欧美| 国产精品嫩草影院av在线观看| 中文欧美无线码| 国产又黄又爽又无遮挡在线| 最近视频中文字幕2019在线8| 免费av不卡在线播放| 国产成人aa在线观看| www日本黄色视频网| 成人美女网站在线观看视频| 青青草视频在线视频观看| 在线观看免费视频日本深夜| 免费在线观看成人毛片| 日韩制服骚丝袜av| 久久久欧美国产精品| 性插视频无遮挡在线免费观看| 国产免费一级a男人的天堂| 久久草成人影院| 亚洲自偷自拍三级| 国产成人freesex在线| 午夜福利视频1000在线观看| 亚洲自拍偷在线| 真实男女啪啪啪动态图| 国产精品精品国产色婷婷| 精品不卡国产一区二区三区| a级毛色黄片| 99久久中文字幕三级久久日本| 国产精品不卡视频一区二区| 亚洲va在线va天堂va国产| 人人妻人人澡欧美一区二区| av在线老鸭窝| 免费观看人在逋| 午夜福利在线观看免费完整高清在 | 国产精品一区二区性色av| 亚洲精品456在线播放app| 国产精品人妻久久久影院| 国产探花极品一区二区| 国产片特级美女逼逼视频| av国产免费在线观看| 国产亚洲av片在线观看秒播厂 | 直男gayav资源| 日日撸夜夜添| 国产av麻豆久久久久久久| 日韩强制内射视频| 久久久久国产网址| 亚洲va在线va天堂va国产| 国产一级毛片七仙女欲春2| 久久精品夜色国产| 欧美又色又爽又黄视频| 欧美3d第一页| 看非洲黑人一级黄片| 中文字幕免费在线视频6| 日韩精品有码人妻一区| 乱码一卡2卡4卡精品| 久久人人爽人人爽人人片va| 看黄色毛片网站| 欧美一区二区亚洲| 国产白丝娇喘喷水9色精品| 男人的好看免费观看在线视频| 在线观看av片永久免费下载| 在线免费观看不下载黄p国产| 亚洲欧美日韩无卡精品| 国产片特级美女逼逼视频| 久久热精品热| 成年免费大片在线观看| 99久久久亚洲精品蜜臀av| 国产视频内射| 亚洲欧洲国产日韩| 国产色婷婷99| 我要搜黄色片| 在线观看66精品国产| 精品欧美国产一区二区三| 午夜福利在线在线| 高清毛片免费看| 床上黄色一级片| 成人永久免费在线观看视频| 中国美女看黄片| 日韩欧美 国产精品| 国产v大片淫在线免费观看| 我的女老师完整版在线观看| 亚洲图色成人| 丰满人妻一区二区三区视频av| 99精品在免费线老司机午夜| 乱人视频在线观看| 国产激情偷乱视频一区二区| 免费大片18禁| 国产高清不卡午夜福利| 亚洲欧美中文字幕日韩二区| 人人妻人人澡欧美一区二区| 成年av动漫网址| 国产亚洲91精品色在线| 久久久久免费精品人妻一区二区| 国产精品久久视频播放| 又粗又爽又猛毛片免费看| 极品教师在线视频| 97超视频在线观看视频| 欧美高清成人免费视频www| 亚洲aⅴ乱码一区二区在线播放| 日本-黄色视频高清免费观看| 能在线免费观看的黄片| 亚洲欧美精品专区久久| 又爽又黄a免费视频| 亚洲国产欧美人成| 婷婷亚洲欧美| 毛片女人毛片| 日本与韩国留学比较| 精品人妻一区二区三区麻豆| 长腿黑丝高跟| 亚洲av中文字字幕乱码综合| 欧美一区二区亚洲| av天堂中文字幕网| 乱系列少妇在线播放|