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

    墊條對(duì)纏繞管換熱器殼側(cè)性能的影響研究

    2016-12-23 01:27:08文鍵童欣楊輝著王斯民厲彥忠
    關(guān)鍵詞:纏繞管換熱器軸向

    文鍵,童欣,楊輝著,王斯民,厲彥忠

    (1.西安交通大學(xué)能源與動(dòng)力學(xué)院,710049,西安;2.西安交通大學(xué)化學(xué)工程與技術(shù)學(xué)院,710049,西安)

    ?

    墊條對(duì)纏繞管換熱器殼側(cè)性能的影響研究

    文鍵1,童欣1,楊輝著1,王斯民2,厲彥忠1

    (1.西安交通大學(xué)能源與動(dòng)力學(xué)院,710049,西安;2.西安交通大學(xué)化學(xué)工程與技術(shù)學(xué)院,710049,西安)

    目前針對(duì)纏繞管換熱器的數(shù)值模擬研究中,大多為了簡(jiǎn)化模型而省略了管層間的墊條,造成模擬結(jié)果偏離實(shí)際情況。通過CFD數(shù)值模擬,著重研究了纏繞管換熱器中墊條的存在及不同的排布方式和數(shù)量對(duì)換熱器殼側(cè)的流動(dòng)、傳熱性能的影響。計(jì)算結(jié)果表明,墊條的存在不僅影響流體軸向流速,并且在徑向截面上更易形成渦旋而加強(qiáng)流體擾動(dòng)。在相同進(jìn)口工況下,相比無(wú)墊條模型,墊條對(duì)齊排布時(shí)Nu下降了1.6%~2.5%,換熱器單位長(zhǎng)度壓降ΔPl增加了14.5%~17.0%,強(qiáng)化傳熱性能指標(biāo)FP減小了5.7%~5.9%;錯(cuò)開排布時(shí)Nu上升了17.0%~18.1%,ΔPl增加了38.4%~39.9%,FP增加了13.7%~15.4%。兩種墊條排布方式中錯(cuò)開排布能夠在更大程度上強(qiáng)化殼側(cè)的湍流,從而強(qiáng)化傳熱。當(dāng)墊條數(shù)量多于12,依次為13、14、15時(shí),對(duì)換熱器殼側(cè)的性能影響不大。所做研究為纏繞管換熱器的數(shù)值模擬提供了新的思路,也為纏繞管換熱器的性能優(yōu)化提供了理論參考。

    纏繞管換熱器;墊條;湍流強(qiáng)化;壓降

    纏繞管換熱器是一種典型的管殼式換熱器,它由換熱盤管纏繞中心筒組成,相鄰兩層管螺旋方向相反。換熱管內(nèi)流體按照螺旋形式流動(dòng),從而出現(xiàn)二次環(huán)流強(qiáng)化換熱,殼側(cè)流體在換熱管外和內(nèi)外殼體間的空隙交叉流動(dòng),不斷分離匯合,形成旺盛湍流。纏繞管換熱器以其結(jié)構(gòu)緊湊、耐高壓、換熱效率高和較低的投入費(fèi)用等優(yōu)點(diǎn)在核能、低溫甲醇洗工業(yè)以及LNG液化天然氣中得到廣泛應(yīng)用[1]。Neeraas等對(duì)殼側(cè)的換熱性能進(jìn)行了數(shù)百種工況的實(shí)驗(yàn)研究,使用了多種實(shí)驗(yàn)工質(zhì),研究了氣液?jiǎn)蜗鄵Q熱以及兩相混合換熱[2-3],為其他學(xué)者的相關(guān)研究提供了可靠的實(shí)驗(yàn)數(shù)據(jù)參考,并且推薦了適用于各工況的數(shù)據(jù)處理方法和經(jīng)驗(yàn)公式。Lu等以空氣為工質(zhì)實(shí)驗(yàn)測(cè)試和數(shù)值模擬了一套3層的纏繞管換熱器在恒定功率電加熱條件下的殼側(cè)流動(dòng)及傳熱特性[4-5],對(duì)比研究了定壁溫、定熱流邊界條件以及考慮管內(nèi)流動(dòng)的傳熱耦合邊界條件對(duì)計(jì)算的影響,在多數(shù)情況下用定壁溫邊界條件進(jìn)行計(jì)算,得出了努塞爾數(shù)Nu及摩擦因子f關(guān)于雷諾數(shù)Re的經(jīng)驗(yàn)關(guān)聯(lián)式。Zeng等對(duì)纏繞管換熱器進(jìn)行了數(shù)值建模,分析了換熱器的各幾何參數(shù)對(duì)流動(dòng)及換熱的影響[6],指出了RNGk-ε湍流模型在計(jì)算螺旋繞流時(shí)有更高的計(jì)算精度,并根據(jù)所得模擬結(jié)果對(duì)各幾何參數(shù)變化的影響大小進(jìn)行了評(píng)估,得出了Nu、f關(guān)于若干影響較大參數(shù)的經(jīng)驗(yàn)關(guān)聯(lián)式。文獻(xiàn)[7-9]利用對(duì)稱和周期性邊界條件對(duì)纏繞管換熱器數(shù)值模型進(jìn)行了簡(jiǎn)化,計(jì)算得出了殼側(cè)流體通道內(nèi)制冷劑的壓力、速度及溫度分布,研究了纏繞管換熱器的各幾何參數(shù)對(duì)殼程換熱性能的影響,并提出了改善換熱性能的建議。但是,目前多數(shù)研究中均為簡(jiǎn)化模型而忽略了管層間墊條的影響,為了使數(shù)值計(jì)算更接近實(shí)際情況,本文重點(diǎn)研究在相同工況下墊條、墊條排布方式、墊條數(shù)量對(duì)纏繞管換熱器殼側(cè)流動(dòng)和傳熱性能的影響。

    1 纏繞管換熱器的物理模型及數(shù)值模擬方法

    1.1 物理模型的建立

    纏繞管換熱器計(jì)算模型軸測(cè)如圖1所示,建立了有墊條和無(wú)墊條的物理模型,其中換熱管層數(shù)為兩層,管層間的墊條采用工業(yè)上較為常用的圓柱形金屬絲,換熱器的軸向剖面視圖如圖2所示。本文采用的流體通道的當(dāng)量直徑為4倍殼側(cè)流通體積與傳熱表面積之比[10]

    (1)

    (a)無(wú)墊條 (b)有墊條圖1 纏繞管換熱器計(jì)算模型軸測(cè)圖

    圖2 纏繞管換熱器的剖面視圖

    1.2 基本方程和數(shù)值方法

    經(jīng)過重整化群處理的RNGk-ε湍流模型能更好處理高應(yīng)變率及流線彎曲程度較大的流動(dòng),從而提高螺旋流的預(yù)測(cè)精度[6]。本文采用RNGk-ε湍流模型,基本方程包括質(zhì)量、動(dòng)量、能量方程,表達(dá)式如下

    連續(xù)性方程

    ui)=0

    (2)

    動(dòng)量方程

    (3)

    能量方程

    (4)

    湍動(dòng)能方程

    ε

    (5)

    湍動(dòng)能耗散方程

    (6)

    模型計(jì)算按照LNG液化時(shí)的冷卻工況設(shè)置,此工況下殼側(cè)的換熱為液態(tài)乙烷單相換熱,假定其物性恒定;流體入口為速度進(jìn)口,乙烷進(jìn)口溫度為143 K,出口為壓力出口邊界條件;管壁采用定壁溫邊界條件[4],溫度設(shè)為220 K;殼體和墊條均按無(wú)滑移絕熱壁面處理。壁面函數(shù)為默認(rèn)的標(biāo)準(zhǔn)壁面函數(shù),壓力-速度耦合采用SIMPLE算法,動(dòng)量和能量方程均為二階迎風(fēng)格式,各項(xiàng)殘差均設(shè)為1×10-6。

    模型結(jié)構(gòu)由SolidWorks建立,網(wǎng)格采用ICEM生成的非結(jié)構(gòu)化網(wǎng)格,模型進(jìn)出口均設(shè)置了延長(zhǎng)段以防止回流的影響,網(wǎng)格示意圖如圖3所示。為保證計(jì)算結(jié)果準(zhǔn)確,對(duì)無(wú)墊條模型進(jìn)行了網(wǎng)格無(wú)關(guān)性驗(yàn)證,結(jié)果如圖4所示。由圖4可知,在網(wǎng)格數(shù)大于2 365 999之后,換熱器的努塞爾數(shù)Nu的變化小于0.9%,換熱器單位長(zhǎng)度壓降ΔPl小于4.5%,模型的網(wǎng)格數(shù)為2 365 999,本文中所有模型的網(wǎng)格數(shù)為2 365 999~5 877 031。

    圖3 纏繞管換熱器網(wǎng)格示意圖

    圖4 兩層無(wú)墊條模型網(wǎng)格無(wú)關(guān)性驗(yàn)證

    1.3 數(shù)學(xué)模型的實(shí)驗(yàn)驗(yàn)證

    為了保證數(shù)值計(jì)算的準(zhǔn)確,本文對(duì)文獻(xiàn)[2]中的實(shí)驗(yàn)進(jìn)行了數(shù)值模擬,實(shí)驗(yàn)工況為甲烷單相換熱,模擬結(jié)果和實(shí)驗(yàn)數(shù)據(jù)對(duì)比如圖5所示。計(jì)算所得的Nu與實(shí)驗(yàn)值的平均誤差為9.9%,換熱管單位長(zhǎng)度壓降ΔPm的平均誤差為12.4%。除了不可避免的計(jì)算誤差之外,數(shù)值模擬與實(shí)驗(yàn)值的誤差主要來(lái)源于對(duì)模型的簡(jiǎn)化,Nu較實(shí)驗(yàn)值偏高主要是因?yàn)閿?shù)值模型將內(nèi)外殼體設(shè)置為絕熱壁面,不計(jì)漏熱損失。數(shù)值模型采用定物性計(jì)算,不考慮氣體被加熱后密度的降低所造成的ΔPm值偏高,實(shí)驗(yàn)數(shù)據(jù)測(cè)量中也會(huì)帶來(lái)一定誤差,因此數(shù)值模擬結(jié)果是正確的。

    (a)Nu

    (b)ΔPm圖5 數(shù)值模擬與實(shí)驗(yàn)值對(duì)比

    2 數(shù)值計(jì)算

    2.1 計(jì)算模型的幾何參數(shù)

    纏繞管換熱器墊條的數(shù)量需根據(jù)換熱器的尺寸來(lái)進(jìn)行選取,選取原則如表1所示[11]。為研究墊條的存在及排布方式對(duì)換熱器性能的影響,

    分別對(duì)圖

    6所示的墊條對(duì)齊及錯(cuò)開排布方式進(jìn)行建模計(jì)算,本文采用的計(jì)算模型的幾何參數(shù)如表2所示。計(jì)算模型的外層纏繞直徑均小于500 mm,每層換熱管的墊條總數(shù)取為12。

    (a)對(duì)齊排布 (b)錯(cuò)開排布圖6 兩層纏繞管模型墊條排布方式俯視圖

    參數(shù)取值換熱管外徑Dt/mm12纏繞層數(shù)N2第1、2層換熱管數(shù)4、5軸向管間距l(xiāng)/mm4墊條直徑B/mm3中心筒體直徑Dc/mm140外殼直徑Ds/mm206換熱部分高度H/mm320

    2.2 數(shù)據(jù)處理

    纏繞管換熱器的常用最小流通面積公式為[12]

    (7)

    殼側(cè)最大流速vmax可用質(zhì)量守恒求得,殼側(cè)雷諾數(shù)Re的計(jì)算式為

    (8)

    熱力計(jì)算式為

    (9)

    (10)

    (11)

    式中:N為換熱管層數(shù);B為墊條厚度;ρ為流體密度;Dh為當(dāng)量直徑;μ為動(dòng)力黏度;h為傳熱系數(shù);m為殼側(cè)流體質(zhì)量流量;cp為比定壓熱容;T為溫度;Δtm為對(duì)數(shù)平均溫差;λ為導(dǎo)熱系數(shù)。

    換熱器單位長(zhǎng)度壓降和摩擦因子為

    (12)

    (13)

    式中:ΔP為進(jìn)出口壓差;L為換熱器有效換熱段的長(zhǎng)度。

    3 計(jì)算結(jié)果分析與討論

    3.1 殼側(cè)流場(chǎng)分析

    3.1.1 墊條及排布方式對(duì)殼側(cè)軸向流動(dòng)的影響 纏繞管換熱器的殼側(cè)流動(dòng)復(fù)雜,流體在換熱管和殼體的間隙中不斷地分離匯合,湍流旺盛。流體在殼側(cè)進(jìn)行軸向流動(dòng)時(shí)主要分為:沿相鄰管層和殼體間空隙的軸向主流流動(dòng);同層管排間的低流速交叉流動(dòng)。當(dāng)存在墊條時(shí),墊條與換熱管之間沿軸向的狹窄區(qū)域易形成流體滯留區(qū)域和渦旋,使得墊條所在區(qū)域的軸向主流流速降低,并且造成壓降的提升及換熱能力的下降,墊條截面剖視圖如圖7所示。相比無(wú)墊條模型,有墊條時(shí)殼側(cè)的流動(dòng)區(qū)域變窄,會(huì)使得軸向流速有所提升。

    (a)對(duì)齊排布 (b)錯(cuò)開排布圖7 墊條截面剖視圖

    選取位于墊條附近的無(wú)墊條軸向截面進(jìn)行研究,當(dāng)進(jìn)口流速vin=1.5 m/s時(shí)截面速度和溫度分布云圖如圖8、9所示。

    (a)無(wú)墊條(b)對(duì)齊排布(c)錯(cuò)開排布圖8 vin=1.5 m/s時(shí)截面速度分布云圖

    (a)無(wú)墊條(b)對(duì)齊排布(c)錯(cuò)開排布圖9 vin=1.5 m/s時(shí)截面溫度分布云圖

    由圖8可知,相比無(wú)墊條模型,對(duì)齊和錯(cuò)開排布模型的管層間流速均有增加,其中錯(cuò)開排布模型的增幅更加明顯。墊條使殼側(cè)流體區(qū)域變小軸向流速得到提升,并且排布方式不同對(duì)軸向流速的提升效果也不同。對(duì)齊排布時(shí)墊條處低流速區(qū)域和無(wú)墊條處的高流速區(qū)域沿徑向依次對(duì)應(yīng),無(wú)加強(qiáng)流體擾動(dòng)的趨勢(shì);錯(cuò)開排布時(shí)墊條處的低流速區(qū)域?qū)?yīng)相鄰無(wú)墊條處的高流速區(qū)域,低流速區(qū)域的流體因黏性力作用流向高流速區(qū)域,增強(qiáng)了管排間交叉流動(dòng),減薄了換熱管壁的邊界層,強(qiáng)化了湍流。由圖9可知,對(duì)齊排布模型和無(wú)墊條模型的溫度分布差別不大,而錯(cuò)開排布模型的流體邊界層厚度比前兩者更薄,而且流體在流經(jīng)相同的管排數(shù)后溫度的變化幅度更大。這說明對(duì)齊排布模型軸向流速的小幅提升對(duì)殼側(cè)的換熱效果影響不大,而錯(cuò)開排布模型在軸向流速及換熱能力方面均有較大程度的提升。

    3.1.2 墊條及排布方式對(duì)殼側(cè)徑向流動(dòng)的影響 流體沿徑向截面做湍流旺盛的螺旋流動(dòng),各模型在徑向z=150 mm、vin=1.5 m/s時(shí)的截面速度分布云圖和流線圖如圖10、11所示。

    (a)無(wú)墊條 (b)對(duì)齊排布 (c)錯(cuò)開排布圖10 z=150 mm、vin=1.5 m/s時(shí)殼側(cè)徑向截面速度云圖

    (a)無(wú)墊條 (b)對(duì)齊排布 (c)錯(cuò)開排布 圖11 z=150 mm、vin=1.5 m/s時(shí)殼側(cè)局部流線分布圖

    由圖10可知,與無(wú)墊條模型相比,對(duì)齊排布模型的整體速度分布相差不大,低流速區(qū)域相對(duì)較多,而錯(cuò)開排布模型的徑向流速明顯增大。由圖11可知:無(wú)墊條模型有兩個(gè)緊挨的渦旋正在形成;在對(duì)齊排布模型中的相同區(qū)域也有渦旋形成并且在墊條附近存在另一正在耗散的渦旋;在錯(cuò)開排布模型中渦旋數(shù)量達(dá)到3個(gè),2個(gè)正在形成而1個(gè)正在耗散。這說明在墊條附近由于流體脫離壁面而更易形成渦旋,并且錯(cuò)開排布的渦旋數(shù)量多于對(duì)齊排布。沿徑向旋轉(zhuǎn)的渦旋能夠加強(qiáng)流體擾動(dòng)和強(qiáng)化換熱,當(dāng)流線匯集形成渦旋時(shí)也會(huì)使得附近區(qū)域的軸向流速進(jìn)一步提高。錯(cuò)開模型中有相對(duì)較多正在形成的渦旋,對(duì)軸向流速的提升作用更為顯著。

    由模擬計(jì)算可見:墊條對(duì)齊排布對(duì)殼側(cè)的流動(dòng)影響相對(duì)較小;錯(cuò)開排布時(shí)殼側(cè)流體軸向和徑向的流動(dòng)速度都有較為明顯的提升,并且湍流強(qiáng)度也得到較大幅度的增強(qiáng)。

    3.2 墊條及排布方式對(duì)殼側(cè)換熱性能的影響

    無(wú)墊條、墊條對(duì)齊和錯(cuò)開排布模型的努塞爾數(shù)Nu和換熱器單位長(zhǎng)度壓降ΔPl隨vin的變化如圖12、13所示。

    圖12 Nu隨vin的變化趨勢(shì)

    圖13 ΔPl隨vin的變化趨勢(shì)

    由圖12可知,對(duì)齊排布模型與無(wú)墊條模型相比Nu下降1.3%~2.5%,而錯(cuò)開排布模型的Nu在同工況下比無(wú)墊條高17.1%~18.1%,兩墊條模型中錯(cuò)開排布模型Nu比對(duì)齊排布高18.6%~20.4%。墊條與換熱管接觸的狹小區(qū)域會(huì)在軸向產(chǎn)生流動(dòng)死區(qū),帶來(lái)壓降的上升并且減弱換熱效果,但是墊條亦可在一定程度上起到加快流速和增強(qiáng)湍流的作用。錯(cuò)開排布模型的湍流強(qiáng)化作用大于死區(qū)渦旋的弱化作用,所以其Nu最高;對(duì)齊模型的強(qiáng)弱化作用基本相互抵消,Nu與無(wú)墊條模型基本持平,這與流場(chǎng)分析中的結(jié)果相符。由圖13可知,兩墊條模型的殼側(cè)換熱器單位長(zhǎng)度壓降ΔPl均比無(wú)墊條模型高,對(duì)齊模型高14.6%~17.0%,錯(cuò)開模型高38.4%~39.9%;兩墊條模型中錯(cuò)開排布模型ΔPl較高,比對(duì)齊模型高18.4%~21.7%。

    3.3 墊條數(shù)量對(duì)殼側(cè)性能的影響

    根據(jù)表1所提墊條數(shù)量選取原則,本文計(jì)算模型的每層合理墊條數(shù)為12~15,為研究墊條數(shù)量對(duì)殼側(cè)性能的影響,對(duì)墊條數(shù)量S為12、13、14、15的模型在相同進(jìn)口工況下進(jìn)行建模計(jì)算,墊條的排布方式為換熱性能較高的錯(cuò)開排布。在不同進(jìn)口流速下不同墊條數(shù)量模型的Nu、ΔPl值對(duì)比如圖14、15所示。

    圖14 不同墊條數(shù)量Nu值對(duì)比

    圖15 不同墊條數(shù)量ΔPl值對(duì)比

    由圖14、15可知,在各進(jìn)口工況下,相對(duì)S=12時(shí)的模型,當(dāng)S=13,14時(shí)Nu變化為-1.0%~0.5%,當(dāng)S增加到15時(shí)Nu下降了1.9%~3.0%。壓降方面各墊條數(shù)量模型差別不大,在較高進(jìn)口速度vin=3 m/s工況下,S=15時(shí)模型的ΔPl下降了5.3%。在相同排布情況下,墊條數(shù)量在一定范圍的改變對(duì)殼側(cè)流動(dòng)和換熱的影響均較小,但過多的墊條數(shù)對(duì)殼側(cè)的換熱有不利影響,因此在墊條數(shù)量的選取時(shí)應(yīng)滿足強(qiáng)度和安裝要求。

    3.4 強(qiáng)化傳熱性能指標(biāo)PEC評(píng)價(jià)

    PEC(performance evaluation criteria)是一種強(qiáng)化傳熱性能的評(píng)價(jià)方法[13]。強(qiáng)化換熱器中流體的流動(dòng)及傳熱方法很多,每種方法在強(qiáng)化對(duì)流換熱的同時(shí)流動(dòng)阻力等功耗也會(huì)相應(yīng)增加,這是強(qiáng)化傳熱中的不利因素。因此,需要提出一個(gè)能涵蓋強(qiáng)化傳熱和阻力增加這兩種因素的評(píng)價(jià)指標(biāo),PEC是其中一種最常用的方法,定義為

    (14)

    式中:Nu、f分別為無(wú)墊條時(shí)殼側(cè)通道內(nèi)的傳熱和阻力參數(shù)。若FP>1.0則表明傳熱強(qiáng)化效果大于流阻增加,當(dāng)墊條數(shù)為12時(shí),兩種墊條排布方式對(duì)應(yīng)的FP如圖16所示。在各進(jìn)口流速下,墊條對(duì)齊排布模型的FP始終小于1.0,即傳熱強(qiáng)化效果小于流阻的增加;錯(cuò)開排布模型則相反,FP始終大于1.0,表明相比于無(wú)墊條模型,錯(cuò)開排布模型的換熱性能更加優(yōu)越。

    圖16 兩種墊條排布方式對(duì)應(yīng)的FP

    4 結(jié) 論

    針對(duì)纏繞管換熱器,對(duì)墊條、墊條排布方式及數(shù)量對(duì)殼側(cè)流動(dòng)換熱性能的影響進(jìn)行了數(shù)值模擬研究,得出以下結(jié)論。

    (1)墊條的存在對(duì)殼側(cè)的流動(dòng)和換熱有顯著的影響,不考慮墊條的換熱器簡(jiǎn)化模型會(huì)帶來(lái)較大的計(jì)算誤差。尤其是墊條錯(cuò)開排布時(shí),軸向高流速區(qū)域?qū)?yīng)鄰層墊條處的低流速區(qū)域,對(duì)流體具有一定擾流作用,且在徑向截面上,因流體在墊條附近壁面脫離的影響,更易形成渦旋而加強(qiáng)流體擾動(dòng)。

    (2)在相同工況下,墊條錯(cuò)開排布使得換熱器的壓降增大,但換熱性能更高。錯(cuò)開排布模型的Nu比無(wú)墊條高17.1%~18.1%,ΔPl高38.4%~39.9%。而且,錯(cuò)開排布的強(qiáng)化傳熱性能指標(biāo)FP較高,始終大于1.0,說明其傳熱強(qiáng)化效果大于流阻增加,有利于殼側(cè)傳熱性能提升。

    (3)墊條數(shù)量的小幅度變化對(duì)殼側(cè)的流動(dòng)和換熱性能影響不大,但當(dāng)墊條數(shù)量過多時(shí)對(duì)換熱有不利影響。在計(jì)算工況下,當(dāng)墊條數(shù)量由12增至15時(shí),Nu下降1.9%~3.0%,因此在墊條的數(shù)量選擇時(shí)滿足強(qiáng)度和安裝要求即可,數(shù)量不宜過多。

    [1] 李亞飛, 毛文睿, 張燦燦. 纏繞管換熱器的發(fā)展與應(yīng)用 [J]. 廣東化工, 2014, 46(1): 131-132. LI Yafei, MAO Wenrui, ZHANG Cancan. Development and applications of helically coiled heat exchanger [J]. Guangdong Chemical Industry, 2014, 46(1): 131-132.

    [2] NEERAAS B O, FREDHEIM A O. Experimental shell-side heat transfer and pressure drop in gas flow for spiral-wound LNG heat exchanger [J]. International Journal of Heat and Mass Transfer, 2004, 47(2): 353-361.

    [3] NEERAAS B O, FREDHEIM A O. Experimental data and model for heat transfer, in liquid falling film flow on shell-side, for spiral-wound LNG heat exchanger [J]. International Journal of Heat and Mass Transfer, 2004, 47(14): 3565-3572.

    [4] LU Xing, DU Xueping, ZENG Min. Shell-side thermal-hydraulic performances of multilayer spiral-wound heat exchangers under different wall thermal boundary conditions [J]. Applied Thermal Engineering, 2014, 70(2): 1216-1227.

    [5] LU Xing. DU Xueping, ZHANG Sen. Experimental and numerical investigation on shell-side performance of multilayer spiral-wound heat exchangers [J]. Chemical Engineering Transactions, 2013, 35: 445-450.

    [6] ZENG Min, ZHANG Gaopeng, LI Ye. Geometrical parametric analysis of flow and heat transfer in the shell side of a spiral-wound heat exchanger [J]. Heat Transfer Engineering, 2014, 36(9): 790-805.

    [7] WU Zhiyong, CHEN Jie, PU Hui. Numerical simulation of superheated flow of refrigerant at shell side of spiral wound heat exchanger [J]. Gas & Heat, 2014, 34(8): 6-11.

    [8] WU Zhiyong, CAI Weihua, QIU Guodong. Prediction of mass transfer time relaxation parameter for boiling simulation on the shell-side of LNG spiral wound heat exchanger [J]. Advances in Mechanical Engineering, 2014(4): 1-11.

    [9] 魏江濤, 曾桃. 纏繞管換熱器殼程傳熱特性的數(shù)值研究 [J]. 化學(xué)工程與裝備, 2015(2): 21-23. WEI Jiangtao, ZENG Tao. Numerical investigation of spiral wound heat exchanger shell side heat transfer [J]. Chemical Engineer & Equipment, 2015(2): 21-23.

    [10]于清野. 纏繞管式換熱器計(jì)算方法研究 [D]. 大連: 大連理工大學(xué), 2011: 19-20.

    [11]王撫華. 化工實(shí)用專題設(shè)計(jì)手冊(cè): 上冊(cè) [M]. 北京: 學(xué)苑出版社, 2002: 183-187

    [12]FREDHEIM A O. Thermal design of coil-wound LNG heat exchangers, shell-side heat transfer and pressure drop [D]. Trondheim, Norway: Norwegian Institute of Technology, 1994: 39-42.

    [13]馬飛. 纏繞管換熱器傳熱數(shù)值模擬 [D]. 鄭州: 鄭州大學(xué), 2014: 18-19.

    [本刊相關(guān)文獻(xiàn)鏈接]

    王哲,厲彥忠,李正宇,等.液氦/超流氦制冷系統(tǒng)負(fù)壓換熱器仿真及優(yōu)化設(shè)計(jì).2016,50(8):143.[doi:10.7652/xjtuxb 201608023]

    唐上朝,胡浩威,牛東,等.大量不凝性氣體存在時(shí)不同潤(rùn)濕性管束對(duì)流冷凝傳熱實(shí)驗(yàn)研究.2016,50(5):24-31.[doi:10.7652/xjtuxb201605004]

    吳毅,胡東帥,王明坤,等.一種新型的跨臨界CO2儲(chǔ)能系統(tǒng).2016,50(3):45-49.[doi:10.7652/xjtuxb201603007]

    王斯民,王萌萌,顧昕,等.換熱器的理論強(qiáng)化傳熱評(píng)價(jià)準(zhǔn)則研究.2016,50(1):1-7.[doi:10.7652/xjtuxb201601001]

    姚爾人,王煥然,席光.一種壓縮空氣儲(chǔ)能與內(nèi)燃機(jī)技術(shù)耦合的冷熱電聯(lián)產(chǎn)系統(tǒng).2016,50(1):22-27.[doi:10.7652/xjtuxb201601004]

    粘權(quán)鑫,郭少龍,方文振,等.液化天然氣浸沒燃燒式氣化器數(shù)值模擬方法研究.2016,50(1):67-71.[doi:10.7652/xjtuxb201601011]

    楊輝著,文鍵,童欣,等.板翅式換熱器鋸齒型翅片參數(shù)的遺傳算法優(yōu)化研究.2015,49(12):90-96.[doi:10.7652/xjtuxb 201512015]

    王斯民,王萌萌,顧昕,等.螺旋折流板換熱器結(jié)構(gòu)參數(shù)多目標(biāo)優(yōu)化的數(shù)值模擬.2015,49(11):14-19.[doi:10.7652/xjtuxb201511003]

    季家東,葛培琪,畢文波.換熱器內(nèi)彈性管束流體組合誘導(dǎo)振動(dòng)響應(yīng)的數(shù)值分析.2015,49(9):24-29.[doi:10.7652/xjtuxb201509005]

    (編輯 趙煒 杜秀杰)

    Effects of Space Bars on the Shell-Side Performance of Spiral-Wound Heat Exchanger

    WEN Jian1,TONG Xin1,YANG Huizhu1,WANG Simin2,LI Yanzhong1

    (1. School of Energy and Power Engineering, Xi’an Jiaotong University, Xi’an 710049, China; 2. School of Chemical Engineering and Technology, Xi’an Jiaotong University, Xi’an 710049, China)

    So far in the simulation research of spiral-wound heat exchanger the space bars are rarely considered, the deviation between simulation results and actual situation is thus enlarged. This paper mainly studies the effects of the number and arrangement of space bars on the shell side’s thermal and hydraulic performance with numerical simulation method. Numerical results showed that the space bars affect the axial flow rate of shell side and cause more vortexes in radial sections. Compared with the calculation results of the models without space bars under the same inlet condition the Nusselt numberNuof lined up-space bar models was decreased by 1.6%-2.5%; the pressure drop per heat exchanger length ΔPlwas increased by 14.5%-17.0%; the performance evaluation criteriaFPwas decreased by 5.7%-5.9%; theNuof staggered-space bar models was increased by 17.0%-18.1%; ΔPlwas increased by 38.4%-39.9% andFPwas increased by 13.7%-15.4%. Moreover, the staggered arrangement of space bars can further enhance the turbulence of shell side, compared with the lined up arrangement theNuof staggered arrangement was increased and ΔPlwas increased accordingly in the same working conditions. When the number of space bars was more than 12 (13, 14 and 15), the number of space bars had little effect on the shell side performance. This paper provides a new approach for spiral-wound heat exchanger simulation and offers a theoretical reference for the performance optimization.

    spiral-wound heat exchanger; space bar; turbulence enhancement; pressure drop

    2016-05-17。 作者簡(jiǎn)介:文鍵(1976—),女,博士,副教授;王斯民(通信作者),男,副教授。 基金項(xiàng)目:國(guó)家自然科學(xué)基金資助項(xiàng)目(51676146);中央高校基本科研業(yè)務(wù)費(fèi)專項(xiàng)資金資助項(xiàng)目。

    時(shí)間:2016-09-14

    10.7652/xjtuxb201611003

    TK124

    A

    0253-987X(2016)11-0014-07

    網(wǎng)絡(luò)出版地址:http:∥www.cnki.net/kcms/detail/61.1069.T.20160914.1807.010.html

    猜你喜歡
    纏繞管換熱器軸向
    大型立式單級(jí)引黃離心泵軸向力平衡的研究
    帶凹槽纏繞管強(qiáng)化傳熱數(shù)值計(jì)算研究
    ASM-600油站換熱器的國(guó)產(chǎn)化改進(jìn)
    能源工程(2021年3期)2021-08-05 07:26:14
    基于停留時(shí)間分布的纏繞管內(nèi)二次流研究
    翅片管式換熱器的傳熱研究進(jìn)展
    制冷(2019年2期)2019-12-09 08:10:36
    荒銑加工軸向切深識(shí)別方法
    一種纏繞管硬度強(qiáng)化設(shè)備
    大型纏繞管式換熱器關(guān)鍵技術(shù)發(fā)展分析
    河南科技(2016年19期)2016-03-13 19:41:44
    310S全焊式板殼型換熱器的焊接
    焊接(2016年10期)2016-02-27 13:05:35
    微小型薄底零件的軸向車銑實(shí)驗(yàn)研究
    国产精品,欧美在线| 国产午夜福利久久久久久| 欧美日韩精品网址| 国产毛片a区久久久久| 在线观看舔阴道视频| 亚洲一区二区三区不卡视频| 欧美日韩国产亚洲二区| 国产精品一区二区三区四区免费观看 | 午夜视频精品福利| 午夜免费观看网址| 偷拍熟女少妇极品色| 亚洲精品久久国产高清桃花| 国产成人av教育| 亚洲国产精品久久男人天堂| 免费看光身美女| 麻豆久久精品国产亚洲av| 精华霜和精华液先用哪个| 综合色av麻豆| 很黄的视频免费| 亚洲在线观看片| 男女床上黄色一级片免费看| 美女cb高潮喷水在线观看 | 久久久久久久久久黄片| 久久精品91无色码中文字幕| 国产精品久久电影中文字幕| 国产免费男女视频| 国产乱人视频| 99精品欧美一区二区三区四区| 亚洲黑人精品在线| 美女扒开内裤让男人捅视频| 精品福利观看| 久久久成人免费电影| 国产黄色小视频在线观看| 熟女人妻精品中文字幕| 国产成人精品久久二区二区免费| 欧美又色又爽又黄视频| 黄色视频,在线免费观看| 国产淫片久久久久久久久 | 十八禁人妻一区二区| 国产精品久久电影中文字幕| 久久国产精品影院| 国产欧美日韩一区二区三| 免费大片18禁| 性欧美人与动物交配| 看黄色毛片网站| 黄色 视频免费看| 黄片小视频在线播放| 亚洲国产欧美一区二区综合| 中文在线观看免费www的网站| 手机成人av网站| 久久国产乱子伦精品免费另类| 亚洲精品国产精品久久久不卡| 免费观看的影片在线观看| 欧美中文综合在线视频| 日韩精品中文字幕看吧| 久久久成人免费电影| 免费人成视频x8x8入口观看| 久久伊人香网站| 国产精品国产高清国产av| 久久久国产欧美日韩av| а√天堂www在线а√下载| 国产亚洲欧美98| 国产精品野战在线观看| 欧美又色又爽又黄视频| 男女之事视频高清在线观看| 欧美丝袜亚洲另类 | 我要搜黄色片| 在线看三级毛片| 真人做人爱边吃奶动态| 国产男靠女视频免费网站| 久久久久久大精品| www国产在线视频色| 亚洲国产日韩欧美精品在线观看 | 欧美日韩中文字幕国产精品一区二区三区| 精品一区二区三区四区五区乱码| 成人精品一区二区免费| 国产野战对白在线观看| 亚洲欧美一区二区三区黑人| 成年女人毛片免费观看观看9| 国产亚洲精品一区二区www| 此物有八面人人有两片| 亚洲av成人av| 午夜激情福利司机影院| 好看av亚洲va欧美ⅴa在| 婷婷精品国产亚洲av在线| 露出奶头的视频| bbb黄色大片| 真人做人爱边吃奶动态| 亚洲国产精品999在线| 叶爱在线成人免费视频播放| 性欧美人与动物交配| 久久精品亚洲精品国产色婷小说| 国产精品香港三级国产av潘金莲| 成人无遮挡网站| 久久久久精品国产欧美久久久| 国产精品99久久99久久久不卡| 国产精品香港三级国产av潘金莲| 亚洲国产欧美一区二区综合| 午夜精品一区二区三区免费看| 全区人妻精品视频| 亚洲精品粉嫩美女一区| 在线看三级毛片| 91九色精品人成在线观看| 99国产极品粉嫩在线观看| 中文字幕熟女人妻在线| 亚洲第一欧美日韩一区二区三区| 一级黄色大片毛片| 精品乱码久久久久久99久播| 18禁黄网站禁片午夜丰满| 国产成人精品久久二区二区免费| 国产一区在线观看成人免费| 18禁裸乳无遮挡免费网站照片| 久久午夜综合久久蜜桃| 久久久久久人人人人人| 黄色视频,在线免费观看| 首页视频小说图片口味搜索| 国产乱人伦免费视频| 久久99热这里只有精品18| 天天躁日日操中文字幕| 级片在线观看| 91字幕亚洲| 一个人免费在线观看的高清视频| 少妇的丰满在线观看| 岛国在线观看网站| 热99在线观看视频| 禁无遮挡网站| 99久久99久久久精品蜜桃| 91av网一区二区| 亚洲黑人精品在线| 人妻丰满熟妇av一区二区三区| 老熟妇仑乱视频hdxx| 国产成人福利小说| 国产精品亚洲av一区麻豆| 亚洲va日本ⅴa欧美va伊人久久| 亚洲成人中文字幕在线播放| 黄片大片在线免费观看| 欧美高清成人免费视频www| 日本与韩国留学比较| 亚洲七黄色美女视频| 成年女人毛片免费观看观看9| 国产aⅴ精品一区二区三区波| 老司机在亚洲福利影院| 亚洲成人中文字幕在线播放| 一个人免费在线观看的高清视频| 亚洲国产精品久久男人天堂| 欧美一区二区精品小视频在线| 亚洲中文av在线| 国产麻豆成人av免费视频| 9191精品国产免费久久| 美女免费视频网站| 欧美乱妇无乱码| 免费搜索国产男女视频| 欧美一级a爱片免费观看看| 18禁观看日本| 狠狠狠狠99中文字幕| 国产久久久一区二区三区| 成人特级黄色片久久久久久久| 88av欧美| 国产亚洲欧美98| 激情在线观看视频在线高清| 国产欧美日韩精品亚洲av| 国产精品美女特级片免费视频播放器 | 国产精品久久久av美女十八| 免费观看的影片在线观看| 99久久精品一区二区三区| 人人妻,人人澡人人爽秒播| 观看美女的网站| 成人18禁在线播放| 好男人电影高清在线观看| 国内毛片毛片毛片毛片毛片| 天堂影院成人在线观看| 国产精品综合久久久久久久免费| 欧美zozozo另类| 久久午夜综合久久蜜桃| 最好的美女福利视频网| 这个男人来自地球电影免费观看| 真人做人爱边吃奶动态| 黄色女人牲交| 国产成人一区二区三区免费视频网站| 天天一区二区日本电影三级| 男女床上黄色一级片免费看| 亚洲第一欧美日韩一区二区三区| 精品国产超薄肉色丝袜足j| 麻豆国产97在线/欧美| 97超级碰碰碰精品色视频在线观看| 亚洲欧美日韩东京热| 他把我摸到了高潮在线观看| 日本撒尿小便嘘嘘汇集6| 成人国产综合亚洲| 国产三级黄色录像| 亚洲国产精品久久男人天堂| 久久欧美精品欧美久久欧美| 嫩草影院入口| 国产伦精品一区二区三区四那| 村上凉子中文字幕在线| 日本一二三区视频观看| 亚洲成av人片在线播放无| 草草在线视频免费看| 最近在线观看免费完整版| 99久久精品一区二区三区| 亚洲欧美日韩无卡精品| 两性夫妻黄色片| 黑人操中国人逼视频| 一区福利在线观看| 午夜激情福利司机影院| 精品不卡国产一区二区三区| 舔av片在线| 亚洲欧洲精品一区二区精品久久久| 91字幕亚洲| 一个人观看的视频www高清免费观看 | 国产精品电影一区二区三区| netflix在线观看网站| 欧美色视频一区免费| 亚洲成a人片在线一区二区| 中文字幕av在线有码专区| 老熟妇乱子伦视频在线观看| 热99在线观看视频| 欧美zozozo另类| 桃色一区二区三区在线观看| 天堂网av新在线| 熟女电影av网| 午夜免费观看网址| 麻豆av在线久日| 午夜福利成人在线免费观看| 国产午夜精品论理片| 国产野战对白在线观看| 嫩草影院精品99| 精品一区二区三区四区五区乱码| 精品国内亚洲2022精品成人| 老汉色∧v一级毛片| 欧美黄色片欧美黄色片| 午夜免费激情av| 最近最新中文字幕大全电影3| 国产亚洲av高清不卡| 一进一出好大好爽视频| 嫩草影院入口| 丰满人妻熟妇乱又伦精品不卡| 亚洲av成人不卡在线观看播放网| 两性夫妻黄色片| 久久久久国产精品人妻aⅴ院| 麻豆国产97在线/欧美| 亚洲国产精品久久男人天堂| 亚洲自拍偷在线| 青草久久国产| 国语自产精品视频在线第100页| 久久久国产精品麻豆| 日韩欧美国产在线观看| 亚洲成人久久性| 九九热线精品视视频播放| 国产亚洲精品久久久com| 日本一本二区三区精品| 国产极品精品免费视频能看的| 国产高清激情床上av| 91久久精品国产一区二区成人 | 一夜夜www| 亚洲第一欧美日韩一区二区三区| 欧美乱色亚洲激情| 亚洲成人精品中文字幕电影| 日本免费a在线| 国产v大片淫在线免费观看| 看片在线看免费视频| 999久久久国产精品视频| 久久这里只有精品中国| 国产成人精品久久二区二区91| aaaaa片日本免费| 欧美日韩综合久久久久久 | 免费av不卡在线播放| 亚洲性夜色夜夜综合| 一进一出抽搐动态| 两个人的视频大全免费| 成年女人永久免费观看视频| 久久午夜综合久久蜜桃| 美女午夜性视频免费| 老鸭窝网址在线观看| 国产成人精品久久二区二区免费| 热99在线观看视频| 一卡2卡三卡四卡精品乱码亚洲| 国内精品一区二区在线观看| 我要搜黄色片| 亚洲午夜理论影院| 国产精品香港三级国产av潘金莲| 久久这里只有精品19| 美女被艹到高潮喷水动态| 男女午夜视频在线观看| 日本撒尿小便嘘嘘汇集6| 欧美激情久久久久久爽电影| 亚洲欧美一区二区三区黑人| 亚洲在线自拍视频| 日韩欧美三级三区| 法律面前人人平等表现在哪些方面| 一区二区三区激情视频| 国产熟女xx| 亚洲熟妇中文字幕五十中出| 伦理电影免费视频| 神马国产精品三级电影在线观看| 黄色片一级片一级黄色片| 动漫黄色视频在线观看| 久久久久久久久久黄片| 国产精品 欧美亚洲| 2021天堂中文幕一二区在线观| 免费人成视频x8x8入口观看| xxxwww97欧美| 午夜两性在线视频| 亚洲av熟女| 男女那种视频在线观看| 日本 av在线| 免费在线观看亚洲国产| 成人特级黄色片久久久久久久| 69av精品久久久久久| 黄色丝袜av网址大全| 99re在线观看精品视频| 午夜精品久久久久久毛片777| 亚洲av五月六月丁香网| 视频区欧美日本亚洲| 亚洲在线观看片| 在线观看免费视频日本深夜| 亚洲aⅴ乱码一区二区在线播放| 久久久久性生活片| 特级一级黄色大片| 99热只有精品国产| 成年女人毛片免费观看观看9| 国产蜜桃级精品一区二区三区| 久久草成人影院| 精品电影一区二区在线| 久久久久久久精品吃奶| 免费无遮挡裸体视频| 男女床上黄色一级片免费看| 午夜精品一区二区三区免费看| 欧美+亚洲+日韩+国产| 99精品久久久久人妻精品| 成人精品一区二区免费| 热99在线观看视频| 校园春色视频在线观看| 人妻夜夜爽99麻豆av| 美女cb高潮喷水在线观看 | 免费大片18禁| 欧美成人一区二区免费高清观看 | 久久久国产欧美日韩av| 听说在线观看完整版免费高清| xxxwww97欧美| 波多野结衣巨乳人妻| 亚洲成人久久性| 国产精品精品国产色婷婷| 久久午夜亚洲精品久久| 男人和女人高潮做爰伦理| 久久久久精品国产欧美久久久| 亚洲欧美精品综合一区二区三区| 国产私拍福利视频在线观看| 精品免费久久久久久久清纯| 成年女人毛片免费观看观看9| 午夜福利在线观看免费完整高清在 | 俄罗斯特黄特色一大片| 欧美高清成人免费视频www| 亚洲av成人av| 级片在线观看| 美女黄网站色视频| 精品久久久久久成人av| 亚洲美女视频黄频| 成人三级做爰电影| 男人舔奶头视频| 在线免费观看不下载黄p国产 | 欧美一区二区国产精品久久精品| 窝窝影院91人妻| 亚洲成人精品中文字幕电影| 欧美一级毛片孕妇| 亚洲中文字幕一区二区三区有码在线看 | 中出人妻视频一区二区| 午夜福利18| 夜夜夜夜夜久久久久| 全区人妻精品视频| 好男人在线观看高清免费视频| 国产精品亚洲美女久久久| 亚洲九九香蕉| 国内精品久久久久精免费| 一二三四社区在线视频社区8| 国产 一区 欧美 日韩| 少妇人妻一区二区三区视频| 狂野欧美激情性xxxx| 精品国产亚洲在线| 日韩精品青青久久久久久| 99久久无色码亚洲精品果冻| 国产精品野战在线观看| 国产伦精品一区二区三区视频9 | 久久久久久国产a免费观看| 制服人妻中文乱码| 五月伊人婷婷丁香| 蜜桃久久精品国产亚洲av| 亚洲中文字幕一区二区三区有码在线看 | 久久热在线av| 国产真人三级小视频在线观看| 欧美极品一区二区三区四区| 老熟妇仑乱视频hdxx| 国产成人av教育| 国产精品1区2区在线观看.| 亚洲七黄色美女视频| 日韩精品中文字幕看吧| 极品教师在线免费播放| 99国产精品一区二区三区| 日韩三级视频一区二区三区| 窝窝影院91人妻| 欧美日韩中文字幕国产精品一区二区三区| 天堂影院成人在线观看| 老汉色av国产亚洲站长工具| 午夜福利成人在线免费观看| 久久精品国产清高在天天线| 天堂网av新在线| 搡老岳熟女国产| 亚洲av电影不卡..在线观看| 日本 av在线| 非洲黑人性xxxx精品又粗又长| 男人舔女人的私密视频| 高清毛片免费观看视频网站| 老汉色∧v一级毛片| 最近在线观看免费完整版| 三级男女做爰猛烈吃奶摸视频| 一本久久中文字幕| 国产伦一二天堂av在线观看| 国产成人精品久久二区二区91| 婷婷精品国产亚洲av| 国产精品爽爽va在线观看网站| 欧美日韩国产亚洲二区| 午夜免费成人在线视频| 欧美日韩亚洲国产一区二区在线观看| 久久精品91蜜桃| 亚洲色图 男人天堂 中文字幕| 免费无遮挡裸体视频| 岛国视频午夜一区免费看| 级片在线观看| 男人舔女人下体高潮全视频| 亚洲五月婷婷丁香| 国产亚洲欧美98| 一a级毛片在线观看| 操出白浆在线播放| 国产成+人综合+亚洲专区| 亚洲国产精品sss在线观看| 九色成人免费人妻av| 亚洲狠狠婷婷综合久久图片| 精品一区二区三区四区五区乱码| 1000部很黄的大片| 成年免费大片在线观看| 一进一出抽搐gif免费好疼| 精品欧美国产一区二区三| 老熟妇仑乱视频hdxx| a级毛片在线看网站| 日韩国内少妇激情av| 免费在线观看影片大全网站| 久99久视频精品免费| 美女被艹到高潮喷水动态| 精品久久久久久久毛片微露脸| 两人在一起打扑克的视频| 国产高潮美女av| 真人一进一出gif抽搐免费| 精品99又大又爽又粗少妇毛片 | 欧美性猛交黑人性爽| 日本一本二区三区精品| 久久久久亚洲av毛片大全| 日韩高清综合在线| 日日夜夜操网爽| 免费av毛片视频| 午夜精品久久久久久毛片777| 亚洲专区国产一区二区| 亚洲精品美女久久av网站| 日韩欧美国产在线观看| 亚洲av成人不卡在线观看播放网| 香蕉国产在线看| 亚洲av成人一区二区三| 一二三四社区在线视频社区8| 中出人妻视频一区二区| 欧美又色又爽又黄视频| 国产精品久久久久久人妻精品电影| 欧美高清成人免费视频www| 久久这里只有精品中国| 欧美一级毛片孕妇| 日本成人三级电影网站| 一个人观看的视频www高清免费观看 | 黄色视频,在线免费观看| 可以在线观看毛片的网站| 亚洲精品色激情综合| 午夜成年电影在线免费观看| 国产成人福利小说| 国产精品乱码一区二三区的特点| 色精品久久人妻99蜜桃| 在线a可以看的网站| 久99久视频精品免费| 午夜福利18| 美女高潮的动态| 国内揄拍国产精品人妻在线| 欧美极品一区二区三区四区| 长腿黑丝高跟| 天天一区二区日本电影三级| 搡老妇女老女人老熟妇| 中文亚洲av片在线观看爽| 在线免费观看不下载黄p国产 | 母亲3免费完整高清在线观看| 国产av麻豆久久久久久久| 亚洲av第一区精品v没综合| 亚洲中文字幕日韩| 国产精品久久视频播放| 伊人久久大香线蕉亚洲五| 久久久久国产一级毛片高清牌| 国产av在哪里看| 伦理电影免费视频| 国产午夜精品久久久久久| 啦啦啦韩国在线观看视频| 亚洲专区国产一区二区| 一个人看视频在线观看www免费 | 成人一区二区视频在线观看| 久久香蕉精品热| 午夜福利免费观看在线| 不卡av一区二区三区| 深夜精品福利| 免费看光身美女| 亚洲成av人片免费观看| 全区人妻精品视频| 欧美黄色淫秽网站| 国产91精品成人一区二区三区| 日本在线视频免费播放| 久久久久久大精品| 白带黄色成豆腐渣| 亚洲国产欧美一区二区综合| 日韩精品青青久久久久久| 国产成人一区二区三区免费视频网站| 九九久久精品国产亚洲av麻豆 | 国产乱人伦免费视频| 麻豆国产av国片精品| 国产精品综合久久久久久久免费| 欧美另类亚洲清纯唯美| www日本黄色视频网| 久久久久久久久免费视频了| 久久中文字幕一级| 欧美乱码精品一区二区三区| 在线观看午夜福利视频| 亚洲真实伦在线观看| 国产美女午夜福利| 久久99热这里只有精品18| 国产激情偷乱视频一区二区| 亚洲无线在线观看| 精品不卡国产一区二区三区| 岛国在线免费视频观看| 久久久久久人人人人人| 法律面前人人平等表现在哪些方面| 成年女人永久免费观看视频| 午夜日韩欧美国产| 国产精品久久视频播放| 亚洲精品一区av在线观看| 国产69精品久久久久777片 | 亚洲精品美女久久av网站| 国产成人av激情在线播放| 搡老岳熟女国产| 色综合欧美亚洲国产小说| 麻豆久久精品国产亚洲av| 免费观看的影片在线观看| 国产精品国产高清国产av| 99国产精品一区二区蜜桃av| 99久久精品一区二区三区| 九九热线精品视视频播放| 国产高清三级在线| 最新中文字幕久久久久 | 精品午夜福利视频在线观看一区| 久久久久久人人人人人| 午夜激情福利司机影院| 变态另类成人亚洲欧美熟女| 18禁裸乳无遮挡免费网站照片| 中文字幕高清在线视频| 视频区欧美日本亚洲| 法律面前人人平等表现在哪些方面| 热99在线观看视频| 国产美女午夜福利| 日韩精品青青久久久久久| 性色av乱码一区二区三区2| 色综合婷婷激情| 欧美黑人欧美精品刺激| 高潮久久久久久久久久久不卡| 久久精品人妻少妇| 最新在线观看一区二区三区| 搡老岳熟女国产| 免费在线观看亚洲国产| 久久国产乱子伦精品免费另类| 熟女人妻精品中文字幕| 日韩欧美国产一区二区入口| 毛片女人毛片| av福利片在线观看| 一个人免费在线观看的高清视频| 成人av一区二区三区在线看| 国产成人系列免费观看| 亚洲av第一区精品v没综合| 少妇熟女aⅴ在线视频| 99热只有精品国产| 一级毛片女人18水好多| 热99在线观看视频| 久久久精品大字幕| 国产精品九九99| 欧美一级a爱片免费观看看| 国产亚洲欧美在线一区二区| 久久久久久久久免费视频了| 天天添夜夜摸| www日本在线高清视频| 黄频高清免费视频| 黄色日韩在线| 偷拍熟女少妇极品色| 曰老女人黄片| 免费看美女性在线毛片视频| 一个人看的www免费观看视频| 美女cb高潮喷水在线观看 | 日本成人三级电影网站| 久9热在线精品视频| 免费看光身美女| 亚洲av熟女| 99视频精品全部免费 在线 | 欧美最黄视频在线播放免费| 日韩人妻高清精品专区| 亚洲第一欧美日韩一区二区三区| 亚洲av日韩精品久久久久久密| 亚洲 欧美 日韩 在线 免费| 少妇裸体淫交视频免费看高清|