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

    吸積盤驅(qū)動噴流模型研究進展

    2020-01-18 01:21:10王薦釗錢宇梁仲佳勇
    天文學進展 2019年4期
    關(guān)鍵詞:磁力線噴流天體

    王薦釗,錢宇梁,仲佳勇

    (北京師范大學天文系,北京100875)

    1 引 言

    高速準直的噴流是非常普遍和重要的天文觀測現(xiàn)象之一,它們是來自活動星系核(active galactic nucleus, AGN)、X 射線雙星、原初恒星體(young stellar object, YSO)等天體中心的高強度等離子體流。這些高速旋轉(zhuǎn)并吸收周圍環(huán)境中物質(zhì)的致密中心天體包括原初恒星、中子星、恒星質(zhì)量黑洞和星系中央的超大質(zhì)量黑洞等。自M87 星系中心光學噴流被發(fā)現(xiàn)以來,噴流的產(chǎn)生、加速和準直機制一直是天體物理研究的前沿熱點之一[1]。人們普遍認為,這種大尺度噴流與致密天體周圍的吸積盤有關(guān)。吸積盤是存在于中心天體周圍含有吸積流的轉(zhuǎn)動氣體盤,它們會導致或激發(fā)各種天文現(xiàn)象,如恒星形成、高能輻射和天體噴流等[2]。半個世紀以來,人們建立了四種重要的吸積盤模型。1973 年Shakura 和Sunyaev[3]建立了幾何薄、光學厚的標準薄盤模型,即SSD (Shakura Sunyaev disk)模型。然而通過這個模型計算出的盤溫度偏低,無法解釋X 射線波段到極紫外波段輻射的觀測事實。為解決SSD 模型的不穩(wěn)定性問題以及高能輻射問題,1976 年Shapiro等人[4]提出了幾何薄、光學薄的SLE (Shapiro-Lightman-Eardley)盤模型。然而這種盤因為熱不穩(wěn)定性而不能真實存在。1988 年Abramowicz 等人[5]建立了幾何厚、光學厚的輻射壓主導的細盤(slim)模型。為了解釋致密天體系統(tǒng)的硬X 射線來源,Narayan 和Abramowicz 等人[6]又建立了徑移主導的吸積流(advection-dominated accretion flow, ADAF)模型。然而對于真實的復雜天體環(huán)境,人們通常會利用雙模式解來描述黑洞X 射線雙星和AGN 等天體系統(tǒng)的吸積盤。了解中心天體及其周圍的吸積盤,有助于人們理解天體噴流和出流的形成機制及其傳輸演化過程。

    雖然在過去幾十年里科學家們提出了很多模型,但從致密天體系統(tǒng)產(chǎn)生強相對論性噴流的具體機制目前還不清楚[7]。最近的研究表明,高速噴流與連接吸積盤的扭曲磁場結(jié)構(gòu)有關(guān),主要有兩種機制:磁流體機制和坡印亭通量機制(Poynting flux regime)。在磁流體機制中,能量和角動量被電磁場和物質(zhì)動能通量所攜帶,也就是說,這種機制中物質(zhì)起主導作用,電磁場起約束作用。這種機制主要與原初恒星體產(chǎn)生噴流過程有關(guān)。坡印亭通量機制中,能量和角動量主要被電磁場所攜帶,坡印亭能流起主導作用,物質(zhì)成分很少。這種機制主要與微類星體噴流和γ射線暴過程有關(guān)[8]。這兩種機制分別對應BP (Blandford-Payne)和BZ (Blandford-Znajek)過程[9]。雖然在噴流形成中處于主導地位的機制有待進一步研究,但連接黑洞和吸積盤系統(tǒng)的大尺度磁場驅(qū)動了噴流這一觀點已被大多數(shù)學者所認同。這樣的噴流速度可以在每秒幾十千米到接近光速之間變化。噴流的總能量很大部分是在宿主形成過程中釋放的引力能,且該引力能大到足以使噴流脫離宿主致密天體。本文列舉了一些典型的吸積供能天體,簡述了當前主流的幾種吸積盤噴流模型,并介紹了近年來計算機模擬吸積盤驅(qū)動噴流以及實驗室相關(guān)模擬噴流的最新研究進展。

    2 吸積盤噴流模型

    2.1 BZ 過程與BP 過程

    1977 年Blandford 和Znajek[10]首先提出,活動星系核噴流受帶有大尺度磁場的旋轉(zhuǎn)黑洞驅(qū)動(BZ 過程),吸積系統(tǒng)通過大尺度磁場提取黑洞轉(zhuǎn)動能。BZ 過程中,黑洞視界面附近的大尺度磁場被中心克爾黑洞拖曳而隨黑洞一起繞轉(zhuǎn)。大尺度磁場的另一端與遙遠的天體物理負載相連,天體物理負載可以是外部載入的(被磁力線束縛的)沿磁力線螺旋運動的等離子體,也可以是與磁力線另一端“凍結(jié)”在一起運動的等離子體物質(zhì)。當外部載入物質(zhì)與磁力線存在角速度差時,磁力線與載入物質(zhì)間會產(chǎn)生力矩作用,進而傳遞角動量和能量。這里的磁力線實質(zhì)上是旋轉(zhuǎn)黑洞與載入物質(zhì)間傳遞角動量和能量的媒介。一般認為,磁力線角速度必須保持恒定,在離黑洞旋轉(zhuǎn)軸足夠遠的地方,磁力線上某一點的轉(zhuǎn)速將超過光速,而沿磁力線運動的等離子體運動速度不能超過光速,因而等離子體不可能沿磁力線一直運動,當?shù)入x子體的速度被加速到接近光速時,它就會脫離磁力線的束縛。實際上,當?shù)入x子體的動能密度大于當?shù)卮蟪叨却艌龅拇拍苊芏葧r,磁場就不能再束縛住等離子體。等離子體的這一加速過程,就像串在繞軸旋轉(zhuǎn)的木棍上的珠子被甩出去一樣[11]。最初人們用BZ 過程成功地解釋了類星體和活動星系核中磁場提取黑洞旋轉(zhuǎn)能的過程,它們最多可提取31%的旋轉(zhuǎn)能,這可謂是一種有效的提能機制。同時,用BZ 過程可以很好地解釋觀測上的相對論性噴流、磁場對噴流的準直效應以及γ射線暴的快速光變等問題[12]。

    除BZ 過程中的磁場位形外,還存在著連接吸積盤和外部天體負載的大尺度有序磁場。1982 年Blandford 和Payne[13]假設(shè)大尺度磁場“凍結(jié)”在吸積盤上。他們提出,隨著吸積盤以開普勒角速度轉(zhuǎn)動,物質(zhì)的外流會受磁場離心作用驅(qū)動,吸積系統(tǒng)通過大尺度磁場提取吸積物質(zhì)的引力能。若盤平面法向與極向磁場的夾角大于臨界值,就會有物質(zhì)受離心作用而流出(BP 過程)。類似于BZ 過程,吸積盤作為等效電路的“電源”,向外部負載輸出能量。同時,吸積盤也是噴流物質(zhì)載入的重要來源,吸積盤上的物質(zhì)載入到大尺度磁場并通過磁場加速,因而大尺度磁場除了向外輸出坡印亭能流外,還向從吸積盤載入并沿磁力線加速的等離子體轉(zhuǎn)移能量和角動量。不同于BZ 過程,BP 過程產(chǎn)生的噴流攜帶有大量重子物質(zhì),這一現(xiàn)象可以較好地解釋低速外流的產(chǎn)生機制[14]。同時,BP 過程帶來的物質(zhì)外流會導致盤光度明顯下降,此外,噴流出現(xiàn)時盤的X 射線輻射通量也會隨之減小[15]。

    2.2 BZ-MC-BP 過程中的磁場結(jié)構(gòu)

    除了BZ 過程和BP 過程,旋轉(zhuǎn)黑洞與它周圍吸積盤的耦合過程,即磁耦合(magnetic coupling, MC)過程也引起了天文學家的注意。通過MC 過程,能量和角動量可以從旋轉(zhuǎn)黑洞被傳送到周圍的盤中,并在盤吸積過程中起重要作用。這里將BZ 過程、MC 過程和BP過程結(jié)合起來,以解釋活動星系核和黑洞雙星中噴流的產(chǎn)生機制(以下稱作BZ-MC-BP 模型,見圖1)。該模型與觀測到的恒星級黑洞和超大質(zhì)量黑洞噴流的結(jié)構(gòu)相一致。

    圖1 BZ-MC-BP 模型磁場結(jié)構(gòu)示意圖[9]

    當磁場既包含極向成分,又包含環(huán)向成分,并且環(huán)向分量很強,以至于磁力線環(huán)繞著其自身時,磁場結(jié)構(gòu)會產(chǎn)生螺旋不穩(wěn)定性[16]。在某種條件下,BZ 過程的螺旋不穩(wěn)定性(screwinstability in the BZ process, SIBZ)可以與MC 過程的螺旋不穩(wěn)定性(screw instability in the MC process, SIMC)共存,因此借助磁場螺旋不穩(wěn)定性把三種過程考慮到一個模型,即BZ-MC-BP 模型中,其磁場結(jié)構(gòu)如圖1 所示。圖1 中,θs表示黑洞視界上開放磁力線與閉合磁力線的臨界角度;θL是閉合磁力線的最小邊界角度,這里取θL= 0.45 π;rSIMC是通過SIMC 過程限制的磁場臨界半徑;rout表示盤的外邊界半徑;rms是吸積盤內(nèi)邊界到黑洞中心的距離;rSIBZ表示由SIBZ 過程限制的臨界半徑;αFL表示磁力線極向分量與盤法向的夾角,只有當該夾角大于某臨界值時,物質(zhì)才可能被拋出。區(qū)域Ⅰ和區(qū)域Ⅱ分別為視界和吸積盤上的開放磁力線區(qū)域,它們分別與BZ 和BP 過程相關(guān)。區(qū)域Ⅲ是連接黑洞和盤的閉合磁力線區(qū)域,它對應MC 過程。區(qū)域Ⅰ被限制在θs的范圍內(nèi),區(qū)域Ⅱ被限制在rSIMC~rout之間,區(qū)域Ⅲ被限制在視界線處的θs~θL之間和盤上的rms~rSIMC之間。

    2.3 吸積噴流模型

    經(jīng)過多年的完善,BZ 過程和BP 過程已經(jīng)可以作為通用模型解釋許多天體物理現(xiàn)象。然而隨著觀測手段和觀測精度的提升,單一的模型并不能解釋觀測的結(jié)果,并且計算機模擬結(jié)果表明,兩個過程都在驅(qū)動噴流過程中扮演了重要角色,因此Meier 等人[17]提出了改進的BZ 模型,即大尺度磁場穿過能層外的吸積盤和能層內(nèi)的旋轉(zhuǎn)等離 子體驅(qū)動噴流的模型。為了進一步解決黑洞X 射線雙星多波段輻射等問題,Yuan 等人[18]基于BP 過程提出了吸積噴流模型。這種噴流與BP 過程產(chǎn)生的噴流類似,也是準相對論性的物質(zhì)主導的噴流。由于它是從盤中汲取能量,因此產(chǎn)生這種噴流甚至不需要一個旋轉(zhuǎn)的黑洞[19]。吸積噴流模型主要由三部分組成:外部的SSD 盤、內(nèi)部的熱吸積流和噴流。一部分吸積流在臨近黑洞的地方轉(zhuǎn)變成了垂直于盤的噴流,并產(chǎn)生紅外和射電波段的輻射;外部SSD 盤產(chǎn)生紫外輻射;內(nèi)部的熱吸積流則產(chǎn)生強的X 射線輻射。這種熱吸積流有兩種存在形式,即低于臨界吸積率的徑移主導的吸積流(advection-dominated accretion flow, ADAF)[20]和高于臨界吸積率的發(fā)光熱吸積流(luminous hot accretion flow, LHAF)[21]。LHAF 按照吸積率的高低可分為Ⅰ型和Ⅱ型。人們通過計算發(fā)現(xiàn),熱吸積流具體以哪種形式存在取決于星風指數(shù)s、粘滯系數(shù)α、等離子體參數(shù)β、湍流耗散系數(shù)δ以及描述噴流與吸積流耦合的參數(shù)η。天鵝座X-1是一個很好的例子,該系統(tǒng)有著很高的吸積率和相對低的粘滯系數(shù)。按照吸積噴流模型,其熱流是Ⅱ型LHAF 形式的,這也與觀測結(jié)果相吻合。相對于普通的吸積盤-噴流系統(tǒng),M87在某些波段(例如射電波段)有著明顯的結(jié)構(gòu),能夠清晰地觀測到其射電核心及其兩側(cè)的噴流,這對于研究各個理論模型的正確性有著重要意義[22]。于照蘢[23]利用耦合的ADAF 噴流模型擬合了M87 的射電和X 射線譜,發(fā)現(xiàn)M87 的X 射線輻射主要來自噴流的同步輻射,而不是ADAF。當然,在噴流的產(chǎn)生和準直性傳播過程中,磁場的作用仍然是不可忽視的。磁場控制著吸積盤的角動量,并為噴流的產(chǎn)生及傳播提供動力[24]。貫穿盤的磁場與盤中的等離子體相作用,產(chǎn)生沿著盤面的磁流體動力學(magnetohydrodynamics, MHD) 坡印亭通量,這種能量通量可以進一步轉(zhuǎn)化為動能,產(chǎn)生噴流。由于噴流密度比盤密度低很多,故在此模型中噴流可以達到很高的速度。

    3 吸積盤噴流天體分類

    3.1 原初恒星體

    原恒星是在星際巨分子云緩慢收縮時誕生的早期恒星。原恒星外圍有一個包圍著它的氣體盤,即原行星盤。按多波段光譜能量分布來分類,原初恒星體主要有紅外類型Ⅰ, Ⅱ,Ⅲ, Ⅳ四種類型。這些類型與年輕恒星的演化階段密切相關(guān)。紅外類型Ⅰ包含一顆原恒星,紅外類型Ⅱ與經(jīng)典金牛座T 型星相關(guān),紅外類型Ⅲ與弱線金牛座T 型星相關(guān),紅外類型Ⅳ對應零齡主序星。原行星盤存在于前三類天體系統(tǒng)中,盤上不同區(qū)域?qū)煌墓庾V區(qū)域。與微類星體和活動星系核的盤內(nèi)區(qū)域主要受輻射壓力不同,原行星盤內(nèi)區(qū)域主要受氣體壓力,能量可能從類似于紅外Ⅱ型的經(jīng)典金牛座T 型星系統(tǒng)中通過磁重聯(lián)被提取出來[25]。

    Ellerbroek 等人[26]通過甚大陣(Very Large Array, VLA)發(fā)現(xiàn)了一個擁有復雜吸積盤噴流系統(tǒng)的高吸積率中等質(zhì)量原初恒星體08576nr292。它的星周盤結(jié)構(gòu)很可能是雙峰發(fā)射線產(chǎn)生的原因,MgⅠ和FeⅠ通常有著小的譜峰分裂,這與開普勒盤外圍結(jié)構(gòu)相一致。而HeⅠ和FeⅡ通常有著明顯的譜峰分裂,它們可能產(chǎn)生于盤內(nèi)旋轉(zhuǎn)速度更快的區(qū)域。而一些單峰發(fā)射線則與盤上方旋轉(zhuǎn)的盤風以及準直的噴流有關(guān)(見圖2)。

    圖2 原初恒星體吸積盤-噴流系統(tǒng)示意圖[26]

    3.2 X射線雙星

    X 射線雙星通常由一顆普通恒星和一顆致密天體(如中子星或黑洞)組成。低質(zhì)量X 射線雙星中,伴星通過洛希瓣溢流向致密星傳遞物質(zhì)。一些大質(zhì)量X 射線雙星主要通過致密星周圍的吸積盤傳遞物質(zhì)和角動量。在一些X 射線雙星系統(tǒng)中,如果致密星為黑洞,則該系統(tǒng)被稱為黑洞X 射線雙星。在觀測上,同一個黑洞雙星系統(tǒng)常常表現(xiàn)出不同的態(tài),如硬態(tài)、軟態(tài)和甚高態(tài)。硬態(tài)時系統(tǒng)常常伴隨著準相對論性的連續(xù)性噴流,而在黑洞從硬態(tài)轉(zhuǎn)變到軟態(tài)的過程中,系統(tǒng)則會產(chǎn)生相對論性的間歇性噴流[27]。黑洞X 射線雙星的噴流產(chǎn)生于黑洞附近。在靠近黑洞轉(zhuǎn)軸的區(qū)域,磁場呈現(xiàn)出大尺度上規(guī)則的極向結(jié)構(gòu),并被限制在一個漏斗狀的結(jié)構(gòu)中。這一結(jié)論被MHD 模擬結(jié)果所證實。圖3 為吸積盤-噴流系統(tǒng)的三維密度圖[28]。在這個模型中,反向噴流通過磁場線與視界相連,也就是說,噴流通過BZ 機制來提取黑洞旋轉(zhuǎn)能。這種噴流被稱為相對論性噴流(相對論因子γ>10)。這樣的結(jié)構(gòu)也對噴流的準直起到了重要作用。同時,隨著時間的演化,盤的傾斜程度不斷變化,而且噴流的角度也隨著盤一起變化。這也與黑洞X 射線雙星的觀測事實相符合。

    圖3 黑洞吸積盤-噴流系統(tǒng)三維MHD 模擬密度圖[28]

    3.3 活動星系核

    活動星系核是活動星系的核心,其光度比普通星系強得多,可以在很多波段發(fā)射出很強的電磁波。黑洞周圍氣體的吸積是活動星系核高速噴流的驅(qū)動來源。黑洞-吸積盤模型是被廣泛承認的活動星系核標準模型。根據(jù)該模型,活動星系核的中心是一個超大質(zhì)量黑洞,黑洞周圍的氣體在引力的作用下朝黑洞下落。由于吸積氣體具有角動量,因此,它們便形成了一個圍繞黑洞的吸積盤。吸積盤中的氣體會被加熱到很高的溫度,并不斷損失角動量,因此,它們逐漸下落到黑洞中央,同時釋放出巨大的引力能,并以電磁波的形式輻射出來,從而在吸積盤兩側(cè)沿法線方向形成高速噴流。在此過程中,磁場扮演了很重要的角色。當噴流的方向與觀測者視線方向夾角很小時,觀測者還能觀測到視超光速現(xiàn)象。

    Kato 和Mineshige[29]通過MHD 模擬發(fā)現(xiàn),黑洞吸積盤系統(tǒng)中的噴流并不是單一成分,磁場結(jié)構(gòu)也不是單一成分。如圖4 所示,在黑洞視界附近,纏繞著很強的極向磁場,在極向磁場之外的內(nèi)盤區(qū)域又包裹著環(huán)形磁場。除了吸積盤軸附近坡印亭能流主導的相對論性噴流外,他們還發(fā)現(xiàn)了起源于吸積盤的物質(zhì)主導的準相對論性噴流,這種模型被稱為磁塔模型。觀測結(jié)果表明,這種模型或許能更好地解釋噴流多波段觀測結(jié)果。

    圖4 磁塔模型示意圖[29]

    Sadowski 和Narayan[30]利用輻射磁流體力學程序KORAL 模擬了無旋轉(zhuǎn)超大質(zhì)量黑洞的吸積過程。由于黑洞沒有旋轉(zhuǎn),因此該模擬產(chǎn)生的噴流與BZ 過程無直接關(guān)系。圖5 為高吸積率和幾何厚盤情況下吸積盤系統(tǒng)的密度和亮度分布圖。吸積流以幾何厚盤的形式存在,圍繞著系統(tǒng)旋轉(zhuǎn)軸有兩個對稱的漏斗狀結(jié)構(gòu),其中的大量氣體使之達到光學厚水平。盤中磁旋轉(zhuǎn)不穩(wěn)定性帶來的能量耗散使氣體加熱,并使其在氣體冷卻過程中產(chǎn)生輻射。來自盤上的輻射首先流入漏斗狀結(jié)構(gòu)中,所產(chǎn)生的輻射壓沿軸線加速其中的物質(zhì)至準相對論性,形成噴流,同時漏斗狀結(jié)構(gòu)也逐漸變得透明。從能量上來看,模擬結(jié)果與觀測一致。因為從特定方向可以觀測到很強的超愛丁頓光度現(xiàn)象,所以這個模擬可以很好地解釋觀測到的超亮X射線源等問題。

    圖5 無旋黑洞吸積盤系統(tǒng)密度和亮度分布圖[30]

    3.4 雙中子星

    雙中子星系統(tǒng)或中子星-恒星質(zhì)量黑洞系統(tǒng)中,天體的并合可導致形成環(huán)繞中心黑洞旋轉(zhuǎn)的吸積盤[31]。Kawamura 等人[32]通過相對論性模擬程序LORENE 模擬了質(zhì)量比為0.8 的兩顆中子星的并合過程。圖6a)為中子雙星并合產(chǎn)生大質(zhì)量黑洞后26 ms 時形成的吸積盤結(jié)構(gòu)圖;圖6b)為兩顆等質(zhì)量中子星并合后35.1 ms 時的磁場結(jié)構(gòu)圖,其中,黃色和藍色分別代表密度為108g·cm?3和1010g·cm?3的等密度面,纏卷的彩色線代表磁場線。由圖6b)可以看出盤中的環(huán)形磁場以及極向的漏斗狀結(jié)構(gòu),這與BZ 模型相符。但是,雖然其中有準直性結(jié)構(gòu),但并沒有相對論性的出流,這可能是由于在模擬中未考慮雙星并合時開爾文-亥姆霍茲不穩(wěn)定性(Kelvin-Helmholtz instability, KHI)和并合后磁旋轉(zhuǎn)不穩(wěn)定性(magneto-rotational instability, MRI) 對磁場的放大作用,導致磁場強度不足以觸發(fā)BZ 過程驅(qū)動噴流;也可能是由于模擬時間在毫秒量級,噴流還沒有形成。

    圖6 雙中子星并合密度和磁場模擬結(jié)果[32]

    4 吸積盤噴流的實驗室模擬

    實驗室天體物理作為一門新興的交叉學科,比傳統(tǒng)天體物理有著獨特的優(yōu)勢。隨著實驗室天體物理的發(fā)展和高能量密度設(shè)備的不斷進步,我們已經(jīng)可以在實驗室中利用等離子體的一系列無量綱參數(shù)更加真實地模擬天體物理現(xiàn)象。目前,主要模擬手段有激光驅(qū)動等離子體、Z 箍縮(Z-Pinch)和流體動力學等,在這里我們將主要介紹前兩種方式。

    4.1 激光驅(qū)動等離子體

    由于距離和分辨率的限制,我們對天體的吸積盤了解甚少。但是通過在實驗室恰當排列一組激光器,我們可以在某些方面模擬天體的吸積盤。盤的較差自轉(zhuǎn)、湍動、不規(guī)則的粘滯力等特性共同影響著盤的物理特征。Ryutov[33]利用強激光,制造了接近真實天體情況的高雷諾數(shù)和磁雷諾數(shù)的旋轉(zhuǎn)準平面盤以及漏斗狀極向出流,如圖7a)所示。實驗采用碳氫靶,9 路參數(shù)完全相同的激光同時照射靶,單束能量約為30 J。為了得到帶有角動量的旋轉(zhuǎn)吸積盤,靶的法線方向與半徑方向有一個角度θ。等離子體流在中心區(qū)域重疊,形成旋轉(zhuǎn)的盤,對于更大的θ角,可能會形成環(huán)形的結(jié)構(gòu)。在盤旋轉(zhuǎn)的同時,等離子體團受到壓力的作用開始沿垂直盤面的方向運動,形成如圖7b)所示的一對反向出流。

    值得注意的是,粘滯性對于吸積盤系統(tǒng)有重要影響,低粘度會導致向中心運動的氣體角頻率增加,而高粘度則會導致系統(tǒng)向剛體運動轉(zhuǎn)變,即中心處離心力變小,氣體被推向中心而形成盤狀結(jié)構(gòu)。實驗得到了大于104的雷諾數(shù),這意味著剪切湍流將影響角動量的徑向輸運。通過引入一個極向磁場以及盤中剪切流的放大可以獲得10~100 的磁雷諾數(shù),這對于盤的動力學有著極大影響。通過此實驗發(fā)現(xiàn),雷諾數(shù)和磁雷諾數(shù)可以在很大范圍內(nèi)變化,通過在盤中增加較重的“雜質(zhì)”可以有效降低盤的溫度,從而減小盤的粘性,得到更大的雷諾數(shù)和磁雷諾數(shù)。

    圖7 實驗陣列的幾何結(jié)構(gòu)[33]

    4.2 Z 箍縮驅(qū)動噴流

    與激光驅(qū)動等離子體類似,用Z 箍縮方法也可以驅(qū)動多團等離子體流,并且可以很好地控制多團等離子體的旋轉(zhuǎn)。類似于Ryutov 的激光排列,Bennett 等人[34]利用Z 箍縮方法在MAGPIE 裝置上(電流強度為1.4 MA,持續(xù)時間為250 ns)模擬了旋轉(zhuǎn)的吸積盤結(jié)構(gòu)。實驗設(shè)置如圖8 所示,上下兩個電極之間連接了8 根或16 根導線,組成半徑為2 mm 的圓柱形陣列。該陣列在通電時可以產(chǎn)生全局尺度上的環(huán)向磁場,并在洛倫茲力作用下向內(nèi)箍縮。同時,他們還添加了一個徑向磁場,離化的等離子體在該磁場作用下偏離軸心θ角度向中心運動,最終形成旋轉(zhuǎn)的吸積盤結(jié)構(gòu)。實驗使用的導線為長度4 mm,直徑40μm 的鋁導線。除自發(fā)光成像和分幅相機外,實驗還通過湯姆孫散射探針光診斷等離子體團旋轉(zhuǎn)速度。通過該實驗,他們成功地觀測到了210 ns 時的一個直徑約3 mm 的環(huán)狀結(jié)構(gòu)(如圖8c)所示),并測得等離子體密度大于1019cm?3,盤的旋轉(zhuǎn)可持續(xù)約150 ns。通過測量湯姆孫散射譜線輪廓展寬和中心波長的移動,他們測得等離子體團徑向運動速度約為100 km·s?1,盤旋轉(zhuǎn)速度約為60 km·s?1(馬赫數(shù)約為2),并通過計算得到離子溫度約為60 eV。這個實驗為我們研究可壓縮理想磁流體奠定了基礎(chǔ)。同時,由于盤可存在較長時間,因此,我們可以更直接地研究盤中的不穩(wěn)定性和湍流等的產(chǎn)生和傳播過程。此實驗雷諾數(shù)可達105,磁雷諾數(shù)可達10~100,這將為我們研究MRI 引起的角動量輸運以及諸多磁場主導的作用提供幫助。

    圖8 用Z 箍縮方式模擬盤結(jié)構(gòu)實驗排布及結(jié)果示意圖[34]

    5 展 望

    天體吸積盤-噴流系統(tǒng)是宇宙中廣泛存在的天體系統(tǒng),不同的系統(tǒng)有著不同的形態(tài)結(jié)構(gòu)、磁場結(jié)構(gòu)和形成機制等。用BZ-MC-BP 模型可以基本解釋多數(shù)吸積天體系統(tǒng)的出流問題,但有的問題該模型無法解釋,其他模型也各有利弊。對于盤狀結(jié)構(gòu)的形成原因,磁場對于盤的旋轉(zhuǎn)和形態(tài)的影響,以及噴流的產(chǎn)生原因等還需要進一步討論。利用實驗室天體物理方法,人們可以在短時標小尺度上模擬出盤狀結(jié)構(gòu),從而可以近距離可控地研究盤的形成、盤內(nèi)不穩(wěn)定性、噴流的產(chǎn)生、準直性傳播等問題。近年來,中國實驗室天體物理的研究也在蓬勃發(fā)展。隨著神光三號等大型激光器的建成,在Ryutov 的激光陣列和Bennett 等人的Z 箍縮導線陣列排布的基礎(chǔ)上,人們可以利用中國高能量密度設(shè)備來模擬多團等離子體“首尾”碰撞形成盤狀結(jié)構(gòu)的過程,從而使這類吸積盤噴流實驗有望實現(xiàn)。通過多束激光與V 形靶作用產(chǎn)生多團準直性較好的等離子體,并使等離子體團向中心運動,最終可使其因碰撞擠壓形成類吸積盤結(jié)構(gòu)。希望這種盤狀結(jié)構(gòu)可以通過尺度變換與特定天體的吸積盤相對應,以進一步解決吸積盤的吸積過程、形成機制、盤內(nèi)磁場結(jié)構(gòu)、磁場在吸積過程中扮演的角色以及磁場對于出流的影響等問題。

    猜你喜歡
    磁力線噴流天體
    多根無限長恒定電流形成的磁力線方程研究
    一種改善磁瓦波形結(jié)構(gòu)的模具設(shè)計方法
    用Maltab實現(xiàn)恒定磁場的可視化
    太陽系中的小天體
    太空探索(2020年10期)2020-10-22 03:59:40
    “慧眼”發(fā)現(xiàn)迄今距離黑洞最近的高速噴流
    測量遙遠天體的秘籍
    一分鐘認識深空天體
    磁致炫光研究
    噴流干擾氣動熱數(shù)值模擬的若干影響因素
    耀變體噴流高能電子譜的形成機制
    丰满少妇做爰视频| 午夜免费男女啪啪视频观看| 真实男女啪啪啪动态图| 精品久久久久久久久久久久久| 人人妻人人澡人人爽人人夜夜 | 小蜜桃在线观看免费完整版高清| 尤物成人国产欧美一区二区三区| 久久精品久久久久久久性| 欧美3d第一页| 亚洲美女搞黄在线观看| 精品久久久久久成人av| 亚洲精品一二三| 大话2 男鬼变身卡| 丰满人妻一区二区三区视频av| 男女国产视频网站| 一本一本综合久久| 精品一区二区三卡| 精品人妻偷拍中文字幕| av一本久久久久| 精品欧美国产一区二区三| 极品少妇高潮喷水抽搐| 看十八女毛片水多多多| 免费播放大片免费观看视频在线观看| 免费大片18禁| 久久久久久久大尺度免费视频| 国产综合懂色| 三级男女做爰猛烈吃奶摸视频| 91午夜精品亚洲一区二区三区| 我的女老师完整版在线观看| 99久久精品国产国产毛片| 欧美成人午夜免费资源| 黄色配什么色好看| 日韩不卡一区二区三区视频在线| 青春草亚洲视频在线观看| 人妻少妇偷人精品九色| kizo精华| 国产真实伦视频高清在线观看| 国产亚洲午夜精品一区二区久久 | 久久精品国产亚洲av天美| 成人高潮视频无遮挡免费网站| 亚洲精品乱码久久久久久按摩| 美女黄网站色视频| 男人爽女人下面视频在线观看| 最近最新中文字幕大全电影3| 看免费成人av毛片| 免费观看a级毛片全部| 国产综合精华液| 美女高潮的动态| 免费av观看视频| 中文字幕av成人在线电影| 精品久久久久久久久av| 国产视频内射| 少妇熟女aⅴ在线视频| 亚洲国产精品sss在线观看| 一个人免费在线观看电影| 天堂俺去俺来也www色官网 | av免费在线看不卡| av国产久精品久网站免费入址| 晚上一个人看的免费电影| 99re6热这里在线精品视频| 日韩精品青青久久久久久| 精品久久久久久久久av| 色视频www国产| 久久久久免费精品人妻一区二区| freevideosex欧美| av在线播放精品| 国产精品蜜桃在线观看| 国内精品一区二区在线观看| 亚洲成人中文字幕在线播放| 校园人妻丝袜中文字幕| 国产一级毛片在线| av天堂中文字幕网| 在线免费观看的www视频| 大香蕉97超碰在线| a级一级毛片免费在线观看| 久久久久久九九精品二区国产| 亚洲av福利一区| 久久久久网色| 亚洲av中文字字幕乱码综合| 午夜福利在线在线| 亚洲欧美日韩卡通动漫| 久久久久久久久久久丰满| 国产高清不卡午夜福利| 亚洲精品色激情综合| 精品欧美国产一区二区三| 国产欧美日韩精品一区二区| 成人毛片a级毛片在线播放| 2018国产大陆天天弄谢| 国产视频内射| 三级经典国产精品| 欧美精品国产亚洲| av女优亚洲男人天堂| 久久这里有精品视频免费| 久久热精品热| 日日干狠狠操夜夜爽| 午夜福利成人在线免费观看| 久久久久久久久久人人人人人人| 国产精品国产三级国产av玫瑰| 欧美xxxx黑人xx丫x性爽| 午夜久久久久精精品| 日本免费在线观看一区| 国产亚洲最大av| 青春草国产在线视频| 亚洲在久久综合| 黄片wwwwww| 六月丁香七月| 日韩一区二区视频免费看| 亚洲精品乱码久久久久久按摩| 国产亚洲最大av| 超碰97精品在线观看| 岛国毛片在线播放| 国产精品人妻久久久影院| 久久久久网色| 日本熟妇午夜| 日韩欧美三级三区| a级毛色黄片| 亚洲人成网站在线播| 国产精品一区二区三区四区免费观看| 九九久久精品国产亚洲av麻豆| 亚洲国产精品成人综合色| 亚洲美女搞黄在线观看| 在线观看av片永久免费下载| 久久精品久久精品一区二区三区| av线在线观看网站| 免费观看av网站的网址| 国产精品精品国产色婷婷| 久热久热在线精品观看| 亚洲成人av在线免费| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 视频中文字幕在线观看| 色播亚洲综合网| 91精品伊人久久大香线蕉| 一本一本综合久久| 成人亚洲精品av一区二区| 成人欧美大片| 在线天堂最新版资源| 亚洲精品成人久久久久久| 又粗又硬又长又爽又黄的视频| 在线观看美女被高潮喷水网站| 18禁裸乳无遮挡免费网站照片| 99热全是精品| 丝袜喷水一区| 国产av在哪里看| 丝袜美腿在线中文| 午夜久久久久精精品| 久久热精品热| 六月丁香七月| 亚洲三级黄色毛片| 免费看a级黄色片| 日日摸夜夜添夜夜爱| 亚洲天堂国产精品一区在线| 纵有疾风起免费观看全集完整版 | 极品少妇高潮喷水抽搐| 免费看a级黄色片| 女人十人毛片免费观看3o分钟| 亚洲久久久久久中文字幕| 日本-黄色视频高清免费观看| 777米奇影视久久| 三级毛片av免费| 免费看光身美女| 99久久中文字幕三级久久日本| 婷婷色综合www| 少妇的逼水好多| 五月伊人婷婷丁香| 乱系列少妇在线播放| 特级一级黄色大片| 国产一区亚洲一区在线观看| 蜜桃亚洲精品一区二区三区| 国产女主播在线喷水免费视频网站 | 成人二区视频| 三级男女做爰猛烈吃奶摸视频| 亚洲18禁久久av| 最近2019中文字幕mv第一页| 老司机影院毛片| 国国产精品蜜臀av免费| 一级av片app| 人体艺术视频欧美日本| 天美传媒精品一区二区| 国产探花在线观看一区二区| 男人舔女人下体高潮全视频| 好男人在线观看高清免费视频| 80岁老熟妇乱子伦牲交| 看非洲黑人一级黄片| 毛片女人毛片| 伊人久久精品亚洲午夜| 国产一区二区亚洲精品在线观看| 欧美 日韩 精品 国产| 永久网站在线| 免费大片18禁| 日本一本二区三区精品| 亚洲真实伦在线观看| 大片免费播放器 马上看| 亚洲丝袜综合中文字幕| 亚洲伊人久久精品综合| 国产美女午夜福利| av.在线天堂| 舔av片在线| 亚洲图色成人| 欧美日韩在线观看h| 久久国内精品自在自线图片| 久久鲁丝午夜福利片| 亚洲精品色激情综合| 又爽又黄无遮挡网站| 高清午夜精品一区二区三区| 91在线精品国自产拍蜜月| 黄色欧美视频在线观看| av国产免费在线观看| 大香蕉97超碰在线| 成人亚洲精品一区在线观看 | av网站免费在线观看视频 | a级一级毛片免费在线观看| 久久精品国产自在天天线| 久久久久网色| 国产av码专区亚洲av| 91久久精品国产一区二区三区| 干丝袜人妻中文字幕| 日本-黄色视频高清免费观看| 日韩强制内射视频| 狠狠精品人妻久久久久久综合| 国产有黄有色有爽视频| 日日摸夜夜添夜夜爱| 日韩一区二区视频免费看| 日韩不卡一区二区三区视频在线| 一级黄片播放器| 色网站视频免费| 噜噜噜噜噜久久久久久91| 国产在线一区二区三区精| 婷婷色麻豆天堂久久| 午夜福利在线观看吧| 99久久精品一区二区三区| 成人国产麻豆网| 51国产日韩欧美| 99re6热这里在线精品视频| 亚洲国产欧美在线一区| 午夜福利成人在线免费观看| 人人妻人人澡欧美一区二区| 精品欧美国产一区二区三| 国产成人精品福利久久| 青青草视频在线视频观看| 天堂网av新在线| 日本色播在线视频| 韩国av在线不卡| av.在线天堂| 91在线精品国自产拍蜜月| 99久久精品一区二区三区| 菩萨蛮人人尽说江南好唐韦庄| 最近2019中文字幕mv第一页| 精品一区二区三区人妻视频| 91精品伊人久久大香线蕉| 亚洲精华国产精华液的使用体验| 日韩欧美三级三区| 日本三级黄在线观看| 国产黄色免费在线视频| 成人性生交大片免费视频hd| 欧美日韩一区二区视频在线观看视频在线 | 99久久精品热视频| 欧美xxxx性猛交bbbb| 国产高清不卡午夜福利| 少妇猛男粗大的猛烈进出视频 | 亚洲精品乱码久久久久久按摩| 久久韩国三级中文字幕| 欧美最新免费一区二区三区| 国产一区二区亚洲精品在线观看| 伊人久久国产一区二区| 亚洲欧美成人综合另类久久久| 国产麻豆成人av免费视频| 老司机影院毛片| 综合色丁香网| 在线播放无遮挡| 亚洲av成人av| 观看免费一级毛片| 国产精品久久久久久精品电影| 午夜激情欧美在线| 精品人妻偷拍中文字幕| 高清欧美精品videossex| 大香蕉久久网| 搡老妇女老女人老熟妇| 国产亚洲精品av在线| 国产精品一区二区三区四区久久| 国产成人aa在线观看| 国产成人一区二区在线| av国产免费在线观看| 青青草视频在线视频观看| 女的被弄到高潮叫床怎么办| 国产一级毛片七仙女欲春2| 午夜福利视频精品| 女人久久www免费人成看片| 亚洲av电影在线观看一区二区三区 | 日本三级黄在线观看| 2021天堂中文幕一二区在线观| 久久久久九九精品影院| 久久久久久久久久人人人人人人| 一区二区三区四区激情视频| 欧美+日韩+精品| a级毛片免费高清观看在线播放| 久久精品国产鲁丝片午夜精品| 免费看a级黄色片| 能在线免费观看的黄片| 国产女主播在线喷水免费视频网站 | 成年女人在线观看亚洲视频 | 美女内射精品一级片tv| 久久鲁丝午夜福利片| 听说在线观看完整版免费高清| 午夜久久久久精精品| 欧美bdsm另类| 女的被弄到高潮叫床怎么办| 欧美日韩精品成人综合77777| 简卡轻食公司| av.在线天堂| 亚洲精品影视一区二区三区av| 久久久欧美国产精品| 又粗又硬又长又爽又黄的视频| 蜜桃久久精品国产亚洲av| 久久99精品国语久久久| 少妇熟女欧美另类| 国产av在哪里看| 欧美xxxx性猛交bbbb| 永久免费av网站大全| 国产亚洲精品久久久com| 欧美激情国产日韩精品一区| 成人特级av手机在线观看| 精品国内亚洲2022精品成人| 久久人人爽人人爽人人片va| 精品国内亚洲2022精品成人| 国产成人福利小说| 免费人成在线观看视频色| av免费在线看不卡| 国产高清三级在线| 亚洲国产欧美在线一区| 国产淫语在线视频| 国产国拍精品亚洲av在线观看| 亚洲美女搞黄在线观看| 九九久久精品国产亚洲av麻豆| 又爽又黄无遮挡网站| 免费高清在线观看视频在线观看| 国产淫片久久久久久久久| 欧美xxⅹ黑人| 亚洲精品日本国产第一区| 午夜福利在线观看免费完整高清在| 日韩av在线大香蕉| 人妻制服诱惑在线中文字幕| 日韩av在线大香蕉| 亚洲经典国产精华液单| 亚洲精品影视一区二区三区av| 久久精品国产亚洲av涩爱| 久久久精品欧美日韩精品| 美女被艹到高潮喷水动态| 亚洲国产欧美在线一区| 成人欧美大片| 狂野欧美激情性xxxx在线观看| 联通29元200g的流量卡| 一级a做视频免费观看| 天堂网av新在线| 亚洲精品色激情综合| 国产 一区 欧美 日韩| 亚洲国产av新网站| videos熟女内射| 久久久成人免费电影| 久久久久久久国产电影| 日韩制服骚丝袜av| 乱码一卡2卡4卡精品| 九九久久精品国产亚洲av麻豆| 晚上一个人看的免费电影| 日本色播在线视频| 免费黄网站久久成人精品| 精品久久久久久久久亚洲| 国产精品伦人一区二区| 最新中文字幕久久久久| 精品人妻一区二区三区麻豆| 亚洲一级一片aⅴ在线观看| 老师上课跳d突然被开到最大视频| 3wmmmm亚洲av在线观看| 欧美成人午夜免费资源| 最近最新中文字幕免费大全7| 熟妇人妻不卡中文字幕| 国产精品久久久久久精品电影| 午夜精品国产一区二区电影 | av免费在线看不卡| 黄色欧美视频在线观看| 色综合站精品国产| 久久亚洲国产成人精品v| 我要看日韩黄色一级片| 麻豆成人av视频| 嫩草影院新地址| 免费大片黄手机在线观看| 欧美zozozo另类| 国产成人精品福利久久| 最后的刺客免费高清国语| 国产av码专区亚洲av| 久久国内精品自在自线图片| 22中文网久久字幕| 国产中年淑女户外野战色| 亚洲精品视频女| 非洲黑人性xxxx精品又粗又长| 久久久亚洲精品成人影院| 亚洲精品成人av观看孕妇| 人人妻人人澡欧美一区二区| 国产精品蜜桃在线观看| 国产精品久久久久久精品电影| 亚洲图色成人| 精品午夜福利在线看| 免费少妇av软件| 麻豆国产97在线/欧美| 秋霞伦理黄片| 青青草视频在线视频观看| 日本熟妇午夜| 一区二区三区免费毛片| 成年女人看的毛片在线观看| 国产麻豆成人av免费视频| 国产av在哪里看| 免费看日本二区| 久久精品久久久久久噜噜老黄| 欧美成人精品欧美一级黄| 国内揄拍国产精品人妻在线| 国产乱人偷精品视频| 99久久中文字幕三级久久日本| 中文在线观看免费www的网站| 久久久久久伊人网av| 国产成人精品一,二区| 又粗又硬又长又爽又黄的视频| 午夜精品一区二区三区免费看| av又黄又爽大尺度在线免费看| 一级爰片在线观看| or卡值多少钱| 不卡视频在线观看欧美| 精品久久久久久久末码| 欧美xxⅹ黑人| 一级黄片播放器| 日本色播在线视频| 中文字幕人妻熟人妻熟丝袜美| 国产 亚洲一区二区三区 | 午夜精品国产一区二区电影 | 床上黄色一级片| 国产午夜精品一二区理论片| 日本免费a在线| 国产成人精品福利久久| 国产亚洲av嫩草精品影院| 亚洲国产高清在线一区二区三| 国内精品宾馆在线| 亚洲精品久久久久久婷婷小说| 麻豆乱淫一区二区| eeuss影院久久| 极品少妇高潮喷水抽搐| 国产一区二区亚洲精品在线观看| 日韩成人伦理影院| 欧美成人a在线观看| 午夜福利在线观看免费完整高清在| 亚洲av免费高清在线观看| 亚洲激情五月婷婷啪啪| 亚洲精品456在线播放app| 成人国产麻豆网| av在线老鸭窝| 欧美日本视频| 观看美女的网站| 人人妻人人澡人人爽人人夜夜 | 国产精品99久久久久久久久| 好男人视频免费观看在线| 国产精品99久久久久久久久| 国产午夜精品久久久久久一区二区三区| 性插视频无遮挡在线免费观看| 成年女人在线观看亚洲视频 | 精品一区在线观看国产| 午夜精品一区二区三区免费看| 特大巨黑吊av在线直播| 在线观看av片永久免费下载| 禁无遮挡网站| 国内精品一区二区在线观看| 欧美一区二区亚洲| 看黄色毛片网站| 国产男人的电影天堂91| 亚洲天堂国产精品一区在线| 亚洲av成人精品一二三区| 91精品国产九色| 亚洲成人精品中文字幕电影| 久久精品夜色国产| 中文天堂在线官网| 五月天丁香电影| 一级毛片黄色毛片免费观看视频| 视频中文字幕在线观看| 久久久久免费精品人妻一区二区| 高清av免费在线| 精品国内亚洲2022精品成人| 最近最新中文字幕免费大全7| 欧美极品一区二区三区四区| 日本与韩国留学比较| 麻豆精品久久久久久蜜桃| 五月伊人婷婷丁香| 国产成人精品一,二区| 国产精品精品国产色婷婷| 人人妻人人澡人人爽人人夜夜 | 国产欧美日韩精品一区二区| 日日干狠狠操夜夜爽| 国产精品久久久久久久电影| 人妻系列 视频| 日韩精品有码人妻一区| 日韩视频在线欧美| 久久99热这里只频精品6学生| 人人妻人人澡人人爽人人夜夜 | 黄色日韩在线| 熟女人妻精品中文字幕| 日韩伦理黄色片| 中文字幕久久专区| 午夜亚洲福利在线播放| 国产一级毛片七仙女欲春2| 色视频www国产| 一级毛片 在线播放| 国产精品人妻久久久久久| 丝袜美腿在线中文| 国产精品综合久久久久久久免费| 亚洲成人一二三区av| 丰满人妻一区二区三区视频av| 日韩一区二区视频免费看| 国产高清有码在线观看视频| 天堂中文最新版在线下载 | 国产91av在线免费观看| 麻豆国产97在线/欧美| 国内精品一区二区在线观看| 国产老妇伦熟女老妇高清| 久久国产乱子免费精品| 亚洲人成网站高清观看| 超碰av人人做人人爽久久| 免费av不卡在线播放| 午夜激情欧美在线| 亚洲国产精品专区欧美| 在线播放无遮挡| 能在线免费看毛片的网站| 午夜福利在线在线| 国产综合懂色| 国产综合精华液| 免费观看av网站的网址| 一级av片app| 精品午夜福利在线看| 一个人看视频在线观看www免费| 精品久久久精品久久久| 日本与韩国留学比较| 丝袜喷水一区| 日本av手机在线免费观看| 大话2 男鬼变身卡| 中文字幕亚洲精品专区| 国产成人福利小说| 亚洲av二区三区四区| 少妇人妻精品综合一区二区| 美女国产视频在线观看| 久久亚洲国产成人精品v| 午夜免费激情av| 国产黄频视频在线观看| 亚洲人成网站高清观看| 最近中文字幕高清免费大全6| 欧美日韩一区二区视频在线观看视频在线 | 亚洲欧美日韩卡通动漫| 午夜福利高清视频| 亚洲久久久久久中文字幕| 久久久久久久午夜电影| 色网站视频免费| 美女高潮的动态| 国产精品一区二区三区四区久久| 欧美人与善性xxx| 久久久色成人| 日韩精品有码人妻一区| 国产又色又爽无遮挡免| 亚洲va在线va天堂va国产| 免费看a级黄色片| 国产成人a区在线观看| 国产高清有码在线观看视频| 天天躁夜夜躁狠狠久久av| 直男gayav资源| 欧美性猛交╳xxx乱大交人| 国产熟女欧美一区二区| 干丝袜人妻中文字幕| 天堂√8在线中文| 网址你懂的国产日韩在线| 秋霞在线观看毛片| 在线免费观看的www视频| 久久鲁丝午夜福利片| 亚洲欧美一区二区三区国产| 亚洲av男天堂| 欧美人与善性xxx| 在线免费十八禁| 免费播放大片免费观看视频在线观看| 国产淫语在线视频| 18禁在线无遮挡免费观看视频| 久久久亚洲精品成人影院| 日韩欧美国产在线观看| 一级二级三级毛片免费看| 寂寞人妻少妇视频99o| 国产91av在线免费观看| 嫩草影院精品99| 国产一级毛片七仙女欲春2| 中文字幕免费在线视频6| or卡值多少钱| 大话2 男鬼变身卡| 丝瓜视频免费看黄片| 搡老乐熟女国产| 寂寞人妻少妇视频99o| 高清毛片免费看| 偷拍熟女少妇极品色| 日本黄大片高清| 高清毛片免费看| 超碰av人人做人人爽久久| 国产精品综合久久久久久久免费| 99热这里只有精品一区| 我要看日韩黄色一级片| av在线蜜桃| 精品亚洲乱码少妇综合久久| 亚洲av一区综合| ponron亚洲| 观看免费一级毛片| 欧美激情在线99| av在线观看视频网站免费| 国产成人精品一,二区| 性插视频无遮挡在线免费观看| 日日啪夜夜撸| 国产探花在线观看一区二区| videossex国产| 婷婷色综合大香蕉|