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

    擋板長度對輻流式沉淀池流場與污泥濃度場影響的數(shù)值模擬

    2016-08-04 06:47:25魏文禮李盼盼洪云飛劉玉玲
    關(guān)鍵詞:沉淀池數(shù)值模擬

    魏文禮,李盼盼,洪云飛,劉玉玲

    (西安理工大學(xué) 陜西省西北旱區(qū)生態(tài)水利工程重點(diǎn)實(shí)驗(yàn)室,陜西 西安 710048)

    ?

    擋板長度對輻流式沉淀池流場與污泥濃度場影響的數(shù)值模擬

    魏文禮,李盼盼,洪云飛,劉玉玲

    (西安理工大學(xué) 陜西省西北旱區(qū)生態(tài)水利工程重點(diǎn)實(shí)驗(yàn)室,陜西 西安 710048)

    [摘要]【目的】 研究進(jìn)水口處垂直擋板長度對沉淀池內(nèi)速度場和污泥質(zhì)量濃度場分布的影響,為實(shí)際工程中沉淀池的優(yōu)化設(shè)計提供參考。【方法】 采用簡化的多相流混合(Mixture)模型,對輻流式沉淀池內(nèi)液固兩相流水力特性進(jìn)行二維數(shù)值模擬,選取RNG k-ε紊流模型封閉時均流方程,采用有限體積法對微分方程進(jìn)行離散,使用壓力隱式算子分裂PISO(Pressure-Implicit with Splitting of Operators)算法進(jìn)行速度與壓力耦合求解。然后設(shè)初始時刻沉淀池內(nèi)為清水,進(jìn)水為含一定質(zhì)量濃度污泥的污水,在此條件下模擬并分析不同擋板長度時2個模型沉淀池內(nèi)速度場和污泥質(zhì)量濃度場的變化規(guī)律。【結(jié)果】 沉淀池內(nèi)有1個大回流區(qū)和2個小回流區(qū),長擋板模型的回流區(qū)比短檔板模型小;沉淀池內(nèi)污泥分布呈現(xiàn)分層現(xiàn)象,且長擋板模型內(nèi)污泥到達(dá)出口的時間比短檔板模型長;長擋板模型的小流速區(qū)域比短檔板模型大,更有利于污泥的沉降?!窘Y(jié)論】 輻流式沉淀池進(jìn)口處設(shè)置長擋板有利于提高沉淀池的運(yùn)行效率。

    [關(guān)鍵詞]沉淀池;擋板長度;數(shù)值模擬;速度場;污泥質(zhì)量濃度場

    沉淀池是污水處理過程中的一種重要構(gòu)筑物,輻流式沉淀池由于其排泥性能好、出水水質(zhì)高、占地面積少而在現(xiàn)代污水處理廠里被廣泛采用。沉淀池的容量和處理能力受諸多因素的影響,如固體濃度、池子形狀、沉淀機(jī)理、進(jìn)出口設(shè)計、擋板的位置和承載率等。

    近年來,國內(nèi)外研究者對沉淀池做了大量的研究,其中一部分采用計算流體力學(xué)(CFD)技術(shù)來模擬沉淀池的運(yùn)行情況。Krebs等[1-2]和Lyn等[3]人利用標(biāo)準(zhǔn)k-ε紊流模型和不同的污泥顆粒運(yùn)動方程[4-5]對平流式二沉池進(jìn)行了數(shù)值模擬,得到了池內(nèi)流速分布和污泥質(zhì)量濃度的分布情況。Tamayol等[6]研究發(fā)現(xiàn),擋板的最好位置是在回流區(qū)域,可用來破壞回流區(qū),從而可以提高初沉池的運(yùn)行效率。Razmi等[7]用試驗(yàn)和數(shù)值模擬2種方法研究了初沉池中擋板位置的影響。Liu等[8]用2D激光多普測速儀對一個矩形初沉池中的流場進(jìn)行測量,然后用研究結(jié)果優(yōu)化沉淀池的設(shè)計參數(shù),表明擋板附近的流速梯度很大,且流速梯度隨著水流變均勻而漸漸減小。秦博[9]通過在輻流式沉淀池內(nèi)設(shè)置斜管,對輻流式斜管沉淀池處理生活污水的效率進(jìn)行試驗(yàn)研究,基本掌握了斜管沉淀池的構(gòu)造功能及沉淀原理,并通過改變運(yùn)行因素驗(yàn)證了斜管沉淀池的處理效率。劉玉玲等[10-11]采用兩相流混合模型和RNGk-ε紊流模型,對輻流式二次沉淀池液固兩相流的力學(xué)特性進(jìn)行了數(shù)值模擬,得到了速度場、紊動能和污泥質(zhì)量濃度等參量的空間分布規(guī)律。蔡金傍等[12]采用標(biāo)準(zhǔn)k-ε紊流模型,對平流式沉淀池進(jìn)行了數(shù)值模擬,分析了不同擋板位置、進(jìn)水流速以及沉淀池長高比對沉淀池內(nèi)流速場和濃度場的影響,并計算分析了平流式沉淀池的過流曲線。曾光明等[13]利用渦量-流函數(shù)法建立了二維沉淀池速度場模型的控制方程,并用有限差分法進(jìn)行求解,利用二維濃度遷移方程對沉淀池濃度場進(jìn)行計算。蔣成義等[14]采用多相流歐拉模型、標(biāo)準(zhǔn)k-ε湍流模型和Boussinesq假設(shè)研究了二沉池中活性污泥和表面散熱引起的異重流。屈強(qiáng)等[15]利用改進(jìn)的RNGk-ε湍流模型和簡化的多相流Mixture模型,對輻流式二沉池內(nèi)的速度場和污泥質(zhì)量濃度場進(jìn)行了數(shù)值模擬。金光等[16]利用計算流體動力學(xué)分析軟件建立可以實(shí)現(xiàn)混凝沉淀的新型輻流式沉淀池數(shù)值模型,并對其可行性進(jìn)行了數(shù)值模擬研究,表明在傳統(tǒng)輻流式沉淀池內(nèi)增設(shè)強(qiáng)化擾流和改變水流方向的反應(yīng)罩和導(dǎo)流板,可以使輻流式沉淀池內(nèi)部流場紊流強(qiáng)度明顯增強(qiáng)。綜上,目前有關(guān)沉淀池的研究較多,但關(guān)于輻流式沉淀池進(jìn)口擋板長度對流場和污泥濃度場影響的研究很少。

    為此,本研究采用數(shù)值模擬的方法,分別建立進(jìn)水口處設(shè)一短擋板和一長擋板的二維模型,采用FLUENT軟件進(jìn)行固液兩相流數(shù)值模擬。然后通過分析輻流式沉淀池內(nèi)固液兩相流的速度場和污泥濃度場,對進(jìn)水口處不同垂直擋板長度對沉淀池運(yùn)行效率的影響進(jìn)行研究,以期為實(shí)際工程中沉淀池的優(yōu)化設(shè)計提供參考。

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

    多相流模型可用來模擬沉淀池中的流場。Mixture模型[17]是一種簡化的多相流模型,其利用小空間尺度上的局部平衡來求解混合相的動量、連續(xù)性,以及第二相的體積分率、滑移速度和漂移速度??紤]該模型具有模型簡單、計算量小且結(jié)果較為可靠等特點(diǎn),這里采用Mixture模型。

    Mixture模型連續(xù)性方程為:

    (1)

    (2)

    (3)

    Mixture模型動量方程可以通過對所有相各自的動量方程求和來獲得。其形式如下:

    (4)

    (5)

    vdr,k=vk-vm。

    (6)

    式中:p為壓強(qiáng),μm為混合黏度,g為重力加速度,F(xiàn)為體積力,vdr,k為第二相k的飄移速度。

    滑流速度被定義為第二相(p)的速度相對于主相(q)的速度,有:

    vqp=vp-vq。

    (7)

    式中:vqp為滑移速度,vp為p相的速度,vq為q相的速度。

    滑流速度和飄移速度的關(guān)系為:

    (8)

    式中:vdr,p為p相的飄移速度,vqk表示第k相的滑移速度。

    從第二相p的連續(xù)方程,可以得到第二相p相的體積分?jǐn)?shù)方程為:

    (9)

    式中:αp表示p相的體積分?jǐn)?shù),ρp表示p相的密度。

    利用改進(jìn)的RNGk-ε紊流模型[18],可以通過修正湍動黏度,在考慮平均流動中的旋轉(zhuǎn)及旋流流動情況的基礎(chǔ)上,提高計算精度并可以正確處理近壁區(qū)域,可以更好地處理高應(yīng)變率及流線彎曲程度較大的流動。其紊動能k和耗散率ε的運(yùn)輸方程分別為:

    (10)

    (11)

    以上方程構(gòu)成了求解流場分布規(guī)律的封閉方程組,根據(jù)實(shí)際工況施加相應(yīng)的邊界條件后,就構(gòu)成該方程組的定解問題。

    2建模和邊界條件

    2.1沉淀池模型

    以中心進(jìn)水的輻流式沉淀池為研究對象,其處理能力大約是150 000 m3/d,池子容積大約是2 960 m3。該沉淀池中心設(shè)環(huán)形進(jìn)水口進(jìn)水,四周是圓形圍堰,出口為設(shè)置在圍堰內(nèi)側(cè)附近的一個環(huán)形出水槽,沉淀池中經(jīng)沉淀過后的清水溢流到出水槽中排出。沉淀池底部污泥斗的坡度為12.5°,沉淀池進(jìn)水口圓環(huán)半徑r=3.75 m。

    輻流式沉淀池內(nèi)的流動屬于三維運(yùn)動,考慮到輻流式沉淀池的幾何對稱性,為了降低網(wǎng)格數(shù),減少計算時間,忽略環(huán)向流動,特簡化為僅有軸向和徑向的二維運(yùn)動,在此基礎(chǔ)上建立2個二維模型:一個是進(jìn)水口后端設(shè)有一小段垂直擋板,高1.3 m;另一個是進(jìn)水口后端設(shè)有一長擋板,擋板走勢跟外墻相同,高約3 m,兩模型其他尺寸都一樣,其幾何尺寸如圖1所示。為敘述方便,下文將2種模型分別用a模型和b模型表示:a模型代表進(jìn)口處垂直擋板為短擋板的模型,b模型代表進(jìn)口處垂直擋板為長擋板的模型。

    圖 1不同擋板長度的二維模型的計算區(qū)域和尺寸

    Fig.1Computational region and size for two-dimensional model of different baffle lengths

    2.2網(wǎng)格劃分

    采用GAMBIT軟件建模,使用結(jié)構(gòu)化網(wǎng)格對整體區(qū)域進(jìn)行劃分。由于進(jìn)水口和出水口都在沉淀池上部,為了提高計算精度,故對這部分網(wǎng)格在垂直方向進(jìn)行適度加密處理。2種模型網(wǎng)格劃分方法一致,劃分結(jié)果如圖2所示,網(wǎng)格總單元數(shù)為7 444個。

    圖 2 二維計算區(qū)域的網(wǎng)格圖Fig.2 Computational mesh for two-dimensional region

    2.3邊界條件及求解方法

    本次試驗(yàn)?zāi)M了沉淀池中的固液兩相流,主相為水,密度為ρ=1 000.35 kg/m3,動力黏度μ動=0.001 005 Pa·s;次相為污泥,密度為ρ=1 051 kg/m3,時均黏度μ時=0.020 01 Pa·s,并且假設(shè)固體為均勻顆粒,顆粒直徑為d=0.001 m,進(jìn)水中污泥的質(zhì)量濃度C=35.734 kg/m3。

    進(jìn)口邊界條件給定速度值,固液兩相的流速相同,v=0.085 m/s;出水口邊界條件采用壓力出流,相對壓強(qiáng)為0。池底和邊壁(包括擋板)為固體壁面,

    使用標(biāo)準(zhǔn)壁面函數(shù)法確定,頂面采用“剛蓋”假定。

    假設(shè)初始時刻沉淀池中充滿清水,污泥質(zhì)量濃度為0,由于沉淀池中污泥質(zhì)量濃度隨計算時間的增加而變化,故全場采用非恒定流進(jìn)行計算。采用簡化的Mixture多相流模型模擬混合流場,用RNGk-ε模型封閉時均流方程,控制方程的離散采用有限體積法,速度與壓力耦合求解時使用壓力隱式算子分裂PISO算法,紊動能、耗散率、動能方程均采用二階迎風(fēng)(Second Order Upwind)離散格式。

    本次試驗(yàn)是為了研究污泥在沉淀池中的遷移規(guī)律和質(zhì)量濃度場隨時間的變化,不依計算過程中進(jìn)出口流量相等為計算平衡的依據(jù),取時間步長=0.005 s,計算總時間t=1 000 s,并且2種模型的各參數(shù)和邊界條件完全相同。

    3結(jié)果與分析

    3.1流線分布

    沉淀池中的流線隨時間的分布如圖3所示。

    圖3表明,隨著時間的推移,在模擬時間為300 s左右時,顯示在擋板附近有1個小回流區(qū),隨著模擬時間增加,回流區(qū)向后部推移;模擬時間為500 s時,在池子中部形成1個大回流區(qū);模擬時間等于700 s時,沉淀池底部開始產(chǎn)生1個小回流區(qū),到 1 000 s左右,池子底部形成2個小回流區(qū),貼近池子底板,再往后2個模型的流線圖均基本同1 000 s時一樣,池中一直保持3個回流區(qū)不變。之前的文獻(xiàn)[10,15]提到沉淀池中是2個回流區(qū),這可能是沉淀池體型不同所致,也可能是文獻(xiàn)[10,15]中模擬計算時間不夠的緣故?;亓鲄^(qū)基本上都是死區(qū),回流區(qū)的存在降低了沉淀池的有效利用容積,是影響沉淀池去除率的一個重要因素。

    流線的疏密代表了流速的大小,擋板附近和池子中部回流區(qū)流線較密,說明該處流速較大。水面附近和池子后區(qū)流線稀疏,說明該處流速較低,這有利于固體顆粒沉降到池底。

    比較a、b 2個模型,流線的變化規(guī)律基本一致,回流區(qū)的位置和個數(shù)也相近,說明擋板的長短并不能減少回流區(qū)的個數(shù)。但是每個時刻模型b中的回流區(qū)面積比a模型小很多,比如模擬時間為300和500 s時,a中的回流區(qū)面積接近b回流區(qū)面積的2倍;模擬時間等于700 s時,a中的回流區(qū)面積約占整個池子面積的3/4,而此時b的回流區(qū)只占池子面積的1/2左右。回流區(qū)面積減少,有效流動區(qū)域增大,就會相應(yīng)提高沉淀池的運(yùn)行效率。因此,進(jìn)口處長擋板與短擋板相比,其可以大大減小回流區(qū)面積,提高沉淀池運(yùn)行效率。但如何減少回流區(qū)的個數(shù),今后還要進(jìn)行更多的研究。

    3.2污泥質(zhì)量濃度的分布

    污泥在沉淀池中的遷移和質(zhì)量濃度分布如圖4所示。

    圖 4 2個模型不同時刻污泥質(zhì)量濃度的分布Fig.4 Sediment mass concentration distribution of the two models at different times

    圖4顯示,模擬開始時污泥從擋板下沿向池內(nèi)擴(kuò)散,到模擬時間為300 s時,池內(nèi)存在著清晰的泥水分界面,在此界面以上,污泥質(zhì)量濃度為0,為清水區(qū),此時整個沉淀池污泥質(zhì)量質(zhì)量濃度很小,不超過137 g/L;隨時間增長,污泥開始在池底沉淀下來,形成成層沉淀區(qū),在模擬時間為500 s時,成層沉淀區(qū)的污泥質(zhì)量濃度為274~343 g/L;到模擬時間為700 s時,可清晰看到池底的污泥壓縮層,壓縮層的污泥質(zhì)量濃度大于480 g/L,此時污泥在沉淀池內(nèi)的分布規(guī)律為:上面是清水區(qū),泥水分界面以下是污泥質(zhì)量濃度逐漸升高的絮凝區(qū),接下來是污泥質(zhì)量濃度基本一致的成層沉淀區(qū),池底有一層較薄的污泥壓縮區(qū)。隨著模擬時間增加,壓縮區(qū)厚度越來越大,該區(qū)域污泥的質(zhì)量濃度也逐漸增大。污水區(qū)中絮凝區(qū)的范圍最大,且不斷向清水區(qū)擴(kuò)展。成層沉淀區(qū)和壓縮區(qū)之間也存在明顯的壓縮界面,最后壓縮區(qū)的高度大于成層沉淀區(qū)的高度。污泥區(qū)逐漸擴(kuò)大,清水區(qū)逐漸縮小,到最后污泥擴(kuò)散到出口處,充滿整個沉淀池,此時出流中開始夾帶污泥,如果不及時對池中污泥進(jìn)行清理,沉淀池的處理能力就會逐漸喪失。

    從污泥質(zhì)量濃度分布圖(圖4)可以看出,2種模型的污泥壓縮區(qū)都集中在左下角的污泥斗中,從污泥斗往外,污泥質(zhì)量濃度逐漸減小。這非常符合實(shí)際情況,證明了數(shù)值模擬的正確性?,F(xiàn)實(shí)中的沉淀池在污泥斗處設(shè)有污泥出口,污泥堆積到一定程度時,用專門的設(shè)備將污泥斗中的污泥排出,才可使沉淀池能夠循環(huán)運(yùn)行。

    比較a、b 2個模型,其污泥的運(yùn)移規(guī)律相同,但質(zhì)量濃度分布存在差異。a模型中污泥擴(kuò)散明顯較b快很多:模擬時間為700 s時,a模型中污泥已經(jīng)擴(kuò)散到出口處,污泥開始從出口流出,沉淀池處理能力喪失,而b模型中污泥在模擬時間為1 000 s時污泥還未到達(dá)出口,出流仍然是清水,沉淀池依然運(yùn)行正常。由此可見,進(jìn)口處長擋板能明顯提高沉淀池的污泥處理能力,延長沉淀池運(yùn)行周期,提高去除率,降低能耗。

    3.3流速分布

    流速云圖能清晰地顯示出流速場的分布情況,模擬時間為300和500 s時的流速場分布如圖5所示。由圖5可以發(fā)現(xiàn),a、b兩模型流速場的差異也很明顯。不同的顏色深度代表了不同的流速大?。簭陌咨珴u變到黑色,流速逐漸增大,顏色深度變化越快的地方,表示流速梯度越大,沉淀池內(nèi)流速的分布范圍為0.016~0.250 m/s。從圖5可以看出,沉淀池內(nèi)流速場的分布規(guī)律為:大流速區(qū)域分布在沉淀池前部,進(jìn)口附近流速最大,在擋板的束縛下,高速水流流向沉淀池下部,然后從擋板下沿開始向池中擴(kuò)散。云圖中顏色深度由黑色逐漸向白色轉(zhuǎn)變,表示流速逐漸減小。沉淀池中后部分基本為白色,表明流場較均勻,該區(qū)域流速大小為0.016 m/s左右,小于0.03 m/s,這樣的流速有利于固體顆粒的沉降。

    圖 52個模型不同時刻的流速云圖

    Fig.5Velocity magnitude contour of the two models at different times

    從圖5還可以看出,b模型長擋板的束流效果比a模型短擋板好:模擬時間為300 s和500 s時,b模型的流速小于0.03 m/s區(qū)域明顯大于a模型,說明b模型小流速區(qū)域比a模型大,這更有利于污泥的沉降。另外從圖5左邊2幅圖可以看到,a模型中流速為0.03~0.063 m/s的區(qū)域處于沉淀池頂部,這樣的流場會將小顆粒污泥卷揚(yáng)到沉淀池上部,最后隨上部水流從出口流出,使沉淀效果降低。500 s時,b模型左下角污泥斗處顏色變成白色,表明此處流速降到0.016 m/s左右,這有利于污泥在污泥斗中的沉淀。其他時刻2個模型的流場圖對比效果與以上結(jié)果雷同,不再贅述。

    4結(jié)論

    本研究采用數(shù)值模擬方法,模擬了輻流式沉淀池內(nèi)流場的變化和污泥的運(yùn)移規(guī)律,通過對不同模擬時間的計算分析,可以清楚看到流場中回流區(qū)的產(chǎn)生、變化,最后在流場內(nèi)形成大小不等的3個回流區(qū);模擬演示了沉淀池內(nèi)污泥的遷移過程,污泥在沉淀池內(nèi)的分布表現(xiàn)出分層現(xiàn)象:從上往下,污泥濃度遞增,最底層出現(xiàn)污泥壓縮區(qū)。這些直觀的顯示在實(shí)驗(yàn)室內(nèi)很難達(dá)到,體現(xiàn)了數(shù)值模擬的優(yōu)勢。

    另外,本研究通過在進(jìn)水口處設(shè)置不同長度的擋板,模擬并分析進(jìn)水口處擋板不同長度對流場的影響。從模擬結(jié)果可以看出:進(jìn)水口處設(shè)置長擋板可以明顯減小流場內(nèi)的回流面積,增大有效流動區(qū)域,提高沉淀池的污泥處理能力,延長沉淀池運(yùn)行周期,提高去除率,改善出流水質(zhì),這一結(jié)論可以為沉淀池的設(shè)計和運(yùn)行提供參考。為了進(jìn)一步減小沉淀池內(nèi)的回流區(qū)個數(shù)和回流區(qū)大小,今后還需結(jié)合試驗(yàn)和數(shù)值模擬,進(jìn)行更多的相關(guān)研究。

    [參考文獻(xiàn)]

    [1]Krebs P,Armbruster M,Rodi W.Laboratory experi-ments of buoyancy-influenced flow in clarifiers [J].Journal of Hydraulic Research,1998,36(5):831-851.

    [2]Krebs P,Stamou A I,Garcia-Heras J L.Influence of inlet and outlet configuration on the flow in secondary clarifiers [J].Water Science and Technology,1996,34(5/6):1-9.

    [3]Lyn D A,Stamou A I,Rodi W.Density currents and shear-induced flocculation in sedimentation tanks [J].Journal of Hydraulic Engineering,1992,118(6):849-867.

    [4]Takacs I,Patry G G,Nolasco D.A dynamic model of the clarification-thickening process [J].Water Research,1991,25(10):1263-1271.

    [5]Zhou S,McCorquodale J A.Modeling of rectangular settling ta-nks [J].Journal of Hydraulic Engineering,1992,118(10):1391-1405.

    [6]Tamayol A,F(xiàn)iroozabadi B,Ahmadi G.Effects of inlet position and baffle configuration on hydraulic performance of primary settling tanks [J].Journal of Hydraulic Engineering,2008,134(7):1004-1009.

    [7]Razmi A,F(xiàn)iroozabadi B,Ahmadi G.Experimental and numerical approach to enlargement of performance of primary settling tanks [J].Journal of Applied Fluid Mechanics,2009,2(1):1-12.

    [8]Liu B,Ma J,Luo J,et al.Two-dimensional LDV measurement,modeling,and optimal design of rectangular primary settling tanks [J].Journal of Environmental Engineering,2010,136(5):501-507.

    [9]秦博.輻流式斜管沉淀池處理生活污水的實(shí)驗(yàn)研究 [J].環(huán)境科學(xué)與管理,2012,37(4):82-85.

    Qin B.Research on domestic wastewater treatment with radial flow inclined tube settling tank [J].Environmental Science and Management,2012,37(4):82-85.

    [10]劉玉玲,張沛,魏文禮.輻流式沉淀池液固兩相流力學(xué)特性二維數(shù)值模擬 [J].武漢大學(xué)學(xué)報(工學(xué)版),2013,46(4):410-412.

    Liu Y L,Zhang P,Wei W L.Two-dimensional numerical simulation of properties of liquid-solid two-phase flow in a circular secondary clarifier [J].Engineering Journal of Wuhan University,2013,46(4):410-412.

    [11]劉玉玲,張沛,魏文禮,等.輻流式沉淀池液固兩相流力學(xué)特性三維數(shù)值模擬 [J].水資源與水工程學(xué)報,2013,24(4):25-27.

    Liu Y L,Zhang P,Wei W L,et al.Numerical simulation of mechanical property of solid-liquid two-phase turbulent flow in a secondary sedimen-tation tank of radial flow [J].Journal of Water Resources & Water Engineering,2013,24(4):25-27.

    [12]蔡金傍,朱亮,段祥寶.平流式流沉池數(shù)值模擬分析 [J].河海大學(xué)學(xué)報(自然科學(xué)版),2004,32(1):27-31.

    Cai J B,Zhu L,Duan X B.Numerical simulation for rectangular settling tanks [J].Journal of Hohai Unversity(Natural Science),2004,32(1):27-31.

    [13]曾光明,葛衛(wèi)華,秦肖生,等.污水廠二維沉淀池水流和懸浮物運(yùn)動數(shù)值模擬 [J].中國環(huán)境科學(xué),2002,22(4):338-341.

    Zeng G M,Ge W H,Qin X S,et al.Numerical modeling on the movement of water and SS in two-dimensional sedimentation tanks of sewage factory [J].China Environmental Science,2002,22(4):338-341.

    [14]蔣成義,吳春篤,黃衛(wèi)東,等.輻流式二沉池中異重流的計算流體力學(xué)模型研究 [J].水科學(xué)進(jìn)展,2007,18(6):846-852.

    Jiang C Y,Wu C D,Huang W D,et al.Computational fluid dynamic model for density currents in activated sludge secondary clarifier [J].Advances in Water Science,2007,18(6):846-852.

    [15]屈強(qiáng),馬魯銘,王紅武.輻流式二沉池固液兩相流數(shù)值模擬 [J].同濟(jì)大學(xué)學(xué)報(自然科學(xué)版),2006,34(9):1212-1216.

    Qu Q,Ma L M,Wang H W.Solid-liquid two-phase flow simulation in a circular secondary clarifier [J].Journal of Tongji Unversity(Natural Science Edition),2006,34(9):1212-1216.

    [16]金光,張保生,張菊,等.實(shí)現(xiàn)輻流式沉淀池混凝沉淀的數(shù)值模擬研究 [J].人民黃河,2011,33(5):41-45.

    Jin G,Zhang B S,Zhang J,et al.Realization of radial flow sedimentation tank coagulation precipitation numerical simulation research [J].Yellow River,2011,33(5):41-45.

    [17]戴會超,槐文信,吳玉林,等.水利水電工程水流精細(xì)模擬理論與應(yīng)用 [M].北京:科學(xué)出版社,2008.

    Dai H C,Huai W X,Wu Y L,et al.Theories and applications of flow refined simulation for hydraulic and hydropower engineering [M].Beijing:Science and Technology Press,2008.

    [18]魏文禮,戴會超.紊流模型理論及工程應(yīng)用 [M].西安:陜西科學(xué)技術(shù)出版社,2006.

    Wei W L,Dai H C.Turbulence model theory and engineering applications [M].Xi’an:Shaanxi Science and Technology Press,2006.

    DOI:網(wǎng)絡(luò)出版時間:2016-06-0816:2110.13207/j.cnki.jnwafu.2016.07.032

    [收稿日期]2014-11-03

    [基金項目]國家自然科學(xué)基金項目(51578452);陜西省工業(yè)科技攻關(guān)項目(2016GY-180);陜西省科學(xué)研究計劃項目(2014K15-03-05)

    [作者簡介]魏文禮(1965-),男,陜西大荔人,教授,博士,主要從事環(huán)境水力學(xué)、水污染控制理論與技術(shù)研究。 E-mail:wei_wenli@126.com

    [中圖分類號]X703

    [文獻(xiàn)標(biāo)志碼]A

    [文章編號]1671-9387(2016)07-0228-07

    Influence of baffle length on flow and sludge concentration fields in a radial sedimentation tank

    WEI Wenli,LI Panpan,HONG Yunfei,LIU Yuling

    (StateKeyLaboratoryofEco-HydraulicEngineeringinShaanxi,Xi’anUniversityofTechnology,Xi’an,Shaanxi710048,China)

    Abstract:【Objective】 The influence of inlet vertical baffle on distribution of flow field and sludge mass concentration field was studied in a sedimentation tank to provide reference for optimization design of actual sedimentation tanks.【Method】 The simplified Mixture multiphase model was used to study the hydraulic characteristics of solid-liquid two-phase flow in a radial sedimentation tank by 2D numerical simulation.The RNG k-ε turbulent model was used to close the time-averaged flow equations.The finite volume method was used to discretize the differential equations,and the coupling velocity and pressure equations were solved with the pressure-implicit method of splitting operators (PISO) algorithm.At the initial time,the sedimentation tank was filled with stationary clear water and sewage ran through inlet.The variation of flow field and sludge mass concentration field in different times was then studied.【Result】 There were three recirculation zones in the sedimentation tank,a big one and two small ones.The model with a longer feed flow baffle caused a smaller recirculation zone.The distribution of sludge showed stratified phenomenon in the sedimentation tank.Sludge took a longer time to arrive at export in the model with a longer feed flow baffle.The region of fluid particles with a smaller velocity was bigger in the model with a longer feed flow baffle,which was beneficial to sludge sedimentation.【Conclusion】 The sedimentation tank with a longer feed flow baffle can produce better running efficiency.

    Key words:sedimentation tank;baffle length;numerical simulation;flow field;sludge mass concentration field

    網(wǎng)絡(luò)出版地址:http://www.cnki.net/kcms/detail/61.1390.S.20160608.1621.064.html

    猜你喜歡
    沉淀池數(shù)值模擬
    新型異形橢圓沉淀池的CFD數(shù)值模擬研究
    污水處理沉淀池進(jìn)水系統(tǒng)的改進(jìn)運(yùn)用
    長興水廠沉淀池反應(yīng)區(qū)浮沫現(xiàn)象的探討和解決方案
    水廠沉淀池改造前后排泥水含固率研究
    生物化工(2021年3期)2021-07-10 09:53:44
    污水處理廠沉淀池剖析——以烏魯木齊某污水處理廠為例
    張家灣煤礦巷道無支護(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ī)的流場數(shù)值模擬研究
    科技視界(2016年22期)2016-10-18 14:53:19
    国产在线精品亚洲第一网站| 亚洲精品国产区一区二| 午夜精品在线福利| 日本a在线网址| 欧美日韩乱码在线| 日韩欧美三级三区| 在线永久观看黄色视频| 亚洲欧美日韩无卡精品| 亚洲 国产 在线| 国产视频一区二区在线看| 欧美乱色亚洲激情| 日本在线视频免费播放| 国产亚洲av嫩草精品影院| а√天堂www在线а√下载| 欧美日韩中文字幕国产精品一区二区三区 | 少妇裸体淫交视频免费看高清 | 久久九九热精品免费| 黄色视频,在线免费观看| 日韩三级视频一区二区三区| 黄片小视频在线播放| 在线视频色国产色| 国产精品久久久久久亚洲av鲁大| 国产一卡二卡三卡精品| 变态另类成人亚洲欧美熟女 | 亚洲第一av免费看| 男人舔女人的私密视频| 精品国产亚洲在线| 国产精品1区2区在线观看.| 老司机深夜福利视频在线观看| 男男h啪啪无遮挡| 高清在线国产一区| 久久久国产成人精品二区| 久久这里只有精品19| 好男人在线观看高清免费视频 | 黑丝袜美女国产一区| 国产三级黄色录像| 国产精品,欧美在线| 日韩精品免费视频一区二区三区| 99久久久亚洲精品蜜臀av| 国产亚洲av高清不卡| 欧美成人午夜精品| 精品国产亚洲在线| 一区二区三区高清视频在线| 久久婷婷成人综合色麻豆| 黄网站色视频无遮挡免费观看| 亚洲va日本ⅴa欧美va伊人久久| 真人做人爱边吃奶动态| 国产在线精品亚洲第一网站| 最新美女视频免费是黄的| 成熟少妇高潮喷水视频| 亚洲一区高清亚洲精品| 三级毛片av免费| 在线观看免费午夜福利视频| 中亚洲国语对白在线视频| 丝袜美腿诱惑在线| 19禁男女啪啪无遮挡网站| 日本精品一区二区三区蜜桃| 黄片小视频在线播放| 国产精品免费一区二区三区在线| 国产免费av片在线观看野外av| 国产一区在线观看成人免费| 最好的美女福利视频网| 伊人久久大香线蕉亚洲五| 国产欧美日韩一区二区三区在线| 亚洲 欧美 日韩 在线 免费| 亚洲三区欧美一区| 久久精品91蜜桃| 亚洲 欧美 日韩 在线 免费| 一区二区三区激情视频| 国产视频一区二区在线看| 在线观看免费视频日本深夜| 午夜激情av网站| 又黄又粗又硬又大视频| 亚洲一区二区三区不卡视频| 女人爽到高潮嗷嗷叫在线视频| 亚洲国产欧美网| 国产成人精品久久二区二区免费| 99国产精品99久久久久| 欧美久久黑人一区二区| 级片在线观看| 欧美国产日韩亚洲一区| 日本在线视频免费播放| 成人亚洲精品av一区二区| 婷婷丁香在线五月| 亚洲avbb在线观看| 精品日产1卡2卡| 欧美日韩亚洲国产一区二区在线观看| 久久精品国产99精品国产亚洲性色 | 久久中文字幕人妻熟女| 久久久久久久午夜电影| 国产一区二区激情短视频| 精品久久久久久久毛片微露脸| 成人亚洲精品一区在线观看| 精品国产乱码久久久久久男人| 国产亚洲精品久久久久5区| 一边摸一边抽搐一进一小说| 成人手机av| 真人一进一出gif抽搐免费| av片东京热男人的天堂| 久久午夜综合久久蜜桃| 国产高清videossex| 午夜成年电影在线免费观看| 亚洲男人的天堂狠狠| 国产激情欧美一区二区| 91在线观看av| 女性被躁到高潮视频| 欧美日韩黄片免| 午夜视频精品福利| 亚洲av成人不卡在线观看播放网| 免费少妇av软件| 日本 av在线| 亚洲激情在线av| 啦啦啦 在线观看视频| 亚洲中文av在线| 亚洲人成77777在线视频| 国产av精品麻豆| 少妇粗大呻吟视频| 嫩草影视91久久| av欧美777| 久久香蕉国产精品| av在线天堂中文字幕| 久久人妻熟女aⅴ| 国产成人精品无人区| 啪啪无遮挡十八禁网站| 一本久久中文字幕| 日本撒尿小便嘘嘘汇集6| 熟女少妇亚洲综合色aaa.| 黑人欧美特级aaaaaa片| 免费一级毛片在线播放高清视频 | 免费高清视频大片| 91成年电影在线观看| 两个人视频免费观看高清| www.www免费av| 后天国语完整版免费观看| 99久久国产精品久久久| 久久草成人影院| 日韩精品免费视频一区二区三区| 91老司机精品| 国产在线精品亚洲第一网站| 精品国产乱子伦一区二区三区| 亚洲国产欧美网| 亚洲色图av天堂| 国产亚洲精品久久久久久毛片| 亚洲五月色婷婷综合| 十八禁人妻一区二区| 这个男人来自地球电影免费观看| 99在线视频只有这里精品首页| 亚洲欧美激情综合另类| 亚洲天堂国产精品一区在线| 欧美中文综合在线视频| 操出白浆在线播放| av福利片在线| 欧美丝袜亚洲另类 | 精品国产乱子伦一区二区三区| 女人爽到高潮嗷嗷叫在线视频| 亚洲精品一卡2卡三卡4卡5卡| 久久久久久免费高清国产稀缺| 一区二区三区高清视频在线| 亚洲国产精品sss在线观看| 两个人免费观看高清视频| 国产av在哪里看| 国产精品久久久人人做人人爽| 90打野战视频偷拍视频| 一区二区三区精品91| 99久久综合精品五月天人人| 免费一级毛片在线播放高清视频 | 满18在线观看网站| 一区在线观看完整版| 多毛熟女@视频| 黑人操中国人逼视频| 亚洲成人久久性| 国产精品久久久人人做人人爽| 久久精品国产亚洲av高清一级| 久久久精品国产亚洲av高清涩受| 最新美女视频免费是黄的| 亚洲色图av天堂| 无遮挡黄片免费观看| 露出奶头的视频| 亚洲成国产人片在线观看| 97人妻天天添夜夜摸| 91精品三级在线观看| 日韩一卡2卡3卡4卡2021年| 亚洲国产精品合色在线| АⅤ资源中文在线天堂| 成人亚洲精品av一区二区| 国产一级毛片七仙女欲春2 | 18禁黄网站禁片午夜丰满| 成人永久免费在线观看视频| 丰满的人妻完整版| 欧美乱妇无乱码| 中文字幕高清在线视频| 亚洲视频免费观看视频| 免费不卡黄色视频| 国产精品亚洲一级av第二区| 日韩三级视频一区二区三区| 法律面前人人平等表现在哪些方面| 亚洲欧美激情在线| 视频区欧美日本亚洲| 精品日产1卡2卡| 亚洲 欧美 日韩 在线 免费| 亚洲欧美一区二区三区黑人| 久久青草综合色| 亚洲精品一区av在线观看| 国产精品 国内视频| 啦啦啦观看免费观看视频高清 | 亚洲国产精品久久男人天堂| www.自偷自拍.com| a在线观看视频网站| 国产亚洲精品av在线| 国产又色又爽无遮挡免费看| 热99re8久久精品国产| 搞女人的毛片| 久久精品影院6| 人妻久久中文字幕网| 麻豆久久精品国产亚洲av| 一区福利在线观看| 丝袜美腿诱惑在线| 国产精品久久久久久亚洲av鲁大| 欧美黑人精品巨大| 91av网站免费观看| 国产av又大| 99久久综合精品五月天人人| 日韩大尺度精品在线看网址 | 村上凉子中文字幕在线| 精品国产美女av久久久久小说| 满18在线观看网站| 久久精品国产99精品国产亚洲性色 | 熟女少妇亚洲综合色aaa.| 又黄又粗又硬又大视频| 亚洲专区中文字幕在线| 99国产精品免费福利视频| 欧美成狂野欧美在线观看| 精品久久久久久成人av| 成人国产综合亚洲| 国产精品香港三级国产av潘金莲| 国产精品免费一区二区三区在线| 国产精华一区二区三区| 免费久久久久久久精品成人欧美视频| 久久伊人香网站| 午夜精品久久久久久毛片777| 久久精品成人免费网站| 久久久久久久久久久久大奶| 日日干狠狠操夜夜爽| 视频区欧美日本亚洲| 亚洲成人国产一区在线观看| 久久久水蜜桃国产精品网| av在线天堂中文字幕| 亚洲av第一区精品v没综合| 免费女性裸体啪啪无遮挡网站| 免费人成视频x8x8入口观看| bbb黄色大片| 国产精品美女特级片免费视频播放器 | 一二三四社区在线视频社区8| 一进一出抽搐gif免费好疼| 成人特级黄色片久久久久久久| 成年版毛片免费区| 国产成人欧美| 91字幕亚洲| 女人被躁到高潮嗷嗷叫费观| 黄色 视频免费看| 亚洲精品在线观看二区| 女人高潮潮喷娇喘18禁视频| 亚洲第一青青草原| 亚洲精品粉嫩美女一区| 亚洲五月婷婷丁香| 亚洲最大成人中文| 亚洲av日韩精品久久久久久密| 亚洲国产精品久久男人天堂| 最好的美女福利视频网| 别揉我奶头~嗯~啊~动态视频| 久久久久国产精品人妻aⅴ院| 岛国视频午夜一区免费看| 国产区一区二久久| 久久久久久免费高清国产稀缺| 久久中文字幕一级| 精品国内亚洲2022精品成人| 日韩欧美一区二区三区在线观看| 久久精品国产99精品国产亚洲性色 | 国产男靠女视频免费网站| 久久精品91无色码中文字幕| 精品久久久精品久久久| 国产精品二区激情视频| 亚洲avbb在线观看| 国产精品久久久久久亚洲av鲁大| 高清毛片免费观看视频网站| 亚洲va日本ⅴa欧美va伊人久久| 日本精品一区二区三区蜜桃| 熟女少妇亚洲综合色aaa.| 国产精品综合久久久久久久免费 | 高潮久久久久久久久久久不卡| 国产亚洲精品第一综合不卡| 男人操女人黄网站| 麻豆成人av在线观看| 亚洲性夜色夜夜综合| www.熟女人妻精品国产| 亚洲欧美日韩无卡精品| 无人区码免费观看不卡| 亚洲 欧美 日韩 在线 免费| 国产成人精品久久二区二区91| 可以在线观看毛片的网站| 免费少妇av软件| 99精品久久久久人妻精品| 国产精品1区2区在线观看.| 俄罗斯特黄特色一大片| 成人精品一区二区免费| 99国产精品免费福利视频| 欧美中文综合在线视频| 国产亚洲精品第一综合不卡| 国产国语露脸激情在线看| 又紧又爽又黄一区二区| 日本a在线网址| 人妻丰满熟妇av一区二区三区| 亚洲国产精品999在线| 别揉我奶头~嗯~啊~动态视频| 少妇 在线观看| 久久中文字幕人妻熟女| 麻豆一二三区av精品| 在线免费观看的www视频| 大码成人一级视频| av中文乱码字幕在线| 99在线视频只有这里精品首页| 午夜激情av网站| 久久久精品欧美日韩精品| 欧美久久黑人一区二区| 亚洲视频免费观看视频| 成人亚洲精品av一区二区| 欧美+亚洲+日韩+国产| 给我免费播放毛片高清在线观看| 一二三四在线观看免费中文在| 黄频高清免费视频| 老司机深夜福利视频在线观看| 亚洲精品粉嫩美女一区| 欧美最黄视频在线播放免费| 美女高潮喷水抽搐中文字幕| 91大片在线观看| 99re在线观看精品视频| 大码成人一级视频| 欧美在线黄色| 国产区一区二久久| 男女午夜视频在线观看| 性欧美人与动物交配| 日本免费一区二区三区高清不卡 | 亚洲av成人不卡在线观看播放网| 久久午夜亚洲精品久久| 久久精品91蜜桃| 啪啪无遮挡十八禁网站| 在线十欧美十亚洲十日本专区| 国产一区二区三区在线臀色熟女| 亚洲 欧美一区二区三区| 亚洲激情在线av| 国产精品久久电影中文字幕| 男女下面进入的视频免费午夜 | 亚洲第一电影网av| 在线观看一区二区三区| 一卡2卡三卡四卡精品乱码亚洲| 国产精品久久久久久亚洲av鲁大| 热99re8久久精品国产| 日本免费一区二区三区高清不卡 | 女人高潮潮喷娇喘18禁视频| 日本撒尿小便嘘嘘汇集6| 少妇粗大呻吟视频| 欧美乱妇无乱码| 亚洲专区国产一区二区| 视频在线观看一区二区三区| 国产精品98久久久久久宅男小说| 免费人成视频x8x8入口观看| 亚洲色图av天堂| 最近最新中文字幕大全免费视频| 亚洲自偷自拍图片 自拍| av有码第一页| 制服人妻中文乱码| 一区二区日韩欧美中文字幕| 国产精品久久电影中文字幕| 男女做爰动态图高潮gif福利片 | 天天躁狠狠躁夜夜躁狠狠躁| 亚洲一区高清亚洲精品| 国产野战对白在线观看| 两人在一起打扑克的视频| 日韩欧美国产一区二区入口| 日本撒尿小便嘘嘘汇集6| 国产av在哪里看| 精品国产一区二区久久| 久久欧美精品欧美久久欧美| 国内精品久久久久久久电影| 波多野结衣高清无吗| 国产麻豆成人av免费视频| 久久婷婷成人综合色麻豆| 黄片大片在线免费观看| 中文亚洲av片在线观看爽| 午夜福利在线观看吧| 精品福利观看| 日韩大码丰满熟妇| 亚洲三区欧美一区| 久久九九热精品免费| 久久香蕉精品热| av超薄肉色丝袜交足视频| 美女国产高潮福利片在线看| 久久亚洲真实| 国产精品98久久久久久宅男小说| 国产亚洲av嫩草精品影院| 亚洲成人免费电影在线观看| 成人手机av| 日本黄色视频三级网站网址| 成年人黄色毛片网站| 久久香蕉精品热| 国产成人影院久久av| 国产精品一区二区精品视频观看| 夜夜爽天天搞| 亚洲欧美激情在线| 后天国语完整版免费观看| 日韩免费av在线播放| 日本vs欧美在线观看视频| 神马国产精品三级电影在线观看 | 国产精华一区二区三区| 不卡一级毛片| 好男人在线观看高清免费视频 | 午夜福利,免费看| 最新美女视频免费是黄的| 美女高潮喷水抽搐中文字幕| 一区二区三区国产精品乱码| 色综合站精品国产| 精品卡一卡二卡四卡免费| 国产av精品麻豆| 欧美乱色亚洲激情| 国产成人精品久久二区二区91| 手机成人av网站| 老鸭窝网址在线观看| 女同久久另类99精品国产91| 成人手机av| 一区二区三区高清视频在线| 免费在线观看亚洲国产| 亚洲成人精品中文字幕电影| 国产欧美日韩一区二区三区在线| 中国美女看黄片| 亚洲成人免费电影在线观看| 在线观看日韩欧美| 亚洲精品一区av在线观看| 女人精品久久久久毛片| 午夜福利欧美成人| 日本三级黄在线观看| 国产精品久久久人人做人人爽| 国产av一区在线观看免费| 18禁国产床啪视频网站| 中文字幕久久专区| 久久精品人人爽人人爽视色| av电影中文网址| 亚洲九九香蕉| 女性被躁到高潮视频| 免费人成视频x8x8入口观看| 亚洲三区欧美一区| 免费在线观看黄色视频的| 午夜免费成人在线视频| 老司机午夜福利在线观看视频| 日本 av在线| 亚洲精品国产色婷婷电影| 欧美激情 高清一区二区三区| 91老司机精品| 欧美日韩中文字幕国产精品一区二区三区 | 久久精品91蜜桃| 99香蕉大伊视频| 亚洲 国产 在线| aaaaa片日本免费| 精品久久久久久,| 精品日产1卡2卡| 一区二区日韩欧美中文字幕| 高潮久久久久久久久久久不卡| 女性生殖器流出的白浆| 欧美绝顶高潮抽搐喷水| 一级片免费观看大全| 亚洲电影在线观看av| 中文字幕人妻熟女乱码| 午夜福利18| 在线播放国产精品三级| 亚洲久久久国产精品| 国内精品久久久久精免费| 露出奶头的视频| 天天一区二区日本电影三级 | 老汉色∧v一级毛片| 国产极品粉嫩免费观看在线| 亚洲第一av免费看| 久久影院123| 亚洲熟妇熟女久久| 无限看片的www在线观看| 日本 欧美在线| 亚洲 欧美 日韩 在线 免费| 精品一区二区三区四区五区乱码| 丁香欧美五月| 成人18禁高潮啪啪吃奶动态图| 日韩三级视频一区二区三区| 色播在线永久视频| 狠狠狠狠99中文字幕| 久久久久国产精品人妻aⅴ院| 亚洲人成电影观看| 纯流量卡能插随身wifi吗| 在线av久久热| 岛国视频午夜一区免费看| 在线国产一区二区在线| 久久人妻熟女aⅴ| 激情在线观看视频在线高清| 大香蕉久久成人网| 伊人久久大香线蕉亚洲五| 桃色一区二区三区在线观看| 中文字幕人成人乱码亚洲影| 国产av一区在线观看免费| 天堂动漫精品| 9191精品国产免费久久| 午夜亚洲福利在线播放| 亚洲熟女毛片儿| 麻豆av在线久日| 好男人电影高清在线观看| e午夜精品久久久久久久| 免费少妇av软件| 久久午夜综合久久蜜桃| 日韩三级视频一区二区三区| 亚洲一区二区三区不卡视频| 色播亚洲综合网| 久久香蕉精品热| 免费在线观看黄色视频的| 美女免费视频网站| 久久午夜综合久久蜜桃| 精品久久久久久久毛片微露脸| 九色亚洲精品在线播放| 久久午夜亚洲精品久久| 久久久国产欧美日韩av| 法律面前人人平等表现在哪些方面| 国产91精品成人一区二区三区| 日韩大尺度精品在线看网址 | 一级a爱视频在线免费观看| 99香蕉大伊视频| 国产又色又爽无遮挡免费看| 日韩大尺度精品在线看网址 | 国产99白浆流出| e午夜精品久久久久久久| 久久精品aⅴ一区二区三区四区| 中文字幕久久专区| 免费观看人在逋| 久久久久国产精品人妻aⅴ院| 国产精华一区二区三区| 午夜免费成人在线视频| 多毛熟女@视频| 亚洲精品美女久久av网站| 国产伦一二天堂av在线观看| 久久久久国内视频| 久久久久亚洲av毛片大全| 亚洲avbb在线观看| 国产午夜精品久久久久久| 搞女人的毛片| 午夜福利视频1000在线观看 | 国产激情久久老熟女| 亚洲成av片中文字幕在线观看| 91在线观看av| 亚洲国产欧美一区二区综合| 国产一级毛片七仙女欲春2 | 日日爽夜夜爽网站| 免费观看人在逋| 变态另类成人亚洲欧美熟女 | 天天躁狠狠躁夜夜躁狠狠躁| 黄片播放在线免费| 叶爱在线成人免费视频播放| 国产99白浆流出| 中出人妻视频一区二区| 亚洲av成人av| 午夜福利影视在线免费观看| 午夜福利在线观看吧| 久久欧美精品欧美久久欧美| 欧洲精品卡2卡3卡4卡5卡区| 久久精品影院6| 超碰成人久久| 如日韩欧美国产精品一区二区三区| 国产三级在线视频| 中文字幕最新亚洲高清| 叶爱在线成人免费视频播放| 精品一区二区三区视频在线观看免费| 女人被狂操c到高潮| 精品乱码久久久久久99久播| 午夜日韩欧美国产| 色播亚洲综合网| 亚洲专区中文字幕在线| 色尼玛亚洲综合影院| av有码第一页| 老熟妇乱子伦视频在线观看| 欧美色视频一区免费| 色在线成人网| 午夜免费鲁丝| 欧美成人性av电影在线观看| 麻豆久久精品国产亚洲av| 国产亚洲欧美精品永久| 国产三级黄色录像| 九色亚洲精品在线播放| 女人精品久久久久毛片| 天堂√8在线中文| 激情在线观看视频在线高清| aaaaa片日本免费| 99热只有精品国产| 成人永久免费在线观看视频| 亚洲色图 男人天堂 中文字幕| 日韩 欧美 亚洲 中文字幕| 免费在线观看日本一区| 国产成人精品无人区| 国产成人欧美| √禁漫天堂资源中文www| 亚洲av日韩精品久久久久久密| 日韩中文字幕欧美一区二区| 狂野欧美激情性xxxx| 久久久水蜜桃国产精品网| 国产精品一区二区精品视频观看| aaaaa片日本免费| 亚洲男人天堂网一区| 99国产精品一区二区三区| 麻豆国产av国片精品| 亚洲第一av免费看| 91精品国产国语对白视频| 亚洲精华国产精华精| 村上凉子中文字幕在线| 午夜免费鲁丝|