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

    Radio frequency conditioning of an S-band accelerating structure prototype for compact proton therapy facility

    2021-07-02 09:31:56YuZhangWenChengFangXiaoXiaHuangJianHaoTanShaoPengZhongChengChengXiaoJunQiangZhangChengWangYuQingYangZhenTangZhao
    Nuclear Science and Techniques 2021年6期

    Yu Zhang· Wen-Cheng Fang · Xiao-Xia Huang· Jian-Hao Tan ·Shao-Peng Zhong · Cheng-Cheng Xiao · Jun-Qiang Zhang· Cheng Wang·Yu-Qing Yang · Zhen-Tang Zhao,

    Abstract The development of a high-gradient accelerating structure is underway to construct a compact proton linear accelerator for cancer treatment. Extensive experiments and numerous studies are being conducted to develop compact linear accelerators for proton therapy. Optimization of the electromagnetic and mechanical design has been performed to simplify the manufacturing process and reduce costs. A novel high-gradient structure with a low relativistic proton velocity (β), v/c = 0.38, was designed,fabricated, and tested at high power. The first full-scale prototype was also successfully tested with high radio frequency(RF)power,a repetition rate of 50 Hz,and pulse length of 3 μs to reach a high-gradient of 46 MV/m using a 50 MW S-band klystron power supply obtained from the Shanghai Soft X-ray Free Electron Laser Facility. This is the first high-power test in China,which is in line with the expected experimental goal. This study presents preliminary high-power testing of S-band standing wave accelerating structures with 11 cells.This work aims to verify the feasibility of using a high-gradient RF accelerating structure in compact proton therapy facilities. The cold test of the prototype cavity was completed in advance. Details of the high-power RF test setup, the process of RF conditioning, and the high-power results are described.

    Keywords S-band·Standing wave(SW)·High-gradient·Accelerating structure · RF conditioning

    1 Introduction

    Proton therapy, also known as proton beam radiotherapy, is an essential component of particle radiotherapy.Owing to the excellent protection, it provides to normal tissues and its biological effects near the Bragg peak,proton therapy has gained the approval of an increasing number of doctors and is favored by patients[1]. In recent years, many proton therapy centers have been established globally [2, 3]. As of February 2021, there are approximately 111 facilities dedicated to proton or heavy-ion therapy worldwide. Proton therapy is the primary method of heavy-ion therapy, which accounts for approximately 89% of heavy-ion therapy. There are approximately 37 particle therapy facilities under construction, and approximately 26 are planned to be built-in future [1]. According to the current requirements for single-room facilities,ultrahigh dose rate (FLASH) radiotherapy [4], and proton computed tomography (CT) [5], research on compact and efficient accelerators for proton therapy equipment has gained considerable attention.In a single-room facility,the linear accelerator (linac) is directly installed on a rotating gantry, which makes the proton therapy facility more compact and cost efficient. In clinical practice, the accuracy of proton therapy is essential to the quality of treatment [6, 7]. The accuracy can be improved using proton CT, which solves the error of conversion from traditional CT to proton therapy. The minimum beam energy required for the beam to pass through the human body is 350 MeV.A section of the booster can be connected to the existing proton therapy facility,and the high-gradient linac can be used to achieve an energy upgrade.FLASH therapy requires ultra-high dose irradiation in a short time.A single short pulse (<100 ms) is delivered at an ultra-high dose rate that exceeds 100 Gy/s[8],which significantly exceeds the conventional dose rate of 0.05 Gy/s; only a linac can achieve such rapid energy modulation. The most representative devices at present are IMPULSE[9],TULIP[10],and ProBe[11].Further research will make proton therapy facilities more popular and preferential, which will significantly promote proton therapy development.

    The current development of high-gradient radio frequency(RF)technology can contribute to achieving a more compact and efficient proton therapy device. Compared with synchrotron radiation accelerators and cyclotrons,linacs can easily meet the requirements of compact proton therapy facilities. In the early 1990s, Hamm et al. first proposed the use of an S-band linac for proton therapy[12, 13]. The TERA Foundation proposed and designed a proton therapy facility based on a high-frequency linac[14, 15] and developed a scheme based on a high-gradient linac, which combines a cyclotron used as an injector with a high-gradient linac used as a booster [16, 17]. The main component of the entire accelerator facility is the linear accelerating structure with a high accelerating gradient,which can be used to make the linac more compact and cheaper [18].

    Research and development of high-gradient RF technology is essential to develop proton therapy technology.The Shanghai Advanced Proton Therapy Device(APTRON) is the first hospital-based proton therapy facility in China[19,20];however,its design is based on a synchrotron radiation accelerator.We conducted a series of investigations and designed an all-linac proton accelerator based on a high-gradient linear accelerating structure [21].The energy of the S-band proton therapy linac that we designed is adjustable from 70 to 250 MeV; the average current intensity is 20 nA, and the pulse width is 10 μs with a repetition frequency of 1000 Hz. The linac consists of an ion source, the radio frequency quadrupole (RFQ),the side-coupled drift Tube Linac (SCDTL), the main accelerating structure (2856 MHz), permanent magnet quadrupoles (PMQs), and a power source. The overall length of the accelerator is less than 7 m.The designed alllinac proton therapy facility scheme is shown in Fig. 1.For the design of the main accelerating structure, we chose a relatively high operating frequency (up to the S-band) to reduce the cavity size. However, the higher the operating frequency, the smaller is the radial size of the accelerating cavity, and the greater is the manufacturing difficulty. For this,we chose a standing wave(SW)accelerating structure with an operating frequency in the S-band. In addition to considering the breakdown field strength when selecting the accelerating electric field gradient, the cavity consumption must be considered for the design of the accelerating structure, which is applied to a proton linac with a relatively large duty cycle.This accelerator is characterized by its compact structure and low cost. The proton energy must also be adjusted by turning off a part of the highfrequency power source or changing the input power. The linac uses a series of PMQs to focus the beam size and reduce the emittance growth; the PMQs are also used to reduce the overall size and the subsequent operation energy consumption.

    The development of the main accelerating structure of the linac is complete, and the optimization of the accelerating structure, processing and manufacturing, low-power testing, and tuning is complete as well and is described in detail in Ref. [21]. Compared with the backward traveling wave (BTW) structure, we adopt a SW accelerating structure; the operating model of the BTW structure is a 5π/6 mode, while that of the accelerating structure introduced in this study is a π mode, which has a high shunt impedance. For low-beta accelerating structures, it is difficult to control the traveling wave phase and synchronize the proton beam with the electromagnetic field. This study describes a high-power test of a high-gradient accelerating structure. A high-power test was conducted on the Shanghai Soft X-ray Free Electron Laser Facility(SXFEL)facility platform. The high-gradient experiment on the accelerating structure is essential in the development of compact proton therapy equipment. The purpose of highpower testing is to verify the performance of all aspects of the accelerating structure.

    After a series of processing, the high-gradient accelerating structure retains a large amount of residual adsorbed gas inside the accelerating structure. There are also some raised particles on the surface of the cavity, which will become a field emission source and affect the performance of the accelerating structure. Therefore, before the highgradient accelerating structure is applied to the actual project,it needs to go through high-power testing to reduce the particle emission sources in the RF cavity.This enables the high-gradient accelerating structure to operate at a higher electric field. The purpose of RF conditioning is to enable the particle accelerating cavity to operate stably at a high input power, and it does not produce RF discharge in the required vacuum range.RF conditioning has four main functions [22, 23]:

    Fig. 1 (Color online) Layout of all-linac scheme

    (1) Remove the remaining adsorbed gas in the accelerating cavity.

    (2) Remove the dielectric impurities on the internal surface of the accelerating cavity.

    (3) Burn whiskers or micro-protrusions on the inner surface of the accelerating cavity and smoothen the surface roughness, pits, and mountains.

    (4) Detect the location and frequency of ‘‘breakdowns’’as well as analyse and find possible processing defects.

    In this study,we will primarily introduce the high-power test of the S-band SW accelerating structure conducted on the SXFEL facility platform. We present the design and installation of the high-power RF conditioning platform and the results and analysis of RF conditioning.

    2 High-power test platform

    The designed accelerating structure is a single periodic S-band SW structure.The accelerating structure consists of π-mode standing wave (SW) cavities and a circular waveguide fed with the TM01mode because the disk-loaded waveguide cavity of the π-mode operating frequency has the highest shunt impedance among all operating modes. However, it was previously observed in a cold test that the single periodic accelerating structure with a πmode operating frequency is unstable [21]. This instability indicates that the mode interval between adjacent operating frequencies is too small, and the mode jumps are likely to occur during the actual operation of the accelerating structure. Furthermore, the field distribution will experience interference from adjacent modes.

    This can be theoretically analyzed. In the SW accelerating structure,when the structure resonates in the π-mode,the forward and backward spatial harmonics merge. The amplitude of the composite field then reaches twice the amplitude of a single harmonic,which greatly increases the shunt impedance of the structure. According to resonant circuit theory, the dispersion relation of a single periodic SW accelerating structure can be obtained as follows:

    This is the narrowest mode interval. From the above analysis, it can be concluded that when the selected operating mode of the accelerating structure is the π-mode, it has a higher shunt impedance, but is less stable [24].

    Although the accelerating structure is unstable, we completed the low-power test and subsequent high-power test after careful calculation, simulation, and manufacturing.The experimental results are consistent with the design goal. Although the reflected power is large, the circulator can ensure that the reflected power will not return to the klystron; thereby preventing damage to the klystron. This verifies the feasibility of this accelerating structure; however, owing to its instability, it cannot be used in practical projects. We subsequently developed a biperiodic SW accelerating structure with an operating mode of π/2,which included the design of the accelerating structure and the optimization of physical quantities. The specific content is not included here. It is used to study the accelerating structure of the linac of a compact proton therapy facility.A low-β single periodic accelerating structure with an expected accelerating gradient of 50 MV/m was designed.The prototype,fabrication,low-power test,and tuning was completed for this design in 2020. The low-power test results after tuning at room temperature (25.8 °C) and standard atmospheric pressure are as follows:the operating frequency of the accelerating structure is 2852.88 MHz,the S parameter at this frequency is - 6.57 dB, and the standing wave rate (SWR) is 2.78. The accelerating structure is in a state of under-coupling, as shown in the Smith chart of the previous microwave measurement in Fig. 2.Although the SWR is relatively large(higher power reflection), there is a circulator between the accelerating structure and the klystron during the high-power test, and the reflected power will not cause damage to the klystron.The final measurement results show that the S-band SW accelerating structure is suitable for high-power testing.

    Fig. 2 (Color online) Smith chart

    The accelerating structure is part of our research on a high-performance proton therapy facility,which is a singleroom proton therapy facility with the unique feature of a linac mounted on a rotating gantry [14]. From this perspective,the high-gradient structure significantly improves the compactness and lightness of the rotating mechanism,while also benefitting the cost and reliability [25]. The high-power test was completed in January 2021. The following describes the layout and feasibility of the highpower testing.

    The first prototype structure was subject to a high-power test using an S-band 50 MW klystron power supply from the SXFEL. According to our calculations, if the accelerating gradient of the accelerating structure is 50 MV/m,the power input to the accelerating structure is approximately 9.15 MW. The formula is as follows:

    where ZTT is the effective shunt impedance,ZTT =59.983 MΩ/m,Plis the power loss per unit length,T is the transit time factor of the accelerating structure,T = 0.846, and E0T is the accelerating gradient design value of the accelerating structure, E0T =50 MV/m.

    On the SXFEL platform, the machine time is limited because there are other accelerating tubes that require highpower testing. Considering the power transmission attenuation, a quarter of the power of the 50 MW klystron(approximately 12.5 MW) can be distributed to the accelerating structure through a 6 dB power divider, which can meet the requirement of the accelerating structure for highpower testing. The remaining power of the klystron(37.5 MW) is allocated to another accelerating structure.Moreover, high-power testing of the proton accelerator tube is installed on the SXFEL platform using SXFEL’s existing technical conditions and equipment and the waveguide system.The original facility uses a 6 dB power divider; and this high-power test did not damage the original facility. After the high-power testing is complete,the normal operation of the SXFEL must be restored.More details on the SXFEL facility can be found in Ref. [26].The layout of the test platform of the SXFEL facility is shown in Fig. 3;1 indicates a klystron,2 is a laser device,3 is an electron gun, and 4 is an accelerating structure. The accelerating structure is installed in the position of the electron gun. During the installation procedure, nitrogen gas flow, from the inside to the outside of the system,prevented any contamination from dust [14]. Next, weintroduce the process of installing this accelerating structure in the SXFEL facility in detail.

    Fig.3 (Color online)Layout of the test platform of the SXFEL facility

    The installation of the S-band accelerating structure began on December 3,2020.To ensure the vacuum degree of the entire system, each component was individually installed and vacuumed, and leak checks were performed during the installation process.Gaskets are required to seal the splicing between the waveguide and accelerating structure;however,the entire device still leaked during the initial installation process. This occurred because the gaskets were not annealed, and the annealed gaskets were fitted at a later time. The installation process proceeded smoothly and required approximately 2 weeks.

    Before being installed on the high-power test platform,the accelerating structure’s RF parameters must be measured for the klystron’s parameter setting. The RF parameters were measured in the SXFEL tunnel. The device diagram for RF parameter measurement is shown in Fig. 4.

    The molecular pump is used for vacuuming, and the vacuum degree is 5.49×10-7mbar.A network analyzer is used for the RF measurement. The designed accelerating structure adopts a cooling method of four-in and four-out water circulation with built-in water pipes, which is a traditional internal water-cooling method. We installed an adjustable water tank in the SXFEL tunnel to control the temperature of the accelerating structure. The operating frequency of the accelerating structure refers to the frequency in a vacuum environment and at a constant temperature. Changes in temperature, relative humidity, and air pressure are considered for the microwave measurement, and the frequency needs to be corrected. Frequency correction primarily includes two quantities: pressure and temperature [24, 27]. Equation (4) describes the pressure correction:where Δfprepresents the change in pressure with respect to frequency,Pdis the pressure of dry gas,Pwis the pressure of water vapor, and T is the absolute temperature.Assuming f0=2856 MHz, Pd=760 Torr, T =298.95 K,a humidity of 77%, and Pw= 17.54 Torr, then Δfp=-0.958 MHz.

    Fig. 4 (Color online) RF parameter measurement device diagram of the prototype cavity

    Changes in the cooling water temperature and cavity working temperature may change the temperature of the accelerating structure. Equation (5) describes the temperature correction:

    The influence of temperature on the frequency is related to the sensitivity of the structure size. The accelerating structure is made of copper, and g is the corresponding expansion coefficient,g =1.77×10-5/°C.In the formula,ΔfTrepresents the change in temperature with respect to frequency.The measuring temperature is 298.95 K,and the working temperature is 293.15 K; thus,ΔfT=fT-f0=- 0.293 MHz [24, 27]. The frequency is increased by 1.251 MHz, and the corrected frequency is(2852.88 - (- 0.958) - (- 0.293))

    MHz = 2854.131 MHz.

    Considering the change in the working environment of the accelerating structure, the working temperature can be changed by adjusting the temperature of the water tank.Changes in the working temperature of the accelerating structure will result in changes in the accelerating structure operating frequency.We measured the operating frequency of the accelerating structure at different working temperatures. We constructed the temperature-frequency curve from the experiment and compared it with the theoretical value, as shown in Fire 5. The theoretical value is consistent with the measured value. It can be observed from Fig. 5 that the operating frequency of the accelerating structure can be increased by reducing the temperature of the cooling water.However, to achieve the ideal operating frequency,the temperature of the cooling water needs to be reduced below 0 °C, which will freeze the cooling water.

    Fig. 5 Temperature-frequency curve of the experiment

    After the installation in the SXFEL tunnel, we adjusted the temperature of the cooling water to a suitable temperature of 20 °C. The measured operating frequency was 2854.23 MHz; the difference of 100 kHz from the theoretical value may be owing to other errors, which is consistent with the law. The measured input reflection coefficient at this operating frequency is - 6.2 dB; therefore,the SW ratio is too large,which is an irreparable error caused during processing. This causes the reflected power to increase, and the considerable reflected power may damage the klystron. To solve this problem, two solutions are proposed: the operating frequency of the accelerating structure may be increased by reducing the temperature of the cooling water,alternatively,the klystron frequency can be changed so that it is consistent with the operating frequency of the accelerating structure. Therefore, the reflection is reduced to a level that does not damage the klystron. Considering that the prototype is only used to research RF technology and high-power performance, it can be used for beam testing under the correct conditions.Finally,we adopt a compromised solution by changing the frequency of the klystron to 2854.23 MHz. This is consistent with the operating frequency of the accelerating structure, which can reduce the reflection in the system.The large reflection coefficient may cause excessive power reflection. In the high-power test installation, a 4-port ferrite phase shifter circulator is installed in front of the accelerating structure to recover the reflected power,which can prevent damage to the klystron. The circulator is produced by AFT Microwave (Backnang, Germany). For further details, please refer to the installation diagram of the high-power test platform in Fig. 6b.

    Fig. 6 (Color online) Site installation diagram and schematic diagram of the S-band accelerating structure on the high-power test platform. a On-site installation diagram of the S-band accelerating structure on the high-power test platform. b Schematic of the S-band accelerating structure high-power test layout

    The on-site installation diagram of the S-band accelerating structure after the preliminary installation on the high-power test platform is shown in Fig. 6a. The scheme of the S-band accelerating structure high-power test layout is shown in Fig. 6b.In this figure,the ion pump is used to ensure that the vacuum of the entire system meets the experimental standards, and the Faraday cup is used to monitor the dark current. A power divider is installed behind the klystron to distribute a quarter of the power source to the accelerating structure, and the remaining power is distributed to other test equipment. The 4-port circulator behind the ceramic window is used to recover the reflected power and protect the klystron from breakdown due to excessive reflected power. We use two directional couplers, DC1 and DC2, to measure the input waves and the reflected and transmitted waves. The directional couplers were calibrated in the laboratory using waveguidecoaxial transitions and a network analyzer. The power transmission direction is from port 1 to port 2 of the directional coupler. Port 1 is upstream and is close to the direction of the klystron;port 2 is downstream and is close to the direction of the accelerator tube.Ports 4 and 1 are on the same side, as are ports 3 and 2. The transmission waveguide coefficient S11 of the directional coupler is -25.91 dB at the operating frequency, S31 and S41 are -50.40 dB and - 76.76 dB, respectively, and the isolation is approximately - 26.27 dB, which meets the measurement isolation requirements of the directional coupler.Moreover,a ceramic window is inserted to avoid exposing the klystron window to the cavity vacuum during installation and then to under conditioning.The purpose of passing sulfur hexafluoride (SF6) gas is to ensure the working requirements of the circulator because it cannot operate in a vacuum.On the last day of 2020,we conducted a series of inspections and adjustments for the entire experimental device. The high-power test officially started on January 5, 2021 and ended on January 15, 2021. The entire conditioning procedure is semi-automatic. Next, we introduce the entire conditioning process in detail.

    3 Preparation for RF conditioning

    The power source is an RF unit based on a solid-state modulator, and a 50 MW S-band high-power pulse klystron (E3730A, Toshiba, Tokyo, Japan) was used for the RF conditioning. The klystron facility is shown in Fig. 7a,and it is located in the technology corridor. After passing the 6 dB power divider, the power allocated from the klystron is approximately 12 MW. We measured the input power at the entrance of the accelerating structure with a power meter and calculated the power attenuation of the fixed attenuator, directional coupler, and cable. The input power of the accelerating structure was calculated to be approximately 9.9 MW. A detailed calculation process is presented in Sect. 2. This result exceeds our experimental goal. The purpose of RF conditioning is to achieve the nominal parameters in terms of the input power at the input coupler (9.15 MW, corresponding to an accelerating gradient of 50 MV/m), the repetition rate (50 Hz), and pulse length (3 μs, corresponding to one filling time).

    Fig. 7 (Color online) a 50 MW S-band high-power pulse klystron.b Input and output characteristic curves of the klystron

    The relevant measurement personnel of the SXFEL calibrated the power measurement values. The following section briefly introduces the performance of the klystron.The input and output characteristic curves of the klystron are shown in Fig. 7b. The maximum value of the charging voltage of the modulator should not exceed 45 kV. All the data in this figure were measured by the relevant technical personnel of the SXFEL.The charging voltage value of the modulator is manually changed, and the output pulse and output voltage values of the klystron are recorded on the control panel.Figure 7b is a critical reference basis for RF conditioning. In the RF conditioning control interface in the central control room, the charging voltage of the modulator can be changed to complete the high-power test.The output power data of the klystron can be obtained according to the input and output characteristic curve of the klystron.Finally, the accelerating gradient is calculated by calculating the power that is input to the accelerating structure after attenuation. The charge voltage of the modulator gradually increases during the conditioning process; therefore, the power that is input to the accelerating structure increases,and the accelerating gradient also increases accordingly.

    The input signal, reflected signal, and pickup signal of the RF pulse are shown in Fig. 8.In Fig. 8,the red and blue lines indicate the voltage and RF power, respectively. In the graph of the reflected signal,the blue line is lower than the red line, which indicates that the accelerating structure is under-coupled. If the accelerating structure is in a critically coupled state, the red and blue lines would coincide,and if the accelerating structure is in an over-coupled state,the blue line would be higher than the red line. It can be observed from the shape of the reflected signal that this accelerating structure is under-coupled. This result is consistent with our previous microwave measurement results in Ref. [21].

    4 Process and results of RF conditioning

    The entire conditioning procedure is semi-automatic,and the modulator charge voltage may be switched off(which is what we often call the breakdown phenomenon)by an artificial breakdown made by local operators or a broken ion pump vacuum, which causes the pressure to exceed the threshold corresponding to 1 × 10-8mbar or causes the monitored reflection signal of the klystron to be more than twice the normal value. In our conditioning procedure, after a breakdown phenomenon occurs, the conditioning procedure will activate the automatic protection mechanism, whether the phenomenon is due to a vacuum drop, caused by breakdown, or the sparking phenomenon, caused by exceeding the reflected power threshold. First, the modulator’s high voltage is automatically decreased by 2 kV, the corresponding input power will decrease, and the high voltage of the modulator is slowly increased by 0.01 kV every 4 s on average. Thisprocess lasts for several minutes. After 30 s of stable operation under the current input power, the RF high voltage is increased by 0.01 kV. The entire RF conditioning process lasts for approximately 240 h.

    Fig. 8 (Color online) The shapes of different RF pulses for a the input signal, b the reflected signal, and c the pickup signal from the cavity.The red and blue lines indicate the voltage and RF power,respectively

    Fig. 9 (Color online) Pulse length

    We initially set the repetition frequency to 50 Hz and the RF pulse length to 2.5 μs during conditioning. Subsequently,the charging voltage of the modulator is gradually increased.After approximately 200 h,the charging voltage of the modulator rises to 40 kV. Before increasing the RF pulse length, we ran the system stably for 10 h with the maximum charge voltage of the modulator, and no breakdown event occurred. After increasing the pulse length to 3 μs (as shown in Fig. 9), the charging voltage of the modulator increases to 40 kV after a short time.

    Fig. 10 Trend of modulator high voltage change

    Fig. 11 Pulse current of the klystron measured by the oscilloscope

    Fig. 12 (Color online) Snapshot of the actual input power value of the accelerator tube measured by the power meter

    Table 1 Attenuation of each component

    The high-voltage signal of the modulator is recorded,and the high-voltage trend during the entire high-power experiment is shown in Fig. 10. As can be observed in the figure, in the early stage of the high-power testing, the modulator’s high voltage rises slowly,and breakdowns are frequent. This occurs because the proton-accelerating structure has recently left the factory, which indicates that there may be more burrs or residual gas inside the cavity.The intermediate interruption is owing to the installation of other equipment by the SXFEL equipment,and not because the equipment is malfunctioning. The high voltage rise of the modulator in the later stage is faster because the cavity inner surface burr decreases,and the internal residual gas is reduced after the previous RF conditioning; this improves the working environment in the cavity. Therefore, the number of breakdowns in this stage is reduced.Finally,the system runs stably while maintaining a very low breakdown rate.

    The pulse current of the klystron measured by the oscilloscope is shown in Fig. 11. In this figure, we can observe that the measured minimum value is - 35 V, and the klystron pulse current is calculated to be 350 A. The output pulse current value of the klystron is verified. This value is obtained based on the current transformer monitor with a nominal ratio of 10:1; the related theory can be found in Ref. [28].

    Fig. 13 (Color online) Four signals collected by an oscilloscope in the case of an approximate input power of 9.9 MW and pulse length of 3 μs pulse length. a Input and reflected power signals of the accelerating structure,output power signal of the klystron,and pickup signal. b Input and reflected power signals of the accelerating structure, dark current signal, and pickup signal

    We use a power meter to measure the input power at the entrance of the proton-accelerating structure.The measurement result is shown in Fig. 12,and the value is 14.15 dBm.The attenuations of the directional coupling cavity, fixed attenuator,and cable are shown in Table 1.The input power at the entrance of the proton-accelerating structure is Pin.From the formula 10log10Pin(mW)=^P(dBm), Pincan be calculated to be 9.9 MW.From this value,the input power atthe entrance of the proton-accelerating structure is obtained as 9.9 MW.

    Table 2 Design and measurement values of S-band standing wave accelerating structure

    The incident power at the entrance of the proton-accelerating structure is Pin, part of which enters the accelerating structure, Pcav; some power is reflected as reverse power, Prev. The incident power is the sum of the cavity and reverse power, Pin=Pcav+Prev. The incident power enters the accelerating cavity through a coupling hole with the waveguide cavity coupling parameter βc. The reverse power is given by [29]

    The effective input power entering the acceleration cavity can be calculated as Pcav=Pin-Prev=7.7 MW.According to Eq. (3), the actual accelerating gradient can be calculated to be 46 MV/m.

    We increased the input power to a maximum of 9.9 MW, which theoretically means that the accelerating gradient of the S-band accelerating structure has reached the expected goal. However, owing to processing errors,the SWR and reflected power are too large, which eventually causes the accelerating gradient to be slightly different from the theoretical value. This indicates that the proposed high-gradient accelerating structure for a compact proton therapy facility was successfully manufactured,and the measurement was completed. This prototype provides valuable experience for future proton therapy facilities. Details on the process and results of the high-power test will be presented next.

    For an input power of approximately 9.9 MW and a pulse length of 3 μs,we use an oscilloscope to measure the input and reflected power signals of the accelerating structure, the output power signal of the klystron and the pickup signal. The four signals are shown in Fig. 13a. It can also be observed from the reflected signal measured in Fig. 13a that the accelerating structure is under-coupled.An oscilloscope was used to measure the dark current onsite; however, it can be observed from Fig. 13b that no dark current signal is detected. This is because the beam aperture of the accelerating structure is too small to observe the dark current.

    To date,all the tests have been completed.According to the tested S-band SW accelerating structure results, all the design goals have been achieved. Table 2 lists the design targets and measured values. The measurement result shows that the effective input power is 7.7 MW.According to the above formula (3), the actual accelerating gradient can be calculated as 46 MV/m.

    5 Discussion

    Initially,an accelerating structure with a single periodic operating frequency in the π-mode was chosen because it corresponded to the highest shunt impedance. However, it was observed that the operation of the single periodic accelerating structure was unstable. If the mode interval is too small, it would cause jumps between adjacent modes.We will subsequently conduct an investigation on the biperiodic SW accelerating structure for proton therapy facilities, including the side-coupled SW accelerating structure and the axis-coupled accelerating structure. The Institute of Nuclear Science and Technology, Sichuan University, China, proposed the construction of a 3 GHz side-coupled cavity linear accelerator(SCL)to reaccelerate the 26 MeV proton beam extracted from the CS-30 cyclotron to 120 MeV [30]. However, the ultimate design goal of the accelerating structure that we designed is to increase the energy of the proton beam to 350 MeV to meet the proton energy required by proton CT.

    6 Conclusion

    A prototype of the S-band accel erating structure was successfully designed,fabricated,and tuned.Subsequently,the high-power test platform was set up in the tunnel of the SXFEL. After RF conditioning for approximately 240 h,the prototype structure was completed, and a high-power test was performed with an input power of 9.9 MW,under the condition of a 50 Hz repetition rate and 3 μs pulse length. The results of the high-power test mostly verified the design of the S-band high-gradient technology. More importantly, the results demonstrated the strong feasibility of developing compact proton therapy facilities based on the high-gradient linac.The development of an S-band SW accelerating structure prototype provides valuable experience and a solid foundation for the research, development,and construction of critical technologies for future proton therapy facilities.

    Author contributionsAll authors contributed to the study conception and design. Material preparation, data collection, and analysis were performed by Yu Zhang, Wen-Cheng Fang, Xiao-Xia Huang,Jian-Hao Tan, Shao-Peng Zhong, Cheng-Cheng Xiao, Jun-Qiang Zhang,Cheng Wang,Yu-Qing Yang,and Zhen-Tang Zhao.The first draft of the manuscript was written by Yu Zhang and all authors commented on previous versions of the manuscript. All authors read and approved the final manuscript.

    高清黄色对白视频在线免费看| 亚洲成人手机| 每晚都被弄得嗷嗷叫到高潮| 久久久久精品人妻al黑| 18在线观看网站| 婷婷丁香在线五月| 下体分泌物呈黄色| 国产又色又爽无遮挡免费看| 51午夜福利影视在线观看| 亚洲伊人久久精品综合| 国产在视频线精品| 精品亚洲成国产av| 成在线人永久免费视频| 99久久精品国产亚洲精品| www.精华液| 亚洲综合色网址| 久久影院123| 国产老妇伦熟女老妇高清| 国产日韩欧美亚洲二区| 视频区欧美日本亚洲| 中国美女看黄片| 久久久久久亚洲精品国产蜜桃av| av电影中文网址| 亚洲国产成人一精品久久久| 一个人免费看片子| 80岁老熟妇乱子伦牲交| 80岁老熟妇乱子伦牲交| 久久久水蜜桃国产精品网| 久久国产精品大桥未久av| 国产成人系列免费观看| 超碰成人久久| 欧美精品亚洲一区二区| 亚洲精品成人av观看孕妇| 国产精品久久久久久精品古装| 久久精品国产亚洲av香蕉五月 | 精品亚洲成国产av| 欧美老熟妇乱子伦牲交| 亚洲精品国产精品久久久不卡| 18在线观看网站| 99热网站在线观看| 99精品久久久久人妻精品| 欧美成人免费av一区二区三区 | 欧美日韩福利视频一区二区| 黄片大片在线免费观看| 欧美成狂野欧美在线观看| 中文字幕人妻熟女乱码| 亚洲精品国产色婷婷电影| 777米奇影视久久| 黄色a级毛片大全视频| 搡老熟女国产l中国老女人| 岛国在线观看网站| 成人特级黄色片久久久久久久 | 黄片小视频在线播放| 午夜福利在线免费观看网站| 一夜夜www| 国产精品亚洲av一区麻豆| 欧美乱码精品一区二区三区| 亚洲国产看品久久| 国产黄频视频在线观看| 黄色毛片三级朝国网站| 精品亚洲成国产av| 亚洲精华国产精华精| 欧美中文综合在线视频| 丰满少妇做爰视频| 中文字幕制服av| 久久毛片免费看一区二区三区| 19禁男女啪啪无遮挡网站| 国产极品粉嫩免费观看在线| 一区二区三区国产精品乱码| 精品亚洲乱码少妇综合久久| 少妇裸体淫交视频免费看高清 | 满18在线观看网站| 国产aⅴ精品一区二区三区波| 欧美成人午夜精品| 国产精品.久久久| 人人妻人人澡人人看| 99久久99久久久精品蜜桃| 新久久久久国产一级毛片| 国产有黄有色有爽视频| 夫妻午夜视频| 欧美亚洲日本最大视频资源| 午夜福利在线免费观看网站| 亚洲精品国产一区二区精华液| 亚洲午夜精品一区,二区,三区| 美女午夜性视频免费| 欧美精品人与动牲交sv欧美| 我的亚洲天堂| 免费看a级黄色片| 男女之事视频高清在线观看| 国产高清国产精品国产三级| 999久久久国产精品视频| 国产免费现黄频在线看| 亚洲国产欧美在线一区| 丰满饥渴人妻一区二区三| 中文字幕人妻丝袜一区二区| 亚洲成a人片在线一区二区| 久久精品亚洲精品国产色婷小说| 国产av一区二区精品久久| 久久久久国内视频| 免费一级毛片在线播放高清视频 | 久久久精品94久久精品| 每晚都被弄得嗷嗷叫到高潮| 免费观看av网站的网址| 大片电影免费在线观看免费| 精品一区二区三区视频在线观看免费 | 午夜成年电影在线免费观看| 亚洲精品国产区一区二| 国产av精品麻豆| 久久久久国产一级毛片高清牌| 国产精品久久久久久精品电影小说| 另类亚洲欧美激情| bbb黄色大片| 搡老熟女国产l中国老女人| 欧美在线黄色| 免费观看a级毛片全部| xxxhd国产人妻xxx| 亚洲专区国产一区二区| 免费在线观看黄色视频的| 日韩有码中文字幕| 高清av免费在线| 国产精品99久久99久久久不卡| 国产91精品成人一区二区三区 | 又紧又爽又黄一区二区| 人人妻人人澡人人爽人人夜夜| 一本—道久久a久久精品蜜桃钙片| www日本在线高清视频| 午夜福利乱码中文字幕| 国产av又大| 丝袜美足系列| 日韩欧美一区视频在线观看| 精品熟女少妇八av免费久了| 搡老熟女国产l中国老女人| 亚洲国产看品久久| 欧美日韩福利视频一区二区| 激情在线观看视频在线高清 | 欧美成狂野欧美在线观看| 精品高清国产在线一区| 人妻 亚洲 视频| 亚洲色图 男人天堂 中文字幕| 亚洲av片天天在线观看| 国产免费av片在线观看野外av| 香蕉国产在线看| 色综合婷婷激情| tube8黄色片| 一级,二级,三级黄色视频| 色视频在线一区二区三区| 嫩草影视91久久| 国产精品av久久久久免费| 亚洲精品国产一区二区精华液| 建设人人有责人人尽责人人享有的| 亚洲第一青青草原| 1024视频免费在线观看| 中国美女看黄片| 啦啦啦 在线观看视频| 人成视频在线观看免费观看| 日本wwww免费看| 一级毛片电影观看| 亚洲国产欧美日韩在线播放| 国产成人精品久久二区二区91| 欧美大码av| 嫩草影视91久久| 国产成人精品久久二区二区91| 欧美久久黑人一区二区| 欧美av亚洲av综合av国产av| 69精品国产乱码久久久| 成年版毛片免费区| 久久久国产精品麻豆| 欧美一级毛片孕妇| 超碰97精品在线观看| 日韩视频一区二区在线观看| 久久人妻av系列| 丝瓜视频免费看黄片| 精品国产乱码久久久久久男人| 日本一区二区免费在线视频| 国产野战对白在线观看| 每晚都被弄得嗷嗷叫到高潮| 久久中文看片网| 国产成人精品无人区| 亚洲免费av在线视频| 国产1区2区3区精品| 久久久久久免费高清国产稀缺| 国产成人系列免费观看| 18禁美女被吸乳视频| 欧美黄色淫秽网站| 亚洲色图综合在线观看| 亚洲精品国产区一区二| 高清视频免费观看一区二区| 午夜福利在线观看吧| 多毛熟女@视频| 天天躁日日躁夜夜躁夜夜| 国产欧美日韩一区二区三区在线| 无遮挡黄片免费观看| www.999成人在线观看| 超色免费av| 久久国产精品人妻蜜桃| 一级黄色大片毛片| 久久久精品免费免费高清| 18禁美女被吸乳视频| 亚洲avbb在线观看| 欧美日韩国产mv在线观看视频| 12—13女人毛片做爰片一| 18禁国产床啪视频网站| 久久精品91无色码中文字幕| 午夜成年电影在线免费观看| a级毛片黄视频| 国产日韩欧美在线精品| 免费女性裸体啪啪无遮挡网站| 日本wwww免费看| 国产成人精品久久二区二区免费| 久久久久久人人人人人| 久久精品国产亚洲av香蕉五月 | 亚洲第一青青草原| 亚洲五月婷婷丁香| 亚洲欧美日韩高清在线视频 | 嫁个100分男人电影在线观看| 欧美中文综合在线视频| 天堂俺去俺来也www色官网| av有码第一页| 中文字幕色久视频| 国产精品一区二区精品视频观看| 日韩三级视频一区二区三区| 操美女的视频在线观看| 亚洲色图av天堂| 午夜福利,免费看| 欧美成狂野欧美在线观看| 黄色怎么调成土黄色| 免费看a级黄色片| 亚洲国产欧美一区二区综合| 午夜免费鲁丝| 欧美+亚洲+日韩+国产| 久久人妻福利社区极品人妻图片| 成人18禁在线播放| 国产精品熟女久久久久浪| 80岁老熟妇乱子伦牲交| 99热国产这里只有精品6| 亚洲自偷自拍图片 自拍| 久热这里只有精品99| 国产av又大| 成人永久免费在线观看视频 | 狂野欧美激情性xxxx| 欧美精品亚洲一区二区| 亚洲精品美女久久久久99蜜臀| 女人精品久久久久毛片| 一本久久精品| 精品国产一区二区三区四区第35| 免费看十八禁软件| 久久人人爽av亚洲精品天堂| 999久久久国产精品视频| 香蕉丝袜av| 人妻久久中文字幕网| 亚洲午夜精品一区,二区,三区| 在线观看一区二区三区激情| 国产成人精品无人区| 婷婷丁香在线五月| 亚洲色图av天堂| 1024香蕉在线观看| 欧美成狂野欧美在线观看| 久久99热这里只频精品6学生| 欧美日韩亚洲高清精品| 色94色欧美一区二区| 国产成人欧美在线观看 | 日韩免费高清中文字幕av| 国产精品香港三级国产av潘金莲| 母亲3免费完整高清在线观看| 国产亚洲欧美精品永久| 黄色视频在线播放观看不卡| 人妻 亚洲 视频| 亚洲国产精品一区二区三区在线| 老司机福利观看| 亚洲一区中文字幕在线| 久久久欧美国产精品| 亚洲熟女精品中文字幕| 啦啦啦视频在线资源免费观看| 日本五十路高清| 男女免费视频国产| 在线播放国产精品三级| 超碰97精品在线观看| 久久精品人人爽人人爽视色| 成人国产av品久久久| 午夜福利乱码中文字幕| 国产精品成人在线| 真人做人爱边吃奶动态| 日韩欧美三级三区| 国产片内射在线| 国产精品久久久av美女十八| 色婷婷av一区二区三区视频| 欧美+亚洲+日韩+国产| tocl精华| 另类精品久久| 一夜夜www| 国产aⅴ精品一区二区三区波| 欧美日韩黄片免| 成年人免费黄色播放视频| 欧美日韩一级在线毛片| 午夜91福利影院| 国产xxxxx性猛交| 嫁个100分男人电影在线观看| 91精品国产国语对白视频| 少妇精品久久久久久久| 高清毛片免费观看视频网站 | 91麻豆精品激情在线观看国产 | 久久精品熟女亚洲av麻豆精品| 又大又爽又粗| 丁香六月欧美| 国产日韩欧美亚洲二区| 国产av精品麻豆| 国产亚洲欧美精品永久| 999久久久精品免费观看国产| 中文字幕人妻丝袜制服| 国产高清国产精品国产三级| 丁香六月天网| 日本av免费视频播放| 不卡av一区二区三区| 18禁国产床啪视频网站| 18禁裸乳无遮挡动漫免费视频| 日韩一区二区三区影片| 国产精品av久久久久免费| 亚洲欧洲日产国产| 国产成人欧美| 人人妻人人澡人人爽人人夜夜| 精品欧美一区二区三区在线| 在线观看66精品国产| 亚洲av美国av| 在线av久久热| 国产精品一区二区精品视频观看| 99在线人妻在线中文字幕 | 天堂中文最新版在线下载| 精品午夜福利视频在线观看一区 | 欧美+亚洲+日韩+国产| 99精品欧美一区二区三区四区| 一级片免费观看大全| 一级毛片电影观看| 黑人操中国人逼视频| 午夜两性在线视频| 国产xxxxx性猛交| 性高湖久久久久久久久免费观看| 1024视频免费在线观看| 一区二区av电影网| 一级毛片电影观看| 热99国产精品久久久久久7| 黄色 视频免费看| 啦啦啦免费观看视频1| 十八禁网站网址无遮挡| 精品免费久久久久久久清纯 | 欧美日韩视频精品一区| 欧美性长视频在线观看| 国产三级黄色录像| 亚洲欧洲精品一区二区精品久久久| videos熟女内射| 久久午夜亚洲精品久久| 久久精品亚洲熟妇少妇任你| 日韩三级视频一区二区三区| 亚洲精品国产色婷婷电影| www.自偷自拍.com| a在线观看视频网站| 精品国产国语对白av| 国产91精品成人一区二区三区 | 欧美 日韩 精品 国产| 狠狠婷婷综合久久久久久88av| 亚洲欧洲精品一区二区精品久久久| 婷婷成人精品国产| 免费高清在线观看日韩| 老司机午夜十八禁免费视频| 成年版毛片免费区| 美女高潮到喷水免费观看| 黑人猛操日本美女一级片| 成年动漫av网址| 国产亚洲欧美精品永久| 欧美变态另类bdsm刘玥| 久久影院123| 黄片小视频在线播放| 两个人看的免费小视频| 久久久水蜜桃国产精品网| av国产精品久久久久影院| 亚洲精品美女久久久久99蜜臀| 久久99热这里只频精品6学生| 一二三四社区在线视频社区8| 天堂8中文在线网| 天堂动漫精品| 午夜福利视频精品| 叶爱在线成人免费视频播放| 在线十欧美十亚洲十日本专区| 国产福利在线免费观看视频| 久久国产精品影院| 超色免费av| 久久久久久免费高清国产稀缺| 欧美成狂野欧美在线观看| 老汉色∧v一级毛片| 一区在线观看完整版| 欧美 日韩 精品 国产| 国产精品美女特级片免费视频播放器 | 国产有黄有色有爽视频| 国产老妇伦熟女老妇高清| 久久精品国产亚洲av高清一级| 女性被躁到高潮视频| 国产成人精品无人区| 97人妻天天添夜夜摸| 色精品久久人妻99蜜桃| 久久久精品区二区三区| 一进一出抽搐动态| 91成人精品电影| 婷婷丁香在线五月| 国产欧美日韩一区二区三| 久久ye,这里只有精品| 亚洲,欧美精品.| 日本av手机在线免费观看| 国产精品九九99| 精品卡一卡二卡四卡免费| 丝袜美腿诱惑在线| 五月天丁香电影| 中亚洲国语对白在线视频| 国产日韩欧美视频二区| 日韩有码中文字幕| 一个人免费看片子| 久久久精品94久久精品| 91精品三级在线观看| 久久 成人 亚洲| 一个人免费在线观看的高清视频| 免费日韩欧美在线观看| 国产黄色免费在线视频| 欧美日韩亚洲高清精品| 成人国产av品久久久| 可以免费在线观看a视频的电影网站| 99re在线观看精品视频| 亚洲精品国产色婷婷电影| 日本精品一区二区三区蜜桃| 色综合婷婷激情| 免费观看av网站的网址| 一级a爱视频在线免费观看| 一区二区三区精品91| 女人被躁到高潮嗷嗷叫费观| 国产区一区二久久| 80岁老熟妇乱子伦牲交| 精品少妇久久久久久888优播| 久久久国产精品麻豆| 久久精品亚洲精品国产色婷小说| 母亲3免费完整高清在线观看| 中文亚洲av片在线观看爽 | 成人精品一区二区免费| 亚洲七黄色美女视频| 欧美精品亚洲一区二区| 久久九九热精品免费| 自线自在国产av| 可以免费在线观看a视频的电影网站| 午夜福利乱码中文字幕| 色视频在线一区二区三区| 国产免费现黄频在线看| 在线 av 中文字幕| 两个人免费观看高清视频| 中文字幕色久视频| 1024香蕉在线观看| 亚洲中文日韩欧美视频| 最近最新中文字幕大全免费视频| 欧美成人午夜精品| 精品国产一区二区久久| 国产淫语在线视频| 青青草视频在线视频观看| av有码第一页| 欧美精品啪啪一区二区三区| 国产黄频视频在线观看| 十八禁高潮呻吟视频| 大片电影免费在线观看免费| 日本精品一区二区三区蜜桃| 亚洲欧美一区二区三区久久| 天天操日日干夜夜撸| 亚洲三区欧美一区| 精品久久久久久电影网| 天天影视国产精品| 18在线观看网站| 欧美日韩黄片免| 老司机午夜十八禁免费视频| 日韩欧美一区视频在线观看| 国产av国产精品国产| 免费观看人在逋| 麻豆成人av在线观看| 亚洲五月婷婷丁香| 97在线人人人人妻| 一级毛片电影观看| 女人久久www免费人成看片| 日韩免费高清中文字幕av| 久久中文字幕一级| 亚洲av片天天在线观看| 成年人免费黄色播放视频| 搡老乐熟女国产| 日韩精品免费视频一区二区三区| 可以免费在线观看a视频的电影网站| 韩国精品一区二区三区| 午夜精品久久久久久毛片777| 国产亚洲一区二区精品| 精品熟女少妇八av免费久了| 国产精品麻豆人妻色哟哟久久| 中文字幕最新亚洲高清| 91成年电影在线观看| 色老头精品视频在线观看| 欧美成狂野欧美在线观看| 国产日韩欧美亚洲二区| av天堂久久9| 亚洲av电影在线进入| 97人妻天天添夜夜摸| 美女主播在线视频| 精品国产一区二区三区四区第35| 精品国产乱子伦一区二区三区| 亚洲性夜色夜夜综合| 欧美激情久久久久久爽电影 | videos熟女内射| 在线播放国产精品三级| 亚洲av电影在线进入| 欧美精品亚洲一区二区| 国产成人av教育| 午夜福利乱码中文字幕| 在线观看免费高清a一片| 天天操日日干夜夜撸| 成人av一区二区三区在线看| 成人三级做爰电影| 岛国毛片在线播放| 手机成人av网站| 亚洲天堂av无毛| 波多野结衣av一区二区av| 最新美女视频免费是黄的| 亚洲性夜色夜夜综合| 国产亚洲av高清不卡| 菩萨蛮人人尽说江南好唐韦庄| 少妇被粗大的猛进出69影院| 亚洲国产欧美网| 亚洲精品自拍成人| 午夜免费成人在线视频| 制服人妻中文乱码| 99精国产麻豆久久婷婷| 久久国产精品影院| 日本撒尿小便嘘嘘汇集6| 亚洲视频免费观看视频| 一区二区三区乱码不卡18| 视频在线观看一区二区三区| 日韩欧美国产一区二区入口| 亚洲精品国产色婷婷电影| 久久国产精品大桥未久av| 在线观看免费视频日本深夜| 女人精品久久久久毛片| 欧美黑人欧美精品刺激| 在线十欧美十亚洲十日本专区| 波多野结衣一区麻豆| 美女午夜性视频免费| 啦啦啦在线免费观看视频4| 男女下面插进去视频免费观看| 中文字幕制服av| √禁漫天堂资源中文www| 国产成人精品在线电影| 女同久久另类99精品国产91| 久久久久久亚洲精品国产蜜桃av| 欧美黑人欧美精品刺激| 在线 av 中文字幕| 国产精品美女特级片免费视频播放器 | 黄色片一级片一级黄色片| 精品人妻在线不人妻| 大片免费播放器 马上看| 一级,二级,三级黄色视频| 人人妻,人人澡人人爽秒播| 久久久久精品国产欧美久久久| 黄色视频,在线免费观看| 亚洲色图av天堂| 亚洲人成伊人成综合网2020| 国产精品免费一区二区三区在线 | 国产精品99久久99久久久不卡| 久久久精品国产亚洲av高清涩受| 91麻豆精品激情在线观看国产 | 男男h啪啪无遮挡| 国产黄色免费在线视频| 波多野结衣av一区二区av| 一级毛片女人18水好多| 母亲3免费完整高清在线观看| av线在线观看网站| 国产精品一区二区免费欧美| 大型av网站在线播放| 大型黄色视频在线免费观看| 成人影院久久| 国产伦人伦偷精品视频| 在线亚洲精品国产二区图片欧美| 夫妻午夜视频| 欧美日韩亚洲综合一区二区三区_| 亚洲精品国产区一区二| 国产精品免费大片| 99精国产麻豆久久婷婷| 婷婷成人精品国产| 法律面前人人平等表现在哪些方面| 午夜福利在线免费观看网站| 精品国产超薄肉色丝袜足j| av又黄又爽大尺度在线免费看| 国产一区二区激情短视频| 国产视频一区二区在线看| 亚洲第一欧美日韩一区二区三区 | 在线播放国产精品三级| 久久久水蜜桃国产精品网| 91九色精品人成在线观看| 人成视频在线观看免费观看| 亚洲欧美激情在线| 精品熟女少妇八av免费久了| 91成人精品电影| 免费久久久久久久精品成人欧美视频| 日韩欧美一区视频在线观看| 欧美久久黑人一区二区| 亚洲精品国产精品久久久不卡| 美女国产高潮福利片在线看| 亚洲伊人久久精品综合| 99国产精品免费福利视频| 水蜜桃什么品种好| 夜夜夜夜夜久久久久| 欧美日本中文国产一区发布| 在线观看66精品国产| 可以免费在线观看a视频的电影网站| 欧美激情极品国产一区二区三区| 国产亚洲av高清不卡| 777米奇影视久久| 露出奶头的视频| 久久婷婷成人综合色麻豆| 黄色丝袜av网址大全|