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

    正丙醇黏度方程研究

    2017-11-22 10:19:59孟現(xiàn)陽(yáng)于野吳江濤
    化工學(xué)報(bào) 2017年11期
    關(guān)鍵詞:正丙醇偏差黏度

    孟現(xiàn)陽(yáng),于野,吳江濤

    (西安交通大學(xué)能源與動(dòng)力工程學(xué)院,熱流科學(xué)與工程教育部重點(diǎn)實(shí)驗(yàn)室,陜西 西安 710049)

    正丙醇黏度方程研究

    孟現(xiàn)陽(yáng),于野,吳江濤

    (西安交通大學(xué)能源與動(dòng)力工程學(xué)院,熱流科學(xué)與工程教育部重點(diǎn)實(shí)驗(yàn)室,陜西 西安 710049)

    正丙醇被廣泛應(yīng)用于化工、醫(yī)藥和農(nóng)業(yè)等領(lǐng)域,可靠的熱物性數(shù)據(jù)是理論研究和工程應(yīng)用的基礎(chǔ)。黏度作為流體重要的熱物理性質(zhì)之一,在生產(chǎn)流程的設(shè)計(jì)與優(yōu)化和產(chǎn)品質(zhì)量監(jiān)測(cè)等方面具有重要意義。通過(guò)實(shí)驗(yàn)測(cè)量得到的黏度數(shù)據(jù)一般離散分布于不同的熱力學(xué)狀態(tài),在工程中難以直接應(yīng)用,因此,建立可靠的黏度方程以實(shí)現(xiàn)任意狀態(tài)黏度的計(jì)算非常必要。經(jīng)文獻(xiàn)調(diào)研發(fā)現(xiàn),目前缺少寬范圍、高精度的正丙醇黏度方程,基于國(guó)內(nèi)外已發(fā)表的實(shí)驗(yàn)數(shù)據(jù)開(kāi)發(fā)了正丙醇高精度黏度方程。該方程不僅能夠很好地復(fù)現(xiàn)實(shí)驗(yàn)數(shù)據(jù),而且具有合理的外推性,其溫度適用范圍為150~400 K,壓力適用范圍為0.1~120 MPa,方程在計(jì)算范圍內(nèi)估計(jì)不確定度為1%~5%。

    正丙醇;黏度;模型;熱力學(xué)

    引 言

    正丙醇被廣泛應(yīng)用于化工、醫(yī)藥和農(nóng)業(yè)等領(lǐng)域。在化工生產(chǎn)中,可以作為原料參與乙酸丙酯的生產(chǎn);在醫(yī)藥生產(chǎn)中,可以用其制取黏合止血?jiǎng)?、紅霉素等;在農(nóng)業(yè)中,其衍生物可以用于制取菌達(dá)滅和滅草蜢等,此外其衍生出的醇醚類(lèi)物質(zhì)可以作為添加劑加入飼料中。正丙醇在不同領(lǐng)域的應(yīng)用離不開(kāi)準(zhǔn)確的熱物性數(shù)據(jù),黏度作為其中一種重要的熱物性參數(shù),對(duì)于生產(chǎn)流程的設(shè)計(jì)與優(yōu)化和產(chǎn)品質(zhì)量監(jiān)測(cè)等具有重要意義,因此對(duì)黏度進(jìn)行研究有利于正丙醇在工程中的實(shí)際應(yīng)用。雖然每年都有大量的黏度實(shí)驗(yàn)數(shù)據(jù)發(fā)表,但是這些數(shù)據(jù)一般離散分布于不同的溫度和壓力下,在工業(yè)應(yīng)用中難以直接應(yīng)用,因此有必要建立可靠的黏度方程以實(shí)現(xiàn)數(shù)據(jù)的連續(xù)計(jì)算。本文基于文獻(xiàn)中公開(kāi)發(fā)表的黏度實(shí)驗(yàn)數(shù)據(jù),采用多參數(shù)黏度模型,開(kāi)發(fā)了高精度、寬范圍的正丙醇黏度方程。

    1 黏度實(shí)驗(yàn)數(shù)據(jù)

    通過(guò)系統(tǒng)全面的文獻(xiàn)調(diào)研,對(duì)正丙醇黏度實(shí)驗(yàn)數(shù)據(jù)進(jìn)行了廣泛搜集,共得到132篇正丙醇黏度實(shí)驗(yàn)數(shù)據(jù)文獻(xiàn),共包含538個(gè)實(shí)驗(yàn)數(shù)據(jù),其中7篇文獻(xiàn)包含高壓狀態(tài)下的實(shí)驗(yàn)數(shù)據(jù),3篇文獻(xiàn)包含 0.1 MPa下氣相實(shí)驗(yàn)數(shù)據(jù),其余數(shù)據(jù)均為常壓液相實(shí)驗(yàn)數(shù)據(jù)。然而并非所有數(shù)據(jù)都可以用于方程開(kāi)發(fā),不同作者、不同年代、不同的實(shí)驗(yàn)方法得到的實(shí)驗(yàn)數(shù)據(jù)相互之間存在較大差異,因此需要對(duì)實(shí)驗(yàn)數(shù)據(jù)進(jìn)行有效評(píng)價(jià),篩選出可靠的黏度實(shí)驗(yàn)數(shù)據(jù)用于方程的開(kāi)發(fā)。本文將搜集到的黏度實(shí)驗(yàn)數(shù)據(jù)分為了兩類(lèi):第一類(lèi)是基本數(shù)據(jù),這類(lèi)數(shù)據(jù)用于方程的開(kāi)發(fā);第二類(lèi)數(shù)據(jù)主要用于方程的檢驗(yàn)。數(shù)據(jù)的篩選遵循以下原則[1]:(1)文獻(xiàn)中黏度實(shí)驗(yàn)測(cè)量方法具備完善的工作方程,并且其中的相關(guān)參數(shù)能夠通過(guò)標(biāo)定或者實(shí)驗(yàn)測(cè)量得到;(2)文獻(xiàn)中應(yīng)包含實(shí)驗(yàn)所用樣品的純度信息;(3)文獻(xiàn)中應(yīng)給出實(shí)驗(yàn)數(shù)據(jù)的不確定度分析。此外,還需考慮實(shí)驗(yàn)數(shù)據(jù)的溫度、壓力范圍,應(yīng)盡可能拓展方程的適用范圍。

    按照以上數(shù)據(jù)篩選原則,對(duì)132篇文獻(xiàn)中正丙醇的黏度實(shí)驗(yàn)數(shù)據(jù)進(jìn)行了詳細(xì)分析和評(píng)價(jià)。在常壓實(shí)驗(yàn)數(shù)據(jù)篩選中,選取了溫度覆蓋范圍較廣、精度較高的6 篇文獻(xiàn)實(shí)驗(yàn)數(shù)據(jù)作為基本數(shù)據(jù)[2-7]。其中,Komarenko等[2]的實(shí)驗(yàn)溫度范圍為153~298 K,其最低溫度接近三相點(diǎn)溫度。Kumagai等[3]、Tu等[4]、Saleh等[5]、Yang等[6]和 Pang等[7]采用了毛細(xì)管黏度計(jì),測(cè)量溫度范圍為273~363 K,測(cè)量不確定度大部分小于 1%。高壓黏度實(shí)驗(yàn)數(shù)據(jù)[8-13],Weber[8]的實(shí)驗(yàn)方法為滾球黏度計(jì),實(shí)驗(yàn)壓力范圍為 0.1~49.1 MPa,溫度范圍為 273~373 K,不確定度為1.5%。Tanaka等[9]的實(shí)驗(yàn)方法為落體黏度計(jì),實(shí)驗(yàn)壓力范圍為0.1~117.8 MPa,溫度范圍為283~323 K,不確定度為2%。Papaioannou等[10]、Assael等[11]的實(shí)驗(yàn)方法為落體黏度計(jì),實(shí)驗(yàn)壓力范圍分別為0.1~71.75 MPa和0.1~51.8 MPa,溫度為298 K,不確定度小于2.5%。Papaioannou等[12]的實(shí)驗(yàn)方法為振動(dòng)弦黏度計(jì),實(shí)驗(yàn)壓力范圍為0.1~27.86 MPa,溫度范圍為 295~328 K,不確定度為 0.5%。Baylaucq等[13]的實(shí)驗(yàn)方法為落體黏度計(jì),實(shí)驗(yàn)壓力范圍為0.1~100 MPa,溫度范圍為293~353 K,不確定度為 2%。通過(guò)對(duì)高壓黏度數(shù)據(jù)的對(duì)比,并結(jié)合方程的適用范圍,本文選用以上6篇文獻(xiàn)高壓實(shí)驗(yàn)數(shù)據(jù)作為基本數(shù)據(jù)。最終用于正丙醇黏度方程開(kāi)發(fā)的實(shí)驗(yàn)數(shù)據(jù)包括12篇文獻(xiàn),共計(jì)178個(gè)數(shù)據(jù)點(diǎn),表1給出了基本實(shí)驗(yàn)數(shù)據(jù)相關(guān)信息。

    表1 正丙醇黏度基本實(shí)驗(yàn)數(shù)據(jù)Table 1 Primary data used in developing viscosity correlation of n-propanol

    表2 密度方程式(1)~式(4)的系數(shù)Table 2 Coefficients for Eq. (1) to Eq.(4)

    2 黏度方程形式

    根據(jù)現(xiàn)有黏度理論研究和已有黏度方程開(kāi)發(fā)經(jīng)驗(yàn)[19-24],純質(zhì)黏度的計(jì)算以溫度和密度為自變量,并由四項(xiàng)互相獨(dú)立的項(xiàng)組成。

    式中,η0(T)為零密度項(xiàng);η1(T)ρ為初始密度依賴(lài)項(xiàng),通常認(rèn)為在0.2 MPa內(nèi)黏度和密度之間是線(xiàn)性關(guān)系[23];Δηc(ρ,T)為臨界增強(qiáng)黏度項(xiàng),作用是修正黏度在臨界點(diǎn)附近區(qū)域的奇異現(xiàn)象,已有文獻(xiàn)研究表明[19-24],對(duì)于黏度僅在臨界點(diǎn)溫度1%~2%范圍內(nèi),Δηc/η可能超過(guò)0.01。此外,鑒于所搜集的黏度實(shí)驗(yàn)數(shù)據(jù)中沒(méi)有近臨界區(qū)的實(shí)驗(yàn)數(shù)據(jù),無(wú)法開(kāi)發(fā)準(zhǔn)確的臨界增強(qiáng)型表達(dá)式,參考已有文獻(xiàn)研究,本文采用了同樣的辦法,將該項(xiàng)設(shè)為零。Δη(ρ,T)為剩余黏度項(xiàng),該項(xiàng)幾乎沒(méi)有理論指導(dǎo),其形式主要是參考已有黏度模型并結(jié)合實(shí)驗(yàn)數(shù)據(jù)得到。

    2.1 零密度項(xiàng)和初始密度依賴(lài)項(xiàng)

    該項(xiàng)的主要理論基礎(chǔ)為分子運(yùn)動(dòng)理論,為了更準(zhǔn)確地復(fù)現(xiàn)實(shí)驗(yàn)數(shù)據(jù),本文采用了Chung等[25]提出的修正模型,該模型引入了參數(shù)Fc用于修正分子結(jié)構(gòu)、分子極性和締合性質(zhì)對(duì)正丙醇黏度的影響。零密度下的黏度計(jì)算方程如式(6)~式(10)所示。

    式(6)~式(10)中需要確定的關(guān)鍵參數(shù)為勢(shì)能尺度特征參數(shù)σ和勢(shì)能能量特征參數(shù)ε/kB,確定方法分為兩種,一種為關(guān)聯(lián)臨界參數(shù)[25],另一種為關(guān)聯(lián)實(shí)驗(yàn)數(shù)據(jù)進(jìn)行求解,本文采用了后者。因沒(méi)有搜集到正丙醇在極低密度下的黏度實(shí)驗(yàn)數(shù)據(jù),故無(wú)法通過(guò)推導(dǎo)得到零密度下氣體的黏度值,本文將η0(T)+η1(T)ρ作為低密度黏度項(xiàng)與Titani[26]的實(shí)驗(yàn)數(shù)據(jù)相關(guān)聯(lián),并結(jié)合Golubev等[27]的密度實(shí)驗(yàn)數(shù)據(jù),回歸得到σ和ε/kB。式(6)~式(10)中的相關(guān)參數(shù)、系數(shù)和指數(shù)見(jiàn)表3。

    文獻(xiàn)研究表明,無(wú)論是極性還是非極性物質(zhì),在一定溫度范圍內(nèi),低壓下氣體的黏度會(huì)沿等溫線(xiàn)隨密度的增加而降低,這種現(xiàn)象可以通過(guò)第二黏度維里系數(shù)進(jìn)行量化說(shuō)明,其主要理論依據(jù)為Rainwater and Friend 理論[28-30],并可由式(1)關(guān)聯(lián)得到。

    為了方便與實(shí)驗(yàn)數(shù)據(jù)進(jìn)行對(duì)比,引入量綱1第二黏度virial系數(shù)B*η(T),該系數(shù)與溫度的關(guān)系式由文獻(xiàn)[31]給出,其溫度適用范圍為 0.3

    式中,NA為Avogadro常數(shù),式(13)中的相關(guān)系數(shù)見(jiàn)表3。

    2.2 剩余黏度項(xiàng)

    剩余黏度項(xiàng)形式的選擇,需要考慮黏度在高密度區(qū)域的數(shù)量級(jí)變化特性以及方程的外推性等因素,結(jié)合正丙醇的黏度特性并參考已有黏度方程研究[19-24],最終確定剩余黏度項(xiàng)的形式如下

    表3 式(6)~式(13)中的系數(shù)Table 3 Coefficients for Eq. (6) to Eq.(13)

    表4 剩余黏度項(xiàng)系數(shù)和指數(shù)Table 4 Coefficients and parameters for Eq. (14)

    其中,Δη(ρ,T)的單位為 μPa?s,ρr=ρ/ρc,Tr=T/Tc。由于基本實(shí)驗(yàn)數(shù)據(jù)的不確定度不盡相同,因此需要對(duì)不同的實(shí)驗(yàn)數(shù)據(jù)進(jìn)行權(quán)重設(shè)置,初始權(quán)重采用不確定度平方分之一進(jìn)行設(shè)置,即w=1/(ηu)2。對(duì)于不確定度不明確的黏度實(shí)驗(yàn)數(shù)據(jù),參考其實(shí)驗(yàn)方法、發(fā)表年代和樣品純度信息等進(jìn)行初始權(quán)重設(shè)置。通過(guò)合理的權(quán)重設(shè)置,最終得到的擬合結(jié)果如表4所示。

    表5給出了黏度文獻(xiàn)實(shí)驗(yàn)數(shù)據(jù)與開(kāi)發(fā)的方程計(jì)算值偏差情況,方程很好地復(fù)現(xiàn)了實(shí)驗(yàn)數(shù)據(jù),其基本實(shí)驗(yàn)數(shù)據(jù)與方程計(jì)算值的平均偏差(Bias)為0.14%,平均絕對(duì)偏差(AAD)為 0.74%(Bias和AAD計(jì)算方法與密度一樣),最大偏差為-5.02%。

    3 實(shí)驗(yàn)數(shù)據(jù)分析

    下述實(shí)驗(yàn)數(shù)據(jù)分析來(lái)自于文獻(xiàn)[2-13,33-50]。圖1給出了0.1 MPa下溫度范圍為153~363 K的液相黏度實(shí)驗(yàn)數(shù)據(jù)與方程的計(jì)算值之間的偏差,圖2給出在壓力范圍為 0.1~117.8 MPa,溫度范圍為273~373 K下的實(shí)驗(yàn)數(shù)據(jù)與方程計(jì)算值間的偏差隨密度變化的分布,從圖中可以看出偏差基本在±2.5%以?xún)?nèi),且均勻分布在正負(fù)兩側(cè),沒(méi)有明顯的趨勢(shì)走向。

    表5 正丙醇黏度方程相關(guān)統(tǒng)計(jì)指標(biāo)Table 5 Evaluation of n-propanol viscosity correlation for primary data

    在常壓狀態(tài)下,Komarenko等[2]的實(shí)驗(yàn)數(shù)據(jù)的溫度范圍為153~298 K,黏度在低溫區(qū)域隨著密度的增加會(huì)呈現(xiàn)出數(shù)量級(jí)的變化,方程很好地復(fù)現(xiàn)了這種特性,實(shí)驗(yàn)數(shù)據(jù)與計(jì)算值的Bias為0.17%,AAD為0.31%;Yang等[6]的實(shí)驗(yàn)數(shù)據(jù)與方程的計(jì)算值的AAD為1.18%,與其不確定度相當(dāng);Pang等[7]的實(shí)驗(yàn)數(shù)據(jù)與方程的計(jì)算值的 AAD為 0.20%,在其不確定度0.4%以?xún)?nèi)。在高壓黏度實(shí)驗(yàn)數(shù)據(jù)中,Tanaka等[9]的實(shí)驗(yàn)數(shù)據(jù)與方程的計(jì)算值的AAD為0.89%,在其不確定度 2%內(nèi);Baylaucq等[13]的實(shí)驗(yàn)數(shù)據(jù)與方程的計(jì)算值的AAD為1.69%,在其不確定度2%內(nèi)。根據(jù)實(shí)驗(yàn)數(shù)據(jù)范圍,最終方程可適用于溫度范圍為150~400 K,壓力范圍為0.1~120 MPa的狀態(tài)區(qū)域內(nèi)。150~373 K、0.1 MPa范圍內(nèi)方程的不確定度為1%,273~373 K、壓力120 MPa以?xún)?nèi)方程的不確定度估計(jì)為 2%,其他區(qū)域根據(jù)范圍不同估計(jì)不確定度為2%~5%。

    圖1 正丙醇常壓下液相黏度實(shí)驗(yàn)數(shù)據(jù)與計(jì)算值偏差Fig.1 Percentage deviations of primary data at 0.1 MPa for n-propanol

    圖2 正丙醇高壓下液相黏度實(shí)驗(yàn)數(shù)據(jù)與計(jì)算值偏差Fig.2 Percentage deviations of primary data at high pressure for n-propanol

    圖3為方程在不同等溫線(xiàn)下黏度與密度的曲線(xiàn),由圖可以發(fā)現(xiàn)黏度的變化曲線(xiàn)平滑,為混合物黏度的預(yù)測(cè)奠定了基礎(chǔ)。表6給出了不同溫度和密度下方程的計(jì)算值,以便實(shí)際應(yīng)用中方程的檢驗(yàn)。

    在第二類(lèi)數(shù)據(jù)中,張建侯等[33]測(cè)量了298.15 K下壓力范圍為0.1~29.5 MPa的正丙醇黏度,其實(shí)驗(yàn)數(shù)據(jù)與計(jì)算值的Bias為1.50%,AAD為1.54%,雖然偏差不算太大,但該溫度下偏差隨壓力變化存在系統(tǒng)偏差。在剩余的第二類(lèi)數(shù)據(jù)中,部分?jǐn)?shù)據(jù)與其他數(shù)據(jù)有著明顯的偏差,例如 Mokhtarani等[34]采用滾球式黏度計(jì)測(cè)量了常壓下溫度范圍為 283~333 K的正丙醇黏度,其實(shí)驗(yàn)數(shù)據(jù)與計(jì)算值的偏差從 7.6%增大到 17.6%;Kermanpour等[35]采用了毛細(xì)管黏度計(jì)測(cè)量了溫度范圍為293~333 K的正丙醇黏度,其實(shí)驗(yàn)數(shù)據(jù)與計(jì)算值的最大偏差為-52%,結(jié)果明顯存在錯(cuò)誤。為了驗(yàn)證方程的合理性,對(duì)于第二類(lèi)數(shù)據(jù)首先排除了與其他實(shí)驗(yàn)數(shù)據(jù)有明顯偏差的文獻(xiàn)實(shí)驗(yàn)數(shù)據(jù),并選擇包含5個(gè)以上實(shí)驗(yàn)數(shù)據(jù)的文獻(xiàn)[36-53]與方程計(jì)算值進(jìn)行了對(duì)比,結(jié)果如圖4所示,第二類(lèi)數(shù)據(jù)偏差基本在±5%以?xún)?nèi),整體AAD為1.99%,Bias為-0.28%,證明所擬合的正丙醇黏度方程能夠很好地復(fù)現(xiàn)實(shí)驗(yàn)數(shù)據(jù)。

    圖3 正丙醇黏度沿不同等溫線(xiàn)隨密度的變化曲線(xiàn)Fig.3 Viscosity of n-propanol as a function of density along isotherms

    表6 正丙醇黏度方程計(jì)算驗(yàn)證值Table 6 Sample points for computer verification of correlating equations

    4 結(jié) 論

    圖4 第二類(lèi)數(shù)據(jù)0.1 MPa下液相黏度實(shí)驗(yàn)值與計(jì)算值偏差Fig.4 Percentage deviations of the secondary data

    本文對(duì)正丙醇公開(kāi)發(fā)表的黏度文獻(xiàn)實(shí)驗(yàn)數(shù)據(jù)進(jìn)行全面搜集和整理,對(duì)文獻(xiàn)中的黏度實(shí)驗(yàn)數(shù)據(jù)進(jìn)行了評(píng)價(jià)和篩選,利用篩選出的基本數(shù)據(jù)開(kāi)發(fā)了寬范圍高精度正丙醇黏度方程,基本實(shí)驗(yàn)數(shù)據(jù)與方程計(jì)算值的平均絕對(duì)偏差為 0.74%。方程的溫度適用范圍為 150~400 K,壓力適用范圍為 0.1~120 MPa,在適用范圍內(nèi)估計(jì)不確定度為 1%~5%。本文所開(kāi)發(fā)的正丙醇黏度方程不僅很好地復(fù)現(xiàn)了實(shí)驗(yàn)數(shù)據(jù),而且具有合理的外推特性。

    符 號(hào) 說(shuō) 明

    Bη——第二黏度 virial系數(shù),m3?kg-1

    M——物質(zhì)的摩爾質(zhì)量,g?mol-1

    NA——Avogadro常數(shù)

    p——壓力,MPa

    T——熱力學(xué)溫度,K

    T*——對(duì)比溫度

    Tc——臨界溫度,K

    Tr——臨界對(duì)比溫度

    vc——臨界比體積,cm3?mol-1

    η——?jiǎng)恿︷ざ龋蘌a·s

    ηc——臨界黏度增強(qiáng)項(xiàng),μPa·s

    η0——稀薄氣體黏度,μPa·s

    η1——密度線(xiàn)性相關(guān)項(xiàng),μPa·s·m3·kg-1

    Δη——剩余黏度,μPa·s

    κ——經(jīng)驗(yàn)性締合參數(shù)

    μ——偶極矩,D

    μr——對(duì)比偶極矩

    ρ——密度,kg?m-3

    ρc——臨界密度,kg?m-3

    ρr——臨界對(duì)比密度

    Ω*——分子碰撞積分

    ω——偏心因子

    下角標(biāo)

    exp ——實(shí)驗(yàn)值

    cal ——計(jì)算值

    [1] ASSAEL M J, RAMIRES M, NIETO DE CASTRO C A,et al.Benzene: a further liquid thermal conductivity standard[J]. Journal of Physical and Chemical Reference Data, 1990, 19(1): 113-117.

    [2] KOMARENKO V G, MANZHELII V G, RADTSIG A V. Density and viscosity of normal alcohols at low temperatures[J]. Ukrainian Physical Journal, 1967, 12: 676.

    [3] KUMAGAI A, YOKOYAMA C. Liquid viscosity of binary mixtures of methanol with ethanol and 1-propanol from 273.15 to 333.15 K[J].International Journal of Thermophysics, 1998, 19(1): 3-13.

    [4] TU C H, LIU C Y, WANG W F,et al. Volumetric and viscometric properties of binary mixtures of aliphatic alcohols (C1-C4) with nitroethane from 293.15 K to 313.15 K[J]. Journal of Chemical and Engineering Data, 2000, 45(3): 450-456.

    [5] SALEH M A, AKHTAR S, BEGUM S,et al. Density and viscosity of 1-alkanols[J]. Physics and Chemistry of Liquids, 2004, 42(6):615-623.

    [6] YANG C S, LAI H X, LIU Z G,et al. Density and viscosity of binary mixtures of diethyl carbonate with alcohols at (293.15 to 363.15) K and predictive results by UNIFAC-VISCO group contribution method[J]. Journal of Chemical and Engineering Data, 2006, 51(4):1345-1351.

    [7] PANG F M, SENG C E, TENG T T,et al. Densities and viscosities of aqueous solutions of 1-propanol and 2-propanol at temperatures from 293.15 K to 333.15 K[J]. Journal of Molecular Liquids, 2007, 136(1):71-78.

    [8] WEBER W. über die druckabh?ngigkeit der viskosit?t von alkohol-wasser-gemischen[J]. Rheologica Acta, 1975, 14(11):1012-1025.

    [9] TANAKA Y, MATSUDA Y, FUJIWARA H,et al. Viscosity of(water+alcohol) mixtures under high pressure[J]. International Journal of Thermophysics, 1987, 8(2): 147-163.

    [10] PAPAIOANNOU D, BRIDAKIS M, PANAYIOTOU C G. Excess dynamic viscosity and excess volume ofn-butylamine + 1-alkanol mixtures at moderately high pressures[J]. Journal of Chemical and Engineering Data, 1993, 38(3): 370-378.

    [11] ASSAEL M J, POLIMATIDOU S K. Measurements of the viscosity of alcohols in the temperature range 290—340 K at pressures up to 30 MPa[J]. International Journal of Thermophysics, 1994, 15(1): 95-107.

    [12] PAPAIOANNOU D, PANAYIOTOU C. Viscosity of alkanol + alkane mixtures at moderately high pressures[J]. Journal of Chemical and Engineering Data, 1994, 39(3): 463-466.

    [13] BAYLAUCQ A, WATSON G, ZéBERG-MIKKELSEN C K,et al.Dynamic viscosity of the binary system 1-propanol+ toluene as a function of temperature and pressure[J]. Journal of Chemical and Engineering Data, 2009, 54(9): 2715-2721.

    [14] ZAKHARYAEV Z R. Experimental observation ofpVTproperties of 1-propanol[J]. Inzh.-Fiz. Zh, 1982, 43(3): 796-798.

    [15] Zú?IGA-MORENO A, GALICIA-LUNA L A. Densities of 1-propanol and 2-propanolviaa vibrating tube densimeter from 313 to 363 K and up to 25 MPa[J]. Journal of Chemical and Engineering Data, 2002, 47(2): 155-160.

    [16] KITAJIMA H, KAGAWA N, ENDO H,et al. Isochoric heat capacities of alkanols and their aqueous mixtures[J]. Journal of Chemical and Engineering Data, 2003, 48(6): 1583-1586.

    [17] WATSON G, LAFITTE T, ZéBERG-MIKKELSEN C K,et al.Volumetric and derivative properties under pressure for the system 1-propanol + toluene: a discussion of PC-SAFT and SAFT-VR[J].Fluid Phase Equilibria, 2006, 247(1): 121-134.

    [18] DáVILA M J, ALCALDE R, ATILHAN M,et al.pρTmeasurements and derived properties of liquid 1-alkanols[J]. Journal of Chemical Thermodynamics, 2012, 47(3): 241–259.

    [19] VOGEL E, HERRMANN S. New Formulation for the viscosity of propane[J]. Journal of Physical and Chemical Reference Data, 2016,45(4): 043103.

    [20] CAO F L, MENG X Y, WU J T,et al. Reference correlation of the viscosity of meta-xylene from 273 to 673 K and up to 200 MPa[J].Journal of Physical and Chemical Reference Data, 2016, 45(1):013103.

    [21] VOGEL E, SPAN R, HERRMANN S. Reference correlation for the viscosity of ethane[J]. Journal of Physical and Chemical Reference Data, 2015, 44(4): 043101.

    [22] BALOGUN B, RIESCO N, VESOVIC V. Reference correlation of the viscosity of para-xylene from the triple point to 673 K and up to 110 MPa[J]. Journal of Physical and Chemical Reference Data, 2015,44(1): 013103.

    [23] AVGERI S, ASSAEL M J, HUBER M L,et al. Reference correlation of the viscosity of toluene from the triple point to 675 K and up to 500 MPa[J]. Journal of Physical and Chemical Reference Data, 2015,44(3): 033101.

    [24] AVGERI S, ASSAEL M J, HUBER M L,et al. Reference correlation of the viscosity of benzene from the triple point to 675 K and up to 300 MPa[J]. Journal of Physical and Chemical Reference Data, 2014,43(3): 033103.

    [25] CHUNG T H, AJLAN M, LEE L L,et al. Generalized multiparameter correlation for nonpolar and polar fluid transport properties[J].Industrial and Engineering Chemistry Research, 1988, 27(4):671-679.

    [26] TITANI T. The viscosity of vapours of organic compounds(PartⅢ)[J]. Bulletin of the Chemical Society of Japan, 1933, 8(9): 255-276.

    [27] GOLUBEV I F, VASIL'KOVSKAYA T N, ZOLIN V S. Density of propyl and isopropyl alcohols at various temperatures and pressures[J]. Trudy GIAP, 1979, 54(5): 5-15.

    [28] RAINWATER J C, FRIEND D G. Second viscosity and thermal-conductivity virial coefficients of gases: extension to low reduced temperature[J]. Physical Review A, 1987, 36(8): 4062.

    [29] RAINWATER J C. On the phase space subdivision of the second virial coefficient and its consequences for kinetic theory[J]. The Journal of Chemical Physics, 1984, 81(1): 495-510.

    [30] FRIEND D G, RAINWATER J C. Transport properties of a moderately dense gas[J]. Chemical Physics Letters, 1984, 107(6):590-594.

    [31] VOGEL E, KUECHENMEISTER C, BICH E,et al. Reference correlation of the viscosity of propane[J]. Journal of Physical and Chemical Reference Data, 1998, 27(5): 947-970.

    [32] POLING B E, PRAUSNITZ J M, O'CONNELL J P. The properties of gases and liquids[M]. 5th ed. New York: McGraw-Hill, 2001.

    [33] 張建侯, 劉洪勤. 有機(jī)物質(zhì)粘度的加壓測(cè)定和高壓液體粘度與密度的關(guān)聯(lián)方程[J]. 化工學(xué)報(bào), 1991, 42(3): 269-277.ZHANG J H, LIU H Q. Determination of viscosity under pressure for organic substances and its correlation with density under high pressure[J]. Journal of Chemical Industry and Engineering(China),1991, 42(3): 269-277.

    [34] MOKHTARANI B, SHARIFI A, MORTAHEB H R,et al. Density and viscosity of 1-butyl-3-methylimidaz-olium nitrate with ethanol,1-propanol, or 1-butanol at several temperatures[J]. The Journal of Chemical Thermodynamics, 2009, 41(12): 1432-1438.

    [35] KERMANPOUR F, NIAKAN H, SHARIFI T. Density and viscosity measurements of binary alkanol mixtures from (293.15 to 333.15) K at atmospheric pressure[J]. Journal of Chemical and Engineering Data,2013, 58(5): 1086-1091.

    [36] MIKHAIL S Z, KIMEL W R. Densities and viscosities of 1-propanol-water mixtures[J]. Journal of Chemical and Engineering Data, 1963, 8(3): 323-328.

    [37] RAUF M A, STEWART G H, FARHATAZIZ. Viscosities and densities of binary mixtures of 1-alkanols from 15 to 55℃[J]. Journal of Chemical and Engineering Data, 1983, 28(3): 324-328.

    [38] PAEZ S, CONTRERAS M. Densities and viscosities of binary mixtures of 1-propanol and 2-propanol with acetonitrile[J]. Journal of Chemical and Engineering Data, 1989, 34(4): 455-459.

    [39] GARCIA B, HERRERA C, LEAL J M. Shear viscosities of binary liquid mixtures: 2-pyrrolidinone with 1-alkanols[J]. Journal of Chemical and Engineering Data, 1991, 36(3): 269-274.

    [40] SALEH M A, UDDIN M H. Viscosity and thermodynamics of viscous flow of the systems, 1-propanol + formamide, +N-methylformamide, +N,N-dimethylformamide[J]. Physics and Chemistry of Liquids, 1999, 37(6): 701-718.

    [41] SALEH M A, BEGUM S, UDDIN M H. Viscosity and thermodynamics of viscous flow of 1-propanol with aniline,N-methylaniline andN,N-dimethylaniline[J]. Physics and Chemistry of Liquids, 2001, 39(4): 465-480.

    [42] BHUIYAN M, FERDAUSH J, UDDIN M. Densities and viscosities of binary mixtures of {dimethylsulfoxide+ aliphatic lower alkanols(C1—C3)} at temperatures fromT=303.15 K toT=323.15 K[J]. The Journal of Chemical Thermodynamics, 2007, 39(5): 675-683.

    [43] NAIN A K. Ultrasonic and viscometric studies of molecular interactions in binary mixtures of formamide with ethanol, 1-propanol,1, 2-ethanediol and 1, 2-propanediol at different temperatures[J].Journal of Molecular Liquids, 2008, 140(1): 108-116.

    [44] DOMANSKA U, LASKOWSKA M. Temperature and composition dependence of the density and viscosity of binary mixtures of{1-butyl-3-methylimidazolium thiocyanate+ 1-alcohols}[J]. Journal of Chemical and Engineering Data, 2009, 54(7): 2113-2119.

    [45] RANJBAR S, FAKHRI K, GHASEMI J B. Densities and viscosities of (1-propanol + 1,2-dichloroethane), (1-propan-ol + benzaldehyde),(benzaldehyde + 1,2-dichloroethane), and (1-propanol +1,2-dich-loroethane + benzaldehyde) mixtures fromT=288.15 K to 313.15 K[J]. Journal of Chemical and Engineering Data, 2009, 5(12):3284-3290.

    [46] YAN J H, DAI L Y, WANG X Z,et al. Densities and viscosities of binary mixtures of cyclopropanecarboxylic acid with methanol,ethanol, propan-1-ol, and butan-1-ol at different temperatures[J].Journal of Chemical and Engineering Data, 2009, 54(3):1147-1152.

    [47] HAN S, ZHU C, MA Y. Densities and viscosities of naphthalen-1-ol,naphthalen-2-ol, and 1-aminonaphthalene in the solvents of ethanol,propan-1-ol, and butan-1-ol[J]. Journal of Chemical and Engineering Data, 2012, 57(5): 1396-1403.

    [48] QIAN W, XU Y, ZHU H,et al. Properties of pure 1-methylimidazolium acetate ionic liquid and its binary mixtures with alcohols[J]. The Journal of Chemical Thermodynamics, 2012, 49(3):87-94.

    [49] XU Y, YAO J, WANG C,et al. Density, viscosity, and refractive index properties for the binary mixtures ofn-butylammonium acetate ionic liquid+ alkanols at several temperatures[J]. Journal of Chemical and Engineering Data, 2012, 57(2): 298-308.

    [50] ESTRADA-BALTAZAR A, IGLESIAS-SILVA G, CABALLEROCERóN C. Volumetric and transport properties of binary mixtures ofn-octane + ethanol, + 1-propanol,+ 1-butanol, and + 1-pentanol from(293.15 to 323.15) K at atmospheric pressure[J]. Journal of Chemical and Engineering Data, 2013, 58(12): 3351-3363.

    [51] KIJEV?ANIN M L, ?IVKOVI? E M, DJORDJEVI? B D,et al.Experimental determination and modeling of excess molar volumes,viscosities and refractive indices of the binary systems(pyridine+1-propanol, + 1,2-propanediol, + 1,3-propanediol, and +glycerol). New UNIFAC-VISCO parameters determination[J]. The Journal of Chemical Thermodynamics, 2013, 56(1): 49-56.

    [52] KERMANPOUR F, SHARIFI T. Measurement and correlation of the excess properties of ternary mixture of {x1[Hmim][BF4] +x21-propanol +x32-propanol} at different temperatures[J]. Journal of Chemical and Engineering Data, 2014, 59(6): 1922-1929.

    [53] ZHU C, HAN S, LIU J,et al. Densities and viscosities of diaminotoluene with water, ethanol, propan-1-ol, and butan-1-ol from(293.15 to 333.15) K[J]. Journal of Chemical and Engineering Data,2014, 59(3): 880-889.

    date:2017-03-30.

    MENG Xianyang, associate professor,xymeng@ mail.xjtu.edu.cn

    supported by the National Natural Science Foundation of China (51676159) and the Natural Science Basic Research Plan in Shaanxi Province of China (2015JM5214).

    Correlation of viscosity ofn-propanol

    MENG Xianyang, YU Ye, WU Jiangtao
    (Key Laboratory of Thermo-Fluid Science and Engineering of Ministry of Education,School of Energy and Power Engineering,Xi’an Jiaotong University,Xi’an710049,Shaanxi,China)

    n-Propanol is widely used in different industries such as pharmaceutical manufacturing and agriculture.The thermophysical properties are the fundamentals for those processes. Viscosity is an important property which influences the design of process equipment, product quality monitoring and so on. However, the experimental data of viscosity are often obtained with discrete points. Correlating equations with these experimental data are prerequisite to use. To author’s best knowledge, there is no such a viscosity correlation forn-propanol available in wide temperature and pressure ranges. The aim of this work is to develop the viscosity correlation ofn-propanol with high accuracy, which can be used over the whole thermodynamic surface. Based on the selected reliable experimental data, the correlation for the viscosity ofn-propanol was presented. The extrapolation behavior of the correlation is reasonable, and the temperature is valid from 150 to 400 K at pressures up to 120 MPa. The estimated uncertainties of this correlation are from 1% to 5%.

    n-propanol; viscosity; model; thermodynamics

    TK 121

    A

    0438—1157(2017)11—4035—08

    10.11949/j.issn.0438-1157.20170328

    2017-03-30收到初稿,2017-07-17收到修改稿。

    聯(lián)系人及第一作者:孟現(xiàn)陽(yáng)(1978—),男,副教授。

    國(guó)家自然科學(xué)基金項(xiàng)目(51676159);陜西省自然科學(xué)基礎(chǔ)研究計(jì)劃項(xiàng)目(2015JM5214)。

    猜你喜歡
    正丙醇偏差黏度
    清香型白酒醛雜味貢獻(xiàn)成分定量分析研究
    釀酒科技(2023年11期)2023-12-19 12:43:18
    氣相色譜法測(cè)定白酒中正丙醇含量的研究
    釀酒科技(2023年9期)2023-11-25 04:58:46
    苯-正丙醇-三辛基甲基醋酸銨物系等壓氣液相平衡數(shù)據(jù)的測(cè)定及關(guān)聯(lián)
    石油化工(2020年8期)2020-09-15 06:33:14
    如何走出文章立意偏差的誤區(qū)
    兩矩形上的全偏差
    超高黏度改性瀝青的研發(fā)與性能評(píng)價(jià)
    上海公路(2019年3期)2019-11-25 07:39:30
    響應(yīng)面法優(yōu)化蠟質(zhì)芽孢桿菌低產(chǎn)正丙醇發(fā)酵工藝
    水的黏度的分子動(dòng)力學(xué)模擬
    關(guān)于均數(shù)與偏差
    SAE J300新規(guī)格增加了SAE 8和SAE 12兩種黏度級(jí)別
    久久97久久精品| 久久久久国产网址| 国产精品福利在线免费观看| 2018国产大陆天天弄谢| 久久国内精品自在自线图片| 在线观看人妻少妇| 午夜福利在线观看免费完整高清在| 少妇被粗大猛烈的视频| 亚洲精品久久午夜乱码| 亚洲综合精品二区| 精品一区二区三区视频在线| 国产精品熟女久久久久浪| 国产在视频线精品| 亚洲国产色片| 青春草国产在线视频| 久久这里有精品视频免费| 日本黄大片高清| 美女国产视频在线观看| 免费观看的影片在线观看| 日日啪夜夜爽| av.在线天堂| 亚洲国产精品专区欧美| av卡一久久| 欧美日韩亚洲高清精品| 三级经典国产精品| 色吧在线观看| 一本一本综合久久| av专区在线播放| av.在线天堂| 我的女老师完整版在线观看| 97热精品久久久久久| 亚洲人成网站在线播| 看免费成人av毛片| 亚洲欧美成人综合另类久久久| av国产精品久久久久影院| 国产无遮挡羞羞视频在线观看| 欧美成人一区二区免费高清观看| 久久精品国产亚洲av天美| 亚洲在久久综合| 日本与韩国留学比较| 国产国拍精品亚洲av在线观看| 色婷婷av一区二区三区视频| 亚洲色图综合在线观看| 国产91av在线免费观看| 大香蕉久久网| 少妇被粗大猛烈的视频| 欧美成人一区二区免费高清观看| 亚洲人成网站高清观看| 我的老师免费观看完整版| 我的女老师完整版在线观看| 精品一区二区三区视频在线| 尾随美女入室| 日韩伦理黄色片| 免费观看无遮挡的男女| 亚洲精品一区蜜桃| 身体一侧抽搐| 久久久久人妻精品一区果冻| 热re99久久精品国产66热6| 尾随美女入室| 国产精品av视频在线免费观看| 亚洲欧美日韩卡通动漫| av黄色大香蕉| 天堂中文最新版在线下载| 免费不卡的大黄色大毛片视频在线观看| 最后的刺客免费高清国语| 国产亚洲av片在线观看秒播厂| 欧美最新免费一区二区三区| 国产精品一及| 欧美精品亚洲一区二区| 99国产精品免费福利视频| 精品人妻视频免费看| 久久综合国产亚洲精品| 久久久久久久国产电影| 插逼视频在线观看| 精品久久国产蜜桃| 妹子高潮喷水视频| 久久女婷五月综合色啪小说| 在线观看三级黄色| 91精品伊人久久大香线蕉| 秋霞在线观看毛片| 又爽又黄a免费视频| 三级经典国产精品| 18禁裸乳无遮挡动漫免费视频| 高清欧美精品videossex| 亚洲人成网站高清观看| 欧美亚洲 丝袜 人妻 在线| 一级二级三级毛片免费看| 欧美3d第一页| 99re6热这里在线精品视频| 人妻夜夜爽99麻豆av| 91久久精品电影网| 看非洲黑人一级黄片| 欧美一区二区亚洲| 熟妇人妻不卡中文字幕| 欧美xxⅹ黑人| 久久ye,这里只有精品| 午夜激情福利司机影院| 日韩,欧美,国产一区二区三区| 最黄视频免费看| 国产成人a∨麻豆精品| 亚洲图色成人| 身体一侧抽搐| 我要看日韩黄色一级片| 亚洲人成网站在线播| 精品国产一区二区三区久久久樱花 | 亚洲欧美精品专区久久| 一区二区三区免费毛片| 久久婷婷青草| 亚洲熟女精品中文字幕| 一区二区三区免费毛片| 久久婷婷青草| 久久国产乱子免费精品| 国产成人免费无遮挡视频| 欧美成人午夜免费资源| 男的添女的下面高潮视频| 国产亚洲午夜精品一区二区久久| 国产亚洲一区二区精品| 美女脱内裤让男人舔精品视频| 久久久久视频综合| 国产亚洲一区二区精品| 国产白丝娇喘喷水9色精品| 亚洲av男天堂| 亚洲精品日本国产第一区| 国产在线男女| 男女下面进入的视频免费午夜| 97热精品久久久久久| 日韩视频在线欧美| 国产 一区精品| 国产免费视频播放在线视频| 精品亚洲成a人片在线观看 | 蜜桃久久精品国产亚洲av| 最近中文字幕2019免费版| 精品午夜福利在线看| 狂野欧美白嫩少妇大欣赏| 久久久久久久久久成人| 九九久久精品国产亚洲av麻豆| 国产亚洲精品久久久com| 80岁老熟妇乱子伦牲交| 国产真实伦视频高清在线观看| 多毛熟女@视频| 春色校园在线视频观看| 欧美3d第一页| 亚洲av成人精品一区久久| 国产在线视频一区二区| 成人毛片a级毛片在线播放| av免费在线看不卡| 欧美zozozo另类| 建设人人有责人人尽责人人享有的 | 秋霞伦理黄片| 成人国产麻豆网| 18禁裸乳无遮挡免费网站照片| 日韩电影二区| 在线观看一区二区三区| 成年人午夜在线观看视频| 97超碰精品成人国产| 欧美bdsm另类| 国产精品99久久99久久久不卡 | 午夜免费男女啪啪视频观看| 美女主播在线视频| 亚洲国产精品成人久久小说| www.av在线官网国产| 一本—道久久a久久精品蜜桃钙片| 国产 精品1| 亚洲在久久综合| 国产精品一及| 久久人人爽人人片av| 色综合色国产| av一本久久久久| 在线亚洲精品国产二区图片欧美 | 一区二区三区四区激情视频| 久久精品熟女亚洲av麻豆精品| av播播在线观看一区| 18禁在线无遮挡免费观看视频| 男人添女人高潮全过程视频| 搡女人真爽免费视频火全软件| 在现免费观看毛片| 人妻一区二区av| 亚洲怡红院男人天堂| 日韩亚洲欧美综合| 少妇高潮的动态图| 激情 狠狠 欧美| 老熟女久久久| 自拍偷自拍亚洲精品老妇| 91精品国产国语对白视频| 少妇精品久久久久久久| 国产一区二区三区av在线| 午夜福利影视在线免费观看| 日本黄色日本黄色录像| 一级a做视频免费观看| 久热久热在线精品观看| 制服丝袜香蕉在线| 色视频在线一区二区三区| 丝袜喷水一区| 欧美bdsm另类| 在线免费观看不下载黄p国产| 黄色配什么色好看| av又黄又爽大尺度在线免费看| 亚洲精品成人av观看孕妇| 看十八女毛片水多多多| 亚洲国产精品一区三区| 国产色婷婷99| 色婷婷av一区二区三区视频| 建设人人有责人人尽责人人享有的 | www.色视频.com| 国产 精品1| 国产精品一区二区三区四区免费观看| 久久影院123| 亚洲精品一区蜜桃| 免费观看a级毛片全部| 久久久色成人| 欧美激情极品国产一区二区三区 | 伊人久久精品亚洲午夜| 免费av中文字幕在线| 中文字幕亚洲精品专区| 免费高清在线观看视频在线观看| 日韩欧美精品免费久久| 久久久久国产网址| 国产精品精品国产色婷婷| 精品人妻一区二区三区麻豆| 亚洲天堂av无毛| 成年人午夜在线观看视频| 黄色一级大片看看| 深爱激情五月婷婷| 亚洲丝袜综合中文字幕| 男女下面进入的视频免费午夜| 性色avwww在线观看| 精品久久久噜噜| 亚洲欧美精品自产自拍| 国产免费视频播放在线视频| 久热这里只有精品99| 草草在线视频免费看| 久久人人爽人人片av| 99视频精品全部免费 在线| 97在线视频观看| 亚洲精品456在线播放app| 亚洲最大成人中文| 3wmmmm亚洲av在线观看| 国精品久久久久久国模美| 精品人妻一区二区三区麻豆| 日韩中字成人| 成年免费大片在线观看| 91久久精品电影网| 国产精品久久久久成人av| 免费看光身美女| 免费久久久久久久精品成人欧美视频 | 日本黄色日本黄色录像| 舔av片在线| 国产一区有黄有色的免费视频| 日本色播在线视频| 日韩一区二区三区影片| 成人毛片a级毛片在线播放| 六月丁香七月| 黄色配什么色好看| 国产精品女同一区二区软件| 亚洲av中文字字幕乱码综合| 国产亚洲欧美精品永久| 精品亚洲成国产av| 精品一区二区免费观看| 欧美日韩亚洲高清精品| 亚洲欧美日韩无卡精品| 99热6这里只有精品| 精品熟女少妇av免费看| 亚洲高清免费不卡视频| 精品酒店卫生间| 丰满乱子伦码专区| 日韩电影二区| 亚洲av福利一区| 欧美极品一区二区三区四区| 日本色播在线视频| 久久人人爽人人片av| 久久人妻熟女aⅴ| 亚洲,一卡二卡三卡| 久久久欧美国产精品| 各种免费的搞黄视频| 午夜精品国产一区二区电影| 亚洲欧美精品专区久久| 99热这里只有是精品在线观看| 男女免费视频国产| 99久久精品国产国产毛片| www.av在线官网国产| 青春草亚洲视频在线观看| 中文字幕久久专区| 国产欧美另类精品又又久久亚洲欧美| 亚洲国产精品国产精品| 一区二区av电影网| 久久久久久久大尺度免费视频| 亚洲av中文av极速乱| 国产精品嫩草影院av在线观看| 久久久久久久久久久免费av| 卡戴珊不雅视频在线播放| 日日摸夜夜添夜夜添av毛片| 狂野欧美激情性xxxx在线观看| 亚洲av成人精品一区久久| 草草在线视频免费看| 男人狂女人下面高潮的视频| 国产午夜精品久久久久久一区二区三区| 人妻少妇偷人精品九色| 成人毛片a级毛片在线播放| 最近手机中文字幕大全| 国产高清不卡午夜福利| 久久午夜福利片| 国产一区二区三区综合在线观看 | 99精国产麻豆久久婷婷| 久久午夜福利片| 久久鲁丝午夜福利片| 成人影院久久| 亚洲欧美成人精品一区二区| 成人特级av手机在线观看| 夜夜爽夜夜爽视频| 免费不卡的大黄色大毛片视频在线观看| 九九爱精品视频在线观看| 一区在线观看完整版| 少妇熟女欧美另类| 婷婷色麻豆天堂久久| 男人爽女人下面视频在线观看| 国产av精品麻豆| 国产成人a∨麻豆精品| 久久鲁丝午夜福利片| 免费大片18禁| 国产伦在线观看视频一区| 免费av中文字幕在线| 久久国产精品男人的天堂亚洲 | 免费看不卡的av| 国产亚洲欧美精品永久| 亚洲丝袜综合中文字幕| 夫妻午夜视频| 国产伦精品一区二区三区四那| xxx大片免费视频| 乱码一卡2卡4卡精品| 亚洲国产av新网站| 91久久精品国产一区二区三区| 一级爰片在线观看| 国产男人的电影天堂91| 免费少妇av软件| 我的女老师完整版在线观看| av在线老鸭窝| 777米奇影视久久| 国产 一区 欧美 日韩| 国产精品秋霞免费鲁丝片| 亚洲国产成人一精品久久久| 日韩精品有码人妻一区| 免费看日本二区| 街头女战士在线观看网站| 少妇 在线观看| 久久热精品热| 狂野欧美白嫩少妇大欣赏| 男的添女的下面高潮视频| 久久久久网色| 天天躁夜夜躁狠狠久久av| 国产淫片久久久久久久久| 国产男女超爽视频在线观看| 丰满少妇做爰视频| 熟女电影av网| 老司机影院成人| 在线看a的网站| 大香蕉97超碰在线| 边亲边吃奶的免费视频| 国产精品99久久99久久久不卡 | 制服丝袜香蕉在线| 在线免费十八禁| 91久久精品电影网| 51国产日韩欧美| av网站免费在线观看视频| 在线 av 中文字幕| 欧美精品人与动牲交sv欧美| 狂野欧美白嫩少妇大欣赏| 男女啪啪激烈高潮av片| 三级国产精品片| av在线观看视频网站免费| 亚洲精品第二区| 日日啪夜夜撸| 干丝袜人妻中文字幕| 国产精品成人在线| 极品教师在线视频| 亚洲欧洲国产日韩| 国产亚洲av片在线观看秒播厂| av免费在线看不卡| 免费少妇av软件| 777米奇影视久久| 亚洲三级黄色毛片| 欧美国产精品一级二级三级 | 国产精品国产三级专区第一集| 亚洲精品第二区| 午夜福利在线在线| 干丝袜人妻中文字幕| 成人黄色视频免费在线看| 美女视频免费永久观看网站| 亚洲欧美成人综合另类久久久| 国产久久久一区二区三区| 日韩精品有码人妻一区| 中文在线观看免费www的网站| 一本色道久久久久久精品综合| 乱码一卡2卡4卡精品| 免费不卡的大黄色大毛片视频在线观看| 又黄又爽又刺激的免费视频.| 亚洲国产欧美在线一区| av专区在线播放| 97在线视频观看| 欧美高清成人免费视频www| 久久久久人妻精品一区果冻| 国产在线视频一区二区| 91精品伊人久久大香线蕉| 老司机影院毛片| 日韩免费高清中文字幕av| 永久网站在线| 51国产日韩欧美| 欧美亚洲 丝袜 人妻 在线| 欧美3d第一页| 搡女人真爽免费视频火全软件| 亚洲av在线观看美女高潮| 亚洲av男天堂| 一级毛片aaaaaa免费看小| 欧美国产精品一级二级三级 | 亚洲精品aⅴ在线观看| 嫩草影院新地址| 18+在线观看网站| 精品久久国产蜜桃| 国产伦在线观看视频一区| 日本爱情动作片www.在线观看| 五月开心婷婷网| 男男h啪啪无遮挡| av福利片在线观看| 国产老妇伦熟女老妇高清| 国内少妇人妻偷人精品xxx网站| 久久人人爽人人片av| 国产一区二区三区av在线| 亚洲国产成人一精品久久久| 色吧在线观看| 成人午夜精彩视频在线观看| 久久久久久久精品精品| 色网站视频免费| 国产精品熟女久久久久浪| 国产乱人视频| 亚洲av综合色区一区| 少妇被粗大猛烈的视频| 国产 一区 欧美 日韩| 高清av免费在线| 亚洲精华国产精华液的使用体验| av国产精品久久久久影院| 欧美最新免费一区二区三区| 久久久久网色| 国产精品成人在线| 亚洲图色成人| 2018国产大陆天天弄谢| 国产伦在线观看视频一区| 国产免费视频播放在线视频| 欧美xxⅹ黑人| 久久久久性生活片| 日韩精品有码人妻一区| 日本欧美视频一区| 久久久久人妻精品一区果冻| 国产精品.久久久| 美女cb高潮喷水在线观看| 久久热精品热| 晚上一个人看的免费电影| 99热网站在线观看| 久热这里只有精品99| 国产综合精华液| 国产白丝娇喘喷水9色精品| 国产亚洲精品久久久com| 91午夜精品亚洲一区二区三区| 热re99久久精品国产66热6| 日韩精品有码人妻一区| 日韩视频在线欧美| 在线观看美女被高潮喷水网站| 看十八女毛片水多多多| 精品一区二区三区视频在线| 91aial.com中文字幕在线观看| 亚洲美女搞黄在线观看| 我要看日韩黄色一级片| 成人高潮视频无遮挡免费网站| 美女内射精品一级片tv| 欧美精品人与动牲交sv欧美| 国产精品久久久久久久久免| 18禁裸乳无遮挡动漫免费视频| 国产在线男女| 亚洲精品国产成人久久av| 久久久精品免费免费高清| 男女边吃奶边做爰视频| h视频一区二区三区| 在线观看三级黄色| 最近最新中文字幕大全电影3| 亚洲伊人久久精品综合| 欧美极品一区二区三区四区| 哪个播放器可以免费观看大片| 国产亚洲一区二区精品| 婷婷色麻豆天堂久久| 国产黄色视频一区二区在线观看| 天天躁夜夜躁狠狠久久av| 国产美女午夜福利| 永久免费av网站大全| 毛片女人毛片| 精品国产一区二区三区久久久樱花 | 嫩草影院新地址| 99国产精品免费福利视频| 18禁在线播放成人免费| 日韩成人伦理影院| 激情 狠狠 欧美| 蜜桃在线观看..| 国产亚洲最大av| 亚洲人成网站在线观看播放| 777米奇影视久久| 精品久久久噜噜| 亚洲av福利一区| 欧美97在线视频| 在线观看三级黄色| 国产免费又黄又爽又色| xxx大片免费视频| 亚洲三级黄色毛片| 免费久久久久久久精品成人欧美视频 | 久久人妻熟女aⅴ| 人妻系列 视频| 久久99热这里只有精品18| 久热这里只有精品99| 国精品久久久久久国模美| 国产成人精品久久久久久| 国产精品一二三区在线看| 99热网站在线观看| 国产精品不卡视频一区二区| 大片免费播放器 马上看| 在线观看av片永久免费下载| 欧美高清性xxxxhd video| 欧美bdsm另类| 国产欧美日韩精品一区二区| 91精品国产国语对白视频| 亚洲精品aⅴ在线观看| 特大巨黑吊av在线直播| 国产男人的电影天堂91| 看免费成人av毛片| 久久久久精品性色| 日韩伦理黄色片| 一级毛片aaaaaa免费看小| 久久6这里有精品| 在线看a的网站| 蜜桃亚洲精品一区二区三区| 国产精品人妻久久久影院| 久热这里只有精品99| 男男h啪啪无遮挡| 久久久亚洲精品成人影院| 亚洲精品乱码久久久v下载方式| 国产精品欧美亚洲77777| 国产精品一区二区三区四区免费观看| 国产 精品1| 搡女人真爽免费视频火全软件| 亚洲av中文字字幕乱码综合| 18禁在线无遮挡免费观看视频| 亚洲一区二区三区欧美精品| 新久久久久国产一级毛片| 亚洲内射少妇av| 女性生殖器流出的白浆| 97精品久久久久久久久久精品| 国产精品秋霞免费鲁丝片| 国产成人午夜福利电影在线观看| 性色avwww在线观看| av播播在线观看一区| 一级爰片在线观看| 亚洲色图综合在线观看| 久久热精品热| 又爽又黄a免费视频| 成人18禁高潮啪啪吃奶动态图 | 夜夜爽夜夜爽视频| 久久6这里有精品| 天天躁夜夜躁狠狠久久av| 国产在线一区二区三区精| 成人18禁高潮啪啪吃奶动态图 | 一级毛片aaaaaa免费看小| 91精品伊人久久大香线蕉| 午夜日本视频在线| 国产 一区 欧美 日韩| 综合色丁香网| 91久久精品电影网| 99久久人妻综合| 99九九线精品视频在线观看视频| 亚洲电影在线观看av| 欧美 日韩 精品 国产| 夫妻午夜视频| 99精国产麻豆久久婷婷| 亚洲av综合色区一区| 九色成人免费人妻av| 免费av不卡在线播放| 久久人人爽人人片av| 久久久久视频综合| 国产在线免费精品| 久久人人爽人人片av| 免费不卡的大黄色大毛片视频在线观看| 伊人久久精品亚洲午夜| 国产精品99久久99久久久不卡 | 精品国产三级普通话版| 97超碰精品成人国产| 毛片女人毛片| 在线观看一区二区三区| 日韩在线高清观看一区二区三区| 亚洲av免费高清在线观看| 男女下面进入的视频免费午夜| 欧美高清成人免费视频www| av不卡在线播放| 日韩精品有码人妻一区| 久久亚洲国产成人精品v| 国内精品宾馆在线| 精品国产三级普通话版| 亚洲无线观看免费| 亚洲国产精品成人久久小说| 亚洲综合色惰| 日本黄色日本黄色录像| 国产中年淑女户外野战色| 中国三级夫妇交换| 黄色一级大片看看| 麻豆乱淫一区二区| 欧美精品国产亚洲| 国产白丝娇喘喷水9色精品| 性色av一级| 一级a做视频免费观看| 国产有黄有色有爽视频| 香蕉精品网在线| 成人影院久久|