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

    激光等離子體中高能電子各向異性壓強(qiáng)的粒子模擬?

    2017-08-09 07:34:28王宬朕1董全力1劉蘋1吳奕瑩1盛政明2張杰2
    物理學(xué)報(bào) 2017年11期
    關(guān)鍵詞:高能張量動(dòng)量

    王宬朕1) 董全力1)3)? 劉蘋1) 吳奕瑩1) 盛政明2)3) 張杰2)3)

    1)(魯東大學(xué)物理與光電工程學(xué)院,煙臺(tái) 264000)

    2)(上海交通大學(xué)物理與天文學(xué)院,上海 200240)

    3)(上海交通大學(xué)IFSA協(xié)同創(chuàng)新中心,上海 200240)

    激光等離子體中高能電子各向異性壓強(qiáng)的粒子模擬?

    王宬朕1) 董全力1)3)? 劉蘋1) 吳奕瑩1) 盛政明2)3) 張杰2)3)

    1)(魯東大學(xué)物理與光電工程學(xué)院,煙臺(tái) 264000)

    2)(上海交通大學(xué)物理與天文學(xué)院,上海 200240)

    3)(上海交通大學(xué)IFSA協(xié)同創(chuàng)新中心,上海 200240)

    (2017年3月3日收到;2017年4月23日收到修改稿)

    直接驅(qū)動(dòng)慣性約束聚變(ICF)的實(shí)現(xiàn)需要對(duì)靶丸進(jìn)行嚴(yán)格的對(duì)稱壓縮,以達(dá)到自持熱核反應(yīng)(點(diǎn)火)所需的條件.快點(diǎn)火方案的應(yīng)用降低了對(duì)靶丸壓縮對(duì)稱性以及驅(qū)動(dòng)能量的要求,但壓縮及核反應(yīng)過(guò)程中良好的靶丸對(duì)稱性無(wú)疑有助于核反應(yīng)增益的提高.本文研究了快點(diǎn)火方案中高能電子注入高密等離子體后導(dǎo)致的各向異性電子的壓強(qiáng)張量.這一現(xiàn)象存在于ICF快點(diǎn)火方案中的高能電子束“點(diǎn)火”及核反應(yīng)階段.鑒于高能電子加熱離子過(guò)程以及靶丸核反應(yīng)自持燃燒過(guò)程的時(shí)間較長(zhǎng),高密靶核會(huì)由于超高的各向異性壓強(qiáng)的作用破壞高密靶丸的對(duì)稱性,降低核燃料密度,進(jìn)而降低了核燃料燃燒效率以及核反應(yīng)增益.

    各向異性壓強(qiáng),粒子模擬,快點(diǎn)火

    1 引言

    慣性約束聚變中心點(diǎn)火方案是指在非常短的時(shí)間尺度內(nèi)利用高能量高功率密度的脈沖激光均勻照射氘氚燃料靶丸[1],該過(guò)程產(chǎn)生的靶面物質(zhì)等離子體外向噴射的反沖作用將靶內(nèi)物質(zhì)快速地壓縮至超高密度(約為1024cm?3)、超高溫度(約為108K)和超高壓強(qiáng)(約為1010atm)狀態(tài)[2],同時(shí)利用物質(zhì)自身的慣性將該狀態(tài)維持一定的時(shí)間(約為10?9s)[3],使之達(dá)到自持熱核反應(yīng)所需條件,業(yè)內(nèi)稱之為“點(diǎn)火”[4].為了能夠?qū)型桧樌貕嚎s至以上狀態(tài),需要在壓縮過(guò)程中嚴(yán)格保持靶丸的對(duì)稱性.影響靶丸對(duì)稱壓縮的初始性條件有很多,包括驅(qū)動(dòng)光源的對(duì)稱性、作為熱核燃料靶丸的球形度、壁厚的均勻性、表面粗糙度、殼層密度分布、殼層材料的組成和摻雜水平以及靶丸裝配偏心程度等.在慣性約束聚變(ICF)的內(nèi)爆實(shí)驗(yàn)中,通常對(duì)稱布置多路高能量、高功率密度的脈沖激光[5,6],同時(shí)照射到球形氘氚靶丸上.靶壓縮對(duì)稱性對(duì)激光輻照的空間分布極為敏感.驅(qū)動(dòng)光源光束強(qiáng)度的不均勻性易激發(fā)壓縮靶面的流體力學(xué)穩(wěn)定性,如Richtmyer-Meshkov(RM)不穩(wěn)定性等[7,8].為實(shí)現(xiàn)高效的對(duì)稱內(nèi)爆和高增益,其驅(qū)動(dòng)輻照不均勻性必須控制到1%以下[2,5].靶丸表面的粗糙度以及小尺度的靶面變形也會(huì)導(dǎo)致內(nèi)爆壓縮過(guò)程中流體力學(xué)不穩(wěn)定性的出現(xiàn)與增長(zhǎng)[5,9?11].瑞利-泰勒不穩(wěn)定性是影響激光內(nèi)爆對(duì)稱壓縮靶丸的重要問(wèn)題之一,容易發(fā)生在內(nèi)爆壓縮的加速階段和減速階段[12].它能夠使靶丸初始的不均勻特征在壓縮過(guò)程中呈指數(shù)增長(zhǎng),從而極大地降低靶丸結(jié)構(gòu)的對(duì)稱和完整,破壞性很大[13,14].在美國(guó)國(guó)家點(diǎn)火裝置上進(jìn)行的慣性約束聚變氘氚內(nèi)爆實(shí)驗(yàn)中,測(cè)量到的中子產(chǎn)額和熱斑壓強(qiáng)明顯低于模擬數(shù)值,其原因就是因?yàn)榘型璧拿婷芏炔痪鶆蛐云茐牧藟嚎s殼層的穩(wěn)定性,并且極大地影響了內(nèi)爆動(dòng)能轉(zhuǎn)化為熱斑內(nèi)能[15].

    針對(duì)中心點(diǎn)火方案中出現(xiàn)的流體力學(xué)不穩(wěn)定性導(dǎo)致的靶丸壓縮問(wèn)題,業(yè)內(nèi)提出了將“壓縮”和“點(diǎn)火”兩個(gè)過(guò)程分開(kāi)進(jìn)行的“快點(diǎn)火”技術(shù)方案[16].在壓縮階段利用納秒長(zhǎng)脈沖激光束對(duì)靶丸進(jìn)行高度對(duì)稱壓縮,當(dāng)聚變?nèi)剂媳痪鶆驂嚎s到最大密度時(shí),將一束超短脈沖強(qiáng)激光聚焦在靶丸表面(光強(qiáng)>1020W/cm2),極高的有質(zhì)動(dòng)力在靶丸表面的等離子體臨界密度面上“打洞”,并將臨界密度面向靶芯的高密核處擠壓.在這個(gè)過(guò)程中產(chǎn)生的大量的能量高達(dá)MeV的電子穿透臨界密度面射入高密核,并將能量傳輸給核燃料離子,使其溫度迅速升至點(diǎn)火所要求的5—10 keV的高溫,實(shí)現(xiàn)對(duì)靶丸的快速點(diǎn)火.由于ICF快點(diǎn)火方案降低了對(duì)靶設(shè)計(jì)、光束質(zhì)量以及輻照均勻性的嚴(yán)格要求,因而很快成為ICF研究中的一個(gè)重要方向.科學(xué)家們對(duì)快點(diǎn)火方案涉及的諸多物理和技術(shù)問(wèn)題進(jìn)行了理論與分解實(shí)驗(yàn)研究,并取得了較好的研究結(jié)果.比如,針對(duì)“點(diǎn)火”過(guò)程需要的超短脈沖強(qiáng)激光及其所產(chǎn)生的高能電子在高密核外圍等離子體中存在各種傳輸不穩(wěn)定性,通過(guò)采用錐球靶,吳鳳娟等[17]對(duì)強(qiáng)激光與錐型三明治結(jié)構(gòu)靶相互作用快電子束的產(chǎn)生與傳輸進(jìn)行了粒子模擬,并且在相同激光參數(shù)下與其他錐型靶的作用結(jié)果進(jìn)行了比較.采用錐型靶的目的就是讓點(diǎn)火激光脈沖通過(guò)錐殼屏蔽的空間直接與壓縮后的高密度靶核相互作用,產(chǎn)生高能電子實(shí)施點(diǎn)火,避免了激光以及高能電子在較低密等離子體區(qū)的傳輸不穩(wěn)定性,提高了靶丸核反應(yīng)效率與中子產(chǎn)額[18].2005年,谷渝秋等[19]利用3 TW飛秒激光器與固體靶相互作用,研究了超熱電子構(gòu)成各向異性的能量分布,而針對(duì)快點(diǎn)火中相對(duì)論電子束能量沉積的問(wèn)題,周滄濤等[20,21]在物理建模、程序研制以及數(shù)值模擬研究等方面做了大量工作,同時(shí)在激光驅(qū)動(dòng)強(qiáng)流電子束的產(chǎn)生和控制方面做了細(xì)致的研究.

    快點(diǎn)火方案中,高能電子束注入到高密靶核里,需要有效地將其能量淀積給燃料離子,同時(shí)也給高密靶丸帶來(lái)了對(duì)稱性問(wèn)題.由于高能電子僅從一個(gè)方向注入,其所造成的壓強(qiáng)張量必然是各向異性的.在電子將能量轉(zhuǎn)移給燃料離子以及核反應(yīng)持續(xù)發(fā)生這一較長(zhǎng)時(shí)間內(nèi),電子各向異性壓強(qiáng)張量預(yù)期將通過(guò)影響核燃料的空間分布,進(jìn)而影響到核燃料的燃燒效率及核反應(yīng)過(guò)程的增益.但之前的研究顯然對(duì)此效應(yīng)影響點(diǎn)火過(guò)程及結(jié)果缺乏評(píng)估與研究[22].本文利用一維粒子模擬(PIC),研究了激光等離子體相互作用過(guò)程中產(chǎn)生的高能電子[23]及其對(duì)等離子體內(nèi)部產(chǎn)生的超高壓強(qiáng)的各向異性特征.由于目前人們通過(guò)在快點(diǎn)火方案壓縮階段結(jié)束后引入外加磁場(chǎng)來(lái)引導(dǎo)或者束縛高能電子對(duì)高密壓縮核的注入及能量淀積[24],本文亦研究了背景磁場(chǎng)對(duì)各向異性電子壓強(qiáng)張量的影響.

    2 基本方程

    壓強(qiáng)張量的定義式為[25?30]

    (2)式中的?···?表示對(duì)其中的物理量進(jìn)行加權(quán)平均.所以靶物質(zhì)內(nèi)某點(diǎn)處的壓強(qiáng)狀態(tài)在直角坐標(biāo)系下可以表示為張量:

    壓強(qiáng)張量有6個(gè)獨(dú)立分量.物體在壓強(qiáng)張量作用下的形變分為體積變化和形狀變化兩部分,其中體積變化取決于各向同性壓強(qiáng),形狀變化取決于各向異性壓強(qiáng).我們將壓強(qiáng)張量Pe進(jìn)行分解,取三個(gè)主軸Pexx,Peyy,Pezz中數(shù)值最小的分量記為標(biāo)量壓強(qiáng)Pem,構(gòu)成各向同性壓強(qiáng)張量;用壓強(qiáng)張量減去各向同性壓強(qiáng)張量之后得到的就是導(dǎo)致形變的各向異性壓強(qiáng)張量,記作πeij.于是有

    若系統(tǒng)處在熱平衡狀態(tài),人們可以利用理想氣體狀態(tài)方程計(jì)算各向同性壓強(qiáng):

    其中Pe為電子的標(biāo)量壓強(qiáng),ne為電子數(shù)密度,kB=1.38×10?23J/K為玻爾茲曼常數(shù),Te=mev2/(3kB)為電子溫度.多粒子系統(tǒng)通過(guò)粒子間的相互作用實(shí)現(xiàn)能量和動(dòng)量的交換,最終達(dá)到熱平衡態(tài)[29,31?35].可以大致估算一下激光等離子體需要經(jīng)過(guò)多長(zhǎng)時(shí)間能達(dá)到平衡.電子的碰撞頻率νee=2.91×10?6nelnΛT?3/2es?1,其中ne=1027m?3,lnΛ=10,Te=108K,可計(jì)算出電子的碰撞頻率大約為1010s?1.若我們?cè)O(shè)置的激光脈寬為1 ns,則電子碰撞一次大概需要經(jīng)歷1/10個(gè)激光脈沖的時(shí)間,即在激光脈沖時(shí)間內(nèi),電子可運(yùn)動(dòng)至平衡狀態(tài).但實(shí)驗(yàn)表明較高能量的激光可能會(huì)產(chǎn)生部分高能電子.這部分電子可占總數(shù)的1/10,其溫度可達(dá)幾百keV,熱化過(guò)程比較慢,且其能量空間分布亦不均勻,必須直接利用壓強(qiáng)張量定義式對(duì)各向異性壓強(qiáng)張量πe進(jìn)行計(jì)算.在本文中我們利用一維PIC模擬產(chǎn)生各向異性分布的高能電子并探討其壓強(qiáng)張量Pe的特征.

    3 PIC模擬結(jié)果

    3.1 PIC模擬參數(shù)設(shè)置

    模擬中我們對(duì)激光場(chǎng)的部分參數(shù)進(jìn)行了如下設(shè)置:無(wú)量綱化的電場(chǎng)振幅峰值為2.7;激光波長(zhǎng)為1μm;激光場(chǎng)的功率密度I=1.0×1019W/cm2;脈沖包絡(luò)為高斯型,脈寬τ=30τ0≈100 fs,τ0為激光周期;激光的入射角度為45°;偏振方向?yàn)镻偏振;初始的電子溫度設(shè)置為2.0 eV.圖1(a)中紅線給出了靶物質(zhì)的初始電子密度分布,橫坐標(biāo)是以波長(zhǎng)λ為單位的空間長(zhǎng)度L,右側(cè)縱坐標(biāo)是以臨界密度nc歸一化的電子密度.橫坐標(biāo)8λ—16λ是電子密度分布為e指數(shù)增長(zhǎng)區(qū)域,16λ—20λ是電子密度為5nc且均勻分布的高密區(qū)域,0—8λ和20λ—28λ是密度為零的真空區(qū)域,臨界密度nc的位置在13.5λ處.考慮到激光等離子體中存在極強(qiáng)的自生磁場(chǎng),而且在ICF快點(diǎn)火方案壓縮階段結(jié)束后,可利用外加磁場(chǎng)進(jìn)行高能電子注入引導(dǎo)或者約束提高其能量向燃料離子的淀積效率.我們還將討論背景磁場(chǎng)對(duì)各向異性壓強(qiáng)特征的影響.

    3.2 無(wú)背景磁場(chǎng)條件下各向異性壓強(qiáng)的特征

    圖1(網(wǎng)刊彩色)(a)電子密度分布;(b)電子動(dòng)量分布Fig.1.(color online)(a)Electron density distribution;(b)electron momentum distribution.

    圖1 是在激光與靶物質(zhì)相互作用過(guò)程中,電子密度(圖1(a))以及電子動(dòng)量(圖1(b))的時(shí)空分布.由于在臨界密度處激光場(chǎng)強(qiáng)梯度最大,有質(zhì)動(dòng)力一方面導(dǎo)致電子密度凹坑的形成,另一方面將電子加速至超高能量.當(dāng)這些高能電子進(jìn)入到高密區(qū)域,激發(fā)了高密區(qū)域等離子體的振蕩,如圖1(a)中黑線表示的75τ0時(shí)刻的電子密度分布(右側(cè)縱坐標(biāo))所示.圖1(a)中的電子密度的時(shí)空演化(左側(cè)縱坐標(biāo),時(shí)間t)更加清晰地體現(xiàn)出這一點(diǎn).從圖1(b)的電子動(dòng)量的時(shí)空分中可以看出,電子動(dòng)量空間分布具有一定的周期性,這顯示了高能電子主要是由激光有質(zhì)動(dòng)力加速而成,這些高能電子進(jìn)入到高密區(qū)域后保留了這一群特征.同時(shí)隨著激光作用時(shí)間的延續(xù),高密區(qū)域中高能電子的動(dòng)量以及數(shù)目明顯地增大、增多,其加速機(jī)制獲取的周期性特征在100τ0后逐漸消失.

    圖2(網(wǎng)刊彩色)(a)各向同性壓強(qiáng)Pe;(b)各向同性壓強(qiáng)pem;(c)沿X方向上電子動(dòng)量pe;(d)電子的壓強(qiáng)張量PexxFig.2.(color online)(a)Isotropic pressure Pe;(b)isotropic pressure Pem;(c)electron momentum along the X direction pe;(d)pressure tensor of electron Pexx.

    圖2 (a)和圖2(b)是在16λ—19λ的高密區(qū)域中利用(5)式求得的各向同性壓強(qiáng)Pe和利用壓強(qiáng)張量分量求得的各向同性壓強(qiáng)Pem的時(shí)空分布圖.從圖2(a)中可以看出,各向同性壓強(qiáng)隨時(shí)間變化得比較平緩,在空間分布上也呈現(xiàn)出較為均勻的態(tài)勢(shì).但圖2(b)呈現(xiàn)出與圖2(a)不同的趨勢(shì),Pem的時(shí)空演化有較大振蕩,而在數(shù)值上也明顯比Pe要小1/20之多.從圖2(c)在X方向的電子動(dòng)量的時(shí)空分布可知,在50τ0時(shí)刻以后,高密區(qū)域才開(kāi)始有高能電子進(jìn)入,與此同時(shí)圖2(d)中電子的壓強(qiáng)張量分量Pexx也開(kāi)始逐步增大,并且可以發(fā)現(xiàn)壓強(qiáng)張量分量Pexx隨著X方向的電子動(dòng)量pe的變化而變化,二者幾乎是同步進(jìn)行的,表明在高密區(qū)域的壓強(qiáng)張量是由于電子的進(jìn)入而導(dǎo)致的.壓強(qiáng)張量是各向同性和各向異性壓強(qiáng)之和,當(dāng)電子達(dá)到熱平衡狀態(tài)之后,則表現(xiàn)出各向同性壓強(qiáng)的特征及性質(zhì),但高能電子卻由于短時(shí)間內(nèi)無(wú)法達(dá)到熱平衡態(tài),主要導(dǎo)致各向異性壓強(qiáng).

    圖3(a)給出了65τ0時(shí)刻各向異性壓強(qiáng)張量πe的6個(gè)獨(dú)立分量,可見(jiàn)6個(gè)壓強(qiáng)分量無(wú)論是在大小、方向,還是在變化趨勢(shì)上均有較大差異.最明顯的差異是各個(gè)方向上的壓強(qiáng)分量大小不盡相同.πexx是其中最大的各向異性壓強(qiáng)分量值,而πezz值最小.πexx,πeyy和πexy明顯地以不同的振動(dòng)幅度隨著空間變化而振蕩,其原因是電子被激光有質(zhì)動(dòng)力加速的過(guò)程中,產(chǎn)生的電子動(dòng)量具有一定的空間周期性.所以高密區(qū)域內(nèi)各向異性壓強(qiáng)的產(chǎn)生歸因于空間分布不均勻的高能電子群.

    圖3(網(wǎng)刊彩色)(a)65τ0時(shí)刻各向異性壓強(qiáng)張量πe;(b)75τ0(紅色)和85τ0(藍(lán)色)時(shí)刻的πexx,虛線是低能電子的貢獻(xiàn),點(diǎn)線是高能電子的貢獻(xiàn),實(shí)線是總的電子貢獻(xiàn)值;(c)低能電子對(duì)πexx的貢獻(xiàn);(d)高能電子對(duì)πexx的貢獻(xiàn)Fig.3.(color online)(a)Anisotropic pressure πeon the moment of 65τ0;(b)πexxon the moment of 75τ0(red)and 85τ0(blue),the contribution of thermal electronsis dashed line,the contribution of hot electrons is dotted line,the contribution of electron is real line;(c)the contribution of thermal electrons to πexx;(d)the contribution of hot electrons to πexx.

    圖3 (b)是在75τ0(紅色線條)和85τ0(藍(lán)色線條)兩個(gè)時(shí)刻的各向異性壓強(qiáng)πexx隨空間的變化情況.圖中虛線表示低能電子對(duì)πexx的貢獻(xiàn),點(diǎn)線表示高能電子對(duì)πexx的貢獻(xiàn),實(shí)線則表示高密區(qū)域內(nèi)總的各向異性壓強(qiáng)πexx值.從圖3(b)可以更清晰地看出,高能電子對(duì)πexx的貢獻(xiàn)遠(yuǎn)超出低能電子,并且在總的各向異性壓強(qiáng)πexx中占據(jù)主導(dǎo)地位.圖3(c)和圖3(d)分別是高密區(qū)域中高能電子和低能電子對(duì)各向異性壓強(qiáng)πexx貢獻(xiàn)的時(shí)空演化圖.隨著作用時(shí)間的延長(zhǎng),低能電子對(duì)各向異性壓強(qiáng)的貢獻(xiàn)逐漸減少,最終趨于一個(gè)較小且穩(wěn)定的數(shù)值范圍.從數(shù)目上看,高密區(qū)域中的高能電子數(shù)目遠(yuǎn)遠(yuǎn)小于低能電子.在65τ0時(shí)刻以前高能電子所占比例不到總電子數(shù)的5%,75τ0時(shí)刻為止所占比例不超過(guò)15%,而85τ0時(shí)刻及以后高能電子所占比例基本趨于一個(gè)穩(wěn)定值,大約為總電子數(shù)的20%.即便數(shù)目比例如此小,高能電子卻在各向異性壓強(qiáng)中占據(jù)了大約80%—90%的貢獻(xiàn)值,是各向異性壓強(qiáng)的主要貢獻(xiàn)者.這說(shuō)明當(dāng)激光等離子體相互作用過(guò)程中,產(chǎn)生的大量低能電子在經(jīng)歷一段時(shí)間之后會(huì)趨于熱平衡狀態(tài),表現(xiàn)出較強(qiáng)的各向同性,而高能電子則可能由于熱化過(guò)程較慢,并未達(dá)到一個(gè)熱平衡態(tài),導(dǎo)致其各向異性較為突出.在直接驅(qū)動(dòng)慣性約束聚變實(shí)驗(yàn)中靶丸如此巨大的各向異性超高壓,主要是由激光與靶丸作用過(guò)程中產(chǎn)生的各向異性的高能電子造成的.同時(shí)我們也注意到85τ0時(shí)刻的πexx值比75τ0時(shí)刻的πexx值要大得多,并且高能電子的貢獻(xiàn)比值也進(jìn)一步增大了,這是因?yàn)檫M(jìn)入高密區(qū)域內(nèi)的高能電子的數(shù)目以及動(dòng)量都有所增大.

    3.3 背景磁場(chǎng)對(duì)各向異性壓強(qiáng)的影響

    激光等離子體內(nèi)自生磁場(chǎng)是相當(dāng)復(fù)雜的,本文僅討論當(dāng)加入大小0.1 T、垂直激光入射方向(模擬設(shè)置為Z軸正方向)的背景磁場(chǎng)的情形下高密區(qū)域內(nèi)各向異性壓強(qiáng)的特征.

    圖4(a)是60τ0時(shí)刻等離子體高密區(qū)域電子相空間圖,給出了背景磁場(chǎng)對(duì)高能電子動(dòng)量分布的影響.選取相同時(shí)刻(60τ0)電子在高密區(qū)域沿X方向的動(dòng)量分布,紅色點(diǎn)表示無(wú)背景磁場(chǎng)時(shí)電子的動(dòng)量,藍(lán)色點(diǎn)則是存在背景磁場(chǎng)時(shí)的電子動(dòng)量.可以看出,在相同作用時(shí)間的情況下,無(wú)背景磁場(chǎng)約束的高能電子保持了激光有質(zhì)動(dòng)力加速特征,在整個(gè)區(qū)域內(nèi)的分布具有一定的空間周期性;而背景磁場(chǎng)的洛倫茲力卻束縛了高能電子,將其阻滯在高密區(qū)域的前端,僅使其中很少一部分能進(jìn)入深層區(qū)域.

    圖4 (網(wǎng)刊彩色)(a)60τ0時(shí)刻電子動(dòng)量分布;(b)65τ0時(shí)刻各向異性壓強(qiáng)張量πe;(c)低能電子對(duì)πexx的貢獻(xiàn);(d)高能電子對(duì)πexx的貢獻(xiàn)Fig.4.(color online)(a)Electron momentum distribution on the moment of 60τ0;(b)anisotropic pressure πeon the moment of 65τ0;(c)the contribution of thermal electrons to πexx;(d)the contribution of hot electrons to πexx.

    圖4(b)是65τ0時(shí)刻存在背景磁場(chǎng)時(shí)各向異性壓強(qiáng)張量πe的6個(gè)分量.從6個(gè)壓強(qiáng)分量的整體變化來(lái)看,振蕩幅度相差較大,而且各個(gè)壓強(qiáng)分量值主要集中在高密區(qū)域的前端,深層區(qū)域內(nèi)的各個(gè)壓強(qiáng)分量值都比較小,基本都是此區(qū)域內(nèi)的低能電子的貢獻(xiàn).與無(wú)背景磁場(chǎng)條件下的各向異性壓強(qiáng)張量分量圖3(a)進(jìn)行比較可以看出,在正X方向上的壓強(qiáng)πexx比無(wú)背景磁場(chǎng)時(shí)的數(shù)值小5—7倍之多,且隨著空間改變整體呈減小的趨勢(shì);在正Y方向上的壓強(qiáng)πeyy也隨著空間改變整體呈減小的趨勢(shì);但與背景磁場(chǎng)平行的正Z方向上的壓強(qiáng)πezz其大小及方向與無(wú)背景磁場(chǎng)條件下的πezz相差不大明顯.兩圖的比較顯示,在Z方向上的背景磁場(chǎng)對(duì)各向異性壓強(qiáng)分量πexx,πeyy和πexy的影響較為明顯,而對(duì)πezz的影響較小.圖4(c)和圖4(d)是存在背景磁場(chǎng)時(shí)高密區(qū)域內(nèi)低能電子和高能電子分別對(duì)各向異性壓強(qiáng)πexx的貢獻(xiàn)情況,可以看出,低能電子對(duì)πexx的貢獻(xiàn)比高能電子的貢獻(xiàn)小了近10倍之多,但隨著作用時(shí)間的延長(zhǎng),其貢獻(xiàn)逐漸減小,并趨于一個(gè)較穩(wěn)定的值;而高能電子的貢獻(xiàn)則主要集中在高密區(qū)域的前端,并隨著空間的改變壓強(qiáng)值逐漸減小.

    由高能電子引起的高密區(qū)域內(nèi)的各向異性壓強(qiáng)會(huì)導(dǎo)致對(duì)稱壓縮的高密靶丸在超高壓強(qiáng)的作用下體積膨脹增大,進(jìn)而引起靶丸密度的降低.我們假設(shè)在高能電子注入之初為等壓膨脹,當(dāng)5.0×108atm的壓強(qiáng)作用在半徑為100μm,密度達(dá)到300 g/cm3的高密靶核時(shí),靶丸的半徑會(huì)在2 ps內(nèi)迅速?gòu)脑瓉?lái)的100μm增大到150μm,其密度也會(huì)由初始的300 g/cm3減小到100 g/cm3以下,如此低的靶丸密度將會(huì)明顯地降低核燃料燃燒效率以及核反應(yīng)增益.

    4 總結(jié)

    本文利用一維PIC模擬程序,研究了激光與等離子體相互作用過(guò)程中產(chǎn)生的高能電子傳輸?shù)礁呙芏葏^(qū)域后導(dǎo)致的壓強(qiáng)張量特征改變.研究表明,等離子體高密區(qū)域內(nèi)各向異性壓強(qiáng)主要是由高能電子造成的.本文也研究了背景磁場(chǎng)對(duì)各向異性壓強(qiáng)的影響,研究結(jié)果發(fā)現(xiàn),在垂直于背景磁場(chǎng)方向上,因電子被磁場(chǎng)束縛,其各向異性壓強(qiáng)張量分量的數(shù)值要大于平行于背景磁場(chǎng)方向上的壓強(qiáng)分量;垂直方向上壓強(qiáng)分量的空間變化也比平行分量要?jiǎng)×?在快點(diǎn)火方案中,在靶丸壓縮結(jié)束后,高能電子向燃料離子的能量傳遞過(guò)程以及隨后的核反應(yīng)過(guò)程非一蹴而就,而是需要相對(duì)較長(zhǎng)的時(shí)間.在這段時(shí)間內(nèi),流體力學(xué)效應(yīng)比較明顯,高能電子的各向異性壓強(qiáng)勢(shì)必會(huì)影響到燃料的空間分布,進(jìn)而影響燃燒效率以及中子的產(chǎn)額.本文對(duì)高能電子各向異性壓強(qiáng)張量的研究意在提醒人們關(guān)注這一點(diǎn).

    [1]Craxton R S,Anderson K S,Boehly T R,Goncharov V N,Harding D R,Knauer J P,McCrory R L,McKenty P W,Meyerhofer D D,Myatt J F,Schmitt A J,Sethian J D,Short R W,Skupsky S,Theobald W,Kruer W L,Tanaka K,Betti R,Collins T J B,Delettrez J A,Hu S X,Marozas J A,Maximov A V,Michel D T,Radha P B,Regan S P,Sangster T C,Seka W,Solodov A A,Soures J M,Stoeckl C,Zuegel J D 2015 Phys.Plasmas 22 110501

    [2]Lindl J 1995 Phys.Plasmas 2 3933

    [3]Drake R P 2006 High-Energy-Density Physics:Fundamentals,Inertial Fusion,and Experimental Astrophysics(1st Ed.)(Berlin:Springer Science&Business Media)pp392–419

    [4]McCrory R L,Meyerhofer D D,Betti R,Craxton R S,Delettrez J A,Edgell D H,Glebov V Yu,Goncharov V N,Harding D R,Jacobs-Perkins D W,Knauer J P,Marshall F J,McKenty P W,Radha P B,Regan S P,Sangster T C,Seka W,Short R W,Skupsky S,Smalyuk V A,Soures J M,Stoeckl C,Yaakobi B,Shvarts D,Frenje J A,Li C K,Petrasso R D,Séguin F H 2008 Phys.Plasmas 15 055503

    [5]Rosen M D 1999 Phys.Plasmas 6 1690

    [6]Bodner S E,Colombant D G,Gardner J H,Lehmberg R H,Obenschain S P,Phillips L,Schmitt A J,Sethian J D 1998 Phys.Plasmas 5 1901

    [7]Sharp D H 1984 Physica D 12 3IN111

    [8]Brouillette M 2002 Annu.Rev.Fluid Mech.34 445

    [9]Wesson J,Campbell D J 2011 Tokamaks(4th Ed.)(Oxford:Oxford University Press)pp356–358

    [10]Li C K,Séguin F H,Frenje J A,Petrasso R D,Delettrez JA,McKenty P W,Sangster T C,Keck R L,Soures J M,Marshall F J,Meyerhofer D D,Goncharov V N,Knauer J P,Radha P B,Regan S P,Seka W 2004 Phys.Rev.Lett.92 205001

    [11]Shigemori K,Azechi H,Nakai M,Honda M,Meguro K,Miyanaga N,Takabe H,Mima K 1997 Phys.Rev.Lett.78 250

    [12]Honda M,Mima K,Shigemori K,Nakai M,Azechi H,Nishiguchi A 1999 Fusion Eng.Des.44 205

    [13]Lindl J D,McCrory R L,Campbell E M 1992 Phys.Today 45 32

    [14]Wouchuk J G 2001 Phys.Rev.E 63 056303

    [15]Gu J F,Dai Z S,Fan Z F,Zou S Y,Ye W H,Pei W B,Zhu S P 2014 Phys.Plasmas 21 012704

    [16]Tabak M,Hammer J,Glinsky M E,Kruer W L,Wilks S C,Woodworth J,Campbell E M,Perry M D 1994 Phys.Plasmas 1 1626

    [17]Wu F J,Zhou W M,Shan L Q,Li F,Liu D X,Zhang Z M,Li B Y,Bi B,Wu B,Wang W W,Zhang F,Gu Y Q,Zhang B H 2014 Acta Phys.Sin.63 94101(in Chinese)[吳鳳娟,周維民,單連強(qiáng),李芳,劉東曉,張智猛,李博原,畢碧,伍波,王為武,張鋒,谷渝秋,張保漢2014物理學(xué)報(bào)63 94101]

    [18]Kodama R,Norreys P A,Mima K,Dangor A E,Evans R G,Fujita H,Kitagawa Y,Krushelnick K,Miyakoshi T,Miyanaga N,Norimatsu T,Rose S J,Shozaki T,Shigemori K,Sunahara A,Tampo M,Tanaka K A,Toyama Y,Yamanaka T,Zepf M 2001 Nature 412 798

    [19]Gu Y Q,Cai D F,Zheng Z J,Yang X D,Zhou W M,Jiao C Y,Chen H,Wen T S,Chunyu S T 2005 Acta Phys.Sin.54 186(in Chinese)[谷渝秋,蔡達(dá)鋒,鄭志堅(jiān),楊向東,周維民,焦春曄,陳豪,溫天舒,淳于書泰2005物理學(xué)報(bào)54 186]

    [20]Wu S Z,Zhang H,Zhou C T,Wu J F,Cai H B,Cao L H,He M Q,Zhu S P,He X T 2015 High Power Laser and Particle Beams 27 77(in Chinese)[吳思忠,張華,周滄濤,吳俊峰,蔡洪波,曹莉華,何民卿,朱少平,賀賢土2015強(qiáng)激光與粒子束27 77]

    [21]Cai H B,Zhou C T,Jia Q,Wu S Z,He M Q,Cao L H,Chen M,Zhang H,Liu J,Zhu S P,He X T 2015 High Power Laser and Particle Beams 27 27032001(in Chinese)[蔡洪波,周滄濤,賈青,吳思忠,何民卿,曹莉華,陳默,張華,劉杰,朱少平,賀賢土2015強(qiáng)激光與粒子束27 27032001]

    [22]Zhang J 1999 Physics 28 0(in Chinese)[張杰1999物理28 0]

    [23]Cassak P A,Baylor R N,Fermo R L,Beidler M T,Shay M A,Swisdak M,Drake J F,Karimabadi H 2015 Phys.Plasmas 22 020705

    [24]Wang W M,Gibbon P,Sheng Z M,Li Y T 2015 Phys.Rev.Lett.114 015001

    [25]Divin A,Markidis S,Lapenta G,Semenov V S,Erkaev N V,Biernat H K 2010 Phys.Plasmas 17 122102

    [26]Hoshino M 2005 J.Geophys.Res.110 A10215

    [27]Liu C,Fox W,Bhattacharjee A 2015 Phys.Plasmas 22 053302

    [28]Wan W G,Lapenta G 2008 Phys.Rev.Lett.101 015001

    [29]Yin L,Winske D,Gary S P,Birn J 2001 J.Geophys.Res.106 10761

    [30]Wang L,Hakim A H,Bhattacharjee A,Germaschewski K 2015 Phys.Plasmas 22 012108

    [31]Mottez F 2004 Ann.Geophys.22 3033

    [32]Heinz H,Paul W,Binder K 2005 Phys.Rev.E 72 066704

    [33]Cai H S,Li D 2009 Phys.Plasmas 16 052107

    [34]Le A,Daughton W,Karimabadi H,Egedal J 2016 Phys.Plasmas 23 032114

    [35]Yin L,Winske D 2003 Phys.Plasmas 10 1595

    PACS:52.27.–h,52.38.–r,52.65.–y,52.65.RrDOI:10.7498/aps.66.115203

    Particle simulation study on anisotropic pressure of electrons in laser-produced plasma interaction?

    Wang Cheng-Zhen1)Dong Quan-Li1)3)?Liu Ping1)Wu Yi-Ying1)Sheng Zheng-Ming2)3)Zhang Jie2)3)
    1)(School of Physics and Optoelectronic Engineerings,Ludong University,Yantai 264000,China)
    2)(Department of Physics and Astronomy,Shanghai Jiao Tong University,Shanghai 200240,China)
    3)(IFSA Collaborative Innovation Center,Shanghai Jiao Tong University,Shanghai 200240,China)

    3 March 2017;revised manuscript

    23 April 2017)

    Direct-drive inertial con fi nement fusion(ICF)requires a symmetric compression of the fuel target to achieve physical conditions for the ignition.The fast ignition scheme reduces the symmetry requirements for the target compression and the necessary driving energy,but symmetrically compressed target will certainly help improve the efficiency of the nuclear fuel burning.In this paper,with the particle-in-cell(PIC)simulation method,characteristics of the anisotropic pressure tensor of hot electrons are reported for the ultra intense laser pulse interaction with over dense plasmas,which mimics the scenario of the last stage when hot electrons are utilized to ignite the compressed fuel core in the ICF fast ignition scheme.A large number of hot electrons can stimulate pressure oscillations in the high density plasma.As the component parallel to the electron velocity dominates the pressure tensor,the electron density distribution perturbation propagates rapidly in this direction.In order to keep those hot electrons in the high density fuel plasma core for a period long enough for them to deposit energy and momentum,a magnetic fi eld perpendicular to the electron velocity is used.The PIC simulation results indicate that the hot electrons can be trapped by the magnetic fi eld,and the components of the anisotropic pressure tensor related to the parallel direction are signi fi cantly a ff ected,thereby producing a high peak near the incidence surface.Since it is a relatively long process for the energy transfer from electrons to fuel ions and the nuclear interaction to be completed,the fl uid e ff ects take their roles in the fuel target evolution.The anisotropic electron pressure will deteriorate the fuel core symmetry,reduce the density,and achieve a lower efficiency of nuclear fuel burning and a lower gain of nuclear reaction than expected.The e ff ects of the hot electron anisotropic pressure tensor in the fast ignition scheme should be considered as a factor in experiments where the nuclear reaction gain is measured to be much lower than the theoretical prediction.

    anisotropic electron pressure tensor,particle-in-cell simulation,fast ignition scheme

    10.7498/aps.66.115203

    ?國(guó)家自然科學(xué)基金(批準(zhǔn)號(hào):11674146,11274152)資助的課題.

    ?通信作者.E-mail:qldong@aphy.iphy.ac.cn

    ?2017中國(guó)物理學(xué)會(huì)Chinese Physical Society

    http://wulixb.iphy.ac.cn

    *Project supported by the National Natural Science Foundation of China(Grant Nos.11674146,11274152).

    ?Corresponding author.E-mail:qldong@aphy.iphy.ac.cn

    猜你喜歡
    高能張量動(dòng)量
    動(dòng)量守恒定律在三個(gè)物體系中的應(yīng)用
    前方高能!戰(zhàn)機(jī)懟瞼
    偶數(shù)階張量core逆的性質(zhì)和應(yīng)用
    四元數(shù)張量方程A*NX=B 的通解
    搞笑秀
    意林彩版(2020年3期)2020-11-22 20:54:12
    應(yīng)用動(dòng)量守恒定律解題之秘訣
    動(dòng)量相關(guān)知識(shí)的理解和應(yīng)用
    《高能少年團(tuán)》少年 未來(lái)可期
    擴(kuò)散張量成像MRI 在CO中毒后遲發(fā)腦病中的應(yīng)用
    Duang!6·18巾幗館前方高能
    海峽姐妹(2015年7期)2015-02-27 15:11:35
    最黄视频免费看| 国产又爽黄色视频| 成人亚洲精品一区在线观看| 中文字幕人妻熟女乱码| 日日爽夜夜爽网站| 欧美日韩黄片免| 一区二区日韩欧美中文字幕| 久久人妻熟女aⅴ| 黄色片一级片一级黄色片| 人人澡人人妻人| 亚洲欧美中文字幕日韩二区| 天天影视国产精品| 男女之事视频高清在线观看 | 久热爱精品视频在线9| av在线老鸭窝| 免费女性裸体啪啪无遮挡网站| 国产精品秋霞免费鲁丝片| av一本久久久久| 男人爽女人下面视频在线观看| 国产野战对白在线观看| 中文字幕另类日韩欧美亚洲嫩草| 国产一区亚洲一区在线观看| 中文乱码字字幕精品一区二区三区| 无遮挡黄片免费观看| 天天添夜夜摸| 各种免费的搞黄视频| av国产久精品久网站免费入址| 欧美日韩亚洲高清精品| 性高湖久久久久久久久免费观看| 在线观看免费午夜福利视频| 精品国产一区二区三区四区第35| 大陆偷拍与自拍| 别揉我奶头~嗯~啊~动态视频 | 韩国高清视频一区二区三区| 久久精品国产亚洲av涩爱| 丝袜人妻中文字幕| 久久久久精品国产欧美久久久 | 美女视频免费永久观看网站| 国产精品二区激情视频| 婷婷色综合大香蕉| 99国产精品99久久久久| 精品亚洲成a人片在线观看| 国产三级黄色录像| 老司机午夜十八禁免费视频| 91精品伊人久久大香线蕉| 中文字幕制服av| 国产一区二区三区综合在线观看| 午夜老司机福利片| 人妻人人澡人人爽人人| 免费黄频网站在线观看国产| 天天操日日干夜夜撸| av又黄又爽大尺度在线免费看| 亚洲精品国产av成人精品| 97在线人人人人妻| 日本a在线网址| 一本—道久久a久久精品蜜桃钙片| 丝袜美足系列| 欧美av亚洲av综合av国产av| 亚洲av欧美aⅴ国产| 亚洲精品国产一区二区精华液| 免费日韩欧美在线观看| 亚洲第一av免费看| 亚洲精品自拍成人| 我的亚洲天堂| 丝袜在线中文字幕| 午夜福利在线免费观看网站| 极品人妻少妇av视频| 亚洲中文日韩欧美视频| 国产精品久久久av美女十八| 日韩av不卡免费在线播放| 国产亚洲午夜精品一区二区久久| 丝袜脚勾引网站| 精品视频人人做人人爽| 一级片免费观看大全| 国产精品三级大全| 欧美日韩黄片免| 啦啦啦 在线观看视频| 久久精品国产a三级三级三级| 性高湖久久久久久久久免费观看| 精品欧美一区二区三区在线| 丝袜脚勾引网站| 另类亚洲欧美激情| 男女边摸边吃奶| 激情五月婷婷亚洲| 国产亚洲一区二区精品| 成人亚洲精品一区在线观看| 首页视频小说图片口味搜索 | 亚洲人成电影观看| 人人澡人人妻人| 欧美成人精品欧美一级黄| 亚洲av美国av| 国产精品av久久久久免费| 成人国产av品久久久| 欧美少妇被猛烈插入视频| 狂野欧美激情性bbbbbb| 亚洲成人国产一区在线观看 | 亚洲精品国产区一区二| 亚洲精品日韩在线中文字幕| 精品一区在线观看国产| 热re99久久精品国产66热6| 中文字幕人妻丝袜一区二区| 七月丁香在线播放| 精品少妇一区二区三区视频日本电影| 一二三四社区在线视频社区8| 久9热在线精品视频| 国产免费视频播放在线视频| 精品亚洲成a人片在线观看| 亚洲七黄色美女视频| 成人三级做爰电影| 9热在线视频观看99| 久久毛片免费看一区二区三区| 久久免费观看电影| av天堂久久9| 超碰97精品在线观看| 欧美日本中文国产一区发布| 首页视频小说图片口味搜索 | 久9热在线精品视频| 精品高清国产在线一区| 日韩 欧美 亚洲 中文字幕| 亚洲av欧美aⅴ国产| 欧美大码av| 亚洲精品国产区一区二| 国产伦人伦偷精品视频| 欧美大码av| 女人被躁到高潮嗷嗷叫费观| 亚洲av电影在线进入| 成人国语在线视频| 老司机在亚洲福利影院| 久久鲁丝午夜福利片| 日韩一卡2卡3卡4卡2021年| 2021少妇久久久久久久久久久| 在线观看一区二区三区激情| av网站免费在线观看视频| 国产视频首页在线观看| 亚洲av成人不卡在线观看播放网 | 亚洲精品乱久久久久久| 免费看av在线观看网站| 成在线人永久免费视频| 亚洲精品国产色婷婷电影| 老司机靠b影院| 人体艺术视频欧美日本| 欧美另类一区| 午夜精品国产一区二区电影| 国产福利在线免费观看视频| 久久国产精品人妻蜜桃| 无遮挡黄片免费观看| 国产精品 国内视频| 中国国产av一级| 黄色 视频免费看| 国产精品一区二区在线观看99| 亚洲av成人精品一二三区| 日韩免费高清中文字幕av| 日韩电影二区| 最黄视频免费看| 亚洲欧美清纯卡通| 亚洲第一青青草原| 在线av久久热| 丝袜在线中文字幕| 国产成人精品久久久久久| 久久午夜综合久久蜜桃| 在线亚洲精品国产二区图片欧美| av网站免费在线观看视频| 欧美日韩视频高清一区二区三区二| 久久久久久亚洲精品国产蜜桃av| 国产精品国产三级国产专区5o| 黄色视频不卡| 国产亚洲午夜精品一区二区久久| 午夜精品国产一区二区电影| 一级片免费观看大全| 免费在线观看日本一区| 精品人妻在线不人妻| 欧美精品一区二区免费开放| 老司机影院成人| 亚洲第一青青草原| 精品欧美一区二区三区在线| 啦啦啦啦在线视频资源| 大话2 男鬼变身卡| 亚洲中文日韩欧美视频| 熟女av电影| 精品国产一区二区久久| 五月天丁香电影| 精品少妇黑人巨大在线播放| 一级毛片 在线播放| 精品福利观看| 精品久久久久久久毛片微露脸 | 啦啦啦中文免费视频观看日本| 国产一区有黄有色的免费视频| 亚洲人成电影观看| 人妻人人澡人人爽人人| 一本一本久久a久久精品综合妖精| 久久这里只有精品19| 成年人午夜在线观看视频| 亚洲,欧美,日韩| 国产精品亚洲av一区麻豆| 亚洲三区欧美一区| 精品国产乱码久久久久久男人| av天堂久久9| 亚洲国产精品国产精品| 日韩大码丰满熟妇| 免费在线观看视频国产中文字幕亚洲 | 人妻人人澡人人爽人人| 男女床上黄色一级片免费看| 色综合欧美亚洲国产小说| 天堂俺去俺来也www色官网| 99精国产麻豆久久婷婷| 精品少妇久久久久久888优播| 无遮挡黄片免费观看| 久久99精品国语久久久| 91九色精品人成在线观看| 另类亚洲欧美激情| 自拍欧美九色日韩亚洲蝌蚪91| 蜜桃国产av成人99| av网站免费在线观看视频| 啦啦啦视频在线资源免费观看| 久久久国产欧美日韩av| 国产精品一国产av| 丝袜脚勾引网站| 看十八女毛片水多多多| 成人影院久久| 在线观看人妻少妇| 别揉我奶头~嗯~啊~动态视频 | 两个人免费观看高清视频| 亚洲欧洲国产日韩| 国产精品国产三级国产专区5o| 99精国产麻豆久久婷婷| 丁香六月欧美| 亚洲视频免费观看视频| 777米奇影视久久| 亚洲精品中文字幕在线视频| 欧美精品人与动牲交sv欧美| 又黄又粗又硬又大视频| 日韩,欧美,国产一区二区三区| 亚洲国产欧美在线一区| 嫩草影视91久久| 日韩av不卡免费在线播放| 精品卡一卡二卡四卡免费| 国产不卡av网站在线观看| av片东京热男人的天堂| 好男人电影高清在线观看| 天天操日日干夜夜撸| 热re99久久国产66热| 亚洲精品一二三| 日韩大片免费观看网站| 嫩草影视91久久| 日韩大码丰满熟妇| 曰老女人黄片| 一级毛片黄色毛片免费观看视频| 中文字幕av电影在线播放| av欧美777| 欧美日韩视频精品一区| 手机成人av网站| 性高湖久久久久久久久免费观看| 校园人妻丝袜中文字幕| 蜜桃在线观看..| 在线av久久热| 十八禁高潮呻吟视频| 国产精品偷伦视频观看了| 老司机亚洲免费影院| 国产欧美亚洲国产| 精品第一国产精品| 日韩欧美一区视频在线观看| 十分钟在线观看高清视频www| 激情视频va一区二区三区| 高清黄色对白视频在线免费看| 国产成人a∨麻豆精品| 国产又色又爽无遮挡免| 多毛熟女@视频| 人人妻,人人澡人人爽秒播 | 亚洲三区欧美一区| 男女边摸边吃奶| 亚洲自偷自拍图片 自拍| www.av在线官网国产| 午夜久久久在线观看| 涩涩av久久男人的天堂| 亚洲人成电影免费在线| 性色av一级| 大话2 男鬼变身卡| 我要看黄色一级片免费的| 亚洲精品美女久久av网站| 国产精品三级大全| 一区福利在线观看| 一边亲一边摸免费视频| 日韩视频在线欧美| 国产一区二区在线观看av| 妹子高潮喷水视频| 国产黄频视频在线观看| 久久国产精品影院| 精品福利永久在线观看| 亚洲一区中文字幕在线| 国产免费福利视频在线观看| av电影中文网址| 欧美日本中文国产一区发布| 制服人妻中文乱码| av视频免费观看在线观看| 日韩中文字幕视频在线看片| 极品少妇高潮喷水抽搐| 亚洲欧美成人综合另类久久久| 妹子高潮喷水视频| 国产精品成人在线| 水蜜桃什么品种好| 亚洲国产精品一区二区三区在线| 99国产精品一区二区三区| 午夜福利在线免费观看网站| 色播在线永久视频| 性高湖久久久久久久久免费观看| 日日摸夜夜添夜夜爱| 亚洲视频免费观看视频| 国产视频首页在线观看| 久久久欧美国产精品| xxx大片免费视频| 亚洲av成人不卡在线观看播放网 | 少妇猛男粗大的猛烈进出视频| 熟女av电影| 欧美另类一区| 另类亚洲欧美激情| 欧美精品av麻豆av| 午夜福利一区二区在线看| 国产男人的电影天堂91| 少妇精品久久久久久久| 可以免费在线观看a视频的电影网站| 亚洲专区中文字幕在线| av福利片在线| 91字幕亚洲| 美女主播在线视频| 18禁观看日本| 一区二区av电影网| 亚洲少妇的诱惑av| 老司机午夜十八禁免费视频| 十八禁高潮呻吟视频| 国产成人av教育| 美女视频免费永久观看网站| 国产精品二区激情视频| 亚洲情色 制服丝袜| 97在线人人人人妻| 欧美黄色淫秽网站| 在线观看一区二区三区激情| 777米奇影视久久| 波多野结衣av一区二区av| 日韩制服骚丝袜av| 亚洲 国产 在线| 99精品久久久久人妻精品| 亚洲精品中文字幕在线视频| 观看av在线不卡| 丰满饥渴人妻一区二区三| 视频在线观看一区二区三区| 久久国产精品影院| 一区二区三区精品91| 9热在线视频观看99| 99国产精品一区二区蜜桃av | 久久九九热精品免费| 97精品久久久久久久久久精品| 美女大奶头黄色视频| 曰老女人黄片| 中国美女看黄片| 亚洲五月色婷婷综合| 1024视频免费在线观看| 建设人人有责人人尽责人人享有的| av网站免费在线观看视频| a级毛片黄视频| 大话2 男鬼变身卡| 精品福利观看| 十八禁人妻一区二区| 日韩一区二区三区影片| 久久精品亚洲av国产电影网| 亚洲av综合色区一区| 在线观看免费日韩欧美大片| 老司机靠b影院| 久9热在线精品视频| 悠悠久久av| 99香蕉大伊视频| 午夜视频精品福利| 老司机午夜十八禁免费视频| 无遮挡黄片免费观看| 狂野欧美激情性xxxx| 午夜福利视频在线观看免费| xxx大片免费视频| 欧美97在线视频| 2018国产大陆天天弄谢| 亚洲人成网站在线观看播放| 久久影院123| 亚洲国产精品999| 在线观看国产h片| 国产欧美日韩一区二区三 | 久久 成人 亚洲| 19禁男女啪啪无遮挡网站| 久久 成人 亚洲| 国产成人精品在线电影| 我要看黄色一级片免费的| 51午夜福利影视在线观看| 国产精品久久久久成人av| 中文字幕另类日韩欧美亚洲嫩草| 一边摸一边抽搐一进一出视频| 亚洲av欧美aⅴ国产| 国产黄色视频一区二区在线观看| 婷婷色麻豆天堂久久| 秋霞在线观看毛片| 麻豆国产av国片精品| 国产高清视频在线播放一区 | 亚洲精品久久午夜乱码| 黄色一级大片看看| 国语对白做爰xxxⅹ性视频网站| 叶爱在线成人免费视频播放| 女性被躁到高潮视频| 久久久久久久久久久久大奶| 亚洲五月色婷婷综合| 国产精品偷伦视频观看了| 蜜桃国产av成人99| 国产成人一区二区三区免费视频网站 | 亚洲中文av在线| 好男人视频免费观看在线| 一级毛片黄色毛片免费观看视频| 下体分泌物呈黄色| 91国产中文字幕| 少妇粗大呻吟视频| 桃花免费在线播放| 国产淫语在线视频| 欧美久久黑人一区二区| 国产一区二区激情短视频 | 亚洲成人免费av在线播放| 成人国语在线视频| 考比视频在线观看| 国产99久久九九免费精品| 色94色欧美一区二区| netflix在线观看网站| 国产伦人伦偷精品视频| 侵犯人妻中文字幕一二三四区| 国产精品成人在线| 亚洲图色成人| 国产精品人妻久久久影院| 两个人免费观看高清视频| 欧美久久黑人一区二区| 日韩熟女老妇一区二区性免费视频| 91九色精品人成在线观看| 欧美日韩综合久久久久久| 国精品久久久久久国模美| 国产91精品成人一区二区三区 | 亚洲成av片中文字幕在线观看| 成人三级做爰电影| videos熟女内射| 久久这里只有精品19| 在线亚洲精品国产二区图片欧美| 18在线观看网站| 下体分泌物呈黄色| 亚洲国产精品一区二区三区在线| 精品久久久精品久久久| 夫妻性生交免费视频一级片| 久久 成人 亚洲| 欧美日本中文国产一区发布| 久久人人爽av亚洲精品天堂| 国产不卡av网站在线观看| 免费av中文字幕在线| 久热这里只有精品99| av网站在线播放免费| 日本欧美视频一区| 精品人妻在线不人妻| 18在线观看网站| av电影中文网址| 久久人妻福利社区极品人妻图片 | 午夜福利一区二区在线看| 午夜久久久在线观看| h视频一区二区三区| 国产精品欧美亚洲77777| 国产真人三级小视频在线观看| 一本—道久久a久久精品蜜桃钙片| 亚洲av成人不卡在线观看播放网 | 少妇人妻 视频| 欧美 日韩 精品 国产| 日本av免费视频播放| 精品少妇一区二区三区视频日本电影| 久久女婷五月综合色啪小说| 国产精品免费大片| 日韩精品免费视频一区二区三区| 精品一区二区三区四区五区乱码 | 日韩精品免费视频一区二区三区| 国产精品三级大全| 考比视频在线观看| 国产视频首页在线观看| 国产91精品成人一区二区三区 | 国产精品偷伦视频观看了| 亚洲精品成人av观看孕妇| 欧美亚洲 丝袜 人妻 在线| 一区二区三区激情视频| 青草久久国产| 亚洲欧美日韩另类电影网站| 婷婷色麻豆天堂久久| 777久久人妻少妇嫩草av网站| 精品国产乱码久久久久久男人| 色播在线永久视频| 国产三级黄色录像| 日日夜夜操网爽| 在线精品无人区一区二区三| 国产一区二区 视频在线| 亚洲精品一二三| 制服人妻中文乱码| 亚洲第一青青草原| 久久国产精品人妻蜜桃| 久久亚洲国产成人精品v| 精品一区在线观看国产| 国产精品国产av在线观看| 国产成人啪精品午夜网站| 亚洲国产看品久久| 欧美 亚洲 国产 日韩一| 免费女性裸体啪啪无遮挡网站| 一本大道久久a久久精品| 中文字幕另类日韩欧美亚洲嫩草| 人妻人人澡人人爽人人| 中文字幕另类日韩欧美亚洲嫩草| 中文字幕精品免费在线观看视频| 亚洲九九香蕉| 欧美大码av| 9色porny在线观看| 蜜桃国产av成人99| 亚洲国产欧美一区二区综合| 水蜜桃什么品种好| 中文字幕人妻丝袜制服| 久久 成人 亚洲| 五月开心婷婷网| 亚洲国产欧美网| 性色av乱码一区二区三区2| 欧美精品亚洲一区二区| 国产女主播在线喷水免费视频网站| 婷婷色综合www| 国产福利在线免费观看视频| 99热全是精品| a级毛片在线看网站| 国产日韩欧美视频二区| 欧美日韩视频精品一区| 亚洲av在线观看美女高潮| 777米奇影视久久| 人妻 亚洲 视频| 午夜激情久久久久久久| h视频一区二区三区| 国产激情久久老熟女| 国产淫语在线视频| 国产成人啪精品午夜网站| 午夜福利乱码中文字幕| 国产精品 欧美亚洲| 少妇被粗大的猛进出69影院| 男女国产视频网站| 国产av一区二区精品久久| 伦理电影免费视频| 亚洲国产精品999| 脱女人内裤的视频| av有码第一页| 午夜福利视频精品| 国产av精品麻豆| av电影中文网址| 久久精品国产亚洲av高清一级| 手机成人av网站| 十八禁人妻一区二区| 婷婷成人精品国产| 美女主播在线视频| 老熟女久久久| 狂野欧美激情性xxxx| 久久亚洲精品不卡| 久久影院123| 亚洲国产日韩一区二区| 精品卡一卡二卡四卡免费| 狠狠婷婷综合久久久久久88av| 热99国产精品久久久久久7| cao死你这个sao货| 亚洲三区欧美一区| av福利片在线| 91成人精品电影| 91老司机精品| 亚洲欧美一区二区三区黑人| 三上悠亚av全集在线观看| 免费久久久久久久精品成人欧美视频| 视频区欧美日本亚洲| 国产视频一区二区在线看| 国产97色在线日韩免费| 亚洲精品美女久久久久99蜜臀 | 精品国产乱码久久久久久男人| 亚洲av综合色区一区| 一边亲一边摸免费视频| 久久女婷五月综合色啪小说| 日韩一区二区三区影片| 可以免费在线观看a视频的电影网站| 欧美日本中文国产一区发布| 九色亚洲精品在线播放| 精品一区二区三区四区五区乱码 | 日本猛色少妇xxxxx猛交久久| 90打野战视频偷拍视频| 久久人人爽人人片av| 日韩大片免费观看网站| 国产欧美亚洲国产| 亚洲五月色婷婷综合| 爱豆传媒免费全集在线观看| 国产在线免费精品| 99热全是精品| 91麻豆av在线| 亚洲成色77777| 热re99久久精品国产66热6| 成人三级做爰电影| 免费高清在线观看日韩| 欧美日韩黄片免| www.精华液| 99精国产麻豆久久婷婷| 美女视频免费永久观看网站| 免费av中文字幕在线| 少妇裸体淫交视频免费看高清 | 国产av一区二区精品久久| 一本综合久久免费| 男人操女人黄网站| 久久亚洲精品不卡| 国语对白做爰xxxⅹ性视频网站| 麻豆国产av国片精品| 国产精品久久久人人做人人爽| 欧美人与性动交α欧美软件| 女警被强在线播放| 1024香蕉在线观看| 亚洲一区中文字幕在线| e午夜精品久久久久久久| 国产欧美亚洲国产| 亚洲三区欧美一区|