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

    非共沸混合工質(zhì)噴射器內(nèi)部流動特性

    2020-05-22 09:28:08楊孟柯黃鯉生

    楊孟柯, 王 林, 黃鯉生

    (河南科技大學 土木工程學院,河南 洛陽 471023)

    0 引言

    噴射制冷機具有體積小、機械能消耗少和低溫熱源驅(qū)動等優(yōu)點,但由于噴射制冷系統(tǒng)壓縮比較低,造成制冷溫度較高,制冷效率較低,成為制約噴射制冷技術(shù)推廣的關(guān)鍵因素。低壓縮比不僅要求較高蒸發(fā)溫度和較低冷凝溫度的運行條件,還限制高壓縮比制冷劑應(yīng)用于噴射制冷。噴射器壓縮比增加勢必會引起噴射因數(shù)下降,進而降低系統(tǒng)性能系數(shù),可見,壓縮比與噴射因數(shù)對噴射制冷循環(huán)的影響是相互制約、此消彼長的關(guān)系。

    為了克服噴射制冷系統(tǒng)低壓縮比的局限性,文獻[1-2]嘗試將非共沸混合制冷劑應(yīng)用于噴射制冷系統(tǒng)中,結(jié)果表明: 選擇合適的混合工質(zhì),可以使噴射器在較低壓縮比下制取較低溫度的冷量。文獻[3]采用R134a/R143a非共沸混合工質(zhì)對噴射膨脹制冷系統(tǒng)進行了研究并得出:當混合工質(zhì)中R134a的質(zhì)量分數(shù)為0.9時,該系統(tǒng)的制冷性能系數(shù)(coefficient of performance,COP)比純工質(zhì)時最大提高了10.47%。文獻[4]研究了R290/R600a非共沸混合制冷劑的噴射器強化壓縮制冷循環(huán),利用混合物組分變化特性,有效降低了噴射器壓縮比,進而提高了系統(tǒng)性能。文獻[5]分析了純R245fa制冷劑、純R601a制冷劑及R245fa/R601a混合制冷劑對串聯(lián)式冷電聯(lián)產(chǎn)噴射制冷系統(tǒng)的熱力學性能影響,并通過多目標優(yōu)化找出最佳組分配比。為了能夠在相同冷凝溫度下獲取更低的制冷溫度,文獻[6]在傳統(tǒng)噴射制冷循環(huán)中增設(shè)了一個噴射器和發(fā)生器,設(shè)計出一種雙溫熱源雙噴射制冷系統(tǒng),從而提升了系統(tǒng)的整體壓縮比,在同等條件下,該系統(tǒng)的蒸發(fā)溫度比傳統(tǒng)噴射制冷系統(tǒng)低10~15 ℃。文獻[7-9]將非共沸混合工質(zhì)依次應(yīng)用于新提出的一級分凝和二級分凝的自復(fù)疊噴射制冷系統(tǒng)中,可以在噴射器較低壓縮比下制取-40 ℃的冷量。

    噴射器作為噴射制冷系統(tǒng)的關(guān)鍵部件,其性能的優(yōu)劣直接影響整個系統(tǒng)的運行狀況。但由于噴射器自身獨特的結(jié)構(gòu)特征,使噴射器內(nèi)部的流動特性比較復(fù)雜,如跨音速流動、壅塞流動、激波和混合邊界層等。文獻[10]研究了亞音速和超音速噴射器的流動特性。文獻[11]提出了一種預(yù)測噴射器性能的臨界圓模型,該模型比許多一維建模方法更簡單,并能更準確地預(yù)測臨界工作模式噴射器的性能。文獻[12]提出用復(fù)合壅塞理論解釋超音速噴射器中流體的壅塞流動現(xiàn)象。文獻[13]研究了激波結(jié)構(gòu)對引射性能和臨界背壓的影響。文獻[14]研究了邊界層分離的原因、分離區(qū)域的確定以及邊界層分離對噴射器性能的影響。文獻[15]研究了混合邊界層的生長對噴射器性能的影響。

    只有對噴射器內(nèi)部的流動特性有了深刻的了解,才能更快地找到提升噴射器引射性能的方法,進而提升整個制冷系統(tǒng)的性能。目前,已有的研究工作集中于混合工質(zhì)在制冷循環(huán)系統(tǒng)中的系統(tǒng)特性,而對混合工質(zhì)噴射器內(nèi)部流體的工作特性研究尚未成熟。本文在對純工質(zhì)噴射器現(xiàn)有研究的基礎(chǔ)上,選擇非共沸混合工質(zhì)R32/R245fa作為噴射器內(nèi)部流體,對非共沸混合工質(zhì)噴射器的內(nèi)部流動特性進行初步研究,捕捉非共沸混合工質(zhì)噴射器內(nèi)部的激波現(xiàn)象,探尋非共沸混合工質(zhì)噴射器內(nèi)部流體的速度場、壓力場、溫度場及濃度場的變化規(guī)律,分析非共沸混合工質(zhì)噴射器內(nèi)部流體的傳熱傳質(zhì)過程,可為后期混合工質(zhì)噴射器結(jié)構(gòu)尺寸及運行工況優(yōu)化提供一定的理論指導(dǎo)。

    1 混合工質(zhì)噴射器尺寸設(shè)計計算

    圖1 噴射器結(jié)構(gòu)圖

    為了研究混合工質(zhì)噴射器流動特性,首先,必須要設(shè)計噴射器的結(jié)構(gòu)尺寸,本文運用氣體動力學函數(shù)法對噴射器結(jié)構(gòu)進行設(shè)計計算。氣體動力學函數(shù)法是基于氣體動力學理論,引入折算等熵速度建立相對溫度、相對壓力及相對比容,根據(jù)守恒定律及自由流束理論,計算噴射因數(shù)及相關(guān)尺寸,該計算方法可靠性較高。目前,大多數(shù)研究者運用該方法計算純工質(zhì)噴射器結(jié)構(gòu)尺寸,在公開文獻中尚未見到將氣體動力學函數(shù)法應(yīng)用于混合工質(zhì)噴射器結(jié)構(gòu)的設(shè)計計算。下面將基于索科洛夫氣體動力學函數(shù)法建立混合工質(zhì)噴射器動力學模型,計算噴射因數(shù),并設(shè)計混合工質(zhì)噴射器各個部位的結(jié)構(gòu)尺寸。噴射器結(jié)構(gòu)如圖1所示。

    將折算等熵速度應(yīng)用于噴射器中混合室的動量守恒方程,再結(jié)合質(zhì)量守恒方程,推導(dǎo)出噴射器噴射因數(shù)[16]:

    (1)

    其中:Κ為速度系數(shù);ag*、ae*和ac*分別為工作流體、引射流體及出口流體的臨界速度;λge、λe2和λc3分別為2-2截面上的工作流體、2-2截面上的引射流體及3-3截面上混合流體的折算等熵速度。

    K1=φ1φ2φ3;

    (2)

    K2=φ2φ3φ4;

    (3)

    (4)

    (5)

    其中:φ1、φ2、φ3和φ4分別為噴嘴、混合室、擴壓室及混合室入口的速度系數(shù);β為混合室漸縮段進出口面積比;Pg、Pe和Pc分別為工作流體、引射流體及出口流體的壓力;Πc2、Πe2和Πc3分別為2-2截面上的混合流體、2-2截面上的引射流體及3-3截面上的混合流體的相對壓力;Πg*和Πe*分別為工作流體與引射流體的臨界相對壓力;kg和ke分別為工作流體與引射流體的絕熱指數(shù);qge和qe2分別為2-2截面上工作流體與引射流體的折算質(zhì)量流量。

    在混合室某一位置(假設(shè)圖1中s-s截面),引射流體的速度達到臨界流速,即噴射器中出現(xiàn)第二臨界狀態(tài),此時相對應(yīng)的噴射因數(shù)為:

    (6)

    其中:μng為噴射器中出現(xiàn)第二臨界狀態(tài)時的噴射因數(shù);kc為混合流體的絕熱指數(shù);Πc*為混合流體的臨界相對壓力;qgs和qc3分別為s-s截面上工作流體的折算質(zhì)量流量和3-3截面上混合流體的折算質(zhì)量流量。

    根據(jù)上述計算公式并添加組分守恒定律方程,對R32/R245fa非共沸混合工質(zhì)噴射器進行幾何尺寸計算。計算典型工況為:工作流體,溫度363 K,壓力1 153 kPa,工作流體中R32的質(zhì)量分數(shù)為0.05;引射流體,溫度283 K,壓力134 kPa,引射流體中R32的質(zhì)量分數(shù)為0.20;混合流體,噴射器出口溫度309 K,壓力273 kPa。在MATLAB軟件編程環(huán)境中調(diào)用REFPROP中R32/R245fa混合工質(zhì)物性程序,編寫計算程序,得出最大噴射因數(shù)為0.33。根據(jù)所得的最大噴射因數(shù)確定R32/R245fa非共沸混合工質(zhì)噴射器各部位的尺寸,結(jié)果見表1,以此尺寸為基礎(chǔ),對混合工質(zhì)噴射器進行數(shù)值模擬,以研究噴射器內(nèi)部流態(tài)。

    表1 R32/R245fa非共沸混合工質(zhì)噴射器結(jié)構(gòu)尺寸 mm

    2 模擬計算

    相比于噴嘴出口的超音速流動,工作流體和引射流體入口速度較小,可以忽略不計,因此可以將引射流體入口簡化為軸向環(huán)形入口[17],與工作流體同側(cè)。又由于此時噴射器關(guān)于軸心完全對稱,文獻[18]對噴射器二維對稱軸模型和三維模型進行了比較,比較結(jié)果基本一致,可再次將其簡化為以軸心為橫坐標的二維模型。在保證計算精度的前提下,為更一步縮減計算時間,對二維軸對稱模型計算一半即可。

    根據(jù)噴射器幾何結(jié)構(gòu)尺寸,利用ICEM軟件建立二維噴射器模型并進行結(jié)構(gòu)性網(wǎng)格的劃分,而網(wǎng)格的質(zhì)量直接影響數(shù)值模擬計算的精確度和所需時間,因此在實際模擬計算前,需對網(wǎng)格無關(guān)性進行檢驗以確保模擬結(jié)果具有準確性。將噴嘴出口中心及擴壓室入口中心作為觀測點,設(shè)置不同的網(wǎng)格數(shù)量,觀察這兩個位置壓力和速度的變化,結(jié)果如圖2所示。圖2a中網(wǎng)格數(shù)在4×104以內(nèi)時,兩個觀測點的參數(shù)均變化較大,噴嘴出口壓力隨網(wǎng)格數(shù)的增加而降低,噴嘴出口速度隨網(wǎng)格數(shù)的增加而增加。圖2b中擴壓室入口壓力隨網(wǎng)格數(shù)的增加而增加,擴壓室入口速度隨網(wǎng)格數(shù)的增加而降低。當網(wǎng)格數(shù)達到4×104以后,噴嘴出口及擴壓室入口的壓力與速度基本維持不變。因此,在數(shù)值模擬研究中,網(wǎng)格數(shù)設(shè)置在4×104左右,網(wǎng)格質(zhì)量保持在0.9以上。同時,為了更準確地捕捉內(nèi)部流體的激波現(xiàn)象,對流動劇烈的地方進行適當?shù)木W(wǎng)格加密處理,噴射器壁面也進行了逐層加密。

    (a) 噴嘴出口中心

    (b) 擴壓室入口中心

    圖2 兩個觀測點的壓力與速度隨網(wǎng)格數(shù)的變化圖

    本文利用Fluent 19.0軟件進行求解計算,選用基于壓力的二維軸對稱求解器,計算過程中要遵守質(zhì)量、動量、能量及組分守恒控制方程。湍流模型采用標準k-epsilion模型,近壁處選用標準壁面函數(shù)法。對噴射器入口采用壓力入口邊界條件,出口采用壓力出口邊界條件,壁面采用絕熱無滑移壁面邊界條件?;旌瞎べ|(zhì)噴射器邊界條件見表2。Fluent 19.0軟件采用有限體積法對控制方程進行離散,求解方法中壓力-速度耦合方法選用Simple算法,梯度項選用基于最小二乘單元,壓力項選用二階格式,其余項均采用二階迎風格式?;旌瞎べ|(zhì)物性通過調(diào)用NIST multispecies real gas model獲得,該材料模型只能限制在單相狀態(tài)下使用,為了保證噴射器內(nèi)部流體全為氣態(tài),進出口溫度均設(shè)5 K的過熱度[19]。

    表2 混合工質(zhì)噴射器邊界條件

    為了驗證模型的準確性與可靠性,以文獻[20]的試驗數(shù)據(jù)作為參考值,對該文獻所提供的噴射器結(jié)構(gòu)尺寸建立模型并進行模擬計算。相同結(jié)構(gòu)尺寸與工況下,模擬計算結(jié)果與試驗結(jié)果如表3所示。由表3可知:模擬計算結(jié)果與試驗結(jié)果誤差均在5%以內(nèi),在允許的誤差范圍以內(nèi)。

    表3 模擬計算結(jié)果與試驗結(jié)果對照表

    3 結(jié)果分析

    圖3為混合工質(zhì)(R32/R245fa)噴射器內(nèi)部速度云圖。由圖3可以觀察到:較低速度的工作流體由拉伐爾噴嘴進入噴射器,在噴嘴漸縮段速度迅速增加,在噴嘴出口出現(xiàn)了鉆石激波鏈,激波核心沿噴嘴下游方向逐漸衰減,消失在混合室等截面前段,該激波鏈的出現(xiàn)會引起能量損失。到擴壓室入口時,充分混合的流體速度突然降低,出現(xiàn)正激波,流體突然由超音速流動降至亞音速流動,動力能轉(zhuǎn)化為壓力能,流體速度在擴壓室內(nèi)繼續(xù)降低,引起壓力繼續(xù)增加,最終噴射器出口的壓力大于出口背壓,流體順利從噴射器出口流出。圖4為混合工質(zhì)噴射器的馬赫數(shù)沿噴射器軸向分布圖。由圖4可知:在噴嘴喉部,馬赫數(shù)增加到1.0,工作流體由亞音速增至音速,進入噴嘴漸擴段速度繼續(xù)增加;在噴嘴出口,馬赫數(shù)增加到2.0,壓力能轉(zhuǎn)化為動能;在噴嘴出口下游,馬赫數(shù)繼續(xù)增加,速度持續(xù)攀升,馬赫數(shù)最大達到2.5,工作流體在噴嘴出口迅速膨脹形成高速射流以卷吸低速的引射流體。兩股流體在混合過程中,馬赫數(shù)從噴嘴出口至混合室等截面前段,上下波動幅度沿噴射器軸向逐漸減小。在整個過程中,噴嘴出口的高速射流不斷膨脹壓縮形成了連續(xù)激波,而噴嘴出口的工作流體與引射流體速度差較大,在工作流體與引射流體接觸面形成了速度邊界層,高速的工作流體與低速引射流體在混合邊界層產(chǎn)生劇烈碰撞,兩股流體不斷進行動量交換,工作流體將自身部分動量轉(zhuǎn)移給引射流體,兩股流體相互摻混。隨著流體向混合室涌入,速度邊界層逐漸加厚,兩股流體混合越來越均勻,徑向速度梯度逐漸縮小,在混合室等截面段達到充分混合。從圖4中還可以看出:整個混合室流體的馬赫數(shù)始終大于1.3,表明流體一直處于超音速流動,到擴壓室出口馬赫數(shù)降低至0.2,因此,能量損失較小。由計算得知:與傳統(tǒng)的純工質(zhì)噴射器相比,混合工質(zhì)噴射器具有更好的性能。相同進出口溫度下,混合工質(zhì)(R32/R245fa)進入噴射器的工作流體質(zhì)量流量為0.196 kg/s,卷吸引射流體的質(zhì)量流量為0.067 kg/s,噴射因數(shù)達到了0.34;而采用純工質(zhì)(R245fa)的噴射器噴射因數(shù)為0.19。混合工質(zhì)(R32/R245fa)噴射器的噴射因數(shù)較純工質(zhì)(R245fa)噴射器的噴射因數(shù)提高了78.95%,與采用氣體動力學函數(shù)法所獲得的混合工質(zhì)噴射器的噴射系統(tǒng)相比,采用模擬方法所獲得的混合工質(zhì)(R32/R245fa)噴射器的噴射因數(shù)也提高了3.03%。

    圖3 混合工質(zhì)噴射器內(nèi)部速度云圖

    圖4 混合工質(zhì)噴射器的馬赫數(shù)沿噴射器軸向分布圖

    圖5為混合工質(zhì)(R32/R245fa)噴射器內(nèi)部壓力云圖。由圖5可以看出:壓力從噴嘴入口到噴嘴出口變化巨大,并且在噴嘴出口下游出現(xiàn)一條明顯的鉆石激波鏈。圖6為混合工質(zhì)噴射器的靜壓沿噴射器軸向分布圖。由圖6可以看出:高壓狀態(tài)的工作流體經(jīng)過噴嘴后壓力驟降,工作流體經(jīng)過噴嘴后壓力由1 153 kPa降低至117 kPa,在噴嘴出口處壓力繼續(xù)降低形成真空域。計算結(jié)果顯示噴射器中工作流體噴嘴出口下游與引射流體入口的最高壓差為78 kPa,在這個壓差作用下引射流體被抽吸進入吸入室。超音速的工作流體因在噴嘴漸擴段出口流道突然擴大,在噴嘴出口出現(xiàn)斜激波,斜激波到混合邊界層被反射為膨脹波,膨脹波到混合邊界層被反射為壓縮波,經(jīng)過多次反射而產(chǎn)生了鉆石激波鏈,激波的能量不斷耗散,直至激波消失殆盡。從圖6還可以明顯看出:噴嘴出口壓力出現(xiàn)了一段連續(xù)波動。工作流體與引射流體在混合室等截面段以等壓的形式進行混合,并繼續(xù)向前流動,充分混合的流體進入擴壓室時,因壓力先小幅下降再突然升高而出現(xiàn)正激波,該現(xiàn)象有助于噴射器內(nèi)部流體克服出口背壓,混合后流體進入擴壓室壓力進一步提升,直至達到噴射器出口背壓。

    圖7為混合工質(zhì)(R32/R245fa)噴射器內(nèi)部溫度云圖。由圖7可以看出:流體的溫度經(jīng)過噴嘴后大幅降低,并且在噴嘴出口下游形成了兩個較大的低溫區(qū)。圖8為混合工質(zhì)噴射器的溫度沿噴射器軸向分布圖。由圖8可以看出:混合工質(zhì)高溫高壓的工作流體經(jīng)過拉伐爾噴嘴,溫度急劇下滑,對應(yīng)速度急速上升,工作流體的內(nèi)能轉(zhuǎn)化為動能。噴嘴出口工作流體溫度繼續(xù)下降形成低溫區(qū),引射流體流入吸入室,在激波的作用下,工作流體與引射流體在溫度邊界層內(nèi)開始混合,兩股流體能量與動量相互傳遞,溫度也隨著波動升高?;旌鲜覞u縮段,由于激波強度逐漸減弱,兩股流體溫度梯度逐漸縮小;混合室等截面段,隨著激波的消失,兩股流體混合充分,溫度達到一致。從圖8中還可以看出:在混合室等截面段噴射器內(nèi)部流體溫度維持在322 K幾乎不變。混合后的流體在擴壓室入口因正激波作用,溫度出現(xiàn)先小幅降低后大幅提升的現(xiàn)象。進入擴壓室后,速度逐漸降低,流體的動能轉(zhuǎn)化為內(nèi)能,因此,溫度繼續(xù)緩慢提升。

    圖5 混合工質(zhì)噴射器內(nèi)部壓力云圖

    圖6 混合工質(zhì)噴射器的靜壓沿噴射器軸向分布圖

    圖7 混合工質(zhì)噴射器內(nèi)部溫度云圖

    圖8 混合工質(zhì)噴射器的溫度沿噴射器軸向分布圖

    圖9和圖10呈現(xiàn)了混合工質(zhì)(R32/R245fa)噴射器內(nèi)部R32的質(zhì)量分數(shù)分布云圖及噴射器內(nèi)流體R32質(zhì)量分數(shù)的變化。從圖9中可以看出:噴嘴出口處工作流體中R32的質(zhì)量分數(shù)與入口相同,為0.05;而引射流體中R32的質(zhì)量分數(shù)為0.20。噴嘴出口的高速工作流體卷吸引射流體后,因兩股流體存在質(zhì)量分數(shù)差,在工作流體與引射流體接觸面上,出現(xiàn)濃度邊界層,工作流體與引射流體會在此發(fā)生組分遷移,引射流體中較高質(zhì)量分數(shù)的R32向工作流體遷移,同樣,工作流體中較高質(zhì)量分數(shù)的R245fa向引射流體遷移。在混合室入口之前,激波核心內(nèi)工作流體R32的質(zhì)量分數(shù)不變。由圖10可知:自噴射器混合室入口開始,工作流體與引射流體中R32的質(zhì)量分數(shù)出現(xiàn)劇烈變化,兩股流體迅速發(fā)生組分傳遞。混合室漸縮段濃度邊界層不斷加厚,兩股流體進一步混合,濃度梯度逐漸縮小,混合室等截面段隨著激波的消散,濃度邊界層將被破壞,更加速了工作流體與引射流體的傳質(zhì)過程,質(zhì)量分數(shù)逐漸趨于一致。從圖9中可以看出:在擴壓室入口R32已混合充分,組分傳遞過程完成。擴壓室內(nèi)混合工質(zhì)(R32/R245fa)以穩(wěn)定的質(zhì)量分數(shù)流出噴射器。由于混合工質(zhì)中只有兩種組分,這兩種組分需相互平衡,R32與R245fa的質(zhì)量分數(shù)分布必是此消彼長。因此,只需分析一種組分變化即可。

    圖9 混合工質(zhì)噴射器內(nèi)部R32的質(zhì)量分數(shù)分布云圖

    圖10 混合工質(zhì)噴射器內(nèi)R32質(zhì)量分數(shù)的變化圖

    綜上所述,噴嘴出口出現(xiàn)了一條鉆石激波鏈,該斜激波鏈雖有助于工作流體與引射流體的混合,但也消耗了流體的部分能量,導(dǎo)致流體混合后能量降低。當擴壓室入口出現(xiàn)正激波,該正激波的出現(xiàn)有助于噴射器內(nèi)部流體克服出口背壓,幫助流體順利從噴射器出口流出,因此,混合工質(zhì)噴射器在雙激波的狀態(tài)下工作性能較佳。工作流體與引射流體在混合室等截面段近似以等溫等壓形式混合,減少了混合過程的能量損失。對于混合工質(zhì)(R32/R245fa)自復(fù)疊噴射制冷系統(tǒng)而言,混合工質(zhì)(R32/R245fa)噴射器的噴射因數(shù)增加,一方面會使制冷效率提升,另一方面也會使制冷系統(tǒng)獲得更低的制冷溫度。

    4 結(jié)論

    (1)混合工質(zhì)(R32/R245fa)噴射器性能優(yōu)于純工質(zhì)(R245fa)噴射器性能。在噴射器壓縮比不高時,混合工質(zhì)噴射器的能量損失較小,噴射因數(shù)高于純工質(zhì)噴射器。在給定工況下,噴射器噴嘴出口下游馬赫數(shù)高達2.5,在整個混合室內(nèi)部馬赫數(shù)始終大于1.3,流體一直處于超音速流動,在擴壓室出口馬赫數(shù)降低至0.2。混合工質(zhì)(R32/R245fa)噴射器的噴射因數(shù)較純工質(zhì)(R245fa)噴射器提高了78.95%。

    (2)混合工質(zhì)(R32/R245fa)噴射器的組分遷移具有顯著特征。工作流體和引射流體的接觸面上出現(xiàn)濃度邊界層。在混合室漸縮段前段,激波核心內(nèi)的工作流體R32的質(zhì)量分數(shù)一直不變,但在進入混合室以后,工作流體與引射流體中R32和R245fa在濃度邊界層迅速發(fā)生組分傳遞。

    热re99久久国产66热| 欧美乱码精品一区二区三区| 久久香蕉精品热| 新久久久久国产一级毛片| 免费av毛片视频| 亚洲成av片中文字幕在线观看| 国产av精品麻豆| 日韩大码丰满熟妇| 午夜久久久在线观看| 看片在线看免费视频| 手机成人av网站| 久久久久久免费高清国产稀缺| 亚洲第一av免费看| 俄罗斯特黄特色一大片| 热99re8久久精品国产| 久久午夜综合久久蜜桃| 日韩精品中文字幕看吧| 国产深夜福利视频在线观看| 国产av一区二区精品久久| 91成年电影在线观看| av超薄肉色丝袜交足视频| 操出白浆在线播放| 久久人妻av系列| 亚洲欧洲精品一区二区精品久久久| 国产亚洲欧美精品永久| 成人三级做爰电影| 日韩欧美在线二视频| 丰满迷人的少妇在线观看| xxx96com| 精品日产1卡2卡| 午夜老司机福利片| 成年人黄色毛片网站| 亚洲精华国产精华精| 97超级碰碰碰精品色视频在线观看| 搡老乐熟女国产| 超碰成人久久| 久久人人精品亚洲av| 天堂动漫精品| 国产精品偷伦视频观看了| 两人在一起打扑克的视频| 看免费av毛片| 亚洲中文字幕日韩| av在线播放免费不卡| 亚洲性夜色夜夜综合| 国产精品久久久久久人妻精品电影| 香蕉国产在线看| 性欧美人与动物交配| 成人国语在线视频| 亚洲成人免费av在线播放| 日韩一卡2卡3卡4卡2021年| 欧美成人性av电影在线观看| 亚洲性夜色夜夜综合| 久久精品影院6| 国产精品国产av在线观看| √禁漫天堂资源中文www| 熟女少妇亚洲综合色aaa.| 久久国产精品人妻蜜桃| 在线十欧美十亚洲十日本专区| 精品一区二区三区视频在线观看免费 | 国产精品电影一区二区三区| 水蜜桃什么品种好| 男人舔女人下体高潮全视频| 中文字幕人妻丝袜制服| tocl精华| 国产高清videossex| 久久国产亚洲av麻豆专区| 亚洲欧美激情在线| 他把我摸到了高潮在线观看| cao死你这个sao货| 亚洲第一欧美日韩一区二区三区| 亚洲国产欧美网| 每晚都被弄得嗷嗷叫到高潮| 久热这里只有精品99| 亚洲狠狠婷婷综合久久图片| 中文字幕人妻熟女乱码| 高潮久久久久久久久久久不卡| 国产精品乱码一区二三区的特点 | 亚洲精品粉嫩美女一区| 这个男人来自地球电影免费观看| 成年女人毛片免费观看观看9| 真人做人爱边吃奶动态| 国产精品一区二区精品视频观看| 99re在线观看精品视频| 黄频高清免费视频| www国产在线视频色| 麻豆一二三区av精品| 亚洲五月色婷婷综合| 在线观看一区二区三区激情| 久久中文字幕人妻熟女| 国产成人av教育| 精品久久久久久久久久免费视频 | 成在线人永久免费视频| 三级毛片av免费| 国产精品永久免费网站| 男男h啪啪无遮挡| 国产精品1区2区在线观看.| 啦啦啦免费观看视频1| 国产一区二区三区在线臀色熟女 | 日韩 欧美 亚洲 中文字幕| 亚洲av日韩精品久久久久久密| 国产精华一区二区三区| 久久青草综合色| 日韩中文字幕欧美一区二区| 久热这里只有精品99| 嫩草影院精品99| 国产熟女午夜一区二区三区| 亚洲欧美日韩高清在线视频| 好男人电影高清在线观看| 男女床上黄色一级片免费看| 一进一出抽搐动态| 精品一区二区三区av网在线观看| 午夜两性在线视频| 黄片小视频在线播放| 又大又爽又粗| 亚洲情色 制服丝袜| 国产精品一区二区精品视频观看| 日韩欧美一区视频在线观看| 久久人妻熟女aⅴ| 嫁个100分男人电影在线观看| 欧洲精品卡2卡3卡4卡5卡区| 热99国产精品久久久久久7| 国产99白浆流出| 午夜a级毛片| 少妇裸体淫交视频免费看高清 | 色婷婷久久久亚洲欧美| 久久久国产欧美日韩av| 啪啪无遮挡十八禁网站| 色综合欧美亚洲国产小说| 热re99久久国产66热| 国产99白浆流出| 欧美黑人精品巨大| 亚洲国产欧美一区二区综合| 久久国产精品人妻蜜桃| 少妇粗大呻吟视频| 极品教师在线免费播放| 一级片'在线观看视频| 韩国av一区二区三区四区| 欧美一区二区精品小视频在线| 午夜福利欧美成人| 免费日韩欧美在线观看| 欧美激情高清一区二区三区| 午夜福利一区二区在线看| 人人妻,人人澡人人爽秒播| 女人被躁到高潮嗷嗷叫费观| 不卡一级毛片| 欧美+亚洲+日韩+国产| 天天躁夜夜躁狠狠躁躁| 午夜福利欧美成人| 91成年电影在线观看| 日韩视频一区二区在线观看| 99热只有精品国产| 成年人免费黄色播放视频| 国产成人精品在线电影| 在线天堂中文资源库| 免费高清视频大片| 宅男免费午夜| 亚洲精品国产精品久久久不卡| 亚洲av电影在线进入| 国产精品偷伦视频观看了| 国产成人欧美在线观看| 欧美久久黑人一区二区| 国产99白浆流出| 日本精品一区二区三区蜜桃| 亚洲一卡2卡3卡4卡5卡精品中文| 夜夜躁狠狠躁天天躁| 亚洲欧美激情在线| 9热在线视频观看99| 久久久国产欧美日韩av| 国产免费现黄频在线看| 日韩一卡2卡3卡4卡2021年| 韩国av一区二区三区四区| 乱人伦中国视频| 日本黄色日本黄色录像| 欧美日韩一级在线毛片| 久久人人爽av亚洲精品天堂| 亚洲精品一卡2卡三卡4卡5卡| 久久狼人影院| 国产在线观看jvid| 午夜亚洲福利在线播放| 中出人妻视频一区二区| 曰老女人黄片| 午夜a级毛片| 日本a在线网址| 国产单亲对白刺激| 久久久久国产精品人妻aⅴ院| 亚洲成人免费av在线播放| 国产一区二区激情短视频| 国产国语露脸激情在线看| 国产成人免费无遮挡视频| 夜夜看夜夜爽夜夜摸 | 欧美日韩中文字幕国产精品一区二区三区 | 亚洲av电影在线进入| 无人区码免费观看不卡| 国产精品久久久久成人av| 人人妻,人人澡人人爽秒播| 精品国产乱码久久久久久男人| 热99re8久久精品国产| 中文字幕精品免费在线观看视频| 岛国视频午夜一区免费看| 真人一进一出gif抽搐免费| 50天的宝宝边吃奶边哭怎么回事| 久久久久亚洲av毛片大全| 91av网站免费观看| 黑人猛操日本美女一级片| 色综合婷婷激情| 九色亚洲精品在线播放| 婷婷精品国产亚洲av在线| 国产精品久久久人人做人人爽| 成人免费观看视频高清| 亚洲性夜色夜夜综合| videosex国产| 久久99一区二区三区| 国产精品偷伦视频观看了| 亚洲成国产人片在线观看| 日日夜夜操网爽| 欧美日韩精品网址| 一区二区三区激情视频| 亚洲国产精品999在线| 亚洲精品美女久久av网站| 日本免费a在线| 欧美乱妇无乱码| www.精华液| 国产成人免费无遮挡视频| 亚洲av五月六月丁香网| 国产av在哪里看| 精品一品国产午夜福利视频| 99精国产麻豆久久婷婷| 性少妇av在线| 黄色成人免费大全| 亚洲中文av在线| 桃色一区二区三区在线观看| 国产精品99久久99久久久不卡| 国产成年人精品一区二区 | 亚洲专区中文字幕在线| 国产精品av久久久久免费| 久久香蕉激情| 中文字幕高清在线视频| 国产在线观看jvid| 一边摸一边抽搐一进一小说| 亚洲色图综合在线观看| 精品人妻1区二区| 亚洲七黄色美女视频| 婷婷丁香在线五月| 久久婷婷成人综合色麻豆| 色综合站精品国产| 久久精品国产99精品国产亚洲性色 | 成人影院久久| 欧美亚洲日本最大视频资源| 叶爱在线成人免费视频播放| www.精华液| 亚洲一区二区三区色噜噜 | av天堂久久9| 亚洲国产欧美日韩在线播放| 狠狠狠狠99中文字幕| 久久热在线av| 女人高潮潮喷娇喘18禁视频| 男人的好看免费观看在线视频 | 欧美日韩中文字幕国产精品一区二区三区 | 久久久久国内视频| 免费在线观看日本一区| av在线播放免费不卡| 免费看十八禁软件| 国产精品野战在线观看 | 日韩欧美在线二视频| 久久婷婷成人综合色麻豆| 男女之事视频高清在线观看| 午夜成年电影在线免费观看| 欧美日韩亚洲高清精品| 日韩视频一区二区在线观看| 亚洲欧美精品综合久久99| 亚洲熟女毛片儿| 国产一区二区三区综合在线观看| 国产区一区二久久| 美女国产高潮福利片在线看| 国产欧美日韩精品亚洲av| 精品久久久久久电影网| 精品少妇一区二区三区视频日本电影| av欧美777| 亚洲国产欧美一区二区综合| 一级毛片女人18水好多| 91字幕亚洲| 99国产综合亚洲精品| 国产精品久久电影中文字幕| 亚洲一区二区三区色噜噜 | 亚洲精品av麻豆狂野| 亚洲国产精品999在线| 国产成人影院久久av| 久久久久精品国产欧美久久久| 一进一出好大好爽视频| 久久久久九九精品影院| 9191精品国产免费久久| 久久久国产欧美日韩av| 国产av又大| netflix在线观看网站| 国产xxxxx性猛交| 亚洲五月婷婷丁香| 国产1区2区3区精品| 午夜a级毛片| 如日韩欧美国产精品一区二区三区| 免费一级毛片在线播放高清视频 | av网站免费在线观看视频| 亚洲欧美日韩高清在线视频| 久久久国产欧美日韩av| 波多野结衣一区麻豆| 中文字幕高清在线视频| 亚洲av美国av| 高清在线国产一区| 久久欧美精品欧美久久欧美| www.www免费av| netflix在线观看网站| 免费在线观看完整版高清| 国产麻豆69| 亚洲国产精品合色在线| 琪琪午夜伦伦电影理论片6080| 午夜久久久在线观看| 久久久久久人人人人人| 国产色视频综合| 欧美国产精品va在线观看不卡| 欧美激情极品国产一区二区三区| 国产亚洲欧美精品永久| 女人高潮潮喷娇喘18禁视频| 大型av网站在线播放| 国产精品一区二区免费欧美| 韩国精品一区二区三区| 精品福利永久在线观看| 国产精品电影一区二区三区| 波多野结衣高清无吗| 国产精品 国内视频| 美女福利国产在线| 少妇被粗大的猛进出69影院| 欧美一级毛片孕妇| 12—13女人毛片做爰片一| 狠狠狠狠99中文字幕| 久久人人爽av亚洲精品天堂| 婷婷精品国产亚洲av在线| 亚洲五月婷婷丁香| 手机成人av网站| 19禁男女啪啪无遮挡网站| 国产精品久久电影中文字幕| 婷婷精品国产亚洲av在线| 午夜精品在线福利| 欧美不卡视频在线免费观看 | 久久久久久久久中文| 亚洲成人国产一区在线观看| 美女扒开内裤让男人捅视频| 亚洲伊人色综图| 后天国语完整版免费观看| 欧美不卡视频在线免费观看 | 国产伦一二天堂av在线观看| 亚洲久久久国产精品| 亚洲一区高清亚洲精品| 色播在线永久视频| 国产精品综合久久久久久久免费 | 日韩欧美一区视频在线观看| 嫩草影院精品99| 视频区图区小说| 窝窝影院91人妻| 老熟妇乱子伦视频在线观看| 亚洲av成人不卡在线观看播放网| 亚洲精品一二三| 日韩人妻精品一区2区三区| 国产亚洲欧美精品永久| 亚洲自偷自拍图片 自拍| 日韩免费av在线播放| 极品教师在线免费播放| 女人精品久久久久毛片| 一本大道久久a久久精品| 美女大奶头视频| 在线看a的网站| 高清在线国产一区| 亚洲五月色婷婷综合| 日韩欧美三级三区| 久久久国产成人免费| 搡老熟女国产l中国老女人| 亚洲国产精品sss在线观看 | 日本一区二区免费在线视频| 一级毛片精品| 午夜精品国产一区二区电影| 国产一区二区在线av高清观看| 国产色视频综合| 久久午夜综合久久蜜桃| 国产99白浆流出| 麻豆av在线久日| 久久精品国产综合久久久| 亚洲第一青青草原| 男女高潮啪啪啪动态图| 美女午夜性视频免费| 午夜老司机福利片| 成人av一区二区三区在线看| 国产精华一区二区三区| 老熟妇乱子伦视频在线观看| 免费女性裸体啪啪无遮挡网站| 国产成人av激情在线播放| 亚洲第一av免费看| 18禁黄网站禁片午夜丰满| 亚洲人成网站在线播放欧美日韩| 叶爱在线成人免费视频播放| 精品卡一卡二卡四卡免费| av视频免费观看在线观看| 国产精品秋霞免费鲁丝片| 亚洲中文字幕日韩| 欧美日韩亚洲综合一区二区三区_| 精品第一国产精品| 免费av毛片视频| 亚洲欧美激情综合另类| 妹子高潮喷水视频| 国产成人精品无人区| 国产精品野战在线观看 | 两人在一起打扑克的视频| 国产精品电影一区二区三区| 老司机靠b影院| 国产免费现黄频在线看| 亚洲五月婷婷丁香| 最近最新中文字幕大全电影3 | 51午夜福利影视在线观看| 老司机午夜福利在线观看视频| 高清av免费在线| 夜夜看夜夜爽夜夜摸 | 日韩有码中文字幕| 亚洲国产欧美网| 涩涩av久久男人的天堂| 亚洲欧美一区二区三区黑人| 麻豆av在线久日| 国产人伦9x9x在线观看| 岛国在线观看网站| 在线观看免费视频网站a站| 久久亚洲真实| 丰满迷人的少妇在线观看| 纯流量卡能插随身wifi吗| 亚洲性夜色夜夜综合| 成人影院久久| 嫩草影视91久久| 亚洲av熟女| 9色porny在线观看| 男人舔女人的私密视频| 一夜夜www| 日韩大尺度精品在线看网址 | 久久人人精品亚洲av| 桃红色精品国产亚洲av| 日日摸夜夜添夜夜添小说| 国产高清激情床上av| 在线观看免费午夜福利视频| 国产精品久久久久成人av| 丝袜美腿诱惑在线| 国产成人免费无遮挡视频| 美女高潮到喷水免费观看| 免费一级毛片在线播放高清视频 | 久久久精品国产亚洲av高清涩受| 午夜精品国产一区二区电影| 性少妇av在线| 91成人精品电影| 成人黄色视频免费在线看| 亚洲一区二区三区色噜噜 | 国产精品爽爽va在线观看网站 | 黄色怎么调成土黄色| 日韩免费av在线播放| 大型av网站在线播放| 黑人欧美特级aaaaaa片| 亚洲专区中文字幕在线| 亚洲狠狠婷婷综合久久图片| x7x7x7水蜜桃| 亚洲专区国产一区二区| 亚洲欧美日韩另类电影网站| 精品一区二区三卡| 嫩草影院精品99| 久久精品亚洲av国产电影网| 国产精品电影一区二区三区| 伊人久久大香线蕉亚洲五| 亚洲av第一区精品v没综合| 久久人人爽av亚洲精品天堂| 欧美亚洲日本最大视频资源| 琪琪午夜伦伦电影理论片6080| 成人亚洲精品av一区二区 | 久久久久久久久久久久大奶| 日韩国内少妇激情av| 少妇裸体淫交视频免费看高清 | 91精品三级在线观看| 久久影院123| 久热这里只有精品99| 正在播放国产对白刺激| 亚洲精品美女久久久久99蜜臀| 国产99久久九九免费精品| 国产三级黄色录像| 亚洲av成人一区二区三| 国产黄a三级三级三级人| 男男h啪啪无遮挡| 亚洲中文av在线| 国产成人精品久久二区二区免费| 亚洲av电影在线进入| 91九色精品人成在线观看| 51午夜福利影视在线观看| 两性午夜刺激爽爽歪歪视频在线观看 | 国产单亲对白刺激| 亚洲欧美日韩无卡精品| av网站在线播放免费| 国产av在哪里看| 十八禁人妻一区二区| 午夜日韩欧美国产| 国产精品综合久久久久久久免费 | 麻豆av在线久日| 亚洲 欧美一区二区三区| 亚洲第一av免费看| 国产色视频综合| 美国免费a级毛片| 日韩欧美三级三区| 免费久久久久久久精品成人欧美视频| 欧美成人午夜精品| 日韩欧美一区视频在线观看| 国产麻豆69| 欧美最黄视频在线播放免费 | 午夜视频精品福利| 18禁美女被吸乳视频| 午夜免费鲁丝| 日本 av在线| 日韩有码中文字幕| 色哟哟哟哟哟哟| 久久久久久久久久久久大奶| 热re99久久国产66热| 他把我摸到了高潮在线观看| 欧美成人免费av一区二区三区| 亚洲av五月六月丁香网| 国产亚洲精品久久久久5区| 在线国产一区二区在线| 少妇被粗大的猛进出69影院| 亚洲欧美精品综合久久99| 久久精品91蜜桃| 高清黄色对白视频在线免费看| 热re99久久国产66热| 在线天堂中文资源库| 亚洲精品国产一区二区精华液| 美国免费a级毛片| 亚洲成av片中文字幕在线观看| 久久久久国产一级毛片高清牌| 午夜精品国产一区二区电影| 一个人观看的视频www高清免费观看 | 亚洲精品在线观看二区| 精品熟女少妇八av免费久了| 丝袜在线中文字幕| 日韩人妻精品一区2区三区| 老司机深夜福利视频在线观看| 日本黄色日本黄色录像| 亚洲人成电影免费在线| 久久精品影院6| 亚洲第一欧美日韩一区二区三区| 12—13女人毛片做爰片一| 丝袜美腿诱惑在线| 搡老乐熟女国产| 午夜福利一区二区在线看| 亚洲一卡2卡3卡4卡5卡精品中文| 在线观看www视频免费| 亚洲五月婷婷丁香| 亚洲av五月六月丁香网| 91成年电影在线观看| 中文字幕最新亚洲高清| 国产精品影院久久| 久久午夜亚洲精品久久| 不卡一级毛片| 一区在线观看完整版| 露出奶头的视频| 国产亚洲欧美精品永久| 91字幕亚洲| 亚洲aⅴ乱码一区二区在线播放 | 亚洲情色 制服丝袜| 搡老岳熟女国产| 国产精品偷伦视频观看了| 国产亚洲精品一区二区www| 99精品欧美一区二区三区四区| 桃色一区二区三区在线观看| 欧美亚洲日本最大视频资源| 在线看a的网站| 国产激情欧美一区二区| 久久中文看片网| 欧美日韩亚洲国产一区二区在线观看| 欧美+亚洲+日韩+国产| 麻豆国产av国片精品| 色婷婷久久久亚洲欧美| 婷婷六月久久综合丁香| 极品人妻少妇av视频| 18禁观看日本| 久久精品亚洲av国产电影网| 日韩免费av在线播放| e午夜精品久久久久久久| 巨乳人妻的诱惑在线观看| 最新美女视频免费是黄的| 在线天堂中文资源库| 91九色精品人成在线观看| 久久精品成人免费网站| 欧美亚洲日本最大视频资源| av片东京热男人的天堂| 一进一出好大好爽视频| 精品国产国语对白av| 在线观看午夜福利视频| av天堂在线播放| 国产精品一区二区在线不卡| 99精国产麻豆久久婷婷| 亚洲欧美日韩另类电影网站| 免费在线观看亚洲国产| 亚洲视频免费观看视频| 丝袜美足系列| 中亚洲国语对白在线视频| 国产深夜福利视频在线观看| 久久久国产一区二区| 男女做爰动态图高潮gif福利片 | 一本大道久久a久久精品| 丝袜人妻中文字幕| 婷婷精品国产亚洲av在线| 亚洲国产中文字幕在线视频| 在线天堂中文资源库| 9热在线视频观看99| 69精品国产乱码久久久| 色在线成人网| 精品一区二区三区av网在线观看| 亚洲九九香蕉| 热re99久久精品国产66热6| 三级毛片av免费|