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

    強(qiáng)非共沸工質(zhì)R1234yf/R170/R14系氣液相平衡的模擬

    2018-03-27 09:17:11
    制冷學(xué)報(bào) 2018年2期
    關(guān)鍵詞:狀態(tài)方程制冷系統(tǒng)工質(zhì)

    (上海理工大學(xué)制冷與低溫工程研究所 上海 200093)

    自動(dòng)復(fù)疊制冷系統(tǒng)在小型制冷裝置、天然氣液化中具有較強(qiáng)的優(yōu)勢(shì),能制取氮?dú)庖夯瘻囟?7 K到常規(guī)單機(jī)壓縮制冷溫度230 K,常被用于氣體液化、半導(dǎo)體加工、低溫醫(yī)學(xué)和低溫生物等領(lǐng)域。國(guó)內(nèi)外對(duì)自復(fù)疊節(jié)流制冷系統(tǒng)的研究較多,從兩級(jí)到五級(jí),從二元混合制冷劑到五元混合制冷劑。一般而言,200~230 K采用兩級(jí)自復(fù)疊制冷系統(tǒng),160~200 K采用三級(jí)自復(fù)疊制冷系統(tǒng),120~160 K采用四級(jí)自復(fù)疊系統(tǒng),70~120 K采用采用五級(jí)自復(fù)疊系統(tǒng)[1]。M. Ruhemann[2]將氟利昂制冷劑應(yīng)用于自動(dòng)復(fù)疊制冷系統(tǒng),采用R13/R22混合制冷劑取得了200 K的低溫。A. P. Kleemenko[3]采用摩爾比為65%/20%/15%的甲烷、乙烷和正丁烷的混合制冷劑通過(guò)單級(jí)壓縮兩級(jí)分離獲得的低溫成功液化了天然氣,該制冷循環(huán)具有很高的熱力學(xué)效率和運(yùn)行可靠性,因此成為單流混合工質(zhì)復(fù)疊循環(huán)(mixed refrigerant cascade cycle,MRC)。A. Fuderer[4]采用質(zhì)量分?jǐn)?shù)為20%/80%的R50/R12混合制冷劑成功獲取了117 K的低溫。此后,眾多研究者不斷改進(jìn)自復(fù)疊系統(tǒng),不斷提高制冷效率,能獲取的溫度也越來(lái)越低[5-8]。李文林等[9-10]展開(kāi)了對(duì)自動(dòng)復(fù)疊制冷系統(tǒng)的研究,以R12/R13作為混合制冷劑,在70 L的低溫箱中得到了218 K的溫度,熱力學(xué)效率為11.4%。Luo Ercang等[11-14]對(duì)混合制冷劑做了大量的研究工作,包括混合制冷劑熱物性計(jì)算、系統(tǒng)優(yōu)化和氣液相平衡等方面。王國(guó)棟等[15-16]采用以R22/R23為混合制冷劑、單級(jí)壓縮、單級(jí)分凝的自動(dòng)復(fù)疊制冷系統(tǒng),在4 h內(nèi)成功將150 L的低溫箱內(nèi)溫度降低到213 K以下,此后又以R134a/R23/R14為混合制冷劑,將冷箱溫度降低到173 K,制冷量為38 W。此外,還對(duì)R600a/R23[17]、R600a/R23/R14/R740[18]、R245fa/R600a/R508B/R14[19]等混合制冷劑的自復(fù)疊系統(tǒng)做了相關(guān)實(shí)驗(yàn)研究,并對(duì)系統(tǒng)進(jìn)行了改進(jìn),分析了系統(tǒng)降溫特性、混合物分離匯合特性以及混合物的優(yōu)化配比等多方面內(nèi)容。隨著人們對(duì)全球變暖、溫室效應(yīng)、臭氧層破壞等世界性環(huán)境問(wèn)題的日益關(guān)注,制冷行業(yè)使用的制冷工質(zhì)也加快了更新?lián)Q代的步伐。自復(fù)疊系統(tǒng)中使用的制冷劑大多為氯氟烴(CFCs)和氫氯氟烴(HCFCs)類制冷工質(zhì),它們或?qū)Υ髿獬粞鯇佑衅茐淖饔?,或具有很?qiáng)的溫室效應(yīng),或兩者兼有。具有高全球變暖潛能值(GWP)和消耗臭氧潛能值(ODP)的制冷劑都將被逐步淘汰,新型環(huán)保制冷劑將逐步出現(xiàn)。

    自復(fù)疊制冷系統(tǒng)結(jié)構(gòu)相對(duì)簡(jiǎn)單,但混合制冷劑的工作機(jī)理復(fù)雜,涉及到很多理論基礎(chǔ)知識(shí),例如工質(zhì)的分離與混合、混合工質(zhì)傳熱、氣液相平衡等。其中氣液相平衡的研究尤為重要,它是對(duì)混合物基本物性的研究,以確定最佳工作狀態(tài)、最佳配比等,影響整個(gè)系統(tǒng)的可靠性及運(yùn)行效率。已有許多學(xué)者對(duì)混合制冷劑的氣液相平衡進(jìn)行了實(shí)驗(yàn)和模擬研究,但多數(shù)集中在近共沸工質(zhì),針對(duì)適用于多級(jí)自復(fù)疊系統(tǒng)的強(qiáng)非共沸工質(zhì)的氣液相平衡研究較少。Hu Peng等[20-25]對(duì)多組混合工質(zhì)進(jìn)行了氣液相平衡研究并取得了一些相平衡數(shù)據(jù),包括R1234yf/R600a[20]、R1234yf/R152a[21],R1234yf/R161[22]、R134a/R600a/R1234yf[23]、R134a/R227ea和R143 a/R1234yf[24],并通過(guò)實(shí)驗(yàn)數(shù)據(jù)結(jié)合理論模型對(duì)相平衡計(jì)算中的重要參數(shù)——二元相互作用系數(shù)kij進(jìn)行了歸納總結(jié)[25]。J. Lia等[26]利用狀態(tài)方程法對(duì)R125/R134a/R143a系進(jìn)行了模擬研究,指出狀態(tài)方程法能很好的預(yù)測(cè)混合物相平衡特性。Gong Maoqiong等[27-31]針對(duì)含R290的混合物進(jìn)行了相關(guān)理論和實(shí)驗(yàn)研究。R. Budinsky等[32]采用Gibbs系綜Monto Carlo(GEMC)模擬方法對(duì)兩組混合物(R134a/R125和R134a/R32)進(jìn)行了計(jì)算,計(jì)算結(jié)果與文獻(xiàn)數(shù)據(jù)具有較好的一致性,并指出GEMC與以熱力學(xué)為基礎(chǔ)的Wilson和UNIFAC方法的預(yù)測(cè)精度相當(dāng)。Q. N. Ho等[33-34]搭建了相關(guān)實(shí)驗(yàn)臺(tái),對(duì)二元混合物R1270/R134a和碳?xì)涔べ|(zhì)對(duì)R1270/R290的氣液相平衡數(shù)據(jù)進(jìn)行了實(shí)驗(yàn),發(fā)現(xiàn)實(shí)驗(yàn)數(shù)據(jù)與采用Peng-Robinson狀態(tài)方程(PR-EoS)結(jié)合Wong-Sandler(WS)混合法則的模擬結(jié)果吻合較好。以上學(xué)者多以實(shí)驗(yàn)的方法對(duì)混合物的氣液相平衡進(jìn)行研究。雖然實(shí)驗(yàn)方法是獲取氣液相平衡數(shù)據(jù)的第一手段,但僅僅通過(guò)實(shí)驗(yàn)無(wú)法取得大量的相平衡數(shù)據(jù),因此通過(guò)有限的實(shí)驗(yàn)數(shù)據(jù)發(fā)展理論模型,提高預(yù)測(cè)精度,也是獲取相平衡數(shù)據(jù)的重要途徑[35]。

    綜合考慮各種制冷劑的物性參數(shù),并結(jié)合三級(jí)自復(fù)疊制冷系統(tǒng)的特點(diǎn)和溫度范圍,本文對(duì)適用于三級(jí)自復(fù)疊制冷系統(tǒng)的R1234yf/R170/R14系混合制冷劑的氣液相平衡數(shù)據(jù)進(jìn)行了模擬計(jì)算。其中R1234yf是近年來(lái)被廣泛使用新型制冷劑,且GWP<1,ODP=0,具有非常短的大氣壽命(0.029 年),被公認(rèn)為是很有發(fā)展前景的環(huán)保替代制冷劑。R170為碳?xì)渲评鋭?,是天然工質(zhì),同時(shí)具備節(jié)能與環(huán)保的兩大優(yōu)勢(shì)。R14的ODP=0,但GWP相對(duì)較高,雖然不是很理想的選擇,但在目前情況下很難找到其對(duì)應(yīng)溫度段更好的替代工質(zhì)。在理論模型的選擇方面,超額吉布斯自由能-狀態(tài)方程GE-EoS是繼傳統(tǒng)的狀態(tài)方程法和活度系數(shù)法之后預(yù)測(cè)相平衡的一個(gè)新思路,此模型結(jié)合了狀態(tài)方程法和活度系數(shù)法在相平衡預(yù)測(cè)方面的優(yōu)點(diǎn),并將極性體系預(yù)測(cè)能力非常強(qiáng)的活度系數(shù)模型直接應(yīng)用于狀態(tài)方程法的相平衡預(yù)測(cè)中。本文基于PR狀態(tài)方程與WS混合法則,結(jié)合PSRK(Predictive-Soave-Redlich-Kwong)方程中使用的UNIFAC基團(tuán)貢獻(xiàn)法構(gòu)建了混合物氣液相平衡預(yù)測(cè)模型(PRWS-UNIFAC-PSRK)。利用此模型計(jì)算了二元混合制冷劑R161/R1234yf及三元混合制冷劑R32/R125/R134a的氣液相平衡數(shù)據(jù),并與相關(guān)文獻(xiàn)的實(shí)驗(yàn)數(shù)據(jù)進(jìn)行對(duì)比,驗(yàn)證了模型的精確度。在此基礎(chǔ)上,對(duì)還未見(jiàn)有相平衡實(shí)驗(yàn)數(shù)據(jù)報(bào)道的R1234yf/R170/R14系混合物進(jìn)行了模擬計(jì)算和分析,并給出了相關(guān)的相平衡圖。這些相平衡特性對(duì)混合物的優(yōu)化分離及相關(guān)自動(dòng)復(fù)疊系統(tǒng)設(shè)計(jì)及改進(jìn)都具有指導(dǎo)和借鑒意義。

    1 方法及驗(yàn)證

    1.1 模型選擇

    流體介質(zhì)最基本的平衡物性參數(shù)即為壓力-比容-溫度(pvt)三者之間的關(guān)系,是研究流體其它熱物性參數(shù)的基礎(chǔ)。而體現(xiàn)流體pvt性質(zhì)的方程式就是狀態(tài)方程,因而對(duì)流體熱物性的研究最終可歸為對(duì)狀態(tài)方程的獲取。其中,兩參數(shù)的立方型狀態(tài)方程因其參數(shù)少,形式簡(jiǎn)單,計(jì)算精度高而備受研究人員的重視。應(yīng)用狀態(tài)方程法預(yù)測(cè)混合物氣液相平衡特性,需要結(jié)合合適的混合法則,目前大多數(shù)文獻(xiàn)中使用的模擬方法為van der Waals(vdW)法則,但是其在預(yù)測(cè)非理性體系時(shí)的精度并不高。為了克服這一缺陷,采用混合法則結(jié)合立方形狀態(tài)方程及超額吉布斯自由能越來(lái)越多的在混合物氣液相平衡的預(yù)測(cè)中使用。目前來(lái)看,還沒(méi)有采用超額吉布斯自由能-狀態(tài)方程模型(GE-EoS)對(duì)強(qiáng)非共沸的三元混合工質(zhì)的氣液相平衡特性的研究報(bào)告。本文利用PR狀態(tài)方程和WS混合法則,結(jié)合PSRK方程中使用的UNIFAC活度系數(shù)模型,編譯了相關(guān)模擬程序,對(duì)多組二元及三元混合物的氣液相平衡數(shù)據(jù)進(jìn)行了計(jì)算,并與參考文獻(xiàn)數(shù)據(jù)和美國(guó)國(guó)家標(biāo)準(zhǔn)與技術(shù)研究院(NIST)數(shù)據(jù)庫(kù)進(jìn)行對(duì)比,驗(yàn)證了模型的準(zhǔn)確性。在此基礎(chǔ)上,對(duì)強(qiáng)非共沸工質(zhì)R1234yf/R170/R14系的氣液相平衡特性進(jìn)行了模擬研究。

    PR方程[36]的形式為:

    (1)

    (2)

    (3)

    kω=0.374 64+1.542 26ω-0.269 92ω2

    (4)

    表1 純工質(zhì)基本物性參數(shù)Tab.1 Basic physical parameters of pure components

    WS混合法則[37]的基本形式為:

    (5)

    (6)

    (7)

    (8)

    由PR方程和WS混合規(guī)則推導(dǎo)所得逸度系數(shù)求解方程為:

    (9)

    (10)

    (11)

    (12)

    (13)

    活度系數(shù)γi由改進(jìn)的PSRK中使用的UNIFAC模型[38]計(jì)算得到:

    lnγi=lnγiC+lnγiR

    (14)

    (15)

    (16)

    (17)

    (18)

    (19)

    式中:Ni為混合物2組分的分子數(shù);N為混合物總的分子數(shù)。用于WS混合法則中相關(guān)參數(shù)計(jì)算的各組分二元相互作用系數(shù)kij列于表2,它們由表3中的文獻(xiàn)數(shù)據(jù)回歸擬合得到。

    表2 二元混合物的相互作用系數(shù)kijTab.2 Regressed kij from binary system

    表3 混合物文獻(xiàn)匯總Tab.3 Literature summary for mixtures

    由于R1234yf/R170/R14系中各二元體系的氣液相平衡數(shù)據(jù)未見(jiàn)有公開(kāi)報(bào)道,因而相互作用系數(shù)采用Gong Maoqiong等[39-40]擬合的數(shù)值。用于活度系數(shù)和吉布斯自由能計(jì)算的UNIFAC-PSRK模型相關(guān)基團(tuán)參數(shù)列于表4和表5?;鶊F(tuán)交互作用系數(shù)Amk是基團(tuán)k和基團(tuán)m之間相互作用能與兩個(gè)m(k)基團(tuán)之間相互作用能差異的度量(Amk≠Akm),例如基團(tuán)m為CH2,基團(tuán)k為CF2時(shí),Amk=42.257,Akm=-7.474。雖然表5中每一個(gè)主基團(tuán)下子基團(tuán)的Rk和Qk的值不等,但在同一主基團(tuán)中所有子基團(tuán)的交互作用系數(shù)是相同的,因而絕大部分的HFCs和HFO制冷劑的相平衡數(shù)據(jù)都可以通過(guò)三個(gè)主基團(tuán)CH2、CF2和F的9個(gè)交互作用系數(shù)進(jìn)行預(yù)測(cè)。本文參考吳獻(xiàn)忠等[41]擬合的基團(tuán)交互作用系數(shù),對(duì)多組二元及三元混合工質(zhì)進(jìn)行氣液相平衡預(yù)測(cè)。

    表4 UNIFAC模型中的基團(tuán)相互作用系數(shù)Amk和Akm[41]

    表5 基團(tuán)體積參數(shù)Rk和表面積參數(shù)Qk[41]

    1.2 模型驗(yàn)證

    為驗(yàn)證模型計(jì)算精度,本文選取已有文獻(xiàn)實(shí)驗(yàn)數(shù)據(jù)的二元混合工質(zhì)R161/R1234yf進(jìn)行模擬計(jì)算,并將模擬結(jié)果與文獻(xiàn)實(shí)驗(yàn)數(shù)據(jù)和美國(guó)國(guó)家標(biāo)準(zhǔn)與技術(shù)研究院(NIST)制定的REFPROP9.0軟件數(shù)據(jù)庫(kù)進(jìn)行對(duì)比分析,如表6所示。

    圖1所示為R161/R1234yf系從283.15~323.15 K溫度范圍內(nèi)壓力與組分關(guān)系的對(duì)比,其中五角星和圓圈分別為實(shí)驗(yàn)測(cè)量的露點(diǎn)壓力值和泡點(diǎn)壓力值,實(shí)線和虛線分別為模型計(jì)算結(jié)果的壓力露點(diǎn)線和泡點(diǎn)線??梢钥闯觯瑢?shí)驗(yàn)數(shù)據(jù)點(diǎn)與計(jì)算結(jié)果的泡點(diǎn)和露點(diǎn)線吻合較好。圖2所示為模型計(jì)算結(jié)果與實(shí)驗(yàn)數(shù)據(jù)的偏差,壓力計(jì)算結(jié)果與實(shí)驗(yàn)值的相對(duì)誤差基本在±1.5%以內(nèi),氣相組分質(zhì)量分?jǐn)?shù)的絕對(duì)誤差均在±0.02內(nèi),表明模型對(duì)二元混合工質(zhì)氣液相平衡特性的預(yù)測(cè)精度較高,滿足實(shí)際工程應(yīng)用的需求。

    圖1 R161(1)/R1234yf(2)系計(jì)算值與實(shí)驗(yàn)值的壓力-組分關(guān)系Fig.1 Pressure-composition relationship between calculated results and experimental data for the R161(1)/R1234yf(2) system

    圖2 R161(1)/R1234yf(2)系計(jì)算值與實(shí)驗(yàn)值的誤差關(guān)系Fig.2 Deviation between calculated results and experimental data for the R161(1)/R1234yf(2) system

    圖3 R161(1)/R1234yf(2)系計(jì)算值與REFPROP9.0數(shù)據(jù)庫(kù)的壓力-組分關(guān)系Fig.3 Pressure-composition relationship between calculated results and REFPROP9.0 database for the R161(1)/R1234yf(2) system

    將模擬結(jié)果與REFPROP9.0數(shù)據(jù)庫(kù)進(jìn)行對(duì)比,如圖3和圖4所示。壓力計(jì)算值與數(shù)據(jù)庫(kù)的相對(duì)誤差都在±5%以內(nèi),氣相組分質(zhì)量分?jǐn)?shù)的絕對(duì)誤差基本在±0.012 以內(nèi),表明模型的計(jì)算精度較高。觀察圖2可知,當(dāng)R161的質(zhì)量分?jǐn)?shù)小于0.5時(shí),壓力偏差主要表現(xiàn)為正偏差,即實(shí)驗(yàn)值大于計(jì)算值,而氣相組分質(zhì)量分?jǐn)?shù)主要表現(xiàn)為負(fù)偏差,即實(shí)驗(yàn)值小于計(jì)算值;當(dāng)R161的質(zhì)量分?jǐn)?shù)大于0.5時(shí),變化趨勢(shì)剛好與上述情況相反。分析其原因如下:1)參考文獻(xiàn)本身的測(cè)量誤差;2)由于WS混合規(guī)則直接采用由低壓氣液平衡數(shù)據(jù)得到的活度系數(shù)模型參數(shù)和二元相互作用系數(shù)kij直接推算高壓區(qū)氣液相平衡,雖有良好的溫度和壓力外推性能,但依然會(huì)有一定的誤差;3)采用PRWS-UNIFAC-PSRK模型對(duì)混合物進(jìn)行氣液相平衡預(yù)測(cè)時(shí),需要建立在對(duì)混合工質(zhì)基團(tuán)分析的基礎(chǔ)上,而UNIFAC基團(tuán)貢獻(xiàn)法主要應(yīng)用于液相逸度的計(jì)算,氣相逸度的計(jì)算還是采用狀態(tài)方程,而狀態(tài)方程在預(yù)測(cè)大分子化合物時(shí)比較困難,因而會(huì)出現(xiàn)負(fù)偏差的情況。因此,采用更高精度的狀態(tài)方程以及對(duì)UNIFAC基團(tuán)貢獻(xiàn)法及參數(shù)進(jìn)行更深入全面的分析研究,將會(huì)提高模型的預(yù)測(cè)精度和適用范圍。

    表6 溫度為283.15~323.15 K時(shí)R161(1)/R1234yf(2)系氣液相平衡數(shù)據(jù)

    續(xù)表6

    圖4 R161(1)/R1234yf(2)系計(jì)算值與REFPROP9.0數(shù)據(jù)庫(kù)的誤差關(guān)系Fig.4 Deviation between calculated results and REFPROP9.0 database for the R161(1)/R1234yf(2) system

    圖5 R32(1)/R125(2)/R134a(3)系計(jì)算值與實(shí)驗(yàn)值的壓力偏差關(guān)系Fig.5 Deviation of the pressures between the experimental data and the calculated data for the R32(1)/R125(2)/R134a(3) system

    圖6 R32(1)/R125(2)/R134a(3)系計(jì)算值與實(shí)驗(yàn)值的組分質(zhì)量分?jǐn)?shù)偏差關(guān)系Fig.6 Deviation of the vapor phase mass fraction between the experimental data and the calculated data for the R32(1)/R125(2)/R134a(3) system

    為了驗(yàn)證模型對(duì)三元混合物氣液相平衡的預(yù)測(cè)精度,對(duì)三元混合物R32/R125/R134a在組分質(zhì)量配比分別0.272 1/0.126 8/0.601 1、0.191 0/0.433 2/0.375 8、0.097 6/0.705 7/0.196 7下的氣液相平衡特性進(jìn)行計(jì)算,并與文獻(xiàn)[45]的實(shí)驗(yàn)數(shù)據(jù)進(jìn)行對(duì)比,對(duì)比結(jié)果列于表7。圖5所示為壓力計(jì)算值與實(shí)驗(yàn)值的誤差關(guān)系,圖6所示為氣相組分質(zhì)量分?jǐn)?shù)計(jì)算值與實(shí)驗(yàn)值的誤差。由圖5和圖6可知,氣相組分和壓力的模擬結(jié)果與實(shí)驗(yàn)數(shù)據(jù)的偏差基本在±0.02和±4%以內(nèi),表明模型在預(yù)測(cè)三元混合物的時(shí)候亦具有較高的精度,滿足實(shí)際工程應(yīng)用的精度要求。雖然模型具有較好的預(yù)測(cè)效果,但不論是以無(wú)窮壓力或零壓力為參考態(tài)的GE-EoS模型在計(jì)算和預(yù)測(cè)非對(duì)稱體系的氣液相平衡時(shí)都會(huì)有一定的偏差,因此對(duì)GE-EoS模型進(jìn)行進(jìn)一步的擴(kuò)展研究,如對(duì)UNIFAC的參數(shù)表進(jìn)行修補(bǔ)和增訂,以期GE-EoS模型在計(jì)算非對(duì)稱體系和消除組合項(xiàng)差異時(shí)具有更好的預(yù)測(cè)效果,將會(huì)是下一步工作的重點(diǎn)。

    2 R1234yf/R170/R14系等壓氣液相平衡研究

    通過(guò)前文對(duì)二元混合物R161/R1234yf和三元混合物R32/R125/R134a的氣液相平衡特性的研究,對(duì)比計(jì)算結(jié)果與實(shí)驗(yàn)數(shù)據(jù),發(fā)現(xiàn)計(jì)算值與實(shí)驗(yàn)值具有很好的一致性。雖然在三元混合物的計(jì)算中直接采用二元體系中的相互作用系數(shù),但模擬結(jié)果具有很高的可信度,表明模型在預(yù)測(cè)二元及三元混合制冷劑氣液相平衡時(shí)具有較好的效果。在此基礎(chǔ)上,利用PRWS-UNIFAC-PSRK模型對(duì)強(qiáng)非共沸工質(zhì)R1234yf/R170/R14系的氣液相平衡特性進(jìn)行模擬研究。

    圖7 0.2 MPa壓力下R1234yf/R170/R14系的三維立體相圖Fig.7 Three-dimensional phase equilibria diagrams of R1234yf/R170/R14 ternary system at 0.2 MPa

    圖7~圖10分別為R1234yf/R170/R14系在壓力為0.2、0.5、1.5、2.2 MPa下的三維立體相圖,其中靠上一面為露點(diǎn)面。露點(diǎn)面以上部分為氣相區(qū);靠下一面為泡點(diǎn)面,泡點(diǎn)面以下部分為液相區(qū);泡點(diǎn)面和露點(diǎn)面之間的區(qū)域?yàn)闅庖簝上鄥^(qū)。從圖中可以很直觀的看到混合物組分、溫度、壓力及泡露點(diǎn)之間關(guān)系,紅色越深,表明溫度值越高。由于R1234yf、R170和R14三種制冷劑兩兩之間互為強(qiáng)非共沸工質(zhì),因而可以看到,泡點(diǎn)面與露點(diǎn)面之間的空間很大,即氣液兩相區(qū)范圍很大,同時(shí)意味著混合物具有很大的滑移溫度。通過(guò)數(shù)據(jù)篩選,混合物質(zhì)量分?jǐn)?shù)比為0.4/0.2/0.4時(shí)體系的溫度滑移現(xiàn)象最為明顯,最大的滑移溫度達(dá)到72.5 K。對(duì)比圖7和圖10可發(fā)現(xiàn),隨著系統(tǒng)壓力的增大,泡點(diǎn)溫度和露點(diǎn)溫度也隨之上升,且R1234yf組分的質(zhì)量分?jǐn)?shù)越大,泡露點(diǎn)溫度也越大。

    表7 R32(1)/R125(2)/R134a(3)系氣液相平衡數(shù)據(jù)[45]

    圖8 0.5 MPa壓力下R1234yf/R170/R14系的三維立體相圖Fig.8 Three-dimensional phase equilibria diagrams of R1234yf/R170/R14 ternary system at 0.5 MPa

    圖9 1.5 MPa壓力下R1234yf/R170/R14系的三維立體相圖Fig.9 Three-dimensional phase equilibria diagrams of R1234yf/R170/R14 ternary system at 1.5 MPa

    圖10 2.2 MPa壓力下R1234yf/R170/R14系的三維立體相圖Fig.10 Three-dimensional phase equilibria diagrams of R1234yf/R170/R14 ternary system at 2.2 MPa

    在三級(jí)自復(fù)疊系統(tǒng)中,每級(jí)循環(huán)中兩組分的標(biāo)準(zhǔn)沸點(diǎn)差值在40~80 K較合適,通過(guò)分析R1234yf/R170/R14系三維立體相圖,可大致判斷適用于實(shí)際工程應(yīng)用的工況,對(duì)制冷裝置的設(shè)計(jì)和改進(jìn)、系統(tǒng)中工質(zhì)組分的配比、混合物的優(yōu)化分離、合適的運(yùn)行溫度及壓力范圍等具有很好的指導(dǎo)及參考價(jià)值,以提高系統(tǒng)運(yùn)行效率。但由于還未見(jiàn)有相關(guān)實(shí)驗(yàn)數(shù)據(jù)的報(bào)道,本文的計(jì)算結(jié)果還無(wú)法進(jìn)行實(shí)驗(yàn)對(duì)比驗(yàn)證,只能作為一種預(yù)測(cè)。未來(lái)將會(huì)對(duì)R1234yf/R170/R14系混合物的氣液相平衡特性進(jìn)行實(shí)驗(yàn)研究,提供更為準(zhǔn)確的氣液相平衡數(shù)據(jù)。

    3 結(jié)論

    本文利用混合工質(zhì)氣液相平衡計(jì)算模型(PRWS-UNIFAC-PSRK)對(duì)溫度范圍283.15~323.15 K和壓力范圍0.44~1.74 MPa下二元體系R161/R1234yf,以及200~350 K溫度下三元體系R32/R125/R134a的氣液相平衡特性進(jìn)行了研究,并對(duì)比了計(jì)算結(jié)果與參考文獻(xiàn)的實(shí)驗(yàn)數(shù)據(jù)。二元混合物R161/R1234yf的壓力計(jì)算值與實(shí)驗(yàn)值相對(duì)誤差基本在±1.5%以內(nèi),氣相組分質(zhì)量分?jǐn)?shù)的絕對(duì)誤差均在±0.02內(nèi);而三元混合物R32/R125/R134a的壓力及質(zhì)量分?jǐn)?shù)偏差則分別在±4%和±0.02以內(nèi),說(shuō)明模型對(duì)二元混合工質(zhì)氣液相平衡特性的預(yù)測(cè)精度較高,滿足實(shí)際工程應(yīng)用的需求。

    利用文獻(xiàn)已有的相互作用系數(shù),對(duì)強(qiáng)非共沸工質(zhì)R1234yf/R170/R14系在0.2、0.5、1.5、2.2 MPa下的等壓氣液相平衡特性進(jìn)行模擬研究,并構(gòu)建了相關(guān)的三維立體相圖。混合物質(zhì)量分?jǐn)?shù)比為0.4/0.2/0.4時(shí)體系的溫度滑移現(xiàn)象最明顯,最大滑移溫度達(dá)72.5 K。隨著系統(tǒng)壓力的增大,泡點(diǎn)溫度和露點(diǎn)溫度也隨之上升,且R1234yf組分的質(zhì)量分?jǐn)?shù)越大,泡點(diǎn)溫度和露點(diǎn)溫度也越大。模擬結(jié)果對(duì)相關(guān)制冷系統(tǒng)設(shè)計(jì)及優(yōu)化有一定的參考價(jià)值。采用多參數(shù)狀態(tài)方程及改進(jìn)活度系數(shù)模型,可進(jìn)一步提高模型的預(yù)測(cè)精度。

    [1] 芮勝軍.帶兩路旁通的三級(jí)自動(dòng)復(fù)疊制冷系統(tǒng)實(shí)驗(yàn)研究[D].上海:上海理工大學(xué),2014. (RUI Shengjun. Experimental study on three stage cascade refrigeration system with two bypass[D].Shanghai: University of Shanghai for Science and Technology, 2014.)

    [2] RUHEMANN M. The separation of gases[M]. 2nd ed. London: Oxford University Press,1949.

    [3] KLEEMENKO A P. One-flow cascade cycle[C]//Proceedings of 10th International Congress of Refrigeration.Copenhagen, 1959.

    [4] FUDERER A. Compression process for refrigeration:US3203194[P]. 1965-08-31.

    [5] LITTLE W A, SAPOZHNIKOV I. Low cost cryocoolers for cryoelectronics[J]. Cryocooler, 1997, 9: 509-513.

    [6] NAYAK H G, VENKATARATHNAM G. Performance of an auto refrigerant cascade refrigerator operating in gas refrigerant supply (GRS) mode with nitrogen hydrocarbon and argon-hydrocarbon refrigerants[J]. Cryogenics, 2009, 49(7):350-359.

    [7] NARASIMHAN N L, VENKATARATHNAM G. Effect of mixture composition and hardware on the performance of a single stage JT refrigerator[J]. Cryogenics, 2011, 51(8): 446-451.

    [8] REDDY K R, MURTHY S S, VENKATARATHNAM G. Relationship between the cool down characteristics of J-T refrigerators and mixture composition[J]. Cryogenics, 2010, 50(6/7): 421-425.

    [9] 李文林.混合制冷劑的研究與發(fā)展[J].流體工程,1986,15(4):60-64. (LI Wenlin. Research and development of mixed refrigerants[J]. Fluid Machinery, 1986,15(4):60-64.)

    [10] 李文林,華小龍. 混合工質(zhì)制冷循環(huán)的試驗(yàn)研究[J]. 流體工程,1987,9(5):49-54. (LI Wenlin, HUA Xiaolong. Experimental study on refrigeration cycle of mixed refrigerant[J]. Fluid Machinery, 1987, 9(5):49-54.)

    [11] LUO Ercang, GONG Maoqiong, ZHOU Yuan, et al. Experimental investigation of a mixed refrigerant T-T cryocooler operating from 30 to 60 K[J]. Advanced in Cryogenic Engineering, 2000, 45:315-322.

    [12] WU J F, GONG M Q, LIU J L, et al. A new type mixture refrigeration auto-cascade cycle with partial condensation and separation reflux exchanger and its preliminary experimental test[J]. Advances in Cryogenic Engineering, 2002, 47:887-892.

    [13] GONG Maoqiong, WU Jianfeng, LUO Ercang, et al. Further development of the mixture refrigeration cycle with a dephlegmation separator[J]. Cryocoolers, 2003, 12:603-608.

    [14] GONG Maoqiong, WU Jianfeng, LUO Ercang, et al. Study of the single-stage mixed-gases refrigeration cycle for cooling temperature-distributed heat loads[J]. International Journal of Thermal Sciences, 2004, 43(1): 31-41.

    [15] 王國(guó)棟.非共沸制冷劑自復(fù)疊制冷系統(tǒng)特性及應(yīng)用研究[D]. 上海:上海理工大學(xué), 2004. (WANG Guodong. Study on the characteristics and application of non azeotropic refrigerant auto cascade refrigeration system[D]. Shanghai: University of Shanghai for Science and Technology, 2002.)

    [16] 陸向陽(yáng).三級(jí)自動(dòng)復(fù)疊制冷循環(huán)的實(shí)驗(yàn)研究[D]. 上海:上海理工大學(xué), 2002. (LU Xiangyang. Experimental study on three stage auto cascade refrigeration cycle[D]. Shanghai: University of Shanghai for Science and Technology, 2004.)

    [17] 荊磊. R600a/R23自動(dòng)復(fù)疊制冷系統(tǒng)及其應(yīng)用[D]. 上海:上海理工大學(xué), 2006. (JING Lei. Auto cascade refrigeration system and application of R600a/R23[D]. Shanghai: University of Shanghai for Science and Technology, 2006.)

    [18] 劉紅紹. 單級(jí)自動(dòng)復(fù)疊制冷系統(tǒng)的實(shí)驗(yàn)研究[D]. 上海:上海理工大學(xué), 2008. (LIU Hongshao. Experimental study on single stage auto cascade refrigeration system[D]. Shanghai: University of Shanghai for Science and Technology, 2008.)

    [19] 朱軍韜,張華. 混合工質(zhì)Linde-Hampson 制冷系統(tǒng)的實(shí)驗(yàn)研究[J].制冷學(xué)報(bào),2012, 33(2):18-23. (ZHU Juntao, ZHANG Hua. The experimental research on refrigerant mixtures of a Linde-Hampson refrigeration system[J]. Journal of Refrigeration, 2012, 33(2):18-23.)

    [20] HU Peng, CHEN Longxiang, ZHU Wanbao, et al. Isothermal VLE measurements for the binary mixture of 2,3,3,3-tetrafluoroprop-1-ene(HFO-1234yf) + 1,1-difluoroethane(HFC-152a)[J]. Fluid Phase Equilibria, 2014, 373:80-83.

    [21] HU Peng, CHEN Longxiang, CHEN Zeshao. Vapor-liquid equilibria for binary system of 2,3,3,3-tetrafluoroprop-1-ene (HFO-1234yf) + isobutane (HC-600a)[J]. Fluid Phase Equilibria, 2014, 365:1-4.

    [22] CHEN Longxiang, HU Peng, ZHU Wanbao, et al. Vapor-liquid equilibria of fluoroethane (HFC-161) + 2,3,3,3-tetrafluoroprop-1-ene (HFO-1234yf)[J]. Fluid Phase Equilibria, 2015, 392:19-23.

    [23] HU Peng, ZHU Wanbao, CHEN Longxiang, et al. Vapor-liquid equilibria measurements of 1,1,1,2-tetrafluoroethane (HFC-134a) + 2,3,3,3-tetrafluoroprop-1-ene (HFO-1234yf) + isobutane (HC-600a) ternary system[J]. Fluid Phase Equilibria, 2016, 414:111-116.

    [24] HU Peng, CHEN Longxiang, CHEN Zeshao. Vapor-liquid equilibria for the 1,1,1,2-tetrafluoroethane (HFC-134a) + 1,1,1,2,3,3,3-heptafluoropropane (HFC-227ea) and 1,1,1-trifluoroethane (HFC-143a) + 2,3,3,3-tetrafluoroprop-1-ene (HFO-1234yf) systems[J]. Fluid Phase Equilibria, 2013, 360:293-297.

    [25] 胡芃, 高敬軒, 陳龍祥, 等. PR方程結(jié)合vdW混合法則推算二元及三元HFC/HC混合工質(zhì)汽液相平衡性質(zhì)[J].化工學(xué)報(bào), 2012, 63(2):350-355. (HU Peng, GAO Jingxuan, CHEN Longxiang, et al. Prediction of vapor-liquid equilibria for binary and ternary HFC/HC mixtures using PR equation of state and vdW mixing rule[J]. CIESC Journal, 2012, 63(2):350-355.)

    [26] LIA J, TILLNER-ROTHB R, SATOA H, et al. Equation of state for hydrofluorocarbon refrigerant mixtures of HFC-125/143a, HFC-125/134a, HFC-134a/143a and HFC-125/134a/143a[J]. Fluid Phase Equilibria, 1999, 161(2):225-239.

    [27] GONG Maoqiong, ZHAO Yanxing, DONG Xueqiang. Measurements of isothermal (vapor + liquid) equilibria for the (propane + cis-1,3,3,3-tetrafluoropropene) system at temperatures from (253.150 to 293.150)K[J]. The Journal of Chemical Thermodynamics, 2016, 98:319-323.

    [28] GONG Maoqiong, GUO Hao, DONG Xueqiang. (Vapor + liquid) phase equilibria measurements for {trifluoroiodomethane (R13I1) + propane (R290)} from T=(258.150 to 283.150)K[J]. The Journal of Chemical Thermodynamics, 2014, 79:167-170.

    [29] LIM J S,PARK J Y, KANG J W, et al. Measurement of vapor-liquid equilibria for the binary systems of propane + 1,1,1,2-tetrafluoroethane and 1,1,1-trifluoroethane + propane at various temperatures[J]. Fluid Phase Equilibria, 2006, 243(1/2):57-63.

    [30] JU M, YUN Y,SHINM S, et al. (Vapour + liquid) equilibria of the {trifluoromethane (HFC-23) + propane} and {trifluoromethane (HFC-23) + n-butane} systems[J]. The Journal of Chemical Thermodynamics, 2009, 41(12):1339-1342.

    [31] 陳秀萍, 祁影霞, 趙勝喜, 等. 新型二元混合制冷劑(R290 + R227ea) 氣液相平衡研究[J]. 制冷學(xué)報(bào), 2014,35(5):94-100. (CHEN Xiuping, QI Yingxia, ZHAO Shengxi, et al. Vapor-liquid equilibrium properties of new binary mixture refrigerant (R290 + R227ea)[J]. Journal of Refrigeration, 2014, 35(5):94-100.)

    [32] BUDINSKY R, VACEK V, LISAL M. Vapor-liquid equilibria of alternative refrigerants and their binaries by molecular simulations employing the reaction Gibbs ensemble Monte Carlo method[J]. Fluid Phase Equilibria, 2004, (222/223):213-220.

    [33] HO Q N, YOO K S, LEE B G. Measurement of vapor-liquid equilibria for the binary mixture of propylene (R1270) + propane (R290)[J]. Fluid Phase Equilibria, 2006, 245(1):63-70.

    [34] HO Q N, LEE B G, PARK J Y, et al. Measurement of vapor-liquid equilibria for the binary mixture of propylene (R-1270) + 1,1,1,2-tetrafluoroethane (HFC-134a)[J]. Fluid Phase Equilibria, 2004, 225:125-132.

    [35] 公茂瓊, 吳劍鋒, 羅二倉(cāng). 深冷混合工質(zhì)節(jié)流制冷原理及應(yīng)用[M].北京:中國(guó)科學(xué)技術(shù)出版社, 2014. (GONG Maoqiong, WU Jianfeng, LUO Ercang. The principle and application of throttling refrigeration for cryogenic mixtures[M]. Beijing: China Science and Technology Press, 2014.)

    [36] PENG D Y,ROBINSON D B. A new two-constant equation of state[J]. Industrial & Engineering Chemistry Fundamentals, 1976, 15(1):59-64.

    [37] WONG D S H, SANDLER S. A theoretically correct mixing rule for cubic equations of state[J]. Aiche Journal, 1992, 38(5):671-680.

    [38] ORBEY H, SANDLER S I. Modelling vapour-liquid equilibria: cubic equation of state and their mixing rules[J]. Cambridge: Cambridge University Press, 1998.

    [39] ZHANG Haiyang, GONG Maoqiong, LI Huiya, et al. A simple model for temperature-independent kij of the PR-vdW model for mixtures containing HCs, HFCs, PFCs, HFOs, CO2, RE170 and R13I1[J]. Fluid Phase Equilibria, 2016, 425:374-384.

    [40] ZHAO Yanxing, GONG Maoqiong, DONG Xueqiang, et al. Prediction of ternary azeotropic refrigerants with a simple method[J]. Fluid Phase Equilibria, 2016, 425:72-83.

    [41] 吳獻(xiàn)忠, 崔曉鈺, 李美玲. UNIFAC基團(tuán)貢獻(xiàn)法預(yù)測(cè)混合制冷劑的氣液相平衡[J].化工學(xué)報(bào),2005, 56(10):1832-1836. (WU Xianzhong, CUI Xiaoyu, LI Meiling. Estimation of refrigerant mixture phase equilibria with UNIFAC model[J]. CIESC Journal, 2005, 56(10):1832-1836.)

    [42] MIYAUCHI H, YASUDA K, MATSUMOTO Y,et al. Isothermal phase equilibria for the (HFC-32 + HFC-134a) mixed-gas hydrate system[J]. Journal of Chemical Thermodynamics, 2012, 47:1-5.

    [43] NICOLA G D, GIULIANI G,PASSERINIG, et al. Vapor-Liquid-Equilibria (VLE) properties of R-32+R-134a system derived from isochoric measurements[J]. Fluid Phase Equilibria, 1998, 153(1):143-165.

    [44] SHIMAWAKI S, FUJII K, HIGASHI Y. Precise measurements of the vapor-liquid equilibria (VLE) of HFC-32/134a mixtures using a new apparatus[J]. International Journal of Thermophysics, 2002, 23(3):801-808.

    [45] KATO R, NISHIUMI H. Vapor-liquid equilibria and critical loci of binary and ternary systems composed of CH2F2, C2HF5and C2H2F4[J]. Fluid Phase Equilibria, 2006, 249:140-146.

    [46] NAGEL M, BIER K. Vapor-liquid equilibria of ternary mixtures of the refrigerants R32, R125and R134a[J]. International Journal of Refrigeration, 1995, 18(8):534-543.

    [47] HU Xiaozhen, YANG Tao, MENG Xianyang, et al. Vapor liquid equilibria measurements for difluoromethane(R32)+2,3,3,3-tetrafluoroprop-1-ene(R1234yf) and fluoroethane(R161) + 2,3,3,3-tetrafluoroprop-1-ene(R1234yf)[J]. Fluid Phase Equilibria, 2017, 438:10-17.

    猜你喜歡
    狀態(tài)方程制冷系統(tǒng)工質(zhì)
    海洋溫差能發(fā)電熱力循環(huán)系統(tǒng)的工質(zhì)優(yōu)選
    R290/R170單級(jí)壓縮回?zé)嶂评湎到y(tǒng)模擬研究
    LKP狀態(tài)方程在天然氣熱物性參數(shù)計(jì)算的應(yīng)用
    煤氣與熱力(2021年6期)2021-07-28 07:21:30
    采用R1234ze(E)/R245fa的非共沸混合工質(zhì)有機(jī)朗肯循環(huán)系統(tǒng)實(shí)驗(yàn)研究
    采用二元非共沸工質(zhì)的有機(jī)朗肯循環(huán)熱力學(xué)分析
    基于隨機(jī)與區(qū)間分析的狀態(tài)方程不確定性比較
    若干低GWP 純工質(zhì)在空調(diào)系統(tǒng)上的應(yīng)用分析
    奧迪A6L車制冷系統(tǒng)工作異常
    用狀態(tài)方程模擬氨基酸水溶液的熱力學(xué)性質(zhì)
    制冷系統(tǒng)節(jié)電控制改造的研究
    河南科技(2014年10期)2014-02-27 14:09:20
    欧美+亚洲+日韩+国产| 青草久久国产| svipshipincom国产片| 在线播放无遮挡| 999久久久精品免费观看国产| 在线观看一区二区三区| 91在线精品国自产拍蜜月 | 在线观看66精品国产| 此物有八面人人有两片| 俄罗斯特黄特色一大片| 欧美日韩一级在线毛片| 国产欧美日韩精品一区二区| 一进一出抽搐gif免费好疼| 国产蜜桃级精品一区二区三区| 国产高清有码在线观看视频| 亚洲精品456在线播放app | 99热这里只有是精品50| 三级男女做爰猛烈吃奶摸视频| 一a级毛片在线观看| 香蕉丝袜av| av片东京热男人的天堂| 成人永久免费在线观看视频| 国产真人三级小视频在线观看| 国产精品爽爽va在线观看网站| 成人性生交大片免费视频hd| 午夜免费男女啪啪视频观看 | 国产97色在线日韩免费| 亚洲va日本ⅴa欧美va伊人久久| 女人高潮潮喷娇喘18禁视频| 国产黄色小视频在线观看| 免费在线观看成人毛片| 一级毛片高清免费大全| 午夜免费成人在线视频| 99国产综合亚洲精品| 99热6这里只有精品| 国产av麻豆久久久久久久| 三级国产精品欧美在线观看| 国产一区二区激情短视频| 真人一进一出gif抽搐免费| 国产又黄又爽又无遮挡在线| 国内精品美女久久久久久| 国产精品久久久人人做人人爽| 久久久久久久精品吃奶| 天堂影院成人在线观看| 麻豆久久精品国产亚洲av| 美女高潮喷水抽搐中文字幕| 久9热在线精品视频| 精品一区二区三区av网在线观看| 国产高清视频在线播放一区| 露出奶头的视频| 亚洲人成伊人成综合网2020| 国产探花极品一区二区| 热99在线观看视频| 毛片女人毛片| 天堂av国产一区二区熟女人妻| 我的老师免费观看完整版| 天堂动漫精品| 18禁黄网站禁片午夜丰满| 看片在线看免费视频| 欧美一区二区精品小视频在线| 亚洲内射少妇av| tocl精华| 国产aⅴ精品一区二区三区波| 日本成人三级电影网站| 中文字幕人妻丝袜一区二区| 国产成人av教育| a级一级毛片免费在线观看| 国产亚洲精品一区二区www| 91字幕亚洲| 老汉色av国产亚洲站长工具| xxxwww97欧美| 又黄又爽又免费观看的视频| 午夜免费观看网址| 亚洲成人久久爱视频| 中文字幕av在线有码专区| 又粗又爽又猛毛片免费看| 成人特级av手机在线观看| 国产成人av教育| 国产色婷婷99| 精品电影一区二区在线| 我要搜黄色片| 欧美国产日韩亚洲一区| 中文亚洲av片在线观看爽| 国产精品,欧美在线| 国产精品久久视频播放| 成人欧美大片| 欧美+亚洲+日韩+国产| 听说在线观看完整版免费高清| 两人在一起打扑克的视频| 亚洲一区二区三区不卡视频| 久久国产乱子伦精品免费另类| 国产伦人伦偷精品视频| 午夜久久久久精精品| 亚洲电影在线观看av| 亚洲18禁久久av| 久9热在线精品视频| 看片在线看免费视频| 精品人妻1区二区| 国产精品女同一区二区软件 | 三级毛片av免费| 在线观看舔阴道视频| 精品久久久久久,| 亚洲av第一区精品v没综合| 亚洲久久久久久中文字幕| 欧美日韩综合久久久久久 | 国产精品久久久久久久电影 | netflix在线观看网站| 国产av在哪里看| 一个人免费在线观看电影| 欧美极品一区二区三区四区| 一区二区三区激情视频| 国产精品爽爽va在线观看网站| 夜夜夜夜夜久久久久| 天天一区二区日本电影三级| 美女高潮喷水抽搐中文字幕| 日本撒尿小便嘘嘘汇集6| 天天一区二区日本电影三级| 波多野结衣高清作品| 国产美女午夜福利| 岛国在线观看网站| 日本精品一区二区三区蜜桃| 桃色一区二区三区在线观看| 欧美日韩综合久久久久久 | 每晚都被弄得嗷嗷叫到高潮| 亚洲性夜色夜夜综合| 免费看日本二区| 激情在线观看视频在线高清| 久久精品国产清高在天天线| 亚洲成a人片在线一区二区| 天天一区二区日本电影三级| 老汉色av国产亚洲站长工具| 人妻久久中文字幕网| 国产精品亚洲一级av第二区| av在线天堂中文字幕| 人妻久久中文字幕网| 久久国产精品影院| 亚洲无线在线观看| 精品日产1卡2卡| 麻豆成人av在线观看| 中文字幕人妻丝袜一区二区| 亚洲欧美精品综合久久99| 三级男女做爰猛烈吃奶摸视频| 国产高清视频在线观看网站| 国产v大片淫在线免费观看| 久久国产精品影院| 久久久久久久午夜电影| www.熟女人妻精品国产| 精品电影一区二区在线| 国产精品久久视频播放| 91九色精品人成在线观看| 国产不卡一卡二| 午夜免费激情av| 夜夜躁狠狠躁天天躁| 91字幕亚洲| 国产高清有码在线观看视频| 美女被艹到高潮喷水动态| 村上凉子中文字幕在线| 一本一本综合久久| 欧美绝顶高潮抽搐喷水| 欧美黄色片欧美黄色片| 欧美中文日本在线观看视频| 亚洲一区二区三区色噜噜| 丰满乱子伦码专区| 国产精品 欧美亚洲| 精品一区二区三区人妻视频| 亚洲av第一区精品v没综合| 九色国产91popny在线| 日韩欧美国产一区二区入口| 性色avwww在线观看| 成人欧美大片| 日韩欧美一区二区三区在线观看| 少妇人妻精品综合一区二区 | 日本黄色视频三级网站网址| 天天添夜夜摸| xxx96com| 久久精品国产自在天天线| 国产精华一区二区三区| 两性午夜刺激爽爽歪歪视频在线观看| 五月玫瑰六月丁香| 欧美另类亚洲清纯唯美| 国产精品影院久久| 午夜免费激情av| 在线播放无遮挡| 国产免费av片在线观看野外av| 欧美日韩福利视频一区二区| 国产精品1区2区在线观看.| 床上黄色一级片| 女人被狂操c到高潮| 亚洲人成网站在线播放欧美日韩| 免费大片18禁| 啦啦啦观看免费观看视频高清| 男女午夜视频在线观看| 国产精品嫩草影院av在线观看 | 操出白浆在线播放| 小说图片视频综合网站| 国产黄色小视频在线观看| 两个人看的免费小视频| 五月伊人婷婷丁香| 无人区码免费观看不卡| 精品一区二区三区视频在线观看免费| 观看美女的网站| 国产 一区 欧美 日韩| 国产精品久久久久久亚洲av鲁大| 久久九九热精品免费| 国产伦一二天堂av在线观看| 18禁国产床啪视频网站| 91麻豆av在线| 日韩欧美 国产精品| 亚洲欧美一区二区三区黑人| 波多野结衣巨乳人妻| 国产精品电影一区二区三区| 99久久精品热视频| 又爽又黄无遮挡网站| 成年女人看的毛片在线观看| 久久亚洲真实| 欧美+日韩+精品| 在线免费观看的www视频| 亚洲成a人片在线一区二区| 欧美区成人在线视频| 欧美性猛交╳xxx乱大交人| 中文字幕人妻熟人妻熟丝袜美 | 99热这里只有精品一区| 午夜福利在线在线| 99在线视频只有这里精品首页| 伊人久久精品亚洲午夜| 90打野战视频偷拍视频| 精品一区二区三区av网在线观看| 久久精品国产亚洲av香蕉五月| 国产中年淑女户外野战色| 一个人看的www免费观看视频| 女生性感内裤真人,穿戴方法视频| 欧美另类亚洲清纯唯美| 久久人人精品亚洲av| 黄片大片在线免费观看| 老司机在亚洲福利影院| 亚洲av一区综合| 国产精品一区二区免费欧美| av欧美777| 一级黄片播放器| 亚洲成人免费电影在线观看| 亚洲av熟女| a级一级毛片免费在线观看| 亚洲美女视频黄频| 中出人妻视频一区二区| 午夜福利在线在线| 亚洲成人久久爱视频| 又紧又爽又黄一区二区| 两个人视频免费观看高清| 免费一级毛片在线播放高清视频| 美女高潮喷水抽搐中文字幕| 在线播放国产精品三级| 久久精品国产自在天天线| 黄色片一级片一级黄色片| 三级男女做爰猛烈吃奶摸视频| 极品教师在线免费播放| 天堂av国产一区二区熟女人妻| 国产精品久久久久久精品电影| 很黄的视频免费| 夜夜看夜夜爽夜夜摸| 两个人视频免费观看高清| 麻豆一二三区av精品| 香蕉av资源在线| 性欧美人与动物交配| 12—13女人毛片做爰片一| 最近最新免费中文字幕在线| 亚洲欧美日韩东京热| 免费高清视频大片| 国产 一区 欧美 日韩| 3wmmmm亚洲av在线观看| 性欧美人与动物交配| 老熟妇仑乱视频hdxx| 一区二区三区免费毛片| 少妇高潮的动态图| 国产黄色小视频在线观看| 天堂av国产一区二区熟女人妻| 久久久久久久久大av| 99久国产av精品| 亚洲av成人不卡在线观看播放网| 99精品久久久久人妻精品| 亚洲专区国产一区二区| 精品免费久久久久久久清纯| 成人永久免费在线观看视频| 亚洲狠狠婷婷综合久久图片| 亚洲av中文字字幕乱码综合| 老汉色av国产亚洲站长工具| 精品一区二区三区av网在线观看| 97超视频在线观看视频| 天堂影院成人在线观看| 国产精品久久久久久精品电影| 国产v大片淫在线免费观看| 日韩精品青青久久久久久| 禁无遮挡网站| 熟女电影av网| 中国美女看黄片| 欧美精品啪啪一区二区三区| 18禁美女被吸乳视频| 午夜日韩欧美国产| 国产精品日韩av在线免费观看| 久久伊人香网站| 美女免费视频网站| 国产v大片淫在线免费观看| 国产精品影院久久| 国内精品久久久久久久电影| 波多野结衣高清无吗| 精品一区二区三区视频在线 | 亚洲人与动物交配视频| 久久久久久人人人人人| 最近在线观看免费完整版| 小说图片视频综合网站| 一边摸一边抽搐一进一小说| 日本成人三级电影网站| 午夜老司机福利剧场| 欧美激情久久久久久爽电影| 欧美成人性av电影在线观看| 亚洲午夜理论影院| 国产国拍精品亚洲av在线观看 | 国产一区在线观看成人免费| 老司机福利观看| 亚洲自拍偷在线| 国产麻豆成人av免费视频| 国内精品美女久久久久久| 久久久久九九精品影院| 精品福利观看| 日本 欧美在线| 国产精品99久久久久久久久| 叶爱在线成人免费视频播放| 在线播放无遮挡| 国产精品日韩av在线免费观看| 人妻丰满熟妇av一区二区三区| 免费电影在线观看免费观看| 最新美女视频免费是黄的| 九九久久精品国产亚洲av麻豆| 成熟少妇高潮喷水视频| 国内揄拍国产精品人妻在线| 日本三级黄在线观看| 国产一区在线观看成人免费| 亚洲人成网站在线播| 丁香欧美五月| 亚洲第一欧美日韩一区二区三区| 美女免费视频网站| 亚洲av电影不卡..在线观看| 最新中文字幕久久久久| 熟女少妇亚洲综合色aaa.| 精品久久久久久,| 国产亚洲欧美98| 亚洲av成人不卡在线观看播放网| 日本撒尿小便嘘嘘汇集6| 男人舔奶头视频| 久久久色成人| 久久久久久人人人人人| 热99re8久久精品国产| 91久久精品电影网| 午夜激情福利司机影院| 久久久久亚洲av毛片大全| 亚洲人与动物交配视频| 午夜免费成人在线视频| 女人被狂操c到高潮| 亚洲乱码一区二区免费版| x7x7x7水蜜桃| 久久人人精品亚洲av| 国产亚洲欧美98| 人人妻人人澡欧美一区二区| 亚洲精品456在线播放app | 伊人久久大香线蕉亚洲五| 99精品久久久久人妻精品| 一区二区三区高清视频在线| 99精品久久久久人妻精品| 久久久久久人人人人人| 国产中年淑女户外野战色| 波多野结衣高清无吗| 国内精品久久久久久久电影| 久久人人精品亚洲av| 成人鲁丝片一二三区免费| 久久国产乱子伦精品免费另类| 亚洲一区高清亚洲精品| 啦啦啦免费观看视频1| 亚洲第一欧美日韩一区二区三区| 久久精品夜夜夜夜夜久久蜜豆| 亚洲男人的天堂狠狠| 久久精品国产自在天天线| 熟女少妇亚洲综合色aaa.| 久久精品夜夜夜夜夜久久蜜豆| 亚洲人成网站在线播| av欧美777| svipshipincom国产片| 久久久成人免费电影| 国产精品永久免费网站| 宅男免费午夜| 一级毛片女人18水好多| 久久久久久国产a免费观看| 亚洲中文日韩欧美视频| 国产精品久久视频播放| 欧美日韩综合久久久久久 | 精品乱码久久久久久99久播| 高潮久久久久久久久久久不卡| www.999成人在线观看| 欧美一级a爱片免费观看看| 亚洲av成人av| 日韩欧美在线乱码| 亚洲欧美一区二区三区黑人| 欧美最新免费一区二区三区 | 很黄的视频免费| 日本与韩国留学比较| 欧美中文日本在线观看视频| 国内揄拍国产精品人妻在线| 国产 一区 欧美 日韩| 日本精品一区二区三区蜜桃| 成人精品一区二区免费| 超碰av人人做人人爽久久 | 综合色av麻豆| 欧美+日韩+精品| 69av精品久久久久久| 国产精品99久久99久久久不卡| 性色avwww在线观看| 午夜福利在线观看免费完整高清在 | 丁香欧美五月| 国产精品永久免费网站| 波多野结衣高清作品| 免费一级毛片在线播放高清视频| 欧美一区二区亚洲| 国产精品国产高清国产av| 亚洲五月天丁香| 夜夜看夜夜爽夜夜摸| 国产淫片久久久久久久久 | 最新中文字幕久久久久| 久久国产精品人妻蜜桃| 女生性感内裤真人,穿戴方法视频| 9191精品国产免费久久| 狠狠狠狠99中文字幕| 成人国产综合亚洲| 国产亚洲精品久久久com| 国产免费一级a男人的天堂| 午夜激情福利司机影院| 3wmmmm亚洲av在线观看| 观看免费一级毛片| 丰满乱子伦码专区| 成人高潮视频无遮挡免费网站| 欧美在线一区亚洲| 成人特级黄色片久久久久久久| 好男人在线观看高清免费视频| 国产爱豆传媒在线观看| 久久精品国产99精品国产亚洲性色| 欧美一级a爱片免费观看看| 搡老岳熟女国产| 成人18禁在线播放| 国产成年人精品一区二区| 国产伦精品一区二区三区四那| 最近视频中文字幕2019在线8| 特级一级黄色大片| www.www免费av| 成人av一区二区三区在线看| 亚洲一区高清亚洲精品| www.色视频.com| 黑人欧美特级aaaaaa片| 成年女人永久免费观看视频| 国产探花极品一区二区| 十八禁人妻一区二区| 人妻夜夜爽99麻豆av| 亚洲人成伊人成综合网2020| 最新中文字幕久久久久| 人妻丰满熟妇av一区二区三区| 一个人免费在线观看的高清视频| 亚洲国产精品合色在线| 成年人黄色毛片网站| 日韩精品青青久久久久久| 中出人妻视频一区二区| 97超视频在线观看视频| 亚洲午夜理论影院| 亚洲精品影视一区二区三区av| 女同久久另类99精品国产91| 99精品久久久久人妻精品| 村上凉子中文字幕在线| 天堂影院成人在线观看| 久久久久久大精品| 又黄又爽又免费观看的视频| 欧美乱妇无乱码| av天堂在线播放| 亚洲av不卡在线观看| 久久久久久久久久黄片| 亚洲国产日韩欧美精品在线观看 | 每晚都被弄得嗷嗷叫到高潮| 欧美3d第一页| 久久中文看片网| 看黄色毛片网站| 欧美zozozo另类| 国产极品精品免费视频能看的| 亚洲人成伊人成综合网2020| 搡女人真爽免费视频火全软件 | 国产精品一及| 免费av观看视频| 69av精品久久久久久| 不卡一级毛片| 国产成人欧美在线观看| 美女黄网站色视频| 搡老熟女国产l中国老女人| 欧美成人a在线观看| 两个人的视频大全免费| 狂野欧美激情性xxxx| 国产精品,欧美在线| 日韩欧美在线乱码| 一区二区三区免费毛片| 男人舔奶头视频| 内地一区二区视频在线| www.色视频.com| 精品久久久久久久人妻蜜臀av| 最近最新中文字幕大全免费视频| 色哟哟哟哟哟哟| 国产成人欧美在线观看| 88av欧美| 国产免费男女视频| 青草久久国产| 亚洲成人精品中文字幕电影| xxx96com| 国产伦人伦偷精品视频| 欧美日韩黄片免| 欧美丝袜亚洲另类 | 丁香欧美五月| 高潮久久久久久久久久久不卡| 免费大片18禁| 国产午夜福利久久久久久| 美女高潮的动态| 亚洲,欧美精品.| 久久精品国产综合久久久| 最近最新免费中文字幕在线| 日韩精品青青久久久久久| 国产成人啪精品午夜网站| 床上黄色一级片| 欧美激情久久久久久爽电影| 成人性生交大片免费视频hd| 日韩欧美精品v在线| 校园春色视频在线观看| 国产国拍精品亚洲av在线观看 | 深夜精品福利| 国产免费男女视频| 亚洲最大成人手机在线| av中文乱码字幕在线| 丰满人妻熟妇乱又伦精品不卡| 久久久久久大精品| 日本成人三级电影网站| 成年免费大片在线观看| 99热只有精品国产| 特级一级黄色大片| 国产欧美日韩一区二区精品| 国产不卡一卡二| 日日摸夜夜添夜夜添小说| 一个人观看的视频www高清免费观看| 12—13女人毛片做爰片一| 欧美极品一区二区三区四区| 国产三级黄色录像| 国产精品99久久99久久久不卡| 热99re8久久精品国产| 色综合亚洲欧美另类图片| 在线观看舔阴道视频| 男插女下体视频免费在线播放| 日韩欧美精品v在线| 亚洲欧美日韩高清在线视频| 欧美zozozo另类| 日本免费一区二区三区高清不卡| 精品久久久久久成人av| 淫秽高清视频在线观看| 在线a可以看的网站| 国产亚洲欧美98| 中文字幕人妻丝袜一区二区| 色老头精品视频在线观看| 国产美女午夜福利| 欧美午夜高清在线| 搡老熟女国产l中国老女人| 国产亚洲精品久久久com| 亚洲成人中文字幕在线播放| 丁香六月欧美| 一区二区三区免费毛片| 国产野战对白在线观看| 精品一区二区三区av网在线观看| 女人十人毛片免费观看3o分钟| 在线观看av片永久免费下载| 亚洲不卡免费看| 噜噜噜噜噜久久久久久91| 亚洲人成网站在线播放欧美日韩| www国产在线视频色| 有码 亚洲区| 麻豆成人av在线观看| 欧美3d第一页| 国产黄a三级三级三级人| 波野结衣二区三区在线 | 日韩人妻高清精品专区| 天堂影院成人在线观看| 国产黄色小视频在线观看| 国产精品综合久久久久久久免费| 亚洲第一电影网av| 在线观看美女被高潮喷水网站 | 亚洲成人中文字幕在线播放| 丁香六月欧美| a在线观看视频网站| 日韩中文字幕欧美一区二区| 婷婷丁香在线五月| 淫秽高清视频在线观看| 国产午夜福利久久久久久| 精品久久久久久成人av| 欧美性猛交╳xxx乱大交人| 精品一区二区三区人妻视频| 成熟少妇高潮喷水视频| 少妇的逼好多水| 在线观看免费视频日本深夜| 99热6这里只有精品| 久久精品国产清高在天天线| 亚洲最大成人中文| 国产成人aa在线观看| 久久精品国产自在天天线| 在线播放无遮挡| 好男人在线观看高清免费视频| 久久婷婷人人爽人人干人人爱| 国产午夜精品论理片| 一本一本综合久久| 美女免费视频网站| 色综合婷婷激情|