• <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
    亚洲av一区综合| 夜夜爽天天搞| 国产蜜桃级精品一区二区三区| 熟女电影av网| 真实男女啪啪啪动态图| 国产人妻一区二区三区在| 国产精品一区二区三区四区久久| 国产一级毛片七仙女欲春2| 中文字幕久久专区| 国产中年淑女户外野战色| 国产精品自产拍在线观看55亚洲| 亚洲乱码一区二区免费版| 久久久久国产精品人妻aⅴ院| 亚洲第一电影网av| 99在线视频只有这里精品首页| 国产亚洲欧美在线一区二区| 此物有八面人人有两片| 欧美乱色亚洲激情| 中出人妻视频一区二区| 精品久久久久久久末码| www.999成人在线观看| 好男人电影高清在线观看| 欧美激情国产日韩精品一区| 亚洲男人的天堂狠狠| 成年人黄色毛片网站| 国产一区二区三区在线臀色熟女| 少妇被粗大猛烈的视频| 全区人妻精品视频| 午夜免费男女啪啪视频观看 | 在线a可以看的网站| 国产成人a区在线观看| xxxwww97欧美| 精品久久久久久久末码| 精品久久久久久久久av| 搡老岳熟女国产| 亚洲中文字幕日韩| 午夜日韩欧美国产| 亚洲久久久久久中文字幕| 麻豆国产av国片精品| 免费观看人在逋| 99热只有精品国产| 日本五十路高清| 中文资源天堂在线| 老司机福利观看| 欧美成狂野欧美在线观看| 我要看日韩黄色一级片| 国产免费av片在线观看野外av| 每晚都被弄得嗷嗷叫到高潮| 国产伦一二天堂av在线观看| 久久久久精品国产欧美久久久| 少妇人妻精品综合一区二区 | 婷婷色综合大香蕉| 久久精品91蜜桃| 91狼人影院| 久久国产精品人妻蜜桃| 日本撒尿小便嘘嘘汇集6| 成人无遮挡网站| 90打野战视频偷拍视频| 国产免费av片在线观看野外av| 久久精品人妻少妇| 欧美最新免费一区二区三区 | 亚洲第一区二区三区不卡| 亚洲成av人片在线播放无| 精品人妻一区二区三区麻豆 | 真人做人爱边吃奶动态| 亚洲av不卡在线观看| 国产精品嫩草影院av在线观看 | 神马国产精品三级电影在线观看| 亚洲色图av天堂| 免费一级毛片在线播放高清视频| 成人亚洲精品av一区二区| 亚洲人成电影免费在线| 亚洲精华国产精华精| 最近视频中文字幕2019在线8| avwww免费| 久久久久亚洲av毛片大全| 欧美日韩综合久久久久久 | 国产美女午夜福利| 色精品久久人妻99蜜桃| 99国产综合亚洲精品| 精品午夜福利在线看| 亚洲成人中文字幕在线播放| 两个人的视频大全免费| 国产白丝娇喘喷水9色精品| netflix在线观看网站| 观看免费一级毛片| 亚洲第一欧美日韩一区二区三区| 久久久国产成人免费| 亚洲精品影视一区二区三区av| 综合色av麻豆| 亚洲,欧美精品.| 中文字幕人成人乱码亚洲影| 成人毛片a级毛片在线播放| 变态另类成人亚洲欧美熟女| 国产亚洲精品av在线| 久久久久久久亚洲中文字幕 | av在线观看视频网站免费| 夜夜看夜夜爽夜夜摸| 黄色一级大片看看| 免费一级毛片在线播放高清视频| 国产在视频线在精品| 国产日本99.免费观看| 亚洲成av人片免费观看| 亚洲黑人精品在线| 欧美一区二区亚洲| 国产一区二区亚洲精品在线观看| 久久婷婷人人爽人人干人人爱| 午夜影院日韩av| 亚洲专区中文字幕在线| 亚洲 欧美 日韩 在线 免费| 国产毛片a区久久久久| 国产精品人妻久久久久久| 久久久精品大字幕| 亚洲国产精品久久男人天堂| 亚洲成a人片在线一区二区| 国产真实乱freesex| 99riav亚洲国产免费| 欧美在线一区亚洲| 深爱激情五月婷婷| 国产一区二区激情短视频| 我的女老师完整版在线观看| 最新中文字幕久久久久| 亚洲乱码一区二区免费版| 欧美+亚洲+日韩+国产| 国产精品一区二区性色av| 精品不卡国产一区二区三区| 成人国产综合亚洲| 久久国产精品影院| 一本综合久久免费| 在线国产一区二区在线| 亚洲三级黄色毛片| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | av黄色大香蕉| 色av中文字幕| 久久草成人影院| 亚洲欧美激情综合另类| bbb黄色大片| 国产欧美日韩一区二区精品| 日本一二三区视频观看| 午夜福利在线观看吧| 亚洲av熟女| 我要搜黄色片| av在线天堂中文字幕| 嫩草影院精品99| 欧美日韩中文字幕国产精品一区二区三区| 欧美精品国产亚洲| 黄片小视频在线播放| 香蕉av资源在线| 欧美午夜高清在线| 色综合亚洲欧美另类图片| 久久这里只有精品中国| 精品久久久久久久久亚洲 | 91午夜精品亚洲一区二区三区 | 一区二区三区激情视频| h日本视频在线播放| 伦理电影大哥的女人| 日韩精品中文字幕看吧| 在线十欧美十亚洲十日本专区| 国产午夜精品论理片| 在线观看舔阴道视频| 中文字幕av在线有码专区| 99久久精品一区二区三区| 亚洲欧美日韩高清在线视频| 国产高清三级在线| 99热这里只有是精品50| 日本一本二区三区精品| 亚洲人成网站高清观看| 日本一本二区三区精品| 精品久久久久久久久久免费视频| 九色国产91popny在线| 国产大屁股一区二区在线视频| 又爽又黄a免费视频| 亚洲第一电影网av| 国产大屁股一区二区在线视频| 观看免费一级毛片| 亚洲av成人精品一区久久| 日本一本二区三区精品| 国产色婷婷99| 成人欧美大片| 中文字幕高清在线视频| 99久久99久久久精品蜜桃| 搡老熟女国产l中国老女人| 欧美+日韩+精品| 欧美精品啪啪一区二区三区| 婷婷亚洲欧美| av天堂中文字幕网| 国产精品久久久久久久电影| 999久久久精品免费观看国产| 国产一区二区三区视频了| 国产精品一区二区性色av| 国产爱豆传媒在线观看| 神马国产精品三级电影在线观看| av在线观看视频网站免费| 少妇人妻一区二区三区视频| 欧美在线黄色| 热99在线观看视频| 他把我摸到了高潮在线观看| 欧美黄色片欧美黄色片| 长腿黑丝高跟| 90打野战视频偷拍视频| 在线观看美女被高潮喷水网站 | 亚洲久久久久久中文字幕| 级片在线观看| 欧美激情久久久久久爽电影| 国产伦精品一区二区三区视频9| 日韩欧美三级三区| 久久性视频一级片| 国产欧美日韩一区二区精品| 51国产日韩欧美| 国产精品亚洲av一区麻豆| 精品99又大又爽又粗少妇毛片 | 亚洲男人的天堂狠狠| 天堂影院成人在线观看| 国产色婷婷99| 国产精品不卡视频一区二区 | 欧美性猛交黑人性爽| 免费看美女性在线毛片视频| 51午夜福利影视在线观看| 国产精品久久视频播放| 男人和女人高潮做爰伦理| 无人区码免费观看不卡| 国产单亲对白刺激| 久久精品国产亚洲av天美| 亚洲人成网站高清观看| 中文字幕高清在线视频| 国产中年淑女户外野战色| 成年女人毛片免费观看观看9| 中文字幕免费在线视频6| 又黄又爽又刺激的免费视频.| 欧美最黄视频在线播放免费| 在线国产一区二区在线| 成人高潮视频无遮挡免费网站| 一级黄色大片毛片| 国产高清视频在线播放一区| 国产av在哪里看| 3wmmmm亚洲av在线观看| 一进一出好大好爽视频| 女生性感内裤真人,穿戴方法视频| 亚洲av一区综合| 一边摸一边抽搐一进一小说| 亚洲欧美日韩高清专用| 两人在一起打扑克的视频| 亚洲人成网站在线播| 国内精品久久久久精免费| 久久久久免费精品人妻一区二区| 亚洲经典国产精华液单 | 国产精品久久久久久人妻精品电影| 亚洲五月婷婷丁香| 少妇人妻精品综合一区二区 | 夜夜爽天天搞| 国产真实乱freesex| 如何舔出高潮| 高潮久久久久久久久久久不卡| 亚洲成av人片在线播放无| 亚洲 国产 在线| 久久这里只有精品中国| 国产成人福利小说| 亚洲欧美激情综合另类| 欧美激情国产日韩精品一区| 少妇被粗大猛烈的视频| 成年人黄色毛片网站| 丁香六月欧美| 3wmmmm亚洲av在线观看| 琪琪午夜伦伦电影理论片6080| 欧美在线一区亚洲| 日韩欧美在线二视频| 日日夜夜操网爽| 亚洲专区中文字幕在线| 少妇裸体淫交视频免费看高清| 91麻豆av在线| 国产精品综合久久久久久久免费| 综合色av麻豆| 欧美3d第一页| 亚洲成人免费电影在线观看| 国产一区二区三区在线臀色熟女| 久久草成人影院| 国语自产精品视频在线第100页| 最新中文字幕久久久久| 午夜老司机福利剧场| 搡老熟女国产l中国老女人| 亚洲中文字幕一区二区三区有码在线看| 十八禁网站免费在线| 亚洲人与动物交配视频| 一本综合久久免费| 欧美+亚洲+日韩+国产| 最后的刺客免费高清国语| 美女被艹到高潮喷水动态| av在线观看视频网站免费| 最近最新免费中文字幕在线| 亚洲av成人精品一区久久| 1000部很黄的大片| 夜夜爽天天搞| 看免费av毛片| 村上凉子中文字幕在线| 男人的好看免费观看在线视频| 国产精品野战在线观看| 亚洲真实伦在线观看| 国产精品精品国产色婷婷| 99久久成人亚洲精品观看| 搡老岳熟女国产| 99久久精品国产亚洲精品| 久久久久久九九精品二区国产| 亚洲欧美清纯卡通| 看十八女毛片水多多多| 久久亚洲真实| 国产精品98久久久久久宅男小说| 日本 欧美在线| 亚洲,欧美,日韩| 免费av毛片视频| 国内精品久久久久久久电影| 国产精品野战在线观看| 久久精品国产清高在天天线| 性插视频无遮挡在线免费观看| 国产 一区 欧美 日韩| 每晚都被弄得嗷嗷叫到高潮| 精品免费久久久久久久清纯| 亚洲人成网站高清观看| av国产免费在线观看| 国产人妻一区二区三区在| 欧美精品啪啪一区二区三区| 哪里可以看免费的av片| 色综合欧美亚洲国产小说| 九九在线视频观看精品| 一级毛片久久久久久久久女| 欧美国产日韩亚洲一区| 国产免费一级a男人的天堂| 亚州av有码| 人妻久久中文字幕网| av在线蜜桃| 欧美性猛交╳xxx乱大交人| 欧美国产日韩亚洲一区| 91字幕亚洲| 国产黄色小视频在线观看| 高清在线国产一区| 嫩草影视91久久| 亚洲av一区综合| 一本精品99久久精品77| 村上凉子中文字幕在线| 免费观看人在逋| 在线国产一区二区在线| 九九久久精品国产亚洲av麻豆| 一级作爱视频免费观看| 一个人看视频在线观看www免费| 97超视频在线观看视频| 少妇人妻精品综合一区二区 | 中文字幕高清在线视频| 久久婷婷人人爽人人干人人爱| 麻豆国产av国片精品| 看十八女毛片水多多多| 一级毛片久久久久久久久女| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 国产精品久久久久久精品电影| 久久久久性生活片| 国产精品一区二区性色av| 久久精品国产亚洲av涩爱 | 99国产精品一区二区三区| 亚洲无线在线观看| 一区二区三区免费毛片| 成人一区二区视频在线观看| avwww免费| 日韩有码中文字幕| 国产白丝娇喘喷水9色精品| 哪里可以看免费的av片| 欧洲精品卡2卡3卡4卡5卡区| 91字幕亚洲| 一级作爱视频免费观看| 国产高清三级在线| 久久人人精品亚洲av| 日韩欧美精品免费久久 | 美女xxoo啪啪120秒动态图 | 国产精品一及| 日韩欧美在线二视频| 精品一区二区三区视频在线观看免费| 丰满乱子伦码专区| 久久精品国产亚洲av天美| 丁香欧美五月| 别揉我奶头~嗯~啊~动态视频| 欧美乱色亚洲激情| 久久精品国产清高在天天线| 亚洲欧美激情综合另类| 日韩中文字幕欧美一区二区| 一夜夜www| 老司机深夜福利视频在线观看| 国内毛片毛片毛片毛片毛片| 黄色一级大片看看| www.www免费av| 亚洲成人免费电影在线观看| 麻豆成人午夜福利视频| 久久久久亚洲av毛片大全| 欧美色欧美亚洲另类二区| 亚洲乱码一区二区免费版| 欧美中文日本在线观看视频| 一a级毛片在线观看| 国产精品久久电影中文字幕| 色吧在线观看| 可以在线观看毛片的网站| 国产亚洲欧美在线一区二区| 色5月婷婷丁香| 亚洲av第一区精品v没综合| 欧美激情久久久久久爽电影| 最后的刺客免费高清国语| 日本黄色视频三级网站网址| 成人性生交大片免费视频hd| 国产精品久久久久久亚洲av鲁大| 97超级碰碰碰精品色视频在线观看| 亚洲自拍偷在线| 成人亚洲精品av一区二区| 久久精品国产亚洲av香蕉五月| 中文字幕人妻熟人妻熟丝袜美| 99热只有精品国产| 亚洲av中文字字幕乱码综合| 少妇裸体淫交视频免费看高清| 欧美精品啪啪一区二区三区| 日韩高清综合在线| 成人av一区二区三区在线看| 国产视频内射| 久久精品夜夜夜夜夜久久蜜豆| 日本三级黄在线观看| 色综合亚洲欧美另类图片| 又黄又爽又刺激的免费视频.| 国产国拍精品亚洲av在线观看| 赤兔流量卡办理| 国产精品女同一区二区软件 | avwww免费| 午夜久久久久精精品| 波多野结衣高清作品| 天堂av国产一区二区熟女人妻| 听说在线观看完整版免费高清| 免费看光身美女| 亚洲精华国产精华精| 久久精品91蜜桃| 欧美日韩福利视频一区二区| 无遮挡黄片免费观看| 亚洲久久久久久中文字幕| 宅男免费午夜| 成人精品一区二区免费| 狠狠狠狠99中文字幕| av福利片在线观看| 简卡轻食公司| 少妇人妻一区二区三区视频| 国产av一区在线观看免费| 亚洲自偷自拍三级| netflix在线观看网站| 一区二区三区四区激情视频 | 欧美色欧美亚洲另类二区| 亚洲精品成人久久久久久| 国产精品久久电影中文字幕| 国产高清视频在线观看网站| 看黄色毛片网站| а√天堂www在线а√下载| 亚洲精品一卡2卡三卡4卡5卡| 午夜福利视频1000在线观看| 亚洲 国产 在线| 亚洲第一欧美日韩一区二区三区| 日韩欧美一区二区三区在线观看| 亚洲一区二区三区不卡视频| 免费搜索国产男女视频| av在线蜜桃| 97超级碰碰碰精品色视频在线观看| 国产一区二区在线观看日韩| 国产综合懂色| 色综合站精品国产| 国产精品一区二区三区四区免费观看 | 欧美一区二区精品小视频在线| 桃色一区二区三区在线观看| 一区福利在线观看| 亚洲在线自拍视频| 女生性感内裤真人,穿戴方法视频| 亚洲欧美日韩高清在线视频| 亚洲18禁久久av| 国产色婷婷99| 国产欧美日韩精品亚洲av| 夜夜躁狠狠躁天天躁| a级毛片a级免费在线| 亚洲欧美日韩东京热| 18禁裸乳无遮挡免费网站照片| 在线观看av片永久免费下载| 午夜福利成人在线免费观看| 亚洲18禁久久av| 91九色精品人成在线观看| 欧美成人a在线观看| av在线蜜桃| 美女高潮的动态| 色综合亚洲欧美另类图片| 成人欧美大片| 亚洲欧美日韩卡通动漫| 色哟哟哟哟哟哟| 最新在线观看一区二区三区| 国产91精品成人一区二区三区| 亚洲av中文字字幕乱码综合| 欧美在线一区亚洲| 一区二区三区高清视频在线| 特大巨黑吊av在线直播| 亚洲国产欧洲综合997久久,| 国产免费男女视频| 久久性视频一级片| 午夜影院日韩av| 国产在视频线在精品| 丰满乱子伦码专区| 亚洲av二区三区四区| 成人av在线播放网站| 亚洲成人中文字幕在线播放| 亚洲在线观看片| 亚洲av一区综合| 亚洲成人免费电影在线观看| xxxwww97欧美| 亚洲综合色惰| 国产精品不卡视频一区二区 | 91午夜精品亚洲一区二区三区 | .国产精品久久| 一进一出抽搐动态| 久久这里只有精品中国| 国产激情偷乱视频一区二区| 成人午夜高清在线视频| 国产高潮美女av| 免费看日本二区| 亚洲自偷自拍三级| 伊人久久精品亚洲午夜| 国产欧美日韩一区二区精品| 日韩中字成人| 中文字幕久久专区| 欧美zozozo另类| 国产精品综合久久久久久久免费| 亚洲avbb在线观看| 久久精品综合一区二区三区| 亚洲av二区三区四区| 国产精品美女特级片免费视频播放器| 在线观看一区二区三区| 成人欧美大片| 欧美黑人巨大hd| 成人国产一区最新在线观看| 亚洲第一电影网av| 听说在线观看完整版免费高清| 女同久久另类99精品国产91| 色视频www国产| 欧美黄色片欧美黄色片| 国产探花极品一区二区| 麻豆国产av国片精品| eeuss影院久久| 国产 一区 欧美 日韩| 国产av不卡久久| 色综合婷婷激情| 久久这里只有精品中国| 日本五十路高清| 国产精品亚洲av一区麻豆| 91av网一区二区| 亚洲精品456在线播放app | 中文字幕av成人在线电影| 精品一区二区三区视频在线| 日本五十路高清| 婷婷丁香在线五月| 香蕉av资源在线| 国内精品久久久久精免费| 天堂√8在线中文| 免费av观看视频| 一级毛片久久久久久久久女| 欧美日韩瑟瑟在线播放| 欧美成人一区二区免费高清观看| 亚洲国产高清在线一区二区三| 久久久色成人| 国产精品综合久久久久久久免费| 国产成人aa在线观看| 欧美日本亚洲视频在线播放| 久久久国产成人精品二区| 在线观看一区二区三区| 国产大屁股一区二区在线视频| 少妇的逼好多水| 国产熟女xx| 免费人成视频x8x8入口观看| 亚洲av二区三区四区| 给我免费播放毛片高清在线观看| 12—13女人毛片做爰片一| 99热这里只有是精品在线观看 | 一级av片app| 国产精品久久久久久精品电影| 精品国内亚洲2022精品成人| 九色成人免费人妻av| 色尼玛亚洲综合影院| 99久久99久久久精品蜜桃| 成人特级av手机在线观看| 欧美日本视频| 午夜亚洲福利在线播放| 亚洲最大成人中文| 性插视频无遮挡在线免费观看| 夜夜看夜夜爽夜夜摸| 狂野欧美白嫩少妇大欣赏| 美女大奶头视频| 亚洲第一电影网av| 国产成人福利小说| 国产精品99久久久久久久久| 18美女黄网站色大片免费观看| 精品一区二区三区视频在线观看免费| 久99久视频精品免费| 久久久成人免费电影| 亚洲无线观看免费| 麻豆国产av国片精品| www.色视频.com| 国产探花极品一区二区| 亚洲人与动物交配视频| 亚洲av美国av| h日本视频在线播放| 亚洲一区高清亚洲精品| 国产av不卡久久| 欧美日韩综合久久久久久 | 性色avwww在线观看| www日本黄色视频网| 免费av观看视频| 亚洲欧美日韩高清专用| 男女做爰动态图高潮gif福利片| 91在线观看av| 亚洲国产色片| 国产精品一及| 伦理电影大哥的女人| 最新在线观看一区二区三区| 亚洲精品粉嫩美女一区|