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

    Investigation of self-generated magnetic field and dynamics of a pulsed plasma flow

    2022-06-01 07:55:56AigerimTAZHENMerlanDOSBOLAYEVandTlekkabulRAMAZANOV
    Plasma Science and Technology 2022年5期

    Aigerim TAZHEN, Merlan DOSBOLAYEV and Tlekkabul RAMAZANOV

    Institute of Experimental and Theoretical Physics,Al-Farabi Kazakh National University,Almaty 050040,Kazakhstan

    Abstract Due to the growing interest in studying the compression and disruption of the plasma filament in magnetic fusion devices and Z-pinches,this work may be important for new developments in the field of controlled thermonuclear fusion.Recently,on a coaxial plasma accelerator,we managed to obtain the relatively long-lived (~300 μs) plasma filaments with its self-magnetic field.This was achieved after modification of the experimental setup by using high-capacitive and lowinductive energy storage capacitor banks, as well as electrical cables with low reactive impedance.Furthermore,we were able to avoid the reverse reflection of the plasma flux from the end of the plasma accelerator by installing a special plasma-absorbing target.Thus, these constructive changes of the experimental setup allowed us to investigate the physical properties of the plasma filament by using the comprehensive diagnostics including Rogowski coil,magnetic probes, and Faraday cup.As a result, such important plasma parameters as density of ions and temperature of electrons in plasma flux, time dependent plasma filament’s azimuthal magnetic field were measured in discharge gap and at a distance of 23.5 cm from the tip of the cathode.In addition, the current oscillograms and I-V characteristics of the plasma accelerator were obtained.In the experiments, we also observed the charge separation during the acceleration of plasma flow via oscillograms of electron and ion beam currents.

    Keywords:coaxial plasma accelerator,pulsed plasma flow,self magnetic field,electron density,electron temperature, Faraday cup, magnetic probe

    1.Introduction

    It is well known that during edge localized mode instabilities in magnetic fusion devices, the disruption of the plasma filament and the production of the heat and particle fluxes occur on the wall of the reactors.As a result, the plasmafacing material surface is damaged and dust particles are created.All of these events negatively affect the safety and normal operation of the reactor.As an example,dust ejection into the plasma core volume leads to the contamination of the discharge and may contribute to tritium safety issues [1].

    In studies [2-5], the disruption and plasma-wall interaction simulation experiments were performed at the plasma accelerators.As a potential candidate material for the first wall,the graphite was chosen in[2,3].The graphite plate was irradiated with the high-energy plasma flux generated in a coaxial plasma accelerator.In work [4] a quasi-stationary pulsed plasma accelerator was used to provide high heat plasma flux, which is expected during operation of an international thermonuclear experimental reactor,and to study the damage of divertor plasma-facing components with a combination of the pulsed and quasi-stationary heat fluxes.In that work the tungsten target was used as a candidate material of the divertor.In work [5] a pulsed coaxial plasma gun was assembled to investigate the dependence plasma flow structure on initial operating working gas pressure at‘continuously filled’regime.In studies[6-8]new design and construction of pulsed plasma accelerators intended for creation of dense pinch plasma focus discharges were presented.Moreover, in this work, the impact of the electrode system on the neutron yield, x-ray radiation in a dense pinch plasma focus was investigated.

    As mentioned above, one of the applications of plasma accelerators is the experimental simulation of the plasma filament compression and disruption in magnetic fusion devices and Z-pinches [9, 10].For these purposes, recently we have modified the coaxial plasma accelerator by installing more capacitive and low-inductive energy storage capacitors and also by using electrical coupling cables with low reactive impedance.In addition, the avoidance of reverse reflected plasma flow from the end of the accelerator is taken into account.It was achieved by using a special plasma-absorbing target installed inside the vacuum chamber.This constructive modification allows obtaining plasma flow with its selfmagnetic field and the relatively long body filament.Those changes give the opportunity to experimentally simulate the processes occurring in the plasma filament.The main purpose of this work is to conduct complex plasma diagnostics for studying the physical properties and plasma flow formation dynamics in the channel of modified plasma accelerator.There is a wide spectrum of diagnostic methods for the investigation of the plasma filament and its basic parameters in CPA at extreme conditions (at high temperatures and plasma densities, at high frequency oscillations created in an accelerator by pinch discharges).Of these, some of the frequently used electrical signal measurement-based techniques for plasma diagnostics such as Rogowski coil, Faraday cup,and magnetic probes, were utilized.Furthermore, high-speed imaging techniques, for visual observation of the accelerated plasma, were used.

    The work novelty is summarized as follows: in comparison with related former researches, this work considers the most important physical properties of plasma flow and the physical processes, such as charge separation (electrons and ions) in plasma flow with FC, compression and disruption of plasma filaments using frames synchronized with magnetic probe signals.

    2.Experimental setup

    A principle scheme of the coaxial plasma accelerator is shown in figure 1.An experimental setup consists of two concentric(coaxial) copper electrodes: an inner anode and an outer cathode.Anode and cathode are isolated by a dielectric material.The operational principle of experimental setup is based on the acceleration of plasma generated by ionization of gas filling in interelectrode space, when a high-voltage is supplied to the electrodes.The acceleration takes place due to the J×B force acting on plasma.A vacuum arrester was used for the transfer of capacitor-stored energy to the gas between two electrodes.The installation design and its operational principle are described in detail elsewhere[11-14].

    The power system of the CPA includes a low-inductive capacitor banks with five parallel-connected capacitors of the KPIMK-8-288 type (288 μF).The total capacitance is 1.44 mF and maximum energy stored on the capacitors at a maximum charging operating voltage of 6 kV is W = 26 kJ.

    Figure 1.A principle scheme of the coaxial plasma accelerator.

    3.Diagnostics techniques and apparatus

    The current flowing through the electric circuit of the CPA was measured by a Rogowski coil [15] with the ~75 mm diameter loop and L = 27.5 μH inductance.

    The magnitudes of magnetic fields generated by pulsed plasma flow in space between two concentric electrodes and at the out of an accelerator were measured using two magnetic probes (No.1 and No.2 in figure 1).The dimension and characteristics of each probe are as follows: number of turns n = 7, loop radius r = 1.1 mm, loop cross-sectional area A = 3.8 mm2, time constant τ = L/R = 2.3 ns [16].The winding ends of magnetic probes were twisted and brought out.Then they were connected to the BNC cable attached to the end of glass tube.The impedance of this cable is 50 Ω.The BNC cable was used to connect the twisted pair terminals of magnetic probes to the oscilloscope through a passive RC integrator circuit (RC = 100 μs).

    The formation, dynamics and development of pulsed plasma flow inside the CPA vacuum chamber were captured by the high-speed Phantom VEO710S CMOS camera [17].The plasma flow velocity was measured by analyzing the obtained images.The frequency and exposure time of the high-speed camera were 280 000 frames s-1and 3.1 μs respectively.

    The Faraday cup with 160 μm orifice [18] was used to measure the electron and ion currents,energy and densities of charged particles of plasma flow.

    All signals were recorded on the digital storage oscilloscope LeCroy WaveJet 354 A, 500 MHz/4-channel.Also,Berkeley Nucleonics pulse generator (model 577) with an external synchronizing plasma-created signal was used to trigger the scope and CMOS camera in the external trigger mode.

    All the experiments were carried out at voltages varied between 2 and 6 kV and at a plasma-forming gas pressures varied between 0.03 and 11 Torr.Two types of gases such as hydrogen, residual air were used as a plasma source.An electrical discharge was initiated by a high-voltage vacuum gap with a residual pressure of ~10-2Torr.

    4.Results and discussion

    This section presents the results of plasma diagnostics using the methods described above:(i)Rogowski coil,(ii)magnetic probe and Faraday cup, (iii) high-speed CMOS camera.

    For controlling the operation of plasma accelerators, it is important to know the value of the current flowing through its driving circuit.In this regard,the typical current oscillograms and I-V characteristic of the CPA,obtained for the two types of working gas, shown in figures 2(a) and (b).These discharge current oscillograms are in well agreement with the oscillograms presented in other studies[19,20]and represent a rapidly damping waveform.This indicates a low inductance of the system and demonstrates the fact of effective transfer of the energy stored in the capacitors.The use of air as the working gas for plasma generation is most convenient for conducting preliminary experiments.

    In this work,with a 5 kV charging voltage and 60 mTorr gas pressure, we could reach the maximum current values of 80 kA for the air and 65 kA for the hydrogen, respectively.From obtained I-V characteristics of the CPA (figure 2(b))one can see that the maximum current linearly increases with increasing the charging voltage supplied to the capacitor banks.We explain this by the insignificant contribution of the pulsed plasma active resistance and the absence of current skipping.In this case,it can be considered that the maximum current Imaxflowing through the CPA driving circuit is proportional to the applied voltage U.A relevant I-V characteristic can be also found in [21].

    Figure 2.(a) The typical current waveforms and (b) the current dependence on the voltages supplied to the capacitor banks (I-V characteristics of the CPA).

    Also, the dependence of the discharge current on the hydrogen gas pressure was obtained in this work(figure 3(a)).This graph shows that the current initially increases to its absolute value,and then it gradually decreases with increasing the pressure of hydrogen in the vacuum chamber.We explain this result by increasing the inductance and inductive resistance of the CPA circuit, which occurs while forming longlived plasma filaments.The absolute value of the discharge current was measured at pressure of 1.2 Torr.A relevant dependence was also presented in studies [22, 23].The dependence of the plasma filament’s self-magnetic field on the hydrogen pressure was shown in figure 3(b).As can be seen in figure 3(b), the magnitude of the magnetic field decreased more than twofold with tenfold increasing gas pressure.At higher gas pressures of 2-4 Torr, the magnetic field magnitude almost did not change.

    Figure 3.The dependence of the discharge current (a) and the self-generated magnetic field of plasma flow (b) on the gas pressure.

    Thus,the experimentally confirmed optimal gas pressure in the CPA is about 200 mTorr, because at this pressure it allowed us to achieve the maximum value of the plasma filament’s self-magnetic field.In this regard, the experiments were performed at pressures at which the maximum plasma compression (Z-pinch) was achieved.Consequently, these provide us for steady focusing of the plasma flux to the axis of the vacuum chamber.

    It can be seen from the obtained results shown in figure 3,that the absolute values of the magnetic field and the discharge current are observed at various gas pressures in the vacuum chamber.Therefore, a small change in the discharge current leads to a significant change in the plasma filament’s self-magnetic field.

    The typical images of the formation of plasma flow in hydrogen gas, obtained in one pulse shot at voltages of 4 kV and gas pressure of 2 Torr, through the quartz viewing side window of the vacuum chamber, are shown in figure 4(b).It should also be noted that the time evolution of the images matches with the time dependence of the magnetic probe signal shown in figure 4(a) (the CMOS camera was synchronized with the discharge in the experiment).The typical magnetic probe signals as in figure 4(a) were presented in [24].

    Figure 4.(a) The typical oscillograms of the plasma filament’s selfmagnetic field and (b) the instaneous images of the formation of pulsed plasma flow (side view).

    Figure 4(b) shows the plasma flux formation and acceleration frames recorded by using a high-speed camera.The camera frequency, resolution and exposure time are 280 000 frames s-1,256×64,3.1 μs.These frames were registered in a single pulse, but at different time moments.In figure 4(b),plasma acceleration direction is from the left to the right.The time moment of the vacuum arc discharge ignition between the coaxial cylindrical electrodes was fixed as the initial reference point.A high-current arc discharge produces strong self-induced magnetic field.Thus, plasma is accelerated toward the cathode tip by electromagnetic forces.As can be seen in figure 4(b)I, the head of the accelerated plasma flux was registered at a time of 17.8 μs.At the same time, as can be seen from the oscillograms, the magnetic probe No.1,placed between inner and outer electrodes registered the magnitude of the magnetic fields with an absolute value of the order of ~0.3 T, while the magnetic probe No.2 did not register this field.The magnetic probe No.2 was shielded by a ferromagnetic material with a 1 mm side slit.Respectively,it was sensitive only to the magnetic field created by the plasma flux current, flowing via the area where the magnetic probe No.2 is located.Then, in the subsequent frame, at t = 21.4 μs, the visible region of the high-speed camera is completely filled with bright plasma, and plasma flow with a duration of more than ~70 μs is formed (figure 4(b)II).At this time moment, plasma flow passes the vacuum chamber cross section, where the magnetic probe No.2 is mounted,and the voltage is gradually induced at its winding ends.Thus,the magnetic field first increases to a maximum value of~0.6 T, and then decreases.The magnetic field fluctuations are observed in the CPA at relatively high gas pressures, but they have a tendency to decay as the pressure increases.In this case,it can be assumed that instability occurs in the CPA due to the anomalous plasma resistance,and subfilament-type structures are formed.Therefore,the voltage decreasing from its absolute value on the probe No.2 and its fluctuation are explained with the plasma flux compression and disruption with its self-magnetic field.Thus, in one pulse, the plasma filament is disrupted into multiple sub-filaments.These indicate the creation of an inhomogeneous plasma fluxes.The estimated from two consecutive frames first sub-filament duration and directional velocity are ~70 μs, ~25.5 km s-1(figure 4(b)I).The following sub-filament was recorded after 96.4 μs (figure 4(b)III).The frame shown in figure 4(b)IV(140-160 μs) corresponds to the third peak of the magnetic probe signal No.2 with amplitude of 0.1 T (figure 4(a)).

    Figures 5(a) and (b) show the measured values of the magnetic field at a distance of 23.5 cm from the tip of the cathode and in space between the electrodes as a function of the charging voltage supplied to the capacitor banks.One can see that the value of the magnetic field in both cases increases proportionally with the voltage.Furthermore, the values of the magnetic field are on average 30% higher when the pure hydrogen is used as the plasma-forming gas than in case of the residual air.The physical interpretation of this phenomenon has not yet been revealed.

    Figure 5.The dependence of the plasma filament’s self-magnetic field on the voltage supplied to the capacitor banks: (a) at a distance of 23.5 cm from the cathode tip and (b) in the space between the electrodes.

    Figure 6.The oscillograms of magnetic probes obtained at gas pressures: (a) p = 60 mTorr, (b) p = 11 Torr.

    As can be seen from the oscillogram in figure 6(b), the increase in pressure up to 11 Torr leads to a decrease in amplitude of the No.2 magnetic probes signal.Note that probe No.2 was positioned at a 23.5 cm distance from the cathode tip, due to the fact that at this position the vacuum chamber-viewing window is located.Therefore, this distance is the most convenient and closest for observation by the high-speed camera.Moreover, the signal from the magnetic probe in interelectrode space repeats the current fluctuation in the accelerator circuit in comparison with figure 5(a)obtained at p = 60 mTorr.This indicates that,when the gas pressure is high, the resistance of the medium becomes larger, and as a result, the flow rate decreases, the primary plasma sub-filament is concentrated in interelectrode space.We also assume that this is due to an increase in mass of the plasma bullet,which can also affect the acceleration of plasma flow.One can see from figure 6(b), at higher pressures in the vacuum chamber, both positive and negative peaks appear in the oscillogram of the magnetic field.

    Figure 7.The ion and electron current waveforms.

    This fact shows the separation of plasma charges.The charge separation in plasma flow during the accelerating along the axis of the vacuum chamber can be observed in the FC signal also.When the FC collector voltage is zero, two peaks of opposite polarity are observed on the oscillogram(see figure 7).This happens when the collector of a FC is bombarded first by the electron and then by the ion beams,due to the center of a mass shift in space.We assume that this is firstly related to the principle of plasma flow acceleration in this type of accelerators and, secondly, to the collision processes between the particles: generally between directed ions and thermal neutrals.When the pressure of the gas in the vacuum chamber increases, the number of neutral particles increases too,resulting in an increased frequency of collisions of plasma charges with neutral particles.So, when neutrals collide frequently with ions, the ions begin to be slowed down.During this,the electrons have enough time to separate from the ions at a certain distance along the axis of acceleration.

    The plasma accelerators are known as potential sources of a wide variety of radiation such as electrons,ions,neutrons and x-rays,which has great interest for different applications.It is very important to investigate the ion beams and its general characteristics, such as energy, density, and uniformity for material processing, surface modification, ion implantation etc.In this regard,in this work we fabricated and employed FC for measuring the density of ions in plasma flux.In the experiment, FC was mounted in the path of the plasma flux.The resulting current of accelerated plasma particles absorbed by the collector was measured using an oscilloscope.The typical oscillogram at zero voltage on the collector, as well as at a potential of ±40 V is shown in figure 7 as an example.

    At zero voltage on the oscillogram one can observe peaks of negative and positive polarity, which is associated with the flux of electrons and ions with a time difference(70-75 μs),and indicates the separation of charge fluxes in plasma.After furthermore analysis of this result, we concluded, that for such a long period of time charge separation is not probable.Because time difference estimation is approximate and is obtained only on the waveform at zero voltage.It should be noted,that charge separation occurs, but only in a short time interval.Therefore,the ion beam current can be appeared earlier than the estimated time.We cannot see this from the FC oscillograms,obtained at zero collector potential,because the ion current can be very low and it can also be dampened by the electron current.The charge fluxes separation in plasma flow occurs in 5 μs, which is evaluated by applying a certain value of negative and positive potential to the collector.As the collector potential increases,for example, positive potential (+40 V), electrons will be attracted and ions will be rejected.In that case, the amplitude of the negative polarity current will increase on the oscillogram.When the potential is negative(-40 V),the amplitude of the current of positive polarity, i.e.ions, increases.The electron and ion currents do not flow to the collector simultaneously at the same collector voltages, as seen in figure 7.The time difference between two current signals is 5 μs.This can be explained by the fact that the collector adsorbs electrons and ions which are separated in space.

    Thus, the dependence of electron beam current on the collector voltage was obtained.The obtained I-V characteristic will be the same as that of the single planar probe,therefore for calculations of the electron temperature,we will use the expression (1), derived from the theory and assumption of the single probe that the particle energy distribution must be Maxwellian

    where e-lectrical charge, k-Boltzmann’s constant,Ie-the current of electrons,V-probe potential.

    The thermal velocity of ions is required to know to calculate their density.This can be estimated using Bohm criterion as follows:

    whereMi-ion mass,Te-electron temperature.Then, the density of ions is calculated by the following expression:

    whereIi-ion saturation current,S-area of the orifice,qelectrical charge.

    In our case, we assume that the plasma particle velocity distribution is Maxwellian, since the directional velocity of plasma flow calculated by the high-speed camera is 2.5×104m s-1,which is much less than the thermal velocity of electrons.The temperature of the electrons in this work was estimated from the slope of the Faraday cup I-V characteristic plotted on a logarithmic scale and obtained by averaging the currents measured in five repeated pulses at the same collector voltage.The averaging was performed,due to the operation of the experimental setup in the pulse mode.Thus,the calculated electron temperatures for the two types of the plasma-forming gases:air and hydrogen were 40 eV.However,as can be seen in figures 8(a) and (b), in our case, the I-V characteristic of the CPA deviates strongly from the exponential one.This can be explained by the fact that experimentally it is very difficult to obtain a real I-V characteristic, since plasma flows, generated in the plasma accelerator, are inhomogeneous.In addition, it is necessary to take into account the pulsed mode of operation of the accelerator,which limits the measurements in a large range and quantity.Usually, measurement errors with electric probes do not exceed 25% [25].In our case, the estimate of the error contribution to this is no more than 10%,which is due to the operation in the pulsed mode and the small number of repetitive measurements.Consequently,considering this contribution, we can judge that the total measurement error is 35%.

    The calculations of the electron temperature and ion density for two types of plasma-forming gases are shown in table 1.These results were compared with other works[15, 26] and are in good agreement.

    From the results of calculations,the following conclusion can be done.The temperature of electrons in plasma flow for hydrogen and residual air is practically the same.The density of ions in plasma flow,generated in residual gas is higher than that in hydrogen, which is possibly associated with multiple ionization at a higher gas breakdown voltage.

    Plasma flows, generated in high-current pulsed systems,are themselves complex objects of investigation due to the non-homogeneity of the fields and currents created by accelerating plasma.Therefore, the complex plasma diagnostics allowed us to establish representation of the formation of plasma flows in the CPA channel.In particular, it allowed us to observe such phenomena as the plasma filament disruption, the charge (electrons and ions) separation in plasma flow, even with a limited number of measurement points.These phenomena were investigated depending on the external conditions of experiments also.As a result, the optimal parameters of the experimental setup were determined: working gas pressure of 200 mTorr, capacitor bank charge voltage of 5 kV.

    Figure 8.I-V characteristics of electron currents on a logarithmic scale for two different gases: (a) hydrogen and (b) air.

    Table 1.The calculation results.

    5.Conclusion

    As a result of the complex diagnostics, the parameters of the external electrical circuit of the CPA and the energetic,dynamic characteristics of the plasma filament were obtained.From the results obtained, it was revealed that the current curves represent rapidly damped aperiodic oscillations with a small number of half-periods.Such a condition is advantageous and indicates the low value of inductive resistance of the circuit, including connecting wires, capacitors and vacuum arrester.One can conclude from this that the stored energy is effectively transported into the discharge.Also, it was found that the current depends linearly on the voltage supplied to the capacitor banks, which is explained by the high conductivity and the low resistance of the plasma filament.In the experiment,the absolute current values of 80 and 65 kA for hydrogen and air plasma-forming gases were obtained at a maximum voltage of 5 kV and working gas pressure of 60 mTorr.With increasing gas pressure, it was found that the current amplitudes decrease from the maximum value, which is explained by an increase in inductance associated with the formation of long-lived plasma filaments.Based on the results obtained by the high-speed CMOS camera, the plasma filament disruption was also observed in the experiments, which explains a quasiperiodic ejection of plasma from space between the cathode and anode.The peaks observed in the oscillograms of the magnetic probe also indicate this phenomenon.The electron temperature and ion density in plasma flow were measured by using FC.It is shown that in the dynamic plasma, even when the collector potential is zero,currents are generated at the collector,which is explained by bombarding of the collector surface by charged particles.Therefore, in the calculations, the parasitic currents are extracted from the major currents measured at various collector bias voltages.As a result,we eliminated the influence of background currents on the error of the measurement.At the conclusion of this work, we achieved the following purposes: a physical representation of plasma flow formation in the CPA channel as a function of working gas pressure and the capacitor bank charging voltage has been established; the optimal experimental conditions were determined at which the maximum plasma filament compression can be obtained; charge fluxes separation in space relative to the common center of mass with a time difference of ~5 μs was observed by using FC; the oscillograms of magnetic probes signals and plasma flow acceleration frames simultaneously confirm plasma filament disruption.Especially it is noticeable at high pressures,which is related with an increase in the medium resistance and collision frequency of ions with neutrals.In the future, we will use these results, the new insights we have gained, and proven diagnostic methods in our fusion research work.In particular,to further modify and investigate the processes of plasma flow interaction with materials that are candidates for the fusion reactor’s first-wall,and to improve the experimental setup for producing homogeneous plasma flows with the best energetic parameters.

    Acknowledgments

    This work was supported by the Ministry of Education and Science of the Republic of Kazakhstan (IRN AP08053373).

    超碰97精品在线观看| 免费在线观看日本一区| 久久中文字幕一级| 麻豆成人av在线观看| 韩国av一区二区三区四区| 午夜免费鲁丝| 精品第一国产精品| 亚洲成国产人片在线观看| 少妇的丰满在线观看| 免费不卡黄色视频| 操美女的视频在线观看| 熟女少妇亚洲综合色aaa.| 97超级碰碰碰精品色视频在线观看| www.www免费av| 成年人免费黄色播放视频| 亚洲欧美日韩另类电影网站| 人妻丰满熟妇av一区二区三区| 女人被狂操c到高潮| 午夜精品在线福利| 亚洲七黄色美女视频| 日韩精品免费视频一区二区三区| 精品高清国产在线一区| 日本撒尿小便嘘嘘汇集6| 90打野战视频偷拍视频| 色哟哟哟哟哟哟| 女人爽到高潮嗷嗷叫在线视频| 日本wwww免费看| 国产欧美日韩一区二区精品| 又大又爽又粗| 亚洲男人天堂网一区| 国产xxxxx性猛交| 黑人巨大精品欧美一区二区mp4| 亚洲一区中文字幕在线| 精品午夜福利视频在线观看一区| 中出人妻视频一区二区| 桃红色精品国产亚洲av| 18禁裸乳无遮挡免费网站照片 | 精品一品国产午夜福利视频| 在线观看免费视频网站a站| avwww免费| 欧美精品啪啪一区二区三区| 午夜成年电影在线免费观看| 国产片内射在线| 多毛熟女@视频| 免费一级毛片在线播放高清视频 | 精品久久久久久电影网| av天堂在线播放| 国产精品野战在线观看 | 久久 成人 亚洲| 久久伊人香网站| 精品一区二区三区四区五区乱码| 中亚洲国语对白在线视频| 免费观看人在逋| 一区福利在线观看| 少妇粗大呻吟视频| 日韩免费高清中文字幕av| 亚洲av成人一区二区三| 欧美最黄视频在线播放免费 | 男人的好看免费观看在线视频 | 嫩草影院精品99| 琪琪午夜伦伦电影理论片6080| 国产一区二区激情短视频| 日本黄色日本黄色录像| 中文字幕人妻丝袜一区二区| 每晚都被弄得嗷嗷叫到高潮| 久久 成人 亚洲| 一级毛片高清免费大全| 怎么达到女性高潮| 国产aⅴ精品一区二区三区波| 老司机午夜福利在线观看视频| 国产欧美日韩综合在线一区二区| 亚洲五月婷婷丁香| 一进一出好大好爽视频| 99re在线观看精品视频| 亚洲国产精品999在线| 国产真人三级小视频在线观看| 99国产精品99久久久久| 中出人妻视频一区二区| 熟女少妇亚洲综合色aaa.| 91国产中文字幕| 国产黄色免费在线视频| 夜夜爽天天搞| 老鸭窝网址在线观看| 日韩欧美三级三区| 欧美激情 高清一区二区三区| 国产主播在线观看一区二区| 日韩成人在线观看一区二区三区| 18禁观看日本| 大型黄色视频在线免费观看| 成人特级黄色片久久久久久久| 国产精品亚洲一级av第二区| 亚洲欧美日韩另类电影网站| 久热爱精品视频在线9| www.精华液| 两性夫妻黄色片| 国产精品一区二区三区四区久久 | 超碰97精品在线观看| 99久久人妻综合| 涩涩av久久男人的天堂| 国产又色又爽无遮挡免费看| 久久天堂一区二区三区四区| 在线十欧美十亚洲十日本专区| 亚洲午夜理论影院| 色综合婷婷激情| 9色porny在线观看| 真人做人爱边吃奶动态| 亚洲成国产人片在线观看| 免费搜索国产男女视频| 最近最新免费中文字幕在线| www.www免费av| 久久精品亚洲av国产电影网| 欧美激情极品国产一区二区三区| 在线观看一区二区三区激情| www.自偷自拍.com| av免费在线观看网站| 国产精品亚洲一级av第二区| 中文亚洲av片在线观看爽| 日韩欧美国产一区二区入口| 精品国产一区二区三区四区第35| 亚洲三区欧美一区| 成人影院久久| 国产99久久九九免费精品| 国产区一区二久久| www国产在线视频色| 婷婷丁香在线五月| 丝袜人妻中文字幕| 成人免费观看视频高清| 不卡av一区二区三区| 久久久国产一区二区| 正在播放国产对白刺激| 人人妻人人添人人爽欧美一区卜| 免费在线观看日本一区| 免费不卡黄色视频| 男人的好看免费观看在线视频 | 欧美亚洲日本最大视频资源| 一级黄色大片毛片| 国产成人免费无遮挡视频| 国产精品野战在线观看 | 日韩免费av在线播放| 91字幕亚洲| 日韩中文字幕欧美一区二区| 亚洲欧美精品综合一区二区三区| xxx96com| 黄色丝袜av网址大全| 亚洲中文日韩欧美视频| 一级毛片女人18水好多| 精品久久久精品久久久| 成年版毛片免费区| 999久久久国产精品视频| 亚洲一区二区三区不卡视频| 亚洲 欧美一区二区三区| 别揉我奶头~嗯~啊~动态视频| 久久久久久亚洲精品国产蜜桃av| 操出白浆在线播放| 一个人观看的视频www高清免费观看 | 国产精品99久久99久久久不卡| 一级毛片精品| 国产精品久久久久成人av| 亚洲成人免费电影在线观看| 亚洲精品中文字幕在线视频| 99国产极品粉嫩在线观看| 美女午夜性视频免费| 天天躁夜夜躁狠狠躁躁| 身体一侧抽搐| 亚洲av成人不卡在线观看播放网| 亚洲国产毛片av蜜桃av| 日韩精品免费视频一区二区三区| 亚洲 国产 在线| 国产成人欧美| 黄色 视频免费看| 国产精品一区二区精品视频观看| 久久精品亚洲精品国产色婷小说| 亚洲精品成人av观看孕妇| 一进一出抽搐gif免费好疼 | 18禁裸乳无遮挡免费网站照片 | 精品免费久久久久久久清纯| 亚洲av第一区精品v没综合| 大型av网站在线播放| 日韩欧美免费精品| av超薄肉色丝袜交足视频| 成人黄色视频免费在线看| 少妇的丰满在线观看| 亚洲精华国产精华精| 久久香蕉激情| www国产在线视频色| 国产aⅴ精品一区二区三区波| 夜夜躁狠狠躁天天躁| 国产熟女xx| 一级黄色大片毛片| 国产伦一二天堂av在线观看| 久久国产精品影院| 99久久99久久久精品蜜桃| xxxhd国产人妻xxx| 午夜日韩欧美国产| www.熟女人妻精品国产| 搡老乐熟女国产| 搡老熟女国产l中国老女人| 一级黄色大片毛片| 亚洲va日本ⅴa欧美va伊人久久| 国产99白浆流出| 欧美av亚洲av综合av国产av| 日韩欧美国产一区二区入口| 99热只有精品国产| 精品免费久久久久久久清纯| 狠狠狠狠99中文字幕| 亚洲黑人精品在线| 成人亚洲精品一区在线观看| 日韩人妻精品一区2区三区| 日韩欧美免费精品| 午夜福利在线观看吧| 女警被强在线播放| 人人妻人人澡人人看| 国产精品 欧美亚洲| 日本三级黄在线观看| 亚洲 欧美一区二区三区| 91老司机精品| 怎么达到女性高潮| 亚洲国产精品合色在线| 亚洲国产欧美日韩在线播放| av免费在线观看网站| 精品国产亚洲在线| 婷婷精品国产亚洲av在线| 精品乱码久久久久久99久播| 90打野战视频偷拍视频| 他把我摸到了高潮在线观看| 亚洲av电影在线进入| 国产一区二区在线av高清观看| 精品国产乱码久久久久久男人| 叶爱在线成人免费视频播放| 精品国产乱码久久久久久男人| 久久久久久久午夜电影 | 色尼玛亚洲综合影院| 9色porny在线观看| 国产精品久久久久成人av| 一边摸一边做爽爽视频免费| 日韩成人在线观看一区二区三区| 19禁男女啪啪无遮挡网站| 国产成年人精品一区二区 | 久9热在线精品视频| 电影成人av| 淫秽高清视频在线观看| av超薄肉色丝袜交足视频| 中文字幕色久视频| 99久久综合精品五月天人人| 中文字幕av电影在线播放| 国产精品免费视频内射| 精品福利观看| 精品久久久精品久久久| 12—13女人毛片做爰片一| 免费少妇av软件| 亚洲av成人av| 久久亚洲真实| 久久精品国产99精品国产亚洲性色 | 亚洲av日韩精品久久久久久密| 18禁国产床啪视频网站| 亚洲精品成人av观看孕妇| 涩涩av久久男人的天堂| 国产精品久久久人人做人人爽| 午夜a级毛片| 久久中文字幕人妻熟女| 欧美精品一区二区免费开放| 成熟少妇高潮喷水视频| 看免费av毛片| 国产av一区在线观看免费| 免费日韩欧美在线观看| 一本大道久久a久久精品| 超碰97精品在线观看| 国产又色又爽无遮挡免费看| 成在线人永久免费视频| 精品国产美女av久久久久小说| 少妇裸体淫交视频免费看高清 | 最好的美女福利视频网| 日韩欧美三级三区| 欧美成人性av电影在线观看| 99香蕉大伊视频| 精品国产一区二区久久| www.精华液| 午夜久久久在线观看| 91九色精品人成在线观看| 女人爽到高潮嗷嗷叫在线视频| 精品一区二区三区av网在线观看| 黄色成人免费大全| 精品一区二区三卡| 嫩草影院精品99| 亚洲一卡2卡3卡4卡5卡精品中文| 在线观看免费日韩欧美大片| 免费在线观看亚洲国产| 亚洲欧洲精品一区二区精品久久久| 丰满人妻熟妇乱又伦精品不卡| 久久久久精品国产欧美久久久| 国产成+人综合+亚洲专区| 热re99久久精品国产66热6| 两个人免费观看高清视频| 欧美日韩乱码在线| 久久香蕉激情| 亚洲情色 制服丝袜| 亚洲成人免费电影在线观看| 婷婷六月久久综合丁香| 成人免费观看视频高清| 久久性视频一级片| 亚洲精品美女久久av网站| 高清毛片免费观看视频网站 | 成人三级黄色视频| 中文字幕另类日韩欧美亚洲嫩草| 亚洲精品国产精品久久久不卡| 又大又爽又粗| 黄色女人牲交| 成年人免费黄色播放视频| 黄片大片在线免费观看| 亚洲欧洲精品一区二区精品久久久| 一个人观看的视频www高清免费观看 | 性少妇av在线| 国产精品日韩av在线免费观看 | 新久久久久国产一级毛片| 宅男免费午夜| 色哟哟哟哟哟哟| 成人三级黄色视频| 国产精品免费一区二区三区在线| 亚洲九九香蕉| 久久99一区二区三区| 亚洲成国产人片在线观看| 欧美国产精品va在线观看不卡| 在线观看一区二区三区| 久久性视频一级片| 亚洲精品av麻豆狂野| 国产在线观看jvid| 国产伦人伦偷精品视频| 精品久久蜜臀av无| 黄色成人免费大全| 成人亚洲精品一区在线观看| 欧美成人免费av一区二区三区| 国产99白浆流出| 久久中文字幕人妻熟女| 亚洲第一欧美日韩一区二区三区| 亚洲精品成人av观看孕妇| 波多野结衣高清无吗| 在线天堂中文资源库| 嫩草影院精品99| 日本wwww免费看| 欧美日韩国产mv在线观看视频| 国产精品综合久久久久久久免费 | 国产精品 欧美亚洲| xxxhd国产人妻xxx| 好男人电影高清在线观看| 亚洲自拍偷在线| 久久精品国产亚洲av高清一级| av超薄肉色丝袜交足视频| 午夜久久久在线观看| 亚洲国产精品合色在线| 日本 av在线| 人成视频在线观看免费观看| 久久久久国产一级毛片高清牌| 成年女人毛片免费观看观看9| 人妻丰满熟妇av一区二区三区| 亚洲人成网站在线播放欧美日韩| 丰满的人妻完整版| 99国产精品一区二区蜜桃av| 久久精品国产99精品国产亚洲性色 | 美国免费a级毛片| 成年版毛片免费区| 黄色女人牲交| av网站免费在线观看视频| 中出人妻视频一区二区| 老司机深夜福利视频在线观看| 国产av在哪里看| 欧美激情高清一区二区三区| 久久精品国产亚洲av高清一级| 在线观看午夜福利视频| av天堂久久9| 麻豆久久精品国产亚洲av | 成人免费观看视频高清| 午夜老司机福利片| 日韩一卡2卡3卡4卡2021年| 亚洲avbb在线观看| 色老头精品视频在线观看| 一边摸一边做爽爽视频免费| 亚洲激情在线av| 高清毛片免费观看视频网站 | e午夜精品久久久久久久| 狂野欧美激情性xxxx| 欧美大码av| netflix在线观看网站| 国产亚洲精品综合一区在线观看 | 日本免费一区二区三区高清不卡 | 又黄又爽又免费观看的视频| 新久久久久国产一级毛片| 国产一区二区三区在线臀色熟女 | 91老司机精品| 国产99久久九九免费精品| 一边摸一边抽搐一进一出视频| 国产无遮挡羞羞视频在线观看| 亚洲熟妇熟女久久| 在线观看免费午夜福利视频| 97超级碰碰碰精品色视频在线观看| 亚洲av片天天在线观看| 久久青草综合色| 日韩中文字幕欧美一区二区| 在线视频色国产色| 亚洲欧美一区二区三区黑人| 久久性视频一级片| 夜夜夜夜夜久久久久| 亚洲中文日韩欧美视频| 咕卡用的链子| 亚洲一区二区三区色噜噜 | 丝袜在线中文字幕| 国产真人三级小视频在线观看| 欧美日韩中文字幕国产精品一区二区三区 | 色综合婷婷激情| 精品电影一区二区在线| 欧美日韩中文字幕国产精品一区二区三区 | 日本 av在线| 亚洲色图 男人天堂 中文字幕| 男人的好看免费观看在线视频 | 天堂中文最新版在线下载| 成人三级做爰电影| 最新美女视频免费是黄的| 亚洲中文字幕日韩| 新久久久久国产一级毛片| 满18在线观看网站| 中出人妻视频一区二区| 变态另类成人亚洲欧美熟女 | 国产欧美日韩一区二区三区在线| 18禁观看日本| 久久中文字幕一级| 亚洲七黄色美女视频| 国产极品粉嫩免费观看在线| 欧美一级毛片孕妇| 激情在线观看视频在线高清| 一区福利在线观看| 天天添夜夜摸| 国产熟女午夜一区二区三区| 精品久久蜜臀av无| a级毛片黄视频| 亚洲精品一二三| 精品一区二区三区四区五区乱码| 在线观看www视频免费| 757午夜福利合集在线观看| 色在线成人网| 91国产中文字幕| 精品福利观看| 国产一区在线观看成人免费| 三级毛片av免费| 热re99久久精品国产66热6| 电影成人av| 热re99久久国产66热| 长腿黑丝高跟| 国产高清videossex| 人人妻,人人澡人人爽秒播| 欧美日韩一级在线毛片| 国产精品永久免费网站| 国产成人欧美| 91av网站免费观看| 黄色毛片三级朝国网站| 久久中文看片网| 亚洲精品国产一区二区精华液| 欧美+亚洲+日韩+国产| 亚洲一区二区三区欧美精品| 欧美国产精品va在线观看不卡| 色综合站精品国产| 亚洲欧洲精品一区二区精品久久久| 精品人妻1区二区| 如日韩欧美国产精品一区二区三区| 亚洲国产精品合色在线| 欧美日韩亚洲综合一区二区三区_| 国产aⅴ精品一区二区三区波| 久久国产精品男人的天堂亚洲| 久久久精品国产亚洲av高清涩受| 变态另类成人亚洲欧美熟女 | 满18在线观看网站| netflix在线观看网站| 亚洲第一青青草原| 中出人妻视频一区二区| 国产三级在线视频| www国产在线视频色| 日韩精品青青久久久久久| 欧美人与性动交α欧美精品济南到| 热99国产精品久久久久久7| 成在线人永久免费视频| 人妻丰满熟妇av一区二区三区| av视频免费观看在线观看| 自拍欧美九色日韩亚洲蝌蚪91| 婷婷精品国产亚洲av在线| 一进一出抽搐动态| 女人被躁到高潮嗷嗷叫费观| 亚洲人成电影观看| 欧美日韩av久久| 一个人观看的视频www高清免费观看 | 亚洲精品一区av在线观看| 国产精品美女特级片免费视频播放器 | 黄片大片在线免费观看| 成年版毛片免费区| cao死你这个sao货| 欧美大码av| 亚洲精品在线观看二区| 久久这里只有精品19| 在线观看www视频免费| 久久亚洲真实| 亚洲av五月六月丁香网| 午夜亚洲福利在线播放| 亚洲午夜理论影院| 亚洲精品美女久久av网站| www.精华液| 天天躁狠狠躁夜夜躁狠狠躁| 久久香蕉激情| 亚洲一区二区三区色噜噜 | 美女扒开内裤让男人捅视频| 久久精品91蜜桃| 日韩人妻精品一区2区三区| 日日夜夜操网爽| 午夜福利一区二区在线看| 精品午夜福利视频在线观看一区| 神马国产精品三级电影在线观看 | √禁漫天堂资源中文www| 无遮挡黄片免费观看| 99热国产这里只有精品6| 12—13女人毛片做爰片一| 亚洲五月婷婷丁香| 欧美精品亚洲一区二区| 久久99一区二区三区| 久久精品亚洲精品国产色婷小说| 国产亚洲精品久久久久5区| 国产精品野战在线观看 | 久久人妻av系列| 午夜影院日韩av| 国产亚洲欧美98| 精品人妻1区二区| 在线观看免费日韩欧美大片| 亚洲专区国产一区二区| 18禁美女被吸乳视频| 精品久久久精品久久久| 99精国产麻豆久久婷婷| 欧美日韩国产mv在线观看视频| 国产一区二区三区在线臀色熟女 | 欧美 亚洲 国产 日韩一| 成在线人永久免费视频| av国产精品久久久久影院| 91精品三级在线观看| 日韩大尺度精品在线看网址 | 国产精品自产拍在线观看55亚洲| 99久久99久久久精品蜜桃| 琪琪午夜伦伦电影理论片6080| 在线观看免费高清a一片| 亚洲精品粉嫩美女一区| 中文字幕av电影在线播放| 亚洲精品粉嫩美女一区| 啦啦啦在线免费观看视频4| 两个人免费观看高清视频| 熟女少妇亚洲综合色aaa.| 亚洲成人精品中文字幕电影 | 免费在线观看视频国产中文字幕亚洲| 久久中文字幕人妻熟女| 久久精品亚洲精品国产色婷小说| 18禁美女被吸乳视频| 亚洲午夜理论影院| 亚洲九九香蕉| 久久国产亚洲av麻豆专区| 免费高清视频大片| 国产99久久九九免费精品| a级片在线免费高清观看视频| 成人18禁高潮啪啪吃奶动态图| av欧美777| bbb黄色大片| 动漫黄色视频在线观看| 在线观看免费日韩欧美大片| 黄色成人免费大全| 亚洲国产欧美一区二区综合| 中文字幕高清在线视频| 久久久久久大精品| 亚洲国产欧美网| 露出奶头的视频| 搡老熟女国产l中国老女人| 国内毛片毛片毛片毛片毛片| 不卡av一区二区三区| 十八禁网站免费在线| 国产在线精品亚洲第一网站| 国产精品一区二区三区四区久久 | 亚洲一区二区三区色噜噜 | 亚洲一区高清亚洲精品| 一级a爱片免费观看的视频| 老司机福利观看| 国产一卡二卡三卡精品| www国产在线视频色| 亚洲片人在线观看| 丰满人妻熟妇乱又伦精品不卡| 国产av一区二区精品久久| xxx96com| 十八禁网站免费在线| 精品卡一卡二卡四卡免费| 亚洲av第一区精品v没综合| 亚洲aⅴ乱码一区二区在线播放 | 99国产精品一区二区三区| 精品久久久久久久毛片微露脸| 国产免费男女视频| 免费看十八禁软件| 麻豆成人av在线观看| 久久午夜亚洲精品久久| 国产无遮挡羞羞视频在线观看| 欧美日韩一级在线毛片| 日韩免费av在线播放| 女人精品久久久久毛片| 国产又爽黄色视频| 亚洲精品国产一区二区精华液| 婷婷六月久久综合丁香| 精品久久久久久久久久免费视频 | 国产欧美日韩综合在线一区二区| 免费看a级黄色片| 人妻久久中文字幕网| 精品国产美女av久久久久小说| 怎么达到女性高潮| 午夜福利在线免费观看网站| 在线观看舔阴道视频| 国产成人精品无人区| 夜夜看夜夜爽夜夜摸 | 91在线观看av| 免费人成视频x8x8入口观看|