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

    Preliminary study on detection and cleaning of parasitic bunches

    2021-11-13 01:30:38JunHeYanFengSuiYanHuaLuDiYinZheDuanSaiKeTianDeChongZhuYingZhaoJunHuiYueJianSheCao
    Nuclear Science and Techniques 2021年10期

    Jun He· Yan-Feng Sui · Yan-Hua Lu · Di Yin· Zhe Duan ·Sai-Ke Tian· De-Chong Zhu· Ying Zhao· Jun-Hui Yue ·Jian-She Cao

    Abstract To establish a nuclear resonant scattering beamline at the High Energy Photon Source (HEPS), it is essential to develop tools for detection and cleaning of parasitic bunches, for meeting the stringent demands on bunch purity. To this end, a novel time-correlated singlephoton counting system was implemented at the electron storage ring of the Beijing Electron-Positron Collider II(BEPCII). The purity deterioration process over a weeklong operation was recorded by the system. In this study,the mechanism of impurity growth was analyzed by numerical methods and validated on measurements. The agreement between the experimental results and the calculation was fairly good. Two main sources of parasitic bunches,pre-accelerators and the Touschek scattering were confirmed. A bunch-cleaning technique, based on a sinusoidal signal mixed with a pseudo-square wave, was also developed and implemented, and its capability to improve the bunch purity to the level of 10-7 was experimentally demonstrated. We present the experimental setup, principle,and measurement results of a system for detection and cleaning of parasitic bunches.

    Keywords Parasitic bunches · Impurity growth · High energy photon source (HEPS)

    1 Introduction

    In synchrotron radiation light sources, bunch purity,defined as‘‘the ratio of the number of electrons in parasitic bunches to that in the main bunches’’, must be kept sufficiently low in some time-resolved experiments [1]. Bunch purity is 10-4for most third-generation light sources without a bunch-cleaning system [2-6]; however, nuclear resonant scattering spectroscopy experiments, unique in high-energy synchrotron radiation light sources, demand much higher bunch purities, reaching those of a singleelectron-level, because the ratio of the non-resonant signal to the delayed resonant signal is above 106, even for highresolution monochromators [7-10]. Achieving singleelectron-level purity is challenging; for that, the dynamic range of the detection system should be extended to 1010.At the same time,a system for cleaning parasitic bunches is compulsory, because purity is likely to worsen gradually.Impurities proliferate owing to the injector chain and the Touschek scattering in the storage ring [9-11].

    The 6-GeV-energy High Energy Photon Source(HEPS)with a circumference of 1360.4 m and a designed natural emittance of 34 pm·rad has been under construction in Beijing [12]. As shown in Fig. 1, there are two filling patterns with bunch numbers of 680 and 63, respectively.The single-bunch charges are 1.33 nC and 14 nC,while the total current is 200 mA. The bunch spacings are 6 ns and 72 ns. The lifetime of the nuclear resonance transition(14.4 keV) of the iron-57 isotope is approximately 141 ns.The bunch spacing of 72 ns in the high-charge mode at the HEPS is slightly longer than a half of the lifetime. There are two beamline installations that demand high bunch purity in 14 beamlines in Phase I of the HEPS project.The Structural Dynamics Beamline (SDB) is intended toperform laser pump X-ray probe experiments with the required purity of 10-6. The Hard X-ray High Energy Resolution Spectroscopy Beamline explores nuclear resonant scattering spectroscopy experiments[13]and requires bunch purities of 10-9in the pattern of 63 bunches.Related beam diagnostic methods should be applied prior to the start-up of the new photon source HEPS. For this purpose,systems for bunch purity measurement and cleaning were installed in the storage ring of the Beijing Electron Positron Collider (BEPCII), which also serves as a synchrotron radiation light source parasitically during high-energy physics runs (the HEP mode), as well as in the dedicated synchrotron radiation mode (the SR mode). For bunch purity studies, all of the results were acquired in the SR mode. Table 1 lists the main parameters of BEPCII in theSR mode. The RF frequency was 499.8 MHz and the minimal bucket spacing was approximately 2 ns. The bunch spacing during a typical operation varied across 4,6,and 8 ns.Normally,only 1/4 of all the buckets were filled.

    Table 1 Main nominal parameters of the BEPCII

    Detection of parasitic bunches utilized the time-correlated single-photon counting (TCSPC) [14] technique while the cleaning system used the RF Knock Out(RFKO)[15, 16] technique. For the detection of parasitic bunches,the arrival time of single photons emitted by the electron bunches was measured relative to a clock pulse that was synchronized with the bunch revolution frequency via the storage-ring RF system. The TCSPC system generated a histogram of binned arrival time measurements that statistically represented the relative time distribution of electrons in the storage ring. For the cleaning of parasitic bunches, the conventional approach is to stimulate all of the bunches at the betatron frequency of parasitic bunches.The beam impedance of the vacuum chamber induces a shift in the betatron frequency, which depends on the single-bunch current, which has a significantly different betatron frequency (several kHz) from that of the main bunches. The amplitude of the transverse oscillations of parasitic bunches is much larger because of resonance.Consequently, the parasitic branches are filtered out by a scraper. The conventional approach was capable of achieving the required level of purity, but was difficult to tune and it was very sensitive to betatron tune shifts induced by the thermal drift of the storage ring.In addition,even for a properly tuned system, it induced significant main beam losses and required several minutes to complete a cleaning cycle. By contrast, the RFKO approach only required dozens of milliseconds. The RFKO method only stimulated the parasitic bunches and had little effect on the main bunches. In Sect. 2, we briefly describe the mechanism of the bunch purity deterioration, and the results of purity on BEPCII are analyzed using T. Obina’s method[11].The Touschek impurity growth rate is calculated,and the linac impurity growth rate is derived from the difference between the decay and top-up modes. The top-up injection control refilled the smallest bunch in sequence until an additional current of 1 mA was added when the total current dropped to 249.5 mA. The bunches were not refilled during the 4-h-long decay mode operation. Section 3 presents the hardware and the results for the system that was used to clean the parasitic bunches. The conclusions and future improvement approaches are discussed in Sect. 4.

    2 Detection of parasitic bunches

    2.1 Theory of the bunch purity deterioration

    As mentioned before, parasitic electrons in the storage ring could be inherited from pre-accelerators. The sources could be the electron gun, the dark current from the buncher in the linac, or the Touschek scattering in the booster. Parasitic electrons could also be produced owing to the Touschek scattering in the main bunches of the storage ring [9-11]. The contribution ratios of these are different for different machines. At the BEPCII, a full energy linac is injected directly into the storage ring, and we hope to determine the relative ratios of parasitic electron contributions caused by the linac beam and the Touschek scattering in the storage ring.

    2.1.1 LINAC

    The limit of refill f of the storage ring is 0.5 mA /250 mA = 0.2%. When the number of electrons in the main bunch N0drops to(1 - f)N0,the linac fills the main bunch to N0during the top-up mode. The refilling time dt and lifetime τ satisfy the following equation [10]:

    2.1.1.1 Touschek scattering The trajectories of electrons in the longitudinal phase space are described by the synchrotron oscillation equation [17]. When the momentum deviation is small, an electron undergoes a damping oscillation around the synchronous phase. An electron that gains a positive momentum larger than the height of the RF bucket exits the main bunch. The trajectories of the electrons in the longitudinal phase space are shown in Fig. 2.The electrons with positive momenta larger than the RF bucket height move toward the backward buckets + 1, + 2, + 3 (+ represents backward buckets).Without the synchrotron radiation,the electrons thrown out from the 0thbucket are never captured by the following buckets. However, radiation damping produces a momentum window Δpion the i-th bucket. Electrons with momenta corresponding to this window can be captured by the + 1, + 2, + 3 buckets.The increasing rate of parasitic bunches owing to the Touschek scattering and recaptures in subsequent buckets is given by [10, 11]

    Fig. 2 (Color online) Trajectories of electrons that are captured by the following buckets in the longitudinal phase space

    The electrons in the parasitic bunches can be lost owing to the beam-gas or the Touschek scattering. Therefore, the growth rate of the electrons in a parasitic bunch is given by

    2.2 Experimental setup

    The purity measurement system is shown schematically in Fig. 3a. The visible light that came from a bending magnet was sent to a photon detector(Microchannel Plate-Photomultiplier Tube MCP-PMT, China Electronics Technology Group Corporation or PMA Hybrid 400,Picoquant [19]) after passing through a full reflection mirror, a focus lens, a variable attenuator, a grating spectrometer, and two adjustable apertures. It shared a photon beam transport system with a beam-profile measurement system [20]. Two slits and neutral density (ND) filters located in front of the photon detector were used for reducing the photon flux down to the level of one photon detection per approximately dozens of revolutions of the beam. A grating spectrometer, for rejecting the stray light,was set on the optical line to the detector. As presented in Fig. 3b, the output signal of the detector was sent to a multichannel picosecond event timer (HydraHarp 400,Picoquant or Constant Fraction Discriminator CFD,ORTEC935, time to amplitude converter TAC, ORTEC 567, and multi-channel analyzer MCA, ORTECTRUMPPCI-2 k) modules [21], where the time difference between the trigger signal and the photon detector signal was measured using a time-to-digital converter (TDC). The detector recorded the distribution of photons that corresponded to the distribution of electrons over time in the revolution period.The reference signal was provided by the event generator / event receiver from the timing system,which provided a revolution period signal of 1.2433 MHz in the SR mode. A digital delay generator (DG645, SRS)was used for adjusting the delay between the detection of photons and the reference signal.

    Originally, the detector electronics were composed of a nuclear instrumentation module (NIM), and the reason for the recent upgrade to the HydraHarp module was a short dead time of 80 ns, allowing a maximal count rate of 12.5 × 106s-1. The count rates were 10 counts/s for background noise (without synchrotron light) and 106counts/s for purity measurements.The maximal number of time bins was 216= 65,536, which allowed the minimal time resolution of 16 ps for a revolution time of 800 ns.The count depth per time bin was 32 bits,corresponding to 232= 4,294,967,296 increments,which was very useful for high-purity measurements, compared with the Pico-Harp300 module with the storage capability of 16 bits[19].

    2.3 Purity measurements

    Figure 3c shows the purity measurement results for the top-up operation at the BEPCII. The filling pattern was designed for pump-probe experiments[22].The reading of the bunch current monitor (BCM) is shown in the inset of Fig. 3c; this provided the filling pattern information for operators. The total beam current of 250 mA was distributed equally across a single isolated bunch that was located at the center of a 200-ns-wide gap and another 99 bunches in the bunch train,with the bunch spacing of 6 ns.

    The single-photon responses of the two detectors are shown in Fig. 4a. The rise time (10% to 90%) was 380 psfor the hybrid detector,and the response time of the MCPPMT was at the same level but had a different tail. The single-bunch purity measurement results are shown in Fig. 5,both with a tail of several tens of nanoseconds.The bunch root mean square (RMS) length was approximately 50 ps, as measured by a streak camera. The full width at half maximum (FWHM) values of the bunch responses were 427 ps and 128 ps, respectively. The difference was mainly caused by the NIM’s response being much slower than that of HydraHarp. The width of a single bunch is determined by the instrument response function (IRF) of the system and the time structure of the bunch photons.The parasitic bunch signals are very easily smeared by spurious signals induced by the main bunch,especially for detectors with fast responses and high photon-electron conversion gains. This effect is alleviated for slow detectors, such as avalanche photo diodes(APDs).There are several different names and possible reasons for these spurious signals,such as ghost reflection[2]and after-pulsing[14].The details of the hybrid detector’s measurement results near the main signal are shown in Fig. 4b.There are four spurious peaks,located at 0.93 ns, 1.34 ns, 2.94 ns, and 6.81 ns after the main bunch, compared with the MCP-PMT detector that revealed peaks at 4 ns and 26 ns, respectively (data not shown). Several reasons explain why these spurious peaks do not correspond to real electrons. First, they appear outside stable RF phase regions. Second, their characteristics do not depend on injections,and the ratio of spurious to the main peaks sometimes shifts slightly depending on the photon flux but not by the time after an injection, as shown in Fig. 6a.Third,they are not affected by the RFKO system.Finally,different labs have reported similar results,although different detectors yielded spurious peaks at different positions [2, 11]. Yet, baseline counts increase,hindering the accurate evaluation of parasitic branches.Fortunately, none of the spurious peaks were at the locations of the parasitic bunches for the PMA hybrid detector,as shown in Fig. 4b.

    Fig. 4 (Color online) a Single-photon responses of the two detectors. b Spurious peaks of the PMA hybrid detector

    Fig. 5 (Color online) Single-bunch purity measurement results for (a) the MCP-PMT + NIM modules, and (b) Hybrid + HydraHarp

    Figure 6a shows the bunch purity measurement results at different times during a week-long top-up operation.The main bunch (No. 0, the first bunch in Fig. 3c) located at 635.2 ns is in the center of a 200-ns-wide gap. Besides the parasitic bunches located at 635.2 + 2 × i ns, there are three spurious peaks located at 636 ns, 638.5 ns, and 642 ns. These spurious peaks could not be cleaned by the bunch purification system. Therefore, they do not correspond to real electrons. One possible explanation for the appearance of these peaks is the inner reflection of the photon detector. We excluded the possibility of reflection photons created in the photon transport process. Further studies need to be done for confirming the origin of these spurious peaks.As shown in Fig. 6a,the counts associated with parasitic bunches increase gradually, but those associated with the spurious peaks do not. The ratios of the numbers of electrons in different parasitic branches to that in the main bunch are shown in Fig. 6b. The experimental results confirmed the Obina and Keil theory [10, 11].Impurity growth was never observed in forward (No. -1)bunches. Impurity growth rates were different for parasitic bunches, and were mainly determined by the different momentum gain probabilities through the Touschek scattering.The whole-bunch purity deterioration could be seen clearly throughout the week-long measurement. Equation (6)suggests that purity approaches equilibrium with a rise time of τg. The typical gas lifetime τgfor the BEPCII was 6.77 h, while the elastic and inelastic scattering gas lifetimes were 14.32 h and 12.84 h, respectively, for the pressure of 1.3 × 10-8Pa.

    Fig. 6 (Color online) a The bunch purity measurement results for a week-long operation,showing the details of a 14-ns-wide span around the main bunch, located at the center of a 200-ns-wide gap.b Evolution of the different parasitic bunches’ purities (ratios of the numbers of electrons in different parasitic branches to that in the main bunch) during the week-long operation

    To elucidate the contributions of the injection and the Touschek scattering to the ratios,the average purities of the 99 bunches in the bunch train in the top-up and decay modes are shown in Fig. 7a. (Figure 6b depicts the purity of a single isolated bunch.)The average purities in the 4-hlong decay mode were compared with that of the normal top-up operation under the same filling pattern.

    The main bunch spacing was 6 ns, and there were two possible parasitic bunches between the two main bunches,each of which could be the contribution of two forward main bunches. The bunch purity (+ 2 ns) corresponded to C1T+ C4T, while the bunch purity (+ 4 ns) corresponded to C2T+ C5T. According to Eq. (6), we obtained Eqs. (7)and (10):

    3 Cleaning of parasitic bunches

    3.1 The principle and experimental setup

    The bunch-cleaning technique developed at SPring-8[15] is shown schematically in Fig. 8. A sinusoidal signal comprising the Nthharmonics of the revolution frequency and the vertical betatron frequency as f = N × frev+ fβyis mixed with a pseudo-square wave that defines the bunch fill pattern. The main bunches, which are represented by blue dots, are located at the middle of the rise or fall edge(also the zero-crossing point) of the pseudo-square signal,and all the other bunches (red dots) at both positive and negative polarities of the square signal are cleaned. The mixed signal is amplified and sent to a kicker to stimulate a beam. If the amplitude of the excited oscillations is sufficiently large, parasitic bunches are lost in a vacuum chamber. The main apparatus of the bunch-cleaning system, including the timing system, the kickers, and the amplifiers, is shared with a transverse bunch-by-bunch feedback system.

    Fig. 7 (Color online) a Purity results for the top-up and decay modes, with the fill pattern being the same as in Fig. 1c. b Fit results for Pu2nsTop-up, using Eq. (6)

    Fig. 8 (Color online) Schematic of the bunch-cleaning system

    The bunch-cleaning system in the BEPCII consists of an arbitrary waveform generator (UHFAWG, 1.8 GSa/s, 14 bit,600 MHz,Zurich Instruments;the minimal time step of programming is 555 ps), a function/arbitrary waveform generator (33250A, Agilent), a digital delay generator(DG645), a solid-state amplifier (Amplifier Research 10 kHz-250 MHz 75 W),and a 60 cm strip-line transverse feedback kicker.

    3.2 Bunch-cleaning experiments

    Fig. 9 (Color online) a-b Typical signals generated by the bunchcleaning system. c BPM (upper trace) and the strip-line kicker downstream port (lower trace) signals of the cleaning system.d Detailed view of the rise edge of the excitation signal. The arrow shows the safe point where the main bunch is located. The results in Fig. 9a and Fig. 9b were obtained in the single-shot mode, while the results in Fig. 9c and Fig. 9d were obtained in the infinite persistence mode

    Figures 9a and 9b shows the typical signals generated by the bunch-cleaning system. Curve 1 in red shows the output of the UHFAWG arbitrary waveform generator,while Curve 2 in blue is a sinusoidal wave with a frequency of frev+ fβy, produced by the Agilent arbitrary waveform generator. Curve 3 in green shows the waveform after mixing.(The colors of the curves in Fig. 9 are the same as in Fig. 8). Curve 4 is not a real signal but the result of a mathematical operation(multiplication)on Curves 1 and 2.The maximal output Vpp of the UHFAWG was 750 mV,so there is a factor of 3 in Curve 4.To show this more clearly,the curves were translated along the vertical axis.Figure 9b shows the result of the scope in the single-shot mode, and Figs. 9c and 9d shows the infinite persistence modes. The time scales were 100 ns and 4 ns per division,respectively.Figure 9c shows the BPM (upper trace) and the strip-line kicker (lower trace) signals of the cleaning system. The BPM signal was used for monitoring the effect of cleaning,and the signal originating from the downstream port of the strip-line kicker was used for observing the temporal relations between the excitation signal and the bunches.First,we set the frequency of the sinusoidal wave far away from the resonance frequency,so the bunches located at the flat part of the kicker signal would not be cleaned, as shown in Fig. 9c. Several bunches were located at the flat of the square wave, and we then changed the delay to ensure that the main bunch was located at the lowest point,as indicated by the arrow in Fig. 9d. After finding the proper timing set, we changed the frequency of the sinusoidal wave to the right resonance frequency, which was cleaned in less than 100 ms.The lost electrons in the main bunch were ignored, and the perturbation was very small,demonstrating the possibility of using the system transparently during user operation.

    After confirming the functionality of the bunch-cleaning system,a special fill pattern was designed for testing.Five bunches with a spacing of 100 ns shared a total current of 15 mA equally, and the bunches were not refilled during the entire test. Figure 10a shows the purity results before and after the cleaning system came into effect. Purity before cleaning was measured 10 min after the injection,and the purity measurement process took several minutes.Figure 10b shows the details of the remaining bunch.Only the third main bunch survived. Table 3 shows the purity results after and before the cleaning system was used. All of the parasitic bunch purities after cleaning reached 10-7.The parasitic bunches after No. + 2 were submerged in noise. The purity of the theoretical value given by CiTin Table 2 is shown in the last column of Table 3.The purity results show that the cleaning system works well and the sensitivity of the purity measurement system is 10-7, limited by the photon detector’s characteristics. Detection of visible light suffered from a relatively high background caused by the detector itself or by stray light. Both the spurious peaks and dark counts increased the uncertainty of the purity. Using the method mentioned in Sect. 2.3, the uncertainty was in the 8 × 10-8-9 × 10-8range, much better than that for the normal operation.

    The absolute uncertainty (or sensitivity) of the system was difficult to evaluate, but the factors affecting it could be identified, as follows. The impulse response of the photon detector, with spurious peak counts much higher than normal noise counts, is expected to strongly deteriorate the uncertainty. The effective dynamic range of the system, which is determined by the minimum across all components, is smaller for the photonic detector than for the electronic one for the purity system, and count rate above 107is expected to trigger the PMA Hybrid 400 detector self-protection.The overall current for the normal operation was 17 (20 for the bunch amount) times higher than that in Fig. 10, and the uncertainty was 70 times worse. The time spent on purity measurements is typically 4-8 min, and purity increases to the 9.8 × 10-7-19.6 × 10-7range, according to the growth rate C1T-= 2.64 × 10-7in Table 2. The discrimination threshold setting of the electronic system also plays an important role. Quantification using the ratio of average background counts to the main bunch counts only provides a reference,but it is a good approximation.

    Fig. 10 (Color online) a Purity results, after and before the cleaning system operation. The purity results before cleaning were obtained in 10 min after the injection. b Details around the remaining bunch

    Table 3 Purity results for different bunches, before and after the cleaning procedure

    4 Conclusion and outlook

    A novel bunch-purity measurement system was implemented at the BEPCII, and we successfully tested the HydraHarp module combined with the PMA hybrid photon detector to evaluate the growth of impurities. HydraHarp 400, with its 32-bit histogram capability, 1 ps resolution,and the acquisition rate of 107counts/s appears to be an excellent instrument adapted for purity measurements.The growth of impurities was recorded by the TCSPC system over a week-long operation. Two sources of parasitic bunch populations have been identified: the linac and the Touschek scattering of the main bunch. The parasitic bunches exhibited different growth rates, which has been confirmed by theoretical calculations and measurements.The difference between the growth rates of the decay and top-up modes indicates that the contributions to the linac were C1inj≈ 10-5/h and C2inj≤10-6/h. The transfer fraction of the linac was k ≈10-6. The measurements showed that the bunch purity gradually worsened after the first injection, approaching equilibrium after a time that was several-fold the gas lifetime τg. A bunch-cleaning system based on the RFKO technique was successfully tested,and the achieved purity was 10-7.Overall,the novel procedure was much easier to use than the previous one,and allowed very fast cleaning of bunches.The purification procedure only took dozens of milliseconds compared with dozens of seconds required by the traditional method. The effect on the main bunches was much smaller.

    The sensitivity of the purity measurement system was better than 10-7in the best case, and was limited by the dynamic range and noise counts of the detector. An APD detector that is only sensitive to X-rays and a gated device that could provide a better contrast for main bunches relative to parasitic bunches, are under development. The APD detector will be helpful for improving the sensitivity by decreasing the counts associated with spurious peaks and stray photons. A gated device or a fast light shutter system,such as the one using Pockels cells,will extend the dynamic range of the system, which is considered to be a necessary condition for achieving purities on the order of 10-9. The purity measurement and purification system for the HEPS is currently under construction.

    AcknowledgementsThe authors would like to thank Joachim Keil and Jiu-Qing Wang for their helpful discussions.

    Author contributionsAll authors contributed to the study conception and design. Material preparation, data collection and analysis were performed by Jun He and Yan-Feng Sui. The first draft of the manuscript was written by Jun He and all authors commented on previous versions of the manuscript. All authors read and approved the final manuscript.

    亚洲一区中文字幕在线| 一个人免费在线观看的高清视频| 免费人成视频x8x8入口观看| 三级毛片av免费| 成人亚洲精品av一区二区 | 欧美一区二区精品小视频在线| 欧美av亚洲av综合av国产av| 91精品三级在线观看| 精品免费久久久久久久清纯| 女警被强在线播放| 成人影院久久| 老司机在亚洲福利影院| 韩国精品一区二区三区| 在线观看免费视频日本深夜| 亚洲成人免费av在线播放| 美女高潮喷水抽搐中文字幕| 久久精品91无色码中文字幕| av网站在线播放免费| 一级毛片高清免费大全| 不卡一级毛片| 一进一出好大好爽视频| 免费女性裸体啪啪无遮挡网站| 久久午夜亚洲精品久久| 制服人妻中文乱码| 老熟妇乱子伦视频在线观看| 日韩成人在线观看一区二区三区| 男女床上黄色一级片免费看| 在线看a的网站| 欧美日韩福利视频一区二区| 亚洲欧美精品综合久久99| 精品国产乱子伦一区二区三区| 老司机亚洲免费影院| 麻豆久久精品国产亚洲av | 国产精品亚洲av一区麻豆| 久久人人爽av亚洲精品天堂| 亚洲专区国产一区二区| 国产不卡一卡二| 成人亚洲精品一区在线观看| 啪啪无遮挡十八禁网站| 久久久国产成人精品二区 | 中文字幕人妻丝袜一区二区| 一级片免费观看大全| 久久久精品国产亚洲av高清涩受| 亚洲精品久久午夜乱码| 亚洲精品成人av观看孕妇| 天天躁狠狠躁夜夜躁狠狠躁| 淫秽高清视频在线观看| 热99re8久久精品国产| 久久久久精品国产欧美久久久| 日日爽夜夜爽网站| 国产精品久久久久成人av| 亚洲精品一卡2卡三卡4卡5卡| 国产av在哪里看| 天堂√8在线中文| 久久精品成人免费网站| 18禁裸乳无遮挡免费网站照片 | 麻豆一二三区av精品| 99国产精品一区二区三区| 午夜精品国产一区二区电影| 亚洲专区国产一区二区| 一边摸一边做爽爽视频免费| 天堂影院成人在线观看| 麻豆国产av国片精品| 国产伦人伦偷精品视频| 亚洲人成77777在线视频| 国产一区二区激情短视频| 欧美黑人欧美精品刺激| 精品乱码久久久久久99久播| 久久精品国产亚洲av香蕉五月| 一区福利在线观看| 99久久人妻综合| 9热在线视频观看99| 女性被躁到高潮视频| 亚洲成人免费电影在线观看| 他把我摸到了高潮在线观看| 久久久久久久精品吃奶| 超色免费av| 精品日产1卡2卡| 亚洲av片天天在线观看| 国产熟女xx| 99久久久亚洲精品蜜臀av| 99国产综合亚洲精品| 18禁裸乳无遮挡免费网站照片 | 欧美日韩av久久| 97碰自拍视频| av有码第一页| 最近最新中文字幕大全免费视频| 国产97色在线日韩免费| 亚洲国产欧美网| 亚洲欧美日韩另类电影网站| 变态另类成人亚洲欧美熟女 | 精品久久久精品久久久| 精品国产国语对白av| 国产精品一区二区在线不卡| 久久精品成人免费网站| tocl精华| 色综合欧美亚洲国产小说| 国产精品永久免费网站| xxxhd国产人妻xxx| 搡老熟女国产l中国老女人| 亚洲人成网站在线播放欧美日韩| 好男人电影高清在线观看| 久久精品aⅴ一区二区三区四区| 两个人免费观看高清视频| 亚洲成人国产一区在线观看| 免费高清视频大片| 成人免费观看视频高清| 99精国产麻豆久久婷婷| 高清黄色对白视频在线免费看| 欧美成狂野欧美在线观看| 丰满人妻熟妇乱又伦精品不卡| 欧美日韩瑟瑟在线播放| av欧美777| 国产亚洲精品久久久久5区| 国产精品秋霞免费鲁丝片| 男男h啪啪无遮挡| 国产男靠女视频免费网站| 国产精品av久久久久免费| 国产单亲对白刺激| 9热在线视频观看99| 亚洲国产精品999在线| 午夜福利在线免费观看网站| 一进一出抽搐动态| 国产精品国产高清国产av| av视频免费观看在线观看| 淫秽高清视频在线观看| 波多野结衣一区麻豆| 亚洲精品国产色婷婷电影| 亚洲欧美日韩高清在线视频| 每晚都被弄得嗷嗷叫到高潮| 88av欧美| 久久午夜亚洲精品久久| 丰满迷人的少妇在线观看| 在线观看66精品国产| 欧美在线黄色| a级毛片在线看网站| 校园春色视频在线观看| 精品福利观看| 亚洲美女黄片视频| 美女大奶头视频| 国产av精品麻豆| 很黄的视频免费| 看黄色毛片网站| 精品一区二区三区视频在线观看免费 | 国产人伦9x9x在线观看| 亚洲成国产人片在线观看| 首页视频小说图片口味搜索| 亚洲国产中文字幕在线视频| 男女高潮啪啪啪动态图| 高潮久久久久久久久久久不卡| 午夜成年电影在线免费观看| 亚洲欧美日韩高清在线视频| 久久精品亚洲熟妇少妇任你| 成人永久免费在线观看视频| 国产精品一区二区精品视频观看| 黄色 视频免费看| 亚洲成人免费av在线播放| a级片在线免费高清观看视频| 在线观看免费日韩欧美大片| 久久国产亚洲av麻豆专区| 欧美日韩精品网址| 亚洲精品在线美女| 色精品久久人妻99蜜桃| 99精品久久久久人妻精品| 两人在一起打扑克的视频| 法律面前人人平等表现在哪些方面| 黑人操中国人逼视频| 欧美日韩亚洲高清精品| 精品人妻在线不人妻| 最近最新中文字幕大全免费视频| 精品国产一区二区久久| www.熟女人妻精品国产| 韩国av一区二区三区四区| 久久国产精品影院| 精品国产超薄肉色丝袜足j| 午夜影院日韩av| 涩涩av久久男人的天堂| 88av欧美| 午夜福利在线观看吧| 欧美人与性动交α欧美精品济南到| 精品久久久久久电影网| 成人亚洲精品av一区二区 | 中文字幕人妻丝袜制服| av在线天堂中文字幕 | 免费高清在线观看日韩| 长腿黑丝高跟| 午夜亚洲福利在线播放| 在线观看午夜福利视频| 麻豆国产av国片精品| 久久久久久亚洲精品国产蜜桃av| 丰满人妻熟妇乱又伦精品不卡| 欧美日韩国产mv在线观看视频| 精品少妇一区二区三区视频日本电影| 变态另类成人亚洲欧美熟女 | 国产片内射在线| 午夜激情av网站| 午夜免费观看网址| 国产国语露脸激情在线看| 色综合婷婷激情| 日韩成人在线观看一区二区三区| 黄片大片在线免费观看| 免费观看精品视频网站| 巨乳人妻的诱惑在线观看| av网站免费在线观看视频| av天堂久久9| 超碰97精品在线观看| 久久久久久久午夜电影 | 亚洲国产欧美日韩在线播放| 欧美日韩中文字幕国产精品一区二区三区 | 九色亚洲精品在线播放| 一级黄色大片毛片| 久久精品国产99精品国产亚洲性色 | 国产精品免费视频内射| 在线观看免费午夜福利视频| 免费一级毛片在线播放高清视频 | 80岁老熟妇乱子伦牲交| 大型av网站在线播放| 精品欧美一区二区三区在线| 亚洲精品国产区一区二| 露出奶头的视频| 久9热在线精品视频| 亚洲精品av麻豆狂野| 男女午夜视频在线观看| 超色免费av| 国产成人免费无遮挡视频| 免费看a级黄色片| 老鸭窝网址在线观看| 色尼玛亚洲综合影院| 亚洲 欧美 日韩 在线 免费| 国产伦人伦偷精品视频| 亚洲在线自拍视频| 天天躁狠狠躁夜夜躁狠狠躁| 69精品国产乱码久久久| 国产精品一区二区精品视频观看| 天天添夜夜摸| 久久久久久久久免费视频了| 国产真人三级小视频在线观看| 久久久久九九精品影院| 韩国av一区二区三区四区| 一区二区日韩欧美中文字幕| aaaaa片日本免费| 亚洲欧美精品综合一区二区三区| 亚洲性夜色夜夜综合| 高清av免费在线| 少妇粗大呻吟视频| 色精品久久人妻99蜜桃| 久久精品亚洲精品国产色婷小说| 长腿黑丝高跟| 久久久久久亚洲精品国产蜜桃av| 超碰97精品在线观看| 日本免费a在线| 丰满迷人的少妇在线观看| 久久久久久大精品| 免费在线观看影片大全网站| 三上悠亚av全集在线观看| 变态另类成人亚洲欧美熟女 | 精品久久久久久电影网| 女人精品久久久久毛片| 久久久国产成人精品二区 | 少妇粗大呻吟视频| 国产成人系列免费观看| 国产免费av片在线观看野外av| 亚洲精品美女久久久久99蜜臀| 国产免费现黄频在线看| 交换朋友夫妻互换小说| 一区二区三区精品91| aaaaa片日本免费| 日本黄色视频三级网站网址| 国产成人精品无人区| 精品国内亚洲2022精品成人| 国产精品国产av在线观看| 国产精品秋霞免费鲁丝片| 一级黄色大片毛片| 午夜福利在线免费观看网站| 久久热在线av| 男女做爰动态图高潮gif福利片 | av网站在线播放免费| ponron亚洲| 久久精品成人免费网站| 热99国产精品久久久久久7| 久久香蕉国产精品| 黑人欧美特级aaaaaa片| 国产高清激情床上av| 午夜老司机福利片| 色哟哟哟哟哟哟| 啪啪无遮挡十八禁网站| 高清欧美精品videossex| 亚洲av成人av| av天堂久久9| 日韩成人在线观看一区二区三区| 国产成人一区二区三区免费视频网站| 悠悠久久av| 母亲3免费完整高清在线观看| 色播在线永久视频| 一二三四社区在线视频社区8| 一级片'在线观看视频| 久久欧美精品欧美久久欧美| 久9热在线精品视频| 亚洲精品一卡2卡三卡4卡5卡| 亚洲人成77777在线视频| xxx96com| 久久中文看片网| 老司机福利观看| 人人妻,人人澡人人爽秒播| 中文欧美无线码| 色老头精品视频在线观看| 嫩草影院精品99| 亚洲在线自拍视频| 在线免费观看的www视频| 99热国产这里只有精品6| 99久久99久久久精品蜜桃| 免费在线观看完整版高清| 在线观看免费午夜福利视频| 欧美精品啪啪一区二区三区| 99精国产麻豆久久婷婷| 中文字幕高清在线视频| 天堂√8在线中文| 亚洲欧洲精品一区二区精品久久久| 亚洲自拍偷在线| 日本一区二区免费在线视频| 纯流量卡能插随身wifi吗| 久久久久国内视频| 欧美精品亚洲一区二区| 久久久久久久久久久久大奶| 老司机在亚洲福利影院| 亚洲aⅴ乱码一区二区在线播放 | 亚洲人成电影观看| 黄色成人免费大全| 亚洲黑人精品在线| www日本在线高清视频| 91成年电影在线观看| 精品国内亚洲2022精品成人| 日本黄色日本黄色录像| av国产精品久久久久影院| 欧美激情 高清一区二区三区| 成熟少妇高潮喷水视频| 国产又色又爽无遮挡免费看| 99久久综合精品五月天人人| 无遮挡黄片免费观看| 热99re8久久精品国产| 亚洲一区二区三区色噜噜 | 久久精品亚洲精品国产色婷小说| 丝袜美腿诱惑在线| 在线观看舔阴道视频| 国产一区二区在线av高清观看| 国产aⅴ精品一区二区三区波| 日韩 欧美 亚洲 中文字幕| 女人精品久久久久毛片| 欧美日本亚洲视频在线播放| 首页视频小说图片口味搜索| 久久国产乱子伦精品免费另类| 日本wwww免费看| 悠悠久久av| 午夜精品在线福利| 免费在线观看影片大全网站| 精品熟女少妇八av免费久了| 纯流量卡能插随身wifi吗| 欧美色视频一区免费| 亚洲熟女毛片儿| 日韩人妻精品一区2区三区| 免费在线观看视频国产中文字幕亚洲| 国产有黄有色有爽视频| 美国免费a级毛片| 亚洲性夜色夜夜综合| 国产成人精品久久二区二区免费| 久久精品人人爽人人爽视色| 在线观看免费高清a一片| 亚洲熟妇熟女久久| 一级毛片精品| 久久久久久久午夜电影 | 欧美激情高清一区二区三区| 精品久久久精品久久久| 日本免费a在线| 老司机深夜福利视频在线观看| 黑丝袜美女国产一区| 麻豆一二三区av精品| 亚洲成av片中文字幕在线观看| 欧洲精品卡2卡3卡4卡5卡区| 久热这里只有精品99| 99久久综合精品五月天人人| 国产主播在线观看一区二区| svipshipincom国产片| avwww免费| 美女扒开内裤让男人捅视频| av视频免费观看在线观看| 法律面前人人平等表现在哪些方面| 午夜福利欧美成人| 午夜久久久在线观看| avwww免费| √禁漫天堂资源中文www| 男人舔女人下体高潮全视频| 午夜成年电影在线免费观看| 少妇被粗大的猛进出69影院| 成人影院久久| 91麻豆精品激情在线观看国产 | 淫秽高清视频在线观看| 亚洲欧美一区二区三区久久| 一夜夜www| 在线十欧美十亚洲十日本专区| 亚洲狠狠婷婷综合久久图片| 老司机深夜福利视频在线观看| 精品免费久久久久久久清纯| 亚洲av五月六月丁香网| 电影成人av| 亚洲欧美精品综合一区二区三区| 一区二区三区精品91| 精品一区二区三区视频在线观看免费 | 色综合欧美亚洲国产小说| 国产又爽黄色视频| svipshipincom国产片| 亚洲视频免费观看视频| 久久人人爽av亚洲精品天堂| 欧美大码av| 欧美黑人欧美精品刺激| 国产在线精品亚洲第一网站| 国产精品九九99| 在线观看免费高清a一片| 一个人免费在线观看的高清视频| 欧美最黄视频在线播放免费 | 丰满迷人的少妇在线观看| 岛国视频午夜一区免费看| 黄片播放在线免费| 国产黄a三级三级三级人| av片东京热男人的天堂| 在线国产一区二区在线| 成人特级黄色片久久久久久久| 亚洲熟妇熟女久久| 黄片大片在线免费观看| 成人亚洲精品av一区二区 | 天天添夜夜摸| 欧美成人免费av一区二区三区| 色老头精品视频在线观看| 国产免费av片在线观看野外av| 一边摸一边抽搐一进一小说| 婷婷丁香在线五月| 免费少妇av软件| 美女扒开内裤让男人捅视频| 亚洲第一青青草原| 可以在线观看毛片的网站| 欧美日韩中文字幕国产精品一区二区三区 | 多毛熟女@视频| av视频免费观看在线观看| 丰满迷人的少妇在线观看| 久久久水蜜桃国产精品网| 一级毛片精品| 日本免费a在线| 亚洲色图av天堂| 1024视频免费在线观看| 妹子高潮喷水视频| 无人区码免费观看不卡| 久久人人精品亚洲av| 可以免费在线观看a视频的电影网站| 日韩 欧美 亚洲 中文字幕| 伦理电影免费视频| 亚洲精品国产精品久久久不卡| 国产成人精品久久二区二区免费| 狠狠狠狠99中文字幕| 变态另类成人亚洲欧美熟女 | 91老司机精品| 免费高清在线观看日韩| www日本在线高清视频| 视频在线观看一区二区三区| 999久久久精品免费观看国产| 别揉我奶头~嗯~啊~动态视频| 亚洲专区中文字幕在线| 大码成人一级视频| 国产免费av片在线观看野外av| 50天的宝宝边吃奶边哭怎么回事| 久久久久久大精品| 久久久久国产精品人妻aⅴ院| a级毛片在线看网站| 变态另类成人亚洲欧美熟女 | 淫妇啪啪啪对白视频| 日韩三级视频一区二区三区| 免费一级毛片在线播放高清视频 | 中文字幕人妻熟女乱码| 欧美成人午夜精品| 久久午夜亚洲精品久久| 久久久久亚洲av毛片大全| 欧美日韩福利视频一区二区| 成人av一区二区三区在线看| 日韩一卡2卡3卡4卡2021年| 搡老岳熟女国产| 亚洲精品国产色婷婷电影| 99在线人妻在线中文字幕| 咕卡用的链子| 久久精品成人免费网站| 亚洲全国av大片| 可以免费在线观看a视频的电影网站| 欧美午夜高清在线| 在线观看免费高清a一片| 欧美成人性av电影在线观看| 精品国产超薄肉色丝袜足j| 热99re8久久精品国产| av在线播放免费不卡| 丰满迷人的少妇在线观看| 亚洲人成77777在线视频| 成年版毛片免费区| 中文字幕人妻熟女乱码| 欧美在线黄色| 亚洲在线自拍视频| 国产伦人伦偷精品视频| 亚洲成人免费电影在线观看| 久久这里只有精品19| 日韩 欧美 亚洲 中文字幕| 淫妇啪啪啪对白视频| 国内毛片毛片毛片毛片毛片| 日韩中文字幕欧美一区二区| 国产亚洲欧美98| 精品国产美女av久久久久小说| 美女福利国产在线| 最近最新免费中文字幕在线| 老熟妇仑乱视频hdxx| 人妻久久中文字幕网| 亚洲成av片中文字幕在线观看| 国产亚洲精品一区二区www| 18美女黄网站色大片免费观看| 欧美日韩亚洲高清精品| 熟女少妇亚洲综合色aaa.| 99国产精品免费福利视频| 搡老熟女国产l中国老女人| 久久久久国内视频| 国产男靠女视频免费网站| 久久伊人香网站| 琪琪午夜伦伦电影理论片6080| 欧美日韩瑟瑟在线播放| 黄片大片在线免费观看| 波多野结衣av一区二区av| 亚洲情色 制服丝袜| 大陆偷拍与自拍| 一级a爱视频在线免费观看| 一级作爱视频免费观看| 美女午夜性视频免费| xxxhd国产人妻xxx| 国内毛片毛片毛片毛片毛片| 亚洲精品在线美女| 性欧美人与动物交配| 深夜精品福利| 欧美人与性动交α欧美软件| 男人操女人黄网站| 国产不卡一卡二| 99精品在免费线老司机午夜| 啪啪无遮挡十八禁网站| www.自偷自拍.com| 免费av中文字幕在线| 露出奶头的视频| 精品少妇一区二区三区视频日本电影| 精品无人区乱码1区二区| 在线永久观看黄色视频| 色播在线永久视频| 久久久久九九精品影院| 精品久久久久久,| 人人澡人人妻人| 免费观看人在逋| 老汉色∧v一级毛片| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲精品美女久久久久99蜜臀| 女性被躁到高潮视频| 99在线视频只有这里精品首页| 成人三级黄色视频| 99精品在免费线老司机午夜| 老鸭窝网址在线观看| 男人舔女人的私密视频| 日本三级黄在线观看| 色综合站精品国产| 一级黄色大片毛片| 成人国语在线视频| 日韩三级视频一区二区三区| av免费在线观看网站| 精品一区二区三区四区五区乱码| 亚洲视频免费观看视频| 国产免费男女视频| 后天国语完整版免费观看| 久久狼人影院| 看黄色毛片网站| 久久精品成人免费网站| 在线视频色国产色| 一级毛片精品| 久久久久国内视频| 亚洲国产精品999在线| 可以在线观看毛片的网站| 麻豆av在线久日| av中文乱码字幕在线| 夜夜躁狠狠躁天天躁| 69av精品久久久久久| 波多野结衣高清无吗| 久久热在线av| 69av精品久久久久久| 国产精品影院久久| 9色porny在线观看| cao死你这个sao货| 国产av在哪里看| 精品国产乱子伦一区二区三区| 母亲3免费完整高清在线观看| 国产av在哪里看| 亚洲精品久久午夜乱码| 色精品久久人妻99蜜桃| 波多野结衣av一区二区av| 窝窝影院91人妻| 亚洲精华国产精华精| 精品人妻1区二区| 精品国产乱子伦一区二区三区| 国产精品 国内视频| 亚洲美女黄片视频| 窝窝影院91人妻| 99riav亚洲国产免费| 美女午夜性视频免费| 99久久99久久久精品蜜桃| 国产av又大| 在线永久观看黄色视频| 99国产精品一区二区三区| 国产欧美日韩综合在线一区二区| 亚洲成人国产一区在线观看| 日韩精品中文字幕看吧| av片东京热男人的天堂|