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

    基于深度學習的鉆孔輻射壓離子加速建模*

    2023-10-06 07:03:56張普渡王偉權(quán)李哲民張資旋王葉晨周泓宇銀燕
    物理學報 2023年18期
    關(guān)鍵詞:模型

    張普渡 王偉權(quán)? 李哲民 張資旋 王葉晨 周泓宇 銀燕?

    1) (國防科技大學理學院物理系,長沙 410073)

    2) (國防科技大學理學院數(shù)學系,長沙 410073)

    超短超強激光脈沖與固體靶相互作用可通過鉆孔輻射壓加速機制產(chǎn)生百MeV 量級的離子束,離子束的品質(zhì)強烈依賴于激光和靶的作用參量.本文以近400 組激光驅(qū)動固體靶的粒子模擬結(jié)果作為數(shù)據(jù)集,以激光強度、靶密度、靶厚和離子質(zhì)量作為輸入?yún)⒘?基于全連接神經(jīng)網(wǎng)絡建立了一個離子峰值能量和截止能量連續(xù)映射模型.該模型用較為稀疏的參量取值獲得了較大參量范圍內(nèi)的分析結(jié)果,大大減少了多維參量大范圍掃參的計算量.基于連續(xù)映射模型的結(jié)果,得到了鉆孔輻射壓加速機制下離子峰值能量的修正公式和截止能量的擬合公式,可為激光離子加速的實驗設計提供重要參考.

    1 引言

    隨著超短超強激光脈沖與等離子體相互作用研究的不斷深入,基于超強激光驅(qū)動的離子加速逐漸成為激光等離子體領(lǐng)域的重要研究方向,在離子束驅(qū)動快點火[1,2]、束靶中子源[3,4]、溫稠密物質(zhì)加熱[5-7]、質(zhì)子照相[8,9]等領(lǐng)域得到了廣泛研究與應用.常見的激光驅(qū)動的離子加速機制包含輻射壓加速(radiation pressure acceleration,RPA)[10,11]、靶背法向鞘層加速(target normal sheath acceleration,TNSA)[12,13]、靶破燒蝕加速[14,15]、無碰撞激波加速[16,17]等.其中,RPA 機制由于其能量轉(zhuǎn)化效率高、離子束單能性好以及可加速重離子等特點,在實驗上備受關(guān)注.根據(jù)固體靶厚度的不同,RPA機制可分為薄靶條件下的光帆加速(light-sail RPA,LS-RPA)[18,19]和厚靶條件下的鉆孔加速(hole-boring RPA,HB-RPA)[20,21].在LS-RPA 機制下,離子與電子將被壓縮為“等離子體飛鏡(plasma mirror)”結(jié)構(gòu)被激光脈沖整體推動加速,但這種機制對激光脈沖對比度要求較高,且靶厚度和激光脈沖強度之間有嚴格的匹配關(guān)系[18,19,22].而對于HB-RPA 機制,激光和靶的參數(shù)要求相對寬松,更容易在實驗中實現(xiàn)[23],且可在較高激光強度下對重離子實現(xiàn)有效加速[24].

    在基于HB-RPA 機制的數(shù)值模擬和實驗設計中,分析研究對象與參量的關(guān)聯(lián)時需進行參數(shù)掃描.然而,激光和靶參數(shù)的選取范圍較為寬泛,通過開展大量實驗或數(shù)值模擬進行線性掃參的代價過于高昂.同時,不同掃描參量與目標參量間可能存在一定映射關(guān)聯(lián),而通過線性掃參只能得到有限的孤立數(shù)據(jù)點,基于有限的數(shù)據(jù)建立較為可靠的連續(xù)映射關(guān)系存在較大困難.因此,有必要尋找新的研究參量映射規(guī)律的手段.近幾年內(nèi),依賴于數(shù)據(jù)學習的神經(jīng)網(wǎng)絡算法[25,26]被逐漸應用在科研領(lǐng)域,例如托卡馬克裝置中磁場的控制優(yōu)化[27]、慣性約束聚變實驗結(jié)果對實驗參數(shù)的敏感度研究[28]、電子加速能譜的自動識別與分類[29]、高輻射環(huán)境下測得的電子能譜圖的去噪分析[30]、預測粒子加速器內(nèi)粒子束縱向相空間演化趨勢[31]等.在眾多類型的神經(jīng)網(wǎng)絡算法中,卷積神經(jīng)網(wǎng)絡(convolutional neural network,CNN)主要面向二維矩陣數(shù)據(jù)的處理,在圖像處理和識別等領(lǐng)域應用廣泛[30].自組織映射神經(jīng)網(wǎng)絡(self-organizing map,SOM)的訓練過程為無監(jiān)督的競爭學習,適用于高維數(shù)據(jù)的聚類分析[29].全連接構(gòu)型人工神經(jīng)網(wǎng)絡(fully connected neural network,FCNN)可以建立起數(shù)個變量關(guān)于單一映射量的連續(xù)模型,適用于數(shù)據(jù)回歸類問題的建模分析.Djordjevi?等[32]基于FCNN 神經(jīng)網(wǎng)絡算法對TNSA 加速機制的多維參數(shù)進行了建模分析,得到離子截止能量的連續(xù)映射模型,并給出TNSA 機制下離子能量增大速率的連續(xù)映射模型[33].

    本文所研究的HB-RPA 機制下的參數(shù)映射也是一個數(shù)據(jù)回歸類問題,采用FCNN 神經(jīng)網(wǎng)絡較為適合.首先將有限的數(shù)值模擬結(jié)果整理為數(shù)據(jù)集輸入FCNN,然后在數(shù)據(jù)迭代的過程中,通過神經(jīng)元權(quán)重的反向傳遞,把目標參量關(guān)于輸入?yún)⒘康淖兓卣饔涗浽谏窠?jīng)網(wǎng)絡中,從而得到連續(xù)映射模型.這一所得模型可較好反映出掃參空間中所有數(shù)據(jù)點的映射結(jié)果,并可在訓練范圍外進行推廣,進而對更大范圍的結(jié)果進行預測.本文利用FCNN建立了HB-RPA 機制下離子峰值能量和截止能量關(guān)于激光強度、靶密度、靶厚度和離子質(zhì)量這4 個輸入?yún)⒘康倪B續(xù)映射模型.根據(jù)該模型與HB-RPA的理論公式及PIC 數(shù)值模擬結(jié)果的對比結(jié)果,對離子峰值能量公式進行修正,同時給出離子截止能量的經(jīng)驗公式,實現(xiàn)了對HB-RPA 機制下離子峰值能量和截止能量的定量描述.

    2 物理模型

    當激光強度與靶厚度近似滿足a<πdn/(ncλ0)關(guān)系時,HB-RPA 機制將占主導地位[34],其中,a為歸一化激光強度,與激光強度I0的關(guān)系滿足a=為激光波長,d為靶厚度為等離子體臨界密度,me為電子質(zhì)量,ω為激光頻率,e為電子電荷量.一維(one-dimension,1D)情況下,一束強度恒定為I0的圓極化激光正入射到均勻密度為n的等離子體,光束在等離子體表面被完全反射.由穩(wěn)態(tài)加速離子所滿足的動量平衡方程,可以得到實驗室系中相對論條件下動量平衡的穩(wěn)態(tài)離子能量表達式[20]:

    本文通過PIC 程序EPOCH[35]對HB-RPA 過程進行了約400 次1D 數(shù)值模擬.在模擬中改變激光脈沖的歸一化激光強度a,等離子體靶厚度d,等離子體靶密度n和離子質(zhì)量mi,模擬參數(shù)選取范圍均滿足HB-RPA 離子加速機制要求.由于靶厚度d是決定HB-RPA 加速機制是否起效的重要參量,將d作為輸入?yún)⒘考尤氲缴窠?jīng)網(wǎng)絡的訓練中,可以更好地保證訓練所得模型的可靠性.在1D 數(shù)值模擬中,激光時間脈沖設置為梯形包絡,如下式所示:

    其中T0為激光周期,I0為激光峰值強度.等離子體為單一離子組分的平面靶,靶離子為完全電離.模擬盒子范圍為 0 ≤x≤90λ0,靶前表面位置處于10λ0處.激光脈沖由左邊界入射.模擬空間分辨率為 ?x=0.01λ0,每個網(wǎng)格內(nèi)各組分的宏粒子數(shù)目為50 個.

    圖1(a)和圖1(c) 分別為He2+離子和C6+離子的x-v相空間圖,其中,v為離子速度.x-v相空間圖展現(xiàn)了HB-RPA 機制的典型振蕩特征[20,21].在HB-RPA 機制中,激光推動電子形成壓縮層,由此產(chǎn)生的電荷分離場把大量靶內(nèi)離子加速到穩(wěn)態(tài)速度,如紅色虛線框所示.加速電場的峰值位于電子壓縮層的前表面,隨時間振蕩[20],離子在振蕩電場作用下形成形似魚骨的x-v相空間結(jié)構(gòu).少量滯后的離子被加速電場的上升沿推動至峰值電場區(qū)域,進而被加速至更高能量,如藍色虛線框所示.圖1(b)和圖1(d)分別為He2+離子和C6+離子在t=70T0時刻的能譜圖,此時離子相空間、離子峰值能量和最大離子能量都已趨于穩(wěn)定.在能譜圖上略去了能量位于能譜中最低0.5%部分的低能離子,呈現(xiàn)清晰的準單能結(jié)構(gòu).對于每組數(shù)值模擬,通過診斷離子的x-v相空間分布來確定該組模擬中HB-RPA 加速機制占主導地位,通過離子能譜圖統(tǒng)計離子束的峰值能量和截止能量.舍棄掉不符合HB-RPA 加速機制的數(shù)據(jù)點,最終得到400 組數(shù)據(jù)集樣本數(shù)據(jù)點.

    圖1 在 t=70T0 時刻,(a) He2+離子的 x-v 相空間分布,(b) He2+離子的能譜分布,(c) C6+離子的 x-v 相空間分布,(d) C6+離子的能譜分布.其中,激光強度 a=100 ,靶密度 n=25nc ,靶厚 d=25λ0 ;紅色虛線框內(nèi)為峰值能量所對應的離子,藍色虛線框內(nèi)為截止能量所對應的離子Fig.1.At t=70T0 ,(a) the x-v phase space diagram of He2+,(b) the energy spectrum of He2+,(c) the x-v phase space diagram of C6+,(d) the energy spectrum of C6+,with d=25λ0 ,n=25nc and a=100 .The red dashed box circles the ions with peak energy,and the blue dashed box circles the ions with cutoff energy.

    3 數(shù)據(jù)準備和搭建神經(jīng)網(wǎng)絡

    3.1 數(shù)據(jù)集生成

    整體數(shù)據(jù)集的采樣采用系統(tǒng)抽樣方式,在一定的參數(shù)范圍內(nèi)均勻間隔取值,通過離子種類將數(shù)據(jù)集劃分為質(zhì)子集、He2+集、C6+集、O8+集和測試集共5 個集合.其中,質(zhì)子集合經(jīng)過數(shù)據(jù)篩選后總共記錄了231 次數(shù)值模擬結(jié)果,具體的輸入?yún)?shù)取值為: 靶厚d在 [7.5λ0,35λ0] 范圍內(nèi)間隔 2.5λ0均勻取值,靶密度n在 [12.5nc,40nc] 范圍內(nèi)間隔 2.5nc均勻取值,激光強度分別取值4×1021,6×1021,8×1021,1×1022,2×1022,3×1022W/cm2.He2+離子集合包含49 次數(shù)值模擬的結(jié)果,C6+離子與O8+離子分別包含50 次數(shù)值模擬的結(jié)果,3 種離子的靶厚d取值為 15λ0和 25λ0,靶密度n在 [10nc,30nc] 范圍內(nèi)間隔 5nc均勻取值,激光強度取值為6×1021,8×1021,1×1022,2×1022,3×1022W/cm2.最后一個集合為獨立于其他數(shù)據(jù)的測試集,離子種類包括D+,T+,He2+,B4+,C6+,F9+,Mg12+和Al13+離子,在每組數(shù)值模擬中,模擬參數(shù)d,n和a均獨立取值,共20 個數(shù)據(jù)點.記錄每組模擬結(jié)果的離子峰值能量Ep與截止能量Em作為訓練參量,具體取值分布如表1 所列.圖2(a)–(d)更直觀地給出了4 種離子的參數(shù)取值分布,并用散點顏色標注了該數(shù)據(jù)點所對應的離子峰值能量.為控制輸入?yún)⒘康臄?shù)量級相近,靶厚參數(shù)d以λ0歸一化,靶密度參數(shù)n以臨界密度nc歸一化,激光強度以歸一化激光強度a輸入,離子質(zhì)量以原子質(zhì)量數(shù)A輸入.

    表1 數(shù)值模擬數(shù)據(jù)集參數(shù)取值分布Table 1. Simulation datasets prepared for neural network training.

    圖2 離子峰值能量隨輸入?yún)?shù) d ,n 和 a 分布的三維散點圖 (a) 質(zhì)子;(b) He2+ 離子;(c) C6+ 離子;(d) O8+離子Fig.2.Scatter plot depiction of ion peak energy data ensembles as a function of input parameters d,n and a: (a) Proton case;(b) He2+,(c) C6+;(d) O8+.

    3.2 神經(jīng)網(wǎng)絡搭建

    本文基于Pytorch 中的nn 模塊來搭建深度神經(jīng)網(wǎng)絡,所采用的網(wǎng)絡類型為FCNN,通過測試與對比訓練結(jié)果對神經(jīng)網(wǎng)絡結(jié)構(gòu)與超參數(shù)的選取進行調(diào)整,最終確定神經(jīng)網(wǎng)絡結(jié)構(gòu)如圖3 所示.經(jīng)過篩選和后處理的樣本數(shù)據(jù)被輸入到輸入層,4 個神經(jīng)元分別對應輸入?yún)?shù)d,n,a,A.經(jīng)過6 個隱藏層后,結(jié)果收斂至輸出層的單一輸出量.Ep與Em采用相同的神經(jīng)網(wǎng)絡結(jié)構(gòu).

    圖3 訓練所使用的神經(jīng)網(wǎng)絡結(jié)構(gòu)Fig.3.Neural network architecture used in the following training.

    神經(jīng)網(wǎng)絡訓練參數(shù)的選取與所得到代理模型(surrogates ensemble,SE)[36]的泛化能力密切相關(guān).我們所使用的神經(jīng)網(wǎng)絡采用均方方差損失(mean square error,MSE)函數(shù),采用nn 模塊默認的L2正則化,其中L2=2×10-5.神經(jīng)元激活函數(shù)選取為Sigmoid 與PReLU,這兩種激活函數(shù)皆為單調(diào)函數(shù),模型映射曲線在有限的預測范圍內(nèi)保持緩慢平滑變化,這一特征與HB-RPA 理論公式相契合.訓練結(jié)果表明,我們所確定的神經(jīng)網(wǎng)絡結(jié)構(gòu)可使模型收斂至接近理論預測結(jié)果,同時保持了較好的平滑程度和泛化性,這使得分析參量的訓練和預測結(jié)果都具有一定的可靠性.

    3.3 訓練與測試

    針對模擬數(shù)據(jù)所展開的訓練類型均為監(jiān)督學習,在給定輸入數(shù)據(jù)上訓練神經(jīng)網(wǎng)絡,訓練過程中優(yōu)化器選擇基于隨機梯度下降的Adam 優(yōu)化器,損失函數(shù)選擇MSE.本文將數(shù)據(jù)集劃分為380 組數(shù)據(jù)的訓練集和20 組數(shù)據(jù)的測試集,訓練集與測試集皆為固定.訓練數(shù)據(jù)在神經(jīng)網(wǎng)絡上迭代次數(shù)固定為8000 次,學習率固定為1×10-4,損失函數(shù)誤差下降至1×10-4時,模型趨于收斂,且平滑程度與擬合程度良好.

    4 研究結(jié)果

    4.1 神經(jīng)網(wǎng)絡集合

    在初始化權(quán)值的訓練條件下,每次的訓練結(jié)果會略有不同.由于本文所使用的參數(shù)空間相對稀疏,每次訓練所得模型的差異將會更加明顯.為此,本文對10 個神經(jīng)網(wǎng)絡同時進行訓練,得到一個神經(jīng)網(wǎng)絡集合,并根據(jù)模型收斂時的損失函數(shù)取值對訓練結(jié)果做加權(quán)平均,從而生成訓練結(jié)果的代理模型.取每個模型最后2000 次迭代的損失函數(shù)均值的倒數(shù)作為該模型的權(quán)值,對所有神經(jīng)網(wǎng)絡模型加權(quán)平均,由模型權(quán)值進行歸一化,即可得到

    式中,代理模型SE 是所有訓練模型Si的均值;wi為每個訓練模型的權(quán)值,為最后2000 次迭代損失函數(shù)均值的倒數(shù);x為輸入?yún)?shù)的取值;N=10 表示神經(jīng)網(wǎng)絡模型訓練個數(shù).為了衡量代理模型在一定取值范圍內(nèi)的可變性,定義了標準差σ:

    4.2 代理模型映射

    利用神經(jīng)網(wǎng)絡訓練結(jié)果和加權(quán)所得的代理模型,在一定范圍內(nèi)的參數(shù)空間上對輸入?yún)?shù)進行連續(xù)掃參,可得到目標參量關(guān)于不同輸入?yún)?shù)的映射關(guān)系,并與測試集數(shù)據(jù)點進行比照.

    圖4 為根據(jù)神經(jīng)網(wǎng)絡代理模型所繪出的離子峰值能量與截止能量關(guān)于靶密度n和激光強度a的二維連續(xù)映射圖.從圖4(a) 可看出,當激光強度一定時,峰值能量Ep隨靶密度的增大而減小;當靶密度固定時,峰值能量Ep隨激光強度的增大而增大.當激光強度與靶密度滿足保持不變時,離子峰值能量接近為定值,峰值能量Ep映射結(jié)果與(1)式理論結(jié)果符合很好.同時,圖4(b)給出的截止能量Em隨靶密度增大而減小,隨激光強度增大而增大,與峰值能量Ep呈現(xiàn)出一致的映射規(guī)律.

    圖4 He2+離子 (a) 峰值能量 Ep 與 (b) 截止能量 Em 的二維連續(xù)映射圖,其中靶厚 d=15λ0 ,參數(shù)映射范圍為63.2 ≤a ≤173.2和10nc ≤n ≤40ncFig.4.Two-dimensional continuous mapping of peak energy Ep (a) and cutoff energy Em (b) for He2+ over 63.2 ≤a ≤173.2 and 10nc ≤n ≤40nc with d=15λ0 .

    圖5 為離子峰值能量代理模型關(guān)于單個輸入?yún)⒘康挠成淝€和峰值能量理論曲線(由(1)式給出)的比照結(jié)果.圖5(a) 給出了He2+離子峰值能量隨靶密度n變化的連續(xù)映射曲線,訓練參數(shù)取值范圍為 10nc≤n≤30nc,預測參數(shù)取值范圍為30nc≤n≤50nc.由圖5(a) 所示,訓練參數(shù)取值范圍內(nèi),代理模型與理論曲線和測試集數(shù)據(jù)吻合程度相當好;在預測參數(shù)取值范圍內(nèi),模型映射值與測試集數(shù)據(jù)偏差小于10%,模型的可靠性可以得到保證.而當n>50nc時,模型映射值與理論值相差超過10%,模型的可靠性不好,所以這里沒有給出映射結(jié)果.若增大訓練參數(shù)取值范圍,則可在更大范圍內(nèi)保證模型預測的可靠性.

    圖5 (a) He2+離子 峰值 能量 Ep 關(guān)于靶密度 n 的映 射曲線,其中靶厚 d=15λ0 ,激光強度 a=100 ;(b) He2+離 子峰 值能量Ep關(guān)于激光強度 a 的映射曲線,其中靶厚 d=15λ0 ,靶 密度 n=20nc ;(c) He2+離子峰值 能量 Ep 關(guān)于 靶厚 d 的映射曲 線,其中 靶密度 n=20nc ,激光強度 a=100 ;(d)離子峰值 能量 Ep 關(guān)于離子質(zhì)量數(shù) A 的映 射曲 線,其中靶厚 d=15λ0 ,靶密度 n=20nc,激光強度 a=100 .圖中紅色實線為代理模型的映射曲線,模型標準差用紅色色塊填充于代理模型曲線兩側(cè).藍色虛線給出的是根據(jù)HB-RPA 機制理論公式所繪出的離子峰值能量隨離子密度變化的曲線,黑色實心數(shù)據(jù)點為測試集數(shù)據(jù)點,參數(shù)取值范圍與模型預測范圍用黑色虛線分隔Fig.5.(a) Parameter scan of He2+ peak energy Ep over target density n with d=15λ0 ,a=100 ;(b) parameter scan of He2+peak energy Ep over laser intensity a with d=15λ0 ,n=20nc ;(c) parameter scan of He2+ peak energy Ep over target thickness d with a=100 ,n=20nc ;(d) parameter scan of ion peak energy Ep over ion mass number A with a=100 ,n=20nc and d=15λ0 .The SE mapping curves are drawn with red solid,the theoretical curves are drawn with dashed,the red filled region indicates the standard deviation,and the untrained data from test subset are drawn with black dot.

    圖5(b)給出了He2+離子峰值能量與激光強度a的模型連續(xù)映射曲線.激光強度的訓練參數(shù)取值范圍為 77.5

    4.3 離子的峰值能量和截止能量

    4.3.1 峰值能量修正

    圖5(d) 展現(xiàn)出映射曲線與理論曲線存在偏差,離子峰值能量在質(zhì)量較低時略高于理論值,而在質(zhì)量較高時略低于理論值,模型曲線與理論曲線于離子質(zhì)量數(shù)A=5 附近相交.對比了更大范圍的數(shù)值模擬結(jié)果與理論值,發(fā)現(xiàn)這一偏差仍然存在.結(jié)合代理模型所給出的模型和測試集所得模擬數(shù)據(jù),對原理論公式進行修正,修正后的公式可描述離子質(zhì)量對離子峰值能量的影響:

    式中,A為離子質(zhì)量數(shù).經(jīng)過修正的模型曲線與測試集數(shù)據(jù)點對比如圖6 所示.

    圖6 在 d=15λ0 ,n=20nc 和 a=100 下,離子峰值 能量隨離子質(zhì)量的變化Fig.6.Ion peak energy scan with ion mass number A with d=15λ0,n=20nc and a=100 .

    如圖6 所示,當質(zhì)量數(shù)A>5 時,數(shù)值模擬結(jié)果與理論曲線開始出現(xiàn)偏移.當質(zhì)量數(shù)A>10 時,模擬結(jié)果與理論結(jié)果差異逐漸增大,誤差率接近10%左右.當采用基于離子質(zhì)量數(shù)的修正公式時,理論結(jié)果與數(shù)值模擬結(jié)果及模型映射結(jié)果符合很好.對于質(zhì)量數(shù)A>30 的中高Z離子,其較為復雜的電離過程將對離子的加速過程產(chǎn)生影響.在電離度較低的情況下,重離子將滯后于電子被推動形成的鉆孔活塞結(jié)構(gòu),此時HB-RPA 模型不再適用,離子峰值能量將顯著低于理論值[21].因此,這里基于映射結(jié)果給出的修正模型僅適用于A<30 的靶材料.在質(zhì)量數(shù)A>30 時,需增大數(shù)據(jù)集的質(zhì)量數(shù)取值范圍,并考慮靶材料的離化過程對加速過程的影響.

    4.3.2 截止能量映射

    當離子截止能量作為目標參量時,重復與離子峰值能量相同的訓練和掃參操作,可以得到連續(xù)映射曲線,如圖7 所示.

    圖7 (a) He2+離子截止能量 Em 關(guān)于靶密度 n 的映射曲線,其中靶厚 d=15λ0 ,激光強度 a=100 ;(b) He2+離子截止能量Em關(guān)于激光強度 a 的映射曲線,其中靶厚 d=15λ0 ,靶密度 n=20nc ;(c) He2+離子截止能量 Em 關(guān)于靶厚 d 的映射曲線,其中靶密度 n=20nc ,激 光強度 a=100 ;(d) 離子截止能量 Em 關(guān)于離子質(zhì)量數(shù) A 的映射曲線,其中靶厚 d=15λ0 ,靶密度 n=20nc,激光強度a=100Fig.7.(a) Parameter scan of He2+ cutoff energy Em over target density n with d=15λ0 ,a=100 ;(b) parameter scan of He2+cutoff energy Em over laser intensity a with d=15λ0 ,n=20nc ;(c) parameter scan of He2+ cutoff energy Em over target thickness d with a=100 ,n=20nc ;(d) parameter scan of ion cutoff energy Em over ion mass number A with a=100,n=20nc and d=15λ0 .

    圖7 給出了He2+離子的截止能量關(guān)于靶密度n、激光強度、靶厚和離子質(zhì)量的映射關(guān)系.結(jié)果表明,離子截止能量隨n的增大而減小,隨激光強度的增強和離子質(zhì)量的增大而上升,而與靶厚無明顯關(guān)聯(lián).離子截止能量的映射結(jié)果與離子峰值能量的映射結(jié)果呈現(xiàn)出一致的單調(diào)性變化規(guī)律.因此,將結(jié)合 離子峰值能量Ep與離子截止能量Em的比值映射模型給出HB-RPA 機制下離子截止能量的經(jīng)驗公式.

    4.3.3 能量比值關(guān)系

    將每組模擬截止能量Em和峰值能量Ep的比值記為k=Em/Ep.將k作為建模參量,神經(jīng)網(wǎng)絡結(jié)構(gòu)與超參數(shù)設置均與訓練Ep和Em時保持一致,損失函數(shù)閾值設為2×10-4,重復與Ep,Em類似的訓練和代理模型加權(quán)平均操作,得到k關(guān)于不同輸入?yún)⒘康挠成淠P?

    圖8 為質(zhì)子、He2+、C6+、O8+四種離子的截止能量Em與峰值能量Ep的比值k的二維連續(xù)映射圖.由圖8(a) 所示,質(zhì)子的能量比值k隨靶密度n增大而增大,隨激光強度a增大而減小,k值在較小范圍內(nèi)緩慢變化且近似均勻.在a<90 ,n>25nc的部分區(qū)域,k變化幅度增大,呈現(xiàn)出一定非線性變化趨勢.圖8(b)–(d)顯示,其他3 種離子的k值變化趨勢與質(zhì)子相同,且k值的變化范圍隨離子質(zhì)量數(shù)A的增大而增大.

    圖8 (a) 質(zhì)子,(b) He2+,(c) C6+和 (d) O8+的能量比值 k=Em/Ep 的二維連續(xù)映射圖,其中靶厚 d=20λ0,參 數(shù)映射范圍 為63.2 ≤a ≤173.2 和10nc ≤n ≤30ncFig.8.Two-dimensional continuous mapping of k=Em/Ep over 63.2 ≤a ≤173.2 and 10nc ≤n ≤30nc with d=20λ0 :(a) Proton case;(b) He2+;(c) C6+;(d) O8+.

    圖9 為能量比值k隨激光強度和靶密度變化的映射曲線.圖9(a)顯示,k隨著激光強度a的增大近似以線性規(guī)律減小,其斜率隨靶密度增大而增大.圖9(b)顯示,當n<25nc時,k隨著靶密度n的增大而近似線性增大;當靶密度n>25nc時,k的上升幅度顯著變大,呈現(xiàn)非線性變化趨勢.前面圖1 中關(guān)于離子相空間的分析已給出: 峰值能量對應的離子主要在電荷分離場的峰值區(qū)域內(nèi)被加速;而截止能量對應的離子先被電荷分離場的上升沿預加速,再進入峰值區(qū)域被加速,所以具有更高的能量.當激光強度較大或靶密度較低時,電荷分離場峰值強度增加,導致離子峰值能量增大,截止能量與峰值能量差距減小.當激光強度較小或靶密度較高時,電荷分離場峰值強度降低,離子受電荷分離場上升沿加速所獲得的能量占比增大,導致截止能量與峰值能量差距增大.綜合上述映射模型特征和物理過程分析,給出了能量比值k關(guān)于輸入?yún)⒘康臄M合公式:

    圖9 C6+離子能量比值 k 隨(a)激光強度和(b)靶密度的映射曲線(實線)和擬合公式結(jié)果(虛線),其中參數(shù)范圍分別為63.2 ≤a ≤173.2,10nc ≤n ≤30nc ,靶厚d=15λ0Fig.9.Parameter scan (solid) and fitted formula (dashed) of k over (a) 63.2 ≤a ≤173.2 and (b) 10nc ≤n ≤30nc for C6+ with d=15λ0.

    在n<25nc時,(5)式中第1 項描述了k關(guān)于激光強度a的線性變化趨勢,以及靶密度n和離子質(zhì)量數(shù)A對激光強度映射的影響;第2 項和第3 項反映了靶密度n和離子質(zhì)量數(shù)A對k值整體取值分布的影響.n>25nc時,第1 項與第2 項中靶密度n增加了高次項,描述了在較大密度下能量比值關(guān)于n變化的非線性趨勢,此時上升幅度增大.離子截止能量Em可以由下式給出:

    5 結(jié)論與展望

    本文基于FCNN 算法對粒子模擬所生成的數(shù)據(jù)集展開了分析,建立了HB-RPA 機制下離子峰值能量和截止能量關(guān)于激光強度、靶密度、靶厚度和離子質(zhì)量的連續(xù)映射模型.該模型在與測試集數(shù)據(jù)和理論公式的比對中表現(xiàn)良好,且具有一定參數(shù)范圍內(nèi)的預測能力.研究發(fā)現(xiàn),當離子質(zhì)量數(shù)較大時,HB-RPA 機制下的峰值能量公式與數(shù)值模擬結(jié)果存在差異.基于映射模型和數(shù)值模擬結(jié)果,給出了離子峰值能量的修正公式和離子截止能量的經(jīng)驗公式.在基于HB-RPA 的激光加速離子實驗設計中,我們的映射模型和經(jīng)驗公式可為激光和靶參數(shù)的選取提供重要參考.本工作也可以拓展到更多樣化的目標參量,對離子束產(chǎn)額、單能性和準直度等束流品質(zhì)進行分析.

    猜你喜歡
    模型
    一半模型
    一種去中心化的域名服務本地化模型
    適用于BDS-3 PPP的隨機模型
    提煉模型 突破難點
    函數(shù)模型及應用
    p150Glued在帕金森病模型中的表達及分布
    函數(shù)模型及應用
    重要模型『一線三等角』
    重尾非線性自回歸模型自加權(quán)M-估計的漸近分布
    3D打印中的模型分割與打包
    长腿黑丝高跟| 丝袜人妻中文字幕| 中文字幕另类日韩欧美亚洲嫩草| 国产精品av久久久久免费| av在线天堂中文字幕 | 国产av又大| 久久久久国产一级毛片高清牌| 精品国产乱子伦一区二区三区| 激情在线观看视频在线高清| 午夜福利欧美成人| 日本a在线网址| 热99re8久久精品国产| 91大片在线观看| 国产精品久久久久成人av| 麻豆久久精品国产亚洲av | 久久久久九九精品影院| 久久天躁狠狠躁夜夜2o2o| 国产色视频综合| 国产成人av教育| 成年女人毛片免费观看观看9| 色播在线永久视频| 亚洲成人久久性| 美女高潮喷水抽搐中文字幕| 亚洲av第一区精品v没综合| 丝袜美腿诱惑在线| 亚洲精品一二三| 亚洲,欧美精品.| 久久精品影院6| 日韩视频一区二区在线观看| 免费人成视频x8x8入口观看| 91国产中文字幕| 高清av免费在线| 精品国产乱子伦一区二区三区| 中文亚洲av片在线观看爽| 久久精品aⅴ一区二区三区四区| 国产成人精品久久二区二区91| 极品人妻少妇av视频| 成人精品一区二区免费| 大型黄色视频在线免费观看| 精品久久久久久久毛片微露脸| 免费在线观看黄色视频的| 亚洲精品久久午夜乱码| 国产深夜福利视频在线观看| 亚洲九九香蕉| 欧美丝袜亚洲另类 | 侵犯人妻中文字幕一二三四区| 国产乱人伦免费视频| 国产99白浆流出| 国产野战对白在线观看| 美女高潮喷水抽搐中文字幕| 久99久视频精品免费| 日日爽夜夜爽网站| 日韩免费av在线播放| 亚洲精品国产精品久久久不卡| 麻豆国产av国片精品| 国产成人av激情在线播放| 动漫黄色视频在线观看| 19禁男女啪啪无遮挡网站| 欧美日韩视频精品一区| 麻豆一二三区av精品| 亚洲全国av大片| 国产午夜精品久久久久久| 久9热在线精品视频| 黄色怎么调成土黄色| 亚洲精品在线美女| 色综合站精品国产| 精品无人区乱码1区二区| 欧美中文综合在线视频| 国产亚洲欧美精品永久| av天堂久久9| 亚洲中文日韩欧美视频| 美女国产高潮福利片在线看| 亚洲专区中文字幕在线| 一进一出好大好爽视频| 亚洲精品粉嫩美女一区| 国产精品爽爽va在线观看网站 | 好看av亚洲va欧美ⅴa在| 亚洲精品一卡2卡三卡4卡5卡| 波多野结衣高清无吗| 日韩欧美一区二区三区在线观看| 国产亚洲欧美在线一区二区| 91九色精品人成在线观看| 99在线人妻在线中文字幕| 国产伦人伦偷精品视频| 在线观看免费午夜福利视频| e午夜精品久久久久久久| 最近最新中文字幕大全电影3 | 18禁黄网站禁片午夜丰满| 三级毛片av免费| 69精品国产乱码久久久| 女人被狂操c到高潮| 国内久久婷婷六月综合欲色啪| 国产av一区在线观看免费| 香蕉国产在线看| 日本精品一区二区三区蜜桃| 亚洲成人精品中文字幕电影 | av欧美777| 欧美黑人欧美精品刺激| 国产亚洲精品久久久久5区| 又大又爽又粗| 免费不卡黄色视频| 日韩欧美在线二视频| 91老司机精品| 老司机深夜福利视频在线观看| 国产成+人综合+亚洲专区| 亚洲人成网站在线播放欧美日韩| 国产av在哪里看| 亚洲视频免费观看视频| 99国产综合亚洲精品| 久久精品成人免费网站| 狠狠狠狠99中文字幕| 又大又爽又粗| 视频区图区小说| 精品一区二区三区四区五区乱码| 欧美一级毛片孕妇| 人妻丰满熟妇av一区二区三区| 又黄又爽又免费观看的视频| 丁香欧美五月| av有码第一页| 国产精品av久久久久免费| 日本免费a在线| 日韩高清综合在线| av片东京热男人的天堂| 一级片'在线观看视频| 一级毛片高清免费大全| 午夜91福利影院| 色播在线永久视频| 久久香蕉国产精品| 成在线人永久免费视频| 国产精品偷伦视频观看了| 可以免费在线观看a视频的电影网站| 丰满迷人的少妇在线观看| 亚洲欧美一区二区三区久久| 女人被躁到高潮嗷嗷叫费观| 亚洲成人国产一区在线观看| 五月开心婷婷网| 亚洲aⅴ乱码一区二区在线播放 | www.999成人在线观看| 午夜亚洲福利在线播放| 中文字幕人妻熟女乱码| 亚洲精品国产一区二区精华液| 精品国产国语对白av| 亚洲 欧美 日韩 在线 免费| 亚洲少妇的诱惑av| 成年人免费黄色播放视频| 电影成人av| 久久久久国产一级毛片高清牌| 看黄色毛片网站| 大型黄色视频在线免费观看| 老鸭窝网址在线观看| 久久影院123| 伊人久久大香线蕉亚洲五| 亚洲精品美女久久av网站| 免费观看人在逋| 亚洲九九香蕉| 桃红色精品国产亚洲av| 精品午夜福利视频在线观看一区| 高潮久久久久久久久久久不卡| 美国免费a级毛片| 三上悠亚av全集在线观看| av在线播放免费不卡| 又大又爽又粗| 欧美乱色亚洲激情| 久久精品影院6| 精品国产一区二区久久| 免费观看精品视频网站| 岛国在线观看网站| 久久草成人影院| 亚洲va日本ⅴa欧美va伊人久久| 成人三级黄色视频| 在线看a的网站| 亚洲avbb在线观看| 亚洲午夜理论影院| 天堂俺去俺来也www色官网| 免费av毛片视频| a级毛片在线看网站| 欧美日韩福利视频一区二区| 精品国产国语对白av| 日本五十路高清| 欧美成人性av电影在线观看| 老司机午夜福利在线观看视频| 亚洲一区中文字幕在线| 国产亚洲精品第一综合不卡| 黄色视频,在线免费观看| 成熟少妇高潮喷水视频| 麻豆成人av在线观看| 亚洲激情在线av| 午夜精品在线福利| 中文字幕av电影在线播放| 国产精品一区二区精品视频观看| 无遮挡黄片免费观看| 国产亚洲精品久久久久5区| 熟女少妇亚洲综合色aaa.| 亚洲va日本ⅴa欧美va伊人久久| 日韩成人在线观看一区二区三区| 丝袜美腿诱惑在线| 久久久精品国产亚洲av高清涩受| 国产成人av激情在线播放| 曰老女人黄片| 99国产极品粉嫩在线观看| 中文字幕人妻丝袜一区二区| 国产成人影院久久av| 黑人猛操日本美女一级片| 日韩国内少妇激情av| 亚洲精品一二三| 色婷婷av一区二区三区视频| 超碰97精品在线观看| 久久人妻福利社区极品人妻图片| 又黄又粗又硬又大视频| 婷婷六月久久综合丁香| 亚洲五月婷婷丁香| 最近最新中文字幕大全电影3 | 亚洲欧美日韩无卡精品| 日韩有码中文字幕| 亚洲精品一卡2卡三卡4卡5卡| 亚洲中文av在线| av片东京热男人的天堂| 我的亚洲天堂| 久9热在线精品视频| www国产在线视频色| 啦啦啦在线免费观看视频4| 超碰97精品在线观看| 亚洲第一av免费看| 欧美黑人精品巨大| 天堂中文最新版在线下载| 中出人妻视频一区二区| 看免费av毛片| 日韩国内少妇激情av| 老熟妇仑乱视频hdxx| 亚洲av成人av| 青草久久国产| 久久国产精品人妻蜜桃| 亚洲欧美精品综合一区二区三区| 亚洲视频免费观看视频| 男人的好看免费观看在线视频 | av网站在线播放免费| 丰满饥渴人妻一区二区三| 国产亚洲精品久久久久5区| 精品国产乱子伦一区二区三区| 国产野战对白在线观看| 亚洲人成电影免费在线| 麻豆成人av在线观看| 国产有黄有色有爽视频| 看黄色毛片网站| 亚洲视频免费观看视频| 熟女少妇亚洲综合色aaa.| 伊人久久大香线蕉亚洲五| www.www免费av| 一本大道久久a久久精品| 精品一品国产午夜福利视频| 日韩大码丰满熟妇| 夫妻午夜视频| 麻豆成人av在线观看| 国产三级黄色录像| 国产国语露脸激情在线看| 伦理电影免费视频| 韩国av一区二区三区四区| 亚洲av五月六月丁香网| 色在线成人网| 天堂动漫精品| 如日韩欧美国产精品一区二区三区| 最近最新中文字幕大全电影3 | 欧美激情 高清一区二区三区| 久久久久久久久久久久大奶| 黑人猛操日本美女一级片| 99久久久亚洲精品蜜臀av| 久久人妻熟女aⅴ| 欧美中文综合在线视频| 日韩精品中文字幕看吧| 一级a爱片免费观看的视频| 男人操女人黄网站| 亚洲欧美日韩另类电影网站| 国产成人精品无人区| 午夜日韩欧美国产| 在线天堂中文资源库| 级片在线观看| 男女做爰动态图高潮gif福利片 | 黄色丝袜av网址大全| 91成人精品电影| 久9热在线精品视频| 亚洲专区国产一区二区| 国产精品久久久久久人妻精品电影| 精品卡一卡二卡四卡免费| 欧美黑人欧美精品刺激| 国内久久婷婷六月综合欲色啪| 亚洲欧美日韩无卡精品| av片东京热男人的天堂| 亚洲 欧美 日韩 在线 免费| av有码第一页| 人人妻人人澡人人看| 搡老熟女国产l中国老女人| av视频免费观看在线观看| 青草久久国产| 最新美女视频免费是黄的| 超碰成人久久| 天天躁狠狠躁夜夜躁狠狠躁| 久久精品人人爽人人爽视色| 一区在线观看完整版| 国产激情欧美一区二区| 涩涩av久久男人的天堂| 夜夜爽天天搞| 久久午夜综合久久蜜桃| 悠悠久久av| 午夜福利在线免费观看网站| 极品人妻少妇av视频| 手机成人av网站| 久久中文字幕人妻熟女| 在线观看舔阴道视频| 日韩三级视频一区二区三区| 91老司机精品| av在线播放免费不卡| 在线观看免费日韩欧美大片| www.熟女人妻精品国产| 黄网站色视频无遮挡免费观看| 日本黄色日本黄色录像| 性少妇av在线| 欧美乱色亚洲激情| av视频免费观看在线观看| 日本一区二区免费在线视频| 热re99久久国产66热| 美女午夜性视频免费| 桃色一区二区三区在线观看| 国产男靠女视频免费网站| bbb黄色大片| 免费在线观看影片大全网站| 这个男人来自地球电影免费观看| 午夜日韩欧美国产| 男女做爰动态图高潮gif福利片 | 欧美色视频一区免费| 一边摸一边做爽爽视频免费| 亚洲一区二区三区欧美精品| 欧美国产精品va在线观看不卡| 88av欧美| 欧美人与性动交α欧美精品济南到| 亚洲 欧美一区二区三区| ponron亚洲| 成人亚洲精品一区在线观看| 十分钟在线观看高清视频www| 国产精品国产高清国产av| av视频免费观看在线观看| 黑人欧美特级aaaaaa片| 国产高清激情床上av| 操出白浆在线播放| 高清欧美精品videossex| 深夜精品福利| 伦理电影免费视频| 亚洲成人国产一区在线观看| 9色porny在线观看| 一二三四社区在线视频社区8| 久久热在线av| 88av欧美| 欧美最黄视频在线播放免费 | 大型av网站在线播放| 日日摸夜夜添夜夜添小说| 亚洲精品久久成人aⅴ小说| 精品欧美一区二区三区在线| 精品国产一区二区三区四区第35| 亚洲精品国产一区二区精华液| 日韩 欧美 亚洲 中文字幕| 老司机午夜福利在线观看视频| 亚洲美女黄片视频| 国产成人欧美在线观看| 亚洲国产看品久久| 热99re8久久精品国产| 亚洲自偷自拍图片 自拍| 国产精品国产av在线观看| 大码成人一级视频| 欧美黑人精品巨大| 性少妇av在线| 欧美不卡视频在线免费观看 | 亚洲av电影在线进入| 午夜成年电影在线免费观看| 嫩草影视91久久| 午夜福利在线观看吧| 久久 成人 亚洲| 99久久人妻综合| 国产成人精品在线电影| 美女 人体艺术 gogo| 激情视频va一区二区三区| 十八禁人妻一区二区| 一边摸一边做爽爽视频免费| 精品国产乱码久久久久久男人| 淫妇啪啪啪对白视频| 老司机靠b影院| 亚洲国产精品合色在线| av在线天堂中文字幕 | 曰老女人黄片| av有码第一页| 国产精品一区二区免费欧美| 亚洲国产看品久久| 免费不卡黄色视频| 999精品在线视频| 久久精品91无色码中文字幕| 在线观看舔阴道视频| 黄色女人牲交| 久久这里只有精品19| ponron亚洲| 中文欧美无线码| 亚洲中文字幕日韩| 亚洲精华国产精华精| 亚洲午夜精品一区,二区,三区| 久久中文字幕人妻熟女| 国产精品野战在线观看 | xxxhd国产人妻xxx| 亚洲精品一二三| 法律面前人人平等表现在哪些方面| 欧美av亚洲av综合av国产av| 免费不卡黄色视频| 亚洲av成人不卡在线观看播放网| 高清毛片免费观看视频网站 | 一二三四社区在线视频社区8| 久久午夜综合久久蜜桃| 国产真人三级小视频在线观看| 午夜老司机福利片| 亚洲精品在线美女| 久久久久久久久久久久大奶| 欧美国产精品va在线观看不卡| 叶爱在线成人免费视频播放| 一级毛片高清免费大全| 久久久水蜜桃国产精品网| 亚洲精品国产色婷婷电影| 村上凉子中文字幕在线| 真人一进一出gif抽搐免费| 精品福利观看| 免费少妇av软件| 欧美午夜高清在线| 精品无人区乱码1区二区| 免费在线观看视频国产中文字幕亚洲| 黄色视频不卡| 欧美久久黑人一区二区| av欧美777| 亚洲一卡2卡3卡4卡5卡精品中文| 亚洲专区国产一区二区| 国产一区二区三区在线臀色熟女 | 无限看片的www在线观看| 男女下面进入的视频免费午夜 | 最新美女视频免费是黄的| 国产精品久久久av美女十八| 黄色片一级片一级黄色片| 国产视频一区二区在线看| 精品国产乱码久久久久久男人| 免费一级毛片在线播放高清视频 | 午夜福利一区二区在线看| bbb黄色大片| 亚洲欧美激情在线| 69av精品久久久久久| 脱女人内裤的视频| 亚洲国产看品久久| 国产精品永久免费网站| 一二三四社区在线视频社区8| 精品高清国产在线一区| 在线看a的网站| 国产精品免费视频内射| 制服诱惑二区| 神马国产精品三级电影在线观看 | 男女做爰动态图高潮gif福利片 | 两性午夜刺激爽爽歪歪视频在线观看 | 伊人久久大香线蕉亚洲五| 国产熟女午夜一区二区三区| 免费看a级黄色片| 欧美另类亚洲清纯唯美| 亚洲第一青青草原| 色综合婷婷激情| 狠狠狠狠99中文字幕| 国产精品亚洲一级av第二区| 久99久视频精品免费| 无遮挡黄片免费观看| 久久伊人香网站| 天堂动漫精品| 国产成+人综合+亚洲专区| 午夜亚洲福利在线播放| 搡老熟女国产l中国老女人| 亚洲av电影在线进入| 久久国产亚洲av麻豆专区| 午夜福利,免费看| 午夜成年电影在线免费观看| 国产真人三级小视频在线观看| 精品久久蜜臀av无| 中文亚洲av片在线观看爽| 1024视频免费在线观看| 三上悠亚av全集在线观看| 国产单亲对白刺激| 欧美激情 高清一区二区三区| 窝窝影院91人妻| 丝袜美腿诱惑在线| 在线观看免费视频日本深夜| 三级毛片av免费| 99re在线观看精品视频| 精品一区二区三区av网在线观看| 黄色视频,在线免费观看| 97人妻天天添夜夜摸| 久久久国产成人免费| 亚洲自拍偷在线| 男女高潮啪啪啪动态图| 国产一区二区在线av高清观看| 国产精品久久久久成人av| 1024视频免费在线观看| 嫩草影视91久久| 成熟少妇高潮喷水视频| 99国产精品免费福利视频| 搡老熟女国产l中国老女人| av超薄肉色丝袜交足视频| 女人被躁到高潮嗷嗷叫费观| 国产成人欧美在线观看| 亚洲av美国av| 91九色精品人成在线观看| 中亚洲国语对白在线视频| 丝袜在线中文字幕| 亚洲第一欧美日韩一区二区三区| 欧美黑人欧美精品刺激| 一本大道久久a久久精品| 日韩高清综合在线| 精品人妻1区二区| 国产成人啪精品午夜网站| 国产一区二区三区综合在线观看| 久久精品亚洲熟妇少妇任你| 国产精品成人在线| 波多野结衣av一区二区av| 成人18禁高潮啪啪吃奶动态图| 淫秽高清视频在线观看| 亚洲精品一卡2卡三卡4卡5卡| 免费在线观看黄色视频的| 国产不卡一卡二| 欧美最黄视频在线播放免费 | 丁香六月欧美| 亚洲一码二码三码区别大吗| 少妇粗大呻吟视频| 熟女少妇亚洲综合色aaa.| 国产一区二区三区综合在线观看| 在线十欧美十亚洲十日本专区| 亚洲五月婷婷丁香| 波多野结衣av一区二区av| 9色porny在线观看| 国产伦人伦偷精品视频| 身体一侧抽搐| 99久久综合精品五月天人人| 欧美日本中文国产一区发布| 又黄又粗又硬又大视频| 制服诱惑二区| 美女午夜性视频免费| 亚洲 欧美一区二区三区| 国产成人精品久久二区二区91| 丝袜在线中文字幕| 欧美日韩亚洲国产一区二区在线观看| 亚洲精品粉嫩美女一区| 成年版毛片免费区| 国产国语露脸激情在线看| 精品人妻在线不人妻| 看黄色毛片网站| 一边摸一边抽搐一进一小说| 亚洲全国av大片| 国产精品永久免费网站| 成在线人永久免费视频| 日本黄色视频三级网站网址| 成年女人毛片免费观看观看9| 亚洲精品美女久久久久99蜜臀| 国产精品久久久久成人av| 国产av又大| 丁香欧美五月| 中文字幕av电影在线播放| 电影成人av| 国产精品影院久久| 最新美女视频免费是黄的| 欧美日韩乱码在线| 悠悠久久av| www.精华液| 国产一区二区三区在线臀色熟女 | 亚洲欧美激情综合另类| 久久中文看片网| 欧美性长视频在线观看| av福利片在线| 少妇被粗大的猛进出69影院| 香蕉久久夜色| 国产av又大| 在线观看免费视频网站a站| 成人精品一区二区免费| 日韩大码丰满熟妇| 在线观看一区二区三区激情| 精品国产乱子伦一区二区三区| 99国产精品一区二区蜜桃av| 伊人久久大香线蕉亚洲五| 国产高清国产精品国产三级| 日本免费a在线| 国产高清视频在线播放一区| www.精华液| 亚洲性夜色夜夜综合| 妹子高潮喷水视频| 欧美日韩亚洲国产一区二区在线观看| videosex国产| 久久精品91蜜桃| 欧美日韩亚洲国产一区二区在线观看| 久久久久久大精品| 久久精品91蜜桃| 欧美另类亚洲清纯唯美| 亚洲狠狠婷婷综合久久图片| 久久久久久免费高清国产稀缺| 高清av免费在线| 亚洲狠狠婷婷综合久久图片| 国产无遮挡羞羞视频在线观看| 久久久国产一区二区| 精品久久久久久久毛片微露脸| 动漫黄色视频在线观看| 熟女少妇亚洲综合色aaa.| 亚洲人成77777在线视频| 精品国产超薄肉色丝袜足j| 美女扒开内裤让男人捅视频| 午夜福利免费观看在线| 日本黄色视频三级网站网址| 高清av免费在线| 99精国产麻豆久久婷婷| 999精品在线视频| 国内久久婷婷六月综合欲色啪| 久久久国产精品麻豆| 在线十欧美十亚洲十日本专区| 中国美女看黄片|