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

    鋇和銫釋放的電離層擾動效應(yīng)對比

    2020-02-18 03:18:16朱肖麗胡耀垓趙正予張援農(nóng)
    物理學(xué)報 2020年2期
    關(guān)鍵詞:云團(tuán)電離層風(fēng)場

    朱肖麗 胡耀垓 趙正予 張援農(nóng)

    (武漢大學(xué)電子信息學(xué)院,電離層實驗室,武漢 430079)

    堿金屬或堿土金屬在電離層釋放后,迅速在太陽輻射作用下發(fā)生光電離,產(chǎn)生正離子和電子,形成人工等離子體云團(tuán).本文基于三維雙成分流體模型,考慮釋放區(qū)域水平風(fēng)場的影響,探討了鋇和銫在電離層釋放后的時空演化規(guī)律,并對鋇和銫的電離層擾動效應(yīng)進(jìn)行了對比.模擬結(jié)果表明,不考慮中性風(fēng)場時,生成的等離子體云團(tuán)逐漸沿磁場被拉伸成橢球形結(jié)構(gòu),同時,膨脹的等離子體云會推開背景氧離子,在釋放中心形成氧離子密度空洞,并在兩側(cè)產(chǎn)生兩個對稱的密度尖峰;水平風(fēng)場的存在會使得生成的離子云逆風(fēng)側(cè)的密度梯度變陡,釋放物質(zhì)對背景氧離子的擾動也更大;對比鋇與銫的釋放結(jié)果發(fā)現(xiàn),由于銫的擴(kuò)散系數(shù)較小,鋇云的膨脹更為迅速,Ba+云團(tuán)的覆蓋區(qū)域更廣;而由于光電離率較大,釋放相同質(zhì)量下銫的離子產(chǎn)率更高;此外,Cs+的掃雪機(jī)效應(yīng)比Ba+掃雪機(jī)更強(qiáng),氧離子密度空穴和凸起處的擾動也更大.

    1 引 言

    20世紀(jì)60年代,人類首次發(fā)現(xiàn)火箭發(fā)射期間對電離層電子密度的擾動現(xiàn)象.此后,美歐等國開展了大量化學(xué)物質(zhì)主動釋放擾動電離層的試驗[1?4],包括電離層電子密度耗空類物質(zhì)和電子密度增強(qiáng)類物質(zhì),其中,電離層電子密度增強(qiáng)試驗通常的釋放物質(zhì)為堿金屬或堿土金屬,如 Ba,Cs,Li,Na,Ca等.

    通過空間物理主動釋放試驗、實驗室研究以及數(shù)值模擬研究等對人工等離子體云團(tuán)釋放后的演化特性進(jìn)行研究已經(jīng)持續(xù)了很多年[5?8].早在1967年,Haerendel等[9]就通過簡化的低密度擾動模型對人工等離子體云進(jìn)行了初步定性探討.1988 年,Schunk和Szuszczewicz[10]首先利用一維Vlasov-Poisson模型研究了高密度 Ba+,Li+和Ba+-Li+等離子體云向低密度O+背景等離子體的無碰撞擴(kuò)展,以闡明早期等離子體膨脹的離子特征.隨后,Mitchell等[11]利用二維靜電模型對等離子體云沿磁場方向和垂直磁場的運動進(jìn)行了理論和數(shù)值研究;為了全面地描述等離子體的膨脹和三維運動,更為精細(xì)的三維模型被構(gòu)建,以研究背景中性風(fēng)場、電磁場作用、粒子間碰撞效應(yīng)和慣性作用影響下等離子體云團(tuán)的膨脹特性[12?15].

    化學(xué)釋放人工擾動電離層的另一類物質(zhì)為中性氣體,從20世紀(jì)60年代開始,國外開展了大量中性物質(zhì)釋放形成電離層電子密度“空洞”的空間主動試驗[16?18].基于試驗觀測結(jié)果,很多學(xué)者也從理論和數(shù)值模擬的角度對中性氣體釋放擾動電離層的機(jī)理及一些觀測到的效應(yīng)進(jìn)行了理論解釋[19?22].近年來,國內(nèi)也開展了中性物質(zhì)釋放人工干擾電離層的相關(guān)研究.黃文耿和古士芬[23]、黃勇等[24]、胡耀垓等[25,26]和汪四成等[27]基于二維動力學(xué)模型,綜合考慮了中性氣體熱擴(kuò)散、等離子體雙極擴(kuò)散、以及離子化學(xué)反應(yīng)等過程,數(shù)值模擬了多種中性氣體釋放后對電離層的擾動特性,并比較了不同釋放條件下的電離層擾動結(jié)果;趙海生等[28]建立了更為精細(xì)的模型,考慮了熱層風(fēng)場、釋放物初速度及運載器飛行姿態(tài)等參量對中性擴(kuò)散過程的影響,并引入地磁傾角和沿場擴(kuò)散項,將二維動力學(xué)模型擴(kuò)展到了三維.

    相比于中性物質(zhì)的釋放,國內(nèi)對堿土金屬類物質(zhì)釋放的研究較少.Li和Xu[29]研究了在擴(kuò)散、碰撞、氧化和光電離作用下不同高度中性鋇云的自擴(kuò)散過程.胡耀垓等[30]在中性鋇云自擴(kuò)散模型的基礎(chǔ)上,考慮鋇原子的氧化和光電離損耗,探討了釋放早期中性鋇云形態(tài)、亮度分布以及釋放區(qū)域電子密度分布.謝良海[31]利用美國BATS-R-US代碼,基于多成分磁流體模型,模擬了鋇的釋放效應(yīng)及其與環(huán)境等離子體的相互作用過程.

    開展空間化學(xué)物質(zhì)主動釋放試驗,觀測釋放物質(zhì)與背景電離層的相互作用過程,對于相關(guān)空間物理問題(如磁重聯(lián)、等離子體不穩(wěn)定性、等離子體波等)的研究,高層大氣風(fēng)場和電磁場的測量,實現(xiàn)人工電離層變態(tài),影響短波通信及衛(wèi)星通信等都具有非常重要的意義[1,4,7,9,29,32].金屬鋇由于質(zhì)量輕、電離電位低、易氣化以及便于觀測等優(yōu)點,成為空間物理主動試驗中最常用的釋放物.但是,鋇在電離層的電離依賴于太陽光的作用,在黑暗條件下無法進(jìn)行光電離,而銫作為堿金屬族中電離勢最低的元素,除了光電離外,銫在無光環(huán)境下也能通過自身熱電離產(chǎn)生電子,打破了光照條件對釋放試驗的限制,早期的空間主動試驗也常用銫作為釋放物[33?35].相比鋇而言,銫釋放的模擬和試驗研究都尚不充分,也鮮見報道.開展鋇和銫釋放效應(yīng)的對比研究,比較不同人工等離子體云團(tuán)的演化規(guī)律及其對背景電離層的擾動特性,對于空間主動試驗中釋放物的選擇具有一定的指導(dǎo)意義.

    本文基于堿金屬在電離層釋放的三維雙成分流體模型,從粒子的連續(xù)性方程、動量方程以及電流平衡方程出發(fā),考慮了背景電磁場力、各種粒子成分間的碰撞及背景水平風(fēng)場作用,模擬了鋇和銫釋放后等離子體云團(tuán)和背景主要粒子的時空演化規(guī)律.

    2 理論模型

    由于釋放物質(zhì)光電離產(chǎn)生的離子和背景電離層離子成分之間具有不同的空間分布,我們將不同離子分開處理.此外,在電離層F層,背景離子成分主要為氧離子,因此這里考慮金屬離子和氧離子的雙成分模型.忽略由于化學(xué)反應(yīng)引起的背景氧離子的密度變化,并假設(shè)離子云的產(chǎn)生僅來源于光電離,釋放區(qū)域等離子體滿足的一般方程如下:

    其中下標(biāo)1,2,e和s分別表示氧離子、金屬離子、電子和釋放的中性物質(zhì).n和u分別表示各粒子的數(shù)密度和漂移速度.(1)式和(2)式是離子的連續(xù)性方程,σ 是離子云的光電離產(chǎn)生率.假設(shè)光電離產(chǎn)生率為離子的凈生成速率,即忽略電荷交換和復(fù)合效應(yīng).在電離層F區(qū),準(zhǔn)中性條件是滿足的(方程 (3)).離子的動量方程由 (4)式給出,其中,nj,mj,uj,ej和pj分別表示離子 j的數(shù)密度、質(zhì)量、漂移速度、電荷量和壓強(qiáng);Dj/Dt=?/?t+uj·?表示遷移導(dǎo)數(shù);E和 B是電場和磁場;δ Mj/δt 代表動量變化;? pj=?kTjnj為壓強(qiáng)梯度項;J 代表總電流密度.

    氧離子和金屬離子的動量變化分別為:

    νjα是粒子間的碰撞頻率(α,s表示帶電粒子和中性粒子),us是中性原子的速度.除了由于碰撞帶來的動量變化,離子2還應(yīng)包括源項引起的動量變化.由于地磁場的束縛作用,在垂直于磁場的方向上帶電粒子之間的速度差很小,與之有關(guān)的碰撞項可以忽略[14],同時忽略電子慣性,將離子慣性項作為小的修正項,并假設(shè)地磁場恒定,求解動量方程可以得到各粒子漂移速度的表達(dá)式[14]:

    其中 z=B/|B| 表示沿磁場方向的單位向量.Ωj=ejB/mj是回旋頻率.aj和sj的表達(dá)式如下:

    電場由公式 E=E0-?φ 計算,其中 E0是背景恒定電場,φ是擾動靜電勢,定義有效電勢Φ=φ-(kTe/e)ln(ne),根據(jù)電荷守恒,可以得到關(guān)于有效電勢的方程如下:

    其中σ⊥為 Pedersen 電導(dǎo) 率,σ//為平行電導(dǎo)率,Q為總的電荷量[14].

    釋放中性云團(tuán)粒子的密度和速度由以下方程給出[27]:

    其中 ns是中性云的數(shù)密度;us是中性云團(tuán)速度;σ表示光電離率;kc是氧化反應(yīng)速率(通常在電離層、磁層高度處 σ ?kcnO2,可不考慮氧化損耗);nO2是氧氣數(shù)密度;D=kT/(mνD) 是擴(kuò)散系數(shù),νD為與擴(kuò)散相關(guān)的碰撞頻率;(x0,y0,z0) 為云團(tuán)初始釋放中心;ux,uy和uz分別是云團(tuán)中心在三個方向上的速度分量;r0是云團(tuán)初始半徑;N0是釋放粒子總數(shù);ε 表示形狀因子.擴(kuò)散系數(shù)D可由以下公式計算[36]:

    其中,nO為氧原子數(shù)密度,為約化質(zhì)量,ms,mO,rs和rO分別為金屬原子和氧原子的質(zhì)量和原子半徑.在電離層F區(qū),中性云會在背景粒子的碰撞作用影響下做減速運動,其速度隨時間的變化可近似表示為

    這里 u0為云團(tuán)初始釋放速度,un為中性風(fēng)場速度,νs為云團(tuán)減速運動的阻尼系數(shù),可根據(jù)彈性碰撞理論求得[37]:

    考慮背景風(fēng)場作用,對于有初速的釋放情形,疊加的移動量如下:

    對上述(1)—(17)式進(jìn)行數(shù)值求解,即可得到任意時刻釋放區(qū)域粒子的空間分布.

    3 數(shù)值模擬結(jié)果及討論

    3.1 仿真算法設(shè)計

    利用有限差分方法求解連續(xù)性方程(1)和(2),方程(12)是橢圓方程,利用交替方向隱式算法求解.時間步長取 ? t=0.01 s,空間步長?x=?y=?z=1 km,基于MATLAB軟件實現(xiàn)了人工等離子體云演化的模擬程序,數(shù)值仿真算法流程如圖1所示.

    圖1 仿真算法流程圖Fig.1.Flow chart of simulation algorithm.

    背景大氣密度、電離層粒子數(shù)密度、溫度、磁場強(qiáng)度等初始條件可由MSIS-E-90大氣模型、電離層IRI-2016模型和磁場模型DGRF/IGRF 1945-2020獲得.所有邊界條件均采用等值外推形式.主要仿真參數(shù)如表1所列.

    表1 主要仿真參數(shù)表Table 1.The main simulation parameters.

    3.2 數(shù)值模擬結(jié)果及討論

    為探討金屬原子在電離層釋放后產(chǎn)生的擾動特性和釋放區(qū)域主要粒子的時空演化規(guī)律,分別模擬了10 kg鋇和銫在300 km高度處的釋放,釋放云團(tuán)的初始特征半徑均取 2 km,計算域為 X=Y=Z=[–25 25]km.所有的剖面圖都是經(jīng)過釋放中心的平面,圖中的等值圖代表粒子數(shù)密度分布,帶電粒子的速度用矢量場表示(圖中黑色箭頭).

    3.2.1 無中性風(fēng)場時

    圖2和圖3描述了 300 km高度釋放 10 kg Ba 原子,釋放后 5,30和200 s鋇離子及背景氧離子的數(shù)密度分布.不考慮背景風(fēng)速時,鋇離子云的早期密度分布為球?qū)ΨQ的,由于磁場的存在,鋇離子云在垂直磁場方向上的膨脹受到束縛(圖2),而在沿磁場方向,由于鋇離子云在密度梯度作用下的運動不受限制,離子云團(tuán)逐漸沿著磁場方向被拉伸,逐漸變成橢球狀結(jié)構(gòu) (圖3).由于碰撞作用,鋇離子和氧離子的動量相互耦合,鋇離子平行于磁場的動能傳遞給氧離子,促使氧離子沿著磁場向兩側(cè)運動,不斷將中心處的氧離子往外傳送,形成氧離子密度“空穴”;O+到達(dá)兩側(cè)后,又受到背景熱壓梯度的作用,阻止了其向外的膨脹,最終造成氧離子在兩側(cè)的堆積,產(chǎn)生兩個密度尖峰,即所謂的“掃雪機(jī)效應(yīng)”[14],氧離子空穴處的形態(tài)分布與鋇離子云團(tuán)形態(tài)是一致的.5 s時,空穴區(qū)的氧離子密度比背景降低了27.6%,密度凸起處比背景高14.3%,鋇離子的數(shù)密度峰值達(dá)到 1.332×107cm–3.

    圖2 無中性風(fēng)場時,300 km 高度釋放 10 kg 鋇后鋇離子和氧離子的離子數(shù)密度分布 (x-y 平面) (a) Ba+,t=5 s;(b) Ba+,t=30 s;(c) Ba+,t=200 s;(d) O+,t=5 s;(e) O+,t=30 s;(f) O+,t=200 sFig.2.Density distribution of Ba+ and O+ (in x-y plane) after 10 kg barium released at 300 km while no neutral wind is considered:(a) Ba+,t=5 s;(b) Ba+,t=30 s;(c) Ba+,t=200 s;(d) O+,t=5 s;(e) O+,t=30 s;(f) O+,t=200 s.

    圖3 無中性風(fēng)場時,300 km 高度釋放 10 kg 鋇后鋇離子和氧離子的粒子數(shù)密度分布 (x-z 平面) (a) O+,t=5 s;(b) O+,t=30 s;(c) O+,t=200 s;(d) Ba+,t=5 s;(e) Ba+,t=30 s;(f) Ba+,t=200 sFig.3.Density distribution of Ba+ and O+ (in x-z plane) after 10 kg barium released at 300 km while no neutral wind is considered:(a) O+,t=5 s;(b) O+,t=30 s;(c) O+,t=200 s;(d) Ba+,t=5 s;(e) Ba+,t=30 s;(f) Ba+,t=200 s.

    3.2.2 考慮中性風(fēng)場的情況

    圖4給出了x方向存在1 km/s的中性風(fēng)時,300 km高度釋放10 kg鋇的模擬結(jié)果(x-z平面),釋放中心位于 (–15,0,0) km 處,由于帶電粒子的密度分布在空間上是柱對稱的,這里不再給出x-y平面的圖像.與無風(fēng)場的結(jié)果相比,中性風(fēng)場的存在使得釋放的中性原子有了一個水平方向的運動,生成的等離子體云團(tuán)不再是對稱的橢球狀結(jié)構(gòu),而是一個逆風(fēng)側(cè)的密度梯度“較陡”,順風(fēng)側(cè)的密度梯度較為“平緩”的不對稱結(jié)構(gòu).比較圖3和圖4可以發(fā)現(xiàn),中性風(fēng)場存在時,釋放相同質(zhì)量鋇原子在早期對背景氧離子的密度擾動更大,由于中性云的移動,擾動區(qū)域也有所增加.5 s 時,空穴區(qū)的氧離子密度比背景降低30.7%,密度凸起處比背景高17.7%,由于分布區(qū)域變廣,鋇離子的數(shù)密度峰值略微降低,為 1.223×107cm–3.

    圖4 存在 x 方向大小為 1 km/s的中性風(fēng)時,300 km 高度釋放 10 kg 鋇后鋇離子和氧離子的粒子數(shù)密度分布 (x-z 平面) (a) O+,t=5 s;(b) O+,t=30 s;(c) O+,t=200 s;(d) Ba+,t=5 s;(e) Ba+,t=30 s;(f) Ba+,t=200 sFig.4.Density distribution of Ba+ and O+ (in x-z plane) after 10 kg barium released at 300 km with a neutral wind of 1 km/s in the x direction:(a) O+,t=5 s;(b) O+,t=30 s;(c) O+,t=200 s;(d) Ba+,t=5 s;(e) Ba+,t=30 s;(f) Ba+,t=200 s.

    圖5中右側(cè)的綠色球表示中性鋇原子,藍(lán)色部分代表鋇離子,雖然鋇離子的運動在垂直磁場方向上受到磁場的束縛,但中性鋇原子不受磁場力的約束,可以自由地穿過磁場線,因此鋇離子云團(tuán)與中性云團(tuán)會慢慢分離,如圖5 所示,30 s時離子云和中性云的分離已經(jīng)較為明顯.中性鋇原子在背景中性風(fēng)的作用下向x軸正方向移動,由于移動過程中持續(xù)的光電離作用,鋇中性云運動后方形成了拉長

    圖5 鋇中性云團(tuán) (綠色)和離子云團(tuán) (藍(lán)色)在釋放后30 s時的三維分布示意圖Fig.5.Three-dimensional density distribution of barium neutral cloud (green sphere) and ion cloud (blue sphere) at 30 s after release.

    的離子結(jié)構(gòu)(藍(lán)色部分),較早產(chǎn)生的鋇離子已經(jīng)沿磁場方向被拉伸(藍(lán)色球左側(cè)沿z方向拉伸的區(qū)域).

    3.2.3 銫釋放結(jié)果

    10 kg Cs原子釋放后 5,30和200 s銫離子及背景氧離子密度分布的演化過程如圖6所示,除釋放物質(zhì)外,其他釋放參數(shù)與 3.2.2節(jié)中一致.雖然銫由于電離勢較低容易發(fā)生熱電離,但在300 km高度及有日照的條件下,其自身的熱電離與光電離相比是微不足道的,因此這里沒有對銫的熱電離過程進(jìn)行探討.定性地說,Cs+和 Ba+的膨脹特性以及對背景O+的擾動效應(yīng)是類似的,雖然銫比鋇原子的質(zhì)量小,但銫原子的極化率較鋇原子要大,因而銫的擴(kuò)散系數(shù)較小,鋇云的膨脹更為迅速,鋇離子云團(tuán)的覆蓋區(qū)域更廣;同時,由于銫的光電離率較大,釋放相同質(zhì)量下銫的離子產(chǎn)率更高,5 s時銫離子云的峰值數(shù)密度達(dá)到 2.248×107cm–3,是相同條件下產(chǎn)生的鋇離子云密度的近兩倍;此外,Cs+-O+的碰撞頻率較 Ba+-O+更大,Cs+的掃雪機(jī)效應(yīng)比Ba+掃雪機(jī)更強(qiáng),氧離子密度空穴和凸起更大,30 s時鋇和銫釋放產(chǎn)生的背景氧離子的最大擾動分別為74.2%和75.1%.

    圖7給出了背景中性風(fēng)場存在下10 kg銫和鋇在300 km高度釋放后,生成的等離子體云團(tuán)密度的最大值以及背景氧離子密度凸起最高點的值隨時間的變化情況對比.釋放初期,伴隨著持續(xù)的光電離,鋇離子和銫離子數(shù)密度持續(xù)增加,之后由于擴(kuò)散作用逐漸下降,銫離子的數(shù)密度最大值比鋇離子大,但隨時間的變化趨勢是一致的;隨著等離子體云團(tuán)密度的增加,增強(qiáng)的密度梯度使得越來越多的氧離子被推開,兩側(cè)密度凸起處的值不斷增大,隨后,密度梯度開始下降,氧離子擾動也逐漸回復(fù).

    圖6 存在 x 方向大小為 1 km/s 的中性風(fēng)時,300 km 高度釋放 10 kg 銫的粒子數(shù)密度分布 (x-z 平面) (a) O+,t=5 s;(b) O+,t=30 s;(c) O+,t=200 s;(d) Cs+,t=5 s;(e) Cs+,t=30 s;(f) Cs+,t=200 sFig.6.Density distribution of Cs+ and O+ (in x-z plane) after 10 kg cesium released at 300 km with a neutral wind of 1 km/s in the x direction:(a) O+,t=5 s;(b) O+,t=30 s;(c) O+,t=200 s;(d) Cs+,t=5 s;(e) Cs+,t=30 s;(f) Cs+,t=200 s.

    圖7 生成的等離子體云團(tuán)的密度最大值(a)和背景氧離子的最大擾動值(b)隨時間的變化Fig.7.The maximum density of artificial plasma cloud (a) and the maximum disturbance of background oxygen ion (a) versus time.

    4 結(jié) 論

    采用三維雙成分的流體模型研究了電離層人工等離子體云的演化,分別給出了鋇在有無中性風(fēng)場時的模擬結(jié)果,并對存在背景風(fēng)場時鋇和銫釋放的電離層擾動效應(yīng)進(jìn)行了對比.基于已有的模擬結(jié)果,得到以下結(jié)論:1)鋇和銫在電離層釋放后迅速向四周膨脹,伴隨著太陽輻射作用下的光電離,生成的等離子體云在磁場的束縛下逐漸由內(nèi)向外呈橢球形結(jié)構(gòu),同時,由于離子間的動量耦合,膨脹的等離子體云會推開背景氧離子,在釋放中心形成氧離子密度空洞,并在兩側(cè)產(chǎn)生兩個氧離子密度尖峰,背景O+空洞的形狀類似于生成的人工等離子體云團(tuán)的形狀;2)在不存在中性風(fēng)的情況下,等離子體云團(tuán)以平行磁場的運動為主,考慮背景水平風(fēng)場時,云和背景受擾區(qū)域除了沿磁場方向的運動外,還會有沿著風(fēng)速方向的移動,逆風(fēng)側(cè)的等離子體云團(tuán)密度梯度變陡,而順風(fēng)側(cè)密度梯度較為平緩,雖然離子云在橫向上的運動會受到磁場的約束,但中性云團(tuán)可以自由地穿過磁場線,因此離子云與中性云會慢慢分離;3)對比鋇和銫的模擬結(jié)果發(fā)現(xiàn),定性地說,Cs+和 Ba+的膨脹特性以及對背景O+的擾動效應(yīng)是類似的,由于銫的擴(kuò)散系數(shù)較小,鋇云的膨脹更為迅速,鋇離子云團(tuán)的覆蓋區(qū)域更廣,同時,由于銫的光電離率較大,釋放相同質(zhì)量下銫的離子產(chǎn)率更高,此外,Cs+-O+的碰撞頻率較Ba+-O+更大,Cs+的掃雪機(jī)效應(yīng)比Ba+掃雪機(jī)更強(qiáng),造成的氧離子密度擾動也更大.

    本文僅利用簡化的模型進(jìn)行了研究,在下一步的工作中,可以綜合考慮地磁傾角、磁偏角、背景電離層密度梯度、溫度梯度等因素對模型加以完善;此外,對不同釋放條件以及不同地球物理條件包括不同太陽和地磁活動、季節(jié)以及地域等條件下人工等離子體云團(tuán)的形成和演化也可以進(jìn)行更加深入、系統(tǒng)的研究.

    猜你喜歡
    云團(tuán)電離層風(fēng)場
    基于FLUENT的下?lián)舯┝魅S風(fēng)場建模
    巴蜀少年齊上一堂云團(tuán)課
    少先隊活動(2021年3期)2021-12-04 13:08:26
    一種電離層TEC格點預(yù)測模型
    Kalman濾波估算電離層延遲的一種優(yōu)化方法
    “最美風(fēng)場”的贏利法則
    能源(2017年8期)2017-10-18 00:47:39
    面陣探測下的污染云團(tuán)紅外光譜仿真?
    電離層對中高軌SAR影響機(jī)理研究
    葡萄牙現(xiàn)明亮橘色云團(tuán)似握著火球的拳頭
    潤·文摘(2016年4期)2016-07-13 04:13:03
    側(cè)向風(fēng)場中無人機(jī)的飛行研究
    Linux Shell語言在電離層解算中的應(yīng)用
    制服人妻中文乱码| av片东京热男人的天堂| 中亚洲国语对白在线视频| 亚洲av日韩在线播放| 国产成人精品无人区| 夜夜爽天天搞| 三级毛片av免费| 欧美日韩视频精品一区| 国产一区二区在线观看av| 大片电影免费在线观看免费| 国产成人精品久久二区二区免费| 亚洲欧美日韩另类电影网站| 成年版毛片免费区| 如日韩欧美国产精品一区二区三区| 国产一区二区三区在线臀色熟女 | 亚洲免费av在线视频| 免费人妻精品一区二区三区视频| 交换朋友夫妻互换小说| 国产精品av久久久久免费| 亚洲黑人精品在线| 久久久久国内视频| 制服诱惑二区| 国产国语露脸激情在线看| 婷婷丁香在线五月| 搡老乐熟女国产| 久久 成人 亚洲| 在线观看免费视频日本深夜| 久久av网站| 一级片免费观看大全| 成人国产av品久久久| 一个人免费看片子| 日韩欧美免费精品| 大码成人一级视频| 亚洲精品一卡2卡三卡4卡5卡| 丁香欧美五月| 国产精品久久久久久人妻精品电影 | 成人国产av品久久久| 国产一区二区三区综合在线观看| 18禁黄网站禁片午夜丰满| 制服人妻中文乱码| 国产高清视频在线播放一区| 视频区欧美日本亚洲| 一区在线观看完整版| 91大片在线观看| 国产精品1区2区在线观看. | 人人妻人人澡人人爽人人夜夜| 亚洲一卡2卡3卡4卡5卡精品中文| 黄色怎么调成土黄色| 成人特级黄色片久久久久久久 | av网站在线播放免费| 中文字幕精品免费在线观看视频| 男女无遮挡免费网站观看| 久久人妻福利社区极品人妻图片| 777米奇影视久久| 欧美人与性动交α欧美软件| 在线观看免费视频网站a站| 天天操日日干夜夜撸| 老司机福利观看| 久久精品亚洲熟妇少妇任你| 久久久精品免费免费高清| cao死你这个sao货| 国产日韩一区二区三区精品不卡| 少妇粗大呻吟视频| 成在线人永久免费视频| 亚洲国产欧美在线一区| 欧美精品亚洲一区二区| 美女高潮喷水抽搐中文字幕| 欧美日韩一级在线毛片| 久久精品国产综合久久久| 乱人伦中国视频| 国产成人精品在线电影| 狂野欧美激情性xxxx| 91九色精品人成在线观看| 国产精品影院久久| 亚洲欧美一区二区三区黑人| 色尼玛亚洲综合影院| 激情在线观看视频在线高清 | 国产免费福利视频在线观看| 法律面前人人平等表现在哪些方面| a级毛片黄视频| 亚洲国产看品久久| 亚洲av国产av综合av卡| 亚洲欧美日韩高清在线视频 | 日韩有码中文字幕| 99久久99久久久精品蜜桃| 日韩免费av在线播放| 91精品三级在线观看| 国产欧美日韩一区二区精品| 免费一级毛片在线播放高清视频 | 国产激情久久老熟女| 欧美日韩av久久| 另类亚洲欧美激情| 别揉我奶头~嗯~啊~动态视频| 啦啦啦在线免费观看视频4| 老熟女久久久| 亚洲av第一区精品v没综合| 亚洲国产欧美在线一区| 精品国产超薄肉色丝袜足j| 欧美激情高清一区二区三区| 女人精品久久久久毛片| 黄色成人免费大全| 国产精品.久久久| 日韩有码中文字幕| 久久精品亚洲熟妇少妇任你| 欧美黑人欧美精品刺激| 老司机午夜十八禁免费视频| 老司机在亚洲福利影院| 国产精品久久久人人做人人爽| 丝袜人妻中文字幕| 50天的宝宝边吃奶边哭怎么回事| 国产精品影院久久| 男女边摸边吃奶| av天堂久久9| 国产精品免费视频内射| 欧美日本中文国产一区发布| 在线av久久热| 黄色视频在线播放观看不卡| 男女床上黄色一级片免费看| 黄色片一级片一级黄色片| av国产精品久久久久影院| 亚洲国产成人一精品久久久| 狠狠狠狠99中文字幕| 亚洲欧洲日产国产| 日韩免费av在线播放| 黄色视频,在线免费观看| 色综合欧美亚洲国产小说| 麻豆成人av在线观看| 国产在线视频一区二区| 成人三级做爰电影| 女人被躁到高潮嗷嗷叫费观| 亚洲欧美一区二区三区黑人| 成人18禁高潮啪啪吃奶动态图| 亚洲精品中文字幕在线视频| 国产亚洲欧美在线一区二区| 最近最新免费中文字幕在线| 后天国语完整版免费观看| 日韩大片免费观看网站| 美女扒开内裤让男人捅视频| 亚洲av成人一区二区三| 建设人人有责人人尽责人人享有的| 精品亚洲成a人片在线观看| 在线天堂中文资源库| 一边摸一边抽搐一进一小说 | 欧美另类亚洲清纯唯美| 国产成人免费无遮挡视频| 亚洲熟女毛片儿| 日韩视频在线欧美| 久久国产精品大桥未久av| 国产熟女午夜一区二区三区| 久久午夜综合久久蜜桃| 亚洲欧美精品综合一区二区三区| 国产欧美日韩精品亚洲av| 极品人妻少妇av视频| 久久久久网色| 中文字幕人妻丝袜一区二区| 少妇的丰满在线观看| 777久久人妻少妇嫩草av网站| 欧美老熟妇乱子伦牲交| 菩萨蛮人人尽说江南好唐韦庄| 免费日韩欧美在线观看| 国产在线免费精品| 免费av中文字幕在线| 国产无遮挡羞羞视频在线观看| 女性生殖器流出的白浆| 午夜激情久久久久久久| 久久久久久人人人人人| 麻豆av在线久日| 脱女人内裤的视频| 久热这里只有精品99| 一区在线观看完整版| 国产在视频线精品| 天堂动漫精品| 亚洲国产中文字幕在线视频| 日韩欧美国产一区二区入口| 久久99一区二区三区| 色尼玛亚洲综合影院| 欧美亚洲日本最大视频资源| 亚洲性夜色夜夜综合| 91麻豆av在线| 丝袜在线中文字幕| 免费女性裸体啪啪无遮挡网站| 欧美成人免费av一区二区三区 | 国产深夜福利视频在线观看| 久久久久国内视频| 亚洲va日本ⅴa欧美va伊人久久| 91麻豆精品激情在线观看国产 | 国产av又大| 欧美乱码精品一区二区三区| 欧美黑人欧美精品刺激| 在线av久久热| 国产欧美日韩一区二区三区在线| 亚洲国产精品一区二区三区在线| 日日摸夜夜添夜夜添小说| 国产一区有黄有色的免费视频| 在线观看免费日韩欧美大片| 久久精品亚洲熟妇少妇任你| 久久精品人人爽人人爽视色| 91大片在线观看| 97在线人人人人妻| 午夜福利在线观看吧| 一级毛片精品| 欧美午夜高清在线| 国产精品秋霞免费鲁丝片| 国产精品98久久久久久宅男小说| 视频区欧美日本亚洲| 亚洲国产中文字幕在线视频| 亚洲美女黄片视频| 国产精品成人在线| 成年人午夜在线观看视频| 久久久水蜜桃国产精品网| 亚洲精品乱久久久久久| 中文字幕制服av| 亚洲黑人精品在线| 黑人操中国人逼视频| 国产免费av片在线观看野外av| 三级毛片av免费| 一二三四在线观看免费中文在| 欧美激情极品国产一区二区三区| 正在播放国产对白刺激| 国产一区有黄有色的免费视频| 香蕉国产在线看| 欧美国产精品va在线观看不卡| 国产欧美日韩精品亚洲av| 成人三级做爰电影| 丁香六月天网| 后天国语完整版免费观看| 在线 av 中文字幕| 国产一区二区三区综合在线观看| 久久ye,这里只有精品| 精品亚洲乱码少妇综合久久| 中国美女看黄片| 香蕉久久夜色| 19禁男女啪啪无遮挡网站| √禁漫天堂资源中文www| 男人舔女人的私密视频| 在线观看免费日韩欧美大片| 亚洲va日本ⅴa欧美va伊人久久| 国产一区二区三区综合在线观看| 午夜成年电影在线免费观看| 国产亚洲精品第一综合不卡| 99香蕉大伊视频| 国产精品国产高清国产av | 精品一区二区三区四区五区乱码| 丰满饥渴人妻一区二区三| 狠狠婷婷综合久久久久久88av| 国产精品国产av在线观看| 久9热在线精品视频| 五月天丁香电影| 19禁男女啪啪无遮挡网站| 国产精品欧美亚洲77777| 成人av一区二区三区在线看| 嫁个100分男人电影在线观看| 一本综合久久免费| 精品久久久久久久毛片微露脸| 久久人妻av系列| 欧美精品高潮呻吟av久久| 亚洲色图 男人天堂 中文字幕| 国产精品国产高清国产av | 国产极品粉嫩免费观看在线| 麻豆乱淫一区二区| 大型av网站在线播放| 中文字幕av电影在线播放| 国产欧美日韩一区二区三| 亚洲国产av新网站| 久久久久久人人人人人| 国产av又大| 老司机在亚洲福利影院| 满18在线观看网站| av电影中文网址| 十八禁网站免费在线| 欧美国产精品一级二级三级| avwww免费| 精品一区二区三区四区五区乱码| 999久久久国产精品视频| 麻豆av在线久日| 老司机深夜福利视频在线观看| 亚洲精品美女久久久久99蜜臀| 国产亚洲欧美精品永久| 国产日韩欧美在线精品| 一本色道久久久久久精品综合| 黄色 视频免费看| 男女下面插进去视频免费观看| 亚洲av成人一区二区三| 亚洲中文av在线| 国产人伦9x9x在线观看| 国产不卡一卡二| 久久亚洲真实| 亚洲情色 制服丝袜| 麻豆乱淫一区二区| a级片在线免费高清观看视频| 国产精品偷伦视频观看了| 高清在线国产一区| 悠悠久久av| 国产精品自产拍在线观看55亚洲 | 久久久国产一区二区| av不卡在线播放| 1024香蕉在线观看| 操出白浆在线播放| 一级毛片精品| 日本欧美视频一区| 18在线观看网站| 欧美日韩亚洲国产一区二区在线观看 | 日本一区二区免费在线视频| 久久精品国产综合久久久| 精品一品国产午夜福利视频| 欧美成狂野欧美在线观看| 在线观看舔阴道视频| 国产淫语在线视频| 免费人妻精品一区二区三区视频| 免费在线观看日本一区| videos熟女内射| 久久久久国内视频| 亚洲成人手机| 成年动漫av网址| 亚洲人成电影免费在线| 女人被躁到高潮嗷嗷叫费观| 别揉我奶头~嗯~啊~动态视频| 黄色成人免费大全| 视频在线观看一区二区三区| 一级a爱视频在线免费观看| 日韩视频一区二区在线观看| 热re99久久国产66热| 久久久水蜜桃国产精品网| 两个人看的免费小视频| 高清视频免费观看一区二区| 久久国产精品人妻蜜桃| 国产伦理片在线播放av一区| 国产精品一区二区在线不卡| 他把我摸到了高潮在线观看 | 精品人妻熟女毛片av久久网站| 最新美女视频免费是黄的| 国产免费av片在线观看野外av| 午夜激情久久久久久久| 一区福利在线观看| 亚洲精品国产一区二区精华液| 国产一区有黄有色的免费视频| 亚洲成a人片在线一区二区| 丝瓜视频免费看黄片| 欧美人与性动交α欧美精品济南到| 一区二区三区激情视频| 青草久久国产| 国产成人欧美| 香蕉国产在线看| 激情视频va一区二区三区| 日韩大码丰满熟妇| 精品一区二区三区视频在线观看免费 | 美女视频免费永久观看网站| tocl精华| 老司机靠b影院| 国产黄频视频在线观看| 999精品在线视频| 欧美成狂野欧美在线观看| 黄网站色视频无遮挡免费观看| 一边摸一边抽搐一进一小说 | 另类亚洲欧美激情| 成人18禁高潮啪啪吃奶动态图| 国产不卡av网站在线观看| 多毛熟女@视频| 久久精品aⅴ一区二区三区四区| 亚洲avbb在线观看| 国产精品国产av在线观看| 男女高潮啪啪啪动态图| 精品乱码久久久久久99久播| 9色porny在线观看| 中文字幕人妻丝袜一区二区| 真人做人爱边吃奶动态| 欧美性猛交╳xxx乱大交人| 国产成人精品久久二区二区免费| 欧美性猛交╳xxx乱大交人| 国产熟女xx| 国产精品女同一区二区软件 | 午夜免费激情av| 久久国产精品影院| 亚洲成av人片免费观看| 国产高潮美女av| 国产精品美女特级片免费视频播放器 | 国产成人精品无人区| 又黄又爽又免费观看的视频| 色精品久久人妻99蜜桃| 国产精品久久久久久久电影 | 脱女人内裤的视频| 亚洲无线在线观看| 亚洲国产日韩欧美精品在线观看 | 国产69精品久久久久777片 | 啦啦啦观看免费观看视频高清| 国产高清视频在线观看网站| 国产激情偷乱视频一区二区| av在线天堂中文字幕| 久久精品亚洲精品国产色婷小说| 精品免费久久久久久久清纯| 国产黄片美女视频| 又大又爽又粗| 动漫黄色视频在线观看| 99在线视频只有这里精品首页| 一进一出抽搐gif免费好疼| 国产成人aa在线观看| 日韩欧美国产一区二区入口| 久久久国产成人精品二区| 精华霜和精华液先用哪个| 欧美极品一区二区三区四区| 久久天堂一区二区三区四区| 看免费av毛片| 天天添夜夜摸| 欧美激情久久久久久爽电影| 99riav亚洲国产免费| 国产午夜福利久久久久久| 欧美乱妇无乱码| 国产高清视频在线播放一区| 97超级碰碰碰精品色视频在线观看| 日韩成人在线观看一区二区三区| 亚洲国产中文字幕在线视频| 岛国视频午夜一区免费看| 俄罗斯特黄特色一大片| 欧美乱色亚洲激情| 69av精品久久久久久| 国产精品国产高清国产av| 亚洲精华国产精华精| 日本成人三级电影网站| 久久九九热精品免费| 国内精品美女久久久久久| 久久久色成人| av天堂在线播放| 最新在线观看一区二区三区| 五月伊人婷婷丁香| 国产综合懂色| 琪琪午夜伦伦电影理论片6080| 国产亚洲精品综合一区在线观看| 18禁美女被吸乳视频| 日韩精品中文字幕看吧| 亚洲五月婷婷丁香| 嫩草影视91久久| 成人鲁丝片一二三区免费| 丁香欧美五月| 校园春色视频在线观看| 真人一进一出gif抽搐免费| 变态另类成人亚洲欧美熟女| 精品一区二区三区视频在线 | 俄罗斯特黄特色一大片| 精品福利观看| 禁无遮挡网站| 神马国产精品三级电影在线观看| 色综合欧美亚洲国产小说| 久久久成人免费电影| svipshipincom国产片| 日本 av在线| 国产v大片淫在线免费观看| 亚洲中文字幕一区二区三区有码在线看 | 黄色片一级片一级黄色片| 久久精品91蜜桃| 99国产精品一区二区蜜桃av| 一级毛片精品| 亚洲美女视频黄频| 很黄的视频免费| 精品国产超薄肉色丝袜足j| 免费搜索国产男女视频| 欧美高清成人免费视频www| www.999成人在线观看| 99久久综合精品五月天人人| 国产高潮美女av| www.www免费av| 成人国产综合亚洲| 久久精品综合一区二区三区| 色综合欧美亚洲国产小说| 大型黄色视频在线免费观看| 天天躁日日操中文字幕| 日本熟妇午夜| 在线观看午夜福利视频| 男人舔女人的私密视频| 成人三级做爰电影| 波多野结衣高清作品| netflix在线观看网站| 久久久久久久精品吃奶| 国内精品美女久久久久久| 在线看三级毛片| 午夜a级毛片| 亚洲在线观看片| a级毛片a级免费在线| 熟妇人妻久久中文字幕3abv| 俺也久久电影网| 久久久久久九九精品二区国产| 最新在线观看一区二区三区| 欧美午夜高清在线| 国产成人精品久久二区二区免费| 免费人成视频x8x8入口观看| 国产极品精品免费视频能看的| 国内精品久久久久久久电影| 国内揄拍国产精品人妻在线| 脱女人内裤的视频| 91在线精品国自产拍蜜月 | 丰满人妻一区二区三区视频av | 成年版毛片免费区| 亚洲国产色片| 国产一区在线观看成人免费| 两个人看的免费小视频| 国产精品女同一区二区软件 | 亚洲欧美激情综合另类| 婷婷六月久久综合丁香| 看片在线看免费视频| 国产成人精品久久二区二区91| 欧美三级亚洲精品| 亚洲一区高清亚洲精品| 黄片小视频在线播放| 波多野结衣巨乳人妻| 两个人视频免费观看高清| 熟女电影av网| 两人在一起打扑克的视频| 国产精品久久久人人做人人爽| 日日摸夜夜添夜夜添小说| av视频在线观看入口| 色哟哟哟哟哟哟| 夜夜夜夜夜久久久久| 色综合亚洲欧美另类图片| 波多野结衣巨乳人妻| 亚洲国产精品成人综合色| 亚洲成人久久性| 中出人妻视频一区二区| 十八禁网站免费在线| 超碰成人久久| 青草久久国产| 美女大奶头视频| 久久午夜综合久久蜜桃| 俄罗斯特黄特色一大片| 日本黄色视频三级网站网址| e午夜精品久久久久久久| 亚洲激情在线av| 9191精品国产免费久久| 性色avwww在线观看| av天堂在线播放| 天堂影院成人在线观看| 久久中文字幕人妻熟女| 久久精品91蜜桃| 亚洲一区高清亚洲精品| 男女床上黄色一级片免费看| 99国产精品一区二区三区| 一本一本综合久久| 巨乳人妻的诱惑在线观看| 中文字幕久久专区| 亚洲aⅴ乱码一区二区在线播放| 国产一区二区三区视频了| 免费看日本二区| 亚洲av美国av| 国产成+人综合+亚洲专区| 给我免费播放毛片高清在线观看| 女人高潮潮喷娇喘18禁视频| 久久国产精品人妻蜜桃| 一个人免费在线观看电影 | 国产精品亚洲av一区麻豆| 精品国产美女av久久久久小说| 成人精品一区二区免费| 色精品久久人妻99蜜桃| 女人高潮潮喷娇喘18禁视频| 国产高清三级在线| 91av网站免费观看| 精品一区二区三区视频在线观看免费| 免费看日本二区| 亚洲精品美女久久av网站| 欧美不卡视频在线免费观看| 午夜激情欧美在线| 国产精品一及| 两性午夜刺激爽爽歪歪视频在线观看| 精品一区二区三区视频在线 | 午夜福利在线观看吧| 人人妻,人人澡人人爽秒播| 高清在线国产一区| 久久中文字幕一级| 999精品在线视频| 精品99又大又爽又粗少妇毛片 | av天堂中文字幕网| 国产亚洲欧美98| 美女高潮喷水抽搐中文字幕| 亚洲最大成人中文| 国产精品美女特级片免费视频播放器 | 两人在一起打扑克的视频| 欧美黄色片欧美黄色片| 亚洲自偷自拍图片 自拍| 欧美成狂野欧美在线观看| 亚洲国产看品久久| 中文字幕av在线有码专区| 俺也久久电影网| 97人妻精品一区二区三区麻豆| 国产精品电影一区二区三区| 国产精品久久视频播放| 亚洲成人精品中文字幕电影| 99精品在免费线老司机午夜| 最近最新免费中文字幕在线| 亚洲乱码一区二区免费版| 免费观看精品视频网站| 91久久精品国产一区二区成人 | 色综合亚洲欧美另类图片| 国产伦在线观看视频一区| 欧美乱妇无乱码| 久久精品夜夜夜夜夜久久蜜豆| 人人妻人人看人人澡| a在线观看视频网站| 五月伊人婷婷丁香| 中文字幕高清在线视频| 又大又爽又粗| 巨乳人妻的诱惑在线观看| 中文字幕高清在线视频| 久久性视频一级片| 男女下面进入的视频免费午夜| 黑人操中国人逼视频| 亚洲欧美日韩高清专用| av欧美777| 国产亚洲精品一区二区www| 午夜精品久久久久久毛片777| 丰满的人妻完整版| 亚洲va日本ⅴa欧美va伊人久久| 日韩有码中文字幕| 国产精品免费一区二区三区在线| 又黄又爽又免费观看的视频| 亚洲人成伊人成综合网2020| 波多野结衣高清作品| 亚洲专区国产一区二区| 精品久久蜜臀av无|