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

    基于CFD模擬及PIV技術的魚池流場分析與優(yōu)化

    2022-06-06 07:08:14張成林張宇雷鄭躍平
    漁業(yè)現(xiàn)代化 2022年2期
    關鍵詞:池壁進水口魚池

    張成林,張宇雷,吳 凡,吳 鎧,鄭躍平

    (1 中國水產(chǎn)科學研究院漁業(yè)機械儀器研究所,農(nóng)業(yè)農(nóng)村部漁業(yè)裝備與工程技術重點實驗室,上海 200092;2 青島海洋科學與技術國家實驗室深藍漁業(yè)工程裝備聯(lián)合實驗室,山東青島 266237;3 上海市水生野生動植物保護研究中心,上海 200092)

    魚池作為水產(chǎn)養(yǎng)殖尤其是循環(huán)水養(yǎng)殖系統(tǒng)中最基本也是極為重要的一個生產(chǎn)單元,其形狀、尺寸、結構對養(yǎng)殖產(chǎn)量、養(yǎng)殖效率及后期維護費用有著很大的影響,尤其是魚池的水流特性,直接關系到魚類生長的最佳流態(tài),因此,通過現(xiàn)代科技手段對其流場進行詳細分析測量,探索優(yōu)化流動環(huán)境是循環(huán)水高密度養(yǎng)殖的前提條件。

    魚池設計主要考慮建設成本、空間利用、水質(zhì)維護及便于管理等因素[1],常用形狀有圓形、矩形(方切角)等。圓形魚池水體流動情況較好,方便自行集排污,但土地利用率不高且流速不能準確預測[2-3];矩形魚池易于建設,土地利用率高,但池中會出現(xiàn)死水區(qū),水體交換不均勻,造成局部缺氧或代謝廢物集中于死水區(qū)[4]。由此可見,如何準確地設計魚池的形狀及結構,使其內(nèi)部水流流速分布均勻、水體交換完全、不存在運動死角、集排污效果好,對陸基高效全循環(huán)水養(yǎng)殖工程有重要意義。

    目前,國內(nèi)對于魚池內(nèi)水體流態(tài)的研究仍較少,主要是缺少研究方法和手段。國外方面,研究主要集中在對于魚池流態(tài)流場的描述方面。Joan等[5]使用ADV(Acoustic Doppler Velocimeter)測量了魚池內(nèi)不同深度,不同半徑處的水流切向速度,提出魚池內(nèi)的水體是由自由渦和強制渦組成的組合渦形態(tài)。湛含輝等[6]研究提出利用二次流理論來實現(xiàn)池內(nèi)顆粒污染物的快速收集,但是由于涉及邊界層、渦流現(xiàn)象等問題,在定量研究方面存在較大的難度。此外,在魚池進水方式及流速控制的調(diào)研中,有研究表明,魚池內(nèi)部水體的旋轉速度從池壁到池中心,從池表面到池底都要盡可能地保持一致,它可保證魚池完成自清洗,但速度不應過大,避免能量浪費。水流速度為魚體長度的0.5~2.0倍/s,有利于魚類生長,要求推動可沉淀固體顆粒物到池中心的速度應大于15~30 cm/s[7-9]。由此可見,在魚池的設計中,除結構參數(shù)外,應著重研究進水口位置的布置及進水流速的控制[10-12]。

    魚池形狀及結構的研究可分為傳統(tǒng)經(jīng)驗方法和數(shù)學方法兩類。傳統(tǒng)經(jīng)驗方法主要依靠常年積累的養(yǎng)殖經(jīng)驗,該方法存在一定的主觀性,多限于小范圍的定性研究;數(shù)學方法主要是指利用物理數(shù)學模型進行計算機模擬的方法,習慣上稱為計算流體力學方法(Computational Fluid Dynamics,CFD),是當前的一種非常有效的研究方法,其在水體的流動特性研究方面得到了廣泛的應用[13-15]。粒子圖像測速(Particle Image Velocimetry,PIV)技術作為一種先進的全流場、瞬態(tài)、無接觸測量技術,可將測得的數(shù)據(jù)和控制相間動力學特性的物理機理聯(lián)系起來[16],對部分通過CFD方法得出的模擬結果進行流態(tài)驗證。在CFD方法與PIV技術相結合解決實際問題方面,國內(nèi)外尚無將其應用到魚池相關問題的研究上,但在人工魚礁水動力學及螺旋分離器水流動特性研究方面均有涉及[17-18],其采用三維湍流模擬技術為手段,開展了水動力學的數(shù)值模擬研究,同時與動力水槽試驗,特別是PIV技術的實測結果進行對比分析,使兩類方法得到相互驗證。

    以魚池作為研究對象,采用CFD模擬與PIV技術相結合的手段,研究在不同循環(huán)率及進水流速下的魚池流場特性,同時進一步研究進水口位置變化對魚池流場的影響,以期為今后魚池的結構建造與優(yōu)化提供研究手段及理論依據(jù)。

    1 魚池結構與參數(shù)設計

    1.1 魚池結構

    本研究選用的魚池需符合:滿足魚類快速、健康生長需求;便于實現(xiàn)顆粒污染物的快速收集,保證養(yǎng)殖水質(zhì);滿足操作性和經(jīng)濟性需求。以中國水產(chǎn)科學研究院漁業(yè)機械儀器研究所“農(nóng)業(yè)農(nóng)村部漁業(yè)裝備與工程技術重點實驗室”中的高密度循環(huán)水養(yǎng)殖魚池為研究原型,其主要由3部分組成,分別為進水口、池體和出水口,主要參數(shù)為:池體直徑2.0 m,水深0.8 m,魚池總水體量約為2.5 m3,進水口中心線距離魚池池壁0.1 m,魚池池底坡度為3 %,即魚池中心水深0.83 m,魚池出水口直徑為50 mm。固液混合液體以一定的初速度從進水口流入池體,而后經(jīng)出水口排出。幾何模型及尺寸如圖1所示。

    圖1 魚池幾何模型

    1.2 參數(shù)設計

    選擇對魚池的循環(huán)率、進水口流速和進水口位置這3個關鍵設計參數(shù)進行模擬優(yōu)化。根據(jù)Davidson等國內(nèi)外相關學者對魚池各項參數(shù)的研究與建議[19-21]以及循環(huán)水養(yǎng)殖系統(tǒng)的運行經(jīng)驗,初始運行參數(shù)設計如下:循環(huán)率設置2個梯度,分別為1次/h和2 次/h,對應循環(huán)量(Q)為800 L/h和1 200 L/h;魚池進水速度(Vin)設置3個梯度,分別為1 m/s、2 m/s、3 m/s;進水口位置(H)設置為水面以下10 cm和30 cm。對各組參數(shù)分別進行計算機模擬研究與分析,具體參數(shù)詳見表1。

    表1 魚池結構及運行參數(shù)

    2 數(shù)值模擬與PIV流場測試

    2.1 數(shù)值模擬

    2.1.1 數(shù)學模型

    用ANSYS Fluent 16.1軟件和聯(lián)想服務器進行數(shù)值模擬??紤]到魚池流場特性,首先可將魚池流場視為等溫、不可壓縮三維流場;其次,魚池出水流量由進水流量控制,當魚池流場發(fā)展充分后,在進水流量不變的前提下,魚池中各點的速度、壓力等參數(shù)為常量,即其流場流動狀態(tài)不隨時間變化而變化,可按定常流動求解。故選擇非耦合、隱式定常的3D求解器,采用N-S方程配合標準k-ε雙方程紊流模型對魚池進行三維數(shù)值模擬[22]。運用SIMPLE算法求解壓力-速度耦合,動量、能量、組分、湍動能和湍流耗散率的離散格式均取一階迎風差分格式,壁面附近采用標準壁面函數(shù)方法處理。設定殘差監(jiān)視器,當計算殘差小于10-5時認為計算結果達到收斂。

    標準k-ε模型是由 Launder 等[23]于 1972 年提出的,其雙方程模型是通過補充微分方程的方法來獲得紊流尺度的,用來封閉雷諾方程的微分方程數(shù)如下[13,24]:

    GK-ρε

    (1)

    (2)

    式中:GK為平均速度梯度引起的湍動能k的產(chǎn)生項,定義為:

    (3)

    式中:t為時間,s;xi為橫坐標量,m;xj為縱坐標量,m;k為湍動能,m2/s2;ε為擴散率,m2/s;μ為水的運動黏性系數(shù),m2/s;μt為水的渦黏系數(shù),m2/s;σK和σε分別表示與湍動能k和耗散率ε對應的普朗特常數(shù),C1ε和C2ε均為經(jīng)驗常數(shù)。根據(jù)經(jīng)驗公式模型的相關參數(shù)取值如下:C1ε=1.44,C2ε=1.92,μt=0.09,σK=1.0,σε=1.3。

    2.1.2 網(wǎng)格劃分與邊界條件

    利用GAMBIT建模平臺構建魚池三維模型,并結合魚池模型形狀進行網(wǎng)格劃分及邊界條件設定。在網(wǎng)格劃分中,高質(zhì)量的網(wǎng)格是完成模擬的前提,但同時又要避免因網(wǎng)格過密造成計算量過大。通常,網(wǎng)格質(zhì)量以滿足要求為準,即隨著網(wǎng)格數(shù)目的增加,計算結果不再有顯著的變化即可[25]。根據(jù)上述原則,本試驗模型主要劃分為六面體結構性網(wǎng)格,并將進水口和出水口區(qū)域進行網(wǎng)格加密操作,同時進行了網(wǎng)格獨立性檢驗,最后確定的整體計算域網(wǎng)格數(shù)量為195 483個,具體劃分見圖2;選擇非耦合、隱式定常的3D求解器,運用SIMPLE算法求解壓力-速度耦合,動量、能量、組分、湍動能和湍流耗散率的離散格式均取一階迎風差分格式,壁面附近采用標準壁面函數(shù)方法處理。在邊界條件設定中,進水口選擇速度進口,出水口設定為壓力出口。

    圖2 網(wǎng)格劃分

    2.2 PIV流場測試

    2.2.1 試驗系統(tǒng)

    為準確獲取不同水流工況下魚池的瞬時流場圖像,搭建了流場測試試驗系統(tǒng)。該試驗系統(tǒng)主要包括光學補償魚池、水泵、流量計、控制裝置和PIV系統(tǒng),PIV系統(tǒng)主要部件及參數(shù)為:LaVision激光器,激光能量為2×135 mJ,最大頻率為15 Hz;LaVision VC10-0090 CCD相機,分辨率為 1 936×1 216,最大頻率為20 Hz;高分辨率8通道同步器。

    2.2.2 匹配原則

    PIV粒子圖像測速主要基于基本模式匹配[26-27],其通過激光光源照射所測流場區(qū)域,形成流場光照平面,同時利用高速CCD攝像機同步獲取示蹤粒子的運動軌跡圖像,并記錄相鄰兩張圖像I1和I2的時間間隔t,通過I1和I2中同一示蹤粒子的位移與時間t的比值計算速度場,位移方向即為速度矢量。大量粒子的流動成像即可形成速度矢量場。其基本計算公式如下:

    (4)

    式中:m,n為對應粒子數(shù)。

    在魚池模型匹配方面,考慮到PIV拍攝鏡頭的大小,要對研究對象魚池進行流場實測可行性較差,故采用制作縮尺模型的方法來拍攝實際流場??紤]到魚池流場特性,動力源主要來自進水口,依靠水體黏性力帶動魚池水體整體流動。而且本試驗模型進水口位于液面以下,基本不會在水面產(chǎn)生波浪,故選擇雷諾數(shù)相似能夠保證魚池模型與原型的流動相似性,亦可更加準確地反映小型流場的流動特性。

    2.2.3 試驗方法

    水源為自來水,示蹤粒子選擇聚苯乙烯顆粒,直徑在60~100 μm之間,示蹤粒子的數(shù)目以保證魚池被測區(qū)域內(nèi)粒子分布均勻為宜。試驗系統(tǒng)搭建完成以后,連接好魚池模型的循環(huán)管路,向魚池內(nèi)注入自來水并加入示蹤粒子,然后水泵通電,通過流量計將進水速度調(diào)整到所需數(shù)值,使魚池水體流動一段時間,流態(tài)基本穩(wěn)定后調(diào)節(jié)攝像機鏡頭,使焦距對準中心點,使鏡頭軸線與反應器成90°角,待光線最佳、流態(tài)穩(wěn)定后打開攝像系統(tǒng),通過調(diào)整攝像機的光圈、快門速度、曝光強度、焦距等參數(shù),使拍攝效果達到最佳,在確認拍攝效果可以后將影像保存至硬盤中。對保存的影像進行處理,用計算機軟件提取出圖片中的示蹤粒子運動軌跡并跟基于數(shù)值模擬的魚池流態(tài)進行擬合。

    試驗使用PIV系統(tǒng)二維成像功能,根據(jù)激光照射及光學補償?shù)葘嶋H測量條件,共拍攝魚池內(nèi)4個不同深度水平面上的流速分布情況,以池底為基準分別為:h=140 mm、120 mm、80 mm和10 mm,通過對原始拍攝圖像的分析,得到魚池速度矢量場、顆粒物運動狀態(tài)及魚池渦量場等。

    3 結果與討論

    3.1 模擬結果與分析

    3.1.1 速度分布

    在不同的進水條件下,速度變化規(guī)律類似:水流速度從池壁到魚池中心呈減小趨勢,在中間半徑處達到最小值,而后朝著魚池中心逐漸增加。但在較低的進水流速下,魚池池壁與魚池中心的水流速度變化較小。

    在魚池直徑方向的剖面上,速度分布斜率持續(xù)變化,在靠近池壁處達到最大值,而在魚池中心則為負值。速度斜率在距離魚池中心1/3處有較大提升,尤其是在進水速度較大的試驗條件下,此現(xiàn)象與圓形魚池二次流的相關研究相吻合[6,28]。

    在不同進水條件下,魚池沿直徑方向的速度分布如圖3、圖4所示。

    圖3 不同進水條件下魚池在直徑軸上的速度分布(H=10 cm)

    圖4 不同進水條件下魚池在直徑軸上的速度分布(H=30 cm)

    在相同的進水流量下,進水速度越快,魚池整體流速分布越高,在池壁和魚池中心處尤為明顯,而進水口位置對魚池內(nèi)部流速分布影響差異性不顯著(P>0.05)。

    在相同的進水速度下,進水量越大(循環(huán)率高),魚池整體流速分布越高,且進水量大使得魚池中心流速顯著高于靠近魚池池壁附近的流速,但此現(xiàn)象在進水量小的試驗條件下表現(xiàn)不明顯。同樣,進水口位置對魚池內(nèi)部流速分布影響差異性不顯著(P>0.05)。

    不同的進水口位置,會導致靠近魚池池壁的流速有所不同,較深的進水口位置,其池壁附近的流速相對較快,但對魚池整體流速分布影響甚微。

    魚池左右兩個半徑方向的流速分布并不完全對稱,主要是受到進水口的影響,有進水口分布的一側速度要高,但其差異在實驗允許的誤差范圍之內(nèi),并不影響數(shù)值模擬結果,此結論與國外相關研究相吻合[29]。

    3.1.2 流線分布與顆粒物運動軌跡

    根據(jù)數(shù)值模擬的流場狀態(tài)及流速分布情況,循環(huán)量Q為1 200 L/h時,魚池流場的穩(wěn)定性、均勻性及整體速度分布狀態(tài)較好,以水體循環(huán)量Q為1 200 L/h時各工況下的流態(tài)模擬情況做主要分析。

    當進水口位置在水面以下10 cm時,隨時間發(fā)展的不同進水速度的魚池流場流線圖和固體顆粒物運動軌跡圖(顆粒物直徑設定為1 mm,密度為1.19 g/cm3,從進水口釋放)見圖5、圖6。

    注:左1 m/s,中2 m/s,右3 m/s

    注:左1 m/s,中2 m/s,右3 m/s

    顆粒物運動軌跡與水流發(fā)展并不同步。在循環(huán)率不變且進水口位置固定的前提下,隨著進水口水流流速的增大(進水口直徑變小),魚池中水體旋轉速度變快,水流更容易也更迅速地流向出水口,但同時,由于快速旋轉帶來的離心力增大,顆粒物隨著進水流速的增加,越來越集中的在池底做旋轉運動,無法順利運動到池底中心排出。由此可見,進水速度增大可以加快魚池中水體旋轉速度,但無法促進魚池快速排污;當進水口位置在水面以下30 cm時,除進水口附近及魚池中心處速度稍快外,魚池整體流場流線圖與進水口位置在水下10 cm處并無顯著區(qū)別,但進水口位置遠離水面會使魚池流場變動復雜,顆粒物的運動軌跡更為發(fā)散。

    綜上所述,水體循環(huán)量Q為1 200 L/h時,不同進水條件下,顆粒物運動軌跡和水流運動軌跡并不完全重合;進水流速的增加能夠提高水體混合度,但并不會促進顆粒物的有效排出,反而會因離心力的增大產(chǎn)生相反的結果;在進水口位置在水面以下10 cm、進水流速1 m/s的水利條件下,魚池內(nèi)的水體流態(tài)無論從穩(wěn)定性、排污性能還是節(jié)能方面都更符合實際使用需要。

    3.2 PIV試驗結果與對比分析

    以水體循環(huán)量Q為1 200 L/h、進水口位置在水面以下10 cm、進水流速1 m/s水利條件的縮尺模型為例,進行PIV實測,使用總物質(zhì)平均流速(Vmass average,即所有單位質(zhì)量點的流速的平均值)作為流場描述的特征值。圖7為不同水層位置CFD模擬速度與PIV結果的對比,從圖中可以看出,CFD預測結果與試驗數(shù)據(jù)吻合,誤差基本保持在10%以內(nèi),CFD模擬圖可以較好地反映了魚池的實際流態(tài)。

    以h=100 mm水層為例,CFD模擬結果顯示,進水口附近(左下方)水流流速最快,對應流速約為35.9 cm/s,隨著水流遠離進水口,流速逐漸降低,在靠近進水口背面位置,達到30.5 cm/s左右。該趨勢與PIV實測結果基本一致,分別為38.1 cm/s和34.2 cm/s。另外,模擬結果顯示,沿著半徑方向,水流流速從池壁至半徑中心位置逐步遞減,在半徑中心處降至最低,而后逐漸升高,該現(xiàn)象也與PIV實測結果一致。但在以h=10 mm水層上,CFD模擬的總物質(zhì)平均流速與PIV實測結果誤差相對較大,尤其是進水口和魚池中心處,最大誤差可達21%,究其原因,是該水層靠近魚池池壁及底部邊界層,受邊界層影響較大,在水流速度較大處影響更為顯著。總體考慮,由于研究對象是流體流場,其受到的影響因素眾多,現(xiàn)有的理論基礎、技術條件和手段在這方面尚無法做到精準化模擬,本試驗PIV實測與CFD模擬的總體趨勢一致,產(chǎn)生的誤差亦在可接受的范圍內(nèi),因此,可以認為在流態(tài)描述上,CFD可以準確地模擬出了魚池實際流態(tài)[30]。

    a:h=140 mm,b:h=100 mm,c:h=80 mm,d:h=10 mm

    4 結論

    通過設置不同流量、不同進水速度和不同進水高度,分析了基于CFD模擬的圓形養(yǎng)殖池的徑向速度分布以及基于PIV實測的魚池實際流場。

    (1)決定圓形魚池流場的除魚池幾何參數(shù)外,還有進水流量、進水速度及進水高度,其中進水流量和速度是主要影響因素,這些因素決定了魚池整體流場及流速分布。圓形魚池可以通過擬合進水流量、進水口直徑和進水高度來獲得最佳流速和適宜的流場。

    (2)在圓形魚池直徑2.0 m、水深0.8 m、水體量約為2.4 m3的幾何尺度下,水體循環(huán)量Q為1 200 L/h時(循環(huán)量2 次/h),在進水口位置在水面以下10 cm、進水流速1 m/s的水利條件,魚池內(nèi)的水體流態(tài)無論從穩(wěn)定性、排污性能還是節(jié)能方面都更符合實際使用需要。

    (3)基于CFD技術的流場分析可以準確地反映圓形魚池內(nèi)部流場及速度分布,并能進一步地研究魚池固液兩相流,掌握顆粒物運動規(guī)律。采用CFD技術與PIV技術相結合的手段,可以得出不同魚池形狀及不同長深比、池底坡度等因素下的流場特性,既簡化了試驗手段,又具備實際工程意義,可為今后魚池的結構建造與優(yōu)化提供研究手段及理論依據(jù)。

    猜你喜歡
    池壁進水口魚池
    筏板基礎-敞口水池池壁豎向最大彎矩影響因素分析
    特種結構(2023年6期)2023-12-31 04:35:52
    設計魚池圖
    池壁冷卻風在玻璃纖維窯爐中的應用
    江蘇建材(2023年5期)2023-06-06 20:31:44
    設計魚池圖
    水庫臥管涵管進水口現(xiàn)狀分析及改造型式研究★
    山西建筑(2019年14期)2019-08-17 08:13:44
    兩河口電站進水口穩(wěn)定性分析
    水電站設計(2018年1期)2018-04-12 05:32:09
    魚池這樣挖 養(yǎng)魚能高產(chǎn)
    深降幅水電站進水口分層取水口設計
    鱷魚池
    壁面溫差對鋼筋混凝土圓形水池池壁結構的影響
    亚洲片人在线观看| 国产不卡一卡二| 一a级毛片在线观看| 草草在线视频免费看| 一进一出好大好爽视频| 色综合站精品国产| 91在线精品国自产拍蜜月 | 一进一出抽搐gif免费好疼| 午夜日韩欧美国产| 国产亚洲精品久久久久久毛片| 亚洲电影在线观看av| 欧美日韩精品网址| 俺也久久电影网| 白带黄色成豆腐渣| xxx96com| 欧美黄色淫秽网站| 欧美丝袜亚洲另类 | 此物有八面人人有两片| 日韩欧美国产一区二区入口| 久久亚洲精品不卡| 床上黄色一级片| 午夜亚洲福利在线播放| a级毛片a级免费在线| cao死你这个sao货| 欧美乱妇无乱码| 国产成人aa在线观看| 精品久久蜜臀av无| 国产日本99.免费观看| 国产真实乱freesex| 国产亚洲欧美98| 国产成人欧美在线观看| 国产综合懂色| 欧美乱色亚洲激情| 中文在线观看免费www的网站| 欧美日韩一级在线毛片| 88av欧美| 香蕉久久夜色| 亚洲男人的天堂狠狠| 久99久视频精品免费| 亚洲无线在线观看| 国产精品av视频在线免费观看| 日韩精品青青久久久久久| 91九色精品人成在线观看| 夜夜躁狠狠躁天天躁| 午夜免费激情av| 搡老妇女老女人老熟妇| or卡值多少钱| 日韩欧美精品v在线| 成人国产一区最新在线观看| 男女做爰动态图高潮gif福利片| 黄色成人免费大全| 天天躁狠狠躁夜夜躁狠狠躁| 97碰自拍视频| 精品国产乱子伦一区二区三区| 国产1区2区3区精品| 两性午夜刺激爽爽歪歪视频在线观看| 在线播放国产精品三级| av黄色大香蕉| 亚洲av成人精品一区久久| 精品日产1卡2卡| 人妻久久中文字幕网| 窝窝影院91人妻| 欧美日韩精品网址| 亚洲欧美日韩无卡精品| 久久精品国产清高在天天线| 97超视频在线观看视频| 69av精品久久久久久| 一级作爱视频免费观看| 国产精品久久电影中文字幕| 99热这里只有是精品50| 在线十欧美十亚洲十日本专区| 无限看片的www在线观看| 日本黄色视频三级网站网址| 中文字幕人成人乱码亚洲影| 亚洲成人免费电影在线观看| 白带黄色成豆腐渣| 亚洲av成人一区二区三| 亚洲 欧美一区二区三区| 亚洲第一电影网av| 亚洲精品456在线播放app | 国产精品一及| 黄频高清免费视频| 久久性视频一级片| 免费搜索国产男女视频| 成人特级av手机在线观看| 国产av不卡久久| 亚洲美女视频黄频| 国产精品99久久99久久久不卡| 波多野结衣巨乳人妻| 最新美女视频免费是黄的| 久久人妻av系列| 精品久久蜜臀av无| 国产精品 欧美亚洲| 淫妇啪啪啪对白视频| 白带黄色成豆腐渣| 色哟哟哟哟哟哟| 男人舔女人下体高潮全视频| 一个人免费在线观看电影 | 熟女人妻精品中文字幕| 给我免费播放毛片高清在线观看| 国产高清有码在线观看视频| 欧美激情在线99| 嫩草影院入口| 国产亚洲av嫩草精品影院| 免费在线观看成人毛片| 99热精品在线国产| 欧美日韩综合久久久久久 | av黄色大香蕉| 韩国av一区二区三区四区| 九九在线视频观看精品| 女人高潮潮喷娇喘18禁视频| 一个人免费在线观看的高清视频| 叶爱在线成人免费视频播放| 国产一区二区在线观看日韩 | 久久这里只有精品中国| 欧美av亚洲av综合av国产av| 免费无遮挡裸体视频| 亚洲av五月六月丁香网| 亚洲国产日韩欧美精品在线观看 | 欧美最黄视频在线播放免费| 国产精品电影一区二区三区| 蜜桃久久精品国产亚洲av| 国产精华一区二区三区| 91av网站免费观看| 午夜福利免费观看在线| 国产精品九九99| 久久天堂一区二区三区四区| 在线永久观看黄色视频| 国产真人三级小视频在线观看| 男人舔奶头视频| 丁香欧美五月| 99国产极品粉嫩在线观看| 亚洲精品在线观看二区| 久久久国产欧美日韩av| 女人高潮潮喷娇喘18禁视频| 国产精品香港三级国产av潘金莲| 少妇人妻一区二区三区视频| 国产精品影院久久| av中文乱码字幕在线| 成人国产一区最新在线观看| 婷婷精品国产亚洲av在线| 久久欧美精品欧美久久欧美| 国产三级在线视频| 亚洲乱码一区二区免费版| 97人妻精品一区二区三区麻豆| avwww免费| 一二三四社区在线视频社区8| 亚洲无线在线观看| 日本免费a在线| 噜噜噜噜噜久久久久久91| 亚洲乱码一区二区免费版| 五月伊人婷婷丁香| 亚洲片人在线观看| 成年女人看的毛片在线观看| 午夜视频精品福利| 日韩免费av在线播放| 欧美不卡视频在线免费观看| 亚洲午夜理论影院| 韩国av一区二区三区四区| 国产成+人综合+亚洲专区| 真人做人爱边吃奶动态| 熟妇人妻久久中文字幕3abv| 国产精品 欧美亚洲| 国产人伦9x9x在线观看| 天堂网av新在线| 久久久精品大字幕| 亚洲av中文字字幕乱码综合| 特大巨黑吊av在线直播| 国产高清videossex| 国产av麻豆久久久久久久| 少妇熟女aⅴ在线视频| 日韩精品青青久久久久久| 久久久水蜜桃国产精品网| 日本a在线网址| av女优亚洲男人天堂 | 成人av一区二区三区在线看| 国产精品久久久人人做人人爽| 国产高清有码在线观看视频| e午夜精品久久久久久久| 成人一区二区视频在线观看| 免费av不卡在线播放| 国产亚洲精品一区二区www| 免费在线观看视频国产中文字幕亚洲| 日本熟妇午夜| 婷婷精品国产亚洲av| 国产99白浆流出| 亚洲专区中文字幕在线| 日本与韩国留学比较| 久久精品夜夜夜夜夜久久蜜豆| 亚洲av免费在线观看| 色噜噜av男人的天堂激情| 婷婷六月久久综合丁香| 欧美黄色片欧美黄色片| 久久精品综合一区二区三区| 国产精品久久久久久精品电影| 国内揄拍国产精品人妻在线| 黄片小视频在线播放| 麻豆国产97在线/欧美| 精品国产乱码久久久久久男人| 国产成人精品无人区| 国产一区二区三区在线臀色熟女| 宅男免费午夜| 五月伊人婷婷丁香| 色吧在线观看| 综合色av麻豆| 波多野结衣高清作品| 一a级毛片在线观看| 日韩 欧美 亚洲 中文字幕| 国产精品久久久久久精品电影| 美女高潮喷水抽搐中文字幕| 亚洲美女黄片视频| 色吧在线观看| 嫩草影院精品99| 色播亚洲综合网| 1024手机看黄色片| 91麻豆av在线| 国产成人av激情在线播放| 中文资源天堂在线| 老熟妇乱子伦视频在线观看| 99久久久亚洲精品蜜臀av| 美女大奶头视频| www.999成人在线观看| 制服丝袜大香蕉在线| 久久性视频一级片| 免费搜索国产男女视频| 熟女人妻精品中文字幕| 成人三级做爰电影| 欧美3d第一页| 亚洲国产精品久久男人天堂| 欧美绝顶高潮抽搐喷水| 91字幕亚洲| 色在线成人网| 88av欧美| 免费观看人在逋| 久久久久国内视频| 色视频www国产| 国产av一区在线观看免费| 国产探花在线观看一区二区| 午夜久久久久精精品| 精品国内亚洲2022精品成人| 国产熟女xx| 午夜两性在线视频| 无遮挡黄片免费观看| www.自偷自拍.com| 欧美日本亚洲视频在线播放| 亚洲,欧美精品.| 亚洲18禁久久av| 制服丝袜大香蕉在线| 久久久色成人| 午夜精品久久久久久毛片777| 久久这里只有精品19| 国产v大片淫在线免费观看| 亚洲在线观看片| 熟妇人妻久久中文字幕3abv| 国产99白浆流出| 色噜噜av男人的天堂激情| 淫秽高清视频在线观看| 母亲3免费完整高清在线观看| 亚洲五月天丁香| 热99re8久久精品国产| 午夜成年电影在线免费观看| 亚洲av成人一区二区三| 在线观看免费视频日本深夜| 国产精品98久久久久久宅男小说| 亚洲专区中文字幕在线| 国语自产精品视频在线第100页| 欧美日韩黄片免| 国产成人系列免费观看| 免费人成视频x8x8入口观看| 国产人伦9x9x在线观看| 我的老师免费观看完整版| 国产午夜精品久久久久久| 午夜激情欧美在线| 国产伦在线观看视频一区| 欧美午夜高清在线| 久久性视频一级片| 一级毛片高清免费大全| 人人妻,人人澡人人爽秒播| а√天堂www在线а√下载| 99热这里只有精品一区 | 午夜福利欧美成人| 天天躁日日操中文字幕| 最新美女视频免费是黄的| 三级国产精品欧美在线观看 | 精品久久久久久久人妻蜜臀av| 婷婷精品国产亚洲av在线| 男女之事视频高清在线观看| 在线播放国产精品三级| 日本黄大片高清| 精品久久久久久久毛片微露脸| 两人在一起打扑克的视频| 国内精品一区二区在线观看| av欧美777| 国产精品香港三级国产av潘金莲| 久久香蕉国产精品| 熟女少妇亚洲综合色aaa.| 国产亚洲av高清不卡| 女生性感内裤真人,穿戴方法视频| 美女扒开内裤让男人捅视频| 国产亚洲精品久久久久久毛片| 免费高清视频大片| 亚洲第一欧美日韩一区二区三区| 精品不卡国产一区二区三区| 国产私拍福利视频在线观看| 日韩高清综合在线| 亚洲人与动物交配视频| 最近在线观看免费完整版| 国产亚洲精品av在线| av国产免费在线观看| 午夜福利视频1000在线观看| 在线国产一区二区在线| 午夜a级毛片| 麻豆国产av国片精品| 91在线观看av| 国产蜜桃级精品一区二区三区| 99精品久久久久人妻精品| 欧美绝顶高潮抽搐喷水| 国产精品免费一区二区三区在线| 无人区码免费观看不卡| 白带黄色成豆腐渣| 九色成人免费人妻av| 欧美丝袜亚洲另类 | 国产亚洲av嫩草精品影院| 国产乱人伦免费视频| 在线免费观看的www视频| 国产精华一区二区三区| 国产欧美日韩精品一区二区| 久久精品人妻少妇| 男女之事视频高清在线观看| 色尼玛亚洲综合影院| 久久香蕉精品热| 国产成人影院久久av| 中文字幕精品亚洲无线码一区| 黑人操中国人逼视频| 亚洲自拍偷在线| 麻豆av在线久日| 国产伦精品一区二区三区四那| 午夜福利18| 日本 欧美在线| 国产 一区 欧美 日韩| 动漫黄色视频在线观看| 国产亚洲欧美在线一区二区| 美女高潮的动态| 熟女电影av网| 亚洲av美国av| 精品一区二区三区视频在线 | 免费看日本二区| 婷婷精品国产亚洲av| 国产成人福利小说| 亚洲欧美精品综合久久99| 熟女电影av网| 91av网站免费观看| 精品一区二区三区av网在线观看| 精品久久久久久,| 久久久色成人| 岛国在线观看网站| 天堂√8在线中文| 欧美极品一区二区三区四区| 三级国产精品欧美在线观看 | 手机成人av网站| 久久久久久久午夜电影| 99热这里只有精品一区 | 偷拍熟女少妇极品色| 男女那种视频在线观看| 又粗又爽又猛毛片免费看| 成人特级av手机在线观看| 日韩有码中文字幕| 亚洲欧美精品综合一区二区三区| 天堂动漫精品| 欧美一级毛片孕妇| 一二三四在线观看免费中文在| 日本五十路高清| 成人午夜高清在线视频| 97超级碰碰碰精品色视频在线观看| 免费在线观看视频国产中文字幕亚洲| 亚洲av美国av| 亚洲成a人片在线一区二区| 国产精品98久久久久久宅男小说| 国产三级黄色录像| 成年女人毛片免费观看观看9| 亚洲 欧美 日韩 在线 免费| 亚洲av第一区精品v没综合| 成人一区二区视频在线观看| 国产主播在线观看一区二区| 精品一区二区三区四区五区乱码| 一二三四在线观看免费中文在| 亚洲性夜色夜夜综合| 亚洲人成网站高清观看| 午夜福利欧美成人| 非洲黑人性xxxx精品又粗又长| 国产精品亚洲一级av第二区| 午夜亚洲福利在线播放| 窝窝影院91人妻| 国产成人精品久久二区二区免费| 男人舔奶头视频| 国产精品久久久久久人妻精品电影| 人妻丰满熟妇av一区二区三区| 村上凉子中文字幕在线| 俺也久久电影网| 国内精品一区二区在线观看| 国产成年人精品一区二区| 两性夫妻黄色片| 亚洲成av人片在线播放无| 在线免费观看的www视频| 美女高潮的动态| 一级黄色大片毛片| 黄色视频,在线免费观看| 国产又黄又爽又无遮挡在线| 人人妻人人澡欧美一区二区| 免费看十八禁软件| 特级一级黄色大片| 久久精品影院6| 亚洲avbb在线观看| 成人av一区二区三区在线看| 中国美女看黄片| 国产精品亚洲美女久久久| 国产精品一及| 又爽又黄无遮挡网站| 国产99白浆流出| 亚洲国产精品成人综合色| 美女被艹到高潮喷水动态| 色综合亚洲欧美另类图片| 日本黄色片子视频| 夜夜看夜夜爽夜夜摸| 成人午夜高清在线视频| 丁香六月欧美| 18禁国产床啪视频网站| x7x7x7水蜜桃| 亚洲熟女毛片儿| 一个人免费在线观看电影 | 美女黄网站色视频| 校园春色视频在线观看| 99国产极品粉嫩在线观看| 成年人黄色毛片网站| 日本黄色视频三级网站网址| 18禁黄网站禁片午夜丰满| 国产野战对白在线观看| 午夜成年电影在线免费观看| 男女床上黄色一级片免费看| 欧美高清成人免费视频www| 国产伦人伦偷精品视频| 国产精品99久久99久久久不卡| 桃红色精品国产亚洲av| 久久久久久人人人人人| 女人被狂操c到高潮| 亚洲片人在线观看| 高潮久久久久久久久久久不卡| 亚洲熟妇熟女久久| 韩国av一区二区三区四区| 亚洲av熟女| 国产精品1区2区在线观看.| 国产亚洲av嫩草精品影院| 999久久久国产精品视频| 国产日本99.免费观看| 99re在线观看精品视频| 国产激情久久老熟女| 亚洲美女黄片视频| 99久久精品国产亚洲精品| 久久香蕉精品热| 日日干狠狠操夜夜爽| 中文字幕高清在线视频| 99精品欧美一区二区三区四区| 国产精品永久免费网站| 午夜福利欧美成人| 女人被狂操c到高潮| 国产av在哪里看| 国产麻豆成人av免费视频| 国产男靠女视频免费网站| 一卡2卡三卡四卡精品乱码亚洲| 精品一区二区三区av网在线观看| www日本黄色视频网| 狂野欧美激情性xxxx| 老熟妇仑乱视频hdxx| 51午夜福利影视在线观看| 男插女下体视频免费在线播放| 欧美乱色亚洲激情| 久久香蕉国产精品| 国产精品野战在线观看| 精品久久久久久久毛片微露脸| 久久午夜综合久久蜜桃| 久久亚洲精品不卡| 亚洲精品456在线播放app | 99re在线观看精品视频| 久久久久国产一级毛片高清牌| 在线视频色国产色| 99久国产av精品| 久久久国产欧美日韩av| 18禁黄网站禁片午夜丰满| 老熟妇仑乱视频hdxx| 久久久色成人| 久久精品91蜜桃| 色精品久久人妻99蜜桃| 不卡av一区二区三区| 视频区欧美日本亚洲| 99热精品在线国产| 十八禁网站免费在线| 九色国产91popny在线| 草草在线视频免费看| 久久久久精品国产欧美久久久| 国产午夜精品久久久久久| 不卡av一区二区三区| 免费观看的影片在线观看| 欧美又色又爽又黄视频| 国产激情久久老熟女| 天堂√8在线中文| 亚洲自偷自拍图片 自拍| 高清毛片免费观看视频网站| 国产av在哪里看| 亚洲第一电影网av| 久久久久久大精品| 亚洲色图av天堂| 好男人在线观看高清免费视频| 久久久国产欧美日韩av| 国产高清三级在线| 淫秽高清视频在线观看| 白带黄色成豆腐渣| 国产美女午夜福利| 热99在线观看视频| 一级毛片高清免费大全| 嫩草影院入口| 此物有八面人人有两片| www.熟女人妻精品国产| 噜噜噜噜噜久久久久久91| 白带黄色成豆腐渣| 女警被强在线播放| 黄色女人牲交| 亚洲午夜理论影院| 99riav亚洲国产免费| 婷婷精品国产亚洲av在线| 久久久久精品国产欧美久久久| 九九热线精品视视频播放| 亚洲精品一卡2卡三卡4卡5卡| 欧美日韩一级在线毛片| 亚洲片人在线观看| 一区二区三区国产精品乱码| 69av精品久久久久久| 久久久久久九九精品二区国产| 99久久精品国产亚洲精品| 国产私拍福利视频在线观看| 亚洲无线观看免费| 毛片女人毛片| 亚洲成人精品中文字幕电影| 精品一区二区三区av网在线观看| 51午夜福利影视在线观看| 999精品在线视频| 色综合欧美亚洲国产小说| 变态另类丝袜制服| 精品福利观看| 一进一出抽搐动态| 日韩欧美一区二区三区在线观看| 看黄色毛片网站| 嫩草影院入口| 人人妻,人人澡人人爽秒播| 午夜久久久久精精品| 久久中文字幕一级| 欧美三级亚洲精品| 婷婷丁香在线五月| 亚洲中文av在线| av天堂中文字幕网| 99热6这里只有精品| 国产精品久久久人人做人人爽| 91在线观看av| 给我免费播放毛片高清在线观看| 亚洲成人中文字幕在线播放| 窝窝影院91人妻| 亚洲五月天丁香| 又粗又爽又猛毛片免费看| 黄色丝袜av网址大全| 综合色av麻豆| 全区人妻精品视频| 女人被狂操c到高潮| 久久精品影院6| 精品福利观看| 桃色一区二区三区在线观看| 五月玫瑰六月丁香| 亚洲欧美日韩东京热| 99在线视频只有这里精品首页| 啦啦啦免费观看视频1| 精品国内亚洲2022精品成人| 老鸭窝网址在线观看| 精品人妻1区二区| 国产成人av激情在线播放| 欧美zozozo另类| 97碰自拍视频| 黑人欧美特级aaaaaa片| 成人午夜高清在线视频| 日韩欧美国产一区二区入口| 两个人看的免费小视频| 男女之事视频高清在线观看| 18禁观看日本| 国产又黄又爽又无遮挡在线| 狠狠狠狠99中文字幕| 一进一出抽搐gif免费好疼| 亚洲第一电影网av| 制服丝袜大香蕉在线| 午夜视频精品福利| av黄色大香蕉| 精品国产乱码久久久久久男人| 国产一区二区三区视频了| 国产成人啪精品午夜网站| 久久中文字幕人妻熟女| 不卡一级毛片| 日本三级黄在线观看| 黄色视频,在线免费观看| 黄色 视频免费看| 欧美中文综合在线视频| 中文字幕熟女人妻在线| 亚洲国产高清在线一区二区三| 人妻久久中文字幕网| 噜噜噜噜噜久久久久久91| 不卡一级毛片| 国内久久婷婷六月综合欲色啪| 欧美日韩乱码在线| 好男人电影高清在线观看| 久久久久久久久免费视频了|