• <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í)別
    欧美精品高潮呻吟av久久| 亚洲一码二码三码区别大吗| 电影成人av| 一区二区三区精品91| 国产欧美日韩综合在线一区二区| 亚洲熟女精品中文字幕| 国产男人的电影天堂91| 亚洲欧洲日产国产| 丝袜美足系列| 欧美日韩黄片免| 18禁观看日本| 青春草视频在线免费观看| 国产又色又爽无遮挡免| 久久人妻熟女aⅴ| 日韩一本色道免费dvd| 欧美变态另类bdsm刘玥| 色网站视频免费| 美女大奶头黄色视频| 欧美日韩福利视频一区二区| 多毛熟女@视频| 日日爽夜夜爽网站| 侵犯人妻中文字幕一二三四区| 天天躁夜夜躁狠狠躁躁| 久久亚洲国产成人精品v| 捣出白浆h1v1| 亚洲精品久久久久久婷婷小说| 国产成人91sexporn| 亚洲精品国产色婷婷电影| 91老司机精品| 亚洲欧美精品综合一区二区三区| 97在线人人人人妻| 国产亚洲欧美在线一区二区| 如日韩欧美国产精品一区二区三区| 爱豆传媒免费全集在线观看| 99久久综合免费| 97精品久久久久久久久久精品| 欧美亚洲日本最大视频资源| 亚洲人成网站在线观看播放| 欧美激情 高清一区二区三区| 成在线人永久免费视频| 汤姆久久久久久久影院中文字幕| 国产亚洲欧美精品永久| 国产成人精品久久二区二区免费| 国产一级毛片在线| 免费黄频网站在线观看国产| 国产一区二区激情短视频 | 久久国产精品影院| 国产一区有黄有色的免费视频| 亚洲av美国av| 男女之事视频高清在线观看 | 无限看片的www在线观看| 亚洲成av片中文字幕在线观看| 又粗又硬又长又爽又黄的视频| 青草久久国产| 女人高潮潮喷娇喘18禁视频| 久久精品国产综合久久久| 黄色片一级片一级黄色片| 国产成人精品久久久久久| 丰满迷人的少妇在线观看| 黄色 视频免费看| 精品高清国产在线一区| 18在线观看网站| 午夜福利乱码中文字幕| 秋霞在线观看毛片| 欧美 亚洲 国产 日韩一| 日本猛色少妇xxxxx猛交久久| 亚洲国产av影院在线观看| 亚洲精品久久久久久婷婷小说| 久久人妻熟女aⅴ| 国产成人一区二区在线| 亚洲伊人色综图| 超碰成人久久| 久久久久久久久免费视频了| 又大又黄又爽视频免费| 国产福利在线免费观看视频| 久久久久精品人妻al黑| 老熟女久久久| 中文字幕亚洲精品专区| 99九九在线精品视频| 国产99久久九九免费精品| 在线av久久热| 中文欧美无线码| 亚洲精品久久成人aⅴ小说| 黄色视频在线播放观看不卡| 欧美久久黑人一区二区| 大陆偷拍与自拍| 午夜福利,免费看| 最近中文字幕2019免费版| 精品亚洲乱码少妇综合久久| 精品国产国语对白av| 久久综合国产亚洲精品| 欧美黄色淫秽网站| 在线观看免费日韩欧美大片| 国产精品一国产av| 国产国语露脸激情在线看| 国产亚洲欧美精品永久| www.精华液| 中国国产av一级| 国产免费视频播放在线视频| av国产久精品久网站免费入址| 国产精品免费视频内射| 啦啦啦中文免费视频观看日本| 午夜久久久在线观看| 色婷婷av一区二区三区视频| 一本—道久久a久久精品蜜桃钙片| 日本wwww免费看| 50天的宝宝边吃奶边哭怎么回事| 黄色a级毛片大全视频| 国产精品一区二区在线不卡| 免费看十八禁软件| 热99久久久久精品小说推荐| 超碰成人久久| 每晚都被弄得嗷嗷叫到高潮| 女人爽到高潮嗷嗷叫在线视频| 亚洲国产最新在线播放| 精品久久久精品久久久| 别揉我奶头~嗯~啊~动态视频 | 国产极品粉嫩免费观看在线| 又粗又硬又长又爽又黄的视频| a级毛片在线看网站| 叶爱在线成人免费视频播放| 黄频高清免费视频| 日韩熟女老妇一区二区性免费视频| 91麻豆精品激情在线观看国产 | 成人亚洲精品一区在线观看| 丁香六月天网| 男人爽女人下面视频在线观看| 中文字幕av电影在线播放| 日韩大片免费观看网站| 亚洲天堂av无毛| 性色av乱码一区二区三区2| 亚洲精品久久成人aⅴ小说| 成年人午夜在线观看视频| 精品高清国产在线一区| 夫妻午夜视频| 最新在线观看一区二区三区 | 精品免费久久久久久久清纯 | h视频一区二区三区| 日韩精品免费视频一区二区三区| av在线老鸭窝| 91成人精品电影| www.自偷自拍.com| 伊人久久大香线蕉亚洲五| 女人被躁到高潮嗷嗷叫费观| 无遮挡黄片免费观看| 无遮挡黄片免费观看| 亚洲国产毛片av蜜桃av| 亚洲欧美色中文字幕在线| 美女午夜性视频免费| 久久亚洲精品不卡| 午夜免费男女啪啪视频观看| 高清视频免费观看一区二区| 精品少妇黑人巨大在线播放| 黄片小视频在线播放| 日韩精品免费视频一区二区三区| 波多野结衣av一区二区av| 午夜福利乱码中文字幕| 国产精品香港三级国产av潘金莲 | 亚洲国产av新网站| 菩萨蛮人人尽说江南好唐韦庄| 乱人伦中国视频| 在线 av 中文字幕| 天天躁日日躁夜夜躁夜夜| 国产精品久久久久成人av| 亚洲自偷自拍图片 自拍| 美女福利国产在线| av在线播放精品| 免费日韩欧美在线观看| 又黄又粗又硬又大视频| 欧美国产精品va在线观看不卡| 黄色视频在线播放观看不卡| 午夜视频精品福利| 18禁观看日本| 在线亚洲精品国产二区图片欧美| 黑人巨大精品欧美一区二区蜜桃| 亚洲伊人色综图| 女人高潮潮喷娇喘18禁视频| 欧美+亚洲+日韩+国产| 国产亚洲av片在线观看秒播厂| 国产精品香港三级国产av潘金莲 | 亚洲中文av在线| 国产高清videossex| 色视频在线一区二区三区| 99国产精品一区二区三区| 侵犯人妻中文字幕一二三四区| 国产老妇伦熟女老妇高清| 九色亚洲精品在线播放| 9191精品国产免费久久| 亚洲一区中文字幕在线| 久久久久网色| 在线亚洲精品国产二区图片欧美| 久久久久久久国产电影| av在线app专区| 亚洲成人免费av在线播放| 在线观看免费视频网站a站| 日韩视频在线欧美| 天天躁狠狠躁夜夜躁狠狠躁| 国产一卡二卡三卡精品| 19禁男女啪啪无遮挡网站| 国产伦理片在线播放av一区| 伊人久久大香线蕉亚洲五| 老熟女久久久| 男女床上黄色一级片免费看| 久久人妻福利社区极品人妻图片 | 9191精品国产免费久久| 欧美日韩精品网址| 2021少妇久久久久久久久久久| 少妇人妻久久综合中文| 高潮久久久久久久久久久不卡| 一边摸一边做爽爽视频免费| 赤兔流量卡办理| 麻豆av在线久日| 久久影院123| 精品少妇一区二区三区视频日本电影| 亚洲欧美色中文字幕在线| 国产精品久久久久久人妻精品电影 | 国产一区二区三区综合在线观看| 最新在线观看一区二区三区 | 欧美在线一区亚洲| 一级毛片女人18水好多 | 国产片内射在线| 国产在线一区二区三区精| 久久精品成人免费网站| 在线 av 中文字幕| 丝袜美腿诱惑在线| 免费一级毛片在线播放高清视频 | 精品人妻一区二区三区麻豆| 欧美人与性动交α欧美精品济南到| 国产片特级美女逼逼视频| 亚洲,欧美,日韩| 婷婷丁香在线五月| 国产亚洲精品第一综合不卡| 国产欧美日韩综合在线一区二区| 激情视频va一区二区三区| 亚洲欧美一区二区三区黑人| av国产久精品久网站免费入址| 熟女av电影| 久久久久视频综合| av天堂在线播放| 国产精品成人在线| 中文乱码字字幕精品一区二区三区| www.av在线官网国产| 国产精品欧美亚洲77777| 人人妻,人人澡人人爽秒播 | 亚洲成人免费电影在线观看 | 亚洲欧美激情在线| 伊人亚洲综合成人网| av国产久精品久网站免费入址| 久久av网站| 一级片'在线观看视频| 91字幕亚洲| 国产成人影院久久av| 免费在线观看视频国产中文字幕亚洲 | 国产伦理片在线播放av一区| 黑人猛操日本美女一级片| 精品亚洲成a人片在线观看| 欧美精品一区二区大全| 精品久久久精品久久久| 久久女婷五月综合色啪小说| 国精品久久久久久国模美| 国产精品国产av在线观看| 国产97色在线日韩免费| 人人妻,人人澡人人爽秒播 | 在线观看www视频免费| 大型av网站在线播放| 热re99久久国产66热| 亚洲第一青青草原| 国产视频首页在线观看| 天天影视国产精品| 日本vs欧美在线观看视频| 亚洲精品中文字幕在线视频| 亚洲av国产av综合av卡| 手机成人av网站| 人体艺术视频欧美日本| 亚洲国产欧美网| av有码第一页| 欧美另类一区| 午夜福利乱码中文字幕| 精品一区二区三区四区五区乱码 | 亚洲av电影在线进入| 人人澡人人妻人| 天天影视国产精品| 各种免费的搞黄视频| 国产爽快片一区二区三区| 在线观看免费日韩欧美大片| 欧美成人午夜精品| 亚洲av国产av综合av卡| 9热在线视频观看99| 叶爱在线成人免费视频播放| 亚洲欧美一区二区三区久久| 老司机影院毛片| 国产99久久九九免费精品| 18禁国产床啪视频网站| 国产高清国产精品国产三级| 国产精品国产av在线观看| 精品久久蜜臀av无| 97在线人人人人妻| 菩萨蛮人人尽说江南好唐韦庄| 丝袜在线中文字幕| 午夜av观看不卡| 又紧又爽又黄一区二区| 中国美女看黄片| 国产片特级美女逼逼视频| 香蕉丝袜av| 丝袜美足系列| 午夜福利在线免费观看网站| 精品卡一卡二卡四卡免费| 在线天堂中文资源库| 久久热在线av| 另类亚洲欧美激情| 伦理电影免费视频| 久久久精品免费免费高清| xxx大片免费视频| 精品人妻一区二区三区麻豆| 免费观看a级毛片全部| 欧美日韩国产mv在线观看视频| 欧美变态另类bdsm刘玥| 婷婷丁香在线五月| 欧美日韩亚洲高清精品| 久久天躁狠狠躁夜夜2o2o | 国产黄频视频在线观看| 手机成人av网站| 夫妻午夜视频| 两个人免费观看高清视频| 啦啦啦在线观看免费高清www| 丁香六月天网| 亚洲精品乱久久久久久| 国产精品久久久av美女十八| 777米奇影视久久| 无遮挡黄片免费观看| 少妇人妻久久综合中文| 精品亚洲成a人片在线观看| 国产成人av教育| 国产欧美日韩一区二区三区在线| 女警被强在线播放| 老熟女久久久| 桃花免费在线播放| 少妇精品久久久久久久| 亚洲综合色网址| 国产精品免费视频内射| 国产91精品成人一区二区三区 | 日本wwww免费看| 纯流量卡能插随身wifi吗| 亚洲色图 男人天堂 中文字幕| 黄网站色视频无遮挡免费观看| 婷婷色麻豆天堂久久| 亚洲专区国产一区二区| 久久99精品国语久久久| 婷婷成人精品国产| 精品熟女少妇八av免费久了| 久久久久久久大尺度免费视频| 亚洲熟女精品中文字幕| 亚洲精品国产区一区二| 中文字幕人妻熟女乱码| 日韩制服丝袜自拍偷拍| 精品一区二区三区四区五区乱码 | 国产一区二区三区av在线| 激情视频va一区二区三区| 大香蕉久久成人网| 免费在线观看影片大全网站 | 国产深夜福利视频在线观看| 日本vs欧美在线观看视频| 伊人久久大香线蕉亚洲五| 久久精品久久久久久噜噜老黄| 久久久久久人人人人人| 日本wwww免费看| av线在线观看网站| 成人国语在线视频| 亚洲欧美色中文字幕在线| 狂野欧美激情性xxxx| 丁香六月天网| 一级片'在线观看视频| 男女无遮挡免费网站观看| 可以免费在线观看a视频的电影网站| 在线观看人妻少妇| 免费不卡黄色视频| 亚洲成人免费av在线播放| 日韩大片免费观看网站| 精品一区二区三卡| 成人国产av品久久久| 精品国产国语对白av| 大陆偷拍与自拍| 一本综合久久免费| cao死你这个sao货| 一级黄片播放器| 国产精品国产三级专区第一集| 午夜两性在线视频| 69精品国产乱码久久久| 免费在线观看影片大全网站 | 亚洲成人免费电影在线观看 | 国产亚洲av片在线观看秒播厂| 国产精品亚洲av一区麻豆| 搡老岳熟女国产| 大片电影免费在线观看免费| 国产欧美日韩一区二区三 | 日日爽夜夜爽网站| 成人亚洲欧美一区二区av| 最近手机中文字幕大全| 性色av乱码一区二区三区2| 在现免费观看毛片| 久久精品成人免费网站| 欧美亚洲日本最大视频资源| 亚洲国产av新网站| 国产在线观看jvid| 男女高潮啪啪啪动态图| 一二三四社区在线视频社区8| 一区二区三区四区激情视频| 成在线人永久免费视频| av不卡在线播放| 精品免费久久久久久久清纯 | 欧美人与性动交α欧美精品济南到| 女人被躁到高潮嗷嗷叫费观| 性少妇av在线| 男女无遮挡免费网站观看| www.熟女人妻精品国产| 老汉色av国产亚洲站长工具| 中文字幕高清在线视频| 99国产综合亚洲精品| 国产精品熟女久久久久浪| 亚洲黑人精品在线| 激情视频va一区二区三区| 久9热在线精品视频| 成人国产一区最新在线观看 | 久久天堂一区二区三区四区| 免费在线观看日本一区| 人体艺术视频欧美日本| 97人妻天天添夜夜摸| 黄片小视频在线播放| 丰满迷人的少妇在线观看| 欧美人与善性xxx| 久久国产亚洲av麻豆专区| 国产免费现黄频在线看| 五月天丁香电影| 色精品久久人妻99蜜桃| 免费在线观看影片大全网站 | 老汉色∧v一级毛片| 亚洲欧洲国产日韩| 国产免费又黄又爽又色| 欧美大码av| 在线天堂中文资源库| 99久久99久久久精品蜜桃| 亚洲欧美一区二区三区黑人| 麻豆乱淫一区二区| 人人妻人人澡人人看| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲精品国产色婷婷电影| 99国产精品免费福利视频| 亚洲欧美日韩高清在线视频 | 精品国产超薄肉色丝袜足j| 超色免费av| 精品国产超薄肉色丝袜足j| 久久99精品国语久久久| 国产精品国产三级国产专区5o| 国产在线免费精品| 在线观看免费高清a一片| 国产免费现黄频在线看| 大片免费播放器 马上看| 波野结衣二区三区在线| 丰满少妇做爰视频| 免费高清在线观看日韩| 国语对白做爰xxxⅹ性视频网站| 黑丝袜美女国产一区| 人妻人人澡人人爽人人| 香蕉丝袜av| 亚洲自偷自拍图片 自拍| 男女免费视频国产| 亚洲自偷自拍图片 自拍| 美女脱内裤让男人舔精品视频| 脱女人内裤的视频| 19禁男女啪啪无遮挡网站| 在线观看免费高清a一片| 性高湖久久久久久久久免费观看| 看免费av毛片| 久久久国产精品麻豆| 啦啦啦中文免费视频观看日本| 人人妻人人爽人人添夜夜欢视频| 久久久精品94久久精品| 男女无遮挡免费网站观看| 日本猛色少妇xxxxx猛交久久| h视频一区二区三区| 色综合欧美亚洲国产小说| 精品视频人人做人人爽| 免费久久久久久久精品成人欧美视频| 新久久久久国产一级毛片| 亚洲国产av影院在线观看| 亚洲精品一区蜜桃| 日韩伦理黄色片| 18禁观看日本| 欧美日韩一级在线毛片| 每晚都被弄得嗷嗷叫到高潮| 欧美日韩视频精品一区| 国产日韩欧美视频二区| 欧美精品一区二区免费开放| 麻豆乱淫一区二区| 一区二区三区精品91| 国产日韩欧美视频二区| 国产精品偷伦视频观看了| 久热爱精品视频在线9| 日韩大码丰满熟妇| 超碰97精品在线观看| 一区在线观看完整版| 另类精品久久| 婷婷色综合大香蕉| 嫩草影视91久久| 精品久久久久久电影网| 丝瓜视频免费看黄片| 日本wwww免费看| 99久久人妻综合| 亚洲国产av影院在线观看| 人人妻,人人澡人人爽秒播 | 国语对白做爰xxxⅹ性视频网站| 中文字幕另类日韩欧美亚洲嫩草| 欧美乱码精品一区二区三区| 国产三级黄色录像| 美女福利国产在线| 久久女婷五月综合色啪小说| 99九九在线精品视频| 午夜福利在线免费观看网站| 99九九在线精品视频| 精品久久久久久久毛片微露脸 | 伦理电影免费视频| 考比视频在线观看| 久久久久国产一级毛片高清牌| 免费少妇av软件| 欧美少妇被猛烈插入视频| 亚洲视频免费观看视频| 搡老乐熟女国产| 人人妻人人爽人人添夜夜欢视频| 一级毛片女人18水好多 | 色视频在线一区二区三区| 另类亚洲欧美激情| 制服人妻中文乱码| 免费在线观看完整版高清| 久久久久精品国产欧美久久久 | 91成人精品电影| av片东京热男人的天堂| 国产无遮挡羞羞视频在线观看| 婷婷色麻豆天堂久久| 亚洲欧美一区二区三区国产| 又大又爽又粗| 国产精品国产av在线观看| 欧美大码av| www.999成人在线观看| 中文字幕av电影在线播放| 人人妻人人澡人人爽人人夜夜| 国产免费一区二区三区四区乱码| 亚洲九九香蕉| 免费日韩欧美在线观看| 国产成人精品久久久久久| 女性被躁到高潮视频| 亚洲七黄色美女视频| 国产日韩欧美视频二区| 首页视频小说图片口味搜索 | 99九九在线精品视频| 免费高清在线观看视频在线观看| 91成人精品电影| 国产成人av激情在线播放| 中文字幕另类日韩欧美亚洲嫩草| 亚洲精品中文字幕在线视频| 大片电影免费在线观看免费| 五月开心婷婷网| 一级a爱视频在线免费观看| 精品久久久久久电影网| 国产伦理片在线播放av一区| 不卡av一区二区三区| 成年人黄色毛片网站| 亚洲欧美色中文字幕在线| 免费在线观看日本一区| av网站免费在线观看视频| 亚洲专区中文字幕在线| 精品一区在线观看国产| 国产有黄有色有爽视频| 久久精品久久久久久久性| 中文精品一卡2卡3卡4更新| 美女高潮到喷水免费观看| 亚洲专区国产一区二区| 久久久久久亚洲精品国产蜜桃av| 亚洲精品美女久久久久99蜜臀 | 成年人黄色毛片网站| 国产精品久久久av美女十八| 高清不卡的av网站| 免费看十八禁软件| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲精品av麻豆狂野| 日韩av在线免费看完整版不卡| 丁香六月天网| av国产精品久久久久影院| 久久免费观看电影| 亚洲伊人色综图| 天堂8中文在线网| 欧美精品一区二区大全| 女人高潮潮喷娇喘18禁视频| 高潮久久久久久久久久久不卡| 色网站视频免费| 亚洲欧美一区二区三区国产| 成人黄色视频免费在线看| 十八禁高潮呻吟视频| 亚洲av成人精品一二三区| 少妇精品久久久久久久| 久久久精品94久久精品| www.999成人在线观看| 亚洲国产欧美网| 日本a在线网址| 日韩av不卡免费在线播放| 久久久久久久大尺度免费视频| 久久毛片免费看一区二区三区| 尾随美女入室| av网站在线播放免费| 国产老妇伦熟女老妇高清| 人成视频在线观看免费观看| 最近手机中文字幕大全| 最近中文字幕2019免费版| 一区二区三区乱码不卡18| 久久精品国产亚洲av涩爱| 一级毛片女人18水好多 |