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

    湍流強(qiáng)化傳熱反饋控制機(jī)理

    2011-07-07 15:03:08潘利生史維秀汪健生
    關(guān)鍵詞:塞爾邊界層流向

    潘利生,史維秀,汪健生

    (天津大學(xué)機(jī)械工程學(xué)院,天津 300072)

    湍流強(qiáng)化傳熱反饋控制機(jī)理

    潘利生,史維秀,汪健生

    (天津大學(xué)機(jī)械工程學(xué)院,天津 300072)

    對(duì)流換熱與湍流結(jié)構(gòu)密切相關(guān),運(yùn)用大渦模擬方法,在湍流流場(chǎng)中進(jìn)行強(qiáng)化傳熱反饋控制,分別采用3種反饋控制速度,研究反饋控制對(duì)傳熱的影響機(jī)理和效果,并建立目標(biāo)函數(shù),研究平均強(qiáng)化傳熱效果.結(jié)果表明:在湍流流場(chǎng)中輸入適當(dāng)反饋控制,強(qiáng)化傳熱效果明顯;當(dāng)α為5.0、10.0和15.0時(shí),努塞爾數(shù)最大增幅分別為8.78%、11.17%和14.08%,平均增幅分別為4.80%、6.88%和9.09%;當(dāng)α為15.0時(shí),y+=5處流場(chǎng)溫度比未輸入反饋控制的流場(chǎng)高10,℃左右.

    反饋控制;強(qiáng)化傳熱;大渦模擬;湍流擬序結(jié)構(gòu)

    反饋控制能夠以很小的能量輸入而達(dá)到對(duì)湍流的完全控制,因此備受關(guān)注.文獻(xiàn)[1-3]研究表明,在近壁面區(qū)域進(jìn)行反饋控制,最多可以使壁面摩擦阻力降低25%.在中等雷諾數(shù)及以下,采用直接數(shù)值模擬和大渦模擬進(jìn)行反饋控制的研究已取得了顯著成果[4].

    Choi等[1]于1994年應(yīng)用一種所謂主動(dòng)取消法(active cancellation),依據(jù)近壁面區(qū)域的法向速度(y+=10),在壁面進(jìn)行吹/吸操作而使得壁面摩擦阻力降低了25%.Lee等[2]于1998年提出了一種新的以壁面變量為參照的亞優(yōu)法則,這些變量為展向壁面剪切應(yīng)力和壁面壓力波動(dòng).同時(shí)Lee等[5]于1997年引入了神經(jīng)網(wǎng)絡(luò)并在直接數(shù)值模擬槽道時(shí)獲得了20%的阻力降低.Endo等[6]于2000年對(duì)槽道進(jìn)行了直接數(shù)值模擬,所采用的控制方法是通過(guò)布置一系列離散的壁面?zhèn)鞲衅骱涂梢宰冃蔚目刂破?,發(fā)現(xiàn)這種真實(shí)的控制單元可以推動(dòng)阻力在實(shí)際中的降低.隨著直接數(shù)值模擬應(yīng)用到高雷諾數(shù)的情況,Iwamoto等[7]于2002年認(rèn)為無(wú)論是主動(dòng)取消法還是亞優(yōu)化控制模式在高雷諾數(shù)的槽道流動(dòng)中所起的作用不大.

    MEMS(micro-electro-mechanical systems)的一個(gè)重要應(yīng)用舞臺(tái)就是流動(dòng)控制,MEMS傳感器和控制器為湍流擬序結(jié)構(gòu)的控制創(chuàng)造了機(jī)會(huì)[8].雖然邊界層附近的湍流結(jié)構(gòu)非常微小,但是隨著MEMS的發(fā)展,構(gòu)建亞毫米級(jí)的流動(dòng)傳感器是完全可能的.Lundell[9]對(duì)傳感器和執(zhí)行器的結(jié)構(gòu)布置進(jìn)行研究,采用交叉結(jié)構(gòu)來(lái)優(yōu)化控制.

    眾多學(xué)者應(yīng)用反饋控制,在湍流減阻方面取得了不少成果,對(duì)線性反饋控制[10]的研究也逐漸增多,盡管湍流強(qiáng)化傳熱的技術(shù)[11]仍在不斷進(jìn)步,但較少關(guān)注反饋控制對(duì)傳熱的影響,也較少將反饋控制應(yīng)用于強(qiáng)化傳熱方面.筆者基于強(qiáng)化傳熱的目的,應(yīng)用大渦模擬方法,研究反饋控制湍流強(qiáng)化傳熱的機(jī)理和效果.

    1 數(shù)值計(jì)算模型及控制方程

    湍流流場(chǎng)為長(zhǎng)方形槽道,見(jiàn)圖1.槽道中流體沿x正方向流動(dòng).控制方程為不可壓縮Navier-Stokes方程和連續(xù)性方程.流向(x向)和展向(z向)為周期性邊界條件,頂面和地面為無(wú)滑移壁面.底面熱邊界條件為定溫.

    流向、展向和法向的計(jì)算區(qū)域分別為2πδ、πδ和2δ.其中δ為槽道的半高度.?dāng)?shù)值模擬計(jì)算在定流量情況下進(jìn)行,雷諾數(shù)為160.本文以后提到的()+是基于壁面剪切速度uτ和運(yùn)動(dòng)黏度v的無(wú)量綱數(shù).

    流向、展向和法向3個(gè)方向的網(wǎng)格數(shù)都為72.法向采用非均勻網(wǎng)格,第1層網(wǎng)格在距離壁面y+=0.46的位置.?dāng)?shù)值計(jì)算時(shí)間步長(zhǎng)選擇為Δt+=0.25.反饋控制輸入前的流場(chǎng)為由LES獲得的充分發(fā)展湍流流場(chǎng).

    圖1 計(jì)算模型示意Fig.1 Schematic diagram of simulation model

    大渦模擬將湍流劃分成大尺度渦和小尺度渦.大尺度渦直接進(jìn)行求解,而小尺度渦采用模型來(lái)求解.在數(shù)值計(jì)算中,采用Deardorff box濾波模型來(lái)劃分大尺度渦和小尺度渦的湍流參數(shù).對(duì)于不可壓縮流體,大渦模擬中的質(zhì)量守恒方程、動(dòng)量守恒方程和能量守恒方程如下.

    質(zhì)量守恒方程

    動(dòng)量守恒方程

    能量守恒方程

    2 計(jì)算結(jié)果與分析

    2.1 湍流擬序結(jié)構(gòu)

    運(yùn)用大渦模擬對(duì)計(jì)算區(qū)域進(jìn)行計(jì)算,獲得充分發(fā)展湍流流場(chǎng).充分發(fā)展的湍流,在近壁面區(qū)域,猝發(fā)過(guò)程伴隨著低速流體的堆積,逐漸形成低速流體向外噴射和高速流體向內(nèi)掃掠,低速流體和高速流體的間歇出現(xiàn)形成了條紋結(jié)構(gòu).圖2為y+=10的法向截面的條紋結(jié)構(gòu).

    圖2 y+=10處的條紋結(jié)構(gòu)Fig.2 Stripe structure at y+=10

    近壁面區(qū)準(zhǔn)流向渦被證明對(duì)湍流速度邊界層和熱邊界層起主要作用,對(duì)外區(qū)流場(chǎng)影響很大.圖3所示為某流向截面的流向渦結(jié)構(gòu),從圖中可以看出,在近壁面區(qū)域存在著大量流向渦,這些渦結(jié)構(gòu)數(shù)量上多于核心區(qū),對(duì)近壁面區(qū)域的阻力和換熱影響很大.為了達(dá)到強(qiáng)化傳熱的目的,必須對(duì)近壁面區(qū)域的流向渦加以研究和控制.

    圖3 某流向截面的流向渦結(jié)構(gòu)Fig.3 Streamwise vorticities at a streamwise cross section

    2.2 反饋控制模式

    Choi等[1]于1994年提出法向速度控制方法(V-control),通過(guò)在y+=10~15處監(jiān)測(cè)流場(chǎng),在底面進(jìn)行吹/吸控制,能夠達(dá)到減弱壁面渦結(jié)構(gòu)的作用.

    經(jīng)過(guò)充分發(fā)展湍流數(shù)值計(jì)算,在y+=10處可以很好地監(jiān)測(cè)到近壁面區(qū)域的第1層準(zhǔn)流向渦結(jié)構(gòu).故在y+=10處設(shè)置監(jiān)測(cè)面,監(jiān)測(cè)此面內(nèi)的法向速度和展向速度,然后將監(jiān)測(cè)結(jié)果返回到控制程序,進(jìn)行判斷,若監(jiān)測(cè)結(jié)果顯示近壁面區(qū)域存在準(zhǔn)流向渦結(jié)構(gòu),則在壁面采取V-control方法進(jìn)行控制操作,見(jiàn)圖4.在V-control中,控制速度表達(dá)式為

    分別取α為5.0、10.0和15.0進(jìn)行數(shù)值計(jì)算.然后對(duì)比不同控制速度下的底面?zhèn)鳠崆闆r.

    圖4 反饋控制示意Fig.4 Sketch map of feedback control method

    2.3 強(qiáng)化傳熱機(jī)理

    在湍流流場(chǎng)中,近壁面區(qū)域基本由相對(duì)強(qiáng)的準(zhǔn)流向渦所占據(jù),這些渦與近壁面區(qū)域流體的噴出和掃掠過(guò)程緊密相關(guān),并且是近壁區(qū)雷諾應(yīng)力產(chǎn)生的根源,緩沖層的流向渦在低速條紋形成中起明顯的作用.準(zhǔn)流向渦常常以非對(duì)稱反轉(zhuǎn)渦對(duì)出現(xiàn),有時(shí)也呈單個(gè)出現(xiàn).

    在近壁面區(qū)域,由于準(zhǔn)流向渦的作用,不斷地有流體從壁面區(qū)域向外噴射,有流體從外區(qū)下掃.流體向外噴射減小近壁面區(qū)域的速度梯度和溫度梯度,流體下掃有利于增大近壁面區(qū)域速度梯度和溫度梯度,減薄速度邊界層和溫度邊界層,有利于傳熱的加強(qiáng).故在近壁面區(qū)域設(shè)置監(jiān)測(cè)面監(jiān)測(cè)流場(chǎng),當(dāng)監(jiān)測(cè)到的信息顯示近壁面區(qū)域存在準(zhǔn)流向渦并將形成流體向外噴射時(shí),即時(shí)采取控制,抑制流體向上噴出,盡量形成下掃,減薄速度邊界層和溫度邊界層,以達(dá)到強(qiáng)化傳熱的目的.

    2.4 反饋控制對(duì)傳熱的效果分析圖5為同一時(shí)刻,y+=10處溫度場(chǎng)云圖.從圖中可以看出,在反饋控制輸入?yún)^(qū)域之前,溫度場(chǎng)基本相同,在反饋控制區(qū)域輸入控制的流場(chǎng),換熱得到加強(qiáng),溫度場(chǎng)溫度顯著提高,流場(chǎng)最高溫度較未輸入反饋控制的流場(chǎng)高10,℃左右.

    圖5 y+=10處的溫度分布云圖Fig.5 Contour of temperature at y+=10

    圖6 為壁面努塞爾數(shù)隨無(wú)量綱時(shí)間的變化情況.從圖6中可以看出,通過(guò)在充分發(fā)展湍流近壁面區(qū)域輸入反饋控制,壁面努塞爾數(shù)顯著提高:

    在t+=0時(shí)輸入反饋控制后,努塞爾數(shù)突然增大,當(dāng)α=5.0時(shí),增幅可達(dá)25.87%,當(dāng)α=10.0時(shí),增幅可達(dá)32.27%,當(dāng)α=15.0時(shí),增幅可達(dá)38.98%.在近壁面流場(chǎng)中輸入法向反饋控制,即使是很小的能量輸入,也能大幅度減薄溫度邊界層,從而在輸入反饋控制之初,強(qiáng)化傳熱效果異常明顯.

    t+增加到100以后,強(qiáng)化傳熱效果趨于平緩,努塞爾數(shù)比隨時(shí)間呈現(xiàn)波浪型變化,當(dāng)α=5.0時(shí),最大增幅可達(dá)8.78%,當(dāng)α=10.0時(shí),最大增幅可達(dá)11.17%,當(dāng)α=15.0時(shí),最大增幅可達(dá)14.08%.

    圖6 壁面努塞爾數(shù)變化Fig.6 Variation of Nusselt number ratio at bottom

    3種輸入反饋控制的情況下,努塞爾數(shù)隨時(shí)間的變化趨勢(shì)相似.由于在近壁面輸入反饋控制的能量非常小,對(duì)流場(chǎng)的整體流動(dòng)狀況影響很小,只是減薄了近壁面區(qū)域邊界層厚度,強(qiáng)化了傳熱,故3種反饋控制下,努塞爾數(shù)的變化趨勢(shì)相同.

    湍流流場(chǎng)和傳熱效果隨時(shí)間不斷變化,經(jīng)過(guò)反饋控制,強(qiáng)化傳熱的效果隨時(shí)間不盡相同.故最大強(qiáng)化傳熱數(shù)值并不能完全反映綜合強(qiáng)化傳熱效果,有必要研究在計(jì)算時(shí)間范圍內(nèi),強(qiáng)化傳熱的平均值.為研究當(dāng)反饋控制作用達(dá)到平穩(wěn)后,反饋控制對(duì)強(qiáng)化傳熱的平均效果,建立目標(biāo)函數(shù)為

    式中:h為對(duì)流換熱系數(shù);λ為流體導(dǎo)熱系數(shù);1t+為初始計(jì)算時(shí)間;2t+為終止計(jì)算時(shí)間;Nu為輸入反饋控制的努塞爾數(shù);0Nu為未輸入反饋控制的努塞爾數(shù).

    鑒于在反饋控制之初,反饋控制效果不穩(wěn)定,努塞爾數(shù)變化劇烈,為研究反饋控制的平均傳熱效果,將t+=100之前視為反饋控制非平穩(wěn)區(qū),重點(diǎn)研究t+在100~500范圍內(nèi)Nu/Nu0的平均值.應(yīng)用式(6)對(duì)3種反饋控制作用下數(shù)據(jù)進(jìn)行處理.計(jì)算結(jié)果表明,t+在100~500范圍內(nèi),當(dāng)α=5.0時(shí),努塞爾數(shù)平均增幅為4.80%,當(dāng)α=10.0時(shí),努塞爾數(shù)平均增幅為6.88%,當(dāng)α=15.0時(shí),努塞爾數(shù)平均增幅為9.09%.

    表1給出了3種反饋控制情況下,努塞爾數(shù)的最大增幅和平均增幅.

    表1 不同反饋控制速度時(shí)的強(qiáng)化傳熱效果對(duì)比Tab.1 Comparison of heat transfer enhancement effects with different control velocities

    3 結(jié) 論

    (1)在近壁面區(qū)域輸入適當(dāng)反饋控制,可以顯著增強(qiáng)換熱.

    (2)當(dāng)α=5.0時(shí),最大增幅可達(dá)8.78%;當(dāng)α= 10.0時(shí),最大增幅可達(dá)11.17%;當(dāng)α=15.0時(shí),最大增幅可達(dá)14.08%.

    (3)引入目標(biāo)函數(shù),研究平均強(qiáng)化傳熱效果,t+在100~500范圍內(nèi),當(dāng)α=5.0時(shí),努塞爾數(shù)平均增幅為4.80%;當(dāng)α=10.0時(shí),努塞爾數(shù)平均增幅為6.88%;當(dāng)α=15.0時(shí),努塞爾數(shù)平均增幅為9.09%.

    [1] Choi H,Moin P,Kim J. Active turbulence control for drag reduction in wall bounded flows[J]. Fluid Mech,1994,262:75-110.

    [2] Lee C,Kim J,Choi H. Suboptimal control of turbulent channel flow for drag reduction[J]. Fluid Mech,1998,358:245-258.

    [3] Min T,Kim J. Effects of hydrophobic surface on skinfriction drag[J]. Physics of Fluids,2004,16(7):55-58.

    [4] Kim J. Control of turbulent boundary layers[J]. Physics of Fluids,2003,15(5):1093-1105.

    [5] Lee C,Kim J,Babcock R. Application of neural networks to turbulence control for drag reduction[J]. Physics of Fluids,1997,9(6):1740-1747.

    [6] Endo t,Kasagi N,Suzuki Y. Feedback control of wall turbulence with wall deformation[J]. Heat Fluid Flow,2000,21(5):568-575.

    [7] Iwamoto K,Suzuki Y,Kasagi N. Reynolds number

    effect on wall turbulence:Toward effective feedback

    control[J]. Heat Fluid Flow,2002,23(5):678-689.[8] Lofdahl L˙˙,Gad-el-Hak M. MEMS applications in turbulence and flow control[J]. Progress in Aerospace Sciences,1999,35(2):101-203.

    [9] Lundell F. Reactive control of transition induced by freestream turbulence:An experimental demonstration[J]. Fluid Mech,2007,585:41-71.

    [10] Monokrousos A,Brandt L,Schlatter P,et al. DNS and LES of estimation and control of transition in boundary layers subject[J]. Heat and Fluid Flow,2008,29(3):841-855.

    [11] Yang Yue-Tzu,Hwang Ming-Lu. Numerical simulation of turbulent fluid flow and heat transfer characteristics in heat exchangers fitted with porous media[J]. Heat and Mass Transfer,2009,52(13/14):2956-2965.

    Feedback Control Mechanism for Turbulent Heat Transfer Enhancement

    PAN Li-sheng,SHI Wei-xiu,WANG Jian-sheng
    (School of Mechanical Engineering,Tianjin University,Tianjin 300072,China)

    Convective heat transfer is closely related to turbulence structure. Using large eddy simulation(LES),feedback control was executed in turbulent flow field with three kinds of control velocities,mechanism and effect of feedback control upon heat transfer were studied,and an objective function was established to study the average effect on heat transfer enhancement. The results show that,inputting feedback control in turbulent flow field,the heat transfer is enhanced remarkably. When α is 5.0,10.0 and 15.0,the largest enhancements of Nusselt number are 8.78%,11.17% and 14.08%,respectively,and the averages of which are 4.80%,6.88% and 9.09% respectively. When α is 15.0,the highest temperature is about 10 ℃ higher than that without control at y+=5.

    feedback-control;heat transfer enhancement;large eddy simulation;turbulent coherent structure

    TK124

    A

    0493-2137(2011)07-0655-04

    2010-04-02;

    2010-06-01.

    國(guó)家自然科學(xué)基金資助項(xiàng)目(50976079);國(guó)家高技術(shù)研究發(fā)展計(jì)劃(863計(jì)劃)資助項(xiàng)目(2009AA05Z431).

    潘利生(1982— ),男,博士研究生.

    潘利生,Lisheng_P@163.com.

    猜你喜歡
    塞爾邊界層流向
    如果地球被我們吃掉了
    小溪??!流向遠(yuǎn)方
    井岡教育(2020年6期)2020-12-14 03:04:42
    基于HIFiRE-2超燃發(fā)動(dòng)機(jī)內(nèi)流道的激波邊界層干擾分析
    十大漲幅、換手、振副、資金流向
    淺談點(diǎn)-連式ATP系統(tǒng)在埃塞爾比亞輕軌中的應(yīng)用
    流向逆轉(zhuǎn)的啟示
    一類具有邊界層性質(zhì)的二次奇攝動(dòng)邊值問(wèn)題
    非特征邊界的MHD方程的邊界層
    鄭州市春季邊界層風(fēng)氣候變化研究
    河南科技(2014年23期)2014-02-27 14:19:08
    秋天的流向(組詩(shī))
    鴨綠江(2013年11期)2013-03-11 19:42:04
    最新中文字幕久久久久| 欧美成人午夜免费资源| 国产精品人妻久久久久久| 人体艺术视频欧美日本| 成人国产麻豆网| www.熟女人妻精品国产 | 边亲边吃奶的免费视频| www.色视频.com| 成人18禁高潮啪啪吃奶动态图| 丰满少妇做爰视频| 如何舔出高潮| 久久鲁丝午夜福利片| 日韩一本色道免费dvd| 五月伊人婷婷丁香| 观看美女的网站| a级毛片黄视频| 高清在线视频一区二区三区| 伦理电影大哥的女人| 一区二区三区精品91| 国产国拍精品亚洲av在线观看| 国产成人a∨麻豆精品| 久久人人爽人人片av| 国产日韩欧美在线精品| 免费高清在线观看日韩| www.色视频.com| 亚洲,一卡二卡三卡| 人妻系列 视频| 国产色爽女视频免费观看| www.av在线官网国产| 午夜日本视频在线| 搡老乐熟女国产| 国产老妇伦熟女老妇高清| 97在线人人人人妻| 成人亚洲欧美一区二区av| 国产精品国产三级专区第一集| 在线观看三级黄色| 亚洲精品日本国产第一区| 精品人妻一区二区三区麻豆| 一级毛片黄色毛片免费观看视频| 成人国产av品久久久| 久久精品国产自在天天线| 青春草亚洲视频在线观看| 中文乱码字字幕精品一区二区三区| 久久精品国产a三级三级三级| 在线天堂中文资源库| freevideosex欧美| 制服人妻中文乱码| 国产乱来视频区| 日韩一区二区视频免费看| 成年美女黄网站色视频大全免费| 亚洲综合色网址| videossex国产| 天天躁夜夜躁狠狠躁躁| 国产男女超爽视频在线观看| 一级毛片 在线播放| 久久国产亚洲av麻豆专区| 性色av一级| 少妇被粗大的猛进出69影院 | 精品少妇内射三级| 久久久久久久国产电影| 人人妻人人澡人人爽人人夜夜| 日本色播在线视频| 91精品伊人久久大香线蕉| 亚洲成人一二三区av| 少妇的逼水好多| av有码第一页| 视频在线观看一区二区三区| 综合色丁香网| 纯流量卡能插随身wifi吗| 亚洲国产精品999| 亚洲精品久久久久久婷婷小说| 国产精品三级大全| 黄色视频在线播放观看不卡| 国产精品不卡视频一区二区| 亚洲熟女精品中文字幕| 久久韩国三级中文字幕| 男女午夜视频在线观看 | 亚洲欧美一区二区三区国产| 久久久精品免费免费高清| 国产成人精品久久久久久| 成人国语在线视频| 大香蕉97超碰在线| 久久久欧美国产精品| 久久久久久久久久久久大奶| 狂野欧美激情性xxxx在线观看| a级毛片黄视频| 亚洲精品久久成人aⅴ小说| 亚洲精品aⅴ在线观看| 久久人人爽人人爽人人片va| 97精品久久久久久久久久精品| 午夜视频国产福利| 成人手机av| 两个人看的免费小视频| 中文字幕亚洲精品专区| 亚洲成人av在线免费| 日本wwww免费看| av国产久精品久网站免费入址| 伊人久久国产一区二区| 午夜福利网站1000一区二区三区| 精品第一国产精品| 亚洲国产成人一精品久久久| 国产精品.久久久| 亚洲欧美成人综合另类久久久| 国产精品久久久久久av不卡| 亚洲伊人色综图| 成人亚洲欧美一区二区av| 我的女老师完整版在线观看| 日本欧美国产在线视频| 精品人妻在线不人妻| 午夜日本视频在线| 国产综合精华液| 国产日韩欧美亚洲二区| 日本免费在线观看一区| 国产成人精品婷婷| 少妇熟女欧美另类| 亚洲,欧美,日韩| 国产熟女午夜一区二区三区| 久久精品国产综合久久久 | 欧美另类一区| 亚洲国产成人一精品久久久| 波野结衣二区三区在线| 亚洲精品乱久久久久久| 成年人午夜在线观看视频| 亚洲国产看品久久| 街头女战士在线观看网站| 国产成人av激情在线播放| 免费看av在线观看网站| 欧美日韩综合久久久久久| 亚洲av.av天堂| 免费日韩欧美在线观看| 国产色爽女视频免费观看| 成人毛片60女人毛片免费| 中文精品一卡2卡3卡4更新| 色婷婷久久久亚洲欧美| 考比视频在线观看| av在线观看视频网站免费| 黄色怎么调成土黄色| 亚洲精品国产av蜜桃| 亚洲成色77777| 亚洲一级一片aⅴ在线观看| 色婷婷久久久亚洲欧美| 热99久久久久精品小说推荐| 亚洲欧洲国产日韩| 国产片特级美女逼逼视频| 午夜免费观看性视频| 国产黄色免费在线视频| 99九九在线精品视频| 又黄又爽又刺激的免费视频.| 欧美日韩综合久久久久久| 国产男人的电影天堂91| 国产69精品久久久久777片| 久久鲁丝午夜福利片| 少妇人妻久久综合中文| 国产 精品1| a级片在线免费高清观看视频| 制服诱惑二区| 99久国产av精品国产电影| 成人黄色视频免费在线看| 各种免费的搞黄视频| 男人爽女人下面视频在线观看| 考比视频在线观看| 免费在线观看黄色视频的| 精品久久久精品久久久| 午夜福利,免费看| 免费在线观看黄色视频的| 一级毛片我不卡| 久久狼人影院| 看非洲黑人一级黄片| 日韩欧美精品免费久久| a级片在线免费高清观看视频| 观看av在线不卡| 9191精品国产免费久久| 国产精品久久久av美女十八| 街头女战士在线观看网站| 国产免费视频播放在线视频| 久久精品aⅴ一区二区三区四区 | 日韩av免费高清视频| 久久久久国产精品人妻一区二区| 国产白丝娇喘喷水9色精品| 一本色道久久久久久精品综合| 最近的中文字幕免费完整| 久久久久精品人妻al黑| 国产男女内射视频| 午夜精品国产一区二区电影| 欧美精品一区二区大全| 丝袜美足系列| 亚洲国产毛片av蜜桃av| 亚洲av中文av极速乱| 国产免费视频播放在线视频| 国产精品不卡视频一区二区| 最近最新中文字幕免费大全7| 亚洲欧美日韩另类电影网站| 亚洲第一区二区三区不卡| 两个人看的免费小视频| 亚洲av在线观看美女高潮| 国产永久视频网站| 黄色毛片三级朝国网站| 国产精品久久久av美女十八| freevideosex欧美| 高清欧美精品videossex| 久久久欧美国产精品| 99热国产这里只有精品6| 国产亚洲精品久久久com| 久久亚洲国产成人精品v| 日韩不卡一区二区三区视频在线| 欧美人与善性xxx| 亚洲精品久久午夜乱码| 国产色爽女视频免费观看| 日产精品乱码卡一卡2卡三| 日日爽夜夜爽网站| 精品亚洲乱码少妇综合久久| 国产黄色免费在线视频| 久久韩国三级中文字幕| 亚洲,欧美精品.| 少妇人妻久久综合中文| 另类精品久久| 久久久久久久久久久久大奶| 一区二区三区乱码不卡18| kizo精华| 爱豆传媒免费全集在线观看| 久热久热在线精品观看| 寂寞人妻少妇视频99o| 免费av中文字幕在线| 久久久久久久大尺度免费视频| 国产精品国产av在线观看| 亚洲av成人精品一二三区| 久久久久精品性色| 国产av精品麻豆| 久久ye,这里只有精品| 成年av动漫网址| 国产成人精品无人区| 久久97久久精品| 欧美97在线视频| 欧美人与性动交α欧美精品济南到 | 免费播放大片免费观看视频在线观看| 中文字幕av电影在线播放| 伊人久久国产一区二区| 午夜av观看不卡| 久久午夜福利片| 两个人免费观看高清视频| 在线观看免费日韩欧美大片| 精品久久久精品久久久| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 国产成人91sexporn| 三级国产精品片| 国产毛片在线视频| 国产成人精品婷婷| 韩国精品一区二区三区 | 高清毛片免费看| 这个男人来自地球电影免费观看 | 中国三级夫妇交换| 自线自在国产av| 国产黄色免费在线视频| 人妻少妇偷人精品九色| 午夜福利在线观看免费完整高清在| 久久精品久久精品一区二区三区| 宅男免费午夜| 欧美精品国产亚洲| 亚洲少妇的诱惑av| 最新中文字幕久久久久| 黄色 视频免费看| 母亲3免费完整高清在线观看 | 在线免费观看不下载黄p国产| 中文字幕另类日韩欧美亚洲嫩草| 国产成人av激情在线播放| 国产淫语在线视频| 黑人高潮一二区| 三上悠亚av全集在线观看| 丝袜人妻中文字幕| 免费日韩欧美在线观看| 久久精品人人爽人人爽视色| 久久热在线av| 国产精品不卡视频一区二区| 丰满乱子伦码专区| 久久这里有精品视频免费| 我的女老师完整版在线观看| 国产在线一区二区三区精| 蜜臀久久99精品久久宅男| 亚洲四区av| 色94色欧美一区二区| 亚洲情色 制服丝袜| 国产免费又黄又爽又色| 18禁裸乳无遮挡动漫免费视频| 欧美日韩亚洲高清精品| 免费观看a级毛片全部| 亚洲精品第二区| 免费观看在线日韩| 亚洲欧美清纯卡通| 国产成人精品在线电影| 亚洲 欧美一区二区三区| 国产精品麻豆人妻色哟哟久久| 永久免费av网站大全| 亚洲天堂av无毛| 美女国产高潮福利片在线看| 亚洲,欧美,日韩| 最新中文字幕久久久久| 高清不卡的av网站| 亚洲国产欧美日韩在线播放| 久久久久久久久久成人| 成人国产麻豆网| 亚洲国产精品成人久久小说| 色94色欧美一区二区| 高清av免费在线| 亚洲精品久久午夜乱码| 一区二区日韩欧美中文字幕 | 国产日韩一区二区三区精品不卡| 久久久久久久亚洲中文字幕| 亚洲熟女精品中文字幕| 久久精品国产综合久久久 | 国产成人精品婷婷| 亚洲欧美色中文字幕在线| 亚洲伊人久久精品综合| 高清毛片免费看| 美女主播在线视频| 欧美 日韩 精品 国产| 色视频在线一区二区三区| 久久狼人影院| 国产精品三级大全| 人妻 亚洲 视频| 久久人人爽人人爽人人片va| 久久久a久久爽久久v久久| 国产熟女午夜一区二区三区| 一级a做视频免费观看| 伊人久久国产一区二区| 最新中文字幕久久久久| 成年av动漫网址| av在线观看视频网站免费| 久久久国产一区二区| 少妇 在线观看| 成年人午夜在线观看视频| 热99久久久久精品小说推荐| 久久精品久久久久久噜噜老黄| 国产熟女午夜一区二区三区| 亚洲综合精品二区| 两个人免费观看高清视频| 波多野结衣一区麻豆| 日韩免费高清中文字幕av| 两个人免费观看高清视频| 五月伊人婷婷丁香| 99久久精品国产国产毛片| 久久久欧美国产精品| 免费播放大片免费观看视频在线观看| 九色亚洲精品在线播放| 免费av不卡在线播放| 啦啦啦在线观看免费高清www| 久久久久国产精品人妻一区二区| 国产男女超爽视频在线观看| 制服人妻中文乱码| 黑人巨大精品欧美一区二区蜜桃 | 黄色毛片三级朝国网站| 久热这里只有精品99| 久久99精品国语久久久| 人妻人人澡人人爽人人| 色5月婷婷丁香| 精品一区二区三区四区五区乱码 | 国产精品嫩草影院av在线观看| 伦精品一区二区三区| 成人漫画全彩无遮挡| 男人添女人高潮全过程视频| 久久久久视频综合| 久久精品久久久久久久性| 18在线观看网站| 在线免费观看不下载黄p国产| 美女脱内裤让男人舔精品视频| 精品一区在线观看国产| 久久精品国产亚洲av天美| 国产极品粉嫩免费观看在线| 国产乱来视频区| 亚洲精品久久久久久婷婷小说| 多毛熟女@视频| 免费少妇av软件| 欧美精品一区二区大全| 久久久久精品性色| 亚洲成色77777| 99久久综合免费| 久久国产精品男人的天堂亚洲 | 国产亚洲av片在线观看秒播厂| 日本vs欧美在线观看视频| 免费黄色在线免费观看| 97超碰精品成人国产| av一本久久久久| 一个人免费看片子| 香蕉丝袜av| 午夜91福利影院| 免费在线观看完整版高清| 十分钟在线观看高清视频www| 丁香六月天网| 国产国语露脸激情在线看| 九九爱精品视频在线观看| 亚洲av国产av综合av卡| 中文乱码字字幕精品一区二区三区| 在线看a的网站| 少妇猛男粗大的猛烈进出视频| 日韩视频在线欧美| 2018国产大陆天天弄谢| 精品人妻一区二区三区麻豆| 在线观看三级黄色| 在线观看人妻少妇| 国产国语露脸激情在线看| 成年av动漫网址| 91国产中文字幕| 国产精品久久久久久精品古装| 人成视频在线观看免费观看| 一边摸一边做爽爽视频免费| 超色免费av| 国产乱来视频区| 22中文网久久字幕| 激情视频va一区二区三区| 国产女主播在线喷水免费视频网站| 国产 精品1| av.在线天堂| 久久久久网色| 亚洲成人av在线免费| 夜夜骑夜夜射夜夜干| 99国产精品免费福利视频| 97在线视频观看| 日韩欧美一区视频在线观看| 亚洲国产欧美在线一区| 在线观看免费视频网站a站| 国产精品国产三级国产av玫瑰| 天天影视国产精品| 我的女老师完整版在线观看| 国产色婷婷99| 日日撸夜夜添| 久久狼人影院| 自拍欧美九色日韩亚洲蝌蚪91| 秋霞伦理黄片| 2018国产大陆天天弄谢| 尾随美女入室| 国产成人91sexporn| 国产一区有黄有色的免费视频| 国产熟女午夜一区二区三区| 97在线人人人人妻| 视频区图区小说| av在线app专区| 免费播放大片免费观看视频在线观看| 丝袜喷水一区| 国产成人91sexporn| a级毛色黄片| 免费观看a级毛片全部| 99热这里只有是精品在线观看| 欧美日韩视频精品一区| 妹子高潮喷水视频| 国产一区二区在线观看av| 女的被弄到高潮叫床怎么办| 91aial.com中文字幕在线观看| 亚洲精品国产av蜜桃| 久久99蜜桃精品久久| 亚洲精品国产av蜜桃| 国产精品无大码| 少妇人妻精品综合一区二区| 最近中文字幕高清免费大全6| 赤兔流量卡办理| 国产成人av激情在线播放| 日韩在线高清观看一区二区三区| 久久这里只有精品19| 超碰97精品在线观看| 国产极品天堂在线| 精品卡一卡二卡四卡免费| kizo精华| 极品人妻少妇av视频| 日本欧美国产在线视频| 久久久久久久久久人人人人人人| 欧美国产精品一级二级三级| 亚洲国产看品久久| freevideosex欧美| av在线播放精品| 日韩精品免费视频一区二区三区 | 国产国语露脸激情在线看| 中文乱码字字幕精品一区二区三区| 日韩欧美一区视频在线观看| 国产一区二区激情短视频 | 有码 亚洲区| 亚洲欧美一区二区三区国产| 婷婷色综合大香蕉| 亚洲精品第二区| 欧美日韩综合久久久久久| 婷婷色综合www| 国产成人精品在线电影| 亚洲国产色片| 中文字幕人妻丝袜制服| a级毛色黄片| 免费高清在线观看视频在线观看| 久久人人97超碰香蕉20202| 婷婷色综合大香蕉| 18禁动态无遮挡网站| 一本—道久久a久久精品蜜桃钙片| 欧美成人精品欧美一级黄| 亚洲一级一片aⅴ在线观看| 高清欧美精品videossex| 女人久久www免费人成看片| 高清欧美精品videossex| 免费观看av网站的网址| 青春草国产在线视频| 日韩一本色道免费dvd| 18禁动态无遮挡网站| 国产一区二区在线观看日韩| 在线观看一区二区三区激情| 男女国产视频网站| 欧美激情 高清一区二区三区| 美女内射精品一级片tv| 成年人免费黄色播放视频| 国产亚洲精品久久久com| 久久99精品国语久久久| 亚洲国产精品国产精品| 五月玫瑰六月丁香| 亚洲成av片中文字幕在线观看 | www.色视频.com| 国产老妇伦熟女老妇高清| 国产淫语在线视频| 日韩成人av中文字幕在线观看| 一级毛片黄色毛片免费观看视频| 欧美 亚洲 国产 日韩一| 大香蕉久久成人网| 另类精品久久| 久久久久视频综合| 亚洲欧美日韩卡通动漫| 熟女电影av网| 伦理电影免费视频| 久久女婷五月综合色啪小说| 午夜激情av网站| 午夜福利网站1000一区二区三区| 飞空精品影院首页| 中文字幕另类日韩欧美亚洲嫩草| 国产片特级美女逼逼视频| 亚洲色图 男人天堂 中文字幕 | 日韩一本色道免费dvd| 国产爽快片一区二区三区| 亚洲内射少妇av| 久久毛片免费看一区二区三区| 精品少妇内射三级| 久久av网站| 久久97久久精品| 波野结衣二区三区在线| 少妇熟女欧美另类| 1024视频免费在线观看| 亚洲伊人色综图| 一级毛片 在线播放| 中国美白少妇内射xxxbb| 看非洲黑人一级黄片| 男人爽女人下面视频在线观看| 香蕉丝袜av| 亚洲成人一二三区av| 中文字幕人妻熟女乱码| 男女无遮挡免费网站观看| 亚洲av.av天堂| 色5月婷婷丁香| 欧美 亚洲 国产 日韩一| 亚洲成色77777| 亚洲精品第二区| 国产成人免费观看mmmm| 亚洲三级黄色毛片| 卡戴珊不雅视频在线播放| 91午夜精品亚洲一区二区三区| 麻豆精品久久久久久蜜桃| 18禁裸乳无遮挡动漫免费视频| 亚洲国产精品999| 国产视频首页在线观看| 香蕉丝袜av| 亚洲精品日本国产第一区| 日韩伦理黄色片| 80岁老熟妇乱子伦牲交| 青青草视频在线视频观看| 男女免费视频国产| 亚洲四区av| 看非洲黑人一级黄片| 国产av一区二区精品久久| 免费大片18禁| 国产女主播在线喷水免费视频网站| av不卡在线播放| 日韩一本色道免费dvd| 97在线人人人人妻| 一级毛片 在线播放| 亚洲在久久综合| 亚洲精品久久午夜乱码| 欧美精品人与动牲交sv欧美| 欧美成人午夜免费资源| 成人国产av品久久久| 国产欧美另类精品又又久久亚洲欧美| 美女视频免费永久观看网站| 国产一区二区在线观看av| 尾随美女入室| 日日爽夜夜爽网站| 亚洲精品乱码久久久久久按摩| 丝袜在线中文字幕| 97在线视频观看| 国产av码专区亚洲av| 九色亚洲精品在线播放| 女人精品久久久久毛片| 亚洲成av片中文字幕在线观看 | 日韩三级伦理在线观看| 亚洲成av片中文字幕在线观看 | 在线观看美女被高潮喷水网站| 中文字幕另类日韩欧美亚洲嫩草| 国产精品蜜桃在线观看| 精品一区二区免费观看| 免费在线观看完整版高清| 日本欧美国产在线视频| 久久99蜜桃精品久久| 精品视频人人做人人爽| 亚洲欧洲日产国产| 全区人妻精品视频| 中文字幕亚洲精品专区| 爱豆传媒免费全集在线观看| 伦理电影大哥的女人| a级毛色黄片| 日本色播在线视频| 考比视频在线观看| 色视频在线一区二区三区| 精品人妻偷拍中文字幕| 亚洲av在线观看美女高潮| 国产成人免费观看mmmm| 免费大片黄手机在线观看| 成人漫画全彩无遮挡| 免费高清在线观看日韩|