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

    換熱管表面的分形表征及積灰特性數(shù)值模擬

    2024-11-27 00:00:00周友行謝寶安高騰騰楊文佳宮天宇彭麗
    中國粉體技術(shù) 2024年6期
    關(guān)鍵詞:數(shù)值模擬

    摘要:【目的】模擬換熱器管束的實(shí)際工況,研究不同粒徑飛灰顆粒在粗糙管束表面的沉積特性?!痉椒ā炕诜中卫碚?,通過改進(jìn)的Weierstrass-Mandelbrot函數(shù)建立不同粗糙程度的管束表面模型,使用Fluent軟件,結(jié)合用戶自定義函數(shù),分析表面形貌對(duì)流體流動(dòng)的影響以及不同粒徑下的顆粒沉積與碰撞特性?!窘Y(jié)果】換熱器管束表面粗糙度對(duì)壁面附近流速和湍流強(qiáng)度影響顯著,粗糙表面會(huì)使湍流強(qiáng)度增大、流體速度降低,進(jìn)一步加快顆粒沉積;流體壓降隨著管束表面粗糙度的增大而增大;相對(duì)于光滑管束表面,粗糙表面通過渦流卷吸作用增強(qiáng)對(duì)顆粒的捕獲效果,導(dǎo)致顆粒具有更高的沉積率和壁面碰撞概率。【結(jié)論】飛灰顆粒的沉積與管束粗糙表面的形成具有正反饋效應(yīng),揭示表面粗糙度與顆粒沉積特性之間的內(nèi)在關(guān)聯(lián)。

    關(guān)鍵詞:分形理論;飛灰沉積;數(shù)值模擬;氣固兩相流

    中圖分類號(hào):TB44文獻(xiàn)標(biāo)志碼:A

    引用格式:

    周友行,謝寶安,高騰騰,等. 換熱管表面的分形表征及積灰特性數(shù)值模擬[J]. 中國粉體技術(shù),2024,30(6):85-96.

    ZHOU Youhang,XIE Baoan,GAO Tengteng,et al. Fractal characterization of heat exchange tube surfaces and numerical simu?lation on ash deposition characteristics[J]. China Powder Science and Technology,2024,30(6):85?96.

    燃煤在燃燒的過程中會(huì)產(chǎn)生煙氣,煙氣中攜帶著大量的熱量,在鍋爐中使用換熱設(shè)備能夠有效提高能源利用率,節(jié)約煤炭資源,然而,煤中含有多種礦物成分,當(dāng)煤粉在爐膛中燃燒時(shí),煤中所含的礦物質(zhì)會(huì)轉(zhuǎn)化為飛灰,飛灰在煙氣的帶動(dòng)下會(huì)沉積在換熱管道表面。一方面,沉積在管道表面的灰層具有很強(qiáng)的絕熱性,幾毫米厚的灰層就能導(dǎo)致傳熱性能降低約25%;另一方面,飛灰中的硫化物會(huì)腐蝕金屬管表面,致使管壁變薄而引起爆管,降低了設(shè)備的安全性[1]。對(duì)此,國內(nèi)外學(xué)者開展了大量相關(guān)研究,用于解釋飛灰顆粒的沉積機(jī)制并預(yù)測(cè)不同條件下?lián)Q熱管表面的顆粒沉積行為。

    An等[2]設(shè)計(jì)了一種在管后設(shè)置附有小圓柱體的混合結(jié)構(gòu),能夠有效降低管束表面顆粒沉積率,但同時(shí)會(huì)降低傳熱性能。Aarnes等[3]通過模擬管束上的顆粒沖擊,發(fā)現(xiàn)隨著雷諾數(shù)增大,背風(fēng)面的顆粒撞擊率逐漸提高。彭巖等[4]通過實(shí)驗(yàn)研究發(fā)現(xiàn),流速的增加會(huì)減少管道迎風(fēng)面的積灰量,但同時(shí)會(huì)增加背風(fēng)面的積灰量。Yang等[5]結(jié)合慣性沖擊、熱泳、冷凝3種沉積機(jī)制,建立了動(dòng)態(tài)計(jì)算流體力學(xué)模型來預(yù)測(cè)煤粉燃燒過程中的飛灰沉積,發(fā)現(xiàn)慣性沖擊和熱泳對(duì)飛灰顆粒沉積起主導(dǎo)作用。Lin等[6]研究了堿金屬富集煤粉燃燒后的積灰現(xiàn)象,發(fā)現(xiàn)隨著入口煙氣溫度的升高,飛灰顆粒的沉積率也隨之增大。袁寶強(qiáng)等[7]采用計(jì)算流體力學(xué)-離散相模型和動(dòng)網(wǎng)格方法,揭示了中低溫順排管束迎風(fēng)面飛灰沉積特性,并描述了其沉積形貌。唐嬋等[8]模擬了管束微小顆粒的沉積,發(fā)現(xiàn)顆粒粒徑較小時(shí),流場(chǎng)和溫度分布對(duì)其運(yùn)動(dòng)軌跡影響顯著。

    以上對(duì)于換熱器管束積灰特性的研究大多著重于光滑表面,對(duì)影響顆粒沉積的因素考慮不夠全面,在實(shí)際應(yīng)用過程中,換熱器管道受到來自飛灰顆粒的沉積、磨損、腐蝕等影響,表面會(huì)變得越來越粗糙。換熱器管道表面形貌也是影響顆粒沉積的一個(gè)重要條件,但目前有關(guān)壁面粗糙管束的積灰特性研究比較欠缺。

    本文中基于分形理論,通過改進(jìn)的Weierstrass-Mandelbrot(W-M)分形函數(shù)模擬不同復(fù)雜程度的換熱表面,分析換熱器管道表面形貌對(duì)流體流動(dòng)以及飛灰顆粒沉積、碰撞的影響。通過深入了解飛灰顆粒在實(shí)際工況下的沉積機(jī)制,有助于預(yù)測(cè)顆粒沉積速率和沉積分布,從而制定相應(yīng)的維護(hù)策略,同時(shí)為氣固兩相流動(dòng)及顆粒沉積理論提供新的數(shù)據(jù)支撐和理論基礎(chǔ)。

    1計(jì)算模型

    1.1湍流模型

    本文中氣流速度較低,流體被假設(shè)為不可壓縮湍流。流體區(qū)域計(jì)算基于質(zhì)量、動(dòng)量和能量守恒方程。相比標(biāo)準(zhǔn)湍動(dòng)能-湍流耗散率(k-ε)模型,重整化群(renormalization group,RNG)k-ε模型考慮了渦流對(duì)湍流的影響,計(jì)算結(jié)果更加精確、可靠,因此,本研究中選用RNG k-ε模型模擬湍流流動(dòng),湍動(dòng)能k與湍流耗散率ε的方程分別如下

    式中:ρ為流體密度;t為時(shí)間;xi和xj為空間坐標(biāo);μeff為有效黏度;ui為流體的速度分量;Gk為由平均

    速度梯度引起的湍動(dòng)能k的產(chǎn)生項(xiàng);αk和αε分別為湍動(dòng)能k和湍流耗散率ε的模型系數(shù);C1(*)ε和C2ε為經(jīng)

    驗(yàn)常數(shù)。

    1.2離散相模型

    由于離散相體積分?jǐn)?shù)遠(yuǎn)小于10%,因此采用離散相模型追蹤顆粒軌跡,用歐拉-拉格朗日方法來描述流體與顆粒運(yùn)動(dòng)。在本文中,忽略顆粒間相互碰撞和顆粒對(duì)流場(chǎng)的影響,假設(shè)粒子是球型并且不可旋轉(zhuǎn)?;谂nD第二定律的顆粒運(yùn)動(dòng)方程為

    式中:vp為顆粒速度;v為流體速度;τr為顆粒的弛豫時(shí)間;g為重力加速度;ρp為顆粒密度;ρg為氣體密度;mp為顆粒質(zhì)量;Fx為作用在單位質(zhì)量顆粒上的其他作用力,包括壓力梯度力、薩夫曼升力、附加質(zhì)量力、熱泳力等。

    飛灰顆粒的主要沉積機(jī)制有慣性撞擊、熱泳沉積、冷凝沉積、化學(xué)沉積[9]。處于中低溫條件下的固化飛灰顆粒主要是因慣性撞擊而沉積。薩夫曼升力和熱泳力對(duì)飛灰顆粒的沉積有一定的影響,而巴塞特力與壓力梯度力等對(duì)顆粒沉積的作用可以忽略不計(jì)[10]。鑒于此,本文中僅考慮重力、曳力、薩夫曼升力和熱泳力。薩夫曼升力是當(dāng)流體區(qū)域存在速度梯度變化時(shí),流場(chǎng)施加在顆粒上的剪切作用力,計(jì)算公式為

    式中:FS為薩夫曼升力;dp為顆粒的等效直徑;μ為流體的動(dòng)力黏度為流速v關(guān)于垂直流動(dòng)方向的

    坐標(biāo)y的梯度。

    熱泳力Ft是由于溫度梯度引起的微小顆?;蚍肿釉诹黧w中移動(dòng)的力,可表示為

    Ft=-6πμdpf(Kn)?T,(5)

    式中:?T為氣體溫差;f(Kn)為與克努森數(shù)Kn相關(guān)的修正函數(shù),其計(jì)算公式為

    式中:λ為氣體導(dǎo)熱率;λp為顆粒導(dǎo)熱率;Ttm、Tm、Tt分別為熱速度系數(shù)、黏性滑動(dòng)系數(shù)、溫度躍變系數(shù),分別取0.461、1.19、3.32;a3為二次滑動(dòng)系數(shù),取2.4;Kn為克努森數(shù),計(jì)算公式如下:

    (7)

    式中,l為空氣分子平均自由程。

    1.3顆粒沉積模型

    當(dāng)飛灰顆粒與壁面接觸時(shí),須要選用合適的沉積模型判定顆粒狀態(tài)。Kleinhans等[11]的研究工作表明,換熱器或省煤器等在中低溫(Tlt;500℃)工作條件下,應(yīng)用基于能量守恒和臨界速度理論的沉積模型預(yù)測(cè)結(jié)果較好。該模型通過臨界沉積速度vcr判斷顆粒狀態(tài),當(dāng)顆粒對(duì)壁面的法向碰撞速度vn低于或等于臨界沉積速度vcr時(shí),顆粒將粘附在壁面,否則顆粒回彈[12]。顆粒沉積模型可表示如下:

    pst={(8)

    式中,pst為顆粒附壁概率,

    vcr=(9)

    式中:R為恢復(fù)系數(shù);E為有效剛度,計(jì)算公式如下

    E=0.51(10)

    (11)

    (12)

    式中:Ys和Yp分別為壁面和顆粒的等效楊氏模量;σ s和σp為壁面和顆粒的泊松比;Es和Ep分別為壁面材料剛度和顆粒材料剛度。

    本文中所用飛灰顆粒材料設(shè)置為K2 SO4,表面材料為結(jié)構(gòu)鋼。顆粒和換熱管壁面的具體參數(shù)如表1所示。

    1.4表面形貌模型

    分形理論是一種數(shù)學(xué)理論,用于描述和分析那些在不同尺度上表現(xiàn)出自相似特征的復(fù)雜、不規(guī)則形態(tài),通過分形維數(shù)量化其復(fù)雜度。在工程領(lǐng)域,分形理論用于分析復(fù)雜流體力學(xué)問題,如湍流、顆粒沉積、燃燒過程等,它提供了一種新的視角來理解和模擬這些復(fù)雜現(xiàn)象。Majumdar等[13]提出的W-M分形函數(shù)滿足連續(xù)性、處處不可微、自仿射性等表面特征,因此,本文中將W-M分形函數(shù)引入到換熱器表面模擬中。為了更好地模擬換熱管表面粗糙圓形輪廓,我們對(duì)方程進(jìn)行變形并且加入隨機(jī)數(shù),使其具有隨機(jī)的相位,改進(jìn)W-M方程為

    式中:Z為表面粗糙度輪廓;R0為圓管半徑;G為尺度系數(shù);D為表面分形維數(shù),1lt;Dlt;2;γ為常數(shù),通常取1.5;θ為相位偏移量;x0為θ=0處分形表面波動(dòng)的初始位置的取值;n 1為截止頻率;φ為[0,2π]范圍內(nèi)的隨機(jī)數(shù)。

    本文中表面形貌參數(shù)采用表面粗糙度,用來描述表面幾何特征,均方根粗糙度Rq公式如下

    |Zi|,(14)

    式中:N為測(cè)量點(diǎn)的總數(shù)量;i為測(cè)量點(diǎn);Zi為測(cè)量輪廓相對(duì)于平均線的高度。

    依據(jù)公式(13),使用MATLAB軟件編寫相關(guān)程序建立模型?;诂F(xiàn)有W-M分形函數(shù)法對(duì)粗糙表面模擬的研究,本文中參考文獻(xiàn)[14]將表面分形維數(shù)設(shè)置為D=1.205,圓管半徑為R0=20 mm,為了建立不同程度的粗糙表面,分別取G=1×10-6、1×10-8、1×10-10。依據(jù)公式(14),編寫粗糙度計(jì)算程序,計(jì)算得到粗糙表面的粗糙度分別為0.25、0.60、1.50 mm,同時(shí)定義光滑壁面粗糙度為0,根據(jù)粗糙度從小到大的順序?qū)⑵浞譃?種案例(案例1—案例4),通過MATLAB軟件生成坐標(biāo)點(diǎn)數(shù)據(jù),將坐標(biāo)點(diǎn)數(shù)據(jù)導(dǎo)入U(xiǎn)G軟件得到實(shí)體模型,并對(duì)邊界光順處理。

    本文中共對(duì)比4種不同粗糙度下的換熱管束表面流動(dòng)與積灰特性,在x、y平面中圓管外輪廓如圖1所示。

    2物理模型

    2.1管式換熱器的幾何結(jié)構(gòu)

    以典型管式換熱器作為研究對(duì)象,多項(xiàng)研究表明,由于管排式換熱器具有周期性和對(duì)稱性的特征,其二維模型也能得到可靠結(jié)果,并且能夠降低大量計(jì)算成本,因此,將具有5排管束的管式換熱器的二維通道模型作為計(jì)算域,邊界條件設(shè)置為周期性邊界條件。具體模型如圖2所示。

    入口區(qū)域長度為200 mm,以確保流體速度的均勻性,出口區(qū)域長度延伸至280 mm,避免氣體回流。計(jì)算域總長為840 mm。

    2.2求解方式與模擬參數(shù)

    數(shù)值模擬求解采用軟件ANSYS Fluent18.0。左側(cè)邊界設(shè)置為速度入口,右側(cè)邊界設(shè)置為壓力出口,y軸方向兩側(cè)邊界面設(shè)置為周期性邊界條件,管道表面設(shè)置為無滑移壁面條件。其他數(shù)值模擬參數(shù)見表2。

    通過Fluent軟件編寫自定義函數(shù)(user-defined function,UDF)建立顆粒沉積模型,當(dāng)顆粒碰撞管壁時(shí),用顆粒沉積模型來判定粒子狀態(tài),方程組計(jì)算求解采用耦合式求解器。控制方程的對(duì)流項(xiàng)采用二階迎風(fēng)格式,擴(kuò)散項(xiàng)采用中心差分格式。設(shè)定當(dāng)各個(gè)方程殘差小于1×10-4時(shí),模型收斂。本文中每個(gè)案例釋放2萬個(gè)顆粒,模擬顆粒粒徑為1~20 μm,粒徑分布擬合為羅辛-拉姆勒分布,使用隨機(jī)游走模型來描述顆粒在湍流中的擴(kuò)散。

    2.3網(wǎng)格劃分與無關(guān)性驗(yàn)證

    采用前處理軟件ICEM對(duì)計(jì)算域進(jìn)行離散化,由于壁面結(jié)構(gòu)比較復(fù)雜,因此使用非結(jié)構(gòu)化單元進(jìn)行網(wǎng)格劃分,并且對(duì)近壁面網(wǎng)格進(jìn)行了局部加密,如圖3所示。

    網(wǎng)格數(shù)的增加有助于提高計(jì)算結(jié)果的精度,然而,隨著網(wǎng)格數(shù)量的增加,計(jì)算規(guī)模和時(shí)間成本也會(huì)顯著增加,每個(gè)案例共劃分4種不同數(shù)量的網(wǎng)格進(jìn)行網(wǎng)格獨(dú)立性驗(yàn)證。以案例4為例,采用4種網(wǎng)格方案(網(wǎng)格a、b、c、d)進(jìn)行檢驗(yàn),網(wǎng)格數(shù)量分別為78563、126524、264839、397433。在相同初始條件下,以表面換熱系數(shù)h和壓降Δp作為監(jiān)測(cè)參數(shù),驗(yàn)證網(wǎng)格獨(dú)立性,具體數(shù)值見表3,從中可以看出網(wǎng)格c與網(wǎng)格d的計(jì)算結(jié)果最大偏差小于1%,綜合考慮計(jì)算時(shí)間和計(jì)算精度的權(quán)衡,最終確定案例4的網(wǎng)格數(shù)量為264839。同理,得到案例1、案例2和案例3劃分的網(wǎng)格數(shù)量分別為118355、234505、208827。

    3仿真結(jié)果及分析

    3.1表面形貌對(duì)氣體流動(dòng)的影響

    在氣固兩相流中,流體狀態(tài)會(huì)對(duì)顆粒運(yùn)動(dòng)產(chǎn)生至關(guān)重要的影響,因此,通過對(duì)流體流動(dòng)進(jìn)行分析,能夠更好地解釋顆粒在模型中的運(yùn)動(dòng)特性。研究表明,復(fù)雜的表面形貌將會(huì)對(duì)其附近的流場(chǎng)產(chǎn)生擾動(dòng)[15]。下面將從壓力、速度、流動(dòng)結(jié)構(gòu)等角度探究不同粗糙程度對(duì)流場(chǎng)的影響。

    當(dāng)入口處流體速度為5 m/s時(shí),對(duì)數(shù)值模擬結(jié)果進(jìn)行后處理分析,不同管壁表面粗糙度的流場(chǎng)壓力分布如圖4所示。由圖可見,隨著表面粗糙度的增加,換熱管的入口壓力逐漸增大,這是由于粗糙度增大導(dǎo)致流體摩擦增大,煙氣在流動(dòng)過程中會(huì)損失更多的能量。

    不同管壁表面粗糙度的壓降變化如圖5所示。由圖可見,壓降隨著管壁表面粗糙度增加而增大。這是因?yàn)?,管壁表面粗糙度的增加?huì)導(dǎo)致流體流動(dòng)所要克服的內(nèi)摩擦力增大,從而引起壓力損失。較大的壓降會(huì)導(dǎo)致流體流速減緩,顆粒動(dòng)能減小,特別是在近壁面位置,較低的流速會(huì)導(dǎo)致顆粒更容易沉積。

    圖6展示了對(duì)于不同粗糙程度的表面,氣流在x=0處近壁面速度分布隨壁面距離的變化。在Rq=0 mm的光滑表面上,流體速度迅速達(dá)到最大值并趨于穩(wěn)定。隨著粗糙度的增加,速度分布曲線出現(xiàn)明顯變化。在Rq=1.50 mm的粗糙表面上,流體速度在靠近壁面的區(qū)域顯著降低,這意味著顆粒更容易沉積在壁面,從而具有更高的沉積幾率。

    圖7展示了對(duì)于不同粗糙程度的表面,氣流在x=0處近壁面湍動(dòng)能分布隨壁面距離的變化。由圖可見,粗糙度對(duì)湍動(dòng)能的影響只發(fā)生在近壁區(qū)域,遠(yuǎn)離壁面的區(qū)域幾乎不受影響,這與Marchis等[16]的研究結(jié)果一致。隨著粗糙度的增加,湍動(dòng)能的峰值逐漸上升,較大的湍動(dòng)能會(huì)引發(fā)顆粒的劇烈擾動(dòng)。對(duì)于Rq=1.50 mm的粗糙表面,湍動(dòng)能峰值顯著高于其他表面,表明粗糙表面增強(qiáng)了湍流強(qiáng)度和能量的耗散,這意味著較大的湍流強(qiáng)度會(huì)增加顆粒在近壁面區(qū)域的碰撞次數(shù),進(jìn)而增加顆粒沉積率。

    圖8展示了粗糙度為0的光滑表面與粗糙度為1.5 mm的粗糙表面在xy平面內(nèi)的近壁面區(qū)域流線形態(tài)。由圖可見,光滑表面相對(duì)平坦,對(duì)流體的干擾較小,而粗糙表面則因突起和凹陷的存在,會(huì)顯著影響流體流動(dòng),導(dǎo)致流體在接近表面時(shí)產(chǎn)生不規(guī)則的變化和擾動(dòng)。當(dāng)粗糙表面上的波峰之間的距離足夠大時(shí),流體在波峰的迎風(fēng)面和背風(fēng)面之間形成渦旋結(jié)構(gòu),這是由迎風(fēng)面的壓力變化和背風(fēng)面的流體回流造成的。渦旋結(jié)構(gòu)不僅使流體局部速度和壓力發(fā)生變化,而且可能對(duì)傳熱和質(zhì)量傳遞過程產(chǎn)生重要影響。粗糙表面通過渦流卷吸作用會(huì)增強(qiáng)對(duì)顆粒的捕獲效果,導(dǎo)致顆粒更容易沉積在輪廓表面。

    3.2表面形貌對(duì)顆粒沉積率的影響

    壁面捕獲效率是指在氣固兩相流中,顆粒被壁面捕獲的效率或比例。通常用捕獲到的顆粒數(shù)與進(jìn)入近壁區(qū)域的顆??倲?shù)之比來表示,反映了壁面對(duì)顆粒沉積的影響,高壁面捕獲效率意味著更多的顆粒在接觸壁面后被沉積。在本文中,壁面捕獲效率Pr被定義為

    (15)

    式中:Nd為沉積在管表面的顆粒數(shù)量;Nc為顆粒撞擊壁面的次數(shù)。

    顆粒沉積率是指單位時(shí)間內(nèi)顆粒在管道壁面上沉積的質(zhì)量或數(shù)量,粗糙表面會(huì)導(dǎo)致顆粒更容易沉積,因?yàn)楸砻娴陌纪菇Y(jié)構(gòu)會(huì)產(chǎn)生渦旋和湍流,減慢顆粒的運(yùn)動(dòng)速度,增加顆粒與壁面的接觸時(shí)間和頻率。沉積率Pd表示為

    (16)

    式中,Nt為追蹤到的顆??倲?shù)。

    為了探究不同表面形貌下的管束積灰特性,分別模擬了粗糙度為0、0.25、0.60、1.50 mm 4種不同粗糙程度的管束,入口煙氣速度為5 m/s,顆粒注入速度與煙氣速度一致。

    圖9所示為不同管束下顆粒沉積率并與Li等[12]的研究數(shù)據(jù)進(jìn)行對(duì)比。由圖可見,本文中的結(jié)果與Li等的結(jié)果變化趨勢(shì)一致,表明了數(shù)值模型的可靠性。隨著管道表面粗糙程度的增加,各個(gè)粒徑顆粒沉積率均有不同程度的增加,當(dāng)Rq=1.50 mm時(shí),各粒徑顆粒沉積率最大。出現(xiàn)這種現(xiàn)象的主要原因有2個(gè):一方面,粗糙表面會(huì)增加流體中顆粒的運(yùn)動(dòng)路徑和強(qiáng)度,從而提高顆粒在流體中與管道壁面碰撞的頻率,導(dǎo)致沉積率增大;另一方面,粗糙表面形成的局部渦流會(huì)影響近壁面湍流運(yùn)動(dòng),顆粒會(huì)在流體攜帶作用下進(jìn)入湍流渦旋。渦旋區(qū)域流速較低,顆粒受到運(yùn)動(dòng)阻力導(dǎo)致速度降低,因此顆粒會(huì)以較慢的速度撞擊管壁并且有較大的可能性產(chǎn)生沉積。

    圖10為不同管束下的壁面捕獲效率。由圖可見,隨著顆粒粒徑的增大,壁面捕獲效率隨之降低。對(duì)比Rq=0mm的光滑表面,Rq=1.50mm的粗糙表面捕獲效率明顯提升,說明表面粗糙度對(duì)顆粒沉積有顯著影響。其中粒徑為1~3μm的小顆粒表現(xiàn)出較高的壁面捕獲率,這是由于小粒徑顆粒慣性較小,更容易受流體影響,導(dǎo)致其更容易沉積在表面。對(duì)于較大顆粒(粒徑為10~20μm),捕獲效率隨著粗糙度的增加變化不大。這是因?yàn)檩^大粒徑顆粒具有更大的慣性,盡管粗糙度增加會(huì)影響流體運(yùn)動(dòng),但顆粒更大的慣性使其更容易脫離流體束縛,減少了顆粒沉積在壁面的機(jī)會(huì),從而導(dǎo)致捕獲效率降低。

    3.3表面形貌對(duì)顆粒碰壁的影響

    碰撞提升率是用于衡量粗糙壁面與光滑壁面之間的碰撞次數(shù)增加的比例,這個(gè)率值反映了粗糙表面結(jié)構(gòu)對(duì)顆粒沉積行為的影響。為了便于分析,定義碰撞提升率為

    ×100%,(17)

    式中:η為碰撞提升率;Nr表示粗糙壁碰撞次數(shù);Ns表示光滑壁碰撞次數(shù)。

    慣性撞擊是中低溫狀態(tài)下顆粒沉積的主要機(jī)制,因此顆粒碰撞壁面的次數(shù)也是影響顆粒沉積率的重要因素。粗糙壁面影響顆粒對(duì)壁面碰撞次數(shù)的主要因素有2種:首先,粗糙表面凹凸不平的輪廓相對(duì)增大了管束表面的面積,從而增大了顆粒與壁面的撞擊概率,其次,近壁面的渦旋結(jié)構(gòu)會(huì)吸引顆粒與壁面發(fā)生碰撞[17]。

    使用UDF壁面函數(shù)記錄顆粒碰壁次數(shù),輸出數(shù)據(jù)并結(jié)合公式(17)計(jì)算碰撞提升率。為了更直觀地觀察不同管壁表面粗糙度對(duì)顆粒行為的影響,圖11對(duì)比了不同粒徑的顆粒在光滑表面(Rq=0 mm)和粗糙表面(Rq=1.50 mm)時(shí)的管壁碰撞次數(shù)與碰撞提升率。由圖可見,粗糙表面的顆粒碰撞次數(shù)均高于光滑表面的,說明管壁表面粗糙度對(duì)顆粒碰撞行為具有明顯的影響。隨著顆粒粒徑的增加,總體碰撞提升率呈現(xiàn)先增大再減小的趨勢(shì),其中顆粒粒徑為1~3 μm的小顆粒碰撞提升率較低,這是由于小顆粒具有較小的慣性,使得其碰撞行為主要受到渦旋作用的影響,管壁表面粗糙度對(duì)其碰撞提升作用相對(duì)有限。顆粒平均粒徑為10 μm的顆粒具有最高的碰撞提升率,而顆粒粒徑為10~20 μm的大顆粒管壁碰撞次數(shù)與碰撞提升率開始下降,這是由于大粒徑顆粒慣性較大,這些顆粒更容易反彈,且受到渦旋作用的影響相對(duì)有限,因此碰撞提升率有所降低。

    4結(jié)論

    本文中基于分形理論建立了不同形貌的換熱管壁模型,采用計(jì)算流體力學(xué)方法結(jié)合臨界沉積速度模型,探究不同粗糙度的換熱管壁面對(duì)煙氣流動(dòng)、顆粒碰撞和顆粒沉積的影響,揭示表面粗糙度與顆粒沉積特性之間的內(nèi)在關(guān)聯(lián),為換熱器管束表面優(yōu)化設(shè)計(jì)提供指導(dǎo)依據(jù),同時(shí)為氣固兩相流動(dòng)以及顆粒沉積理論提供新的數(shù)據(jù)支撐和理論基礎(chǔ)。

    1)隨著換熱管表面粗糙度的增加,管道內(nèi)流體的壓降顯著增大,粗糙表面引起的渦流會(huì)使流體與壁面的摩擦阻力增大,導(dǎo)致流動(dòng)過程中的能量損失增大。

    2)換熱管壁表面粗糙度對(duì)流體流速和湍動(dòng)能的影響較為顯著,換熱管壁面粗糙度Rq=1.50 mm時(shí),較大粗糙度表面會(huì)明顯增強(qiáng)壁面附近的湍流效應(yīng),導(dǎo)致湍動(dòng)能增加和流速減小。

    3)隨著換熱管表面粗糙度的增加,顆粒的沉積率與捕獲效率顯著提升,其中粒徑為1~3 μm的小顆粒因慣性較小而表現(xiàn)出較高的沉積率與壁面捕獲率,粒徑為10~20 μm的顆粒慣性較大,容易脫離流體流動(dòng),導(dǎo)致其沉積率和捕獲效率受粗糙度的影響較小。

    4)飛灰顆粒的沉積與管束粗糙表面的形成具有正反饋效應(yīng)。顆粒沉積在管束表面,促使管束表面變得更加粗糙,而表面粗糙程度的增大,會(huì)進(jìn)一步提高顆粒沉積率。

    利益沖突聲明(Conflict of Interests)

    所有作者聲明不存在利益沖突。

    All authors declare no relevant conflict of interests.

    作者貢獻(xiàn)(Authors’Contributions)

    周友行、謝寶安和高騰騰參與論文構(gòu)思和寫作,楊文佳、宮天宇和彭麗參與論文審閱和修改。所有作者均閱讀并同意了最終稿件的提交。

    ZHOUYouhang,XIEBaoan,andGAO Tengtengcontributed to theconceptualizationand writingof the manuscript. YANG Wenjia,GONG Tianyu,and PENG Li participated in the review and revision. All authors have read the final version of the paper and consented to its submission.

    參考文獻(xiàn)(References)

    [1]馬廣興,孫煜光,杜聰,等. 空氣源熱泵換熱器沙塵沖蝕磨損模擬研究[J]. 流體機(jī)械,2021,49(12):62-69.

    MA GX,SUN YG,DU C,et al. Numerical simulation on sand erosion wear of heat exchanger of airsource heat pump[J]. Fluid Machinery,2021,49(12):62-69.

    [2]AN R,ZHANG XB. A strategy of mixed tube bundle to improve anti-deposition performance and heat transfer[J]. Interna-tional Communications in Heat and Mass Transfer,2024,154:107450.

    [3]AARNES RJ,HAUGEN EN,ANDERSSON IH. Inertial particle impaction on acylinder in turbulent cross-flow at modest reynolds numbers[J]. International Journal of Multiphase Flow,2018,111:53-61.

    [4]彭巖,趙欽新,王為術(shù). 煙氣橫掠管束外表面的積灰特性[J]. 化工學(xué)報(bào),2018,69(12):5034-5041.

    PENG Y,ZHAO QX,WANG WS. Ash deposition characteristics of external surface of flue gas swept across tube bundles[J]. CIESC Journal,2018,69(12):5034-5041.

    [5]YANG X,LIN M,MA L,et al. Understanding the ash deposition formation in zhundong lignite combustion through dynamic cfd modelling analysis[J]. Fuel,2017,194:533-543.

    [6]LIN M,WANG ST,ZHAI ZD,et al. Unsteady cfd simulation on ash particle deposition and removal characteristics in tube banks:focusing on particle diameter,flow velocity,and temperature[J]. Journal of the Energy Institute,2020,93(4):1481-1494.

    [7]袁寶強(qiáng),趙哲華,王湛,等. 基于動(dòng)網(wǎng)格的飛灰沉積特性數(shù)值模擬[J]. 工程熱物理學(xué)報(bào),2019,40(6):1332-1337. YUAN BQ,ZHAO ZH,WANG Z,etal. CFD analysis of ash deposition performance based on dynamics mesh[J]. Journal of Engineering Thermophysics,2019,40(6):1332-1337.

    [8]唐嬋,張靖周. 微小顆粒橫掠圓管的數(shù)值分析[J]. 熱科學(xué)與技術(shù),2015,14(6):497-504.

    TANGC,ZHANGJZ. Numerical analysis of tiny particles sweep round tube bundles[J]. Journal of Thermal Science and Technology,2015,14(6):497-504.

    [9]LIANG YT,LIJ B,LONG XF,etal. A numerical simulation study of ash deposition in acirculating fluidized bed during zhundong lignite combustion[J]. Fuel,2023,333(P2).

    [10]AN R,ZHANG XB. Dynamic evaluation of ash deposition behavior and parametric study using an inverse distance weighted interpolation dynamic mesh method[J]. Powder Technology,2024,437:119569.

    [11]KLEINHANS U,WIELAND C,F(xiàn)RANDSEN FJ,et al. Ash formation and deposition in coal and biomass fired combustionsystems:progress and challenges in the field of ash particle sticking and rebound behavior[J]. Progress in Energy and Com-bustion Science,2018,68:65-168.

    [12]LIJ B,DU WJ,LIN C. Numerical simulation and experiment of gas-solid two phase flow and ash deposition on anovel heat transfer surface[J]. Applied Thermal Engineering,2017,113:1033-1046.

    [13]MAJUMDAR A,TIEN CL. Fractal characterization and simulation of rough surfaces[J]. Wear,1990,136(2):313-327.

    [14]FU Y,ZHENG ZY,XIAO R,et al. Comparison of two fractal interpolation methods[J]. Physica A:Statistical Mechanics and its Applications,2017,469:563-571.

    [15]LI SL,YANG XL,JIN GD,et al. Wall-resolved large-eddy simulation of turbulent channel flows with rough walls[J]. Theoretical and Applied Mechanics Letters,2021,11(1):32-36.

    [16]MARCHIS DM,MILICI B,NAPOLI E. Numerical observations of turbulence structure modification in channel flow over 2D and 3D rough walls[J]. International Journal of Heat and Fluid Flow,2015,56:108-123.

    [17]洪文鵬,王欣. 粗糙流道微細(xì)顆粒物沉積特性的數(shù)值研究[J]. 中國電機(jī)工程學(xué)報(bào),2018,38(13):3889-3895,4030. HONG WP,WANG X. Numerical investigation on fine particles deposition in channels with rough surface[J]. Proceedings of the CSEE,2018,38(13):3889-3895,4030.

    Fractal characterization of heat exchange tube surfaces and numerical simulation on ash deposition characteristics

    ZHOU Youhang,XIE Baoan,GAO Tengteng,YANG Wenjia,GONG Tianyu,PENG Li School of Mechanical Engineering and Mechanics,University of Xiangtan,Xiangtan 411105,China

    Abstract

    ObjectiveThe surfaces of heat exchanger tube bundles become increasingly rougher due to prolonged use. To explore the deposi?tion characteristics of fly ash particles on heat exchanger tube surfaces with varying roughness and the mechanisms affecting this process,the deposition behavior of fly ash particles of different sizes under actual working conditions is simulated. The study aims to reveal the intrinsic relationship between particle deposition and surface roughness.

    MethodsTube bundle surface models with different roughness levels were developed using fractal theory and the improved Weierstrass-Mandelbrot(W-M)function. Thesurface models were generated usingMatlab,and mesh generation was con?ducted using ICEM software. Numerical simulations were performed in Fluent,combined with auser-defined function(UDF). The renormalization group(RNG)k-epsilon(k-ε)turbulence model was employed to simulate the incompressible turbulent flow,and the computational fluid dynamics-discrete phase model(CFD-DPM)was applied to track and analyze the motion of fly ash particles of different sizes in gas-solid two-phase flow. Particle deposition behavior was modeled using energy conservation and critical velocity theory.

    Results and DiscussionIncreasing surface roughness of the heat exchanger tube led to asignificant increase in both inlet pres?sure and pressure drop,indicating increased energy loss during fluid flow. This was attributed to the enhanced internal friction caused by highersurface roughness. Greatersurface roughnessalsointensified fluid velocityand turbulentkineticenergy. Reducing near-wall fluid velocity would in turn promote particle deposition. Rough surface peaks influenced fluid flow,generat?ing vortex structures on both windward and leeward sides,which further enhanced fluid-wall interactions. As surface roughness increased,the particle deposition rate rose significantly,especially for small particles with diameters of 1-3 μm. These particles showed higher deposition rates and wall capture efficiency due to their low inertia. As particle size increased,the number of particle-wall collisions and the overall collision enhancement rate first increased and then decreased. Compared to smooth sur?faces,the rough surfaces with an Rqof1. 50mm exhibited ahigher number of particle collisions,indicating the significant impact of surface roughness on particle collision behavior.

    ConclusionThe deposition of fly ash particles and the formation of rough surfaces on tube bundles exhibit apositive feedback loop. Particle deposition gradually makes the tube wall surfaces rougher,which in turn further accelerates particle deposition. Vortices generated by rough surfaces increase fluid-wall friction,resulting in greater energy loss,increased turbulent kinetic energy,and reduced flow velocity,thereby intensifying pipe wear and corrosion,ultimately shortening the equipment’s lifes?pan. Understanding the deposition mechanisms of fly ash particles under actual operating conditions allows for the prediction of deposition rates and distribution. This study provides atheoretical basis for optimizing tube bundle design and maintenance strat?egies,and offers new data and theoretical support for gas-solid two-phase flow and particle deposition theory.

    Keywords:fractal theory;fly ash deposition;numerical simulation;gas-solid two-phase flow

    (責(zé)任編輯:趙雁)

    猜你喜歡
    數(shù)值模擬
    基于AMI的雙色注射成型模擬分析
    錐齒輪精密冷擺輾成形在“材料成型數(shù)值模擬”課程教學(xué)中的應(yīng)用
    基于氣象信息及風(fēng)場(chǎng)信息的風(fēng)機(jī)輪轂處風(fēng)速預(yù)測(cè)
    鉆孔灌注樁樁底沉渣對(duì)樁體承載特性影響的模擬分析
    西南地區(qū)氣象資料測(cè)試、預(yù)處理和加工研究報(bào)告
    科技資訊(2016年18期)2016-11-15 08:01:18
    張家灣煤礦巷道無支護(hù)條件下位移的數(shù)值模擬
    科技視界(2016年18期)2016-11-03 23:14:27
    張家灣煤礦開切眼錨桿支護(hù)參數(shù)確定的數(shù)值模擬
    科技視界(2016年18期)2016-11-03 22:57:21
    跨音速飛行中機(jī)翼水汽凝結(jié)的數(shù)值模擬研究
    科技視界(2016年18期)2016-11-03 20:38:17
    姚橋煤礦采空區(qū)CO2防滅火的數(shù)值模擬分析
    雙螺桿膨脹機(jī)的流場(chǎng)數(shù)值模擬研究
    科技視界(2016年22期)2016-10-18 14:53:19
    国产亚洲欧美精品永久| 如何舔出高潮| 国产高清国产精品国产三级| 男女边摸边吃奶| 国产精品久久久久久久久免| 国产欧美另类精品又又久久亚洲欧美| 亚洲丝袜综合中文字幕| 狂野欧美激情性xxxx在线观看| 久久婷婷青草| 日本av手机在线免费观看| 不卡视频在线观看欧美| 最近手机中文字幕大全| 国产av精品麻豆| 成人二区视频| 日韩强制内射视频| 午夜日本视频在线| 精品99又大又爽又粗少妇毛片| 狂野欧美激情性xxxx在线观看| 日本av免费视频播放| 亚洲无线观看免费| 日本黄色片子视频| 18禁在线无遮挡免费观看视频| 69精品国产乱码久久久| 久久久国产精品麻豆| 久久狼人影院| 如何舔出高潮| 在线免费观看不下载黄p国产| 中文字幕av电影在线播放| 日韩中文字幕视频在线看片| av播播在线观看一区| 美女脱内裤让男人舔精品视频| 久热这里只有精品99| 国产91av在线免费观看| 亚洲国产色片| 成人无遮挡网站| 欧美日韩视频高清一区二区三区二| 女人精品久久久久毛片| 最新的欧美精品一区二区| .国产精品久久| 亚洲人成网站在线观看播放| 少妇精品久久久久久久| 美女视频免费永久观看网站| 麻豆成人午夜福利视频| 欧美日韩视频高清一区二区三区二| 大码成人一级视频| 成人国产麻豆网| 免费观看a级毛片全部| 男女免费视频国产| 亚洲性久久影院| 国产日韩欧美亚洲二区| 色婷婷久久久亚洲欧美| av天堂久久9| 我的女老师完整版在线观看| 成人毛片a级毛片在线播放| 欧美日本中文国产一区发布| 免费少妇av软件| 国产精品一二三区在线看| 欧美日韩国产mv在线观看视频| 男女国产视频网站| 国产极品粉嫩免费观看在线 | 欧美日韩在线观看h| 国产成人一区二区在线| 乱人伦中国视频| 欧美3d第一页| 黄色视频在线播放观看不卡| 久久韩国三级中文字幕| 午夜免费鲁丝| 美女主播在线视频| 国产一区二区三区av在线| 高清在线视频一区二区三区| 少妇 在线观看| 久久6这里有精品| 欧美日韩精品成人综合77777| 国产精品.久久久| 国产高清三级在线| 日韩一区二区三区影片| 亚洲综合精品二区| 观看美女的网站| 少妇人妻久久综合中文| 欧美亚洲 丝袜 人妻 在线| 精品国产一区二区三区久久久樱花| 欧美97在线视频| 亚洲国产成人一精品久久久| 一本久久精品| 久久韩国三级中文字幕| 午夜免费鲁丝| 亚州av有码| 一级二级三级毛片免费看| 国产精品一区二区在线观看99| 亚洲国产日韩一区二区| 国产成人免费无遮挡视频| 精品久久久久久电影网| 亚洲美女视频黄频| av播播在线观看一区| 高清黄色对白视频在线免费看 | 久久久久视频综合| 一区二区三区精品91| 日韩av在线免费看完整版不卡| 在线看a的网站| 国产男女超爽视频在线观看| 一个人看视频在线观看www免费| 少妇被粗大的猛进出69影院 | 色5月婷婷丁香| 日韩视频在线欧美| 在线观看一区二区三区激情| 久久精品国产a三级三级三级| 国产无遮挡羞羞视频在线观看| 一级黄片播放器| 精品人妻熟女av久视频| 国内精品宾馆在线| 女性被躁到高潮视频| 久久午夜福利片| 卡戴珊不雅视频在线播放| 成人特级av手机在线观看| 日韩三级伦理在线观看| 各种免费的搞黄视频| 免费播放大片免费观看视频在线观看| 亚洲精品久久久久久婷婷小说| 欧美xxⅹ黑人| 国产黄片美女视频| 亚洲经典国产精华液单| 日本黄色日本黄色录像| 色5月婷婷丁香| 秋霞伦理黄片| 精品亚洲成国产av| 高清在线视频一区二区三区| av线在线观看网站| 久久久久久人妻| 99热6这里只有精品| 热re99久久国产66热| 国产av国产精品国产| 国产黄片视频在线免费观看| 青春草亚洲视频在线观看| 高清av免费在线| 国产成人免费观看mmmm| 亚洲欧美成人综合另类久久久| 精华霜和精华液先用哪个| 久久99蜜桃精品久久| av播播在线观看一区| 亚洲精品日韩在线中文字幕| 国产精品福利在线免费观看| 在线天堂最新版资源| 国产熟女欧美一区二区| 不卡视频在线观看欧美| 一个人看视频在线观看www免费| 中国美白少妇内射xxxbb| h视频一区二区三区| 欧美精品国产亚洲| 日韩中字成人| 国产精品偷伦视频观看了| 91午夜精品亚洲一区二区三区| 午夜福利网站1000一区二区三区| 制服丝袜香蕉在线| 伊人亚洲综合成人网| 一本色道久久久久久精品综合| 欧美性感艳星| 国精品久久久久久国模美| 王馨瑶露胸无遮挡在线观看| 五月天丁香电影| 99热这里只有是精品在线观看| 国产一级毛片在线| 自线自在国产av| 丁香六月天网| h视频一区二区三区| 熟女av电影| 欧美+日韩+精品| www.av在线官网国产| 精品少妇黑人巨大在线播放| 性高湖久久久久久久久免费观看| 伦理电影大哥的女人| 久久女婷五月综合色啪小说| 亚洲欧美精品自产自拍| 交换朋友夫妻互换小说| 国产精品国产三级专区第一集| 欧美日韩视频精品一区| 日本wwww免费看| 久久久久视频综合| 黑人猛操日本美女一级片| 成人漫画全彩无遮挡| 黄色怎么调成土黄色| 午夜福利视频精品| 王馨瑶露胸无遮挡在线观看| 五月天丁香电影| 午夜激情福利司机影院| 妹子高潮喷水视频| 日日啪夜夜爽| 80岁老熟妇乱子伦牲交| 日本黄色日本黄色录像| 欧美三级亚洲精品| 久久久久久久久久久久大奶| 久久99蜜桃精品久久| 又大又黄又爽视频免费| 亚洲综合色惰| 黑人高潮一二区| 国产视频首页在线观看| 亚洲精品久久午夜乱码| 久久午夜福利片| 人妻少妇偷人精品九色| 亚洲成人一二三区av| 日日啪夜夜爽| 久热这里只有精品99| 尾随美女入室| 久久人妻熟女aⅴ| 免费高清在线观看视频在线观看| 少妇人妻一区二区三区视频| 3wmmmm亚洲av在线观看| 最近最新中文字幕免费大全7| 有码 亚洲区| 黄色怎么调成土黄色| 婷婷色综合www| 亚洲欧美中文字幕日韩二区| 超碰97精品在线观看| 熟妇人妻不卡中文字幕| 欧美一级a爱片免费观看看| 内射极品少妇av片p| 精品酒店卫生间| 国产有黄有色有爽视频| 亚洲av国产av综合av卡| 日韩不卡一区二区三区视频在线| 久久久久久久久大av| 黑人猛操日本美女一级片| 免费观看无遮挡的男女| 国产老妇伦熟女老妇高清| 天堂俺去俺来也www色官网| 91成人精品电影| 少妇的逼水好多| 亚洲国产精品999| 亚洲经典国产精华液单| 天天操日日干夜夜撸| 欧美亚洲 丝袜 人妻 在线| 亚洲国产精品专区欧美| 人人澡人人妻人| 中文字幕精品免费在线观看视频 | 日韩熟女老妇一区二区性免费视频| 永久免费av网站大全| 在线观看免费日韩欧美大片 | av免费观看日本| 大码成人一级视频| 久久97久久精品| 中文字幕精品免费在线观看视频 | 亚洲欧美精品专区久久| 爱豆传媒免费全集在线观看| 免费大片18禁| 精品一区二区三区视频在线| 又爽又黄a免费视频| 中文资源天堂在线| 久久人妻熟女aⅴ| 婷婷色综合大香蕉| 大片电影免费在线观看免费| 亚洲色图综合在线观看| 狂野欧美激情性bbbbbb| 欧美性感艳星| 男人和女人高潮做爰伦理| 久久久久久久国产电影| 国产伦精品一区二区三区四那| 亚洲av二区三区四区| 国产午夜精品久久久久久一区二区三区| 久久久久视频综合| 在线观看免费高清a一片| 嘟嘟电影网在线观看| 国产在线视频一区二区| 国产黄色免费在线视频| 亚洲国产精品专区欧美| 9色porny在线观看| 在线 av 中文字幕| 国产爽快片一区二区三区| 日韩精品免费视频一区二区三区 | 精品国产露脸久久av麻豆| 亚洲精品视频女| 一级毛片 在线播放| 9色porny在线观看| 大香蕉97超碰在线| 另类精品久久| 少妇被粗大的猛进出69影院 | 亚洲天堂av无毛| 少妇人妻久久综合中文| 中国国产av一级| 亚洲av免费高清在线观看| 欧美 日韩 精品 国产| 免费人成在线观看视频色| 黄色日韩在线| 精华霜和精华液先用哪个| 激情五月婷婷亚洲| 性色avwww在线观看| 啦啦啦啦在线视频资源| 日韩电影二区| 日韩av在线免费看完整版不卡| 日本爱情动作片www.在线观看| 久久久久久久久久久丰满| 国产免费一级a男人的天堂| 亚洲情色 制服丝袜| 丝瓜视频免费看黄片| 纵有疾风起免费观看全集完整版| 黄色怎么调成土黄色| 乱码一卡2卡4卡精品| 夫妻性生交免费视频一级片| 中文欧美无线码| 国国产精品蜜臀av免费| 热re99久久国产66热| 国产精品国产三级国产专区5o| 18+在线观看网站| 男的添女的下面高潮视频| 99热网站在线观看| 久久99热这里只频精品6学生| 三级国产精品片| 人人妻人人看人人澡| 三级国产精品欧美在线观看| 日韩中字成人| 少妇裸体淫交视频免费看高清| 国产老妇伦熟女老妇高清| 老女人水多毛片| 久久女婷五月综合色啪小说| 哪个播放器可以免费观看大片| 97在线视频观看| 简卡轻食公司| 成人综合一区亚洲| 在线天堂最新版资源| 日韩不卡一区二区三区视频在线| 9色porny在线观看| 一级黄片播放器| 久久婷婷青草| 寂寞人妻少妇视频99o| 草草在线视频免费看| 色视频在线一区二区三区| 国产成人aa在线观看| 赤兔流量卡办理| 欧美精品人与动牲交sv欧美| 99热网站在线观看| av卡一久久| 免费黄网站久久成人精品| www.av在线官网国产| 国产一区二区三区av在线| 亚洲国产精品一区三区| 麻豆精品久久久久久蜜桃| 在线观看免费视频网站a站| 国产永久视频网站| 中国三级夫妇交换| 色婷婷av一区二区三区视频| 最近中文字幕2019免费版| 丝袜喷水一区| 日本欧美国产在线视频| 亚洲av国产av综合av卡| 你懂的网址亚洲精品在线观看| 亚洲天堂av无毛| 亚洲情色 制服丝袜| 狂野欧美白嫩少妇大欣赏| 久久毛片免费看一区二区三区| 中文天堂在线官网| 少妇的逼好多水| 在现免费观看毛片| 欧美一级a爱片免费观看看| 国产有黄有色有爽视频| av免费在线看不卡| 一级毛片我不卡| 高清在线视频一区二区三区| 国产男女超爽视频在线观看| 国产 一区精品| 视频中文字幕在线观看| 人妻系列 视频| 日日摸夜夜添夜夜添av毛片| 日产精品乱码卡一卡2卡三| 久久影院123| 性色av一级| 国产在线视频一区二区| 精品国产乱码久久久久久小说| 国产精品久久久久久av不卡| 亚洲av在线观看美女高潮| 免费av不卡在线播放| 日日撸夜夜添| 日韩免费高清中文字幕av| 免费观看无遮挡的男女| 纯流量卡能插随身wifi吗| h日本视频在线播放| 久久久久久久久大av| 亚洲精品国产av成人精品| tube8黄色片| 亚洲欧美精品自产自拍| 欧美精品国产亚洲| 久久精品夜色国产| 日韩成人伦理影院| 最近最新中文字幕免费大全7| 高清欧美精品videossex| av在线老鸭窝| 交换朋友夫妻互换小说| 青春草亚洲视频在线观看| 欧美 日韩 精品 国产| 人体艺术视频欧美日本| 国产一区亚洲一区在线观看| 亚洲国产精品成人久久小说| 久久婷婷青草| 免费高清在线观看视频在线观看| 偷拍熟女少妇极品色| 在线亚洲精品国产二区图片欧美 | 欧美精品高潮呻吟av久久| 高清黄色对白视频在线免费看 | 自线自在国产av| 在线天堂最新版资源| 大片电影免费在线观看免费| 日韩 亚洲 欧美在线| 欧美高清成人免费视频www| av播播在线观看一区| 国产高清有码在线观看视频| 纯流量卡能插随身wifi吗| 99视频精品全部免费 在线| 各种免费的搞黄视频| 我要看黄色一级片免费的| av专区在线播放| 日韩av在线免费看完整版不卡| 在线观看国产h片| 免费看日本二区| 国产精品.久久久| 日韩三级伦理在线观看| 在现免费观看毛片| 日韩欧美一区视频在线观看 | 在现免费观看毛片| 欧美一级a爱片免费观看看| 中文天堂在线官网| 伦理电影免费视频| 欧美成人精品欧美一级黄| 免费av不卡在线播放| 久久久久视频综合| 少妇人妻精品综合一区二区| 亚洲无线观看免费| 一级毛片 在线播放| 麻豆成人午夜福利视频| 蜜桃在线观看..| 欧美精品一区二区免费开放| 六月丁香七月| 美女中出高潮动态图| 欧美精品一区二区免费开放| 卡戴珊不雅视频在线播放| 在线观看一区二区三区激情| 亚洲av中文av极速乱| 视频中文字幕在线观看| 又爽又黄a免费视频| 久久热精品热| 中文字幕av电影在线播放| tube8黄色片| 这个男人来自地球电影免费观看 | 欧美高清成人免费视频www| 五月伊人婷婷丁香| 一本久久精品| 国产又色又爽无遮挡免| 日韩伦理黄色片| 五月天丁香电影| 久久av网站| 亚洲av.av天堂| 狂野欧美白嫩少妇大欣赏| 国产精品不卡视频一区二区| 五月伊人婷婷丁香| 人妻人人澡人人爽人人| 午夜视频国产福利| 国产熟女午夜一区二区三区 | 最近中文字幕高清免费大全6| 中文在线观看免费www的网站| 国产极品天堂在线| 亚洲婷婷狠狠爱综合网| h视频一区二区三区| 日本午夜av视频| 国产片特级美女逼逼视频| 99re6热这里在线精品视频| 日韩中文字幕视频在线看片| 建设人人有责人人尽责人人享有的| 五月开心婷婷网| 亚洲国产欧美在线一区| 成人二区视频| av.在线天堂| 人人澡人人妻人| 99热这里只有是精品50| 国产黄频视频在线观看| 高清欧美精品videossex| 色5月婷婷丁香| 亚洲精品乱久久久久久| 建设人人有责人人尽责人人享有的| 久久女婷五月综合色啪小说| 国产亚洲最大av| 乱码一卡2卡4卡精品| 自拍偷自拍亚洲精品老妇| 午夜老司机福利剧场| 在线播放无遮挡| 久久久久久伊人网av| 狂野欧美白嫩少妇大欣赏| 三级经典国产精品| 少妇人妻 视频| 亚洲四区av| 欧美日本中文国产一区发布| 大码成人一级视频| 久久精品国产鲁丝片午夜精品| freevideosex欧美| 日本欧美国产在线视频| 大香蕉久久网| 伦理电影免费视频| 久久久久久久亚洲中文字幕| 精华霜和精华液先用哪个| 亚洲欧洲日产国产| 亚洲国产成人一精品久久久| 2022亚洲国产成人精品| av不卡在线播放| 伦理电影大哥的女人| 18禁在线播放成人免费| 日韩电影二区| av免费在线看不卡| 纯流量卡能插随身wifi吗| 日日啪夜夜撸| 老司机影院成人| 免费av不卡在线播放| 久久精品久久精品一区二区三区| 99re6热这里在线精品视频| 午夜激情久久久久久久| 久久精品国产亚洲av天美| 国产精品无大码| 亚洲欧美日韩另类电影网站| 成年美女黄网站色视频大全免费 | 有码 亚洲区| videos熟女内射| 国产欧美日韩综合在线一区二区 | 成人综合一区亚洲| 视频区图区小说| 在线观看www视频免费| 22中文网久久字幕| 丰满少妇做爰视频| 国产精品一区二区三区四区免费观看| 国产视频首页在线观看| 免费少妇av软件| 七月丁香在线播放| 蜜桃久久精品国产亚洲av| 国产真实伦视频高清在线观看| 水蜜桃什么品种好| 精品国产一区二区三区久久久樱花| 国产黄片视频在线免费观看| 另类亚洲欧美激情| 九色成人免费人妻av| 观看免费一级毛片| 国产伦精品一区二区三区视频9| 国产精品久久久久久精品古装| 久久久午夜欧美精品| 欧美日韩av久久| 国产av码专区亚洲av| 在线观看免费高清a一片| 亚洲国产精品999| av一本久久久久| 精品一区二区三卡| 亚洲av男天堂| 免费黄网站久久成人精品| 国产熟女午夜一区二区三区 | 美女内射精品一级片tv| 亚洲久久久国产精品| 丰满迷人的少妇在线观看| 精品视频人人做人人爽| 一本色道久久久久久精品综合| 亚洲美女搞黄在线观看| 五月伊人婷婷丁香| 男人狂女人下面高潮的视频| freevideosex欧美| 久久久久国产精品人妻一区二区| 国产女主播在线喷水免费视频网站| 插阴视频在线观看视频| av天堂中文字幕网| 在线观看免费高清a一片| 成年人免费黄色播放视频 | 成年人午夜在线观看视频| 亚洲第一区二区三区不卡| 午夜91福利影院| h日本视频在线播放| 亚洲精品久久午夜乱码| 在线播放无遮挡| 国国产精品蜜臀av免费| h视频一区二区三区| 亚洲内射少妇av| 青春草视频在线免费观看| 大片电影免费在线观看免费| 夜夜爽夜夜爽视频| 亚洲一级一片aⅴ在线观看| 国产日韩欧美亚洲二区| 国产成人91sexporn| 日日啪夜夜爽| 国产一区二区在线观看av| 久久精品国产a三级三级三级| 少妇人妻 视频| 老熟女久久久| 搡女人真爽免费视频火全软件| 特大巨黑吊av在线直播| 黄片无遮挡物在线观看| 99热这里只有是精品在线观看| av在线观看视频网站免费| 国产乱来视频区| 亚洲精品乱码久久久久久按摩| 亚洲欧美中文字幕日韩二区| 91久久精品电影网| 黄色视频在线播放观看不卡| av网站免费在线观看视频| 久久久久精品性色| av网站免费在线观看视频| 日韩 亚洲 欧美在线| 国产黄色免费在线视频| 日韩一区二区三区影片| 久久99精品国语久久久| 日本黄色片子视频| 在线观看免费日韩欧美大片 | 在线 av 中文字幕| 寂寞人妻少妇视频99o| 国产精品熟女久久久久浪| 男女免费视频国产| 99久国产av精品国产电影| 日本黄大片高清| 曰老女人黄片| 男人舔奶头视频| 天堂俺去俺来也www色官网| 如日韩欧美国产精品一区二区三区 | 青春草视频在线免费观看| 国产男人的电影天堂91| 日韩欧美一区视频在线观看 | 国产高清三级在线| 在线观看三级黄色| 国产伦精品一区二区三区四那|