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

    螺旋波等離子體原型實驗裝置中天線的優(yōu)化設(shè)計與功率沉積*

    2019-10-25 06:58:00平蘭蘭張新軍楊樺徐國盛萇磊吳東升呂虹鄭長勇彭金花金海紅何超甘桂華
    物理學(xué)報 2019年20期
    關(guān)鍵詞:拋物面等離子體徑向

    平蘭蘭 張新軍 楊樺 徐國盛 萇磊 吳東升 呂虹 鄭長勇 彭金花 金海紅 何超 甘桂華

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

    2) (中國科學(xué)技術(shù)大學(xué),合肥 230026)

    3) (安徽建筑大學(xué)電子與信息工程學(xué)院,合肥 230601)

    4) (四川大學(xué)空天科學(xué)與工程學(xué)院,成都 610065)

    5) (國防科技大學(xué)電子對抗學(xué)院,合肥 230036)

    近年來,螺旋波等離子體源在核聚變條件下等離子體與材料的相互作用方面具有十分重要的研究意義.本文對高密度螺旋波等離子體原型實驗裝置(helicon physics prototype experiment,HPPX)中的螺旋波系統(tǒng)發(fā)射天線進行了優(yōu)化設(shè)計,利用HELIC程序,對螺旋波的耦合和功率沉積與天線的類型、天線長度、運行頻率等關(guān)鍵參數(shù)的相互關(guān)系作了數(shù)值模擬計算,進而給出了天線的最佳天線結(jié)構(gòu)和物理尺寸.同時還分析了靜磁場強度和軸心等離子體密度對功率沉積及其分布的影響,發(fā)現(xiàn)等離子體對螺旋波的吸收功率在多個靜磁場和軸心密度處有不同的峰值功率點,且整體的耦合趨勢隨靜磁場增大呈上升的趨勢,而隨軸心密度增大是下降的趨勢; 并根據(jù)這些仿真結(jié)果深入探討了螺旋波等離子體的電離機制.為了進一步研究HPPX裝置中螺旋波與等離子體耦合的問題,還給出了等離子體放電時的感生電磁場與電流密度的分布情況.本文工作可為HPPX裝置上螺旋波天線設(shè)計和相關(guān)物理實驗提供理論依據(jù).

    1 引 言

    螺旋波是由射頻(radio frequency,RF)天線激發(fā)的一種在高電導(dǎo)率介質(zhì)中傳播的低頻電磁波,使用螺旋波電離可以產(chǎn)生高密度的等離子體.在軸向均勻磁場條件下,螺旋波等離子體可通過頻率介于離子和電子回旋頻率間的射頻波激發(fā)石英管中的介質(zhì)氣體得到[1].在0.1 Pa的低氣壓下等離子體密度可達1013-1014cm-3量級,并且其電離效率高達100%,這是迄今采用人工方法在低氣壓下所能獲得的最大等離子體密度[2].螺旋波等離子體因其高密度、高效率、均勻平穩(wěn)、無電極污染、粒子能量可控等優(yōu)點,廣泛應(yīng)用于薄膜沉積[3]、等離子體火箭推進[4]、半導(dǎo)體刻蝕[5]和核聚變[6]等領(lǐng)域.

    螺旋波等離子體的電離機制十分復(fù)雜,其核心是波-粒能量耦合機制,自1960年Aigrain[7]首次提出螺旋波概念,1970年Boswell[8]首次實現(xiàn)螺旋波等離子體放電以來,人們對螺旋波產(chǎn)生高密度的等離子體的物理機理、螺旋波能量耦合到等離子體中的物理方式始終沒有給出最終的解釋.1991年,Chen[9]利用朗道阻尼機制導(dǎo)出了均勻有界等離子體中螺旋波的色散關(guān)系,認(rèn)為電子通過朗道阻尼進行加熱具有很高的吸收效率.1999年,Chen和Blackwell[10]又提出了螺旋波放電中朗道阻尼的限制.朗道阻尼可以解釋在一般的放電實驗中的螺旋波等離子體的能量沉積,但在較大尺度實驗裝置中或者在天線長度大于半個波長的情況下,朗道阻尼沉積將不再成立.通過深入研究發(fā)現(xiàn),在螺旋波放電能量沉積過程中,存在著螺旋波之外的一種靜電性質(zhì)的表面波Trivelpiece-Gould (TG)波[11].Shamrai和Taranov[12]認(rèn)為放電過程中RF功率通過兩種通道轉(zhuǎn)移給等離子體,即體吸收的螺旋波(H波)和表面吸收的TG波實現(xiàn)的,這兩種波在等離子體柱徑向位置上進行耦合,而大部分的RF功率通過強阻尼的TG波轉(zhuǎn)移給等離子體.到目前為止,有關(guān)螺旋波高功率耦合到等離子體中的機理和能量沉積的機制仍在不斷探索中.因此,磁場對兩種波的影響需進一步討論,但TG波能量沉積機理的假設(shè)已逐漸成為主流共識[13].

    近年來,螺旋波等離子體的放電模擬取得了較大發(fā)展,促進了人們對螺旋波放電機理的認(rèn)識.如:由Chen和Arnush[14-16]開發(fā)的HELIC代碼,通過電磁場理論分析計算了螺旋波和TG波的能量沉積及其電場、磁場和電流的分布情況; 由Mouzouris和Scharer[17]開發(fā)的ANTENA2代碼利用碰撞機制和朗道阻尼加熱機制,計算了等離子體中的電磁場和功率沉積; 2012年,Melazzi等[18]開發(fā)了SPIREs代碼,采用FDTD算法對圓柱等離子體中的電磁場問題和功率沉積問題進行了快速求解; 2015年,Melazzi和Lancellotti[19]開發(fā)了ADAMANT代碼,該代碼利用MOM算法研究了天線幾何結(jié)構(gòu)及其在功率沉積方面的影響.國內(nèi)方面,成玉國等[20]考察了徑向電子數(shù)密度非均勻分布條件下,不同強度磁場對能量沉積的影響; 楊雄等[21]建立三維條件下螺旋波放電模型,充分考慮了電離過程中的粒子碰撞和電化學(xué)反應(yīng),利用漂移-擴散輸運方程計算等離子體中各粒子密度及電子溫度;萇磊等[22,23]針對徑向密度梯度對于螺旋波等離子體中電磁波和能量吸收進行了數(shù)值模擬.

    HELIC代碼常用于射頻等離子體源設(shè)計和實驗中,該程序利用特定的邊界條件求解6個徑向耦合微分方程,得到兩個獨立的波-螺旋波和TG波.這比將非均勻等離子體劃分成層并在每個界面匹配邊界條件(如ANTENA2代碼)要快得多,而且改變天線的幾何形狀也特別容易.本文首次采用HELIC代碼對高密度螺旋波等離子體原型實驗裝置(helicon physics prototype experiment,HPPX)中等離子體放電現(xiàn)象進行數(shù)值模擬,分析了HPPX裝置中螺旋波天線種類、天線長度及運行頻率等參量對功率沉積的影響,從而得到了發(fā)射天線的最佳結(jié)構(gòu)和尺寸.眾所周知,忽略電子質(zhì)量me的簡單螺旋波理論,只有當(dāng)軸向電場Ez消失時才成立,而實際實驗中并未滿足這一條件.通過對冷等離子體方程中有限的me和Ez進行精確求解,預(yù)測出了TG波的存在,這極大改變了人們對等離子體放電的認(rèn)識.因此,在分析靜磁場大小和軸心等離子體密度對波-粒耦合能量的影響時,本文對等離子體中的螺旋波和TG波的功率沉積機理做了進一步探討,進而,研究了等離子體放電時的感生電磁場與電流密度等參數(shù).本文所研究的成果對HPPX裝置上螺旋波天線設(shè)計和今后相關(guān)物理實驗開展具有一定的指導(dǎo)意義.

    2 理論公式與計算模型

    2.1 理論分析

    為了研究RF天線產(chǎn)生的螺旋波耦合到等離子體中的功率沉積,從麥克斯韋方程出發(fā),得到電場E和磁場B滿足的法拉第定律和安培定律:

    其中,

    式中下標(biāo)α表示粒子種類(如電子和離子),ω為波的頻率,ωpα是等離子體頻率,ωcα是等離子體回旋頻率,υα為粒子間的有效碰撞頻率.考慮在圓柱體坐標(biāo) (r,φ,z)下,外加靜磁場B0與z軸平行,等離子體的擾動與軸向波數(shù)kz和周向波數(shù)m有關(guān),可表示為exp[i(kzz+mφ - ωt)],為了方便,使用k=kz,采用 ?/?z→ik ,?/??→im ,?/?t→iω 進行傅里葉變換,(1)式和(2)式可轉(zhuǎn)化成以下6個微分方程組形式[14]:

    其中k0=ω/c為真空中的波數(shù); Er,Eφ,Ez,Br,Bφ,Bz分別為徑向、角向、軸向感生電磁場.公式(7)-(12)給出了螺旋波等離子體中波場理論模型,對上述6個方程進行編程求解,生成基函數(shù)br,bz,er和ez.如果等離子體是均勻的,在圓柱體坐標(biāo)中,基函數(shù)是貝塞爾函數(shù),即bn,z(r)=Jm(Tnr) (n=1,2),根據(jù)螺旋波理論有

    其中R,L=(S ± D)/2,根據(jù)邊界條件,對于T2值,方程(13)有兩個根,一個是位于軸上頻率較低的H波,另一個是位于等離子體邊界處頻率較高的TG波.當(dāng)?shù)入x子體非均勻時,靠近原點處表現(xiàn)為貝塞爾函數(shù)Jm(r)的波被識別為H波,另一個根為TG波.

    假設(shè)天線是一個半徑r=Rb的無限薄的薄片,厚度可忽略不計,其上的表面電流J滿足 ?·J=0 ,因此傳送到天線上的功率為

    從時間平均能量流的角度來看,天線終端的輸入功率 Pant主要由天線材料中的電阻損耗 Pcopper,等離子體吸收的功率 Pabs,以及天線輻射到真空中的功率損耗 Pspace組成,即

    本文所研究的螺旋波放電是在密封的真空腔內(nèi).(15)式中最后一項為零,因為任何逃逸粒子都會被這個邊界條件反射回來.此外,螺旋波天線被認(rèn)為是理想的導(dǎo)體,因此 Pcopper=0.在這些條件下,天線輸入的時間平均功率等于等離子體吸收的時間平均功率,即 Pant=Pabs.

    對于等離子體的吸收功率,定義特定的等離子體功率譜函數(shù) Sp(k′)[14],

    等離子體電流密度可用冷等離子體張量形式給出[15]:

    這里等離子體電場及電流密度是由頻譜為Kφ(k′)=δ(k-k′)的天線激勵,由于 Sp(k′) 是一個偶函數(shù),可以用天線功率譜函數(shù)pA(k)=|Kφ(k)|2+|Kφ(-k)|2簡化計算.在k到k+dk,及-k到-(k+dk)之間,對于電流為1 A的等離子體的功率譜函數(shù) Pp(k) ,有

    則總吸收功率為

    因此當(dāng)電流為1 A時,等離子體吸收功率數(shù)值上等于等離子體負(fù)載阻抗 Rp=Pabs/I2.

    2.2 計算模型

    HPPX是2016年中國科學(xué)院合肥等離子體研究所自主研發(fā)的螺旋波放電裝置.該裝置采用高功率射頻源和模塊化的設(shè)計,束流密度高達1023m-2·s-1量級,可產(chǎn)生高密度高能量等離子體,是用于研究核聚變的新型材料裝置(圖1).該裝置總長度為4 m,分為4個腔室,其徑向中心處是等離子體,等離子體外有個法拉第屏蔽罩,纏繞在屏蔽罩外是螺旋波天線,放置在等離子體軸向中心處(z=0 m),用以產(chǎn)生螺旋波來激發(fā)等離子體,最外面是一個真空腔.等離子體參數(shù)和邊界條件如表1所列.

    圖1 HPPX裝置結(jié)構(gòu)圖Fig.1.The structure of HPPX.

    表1 HPPX裝置參數(shù)及等離子體參數(shù)Table 1.HPPX device parameters and plasma parameters.

    磁場密度沿軸向為均勻的,粒子的密度分布沿軸向也是均勻的,徑向密度可以由下面函數(shù)確定:

    其中s和t是常量,fa表示r=Ra處的相對密度n/ne0,ne0是等離體中心處的密度.如果fa=0,則w=Ra,密度函數(shù)就只有兩個參數(shù),如果fa≠0,密度函數(shù)可以設(shè)置為各種密度形式,例如拋物面函數(shù)和高斯函數(shù)等類型.

    3 結(jié)果與討論

    本文基于上述的理論公式和計算模型,對HPPX裝置進行數(shù)值模擬.構(gòu)建2種徑向密度結(jié)構(gòu):拋物面分布(s=2,t=1,fa=0.1)和高斯分布(s=2,t=6,fa=0.01),如圖2所示.這兩種密度分布雖然在軸向上有相同的峰值,但徑向上卻有不同的密度梯度,這對螺旋波與TG波的非共振模式轉(zhuǎn)換和能量吸收有著非常重要的影響.在HPPX裝置的相關(guān)物理實驗中可以通過調(diào)節(jié)磁場位型等參數(shù),獲得所需要的等離子體徑向密度分布,文獻[24]顯示,較強的磁場對等離子體密度產(chǎn)生較好的約束,其特征半徑較小.將朗繆爾探針測得的數(shù)據(jù)進行擬合便可確定實驗中加載的等離子體密度結(jié)構(gòu).拋物面分布和高斯分布是螺旋波等離子體源實驗中非常重要的徑向密度結(jié)構(gòu),用以分析天線與等離子體間的功率耦合過程,它對HPPX這種大功率的螺旋波等離子體裝置的實驗設(shè)計,有著非常重要的指導(dǎo)價值.

    圖2 等離子體徑向密度分布 (a) 拋物面分布; (b)高斯分布Fig.2.Radial profiles of plasma density:(a) Parabolic density profile; (b) Gaussian density profile.

    螺旋波屬于哨聲波,對于絕緣介質(zhì)為邊界、半徑為Ra的等離子體中產(chǎn)生的m=1螺旋波,有如下關(guān)系成立[25]:

    由此可見,對于HPPX裝置中天線的設(shè)計除了與等離子體的半徑有關(guān),還與RF源運行頻率f,電子密度ne,磁場強度B等參數(shù)有關(guān).接下來,本文將從以下這幾個方面來分析.

    3.1 螺旋波天線種類

    螺旋波通過TG-H耦合模式將能量傳輸給電子,要求天線能很好地將射頻源供給的能量耦合傳遞給螺旋波,因此天線的選擇至關(guān)重要.目前,用于直線磁約束產(chǎn)生等離子體的螺旋波天線主要有3種,即Nagoya III型[26]、Boswell型[27]和Half helix型[28].這3種天線可有效激發(fā)m=+1的波動,其能量沉積主要是角向模式[29].因其結(jié)構(gòu)簡單,可激發(fā)高效螺旋波,因此逐漸應(yīng)用于各類螺旋波等離子體實驗裝置中.從圖3(a)中可以看出,在拋物面密度分布下,Half helix型天線在等離子體中心及邊緣處較Nagoya III型、Boswell型天線有較高的徑向相對吸收功率,且Nagoya III型天線和Boswell型天線在等離子體中耦合的效果差不多.沿軸向的相對吸收功率(圖3(c))更能明顯地看出各種天線的吸收功率大小,Half helix型天線產(chǎn)生的電磁波在等離子體中的相對吸收功率最高,其次是Nagoya III型天線,最差的是Boswell型天線.

    圖3(b)和圖3(d)是高斯密度分布下的相對吸收功率,從徑向及軸向上都可以看出,吸收效果最好的是Half helix型天線,最差的仍是Boswell型天線.而且在這兩種密度分布下,除了Half helix型天線外,其余兩種天線都關(guān)于天線的中心處對稱.這種非軸對稱的能量耦合主要與Half helix天線自身結(jié)構(gòu)和靜磁場的方向有關(guān),在很多螺旋波設(shè)備中[23]都被觀察到過,且高斯密度分布下的徑向相對功率沉積高于拋物面下的.

    通過以上分析,可以得出,Half helix天線產(chǎn)生的是非軸對稱的射頻能量耦合,較其他兩種天線能激發(fā)產(chǎn)生出更高的電子密度.因此,在HPPX裝置中選用Half helix型天線作為激勵天線.

    圖3 3種典型的螺旋波天線的徑向(z=0.2 m)和軸向(r=0.02 m)相對吸收功率 (a)拋物面密度分布下3種天線的徑向相對吸收功率; (b)高斯密度分布下3種天線的徑向相對吸收功率; (c)拋物面密度分布下3種天線的軸向相對吸收功率; (d)高斯密度分布下3種天線的軸向相對吸收功率Fig.3.Relative power absorption in radial (z=0.2 m) and axial (r=0.02 m) directions for three typical helicon wave antennas:(a) Radial relative absorption power of three antennas under parabolic density distribution; (b) radial relative absorption power of three antennas under Gaussian density distribution; (c) axial relative absorption of power of three antennas under parabolic density distribution; (d) axial relative absorption of power of three antennas under Gaussian density distribution.

    3.2 天線長度

    天線的長度也是HPPX裝置中天線設(shè)計的一項重要指標(biāo).在密度均勻的等離子體中,軸向螺旋波的波長由(23)式得到:λz≈0.1862 m,由文獻[30]可知,天線和螺旋模式耦合較好時,kz≈π/LA,3π/LA等,對應(yīng)的波長為λz≈2LA,2LA/3,2LA/5等.根據(jù)HPPX裝置尺寸,取λz≈2LA/5,則天線長度為LA≈5λz/2≈0.4655 m,考察不同天線長度下的吸收功率,如圖4所示.

    從圖4(a)可以看出,在拋物面密度分布下等離子體在中心及邊緣的相對吸收功率比中間的要高很多,并且隨著天線長度的增大,徑向相對吸收功率也慢慢變大,當(dāng)天線長度等于0.4 m時,相對吸收功率達到最大,但之后天線長度再增大,相對吸收功率卻隨之減小.為了更清楚地看清天線長度對吸收功率的影響,利用公式將每一個天線長度對應(yīng)的相對吸收功率沿徑向進行積分得到總的徑向吸收功率(圖4(c)),可明顯看出天線長度在0.4 m時吸收功率最大,基本上和計算所得的理論值吻合.

    在高斯密度分布下,等離子體的徑向相對吸收功率在等離子體中心處且LA=0.5 m時較強些(圖4(b)).同樣利用公式可得到圖4(d),從圖4(d)中可以更加清晰地看出天線長度在0.5 m時吸收功率最好.總體來說,高斯分布下的吸收功率比拋物線分布時的吸收功率要好,但在實際HPPX實驗中,常用的是拋物面形式的密度分布,因此HPPX裝置中天線長度選擇0.4 m.

    3.3 運行頻率

    天線的運行頻率是一個非常重要的參數(shù),它影響著天線發(fā)射出的螺旋波在等離子體中的耦合效果.選擇6種典型的射頻頻率仿真對比它們的相對吸收功率(圖5).從圖5(a)可以看出,隨著頻率的增大,相對吸收功率在等離子體邊緣和中心處都逐漸減小,而中心處的相對吸收功率降低更為明顯,說明TG波比H波在高頻處能吸收更多的功率,這可能是由于運行頻率增加導(dǎo)致快速的電子運動和強烈的靜電加熱造成的.減少徑向壓力,使碰撞阻尼增大從而增加電子自由路徑,可能有助于將能量吸收從等離子體邊緣重新分配到中心處.圖5(b)顯示,隨著頻率增加,總體上等離子體中心處的徑向吸收明顯降低,而邊緣處變化不大,這是由于高斯型密度分布下邊緣處密度梯度變化較拋物面型分布平緩許多.從這2幅圖中可以看出,并不是運行頻率越高,相對吸收功率就越高.在徑向吸收功率上,f=13.56 MHz時的吸收最好,之后隨著運行頻率的增加,相對吸收功率下降得很快,這從圖5(c)和(d)中也可以明顯地看出來.因此,在HPPX實驗中13.56 MHz是螺旋波等離子體常用的頻率[31].再次,從圖5中可以看出高斯分布下的等離子體的相對吸收功率要比拋物面時的大,這和前面提到的基本一致.

    圖4 不同天線長度下的徑向(z=0.2 m)吸收功率 (a)拋物面密度分布下天線長度對徑向相對吸收功率的影響; (b)高斯密度分布下天線長度對徑向相對吸收功率的影響; (c)拋物面密度分布下徑向吸收功率隨天線長度的變化曲線; (d)高斯密度分布下徑向吸收功率隨天線長度的變化曲線Fig.4.Relative power absorption in radial (z=0.2 m) directions for different antenna lengths:(a) Effect of antenna length on radial relative absorption power under parabolic density distribution; (b) effect of antenna length on radial relative absorption power under Gaussian density distribution; (c) radial relative absorption power of different antennas lengths under parabolic density distribution; (d) radial relative absorption power of different antennas lengths under Gaussian density distribution.

    圖5 不同運行頻率下的徑向(z=0.2 m)吸收功率 (a)拋物面密度分布下天線運行頻率對徑向相對吸收功率的影響; (b)高斯密度分布下天線運行頻率對徑向相對吸收功率的影響; (c) 拋物面密度分布下徑向吸收功率隨運行頻率的變化曲線; (d)高斯密度分布下徑向吸收功率隨運行頻率的變化曲線Fig.5.Relative power absorption in radial (z=0.2 m) directions for various operating frequencies:(a) Effect of various operating frequencies on radial relative absorption power under parabolic density distribution; (b) effect of various operating frequencies on radial relative absorption power under Gaussian density distribution; (c) radial relative absorption power of various operating frequencies under parabolic density distribution; (d) radial relative absorption power of various operating frequencies under Gaussian density distribution.

    3.4 靜磁場強度

    等離子體中引入外部約束靜磁場,使得螺旋波深入等離子體柱中傳播,獲得了高的功率耦合效率.靜磁場大小的改變對帶電粒子的運動軌跡產(chǎn)生影響,改變了各種粒子間的碰撞頻率,從而致使等離子體介電張量隨之變化.同時靜磁場對電子與離子繞磁力線的回旋尺度的影響導(dǎo)致帶電粒子的能量吸收區(qū)域發(fā)生改變,從而產(chǎn)生波在等離子體中能量沉積的不均勻性問題,可見磁場對等離子體的影響比較大.設(shè)定磁場值從100 Gs (1 Gs=10-4T)到1000 Gs之間變化,采樣間隔為50 Gs,計算得到等離子體的相對徑向吸收功率(圖6).

    不同的磁場條件下,螺旋波和TG波的徑向波長不同且所受阻尼不同,其在等離子體中沿徑向的傳播距離或者功率耦合范圍就會有較大的差異.在磁場很低時,由于磁場太小,從色散關(guān)系上分析TG波被禁止傳播,此時不存在TG波(圖6(a)).隨著磁場強度增大,在等離子體邊緣處的相對吸收功率迅速增大,這是由于在強磁場下TG波所受的阻尼較大,徑向波長較短,很難在等離子體柱中傳播,所以TG波影響范圍由徑向全場縮減至邊界的極小范圍[20],同時能量的吸收逐漸集中于等離子體-真空邊界處,說明外加穩(wěn)恒磁場較大不利于TG波的深入傳播,較高的磁場限制TG波的傳播,從而造成徑向能量分布的嚴(yán)重不均.與拋物面分布不同的是,在高斯分布中(圖6(b)),H波所受阻尼小,在等離子體內(nèi)的穿透能力強,可沿徑向傳播距離遠(yuǎn),等離子體中的功率沉積主要發(fā)生在等離子體中心位置.隨著磁場的增大,TG波和H波耦合模式發(fā)生轉(zhuǎn)變,TG波會被抑制,H波可以在磁化等離子體柱內(nèi)部進行傳播,將能量耦合給等離子體.

    從圖6(c)和圖6(d)中可以看出,在某些磁場強度處有一些功率峰值點,說明該處天線和等離子體有較強的耦合,而總體的耦合程度在上升.在HPPX裝置的實驗設(shè)計中,可以參考仿真結(jié)果,選擇合適的磁場值,獲得更高的等離子體吸收功率.

    圖6 不同靜磁場下的徑向(z=0.2 m)吸收功率 (a)拋物面密度分布下磁場對徑向相對吸收功率的影響; (b)高斯密度分布下磁場對徑向相對吸收功率的影響; (c)拋物面密度分布下徑向吸收功率隨磁場強度的變化曲線; (d)高斯密度分布下徑向吸收功率隨磁場強度的變化曲線Fig.6.Relative power absorption in radial (z=0.2 m) directions for various static magnetic intensity:(a) Effect of various magnetic intensity on radial relative absorption power under parabolic density distribution; (b) effect of various magnetic intensity on radial relative absorption power under Gaussian density distribution; (c) radial relative absorption power of various magnetic intensity under parabolic density distribution; (d) radial relative absorption power of various magnetic intensity under Gaussian density distribution.

    3.5 軸心等離子體密度

    同樣,徑向等離子體密度的非均勻性導(dǎo)致徑向擾動電場的產(chǎn)生,造成帶電粒子的徑向漂移,從而改變了各種粒子間的碰撞頻率,進而使等離子體色散方程中介電張量也隨之改變,最終對等離子體的徑向功率沉積產(chǎn)生影響.設(shè)定密度值從1×1012到1×1013cm-3之間變化,每隔0.3×1012cm-3進行一次采樣,計算得到等離子體吸收功率(圖7).在拋物面密度分布下(圖7(a)),當(dāng)密度較低時,等離子體中心處的吸收功率較低,大部分功率集中在等離子體邊緣部分.這是由于TG波強烈的邊緣加熱,TG波是短波長準(zhǔn)靜電波,當(dāng)波在等離子體中向中心運動時,TG迅速衰減,而H波只有微弱的阻尼,能夠深入等離子體中心處,因此邊緣處H波功率轉(zhuǎn)換到TG波中.但是隨著密度逐漸變大,等離子體中心的吸收有所下降.在高斯密度分布下(圖7(b)),密度低時等離子體中心處的相對吸收功率較高,但隨著密度的逐漸增大,中心處的吸收功率也慢慢變小,邊緣處變化不大明顯.

    圖7 不同密度下的徑向(z=0.2 m)吸收功率 (a)拋物面密度分布下密度對徑向相對吸收功率的影響; (b)高斯密度分布下密度對徑向相對吸收功率的影響; (c)拋物面密度分布下徑向吸收功率隨密度大小的變化曲線; (d)高斯密度分布下徑向吸收功率隨密度大小的變化曲線Fig.7.Relative power absorption in radial (z=0.2 m) directions for various density:(a) Effect of various density on radial relative absorption power under parabolic density distribution; (b) effect of various density on radial relative absorption power under Gaussian density distribution; (c) radial relative absorption power of various density under parabolic density distribution; (d) radial relative absorption power of various density under Gaussian density distribution.

    從圖7(c)和(d)可以看出,在徑向上某些密度處有一些功率峰值點,而總體的耦合程度在下降,這意味著并不是密度越大,等離子體吸收波的情況越好,而是在某一特定的密度下吸收效果最好.同樣地,等離子體在高斯分布下比拋物面分布下有著更高的相對吸收功率.因此通過對等離子體中不同密度下的吸收功率模擬,能為HPPX實驗指導(dǎo)提供良好的理論依據(jù).

    3.6 感生電磁場與電流密度

    為了進一步研究HPPX裝置中螺旋波與等離子體的耦合,圖8給出了螺旋波等離子體的電場強度、感應(yīng)磁場及電流密度的分布.從圖8(a)中可以看出,在r < 0.07 m時拋物面下的電場強度和高斯分布下的電場強度相差不大,越靠近邊緣處高斯分布下的電場值越大,這是因為兩種密度分布下等離子體邊緣處電子密度較小,邊緣處趨膚效應(yīng)較小,使得電子與離子碰撞幾率增加,相對于邊緣處的電阻增大,必然造成邊緣處負(fù)載電壓增大,進而感生電場顯著增大.總體來說高斯分布下的電場強度徑向分布較拋物面密度分布下的要大.從圖8(b)和圖8(c)中可以看出在,等離子體中心處高斯分布下的感生磁場和電流密度比拋物面下的分布大很多,而靠近邊緣處的磁場和電流逐漸變小,這是因為從密度分布位型上看,高斯分布時中心處電子密度最高,H波能夠深入中心處傳播引起電子運動速度加大進而導(dǎo)致電流密度增大,從中心至r=0.05 m處高斯分布下密度梯度變化較大,由電流密度公式J=eneve(其中e為電子電荷量,ne為電子密度,ve為電子速度)可以看出,隨著密度的降低,磁感應(yīng)強度及電流大小逐漸變小的趨勢較為明顯,拋物面分布情況下,磁感應(yīng)強度及電流變化趨勢相對平緩些,說明場能量較均勻地分布在等離子體中.

    圖8 不同密度分布下的等離子體徑向(z=0.2 m)感應(yīng)電場、感應(yīng)磁場及電流密度分布Fig.8.Radial profiles (z=0.2 m) of wave electric field,magnetic field and current density in parabolic density profile and Gaussian density profile.

    4 總 結(jié)

    本文針對最新高密度螺旋波等離子體原型實驗裝置HPPX,基于HELIC軟件計算了在拋物面密度和高斯密度分布下的天線類型、天線長度、運行頻率、靜磁場強度、軸心等離子體密度等參數(shù)對螺旋波的功率沉積的影響.計算結(jié)果表明:

    1) Half helix螺旋波天線能夠產(chǎn)生非軸對稱的射頻能量耦合,較Nagoya III型天線和Boswell型天線可激發(fā)出更高的電子密度; 通過理論分析及仿真計算,得出天線長度為0.4 m,運行頻率為13.56 MHz時等離子體與天線有著較好的能量耦合,從而得到了HPPX裝置中螺旋波系統(tǒng)發(fā)射天線的最優(yōu)設(shè)計;

    2)外加穩(wěn)恒的靜磁場時,螺旋波天線在等離子體中激發(fā)出兩種不同的性質(zhì)的波,即H波和TG波,在拋物面密度分布下,隨著磁場的增大,TG波所受的阻尼大,很難在等離子體柱中傳播,能量的吸收逐漸集中于等離子體-真空邊界處,說明較高的磁場限制了TG波向中心處的傳播; 而在高斯分布中,隨著磁場增加,TG波會被抑制,H波所受的阻尼小,可以在磁化等離子體柱內(nèi)部進行傳播,將能量耦合給等離子體; 通過對徑向積分計算可以看出,吸收功率在多個磁場處有不同的峰值功率點,且整體的耦合趨勢是在上升的;

    3)當(dāng)軸心等離子體密度增加時,在拋物面密度和高斯密度分布下,H波所受的阻尼增大,限制了H波的傳播,H波的功率轉(zhuǎn)換到TG波中; 同樣,通過積分可發(fā)現(xiàn)吸收功率在多個密度處有不同的峰值功率點,但整體的耦合趨勢卻是在下降的;

    4)總體上高斯密度分布下感應(yīng)電場、感應(yīng)磁場及電流大小較拋物面密度時要大,但由于拋物面密度梯度變化較小,等離子體中各場值和電流值較為平緩,所以能量較均勻地分布在等離子體中,此特點在之前計算功率沉積時已有體現(xiàn).

    以上這些分析可為HPPX裝置的螺旋波天線設(shè)計及相關(guān)物理實驗提供參考,同時也能為螺旋波等離子體的放電機理研究提供一定的技術(shù)支持.

    猜你喜歡
    拋物面等離子體徑向
    薯片為什么是這個形狀
    FAST照明口徑分析*
    淺探徑向連接體的圓周運動
    RN上一類Kirchhoff型方程徑向?qū)ΨQ正解的存在性
    連續(xù)磁活動對等離子體層演化的影響
    基于低溫等離子體修飾的PET/PVC浮選分離
    基于PID+前饋的3MN徑向鍛造機控制系統(tǒng)的研究
    重型機械(2020年3期)2020-08-24 08:31:40
    一類無窮下級整函數(shù)的Julia集的徑向分布
    等離子體種子處理技術(shù)介紹
    復(fù)合拋物面-漸開線聚光的悶曬式太陽能熱水器
    女性生殖器流出的白浆| 国产高清激情床上av| 视频区欧美日本亚洲| 久久久久国内视频| 久久ye,这里只有精品| 天天躁日日躁夜夜躁夜夜| 女人被躁到高潮嗷嗷叫费观| 十八禁网站网址无遮挡| 大香蕉久久成人网| www.熟女人妻精品国产| 精品少妇一区二区三区视频日本电影| 午夜福利一区二区在线看| 国产精品av久久久久免费| 最黄视频免费看| 欧美久久黑人一区二区| 亚洲国产欧美日韩在线播放| 国产在线免费精品| 亚洲人成电影免费在线| 一区二区三区国产精品乱码| 国产精品影院久久| 亚洲九九香蕉| 大片电影免费在线观看免费| 777久久人妻少妇嫩草av网站| 国产一区有黄有色的免费视频| 成在线人永久免费视频| 国产一区二区在线观看av| 免费一级毛片在线播放高清视频 | 18禁美女被吸乳视频| 国产麻豆69| 亚洲精品av麻豆狂野| 亚洲精品在线美女| 久久久久精品人妻al黑| 亚洲人成电影免费在线| 精品久久久精品久久久| 日本欧美视频一区| 中亚洲国语对白在线视频| 国产97色在线日韩免费| av不卡在线播放| 色在线成人网| 久久久久久久久久久久大奶| 中文字幕最新亚洲高清| 亚洲三区欧美一区| 精品人妻熟女毛片av久久网站| 如日韩欧美国产精品一区二区三区| 最近最新中文字幕大全电影3 | 国产在线一区二区三区精| 在线天堂中文资源库| 成年人黄色毛片网站| av一本久久久久| 不卡av一区二区三区| 天天影视国产精品| 在线看a的网站| 狂野欧美激情性xxxx| 丰满少妇做爰视频| 女人高潮潮喷娇喘18禁视频| 精品国产一区二区三区久久久樱花| 精品国产一区二区久久| 又黄又粗又硬又大视频| 国产主播在线观看一区二区| 嫁个100分男人电影在线观看| 正在播放国产对白刺激| 12—13女人毛片做爰片一| 深夜精品福利| 亚洲五月婷婷丁香| 日日摸夜夜添夜夜添小说| 女人被躁到高潮嗷嗷叫费观| 老司机影院毛片| 午夜老司机福利片| 国精品久久久久久国模美| 飞空精品影院首页| 亚洲人成77777在线视频| 久久这里只有精品19| 色综合婷婷激情| 成人国产一区最新在线观看| 国产免费现黄频在线看| 成人三级做爰电影| 天堂俺去俺来也www色官网| 亚洲成人国产一区在线观看| 欧美人与性动交α欧美精品济南到| 丁香六月欧美| 国产日韩欧美亚洲二区| 69av精品久久久久久 | 十八禁网站免费在线| 一区在线观看完整版| 高潮久久久久久久久久久不卡| 激情视频va一区二区三区| 亚洲第一青青草原| bbb黄色大片| 亚洲中文日韩欧美视频| 午夜福利欧美成人| 香蕉国产在线看| 脱女人内裤的视频| 亚洲美女黄片视频| 母亲3免费完整高清在线观看| 成人18禁在线播放| 999精品在线视频| 国产三级黄色录像| 国产成人欧美在线观看 | 日韩欧美一区视频在线观看| 在线观看免费日韩欧美大片| 成人精品一区二区免费| 十八禁高潮呻吟视频| 俄罗斯特黄特色一大片| 亚洲欧洲日产国产| 亚洲精品中文字幕在线视频| 老司机影院毛片| 可以免费在线观看a视频的电影网站| 热99国产精品久久久久久7| 午夜福利乱码中文字幕| 亚洲天堂av无毛| 日日爽夜夜爽网站| 亚洲欧美色中文字幕在线| 久久99一区二区三区| 欧美另类亚洲清纯唯美| 动漫黄色视频在线观看| 免费高清在线观看日韩| 97在线人人人人妻| 交换朋友夫妻互换小说| 精品国产乱子伦一区二区三区| 悠悠久久av| 国产一区二区三区在线臀色熟女 | 性少妇av在线| 97人妻天天添夜夜摸| 蜜桃在线观看..| 欧美成狂野欧美在线观看| 久久 成人 亚洲| av线在线观看网站| 老鸭窝网址在线观看| 国产日韩一区二区三区精品不卡| 日韩制服丝袜自拍偷拍| 国产精品久久久久久人妻精品电影 | 国产激情久久老熟女| 中文字幕av电影在线播放| 韩国精品一区二区三区| 亚洲精品美女久久久久99蜜臀| 91国产中文字幕| 日日夜夜操网爽| 十八禁高潮呻吟视频| 99国产精品一区二区蜜桃av | 亚洲精品国产区一区二| 一夜夜www| √禁漫天堂资源中文www| 中文亚洲av片在线观看爽 | 免费在线观看黄色视频的| 一本综合久久免费| 国产xxxxx性猛交| 亚洲一区二区三区欧美精品| 黄色怎么调成土黄色| 伊人久久大香线蕉亚洲五| 国内毛片毛片毛片毛片毛片| 久久人妻av系列| 国产欧美日韩一区二区三区在线| 老鸭窝网址在线观看| 亚洲欧美日韩高清在线视频 | 国产精品 欧美亚洲| 日本黄色视频三级网站网址 | 极品人妻少妇av视频| 成年人午夜在线观看视频| 国产91精品成人一区二区三区 | h视频一区二区三区| 国产av又大| 在线观看舔阴道视频| 久久精品熟女亚洲av麻豆精品| 十八禁人妻一区二区| 中文字幕av电影在线播放| 国产日韩欧美在线精品| 日韩有码中文字幕| av免费在线观看网站| 亚洲成人国产一区在线观看| av又黄又爽大尺度在线免费看| 久久国产精品影院| 天天添夜夜摸| 制服人妻中文乱码| 天天躁狠狠躁夜夜躁狠狠躁| 日本av免费视频播放| 丰满少妇做爰视频| 99国产精品免费福利视频| 亚洲国产av影院在线观看| 免费日韩欧美在线观看| 女人久久www免费人成看片| 狂野欧美激情性xxxx| 成人特级黄色片久久久久久久 | 午夜免费鲁丝| 精品一品国产午夜福利视频| 日本wwww免费看| 五月天丁香电影| 亚洲男人天堂网一区| 精品一区二区三区四区五区乱码| 日韩一区二区三区影片| 中文字幕人妻丝袜制服| xxxhd国产人妻xxx| 国产免费av片在线观看野外av| 国产精品美女特级片免费视频播放器 | 国产三级黄色录像| 成人影院久久| 99re在线观看精品视频| 菩萨蛮人人尽说江南好唐韦庄| 欧美日韩精品网址| 90打野战视频偷拍视频| 久久国产精品人妻蜜桃| 国产xxxxx性猛交| 麻豆国产av国片精品| 两性午夜刺激爽爽歪歪视频在线观看 | 999久久久精品免费观看国产| 男女床上黄色一级片免费看| 一本—道久久a久久精品蜜桃钙片| 亚洲精品成人av观看孕妇| 欧美日韩中文字幕国产精品一区二区三区 | 老司机亚洲免费影院| 男女床上黄色一级片免费看| 亚洲第一欧美日韩一区二区三区 | 中国美女看黄片| 热re99久久国产66热| 久热这里只有精品99| 国产精品秋霞免费鲁丝片| 变态另类成人亚洲欧美熟女 | 黑人巨大精品欧美一区二区mp4| 大香蕉久久网| 国产在线精品亚洲第一网站| 久热这里只有精品99| 亚洲国产av影院在线观看| 亚洲欧美日韩另类电影网站| 男女免费视频国产| 色播在线永久视频| 天堂动漫精品| 久久久国产一区二区| 亚洲全国av大片| 1024视频免费在线观看| 亚洲自偷自拍图片 自拍| 黄片小视频在线播放| 桃花免费在线播放| 午夜精品国产一区二区电影| 日本av免费视频播放| 五月天丁香电影| 亚洲中文av在线| 午夜成年电影在线免费观看| 在线观看一区二区三区激情| 精品国产一区二区久久| 每晚都被弄得嗷嗷叫到高潮| 嫁个100分男人电影在线观看| 亚洲专区字幕在线| 精品久久久久久久毛片微露脸| 成人18禁高潮啪啪吃奶动态图| 欧美av亚洲av综合av国产av| aaaaa片日本免费| 亚洲欧美激情在线| 国产成人精品无人区| 国产av国产精品国产| 两个人免费观看高清视频| 老司机午夜十八禁免费视频| 久久久久久久大尺度免费视频| 纵有疾风起免费观看全集完整版| 亚洲第一av免费看| av福利片在线| 亚洲精品中文字幕在线视频| 桃红色精品国产亚洲av| 欧美日韩精品网址| 最近最新中文字幕大全免费视频| 午夜91福利影院| 久久中文看片网| 亚洲成国产人片在线观看| 999久久久精品免费观看国产| 69av精品久久久久久 | 色精品久久人妻99蜜桃| 欧美人与性动交α欧美软件| 久久中文字幕一级| 王馨瑶露胸无遮挡在线观看| 精品少妇黑人巨大在线播放| 亚洲精品国产色婷婷电影| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲成a人片在线一区二区| 在线观看www视频免费| 久久精品国产亚洲av高清一级| 中文字幕高清在线视频| 亚洲熟妇熟女久久| 怎么达到女性高潮| 黄片小视频在线播放| 欧美另类亚洲清纯唯美| 日本欧美视频一区| 18禁裸乳无遮挡动漫免费视频| 丝袜在线中文字幕| 国产黄频视频在线观看| 久久国产精品影院| 久久天躁狠狠躁夜夜2o2o| 超碰成人久久| 人人妻人人澡人人爽人人夜夜| 精品第一国产精品| 国产精品免费一区二区三区在线 | 一级毛片女人18水好多| 亚洲,欧美精品.| 老司机福利观看| 免费看a级黄色片| 欧美日韩一级在线毛片| 亚洲精品中文字幕在线视频| 欧美成狂野欧美在线观看| 一二三四在线观看免费中文在| 久久人妻熟女aⅴ| 亚洲中文日韩欧美视频| 视频区欧美日本亚洲| 精品一区二区三区av网在线观看 | 少妇的丰满在线观看| 纯流量卡能插随身wifi吗| 国产黄频视频在线观看| 建设人人有责人人尽责人人享有的| 啪啪无遮挡十八禁网站| 亚洲色图av天堂| 一夜夜www| 丰满少妇做爰视频| 99精国产麻豆久久婷婷| 午夜91福利影院| 五月开心婷婷网| 国产精品一区二区免费欧美| 国产伦理片在线播放av一区| av天堂久久9| 一区福利在线观看| 人妻久久中文字幕网| 亚洲精品在线观看二区| 午夜精品久久久久久毛片777| 黑人巨大精品欧美一区二区mp4| 天天操日日干夜夜撸| 黑人操中国人逼视频| 美女高潮到喷水免费观看| 亚洲欧美日韩高清在线视频 | 久久精品国产综合久久久| 美女高潮到喷水免费观看| 高潮久久久久久久久久久不卡| 宅男免费午夜| 国产一区二区 视频在线| 久久亚洲精品不卡| 成人18禁在线播放| 在线永久观看黄色视频| 成人影院久久| 黄片小视频在线播放| 久久久久精品国产欧美久久久| 美女午夜性视频免费| 女人精品久久久久毛片| 十八禁高潮呻吟视频| 人人妻人人澡人人看| 中文字幕制服av| 欧美乱妇无乱码| 国产亚洲一区二区精品| 一区二区三区精品91| 777久久人妻少妇嫩草av网站| 满18在线观看网站| 老司机在亚洲福利影院| 1024视频免费在线观看| 亚洲三区欧美一区| 香蕉国产在线看| 精品一区二区三区视频在线观看免费 | 久久久久久免费高清国产稀缺| 麻豆av在线久日| 亚洲午夜精品一区,二区,三区| 天堂中文最新版在线下载| 国产亚洲精品一区二区www | 亚洲综合色网址| 亚洲成人免费av在线播放| 亚洲伊人久久精品综合| 男女午夜视频在线观看| 桃花免费在线播放| 亚洲第一av免费看| 国产精品免费大片| 久热这里只有精品99| 国产精品免费大片| 久久天堂一区二区三区四区| 久久久久久久久免费视频了| 99re6热这里在线精品视频| 国产成人av激情在线播放| 午夜福利,免费看| 欧美大码av| 国产精品国产av在线观看| 夫妻午夜视频| 欧美激情 高清一区二区三区| av一本久久久久| 久久久欧美国产精品| 午夜成年电影在线免费观看| 国产一卡二卡三卡精品| 精品国产乱子伦一区二区三区| 91精品三级在线观看| xxxhd国产人妻xxx| 成人亚洲精品一区在线观看| 午夜福利影视在线免费观看| 免费观看av网站的网址| 五月天丁香电影| 天天操日日干夜夜撸| 国产精品99久久99久久久不卡| 一级,二级,三级黄色视频| 侵犯人妻中文字幕一二三四区| 母亲3免费完整高清在线观看| 亚洲欧美精品综合一区二区三区| 国产精品 国内视频| 国产精品99久久99久久久不卡| 色精品久久人妻99蜜桃| 精品国产超薄肉色丝袜足j| 制服诱惑二区| 嫩草影视91久久| 日韩一卡2卡3卡4卡2021年| videos熟女内射| 可以免费在线观看a视频的电影网站| 国产又爽黄色视频| 国产精品久久电影中文字幕 | 午夜福利免费观看在线| av电影中文网址| 日本a在线网址| 欧美国产精品va在线观看不卡| cao死你这个sao货| 中文亚洲av片在线观看爽 | 国产精品电影一区二区三区 | av片东京热男人的天堂| 欧美av亚洲av综合av国产av| 夫妻午夜视频| 天堂动漫精品| 王馨瑶露胸无遮挡在线观看| 狠狠狠狠99中文字幕| 在线 av 中文字幕| a在线观看视频网站| 超碰97精品在线观看| 精品熟女少妇八av免费久了| 男女午夜视频在线观看| 国产精品免费大片| 亚洲男人天堂网一区| 亚洲色图综合在线观看| 露出奶头的视频| 又黄又粗又硬又大视频| 最近最新免费中文字幕在线| 蜜桃国产av成人99| aaaaa片日本免费| 丁香欧美五月| 色尼玛亚洲综合影院| 久久天堂一区二区三区四区| 国产av国产精品国产| 免费高清在线观看日韩| 成年人午夜在线观看视频| 久久精品91无色码中文字幕| 精品国产超薄肉色丝袜足j| 中文字幕人妻熟女乱码| 午夜免费成人在线视频| 国产激情久久老熟女| 亚洲成人国产一区在线观看| 国产1区2区3区精品| 久久精品成人免费网站| 亚洲精品中文字幕在线视频| 亚洲伊人久久精品综合| 日韩精品免费视频一区二区三区| 久久人人爽av亚洲精品天堂| 女性生殖器流出的白浆| 欧美国产精品va在线观看不卡| 久久天躁狠狠躁夜夜2o2o| 999久久久精品免费观看国产| 99热国产这里只有精品6| 91成人精品电影| 欧美另类亚洲清纯唯美| 男女床上黄色一级片免费看| 丰满少妇做爰视频| 俄罗斯特黄特色一大片| 精品福利观看| 69精品国产乱码久久久| 国产日韩欧美亚洲二区| 成人18禁在线播放| 波多野结衣一区麻豆| 美女视频免费永久观看网站| av片东京热男人的天堂| 人妻 亚洲 视频| 黄色怎么调成土黄色| 日韩熟女老妇一区二区性免费视频| 一本大道久久a久久精品| 交换朋友夫妻互换小说| 亚洲精品国产区一区二| 国产97色在线日韩免费| 欧美黑人精品巨大| 日日夜夜操网爽| 精品国产乱码久久久久久男人| 欧美精品亚洲一区二区| 波多野结衣一区麻豆| 女人高潮潮喷娇喘18禁视频| 美女扒开内裤让男人捅视频| 美女高潮喷水抽搐中文字幕| 欧美日韩国产mv在线观看视频| 欧美精品人与动牲交sv欧美| 欧美 日韩 精品 国产| 亚洲免费av在线视频| 日韩欧美一区视频在线观看| 99精品在免费线老司机午夜| 日本a在线网址| 黄网站色视频无遮挡免费观看| 中文字幕人妻熟女乱码| 日本欧美视频一区| 精品一区二区三区av网在线观看 | 日日夜夜操网爽| 免费女性裸体啪啪无遮挡网站| 又紧又爽又黄一区二区| 日本一区二区免费在线视频| 天天影视国产精品| 黄色a级毛片大全视频| 精品国产一区二区三区久久久樱花| 97人妻天天添夜夜摸| 免费观看a级毛片全部| 精品亚洲乱码少妇综合久久| 国产激情久久老熟女| 老汉色av国产亚洲站长工具| 久久精品亚洲精品国产色婷小说| 99国产综合亚洲精品| 少妇 在线观看| 黄色片一级片一级黄色片| 在线观看一区二区三区激情| 天堂8中文在线网| 国产av精品麻豆| 少妇 在线观看| 黄色片一级片一级黄色片| av在线播放免费不卡| 久久99热这里只频精品6学生| 日韩欧美免费精品| 亚洲精品成人av观看孕妇| 久久久精品国产亚洲av高清涩受| 两个人免费观看高清视频| 免费一级毛片在线播放高清视频 | 12—13女人毛片做爰片一| 成年女人毛片免费观看观看9 | 日韩欧美免费精品| 成人特级黄色片久久久久久久 | 亚洲精品成人av观看孕妇| 国产成人av教育| 久热这里只有精品99| 国产1区2区3区精品| 人成视频在线观看免费观看| 久久久久久久国产电影| 法律面前人人平等表现在哪些方面| 国产在线免费精品| 欧美日韩精品网址| 无人区码免费观看不卡 | 中文字幕最新亚洲高清| 黑丝袜美女国产一区| 国产精品亚洲一级av第二区| 精品一区二区三区四区五区乱码| 在线永久观看黄色视频| 午夜91福利影院| 国产成人影院久久av| 免费看a级黄色片| 成人18禁高潮啪啪吃奶动态图| 丝袜美腿诱惑在线| 热re99久久精品国产66热6| 18禁裸乳无遮挡动漫免费视频| 人妻一区二区av| 免费在线观看视频国产中文字幕亚洲| 免费人妻精品一区二区三区视频| 无遮挡黄片免费观看| 黄色a级毛片大全视频| 亚洲精品中文字幕一二三四区 | 欧美黄色淫秽网站| 中文字幕人妻丝袜一区二区| 女性被躁到高潮视频| 亚洲天堂av无毛| 精品一区二区三卡| 国产1区2区3区精品| a级毛片在线看网站| 亚洲精品在线观看二区| 成人18禁高潮啪啪吃奶动态图| 91av网站免费观看| 国产成人欧美| www.精华液| 亚洲人成77777在线视频| 亚洲欧美激情在线| 狂野欧美激情性xxxx| 捣出白浆h1v1| 欧美精品亚洲一区二区| 性少妇av在线| 日本精品一区二区三区蜜桃| a级毛片黄视频| 午夜老司机福利片| kizo精华| 99国产精品一区二区三区| 美女主播在线视频| 两个人免费观看高清视频| 黑人操中国人逼视频| 国产色视频综合| 亚洲国产毛片av蜜桃av| 高潮久久久久久久久久久不卡| cao死你这个sao货| 老熟女久久久| av天堂久久9| 18禁美女被吸乳视频| 日韩免费av在线播放| 亚洲精品国产色婷婷电影| 精品国产乱码久久久久久男人| 蜜桃国产av成人99| www.精华液| 久久中文字幕人妻熟女| 久久婷婷成人综合色麻豆| 国产激情久久老熟女| 天天添夜夜摸| 亚洲精品中文字幕在线视频| 欧美中文综合在线视频| 国产成人精品在线电影| 激情视频va一区二区三区| 人妻久久中文字幕网| 露出奶头的视频| 久久久久久免费高清国产稀缺| 熟女少妇亚洲综合色aaa.| 国产精品一区二区在线观看99| 久久久精品区二区三区| 美女扒开内裤让男人捅视频| 天天添夜夜摸| 亚洲国产欧美日韩在线播放| 丰满人妻熟妇乱又伦精品不卡| 欧美激情高清一区二区三区| 在线亚洲精品国产二区图片欧美| 中文字幕最新亚洲高清| 亚洲七黄色美女视频| 午夜精品国产一区二区电影| 制服人妻中文乱码| 大型av网站在线播放| 欧美+亚洲+日韩+国产| 人人澡人人妻人| av超薄肉色丝袜交足视频| 午夜免费成人在线视频| 久久青草综合色|