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

    基于二次攜帶現(xiàn)象的折流板除霧性能仿真優(yōu)化

    2021-12-07 08:13:20劉定平張宇駿
    關(guān)鍵詞:流板液膜液滴

    劉定平 張宇駿

    (華南理工大學(xué) 電力學(xué)院∥廣東省能源高效低污染轉(zhuǎn)化工程技術(shù)研究中心,廣東 廣州510640)

    濕法煙氣脫硫(WFGD)系統(tǒng)由于效率高、成本低、運(yùn)行穩(wěn)定,廣泛應(yīng)用于火電、化工、水泥、陶瓷等行業(yè)。然而,經(jīng)WFGD系統(tǒng)處理后的煙氣會(huì)攜帶大量酸性液滴。這些液滴不僅會(huì)對(duì)系統(tǒng)下游設(shè)備造成腐蝕和堵塞,而且排至大氣后會(huì)形成“煙囪雨”現(xiàn)象,對(duì)環(huán)境造成嚴(yán)重污染[1-2]。因此,開(kāi)發(fā)性能優(yōu)良的除霧器,解決“煙囪雨”問(wèn)題是當(dāng)前研究的熱點(diǎn)。

    CFD數(shù)值模擬是研究除霧器性能的重要手段。許多學(xué)者[3-7]運(yùn)用數(shù)值模擬對(duì)除霧器結(jié)構(gòu)及性能進(jìn)行了優(yōu)化研究,得出的結(jié)論對(duì)提高除霧器性能具有一定指導(dǎo)意義。但是,由于其數(shù)學(xué)模型的局限性,上述研究均側(cè)重于結(jié)構(gòu)參數(shù)對(duì)液滴捕集規(guī)律的探究,很少關(guān)注被捕集液滴在壁面形成液膜后,在高速氣流作用下液滴飛濺和液膜分離等造成的液滴二次攜帶現(xiàn)象。這也使得研究結(jié)論總是一致的表明:更高的入口流速有更佳的除霧效率[8]。而在實(shí)際應(yīng)用中,當(dāng)入口流速超過(guò)臨界速度時(shí)會(huì)發(fā)生二次攜帶現(xiàn)象[9],導(dǎo)致除霧器性能惡化,使優(yōu)化后的除霧器在高流速工況下失效。

    針對(duì)上述問(wèn)題,本文在傳統(tǒng)歐拉-拉格朗日方法的基礎(chǔ)上耦合歐拉壁膜(EWF)模型,建立了基于二次攜帶的除霧模型。以實(shí)驗(yàn)參數(shù)[10]作為模擬參數(shù),對(duì)比模擬結(jié)果與實(shí)驗(yàn)結(jié)果驗(yàn)證了所構(gòu)建模型的可靠性;通過(guò)仿真研究了折流板在不同葉片板間距和轉(zhuǎn)折角度情況下的除霧效率、臨界速度及壓降等參數(shù)的變化規(guī)律,并進(jìn)行了流場(chǎng)與結(jié)構(gòu)優(yōu)化,得到了除霧器的結(jié)構(gòu)優(yōu)化參數(shù)。

    1 物理模型

    本次研究選取折流板除霧器的單一通道作為研究對(duì)象。模型為三維結(jié)構(gòu),平面結(jié)構(gòu)特征如圖1所示,主要參數(shù)如表1所示,寬度為150 mm。為了便于模型驗(yàn)證,結(jié)構(gòu)參數(shù)與文獻(xiàn)[10]保持一致。

    圖1 折流板除霧器結(jié)構(gòu)Fig.1 Structural of wave-plate demister

    表1 折流板主要尺寸參數(shù)[10]Table 1 Dimensional parameters of wave-plate demister[10]

    本文應(yīng)用ICEM CFD軟件對(duì)幾何模型進(jìn)行網(wǎng)格劃分,采用結(jié)構(gòu)型六面體網(wǎng)格。為了保證計(jì)算精度的同時(shí)避免過(guò)大的計(jì)算量,在進(jìn)行網(wǎng)格無(wú)關(guān)性驗(yàn)證后,確定了最佳網(wǎng)格劃分方案:最大網(wǎng)格尺寸為0.2 mm,第1層邊界層高度及膨脹比率分別為0.06 mm和1.025,最終網(wǎng)格數(shù)量為356 178。

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

    2.1 模型假設(shè)

    本文在傳統(tǒng)歐拉-拉格朗日法的基礎(chǔ)上耦合EWF模型構(gòu)建基于二次攜帶現(xiàn)象的折流板除霧數(shù)學(xué)模型:采用歐拉法(Realizablek-ε湍流模型)計(jì)算煙氣在除霧器內(nèi)的流動(dòng);采用拉格朗日法(離散相模型)計(jì)算液滴的軌跡;采用EWF模型計(jì)算液膜發(fā)展。仿真基于以下假設(shè)[6,11-12]:

    (1)除霧器理想工況下煙氣各參數(shù)不隨時(shí)間變化,因此將煙氣流動(dòng)視為定常流動(dòng);

    (2)煙氣在除霧器內(nèi)的流速低于0.1個(gè)馬赫數(shù),可視為不可壓縮流。選用空氣模擬煙氣;

    (3)將液滴視為均質(zhì)的硬球模型,酸性液滴簡(jiǎn)化為水滴,忽略液滴間的相互作用力,即液滴間無(wú)摩擦、碰撞及破碎的發(fā)生。

    2.2 控制方程

    2.2.1 除霧器內(nèi)液滴運(yùn)動(dòng)的控制方程

    液滴在流場(chǎng)中運(yùn)動(dòng)主要受各種力支配,如曳力、重力等。因此,液滴軌跡可通過(guò)牛頓第二定律求解。離散相模型中,液滴控制方程[13]如下:

    (1)

    式中,u和ρ分別為煙氣的流速和密度;up、ρp和dp分別為液滴的速度、密度和粒徑;μ為煙氣的動(dòng)力粘度;CD和Rer分別為液滴的曳力系數(shù)和相對(duì)雷諾數(shù);F為作用于液滴上的附加力,由于ρp/ρ?1,因此考慮Saffman升力[14]。

    離散相計(jì)算不僅要考慮煙氣流動(dòng)對(duì)液滴的作用,還應(yīng)考慮液滴運(yùn)動(dòng)對(duì)流場(chǎng)的作用。為了貼合實(shí)際流動(dòng),本文對(duì)連續(xù)相(煙氣)和離散相(液滴)進(jìn)行雙向耦合計(jì)算。

    在雙向耦合中[15],連續(xù)相與離散相交替迭代計(jì)算,直至兩相方程收斂。在此過(guò)程中,連續(xù)相傳遞給離散相的動(dòng)量值可通過(guò)計(jì)算離散相的動(dòng)量變化量獲得,并在隨后的連續(xù)相動(dòng)量守恒方程中以源項(xiàng)的形式反映。在本文中,液滴的動(dòng)量變化值Fk為

    (2)

    2.2.2 除霧器內(nèi)液膜運(yùn)動(dòng)的控制方程

    液膜在壁面的發(fā)展及運(yùn)動(dòng)是通過(guò)EWF模型控制。在模型假設(shè)中,相比于壁面曲率液膜的厚度足夠小,因此可忽略液膜物性在沿厚度方向的變化。

    液膜的質(zhì)量守恒方程[16]如下:

    (3)

    式中:ρl為液膜的密度;h為液膜的厚度;s為表面梯度算子;Vl為平均液膜速度;為由于液滴捕集、液滴撞擊液膜產(chǎn)生飛濺、液膜分離等現(xiàn)象引起的單位壁面液膜質(zhì)量變化的源項(xiàng)。

    液膜的動(dòng)量守恒方程[17]如下:

    (4)

    (5)

    (6)

    2.2.3 液滴飛濺機(jī)制

    以往的研究大多只關(guān)注除霧器葉片對(duì)液滴的捕集特性,少有考慮被捕集液滴在壁面的成膜特性和液滴對(duì)液膜的慣性碰撞。當(dāng)液滴撞擊液膜時(shí),液滴可以發(fā)生粘附、鋪展或飛濺,這均取決于液滴撞擊液膜的碰撞能量E,可表達(dá)為[19-20]:

    (7)

    式中:Vp,n為液滴相對(duì)于被撞擊壁面的法向速度;σ為液膜的表面張力;雷諾數(shù)Re=ρpVp,ndp/μ。

    當(dāng)E<16時(shí),液滴在撞擊壁面后將以較為完整的球型粘附在液膜;當(dāng)1657.7時(shí),液滴撞擊將會(huì)發(fā)生飛濺,并生成新的液滴進(jìn)入流場(chǎng)。在上述機(jī)理中,粘附和鋪展控制液膜的發(fā)展,飛濺則是造成二次攜帶現(xiàn)象的重要原因。

    2.2.4 液膜分離機(jī)制

    除了液滴飛濺,液膜的分離也是造成除霧器二次攜帶現(xiàn)象的重要原因。在重力、煙氣剪切力等驅(qū)動(dòng)力作用下,當(dāng)液膜流經(jīng)壁面拐角處,只要液膜運(yùn)動(dòng)的慣性足以克服維持其膜狀的表面張力,新的液滴就會(huì)在拐角邊緣被分離出,并隨煙氣帶離流道,如圖2所示。

    圖2 剪切力作用下液膜在轉(zhuǎn)角處分離示意圖[21]Fig.2 Schematic of shear-driven film separation atcorners[21]

    液膜的分離需要同時(shí)滿足兩個(gè)條件[22]:轉(zhuǎn)折角θ足夠大;液膜的韋伯?dāng)?shù)Wef大于發(fā)生分離的臨界韋伯?dāng)?shù),即同時(shí)滿足式(8)和(9):

    θ>θcritic

    (8)

    Wef>Wef,critic

    (9)

    2.3 邊界條件與計(jì)算方法

    本文采用歐拉-拉格朗日法耦合EWF模型對(duì)折流板除霧器的除霧性能進(jìn)行仿真。按照文獻(xiàn)[10]的實(shí)驗(yàn)參數(shù)對(duì)仿真進(jìn)行設(shè)定。

    連續(xù)相關(guān)鍵參數(shù)的設(shè)置如下:煙氣使用密度為1.204 kg/m3和動(dòng)力粘度為1.825×10-5kg/(m·s)的不可壓縮空氣替代。采用速度進(jìn)口,入口速度U設(shè)定在3~7 m/s,方向垂直于煙道入口;采用壓力出口,出口表壓取0(操作壓力為10 135 Pa)。

    離散相關(guān)鍵參數(shù)的設(shè)置如下:液滴取密度為998.2 kg/m3、動(dòng)力粘度為1.002×10-3kg/(m·s)、表面張力為0.072 7 N/m的水替代。假定霧滴的噴射速度與煙氣流速相同,沿?zé)煔馊肟诿娣ㄏ蚓鶆蛏淙耄|(zhì)量流量Q=0.022 62UD,單位為kg/s。根據(jù)文獻(xiàn)[10]的實(shí)驗(yàn)參數(shù),入射液滴粒徑分布遵循Rosin-Rammler[23]分布,粒徑范圍集中在2~80 μm,平均粒徑為39.35 μm,分布系數(shù)為2.8。此外,考慮到湍流對(duì)液滴運(yùn)動(dòng)的影響,采用隨機(jī)游走模型[24]對(duì)液滴的運(yùn)動(dòng)軌跡進(jìn)行修正。

    仿真采用SIMPLE算法對(duì)壓力場(chǎng)與速度場(chǎng)進(jìn)行耦合。動(dòng)量方程、煙氣湍流動(dòng)能和湍流耗散率等均采用二階迎風(fēng)格式離散。液膜計(jì)算的時(shí)間步長(zhǎng)設(shè)為0.000 1 s,每計(jì)算液膜方程10次更新一次離散相方程。在計(jì)算過(guò)程中監(jiān)控入口面壓力、出口面流速及液膜厚度的變化,當(dāng)收斂殘差低于10-4且上述3個(gè)監(jiān)控值穩(wěn)定后,計(jì)算視為收斂。

    2.4 除霧器性能評(píng)價(jià)指標(biāo)

    文中研究不同葉片板間距和轉(zhuǎn)折角度情況下折流板的除霧性能。主要的評(píng)價(jià)指標(biāo)有除霧效率、臨界速度和進(jìn)出口壓降。

    除霧效率是表征除霧器除霧性能最關(guān)鍵的指標(biāo)。它直接反映除霧器脫去煙氣中液滴的能力。本文中,除霧效率η定義為

    (10)

    式中,Qin為除霧器進(jìn)口處液滴質(zhì)量流量,Qout為出口處液滴質(zhì)量流量。

    臨界速度是表征除霧器除霧性能穩(wěn)定性的重要指標(biāo)。通常來(lái)說(shuō),除霧效率會(huì)隨著煙氣入口流速的增大而提高。但由于煙氣對(duì)已捕集液滴的再攜帶作用,當(dāng)入口流速大于一定值時(shí),除霧效率將開(kāi)始下降。這時(shí),對(duì)應(yīng)于最高除霧效率的入口流速即為臨界速度。臨界速度越小,說(shuō)明除霧器穩(wěn)定性越差。

    鑒于在實(shí)際研究中往往只能得到臨界速度的近似值[25],因此本文取仿真效率曲線開(kāi)始下降前最高點(diǎn)對(duì)應(yīng)的入口流速值為臨界速度的近似值。

    除霧器進(jìn)出口壓降是評(píng)判除霧器除霧性能重要經(jīng)濟(jì)指標(biāo)。越大的壓降值意味著除霧系統(tǒng)所需的能耗越高,因此在除霧器設(shè)計(jì)及優(yōu)化時(shí)要在保證除霧效率的情況下盡可能降低除霧器壓降。本文中,除霧器進(jìn)出口壓降ΔP定義為

    ΔP=Pin-Pout

    (11)

    式中,Pin為除霧器入口面壓力,Pout為除霧器出口面壓力。

    2.5 模型驗(yàn)證

    為了驗(yàn)證模型的可靠性,以文獻(xiàn)[10]折流板除霧器為原型,通過(guò)改變煙氣入口流速得到不同工況下折流板的除霧效率和進(jìn)出口壓降,將文獻(xiàn)[10]中相同工況下的實(shí)驗(yàn)數(shù)據(jù)與文中仿真結(jié)果進(jìn)行比較。

    圖3為在不同入口流速下折流板除霧效率模擬結(jié)果和文獻(xiàn)[10]中實(shí)驗(yàn)及模擬結(jié)果的對(duì)比。從圖中可以看出,傳統(tǒng)的歐拉-拉格朗日方法未能預(yù)測(cè)出由于二次攜帶現(xiàn)象導(dǎo)致的除霧效率下降,而引入EWF模型改良后的歐拉-拉格朗日方法不僅能較為準(zhǔn)確的預(yù)測(cè)除霧效率隨入口流速的變化關(guān)系,而且預(yù)測(cè)誤差小,平均相對(duì)誤差僅為2.21%。

    如圖3所示,除霧效率隨著入口煙氣流速變化呈現(xiàn)先升高后下降的趨勢(shì),存在臨界速度。在U=3~6 m/s的變化范圍內(nèi),除霧效率隨著入口流速的增大而提升;當(dāng)入口流速達(dá)到U=6 m/s時(shí),除霧效率達(dá)到最大;當(dāng)入口流速為7 m/s時(shí),除霧效率開(kāi)始降低,這是由于除霧器內(nèi)已被捕集的液滴在高速氣流的作用下再次進(jìn)入煙氣,并隨煙氣流出除霧器,造成了二次攜帶現(xiàn)象。

    圖3 變流速下除霧效率模擬結(jié)果與文獻(xiàn)[10]結(jié)果對(duì)比Fig.3 Comparison of simulation results of demisting efficiency with the results of reference[10] under variable flow rate

    仿真計(jì)算結(jié)果表明,造成二次攜帶現(xiàn)象的主要原因是轉(zhuǎn)折角處煙氣對(duì)液膜的剪切分離作用。雖然液滴在流場(chǎng)中具有較高的速度,但其相對(duì)于被撞擊壁面的法向速度較小,難以發(fā)生飛濺。因此在本模型下,液滴大多以粘附、鋪展的方式形成液膜,飛濺可忽略不計(jì)。

    文獻(xiàn)[10]中實(shí)驗(yàn)所得和模型模擬所得的除霧器進(jìn)出口壓降比較見(jiàn)圖4。從圖中可以看出,模擬計(jì)算得到的進(jìn)出口壓降與文獻(xiàn)[10]中實(shí)驗(yàn)所測(cè)的變化趨勢(shì)基本一致,且數(shù)值十分接近,平均相對(duì)誤差僅為6.46%。綜合除霧效率和進(jìn)出口壓降的對(duì)比效果可得,文中采用的數(shù)值模型能夠較為準(zhǔn)確的預(yù)測(cè)折流板的除霧性能。

    圖4 變流速下進(jìn)出口壓降模擬結(jié)果與文獻(xiàn)[10]實(shí)驗(yàn)結(jié)果對(duì)比Fig.4 Comparison of simulation results of the pressure drop with experimental results in reference [10] under variable flow rate

    3 數(shù)值計(jì)算與分析

    3.1 板間距對(duì)除霧性能的影響

    為研究葉片板間距對(duì)折流板除霧性能的影響,在其他條件不變的情況下,選取D=20,25,30,35,40 mm 5種不同板間距進(jìn)行比較研究。其中,轉(zhuǎn)折角α=40°,直流段長(zhǎng)度H1=38 mm,捕集段長(zhǎng)度H2=49.6 mm。

    3.1.1 板間距對(duì)除霧效率的影響

    圖5示出了在不同入口流速、不同板間距下折流板除霧效率的變化曲線。由圖5可以看出,在同一入口流速條件下,減小板間距能有效提高折流板的除霧效率。當(dāng)D=20 mm時(shí),折流板在入口流速U的全范圍內(nèi)均具有最高的除霧效率,且除霧效率隨著U的增大而平穩(wěn)提高;在U=7 m/s時(shí)除霧效率達(dá)到最高值,為92.99%。

    圖5 不同流速和不同板間距下除霧器的效率Fig.5 Efficiency of demister under different flow rates and different plate spacing

    經(jīng)分析主要有兩方面原因:一方面,折流板是通過(guò)改變煙氣流動(dòng)方向使具有較大慣性的液滴撞擊捕集葉片實(shí)現(xiàn)除霧的,而較小的葉片間距縮短了慣性液滴撞擊捕集面的運(yùn)動(dòng)距離,提高了液滴被捕集的概率;另一方面,更高的入口流速意味著液滴具有更大的慣性。這使得液滴在氣流突然改變流向的情況下依舊能夠以慣性保持原有的運(yùn)動(dòng)軌跡,利于液滴的分離。

    然而,除霧效率并非總隨著入口流速的增大而提高。從圖5可以看出,當(dāng)D=30,35,40 mm時(shí),折流板在較高流速下均發(fā)生不同程度的二次攜帶現(xiàn)象,極大降低了除霧效率。以D=35 mm為例,隨著U由3 m/s提高至7 m/s,除霧效率先是穩(wěn)定提升,由59.44%提升至68.17%,在其臨界速度5 m/s時(shí)達(dá)到最大,之后便急速惡化,由最高的74.19%減低至51.41%和39.66%。

    此外,值得注意的是,除霧器的臨界速度并非為定值,而是隨著板間距的增大而減小。對(duì)于D=20,25 mm,臨界速度大于或等于7 m/s;對(duì)于D=30 mm,臨界速度為6 m/s;對(duì)于D=35,40 mm,臨界速度為5 m/s。臨界速度越低,說(shuō)明除霧器對(duì)工況的適應(yīng)性越差,在煙氣量發(fā)生波動(dòng)的工況下更易發(fā)生二次攜帶現(xiàn)象,導(dǎo)致除霧性能惡化。

    圖6為3種不同板間距結(jié)構(gòu)在入口流速為7 m/s的工況下的流場(chǎng)速度分布云圖對(duì)比。結(jié)果表明,煙氣流速在流經(jīng)轉(zhuǎn)角處具有最大值,且葉片板間距越大,轉(zhuǎn)折角處流速越高。如圖6所示,當(dāng)D=20 mm時(shí),轉(zhuǎn)角處最大流速為12.06 m/s;當(dāng)D=30 mm時(shí),轉(zhuǎn)角處最大流速為13.15 m/s;當(dāng)D=40 mm時(shí),轉(zhuǎn)角處最大流速為15.45 m/s。轉(zhuǎn)折角處煙氣的流速越高,流經(jīng)該處的液膜受煙氣剪切力就越大,具有動(dòng)量也就越大。當(dāng)液膜的流動(dòng)慣性力足以克服表面張力對(duì)其膜狀形態(tài)的保持時(shí),已被捕集在葉片上的液滴就會(huì)重新卷入主流煙氣并夾帶出除霧器,降低除霧效率。這解釋了在大葉片間距和高入口流速的交互作用下折流板除霧效率惡化的原因。

    圖6 3種板間距在U=7 m/s時(shí)的速度分布云圖對(duì)比Fig.6 Comparison of velocity contour of the three kinds of plate spacing at U=7 m/s

    3.1.2 板間距對(duì)進(jìn)出口壓降的影響

    圖7示出了不同入口流速、不同板間距下折流板進(jìn)出口壓降的變化曲線。從圖中可知,板間距越大,除霧器進(jìn)出口壓降就越高。其中當(dāng)D=50 mm,除霧器進(jìn)出口壓降最大,介于11.55~55.76 Pa。這主要是由于在葉片板間距增大后,當(dāng)煙氣流經(jīng)轉(zhuǎn)折角發(fā)生偏折時(shí),局部渦旋增強(qiáng),部分區(qū)域產(chǎn)生回流區(qū)所致。從圖6可以明顯看出,相較于大板間距,狹小流道更貼合煙氣主流區(qū)輪廓,使得流場(chǎng)中回流區(qū)域更小,降低了煙氣流動(dòng)過(guò)程的阻力。

    圖7 不同流速和不同板間距下除霧器進(jìn)出口的壓降Fig.7 Pressure drop at the inlet and outlet of the demister under different flow rates and plate spacing

    綜上所述,為了獲得理想的除霧性能,板間距優(yōu)化值可取D=20 mm。除霧效率相對(duì)于優(yōu)化前平均提升了16.38%,壓降平均降低了9.74%。

    3.2 轉(zhuǎn)折角度對(duì)除霧器性能的影響

    為研究葉片轉(zhuǎn)折角對(duì)折流板除霧性能的影響,在其他條件不變的情況下,選取α=46°,42°,40°,38°,34°5種不同轉(zhuǎn)折角情況進(jìn)行比較研究。其中,板間距D=30 mm,直流段長(zhǎng)度H1=38 mm,捕集段長(zhǎng)度H2=49.6 mm。

    3.2.1 轉(zhuǎn)折角度對(duì)除霧效率的影響

    圖8示出了在不同入口流速、不同轉(zhuǎn)折角情況下折流板除霧效率的變化曲線。可以看出,葉片轉(zhuǎn)折角和煙氣入口流速的交互作用對(duì)除霧效率的影響非常明顯。隨著葉片轉(zhuǎn)折角的增大,除霧器臨界速度不斷下降。當(dāng)α=34°,38°時(shí),臨界速度大于或等于7 m/s;當(dāng)α=40°,42°時(shí),臨界速度為6 m/s;當(dāng)α=46°時(shí),臨界速度為5 m/s。在入口流速低于臨界流速情況下,增大轉(zhuǎn)折角能夠有效提高除霧效率。但是,入口流速一旦高于臨界流速,增大轉(zhuǎn)折角則會(huì)加劇二次攜帶現(xiàn)象的發(fā)生。以U=7 m/s為例,當(dāng)α由40°增大至42°和46°時(shí),除霧效率由80.08%分別下降至45.35%和7.71%,極大惡化了折流板的除霧性能。

    圖8 不同流速、不同轉(zhuǎn)折角下除霧器的效率Fig.8 Efficiency of the demister at different flow rates and angles

    圖9為3種不同葉片轉(zhuǎn)折角在入口流速為7 m/s的工況下的流場(chǎng)速度分布云圖對(duì)比。從圖中可以看出,增大葉片轉(zhuǎn)折角增強(qiáng)了主流煙氣流動(dòng)的偏折程度,并且減小了捕集段流域法向的流通面積,提高了局部流動(dòng)速度,這都是有利于折流板對(duì)液滴的捕集。然而,越高的局部流速意味著液膜所受的剪切力越大,具有的流動(dòng)慣性越大,流經(jīng)轉(zhuǎn)折角處分離出液滴的可能性也越大。

    圖9 3種轉(zhuǎn)折角在U=7 m/s時(shí)的速度分布云圖對(duì)比Fig.9 Comparison of velocity contour of the three kinds of angles at U=7 m/s

    3.2.2 轉(zhuǎn)折角度對(duì)進(jìn)出口壓降的影響

    圖10為不同入口流速、不同轉(zhuǎn)折角下折流板進(jìn)出口壓降的變化曲線。計(jì)算結(jié)果顯示,轉(zhuǎn)折角越大,除霧器進(jìn)出口壓降越大,除霧系統(tǒng)所需的能耗越高。這是由于越大的轉(zhuǎn)折角會(huì)使煙氣流場(chǎng)更為曲折,導(dǎo)致流場(chǎng)內(nèi)回流區(qū)加劇,增大局部阻力。從數(shù)值上看,當(dāng)α由40°分別增大至42°和46°時(shí),除霧器進(jìn)出口壓降值平均上升了21.29%和80.90%,壓力損失增加明顯。

    圖10 不同流速、不同轉(zhuǎn)折角下除霧器進(jìn)出口的壓降Fig.10 Pressure drop at the inlet and outlet of the demister under different flow rates and angles

    綜上所述,增大葉片轉(zhuǎn)折角度雖然能在低入口流速下有效提高折流板的除霧效率,但入口流速一旦高于臨界速度,則除霧效率會(huì)急劇降低,難以滿足實(shí)際工程需要。此外,增大葉片轉(zhuǎn)折角會(huì)大幅提高除霧系統(tǒng)所需的能耗。因此,綜合考慮除霧效率的穩(wěn)定性壓降變化,原模型轉(zhuǎn)折角α=40°為最優(yōu)值。

    4 結(jié)論

    本文以折流板除霧器為研究對(duì)象,采用數(shù)值模擬的方法,在傳統(tǒng)歐拉-拉格朗日方法的基礎(chǔ)上耦合EWF模型,構(gòu)建了基于二次攜帶的折流板除霧性能的預(yù)測(cè)模型。從葉片板間距和轉(zhuǎn)折角兩個(gè)方面出發(fā),探究不同流速下折流板除霧性能的變化規(guī)律,得到以下結(jié)論:

    1)將描述液膜生成及運(yùn)動(dòng)特性的EWF模型與傳統(tǒng)歐拉-拉格朗日法進(jìn)行耦合,建立了考慮二次攜帶現(xiàn)象的除霧模型。該模型預(yù)測(cè)曲線變化與實(shí)驗(yàn)結(jié)果一致,模型可行;

    2)折流板除霧過(guò)程中二次攜帶現(xiàn)象的發(fā)生主要是由流道轉(zhuǎn)角處局部煙氣流速過(guò)大造成液膜分離導(dǎo)致。增大葉片板間距和轉(zhuǎn)折角均會(huì)提升該處流速,使除霧器臨界速度下降,加劇液膜分離。

    3)在低于臨界速度時(shí),除霧效率隨板間距增大而減少,隨轉(zhuǎn)折角增大而增大;在高于臨界速度時(shí),除霧效率隨板間距增大而減小,隨轉(zhuǎn)折角增大而降低。

    4)模型板間距優(yōu)化值取20 mm,轉(zhuǎn)折角度取40°時(shí),除霧效率可提升16.38%,壓降降低9.74%。該結(jié)果對(duì)防止除霧器二次攜帶及提高除霧性能具有指導(dǎo)作用。

    猜你喜歡
    流板液膜液滴
    考慮軸彎曲的水潤(rùn)滑軸承液膜建模方法
    高空高速氣流下平板液膜流動(dòng)與破裂規(guī)律
    結(jié)構(gòu)參數(shù)對(duì)交錯(cuò)百葉折流板管殼式換熱器性能影響的研究
    新型錐形孔折流板管殼式換熱器數(shù)值模擬研究
    液膜破裂對(duì)PCCS降膜的影響*
    液滴間相互碰撞融合與破碎的實(shí)驗(yàn)研究
    噴淋液滴在空氣環(huán)境下的運(yùn)動(dòng)特性
    螺旋折流板換熱器在高溫高壓環(huán)境中的應(yīng)用研究
    豎直窄矩形通道內(nèi)彈狀流中液膜特性研究
    氣井多液滴攜液理論模型研究
    河南科技(2013年10期)2013-08-12 05:55:10
    欧美性长视频在线观看| 国产麻豆69| 视频区图区小说| 精品一区二区三区四区五区乱码 | 精品第一国产精品| 日韩av免费高清视频| 国产黄色免费在线视频| 99精国产麻豆久久婷婷| 久久国产精品大桥未久av| 国产爽快片一区二区三区| 久久综合国产亚洲精品| 国产精品人妻久久久影院| 如日韩欧美国产精品一区二区三区| 国产亚洲精品久久久久5区| av视频免费观看在线观看| 丝瓜视频免费看黄片| 国产精品99久久99久久久不卡| 真人做人爱边吃奶动态| 久久久精品94久久精品| 亚洲av男天堂| 大香蕉久久成人网| 叶爱在线成人免费视频播放| 人人妻人人爽人人添夜夜欢视频| tube8黄色片| 人人妻人人爽人人添夜夜欢视频| www.av在线官网国产| 91精品三级在线观看| videosex国产| 欧美激情极品国产一区二区三区| 亚洲国产精品999| 国产不卡av网站在线观看| 欧美人与善性xxx| 亚洲人成网站在线观看播放| 免费少妇av软件| 午夜免费成人在线视频| 国产精品久久久久久精品古装| 老司机在亚洲福利影院| 日本欧美国产在线视频| 天堂俺去俺来也www色官网| 黑人巨大精品欧美一区二区蜜桃| 在线观看一区二区三区激情| 黄片小视频在线播放| 18禁国产床啪视频网站| 日韩精品免费视频一区二区三区| 熟女av电影| 色94色欧美一区二区| 婷婷色麻豆天堂久久| 日本色播在线视频| 一本久久精品| 一级a爱视频在线免费观看| 国产免费视频播放在线视频| 一本综合久久免费| 精品国产国语对白av| 99热全是精品| 丝袜喷水一区| 男人添女人高潮全过程视频| 香蕉丝袜av| 午夜av观看不卡| 一本—道久久a久久精品蜜桃钙片| www.熟女人妻精品国产| 一区二区av电影网| 亚洲欧洲国产日韩| 少妇人妻 视频| 成人三级做爰电影| 亚洲精品国产一区二区精华液| 中国国产av一级| 曰老女人黄片| 日韩中文字幕欧美一区二区 | 成年人黄色毛片网站| 如日韩欧美国产精品一区二区三区| 性少妇av在线| av一本久久久久| 日本av手机在线免费观看| 男男h啪啪无遮挡| 一二三四在线观看免费中文在| 国产成人系列免费观看| 国产伦人伦偷精品视频| 久久ye,这里只有精品| 97精品久久久久久久久久精品| 日韩欧美一区视频在线观看| 在线观看免费日韩欧美大片| 久久国产亚洲av麻豆专区| 在线亚洲精品国产二区图片欧美| 国产欧美亚洲国产| 欧美日韩视频精品一区| 一本久久精品| 国产又色又爽无遮挡免| 国产精品偷伦视频观看了| 成人18禁高潮啪啪吃奶动态图| 国产亚洲av高清不卡| 久久av网站| 欧美日韩视频精品一区| 美国免费a级毛片| 91成人精品电影| 久久人妻福利社区极品人妻图片 | 欧美变态另类bdsm刘玥| 麻豆国产av国片精品| 91精品国产国语对白视频| 亚洲精品国产av蜜桃| 人人妻人人澡人人爽人人夜夜| 色婷婷av一区二区三区视频| 曰老女人黄片| 成年人免费黄色播放视频| 18禁黄网站禁片午夜丰满| 精品高清国产在线一区| 精品视频人人做人人爽| 香蕉国产在线看| 看免费av毛片| 搡老岳熟女国产| 亚洲国产欧美一区二区综合| 丰满人妻熟妇乱又伦精品不卡| 建设人人有责人人尽责人人享有的| 精品少妇黑人巨大在线播放| 亚洲伊人久久精品综合| 久久久久国产精品人妻一区二区| 日韩,欧美,国产一区二区三区| 国产极品粉嫩免费观看在线| 高清欧美精品videossex| 天天影视国产精品| 在线观看免费日韩欧美大片| 国产成人免费无遮挡视频| 中文乱码字字幕精品一区二区三区| 欧美 亚洲 国产 日韩一| 女性被躁到高潮视频| 国产片特级美女逼逼视频| 久久久精品区二区三区| 十八禁高潮呻吟视频| 丝袜喷水一区| 国产成人一区二区在线| www.999成人在线观看| 日韩熟女老妇一区二区性免费视频| 国产免费福利视频在线观看| 亚洲国产精品一区二区三区在线| 老鸭窝网址在线观看| 日本vs欧美在线观看视频| 一级毛片我不卡| 美女午夜性视频免费| 国产成人精品久久久久久| 国产精品欧美亚洲77777| 又紧又爽又黄一区二区| 尾随美女入室| 18禁裸乳无遮挡动漫免费视频| 黄色视频在线播放观看不卡| 亚洲国产日韩一区二区| 自拍欧美九色日韩亚洲蝌蚪91| 女性被躁到高潮视频| 亚洲精品一二三| 美女主播在线视频| 午夜福利影视在线免费观看| h视频一区二区三区| 1024视频免费在线观看| 脱女人内裤的视频| 国产精品麻豆人妻色哟哟久久| 久久天躁狠狠躁夜夜2o2o | 免费高清在线观看日韩| 丰满人妻熟妇乱又伦精品不卡| 蜜桃国产av成人99| 中文字幕人妻丝袜制服| 天堂8中文在线网| 国产在线一区二区三区精| 男女下面插进去视频免费观看| 国产片内射在线| 999精品在线视频| 人人妻人人澡人人爽人人夜夜| 啦啦啦在线观看免费高清www| 亚洲欧美一区二区三区国产| 女性生殖器流出的白浆| 夫妻午夜视频| 日本五十路高清| 成年美女黄网站色视频大全免费| 亚洲av男天堂| 男人爽女人下面视频在线观看| 日韩精品免费视频一区二区三区| 又大又黄又爽视频免费| 亚洲自偷自拍图片 自拍| 一级毛片电影观看| 亚洲 国产 在线| 99国产综合亚洲精品| 欧美国产精品一级二级三级| 韩国精品一区二区三区| 麻豆av在线久日| 久久99热这里只频精品6学生| 在线观看人妻少妇| 人体艺术视频欧美日本| videosex国产| 亚洲 国产 在线| 男人舔女人的私密视频| 女人久久www免费人成看片| 中文欧美无线码| 男女边吃奶边做爰视频| 黑人猛操日本美女一级片| 免费久久久久久久精品成人欧美视频| 亚洲国产中文字幕在线视频| 欧美中文综合在线视频| 亚洲欧洲日产国产| 国产精品三级大全| 欧美国产精品一级二级三级| 搡老乐熟女国产| 国产精品一区二区免费欧美 | 黄色片一级片一级黄色片| 黄色片一级片一级黄色片| 亚洲av成人精品一二三区| 麻豆国产av国片精品| 日本色播在线视频| 女警被强在线播放| 久久久久久亚洲精品国产蜜桃av| 欧美黄色淫秽网站| 久久久国产一区二区| 久久久精品94久久精品| 午夜免费鲁丝| 国产av国产精品国产| 啦啦啦视频在线资源免费观看| 亚洲国产精品一区三区| 精品少妇久久久久久888优播| 亚洲精品久久午夜乱码| 中文字幕高清在线视频| 国产精品久久久久成人av| 亚洲精品一卡2卡三卡4卡5卡 | 精品一区在线观看国产| 日本av手机在线免费观看| 美女福利国产在线| 欧美精品啪啪一区二区三区 | 亚洲欧洲日产国产| 欧美日韩视频高清一区二区三区二| 亚洲欧美日韩高清在线视频 | 一级片免费观看大全| 777米奇影视久久| av电影中文网址| 国产亚洲精品第一综合不卡| 老司机影院毛片| 精品久久蜜臀av无| 在现免费观看毛片| 国产又色又爽无遮挡免| 久久天躁狠狠躁夜夜2o2o | 亚洲欧美色中文字幕在线| 高清黄色对白视频在线免费看| 色综合欧美亚洲国产小说| 王馨瑶露胸无遮挡在线观看| 亚洲少妇的诱惑av| 精品免费久久久久久久清纯 | 热99国产精品久久久久久7| 久久国产精品男人的天堂亚洲| 免费女性裸体啪啪无遮挡网站| 操美女的视频在线观看| 老司机亚洲免费影院| 80岁老熟妇乱子伦牲交| 日韩大码丰满熟妇| 日本av手机在线免费观看| 成人国产一区最新在线观看 | 一本综合久久免费| 在线观看免费视频网站a站| 最新的欧美精品一区二区| 9色porny在线观看| 嫩草影视91久久| 久久人妻福利社区极品人妻图片 | 国产欧美日韩一区二区三 | 一级片免费观看大全| 啦啦啦中文免费视频观看日本| 亚洲国产精品国产精品| 中文字幕人妻丝袜制服| 亚洲综合色网址| 曰老女人黄片| 成年动漫av网址| 亚洲av日韩在线播放| 18禁裸乳无遮挡动漫免费视频| 亚洲国产欧美网| 国产亚洲一区二区精品| 成人午夜精彩视频在线观看| 欧美激情高清一区二区三区| 别揉我奶头~嗯~啊~动态视频 | 欧美精品av麻豆av| 伦理电影免费视频| 我要看黄色一级片免费的| 亚洲精品一区蜜桃| 老司机影院毛片| 日韩伦理黄色片| 欧美日韩视频高清一区二区三区二| 久久久久久人人人人人| 在线看a的网站| 黑人欧美特级aaaaaa片| 亚洲视频免费观看视频| 国产精品 欧美亚洲| 亚洲av成人不卡在线观看播放网 | 国产人伦9x9x在线观看| avwww免费| 日韩一本色道免费dvd| 欧美日韩国产mv在线观看视频| 亚洲国产精品999| 美女扒开内裤让男人捅视频| 欧美av亚洲av综合av国产av| 十分钟在线观看高清视频www| 婷婷成人精品国产| 两性夫妻黄色片| 手机成人av网站| 久久人妻福利社区极品人妻图片 | 久久人妻熟女aⅴ| 午夜福利视频精品| 在线天堂中文资源库| 亚洲精品国产色婷婷电影| 欧美97在线视频| 桃花免费在线播放| 女人高潮潮喷娇喘18禁视频| 日本一区二区免费在线视频| 老司机靠b影院| 男女边摸边吃奶| 亚洲精品日本国产第一区| 精品久久久久久久毛片微露脸 | 亚洲自偷自拍图片 自拍| 性高湖久久久久久久久免费观看| 一级毛片女人18水好多 | 午夜福利视频精品| 女性生殖器流出的白浆| 国产av国产精品国产| 人妻一区二区av| 亚洲国产欧美在线一区| 国产真人三级小视频在线观看| e午夜精品久久久久久久| 欧美成人精品欧美一级黄| 视频在线观看一区二区三区| 一级黄色大片毛片| 国产黄色免费在线视频| 秋霞在线观看毛片| 国产成人影院久久av| 黄色视频在线播放观看不卡| 日韩人妻精品一区2区三区| 人妻人人澡人人爽人人| 亚洲九九香蕉| 99久久99久久久精品蜜桃| 老汉色av国产亚洲站长工具| 我的亚洲天堂| 80岁老熟妇乱子伦牲交| 亚洲精品国产色婷婷电影| 午夜免费观看性视频| 精品欧美一区二区三区在线| 欧美在线黄色| 免费久久久久久久精品成人欧美视频| 一区福利在线观看| 日本猛色少妇xxxxx猛交久久| 亚洲第一av免费看| 大码成人一级视频| 国产日韩欧美在线精品| 大片电影免费在线观看免费| 日韩 欧美 亚洲 中文字幕| 狂野欧美激情性xxxx| 男女午夜视频在线观看| 亚洲精品乱久久久久久| 免费少妇av软件| 日韩中文字幕视频在线看片| 欧美av亚洲av综合av国产av| 亚洲国产欧美在线一区| 一区福利在线观看| 激情五月婷婷亚洲| 宅男免费午夜| 午夜免费鲁丝| 五月开心婷婷网| 丁香六月欧美| 人人妻人人澡人人爽人人夜夜| 免费一级毛片在线播放高清视频 | 黄片小视频在线播放| 国产人伦9x9x在线观看| 一本一本久久a久久精品综合妖精| www.熟女人妻精品国产| 黄片小视频在线播放| 国产男女内射视频| 自线自在国产av| 国产av一区二区精品久久| 男女午夜视频在线观看| 国产福利在线免费观看视频| svipshipincom国产片| 亚洲精品中文字幕在线视频| 国产91精品成人一区二区三区 | 老司机亚洲免费影院| 午夜福利一区二区在线看| 一级毛片 在线播放| 久久久精品区二区三区| 满18在线观看网站| 高清视频免费观看一区二区| 免费在线观看完整版高清| 大话2 男鬼变身卡| 国产激情久久老熟女| 欧美人与性动交α欧美软件| 午夜免费成人在线视频| 视频区欧美日本亚洲| 国产一区亚洲一区在线观看| av国产精品久久久久影院| 天天影视国产精品| 桃花免费在线播放| 亚洲国产最新在线播放| 精品国产一区二区三区四区第35| 少妇粗大呻吟视频| 亚洲,欧美,日韩| 操美女的视频在线观看| 男人舔女人的私密视频| 在线观看免费日韩欧美大片| 男人爽女人下面视频在线观看| 少妇裸体淫交视频免费看高清 | 国产成人欧美在线观看 | 亚洲国产av影院在线观看| 晚上一个人看的免费电影| 91国产中文字幕| 一二三四在线观看免费中文在| cao死你这个sao货| 丰满少妇做爰视频| 国产高清不卡午夜福利| 啦啦啦视频在线资源免费观看| 国产一区二区三区综合在线观看| 青春草亚洲视频在线观看| 久久久精品免费免费高清| 免费久久久久久久精品成人欧美视频| 国产日韩欧美在线精品| 亚洲精品乱久久久久久| 新久久久久国产一级毛片| 亚洲欧洲国产日韩| 秋霞在线观看毛片| 99精国产麻豆久久婷婷| 欧美国产精品一级二级三级| 亚洲av在线观看美女高潮| 久久人妻熟女aⅴ| 黄色一级大片看看| 亚洲欧美成人综合另类久久久| 亚洲av片天天在线观看| 国产免费又黄又爽又色| 人人妻人人澡人人看| 亚洲精品成人av观看孕妇| 午夜福利免费观看在线| 99国产综合亚洲精品| 国产熟女欧美一区二区| 熟女少妇亚洲综合色aaa.| 精品国产超薄肉色丝袜足j| 黄色视频不卡| 伦理电影免费视频| 亚洲少妇的诱惑av| 国产熟女午夜一区二区三区| 亚洲精品成人av观看孕妇| 久久久亚洲精品成人影院| 国产成人免费无遮挡视频| 中文字幕人妻丝袜一区二区| 午夜久久久在线观看| 韩国精品一区二区三区| 国产黄色免费在线视频| 咕卡用的链子| 国产一级毛片在线| 亚洲精品久久成人aⅴ小说| 色精品久久人妻99蜜桃| 丝袜美足系列| 黄片小视频在线播放| 自拍欧美九色日韩亚洲蝌蚪91| 男人舔女人的私密视频| 午夜av观看不卡| 日本欧美国产在线视频| 夫妻性生交免费视频一级片| av国产精品久久久久影院| 另类亚洲欧美激情| 亚洲成国产人片在线观看| 伊人久久大香线蕉亚洲五| 国产一区二区三区av在线| 国产精品免费视频内射| 大型av网站在线播放| 99国产精品一区二区三区| 久久精品人人爽人人爽视色| 久久精品久久精品一区二区三区| 91精品三级在线观看| 50天的宝宝边吃奶边哭怎么回事| 蜜桃在线观看..| 一级毛片 在线播放| 国产一级毛片在线| 亚洲国产欧美一区二区综合| 欧美性长视频在线观看| 女性生殖器流出的白浆| 亚洲精品久久午夜乱码| 熟女av电影| www.精华液| 日本欧美视频一区| 99热网站在线观看| www.自偷自拍.com| 两性夫妻黄色片| 亚洲精品第二区| 国产精品av久久久久免费| av在线播放精品| 岛国毛片在线播放| 丁香六月欧美| 亚洲精品美女久久av网站| 中文字幕色久视频| 亚洲精品国产av蜜桃| av福利片在线| 老熟女久久久| 七月丁香在线播放| 欧美+亚洲+日韩+国产| 亚洲熟女精品中文字幕| 国产亚洲一区二区精品| 国产女主播在线喷水免费视频网站| 9191精品国产免费久久| 国产精品九九99| 好男人视频免费观看在线| 国产av精品麻豆| 人妻人人澡人人爽人人| 女人久久www免费人成看片| 黄色怎么调成土黄色| 亚洲精品美女久久av网站| 男男h啪啪无遮挡| 色94色欧美一区二区| 脱女人内裤的视频| 亚洲精品美女久久久久99蜜臀 | 午夜免费成人在线视频| 18禁黄网站禁片午夜丰满| 男人舔女人的私密视频| 亚洲欧美成人综合另类久久久| 99久久人妻综合| 精品熟女少妇八av免费久了| 啦啦啦啦在线视频资源| 婷婷色综合大香蕉| 欧美成人精品欧美一级黄| 久久久亚洲精品成人影院| 欧美精品亚洲一区二区| 纵有疾风起免费观看全集完整版| 一本大道久久a久久精品| 国产女主播在线喷水免费视频网站| 在线 av 中文字幕| 亚洲久久久国产精品| √禁漫天堂资源中文www| 国产精品久久久av美女十八| 色婷婷av一区二区三区视频| 婷婷色av中文字幕| 最新的欧美精品一区二区| 免费看av在线观看网站| 女警被强在线播放| 伊人久久大香线蕉亚洲五| 欧美日韩福利视频一区二区| 成年av动漫网址| 日本wwww免费看| 国产高清视频在线播放一区 | 波多野结衣一区麻豆| 一二三四在线观看免费中文在| 叶爱在线成人免费视频播放| 高潮久久久久久久久久久不卡| 人妻 亚洲 视频| 欧美日韩亚洲国产一区二区在线观看 | 午夜视频精品福利| 香蕉丝袜av| 欧美 亚洲 国产 日韩一| 亚洲九九香蕉| 亚洲国产欧美网| 成人免费观看视频高清| 丝袜美足系列| 黄频高清免费视频| 夜夜骑夜夜射夜夜干| 精品免费久久久久久久清纯 | 久久久久久人人人人人| 热re99久久国产66热| 亚洲人成77777在线视频| 五月天丁香电影| 中国美女看黄片| 九色亚洲精品在线播放| 男女免费视频国产| a级毛片黄视频| 欧美在线一区亚洲| 丝袜脚勾引网站| 亚洲黑人精品在线| 老司机深夜福利视频在线观看 | 亚洲中文av在线| 丝袜美足系列| 丝袜美腿诱惑在线| 我要看黄色一级片免费的| a级毛片在线看网站| 狂野欧美激情性bbbbbb| 久久精品国产亚洲av高清一级| 日韩免费高清中文字幕av| 精品国产乱码久久久久久男人| 老司机影院成人| 肉色欧美久久久久久久蜜桃| 老司机深夜福利视频在线观看 | 91国产中文字幕| 亚洲 国产 在线| 又粗又硬又长又爽又黄的视频| 精品一品国产午夜福利视频| 夫妻午夜视频| 亚洲第一青青草原| 99久久99久久久精品蜜桃| 国产免费福利视频在线观看| 久久久久网色| 大型av网站在线播放| 国产日韩欧美视频二区| 精品一区二区三区四区五区乱码 | 日韩熟女老妇一区二区性免费视频| 丝袜美足系列| 日本a在线网址| 亚洲熟女毛片儿| 久久久精品国产亚洲av高清涩受| 精品亚洲乱码少妇综合久久| 两人在一起打扑克的视频| 日韩欧美一区视频在线观看| xxxhd国产人妻xxx| 国产av精品麻豆| 欧美性长视频在线观看| 久久久久久久久久久久大奶| 悠悠久久av| 777米奇影视久久| 男女床上黄色一级片免费看| 丰满少妇做爰视频| 肉色欧美久久久久久久蜜桃| 在线 av 中文字幕| 久久九九热精品免费| 99久久人妻综合| 精品久久蜜臀av无| 久久毛片免费看一区二区三区| 日本av手机在线免费观看| 中文字幕另类日韩欧美亚洲嫩草| 精品福利观看| 欧美人与性动交α欧美精品济南到| 老熟女久久久| 国产精品欧美亚洲77777| av有码第一页| a级毛片黄视频| 日韩伦理黄色片| 国产真人三级小视频在线观看|