• <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
    欧美绝顶高潮抽搐喷水| 国产精品 欧美亚洲| 久久精品夜夜夜夜夜久久蜜豆| 九色国产91popny在线| 国产av不卡久久| 好看av亚洲va欧美ⅴa在| 最近最新免费中文字幕在线| 黄色日韩在线| 制服丝袜大香蕉在线| 日韩人妻高清精品专区| 一进一出抽搐gif免费好疼| 国产av在哪里看| 黄片小视频在线播放| 久久香蕉国产精品| 99国产精品一区二区蜜桃av| 十八禁人妻一区二区| 成人av在线播放网站| 国产精品亚洲av一区麻豆| 国产在视频线在精品| 国产精品电影一区二区三区| 色老头精品视频在线观看| 91久久精品电影网| 国产av在哪里看| 欧美日本亚洲视频在线播放| 国产精品久久久久久精品电影| 他把我摸到了高潮在线观看| 91av网一区二区| 人妻夜夜爽99麻豆av| 女同久久另类99精品国产91| 国产精品99久久99久久久不卡| 网址你懂的国产日韩在线| 1000部很黄的大片| 在线观看免费视频日本深夜| 在线观看av片永久免费下载| 日韩欧美 国产精品| 亚洲成人久久性| 国产精品电影一区二区三区| a在线观看视频网站| 色在线成人网| 欧美乱色亚洲激情| 久99久视频精品免费| 久久香蕉国产精品| 三级国产精品欧美在线观看| 成熟少妇高潮喷水视频| 亚洲内射少妇av| 久久精品91蜜桃| bbb黄色大片| 亚洲真实伦在线观看| 久久草成人影院| 狂野欧美白嫩少妇大欣赏| 欧美中文日本在线观看视频| 久久伊人香网站| 亚洲在线观看片| 国产综合懂色| 久久人妻av系列| 三级毛片av免费| 欧美成人性av电影在线观看| 男女之事视频高清在线观看| 欧美色视频一区免费| 精品午夜福利视频在线观看一区| 午夜精品一区二区三区免费看| 一个人看视频在线观看www免费 | 最新中文字幕久久久久| 色吧在线观看| 哪里可以看免费的av片| 精品久久久久久久毛片微露脸| 午夜福利免费观看在线| 丰满乱子伦码专区| 久久久久久国产a免费观看| 乱人视频在线观看| 一个人免费在线观看电影| 国产精品三级大全| 超碰av人人做人人爽久久 | 少妇裸体淫交视频免费看高清| 欧美一区二区亚洲| 最近视频中文字幕2019在线8| 少妇裸体淫交视频免费看高清| 丝袜美腿在线中文| 欧美在线一区亚洲| 日韩高清综合在线| 无人区码免费观看不卡| 中出人妻视频一区二区| 中出人妻视频一区二区| 日本在线视频免费播放| 日日夜夜操网爽| 又紧又爽又黄一区二区| 最好的美女福利视频网| 88av欧美| 欧美成人性av电影在线观看| 国产成人av教育| 国产真实乱freesex| 亚洲精品456在线播放app | 黄片大片在线免费观看| 精品久久久久久久末码| 午夜精品在线福利| 欧美中文日本在线观看视频| 亚洲av成人精品一区久久| 神马国产精品三级电影在线观看| 国产av一区在线观看免费| 亚洲精品影视一区二区三区av| 久久九九热精品免费| 久久中文看片网| 在线十欧美十亚洲十日本专区| 国产精品乱码一区二三区的特点| 91麻豆av在线| 母亲3免费完整高清在线观看| 欧美成人性av电影在线观看| 久久人妻av系列| 亚洲七黄色美女视频| 国产精品电影一区二区三区| 69av精品久久久久久| 一进一出抽搐gif免费好疼| 亚洲成av人片在线播放无| 99久久久亚洲精品蜜臀av| 成人三级黄色视频| 色老头精品视频在线观看| 两人在一起打扑克的视频| 亚洲欧美激情综合另类| www.www免费av| 欧美性猛交黑人性爽| 欧美又色又爽又黄视频| 一进一出抽搐gif免费好疼| 九九在线视频观看精品| 欧美不卡视频在线免费观看| 在线观看舔阴道视频| 久久精品亚洲精品国产色婷小说| 真人做人爱边吃奶动态| 国产亚洲精品综合一区在线观看| 麻豆一二三区av精品| 亚洲午夜理论影院| 欧美+日韩+精品| 亚洲最大成人中文| 老熟妇乱子伦视频在线观看| 欧美不卡视频在线免费观看| 亚洲色图av天堂| 久久久久国产精品人妻aⅴ院| 国产成+人综合+亚洲专区| 少妇高潮的动态图| 中亚洲国语对白在线视频| 老鸭窝网址在线观看| 国产精品久久视频播放| 嫁个100分男人电影在线观看| 亚洲精品久久国产高清桃花| 国产精品久久视频播放| 毛片女人毛片| 禁无遮挡网站| 午夜精品在线福利| 美女cb高潮喷水在线观看| 18禁裸乳无遮挡免费网站照片| 亚洲精品久久国产高清桃花| 俺也久久电影网| 欧美成人一区二区免费高清观看| 毛片女人毛片| 三级毛片av免费| 午夜激情福利司机影院| 成年免费大片在线观看| 老司机福利观看| 99久久成人亚洲精品观看| 中文字幕人妻熟人妻熟丝袜美 | 久久国产乱子伦精品免费另类| 亚洲成人精品中文字幕电影| 午夜福利欧美成人| 日韩欧美一区二区三区在线观看| 长腿黑丝高跟| 免费大片18禁| 母亲3免费完整高清在线观看| 男人的好看免费观看在线视频| 久久久国产成人免费| 久久久国产成人精品二区| 高清在线国产一区| 99热精品在线国产| 韩国av一区二区三区四区| 成人国产综合亚洲| 成人av一区二区三区在线看| 成人国产一区最新在线观看| 午夜福利18| 国产免费av片在线观看野外av| АⅤ资源中文在线天堂| 在线观看美女被高潮喷水网站 | 首页视频小说图片口味搜索| 精品久久久久久成人av| 怎么达到女性高潮| 国产熟女xx| 久久久久国产精品人妻aⅴ院| 成年女人看的毛片在线观看| 久久人人精品亚洲av| 久久亚洲精品不卡| 免费电影在线观看免费观看| 色综合站精品国产| 一个人观看的视频www高清免费观看| 高清毛片免费观看视频网站| 国产麻豆成人av免费视频| 色综合欧美亚洲国产小说| 国产aⅴ精品一区二区三区波| 国产精品电影一区二区三区| 香蕉久久夜色| 国产精品亚洲一级av第二区| 国产爱豆传媒在线观看| 18禁黄网站禁片免费观看直播| 色吧在线观看| 欧美性感艳星| 亚洲五月天丁香| 午夜精品久久久久久毛片777| 97超级碰碰碰精品色视频在线观看| 亚洲国产中文字幕在线视频| 欧美最黄视频在线播放免费| 在线十欧美十亚洲十日本专区| 在线观看免费视频日本深夜| 亚洲真实伦在线观看| 国产精品久久久人人做人人爽| 亚洲精品色激情综合| 亚洲,欧美精品.| 国产精品,欧美在线| 日本 av在线| 搡女人真爽免费视频火全软件 | 色av中文字幕| 亚洲avbb在线观看| 禁无遮挡网站| 欧美日韩一级在线毛片| x7x7x7水蜜桃| 99久国产av精品| 亚洲国产精品久久男人天堂| 国产在视频线在精品| 啦啦啦观看免费观看视频高清| 美女被艹到高潮喷水动态| 国内精品久久久久精免费| 日韩大尺度精品在线看网址| 美女被艹到高潮喷水动态| 免费av观看视频| 国产精品亚洲美女久久久| 无限看片的www在线观看| 欧美日本视频| 内射极品少妇av片p| 国产精品98久久久久久宅男小说| 18+在线观看网站| 免费一级毛片在线播放高清视频| 久久久久久久久中文| 人人妻人人看人人澡| 熟妇人妻久久中文字幕3abv| 老熟妇乱子伦视频在线观看| 小说图片视频综合网站| 国产精品久久久久久久电影 | 69人妻影院| 亚洲欧美日韩卡通动漫| 国产精华一区二区三区| 久9热在线精品视频| 狂野欧美白嫩少妇大欣赏| 亚洲成人久久爱视频| 成年免费大片在线观看| 欧美成人免费av一区二区三区| 精品99又大又爽又粗少妇毛片 | 成人特级黄色片久久久久久久| 中文亚洲av片在线观看爽| 欧美一区二区精品小视频在线| 女人高潮潮喷娇喘18禁视频| 91麻豆精品激情在线观看国产| 免费高清视频大片| 九九热线精品视视频播放| 19禁男女啪啪无遮挡网站| 在线观看美女被高潮喷水网站 | 啦啦啦免费观看视频1| 亚洲专区中文字幕在线| 国产高潮美女av| 国产免费一级a男人的天堂| 久久国产精品影院| 欧美精品啪啪一区二区三区| 国产成人a区在线观看| 制服人妻中文乱码| 国产伦精品一区二区三区四那| netflix在线观看网站| netflix在线观看网站| 国产乱人伦免费视频| 精品一区二区三区视频在线观看免费| 每晚都被弄得嗷嗷叫到高潮| 午夜福利高清视频| 神马国产精品三级电影在线观看| 国产一区二区在线观看日韩 | а√天堂www在线а√下载| 欧美黑人巨大hd| 成人无遮挡网站| av片东京热男人的天堂| 丁香六月欧美| 久久香蕉精品热| 亚洲精品日韩av片在线观看 | 99在线视频只有这里精品首页| 国产精品一区二区免费欧美| 少妇高潮的动态图| 91av网一区二区| 欧美性感艳星| 亚洲人成伊人成综合网2020| 成人欧美大片| 日本免费一区二区三区高清不卡| 国产av不卡久久| 久99久视频精品免费| 国产 一区 欧美 日韩| 伊人久久精品亚洲午夜| 午夜福利在线在线| 欧美性猛交黑人性爽| 免费在线观看成人毛片| 高潮久久久久久久久久久不卡| 两性午夜刺激爽爽歪歪视频在线观看| 人人妻,人人澡人人爽秒播| 成年女人永久免费观看视频| 1024手机看黄色片| 99热精品在线国产| 国产精品综合久久久久久久免费| 久久精品国产清高在天天线| 午夜老司机福利剧场| 国产 一区 欧美 日韩| 中文亚洲av片在线观看爽| 高清日韩中文字幕在线| 免费看十八禁软件| 国产av一区在线观看免费| 欧美不卡视频在线免费观看| av国产免费在线观看| 久久6这里有精品| 精品久久久久久成人av| 在线观看舔阴道视频| 欧美日本亚洲视频在线播放| 在线观看美女被高潮喷水网站 | 精品久久久久久久毛片微露脸| 麻豆成人午夜福利视频| 国产精品一区二区免费欧美| 欧美日韩黄片免| 波多野结衣高清作品| 精品欧美国产一区二区三| 亚洲av成人精品一区久久| 亚洲中文字幕一区二区三区有码在线看| 午夜免费观看网址| 男插女下体视频免费在线播放| 黄色丝袜av网址大全| 国产乱人伦免费视频| 日本 av在线| 久久久久国内视频| 亚洲天堂国产精品一区在线| 在线a可以看的网站| 国产私拍福利视频在线观看| 亚洲精品影视一区二区三区av| 日韩欧美在线二视频| 三级男女做爰猛烈吃奶摸视频| 91麻豆精品激情在线观看国产| 性色avwww在线观看| 国产伦一二天堂av在线观看| 婷婷精品国产亚洲av在线| 一个人免费在线观看的高清视频| 两人在一起打扑克的视频| 无限看片的www在线观看| 五月伊人婷婷丁香| 亚洲aⅴ乱码一区二区在线播放| 蜜桃亚洲精品一区二区三区| 免费av观看视频| 国产真实伦视频高清在线观看 | 免费av观看视频| 2021天堂中文幕一二区在线观| 少妇裸体淫交视频免费看高清| 欧美日韩亚洲国产一区二区在线观看| h日本视频在线播放| 88av欧美| 欧美日韩福利视频一区二区| 国产av麻豆久久久久久久| 久久人人精品亚洲av| 69av精品久久久久久| 久久伊人香网站| 在线天堂最新版资源| 亚洲18禁久久av| 亚洲人成电影免费在线| 亚洲午夜理论影院| 亚洲精品日韩av片在线观看 | 亚洲成a人片在线一区二区| 少妇熟女aⅴ在线视频| 久久精品夜夜夜夜夜久久蜜豆| 性色avwww在线观看| 99热精品在线国产| 人人妻,人人澡人人爽秒播| 18禁美女被吸乳视频| 蜜桃久久精品国产亚洲av| 在线天堂最新版资源| АⅤ资源中文在线天堂| 国内精品久久久久久久电影| 无限看片的www在线观看| 国内揄拍国产精品人妻在线| 欧美成人一区二区免费高清观看| 黑人欧美特级aaaaaa片| ponron亚洲| 欧美一区二区亚洲| 精品午夜福利视频在线观看一区| 国产精品99久久99久久久不卡| 亚洲国产欧洲综合997久久,| 久9热在线精品视频| 中文字幕av在线有码专区| 国产三级在线视频| 两个人视频免费观看高清| 欧美国产日韩亚洲一区| 不卡一级毛片| 成人特级av手机在线观看| 女人被狂操c到高潮| 中文字幕高清在线视频| 亚洲自拍偷在线| 成人无遮挡网站| 国产欧美日韩精品一区二区| 一级a爱片免费观看的视频| 久久久久九九精品影院| 免费在线观看日本一区| 日本与韩国留学比较| 中国美女看黄片| 可以在线观看毛片的网站| 国产男靠女视频免费网站| а√天堂www在线а√下载| 中文在线观看免费www的网站| 老汉色∧v一级毛片| 国内毛片毛片毛片毛片毛片| 国产精品 欧美亚洲| 999久久久精品免费观看国产| 高清日韩中文字幕在线| 久久久国产成人免费| 欧美乱码精品一区二区三区| 国产蜜桃级精品一区二区三区| 国产精品嫩草影院av在线观看 | 午夜精品久久久久久毛片777| 一个人看视频在线观看www免费 | 狂野欧美激情性xxxx| 午夜两性在线视频| 90打野战视频偷拍视频| 成人亚洲精品av一区二区| 不卡一级毛片| 午夜亚洲福利在线播放| 真实男女啪啪啪动态图| 女人高潮潮喷娇喘18禁视频| 久久精品91无色码中文字幕| 免费一级毛片在线播放高清视频| 亚洲成人中文字幕在线播放| 中文字幕人成人乱码亚洲影| 亚洲精品影视一区二区三区av| 可以在线观看毛片的网站| 国产真实伦视频高清在线观看 | 欧美不卡视频在线免费观看| 波多野结衣高清作品| 久久久国产成人精品二区| 变态另类成人亚洲欧美熟女| a级毛片a级免费在线| 亚洲欧美激情综合另类| 成人精品一区二区免费| 国产伦精品一区二区三区视频9 | 高潮久久久久久久久久久不卡| 丁香欧美五月| 好男人在线观看高清免费视频| 久久久久国产精品人妻aⅴ院| 国产一区二区在线av高清观看| 欧美另类亚洲清纯唯美| 国产极品精品免费视频能看的| 一个人免费在线观看电影| 国产高清三级在线| 国产成人系列免费观看| 久久久色成人| 国产精品一区二区三区四区免费观看 | 亚洲av不卡在线观看| 91字幕亚洲| 最后的刺客免费高清国语| 中文字幕av在线有码专区| 国产免费av片在线观看野外av| 日本免费一区二区三区高清不卡| 亚洲欧美日韩卡通动漫| 日日夜夜操网爽| 一进一出好大好爽视频| 首页视频小说图片口味搜索| 我的老师免费观看完整版| 日本a在线网址| 亚洲av中文字字幕乱码综合| 国产免费av片在线观看野外av| 桃色一区二区三区在线观看| 欧美区成人在线视频| 亚洲av电影在线进入| 男女那种视频在线观看| 亚洲国产精品999在线| 午夜精品在线福利| 亚洲色图av天堂| 中文字幕熟女人妻在线| 午夜久久久久精精品| 久久国产乱子伦精品免费另类| 久久久色成人| 亚洲色图av天堂| 乱人视频在线观看| 观看免费一级毛片| 一区福利在线观看| 18禁黄网站禁片午夜丰满| 国产黄片美女视频| 国产精品自产拍在线观看55亚洲| 看免费av毛片| 国产免费av片在线观看野外av| 老司机午夜福利在线观看视频| 在线观看av片永久免费下载| 久久久久久久精品吃奶| 小蜜桃在线观看免费完整版高清| 国产私拍福利视频在线观看| 欧美成人免费av一区二区三区| 伊人久久大香线蕉亚洲五| 桃红色精品国产亚洲av| 岛国在线观看网站| 色综合站精品国产| 真实男女啪啪啪动态图| 少妇丰满av| 欧美中文日本在线观看视频| av在线蜜桃| 国产精品 国内视频| 搡老熟女国产l中国老女人| 国产一区二区亚洲精品在线观看| 中文字幕人妻熟人妻熟丝袜美 | 中文字幕精品亚洲无线码一区| 麻豆成人av在线观看| 亚洲七黄色美女视频| 欧美高清成人免费视频www| 一二三四社区在线视频社区8| 国产久久久一区二区三区| 国产真实伦视频高清在线观看 | 99热这里只有是精品50| 亚洲成人中文字幕在线播放| 国产欧美日韩一区二区精品| 看黄色毛片网站| 3wmmmm亚洲av在线观看| 一进一出抽搐动态| 香蕉久久夜色| 国产私拍福利视频在线观看| 日韩欧美在线二视频| 久久国产乱子伦精品免费另类| 操出白浆在线播放| 国产一区二区三区视频了| 十八禁人妻一区二区| 成人国产综合亚洲| 99国产综合亚洲精品| 亚洲中文日韩欧美视频| 亚洲精品国产精品久久久不卡| 国产色爽女视频免费观看| 男人舔女人下体高潮全视频| 欧美中文日本在线观看视频| 国产精品久久久久久亚洲av鲁大| 1024手机看黄色片| 精品一区二区三区视频在线观看免费| 亚洲黑人精品在线| 两性午夜刺激爽爽歪歪视频在线观看| 午夜福利免费观看在线| 美女 人体艺术 gogo| 久久精品91蜜桃| www.www免费av| 久久香蕉精品热| 一区二区三区免费毛片| 手机成人av网站| 日本免费a在线| 久久久久免费精品人妻一区二区| 一个人免费在线观看的高清视频| 脱女人内裤的视频| 亚洲18禁久久av| 欧美大码av| 亚洲美女视频黄频| 亚洲久久久久久中文字幕| 国产真实伦视频高清在线观看 | 女人高潮潮喷娇喘18禁视频| 国产中年淑女户外野战色| 中文字幕av成人在线电影| 欧美色视频一区免费| 夜夜看夜夜爽夜夜摸| 嫁个100分男人电影在线观看| 亚洲av二区三区四区| 欧美成人a在线观看| 亚洲色图av天堂| 欧美成人性av电影在线观看| 日韩国内少妇激情av| 18禁美女被吸乳视频| 一个人看视频在线观看www免费 | 精品久久久久久久人妻蜜臀av| 国产高清视频在线观看网站| 久久国产精品人妻蜜桃| 99久久精品热视频| 亚洲第一电影网av| 亚洲激情在线av| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 一区二区三区国产精品乱码| 麻豆国产97在线/欧美| 高潮久久久久久久久久久不卡| 波多野结衣高清作品| 看免费av毛片| 国产一区二区在线av高清观看| 欧美最黄视频在线播放免费| 成年女人永久免费观看视频| 国产欧美日韩精品一区二区| 亚洲精品国产精品久久久不卡| 亚洲欧美日韩高清在线视频| 长腿黑丝高跟| 日日夜夜操网爽| 久久精品综合一区二区三区| 中文字幕av在线有码专区| 亚洲狠狠婷婷综合久久图片| 脱女人内裤的视频| 国产一区二区激情短视频| 午夜免费男女啪啪视频观看 | 亚洲内射少妇av| 国产伦精品一区二区三区四那| 久久久久久久午夜电影| 亚洲国产欧美网| 少妇裸体淫交视频免费看高清| 日本熟妇午夜| 夜夜爽天天搞| 99视频精品全部免费 在线| 禁无遮挡网站| 免费在线观看影片大全网站| 日韩亚洲欧美综合| 亚洲精品亚洲一区二区| e午夜精品久久久久久久| 亚洲成人中文字幕在线播放| 国产精品久久电影中文字幕| 精品国产美女av久久久久小说| 国产高清三级在线| 99在线人妻在线中文字幕| 波多野结衣巨乳人妻| 亚洲性夜色夜夜综合| 欧美日韩瑟瑟在线播放|