• <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
    波多野结衣高清无吗| 国内精品一区二区在线观看| 色哟哟哟哟哟哟| 无遮挡黄片免费观看| 国产成人精品久久二区二区91| 91av网站免费观看| 可以免费在线观看a视频的电影网站| 精品国内亚洲2022精品成人| 精品国产亚洲在线| 久久精品aⅴ一区二区三区四区| 午夜亚洲福利在线播放| 国产成人欧美在线观看| 女人爽到高潮嗷嗷叫在线视频| 悠悠久久av| 在线永久观看黄色视频| 亚洲性夜色夜夜综合| 成人手机av| 欧美国产日韩亚洲一区| 成人国产综合亚洲| 日韩精品青青久久久久久| 男女床上黄色一级片免费看| 成年人黄色毛片网站| 日本在线视频免费播放| 久久香蕉精品热| 婷婷六月久久综合丁香| 中亚洲国语对白在线视频| 欧美日韩精品网址| 色综合亚洲欧美另类图片| 亚洲aⅴ乱码一区二区在线播放 | 国产免费男女视频| 国产又色又爽无遮挡免费看| 国产精品 国内视频| 婷婷丁香在线五月| 午夜激情av网站| www日本黄色视频网| 91成年电影在线观看| 操出白浆在线播放| 一级毛片高清免费大全| 国产男靠女视频免费网站| 久久久国产精品麻豆| 亚洲成人久久爱视频| 国产探花在线观看一区二区| 亚洲 国产 在线| 国产亚洲欧美在线一区二区| 久久久久性生活片| 国产精品国产高清国产av| 97人妻精品一区二区三区麻豆| 极品教师在线免费播放| 亚洲avbb在线观看| 日韩精品中文字幕看吧| 草草在线视频免费看| 欧美在线一区亚洲| 国产久久久一区二区三区| 激情在线观看视频在线高清| 国产精品久久久久久精品电影| 母亲3免费完整高清在线观看| 精品日产1卡2卡| 在线十欧美十亚洲十日本专区| 日韩成人在线观看一区二区三区| 亚洲av熟女| 精品午夜福利视频在线观看一区| 在线观看66精品国产| 国产免费av片在线观看野外av| 色精品久久人妻99蜜桃| 50天的宝宝边吃奶边哭怎么回事| 久久欧美精品欧美久久欧美| 亚洲欧美一区二区三区黑人| 久久久国产成人精品二区| 欧美乱码精品一区二区三区| 美女午夜性视频免费| 狠狠狠狠99中文字幕| 亚洲欧美精品综合久久99| 成人高潮视频无遮挡免费网站| 搡老岳熟女国产| 51午夜福利影视在线观看| www日本在线高清视频| 精品人妻1区二区| 99国产综合亚洲精品| 欧美黄色淫秽网站| 亚洲色图av天堂| 国内揄拍国产精品人妻在线| 亚洲五月婷婷丁香| 啦啦啦观看免费观看视频高清| 999精品在线视频| 欧美黑人精品巨大| 亚洲,欧美精品.| 国产av不卡久久| 亚洲无线在线观看| av视频在线观看入口| 丝袜人妻中文字幕| 又爽又黄无遮挡网站| 波多野结衣巨乳人妻| 亚洲专区字幕在线| 小说图片视频综合网站| 伊人久久大香线蕉亚洲五| 色综合亚洲欧美另类图片| 亚洲九九香蕉| 午夜免费激情av| www国产在线视频色| 久久 成人 亚洲| 又紧又爽又黄一区二区| 亚洲自拍偷在线| 国产一区二区在线av高清观看| 美女 人体艺术 gogo| netflix在线观看网站| 国产三级在线视频| 在线观看午夜福利视频| 看黄色毛片网站| 女人高潮潮喷娇喘18禁视频| 最近最新免费中文字幕在线| 国产成人av激情在线播放| 我的老师免费观看完整版| 久久香蕉国产精品| 中出人妻视频一区二区| 精品不卡国产一区二区三区| 久久香蕉国产精品| 精品日产1卡2卡| 久久国产精品影院| 国产精品日韩av在线免费观看| 午夜激情福利司机影院| 国产精品一区二区三区四区免费观看 | 亚洲成a人片在线一区二区| 成人永久免费在线观看视频| 欧美午夜高清在线| 国产成人一区二区三区免费视频网站| 亚洲自偷自拍图片 自拍| 男人舔女人下体高潮全视频| 久9热在线精品视频| 中文字幕高清在线视频| 国产亚洲精品第一综合不卡| 免费看十八禁软件| 国产成人精品久久二区二区免费| 国产亚洲精品第一综合不卡| 一本久久中文字幕| 搡老熟女国产l中国老女人| 免费搜索国产男女视频| 欧美最黄视频在线播放免费| 搡老岳熟女国产| 1024视频免费在线观看| 精品久久久久久久人妻蜜臀av| 午夜影院日韩av| 亚洲国产精品成人综合色| 午夜久久久久精精品| 日韩欧美 国产精品| 听说在线观看完整版免费高清| 精品欧美一区二区三区在线| 亚洲精品中文字幕在线视频| 精品免费久久久久久久清纯| 脱女人内裤的视频| 免费看a级黄色片| 亚洲aⅴ乱码一区二区在线播放 | 日本a在线网址| 日本三级黄在线观看| 动漫黄色视频在线观看| 少妇裸体淫交视频免费看高清 | 亚洲精品一区av在线观看| 免费在线观看影片大全网站| 99久久精品热视频| 大型av网站在线播放| 最新在线观看一区二区三区| 色综合婷婷激情| 12—13女人毛片做爰片一| 国产精品免费一区二区三区在线| 亚洲精品中文字幕一二三四区| xxxwww97欧美| 麻豆一二三区av精品| 日本黄大片高清| 黄片小视频在线播放| bbb黄色大片| 一级毛片高清免费大全| 麻豆久久精品国产亚洲av| 欧美日韩中文字幕国产精品一区二区三区| 国产成人一区二区三区免费视频网站| 又紧又爽又黄一区二区| videosex国产| 中文在线观看免费www的网站 | 成人三级做爰电影| 午夜视频精品福利| 99久久无色码亚洲精品果冻| 国产精品永久免费网站| 亚洲va日本ⅴa欧美va伊人久久| 欧美+亚洲+日韩+国产| 久久香蕉国产精品| 欧美中文日本在线观看视频| 亚洲精品中文字幕在线视频| 啪啪无遮挡十八禁网站| 国产精品日韩av在线免费观看| 久久久久久亚洲精品国产蜜桃av| 欧美三级亚洲精品| 国产精品av久久久久免费| 日韩欧美在线二视频| 亚洲成人久久性| 嫩草影院精品99| 日本在线视频免费播放| 97碰自拍视频| 成年人黄色毛片网站| 最好的美女福利视频网| 国产欧美日韩一区二区三| 麻豆久久精品国产亚洲av| 国产一区二区在线av高清观看| 国产精品美女特级片免费视频播放器 | 久久中文字幕一级| 这个男人来自地球电影免费观看| 一本一本综合久久| 脱女人内裤的视频| 欧美精品亚洲一区二区| 一区福利在线观看| 亚洲国产欧美一区二区综合| 中文字幕最新亚洲高清| 日本黄色视频三级网站网址| 免费电影在线观看免费观看| 欧美成人性av电影在线观看| 97超级碰碰碰精品色视频在线观看| 亚洲天堂国产精品一区在线| 久9热在线精品视频| 午夜日韩欧美国产| 女人高潮潮喷娇喘18禁视频| 久久亚洲精品不卡| 亚洲黑人精品在线| 人妻久久中文字幕网| 欧美不卡视频在线免费观看 | 91麻豆av在线| 国产伦在线观看视频一区| 精品久久久久久,| 日日干狠狠操夜夜爽| 国产精品美女特级片免费视频播放器 | 国产区一区二久久| 无限看片的www在线观看| 搡老岳熟女国产| 制服诱惑二区| 国语自产精品视频在线第100页| 亚洲国产欧美人成| 九色国产91popny在线| 亚洲国产精品999在线| 欧美最黄视频在线播放免费| 国产精品一及| 一卡2卡三卡四卡精品乱码亚洲| 日韩欧美 国产精品| 母亲3免费完整高清在线观看| 成人精品一区二区免费| 中文字幕人妻丝袜一区二区| www.熟女人妻精品国产| 日韩成人在线观看一区二区三区| bbb黄色大片| 小说图片视频综合网站| 色播亚洲综合网| 亚洲成av人片在线播放无| 99热只有精品国产| 狂野欧美激情性xxxx| 国产精品久久久久久久电影 | 欧美日韩中文字幕国产精品一区二区三区| 国产精品 欧美亚洲| 人人妻人人澡欧美一区二区| 亚洲人成电影免费在线| 国产一区二区三区视频了| 免费人成视频x8x8入口观看| 国产伦一二天堂av在线观看| 国产一区二区激情短视频| 亚洲欧洲精品一区二区精品久久久| 精品人妻1区二区| 一本一本综合久久| 日韩中文字幕欧美一区二区| 国模一区二区三区四区视频 | 国产一区二区在线观看日韩 | 亚洲性夜色夜夜综合| 黄色a级毛片大全视频| 日韩 欧美 亚洲 中文字幕| 日韩有码中文字幕| 亚洲七黄色美女视频| 免费在线观看日本一区| 深夜精品福利| 男女那种视频在线观看| www日本黄色视频网| 99热6这里只有精品| 欧美黑人精品巨大| 久久这里只有精品中国| 中出人妻视频一区二区| 午夜两性在线视频| 母亲3免费完整高清在线观看| 黑人操中国人逼视频| 国产91精品成人一区二区三区| aaaaa片日本免费| 在线播放国产精品三级| 可以在线观看的亚洲视频| 国产免费男女视频| 亚洲精品久久成人aⅴ小说| 18美女黄网站色大片免费观看| 国产探花在线观看一区二区| 97碰自拍视频| 午夜亚洲福利在线播放| 午夜成年电影在线免费观看| 国产视频一区二区在线看| 哪里可以看免费的av片| 又爽又黄无遮挡网站| 精品少妇一区二区三区视频日本电影| 欧美黑人欧美精品刺激| 欧美日韩乱码在线| 欧美国产日韩亚洲一区| 悠悠久久av| 免费观看人在逋| 亚洲男人的天堂狠狠| 99热只有精品国产| 中文字幕精品亚洲无线码一区| 两性午夜刺激爽爽歪歪视频在线观看 | 欧美黄色片欧美黄色片| 99久久精品国产亚洲精品| 一本大道久久a久久精品| 女生性感内裤真人,穿戴方法视频| 亚洲精品国产一区二区精华液| 黄色 视频免费看| www.www免费av| 老司机在亚洲福利影院| 国产精品久久久久久人妻精品电影| 特级一级黄色大片| 伦理电影免费视频| 日韩欧美一区二区三区在线观看| 欧洲精品卡2卡3卡4卡5卡区| 久久久国产欧美日韩av| 国产亚洲欧美98| 18禁观看日本| 亚洲欧美日韩高清专用| 欧美日韩黄片免| 俺也久久电影网| 黄色成人免费大全| 欧美成人一区二区免费高清观看 | 国内揄拍国产精品人妻在线| 国产激情久久老熟女| 色在线成人网| 久久草成人影院| 亚洲美女黄片视频| 很黄的视频免费| 国产午夜福利久久久久久| 91av网站免费观看| 欧美性长视频在线观看| 欧美zozozo另类| 老司机深夜福利视频在线观看| 极品教师在线免费播放| 亚洲专区字幕在线| а√天堂www在线а√下载| 中文资源天堂在线| 国产精品一区二区免费欧美| 黄片小视频在线播放| 一本大道久久a久久精品| 97超级碰碰碰精品色视频在线观看| 日韩免费av在线播放| 久久久国产成人免费| 动漫黄色视频在线观看| 中文字幕精品亚洲无线码一区| 人妻夜夜爽99麻豆av| 美女黄网站色视频| 午夜精品一区二区三区免费看| 久久精品国产亚洲av香蕉五月| 国产成人啪精品午夜网站| 男人舔女人下体高潮全视频| 欧美激情久久久久久爽电影| 亚洲欧美一区二区三区黑人| 老熟妇乱子伦视频在线观看| 夜夜爽天天搞| 变态另类成人亚洲欧美熟女| 99riav亚洲国产免费| 国产精品免费视频内射| 18禁裸乳无遮挡免费网站照片| 婷婷精品国产亚洲av| 亚洲国产欧美一区二区综合| 欧美成人午夜精品| 午夜福利视频1000在线观看| 悠悠久久av| 成人av在线播放网站| 人妻夜夜爽99麻豆av| 久9热在线精品视频| 老司机靠b影院| 亚洲一卡2卡3卡4卡5卡精品中文| 午夜福利高清视频| 久久久精品国产亚洲av高清涩受| 国产三级在线视频| 老司机福利观看| 国产精品香港三级国产av潘金莲| 欧美国产日韩亚洲一区| 免费在线观看完整版高清| 99热这里只有是精品50| 一二三四社区在线视频社区8| 人妻丰满熟妇av一区二区三区| 深夜精品福利| 黄色毛片三级朝国网站| 国产精品精品国产色婷婷| 女人爽到高潮嗷嗷叫在线视频| 免费看美女性在线毛片视频| 亚洲欧美日韩东京热| а√天堂www在线а√下载| 又黄又爽又免费观看的视频| 男女视频在线观看网站免费 | 99riav亚洲国产免费| 精品久久久久久久人妻蜜臀av| 黄色女人牲交| 亚洲免费av在线视频| 欧美日本亚洲视频在线播放| 欧美一级毛片孕妇| 久久婷婷人人爽人人干人人爱| 伦理电影免费视频| 亚洲精品美女久久久久99蜜臀| 精品国产乱子伦一区二区三区| 十八禁网站免费在线| 国产99久久九九免费精品| 国产成人一区二区三区免费视频网站| 日韩欧美国产一区二区入口| 国产69精品久久久久777片 | 免费搜索国产男女视频| 亚洲成人免费电影在线观看| av福利片在线| 久久人人精品亚洲av| 久久久国产成人精品二区| 美女午夜性视频免费| 亚洲国产精品久久男人天堂| 国产成年人精品一区二区| 国产探花在线观看一区二区| 国产精品自产拍在线观看55亚洲| 757午夜福利合集在线观看| 欧美zozozo另类| 一级黄色大片毛片| 制服诱惑二区| 99国产极品粉嫩在线观看| 又爽又黄无遮挡网站| av在线播放免费不卡| 亚洲第一电影网av| 哪里可以看免费的av片| 久久久久免费精品人妻一区二区| 色播亚洲综合网| 国产亚洲精品久久久久久毛片| 亚洲欧美日韩高清专用| 女人高潮潮喷娇喘18禁视频| 一卡2卡三卡四卡精品乱码亚洲| 久久伊人香网站| 国产私拍福利视频在线观看| 丁香六月欧美| 性色av乱码一区二区三区2| 国产精品一区二区三区四区免费观看 | 精品欧美一区二区三区在线| 女生性感内裤真人,穿戴方法视频| 色综合欧美亚洲国产小说| 日韩欧美在线二视频| 村上凉子中文字幕在线| 成年人黄色毛片网站| 在线观看66精品国产| 男女之事视频高清在线观看| 看黄色毛片网站| 九色成人免费人妻av| 一本精品99久久精品77| 国产高清激情床上av| 国产精品亚洲av一区麻豆| 97碰自拍视频| av福利片在线| 三级国产精品欧美在线观看 | 99国产综合亚洲精品| 99国产精品一区二区三区| 麻豆av在线久日| 国产一区二区三区视频了| 99热6这里只有精品| 啦啦啦免费观看视频1| 久久伊人香网站| 亚洲精品美女久久久久99蜜臀| 欧美高清成人免费视频www| 国产亚洲精品久久久久5区| 手机成人av网站| 午夜免费激情av| 成人国产一区最新在线观看| 久久久久国产一级毛片高清牌| 国产又色又爽无遮挡免费看| 淫妇啪啪啪对白视频| 久久久国产精品麻豆| 在线十欧美十亚洲十日本专区| 久久精品91无色码中文字幕| 99久久久亚洲精品蜜臀av| 国产97色在线日韩免费| 久久婷婷人人爽人人干人人爱| 啦啦啦观看免费观看视频高清| 国产精品精品国产色婷婷| 欧美日韩精品网址| tocl精华| 国产亚洲精品第一综合不卡| 亚洲欧洲精品一区二区精品久久久| 19禁男女啪啪无遮挡网站| 熟妇人妻久久中文字幕3abv| 男人的好看免费观看在线视频 | www.精华液| 深夜精品福利| 亚洲专区字幕在线| 日本一二三区视频观看| 男女午夜视频在线观看| 日韩欧美国产在线观看| 亚洲狠狠婷婷综合久久图片| 亚洲精品一卡2卡三卡4卡5卡| 成人av在线播放网站| 97人妻精品一区二区三区麻豆| 最新美女视频免费是黄的| 久久久久久大精品| 欧美激情久久久久久爽电影| 99在线视频只有这里精品首页| 亚洲真实伦在线观看| 丰满人妻熟妇乱又伦精品不卡| 麻豆国产97在线/欧美 | 最近视频中文字幕2019在线8| 欧美激情久久久久久爽电影| 亚洲一卡2卡3卡4卡5卡精品中文| 亚洲自拍偷在线| 制服人妻中文乱码| 国产一区二区三区在线臀色熟女| 国产精品98久久久久久宅男小说| 欧美日韩福利视频一区二区| 欧美又色又爽又黄视频| 久久久久性生活片| 成人三级黄色视频| xxx96com| 一个人免费在线观看电影 | 亚洲国产欧洲综合997久久,| 久久精品亚洲精品国产色婷小说| 亚洲精华国产精华精| 亚洲 国产 在线| 国产激情久久老熟女| 床上黄色一级片| 日韩国内少妇激情av| 亚洲精品久久国产高清桃花| 精品人妻1区二区| 精品久久久久久久久久久久久| 欧美一区二区国产精品久久精品 | 黄色视频,在线免费观看| 国产精品98久久久久久宅男小说| 99久久国产精品久久久| 久久这里只有精品中国| 看黄色毛片网站| 18美女黄网站色大片免费观看| 波多野结衣高清无吗| 黄色丝袜av网址大全| or卡值多少钱| 亚洲人成网站在线播放欧美日韩| 在线看三级毛片| 日韩欧美一区二区三区在线观看| netflix在线观看网站| www.www免费av| 男女午夜视频在线观看| 欧美久久黑人一区二区| 精品免费久久久久久久清纯| 亚洲成人久久性| 熟女少妇亚洲综合色aaa.| 亚洲熟妇熟女久久| 欧美精品啪啪一区二区三区| 久9热在线精品视频| 男人舔女人下体高潮全视频| 精品久久久久久久人妻蜜臀av| 日本 欧美在线| 又粗又爽又猛毛片免费看| 国产精品久久电影中文字幕| 欧美性猛交黑人性爽| av免费在线观看网站| 中文字幕精品亚洲无线码一区| 色综合亚洲欧美另类图片| 中文资源天堂在线| 欧美成狂野欧美在线观看| 亚洲欧美日韩无卡精品| 50天的宝宝边吃奶边哭怎么回事| 大型av网站在线播放| 18禁观看日本| 国产人伦9x9x在线观看| 国内毛片毛片毛片毛片毛片| 日本黄大片高清| 欧美极品一区二区三区四区| 亚洲国产欧美人成| 黑人巨大精品欧美一区二区mp4| 免费在线观看成人毛片| 1024香蕉在线观看| 欧美绝顶高潮抽搐喷水| 久9热在线精品视频| 一区福利在线观看| 久久精品国产综合久久久| 毛片女人毛片| 午夜a级毛片| 久久久久久久久中文| 男人的好看免费观看在线视频 | 91老司机精品| 狠狠狠狠99中文字幕| 18禁观看日本| 制服人妻中文乱码| 狠狠狠狠99中文字幕| 91字幕亚洲| 伦理电影免费视频| 黄频高清免费视频| 久久久久久九九精品二区国产 | АⅤ资源中文在线天堂| 手机成人av网站| aaaaa片日本免费| 亚洲国产欧洲综合997久久,| 国产乱人伦免费视频| 亚洲成a人片在线一区二区| 亚洲熟妇中文字幕五十中出| 国语自产精品视频在线第100页| 欧美色欧美亚洲另类二区| 久久久久久久久久黄片| 亚洲国产精品合色在线| 很黄的视频免费| 十八禁网站免费在线| 亚洲精品一卡2卡三卡4卡5卡| 最新在线观看一区二区三区| 九色成人免费人妻av| 在线观看免费视频日本深夜| 国产精品亚洲av一区麻豆| 女生性感内裤真人,穿戴方法视频| 精品少妇一区二区三区视频日本电影| 久久精品国产99精品国产亚洲性色| 亚洲熟妇熟女久久| 小说图片视频综合网站| 特大巨黑吊av在线直播| 中文字幕熟女人妻在线| 18禁黄网站禁片午夜丰满| 成人国产综合亚洲| 成年人黄色毛片网站| 国产成人系列免费观看|