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

    EAST 上離子回旋波與中性束注入?yún)f(xié)同加熱產(chǎn)生的高能粒子分布及輸運(yùn)研究*

    2023-11-24 05:05:42張偉張新軍劉魯南朱光輝楊樺張華朋鄭藝峰何開洋黃娟
    物理學(xué)報(bào) 2023年21期
    關(guān)鍵詞:限制器高能共振

    張偉 張新軍? 劉魯南 朱光輝 楊樺 張華朋 鄭藝峰 何開洋 黃娟

    1) (中國(guó)科學(xué)院合肥物質(zhì)科學(xué)研究院,等離子體物理研究所,合肥 230031)

    2) (深圳大學(xué)物理與光電工程學(xué)院,深圳 518060)

    3) (中國(guó)科學(xué)技術(shù)大學(xué)核科學(xué)技術(shù)學(xué)院,合肥 230026)

    在磁約束聚變等離子體中,離子回旋共振加熱(ICRF)與中性束注入(NBI)是兩種主要的加熱方法.它們的協(xié)同加熱一直都是聚變領(lǐng)域研究的重點(diǎn).本文首先闡明了ICRF 高次諧波加熱以及ICRF 與NBI 協(xié)同加熱的基本原理.通過(guò)EAST 托卡馬克上實(shí)驗(yàn)和相應(yīng)的TRANSP 模擬,發(fā)現(xiàn)了ICRF 與NBI 的協(xié)同加熱不僅可以顯著提高等離子體參數(shù)(極向比壓、等離子體儲(chǔ)能、離子溫度、中子產(chǎn)額等),而且能產(chǎn)生大量高能粒子,形成高能粒子尾巴.例如,1 MW 的ICRF 三次諧波可將初始能量為60 keV 的NBI 高能氘離子加速至600 keV.通過(guò)改變氫少子含量、提高ICRF 和NBI 加熱功率、使用ICRF 在軸加熱、優(yōu)化NBI 注入角度等,可以有效地提高協(xié)同加熱效率以及高能粒子的能量.進(jìn)一步地,將協(xié)同加熱產(chǎn)生的高能粒子分布代入粒子軌道程序中,計(jì)算了高能粒子的輸運(yùn)以及其在第一壁上的損失.結(jié)果表明,損失的高能粒子的初始位置位于低場(chǎng)側(cè),且損失軌道大部分為捕獲粒子軌道.高能粒子損失位置主要位于主限制器以及ICRF 和低雜波限制器的中上平面.這些損失的高能粒子被認(rèn)為是造成限制器上熱斑的主要原因之一.

    1 引言

    自20 世紀(jì)90 年代以來(lái),離子回旋共振加熱(ICRF)與中性束注入(NBI)的協(xié)同加熱一直都是磁約束聚變研究的重點(diǎn).在世界著名的磁約束聚變裝置如JET[1,2],ASDEX Upgrade[3,4],DIII-D[5],TEXTOR[6],JT-60[7-9]上,都開展了大量的協(xié)同加熱的實(shí)驗(yàn)和理論研究.比如,日本JT-60 托卡馬克裝置上的實(shí)驗(yàn)表明ICRF 和NBI 的協(xié)同加熱可以產(chǎn)生高能離子尾巴,并將等離子體能量約束時(shí)間提高到了只有NBI 或ICRF 單獨(dú)加熱時(shí)的3倍[7].歐盟JET 托卡馬克上的實(shí)驗(yàn)表明NBI 的高能氘離子與ICRF 二次[2,10]或三次[11,12]諧波的協(xié)同可以產(chǎn)生能量達(dá)2 MeV 的高能氘離子.最近的JET 實(shí)驗(yàn)發(fā)現(xiàn)可以將NBI 產(chǎn)生的氘高能粒子視為第三種離子,通過(guò)ICRF 將其在D-3He[13,14]或H-D[1,15]等離子體中的離子-離子混雜共振處進(jìn)一步加速,產(chǎn)生協(xié)同效應(yīng).德國(guó)ASDEX Upgrade 托卡馬克上的實(shí)驗(yàn)和模擬表明NBI 與ICRF 的協(xié)同加熱可以將NBI約50 keV 的高能離子能量提高至500 keV[4,16].

    ICRF 與NBI 的協(xié)同加熱不僅可以大幅提高注入的NBI 高能粒子能量,也能提高等離子體的參數(shù),包括離子溫度、儲(chǔ)能、極向比壓、中子產(chǎn)額等.在JET 的DT 實(shí)驗(yàn)中,通過(guò)使用5.5 MW ICRF與23.5 MW NBI 獲得了36 keV 的離子溫度,這相對(duì)于單純使用29 MW NBI 獲得的離子溫度提高了10 keV[17].EAST 托卡馬克上的實(shí)驗(yàn)也表明,通過(guò)ICRF 和NBI 的協(xié)同加熱,可以大幅提高等離子體的極向比壓、等離子體儲(chǔ)能、離子溫度和中子產(chǎn)額[18].特別是在歐盟JET 和美國(guó)TFTR 裝置上,都曾廣泛使用了ICRF 與NBI 的協(xié)同加熱來(lái)提高等離子體的性能,并分別產(chǎn)生了15 MW 和10 MW 的聚變輸出功率[19,20].

    本文重點(diǎn)分析EAST 托卡馬克上ICRF 與NBI協(xié)同加熱產(chǎn)生的高能粒子分布及輸運(yùn).通過(guò)協(xié)同加熱產(chǎn)生的大量高能粒子,可以有效地提高等離子體性能和聚變反應(yīng)率.然而,未經(jīng)充分慢化的高能粒子損失則可能對(duì)裝置器壁造成損傷.因此,理解這些高能粒子的行為對(duì)提高未來(lái)聚變堆的聚變輸出功率、實(shí)現(xiàn)高參數(shù)等離子體的安全運(yùn)行尤為重要.

    2 協(xié)同加熱基本原理

    由于有限拉莫爾半徑效應(yīng),ICRF 的電場(chǎng)可有效加速位于高次諧波共振層附近的NBI 高能離子,產(chǎn)生協(xié)同效應(yīng).協(xié)同加熱利用ICRF 快波電場(chǎng)梯度來(lái)加速NBI 高能離子,其基本物理圖像如圖1 所示.當(dāng)高能離子能量越大(即回旋半徑越大)時(shí),其經(jīng)歷的電場(chǎng)梯度越大,從ICRF 電場(chǎng)中獲得的能量也越大,進(jìn)而更容易被加速.當(dāng)回旋軌道直徑等于半波長(zhǎng)時(shí),高能粒子可以從ICRF 電場(chǎng)中獲得最大的能量.通過(guò)計(jì)算可知,1 MeV 的氘離子其回旋軌道直徑(7.2 cm)仍遠(yuǎn)小于半波長(zhǎng)(dm 到m 量級(jí)).而目前聚變裝置中能產(chǎn)生的高能粒子最大能量也在MeV 量級(jí).因此,一般認(rèn)為能量越大的離子越容易被ICRF 高次諧波加熱.

    圖1 (a) 捕獲離子在離子回旋共振加熱層被加速示意圖;(b) 高次諧波加熱基本原理圖Fig.1.(a) Demonstration of trapped ion acceleration at the ion cyclotron resonance layer;(b) demonstration of basic mechanism of high harmonic ICRF heating.

    ICRF 與NBI 協(xié)同加熱產(chǎn)生的高能離子分布(fi)可以用Fokker-Planck 方程來(lái)計(jì)算[21]:

    其中C(fi),S(v)和L(v) 分別為碰撞項(xiàng)、源項(xiàng)和損失項(xiàng).兩種加熱的協(xié)同作用則是通過(guò)波與粒子相互作用項(xiàng)來(lái)實(shí)現(xiàn).這里v為離子速度;DRF是波準(zhǔn)線性算法,與ICRF 的電場(chǎng)、垂直波矢、能量密度以及共振離子的回旋頻率和垂直速度相關(guān).本文中使用的TRANSP 集成模擬,通過(guò)TORIC(ICRF 程序)、NUBEAM(NBI 程序)和FPP(計(jì)算Fokker-Planck 方程模塊)的自洽迭代,可以很好地計(jì)算協(xié)同加熱產(chǎn)生的高能離子分布.

    ICRF 與NBI 的協(xié)同加熱一般要求選定的高次諧波共振層位于磁軸附近,這樣才能使NBI 高能粒子盡可能多地位于共振層的多普勒展寬內(nèi),并獲得加速.ICRF 波與粒子共振關(guān)系為

    其中ω是ICRF 天線加熱頻率,ωcj是共振粒子回旋頻率,k//和v//是波數(shù)和平行速度,n則代表了諧波數(shù).對(duì)于D 離子,常用的高次諧波加熱為二次(n=2)和三次(n=3)諧波加熱.從波與粒子共振關(guān)系,可得諧波共振的位置為

    其中B0和R0分別為磁軸處的磁場(chǎng)和大半徑,qi和mi分別為帶電粒子的電荷和質(zhì)量.因此,在已知天線加熱頻率的前提下,離子的回旋共振層可以通過(guò)環(huán)向磁場(chǎng)來(lái)確定.對(duì)于NBI 高能氘離子來(lái)說(shuō),還會(huì)存在較大的多普勒展寬,位于其內(nèi)的高能離子都能獲得有效加速.其中,為背景離子的速度,ωci為離子回旋頻率.

    基于(3)式,可得出EAST 的氘離子諧波共振層位置隨磁場(chǎng)的變化,如圖2(a)所示.為使諧波共振層位于磁軸(R0=1.9 m)處,二次到五次諧波加熱所需的環(huán)向磁場(chǎng)分別為Bt=2.43,1.62,1.21 0.97 T.EAST 常用的磁場(chǎng)為2.4—2.5 T,因而氘離子二次諧波是最常用的高次諧波加熱模式.然,種加熱模式是相互競(jìng)爭(zhēng)的.一般說(shuō)來(lái),當(dāng)氫少子的含量大于1%時(shí),氫少子加熱模式占主導(dǎo),ICRF電場(chǎng)主要加速氫離子;當(dāng)氫少子含量少于1%時(shí),氘離子的二次諧波加熱模式占主導(dǎo),ICRF 電場(chǎng)主要加速背景和注入的氘離子.氫高能離子相對(duì)氘高能離子在慢化過(guò)程中更易加熱電子.

    圖2 (a) 高次諧波對(duì)應(yīng)的環(huán)向磁場(chǎng).其中n 來(lái)自方程(2),代表第n 次諧波加熱,虛線為EAST 磁軸R0=1.92 m 處;(b) 當(dāng)環(huán)向磁場(chǎng)Bt=1.63 T時(shí)n=2—4 次諧波共振所在位置;(c) EAST 上四條NBI 束線和兩個(gè)ICRF 天線的環(huán)向位置Fig.2.(a) Toroidal magnetic field for high harmonic ICRF heating.The integer n comes from Eq.(2) and represents the nth high harmonic heating.The vertical dashed line represents the magnetic axis at R0=1.92m;(b) high harmonic resonance positions for n=2-4 with Bt=1.63 T.(c) toroidal locations of the four NBI linesand two ICRF antennas.

    當(dāng)磁場(chǎng)降低至1.6 T 時(shí),此時(shí)氘離子三次諧波位于磁軸,而二次和四次諧波共振層位于裝置外,如圖2(b).因而不管主等離子中氫離子含量如何,氘離子的三次諧波一直占主導(dǎo),這樣可以使三次諧波加熱最大化.但是在EAST 上,使用低縱場(chǎng)容易而,氘離子的二次諧波與H 離子的基頻(n=1)共振位置是相同的,在存在氫少子的情況下,這兩導(dǎo)致等離子體的約束變差,放電不穩(wěn)定,因此不經(jīng)常使用.值得一提的是,在ASDEX-U 等裝置中還嘗試過(guò)改變天線頻率和環(huán)向磁場(chǎng),使二次諧波和三次諧波共振層同時(shí)位于裝置內(nèi),實(shí)現(xiàn)了多種諧波同時(shí)與NBI 協(xié)同加熱.

    EAST 目前共有兩個(gè)ICRF 天線,最大耦合功率共為4 MW.配有四條NBI 束線,源功率共8 MW,注入等離子體的最大功率共為4—5 MW.這些加熱系統(tǒng)的環(huán)向位置如圖2(c)所示.此外,EAST 還有多項(xiàng)高能粒子診斷,包括中子發(fā)射譜儀(NES)、中子飛行時(shí)間譜儀(TOFED)、快離子Dα光譜儀(FIDA)和高能粒子損失探針(FILD).因此,EAST 具備了開展ICRF 與NBI 協(xié)同加熱的實(shí)驗(yàn)條件.

    3 協(xié)同加熱產(chǎn)生的高能粒子分布

    在EAST 上,開展了大量的ICRF 二次或三次諧波與NBI 協(xié)同加熱的實(shí)驗(yàn),重點(diǎn)研究了協(xié)同加熱產(chǎn)生的高能粒子分布和輸運(yùn),并取得了一系列成果.研究發(fā)現(xiàn),ICRF 電場(chǎng)極大地提高了NBI 注入的高能粒子能量,并產(chǎn)生高能氘離子尾巴.圖3 給出了ICRF 二次諧波與NBI 協(xié)同加熱的一些實(shí)驗(yàn)設(shè)置和結(jié)果.結(jié)果表明,通過(guò)1.5 MW 的ICRF 和1.0 MW 的NBI 協(xié)同加熱,可以將等離子體極向比壓、等離子體儲(chǔ)能、離子溫度和中子產(chǎn)額分別提高約35%,33%,22%和80%.實(shí)驗(yàn)中還發(fā)現(xiàn)可以通過(guò)降低少子含量、優(yōu)化共振層位置、提高ICRF 功率、提高NBI 束壓、改變NBI 注入角度等方法來(lái)增大NBI 高能離子吸收的ICRF 波能量,有效提高ICRF 與NBI 的協(xié)同效應(yīng),增大高能粒子能量,并獲得更高參數(shù)等離子體.

    圖3 (a) ICRF 二次諧波與NBI 協(xié)同加熱的實(shí)驗(yàn)設(shè)置;(b) 協(xié)同加熱對(duì)離子溫度的影響;(c) 協(xié)同加熱對(duì)中子產(chǎn)額的影響.圖中NBI 的功率設(shè)置成兩個(gè)脈沖是為了獲得更多的協(xié)同加熱和NBI 單獨(dú)加熱的數(shù)據(jù).設(shè)置NBI blips 是為了獲得FIDA 的診斷數(shù)據(jù)Fig.3.(a) Experimental setups of the ICRF second harmonic heating and NBI synergetic heating;(b),(c) influences of synergetic heating on the (b) ion temperature and (c) neutron yield.The NBI heating power is set with two separate pulses in order to obtain more data for cases with synergetic heating and NBI heating alone.The NBI blips is set to facilitate the FIDA diagnostics.

    在實(shí)驗(yàn)中,通過(guò)NES 和FIDA 測(cè)得了高能中子關(guān)于能量的分布,并用FIDA 測(cè)得了D 離子在低能段(10—80 keV)的能譜.從這些實(shí)驗(yàn)可知,在有協(xié)同加熱時(shí),高能中子產(chǎn)生了能量大于3 MeV的高能尾巴.然而,高能粒子的位置和速度分布等卻無(wú)法從實(shí)驗(yàn)診斷中獲得.因此,使用實(shí)驗(yàn)測(cè)得的電流、環(huán)電壓、磁場(chǎng)以及電子、離子密度和溫度等作為TRANSP 的輸入,并計(jì)算了不同ICRF 諧波和NBI 協(xié)同加熱產(chǎn)生的高能粒子分布.不管是對(duì)二次諧波加熱還是三次諧波加熱模式,都確保模擬得到的中子產(chǎn)額和實(shí)驗(yàn)測(cè)量在定量上保持一致.

    在TRANSP 模擬中,首先計(jì)算了在協(xié)同加熱條件下ICRF 波場(chǎng)對(duì)等離子體中各種帶電粒子的加熱情況,結(jié)果如圖4(a)所示.在三次諧波加熱中,ICRF 的能量基本被NBI 高能氘離子吸收了.NBI 高能離子吸收的ICRF 功率密度如圖4(b)所示.由于約80%的NBI 高能離子能量在66 keV,這些高能離子會(huì)使共振層形成比較大的多普勒展寬,并在芯部形成較寬的波能量沉積區(qū)域.而對(duì)于背景的氘離子,雖然其共振層也在磁軸處,但由于其1 keV 的能量導(dǎo)致其回旋半徑和多普勒展寬都較小,因而在有NBI 高能離子競(jìng)爭(zhēng)加熱的情況下,其吸收的波能量并不顯著.此外,氫少子其共振層不在裝置內(nèi),故其吸收的波能量可忽略不計(jì).背景電子則通過(guò)朗道加熱(電場(chǎng)平行分量)和磁泵加熱(平行磁場(chǎng)梯度)來(lái)吸收很少的波能量.

    圖4 (a) 不同粒子吸收的ICRF 三次諧波能量,可以看出ICRF 能量主要被NBI 高能氘離子吸收;(b) NBI 高能離子吸收的三次諧波能量密度Fig.4.(a) ICRF third harmonic heating power absorbed by different charged particles.The dominant ICRF power is absorbed by NBI fast Deuterium ions.(b) The ICRF power density absorbed by NBI fast Deuterium ions.

    基于ICRF 三次諧波與NBI 協(xié)同加熱實(shí)驗(yàn)數(shù)據(jù)(#102242),用TRANSP 模擬得到的氘高能粒子在能量和投擲角上的分布如圖5 所示.結(jié)果表明,相對(duì)于2 MW NBI 加熱,2 MW NBI 與1 MW ICRF 的協(xié)同加熱可以大幅提高氘離子分布的高能尾巴,甚至將初始60 keV 的NBI 高能氘離子的能量提升至約600 keV.這些被加速的離子投擲角的主要區(qū)間為v///v=0.2—0.4.在離子能量E< 180 keV 時(shí),該投擲角區(qū)間的離子主要為捕獲粒子.而當(dāng)E> 320 keV 時(shí),該投擲角區(qū)間的離子則基本為通行粒子.對(duì)于捕獲離子,一般是當(dāng)其反彈點(diǎn)位于共振層(包括多普勒展寬)時(shí),其被波電場(chǎng)加速.對(duì)于通行離子,當(dāng)其軌道穿過(guò)共振層被波電場(chǎng)加速,且在極向上越靠近中平面時(shí)獲得的加速也越大.這是因?yàn)樽畲蟮牟üβ食练e位置在中平面附近的區(qū)域,如圖4 所示.

    圖5 (a1) 2 MW NBI 和(a2) 2 MW NBI+1 MW ICRF 協(xié)同加熱下的高能氘離子在能量-投擲角上的二維分布以及(a3) 高能氘離子關(guān)于能量的一維分布,其中F 是高能氘離子的密度,E 是能量;計(jì)算這些分布時(shí)對(duì)不同徑向位置做了平均;(b1)—(b3) 2 MW NBI+1 MW ICRF 協(xié)同加熱下高能粒子在不同徑向位置的分布Fig.5.(a1)2 MW NBI and (a2)2 MW NBI +1MW ICRF synergetic heating induced fast Deuterium ion distribution as a function of the energy andpitch angle and(a3)1D distribution as a function of energy.(b1)-(b3)2 MWNBI+1 MW ICRF synergetic heating induced fast Deuterium ion distribution at different radial locations.

    通過(guò)TRANSP 計(jì)算的高能粒子分布如圖5 所示.其中,圖5(a1)—(a3)中高能粒子的分布在徑向上做了一個(gè)平均.為了弄清在不同徑向位置的分布,分別計(jì)算了不同徑向位置的高能氘離子在能量-投擲角空間的分布,如圖5(b1)—(b3).顯然,在越靠近等離子體芯部(如ρ=0.1)的位置,高能氘離子被加速的份額越大,且其在投擲角上的分布更寬.這說(shuō)明不僅是反彈點(diǎn)位于共振層的捕獲離子獲得了加速,而且與磁場(chǎng)同向的通行粒子由于其軌道可以長(zhǎng)時(shí)間停留在共振層內(nèi),因而也獲得了加速.而在靠近等離子邊界(如ρ=0.7)的位置,一方面由于在這個(gè)徑向區(qū)間的離子數(shù)相對(duì)較少,另一方面由于捕獲離子其反彈點(diǎn)不能經(jīng)過(guò)共振層,或者捕獲和通行離子穿過(guò)共振層的軌道遠(yuǎn)離中平面,因而被加速的高能粒子較少.

    通過(guò)對(duì)高能離子分布作進(jìn)一步的統(tǒng)計(jì)分析,可以得到高能氘離子關(guān)于投擲角和徑向位置的分布,如圖6 所示.該結(jié)果進(jìn)一步證實(shí)了ICRF 加速的NBI 高能離子主要集中在芯部,且加速的離子的投擲角多為v///v=0.2—0.4.總體來(lái)說(shuō),由于初始的NBI 高能離子主要為通行離子,所以NBI 與ICRF協(xié)同加熱后的高能粒子也以通行離子為主.

    圖6 NBI 單獨(dú)加熱(藍(lán)色)及NBI 與ICRF 協(xié)同加熱(紅色)情況下產(chǎn)生的高能氘離子關(guān)于(a) 徑向和(b) 投擲角的分布.圖中藍(lán)色和紅色重疊區(qū)域顯示為深紅色Fig.6.NBI heating alone (blue) and NBI-ICRF synergetic heating (red) induced fast Deuterium ion distribution as a function of(a) radial position and (b) pitch angle.The dark red color in the figure is due to the overlap of the blue and red colors.

    進(jìn)一步地,可以計(jì)算得到不同能量段的高能氘離子在磁矩-環(huán)向動(dòng)量上的分布,從而理解產(chǎn)生與損失的是哪種類型的離子,結(jié)果如圖7 所示.從圖7 可以看出,被約束或損失的高能粒子可以分為捕獲、同向通行(與等離子體電流同向)和反向通行三大類.顯然,對(duì)E=60 keV 和100 keV 附近的高能粒子,其都在捕獲和同向通行區(qū)間.由于NBI單獨(dú)加熱產(chǎn)生的離子最大能量為80 keV,所以E>80 keV (如圖7(b)中的E=100 keV)的高能粒子只有NBI 和ICRF 協(xié)同加熱才會(huì)產(chǎn)生.從結(jié)果可以看出,協(xié)同加熱對(duì)捕獲和通行離子都有加速作用.而對(duì)于能量更高的離子(如E> 400 keV),其只落在同向通行區(qū)間,即協(xié)同加熱只加速了通行離子.這是因?yàn)楦吣芰康碾x子更容易穿透環(huán)向磁場(chǎng)產(chǎn)生的磁阱而成為通行粒子.

    4 高能粒子輸運(yùn)與損失

    通過(guò)TRANSP 計(jì)算得到高能粒子分布后,可以將其代入粒子軌道程序如SOFT[22]和ISSDE[23]中做進(jìn)一步的計(jì)算,從而分析高能粒子的軌道及其在第一壁上的損失.在該計(jì)算中,考慮了包括主限制器以及ICRF 和低雜波天線限制器的第一壁結(jié)構(gòu).損失機(jī)制包含了碰撞和環(huán)向波紋場(chǎng),模擬時(shí)間為0.5 s 時(shí).計(jì)算結(jié)果表明,NBI 與ICRF 協(xié)同加熱相對(duì)于NBI 單獨(dú)加熱產(chǎn)生了更多的高能粒子,且引起的高能粒子損失也更大.NBI 單獨(dú)加熱可引起約3%的高能粒子損失,而NBI 與ICRF 協(xié)同加熱(H 少子含量為1%)可引起約4%的高能粒子損失.由于主限制器以及ICRF 與LH 天線限制器最靠近等離子體,因而高能粒子也主要損失在這些第一壁部件上.大部分高能粒子損失在中平面上,而剩余的少部分高能粒子則損失在第一壁的中上平面上,如圖8(a)所示.這些損失高能粒子的初始位置位于低場(chǎng)側(cè),損失軌道大部分為捕獲粒子軌道.

    圖8 NBI 與ICRF 協(xié)同加熱產(chǎn)生的高能氘離子 (a)初始位置(藍(lán)色)與第一壁損失(紅色);(b) 在第一壁上的能量沉積;(c) 實(shí)驗(yàn)上觀測(cè)到的 Dα 信號(hào).其中,在SOFT 和ISSDE 模擬計(jì)算中都只考慮了碰撞和環(huán)向波紋場(chǎng)Fig.8.(a) Initial position (blue) and first wall loss (red) of the NBI-ICRF synergetic heating induced fast Deuterium ions;(b) power deposition of the fast ions on the first wall and its comparison with (c) observed Dα signal in experiments.Here,both collision and toroidal field ripple are considered in the SOFT and ISSDE calculations.

    從三維粒子軌道程序ISSDE 的模擬結(jié)果(圖8(b))也可以看出,高能粒子主要損失在中平面上,且會(huì)在限制器的中平面位置引起明顯的熱斑.高能粒子在主限制器和ICRF 限制器上產(chǎn)生的能量沉積與EAST 上由可見光相機(jī)觀測(cè)得實(shí)驗(yàn)結(jié)果(圖8(c))一致.由于主限制器相對(duì)于I 窗口的ICRF 天線限制器更靠近等離子體,因而主限制器上產(chǎn)生的熱斑更大.而由于等離子體電流為逆時(shí)針?lè)较?ICRF 天線左限制器上高能粒子的能量沉積大于右限制器上的能量沉積.

    除了環(huán)向波紋場(chǎng),磁流體不穩(wěn)定性也是引起高能粒子損失的重要原因.為了理解它對(duì)高能粒子損失軌道的影響,使用了粒子導(dǎo)心軌道程序ORBIT[24]對(duì)無(wú)磁擾動(dòng)、有磁擾動(dòng)和同時(shí)存在磁擾動(dòng)和波紋場(chǎng)這三種情況進(jìn)行了計(jì)算.在該計(jì)算中,高能氘離子的能量為60 keV,初始位置為0.8ψ,投擲角(v///v)為0.3.其中ψ為歸一化磁通,v//為平行于磁力線的速度.磁擾動(dòng)使用了理想磁島模型[25],磁擾動(dòng)幅度為 2×10-4,頻率為10 kHz,模數(shù)為(2,1)模.計(jì)算結(jié)果表明,該捕獲離子在磁擾動(dòng)的影響下,其粒子軌道會(huì)在徑向上發(fā)生向外的偏移,并最終損失在第一壁.軌道偏移最大的地方為離子進(jìn)動(dòng)到捕獲軌道反彈點(diǎn)時(shí).當(dāng)磁擾動(dòng)和環(huán)向波紋場(chǎng)同時(shí)存在時(shí),粒子軌道的徑向偏移量更大,粒子損失得越快(如圖9).從高能粒子損失方面來(lái)說(shuō),這兩種機(jī)制也存在協(xié)同作用.因此,在磁約束聚變等離子體中,當(dāng)環(huán)向電流線圈及其產(chǎn)生的環(huán)向波紋場(chǎng)不可改變時(shí),應(yīng)盡量減小磁流體不穩(wěn)定性(包括各種高能粒子模)及其引起的磁擾動(dòng),從而減小高能粒子損失.

    圖9 單個(gè)高能D 離子(E=60 keV)在(a)無(wú)磁擾動(dòng),(b)有磁擾動(dòng),(c)有磁擾動(dòng)和環(huán)向波紋場(chǎng)時(shí)的粒子導(dǎo)心軌道Fig.9.Guiding center orbit of a fast Deuterium ion with energy of E=60 keV in the presence of (a) no magnetic field perturbation,(b) magnetic field perturbation,(c) both magnetic field perturbation and toroidal magnetic field ripple.

    5 結(jié)論

    ICRF 與NBI 不僅是現(xiàn)有磁約束聚變裝置廣泛使用的兩種加熱方法,也是未來(lái)聚變堆將使用的主要加熱方法.ICRF 與NBI 的協(xié)同加熱可以產(chǎn)生大量的高能粒子,這些高能粒子在提高等離子體參數(shù)的同時(shí)也能提高聚變反應(yīng)率和中子產(chǎn)額.然而,未經(jīng)充分慢化的高能粒子損失會(huì)對(duì)裝置器壁造成損傷.因此,協(xié)同加熱條件下高能粒子的產(chǎn)生和輸運(yùn)一直都是磁約束核聚變研究的重點(diǎn).

    本文首先對(duì)ICRF 與NBI 協(xié)同加熱的基本原理做了介紹.由于有限拉莫爾半徑效應(yīng),ICRF 波的電場(chǎng)梯度可有效加速位于高次諧波共振層的NBI高能離子,從而實(shí)現(xiàn)協(xié)同加熱.EAST 實(shí)驗(yàn)表明,通過(guò)1.5 MW 的ICRF 和1.0 MW 的NBI 協(xié)同加熱,可以將等離子體極向比壓、等離子體儲(chǔ)能、離子溫度和中子產(chǎn)額分別提高約35%,33%,22%和80%.同時(shí),協(xié)同加熱也大量產(chǎn)生了能量在約100 keV 量級(jí)的高能粒子,形成了高能氘離子尾巴.

    通過(guò)基于實(shí)驗(yàn)的TRANSP 模擬可知,初始能量為60 keV 的NBI 高能氘離子可被協(xié)同加熱提高至600 keV.協(xié)同加熱下加速的高能離子主要集中在等離子體芯部,且投擲角主要為v///v=0.2—0.4.被加速后的離子在E< 400 keV 時(shí)既有捕獲離子,又有同向通行離子.而被加速后的離子在E>400 keV 時(shí)則只有同向通行離子.這是因?yàn)閷?duì)于投擲角相同的離子來(lái)說(shuō),當(dāng)粒子能量越大時(shí),其平行方向的速度也越大,粒子更容易穿透磁阱而成為通行離子.研究中還發(fā)現(xiàn),當(dāng)EAST 使用常用的環(huán)向磁場(chǎng)(Bt=2.5 T)時(shí),氘離子的二次諧波共振層和氫離子的基頻共振層重疊.此時(shí)當(dāng)氫少子含量越低時(shí),ICRF 的二次諧波與NBI 的協(xié)同加熱效果越強(qiáng),所產(chǎn)生的高能粒子能量也越大.而當(dāng)EAST 降低磁場(chǎng)至Bt=1.7 T 時(shí),氘離子的三次諧波共振層位于磁軸附近,氫離子的基頻共振層位于托卡馬克外.此時(shí)氫少子含量對(duì)ICRF 的三次諧波與NBI的協(xié)同加熱不產(chǎn)生影響.除了改變氫少子含量外,提高ICRF 和NBI 加熱功率、使用ICRF 在軸加熱、優(yōu)化NBI 注入角度等都可以進(jìn)一步提高協(xié)同加熱以及氘離子分布的高能尾巴.

    進(jìn)一步地,我們將TRANSP 計(jì)算得到高能粒子分布代入到粒子軌道程序ORBIT 和ISSDE 中,計(jì)算了高能粒子的輸運(yùn)和損失.研究表明,NBI 與ICRF 協(xié)同加熱產(chǎn)生的更多的高能粒子,因此在考慮存在環(huán)向波紋場(chǎng)和碰撞的情況下,其引起的高能粒子損失也越多.而當(dāng)進(jìn)一步考慮磁流體不穩(wěn)定性時(shí),高能粒子向外偏移的軌道更大,輸運(yùn)更強(qiáng),損失的高能粒子越多.這些損失高能粒子的初始位置都位于低場(chǎng)側(cè),且損失軌道大部分為捕獲粒子軌道.這是因?yàn)槲挥诘蛨?chǎng)側(cè)的粒子才能落到磁阱中成為捕獲粒子.相對(duì)于通行粒子,位于低場(chǎng)側(cè)的捕獲粒子會(huì)在反彈點(diǎn)附近平行速度變得非常小,因而受環(huán)向波紋場(chǎng)(ripple)和碰撞的影響越大,向外的軌道漂移也更大,更容易損失掉.計(jì)算結(jié)果和實(shí)驗(yàn)結(jié)果都表明,高能粒子主要損失在主限制器以及ICRF和低雜波限制器的中上平面.這主要是由于這些限制器的中平面最靠近等離子體,因而最容易受到高能粒子的轟擊.其次,由于等離子體電流從裝置頂部俯視為逆時(shí)針?lè)较?因而同向捕獲粒子軌道的外半圈在極向上是向下的,捕獲粒子更容易落到限制器的中上平面.這些損失的高能粒子是造成限制器上熱斑的主要原因之一.因此,為了保障裝置的安全穩(wěn)定運(yùn)行,亟需減小高能粒子的徑向輸運(yùn)和損失.在未來(lái)的研究中,將進(jìn)一步在實(shí)驗(yàn)和理論上探索控制高能粒子分布以及減小高能粒子損失的方法.

    猜你喜歡
    限制器高能共振
    海上風(fēng)電工程彎曲限制器受力特性數(shù)值模擬研究
    前方高能!戰(zhàn)機(jī)懟瞼
    搞笑秀
    意林彩版(2020年3期)2020-11-22 20:54:12
    電梯或起重機(jī)極限位置限制器的可靠性分析
    安然 與時(shí)代同頻共振
    《高能少年團(tuán)》少年 未來(lái)可期
    選硬人打硬仗——紫陽(yáng)縣黨建與脫貧同頻共振
    新型三階TVD限制器性能分析
    CTA 中紡院+ 化纖聯(lián)盟 強(qiáng)強(qiáng)聯(lián)合 科技共振
    隨車起重機(jī)力矩限制器的振動(dòng)設(shè)計(jì)
    專用汽車(2015年1期)2015-03-01 04:05:29
    欧美+亚洲+日韩+国产| 日韩成人在线观看一区二区三区| av线在线观看网站| 极品少妇高潮喷水抽搐| 午夜免费成人在线视频| 人人妻,人人澡人人爽秒播| 久久久国产欧美日韩av| 免费黄频网站在线观看国产| 嫩草影视91久久| 成人av一区二区三区在线看| 999久久久精品免费观看国产| av免费在线观看网站| 美女视频免费永久观看网站| 久久久久久免费高清国产稀缺| 亚洲成人免费电影在线观看| 九色亚洲精品在线播放| 亚洲欧美激情在线| 韩国av一区二区三区四区| 国产成+人综合+亚洲专区| 久久精品亚洲熟妇少妇任你| 亚洲成人免费电影在线观看| 欧美+亚洲+日韩+国产| 国产淫语在线视频| 天堂√8在线中文| 亚洲精品中文字幕一二三四区| 亚洲精品自拍成人| 香蕉国产在线看| 夜夜夜夜夜久久久久| 很黄的视频免费| 亚洲av美国av| 日本精品一区二区三区蜜桃| 极品教师在线免费播放| 满18在线观看网站| 国产精品久久视频播放| 国产欧美日韩一区二区精品| 精品熟女少妇八av免费久了| 成年女人毛片免费观看观看9 | 国产亚洲欧美精品永久| 正在播放国产对白刺激| 亚洲 欧美一区二区三区| 人妻丰满熟妇av一区二区三区 | 日韩制服丝袜自拍偷拍| 国产乱人伦免费视频| 亚洲精品中文字幕一二三四区| 亚洲aⅴ乱码一区二区在线播放 | 欧美精品亚洲一区二区| 亚洲欧美日韩另类电影网站| a级毛片黄视频| 久久久精品区二区三区| 两性午夜刺激爽爽歪歪视频在线观看 | 十八禁人妻一区二区| 一边摸一边抽搐一进一小说 | 亚洲中文日韩欧美视频| 国产日韩一区二区三区精品不卡| 亚洲欧美日韩高清在线视频| av电影中文网址| 999精品在线视频| 久久久国产欧美日韩av| 亚洲精品美女久久久久99蜜臀| 美女午夜性视频免费| 久久狼人影院| 99在线人妻在线中文字幕 | 777米奇影视久久| 搡老乐熟女国产| 高潮久久久久久久久久久不卡| 两个人看的免费小视频| netflix在线观看网站| 成年人午夜在线观看视频| 激情视频va一区二区三区| 好男人电影高清在线观看| avwww免费| 亚洲精品美女久久av网站| 俄罗斯特黄特色一大片| 啪啪无遮挡十八禁网站| 又大又爽又粗| 国产国语露脸激情在线看| 黑人猛操日本美女一级片| 日韩精品免费视频一区二区三区| 久久中文看片网| 美女高潮到喷水免费观看| 国产精品九九99| 精品国产乱码久久久久久男人| 久久99一区二区三区| 久99久视频精品免费| 欧美成人免费av一区二区三区 | netflix在线观看网站| 一区二区三区激情视频| av线在线观看网站| www日本在线高清视频| av网站在线播放免费| 精品高清国产在线一区| 女人爽到高潮嗷嗷叫在线视频| 国产一区二区三区视频了| 一级a爱片免费观看的视频| 国产精品香港三级国产av潘金莲| 在线天堂中文资源库| 国产欧美日韩精品亚洲av| 精品高清国产在线一区| 国产精品乱码一区二三区的特点 | 69av精品久久久久久| 免费观看人在逋| 高清毛片免费观看视频网站 | 黄色女人牲交| 国产91精品成人一区二区三区| 欧美黑人欧美精品刺激| 美女高潮到喷水免费观看| 男女之事视频高清在线观看| 欧美日韩视频精品一区| 欧美黄色淫秽网站| www.精华液| 怎么达到女性高潮| 美女福利国产在线| 女同久久另类99精品国产91| 精品亚洲成a人片在线观看| 国产91精品成人一区二区三区| 免费观看人在逋| 人妻久久中文字幕网| 亚洲一区中文字幕在线| 日本黄色日本黄色录像| 可以免费在线观看a视频的电影网站| 欧美激情久久久久久爽电影 | 亚洲国产精品sss在线观看 | 国产欧美日韩一区二区精品| 99国产精品一区二区三区| 久热这里只有精品99| a在线观看视频网站| 免费在线观看日本一区| 乱人伦中国视频| 狠狠婷婷综合久久久久久88av| 亚洲精品国产一区二区精华液| 中文亚洲av片在线观看爽 | 欧美 亚洲 国产 日韩一| 一区二区日韩欧美中文字幕| 欧美日韩福利视频一区二区| 亚洲午夜精品一区,二区,三区| 久久中文看片网| 在线观看免费午夜福利视频| 搡老熟女国产l中国老女人| 国产xxxxx性猛交| 免费看十八禁软件| 国产成人免费无遮挡视频| 在线视频色国产色| 18禁裸乳无遮挡免费网站照片 | 丰满的人妻完整版| x7x7x7水蜜桃| 亚洲中文字幕日韩| 日韩欧美免费精品| 一级毛片女人18水好多| 精品视频人人做人人爽| 无人区码免费观看不卡| 在线十欧美十亚洲十日本专区| 两个人免费观看高清视频| 久久性视频一级片| 757午夜福利合集在线观看| 91麻豆精品激情在线观看国产 | 91精品国产国语对白视频| 精品亚洲成国产av| 久久九九热精品免费| 在线观看免费高清a一片| 天天躁狠狠躁夜夜躁狠狠躁| 欧美日韩成人在线一区二区| 99国产精品99久久久久| 另类亚洲欧美激情| 99久久精品国产亚洲精品| 成年版毛片免费区| 露出奶头的视频| 亚洲综合色网址| 超碰97精品在线观看| 大陆偷拍与自拍| 校园春色视频在线观看| 日本黄色视频三级网站网址 | 国产精品99久久99久久久不卡| 99re6热这里在线精品视频| 老司机午夜十八禁免费视频| 亚洲成人免费电影在线观看| а√天堂www在线а√下载 | 国产精品自产拍在线观看55亚洲 | 高清视频免费观看一区二区| av网站在线播放免费| 精品一区二区三区视频在线观看免费 | 久久亚洲精品不卡| 国产麻豆69| 97人妻天天添夜夜摸| 久久久久精品人妻al黑| 亚洲精品国产一区二区精华液| 欧美亚洲日本最大视频资源| 在线看a的网站| 黄色怎么调成土黄色| 亚洲美女黄片视频| 高清毛片免费观看视频网站 | 国产精品一区二区精品视频观看| 国产av精品麻豆| 国产精品欧美亚洲77777| 久久久久久人人人人人| 一夜夜www| 精品一区二区三区av网在线观看| 人人妻,人人澡人人爽秒播| 精品视频人人做人人爽| 99久久精品国产亚洲精品| 一级黄色大片毛片| 少妇裸体淫交视频免费看高清 | 在线国产一区二区在线| 欧美黄色片欧美黄色片| 可以免费在线观看a视频的电影网站| 日本黄色视频三级网站网址 | 亚洲欧美日韩另类电影网站| 99国产综合亚洲精品| 久久ye,这里只有精品| 国产在线一区二区三区精| a级毛片在线看网站| 黄色怎么调成土黄色| 免费少妇av软件| 成人黄色视频免费在线看| 99国产精品一区二区三区| 99re6热这里在线精品视频| 可以免费在线观看a视频的电影网站| 精品一区二区三区视频在线观看免费 | 人人妻人人爽人人添夜夜欢视频| 中文字幕色久视频| 亚洲av成人不卡在线观看播放网| 国产黄色免费在线视频| 亚洲自偷自拍图片 自拍| 天堂动漫精品| www.精华液| 欧美亚洲日本最大视频资源| 成人免费观看视频高清| 午夜精品在线福利| 欧美黑人精品巨大| 国产成人精品无人区| 国产国语露脸激情在线看| 亚洲 欧美一区二区三区| 亚洲 国产 在线| 黄频高清免费视频| 欧美丝袜亚洲另类 | 久久香蕉精品热| 91成年电影在线观看| 搡老岳熟女国产| 精品少妇久久久久久888优播| 久久午夜综合久久蜜桃| 99久久精品国产亚洲精品| 国产欧美日韩综合在线一区二区| 日韩有码中文字幕| 久久性视频一级片| 男人的好看免费观看在线视频 | 黑人猛操日本美女一级片| 亚洲精品国产色婷婷电影| 怎么达到女性高潮| 啦啦啦在线免费观看视频4| av免费在线观看网站| 美女国产高潮福利片在线看| 精品福利观看| 在线观看免费视频日本深夜| 在线观看www视频免费| ponron亚洲| 成年动漫av网址| 日本精品一区二区三区蜜桃| 国产亚洲精品久久久久5区| 免费观看a级毛片全部| 涩涩av久久男人的天堂| 十八禁高潮呻吟视频| 在线十欧美十亚洲十日本专区| 老司机亚洲免费影院| 黑人巨大精品欧美一区二区mp4| 国产精品二区激情视频| 免费女性裸体啪啪无遮挡网站| 国产极品粉嫩免费观看在线| 怎么达到女性高潮| 午夜两性在线视频| 天天影视国产精品| 搡老岳熟女国产| 岛国毛片在线播放| 又大又爽又粗| 高清av免费在线| 又黄又粗又硬又大视频| 亚洲专区国产一区二区| 午夜精品国产一区二区电影| 999精品在线视频| 亚洲专区中文字幕在线| 午夜福利欧美成人| 国产精品av久久久久免费| 国产伦人伦偷精品视频| 精品人妻熟女毛片av久久网站| 午夜福利免费观看在线| 男女免费视频国产| 亚洲五月婷婷丁香| 99久久国产精品久久久| 久久久国产一区二区| 日韩欧美在线二视频 | 欧美av亚洲av综合av国产av| 国产真人三级小视频在线观看| 成年动漫av网址| 国产精品二区激情视频| 国产欧美亚洲国产| 最近最新免费中文字幕在线| 亚洲精品国产一区二区精华液| 麻豆乱淫一区二区| 国产精品国产av在线观看| 男人的好看免费观看在线视频 | 水蜜桃什么品种好| 激情视频va一区二区三区| 丰满的人妻完整版| 一二三四社区在线视频社区8| 免费看十八禁软件| av天堂久久9| 侵犯人妻中文字幕一二三四区| 久久精品国产清高在天天线| 交换朋友夫妻互换小说| 国产成人影院久久av| 国产精品秋霞免费鲁丝片| 精品久久久久久久久久免费视频 | 国产色视频综合| 精品久久久久久电影网| 亚洲国产欧美网| 99国产极品粉嫩在线观看| 91成年电影在线观看| 久久精品国产99精品国产亚洲性色 | 身体一侧抽搐| 日韩精品免费视频一区二区三区| 男女之事视频高清在线观看| 欧美另类亚洲清纯唯美| 久久草成人影院| xxx96com| 欧美日韩av久久| 波多野结衣av一区二区av| 国产精品九九99| 91国产中文字幕| 在线av久久热| 91老司机精品| 国产精品一区二区精品视频观看| av在线播放免费不卡| 丁香欧美五月| 另类亚洲欧美激情| 51午夜福利影视在线观看| 精品无人区乱码1区二区| 久久国产精品人妻蜜桃| av不卡在线播放| 亚洲欧美日韩另类电影网站| videos熟女内射| 俄罗斯特黄特色一大片| 国产区一区二久久| 午夜成年电影在线免费观看| 免费在线观看日本一区| 欧美日韩av久久| 亚洲美女黄片视频| 亚洲七黄色美女视频| 两个人看的免费小视频| ponron亚洲| 91av网站免费观看| 午夜成年电影在线免费观看| 亚洲专区国产一区二区| 亚洲美女黄片视频| 丝袜美腿诱惑在线| av片东京热男人的天堂| 久久国产亚洲av麻豆专区| 黑人巨大精品欧美一区二区蜜桃| 久久精品人人爽人人爽视色| 超色免费av| 制服诱惑二区| 国产精品欧美亚洲77777| 又黄又爽又免费观看的视频| 两个人免费观看高清视频| 国产精品成人在线| 制服人妻中文乱码| 啦啦啦视频在线资源免费观看| 国产亚洲一区二区精品| 黄色女人牲交| 交换朋友夫妻互换小说| 色综合欧美亚洲国产小说| 国产免费男女视频| 女人高潮潮喷娇喘18禁视频| 91av网站免费观看| 色婷婷av一区二区三区视频| 黄色毛片三级朝国网站| 亚洲精品国产区一区二| a级片在线免费高清观看视频| 亚洲精品自拍成人| av天堂久久9| 一级a爱视频在线免费观看| 国产一区在线观看成人免费| 欧美激情高清一区二区三区| 国产欧美日韩综合在线一区二区| xxxhd国产人妻xxx| 亚洲精品中文字幕在线视频| 亚洲久久久国产精品| 亚洲熟女精品中文字幕| 欧美久久黑人一区二区| 成人av一区二区三区在线看| 18禁国产床啪视频网站| 亚洲av日韩在线播放| 国产亚洲欧美98| 久久精品亚洲熟妇少妇任你| 黄色女人牲交| 99香蕉大伊视频| 国产精品电影一区二区三区 | av在线播放免费不卡| 亚洲精品国产精品久久久不卡| 国产av又大| 王馨瑶露胸无遮挡在线观看| 国产淫语在线视频| 波多野结衣一区麻豆| 亚洲av片天天在线观看| 嫩草影视91久久| 国产精品九九99| 女人精品久久久久毛片| 欧美日韩瑟瑟在线播放| 国产av一区二区精品久久| 后天国语完整版免费观看| 久久青草综合色| 国产精品1区2区在线观看. | 交换朋友夫妻互换小说| 丰满迷人的少妇在线观看| 成人国产一区最新在线观看| 亚洲精品美女久久久久99蜜臀| 美女高潮到喷水免费观看| 午夜精品久久久久久毛片777| 正在播放国产对白刺激| 成年人午夜在线观看视频| 国产在线观看jvid| av超薄肉色丝袜交足视频| 夜夜躁狠狠躁天天躁| 国产精品久久久av美女十八| 国产精品久久视频播放| 啦啦啦视频在线资源免费观看| 亚洲成人免费av在线播放| 国产午夜精品久久久久久| a在线观看视频网站| 两个人看的免费小视频| 国产99白浆流出| 男女午夜视频在线观看| 99久久人妻综合| 黄频高清免费视频| 国产成人啪精品午夜网站| 黑人巨大精品欧美一区二区mp4| 欧美激情极品国产一区二区三区| 国产97色在线日韩免费| 在线观看免费视频日本深夜| 亚洲中文av在线| 国产亚洲欧美精品永久| 日本黄色视频三级网站网址 | 性少妇av在线| 免费在线观看视频国产中文字幕亚洲| 国产精品国产av在线观看| 久久久久久亚洲精品国产蜜桃av| 在线观看免费视频日本深夜| 欧美在线黄色| 亚洲精品中文字幕一二三四区| 成人影院久久| 伊人久久大香线蕉亚洲五| 亚洲国产精品一区二区三区在线| 搡老乐熟女国产| 亚洲 国产 在线| 久久精品国产亚洲av高清一级| 亚洲欧洲精品一区二区精品久久久| 国产精品久久久久久精品古装| 国产精品自产拍在线观看55亚洲 | 国产精品影院久久| 满18在线观看网站| 国产99白浆流出| 高清视频免费观看一区二区| 久久这里只有精品19| 美女国产高潮福利片在线看| 天堂动漫精品| 免费不卡黄色视频| 亚洲五月色婷婷综合| 岛国毛片在线播放| 成熟少妇高潮喷水视频| 亚洲国产欧美日韩在线播放| 91av网站免费观看| 精品第一国产精品| 搡老岳熟女国产| 亚洲欧美激情综合另类| 精品国产亚洲在线| 亚洲av欧美aⅴ国产| 一本大道久久a久久精品| 亚洲精品国产精品久久久不卡| 国产片内射在线| 女人被躁到高潮嗷嗷叫费观| 亚洲欧美一区二区三区黑人| 人人妻人人添人人爽欧美一区卜| 久久久国产欧美日韩av| 欧美精品一区二区免费开放| av免费在线观看网站| 亚洲七黄色美女视频| 精品卡一卡二卡四卡免费| 国产欧美日韩精品亚洲av| 国产片内射在线| 精品国产一区二区三区四区第35| 国产精品成人在线| 激情在线观看视频在线高清 | 50天的宝宝边吃奶边哭怎么回事| 两个人免费观看高清视频| 久久中文字幕人妻熟女| 久久影院123| 一区二区日韩欧美中文字幕| 狂野欧美激情性xxxx| 十八禁高潮呻吟视频| 欧美乱妇无乱码| 欧美激情久久久久久爽电影 | 久久影院123| 丁香欧美五月| 久久午夜亚洲精品久久| 成人影院久久| 黄色女人牲交| 久久精品国产99精品国产亚洲性色 | 欧美日韩黄片免| 桃红色精品国产亚洲av| 国产精品一区二区精品视频观看| 午夜福利免费观看在线| 搡老乐熟女国产| 国产成人系列免费观看| 波多野结衣一区麻豆| 老司机在亚洲福利影院| 久久久国产精品麻豆| 成人18禁高潮啪啪吃奶动态图| 日韩欧美国产一区二区入口| 中国美女看黄片| 午夜福利影视在线免费观看| 亚洲精品美女久久久久99蜜臀| www.精华液| 久久久国产欧美日韩av| 亚洲人成电影免费在线| 一边摸一边做爽爽视频免费| 久久人妻熟女aⅴ| 啦啦啦视频在线资源免费观看| 一进一出好大好爽视频| 久久久久国产一级毛片高清牌| 18禁观看日本| 男人的好看免费观看在线视频 | 亚洲五月婷婷丁香| 亚洲专区国产一区二区| 99国产精品一区二区蜜桃av | videosex国产| 搡老乐熟女国产| 日本欧美视频一区| 三级毛片av免费| 久久久久久人人人人人| a在线观看视频网站| 国产亚洲av高清不卡| 午夜福利在线观看吧| 国产乱人伦免费视频| 人成视频在线观看免费观看| 曰老女人黄片| 亚洲aⅴ乱码一区二区在线播放 | 99国产综合亚洲精品| 精品高清国产在线一区| 一二三四社区在线视频社区8| 丰满人妻熟妇乱又伦精品不卡| 99re在线观看精品视频| 婷婷丁香在线五月| 亚洲国产毛片av蜜桃av| 日本撒尿小便嘘嘘汇集6| 亚洲伊人色综图| 免费不卡黄色视频| 女人高潮潮喷娇喘18禁视频| 热99re8久久精品国产| 身体一侧抽搐| 黄片大片在线免费观看| 亚洲五月色婷婷综合| 亚洲一区二区三区欧美精品| 成年版毛片免费区| 狠狠狠狠99中文字幕| 高清欧美精品videossex| 国产一区二区三区在线臀色熟女 | 丁香六月欧美| 午夜福利免费观看在线| 咕卡用的链子| 黄频高清免费视频| 免费日韩欧美在线观看| 亚洲av第一区精品v没综合| av电影中文网址| 丝袜美腿诱惑在线| 美国免费a级毛片| 99精国产麻豆久久婷婷| 看黄色毛片网站| 久久久国产精品麻豆| 在线国产一区二区在线| 男女之事视频高清在线观看| 热99国产精品久久久久久7| 满18在线观看网站| 欧美不卡视频在线免费观看 | 免费日韩欧美在线观看| 男女之事视频高清在线观看| 亚洲自偷自拍图片 自拍| 性色av乱码一区二区三区2| 99国产精品一区二区三区| 亚洲精品中文字幕一二三四区| 亚洲va日本ⅴa欧美va伊人久久| 精品国产乱码久久久久久男人| 日韩免费高清中文字幕av| 欧美日韩乱码在线| 国产精品乱码一区二三区的特点 | 欧美乱色亚洲激情| 亚洲av电影在线进入| 女人爽到高潮嗷嗷叫在线视频| 免费观看人在逋| 老司机在亚洲福利影院| 最近最新免费中文字幕在线| 亚洲精品一卡2卡三卡4卡5卡| 亚洲aⅴ乱码一区二区在线播放 | 99国产综合亚洲精品| 精品久久久精品久久久| 香蕉丝袜av| 欧美 亚洲 国产 日韩一| 天堂中文最新版在线下载| 又黄又粗又硬又大视频| 亚洲精品中文字幕一二三四区| 国产精品乱码一区二三区的特点 | 在线观看舔阴道视频| 一个人免费在线观看的高清视频| 看片在线看免费视频| 亚洲中文av在线| 99久久99久久久精品蜜桃| 久久中文字幕人妻熟女| 99国产精品一区二区三区| 亚洲aⅴ乱码一区二区在线播放 |