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

    螺旋管結(jié)構(gòu)對(duì)超臨界甲烷換熱特性的影響分析

    2023-10-14 13:41:28李選平鄭文科田中允姜益強(qiáng)
    煤氣與熱力 2023年9期
    關(guān)鍵詞:升角螺旋管壓力降

    李選平, 鄭文科, 田中允, 姜益強(qiáng), 孫 澄

    (1.哈爾濱工業(yè)大學(xué)建筑學(xué)院,黑龍江哈爾濱150001;2.寒地城鄉(xiāng)人居環(huán)境科學(xué)與技術(shù)工業(yè)和信息化部重點(diǎn)實(shí)驗(yàn)室,黑龍江哈爾濱150090)

    1 概述

    我國(guó)海上天然氣資源十分豐富,為了實(shí)現(xiàn)長(zhǎng)久穩(wěn)定的天然氣供給,需要對(duì)儲(chǔ)量豐富的海上天然氣進(jìn)行開發(fā)利用[1]。為了方便運(yùn)輸,在海上天然氣的開采過程中,需要對(duì)天然氣進(jìn)行液化工作[2]。天然氣的液化工藝包括階式制冷循環(huán)工藝、混合制冷劑制冷液化工藝以及膨脹機(jī)制冷循環(huán)液化工藝[3]。其中,繞管式換熱器為混合制冷劑制冷液化工藝的主要設(shè)備,一般由3股流體通道構(gòu)成,分別為天然氣、液態(tài)混合制冷劑及氣態(tài)混合制冷劑,其中天然氣一般處于超臨界狀態(tài)[4]。我國(guó)的大型繞管式換熱器主要依靠進(jìn)口,因此為了打破國(guó)外技術(shù)壟斷,對(duì)超臨界天然氣在螺旋管內(nèi)的流動(dòng)換熱特性進(jìn)行研究具有重大意義。

    國(guó)內(nèi)外學(xué)者對(duì)超臨界流體流動(dòng)傳熱特性進(jìn)行了相關(guān)研究。Bai等[5]基于PCHE數(shù)值模型,探究了熱流密度和壓力對(duì)正弦波紋通道內(nèi)超臨界液化天然氣熱工水力特性的影響。Wang等[6]對(duì)超臨界水在傾斜角度為25°的光滑管內(nèi)的傳熱進(jìn)行了實(shí)驗(yàn)和數(shù)值研究。實(shí)驗(yàn)結(jié)果表明,質(zhì)量流量的增加有利于傳熱。隨著壓力增加,物理性質(zhì)變化對(duì)傳熱的影響變得平緩。Wang等[7]采用直接模擬的方法,研究了在混合對(duì)流和強(qiáng)制對(duì)流條件下,超臨界CO2在半圓形管道中的傳熱過程,并探究了熱物性變化對(duì)傳熱性能的影響。李仲珍等[8]利用數(shù)值模擬的方法對(duì)超臨界液化天然氣在不同結(jié)構(gòu)的豎直管內(nèi)的強(qiáng)制對(duì)流換熱進(jìn)行相關(guān)研究,并驗(yàn)證了甲烷代替LNG進(jìn)行數(shù)值模擬的準(zhǔn)確性。蔡浩飛[9]利用數(shù)值模擬及實(shí)驗(yàn)研究相結(jié)合的方式,探究了超臨界CO2在毛細(xì)管換熱器內(nèi)的冷卻傳熱機(jī)理,并針對(duì)變物性換熱器提出優(yōu)化準(zhǔn)則。楊梅[10]利用數(shù)值模擬的方法,對(duì)超臨界CO2在不同形式的螺旋管及直管內(nèi)對(duì)流換熱性能展開研究。于博文等[11]采用數(shù)值模擬的方法對(duì)超臨界CO2/丙烷混合物在水平直管內(nèi)的流動(dòng)換熱特性展開研究,重點(diǎn)探究了工質(zhì)配比濃度、運(yùn)行壓力、質(zhì)量流量及熱流密度對(duì)換熱性能的影響。關(guān)于超臨界甲烷的研究,張羽楠[12]對(duì)超臨界甲烷在水平直管及U形管中的流動(dòng)換熱性能展開研究,利用數(shù)值模擬重點(diǎn)分析了質(zhì)量流量、熱流密度、壓力及重力對(duì)流動(dòng)換熱性能的影響規(guī)律。孫鵬等[13]運(yùn)用相同的研究方法研究了彎管曲率半徑以及U形管道不同布置形式對(duì)超臨界甲烷在U形管內(nèi)的流動(dòng)換熱特性的影響。周昭勇等[14]針對(duì)超臨界甲烷在水平圓管中的冷卻傳熱的影響因素進(jìn)行相關(guān)研究,重點(diǎn)探索了質(zhì)量流量、熱流密度及不同入口壓力對(duì)傳熱系數(shù)的影響。

    此外,關(guān)于流體在螺旋管內(nèi)流動(dòng)傳熱特性方面,國(guó)內(nèi)外也有諸多學(xué)者對(duì)此展開相關(guān)研究。Zhai等[15]對(duì)螺旋管內(nèi)納米流體的流動(dòng)和傳熱特性進(jìn)行了實(shí)驗(yàn)研究,探究了螺距、旋轉(zhuǎn)角度和納米粒子質(zhì)量分?jǐn)?shù)對(duì)流動(dòng)和傳熱性能的影響。Dong等[16]對(duì)高黏度熱油和水在螺旋扭曲管換熱器中的傳熱和流阻特性進(jìn)行了實(shí)驗(yàn)研究,并將螺旋扭曲管換熱器與普通管和螺旋槽管管殼式換熱器的傳熱和流阻性能進(jìn)行比較。Yu等[17]分析了螺旋翅片管內(nèi)低雷諾數(shù)流動(dòng)的速度和溫度分布的解析解,并提出了最佳翅片數(shù)量和最佳扭轉(zhuǎn)比之間的關(guān)系。刁林[18]對(duì)超臨界狀態(tài)下的CO2在直管及螺旋管內(nèi)的流動(dòng)換熱特性進(jìn)行研究,主要探索了不同通道傾角、熱流密度、重力水平、操作壓力等對(duì)平均表面?zhèn)鳠嵯禂?shù)的影響。馮璐璐等[19]基于Realizablek-ε湍流模型及第一類熱邊界條件下的薄壁熱阻模型,采用CFD模擬技術(shù)對(duì)螺旋彎管內(nèi)水的流動(dòng)與換熱特性進(jìn)行數(shù)值模擬,探究了螺旋彎管內(nèi)流動(dòng)換熱特性。趙敏等[20]采用數(shù)值模擬的方法研究了兩圈纏繞管內(nèi)部超臨界甲烷的流動(dòng)換熱特性,重點(diǎn)分析了流體物性變化和離心力對(duì)流動(dòng)換熱特性的影響。崔海亭等[21]利用Fluent軟件對(duì)超臨界CO2在水平螺旋管內(nèi)的流動(dòng)換熱性能進(jìn)行模擬研究,重點(diǎn)研究了結(jié)構(gòu)參數(shù)及有無重力對(duì)流動(dòng)換熱特性的影響。

    綜上,目前對(duì)以上兩方面的研究方法較為完善,但對(duì)超臨界甲烷在螺旋管內(nèi)流動(dòng)換熱特性的研究仍需探索。本文將超臨界甲烷作為螺旋管內(nèi)的工質(zhì),利用數(shù)值模擬的方法探究不同結(jié)構(gòu)參數(shù)下超臨界甲烷在螺旋管內(nèi)的流動(dòng)換熱特性。主要探究螺旋管內(nèi)直徑、螺旋直徑及螺旋升角對(duì)超臨界甲烷在螺旋管內(nèi)流動(dòng)換熱性能的影響規(guī)律,為我國(guó)大型繞管式換熱器的設(shè)計(jì)研究提供參考。

    2 計(jì)算模型建立

    2.1 超臨界甲烷物性分析

    本文選取甲烷為工作流體,其臨界溫度和臨界壓力分別為190.56 K和4.599 2 MPa。通過NIST軟件查詢壓力為7 MPa下的甲烷物性參數(shù),甲烷物性參數(shù)隨溫度變化見圖1。不同壓力下的超臨界甲烷,其物性均會(huì)在某溫度點(diǎn)產(chǎn)生突變,稱此溫度點(diǎn)為擬臨界溫度。由圖1可知,在壓力恒定的情況下,比定壓熱容隨著溫度升高呈現(xiàn)出先增加后降低再緩慢上升的變化趨勢(shì),且在擬臨界溫度附近達(dá)到峰值。密度隨著溫度升高不斷減小,并在擬臨界溫度附近陡降。熱導(dǎo)率和黏度均隨著溫度升高呈現(xiàn)先降低后增加的變化趨勢(shì),且均在擬臨界溫度附近達(dá)到最低值。

    2.2 物理模型

    繞管式換熱器內(nèi)由多根螺旋管組成,纏繞結(jié)構(gòu)復(fù)雜且實(shí)際尺寸較大,綜合計(jì)算成本和計(jì)算精度,對(duì)大型繞管式換熱器管側(cè)流動(dòng)換熱特性的研究可以簡(jiǎn)化為對(duì)半圈螺旋管內(nèi)流動(dòng)換熱特性的研究。為了探究螺旋管結(jié)構(gòu)參數(shù)對(duì)超臨界甲烷在螺旋管內(nèi)流動(dòng)換熱特性的影響,本文構(gòu)建了不同結(jié)構(gòu)的模型進(jìn)行模擬計(jì)算,物理模型結(jié)構(gòu)參數(shù)見表1。

    為了確保模擬的準(zhǔn)確性,將模型分為發(fā)展段和測(cè)量段兩個(gè)部分,發(fā)展段長(zhǎng)為800 mm,保證流體流動(dòng)穩(wěn)定且接近真實(shí)流動(dòng)工況。物理模型幾何結(jié)構(gòu)見圖2,其中R為螺旋半徑(是螺旋直徑的1/2)。螺旋升角記為β,指螺旋管管壁上任何一點(diǎn)切線與螺旋管纏繞軸線的垂直平面之間的夾角。

    2.3 控制方程

    超臨界狀態(tài)下螺旋管內(nèi)流動(dòng)與換熱的模擬是基于連續(xù)性方程、動(dòng)量守恒方程、能量守恒方程以及湍流模型的求解過程,從而獲得螺旋管內(nèi)流動(dòng)中各相關(guān)參數(shù)的變化情況。連續(xù)性方程、動(dòng)量守恒方程、能量守恒方程及湍流模型的數(shù)學(xué)描述見文獻(xiàn)[7]。

    圖1 壓力為7 MPa時(shí)超臨界甲烷物性參數(shù)隨溫度變化

    表1 物理模型結(jié)構(gòu)參數(shù)

    圖2 物理模型幾何結(jié)構(gòu)

    2.4 數(shù)值方法

    本文采用Fluent軟件對(duì)超臨界下甲烷在螺旋管內(nèi)的流動(dòng)換熱特性進(jìn)行模擬。利用Gambit軟件對(duì)計(jì)算流體域進(jìn)行網(wǎng)格劃分,由于近壁面處溫度、速度以及物性變化較為劇烈,為精確捕捉近壁面處流體的流動(dòng)參數(shù),需要對(duì)近壁面處進(jìn)行網(wǎng)格加密處理。

    其中,內(nèi)直徑為10 mm的螺旋管網(wǎng)格尺寸設(shè)置如下:全局網(wǎng)格尺寸為0.6 mm,邊界層網(wǎng)格首層厚度為0.005 mm,邊界層網(wǎng)格有15層,增長(zhǎng)比例設(shè)置為1.2,螺旋管端面網(wǎng)格及邊界層局部放大圖見圖3。螺旋管入口采用質(zhì)量流率入口邊界,質(zhì)量流率設(shè)置為300 kg/(m2·s),入口壓力為7 MPa。入口邊界設(shè)置中的湍流指定方法選擇湍流強(qiáng)度及水力直徑,其中湍流強(qiáng)度為5%,水力直徑為實(shí)際螺旋管內(nèi)直徑。螺旋管出口采用自由邊界出口,壁面均采用定熱流無滑移光滑壁面邊界條件,不考慮壁面厚度,壁面定熱流密度設(shè)置為-20 kW/m2,其中負(fù)號(hào)表示甲烷向外界環(huán)境散熱,為冷卻過程。

    圖3 螺旋管端面網(wǎng)格及邊界層網(wǎng)格局部放大圖

    此模擬研究中,工質(zhì)物性設(shè)置選取變物性模型。壓力-速度耦合選用SIMPLEC算法,能量方程、動(dòng)量方程、湍流強(qiáng)度及湍動(dòng)耗散率選擇較為嚴(yán)格的二階迎風(fēng)格式。為了保證求解過程的穩(wěn)定性,在本次計(jì)算中除了能量的松弛因子設(shè)置為1外,其余松弛因子均設(shè)置為0.1。由于本文數(shù)值模擬為穩(wěn)態(tài)模擬,因此初始值的設(shè)定采用混合初始化的方式。當(dāng)能量、質(zhì)量、速度的殘差數(shù)量級(jí)小于10-6時(shí),認(rèn)定此次模擬計(jì)算收斂。

    2.5 模型驗(yàn)證及數(shù)據(jù)計(jì)算

    模擬得到的數(shù)據(jù)有螺旋管壁面溫度、螺旋管流體溫度、總壓力降,通過計(jì)算得到螺旋管的表面?zhèn)鳠嵯禂?shù)與摩擦壓力降。本文的表面?zhèn)鳠嵯禂?shù)均指測(cè)量段的內(nèi)表面?zhèn)鳠嵯禂?shù),壓力降也均指測(cè)量段的壓力降。

    表面?zhèn)鳠嵯禂?shù)表征超臨界甲烷在螺旋管內(nèi)的換熱特性,計(jì)算式為:

    (1)

    式中h——表面?zhèn)鳠嵯禂?shù),W/(m2·K)

    q——熱流密度,W/m2

    Tw——螺旋管壁面平均溫度,K

    Tb——螺旋管流體的平均溫度,K

    螺旋管內(nèi)的總壓力降可根據(jù)Fluent軟件直接讀取。在本研究中,總壓力降包括摩擦壓力降、加速壓力降以及重力壓力降。其中,加速壓力降產(chǎn)生的主要原因是超臨界流體的密度隨溫度和壓力而變化,其計(jì)算式為:

    (2)

    式中 Δpac——加速壓力降,Pa

    ρout——螺旋管測(cè)量段出口處平均密度,kg/m3

    vout——螺旋管測(cè)量段出口處平均流速,m/s

    ρin——螺旋管測(cè)量段入口處平均密度,kg/m3

    vin——螺旋管測(cè)量段入口處平均流速,m/s

    重力壓力降產(chǎn)生的主要原因是流道內(nèi)不同高度處流體位能不同而引起靜壓的變化,其計(jì)算式為:

    (3)

    式中 Δpel——重力壓力降,Pa

    g——重力加速度,m/s2

    hout——螺旋管測(cè)量段出口處的水平高度,m

    hin——螺旋管測(cè)量段入口處的水平高度,m

    摩擦壓力降產(chǎn)生的主要原因是管道壁面摩擦產(chǎn)生壓頭損失。摩擦壓力降計(jì)算式為:

    Δpf=Δp-Δpac-Δpel

    (4)

    式中 Δpf——摩擦壓力降,Pa

    Δp——螺旋管內(nèi)總壓力降,Pa

    為了驗(yàn)證所建模型與研究方法的準(zhǔn)確性及適用性,本文在同等條件下將數(shù)值模擬結(jié)果與文獻(xiàn)[22]中壓力為9 MPa、質(zhì)量流率為239 kg/(m2·s)、熱流密度為13.3 kW/m2的工況下的實(shí)驗(yàn)數(shù)據(jù)進(jìn)行對(duì)比,對(duì)比結(jié)果見圖4。圖4中流體溫度為測(cè)量段流體平均溫度,實(shí)驗(yàn)數(shù)據(jù)與模擬數(shù)據(jù)的平均相對(duì)誤差在5%以內(nèi),驗(yàn)證了模型的適用性。

    圖4 實(shí)驗(yàn)數(shù)據(jù)與模擬數(shù)據(jù)的對(duì)比

    數(shù)值模擬結(jié)果會(huì)受到網(wǎng)格數(shù)量的影響,因此為了保證計(jì)算精準(zhǔn)性及計(jì)算效率并盡可能節(jié)約計(jì)算資源,需要進(jìn)行網(wǎng)格無關(guān)性驗(yàn)證。本文以1號(hào)物理模型為基準(zhǔn)進(jìn)行網(wǎng)格無關(guān)性驗(yàn)證,共劃分4種不同數(shù)量的網(wǎng)格,分別為363 022、600 894、1 102 244、1 681 768。在網(wǎng)格無關(guān)性驗(yàn)證過程中,需要保證除網(wǎng)格不同外其他模擬條件完全一致,最終通過計(jì)算得到測(cè)量段的表面?zhèn)鳠嵯禂?shù),見表2。由表2可知,以方案4為基準(zhǔn),方案1、方案2與方案4下的表面?zhèn)鳠嵯禂?shù)相差較大,而方案3與方案4下的表面?zhèn)鳠嵯禂?shù)相差很小。因此為了節(jié)約計(jì)算資源并保證計(jì)算效率,最終選擇網(wǎng)格數(shù)量為1 102 244的方案3進(jìn)行模擬計(jì)算分析。

    表2 網(wǎng)格無關(guān)性驗(yàn)證得到測(cè)量段的表面?zhèn)鳠嵯禂?shù)

    3 模擬結(jié)果及分析

    本文將重點(diǎn)研究螺旋管內(nèi)直徑、螺旋直徑及螺旋升角對(duì)超臨界甲烷在螺旋管內(nèi)的表面?zhèn)鳠嵯禂?shù)以及摩擦壓力降的影響。甲烷入口溫度范圍為200~300 K,入口壓力為7 MPa,質(zhì)量流率為300 kg/(m2·s-1)。螺旋管內(nèi)直徑為6、10、14 mm,螺旋直徑為1 000、2 000、3 000 mm,螺旋升角為2°、6°、10°。當(dāng)研究某個(gè)變量的影響時(shí),其他變量按此取值。但在探究螺旋管內(nèi)直徑與螺旋直徑對(duì)螺旋管表面?zhèn)鳠嵯禂?shù)的影響規(guī)律時(shí),將入口溫度上限提升至360 K。

    3.1 螺旋管內(nèi)直徑對(duì)換熱性能的影響

    ① 不同螺旋管內(nèi)直徑下,表面?zhèn)鳠嵯禂?shù)隨入口溫度的變化

    對(duì)螺旋直徑為2 000 mm,螺旋升角為10°,螺旋管內(nèi)直徑分別為6、10、14 mm下的超臨界甲烷在螺旋管內(nèi)流動(dòng)換熱特性進(jìn)行研究。圖5為螺旋管內(nèi)的超臨界甲烷在不同螺旋管內(nèi)直徑下的表面?zhèn)鳠嵯禂?shù)隨入口溫度變化曲線。由圖5可知,在螺旋管內(nèi)直徑相同的情況下,表面?zhèn)鳠嵯禂?shù)隨著入口溫度升高呈現(xiàn)先增加后降低的變化規(guī)律。認(rèn)為產(chǎn)生上述現(xiàn)象的原因是超臨界甲烷的比定壓熱容會(huì)隨溫度的增加先增大后減小并在擬臨界溫度處達(dá)到峰值,而其密度會(huì)隨溫度增加不斷減小。由于超臨界甲烷物性的變化,且比定壓熱容對(duì)換熱性能的影響程度高于密度對(duì)換熱性能的影響程度,從而引發(fā)上述現(xiàn)象。

    圖5 不同螺旋管內(nèi)直徑下表面?zhèn)鳠嵯禂?shù)隨入口溫度的變化

    當(dāng)入口溫度范圍為200~210 K時(shí),表面?zhèn)鳠嵯禂?shù)與螺旋管內(nèi)直徑呈正相關(guān)關(guān)系。當(dāng)入口溫度范圍為220~360 K時(shí),表面?zhèn)鳠嵯禂?shù)與螺旋管內(nèi)直徑呈負(fù)相關(guān)關(guān)系。認(rèn)為出現(xiàn)小部分溫度區(qū)間存在表面?zhèn)鳠嵯禂?shù)與螺旋管內(nèi)直徑呈正相關(guān)關(guān)系這一現(xiàn)象,主要原因是入口邊界設(shè)置為質(zhì)量流率入口,導(dǎo)致大管徑的管道質(zhì)量流量相應(yīng)較大,縮短了不同管徑之間表面?zhèn)鳠嵯禂?shù)的差距。

    隨著螺旋管內(nèi)直徑減小,表面?zhèn)鳠嵯禂?shù)達(dá)到峰值時(shí)的入口溫度有所增大且所達(dá)到的峰值越大,表面?zhèn)鳠嵯禂?shù)達(dá)到峰值時(shí)的入口溫度分別為210、215、225 K。認(rèn)為發(fā)生上述現(xiàn)象的原因是小管徑管道具有較強(qiáng)的流動(dòng)換熱性能,導(dǎo)致相同的入口溫度下小管徑管道的出口溫度較低。其中,螺旋管內(nèi)直徑為14 mm時(shí)表面?zhèn)鳠嵯禂?shù)峰值為4 078.87 W/(m2·K-1),螺旋管內(nèi)直徑為10 mm時(shí)表面?zhèn)鳠嵯禂?shù)峰值比14 mm時(shí)增長(zhǎng)了5.75%,螺旋管內(nèi)直徑為6 mm時(shí)表面?zhèn)鳠嵯禂?shù)峰值比10 mm時(shí)增長(zhǎng)了6.88%。

    綜上所述,螺旋管內(nèi)直徑對(duì)螺旋管流動(dòng)換熱性能有明顯影響,具體表現(xiàn)為超臨界甲烷在螺旋管內(nèi)流動(dòng)換熱性能會(huì)隨著螺旋管內(nèi)直徑減小而增強(qiáng)。這主要是因?yàn)楫?dāng)螺旋管內(nèi)直徑減小時(shí),表面張力增加導(dǎo)致螺旋管內(nèi)的冷凝換熱被強(qiáng)化。

    ② 不同螺旋管內(nèi)直徑下,比摩擦壓力降隨入口溫度的變化

    將單位長(zhǎng)度摩擦壓力降稱為比摩擦壓力降。圖6為螺旋管內(nèi)的超臨界甲烷在不同螺旋管內(nèi)直徑下的比摩擦壓力降隨入口溫度變化曲線。由圖6可知,在螺旋管內(nèi)直徑相同的情況下,比摩擦壓力降隨著入口溫度升高而升高。此外,比摩擦壓力降隨著螺旋管內(nèi)直徑減小而逐漸增大。當(dāng)入口溫度為300 K時(shí),螺旋管內(nèi)直徑為14 mm時(shí)的比摩擦壓力降為656.80 Pa/m,螺旋管內(nèi)直徑為10 mm時(shí)的比摩擦壓力降比14 mm時(shí)增長(zhǎng)了28.09%,螺旋管內(nèi)直徑為6 mm時(shí)的比摩擦壓力降比10 mm時(shí)增長(zhǎng)了44.04%。上述現(xiàn)象表明螺旋管內(nèi)直徑對(duì)螺旋管內(nèi)比摩擦壓力降影響較大。

    圖6 不同螺旋管內(nèi)直徑下比摩擦壓力降隨入口溫度的變化

    由于本文中的入口邊界選取質(zhì)量流率入口,使得管內(nèi)流速基本與管徑無關(guān),因此排除管徑變小導(dǎo)致管內(nèi)流速增大這一原因。引發(fā)上述現(xiàn)象的根本原因是螺旋管內(nèi)直徑減小使速度梯度有所增加,根據(jù)牛頓內(nèi)摩擦定律,速度梯度與黏性剪切力成正比。隨著螺旋管內(nèi)直徑減小,管內(nèi)流體與壁面間的黏性剪切力有所增加,導(dǎo)致比摩擦壓力降增大。

    為了更清晰地分析螺旋管內(nèi)的對(duì)流換熱過程,本文給出了螺旋直徑為2 000 mm、螺旋升角為10°、螺旋管內(nèi)直徑為10 mm的物理模型,在入口溫度為300 K工況下的螺旋管內(nèi)流體溫度及管壁溫度的沿程變化曲線,見圖7。其中,管程為位于測(cè)量段中沿管子軸線距測(cè)量段入口的距離。由圖7可知,管內(nèi)流體溫度從292.6 K降至269.0 K,管壁溫度從282.5 K降至260.3 K。從測(cè)量段入口到管道出口,管壁溫度與對(duì)應(yīng)的流體溫度的差值由10.1 K降低至8.7 K。

    圖7 螺旋管內(nèi)流體溫度及管壁溫度的沿程變化曲線

    3.2 螺旋直徑對(duì)換熱性能的影響

    對(duì)螺旋管內(nèi)直徑為10 mm,螺旋升角為10°,螺旋直徑分別為1 000、2 000、3 000 mm下的超臨界甲烷在螺旋管內(nèi)流動(dòng)換熱特性進(jìn)行研究。圖8為螺旋管內(nèi)的超臨界甲烷在不同螺旋直徑下的表面?zhèn)鳠嵯禂?shù)隨入口溫度變化曲線。分析圖8可知,在螺旋直徑相同的情況下,表面?zhèn)鳠嵯禂?shù)隨著入口溫度升高呈現(xiàn)先增加并達(dá)到峰值后降低的變化規(guī)律。隨著螺旋直徑增大,表面?zhèn)鳠嵯禂?shù)達(dá)到峰值時(shí)的入口溫度分別為210、215、220 K,峰值有所減小但幅度不大,其相對(duì)應(yīng)的峰值分別為4 389.53、4 313.36、4 270.27 W/(m2·K)。當(dāng)入口溫度范圍為220~290 K及350~360 K時(shí),表面?zhèn)鳠嵯禂?shù)與螺旋直徑呈正相關(guān)關(guān)系。當(dāng)入口溫度范圍為200~210 K及305~325 K時(shí),表面?zhèn)鳠嵯禂?shù)與螺旋直徑呈負(fù)相關(guān)關(guān)系。因?yàn)楸疚难芯康穆菪睆脚c管徑的比值過大,導(dǎo)致螺旋管的曲率效應(yīng)不明顯,使得螺旋直徑對(duì)超臨界甲烷螺旋管內(nèi)的表面?zhèn)鳠嵯禂?shù)影響規(guī)律不清晰。

    圖8 不同螺旋直徑下超臨界甲烷表面?zhèn)鳠嵯禂?shù)隨入口溫度的變化

    圖9為螺旋管內(nèi)的超臨界甲烷在不同螺旋直徑下的比摩擦壓力降隨入口溫度的變化曲線。由圖9可知,在螺旋管內(nèi)直徑、螺旋升角相同的情況下,入口溫度為200 K時(shí),3種螺旋直徑(分別為1 000、2 000、3 000 mm)下的比摩擦壓力降基本一致,而后隨著入口溫度升高不斷增加。此外,比摩擦壓力降隨著螺旋直徑增大逐漸降低,且螺旋直徑越小比摩擦壓力降的增長(zhǎng)率越大。當(dāng)入口溫度為300 K時(shí),螺旋直徑為3 000 mm的比摩擦壓力降為524.32 Pa/m,螺旋直徑為2 000 mm比3 000 mm時(shí)增長(zhǎng)了60.46%,螺旋直徑為1 000 mm比2 000 mm時(shí)增長(zhǎng)了66.01%。上述現(xiàn)象表明螺旋直徑對(duì)超臨界甲烷在螺旋管內(nèi)的比摩擦壓力降影響較大。其主要原因是由于螺旋直徑減小,離心力增強(qiáng),導(dǎo)致管內(nèi)二次流和湍流強(qiáng)度增加,進(jìn)而流動(dòng)阻力有所增大,比摩擦壓力降增大。

    圖9 不同螺旋直徑下比摩擦壓力降隨入口溫度的變化

    3.3 螺旋升角對(duì)換熱性能的影響

    對(duì)螺旋管內(nèi)直徑為10 mm,螺旋直徑為2 000 mm,螺旋升角分別為2°、6°、10°下的超臨界甲烷在螺旋管內(nèi)的流動(dòng)換熱特性進(jìn)行研究。圖10為螺旋管內(nèi)的超臨界甲烷在不同螺旋升角下的表面?zhèn)鳠嵯禂?shù)隨入口溫度的變化曲線。由圖10可知,不同螺旋升角下的表面?zhèn)鳠嵯禂?shù)基本一致,其中螺旋升角為2°的表面?zhèn)鳠嵯禂?shù)變化曲線中的突點(diǎn)認(rèn)為是模擬過程中的壞點(diǎn)。隨著螺旋升角增大,表面?zhèn)鳠嵯禂?shù)峰值分別為4 316.76、4 316.08、4 313.36 W/(m2·K)。因此認(rèn)為螺旋升角對(duì)超臨界甲烷在螺旋管內(nèi)的表面?zhèn)鳠嵯禂?shù)影響較小。

    圖10 不同螺旋升角下表面?zhèn)鳠嵯禂?shù)隨入口溫度的變化

    圖11為螺旋管內(nèi)的超臨界甲烷在不同螺旋升角下的比摩擦壓力降隨入口溫度變化曲線。由圖11可知,在其他條件相同的情況下,隨著入口溫度升高,不同螺旋升角的比摩擦壓力降間的差距緩慢降低。當(dāng)入口溫度為300 K時(shí),螺旋升角為2°的比摩擦壓力降為858.90 Pa/m,螺旋升角為6°比2°降低了1.60%,螺旋升角為10°比6°降低了0.48%。由上述現(xiàn)象可知,螺旋升角對(duì)超臨界甲烷在螺旋管內(nèi)的比摩擦壓力降有影響但影響效果不明顯。這是因?yàn)槁菪堑淖兓秶^小(2°~10°),使得離心力與重力之間的比例關(guān)系變化不明顯。

    圖11 不同螺旋升角下比摩擦壓力降隨入口溫度的變化

    3.4 出口橫截面速度及溫度分布

    ① 不同螺旋管內(nèi)直徑下的出口橫截面速度及溫度分布

    對(duì)入口溫度為300 K,螺旋直徑為2 000 mm,螺旋升角為10°,螺旋管內(nèi)直徑分別為6、10、14 mm下的出口橫截面速度及溫度分布進(jìn)行分析。圖12為不同螺旋管內(nèi)直徑下出口橫截面速度及流線分布。由圖12可知,出口橫截面速度呈現(xiàn)出中心大、四周小的分布情況,同時(shí)在管道兩側(cè)均有2個(gè)二次渦對(duì)稱分布。隨著螺旋管內(nèi)直徑增大,出口橫截面的平均流速逐漸增大,二次流強(qiáng)度逐漸增強(qiáng),但螺旋管內(nèi)直徑對(duì)二次渦結(jié)構(gòu)以及速度分布形式影響不大。

    圖12 不同螺旋管內(nèi)直徑下出口橫截面速度及流線分布(軟件截圖)

    不同螺旋管內(nèi)直徑下出口橫截面溫度分布見圖13。由圖13可知,管道外側(cè)溫度要小于管道內(nèi)側(cè)溫度。此外,隨著螺旋管內(nèi)直徑增大,管道內(nèi)的高溫流體所占橫截面面積越大且越貼近壁面,且橫截面最高溫度由250 K增長(zhǎng)至280 K。

    ② 不同螺旋直徑下的出口橫截面速度及溫度分布

    對(duì)入口溫度為300 K,螺旋管內(nèi)直徑為10 mm,螺旋升角為10°,螺旋直徑分別為1 000、2 000、3 000 mm下的出口橫截面速度及溫度分布進(jìn)行分析。圖14為不同螺旋直徑下出口橫截面速度及流線分布。由圖14可知,隨著螺旋直徑減小,出口橫截面流速增大且橫截面流速的分層現(xiàn)象越明顯。此外,管道內(nèi)二次流強(qiáng)度隨著螺旋直徑增大逐漸減弱且二次渦結(jié)構(gòu)有所改變。上述現(xiàn)象說明,隨著螺旋直徑減小,管內(nèi)流體的湍流強(qiáng)度不斷增強(qiáng)。因?yàn)樵谶\(yùn)行參數(shù)相同的情況下,螺旋直徑越小,螺旋管產(chǎn)生的離心力越大。

    圖13 不同螺旋管內(nèi)直徑下出口橫截面溫度分布(軟件截圖)

    圖15為不同螺旋直徑下出口橫截面溫度分布。由圖15可知,隨著螺旋直徑增大,管內(nèi)高溫流體所占橫截面面積逐漸減少且逐漸遠(yuǎn)離壁面,同時(shí)橫截面溫差由17 K降低至12 K。分析認(rèn)為管內(nèi)二次流強(qiáng)度減弱,導(dǎo)致管道出口橫截面溫度分布產(chǎn)生上述變化。

    圖14 不同螺旋直徑下出口橫截面速度及流線分布(軟件截圖)

    ③ 不同螺旋升角的出口橫截面速度及溫度分布

    對(duì)入口溫度為300 K,螺旋管內(nèi)直徑為10 mm,螺旋直徑為2 000 mm,螺旋升角分別為2°、6°、10°下的出口橫截面速度及溫度分布進(jìn)行分析。圖16、17分別為不同螺旋升角下出口橫截面速度與流線分布及出口橫截面溫度分布。由圖16、17可知,隨著螺旋升角變化,二次渦結(jié)構(gòu)及二次流強(qiáng)度、出口橫截面的流速分布以及溫度分布均無明顯變化。由此可見,螺旋升角變化不會(huì)對(duì)螺旋管內(nèi)流動(dòng)換熱性能產(chǎn)生過大影響。這是因?yàn)槁菪堑淖兓秶^小(2°~10°),使得離心力與重力之間的比例關(guān)系變化不明顯,所以發(fā)生上述現(xiàn)象。

    圖15 不同螺旋直徑下出口橫截面溫度分布(軟件截圖)

    圖16 不同螺旋升角下出口橫截面速度及流線分布(軟件截圖)

    圖17 不同螺旋升角下出口橫截面溫度分布(軟件截圖)

    4 結(jié)論

    ① 螺旋直徑及螺旋升角相同時(shí),隨著螺旋管內(nèi)直徑減小,表面?zhèn)鳠嵯禂?shù)達(dá)到峰值時(shí)的入口溫度及峰值、比摩擦壓力降均增大。

    ② 螺旋直徑及螺旋升角相同時(shí),隨著螺旋管內(nèi)直徑增大,出口橫截面的平均流速逐漸增大,管道內(nèi)的高溫流體所占橫截面面積越大且越貼近壁面,且橫截面最高溫度增大。

    ③ 螺旋管內(nèi)直徑及螺旋升角相同時(shí),隨著螺旋直徑增大,管內(nèi)表面?zhèn)鳠嵯禂?shù)峰值小幅度降低,不同入口溫度范圍表面?zhèn)鳠嵯禂?shù)與螺旋直徑呈現(xiàn)不同變化規(guī)律。此外,比摩擦壓力降隨著螺旋直徑增大而逐漸降低。

    ④ 螺旋管內(nèi)直徑及螺旋升角相同時(shí),隨著螺旋直徑減小,出口橫截面流速增大且橫截面流速的分層現(xiàn)象越明顯。隨著螺旋直徑增大,管內(nèi)高溫流體所占橫截面面積逐漸減少且逐漸遠(yuǎn)離壁面,同時(shí)橫截面溫差降低。

    ⑤ 螺旋升角對(duì)超臨界甲烷在螺旋管內(nèi)的流動(dòng)換熱性能影響較小。

    猜你喜歡
    升角螺旋管壓力降
    陽(yáng)極壓力降在PEMFC故障診斷中的應(yīng)用
    螺旋管內(nèi)氣液分離仿真分析*
    機(jī)載火焰抑制器流通性能研究
    基于ADAMS的采煤機(jī)截割機(jī)構(gòu)運(yùn)動(dòng)分析和結(jié)構(gòu)優(yōu)化
    螺旋升角對(duì)楔橫軋?zhí)菪温菁y軸成形的影響
    螺紋升角
    螺旋管預(yù)精焊技術(shù)在國(guó)內(nèi)的應(yīng)用
    焊接(2016年4期)2016-02-27 13:04:24
    氣體熱載體干餾爐內(nèi)壓力降的研究
    雙螺母滾珠絲杠副預(yù)緊力的研究
    船用膜式螺旋管換熱器傳熱和流動(dòng)特性研究
    久久精品影院6| 久久99热这里只有精品18| 麻豆乱淫一区二区| 久久人人精品亚洲av| 老司机影院成人| 搡女人真爽免费视频火全软件 | 美女被艹到高潮喷水动态| 久久韩国三级中文字幕| 亚洲欧美成人精品一区二区| 老司机影院成人| 国产国拍精品亚洲av在线观看| 亚洲综合色惰| 白带黄色成豆腐渣| 简卡轻食公司| 欧美性感艳星| 精华霜和精华液先用哪个| 午夜视频国产福利| 欧美极品一区二区三区四区| 国产成人影院久久av| 国产精品三级大全| 三级国产精品欧美在线观看| av在线老鸭窝| 国产av不卡久久| 美女被艹到高潮喷水动态| 日本精品一区二区三区蜜桃| 人人妻人人看人人澡| 亚洲四区av| 久久久成人免费电影| 亚洲真实伦在线观看| 亚洲,欧美,日韩| 精品一区二区三区视频在线| 国产成年人精品一区二区| 91久久精品电影网| 亚洲久久久久久中文字幕| 99热全是精品| 在线观看一区二区三区| 精品少妇黑人巨大在线播放 | 成人毛片a级毛片在线播放| 国产91av在线免费观看| 男人的好看免费观看在线视频| 蜜臀久久99精品久久宅男| 久久午夜亚洲精品久久| 国内精品美女久久久久久| 少妇的逼水好多| 18禁裸乳无遮挡免费网站照片| 少妇人妻一区二区三区视频| 国产综合懂色| 日韩亚洲欧美综合| 国产精品1区2区在线观看.| 国产真实乱freesex| 国产一区二区三区在线臀色熟女| 亚洲真实伦在线观看| 精品久久久久久久人妻蜜臀av| 麻豆国产av国片精品| 亚洲最大成人av| 亚洲国产高清在线一区二区三| 亚洲国产精品久久男人天堂| 精品人妻偷拍中文字幕| 长腿黑丝高跟| 欧美成人免费av一区二区三区| 久久精品综合一区二区三区| 中文字幕av成人在线电影| 99久国产av精品| 久久九九热精品免费| 天天躁夜夜躁狠狠久久av| 亚洲最大成人中文| 我要看日韩黄色一级片| 无遮挡黄片免费观看| av在线老鸭窝| 久久久国产成人免费| 亚洲一区二区三区色噜噜| a级一级毛片免费在线观看| 99热6这里只有精品| 狠狠狠狠99中文字幕| 综合色av麻豆| 久99久视频精品免费| 黄色日韩在线| 亚洲三级黄色毛片| 日韩欧美精品免费久久| 大香蕉久久网| 天堂网av新在线| 久久天躁狠狠躁夜夜2o2o| 色噜噜av男人的天堂激情| 在线a可以看的网站| 看片在线看免费视频| 在线免费观看不下载黄p国产| 国产一区二区三区在线臀色熟女| 99视频精品全部免费 在线| 免费一级毛片在线播放高清视频| 精品无人区乱码1区二区| 一个人观看的视频www高清免费观看| 国产女主播在线喷水免费视频网站 | 国产探花极品一区二区| 99久久精品国产国产毛片| 少妇猛男粗大的猛烈进出视频 | 亚洲久久久久久中文字幕| 中国美女看黄片| 国产一区二区在线观看日韩| 国产成人一区二区在线| 一进一出抽搐动态| videossex国产| 深爱激情五月婷婷| 91av网一区二区| av.在线天堂| 在线国产一区二区在线| 久久久久久久久久久丰满| 欧美精品国产亚洲| 亚洲美女视频黄频| 成人综合一区亚洲| 最近中文字幕高清免费大全6| 老熟妇乱子伦视频在线观看| 男女做爰动态图高潮gif福利片| 国语自产精品视频在线第100页| 性色avwww在线观看| 色5月婷婷丁香| 天美传媒精品一区二区| 人人妻人人澡欧美一区二区| 免费观看的影片在线观看| 国产成年人精品一区二区| 欧美一区二区精品小视频在线| 免费观看精品视频网站| 成人av在线播放网站| 成年版毛片免费区| 尤物成人国产欧美一区二区三区| 国产精品嫩草影院av在线观看| 久久久久久久午夜电影| 一级黄片播放器| 中文字幕久久专区| 在现免费观看毛片| 国产一区二区三区av在线 | 听说在线观看完整版免费高清| 国产伦精品一区二区三区四那| 日韩大尺度精品在线看网址| 国产黄色视频一区二区在线观看 | 午夜爱爱视频在线播放| 国内揄拍国产精品人妻在线| 你懂的网址亚洲精品在线观看 | 国产淫片久久久久久久久| 欧美性感艳星| 色综合亚洲欧美另类图片| 亚洲,欧美,日韩| 极品教师在线视频| 美女内射精品一级片tv| 一夜夜www| 色综合亚洲欧美另类图片| 日韩欧美精品v在线| 久久精品91蜜桃| 亚洲中文日韩欧美视频| 久久久久精品国产欧美久久久| 亚洲欧美清纯卡通| 亚洲无线在线观看| 国产免费一级a男人的天堂| 国产精品野战在线观看| 嫩草影院新地址| 久久精品夜夜夜夜夜久久蜜豆| 我的老师免费观看完整版| 免费观看精品视频网站| 国产精品永久免费网站| 成人二区视频| eeuss影院久久| 国产探花极品一区二区| 听说在线观看完整版免费高清| 九色成人免费人妻av| 国产黄a三级三级三级人| 日日啪夜夜撸| 一个人看的www免费观看视频| 国产大屁股一区二区在线视频| 精品久久久久久成人av| 91在线精品国自产拍蜜月| 国产蜜桃级精品一区二区三区| 一区福利在线观看| 少妇的逼水好多| 听说在线观看完整版免费高清| 亚洲自偷自拍三级| 精品午夜福利在线看| 高清日韩中文字幕在线| 国产免费男女视频| 麻豆国产av国片精品| 亚洲性夜色夜夜综合| 亚洲精品日韩在线中文字幕 | 黄色视频,在线免费观看| 亚洲熟妇中文字幕五十中出| 亚洲综合色惰| 久久久久久久久大av| 精品人妻熟女av久视频| 日日干狠狠操夜夜爽| 久久久久久大精品| 欧美成人一区二区免费高清观看| 校园人妻丝袜中文字幕| 一个人免费在线观看电影| 嫩草影视91久久| 精品国内亚洲2022精品成人| 直男gayav资源| 精品久久国产蜜桃| 又粗又爽又猛毛片免费看| 午夜日韩欧美国产| 欧美潮喷喷水| 欧美日本视频| 亚洲国产日韩欧美精品在线观看| 亚洲成av人片在线播放无| 两性午夜刺激爽爽歪歪视频在线观看| 99精品在免费线老司机午夜| 特大巨黑吊av在线直播| 99久国产av精品| 午夜福利成人在线免费观看| 国产精品久久久久久av不卡| 日本爱情动作片www.在线观看 | 国产一区二区三区av在线 | 日韩欧美在线乱码| 老师上课跳d突然被开到最大视频| 日本在线视频免费播放| 国产单亲对白刺激| 乱系列少妇在线播放| 一级黄片播放器| 简卡轻食公司| 国产中年淑女户外野战色| or卡值多少钱| 久久综合国产亚洲精品| 国产高清激情床上av| 在线观看一区二区三区| 国产伦精品一区二区三区视频9| avwww免费| 国产免费一级a男人的天堂| 免费观看在线日韩| 亚洲电影在线观看av| 婷婷精品国产亚洲av在线| 国产爱豆传媒在线观看| 午夜福利在线在线| 亚洲乱码一区二区免费版| 日本色播在线视频| 此物有八面人人有两片| 亚洲成人久久性| 一区二区三区高清视频在线| 天天躁夜夜躁狠狠久久av| 久久欧美精品欧美久久欧美| 12—13女人毛片做爰片一| 欧美日韩一区二区视频在线观看视频在线 | 岛国在线免费视频观看| 久久久色成人| 午夜免费激情av| 99九九线精品视频在线观看视频| 天天躁日日操中文字幕| 欧美成人一区二区免费高清观看| 亚洲熟妇熟女久久| 一级毛片电影观看 | 毛片一级片免费看久久久久| 性插视频无遮挡在线免费观看| 1024手机看黄色片| 亚洲成人久久爱视频| 99久久精品国产国产毛片| 日韩 亚洲 欧美在线| 成人av在线播放网站| 国产精品久久久久久久久免| 给我免费播放毛片高清在线观看| 久久午夜亚洲精品久久| 国产精品嫩草影院av在线观看| 一级黄片播放器| 亚洲第一电影网av| 亚洲人成网站在线观看播放| 国产免费男女视频| 亚洲中文日韩欧美视频| 国产精品一区二区免费欧美| 国产在线男女| 国产激情偷乱视频一区二区| 看黄色毛片网站| 精品一区二区三区av网在线观看| 中文字幕久久专区| 在线国产一区二区在线| 免费看光身美女| 久久精品91蜜桃| 亚洲国产高清在线一区二区三| 日韩,欧美,国产一区二区三区 | 九色成人免费人妻av| 麻豆国产97在线/欧美| 听说在线观看完整版免费高清| 国产成人a∨麻豆精品| 免费人成视频x8x8入口观看| 99久久精品国产国产毛片| 老司机福利观看| 国产精品久久久久久精品电影| 97超视频在线观看视频| 美女xxoo啪啪120秒动态图| 99国产精品一区二区蜜桃av| 99国产精品一区二区蜜桃av| 亚洲五月天丁香| 插阴视频在线观看视频| 一级毛片电影观看 | 久久精品国产亚洲av涩爱 | 久久久成人免费电影| 变态另类丝袜制服| 国内精品宾馆在线| 精品欧美国产一区二区三| 日韩 亚洲 欧美在线| 成年免费大片在线观看| 一进一出抽搐动态| 人妻少妇偷人精品九色| 男人狂女人下面高潮的视频| 亚洲图色成人| 欧美最黄视频在线播放免费| 色av中文字幕| 午夜激情欧美在线| 亚洲av电影不卡..在线观看| 亚洲在线观看片| 18+在线观看网站| .国产精品久久| 久久精品人妻少妇| avwww免费| 一级毛片电影观看 | 亚洲精品色激情综合| 久久精品综合一区二区三区| 女生性感内裤真人,穿戴方法视频| 欧美zozozo另类| 国产成人aa在线观看| 别揉我奶头 嗯啊视频| 亚洲中文日韩欧美视频| 久久久久国产精品人妻aⅴ院| 国产黄片美女视频| 日韩人妻高清精品专区| avwww免费| АⅤ资源中文在线天堂| 精品久久久久久久末码| 国产精品一区www在线观看| 熟女人妻精品中文字幕| 国产精品伦人一区二区| 国产精品不卡视频一区二区| 国产欧美日韩一区二区精品| 日本三级黄在线观看| 亚洲,欧美,日韩| 欧美日韩一区二区视频在线观看视频在线 | 国产精品久久电影中文字幕| 18+在线观看网站| 黄色一级大片看看| 国产极品精品免费视频能看的| 成人欧美大片| h日本视频在线播放| 国产高潮美女av| 成人三级黄色视频| eeuss影院久久| 国产成人福利小说| 18禁在线无遮挡免费观看视频 | 在线播放国产精品三级| 乱码一卡2卡4卡精品| 99久久九九国产精品国产免费| 亚洲激情五月婷婷啪啪| 美女黄网站色视频| 国产精品精品国产色婷婷| 丰满乱子伦码专区| 一边摸一边抽搐一进一小说| 一卡2卡三卡四卡精品乱码亚洲| 精品少妇黑人巨大在线播放 | 在线观看午夜福利视频| 色av中文字幕| av天堂中文字幕网| 久久精品国产亚洲av香蕉五月| 国产高清有码在线观看视频| 国产美女午夜福利| 国产精品野战在线观看| 国产三级中文精品| 亚洲七黄色美女视频| 久久久久国内视频| 色哟哟·www| 卡戴珊不雅视频在线播放| 日日干狠狠操夜夜爽| 久久久久性生活片| 亚洲精华国产精华液的使用体验 | 国产亚洲91精品色在线| 最近最新中文字幕大全电影3| 亚洲三级黄色毛片| 中文字幕av在线有码专区| 97超碰精品成人国产| 亚洲精华国产精华液的使用体验 | 亚洲18禁久久av| 一个人免费在线观看电影| 秋霞在线观看毛片| 国产精品永久免费网站| 熟妇人妻久久中文字幕3abv| 91久久精品电影网| 亚洲av成人av| 九九热线精品视视频播放| 国产高清视频在线观看网站| 一区二区三区四区激情视频 | 欧美成人免费av一区二区三区| 国产探花在线观看一区二区| 国产成人一区二区在线| 身体一侧抽搐| 午夜福利在线在线| 十八禁网站免费在线| 亚洲va在线va天堂va国产| 欧美在线一区亚洲| 亚洲国产色片| 又黄又爽又刺激的免费视频.| 一个人看的www免费观看视频| 婷婷精品国产亚洲av在线| 2021天堂中文幕一二区在线观| 亚洲国产欧洲综合997久久,| 久久天躁狠狠躁夜夜2o2o| 干丝袜人妻中文字幕| 久久人人爽人人片av| 成人综合一区亚洲| 欧美+日韩+精品| 91午夜精品亚洲一区二区三区| 国产精品国产三级国产av玫瑰| av国产免费在线观看| 国产精品国产三级国产av玫瑰| 一级毛片我不卡| 亚洲不卡免费看| 亚洲三级黄色毛片| 可以在线观看毛片的网站| 亚洲色图av天堂| 亚洲精品在线观看二区| 哪里可以看免费的av片| 一夜夜www| 秋霞在线观看毛片| 亚洲丝袜综合中文字幕| 日本熟妇午夜| 欧美最新免费一区二区三区| 日本与韩国留学比较| ponron亚洲| 天堂动漫精品| 亚洲成人久久爱视频| 国产 一区 欧美 日韩| 欧美日本视频| 精品福利观看| 日韩人妻高清精品专区| 九九久久精品国产亚洲av麻豆| av在线天堂中文字幕| 日韩欧美免费精品| 久久九九热精品免费| 精品人妻一区二区三区麻豆 | 亚洲aⅴ乱码一区二区在线播放| www日本黄色视频网| 69人妻影院| 亚洲性夜色夜夜综合| 俺也久久电影网| 国产精品久久久久久av不卡| 日韩高清综合在线| 久久午夜亚洲精品久久| 国产 一区 欧美 日韩| 人人妻人人看人人澡| 人人妻人人看人人澡| 91午夜精品亚洲一区二区三区| 亚洲熟妇熟女久久| 日韩亚洲欧美综合| 中文资源天堂在线| 亚洲中文字幕一区二区三区有码在线看| 日韩av不卡免费在线播放| 国产精品一区二区三区四区久久| 日日摸夜夜添夜夜爱| 精品久久久久久成人av| 国产伦一二天堂av在线观看| 国内少妇人妻偷人精品xxx网站| 天堂动漫精品| 热99re8久久精品国产| 日韩欧美在线乱码| 人人妻,人人澡人人爽秒播| 中国国产av一级| 少妇熟女欧美另类| 91在线观看av| 2021天堂中文幕一二区在线观| 天美传媒精品一区二区| 久久草成人影院| 亚洲中文字幕日韩| 91久久精品国产一区二区三区| 美女cb高潮喷水在线观看| 老师上课跳d突然被开到最大视频| 精品少妇黑人巨大在线播放 | 欧美人与善性xxx| 3wmmmm亚洲av在线观看| 亚洲国产精品成人综合色| 成人av在线播放网站| 日本熟妇午夜| 久久人人爽人人爽人人片va| 免费看日本二区| 亚洲第一区二区三区不卡| 深爱激情五月婷婷| 欧美国产日韩亚洲一区| 久久久久久九九精品二区国产| 亚洲av一区综合| 中文字幕人妻熟人妻熟丝袜美| 欧美又色又爽又黄视频| 中国美女看黄片| 国产亚洲精品综合一区在线观看| 久久99热这里只有精品18| 蜜桃亚洲精品一区二区三区| 国产高潮美女av| 欧美激情在线99| av黄色大香蕉| 最好的美女福利视频网| 六月丁香七月| 亚洲精华国产精华液的使用体验 | 日本黄色视频三级网站网址| 嫩草影院精品99| 久久人妻av系列| 一区福利在线观看| 国产美女午夜福利| 中文资源天堂在线| 热99在线观看视频| 秋霞在线观看毛片| 成人av一区二区三区在线看| 卡戴珊不雅视频在线播放| 国产精品乱码一区二三区的特点| 99在线人妻在线中文字幕| 亚洲美女搞黄在线观看 | 综合色丁香网| 日日摸夜夜添夜夜添av毛片| 国国产精品蜜臀av免费| 久久久a久久爽久久v久久| 91av网一区二区| 成人av在线播放网站| av福利片在线观看| 俄罗斯特黄特色一大片| 又黄又爽又刺激的免费视频.| 精品久久久久久久久亚洲| aaaaa片日本免费| 两个人的视频大全免费| 毛片一级片免费看久久久久| 中文资源天堂在线| 亚洲熟妇熟女久久| 国产在视频线在精品| 亚洲成人中文字幕在线播放| 免费看av在线观看网站| 免费在线观看影片大全网站| 高清毛片免费观看视频网站| 亚洲,欧美,日韩| 欧美性猛交黑人性爽| 国产aⅴ精品一区二区三区波| 少妇被粗大猛烈的视频| 国产视频内射| 婷婷精品国产亚洲av在线| 在现免费观看毛片| 国产精品1区2区在线观看.| 综合色丁香网| 不卡视频在线观看欧美| 久久精品综合一区二区三区| 国产精品一区二区免费欧美| 久久99热这里只有精品18| 国产精品永久免费网站| av国产免费在线观看| 久久精品国产亚洲网站| 日韩欧美 国产精品| 久久久久九九精品影院| 免费无遮挡裸体视频| 综合色丁香网| 97超碰精品成人国产| 久久久精品94久久精品| 成人美女网站在线观看视频| 亚洲欧美成人综合另类久久久 | 日韩欧美国产在线观看| 中文字幕av成人在线电影| 毛片一级片免费看久久久久| 一个人看视频在线观看www免费| 国产一区二区三区在线臀色熟女| 免费观看精品视频网站| 久久久久久久久久黄片| 在线播放国产精品三级| 色综合亚洲欧美另类图片| 亚洲人成网站在线播| 国产一区二区三区av在线 | 国模一区二区三区四区视频| 国产黄片美女视频| 狠狠狠狠99中文字幕| 天天一区二区日本电影三级| 有码 亚洲区| 啦啦啦韩国在线观看视频| 深夜精品福利| 精品久久久久久久久久免费视频| 精品一区二区免费观看| 亚洲图色成人| 久久精品国产清高在天天线| 亚洲精品一卡2卡三卡4卡5卡| 在线观看av片永久免费下载| 禁无遮挡网站| 日本成人三级电影网站| 午夜老司机福利剧场| 一个人看的www免费观看视频| 校园人妻丝袜中文字幕| 色5月婷婷丁香| 午夜亚洲福利在线播放| 亚洲av二区三区四区| 18禁裸乳无遮挡免费网站照片| 免费av不卡在线播放| 一边摸一边抽搐一进一小说| 在线观看66精品国产| 久久天躁狠狠躁夜夜2o2o| 成年版毛片免费区| 色播亚洲综合网| 亚洲经典国产精华液单| 亚洲五月天丁香| 国产精品1区2区在线观看.| 男人的好看免费观看在线视频| 精品不卡国产一区二区三区| 男人的好看免费观看在线视频| 亚洲人成网站在线播放欧美日韩| 97热精品久久久久久| 在线a可以看的网站| 一级毛片我不卡| 高清午夜精品一区二区三区 | 天堂√8在线中文| 成人综合一区亚洲| 亚洲内射少妇av| 亚洲av中文字字幕乱码综合| 久久久精品94久久精品| 真人做人爱边吃奶动态| 国产精品日韩av在线免费观看| 国产精品野战在线观看| 国产精品一区二区三区四区久久| 成人美女网站在线观看视频| 久久精品影院6| 免费一级毛片在线播放高清视频| 精品久久国产蜜桃| 精品久久久久久久人妻蜜臀av| 久久午夜福利片| 小说图片视频综合网站| 久久国内精品自在自线图片| 男女啪啪激烈高潮av片|