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

    換熱器螺旋管冷凝換熱的數(shù)值模擬與實(shí)驗(yàn)研究

    2023-12-14 13:03:30徐裕民孫偉珂于哲林洪俊陳學(xué)永李西兵
    科學(xué)技術(shù)與工程 2023年32期
    關(guān)鍵詞:螺旋管制冷劑冷凝

    徐裕民, 孫偉珂, 于哲, 林洪俊, 陳學(xué)永, 李西兵

    (福建農(nóng)林大學(xué)機(jī)電工程學(xué)院, 福州 350002)

    隨著工業(yè)化進(jìn)程的加快,能源問(wèn)題也越來(lái)越突出,中國(guó)在第75屆聯(lián)合國(guó)大會(huì)上宣布,力爭(zhēng)2030年前二氧化碳排放達(dá)到峰值,努力爭(zhēng)取2060年前實(shí)現(xiàn)碳中和目標(biāo)[1]。在兼顧可靠性與碳中和兩個(gè)目標(biāo)的要求下,提高能效,注重節(jié)能技術(shù)的開(kāi)發(fā)、推廣及使用。換熱器是一種應(yīng)用非常廣泛的熱交換器,在能源產(chǎn)業(yè)中有著舉足輕重的地位,換熱器在整個(gè)生產(chǎn)工藝總設(shè)備中占據(jù)30%以上。具有單位面積換熱量大,傳熱效率高,操作壓力大,補(bǔ)償性能好等優(yōu)點(diǎn),現(xiàn)廣泛應(yīng)用于制冷、冶金和化工等領(lǐng)域[2]。因此加強(qiáng)不同參數(shù)對(duì)螺旋管傳熱的影響研究對(duì)于提高換熱器的節(jié)能效果非常重要。

    對(duì)于換熱器的冷凝換熱,中外學(xué)者開(kāi)展了大量工作。楊英英等[3]采用R32為工質(zhì)對(duì)水平管進(jìn)行可視化冷凝換熱實(shí)驗(yàn)研究,結(jié)果表明,El Hajal J流型圖對(duì)小管徑內(nèi)冷凝流型的預(yù)測(cè)誤差較大。Mozafari等[4]開(kāi)展了螺旋管換熱器內(nèi)R-600a的冷凝特性實(shí)驗(yàn)研究,結(jié)果表明,在每一傾角下,換熱系數(shù)隨質(zhì)量流量和平均蒸汽質(zhì)量的增加而增加。Macdonald等[5]對(duì)丙烷在水平管內(nèi)的冷凝傳熱進(jìn)行了實(shí)驗(yàn),結(jié)果表明,傳熱系數(shù)隨著質(zhì)量流量的增加而增加,隨著飽和溫度的降低而降低。Salimpour等[6]討論了結(jié)構(gòu)參數(shù)對(duì)螺旋盤管中R404A冷凝傳熱的影響,結(jié)果表明,在低蒸汽質(zhì)量下,傳熱系數(shù)隨著盤管直徑的減小而增加,隨著盤管間距的增加而減小。Yu等[7]以碳?xì)漕惢旌瞎べ|(zhì)作為制冷劑,使用實(shí)驗(yàn)和模擬相結(jié)合的研究方法對(duì)螺旋管內(nèi)的凝結(jié)換熱開(kāi)展研究,提出了一種適用于碳?xì)漕惢旌瞎べ|(zhì)的冷凝換熱關(guān)聯(lián)式。Zhu等[8]對(duì)光管和螺旋槽管兩種結(jié)構(gòu)中超臨界CO2的換熱系數(shù)進(jìn)行了實(shí)驗(yàn)研究,結(jié)果表明,螺旋槽管具有較好的換熱效率。趙進(jìn)元[9]對(duì)螺旋管內(nèi)二次流的換熱特性進(jìn)行研究,最終得到了螺旋管內(nèi)強(qiáng)化換熱技術(shù)的方法。Koluib[10]對(duì)螺旋管換熱器R407-C熱流體工質(zhì)的冷凝換熱特性進(jìn)行實(shí)驗(yàn)研究,結(jié)果表明,R407-C的質(zhì)量流量與換熱系數(shù)成正比,隨著R407-C飽和溫度的降低,R407-C側(cè)換熱系數(shù)和整體換熱系數(shù)都相應(yīng)增加。曹佳銘等[11]對(duì)一種纏繞管換熱器進(jìn)行結(jié)構(gòu)優(yōu)化與實(shí)驗(yàn)研究,提出一整套改進(jìn)的計(jì)算模型,結(jié)果表明,改進(jìn)后的模型傳熱誤差最大不超過(guò)3.6%,壓降誤差最大不超過(guò)8.5%,說(shuō)明改進(jìn)的計(jì)算模型具有較高的準(zhǔn)確度。

    綜上可知,關(guān)于水平管和螺旋管冷凝傳熱的實(shí)驗(yàn)研究很多,而對(duì)換熱器螺旋管的數(shù)值模擬研究很少,對(duì)螺旋管在制冷和空調(diào)熱水器系統(tǒng)中的推廣使用而言,現(xiàn)有研究還遠(yuǎn)遠(yuǎn)不足。為此,對(duì)R22制冷劑在螺旋管中冷凝換熱進(jìn)行數(shù)值模擬,分析換熱器螺旋管的流場(chǎng)分布特性,研究流體流速和飽和溫度對(duì)螺旋管內(nèi)換熱性能的影響,并在企業(yè)內(nèi)的換熱器測(cè)試實(shí)驗(yàn)室進(jìn)行實(shí)驗(yàn),對(duì)數(shù)值模擬的準(zhǔn)確性進(jìn)行驗(yàn)證。為螺旋管換熱器的設(shè)計(jì)優(yōu)化以及空調(diào)熱水器一體機(jī)的節(jié)能損耗研究提供一定的科學(xué)依據(jù)。

    1 數(shù)學(xué)模擬方法

    1.1 數(shù)學(xué)模型

    1.1.1 多相流模型控制方程

    流體體積法(volume of fluid,VOF)通常可用于跟蹤氣液界面的運(yùn)動(dòng),是一種簡(jiǎn)單有效的自由邊界數(shù)值處理方法,因此選用VOF模型作為模擬設(shè)置中的多相流模型。FLUENT18.0中的VOF模型能夠用于模擬不同類型管道的冷凝和流動(dòng)[12]。對(duì)于多相流,在VOF模型中,αl和αv分別為液體和氣體的體積分?jǐn)?shù)。液體和氣體的總體積分?jǐn)?shù)應(yīng)限制為1,即

    αl+αv=1

    (1)

    其他熱性質(zhì)(如黏度和導(dǎo)熱率)也通過(guò)上述方程所示的算術(shù)平均值來(lái)計(jì)算。能量和溫度計(jì)算公式分別為

    (2)

    (3)

    式(2)中:El為液體內(nèi)能量;Ev為氣體內(nèi)能量;ρv為氣體密度,kg/m3;ρl為液體密度,kg/m3;Tl為液體的溫度,K;Tv為氣體的溫度,K;密度ρ的計(jì)算公式為

    ρ=αlρl+αvρv

    (4)

    包括連續(xù)性、動(dòng)量、能量守恒在內(nèi)的控制方程如下[13]。

    (1)連續(xù)性方程。

    (5)

    式(5)中:t為時(shí)間,s;?為哈密爾頓算子;u為速度,m/s;Sl為從相變模型獲得的液相質(zhì)量源項(xiàng),kg/(m3s)。

    (6)

    (2)動(dòng)量守恒方程。

    ρg+Fσ

    (7)

    式(7)中:p為壓力,Pa;ɑ為相體積分?jǐn)?shù);μ為物體黏度,Pa·s;g為重力加速度,m/s2;Fσ為表面張力。由連續(xù)表面力模型計(jì)算得出[14],其計(jì)算公式為

    (8)

    式(8)中:σ為表面張力系數(shù);κ為曲率,其計(jì)算公式為

    (9)

    (3)能量守恒方程。

    (10)

    式(10)中:E氣體和液體的總能量;T為氣體和液體的總溫度;λe為氣液相有效熱傳導(dǎo)系數(shù);hlg為汽化潛熱。

    1.1.2 相變模型

    在數(shù)值研究中,Lee相變模型同時(shí)考慮了傳熱和傳質(zhì)。質(zhì)量源項(xiàng)由Lee相變模型[15]計(jì)算。Lee相變模型與VOF模型相結(jié)合已廣泛應(yīng)用于管內(nèi)流動(dòng)冷凝的建模,VOF模型與Lee模型相結(jié)合用于模擬,以獲得流型和傳熱系數(shù)。Lee模型中的傳質(zhì)可表示為

    (11)

    (12)

    式中:Tsat為飽和溫度;Tl為液體溫度;Tv為氣體溫度;Sv為蒸發(fā)中的傳質(zhì);Sl為冷凝中的傳質(zhì);r為松弛因子;ρl和ρv分別為液相和氣相密度,eva為蒸發(fā),con為冷凝。

    Lee相變模型中的松弛因子r對(duì)模擬結(jié)果有較大的影響[16]。因此,為了提高收斂速度和精度,在冷凝換熱的數(shù)值模擬中,r值取103。

    1.1.3 湍流模型

    使用可以旋流校正的“Realizablek-ε模型”對(duì)湍流條件下?lián)Q熱器螺旋管的管內(nèi)流場(chǎng)進(jìn)行數(shù)值仿真能更好的適應(yīng)曲率變化,而且可以加快收斂速度?!癛ealizablek-ε模型”的湍動(dòng)能k的表達(dá)式為[17-18]

    (13)

    式(13)中:σk為湍動(dòng)能的湍流普朗特?cái)?shù);xi為X軸方向坐標(biāo)分量;ui、uj為X、Y軸方向速度分量。。

    湍流耗散率ε的表達(dá)式為

    (14)

    式(14)中:C1、C2為經(jīng)驗(yàn)常數(shù);xj為Y軸方向坐標(biāo)分量;σε為耗散率的湍流普朗特?cái)?shù)。

    湍流黏性系數(shù)μt可表示為

    (15)

    式(15)中:Cμ為常數(shù)。

    1.2 物理模型

    采用Solidworks2018三維軟件構(gòu)建了換熱器螺旋管的物理模型,重要的結(jié)構(gòu)參數(shù)主要有換熱器螺旋管的螺旋直徑D、螺距P和管道直徑d,如圖1所示。利用ICEM CFD18.0軟件對(duì)換熱器螺旋管物理模型進(jìn)行網(wǎng)格劃分,換熱器螺旋管的進(jìn)口端面的網(wǎng)格,如圖2所示,通過(guò)FLUENT18.0軟件對(duì)換熱器螺旋管的物理模型進(jìn)行數(shù)值仿真。

    圖1 螺旋管三維結(jié)構(gòu)示意圖Fig.1 Schematic diagram of three dimensional structure of spiral tube

    圖2 進(jìn)口面網(wǎng)格示意圖Fig.2 Grid diagram of inlet surface

    1.3 數(shù)值方法

    數(shù)值模擬中考慮了重力的影響,重力加速度方向是沿Z軸負(fù)方向。換熱器螺旋管的管內(nèi)工質(zhì)為R22制冷劑,計(jì)算域內(nèi)的工質(zhì)物性可由REFPROP9.1確定。在螺旋管的進(jìn)口處設(shè)置速度入口邊界條件,螺旋管進(jìn)口溫度設(shè)置為323.15 K;在螺旋管的出口處設(shè)置壓力出口邊界條件;設(shè)置入口螺旋管內(nèi)的氣相體積分?jǐn)?shù)為1;換熱器螺旋管的壁面設(shè)置為無(wú)滑移壁面,材料為不銹鋼,由于定壁溫或者定熱流等壁面邊界條件對(duì)壓降和換熱無(wú)影響[19],因此螺旋管壁面設(shè)置為定壁溫?zé)o滑移壁面,壁面溫度設(shè)置為313.15 K,與飽和溫度差值10 K,輸入相應(yīng)的水力直徑和湍流強(qiáng)度。

    在FLUENT軟件中,采用Transient求解計(jì)算器,Explicit格式離散求解。采用Implicity Body Force處理提高計(jì)算收斂性。湍流模型選用Realizablek-ε模型,選擇Scalable Wall Functions來(lái)處理湍流模型的邊界層。求解計(jì)算過(guò)程采用SIMPLE算法;梯度離散項(xiàng)采用Least Squares Cell Based;壓力離散項(xiàng)采用PRESTO!;體積分?jǐn)?shù)項(xiàng)采用Geo-Reconstruction;動(dòng)量方程、能量方程、湍流動(dòng)能方程和耗散率方程采用二階迎風(fēng)格式;松弛因子不變,使用默認(rèn)值,對(duì)管道出口處氣相體積分?jǐn)?shù)、進(jìn)出口質(zhì)量流量進(jìn)行監(jiān)控,保證數(shù)值模擬達(dá)到收斂。

    2 數(shù)值模擬結(jié)果與分析

    2.1 網(wǎng)格無(wú)關(guān)性驗(yàn)證

    網(wǎng)格無(wú)關(guān)性驗(yàn)證將選取結(jié)構(gòu)參數(shù)為螺旋直徑D=300 mm、螺距P=19.52 mm、管道直徑d=9.52 mm的換熱器螺旋管,對(duì)螺旋管的物理模型采用O型結(jié)構(gòu)化網(wǎng)格劃分,并且對(duì)邊界層處的網(wǎng)格進(jìn)行加密,表1為相應(yīng)的網(wǎng)格設(shè)置參數(shù)與網(wǎng)格數(shù)量。

    表1 網(wǎng)格劃分參數(shù)及輸出網(wǎng)格數(shù)Table 1 Mesh partition parameters and output mesh number

    網(wǎng)格無(wú)關(guān)性驗(yàn)證基于FLUENT18.0軟件進(jìn)行分析,在軟件中設(shè)置條件為:工質(zhì)為水,進(jìn)口速度為1 m/s,進(jìn)口溫度為298.15 K,在恒壁溫的條件下,壁面溫度為333.15 K。換熱系數(shù)和進(jìn)出口壓降隨網(wǎng)格數(shù)量的變化分別如圖3(a)和圖3(b)所示。

    圖3 換熱系數(shù)和進(jìn)出口總壓降隨網(wǎng)格數(shù)量的變化示意圖Fig.3 Schematic diagram of wall heat transfer coefficient and total pressure drop at inlet and outlet with the number of grids

    數(shù)值仿真結(jié)果表明,換熱器螺旋管的換熱系數(shù)、進(jìn)出口壓降隨網(wǎng)格數(shù)量增加而減少;當(dāng)網(wǎng)格數(shù)量為4 166 400的情況下,換熱器螺旋管的換熱系數(shù)和進(jìn)出口壓降的變化趨于平穩(wěn),將不隨網(wǎng)格的變化而變化。因此從計(jì)算精度和效率兩方面考慮,采用第4組劃分網(wǎng)格的方法,這種劃分網(wǎng)格的方法可以滿足網(wǎng)格無(wú)關(guān)性驗(yàn)證的要求,對(duì)其它螺旋管進(jìn)行網(wǎng)格劃分時(shí),均可使用類似的方法計(jì)算。

    2.2 各參數(shù)沿管截面分布情況

    2.2.1 螺旋管內(nèi)氣液相體積分布情況

    圖4為流體在螺旋管內(nèi)沿著軸向方向上,是ZX截面的氣液相分布云圖。

    1#~22#為螺旋管入口開(kāi)始每180°的截面

    R22制冷劑在螺旋管內(nèi)的冷凝過(guò)程中,隨著氣態(tài)R22制冷劑的不斷凝結(jié),冷凝液逐漸積累,管內(nèi)的冷凝流型開(kāi)始改變。從圖4可以看出,沿管道方向冷凝液越來(lái)越多,主要流型為分層流。因?yàn)槿肟谶吔鐥l件為純氣相,冷凝時(shí)間非常短暫,因此,入口階段流型為霧狀流。隨著氣態(tài)R22釋放潛熱,在管壁附近持續(xù)冷凝,在剪切力的影響下,冷凝液被管道中心流速較快的氣體吹起,由于剪切力大于重力,因此,冷凝液沿著螺旋管外壁流動(dòng),形成半環(huán)狀流,隨著冷凝過(guò)程的繼續(xù)進(jìn)行,流型逐漸由半環(huán)狀流轉(zhuǎn)化成分層流,在分層流流型中,由于液相密度較大,速度也相應(yīng)減小,流型在重力的作用下,冷凝液聚集在管道底部。此外,還可以看出氣液界面略微傾斜,這是由于離心力導(dǎo)致的。

    2.2.2 螺旋管內(nèi)溫度場(chǎng)分布情況

    圖5為流體在螺旋管內(nèi)沿著軸向方向上并且是ZX截面的溫度場(chǎng)分布云圖。

    1#~22#為螺旋管入口開(kāi)始每180°的截面

    將螺旋管的溫度場(chǎng)與氣液相分布云圖相對(duì)比,螺旋管管內(nèi)溫度場(chǎng)與氣液相分布場(chǎng)情況對(duì)應(yīng)良好,從圖5可以看出,入口處全部為高溫的氣態(tài)工質(zhì),管壁溫度低于制冷劑溫度發(fā)生冷凝,冷凝換熱發(fā)生不久后,螺旋管外壁開(kāi)始出現(xiàn)一層液膜,因此管壁附近溫度較低。隨著氣態(tài)R22制冷劑在螺旋管內(nèi)發(fā)生持續(xù)性的冷凝,在重力和離心力的雙重影響下,管道上部的冷凝液開(kāi)始流向管道的外側(cè)底部,并且冷凝液逐漸增多,導(dǎo)致管底的溫度比較低。因此,與管道頂部相比,管道底部有很明顯的溫度邊界層。

    2.2.3 螺旋管內(nèi)速度場(chǎng)分布情況

    圖6為流體在螺旋管內(nèi)沿著軸向方向上并且是ZX截面的溫度場(chǎng)分布云圖。

    1#~22#為螺旋管入口開(kāi)始每180°的截面

    從圖6中可以看出,流體速度沿著管道流動(dòng)方向逐漸減小,這是因?yàn)殡S著冷凝換熱地進(jìn)行,冷凝液逐漸增多,流體密度逐漸增大,整體流速也就逐漸減小,并且液相部分的速度小于氣相部分的速度。此外,在螺旋管的上端和下端流體的流速則是近似對(duì)稱的,并且速度分布沒(méi)有明顯的變化,由于離心力和二次回流的作用,截面外側(cè)的流速要大于內(nèi)側(cè)的流速,外側(cè)管壁附近具有較大的速度梯度,管內(nèi)側(cè)的速度等值線呈以內(nèi)側(cè)頂點(diǎn)為中心的小的凸起狀。

    2.3 螺旋管內(nèi)冷凝換熱影響因素

    2.3.1 流體流速對(duì)換熱系數(shù)的影響

    圖7為在飽和溫度323.15 K下,螺旋管內(nèi)不同流體流速6、8、10 m/s時(shí)冷凝換熱系數(shù)隨干度的變化趨勢(shì)圖。從圖7中可以看出,在飽和溫度和干度一樣的條件下,流體流速越高,冷凝換熱系數(shù)越大。因?yàn)殡S著流體流速增高,氣液界面的剪切力越大,增加了界面擾動(dòng),促進(jìn)了液膜中的湍流強(qiáng)度,強(qiáng)化了氣相和液相之間的換熱,使得冷凝換熱系數(shù)增大。此外,在相同速度的條件下,干度越大,換熱系數(shù)也相應(yīng)增大。因?yàn)楦啥仍黾?說(shuō)明螺旋管道的截面有很大一部分為氣相,液相部分占據(jù)的較少,冷凝液的厚度隨著干度的增加而減小,其中的熱阻也相應(yīng)減小,因此使得換熱系數(shù)增加。

    Tsat為飽和溫度;v為流體速度

    2.3.2 飽和溫度對(duì)換熱系數(shù)的影響

    圖8為流體流速在8 m/s下,螺旋管內(nèi)不同飽和溫度303.15、313.15、323.15 K時(shí)冷凝換熱系數(shù)隨干度的變化趨勢(shì)圖。從圖8中可以看出,飽和溫度對(duì)冷凝換熱系數(shù)有很大的影響,隨著飽和溫度降低,冷凝換熱系數(shù)相應(yīng)增大,因?yàn)樵诓煌娘柡蜏囟葧r(shí)R22制冷劑有著不同的物理性質(zhì),飽和溫度越低時(shí),R22制冷劑釋放出的潛熱越高,換熱系數(shù)越大。另外,隨著飽和溫度的增高,液相密度減小,而氣相密度增大,氣液相密度的改變會(huì)造成氣體空隙率降低,讓管道截面上更大一部分被液相占據(jù),產(chǎn)生非常大熱阻,從而使得冷凝換熱系數(shù)降低。R22制冷劑在303.15、313.15、323.15 K時(shí)的熱力學(xué)性質(zhì)使用物性軟件REFPROP9.1查詢,R22制冷劑的物性參數(shù)如表2所示,汽化潛熱由標(biāo)準(zhǔn)狀態(tài)焓與摩爾質(zhì)量計(jì)算,計(jì)算公式為[20]

    表2 液相R22和氣相R22在各飽和溫度下的熱物性參數(shù)Table 2 Thermophysical parameters of liquid phase R22 and gas phase R22 at various saturation temperatures

    圖8 飽和溫度對(duì)冷凝換熱系數(shù)的影響Fig.8 Influence of saturation temperature on heat transfer coefficient of condensation

    (16)

    式(16)中:L為汽化潛熱,J/kg;H為標(biāo)準(zhǔn)狀態(tài)焓,J/kmol;M為摩爾質(zhì)量,kg/kmol。

    從表2中還可以看出,在液相狀態(tài)下,隨著飽和溫度增高,液體的黏度降低,使得液體的導(dǎo)熱能力升高,有助于提高螺旋管的換熱性能,但是液相熱阻對(duì)換熱系數(shù)影響更大。因此,在飽和溫度升高時(shí),整體的冷凝換熱系數(shù)降低。

    3 實(shí)驗(yàn)研究

    3.1 實(shí)驗(yàn)系統(tǒng)原理

    本次實(shí)驗(yàn)為螺旋管管內(nèi)冷凝傳熱實(shí)驗(yàn),實(shí)驗(yàn)的主要目的是為了對(duì)不同因素下螺旋管內(nèi)冷凝換熱系數(shù)進(jìn)行測(cè)量,從而對(duì)數(shù)值模擬結(jié)果進(jìn)行驗(yàn)證。采用R22制冷劑作為冷凝換熱的工質(zhì),在常壓的工況下運(yùn)行。將不銹鋼制作的螺旋管作為研究對(duì)象,實(shí)驗(yàn)系統(tǒng)原理如圖9所示,主要由R22制冷劑回路、冷卻水回路和相關(guān)測(cè)試裝置以及數(shù)據(jù)測(cè)量系統(tǒng)組成。

    圖9 實(shí)驗(yàn)原理圖Fig.9 Experimental schematic diagram

    制冷劑系統(tǒng)主要包括壓縮機(jī)、循環(huán)泵、螺旋管實(shí)驗(yàn)段、儲(chǔ)液罐、過(guò)濾閥、流量計(jì)、膨脹閥、蒸發(fā)器、氣液分離器和穩(wěn)壓器,流量計(jì)置于過(guò)濾閥與膨脹閥之間對(duì)制冷劑的流量進(jìn)行測(cè)量。壓縮機(jī)排出的高溫、高壓的氣態(tài)R22制冷劑進(jìn)入換熱器螺旋管實(shí)驗(yàn)段,與注入生活水箱中的冷水進(jìn)行換熱,發(fā)生冷凝產(chǎn)生冷凝液流進(jìn)儲(chǔ)液罐中,被干燥、計(jì)量后經(jīng)膨脹閥節(jié)流進(jìn)入蒸發(fā)器,吸收周圍環(huán)境的熱量,再經(jīng)過(guò)氣液分離器,變?yōu)闅鈶B(tài)制冷劑后被壓縮機(jī)吸入,進(jìn)入下一次循環(huán)。同時(shí),在螺旋管實(shí)驗(yàn)段的進(jìn)出口分別插入PT100溫度傳感器和壓力傳感器對(duì)R22的溫度和壓力進(jìn)行測(cè)量。并且通過(guò)流量計(jì)測(cè)量R22體積流量計(jì)算速度。

    實(shí)驗(yàn)整體裝置如圖10所示,實(shí)驗(yàn)將選取結(jié)構(gòu)參數(shù)為螺旋直徑D=300 mm、螺距P=19.52 mm、管道直徑d=9.52 mm的換熱器螺旋管。為了降低環(huán)境溫度所產(chǎn)生的誤差,所有實(shí)驗(yàn)設(shè)備和管道表面使用了隔熱保溫措施。

    圖10 實(shí)驗(yàn)實(shí)物圖Fig.10 Actual picture of experiment

    3.2 實(shí)驗(yàn)數(shù)據(jù)處理

    實(shí)驗(yàn)過(guò)程中保持螺旋管外壁面溫度為恒定值,通過(guò)調(diào)節(jié)流量大小,再基于螺旋管進(jìn)出口溫度和壓力,計(jì)算出螺旋管的壁面平均熱流密度q,在相變情況下q的計(jì)算公式為[21]

    (17)

    式(17)中:m為質(zhì)量,kg;l為螺旋管長(zhǎng)度,mm;ΔH為螺旋管進(jìn)出口焓差,J/kmol。

    螺旋管實(shí)驗(yàn)段進(jìn)口干度可表示為

    (18)

    螺旋管實(shí)驗(yàn)段出口干度可表示為

    (19)

    式中:xin、xout分別為進(jìn)口和出口的干度;U為加熱電壓,V;I為加熱電流,A;Q為冷凝放熱量,W;hs為R22處于飽和液態(tài)時(shí)的比焓,J/kg;hin為預(yù)熱前過(guò)冷液態(tài)時(shí)的比焓,J/kg;mr為制冷劑測(cè)質(zhì)量流量,kg/h。

    R22制冷劑的平均干度可表示為

    (20)

    R22制冷劑的飽和溫度可表示為

    (21)

    式(21)中:Tin、Tout分別為R22制冷劑進(jìn)出口溫度,K。

    螺旋管內(nèi)冷凝換熱系數(shù)的計(jì)算公式為[22-23]

    (22)

    式(22)中:hex為R22制冷劑的冷凝換熱系數(shù),W/(m2·K);δw為管壁厚度,mm;Tw為壁面溫度,K;λ為螺旋管的導(dǎo)熱系數(shù),W/(m·K)。

    3.3 數(shù)值模擬結(jié)果與實(shí)驗(yàn)結(jié)果對(duì)比分析

    將模擬數(shù)據(jù)與實(shí)驗(yàn)數(shù)據(jù)進(jìn)行對(duì)比,并判斷數(shù)值模擬模型的合理性。

    3.3.1 不同流體流速下冷凝換熱模擬數(shù)據(jù)與實(shí)驗(yàn)數(shù)據(jù)對(duì)比

    圖11描述了流體流速分別為6、8、10 m/s,飽和溫度323.15 K的條件下,數(shù)值模擬和實(shí)驗(yàn)得到的冷凝換熱系數(shù)對(duì)比情況。從圖11中可以看出,流體流速模擬數(shù)據(jù)和實(shí)驗(yàn)數(shù)據(jù)對(duì)冷凝換熱換熱系數(shù)的影響的趨勢(shì)較為一致,冷凝換熱系數(shù)在模擬數(shù)據(jù)整體上比實(shí)驗(yàn)數(shù)據(jù)偏高,流體流速為6 m/s時(shí),最大誤差在11%,平均相對(duì)誤差為7.02%,流體流速為8 m/s時(shí),最大誤差在10.8%,平均相對(duì)誤差為5%,流體流速為10 m/s時(shí),最大誤差在7.02%,平均相對(duì)誤差為4.7%,流速越大時(shí),誤差較小,因?yàn)樵诹黧w流速較大時(shí),傳熱傳質(zhì)得到強(qiáng)化,換熱的損失量較小。

    圖11 不同流體流速時(shí)冷凝換熱系數(shù)模擬值與實(shí)驗(yàn)值的對(duì)比Fig.11 Comparison of simulated and experimental values of condensation heat transfer coefficient at different flow rates

    3.3.2 不同飽和溫度下冷凝換熱模擬數(shù)據(jù)與實(shí)驗(yàn)數(shù)據(jù)對(duì)比

    圖12給出了在流體流速8 m/s,飽和溫度分別為303.15、313.15、323.15 K時(shí),數(shù)值模擬和實(shí)驗(yàn)得到的冷凝換熱系數(shù)對(duì)比情況。模擬數(shù)據(jù)和實(shí)驗(yàn)數(shù)據(jù)均表明,R22在303.15 K時(shí)的冷凝換熱系數(shù)整體上高于323.15 K的冷凝換熱系數(shù)。飽和溫度為303.15K時(shí),最大誤差在3.26%,平均相對(duì)誤差為1.84%,飽和溫度為313.15 K時(shí),最大誤差在6.05%,平均相對(duì)誤差為3.1%,飽和溫度為323.15 K時(shí),最大誤差在8.11%,平均相對(duì)誤差為3.2%。將模擬數(shù)據(jù)和實(shí)驗(yàn)數(shù)據(jù)進(jìn)行對(duì)比,產(chǎn)生的誤差在10%以內(nèi),說(shuō)明了所建模型的準(zhǔn)確性。此外,實(shí)驗(yàn)和模擬的結(jié)果都說(shuō)明,在流體流速一定時(shí),增加干度對(duì)于提高冷凝換熱系數(shù)有十分明顯的作用。

    根據(jù)對(duì)換熱器螺旋管所做的數(shù)值模擬以及實(shí)驗(yàn)研究可知,對(duì)數(shù)值模擬結(jié)果與實(shí)驗(yàn)結(jié)果進(jìn)行對(duì)比,在不同流體流速時(shí)冷凝換熱系數(shù)的模擬數(shù)據(jù)與實(shí)驗(yàn)數(shù)據(jù)之間的相對(duì) 誤差為3%~11%,在不同飽和溫度時(shí)冷凝換熱系數(shù)的模擬數(shù)據(jù)與實(shí)驗(yàn)數(shù)據(jù)之間的相對(duì)誤差為3%~8%,誤差工程計(jì)算的允許范圍之內(nèi),說(shuō)明數(shù)值模擬方法和結(jié)果是合理的。通過(guò)對(duì)流體流速和飽和溫度模擬數(shù)據(jù)和實(shí)驗(yàn)數(shù)據(jù)進(jìn)行對(duì)比,發(fā)現(xiàn)模擬數(shù)據(jù)均高于實(shí)驗(yàn)數(shù)據(jù),產(chǎn)生誤差的原因包括:在數(shù)值模擬中假設(shè)管道內(nèi)壁為光滑的狀態(tài),但在實(shí)驗(yàn)當(dāng)中管道內(nèi)壁有一定的粗糙度,使工質(zhì)的湍動(dòng)能發(fā)生耗散,導(dǎo)致在實(shí)驗(yàn)過(guò)程中管內(nèi)工質(zhì)的實(shí)際速度低于設(shè)置的進(jìn)口速度;模擬中假設(shè)溫度是不發(fā)生變化的,而在實(shí)驗(yàn)過(guò)程中,溫度是發(fā)生變化的,工質(zhì)的熱物性也會(huì)產(chǎn)生變化。

    4 結(jié)論

    對(duì)R22制冷劑在螺旋管中冷凝換熱進(jìn)行了數(shù)值模擬和實(shí)驗(yàn)研究,分析了換熱器螺旋管的流場(chǎng)分布特性,研究了流體流速和飽和溫度對(duì)螺旋管內(nèi)換熱性能的影響,最終通過(guò)實(shí)驗(yàn)對(duì)數(shù)值模擬的準(zhǔn)確性進(jìn)行了驗(yàn)證。得出如下主要結(jié)論。

    (1)通過(guò)數(shù)值仿真分析,螺旋管管內(nèi)溫度場(chǎng)、速度場(chǎng)和氣液相分布場(chǎng)的情況對(duì)應(yīng)良好

    (2)通過(guò)對(duì)換熱器螺旋管的數(shù)值仿真分析,增大流體流速,可以提高換熱器螺旋管冷凝換熱性能,而飽和溫度越低,螺旋管的冷凝換熱性能越好。同時(shí),增大螺旋管的干度能夠加強(qiáng)螺旋管的換熱系數(shù)。

    (3)根據(jù)對(duì)螺旋管管內(nèi)的冷凝傳熱實(shí)驗(yàn),結(jié)果表明:在不同流體流速時(shí)冷凝換熱系數(shù)的模擬數(shù)據(jù)與實(shí)驗(yàn)數(shù)據(jù)之間的相對(duì)誤差為3%~11%,在不同飽和溫度時(shí)冷凝換熱系數(shù)的模擬數(shù)據(jù)與實(shí)驗(yàn)數(shù)據(jù)之間的相對(duì)誤差為3%~8%,誤差工程計(jì)算的允許范圍之內(nèi),說(shuō)明數(shù)值模擬方法和結(jié)果是合理的。

    猜你喜歡
    螺旋管制冷劑冷凝
    揚(yáng)州市出臺(tái)“制冷劑六條”提升制冷劑安全管理水平
    洗脫苯粗苯冷凝冷卻器優(yōu)化改造
    山東冶金(2019年6期)2020-01-06 07:46:14
    螺旋管內(nèi)氣液分離仿真分析*
    近共沸制冷劑R134a/R1234yf的PVTx性質(zhì)的實(shí)驗(yàn)研究
    硫磺回收裝置冷凝冷卻器改造
    螺旋管預(yù)精焊技術(shù)在國(guó)內(nèi)的應(yīng)用
    焊接(2016年4期)2016-02-27 13:04:24
    船用膜式螺旋管換熱器傳熱和流動(dòng)特性研究
    壓縮機(jī)級(jí)間冷卻冷凝析水量計(jì)算
    霍尼韋爾新型制冷劑獲歐盟科研機(jī)構(gòu)認(rèn)可
    美國(guó)擬對(duì)中國(guó)制冷劑產(chǎn)品征收反補(bǔ)貼稅
    丝袜在线中文字幕| 人成视频在线观看免费观看| 99热国产这里只有精品6| √禁漫天堂资源中文www| 男男h啪啪无遮挡| 久久久久国产精品人妻一区二区| 午夜视频国产福利| 午夜日本视频在线| 亚洲国产精品一区三区| 国产精品国产三级国产专区5o| 日本91视频免费播放| 国产精品.久久久| 日韩三级伦理在线观看| 99久久精品国产国产毛片| 亚洲图色成人| 你懂的网址亚洲精品在线观看| 这个男人来自地球电影免费观看 | 一区二区三区乱码不卡18| 婷婷色av中文字幕| 国产一区二区激情短视频 | 久久精品熟女亚洲av麻豆精品| 青春草视频在线免费观看| 免费少妇av软件| 国产又色又爽无遮挡免| 女人久久www免费人成看片| 日本av免费视频播放| av天堂久久9| 久久久久国产精品人妻一区二区| av视频免费观看在线观看| 男女下面插进去视频免费观看 | 久久久欧美国产精品| 天天操日日干夜夜撸| 亚洲精品乱久久久久久| 午夜激情av网站| 美女中出高潮动态图| 最近中文字幕2019免费版| 一级片免费观看大全| 日韩人妻精品一区2区三区| 大香蕉久久成人网| 一区二区三区精品91| 日韩不卡一区二区三区视频在线| 亚洲av福利一区| 国产男女内射视频| 亚洲av福利一区| 成人18禁高潮啪啪吃奶动态图| 在线观看免费高清a一片| 亚洲国产精品成人久久小说| 免费看av在线观看网站| 亚洲第一av免费看| 久久精品人人爽人人爽视色| 国产精品国产三级国产av玫瑰| 一区二区三区乱码不卡18| 极品少妇高潮喷水抽搐| 高清黄色对白视频在线免费看| 亚洲国产欧美在线一区| 9热在线视频观看99| 日韩精品免费视频一区二区三区 | 中文字幕av电影在线播放| 人妻少妇偷人精品九色| 国产日韩欧美亚洲二区| 久久国产精品男人的天堂亚洲 | 五月玫瑰六月丁香| 9色porny在线观看| 精品卡一卡二卡四卡免费| 国产精品久久久久久精品电影小说| 亚洲成人一二三区av| 亚洲三级黄色毛片| 性高湖久久久久久久久免费观看| 狠狠婷婷综合久久久久久88av| 黄片无遮挡物在线观看| 在线观看三级黄色| 91在线精品国自产拍蜜月| 丰满乱子伦码专区| 日韩成人av中文字幕在线观看| 91在线精品国自产拍蜜月| 国产精品熟女久久久久浪| 欧美精品亚洲一区二区| 国产无遮挡羞羞视频在线观看| 乱码一卡2卡4卡精品| 亚洲国产精品999| 22中文网久久字幕| 极品人妻少妇av视频| 亚洲成人手机| 日日啪夜夜爽| 纵有疾风起免费观看全集完整版| 亚洲国产av新网站| 欧美国产精品一级二级三级| 黄色毛片三级朝国网站| 最近2019中文字幕mv第一页| 99热这里只有是精品在线观看| 97在线人人人人妻| 亚洲国产精品999| 久久久久久久久久久久大奶| av又黄又爽大尺度在线免费看| 亚洲国产精品成人久久小说| 新久久久久国产一级毛片| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 国产精品蜜桃在线观看| 男人添女人高潮全过程视频| 成人亚洲欧美一区二区av| 亚洲av国产av综合av卡| 成人国语在线视频| 久久久久久久精品精品| 国产老妇伦熟女老妇高清| 三上悠亚av全集在线观看| 国产精品不卡视频一区二区| 久久女婷五月综合色啪小说| 男女啪啪激烈高潮av片| 一级毛片 在线播放| 日韩成人av中文字幕在线观看| 国产片特级美女逼逼视频| 伊人久久国产一区二区| 亚洲精品美女久久久久99蜜臀 | 在线观看www视频免费| 亚洲在久久综合| 久久精品久久久久久噜噜老黄| 久久午夜福利片| 老女人水多毛片| 寂寞人妻少妇视频99o| 91aial.com中文字幕在线观看| 男女边吃奶边做爰视频| 日韩免费高清中文字幕av| 欧美成人精品欧美一级黄| 久久久久精品人妻al黑| xxx大片免费视频| 一二三四在线观看免费中文在 | 夫妻午夜视频| 天堂俺去俺来也www色官网| 亚洲精品美女久久久久99蜜臀 | 久久人人爽人人片av| 亚洲成av片中文字幕在线观看 | 国产男女内射视频| 久久久久久人人人人人| 满18在线观看网站| 久久精品久久精品一区二区三区| 亚洲第一av免费看| 国产精品偷伦视频观看了| 亚洲精品久久午夜乱码| 涩涩av久久男人的天堂| 亚洲,一卡二卡三卡| 18禁动态无遮挡网站| 久久久久国产网址| 精品久久久久久电影网| 久久综合国产亚洲精品| 人妻系列 视频| 久久亚洲国产成人精品v| 国产精品偷伦视频观看了| 日韩成人av中文字幕在线观看| 免费黄网站久久成人精品| 国产高清不卡午夜福利| 日本猛色少妇xxxxx猛交久久| 国产精品久久久久久精品电影小说| 久热这里只有精品99| 国产精品国产三级专区第一集| 国产综合精华液| 美女视频免费永久观看网站| 99国产综合亚洲精品| 亚洲一级一片aⅴ在线观看| 欧美 亚洲 国产 日韩一| 国产淫语在线视频| av免费观看日本| 成人无遮挡网站| 国产日韩欧美在线精品| 欧美精品一区二区大全| 久久青草综合色| 亚洲av在线观看美女高潮| 99视频精品全部免费 在线| 欧美日韩成人在线一区二区| 亚洲av成人精品一二三区| 夜夜骑夜夜射夜夜干| 国产又色又爽无遮挡免| 中国三级夫妇交换| 亚洲精品456在线播放app| 99视频精品全部免费 在线| 欧美另类一区| 国产亚洲av片在线观看秒播厂| 美女大奶头黄色视频| 日韩成人av中文字幕在线观看| 精品亚洲成a人片在线观看| 国产精品 国内视频| 亚洲色图综合在线观看| 赤兔流量卡办理| 午夜日本视频在线| 在线观看美女被高潮喷水网站| 如何舔出高潮| 男女啪啪激烈高潮av片| 91精品国产国语对白视频| 巨乳人妻的诱惑在线观看| 免费高清在线观看视频在线观看| 人妻人人澡人人爽人人| 捣出白浆h1v1| 高清av免费在线| 激情视频va一区二区三区| 亚洲国产日韩一区二区| 97在线视频观看| 九九爱精品视频在线观看| 午夜福利网站1000一区二区三区| 亚洲成人手机| 亚洲精品美女久久久久99蜜臀 | 国产成人一区二区在线| 在线观看三级黄色| 亚洲图色成人| 亚洲熟女精品中文字幕| 婷婷色综合www| 欧美激情 高清一区二区三区| 少妇人妻久久综合中文| 免费观看无遮挡的男女| 亚洲三级黄色毛片| 亚洲综合色惰| 欧美3d第一页| 久久久亚洲精品成人影院| 亚洲,一卡二卡三卡| 精品一区二区三卡| 日本色播在线视频| 亚洲五月色婷婷综合| 90打野战视频偷拍视频| 亚洲国产精品一区三区| 精品国产一区二区三区四区第35| 成人黄色视频免费在线看| 免费日韩欧美在线观看| 伦理电影大哥的女人| 大片电影免费在线观看免费| 黑人巨大精品欧美一区二区蜜桃 | 麻豆精品久久久久久蜜桃| 夫妻性生交免费视频一级片| 国产片内射在线| av黄色大香蕉| 久久精品国产自在天天线| 曰老女人黄片| 丝瓜视频免费看黄片| 大陆偷拍与自拍| 美女xxoo啪啪120秒动态图| 新久久久久国产一级毛片| 大香蕉久久成人网| 狂野欧美激情性bbbbbb| 亚洲欧美精品自产自拍| 久久 成人 亚洲| 成人午夜精彩视频在线观看| 我的女老师完整版在线观看| 国产成人欧美| 蜜臀久久99精品久久宅男| 日韩免费高清中文字幕av| 十八禁网站网址无遮挡| 国产精品一二三区在线看| 插逼视频在线观看| 视频中文字幕在线观看| 欧美变态另类bdsm刘玥| av播播在线观看一区| 一本久久精品| 色94色欧美一区二区| 制服人妻中文乱码| 亚洲丝袜综合中文字幕| 七月丁香在线播放| 国产深夜福利视频在线观看| 看非洲黑人一级黄片| av女优亚洲男人天堂| 尾随美女入室| 成人国产av品久久久| 综合色丁香网| 宅男免费午夜| 男女边吃奶边做爰视频| 亚洲三级黄色毛片| 少妇的丰满在线观看| 中文字幕制服av| 视频在线观看一区二区三区| 最近中文字幕高清免费大全6| 国产高清三级在线| 一级,二级,三级黄色视频| 一二三四中文在线观看免费高清| 亚洲,欧美,日韩| 2018国产大陆天天弄谢| 涩涩av久久男人的天堂| 精品久久久久久电影网| 九九爱精品视频在线观看| 在线观看免费日韩欧美大片| 97精品久久久久久久久久精品| av.在线天堂| 啦啦啦视频在线资源免费观看| 宅男免费午夜| 国产高清国产精品国产三级| 亚洲人成网站在线观看播放| 国产亚洲欧美精品永久| 久久国产精品男人的天堂亚洲 | 久热这里只有精品99| 日韩大片免费观看网站| av黄色大香蕉| 日本vs欧美在线观看视频| 性色av一级| 精品福利永久在线观看| 成年人午夜在线观看视频| 亚洲经典国产精华液单| 久久韩国三级中文字幕| 少妇高潮的动态图| 日日摸夜夜添夜夜爱| 亚洲,欧美精品.| 精品久久国产蜜桃| 国产免费又黄又爽又色| 免费不卡的大黄色大毛片视频在线观看| 午夜精品国产一区二区电影| 成人手机av| 97人妻天天添夜夜摸| 欧美日韩一区二区视频在线观看视频在线| 少妇猛男粗大的猛烈进出视频| 尾随美女入室| 久久综合国产亚洲精品| 人妻一区二区av| 精品久久久久久电影网| 亚洲五月色婷婷综合| 成人毛片60女人毛片免费| 亚洲国产最新在线播放| 99国产综合亚洲精品| 少妇猛男粗大的猛烈进出视频| 久久久久久久久久成人| 制服人妻中文乱码| 国产精品不卡视频一区二区| 久久影院123| 80岁老熟妇乱子伦牲交| 免费久久久久久久精品成人欧美视频 | 在线天堂中文资源库| 久久ye,这里只有精品| 日韩电影二区| 日本猛色少妇xxxxx猛交久久| 国产精品久久久av美女十八| 亚洲成国产人片在线观看| 国产深夜福利视频在线观看| 国产一区二区在线观看av| 午夜精品国产一区二区电影| 国产极品粉嫩免费观看在线| 亚洲国产成人一精品久久久| 精品国产乱码久久久久久小说| 成人二区视频| 国产精品偷伦视频观看了| 免费黄色在线免费观看| 国产极品天堂在线| 麻豆精品久久久久久蜜桃| 亚洲精品乱码久久久久久按摩| a级毛片黄视频| 中文字幕另类日韩欧美亚洲嫩草| 一区二区三区乱码不卡18| 亚洲成色77777| 国产1区2区3区精品| 99国产精品免费福利视频| 青春草亚洲视频在线观看| 日韩一区二区视频免费看| 亚洲综合色网址| 一个人免费看片子| 精品午夜福利在线看| 一个人免费看片子| 啦啦啦在线观看免费高清www| 欧美精品高潮呻吟av久久| 亚洲国产av新网站| 精品亚洲成国产av| 一级毛片 在线播放| 97在线人人人人妻| 少妇的丰满在线观看| 丝袜喷水一区| 各种免费的搞黄视频| 国产亚洲精品第一综合不卡 | 久久毛片免费看一区二区三区| 成人国语在线视频| 欧美97在线视频| 亚洲精品色激情综合| 一个人免费看片子| 国产国拍精品亚洲av在线观看| 亚洲国产毛片av蜜桃av| 黑人高潮一二区| 国产高清三级在线| 2018国产大陆天天弄谢| 少妇的逼水好多| 啦啦啦在线观看免费高清www| 巨乳人妻的诱惑在线观看| 亚洲av福利一区| 夜夜爽夜夜爽视频| 日本猛色少妇xxxxx猛交久久| 日韩电影二区| 成人影院久久| 91精品三级在线观看| 精品少妇久久久久久888优播| 亚洲精品美女久久久久99蜜臀 | 热99久久久久精品小说推荐| 在线观看国产h片| 蜜桃在线观看..| 国产视频首页在线观看| 亚洲,欧美,日韩| 亚洲精品久久午夜乱码| 色94色欧美一区二区| 好男人视频免费观看在线| 成人亚洲欧美一区二区av| 性色avwww在线观看| 国产有黄有色有爽视频| 夜夜骑夜夜射夜夜干| www.av在线官网国产| 中文字幕最新亚洲高清| 国产探花极品一区二区| 嫩草影院入口| 国产精品久久久久久久电影| 成年人午夜在线观看视频| 交换朋友夫妻互换小说| 国产乱来视频区| 亚洲激情五月婷婷啪啪| 亚洲av综合色区一区| 哪个播放器可以免费观看大片| 国产一级毛片在线| 亚洲欧洲国产日韩| 91成人精品电影| 在线观看三级黄色| 久久久欧美国产精品| 亚洲国产精品一区三区| 黄色一级大片看看| av电影中文网址| 日本黄大片高清| a 毛片基地| 亚洲一码二码三码区别大吗| 亚洲成人手机| 国产在线免费精品| 两个人看的免费小视频| 日韩视频在线欧美| 日韩大片免费观看网站| 国产深夜福利视频在线观看| 亚洲精品视频女| 国内精品宾馆在线| 少妇高潮的动态图| www.熟女人妻精品国产 | 亚洲一级一片aⅴ在线观看| 亚洲国产成人一精品久久久| 婷婷色麻豆天堂久久| 精品一区二区三区视频在线| 全区人妻精品视频| 男女国产视频网站| 熟妇人妻不卡中文字幕| 免费人成在线观看视频色| av又黄又爽大尺度在线免费看| 纯流量卡能插随身wifi吗| freevideosex欧美| a 毛片基地| 欧美精品人与动牲交sv欧美| 久久国内精品自在自线图片| 亚洲人成77777在线视频| 观看美女的网站| 久久久久久久亚洲中文字幕| 69精品国产乱码久久久| 午夜激情av网站| 成人毛片60女人毛片免费| 一本大道久久a久久精品| 国产白丝娇喘喷水9色精品| 最近中文字幕高清免费大全6| 日韩成人伦理影院| 多毛熟女@视频| 男人添女人高潮全过程视频| 女人精品久久久久毛片| 国产深夜福利视频在线观看| 看非洲黑人一级黄片| 秋霞在线观看毛片| 日韩av在线免费看完整版不卡| 在现免费观看毛片| 成年美女黄网站色视频大全免费| 国语对白做爰xxxⅹ性视频网站| 在线天堂最新版资源| 最新的欧美精品一区二区| 国产成人免费无遮挡视频| 亚洲五月色婷婷综合| 母亲3免费完整高清在线观看 | 国产精品久久久久久久电影| av卡一久久| av有码第一页| 久久韩国三级中文字幕| 亚洲av电影在线进入| 久久 成人 亚洲| 亚洲经典国产精华液单| 9191精品国产免费久久| 90打野战视频偷拍视频| 亚洲欧美清纯卡通| 寂寞人妻少妇视频99o| 午夜福利影视在线免费观看| 你懂的网址亚洲精品在线观看| 久久人人爽av亚洲精品天堂| 久久av网站| 在线观看美女被高潮喷水网站| 国产精品不卡视频一区二区| 国产日韩一区二区三区精品不卡| 宅男免费午夜| 午夜免费鲁丝| av有码第一页| 99久久综合免费| 欧美日韩综合久久久久久| 王馨瑶露胸无遮挡在线观看| 国产高清不卡午夜福利| 午夜福利影视在线免费观看| 欧美xxⅹ黑人| 国产 精品1| 1024视频免费在线观看| 欧美精品亚洲一区二区| 成人毛片a级毛片在线播放| 亚洲综合精品二区| 欧美xxⅹ黑人| 国产片特级美女逼逼视频| 日韩一本色道免费dvd| 亚洲欧美中文字幕日韩二区| 亚洲一区二区三区欧美精品| 亚洲欧美中文字幕日韩二区| 欧美精品人与动牲交sv欧美| 51国产日韩欧美| av网站免费在线观看视频| 热99国产精品久久久久久7| 国产成人精品无人区| 水蜜桃什么品种好| 精品熟女少妇av免费看| 免费日韩欧美在线观看| 天天躁夜夜躁狠狠久久av| 亚洲婷婷狠狠爱综合网| 性高湖久久久久久久久免费观看| 啦啦啦视频在线资源免费观看| 久久毛片免费看一区二区三区| 久久久国产一区二区| 久久精品aⅴ一区二区三区四区 | 亚洲精品,欧美精品| 韩国av在线不卡| 亚洲精品日韩在线中文字幕| 欧美成人午夜精品| 久久久久久久久久成人| 美女xxoo啪啪120秒动态图| 亚洲精品一二三| 搡老乐熟女国产| 国产探花极品一区二区| 国产av国产精品国产| 伦理电影大哥的女人| 丰满少妇做爰视频| 中文字幕另类日韩欧美亚洲嫩草| 一本久久精品| 成人毛片60女人毛片免费| 国产日韩欧美视频二区| 91午夜精品亚洲一区二区三区| 久久99热6这里只有精品| 男人操女人黄网站| 久久精品国产鲁丝片午夜精品| 婷婷色综合大香蕉| 欧美成人午夜免费资源| 日韩中文字幕视频在线看片| 赤兔流量卡办理| 91在线精品国自产拍蜜月| 亚洲成人av在线免费| 在现免费观看毛片| 汤姆久久久久久久影院中文字幕| 亚洲国产欧美日韩在线播放| 国产一区二区在线观看av| 一级爰片在线观看| 有码 亚洲区| 熟女av电影| 青春草亚洲视频在线观看| 岛国毛片在线播放| 久久热在线av| 少妇高潮的动态图| 一本大道久久a久久精品| 大话2 男鬼变身卡| 少妇 在线观看| 久久99蜜桃精品久久| 亚洲美女搞黄在线观看| 亚洲一级一片aⅴ在线观看| 久久精品aⅴ一区二区三区四区 | 精品一区二区免费观看| 久久精品熟女亚洲av麻豆精品| 在线观看美女被高潮喷水网站| 99精国产麻豆久久婷婷| 亚洲人与动物交配视频| 五月伊人婷婷丁香| 又黄又爽又刺激的免费视频.| 午夜久久久在线观看| 深夜精品福利| 咕卡用的链子| 亚洲精品日韩在线中文字幕| 中文字幕免费在线视频6| 爱豆传媒免费全集在线观看| 蜜桃在线观看..| 亚洲成人手机| 在线看a的网站| 热re99久久精品国产66热6| 亚洲av.av天堂| 国产不卡av网站在线观看| 嫩草影院入口| 亚洲一码二码三码区别大吗| 超碰97精品在线观看| 亚洲av男天堂| 欧美变态另类bdsm刘玥| 91精品国产国语对白视频| 天天操日日干夜夜撸| 女人被躁到高潮嗷嗷叫费观| 国产精品人妻久久久影院| 国产精品人妻久久久久久| av卡一久久| 男女啪啪激烈高潮av片| 久久久久久人妻| 久久久国产一区二区| 日本黄色日本黄色录像| 在现免费观看毛片| 亚洲熟女精品中文字幕| av在线老鸭窝| 两个人看的免费小视频| 最近中文字幕2019免费版| 久久人人爽人人片av| 日本vs欧美在线观看视频| 国产又色又爽无遮挡免| 精品卡一卡二卡四卡免费| 日本-黄色视频高清免费观看| 91aial.com中文字幕在线观看| 激情视频va一区二区三区| 99久久中文字幕三级久久日本| 久久亚洲国产成人精品v| 久久久久久久精品精品| 日本黄大片高清| 国产精品蜜桃在线观看| 又粗又硬又长又爽又黄的视频| 91精品三级在线观看| 午夜福利影视在线免费观看| 一区二区三区精品91| av在线老鸭窝|