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

    Line identification of boron and nitrogen emissions in extreme- and vacuumultraviolet wavelength ranges in the impurity powder dropping experiments of the Large Helical Device and its application to spectroscopic diagnostics

    2021-08-05 08:29:06TetsutarouOISHI大石鉄太郎NaokoASHIKAWA蘆川直子FedericoNESPOLISuguruMASUZAKI増崎貴MamoruSHOJI莊司主EricGILSONRobertLUNSFORDShigeruMORITA森田繁MotoshiGOTO後藤基志YasukoKAWAMOTO川本靖子ChihiroSUZUKI鈴木千尋ZhenSUN孫震AlexNAGYDavidGATESandTomohiroMORIS
    Plasma Science and Technology 2021年8期
    關(guān)鍵詞:森田直子

    Tetsutarou OISHI (大石鉄太郎), Naoko ASHIKAWA (蘆川直子),Federico NESPOLI, Suguru MASUZAKI (増崎貴), Mamoru SHOJI (莊司主),Eric P GILSON, Robert LUNSFORD, Shigeru MORITA (森田繁),Motoshi GOTO (後藤基志), Yasuko KAWAMOTO (川本靖子),Chihiro SUZUKI(鈴木千尋),Zhen SUN(孫震),Alex NAGY,David A GATES and Tomohiro MORISAKI (森崎友宏)

    1 National Institute for Fusion Science,National Institutes of Natural Sciences,322-6 Oroshi-cho,Toki 509-5292, Gifu, Japan

    2 Department of Fusion Science, The Graduate University for Advanced Studies, SOKENDAI, 322-6 Oroshi-cho, Toki 509-5292, Gifu, Japan

    3 Princeton Plasma Physics Laboratory, 100 Stellarator Road, Princeton, NJ 08540, United States of America

    Abstract An impurity powder dropper was installed in the 21st campaign of the Large Helical Device experiment (Oct.2019–Feb.2020) under a collaboration between the National Institute for Fusion Science and the Princeton Plasma Physics Laboratory for the purposes of real-time wall conditioning and edge plasma control.In order to assess the effective injection of the impurity powders,spectroscopic diagnostics were applied to observe line emission from the injected impurity.Thus,extreme-ultraviolet (EUV) and vacuum-ultraviolet (VUV) emission spectra were analyzed to summarize observable impurity lines with B and BN powder injection.Emission lines released from B and N ions were identified in the EUV wavelength range of 5–300 ? measured using two grazing incidence flat-field EUV spectrometers and in the VUV wavelength range of 300–2400 ? measured using three normal incidence 20 cm VUV spectrometers.BI–BV and NIII–NVII emission lines were identified in the discharges with the B and BN powder injection, respectively.Useful B and N emission lines which have large intensities and are isolated from other lines were successfully identified as follows: BI (1825.89, 1826.40) ? (blended), BII 1362.46 ?, BIII (677.00, 677.14,677.16) ? (blended), BIV 60.31 ?, BV 48.59 ?, NIII (989.79, 991.51, 991.58) ? (blended), NIV 765.15 ?,NV(209.27,209.31)?(blended),NVI 1896.80 ?,and NVII 24.78 ?.Applications of the line identifications to the advanced spectroscopic diagnostics were demonstrated,such as the vertical profile measurements for the BV and NVII lines using a space-resolved EUV spectrometer and the ion temperature measurement for the BII line using a normal incidence 3 m VUV spectrometer.

    Keywords:plasma spectroscopy,extreme ultraviolet, vacuum ultraviolet, magnetically confined fusion, impurity seeding, wall conditioning

    1.Introduction

    Impurity injection experiments using submillimeter-sized powders which consist of a simple substance or chemical compound of low-Zelements,such as Li,B,or N,have been conducted in several magnetically confined torus plasma devices for the purposes of real-time wall conditioning or edge plasma control.According to the results reported in tokamak devices,recycling of neutral H particles was suppressed by B or BN powder injection, which indicated that the injected B impurities contributed a real-time boronization (B wall coating) [1–3].Stabilization of the edge localized mode through modification of edge plasma profiles by injection of Li or B powders has also been reported [4, 5].Establishment of this kind of real-time plasma control method becomes more important in plasma confinement experiments with longer discharge pulses that progress toward the realization of thermonuclear fusion reactors with steady-state operation.Thus,the stellarator devices can provide an attractive research platform of investigation on behaviors and roles of injected impurities in steady-state plasma operations.In the Large Helical Device(LHD),an impurity powder dropper(IPD)was installed in the 21st campaign of the experiment (Oct.2019–Feb.2020) under a collaboration between the National Institute for Fusion Science and the Princeton Plasma Physics Laboratory [6].In order to assess the effective injection of the impurity powders, spectroscopic diagnostics were applied to observe line emission from the injected impurity.In the present study,extreme-ultraviolet(EUV)and vacuum-ultraviolet(VUV)emission spectra are analyzed to summarize observable impurity lines with B and BN powder injection.Based on the line identification, applications of the observed line to advanced spectroscopic diagnostics are demonstrated such as the spatial profile measurement of the emission intensities or the ion temperature measurement for the injected impurity ions.

    2.Spectroscopic diagnostics for the IPD experiment in the LHD

    The LHD is a superconducting plasma confinement device which employs a heliotron magnetic configuration.The plasma in the standard configuration has major and averaged minor radii of 3.6 and 0.64 m,respectively[7].The coil system consists of a set of two continuous superconducting helical coils with poloidal and toroidal pitch numbers of 2 and 10, respectively, and also three pairs of superconducting poloidal coils.The observation geometry of the spectroscopic diagnostics in the LHD, which includes two flat-field grazing incidence EUV spectrometers(denoted as‘EUV Short’[8]and‘EUV Long’[9]),three 20 cm normal incidence VUV spectrometers(denoted as‘VUV 109L’,‘VUV 106R’, and ‘VUV 102R’ [10]), a space-resolved EUV spectrometer(denoted as‘EUV Short2’[11]),and a 3 m normal incidence VUV spectrometer (denoted as ‘VUV 3 m’ [12]), is illustrated in figure 1.Figure 1(a)shows a top view of the optical axes of the spectrometers and vacuum magnetic surfaces with a magnetic axis positionRaxof 3.6 m.

    All of the spectroscopic diagnostics systems are located at an outer port of the LHD which is called‘10-O’.The IPD is installed at an upper port called‘2.5-U’,which is separated by 90° clockwise in the toroidal direction from the 10-O port.The detailed location of IPD installation was decided based on a simulation study on powder trajectory calculated with the coupled EMC3-EIRENE and DUSTT codes [13].The powders are dropped into the plasma by vibrating the powder reservoirs by piezo blades with an injection rate which can be scanned from a few 10 μg s?1to 900 mg s?1[14].The IPD is equipped with an accelerometer to measure the amplitude of the vibration of the powder reservoirs.A flowmeter consisting of a collimated beam of visible light crossing a guide tube through which the powders fall, and a light detector, is installed in the IPD to measure the injection rate of the powders.The injection rate is obtained by measuring the attenuation of the visible light in the flowmeter.The diameters of B and BN powder particles are roughly 150 μm and 60 μm,respectively.In the powder dropping experiments described in the next chapter, three kinds of spectroscopic measurements are conducted as follows.

    (1)EUV and VUV wavelength spectral measurements with high time resolution: the spectrometers of EUV Short,EUV Long, VUV 109L, VUV106R, and VUV102 cover the wavelength ranges of 5–60 ?, 100–300 ?,300–1050 ?,970–1870 ?,and 1510–2400 ?,respectively.The vertical observation range on a plasma cross section including optical axes of these spectrometers is illustrated in figure 1(b).CCD detectors (1024×256 pixels, pixel size 26×26 μm2, Andor DO420-BN) are placed at the positions of the exit slits of the spectrometers.A CCD data acquisition operational mode applied in this experiment is called‘full-binning’mode in which all CCD pixels aligned in the vertical direction are replaced by a single channel and the vertical spatial resolution is entirely eliminated.A time resolution for the spectral measurements is 5 ms with the full-binning data acquisition mode.

    (2)Emission intensity profile measurement with spaceresolved EUV spectroscopy: EUV Short2 has a spatial resolution in the vertical direction by mounting a spaceresolved slit at the position of the entrance slit.A CCD detector (1024×255 pixels, pixel size 26×26 μm2,Andor DO920P-BN) is placed at the positions of the exit slits of the spectrometer.The long axis of the detector plane is directed perpendicular to the midplane of the LHD plasma and the spatial observation range is set to project the emission profile of the upper half of the plasma at the horizontally elongated poloidal cross section on the detector.Five-pixel binning is applied to the long axis of the detector with 1024 pixels.Thus,the observation range is resolved by 204 viewing chords in the vertical direction.The observation range covers approximately 500 mm in the vertical direction with the spatial resolution of 2.6 mm around the plasma center as illustrated in figure 1(c).The short axis of the detector is directed to the wavelength dispersion direction and twopixel binning is applied.The temporal resolution is 100 ms under the operational conditions employed in this study.The wavelength range is set to cover 40–58 ? to observe vertical intensity profiles of BV 48.59 ? and NVII 24.78×2 ? simultaneously.

    Figure 1.Observation geometry of EUV and VUV spectroscopy diagnostics in the LHD.(a) Top view of optical axes of two grazing incidence flat-field EUV spectrometers (‘EUV Short’ and ‘EUV Long’), three normal incidence 20 cm VUV spectrometers (‘VUV 109L’,‘VUV 106R’, and ‘VUV 102R’), a space-resolved EUV spectrometer (‘EUV Short2’), and a normal incidence 3 m VUV spectrometer(‘VUV 3 m’) shown together with magnetic surfaces (Rax=3.6 m).(b) Plasma cross section including the vertical observation range of(b) EUV Short, EUV Long, VUV 109L, VUV 106R, and VUV 102R, and (c) EUV Short2 and VUV 3 m.

    (3)Impurity ion temperature measurement with high spectral resolution VUV spectroscopy: a 3 m normal incidence VUV spectrometer has been developed to measure the radial distribution of VUV lines in the wavelength range of 300–3200 ? in the edge plasmas of LHD.A CCD detector (1024×1024 pixels, pixel size 13×13 μm2,Andor DO934P-BN)is placed at the positions of the exit slits of the spectrometer.The high spectral resolution of the spectroscopic system with the wavelength dispersion of 0.037 ?/CCD pixel enables us to measure the Doppler profiles of impurity line spectra precisely.The vertical observation range covers the entire vertical height of the LHD plasmas at the horizontally elongated poloidal cross section as illustrated in figure 1(c).The temporal resolution is 20 ms with a full-binning CCD data acquisition mode.The line shape of the wavelength spectrum has a Gaussian profile if the ions are assumed to have a Maxwellian velocity distribution.The wavelength range is set to cover 1350–1387 ? to observe the spectral shape of BII 1362.46 ?.

    3.EUV and VUV spectra of boron and nitrogen emissions

    3.1.Spectra in the boron powder dropping experiments

    Figure 2.Typical waveform of the B powder dropping experiment in the LHD: (a) heating power of ECH and NBI, and gas-puffing duration of H2, (b) amplitude of accelerometer signal Vacc and attenuation of flowmeter signal ΔF/Fof the IPD,(c)central electron temperature Te0, (d) edge electron temperature Te(a99), (e) lineaveraged electron density ne, (f) edge electron density ne(a99),(g) total radiated power Prad, and (h) plasma stored energy Wp.

    Figure 2 shows a typical waveform of a B powder dropping experiment in LHD withRaxof 3.6 m and toroidal magnetic fieldBtof 2.75 T in the counter-clockwise direction from the top view.In this figure, the temporal evolution of (a) the heating power of electron cyclotron heating (ECH) and neutral beam injection (NBI), and gas-puffing duration of H2,(b) the amplitude of the accelerometer signalVaccand attenuation of flowmeter signal ΔF/Fof the IPD, (c) the central electron temperatureTe0, (d) the edge electron temperatureTe(a99), (e) the line-averaged electron densityne,(f) the edge electron densityne(a99), (g) the total radiated powerPrad, and (h) the plasma stored energyWpare plotted together.The plasma edgea99was defined as the effective minor radius in which 99% of electron stored energy was enclosed [15].The plasma was initiated by ECH and further heated by NBI.H2gas was puffed from 3.1 to 5.3 s to control the electron density.A voltage is applied to the IPD feeder piezoelectric blades from 3.0 to 3.5 s as indicated in the signal ofVaccin figure 2(b).As shown in the figures, signals of plasma response to the powder injection are observed at 4.0 s such as an increase inne,ne(a99),andWp,a slight decrease inTe(a99), and positive spikes inPrad.It has already been confirmed that a typical time delay is about 1 s from the onset of the voltage application to the IPD feeder piezoelectric blade to the response by the plasma [6].

    Figure 3 shows EUV/VUV spectra including B line emissions.The spectral data were averaged over 3.5–3.6 s(blue dashed lines, before B powders reach the plasma) and 5.0–5.1 s (red solid lines, after B powders reach the plasma).Identification of B line emissions were conducted based on the wavelength database of NIST [16] for BI–BV lines.The wavelengths of the blended lines are noted in parentheses.

    As shown in figure 3, clear emission lines of BIV 60.31 ? and BV 38.87 ?, 41.00 ?, and 48.59 ? were observed in the 5–60 ? wavelength range measured by the ‘EUV Short’spectrometer.The 100–300 ? wavelength range measured by the‘EUV Long’ spectrometer included higher-order emissions of BIV and BV with relatively strong intensities.In the 300–1050 ? wavelength range measured by the ‘VUV 109L’ spectrometer,many weak lines such as BII (882.50, 882.55, 882.68) ?, BIII(465.55,465.62)?,(510.77,510.85,510.86)?,(677.00,677.14,677.16) ?, and (758.48, 758.67) ?, BIV (384.95, 384.98,385.03) ? and (406.78, 406.81, 406.87) ?, and BV (749.41,749.74, 749.84) ? were observed.In the 970–1870 ? wavelength range measured by the ‘VUV 106R’ spectrometer, the most remarkable lines were BII 1362.46 ? and (1623.60,1623.79, 1623.95, 1624.02, 1624.17, 1624.38) ?.In the 1510–2400 ? wavelength range measured by the ‘VUV 102R’spectrometer,BI(1825.89,1826.40)? and BII 1842.82 ? were observed.At present, the line identifications in the wavelength range of 60–100 ? have not been conducted because this wavelength range was not measured in the discharges which were presented in this paper.We can control the observable wavelength range of the EUV Short spectrometer from 5–60 ? to the longer wavelength range, 40–130 ? at maximum, or that of the EUV Long spectrometer from 100–300 ? to the shorter wavelength range, 60–220 ? at minimum.Therefore, it is possible to measure 60–100 ? by optimizing the wavelength ranges of these two spectrometers in future experiments.

    Figure 3.EUV/VUV spectra including B line emissions observed in a discharge with B powders in(a)5–60 ? measured by the‘EUV Short’spectrometer, (b) 100–300 ? measured by the ‘EUV Long’ spectrometer, (c) 300–1050 ? measured by the ‘VUV 109L’ spectrometer,(d)970–1870 ? measured by the‘VUV 106R’spectrometer,and(e)1510–2400 ? measured by the‘VUV 102R’spectrometer.The spectra were averaged over 3.5–3.6 s (blue, before B powders reach the plasma) and 5.0–5.1 s (red, after B powders reach the plasma).

    Figure 4.Temporal evolutions of BI–BV intensities observed in a discharge with B powders which were integrated over wavelength ranges of (a) 1824.26–1827.82 ? for BI (1825.89, 1826.40) ?,(b) 1360.92–1364.46 ? for BII 1362.46 ?, (c) 674.96–678.34 ? for BIII (677.00, 677.14, 677.16) ?, (d) 120.31–120.94 ? for BIV 60.31×2 ?, and (e) 48.50–48.78 ? for BV 48.59 ?.

    Figure 4 shows temporal evolutions of BI–BV emission intensities observed in a discharge with B powders of which the wavelength spectra were integrated over the wavelength ranges as follows:(a)1824.26–1827.82 ? for BI(1825.89,1826.40)?,(b)1360.92–1364.46 ? for BII 1362.46 ?,(c)674.96–678.34 ? for BIII (677.00, 677.14, 677.16) ?, (d) 120.31–120.94 ? for BIV 60.31×2 ?, and (e) 48.50–48.78 ? for BV 48.59 ?.Background levels of the signals, which mainly consist of bremsstrahlung emission, were subtracted from the raw signals to obtain the intensities of the line emission.The values of ionization potentialEiof BI–BV are 8,25,38,259,and 340 eV,respectively.The electron temperature at the location from which each line is emitted has a value related toEiof each line.Therefore,the BI,BII,and BIII lines are emissions from B ions distributed in the edge plasmas outside the last closed flux surface(LCFS)or divertor plasmas with low electron temperature.On the other hand, the BIV and BV emission lines are considered to originate inside the LCFS.As shown in figure 4,the BIV and BV signals exhibited steep spikes in 4.0–4.5 s coincident with the spikes inPradfollowed by a gradual increase in the signals of all charge states from 4.7 s together with increases inPrad,ne, andne(a99), and a slight decrease inWp.Even though the reason why BIV and BV responded to the powder injection earlier than BI–BIII has not been clarified yet,one possible interpretation is that powders first penetrated into the plasmas through the LCFS and then ablated; thereafter, B ions were transported toward the edge plasmas.The detailed mechanisms of the powder ablation and ion transport processes should be investigated in future studies for a comprehensive understanding of the behaviors of the injected B ions.The deposition depth and ablation profile of the powders are very important information for the impurity powder injection experiment.For more detailed discussion, combination of the line intensity measurement together with the spatial profile measurement,which will be described in section 4.1,is required.

    3.2.Spectra in the boron nitride powder dropping experiments

    Figure 5.Typical waveform of a BN powder dropping experiment in LHD: (a) heating power of ECH and NBI, and gas-puffing duration of H2, (b) amplitude of accelerometer signal Vacc and attenuation of flowmeter signal ΔF/F of the IPD, (c) central electron temperature Te0,(d)edge electron temperature Te(a99),(e)line-averaged electron density ne, (f)edge electron density ne(a99),(g)total radiated power Prad, and (h) plasma stored energy Wp.

    Figure 5 shows a typical waveform of a BN powder dropping experiment in the LHD.The physical quantities plotted in this figure, the magnetic configuration, the heating and fueling patterns,and the onset and the duration of the voltage applied to the IPD feeder piezoelectric blades are the same as those of the discharge in figure 2.As shown in the figures, signals of the plasma response to the powder injection are observed at 4.2 s such as the increase inne,ne(a99),Prad, andWp.

    Figure 6 shows EUV/VUV spectra including B and N line emissions.The spectral data were averaged over 3.5–3.6 s(blue dashed lines,before BN powders reach the plasma)and 4.5–4.6 s (red solid lines, after BN powders reach the plasma).

    Figure 6.EUV/VUV spectra including B and N line emissions observed in a discharge with BN powders in (a) 5–60 ? measured by the‘EUV Short’ spectrometer, (b) 100–300 ? measured by the ‘EUV Long’ spectrometer, (c) 300–1050 ? measured by the ‘VUV 109L’spectrometer,(d)970–1870 ? measured by the‘VUV 106R’spectrometer,and(e)1510–2400 ? measured by the‘VUV 102R’spectrometer.The spectra were averaged over 3.5–3.6 s (blue, before BN powders reach the plasma) and 4.5–4.6 s (red, after BN powders reach the plasma).

    Identification of N line emissions was conducted for NIII–NVII lines.As shown in figure 6,clear emission lines of NVI 28.79 ? and 29.08 ? and NVII 20.91 ? and 24.78 ? were observed in 5–60 ? measured by the ‘EUV Short’spectrometer.The 100–300 ? wavelength range measured by the‘EUV Long’spectrometer included peaks of NV(209.27,209.31) ? and (247.56, 247.71, 247.72) ?.In the 300–1050 ? wavelength range measured by the‘VUV 109L’spectrometer, NIII (989.79, 991.51, 991.58) ? and NIV 765.15 ? were clearly observed.In order to consider whether the NIV 765.15 ? are blended with NIII 763.34 ? and 764.36 ? or not, we compared the temporal evolution of the intensity of the NIV 765.15 ? to those of NIV (921.99,922.52, 923.06, 923.22, 923.68, 924.28) ? observed in the 300–1050 ? wavelength range and NIV 1718.55 ? observed in the wavelength ranges of both 970–1870 ? measured by the‘VUV 106R’spectrometer and 1510–2400 ? measured by the‘VUV 102R’spectrometer.Since the temporal evolutions of these NIV lines were quite similar to each other,we judged that the NIV 765.15 ? has a low possibility to be blended with other lines.In the 970–1870 ? wavelength range, NV(1238.82, 1242.80) ? were also observed.In the 1510–2400 ? wavelength range, NVI 1896.80 ? and NVI(1907.30, 1907.90) ? were observed with relatively strong intensities.

    Figure 7 shows temporal evolutions of BI–BV and NIII–NVII emission intensities observed in a discharge with BN powders.The wavelength range of signal integration for B lines in figures 7(a)–(e) is the same as shown in figure 4.Increases in the signals of BIV and BV from 4.2 s were faster than those of BI, BII, and BIII, which is similar to the behavior observed in the B powder dropping discharge as explained in relation to figure 4.The signals for N lines were integrated over wavelength ranges of(f)990.78–993.35 ? for NIII (989.79, 991.51, 991.58) ?, (g) 763.23–766.58 ? for NIV 765.15 ?,(h)208.80–209.59 ? for NV(209.27,209.31)?, (i) 1894.29–1897.83 ? for NVI 1896.80 ?, and (j)24.65–24.85 ? for NVII 24.78 ?.TheEivalues of NIII–NVII are 47, 77, 98, 552, and 667 eV, respectively.Therefore, the NIII, NIV, and NV lines are emissions from N ions distributed in the edge plasmas outside LCFS while the NVI and NVII emission lines originate inside LCFS.The sharpest increases of the signals seem to be in NIV and NV compared to those in NIII, NVI, and NVII.The sharp increases in NIV and NV signals suggested that the powders might have been ablated outside LCFS followed by the transport of B and N ions to both the higher and the lower sides of the electron temperature.The spike signals in BII and NIII which appeared at 3.0 s in figure 7 are not due to the powder injection.ECH was injected to initiate the discharge at 3.0 s,and then the BII and NIII signals had onset at the same time as the plasma initiation.The emissions from impurity ions in lower charge states such as BII and NIII were observed even in the discharges without the powder injection; therefore, the emission source at this timing was considered to be the impurity ions which were deposited on the surface of the plasma-facing components.Note that the powder size in this BN powder dropping experiment was 60 μm, which was smaller than the B powder dropping experiment with the powder size of 150 μm.If the locations of the powder ablation were actually different between the B and the BN powder dropping experiments analyzed in this paper,the difference in the powder sizes could be one of the reasons for the difference in the locations of the powder ablation.

    Figure 7.Temporal evolutions of BI–BV and NIII-NVII intensities observed in a discharge with BN powders.The wavelength range of signal integration for B lines in (a)–(e) are the same as shown in figure 4.The signals for N lines were integrated over wavelength ranges of (f) 990.78–993.35 ? for NIII (989.79, 991.51, 991.58) ?,(g) 763.23–766.58 ? for NIV 765.15 ?, (h) 208.80–209.59 ? for NV(209.27,209.31)?,(i)1894.29–1897.83 ? for NVI 1896.80 ?,and (j) 24.65–24.85 ? for NVII 24.78 ?.

    Table 1.Useful B and N lines for spectroscopic diagnostics observed in this study.Wavelengths of blended lines are enclosed in parentheses.

    Further studies of the spatial profiles of the emission lines are necessary to explain the difference or similarity in the temporal evolution of emissions from each charge state.Therefore,it is important for the profile measurements to find and summarize lines which are isolated from other lines in the wavelength and have substantially large intensities.The B and N lines observed in this study are summarized in table 1 together withEifor each charge state as useful tools for further spectroscopic diagnostics.Some of the lines are blended with each other;thus,the wavelengths of the blended lines are enclosed in parentheses in the table.

    4.Application of identified lines to advanced spectroscopy

    4.1.Emission intensity profile measurement of BV and NVII lines with space-resolved EUV spectroscopy

    As stated in the previous chapter,the emission lines identified in this study can be used for spatial profile measurements using space-resolved EUV and VUV spectrometers which have been developed in the LHD [11, 12, 17–19].Figure 8 shows the temporal evolution of vertical profiles of emission intensities of (a)–(c) BV 48.59 ? and (d)–(f) NVII 24.78×2 ? measured using the ‘EUV Short2’ spectrometer in a discharge with BN powders.In this profile measurement,the observation chords are aligned vertically in the observation range as illustrated in figure 1(c).The horizontal axis in figure 8 is the position of each observation chord at the major radius ofR=3600 mm.The position of the LCFS in the vacuum magnetic field is also indicated.The BV emission profile has a peak atZ=450 mm as shown in figures 8(a)–(c)while the NVII emission profile has a peak atZ=420 mm.Both of them are located inside the position of LCFS,Z=460 mm.It is also expected that the peak of NVII is located inside that of BV because theEiof NVII (667 eV) is larger than theEiof BV(340 eV).At present,these data of the spatial profiles are line-integral values along the observation chords.Some form of inversion technique such as Abel inversion will be applied in a future analysis to obtain the spatial profiles as local values.Moreover,if the spatial profile measurement can be applied to not only BV and NVII but also various charge states, it should contribute a comprehensive understanding of the spatio-temporal evolution of injected impurities including the powder ablation.

    4.2.B+ ion temperature measurement with high spectral resolution VUV spectroscopy

    The Doppler broadening of impurity emission lines provides information regarding the ion temperature.Figure 9 shows ion temperature measurements of B+ions from Doppler broadening of the BII 1362.46 ? line using a 3 m normal incidence VUV spectrometer.The vertical observation range of the spectrometer covers the full height of the LHD plasmas at the horizontally elongated poloidal cross section as illustrated in figure 1(c).The temporal resolution is 20 ms with a full-binning CCD data acquisition mode.A wavelength spectrum of the BII line which is averaged over a time range of 5.3–5.4 s in a B powder dropping discharge, as shown in figure 2, is plotted in figure 9(a) together with a Gaussian function as a fitting curve.The instrumental function of the spectrometer, ΔINST,was 0.154 ? for the entrance slit width of 50 μm which was employed in the diagnostics.It was determined as the full width at half maximum(FWHM)of the Gaussian-fitted HgI 2536.52 ? spectral profile from a lowpressure Hg lamp (HAMAMATSU L937-02).Since the instrumental function has a comparable value to the observed broadening of the BII spectrum in FWHM,Δobs,for example 0.210 ? as shown in figure 9(a), the signal to noise ratio should be improved as much as possible by averaging the signals over some time ranges to confirm that the observed broadening is significantly larger than the instrumental function.

    Figure 8.Temporal evolution of vertical profiles of emission intensities of(a)–(c)BV 48.59 ? and(d)–(f)NVII 24.78×2 ? measured using the ‘EUV Short2’ spectrometer in a discharge with BN powders.The position of LCFS in the vacuum magnetic field is also indicated.

    Figure 9.Ion temperature measurement of B+ ions using a 3 m normal incidence VUV spectrometer.(a) Wavelength spectrum of BII 1362.46 ? together with a Gaussian function as a fitting curve.Temporal evolutions of (b) edge electron temperature Te(a99) and edge electron density ne(a99), and (c) the temperature of B+ ions,TiBII, and the emission intensity of the BII line.(d) Edge electron density dependence of TiBII.

    The ion temperatureTiin eV is provided byTi=1.68×108M(ΔFWHM/λ0)2, whereMis the atomic mass number,ΔFWHMis the Doppler width at FWHM, and λ0is the central wavelength.As a result derived from the spectrum shown in figure 9(a) withM=10.81, ΔFWHM=(?)1/2=0.143 ?, and λ0=1362.46 ?, the ion temperature of B+ions,TiBII,of 19.9 eV was obtained.The error inTiBII,δTiBII,can be derived by considering a propagation of the fitting error in ΔFWHM, δΔFWHM, which is provided by δTiBII=2TiBIIδΔFWHM/ΔFWHM.In the case in figure 9(a), δTiBII= 5.0 eV was obtained with δΔFWHM=0.018 ?.Because of the lowEivalue,25 eV,of B+ions,theTiBIIvalue obtained here is the ion temperature of B+ions which are distributed in the edge plasmas where the electron temperature is low.Thus, comparison betweenTiBIIand edge plasma parameters contributes to the characterization of the behavior of B+ions.Temporal evolutions ofTe(a99) andne(a99) are shown in figure 9(b) as well asTiBIIand the emission intensity of the BII line in figure 9(c).In this discharge,TiBIIstarted to increase slightly from 5.15 s,whereasne(a99)decreased from 5.35 s.At present,it is not clear whether the time difference between the onset of the increase inTiBIIand the onset of the decrease inne(a99) has some physical meaning or not.The edge electron density dependence ofTiBIIis plotted in figure 9(d)for four discharges with B powder dropping,three discharges with BN powder dropping, and one discharge without powder dropping.As shown in figure 9(d),TiBIIranges around the ionization potential and has a weak negative correlation with the edge electron density.Moreover, in a previous study, a negative correlation between the edge electron density and the impurity ion temperature has been found for CII, CIII,and CIV emission lines released from intrinsic C impurity ions which originate from the sputtering of the C divertor plates in LHD [20].Further investigation of this kind of relationship among edge plasma parameters can contribute to the validation of simulations of the edge plasmas such as three-dimensional impurity transport calculations using EMC3-EIRENE code[21].In addition,the observable wavelength range of the 3 m normal incidence VUV spectrometer,300–3200 ?,covers some N lines such as NIV 765.15 ?, NV (1238.82, 1242.80) ?, and NVI 1896.80 ? which have relatively large intensities as shown in figure 6.They should be the next targets of the impurity ion temperature measurement.

    5.Summary

    An IPD was installed in the LHD under a collaboration between the National Institute for Fusion Science and the Princeton Plasma Physics Laboratory.In order to assess the effectiveness of the injection of B or BN impurity powders,spectroscopic diagnostics were applied to observe line emission from injected impurities.Emission lines released from B and N ions were identified in the EUV wavelength range of 5–300 ? measured using two grazing incidence flat-field EUV spectrometers and in the VUV wavelength range of 300–2400 ? measured using three normal incidence 20 cm VUV spectrometers.BI–BV and NIII–NVII emission lines were identified in discharges with B and BN powder injection,respectively.Useful B and N emission lines which have large intensities and are isolated from other lines were successfully identified as follows: BI (1825.89, 1826.40) ?(blended), BII 1362.46 ?, BIII (677.00, 677.14, 677.16) ?(blended), BIV 60.31 ?, BV 48.59 ?, NIII (989.79, 991.51,991.58) ? (blended), NIV 765.15 ?, NV (209.27, 209.31) ?(blended), NVI 1896.80 ?, and NVII 24.78 ?.Applications of the lines identified to advanced spectroscopic diagnostics were demonstrated, such as vertical profile measurements for the BV and NVII lines using a space-resolved EUV spectrometer and ion temperature measurement for the BII line using a normal incidence 3 m VUV spectrometer.

    Acknowledgments

    The authors thank all the members of the LHD experiment group for their cooperation with the LHD operation.This work is partially supported by the Post-CUP program, JSPSCAS Bilateral Joint Research Projects, ‘Control of wall recycling on metallic plasma facing materials in fusion reactor,’ 2019–2022, (No.GJHZ201984), US Department of Energy (No.DE-AC02-09CH11466) with Princeton University, the LHD project financial support (Nos.ULPP010,ULFF022) and JSPS KAKENHI (Nos.17K14426,20K03896).

    猜你喜歡
    森田直子
    聲音的森林
    上便利商店去
    上便利商店去
    Two-point model analysis of SOL plasma in EAST
    挪威的森林第06章 綠茵藏艷(38)
    挪威的森林第06章 綠茵藏艷(16)
    《挪威的森林》中的女性
    森田化學(xué)提升在華六氟磷酸鋰產(chǎn)能
    暈輪效應(yīng)下《挪威的森林》中直子形象
    改良森田治療法對(duì)改善康復(fù)期精神病患者應(yīng)對(duì)方式的成效
    亚洲av一区综合| 国产精品一及| 一级爰片在线观看| 久久久久性生活片| 简卡轻食公司| 国产高清国产精品国产三级 | 亚洲不卡免费看| 亚洲人与动物交配视频| 一区二区三区乱码不卡18| 男女国产视频网站| 日韩欧美三级三区| 日韩av不卡免费在线播放| 熟女电影av网| 女人久久www免费人成看片 | 全区人妻精品视频| 中文字幕久久专区| 久久久久久久久久久丰满| 中文在线观看免费www的网站| 尤物成人国产欧美一区二区三区| 免费黄色在线免费观看| 精品人妻熟女av久视频| 2021少妇久久久久久久久久久| 欧美成人a在线观看| 久久精品影院6| 免费电影在线观看免费观看| 我的女老师完整版在线观看| 色噜噜av男人的天堂激情| 男人舔女人下体高潮全视频| 国产真实伦视频高清在线观看| 日本与韩国留学比较| 亚洲精品乱码久久久久久按摩| 久久久久久久久久成人| 亚洲丝袜综合中文字幕| 国产精品电影一区二区三区| 中文字幕久久专区| av线在线观看网站| 精品国内亚洲2022精品成人| 欧美变态另类bdsm刘玥| 国产探花在线观看一区二区| 国产综合懂色| 丝袜喷水一区| 欧美极品一区二区三区四区| 亚洲美女搞黄在线观看| 嘟嘟电影网在线观看| 热99re8久久精品国产| 欧美xxxx黑人xx丫x性爽| 欧美成人一区二区免费高清观看| videossex国产| 国产在视频线精品| 春色校园在线视频观看| 成年女人永久免费观看视频| 国产精品日韩av在线免费观看| 嫩草影院精品99| 一级二级三级毛片免费看| 国产亚洲最大av| 免费不卡的大黄色大毛片视频在线观看 | 久久热精品热| 国产精品综合久久久久久久免费| 边亲边吃奶的免费视频| 免费观看在线日韩| 麻豆久久精品国产亚洲av| 最近2019中文字幕mv第一页| 丝袜喷水一区| 日日干狠狠操夜夜爽| 亚洲成av人片在线播放无| 日本一本二区三区精品| 午夜久久久久精精品| 午夜精品国产一区二区电影 | 亚洲天堂国产精品一区在线| 国产一区有黄有色的免费视频 | 日韩一区二区三区影片| 非洲黑人性xxxx精品又粗又长| 麻豆av噜噜一区二区三区| 国产中年淑女户外野战色| 国产片特级美女逼逼视频| 日本爱情动作片www.在线观看| 乱码一卡2卡4卡精品| 少妇人妻精品综合一区二区| 熟女人妻精品中文字幕| 简卡轻食公司| 老司机福利观看| 亚洲av成人精品一区久久| 亚洲,欧美,日韩| 岛国毛片在线播放| 国产精品国产高清国产av| 免费无遮挡裸体视频| 欧美激情久久久久久爽电影| av线在线观看网站| 亚洲美女视频黄频| 人妻夜夜爽99麻豆av| 日本五十路高清| 婷婷色综合大香蕉| 日韩大片免费观看网站 | 最近手机中文字幕大全| 一本久久精品| 在现免费观看毛片| 欧美丝袜亚洲另类| 99热6这里只有精品| 最后的刺客免费高清国语| 亚洲精品国产成人久久av| 亚洲av成人av| 国产成年人精品一区二区| 超碰av人人做人人爽久久| 一区二区三区乱码不卡18| 超碰97精品在线观看| 亚洲图色成人| 国产探花在线观看一区二区| 久久久久久久久中文| 国产精品不卡视频一区二区| 亚洲国产精品合色在线| h日本视频在线播放| 乱码一卡2卡4卡精品| 99热这里只有精品一区| 十八禁国产超污无遮挡网站| 亚洲精华国产精华液的使用体验| 日韩国内少妇激情av| 国产一级毛片七仙女欲春2| 国产乱人视频| 男女下面进入的视频免费午夜| 久久人人爽人人爽人人片va| 一个人免费在线观看电影| 免费大片18禁| 人妻少妇偷人精品九色| 中文字幕亚洲精品专区| 免费黄网站久久成人精品| 秋霞伦理黄片| 国产白丝娇喘喷水9色精品| 久久久午夜欧美精品| 亚洲一级一片aⅴ在线观看| 色综合站精品国产| 午夜免费激情av| 日本一二三区视频观看| 欧美xxxx性猛交bbbb| 熟女人妻精品中文字幕| 99热6这里只有精品| 国产日韩欧美在线精品| 国产成人a区在线观看| 中文字幕精品亚洲无线码一区| 国产成人aa在线观看| 美女高潮的动态| 久久久久久久亚洲中文字幕| 国产白丝娇喘喷水9色精品| 夜夜看夜夜爽夜夜摸| 91午夜精品亚洲一区二区三区| 亚洲怡红院男人天堂| 又粗又硬又长又爽又黄的视频| 中国国产av一级| 日韩av在线大香蕉| 禁无遮挡网站| 在线免费十八禁| 日本黄色片子视频| 天堂影院成人在线观看| 欧美潮喷喷水| 久久精品国产亚洲av涩爱| 欧美三级亚洲精品| 美女大奶头视频| 日韩中文字幕视频在线看片| 亚洲精品久久午夜乱码| 中文字幕最新亚洲高清| 国产免费一区二区三区四区乱码| 国产精品久久久久久精品电影小说| 999精品在线视频| av在线app专区| 色5月婷婷丁香| 久久久亚洲精品成人影院| 超色免费av| 中文精品一卡2卡3卡4更新| 校园人妻丝袜中文字幕| 午夜激情av网站| 乱人伦中国视频| 80岁老熟妇乱子伦牲交| 亚洲国产日韩一区二区| 男女啪啪激烈高潮av片| 爱豆传媒免费全集在线观看| 人人澡人人妻人| 美国免费a级毛片| 久久久久久人人人人人| 国产又爽黄色视频| 不卡视频在线观看欧美| 永久网站在线| 国产欧美日韩综合在线一区二区| 又粗又硬又长又爽又黄的视频| 免费女性裸体啪啪无遮挡网站| 亚洲一区二区三区欧美精品| 成人手机av| 国产精品人妻久久久影院| 人人妻人人添人人爽欧美一区卜| 99久国产av精品国产电影| 黄色 视频免费看| 十八禁高潮呻吟视频| 少妇的逼水好多| 亚洲经典国产精华液单| 亚洲精品自拍成人| 成人国语在线视频| 日本爱情动作片www.在线观看| 中文字幕精品免费在线观看视频 | 男女免费视频国产| 五月伊人婷婷丁香| 成人影院久久| videossex国产| 国产成人欧美| 热99久久久久精品小说推荐| 国产熟女欧美一区二区| 亚洲综合色惰| 国产成人免费无遮挡视频| 久久人人爽人人片av| 亚洲第一av免费看| 尾随美女入室| 国产成人精品无人区| 日韩精品有码人妻一区| 日韩大片免费观看网站| 成人漫画全彩无遮挡| 国产 精品1| 精品人妻一区二区三区麻豆| 欧美日本中文国产一区发布| 成人毛片a级毛片在线播放| 大码成人一级视频| 丰满少妇做爰视频| 啦啦啦啦在线视频资源| 久久国产精品男人的天堂亚洲 | 久久人人爽av亚洲精品天堂| 婷婷色综合www| 日韩制服丝袜自拍偷拍| 国产欧美亚洲国产| 大香蕉久久成人网| 国产片内射在线| 欧美日韩av久久| 国产精品久久久av美女十八| 色5月婷婷丁香| 蜜臀久久99精品久久宅男| 国国产精品蜜臀av免费| 欧美 亚洲 国产 日韩一| 纯流量卡能插随身wifi吗| 人妻一区二区av| 亚洲一区二区三区欧美精品| 午夜久久久在线观看| 高清欧美精品videossex| 亚洲美女搞黄在线观看| www.av在线官网国产| 成人黄色视频免费在线看| 91久久精品国产一区二区三区| 国产精品熟女久久久久浪| 在线观看www视频免费| 久久久国产欧美日韩av| 精品国产露脸久久av麻豆| 最黄视频免费看| 人人澡人人妻人| 热re99久久精品国产66热6| 91精品三级在线观看| av播播在线观看一区| 综合色丁香网| 丝袜脚勾引网站| 天天影视国产精品| 99精国产麻豆久久婷婷| 99热网站在线观看| xxxhd国产人妻xxx| 亚洲精品成人av观看孕妇| 成年女人在线观看亚洲视频| 日本黄大片高清| 日韩不卡一区二区三区视频在线| 看免费av毛片| 欧美 日韩 精品 国产| 亚洲国产精品999| 亚洲av日韩在线播放| 精品久久蜜臀av无| 观看美女的网站| 久久精品夜色国产| 久久精品aⅴ一区二区三区四区 | 夜夜爽夜夜爽视频| 精品国产一区二区三区久久久樱花| 亚洲欧美成人综合另类久久久| 日本av手机在线免费观看| 亚洲五月色婷婷综合| 亚洲人成网站在线观看播放| av网站免费在线观看视频| 97在线视频观看| 日本欧美国产在线视频| 国产精品久久久av美女十八| 亚洲av欧美aⅴ国产| 丰满乱子伦码专区| 国产精品一二三区在线看| 亚洲丝袜综合中文字幕| 日韩人妻精品一区2区三区| 国产精品免费大片| 中文乱码字字幕精品一区二区三区| 久久狼人影院| 婷婷色麻豆天堂久久| 国产视频首页在线观看| 亚洲成人一二三区av| 老司机影院毛片| 一本—道久久a久久精品蜜桃钙片| 成年人午夜在线观看视频| 成年动漫av网址| 纵有疾风起免费观看全集完整版| 激情五月婷婷亚洲| 亚洲av综合色区一区| 90打野战视频偷拍视频| 精品久久蜜臀av无| 曰老女人黄片| 高清视频免费观看一区二区| 亚洲精品视频女| 日本av手机在线免费观看| 18禁观看日本| 最近中文字幕2019免费版| 国产亚洲午夜精品一区二区久久| 国产日韩欧美视频二区| 久久国产精品男人的天堂亚洲 | 国产黄色视频一区二区在线观看| 另类亚洲欧美激情| 在线观看免费日韩欧美大片| 日韩成人av中文字幕在线观看| 亚洲国产成人一精品久久久| 国产熟女欧美一区二区| 一级毛片电影观看| 晚上一个人看的免费电影| 免费高清在线观看日韩| 一级a做视频免费观看| 麻豆精品久久久久久蜜桃| 国产不卡av网站在线观看| 免费观看a级毛片全部| 国产亚洲精品久久久com| a级片在线免费高清观看视频| 欧美精品一区二区免费开放| 99久久综合免费| 国产免费又黄又爽又色| 2022亚洲国产成人精品| 亚洲欧美日韩卡通动漫| 男人舔女人的私密视频| 成年美女黄网站色视频大全免费| 久久久精品94久久精品| 国产精品蜜桃在线观看| 最近中文字幕2019免费版| 国产精品欧美亚洲77777| 国产精品久久久久久久久免| 久久精品人人爽人人爽视色| 免费不卡的大黄色大毛片视频在线观看| 日韩一区二区三区影片| 久久久久久久久久久免费av| 十分钟在线观看高清视频www| 最近的中文字幕免费完整| 久久精品熟女亚洲av麻豆精品| 91午夜精品亚洲一区二区三区| 国产av国产精品国产| av在线老鸭窝| 只有这里有精品99| 大片免费播放器 马上看| 久久狼人影院| 超色免费av| 精品久久久精品久久久| 考比视频在线观看| 久久精品国产自在天天线| 国产亚洲精品久久久com| 久久亚洲国产成人精品v| 美女主播在线视频| 婷婷色综合www| 中国美白少妇内射xxxbb| 国产国语露脸激情在线看| 午夜福利影视在线免费观看| 考比视频在线观看| 久久久久精品久久久久真实原创| 热99国产精品久久久久久7| 国产av一区二区精品久久| 日本与韩国留学比较| 久久韩国三级中文字幕| 夜夜爽夜夜爽视频| 欧美精品一区二区免费开放| 咕卡用的链子| a级片在线免费高清观看视频| 99热国产这里只有精品6| 99久国产av精品国产电影| 亚洲性久久影院| 亚洲成国产人片在线观看| 老司机影院成人| 欧美老熟妇乱子伦牲交| 免费在线观看黄色视频的| 2021少妇久久久久久久久久久| av视频免费观看在线观看| 欧美人与性动交α欧美精品济南到 | av不卡在线播放| 赤兔流量卡办理| 哪个播放器可以免费观看大片| 午夜老司机福利剧场| 国产av一区二区精品久久| av福利片在线| 人妻人人澡人人爽人人| 精品人妻在线不人妻| 欧美人与性动交α欧美精品济南到 | 在线天堂中文资源库| 色吧在线观看| 日韩一区二区三区影片| 人人妻人人澡人人爽人人夜夜| 日韩av免费高清视频| 一级黄片播放器| 精品国产一区二区三区久久久樱花| 免费av中文字幕在线| 国产爽快片一区二区三区| 交换朋友夫妻互换小说| 久久 成人 亚洲| 午夜影院在线不卡| 91国产中文字幕| 亚洲高清免费不卡视频| 久热这里只有精品99| av福利片在线| 亚洲欧美中文字幕日韩二区| 久久亚洲国产成人精品v| 爱豆传媒免费全集在线观看| 少妇人妻精品综合一区二区| 性高湖久久久久久久久免费观看| 久久人人爽人人爽人人片va| 天天躁夜夜躁狠狠躁躁| 永久免费av网站大全| 黑人欧美特级aaaaaa片| 涩涩av久久男人的天堂| 亚洲一区二区三区欧美精品| 久久 成人 亚洲| 成年女人在线观看亚洲视频| 免费观看a级毛片全部| 国产免费一区二区三区四区乱码| 久久狼人影院| 少妇被粗大猛烈的视频| 春色校园在线视频观看| 免费在线观看完整版高清| 亚洲伊人久久精品综合| 久久精品夜色国产| 九九在线视频观看精品| 成人影院久久| videos熟女内射| 五月伊人婷婷丁香| 黄色视频在线播放观看不卡| 成年人免费黄色播放视频| 亚洲av男天堂| 啦啦啦啦在线视频资源| 日本黄大片高清| 免费观看无遮挡的男女| 成人午夜精彩视频在线观看| 99久久中文字幕三级久久日本| 欧美日韩亚洲高清精品| 欧美精品人与动牲交sv欧美| 高清不卡的av网站| 国产成人精品一,二区| 最黄视频免费看| 成年人午夜在线观看视频| 精品少妇黑人巨大在线播放| 久久人人爽人人爽人人片va| 天天躁夜夜躁狠狠躁躁| 一级,二级,三级黄色视频| 久久国内精品自在自线图片| 免费大片黄手机在线观看| 夫妻午夜视频| 国产在线一区二区三区精| 只有这里有精品99| 欧美人与性动交α欧美软件 | 观看美女的网站| 高清在线视频一区二区三区| 国产精品国产三级专区第一集| 性色avwww在线观看| 亚洲久久久国产精品| 亚洲欧洲日产国产| 久久精品国产a三级三级三级| 免费人妻精品一区二区三区视频| 在现免费观看毛片| 色网站视频免费| 校园人妻丝袜中文字幕| 精品人妻一区二区三区麻豆| 一区在线观看完整版| 日韩 亚洲 欧美在线| 亚洲第一区二区三区不卡| 久久精品aⅴ一区二区三区四区 | 久久久久久久大尺度免费视频| freevideosex欧美| 久久久久精品性色| av.在线天堂| 久久午夜福利片| 亚洲在久久综合| 伦理电影免费视频| 一级爰片在线观看| 26uuu在线亚洲综合色| 麻豆乱淫一区二区| 男的添女的下面高潮视频| 99久久精品国产国产毛片| 欧美最新免费一区二区三区| 亚洲精品久久午夜乱码| 久久亚洲国产成人精品v| 国产激情久久老熟女| 国产亚洲av片在线观看秒播厂| 2022亚洲国产成人精品| 亚洲精品成人av观看孕妇| 亚洲一区二区三区欧美精品| 亚洲人与动物交配视频| 美女大奶头黄色视频| 婷婷色av中文字幕| 日韩av在线免费看完整版不卡| 高清欧美精品videossex| 久久精品人人爽人人爽视色| 久久久久久久精品精品| 日韩在线高清观看一区二区三区| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲伊人色综图| 日本色播在线视频| 久久精品夜色国产| 国产成人精品无人区| 99九九在线精品视频| 国产白丝娇喘喷水9色精品| 久久99热这里只频精品6学生| 日韩欧美一区视频在线观看| 精品久久久久久电影网| 久久久久久久国产电影| 如日韩欧美国产精品一区二区三区| 狂野欧美激情性bbbbbb| 亚洲欧美一区二区三区黑人 | 两个人看的免费小视频| 赤兔流量卡办理| 一本—道久久a久久精品蜜桃钙片| 国产爽快片一区二区三区| 制服诱惑二区| 热re99久久精品国产66热6| 久久99精品国语久久久| 国产精品国产三级国产专区5o| 欧美精品人与动牲交sv欧美| 大香蕉97超碰在线| 香蕉丝袜av| 最后的刺客免费高清国语| 亚洲性久久影院| 精品第一国产精品| av黄色大香蕉| 国产成人91sexporn| 黑人巨大精品欧美一区二区蜜桃 | 侵犯人妻中文字幕一二三四区| 男女免费视频国产| 高清欧美精品videossex| 香蕉丝袜av| 午夜91福利影院| 熟女电影av网| 久久精品国产综合久久久 | 国产精品女同一区二区软件| 国产又色又爽无遮挡免| 男女国产视频网站| 少妇的逼水好多| 一级毛片电影观看| 国产av一区二区精品久久| 国国产精品蜜臀av免费| 精品视频人人做人人爽| 一级a做视频免费观看| 美女中出高潮动态图| 欧美精品国产亚洲| 精品国产一区二区三区四区第35| 一级片'在线观看视频| 久久精品aⅴ一区二区三区四区 | 久久99蜜桃精品久久| 热re99久久精品国产66热6| 成人黄色视频免费在线看| 18禁动态无遮挡网站| 日本色播在线视频| 久久久久久人人人人人| 免费播放大片免费观看视频在线观看| 自线自在国产av| 搡女人真爽免费视频火全软件| 九九爱精品视频在线观看| 国产国语露脸激情在线看| 久久久久久人妻| 亚洲精品色激情综合| 一级爰片在线观看| 欧美精品国产亚洲| 亚洲国产看品久久| 少妇的逼好多水| 亚洲丝袜综合中文字幕| 免费av中文字幕在线| 尾随美女入室| 黄色怎么调成土黄色| 欧美日本中文国产一区发布| 最近中文字幕高清免费大全6| 99国产精品免费福利视频| 黄色毛片三级朝国网站| 老司机影院毛片| 欧美人与性动交α欧美软件 | 久久国内精品自在自线图片| 男女免费视频国产| 国产精品99久久99久久久不卡 | 国产深夜福利视频在线观看| 国产成人91sexporn| 午夜影院在线不卡| 国产亚洲最大av| 不卡视频在线观看欧美| 久久久久久久精品精品| 国产男女超爽视频在线观看| 亚洲四区av| 日韩一本色道免费dvd| 侵犯人妻中文字幕一二三四区| 日本与韩国留学比较| 男人舔女人的私密视频| 亚洲一级一片aⅴ在线观看| 两个人免费观看高清视频| 日韩成人av中文字幕在线观看| 国产精品.久久久| 咕卡用的链子| 男的添女的下面高潮视频| 久久精品久久久久久久性| 99热国产这里只有精品6| 精品亚洲乱码少妇综合久久| 女人精品久久久久毛片| 久久狼人影院| 国产精品女同一区二区软件| 黄色配什么色好看| 性色avwww在线观看| 成人毛片a级毛片在线播放| 亚洲精品第二区| 亚洲色图 男人天堂 中文字幕 | 久久精品人人爽人人爽视色| 婷婷色综合www| 国产精品三级大全| 宅男免费午夜| 亚洲欧美日韩卡通动漫| 国产av码专区亚洲av| 国产精品一区www在线观看| 91国产中文字幕|