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

    纈沙坦在乙酸乙酯中的溶解度及介穩(wěn)區(qū)

    2016-08-18 06:36:00王清清孫勤楊阿三程榕鄭燕萍浙江工業(yè)大學(xué)化學(xué)工程學(xué)院浙江杭州310014
    化工進(jìn)展 2016年8期
    關(guān)鍵詞:傳熱

    王清清,孫勤,楊阿三,程榕,鄭燕萍(浙江工業(yè)大學(xué)化學(xué)工程學(xué)院,浙江 杭州 310014)

    孫斌,曲藝,楊迪(東北電力大學(xué)能源與動(dòng)力工程學(xué)院,吉林 吉林 132012)

    研究開發(fā)

    纈沙坦在乙酸乙酯中的溶解度及介穩(wěn)區(qū)

    王清清,孫勤,楊阿三,程榕,鄭燕萍
    (浙江工業(yè)大學(xué)化學(xué)工程學(xué)院,浙江 杭州 310014)

    采用靜態(tài)平衡法測(cè)定纈沙坦在乙酸乙酯中278.15~323.15K的溶解度數(shù)據(jù),利用Modified Apelblat、NRTL 和λh方程分別對(duì)溶解度數(shù)據(jù)進(jìn)行了關(guān)聯(lián),并通過van't Hoff分析計(jì)算了纈沙坦在溶解過程中的熱力學(xué)參數(shù)(焓變、熵變和吉布斯自由能)。結(jié)果表明,纈沙坦在乙酸乙酯中的溶解度隨著溫度升高而顯著增大,3個(gè)模型均具有較好的關(guān)聯(lián)性,平均相對(duì)偏差分別為1.03%、3.87%、1.72%,Modified Apelblat方程對(duì)溶解度數(shù)據(jù)關(guān)聯(lián)的效果最好;纈沙坦在乙酸乙酯中的溶解過程為自發(fā)吸熱熵驅(qū)動(dòng)過程,焓變?cè)谌芙膺^程中對(duì)吉布斯自由能貢獻(xiàn)較大。利用激光法測(cè)定了纈沙坦在乙酸乙酯中的超溶解度,研究了不同攪拌速率、降溫速率對(duì)結(jié)晶介穩(wěn)區(qū)的影響。結(jié)果表明,隨飽和溫度的升高,介穩(wěn)區(qū)顯著變寬;攪拌速率越小、降溫速率越快,介穩(wěn)區(qū)越寬,降溫速率的影響相對(duì)較小。

    纈沙坦;乙酸乙酯;溶解度;超溶解度;介穩(wěn)區(qū)

    纈沙坦,白色晶體,化學(xué)名為 N-(1-氧戊基)-N-[[2'-(1H-四唑-5-基)[1,1'-聯(lián)苯]-4-基]甲基]-L-纈氨酸,分子式 C24H29N5O3,相對(duì)分子質(zhì)量435.53,熔點(diǎn)383.15K。纈沙坦屬于非肽類、口服有效的血管緊張素Ⅱ(AT)受體拮抗劑,選擇性作用于Ⅰ型受體(AT1),用于治療各種類型高血壓,具有療效顯著、耐藥性好和成本效益低的優(yōu)點(diǎn),對(duì)需要長(zhǎng)期治療高血壓的患者來說是一個(gè)很好的選擇[1-2]。

    乙酸乙酯是纈沙坦工業(yè)結(jié)晶的首選溶劑,研究結(jié)晶體系的介穩(wěn)區(qū)性質(zhì)對(duì)優(yōu)化結(jié)晶過程和結(jié)晶器設(shè)計(jì)起到重要的作用[3],特別是超溶解度數(shù)據(jù)對(duì)工業(yè)化的實(shí)際生產(chǎn)有著更為重要的意義。纈沙坦在乙酸乙酯中的溶解度已有文獻(xiàn)報(bào)道[4],但還未發(fā)現(xiàn)有關(guān)纈沙坦在乙酸乙酯中超溶解度數(shù)據(jù)的報(bào)道。

    本文通過激光法測(cè)定纈沙坦在乙酸乙酯中的超溶解度,采用靜態(tài)平衡法測(cè)定纈沙坦在乙酸乙酯溶劑中278.15~323.15K的溶解度,研究攪拌速率和降溫速率對(duì)介穩(wěn)區(qū)寬度的影響。利用Modified Apelblat、NRTL和λh方程對(duì)溶解度數(shù)據(jù)進(jìn)行交叉關(guān)聯(lián)驗(yàn)證,并采用van't Hoff分析,計(jì)算溶解過程中的熱力學(xué)參數(shù)(焓變、熵變和吉布斯自由能)。

    1 實(shí) 驗(yàn)

    1.1原料和儀器

    纈沙坦,純度>99%,浙江華海藥業(yè)股份有限公司;乙酸乙酯,分析純,上海凌峰化學(xué)試劑有限公司;溴百里香酚藍(lán)指示劑,灼燒殘?jiān)?.3%,國(guó)藥集團(tuán)化學(xué)試劑有限公司;氦氖激光器,SPL-HN1.5P型,工作波長(zhǎng)632.8nm,杭州譜鐳光電技術(shù)有限公司;低溫恒溫水槽,DC-3010型,寧波天恒儀器廠。

    溶解度測(cè)定實(shí)驗(yàn)裝置由低溫恒溫水槽、250mL三口燒瓶、精密溫度計(jì)和機(jī)械攪拌器組成。超溶解度測(cè)定實(shí)驗(yàn)裝置由自制夾套冷卻結(jié)晶器、低溫恒溫水槽、激光發(fā)射器、激光接收器、精密溫度計(jì)和機(jī)械攪拌器組成。

    1.2纈沙坦溶解度測(cè)定

    纈沙坦在乙酸乙酯中的溶解度采用靜態(tài)平衡法測(cè)定:置已知量的乙酸乙酯于250mL三口燒瓶中,加入過量纈沙坦,在恒定溫度下,充分?jǐn)嚢韬箝L(zhǎng)時(shí)間靜置。移取上層清液精密稱重后加入適量乙醇溶劑,再加入溴百里香酚藍(lán)指示劑,用配好的氫氧化鈉滴定液(0.05mol/L)滴定至藍(lán)色,并將結(jié)果用空白試驗(yàn)進(jìn)行校正。1mL氫氧化鈉滴定液(0.1mol/L)相對(duì)于21.78mgC24H29N5O3(來源于中國(guó)藥典)。依次測(cè)定278.15~323.15K下纈沙坦的溶解度。溶液溫度由恒溫水浴控制,每個(gè)點(diǎn)重復(fù)取樣 3次取平均值。

    纈沙坦在溶劑中的摩爾分?jǐn)?shù)x,用式(1)來計(jì)算。

    式中,mA、mB分別為溶液中纈沙坦和乙酸乙酯的質(zhì)量;MA、MB分別為纈沙坦和乙酸乙酯相對(duì)分子質(zhì)量。

    1.3纈沙坦超溶解度測(cè)定

    采用激光法測(cè)定超溶解度存在一定的滯后性,但由于纈沙坦在乙酸乙酯中的結(jié)晶速率快且實(shí)驗(yàn)中控制的降溫速率較小,因此激光法的滯后性對(duì)實(shí)驗(yàn)結(jié)果的影響較小。

    根據(jù)溶解度數(shù)據(jù),配置成一定溫度下的纈沙坦飽和溶液置于結(jié)晶器中。在攪拌速率恒定時(shí),以一定的降溫速率進(jìn)行冷卻。激光發(fā)射器發(fā)射激光穿透飽和溶液到達(dá)另一端的激光接收器。當(dāng)接收激光的強(qiáng)度發(fā)生突變時(shí),即有晶體產(chǎn)生,記錄下此時(shí)的溫度。則所配的飽和溶液濃度即為纈沙坦此溫度下的超溶解度,與初始飽和溶液的溫度差ΔT即為介穩(wěn)區(qū)寬度。

    2 結(jié)果與討論

    2.1纈沙坦在乙酸乙酯中的溶解度

    雖然LIU等[4]早期采用激光法測(cè)定了纈沙坦在乙酸乙酯中278.15~313.15K的溶解度數(shù)據(jù),但是經(jīng)實(shí)驗(yàn)發(fā)現(xiàn)利用接收器端激光強(qiáng)度達(dá)到最大值來判斷實(shí)驗(yàn)終點(diǎn)一般會(huì)產(chǎn)生延后,導(dǎo)致加入的溶劑質(zhì)量偏大,這樣測(cè)定的溶解度數(shù)據(jù)相比實(shí)際值會(huì)偏小。而且在文中只采用Modified Apelblat對(duì)實(shí)驗(yàn)數(shù)據(jù)進(jìn)行了擬合,得到的均方根偏差(rmsd)為6.0×10-5,計(jì)算其平均相對(duì)偏差(MD)為1.21%,單個(gè)方程的擬合結(jié)果不足以論證實(shí)驗(yàn)的準(zhǔn)確性。所以,本文采用靜態(tài)平衡法重新測(cè)定了纈沙坦在乙酸乙酯溶劑中278.15~323.15K的溶解度。

    圖 1是實(shí)驗(yàn)值與文獻(xiàn)值[4]的對(duì)比結(jié)果。從圖1中可以發(fā)現(xiàn),纈沙坦在乙酸乙酯中的溶解度隨溫度的升高而增大,但增大的幅度不相同,在308.15~323.15K時(shí)變化迅速。實(shí)驗(yàn)值基本與文獻(xiàn)值吻合,但略高于文獻(xiàn)值。

    圖1 纈沙坦在乙酸乙酯中溶解度實(shí)驗(yàn)值與文獻(xiàn)值比較

    2.2溶解度不同模型關(guān)聯(lián)結(jié)果

    纈沙坦溶解度分別采用 Modified Apelblat、NRTL和λh方程進(jìn)行關(guān)聯(lián),其擬合效果可以用平均偏差(σ)及平均相對(duì)偏差(MD)來表示,其定義如式(2)~式(3)。

    式中,xexp、xcal分別代表實(shí)驗(yàn)值和理論值;N為實(shí)驗(yàn)點(diǎn)數(shù)。

    2.2.1Modified Apelblat 方程

    Modified Apelblat 方程廣泛用于溶解度的關(guān)聯(lián)[5-8],如式(4)所示。

    式中,A、B、C為Apelblat關(guān)聯(lián)參數(shù);T為熱力學(xué)溫度,K。

    2.2.2NRTL方程

    要研究NRTL方程對(duì)溶解度的關(guān)聯(lián),首先通過熱力學(xué)理論和假設(shè),可以得到式(5)。

    式中,溶劑的摩爾分?jǐn)?shù) x2=1- x1,g12-g22,

    g21-g11和α1(2α12=α2)1為NRTL參數(shù)。

    2.2.3λh方程

    λh方程也經(jīng)常用于溶解度的關(guān)聯(lián),方程具體形式如式(9)[12]。

    式中,λ和h為方程參數(shù)。

    表1~表3分別為上述3種方程回歸得到的模型參數(shù)及平均相對(duì)偏差,表4為實(shí)驗(yàn)值與模型計(jì)算值的對(duì)比。

    從表中可以看出,3個(gè)方程平均相對(duì)偏差≤3.87%,關(guān)聯(lián)較好。其中Modified Apelblat方程關(guān)聯(lián)效果最好,平均相對(duì)偏差只有1.03%。

    農(nóng)業(yè)技術(shù)推廣人員的專業(yè)能力將會(huì)直接決定我國(guó)農(nóng)業(yè)技術(shù)推廣的效果,在這樣的情況下,要更好的提升農(nóng)業(yè)技術(shù)推廣的作用,對(duì)于農(nóng)業(yè)技術(shù)推廣人員的專業(yè)能力進(jìn)行提升是非常有必要的。政府的相關(guān)部門可以通過定期對(duì)農(nóng)業(yè)技術(shù)推廣人員進(jìn)行培訓(xùn),組織技術(shù)人員下鄉(xiāng)與弄濕實(shí)際溝通等方式,提升其專業(yè)能力。在此基礎(chǔ)上還可以通過編制調(diào)整,更好的引入人才,為農(nóng)業(yè)技術(shù)推廣的開展提供充分支持。

    表1 Modified Apelblat方程回歸的參數(shù)值及平均相對(duì)偏差

    表2 NRTL方程回歸的參數(shù)值及平均相對(duì)偏差

    表3 λh方程回歸的參數(shù)值及平均相對(duì)偏差

    表4 纈沙坦的溶解度數(shù)據(jù)及Modified Apelbla、NRTL和λh方程回歸結(jié)果

    2.3溶解熱力學(xué)計(jì)算

    式中,Tmean為平均熱力學(xué)溫度,本研究中Tmean=299.96K。

    此外,為了比較溶解過程中焓變和熵變對(duì)吉布斯自由能的貢獻(xiàn)大小,引入了ξH和ξS,定義如式(15)、式(16)[15]。

    表5 纈沙坦溶解過程中的熱力學(xué)參數(shù)

    溶解過程是一個(gè)熵驅(qū)動(dòng)過程。ξH>0.58,說明溶解過程中焓變對(duì)吉布斯自由能的貢獻(xiàn)較大,為焓控制過程。

    2.4纈沙坦在乙酸乙酯中結(jié)晶介穩(wěn)區(qū)分析

    介穩(wěn)區(qū)是介于溶解度曲線和超溶解度曲線之間的區(qū)域。介穩(wěn)區(qū)寬度可以用極限過飽和度 ΔCmax或極限溫度過冷度ΔTmax來表示,本研究采用ΔTmax表示。

    對(duì)介穩(wěn)區(qū)寬度的影響因素有很多,如降溫速率、攪拌速率、有無晶種、有無外加場(chǎng)及溶液的性質(zhì)等。實(shí)驗(yàn)中采用控制變量法考察攪拌速率和降溫速率對(duì)介穩(wěn)區(qū)寬度的影響。

    2.4.1攪拌速率對(duì)介穩(wěn)區(qū)的影響

    保持降溫速率為0.2K/min恒定,分別設(shè)定攪拌速率為100r/min、200r/min和300r/min,考察攪拌速率對(duì)介穩(wěn)區(qū)寬度的影響。

    從圖2中可以看出,纈沙坦在乙酸乙酯中的介穩(wěn)區(qū)寬度隨著攪拌速率的增大而變窄,而隨飽和溫度的升高,攪拌速率對(duì)介穩(wěn)區(qū)的影響增大。這可能由于攪拌速率的變大使得纈沙坦分子碰撞的機(jī)會(huì)增多,更利于晶體成核,從而使介穩(wěn)區(qū)變窄。此外,從圖中還知,飽和溫度越高,介穩(wěn)區(qū)越寬。在低溫區(qū)(293.15~313.15K)ΔTmax為 3~8K,高溫區(qū)(313.15~323.15K)ΔTmax為8~12K。

    2.4.2降溫速率對(duì)介穩(wěn)區(qū)的影響

    保持?jǐn)嚢杷俾蕿?00r/min恒定,分別設(shè)定降溫速率為0.2K/min、0.3K/min和0.4K/min,考察降溫速率對(duì)介穩(wěn)區(qū)寬度的影響。

    從圖3可知,纈沙坦在乙酸乙酯中的介穩(wěn)區(qū)寬度隨著降溫速率的增大而變寬。這是由于降溫速率的增大使得溶質(zhì)經(jīng)過成核溫度區(qū)域過短,不利于緩慢穩(wěn)定成核,從而使介穩(wěn)區(qū)有所變寬,但影響并不明顯。

    圖2 不同攪拌速率下纈沙坦在乙酸乙酯中的介穩(wěn)區(qū)

    圖3 不同降溫速率下纈沙坦在乙酸乙酯中的介穩(wěn)區(qū)

    3 結(jié) 論

    (1)纈沙坦在乙酸乙酯中278.15~323.15K下的溶解度隨著溫度的升高而顯著增大,采用Modified Apelblat、NRTL和λh方程進(jìn)行擬合,均具有良好的關(guān)聯(lián)性,其中Modified Apelblat方程關(guān)聯(lián)效果最好,平均相對(duì)偏差僅有1.03%。

    (2)通過van't Hoff分析得到纈沙坦在乙酸乙酯溶解過程中的熱力學(xué)參數(shù):焓變、熵變和吉布斯自由能均為正值,說明此過程為熵驅(qū)動(dòng)自發(fā)吸熱過程。而ξH>0.58說明焓變?cè)谌芙膺^程中對(duì)吉布斯自由能貢獻(xiàn)較大。

    (3)纈沙坦在乙酸乙酯中的結(jié)晶介穩(wěn)區(qū)寬度隨著飽和溫度的升高、降溫速率的增大而變寬,隨著攪拌速率的增大而變窄,其中飽和溫度影響顯著,而降溫速率的影響不太明顯。

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

    Tmean—— 熱力學(xué)平均溫度,K

    Tmi—— 組分i的熔點(diǎn)溫度,K—— 標(biāo)準(zhǔn)摩爾溶解熵,J/(mol·K)

    參考文獻(xiàn)

    [1]NIXON R M,MULLER E,LOWY A,et al. Valsartan vs. other aongiotensinⅡ receptor blockers in the treatment of hypertension:a meta-analytical approach[J]. Int. J. Clin. Pract.,2009,63(5):766-775.

    [2]NEUTEL J,WEBER M,POOL J,et al. Valsartan,a new angotensin Ⅱ antagonist:antihypertensive effects over 24 hours[J]. Clin. Ther.,1997,19(3):447-458.

    [3]馬勇,朱家文,陳蔡,等. 磷酸結(jié)晶介穩(wěn)區(qū)性質(zhì)的研究[J]. 高?;瘜W(xué)工程學(xué)報(bào),2010,24(2):331-335.

    [4]LIU Y,WANG J K,WANG X B,et al. Solubility of valsartan in ethyl acetate + hexane binary mixtures from(278.15 to 313.15)K[J]. J. Chem. Eng. Data,2009,54(4):1412-1414.

    [5]王小敏,汪全義,彭琦,等. 對(duì)乙酰氨基苯甲酸在甲醇+乙醇混合溶劑中溶解度測(cè)定及關(guān)聯(lián)計(jì)算[J]. 化工進(jìn)展,2015,34(3):647-651.

    [6]TSONOPOULOS C. Thermodynamic analysis of the mutual solubilities of normal alkanes and water[J]. Fluid Phase Equilibr.,1999,156(1):21-33.

    [7]MANZUROLA E,APELBLAT A. Solubilities of L-glutamic acid,3-nitrobenzoic acid,p-toluic acid,calcium-L-lactate,calcium gluconate,magnesium-DL-aspartate,and magnesium-L-lactate in water[J]. J. Chem. Thermodyn.,2002,34(7):127-136.

    [8]HUANG C,XIE Z P,XU J C,et al. Experimental and modeling studies on the solubility of d-pantolactone in four pure solvents and ethanol-water mixtures[J]. J. Chem. Eng. Data,2015,60(3):870-875.

    [9]MARTI E E. Crystalline salt forms of valsartan:US 6869970[P]. 2005-03-22.

    [10]宋微微,馬沛生,相政樂. 丁二酸在五種有機(jī)溶劑中溶解度測(cè)定及關(guān)聯(lián)[J]. 高?;瘜W(xué)工程學(xué)報(bào),2007,21(2):341-344.

    [11]李柏春,杜東雪,張文林,等. 氯化鋰、氯化鈉在NMP中溶解度的測(cè)定與關(guān)聯(lián)[J]. 化工學(xué)報(bào),2014,65(12):4664-4668.

    [12]張迪,肖清貴,張炳燭,等. 氯氧化鉍在鹽酸溶液中溶解度的測(cè)定和關(guān)聯(lián)[J]. 化工學(xué)報(bào),2014,65(6):1987-1992.

    [13]SU J H,QIAN C,LUO N Z,et al. Experimental measurement and modeling of the solubility of biotin in six pure solvents at temperatures from 298.15 K to 333.85 K[J]. J. Chem. Eng. Data,2014,59(11):3894-3899.

    [14]SHAKEEL F, HAQ N,SIDDIQUI N A,et al. Solubility and thermodynamic behavior of vanillin in propane-1,2-diol + water cosolvent mixtures at different temperatures[J]. Food Chem.,2015,188:57-61.

    [15]YU C,ZENG Z X,XUE W L. Measurement and correlation of the solubility of hexaquonickel(Ⅱ) bis(ptoluenesul-fonate) in water + ethanol solvents within 288.15-333.15 K[J]. Ind. Eng. Chem. Res.,2015,54(15):3961-3967.

    研究開發(fā)

    納米流體沖擊射流換熱特性實(shí)驗(yàn)

    孫斌,曲藝,楊迪
    (東北電力大學(xué)能源與動(dòng)力工程學(xué)院,吉林 吉林 132012)

    摘要:以納米流體為工質(zhì)對(duì)沖擊射流冷卻系統(tǒng)的綜合性能進(jìn)行實(shí)驗(yàn),主要研究了添加納米顆粒的納米流體與水在不同流速、不同射流高度等條件下沖擊射流的傳熱效率,同時(shí)也對(duì)不同種類的納米流體的換熱效率進(jìn)行了對(duì)比。結(jié)果表明:對(duì)于添加了納米顆粒的沖擊射流冷卻系統(tǒng),傳熱效率得到顯著提高,但當(dāng)質(zhì)量分?jǐn)?shù)達(dá)到0.5%時(shí),傳熱系數(shù)變化不明顯。對(duì)于不同種類的納米流體:Cu-水、Al2O3-水和Al-水納米流體,其中Cu-水的換熱效率最高,存在一個(gè)特定的射流高度,使傳熱系數(shù)達(dá)到最大值。研究結(jié)果對(duì)設(shè)計(jì)制造輕型緊湊的高效換熱器有實(shí)用的工程價(jià)值。

    關(guān)鍵詞:納米粒子;多相流;傳熱;沖擊射流

    中圖分類號(hào):TK 121文獻(xiàn)標(biāo)志碼:A文章編號(hào):1000-6613(2016)08-2334-08

    DOI:10.16085/j.issn.1000-6613.2016.08.05

    Experimental study on heat transfer characteristics of nanofluids impacted jet

    SUN Bin,QU Yi,YANG Di
    (Energy and Power Engineering Institute,Northeast Dianli University,Jilin 132012,Jilin,China)

    Abstract:In this paper,comprehensive performance of impinging jet cooling system heat exchanger was experimentally studied using nanofluids. The heat transfer efficiencies were compared for nanofluids of different flow rates,jet height and types . The results revealed that heat transfer efficiency significantly increased with the introduction of nanofluids in jet,but,when the mass percentage of nanofluids exceeded of 0.5%,the heat transfer coefficient did not change significantly. For different types of nanofluids:Cu- water,Al2O3-water,Al-water nanofluids,the highest heat transfer efficiency was observed for Cu- water. In addition,there was a particular jet height,where the maximum heat transfer coefficient could be reached. The results would be practically valuable in designing and manufacturing light and efficient heat exchanger.

    Key words:nanoparticles;multiphase flow;heat transfer;impinging jet

    換熱設(shè)備是在航天、電子、化工、原子能、動(dòng)力等眾多領(lǐng)域中廣泛使用的一種通用設(shè)備,其主要作用是實(shí)現(xiàn)物料之間的熱量傳遞[1]。沖擊射流是一種極其有效的局部強(qiáng)化換熱方法。由于沖擊射流在換熱面上形成的流動(dòng)邊界層較薄,沖擊射流的傳熱系數(shù)比常規(guī)的對(duì)流傳熱系數(shù)高幾倍甚至一個(gè)數(shù)量級(jí)[2]。

    以往的研究主要針對(duì)沖擊射流的流速、沖擊距離、射流角度、噴嘴孔徑等的運(yùn)行參數(shù)和結(jié)構(gòu)參數(shù)。BARDBURY[3]對(duì)單孔自由射流的流動(dòng)機(jī)理進(jìn)行了探索,研究了單孔自由射流典型的速度分布形式及軸向速度的衰減規(guī)律,實(shí)驗(yàn)結(jié)果表明,在單孔自由射流的中心區(qū),流體的速度高、波動(dòng)大,而且湍流度非常高,因此具有極大的強(qiáng)化換熱潛力。

    收稿日期:2016-01-11;修改稿日期:2016-03-16。
    基金項(xiàng)目:教育部新世紀(jì)優(yōu)秀人才支持計(jì)劃項(xiàng)目(NCET-12- 0727)。
    第一作者及聯(lián)系人:孫斌(1972—),男,博士,教授,主要研究方向?yàn)槎嘞嗔骼碚摷皯?yīng)用。E-mail sunbin@nedu.edu.cn。

    BAONGA等[4]利用激光誘導(dǎo)熒光技術(shù)和圖像處理在圓形水平加熱盤上進(jìn)行自由液體射流沖擊試驗(yàn)測(cè)得的局部溫度和速度曲線,通過求解逆熱傳導(dǎo)問題建立了徑向分布的局部傳熱系數(shù),結(jié)果表明該噴嘴到熱交換面間距對(duì)換熱影響不大,雷諾數(shù)相同的情況下噴嘴直徑的改變使得Nu數(shù)發(fā)生明顯變化,并且液體流速影響局部傳熱系數(shù)。KATTI等[5]研究了圓形直管噴嘴沖擊噴氣到光滑平整表面的局部發(fā)熱理論,研究了局部傳熱噴射到板間距和雷諾數(shù)的影響。近幾年,冷卻射流被應(yīng)用到精密元件及微小空間散熱的問題中。NAPHON等[6]進(jìn)行了關(guān)于個(gè)人計(jì)算機(jī)中央處理器沖擊射流強(qiáng)化傳熱的實(shí)驗(yàn),利用小型矩形翅片板進(jìn)行噴射液體的冷卻實(shí)驗(yàn),實(shí)驗(yàn)證明,冷卻效果確實(shí)要比常規(guī)冷卻方式要好,但是需要較高的能量消耗。GUO等[7]使用微針噴頭進(jìn)行沖擊射流研究強(qiáng)化流動(dòng)沸騰換熱過程,結(jié)果表明該沸騰傳熱可以通過提高增加的總面積及射流沖擊速度來提高傳熱系數(shù)。當(dāng)噴射速度過高時(shí),對(duì)流換熱占主導(dǎo)地位;當(dāng)沖擊速度過高時(shí),滯流區(qū)的熱傳遞與渦擴(kuò)大了散熱面積并增加了紊流,影響微對(duì)流過程,影響換熱。

    納米流體是以一定的方式和比例在基液(水、酒精等)中添加納米級(jí)金屬或非金屬氧化物粒子,形成的一類具有高熱導(dǎo)率的均勻穩(wěn)定的新型傳熱工質(zhì)[8]。納米流體以優(yōu)異的導(dǎo)熱性與穩(wěn)定性,為沖擊射流的工質(zhì)提供了新的選擇[9]。凌智勇等[10]對(duì)納米流體的黏度進(jìn)行了研究;李強(qiáng)等[11]發(fā)現(xiàn),納米流體的熱導(dǎo)率通常比基液高。采用納米流體代替常規(guī)的流體對(duì)具有高熱量輸出的設(shè)備冷卻有特別的效果。

    本實(shí)驗(yàn)在一個(gè)典型的射流冷卻系統(tǒng)中研究沖擊射流的傳熱特性,通過使用不同的納米顆粒制備的納米流體進(jìn)行實(shí)驗(yàn),對(duì)不同納米流體的熱導(dǎo)率進(jìn)行對(duì)比,得到的一些相關(guān)結(jié)果,實(shí)驗(yàn)數(shù)據(jù)和納米流體的傳熱特性進(jìn)行討論,具有較高的工程實(shí)用意義。

    1 實(shí)驗(yàn)系統(tǒng)及流程

    1.1納米流體制備

    實(shí)驗(yàn)采用Cu、Al和Al2O3納米粒子,其基本物性如表1所示,基液為去離子水,采用十二烷基苯磺酸鈉(SDBS)作為分散劑[12]。采用“兩步法”[13],按照1∶1的比例加入納米流體和分散劑,此時(shí)納米流體穩(wěn)定性最佳[14]。配置質(zhì)量分?jǐn)?shù)分別為0.1%、0.2%、0.3%、0.4%、0.5%的Cu-水、Al2O3-水、Al-水納米流體。隨后在超聲波震蕩儀中震蕩60min,得到穩(wěn)定的納米流體。

    表1 納米流體熱物性

    1.2納米流體穩(wěn)定性分析

    采用透射比法分別對(duì)質(zhì)量分?jǐn)?shù) 0.1%、0.2%、0.3%、0.4%、0.5%的Cu-水、Al-水、Al2O3-水納米流體的穩(wěn)定性進(jìn)行分析。透射比法即分光光度計(jì)法,又稱濁度法,是將分散體系經(jīng)過重力沉降或離心沉降等手段處理后,吸取上層清液,在分光光度計(jì)上測(cè)定其透過率或吸收率。一般來說,穩(wěn)定分散的懸浮液,其上層清液的透過率較低,吸收率較高[15]。

    圖1為重力沉降出來后的納米流體,吸取上層清液,在分光光度計(jì)中進(jìn)行透射率測(cè)定,所用儀器為752型分光光度計(jì),所用波長(zhǎng)為540nm,每隔12h測(cè)量并記錄一次數(shù)據(jù),重復(fù)3~4次,測(cè)定結(jié)果如圖2所示。

    采用上述實(shí)驗(yàn)方法,實(shí)驗(yàn)范圍內(nèi)各質(zhì)量分?jǐn)?shù)的Cu-水、Al-水、Al2O3-水納米流體均具有很強(qiáng)的穩(wěn)定性,在較長(zhǎng)時(shí)間內(nèi)不易沉積,適于后期的實(shí)驗(yàn)研究及應(yīng)用。

    2 實(shí)驗(yàn)系統(tǒng)

    實(shí)驗(yàn)系統(tǒng)如圖3所示,由儲(chǔ)液箱、磁力循環(huán)泵、預(yù)熱段、實(shí)驗(yàn)測(cè)試段、冷卻水箱及數(shù)據(jù)采集部分組成。

    實(shí)驗(yàn)流程:將制備好的納米流體放入儲(chǔ)液箱中,由揚(yáng)程為30m的磁力泵提供循環(huán)動(dòng)力,采用渦輪流量變送器測(cè)量流體體積流量。納米流體經(jīng)預(yù)熱段預(yù)熱為恒定25℃后進(jìn)入實(shí)驗(yàn)測(cè)試段進(jìn)行換熱實(shí)驗(yàn)。將4個(gè)熱電阻分別布置在射流入口處、射流出口處和紫銅圓柱中心及二分之一半徑處,從而測(cè)量流體入口溫度、流體出口溫度以及紫銅圓柱的溫度。兩個(gè)功率為 50W 的陶瓷加熱芯對(duì)稱放置在紫銅圓柱二分之一半徑處,用來提供熱量。為降低測(cè)量誤差,紫銅柱體內(nèi)部的陶瓷加熱芯以及熱電阻都涂抹導(dǎo)熱膏減小其傳熱熱阻。紫銅柱體外包裹石棉進(jìn)行保溫,紫銅換熱面進(jìn)行鏡面打磨處理。Rosemount 3051S電容式差壓變送器連接在進(jìn)出口位置用來測(cè)量流體進(jìn)出口兩端的壓力差。壓差計(jì)、流量計(jì)和熱電阻所采集的數(shù)據(jù)通過 USB-4716數(shù)據(jù)采集器輸入計(jì)算機(jī),進(jìn)行實(shí)時(shí)檢測(cè)記錄。

    圖1 3種納米流體的重力沉降觀測(cè)照片

    圖4為實(shí)驗(yàn)測(cè)試段結(jié)構(gòu)示意圖。實(shí)驗(yàn)段采用法蘭盤與長(zhǎng)徑螺栓進(jìn)行連接固定,橡膠墊圈進(jìn)行防水密封,易于拆卸組裝。實(shí)驗(yàn)射流距離分別設(shè)定為3mm、6mm、9mm進(jìn)行對(duì)比實(shí)驗(yàn)。射流進(jìn)口管外側(cè)加工出螺紋,通過螺紋調(diào)整射流距離,既保證了密封性能,也易于實(shí)驗(yàn)參數(shù)的變更。實(shí)驗(yàn)段采用有機(jī)玻璃設(shè)計(jì),可以調(diào)整射流距離以及觀察實(shí)驗(yàn)中流體流動(dòng)情況。

    3 數(shù)據(jù)處理與不確定性分析

    采用牛頓(Newton)冷卻公式,流體流動(dòng)的傳熱系數(shù)可寫為式(1)。

    流體的平均溫度采用式(2)計(jì)算。

    式中,Tin是射流實(shí)驗(yàn)段納米流體進(jìn)口溫度;Tout射流實(shí)驗(yàn)段納米流體出口溫度。由于熱電阻測(cè)點(diǎn)與換熱面具有一定距離,可能導(dǎo)致熱電阻測(cè)溫面與實(shí)際的射流沖擊換熱表面溫度不一致,而且這種溫度不一致性隨著熱流密度的增大而增大。為了減少實(shí)驗(yàn)誤差,應(yīng)考慮紫銅塊熱阻的影響如式(3)、式(4)。

    圖2 3種納米流體的穩(wěn)定性

    圖3 實(shí)驗(yàn)系統(tǒng)結(jié)構(gòu)

    圖4 測(cè)試段結(jié)構(gòu)示意圖

    式中,Tw是換熱面的溫度;T'w是熱紫銅塊的實(shí)測(cè)溫度的平均值;ΔT是換熱面溫度與納米流體平均溫度的溫差。

    努塞爾數(shù)Nu定義見式(5)。

    本實(shí)驗(yàn)使用的不同濃度的納米流體熱導(dǎo)率 knf均由DRE-Ш型熱導(dǎo)率測(cè)定儀測(cè)量。雷諾數(shù)定義如式(6)。

    式中,納米流體的密度ρnf定義如式(7)[16]。)

    由于納米流體的體積分?jǐn)?shù)難以精準(zhǔn)測(cè)定,本實(shí)驗(yàn)由納米流體質(zhì)量分?jǐn)?shù)進(jìn)行計(jì)算,如式(8)。

    為了得到更準(zhǔn)確的實(shí)驗(yàn)結(jié)果,需要對(duì)系統(tǒng)進(jìn)行誤差分析及不確定度評(píng)估[17],使用儀器的參數(shù)及不確定度見表2,變量的不確定度見表3。

    4 實(shí)驗(yàn)結(jié)果及分析

    4.1傳熱特性

    在實(shí)驗(yàn)前,對(duì)試驗(yàn)系統(tǒng)進(jìn)行熱平衡計(jì)算,計(jì)算輸入功率Q與測(cè)試段工質(zhì)獲得的熱功率Ql的偏差。熱平衡偏差如式(9)[18]。

    表2 儀器的參數(shù)及不確定度

    表3 實(shí)驗(yàn)變量及其不確定度

    經(jīng)計(jì)算,熱平衡偏差小于 7%,因此工質(zhì)能夠獲得的總功率非常接近100W。

    4.2納米流體沖擊射流換熱性能分析

    實(shí)驗(yàn)中,分別測(cè)試了質(zhì)量分?jǐn)?shù)分別是 0.1%、0.2%、0.3%、0.4%、0.5%的Cu-水、Al2O3-水、Al-水納米流體的射流傳熱系數(shù),分析了納米粒子體積分?jǐn)?shù)、種類、射流速度以及射流高度對(duì)沖擊射流換熱性能的影響。

    4.2.1質(zhì)量分?jǐn)?shù)對(duì)納米流體換熱的影響

    圖5、圖6和圖7分別給出了3種不同種類的納米流體在不同質(zhì)量分?jǐn)?shù)下的沖擊射流傳熱系數(shù)隨Re數(shù)的變化情況。從圖中可以看出,沖擊射流的傳熱系數(shù)隨著 Re數(shù)的增加而增大;另外,在相同的Re數(shù)情況下,沖擊射流傳熱系數(shù)h隨著濃度的增加而增加。這是因?yàn)樵谝后w中添加納米粒子,增加了液體的表面積和熱容量,在液體中添加納米粒子,增加了液體的熱導(dǎo)率,由于粒子與粒子、粒子與液體及粒子與壁面間的相互作用及碰撞,傳熱效果增強(qiáng)。由圖可以發(fā)現(xiàn),3種納米流體在質(zhì)量分?jǐn)?shù)為0.4%時(shí),傳熱系數(shù)提升最為明顯,而質(zhì)量分?jǐn)?shù)為0.5%時(shí),相對(duì)于質(zhì)量分?jǐn)?shù) 0.4%的納米流體的傳熱系數(shù)變化較小。這是因?yàn)榧{米粒子的質(zhì)量分?jǐn)?shù)過大,會(huì)導(dǎo)致納米流體的黏度過大,射流路徑中的能量衰減加快,出口射流速度相同時(shí),沖擊到換熱面的速度減小,對(duì)換熱效果起負(fù)面作用。另外,當(dāng)納米流體質(zhì)量分?jǐn)?shù)過大時(shí),納米顆粒更易發(fā)生沉淀現(xiàn)象,會(huì)有少部分納米顆粒沉淀到儲(chǔ)液箱底部,導(dǎo)致實(shí)際參與換熱的納米流體質(zhì)量分?jǐn)?shù)下降。另外,在 Re數(shù)較小的區(qū)域,傳熱系數(shù)相差不大,因此低流速下,納米流體對(duì)傳熱系數(shù)提高的效果不明顯;而在較高的流速下,傳熱系數(shù)提升顯著。這是因?yàn)樵诹魉佥^高的情況下,紊流度增加,顯著提升了傳熱系數(shù),所以高速射流應(yīng)用納米流體更為有利。

    實(shí)驗(yàn)中使用了較高比例的分散劑,考慮到分散劑對(duì)去離子水換熱性能的影響現(xiàn)作如下實(shí)驗(yàn)分析:將質(zhì)量分?jǐn)?shù)為0.1%、0.2%、0.3%、0.4%和0.5%SDBS分散劑加入去離子水中制備為射流工質(zhì),進(jìn)行沖擊射流實(shí)驗(yàn),如圖8所示。

    圖5 Cu-水納米流體與水傳熱系數(shù)隨Re數(shù)的變化

    圖6 Al2O3-水納米流體與水傳熱系數(shù)隨Re數(shù)的變化

    圖7 Al-水納米流體與水傳熱系數(shù)隨Re數(shù)的變化

    圖8 只添加分散劑SDBS的去離子水的傳熱系數(shù)隨Re數(shù)的變化

    從圖8可以看出,隨著分散劑質(zhì)量分?jǐn)?shù)的增加,該射流工質(zhì)的傳熱系數(shù)隨之減??;同時(shí),隨著射流流速的增加,傳熱系數(shù)呈近似線性增加。對(duì)于只加入分散劑的去離子水為射流工質(zhì),傳熱系數(shù)較去離子水稍有下降,最大降幅 4%。這是因?yàn)榉稚┑臒釋?dǎo)率較低,使得只添加分散劑的去離子水換熱性能略降低。與文獻(xiàn)[19]的研究結(jié)果一致。

    納米流體作為換熱工質(zhì),結(jié)合沖擊射流高效冷卻方式,可以顯著提高換熱效率。主要原因如下:懸浮顆粒增加了兩相混合物的熱導(dǎo)率;納米顆粒加快了納米顆粒、去離子水和換熱壁面三者之間能量的傳遞速率;由于流速加快,單個(gè)納米顆粒強(qiáng)力沖擊換熱面,使換熱流體更好的接觸換熱壁面,明顯降低了換熱壁面上流體層的厚度。此外,本實(shí)驗(yàn)使用的是圓形漸縮噴嘴,更利于沖擊射流的強(qiáng)化換熱。由于上述原因,導(dǎo)致本實(shí)驗(yàn)中傳熱系數(shù)增加較大,質(zhì)量分?jǐn)?shù)為0.5%的Cu-水納米流體增幅最大,傳熱系數(shù)最大增加54%。文獻(xiàn)[20]也給出了相近的結(jié)果。

    4.2.2不同納米流體對(duì)傳熱系數(shù)的影響

    圖9描述的是質(zhì)量分?jǐn)?shù)為0.4%的Cu-水、Al2O3-水和 Al-水納米流體在不同 Re數(shù)下的射流傳熱系數(shù)。由圖9中實(shí)驗(yàn)數(shù)據(jù)可以看出,添加了納米顆粒的納米流體,傳熱系數(shù)明顯高于去離子水。其中,Cu-水納米流體的傳熱系數(shù)相對(duì)于其他兩種納米流體相對(duì)較高,說明 Cu-水納米流體更適宜應(yīng)用在沖擊射流的實(shí)際應(yīng)用中。這是因?yàn)镃u納米顆粒的熱導(dǎo)率比Al納米顆粒以及 Al2O3納米顆粒的熱導(dǎo)率大,另外本實(shí)驗(yàn)使用的Al2O3納米顆粒粒徑較另外兩種納米顆粒粒徑大,對(duì)傳熱系數(shù)具有一定影響[21]。

    圖9 相同質(zhì)量分?jǐn)?shù)不同種類的納米流體與水傳熱系數(shù)隨Re數(shù)的變化

    4.2.3不同射流參數(shù)對(duì)傳熱系數(shù)的影響

    由于3種納米流體中質(zhì)量分?jǐn)?shù)為0.4%的Cu-水納米流體的換熱效果最好,因此選用Cu-水納米流體進(jìn)行射流距離對(duì)比實(shí)驗(yàn)。圖10給出在3種射流高度下0.4%質(zhì)量分?jǐn)?shù)的Cu-水納米流體傳熱系數(shù)隨Re的變化。從圖10中可以看出,射流高度存在一個(gè)特定的使傳熱系數(shù)達(dá)到最高的值。從實(shí)驗(yàn)結(jié)果發(fā)現(xiàn),在射流高度為3mm時(shí),傳熱系數(shù)較小。這是因?yàn)闆_擊射流的流場(chǎng)決定了換熱效果,而射流高度決定了流場(chǎng)的形態(tài)。當(dāng)射流高度很小的時(shí)候,沖擊在射流面的流體會(huì)有一部分反彈到射流噴嘴處,明顯降低正向射流的流速,使得圖10中射流高度3mm的曲線在高Re數(shù)狀態(tài)下傳熱系數(shù)提升不明顯。隨著射流高度的增加,射流得到充分的發(fā)展,紊流度增加,換熱效果增強(qiáng)。另外,當(dāng)射流高度繼續(xù)增大至9mm時(shí),由于射流高度過大,導(dǎo)致了在此過程中一定的能量損失,使得傳熱系數(shù)不升反降。

    4.3納米流體沖擊射流的流動(dòng)特性

    圖10 不同射流高度下Cu-水納米流體傳熱系數(shù)隨Re數(shù)的變化

    圖11展示了不同納米流體的射流壓力隨Re數(shù)的增大而增大,添加了納米顆粒的納米流體的壓降高于去離子水的壓降,這是由于濃度增加而導(dǎo)致黏度增大等原因,使得壓降也相應(yīng)的有所提高。圖11中,3種納米流體,Al2O3-水納米流體壓降最大,Al-水納米流體和Cu-水納米流體的壓降基本相同,表明添加了 Al2O3顆粒的納米黏度增大較為明顯導(dǎo)致流體流動(dòng)阻力大,因此壓降較大。

    圖12展示了Cu-水納米流體的噴射壓力隨Re數(shù)和濃度的變化。由圖12中可以看出,較去離子水的壓降相比,在相同雷諾數(shù)下,濃度越大,壓降越大。這是由于納米顆粒濃度的增大,使納米流體的黏度增大,管內(nèi)流動(dòng)阻力與噴嘴射流阻力都隨之增大,使壓降增大。

    圖11 不同納米流體的射流壓降隨Re數(shù)的變化

    圖12 不同質(zhì)量分?jǐn)?shù)下Cu-水納米流體的射流壓降隨Re數(shù)的變化

    5 結(jié) 論

    本實(shí)驗(yàn)中,通過利用去離子水、Cu-水、Al2O3-水、Al-水納米流體研究了納米流體沖擊射流的傳熱特性。使用不同種類的納米顆粒以及對(duì)應(yīng)其不同的質(zhì)量分?jǐn)?shù)進(jìn)行實(shí)驗(yàn)研究,Re數(shù)的范圍在 2000~8000,Nu數(shù)的范圍在10~70,經(jīng)實(shí)驗(yàn)測(cè)得數(shù)據(jù),得出以下結(jié)論。

    (1)3種納米流體中,Cu-水納米流體的射流冷卻效果要稍強(qiáng)于Al2O3-水納米流體,Al相對(duì)較差。但是添加了納米顆粒的納米流體的換熱效果都遠(yuǎn)好于去離子水的換熱效果。

    (2)在同一種納米顆粒制備的納米流體中,質(zhì)量分?jǐn)?shù)的增大使傳熱效果有明顯的提升。當(dāng)質(zhì)量分?jǐn)?shù)達(dá)到0.5%時(shí),由于黏度和懸浮性的影響,傳熱效果的變化明顯減小。以 Cu-水納米流體為例,質(zhì)量分?jǐn)?shù)分別為0.1%、0.2%、0.3%、0.4%、0.5%的Cu-水納米流體,傳熱系數(shù)分別提升 12.6%、34.6%、41%、51.2%、54.5%。

    (3)在Re<3000時(shí),納米流體相對(duì)于去離子水的換熱優(yōu)勢(shì)不明顯,所以在實(shí)際應(yīng)用中,應(yīng)保證有較高的流速。

    (4)實(shí)驗(yàn)中測(cè)量了不同射流高度對(duì)傳熱系數(shù)的影響。無論是納米流體還是水,都存在一個(gè)特定的高度,即本實(shí)驗(yàn)中高度為6mm時(shí),傳熱系數(shù)達(dá)到最大值。

    符號(hào)說明

    A—— 換熱面積,m2

    D —— 噴嘴直徑,mm

    H —— 射流距離,mm

    h——傳熱系數(shù),W/(m2·K)

    k——熱導(dǎo)率,W/(m·K)

    Nu—努塞爾數(shù)

    Q—輸入功率,W

    Ql—工質(zhì)獲得的功率,W

    q—熱量,J

    Re——雷諾數(shù)

    T——溫度,℃

    —平均溫度,℃

    ΔT —— 溫差,℃

    u——

    μ——?jiǎng)恿︷ざ?,Pa/s

    ρ—密度,kg/m3

    φ—體積分?jǐn)?shù),%

    下角標(biāo)

    bf—— 基液

    in—— 入口

    nf—— 納米流體

    out——出口

    p—納米顆粒

    w—— 壁面

    [1]WANG J C. L- and U- shaped heat pipes thermal modules with twin fans for cooling of electronic system under variable heat source areas[J]. Heat & Mass Transfer,2014,50(11):1487-1498.

    [2]SUFIAN S F,ABDULLAH M Z,MOHAMED J J. Effect of synchronized piezoelectric fans on microelectronic cooling performance [J]. International Communications in Heat & Mass Transfer,2013,43(4):81-89.

    [3]BRADBURY L J S. The structure of a self-preserving turbulent plane jet[J]. Journal of Fluid Mechanics,1965,23(1):31-64.

    [4]BAONGA J B,LOUAHLIA-GUALOUS H,IMBERT M. Experimental study of the hydrodynamic and heat transfer of free liquid jet impinging a flat circular heated disk[J]. Applied Thermal Engineering,2006,26(11/12):1125-1138.

    [5]KATTI V,PRABHU S V. Experimental study and theoretical analysis of local heat transfer distribution between smooth flat surface and impinging air jet from a circular straight pipe nozzle[J]. International Journal of Heat & Mass Transfer,2008,51:4480-4495.

    [6]NAPHON P,WONGWISES S. Investigation on the jet liquid impingement heat transfer for the central processing unit of personal computers[J]. International Communications in Heat & Mass Transfer,2010,37(7):822-826.

    [7]GUO D,WEI J J,ZHANG Y H. Enhanced flow boiling heat transfer with jet impingement on micro-pin-finned surfaces[J]. Applied Thermal Engineering,2011,31(11):2042-2051.

    [8]CHOI S U S. Enhancing thermal conductivity of fluids with nanoparticles[J]. ASME,1995,231:99-106.

    [9]WANG X,XU X,CHOI S U S. Thermal conductivity of nanoparticle -fluid mixture[J]. Journal of Thermophysics & Heat Transfer,1999,13 (13):474-80.

    [10]凌智勇,鄒濤,丁建寧,等. 納米流體黏度特性[J]. 化工學(xué)報(bào),2012,63(5):1409-1414.

    [11]李強(qiáng),宣益民. 納米流體熱導(dǎo)率的測(cè)量[J]. 化工學(xué)報(bào),2003,54 (1):42-46.

    [12]WUSIMAN K,JEONG H,TULUGAN K,et al. Thermal performance of multi-walled carbon nanotubes (MWCNTs) in aqueous suspensions with surfactants SDBS and SDS[J]. International Communications in Heat & Mass Transfer,2013,41(1):28-33.

    [13]HADDAD Z,ABID C,OZTOP H F,et al. A review on how the researchers prepare their nanofluids[J]. International Journal of Thermal Sciences,2014,76(76):168-189.

    [14]徐小嬌,劉妮,王玉強(qiáng),等. 納米流體懸浮液穩(wěn)定性的最新研究進(jìn)展[J]. 流體機(jī)械,2012,40(10):46-49.

    [15]毋偉,陳建峰,盧壽慈. 超細(xì)粉體表面修飾[M]. 北京:化學(xué)工業(yè)出版,2004:301-301.

    [16]SUNDAR L S,NAIK M T,SHARMA K V,et al. Experimental investigation of forced convection heat transfer and friction factor in a tube with Fe3O4magnetic nanofluid[J]. Experimental Thermal & Fluid Science,2012,37(2):65-71.

    [17]VENKATESHAN S P. Mechanical Measurements[M]. 2nd Edition. New York:John Wiley & Sons,Ltd.,2015:3-45.

    [18]TRCALA M,?ERMáK J. Nonlinear finite element analysis of thermal inertia in heat-balance sap flow measurement[J]. International Journal of Thermal Sciences,2014,76(2):200-207.

    [19]WUSIMAN K,JEONG H,TULUGAN K,et al. Thermal performance of multi-walled carbon nanotubes (MWCNTs) in aqueous suspensions with surfactants SDBS and SDS [J]. International Communications in Heat & Mass Transfer,2013,41 (1):28-33.

    [20]LI Q,XUAN Y,YU F. Experimental investigation of submerged single jet impingement using Cu-water nanofluid[J]. Applied Thermal Engineering,2012,36(2):426-433.

    [21]KELLY K L,CORONADO E,LIN L Z,et al. The optical properties of metal nanoparticles:the influence of size,shape,and dielectric environment[J]. Cheminform,2003,34(16):668-677.

    Solubility and metastable zone of valsartan in ethyl acetate

    WANG Qingqing,SUN Qin,YANG Asan,CHENG Rong,ZHENG Yanping
    (College of Chemical Engineering,Zhejiang University of Technology,Hangzhou 310014,Zhejiang,China)

    The solubility of valsartan in ethyl acetate was measured by a static equilibrium method at temperatures from 278.15K to 323.15K. The solubility data were correlated by the modified Apelblat equation,the NRTL equation,and the λh equation. In addition,the thermodynamic parameters (including the enthalpy,entropy and Gibbs energy) were calculated by the van't Hoff analysis. The results showed that the solubility of valsartan obviously increased with an increase of temperature. The experimental data were well correlated with the above models because the mean deviation values of these three models were 1.03%,3.87%,and 1.7%,respectively. Therefore,the modified Apelblat equation was the best choice in the correlation of the solubility data of valsartan. The values of thermodynamic parameters proved that the dissolution process of valsartan in ethyl acetate was endothermic,spontaneous and entropy-driven. The main contributor to the standard molar Gibbs energy of solution was the enthalpy change during the dissolution process. Furthermore,the supersolubility of valsartan in ethyl acetate was measured at different stirring rates and cooling rates by the laser method. The results indicated that the metastable zone was wider under the condition of low stirring rate,high temperature and large cooling rate. Meanwhile,the temperature effect was remarkable but the cooling rate had little effect on the metastable zone.

    valsartan;ethyl acetate;solubility;correlation;metastable zone

    O 795

    A

    1000-6613(2016)08-2329-05

    10.16085/j.issn.1000-6613.2016.08.04

    2016-01-11;修改稿日期:2016-01-27。

    王清清(1991—),男,碩士研究生,研究方向?yàn)榻Y(jié)晶動(dòng)力

    學(xué)。聯(lián)系人:孫勤,副教授,研究?jī)?nèi)容為結(jié)晶與干燥器設(shè)計(jì)。E-mail qins@zjut.edu.com。

    猜你喜歡
    傳熱
    換熱器強(qiáng)化傳熱方法及研究進(jìn)展
    世界家苑(2018年3期)2018-04-27 11:24:20
    雜談數(shù)學(xué)教學(xué)中的“熱處理”藝術(shù)
    考試周刊(2017年103期)2018-01-31 06:43:45
    干發(fā)酵過程中傳熱傳質(zhì)研究進(jìn)展
    綠色科技(2017年6期)2017-04-20 08:20:51
    不同外墻保溫體系框架柱熱橋二維非穩(wěn)態(tài)傳熱比較
    高爐熱風(fēng)爐理論與設(shè)計(jì)
    冬季供暖工況下混凝土輻射頂板傳熱特性
    管道化熔鹽段傳熱規(guī)律初探
    微通道換熱器的數(shù)值模擬和結(jié)構(gòu)優(yōu)化
    非牛頓流體繞鈍體傳熱的數(shù)值模擬
    高職化工原理課傳熱單元操作課堂教學(xué)改革初探
    亚洲自拍偷在线| 51午夜福利影视在线观看| 精品不卡国产一区二区三区| 精品少妇一区二区三区视频日本电影| 亚洲专区国产一区二区| 韩国av一区二区三区四区| 9191精品国产免费久久| 国内精品久久久久精免费| 夜夜躁狠狠躁天天躁| 国产日韩一区二区三区精品不卡| 9色porny在线观看| 亚洲国产精品久久男人天堂| 日本在线视频免费播放| 九色国产91popny在线| 国产成人欧美在线观看| 久久久精品欧美日韩精品| 久热这里只有精品99| 日韩 欧美 亚洲 中文字幕| 亚洲精品中文字幕在线视频| 亚洲va日本ⅴa欧美va伊人久久| 又大又爽又粗| 狠狠狠狠99中文字幕| 黄频高清免费视频| av视频在线观看入口| 精品国产乱码久久久久久男人| 黄频高清免费视频| 久久亚洲真实| 欧美在线一区亚洲| 欧美激情极品国产一区二区三区| 国产蜜桃级精品一区二区三区| 亚洲精品久久成人aⅴ小说| 一二三四在线观看免费中文在| 91精品国产国语对白视频| 欧美最黄视频在线播放免费| 欧美 亚洲 国产 日韩一| 久热爱精品视频在线9| 少妇 在线观看| 久久久久国产精品人妻aⅴ院| 久久性视频一级片| 搡老妇女老女人老熟妇| 欧美激情高清一区二区三区| 国产精品自产拍在线观看55亚洲| 俄罗斯特黄特色一大片| 嫩草影院精品99| 搡老妇女老女人老熟妇| 两个人视频免费观看高清| 日韩欧美免费精品| 一本综合久久免费| 怎么达到女性高潮| 人妻丰满熟妇av一区二区三区| 亚洲国产中文字幕在线视频| 老司机福利观看| 黑人欧美特级aaaaaa片| 国产人伦9x9x在线观看| 91大片在线观看| 最新美女视频免费是黄的| 欧美国产精品va在线观看不卡| 老汉色∧v一级毛片| 中文亚洲av片在线观看爽| 免费看a级黄色片| 亚洲精品美女久久久久99蜜臀| 欧美日韩中文字幕国产精品一区二区三区 | 国产午夜精品久久久久久| 一区二区日韩欧美中文字幕| 一进一出好大好爽视频| 美女高潮喷水抽搐中文字幕| 99国产综合亚洲精品| 高清毛片免费观看视频网站| 亚洲一卡2卡3卡4卡5卡精品中文| 黄片大片在线免费观看| 久久精品91无色码中文字幕| 制服诱惑二区| av有码第一页| 午夜福利,免费看| 免费av毛片视频| 91精品三级在线观看| avwww免费| 欧美日韩瑟瑟在线播放| 色av中文字幕| 久久久国产成人精品二区| 国产97色在线日韩免费| 性少妇av在线| 久久久久久久久久久久大奶| 免费av毛片视频| 日韩欧美国产一区二区入口| 婷婷精品国产亚洲av在线| 国产精品久久久久久亚洲av鲁大| 一区二区三区激情视频| 国产一级毛片七仙女欲春2 | 午夜久久久在线观看| 国产单亲对白刺激| 国内久久婷婷六月综合欲色啪| 国内久久婷婷六月综合欲色啪| 欧美久久黑人一区二区| 亚洲天堂国产精品一区在线| 人人澡人人妻人| 色播在线永久视频| 日日爽夜夜爽网站| 18禁观看日本| 91成年电影在线观看| 18美女黄网站色大片免费观看| 国产精品秋霞免费鲁丝片| 黄色a级毛片大全视频| 免费无遮挡裸体视频| 国产99白浆流出| 99热只有精品国产| 在线观看午夜福利视频| 我的亚洲天堂| 亚洲 欧美 日韩 在线 免费| 久久精品国产亚洲av香蕉五月| 日韩大码丰满熟妇| 午夜成年电影在线免费观看| 亚洲va日本ⅴa欧美va伊人久久| 97碰自拍视频| 亚洲精华国产精华精| 熟女少妇亚洲综合色aaa.| 久久久国产成人精品二区| 9191精品国产免费久久| 男女午夜视频在线观看| 精品国产乱码久久久久久男人| 日韩一卡2卡3卡4卡2021年| 日日摸夜夜添夜夜添小说| 亚洲av成人不卡在线观看播放网| 一二三四社区在线视频社区8| 十八禁人妻一区二区| www.熟女人妻精品国产| 国产1区2区3区精品| 国产精品免费视频内射| 黄片小视频在线播放| 日本免费a在线| 又大又爽又粗| 国产激情欧美一区二区| 婷婷精品国产亚洲av在线| 亚洲av日韩精品久久久久久密| 别揉我奶头~嗯~啊~动态视频| 啦啦啦韩国在线观看视频| 动漫黄色视频在线观看| 亚洲第一青青草原| 91成年电影在线观看| 日本vs欧美在线观看视频| 亚洲精品av麻豆狂野| 精品国产超薄肉色丝袜足j| 欧美激情极品国产一区二区三区| 在线观看免费午夜福利视频| 丁香六月欧美| 亚洲伊人色综图| 久久天堂一区二区三区四区| 欧美中文日本在线观看视频| 成年人黄色毛片网站| 一进一出抽搐动态| 国产精品久久电影中文字幕| 色综合婷婷激情| 亚洲欧美精品综合一区二区三区| 神马国产精品三级电影在线观看 | 国产一区二区在线av高清观看| 亚洲精品久久国产高清桃花| 欧美中文日本在线观看视频| 国产亚洲精品av在线| 老熟妇仑乱视频hdxx| 国产欧美日韩一区二区三| 日日干狠狠操夜夜爽| 99久久国产精品久久久| 久久精品影院6| 丝袜美足系列| 精品卡一卡二卡四卡免费| 国产男靠女视频免费网站| 色婷婷久久久亚洲欧美| 热re99久久国产66热| 在线观看舔阴道视频| 久久久久久免费高清国产稀缺| 老司机靠b影院| 少妇粗大呻吟视频| 美女国产高潮福利片在线看| 黑人巨大精品欧美一区二区蜜桃| 波多野结衣av一区二区av| 国产精品av久久久久免费| 色精品久久人妻99蜜桃| 国产精品秋霞免费鲁丝片| 国产欧美日韩综合在线一区二区| 欧美日韩亚洲综合一区二区三区_| 琪琪午夜伦伦电影理论片6080| 岛国视频午夜一区免费看| 女生性感内裤真人,穿戴方法视频| 两个人视频免费观看高清| 亚洲人成77777在线视频| 国产欧美日韩一区二区三| 一a级毛片在线观看| videosex国产| 欧美在线一区亚洲| 18美女黄网站色大片免费观看| 桃红色精品国产亚洲av| 黄色片一级片一级黄色片| 日韩视频一区二区在线观看| e午夜精品久久久久久久| 一本综合久久免费| 99国产极品粉嫩在线观看| 免费观看人在逋| 狂野欧美激情性xxxx| 又黄又爽又免费观看的视频| 久久天堂一区二区三区四区| 99香蕉大伊视频| 精品高清国产在线一区| 亚洲自拍偷在线| 国产精品久久久久久人妻精品电影| 国产激情久久老熟女| 美女国产高潮福利片在线看| 国产三级黄色录像| 成年女人毛片免费观看观看9| 免费一级毛片在线播放高清视频 | 精品久久久久久久毛片微露脸| 在线播放国产精品三级| 大码成人一级视频| 午夜免费激情av| 女人被狂操c到高潮| 可以在线观看毛片的网站| 又黄又粗又硬又大视频| 精品国产亚洲在线| 此物有八面人人有两片| 精品人妻在线不人妻| 成人国产一区最新在线观看| 国产乱人伦免费视频| 久久欧美精品欧美久久欧美| 日韩精品免费视频一区二区三区| 19禁男女啪啪无遮挡网站| av在线天堂中文字幕| x7x7x7水蜜桃| 国产av又大| 一级,二级,三级黄色视频| 欧美成狂野欧美在线观看| 午夜久久久在线观看| 91大片在线观看| 久久精品亚洲精品国产色婷小说| 亚洲最大成人中文| 9热在线视频观看99| 黄网站色视频无遮挡免费观看| 亚洲专区中文字幕在线| 免费少妇av软件| 视频区欧美日本亚洲| 香蕉国产在线看| 又黄又爽又免费观看的视频| 国产精品一区二区三区四区久久 | 色在线成人网| 国产精品免费一区二区三区在线| 国产精品,欧美在线| 伦理电影免费视频| 成人亚洲精品av一区二区| 久久久久久久久久久久大奶| 日本一区二区免费在线视频| 久久热在线av| 在线av久久热| 看片在线看免费视频| 精品不卡国产一区二区三区| 国产主播在线观看一区二区| 51午夜福利影视在线观看| 精品久久蜜臀av无| 日日夜夜操网爽| 午夜激情av网站| 人人妻人人澡欧美一区二区 | 久热这里只有精品99| 变态另类丝袜制服| 欧美在线一区亚洲| 91国产中文字幕| 欧美日韩福利视频一区二区| 人妻久久中文字幕网| 国产激情久久老熟女| 黄色毛片三级朝国网站| 亚洲av五月六月丁香网| 国产99白浆流出| 国产亚洲精品第一综合不卡| 99久久久亚洲精品蜜臀av| 国产xxxxx性猛交| 国产真人三级小视频在线观看| 乱人伦中国视频| 亚洲伊人色综图| 国产一区二区三区综合在线观看| 久久精品国产亚洲av高清一级| 日韩一卡2卡3卡4卡2021年| 中文字幕最新亚洲高清| 日韩精品中文字幕看吧| 亚洲精品一卡2卡三卡4卡5卡| 一区二区日韩欧美中文字幕| 日日干狠狠操夜夜爽| 极品人妻少妇av视频| 丰满人妻熟妇乱又伦精品不卡| 日韩视频一区二区在线观看| 满18在线观看网站| 国产麻豆69| 搡老熟女国产l中国老女人| 黄色 视频免费看| 在线永久观看黄色视频| 琪琪午夜伦伦电影理论片6080| 欧美黑人欧美精品刺激| 在线观看午夜福利视频| 久久草成人影院| 91精品国产国语对白视频| 亚洲精品在线美女| 宅男免费午夜| a级毛片在线看网站| 黄色丝袜av网址大全| 麻豆成人av在线观看| 精品乱码久久久久久99久播| 一本久久中文字幕| xxx96com| 色播在线永久视频| 女人被躁到高潮嗷嗷叫费观| 免费av毛片视频| 国产av又大| 亚洲精品国产色婷婷电影| 免费无遮挡裸体视频| 午夜影院日韩av| 久久久国产精品麻豆| 这个男人来自地球电影免费观看| 精品一区二区三区四区五区乱码| 国产亚洲精品一区二区www| 搡老岳熟女国产| 国产午夜精品久久久久久| 一级作爱视频免费观看| 两个人看的免费小视频| 看片在线看免费视频| 国产精品免费视频内射| 亚洲中文字幕日韩| 日韩 欧美 亚洲 中文字幕| 中国美女看黄片| 国产片内射在线| 岛国在线观看网站| 亚洲全国av大片| 91精品三级在线观看| 丝袜美腿诱惑在线| 久久久精品国产亚洲av高清涩受| 国产成人av激情在线播放| 亚洲专区字幕在线| 日韩成人在线观看一区二区三区| 亚洲aⅴ乱码一区二区在线播放 | 欧美日韩亚洲国产一区二区在线观看| 伦理电影免费视频| 精品不卡国产一区二区三区| 亚洲片人在线观看| 最近最新免费中文字幕在线| 国产激情久久老熟女| 一级a爱视频在线免费观看| 国产真人三级小视频在线观看| 丝袜美腿诱惑在线| 国产区一区二久久| 国产成+人综合+亚洲专区| 神马国产精品三级电影在线观看 | 亚洲精品国产色婷婷电影| 欧美不卡视频在线免费观看 | 中文亚洲av片在线观看爽| 日韩欧美三级三区| 日韩大码丰满熟妇| 色av中文字幕| 777久久人妻少妇嫩草av网站| 日韩欧美一区二区三区在线观看| 欧美日韩乱码在线| 美女扒开内裤让男人捅视频| 看片在线看免费视频| 久久香蕉精品热| 亚洲成人免费电影在线观看| 国产一卡二卡三卡精品| 桃红色精品国产亚洲av| 99久久精品国产亚洲精品| 亚洲中文日韩欧美视频| 亚洲成人久久性| 久久中文字幕人妻熟女| 韩国av一区二区三区四区| 12—13女人毛片做爰片一| 久久香蕉国产精品| 一级作爱视频免费观看| 亚洲av成人一区二区三| 免费看a级黄色片| 亚洲av熟女| 中文字幕最新亚洲高清| 妹子高潮喷水视频| 中文字幕另类日韩欧美亚洲嫩草| 美女 人体艺术 gogo| 亚洲国产日韩欧美精品在线观看 | 久久久国产成人精品二区| 国产在线精品亚洲第一网站| 亚洲最大成人中文| 婷婷六月久久综合丁香| 久久久久久久久久久久大奶| 在线播放国产精品三级| 侵犯人妻中文字幕一二三四区| 黄色女人牲交| 制服诱惑二区| 9色porny在线观看| 久久久久国产精品人妻aⅴ院| 亚洲国产毛片av蜜桃av| 人成视频在线观看免费观看| 午夜福利成人在线免费观看| 亚洲专区国产一区二区| 19禁男女啪啪无遮挡网站| 日本一区二区免费在线视频| 香蕉久久夜色| 午夜福利欧美成人| 欧美成狂野欧美在线观看| 久久青草综合色| 亚洲专区国产一区二区| 99国产精品免费福利视频| 午夜久久久在线观看| 国内久久婷婷六月综合欲色啪| 欧美久久黑人一区二区| 精品人妻1区二区| 国产av又大| 色av中文字幕| 亚洲少妇的诱惑av| 午夜亚洲福利在线播放| 久久伊人香网站| 曰老女人黄片| 国产91精品成人一区二区三区| 91国产中文字幕| 超碰成人久久| 欧美一区二区精品小视频在线| 国产一区二区三区综合在线观看| 午夜老司机福利片| 超碰成人久久| 波多野结衣av一区二区av| 国产精品香港三级国产av潘金莲| 日韩大码丰满熟妇| xxx96com| 午夜精品国产一区二区电影| 90打野战视频偷拍视频| 精品欧美一区二区三区在线| 欧美日韩精品网址| 亚洲少妇的诱惑av| 欧美黄色片欧美黄色片| 日韩av在线大香蕉| 久久人人爽av亚洲精品天堂| 99久久久亚洲精品蜜臀av| 亚洲av日韩精品久久久久久密| 黑人操中国人逼视频| 欧美激情极品国产一区二区三区| 午夜免费鲁丝| 国产一区二区在线av高清观看| 高潮久久久久久久久久久不卡| 91麻豆精品激情在线观看国产| 两个人看的免费小视频| 88av欧美| 黑人操中国人逼视频| 嫩草影院精品99| 不卡av一区二区三区| 免费女性裸体啪啪无遮挡网站| 亚洲国产日韩欧美精品在线观看 | 啦啦啦免费观看视频1| 极品人妻少妇av视频| 美女大奶头视频| 成人国语在线视频| 人人妻人人澡欧美一区二区 | 亚洲aⅴ乱码一区二区在线播放 | 性色av乱码一区二区三区2| 看免费av毛片| 欧美日韩亚洲综合一区二区三区_| 精品卡一卡二卡四卡免费| 精品免费久久久久久久清纯| 日韩精品中文字幕看吧| 免费搜索国产男女视频| 国产99白浆流出| 两性夫妻黄色片| 91九色精品人成在线观看| 免费人成视频x8x8入口观看| 亚洲色图av天堂| 两性夫妻黄色片| 亚洲伊人色综图| 一卡2卡三卡四卡精品乱码亚洲| 一进一出好大好爽视频| 麻豆一二三区av精品| 久久人妻福利社区极品人妻图片| 久久国产精品影院| 女警被强在线播放| 国产又色又爽无遮挡免费看| 欧美一级a爱片免费观看看 | 日本精品一区二区三区蜜桃| 一级毛片高清免费大全| 久久久久久久精品吃奶| 久久婷婷成人综合色麻豆| av片东京热男人的天堂| 男人舔女人下体高潮全视频| 中文字幕久久专区| 男女之事视频高清在线观看| 亚洲熟女毛片儿| 国产成人一区二区三区免费视频网站| 男人舔女人下体高潮全视频| 午夜免费鲁丝| 变态另类丝袜制服| 久久精品亚洲熟妇少妇任你| 最近最新中文字幕大全免费视频| 欧美日韩黄片免| 午夜精品久久久久久毛片777| 精品不卡国产一区二区三区| 亚洲狠狠婷婷综合久久图片| 久久人妻av系列| 脱女人内裤的视频| 麻豆一二三区av精品| 久久婷婷人人爽人人干人人爱 | 国产精品亚洲av一区麻豆| 手机成人av网站| 欧美老熟妇乱子伦牲交| 国产亚洲欧美98| 亚洲三区欧美一区| 一二三四社区在线视频社区8| www国产在线视频色| 十八禁人妻一区二区| 老司机深夜福利视频在线观看| 一区二区日韩欧美中文字幕| 国产av又大| 亚洲色图av天堂| 一个人免费在线观看的高清视频| 午夜久久久久精精品| 国产成人系列免费观看| 久久久久久亚洲精品国产蜜桃av| 999精品在线视频| 黄色视频,在线免费观看| 日韩视频一区二区在线观看| 国产熟女午夜一区二区三区| 男男h啪啪无遮挡| 国产91精品成人一区二区三区| 性少妇av在线| 黑人巨大精品欧美一区二区mp4| 国产亚洲精品久久久久5区| 黄片播放在线免费| 人人妻人人澡欧美一区二区 | 狂野欧美激情性xxxx| 999久久久精品免费观看国产| 18禁美女被吸乳视频| 国产成人精品久久二区二区免费| 欧美黄色片欧美黄色片| 国产精品亚洲av一区麻豆| 天天躁狠狠躁夜夜躁狠狠躁| 在线观看日韩欧美| 亚洲专区国产一区二区| 黄频高清免费视频| 看黄色毛片网站| 国产精品日韩av在线免费观看 | 妹子高潮喷水视频| 久久中文字幕一级| 欧美精品啪啪一区二区三区| 亚洲国产欧美一区二区综合| 欧美精品亚洲一区二区| 一个人观看的视频www高清免费观看 | 久9热在线精品视频| 久久婷婷成人综合色麻豆| 久99久视频精品免费| 欧美黑人精品巨大| 欧美久久黑人一区二区| 淫秽高清视频在线观看| 亚洲午夜精品一区,二区,三区| 日韩高清综合在线| 亚洲美女黄片视频| 好男人在线观看高清免费视频 | 色婷婷久久久亚洲欧美| av视频在线观看入口| 国产精品综合久久久久久久免费 | 久久青草综合色| 51午夜福利影视在线观看| 757午夜福利合集在线观看| 老鸭窝网址在线观看| 国产精品亚洲一级av第二区| 色精品久久人妻99蜜桃| 亚洲五月色婷婷综合| 人成视频在线观看免费观看| 曰老女人黄片| 亚洲aⅴ乱码一区二区在线播放 | 亚洲在线自拍视频| 一二三四社区在线视频社区8| 精品高清国产在线一区| 亚洲人成电影观看| 欧美一区二区精品小视频在线| 日韩欧美在线二视频| 亚洲男人天堂网一区| 亚洲精品av麻豆狂野| 午夜福利在线观看吧| 999久久久精品免费观看国产| 老司机午夜福利在线观看视频| 亚洲激情在线av| av天堂久久9| 免费观看人在逋| 一级黄色大片毛片| 成人国语在线视频| 精品国产美女av久久久久小说| 国产伦一二天堂av在线观看| 一级片免费观看大全| 国产精品久久久久久亚洲av鲁大| 国产激情久久老熟女| 国产91精品成人一区二区三区| 两性午夜刺激爽爽歪歪视频在线观看 | 黄色毛片三级朝国网站| 国内毛片毛片毛片毛片毛片| 每晚都被弄得嗷嗷叫到高潮| 国内精品久久久久久久电影| 国产精品免费一区二区三区在线| 亚洲中文日韩欧美视频| 亚洲黑人精品在线| 老司机在亚洲福利影院| 午夜福利在线观看吧| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲天堂国产精品一区在线| 久久精品国产清高在天天线| 制服诱惑二区| 91字幕亚洲| 国产高清videossex| 久久亚洲精品不卡| 久久伊人香网站| 国产亚洲精品久久久久5区| 国产精品国产高清国产av| 久久天躁狠狠躁夜夜2o2o| 伊人久久大香线蕉亚洲五| 午夜精品在线福利| 欧美乱码精品一区二区三区| 久久精品91无色码中文字幕| 国产午夜福利久久久久久| 亚洲成a人片在线一区二区| 精品国产美女av久久久久小说| 波多野结衣巨乳人妻| 欧美激情极品国产一区二区三区| 激情在线观看视频在线高清|