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

    基于MDSC程序的Z箍縮內(nèi)爆單溫和三溫模擬分析*

    2016-04-25 08:46:33王剛?cè)A闞明先
    高壓物理學(xué)報(bào) 2016年1期
    關(guān)鍵詞:狀態(tài)方程套筒等離子體

    楊 龍,王剛?cè)A,闞明先,李 平

    (中國(guó)工程物理研究院流體物理研究所,四川綿陽(yáng) 621999)

    1 引 言

    強(qiáng)電流脈沖驅(qū)動(dòng)的Z箍縮裝置在過(guò)去的幾十年里一直作為一種高功率的X射線裝置被人們研究和使用。特別是在20世紀(jì)90年代以后,脈沖功率技術(shù)的進(jìn)步使Z箍縮科學(xué)得到了迅速發(fā)展,美國(guó)、俄羅斯等國(guó)家都建立起大型的Z箍縮裝置,希望在實(shí)驗(yàn)室內(nèi)創(chuàng)造一個(gè)強(qiáng)X射線輻射源,用以研究原子參數(shù)、輻射對(duì)材料的損傷、輻射輸運(yùn)、高溫高壓物態(tài)方程和不透明度、內(nèi)爆動(dòng)力學(xué)和界面不穩(wěn)定性等物理問(wèn)題。受實(shí)驗(yàn)手段的限制,諸如輻射輸運(yùn)、界面不穩(wěn)定性的研究等需要借助模擬工具,為此國(guó)內(nèi)外許多科研機(jī)構(gòu)各自發(fā)展了二維和三維輻射磁流體力學(xué)程序[1-7],對(duì)理論和實(shí)驗(yàn)研究都起到了很大的作用。進(jìn)入21世紀(jì)以來(lái),Z箍縮科學(xué)的發(fā)展方向轉(zhuǎn)向大質(zhì)量的套筒負(fù)載,不僅可以用于研究物態(tài)方程,還可以用來(lái)形成極端的高溫高壓環(huán)境,從而進(jìn)行磁化套筒慣性約束聚變(Magnetized Liner Inertial Fusion,MagLIF)研究[8-9]。然而,這些實(shí)驗(yàn)研究也需要借助于模擬程序來(lái)進(jìn)行負(fù)載參數(shù)的選取。模擬程序中的輻射磁流體模型有單溫模型和三溫模型,但是到目前為止,還沒有對(duì)這兩種不同的模型進(jìn)行對(duì)比研究。單溫模型和三溫模型之間的差異有多大,為什么會(huì)有這樣的差異等問(wèn)題不得而知。了解兩種模型之間的差異,有助于更好地運(yùn)用輻射磁流體模型。

    本研究通過(guò)改變負(fù)載材料、初始溫度和負(fù)載質(zhì)量,代表性地選取3個(gè)算例,運(yùn)用單溫模型和三溫模型分別進(jìn)行模擬計(jì)算,對(duì)比計(jì)算結(jié)果,分析產(chǎn)生差異的原因。

    2 MDSC程序及物理模型

    2.1 MDSC程序簡(jiǎn)介

    本研究所用的模擬程序?yàn)橹袊?guó)工程物理研究院流體物理研究所開發(fā)的MDSC磁流體程序[4,7]。該程序采用合理優(yōu)化的金屬電阻率模型[10],考慮了熱擴(kuò)散、焦耳熱、磁擴(kuò)散等物理過(guò)程,是二維、單流體的輻射磁流體力學(xué)模擬程序,目前已成功應(yīng)用于Z箍縮絲陣和套筒內(nèi)爆研究[11]以及磁驅(qū)動(dòng)飛片研究[12-13]。此前的MDSC輻射磁流體程序是單溫程序,最近又研發(fā)出基于三溫模型的三溫程序。

    2.2 物理模型

    采用單流體模型,考慮一個(gè)有限長(zhǎng)度的空心圓柱導(dǎo)體在大電流脈沖作用下的融化、汽化、等離子體化,最后在磁場(chǎng)作用下自箍縮的過(guò)程[14]。在負(fù)載導(dǎo)體從固態(tài)發(fā)展到等離子體態(tài)的過(guò)程中,其物質(zhì)成分是非常復(fù)雜的,有原子、離子、電子、光子等,只有當(dāng)它們內(nèi)部或彼此之間達(dá)到了局域熱平衡,溫度才有意義。物質(zhì)由初始的非平衡態(tài)到達(dá)局域熱平衡態(tài)需要一定的能量弛豫時(shí)間,記電子之間、離子之間、電子與離子之間的弛豫時(shí)間分別為τee、τii和τei,則3個(gè)弛豫時(shí)間量級(jí)的比為[15]

    式中:mi、me分別為離子和電子的質(zhì)量。如果等離子體系統(tǒng)的時(shí)間tp滿足τeeτei,則電子與離子之間達(dá)到局域熱平衡狀態(tài),所有物質(zhì)可以用同一個(gè)溫度來(lái)表示,單溫假設(shè)來(lái)源于此。假設(shè)只有一個(gè)溫度的輻射磁流體模型即為單溫模型,其物態(tài)方程取物質(zhì)的原子狀態(tài)方程[4]。如果認(rèn)為離子、電子和光子的溫度不相等,在局域熱平衡下有各自的溫度和壓力,那么在模型的建立過(guò)程中就需要3個(gè)不同的能量方程來(lái)描述離子、電子和光子本身的熱擴(kuò)散過(guò)程以及這些粒子間的熱傳導(dǎo)和熱輻射過(guò)程,此時(shí)的模型就是三溫模型,分別采用離子、電子和光子的狀態(tài)方程[4-5]。

    單溫模型的能量方程為

    三溫模型的能量方程為

    將三溫能量方程(即(3)式)代替單溫能量方程(即(2)式)以及用離子、電子和光子的狀態(tài)方程代替單溫原子狀態(tài)方程是三溫模型和單溫模型的主要差別。

    3 程序檢驗(yàn)

    為了檢驗(yàn)單溫(1T)模型和三溫(3T)模型對(duì)模擬結(jié)果的影響,選取3個(gè)計(jì)算算例[16],對(duì)比單溫程序和三溫程序得到的內(nèi)爆路徑,進(jìn)而分析兩種模型的合理性。計(jì)算程序中,r方向和z方向的網(wǎng)格尺寸均為30 μm,初始時(shí)間步長(zhǎng)為1 ps,時(shí)間步長(zhǎng)隨著計(jì)算需求而變化,狀態(tài)方程庫(kù)采用Sesame數(shù)據(jù)庫(kù)[17],輸入電流波形如圖1所示。表1列出了3個(gè)算例的模型參數(shù),其中mtot、H、D和δ分別為套筒的總質(zhì)量、高度、外徑和厚度,ρ0為初始密度,T0為初始溫度。

    表1 3個(gè)計(jì)算模型的初始參數(shù)Table 1 Initial parameters for three calculation models

    圖1 采用單溫程序和三溫程序計(jì)算得到的3個(gè)計(jì)算模型的自由面運(yùn)動(dòng)軌跡Fig.1 Trajectories of the inner surface calculated by mono-temperature and tri-temperature programs in the three models

    計(jì)算得到的套筒自由面(即套筒內(nèi)表面)的運(yùn)動(dòng)軌跡如圖1所示。模擬計(jì)算結(jié)果列于表2,其中:timp為碰軸時(shí)間,其實(shí)驗(yàn)值取最大輻射功率所對(duì)應(yīng)的時(shí)間,模擬值取自由面運(yùn)動(dòng)到距軸0.1 mm時(shí)的時(shí)間;vmax為套筒自由面的最大速度。從表2可以看出:在3個(gè)計(jì)算模型中,三溫程序得到的內(nèi)爆速度(即vmax)均大于單溫程序得到的內(nèi)爆速度,而碰軸時(shí)間均小于單溫程序得到的結(jié)果;在具有大負(fù)載質(zhì)量的Model 3中,單溫程序和三溫程序所得的碰軸時(shí)間差達(dá)到30 ns,而在Model 1和Model 2中實(shí)驗(yàn)值與模擬計(jì)算值之間存在一定的差別,但是差別較小,表明單溫模型和三溫模型對(duì)計(jì)算結(jié)果產(chǎn)生了影響。

    表2 3個(gè)計(jì)算模型的MDSC程序模擬結(jié)果Table 2 MDSC calculation results of the three models

    4 結(jié)果分析

    4.1 對(duì)溫度的影響

    圖2 由單溫程序和三溫程序得到的Model 3的計(jì)算結(jié)果Fig.2 Calculation results of Model 3 by mono-temperature and tri-temperature programs

    4.2 對(duì)速度的影響

    除此之外,由于狀態(tài)方程的差異,單溫模型下加載面上的壓縮波不會(huì)在物質(zhì)內(nèi)部產(chǎn)生沖擊波,而三溫模型下物質(zhì)內(nèi)部產(chǎn)生了沖擊波。從Model 3中300 ns時(shí)負(fù)載區(qū)的壓力剖面(見圖2(d))可以明顯地看到:三溫模型所得的壓縮波剖面陡峭,已經(jīng)形成了沖擊波;而在單溫模型下的相同時(shí)刻,波剖面變化平緩。當(dāng)沖擊波到達(dá)自由面時(shí),在自由面上產(chǎn)生一層稀薄的先驅(qū)等離子體,該先驅(qū)等離子體的速度很快,促使套筒提前碰軸[11]。

    對(duì)于Model 1和Model 2,由于其初始溫度均為104K,此時(shí)負(fù)載已經(jīng)處于等離子體態(tài),膨脹過(guò)程不明顯,但是從圖3(a)和圖3(b)可以看出,三溫模型仍然保持著比單溫模型更快的膨脹趨勢(shì)。這是由于:(1) 三溫模型中的洛倫茲力做功比單溫模型大,只是不明顯;(2) 先驅(qū)等離子體使三溫模型下的等離子體提前碰軸,由同一時(shí)刻負(fù)載區(qū)的壓力剖面(見圖3(c)和圖3(d))可以看出,負(fù)載區(qū)內(nèi)部形成了沖擊波,沖擊波到達(dá)自由面后形成了稀薄的先驅(qū)等離子體,在圖3(a)和圖3(b)中則表現(xiàn)為58和62 ns時(shí)自由面速度發(fā)生突然轉(zhuǎn)折。

    圖3 由單溫程序和三溫程序得到的Model 1和Model 2的計(jì)算結(jié)果Fig.3 Calculation results of Model 1 and Model 2 by mono-temperature and tri-temperature programs

    4.3 小結(jié)

    從加載面的運(yùn)動(dòng)軌跡來(lái)看,Model 1和Model 2中由單溫模型和三溫模型所得的計(jì)算結(jié)果幾乎是重合的,而Model 3中兩個(gè)模型的計(jì)算結(jié)果相差很大。自由面的運(yùn)動(dòng)軌跡與加載面的運(yùn)動(dòng)軌跡和先驅(qū)等離子體相關(guān)。加載面的運(yùn)動(dòng)軌跡本質(zhì)上反映了洛倫茲力做功,在單溫模型和三溫模型得到的加載面運(yùn)動(dòng)軌跡相差很小時(shí),自由面運(yùn)動(dòng)軌跡也相差很小。在Model 1和Model 2中,單溫模型和三溫模型所得的碰軸時(shí)間分別相差3和1 ns;而在Model 3中由于單溫模型和三溫模型下加載面運(yùn)動(dòng)軌跡的差異顯著,因此碰軸時(shí)間相差較大,達(dá)到30 ns。雖然沖擊波也會(huì)影響自由面的運(yùn)動(dòng)軌跡,但是通過(guò)Model 1和Model 2可以看到,沖擊波的影響不大,因此在單溫模型和三溫模型中洛倫茲力做功的差異才是使碰軸時(shí)間產(chǎn)生差異的最主要原因。

    加載面在一定溫度和壓力下是否膨脹取決于其狀態(tài)方程。根據(jù)離子和原子的狀態(tài)方程,在低溫情況下,離子比原子更易于膨脹和壓縮,從而導(dǎo)致單溫模型和三溫模型計(jì)算得到的負(fù)載加載面運(yùn)動(dòng)軌跡的差異較大。在高溫情況下,由于壓力變小,加載面的運(yùn)動(dòng)軌跡反而趨于一致。研究表明,采用MDSC的單溫程序模擬磁驅(qū)動(dòng)飛片時(shí),所獲得的自由面速度歷史與實(shí)驗(yàn)結(jié)果符合很好[12-13]。因此,可以認(rèn)為:在低溫情況下,單溫模型更趨近于實(shí)際情況,此時(shí)大多數(shù)負(fù)載物質(zhì)并沒有等離子體化,需要考慮其結(jié)合力,這在單溫模型的狀態(tài)方程中已有所考慮,而三溫模型的離子狀態(tài)方程則沒有這種相互作用。

    5 結(jié) 論

    單溫模型和三溫模型的重要區(qū)別在于狀態(tài)方程。低溫下,兩種模型中狀態(tài)方程的區(qū)別很大,導(dǎo)致計(jì)算得到的加載面運(yùn)動(dòng)軌跡產(chǎn)生較大的差異,進(jìn)而使洛倫茲力做功產(chǎn)生差別,直接影響所獲得的動(dòng)能,是造成兩種模型所得碰軸時(shí)間相差很大的最主要原因;高溫下,兩種模型的狀態(tài)方程趨于一致,致使等離子體的膨脹和壓縮過(guò)程以及洛倫茲力做功均趨于一致。此外,三溫模型中更容易形成沖擊波,形成先驅(qū)等離子體,對(duì)于計(jì)算結(jié)果也將產(chǎn)生一定的影響。單溫模型在計(jì)算磁驅(qū)動(dòng)飛片、固體套筒等應(yīng)用中比三溫模型更符合實(shí)驗(yàn)數(shù)據(jù),采用三溫模型模擬固體物質(zhì)時(shí)應(yīng)考慮物質(zhì)變?yōu)榈入x子體的元過(guò)程。

    [1] PETERKIN R E Jr,DEGNAN J H,HUSSEY T W,et al.A long conduction time compact torus plasma opening switch [J].IEEE Trans Plasma Sci,1993,21(5):522-528.

    [2] HELLINGER P,MANGENEY A,MATTHEWS A.Whistler waves in 3D hybrid simulations of quasiperpendicular shocks [J].Geophys Res Lett,1996,23(6):621-624.

    [3] CHITTENDEN J P,LEBEDEV S V,JENNINGS C A,et al.X-ray generation mechanisms in three-dimensional simulations of wire arrayZ-pinches [J].Plasma Phys Control Fusion,2004,46(12B):457-476.

    [4] 闞明先,胡熙靜,王剛?cè)A.二維磁流體力學(xué)ALE數(shù)值模擬 [J].高能量密度物理,2009(2):56-60.

    KAN M X,HU X J,WANG G H.Two dimensional MHD ALE numerical simulation [J].High Energy Density Physics,2009(2):56-60.

    [5] 丁 寧,鄔吉明,戴自換,等.Z箍縮內(nèi)爆的MARED程序數(shù)值模擬分析 [J].物理學(xué)報(bào),2010,59(12):8707-8716.

    DING N,WU J M,DAI Z H,et al.Numerical simulation analysis ofZ-pinch implosion using MARED code [J].Acta Physica Sinica,2010,59(12):8707-8716.

    [6] ROBINSON A C,GARASI C J.Three-dimensionalz-pinch wire array modeling with ALEGRA-HEDP [J].Comput Phys Commun,2004,164(1/2/3):408-413.

    [7] 闞明先,蔣吉昊,王剛?cè)A,等.襯套內(nèi)爆ALE方法二維MHD數(shù)值模擬 [J].四川大學(xué)學(xué)報(bào)(自然科學(xué)版),2007,44(1):91-96.

    KAN M X,JIANG J H,WANG G H,et al.ALE simulation of 2D MHD for liner [J].Journal of Sichuan University (Natural Science Edition),2007,44(1):91-96.

    [8] MARTIN M R,LEMKE R W,MCBRIDE R D,et al.Solid liner implosions onZfor producing multi-megabar,shockless compressions [J].Phys Plasmas,2012,19(5):056310.

    [9] MCBRIDE R D,MARTIN M R,LEMKE R W,et al.Beryllium liner implosion experiments on theZaccelerator in preparation for magnetized liner inertial fusion [J].Phys Plasmas,2013,20(5):056309.

    [10] 闞明先,王剛?cè)A,趙海龍,等.金屬電阻率模型 [J].爆炸與沖擊,2013,33(3):282-286.

    KAN M X,WANG G H,ZHAO H L,et al.Electrical resistivity model for metals [J].Explosion and Shock Waves,2013,33(3):282-286.

    [11] 王剛?cè)A,胡熙靜,孫承緯.雙層鎢襯套Z箍縮內(nèi)爆數(shù)值模擬 [J].高壓物理學(xué)報(bào),2004,18(4):364-367.

    WANG G H,HU X J,SUN C W.Simulation of magnetohydrodynamics for plasma jetting on wire pinch [J].Chinese Journal of High Pressure Physics,2004,18(4):364-367.

    [12] 闞明先,王剛?cè)A,趙海龍,等.磁驅(qū)動(dòng)飛片二維磁流體力學(xué)數(shù)值模擬 [J].強(qiáng)激光與粒子束,2013,25(8):2137-2140.

    KAN M X,WANG G H,ZHAO H L,et al.Two-dimensional magneto-hydrodynamic simulations of magnetically accelerated flyer plates [J].High Power Laser and Particle Beams,2013,25(8):2137-2140.

    [13] 王剛?cè)A,孫承緯,趙劍衡,等.磁驅(qū)動(dòng)平面飛片的一維磁流體力學(xué)計(jì)算 [J].爆炸與沖擊,2008,28(3):261-264.

    WANG G H,SUN C W,ZHAO J H,et al.One-dimensional,magnetohydrodynamic simulations of magnetically driven flyer plates [J].Explosion and Shock Waves,2008,28(3):261-264.

    [14] HAINS M G.A heuristic model of the wire arrayZ-pinch [J].IEEE Trans Plasma Sci,1998,26(4):1275-1281.

    [15] 徐家鸞,金尚憲.等離子體物理學(xué) [M].北京:原子能出版社,1981:37-40.

    XU J L,JIN S X.Plasma physics [M].Beijing:Atomic Energy Press,1981:37-40.

    [16] SANDFORD T W L,ALLSHOUSE G O,MARDER B M,et al.Improved symmetry greatly increases X-ray power from wire-arrayZ-pinches [J].Phys Rev Lett,1996,77(25):5063-5066.

    [17] COOPER N G.An invitation to participate in the LASL equation of state library:LASL-79-62 [R].Los Alamos,USA:Los Alamos Scientific Laboratory,1979.

    [18] LYON S P,JOHNSON J D.T-1 handbook of the SESAME equation of state library:LA-CP-98-100 [R].Los Alamos,USA:Los Alamos Scientific Laboratory,1998.

    猜你喜歡
    狀態(tài)方程套筒等離子體
    套筒灌漿連接密實(shí)性檢測(cè)研究現(xiàn)狀及展望
    LKP狀態(tài)方程在天然氣熱物性參數(shù)計(jì)算的應(yīng)用
    煤氣與熱力(2021年6期)2021-07-28 07:21:30
    連續(xù)磁活動(dòng)對(duì)等離子體層演化的影響
    基于低溫等離子體修飾的PET/PVC浮選分離
    一種尾架套筒自動(dòng)機(jī)械鎖緊機(jī)構(gòu)
    基于隨機(jī)與區(qū)間分析的狀態(tài)方程不確定性比較
    套筒類零件內(nèi)孔精加工工藝分析
    等離子體種子處理技術(shù)介紹
    用狀態(tài)方程模擬氨基酸水溶液的熱力學(xué)性質(zhì)
    鋁帶、箔軋機(jī)上、卸套筒裝置的結(jié)構(gòu)分析
    禁无遮挡网站| 亚洲自偷自拍三级| 亚洲av日韩在线播放| 日韩在线高清观看一区二区三区| 中文字幕人妻熟人妻熟丝袜美| 欧美丝袜亚洲另类| 最近的中文字幕免费完整| 高清欧美精品videossex| 六月丁香七月| 日本爱情动作片www.在线观看| 亚洲熟妇中文字幕五十中出| 久久精品夜夜夜夜夜久久蜜豆| 亚洲欧美精品自产自拍| 国产亚洲精品av在线| 国产探花在线观看一区二区| 久久久久久久国产电影| 中文乱码字字幕精品一区二区三区 | 国产成人一区二区在线| 精品亚洲乱码少妇综合久久| 亚洲成人久久爱视频| videossex国产| 免费观看无遮挡的男女| 激情五月婷婷亚洲| 久久久久久久久大av| av免费观看日本| 五月伊人婷婷丁香| 日产精品乱码卡一卡2卡三| 成人亚洲欧美一区二区av| 成人高潮视频无遮挡免费网站| videos熟女内射| 十八禁国产超污无遮挡网站| 国产亚洲一区二区精品| 成人av在线播放网站| 美女cb高潮喷水在线观看| 亚洲精品影视一区二区三区av| av一本久久久久| 一级毛片电影观看| 九九久久精品国产亚洲av麻豆| 只有这里有精品99| 亚洲av日韩在线播放| 亚洲电影在线观看av| 精品午夜福利在线看| 天堂中文最新版在线下载 | 免费看a级黄色片| 黄色日韩在线| 好男人视频免费观看在线| 国产精品国产三级专区第一集| 国产有黄有色有爽视频| 欧美区成人在线视频| av天堂中文字幕网| 欧美日韩综合久久久久久| 白带黄色成豆腐渣| 国产精品99久久久久久久久| 日韩在线高清观看一区二区三区| 少妇的逼水好多| 观看免费一级毛片| 青春草视频在线免费观看| 国产高清三级在线| 午夜精品一区二区三区免费看| 国产高清三级在线| 天堂av国产一区二区熟女人妻| 99久久精品一区二区三区| 亚洲国产精品国产精品| 女的被弄到高潮叫床怎么办| 99热这里只有精品一区| 国产精品伦人一区二区| 国产精品伦人一区二区| 免费看日本二区| 国产亚洲一区二区精品| 欧美3d第一页| 亚洲久久久久久中文字幕| 欧美 日韩 精品 国产| 超碰97精品在线观看| 久久久久久久大尺度免费视频| 18禁动态无遮挡网站| 三级国产精品片| 久热久热在线精品观看| 亚洲精品456在线播放app| 女人被狂操c到高潮| 国产精品国产三级国产专区5o| 亚洲成人一二三区av| 中文字幕亚洲精品专区| 直男gayav资源| 日韩电影二区| 国产成人午夜福利电影在线观看| 91精品一卡2卡3卡4卡| 免费看不卡的av| 久久午夜福利片| 免费大片18禁| 久久久久久久午夜电影| 我要看日韩黄色一级片| 最近中文字幕高清免费大全6| 91在线精品国自产拍蜜月| 综合色丁香网| 女人十人毛片免费观看3o分钟| 日韩av在线大香蕉| 成人无遮挡网站| 亚洲精品影视一区二区三区av| 18禁裸乳无遮挡免费网站照片| 久久久久九九精品影院| 大香蕉97超碰在线| 久久精品夜色国产| 午夜福利成人在线免费观看| 日日啪夜夜爽| av在线老鸭窝| av专区在线播放| 欧美人与善性xxx| 日本免费在线观看一区| 欧美精品一区二区大全| 日日啪夜夜撸| 精品久久久久久久久亚洲| 亚洲av.av天堂| 国产一区二区亚洲精品在线观看| 欧美性感艳星| 欧美成人a在线观看| 一本久久精品| 少妇熟女欧美另类| 亚洲久久久久久中文字幕| 又黄又爽又刺激的免费视频.| 十八禁国产超污无遮挡网站| 国产探花在线观看一区二区| 亚洲人成网站在线观看播放| 欧美+日韩+精品| 黄色日韩在线| 亚洲国产高清在线一区二区三| 蜜臀久久99精品久久宅男| 久久精品国产自在天天线| 亚洲欧美成人综合另类久久久| 久久久精品94久久精品| 亚洲aⅴ乱码一区二区在线播放| 精品一区二区三区视频在线| 男人狂女人下面高潮的视频| 男女那种视频在线观看| 97超碰精品成人国产| 成人亚洲精品一区在线观看 | 日韩大片免费观看网站| 日韩三级伦理在线观看| 最近视频中文字幕2019在线8| 欧美日韩视频高清一区二区三区二| 亚洲四区av| 51国产日韩欧美| 亚洲精品456在线播放app| 久久久久久久久久久丰满| 国产成人a区在线观看| 亚洲,欧美,日韩| 天天躁夜夜躁狠狠久久av| 亚洲精品乱久久久久久| 日韩强制内射视频| 久久6这里有精品| 亚洲av免费高清在线观看| 精品一区二区三区视频在线| 国产精品久久久久久久久免| 别揉我奶头 嗯啊视频| 国产人妻一区二区三区在| 能在线免费观看的黄片| 欧美激情在线99| 五月天丁香电影| 日韩 亚洲 欧美在线| 亚洲av不卡在线观看| 精品久久久噜噜| 日韩国内少妇激情av| 一级毛片黄色毛片免费观看视频| 亚洲av国产av综合av卡| 亚洲国产av新网站| 精品午夜福利在线看| 日韩av免费高清视频| 波野结衣二区三区在线| 欧美极品一区二区三区四区| 乱码一卡2卡4卡精品| 尤物成人国产欧美一区二区三区| 亚洲av电影不卡..在线观看| 99热这里只有是精品在线观看| av线在线观看网站| 美女高潮的动态| 国产黄色视频一区二区在线观看| 最近手机中文字幕大全| 在现免费观看毛片| 亚洲最大成人手机在线| 精品久久久精品久久久| 亚洲精品第二区| 熟女人妻精品中文字幕| 午夜福利视频1000在线观看| 午夜免费激情av| 欧美成人a在线观看| 乱人视频在线观看| 性插视频无遮挡在线免费观看| 美女内射精品一级片tv| 狂野欧美激情性xxxx在线观看| 一级毛片久久久久久久久女| 卡戴珊不雅视频在线播放| 午夜激情欧美在线| 秋霞在线观看毛片| 久久精品夜色国产| 九九爱精品视频在线观看| 欧美成人精品欧美一级黄| 内地一区二区视频在线| 亚洲av在线观看美女高潮| 免费看av在线观看网站| 久久这里有精品视频免费| 在线观看人妻少妇| 高清在线视频一区二区三区| videos熟女内射| av在线老鸭窝| 99久国产av精品| 又爽又黄无遮挡网站| 99热这里只有是精品50| 熟妇人妻不卡中文字幕| 十八禁国产超污无遮挡网站| 午夜精品一区二区三区免费看| 亚洲av成人精品一区久久| 观看免费一级毛片| 99久久九九国产精品国产免费| 欧美人与善性xxx| 亚洲精品一二三| 亚洲三级黄色毛片| 午夜久久久久精精品| 大又大粗又爽又黄少妇毛片口| 国产成人a区在线观看| 国产成人a区在线观看| av专区在线播放| 亚洲人成网站高清观看| 日韩视频在线欧美| 欧美97在线视频| 国产日韩欧美在线精品| 国内精品宾馆在线| 好男人在线观看高清免费视频| 国产91av在线免费观看| av在线老鸭窝| 免费av观看视频| 插逼视频在线观看| 欧美日韩一区二区视频在线观看视频在线 | 夜夜爽夜夜爽视频| 中文字幕人妻熟人妻熟丝袜美| 秋霞伦理黄片| 成人高潮视频无遮挡免费网站| 午夜福利在线在线| 3wmmmm亚洲av在线观看| 国产精品一区二区三区四区久久| 国产成人精品婷婷| 久久综合国产亚洲精品| 久久97久久精品| 2021少妇久久久久久久久久久| 韩国av在线不卡| 久久精品夜色国产| 联通29元200g的流量卡| 欧美激情国产日韩精品一区| 欧美三级亚洲精品| freevideosex欧美| 日韩三级伦理在线观看| 九色成人免费人妻av| 国产精品女同一区二区软件| 成年免费大片在线观看| 国产欧美日韩精品一区二区| 看免费成人av毛片| 一级毛片黄色毛片免费观看视频| 美女大奶头视频| 日韩三级伦理在线观看| 免费大片18禁| 成年人午夜在线观看视频 | 国产久久久一区二区三区| 在线观看人妻少妇| 国产午夜精品一二区理论片| 亚洲在久久综合| 男人舔奶头视频| 日韩制服骚丝袜av| 久久久久久国产a免费观看| 天美传媒精品一区二区| 精品久久久久久久末码| 最近2019中文字幕mv第一页| 丝袜美腿在线中文| 成人无遮挡网站| 精品亚洲乱码少妇综合久久| 国产在线男女| 看十八女毛片水多多多| av在线观看视频网站免费| 成人特级av手机在线观看| 少妇猛男粗大的猛烈进出视频 | 一级毛片电影观看| 亚洲国产欧美人成| 亚洲婷婷狠狠爱综合网| 国产淫片久久久久久久久| 成人鲁丝片一二三区免费| 国产成人a区在线观看| 国产伦一二天堂av在线观看| 国产精品1区2区在线观看.| 欧美bdsm另类| 久久久久久国产a免费观看| 69av精品久久久久久| 国产乱来视频区| 午夜免费男女啪啪视频观看| 午夜免费激情av| 一级黄片播放器| 国产91av在线免费观看| 日韩欧美一区视频在线观看 | 日韩三级伦理在线观看| 一区二区三区四区激情视频| 成年人午夜在线观看视频 | 国产精品麻豆人妻色哟哟久久 | 在线a可以看的网站| 日韩av在线免费看完整版不卡| 97在线视频观看| 一本久久精品| 18禁裸乳无遮挡免费网站照片| 亚洲成人一二三区av| 美女被艹到高潮喷水动态| 你懂的网址亚洲精品在线观看| 草草在线视频免费看| 丰满乱子伦码专区| 嫩草影院入口| 特大巨黑吊av在线直播| 免费大片黄手机在线观看| 男女啪啪激烈高潮av片| 国产亚洲最大av| 久久久久久久大尺度免费视频| 啦啦啦韩国在线观看视频| 亚洲精品视频女| 夫妻午夜视频| 亚洲欧美日韩卡通动漫| 日日啪夜夜爽| 国产精品一区二区在线观看99 | 高清日韩中文字幕在线| 亚洲精品国产成人久久av| 国语对白做爰xxxⅹ性视频网站| 国产v大片淫在线免费观看| 日韩成人伦理影院| 免费观看av网站的网址| 免费观看a级毛片全部| 久久久久久国产a免费观看| 69av精品久久久久久| 精品亚洲乱码少妇综合久久| 日韩在线高清观看一区二区三区| 亚洲电影在线观看av| 婷婷色综合大香蕉| 在线观看免费高清a一片| 又爽又黄无遮挡网站| 啦啦啦韩国在线观看视频| 国产精品1区2区在线观看.| 亚洲伊人久久精品综合| 国产亚洲最大av| 高清欧美精品videossex| av线在线观看网站| 精品一区二区三区人妻视频| 日韩伦理黄色片| 一级二级三级毛片免费看| 一本一本综合久久| 如何舔出高潮| 日韩一区二区视频免费看| 国产老妇女一区| 日韩精品有码人妻一区| 国产国拍精品亚洲av在线观看| 亚洲国产精品成人久久小说| 最近最新中文字幕大全电影3| 免费看光身美女| 日本与韩国留学比较| 嫩草影院新地址| 少妇熟女欧美另类| 天堂网av新在线| 国产白丝娇喘喷水9色精品| 一级毛片 在线播放| 卡戴珊不雅视频在线播放| 狂野欧美白嫩少妇大欣赏| 国产高清三级在线| 日本-黄色视频高清免费观看| 国产一区二区三区av在线| 国产男人的电影天堂91| 插阴视频在线观看视频| 国产精品伦人一区二区| 亚洲成色77777| 97人妻精品一区二区三区麻豆| 久久久精品94久久精品| 精品久久久精品久久久| 久久久精品免费免费高清| 国产高潮美女av| 欧美精品国产亚洲| 在线免费观看不下载黄p国产| 成年女人在线观看亚洲视频 | 成年女人在线观看亚洲视频 | 99久久人妻综合| 日韩欧美 国产精品| 爱豆传媒免费全集在线观看| 一级毛片我不卡| 美女被艹到高潮喷水动态| 边亲边吃奶的免费视频| 免费av毛片视频| 一二三四中文在线观看免费高清| 国产在视频线在精品| 大陆偷拍与自拍| 在线观看一区二区三区| 毛片女人毛片| 日韩电影二区| 日韩av不卡免费在线播放| 久久精品久久精品一区二区三区| 久久久久性生活片| 国产成人午夜福利电影在线观看| 久久久午夜欧美精品| 午夜福利视频精品| 成年av动漫网址| 国产黄片美女视频| 中文精品一卡2卡3卡4更新| 久久国内精品自在自线图片| 国产亚洲av片在线观看秒播厂 | 三级男女做爰猛烈吃奶摸视频| 男女边吃奶边做爰视频| 国产精品熟女久久久久浪| 国产探花在线观看一区二区| 国产午夜精品一二区理论片| 日韩一本色道免费dvd| 国产精品无大码| 超碰av人人做人人爽久久| 一级黄片播放器| 亚洲av二区三区四区| 国产真实伦视频高清在线观看| 乱人视频在线观看| 极品少妇高潮喷水抽搐| 成人亚洲精品一区在线观看 | 日本色播在线视频| 少妇的逼好多水| 国产精品美女特级片免费视频播放器| 免费看光身美女| 国产不卡一卡二| 啦啦啦中文免费视频观看日本| 欧美潮喷喷水| 午夜福利高清视频| 欧美性感艳星| 国产精品麻豆人妻色哟哟久久 | 国内揄拍国产精品人妻在线| 欧美性感艳星| 久久久久久久亚洲中文字幕| 99视频精品全部免费 在线| av网站免费在线观看视频 | 在线观看免费高清a一片| 99热网站在线观看| kizo精华| 亚洲图色成人| 国产在线男女| 国产成人精品久久久久久| xxx大片免费视频| 婷婷色av中文字幕| 高清视频免费观看一区二区 | 精品99又大又爽又粗少妇毛片| 国产精品一区二区三区四区免费观看| 91午夜精品亚洲一区二区三区| 日韩av在线大香蕉| 午夜福利成人在线免费观看| 欧美激情国产日韩精品一区| 免费高清在线观看视频在线观看| 国产黄片视频在线免费观看| 国产高清三级在线| 午夜福利网站1000一区二区三区| 亚洲av成人av| 日韩亚洲欧美综合| 可以在线观看毛片的网站| 国产极品天堂在线| 欧美成人一区二区免费高清观看| 极品教师在线视频| 国产精品.久久久| 在线免费观看的www视频| 青春草国产在线视频| 91狼人影院| av国产免费在线观看| 免费看美女性在线毛片视频| 国产在线一区二区三区精| 亚洲成人av在线免费| 天堂√8在线中文| 中文字幕亚洲精品专区| 亚洲色图av天堂| 黑人高潮一二区| 亚洲av免费高清在线观看| 成人美女网站在线观看视频| 久久这里有精品视频免费| 国产激情偷乱视频一区二区| 国产伦一二天堂av在线观看| 国产真实伦视频高清在线观看| 一边亲一边摸免费视频| 春色校园在线视频观看| 国产精品麻豆人妻色哟哟久久 | 亚洲av电影在线观看一区二区三区 | 波野结衣二区三区在线| 丰满人妻一区二区三区视频av| 麻豆成人午夜福利视频| 国产伦理片在线播放av一区| 精品久久久久久久久av| 又爽又黄无遮挡网站| av.在线天堂| 一二三四中文在线观看免费高清| 观看免费一级毛片| 国产黄频视频在线观看| 久久国内精品自在自线图片| 国产在线男女| 五月伊人婷婷丁香| 久久韩国三级中文字幕| 免费av观看视频| 免费播放大片免费观看视频在线观看| 好男人在线观看高清免费视频| 免费av不卡在线播放| 午夜激情久久久久久久| a级毛色黄片| 我的老师免费观看完整版| 欧美性猛交╳xxx乱大交人| 99热这里只有精品一区| 精华霜和精华液先用哪个| 亚洲va在线va天堂va国产| 免费看不卡的av| 97在线视频观看| 91精品国产九色| 一级毛片黄色毛片免费观看视频| 国产免费视频播放在线视频 | 日韩,欧美,国产一区二区三区| 搡老乐熟女国产| 午夜福利成人在线免费观看| 午夜福利在线观看免费完整高清在| 亚洲欧美日韩卡通动漫| 亚洲在线观看片| 亚洲精品成人久久久久久| 人妻系列 视频| av福利片在线观看| 久久久精品欧美日韩精品| 国产淫语在线视频| 免费观看a级毛片全部| 国产精品嫩草影院av在线观看| 国产老妇女一区| 久久久精品欧美日韩精品| 国产色爽女视频免费观看| 久久久久久久亚洲中文字幕| 在线观看av片永久免费下载| 日本猛色少妇xxxxx猛交久久| 又爽又黄无遮挡网站| 国产av在哪里看| 亚洲高清免费不卡视频| 可以在线观看毛片的网站| 日日撸夜夜添| 久久精品久久久久久久性| 免费观看av网站的网址| 能在线免费看毛片的网站| 亚洲精华国产精华液的使用体验| 美女国产视频在线观看| 内地一区二区视频在线| 插逼视频在线观看| 舔av片在线| 国产大屁股一区二区在线视频| 边亲边吃奶的免费视频| 久久韩国三级中文字幕| 国产毛片a区久久久久| 国产 亚洲一区二区三区 | 色播亚洲综合网| 亚洲国产精品国产精品| 最新中文字幕久久久久| 人妻少妇偷人精品九色| av在线蜜桃| 91精品伊人久久大香线蕉| 国产成人午夜福利电影在线观看| 男女边吃奶边做爰视频| 亚洲人成网站在线播| 日产精品乱码卡一卡2卡三| 最近2019中文字幕mv第一页| 少妇高潮的动态图| 少妇的逼好多水| 国产淫语在线视频| 色综合亚洲欧美另类图片| 熟妇人妻不卡中文字幕| 男女边摸边吃奶| 欧美日韩国产mv在线观看视频 | 男人舔奶头视频| 亚洲aⅴ乱码一区二区在线播放| kizo精华| 听说在线观看完整版免费高清| 男女那种视频在线观看| 97人妻精品一区二区三区麻豆| 亚洲精品中文字幕在线视频 | av卡一久久| 亚洲天堂国产精品一区在线| 一级毛片 在线播放| 中文字幕免费在线视频6| 婷婷色综合大香蕉| 国产精品美女特级片免费视频播放器| 欧美日韩综合久久久久久| 国产在视频线精品| 亚洲欧美日韩东京热| 乱人视频在线观看| 久久久久久久久久成人| 色综合站精品国产| 欧美变态另类bdsm刘玥| 婷婷六月久久综合丁香| 一级av片app| 久久精品国产自在天天线| 在线免费十八禁| 白带黄色成豆腐渣| 欧美zozozo另类| 久久这里有精品视频免费| 亚洲国产av新网站| 十八禁国产超污无遮挡网站| av播播在线观看一区| 国产淫片久久久久久久久| 国产永久视频网站| 国产精品久久视频播放| 精品国产一区二区三区久久久樱花 | 亚洲一区高清亚洲精品| 日日摸夜夜添夜夜添av毛片| 欧美性猛交╳xxx乱大交人| 一级毛片aaaaaa免费看小| 黄片无遮挡物在线观看| 99久国产av精品国产电影| 91久久精品电影网| 精品久久久久久久末码| 久久久久免费精品人妻一区二区| 欧美一区二区亚洲| 啦啦啦韩国在线观看视频| 亚洲国产av新网站| 国产大屁股一区二区在线视频| 日本色播在线视频| 别揉我奶头 嗯啊视频| 国产真实伦视频高清在线观看| 久久久久性生活片| 日韩中字成人| 日韩一区二区三区影片|