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

    KdV-Burgers方程的一類新本性并行差分格式*

    2023-06-06 07:26:44潘悅悅楊曉忠
    應用數(shù)學和力學 2023年5期
    關(guān)鍵詞:計算精度差分解析

    潘悅悅, 楊曉忠

    (1.華北電力大學 控制與計算機工程學院, 北京 102206;2.華北電力大學 數(shù)理學院 信息與計算研究所, 北京 102206)

    0 引 言

    KdV-Burgers方程是一種帶阻尼的KdV方程模型,也是一個典型的非線性波動方程.其一方面,可以描述同時含有色散、非線性水平對流、黏性流體的淺水長波的傳播過程[1-2];另一方面,可以用來描述同時含有色散和耗散的等離子物理中的一些重要物理現(xiàn)象[3-5].KdV-Burgers方程求解的快速算法一直是計算數(shù)學研究的熱點問題,其高效數(shù)值解法具有重要的科學意義和應用價值.

    隨著越來越多大規(guī)??茖W計算問題的出現(xiàn),人們開始對串行計算程序的合理化和并行化進行深入研究.近年來的大多數(shù)并行差分方法是條件穩(wěn)定的,或者是無條件穩(wěn)定但精度只有一階.為了得到穩(wěn)定性條件更寬松、計算精度更高的并行方法,實現(xiàn)串行差分方法的合理并行化,一種新的高效并行方法是必不可少的.針對拋物型方程,Evans和Abdullah[6]提出了分組顯式(group explicit,GE)的思想,設計出交替分組顯式(AGE)差分格式.通過計算網(wǎng)格區(qū)域中的兩個相鄰點,將得到的隱式方程轉(zhuǎn)化為顯式方程.交替使用這種方法,截斷誤差項可以部分抵消,從而提高了算法的精度.

    隱格式具有良好的穩(wěn)定性,但不適合并行計算.受到構(gòu)造AGE方法的啟示,張寶琳等[7]提出了用Saul’yev非對稱格式構(gòu)造分段隱式的思想,并且恰當?shù)厥褂媒惶婕夹g(shù)建立了多種交替分段顯-隱式(ASE-I)和交替分段Crank-Nicolson(ASC-N)并行方法,得到了穩(wěn)定性和并行性兼顧的研究成果.周毓麟院士[8]將最一般的拋物型方程的顯隱混合格式稱之為具有并行本性的差分格式,對其差分解的存在性、唯一性、收斂性和穩(wěn)定性等理論問題進行了研究,建立了拋物方程并行本性差分方法的基本理論.現(xiàn)在,并行本性差分方法已經(jīng)推廣到很多發(fā)展方程的數(shù)值求解中[9-10].

    王文洽[11]使用一組非對稱格式構(gòu)造了具有本性并行的ASC-N差分格式,求解了KdV方程,給出了線性絕對穩(wěn)定性的理論分析.Yuan等[12]針對非線性拋物型方程組提出了一類無條件穩(wěn)定,空間二階精度的并行差分格式.Borhanifar和Abazari[13]應用無條件穩(wěn)定的差分格式求解了電報方程.數(shù)值算例表明了新方法的穩(wěn)定性和準確性.針對四階拋物型方程,Guo和Liu[14]給出了經(jīng)典的顯式、隱式和C-N格式,主要分析了格式的無條件穩(wěn)定性.Namjoo等[15]提出了廣義Burgers-Fisher方程的3種格式:兩種精確有限差分格式和一種非標準有限差分格式.通過誤差分析,驗證了該格式的精確性和有效性.Xue和Feng[16]建立了求解拋物型方程的有效算法.在第n層使用兩種區(qū)域分解法,兩種方法的解取平均值,得到第n+1層的解.該方法保持了并行性和穩(wěn)定性.并行本性方法既適用于并行計算,又是無條件穩(wěn)定或線性穩(wěn)定的.然而,這些方法誤差稍大,精度較低[17-18].

    受到上述差分方法的啟發(fā),本文結(jié)合經(jīng)典的顯式、隱式和C-N格式,提出了一種求解KdV-Burgers方程的新并行差分方法.MASC-N方法具有更精確的結(jié)果和更寬松的穩(wěn)定性條件.理論分析和數(shù)值算例表明,MASC-N格式線性絕對穩(wěn)定,收斂階為二階.本文格式的計算效率遠高于C-N格式的計算效率,新方法對于求解KdV-Burgers方程是高效合理的.

    1 KdV-Burgers方程的混合交替分段C-N差分格式

    1.1 KdV-Burgers方程

    KdV-Burgers方程的一般形式為[1-2]

    ut+αuux+βuxx+εuxxx=0,L1

    (1)

    其中α>0,β<0,ε>0.初始條件和邊界條件為

    u(x,0)=f(x),L1≤x≤L2,

    (2)

    u(L1,t)=g1(t),u(L2,t)=g2(t), 0

    (3)

    其中L1,L2都是適當大的數(shù).

    1.2 MASC-N并行差分格式的構(gòu)造

    以下是方程(1)的顯格式、隱格式和C-N格式:

    (4)

    (5)

    (6)

    其中

    現(xiàn)在討論KdV-Burgers方程的交替分段并行差分格式,本文的MASC-N格式(圖1)將結(jié)合以上幾種差分格式使用.

    圖1 MASC-N格式構(gòu)造示意圖

    (7)

    式中

    Un,F1,F2是M-3維向量,I是M-3階單位矩陣,A1和A2均是M-3階對角矩陣,且滿足A1+A2=I,A1=diag(θ1,θ2,…,θM-4,θM-3),

    2 MASC-N并行差分方法的數(shù)值分析

    2.1 MASC-N格式解的存在唯一性

    引理1(Kellogg引理)[19]設θ>0,矩陣G+GT是非負定的,則(θI+G)-1存在,并且有

    ‖(θI+G)-1‖2≤θ-1.

    引理2 MASC-N格式的A1G和A2G是非負定矩陣.

    其中c<0.顯然A1G+(A1G)T是弱對角占優(yōu)矩陣,對角線上的元素是非負實數(shù).因此A1G+(A1G)T是非負定矩陣,A2G+(A2G)T也是非負定矩陣,則A1G和A2G是非負定矩陣.

    由初始條件得,U0已知,假設U2n已知,求U2n+1.由式(7)得到求解U2n+1的公式為

    (I+A1G)U2n+1=(I-A2G)U2n+F1.

    (8)

    顯然,等式右邊已知,(I+A1G)-1存在,式(8)有唯一解.同樣地,使用MASC-N格式計算2n+2層的點:

    (I+A2G)U2n+2=(I-A1G)U2n+1+F2.

    (9)

    同理可得式(9)有唯一解,定理得證.

    定理1 KdV-Burgers方程MASC-N格式(7)的解是存在且唯一的.

    2.2 MASC-N格式的線性絕對穩(wěn)定性

    引理3[19]設θ>0,矩陣G+GT是非負定的,則‖(θI-G)(θI+G)-1‖2≤1.

    定理2 KdV-Burgers方程的MASC-N格式(7)是線性絕對穩(wěn)定的.

    2.3 MASC-N格式的收斂性

    古典顯格式為

    古典隱格式為

    由圖1可知,在每段的“內(nèi)點”使用C-N格式,其計算精度為二階.在每段的“內(nèi)邊界點”使用古典顯式和古典隱式,計算精度為O(τ+h2τ+h2).在T1(τ,h)和T2(τ,h)的表達式中,分別包含了絕對值相同但符號相反的項.古典顯式和古典隱式在每一時間層交替使用,部分項將相互抵消.因此,MASC-N格式的計算精度近似為二階.

    定理3 KdV-Burgers方程的MASC-N格式(7)的計算精度為O(h2+τ2),內(nèi)邊界點的計算精度為O(τ+h2τ+h2).

    ‖en+2‖≤‖(I+G2)-1(I-G1)(I+G1)-1(I-G2)en‖+

    ‖(I+G2)-1(I-G1)(I+G1)-1Rn+1‖+‖(I+G2)-1Rn+2‖≤

    C′n+2(τ2+h2).

    定理4 KdV-Burgers方程的MASC-N格式(7)是收斂的.

    3 數(shù) 值 試 驗

    為了驗證理論分析,給出數(shù)值算例證明MASC-N格式的計算精度、穩(wěn)定性、收斂性和并行性.同時給出MASC-N格式與ASE-I格式和ASC-N格式的比較試驗、MASC-N格式擴展到mKdV-Burgers方程的數(shù)值試驗.

    例1 令ε=1,γ=1,μ=1,考慮以下KdV-Burgers方程[20-21]:

    ut+εuux-γuxx+μuxxx=0, -50≤x≤50, 0≤t≤1.

    (10)

    初始條件為

    邊界條件為

    解析解為

    取時間層N=500,空間層M=500,圖2和3分別給出了式(10)的孤立波解從時間t=0到t=1的精確解和MASC-N格式解的波形圖.從圖中可以看到,數(shù)值結(jié)果穩(wěn)定可靠,保持了與精確解幾乎完全一致的波形.

    圖2 精確解的波形 圖3 MASC-N格式解的波形

    表1 兩種格式解相對于解析解的絕對誤差

    以下給出C-N格式和MASC-N格式收斂階的數(shù)值試驗.在點(xi,tn)處,定義誤差E∞(h,τ)、空間收斂階S和時間收斂階T[22-23]:

    令M=25,50,100,200,400,N=500和N=100,200,400,800,1 600,M=100,表2和3分別給出了C-N格式和MASC-N格式的空間和時間精度.由表中數(shù)據(jù)可知,兩種格式的空間和時間精度為二階.因此本文MASC-N 格式和已有的C-N格式有相同的階數(shù),數(shù)值試驗結(jié)果與理論分析相符.

    表2 兩種格式的空間收斂階

    表3 兩種格式的時間收斂階

    定義加速比Sp=T1/Tp和效率Ep=Sp/p[24],T1和Tp分別是C-N和MASC-N格式的運行時間,p是處理器的數(shù)量.使用八核處理器,令N=1 000,M=600,1 100,1 600,2 100,2 600,計算結(jié)果見表4.表4表明當M增加時,兩種格式的運行時間都增加.但與C-N格式相比,MASC-N格式的運行時間增加緩慢.Sp和Ep都逐漸增加.隨著M的增大,MASC-N格式的效率大幅度提高.當M>600后,MASC-N 格式的運行時間比C-N格式節(jié)約56%.

    表4 兩種格式運行時間的比較

    圖4給出了C-N和MASC-N格式的計算時間.可以看到C-N格式的計算時間迅速增加,變化率大.本文構(gòu)造的MASC-N格式的計算時間更少,且變化率小.C-N格式是串行格式,通過LU分解計算5對角方程.但使用MASC-N 格式計算時,全局問題被分解成多個不相關(guān)的問題同步進行計算.程序循環(huán)有效地提高了計算效率,縮短了計算時間.從時間的角度來看,MASC-N格式具有明顯的優(yōu)勢.

    圖4 C-N和MASC-N格式的計算時間 圖5 不同CPU核數(shù)的MASC-N格式計算時間

    圖5給出了MASC-N格式在四核和八核CPU下的計算時間.可以看到隨著空間網(wǎng)格數(shù)的增加,八核CPU的計算時間比四核CPU更少,MASC-N格式并行性明顯.

    為了驗證MASC-N格式的穩(wěn)定性和精度,給出隨時間和空間變化的相對誤差.定義時間層相對誤差和(the sum of relative error for every time level,SRET,ΔSRET)與誤差能量(difference total energy,DTE,ΔDTE):

    取N=200,M=200,圖6表明當N增加時,兩種格式解的ΔSRET迅速減小,且趨向于0.因此,MASC-N方法求解KdV-Burgers方程是計算穩(wěn)定的.

    圖6 C-N和MASC-N格式解的SRET 圖7 C-N和MASC-N格式解的DTE

    取N=500,M=500,從圖7可以看出兩種格式ΔDTE的最大值是5×10-3.表明兩種格式數(shù)值解與解析解非常接近,誤差很小,MASC-N方法有良好的計算精度.

    例2 令L=100,T=1,考慮以下KdV-Burgers方程[25]:

    ut+uux-uxx+uxxx=0.

    (11)

    初邊界值由解析解給出,解析解為

    比較3種并行格式的計算時間.固定時間層N=10 000,計算結(jié)果見圖8和表5.ASC-N格式有最長的運行時間.ASE-I和MASC-N格式運行時間相似且MASC-N格式計算時間略短.

    表5 3種并行格式的計算時間

    圖8 3種并行格式計算時間比較

    比較解析解和ASE-I[9]、ASC-N[11]、MASC-N格式解.令N=210,M=210,t=0.5,表6給出了計算結(jié)果.表中數(shù)據(jù)表明3種格式的數(shù)值解相似,且MASC-N格式解最接近解析解.

    表6 數(shù)值解和解析解

    最后,取計算精度的權(quán)重為0.6,計算時間的權(quán)重為0.4.根據(jù)計算精度和時間比較3種并行格式,見表7.可以看到,MASC-N格式的加權(quán)排序總和數(shù)最小,可知其具有最好的性能和實用性.

    表7 3種并行格式的比較

    例3 令β=-1,ε=-1,考慮以下mKdV-Burgers方程[26-27]:

    ut+u2ux+βuxx+εuxxx=0, -50≤x≤50, 0≤t≤1.

    (12)

    初邊界值由解析解給出,解析解為

    令N=210,M=210,t=0.5,表8為mKdV-Burgers方程解析解和兩種格式數(shù)值解的比較.可以看出C-N格式解和MASC-N格式解幾乎相同.圖9和10分別是C-N和MASC-N格式的絕對誤差,對比發(fā)現(xiàn)兩種差分格式絕對誤差最大值的級數(shù)為10-2.

    表8 數(shù)值解和解析解

    圖9 C-N格式解的絕對誤差 圖10 MASC-N格式解的絕對誤差

    令時間層N=1 000,圖11是C-N和MASC-N格式的運行時間.可以看到C-N格式的運行時間遠大于MASC-N格式的運行時間.與C-N格式相比,MASC-N格式求解mKdV-Burgers方程更為高效.

    圖11 C-N和MASC-N格式的運行時間

    4 結(jié) 論

    KdV-Burgers方程的MASC-N并行差分格式是線性絕對穩(wěn)定的,格式解存在唯一,與C-N格式的收斂階同為二階.在保證計算精度的前提下,MASC-N格式比C-N格式有更高的效率.與ASE-I和ASC-N格式相比,MASC-N并行差分格式對求解KdV-Burgers方程有更好的性能.同時,將MASC-N并行差分格式擴展到求解mKdV-Burgers方程,本文的數(shù)值試驗結(jié)果表明,MASC-N格式的計算效率明顯高于C-N格式的計算效率,MASC-N 格式求解各種類型的并行計算系統(tǒng)是有效的.

    致謝本文作者衷心感謝德國阿爾弗雷德韋格納研究所Sergey Danilov教授和王強博士在寫作過程中提出的寶貴建議和意見;同時衷心感謝華北電力大學國內(nèi)外聯(lián)合培養(yǎng)博士生項目(2020)對本文的資助.

    猜你喜歡
    計算精度差分解析
    數(shù)列與差分
    三角函數(shù)解析式中ω的幾種求法
    睡夢解析儀
    基于SHIPFLOW軟件的某集裝箱船的阻力計算分析
    廣東造船(2018年1期)2018-03-19 15:50:50
    電競初解析
    商周刊(2017年12期)2017-06-22 12:02:01
    相機解析
    單元類型和尺寸對拱壩壩體應力和計算精度的影響
    價值工程(2015年9期)2015-03-26 06:40:38
    基于差分隱私的大數(shù)據(jù)隱私保護
    鋼箱計算失效應變的沖擊試驗
    相對差分單項測距△DOR
    太空探索(2014年1期)2014-07-10 13:41:50
    国产国拍精品亚洲av在线观看| 国产色婷婷99| 久久久精品欧美日韩精品| 美女黄网站色视频| 久久99热这里只有精品18| 久久精品夜夜夜夜夜久久蜜豆| 久久国产乱子免费精品| 日日啪夜夜撸| 日韩欧美 国产精品| 在线观看av片永久免费下载| 日韩人妻高清精品专区| 人人妻人人澡欧美一区二区| 黄片无遮挡物在线观看| 丰满人妻一区二区三区视频av| 亚洲伊人久久精品综合 | 一区二区三区免费毛片| 日本与韩国留学比较| 久久久午夜欧美精品| 国产乱人视频| 国产一区二区三区av在线| 国产 一区精品| 久久亚洲精品不卡| 亚洲一级一片aⅴ在线观看| 国产国拍精品亚洲av在线观看| 国产精品三级大全| 美女xxoo啪啪120秒动态图| 嫩草影院新地址| 男插女下体视频免费在线播放| 1024手机看黄色片| 国语对白做爰xxxⅹ性视频网站| 国产高潮美女av| 国产老妇女一区| 久久久亚洲精品成人影院| 日本五十路高清| 3wmmmm亚洲av在线观看| 国产一区二区在线观看日韩| 亚洲欧洲国产日韩| 久久欧美精品欧美久久欧美| 边亲边吃奶的免费视频| av.在线天堂| 午夜福利网站1000一区二区三区| 成年女人看的毛片在线观看| 久久久久国产网址| 看十八女毛片水多多多| 人妻夜夜爽99麻豆av| av在线蜜桃| 日韩欧美精品v在线| 成人特级av手机在线观看| 毛片女人毛片| 中文字幕人妻熟人妻熟丝袜美| 黄色配什么色好看| 日韩av不卡免费在线播放| 丰满人妻一区二区三区视频av| 好男人视频免费观看在线| 国产精品久久久久久av不卡| 99久久无色码亚洲精品果冻| 人人妻人人澡欧美一区二区| 最新中文字幕久久久久| 亚洲精品色激情综合| 一区二区三区四区激情视频| 99热网站在线观看| 久久久久久久午夜电影| 国产精品久久视频播放| 麻豆精品久久久久久蜜桃| 亚洲国产色片| 国产视频内射| 极品教师在线视频| 最近最新中文字幕大全电影3| 久久人妻av系列| 久久久欧美国产精品| 人妻制服诱惑在线中文字幕| 在线天堂最新版资源| 日韩av不卡免费在线播放| 国产男人的电影天堂91| 一本一本综合久久| 国产免费又黄又爽又色| 亚洲第一区二区三区不卡| 99久久中文字幕三级久久日本| 久久午夜福利片| av线在线观看网站| 老司机影院毛片| 亚洲欧美成人综合另类久久久 | 久久久欧美国产精品| 一区二区三区高清视频在线| 九色成人免费人妻av| 日韩强制内射视频| 国产亚洲av嫩草精品影院| 国产乱来视频区| 观看免费一级毛片| av天堂中文字幕网| 精品一区二区免费观看| 午夜老司机福利剧场| 久久久久久久久久久丰满| 国产黄片美女视频| 久久这里只有精品中国| 97热精品久久久久久| 高清毛片免费看| 亚洲成人av在线免费| 日日撸夜夜添| 国产精华一区二区三区| 99热这里只有是精品50| 中国美白少妇内射xxxbb| 亚洲中文字幕日韩| www.av在线官网国产| 亚洲美女搞黄在线观看| 午夜精品一区二区三区免费看| 欧美日韩国产亚洲二区| 国产高潮美女av| 国产精品一及| 男女边吃奶边做爰视频| 精品一区二区三区视频在线| 国产一区有黄有色的免费视频 | 一卡2卡三卡四卡精品乱码亚洲| 伊人久久精品亚洲午夜| 麻豆久久精品国产亚洲av| 亚洲伊人久久精品综合 | 亚洲精品乱码久久久v下载方式| 麻豆一二三区av精品| 两性午夜刺激爽爽歪歪视频在线观看| 国产极品精品免费视频能看的| 在线观看一区二区三区| 国国产精品蜜臀av免费| 午夜福利视频1000在线观看| 99久国产av精品| 九草在线视频观看| 亚洲综合色惰| 久久久国产成人免费| 97超碰精品成人国产| 五月伊人婷婷丁香| 久久久久九九精品影院| 久久精品熟女亚洲av麻豆精品 | 免费播放大片免费观看视频在线观看 | 男女国产视频网站| 女人被狂操c到高潮| 99久久九九国产精品国产免费| 精品一区二区免费观看| 亚洲国产色片| 伦精品一区二区三区| 在线观看一区二区三区| 99久久精品国产国产毛片| av专区在线播放| 日产精品乱码卡一卡2卡三| 国产黄片视频在线免费观看| 99久久成人亚洲精品观看| 永久网站在线| 我要看日韩黄色一级片| 欧美成人免费av一区二区三区| 国产熟女欧美一区二区| 国产视频首页在线观看| 麻豆精品久久久久久蜜桃| 精品酒店卫生间| 99久久中文字幕三级久久日本| 插阴视频在线观看视频| av专区在线播放| 午夜福利高清视频| 亚洲第一区二区三区不卡| 亚洲欧美精品专区久久| 国产视频首页在线观看| 久久久久性生活片| 晚上一个人看的免费电影| 亚洲自偷自拍三级| 欧美激情在线99| 男人和女人高潮做爰伦理| 国产中年淑女户外野战色| 十八禁国产超污无遮挡网站| 尾随美女入室| 亚洲av日韩在线播放| 欧美日本亚洲视频在线播放| 免费播放大片免费观看视频在线观看 | 老师上课跳d突然被开到最大视频| 亚洲欧美成人综合另类久久久 | 亚洲18禁久久av| 亚洲国产欧洲综合997久久,| 人妻少妇偷人精品九色| 美女被艹到高潮喷水动态| 亚洲成人精品中文字幕电影| 一夜夜www| 99久久成人亚洲精品观看| 精品酒店卫生间| 长腿黑丝高跟| 国产精品嫩草影院av在线观看| 日韩在线高清观看一区二区三区| 中文欧美无线码| 国模一区二区三区四区视频| 九九爱精品视频在线观看| 日韩欧美 国产精品| 亚洲18禁久久av| 男插女下体视频免费在线播放| a级一级毛片免费在线观看| 久久久久久九九精品二区国产| 视频中文字幕在线观看| 欧美丝袜亚洲另类| 天天一区二区日本电影三级| 哪个播放器可以免费观看大片| 日本av手机在线免费观看| 直男gayav资源| 国产黄片美女视频| 成人特级av手机在线观看| 亚洲人与动物交配视频| av黄色大香蕉| 国产女主播在线喷水免费视频网站 | 国产av在哪里看| 日本爱情动作片www.在线观看| 久久精品国产99精品国产亚洲性色| 亚洲成人av在线免费| 日本wwww免费看| 午夜爱爱视频在线播放| 美女被艹到高潮喷水动态| 久久精品国产自在天天线| 两个人视频免费观看高清| 亚洲国产精品成人久久小说| 亚洲国产日韩欧美精品在线观看| 天天躁夜夜躁狠狠久久av| 久久久精品94久久精品| 成人鲁丝片一二三区免费| 色5月婷婷丁香| 国产精品一及| 老司机影院成人| 国产探花极品一区二区| 精品久久久久久久久亚洲| av女优亚洲男人天堂| 国产在线男女| 久久精品夜夜夜夜夜久久蜜豆| 日韩高清综合在线| 久久久欧美国产精品| 村上凉子中文字幕在线| 97在线视频观看| 老司机福利观看| 99久久九九国产精品国产免费| 97在线视频观看| 非洲黑人性xxxx精品又粗又长| 亚洲成人av在线免费| 亚洲精品色激情综合| 免费看av在线观看网站| 国产伦在线观看视频一区| 亚洲精品aⅴ在线观看| 亚洲欧洲日产国产| 亚洲精品一区蜜桃| 欧美+日韩+精品| 22中文网久久字幕| 国模一区二区三区四区视频| 日本三级黄在线观看| 日韩一本色道免费dvd| 精品午夜福利在线看| 最近的中文字幕免费完整| 白带黄色成豆腐渣| 美女黄网站色视频| 国产一级毛片在线| 美女脱内裤让男人舔精品视频| 国产亚洲午夜精品一区二区久久 | 国产欧美另类精品又又久久亚洲欧美| 午夜老司机福利剧场| 中文资源天堂在线| 日韩欧美精品v在线| 在现免费观看毛片| 波多野结衣巨乳人妻| 国产精品,欧美在线| 精品欧美国产一区二区三| 日本熟妇午夜| 欧美bdsm另类| 亚洲成人久久爱视频| 精品久久久久久久人妻蜜臀av| 菩萨蛮人人尽说江南好唐韦庄 | 一级毛片aaaaaa免费看小| 日韩av不卡免费在线播放| 一本久久精品| 黄色欧美视频在线观看| 国内揄拍国产精品人妻在线| 亚洲最大成人手机在线| 国产毛片a区久久久久| 国产精品久久久久久精品电影| 免费看av在线观看网站| 国产精品.久久久| av在线天堂中文字幕| 亚洲欧美精品综合久久99| 国产三级中文精品| 丝袜美腿在线中文| 91精品一卡2卡3卡4卡| 男人舔女人下体高潮全视频| 韩国av在线不卡| 日本熟妇午夜| 成人漫画全彩无遮挡| 中文字幕制服av| 亚洲激情五月婷婷啪啪| 99九九线精品视频在线观看视频| 国产日韩欧美在线精品| av线在线观看网站| 在线观看66精品国产| 国产高清国产精品国产三级 | 天堂√8在线中文| 精品久久久久久久久久久久久| 中文欧美无线码| 日本av手机在线免费观看| 亚洲三级黄色毛片| 日日撸夜夜添| 少妇熟女欧美另类| 麻豆久久精品国产亚洲av| 亚洲五月天丁香| 国产一区有黄有色的免费视频 | 成人三级黄色视频| 国产精品久久久久久av不卡| 亚洲欧美成人精品一区二区| 欧美日韩综合久久久久久| 亚洲精品自拍成人| 九色成人免费人妻av| 我要看日韩黄色一级片| 亚洲美女搞黄在线观看| 欧美日韩在线观看h| 一二三四中文在线观看免费高清| 日韩成人伦理影院| 少妇丰满av| 韩国高清视频一区二区三区| 亚洲成色77777| 天美传媒精品一区二区| 精品久久久久久久久久久久久| eeuss影院久久| 欧美成人a在线观看| 亚洲丝袜综合中文字幕| 黑人高潮一二区| 国产片特级美女逼逼视频| 小说图片视频综合网站| av又黄又爽大尺度在线免费看 | 可以在线观看毛片的网站| 99久久中文字幕三级久久日本| 国模一区二区三区四区视频| 精品久久久久久成人av| 亚洲国产精品合色在线| 国产成人免费观看mmmm| 禁无遮挡网站| 亚洲天堂国产精品一区在线| 久久久午夜欧美精品| 青青草视频在线视频观看| 亚洲av中文字字幕乱码综合| 日本五十路高清| 大香蕉97超碰在线| 精品熟女少妇av免费看| 中国国产av一级| 国产伦精品一区二区三区四那| 亚洲精品自拍成人| 最近的中文字幕免费完整| 三级经典国产精品| av视频在线观看入口| 国产乱人偷精品视频| 亚洲欧洲国产日韩| 一边摸一边抽搐一进一小说| 国产一区有黄有色的免费视频 | 久久久久国产网址| 一级毛片久久久久久久久女| 99热这里只有是精品50| 日韩大片免费观看网站 | 亚洲在久久综合| 99久久中文字幕三级久久日本| 久久精品国产鲁丝片午夜精品| 99久久九九国产精品国产免费| av视频在线观看入口| 熟妇人妻久久中文字幕3abv| 亚洲av福利一区| 天堂网av新在线| 亚洲自拍偷在线| 欧美3d第一页| 国产极品精品免费视频能看的| 美女cb高潮喷水在线观看| 国产免费男女视频| 日本猛色少妇xxxxx猛交久久| 老司机福利观看| ponron亚洲| 欧美色视频一区免费| 精品熟女少妇av免费看| 亚洲av熟女| 国产精品一区二区三区四区免费观看| 一夜夜www| 免费人成在线观看视频色| 亚洲国产欧美在线一区| 别揉我奶头 嗯啊视频| 午夜a级毛片| 国产老妇伦熟女老妇高清| 精品一区二区三区人妻视频| 一边亲一边摸免费视频| 久久久a久久爽久久v久久| 国产成人免费观看mmmm| 亚洲av不卡在线观看| 少妇人妻一区二区三区视频| 毛片女人毛片| 亚洲欧美精品自产自拍| 美女xxoo啪啪120秒动态图| 高清视频免费观看一区二区 | 69人妻影院| 老女人水多毛片| 国产欧美另类精品又又久久亚洲欧美| 少妇熟女欧美另类| 久久精品91蜜桃| 亚洲欧美日韩高清专用| 亚洲av日韩在线播放| 国产免费视频播放在线视频 | 国产毛片a区久久久久| 嫩草影院新地址| 国产淫语在线视频| 91久久精品国产一区二区三区| 在线播放国产精品三级| 久久久成人免费电影| 在线免费十八禁| 亚洲激情五月婷婷啪啪| 国产69精品久久久久777片| av天堂中文字幕网| 午夜免费激情av| 国内精品美女久久久久久| 女人十人毛片免费观看3o分钟| 人妻制服诱惑在线中文字幕| 噜噜噜噜噜久久久久久91| 91精品一卡2卡3卡4卡| 色吧在线观看| 久久久精品大字幕| 高清av免费在线| 观看免费一级毛片| 有码 亚洲区| 蜜桃亚洲精品一区二区三区| 婷婷色麻豆天堂久久 | 中文亚洲av片在线观看爽| 一边摸一边抽搐一进一小说| 久久精品国产亚洲av涩爱| 午夜亚洲福利在线播放| a级一级毛片免费在线观看| 久久久精品欧美日韩精品| 老女人水多毛片| 欧美精品国产亚洲| 日韩一区二区视频免费看| 美女黄网站色视频| 国产欧美日韩精品一区二区| 国产伦精品一区二区三区视频9| 国产 一区 欧美 日韩| ponron亚洲| 国产真实乱freesex| 日日啪夜夜撸| 国产片特级美女逼逼视频| 亚洲精品日韩在线中文字幕| 男插女下体视频免费在线播放| 1000部很黄的大片| 少妇熟女欧美另类| 寂寞人妻少妇视频99o| 在线免费十八禁| 插逼视频在线观看| 国产片特级美女逼逼视频| 日韩人妻高清精品专区| 好男人在线观看高清免费视频| 成人毛片60女人毛片免费| 插阴视频在线观看视频| 永久免费av网站大全| 国产成人freesex在线| 日本av手机在线免费观看| 欧美最新免费一区二区三区| 久久国内精品自在自线图片| 国内精品一区二区在线观看| 日本一二三区视频观看| 日韩人妻高清精品专区| 国产精品久久久久久精品电影小说 | 亚洲中文字幕一区二区三区有码在线看| 日本av手机在线免费观看| 人人妻人人澡人人爽人人夜夜 | 国产成人精品婷婷| 国产成人午夜福利电影在线观看| 综合色av麻豆| 看黄色毛片网站| 狠狠狠狠99中文字幕| 亚洲欧美成人精品一区二区| 最近中文字幕2019免费版| 黄片wwwwww| 成年免费大片在线观看| 久久精品国产自在天天线| 中文字幕精品亚洲无线码一区| 欧美日韩国产亚洲二区| 久久韩国三级中文字幕| 欧美人与善性xxx| 毛片一级片免费看久久久久| 久久鲁丝午夜福利片| 国产一区二区亚洲精品在线观看| 联通29元200g的流量卡| 婷婷六月久久综合丁香| 亚洲aⅴ乱码一区二区在线播放| 日韩国内少妇激情av| 国产精品久久久久久精品电影| 一本久久精品| 国产午夜福利久久久久久| 嫩草影院入口| av女优亚洲男人天堂| 成人二区视频| 国内揄拍国产精品人妻在线| 婷婷色av中文字幕| 久久久久久伊人网av| 亚洲va在线va天堂va国产| 国产白丝娇喘喷水9色精品| 欧美一区二区精品小视频在线| 美女xxoo啪啪120秒动态图| 一级毛片aaaaaa免费看小| 一个人免费在线观看电影| 内地一区二区视频在线| 91aial.com中文字幕在线观看| 亚洲图色成人| 久久精品久久久久久噜噜老黄 | 亚洲人成网站高清观看| 精品午夜福利在线看| 久久精品夜色国产| 免费播放大片免费观看视频在线观看 | 波多野结衣高清无吗| 麻豆av噜噜一区二区三区| 三级国产精品欧美在线观看| 一本久久精品| 久久久久久久久中文| 成年女人永久免费观看视频| 亚洲无线观看免费| 亚洲最大成人av| 亚洲av不卡在线观看| 亚洲av福利一区| 欧美97在线视频| 国产精品女同一区二区软件| 熟妇人妻久久中文字幕3abv| 毛片女人毛片| 亚洲电影在线观看av| 色视频www国产| 神马国产精品三级电影在线观看| 狂野欧美白嫩少妇大欣赏| 免费播放大片免费观看视频在线观看 | 一个人免费在线观看电影| 日日摸夜夜添夜夜添av毛片| 我要看日韩黄色一级片| 色网站视频免费| 好男人视频免费观看在线| 1000部很黄的大片| av女优亚洲男人天堂| 国内精品一区二区在线观看| 中文乱码字字幕精品一区二区三区 | 久久精品国产亚洲av涩爱| 国产精品爽爽va在线观看网站| 国产淫语在线视频| 日日摸夜夜添夜夜爱| 国产探花极品一区二区| 国产一区亚洲一区在线观看| 联通29元200g的流量卡| 一边亲一边摸免费视频| 伦理电影大哥的女人| 91精品伊人久久大香线蕉| 中文天堂在线官网| 2021少妇久久久久久久久久久| 日韩欧美精品v在线| 欧美成人精品欧美一级黄| 亚洲av中文av极速乱| 亚洲不卡免费看| 亚洲av男天堂| 亚洲欧美成人精品一区二区| 99在线人妻在线中文字幕| 日本熟妇午夜| 99在线视频只有这里精品首页| 亚洲成人av在线免费| a级毛色黄片| 神马国产精品三级电影在线观看| 日韩欧美国产在线观看| 久久婷婷人人爽人人干人人爱| a级一级毛片免费在线观看| 卡戴珊不雅视频在线播放| 久久亚洲国产成人精品v| 亚洲精品国产av成人精品| 精品一区二区三区人妻视频| 久久久久久久国产电影| 国产成人福利小说| 欧美日本视频| 亚洲av免费高清在线观看| 国产真实伦视频高清在线观看| 精品99又大又爽又粗少妇毛片| 久久久精品94久久精品| 九九在线视频观看精品| 日韩亚洲欧美综合| 午夜a级毛片| a级一级毛片免费在线观看| 长腿黑丝高跟| 国产毛片a区久久久久| 三级毛片av免费| 日韩精品青青久久久久久| 国产亚洲av片在线观看秒播厂 | 亚洲第一区二区三区不卡| 免费看日本二区| 精品久久国产蜜桃| 国产成人午夜福利电影在线观看| 国产在视频线精品| 国产精品福利在线免费观看| 欧美极品一区二区三区四区| 国产视频首页在线观看| 精品人妻一区二区三区麻豆| 美女xxoo啪啪120秒动态图| 亚洲av日韩在线播放| 免费黄色在线免费观看| 七月丁香在线播放| 尾随美女入室| 国产精品久久久久久av不卡| 日本与韩国留学比较| 男女那种视频在线观看| 国产av不卡久久| 看十八女毛片水多多多| 麻豆成人午夜福利视频| 国产午夜精品久久久久久一区二区三区| 精品国内亚洲2022精品成人| 久久久久久久午夜电影| 三级经典国产精品| 亚洲一区高清亚洲精品| 麻豆国产97在线/欧美| 亚洲精华国产精华液的使用体验| 最近手机中文字幕大全| 久久精品国产自在天天线| 在线播放国产精品三级| 日本av手机在线免费观看| 精品国产一区二区三区久久久樱花 | 两性午夜刺激爽爽歪歪视频在线观看| 免费av毛片视频| 亚洲欧美清纯卡通| 少妇熟女aⅴ在线视频| or卡值多少钱| 国产精品久久久久久久久免|