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

    外日球?qū)拥膶捘芏坞x子及其與湍動(dòng)的耦合作用

    2021-01-30 07:56:12何建森王玲華宗秋剛
    深空探測學(xué)報(bào) 2020年6期
    關(guān)鍵詞:太陽風(fēng)激波星際

    何建森,林 榮,崔 博,王玲華,宗秋剛

    (1.北京大學(xué) 地球與空間科學(xué)院,100871,北京;2.北京大學(xué) 行星與空間科學(xué)研究中心,100871,北京)

    引 言

    太陽系是目前唯一已知的孕育生命和人類文明的星系。在太陽系的四維時(shí)空坐標(biāo)中,當(dāng)今太陽系的空間組成和結(jié)構(gòu)及其時(shí)間演化,是一個(gè)重要的命題(見圖1)。通過探測外太陽系的小天體(沒有明顯的地質(zhì)過程,被認(rèn)為是太陽系的活化石),可以追溯太陽系的形成歷史和早期特征[1-5]。通過探測被太陽風(fēng)所充斥的日球?qū)拥娜S結(jié)構(gòu)及其動(dòng)態(tài)變化,可以認(rèn)識(shí)太陽風(fēng)與局地星際介質(zhì)的相互作用。當(dāng)今的太陽系不再像太陽系的初期形成的行星星云那樣稠密,而是更接近真空狀態(tài),盡管也充斥著太陽風(fēng)等離子體、能量粒子(太陽能量粒子和宇宙線粒子)、塵埃粒子、小天體(小行星和彗星)、矮行星和行星。

    早期的太陽活動(dòng)(比如T-Tauri階段,可能比現(xiàn)在的太陽活動(dòng)強(qiáng)1 000倍)所吹出的太陽風(fēng),或許可以把原初太陽系星云的物質(zhì)清掃輸運(yùn)到很遠(yuǎn)的地方(比如奧特云所在的距離)[6-7]。所以太陽風(fēng)的研究是一個(gè)認(rèn)識(shí)太陽、日球?qū)印⑻栂导捌渑c銀河系的局地星際介質(zhì)相互作用的重要方面。

    人類要全面認(rèn)識(shí)自身所處的太陽系,就需要像環(huán)游地球那樣,周游太陽系開展全面的探測。在太陽系探測遠(yuǎn)景設(shè)想中,認(rèn)為需要在5個(gè)標(biāo)志性的探測能力方面取得重要的突破:寬能段的粒子探測、寬頻段的電磁場測量、寬波段的光學(xué)觀測、大角度的覆蓋、遠(yuǎn)距離的抵達(dá),從而形成全向發(fā)展的探測能力雷達(dá)圖(見圖2)。

    圖1 太陽系的時(shí)空坐標(biāo)系Fig.1 The Solar system in the space-time coordinate system

    圖2 5個(gè)維度的深空探測能力雷達(dá)圖Fig.2 A radar chart illustrating five aspects of deep space exploration capabilities

    在以往的探測中,每次任務(wù)相對于以往的探測任務(wù)在5個(gè)標(biāo)志性的探測能力的某一(某些)方面都會(huì)有突破性,相應(yīng)的也會(huì)帶來科學(xué)上的發(fā)現(xiàn)。Ulysses飛船首次在日球?qū)又酗w躍太陽極區(qū)探測,發(fā)現(xiàn)高緯極區(qū)太陽風(fēng)與低緯黃道面太陽風(fēng)的差異[8-12]。ACE(Advanced Composition Explorer)飛船在日地L1點(diǎn)暈軌道上首次實(shí)現(xiàn)對太陽風(fēng)多種類、多價(jià)態(tài)離子進(jìn)行寬能譜探測,揭示了雙模態(tài)太陽風(fēng)在離子凍結(jié)溫度、不同第一電離勢的離子豐度、超熱離子能譜指數(shù)等方面的差異[13-16]。磁層多尺度星座 (Magnetospheric MultiScale,MMS)星座衛(wèi)星首次實(shí)現(xiàn)電場矢量的測量,研究湍動(dòng)和重聯(lián)中的關(guān)鍵能量轉(zhuǎn)換機(jī)制成為了可能[17-19]。太陽與日光層探測器 (SOlar and Heliospheric Observatory,SOHO)飛船對太陽風(fēng)源區(qū)進(jìn)行寬波段的觀測,發(fā)現(xiàn)了源區(qū)太陽風(fēng)初始外流形成于色球網(wǎng)絡(luò)結(jié)構(gòu)上方的過渡區(qū)高度[20]?!奥眯姓?號(hào)”(Voyager 1)“旅行者2號(hào)”(Voyager 2)飛船首次實(shí)現(xiàn)100 AU以遠(yuǎn)的探測,先后穿越了終止激波、內(nèi)日鞘區(qū)、日球?qū)禹數(shù)汝P(guān)鍵區(qū)域[21-23]。

    針對遠(yuǎn)距離的太陽系邊際探測深空探測任務(wù),由于飛船服役周期長、服役數(shù)量非常有限,如何在有限的發(fā)射機(jī)會(huì)中、在載荷重量和數(shù)據(jù)傳輸率等限制條件下,充分考慮并發(fā)揮上述5個(gè)方面的探測能力,尤為必要。

    1 寬能段粒子相空間密度譜

    寬能段的粒子探測至關(guān)重要。不同來源的粒子具有不同的能量,不同能量的粒子也經(jīng)歷不同的物理過程(見圖3)。源自太陽大氣的太陽風(fēng)等離子體屬于低能粒子,能量在10 keV以下。雖然單個(gè)粒子的能量較低,但是由于相空間密度高,但是單位體積的動(dòng)能和熱能在能量密度中的貢獻(xiàn)都是非常重要。盡管太陽風(fēng)等離子體粒子在低能范疇,然而太陽風(fēng)等離子體的整體流速卻達(dá)幾百千米每秒,遠(yuǎn)大于聲速(在終止激波上游基本小于100 km/s),屬于超聲速范疇。作為一級近似,超聲速的太陽風(fēng)的動(dòng)壓和局地星際介質(zhì)的熱壓之間的平衡假設(shè),決定了日球頂?shù)娜招木嚯x乃至日球?qū)拥拇笮?。超聲速太陽風(fēng)在日球?qū)禹斍靶枰眲p速,從而產(chǎn)生了終止激波。大尺度的終止激波的存在,為粒子經(jīng)過激波加速成為異常宇宙線(太陽系內(nèi)持續(xù)產(chǎn)生的、能量最高的粒子種類)提供了必要的場所。

    圖3 寬能段的日球?qū)与x子相空間密度譜在不同日心距離處的變化[24]Fig.3 The variation of the phase space density spectrum of the heliosphere ions in the wide energy range at different heliocentric distances[24]

    太陽風(fēng)等離子體的速度分布在一定程度上可以用各向異性的麥?zhǔn)戏植紒肀硎?,存在垂直和平行溫度的差異。通過探測分析垂直溫度、平行溫度、密度以及磁場強(qiáng)度隨日心距離的變化,發(fā)現(xiàn)太陽風(fēng)的膨脹不是(雙)絕熱膨脹的,比如質(zhì)子垂直溫度與磁場比值所表征的磁矩是隨距離增加的,導(dǎo)致質(zhì)子投擲角分布不像絕熱機(jī)制預(yù)期的那樣變得越來越趨近平行或者反平行分布[25-28]。太陽風(fēng)質(zhì)子的非絕熱膨脹也導(dǎo)致多方關(guān)系近似的狀態(tài)方程中的多方指數(shù)γ偏離傳統(tǒng)單原子的定壓/定容比熱的比值(-1.67 < γ < -1)。太陽風(fēng)質(zhì)子的非絕熱膨脹是日球?qū)拥闹匾锢憩F(xiàn)象,其背后的加熱機(jī)制(比如湍動(dòng)耗散加熱)是日球?qū)游锢淼年P(guān)鍵科學(xué)問題和深空探測的重要科學(xué)目標(biāo)[29]。要解決/回答該關(guān)鍵問題,需要針對質(zhì)子相空間密度開展高分辨率和高靈敏度的探測,獲得偏離麥?zhǔn)戏植嫉臄_動(dòng)相空間密度的精確分布,并實(shí)施對應(yīng)的場(電場和磁場)的高質(zhì)量的探測,從而能有效分析場-粒子相關(guān)性耦合所導(dǎo)致的能量轉(zhuǎn)換這一無碰撞等離子體湍動(dòng)耗散的關(guān)鍵物理過程[30]。

    針對背景太陽風(fēng)質(zhì)子和氦離子的探測,根據(jù)其能量范圍,采用的經(jīng)典探測儀器是法拉第杯和靜電分析儀。法拉第杯后端的金屬電流收集板具有抗老化長壽命服役、量程大等特點(diǎn),成功用于多項(xiàng)深空探測任務(wù)中[31]。靜電分析儀對低通量粒子具有高靈敏度探測、高角度分辨率等優(yōu)點(diǎn),適合探測向像風(fēng)這樣方向性比較好、數(shù)密度比較稀薄的無碰撞等離子體。法拉第杯和靜電分析儀組合,可以優(yōu)勢互補(bǔ),最大程度發(fā)揮探測太陽風(fēng)等離子體的能力。

    拾起離子也是日球?qū)訉捘芏瘟W拥闹匾M成。“拾起離子”所謂的“拾起”是指:電荷交換后的新生離子在太陽風(fēng)的電磁場中,會(huì)繞著行星際磁場回旋運(yùn)動(dòng),而且發(fā)生電場漂移,從而具有徑向向外的速度分量,與太陽風(fēng)同向運(yùn)動(dòng),被太陽風(fēng)拾起。拾起離子是太陽風(fēng)等離子體與中性原子(中性原子可來源于侵入的星際介質(zhì)流或者太陽系天體釋放逃逸出來的氣體)相互作用交換電荷的重要產(chǎn)物之一(另一產(chǎn)物是能量中性原子)[32]。針對拾起離子成分的探測,能夠探明星際介質(zhì)的(輕)元素組成和豐度(比如H,He,O,Ar,N,Ne等元素)。

    拾起效率隨日心距離而增大:①在內(nèi)日球?qū)?,行星際磁場方向與太陽風(fēng)速度方向近似平行,不利于新生離子被拾起,拾起離子的流速和溫度也遠(yuǎn)低于背景太陽風(fēng)離子的流速和溫度;②在外日球?qū)?,行星際磁場方向與太陽風(fēng)速度方向近似垂直,新生離子容易被拾起,拾起離子的流速與背景太陽風(fēng)流速相當(dāng),拾起離子的等效熱速度由于與太陽風(fēng)流速同數(shù)量級而明顯高于背景太陽風(fēng)離子熱速度。隨著太陽風(fēng)的向外膨脹,其與中性原子發(fā)生電荷交換的次數(shù)越來越多,所產(chǎn)生和攜帶的拾起離子也越來越多,導(dǎo)致拾起離子對外日球?qū)犹栵L(fēng)的改造也越來越明顯,比如對太陽風(fēng)熱壓的貢獻(xiàn)比重變得非常顯著[33]。當(dāng)然拾起離子與背景太陽風(fēng)等離子體的混合,離不開湍動(dòng)的作用。如果沒有湍動(dòng),拾起離子無法在速度相空間中散射,無法有效和背景太陽風(fēng)等離子體混合起來進(jìn)而影響外日球?qū)犹栵L(fēng)的狀態(tài)。有了拾起離子及其受湍動(dòng)的散射熱化,才能導(dǎo)致太陽風(fēng)的微弱減速和持續(xù)加熱。

    散射拾起離子的湍動(dòng)可伴隨背景太陽風(fēng)的湍動(dòng)由拾起離子自激發(fā)產(chǎn)生的波湍動(dòng)。源自星際介質(zhì)的拾起離子自從1985年被動(dòng)態(tài)磁層粒子追蹤探測器(Active Magnetospheric Particle Tracer Explorers,AMPTE)上的超熱能量離子分析儀(SUprathermaL Energy Ionic Charge Analyzer,SULEICA)首次探測到之后[33],其受湍動(dòng)散射形成各向同性殼層分布需要多長時(shí)間(即散射效率)一直備受關(guān)注。經(jīng)過Ulysses飛船上的太陽風(fēng)離子成分譜儀(Solar Wind Ion Composition Spectromete,SWICS)的探測,發(fā)現(xiàn)拾起離子的速度分布沒有預(yù)想的那樣各向同性,而是仍保留各向異性的特征[24],說明湍動(dòng)的散射效率低于預(yù)期,需要更長的散射平均自由程。湍動(dòng)的散射效率取決于湍動(dòng)的狀態(tài):2D結(jié)構(gòu)成分所主導(dǎo)的湍動(dòng)的散射效率低于slab波動(dòng)成分所主導(dǎo)的湍動(dòng)的散射效率[33]。拾起離子自激發(fā)的波動(dòng)主要表現(xiàn)為多普勒頻移的離子回旋波,該波動(dòng)也證明磁強(qiáng)計(jì)可以間接探測拾起離子的活動(dòng)[23]。

    超熱離子是能量介于10 keV和1 MeV之間的離子,是太陽風(fēng)等離子體離子和拾起離子能化的產(chǎn)物,是進(jìn)一步加速成為高能異常宇宙線的種子粒子,在低能粒子和高能粒子之間起到重要的橋梁作用。超熱離子如何由低能離子能化產(chǎn)生,產(chǎn)生后所形成的冪律譜的譜指數(shù)有何特征,譜指數(shù)是否普適(比如相空間密度譜隨速度變化的譜指數(shù)是否在-5附近)[15]?這些是有關(guān)超熱離子的關(guān)鍵問題,也是探測超熱離子的主要?jiǎng)訖C(jī)。針對日球?qū)悠毡榇嬖诘某瑹犭x子,湍動(dòng)(特別是可壓縮的湍動(dòng))被認(rèn)為是低能離子能化形成超熱離子的主要驅(qū)動(dòng)源,湍動(dòng)的磁場和電場導(dǎo)致離子在速度相空間中有效擴(kuò)散,從而填充速度相空間中的高速區(qū)域。所以湍動(dòng)場的測量是認(rèn)識(shí)超熱離子形成之謎的探測基礎(chǔ)。

    2 雙模態(tài)的太陽風(fēng)特征

    太陽系邊際探測能夠在太陽系時(shí)空坐標(biāo)系中獲得對日球?qū)拥膽B(tài)勢感知(見圖4)。像“旅行者1號(hào)”和“旅行者2號(hào)”飛船那樣遠(yuǎn)距離的深空旅行,通過時(shí)間換空間的方式,來獲得日球?qū)訝顟B(tài)變量隨全日心距離的變化剖面[34]。這個(gè)日心距離剖面不是單純某個(gè)時(shí)刻的空間剖面,而是記錄了日球?qū)酉到y(tǒng)對其內(nèi)外邊界活動(dòng)變化的響應(yīng)的信息,是一個(gè)時(shí)空復(fù)合條件下的狀態(tài)感知(見圖5)。所以面對探測數(shù)據(jù)進(jìn)行解讀時(shí),必須區(qū)分出哪些是空間變化(太陽風(fēng)向外膨脹),哪些是時(shí)間變化引起的(太陽活動(dòng)高低年相位和太陽自轉(zhuǎn)效應(yīng))。

    圖4 “旅行者2號(hào)”飛船的日心距離、緯度隨時(shí)間的變化Fig.4 The change of distance and latitude of Voyager 2 as well as sunspot number over time

    為了實(shí)現(xiàn)有效的區(qū)分,離不開聯(lián)合其他觀測計(jì)劃的支持,SOHO飛船、STEREO飛船以及我國的先進(jìn)天基太陽天文臺(tái)衛(wèi)星(Advanced Space-based Solar Observatory,ASO-S)[35]對太陽風(fēng)源區(qū)全球分布的監(jiān)測以及全球源區(qū)隨時(shí)間變化的監(jiān)測,離不開在行星際空間不同日心距離處(比如行星探測衛(wèi)星)對太陽風(fēng)的監(jiān)測。

    圖5 “旅行者1號(hào)”飛船上CRS載荷探測不同能檔高能粒子的微分通量密度隨飛行時(shí)間的變化Fig.5 Time variation of the differential flux density of energetic particles in different energy bands as detected by CRS onboard Voyager 1

    圖6 日球?qū)犹栵L(fēng)及其湍動(dòng)的雙模態(tài)Fig.6 The dual modes of the solar wind and its turbulence in the heliosphere

    關(guān)于日球?qū)犹栵L(fēng)的探測獲得了一個(gè)非?;镜苤匾恼J(rèn)知:雙模態(tài)(bimodal states)。在源區(qū),針對日冕輻射強(qiáng)度的明暗、日冕磁場的拓?fù)湮恍?,可以把太陽風(fēng)源區(qū)分為兩大類(冕洞區(qū)、寧靜區(qū)/活動(dòng)區(qū))。在行星際空間,針對太陽風(fēng)的流速、數(shù)密度和溫度等特征,分為高速太陽風(fēng)和低速太陽風(fēng)[36]。雙模態(tài)的行星際太陽風(fēng)所伴隨的湍動(dòng)經(jīng)常也是具有雙模態(tài)(見圖6):高速流經(jīng)常伴隨強(qiáng)的阿爾芬性湍動(dòng),橫向的擾動(dòng)速度和擾動(dòng)磁場具有很好的相關(guān)性,磁場強(qiáng)度的起伏則不是很大,Els?sser變量有主次之分,主變量的功率譜具有明顯的三段譜的特征(即含能區(qū)、慣性區(qū)、耗散區(qū)),次變量的功率譜在低頻段沒有呈現(xiàn)出明顯的含能區(qū)-慣性區(qū)的拐點(diǎn)特征,可能說明主變量對次變量的串級乃至功率譜輪廓的影響很大;低速流伴隨可壓縮性的湍動(dòng),湍動(dòng)里不時(shí)表現(xiàn)出壓力平衡結(jié)構(gòu)的特征,即相鄰區(qū)域的總壓力(熱壓+磁壓)基本相等[37]。

    高低速流及其伴隨湍動(dòng)的雙模態(tài),說明太陽風(fēng)起源的多樣性與復(fù)雜性,無法用以前的一維流管中阿爾芬波驅(qū)動(dòng)太陽風(fēng)的模型來簡單地統(tǒng)一解釋,需要拓展到三維模型中,而且需要自洽考慮太陽風(fēng)等離子體熱動(dòng)力學(xué)演化和波湍動(dòng)演化相互耦合的過程。無獨(dú)有偶,超熱離子的能譜也具有雙模態(tài)的特征:高速流中的超熱離子能譜偏軟(較少的離子被能化),而低速流中的超熱離子能譜則相對較硬(更多的離子被能化)(見圖7)[38]。所以超熱離子能譜沒有出現(xiàn)有些學(xué)者預(yù)期的冪律譜指數(shù)唯一性的特征。究其原因,高低速流的超熱離子能譜的雙模態(tài)可能源于湍動(dòng)的雙模態(tài),即低速流中的可壓縮湍動(dòng)(相對于高速流中的不可壓湍動(dòng))對于離子能化形成超熱離子更為有效。似乎高速流中的阿爾芬性不可壓湍動(dòng)的串級耗散,只對太陽風(fēng)等離子體的熱離子有明顯的加熱效應(yīng),而對進(jìn)一步能化成超熱離子則顯得有些“無能為力”。搞清楚湍動(dòng)與不同能量粒子之間的能化關(guān)系,是構(gòu)建日球?qū)游锢淼幕弧?/p>

    圖7 寬能段離子的相空間密度譜在低速流和高速流中的差異Fig.7 Difference between the phase space density spectra of the wideenergy-band ions in the low-speed and high-speed flows

    雙模態(tài)的太陽風(fēng)及其湍動(dòng)在外日球?qū)又惺侨绾窝莼??由于相鄰的共轉(zhuǎn)相互作用區(qū)邊界上前向激波和后向激波相互碰撞而滲透到對方區(qū)域中,導(dǎo)致前向激波不僅壓縮本CIR(Co-rotating Interaction Region)前面的低速流區(qū),而且侵入前面的CIR去壓縮其中的高速流區(qū),后向激波不僅壓縮本CIR后面的低速流區(qū),而且侵入后面的CIR去壓縮其中的低速流區(qū)。雙模態(tài)的太陽風(fēng)及其湍動(dòng)是否隨著不同CIR的融合(產(chǎn)生Merged Interaction Region,MIR)而導(dǎo)致界限模糊,目前還缺乏專門的定量研究[34]。雙模態(tài)的太陽風(fēng)及其湍動(dòng)對終止激波及其下游日鞘層的影響是否也有差別?由于超熱離子是加速形成異常宇宙線的種子粒子,超熱離子在高速流和低速流中的能譜又不一樣,這是否意味著低速流中的超熱離子對形成異常宇宙線的貢獻(xiàn)更大?這些問題都是需要通過探測外日球?qū)觼韺?shí)現(xiàn)。

    3 三維日球?qū)咏Y(jié)構(gòu)、拾起離子產(chǎn)生與演化、湍動(dòng)與超熱離子形成

    對比理想的探測能力雷達(dá)圖和現(xiàn)實(shí)中的探測成就,發(fā)現(xiàn)有幾個(gè)可以取得突破的地方。大角度的覆蓋和遠(yuǎn)距離的抵達(dá),這兩個(gè)條件目前沒有同時(shí)滿足?!奥眯姓?號(hào)”和“旅行者2號(hào)”飛船都是飛向日球?qū)禹敱羌夥轿桓浇?,沒有原位探測去驗(yàn)證。

    IBEX衛(wèi)星在地球軌道上測繪了全天空的能量中性原子通量分布,發(fā)現(xiàn)了不曾預(yù)期的一個(gè)近乎封閉的狹窄的條帶狀結(jié)構(gòu)(IBEX-ribbon)[20]。條帶狀的能量中性原子被認(rèn)為是來自日球?qū)禹敻浇鼌^(qū)域,其分布走向被推測是受局地星際磁場所控制。IBEX-ribbon一經(jīng)發(fā)現(xiàn),被之前篤定的半開放的日球?qū)禹斀Y(jié)構(gòu)開始懷疑。在之前認(rèn)為的結(jié)構(gòu)中,日球?qū)禹斘恍途皖愃频厍虻拇艑禹斘恍停栵L(fēng)在尾巴方位與局地星際介質(zhì)混合起來,而不像在鼻尖方位那樣與局地星際介質(zhì)有明顯的界線。

    關(guān)于日球?qū)禹攷缀挝恍偷牟孪?,日球?qū)游锢淼拈_創(chuàng)者之一尤金·帕克(Eugene Parker)曾經(jīng)做過兩頭極端的推測:如果局地星際介質(zhì)的動(dòng)壓遠(yuǎn)大于局地星際介質(zhì)的磁壓,就好比超阿爾芬速的太陽風(fēng)的動(dòng)壓遠(yuǎn)大于太陽風(fēng)的磁壓,這時(shí)日球?qū)禹敃?huì)是一個(gè)半開放的結(jié)構(gòu);反之,如果局地星際介質(zhì)的動(dòng)壓遠(yuǎn)小于其磁壓,則日球?qū)訒?huì)因?yàn)轭愃瓶勾判?yīng)那樣把局地星際磁場排除在外,而形成一個(gè)封閉的準(zhǔn)球狀結(jié)構(gòu)[1]?;贗BEX的能量中性原子的遙感測繪結(jié)果,開始猜測日球?qū)禹斒欠裉幱赑arker猜想的兩頭極端的中間狀態(tài),即尾巴方位既不是嚴(yán)格封閉,但也不是完全開放的。如果認(rèn)為背景太陽風(fēng)和拾起離子之間沒有通過自激發(fā)或背景湍動(dòng)有效地混合起來,而把背景太陽風(fēng)和拾起離子作為兩個(gè)獨(dú)立的成分來建模,輔以適當(dāng)?shù)木值匦请H介質(zhì)條件,會(huì)得到比較新奇的日球?qū)有螤?,形似癟掉的羊角面包狀[5](見圖8)。要解決上述疑問,就要弄清楚日球?qū)拥娜S形狀,需要對日球?qū)游膊亢蜆O區(qū)方位開展原位探測。

    圖8 日球?qū)拥目赡懿煌螤罴捌鋵?yīng)的不同的局地星際介質(zhì)狀態(tài)Fig.8 Possible different shapes of heliosphere and their corresponding different states of local instellar

    拾起離子作為外日球?qū)犹栵L(fēng)熱壓的重要貢獻(xiàn)來源,尚未被系統(tǒng)探測過。因?yàn)槭捌痣x子的存在所導(dǎo)致的熱壓增強(qiáng)、聲速增大是否會(huì)對終止激波的強(qiáng)度有明顯減弱?拾起離子相對于背景太陽風(fēng)的等離子體離子是否更容易在終止激波處受到明顯加速,從而形成異常宇宙線?這些疑問都是“旅行者號(hào)”飛船越過終止激波后留下來的重要問題。拾起離子越過終止激波之后,在日球?qū)忧蕝^(qū)乃至日球?qū)禹敻浇?,與侵入的星際介質(zhì)流發(fā)生電荷交換所產(chǎn)生的能量中性原子對IBEX-ribbon的觀測貢獻(xiàn)多少?

    拾起離子探測器在設(shè)計(jì)視場時(shí)要充分考慮到拾起離子的初始動(dòng)量方向,即星際介質(zhì)來流方向,可以盡可能多地探測被完全拾起前的低能的新生離子。未來要想準(zhǔn)確解釋拾起離子的探測結(jié)果,需要開始對拾起離子的動(dòng)力學(xué)行為(包括局地行為和全局行為)開展相應(yīng)的理論研究,并對以往的拾起離子探測進(jìn)行全過程的數(shù)據(jù)分析。經(jīng)典的拾起離子速度分布模型(V&S模型)假設(shè)拾起離子產(chǎn)生后受湍動(dòng)作用在太陽風(fēng)參考系中快速散射形成各向同性的鞘狀結(jié)構(gòu)?;谠摷僭O(shè)得到的拾起離子速度分布模型能用來擬合相當(dāng)一部分的探測數(shù)據(jù),說明湍動(dòng)的散射效率可能確實(shí)比較高,但也有一部分探測的速度分布偏離模型擬合,出現(xiàn)偏差的物理原因仍待研究。目前關(guān)于外日球?qū)又惺捌痣x子的理論,預(yù)測其數(shù)密度隨徑向距離是r-1衰減?!靶乱曇疤?hào)”(New Horizons)飛船在20 ~ 40 AU附近的探測表明,拾起離子數(shù)密度的徑向衰減速度可能比理論預(yù)測要慢一些(比如r-0.6)。

    利用現(xiàn)有的片段數(shù)據(jù),輔以冪函數(shù)的徑向距離變化進(jìn)行擬合并外推到終止激波上游(假設(shè)為90 AU),得到的推論:拾起離子數(shù)密度/背景太陽風(fēng)數(shù)密度比值從30 AU實(shí)測的0.038上升到90 AU外推預(yù)測的0.14;拾起質(zhì)子溫度/背景太陽風(fēng)質(zhì)子溫度的比值從30 AU實(shí)測的510倍上升到90 AU外推預(yù)測的2 400倍;相應(yīng)的拾起質(zhì)子熱壓/背景太陽風(fēng)質(zhì)子熱壓的比值從30 AU實(shí)測的19倍上升到90 AU外推預(yù)測的350倍;拾起離子的熱壓在20 AU附近開始超過太陽風(fēng)的磁壓(見圖9)。這些推論讓研究者看到拾起離子在外日球?qū)又械娘@著作用和不可替代的角色地位。雖然目前利用“新視野號(hào)”,獲得了外日球?qū)硬糠秩招木嚯x范圍的拾起離子的數(shù)據(jù),但是拾起離子的實(shí)測數(shù)據(jù)在終止激波前后仍然缺失,特別是日球?qū)游膊亢蜆O區(qū)方位沒有深空航天器到達(dá)。目前“新視野號(hào)”上的SWAP載荷(靜電分析儀)的能量上限不到8 keV/q,無法探測完整的He拾起離子的速度分布,在今后的拾起離子載荷研制中,應(yīng)擴(kuò)大能量上限,從而涵蓋He拾起離子的速度分布。

    圖9 “新視野號(hào)”飛船上SWAP載荷測量到的離子計(jì)數(shù)率隨單位電荷能量的變化剖面和根據(jù)“新視野號(hào)”和“旅行者號(hào)”飛船探測,以及模型假設(shè),所畫出的不同壓強(qiáng)隨日心距離的變化以及相互之間大小隨單位電荷能量的變化剖面Fig.9 Profile of ion count rate versus energy per charge measured by SWIP payload onboard New Horizons spacecraft,consisting of solar wind protons,solar wind helium ions,pickup protons and pickup helium ions and changes of different pressures with the heliocentric distance,according to the measurements from New Horizons and Voyage,as well as the model prediction

    在外日球?qū)?,根?jù)拾起過程的物理圖像,拾起離子的速度分布很廣,速度的下限可以低至局地星際介質(zhì)流的流速達(dá)到20多km/s(對應(yīng)能量為幾個(gè)eV),速度的上限可以高至兩倍的太陽風(fēng)流速達(dá)到上千千米每秒(對應(yīng)能量為10 keV量級)。所以未來的探測任務(wù),如果要想完整探測拾起離子的能譜輪廓,需要涵蓋的能量需從幾個(gè)eV到幾十個(gè)keV。拾起離子的探測和背景太陽風(fēng)離子的探測可以整合在一個(gè)儀器上,實(shí)現(xiàn)靜電分析+飛行時(shí)間質(zhì)譜分析的技術(shù)。儀器設(shè)置不同的開口,其中窄的和寬的開口分別接收太陽風(fēng)離子和拾起離子,開口之后是靜電偏置雙層曲面用于選擇不同單位電荷能量的離子,再之后分別是無場飛行路徑和符合探頭部分。

    超熱離子是如何能化(加速)形成的,這個(gè)問題不同的學(xué)者提出了不同的理論機(jī)制,沒有一個(gè)公認(rèn)的答案。日球?qū)又袕V泛存在CIR和MIR所伴隨的激波被認(rèn)為是加速離子的有效場所。因此,有人提出擴(kuò)散激波加速機(jī)制是超熱離子冪律譜的形成原因(見圖10)。但是擴(kuò)散激波加速機(jī)制所形成冪律譜的指數(shù)是可變的,隨激波強(qiáng)度(壓縮比)而變化,這與在壓縮區(qū)里觀測到的冪律譜指數(shù)比較穩(wěn)定的情況不太一致。也有人提出隨機(jī)加速機(jī)制,但是該機(jī)制一般給出的是指數(shù)函數(shù)的能譜形式。多個(gè)磁島之間的磁重聯(lián)也能夠加速離子,但是日球?qū)犹栵L(fēng)中(包括激波下游鞘區(qū))磁重聯(lián)不是很常見,所以磁重聯(lián)加速可能不適用于解釋超熱離子冪律譜的形成機(jī)理。

    圖10 “旅行者號(hào)”飛船在不同年份對超熱離子和能量粒子(宇宙線)的測量能譜和三維磁流體可壓縮湍動(dòng)的模擬結(jié)果顯示出多尺度的壓力平衡結(jié)構(gòu)Fig.10 The energy spectra of super-thermal ions and energetic ions measured by Voyager spacecraft in different years and the simulation results of 3D MHD compressible turbulence showing multi-scale pressure-balanced structures

    從上述幾個(gè)機(jī)制討論中,可以看到湍動(dòng)都是起著非常重要的作用:激波上下游的湍動(dòng)導(dǎo)致離子多次穿越激波發(fā)生擴(kuò)散激波加速;湍動(dòng)場導(dǎo)致場與粒子之間能量交換的隨機(jī)性從而產(chǎn)生隨機(jī)加速;湍動(dòng)串級的間歇性可能有利于相干性結(jié)構(gòu)如磁島的形成,為磁島間磁重聯(lián)創(chuàng)造條件。前面這些機(jī)制沒有考慮能化離子對湍動(dòng)場的反饋?zhàn)饔茫瑢Ω吣芰康碾x子而言,這種忽略是合理的,但是超熱離子由于其對整體的數(shù)密度和熱壓貢獻(xiàn)不能忽略,所以會(huì)反過來影響湍動(dòng)場的可壓縮性,比如原來是壓縮的場會(huì)受到能化的超熱離子增加而變成膨脹的場。一方面,超熱離子中能量較低的離子接著參與局地的膨脹-壓縮的循環(huán),而能量較高的離子則有機(jī)會(huì)在空間中擴(kuò)散開來;另一方面,能量較高的離子又是由能量較低的離子壓縮能化后所形成的。綜合這兩方面的行為,有學(xué)者提出“泵加速”(pump acceleration)機(jī)制,可壓縮湍動(dòng)場的壓縮-膨脹過程類似一個(gè)泵的運(yùn)行,把部分離子能化到更高能量而且有擴(kuò)散自由的狀態(tài)。場和粒子的相互耦合作用(包括粒子對場的反饋?zhàn)饔茫?,這是空間等離子體湍動(dòng)物理研究的重要領(lǐng)域,目前在低能部分已經(jīng)找到充實(shí)的場-粒子相互作用的證據(jù),并得到了湍動(dòng)的耗散率譜。然而要在超熱離子能段,研究場-粒子耦合作用以及其中可能存在的泵加速機(jī)理(包括理論研究、探測設(shè)計(jì)和觀測分析)。

    4 日球?qū)犹綔y

    日球?qū)幼鳛橐粋€(gè)宜居性的星球?qū)樱╤abitable astrosphere),認(rèn)識(shí)日球?qū)拥默F(xiàn)狀乃至其過去與未來,都將為人類在茫茫銀河系和宇宙中探索宜居性星球?qū)犹峁┲匾膮⒖季€索。如同認(rèn)識(shí)地球空間環(huán)境那樣,要對不同方位的磁層進(jìn)行原位探測,磁層頂鼻尖、側(cè)翼邊界層、磁尾等一個(gè)都不能落,這樣才能認(rèn)清地球空間環(huán)境的全貌,摸清地球空間受太陽風(fēng)影響動(dòng)態(tài)演化的機(jī)理。

    日球?qū)禹?,作為從一個(gè)空間(行星際)跨越到另一個(gè)空間(恒星際空間)的分水嶺(見圖11),能夠從不同的方位穿過它,都將是人類探索茫茫宇宙空間的壯舉,也會(huì)帶領(lǐng)領(lǐng)略日球?qū)禹敳煌轿唬ū羌?、尾部、極區(qū))的特色與差異。關(guān)鍵探測要素(特別是拾起離子)的缺位,已經(jīng)成為恒星際空間探路者(“旅行者1號(hào)”和“旅行者2號(hào)”飛船)永遠(yuǎn)的遺憾,但也為后來者繼續(xù)探秘外日球?qū)犹栵L(fēng)與侵入星際風(fēng)的相互作用機(jī)理與效應(yīng)提供了繼續(xù)前行的動(dòng)力。

    圖11 日球?qū)咏Y(jié)構(gòu)示意圖Fig.11 The schematic diagram shows the structure of the heliosphere

    在浩瀚的日球?qū)涌臻g這個(gè)舞臺(tái)上,離子作為主角之一,會(huì)不停地變換角色,比如從拾起離子膨脹冷卻下來與背景太陽風(fēng)離子混合,拾起離子被可壓縮湍動(dòng)能化形成超熱離子,超熱離子在終止激波處進(jìn)一步加速形成異常宇宙線。要追蹤離子在不同角色之間的串位,需要在寬能段范圍里(比如攜帶一套寬能段離子能譜儀:可包括法拉第杯、靜電分析儀、和固體半導(dǎo)體探測器等)測量它們(見表1)。

    表1 針對外日球?qū)犹栵L(fēng)與侵入星際介質(zhì)流的寬能段粒子測量的建議探測載荷和方案Table 1 Possible payloads and schemes suggested for the measurement of wide-energy-band ions in the outer heliosphere and the intrusive interstellar medium flow

    “新視野號(hào)”航天器作為外日球?qū)犹綔y任務(wù)因?yàn)闆]有攜帶磁強(qiáng)計(jì),所以無法真正開展與“場-粒子耦合作用”這一基本物理過程有關(guān)的“超熱離子受可壓縮湍動(dòng)能化加速形成”的前沿課題。因此“探測行星際乃至恒星際中微弱的磁場及其湍動(dòng)”應(yīng)該成為未來太陽系邊際探測任務(wù)設(shè)計(jì)中所面對的眾多探測要素中的必選項(xiàng)。專門針對太陽系邊際開展抵達(dá)探測的設(shè)想,是驅(qū)動(dòng)深空探測能力邁上新臺(tái)階的重要源泉,也被認(rèn)為將是我國建設(shè)航天強(qiáng)國的重要標(biāo)志[39]。

    5 結(jié)束語

    本文圍繞外日球?qū)犹栵L(fēng)及其與侵入星際風(fēng)的作用從目前認(rèn)知、前沿問題、探測建議等3個(gè)方面展開論述,目前認(rèn)知是基于當(dāng)前的探測能力和探測結(jié)果,關(guān)于太陽風(fēng)和星際風(fēng)相互作用的理想探測模式要求能夠同時(shí)探測相互作用前后的所有種類粒子。這些粒子的能量分布范圍很寬,形成寬能段的粒子能譜。對離子而言,在3個(gè)不同的能段包括三種不同的粒子:原初太陽風(fēng)離子的麥?zhǔn)夏茏V、星際風(fēng)拾起離子的能譜平臺(tái)和能譜截?cái)?、超熱離子的冪律能譜。拾起離子和超熱離子的能譜有來自局地星際介質(zhì)侵入流的貢獻(xiàn)。

    給出了面臨的前沿問題有3個(gè)方面:①從未到達(dá)的領(lǐng)域,即黃道面日球?qū)游膊康耐馊涨驅(qū)犹栵L(fēng)以及高緯外日球?qū)犹栵L(fēng)從沒有被探測過,其與日球?qū)颖羌夥轿坏奶栵L(fēng)差異如何不得而知;②曾經(jīng)到達(dá)的領(lǐng)域,但是沒有探測的物理量,例如外日球?qū)犹栵L(fēng)中的拾起離子;③曾經(jīng)到達(dá)并探測,但是形成機(jī)制不明,超熱離子的冪律譜及其雙模態(tài)。

    要回答前沿問題,通過探討嘗試給出建議:①設(shè)計(jì)不同的飛行路徑,朝向不同的方位進(jìn)行探測;②攜帶覆蓋寬能段的離子譜儀器包涵蓋原初太陽風(fēng)等離子體、拾起離子和超熱離子3個(gè)能段;③搭載高靈敏度磁強(qiáng)計(jì)測量弱磁場的可壓縮湍動(dòng)。

    致 謝

    本文撰寫受到北京大學(xué)、中國科學(xué)院國家空間科學(xué)中心、探月與航天工程中心、北京空間飛行器總體設(shè)計(jì)部等多家單位的老師們的幫助,在此對涂傳詒、王赤、李暉、郭孝城、吳偉仁、王倩、龐涪川、康炎、黃江川、孟林智等專家表示感謝!

    猜你喜歡
    太陽風(fēng)激波星際
    通天達(dá)地:穿越星際的通信
    一種基于聚類分析的二維激波模式識(shí)別算法
    多種觀測數(shù)據(jù)驅(qū)動(dòng)的三維行星際太陽風(fēng)MHD模擬
    基于HIFiRE-2超燃發(fā)動(dòng)機(jī)內(nèi)流道的激波邊界層干擾分析
    基于ACE飛船觀測的銀河宇宙線與太陽風(fēng)變化的統(tǒng)計(jì)研究
    星際打劫案
    星際奇葩說
    斜激波入射V形鈍前緣溢流口激波干擾研究
    “穿越星際”去上課
    適于可壓縮多尺度流動(dòng)的緊致型激波捕捉格式
    欧美日韩瑟瑟在线播放| 亚洲国产精品久久男人天堂| 国产成人a区在线观看| 欧美日韩国产亚洲二区| av在线蜜桃| av国产免费在线观看| 欧美一区二区国产精品久久精品| 亚洲欧美清纯卡通| 日韩欧美国产在线观看| 国语自产精品视频在线第100页| 亚洲自拍偷在线| 天堂动漫精品| 久久久久九九精品影院| 久久久久精品国产欧美久久久| 国产精品乱码一区二三区的特点| 精品久久久久久,| 国产精品日韩av在线免费观看| 日韩欧美精品免费久久| 老女人水多毛片| 成人国产综合亚洲| 国产三级在线视频| 在线观看美女被高潮喷水网站| 日韩 亚洲 欧美在线| 久久久久久久亚洲中文字幕| 在线天堂最新版资源| 精品一区二区免费观看| 少妇猛男粗大的猛烈进出视频 | 成人国产一区最新在线观看| 欧美性猛交黑人性爽| 久久国内精品自在自线图片| 又黄又爽又刺激的免费视频.| 变态另类成人亚洲欧美熟女| 久久久久久久久大av| 性插视频无遮挡在线免费观看| 欧美成人免费av一区二区三区| 国产精品亚洲美女久久久| 日韩一本色道免费dvd| 亚洲18禁久久av| 亚洲va在线va天堂va国产| 欧美日韩亚洲国产一区二区在线观看| 一个人看的www免费观看视频| 五月伊人婷婷丁香| 国产精品久久久久久av不卡| 亚洲图色成人| 免费无遮挡裸体视频| 色噜噜av男人的天堂激情| 国产伦人伦偷精品视频| 久久草成人影院| 伦精品一区二区三区| 免费看a级黄色片| 日韩强制内射视频| 欧美3d第一页| 欧美激情久久久久久爽电影| 一a级毛片在线观看| 精品免费久久久久久久清纯| 精品99又大又爽又粗少妇毛片 | 午夜福利视频1000在线观看| 亚洲三级黄色毛片| 久久热精品热| 欧美性感艳星| 成熟少妇高潮喷水视频| 日本成人三级电影网站| 久久中文看片网| 亚洲精品久久国产高清桃花| 天堂网av新在线| 亚洲第一区二区三区不卡| 美女黄网站色视频| 国产欧美日韩一区二区精品| av专区在线播放| 一区二区三区高清视频在线| 嫩草影院精品99| 亚洲欧美激情综合另类| 亚洲va日本ⅴa欧美va伊人久久| 免费看a级黄色片| 亚洲av一区综合| 性欧美人与动物交配| 欧美日韩国产亚洲二区| 97人妻精品一区二区三区麻豆| 国产伦人伦偷精品视频| 精品午夜福利视频在线观看一区| 麻豆成人av在线观看| 丰满的人妻完整版| 最后的刺客免费高清国语| 中文在线观看免费www的网站| av天堂中文字幕网| 波多野结衣高清作品| 亚洲av免费高清在线观看| 免费人成在线观看视频色| 91久久精品电影网| 我要看日韩黄色一级片| 亚洲av电影不卡..在线观看| 夜夜看夜夜爽夜夜摸| 国产亚洲91精品色在线| 在线免费十八禁| 成人一区二区视频在线观看| 神马国产精品三级电影在线观看| 91久久精品国产一区二区三区| 国产精品美女特级片免费视频播放器| av在线亚洲专区| 国产一区二区三区av在线 | 久久久久久久午夜电影| 国产中年淑女户外野战色| 亚洲午夜理论影院| 精品国产三级普通话版| 亚洲黑人精品在线| 啦啦啦观看免费观看视频高清| 中文字幕av成人在线电影| or卡值多少钱| 亚洲avbb在线观看| 国产高清视频在线播放一区| 欧美性猛交╳xxx乱大交人| 色视频www国产| 欧美激情国产日韩精品一区| 亚洲精品一区av在线观看| 嫩草影院精品99| 在线观看一区二区三区| 亚洲四区av| 国产精品电影一区二区三区| 国产三级中文精品| 少妇猛男粗大的猛烈进出视频 | 人妻少妇偷人精品九色| 一级av片app| 亚洲中文日韩欧美视频| 看免费成人av毛片| 成人国产一区最新在线观看| 性插视频无遮挡在线免费观看| 久久国产精品人妻蜜桃| 精品人妻1区二区| 国产黄色小视频在线观看| 午夜福利成人在线免费观看| 亚洲色图av天堂| 国产精品av视频在线免费观看| 天天躁日日操中文字幕| 久久久久久久久中文| 97人妻精品一区二区三区麻豆| 中亚洲国语对白在线视频| 亚洲欧美日韩高清在线视频| 国产一区二区三区av在线 | 亚洲av第一区精品v没综合| 蜜桃久久精品国产亚洲av| 久久香蕉精品热| 日本精品一区二区三区蜜桃| 欧美丝袜亚洲另类 | 俺也久久电影网| 欧美潮喷喷水| 在线观看一区二区三区| 国产亚洲av嫩草精品影院| 日韩一区二区视频免费看| 亚洲精品一区av在线观看| 亚洲国产精品成人综合色| 国产综合懂色| 午夜福利高清视频| 69人妻影院| 又爽又黄a免费视频| 精品一区二区三区人妻视频| or卡值多少钱| 亚洲精品乱码久久久v下载方式| 亚洲人成网站在线播放欧美日韩| 欧美性猛交黑人性爽| 女的被弄到高潮叫床怎么办 | 午夜激情欧美在线| 欧美不卡视频在线免费观看| 国产一区二区在线av高清观看| aaaaa片日本免费| 国产精品一区二区三区四区免费观看 | 久久欧美精品欧美久久欧美| 中国美女看黄片| 免费观看精品视频网站| 亚洲美女视频黄频| 一级黄片播放器| 成人高潮视频无遮挡免费网站| 99热网站在线观看| 成年免费大片在线观看| 免费无遮挡裸体视频| 五月玫瑰六月丁香| 成熟少妇高潮喷水视频| 神马国产精品三级电影在线观看| 美女cb高潮喷水在线观看| netflix在线观看网站| 美女免费视频网站| 日韩亚洲欧美综合| 国产精品亚洲美女久久久| 久久久久久久久大av| 午夜激情福利司机影院| 动漫黄色视频在线观看| 综合色av麻豆| 九色国产91popny在线| 日本五十路高清| 男女视频在线观看网站免费| 日韩大尺度精品在线看网址| 国产真实伦视频高清在线观看 | 97碰自拍视频| 免费搜索国产男女视频| 日韩在线高清观看一区二区三区 | 日本黄色片子视频| 日本黄大片高清| 欧美一区二区亚洲| 亚洲在线自拍视频| 精品久久久久久,| 嫁个100分男人电影在线观看| 成人综合一区亚洲| 成人国产麻豆网| 亚洲国产色片| 国产精品国产高清国产av| 国产高清视频在线观看网站| 中文字幕免费在线视频6| 999久久久精品免费观看国产| 我要搜黄色片| 亚洲欧美日韩高清在线视频| 日韩欧美一区二区三区在线观看| 亚洲人与动物交配视频| 精品一区二区三区视频在线观看免费| 欧美一区二区亚洲| 国产成人a区在线观看| 亚洲美女黄片视频| 久久久久九九精品影院| 欧美不卡视频在线免费观看| 偷拍熟女少妇极品色| 亚洲人成网站高清观看| 五月伊人婷婷丁香| 99热这里只有是精品在线观看| 欧美区成人在线视频| 九九热线精品视视频播放| 九九爱精品视频在线观看| 又粗又爽又猛毛片免费看| 免费人成视频x8x8入口观看| 国产亚洲精品av在线| 一进一出好大好爽视频| 欧美3d第一页| 97碰自拍视频| 亚洲国产色片| 免费观看在线日韩| 天天躁日日操中文字幕| 日韩,欧美,国产一区二区三区 | 日本熟妇午夜| a级毛片a级免费在线| 日韩欧美在线二视频| 日韩欧美 国产精品| 狠狠狠狠99中文字幕| 国产色婷婷99| 一进一出好大好爽视频| 成熟少妇高潮喷水视频| 无遮挡黄片免费观看| 亚洲最大成人手机在线| 波野结衣二区三区在线| 岛国在线免费视频观看| 国产欧美日韩精品一区二区| 久久中文看片网| 极品教师在线视频| 精品人妻1区二区| 一本一本综合久久| 在线观看66精品国产| 中文字幕熟女人妻在线| 看片在线看免费视频| 精品一区二区免费观看| 成人二区视频| 波多野结衣高清无吗| www.www免费av| 搡老熟女国产l中国老女人| 亚洲18禁久久av| 欧美高清性xxxxhd video| av中文乱码字幕在线| 久久久精品大字幕| 亚洲美女视频黄频| 一卡2卡三卡四卡精品乱码亚洲| 在线观看av片永久免费下载| 国产淫片久久久久久久久| 欧美日韩瑟瑟在线播放| 国产精品精品国产色婷婷| av在线老鸭窝| 午夜福利欧美成人| 免费看日本二区| 老熟妇仑乱视频hdxx| 性色avwww在线观看| 久久久久久久久久成人| 国产精品精品国产色婷婷| 成人综合一区亚洲| 91麻豆精品激情在线观看国产| 国产单亲对白刺激| 日本五十路高清| 国产淫片久久久久久久久| 日本成人三级电影网站| 亚洲五月天丁香| 欧美绝顶高潮抽搐喷水| 在线观看美女被高潮喷水网站| 联通29元200g的流量卡| 国产精品一区二区三区四区久久| av黄色大香蕉| 精品人妻1区二区| 午夜精品久久久久久毛片777| 免费av不卡在线播放| 俄罗斯特黄特色一大片| 中文字幕av在线有码专区| 国产成年人精品一区二区| 亚洲avbb在线观看| 一a级毛片在线观看| 狂野欧美白嫩少妇大欣赏| 欧美日韩黄片免| 夜夜爽天天搞| 国产毛片a区久久久久| 久久久精品欧美日韩精品| 日日夜夜操网爽| 日日撸夜夜添| 亚洲av电影不卡..在线观看| 人人妻人人澡欧美一区二区| 久9热在线精品视频| 老熟妇乱子伦视频在线观看| а√天堂www在线а√下载| 人人妻人人看人人澡| 国产精品一区二区三区四区久久| 又爽又黄a免费视频| 在线观看一区二区三区| 精品久久久久久成人av| 亚洲精品日韩av片在线观看| 韩国av在线不卡| 久久热精品热| 亚洲人与动物交配视频| 国产三级中文精品| 一级黄片播放器| 午夜福利高清视频| 国产麻豆成人av免费视频| 国产亚洲精品久久久久久毛片| 夜夜爽天天搞| 免费av不卡在线播放| 亚洲aⅴ乱码一区二区在线播放| 在线国产一区二区在线| 亚洲性夜色夜夜综合| 麻豆一二三区av精品| 深夜精品福利| 十八禁国产超污无遮挡网站| 精品一区二区免费观看| 亚洲avbb在线观看| 国产伦精品一区二区三区视频9| 成年女人毛片免费观看观看9| av中文乱码字幕在线| 国产中年淑女户外野战色| 12—13女人毛片做爰片一| 日本欧美国产在线视频| 99久久久亚洲精品蜜臀av| 美女cb高潮喷水在线观看| 国产亚洲精品久久久久久毛片| 日本 欧美在线| 老熟妇乱子伦视频在线观看| 日本 欧美在线| 成年版毛片免费区| 久久亚洲真实| 国语自产精品视频在线第100页| 亚洲欧美日韩卡通动漫| 亚洲色图av天堂| 欧美在线一区亚洲| 成年版毛片免费区| videossex国产| 看十八女毛片水多多多| 狠狠狠狠99中文字幕| 免费看光身美女| 免费看美女性在线毛片视频| 亚洲欧美日韩东京热| 免费大片18禁| 看免费成人av毛片| 日韩av在线大香蕉| h日本视频在线播放| 狠狠狠狠99中文字幕| 亚洲欧美日韩东京热| 少妇裸体淫交视频免费看高清| 国产亚洲精品av在线| 九九热线精品视视频播放| 国产探花在线观看一区二区| 欧美一区二区国产精品久久精品| 在线观看舔阴道视频| a级毛片免费高清观看在线播放| 可以在线观看毛片的网站| 亚洲精品粉嫩美女一区| 国产在线男女| 久久久久久大精品| 亚洲成av人片在线播放无| 国产精品三级大全| 国产欧美日韩精品一区二区| 琪琪午夜伦伦电影理论片6080| 精品久久久久久久人妻蜜臀av| 国产精品日韩av在线免费观看| 国产av麻豆久久久久久久| 哪里可以看免费的av片| 色综合婷婷激情| a级毛片免费高清观看在线播放| 精品人妻视频免费看| 黄色视频,在线免费观看| 亚洲三级黄色毛片| 国产伦精品一区二区三区视频9| 一卡2卡三卡四卡精品乱码亚洲| www日本黄色视频网| 又爽又黄无遮挡网站| 神马国产精品三级电影在线观看| 桃色一区二区三区在线观看| 日本与韩国留学比较| 一级黄片播放器| 男女啪啪激烈高潮av片| 人妻少妇偷人精品九色| 2021天堂中文幕一二区在线观| 亚洲一区高清亚洲精品| 日韩欧美精品v在线| 性欧美人与动物交配| 国产三级在线视频| 婷婷丁香在线五月| 久久久久久久久大av| 日韩欧美国产在线观看| 国产高清视频在线观看网站| 中文字幕av在线有码专区| 淫秽高清视频在线观看| 欧美+日韩+精品| 男女之事视频高清在线观看| 伦理电影大哥的女人| 少妇猛男粗大的猛烈进出视频 | 赤兔流量卡办理| 久久久久久大精品| 国产亚洲精品综合一区在线观看| 国产精品自产拍在线观看55亚洲| 欧美日韩国产亚洲二区| 欧美激情在线99| 白带黄色成豆腐渣| 99riav亚洲国产免费| 欧美三级亚洲精品| 无人区码免费观看不卡| 国产午夜精品久久久久久一区二区三区 | 久久九九热精品免费| 夜夜看夜夜爽夜夜摸| 午夜精品在线福利| 12—13女人毛片做爰片一| 午夜精品久久久久久毛片777| 12—13女人毛片做爰片一| 美女 人体艺术 gogo| 精品久久久久久,| 亚洲av不卡在线观看| www日本黄色视频网| 精品久久久久久久久亚洲 | 久久久久久九九精品二区国产| 好男人在线观看高清免费视频| 午夜福利在线观看免费完整高清在 | 国产精品免费一区二区三区在线| av专区在线播放| 露出奶头的视频| 国产伦一二天堂av在线观看| 啪啪无遮挡十八禁网站| 不卡一级毛片| xxxwww97欧美| 18禁黄网站禁片免费观看直播| 99热这里只有精品一区| 国产私拍福利视频在线观看| 99热精品在线国产| or卡值多少钱| 成年版毛片免费区| 国产欧美日韩一区二区精品| 久久人人精品亚洲av| 欧美成人免费av一区二区三区| 久久精品夜夜夜夜夜久久蜜豆| 国产在视频线在精品| 国产91精品成人一区二区三区| 99精品在免费线老司机午夜| 99九九线精品视频在线观看视频| 国产精品98久久久久久宅男小说| 午夜精品久久久久久毛片777| 亚洲精品一区av在线观看| 亚洲avbb在线观看| 97超视频在线观看视频| 欧美高清成人免费视频www| 两个人视频免费观看高清| 亚洲人成网站在线播放欧美日韩| 亚洲人与动物交配视频| 午夜福利欧美成人| 国产色婷婷99| 午夜亚洲福利在线播放| 亚洲成人免费电影在线观看| 国产美女午夜福利| 日本色播在线视频| 亚洲人成网站高清观看| 欧美一区二区亚洲| 九九在线视频观看精品| 一进一出抽搐动态| 亚洲精品成人久久久久久| 国产 一区 欧美 日韩| 欧美一区二区国产精品久久精品| 又黄又爽又免费观看的视频| 久久午夜亚洲精品久久| 村上凉子中文字幕在线| 高清在线国产一区| 亚洲精品粉嫩美女一区| 亚洲国产精品合色在线| 久久九九热精品免费| 亚洲精品一区av在线观看| 亚洲精品乱码久久久v下载方式| 免费高清视频大片| 伊人久久精品亚洲午夜| a在线观看视频网站| 国产爱豆传媒在线观看| 欧美极品一区二区三区四区| 88av欧美| 国产免费男女视频| 亚洲最大成人手机在线| 日韩欧美精品v在线| 老师上课跳d突然被开到最大视频| 99在线人妻在线中文字幕| 欧美日韩国产亚洲二区| 亚洲欧美清纯卡通| 神马国产精品三级电影在线观看| 亚洲图色成人| 国产一区二区激情短视频| 黄色配什么色好看| 日本欧美国产在线视频| 亚洲欧美日韩无卡精品| 精品一区二区三区av网在线观看| 成人特级av手机在线观看| 亚洲国产精品合色在线| 在线免费十八禁| 精品久久久久久久久亚洲 | 少妇裸体淫交视频免费看高清| 欧美区成人在线视频| 一级黄片播放器| 国产蜜桃级精品一区二区三区| a级一级毛片免费在线观看| 精品免费久久久久久久清纯| 日本黄大片高清| 成年女人毛片免费观看观看9| 亚洲av美国av| av在线亚洲专区| 俺也久久电影网| 国产男靠女视频免费网站| 亚洲国产精品久久男人天堂| 成人特级av手机在线观看| 一个人看的www免费观看视频| 可以在线观看的亚洲视频| avwww免费| 国产探花在线观看一区二区| 欧美zozozo另类| 老熟妇仑乱视频hdxx| 校园春色视频在线观看| 男女边吃奶边做爰视频| 波野结衣二区三区在线| 亚洲欧美激情综合另类| 久久中文看片网| 亚洲性久久影院| 午夜老司机福利剧场| 97碰自拍视频| 亚洲天堂国产精品一区在线| 欧美一区二区亚洲| 午夜福利在线观看免费完整高清在 | 欧美日韩乱码在线| 日本黄大片高清| 女的被弄到高潮叫床怎么办 | 国产精品乱码一区二三区的特点| 1024手机看黄色片| 国内揄拍国产精品人妻在线| 亚洲中文日韩欧美视频| 日韩高清综合在线| 99国产极品粉嫩在线观看| 老师上课跳d突然被开到最大视频| 美女xxoo啪啪120秒动态图| 久久精品91蜜桃| 91精品国产九色| 国产精品一区二区三区四区久久| 久久久久国产精品人妻aⅴ院| 乱系列少妇在线播放| 久久婷婷人人爽人人干人人爱| 国产亚洲精品综合一区在线观看| 高清毛片免费观看视频网站| 亚洲乱码一区二区免费版| 最近最新中文字幕大全电影3| 午夜亚洲福利在线播放| 一区二区三区高清视频在线| 九九爱精品视频在线观看| 狂野欧美激情性xxxx在线观看| 九九在线视频观看精品| 久久精品国产99精品国产亚洲性色| 亚州av有码| 又黄又爽又免费观看的视频| 国产成人av教育| www.色视频.com| 国产伦在线观看视频一区| 真人一进一出gif抽搐免费| 精华霜和精华液先用哪个| 国产激情偷乱视频一区二区| 久久久久国内视频| 亚洲国产精品sss在线观看| 日韩大尺度精品在线看网址| 夜夜爽天天搞| 2021天堂中文幕一二区在线观| 亚洲在线观看片| 岛国在线免费视频观看| 国产三级中文精品| 亚洲人成网站高清观看| 99国产精品一区二区蜜桃av| 久久欧美精品欧美久久欧美| 乱码一卡2卡4卡精品| 日本a在线网址| 老熟妇仑乱视频hdxx| 久久久久精品国产欧美久久久| 制服丝袜大香蕉在线| 22中文网久久字幕| 久久久久精品国产欧美久久久| 午夜亚洲福利在线播放| 亚洲专区中文字幕在线| 久久久久精品国产欧美久久久| 日本免费a在线| 中文字幕av成人在线电影| 久久人妻av系列| 日日干狠狠操夜夜爽| 精品人妻偷拍中文字幕| 九色国产91popny在线| 欧美在线一区亚洲| av在线观看视频网站免费| 中文字幕久久专区| 成人无遮挡网站| 国产真实伦视频高清在线观看 | 91久久精品国产一区二区成人| 午夜免费成人在线视频| 精品一区二区三区av网在线观看| 久久亚洲真实|