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

    蒸汽浸沒射流與凝結(jié)數(shù)值仿真分析

    2019-08-23 08:00:40汪廣懷孫露王暢郝銳王建軍
    中國艦船研究 2019年4期
    關(guān)鍵詞:冷凝射流蒸汽

    汪廣懷,孫露,王暢,郝銳,王建軍

    1中國艦船研究設(shè)計中心,湖北武漢430064

    2哈爾濱工程大學核科學與技術(shù)學院,黑龍江哈爾濱150001

    0 引 言

    蒸汽浸沒射流是汽水直接接觸冷凝的一種形式,因其良好的傳熱傳質(zhì)特性,在制冷、石油化工、航天航空及核能等工業(yè)領(lǐng)域的應(yīng)用較為廣泛[1]。船用蒸汽動力系統(tǒng)進行蒸汽排放時,也是通過射流冷凝的方式實現(xiàn)高效換熱。高速蒸汽射入冷水而形成的穩(wěn)定連續(xù)蒸汽區(qū)即汽羽,在穩(wěn)定射流中形成的汽羽能夠反映浸沒射流的傳熱特性,所以國內(nèi)外開展了大量的研究工作。曹巖等[2]針對0.2~0.5 MPa飽和蒸汽在20~70℃過冷水中超音速浸沒射流凝結(jié)而形成的汽羽形狀開展了實驗研究,其根據(jù)汽羽膨脹的次數(shù)提出了5類汽羽形狀,得到了汽羽的穿透深度隨著蒸汽壓力增大和過冷水溫度上升而逐漸增加的結(jié)論。Zhou等[3]運用CFD方法對超音速浸沒射流冷凝進行數(shù)值分析,得到了汽羽在凝水中的貫穿長度和膨脹比將隨著環(huán)境壓力的增加而降低的結(jié)論,同時還發(fā)現(xiàn)噴管出口處汽相的體積波動將隨著環(huán)境壓力的下降而上升。袁方等[4]針對音速與超音速蒸汽浸沒射流核心汽羽區(qū)的穿透長度開展實驗研究,得到了汽羽的無量綱穿透長度將隨著蒸汽質(zhì)量流率的增大和過冷水溫度的上升而逐漸增加的結(jié)論,同時發(fā)現(xiàn)超音速蒸汽射流的無量綱穿透長度小于音速蒸汽射流。王超等[5]研究了超音速氣體射流在水中的噴射過程,并基于流體體積函數(shù)(Volume of Fluid,VOF)方法對浸沒射流過程開展二維仿真分析,發(fā)現(xiàn)了射流出口處氣泡生長的周期性及伴隨的脹鼓現(xiàn)象。Isaev等[6]采用修正的SST湍流模型,基于 Zapryagaev等[7]的實驗結(jié)果,精確模擬了浸沒射流汽羽,發(fā)現(xiàn)與實驗結(jié)果高度吻合。

    綜上所述,壓力、凝水溫度、噴管口徑、浸沒深度等因素對浸沒射流冷凝的傳熱傳質(zhì)特性影響較大。目前的研究成果主要與汽羽形態(tài)、長度等空間分布有關(guān),然而,對其形成過程中汽羽內(nèi)部速度場、壓力場等流體細節(jié)卻沒有形成清晰的認識。本文擬基于前期的蒸汽浸沒射流實驗成果,利用CFD仿真軟件ANSYS Fluent進行仿真模擬,分析蒸汽進口壓力、凝水溫度等因素對浸沒射流汽羽的影響,同時研究汽羽內(nèi)部的速度場和壓力場分布情況,用以為浸沒射流的高效換熱應(yīng)用場景設(shè)計提供指導。

    1 CFD模型

    ANSYS Fluent中的Mixture兩相流模型是基于Euler-Euler兩相流方程開發(fā)的應(yīng)用,其遵守的基本控制方程如下。

    1)連續(xù)性方程:

    式中:ρm為混合相密度;t為時間;vm為混合相基于質(zhì)量平均的速度。

    其中,

    式中:k=1,2,…,n,為混合相中的不同組分;αk為k相的體積分數(shù);ρk為k相密度;vk為k相速度。

    2)動量方程:

    式中:P為壓力;μm為混合相的粘度;g為重力加速度;F為體積力;vm,k為k相相對于混合相的漂移速度。

    其中,

    式中:μk為k相的粘度。

    3)漂移速度:

    使用漂移流模型時,1.2節(jié)動量方程的漂移速度可按照如下模型定義:

    式中:vpq為p,q兩相間的滑移速度;ck=αkρk/ρm,為k相對應(yīng)的質(zhì)量份額;vqk為主相q相對于混合相中k相的速度。

    在Fluent中,相間滑移速度的Manninen模型定義如下:

    其中,

    式中:dp為次相氣泡的直徑;μq為主相的粘度。

    基于Schiller-Naumann模型,fdrag的表達式為

    式中,Re為雷諾數(shù)。

    4)能量方程:

    需要說明的是,本文的相間冷凝換熱系數(shù)需通過實驗結(jié)果擬合的經(jīng)驗關(guān)系式進行自定義。

    5)兩熱阻模型:

    在相間傳熱過程中,鑒于相界面兩側(cè)分別對應(yīng)不同的換熱系數(shù),本文將相間換熱模型轉(zhuǎn)換為兩熱阻模型。定義主相(q)與次相(p)之間相界面的溫度為Ts,則相界面?zhèn)鬟f到q相和p相的熱量Qq和Qp分別為:

    式中:hq和hp分別為q相和p相的換熱系數(shù);Ai為相界面面積;Tq和Tp為q相和p相的溫度;m˙pq為由p相轉(zhuǎn)變?yōu)閝相的質(zhì)量速率;Hqs和 Hps分別為q相和p相的的焓值。

    由于相間傳遞的熱量之和為0,則相間溫度為

    相間傳熱的兩熱阻模型如下:

    式中:hpq為p,q兩相的相間傳熱系數(shù)。

    6)冷凝換熱關(guān)系式

    本文利用試驗數(shù)據(jù)擬合得到的冷凝換熱試驗關(guān)系式為

    其中,

    式中:cp,w為水的定壓比熱,kJ/(kg·℃);hfg為蒸汽潛熱;ΔT為冷凝前后水的溫差。

    通過數(shù)據(jù)對比可知,冷凝換熱系數(shù)的計算結(jié)果(式(16))與實驗結(jié)果的最大偏差在±12%以內(nèi)。

    2 幾何模型及仿真模型驗證

    2.1 實驗裝置

    在前期研究工作中,為分析蒸汽浸沒射流與凝結(jié)換熱現(xiàn)象,搭建了如圖1所示的實驗系統(tǒng),該系統(tǒng)主要包括電加熱鍋爐、蒸汽噴嘴、水箱、配置了壓力和溫度探針的三維測量支架、可視化系統(tǒng)和數(shù)據(jù)采集系統(tǒng)[8]。

    圖1 蒸汽浸沒射流實驗系統(tǒng)Fig.1 Experiment system for steam jet in water

    通過調(diào)節(jié)鍋爐功率,即可獲得不同參數(shù)的蒸汽;基于水箱的充水排水系統(tǒng),即可調(diào)節(jié)水箱內(nèi)部的溫度;利用高速攝影儀,即可對不同工況下的射流冷凝現(xiàn)象進行可視化攝像。該系統(tǒng)的實驗參數(shù)如表1所示。

    表1 實驗參數(shù)Table 1 Experimental parameters

    2.2 計算模型

    結(jié)合實驗過程,將實驗裝置簡化為如圖2所示的計算模型。其中,噴管的上端進口設(shè)置為給定壓力邊界的壓力進口,水箱的水面設(shè)置為壓力出口邊界(壓力設(shè)為1個標準大氣壓),水箱的其余壁面則設(shè)置為壁面無滑移邊界條件。

    水箱的幾何尺寸(長×寬×高)為1.0 m×0.5 m×2.0 m,噴管直徑為10 mm,浸沒深度為0.5 m。針對簡化的三維幾何模型,本文采用了結(jié)構(gòu)化網(wǎng)格劃分方法,其空間離散網(wǎng)格如圖3所示,其中噴管中心部分采用O-Block技術(shù)進行了加密處理。

    圖2 計算模型幾何結(jié)構(gòu)示意圖Fig.2 Schematic of geometry structure of calculation model

    圖3 網(wǎng)格劃分示意圖Fig.3 Display of mesh generation

    基于該網(wǎng)格拓撲結(jié)構(gòu),本文將針對4套不同基準尺度的離散網(wǎng)格(基準尺寸為10~40 mm),在入口壓力為0.2 MPa、冷卻劑溫度為25℃的條件下進行蒸汽射流仿真計算,以驗證網(wǎng)格的獨立性。在4套計算網(wǎng)格下,噴管內(nèi)蒸汽的最高射流流速如圖4所示。當網(wǎng)格尺度低于30 mm時,蒸汽射流流速的最大偏差小于0.5%,從而驗證了對應(yīng)網(wǎng)格的獨立性。因此,本文將網(wǎng)格尺度設(shè)為20 mm,而網(wǎng)格整體規(guī)模約為8×105。

    圖4 網(wǎng)格無關(guān)性驗證Fig.4 Verification of grid independence

    蒸汽浸沒射流的直接接觸冷凝是一個瞬態(tài)的物理過程,為保證仿真模擬的計算精度,本文將采用雙精度求解器。同時,仿真計算將考慮采用兩相相間滑移的Mixture兩相流模型,而相間傳熱傳質(zhì)將通過在兩相的體積源項中使用DEFINE_PROPERTY UDF宏添加udf進行自定義。進行壓力—速度的耦合計算時,將采用壓力的隱式算子分割算法(Pressure Implicit with Splitting of Operators,PISO);守恒方程的離散格式中,體積分數(shù)項將采用QUICK格式,其余項采用二階迎風格式。瞬態(tài)計算中:時間步長為10-5s,計算的物理時間為100 ms,每一個計算時間步內(nèi)需要完成50~150次迭代;采用ANSYS Fluent 14.0仿真程序,服務(wù)器的基本配置為 E5 2685 v3(×2)CPU。在此基礎(chǔ)上,完成一個算例約耗時96 h。

    本文需保證連續(xù)性方程的收斂殘差降至10-3量級,動量方程收斂至10-6量級以下,同時需監(jiān)控蒸汽排放管出口速度在每一個瞬態(tài)時間步中迭代的變化,并以此作為下一時間步中迭代計算的輸入。綜合上述收斂評判標準,即可保證在一個瞬態(tài)計算步的迭代滿足收斂精度。

    2.3 仿真模型驗證

    在蒸汽浸沒射流凝結(jié)過程中,蒸汽由噴管射入水中,典型的穩(wěn)定汽羽形狀主要包括錐形和橢球形2種,如圖5所示。其中,圖5(a)~圖5(d)為豎直噴射時不同形狀汽羽的仿真結(jié)果,圖5(e)~圖5(h)為對應(yīng)的實驗結(jié)果[8]。

    在低蒸汽流量、低凝水溫度條件下,蒸汽在噴管出口處充分膨脹,然后進入凝水主流中得到充分冷凝,整個汽羽呈現(xiàn)出如圖5(a)所示的錐形汽羽。隨著蒸汽流量增加,汽羽核心區(qū)的蒸汽無法充分冷凝,故有所膨脹,從而形成如圖5(b)所示的橢球形汽羽。當蒸汽流量進一步增加時,汽羽在凝水中膨脹后收縮,但未能及時冷卻的高溫汽羽將進一步膨脹,形成如圖5(c)所示的雙膨脹收縮形汽羽。在低蒸汽流量、高凝水溫度條件下,由于氣泡產(chǎn)生和脫落,噴射蒸汽將在凝水中呈現(xiàn)如圖5(d)所示的振蕩流形。

    由圖5可知,兩熱阻模型作為相變模型,能夠較好地模擬不同類型的蒸汽射流凝結(jié)實驗,且仿真計算結(jié)果與實驗結(jié)果[8]吻合度較高,驗證了該模型的適用性及準確性。

    3 射流汽羽影響分析

    基于上述仿真模型,本文將進一步分析冷凝溫度、入口壓力等因素對當前實驗條件下浸沒射流汽羽行為的影響。

    圖5 仿真結(jié)果與實驗結(jié)果對比Fig.5 Comparison between simulation results and experimental results

    3.1 汽羽長度確定

    一般而言,在實驗條件下可以通過調(diào)節(jié)可視化圖像中氣、液兩相之間的對比度來確定汽羽長度[9]。以噴管直徑10 mm、浸沒深度0.5 m、水溫25℃、入口壓力0.5 MPa條件下的計算結(jié)果為例,參照文獻[10]截取計算結(jié)果中的汽羽,如圖6所示,橫坐標為試驗裝置中水箱液位的不同高度。本文將以汽相體積份額為0.95作為汽羽長度度量的截斷值,如圖6所示,即以兩相區(qū)域內(nèi)汽相份額為0.95的輪廓線作為形成汽羽的邊界。

    3.2 汽羽形狀對比

    圖6 汽羽長度的計算方法Fig.6 Calculation method of steam plume length

    根據(jù)唐思等[11]的研究,入口壓力和水溫將對蒸汽排入過冷水中的膨脹—壓縮過程及蒸汽的冷凝效率產(chǎn)生影響,從而導致蒸汽汽羽形狀和流場有所區(qū)別。如圖7(a)和圖7(b)所示,蒸汽入口壓力從0.3 MPa變化到0.8 MPa,當水溫為25℃和35℃時,汽羽形狀呈現(xiàn)出與實驗結(jié)果相似的橢球形。

    隨著蒸汽入口處壓力的增加,汽羽的膨脹程度也隨之增加,具體表現(xiàn)為:噴管出口處的汽羽直徑增加,蒸汽的射流速度增加,并導致汽羽長度進一步增加。如圖7(c)所示,在水溫為50℃條件下,當入口壓力由0.3 MPa增加至0.8 MPa時,汽羽形狀逐漸由橢球形轉(zhuǎn)變?yōu)殡p膨脹收縮形。隨著入口壓力的增加,在汽羽內(nèi)部首先發(fā)生第1次膨脹—壓縮過程,但蒸汽凝結(jié)不充分,進而發(fā)生第2次膨脹—壓縮過程,最終呈現(xiàn)出雙膨脹收縮型汽羽。

    圖7 不同溫度、壓力條件下的汽羽形狀對比Fig.7 Comparison of steam plumes under different pressures and temperatures

    3.3 溫度場分布

    徑向溫度場分布是蒸汽射流凝結(jié)的一個重要特征。圖8所示為入口蒸汽壓力為0.5 MPa、冷凝水溫為25℃的工況下,在噴嘴出口不同距離處截面的徑向溫度分布情況。圖8中,x為距離噴管中心軸線的徑向距離,x=0即為噴管軸線處。

    如圖8(a)所示,噴管出口截面上的蒸汽溫度分布與管內(nèi)相似,其中管壁位置出現(xiàn)了較高的汽相溫度分布。這是由于蒸汽噴射到凝水界面時與凝水產(chǎn)生劇烈壓縮,從而導致了對應(yīng)位置處的蒸汽溫度上升。

    如圖8(b)所示,汽羽核心處存在低溫區(qū),由徑向溫度分布曲線可知,蒸汽由于膨脹做功而導致溫度降低。汽羽內(nèi)部沿著徑向溫度是先上升后下降,并逐漸接近凝水溫度。由可壓縮氣體的紊流射流可知,靠近氣液界面處的蒸汽首先處于壓縮狀態(tài),蒸汽的內(nèi)能增加,溫度升高;然后,由于卷吸和凝結(jié)作用,蒸汽溫度逐漸降低,并接近周圍過冷水的溫度。

    圖8 溫度場分布Fig.8 Distribution of temperature field

    在圖8(c)和圖8(d)中,由于汽相受到壓縮,且對應(yīng)位置在壓縮波處,因此均在汽羽軸線達到了峰值溫度。

    3.4 壓力場分布

    圖9所示為10 mm管徑下水溫為50℃時的壓力分布情況。圖9(a)為蒸汽射流流場中心截面的壓力分布,可以發(fā)現(xiàn),蒸汽射流流場中壓力場的分布大致可以分為3個區(qū):對應(yīng)汽羽內(nèi)部汽空間范圍為壓力核心區(qū),其壓力梯度較大;對應(yīng)包裹汽羽汽空間的氣液混合區(qū)為壓力擴散區(qū),其表壓為負值;對應(yīng)汽羽外側(cè)的凝水區(qū)為壓力邊緣區(qū),其接近于環(huán)境壓力。在壓力擴散區(qū),由于蒸汽出現(xiàn)了部分冷凝,故其壓力低于凝水的環(huán)境壓力;在接近液體單相區(qū),其壓力逐漸接近環(huán)境壓力。

    圖9(b)為不同蒸汽壓力進口條件下,距離噴嘴出口0.1 m處的徑向壓力分布圖。由中心位置沿徑向方向,壓力呈現(xiàn)出先降低后升高,然后逐漸接近周圍環(huán)境壓力的變化趨勢。

    圖9(b)中的壓力谷值主要源自蒸汽射流的卷吸過程,卷吸過程會強化蒸汽的凝結(jié)作用,導致局部壓力降低。結(jié)合圖9(a)和圖9(b)可知,壓力場沿噴管軸線呈軸對稱形式,且在中軸線處達到壓力峰值。隨著入口初始壓力的增加,該峰值也隨之增加,且分布曲線的影響范圍逐漸擴大,表現(xiàn)為蒸汽汽羽直徑逐漸增大。另一方面,軸線處的壓力峰值相對于蒸汽入口壓力的壓差較大,表明蒸汽射流在凝水中的傳質(zhì)過程較為強烈,從而導致動量傳遞也較強烈。

    圖9(c)和圖9(d)所示為不同蒸汽壓力進口條件下,沿噴管軸線方向的壓力分布曲線??梢钥闯?,汽羽內(nèi)部存在一定范圍的壓力波動,表現(xiàn)為壓力變化呈現(xiàn)出2個峰值和2~3個谷值。這一結(jié)果表明了汽羽內(nèi)部壓縮波和膨脹波的存在性:首先,蒸汽在噴管出口處發(fā)生膨脹過程,導致壓力過低;然后,進入壓縮過程,導致壓力升高;最后,進入壓縮—膨脹過程。隨著入口壓力的增加,第2次壓縮—膨脹過程產(chǎn)生的壓力波動逐漸增強。同時,通過對比圖9(c)和圖9(d)可知,進口壓力對汽羽內(nèi)部壓縮—膨脹波的波節(jié)數(shù)影響較小,但較大的進口壓力將導致壓縮—膨脹波的波長增加,即壓力波動將傳播更遠,其幅值也將有所增加。

    圖9 10 mm管徑下水溫為50℃時的壓力分布Fig.9 Pressure distribution under 10 mm tube and 50℃condition

    4 結(jié) 論

    基于筆者前期開展的浸沒射流實驗,本文針對浸沒射流冷凝現(xiàn)象開展了CFD仿真分析,通過引入實驗過程中擬合的傳熱經(jīng)驗關(guān)系式,得到了與實驗結(jié)果高度吻合的仿真結(jié)果。本文以汽相體積份額為0.95作為度量汽羽長度的臨界分數(shù),對比研究了不同凝水溫度和進口壓力條件下浸沒射流冷凝汽羽的空間分布、速度場、溫度場及壓力場,充分掌握了這類開放空間純蒸汽浸沒射流的基礎(chǔ)物理現(xiàn)象,得到結(jié)論如下:

    1)隨著蒸汽入口壓力的增加,蒸汽在水域中的膨脹程度隨之增加,產(chǎn)生的汽羽直徑和汽羽長度也逐漸增加。同時,不同溫度條件下,冷凝水箱內(nèi)形成的蒸汽汽羽形狀將有所不同。

    2)蒸汽進口壓力較高時,水溫對汽羽內(nèi)部蒸汽膨脹的影響較小,且不同水溫條件下沿軸線的汽羽速度分布比較接近。

    3)汽羽內(nèi)部的徑向溫度分布將呈現(xiàn)出不同的特征,由于汽羽的膨脹過程,在核心區(qū)將出現(xiàn)局部低溫現(xiàn)象。

    4)蒸汽射流氣羽周圍的壓力場分布主要分為3個區(qū):較大壓力梯度的核心區(qū)、負壓的擴散區(qū)及接近環(huán)境壓力的邊緣區(qū)。

    5)隨著蒸汽初始壓力的增加,核心區(qū)靜壓的峰值逐漸升高,蒸汽的汽羽直徑也隨之增加。

    6)汽羽內(nèi)部沿軸向的壓力分布一般存在2組波峰—波谷,這是因為汽羽內(nèi)部存在壓縮波和膨脹波,蒸汽在噴嘴出口處先發(fā)生膨脹過程,再進行壓縮過程。進口壓力對汽羽內(nèi)部壓縮—膨脹波的波節(jié)數(shù)影響較小,但會導致壓力波的波動幅值和影響范圍有所增加。

    猜你喜歡
    冷凝射流蒸汽
    核電廠蒸汽發(fā)生器一次側(cè)管嘴堵板研發(fā)和應(yīng)用
    深海逃逸艙射流注水均壓過程仿真分析
    低壓天然氣泄漏射流擴散特性研究
    煤氣與熱力(2022年4期)2022-05-23 12:45:00
    洗脫苯粗苯冷凝冷卻器優(yōu)化改造
    山東冶金(2019年6期)2020-01-06 07:46:14
    一種新型蒸汽發(fā)生器結(jié)構(gòu)設(shè)計
    中國核電(2017年1期)2017-05-17 06:10:04
    第一艘蒸汽輪船
    蒸汽閃爆
    硫磺回收裝置冷凝冷卻器改造
    射流齒形噴嘴射流流場與氣動聲學分析
    地鐵站臺活塞風附壁射流起始段的實測和實驗驗證
    制冷學報(2014年3期)2014-03-01 03:07:17
    又紧又爽又黄一区二区| 性高湖久久久久久久久免费观看| 久久精品国产综合久久久| 欧美精品人与动牲交sv欧美| 亚洲成人国产一区在线观看| 日韩制服丝袜自拍偷拍| 欧美性长视频在线观看| 色综合欧美亚洲国产小说| 丰满人妻熟妇乱又伦精品不卡| 日本av免费视频播放| 在线观看免费日韩欧美大片| 日本vs欧美在线观看视频| 欧美日韩亚洲国产一区二区在线观看 | 伊人亚洲综合成人网| 99国产精品一区二区蜜桃av | av欧美777| 一级a爱视频在线免费观看| 免费一级毛片在线播放高清视频 | 性少妇av在线| 国产区一区二久久| 免费在线观看影片大全网站| 日韩欧美国产一区二区入口| 1024视频免费在线观看| 久久人人97超碰香蕉20202| svipshipincom国产片| 人人妻,人人澡人人爽秒播| 免费观看人在逋| 又紧又爽又黄一区二区| 国产人伦9x9x在线观看| 国产精品国产三级国产专区5o| 人妻久久中文字幕网| 亚洲一码二码三码区别大吗| 美女高潮到喷水免费观看| 国产精品 欧美亚洲| 日本黄色日本黄色录像| 巨乳人妻的诱惑在线观看| av在线老鸭窝| 好男人电影高清在线观看| 丝袜在线中文字幕| 悠悠久久av| 亚洲欧美激情在线| 国产在视频线精品| 人人妻人人添人人爽欧美一区卜| 人妻久久中文字幕网| 国产成人av教育| 欧美精品一区二区免费开放| 可以免费在线观看a视频的电影网站| 熟女少妇亚洲综合色aaa.| 国产av国产精品国产| 99国产精品一区二区三区| 国产老妇伦熟女老妇高清| 99久久人妻综合| 久久精品国产亚洲av香蕉五月 | 亚洲中文日韩欧美视频| 9热在线视频观看99| 12—13女人毛片做爰片一| 搡老岳熟女国产| 亚洲精品美女久久av网站| 国产免费福利视频在线观看| 国产区一区二久久| 国产成人a∨麻豆精品| 美女国产高潮福利片在线看| 在线观看www视频免费| 午夜视频精品福利| 亚洲三区欧美一区| 国产成人系列免费观看| 99精品欧美一区二区三区四区| 黑丝袜美女国产一区| h视频一区二区三区| 热99久久久久精品小说推荐| 亚洲欧美精品自产自拍| 国产一区二区三区综合在线观看| 国产三级黄色录像| 欧美在线黄色| 69精品国产乱码久久久| 午夜福利视频精品| 中文字幕人妻丝袜一区二区| 亚洲午夜精品一区,二区,三区| 女性被躁到高潮视频| 日韩中文字幕欧美一区二区| 成年av动漫网址| 啦啦啦免费观看视频1| 黄片小视频在线播放| 91麻豆精品激情在线观看国产 | 国产xxxxx性猛交| 久久久欧美国产精品| 视频区欧美日本亚洲| 成人亚洲精品一区在线观看| 国产又色又爽无遮挡免| 久久天堂一区二区三区四区| 一级黄色大片毛片| 午夜两性在线视频| 国产成人av激情在线播放| 精品欧美一区二区三区在线| av片东京热男人的天堂| 91国产中文字幕| 女人被躁到高潮嗷嗷叫费观| 手机成人av网站| www日本在线高清视频| 免费在线观看视频国产中文字幕亚洲 | 久久久久国产精品人妻一区二区| 成年动漫av网址| 国产高清videossex| 久久国产精品大桥未久av| 午夜福利乱码中文字幕| 啦啦啦中文免费视频观看日本| 欧美国产精品一级二级三级| 黄片大片在线免费观看| 亚洲国产av影院在线观看| 90打野战视频偷拍视频| 香蕉丝袜av| 91字幕亚洲| 婷婷丁香在线五月| 亚洲久久久国产精品| 午夜两性在线视频| 国产成+人综合+亚洲专区| 免费在线观看视频国产中文字幕亚洲 | 色老头精品视频在线观看| 在线观看舔阴道视频| 啦啦啦视频在线资源免费观看| 伦理电影免费视频| 各种免费的搞黄视频| 一本色道久久久久久精品综合| 精品国产乱码久久久久久小说| 可以免费在线观看a视频的电影网站| 丁香六月欧美| 老司机影院成人| 午夜视频精品福利| 欧美日韩国产mv在线观看视频| 一区二区av电影网| 日本av免费视频播放| 亚洲av电影在线观看一区二区三区| 午夜福利影视在线免费观看| 操美女的视频在线观看| 亚洲精品久久午夜乱码| 999久久久精品免费观看国产| 性色av一级| 国产精品久久久久久精品古装| 最近最新免费中文字幕在线| 午夜成年电影在线免费观看| 一边摸一边做爽爽视频免费| 欧美老熟妇乱子伦牲交| 18在线观看网站| 国产亚洲欧美精品永久| 日韩 亚洲 欧美在线| 青春草视频在线免费观看| 成人免费观看视频高清| 中文欧美无线码| 搡老乐熟女国产| 国产视频一区二区在线看| 精品一品国产午夜福利视频| 纵有疾风起免费观看全集完整版| 最近最新中文字幕大全免费视频| 老司机午夜福利在线观看视频 | 欧美激情极品国产一区二区三区| 国内毛片毛片毛片毛片毛片| 成在线人永久免费视频| 亚洲国产av新网站| 午夜日韩欧美国产| 亚洲精品成人av观看孕妇| 少妇被粗大的猛进出69影院| 欧美日韩福利视频一区二区| 日韩免费高清中文字幕av| 精品一区二区三卡| 亚洲天堂av无毛| a级毛片在线看网站| 亚洲五月婷婷丁香| 黄色怎么调成土黄色| 成人av一区二区三区在线看 | 人人妻人人添人人爽欧美一区卜| 婷婷成人精品国产| 香蕉丝袜av| 老司机午夜十八禁免费视频| 日本vs欧美在线观看视频| 18禁观看日本| 欧美精品av麻豆av| 成人手机av| 一区二区三区精品91| 汤姆久久久久久久影院中文字幕| 日本一区二区免费在线视频| 中文精品一卡2卡3卡4更新| 欧美另类一区| 脱女人内裤的视频| 老司机亚洲免费影院| 91精品国产国语对白视频| 国产精品欧美亚洲77777| 亚洲激情五月婷婷啪啪| 亚洲国产欧美日韩在线播放| 欧美 亚洲 国产 日韩一| 亚洲精品av麻豆狂野| 可以免费在线观看a视频的电影网站| 精品国产乱码久久久久久小说| 免费女性裸体啪啪无遮挡网站| 久久中文看片网| 亚洲,欧美精品.| 国产一区二区三区综合在线观看| 999久久久国产精品视频| 亚洲专区中文字幕在线| av线在线观看网站| 少妇精品久久久久久久| 美女午夜性视频免费| a级片在线免费高清观看视频| 少妇被粗大的猛进出69影院| 国产欧美亚洲国产| 香蕉丝袜av| 免费高清在线观看视频在线观看| 亚洲av日韩精品久久久久久密| 国产精品一区二区在线观看99| 国产熟女午夜一区二区三区| 50天的宝宝边吃奶边哭怎么回事| 啦啦啦在线免费观看视频4| 国产精品久久久久久精品电影小说| 亚洲第一av免费看| 在线看a的网站| 大码成人一级视频| 一级毛片精品| 久久久久久久国产电影| 大型av网站在线播放| 欧美一级毛片孕妇| 国产精品熟女久久久久浪| 少妇精品久久久久久久| 欧美亚洲 丝袜 人妻 在线| 久久久精品94久久精品| 国产精品久久久久久精品古装| 成年动漫av网址| 久久香蕉激情| 色婷婷久久久亚洲欧美| 黑人巨大精品欧美一区二区蜜桃| 免费在线观看完整版高清| 亚洲免费av在线视频| 一区二区三区乱码不卡18| 午夜福利视频在线观看免费| 99国产精品一区二区蜜桃av | av免费在线观看网站| av欧美777| 亚洲国产欧美一区二区综合| 我的亚洲天堂| 每晚都被弄得嗷嗷叫到高潮| 国产成人精品无人区| 日韩中文字幕视频在线看片| 成年人午夜在线观看视频| 99国产极品粉嫩在线观看| 国产成人精品在线电影| 国产亚洲午夜精品一区二区久久| 一区二区三区激情视频| 国产一区二区三区综合在线观看| 欧美日韩福利视频一区二区| 丝袜脚勾引网站| tube8黄色片| 欧美精品一区二区大全| av一本久久久久| 999久久久国产精品视频| 男女免费视频国产| 高清在线国产一区| 少妇精品久久久久久久| 亚洲中文av在线| 久久性视频一级片| 午夜福利视频精品| 一边摸一边抽搐一进一出视频| 精品少妇一区二区三区视频日本电影| 久热爱精品视频在线9| 热re99久久国产66热| 精品久久蜜臀av无| 午夜老司机福利片| 三级毛片av免费| 在线av久久热| 黄色a级毛片大全视频| 老司机深夜福利视频在线观看 | 亚洲av日韩精品久久久久久密| 男女边摸边吃奶| 国产激情久久老熟女| www日本在线高清视频| 97在线人人人人妻| 久久人妻福利社区极品人妻图片| 免费在线观看黄色视频的| 国产精品 欧美亚洲| 老汉色∧v一级毛片| 欧美激情 高清一区二区三区| 国产黄频视频在线观看| 亚洲午夜精品一区,二区,三区| 日韩熟女老妇一区二区性免费视频| 亚洲第一欧美日韩一区二区三区 | 国产深夜福利视频在线观看| 国产成人影院久久av| 亚洲av电影在线观看一区二区三区| 超碰成人久久| 国产精品1区2区在线观看. | 午夜激情久久久久久久| 欧美一级毛片孕妇| 啦啦啦视频在线资源免费观看| 黄网站色视频无遮挡免费观看| 久久精品国产亚洲av香蕉五月 | 99香蕉大伊视频| 黄片大片在线免费观看| 母亲3免费完整高清在线观看| 久久精品久久久久久噜噜老黄| 国产精品偷伦视频观看了| www.精华液| 夜夜骑夜夜射夜夜干| 老熟女久久久| 久久亚洲精品不卡| 亚洲色图 男人天堂 中文字幕| 黑人欧美特级aaaaaa片| 中文精品一卡2卡3卡4更新| 欧美日韩精品网址| 我要看黄色一级片免费的| www日本在线高清视频| 国产精品久久久久成人av| 亚洲av电影在线观看一区二区三区| 正在播放国产对白刺激| 亚洲精品国产av成人精品| 久久久久久人人人人人| 欧美变态另类bdsm刘玥| 欧美激情高清一区二区三区| 国产色视频综合| 亚洲精品国产色婷婷电影| av不卡在线播放| 精品亚洲乱码少妇综合久久| 亚洲精品美女久久av网站| 中文字幕人妻丝袜制服| 国产精品秋霞免费鲁丝片| 国产99久久九九免费精品| 一二三四在线观看免费中文在| 好男人电影高清在线观看| av视频免费观看在线观看| 欧美日韩福利视频一区二区| 又黄又粗又硬又大视频| 久久久久精品人妻al黑| av一本久久久久| h视频一区二区三区| 久久久久视频综合| tube8黄色片| 精品久久久久久电影网| 国产野战对白在线观看| 午夜福利在线免费观看网站| 日日爽夜夜爽网站| 国产精品香港三级国产av潘金莲| 最黄视频免费看| 搡老熟女国产l中国老女人| 亚洲久久久国产精品| 一本—道久久a久久精品蜜桃钙片| 国产精品1区2区在线观看. | av网站免费在线观看视频| 亚洲色图综合在线观看| 两个人免费观看高清视频| 青草久久国产| 五月天丁香电影| 午夜福利免费观看在线| 欧美少妇被猛烈插入视频| www.自偷自拍.com| 日韩大码丰满熟妇| 黑人操中国人逼视频| 精品视频人人做人人爽| 纵有疾风起免费观看全集完整版| 亚洲一区中文字幕在线| 亚洲精品粉嫩美女一区| 国产亚洲欧美精品永久| 国产亚洲精品第一综合不卡| 国产色视频综合| 亚洲av片天天在线观看| 国产欧美日韩精品亚洲av| 国产男人的电影天堂91| 欧美+亚洲+日韩+国产| 伊人亚洲综合成人网| 超碰成人久久| 国产成人精品无人区| 精品国产国语对白av| 啦啦啦啦在线视频资源| 精品国产一区二区三区四区第35| 久久精品国产亚洲av高清一级| 久久久久久亚洲精品国产蜜桃av| 人妻人人澡人人爽人人| 国产亚洲精品第一综合不卡| 51午夜福利影视在线观看| 国产精品秋霞免费鲁丝片| 99国产精品一区二区三区| 精品视频人人做人人爽| 高清av免费在线| 国产精品一区二区免费欧美 | 亚洲专区中文字幕在线| 美女脱内裤让男人舔精品视频| e午夜精品久久久久久久| 国产亚洲精品一区二区www | 亚洲av国产av综合av卡| 少妇的丰满在线观看| 久久精品亚洲熟妇少妇任你| 丝袜喷水一区| 国产精品自产拍在线观看55亚洲 | 日本wwww免费看| 丰满人妻熟妇乱又伦精品不卡| 丰满少妇做爰视频| 免费日韩欧美在线观看| tube8黄色片| 国内毛片毛片毛片毛片毛片| 成人国产av品久久久| 成人亚洲精品一区在线观看| 国产av又大| 自拍欧美九色日韩亚洲蝌蚪91| 欧美精品啪啪一区二区三区 | 亚洲欧美色中文字幕在线| 久久中文看片网| 亚洲精品国产av成人精品| 每晚都被弄得嗷嗷叫到高潮| 国产一区二区激情短视频 | 国产精品一区二区精品视频观看| 午夜福利免费观看在线| 久久热在线av| h视频一区二区三区| 18禁裸乳无遮挡动漫免费视频| 搡老熟女国产l中国老女人| 青春草视频在线免费观看| 菩萨蛮人人尽说江南好唐韦庄| 秋霞在线观看毛片| 天天操日日干夜夜撸| 人人妻人人添人人爽欧美一区卜| tube8黄色片| 亚洲天堂av无毛| 欧美日韩亚洲高清精品| 黄片大片在线免费观看| 亚洲精品一卡2卡三卡4卡5卡 | 亚洲伊人久久精品综合| 久久人人爽人人片av| 免费人妻精品一区二区三区视频| 国产老妇伦熟女老妇高清| 女性生殖器流出的白浆| 青青草视频在线视频观看| 可以免费在线观看a视频的电影网站| 黄色视频在线播放观看不卡| 一个人免费在线观看的高清视频 | 三级毛片av免费| 另类精品久久| 99精国产麻豆久久婷婷| xxxhd国产人妻xxx| 免费在线观看完整版高清| 一本色道久久久久久精品综合| 午夜影院在线不卡| e午夜精品久久久久久久| 国产高清视频在线播放一区 | 国产av国产精品国产| 在线永久观看黄色视频| 各种免费的搞黄视频| 丝瓜视频免费看黄片| 国产精品av久久久久免费| 国产男人的电影天堂91| 2018国产大陆天天弄谢| 久久人人爽av亚洲精品天堂| 亚洲专区中文字幕在线| 日本黄色日本黄色录像| 日韩精品免费视频一区二区三区| 亚洲激情五月婷婷啪啪| 黄色a级毛片大全视频| av线在线观看网站| 久久人妻熟女aⅴ| 欧美日韩亚洲国产一区二区在线观看 | 亚洲av成人一区二区三| 新久久久久国产一级毛片| 狠狠狠狠99中文字幕| 中国国产av一级| 三上悠亚av全集在线观看| 国产av一区二区精品久久| 国产亚洲欧美在线一区二区| 亚洲色图 男人天堂 中文字幕| 色视频在线一区二区三区| 久久精品aⅴ一区二区三区四区| 国产成人影院久久av| 夜夜骑夜夜射夜夜干| www.999成人在线观看| 欧美中文综合在线视频| 国产精品麻豆人妻色哟哟久久| 国产av又大| 国产精品免费视频内射| 亚洲精品中文字幕一二三四区 | 在线十欧美十亚洲十日本专区| 国产无遮挡羞羞视频在线观看| 久久精品国产a三级三级三级| 国产一级毛片在线| 亚洲精品美女久久久久99蜜臀| 亚洲国产精品一区三区| 99精品欧美一区二区三区四区| 国产成人免费观看mmmm| 12—13女人毛片做爰片一| 亚洲一区二区三区欧美精品| 50天的宝宝边吃奶边哭怎么回事| 亚洲成人国产一区在线观看| 一级片'在线观看视频| 黑人巨大精品欧美一区二区蜜桃| 欧美精品一区二区免费开放| 国产精品国产三级国产专区5o| 亚洲精品国产精品久久久不卡| 国产精品一区二区免费欧美 | 久久久久久人人人人人| 国产成人影院久久av| 下体分泌物呈黄色| 黄频高清免费视频| netflix在线观看网站| 亚洲精品国产精品久久久不卡| 50天的宝宝边吃奶边哭怎么回事| 91麻豆精品激情在线观看国产 | 人妻人人澡人人爽人人| 精品久久久久久久毛片微露脸 | 亚洲精品一区蜜桃| 麻豆国产av国片精品| 国产一级毛片在线| 久久九九热精品免费| 亚洲一区中文字幕在线| 久久免费观看电影| 精品国产超薄肉色丝袜足j| 日本一区二区免费在线视频| 国产真人三级小视频在线观看| 成人18禁高潮啪啪吃奶动态图| 国产欧美亚洲国产| 欧美精品啪啪一区二区三区 | 亚洲午夜精品一区,二区,三区| 久久性视频一级片| av在线app专区| svipshipincom国产片| 十八禁网站网址无遮挡| 50天的宝宝边吃奶边哭怎么回事| 深夜精品福利| 99国产综合亚洲精品| www.熟女人妻精品国产| 久久狼人影院| 久久久久视频综合| 国产日韩欧美亚洲二区| 久久精品久久久久久噜噜老黄| 婷婷色av中文字幕| 欧美成狂野欧美在线观看| 99热全是精品| 国产黄频视频在线观看| 日韩中文字幕欧美一区二区| 欧美日韩成人在线一区二区| 亚洲精品一卡2卡三卡4卡5卡 | 成人手机av| 每晚都被弄得嗷嗷叫到高潮| 99热国产这里只有精品6| 热re99久久精品国产66热6| 久久九九热精品免费| 亚洲国产欧美日韩在线播放| 精品国产乱码久久久久久小说| 黑人巨大精品欧美一区二区蜜桃| 两人在一起打扑克的视频| 少妇精品久久久久久久| 成人三级做爰电影| 深夜精品福利| 亚洲va日本ⅴa欧美va伊人久久 | 真人做人爱边吃奶动态| 99国产综合亚洲精品| 久久精品久久久久久噜噜老黄| 国产片内射在线| 两个人免费观看高清视频| 99热网站在线观看| 久久久久久久大尺度免费视频| 国产色视频综合| 亚洲精品国产区一区二| 麻豆av在线久日| 亚洲一区中文字幕在线| 中文字幕人妻丝袜制服| 王馨瑶露胸无遮挡在线观看| 亚洲精品国产精品久久久不卡| 精品少妇久久久久久888优播| 国产又色又爽无遮挡免| 天堂8中文在线网| 91精品三级在线观看| 在线观看免费高清a一片| 丝袜美腿诱惑在线| www日本在线高清视频| 午夜福利在线观看吧| 国产精品一二三区在线看| 制服诱惑二区| 天天躁狠狠躁夜夜躁狠狠躁| 久久青草综合色| 久久久久久久久久久久大奶| 精品亚洲乱码少妇综合久久| 欧美久久黑人一区二区| 高清黄色对白视频在线免费看| 亚洲全国av大片| 国产高清国产精品国产三级| 性色av乱码一区二区三区2| 丝袜在线中文字幕| 多毛熟女@视频| 美女脱内裤让男人舔精品视频| 国产精品久久久久久精品古装| 欧美另类亚洲清纯唯美| 久久女婷五月综合色啪小说| 激情视频va一区二区三区| 婷婷丁香在线五月| 亚洲av日韩精品久久久久久密| 日韩 亚洲 欧美在线| 两人在一起打扑克的视频| 精品少妇黑人巨大在线播放| 午夜激情av网站| 女人精品久久久久毛片| 日本精品一区二区三区蜜桃| 欧美精品啪啪一区二区三区 | 久久精品国产综合久久久| 欧美精品一区二区免费开放| 国产淫语在线视频| 99九九在线精品视频| 亚洲熟女精品中文字幕| 久久精品aⅴ一区二区三区四区| 搡老熟女国产l中国老女人| 后天国语完整版免费观看| 久热这里只有精品99| 黄片播放在线免费| 国产av精品麻豆| 国产精品免费视频内射| 日韩有码中文字幕| 欧美日韩视频精品一区| 亚洲精品在线美女| 女人精品久久久久毛片| 午夜免费成人在线视频| 桃花免费在线播放| av欧美777| 欧美中文综合在线视频| 天天躁狠狠躁夜夜躁狠狠躁|