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

    An investigation of the L-shell x-ray conversion efficiency for laser-irradiated tin foils

    2020-05-06 05:59:12DavidBAILIECormacHYLANDRajSINGHStevenWHITEGianlucaSARRIFrancisKEENANDavidRILEYStevenROSEEdwardHILLFeiluWANG王菲鹿DaweiYUAN袁大偉GangZHAO趙剛HuigangWEI魏會岡BoHAN韓波BaoqiangZHU朱寶強JianqiangZHU朱健強andPengqianYANG楊朋千
    Plasma Science and Technology 2020年4期
    關(guān)鍵詞:趙剛王菲大偉

    David BAILIE,Cormac HYLAND,Raj L SINGH,Steven WHITE,Gianluca SARRI,Francis P KEENAN,David RILEY,Steven J ROSE,Edward G HILL,Feilu WANG (王菲鹿),Dawei YUAN (袁大偉),Gang ZHAO (趙剛),Huigang WEI (魏會岡),Bo HAN (韓波),Baoqiang ZHU (朱寶強),Jianqiang ZHU (朱健強) and Pengqian YANG (楊朋千)

    1 Centre for Plasma Physics,School of Mathematics and Physics,Queen’s University Belfast,Belfast,BT7 1NN,United Kingdom

    2 Astrophysics Research Centre,School of Mathematics and Physics,Queen’s University Belfast,Belfast,BT7 1NN,United Kingdom

    3 Plasma Physics Group,Blackett Laboratory,Prince Consort Road,London,SW7 2AZ,United Kingdom

    4 Key Laboratory of Optical Astronomy,National Astronomical Observatories,Chinese Academy of Sciences,Beijing 100012,People’s Republic of China

    5 Department of Astronomy,Beijing Normal University,Beijing 100875,People’s Republic of China

    6 Shanghai Institute of Optics and Fine Mechanics,Chinese Academy of Sciences,Shanghai 201800,People’s Republic of China

    Abstract

    Keywords: conversion efficiency,laser plasma x-ray sources,laboratory astrophysics

    1.Introduction

    Laser plasma x-ray sources have found widespread application in areas such as laboratory astrophysics,inertial fusion physics and warm dense matter research.They have been used as probes,as in x-ray Thomson scattering [1,2] and absorption spectroscopy[3]or as heating sources as in inertial confinement fusion capsules [4],shock drive experiments [5] and laboratory astrophysics experiments,including to create a highly photoionised plasma [6].Many experiments have been performed to date with the aim of producing multi-keV x-ray sources.These experiments include work on single line or line group sources such as K-shell emission lines [7–10].They also include work on multiple line group emission from L-shell transitions[11,12]and broad band M-shell emission [13].A wide variety of elements have been investigated,depending on whether an intense narrow-line source or broad band quasi-continuous source is required.The duration of the x-ray emission,the conversion efficiency of laser energy into x-rays,and the dependence of conversion efficiency on foil thickness have all been investigated.The initial conversion efficiency of laser light into primary x-ray emission is very important in the aforementioned applications.In indirect drive ICF,the compression efficiency will depend on high conversion of the lasers to a broad spectrum of sub-keV x-rays within a Hohlraum[14].For x-ray Thomson scattering applications (e.g.García Saiz et al [15]) the cross sections are low(of order 10-25cm2)and so a high efficiency of narrow bandkeV K-shell line radiation is desirable.For some applications we wish to irradiate a sample with keV x-rays to achieve uniform heating at high density(e.g.Glenzer et al[16]).In this case L-shell emission has proved useful.In addition to this laboratory astrophysics experiments where we wish to preferentially photo-ionise inner shell electrons can make good use of L-shell sources,filtered to remove softer x-rays (e.g.White et al[17]).The more intense we can make such sources,the closer we can come to typical astrophysical photo-ionisation parameters.

    Figure 1.A schematic of the experimental setup,with the Sn foil target placed at the centre of the spherical chamber,(a)looking from the top and(b)from the side.The beam numbers visible from the two directions are marked on the figures.Two curved quartz crystal spectrometers with image plates are on the front and rear sides of the target,at angles of 15°to the horizontal axis of the experiment(indicated by the line WE).Two flat Si crystal spectrometers with CCDs are placed at variable angles in the horizontal plane.

    Figure 2.Typical pulse shape for the nominal 1 ns flat topped pulse.

    Much of the above research has been focused on studies of bulk targets where the emission is measured from the laser irradiated side of the sample.However,in some applications,the experimental geometry may require the utilisation of the x-ray emission out through the non-irradiated side of a thin foil[12,17].This is particularly the case where we are employing the x-ray flux to heat a sample which must then be placed close to the source.In this paper we discuss an experimental setup where we take measurements on both sides of a thin foil target,but concentrate on the non-laser irradiated side.

    Our paper is organised as follows.Details of the experimental setup including diagnostics are provided in section 2,while the data analysis is covered in section 3.Section 4 contains the results and discussion,and finally in section 5 we provide some suggestions for future experimental campaigns,based on our findings.

    2.Experimental setup

    The experiment was conducted at the Shenguang II laser facility [18,19] in November and December 2017.A schematic of our experimental setup is shown in figure 1.A total of 4 overlapping laser beams was used to deliver up to 1050 J of energy at a wavelength of 351 nm on the target.The foils were Sn of thickness 400–800 nm coated onto CH layers of 2–30 μm for mechanical support.Specifically,in our experiment we have investigated the L-shell emission(~3.3–4.4 keV),for irradiation with 3ω(351 nm)laser pulses with both 1 and 2 ns FHWM duration at intensities of~(0.6–3) × 1015W cm-2.In figure 2 we show a typical pulse shape for the 1 ns pulse.Targets were placed at the centre of the spherical chamber,and the diameter of the focal spot of the laser on the target foil varied between 200 and 500 μm,with a flat-top temporal profile.This was the main mechanism for varying the intensity,although we also varied the incident energy for a fixed focus for some of the data,as discussed below.

    L-shell emission spectra between 3.3 and 4.4 keV were recorded using two convex curved quartz (10–10) crystal spectrometers with the signals recorded on image plate.These spectrometers were both placed to view the target from 18.8° above the horizontal.One was set to view the front (laser irradiated) side of the foil at 15° from the horizontal axis of the experiment (indicated by the line WE in

    Figure 3.Sample spectra showing the Sn L-shell emission obtained by(a)the flat Si crystal spectrometer coupled to a CCD for both 30°and 80°viewing,(b)the curved quartz crystal spectrometer coupled to image plate for the same shot as shown in(a)for 30°viewing.The spectra are all corrected for filtering which consisted of 40 μm Ti on both curved spectrometers and 51.5 ± 0.5 μm Al,50 μm Kapton,25 μm Be on the flat Si crystal spectrometers.This error in the Al on the flat Si crystal spectrometers provides the key error in the data and is folded into the error bars of the analysis.

    In figure 3 we plot typical spectra taken with both the CCD,Si crystal and image plate detector based spectrometers,which clearly show the different spectral regions covered.In previous work by White et al [17] the wavelengths of the Sn L-shell spectral features were calculated by comparison with calibration shots using K-shell emission from Cl and K targets,and we used these for the features detected in the present data.The sensitivity of the CCD based spectrometers is relatively straightforward to calibrate,as the detection efficiency of the CCD is known and can be tested with the use of a calibrated55Fe source.This also gives the A–D conversion,i.e.the number of counts per photon absorbed.The crystal reflectivity can be calibrated by the procedure described in Valle Brozas et al [21],and values similar to the calculated ones for Si are obtained.Cross calibration was then possible by placing the CCD based spectrometers at 30° from the target normal,very similar to the angle of the fixed spectrometers.

    Since we know the distances from source to Si crystal and CCD,as well as the calibrations of the CCD and crystal,we can calculate the number of photons emitted per steradian in the spectral range covered by the line group in figure 3(a)(roughly 3.6–3.85 keV).From figure 3(b),we can determine what fraction,F,of the total L-shell emission is represented by this limited spectral range sample and thus estimate the photons per steradian emitted into the whole L-shell spectrum.This fraction was found to be relatively stable from shot-to-shot (0.51 ± 0.03) and not sensitive to the thickness of the Sn layer or intensity.This latter calculation is,however,only strictly valid for data where the CCD spectrometer(set to 30°) is close to the same angle-of-view as the image plate based spectrometers (with an effective angle of 23.9°,implying the cosine differs only by 5.6%).If not the ratio of intensity between the line groups may vary with angle-ofview due to effects of opacity.figure 1),the other to view the rear side of the foil.In each case the combined horizontal and vertical angles meant the spectrometers viewed the foils at an effective angle of 23.9°from the target normal.The resolution of these spectrometers was dominated by the resolution of the image plate detectors,and values of E/ΔE ~250 were achieved.The type of image plate used was 112 μm thick layer of Fuji BAS-SR (density = 3.07 g cm-3),BaF(Br,I):Eu2+,resolution 109 ± 2 μm [20],with the image plate scanner set to 25 μm resolution.In addition to these fixed spectrometers,we employed two spectrometers with flat Si (111) crystals coupled to x-ray CCDs (Andor DX420 and an Andor DX435).These latter two spectrometers recorded a smaller spectral range of L-shell emission but were able to be moved to different viewing angles in the experiment,although they always viewed the targets in the horizontal plane.The Si crystals and CCD detectors have a known efficiency and data from these were used to calibrate the fixed spectrometers as described below.The resolution of the flat crystal spectrometers was dominated by source size and was typically E/ΔE ~ 400–800.

    3.Data analysis

    For the CCD data our procedure is to correct the measured number of counts (Nccd) to account for the energy dependent transmission of the filters and the quantum efficiency of the CCD,again a function of photon energy.By partially shielding the CCD from direct diffraction,we were able to establish that the background level coming from crystal fluorescence was typically 100 counts per pixel,compared to peaks of over 10 000 counts per pixel for the data.In the x-ray regime we know that 3.65 eV of absorbed energy is needed per electron–hole pair [22],and we have calibrations of the number of counts per photon in the keV photon regime from photon counting measurements with K–α sources.The intensity,I,in photons per steradian is given by:

    Figure 4.X-ray conversion yield plotted as a function of flat spectrometer position angle for an average of focal spot sizes 300,400 and 500 μm at different target thicknesses.The dependence on Sn thickness and intensity is weak (here we use data with intensity spread between 3×1014 and 1.3×1015 W cm-2).See text and further data for details on Sn thickness and intensity.The laser pulse duration was 1 ns.

    where Qccdis the quantum efficiency of the CCD (a function of the photon energy,E).T(E)is the filter transmission(again a function of photon energy,E),γ is the number of counts per photon and Ω the effective solid angle of collection of the spectrometer.This is given by Ω =Rcw Dwhere Rcis the integrated reflectivity of the crystal,w the width of the CCD detector and D the distance from source to detector surface.Having determined the photons emitted per steradian we are then in a position to calculate the conversion efficiency to x-rays in the L-shell emission by comparison with the measured laser energy on target.

    4.Results and discussion

    In figure 4 we plot the x-ray conversion yield,into the line group recorded on the CCD based spectrometer,spanning 3.69–3.85 keV.This is presented as a function of angle,defined as being relative to the normal to the rear side of the foil.As expected,there is fall off in yield with increasing angle.We have fitted the data to a simple cosine dependence that might be expected from a black body emitter.As we can see,the cosine fit is consistent with the average to within the statistical error bars at all measured angles.In calculating the total emission efficiency,for subsequent data,we have used the angular dependence of this latter model as a basis to extrapolate from results measured at a particular angle.

    At 30° we have the highest variability with a spread of roughly a factor of two between the highest and lowest conversion yield and an average of 1.11 ± 0.26 × 1012photons/sr/J where the error bar is the standard deviation,which is not untypical of laser-plasma experiments.If we take the cosine fit as giving the variation with angle,we can use this to estimate that the efficiency of conversion to the line group from 3.69 to 3.85 keV is approximately 0.46% ± 0.11%.The statistical spread dominates over expected estimates of the individual error bars for data,which are in turn dominated by an estimated±10%systematic error in the crystal calibration.If we make the approximation of taking the ratio of the first line group to the total from the image plate data and apply it to all angles,we would estimate a total conversion efficiency of around 0.9% ± 0.2%.

    Figure 5.Variation of the x-ray conversion yield with thickness of the Sn foil for pulse duration 1 ns and flat/curved crystal spectrometers at an angle of ≈30°.The crosses are for spot size 300 μm,the solid circles are for spot size 400 μm and the circles for 500 μm focus.Data is for a single line group centred at 3.69–3.85 keV.

    Figure 6.X-ray conversion yield plotted as a function of intensity,with different focal spots(data taken from CCD based spectrometer).The laser pulse duration was 1 ns and the Sn thickness was 400 nm for all shots.Error bars are from estimated uncertainty of ±10% in laser energy and similar error in filter and crystal efficiency.

    The measured variation of the total conversion yield,on the thickness of the Sn foil is plotted in figure 5.An inspection of the figure indicates that whilst there is shot-toshot variation,there is no clear dependence on the thickness for values in the range of 400–800 nm.The average conversion yield overall,for the 400 μm focus,is 1.15 × 1012photons/sr/J with a standard deviation of about 30%.For the same focus,the individual average conversion efficiencies into the single line group,for the three Sn thicknesses agree to better than 10% with each other.The laser energies for this data set were all within 7% of an average of 907 J.

    By utilising focal spot sizes between 200 and 500 μm we can vary the intensity and our results are shown in figure 6.

    Figure 7.Ratio of L-shell emission for front and rear sides of the foil as a function of irradiance.We see that the ratio is close to unity.

    We can see that the yield quickly rises to a peak at approximately 5 × 1014W cm-2but then falls slowly as the intensity increases.This behaviour is expected for the following reasons.The x-ray emission depends on collisional ionisation and excitation processes that are highly temperature dependent.At low irradiance the plasma has not reached a high enough temperature to efficiently emit in the keV range.At higher intensities,the temperature is high enough but the efficiency of inverse bremsstrahlung absorption drops and this is expected to be the dominant absorption mechanism for our short wavelength,long pulse laser parameters [23].Furthermore,as discussed by Riley et al [9] at higher intensity,a higher fraction of the laser energy absorbed is converted to kinetic rather than thermal energy,leading to a slow decrease in conversion efficiency.

    Finally,since we have fixed spectrometers viewing at 23.8° to both the front (i.e.laser irradiated),and rear sides of the foils,we can determine the ratio of emission from either side of the foils and the results are shown in figure 7.We note that although the ratio is close to unity,the rear(non-irradiated)side emission is typically around 10%–15% weaker.This is a plausible result because as the heat front generated by the laser irradiation starts to burn through the Sn layer,the emission generated to the rear side needs to pass through a cooler absorbing plasma.By contrast,on the front side,the emission passes through a hotter,more ionised plasma with a large velocity gradient,reducing opacity to the L-shell emission.

    A summary of the experiment is that we have measured a conversion efficiency of approximately 0.46% for the line group centred around 3.8 keV.Full spectra suggest the conversion overall into L-shell is of order 1%if the emission ratio between the single line group (3.69–3.85 keV) and multi line groups (3.3–4.4 keV) is constant at different viewing angles.However,the ratio between groups may vary with angle as a result of the changing effective opacity.The irradiation was with 1 ns laser pulses in the third harmonic (351 nm) and results are broadly in line with data from previous experiments[12,17].The angular variation appears to be approximated reasonably well by a simple cosine fit.There is no clear evidence for an optimal thickness for the Sn layer between 400 and 800 nm.Furthermore,we see only a modest difference in the ratio of emission from the front to the rear of the foil.

    Figure 8.HYADES simulation for two cases,one with a 400 nm Sn layer,the other 800 nm.In both cases we show only the Sn plasma and not the CH onto which it is coated,and used 40 cells for the Sn layer.The laser intensity is 5 × 1014 W cm-2,the pulse is modelled as a 1 ns full-width-at-half-maximum (FWHM) flat-top pulse,and results are shown at the centre (1.3 ns on figure 2) of the pulse.

    The reason for these findings can be considered.One possibility is that the laser only burns though around 400 nm or less of Sn,so that adding extra thickness does not increase the emission and the cold layer of additional Sn is only a few hundred nm thick,and does not significantly attenuate the emission to the rear side of the foil.This would explain why there is little or no dependence on the thickness of the initial layer.We have tested this hypothesis with hydro-dynamic simulation using the HYADES radiation-hydrodynamic code[24].This has multi-group radiation transport,and we have used it with 50 photon energy groups spaced logarithmically between 0.03 and 20 keV.We have assumed a laser intensity of 5 × 1014W cm-2and a 1 ns flat-top pulse shape.The Sn layer is modelled with 40 cells and the CH layer with 60 cells.As we can see in the figure below,the expansion is expected to exceed the lateral size of the focal spot.For this reason,we have used a 1D spherical geometry with a radius of curvature of 500 μm,which is similar to the focal spot for most data.This means that close to the surface,we have expansion similar to 1D planar but further out,the density is allowed to fall off which is more in keeping with the actual experiment.

    In figure 8 we show the HYADES profile of the Sn layer at the peak of the incident 3ω laser for starting Sn layers of 400 nm and 800 nm.Since the profiles are generated at the same time (the peak of the pulse),the fact that a higher density,cooler region still exists for the 800 nm layer profile shows that burn-through of the layer is not as complete as for the thinner layer,as expected.The 800 nm foil does burnthrough later in the pulse,for example,at the end of the flattopped portion of the pulse profile,the 800 nm Sn layer is all at above 2 keV.Consequently,a limited burn-through is unlikely to be more than a partial explanation as the lack of variation with layer thickness.

    Figure 9.(a) Simulated spectrum for a 400 nm Sn foil case at the peak of the laser pulse,compared to an experimental spectrum.The FLYCHK simulation is initiated at LTE but allowed to evolve with non-LTE time dependent atomic kinetics.The temperature and density as well as plasma size are taken from the HYADES simulation.(b) Simulated total L-shell emission between 3.3 and 4.4 keV.

    In addition to hydrodynamic simulations we have used the FLYCHK code[25]to model the L-shell emission based on the temporal history of the Sn layer.We do this by dividing the plasma into four zones.The reasons for this are,firstly,since there are gradients of density and temperature within the Sn layer,it is not accurate to model the emission simply using the overall averages.On the other hand,modelling 40 separate plasmas,based on the hydro-code cells has the problem that opacity effects would not then be properly accounted for since FLYCHK is a zero-dimensional code.Additionally,the opacity is likely to be strongly affected by Doppler decoupling effects caused by the strong velocity gradient,which we can see from figure 8 is of the order 109s-1.This means that with a single layer,we would overestimate the opacity.Thus division into four zones that are effectively considered as Doppler decoupled is a compromise,based on the estimated line-widths of the strongest lines,and the velocity gradients,that allows us to still account for,in a coarse manner,the gradients in plasma conditions in determining average ionisation and excitation of the L-shell.

    In figure 9 we show two FLYCHK outputs,the first is a simulated integrated spectrum sample for a 400 nm Sn layer compared to a sample experimental spectrum.The features broadly match those seen experimentally,but there are some mismatches in line positions.This is perhaps not so surprising as,in order to keep databases manageable,for higher Z materials,FLYCHK uses schematic energy levels,generated from hydrogenic levels with screening constants,rather than detailed results of relativistic atomic physics calculations.In figure 9(b)we have shown the simulated temporal emission of the total L-shell between 3.3 and 4.4 keV.As we can see,we predict a duration of emission of about 0.3 ns FWHM.During the experiment we were unable to obtain streak camera data,but note that this emission period is shorter than the approximately 0.5 ns seen in previous work at laser wavelength 532 nm [17]with thinner(217 nm)layers.The simulations predict an overall conversion efficiency into L-shell emission of approximately 7%,which is somewhat higher than our experimental observation.We should however caution that the simulations make some approximations in both the hydrodynamic and spectral simulations as described above.

    As noted above,the reason for the lack of dependence on foil thickness appears to be only partially a result of limited burn-through of the Sn layer.Instead the lack of dependence is perhaps more connected to the fact that,even with Doppler decoupling,there is predicted to be a larger opacity for the Sn L-shell lines,with τ ~ 5-10 for the strongest lines.Thus,we might conclude that,for the range of the layer thickness employed here,we create a plasma slab with broadly similar range of conditions and that the amount of material heated is not a controlling factor since emission from the core of the Sn plasma is suppressed by strong opacity anyway.A cosine dependence,is what we would expect from an optically thick slab close to black body conditions.

    In conclusion,we have shown that L-shell Sn can form the basis of an efficient laser-plasma x-ray source in the 3.3–4.4 keV spectral region.As might be expected,the x-ray conversion efficiency of ~1%is higher than generally seen in single K-shell line sources.Modelling of the source for such a high-Z plasma is non-trivial,with radiation transfer likely to play a significant role in the evolution of the plasma density and temperature.The atomic kinetics of the plasma are also strongly dependent on the opacity of the lines,with Doppler decoupling likely to also play a major role.For this reason we use our simulations to confirm the plausibility of the experimental result rather than being definitively predictive.

    5.Future experimental campaigns

    It is evident that more work could be undertaken to explore a wider range of thicknesses for the Sn foil,in particular,to investigate the total emission as the thickness is decreased.Our previous experiments with thinner foils[17]were carried out at 532 nm and hence are not directly comparable.It would also be advantageous to obtain data on the temporal profile of the emission as a function not only of thickness but also intensity.Since one of our primary purposes is to create a keV x-ray drive for photoionisation experiments,it would be useful to deploy additional diagnostics to measure the softer x-ray component,especially through the rear side of the target when the source x-rays have passed through the thick CH layer.The purpose of this layer is to suppress the softer x-ray component and thus create a source with a higher effective radiation temperature[17,26].By controlling the thickness of CH we can,in principle,control the ratio of outer-shell to inner-shell photoionisation.To do this effectively,an array of calibrated and filtered diodes could be used to measure the emitted spectrum.

    Acknowledgments

    This work was supported by the UK Science and Technology Facilities Council,National Natural Science Foundation of China(No.11573040)and Science Challenge Project(No.TZ2016005)and The Royal Society International Exchange (No.IE161039).The authors gratefully acknowledge the expert support from the VULCAN laser facility and from the Shanghai Institute of Laser Plasma.We also thank the staff of the SG-II laser facility.

    猜你喜歡
    趙剛王菲大偉
    《夏日》
    《秋》
    院前急救標準化管理在急性心肌梗死患者中的應用
    張大偉作品
    Enhanced microwave absorption performance of MOF-derived hollow Zn-Co/C anchored on reduced graphene oxide?
    王菲版《我和我的祖國》為什么也那么動聽?
    綠色中國(2019年19期)2019-11-26 07:13:30
    趙剛作品
    王菲與那英的《歲月》
    海峽姐妹(2018年3期)2018-05-09 08:21:05
    心底最大的隱痛
    神奇的邊界線:一不留神就出國
    智慧少年(2017年8期)2018-01-10 21:39:12
    国产av精品麻豆| 伊人久久大香线蕉亚洲五| 日本wwww免费看| 国产在线视频一区二区| 男女床上黄色一级片免费看| 国产欧美亚洲国产| 国产欧美亚洲国产| 国内毛片毛片毛片毛片毛片| 欧美激情高清一区二区三区| 69av精品久久久久久 | 国产av精品麻豆| 捣出白浆h1v1| 永久免费av网站大全| 美女主播在线视频| 老司机午夜十八禁免费视频| 午夜老司机福利片| 免费av中文字幕在线| 黄色a级毛片大全视频| 午夜福利视频精品| 一区二区三区乱码不卡18| 大片免费播放器 马上看| 无限看片的www在线观看| 手机成人av网站| 一级a爱视频在线免费观看| 九色亚洲精品在线播放| 窝窝影院91人妻| 五月天丁香电影| 国产精品一二三区在线看| 日韩制服丝袜自拍偷拍| 无遮挡黄片免费观看| 亚洲男人天堂网一区| 国产在线免费精品| 亚洲一码二码三码区别大吗| 国产成人精品久久二区二区91| 精品免费久久久久久久清纯 | 精品久久久久久电影网| 十八禁网站网址无遮挡| 日韩精品免费视频一区二区三区| 婷婷色av中文字幕| av一本久久久久| 婷婷丁香在线五月| 日本一区二区免费在线视频| 一区在线观看完整版| 久久这里只有精品19| 国产欧美日韩综合在线一区二区| 国产成人一区二区三区免费视频网站| 大片电影免费在线观看免费| 久久热在线av| 在线观看舔阴道视频| 一边摸一边做爽爽视频免费| 亚洲专区中文字幕在线| 人人妻,人人澡人人爽秒播| 欧美变态另类bdsm刘玥| 亚洲精品美女久久久久99蜜臀| 成在线人永久免费视频| 1024视频免费在线观看| 久久香蕉激情| 成年女人毛片免费观看观看9 | 69av精品久久久久久 | 久久亚洲精品不卡| 女性被躁到高潮视频| 欧美日韩亚洲高清精品| av天堂在线播放| av在线app专区| 午夜影院在线不卡| 国产精品.久久久| 制服诱惑二区| 国产高清国产精品国产三级| 男人爽女人下面视频在线观看| 新久久久久国产一级毛片| 亚洲久久久国产精品| 欧美xxⅹ黑人| 两性夫妻黄色片| 免费av中文字幕在线| 天天躁日日躁夜夜躁夜夜| 欧美亚洲日本最大视频资源| 午夜福利视频在线观看免费| 亚洲久久久国产精品| 91av网站免费观看| 久久人妻福利社区极品人妻图片| 日本wwww免费看| 日日夜夜操网爽| 水蜜桃什么品种好| 亚洲精品一区蜜桃| 欧美黑人精品巨大| 中文字幕精品免费在线观看视频| 亚洲一区二区三区欧美精品| 18在线观看网站| 免费av中文字幕在线| 曰老女人黄片| 亚洲欧美色中文字幕在线| 俄罗斯特黄特色一大片| 三级毛片av免费| 国产精品一区二区免费欧美 | 搡老岳熟女国产| 久久毛片免费看一区二区三区| 精品国产一区二区三区久久久樱花| 18禁黄网站禁片午夜丰满| 国产精品二区激情视频| 亚洲精品美女久久久久99蜜臀| 国产亚洲精品久久久久5区| 1024视频免费在线观看| 亚洲精品日韩在线中文字幕| 999精品在线视频| 男女午夜视频在线观看| 欧美黑人欧美精品刺激| 亚洲精品美女久久久久99蜜臀| 国产男女超爽视频在线观看| 老司机午夜十八禁免费视频| 国产1区2区3区精品| 999久久久精品免费观看国产| 亚洲第一青青草原| 久久久久国内视频| 夫妻午夜视频| 天天躁日日躁夜夜躁夜夜| 欧美亚洲日本最大视频资源| 国产淫语在线视频| 狠狠婷婷综合久久久久久88av| 日本wwww免费看| 亚洲中文日韩欧美视频| 成人av一区二区三区在线看 | 久久人妻熟女aⅴ| 亚洲伊人久久精品综合| 午夜免费鲁丝| 国产成+人综合+亚洲专区| 老司机深夜福利视频在线观看 | 午夜91福利影院| 久久精品亚洲av国产电影网| 国产极品粉嫩免费观看在线| 欧美精品av麻豆av| 五月天丁香电影| 自线自在国产av| 热99久久久久精品小说推荐| 亚洲黑人精品在线| 国产成人欧美在线观看 | 精品人妻在线不人妻| 亚洲久久久国产精品| 一二三四社区在线视频社区8| 亚洲中文字幕日韩| 不卡av一区二区三区| 黄色片一级片一级黄色片| 国产激情久久老熟女| 亚洲精品av麻豆狂野| 黑人巨大精品欧美一区二区mp4| 国产真人三级小视频在线观看| 成人免费观看视频高清| 又大又爽又粗| 两人在一起打扑克的视频| 高清在线国产一区| 欧美精品av麻豆av| 最近最新中文字幕大全免费视频| 美女大奶头黄色视频| 国产成人免费观看mmmm| 亚洲国产精品一区三区| 一二三四在线观看免费中文在| 国产日韩欧美在线精品| 欧美日韩中文字幕国产精品一区二区三区 | 国产熟女午夜一区二区三区| 亚洲欧美成人综合另类久久久| 大码成人一级视频| 热re99久久国产66热| 一区二区三区激情视频| 美女午夜性视频免费| 五月开心婷婷网| 又大又爽又粗| www.精华液| 手机成人av网站| 亚洲 国产 在线| 精品少妇一区二区三区视频日本电影| 婷婷丁香在线五月| 亚洲精品国产色婷婷电影| 咕卡用的链子| 黄色 视频免费看| 99国产精品免费福利视频| 中国美女看黄片| 久久这里只有精品19| 高潮久久久久久久久久久不卡| 国产成人a∨麻豆精品| 一本久久精品| 亚洲精品在线美女| 亚洲精品国产av蜜桃| 久久亚洲精品不卡| 在线亚洲精品国产二区图片欧美| 一区二区日韩欧美中文字幕| 少妇裸体淫交视频免费看高清 | 嫩草影视91久久| 99国产精品免费福利视频| 亚洲精品久久久久久婷婷小说| 亚洲avbb在线观看| 国产精品久久久久成人av| 十八禁高潮呻吟视频| 成人免费观看视频高清| 视频区图区小说| 王馨瑶露胸无遮挡在线观看| videos熟女内射| 亚洲精品日韩在线中文字幕| 国产成人a∨麻豆精品| 久久狼人影院| 久久久精品国产亚洲av高清涩受| 水蜜桃什么品种好| 超色免费av| 久久国产精品人妻蜜桃| 啦啦啦免费观看视频1| 如日韩欧美国产精品一区二区三区| 天天躁日日躁夜夜躁夜夜| 久久天堂一区二区三区四区| 精品国产一区二区久久| 97人妻天天添夜夜摸| 国产三级黄色录像| 久久影院123| 新久久久久国产一级毛片| 国产成人精品久久二区二区91| 午夜免费鲁丝| 欧美午夜高清在线| 国产区一区二久久| 黄色 视频免费看| 国产成人精品久久二区二区91| 国产99久久九九免费精品| 亚洲欧美色中文字幕在线| 国产av又大| 我要看黄色一级片免费的| 热99久久久久精品小说推荐| 欧美一级毛片孕妇| 日韩欧美国产一区二区入口| 久久九九热精品免费| 亚洲欧美清纯卡通| 欧美日韩中文字幕国产精品一区二区三区 | 波多野结衣av一区二区av| 美女脱内裤让男人舔精品视频| 日韩制服骚丝袜av| 亚洲 欧美一区二区三区| 丰满人妻熟妇乱又伦精品不卡| 高清欧美精品videossex| 国产又色又爽无遮挡免| 久久久国产成人免费| 精品久久久精品久久久| 亚洲欧美色中文字幕在线| 大香蕉久久成人网| 99国产精品99久久久久| 国产日韩欧美在线精品| 国产成人精品久久二区二区91| 中文精品一卡2卡3卡4更新| 亚洲精品第二区| 亚洲第一青青草原| 国产免费福利视频在线观看| 一进一出抽搐动态| 色视频在线一区二区三区| 亚洲精品第二区| 亚洲男人天堂网一区| 欧美人与性动交α欧美精品济南到| 精品国内亚洲2022精品成人 | 丝袜美腿诱惑在线| 天堂8中文在线网| 在线观看一区二区三区激情| 亚洲国产av影院在线观看| 日韩一区二区三区影片| 91av网站免费观看| 欧美日韩国产mv在线观看视频| 99精品欧美一区二区三区四区| 老司机福利观看| 老熟女久久久| 免费高清在线观看视频在线观看| 国产精品1区2区在线观看. | 日韩 欧美 亚洲 中文字幕| 中文欧美无线码| 国产一区二区 视频在线| 下体分泌物呈黄色| 国产精品一区二区免费欧美 | 性色av乱码一区二区三区2| av片东京热男人的天堂| 国产精品99久久99久久久不卡| 欧美激情 高清一区二区三区| 国产麻豆69| 99国产精品一区二区三区| 国产又爽黄色视频| 日韩大码丰满熟妇| 国精品久久久久久国模美| 婷婷成人精品国产| 在线天堂中文资源库| 亚洲欧美成人综合另类久久久| 国产在视频线精品| 日韩电影二区| 欧美精品一区二区免费开放| 精品少妇内射三级| 国产在线免费精品| xxxhd国产人妻xxx| 欧美老熟妇乱子伦牲交| 男女高潮啪啪啪动态图| 日韩欧美一区视频在线观看| 色综合欧美亚洲国产小说| 久久精品国产a三级三级三级| 汤姆久久久久久久影院中文字幕| 亚洲欧美精品综合一区二区三区| 免费女性裸体啪啪无遮挡网站| 国产成人a∨麻豆精品| 一本一本久久a久久精品综合妖精| 日本猛色少妇xxxxx猛交久久| 91成人精品电影| 女警被强在线播放| 久久久久久久国产电影| 香蕉国产在线看| 亚洲欧美精品综合一区二区三区| 免费女性裸体啪啪无遮挡网站| 国产成人精品久久二区二区免费| 黄色视频不卡| 精品少妇久久久久久888优播| 午夜福利一区二区在线看| 大码成人一级视频| 亚洲国产毛片av蜜桃av| 99久久精品国产亚洲精品| 老汉色∧v一级毛片| 男女边摸边吃奶| 成年人午夜在线观看视频| 国产精品自产拍在线观看55亚洲 | 欧美中文综合在线视频| 精品视频人人做人人爽| 午夜久久久在线观看| 欧美精品高潮呻吟av久久| 女警被强在线播放| 嫩草影视91久久| 久久香蕉激情| 欧美日韩亚洲高清精品| 欧美日韩福利视频一区二区| 午夜视频精品福利| 久久久国产成人免费| 性色av一级| 视频在线观看一区二区三区| 亚洲欧美精品综合一区二区三区| 中亚洲国语对白在线视频| 手机成人av网站| 亚洲欧美一区二区三区黑人| 免费日韩欧美在线观看| 精品免费久久久久久久清纯 | 久久99热这里只频精品6学生| 91麻豆精品激情在线观看国产 | 国产欧美日韩精品亚洲av| 久久午夜综合久久蜜桃| 美女高潮喷水抽搐中文字幕| 欧美久久黑人一区二区| 在线永久观看黄色视频| www.熟女人妻精品国产| 无遮挡黄片免费观看| 欧美成狂野欧美在线观看| 色婷婷久久久亚洲欧美| 国产在线视频一区二区| 中文字幕色久视频| 老熟妇仑乱视频hdxx| 亚洲精品国产色婷婷电影| 国产精品99久久99久久久不卡| 99精品欧美一区二区三区四区| 少妇人妻久久综合中文| 9热在线视频观看99| 久久久久精品人妻al黑| 成人三级做爰电影| 亚洲精品中文字幕在线视频| 精品国产一区二区三区久久久樱花| www.熟女人妻精品国产| 成人国产一区最新在线观看| 免费女性裸体啪啪无遮挡网站| 9色porny在线观看| 亚洲精品国产区一区二| kizo精华| 国产一区二区在线观看av| 99久久综合免费| 99久久人妻综合| 免费看十八禁软件| 国产亚洲精品一区二区www | a级毛片黄视频| 亚洲精品国产精品久久久不卡| 丰满人妻熟妇乱又伦精品不卡| 免费日韩欧美在线观看| 性色av一级| 久久精品aⅴ一区二区三区四区| 少妇人妻久久综合中文| 悠悠久久av| 一级片免费观看大全| 两性午夜刺激爽爽歪歪视频在线观看 | 国产极品粉嫩免费观看在线| 美女国产高潮福利片在线看| 欧美变态另类bdsm刘玥| 欧美日韩一级在线毛片| 一二三四在线观看免费中文在| 国产一区二区三区综合在线观看| 国产精品免费大片| 亚洲av成人不卡在线观看播放网 | 侵犯人妻中文字幕一二三四区| 久久久国产精品麻豆| 欧美另类一区| av天堂在线播放| 丝袜人妻中文字幕| 搡老熟女国产l中国老女人| 他把我摸到了高潮在线观看 | 国产男人的电影天堂91| 悠悠久久av| av超薄肉色丝袜交足视频| 天天躁日日躁夜夜躁夜夜| 深夜精品福利| 欧美另类亚洲清纯唯美| 99久久人妻综合| 久久九九热精品免费| 亚洲国产精品成人久久小说| videosex国产| 搡老乐熟女国产| 日本欧美视频一区| 99热国产这里只有精品6| 侵犯人妻中文字幕一二三四区| 免费女性裸体啪啪无遮挡网站| 中文字幕av电影在线播放| 久久精品国产亚洲av香蕉五月 | 精品久久久久久久毛片微露脸 | 国产三级黄色录像| 人妻人人澡人人爽人人| 99国产精品一区二区蜜桃av | 不卡av一区二区三区| 永久免费av网站大全| 97人妻天天添夜夜摸| 丝袜人妻中文字幕| 亚洲国产成人一精品久久久| 两人在一起打扑克的视频| 无遮挡黄片免费观看| 成人三级做爰电影| 久久中文看片网| 国产99久久九九免费精品| 久久精品人人爽人人爽视色| 超碰成人久久| 视频在线观看一区二区三区| 欧美 日韩 精品 国产| 久久久精品国产亚洲av高清涩受| 搡老岳熟女国产| 夫妻午夜视频| 十八禁网站网址无遮挡| 丰满少妇做爰视频| 美女大奶头黄色视频| 亚洲熟女毛片儿| 中文字幕人妻熟女乱码| 极品少妇高潮喷水抽搐| 欧美中文综合在线视频| 男女床上黄色一级片免费看| 色播在线永久视频| 国产成人精品无人区| 亚洲五月色婷婷综合| 国产色视频综合| 亚洲专区中文字幕在线| av又黄又爽大尺度在线免费看| 一区福利在线观看| 黄片播放在线免费| 老司机影院毛片| 久久精品亚洲熟妇少妇任你| 国产成人影院久久av| 丝袜在线中文字幕| 亚洲免费av在线视频| 美女午夜性视频免费| 免费久久久久久久精品成人欧美视频| 成年动漫av网址| 亚洲av电影在线观看一区二区三区| 99re6热这里在线精品视频| 在线观看免费高清a一片| 免费女性裸体啪啪无遮挡网站| 777久久人妻少妇嫩草av网站| 久久精品亚洲av国产电影网| 久久精品国产a三级三级三级| 大码成人一级视频| 精品乱码久久久久久99久播| 大片免费播放器 马上看| 欧美xxⅹ黑人| 天堂俺去俺来也www色官网| av福利片在线| 波多野结衣av一区二区av| 多毛熟女@视频| 91精品三级在线观看| netflix在线观看网站| 亚洲av日韩在线播放| 亚洲五月色婷婷综合| av福利片在线| 巨乳人妻的诱惑在线观看| 91大片在线观看| 国产亚洲欧美在线一区二区| 在线 av 中文字幕| 亚洲色图 男人天堂 中文字幕| 精品亚洲乱码少妇综合久久| 多毛熟女@视频| 捣出白浆h1v1| 超碰成人久久| 国产黄频视频在线观看| 十八禁高潮呻吟视频| 法律面前人人平等表现在哪些方面 | 国产一区二区三区av在线| 自线自在国产av| 午夜激情久久久久久久| 91精品三级在线观看| 女人爽到高潮嗷嗷叫在线视频| av网站免费在线观看视频| 不卡一级毛片| 久久久久久亚洲精品国产蜜桃av| 他把我摸到了高潮在线观看 | 乱人伦中国视频| 亚洲国产欧美在线一区| 中文字幕人妻丝袜一区二区| 一区二区三区乱码不卡18| 亚洲激情五月婷婷啪啪| 久久人人97超碰香蕉20202| 人妻久久中文字幕网| 人人妻人人爽人人添夜夜欢视频| 美国免费a级毛片| 亚洲精华国产精华精| 老熟妇仑乱视频hdxx| 精品欧美一区二区三区在线| 国产91精品成人一区二区三区 | 亚洲成人手机| 中国美女看黄片| 亚洲一区中文字幕在线| cao死你这个sao货| 岛国在线观看网站| 俄罗斯特黄特色一大片| 黄色怎么调成土黄色| 精品国产一区二区三区四区第35| av片东京热男人的天堂| 好男人电影高清在线观看| 99热国产这里只有精品6| 亚洲 国产 在线| 三级毛片av免费| 十八禁高潮呻吟视频| 国产成人免费观看mmmm| 91成人精品电影| 爱豆传媒免费全集在线观看| 日本欧美视频一区| 国产亚洲欧美精品永久| 美女中出高潮动态图| www.熟女人妻精品国产| 国产一卡二卡三卡精品| 亚洲国产精品一区三区| 男女床上黄色一级片免费看| 亚洲精品久久午夜乱码| 国产欧美亚洲国产| 大片免费播放器 马上看| 美女脱内裤让男人舔精品视频| 夜夜骑夜夜射夜夜干| 久久精品国产亚洲av高清一级| 丝袜脚勾引网站| 视频在线观看一区二区三区| 天堂俺去俺来也www色官网| 美女主播在线视频| 乱人伦中国视频| 久久久久久亚洲精品国产蜜桃av| 国产深夜福利视频在线观看| 亚洲av日韩在线播放| 丝袜脚勾引网站| 性色av乱码一区二区三区2| a 毛片基地| 久久九九热精品免费| 午夜免费成人在线视频| 可以免费在线观看a视频的电影网站| 在线观看一区二区三区激情| 大码成人一级视频| 国产91精品成人一区二区三区 | 正在播放国产对白刺激| 国产亚洲av片在线观看秒播厂| 在线观看免费午夜福利视频| 日韩 亚洲 欧美在线| 最近最新免费中文字幕在线| 女性被躁到高潮视频| 国产一区二区三区综合在线观看| 正在播放国产对白刺激| 亚洲中文av在线| 成年人午夜在线观看视频| 成人手机av| 婷婷成人精品国产| 久久久久国内视频| 国产欧美日韩一区二区三区在线| 国产欧美日韩综合在线一区二区| 亚洲五月色婷婷综合| 国产黄色免费在线视频| 97人妻天天添夜夜摸| 又大又爽又粗| 欧美 日韩 精品 国产| 热99re8久久精品国产| 丝袜喷水一区| www.精华液| 亚洲男人天堂网一区| 少妇猛男粗大的猛烈进出视频| 美女脱内裤让男人舔精品视频| 人妻人人澡人人爽人人| 国产成人精品无人区| 黄网站色视频无遮挡免费观看| 最黄视频免费看| 久久久国产欧美日韩av| 欧美国产精品一级二级三级| 国产精品二区激情视频| 真人做人爱边吃奶动态| 少妇裸体淫交视频免费看高清 | 中文字幕最新亚洲高清| 精品国产乱码久久久久久男人| 国产一级毛片在线| 国产亚洲午夜精品一区二区久久| 国产一区二区三区在线臀色熟女 | 亚洲伊人久久精品综合| 亚洲av美国av| 一边摸一边抽搐一进一出视频| 久久久精品国产亚洲av高清涩受| 亚洲第一青青草原| 五月天丁香电影| 亚洲国产av新网站| 一本一本久久a久久精品综合妖精| 色婷婷久久久亚洲欧美| 三级毛片av免费| 人妻一区二区av| www日本在线高清视频| 久久精品久久久久久噜噜老黄| 亚洲精品一二三| 欧美精品啪啪一区二区三区 | 久久精品亚洲熟妇少妇任你| 亚洲av国产av综合av卡| 伊人久久大香线蕉亚洲五| 成人影院久久| 最近最新中文字幕大全免费视频|