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

    立式水力碎漿機(jī)內(nèi)部漿料流場數(shù)值模擬及其新型槽體結(jié)構(gòu)的研究(一)

    2016-11-16 05:53:27沙九龍
    中國造紙學(xué)報 2016年2期
    關(guān)鍵詞:槽體形槽漿料

    竇 靖 張 放 沙九龍 張 輝

    (南京林業(yè)大學(xué)江蘇省制漿造紙科學(xué)與技術(shù)重點(diǎn)實(shí)驗(yàn)室,江蘇南京,210037)

    ?

    立式水力碎漿機(jī)內(nèi)部漿料流場數(shù)值模擬及其新型槽體結(jié)構(gòu)的研究(一)

    竇靖張放沙九龍張輝*

    (南京林業(yè)大學(xué)江蘇省制漿造紙科學(xué)與技術(shù)重點(diǎn)實(shí)驗(yàn)室,江蘇南京,210037)

    以高效、節(jié)能、運(yùn)行平衡安全為主要依據(jù),在分析國內(nèi)外現(xiàn)有專利和應(yīng)用中的立式水力碎漿機(jī)槽體結(jié)構(gòu)型式基礎(chǔ)上,運(yùn)用FLUENT軟件對O形和D形2類槽體結(jié)構(gòu)中漿濃5%的漿料流場數(shù)值模擬,得到兩者內(nèi)部流場的壓力、流速和湍流強(qiáng)度的分布并進(jìn)行對比分析。結(jié)果表明,傳統(tǒng)O形和D形槽體結(jié)構(gòu)內(nèi)部流場均為垂直環(huán)流和水平旋流合成的漩渦流;傳統(tǒng)O形內(nèi)部流場對稱、運(yùn)行平衡穩(wěn)定,但其高湍流強(qiáng)度區(qū)域局限于下部轉(zhuǎn)子旋轉(zhuǎn)區(qū)域附近,上部和下部流場湍流強(qiáng)度相差1159%,且頂部邊緣處漿料圓周切向速度為0而缺乏碎解動力;D形槽體中心低湍流區(qū)域半徑縮小到O形的13.0%,湍流強(qiáng)度明顯增強(qiáng),但漿料徑向撞擊槽壁致使比O形多造成22.7%的動能浪費(fèi),且內(nèi)部流場不對稱使設(shè)備運(yùn)行缺乏平衡穩(wěn)定性。

    O形立式水力碎漿機(jī);D形立式水力碎漿機(jī);FLUENT軟件;數(shù)值模擬;漿料流場分布

    隨著經(jīng)濟(jì)、科技和社會的不斷發(fā)展,對包裝用紙或紙板的需求不斷增長;加之造紙用纖維資源缺乏,同時倡導(dǎo)節(jié)能、降本、少用化學(xué)品等,致使廢紙回收再利用得到高度重視[1]。碎漿設(shè)備是廢紙漿生產(chǎn)的第一道重要設(shè)備,其基本功能是將漿板或廢紙等碎解成纖維漿料,即使原先交織成紙的纖維最大限度地離解成單根纖維,同時又盡可能地保持纖維的原有形態(tài)和強(qiáng)度[2]。

    目前,廢紙制漿或漿板碎解用碎漿設(shè)備大多采用立式水力碎漿機(jī)。其結(jié)構(gòu)相對簡單、占地面積不等,主要組成部分有槽體、轉(zhuǎn)子、篩板、傳動和支承機(jī)構(gòu)等,其中,槽體和轉(zhuǎn)子結(jié)構(gòu)形式是影響碎漿效率和質(zhì)量、能耗和機(jī)械運(yùn)行穩(wěn)定性的關(guān)鍵因素[3]。本研究主要涉及槽體結(jié)構(gòu)形式,無論是何種形式的槽體,其內(nèi)部的流場均是較復(fù)雜的漩渦流[4-5]。近10多年來,人們越來越多地利用CFD(Computational Fluid Dynamic)技術(shù)對漩渦流場進(jìn)行數(shù)值模擬研究[6-12],但是用于立式水力碎漿機(jī)內(nèi)部流場的數(shù)值模擬未見報道。

    本研究基于國內(nèi)外現(xiàn)有專利技術(shù)、應(yīng)用中的立式水力碎漿機(jī)和碎漿原理分類槽體結(jié)構(gòu)型式,然后選取最具代表性的O形和D形槽體結(jié)構(gòu)為研究對象;基于CFD中的FLUENT軟件對水力碎漿機(jī)中漿濃為5%的紙漿懸浮液流場進(jìn)行了數(shù)值模擬與分析,并對獲得的O形和D形槽體內(nèi)部流場的壓力、流速和湍流強(qiáng)度的分布進(jìn)行對比分析,得出現(xiàn)有立式水力碎漿機(jī)的O形和D形槽體結(jié)構(gòu)在碎漿過程中效率、能耗、安全穩(wěn)定性等方面的差異。

    1 國內(nèi)外現(xiàn)有立式水力碎漿機(jī)槽體結(jié)構(gòu)型式分析

    1.1立式水力碎漿機(jī)槽體結(jié)構(gòu)型式分類

    現(xiàn)有的國內(nèi)外立式水力碎漿機(jī)槽體結(jié)構(gòu)型式大致分為以下4類:傳統(tǒng)O形槽體;在槽體內(nèi)壁設(shè)置垂直向△形擋板、助升推式導(dǎo)流板;D形、σ形、G形槽體;轉(zhuǎn)子位于筒形槽的偏心位置。從內(nèi)部流場特征的本質(zhì)角度,將立式水力碎漿機(jī)上述4類槽體結(jié)構(gòu)型式進(jìn)一步簡化歸納,將設(shè)有擋板和導(dǎo)流板的對稱結(jié)構(gòu)槽體縮歸為O形槽體,將σ形、G形、偏心轉(zhuǎn)子的非對稱結(jié)構(gòu)槽體縮歸為D形槽體。所以,O形槽體和D形槽體為目前最具代表性的2類槽體結(jié)構(gòu)型式。

    1.2國內(nèi)外主要立式水力碎漿機(jī)類型及其特征

    1.2.1圓柱形槽體

    (1)經(jīng)典圓柱形槽體傳統(tǒng)經(jīng)典碎漿機(jī)槽體標(biāo)準(zhǔn)型為圓柱形,也稱O形。其槽體結(jié)構(gòu)具有較大的容積,內(nèi)部物料的碎解主要依賴轉(zhuǎn)子的轉(zhuǎn)動,轉(zhuǎn)子轉(zhuǎn)動輸入的能量主要傳送給水,通過水流非常強(qiáng)烈的攪動、剪切,以及紙料之間、紙料與槽體或轉(zhuǎn)子間相互摩擦、碰撞、撕拉和搓揉等機(jī)械作用,快速完成碎漿功能[13]。

    圖3 不同槽體結(jié)構(gòu)的示意圖

    (2)在圓柱形槽體內(nèi)增設(shè)△形擋板日本相川公司在槽體設(shè)計上采用底部兩角結(jié)構(gòu),垂直設(shè)置△形擋板(見圖1a)[14],使放入的紙料很快形成激烈的回流,并且在最短距離內(nèi)回到轉(zhuǎn)子上。G.E.Brown JR等又增設(shè)了2個擋板,即槽體內(nèi)部共計3個擋板、兩兩呈120°焊在槽壁上[15]。維美德公司在此基礎(chǔ)上進(jìn)一步改進(jìn),在槽體內(nèi)共設(shè)4個△形擋板[16],以改善紙料路徑(見圖1b)。在槽體內(nèi)壁設(shè)置垂直向△形擋板,可以阻撓紙料的水平圓周向流動,進(jìn)而縮短回流中心區(qū)的時間。

    圖1 設(shè)△形擋板圓柱形槽體結(jié)構(gòu)示意圖

    (3)在圓柱形槽體內(nèi)增設(shè)助升推式導(dǎo)流板在1970年Edward H.Cumpston等提出了在槽壁內(nèi)設(shè)置流線形導(dǎo)流板(見圖2)[17],該流線形導(dǎo)流板葉片置于槽體內(nèi)壁工作液面以下。導(dǎo)流板呈45°順著紙料的漩渦方向排布,幫助紙料提升,從而加速紙料的周期循環(huán)運(yùn)動。導(dǎo)流板分布于不同的高度,以滿足不同工作液面的需求。

    圖2 45°流線形助升推式導(dǎo)流板圓柱形 槽體結(jié)構(gòu)示意圖

    1.2.2D形、σ形、G形槽體

    Thermo-Black Clawson公司的D形連續(xù)式水力碎漿機(jī)是對圓形槽體的改進(jìn)(見圖3a)。D形可使物料水平向旋流為非穩(wěn)定的D形湍流,一定程度上產(chǎn)生松散物料塊,使得轉(zhuǎn)子和紙料的接觸更迅速、頻率更高。對使用的圓形槽水力碎漿機(jī),可用擋板將它改造成D形槽,以提高生產(chǎn)能力,但擋板的設(shè)置應(yīng)根據(jù)轉(zhuǎn)子尺寸大小、轉(zhuǎn)速及動力配置等因素進(jìn)行精心設(shè)計,以便達(dá)到最高生產(chǎn)能力。

    原貝洛依特公司設(shè)計的碎漿機(jī)槽體為σ形(見圖3b),目的是改變紙料路徑。Helix碎漿機(jī)為G形槽體(見圖3c),其轉(zhuǎn)子回轉(zhuǎn)形成平滑的水平回流,不僅提高了離解性能,還能降低動力消耗;而且,在垂直方向產(chǎn)生偏心的縱向回流(稱為“G式回流”);當(dāng)廢紙?jiān)贤度氩垠w后就會被迅速有力地吞進(jìn),加速了離解[18]。

    1.2.3偏心轉(zhuǎn)子

    福伊特公司推出的Intensa碎漿機(jī)的轉(zhuǎn)子位于圓柱形槽的偏心位置[19-20](見圖4)。Intensa碎漿機(jī)具有以下特點(diǎn):一個帶有偏心轉(zhuǎn)子的槽;雙錐形槽底;經(jīng)過優(yōu)化的擋板排列,只有1個擋板;經(jīng)過優(yōu)化的轉(zhuǎn)子,改進(jìn)碎漿效果。該碎漿機(jī)效率高于常規(guī)碎漿機(jī),槽內(nèi)紙料流的湍流程度更強(qiáng),可以更有效地實(shí)現(xiàn)纖維分離作用。

    圖4 Intensa碎漿機(jī)槽體內(nèi)轉(zhuǎn)際偏心示意圖

    2 O形和D形槽體結(jié)構(gòu)內(nèi)部流場數(shù)值模擬

    2.1模型的建立

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

    漿料流動假設(shè)為不可壓縮黏性流動,流體被視為連續(xù)介質(zhì),可忽略纖維間的相互作用力。在以上假設(shè)的基礎(chǔ)上,水力碎漿機(jī)中的流體流動過程可以用連續(xù)性方程和Navier-Stokes方程(以下簡稱“N-S方程”)來表述[21]。

    連續(xù)性方程:

    (1)

    N-S方程:

    (2)

    根據(jù)流體流動理論,紙漿懸浮液在水力碎漿機(jī)中的流動是三維湍流問題,由于湍流的復(fù)雜性,通常需要借助合適的湍流模型。在工程中采用較多、較成熟的湍流模型是可實(shí)現(xiàn)性k-ε雙方程模型,其是通過對標(biāo)準(zhǔn)k-ε模型修正而得。由于標(biāo)準(zhǔn)k-ε模型在回流、大曲率過流表面和強(qiáng)旋度等情況下不能很好地預(yù)測湍流特性,因此,在偏離最優(yōu)工況計算時該模型常常失效。由此對標(biāo)準(zhǔn)k-ε模型進(jìn)行修正,最后得到可實(shí)現(xiàn)性k-ε模型,方程如下。

    湍流動能k方程:

    (3)

    湍流耗散率ε方程:

    (4)

    在上述方程中,Gk為平均速度梯度引起的湍動能,Gb為浮力引起的湍動能,YM代表可壓速湍流脈動膨脹對總耗散率的影響。C2和C1ε是常數(shù),分別為C2=1.9、C1ε=1.44;σk和σε分別是湍動能及其耗散率的湍流普特朗數(shù),在FLUENT中,作為常數(shù),分別為σk=1.0、σε=1.2[22-23]。

    2.1.2物理模型

    根據(jù)O形槽體結(jié)構(gòu)的通用尺寸規(guī)范,其體積大約為7~100 m3;選擇大容積的80 m3槽體作為研究對象,其結(jié)構(gòu)尺寸為槽體直徑D1=4.6 m、槽體高度H1=4.6 m、轉(zhuǎn)子直徑d=1.65 m和轉(zhuǎn)子厚度h=0.2 m(見圖5a);轉(zhuǎn)子為應(yīng)用廣泛的伏克斯轉(zhuǎn)子[24-25]。同理,D形槽體體積也為80 m3;其他結(jié)構(gòu)尺寸為槽體最大直徑D2=5 m、槽體高度H2=5 m、D形截面的直線部分長L2=4.4 m和轉(zhuǎn)子中心距D形截面圓心X2=0.625 m(見圖5b);轉(zhuǎn)子結(jié)構(gòu)大小與O形相同。本研究模型均是間歇式敞口水力碎漿機(jī),靜止時漿料液面位于槽體高度的85%。利用Unigraphics NX(UG)軟件分別對O形和D形槽體內(nèi)部流場建立模型。

    圖5 不同槽體尺寸示意圖

    將建好的模型分別導(dǎo)入到ICEM CFD軟件中,對三維實(shí)體模型進(jìn)行非結(jié)構(gòu)網(wǎng)格劃分。由于后期數(shù)值模擬設(shè)置的需要,對槽體結(jié)構(gòu)分塊進(jìn)行網(wǎng)格化,再把網(wǎng)格化好的部件組裝起來(見圖6)。由于立式水力碎漿機(jī)模型較大,整個計算域中除了轉(zhuǎn)子旋轉(zhuǎn)域使用了四面體網(wǎng)格,其余區(qū)域均采用六面體網(wǎng)格,以降低網(wǎng)格數(shù)量,順利進(jìn)行計算機(jī)數(shù)值模擬[26]。

    圖6 不同槽體網(wǎng)格模型示意圖

    2.1.3流體物性及邊界條件的設(shè)置

    在用FLUENT軟件進(jìn)行模擬計算之前,需要對漿料的密度、黏度以及頂層鋪設(shè)的空氣域性質(zhì)等基本參數(shù)進(jìn)行設(shè)置。

    (1)漿料的密度和黏度本研究采用漿濃5%的OCC漿[27]。查閱文獻(xiàn)[28]得到漿料密度ρ及其濃度Cm的關(guān)系滿足以下關(guān)系式:

    (5)

    式中,ρf為絕干纖維密度,1520kg/m3;ρw為水的密度,1000kg/m3;Cm=5%。由式(5)可得漿濃5%OCC漿的密度約為1017kg/m3。

    通過實(shí)驗(yàn)室RST-SST軟固體測試儀測得的漿料黏度μ與剪切速率γ的關(guān)系曲線如圖7所示。根據(jù)源數(shù)據(jù)曲線擬合可得方程μ=32.878γ-0.6193。查閱文獻(xiàn)[29]可得立式水力碎漿機(jī)中轉(zhuǎn)子剪切速率γ與轉(zhuǎn)速N的關(guān)系為γ=KsN(Ks=10,γ的單位為1/s,N的單位為r/s)。綜合上述公式,當(dāng)選擇轉(zhuǎn)子轉(zhuǎn)速N=208r/min[30]時,漿料黏度大約為3.658Pa·s。

    (2)頂層鋪設(shè)空氣域由于設(shè)備是敞口型,因此,需對模型頂部設(shè)置空氣域,空氣的性質(zhì)運(yùn)用FLUENT自帶數(shù)據(jù)庫中空氣的設(shè)置進(jìn)行設(shè)定。

    (3)其他參數(shù)紙漿懸浮液為不可壓縮的非牛頓流體,求解器的類型設(shè)為基于壓力;考慮到紙漿纖維與水完全混合,采用可實(shí)現(xiàn)性k-ε湍流模型,算法選擇SIMPLE進(jìn)行求解,時間設(shè)為穩(wěn)態(tài);考慮重力因素,重力大小近似取9.8N/kg;轉(zhuǎn)子轉(zhuǎn)速設(shè)置為22rad/s;設(shè)備頂部敞口,采用壓力出口的邊界條件,壓力為0.1MPa。

    圖7 漿料黏度與剪切速率關(guān)系曲線圖

    2.2模擬結(jié)果與分析

    在對立式水力碎漿機(jī)進(jìn)行三維建模時,選取槽體底部轉(zhuǎn)子中心為坐標(biāo)原點(diǎn)、槽體高度方向?yàn)閆軸方向作為標(biāo)準(zhǔn)坐標(biāo)系。通過FLUENT軟件對網(wǎng)格化的槽體結(jié)構(gòu)模型進(jìn)行模擬計算后,為了更直觀地體現(xiàn)其內(nèi)部流場的特征,在槽體結(jié)構(gòu)內(nèi)部流場模型選取具有代表性的X-Z平面、Z=0平面和Z=0.50、Z=1.25、Z=2.25、Z=3.25、Z=4.00平面,其中,Z=0.50平面位于轉(zhuǎn)子上方附近平面、Z=4.00平面位于漿料域液面下方附近,同時選取了X-Z平面與Z=0.50、Z=1.25、Z=2.25、Z=3.25、Z=4.00平面的交線,分別命名為Line-1、Line-2、Line-3、Line- 4、Line-5,通過這7個平面和5條交線來描述槽體結(jié)構(gòu)內(nèi)部流場的壓力、速度和湍流強(qiáng)度等特征。

    2.2.1計算殘差

    O形和D形立式水力碎漿機(jī)內(nèi)部流場的殘差監(jiān)測曲線如圖8所示。整個O形和D形水力碎漿機(jī)流場在分別經(jīng)過2059次和4009次迭代計算后,各方程計算結(jié)果殘差都小于所設(shè)置的1×10-3,此時,全局的質(zhì)量、動量能量和標(biāo)量達(dá)到了平衡,所有離散的守恒方程在所有單元中滿足指定的誤差,隨計算不再改變,達(dá)到了收斂狀態(tài)。另外,殘差收斂也驗(yàn)證了網(wǎng)格劃分和湍流模型設(shè)置是合理的。

    圖8 不同槽體的水力碎漿機(jī)內(nèi)部流場殘差監(jiān)測曲線圖

    2.2.2O形立式水力碎漿機(jī)內(nèi)部流場分析

    (1)壓力分布O形槽體結(jié)構(gòu)內(nèi)部流場X-Z平面上的壓力分布如圖9a所示。從圖9a可以看出,從底部到頂端壓力逐漸降低。底部轉(zhuǎn)子葉片附近壓力最大,其次是底部靠近壁面的區(qū)域。漿料在轉(zhuǎn)子葉片處受到的最大壓力為400 Pa左右。

    Z=0平面上的壓力分布如圖9b所示。從圖9b可以看出,從中心到壁面處壓力先是小幅度逐漸降低,在靠近壁面處壓力急劇升高。在壓力最大的轉(zhuǎn)子區(qū)域,轉(zhuǎn)子前端壓力最大,表明漿料在此處受到了最大程度的碎解。

    圖9 O形槽體內(nèi)部流場壓力分布

    O形槽體結(jié)構(gòu)內(nèi)部流場中5條交線沿X軸方向的壓力分布散點(diǎn)圖如圖10(位置為0時位于槽體中心軸上,沿著X軸的負(fù)半軸和正半軸依次向相反方向的壁面靠攏)所示。從圖10可以很明顯地得到5條交線上的壓力分布情況:轉(zhuǎn)子中心區(qū)域(X=0)壓力接近0,沿著X軸方向壓力逐漸降低且為負(fù)壓,在轉(zhuǎn)子葉片外沿處(X=±1)壓力達(dá)到最低值,之后壓力開始逐漸升高且由負(fù)壓逐步轉(zhuǎn)變?yōu)檎龎?在接近壁面處(X=±2)壓力急劇上升,最終在壁面處(X=±2.3)壓力達(dá)到最大值。另外,5條交線中從Line-1至Line-5的最大壓力差逐漸降低,這表明從流場底部到頂部促使?jié){料在水平方向上沿徑向流動的壓力差呈逐漸降低的趨勢。其中,漿料在靠近轉(zhuǎn)子區(qū)域(Line-1)的壁面處受到的最大壓力為120 Pa左右,最大壓力差為140 Pa左右。

    圖10 O形槽體內(nèi)部流場5條交線沿 X軸方向的壓力分布散點(diǎn)圖

    (2)速度分布速度矢量圖能夠直觀表征流體的運(yùn)動方向、速度大小和總體趨勢。O形槽體結(jié)構(gòu)內(nèi)部流場X-Z平面上的速度矢量圖如圖11a所示。由于槽體底部轉(zhuǎn)子高速旋轉(zhuǎn),被轉(zhuǎn)子接觸的漿料由底部中心區(qū)甩向周邊區(qū),然后沿槽壁自下向上流動,再從上部的周邊區(qū)流向中心區(qū),自中心區(qū)下沉至轉(zhuǎn)子,在垂直截面上形成循環(huán)流。底部轉(zhuǎn)子旋轉(zhuǎn)區(qū)域和壁面處漿料總體速度較快,同時在同一水平面上壁面處的速度要大于轉(zhuǎn)子中心區(qū)域的速度,流場頂部速度較中下部而言相對較慢。底部中心轉(zhuǎn)子區(qū)域速度最高達(dá)到15 m/s 左右。

    圖11b從左至右依次為Z=0.50、Z=1.25、Z=2.25、Z=3.25、Z=4.00 5個截面上的速度矢量圖。由此可以看出,槽體內(nèi)上中下各水平面的漿料形成水平向旋轉(zhuǎn)流,旋流方向與轉(zhuǎn)子旋轉(zhuǎn)方向一致。

    圖11 O形槽體內(nèi)部流場速度矢量圖

    圖12 O形槽體流場5條交線不同方向的速度分布示意圖

    軸向速度是反映漿料垂直方向運(yùn)動情況的一個重要依據(jù)。O形槽體結(jié)構(gòu)內(nèi)部流場中5條交線的軸向速度分布情況如圖12a所示。從圖12a可以看到,靠近壁面的漿料軸向速度為正,即在垂直方向上呈現(xiàn)爬升運(yùn)動;靠近轉(zhuǎn)子中心軸的漿料軸向速度為負(fù),即在垂直方向上做下沉運(yùn)動。這與圖11的速度矢量圖規(guī)律一致。

    徑向速度是反映漿料水平方向上徑向運(yùn)動情況的一個重要依據(jù)。5條交線的徑向速度分布散點(diǎn)圖如圖12b所示。Line-3、Line- 4、Line-5與Line-1和Line-2的徑向速度方向相反,即流場下部漿料由中心軸流至壁面,而流場上部漿料則由壁面流向中心軸。

    切向速度是反映漿料水平方向上圓周運(yùn)動情況的一個重要依據(jù)。5條交線的切向速度分布散點(diǎn)圖如圖13所示。從中心軸向壁面處,切向速度先逐漸增大再迅速減小。遠(yuǎn)離轉(zhuǎn)子的漿料(Line-5)在壁面處的切向速度最小值甚至達(dá)到0。因此,流場頂部靠近壁面處的漿料由于動能損失形成的旋流運(yùn)動相對動力不足,向心力顯著降低,但可以通過對槽體結(jié)構(gòu)進(jìn)行優(yōu)化設(shè)計來增強(qiáng)流場頂部壁面附近漿料的動力。

    圖13 O形槽體流場5條交線沿X軸的 切向速度分布情況示意圖

    (3)湍流強(qiáng)度分布立式水力碎漿機(jī)的工作原理是通過水力、轉(zhuǎn)子的直接作用以及紙料之間摩擦作用增強(qiáng)漿料的湍動,從而快速地完成碎漿功能。因此,湍流強(qiáng)度可以作為衡量立式水力碎漿機(jī)碎解效果的直觀指標(biāo)。

    圖14 O形槽體內(nèi)部流場湍流強(qiáng)度分布圖

    O形立式水力碎漿機(jī)內(nèi)部流場X-Z特征平面的湍流強(qiáng)度分布如圖14a所示。圖14b為5個截面上的湍流強(qiáng)度分布圖。由這2個圖可以得到,流場底部轉(zhuǎn)子葉片處湍流強(qiáng)度最高,大約為1180%;流場頂部中心區(qū)域湍流強(qiáng)度最低,大約為21%;同一水平面上壁面處的湍流強(qiáng)度高于中心軸區(qū)域;流場頂部壁面處湍流強(qiáng)度很低,80%左右,低于鄰近區(qū)域的湍流強(qiáng)度,這與上述速度分布情況一致。

    2.2.3D形立式水力碎漿機(jī)內(nèi)部流場分析

    圖15 D形槽體內(nèi)部流場壓力分布圖

    (1)壓力分布D形槽體結(jié)構(gòu)內(nèi)部流場X-Z平面和Z=0平面上的壓力分布如圖15所示。由圖15可以得到,漿料在轉(zhuǎn)子葉片處受到的最大壓力為350 kPa左右,整個流場中壓力最小值為3.5 kPa左右,流場各處壓力均大于O形立式水力碎漿機(jī)內(nèi)部流場各處壓力。

    D形槽體結(jié)構(gòu)內(nèi)部流場中5條交線沿X軸方向的壓力分布散點(diǎn)圖如圖16(X軸負(fù)半軸是靠近D形槽體垂直壁面端)所示。與O形槽體相似的是,D形槽體在壁面處壓力達(dá)到最大值;不同的是,D形槽體結(jié)構(gòu)內(nèi)部流場壓力分布并不以轉(zhuǎn)子中心軸為對稱。

    圖16 D形槽體內(nèi)部流場5條交線沿 X軸方向的壓力分布散點(diǎn)圖

    圖17 D形槽體內(nèi)部流場速度矢量圖

    (2)速度分布D形槽體結(jié)構(gòu)內(nèi)部流場X-Z平面上的速度矢量圖如圖17a所示;圖17b為Z=4.00平面上的速度矢量圖,代表水平面上的速度分布。與O形內(nèi)部流場相比,相同的是,D形內(nèi)部流場垂直面上是循環(huán)流、水平面上是旋流;不同的是,D形槽體結(jié)構(gòu)內(nèi)部流場水平面上的旋流中心不是圓形,而是與槽體相似的D形;由Z=4.00平面上可以較明顯地看出,其旋流中心類似雙旋流,D形內(nèi)部垂直截面上漿料速度矢量方向指向中心軸,向中心軸附近聚攏。

    圖18 D形槽體內(nèi)部流場5條交線 沿X軸方向速度分布散點(diǎn)圖

    D形槽體結(jié)構(gòu)內(nèi)部流場中5條交線的速度分布散點(diǎn)圖如圖18(X軸負(fù)半軸是靠近D形槽體垂直壁面端)所示。同樣,與O形槽體相比,D形槽體內(nèi)部流場速度分布并不以轉(zhuǎn)子中心軸為對稱。另外,Line-1、Line-2和Line-3垂直壁面附近的速度要小于弧形壁面附近的速度。通過分別計算漿料在轉(zhuǎn)子附近最大動能與其在靠近垂直壁面和弧形壁面處動能的差值,可得垂直壁面平均動能損失比弧形壁面平均動能損失高22.7%左右,這是由于漿料撞擊垂直壁面造成能量損失所致?;⌒伪诿骓敹薒ine- 4和Line-5表現(xiàn)出速度下降的趨勢,而垂直壁面頂端Line- 4和Line-5則沒有明顯速度下降現(xiàn)象,這是由于弧形壁面頂端漿料動能損失、向心力降低致使其形成的旋流運(yùn)動相對動力不足,同前述O形槽體。

    圖19 D形槽體內(nèi)部流場X-Z平面湍流強(qiáng)度分布圖

    (3)湍流強(qiáng)度分布D形槽體結(jié)構(gòu)內(nèi)部流場X-Z特征平面的湍流強(qiáng)度分布如圖19(X軸負(fù)半軸是靠近D形槽體垂直壁面端)所示。從圖19可見,垂直壁面頂部邊緣湍流強(qiáng)度減弱的現(xiàn)象較弧形壁面顯著減少,湍流增強(qiáng)。

    D形槽體結(jié)構(gòu)是在O形槽體結(jié)構(gòu)的基礎(chǔ)上演變而來的。圖20對比了D形和O形槽體結(jié)構(gòu)內(nèi)部流場中5條交線上的湍流強(qiáng)度沿X軸的分布情況。從圖20可以很明顯地觀察到,D形結(jié)構(gòu)較O形結(jié)構(gòu)中心低湍流區(qū)域半徑顯著縮小,平均縮小到O形的13.0%左右;這可以有效增強(qiáng)整個流場的平均湍流強(qiáng)度,從而提高碎解效率。

    圖20 不同槽體內(nèi)部流場5條交線的湍流強(qiáng)度對比圖

    3 結(jié) 論

    3.1立式水力碎漿機(jī)槽體內(nèi)部流場是由垂直面和水平面上的流動合成的復(fù)合流。在垂直面上呈現(xiàn)下部由轉(zhuǎn)子到壁面、上部由壁面流向中心軸的環(huán)流;在水平面上則是與轉(zhuǎn)子旋轉(zhuǎn)方向一致的旋流,進(jìn)而使總體流場呈現(xiàn)旋渦流。

    3.2O形槽體內(nèi)部形成對稱的漩渦流流場,設(shè)備運(yùn)行平穩(wěn)且相對安全穩(wěn)定;但是漿料的碎解基本在轉(zhuǎn)子旋轉(zhuǎn)區(qū)域內(nèi)完成,下部區(qū)域最大湍流強(qiáng)度在1180%左右,而上部流場為21%,與上部流場相差1159%,流場上部基本沒有湍動發(fā)生。同時,位于流場頂部邊緣的漿料由于動能損失形成的旋流運(yùn)動相對動力不足,切向速度接近0。

    3.3與O形槽體內(nèi)部流場相比,D形槽體內(nèi)部流場湍流強(qiáng)度明顯增強(qiáng),中心低湍流區(qū)域半徑縮小到O形的13%左右,這是D形改變漿料運(yùn)行路徑致使轉(zhuǎn)子和紙料的接觸更迅速、頻率更高。同時,D形槽體非穩(wěn)定的D形湍流,一定程度上產(chǎn)生松散纖維,漿料湍流的發(fā)生不完全依賴于轉(zhuǎn)子旋轉(zhuǎn)。但與O形槽體相比,D形槽體存在容積小的缺點(diǎn),更為不足的是,D形槽體中有一部分漿料在撞擊槽壁時造成能量浪費(fèi),漿料撞擊垂直壁面較弧形壁面相比動能損失在22.7%左右。同時,從機(jī)械角度來看,D形流使槽體產(chǎn)生的較大振動以及轉(zhuǎn)子受力的不對稱性會影響機(jī)械部件的壽命。

    [1]2014年中國造紙行業(yè)年度研究報告[R/OL]. http://www.chinappi.org/reps/20150512154512227033.html.

    [2]Doshi M R. Paper recycling technology developments[J]. Progress in Paper Recycling, 2003, 12(2): 32.

    [3]陳遠(yuǎn)近. 離心機(jī)械設(shè)備在制漿工序的有效利用[J]. 中國造紙, 2008, 27(2): 69.

    [4]陳克復(fù), 張輝. 制漿造紙機(jī)械與設(shè)備[M]. 北京: 中國輕工業(yè)出版社, 2011.

    [5]陳賢根, 陳啟聲. 從專利文獻(xiàn)看水力碎漿機(jī)轉(zhuǎn)子節(jié)能工藝[J]. 能源與節(jié)能, 2012, 77(2): 38.

    [6]張林近, 葉旭初. 攪拌器內(nèi)湍流場的CFD模擬研究[J]. 南京工業(yè)大學(xué)學(xué)報, 2005, 27(2): 59.

    [7]馬青山, 聶毅強(qiáng), 包青云, 等. 攪拌槽內(nèi)三維流場的數(shù)值模擬[J]. 化工學(xué)報, 2003, 54(5): 612.

    [8]萬富濱, 袁力, 胡知輝. 基于CFD數(shù)值模擬的立式中心攪拌器有限元分析[J]. 石油礦場機(jī)械, 2008, 37(8): 52.

    [9]Bridgeman J, Jefferson B, Parsons S A. Computational Fluid Dynamics Modelling of Flocculation in Water Treatment: A Review[J]. Engineering Applications of Computational Fluid Mechanics, 2009, 3(2): 220.

    [10]Vlaev S D, Nikov I, Martinov M. Shear and Skin Friction on Particles in Power-law Fluids Agitated by Flat-blade and Fluid Foil Impellers[J]. Chemical Engineering Science, 2006, 61(2): 5455.

    [11]Cudmore G C, Holloway A G L, Gerber A G. A Model of Impeller Whirl for Baffled Mixing Vessels[J]. Journal of Fluids and Structures, 2015, 54: 719.

    [12]賴艷萍, 賴維, 盧曉江, 等. 圓柱型水力旋流器篩分紙漿纖維的數(shù)值模擬和實(shí)驗(yàn)研究[J]. 中國造紙學(xué)報, 2009, 24(1): 208.

    [13]陳慶蔚. 當(dāng)代廢紙?zhí)幚砑夹g(shù)[M]. 北京: 中國輕工業(yè)出版社, 1999.

    [14]肖明禎. 日本處理廢紙用水力碎漿機(jī)的幾項(xiàng)技術(shù)改進(jìn)[J]. 紙和造紙, 2000(5): 32.

    [15]G E Brown JR. Apparatus For Defibering Cellulosic Material In Water: US, 3420454[P].1969- 01- 07.

    [16]裴繼誠. 影響水力碎漿機(jī)效應(yīng)的幾個因素[J]. 天津造紙, 2004(1): 34.

    [17]Edward H Cumpston. Pulper: US, 3506202[P]. 1970- 4-14.

    [18]林躍春.立式碎漿機(jī)改進(jìn)的新動態(tài)[J].紙和造紙,2002(1):32.

    [19]Naoyuki Iwashige. Pulper For Pulverizing and Suspending Paper Stock: US, 2006/0273208 A1[P]. 2006-12- 07.

    [20]Wolfgang Mueller. Pulper for Comminuting and Suspending Paper Stock and Its Use: US, 2009/0272828A1[P]. 2009-11-5.

    [21]張健, 陳永利, 王晨, 等. 正向除渣器內(nèi)部纖維漿料流場的數(shù)值模擬[J]. 南京林業(yè)大學(xué)學(xué)報: 自然科學(xué)版, 2013, 37(2): 141.

    [22]周凌九, 胡德義, 王正偉. 可實(shí)現(xiàn)性k-ε模型在水輪機(jī)流場計算中的應(yīng)用[J]. 水動力學(xué)研究與進(jìn)展, 2003, 18(1): 68.

    [23]肖宗亮, 楊健, 林思達(dá), 等. 流漿箱噴嘴及湍動發(fā)生器的內(nèi)流場數(shù)值模擬[J]. 中國造紙學(xué)報, 2005, 20(1): 185.

    [24]余章書. 水力碎漿機(jī)葉輪的技術(shù)改造[J]. 中國造紙, 2015, 34(4): 68.

    [25]OptiSlushTMBale Pulpers: Optimized slushing result for the entire stock line and paper quality[EB/OL]. http://www.valmet.com/zh/services/pulping-and-fiber/recycled-fiber/.

    [26]陳奇峰, 陳克復(fù), 楊仁黨, 等. 基于CFD的中濃紙漿懸浮液在變徑管內(nèi)流動的數(shù)值模擬[J]. 中國造紙學(xué)報, 2004, 19(2): 96.

    [27]詹懷宇. 制漿原理與設(shè)備[M]. 北京: 中國輕工業(yè)出版社, 2010.

    [28]Ishkintana L K, Bennington C P J. Gas Holdup in Pulp Fibre Suspensions: Gas Voidage Profiles in a Batch-operated Sparged Tower[J]. The Canadian Journal of Chemical Engineering, 2015, 93: 486.

    [29]Clara Gomez. Numerical and Experimental Investigation of Macro-scale Mixing Applied to Pulp Fibre Suspensions[D]. Vancouver: The University of British Columbia, 2009: 7.

    [30]Frank Saville, Mark Martinez, James Olson. Variables Affecting Specific Energy Requirements for Repulping[J]. The Canadian Journal of Chemical Engineering, 2015, 93: 486.

    (*E-mail: zhnjfu@163.com)

    (責(zé)任編輯:關(guān)穎)

    Numerical Simulation of Fiber Slurry Flow Field inside the Vertical Hydraulic Pulper and Research on New Type of Its Tank Structure(Ι)

    DOU JingZHANG FangSHA Jiu-longZHANG Hui*

    (JiangsuProvincialKeyLabofPulpandPaperScienceandTechnology,NanjingForestryUniversity,Nanjing,JiangsuProvince, 210037)

    According to the assessment principles for pulpers including high efficiency, low energy consumption and running balance, based on the analysis of the structures of patented and applied vertical hydraulic pulpers domestically and abroad, using Fluent software to simulate the flow field of flow slurry at 5% consistency inside the O-shaped and D-shaped tanks, the distribution of pressure, velocity and turbulence intensity of both flow fields were obtained and analyzed. The results showed that the flow fields inside the O-shaped and D-shaped were vortex-like including a facade circulation and a horizontal swirling;the flow field inside the O-shaped was symmetrical and balanced, but the high-turbulence intensity area was limited to the rotor spinning region, and the turbulence intensity of the upper part was 1159% lower than that of the bottom part, the tangential velocity at the edge of the top of the flow field was zero and it was short of impact power;the flow field inside the D-shaped was with obviously increased turbulence intensity, and at the center of the D-shaped tank, the radius of low-turbulence intensity area was only 13.0% of the O-shaped tank; however, energy consumption inside the D-shaped tank was 22.7% higher caused by radially hitting the wall, and the D-shaped tank was with unbalanced running because of its unsymmetrical flow field.

    O-shaped vertical hydraulic pulper; D-shaped vertical hydraulic pulper; Fluent; numerical simulation; distribution of fiber slurry flow

    2015- 08-10

    南京林業(yè)大學(xué)江蘇省制漿造紙科學(xué)與技術(shù)重點(diǎn)實(shí)驗(yàn)室開放基金(201409);江蘇高校優(yōu)勢學(xué)科建設(shè)工程資助項(xiàng)目(PAPD)。

    竇靖,女,1990年生;在讀碩士研究生;主要從事制漿造紙節(jié)能減排技術(shù)與裝備方面的研究。

    *通信聯(lián)系人:張輝,E-mail:zhnjfu@163.com。

    TS734+.1

    A

    1000- 6842(2016)02- 0034- 09

    猜你喜歡
    槽體形槽漿料
    小型化雙向波束礦用定位終端天線設(shè)計
    隧道式食品油炸機(jī)槽體熱變形研究
    一種用于TPE塑料條的冷卻槽
    退鉻槽有限元分析及結(jié)構(gòu)優(yōu)化
    V形槽位置對偏導(dǎo)射流式伺服閥前置級液流特性的影響
    液壓與氣動(2019年7期)2019-07-18 11:11:46
    退鉻槽有限元分析及結(jié)構(gòu)優(yōu)化
    玻璃漿料鍵合中的孔洞抑制和微復(fù)合調(diào)控
    鋰離子電池漿料的研究
    柱塞泵球面配流盤阻尼槽對流量脈動性的影響
    液壓與氣動(2015年2期)2015-04-16 08:51:21
    V形槽角度快速測量
    丝袜人妻中文字幕| 可以免费在线观看a视频的电影网站| 性色av乱码一区二区三区2| 国产精品国产av在线观看| 婷婷精品国产亚洲av在线| 亚洲专区国产一区二区| 中文字幕av电影在线播放| 久久草成人影院| 中国美女看黄片| 亚洲欧美一区二区三区黑人| av视频免费观看在线观看| 悠悠久久av| 手机成人av网站| 欧美乱码精品一区二区三区| 搡老乐熟女国产| 色综合欧美亚洲国产小说| 99久久国产精品久久久| 十八禁人妻一区二区| 丝袜美腿诱惑在线| 男人舔女人下体高潮全视频| 国产精品综合久久久久久久免费 | 日韩三级视频一区二区三区| 国产又爽黄色视频| 999久久久精品免费观看国产| 一本综合久久免费| 香蕉久久夜色| 亚洲国产欧美网| 最新美女视频免费是黄的| 久久精品国产99精品国产亚洲性色 | 亚洲人成电影免费在线| 久久久水蜜桃国产精品网| 美女高潮喷水抽搐中文字幕| 国产成+人综合+亚洲专区| 欧美久久黑人一区二区| 午夜久久久在线观看| 极品人妻少妇av视频| 午夜精品在线福利| 男人操女人黄网站| 国产av精品麻豆| 琪琪午夜伦伦电影理论片6080| 黄片小视频在线播放| 国产精品免费一区二区三区在线| 国产精品久久久人人做人人爽| 久热这里只有精品99| 美女大奶头视频| 亚洲av成人不卡在线观看播放网| 纯流量卡能插随身wifi吗| 日本黄色视频三级网站网址| 女同久久另类99精品国产91| 手机成人av网站| 热99国产精品久久久久久7| 欧美日韩瑟瑟在线播放| 欧美日韩av久久| 一边摸一边抽搐一进一出视频| 99国产精品一区二区三区| 搡老熟女国产l中国老女人| 久久 成人 亚洲| 亚洲精品av麻豆狂野| 午夜免费鲁丝| av天堂久久9| 两个人看的免费小视频| 老司机福利观看| 91麻豆av在线| av网站在线播放免费| 国产亚洲欧美精品永久| 国产亚洲欧美精品永久| 亚洲男人的天堂狠狠| www.精华液| 一级片'在线观看视频| 熟女少妇亚洲综合色aaa.| 国产精品自产拍在线观看55亚洲| 亚洲精品国产精品久久久不卡| 亚洲第一青青草原| 长腿黑丝高跟| 亚洲狠狠婷婷综合久久图片| 黑丝袜美女国产一区| 俄罗斯特黄特色一大片| 精品久久久久久久毛片微露脸| 日韩高清综合在线| 久久这里只有精品19| 国产精品一区二区在线不卡| www.精华液| 亚洲avbb在线观看| 日韩欧美在线二视频| 老司机午夜福利在线观看视频| 成人国语在线视频| 精品日产1卡2卡| 亚洲第一欧美日韩一区二区三区| x7x7x7水蜜桃| 日本vs欧美在线观看视频| 女警被强在线播放| 99精品久久久久人妻精品| 精品熟女少妇八av免费久了| 久99久视频精品免费| 一区在线观看完整版| 麻豆国产av国片精品| tocl精华| 国产精品秋霞免费鲁丝片| 亚洲成a人片在线一区二区| 午夜成年电影在线免费观看| 真人做人爱边吃奶动态| 国产成人av激情在线播放| 高潮久久久久久久久久久不卡| 黄色毛片三级朝国网站| 一区二区三区国产精品乱码| 成人永久免费在线观看视频| 9191精品国产免费久久| av电影中文网址| www日本在线高清视频| 一级a爱视频在线免费观看| 变态另类成人亚洲欧美熟女 | 国产精品久久久久成人av| 亚洲av五月六月丁香网| 久热爱精品视频在线9| 侵犯人妻中文字幕一二三四区| 少妇粗大呻吟视频| 久久久国产欧美日韩av| 日韩成人在线观看一区二区三区| 国产又色又爽无遮挡免费看| 一区二区日韩欧美中文字幕| 精品乱码久久久久久99久播| 国产一区二区激情短视频| 亚洲av日韩精品久久久久久密| 国产精品自产拍在线观看55亚洲| 国产亚洲精品久久久久5区| 久久久久久久精品吃奶| 久久精品亚洲av国产电影网| 日韩有码中文字幕| www国产在线视频色| 夫妻午夜视频| 国产99久久九九免费精品| 国产成人精品久久二区二区免费| 亚洲一区中文字幕在线| 日韩中文字幕欧美一区二区| 色哟哟哟哟哟哟| 在线视频色国产色| 日韩精品免费视频一区二区三区| 日韩欧美三级三区| 免费高清视频大片| 欧美日韩中文字幕国产精品一区二区三区 | 99精品久久久久人妻精品| 日本vs欧美在线观看视频| 丰满人妻熟妇乱又伦精品不卡| 淫妇啪啪啪对白视频| 9191精品国产免费久久| 水蜜桃什么品种好| 国产无遮挡羞羞视频在线观看| av中文乱码字幕在线| 一区二区日韩欧美中文字幕| 午夜福利,免费看| 亚洲九九香蕉| 十八禁人妻一区二区| 久久香蕉激情| 亚洲成国产人片在线观看| 级片在线观看| 看黄色毛片网站| 久久久久国产精品人妻aⅴ院| 在线观看66精品国产| 纯流量卡能插随身wifi吗| 亚洲欧美日韩无卡精品| 午夜福利一区二区在线看| 男女床上黄色一级片免费看| av在线天堂中文字幕 | 国产麻豆69| 夫妻午夜视频| 窝窝影院91人妻| 欧美乱色亚洲激情| 岛国在线观看网站| 丰满饥渴人妻一区二区三| 久久久久亚洲av毛片大全| 日韩欧美一区视频在线观看| 十八禁网站免费在线| 久久人妻av系列| 国产成人精品久久二区二区免费| 美女国产高潮福利片在线看| 亚洲熟女毛片儿| 美女大奶头视频| 老汉色av国产亚洲站长工具| 51午夜福利影视在线观看| 真人一进一出gif抽搐免费| 亚洲一区高清亚洲精品| 精品第一国产精品| 亚洲色图综合在线观看| 日韩国内少妇激情av| 国产成人欧美| 色在线成人网| 亚洲第一青青草原| 欧美日韩乱码在线| 中国美女看黄片| 少妇的丰满在线观看| 丰满人妻熟妇乱又伦精品不卡| 一级毛片高清免费大全| 国产精品1区2区在线观看.| 亚洲精品在线观看二区| 精品人妻在线不人妻| 久久精品人人爽人人爽视色| 51午夜福利影视在线观看| 精品熟女少妇八av免费久了| 少妇的丰满在线观看| 国产精品电影一区二区三区| 高清欧美精品videossex| 免费在线观看影片大全网站| 男人舔女人下体高潮全视频| 亚洲熟妇熟女久久| 亚洲一区二区三区欧美精品| 99在线视频只有这里精品首页| 女警被强在线播放| 夫妻午夜视频| 色综合婷婷激情| 又紧又爽又黄一区二区| 色在线成人网| 亚洲成av片中文字幕在线观看| 久久精品国产亚洲av高清一级| 99在线视频只有这里精品首页| 在线av久久热| 国产精品一区二区免费欧美| 欧美人与性动交α欧美软件| 三上悠亚av全集在线观看| 久久天堂一区二区三区四区| 午夜亚洲福利在线播放| 日日干狠狠操夜夜爽| 亚洲一区二区三区欧美精品| 精品国产一区二区久久| 高清在线国产一区| 看片在线看免费视频| 国产精品1区2区在线观看.| 999精品在线视频| 欧美在线黄色| 国产欧美日韩精品亚洲av| 国产精品av久久久久免费| 嫁个100分男人电影在线观看| 两人在一起打扑克的视频| 91大片在线观看| 色婷婷久久久亚洲欧美| 日韩免费av在线播放| 男人舔女人下体高潮全视频| 国产精品一区二区三区四区久久 | 亚洲精品美女久久av网站| 亚洲九九香蕉| 午夜精品在线福利| 精品久久久久久电影网| 久久伊人香网站| 久久国产亚洲av麻豆专区| 欧美性长视频在线观看| 少妇粗大呻吟视频| 日本欧美视频一区| 99国产精品免费福利视频| 夜夜爽天天搞| 日韩有码中文字幕| www.精华液| 可以在线观看毛片的网站| 亚洲精品一卡2卡三卡4卡5卡| 日本五十路高清| 天堂动漫精品| 欧美在线一区亚洲| 免费在线观看黄色视频的| 欧美一级毛片孕妇| 久久国产精品人妻蜜桃| 日韩精品中文字幕看吧| 午夜视频精品福利| av在线播放免费不卡| 久久中文字幕人妻熟女| 欧美成狂野欧美在线观看| 首页视频小说图片口味搜索| 精品一区二区三区av网在线观看| 手机成人av网站| 亚洲av成人不卡在线观看播放网| 另类亚洲欧美激情| 精品免费久久久久久久清纯| 一区在线观看完整版| 欧美久久黑人一区二区| 精品福利观看| 波多野结衣av一区二区av| 国产一区二区激情短视频| 18美女黄网站色大片免费观看| 久久久久国产一级毛片高清牌| 麻豆一二三区av精品| 夜夜看夜夜爽夜夜摸 | 老汉色av国产亚洲站长工具| 日韩高清综合在线| 欧美丝袜亚洲另类 | av网站在线播放免费| 看免费av毛片| 精品一区二区三区视频在线观看免费 | 亚洲自偷自拍图片 自拍| 亚洲伊人色综图| 久久中文字幕一级| 视频区图区小说| 夜夜爽天天搞| 精品日产1卡2卡| 女人爽到高潮嗷嗷叫在线视频| 日韩免费av在线播放| 亚洲精品国产精品久久久不卡| 女警被强在线播放| 国产成人免费无遮挡视频| 欧美亚洲日本最大视频资源| 丰满的人妻完整版| 狠狠狠狠99中文字幕| 精品午夜福利视频在线观看一区| 高清欧美精品videossex| 日韩成人在线观看一区二区三区| 欧美激情高清一区二区三区| e午夜精品久久久久久久| 色播在线永久视频| 成人av一区二区三区在线看| 黄片小视频在线播放| 亚洲人成电影观看| 久久午夜综合久久蜜桃| 亚洲欧美精品综合一区二区三区| 国产在线观看jvid| 国产亚洲欧美在线一区二区| 精品国产一区二区久久| 精品高清国产在线一区| 精品一区二区三卡| 天堂中文最新版在线下载| 日本wwww免费看| 国产熟女xx| 嫩草影视91久久| 欧美日本中文国产一区发布| 国产精品成人在线| 久久青草综合色| 夜夜看夜夜爽夜夜摸 | 国产成人精品无人区| 精品久久蜜臀av无| 国产精品乱码一区二三区的特点 | 成人av一区二区三区在线看| 黄色视频不卡| 热99re8久久精品国产| xxx96com| 久久婷婷成人综合色麻豆| 亚洲精品美女久久久久99蜜臀| 在线观看日韩欧美| 又黄又粗又硬又大视频| 午夜福利,免费看| 亚洲人成电影免费在线| 黄色a级毛片大全视频| 国产一区二区三区视频了| 长腿黑丝高跟| 老司机午夜十八禁免费视频| 香蕉久久夜色| 美女大奶头视频| www.www免费av| 精品久久久久久久久久免费视频 | 国产精品影院久久| 美女高潮到喷水免费观看| 18禁黄网站禁片午夜丰满| 成人永久免费在线观看视频| 国产xxxxx性猛交| 午夜福利在线观看吧| 我的亚洲天堂| 不卡一级毛片| 日韩三级视频一区二区三区| 午夜老司机福利片| 9191精品国产免费久久| 成人三级黄色视频| 久久久久久久久久久久大奶| 高清在线国产一区| 精品国产超薄肉色丝袜足j| 亚洲精品美女久久av网站| 无人区码免费观看不卡| 精品国产国语对白av| 国产深夜福利视频在线观看| 满18在线观看网站| 国产成人啪精品午夜网站| 在线观看免费视频日本深夜| 99精品欧美一区二区三区四区| 国产欧美日韩一区二区精品| 免费在线观看日本一区| 午夜免费成人在线视频| 免费女性裸体啪啪无遮挡网站| 黑人猛操日本美女一级片| 亚洲 欧美 日韩 在线 免费| 男女床上黄色一级片免费看| 免费在线观看亚洲国产| 精品福利观看| 国产精品一区二区三区四区久久 | 看黄色毛片网站| 波多野结衣av一区二区av| 母亲3免费完整高清在线观看| 9色porny在线观看| 欧美一区二区精品小视频在线| 国产成人精品在线电影| 两人在一起打扑克的视频| 欧美不卡视频在线免费观看 | www.999成人在线观看| 国产成人一区二区三区免费视频网站| 好看av亚洲va欧美ⅴa在| 欧美av亚洲av综合av国产av| 69av精品久久久久久| 亚洲精品国产色婷婷电影| 免费在线观看影片大全网站| 午夜福利在线免费观看网站| 亚洲男人的天堂狠狠| 久久香蕉精品热| 成人免费观看视频高清| 少妇的丰满在线观看| 亚洲欧美日韩另类电影网站| 十八禁网站免费在线| 性少妇av在线| 亚洲精品一卡2卡三卡4卡5卡| 岛国视频午夜一区免费看| bbb黄色大片| 亚洲国产精品合色在线| 亚洲成av片中文字幕在线观看| 最近最新中文字幕大全免费视频| 老司机午夜福利在线观看视频| 99久久99久久久精品蜜桃| 女人爽到高潮嗷嗷叫在线视频| 黑人猛操日本美女一级片| 欧美日韩亚洲国产一区二区在线观看| 欧洲精品卡2卡3卡4卡5卡区| av网站免费在线观看视频| 日韩 欧美 亚洲 中文字幕| 深夜精品福利| 国产熟女午夜一区二区三区| 无限看片的www在线观看| 男女床上黄色一级片免费看| 国产成人欧美在线观看| 久久精品91蜜桃| 免费观看人在逋| 人人妻人人添人人爽欧美一区卜| 嫩草影视91久久| 午夜免费激情av| 亚洲国产欧美网| 亚洲国产毛片av蜜桃av| 国产成人精品无人区| 国产片内射在线| 美女大奶头视频| 久久性视频一级片| 欧美日韩瑟瑟在线播放| 成年人免费黄色播放视频| 日本三级黄在线观看| 日韩精品青青久久久久久| 国产精品久久视频播放| 亚洲av成人一区二区三| 国产欧美日韩一区二区三| 99精品久久久久人妻精品| 欧美日韩精品网址| 99在线人妻在线中文字幕| av片东京热男人的天堂| 午夜两性在线视频| 午夜福利免费观看在线| 国产精品久久电影中文字幕| www.自偷自拍.com| 黄频高清免费视频| 在线国产一区二区在线| 久久精品国产亚洲av香蕉五月| 日本vs欧美在线观看视频| 97人妻天天添夜夜摸| 国产免费现黄频在线看| av天堂久久9| 亚洲少妇的诱惑av| 久久九九热精品免费| 99国产精品一区二区三区| 成年人免费黄色播放视频| 国产亚洲欧美精品永久| 国产精品影院久久| 国产精华一区二区三区| 老司机深夜福利视频在线观看| 啦啦啦 在线观看视频| 亚洲avbb在线观看| 神马国产精品三级电影在线观看 | 亚洲七黄色美女视频| 在线观看免费高清a一片| 九色亚洲精品在线播放| 日韩成人在线观看一区二区三区| 中文字幕人妻熟女乱码| 亚洲一区中文字幕在线| 免费高清视频大片| 午夜老司机福利片| 热re99久久精品国产66热6| 老鸭窝网址在线观看| 久久国产亚洲av麻豆专区| 宅男免费午夜| 成人国语在线视频| 欧美一区二区精品小视频在线| 久久亚洲真实| 日韩免费高清中文字幕av| 成年人免费黄色播放视频| 成人影院久久| 国产在线观看jvid| 久9热在线精品视频| 国产亚洲av高清不卡| 久久人妻福利社区极品人妻图片| 怎么达到女性高潮| 18禁国产床啪视频网站| 精品日产1卡2卡| 好男人电影高清在线观看| 桃红色精品国产亚洲av| 亚洲午夜精品一区,二区,三区| 女性生殖器流出的白浆| 黄色丝袜av网址大全| 精品一区二区三卡| 女人精品久久久久毛片| 大型黄色视频在线免费观看| 久久九九热精品免费| 女同久久另类99精品国产91| 久久人妻熟女aⅴ| 一级毛片高清免费大全| 男女床上黄色一级片免费看| 午夜久久久在线观看| 国产精品久久久av美女十八| 久久这里只有精品19| 亚洲精品中文字幕在线视频| 麻豆av在线久日| 亚洲午夜精品一区,二区,三区| 久久人妻福利社区极品人妻图片| 99国产精品99久久久久| 丰满的人妻完整版| 可以免费在线观看a视频的电影网站| 亚洲精品中文字幕在线视频| 亚洲国产欧美一区二区综合| 又黄又爽又免费观看的视频| av在线天堂中文字幕 | 中文亚洲av片在线观看爽| 亚洲国产精品一区二区三区在线| 91老司机精品| 日韩精品免费视频一区二区三区| 女人精品久久久久毛片| 精品电影一区二区在线| 国产欧美日韩一区二区精品| 侵犯人妻中文字幕一二三四区| а√天堂www在线а√下载| 亚洲国产欧美网| 日韩三级视频一区二区三区| 精品一区二区三区av网在线观看| 国产精品久久久人人做人人爽| 欧美另类亚洲清纯唯美| 国产又爽黄色视频| 亚洲成国产人片在线观看| 新久久久久国产一级毛片| 精品国产一区二区久久| 欧美一区二区精品小视频在线| 国产熟女午夜一区二区三区| 日本黄色视频三级网站网址| 久久久久国内视频| 黄色成人免费大全| 一级片免费观看大全| 欧美激情久久久久久爽电影 | 日韩av在线大香蕉| 不卡av一区二区三区| 免费女性裸体啪啪无遮挡网站| 亚洲专区字幕在线| 午夜亚洲福利在线播放| 一本综合久久免费| 国产精品av久久久久免费| 日日夜夜操网爽| 亚洲成人国产一区在线观看| 91大片在线观看| 欧美老熟妇乱子伦牲交| 国产aⅴ精品一区二区三区波| 欧美丝袜亚洲另类 | 亚洲自拍偷在线| 国产国语露脸激情在线看| 老熟妇仑乱视频hdxx| 久9热在线精品视频| 一区二区三区激情视频| 十八禁人妻一区二区| 在线观看66精品国产| 亚洲国产欧美一区二区综合| 美女国产高潮福利片在线看| 日日干狠狠操夜夜爽| 超碰成人久久| 亚洲国产精品一区二区三区在线| 欧美人与性动交α欧美软件| 亚洲精品国产精品久久久不卡| 两个人看的免费小视频| 少妇被粗大的猛进出69影院| 日韩一卡2卡3卡4卡2021年| 亚洲精品在线观看二区| 性少妇av在线| 久久久国产欧美日韩av| 少妇的丰满在线观看| 久久国产精品男人的天堂亚洲| 国产亚洲欧美精品永久| 午夜a级毛片| 在线播放国产精品三级| 免费av毛片视频| av天堂久久9| 精品乱码久久久久久99久播| 久久人人爽av亚洲精品天堂| 91精品三级在线观看| 久久人人爽av亚洲精品天堂| 亚洲精品久久成人aⅴ小说| 国产精品 欧美亚洲| 国产欧美日韩一区二区三| 天堂影院成人在线观看| 美女高潮喷水抽搐中文字幕| 久久青草综合色| 国产免费现黄频在线看| 午夜91福利影院| 黄色视频,在线免费观看| 美女福利国产在线| 免费高清视频大片| 国产精品永久免费网站| 欧美黄色片欧美黄色片| 丁香六月欧美| 精品国产乱码久久久久久男人| 欧美日韩瑟瑟在线播放| a级毛片在线看网站| 激情视频va一区二区三区| 免费观看人在逋| 亚洲va日本ⅴa欧美va伊人久久| 亚洲黑人精品在线| 一二三四在线观看免费中文在| 久久久久九九精品影院| 悠悠久久av| e午夜精品久久久久久久| 久久人人精品亚洲av| 欧美日韩亚洲国产一区二区在线观看| 亚洲片人在线观看| 久久影院123| 亚洲成av片中文字幕在线观看| 日日夜夜操网爽| av在线天堂中文字幕 |