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

    不同泥沙粒徑下的分離鰓內(nèi)部流場三維數(shù)值模擬

    2014-07-02 01:38:44陶洪飛邱秀云蘇建江
    水利水電科技進展 2014年5期
    關(guān)鍵詞:水沙含沙量泥沙

    陶洪飛,邱秀云,李 巧,蘇建江

    (1.新疆農(nóng)業(yè)大學(xué)水利與土木工程學(xué)院,新疆烏魯木齊 830052; 2.國電新疆吉林臺水電開發(fā)有限公司,新疆伊犁 835716)

    不同泥沙粒徑下的分離鰓內(nèi)部流場三維數(shù)值模擬

    陶洪飛1,邱秀云1,李 巧1,蘇建江2

    (1.新疆農(nóng)業(yè)大學(xué)水利與土木工程學(xué)院,新疆烏魯木齊 830052; 2.國電新疆吉林臺水電開發(fā)有限公司,新疆伊犁 835716)

    為研究泥沙粒徑對分離鰓的速度場及水沙分離效果的影響,采用層流模型和歐拉模型,運用PC-SIMPLE算法,對分離鰓內(nèi)的水沙兩相流流場進行了三維數(shù)值模擬。計算結(jié)果表明:當泥沙粒徑d=0.0001 mm時,分離鰓內(nèi)的流場與其他泥沙粒徑時有所不同;當0.005 mm≤d≤0.035 mm時,分離鰓內(nèi)不同粒徑泥沙水流的流場相同;當0.015 mm≤d≤0.035 mm時,分離鰓處理泥沙效果最好,且泥沙粒徑越大,去除率越高,水沙分離速度就越快。

    泥沙粒徑;分離鰓;歐拉模型;去除率;數(shù)值模擬

    我國是世界上河流泥沙問題最嚴重的國家之一。據(jù)觀測資料統(tǒng)計,我國西北和華北地區(qū)的大部分河流都具有含沙量高且泥沙顆粒細的特點,這給工業(yè)用水、生活用水以及多沙河流沿岸居住分散的農(nóng)牧民群眾用水造成極大的不便,并對灌區(qū)的滴、噴灌設(shè)備帶來很大的危害,灌溉時常常會發(fā)生噴頭或滴灌頭流道堵塞的問題,嚴重制約了噴、滴灌等先進節(jié)水灌溉技術(shù)的推廣和應(yīng)用。

    分離鰓是新疆農(nóng)業(yè)大學(xué)水利與土木工程學(xué)院水力學(xué)教研室課題組研制的一種新型的水沙兩相流分離裝置,與傳統(tǒng)水沙分離方法相比,分離鰓具有分離速度快、結(jié)構(gòu)簡單、造價低廉及運行費低等特點,同時分離出來的清水和排出的泥沙對人體及周圍環(huán)境無害。因此分離鰓對解決分散居住的農(nóng)牧民生活用水問題具有積極意義,同時分離鰓作為灌區(qū)泥沙前期處理設(shè)施,可為灌區(qū)滴、噴灌提供優(yōu)質(zhì)水源,促進節(jié)水灌溉技術(shù)在多沙河流地區(qū)的推廣應(yīng)用。根據(jù)文獻[1-12]可知,目前分離鰓已做了如下研究:通過系列模型試驗和理論分析,探明了分離鰓的分離機理以及分離鰓內(nèi)水沙分離的水流結(jié)構(gòu),并對其結(jié)構(gòu)進行了優(yōu)化;運用粒子圖像測速技術(shù)測試了分離鰓內(nèi)部流場,并利用流體軟件Fluent對分離鰓內(nèi)部的水沙兩相流流場進行了數(shù)值模擬,獲得了分離鰓內(nèi)部的流場分布;在上述研究的基礎(chǔ)上提出了分離鰓應(yīng)用于實際工程的集成形式,并對此形式進行了室內(nèi)模型集成試驗與室外大比尺試驗。影響分離鰓的水沙分離效果的因素有多種,其中泥沙顆粒的大小對其沉淀效率不容忽視。文獻[13]通過流體軟件Fluent模擬了分離鰓內(nèi)水沙兩相流流場,對計算結(jié)果與粒子圖像測速(particle Image velocimetry,PIV)的測試結(jié)果進行了定性與定量的比較,發(fā)現(xiàn)歐拉模型適合模擬分離鰓內(nèi)的水沙兩相流流場,故本文采用歐拉模型模擬靜水中不同泥沙粒徑下分離鰓內(nèi)部的水沙兩相流流場,分析和討論泥沙粒徑對分離鰓內(nèi)部流場及沉淀效果的影響,為分離鰓的推廣應(yīng)用奠定基礎(chǔ)。

    1 模型建立

    1.1 分離鰓結(jié)構(gòu)

    分離鰓由鰓片與鰓管構(gòu)成,靜水沉降時,分離鰓的上端為開口狀,下端呈封閉狀,該裝置的典型結(jié)構(gòu)特征見圖1。分離鰓的主要特征尺寸有鰓片間距d、鰓片與鰓管寬度側(cè)壁構(gòu)成的α傾斜角、鰓片與鰓管長度側(cè)壁構(gòu)成的β傾斜角、清水流上升通道寬度e及泥沙流下降通道寬度f。采用的分離鰓特征尺寸為:a=10 cm,b=5 cm,c=40 cm,d=5 cm,e=f=1 cm,采用優(yōu)化后的傾斜角α=60°和β=45°。

    圖1 分離鰓結(jié)構(gòu)示意圖

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

    考慮本研究源項為零,歐拉模型的質(zhì)量守恒方程為

    式中:▽是那勃勒算子;vq為第q相的流速;φq為第q相的體積分數(shù);ρq為第q相的物理密度。

    數(shù)值模擬中升力相對于曳力是不重要的,因此動量守恒方程中不包含升力;泥沙(第q相)的密度遠大于水(第p相)的密度,虛擬質(zhì)量力可忽略不計,則動量守恒方程為

    其中

    式中:t為時間;vq為第q相的流速矢量;p′為水與沙共同承擔的壓力;τq為第q相的壓力應(yīng)變張量;Rpq為水與沙之間的作用力;Fq為第q相的外部體積力張量;λq與μq分別為第q相的體積黏度和剪切黏度;I為單位張量。

    Rpq取決于壓力、摩擦、內(nèi)聚力等的影響,其計算公式如下:

    其中

    式中:Kqp為相間動量交換系數(shù);φp為第p相的體積分數(shù);ρp為第p相的物理密度;C為曳力函數(shù);Req為第q相的雷諾數(shù);vr,q為與泥沙相關(guān)的末端流速;vp為第p相的流速矢量;dq為第q相的直徑;A和B為系數(shù);μp為p相的剪切黏度。

    羅菲[9]利用埃施理論說明泥沙顆粒在分離鰓靜水沉降過程中,水沙之間不會相互混摻,水流較穩(wěn)定且流態(tài)為層流,故選擇層流模型與歐拉模型來模擬分離鰓中的流場。層流模型的方程如下:

    式中:v為流速矢量;p″為流體靜壓;μ為流體的動力黏滯系數(shù)。

    1.3 計算方法

    分離鰓的計算區(qū)域和基本控制方程的離散采用有限體積法[13]。考慮減小數(shù)值模擬的計算工作量和保證計算精度,最終采用優(yōu)化后的四面體非結(jié)構(gòu)性網(wǎng)格,網(wǎng)格數(shù)共計81714個。流相的離散格式采用具有良好遷移特性和較好穩(wěn)定性的一階迎風(fēng)格式。對歐拉多相流計算,采用PC-SIMPLE算法來計算離散后的線性代數(shù)方程組,各方程計算精度均為0.001。

    1.4 邊界條件

    a.分離鰓上端邊界條件:分離鰓中泥沙顆粒作靜水沉降運動,分離鰓上端為一水平的自由表面,分離鰓內(nèi)的水頭大小不變,且水面沿y軸負方向無波動現(xiàn)象,可以認為流速及其他各個變量沿y方向的梯度都為零,故可采用“剛蓋假定”。在自由表面上采用對稱邊界條件處理。

    b.壁面邊界條件:分離鰓的固體邊壁包括鰓片、鰓管下端和鰓管外圍邊壁,其邊界條件采用無滑移邊界條件。

    1.5 模型參數(shù)及運行效果

    初始化時,設(shè)置時間步長為0.0001 s,求解時間步數(shù)為1×106,每個時間步數(shù)最多的迭代次數(shù)為200。時間步數(shù)為0~1000時,各方程的殘差線波動很大,1000步以后波動很小;當時間步達到20 000時,殘差線都在數(shù)值0.001以下,說明各方程的計算精度均達到0.001,各方程已經(jīng)收斂;隨著進一步的求解計算,收斂效果更好,此時可適當增大時間步長,以減少計算時間。

    2 監(jiān)測面與工況設(shè)置

    為得到不同泥沙粒徑下分離鰓中的泥沙分布特性,在數(shù)值模擬計算過程中應(yīng)監(jiān)測y方向(高度方向)的平均含沙量與時間的關(guān)系,以便對數(shù)值計算結(jié)果進行分析和對比,具體的斷面位置見圖2。

    圖2 不同位置的監(jiān)測面

    分離鰓的計算區(qū)域包括水和沙兩種介質(zhì),將水定義為主相,泥沙顆粒為第2相,假定泥沙顆粒為球形。研究泥沙粒徑對分離鰓內(nèi)部流場及沉淀效果的影響時,水的密度ρw=998.2 kg/m3,泥沙顆粒的密度ρs=2650 kg/m3,渾水含沙量S=10 kg/m3,它們均為固定值,已有的試驗表明分離鰓在處理細顆粒泥沙時能體現(xiàn)出其優(yōu)越性,故按表1中的不同工況進行數(shù)值模擬計算。

    表1 泥沙粒徑的工況設(shè)置

    3 數(shù)值模擬結(jié)果與分析

    3.1 鰓片表面的流速矢量分布

    限于篇幅,本文只給出了工況1、4、6的歐拉模型數(shù)值模擬結(jié)果(工況2、3、5的計算結(jié)果同工況1、4),即給出這3種工況下鰓片1(靠近分離鰓的自由表面)、鰓片3(位于分離鰓中間)、鰓片5(靠近分離鰓底端)和z方向各斷面的流速矢量分布,通過對比分析確定泥沙粒徑對分離鰓內(nèi)流速分布規(guī)律的影響。

    當計算沉降時間為400 s時,不同泥沙粒徑下鰓片1、3、5的上、下表面流速矢量分布如圖3~8所示。當d=0.035 mm和d=0.01 mm時,同一鰓片的上表面流速分布規(guī)律一致,同樣鰓片下表面的流速分布規(guī)律也一致;當d=0.0001mm時,鰓片上、下表面的流速分布規(guī)律與上述兩種泥沙粒徑下的流速分布規(guī)律有所不同。相同泥沙粒徑下不同鰓片的上表面流速分布規(guī)律有所不同,下表面流速分布規(guī)律也不同,這與鰓片在分離鰓中的位置有關(guān)。

    從圖3可知,當d=0.035 mm和d=0.01 mm時,鰓片1上表面的水流先以傾斜角α沿鰓片短邊向下滑動,當滑至鰓片長邊時以傾斜角β貼壁向下滑動;當d=0.0001 mm時,鰓片上表面的水流完全呈放射狀。從圖4和圖5可看出,當d=0.035 mm和d=0.01 mm時,鰓片3與鰓片5上表面的流速分布規(guī)律同鰓片1上表面,但因泥沙通道中下落水流的影響,從而這兩個鰓片上表面右上角的流速呈向內(nèi)放射狀;當d=0.000 1 mm時,鰓片3與鰓片5上表面的水流在鰓片高端區(qū)域形成順時針漩渦,低端區(qū)域則有部分水流以傾斜角β向下滑動。

    從圖6和圖7可知,當d=0.035 mm和d= 0.01 mm時,鰓片1與鰓片3下表面的水流先以傾斜角α沿鰓片短邊向上流動,當流至鰓片長邊時以傾斜角β貼壁向上流動,但因清水通道中上升水流的影響,從而這兩個鰓片下表面左下角的流速呈向內(nèi)放射狀;當d=0.000 1 mm時,鰓片1與鰓片3下表面的水流在鰓片低端區(qū)域形成順時針的漩渦,高端區(qū)域則有部分水流以傾斜角β向上流動。從圖8中可看出,當d=0.035mm和d=0.01mm時,鰓片5下表面的流速分布規(guī)律與鰓片1、鰓片3的下表面相同,但因鰓片5靠近分離鰓的底端,清水通道中上升水流對其流場沒有影響,故左下角不存在流速呈向內(nèi)放射狀的現(xiàn)象;當d=0.000 1 mm時,鰓片下表面的水流完全呈放射狀。

    圖3 鰓片1上表面流速矢量分布

    圖4 鰓片3上表面流速矢量分布

    圖5 鰓片5上表面流速矢量分布

    圖6 鰓片1下表面流速矢量分布

    圖7 鰓片3下表面流速矢量分布

    圖8 鰓片5下表面流速矢量分布

    圖9 泥沙顆粒受力狀態(tài)

    泥沙顆粒沉降到鰓片長邊時,其受力狀態(tài)如圖9所示。水流阻力FD與摩阻力Ff與泥沙運動方向相反;上舉力FL與斜面支撐力FR與斜面垂直;泥沙顆粒的有效重力W豎直向下。采用泥沙顆?;破饎幽J?當滿足0.707W>FD+Ff時,泥沙顆粒便在鰓片的長邊滑動;當0.01mm≤d≤0.035 mm時,泥沙有效重力為8.48×10-12~3.63×10-10N,鰓片上表面的泥沙顆粒會在有效重力作用下克服阻力,在鰓片上表面滑動至泥沙通道中,對于鰓片2、3、4、5的上表面而言,因同時受到泥沙通道中水流的沖擊影響,從而流場與鰓片1上表面的有所不同;當d=0.0001 mm時,泥沙有效重力為8.48×10-18N,其所受有效重力的大小與其他泥沙粒徑不在同一個量級上,故受流場的影響最大,鰓片上表面的泥沙顆粒,由于有效重力很小,此時上舉力FL占主導(dǎo)地位,同時泥沙通道及清水通道中的水流對小粒徑的運動影響很大,才導(dǎo)致該粒徑下鰓片上表面的流速矢量分布規(guī)律與其他泥沙顆粒不同。根據(jù)連續(xù)原理[2,9]可知,不同泥沙粒徑下鰓片下表面的流場與鰓片上表面的相反。

    圖10 z=0.5 cm斷面流速矢量分布

    圖11 z=2.5 cm斷面流速矢量分布

    圖12 z=4.5 cm斷面流速矢量分布

    3.2 z方向各斷面的流速矢量分布

    當計算沉降時間為400 s時,不同泥沙粒徑下z方向各斷面的流速矢量分布如圖10~12所示。當d=0.035 mm和d=0.01 mm時,z=0.5 cm斷面(清水通道斷面)、z=2.5 cm斷面(分離鰓寬方向的中間斷面)、z=4.5 cm斷面(泥沙下沉通道斷面)相應(yīng)的流速分布規(guī)律均相同。由圖10可知,當d= 0.035 mm和d=0.01mm時,各鰓片下表面的水流沿著各自鰓片下表面向上運動,至鰓片高端匯入清水通道,再一同向上運動;由圖12可知,當d=0.035 mm和d=0.01 mm時,各鰓片上表面的水流沿著各自鰓片上表面向下運動,至鰓片低端匯入泥沙通道,再一同向下運動,這與物理試驗中所觀測到的現(xiàn)象一致[1-9],即在鰓片之間形成了橫向異重流,在整個系統(tǒng)中形成了水沙沿分離鰓邊璧的垂向異重流;而這兩種現(xiàn)象在圖11中當d=0.035 mm和d=0.01 mm得到了完整的體現(xiàn),鰓片上下表面的水流、泥沙通道中的水流及清水通道中的水流互不干擾,各自按著相應(yīng)的運動軌跡流動,從而達到水沙分離的目的。當d=0.0001 mm時,z方向各斷面的流速矢量分布規(guī)律與其他兩個泥沙粒徑的有所不同。由圖10可知,當d=0.0001 mm時,鰓片間大部分的水流斜向下運動,只有少部分水流沿著鰓片下表面向上運動后進入清水通道,分離鰓的右上角與左下角形成了順時針的漩渦;由圖12可知,當d=0.0001 mm時,鰓片間大部分的水流斜向上運動,只有少部分水流沿著鰓片上表面向下運動后進入泥沙通道,分離鰓的右上角與左下角也形成了順時針的漩渦,而這兩種現(xiàn)象在圖11中當d=0.0001 mm時完全體現(xiàn),相鄰鰓片間的水流完全混合,并形成了4個順時針的封閉的漩渦,同時在分離鰓的右上角與左下角也形成了順時針的漩渦,此時水沙不能從分離鰓中分離出來。

    3.3 y方向各斷面的平均含沙量分布

    圖13為不同泥沙粒徑下分離鰓y方向(高度方向)各監(jiān)測斷面的平均含沙量的分布特性。由圖13可知:①當d=0.0001 mm時,平均含沙量的分布規(guī)律與0.005 mm≤d≤0.035 mm時有所不同,平均含沙量隨時間的增大而保持不變,即為一條平行于時間軸的直線。②y=5 cm與y=15 cm斷面靠近分離鰓的底端,當0.005mm≤d≤0.035mm時,平均含沙量隨時間的變化關(guān)系與y=25 cm、y=35 cm兩個段面有所不同,先隨時間的增大而增大,后隨時間的增大而減小,y值越小(越靠近分離鰓底端)水沙分離效果越明顯,這是因為上面的泥沙沉降到底部所引起的;但不同泥沙粒徑下平均含沙量隨時間的增大或減小的幅度不同,泥沙粒徑越大,平均含沙量隨時間的變化幅度就越大,說明泥沙沉降的也就越快,從而分離鰓的沉淀效率也就越高。當0.005 mm≤d≤0.035 mm時,y=25 cm與y=35 cm兩個斷面上的平均含沙量隨時間變化規(guī)律相同,均隨著時間的延長而減小,但不同泥沙粒徑下平均含沙量減小的幅度不同,泥沙粒徑越大,平均含沙量隨時間的減小幅度就越大,同一時間下,平均含沙量就越小,即水沙分離效果就越好,當時間為1000s、d=0.005mm時平均含沙量為9.18 kg/m3,而當d=0.035 mm時,平均含沙量幾乎為零,如圖13(d)所示。

    圖13 y方向各斷面的平均含沙量分布

    3.4 不同泥沙粒徑的去除率分布

    圖14為某一時間y方向不同斷面、不同粒徑的去除率分布,去除率是判斷分離鰓水沙分離效果好壞的一個重要參數(shù),其表達式為

    式中:η為去除率,%;S1為初始時刻某斷面的渾水含沙量,kg/m3;S2為沉降到某時刻該斷面的平均含沙量,kg/m3。當y=15 cm斷面的去除率不小于70%時,可認為分離鰓的水沙分離效果最好,說明該斷面以上的其他斷面去除率必定大于70%,甚至靠近分離鰓頂端的斷面可達到100%,即y=15 cm斷面以上的高含沙渾水已經(jīng)成為低含沙的渾水;反之水沙分離效果則一般。

    從圖14可知:①y=25 cm與y=35 cm兩個斷面的去除率隨著泥沙粒徑的增大而增大,當去除率達到100%時,表明分離鰓中的水沙分離完全;y= 5 cm與y=15 cm兩個斷面因靠近分離鰓的下端,去除率先隨泥沙粒徑的增大而保持不變,都為零,后隨泥沙粒徑的增大而增大,當水沙分離完全時,又保持不變。②當0.0001 mm≤d<0.005 mm時,4個不同斷面的去除率都在23%以下,且d=0.000 1 mm時的去除率全為零,說明分離鰓不能處理這種粒徑的泥沙;當0.005 mm≤d<0.015 mm時,在同一粒徑下,y越大去除率越高,即越靠近分離鰓的上端,其水沙分離效果就越好,就越容易獲得清水,分離鰓處理這樣的泥沙效果一般;當0.015 mm≤d≤0.035 mm時,在同一粒徑下,y越大去除率就越高或為100%,且在同一斷面上相對于其他泥沙粒徑去除率要大很多,說明分離鰓處理這類泥沙效果最好。

    圖14 y方向各斷面不同泥沙粒徑的去除率分布

    4 結(jié) 論

    a.在同一鰓片上,當0.005 mm≤d≤0.035 mm時,鰓片上表面的流速分布規(guī)律一致,同樣鰓片下表面的流速分布規(guī)律也一致;當d=0.0001 mm時,鰓片上、下表面的流速分布規(guī)律則與0.005 mm≤d≤0.035 mm時有所不同。當0.000 1 mm≤d≤0.035 mm時,不同鰓片的流速分布規(guī)律有所不同,這與鰓片在分離鰓中的位置有關(guān)。

    b.當0.005 mm≤d≤0.035 mm時,z=0.5 cm斷面、z=2.5 cm斷面和z=4.5 cm斷面相應(yīng)的流速分布規(guī)律均相同;當d=0.000 1 mm時,z方向這3個斷面的流速分布規(guī)律則與之不同。

    c.當0.015 mm≤d≤0.035 mm時,分離鰓處理泥沙效果最好。

    [1]邱秀云,龔守遠,嚴躍成,等.一種新型水沙分離裝置的研究[J].新疆農(nóng)業(yè)大學(xué)學(xué)報,2007,30(1):68-70. (QIU Xiuyun,GONG Shouyuan,YAN Yaocheng,et al. The research on a new separator device for water-sediment [J].Journal of Xinjiang Agricultural University,2007,30 (1):68-70.(in Chinese))

    [2]朱超,邱秀云,嚴躍成.垂向異重流式水沙分離鰓水沙分離機理淺析[J].水利水電科技進展,2009,29(5): 20-23.(ZHU Chao,QIU Xiuyun,YAN Yaocheng.Analysis of mechanism of vertical-component density flow watersediment separation device[J].Advances in Science and Technology of Water Resources,2009,29(5):20-23.(in Chinese))

    [3]陶洪飛,邱秀云,何照青,等.含沙量對分離鰓的水沙分離影響試驗[J].水電能源科學(xué),2013,31(7):93-95. (TAO Hongfei,QIU Xiuyun,HE Zhaoqing,et al.Sediment concentration on water-sediment separation effect test of separation device[J].Water Resources and Power,2013, 31(7):93-95.(in Chinese))

    [4]朱超,邱秀云,劉艷.垂向異重流式水沙分離鰓鰓片型式對水沙分離的影響研究[J].新疆農(nóng)業(yè)大學(xué)學(xué)報, 2008,31(6):72-75.(ZHU Chao,QIU Xiuyun,LIU Yan. Effects of gill-piece pattern of vertical-component density flow water-sediment separation device on water sand separation[J].JournalofXinjiangAgricultural University,2008,31(6):72-75.(in Chinese))

    [5]陶洪飛,邱秀云,李衛(wèi)鵬,等.分離鰓最優(yōu)鰓片間距的靜水沉降試驗[J].水利水電科技進展,2013,33(4):49-53.(TAO Hongfei,QIU Xiuyun,LI Weipeng,et al.The static water settlement test of optimal gill piece spacing in gill piece separation device[J].Advances in Science and Technology of Water Resources,2013,33(4):49-53.(in Chinese))

    [6]張翔,邱秀云,李琳,等.垂向異重流式分離鰓在動水中的集成試驗研究[J].新疆農(nóng)業(yè)大學(xué)學(xué)報,2010,33 (4):359-362.(ZHANG Xiang,QIU Xiuyun,LI Lin,et al.A study on integrated experiment in dynamic water of separation device with vertical density flow[J].Journal of Xinjiang Agricultural University,2010,33(4):359-362. (in Chinese))

    [7]張翔,邱秀云,李琳,等.垂向異重流式分離鰓在靜水中的集成試驗[J].水 利 水 運 工 程 學(xué) 報,2011(2): 89-93.(ZHANG Xiang,QIU Xiuyun,LI Lin,et al. Integration experiments of vertical-component density flow water-sediment separation device in still water.Hydro-Science and Engineering,2011(2):89-93.(in Chinese))

    [8]嚴躍成,邱秀云,張翔,等.兩相流分離鰓泥沙運動軌跡及加速滑移的力學(xué)分析[J].水利水電科技進展,2011, 31(5):27-29.(YAN Yaocheng,QIU Xiuyun,ZHANG Xiang,et al.Mechanical analysis of sediment motion path and accelerated slip on two-phase flow separation device [J].Advances in Science and Technology of Water Resources,2011,31(5):27-29.(in Chinese))

    [9]羅菲.兩相流分離鰓的水沙流場數(shù)值模擬[D].新疆:新疆農(nóng)業(yè)大學(xué),2011.

    [10]趙麗娜,邱秀云,陳鉑,等.兩相流水沙分離鰓流場的PIV測試及分析[J].新疆農(nóng)業(yè)大學(xué)學(xué)報,2012,35(2): 144-148.(ZHAO LIna,QIU Xiuyun,CHEN Bo,et al.Test and analysis on flow field in a two-phase flow watersediment separation device by piv techniques[J].Journal of Xinjiang Agricultural University,2012,35(2):144-148.(in Chinese))

    [11]陶洪飛,邱秀云,郝永志,等.兩種新型分離鰓沉淀池沉淀試驗及沉降效果研究[J].水電能源科學(xué),2012,30 (8):132-135.(TAOHongfei,QIUXiuyun,HAO Yongzhi,etal.Experimentalresearchoftwonew sedimentation tanks[J].water resources and power,2012, 30(8):132-135.(in Chinese))

    [12]陶洪飛,邱秀云,徐志宗,等.敞開式分離鰓沉淀池試驗[J].水利水電科技進展,2012,32(5):24-28.(TAO Hongfei,QIU Xiuyun,XU Zhizong,et al.An experiment of separation deviceinopensedimentationtank[J]. Advances in Science and Technology of Water Resources, 2012,32(5):24-28.(in Chinese))

    [13]陶洪飛,邱秀云,趙麗娜,等.水沙分離鰓內(nèi)部流場數(shù)值模擬與驗證[J].農(nóng)業(yè)工程學(xué)報,2013,29(17):38-46. (TAO Hongfei,QIU Xiuyun,ZHAO Lina,et al.Numerical simulation and PIV display of internal flow field in gillpiece separation device[J].Transactions of the Chinese Society of Agricultural Engineering,2013,29(17):38-46. (in Chinese))

    Three-dimensional numerical simulation of internal flow field in gill-sheet separation device under different sediment grain size//

    TAO Hongfei1,QIU Xiuyun1,LI Qiao1,SU Jianjiang2
    (1.College of Water Conservancy and Civil Engineering,Xinjiang AgriculturalUniversity,Urumqi 830052,China;2.GuodianXinjiangJilintaiHydropower Development Co.,Ltd.,Yili 835716,China)

    Three-dimensional numerical simulation adopting the laminar model and Euler model were conducted to the flow field of water-sediment two-phase in gill-sheet separation device to study how the sediment grain size effects upon velocity field and water-sediment separation by using PC-SIMPLE algorithm.The results show that the velocity flow is the same in the gill-sheet separation device when sediment grain size ranges from 0.005mm to 0.035mm and when the it is different the sediment grain size is 0.000 1 mm.The best performance is observed when gill-sheet separation device handles sediment grain size that ranges from 0.015mm to 0.035 mm.Additionally,the greater sediment grain size is,the higher removal is. Meanwhile the faster the water-sediment separation is.

    sediment grain size;gill-sheet separation device;Eulerian model;removal;numerical simulation

    TV14;TV131.3

    :A

    :1006-7647(2014)05-0017-07

    10.3880/j.issn.1006-7647.2014.05.004

    2013-0627 編輯:周紅梅)

    國家自然科學(xué)基金(50969009);國家星火計劃(2012GA890002);新疆水利水電工程重點學(xué)科基金(xjslgczdxk20101202)

    陶洪飛(1987—),男,四川南充人,博士研究生,主要從事計算水力學(xué)及河流泥沙工程研究。E-mail:304276290@qq.com

    邱秀云(1955—),女,江蘇泗陽人,教授,主要從事水力學(xué)及河流動力學(xué)研究。E-mail:wlmqqxy@sina.com

    猜你喜歡
    水沙含沙量泥沙
    泥沙做的父親
    0.6 H 層含沙量與垂線平均含沙量代表性探討
    新疆多泥沙河流水庫泥沙處理措施
    大型水利樞紐下游水沙變異特征
    土壤團聚體對泥沙沉降速度的影響
    山區(qū)河流上下雙丁壩回流區(qū)水沙特性淺探
    江西建材(2018年1期)2018-04-04 05:26:28
    泥沙滅火
    兒童繪本(2015年2期)2015-05-25 18:10:15
    羅源灣海洋傾倒區(qū)拋泥過程含沙量增量數(shù)值模擬
    懸移質(zhì)含沙量垂線分布
    水道港口(2014年1期)2014-04-27 14:14:35
    長江上中下游河道水沙特征和水沙關(guān)系
    亚洲精品,欧美精品| 黑人欧美特级aaaaaa片| 欧美精品人与动牲交sv欧美| 高清视频免费观看一区二区| 国产男人的电影天堂91| 欧美成人精品欧美一级黄| 免费在线观看黄色视频的| 波野结衣二区三区在线| 亚洲第一av免费看| 久久人人爽人人片av| 毛片一级片免费看久久久久| 各种免费的搞黄视频| 电影成人av| 看免费成人av毛片| videos熟女内射| 侵犯人妻中文字幕一二三四区| 久久久久久久精品精品| 亚洲美女视频黄频| 男男h啪啪无遮挡| 桃花免费在线播放| 国产一区二区三区av在线| 大陆偷拍与自拍| 19禁男女啪啪无遮挡网站| 建设人人有责人人尽责人人享有的| 亚洲 欧美一区二区三区| 日韩成人av中文字幕在线观看| 汤姆久久久久久久影院中文字幕| 黑丝袜美女国产一区| kizo精华| 午夜影院在线不卡| 天堂俺去俺来也www色官网| 国产日韩欧美视频二区| 亚洲欧美一区二区三区久久| 搡老乐熟女国产| bbb黄色大片| 亚洲国产欧美一区二区综合| 国产精品.久久久| 在线观看www视频免费| 欧美乱码精品一区二区三区| 91精品国产国语对白视频| 免费黄频网站在线观看国产| 成年人午夜在线观看视频| av国产精品久久久久影院| 亚洲国产欧美在线一区| 免费观看a级毛片全部| 国产免费福利视频在线观看| 精品免费久久久久久久清纯 | 99精品久久久久人妻精品| 国产伦理片在线播放av一区| 波多野结衣av一区二区av| 伊人久久大香线蕉亚洲五| 久久97久久精品| 制服诱惑二区| 亚洲一码二码三码区别大吗| 岛国毛片在线播放| 亚洲国产日韩一区二区| 天天躁狠狠躁夜夜躁狠狠躁| 久久久久久久国产电影| 一二三四在线观看免费中文在| 老司机深夜福利视频在线观看 | 国产老妇伦熟女老妇高清| 熟女少妇亚洲综合色aaa.| 国产一区二区三区av在线| 久久久久久久久久久久大奶| 最近中文字幕2019免费版| 久久久久精品人妻al黑| 国产av一区二区精品久久| 涩涩av久久男人的天堂| 麻豆乱淫一区二区| 最近手机中文字幕大全| 久久青草综合色| 男女国产视频网站| av网站在线播放免费| 国产精品亚洲av一区麻豆 | 婷婷色综合www| 久久精品亚洲熟妇少妇任你| 午夜日韩欧美国产| 亚洲国产精品成人久久小说| 丰满饥渴人妻一区二区三| 一边摸一边抽搐一进一出视频| 天天躁日日躁夜夜躁夜夜| 一二三四中文在线观看免费高清| av天堂久久9| 国产高清国产精品国产三级| 又粗又硬又长又爽又黄的视频| 精品一区二区三区四区五区乱码 | 色视频在线一区二区三区| 亚洲av日韩在线播放| 最近中文字幕2019免费版| 你懂的网址亚洲精品在线观看| 夫妻性生交免费视频一级片| 国产精品.久久久| 亚洲男人天堂网一区| 久久久精品免费免费高清| 十八禁高潮呻吟视频| 欧美成人精品欧美一级黄| 韩国av在线不卡| 国产色婷婷99| 韩国av在线不卡| 亚洲欧美中文字幕日韩二区| 最近手机中文字幕大全| 热re99久久国产66热| 制服人妻中文乱码| 午夜久久久在线观看| 日本一区二区免费在线视频| 观看av在线不卡| 亚洲色图 男人天堂 中文字幕| 又大又爽又粗| 电影成人av| 视频在线观看一区二区三区| 在线看a的网站| 亚洲欧美精品自产自拍| 欧美日韩亚洲综合一区二区三区_| 啦啦啦 在线观看视频| 国产伦理片在线播放av一区| 热99久久久久精品小说推荐| 国产xxxxx性猛交| 亚洲精品一区蜜桃| 中文字幕亚洲精品专区| 久久久久久久国产电影| 熟妇人妻不卡中文字幕| 亚洲精品国产av蜜桃| 18禁动态无遮挡网站| 亚洲精品国产av蜜桃| 999久久久国产精品视频| 免费黄网站久久成人精品| 18禁动态无遮挡网站| www.自偷自拍.com| 亚洲伊人色综图| 涩涩av久久男人的天堂| av在线老鸭窝| 欧美av亚洲av综合av国产av | 欧美精品人与动牲交sv欧美| 一级,二级,三级黄色视频| 国产成人精品久久二区二区91 | 午夜老司机福利片| 少妇被粗大猛烈的视频| 美女主播在线视频| 一区在线观看完整版| 两性夫妻黄色片| 国产黄色免费在线视频| 亚洲av欧美aⅴ国产| 亚洲av欧美aⅴ国产| 欧美精品人与动牲交sv欧美| 成人手机av| 在线观看免费午夜福利视频| av又黄又爽大尺度在线免费看| 别揉我奶头~嗯~啊~动态视频 | 一边摸一边抽搐一进一出视频| 欧美日韩av久久| 999精品在线视频| 熟妇人妻不卡中文字幕| 国产精品 欧美亚洲| 中文字幕色久视频| 国产成人免费无遮挡视频| 午夜91福利影院| 老司机影院成人| 久久人人爽人人片av| 午夜精品国产一区二区电影| 免费黄网站久久成人精品| 一个人免费看片子| 秋霞伦理黄片| 久久狼人影院| 国产又色又爽无遮挡免| 国产一区二区三区综合在线观看| 久久人妻熟女aⅴ| 日本色播在线视频| 少妇人妻久久综合中文| 中国三级夫妇交换| 午夜久久久在线观看| 男女国产视频网站| 青青草视频在线视频观看| 久久 成人 亚洲| 男男h啪啪无遮挡| 亚洲一区中文字幕在线| 黄片播放在线免费| 好男人视频免费观看在线| 2021少妇久久久久久久久久久| 亚洲欧美日韩另类电影网站| 午夜激情av网站| 性色av一级| 日本黄色日本黄色录像| 老汉色∧v一级毛片| 国产亚洲av高清不卡| 老汉色∧v一级毛片| 精品亚洲成a人片在线观看| 在线观看www视频免费| 如何舔出高潮| 亚洲精品日本国产第一区| 国产黄色视频一区二区在线观看| 天堂俺去俺来也www色官网| 在线观看一区二区三区激情| 19禁男女啪啪无遮挡网站| 国产精品国产三级专区第一集| 日本wwww免费看| 一边亲一边摸免费视频| 免费观看性生交大片5| 老司机亚洲免费影院| 免费黄网站久久成人精品| 两个人免费观看高清视频| 精品久久蜜臀av无| 久久精品国产亚洲av高清一级| 麻豆av在线久日| 亚洲精品,欧美精品| 日韩中文字幕视频在线看片| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲av中文av极速乱| 久久女婷五月综合色啪小说| 欧美日韩综合久久久久久| 国产成人精品无人区| 亚洲国产欧美一区二区综合| 久久精品亚洲av国产电影网| 成人国产麻豆网| 亚洲伊人久久精品综合| 少妇人妻久久综合中文| 日韩av免费高清视频| 国产视频首页在线观看| 最近中文字幕高清免费大全6| 人人妻人人添人人爽欧美一区卜| 久久久久网色| 女人被躁到高潮嗷嗷叫费观| 男人添女人高潮全过程视频| 99国产精品免费福利视频| 亚洲国产精品一区二区三区在线| 成人国产麻豆网| 欧美精品av麻豆av| av国产精品久久久久影院| 欧美人与性动交α欧美软件| 啦啦啦啦在线视频资源| 国产一区有黄有色的免费视频| 看免费成人av毛片| 丝袜美腿诱惑在线| 成人手机av| 欧美日韩福利视频一区二区| 欧美久久黑人一区二区| 中文字幕av电影在线播放| 亚洲欧美日韩另类电影网站| 午夜福利乱码中文字幕| 欧美精品人与动牲交sv欧美| 精品一区二区免费观看| 久久久久视频综合| 日日摸夜夜添夜夜爱| 亚洲精品第二区| 男女午夜视频在线观看| 亚洲精品久久成人aⅴ小说| 美女高潮到喷水免费观看| 国产成人精品无人区| 日韩精品有码人妻一区| 久久久精品94久久精品| 欧美国产精品va在线观看不卡| 精品国产超薄肉色丝袜足j| 亚洲国产精品国产精品| 久久精品熟女亚洲av麻豆精品| 日本黄色日本黄色录像| 一级,二级,三级黄色视频| 日韩中文字幕视频在线看片| 色视频在线一区二区三区| 无限看片的www在线观看| 建设人人有责人人尽责人人享有的| 亚洲少妇的诱惑av| 人妻一区二区av| 欧美xxⅹ黑人| 亚洲欧美成人综合另类久久久| 国产男人的电影天堂91| 国语对白做爰xxxⅹ性视频网站| 亚洲成人免费av在线播放| 免费女性裸体啪啪无遮挡网站| 亚洲欧洲国产日韩| 天堂8中文在线网| 丝袜脚勾引网站| 极品少妇高潮喷水抽搐| 亚洲国产日韩一区二区| 亚洲一区二区三区欧美精品| 高清av免费在线| 精品少妇久久久久久888优播| av不卡在线播放| 亚洲专区中文字幕在线 | 成年美女黄网站色视频大全免费| 亚洲国产av影院在线观看| 欧美日韩国产mv在线观看视频| 激情视频va一区二区三区| 久久精品熟女亚洲av麻豆精品| av有码第一页| 亚洲美女黄色视频免费看| www.精华液| 日韩欧美精品免费久久| 狂野欧美激情性bbbbbb| 国产日韩欧美亚洲二区| 国产成人一区二区在线| 母亲3免费完整高清在线观看| 成人国产麻豆网| 丰满迷人的少妇在线观看| 美女国产高潮福利片在线看| 亚洲欧美色中文字幕在线| 黄色视频不卡| 一本—道久久a久久精品蜜桃钙片| 一级毛片我不卡| 亚洲人成网站在线观看播放| 免费在线观看黄色视频的| 2021少妇久久久久久久久久久| 亚洲精品一区蜜桃| 日韩熟女老妇一区二区性免费视频| 男女边摸边吃奶| 国产激情久久老熟女| 成年人午夜在线观看视频| 嫩草影视91久久| 国产精品.久久久| 如日韩欧美国产精品一区二区三区| av一本久久久久| 成年人午夜在线观看视频| 涩涩av久久男人的天堂| 亚洲av电影在线进入| 99香蕉大伊视频| 久久人人爽av亚洲精品天堂| 下体分泌物呈黄色| 久久狼人影院| 国产女主播在线喷水免费视频网站| 国产成人91sexporn| 嫩草影院入口| 欧美成人精品欧美一级黄| 十八禁网站网址无遮挡| 欧美亚洲 丝袜 人妻 在线| www.精华液| 久久久久国产精品人妻一区二区| h视频一区二区三区| 久久免费观看电影| 美女扒开内裤让男人捅视频| 超碰成人久久| 国产精品久久久av美女十八| 高清黄色对白视频在线免费看| 亚洲av福利一区| 2021少妇久久久久久久久久久| 无遮挡黄片免费观看| 色94色欧美一区二区| 操出白浆在线播放| 毛片一级片免费看久久久久| 中文字幕制服av| 亚洲中文av在线| 操出白浆在线播放| 亚洲成人手机| 丝瓜视频免费看黄片| 美女高潮到喷水免费观看| 国产一区二区在线观看av| 国产在视频线精品| 免费少妇av软件| 精品免费久久久久久久清纯 | 国产又色又爽无遮挡免| 日韩av在线免费看完整版不卡| svipshipincom国产片| 久久青草综合色| 成年av动漫网址| 精品国产露脸久久av麻豆| 桃花免费在线播放| 免费在线观看视频国产中文字幕亚洲 | 水蜜桃什么品种好| 美女中出高潮动态图| 男的添女的下面高潮视频| 夫妻午夜视频| 欧美成人午夜精品| 一本一本久久a久久精品综合妖精| 国产熟女午夜一区二区三区| 亚洲国产欧美网| 男女午夜视频在线观看| 欧美日韩亚洲综合一区二区三区_| 18禁裸乳无遮挡动漫免费视频| 免费在线观看视频国产中文字幕亚洲 | 十八禁网站网址无遮挡| 中文字幕人妻熟女乱码| 亚洲,一卡二卡三卡| 人人妻,人人澡人人爽秒播 | 国产精品久久久av美女十八| 高清黄色对白视频在线免费看| 午夜久久久在线观看| 精品国产超薄肉色丝袜足j| 成年人免费黄色播放视频| 久久精品久久久久久噜噜老黄| 丰满少妇做爰视频| 天天影视国产精品| 亚洲熟女精品中文字幕| 亚洲成人手机| 午夜福利视频在线观看免费| 久久99一区二区三区| 成人影院久久| 欧美成人精品欧美一级黄| 男女免费视频国产| 一区二区三区四区激情视频| 久久国产精品男人的天堂亚洲| 国产女主播在线喷水免费视频网站| 久久影院123| 热re99久久精品国产66热6| 99re6热这里在线精品视频| 又大又黄又爽视频免费| av卡一久久| 丁香六月天网| 大话2 男鬼变身卡| www.熟女人妻精品国产| 一二三四在线观看免费中文在| 欧美黄色片欧美黄色片| 亚洲国产av新网站| 成年人午夜在线观看视频| 一本一本久久a久久精品综合妖精| bbb黄色大片| 亚洲欧美一区二区三区国产| 亚洲欧美成人精品一区二区| 丝袜人妻中文字幕| 国产99久久九九免费精品| 97在线人人人人妻| 日韩一本色道免费dvd| 少妇猛男粗大的猛烈进出视频| 欧美黑人精品巨大| 日本午夜av视频| 卡戴珊不雅视频在线播放| av线在线观看网站| 国产成人精品久久久久久| 一二三四中文在线观看免费高清| 免费看av在线观看网站| 大香蕉久久网| 久久久精品国产亚洲av高清涩受| 建设人人有责人人尽责人人享有的| 亚洲国产日韩一区二区| 七月丁香在线播放| 黄色视频在线播放观看不卡| 人人妻人人爽人人添夜夜欢视频| 日韩伦理黄色片| 欧美人与性动交α欧美软件| 少妇被粗大的猛进出69影院| 美女主播在线视频| 一本大道久久a久久精品| h视频一区二区三区| 国产黄频视频在线观看| 精品少妇内射三级| 欧美在线黄色| 青草久久国产| 午夜福利视频精品| 国产成人系列免费观看| 亚洲精品日韩在线中文字幕| 黄色怎么调成土黄色| 男女午夜视频在线观看| 国产精品国产三级专区第一集| 中文字幕高清在线视频| 超碰成人久久| 久久精品国产综合久久久| 成年动漫av网址| 免费在线观看黄色视频的| 五月天丁香电影| 中文字幕人妻丝袜一区二区 | 免费日韩欧美在线观看| 亚洲,欧美,日韩| 性少妇av在线| www日本在线高清视频| 19禁男女啪啪无遮挡网站| 免费久久久久久久精品成人欧美视频| 自线自在国产av| av在线观看视频网站免费| 精品国产国语对白av| kizo精华| 日日撸夜夜添| 国产精品国产三级国产专区5o| 免费女性裸体啪啪无遮挡网站| 日本爱情动作片www.在线观看| 久久久国产欧美日韩av| 丝袜人妻中文字幕| 精品国产国语对白av| 午夜免费男女啪啪视频观看| 日韩不卡一区二区三区视频在线| 亚洲五月色婷婷综合| 香蕉丝袜av| 日本av手机在线免费观看| 日韩视频在线欧美| 日韩精品有码人妻一区| 国产精品久久久久久精品电影小说| www.av在线官网国产| 欧美中文综合在线视频| 老熟女久久久| 国产精品一国产av| 一边摸一边抽搐一进一出视频| 免费观看人在逋| 中文字幕高清在线视频| 久久精品熟女亚洲av麻豆精品| av免费观看日本| 一二三四在线观看免费中文在| 美女脱内裤让男人舔精品视频| 中文字幕高清在线视频| 国产在线视频一区二区| 狂野欧美激情性xxxx| 老汉色∧v一级毛片| 成年美女黄网站色视频大全免费| 黑人欧美特级aaaaaa片| 欧美97在线视频| 日韩不卡一区二区三区视频在线| 最近2019中文字幕mv第一页| 色网站视频免费| 9色porny在线观看| avwww免费| 日韩制服骚丝袜av| 国产亚洲一区二区精品| 男女之事视频高清在线观看 | 亚洲成色77777| 久久久久人妻精品一区果冻| 在线观看人妻少妇| 观看av在线不卡| 亚洲欧美一区二区三区国产| 90打野战视频偷拍视频| av天堂久久9| 国产一区有黄有色的免费视频| 欧美人与性动交α欧美精品济南到| 在线 av 中文字幕| 久久久久久久久久久免费av| 精品一区在线观看国产| 少妇精品久久久久久久| 在线观看免费午夜福利视频| 性色av一级| 汤姆久久久久久久影院中文字幕| 欧美日韩av久久| 成年美女黄网站色视频大全免费| 亚洲国产欧美网| 久久精品国产亚洲av涩爱| 欧美中文综合在线视频| 精品亚洲成a人片在线观看| 男人操女人黄网站| 日韩 亚洲 欧美在线| 精品国产乱码久久久久久男人| 久久综合国产亚洲精品| 日韩欧美一区视频在线观看| 日韩一卡2卡3卡4卡2021年| 欧美亚洲日本最大视频资源| 日韩人妻精品一区2区三区| 精品久久久精品久久久| 国产在线视频一区二区| 成人影院久久| 人妻 亚洲 视频| 日韩中文字幕视频在线看片| 成年av动漫网址| 午夜av观看不卡| 国产日韩欧美视频二区| 男的添女的下面高潮视频| 老司机深夜福利视频在线观看 | 2018国产大陆天天弄谢| 亚洲少妇的诱惑av| 少妇人妻 视频| 如何舔出高潮| 国产精品偷伦视频观看了| 777久久人妻少妇嫩草av网站| 欧美日韩亚洲综合一区二区三区_| 高清不卡的av网站| 欧美人与善性xxx| av天堂久久9| 人妻 亚洲 视频| 又粗又硬又长又爽又黄的视频| 欧美精品高潮呻吟av久久| 日韩制服骚丝袜av| 狠狠精品人妻久久久久久综合| 极品少妇高潮喷水抽搐| 最新的欧美精品一区二区| 最近的中文字幕免费完整| 国产福利在线免费观看视频| 精品国产国语对白av| 五月开心婷婷网| 日本av手机在线免费观看| 人妻一区二区av| 久久ye,这里只有精品| 两个人免费观看高清视频| 啦啦啦在线观看免费高清www| 国产精品久久久av美女十八| 成人国产av品久久久| 亚洲 欧美一区二区三区| 欧美黄色片欧美黄色片| 久久影院123| a级片在线免费高清观看视频| 国产成人啪精品午夜网站| 久久狼人影院| svipshipincom国产片| 80岁老熟妇乱子伦牲交| www.熟女人妻精品国产| 午夜福利一区二区在线看| av一本久久久久| 国产一卡二卡三卡精品 | 丝袜脚勾引网站| 日韩av不卡免费在线播放| 一级a爱视频在线免费观看| 麻豆av在线久日| 国产亚洲一区二区精品| 午夜av观看不卡| 黄色 视频免费看| 亚洲一卡2卡3卡4卡5卡精品中文| 纯流量卡能插随身wifi吗| 久久久精品国产亚洲av高清涩受| 97在线人人人人妻| 18禁裸乳无遮挡动漫免费视频| 巨乳人妻的诱惑在线观看| 熟女少妇亚洲综合色aaa.| 考比视频在线观看| 久久热在线av| 亚洲在久久综合| 久久国产亚洲av麻豆专区| 夫妻性生交免费视频一级片| 精品久久久精品久久久| 亚洲人成77777在线视频| 欧美日韩国产mv在线观看视频| 久久人妻熟女aⅴ| 美女国产高潮福利片在线看| 一级,二级,三级黄色视频| 精品国产超薄肉色丝袜足j| 看非洲黑人一级黄片| 精品亚洲乱码少妇综合久久| 桃花免费在线播放| 一区二区日韩欧美中文字幕| 亚洲欧洲精品一区二区精品久久久 | 免费看av在线观看网站| 又大又黄又爽视频免费| 国产亚洲精品第一综合不卡| 久久久久国产一级毛片高清牌| 久久精品aⅴ一区二区三区四区| 黑人猛操日本美女一级片| kizo精华|