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

    HL-2A 裝置上成像型中性粒子分析器的物理設(shè)計(jì)和初步實(shí)驗(yàn)結(jié)果*

    2023-11-24 05:06:06顏筱宇何小斐于利明劉亮陳偉石中兵盧杰魏會(huì)領(lǐng)韓紀(jì)鋒張軼潑鐘武律許敏
    物理學(xué)報(bào) 2023年21期
    關(guān)鍵詞:診斷系統(tǒng)螺距膜片

    顏筱宇 何小斐 于利明? 劉亮 陳偉 石中兵 盧杰魏會(huì)領(lǐng) 韓紀(jì)鋒 張軼潑 鐘武律 許敏

    1) (四川大學(xué)原子核科學(xué)技術(shù)研究所,輻射物理及技術(shù)教育部重點(diǎn)實(shí)驗(yàn)室,成都 610065)

    2) (核工業(yè)西南物理研究院,成都 610041)

    在HL-2A 裝置上發(fā)展了基于硫化鋅銀閃爍體的成像型中性粒子分析器,對(duì)磁約束聚變等離子體中高能量粒子(EP)的分布、能量和螺距角等關(guān)鍵信息,以及EP 與磁流體不穩(wěn)定性之間的相互作用等物理問題進(jìn)行了研究.在中性束注入路徑上逃逸出的具有等離子體中快離子能量和螺距角信息的中性粒子,通過由入射孔和碳微晶體膜片組成的準(zhǔn)直系統(tǒng)后轉(zhuǎn)化為離子,在裝置邊緣磁場(chǎng)中受洛倫茲力偏轉(zhuǎn)而撞擊到閃爍體屏上.通過分析發(fā)光點(diǎn)的位置和光強(qiáng)度,可以推斷出裝置中快離子的位置、能量和螺距角等關(guān)鍵信息.在HL-2A 裝置高能量粒子物理實(shí)驗(yàn)中,通過該診斷和理論計(jì)算初步證實(shí)了長(zhǎng)壽模不穩(wěn)定性是由能量、螺距角和位置分別為E=12.5—32 keV,θ ~ 149.2° (v///v ~ 0.86)和R=170.5—171.5 cm 的芯部快離子激發(fā).

    1 引言

    磁約束聚變裝置中的大功率輔助加熱和電流驅(qū)動(dòng)系統(tǒng)(如中性束注入(NBI)、離子回旋波加熱、電子回旋波加熱/電流驅(qū)動(dòng)和低混雜波電流驅(qū)動(dòng)(LHCD)等),以及氘氚聚變反應(yīng)會(huì)產(chǎn)生大量的能量偏離等離子體溫度的高能量粒子(EP),包括快離子(FI)、快電子和阿爾法(α)粒子等.一方面,EP 可以有效地提高等離子體的溫度,驅(qū)動(dòng)等離子體無(wú)感電流,維持等離子體燃燒等.另一方面,EP會(huì)和磁流體力學(xué)波通過波-粒子相互作用而產(chǎn)生各種高能量粒子的不穩(wěn)定性[1-5].這些不穩(wěn)定性反過來(lái)會(huì)與EP 相互作用導(dǎo)致其發(fā)生顯著的擴(kuò)散、對(duì)流等輸運(yùn)和損失,進(jìn)而引起輔助加熱和電流驅(qū)動(dòng)效率下降,而攜帶著大量能量的損失EP 會(huì)對(duì)裝置第一壁造成破壞,甚至?xí)?dǎo)致其過熱和融化[6].因此,了解EP 在等離子體中的分布、能譜、輸運(yùn)和損失對(duì)理解等離子體加熱、電流驅(qū)動(dòng)、能量、粒子輸運(yùn)以及約束等至關(guān)重要.為了實(shí)現(xiàn)對(duì)熱離子和快離子的深入研究,國(guó)際上已發(fā)展了中性粒子分析器(NPA)[7]、快離子損失探針[8]、快離子α輻射測(cè)量(FiDα)[9]、法拉第筒[10]和湯姆孫相干散射[11]等診斷系統(tǒng).

    NPA 診斷系統(tǒng)的工作原理是使用準(zhǔn)直系統(tǒng)收集通過電荷交換從等離子體中逃逸出來(lái)的攜帶有離子信息的中性粒子,此后利用氣體或者固體電荷剝離器使中性粒子變成帶電離子,最終使用能量分析器等對(duì)離子的能量、質(zhì)量和通量等進(jìn)行分析,得到離子溫度/能量(E)、徑向位置(R)、通量和螺距角(θ,實(shí)驗(yàn)上一般用v///v表示,即 cosθ=v///v,其中v//表示快離子平行于磁場(chǎng)的速度,v表示離子的總速度)等離子的關(guān)鍵信息.自20 世紀(jì)70 年代T-10 托卡馬克[12]發(fā)展NPA 測(cè)量等離子體中的離子溫度和快離子的速度分布函數(shù)以來(lái),NPA 通過對(duì)各種功能模塊的優(yōu)化,已經(jīng)發(fā)展出了測(cè)量離子溫度(及分布)、快離子的分布和能譜、等離子體中粒子的質(zhì)譜等功能.目前,NPA 已經(jīng)作為一種常規(guī)診斷被廣泛應(yīng)用于托卡馬克和仿星器上[13-18].歐洲聯(lián)合環(huán)(JET)裝置上的E//B(粒子分析功能單元中電場(chǎng)(E)平行于磁場(chǎng)(B)) NPA 不僅可以實(shí)現(xiàn)5—740 keV 范圍的快離子能量測(cè)量,還可得到氫-氘-氚粒子的質(zhì)譜[19].在未來(lái)的ITER 裝置上已經(jīng)設(shè)計(jì)了兩套NPA 診斷系統(tǒng),分別用于測(cè)量燃燒等離子體中氘氚離子的能量和質(zhì)譜,以及能量在MeV量級(jí)的α粒子分布函數(shù)和電荷剝離(He1+和He2+的比例)情況等重要信息,并作為關(guān)鍵診斷來(lái)監(jiān)測(cè)燃燒等離子體的性能[20].國(guó)內(nèi)的HL-1M 裝置發(fā)展了一套基于平行板探測(cè)器的NPA,用于測(cè)量等離子體的平均溫度.HL-2A 裝置沿用了平行板NPA,還開發(fā)了一套11 通道的具有空間分辨能力的緊湊型NPA (CP-NPA)[21],可以對(duì)平均離子溫度分布和非熱離子分布等進(jìn)行測(cè)量.除此之外,還在HL-2A 裝置上發(fā)展了一套電流型固態(tài)NPA (ssNPA)對(duì)能量大于22.5 keV 的快離子進(jìn)行了測(cè)量.

    近些年,美國(guó)Heidbrink 教授的EP 物理研究團(tuán)隊(duì)[22,23]研發(fā)了一種基于閃爍體成像探測(cè),并可以放置在裝置真空室內(nèi)部的成像型NPA (INPA),用于測(cè)量快離子的能譜、空間分布和螺距角等關(guān)鍵信息.相對(duì)于放置于裝置真空室外,并額外需要真空維持、電荷剝離和粒子分析和探測(cè)等復(fù)雜功能模塊的傳統(tǒng)NPA,INPA 具有結(jié)構(gòu)和原理簡(jiǎn)單、造價(jià)較低和實(shí)用性強(qiáng)等優(yōu)勢(shì),已在NSTX,ASDEX-U[24]與EAST[25]等裝置上得到了廣泛應(yīng)用并取得了較好的實(shí)驗(yàn)結(jié)果.

    在HL-2A 裝置上已經(jīng)發(fā)展了一套INPA 診斷系統(tǒng)[26],并在相關(guān)高能量粒子物理實(shí)驗(yàn)中取得了初步的實(shí)驗(yàn)結(jié)果.本文第1 部分(引言)簡(jiǎn)要介紹了NPA 診斷系統(tǒng)及發(fā)展.第2 部分介紹INPA 診斷系統(tǒng)的測(cè)量原理、物理設(shè)計(jì)和主要結(jié)構(gòu).第3 部分介紹INPA 診斷系統(tǒng)在HL-2A 裝置上的安裝布局和初步實(shí)驗(yàn)結(jié)果.第4 部分為總結(jié).

    2 INPA 診斷系統(tǒng)的物理設(shè)計(jì)

    INPA 診斷系統(tǒng)固定于HL-2A 裝置的真空室弱場(chǎng)側(cè)內(nèi)壁,其測(cè)量原理如圖1 所示.由粒子入射孔和具有電荷剝離作用的碳微晶體膜片組成準(zhǔn)直器,與NBI 電荷交換后從背景等離子體中逃逸出攜帶有快離子信息的大量中性粒子通過該準(zhǔn)直系統(tǒng),并被碳微晶體膜片剝離電子重新變成帶電離子(D+),透射出剝離膜片的帶電離子在裝置邊緣磁場(chǎng)中受洛倫茲力(f=ev×B,其中,e為單電荷電量,B為裝置邊緣磁場(chǎng)強(qiáng)度)作用偏轉(zhuǎn)180°后打擊到閃爍體上產(chǎn)生熒光.粒子的能量越高,則偏轉(zhuǎn)半徑(rd=,其中,mi表 示 D+的質(zhì)量)越大,在閃爍體上轟擊點(diǎn)與剝離膜片的縱向距離越大,閃爍體上光斑的亮度反映了轟擊到閃爍體上的粒子通量.從6 個(gè)準(zhǔn)直測(cè)量通道透射出的粒子束及其入射軌跡分別用紅色、藍(lán)色、黃色、綠色、紫色和棕色曲線表示,并分別被命名為第1—6 通道.

    圖1 HL-2A 裝置INPA 診斷系統(tǒng)的主要結(jié)構(gòu)及快離子的測(cè)量軌跡示意圖Fig.1.Structure of the INPA and flight trajectories of FIs on HL-2A.

    圖2 給出了包括了 2#NBI 方向和INPA安裝窗口在內(nèi)的HL-2A 裝置部分俯視圖.在本次放電中,HL-2A 的等離子體電流(Ip)方向?yàn)槟鏁r(shí)針,縱向磁場(chǎng)(Bt)方向?yàn)轫槙r(shí)針.圖2 中用兩條粗黑色直線表示 2#NBI 系統(tǒng)的左右兩個(gè)注入分束,而INPA 系統(tǒng)6 個(gè)測(cè)量通道通過準(zhǔn)直系統(tǒng)測(cè)量到的中性粒子主要來(lái)自于NBI 右注入分束.INPA 系統(tǒng)準(zhǔn)直線在裝置法蘭中平面的理想情況下,第1—6通道具體的測(cè)量位置見圖中紅色、藍(lán)色、黃色、綠色、紫色和棕色粗直線與NBI 右注入分束的交叉點(diǎn).INPA 診斷系統(tǒng)6 個(gè)觀測(cè)位置的中心分別位于大半徑R≈ 172.6,170.6,170.9,175.8,180.5和211.8 cm 處.另外,受到窗口位置和可操作空間等實(shí)際限制因素,INPA 系統(tǒng)安裝在HL-2A 裝置的中平面下方,粒子入射孔在裝置中平面下方Z=-10.5 cm處.

    圖2 HL-2A 裝置上2# NBI 的注入路徑和INPA 診斷系統(tǒng)中6 個(gè)測(cè)量通道所對(duì)應(yīng)的觀測(cè)位置Fig.2.Injection path of 2# NBI system and the observed positions for 6 channels of INPA on HL-2A.

    計(jì)算通過6 個(gè)測(cè)量通道進(jìn)入準(zhǔn)直系統(tǒng)的粒子速度方向與測(cè)量位置處磁場(chǎng)方向之間的夾角即可得到快離子的螺距角.圖3 給出了HL-2A 裝置在第38140 次放電(第3 節(jié)的展示放電)實(shí)驗(yàn)中Bt=1.3 T和Ip=(155±3) kA 的實(shí)際情況下,利用EFIT代碼反演的裝置磁位形數(shù)據(jù),計(jì)算得到了理想情況下快離子的產(chǎn)生位置和對(duì)應(yīng)的螺距角信息,其螺距角v///v≈ 0,0.53,0.86,0.98,0.99 和0.93,即θ ≈ 90.0°,121.9°,149.2°,170.9°,172.0°和158.4° (v///v=cosθ).根據(jù)HL-2A 裝置和INPA系統(tǒng)安裝位置的工程信息,以及38140 次放電的實(shí)驗(yàn)參數(shù),可以得知該INPA 系統(tǒng)在理想情況下所測(cè)量粒子的位置范圍和螺距角等信息,如表1 所列.

    表1 INPA 診斷系統(tǒng)的6 個(gè)測(cè)量通道所觀測(cè)粒子的位置和粒子特征信息Table 1.Observed positions and characteristic information of particles from the 6 channels of the INPA system.

    圖3 在HL-2A 裝置極向截面顯示的INPA 系統(tǒng)觀測(cè)到的粒子位置和螺距角Fig.3.Positions and pitch angles of the observed particles from INPA system in the poloidal cross section in HL-2A.

    INPA 診斷系統(tǒng)原理圖(圖1)對(duì)應(yīng)的實(shí)物結(jié)構(gòu)如圖4 所示.INPA 的外部結(jié)構(gòu)為304 不銹鋼制成的長(zhǎng)方形封閉腔體,外部尺寸為10.6 cm ×12.0 cm × 9.5 cm,其封閉結(jié)構(gòu)可以屏蔽等離子體中的各種雜散光和雜質(zhì)粒子等.在INPA 外部腔體上,設(shè)有一個(gè)面向裝置等離子體的直徑φ=3 cm的粒子入射孔,在INPA 腔室內(nèi)入射孔后方3 cm處則水平排布有6 個(gè)間距為1.2 cm 的電荷剝離膜片,從觀測(cè)孔望去從右向左將其編號(hào)為1—6.電荷剝離模塊采用的是10 nm 厚、φ=5 mm 的圓形碳微晶體膜片.為了增加該膜片的耐久性,不僅將膜片固定在內(nèi)/外徑1.0/0.5 cm、厚0.5 mm 的無(wú)氧銅環(huán)上,而且其底部還固定在了每英寸300 格的密銅網(wǎng)上,如圖4(c)所示.碳微晶體膜片的厚度會(huì)使入射粒子的能量產(chǎn)生一定的能量衰減(δE),δE、入射粒子的能量(E)和膜片厚度(d)滿足公式 δE≈4.4×10-2√,其中E和d分別以keV和nm為單位.HL-2A 裝置上NBI 粒子的最高能量約為42 keV,因此當(dāng)膜的厚度為10 nm 時(shí),粒子經(jīng)過剝離膜片后產(chǎn)生的最大能量衰減 δE=2.85 keV.

    圖4 HL-2A 裝置上INPA 診斷系統(tǒng)主要部件的實(shí)物及內(nèi)部布置圖 (a) INPA 的外觀圖;(b)內(nèi)部剝離膜片和閃爍體的布局圖;(c)碳微晶體膜片尺寸和結(jié)構(gòu)Fig.4.External figure and arrangement inside the chamber of INPA diagnostics on HL-2A: (a) External figure;(b) arrangement of carbon microcrystal diaphragm and scintillator inside the chamber: (c) detail structure of carbon microcrystal diaphragm.

    硫化鋅摻雜銀(ZnS(Ag))閃爍體可以發(fā)射波長(zhǎng)峰值為450 nm 的藍(lán)色可見光,且發(fā)光效率約為蒽晶體的300%,光衰減時(shí)間約為200 ns,能夠在中性粒子通量較低時(shí)提高INPA 的信噪比,所以INPA 內(nèi)部顯示粒子轟擊位置和通量等信息的閃爍屏采用了ZnS(Ag) 閃爍體.基于HL-2A 裝置的常規(guī)運(yùn)行磁場(chǎng)強(qiáng)度和約42 keV 的離子能量上限,ZnS(Ag) 閃爍體的尺寸為7.0 cm × 10.0 cm.

    由于直徑為φ=5 mm 的圓形碳微晶體膜片尺寸相對(duì)較大(不能被看作理想的點(diǎn)),因此能量相同、入射方向稍有不同的粒子在入射并擊打到熒光屏的過程中會(huì)存在如圖5(a)所示的發(fā)散情況,即在橫向和縱向形成具有一定寬度的橢圓形轟擊斑.轟擊斑的尺寸隨入射粒子能量以及測(cè)量通道的變化如圖5(b)所示,并將橢圓型光斑在縱向的短軸長(zhǎng)度對(duì)應(yīng)的能量范圍稱作能量分辨率 ΔE.6 個(gè)通道上粒子的落點(diǎn)與剝離膜片中心的距離D(D=2rd)和粒子入射能量的關(guān)系如圖5(c)所示,研究發(fā)現(xiàn)二者成指數(shù)相關(guān).對(duì)2—42 keV 能量范圍的粒子在6 個(gè)通道中的能量分辨率計(jì)算結(jié)果如圖5(d)所示.在HL-2A 的INPA 系統(tǒng)內(nèi),經(jīng)第1 測(cè)量通道對(duì)入射粒子能量分辨率最差,在測(cè)量粒子能量為42 keV 的情況下 ΔE ≈4.5 keV.另外,由于碳微晶體膜片和入射孔的實(shí)際尺寸會(huì)造成準(zhǔn)直系統(tǒng)的觀測(cè)位置偏離理想位置的展寬,如圖2 中理想觀測(cè)粗線兩側(cè)的細(xì)線所示.由此造成的各通道所測(cè)量的最大半徑Rmax、最小半徑Rmin、以及立體張角φ見表2 所列.

    表2 INPA 診斷系統(tǒng)的6 個(gè)通道對(duì)應(yīng)的測(cè)量范圍Table 2.Measurement ranges corresponding to the 6 channels of the INPA diagnostic system.

    圖5 INPA 診斷系統(tǒng)中幾何機(jī)構(gòu)引起的誤差分析 (a)粒子束在磁場(chǎng)中的偏轉(zhuǎn)及在閃爍體上的轟擊斑;(b) INPA 診斷系統(tǒng)6 個(gè)測(cè)量通道的粒子在閃爍體上的落點(diǎn);(c)粒子在閃爍體上的落點(diǎn)位置和入射能量的關(guān)系;(d)能量分辨率與粒子能量的關(guān)系Fig.5.Analysis of errors caused by geometric mechanisms of diagnostic systems: (a) Flight orbits and impact spots of the measured particles on scintillator;(b) positions of impact spots from the particles from 6 channels in INPA;(c) relationship between the position of the particle's landing point on the scintillator and the incident energy;(d) relationship between energy of particles and energy resolution.

    INPA 系統(tǒng)在真空室內(nèi)的實(shí)際安裝位置如圖6(b)所示.由于理想的中平面窗口已經(jīng)被其他診斷占用,INPA 所在直徑為63 cm 的法蘭位于中平面下15 cm 處.經(jīng)過安裝和核實(shí),粒子入射孔位于中平面下10.5 cm 處.在真空室外,通過成像光纖束將真空室內(nèi)閃爍體上的轟擊亮斑圖像信息傳輸給高速相機(jī),如圖6(a)所示.實(shí)驗(yàn)中采用了量子效率峰值與ZnS(Ag)閃爍體發(fā)光波長(zhǎng)峰值較為符合的工業(yè)高速相機(jī)(LUMENERA Lt225)記錄閃爍體的發(fā)光信息,其曝光時(shí)間范圍和最大全畫幅幀頻可以在26μs-500ms和170fps 范圍內(nèi)根據(jù)實(shí)驗(yàn)情況(如 NBI 束線強(qiáng)度、逃逸粒子的通量和閃爍體的發(fā)光強(qiáng)度等)進(jìn)行調(diào)節(jié).

    圖6 HL-2A 裝置上INPA 診斷系統(tǒng)的安裝 (a) INPA 診斷系統(tǒng)在真空室內(nèi)的安裝位置;(b) INPA 診斷系統(tǒng)在真空室外的高速相機(jī)、光纖束和法蘭等Fig.6.Installation of INPA system on HL-2A: (a) Installation of the INPA on the flange inside the vacuum chamber;(b) arrangement of the fast speed camera,light fiber bundle and flange.

    3 INPA 診斷系統(tǒng)的物理實(shí)驗(yàn)結(jié)果及分析

    在HL-2A 裝置中有NBI 加熱提供中性粒子媒介的情況下,INPA 診斷系統(tǒng)被應(yīng)用于快離子與高能量不穩(wěn)定性之間的相互作用研究,確定激發(fā)各種不穩(wěn)定性快離子的特征.在第38140 次放電實(shí)驗(yàn)中,在INPA 診斷系統(tǒng)上得到了有效的實(shí)驗(yàn)數(shù)據(jù),并確定在這次放電實(shí)驗(yàn)中激發(fā)了長(zhǎng)壽模(LLM)不穩(wěn)定性快離子的能量、位置和螺距角等特征信息.

    第38140 次放電中主要的實(shí)驗(yàn)參數(shù)、等離子體運(yùn)行狀態(tài)和快離子激發(fā)的不穩(wěn)定性等見圖7.在這次放電中,等離子體電流和縱向磁場(chǎng)強(qiáng)度分別為Ip=(155 ± 3) kA和Bt=1.3 T,它們隨時(shí)間的演化如圖7(a)中紅色和藍(lán)色曲線所示.1#NBI,2#NBI和LHCD 功率分別為PNBI1=0.42 MW,PNBI2=0.38 MW和PLHCD=0.6 MW,2#NBI 注入粒子的能量Eb≈30—32 keV,其加熱時(shí)序如圖7(b)所示.在NBI 加熱期間的850—1100 ms 和1500—1750 ms 時(shí)間段內(nèi),等離子體處于高約束運(yùn)行模式(H 模),Dα輻射本底突然降低并且出現(xiàn)了邊緣局域模(ELM)引起的周期性振蕩;而在1100—1500 ms期間,等離子體處于低約束運(yùn)行模式(L 模),Dα輻射本底相對(duì)較高,如圖7(c)所示.由于等離子體在H 模運(yùn)行情況下的約束性能優(yōu)于L 模,因此等離子體平均密度(ne)從L 模情況下的1.8 × 1019m-3上升到H 模情況下的2.2×1019—2.5×1019m-3,如圖7(a)中綠色曲線所示.在整個(gè)NBI 加熱的800—1800 ms 期間可在米爾諾夫(Mirnov)磁探針信號(hào)上觀測(cè)到由快離子激發(fā)的高能量不穩(wěn)定性中常見的LLM,LLM 在Mirnov 磁探針信號(hào)上引起的明顯磁擾動(dòng)信號(hào)和其頻率譜圖分別如圖7(d)和圖7(e)所示,LLM 的頻率(fLLM)基本穩(wěn)定在10—12 kHz 之間,并且在L 模運(yùn)行期間LLM 引起的磁擾動(dòng)基本處于較高幅度且飽和狀態(tài);而在H 模運(yùn)行期間,LLM 會(huì)被周期性爆發(fā)的ELM 影響而出現(xiàn)磁擾動(dòng)幅度的明顯減小甚至被致穩(wěn)的狀態(tài),表現(xiàn)在頻率上則是被周期性地打斷和消失.

    圖7 HL-2A 裝置上第38140 次放電的實(shí)驗(yàn)參數(shù)及觀測(cè)到的LLM 不穩(wěn)定性 (a) 等離子體主要放電參數(shù),即 Ip、等離子體平均密度 ne和Bt ;(b) 1#和2# NBI 束線的加熱功率和時(shí)序;(c)氘 α (Dα)輻射信號(hào);(d) Mirnov 磁探針信號(hào)及(e)頻率譜圖Fig.7.Discharge parameters and the observed LLM instabilities in shot 38140 on HL-2A: (a) Main discharge parameters,Ip,lineaveraged electron density ne and Bt ;(b) heating power of 1# and 2# NBI systems and evolution;(c) Dα signal;(d) Mirnov signal and (e) its spectrogram.

    通過調(diào)節(jié)快速相機(jī)記錄了NBI 加熱期間與LLM 引起的與快離子分布等相關(guān)的實(shí)驗(yàn)圖像,如圖8 所示.圖8(a)—(l)給出了快速相機(jī)記錄的在1096—1504 ms 期間內(nèi)與12 個(gè)時(shí)間點(diǎn)相對(duì)應(yīng)的粒子轟擊到閃爍體屏幕上的圖像數(shù)據(jù),其中左上角位置的圓形亮點(diǎn)是等離子體放電發(fā)出的通過入射孔與第1 個(gè)電荷剝離膜片直接被高速相機(jī)記錄的可見光,左上角藍(lán)色圓點(diǎn)為第3 個(gè)剝離膜片的實(shí)際位置.在LLM 擾動(dòng)幅度相對(duì)較強(qiáng)且等離子體約束性能相對(duì)較差的L 模運(yùn)行期間,在閃爍屏右側(cè)觀測(cè)到了與LLM 相關(guān)的明顯的粒子轟擊斑,如圖8(b)—(j)所示.而在LLM 擾動(dòng)幅度較弱的H 模運(yùn)行期間未在閃爍屏上觀測(cè)到轟擊斑,如圖8(a)和圖8(l)所示.在L 模運(yùn)行模式下與LLM 不穩(wěn)定性有關(guān)的快離子可以在INPA 閃爍體上被明顯地觀測(cè)到的原因可能是: 1) H 模運(yùn)行期間裝置對(duì)等離子體整體的約束水平要遠(yuǎn)高于L 模,如ne在H 模運(yùn)行期間高于L 模,裝置對(duì)快離子的約束能力也相對(duì)較高;2)在H 模運(yùn)行期間的高ne情況下,逃逸的中性粒子與等離子體碰撞再次被電離的概率遠(yuǎn)高于L 模運(yùn)行低ne情況;3) LLM 不穩(wěn)定性的強(qiáng)度(Mirnov磁探針的擾動(dòng)幅度)在L 模運(yùn)行模式下要高于H 模運(yùn)行情況,并且在L 模運(yùn)行期間LLM 不會(huì)被頻繁出現(xiàn)的ELM 打斷.在這次實(shí)驗(yàn)中,由于 2#NBI的注入功率和束強(qiáng)度很低,中性粒子束注入等離子體較淺,能為快離子電荷交換逃逸提供的媒介中性粒子較少,因此能打擊到閃爍體上的粒子數(shù)目較少,轟擊斑光強(qiáng)較弱.采用光纖束對(duì)閃爍體閃光圖像進(jìn)行傳輸?shù)姆桨笗?huì)使光通量損失約30%.根據(jù)閃爍體的粒子轟擊閃光效果,高速相機(jī)的曝光時(shí)間和采樣時(shí)間分別被調(diào)整約為30 ms和37 ms.

    圖8 NBI 期間的LLM 引起的粒子轟擊圖像 (a)—(l) LLM 在H 模和L 模運(yùn)行期間的不同時(shí)刻在INPA 閃爍體上觀測(cè)到的粒子轟擊圖像及演化Fig.8.Impact spots caused by LLM instabilities on scintillator screen: (a)-(l) impact light spots of measured particles caused by LLM on scintillator screen in different time during H-and L-mode operation scenarios.

    通過EFIT 文件給出的第38140 次放電過程中在INPA 位置的具體磁場(chǎng)強(qiáng)度和分布情況可以計(jì)算出來(lái)自INPA 診斷系統(tǒng)不同測(cè)量通道(位置)的粒子在閃爍屏上的落點(diǎn)網(wǎng)格,即縱向和橫向分別表示測(cè)量粒子的能量和位置.圖9 給出了在1356 ms 時(shí)刻高速相機(jī)記錄到的閃爍體的粒子轟擊斑圖像(對(duì)應(yīng)于圖8(h))以及經(jīng)過計(jì)算后得到的能量和位置網(wǎng)格.通過對(duì)比發(fā)現(xiàn),LLM 不穩(wěn)定性被激發(fā)期間閃爍屏上的粒子轟擊斑是由INPA 診斷第3 通道測(cè)量到的.根據(jù)亮斑在水平方向的位置和縱向分布范圍可知被測(cè)量粒子主要來(lái)自于大半徑R=170.9 cm 處,其能量范圍和螺距角分別為E≈ 14—30 keV 和149.2° (v///v~ 0.86).因此,LLM 主要由來(lái)自等離子芯部的通行快離子激發(fā).考慮到INPA 診斷的系統(tǒng)誤差問題(圖5(d)),激發(fā)LLM 不穩(wěn)定性的快離子的可能位置和能量范圍分別為R≈ 170.6—171.5 cm和E≈ 12.5—32.0 keV.

    圖9 通過閃爍屏上的轟擊斑位置得到的快離子能量和位置Fig.9.Energy and birthplace of FIs based on the impact light spot on scintillator screen.

    LLM 又被稱為“非共振魚骨?!被颉帮柡汪~骨模”,也是由快粒子與內(nèi)扭曲模通過進(jìn)動(dòng)共振而激發(fā).與典型的魚骨模相比,LLM 一般發(fā)生在安全因子q分布具有反剪切或弱剪切,以及最小安全因子qmin稍大于1 的等離子體中[27,28].LLM 不僅可以引起Mirnov 磁探針的擾動(dòng)(如圖7(d)所示),還可以引起電子溫度和等離子體密度等的擾動(dòng).圖10(a)—(d)給出了具有一定空間分辨率的遠(yuǎn)紅外激光干涉儀診斷系統(tǒng)(FIR1—FIR4)分別在距離等離子體中心弦距d=3.5,10.5,17.5 和24.5 cm 處探測(cè)到的LLM 不穩(wěn)定性引起的密度擾動(dòng).LLM 在d=10.5 cm處引起的擾動(dòng)最強(qiáng),在d=3.5 cm 處較弱,而在d=17.5 和24.5 cm 處LLM 則不會(huì)引起密度擾動(dòng).該現(xiàn)象說(shuō)明LLM 不穩(wěn)定性大致局域在距離等離子體芯部r≈ 10.5 cm 處.

    圖10 通過具有空間分辨率的遠(yuǎn)紅外密度干涉儀的密度擾動(dòng)確定LLM 的局域位置Fig.10.Locations of LLM confirmed by the fluctuations in electron-density by far-infrared laser interferometer with a rough spatial resolution.

    通過計(jì)算得到的LLM 不穩(wěn)定性的頻率(fLab)與等離子體的環(huán)向旋轉(zhuǎn)頻率(ft)以及激發(fā)LLM的快離子的進(jìn)動(dòng)頻率(fp)滿足以下關(guān)系:

    其中,n表示不穩(wěn)定性的環(huán)向模數(shù),對(duì)于屬于內(nèi)扭曲模的LLM,其環(huán)向模數(shù)為1,即fLab=fp+ft.

    粒子的進(jìn)動(dòng)頻率公式[29]為

    通行粒子(κ2>1)的進(jìn)動(dòng)公式中的H(κ,s) 因子可表示為

    這里E(1/κ2)和K(1/κ2) 表示第一類和第二類不完全橢圓積分.

    通過INPA 診斷系統(tǒng),在實(shí)驗(yàn)中得到了激發(fā)LLM 不穩(wěn)定性的快離子能量范圍、螺距角和位置分別為E=12.5—32.0 keV,θ≈149.2°(v///v≈0.86)和R=170.6—171.5 cm.按照通行粒子的進(jìn)動(dòng)公式對(duì)INPA 測(cè)量到的E≈12.5—32.0 keV 快離子的進(jìn)動(dòng)頻率進(jìn)行計(jì)算,取LLM 不穩(wěn)定性局域位置r≈10.5cm處,R0=1.6m,q≈1.05,s=0.25E=12.5-32keV的通行離子的進(jìn)動(dòng)頻率fp≈2.4—6.2 kHz.另外,通過電荷交換譜儀得到了在NBI注入期間等離子體的環(huán)向旋轉(zhuǎn)頻率信息,如圖11所示.圖11(a)給出了在t=1350 ms 附近(接近實(shí)驗(yàn)分析的t=1356 ms 時(shí)刻)等離子體的環(huán)向旋轉(zhuǎn)頻率的空間分布.在LLM 不穩(wěn)定性局域位置處r≈ 10.5 cm (R=170.9 cm)的等離子環(huán)向旋轉(zhuǎn)頻率ft≈8.1 kHz,如圖中虛線所對(duì)應(yīng)的位置.圖11(b)給出了在LLM 不穩(wěn)定性局域位置附近的等離子體環(huán)向旋轉(zhuǎn)頻率隨時(shí)間的演化.ft的變化趨勢(shì)與LLM 頻率變化大體一致,因此,LLM不穩(wěn)定性頻率的變化主要是由ft的變化引起的.

    圖11 (a) 1356 ms 前后等離子體轉(zhuǎn)動(dòng)頻率隨大半徑的變化;(b)大半徑在170 cm 附近,等離子體轉(zhuǎn)動(dòng)頻率隨時(shí)間的變化Fig.11.(a) Variations of plasma rotation frequency with R around 1356 ms;(b) variations of plasma rotation frequency with time for R around 170 cm.

    將上述計(jì)算得到的能量E=12.5—32.0 keV的通行離子的進(jìn)動(dòng)頻率(fp≈2.4—6.2 kHz)和實(shí)驗(yàn)得到的在LLM 不穩(wěn)定性局域位置等離子體環(huán)向旋轉(zhuǎn)頻率(ft=8.1 kHz)代入(1)式可得,fLab范圍在10.5—14.3 kHz,見表3.該結(jié)果與頻譜圖中觀測(cè)到在t=1356 ms 時(shí)刻的LLM 不穩(wěn)定性的頻率fLLM=13.4 kHz 接近.因此,在誤差范圍內(nèi)可以證明實(shí)驗(yàn)中的LLM 不穩(wěn)定性是由能量E=12.5—32.0 keV、螺距角θ≈149.2° (v///v≈ 0.86)和局域位置在R=170.6—171.5 cm 處的通行快離子激發(fā)的.

    表3 通行快離子的理論計(jì)算頻率值與實(shí)驗(yàn)觀測(cè)LLM 不穩(wěn)定性頻率對(duì)比Table 3.Comparisons between the calculated frequency of EIs and fLLM.

    4 總結(jié)與討論

    根據(jù)HL-2A 裝置的運(yùn)行參數(shù)和實(shí)際空間位置等條件,設(shè)計(jì)了一套INPA 診斷系統(tǒng)用于高能量粒子物理研究中快離子的空間分布、能量和螺距角等關(guān)鍵信息的測(cè)量,以及高能量粒子不穩(wěn)定性與快離子相互作用的研究.在HL-2A 裝置物理實(shí)驗(yàn)中,使用INPA 診斷系統(tǒng)對(duì)快離子進(jìn)行了觀測(cè),得到了LLM 不穩(wěn)定性主要由能量、螺距角和位置分別為E=12.5—32.0 keV,θ ≈149.2° (v///v≈ 0.86)和R=170.1—171.5 cm 處的芯部快離子激發(fā),并通過理論計(jì)算對(duì)該結(jié)果進(jìn)行了驗(yàn)證.由于INPA 診斷系統(tǒng)的安裝位置距離NBI 束線較遠(yuǎn),NBI 粒子束在注入到等離子體中所觀測(cè)位置后的衰減比較嚴(yán)重,所能提供給背景等離子體中快離子進(jìn)行電荷交換的中性粒子很少.另外,采用對(duì)光傳輸效率較差的成像光纖束對(duì)閃爍屏上的光斑圖像進(jìn)行傳輸,這些都造成了INPA 系統(tǒng)的有效信號(hào)較弱.目前,HL-2A 裝置上的INPA 系統(tǒng)主要適用于在NBI條件下,對(duì)快粒子約束條件較差的L 模運(yùn)行模式的低密度等離子體進(jìn)行測(cè)量.在此條件下,INPA診斷系統(tǒng)的時(shí)間分別率約為37 ms.在下一輪的實(shí)驗(yàn)中將設(shè)計(jì)專門的圖像傳輸光路提高INPA 診斷系統(tǒng)的信噪比和時(shí)間分辨率.INPA系統(tǒng)在HL-2A裝置上的設(shè)計(jì)和實(shí)驗(yàn)為該診斷在HL-2M 裝置上的發(fā)展積累了重要的物理設(shè)計(jì)、工程安裝和實(shí)驗(yàn)運(yùn)行經(jīng)驗(yàn).根據(jù)HL-2M 裝置高能量粒子物理實(shí)驗(yàn)需求,INPA 診斷系統(tǒng)已經(jīng)開始在HL-2M 裝置上進(jìn)行設(shè)計(jì)和加工.

    猜你喜歡
    診斷系統(tǒng)螺距膜片
    變螺距螺紋數(shù)控車削技巧
    橡膠膜片耐液體性能及等級(jí)評(píng)定
    基于Workbench的膜片聯(lián)軸節(jié)變形分析方法
    區(qū)間軌道電路智能診斷系統(tǒng)的探討
    設(shè)備在線診斷系統(tǒng)在唐鋼的建設(shè)與應(yīng)用
    常用粗牙螺紋螺距簡(jiǎn)易計(jì)算公式及應(yīng)用
    厚度梯度對(duì)鈦制橢球形膜片翻轉(zhuǎn)性能的影響分析
    利用體塊PZT制備膜片式壓電微泵
    雙源CT大螺距掃描在胸部篩查的應(yīng)用
    連鑄板坯質(zhì)量在線診斷系統(tǒng)的應(yīng)用
    新疆鋼鐵(2015年2期)2015-11-07 03:27:52
    叶爱在线成人免费视频播放| 国产精品综合久久久久久久免费 | 久久99一区二区三区| 久久人妻福利社区极品人妻图片| 国产精品二区激情视频| 国产97色在线日韩免费| 免费高清视频大片| av天堂在线播放| 精品一区二区三卡| 久久中文看片网| 午夜亚洲福利在线播放| 国产熟女午夜一区二区三区| 国产熟女午夜一区二区三区| 久久久国产成人精品二区 | 老汉色∧v一级毛片| 老汉色∧v一级毛片| 自线自在国产av| 黑人巨大精品欧美一区二区蜜桃| 黑人巨大精品欧美一区二区蜜桃| 久久香蕉精品热| 变态另类成人亚洲欧美熟女 | 九色亚洲精品在线播放| 三上悠亚av全集在线观看| 黄色女人牲交| 精品久久久久久久毛片微露脸| 一级,二级,三级黄色视频| 午夜老司机福利片| 亚洲免费av在线视频| av网站在线播放免费| 黄色怎么调成土黄色| 欧美激情久久久久久爽电影 | 搡老岳熟女国产| 伊人久久大香线蕉亚洲五| 在线观看一区二区三区激情| 欧美激情 高清一区二区三区| 男女床上黄色一级片免费看| 一区福利在线观看| 精品国产美女av久久久久小说| 一二三四社区在线视频社区8| 国产精品二区激情视频| 精品一品国产午夜福利视频| 波多野结衣av一区二区av| 午夜精品国产一区二区电影| 天天躁夜夜躁狠狠躁躁| 亚洲精品中文字幕在线视频| 999久久久精品免费观看国产| 亚洲自拍偷在线| 国产精品免费一区二区三区在线| 亚洲性夜色夜夜综合| 免费人成视频x8x8入口观看| www.www免费av| 桃红色精品国产亚洲av| 日本wwww免费看| 日本撒尿小便嘘嘘汇集6| 老司机午夜福利在线观看视频| 精品日产1卡2卡| 99久久精品国产亚洲精品| 久久久久久久精品吃奶| 亚洲欧美一区二区三区黑人| 满18在线观看网站| 老鸭窝网址在线观看| 男女下面插进去视频免费观看| 中文字幕另类日韩欧美亚洲嫩草| 91麻豆精品激情在线观看国产 | 久久人妻福利社区极品人妻图片| 女生性感内裤真人,穿戴方法视频| 精品一品国产午夜福利视频| 在线观看一区二区三区| 日韩欧美一区二区三区在线观看| 亚洲自拍偷在线| 久久人人精品亚洲av| 午夜影院日韩av| 亚洲精品一区av在线观看| 国产深夜福利视频在线观看| 亚洲少妇的诱惑av| 午夜影院日韩av| 久久国产精品男人的天堂亚洲| 色在线成人网| 亚洲美女黄片视频| 99re在线观看精品视频| 久久天堂一区二区三区四区| a级毛片在线看网站| 多毛熟女@视频| 亚洲五月色婷婷综合| 国产91精品成人一区二区三区| 国产精品成人在线| 亚洲av成人一区二区三| 美女大奶头视频| 久久精品国产亚洲av高清一级| 亚洲一区高清亚洲精品| 97人妻天天添夜夜摸| 国产在线观看jvid| 69av精品久久久久久| 亚洲国产中文字幕在线视频| 麻豆成人av在线观看| 美女高潮喷水抽搐中文字幕| 99久久国产精品久久久| 日韩国内少妇激情av| 亚洲片人在线观看| 真人一进一出gif抽搐免费| 脱女人内裤的视频| 伊人久久大香线蕉亚洲五| 国产黄a三级三级三级人| 一进一出抽搐动态| 一本综合久久免费| 精品电影一区二区在线| 少妇裸体淫交视频免费看高清 | 国产亚洲精品一区二区www| 免费看a级黄色片| 国产伦人伦偷精品视频| 夜夜夜夜夜久久久久| 91在线观看av| 在线观看一区二区三区激情| 久久久久九九精品影院| 亚洲在线自拍视频| 国产野战对白在线观看| 黄频高清免费视频| 国产精品久久电影中文字幕| 老熟妇仑乱视频hdxx| 淫秽高清视频在线观看| 精品欧美一区二区三区在线| 欧美激情久久久久久爽电影 | 免费看十八禁软件| 在线播放国产精品三级| 午夜视频精品福利| 国产av精品麻豆| 亚洲成国产人片在线观看| 无限看片的www在线观看| 免费少妇av软件| 人成视频在线观看免费观看| 50天的宝宝边吃奶边哭怎么回事| 黑人巨大精品欧美一区二区蜜桃| 国产伦人伦偷精品视频| 午夜久久久在线观看| 亚洲 欧美 日韩 在线 免费| 国产成人精品在线电影| 国产成人啪精品午夜网站| 九色亚洲精品在线播放| 国产精品一区二区三区四区久久 | 国产免费现黄频在线看| 熟女少妇亚洲综合色aaa.| 9热在线视频观看99| 久热这里只有精品99| 国产精品偷伦视频观看了| 亚洲 欧美一区二区三区| 男女做爰动态图高潮gif福利片 | 精品国产一区二区久久| 老司机在亚洲福利影院| 黄片小视频在线播放| 欧美成狂野欧美在线观看| 看免费av毛片| 亚洲少妇的诱惑av| 十八禁人妻一区二区| 精品久久久久久久久久免费视频 | 99国产精品一区二区蜜桃av| 天堂影院成人在线观看| 久久久久亚洲av毛片大全| 亚洲欧美一区二区三区黑人| 亚洲在线自拍视频| 新久久久久国产一级毛片| 日韩av在线大香蕉| 久久热在线av| 亚洲国产欧美日韩在线播放| 在线看a的网站| 成人影院久久| 99国产极品粉嫩在线观看| 狠狠狠狠99中文字幕| 搡老乐熟女国产| 日本免费一区二区三区高清不卡 | 色综合欧美亚洲国产小说| 在线观看66精品国产| 老汉色av国产亚洲站长工具| 亚洲欧美激情综合另类| 亚洲 欧美 日韩 在线 免费| 亚洲七黄色美女视频| 国产成人啪精品午夜网站| 丰满饥渴人妻一区二区三| 成人亚洲精品一区在线观看| 欧美日韩亚洲国产一区二区在线观看| 欧美日韩黄片免| 成人免费观看视频高清| 高清在线国产一区| 91成人精品电影| 十八禁人妻一区二区| 日日摸夜夜添夜夜添小说| 色在线成人网| 久久精品人人爽人人爽视色| 午夜精品在线福利| 成熟少妇高潮喷水视频| 可以在线观看毛片的网站| 国产成+人综合+亚洲专区| 国产欧美日韩一区二区三区在线| 午夜福利欧美成人| 香蕉丝袜av| 在线观看免费午夜福利视频| 伦理电影免费视频| 岛国在线观看网站| 最好的美女福利视频网| 国产伦一二天堂av在线观看| 三上悠亚av全集在线观看| 99精品欧美一区二区三区四区| 国产精品 国内视频| 黄色a级毛片大全视频| www.自偷自拍.com| 少妇裸体淫交视频免费看高清 | 搡老岳熟女国产| 国产精品久久视频播放| 午夜日韩欧美国产| 久久天堂一区二区三区四区| 久久久久久亚洲精品国产蜜桃av| 女同久久另类99精品国产91| 一个人免费在线观看的高清视频| 亚洲欧美日韩高清在线视频| 国产高清激情床上av| 波多野结衣一区麻豆| 亚洲情色 制服丝袜| 久久欧美精品欧美久久欧美| 99riav亚洲国产免费| 人人妻人人添人人爽欧美一区卜| 搡老熟女国产l中国老女人| 国产一区二区三区视频了| 一a级毛片在线观看| 露出奶头的视频| 欧美人与性动交α欧美精品济南到| 在线av久久热| 一级毛片精品| 日本黄色视频三级网站网址| 午夜福利,免费看| 久久精品成人免费网站| 黄色片一级片一级黄色片| 成人18禁在线播放| 免费日韩欧美在线观看| 波多野结衣av一区二区av| 久久久久久久午夜电影 | 动漫黄色视频在线观看| 极品教师在线免费播放| 欧美最黄视频在线播放免费 | 性少妇av在线| 岛国在线观看网站| 欧美日韩中文字幕国产精品一区二区三区 | 深夜精品福利| 成人国语在线视频| 大型黄色视频在线免费观看| 色精品久久人妻99蜜桃| 好看av亚洲va欧美ⅴa在| 欧美老熟妇乱子伦牲交| 亚洲一区高清亚洲精品| 波多野结衣高清无吗| 这个男人来自地球电影免费观看| 亚洲av成人一区二区三| 男女床上黄色一级片免费看| 中文字幕色久视频| 国产高清视频在线播放一区| 国产精品成人在线| 在线观看免费视频日本深夜| 亚洲一码二码三码区别大吗| 脱女人内裤的视频| 欧美午夜高清在线| 欧美另类亚洲清纯唯美| 波多野结衣av一区二区av| 欧美乱妇无乱码| 欧美激情久久久久久爽电影 | av片东京热男人的天堂| 每晚都被弄得嗷嗷叫到高潮| 日韩av在线大香蕉| av电影中文网址| 99国产精品免费福利视频| 亚洲avbb在线观看| 免费看a级黄色片| 后天国语完整版免费观看| 国产成+人综合+亚洲专区| 久久久水蜜桃国产精品网| 国内久久婷婷六月综合欲色啪| 亚洲美女黄片视频| 极品人妻少妇av视频| 欧美日韩黄片免| 亚洲七黄色美女视频| 啦啦啦在线免费观看视频4| av在线播放免费不卡| 久久这里只有精品19| 中文亚洲av片在线观看爽| 亚洲av成人一区二区三| 男女床上黄色一级片免费看| 男女高潮啪啪啪动态图| 精品国产一区二区三区四区第35| av福利片在线| 12—13女人毛片做爰片一| 亚洲五月天丁香| 搡老熟女国产l中国老女人| 欧美激情极品国产一区二区三区| av网站在线播放免费| 精品电影一区二区在线| 久久久国产成人免费| 久久香蕉精品热| 夜夜爽天天搞| 国产精品乱码一区二三区的特点 | 国产区一区二久久| 男人舔女人的私密视频| 91国产中文字幕| 一级a爱视频在线免费观看| 嫩草影院精品99| 国产精品一区二区在线不卡| 午夜福利一区二区在线看| 精品久久久久久电影网| 国产精品综合久久久久久久免费 | 操美女的视频在线观看| 身体一侧抽搐| 两性夫妻黄色片| 欧美日韩av久久| 天堂中文最新版在线下载| 99国产精品99久久久久| 长腿黑丝高跟| av国产精品久久久久影院| 亚洲视频免费观看视频| 久久久国产成人精品二区 | 法律面前人人平等表现在哪些方面| 亚洲少妇的诱惑av| 国产亚洲欧美98| 黑人巨大精品欧美一区二区蜜桃| 中文字幕人妻熟女乱码| 在线观看66精品国产| av在线天堂中文字幕 | 国产精品免费一区二区三区在线| 久久久久久大精品| 欧美一级毛片孕妇| 亚洲一区中文字幕在线| 色在线成人网| 大陆偷拍与自拍| 国产精品永久免费网站| 亚洲国产欧美日韩在线播放| 免费女性裸体啪啪无遮挡网站| 日韩一卡2卡3卡4卡2021年| 日日摸夜夜添夜夜添小说| 亚洲 国产 在线| 欧美激情极品国产一区二区三区| 亚洲伊人色综图| 亚洲性夜色夜夜综合| 国产精品 欧美亚洲| av免费在线观看网站| 欧美老熟妇乱子伦牲交| 午夜精品久久久久久毛片777| 啦啦啦在线免费观看视频4| 亚洲成人久久性| 免费女性裸体啪啪无遮挡网站| 美女扒开内裤让男人捅视频| 国产视频一区二区在线看| 亚洲视频免费观看视频| 欧美精品一区二区免费开放| 国产99久久九九免费精品| 欧美+亚洲+日韩+国产| 黄频高清免费视频| 夫妻午夜视频| 黄色 视频免费看| 久久久精品欧美日韩精品| 色尼玛亚洲综合影院| 欧美日韩亚洲高清精品| 国产精华一区二区三区| 麻豆成人av在线观看| 国产精品久久视频播放| 一区二区三区激情视频| 最近最新免费中文字幕在线| 国产亚洲精品久久久久久毛片| 久久亚洲真实| 黄色成人免费大全| 日韩高清综合在线| 身体一侧抽搐| 久99久视频精品免费| 99久久久亚洲精品蜜臀av| 制服诱惑二区| 男人舔女人的私密视频| 在线观看日韩欧美| 99久久人妻综合| 精品国产一区二区久久| 久久久国产欧美日韩av| av网站在线播放免费| 在线视频色国产色| 高清av免费在线| 国产主播在线观看一区二区| 最新美女视频免费是黄的| 亚洲av美国av| 久久久久久久久免费视频了| 亚洲激情在线av| 丰满迷人的少妇在线观看| 亚洲av美国av| 看免费av毛片| 五月开心婷婷网| 亚洲国产欧美网| 露出奶头的视频| 黑人操中国人逼视频| 日韩视频一区二区在线观看| 精品久久久精品久久久| 69av精品久久久久久| 欧美精品亚洲一区二区| a在线观看视频网站| 精品国产亚洲在线| 别揉我奶头~嗯~啊~动态视频| 免费在线观看影片大全网站| 国产国语露脸激情在线看| 老汉色∧v一级毛片| 午夜福利在线观看吧| 不卡av一区二区三区| 99国产综合亚洲精品| 香蕉久久夜色| av网站免费在线观看视频| 麻豆成人av在线观看| 青草久久国产| 亚洲国产欧美日韩在线播放| aaaaa片日本免费| 亚洲成人精品中文字幕电影 | 动漫黄色视频在线观看| 99热国产这里只有精品6| 成人国产一区最新在线观看| 欧美中文日本在线观看视频| av在线播放免费不卡| 一级片'在线观看视频| 亚洲片人在线观看| 91老司机精品| 日本黄色日本黄色录像| 国产成人精品久久二区二区免费| av在线天堂中文字幕 | 国产精品爽爽va在线观看网站 | 国产亚洲欧美精品永久| 1024香蕉在线观看| 国产成人av激情在线播放| 国产成年人精品一区二区 | 狠狠狠狠99中文字幕| 水蜜桃什么品种好| 国产av又大| 天天躁夜夜躁狠狠躁躁| 桃色一区二区三区在线观看| 精品一品国产午夜福利视频| 久久人人爽av亚洲精品天堂| 91成年电影在线观看| 亚洲自偷自拍图片 自拍| 91av网站免费观看| 亚洲专区国产一区二区| 黄片大片在线免费观看| 伊人久久大香线蕉亚洲五| 精品一区二区三区四区五区乱码| 在线视频色国产色| 日本vs欧美在线观看视频| 女人被躁到高潮嗷嗷叫费观| av在线天堂中文字幕 | 精品国产超薄肉色丝袜足j| 天堂中文最新版在线下载| 黄色怎么调成土黄色| 午夜免费观看网址| 精品久久久久久成人av| 日韩精品免费视频一区二区三区| 日本vs欧美在线观看视频| 老熟妇仑乱视频hdxx| 午夜精品在线福利| 欧美日韩乱码在线| 日日爽夜夜爽网站| 国产精品国产av在线观看| 首页视频小说图片口味搜索| 久久婷婷成人综合色麻豆| av超薄肉色丝袜交足视频| 亚洲va日本ⅴa欧美va伊人久久| 免费高清在线观看日韩| 欧美日韩精品网址| 桃红色精品国产亚洲av| 亚洲精品国产区一区二| 精品国产一区二区久久| 男女下面进入的视频免费午夜 | 国产成人影院久久av| 国产无遮挡羞羞视频在线观看| 操出白浆在线播放| 亚洲成人免费av在线播放| 老司机午夜福利在线观看视频| 国产精品久久久人人做人人爽| 欧美大码av| 日本黄色日本黄色录像| 高清欧美精品videossex| 亚洲国产欧美一区二区综合| 丰满饥渴人妻一区二区三| 热99re8久久精品国产| 国产xxxxx性猛交| 国产片内射在线| 日本黄色日本黄色录像| 欧美 亚洲 国产 日韩一| 日韩高清综合在线| 久久久久亚洲av毛片大全| 国产亚洲欧美精品永久| 亚洲精品美女久久久久99蜜臀| 国产熟女xx| 夜夜爽天天搞| 18禁观看日本| 亚洲男人的天堂狠狠| www.999成人在线观看| 一级毛片高清免费大全| 欧美在线一区亚洲| 亚洲三区欧美一区| 久99久视频精品免费| 免费少妇av软件| 黑人欧美特级aaaaaa片| 99精国产麻豆久久婷婷| 亚洲专区国产一区二区| 久久香蕉激情| 亚洲国产精品999在线| 国产日韩一区二区三区精品不卡| 亚洲av日韩精品久久久久久密| 又黄又爽又免费观看的视频| 欧美乱妇无乱码| 村上凉子中文字幕在线| 两个人看的免费小视频| 五月开心婷婷网| a级毛片黄视频| 99在线视频只有这里精品首页| 成人三级黄色视频| 91老司机精品| 精品久久久久久久毛片微露脸| 欧美精品亚洲一区二区| 99久久国产精品久久久| 久热这里只有精品99| 免费一级毛片在线播放高清视频 | 国产精品 欧美亚洲| 国产成年人精品一区二区 | 日韩大码丰满熟妇| 成人av一区二区三区在线看| 久久草成人影院| 99国产精品免费福利视频| 国产精品久久久久成人av| 成人手机av| 久久久国产欧美日韩av| 精品久久久久久久久久免费视频 | 看片在线看免费视频| 交换朋友夫妻互换小说| 成人av一区二区三区在线看| 高清av免费在线| 久久香蕉国产精品| 九色亚洲精品在线播放| 精品久久久久久电影网| 宅男免费午夜| 看黄色毛片网站| 久久精品成人免费网站| www.自偷自拍.com| 亚洲男人的天堂狠狠| 国产一区二区激情短视频| 欧美日韩精品网址| 别揉我奶头~嗯~啊~动态视频| 国产不卡一卡二| 欧美最黄视频在线播放免费 | 精品久久久精品久久久| 久久性视频一级片| 午夜影院日韩av| 亚洲情色 制服丝袜| 免费女性裸体啪啪无遮挡网站| 午夜两性在线视频| 99久久精品国产亚洲精品| 国产免费av片在线观看野外av| 国产亚洲精品综合一区在线观看 | 黑人巨大精品欧美一区二区蜜桃| 免费久久久久久久精品成人欧美视频| 天堂√8在线中文| 欧美性长视频在线观看| 一区二区三区精品91| 亚洲国产精品999在线| xxxhd国产人妻xxx| av网站免费在线观看视频| 韩国精品一区二区三区| 亚洲男人的天堂狠狠| 久久精品人人爽人人爽视色| 国产视频一区二区在线看| 很黄的视频免费| 久久精品国产清高在天天线| 国产免费男女视频| 麻豆国产av国片精品| 天天影视国产精品| av国产精品久久久久影院| 久久亚洲真实| 免费看a级黄色片| 50天的宝宝边吃奶边哭怎么回事| 免费观看人在逋| 麻豆成人av在线观看| 亚洲人成电影免费在线| 日本五十路高清| 后天国语完整版免费观看| 国产区一区二久久| 一级a爱片免费观看的视频| 欧美成人免费av一区二区三区| 亚洲视频免费观看视频| 亚洲欧洲精品一区二区精品久久久| 女人高潮潮喷娇喘18禁视频| 国产极品粉嫩免费观看在线| 欧美日韩黄片免| 中文字幕人妻熟女乱码| 好男人电影高清在线观看| bbb黄色大片| 国产高清视频在线播放一区| 久久久精品国产亚洲av高清涩受| 免费在线观看影片大全网站| 亚洲人成电影免费在线| 亚洲av五月六月丁香网| 嫁个100分男人电影在线观看| 少妇粗大呻吟视频| 极品人妻少妇av视频| 天天影视国产精品| 国产aⅴ精品一区二区三区波| 黄片小视频在线播放| 亚洲精品久久午夜乱码| 久久久久久久午夜电影 | 精品国产一区二区三区四区第35| 校园春色视频在线观看| 亚洲av美国av| 亚洲精品av麻豆狂野| 少妇裸体淫交视频免费看高清 | 不卡一级毛片| aaaaa片日本免费| 亚洲片人在线观看| 久久青草综合色| 久久人妻熟女aⅴ| 亚洲成人免费电影在线观看| 久久这里只有精品19| 亚洲欧美日韩无卡精品|