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

    雙曲守恒律方程的Lax—Wendroff時間離散WENO格式

    2018-01-09 23:28:33李興華孫陽艾曉輝
    關(guān)鍵詞:高精度

    李興華+孫陽+艾曉輝

    摘 要:雙曲守恒型方程的高精度、高分辨率計算格式的研究一直是計算流體力學(xué)的熱點問題。針對原WENOJS格式分辨率較低和計算量偏大的不足問題,提出利用簡單的重構(gòu)數(shù)值通量的方法以提高計算效率,構(gòu)造了新的簡單限制器的5階迎風(fēng)型WENO格式。通過MATLAB軟件的仿真對LaxWendroff WENOJS格式、LaxWendroff簡單限制器WENO格式、RungeKutta WENOJS格式、RungeKutta簡單限制器的WENO格式的實驗結(jié)果進行了分析,并比較了這四種計算格式的計算效率和計算精度。數(shù)值實驗表明:新格式LaxWendroff簡單限制器WENO格式在保持原WENO分辨率的前提下,計算速度有明顯提高,減少了20%的計算時間。

    關(guān)鍵詞:高精度;WENO;RungeKutta;LaxWendroff;時間離散

    DOI:10.15938/j.jhust.2017.06.026

    中圖分類號: O175

    文獻(xiàn)標(biāo)志碼: A

    文章編號: 1007-2683(2017)06-0134-06

    Abstract:The research of high accuracy and high resolution schemes have been a hot topic in computational mathematics. According to low resolution and large amount of calculation of the original WENOJS scheme, we propose a simple new limiter fifth order upwind WENO scheme to reconstruct the numerical flux of the simple structure to improve the computational efficiency. Compared with other efficient high accuracy schemes such as ENO and WENO, it is shown that the computational cost of this scheme is less than that of WENOJS in the same accuracy. By use of MATLAB software, we compared and analyzed computational efficiencies and computational accuracies of LaxWendroff WENOJS scheme, LaxWendroff simple limiter WENO scheme, RungeKutta simple limiter WENO scheme and RungeKutta WENOJS scheme. The numerical results show that the new LaxWendroff simple limiter WENO scheme can improve the computing speed and reduce the computing time by 20% while maintaining the original WENO resolution.

    Keywords:high accuracy; WENO; RungeKutta; LaxWendroff; time discretization

    0 引 言

    雙曲守恒律方程(組)為科學(xué)理論和工程應(yīng)用研究中一類非常重要的偏微分方程(組)??諝鈩恿W(xué)、爆炸力學(xué)、流體力學(xué)等許多力學(xué)問題的求解都與其密切相關(guān)。由于即使在初始條件充分光滑的條件下,雙曲型守恒律方程的解仍可能出現(xiàn)間斷。因此,為了能夠更加高效地捕捉到間斷,并且避免間斷附近出現(xiàn)數(shù)值振蕩,計算流體力學(xué)及計算數(shù)學(xué)等領(lǐng)域一直致力于研究穩(wěn)定高效的數(shù)值計算方法[1-5]。

    低階精度的離散方法可能會對許多流動結(jié)構(gòu)造成“失真”模擬,不能真實反映真實的流動現(xiàn)象。高精度格式具有較小的耗散誤差和色散誤差,能夠更準(zhǔn)確地模擬雙曲守恒律方程的時空多尺度特性,在同樣網(wǎng)格條件下比低階格式能分辨出更加精細(xì)的流場、捕捉到其更細(xì)微的結(jié)構(gòu)變化。然而,高階精度格式也存在不足:求解在間斷處可能出現(xiàn)Gibbs現(xiàn)象即產(chǎn)生偽物理振蕩,從而導(dǎo)致非線性不穩(wěn)定現(xiàn)象。因此,發(fā)展高分辨率、高精度離散格式是計算流體力學(xué)發(fā)展中的迫切需要,成為計算流體力學(xué)工作者的一大研究方向[6-10]。

    TVD格式的出現(xiàn),標(biāo)志著計算流體力學(xué)步入了高精度計算格式階段,大量高精度計算格式不斷出現(xiàn),而WENO格式由于其良好的計算穩(wěn)定性,已經(jīng)成為一類極為重要的計算格式。眾所周知,經(jīng)典的WENOJS格式存在分辨率較低和計算量偏大兩方面問題,因此計算數(shù)學(xué)工作者嘗試對其改進,構(gòu)造了諸如WENOZ、WENOM等很多改進格式。根據(jù)文獻(xiàn)[11]的基本思想構(gòu)造的簡化限制器的5階WENO格式是對原WENO格式的很好的補充和修正。

    在時間離散方面,本文研究了LaxWendroff型時間離散格式,其具有較RungeKutta方法更為優(yōu)秀的計算效率。本文將WENO格式與LaxWendroff型時間離散相結(jié)合,得到了一種時空同步離散的新數(shù)值方法——基于LaxWendroff型時間離散的WENO格式。本文研究可以說是對高精度數(shù)值方法研究方面的有益補充和探索。

    1 WENO5格式

    針對以上兩種WENO格構(gòu)造不同的時間離散格式,能夠得到不同的數(shù)值方法。針對雙曲守恒律方程給出了具體的格式構(gòu)造過程,主要采用如下介紹的兩種時間離散方法LaxWendroff時間離散和TVD RungeKutta時間離散。下面,我們將介紹這兩種常用的時間離散方法。endprint

    2 時間離散格式

    目前雙曲守恒律方程最為常見的時間離散格式為TVD RungeKutta時間離散格式,也稱為SSP(strong stability preserving)RungeKutta時間離散。當(dāng)一階向前Euler迭代下的空間離散格式為TVD時,所構(gòu)造的全離散格式也TVD的。本文采用的RungeKutta時間離散格式為三階RungeKutta格式(簡記為RK3)[12-16]如下:

    3 結(jié)果與分析

    考慮線性對流方程和非線性方程兩種類型的方程的數(shù)值求解[17-20],利用LaxWendroff WENOJS格式、LaxWendroff簡單限制器WENO格式、RungeKutta WENOJS格式、RungeKutta簡單限制器的WENO格式四種格式進行數(shù)值模擬,研究數(shù)值方法的穩(wěn)定性、收斂性和計算效率等特性,所采用多的通量分裂為LaxFredrichs通量分裂。

    表4和表5分別給出LaxWendroff WENOJS格式、LaxWendroff簡單限制器WENO格式、RungeKutta WENOJS格式、RungeKutta簡單限制器的WENO格式計算結(jié)果。表4為四種格式的L1誤差和精度表,而表5為四種格式的L1誤差和精度表。

    具有如下初值u(x,0)=0.5+sin(πx),計算到t=32π時刻,邊界條件為周期條件,N=400個計算點。此時在1.2~1.3之間的位置產(chǎn)生了一強激波,圖1和圖2為其計算結(jié)果。

    計算結(jié)果中簡單限制器的WENO格式僅僅用到了2-3個計算點就可以識別激波,因此可以斷定簡單限制器的WENO格式的分辨率要高于WENOJS格式,而計算時間上,LaxWendroff時間離散具有一定的優(yōu)勢,WENO LW3和SWENO LW3兩種格式的計算時間為最少。

    4 結(jié) 論

    針對線性和非線性標(biāo)量守恒律方程進行了數(shù)值求解。研究了不同初值條件下的,LaxWendroff WENOJS格式、LaxWendroff簡單限制器WENO格式、RungeKutta WENOJS格式、RungeKutta簡單限制器的WENO格式的計算效率和計算精度。LaxWendroff簡單限制器WENO格式不但形式簡單且便于應(yīng)用到CFD代碼中。算例表明,LaxWendroff WENO格式的分辨率與RungeKutta型WENO格式較為接近,而LaxWendroff簡單限制器WENO格式具有更好的激波及小尺度波分辨率,特別計算速度提高了WENO5的20%以上。本文從分析標(biāo)量雙曲守恒律方程解的性質(zhì)入手,分析了雙曲守恒律方程自身的特點和原有的WENOJS差分格式的優(yōu)缺點后,提出一種簡化限制器的WENO格式,新格式較原格式的分辨率有所提高。LaxWendroff時間離散格式由于其計算效率方面的優(yōu)勢,可以作為RungeKutta方法的一種時間離散改進格式。本文研究了四種不同的數(shù)值方法:LaxWendroff WENOJS格式、LaxWendroff簡單限制器WENO格式、RungeKutta WENOJS格式、RungeKutta簡單限制器的WENO格式計算了一維標(biāo)量守恒型方程,并給出了數(shù)值實驗的結(jié)果。數(shù)值試驗表明這四種方法均出了高精度高分辨的特性,都具有更好的分辨激波捕捉的能力,在激波間斷處不但保持了陡峭的的形狀,而且沒有出現(xiàn)偽物理振蕩。我們發(fā)現(xiàn),對于同樣的問題,LaxWendroff簡單限制器WENO方法能在更短的時間和更小的儲存量。當(dāng)然,從數(shù)值實驗的結(jié)果來看,方法仍有改進的空間。未來需要關(guān)注的主要方面在于:在精度相同的情況下, LaxWendroff型的計算量要小于RungeKutta,因而更具有效率上的優(yōu)勢價值。然而,在推廣到雙曲守恒律方程組時,需要進行張量運算這給編程帶來了一定的困難。下一步作者將試著將推廣到高維守恒律方程方程組問題,現(xiàn)公式推到已經(jīng)完成,正在進行數(shù)值實驗,以得到較優(yōu)結(jié)果。

    參 考 文 獻(xiàn):

    [1] 陸金甫. 偏微分方程數(shù)值解法[M]. 北京:清華大學(xué)出版社, 2004.

    [2] 袁光偉,杭旭登,盛志強,等. 福射擴散計算方法若干研究進展[M]. 計算物理,2009,26(4): 475-500.

    [3] 歐陽頎. 非線性科學(xué)與斑圖動力學(xué)導(dǎo)論[M]. 北京:北京大學(xué)出版社,2010.

    [4] 劉儒勛,舒其望. 計算流體力學(xué)的若干新方法[M]. 北京:科學(xué)出版社,2003.

    [5] WARMING R F. Flux Vector Splitting of the Inviscid GasDynamic Equation with Application of Finite Difference Method [J]. Journal of Computational, 1981,40: 263-293.

    [6] GODONOV S K. A Difference Schemes for Numerical Computation of Discontinuous Solution of Hydrodynamics Equations [J]. Math. Sbornik, 1959,47: 271-306.

    [7] ENGQUIST,OSHER S. Stable and Entropy Satisfying Approximation for Transonic Flow Calculation [J]. Mathematics of Computation, 1980,34: 45-75.

    [8] COCKBURN B,LIN S Y,SHU C W. TVB RungeKutta Local Projection Discontinuous Galerkin Finite Element Method for Conservation Laws III: OneDimensional Systems [J]. J. Comput. Phys.,1989,84: 90-113.endprint

    [9] COCKBURN B,HOU S, SHU C W. TVB RungeKutta Local Projection Discontinuous Galerkin Finite element method for conservation laws IV [J]. Math. Comput.,1990,54: 545-581.

    [10]LIU H,XU H,GONG H. Modeling the Asymmetry in Traffic Flow (b):Macroscopic Approach [J]. Appl. Math. Model,2013,37(22): 9441-9450.

    [11]ZHU J,QIU J X. A New Fifth Order Finite Difference WENO Scheme for Solving Hyperbolic Conservation Laws [J]. J. Comput. Phys.,2016,318: 110-121.

    [12]LIU X, ZHANG S,ZHANG H,SHU C W. A New Class of Central Compact Schemes with Spectrallike Resolution II: Hybrid Weighted Nonlinear Schemes [J]. J. Comput. Phys.,2015,284: 133-154.

    [13]PIROZZOLI S. Conservative Hybrid CompactWENO Schemes for ShockTurbulence Interaction [J]. J. Comput. Phys.,2002,178 (1): 81-117.

    [14]QIU J X, DUMBSER M, SHU C W. The Discontinuous Galerkin Method with LaxWendroff type Time Discretizations [J]. Comput. Methods Appl. Mech. Engrg. 2005, 194: 4528-4543.

    [15]LAX P D. Hyperbolic System of Conservation Laws II [J]. Comm. Pure and Applied Mathematics. 1957,10: 537-566.

    [16]ZHANG X,SHU C W. On MaximumPrincipleSatisfying High Order Schemes for Scalar Conservation Laws [J]. J. Comput. Phys.,2010,229 (9): 3091-3120.

    [17]ZHANG X,SHU C W. MaximumPrincipleSatisfying and PositivityPreserving HighOrder Schemes for Conservation Laws: Survey and New Developments Proc. R. Soc. A [J]. Math. Phys. Eng. Sci.,2011,467(2): 2752-2776.

    [18]LIU Y J,SHU C W,TADMOR E,ZHANG M P. Central Discontinuous Galerkin Methods on Overlapping Cells with A NonOscillatory Hierarchical Reconstruction[J]. SIAM J. Numer. Anal.,2007, 45:2442-2467.

    [19]姚慧麗,宋曉秋,李興華. 一類半線性微分方程的漸進概自守溫和解 [J].哈爾濱理工大學(xué)學(xué)報,2012,17(1): 72-78.

    [20]姚慧麗,卜憲江,宋曉秋. 一類微分方程的指數(shù)增長的溫和漸進概自守解[J]. 哈爾濱理工大學(xué)學(xué)報,2014,19(5):23-26.

    (編輯:王 萍)endprint

    猜你喜歡
    高精度
    基于AD2S1210的高精度RDC電路設(shè)計
    電子測試(2018年23期)2018-12-29 11:11:30
    兩級升降伺服運動系統(tǒng)的高精度位置控制
    基于Niosll高精度超聲波流量計的研究
    電子制作(2018年17期)2018-09-28 01:56:38
    高精度PWM式DAC開發(fā)與設(shè)計
    智富時代(2018年8期)2018-09-28 00:12:34
    高抗擾高精度無人機著艦縱向飛行控制
    船載高精度星敏感器安裝角的標(biāo)定
    基于高精度測角的多面陣航測相機幾何拼接
    高精度免熱處理45鋼的開發(fā)
    山東冶金(2015年5期)2015-12-10 03:27:41
    高精度VLBI技術(shù)在深空探測中的應(yīng)用
    高精度可調(diào)式單端反激開關(guān)電源設(shè)計
    五月伊人婷婷丁香| 丰满少妇做爰视频| 欧美成人午夜精品| 精品一品国产午夜福利视频| 欧美激情高清一区二区三区 | 午夜av观看不卡| 免费av中文字幕在线| 2021少妇久久久久久久久久久| 中文字幕亚洲精品专区| 1024香蕉在线观看| 亚洲一码二码三码区别大吗| 亚洲成国产人片在线观看| av卡一久久| 女人被躁到高潮嗷嗷叫费观| 十八禁网站网址无遮挡| 午夜av观看不卡| 9色porny在线观看| 国产精品一二三区在线看| 久久精品aⅴ一区二区三区四区 | 一区二区三区四区激情视频| 欧美成人午夜精品| 亚洲第一青青草原| 亚洲人成电影观看| 搡老乐熟女国产| 国产黄色视频一区二区在线观看| 成年女人在线观看亚洲视频| 高清黄色对白视频在线免费看| 桃花免费在线播放| av国产久精品久网站免费入址| 免费黄网站久久成人精品| 久久久久精品人妻al黑| 观看av在线不卡| 欧美日韩视频高清一区二区三区二| 一级黄片播放器| 午夜福利一区二区在线看| 成人毛片60女人毛片免费| 亚洲精品乱久久久久久| 欧美 亚洲 国产 日韩一| 丁香六月天网| 熟女电影av网| 爱豆传媒免费全集在线观看| 婷婷色综合www| 极品少妇高潮喷水抽搐| 日韩熟女老妇一区二区性免费视频| 69精品国产乱码久久久| 啦啦啦在线免费观看视频4| 成人二区视频| 亚洲国产日韩一区二区| 久久精品久久久久久噜噜老黄| 久久久久久久精品精品| 国产精品久久久久成人av| 99久久人妻综合| 久久精品人人爽人人爽视色| 亚洲国产av影院在线观看| 亚洲成人手机| 极品少妇高潮喷水抽搐| 免费女性裸体啪啪无遮挡网站| 免费播放大片免费观看视频在线观看| 卡戴珊不雅视频在线播放| 999久久久国产精品视频| 岛国毛片在线播放| 人人妻人人澡人人看| 亚洲精品av麻豆狂野| 女人被躁到高潮嗷嗷叫费观| 美女高潮到喷水免费观看| 日韩不卡一区二区三区视频在线| 久久 成人 亚洲| 欧美 日韩 精品 国产| 久久精品国产亚洲av涩爱| 男人操女人黄网站| 国产精品女同一区二区软件| 蜜桃国产av成人99| 黄色配什么色好看| 亚洲欧洲国产日韩| 亚洲国产欧美网| 亚洲图色成人| 国产熟女欧美一区二区| av又黄又爽大尺度在线免费看| 婷婷色综合www| 国产又爽黄色视频| 久久久久久久久久久久大奶| 制服人妻中文乱码| 国产成人欧美| 国产成人精品在线电影| 亚洲av.av天堂| 午夜91福利影院| 国产亚洲av片在线观看秒播厂| 熟女电影av网| 国产精品三级大全| 麻豆乱淫一区二区| 在线天堂中文资源库| 日韩伦理黄色片| 亚洲婷婷狠狠爱综合网| 波多野结衣av一区二区av| 国产成人精品久久久久久| 又粗又硬又长又爽又黄的视频| 亚洲欧美中文字幕日韩二区| 国产伦理片在线播放av一区| 日韩中字成人| 性色avwww在线观看| 巨乳人妻的诱惑在线观看| 国产精品 国内视频| 亚洲国产精品一区三区| 久久久久久久大尺度免费视频| 国产一区二区三区综合在线观看| 亚洲国产av新网站| 在现免费观看毛片| 国产又色又爽无遮挡免| 制服诱惑二区| 亚洲国产欧美网| 少妇精品久久久久久久| 日韩av不卡免费在线播放| 久久午夜综合久久蜜桃| 午夜精品国产一区二区电影| 你懂的网址亚洲精品在线观看| 精品午夜福利在线看| 亚洲国产精品成人久久小说| 亚洲国产精品999| 99香蕉大伊视频| 久久国产亚洲av麻豆专区| 18禁动态无遮挡网站| 人人妻人人澡人人爽人人夜夜| 少妇被粗大猛烈的视频| 一级毛片我不卡| 日韩电影二区| 青草久久国产| 9191精品国产免费久久| 男女无遮挡免费网站观看| freevideosex欧美| 日日爽夜夜爽网站| 国产日韩欧美在线精品| 在线观看免费日韩欧美大片| 久久韩国三级中文字幕| 欧美在线黄色| 免费人妻精品一区二区三区视频| 国产成人免费无遮挡视频| 在线观看人妻少妇| 日韩av在线免费看完整版不卡| 亚洲情色 制服丝袜| 午夜精品国产一区二区电影| 寂寞人妻少妇视频99o| 91久久精品国产一区二区三区| 啦啦啦啦在线视频资源| 大陆偷拍与自拍| 有码 亚洲区| 天堂俺去俺来也www色官网| 免费大片黄手机在线观看| 亚洲综合色惰| 水蜜桃什么品种好| 久久鲁丝午夜福利片| 日日撸夜夜添| 一边摸一边做爽爽视频免费| 日韩熟女老妇一区二区性免费视频| 成年女人毛片免费观看观看9 | 在线观看国产h片| 国产欧美日韩一区二区三区在线| 熟女电影av网| 青春草亚洲视频在线观看| 久久久久久久国产电影| 大话2 男鬼变身卡| 国产综合精华液| 日韩制服丝袜自拍偷拍| 人成视频在线观看免费观看| 国产综合精华液| 97在线人人人人妻| 精品福利永久在线观看| 亚洲 欧美一区二区三区| 亚洲少妇的诱惑av| 久久免费观看电影| 大片免费播放器 马上看| 成年女人在线观看亚洲视频| av福利片在线| 在线 av 中文字幕| 国产探花极品一区二区| 精品久久蜜臀av无| 超碰97精品在线观看| 国产高清不卡午夜福利| www.熟女人妻精品国产| 亚洲图色成人| 亚洲av国产av综合av卡| 最近中文字幕高清免费大全6| 91成人精品电影| 2022亚洲国产成人精品| 国语对白做爰xxxⅹ性视频网站| 午夜福利乱码中文字幕| 男人操女人黄网站| 中文字幕精品免费在线观看视频| 涩涩av久久男人的天堂| 在线观看免费日韩欧美大片| 日韩免费高清中文字幕av| 国产精品国产av在线观看| 日韩精品有码人妻一区| 久久精品夜色国产| 99精国产麻豆久久婷婷| 成年人免费黄色播放视频| 亚洲一级一片aⅴ在线观看| 亚洲欧美日韩另类电影网站| 叶爱在线成人免费视频播放| 亚洲情色 制服丝袜| 丁香六月天网| 午夜av观看不卡| av不卡在线播放| 亚洲,欧美,日韩| 精品国产一区二区久久| 最近2019中文字幕mv第一页| 久久久国产精品麻豆| 97在线视频观看| 亚洲精品国产av蜜桃| 午夜福利网站1000一区二区三区| 久久韩国三级中文字幕| 18禁国产床啪视频网站| 免费不卡的大黄色大毛片视频在线观看| 麻豆精品久久久久久蜜桃| 亚洲av电影在线进入| 免费av中文字幕在线| av线在线观看网站| 日韩av不卡免费在线播放| 精品少妇内射三级| 熟女av电影| 免费不卡的大黄色大毛片视频在线观看| 热re99久久精品国产66热6| 国产熟女欧美一区二区| 亚洲精品久久午夜乱码| 精品一品国产午夜福利视频| 日本vs欧美在线观看视频| 国产日韩欧美视频二区| 亚洲欧美精品自产自拍| 国精品久久久久久国模美| 亚洲精品久久成人aⅴ小说| 亚洲av成人精品一二三区| 水蜜桃什么品种好| 精品国产乱码久久久久久男人| 制服人妻中文乱码| 欧美97在线视频| 亚洲美女黄色视频免费看| 日韩免费高清中文字幕av| 久久久欧美国产精品| 日韩制服骚丝袜av| 免费看av在线观看网站| 久久人妻熟女aⅴ| 国产一区二区 视频在线| 国产精品久久久久久精品古装| 成年动漫av网址| a级毛片黄视频| 妹子高潮喷水视频| 不卡视频在线观看欧美| 久久精品久久久久久噜噜老黄| videossex国产| 亚洲精品国产av蜜桃| 欧美最新免费一区二区三区| 色婷婷av一区二区三区视频| 咕卡用的链子| av在线app专区| 亚洲欧美一区二区三区黑人 | 天堂俺去俺来也www色官网| 亚洲精品美女久久av网站| 亚洲国产av新网站| 亚洲av免费高清在线观看| 成年动漫av网址| 国产精品国产三级专区第一集| 成人18禁高潮啪啪吃奶动态图| 香蕉精品网在线| 欧美日韩综合久久久久久| 亚洲av福利一区| 精品少妇一区二区三区视频日本电影 | 久久国产亚洲av麻豆专区| 久久久欧美国产精品| 久久久久视频综合| 欧美成人午夜免费资源| 欧美亚洲 丝袜 人妻 在线| 欧美人与性动交α欧美软件| 久久久久久久久久久免费av| 人人妻人人添人人爽欧美一区卜| 国产又色又爽无遮挡免| 亚洲内射少妇av| videos熟女内射| 成人黄色视频免费在线看| 亚洲欧美中文字幕日韩二区| 免费观看在线日韩| 色播在线永久视频| 亚洲男人天堂网一区| 欧美中文综合在线视频| 老女人水多毛片| 看非洲黑人一级黄片| 国产精品国产三级国产专区5o| 日本-黄色视频高清免费观看| 亚洲精品av麻豆狂野| 18禁国产床啪视频网站| 国产激情久久老熟女| 高清在线视频一区二区三区| 午夜av观看不卡| 亚洲精品第二区| 男女免费视频国产| 美国免费a级毛片| 久久久久国产网址| 黄片播放在线免费| 欧美激情 高清一区二区三区| 午夜久久久在线观看| av在线app专区| 久久久久精品久久久久真实原创| 亚洲国产欧美日韩在线播放| 一区二区三区精品91| 国产片特级美女逼逼视频| 免费高清在线观看日韩| 男男h啪啪无遮挡| 亚洲久久久国产精品| 欧美少妇被猛烈插入视频| 波多野结衣av一区二区av| 亚洲国产av影院在线观看| 欧美成人精品欧美一级黄| 18在线观看网站| 国产精品99久久99久久久不卡 | 国产精品人妻久久久影院| 国产深夜福利视频在线观看| 少妇人妻 视频| 国产无遮挡羞羞视频在线观看| 麻豆av在线久日| a级毛片黄视频| 日本-黄色视频高清免费观看| 午夜日韩欧美国产| 成人免费观看视频高清| 日日爽夜夜爽网站| 国产成人精品福利久久| 免费观看性生交大片5| 尾随美女入室| 亚洲情色 制服丝袜| 免费看不卡的av| 亚洲精品国产一区二区精华液| 丁香六月天网| 亚洲成国产人片在线观看| 女人高潮潮喷娇喘18禁视频| 日本wwww免费看| 狠狠精品人妻久久久久久综合| 日本av免费视频播放| 国产精品一二三区在线看| freevideosex欧美| 日韩中文字幕欧美一区二区 | 黄频高清免费视频| 高清欧美精品videossex| 亚洲av电影在线进入| 91精品国产国语对白视频| 捣出白浆h1v1| av国产精品久久久久影院| 亚洲视频免费观看视频| 成人18禁高潮啪啪吃奶动态图| 欧美成人午夜免费资源| 丝袜人妻中文字幕| 黄网站色视频无遮挡免费观看| 美女国产视频在线观看| 亚洲一区二区三区欧美精品| 免费人妻精品一区二区三区视频| 男人爽女人下面视频在线观看| 成人亚洲精品一区在线观看| 免费观看无遮挡的男女| 精品酒店卫生间| 久久 成人 亚洲| www.av在线官网国产| 校园人妻丝袜中文字幕| 黄色配什么色好看| 大片免费播放器 马上看| 久久精品人人爽人人爽视色| 精品卡一卡二卡四卡免费| 黄色一级大片看看| 精品卡一卡二卡四卡免费| 人人妻人人澡人人爽人人夜夜| av卡一久久| 90打野战视频偷拍视频| 久久午夜福利片| 日本黄色日本黄色录像| 黄色毛片三级朝国网站| 久久国产精品大桥未久av| 欧美日韩国产mv在线观看视频| 国产高清不卡午夜福利| 欧美人与性动交α欧美软件| 亚洲精品aⅴ在线观看| 下体分泌物呈黄色| 麻豆av在线久日| 韩国av在线不卡| 久久99精品国语久久久| 国产无遮挡羞羞视频在线观看| 亚洲色图 男人天堂 中文字幕| 超碰97精品在线观看| 波多野结衣一区麻豆| 18禁动态无遮挡网站| 日韩伦理黄色片| 欧美成人午夜免费资源| 亚洲精品乱久久久久久| 成年女人在线观看亚洲视频| 国产人伦9x9x在线观看 | 国产免费视频播放在线视频| 久久久久久久久久人人人人人人| 亚洲四区av| 永久免费av网站大全| 欧美精品亚洲一区二区| 免费高清在线观看视频在线观看| 香蕉精品网在线| 又大又黄又爽视频免费| 免费黄频网站在线观看国产| www.精华液| 午夜免费男女啪啪视频观看| 久久久久久久久免费视频了| av有码第一页| 婷婷色综合www| 久久99一区二区三区| 亚洲国产色片| 大香蕉久久成人网| 欧美日韩亚洲高清精品| 少妇人妻久久综合中文| 免费少妇av软件| av网站在线播放免费| 日韩一区二区视频免费看| 国产视频首页在线观看| 国产一区二区三区综合在线观看| 9191精品国产免费久久| 国产成人精品无人区| 狂野欧美激情性bbbbbb| 成人黄色视频免费在线看| 18禁观看日本| 免费不卡的大黄色大毛片视频在线观看| 日韩视频在线欧美| 久久影院123| 久久久久精品人妻al黑| 午夜91福利影院| 国产成人精品久久久久久| 777久久人妻少妇嫩草av网站| 欧美精品一区二区大全| 精品国产超薄肉色丝袜足j| 国产熟女午夜一区二区三区| 99久久人妻综合| 可以免费在线观看a视频的电影网站 | 欧美变态另类bdsm刘玥| 我的亚洲天堂| 亚洲经典国产精华液单| 欧美 亚洲 国产 日韩一| 电影成人av| 激情视频va一区二区三区| 少妇被粗大猛烈的视频| 午夜福利视频在线观看免费| 大片电影免费在线观看免费| 亚洲第一青青草原| 亚洲成人手机| 边亲边吃奶的免费视频| 日韩中文字幕欧美一区二区 | 欧美亚洲 丝袜 人妻 在线| 香蕉精品网在线| 日本vs欧美在线观看视频| 日韩视频在线欧美| 国产精品女同一区二区软件| 国产成人精品在线电影| 亚洲精品日本国产第一区| 18在线观看网站| 中文字幕亚洲精品专区| 亚洲一级一片aⅴ在线观看| h视频一区二区三区| 制服人妻中文乱码| 亚洲经典国产精华液单| av女优亚洲男人天堂| 久久精品亚洲av国产电影网| 久久久久久人妻| 成年女人毛片免费观看观看9 | 国产无遮挡羞羞视频在线观看| 亚洲美女视频黄频| 男女无遮挡免费网站观看| 成人漫画全彩无遮挡| 国产成人精品福利久久| 99九九在线精品视频| 国产成人精品久久二区二区91 | 欧美成人午夜免费资源| 精品一区二区免费观看| 欧美在线黄色| 综合色丁香网| 欧美激情高清一区二区三区 | 亚洲精华国产精华液的使用体验| 免费观看性生交大片5| 精品国产露脸久久av麻豆| 亚洲欧美一区二区三区国产| av福利片在线| 久久国产精品男人的天堂亚洲| 久久久久人妻精品一区果冻| 久久女婷五月综合色啪小说| 国产 精品1| 欧美日韩视频精品一区| 精品亚洲成国产av| 在线观看国产h片| 成人亚洲精品一区在线观看| 人妻少妇偷人精品九色| 日本午夜av视频| 少妇精品久久久久久久| 成人漫画全彩无遮挡| 亚洲av电影在线进入| 欧美日韩av久久| 日本午夜av视频| 国产一区有黄有色的免费视频| 伦理电影免费视频| 国产精品一区二区在线观看99| 国产精品熟女久久久久浪| 国产黄频视频在线观看| 国产亚洲av片在线观看秒播厂| 精品一区在线观看国产| 精品酒店卫生间| 久久久久网色| 亚洲色图 男人天堂 中文字幕| 老女人水多毛片| 制服人妻中文乱码| 成人亚洲欧美一区二区av| 最近2019中文字幕mv第一页| 人人妻人人澡人人看| 日韩精品免费视频一区二区三区| 日韩成人av中文字幕在线观看| 香蕉丝袜av| 青春草国产在线视频| 国产成人91sexporn| 国产午夜精品一二区理论片| 色播在线永久视频| 亚洲av中文av极速乱| 欧美成人午夜免费资源| 欧美另类一区| 国产白丝娇喘喷水9色精品| 男人爽女人下面视频在线观看| 成人黄色视频免费在线看| tube8黄色片| 精品第一国产精品| 日本免费在线观看一区| 免费看av在线观看网站| 国产一区亚洲一区在线观看| 国产成人一区二区在线| av电影中文网址| 这个男人来自地球电影免费观看 | 人人妻人人澡人人爽人人夜夜| 97人妻天天添夜夜摸| 精品国产乱码久久久久久小说| 久久狼人影院| 少妇猛男粗大的猛烈进出视频| 欧美精品人与动牲交sv欧美| 精品99又大又爽又粗少妇毛片| av网站在线播放免费| 黄色配什么色好看| 国产精品99久久99久久久不卡 | 欧美日韩国产mv在线观看视频| 视频在线观看一区二区三区| 乱人伦中国视频| 亚洲欧洲精品一区二区精品久久久 | 宅男免费午夜| 一区二区三区乱码不卡18| 亚洲精品,欧美精品| 国产日韩欧美亚洲二区| 久久久久人妻精品一区果冻| 香蕉丝袜av| 亚洲精品中文字幕在线视频| 一级片'在线观看视频| 精品酒店卫生间| 国产精品免费视频内射| 婷婷色av中文字幕| 亚洲av福利一区| 九草在线视频观看| 最近2019中文字幕mv第一页| 少妇人妻久久综合中文| 国产毛片在线视频| 国产精品嫩草影院av在线观看| 日韩中字成人| 久久久精品免费免费高清| 久久午夜福利片| 91国产中文字幕| 免费少妇av软件| 免费观看无遮挡的男女| 麻豆av在线久日| 国产亚洲精品第一综合不卡| 综合色丁香网| 精品福利永久在线观看| 国产在线一区二区三区精| 亚洲av电影在线观看一区二区三区| 国产有黄有色有爽视频| 国产精品久久久久久久久免| 亚洲欧洲日产国产| 99国产综合亚洲精品| 婷婷色av中文字幕| 免费久久久久久久精品成人欧美视频| 90打野战视频偷拍视频| 丰满饥渴人妻一区二区三| 久久久久视频综合| 国产成人精品福利久久| 一级黄片播放器| 亚洲一级一片aⅴ在线观看| 欧美日本中文国产一区发布| 一二三四在线观看免费中文在| 女人精品久久久久毛片| 国产一区二区三区综合在线观看| 九色亚洲精品在线播放| 99久久中文字幕三级久久日本| 不卡视频在线观看欧美| 亚洲欧美精品综合一区二区三区 | 亚洲精品国产av蜜桃| 亚洲欧洲精品一区二区精品久久久 | 中文字幕人妻丝袜一区二区 | 久热这里只有精品99| 亚洲欧美一区二区三区久久| 成人黄色视频免费在线看| 精品人妻在线不人妻| 一级毛片我不卡| 国产熟女午夜一区二区三区| av视频免费观看在线观看| 婷婷色综合www| 国产欧美日韩综合在线一区二区| 亚洲精华国产精华液的使用体验| 亚洲四区av| 久久久精品国产亚洲av高清涩受| 欧美黄色片欧美黄色片| 丝袜脚勾引网站| 亚洲伊人久久精品综合| 久久久久精品人妻al黑| 侵犯人妻中文字幕一二三四区| 国产一区亚洲一区在线观看| 成人亚洲精品一区在线观看| 妹子高潮喷水视频| 男女边吃奶边做爰视频|