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

    高超聲速壁湍流入口條件生成方法的比較

    2017-12-25 03:23:31袁湘江朱志斌
    空氣動力學(xué)學(xué)報 2017年6期
    關(guān)鍵詞:馬赫數(shù)波紋湍流

    禹 旻,袁湘江,朱志斌

    (中國航天空氣動力技術(shù)研究院,北京 100074)

    高超聲速壁湍流入口條件生成方法的比較

    禹 旻*,袁湘江,朱志斌

    (中國航天空氣動力技術(shù)研究院,北京 100074)

    給定恰當(dāng)?shù)娜肟跅l件是開展壁湍流數(shù)值模擬的關(guān)鍵問題。選取基于有限差分的直接數(shù)值模擬方法結(jié)合高精度的計算格式,詳細(xì)討論了自然轉(zhuǎn)捩、波紋壁面結(jié)構(gòu)導(dǎo)致的“Bypass”轉(zhuǎn)捩和利用時間發(fā)展湍流場進(jìn)行參數(shù)回收促發(fā)轉(zhuǎn)捩這幾種湍流入口生成方法在高超聲速條件下的可行性,分析了這幾種方法各自存在的優(yōu)點和不足之處。計算結(jié)果表明,Bypass轉(zhuǎn)捩和利用時間發(fā)展流場進(jìn)行參數(shù)回收方法相比自然轉(zhuǎn)捩能快速促發(fā)轉(zhuǎn)捩,但自然轉(zhuǎn)捩得到的湍流場品質(zhì)更好。該項研究為高馬赫數(shù)壁湍流數(shù)值模擬入口條件的選取提供參考依據(jù)。

    高超聲速壁湍流;入口條件;自然轉(zhuǎn)捩;波紋壁面;流場參數(shù)回收

    0 引 言

    湍流的多樣性使得描述一個湍流流動入口是非常困難的。對于湍流的數(shù)值模擬,通常需要在入口給定復(fù)雜的湍流脈動生成方法來模擬真實的湍流場。常用的湍流入口邊界生成方法有自然轉(zhuǎn)捩方法、合成湍流方法和流場參數(shù)回收方法等等[1-3]。

    自然轉(zhuǎn)捩方法最為簡單[4],但計算量大,計算域包含了層流區(qū)、轉(zhuǎn)捩區(qū)和湍流區(qū)。流場參數(shù)回收方法是采用計算本身來產(chǎn)生入口所需的湍流脈動信息,它其中又包括利用時間發(fā)展湍流場方法(流向周期性方法)[5-7]和回收/調(diào)節(jié)方法[8-9],該方法適合的研究對象是不關(guān)心轉(zhuǎn)捩區(qū)的充分發(fā)展湍流。在工程上,經(jīng)常使用在壁面上引入粗糙度、凸起或壁面變形等使流動和幾何參數(shù)發(fā)生改變的方法,對轉(zhuǎn)捩產(chǎn)生影響[10-11]。

    由于壓縮性的影響,高馬赫數(shù)下邊界層流動更為穩(wěn)定,因此在高超聲速下生成湍流場就更為困難。一些常用的生成湍流場的方法在高馬赫數(shù)條件下不能達(dá)到預(yù)期效果。比如回收/調(diào)節(jié)方法,該方法把下游某處速度場、熱力學(xué)參量和各種湍流統(tǒng)計量按照一定的比例縮放后引入入口,但這樣引入的前提是Morkovin假說成立[9],其實該方法在M>5時就不是很有效了。在實際數(shù)值計算中也發(fā)現(xiàn),在低馬赫數(shù)情況下(M=3)使用該方法能產(chǎn)生湍流場,而在高馬赫數(shù)情況下(M=6)就不能激發(fā)出湍流場。

    本文主要針對高馬赫數(shù)流動,對自然轉(zhuǎn)捩、波紋壁面結(jié)構(gòu)和利用時間發(fā)展流場的參數(shù)回收這幾種方法進(jìn)行了數(shù)值計算,并對結(jié)果進(jìn)行了分析和評價。

    1 控制方程和數(shù)值方法

    考慮一般曲線坐標(biāo)系(ξ,η,ζ),在此坐標(biāo)系下三維Navier-Stokes方程形式為:

    采用完全氣體假設(shè),方程的具體表達(dá)式參見文獻(xiàn)[12]。方程(1)的數(shù)值求解采用有限差分方法,通量采用Sterge-Warming矢通量分裂形式,通量離散波紋壁計算時使用五階WENO格式[13],其他算例使用五階迎風(fēng)格式;黏性項采用六階中心差分格式;時間方向采用具有TVD性質(zhì)的三步三階Runge-Kutta方法[14]推進(jìn)求解。

    邊界處理:壁面采用無滑移條件和等溫壁條件;出口采用外推方法;展向采用周期邊界條件。

    2 幾種可壓縮湍流入口條件生成方法

    氣體參數(shù)選取26 km高空處值,來流馬赫數(shù)為6,壁溫為600 K。計算時高空來流速度、密度和溫度值作為相應(yīng)無量綱參考量,長度用L=1 m作為無量綱參考量。

    2.1 自然轉(zhuǎn)捩方法

    在計算域入口存在邊界層問題時,將計算域入口取在上游層流段,根據(jù)Orr-Sommerfeld方程得到不穩(wěn)定特征模態(tài),疊加至層流剖面上,激勵擾動經(jīng)過線性、非線性增長轉(zhuǎn)捩至成為完全發(fā)展的湍流。

    在高馬赫數(shù)情況下,最不穩(wěn)定的是第二模態(tài)二維擾動[15]。為促使擾動快速進(jìn)入非線性增長階段,在計算域入口處加入一個第二模態(tài)二維擾動和一對第一模態(tài)三維擾動,擾動的形式為:

    表1 入口擾動參數(shù)(無量綱量,自然轉(zhuǎn)捩方法)Table 1 Parameters of the disturbances at the inlet (dimensionless quantities,natural transition)

    圖1中黑色實線給出了自然轉(zhuǎn)捩方法得到的物面摩擦系數(shù)分布,虛線是層流和湍流壁面摩擦系數(shù)的理論值[16]。在x=0.55 m左右曲線開始抬升,轉(zhuǎn)捩開始,經(jīng)過一段距離后達(dá)到湍流狀態(tài)。圖2給出了計算到統(tǒng)計定常后速度梯度張量二階不變量Q2=5000的三維等值面,可以看到,計算域入口附近流動是以二維結(jié)構(gòu)為主,隨著向下游的推進(jìn),開始有三維結(jié)構(gòu)出現(xiàn),Λ渦形成,在流向x>0.9 m到達(dá)完全湍流區(qū)。圖3給出了流向為x=0.9 m處的平均速度剖面沿壁面法向的分布,實線為平均速度,三角形為Van Driest變換后的平均速度,虛線分別為壁面律(u+=y+)和對數(shù)律(u+=2.5ln(y+)+5.8),平均速度剖面與對數(shù)率和壁面率符合得很好。在計算域入口處添加上述擾動,數(shù)值模擬獲得了平板邊界層從層流失穩(wěn)轉(zhuǎn)捩到湍流的整個發(fā)展過程。

    2.2 壁面上引入波紋結(jié)構(gòu)

    在結(jié)構(gòu)表面上引入粗糙單元或改變壁面局部結(jié)構(gòu),都可使局部流場狀態(tài)發(fā)生改變,促使流向渦的產(chǎn)生,使流動剪切作用增強,起到促進(jìn)轉(zhuǎn)捩的效果。

    在流向位置為x1=0.45 m至x2=0.5 m之間布置波紋壁面,波紋壁的分布函數(shù)為:

    其中:高度h=2.45 mm(約為0.7倍當(dāng)?shù)剡吔鐚雍穸?,kx=5×2π/(x2-x1),kz=2×2π/lz,lz為展向長度。

    圖4給出了Q2=5000的三維等值面,在下游由于波紋壁面的干擾作用,使得流動剪切作用增強,流場變得不穩(wěn)定,層流快速轉(zhuǎn)捩到湍流,從圖4中可以看出x≈0.6 m左右開始流動結(jié)構(gòu)就比較飽滿了。圖1中實心方形表示的是波紋壁方法得到的物面摩擦系數(shù)曲線,由于波紋壁的存在,能引發(fā)壁面摩擦系數(shù)曲線迅速抬升,在經(jīng)過波紋壁面之后,曲線會緩慢下降,經(jīng)過一段距離后達(dá)到完全湍流狀態(tài)。圖5給出了x=0.68 m處的平均速度剖面,可以看到平均速度剖面同對數(shù)率和壁面率符合得較好。

    2.3 利用時間發(fā)展湍流場進(jìn)行參數(shù)回收方法

    采用時間發(fā)展湍流場作為空間發(fā)展流場的入口條件,預(yù)先模擬的時間發(fā)展湍流場流向采用周期性邊界條件,在進(jìn)行回收時只需用預(yù)先模擬的湍流場的某一時刻的流動信息即可,預(yù)存的數(shù)據(jù)量小。

    首先,計算一個隨時間發(fā)展的湍流場。與2.1節(jié)計算自然轉(zhuǎn)捩時引入擾動的方法類似,初始時刻在整個計算域中引入一些擾動,擾動形式為:

    具體參數(shù)見表2,此處ωi代表擾動的增長率。

    表2 擾動的具體參數(shù)(無量綱量,時間發(fā)展的DNS)Table 2 Parameters of the disturbances at the inlet (dimensionless quantities,temporal DNS)

    圖6給出了壁面摩擦系數(shù)曲線隨時間的分布??梢钥吹絫=0.2 s左右壁面摩擦系數(shù)曲線開始抬升,也就意味著轉(zhuǎn)捩的開始,到t=0.5 s左右達(dá)到峰值進(jìn)入湍流區(qū),到達(dá)湍流區(qū)后系數(shù)曲線的變化趨于平緩且有所下降。圖7為0.8 s時速度梯度張量二階不變量Q2=5000的三維等值面,此時已是充分發(fā)展湍流。取0.8 s時刻的瞬時流場數(shù)據(jù)作為空間發(fā)展湍流流場的入口條件,假定此瞬時湍流流場數(shù)據(jù)為UT(x,y,z),則空間模式DNS入口條件可表示為:

    這里通過時空變換關(guān)系,x=ct,將某一時刻隨時間發(fā)展的湍流場整個計算域的流場信息轉(zhuǎn)換為空間發(fā)展的湍流流場入口處的時間序列。黃章峰等[5]認(rèn)為,壁湍流是具有大尺度相干結(jié)構(gòu)的,而相干結(jié)構(gòu)對湍流的性質(zhì)起著決定性的作用,而相干結(jié)構(gòu)具有某些波的特點,其平均傳播速度約為0.92,在實際計算中也發(fā)現(xiàn)設(shè)定為此對流速度能得到較好的效果,因此取傳播速度c為0.9。

    圖1中的空心圓形為回收方法得到的壁面摩擦系數(shù)沿流向的分布,可看到由于入口條件的人為給定,流動需要經(jīng)過一小段調(diào)整,之后就維持在完全湍流狀態(tài)。圖8為Q2=5000的三維等值面,計算域入口開始就是湍流狀態(tài),而且在向下游的推進(jìn)中能維持住湍流流動。圖9給出了流向為0.55 m處平均速度剖面,同對數(shù)率和壁面率符合得較好。

    3 可壓縮湍流入口條件生成方法評價

    對第2節(jié)中計算出的湍流場進(jìn)行歸納總結(jié),從生成湍流場品質(zhì)和能否快速生成湍流場兩方面來對這幾類入口條件給定方法進(jìn)行評價。首先,圖10分別給出了y+=40時3種方法生成的湍流場速度分布云圖,自然轉(zhuǎn)捩方法生成的湍流場最為均勻飽滿;引入波紋壁面結(jié)構(gòu)促發(fā)轉(zhuǎn)捩生成的湍流場不如自然轉(zhuǎn)捩方法得到的湍流場均勻、飽滿;同樣,使用隨時間發(fā)展湍流場進(jìn)行參數(shù)回收方法時得到的湍流場也不如自然轉(zhuǎn)捩方法得到的湍流場飽滿。

    另從圖1中的壁面摩擦系數(shù)曲線對比圖,可以直觀地看到引入波紋壁面結(jié)構(gòu)方法相比自然轉(zhuǎn)捩方法能更為快速的生成湍流場,不需要經(jīng)歷很長的流向距離(從流動入口位置0.4 m到完全湍流處0.7 m左右)。而采用時間發(fā)展流場進(jìn)行回收作為入口條件的方法,從計算域入口開始就是湍流場,只要經(jīng)過一小段的調(diào)整(入口0.4 m到0.45 m左右)就能獲得一個比較真實的充分發(fā)展湍流,此種方法也是一種快速生成湍流場的方法。

    采用時間發(fā)展流場進(jìn)行回收方法可以推廣到其它外形和不同計算工況。董明等[6]證實此方法是對不同馬赫數(shù)、雷諾數(shù)、壁面條件和典型邊界層類型均可行的有效的空間模式直接數(shù)值模擬入口條件生成方法,在實際數(shù)值計算中也發(fā)現(xiàn),使用文獻(xiàn)中對平均密度和溫度以及脈動量進(jìn)行相應(yīng)修正的方法,可導(dǎo)出不同計算條件下的入口條件。在這里用2.3節(jié)計算得到的隨時間發(fā)展的湍流場作為入口,可以計算出來流馬赫數(shù)7、高空27 km、壁面溫度為700 K時的湍流場。圖11為Q2=5000的三維等值面,可以看到同樣從計算域入口開始就是湍流狀態(tài)。圖12為壁面摩擦系數(shù)沿流向的分布,可看到由于入口條件的影響,流動需要經(jīng)過一段比之前相同計算工況條件下稍長一些的調(diào)整距離之后就維持在完全湍流狀態(tài)??梢?,采用隨時間發(fā)展湍流場進(jìn)行參數(shù)回收方法,可以推廣到更高的馬赫數(shù)。

    4 結(jié) 論

    數(shù)值模擬研究了自然轉(zhuǎn)捩、波紋壁誘發(fā)轉(zhuǎn)捩、利用時間發(fā)展邊界層湍流這三類促發(fā)湍流場生成的入口條件給定方法。研究表明,在高馬赫流動下這三種方法都可以實現(xiàn)湍流場的生成。其中自然轉(zhuǎn)捩方法處理起來最為簡單,得到的湍流場也最為均勻飽滿,但是計算域長,網(wǎng)格量大,不能快速生成湍流場;波紋壁面的引入,能很快誘導(dǎo)湍流生成,是一種能快速生成湍流的方法,但生成的湍流場不如自然轉(zhuǎn)捩方法得到的湍流場均勻飽滿;利用時間發(fā)展流場進(jìn)行參數(shù)回收這種湍流生成方法,從計算域入口就是湍流場,也是一種能有效生成湍流場的方法,該方法還有一個好處就是可由同一預(yù)先模擬的時間發(fā)展湍流場導(dǎo)出不同馬赫數(shù)、雷諾數(shù)、壁面條件所需要的湍流入口條件,這對于要計算很多參數(shù)下的湍流情況是很有用的。

    [1]Lund T S,Wu X H,Squires K D.Generation of turbulent in flow data for spatially-developing boundary layer simulations[J].Journal of Computational Physics,1998,140(2):233-258.

    [2]Sivasubramanian J,Fasel H F.Direct numerical simulation of transition in a sharp cone boundary layer at Mach 6:Fundamental breakdown[J].Journal of Fluid Mechanics,2015,768:175-218.

    [3]Xu S,Martin M P.Assessment of inflow boundary conditions for compressible turbulent boundary layers[J].Physics of Fluids,2004,16(7):2623-2639.

    [4]Batten P,Goldberg U,Chakravarthy S.Interfacing statistical turbulence closures with large-eddy simulation[J].AIAA Journal,2004,42(3):485-492.

    [5]黃章峰,周恒.湍流邊界層的空間模式DNS的入口條件[J].中國科學(xué)(G輯:物理學(xué) 力學(xué) 天文學(xué)),2008,38(3):310-318.

    [6]董明,周恒.超聲速鈍錐湍流邊界層DNS入口邊界條件的研究[J].應(yīng)用數(shù)學(xué)和力學(xué),2008,29(8):893-904.

    [7]逯學(xué)志,黃章峰,羅紀(jì)生.湍流入流條件在強迫轉(zhuǎn)捩研究中的應(yīng)用[J].空氣動力學(xué)學(xué)報,2012,30(6):749-753.

    [8]陳逖,劉衛(wèi)東,范曉檣,等.“回收/調(diào)節(jié)”方法在混合LES/RANS模擬方法中的應(yīng)用[J].航空動力學(xué)報,2011,26(6):1215-1222.

    [9]Xiao X D,Edwards J R,Hassan H A.inflow boundary conditions for hybrid large eddy/Reynolds averaged Navier-Stokes simulations[J].AIAA Journal,2003,41(8):1481-1489.

    [10]Schneider S P.Effects of roughness on hypersonic Boundary-Layer transition[R].AIAA 2007-305,2007.

    [11]Berry S A.Discrete roughness transition for hypersonic flight vehicles[R].AIAA 2007-307,2007.

    [12]劉建新.小攻角鈍錐高超聲速邊界層的擾動演化[D].天津:天津大學(xué),2010.

    [13]Jiang G S,Shu C W.Efficient implementation of weighted ENO Schemes[J].Journal of Computational Physics,1996,126(1):202-228.

    [14]Gottlieb S,Shu C W.Total variation diminishing Runge-Kutta schemes[J].Mathematics of computation,1998,67(221):73-85.

    [15]Yu M,Luo J S.Nonlinear interaction mechanisms of disturbances in supersonic flat-plate boundary layers[J].Science China:Physics Mechanics and Astronomy,2014,57(11):2141-2151.

    [16] White F M.Viscous fluid flow[M].McGraw-Hill Education,2005.

    Assessmentofinflowboundaryconditionsforhypersonicwallboundedturbulentflows

    YU Min*,YUAN Xiangjiang,ZHU Zhibin

    (ChinaAcademyofAerospaceAerodynamics,Beijing100074,China)

    It is critically important to generate turbulent fluctuating inflow boundary conditions for numerical simulation in solving hypersonic turbulent boundary layers.A high-order finite-difference method was applied to investigate the inflow generation methods.Turbulent inflow boundary conditions are broadly classified into natural transition methods,recycling/rescaling-based methods,and synthetic turbulence generators.To investigate the feasibility of inflow generation methods at high Mach number,the simulation of a Mach 6 boundary layer over an isothermal flat plate was carried out.The validity and limitations of different inflow boundary conditions were also discussed.Several feasible methods were studied in great detail including natural transition from laminar to turbulent flows,bypass transition induced by wavy wall configuration,and turbulence generation via reintroducing the flow field obtained from temporal developing turbulent flow.The limitation of the recycling/rescaling method was also proposed in hypersonic flow condition.By analyzing the flow filed and statistical results,we summarized and discussed the advantages and drawbacks about these inflow generation approaches,and further assessed these methods.Compared with natural transition method,bypass transition method and reintroducing method can promote transition to turbulent flow field efficiently.But turbulent field has better quality generated by natural transition route.The paper can be a reference of imposing an inflow boundary condition for wall bounded turbulent flows at high Mach number.

    hypersonic wall bounded turbulent flows; inflow boundary condition; natural transition; wavy wall; reintroducing methods

    0258-1825(2017)06-0772-05

    O357.5

    A

    10.7638/kqdlxxb-2015.0177

    2015-09-10;

    2015-12-15

    國家自然科學(xué)基金(11772316)

    禹旻*(1986-),女,天津人,高級工程師,研究方向:湍流直接數(shù)值模擬,穩(wěn)定性分析.E-mail:yumin8688@126.com

    禹旻,袁湘江,朱志斌.高超聲速壁湍流入口條件生成方法的比較[J].空氣動力學(xué)學(xué)報,2017,35(6):772-776.

    10.7638/kqdlxxb-2015.0177 YU M,YUAN X J,ZHU Z B.Assessment of inflow boundary conditions for hypersonic wall bounded turbulent flows[J].Acta Aerodynamica Sinica,2017,35(6):772-776.

    猜你喜歡
    馬赫數(shù)波紋湍流
    高馬赫數(shù)激波作用下單模界面的Richtmyer-Meshkov不穩(wěn)定性數(shù)值模擬
    爆炸與沖擊(2024年7期)2024-11-01 00:00:00
    一維非等熵可壓縮微極流體的低馬赫數(shù)極限
    基于NACA0030的波紋狀翼型氣動特性探索
    小波紋的童話
    載荷分布對可控擴(kuò)散葉型性能的影響
    重氣瞬時泄漏擴(kuò)散的湍流模型驗證
    為什么水面波紋蕩漾
    學(xué)與玩(2017年5期)2017-02-16 07:06:26
    “青春期”湍流中的智慧引渡(三)
    “青春期”湍流中的智慧引渡(二)
    弱分層湍流輸運特性的統(tǒng)計分析
    99久久中文字幕三级久久日本| 久久午夜亚洲精品久久| 午夜老司机福利剧场| 美女内射精品一级片tv| 久久99热6这里只有精品| 久久久久久久久久久丰满| 欧美日本视频| 日韩精品青青久久久久久| 99久久中文字幕三级久久日本| 成人午夜高清在线视频| 国产真实乱freesex| 老女人水多毛片| 淫妇啪啪啪对白视频| 在线国产一区二区在线| 亚洲精品日韩在线中文字幕 | 美女大奶头视频| 日本 av在线| 我要搜黄色片| 免费看av在线观看网站| 色5月婷婷丁香| 天天躁夜夜躁狠狠久久av| 国产高清视频在线观看网站| 亚洲av熟女| 成年免费大片在线观看| 国产成人一区二区在线| 亚洲激情五月婷婷啪啪| 丰满乱子伦码专区| 99热6这里只有精品| 精品久久久噜噜| 三级经典国产精品| 久久久成人免费电影| 久久综合国产亚洲精品| 久99久视频精品免费| 99热这里只有是精品50| 成人av一区二区三区在线看| 亚洲av成人av| 久久国内精品自在自线图片| 99热这里只有是精品50| 91久久精品电影网| 国国产精品蜜臀av免费| 非洲黑人性xxxx精品又粗又长| 99热这里只有是精品50| 国产精品野战在线观看| 一个人看的www免费观看视频| 3wmmmm亚洲av在线观看| 少妇猛男粗大的猛烈进出视频 | 少妇熟女aⅴ在线视频| 亚洲精品在线观看二区| 少妇裸体淫交视频免费看高清| 亚洲成a人片在线一区二区| 小蜜桃在线观看免费完整版高清| 自拍偷自拍亚洲精品老妇| 一级毛片电影观看 | 亚洲国产欧美人成| 国产麻豆成人av免费视频| 99精品在免费线老司机午夜| 人妻夜夜爽99麻豆av| 97超视频在线观看视频| 国产欧美日韩一区二区精品| 亚洲中文日韩欧美视频| 日本与韩国留学比较| a级毛片免费高清观看在线播放| 国产亚洲av嫩草精品影院| 国产av一区在线观看免费| 欧美不卡视频在线免费观看| 免费av毛片视频| 午夜福利在线观看吧| 蜜桃久久精品国产亚洲av| 久久鲁丝午夜福利片| 国产午夜精品久久久久久一区二区三区 | 国产私拍福利视频在线观看| 性色avwww在线观看| 日韩欧美国产在线观看| 亚洲成人久久性| 天堂动漫精品| 在线免费观看不下载黄p国产| 可以在线观看毛片的网站| 成人av一区二区三区在线看| 成年女人永久免费观看视频| 欧美中文日本在线观看视频| 免费av观看视频| 国产爱豆传媒在线观看| 国产精品精品国产色婷婷| 日本一二三区视频观看| 精品久久久久久久人妻蜜臀av| 精品一区二区三区人妻视频| 久久人人爽人人爽人人片va| 久99久视频精品免费| 午夜视频国产福利| 一级a爱片免费观看的视频| 综合色丁香网| 熟女人妻精品中文字幕| 欧美激情国产日韩精品一区| 国产一区二区在线观看日韩| 日韩欧美免费精品| 亚洲图色成人| 色哟哟哟哟哟哟| 丰满的人妻完整版| 国产亚洲av嫩草精品影院| 黄色欧美视频在线观看| 欧美日韩一区二区视频在线观看视频在线 | 91久久精品电影网| 国国产精品蜜臀av免费| 性欧美人与动物交配| 国产成人一区二区在线| 精品国产三级普通话版| 性色avwww在线观看| 午夜久久久久精精品| 欧美区成人在线视频| 亚洲综合色惰| 三级毛片av免费| 国产在线男女| 3wmmmm亚洲av在线观看| 性色avwww在线观看| 国产成人a区在线观看| 1000部很黄的大片| 日韩欧美免费精品| 久久人妻av系列| 亚洲一区二区三区色噜噜| 久久精品影院6| 日韩 亚洲 欧美在线| 小蜜桃在线观看免费完整版高清| 国产不卡一卡二| 一区二区三区免费毛片| 欧美最黄视频在线播放免费| 精品人妻视频免费看| 69人妻影院| 别揉我奶头 嗯啊视频| 欧美日韩综合久久久久久| 五月玫瑰六月丁香| 欧美3d第一页| 中文亚洲av片在线观看爽| 波多野结衣高清作品| 麻豆成人午夜福利视频| 搞女人的毛片| 国产精品国产三级国产av玫瑰| 麻豆国产av国片精品| 免费一级毛片在线播放高清视频| 色在线成人网| 精品免费久久久久久久清纯| 国产精品永久免费网站| 日韩欧美 国产精品| 国产精品国产三级国产av玫瑰| 国内精品宾馆在线| 淫妇啪啪啪对白视频| 精品一区二区三区视频在线观看免费| 国产男人的电影天堂91| 国产日本99.免费观看| 色5月婷婷丁香| 一级a爱片免费观看的视频| 久久天躁狠狠躁夜夜2o2o| 一本精品99久久精品77| 97碰自拍视频| 天美传媒精品一区二区| 亚洲第一电影网av| 麻豆久久精品国产亚洲av| 国产一区二区激情短视频| av女优亚洲男人天堂| 人人妻人人看人人澡| 精品人妻视频免费看| 日韩av在线大香蕉| 欧美一级a爱片免费观看看| 久久热精品热| 免费搜索国产男女视频| 久久久久久久久中文| 亚洲精品成人久久久久久| av卡一久久| 黄色欧美视频在线观看| 日韩中字成人| 日韩高清综合在线| 91狼人影院| 久久亚洲精品不卡| 少妇人妻精品综合一区二区 | 国产精品亚洲一级av第二区| 精品久久久久久久久亚洲| 国产精品久久久久久av不卡| 成人一区二区视频在线观看| 黄片wwwwww| 午夜老司机福利剧场| 免费人成在线观看视频色| 日本与韩国留学比较| 亚洲高清免费不卡视频| 亚洲中文日韩欧美视频| 久久精品久久久久久噜噜老黄 | 欧美激情在线99| 中文在线观看免费www的网站| 国产女主播在线喷水免费视频网站 | 国产高清三级在线| 国产三级在线视频| 亚洲七黄色美女视频| 午夜久久久久精精品| 中文字幕av成人在线电影| 69人妻影院| 十八禁国产超污无遮挡网站| 亚洲自拍偷在线| 久久午夜福利片| 久久久精品大字幕| 亚洲美女黄片视频| 俺也久久电影网| 日日摸夜夜添夜夜添小说| 日本与韩国留学比较| av卡一久久| 成人精品一区二区免费| av国产免费在线观看| 精品久久久久久久末码| 国产亚洲精品av在线| 22中文网久久字幕| 夜夜爽天天搞| 欧美潮喷喷水| 午夜福利在线在线| 男女视频在线观看网站免费| 亚洲精华国产精华液的使用体验 | 丝袜美腿在线中文| 看片在线看免费视频| 亚洲精品日韩av片在线观看| 国产麻豆成人av免费视频| h日本视频在线播放| 欧美一区二区国产精品久久精品| 国产精品永久免费网站| 又粗又爽又猛毛片免费看| 亚洲av熟女| 国产 一区 欧美 日韩| 蜜臀久久99精品久久宅男| 欧美成人一区二区免费高清观看| 色5月婷婷丁香| 亚洲精品色激情综合| 啦啦啦啦在线视频资源| 身体一侧抽搐| av卡一久久| 免费av观看视频| 美女黄网站色视频| 中文字幕免费在线视频6| 一进一出抽搐gif免费好疼| 久久6这里有精品| 欧美潮喷喷水| 日本精品一区二区三区蜜桃| 色尼玛亚洲综合影院| 99国产极品粉嫩在线观看| 男人的好看免费观看在线视频| 亚洲性久久影院| 久久中文看片网| 最近的中文字幕免费完整| 国产探花极品一区二区| 亚洲国产欧美人成| 国产色婷婷99| 丰满人妻一区二区三区视频av| 啦啦啦啦在线视频资源| 国产一级毛片七仙女欲春2| a级毛片a级免费在线| 51国产日韩欧美| 女人被狂操c到高潮| 色5月婷婷丁香| 一级a爱片免费观看的视频| 丰满乱子伦码专区| 女的被弄到高潮叫床怎么办| 在线a可以看的网站| 中国国产av一级| 成人av在线播放网站| 99久久成人亚洲精品观看| 欧美日韩综合久久久久久| 国产成人91sexporn| 午夜福利视频1000在线观看| 亚洲欧美成人精品一区二区| av国产免费在线观看| 亚洲aⅴ乱码一区二区在线播放| 成人精品一区二区免费| 国产真实伦视频高清在线观看| 日日干狠狠操夜夜爽| 99久国产av精品| 欧美+亚洲+日韩+国产| 天美传媒精品一区二区| 国产在视频线在精品| 国产白丝娇喘喷水9色精品| 九色成人免费人妻av| 51国产日韩欧美| 日韩欧美在线乱码| 一区福利在线观看| 国产美女午夜福利| 国产精品人妻久久久久久| 日韩欧美精品v在线| 久久99热这里只有精品18| 99久久无色码亚洲精品果冻| 男女那种视频在线观看| 欧美不卡视频在线免费观看| 18禁在线播放成人免费| 亚洲国产精品国产精品| 成人av一区二区三区在线看| 午夜精品一区二区三区免费看| 搡老熟女国产l中国老女人| 最新在线观看一区二区三区| 国产色爽女视频免费观看| 亚洲av五月六月丁香网| 久久精品国产亚洲av香蕉五月| 听说在线观看完整版免费高清| 亚洲av美国av| 免费电影在线观看免费观看| 亚洲av第一区精品v没综合| 午夜福利视频1000在线观看| 校园春色视频在线观看| 看黄色毛片网站| av在线蜜桃| 黑人高潮一二区| videossex国产| 亚洲天堂国产精品一区在线| 在线播放国产精品三级| 国语自产精品视频在线第100页| 黑人高潮一二区| 美女内射精品一级片tv| 极品教师在线视频| 亚洲精华国产精华液的使用体验 | 精品久久久久久久久久免费视频| 真人做人爱边吃奶动态| 国产精品久久久久久久电影| 久久99热这里只有精品18| 亚洲高清免费不卡视频| 日本三级黄在线观看| 欧美一级a爱片免费观看看| ponron亚洲| 99热只有精品国产| 美女内射精品一级片tv| 精品久久久噜噜| 成人亚洲欧美一区二区av| 欧美bdsm另类| 午夜精品在线福利| 国产aⅴ精品一区二区三区波| 黄色配什么色好看| av天堂在线播放| 成人午夜高清在线视频| 亚洲av不卡在线观看| 国产色婷婷99| 日本五十路高清| 亚洲综合色惰| 日本免费a在线| 久久九九热精品免费| 亚洲av中文字字幕乱码综合| 男女啪啪激烈高潮av片| 亚州av有码| 99riav亚洲国产免费| 精品99又大又爽又粗少妇毛片| 国产精品不卡视频一区二区| 丰满的人妻完整版| 国产精品一二三区在线看| 两性午夜刺激爽爽歪歪视频在线观看| 国产亚洲91精品色在线| 国产黄色视频一区二区在线观看 | 男女视频在线观看网站免费| 黄色欧美视频在线观看| 无遮挡黄片免费观看| 亚洲欧美清纯卡通| 久久亚洲国产成人精品v| 国产一区二区亚洲精品在线观看| 晚上一个人看的免费电影| 亚洲色图av天堂| 黄片wwwwww| 毛片一级片免费看久久久久| 蜜桃亚洲精品一区二区三区| 日日啪夜夜撸| 一边摸一边抽搐一进一小说| 搡老熟女国产l中国老女人| 精华霜和精华液先用哪个| 日本黄色片子视频| 亚洲电影在线观看av| 亚洲av一区综合| 日日摸夜夜添夜夜添av毛片| 国产成年人精品一区二区| 亚洲人成网站高清观看| 精华霜和精华液先用哪个| 91精品国产九色| 国产午夜福利久久久久久| 俺也久久电影网| 欧美+日韩+精品| av在线天堂中文字幕| 又黄又爽又刺激的免费视频.| 色哟哟哟哟哟哟| 最后的刺客免费高清国语| 国产精品电影一区二区三区| 亚洲国产欧美人成| 成人亚洲精品av一区二区| 日韩精品有码人妻一区| 欧美xxxx性猛交bbbb| 久久久久久九九精品二区国产| 狠狠狠狠99中文字幕| 联通29元200g的流量卡| 日韩欧美一区二区三区在线观看| 亚洲天堂国产精品一区在线| 亚洲欧美日韩无卡精品| 欧美成人免费av一区二区三区| 国产精品一区二区性色av| 国产真实乱freesex| 一个人看的www免费观看视频| 美女内射精品一级片tv| 国产探花极品一区二区| 久久韩国三级中文字幕| 国产精品野战在线观看| 黄色视频,在线免费观看| 精品一区二区三区av网在线观看| 国产国拍精品亚洲av在线观看| 国产v大片淫在线免费观看| 麻豆久久精品国产亚洲av| 69人妻影院| 特级一级黄色大片| 亚洲精品粉嫩美女一区| 久久久a久久爽久久v久久| 国产一区亚洲一区在线观看| 久久草成人影院| 我的老师免费观看完整版| 亚洲成人久久爱视频| 国产精品亚洲一级av第二区| 日日摸夜夜添夜夜爱| 免费av不卡在线播放| 欧美性猛交黑人性爽| 亚洲美女搞黄在线观看 | 欧美丝袜亚洲另类| 亚洲国产精品合色在线| 亚洲av熟女| 天堂网av新在线| eeuss影院久久| 亚洲av电影不卡..在线观看| 亚洲精品在线观看二区| 午夜爱爱视频在线播放| 亚洲精品一区av在线观看| 毛片女人毛片| 少妇熟女aⅴ在线视频| 美女内射精品一级片tv| 一级毛片aaaaaa免费看小| 我的老师免费观看完整版| 深爱激情五月婷婷| 午夜a级毛片| 乱人视频在线观看| 97超级碰碰碰精品色视频在线观看| 国产精品亚洲一级av第二区| 亚洲七黄色美女视频| 国产精品久久久久久久久免| 国产精品1区2区在线观看.| 国产国拍精品亚洲av在线观看| 午夜精品国产一区二区电影 | 国产综合懂色| 99视频精品全部免费 在线| 国产精品国产三级国产av玫瑰| av女优亚洲男人天堂| 国产成人a区在线观看| 亚洲精品亚洲一区二区| 观看美女的网站| 亚洲美女搞黄在线观看 | 噜噜噜噜噜久久久久久91| 亚洲欧美日韩高清专用| 亚洲人成网站高清观看| 婷婷亚洲欧美| 免费在线观看影片大全网站| 日本精品一区二区三区蜜桃| 国产麻豆成人av免费视频| 国产精品乱码一区二三区的特点| 日韩欧美一区二区三区在线观看| 精品欧美国产一区二区三| 亚洲熟妇熟女久久| 久久久色成人| 国产大屁股一区二区在线视频| 欧美日韩精品成人综合77777| 午夜免费激情av| 亚洲丝袜综合中文字幕| 少妇被粗大猛烈的视频| а√天堂www在线а√下载| 国产视频内射| 精品福利观看| 免费高清视频大片| 国产69精品久久久久777片| 亚洲专区国产一区二区| 久久久精品94久久精品| 亚洲四区av| 欧美三级亚洲精品| av卡一久久| 国产伦一二天堂av在线观看| 久久久久久大精品| 淫妇啪啪啪对白视频| 国产一区亚洲一区在线观看| 黑人高潮一二区| 成人二区视频| 久久久久免费精品人妻一区二区| 最新中文字幕久久久久| 女的被弄到高潮叫床怎么办| 亚洲18禁久久av| 久久精品国产亚洲av天美| 亚州av有码| 中文字幕人妻熟人妻熟丝袜美| 女生性感内裤真人,穿戴方法视频| 国产91av在线免费观看| 亚洲自偷自拍三级| av黄色大香蕉| 色综合站精品国产| 久久人妻av系列| or卡值多少钱| 久久人妻av系列| 激情 狠狠 欧美| 国语自产精品视频在线第100页| 午夜亚洲福利在线播放| 高清午夜精品一区二区三区 | av卡一久久| 色哟哟哟哟哟哟| 国模一区二区三区四区视频| 男人舔奶头视频| 在线免费十八禁| 欧美日韩在线观看h| 中文字幕精品亚洲无线码一区| 精华霜和精华液先用哪个| 又爽又黄无遮挡网站| 免费电影在线观看免费观看| av在线老鸭窝| 俄罗斯特黄特色一大片| 日韩欧美在线乱码| 男女边吃奶边做爰视频| 人妻制服诱惑在线中文字幕| 天美传媒精品一区二区| 亚洲第一区二区三区不卡| 午夜福利视频1000在线观看| 一级黄色大片毛片| 国产黄色视频一区二区在线观看 | av视频在线观看入口| 亚洲,欧美,日韩| 久久精品影院6| 啦啦啦观看免费观看视频高清| 精品乱码久久久久久99久播| 日日摸夜夜添夜夜添小说| 亚洲精品粉嫩美女一区| 国产激情偷乱视频一区二区| 亚洲国产精品成人综合色| 欧美激情国产日韩精品一区| 淫秽高清视频在线观看| 精品不卡国产一区二区三区| 久久午夜福利片| 波多野结衣高清无吗| 简卡轻食公司| 日韩一本色道免费dvd| 成年av动漫网址| 欧美日韩精品成人综合77777| 国产亚洲欧美98| 成人三级黄色视频| 午夜激情福利司机影院| 亚洲av.av天堂| 国产又黄又爽又无遮挡在线| 乱人视频在线观看| 国产精品亚洲美女久久久| 99久久精品国产国产毛片| 99视频精品全部免费 在线| 日本撒尿小便嘘嘘汇集6| 九九久久精品国产亚洲av麻豆| 久久久久久久久久久丰满| 中文字幕熟女人妻在线| 亚洲,欧美,日韩| 一本一本综合久久| 九九在线视频观看精品| a级毛色黄片| 又粗又爽又猛毛片免费看| 亚洲第一电影网av| 国产蜜桃级精品一区二区三区| 国产成人aa在线观看| 久久精品国产亚洲av香蕉五月| 亚洲久久久久久中文字幕| 色综合站精品国产| 蜜桃久久精品国产亚洲av| 国产三级中文精品| 三级毛片av免费| 一个人看的www免费观看视频| 午夜免费激情av| 少妇人妻一区二区三区视频| 免费看av在线观看网站| 天堂√8在线中文| 丝袜喷水一区| 午夜亚洲福利在线播放| 亚洲国产高清在线一区二区三| 成人av在线播放网站| 国产精品av视频在线免费观看| 高清毛片免费看| 嫩草影院精品99| 美女被艹到高潮喷水动态| 亚洲婷婷狠狠爱综合网| 亚洲av中文av极速乱| 麻豆乱淫一区二区| 亚洲自偷自拍三级| 欧美最黄视频在线播放免费| 性插视频无遮挡在线免费观看| 国产三级中文精品| 91精品国产九色| 97在线视频观看| 中文字幕av在线有码专区| 久久午夜福利片| 日韩欧美精品v在线| 精品99又大又爽又粗少妇毛片| a级毛片a级免费在线| 国产精品99久久久久久久久| 亚洲精品久久国产高清桃花| 国产精品福利在线免费观看| 国产精品伦人一区二区| 国产av不卡久久| 麻豆一二三区av精品| 欧美成人一区二区免费高清观看| 无遮挡黄片免费观看| 蜜桃久久精品国产亚洲av| 成人漫画全彩无遮挡| 伦精品一区二区三区| 简卡轻食公司| 国产精品一及| 久久久精品大字幕| 毛片一级片免费看久久久久| 在线观看免费视频日本深夜| 久久久久性生活片| 亚洲av免费在线观看| 亚洲av中文字字幕乱码综合| 男人的好看免费观看在线视频| 99视频精品全部免费 在线| 如何舔出高潮| 男女边吃奶边做爰视频| 内射极品少妇av片p| 狂野欧美白嫩少妇大欣赏| 国产男靠女视频免费网站| 免费人成在线观看视频色|