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

    溫度對流場中液滴運動軌跡偏轉(zhuǎn)作用的研究

    2014-02-10 18:14:12唐洪濤耿勝民孟輝波
    天津科技大學(xué)學(xué)報 2014年6期
    關(guān)鍵詞:平均溫度液滴剪力

    唐洪濤,耿勝民,孟輝波

    (1. 天津科技大學(xué)機械工程學(xué)院,天津 300222;2. 沈陽化工大學(xué)機械工程學(xué)院,沈陽 110142)

    溫度對流場中液滴運動軌跡偏轉(zhuǎn)作用的研究

    唐洪濤1,耿勝民1,孟輝波2

    (1. 天津科技大學(xué)機械工程學(xué)院,天津 300222;2. 沈陽化工大學(xué)機械工程學(xué)院,沈陽 110142)

    推導(dǎo)了擴散相液滴運動軌跡模型和溫度場中平均流修正模型,并優(yōu)化了計算方法.直接數(shù)值模擬結(jié)果表明:平均溫度和溫度分布方式不同都會引起溫度剖面和平均流剖面發(fā)生相應(yīng)的變化,但它們的表現(xiàn)形式各不相同.因為溫度場的變化會改變液滴所受到的不平衡剪應(yīng)力,所以溫度不同液滴的運動軌跡也不同,這種變化的動力來源于溫度剖面和平均流剖面發(fā)揮的作用,它顯著地影響液滴的聚結(jié)效率,因此,適當(dāng)?shù)卣{(diào)節(jié)流場溫度可以改善流場的聚結(jié)條件.在層流流場中液滴的運動方式受溫度分布的影響具有一定的隨機性,但隨機變化因素中包含確定性因素.液滴的聚結(jié)效率在一定程度上依賴于流場的平均溫度,但溫度分布方式也是影響聚結(jié)效率的不可忽視的因素,溫度是理論反演計算和簡化工程計算的必要條件.

    溫度;聚結(jié);液滴;平均流剖面;熱導(dǎo)率

    是隨機過程而是在某種范圍內(nèi)可以確定地被描述.他們提出了“雙膜模型”,該模型考慮到了層膜及其受到曳力的液體分別為相互獨立的兩相,二者在兩相的界面處發(fā)生動量傳遞,而且相界面具有流動的特性.計算結(jié)果表明,液滴在傾斜斜板間的聚結(jié)過程在某些操作范圍內(nèi)主要受液力影響.Goldman等[6]的研究結(jié)果表明:球形液滴與傾斜板的接觸不是普通的接觸,而是沿著傾斜壁向下滾動過程中存在“滑動”.Rowley等[7]從理論方面研究了斜板分層器中液滴的運動方式.然而,上述研究建立的理論模型計算結(jié)果與實驗結(jié)果偏差較大[4–5].

    鑒于此種情況,筆者[8–11]在計算分析過程中考慮了斜板流場中的基本流,使得計算結(jié)果比Meon等[4]的計算結(jié)果更加接近實驗結(jié)果,并發(fā)現(xiàn)除了液滴的平動以外由瞬時流場引起的液滴的旋轉(zhuǎn)運動對液滴的聚結(jié)也發(fā)揮重要作用.進一步的研究表明:僅考慮液滴在層膜上的運動方式來判斷聚結(jié)效果是不完全的,液滴在界面波誘發(fā)下產(chǎn)生的運動方式對聚結(jié)的影響也是不可忽視的因素[12].事實上,分層二相流流體流動是以界面波的形式進行的,由于界面波的存在使得流場內(nèi)液滴聚結(jié)更加復(fù)雜多變,它也是導(dǎo)致流場內(nèi)液滴聚結(jié)實驗可重復(fù)性差的原因之一;基本流場物性的變化(密度、黏度以及界面張力等)對分散相液滴在分層界面上的受力特征、運動方式及液滴結(jié)構(gòu)形狀變化具有顯著的影響.該研究分析了改善或抑制液滴的聚結(jié)因素的存在性,揭示了分層二相流流場的復(fù)雜性和多變性[13].

    過去的研究中一直忽視了溫度場對液滴聚結(jié)的影響,其原因在于:溫度與流場中的各種參數(shù)是耦合的,這給數(shù)值模擬帶來極大的困難;溫度引起的復(fù)雜性和隨機性導(dǎo)致數(shù)值模擬困難和實驗可重復(fù)性差.因此,本文統(tǒng)計分析了流場溫度變化和溫度分布變化對聚結(jié)效率的影響,研究了引起液滴運動軌跡的改變和聚結(jié)效率顯著變化的原因,在此基礎(chǔ)上,進一步分析了平均溫度和溫度分布方式對聚結(jié)效率的影響方式.

    1 液滴軌跡理論模型

    1.1 溫度場中平均流修正模型

    本文的計算始于層流模型,計算域如圖1所示,坐標(biāo)x、y和z分別對應(yīng)流向、壁法向和展向.討論的問題是稀溶液的定常流動.Navier-Stokes方程的連續(xù)相的連續(xù)方程、動量方程及能量方程分別表達為

    式中:u、v分別為對應(yīng)流向、壁法向的基本流速度;T為流場溫度;dρ、dμ、dc、dk和cρ、cμ、cc、ck分別為分散相和連續(xù)相的密度、動力黏性系數(shù)、比熱容以及熱傳導(dǎo)系數(shù).邊界條件為

    式中:1T表示上板溫度;0T為下板溫度;h表示上下板間距的一半.

    實際上,由于考慮到溫度對流場的影響,上述方程的各個參數(shù)以及速度應(yīng)該均是溫度的函數(shù),求解比較困難.然而,值得注意的是:在指定的溫度范圍內(nèi),流場受溫度影響的程度主要取決于溫度差cT,如果溫度差變化不大,則密度、黏度、傳熱系數(shù)以及表面張力等參數(shù)變化不大,因此,可以認為密度、黏度等物理化學(xué)參數(shù)基本不變,即可以取流場的平均溫度確定物理化學(xué)參數(shù)值,這樣就可以通過式(1)—式(4)確定定常流場中各物理量的數(shù)值,即流場的初步解.為了研究溫度差cT對流場的影響,在后續(xù)的計算中不斷地對初步解進行修正,使流場的瞬時解真正反映流場的真實解,該方法原則上不失一般性.因此,聯(lián)立

    上述方程并結(jié)合邊界條件可得到流場初步解的表達式:

    基本流方程

    溫度剖面方程

    中間變量表達式

    本文研究的是正戊醇和水的稀混合物系,溫度影響下的各參數(shù)的計算方法如下:正戊醇的密度[14]

    1.2 分散相液滴在連續(xù)相中的運動

    作用在液滴上的力(圖2所示)的平衡表達式為

    選用的流體流動區(qū)域為Stokes滯流區(qū),即10?4<Re< 1,則曳力系數(shù)為

    式中:ξ是阻力關(guān)聯(lián)系數(shù),可根據(jù)文獻[4]的計算方法得到.

    液滴終端速度沉降的平均停留時間為

    通過Stokes方程求解沉降速度

    進料的沉降速度為

    進入通道之中的料液的平均速度為

    其中,進入通道內(nèi)的指定流場體積流量為

    流場體積流量Q=2.614×10-5,m3/h.將式(8)方程代入式(27)中,并令Δp為常數(shù),就可以求出Δp.

    1.3 數(shù)值模擬的方法

    (1)計算域取值為:斜板傾角α=45°,板長L=0.96,m,板寬b=0.473,m,板間距2,h=0.02,m.入口計算域被劃分成101 101×的網(wǎng)格.分散相被認為以當(dāng)量流量[9]的形式均布在網(wǎng)格格點上.

    (2)如果上下板的溫度相同,那么,式(5)、式(6)分別與式(8)、式(9)相同,則整個計算過程相當(dāng)于在平均溫度情況下進行的,該計算比較容易;如果上下板的溫度不同,那么,式(5)、式(6)分別與式(8)、式(9)不同,則首先要取上下板溫度的平均值進行計算,得到的計算初值比較接近于實際計算得到的平均流剖面,將該平均溫度值代入式(10)—式(17)中就得到各個物性參數(shù)值,再由式(5)、式(6)(或式(8)、式(9))以及式(23)—式(27)計算出平均停留時間t、液滴沉降速度、平均流剖面以及壓力差Δp,變化溫度場的平均流修正將依賴于計算入口的初始值.

    (4)流場中液滴運動軌跡的計算方法是:將式(22)代入式(21)中,并將式(21)代入式(19)中,再利用式(20)消去ur形成只含有up和t的等式,然后,將式(19)沿x、y、z方向展開,積分并離散該方程得

    這里,dp0=1.5× 10?3m,該值是取值于常溫20,℃情況下的液滴直徑,ρd0為常溫20,℃情況下分散相液滴

    的密度.

    (6)將步驟(5)計算得到的流場參數(shù)及物性參數(shù)代回式(28),并注意到此時的速度值取值于上一步計算得到的速度值,然后重復(fù)步驟(4)和(5)進行計算,到t=t時迭代終止,就可以得到整個流場的分析數(shù)據(jù)了.另外,液滴在界面處運動的平均速度、相對速度和絕對速度以及液滴在界面處發(fā)生形狀的變化是本文研究必不可少的一部分,有關(guān)它們的定義方法和計算方法參見文獻[12].

    液滴所受不平衡剪力定義為液滴上表面所受到的來自流場的剪力與下表面所受到的剪力之差.它描述的是液滴在界面處的旋轉(zhuǎn)運動和液滴與層膜的相互作用以及液滴與上板的相互作用,包括兩種形式:液滴與層膜之間的連續(xù)相薄膜排出之后液滴所受到的不平衡剪力Fτ(命名為case1);液滴與層膜之間的連續(xù)相薄膜排出之前液滴所受到的不平衡剪力Fτ(命名為case2).有關(guān)它們的計算方法和公式表達參見文獻[12].

    阿基米德數(shù)是密度、液滴直徑和表面張力的綜合表達式,即它反映的是多個物理參數(shù)的綜合效果.阿基米德數(shù)的表達式為

    為了給理論反演計算和工程計算提供簡化的依據(jù)并驗證模型的正確性,本文將計算模擬的結(jié)果與Meon等[2]從實驗擬合得到的無量綱經(jīng)驗方程進行對比,該方程表達式如下:

    2 結(jié)果與討論

    從圖3可以看出:溫度對平均流剖面的影響是很顯著的,溫度高則平均流剖面的平均速度越高,表現(xiàn)為速度峰值越尖;反之,則平均速度越低,表現(xiàn)為速度峰值越平緩.

    從圖4可以看出:T0=20,℃、T1=20,℃和T0=22.5,℃、T1=22.5,℃兩種情況下的平均流剖面表現(xiàn)出類似的特征,即它們的曲線均是對稱的,平均溫度高的平均流剖面尖且其所包含的面積大于溫度低的平均流剖面所包含的面積;當(dāng)T0=20,℃、T1=25,℃時,平均流剖面的峰值處向右偏離中心線,而T0=25,℃、T1=20,℃時的平均流剖面的峰值處向左偏離中心線,而且上述兩種情況的平均流剖面均關(guān)于T0=22.5,℃、T1=22.5,℃情況下的平均流剖面的中心線對稱.可見,平均溫度不同和溫度分布方式不同都會引起平均流剖面發(fā)生相應(yīng)的變化,但它們的表現(xiàn)形式各不相同.

    從圖5可以看出:上下板溫度相同時流場的溫度剖面為水平直線,且平均溫度高的溫度剖面在溫度低的溫度剖面上方;T0=20,℃、T1=25,℃時的溫度剖面與T0=25,℃、T1=20,℃時的溫度剖面關(guān)于T0=22.5,℃且T1=22.5,℃情況下的溫度剖面對稱.因此,溫度剖面受到壁溫的影響也是比較顯著的,它與對應(yīng)的平均流剖面表現(xiàn)出類似的特征.同樣,流場溫度的變化和溫度分布的變化都會引起對應(yīng)的平均流剖面和溫度剖面的變化,該變化將會引起液滴運動軌跡的改變和聚結(jié)效率的顯著變化.

    圖6可視化了不同溫度對相同液滴運動軌跡的影響,該圖跟蹤了網(wǎng)格點初始位置為(0,-45,-45)液滴的運動軌跡,從頁面向里各運動軌跡對應(yīng)溫度分別為:T0=10,℃、T1=10,℃(實線),T0=20,℃、T1=20,℃(虛線),T0=20,℃、T1=25,℃(點劃線)、T0=30,

    ℃、T1=30,℃(點),T0=40,℃、T1=40,℃(雙點劃線).圖6表明,溫度不同液滴的運動軌跡也不同,溫度越高則液滴運動軌跡偏轉(zhuǎn)得越大,溫度不均衡也會引起液滴運動軌跡的偏轉(zhuǎn),這種變化的動力來源于溫度剖面和平均流剖面發(fā)揮的作用,它顯著地影響液滴的聚結(jié)效率,因此,適當(dāng)?shù)卣{(diào)節(jié)流場溫度可以改善流場的聚結(jié)條件.

    圖7描述的是不同溫度下瞬時流場的y–z平面速度矢量圖.從圖7(a)可以看出,均勻溫度影響下的流場的y–z速度矢量圖表現(xiàn)出具有一定的規(guī)律性.而從圖7(b)和圖7(c)可以看出,上下板溫度不同時流場的y–z速度矢量圖變得紊亂,這表明液滴的運動方式發(fā)生了一定的變化,也表明在層流流場中液滴的運動方式受溫度分布的影響具有一定的隨機性,但液滴運動方式的變化要通過溫度剖面和平均流剖面的作用才能改變,因此,溫度引起的液滴運動方式的隨機變化因素中還包含確定性因素.比較圖7(b)和圖7(c)可以看出,T0=20,℃、T1=25,℃時的y–z速度矢量圖比T0=25,℃、T1=20,℃時的y–z速度矢量圖密集,這表明由于上下板溫度不同誘發(fā)了液滴運動軌跡的變化,盡管有時這種變化不是很明顯,但卻足以對液滴的聚結(jié)產(chǎn)生明顯的影響.T0=20,℃、T1=25,℃時的流場聚結(jié)效率要高于T0=25,℃、T1=20,℃時的聚結(jié)效率,其原因在于溫度分布方式不同.T0=25,℃、T1=20,℃時的流場聚結(jié)效率要高于T0=20,℃、T1=20,℃時的聚結(jié)效率,其原因在于T0=25,℃、T1=20,℃時的平均溫度高于T0=20,℃、T1=20,℃時的平均溫度.可見,液滴聚結(jié)效率的高低在一定程度上依賴于流場的平均溫度,平均溫度是影響聚結(jié)效率的重要因素,但溫度分布方式也是影響聚結(jié)效率的不可忽視的因素.

    從圖8可以看出:當(dāng)T=20,℃、dp≤4× 10?3,m時,液滴的平均速度、相對速度和絕對速度隨液滴直徑的增大而增大,此時,液滴的聚結(jié)過程表現(xiàn)為液滴與界面的相互作用,該作用也隨液滴直徑的增大而增

    大;當(dāng)dp>4× 10?3,m時,液滴的相對速度快速減小,相應(yīng)地它的絕對速度會快速地接近平均速度,液滴的運動受到層膜速度影響比較大,液滴與層膜的聚結(jié)方式表現(xiàn)為液滴與上板直接作用,此時液滴直接與上板撞擊并受到上板的摩擦阻滯作用,液滴的聚結(jié)效率也將大幅度提高,因此,大液滴比小液滴聚結(jié)效率高,該結(jié)論從理論上證明了Meon等[2]關(guān)于液滴的聚結(jié)效率隨液滴直徑的增大而增大的結(jié)論;同時,圖8也表明液滴在界面處與層膜作用方式將因液滴直徑的不同而不同.

    另外,從圖8也可以看出,當(dāng)dp≤4× 10?3,m時,T=10,℃的液滴的平均速度、相對速度和絕對速度要小于T=20,℃和T=30,℃的液滴速度的對應(yīng)值,這表明較低平均溫度的流場中的液滴的活動度要低于較高平均溫度的流場中的液滴的活動度,較低平均溫度的流場中的液滴與較高平均溫度的流場中的液滴相比它與層膜相互作用程度較差,在低平均溫度的流場中只有更大的液滴才能直接與上板相互作用而產(chǎn)生聚結(jié).因此,較低平均溫度的流場中的大液滴與上板作用程度比較弱,液滴不易聚結(jié);較高平均溫度的流場中的液滴與上板作用程度比較強,容易聚結(jié).適當(dāng)?shù)靥岣吡鲌龅臏囟扔兄谔岣咭旱蔚木劢Y(jié)效率.

    從圖9可以看出:T=10、20、30,℃時的液滴直徑與case1情況的不平衡剪力曲線幾乎不隨溫度的變化而產(chǎn)生明顯變化,它們的case1情況的不平衡剪力曲線變化主要取決于液滴的直徑變化,液滴所受到的不平衡剪力表現(xiàn)為以case1情況的不平衡剪力為主要方式;當(dāng)dp≤2× 10?3,m時,case2情況的不平衡剪力曲線與case1情況的不平衡剪力曲線比較接近,case2情況的不平衡剪力曲線小于等于零,小液滴表現(xiàn)為順時針旋轉(zhuǎn);當(dāng)2× 10?3m<dp≤3× 10?3m時,case2情況的不平衡剪力曲線開始與case1情況的不平衡剪力曲線分開,液滴所受到的不平衡剪力變得更大,即液滴與層膜作用更加強烈了;當(dāng)3× 10?3m<dp≤4× 10?3m時,case2情況的不平衡剪力曲線快速攀升而轉(zhuǎn)為正值,此時,case2情況的不平衡剪力與case1情況的不平衡剪力在量級上比較接近,此區(qū)域為液滴旋轉(zhuǎn)運動不明朗區(qū)域,在該區(qū)域內(nèi)液滴與層膜之間的連續(xù)相薄膜仍然對液滴的運動方式產(chǎn)生強烈的影響,這表明在dp≤4× 10?3m區(qū)域內(nèi)液滴的聚結(jié)過程主要表現(xiàn)為液滴與層膜的直接作用;當(dāng)dp>4× 10?3m時,液滴所受到的不平衡剪力表現(xiàn)為以case1情況的不平衡剪力為主要方式,但液滴的聚結(jié)過程主要表現(xiàn)為液滴與上板的直接作用.以上結(jié)果從液滴所受到的不平衡剪力角度證明了對圖8分析得到的結(jié)論的正確性.

    另一方面,從圖9還可以看出:T=30,℃時的case2情況的不平衡剪力曲線隨液滴直徑增加變化顯著,T=20,℃時次之,T=10,℃時的case2情況的不平衡剪力曲線隨液滴直徑增加變化不明顯,這也表明較高的流場平均溫度有利于液滴與層膜或上板的相互作用.

    盡管從表面上看case1情況的不平衡剪力曲線幾乎不隨溫度的變化而產(chǎn)生明顯變化,但從case1情況的不平衡剪力曲線放大圖上可以看出,case1情況的不平衡剪力曲線隨溫度變化情況與case2情況的不平衡剪力曲線隨溫度變化情況表現(xiàn)出類似的特征,因此,溫度場的變化會改變液滴的受力方式.

    一般而言,液滴在上板的聚結(jié)效率依賴于液滴群在上板的分布率,液滴群在靠近流場計算入口處且到達上板的分布率越高則聚結(jié)效率越高.為了能從數(shù)值上明確地比較由于液滴在上板的隨機分布不同引起的聚結(jié)效率變化,本文采取了2種方法:截斷積分法和隨機概率期望值法.截斷積分法就是截取

    0≤x≤0.015這個區(qū)間并對圖10中的各個曲線進行積分,然后比較積分值以分析效率變化.隨機概率期望值法就是對圖10中的各個曲線上的隨機點求數(shù)學(xué)期望值的方法.這兩種方法相互結(jié)合使用才能有效地分析效率的變化.

    通過上述2種方法可以看出:平均溫度越高則液滴群在上板的分布率和對應(yīng)的聚結(jié)效率也越高,因此,圖10中T0=40,℃、T1=40,℃時的聚結(jié)效率最高;T0=20,℃且T1=25,℃時的液滴群的分布率要高于T0=25,℃且T1=20,℃時的情況,這表明液滴群的聚結(jié)效率在一定程度上也依賴于溫度分布方式.上下板溫度分布方式在一定程度上可用于改善或抑制流場的聚結(jié)條件,因此,適當(dāng)調(diào)整上下板溫度分布方式可以在一定情況下彌補由于裝備能力受限引起聚結(jié)效率的下降,這就再一次證明了上述結(jié)論的正確性.另一方面,可以看出溫度也是影響聚結(jié)效率和實驗效果及實驗可重復(fù)性的重要因素,這一點在過去的理論研究與實驗研究過程中經(jīng)常被忽視.需要注意的是,上述結(jié)論僅限于本文所取的溫度范圍和實驗范圍,因為更高的平均溫度情況下流場和液滴群的物理化學(xué)性質(zhì)可能發(fā)生顯著地變化,對于更高的平均溫度的情況有待于進一步的研究.

    從圖11可以看出:當(dāng)Ar≥400時,T=10、20、30,℃的阿基米德數(shù)與液滴雷諾數(shù)的對數(shù)曲線(Ar–Rep曲線)與Meon等[2]從實驗得到的經(jīng)驗方程(式(31))描述的曲線結(jié)果吻合較好;當(dāng)Ar<400時,T=20,℃的Ar–Rep曲線與經(jīng)驗方程(式(31))(實驗溫度為20,℃)所確定的曲線間存在誤差,誤差范圍為0.1<Rep<0.8.上述對比證明了本文數(shù)學(xué)模型的正確性.

    值得注意的是,當(dāng)Ar<400時,在相同的阿基米德數(shù)情況下,T=30,℃的液滴雷諾數(shù)要大于T=10,℃和20,℃的液滴雷諾數(shù),T=10,℃的液滴雷諾數(shù)最小,這表明較高的流場平均溫度可以提升小液滴活動程度,使其與層膜作用加強,較低的流場平均溫度使液滴與層膜作用減弱.

    另一方面,就正戊醇–水物系的Ar–Rep曲線而言,當(dāng)Ar≥400時有利于理論反演計算和簡化工程計算,因為此時的Ar–Rep曲線隨溫度變化不明顯;當(dāng)Ar<400時混合物系的Ar–Rep曲線受溫度影響較大,因此,不利于理論反演計算和簡化工程計算.然而,考慮到圖11中的4條曲線相差誤差不大,在一定程度上仍可以利用本文計算結(jié)果進行理論反演計算,但必須明確流場的平均溫度.

    3 結(jié) 論

    (1)流場溫度的變化和溫度分布的變化都會引起對應(yīng)的平均流剖面和溫度剖面的變化,該變化將會引起液滴運動軌跡的改變和聚結(jié)效率的顯著變化.這種變化的動力來源于溫度剖面和平均流剖面發(fā)揮的作用,因此,適當(dāng)?shù)卣{(diào)節(jié)流場溫度可以改善流場的聚結(jié)條件并彌補由于裝備能力受限引起聚結(jié)效率的下降.

    (2)在層流流場中液滴的運動方式受溫度分布的影響具有一定的隨機性,但隨機變化因素中還包含確定性因素.平均溫度是影響聚結(jié)效率的重要因素,而溫度分布方式也是不可忽視的因素.研究溫度場變化對聚結(jié)的影響對更加深入地研究調(diào)節(jié)流場的物理化學(xué)性質(zhì)以加速或減緩聚結(jié)效率的方法成為可能.參考文獻:

    [1] Mizrahi J,Barnea E. Compact settler gives efficient separation of liquid-liquid dispersions[J]. Process Engineering,1973,1(1):60–65.

    [2] Meon W,Blass E. Particle velocity in a stratified twophase flow between inclined parallel plates [J]. Chemical Engineering and Technology,1987,10:281–289.

    [3] Meon W. Plate separators for dispersed liquid-liquid systems:Hydrodynamic coalescence mode[J]. Chemical Engineering Science,1993,48(1):159–168.

    [4] Meon W,Blass E. Drop coalescence on inclined plates in liquids[J]. Chemical Engineering and Technology,1991,14(1):11–19.

    [5] Rommel W,Blass E,Meon W. Plate separators for dispersed liquid-liquid systems:Multiphase flow,droplet coalescence,separation performance and design[J]. Chemical Engineering Science,1992,47(3):555–564.

    [6] Goldman A J,Cox R G,Brenner H. Slow viscous motion of a sphere parallel to a plane wall(I):Motion through a quiescent fluid[J]. Chemical Engineering Science,1967,22(4):637–651.

    [7] Rowley M E,Davies G A. The design of plate separators for the separation of oil-water dispersions[J]. Chemical Engineering Research and Design,1988,66(4):313– 322.

    [8] 唐洪濤,吳劍華. 斜板沉降器的理論板長的研究[J].流體機械,2005,33(10):20–23.

    [9] Tang Hongtao,Wu Jianhua. Efficiency of plate separators for dispersed liquid-liquid systems[J]. Transactions of Tianjin University,2006,12(5):313–318.

    [10] Tang Hongtao,Wang Hui. Characteristic of drop coalescence resting on liquid-liquid interface[J]. Transactions of Tianjin University,2010,16(4):244–250.

    [11] Tang Hongtao,Chen Jianping,Cui Shihai. Coalescence behaviors of drop swarms on liquid-liquid interface[J]. Transactions of Tianjin University,2011,17(2):96–102.

    [12] 唐洪濤,王彥飛. 液滴聚結(jié)過程中隨機因素的分析[J].高?;瘜W(xué)工程學(xué)報,2012,26(6):964–970.

    [13] 唐洪濤,崔世海. 界面物性對液滴聚結(jié)的影響[J]. 化工學(xué)報,2012,63(4):1140–1148.

    [14] 鄢國森,Карапетьянц М Х. 正戊醇及正庚醇的黏度及密度與溫度的關(guān)系[J]. 四川大學(xué)學(xué)報,1959,4:51–53.

    [15] 中國石化集團上海工程有限公司. 化工工藝設(shè)計手冊[M]. 4版. 北京:化學(xué)工業(yè)出版社,2009.

    責(zé)任編輯:常濤

    Temperature Deflection Effect on Drop Trajectory in the Flowing Field

    TANG Hongtao1,GENG Shengmin1,MENG Huibo2
    (1. College of Mechanical Engineering,Tianjin University of Science & Technology,Tianjin 300222,China;2. School of Mechanical Engineering,Shenyang University of Chemical Technology,Shenyang 110142,China)

    A trajectory model of dispersed phase drops and A modified model of basic flow in the temperature field between two inclined parallel plates were derived with the optimized calculation. The analytical results of direct numerical simulation indicate that the difference of both the mean temperature and temperature distribution may lead to changes of mean flow profiles and temperature profiles,but their functions are very different. The trajectories of a drop may vary with the changes of unbalanced shearing force acting on a drop induced by temperature field. The dynamic sources of changeable processes stem from the action of mean flow profiles and temperature profiles,which have a great effect on the coalescence efficiencies. Hence,the coalescence conditions of the flowing field may be improved by the proper adjustment of temperature in the flowing field. Although the flowing field is in the laminar flow state,temperature distribution leads to the fact that in the change of drop motion,there are not only stochastic factors but also determinate factors. The difference in the degree of coalescence efficiency depends on the mean temperature to a certain extent,but the temperature distributions are indispensable to coalescence efficiency. The factors of temperature field are indispensable to sovling the theoretical simplification of the inverse operation and engineering calculation.

    temperature;coalescence;drop;mean flow profile;thermal conductivity

    TB126;TQ021.1

    A

    1672-6510(2014)06-0065-09

    10.13364/j.issn.1672-6510.2014.06.013

    2014–02–05;

    2014–07–07

    天津市自然科學(xué)基金資助項目(09JCYBJC06400);國家自然科學(xué)基金資助項目(21106086)

    唐洪濤(1968—),男,黑龍江海林人,副教授,tanghongtao@tust.edu.cn.

    聚結(jié)理論研究的是自由流場界面的特性和分層體系之間相互作用的關(guān)系[1–2].截至目前,聚結(jié)理論還沒被深入細致地理解和掌握,其原因在于它是一種復(fù)雜的物理化學(xué)過程.因此,聚結(jié)理論還需要引入很多物理化學(xué)原理才能解釋聚結(jié)發(fā)生的細節(jié).

    液滴以什么樣的運動方式在分層二相流體系中運動,是研究聚結(jié)理論領(lǐng)域普遍關(guān)心的問題.Meon等[3–5]的研究表明:在某些限定物理條件下,聚結(jié)不

    猜你喜歡
    平均溫度液滴剪力
    3月熱浪來襲悉尼或迎165年以來新紀(jì)錄
    南方地區(qū)圓拱形和鋸齒形大棚內(nèi)溫度四季差別探究*
    液滴間相互碰撞融合與破碎的實驗研究
    噴淋液滴在空氣環(huán)境下的運動特性
    懸臂箱形截面梁的負剪力滯效應(yīng)
    考慮截面配筋的箱梁剪力滯效應(yīng)分析
    云南保山氣溫變化特征及其均生函數(shù)預(yù)測實驗
    徐州地區(qū)加權(quán)平均溫度模型研究
    箱型梁剪力滯效應(yīng)的解耦求解
    80m簡支鋼箱梁剪力滯系數(shù)研究
    一区二区三区激情视频| av视频在线观看入口| 变态另类成人亚洲欧美熟女| 极品教师在线视频| 嫩草影院入口| 国产一区二区在线观看日韩| 网址你懂的国产日韩在线| 高清在线国产一区| 最近在线观看免费完整版| 国产精品精品国产色婷婷| 无遮挡黄片免费观看| 12—13女人毛片做爰片一| 精品一区二区三区人妻视频| 日韩精品青青久久久久久| 啦啦啦观看免费观看视频高清| 熟女人妻精品中文字幕| 偷拍熟女少妇极品色| 日韩欧美一区二区三区在线观看| 成人av一区二区三区在线看| 午夜日韩欧美国产| 久久久久久久久中文| 亚洲va日本ⅴa欧美va伊人久久| 国产精品人妻久久久久久| 在线播放无遮挡| 亚洲男人的天堂狠狠| 国产69精品久久久久777片| 麻豆av噜噜一区二区三区| 99riav亚洲国产免费| 国产一区二区三区在线臀色熟女| 久久久久久久久大av| 国内少妇人妻偷人精品xxx网站| 一级黄色大片毛片| 国产亚洲精品av在线| 国产伦在线观看视频一区| 国产精品不卡视频一区二区 | 婷婷亚洲欧美| 久久久精品大字幕| 很黄的视频免费| 一边摸一边抽搐一进一小说| 中国美女看黄片| 简卡轻食公司| 国产精品伦人一区二区| 色综合欧美亚洲国产小说| 很黄的视频免费| 狠狠狠狠99中文字幕| 国内揄拍国产精品人妻在线| 搡女人真爽免费视频火全软件 | 一个人免费在线观看的高清视频| 我要看日韩黄色一级片| 亚洲av熟女| 一个人观看的视频www高清免费观看| 一边摸一边抽搐一进一小说| 久99久视频精品免费| 99久久精品国产亚洲精品| 国产69精品久久久久777片| 18美女黄网站色大片免费观看| 欧美性猛交黑人性爽| 国产主播在线观看一区二区| 一区二区三区高清视频在线| 不卡一级毛片| 一级黄色大片毛片| 亚洲狠狠婷婷综合久久图片| 嫩草影视91久久| 又黄又爽又免费观看的视频| 窝窝影院91人妻| 成人鲁丝片一二三区免费| 国产亚洲av嫩草精品影院| 黄色视频,在线免费观看| 亚洲av.av天堂| 亚洲欧美日韩无卡精品| 在线观看66精品国产| 亚洲成av人片在线播放无| 观看免费一级毛片| 麻豆av噜噜一区二区三区| 久久九九热精品免费| 女人十人毛片免费观看3o分钟| 少妇的逼好多水| 色在线成人网| 麻豆久久精品国产亚洲av| 又黄又爽又免费观看的视频| 99精品在免费线老司机午夜| 色尼玛亚洲综合影院| 国产成人av教育| 毛片女人毛片| 观看美女的网站| 十八禁网站免费在线| 亚洲美女视频黄频| 午夜福利欧美成人| 身体一侧抽搐| 看十八女毛片水多多多| 亚洲成av人片在线播放无| 赤兔流量卡办理| 九色国产91popny在线| 久久久国产成人精品二区| a在线观看视频网站| 国产精品久久视频播放| 国产一区二区亚洲精品在线观看| 黄片小视频在线播放| 3wmmmm亚洲av在线观看| 人人妻人人澡欧美一区二区| 三级国产精品欧美在线观看| 看十八女毛片水多多多| 欧美乱色亚洲激情| 国产色婷婷99| 尤物成人国产欧美一区二区三区| 亚洲欧美精品综合久久99| 亚洲人成伊人成综合网2020| 亚洲人与动物交配视频| 久久久色成人| 18禁裸乳无遮挡免费网站照片| 女人十人毛片免费观看3o分钟| 美女xxoo啪啪120秒动态图 | 老司机福利观看| 久久99热6这里只有精品| 午夜激情福利司机影院| 国产美女午夜福利| 亚洲国产色片| 性色av乱码一区二区三区2| 人人妻,人人澡人人爽秒播| 国产亚洲精品av在线| 中文字幕人妻熟人妻熟丝袜美| 精品日产1卡2卡| 亚洲欧美激情综合另类| aaaaa片日本免费| 久久九九热精品免费| 国产亚洲精品久久久com| 久久久久国内视频| 日韩av在线大香蕉| 男女视频在线观看网站免费| 国产伦人伦偷精品视频| 我要看日韩黄色一级片| 嫩草影院入口| 成人午夜高清在线视频| 国产成人影院久久av| 国产精品久久久久久人妻精品电影| 乱人视频在线观看| 丰满的人妻完整版| 精品久久久久久久久亚洲 | 国内精品美女久久久久久| 免费无遮挡裸体视频| 成人高潮视频无遮挡免费网站| 欧美高清成人免费视频www| 啦啦啦韩国在线观看视频| 九色国产91popny在线| www.熟女人妻精品国产| 国产 一区 欧美 日韩| 亚洲久久久久久中文字幕| 国产单亲对白刺激| 美女 人体艺术 gogo| 中文字幕免费在线视频6| 日本黄色片子视频| 欧美中文日本在线观看视频| 亚洲美女搞黄在线观看 | 露出奶头的视频| 麻豆国产97在线/欧美| 国产av不卡久久| 长腿黑丝高跟| 哪里可以看免费的av片| 综合色av麻豆| 噜噜噜噜噜久久久久久91| 九九热线精品视视频播放| 久久久久久久午夜电影| 亚洲熟妇熟女久久| 舔av片在线| 偷拍熟女少妇极品色| 免费观看人在逋| 国产精品女同一区二区软件 | 欧美性猛交黑人性爽| 露出奶头的视频| 黄色视频,在线免费观看| 99国产综合亚洲精品| 免费看美女性在线毛片视频| 国产伦人伦偷精品视频| 国产精品爽爽va在线观看网站| 久久精品综合一区二区三区| 有码 亚洲区| 91在线精品国自产拍蜜月| 欧美成人性av电影在线观看| 91在线精品国自产拍蜜月| 少妇裸体淫交视频免费看高清| 日本免费a在线| 亚洲av第一区精品v没综合| 狠狠狠狠99中文字幕| 男人的好看免费观看在线视频| www.www免费av| 午夜免费激情av| 午夜精品久久久久久毛片777| 18禁黄网站禁片午夜丰满| 午夜精品久久久久久毛片777| www.www免费av| 在现免费观看毛片| 国产美女午夜福利| 国产日本99.免费观看| 夜夜夜夜夜久久久久| 欧美日韩亚洲国产一区二区在线观看| 久久国产精品人妻蜜桃| 国产色婷婷99| netflix在线观看网站| 国产免费男女视频| 黄色配什么色好看| 在线a可以看的网站| 变态另类成人亚洲欧美熟女| 夜夜看夜夜爽夜夜摸| 每晚都被弄得嗷嗷叫到高潮| 欧美又色又爽又黄视频| 搡女人真爽免费视频火全软件 | 欧美zozozo另类| 性插视频无遮挡在线免费观看| 能在线免费观看的黄片| 免费在线观看成人毛片| 国产精品伦人一区二区| 又爽又黄a免费视频| 日本一本二区三区精品| 国产伦精品一区二区三区四那| 亚洲一区二区三区不卡视频| 黄色一级大片看看| 国产老妇女一区| 欧美一区二区精品小视频在线| 2021天堂中文幕一二区在线观| 少妇被粗大猛烈的视频| 色综合亚洲欧美另类图片| 免费高清视频大片| 亚洲性夜色夜夜综合| 国产精品久久久久久人妻精品电影| 午夜福利视频1000在线观看| 中文字幕精品亚洲无线码一区| 99久国产av精品| 国产成人啪精品午夜网站| 九色国产91popny在线| 亚洲国产日韩欧美精品在线观看| 国产精品国产高清国产av| 欧美日韩瑟瑟在线播放| 欧美三级亚洲精品| 亚洲 国产 在线| 男人舔女人下体高潮全视频| 十八禁网站免费在线| 国产一区二区三区在线臀色熟女| 男女之事视频高清在线观看| 亚洲精华国产精华精| bbb黄色大片| 亚洲av不卡在线观看| 日韩大尺度精品在线看网址| 午夜视频国产福利| 欧美日韩瑟瑟在线播放| 91久久精品电影网| 一区福利在线观看| www.999成人在线观看| 亚洲av熟女| 最近在线观看免费完整版| 久久99热6这里只有精品| 久久久久精品国产欧美久久久| 99久久精品一区二区三区| 一个人看视频在线观看www免费| 亚洲真实伦在线观看| 嫩草影院精品99| 免费高清视频大片| 一个人免费在线观看的高清视频| 99久久无色码亚洲精品果冻| 好男人在线观看高清免费视频| 狂野欧美白嫩少妇大欣赏| 国产v大片淫在线免费观看| 午夜久久久久精精品| 亚洲成av人片在线播放无| 国产午夜福利久久久久久| 好男人电影高清在线观看| eeuss影院久久| 亚洲熟妇熟女久久| 69av精品久久久久久| 免费看日本二区| 午夜精品一区二区三区免费看| 欧美激情国产日韩精品一区| 国产午夜福利久久久久久| 久久国产乱子免费精品| 亚洲国产欧洲综合997久久,| 色播亚洲综合网| 免费在线观看亚洲国产| 国产v大片淫在线免费观看| 色综合婷婷激情| 日韩人妻高清精品专区| 国产精品,欧美在线| 欧美zozozo另类| 美女 人体艺术 gogo| 老女人水多毛片| 国产探花极品一区二区| 三级国产精品欧美在线观看| 男人舔女人下体高潮全视频| 99国产综合亚洲精品| 精品欧美国产一区二区三| 亚洲第一欧美日韩一区二区三区| 免费观看精品视频网站| 国语自产精品视频在线第100页| 嫩草影视91久久| 国产精品一区二区三区四区久久| 欧美在线一区亚洲| 国产美女午夜福利| 老女人水多毛片| 欧美绝顶高潮抽搐喷水| 色哟哟哟哟哟哟| 国内揄拍国产精品人妻在线| 美女xxoo啪啪120秒动态图 | 黄色一级大片看看| 男人舔奶头视频| 黄片小视频在线播放| 丝袜美腿在线中文| 热99re8久久精品国产| 久久国产精品人妻蜜桃| 亚洲国产精品合色在线| 成人av一区二区三区在线看| 色哟哟·www| 一区二区三区免费毛片| 首页视频小说图片口味搜索| 可以在线观看毛片的网站| 国产老妇女一区| eeuss影院久久| 超碰av人人做人人爽久久| 好男人在线观看高清免费视频| 亚洲久久久久久中文字幕| 亚洲国产日韩欧美精品在线观看| 日本黄大片高清| 精品一区二区免费观看| 搡老妇女老女人老熟妇| 亚洲人成电影免费在线| 久久人妻av系列| 国产乱人伦免费视频| 一区福利在线观看| 免费看美女性在线毛片视频| 可以在线观看的亚洲视频| av在线老鸭窝| 亚洲色图av天堂| 精品久久久久久,| 亚洲av日韩精品久久久久久密| 嫁个100分男人电影在线观看| 精品99又大又爽又粗少妇毛片 | 看免费av毛片| 身体一侧抽搐| 高清日韩中文字幕在线| 欧美成人性av电影在线观看| 久99久视频精品免费| 精品一区二区三区av网在线观看| 婷婷精品国产亚洲av| 亚洲成人中文字幕在线播放| a级毛片a级免费在线| 麻豆成人午夜福利视频| 天堂影院成人在线观看| 国语自产精品视频在线第100页| 五月伊人婷婷丁香| 一个人免费在线观看的高清视频| 午夜精品在线福利| 老司机福利观看| 国产探花极品一区二区| 国产精品影院久久| 欧美性猛交黑人性爽| 久99久视频精品免费| 怎么达到女性高潮| 国产亚洲精品久久久久久毛片| 精品国产三级普通话版| 久久久久久久精品吃奶| 日韩欧美在线二视频| 国产精品乱码一区二三区的特点| 亚洲av成人不卡在线观看播放网| 国产精品免费一区二区三区在线| 亚洲aⅴ乱码一区二区在线播放| 夜夜夜夜夜久久久久| 国产亚洲精品av在线| 亚洲av成人不卡在线观看播放网| 国产黄a三级三级三级人| 嫁个100分男人电影在线观看| 日本在线视频免费播放| 成人高潮视频无遮挡免费网站| 精华霜和精华液先用哪个| 一夜夜www| 如何舔出高潮| 国内精品一区二区在线观看| 五月伊人婷婷丁香| 在线观看美女被高潮喷水网站 | 亚洲国产日韩欧美精品在线观看| 国产精品女同一区二区软件 | 毛片女人毛片| 久久国产精品人妻蜜桃| 熟女电影av网| 有码 亚洲区| 国产成人aa在线观看| 国产激情偷乱视频一区二区| 国产免费一级a男人的天堂| 天堂动漫精品| 色播亚洲综合网| 中文字幕人妻熟人妻熟丝袜美| 熟女电影av网| 波多野结衣高清作品| av女优亚洲男人天堂| a级一级毛片免费在线观看| 好男人电影高清在线观看| 岛国在线免费视频观看| 校园春色视频在线观看| 亚洲av五月六月丁香网| 精品一区二区免费观看| 桃红色精品国产亚洲av| 久久久久久久亚洲中文字幕 | 禁无遮挡网站| 国产精品电影一区二区三区| 欧美+亚洲+日韩+国产| 俄罗斯特黄特色一大片| 欧美中文日本在线观看视频| 又黄又爽又免费观看的视频| 男女床上黄色一级片免费看| 能在线免费观看的黄片| 搞女人的毛片| 久久精品国产亚洲av香蕉五月| 在线观看午夜福利视频| 夜夜躁狠狠躁天天躁| 又粗又爽又猛毛片免费看| 欧美激情在线99| 免费黄频网站在线观看国产| 1000部很黄的大片| av一本久久久久| 麻豆精品久久久久久蜜桃| 精品酒店卫生间| 丰满人妻一区二区三区视频av| 国产精品99久久99久久久不卡 | 五月玫瑰六月丁香| 国产精品国产三级国产专区5o| 国产成人一区二区在线| 在线观看一区二区三区| 天天躁日日操中文字幕| av在线老鸭窝| 男女啪啪激烈高潮av片| 午夜亚洲福利在线播放| 国产 一区 欧美 日韩| 中文天堂在线官网| 国产精品女同一区二区软件| 欧美成人午夜免费资源| 成年女人看的毛片在线观看| 97在线视频观看| 极品教师在线视频| 亚洲av欧美aⅴ国产| 青青草视频在线视频观看| 亚洲美女搞黄在线观看| 国内精品美女久久久久久| 亚洲av福利一区| 久久久久性生活片| 亚洲人成网站在线观看播放| 日韩欧美 国产精品| 99精国产麻豆久久婷婷| 91aial.com中文字幕在线观看| 在线观看免费高清a一片| 欧美极品一区二区三区四区| 熟女电影av网| av又黄又爽大尺度在线免费看| 国产黄色视频一区二区在线观看| 欧美丝袜亚洲另类| 亚洲精品久久午夜乱码| 日韩免费高清中文字幕av| 少妇的逼水好多| 神马国产精品三级电影在线观看| 国产av国产精品国产| 午夜精品国产一区二区电影 | 日韩在线高清观看一区二区三区| 亚洲av欧美aⅴ国产| 黄色怎么调成土黄色| 青春草国产在线视频| 精品视频人人做人人爽| 亚洲第一区二区三区不卡| 亚洲国产av新网站| 免费黄色在线免费观看| 波野结衣二区三区在线| 亚洲国产av新网站| 午夜激情福利司机影院| 熟女av电影| 久久久久久国产a免费观看| 久久精品久久久久久久性| 永久网站在线| av一本久久久久| 超碰97精品在线观看| 国产毛片a区久久久久| 国产男女内射视频| 人妻系列 视频| 两个人的视频大全免费| 亚洲精品国产成人久久av| 亚洲国产精品国产精品| av.在线天堂| 国产亚洲最大av| 中文字幕制服av| 午夜福利在线观看免费完整高清在| 天美传媒精品一区二区| 蜜桃久久精品国产亚洲av| 日产精品乱码卡一卡2卡三| 女人十人毛片免费观看3o分钟| 午夜精品一区二区三区免费看| 久久久久久九九精品二区国产| 久久6这里有精品| videos熟女内射| 边亲边吃奶的免费视频| 中文精品一卡2卡3卡4更新| 天堂俺去俺来也www色官网| 久久久久久久久久久免费av| av一本久久久久| 国产91av在线免费观看| 欧美高清成人免费视频www| 欧美最新免费一区二区三区| 丰满人妻一区二区三区视频av| 亚洲国产精品成人久久小说| 亚洲精华国产精华液的使用体验| 亚洲欧美成人综合另类久久久| 久久国内精品自在自线图片| 成人高潮视频无遮挡免费网站| 亚洲av免费在线观看| 中文字幕亚洲精品专区| 精品熟女少妇av免费看| 午夜免费鲁丝| 国产亚洲av片在线观看秒播厂| 日韩欧美精品v在线| 一区二区三区四区激情视频| 尤物成人国产欧美一区二区三区| 国产一区二区三区av在线| 久久国内精品自在自线图片| 国产成人精品福利久久| 97超视频在线观看视频| 亚洲av不卡在线观看| 亚洲电影在线观看av| 国精品久久久久久国模美| 在线观看人妻少妇| 看非洲黑人一级黄片| 亚洲精品日本国产第一区| 久久久久久久久久人人人人人人| 精品午夜福利在线看| 国产久久久一区二区三区| 麻豆久久精品国产亚洲av| 男男h啪啪无遮挡| 午夜精品一区二区三区免费看| www.色视频.com| 性色avwww在线观看| 久久鲁丝午夜福利片| 国产男人的电影天堂91| 亚洲图色成人| 亚洲av欧美aⅴ国产| 99久久精品国产国产毛片| 亚洲精品一二三| 国产男人的电影天堂91| 免费看日本二区| 日日啪夜夜撸| 六月丁香七月| www.av在线官网国产| 啦啦啦在线观看免费高清www| 高清在线视频一区二区三区| 午夜福利在线在线| 狂野欧美激情性bbbbbb| 亚洲最大成人中文| 久久人人爽人人爽人人片va| 日日啪夜夜爽| 亚洲成人精品中文字幕电影| 久久久久性生活片| 熟女电影av网| 亚洲va在线va天堂va国产| 国产伦理片在线播放av一区| 噜噜噜噜噜久久久久久91| 熟妇人妻不卡中文字幕| 亚洲欧美成人综合另类久久久| 欧美xxxx性猛交bbbb| 精品人妻偷拍中文字幕| 亚洲精品国产av成人精品| 大陆偷拍与自拍| 精品国产乱码久久久久久小说| 热re99久久精品国产66热6| 伦理电影大哥的女人| 国产黄片视频在线免费观看| 欧美高清成人免费视频www| 99热网站在线观看| 欧美极品一区二区三区四区| 国产高潮美女av| 亚洲av免费高清在线观看| 国产探花在线观看一区二区| 国产欧美另类精品又又久久亚洲欧美| 建设人人有责人人尽责人人享有的 | 99热这里只有精品一区| 亚洲欧美精品自产自拍| 国产在视频线精品| 国产免费一级a男人的天堂| 午夜福利视频1000在线观看| 亚洲精品aⅴ在线观看| 观看美女的网站| 国产男女超爽视频在线观看| 亚洲va在线va天堂va国产| 免费看光身美女| 最后的刺客免费高清国语| 三级经典国产精品| 在线免费十八禁| 日本与韩国留学比较| 国产精品熟女久久久久浪| 亚洲激情五月婷婷啪啪| 亚洲婷婷狠狠爱综合网| 1000部很黄的大片| 国产 精品1| 国内精品美女久久久久久| 国产人妻一区二区三区在| 美女脱内裤让男人舔精品视频| 一本一本综合久久| 22中文网久久字幕| 在线 av 中文字幕| 国产永久视频网站| 国内精品宾馆在线| 国产高潮美女av| 婷婷色综合www| 成人毛片a级毛片在线播放| 狂野欧美白嫩少妇大欣赏| 日韩av不卡免费在线播放| 如何舔出高潮| 免费av毛片视频| 伊人久久国产一区二区| 777米奇影视久久| 老女人水多毛片| 国产男人的电影天堂91| 国产亚洲av片在线观看秒播厂| 中文精品一卡2卡3卡4更新| 亚洲丝袜综合中文字幕| 日韩一本色道免费dvd|