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

    [bmim]Zn2Cl5/NH3溶液的比熱容和過量焓

    2016-05-30 03:29:18陳偉梁世強(qiáng)郭永獻(xiàn)
    山東科學(xué) 2016年2期
    關(guān)鍵詞:離子液體比熱容

    陳偉,梁世強(qiáng),郭永獻(xiàn)

    (1.青島科技大學(xué)機(jī)電工程學(xué)院,山東 青島 266061; 2.中國科學(xué)院工程熱物理研究所,北京 100190)

    ?

    【能源與動(dòng)力】

    [bmim]Zn2Cl5/NH3溶液的比熱容和過量焓

    陳偉1,梁世強(qiáng)2*,郭永獻(xiàn)2

    (1.青島科技大學(xué)機(jī)電工程學(xué)院,山東 青島 266061; 2.中國科學(xué)院工程熱物理研究所,北京 100190)

    摘要:測定了離子液體[bmim]Zn2Cl5在T=(323.15~1 173.15) K 范圍內(nèi)的熱重曲線,結(jié)果顯示[bmim]Zn2Cl5 在T<637.15 K時(shí)具有很高的熱穩(wěn)定性。通過DSC測試得到[bmim]Zn2Cl5的比熱容數(shù)據(jù),在T= (251.15~383.15) K 范圍內(nèi)可以用一個(gè)圓錐曲線很好地?cái)M合。 實(shí)驗(yàn)測定了[bmim]Zn2Cl5(2)+NH3(1)二元體系溶液的摩爾過量焓,其中氨的摩爾分?jǐn)?shù)x1=(0.60~0.95),溫度值為T=288.15 K,303.15 K,318.15 K,333.15 K。采用NRTL模型對過量焓數(shù)據(jù)進(jìn)行擬合,得出二元可調(diào)參數(shù)和非隨機(jī)參數(shù)。過量焓數(shù)據(jù)的測量誤差和最大擬合偏差分別小于 4.8% 和 4.3%。在[bmim]Zn2Cl5比熱容和[bmim]Zn2Cl5/NH3過量焓數(shù)據(jù)的基礎(chǔ)上,計(jì)算了氨質(zhì)量分?jǐn)?shù)w1 = (0 ~ 1)、溫度范圍T= (273.15 ~343.15) K條件下[bmim]Zn2Cl5/NH3溶液的焓,所得焓濃圖對于研究[bmim]Zn2Cl5/NH3吸收式制冷系統(tǒng)性能至關(guān)重要。

    關(guān)鍵詞:比熱容;過量焓;離子液體;[bmim]Zn2Cl5/NH3;NRTL模型

    過去十年中,離子液體型吸收式制冷工質(zhì)對的研究頗受關(guān)注[1],原因在于傳統(tǒng)的工質(zhì)對的某些固有缺陷,如腐蝕、結(jié)晶及毒性等問題[2],可以使用離子液體作為吸收劑解決[3]。

    迄今為止,離子液體型工質(zhì)對的相關(guān)研究已有很多,如美國杜邦公司的Yokozeki等[4],韓國的Kim等[5],西班牙的Martin等[6],國內(nèi)的Chen(陳偉)等[7]、Liang(梁世強(qiáng))等[8]和Zhang(張曉東)等[9]分別研究了各種離子液體應(yīng)用于吸收式制冷系統(tǒng)的可能性。但是能夠在吸收式制冷循環(huán)特性方面超越傳統(tǒng)溴化鋰/水的離子液體型工質(zhì)對一直鮮有報(bào)道。制冷劑本身的性質(zhì)仍然是制約離子液體型工質(zhì)對性能的重要因素。文獻(xiàn)統(tǒng)計(jì)表明,以H2O為制冷劑的離子液體型工質(zhì)對普遍具有相對較好的理論循環(huán)效率,然而,同樣具有良好熱力性質(zhì)的NH3與離子液體組成的很多工質(zhì)對的理論循環(huán)性能卻并不理想[10],原因在于多數(shù)離子液體對NH3的吸收能力有限,導(dǎo)致吸收式制冷系統(tǒng)的高循環(huán)倍率和循環(huán)性能的惡化。

    Chen(陳偉)等[7,11-12]提出采用兩種含金屬離子的離子液體[bmim]Zn2Cl5和[Cnmim]2NiCl4作為NH3的吸收劑,通過測定[bmim]Zn2Cl5/NH3溶液的氣液相平衡性質(zhì),發(fā)現(xiàn)NH3在[bmim]Zn2Cl5中的溶解度要遠(yuǎn)遠(yuǎn)高于常規(guī)的離子液體[11],并且[bmim]Zn2Cl5/NH3的吸收式制冷理論循環(huán)特性非常優(yōu)異[7]。[bmim]Zn2Cl5/NH3溶液的一些熱力學(xué)性質(zhì),如比熱容、過量焓等都是考察該工質(zhì)對應(yīng)用潛力必不可少的基礎(chǔ)物性,本文通過實(shí)驗(yàn)測定[bmim]Zn2Cl5比熱容和[bmim]Zn2Cl5/NH3溶液的過量焓,并分別用絕對溫度的二次函數(shù)和NRTL[13-15]模型進(jìn)行了關(guān)聯(lián)。

    1實(shí)驗(yàn)材料與裝置

    1.1實(shí)驗(yàn)材料

    高純度無水液氨(純度≥99.99%, CAS No. 7664-41-7),購自北京市威特化工有限公司;高純度無水氯化鋅(純度≥99.99%, CAS No. 231-592-0) 購自北京市隆?;び邢薰荆桓呒兌入x子液體[bmim]Cl(純度≥99.9%,C8H15ClN2,CAS No. 231-592-0)購自上海成捷化學(xué)有限公司。將離子液體[bmim]Cl與無水氯化鋅以摩爾比1:2混合,加熱至200 ℃、攪拌均勻并維持48 h即可制得實(shí)驗(yàn)所需的離子液體[bmim]Zn2Cl5。

    1.2實(shí)驗(yàn)裝置和實(shí)驗(yàn)方法

    采用TGA/SDT Q600型熱重分析儀在溫度區(qū)間323.15 K~1 173.15 K對[bmim]Zn2Cl5進(jìn)行熱重分析。

    采用DSC 910S型差式掃描量熱儀在溫度區(qū)間210.15 K~381.15 K測量[bmim]Zn2Cl5的比熱容,操作工況為溫升速度5 K·min-1、氮?dú)怏w積流量40 cm3·min-1。

    采用如圖1所示的實(shí)驗(yàn)裝置在溫度為288.15 K、303.15 K、318.15 K以及333.15 K條件下測定了不同摩爾分?jǐn)?shù)[bmim]Zn2Cl5/NH3溶液過量焓。

    整個(gè)過量焓測量裝置由3個(gè)子系統(tǒng)構(gòu)成:溫度控制系統(tǒng)、[bmim]Zn2Cl5/NH3反應(yīng)器和數(shù)據(jù)采集儀。

    溫度控制系統(tǒng)包括氧化鋁保溫層、溫度傳感器(精度:0.05%,量程:0~450 ℃)、加熱器、風(fēng)扇葉片、風(fēng)扇電機(jī)、絕熱箱、PID溫控器、變壓器、可控電路開關(guān)和鐵絲網(wǎng)。根據(jù)溫度傳感器采集的溫度信號,PID控制器將會調(diào)節(jié)變壓器的輸出電壓以及風(fēng)扇電機(jī)的啟停,從而控制加熱功率和風(fēng)扇轉(zhuǎn)動(dòng),并最終將絕熱箱中的溫度控制在設(shè)定值,系統(tǒng)的控制誤差低于±0.02 K。

    [bmim]Zn2Cl5/NH3反應(yīng)器包括:水浴容器(體積:4 879.6±0.1 mL,質(zhì)量:3 925.5±0.1 g)、液氯鋼瓶(體積:79.64±0.02 mL,質(zhì)量:213.6±0.1 g)、自動(dòng)閥門、高壓反應(yīng)器(體積:216.25±0.05 mL,質(zhì)量351.7±0.1 g)。

    1 NH3鋼瓶; 2 自動(dòng)閥門; 3 高壓反應(yīng)容器; 4 水浴容器; 5 密封墊圈; 6 氧化鋁保溫層 ; 7 Pt100熱電阻; 8 K型熱電偶; 9 溫度傳感器; 10 加熱器; 11 風(fēng)扇葉片; 12 風(fēng)扇電機(jī); 13 絕熱箱; 14 計(jì)算機(jī); 15 數(shù)據(jù)采集系統(tǒng); 16 PID溫控器; 17 變壓器; 18 可控電路開關(guān); 19 鐵絲網(wǎng)圖1 溶液過量焓測量裝置Fig.1 Schematic diagram of the measurement apparatus of vapor excess enthalpies

    數(shù)據(jù)采集系統(tǒng)包括:Pt100熱電阻(精度:0.01%,量程:0~250 ℃)、K型熱電偶(精度:0.02%,量程:0~200 ℃)、安捷倫34970A數(shù)據(jù)采集儀、計(jì)算機(jī)。通過安捷倫采集儀,實(shí)時(shí)觀測壓力容器內(nèi)部溫度、壓力隨時(shí)間的變化曲線,判斷氣液相平衡狀態(tài)。

    實(shí)驗(yàn)流程如下:首先將離子液體裝入高壓反應(yīng)容器中,并用電子秤(精度:±0.001 g)稱量裝入離子液體的質(zhì)量。然后在自動(dòng)閥門開啟的情況下,對NH3鋼瓶抽真空,同時(shí)高壓反應(yīng)容器中也被抽成真空狀態(tài),關(guān)閉自動(dòng)閥門。再將液氨灌入NH3鋼瓶中,灌入液氨的質(zhì)量有一個(gè)精度為0.5%的質(zhì)量流量計(jì)確定。將灌裝好的[bmim]Zn2Cl5/NH3反應(yīng)系統(tǒng)置于水浴中,開啟溫控系統(tǒng),待到系統(tǒng)達(dá)到穩(wěn)態(tài)時(shí),開啟自動(dòng)閥門,使[bmim]Zn2Cl5和NH3在高壓反應(yīng)容器中發(fā)生反應(yīng)。最后,通過數(shù)據(jù)采集系統(tǒng)記錄實(shí)驗(yàn)數(shù)據(jù),一直到整個(gè)系統(tǒng)再次達(dá)到平衡狀態(tài)為止。

    其中,根據(jù)高壓反應(yīng)容器和NH3鋼瓶的體積和溫度,就可以通過氨氣的狀態(tài)方程計(jì)算出反應(yīng)容器中以氣相形式存在的NH3的質(zhì)量,從而確定被吸收的NH3的質(zhì)量和[bmim]Zn2Cl5/NH3溶液的組分。根據(jù)高壓反應(yīng)容器、NH3鋼瓶、水和水浴容器的熱容及溫度變化,即可求出吸收反應(yīng)的放熱量,從而確定 [bmim]Zn2Cl5/NH3溶液的過量焓。

    1.3誤差分析

    隨機(jī)誤差、測試裝置系統(tǒng)誤差以及UNIFAC模型預(yù)測氣液相平衡性質(zhì)[11]是引起過量焓測量誤差的主要因素,溫度、氨摩爾分?jǐn)?shù)和摩爾過量焓等參數(shù)的誤差估算結(jié)果列于表1中,其中,u表示誤差,T為絕對溫度,X為氨摩爾分?jǐn)?shù),HE為溶液摩爾過量焓。

    表1 實(shí)驗(yàn)誤差

    圖2 [bmim]Zn2Cl5的熱重曲線圖Fig.2 TG scan curve of [bmim]Zn2Cl5

    2結(jié)果及分析

    2.1熱重分析

    圖2是[bmim]Zn2Cl5在323.15 K~1173.15 K溫度區(qū)間內(nèi)、在氮?dú)鈿夥罩械臒嶂厍€。溫度為676.15 K時(shí),被測樣品的失重率為2.5%;溫度為774.15 K時(shí),被測樣品的失重率為40%,大致與離子液體中咪唑離子的質(zhì)量分?jǐn)?shù)相符,以上兩個(gè)溫度都是咪唑基離子液體比較典型的揮發(fā)溫度。繼續(xù)加熱樣品至1 043.15 K,其間樣品的質(zhì)量不斷下降并逐漸達(dá)到穩(wěn)定值,最終的失重率為63.6%。結(jié)果表明,在溫度低于676.15 K時(shí),[bmim]Zn2Cl5具有很高的熱穩(wěn)定性。

    2.2差式掃描量熱分析

    圖3是[bmim]Zn2Cl5樣品的比加熱熱流qm和比熱容cp隨溫度T的變化曲線??梢姡燃訜釤崃麟S著溫度的升高而升高,在溫度低于243.15 K時(shí),比加熱熱流的升高速度也是不斷增大的,但是當(dāng)溫度高于243.15 K時(shí),比加熱熱流隨著溫度呈近似線性地升高。在溫度低于243.15 K時(shí),[bmim]Zn2Cl5比熱容隨著溫度的升高而升高,且升高的速度隨著溫度的升高而明顯增大。當(dāng)溫度高于243.15 K時(shí),比熱容先是有一個(gè)較快的下降,然后又緩慢下降且趨于平緩,并在293.15 K左右達(dá)到最低點(diǎn),隨后又隨著溫度的升高而升高,升高的速度隨著溫度的升高亦顯出升高的趨勢。研究發(fā)現(xiàn),如圖4所示,溫度高于251.15 K的比熱容曲線可擬合為溫度的二次函數(shù),表達(dá)式如下:

    cp=2.393 27-0.006 91·T+0.000 011 767·T2

    (1)

    圖3 [bmim]Zn2Cl5樣品的比加熱熱流qm和比熱容cp隨溫度變化曲線Fig.3 DSC scan curve of the variations of heat flow qm and specific heat capacity cp of [bmim]Zn2Cl5 samples with temperature

    圖4 [bmim]Zn2Cl5的比熱容溫度擬合曲線Fig.4 cp-T diagram of [bmim]Zn2Cl5

    2.3溶液過量焓

    根據(jù)Renon等[16]提出的非隨機(jī)的二元液相(NRTL)模型,二元溶液的過量吉布斯自由能GE描述如下:

    (2)

    其中,R為理想氣體常數(shù),x1、x2分別為組分1和組分2的摩爾分?jǐn)?shù),τ12、τ21、G12和G21的表達(dá)式如下:

    (3)

    G12=exp(-ατ12),G21=exp(-ατ21),

    (4)

    其中,gij、gjj分別是組分對ij和jj的相互作用能,α為非隨機(jī)參數(shù)。

    過量吉布斯自由能與組分活度系數(shù)γ的關(guān)系如下:

    (5)

    其中,P為氣壓,i、j和k為組分標(biāo)號。

    因此,各溶液組分的活度系數(shù)表示如下:

    (6)

    (7)

    對[bmim]Zn2Cl5(1)/NH3(2)溶液而言,離子液體[bmim]Zn2Cl5具有不可測量的飽和蒸氣壓,因此,氣相成分即為單一的氨氣,NH3的活度系數(shù)γ2可以計(jì)算如下[17]:

    (8)

    其中,x2表示的是溶液的氨摩爾分?jǐn)?shù),γ2表示的是溶液中NH3的活度系數(shù),V2L表示的是氨的摩爾體積,B22表示的是NH3狀態(tài)方程的第二維里數(shù),P2既是氨的分壓也是總壓,P2S表示的是純氨的蒸氣壓力,當(dāng)溶液的溫度低于NH3的臨界溫度時(shí)[18],

    (9)

    其中,Tr表示的是相對溫度,即當(dāng)前溫度與臨界溫度的比值,當(dāng)溶液溫度高于NH3的臨界溫度時(shí),P2S被定義為NH3在當(dāng)前溫度和臨界體積狀態(tài)時(shí)的氣壓,可以用RK EOS狀態(tài)方程來計(jì)算:

    (10)

    (11)

    b=0.086 64·RTC/PC,

    (12)

    其中,基于溫度變化的溫度項(xiàng)α(T)可表示為:

    (13)

    NH3狀態(tài)方程常數(shù)和臨界參數(shù)TC、VC、PC以及擬合參數(shù)βk見表2。

    表2 NH3狀態(tài)方程常數(shù)和臨界參數(shù)

    [bmim]Zn2Cl5(1)/NH3(2)溶液的NRTL模型組分對ij和jj的相互作用能的定義如下:

    g12-g22=A1+B1T,

    (14)

    g21-g11=A2+B2T,

    (15)

    其中,A1、A2、B1、B2是用于關(guān)聯(lián)擬合的可調(diào)整參數(shù)。

    過量焓HE和過量吉布斯自由能之間還存在著如下的關(guān)系:

    (16)

    因此,過量焓可表示為溶液組分和溫度的函數(shù):

    。 (17)

    采用上述NRTL模型對在溫度為288.15 K、303.15 K、318.15 K以及333.15 K條件下測得的不同摩爾分?jǐn)?shù)的[bmim]Zn2Cl5/NH3溶液過量焓實(shí)驗(yàn)數(shù)據(jù)(見圖5)進(jìn)行關(guān)聯(lián),得出NRTL模型可調(diào)整參數(shù)A1、A2、B1、B2和非隨機(jī)參數(shù)α的關(guān)聯(lián)結(jié)果,列于表3,擬合的相關(guān)系數(shù)高達(dá)0.998。

    圖5 [bmim]Zn2Cl5/NH3溶液T-x-HE圖Fig.5 T-x-HE diagram of [bmim]Zn2Cl5/NH3 solution

    圖6 關(guān)聯(lián)的相對誤差和絕對誤差Fig.6 Absolute and relative deviations between experimental data and NRTL based calculation values

    圖5是[bmim]Zn2Cl5/NH3溶液的T-x-HE圖,符號標(biāo)記表示實(shí)驗(yàn)數(shù)據(jù),實(shí)線表示NRTL模型計(jì)算結(jié)果??梢?,在相同的氨摩爾分?jǐn)?shù)條件下,摩爾過量焓隨著溫度的升高而升高;在溫度一定的條件下,摩爾過量焓隨著氨摩爾分?jǐn)?shù)的升高呈現(xiàn)出一個(gè)先降后升的趨勢。因此,在一定的溫度條件下,都存在著一個(gè)特定的氨摩爾分?jǐn)?shù)使得摩爾過量焓取得最小值。在288.15 K、303.15 K、318.15 K、333.15 K溫度條件下,這個(gè)使得[bmim]Zn2Cl5/NH3溶液摩爾過量焓取得最小值的特定氨摩爾分?jǐn)?shù)分別為0.772、0.774、0.776和0.777,對應(yīng)的最小摩爾過量焓分別為-6 555.7 J·mol-1、-6 707.1 J·mol-1、-6 846.3 J·mol-1和-6 974.7 J·mol-1。

    圖6是NRTL模型關(guān)聯(lián)結(jié)果和實(shí)驗(yàn)測量結(jié)果之間的相對誤差和絕對誤差,絕對誤差被定義為關(guān)聯(lián)結(jié)果減去實(shí)驗(yàn)結(jié)果,相對誤差等于絕對誤差除以實(shí)驗(yàn)值。關(guān)聯(lián)的絕對誤差最大不超過0.15 kJ·mol-1,相對誤差不超過4.0%。該偏差主要有以下幾方面的因素引起:高壓反應(yīng)容器的體積(0.5%)、NH3鋼瓶的體積(0.5%)、[bmim]Zn2Cl5/NH3溶液的體積 (0.2%)、[bmim]Zn2Cl5的質(zhì)量(0.01%)、NH3的質(zhì)量(0.05%)、水的質(zhì)量(0.01%)、水的溫度分布不均勻性(1.4%)、水浴容器溫度分布不均勻性(1.2%)以及UNIFAC模型預(yù)測組分計(jì)算誤差(0.9%)?;谝陨险`差因素可估算出整個(gè)測量和關(guān)聯(lián)總的偏差將在±4.8%以內(nèi),與圖中給出的誤差范圍基本吻合。

    圖7 溫度T= (273.15~343.15) K、x2= (0~1)區(qū)間的[bmim]Zn2Cl5/NH3溶液摩爾過量焓計(jì)算值Fig.7Calculation values for mole excess enthalpies of [bmim]Zn2Cl5/NH3 solution at x2=(0~1) and T=(273.15~343.15) K

    圖7表明[bmim]Zn2Cl5/NH3溶液的摩爾過量焓均為負(fù)值,隨著氨摩爾分?jǐn)?shù)的升高,溶液的摩爾過量焓呈現(xiàn)出先降后升的趨勢,摩爾過量焓的最低值出現(xiàn)在氨摩爾分?jǐn)?shù)為0.78左右;隨著溫度的升高,摩爾過量焓呈現(xiàn)出近似線性上升的趨勢,但是摩爾過量焓的最低值出現(xiàn)在氨摩爾分?jǐn)?shù)基本保持不變時(shí)。

    3結(jié)論

    通過對離子液體[bmim]Zn2Cl5的熱穩(wěn)定性和比熱容、[bmim]Zn2Cl5/NH3溶液的摩爾過量焓的實(shí)驗(yàn)研究,獲得了分析計(jì)算[bmim]Zn2Cl5/NH3工質(zhì)對吸收式制冷循環(huán)特性所必需的基礎(chǔ)數(shù)據(jù)和數(shù)學(xué)模型,并得到如下結(jié)論:

    (1)在溫度低于676.15 K時(shí),[bmim]Zn2Cl5具有很高的熱穩(wěn)定性,適用溫度范圍比傳統(tǒng)工質(zhì)對更寬;

    (2)溫度高于251.15 K的[bmim]Zn2Cl5比熱容曲線可擬合為溫度的二次函數(shù),擬合精度較高;[bmim]Zn2Cl5比熱容比溴化鋰溶液小很多,適合用于開發(fā)吸收式制冷工質(zhì)對;

    (3)[bmim]Zn2Cl5/NH3溶液的摩爾過量焓可以用NRTL模型進(jìn)行關(guān)聯(lián),精度滿足應(yīng)用需要;[bmim]Zn2Cl5/NH3溶液的摩爾過量焓在應(yīng)用范圍內(nèi)均為負(fù)值,符合吸收式制冷工質(zhì)對特征。

    參考文獻(xiàn):

    [1]JAWAHAR C P, SARAVANAN R. Experimental studies on air-cooled NH3-H2O based modified gax absorption cooling system [J], Int J Refrigeration, 2011, 34(3): 658-666.

    [2]EZZINE B N, GARMA R, BELLAGI A. A numerical investigation of a diffusion-absorption refrigeration cycle based on R124-DMAC mixture for solar cooling [J], Energy, 2010, 35(5): 1874-1883.

    [3]ZHANG X D, HU D P. Performance analysis of the single-stage absorption heat transformer using a new working pair composed of ionic liquid and water [J], Appl Therm Eng, 2012, 37(1): 129-135.

    [4]YOKOZEKI A, SHIFLETT M B. Water solubility in ionic liquids and application to absorption cycles [J], Ind Eng Chem Res, 2010, 49 (10): 9496-9503.

    [5]KIM Y J, KIM S, JOSHI Y K, et al. Thermodynamic analysis of an absorption refrigeration system with ionic-liquid/refrigerant mixture as a working fluid [J]. Energy, 2012, 44(1): 1005-1016.

    [7]CHEN W, LIANG S Q, GUO Y X, et al. Thermodynamic analysis of absorption system using [bmim]Zn2Cl5/NH3as working pairs [J]. Energy Conversion and Management, 2014, 85: 13-19.

    [8]LIANG S Q, CHEN W, CHENG K Y, et al. The latent application of ionic liquids in absorption refrigeration [M]//Scott Handy. Application of ionic liquid in science and technology. Croatia: InTech, 2011: 467-494.

    [9]ZHANG X D, HU D P. Performance simulation of the absorption chiller using water and ionic liquid 1-ethyl-3-methylimidazolium dimethylphosphate as the working pair [J]. Applied Thermal Engineering, 2011, 31 (16): 3316-3321.

    [10]YOKOZEKI A, SHIFLETT M B. Vapor-liquid equilibria of ammonia + ionic liquid mixtures [J]. Applied Energy, 2007, 84(12): 1258-1273.

    [11]CHEN W, LIANG S Q, GUO Y X, et al. Investigation on vapor-iquid equilibria for binary systems of metal ion-containing ionic liquid [bmim]Zn2Cl5/NH3by experiment and modified UNIFAC model[J]. Fluid Phase Equilibria, 2013, 360: 1-6.

    [12]LIANG S Q, CHEN W, GUO Y X, et al. Ionic Liquids Facilitate the Development of Absorption Refrigeration[M]// Scott Handy. Ionic Liquids- Current State of the Art. Croatia: InTech, 2015: 623-654.

    [13]ZHANG Y, QUE H L, CHEN C C. Thermodynamic modeling for CO2absorption in aqueous MEA solution with electrolyte NRTL model [J]. Fluid Phase Equilib, 2011,311: 67-75.

    [14]RAYER A V, HENNI A, TONTIWACHWUTHIKUL P. High-pressure solubility of methane (CH4) and ethane (C2H6) in mixed polyethylene glycol dimethyl ethers (Genosorb 1753) and its selectivity in natural gas sweetening operations [J]. Journal of Chemical Engineering Data, 2012, 57(3): 764-775.

    [15]QUE H L, CHEN C C.Thermodynamic modeling of the NH3-CO2-H2O system with electrolyte NRTL model [J]. Industrial & Engineering Chemistry Research, 2011, 50(19): 11406-11421.

    [16]RENON H, PRAUSNITZ J M. Local compositions in thermodynamic excess functions for liquid mixtures [J].AIChE Journal, 1968, 14(1):135-144.

    [17]ISLAM A W, RAHMAN M H. A review of Barker’s activity coefficient method and VLE data reduction [J]. Journal of Chemical Thermodynamics, 2012, 44(1): 31-37.

    [18]李斯特, 程灶亮, 馬永其, 等. 氨工質(zhì)的新簡明狀態(tài)方程 [J]. 工程熱物理學(xué)報(bào), 2000, 21(1): 17-19.

    Heat capacity and excess enthalpy of [bmim]Zn2Cl5/NH3binary system solution

    CHEN Wei1,LIANG Shi-qiang2*,GUO Yong-xian2

    (1. School of Electromechanical Engineering, Qingdao University of Science and Technology, Qingdao 266061, China;2. Institute of Engineering Thermophysics, Chinese Academy of Sciences, Beijing 100190, China)

    Abstract∶We measured thermogravimetric curve for ionic liquid [bmim]Zn2Cl5 at T= (323.15~1 173.15) K. Results show that [bmim]Zn2Cl5 has higher thermal stability at T < 637.15 K. We acquired heat capacity data of [bmim]Zn2Cl5 through DSC scan, which could be well fitted with a conic curve at T= (251.15~383.15) K. We measured mole excess enthalpy of binary system [bmim]Zn2Cl5 (2) + NH3 (1) for ammonia molar fractions x1 = (0.60 ~ 0.95) at T=288.15 K,303.15 K,318.15 K,333.15 K.We employed NRTL model to fit excess enthalpy data and acquired binary adjustable parameters and non-random parameters. Measurement error and maximum fitting deviation of excess enthalpy data are respectively less than 4.8% and 4.3%. Based on heat capacity of [bmim]Zn2Cl5 and excess enthalpy data of [bmim]Zn2Cl5/NH3, we also calculated the enthalpy of [bmim]Zn2Cl5/NH3 solution for ammonia mass fractions of w1 = (0 ~1) and temperature scope T=(273.15 ~ 343.15) K. The obtained enthalpy-concentration diagram is essential to the investigation of the thermodynamic performances of [bmim]Zn2Cl5/NH3 absorption refrigeration system.

    Key words∶heat capacity; excess enthalpy; ionic liquid;[bmim]Zn2Cl5/NH3; NRTL model

    中圖分類號:TB64

    文獻(xiàn)標(biāo)識碼:A

    文章編號:1002-4026(2016)02-0049-08

    作者簡介:陳偉(1986-),男,博士,研究方向?yàn)槲帐街评洹?通訊作者,梁世強(qiáng)。Email:liangsq@iet.cn

    基金項(xiàng)目:國家自然科學(xué)基金(51276180);山東省優(yōu)秀中青年科學(xué)家科研獎(jiǎng)勵(lì)基金(BS2014NJ021)

    收稿日期:2016-01-19

    DOI:10.3976/j.issn.1002-4026.2016.02.010

    猜你喜歡
    離子液體比熱容
    比熱容知識知多少
    話說物質(zhì)的比熱容
    細(xì)說比熱容
    多視角解讀比熱容
    離子液體—基質(zhì)固相分散—超聲霧化—固相萃取結(jié)合高效液相色譜法檢測人參中三嗪類除草劑
    解讀比熱容
    基于離子液體——綠色溶劑的制備和應(yīng)用分析
    科技傳播(2016年11期)2016-07-20 23:32:52
    分離分析中離子液體的應(yīng)用進(jìn)展
    科技傳播(2016年10期)2016-07-15 23:13:53
    離子液體在潤滑行業(yè)中的應(yīng)用及研究進(jìn)展
    潤滑油(2016年1期)2016-07-13 04:01:58
    離子色譜法測污染物中的有機(jī)陽離子
    商(2016年4期)2016-03-24 00:28:51
    久久精品久久久久久久性| 街头女战士在线观看网站| 一本色道久久久久久精品综合| 日韩成人av中文字幕在线观看| 少妇 在线观看| 国产精品久久久av美女十八| 全区人妻精品视频| 免费不卡的大黄色大毛片视频在线观看| 午夜福利视频精品| 黄色 视频免费看| 国产精品国产av在线观看| 熟女电影av网| 各种免费的搞黄视频| 在线观看免费高清a一片| 超碰97精品在线观看| 另类亚洲欧美激情| 只有这里有精品99| 两个人看的免费小视频| 超碰97精品在线观看| 熟女av电影| 中文天堂在线官网| av一本久久久久| 男人爽女人下面视频在线观看| 久久人人97超碰香蕉20202| 国产探花极品一区二区| 男人操女人黄网站| 国产综合精华液| 人妻系列 视频| 国产午夜精品一二区理论片| 一级,二级,三级黄色视频| 国产男女内射视频| 欧美国产精品va在线观看不卡| 精品人妻一区二区三区麻豆| 肉色欧美久久久久久久蜜桃| 国产成人精品婷婷| 99热6这里只有精品| 午夜视频国产福利| 只有这里有精品99| 久久 成人 亚洲| 欧美日韩亚洲高清精品| 久久国产精品大桥未久av| 亚洲国产av影院在线观看| 一二三四中文在线观看免费高清| 国产在线一区二区三区精| 夫妻性生交免费视频一级片| 男女国产视频网站| 麻豆精品久久久久久蜜桃| 考比视频在线观看| 一区在线观看完整版| 视频中文字幕在线观看| 久热久热在线精品观看| 日韩熟女老妇一区二区性免费视频| av免费观看日本| 色视频在线一区二区三区| 亚洲精品第二区| 国产精品国产av在线观看| 亚洲精品自拍成人| 日韩成人av中文字幕在线观看| 亚洲精品中文字幕在线视频| 久久精品人人爽人人爽视色| 国产精品秋霞免费鲁丝片| av卡一久久| 九九爱精品视频在线观看| 少妇高潮的动态图| 中国三级夫妇交换| a 毛片基地| 亚洲人成网站在线观看播放| 2022亚洲国产成人精品| 久久久久人妻精品一区果冻| 国产成人欧美| 熟女电影av网| 自拍欧美九色日韩亚洲蝌蚪91| 在线 av 中文字幕| 男人添女人高潮全过程视频| 国产一区二区激情短视频 | 国产伦理片在线播放av一区| 久久久久精品久久久久真实原创| 久久国内精品自在自线图片| 美女中出高潮动态图| 高清黄色对白视频在线免费看| 国产国拍精品亚洲av在线观看| 中文乱码字字幕精品一区二区三区| 精品国产一区二区三区四区第35| 2022亚洲国产成人精品| 久久国产亚洲av麻豆专区| 高清毛片免费看| 波野结衣二区三区在线| 中文精品一卡2卡3卡4更新| 9热在线视频观看99| 久久久精品免费免费高清| 午夜日本视频在线| 久久97久久精品| 中文字幕免费在线视频6| 看十八女毛片水多多多| 国产欧美日韩综合在线一区二区| 一本大道久久a久久精品| 国产精品久久久av美女十八| 欧美日韩成人在线一区二区| 国产成人免费无遮挡视频| 国产在线一区二区三区精| 中文字幕人妻丝袜制服| 男女边吃奶边做爰视频| 欧美3d第一页| 一级毛片黄色毛片免费观看视频| 欧美少妇被猛烈插入视频| 少妇熟女欧美另类| 欧美激情极品国产一区二区三区 | 天天影视国产精品| 男女啪啪激烈高潮av片| 久久97久久精品| 亚洲国产精品国产精品| 满18在线观看网站| 久久人人爽av亚洲精品天堂| 亚洲国产成人一精品久久久| 永久网站在线| 91久久精品国产一区二区三区| 天堂俺去俺来也www色官网| 精品人妻熟女毛片av久久网站| 免费播放大片免费观看视频在线观看| 99久国产av精品国产电影| 免费人成在线观看视频色| 亚洲精品自拍成人| 美女内射精品一级片tv| 在线观看免费视频网站a站| 亚洲精品日韩在线中文字幕| 欧美性感艳星| 狠狠精品人妻久久久久久综合| 色94色欧美一区二区| 成人国产av品久久久| www.熟女人妻精品国产 | 日韩制服骚丝袜av| 下体分泌物呈黄色| 亚洲国产精品999| 99香蕉大伊视频| 美女内射精品一级片tv| 日韩电影二区| 久久精品国产自在天天线| 精品国产国语对白av| 午夜影院在线不卡| 男人操女人黄网站| 成年人午夜在线观看视频| 赤兔流量卡办理| 香蕉国产在线看| 国产精品久久久久久av不卡| 亚洲国产毛片av蜜桃av| 日韩av不卡免费在线播放| 久久人人爽av亚洲精品天堂| 日日啪夜夜爽| 亚洲一区二区三区欧美精品| 美女内射精品一级片tv| 宅男免费午夜| 我的女老师完整版在线观看| 国产精品一二三区在线看| 大香蕉久久网| 国产激情久久老熟女| 在线观看www视频免费| 性高湖久久久久久久久免费观看| 高清av免费在线| 黄网站色视频无遮挡免费观看| 黄色配什么色好看| 亚洲成国产人片在线观看| 国产成人aa在线观看| 欧美激情国产日韩精品一区| 精品午夜福利在线看| 韩国高清视频一区二区三区| 国产成人精品久久久久久| www.熟女人妻精品国产 | 一级片'在线观看视频| 亚洲精品成人av观看孕妇| 九草在线视频观看| 日韩av免费高清视频| 国产69精品久久久久777片| 久久精品人人爽人人爽视色| 欧美另类一区| 国产欧美日韩一区二区三区在线| 精品福利永久在线观看| 美女内射精品一级片tv| 高清视频免费观看一区二区| 国产一区亚洲一区在线观看| 久久久久久人人人人人| 亚洲精品乱码久久久久久按摩| 久久青草综合色| 国产精品麻豆人妻色哟哟久久| 亚洲精品自拍成人| 亚洲,一卡二卡三卡| 夜夜爽夜夜爽视频| 中文字幕人妻熟女乱码| 亚洲欧美一区二区三区国产| 日韩熟女老妇一区二区性免费视频| 高清av免费在线| 欧美 日韩 精品 国产| 日本-黄色视频高清免费观看| 欧美另类一区| av在线老鸭窝| av播播在线观看一区| 国语对白做爰xxxⅹ性视频网站| 80岁老熟妇乱子伦牲交| 亚洲av.av天堂| 最近的中文字幕免费完整| 国产白丝娇喘喷水9色精品| 国产男人的电影天堂91| 久久久久久人妻| 精品人妻偷拍中文字幕| 亚洲色图综合在线观看| 赤兔流量卡办理| 精品国产一区二区久久| 国产 精品1| 黄色配什么色好看| 亚洲国产精品一区三区| 新久久久久国产一级毛片| 汤姆久久久久久久影院中文字幕| 国内精品宾馆在线| 国产精品.久久久| 日本猛色少妇xxxxx猛交久久| 中文字幕免费在线视频6| 母亲3免费完整高清在线观看 | 夜夜骑夜夜射夜夜干| 亚洲精品一二三| 少妇人妻久久综合中文| 亚洲av综合色区一区| 免费观看av网站的网址| 高清不卡的av网站| 成人午夜精彩视频在线观看| 爱豆传媒免费全集在线观看| 国产精品久久久久久久久免| 亚洲精品日韩在线中文字幕| 婷婷色综合www| 欧美人与善性xxx| 亚洲综合色网址| 国产爽快片一区二区三区| av网站免费在线观看视频| 黄色毛片三级朝国网站| 亚洲性久久影院| 一本色道久久久久久精品综合| 国产精品久久久久久精品古装| 午夜福利在线观看免费完整高清在| 国产成人av激情在线播放| tube8黄色片| 日本午夜av视频| 熟女人妻精品中文字幕| 亚洲欧美一区二区三区黑人 | 久久久久久久大尺度免费视频| av网站免费在线观看视频| 亚洲婷婷狠狠爱综合网| 欧美激情极品国产一区二区三区 | 少妇高潮的动态图| 搡老乐熟女国产| 精品亚洲成国产av| 亚洲美女视频黄频| 国产精品久久久久久久久免| 国产黄频视频在线观看| 日韩 亚洲 欧美在线| 欧美最新免费一区二区三区| av在线播放精品| 在线观看免费视频网站a站| 亚洲伊人久久精品综合| 美女视频免费永久观看网站| 国产精品偷伦视频观看了| 免费日韩欧美在线观看| 免费高清在线观看日韩| 午夜福利视频精品| 久久久久精品人妻al黑| 人人澡人人妻人| 秋霞伦理黄片| 国产一区亚洲一区在线观看| 在线观看人妻少妇| 日本vs欧美在线观看视频| 水蜜桃什么品种好| 国产精品不卡视频一区二区| 1024视频免费在线观看| 日韩大片免费观看网站| 亚洲天堂av无毛| 91精品国产国语对白视频| 夜夜爽夜夜爽视频| 香蕉国产在线看| 美女福利国产在线| 夜夜爽夜夜爽视频| 国产极品天堂在线| 激情五月婷婷亚洲| 国产乱来视频区| 国产色婷婷99| 亚洲av中文av极速乱| 午夜福利视频精品| 巨乳人妻的诱惑在线观看| 久久久亚洲精品成人影院| 国产一区亚洲一区在线观看| 一级黄片播放器| 自线自在国产av| 国产一区二区三区综合在线观看 | 国产成人精品婷婷| av国产久精品久网站免费入址| 婷婷色av中文字幕| kizo精华| 伊人亚洲综合成人网| 十分钟在线观看高清视频www| 精品卡一卡二卡四卡免费| 另类精品久久| 18禁国产床啪视频网站| 超色免费av| 在线观看国产h片| 国产乱来视频区| 考比视频在线观看| 国产精品成人在线| 欧美精品国产亚洲| 大香蕉97超碰在线| 国产乱人偷精品视频| 久久精品国产亚洲av涩爱| 亚洲欧美一区二区三区黑人 | 五月伊人婷婷丁香| 精品一区在线观看国产| 亚洲四区av| 亚洲欧美色中文字幕在线| 黄色一级大片看看| 国产极品天堂在线| 免费av中文字幕在线| 91精品国产国语对白视频| 热re99久久精品国产66热6| 久久av网站| 精品国产乱码久久久久久小说| a 毛片基地| 成年人免费黄色播放视频| 国产麻豆69| 亚洲欧美中文字幕日韩二区| 乱码一卡2卡4卡精品| 中文字幕人妻丝袜制服| av一本久久久久| 韩国av在线不卡| 97超碰精品成人国产| 黄片无遮挡物在线观看| 久久97久久精品| 免费观看av网站的网址| 午夜激情久久久久久久| 激情视频va一区二区三区| 成人毛片60女人毛片免费| 亚洲精品国产av成人精品| 国产男女超爽视频在线观看| av国产久精品久网站免费入址| 女性生殖器流出的白浆| 99久久综合免费| 看免费成人av毛片| 高清不卡的av网站| 成人亚洲精品一区在线观看| 精品国产国语对白av| videosex国产| 黄色毛片三级朝国网站| 五月开心婷婷网| 日韩成人av中文字幕在线观看| 久久av网站| 制服人妻中文乱码| 国产伦理片在线播放av一区| 色视频在线一区二区三区| 欧美精品高潮呻吟av久久| 亚洲一区二区三区欧美精品| 日韩一区二区三区影片| 久久久亚洲精品成人影院| 国产免费视频播放在线视频| 亚洲精品美女久久久久99蜜臀 | 久久久久国产网址| 少妇猛男粗大的猛烈进出视频| 欧美激情极品国产一区二区三区 | 制服丝袜香蕉在线| 老司机影院成人| 国产免费福利视频在线观看| 曰老女人黄片| 天天躁夜夜躁狠狠久久av| 精品一区二区三区视频在线| 天天操日日干夜夜撸| 久久久久网色| 纵有疾风起免费观看全集完整版| 国产黄色视频一区二区在线观看| 日韩av在线免费看完整版不卡| 久久久久国产精品人妻一区二区| 久久免费观看电影| 秋霞在线观看毛片| 国产精品一国产av| 如何舔出高潮| 又大又黄又爽视频免费| 乱人伦中国视频| 亚洲av成人精品一二三区| 国产爽快片一区二区三区| 成人国产麻豆网| 日韩制服丝袜自拍偷拍| 国产一区二区在线观看av| 97在线人人人人妻| 777米奇影视久久| 亚洲欧美色中文字幕在线| 亚洲av福利一区| 精品国产一区二区三区四区第35| 黑人猛操日本美女一级片| 国产成人午夜福利电影在线观看| 国产成人av激情在线播放| 十八禁网站网址无遮挡| 亚洲成人手机| 亚洲第一区二区三区不卡| 亚洲精品日本国产第一区| 成人毛片a级毛片在线播放| 国产精品99久久99久久久不卡 | 精品熟女少妇av免费看| 中文乱码字字幕精品一区二区三区| 天天操日日干夜夜撸| 国产黄色视频一区二区在线观看| 日本与韩国留学比较| 啦啦啦啦在线视频资源| 最黄视频免费看| 韩国精品一区二区三区 | 一区二区日韩欧美中文字幕 | 精品国产一区二区三区四区第35| 男女高潮啪啪啪动态图| 久久精品国产亚洲av涩爱| 高清不卡的av网站| 插逼视频在线观看| 五月天丁香电影| 天天躁夜夜躁狠狠久久av| 麻豆乱淫一区二区| 国产 一区精品| 如何舔出高潮| 桃花免费在线播放| h视频一区二区三区| 搡老乐熟女国产| 91精品三级在线观看| 卡戴珊不雅视频在线播放| 久久久久久久久久人人人人人人| av电影中文网址| 欧美xxⅹ黑人| 一边摸一边做爽爽视频免费| 亚洲第一区二区三区不卡| 免费观看无遮挡的男女| 秋霞伦理黄片| 国产成人91sexporn| 一二三四在线观看免费中文在 | 看非洲黑人一级黄片| 日韩 亚洲 欧美在线| 日韩av在线免费看完整版不卡| 建设人人有责人人尽责人人享有的| 日本色播在线视频| 乱人伦中国视频| 亚洲精品色激情综合| 亚洲国产看品久久| av天堂久久9| 日韩,欧美,国产一区二区三区| 黄色毛片三级朝国网站| 亚洲第一区二区三区不卡| 欧美精品亚洲一区二区| 一区在线观看完整版| 老熟女久久久| 国产男人的电影天堂91| 亚洲,欧美精品.| 啦啦啦中文免费视频观看日本| 国产男女超爽视频在线观看| 亚洲国产精品成人久久小说| av在线观看视频网站免费| 女的被弄到高潮叫床怎么办| 国产精品国产三级国产av玫瑰| 精品一区二区三卡| 伦理电影免费视频| 亚洲精品美女久久av网站| 热re99久久国产66热| 寂寞人妻少妇视频99o| 亚洲精品第二区| 一区在线观看完整版| 97在线人人人人妻| 97在线视频观看| 不卡视频在线观看欧美| 夜夜骑夜夜射夜夜干| 1024视频免费在线观看| 纯流量卡能插随身wifi吗| 成人国语在线视频| 18禁国产床啪视频网站| 99久久综合免费| 成人毛片60女人毛片免费| 99热这里只有是精品在线观看| 美女主播在线视频| 免费黄频网站在线观看国产| 久久久精品94久久精品| 最近中文字幕高清免费大全6| 色94色欧美一区二区| www日本在线高清视频| 超色免费av| 高清毛片免费看| 午夜激情久久久久久久| 在线精品无人区一区二区三| 最近手机中文字幕大全| 欧美亚洲 丝袜 人妻 在线| 在线亚洲精品国产二区图片欧美| 免费不卡的大黄色大毛片视频在线观看| 欧美变态另类bdsm刘玥| 在线观看免费高清a一片| 大话2 男鬼变身卡| 成人二区视频| 爱豆传媒免费全集在线观看| 亚洲,欧美精品.| 王馨瑶露胸无遮挡在线观看| 一级,二级,三级黄色视频| 中文字幕亚洲精品专区| 久久久国产一区二区| 亚洲久久久国产精品| 午夜老司机福利剧场| 1024视频免费在线观看| 黑人巨大精品欧美一区二区蜜桃 | 国产毛片在线视频| 性色avwww在线观看| 国产黄色免费在线视频| 男人舔女人的私密视频| 一级a做视频免费观看| 亚洲久久久国产精品| 国产在视频线精品| 一区二区三区精品91| 97在线人人人人妻| 国产精品一区二区在线不卡| 久久国产精品大桥未久av| 亚洲欧美一区二区三区国产| 亚洲欧洲国产日韩| 午夜免费观看性视频| 亚洲精品久久午夜乱码| 亚洲欧美一区二区三区黑人 | 九草在线视频观看| 美女国产高潮福利片在线看| 最黄视频免费看| 亚洲国产看品久久| 日日摸夜夜添夜夜爱| 亚洲欧美一区二区三区黑人 | 看免费成人av毛片| 亚洲国产精品成人久久小说| 精品国产乱码久久久久久小说| 男人添女人高潮全过程视频| 国精品久久久久久国模美| 天堂中文最新版在线下载| xxxhd国产人妻xxx| 制服丝袜香蕉在线| 如何舔出高潮| 九色成人免费人妻av| 欧美激情极品国产一区二区三区 | 亚洲熟女精品中文字幕| 国产又爽黄色视频| 黄色怎么调成土黄色| 男女免费视频国产| 中文字幕制服av| 免费久久久久久久精品成人欧美视频 | 久久久久精品人妻al黑| 久久久久网色| 国产极品粉嫩免费观看在线| 精品视频人人做人人爽| 午夜福利视频精品| 亚洲精品日韩在线中文字幕| 涩涩av久久男人的天堂| 青春草亚洲视频在线观看| 国产有黄有色有爽视频| 亚洲av综合色区一区| 在线观看三级黄色| 国产一区二区激情短视频 | 久久国产精品男人的天堂亚洲 | 少妇被粗大的猛进出69影院 | 国产又色又爽无遮挡免| 午夜免费鲁丝| 免费av不卡在线播放| 97精品久久久久久久久久精品| 亚洲伊人久久精品综合| 国产毛片在线视频| 老司机影院成人| 久久精品久久久久久噜噜老黄| 80岁老熟妇乱子伦牲交| 日韩欧美精品免费久久| 一区二区三区精品91| 国产成人精品一,二区| 国产av国产精品国产| 天天躁夜夜躁狠狠久久av| 国产乱人偷精品视频| 欧美性感艳星| 高清黄色对白视频在线免费看| 中文字幕制服av| www.熟女人妻精品国产 | 少妇熟女欧美另类| 亚洲成国产人片在线观看| 亚洲精品美女久久av网站| 欧美精品亚洲一区二区| 亚洲av男天堂| 国产探花极品一区二区| 欧美精品亚洲一区二区| 亚洲av男天堂| 久久精品国产亚洲av涩爱| 欧美日韩av久久| 亚洲国产av新网站| 深夜精品福利| www日本在线高清视频| 亚洲欧美成人综合另类久久久| 97超碰精品成人国产| 欧美 亚洲 国产 日韩一| 精品国产乱码久久久久久小说| 国产午夜精品一二区理论片| 婷婷色综合www| 少妇被粗大猛烈的视频| 国产成人午夜福利电影在线观看| 我要看黄色一级片免费的| 国产一区二区在线观看日韩| 国产精品国产av在线观看| 欧美人与性动交α欧美软件 | 亚洲第一av免费看| 毛片一级片免费看久久久久| 国产精品偷伦视频观看了| 精品视频人人做人人爽| 成人无遮挡网站| 国产精品无大码| 日韩欧美一区视频在线观看| 午夜福利影视在线免费观看| av福利片在线| kizo精华| 久久久久网色| 美女国产视频在线观看| 欧美日本中文国产一区发布| 在线免费观看不下载黄p国产| 亚洲精品乱久久久久久| 五月开心婷婷网| 国产69精品久久久久777片| 精品久久国产蜜桃| av.在线天堂|