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

    水平表面氣流剪切作用下的水膜厚度

    2017-11-22 01:28:33冷夢堯常士楠丁亮李曉峰
    航空學(xué)報 2017年2期
    關(guān)鍵詞:水膜雷諾數(shù)氣流

    冷夢堯, 常士楠, 丁亮, 李曉峰

    北京航空航天大學(xué) 航空科學(xué)與工程學(xué)院, 北京 100083

    水平表面氣流剪切作用下的水膜厚度

    冷夢堯, 常士楠*, 丁亮, 李曉峰

    北京航空航天大學(xué) 航空科學(xué)與工程學(xué)院, 北京 100083

    飛機結(jié)冰表面上的液態(tài)水受氣流吹拂作用會發(fā)生向后溢流,從而影響結(jié)冰區(qū)域范圍及防冰系統(tǒng)設(shè)計;為了獲得水膜流動規(guī)律,對水平平板表面上氣流剪切驅(qū)動的水膜流動進行了實驗測量和建模分析。通過水膜流動風(fēng)洞試驗臺產(chǎn)生高速氣流驅(qū)動水膜的流動,使用色散共焦位移計測量同一位置的水膜在不同時刻的厚度變化,結(jié)果表明氣-液界面由底層薄水膜和多種尺度的波動組成,具有變化速度快隨機性強的特點。通過水膜厚度隨氣流速度及水膜雷諾數(shù)的變化規(guī)律,發(fā)現(xiàn)平均水膜厚度與兩者均呈現(xiàn)出單調(diào)非線性的依賴關(guān)系?;诒∷ち鲃永碚摵推骄ず穸葘嶒灲Y(jié)果,提出了高速氣流剪切作用下的氣-液波動界面剪切因子計算式,適用于風(fēng)速17.8~52.2 m/s,水膜雷諾數(shù)26~128之間的平板水膜流動計算。

    飛機防冰; 兩相流; 水膜厚度; 界面剪切力; 建模

    當(dāng)飛機穿過含有過冷水滴的云層時,其迎風(fēng)表面會迅速覆蓋積冰,嚴重威脅飛行安全。在不完全凍結(jié)或濕態(tài)防冰過程中,表面存留的液態(tài)水會在氣流驅(qū)動下向后流動,改變飛機表面水的分布范圍,形成冰角或溢流冰脊;水膜表面上的波動反過來會干擾氣流邊界層的流動,對結(jié)冰或防冰表面的熱量及其分布造成影響。而水膜厚度方向上的導(dǎo)熱對于結(jié)冰預(yù)測或防冰熱流計算也是必不可少的一項。故研究飛機表面水膜流動對于冰形模擬及防冰系統(tǒng)優(yōu)化設(shè)計是十分有意義的。

    在理論和數(shù)值仿真方面,經(jīng)典的Messinger模型[1]首先在飛機結(jié)冰預(yù)測中考慮液態(tài)水的存在,但該模型假設(shè)當(dāng)前控制體內(nèi)未凍結(jié)的水即全部進入下一個控制體。Ai-Khalil等[2]則通過對水膜的穩(wěn)定性分析發(fā)展了防冰表面溢流模型。Myers和Thompson[3-4]基于潤滑近似理論建立水膜流動方程,提高了水膜模型的物理意義,隨后該模型進一步被應(yīng)用于結(jié)冰理論[5-6]、防冰及溢流冰的模擬[7-8]中。Fortin等[9]則根據(jù)機翼表面的近距離攝像觀察結(jié)果,根據(jù)粗糙表面換熱模型研究液態(tài)水帶來的影響。還有部分研究者專注于水膜表面波動的形成和發(fā)展[10-12]。然而這些模型大都是同宏觀結(jié)果進行對比驗證,對于水膜流動過程中的部分物理機制缺乏了解,導(dǎo)致預(yù)測結(jié)果仍有一定誤差。

    近幾年,隨著實驗觀測技術(shù)的發(fā)展,研究飛機表面液態(tài)水流動機理的實驗逐漸增多,例如加速流場中的水珠脫離臨界條件[13],水珠受力情況及運動規(guī)律[14],溪流厚度的空間特征研究[15]等,這些實驗結(jié)果完善了飛機表面水珠及溪流形態(tài)的液態(tài)水流動機制。關(guān)于連續(xù)性水膜的研究,F(xiàn)eo和Tsao[16]為建立明冰縮比理論測量了駐點附近的水膜厚度,并建立其同模型尺寸、液態(tài)水含量、水滴直徑以及風(fēng)速之間的關(guān)系。Muzik等[17]觀察了水膜在翼型表面上的流動特征及破碎現(xiàn)象。Zhang 等[18]使用數(shù)字圖像投影技術(shù)(Digital Image Projection, DIP)測量了翼型表面的水膜厚度,得到水膜厚度約在10~1 000 μm之間,但未建立起有效的水膜厚度與外部條件之間的關(guān)系。

    本文針對高速氣流驅(qū)動的水膜流動進行了測量,并根據(jù)實驗數(shù)據(jù)建立起水膜流動方程。通過搭建特定的水膜觀測實驗臺在平板表面產(chǎn)生連續(xù)性流動的水膜,使用色散共焦位移計記錄一段時間內(nèi)水膜厚度的時域信號,分析水膜形態(tài)隨風(fēng)速及水流量的變化規(guī)律。根據(jù)平均水膜厚度實驗結(jié)果,結(jié)合由Myers模型推導(dǎo)的水平穩(wěn)態(tài)膜厚計算式,擬合出氣-液波動界面剪切系數(shù)計算式,并通過與無量綱水膜厚度文獻結(jié)果對比進行驗證。

    1 水膜流動方程

    已有實驗結(jié)果表明[5,18],飛機表面存在的水膜厚度極薄且流速較低,故而可將水膜簡化為二維不可壓縮層流流動。又由于水膜溫度變化區(qū)間較窄,認為水膜物性參數(shù)(密度和黏度)為常數(shù),假設(shè)表面外法向坐標為n,切向坐標為s,法向速度為v,切向速度為u,水膜連續(xù)性方程及二維不可壓縮穩(wěn)態(tài)Navier-Stokes方程可表達為

    (1)

    (2)

    (3)

    式中:ρw為水的密度;μw為水的黏度;p為水膜內(nèi)壓強;g為重力加速度。由于水膜厚度方向上的尺寸為微米級,相對于流向長度而言,可以忽略不計,即n?s??蓱?yīng)用潤滑近似理論[4]對動量方程考慮以下簡化:① 忽略表面法線方向速度,即認為v=0;② 結(jié)合連續(xù)性方程可忽略速度u的切向變化,即?u/?s=0。簡化后得到動量方程分別為

    (4)

    (5)

    根據(jù)式(5)可知,水膜內(nèi)壓強的法向梯度是由重力造成的,壓強p在水膜厚度方向上成線性分布。分別在氣-液交界面和液-固邊界應(yīng)用以下邊界條件[4]:① 壓強和界面切應(yīng)力連續(xù);② 無滑移平整邊界,對式(4)和式(5)先后求積分可得

    p=p∞(s)+ρwgn(n-h)

    (6)

    (7)

    式中:τ為界面剪切應(yīng)力;h為當(dāng)?shù)厮ず穸?;p∞(s)為氣流邊界層底部的壓強,當(dāng)氣流與水膜流動均處于充分發(fā)展階段,空氣壓力沿切向變化很小[6],可認為?p/?s~0。將速度沿法向積分可求得水流量的計算式為

    (8)

    為便于得到水膜表面剪切力的影響規(guī)律,本實驗關(guān)注在水平平板表面上的水膜流動,重力切向分量ρgs等于零,氣流剪切力成為水膜流動的關(guān)鍵影響因素,得到相應(yīng)的水膜厚度為

    (9)

    分析式(9)可發(fā)現(xiàn),液體黏度由水膜的溫度計算得到,而水流量Q取決于局部撞擊水量及上游流入的水量,故而氣-液界面剪切力τ的大小成為了求解水膜厚度的關(guān)鍵因素。Alzaili和Hammond[5]采用流場求解得到的壁面剪切應(yīng)力作為τ的值,但是顯然氣-液界面上的波動會大大增強兩相之間的相互作用。借鑒氣-液兩相分層流研究的經(jīng)驗,界面剪切力可以表示為[10,19]

    (10)

    式中:fi為界面剪切系數(shù);ρa為空氣密度;Ua為空氣主流速度。界面剪切系數(shù)fi通過實驗結(jié)果測得,并可表達為水膜雷諾數(shù)的函數(shù),即fi=f(Ref)。水膜雷諾數(shù)定義為Ref=Q/Wνw,W為水膜濕潤寬度,νw為液體運動黏度。而C則為氣-液界面波的相位速度。一方面,由于界面波動的速度遠小于主氣流速度,故而在式(10)中往往將其忽略[20-22]。另一方面,氣-液界面波的相位速度同樣取決于風(fēng)速Ua和水膜雷諾數(shù)Ref[23],即C=f(Ua,Ref)。本文在做實驗數(shù)據(jù)處理時亦采用忽略波動速度C的簡化式。

    2 試驗裝置

    水膜厚度測量在北京航空航天大學(xué)防/除冰系統(tǒng)研究實驗室的小型風(fēng)洞中進行,如圖1所示。高速氣流由離心風(fēng)機提供,經(jīng)穩(wěn)定段整流、收縮段加速后進入實驗段中。試驗段承力結(jié)構(gòu)為不銹鋼,底部中間安裝拋光鋁板作為實驗底板,表面粗糙度平均為0.19 μm,平衡接觸角為70°。實驗段頂部及左右側(cè)面使用透明亞克力玻璃板粘接制成,其中頂部中間段設(shè)計可拆卸,方便觀測儀器的安裝和調(diào)節(jié)。實驗段頂部距底板50 mm,橫向?qū)?00 mm,總長0.5 m。風(fēng)速由皮托管風(fēng)速計在位移計測量點附近測得,主氣流的速度范圍為0~55 m/s,橫截面上風(fēng)速均勻區(qū)可達76%。

    觀測所使用的水源為純凈水,由離心水泵驅(qū)動,經(jīng)由供水箱從底部進入試驗段。供水箱使用疏水性的有機玻璃材料制成,有助于穩(wěn)定水流,并增強水箱口的密封性。供水口為距試驗段入口約55 mm的長縫,且其下游側(cè)面邊緣倒有圓角。在試驗底板末端留有出水口,以吸水性材料填充,回收的這部分水流會再進入循環(huán)中。此外,高風(fēng)速大流量試驗中會有小水滴被卷吸脫離水膜表面,并隨氣流從風(fēng)洞擴壓段的出口離開,而這部分水不再進入水流循環(huán)。

    水膜的厚度使用色散共焦位移計(ACR-HNDS100)測量和記錄。該位移計由一個控制器和一支光筆組成,通過感測從物體表面反射回來的單色光波長,將波長換算獲得相應(yīng)距離值,若反射回兩個波長即可得到被測物的厚度值。位移計被安裝在水流出口下游245 mm的位置,試驗過程中與被測物之間的距離保持在其焦距(40 mm)附近,數(shù)據(jù)采集頻率為2 000 Hz,每個試驗點采集20 000組數(shù),每組試驗至少重復(fù)三次。所有試驗均在常壓室溫條件下進行,并待空氣及水流運行穩(wěn)定一段時間后再測量。試驗風(fēng)速范圍為17.8~52.2 m/s,水膜雷諾數(shù)在22~138之間,在此范圍內(nèi)水膜可以均勻地鋪展在鋁板表面。

    圖1 試驗裝置
    Fig.1 Test setup

    3 結(jié)果與分析

    3.1 水膜厚度的瞬態(tài)特性

    試驗狀態(tài)點及典型水膜厚度時域圖分別繪于圖2、圖3和圖4中。盡管水膜厚度時域信號隨風(fēng)速和水流量的變化是一個連續(xù)的過程,仍可參考氣-液兩相分層流動的流域劃分方式[22],將其分為三維波(3D)、滾動波(Roll-wave)以及卷吸(Entrainment)3種不同的流動區(qū)域。三維波區(qū)域內(nèi)的界面波具有振幅小且頻率高的特點,且界面波隨著底部黏性水膜層一起運動,相位速度低。當(dāng)供水流量增大時,黏性底層不足以輸運全部液體,在三維波動之上出現(xiàn)了振幅較大、移動速度更快的滾動波。滾動波由于特殊的運動方式而得名,其瞬態(tài)信號具有相位速度快、隨機性強的特點,水膜表面波的追趕與合并也都發(fā)生在此區(qū)域內(nèi)。卷吸區(qū)發(fā)生在流量大且風(fēng)速高的條件下,波動頂部的液體被氣流卷吸脫離水膜表面。在本文試驗中,風(fēng)速高于36.5 m/s、水膜雷諾數(shù)大于89時能夠觀察氣流中夾帶水滴的現(xiàn)象,然而從水膜厚度時域圖難以區(qū)分滾動波與卷吸區(qū)域。

    本文使用Ishii和Grolmes[24]提出的理論對試驗中發(fā)生的卷吸現(xiàn)象進行簡要分析,圖2展示了室溫條件下本文試驗條件與Ishii和Grolmes的空氣-水兩相流動卷吸判別式的對比。點A實驗條件為Ua=30.5 m/s 且Ref=38,點B實驗條件為Ua=30.5 m/s 且Ref=128。Ishii和Grolmes根據(jù)界面受力分析以及不穩(wěn)定理論建立了卷吸現(xiàn)象的臨界判別式,因此可不受兩相流流道形狀和尺寸的影響。氣液兩相流中的卷吸現(xiàn)象產(chǎn)生的機理按水膜雷諾數(shù)劃分為3種方式:根切(Undercut),滾動波剪切(Roll-wave breakup)以及粗糙湍流(Rough turbulent)。由圖2可知現(xiàn)試驗中的卷吸現(xiàn)象主要為滾動波剪切作用,即當(dāng)氣流剪切力大于水膜表面張力時,滾動波頂部在氣流剪切作用下發(fā)生變形直至斷裂形成水滴,故而主要發(fā)生在高風(fēng)速大流量的試驗點??紤]到水膜對飛機防冰的影響來自存留在表面上的部分,且試驗點大都未達卷吸霧化核心區(qū),故而只關(guān)注表面水膜厚度自身的變化。

    圖2 試驗狀態(tài)點及與文獻[24]卷吸判別式的對比
    Fig.2 Test points and comparison with Ref.[24]entrainment criteria

    圖3 氣流速度U=30.5 m/s,不同供水流量下的水膜厚度時域圖
    Fig.3 Instantaneous thickness of water film with different liquid flow rates under wind speed U=30.5 m/s

    圖4 水膜雷諾數(shù)Re=89,不同風(fēng)速下的水膜厚度時域圖
    Fig.4 Instantaneous thickness of water film with different wind speeds under liquid Reynolds number Re=89

    比較同一風(fēng)速下不同供水流量的水膜厚度時域圖的變化情況(圖3)可以發(fā)現(xiàn),水膜表面滾動波的頻率隨著水膜雷諾數(shù)的增加而逐漸增大,同時基底水膜及滾動波振幅則無明顯變化。相對而言,從圖4可以看出,當(dāng)風(fēng)速發(fā)生變化時,基底水膜厚度及滾動波的頻率和振幅都有較大變化。這表明,高速氣流是水膜表面波動行為的主導(dǎo)因素。由Kelvin-Helmholtz不穩(wěn)定性理論[22,24]可知,當(dāng)上下分層的兩個不同流體之間存在速度差時,流體界面會產(chǎn)生扭曲,兩流體之間發(fā)生位置和動量的交換,即界面波動現(xiàn)象。由于較輕的氣體在上層,較重的液體在下層,當(dāng)速度差較小時界面波動在水膜表面張力和重力的作用下得到抑制,而當(dāng)兩者速度差增大到一定程度后,界面不穩(wěn)定情況加劇,下層液體交換獲得更高的動量,直至克服自身黏性力形成渦旋,即滾動波。若只增大水流量,相間相對速度沒有太大變化時,則界面變形情況較為相似,只是滾動波發(fā)生的頻率提高以承載更多的液體輸運量。

    3.2 平均水膜厚度及預(yù)測模型

    圖5 平均水膜厚度隨風(fēng)速和水膜雷諾數(shù)的變化規(guī)律
    Fig.5 Variation of mean film thickness with different wind speeds and film Reynolds numbers

    另一方面,在風(fēng)速較小的條件下,水膜雷諾數(shù)Ref的變化對平均水膜厚度的影響更大,這是由于氣流對表層水膜的加速作用較弱,上層水膜與底層水膜的速度差小,增加水流量后水膜厚度即迅速增加。當(dāng)風(fēng)速提高后,水膜厚度隨Ref的變化速率隨之降低。而在水膜流動進入卷吸區(qū)時,還可看到水膜厚度變化趨勢發(fā)生輕微的偏轉(zhuǎn),表明卷吸現(xiàn)象進一步削弱了水膜雷諾數(shù)對膜厚的影響。

    fi=0.106Ref-0.133

    (11)

    圖6 界面剪切系數(shù)試驗結(jié)果及擬合曲線
    Fig.6 Test data and correlation for interfacial shearfactor

    值得注意的一點是,本文擬合結(jié)果發(fā)現(xiàn)界面剪切因子隨著水膜雷諾數(shù)的增大而減小,這與水流量的增加引起波動增加從而增強剪切作用的常規(guī)認識相反。隨后可見在分層流Ref=100前后也表現(xiàn)出類似的規(guī)律[19]。這是由于在剪切力計算式中忽略了界面波相位速度的緣故,即水流量增加時,波的相位速度C也會相應(yīng)增大。而通過將風(fēng)速及水流量對波動相位速度的影響分別并入剪切力和界面剪切系數(shù)的計算式中,并不影響最終水膜厚度計算式。

    在關(guān)于飛機表面的水膜流動計算方面,文獻中最終推導(dǎo)得到水平表面水膜厚度計算式的形式大都同式(9)一致[4-5,10-11],而界面剪切力的計算主要應(yīng)用了以下兩種方式:① 近似認為外流場求解得到的壁面摩擦應(yīng)力等于氣-液界面剪切力τ;② 假設(shè)管內(nèi)氣-液兩相流與平板水膜流動規(guī)律相同,使用管內(nèi)流界面應(yīng)力公式計算τ。本文使用這兩種方法得到界面剪切系數(shù)fi,并分別代入水膜計算式中,將預(yù)測結(jié)果與實驗結(jié)果進行了對比,如圖7所示。具體計算方法如下:

    1) 平板邊界層模型[25]

    假設(shè)界面剪切系數(shù)等于半無限大平板湍流邊界層摩擦因子cf,則其計算式為

    fi=cf=0.074(Rea)-0.2

    (12)

    式中:氣流雷諾數(shù)Rea=UaL/νa,其中νa為空氣運動黏度,L為特征長度,此處取值為試驗段高度。

    2) Cheremisinoff分層流模型[19]

    將Ref=100作為小振幅界面波流域及滾動波流域的劃分,可得界面剪切系數(shù)計算式為

    (13)

    從圖7(a)中可以看出,由于水平壁面摩擦系數(shù)顯著小于氣-液波動界面剪切系數(shù),由此計算出來的水膜厚度遠高于試驗值。由于壁面摩擦系數(shù)是在氣流掠過光滑固壁表面上得到的,沒有考慮氣-液界面波動對于空氣流場邊界層的影響,這也符合水膜表面波動能夠增強對氣-液交互作用的結(jié)論。而圖7(b)中表明使用管內(nèi)分層流動經(jīng)驗式計算的界面剪切力同樣偏小,這是由于原試驗中氣流速度較低(<20 m/s),且流道底面的內(nèi)凹形狀令水膜向底部聚集,這樣得到的水膜厚度比平板上的試驗結(jié)果偏大。另外,本文所使用的固壁基底材質(zhì)為親水性的金屬鋁,而文獻中的流道基底材質(zhì)為有機玻璃,表面特性的不同也可能對結(jié)果帶來一定的影響。而使用本文的界面剪切系數(shù)公式計算得到的平均水膜厚度,與試驗結(jié)果相比誤差約在20%以內(nèi)。

    圖7 平均水膜厚度試驗值與各剪切因子模型預(yù)測值的對比
    Fig.7 Comparison of mean film thickness test value with results from different shear factor prediction models

    與經(jīng)典的無量綱水膜厚度經(jīng)驗式的對比,可驗證本文所得的界面剪切力計算方法的正確性。在環(huán)狀氣-液兩相流的研究中,通常認為液體流量,平均膜厚及液膜剪切應(yīng)力三者符合“三角關(guān)系式”,即知道其中任意兩者的值后可確定剩余一項的值。液膜無量綱厚度δ+將平均水膜厚度與剪切應(yīng)力相結(jié)合進行去量綱化,隨后即可表達為無量綱液體流量Ref的函數(shù)。δ+的定義為

    (14)

    (15)

    表1中列出了文獻及本文的無量綱水膜厚度計算模型,各模型的計算值同實驗值的比較結(jié)果如圖8所示。對比發(fā)現(xiàn)使用本文提出的剪切力計算式對平均水膜厚度實驗數(shù)據(jù)進行無量綱化后,與Kosky & Staub根據(jù)水平冷凝環(huán)狀流研究提出的計算模型最為相符。基于逆向豎直環(huán)狀流實驗提出的Hughmark模型在水膜雷諾數(shù)較小的范圍內(nèi)與本試驗結(jié)果符合的較好,但在水流量較大的部分產(chǎn)生明顯差異;而基于豎直向下環(huán)狀流的Asali模型反映出的無量綱水膜厚度隨Ref的變化規(guī)律與試驗結(jié)果相似,但δ+的值整體更小。

    表1無量綱水膜厚度關(guān)系式及應(yīng)用范圍

    Table1Dimensionlessfilmthicknessrelationandtheirapplicationrange

    AuthorApplicationrangeRefRelationKosky&Staub[26]Horizontalcondenseflow13?370δ+=1.316Re0.529fHughmark[27]Upwardflow25?250δ+=0.874Re2/3fAsali[28]Verticalannularflow20?300δ+=0.781Re0.6fPresentworkHorizontalplateflow26?128δ+=1.292Re0.5289f

    圖8 無量綱水膜厚度與各模型預(yù)測結(jié)果的對比
    Fig.8 Comparison of dimensionless water film thickness with results from different prediction models

    4 結(jié) 論

    針對飛機表面結(jié)冰及防冰過程中水膜流動計算和預(yù)測的需求,對高速氣流驅(qū)動的水膜流動進行了試驗研究,分析了水膜厚度瞬態(tài)特征隨氣流速度及水流量的變化規(guī)律,對水膜厚度時域平均值進行統(tǒng)計,校驗了使用平板邊界層模型或管內(nèi)分層流模型的膜厚計算方法,并應(yīng)用新的界面剪切系數(shù)計算式實現(xiàn)平均水膜厚度的準確預(yù)測。研究發(fā)現(xiàn):

    1) 高速氣流驅(qū)動下的水平平板水膜流動,按界面波動特征可分為三維波、滾動波及卷吸區(qū)3種流動區(qū)域,且這3種流域之間存在一個平緩變化的過程。

    2) 平均水膜厚度隨氣流速度的減小或水膜雷諾數(shù)的增大而增大。增大氣流速度時,水膜表面小波動的頻率顯著增加,滾動波的振幅和頻率則有降低;而相同風(fēng)速下,水膜雷諾數(shù)的增大主要表現(xiàn)出滾動波頻率的增大。

    3) 常見的外流場邊界層模型,或低風(fēng)速分層流動中總結(jié)得到的經(jīng)驗關(guān)系式應(yīng)用于氣液兩相界面上出現(xiàn)劇烈的波動時,水膜厚度預(yù)測結(jié)果與實驗數(shù)據(jù)的偏差較大。由此可見,界面剪切應(yīng)力的準確?;谦@得準確水膜厚度的前提和基礎(chǔ)。

    [1] MESSINGER B L. Equilibrium temperature of an un-heated icing surface as a function of air speed[J]. Journal of the Aeronautical Sciences, 1953, 20(1): 29-42.

    [2] AI-KHALIL K M, KEITH T G, DE-WITT K J. Development of an improved model for runback water on aircraft surfaces[J]. Journal of Aircraft, 1994, 31(2): 271-278.

    [3] MYERS T G. Extension to the Messinger model for aircraft icing[J]. AIAA Journal, 2001, 39(2): 211-218.

    [4] MYERS T G, THOMPSON C P. Modeling the flow of water on aircraft in icing conditions[J]. Journal of Aircraft, 1998, 36(6): 1010-1013.

    [5] ALZAILI J, HAMMOND D. Experimental investigation of thin water film stability and its characteristics in SLD icing problem[C]//SAE 2011 International Conference on Aircraft and Engine Icing and Ground Deicing. Chicago: SAE International, 2011.

    [6] DU Y X, GUI Y W, XIAO C H, et al. Investigation on heat transfer characteristics of aircraft icing including runback water[J]. International Journal of Heat and Mass Transfer, 2010, 53(19-20): 3702-3707.

    [7] WRIGHT W B, STRUK P, BARTKUS T,et al. Recent advances in the LEWICE icing model[C]//SAE 2015 International Conference on Icing of Aircraft, Engines, and Structures. Prague: SAE International, 2015.

    [8] HARIRECHE O, VERDIN P, THOMPSON C P, et al. Explicit finite volume modeling of aircraft anti-icing and de-icing[J]. Journal of Aircraft, 2008, 45(6): 1924-1936.

    [9] FORTIN G, LAFORTE J, ILINCA A. Heat and mass transfer during ice accretion on aircraft wings with an improved roughness model[J]. International Journal of Thermal Sciences, 2006, 45(6): 595-606.

    [10] KAREV A R, FARZANEH M, LOZOWSKI E P. Character and stability of a wind-driven supercooled water film on an icing surface—I. Laminar heat transfer[J]. International Journal of Thermal Sciences, 2003, 42(5): 481-498.

    [11] UENO K, FARZANEH M. Linear stability analysis of ice growth under supercooled water film driven by a laminar airflow[J]. Physics of Fluids, 2011, 23(4): 042103.

    [12] WANG G K, ROTHMAYER A P. Thin water films driven by air shear stress through roughness[J]. Computers & Fluids, 2009, 38(2): 235-246.

    [13] WHITE E B, SCHMUCKER J A. A runback criterion for water drops in a turbulent accelerated boundary layer [J]. Journal of Fluids Engineering, 2008, 130(6): 061302.

    [14] 孟繁鑫, 朱光亞, 李榮嘉, 等. 加熱表面水珠運動特性研究[J]. 航空學(xué)報, 2014, 35(5): 1292-1301.

    MENG F X,ZHU G Y, LI R J, et al. Study of water drop motion characteristics on heating surface [J]. Acta Aeronautica et Astronautica Sinica, 2014, 35(5): 1292-1301 (in Chinese).

    [15] MOGHTADERNEJAD S, JADIDI M, NABIL E, et al. Shear driven rivulet dynamics on surfaces with various wettabilities[C]//ASME 2014 International Mechanical Engineering Congress and Exposition. Montreal: American Society of Mechanical Engineers, 2014.

    [16] FEO A, TSAO J. The water film weber number in glaze icing scaling[C]//2007 SAE Aircraft and Engine Icing International Conference. Seville: SAE International, 2007.

    [17] MUZIK T, SAFARIK P, TUCEK A. Analysis of the water film behavior and its breakup on profile using experimental and numerical methods[J]. Journal of Thermal Science, 2014, 23(4): 325-331.

    [18] ZHANG K, WEI T, HU H. An experimental investigation on the surface water transport process over an airfoil by using a digital image projection technique[J]. Experiments in Fluids, 2015, 56(9): 173.

    [19] CHEREMISINOFF N P, DAVIS E J. Stratified turbulent-turbulent gas-liquid flow[J]. AIChE Journal, 1979, 25(1): 48-56.

    [20] ANDRITSOS N, HANRATTY T J. Influence of interfacial waves in stratified gas-liquid flows[J]. AIChE Journal, 1987, 33(3): 444-454.

    [21] PARAS S V, VLACHOS N A, KARABELAS A J. Liquid layer characteristics in stratified—Atomization flow[J]. International Journal of Multiphase Flow, 1994, 20(5): 939-956.

    [22] TZOTZI C, ANDRITSOS N. Interfacial shear stress in wavy stratified gas-liquid flow in horizontal pipes[J]. International Journal of Multiphase Flow, 2013, 54(3): 43-54.

    [23] SETYAWAN A, INDARTO, DEENDARLIANTO. The effect of the fluid properties on the wave velocity and wave frequency of gas-liquid annular two-phase flow in a horizontal pipe[J]. Experimental Thermal and Fluid Science, 2016, 71(4): 25-41.

    [24] ISHII M, GROLMES M A. Inception criteria for droplet entrainment in two-phase concurrent film flow[J]. AIChE Journal, 1975, 21(2): 308-318.

    [25] 吳望一. 流體力學(xué)(下冊)[M]. 北京: 北京大學(xué)出版社, 2004: 370-386.

    WU W Y. Fluiddynamic (Part 2)[M]. Beijing: Peking University Press, 2004: 370-386 (in Chinese).

    [26] KOSKY P G, STAUB F W. Local condensing heat transfer coefficients in the annular flow regime[J]. AIChE Journal, 1971, 17(5): 1037-1043.

    [27] HUGHMARK G A. Film thickness, entrainment, and pressure drop in upward annular and dispersed flow[J]. AIChE Journal, 1973, 19(5): 1062-1065.

    [28] ASALI J C, HANRATTY T T, ANDREUSSI P. Interfacial drag and film height for vertical annular flow[J]. AIChE Journal, 1985, 31(6): 895-902.

    (責(zé)任編輯: 李明敏)

    URL:www.cnki.net/kcms/detail/11.1929.V.20161027.1057.002.html

    Thicknessofwaterfilmdrivenbygasstreamonhorizontalplane

    LENGMengyao,CHANGShinan*,DINGLiang,LIXiaofeng

    SchoolofAeronauticScienceandEngineering,BeihangUniversity,Beijing100083,China

    Liquidwateronthesurfaceofaircraftwillrunbackundertheeffectoftheairflow,resultinginredistributionoficeaccretionandanti-icingheatflux.Experimentalmeasurementandmodelinganalysisareconductedtoinvestigatetheflowbehaviorofshear-drivenwaterfilmonthehorizontalflatsubstrate.Thewaterflowfilmisdriveninawindtunnel,andtheinstantaneousthicknessismeasuredinthesamelocationusingalaserfocusdisplacementmeterbasedonconfocalchromatictechnique.Itisfoundthattheinterfacebetweenthegasandliquidphasesconsistsofunderlyingthinfilmandmultiplescalefluctuations.ThevariationrelationshipofthefilmthicknessbetweenthewindspeedandfilmReynoldsnumberisalsoobtained.Resultsshowthattheaveragefilmthicknessdependsmonotonicallyonthesetwofactors.Basedonfilmflowmodelandexperimentaldata,anewcorrelationforcalculatingtheairshearstressaboveathinfilmisproposedandvalidatedbycomparisonwithpreviousstudies.Thecorrelationcanbeappliedforwaterfilmthicknesscalculationoverarangeofwindspeed(17.8-52.2m/s)andwaterfilmReynoldsnumber(26-128).

    aircraftanti-icing;two-phaseflow;waterfilmthickness;interfacialshearstress;modeling

    2016-08-23;Revised2016-09-18;Accepted2016-10-25;Publishedonline2016-10-271057

    NationalBasicResearchProgramofChina(2015CB755803)

    .E-mailsn_chang@buaa.edu.cn

    2016-08-23;退修日期2016-09-18;錄用日期2016-10-25; < class="emphasis_bold">網(wǎng)絡(luò)出版時間

    時間:2016-10-271057

    www.cnki.net/kcms/detail/11.1929.V.20161027.1057.002.html

    國家“973”計劃 (2015CB755803)

    .E-mailsn_chang@buaa.edu.cn

    冷夢堯, 常士楠, 丁亮, 等. 水平表面氣流剪切作用下的水膜厚度J. 航空學(xué)報,2017,38(2):520696.LENGMY,CHANGSN,DINGL,etal.ThicknessofwaterfilmdrivenbygasstreamonhorizontalplaneJ.ActaAeronauticaetAstronauticaSinica,2017,38(2):520696.

    http://hkxb.buaa.edu.cnhkxb@buaa.edu.cn

    10.7527/S1000-6893.2016.0275

    V211; O359.1

    A

    1000-6893(2017)02-520696-09

    猜你喜歡
    水膜雷諾數(shù)氣流
    巧測水膜張力
    少兒科技(2022年4期)2022-04-14 23:48:10
    氣流的威力
    基于Transition SST模型的高雷諾數(shù)圓柱繞流數(shù)值研究
    濕滑跑道飛機著陸輪胎-水膜-道面相互作用
    固體運載火箭變軌發(fā)動機噴管氣流分離研究
    飛片下的空氣形成的“超強高速氣流刀”
    失穩(wěn)初期的低雷諾數(shù)圓柱繞流POD-Galerkin 建模方法研究
    基于轉(zhuǎn)捩模型的低雷諾數(shù)翼型優(yōu)化設(shè)計研究
    非能動核電站安全殼外壁下降水膜的穩(wěn)定性分析
    民機高速風(fēng)洞試驗的阻力雷諾數(shù)效應(yīng)修正
    av一本久久久久| 精品亚洲成a人片在线观看 | 国产v大片淫在线免费观看| 18禁在线播放成人免费| 最黄视频免费看| 波野结衣二区三区在线| 国产精品麻豆人妻色哟哟久久| 国产无遮挡羞羞视频在线观看| 精品久久久精品久久久| 免费久久久久久久精品成人欧美视频 | 日本黄大片高清| 亚洲久久久国产精品| 天堂中文最新版在线下载| 99久久精品热视频| 久久午夜福利片| 特大巨黑吊av在线直播| 亚洲欧美日韩东京热| 免费高清在线观看视频在线观看| 免费大片18禁| 妹子高潮喷水视频| 久久亚洲国产成人精品v| 男女啪啪激烈高潮av片| 免费久久久久久久精品成人欧美视频 | av黄色大香蕉| 亚洲经典国产精华液单| 高清日韩中文字幕在线| 丰满少妇做爰视频| 国产成人a区在线观看| 最近中文字幕2019免费版| 日韩制服骚丝袜av| 春色校园在线视频观看| 精品一区二区三卡| 亚洲精品视频女| 日韩中文字幕视频在线看片 | 91午夜精品亚洲一区二区三区| 色吧在线观看| 麻豆国产97在线/欧美| 国产男人的电影天堂91| 午夜福利在线在线| 亚洲av国产av综合av卡| 午夜福利网站1000一区二区三区| 我的老师免费观看完整版| 欧美精品人与动牲交sv欧美| 日韩欧美 国产精品| 老司机影院毛片| 久久久欧美国产精品| 免费看日本二区| 嫩草影院入口| 国产一区亚洲一区在线观看| 深夜a级毛片| 亚洲欧美精品自产自拍| 国产精品一区二区三区四区免费观看| h日本视频在线播放| 亚洲av日韩在线播放| 欧美亚洲 丝袜 人妻 在线| 免费播放大片免费观看视频在线观看| 各种免费的搞黄视频| 国产高清国产精品国产三级 | 久热这里只有精品99| 日韩国内少妇激情av| 人体艺术视频欧美日本| 高清黄色对白视频在线免费看 | 亚洲av中文av极速乱| 免费黄网站久久成人精品| 国产 精品1| 欧美成人午夜免费资源| 一级黄片播放器| 国产亚洲av片在线观看秒播厂| 看免费成人av毛片| 日韩av在线免费看完整版不卡| 少妇被粗大猛烈的视频| 日韩强制内射视频| 欧美精品亚洲一区二区| 99久久人妻综合| 国产在线视频一区二区| 亚洲精品视频女| 久久久国产一区二区| 亚洲综合精品二区| 欧美高清性xxxxhd video| 99久久精品一区二区三区| 在线天堂最新版资源| 免费播放大片免费观看视频在线观看| 色哟哟·www| 夜夜骑夜夜射夜夜干| 国产在线一区二区三区精| 日日摸夜夜添夜夜添av毛片| 寂寞人妻少妇视频99o| 欧美日韩视频精品一区| 国产精品一区二区在线不卡| 丝袜脚勾引网站| 亚洲精品一区蜜桃| 久久久久久久久久成人| 在线精品无人区一区二区三 | 一级爰片在线观看| 亚洲精品一区蜜桃| 国产精品麻豆人妻色哟哟久久| 老师上课跳d突然被开到最大视频| 一级av片app| 麻豆精品久久久久久蜜桃| 狂野欧美白嫩少妇大欣赏| 熟女人妻精品中文字幕| 99久久精品国产国产毛片| 精品久久久久久久末码| 国产高清有码在线观看视频| 亚洲av电影在线观看一区二区三区| 嫩草影院新地址| 成人18禁高潮啪啪吃奶动态图 | 久久久国产一区二区| 少妇人妻久久综合中文| 七月丁香在线播放| 久久久成人免费电影| 久久国产乱子免费精品| 日本vs欧美在线观看视频 | 男的添女的下面高潮视频| 亚洲精品亚洲一区二区| 亚洲av中文av极速乱| 久久久久久伊人网av| 最后的刺客免费高清国语| 精品亚洲成国产av| 国产高潮美女av| 丰满乱子伦码专区| 国产欧美日韩一区二区三区在线 | 人人妻人人爽人人添夜夜欢视频 | 大码成人一级视频| 欧美3d第一页| 色5月婷婷丁香| 制服丝袜香蕉在线| 久久99蜜桃精品久久| 亚洲精品久久久久久婷婷小说| 中文天堂在线官网| 国产在线免费精品| 亚洲欧美日韩东京热| 欧美人与善性xxx| 国产毛片在线视频| 赤兔流量卡办理| 成人美女网站在线观看视频| 丝袜喷水一区| 久久精品熟女亚洲av麻豆精品| 国产成人一区二区在线| 免费观看av网站的网址| 欧美精品人与动牲交sv欧美| 久久国产乱子免费精品| 大码成人一级视频| 18禁在线播放成人免费| 在现免费观看毛片| 亚洲高清免费不卡视频| 欧美性感艳星| 久久人人爽人人爽人人片va| 亚洲国产av新网站| 国产一区二区三区av在线| 国产免费一级a男人的天堂| 亚洲在久久综合| 青春草视频在线免费观看| 欧美成人午夜免费资源| 老熟女久久久| 久久精品久久久久久久性| 免费黄频网站在线观看国产| 色网站视频免费| 精品国产乱码久久久久久小说| 麻豆国产97在线/欧美| 最近手机中文字幕大全| 亚洲国产精品成人久久小说| 精品久久久久久久末码| 国产成人a∨麻豆精品| 国产免费一区二区三区四区乱码| 大香蕉97超碰在线| 国产 精品1| 热re99久久精品国产66热6| 六月丁香七月| 亚洲精品第二区| 国产91av在线免费观看| 日韩欧美一区视频在线观看 | 性高湖久久久久久久久免费观看| 成人特级av手机在线观看| 欧美成人精品欧美一级黄| 人体艺术视频欧美日本| 视频中文字幕在线观看| 涩涩av久久男人的天堂| 亚洲天堂av无毛| 免费观看无遮挡的男女| 亚洲人成网站在线观看播放| 日韩中字成人| 在线观看美女被高潮喷水网站| 最近2019中文字幕mv第一页| 高清毛片免费看| 欧美97在线视频| 在线播放无遮挡| 免费av不卡在线播放| 久久99蜜桃精品久久| 777米奇影视久久| 青春草亚洲视频在线观看| 男男h啪啪无遮挡| 在线免费十八禁| 国产精品国产三级国产专区5o| 成年免费大片在线观看| 国产视频内射| 日韩一区二区三区影片| 嘟嘟电影网在线观看| 建设人人有责人人尽责人人享有的 | 大片免费播放器 马上看| 久久99热6这里只有精品| av在线app专区| 国产精品不卡视频一区二区| 久久久久久久久久人人人人人人| 国产极品天堂在线| 十分钟在线观看高清视频www | 最近中文字幕2019免费版| 精品国产露脸久久av麻豆| 夜夜看夜夜爽夜夜摸| 亚洲av中文av极速乱| 国产一区二区三区av在线| 欧美日韩精品成人综合77777| 爱豆传媒免费全集在线观看| 国产精品三级大全| 久久精品久久精品一区二区三区| 久久精品久久久久久久性| 免费看不卡的av| freevideosex欧美| 多毛熟女@视频| 五月天丁香电影| 国产色婷婷99| 身体一侧抽搐| 亚洲国产欧美在线一区| 亚洲精品国产av成人精品| 国产色爽女视频免费观看| 成人亚洲欧美一区二区av| 26uuu在线亚洲综合色| 亚洲高清免费不卡视频| 国产 一区 欧美 日韩| 欧美最新免费一区二区三区| 在线观看一区二区三区| 久久国产亚洲av麻豆专区| 熟妇人妻不卡中文字幕| 人人妻人人澡人人爽人人夜夜| 亚洲国产精品专区欧美| 日韩免费高清中文字幕av| 97超碰精品成人国产| 最近最新中文字幕免费大全7| 国产69精品久久久久777片| 全区人妻精品视频| 欧美zozozo另类| 少妇的逼水好多| 亚洲伊人久久精品综合| 91精品国产九色| 国内少妇人妻偷人精品xxx网站| 亚洲高清免费不卡视频| 伦精品一区二区三区| 国产亚洲91精品色在线| 九九爱精品视频在线观看| 男人舔奶头视频| 我要看黄色一级片免费的| 高清视频免费观看一区二区| 伦理电影免费视频| 少妇精品久久久久久久| 少妇人妻 视频| 97精品久久久久久久久久精品| 成人高潮视频无遮挡免费网站| 2018国产大陆天天弄谢| 亚洲人成网站在线播| 国产乱人视频| 亚洲欧美中文字幕日韩二区| 天堂中文最新版在线下载| 亚洲aⅴ乱码一区二区在线播放| 亚洲人成网站在线观看播放| 免费久久久久久久精品成人欧美视频 | 成人一区二区视频在线观看| 亚洲国产精品专区欧美| 99热这里只有精品一区| 国产 一区精品| 国产人妻一区二区三区在| 色网站视频免费| 久久99热这里只频精品6学生| 精品视频人人做人人爽| 亚洲国产欧美在线一区| 亚洲不卡免费看| 久久99精品国语久久久| 99热全是精品| 午夜福利高清视频| 国产淫片久久久久久久久| 夫妻性生交免费视频一级片| 久久97久久精品| 国产又色又爽无遮挡免| 麻豆国产97在线/欧美| 精品亚洲成国产av| 亚洲av男天堂| 老熟女久久久| 亚洲怡红院男人天堂| 性色av一级| 午夜福利网站1000一区二区三区| 热99国产精品久久久久久7| h日本视频在线播放| 少妇裸体淫交视频免费看高清| 久久久精品免费免费高清| 亚州av有码| 亚洲av在线观看美女高潮| 国产精品无大码| 少妇 在线观看| 高清黄色对白视频在线免费看 | 一个人免费看片子| 七月丁香在线播放| 新久久久久国产一级毛片| 精品国产乱码久久久久久小说| 99热这里只有是精品在线观看| 少妇熟女欧美另类| 永久免费av网站大全| 春色校园在线视频观看| 免费av不卡在线播放| av在线播放精品| 尾随美女入室| 舔av片在线| 国产伦精品一区二区三区四那| 精品午夜福利在线看| av线在线观看网站| 亚洲国产精品专区欧美| 精品人妻一区二区三区麻豆| 亚洲av.av天堂| 熟妇人妻不卡中文字幕| 搡女人真爽免费视频火全软件| 嘟嘟电影网在线观看| av在线观看视频网站免费| 日韩 亚洲 欧美在线| 男女免费视频国产| 男人爽女人下面视频在线观看| 国产男女超爽视频在线观看| 亚洲国产最新在线播放| 熟女电影av网| 大香蕉97超碰在线| 高清黄色对白视频在线免费看 | 久久国产精品男人的天堂亚洲 | 国产欧美日韩精品一区二区| 精品久久国产蜜桃| av视频免费观看在线观看| 欧美精品国产亚洲| 在线观看免费视频网站a站| 亚洲av男天堂| 欧美丝袜亚洲另类| 国产成人一区二区在线| 麻豆成人av视频| 18禁裸乳无遮挡动漫免费视频| 青春草国产在线视频| 色网站视频免费| 日日啪夜夜爽| 一区二区三区精品91| 一区二区三区乱码不卡18| 免费看不卡的av| 亚洲美女黄色视频免费看| 午夜福利影视在线免费观看| 欧美老熟妇乱子伦牲交| 亚洲av成人精品一二三区| 国产人妻一区二区三区在| 看免费成人av毛片| 777米奇影视久久| 国产v大片淫在线免费观看| 亚洲第一区二区三区不卡| 免费大片18禁| 涩涩av久久男人的天堂| 只有这里有精品99| 2021少妇久久久久久久久久久| 在线播放无遮挡| av免费在线看不卡| 久久久久性生活片| 99视频精品全部免费 在线| 在线观看国产h片| 一个人免费看片子| 日韩中字成人| 日韩一区二区三区影片| 欧美日韩在线观看h| 性色av一级| 欧美精品一区二区免费开放| 欧美 日韩 精品 国产| 哪个播放器可以免费观看大片| 亚洲精品456在线播放app| 国产真实伦视频高清在线观看| a级毛色黄片| 国产黄频视频在线观看| 久久久久久久大尺度免费视频| 少妇人妻一区二区三区视频| 国产爽快片一区二区三区| 搡女人真爽免费视频火全软件| 国产久久久一区二区三区| 亚洲伊人久久精品综合| 国产成人a区在线观看| xxx大片免费视频| 大又大粗又爽又黄少妇毛片口| 大码成人一级视频| 80岁老熟妇乱子伦牲交| 18禁在线播放成人免费| 岛国毛片在线播放| 精品久久久精品久久久| 国产女主播在线喷水免费视频网站| 亚洲四区av| 成人毛片60女人毛片免费| 国产淫片久久久久久久久| 22中文网久久字幕| 国产精品国产三级国产专区5o| 自拍偷自拍亚洲精品老妇| 免费看日本二区| 少妇丰满av| 亚洲国产精品国产精品| 在线观看三级黄色| 久久久欧美国产精品| 一级二级三级毛片免费看| 一个人看视频在线观看www免费| 国精品久久久久久国模美| 国产成人精品福利久久| 午夜福利影视在线免费观看| 欧美精品一区二区大全| 国产亚洲一区二区精品| 精华霜和精华液先用哪个| 九九爱精品视频在线观看| 七月丁香在线播放| av专区在线播放| 成人特级av手机在线观看| 亚洲无线观看免费| 91精品伊人久久大香线蕉| 亚洲成人av在线免费| 欧美成人一区二区免费高清观看| 我的老师免费观看完整版| 最近的中文字幕免费完整| 午夜福利在线在线| 久久久久久久国产电影| 国产精品国产三级国产av玫瑰| 欧美97在线视频| av一本久久久久| 久久久久久久久久成人| 好男人视频免费观看在线| 日本-黄色视频高清免费观看| 在线免费十八禁| 2022亚洲国产成人精品| 欧美极品一区二区三区四区| 国产精品国产三级国产av玫瑰| 国产淫语在线视频| 青青草视频在线视频观看| 亚洲中文av在线| 精品久久久精品久久久| 精品人妻偷拍中文字幕| 蜜桃亚洲精品一区二区三区| 七月丁香在线播放| 国产av一区二区精品久久 | 极品教师在线视频| 一级毛片 在线播放| 国产成人精品一,二区| 97在线视频观看| 高清视频免费观看一区二区| 精品久久久久久久末码| 国产高潮美女av| 天堂俺去俺来也www色官网| 日韩中字成人| 国产一区有黄有色的免费视频| 一级a做视频免费观看| 国产毛片在线视频| 久久午夜福利片| 午夜福利网站1000一区二区三区| 久久精品久久久久久久性| 成人国产麻豆网| 国产成人a区在线观看| a级毛片免费高清观看在线播放| 在线精品无人区一区二区三 | 中文字幕av成人在线电影| 免费看av在线观看网站| 美女国产视频在线观看| 80岁老熟妇乱子伦牲交| 一区二区三区四区激情视频| 精品国产露脸久久av麻豆| 人妻少妇偷人精品九色| 国产黄频视频在线观看| 又爽又黄a免费视频| 自拍欧美九色日韩亚洲蝌蚪91 | 日本av免费视频播放| 国产一区有黄有色的免费视频| 久久99蜜桃精品久久| 九九爱精品视频在线观看| av又黄又爽大尺度在线免费看| 午夜免费男女啪啪视频观看| 久久 成人 亚洲| 天堂中文最新版在线下载| 免费人妻精品一区二区三区视频| 亚洲美女搞黄在线观看| 高清av免费在线| 欧美老熟妇乱子伦牲交| 18禁在线无遮挡免费观看视频| 国产精品久久久久久久电影| 少妇熟女欧美另类| 国产精品av视频在线免费观看| 精品一区二区三卡| 人人妻人人添人人爽欧美一区卜 | 久久久久久久精品精品| 夜夜看夜夜爽夜夜摸| 一级黄片播放器| 视频区图区小说| 国产色婷婷99| 久久久久久久精品精品| 亚洲第一区二区三区不卡| 麻豆精品久久久久久蜜桃| 亚洲国产av新网站| 国产av一区二区精品久久 | 日日摸夜夜添夜夜添av毛片| 精品久久久久久久末码| 男女啪啪激烈高潮av片| 国产精品久久久久久精品古装| 国产成人免费观看mmmm| 免费看不卡的av| 天天躁日日操中文字幕| 边亲边吃奶的免费视频| 国产69精品久久久久777片| 日本vs欧美在线观看视频 | 亚洲真实伦在线观看| av线在线观看网站| av一本久久久久| 三级经典国产精品| av女优亚洲男人天堂| 欧美精品一区二区免费开放| 国产高清三级在线| 久久久久人妻精品一区果冻| 大香蕉久久网| 在现免费观看毛片| av免费观看日本| 黄色欧美视频在线观看| 人人妻人人添人人爽欧美一区卜 | 国产成人a∨麻豆精品| 男的添女的下面高潮视频| 中文字幕亚洲精品专区| 国产男女内射视频| 成人免费观看视频高清| 欧美精品亚洲一区二区| 一区二区av电影网| 国产午夜精品一二区理论片| videossex国产| 2022亚洲国产成人精品| 中文在线观看免费www的网站| 我要看日韩黄色一级片| 在线天堂最新版资源| 麻豆乱淫一区二区| 十分钟在线观看高清视频www | 欧美日韩综合久久久久久| 秋霞在线观看毛片| 国产精品一及| 18+在线观看网站| 99久久综合免费| 婷婷色麻豆天堂久久| 尤物成人国产欧美一区二区三区| 中国国产av一级| 一级毛片黄色毛片免费观看视频| 妹子高潮喷水视频| 国产男人的电影天堂91| 黄色一级大片看看| 狂野欧美激情性bbbbbb| 少妇人妻久久综合中文| 中文资源天堂在线| 久久国内精品自在自线图片| 亚洲成人中文字幕在线播放| 黄色欧美视频在线观看| 久久国产亚洲av麻豆专区| 日韩av不卡免费在线播放| 亚洲av在线观看美女高潮| 激情 狠狠 欧美| 亚洲国产高清在线一区二区三| 国产精品一二三区在线看| 99久久中文字幕三级久久日本| 高清视频免费观看一区二区| 久久久色成人| 美女主播在线视频| 亚洲av中文字字幕乱码综合| 男男h啪啪无遮挡| 午夜福利在线在线| 国产精品麻豆人妻色哟哟久久| 久久久久久久精品精品| 交换朋友夫妻互换小说| 久久精品国产a三级三级三级| 国产视频首页在线观看| 最近中文字幕2019免费版| 干丝袜人妻中文字幕| 777米奇影视久久| 国产精品.久久久| 久久人人爽人人片av| 97超视频在线观看视频| 黄色配什么色好看| 免费黄色在线免费观看| 一本—道久久a久久精品蜜桃钙片| 久久精品国产亚洲av涩爱| 毛片女人毛片| 又大又黄又爽视频免费| 欧美高清成人免费视频www| 美女cb高潮喷水在线观看| 国产真实伦视频高清在线观看| 男女无遮挡免费网站观看| 六月丁香七月| 免费人妻精品一区二区三区视频| 欧美一级a爱片免费观看看| 美女福利国产在线 | 精品人妻熟女av久视频| 久久久久久伊人网av| 伊人久久精品亚洲午夜| 久久久久国产精品人妻一区二区| 中国三级夫妇交换| 欧美激情极品国产一区二区三区 | 国产成人91sexporn| 内射极品少妇av片p| 男人爽女人下面视频在线观看| 亚洲国产av新网站| 欧美+日韩+精品| 水蜜桃什么品种好| 最近最新中文字幕免费大全7| 免费大片黄手机在线观看| 亚洲欧美日韩另类电影网站 | 尤物成人国产欧美一区二区三区| 少妇裸体淫交视频免费看高清| 国内少妇人妻偷人精品xxx网站| 国产大屁股一区二区在线视频| 亚洲色图综合在线观看| 日日啪夜夜撸| 亚洲av成人精品一区久久| av卡一久久| 久久国产乱子免费精品| 久久久欧美国产精品|