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

    水力旋流器旋流運動的格子Boltzmann模擬

    2018-01-15 09:24:47胡亮朱衛(wèi)兵張小彬
    哈爾濱工程大學(xué)學(xué)報 2017年12期
    關(guān)鍵詞:切向速度旋流器旋流

    胡亮, 朱衛(wèi)兵, 張小彬

    (哈爾濱工程大學(xué) 航天與建筑工程學(xué)院,黑龍江 哈爾濱 150001)

    利用離心力將兩相分離在工程中已被廣泛采用,如旋風(fēng)除塵器、石油工業(yè)氣液旋流器和航空發(fā)動機動壓式油氣分離器等[1-2]。根據(jù)氣液主相不同,通常稱為旋風(fēng)分離器和水力旋流器。

    Erdal、Hreiz等[3-5]通過實驗探討了結(jié)構(gòu)參數(shù)對水力旋流器流場的影響,結(jié)果表明入口結(jié)構(gòu)及數(shù)量對流場分布均有較大影響。Junwei Yang等[6]建立了水力旋流器液相停留時間預(yù)測模型。研究表明,水力旋流器內(nèi)部流動復(fù)雜,存在三維強旋湍流運動,且軸向存在二次流擾動。在數(shù)值模擬方面,Hreiz等[7]利用Fluent中幾種湍流模型系統(tǒng)模擬了水力旋流器內(nèi)的旋流場,發(fā)現(xiàn)Realizableκ-ε模型和LES模型對于流場細(xì)節(jié)的捕捉還不夠準(zhǔn)確,需要更有效的數(shù)值方法。格子Boltzmann方法(lattice Boltzmann method,LBM)能夠真實的反映流體粒子間的相互作用,并提高流場的計算精度[8],具有描述流場局部二次流現(xiàn)象的潛力,而這類問題的研究可以加深對水力旋流器的理論認(rèn)識,為水力旋流器的研制提供理論依據(jù);針對三維流動問題,LBM演化方程為顯式方程,無需迭代計算,在非穩(wěn)態(tài)計算中效率很高;同時LBM具有易并行的優(yōu)點。Derksen等[9-11]采用LBM對旋風(fēng)分離器氣相旋流場進(jìn)行研究,準(zhǔn)確描述了旋渦中心的旋進(jìn)渦核(precessing vortex core,PVC)現(xiàn)象,即旋進(jìn)渦核現(xiàn)象。上述研究表明,LBM已成功應(yīng)用在旋風(fēng)分離器的流場分析中。與旋風(fēng)分離器相比,由于主相和壁面影響的不同,水力旋流器的流場存在更大的非軸對稱性[12],且流場中心的二次回流區(qū)域也有明顯差異[7]。

    Bhanjee等[13]探討了LBM在水力旋流器模擬中的應(yīng)用方法,結(jié)果表明LBM給出的旋流場流動趨勢與以往認(rèn)知一致,但未對計算結(jié)果做出進(jìn)一步定量分析。本文以文獻(xiàn)[12]的實驗數(shù)據(jù)為基礎(chǔ),進(jìn)一步探討LBM在水力旋流器中的應(yīng)用,通過流場描述、定量分析及對旋進(jìn)渦核現(xiàn)象模擬等,深入了解LBM在此問題上的應(yīng)用前景。

    1 Boltzmann模型

    1.1 格子Boltzmann方法

    本文采用D3Q19單松弛格子Boltzmann模型進(jìn)行數(shù)值模擬,其演化方程即含外力項的格子Boltzmann-BGK方程如下

    fi(r+eiδt,t+δt)-fi(r,t)=δtFi(r,t)-

    (1)

    (2)

    式中:c=δx/δt,δx、δt分別為網(wǎng)格步長和時間步長。

    宏觀密度和速度分別為

    (3)

    (4)

    式中:cs為格子聲速,ωi為權(quán)函數(shù):ωi=1/3(i為原點位置),ωi=1/18(i為6個坐標(biāo)軸方向),ωi=1/36(i為12個對角線方向)。

    1.2 湍流模型

    palabos中湍流模型采用LES中的標(biāo)準(zhǔn)Smagorinsky模型,式(1)中的無量綱弛豫時間τ計算公式為

    (5)

    (6)

    式中:cs=0.14,S為應(yīng)變率張量。

    2 物理模型及網(wǎng)格無關(guān)性驗證

    2.1 物理模型

    采用文獻(xiàn)[12]的實驗?zāi)P?,模型結(jié)構(gòu)與參數(shù)如圖1所示。本文獲取的平均速度分布截面位置均在圖中實驗測試平面上,如圖1所示。

    根據(jù)文獻(xiàn)[12]提供的數(shù)據(jù),模擬選擇了四種不同實驗條件,參數(shù)如表1所示。表中,ν為流體運動粘度,Q為流體體積流量,Uav為分離器筒體平均軸向速度。本文采用速度入口邊界條件,出口為自由出流,壁面邊界采用具有二階精度的Guo格式非平衡外推邊界。

    圖1 物理模型Fig.1 Geometry model

    工況ν×106/(m2·s-1)Q/(m3·s-1)Uav/(m·s-1)Re170.000680.1091514270.003390.5457570310.000630.1029285410.004540.73166855

    2.2 網(wǎng)格無關(guān)性驗證

    由于采用palabos中的分塊網(wǎng)格技術(shù),計算網(wǎng)格量統(tǒng)計只涉及流體域相關(guān)的網(wǎng)格。本文選取三種網(wǎng)格模型對工況3進(jìn)行計算,并對截面X=31.7 cm處的平均切向速度分布進(jìn)行對比,對比結(jié)果如圖2所示。其中速度值進(jìn)行了無量綱處理,下文處理方式一致,不再贅述。

    圖2 截面X=31.7 cm處平均切向速度分布Fig.2 Average tangential velocity profiles at X=31.7 cm

    經(jīng)過對比,1 430萬和2 000萬網(wǎng)格量的模型計算結(jié)果比較接近,最大相差8.4%。考慮到計算時長的影響,本文研究選擇1 430萬網(wǎng)格量的計算模型,其中計算網(wǎng)格步長δx=1.2 mm,時間步長δt=6.23×10-5s。計算在Intel core i7 4770處理器上進(jìn)行,CPU主頻為3.40 GHz、內(nèi)存16 G、8核并行,每個工況計算時長約100 h。

    3 水力旋流器疏場研究與分析

    水力旋流器存在復(fù)雜的三維強旋湍流運動,其主要體現(xiàn)為各分向速度分布以及旋渦形狀變化;其中切向速度對流體粒子所受離心力有直接影響,軸向速度與流體粒子的停留時間呈正比,徑向速度引發(fā)流場的偏擺現(xiàn)象,旋渦形狀則關(guān)系到氣泡的捕獲位置。

    3.1 流場特征分析

    選取工況3進(jìn)行流場特征分析,圖3給出了水力旋流器內(nèi)瞬時速度分布云圖,物理時間取15.22 s。

    由圖3(a)可見,軸向瞬時速度沿軸向劇烈變化,流場分布不連續(xù),且具有非對稱性,其中軸向零速包絡(luò)面為黑色與灰色區(qū)域交界面,該面呈蛇形分布;由圖3(b)可看出,徑向速度在軸向伴有

    正負(fù)交替,說明流動是不穩(wěn)定的,這是造成工程中氣柱左右搖擺的原因;由圖3(c)可見,渦核中心(即切向速度零值點)并非沿中心軸線對稱,而是在軸線附近擺動。流場分析表明,LBM能夠模擬水力旋流器內(nèi)的旋流運動,且捕捉到了二次回流現(xiàn)象和渦核的非穩(wěn)態(tài)特征,流場符合水力旋流器典型的分布特點。

    圖3 工況3瞬時速度分布云圖Fig.3 contours of instantaneous velocity in case 3

    3.2 平均切向速度分布對比

    圖4~7為工況1~4時不同截面的平均切向速度分布。由圖4可見,本文獲得的平均切向速度分布符合蘭金渦(Rankine vortex)運動規(guī)律,即自由渦和強制渦組成的復(fù)合渦運動。中心區(qū)域平均切向速度呈線性變化,兩端靠近壁面區(qū)域平均切向速度逐漸減?。挥捎诒诿孀饔靡约傲鲃雍纳⒌挠绊?,在壁面附近出現(xiàn)較大的速度梯度,且隨截面位置下移,切向速度變小。結(jié)果對比表明,本文模擬結(jié)果與文獻(xiàn)[7]的實驗結(jié)果及模擬結(jié)果基本相符。圖4中還可以看出,旋轉(zhuǎn)中心位置(即切向速度零點)不在軸線上,且不同截面旋轉(zhuǎn)中心的位置發(fā)生偏移,這一現(xiàn)象與文獻(xiàn)[3]結(jié)論吻合。

    通過對比圖5和圖6可以看出,當(dāng)Re數(shù)相近而粘度變化時,模擬結(jié)果與實驗結(jié)果吻合度相差不大,說明在小范圍變化時,粘度對本文采用方法的準(zhǔn)確性影響較小。

    圖4 工況1平均切向速度分布Fig.4 Average tangential velocity profiles in case 1

    圖5 工況2平均切向速度分布Fig.5 Average tangential velocity profiles in case 2

    圖6 工況3平均切向速度分布Fig.6 Average tangential velocity profiles in case 3

    圖7 工況4平均切向速度分布Fig.7 Average tangential velocity profiles in case 4

    在水力旋流器中,平均切向速度最大值在一定程度上反映了流場的旋流特點,通過式(7)可得到平均切向速度最大值處實驗值與模擬值的相對誤差:

    (7)

    式中:WEX為最大切向速度實驗值,WLBM為同一位置的模擬值。表2給出了不同工況切向速度最大值的相對誤差,其中a、b分別表示本文模擬結(jié)果和文獻(xiàn)[7]Realizableκ-ε模型的模擬結(jié)果。

    由表2可見,當(dāng)Re較小時,本文模擬結(jié)果與實驗值的相對誤差均大于文獻(xiàn)[7]模擬結(jié)果的相對誤差;隨著Re增大,兩種模擬結(jié)果的相對誤差值逐漸接近,在部分截面,本文結(jié)果與實驗值的吻合度高于文獻(xiàn)[7]結(jié)果;同時,隨著Re增大,本文模擬結(jié)果與實驗結(jié)果偏差減小,說明LBM更適用于模擬較大

    Re的旋流場。

    表2切向速度最大值相對誤差

    Table2Relativeerrorofmaximumtangentialvelocity

    %

    由圖4~7對比可知:本文模擬得到的各截面平均切向速度分布與實驗結(jié)果趨勢一致,說明LBM能夠定性的模擬水力旋流器中存在的旋流運動;與文獻(xiàn)[7]相比,本文模擬結(jié)果與中心強制旋流區(qū)域的實驗值更接近,體現(xiàn)了LBM作為介觀方法能夠真實反映流體粒子間相互作用的特點。與實驗結(jié)果相比,本文對較小Re旋流場的模擬結(jié)果存在偏差,當(dāng)Re增大后誤差減小,說明本文采用方法對較小Re的流場模擬存在局限性;當(dāng)粘度在小范圍變化時,對本文模擬結(jié)果準(zhǔn)確性影響很小。

    對于靠近壁面的平均切向速度,本文模擬結(jié)果與實驗結(jié)果存在一定差異,且大部分速度值小于實驗值,分析可能因為本文采用的湍流模型為標(biāo)準(zhǔn)Smagorinsky模型,該模型具有各向同性的特點,且在實際應(yīng)用中存在耗散過大的問題[14],因而在邊界處出現(xiàn)模擬結(jié)果與實驗結(jié)果差異明顯的現(xiàn)象;而文獻(xiàn)[7]采用的LES模型為WALE模型(即壁面適應(yīng)局部渦粘模型),WALE模型考慮了旋轉(zhuǎn)張量的影響,在靠近壁面處能獲得較好的渦粘性,從而改善了標(biāo)準(zhǔn)Smagorinsky模型湍流耗散量大的問題[15],所以文獻(xiàn)[7]的LES模擬結(jié)果與實驗值更接近。

    3.3 平均軸向速度分布對比

    圖8~11為工況1~4不同截面平均軸向速度分布。由圖8可見,本文獲得的平均軸向速度分布體現(xiàn)了水力旋流器中沿器壁向下的外旋流運動以及中心軸線附近向上的內(nèi)旋流運動。結(jié)果對比表明,本文模擬結(jié)果與文獻(xiàn)[7]的實驗結(jié)果及模擬結(jié)果趨勢相同。對于靠近壁面的平均軸向速度,與前述切向速度對比結(jié)果類似,本文模擬結(jié)果與實驗結(jié)果存在一定差異。

    圖8 工況1平均軸向速度分布Fig.8 Average axial velocity profiles in case 1

    圖9 工況2平均軸向速度分布Fig.9 Average axial velocity profiles in case 2

    對比圖9和圖10可看出粘度對平均軸向速度模擬效果的影響。通過對比可見,粘度變化時模擬結(jié)果與實驗結(jié)果吻合度相差不大,說明在小范圍變化時,粘度對本文采用方法的準(zhǔn)確性影響較小,與上文結(jié)論一致。

    與切向速度誤差處理類似,表3給出了不同工況軸向速度最大值的相對誤差。

    由表3可以看出,本文平均軸向速度的模擬結(jié)果與實驗結(jié)果存在一定偏差,隨著Re增大,誤差減?。划?dāng)Re為66 855時,誤差值有所增大,但仍小于工況1的誤差;變化趨勢與上文切向速度誤差對比基本相同。

    圖10 工況3平均軸向速度分布Fig.10 Average axial velocity profiles in case 3

    圖11 工況4平均軸向速度分布Fig.11 Average axial velocity profiles in case 4

    X/mm工況1工況2工況3工況4abababab31.744.97.41.96.811.720.733.010.112.546.73.110.81.58.759.423.212.168.616.74.375.725.95.889.915.115.85.025.5

    由圖8~11對比可知:本文模擬得到的各截面平均軸向速度分布與實驗結(jié)果基本吻合,說明LBM能夠模擬水力旋流器中存在的二次回流運動;與文獻(xiàn)[7]的數(shù)值模擬結(jié)果相比,本文在水力旋流器中心二次回流區(qū)域的模擬結(jié)果比實驗結(jié)果更接近,體現(xiàn)了LBM作為介觀方法能夠真實反映流體粒子間相互作用的特點。與實驗結(jié)果相比,本文模擬得到的平均軸向速度分布波動偏大,分析是因為LBM從介觀角度出發(fā),更加強調(diào)對流體粒子間相互作用的描述,模擬得到的湍流脈動強度大于實驗結(jié)果,從而造成平均軸向速度分布的波動大于實驗結(jié)果,這與文獻(xiàn)[9]的結(jié)論一致。

    3.4 平均徑向速度分布

    在水力旋流器內(nèi)徑向速度測量比較困難,文獻(xiàn)[12]也未給出實驗結(jié)果。針對工況3的平均徑向速度分布,圖12給出了本文以及文獻(xiàn)[7]模擬的結(jié)果。

    圖12(a)中平均徑向速度的絕對值隨半徑減小而增大,這與文獻(xiàn)[16]對平均徑向速度分布的描述吻合;對比圖6和圖10可見,徑向速度與軸向速度以及切向速度量級相等,這與文獻(xiàn)[7]結(jié)論相同;同時,由圖12(b)可見,沿軸向平均徑向速度方向是正負(fù)交替變化的,這與文獻(xiàn)[7]結(jié)果(圖12(c))一致。

    圖12 工況3平均徑向速度分布Fig.12 Average radial velocity profiles in case 3

    圖13為工況3平均軸向速度分布云圖,圖13中箭頭指示方向為根據(jù)圖12標(biāo)記的平均徑向速度方向,方向的左右交替造成了渦核的偏擺;平均速度與瞬時速度對流場的影響規(guī)律一致,但與圖3(a)對比可見流場空間形態(tài)是不同的,這也體現(xiàn)了渦核位置的瞬時性。

    3.5 旋進(jìn)渦核現(xiàn)象

    旋進(jìn)渦核普遍存在于旋轉(zhuǎn)流動中,是一種非穩(wěn)態(tài)流動現(xiàn)象,它的存在直接影響到水力旋流器的分離過程。為了更清楚的了解渦核的非穩(wěn)態(tài)變化,針對工況3,在模擬過程中每隔0.3 s記錄X=89.9 cm截面處渦核中心的位置,渦核中心的運動軌跡如圖14所示。

    圖13 平均徑向速度對旋渦的影響Fig.13 Influence of the average radial velocity to vortex

    圖14 水力旋流器PVC現(xiàn)象Fig.14 PVC phenomenon in hydrocyclone

    從圖14可以看出,流場存在明顯的外旋流和內(nèi)旋流,其中外旋流方向向下,內(nèi)旋流方向向上;同時,流場并不是繞著幾何中心旋轉(zhuǎn),呈現(xiàn)明顯的非穩(wěn)態(tài)特征。隨著時間變化,內(nèi)外流動區(qū)域相對固定,但渦核中心的位置持續(xù)變化,且旋轉(zhuǎn)方向與流場旋轉(zhuǎn)方向相反,這與文獻(xiàn)[7]結(jié)論一致。從上述分析可知,LBM能夠定性的模擬水力旋流器存在的旋進(jìn)渦核現(xiàn)象,且渦核的運動趨勢與文獻(xiàn)[7]吻合。

    4 結(jié)論

    1)LBM能夠模擬出水力旋流器中存在的旋流運動,可捕捉到二次回流現(xiàn)象和渦核的非穩(wěn)態(tài)特征,模擬結(jié)果符合水力旋流器典型的分布特點;

    2)LBM模擬獲得的平均軸向速度和平均切向速度分布與實驗結(jié)果基本吻合,對流場中心二次流區(qū)域的模擬結(jié)果比文獻(xiàn)[7]更接近實驗值,體現(xiàn)了LBM對復(fù)雜流動模擬的優(yōu)勢;

    3)模擬可得出平均徑向速度沿軸向的正負(fù)交替變化及流場的旋進(jìn)渦核現(xiàn)象。

    [1] 林基恕, 常春江, 楊九高. 航空發(fā)動機設(shè)計手冊[M]. 北京: 航空工業(yè)出版社, 2002.

    LIN Jishu, CHANG Chunjiang, YANG Jiugao. Aero engine design manual[M]. Beijing: Navigate Industrical Press, 2002.

    [2] GRUSELLE F, STEIMES J, HENDRICK P. Study of a two-phase flow pump and separator system[J]. Journal of engineering for gas turbines and power, 2011, 133(6): 062401.

    [3] ERDAL F M, SHIRAZI S A. Local velocity measurements and computational fluid dynamics (CFD) simulations of swirling flow in a cylindrical cyclone separator[J]. Journal of energy resources technology, 2004, 126(4): 326-333.

    [4] HREIZ R, LAIN R, WU J, et al. On the effect of the nozzle design on the performances of gas-liquid cylindrical cyclone separators[J]. International journal of multiphase flow, 2014, 58: 15-26.

    [5] HREIZ R, GENTRIC C, MIDOUX N, et al. Hydrodynamics and velocity measurements in gas-liquid swirling flows in cylindrical cyclones[J]. Chemical engineering research and design, 2014, 92(11): 2231-2246.

    [6] YANG J, ZHANG X, SHEN G, et al. Modeling the mean residence time of liquid phase in the gas-liquid cyclone[J]. Industrial & engineering chemistry research, 2015, 54(43): 10885-10892.

    [7] HREIZ R, GENTRIC C, MIDOUX N. Numerical investigation of swirling flow in cylindrical cyclones[J]. Chemical engineering research and design, 2011, 89(12): 2521-2539.

    [8] YU H, GIRIMAJI S S, LUO L-S. DNS and LES of decaying isotropic turbulence with and without frame rotation using lattice Boltzmann method[J]. Journal of computational physics, 2005, 209(2): 599-616.

    [9] DERKSEN J, VAN DEN AKKER H. Simulation of vortex core precession in a reverse-flow cyclone[J]. AIChE journal, 2000, 46(7): 1317-1331.

    [10] LI Y, SHOCK R, ZHANG R, et al. Simulation of turbulent flow in a cyclonic separator with Lattice Boltzmann method[C]//36th AIAA Fluid Dynamics Conference and Exhibit. San Francisco, USA, 2006.

    [11] PIRKER S, GONIVA C, KLOSS C, et al. Application of a hybrid Lattice Boltzmann-finite volume turbulence model to cyclone short-cut flow[J]. Powder technology, 2013, 235: 572-580.

    [12] ERDAL F M. Local measurements and computational fluid dynamics simulations in a gas-liquid cylindrical cyclone separator[D]. Tulsa: The University of Tulsa, 2001.

    [13] BHAMJEE M, CONNELL S H, NEL A L. An investigation into the applicability of the Lattice Boltzmann method to modelling of the flow in a hydrocyclone[C]//South African Conference on Computational and Applied Mechanics. [S.l.],2014.

    [14] BEAUBERT F, VIAZZO S. Large eddy simulations of plane turbulent impinging jets at moderate Reynolds numbers[J]. International journal of heat and fluid flow, 2003, 24(4): 512-519.

    [15] NICOUD F, DUCROS F. Subgrid-scale stress modelling based on the square of the velocity gradient tensor[J]. Flow, turbulence and combustion, 1999, 62(3): 183-200.

    [16] 龐學(xué)詩. 水力旋流器技術(shù)與應(yīng)用[M]. 北京:中國石化出版社, 2010.

    PANG Xueshi. Technology and application of hydrocyclone[M]. Beijing: China Petrochemical Press, 2010.

    本文引用格式:

    胡亮, 朱衛(wèi)兵, 張小彬. 水力旋流器旋流運動的格子Boltzmann模擬[J]. 哈爾濱工程大學(xué)學(xué)報, 2017, 38(12): 1864-1871.

    HU Liang, ZHU Weibing, ZHANG Xiaobin. Lattice Boltzmann simulation of swirling flow in hydrocyclone[J]. Journal of Harbin Engineering University, 2017, 38(12): 1864-1871.

    猜你喜歡
    切向速度旋流器旋流
    重介質(zhì)旋流器選煤技術(shù)在我國的創(chuàng)新發(fā)展與應(yīng)用
    旋風(fēng)分離器內(nèi)氣相旋轉(zhuǎn)流不穩(wěn)定性的實驗研究
    雙級徑向旋流器對燃燒性能的影響
    燒嘴旋流器優(yōu)化設(shè)計計算
    雙旋流耦合式旋流反應(yīng)器內(nèi)切向速度分布研究
    流體機械(2020年4期)2020-05-12 09:20:48
    上旋流厭氧反應(yīng)器在造紙廢水處理中的應(yīng)用
    一級旋流偏置對雙旋流杯下游流場的影響
    微氣泡旋流氣浮選設(shè)備的研究與應(yīng)用
    旋流脫氣性能影響因素的CFD模擬
    排水采氣旋流器的固定裝置設(shè)計
    機械工程師(2015年9期)2015-02-26 08:38:11
    天美传媒精品一区二区| kizo精华| 精品少妇久久久久久888优播| 秋霞伦理黄片| 99久国产av精品国产电影| 色吧在线观看| 22中文网久久字幕| 国产精品久久久久成人av| 男女免费视频国产| 中文资源天堂在线| 亚洲第一av免费看| 高清毛片免费看| 日本黄色片子视频| 亚洲欧美日韩另类电影网站| 一本久久精品| 久久综合国产亚洲精品| 伦理电影免费视频| 国产精品女同一区二区软件| 国产精品一区www在线观看| 婷婷色av中文字幕| 亚洲精品一二三| 亚洲成人av在线免费| 在线观看免费日韩欧美大片 | 亚洲精品,欧美精品| 久久久久精品久久久久真实原创| 精品久久国产蜜桃| 精品午夜福利在线看| freevideosex欧美| 大香蕉久久网| 成人综合一区亚洲| 精品亚洲乱码少妇综合久久| 亚洲精品456在线播放app| 在线免费观看不下载黄p国产| 久久久久网色| 国产亚洲精品久久久com| 好男人视频免费观看在线| 人人妻人人添人人爽欧美一区卜| 国产精品一区二区在线不卡| 国产精品免费大片| 简卡轻食公司| 国产高清三级在线| 少妇高潮的动态图| 五月伊人婷婷丁香| 最近中文字幕高清免费大全6| 一本大道久久a久久精品| 欧美日韩亚洲高清精品| 日本wwww免费看| 少妇人妻一区二区三区视频| 成年美女黄网站色视频大全免费 | 久久久久久伊人网av| 男人爽女人下面视频在线观看| 黑人高潮一二区| 中文欧美无线码| 最近中文字幕2019免费版| 午夜免费观看性视频| 少妇 在线观看| 黄色毛片三级朝国网站 | 久久 成人 亚洲| 桃花免费在线播放| 亚洲熟女精品中文字幕| 久久国产精品大桥未久av | 我要看日韩黄色一级片| 99久国产av精品国产电影| 国产精品久久久久久久久免| 日韩成人伦理影院| 日本与韩国留学比较| 搡老乐熟女国产| 国产精品熟女久久久久浪| 熟妇人妻不卡中文字幕| 国产精品一区二区三区四区免费观看| 午夜av观看不卡| 国产男女内射视频| av在线播放精品| 亚洲成人一二三区av| 秋霞伦理黄片| 国内少妇人妻偷人精品xxx网站| 日本vs欧美在线观看视频 | 波野结衣二区三区在线| 视频区图区小说| 偷拍熟女少妇极品色| 男人爽女人下面视频在线观看| 中文字幕免费在线视频6| 高清视频免费观看一区二区| 亚洲激情五月婷婷啪啪| 高清av免费在线| 一级毛片我不卡| 欧美日韩亚洲高清精品| 日本黄色日本黄色录像| 一区在线观看完整版| 啦啦啦视频在线资源免费观看| 建设人人有责人人尽责人人享有的| 国产色婷婷99| 欧美三级亚洲精品| 亚洲av成人精品一二三区| 午夜91福利影院| 夫妻午夜视频| 丰满迷人的少妇在线观看| 国产高清国产精品国产三级| 亚洲av综合色区一区| 蜜桃久久精品国产亚洲av| 欧美变态另类bdsm刘玥| 亚洲精品乱码久久久久久按摩| 色婷婷av一区二区三区视频| 91在线精品国自产拍蜜月| 国产精品国产三级国产专区5o| 久久精品熟女亚洲av麻豆精品| 亚洲成人av在线免费| 亚洲欧美成人综合另类久久久| 一级毛片电影观看| 久热久热在线精品观看| 18禁动态无遮挡网站| 国产69精品久久久久777片| 麻豆精品久久久久久蜜桃| 大片免费播放器 马上看| 国产综合精华液| 国产av国产精品国产| 日韩欧美 国产精品| 日本-黄色视频高清免费观看| 国产白丝娇喘喷水9色精品| 高清视频免费观看一区二区| 99久久人妻综合| 熟妇人妻不卡中文字幕| 秋霞在线观看毛片| 人妻系列 视频| 婷婷色av中文字幕| 99热这里只有是精品在线观看| a 毛片基地| 国产在视频线精品| 精品久久久精品久久久| 99久久综合免费| 国产精品久久久久久av不卡| 视频中文字幕在线观看| 边亲边吃奶的免费视频| 日韩精品免费视频一区二区三区 | 在线观看人妻少妇| 亚洲精品一区蜜桃| 九九爱精品视频在线观看| 亚洲成人一二三区av| 最近的中文字幕免费完整| 亚洲综合精品二区| 亚洲精品国产成人久久av| 国产成人91sexporn| 精品久久久久久久久亚洲| 国产成人精品无人区| 三级经典国产精品| 99re6热这里在线精品视频| 久久这里有精品视频免费| 国产深夜福利视频在线观看| 色婷婷久久久亚洲欧美| 春色校园在线视频观看| 国产成人免费无遮挡视频| 91久久精品国产一区二区三区| 中国三级夫妇交换| 久久久久精品久久久久真实原创| av免费观看日本| 日本与韩国留学比较| 欧美精品一区二区大全| 久久女婷五月综合色啪小说| www.av在线官网国产| 日韩一区二区视频免费看| 亚洲美女搞黄在线观看| 99热这里只有精品一区| av黄色大香蕉| 黑人巨大精品欧美一区二区蜜桃 | 国产av一区二区精品久久| 国产日韩欧美在线精品| 免费人成在线观看视频色| 特大巨黑吊av在线直播| 老司机影院成人| 国产高清国产精品国产三级| 午夜av观看不卡| 中文资源天堂在线| 亚洲国产色片| 成人特级av手机在线观看| 国产黄色免费在线视频| 纯流量卡能插随身wifi吗| 卡戴珊不雅视频在线播放| 十分钟在线观看高清视频www | 三级国产精品片| 人妻系列 视频| 丰满饥渴人妻一区二区三| 亚洲国产欧美日韩在线播放 | 五月开心婷婷网| 亚洲第一区二区三区不卡| 国产永久视频网站| 国产黄频视频在线观看| 国产精品偷伦视频观看了| 日本色播在线视频| 亚洲av.av天堂| 亚洲婷婷狠狠爱综合网| 国产深夜福利视频在线观看| 九九在线视频观看精品| 久久精品国产自在天天线| 天美传媒精品一区二区| 亚洲精品乱码久久久v下载方式| 国产精品一区二区三区四区免费观看| 亚洲伊人久久精品综合| 久久午夜福利片| 国产乱人偷精品视频| 国产精品一区二区三区四区免费观看| www.色视频.com| 国精品久久久久久国模美| 亚洲精品乱码久久久v下载方式| 女性生殖器流出的白浆| 九九久久精品国产亚洲av麻豆| 男男h啪啪无遮挡| 在线亚洲精品国产二区图片欧美 | 国产伦在线观看视频一区| 男女边摸边吃奶| 五月玫瑰六月丁香| 亚洲综合色惰| 日韩制服骚丝袜av| 少妇人妻久久综合中文| 99九九线精品视频在线观看视频| 国产精品久久久久久久久免| 亚洲第一区二区三区不卡| 日韩大片免费观看网站| 午夜激情福利司机影院| 国产成人一区二区在线| 一区二区av电影网| 最近的中文字幕免费完整| 在线精品无人区一区二区三| 最近最新中文字幕免费大全7| 一本久久精品| 国产在线视频一区二区| 各种免费的搞黄视频| 曰老女人黄片| 五月开心婷婷网| 观看免费一级毛片| 一个人免费看片子| 精品人妻偷拍中文字幕| 亚洲va在线va天堂va国产| 精品视频人人做人人爽| 欧美精品一区二区免费开放| 欧美 亚洲 国产 日韩一| 国产精品国产av在线观看| 欧美日韩精品成人综合77777| 六月丁香七月| 国产一区二区三区综合在线观看 | 人妻制服诱惑在线中文字幕| 搡女人真爽免费视频火全软件| 最近的中文字幕免费完整| 亚洲欧美日韩东京热| 久久久久久久久大av| 黄色视频在线播放观看不卡| 欧美激情极品国产一区二区三区 | 高清午夜精品一区二区三区| 搡老乐熟女国产| 久久久久精品性色| 女人久久www免费人成看片| 一本色道久久久久久精品综合| 亚洲精品第二区| 大香蕉久久网| 成人漫画全彩无遮挡| 在线亚洲精品国产二区图片欧美 | 国产探花极品一区二区| 91aial.com中文字幕在线观看| 免费大片18禁| 夜夜骑夜夜射夜夜干| 一个人免费看片子| 日本-黄色视频高清免费观看| 99热国产这里只有精品6| 久久青草综合色| 51国产日韩欧美| 亚洲av不卡在线观看| tube8黄色片| 青春草视频在线免费观看| 赤兔流量卡办理| 亚洲国产精品999| 色网站视频免费| 日韩人妻高清精品专区| 少妇人妻一区二区三区视频| 国产亚洲91精品色在线| av国产久精品久网站免费入址| 日本黄色片子视频| 国产精品麻豆人妻色哟哟久久| 少妇高潮的动态图| 国产精品一区二区在线观看99| 精品人妻偷拍中文字幕| 免费高清在线观看视频在线观看| 最近最新中文字幕免费大全7| 久久午夜福利片| 久久久久久久国产电影| 午夜福利,免费看| 国产精品成人在线| 亚洲av电影在线观看一区二区三区| 亚洲人与动物交配视频| 亚洲丝袜综合中文字幕| 欧美精品人与动牲交sv欧美| 久久99热这里只频精品6学生| 又黄又爽又刺激的免费视频.| av播播在线观看一区| 一区二区三区免费毛片| 97超碰精品成人国产| www.色视频.com| 日韩一区二区三区影片| 黑丝袜美女国产一区| 9色porny在线观看| 一个人免费看片子| 3wmmmm亚洲av在线观看| 下体分泌物呈黄色| 久久久久久久久久人人人人人人| 97在线视频观看| 成人美女网站在线观看视频| 最黄视频免费看| 人妻少妇偷人精品九色| 亚洲中文av在线| 免费av中文字幕在线| 纯流量卡能插随身wifi吗| 亚洲av中文av极速乱| 日韩欧美精品免费久久| 一个人看视频在线观看www免费| av又黄又爽大尺度在线免费看| 丝袜在线中文字幕| 亚洲欧美一区二区三区黑人 | 中文在线观看免费www的网站| 精华霜和精华液先用哪个| 建设人人有责人人尽责人人享有的| av网站免费在线观看视频| av不卡在线播放| 精品国产露脸久久av麻豆| 毛片一级片免费看久久久久| 精品酒店卫生间| 精品少妇久久久久久888优播| 一级毛片aaaaaa免费看小| 精品久久久久久电影网| 嘟嘟电影网在线观看| 日韩一本色道免费dvd| 精品国产国语对白av| 日韩一区二区视频免费看| 日韩中文字幕视频在线看片| 午夜影院在线不卡| xxx大片免费视频| 十八禁高潮呻吟视频 | 亚洲av在线观看美女高潮| 自拍偷自拍亚洲精品老妇| 亚洲四区av| 精品人妻熟女毛片av久久网站| 有码 亚洲区| 国产美女午夜福利| 国产成人精品一,二区| 欧美最新免费一区二区三区| 菩萨蛮人人尽说江南好唐韦庄| 另类精品久久| 免费观看a级毛片全部| 青春草国产在线视频| 日韩av不卡免费在线播放| 丝袜脚勾引网站| 成人毛片60女人毛片免费| 亚洲四区av| 国产av精品麻豆| 热re99久久国产66热| 女人久久www免费人成看片| 亚洲,一卡二卡三卡| 99re6热这里在线精品视频| 国产精品女同一区二区软件| 一个人看视频在线观看www免费| 一级毛片久久久久久久久女| 国产精品人妻久久久久久| 色哟哟·www| 精品视频人人做人人爽| 亚洲无线观看免费| 国产国拍精品亚洲av在线观看| 99久久人妻综合| 国产真实伦视频高清在线观看| 国产69精品久久久久777片| 免费人妻精品一区二区三区视频| 国产精品一二三区在线看| 国产黄色免费在线视频| 免费看av在线观看网站| 日韩伦理黄色片| 色婷婷久久久亚洲欧美| 色网站视频免费| 亚洲,一卡二卡三卡| a级片在线免费高清观看视频| 国产高清三级在线| 色婷婷av一区二区三区视频| 九九爱精品视频在线观看| 国产免费一级a男人的天堂| 久久久国产一区二区| 国产在视频线精品| 涩涩av久久男人的天堂| 你懂的网址亚洲精品在线观看| 日韩三级伦理在线观看| av不卡在线播放| 日韩电影二区| 人妻夜夜爽99麻豆av| 国产一区亚洲一区在线观看| 男女边吃奶边做爰视频| 99久国产av精品国产电影| tube8黄色片| 人体艺术视频欧美日本| 亚洲怡红院男人天堂| 日日摸夜夜添夜夜爱| 久久精品国产鲁丝片午夜精品| 欧美日韩综合久久久久久| 啦啦啦在线观看免费高清www| 国内精品宾馆在线| 日韩电影二区| freevideosex欧美| 91在线精品国自产拍蜜月| 亚洲美女搞黄在线观看| 国产精品一区二区三区四区免费观看| 亚洲美女视频黄频| 久久精品久久精品一区二区三区| 国产深夜福利视频在线观看| 日本-黄色视频高清免费观看| 亚洲精品国产色婷婷电影| 极品教师在线视频| 日本免费在线观看一区| 极品人妻少妇av视频| 黄色视频在线播放观看不卡| 日本黄色日本黄色录像| 最近最新中文字幕免费大全7| 午夜激情福利司机影院| 女人久久www免费人成看片| 国产精品人妻久久久影院| 婷婷色综合www| 国产在线男女| 美女国产视频在线观看| 黑人巨大精品欧美一区二区蜜桃 | 免费黄网站久久成人精品| 欧美日韩在线观看h| 国产69精品久久久久777片| 国产一区有黄有色的免费视频| 伊人亚洲综合成人网| 热99国产精品久久久久久7| 国产有黄有色有爽视频| av卡一久久| 一个人看视频在线观看www免费| 免费大片18禁| 成人综合一区亚洲| 伦理电影免费视频| 蜜桃在线观看..| 亚洲精品国产色婷婷电影| 久热这里只有精品99| 中文精品一卡2卡3卡4更新| 国产一区有黄有色的免费视频| 我的女老师完整版在线观看| 国产在线一区二区三区精| 日韩精品有码人妻一区| 国产深夜福利视频在线观看| 免费不卡的大黄色大毛片视频在线观看| 亚洲美女视频黄频| 97精品久久久久久久久久精品| 国产精品国产三级专区第一集| 嘟嘟电影网在线观看| 日本wwww免费看| 亚洲第一av免费看| 国产日韩一区二区三区精品不卡 | 国产伦在线观看视频一区| 久久青草综合色| 色网站视频免费| 曰老女人黄片| 一区在线观看完整版| 99久久人妻综合| 午夜激情久久久久久久| 国产成人一区二区在线| 精品国产国语对白av| 国产免费一级a男人的天堂| 麻豆乱淫一区二区| 亚洲欧洲国产日韩| 男女无遮挡免费网站观看| 男女边吃奶边做爰视频| 成人美女网站在线观看视频| 亚洲av成人精品一二三区| 老司机影院毛片| 精品熟女少妇av免费看| 9色porny在线观看| 久久国产精品男人的天堂亚洲 | 国产精品嫩草影院av在线观看| 久久久久久久久久成人| 成人特级av手机在线观看| 国产日韩欧美在线精品| 久久久久精品久久久久真实原创| 少妇人妻精品综合一区二区| 亚洲第一区二区三区不卡| 免费观看无遮挡的男女| 国产精品99久久久久久久久| 高清av免费在线| 亚洲精品一区蜜桃| 观看av在线不卡| 亚洲不卡免费看| 久久人人爽人人爽人人片va| 久久久久久久久久成人| 日韩免费高清中文字幕av| 亚洲精品亚洲一区二区| 熟妇人妻不卡中文字幕| 九色成人免费人妻av| 午夜福利,免费看| 精品卡一卡二卡四卡免费| 天天操日日干夜夜撸| 男人狂女人下面高潮的视频| 桃花免费在线播放| 亚洲精品一二三| 嫩草影院入口| 久久久亚洲精品成人影院| 免费看日本二区| 国产精品国产三级国产av玫瑰| 久久久欧美国产精品| 91精品国产国语对白视频| 午夜福利在线观看免费完整高清在| 麻豆乱淫一区二区| 亚洲伊人久久精品综合| 午夜免费男女啪啪视频观看| 大片免费播放器 马上看| av网站免费在线观看视频| 欧美 日韩 精品 国产| 国产精品人妻久久久影院| 人妻一区二区av| 一区二区三区免费毛片| 亚洲av国产av综合av卡| 亚洲精品一区蜜桃| 国产视频首页在线观看| 久久国产精品大桥未久av | 最近最新中文字幕免费大全7| 国产午夜精品一二区理论片| 精品少妇内射三级| 日韩av在线免费看完整版不卡| 国产91av在线免费观看| √禁漫天堂资源中文www| 少妇高潮的动态图| 欧美精品亚洲一区二区| 久久精品久久久久久噜噜老黄| 国产淫语在线视频| 久久精品久久精品一区二区三区| .国产精品久久| 51国产日韩欧美| 一区二区av电影网| 亚洲人成网站在线观看播放| 亚洲婷婷狠狠爱综合网| freevideosex欧美| 丁香六月天网| 久久综合国产亚洲精品| av天堂中文字幕网| 最近手机中文字幕大全| 日韩成人伦理影院| 国产免费福利视频在线观看| 免费观看的影片在线观看| 久久久久久久久久人人人人人人| 国语对白做爰xxxⅹ性视频网站| 插逼视频在线观看| 中文字幕精品免费在线观看视频 | 亚洲精品国产av成人精品| 永久网站在线| 老司机影院成人| 一级毛片aaaaaa免费看小| 国产精品三级大全| 亚洲精品一二三| 青青草视频在线视频观看| 国产探花极品一区二区| 久久韩国三级中文字幕| 99精国产麻豆久久婷婷| 色婷婷久久久亚洲欧美| 日韩一本色道免费dvd| 三级国产精品片| 婷婷色综合大香蕉| 老熟女久久久| 在线观看三级黄色| 99久久精品一区二区三区| av一本久久久久| 亚洲,欧美,日韩| 久久6这里有精品| 黄色日韩在线| 春色校园在线视频观看| 美女大奶头黄色视频| 一本一本综合久久| 男人舔奶头视频| 欧美成人精品欧美一级黄| 熟妇人妻不卡中文字幕| 日韩精品有码人妻一区| 三上悠亚av全集在线观看 | 一本大道久久a久久精品| 日本黄大片高清| 日韩制服骚丝袜av| 成人二区视频| 免费播放大片免费观看视频在线观看| 九九爱精品视频在线观看| 桃花免费在线播放| 国产成人精品一,二区| 日韩av免费高清视频| 日韩电影二区| 国国产精品蜜臀av免费| a级片在线免费高清观看视频| 我要看日韩黄色一级片| 伦精品一区二区三区| 久久97久久精品| 三级经典国产精品| 亚洲国产精品一区三区| 免费少妇av软件| 一级a做视频免费观看| 国产伦精品一区二区三区四那| 久久久久人妻精品一区果冻| 免费看光身美女| 国产免费视频播放在线视频| 精品国产露脸久久av麻豆| av不卡在线播放| 国产成人精品久久久久久| 中文乱码字字幕精品一区二区三区| 亚洲精品aⅴ在线观看| 乱系列少妇在线播放| 这个男人来自地球电影免费观看 | 国产精品免费大片| 亚洲,欧美,日韩| 人妻人人澡人人爽人人| 欧美精品人与动牲交sv欧美| 国产一级毛片在线| 爱豆传媒免费全集在线观看| 在线播放无遮挡| 美女福利国产在线| 亚洲欧美一区二区三区黑人 | 在线播放无遮挡| 青青草视频在线视频观看| 少妇高潮的动态图| 少妇猛男粗大的猛烈进出视频|