高溫金屬熔體黏度突變探索?

商繼祥 趙云波 胡麗娜

(山東大學(xué),材料液固結(jié)構(gòu)演變與加工教育部重點(diǎn)實(shí)驗(yàn)室,濟(jì)南 250061)

1 引 言

高溫金屬熔體黏度作為描述液體動(dòng)力學(xué)性質(zhì)的物理參量,體現(xiàn)了液體結(jié)構(gòu)的演變過(guò)程,始終受到廣泛的關(guān)注.在傳統(tǒng)鑄造工藝中,高溫金屬熔體黏度變化是決定晶態(tài)鑄件成型能力和氣泡或非金屬夾雜物能夠順利排出的關(guān)鍵因素[1].對(duì)于利用超快速冷卻技術(shù)制備的金屬玻璃而言,高溫金屬熔體黏度大小更是決定該合金能否形成金屬玻璃或者玻璃形態(tài)的要素之一.

1985年,美國(guó)科學(xué)家Angell[2]首先用液體的脆性概念表征了黏度隨溫度變化的難易程度.當(dāng)過(guò)冷液體的黏度隨溫度的變化較大時(shí),即液體的結(jié)構(gòu)隨溫度的變化容易發(fā)生變化時(shí),液體是脆性的,在Angell圖上表現(xiàn)為一條曲線;當(dāng)過(guò)冷液體的黏度隨溫度的變化很難發(fā)生改變時(shí),即液體結(jié)構(gòu)不容易發(fā)生變化,為強(qiáng)性液體,在Angell圖上表現(xiàn)為一條直線.借用該過(guò)冷液體的脆性定義,Bian等[3]采用液相線約化溫度來(lái)表示高溫金屬熔體黏度在液相線附近變化的難易程度,提出了過(guò)熱熔體的脆性概念,并發(fā)現(xiàn)該脆性大小與熔體中團(tuán)簇的穩(wěn)定性密切相關(guān).黏度在衡量金屬玻璃形成能力方面也有一定的參考性,玻璃形成能力強(qiáng)的合金體系往往有著較大的黏度值.這是由于大的黏度值可以抑制原子或分子定向移動(dòng)從而抑制晶化過(guò)程,降低了結(jié)晶的可能性,從而有利于金屬玻璃的形成.例如,在共晶凝固的過(guò)程中,當(dāng)溫度降低到共晶點(diǎn)溫度時(shí),液相中會(huì)同時(shí)析出兩種不同成分的固溶體相,直至凝固結(jié)束.兩種不同成分固溶體的形核,以及形核后的晶體長(zhǎng)大過(guò)程都與金屬原子的移動(dòng)速率或擴(kuò)散行為密切相關(guān).當(dāng)金屬原子的遷移因較大的黏滯力受阻后,熔體能在較低的溫度下仍保持原子排列的無(wú)序狀態(tài),從而有利于非晶的形成.以往的研究工作已經(jīng)表明過(guò)熱熔體的脆性系數(shù)與金屬玻璃形成能力的大小存在著負(fù)相關(guān)性[3,4],黏度的變化特征對(duì)合金熔體玻璃形成能力的影響舉足輕重.金屬玻璃形成能力的大小對(duì)拓寬金屬玻璃的應(yīng)用范圍,實(shí)現(xiàn)金屬玻璃從功能材料向結(jié)構(gòu)材料質(zhì)的跨越起著關(guān)鍵性的作用.通常認(rèn)為,高溫金屬熔體的黏度在熔點(diǎn)以上隨著溫度的變化規(guī)律符合Arrhenius公式.但越來(lái)越多的實(shí)驗(yàn)結(jié)果表明,高溫金屬熔體的動(dòng)力學(xué)性質(zhì)存在轉(zhuǎn)折.高溫金屬熔體的黏度在很多情況下不能用同一條Arrhenius曲線進(jìn)行很好的擬合:在特定的溫度處黏度偏離Arrhenius曲線,并在偏離的基礎(chǔ)上繼續(xù)符合另外一條新的Arrhenius曲線.這種黏度突變的行為是液-液相變的直觀體現(xiàn),液-液相變現(xiàn)象的發(fā)生勢(shì)必伴隨著熔體中金屬原子的重組,即熔體中原子或原子團(tuán)簇從一種狀態(tài)轉(zhuǎn)變?yōu)榱硪环N狀態(tài).與固態(tài)相變類似,液-液相變也伴隨著體積、熵或焓的變化.加深液-液相變深層次的研究,對(duì)深入認(rèn)識(shí)液態(tài)金屬的性質(zhì)著重要的理論意義.而對(duì)于該現(xiàn)象的揭示,其動(dòng)力學(xué)性質(zhì)即黏度隨溫度的變化可以給我們提供明顯的論證依據(jù).同時(shí),高溫金屬熔體黏度的突變也極可能是金屬玻璃固體異常弛豫行為的根源.

因此,對(duì)高溫金屬熔體黏度突變的探索研究,不僅有助于認(rèn)識(shí)高溫金屬熔體的動(dòng)力學(xué)性質(zhì)的演變規(guī)律,而且對(duì)它的探索有助于進(jìn)一步認(rèn)識(shí)合金體系的金屬玻璃形成能力和液-液相變規(guī)律,并將高溫合金熔體的動(dòng)力學(xué)性質(zhì)與合金體系的勢(shì)能圖譜、固體的不同弛豫行為相聯(lián)系[3?6].這不僅對(duì)實(shí)際工業(yè)生產(chǎn)具有指導(dǎo)意義,同時(shí)對(duì)理解熔體過(guò)-冷液體-金屬固體的遺傳特征也具有重要的理論價(jià)值.這也是今后金屬熔體動(dòng)力學(xué)性質(zhì)研究的重要發(fā)展方向之一.

2 高溫金屬熔體的黏度測(cè)量原理

黏度的物理意義定義為:在單位速度梯度下,作用在單位接觸面積上的內(nèi)摩擦力.

式中η表示為動(dòng)力學(xué)黏度,F為內(nèi)摩擦力,A為內(nèi)摩擦力與兩液體層接觸面的面積,dv/dy為速度梯度.如圖1所示,設(shè)距離為dy的相鄰液體層以不同的速度v和v+dv的速度流動(dòng),即垂直于剪切力方向產(chǎn)生速度梯度,則在這兩液層之間,將產(chǎn)生一對(duì)大小相等而方向相反的內(nèi)摩擦力,力圖阻止兩液層發(fā)生互相運(yùn)動(dòng),產(chǎn)生的這一阻力就是液體的黏滯力.

圖1 黏度定義示意圖Fig.1.Diagram of viscosity definition.

雖然用于測(cè)量黏度的方法有很多(如毛細(xì)血管法、落體法、旋轉(zhuǎn)法、聲波法等),但適合高溫金屬熔體黏度測(cè)量的方法并不多.這往往是由于高溫金屬熔體的黏度比較小,即使含有稀土元素等大分子量的Fe基、Co基等體系,黏度也不過(guò)只有十幾mPa·s左右,而且高溫金屬熔體黏度的測(cè)量是在高溫下進(jìn)行的,極易被氧化,因此一些常規(guī)的測(cè)量黏度的方法并不適合高溫金屬熔體黏度的測(cè)量.目前用于測(cè)量高溫金屬熔體黏度的方法主要是震蕩法.基于震蕩法的黏度儀有兩種形式,一種是懸掛式[7?12],一種是站立式[13].相比于站立式,懸掛式更容易操作和制作,所以實(shí)驗(yàn)室中測(cè)量高溫金屬熔體的黏度儀主要是懸掛式.懸掛式黏度儀的結(jié)構(gòu)如圖2(a)所示,金屬熔體置于圓柱形坩堝內(nèi),坩堝由懸絲懸掛,懸絲與步進(jìn)電機(jī)相連,單片機(jī)控制步進(jìn)電機(jī)的驅(qū)動(dòng)器,使步進(jìn)電機(jī)轉(zhuǎn)動(dòng),從而使懸絲帶動(dòng)坩堝沿垂直軸旋轉(zhuǎn),在懸絲恢復(fù)力的作用下,坩堝做震蕩運(yùn)動(dòng),由于液體黏滯的作用,震蕩運(yùn)動(dòng)逐漸衰減,但運(yùn)動(dòng)周期保持不變,在懸絲放置一個(gè)反光鏡,用于反射激光,計(jì)算運(yùn)動(dòng)周期和對(duì)數(shù)衰減率.測(cè)量系統(tǒng)如圖2(b)所示,A,B兩點(diǎn)為光敏光電池用于接受反射激光,在開(kāi)始測(cè)量之前需要將反射光電置于A點(diǎn),步進(jìn)電機(jī)轉(zhuǎn)動(dòng)之后,反射激光在A和B兩光電池之間來(lái)回?cái)[動(dòng),光電池通關(guān)光電轉(zhuǎn)換效應(yīng),計(jì)算出激光掃過(guò)AA,AB,BA,BB的時(shí)間,這樣便可以算出周期S和對(duì)數(shù)衰減率σ,再通過(guò)Shvidkovskiy公式[14]計(jì)算出黏度值.

圖2 (a)黏度儀示意圖[15];(b)測(cè)量系統(tǒng)示意圖[15]Fig.2.(a)Schematic of viscometer[15];(b)schematic of measuring system[15].

震蕩運(yùn)動(dòng)的方程為:

式中θ0是前指因子,S為震蕩周期,σ為對(duì)數(shù)衰減率,t為時(shí)間.

Shvidkovskiy公式計(jì)算黏度:

式中I表示轉(zhuǎn)動(dòng)慣量,m表示測(cè)試樣品的質(zhì)量,R表示為坩堝的內(nèi)徑,S0表示空坩堝震蕩運(yùn)動(dòng)周期,σ0表示空坩堝的對(duì)衰減率,W是矯正系數(shù),υ表示為運(yùn)動(dòng)學(xué)黏度,ρ為測(cè)試樣品的密度.

不同的實(shí)驗(yàn)條件和試驗(yàn)方法對(duì)黏度測(cè)量結(jié)果影響較大.以前工作已經(jīng)表明磁場(chǎng)強(qiáng)度、坩堝表面粗糙度、保溫時(shí)間以及放料的多少對(duì)黏度有著重要的影響.Xu等[15]報(bào)道了Ga80Cr20和Ga80Ni20在不同磁場(chǎng)強(qiáng)度下黏度的變化規(guī)律,實(shí)驗(yàn)結(jié)果表明:隨著磁場(chǎng)強(qiáng)度不斷增加,熔體的黏度不斷增加.這是由于高溫金屬熔體的團(tuán)簇帶有一定量的電荷,帶電離子在磁場(chǎng)中會(huì)受到洛倫茲力,洛倫茲力改變團(tuán)簇的運(yùn)動(dòng)路徑,從而使團(tuán)簇的運(yùn)動(dòng)更為復(fù)雜,高溫金屬熔體的黏滯力增加.磁場(chǎng)強(qiáng)度H越大,黏度值也越大.Wu等[16]報(bào)道了測(cè)量黏度時(shí)所用坩堝的表面粗糙程度對(duì)Pb38.1Sn61.9高溫金屬熔體黏度的影響,實(shí)驗(yàn)結(jié)果表明:坩堝的表面粗糙程度對(duì)黏度有一定的影響,且表面粗糙度越大,在同一溫度下測(cè)得的黏度值越大.這是因?yàn)槟Σ亮Φ脑黾訉?dǎo)致對(duì)數(shù)衰減率δ的增加,黏度也相應(yīng)地增大.Sun等[17]的工作也證實(shí)了坩堝表面的粗糙程度對(duì)水的黏度具有正相關(guān)的關(guān)系:表面粗糙程度越大,水的黏度也越大.放料的多少對(duì)黏度影響也比較大[18]:當(dāng)樣品體積較小時(shí),測(cè)量的黏度數(shù)據(jù)大于標(biāo)準(zhǔn)數(shù)值,隨著體積的增大,黏度數(shù)值減小,逐漸接近標(biāo)準(zhǔn)值;當(dāng)體積繼續(xù)增大時(shí),黏度數(shù)值進(jìn)一步降低,低于標(biāo)準(zhǔn)值.這是因?yàn)闃悠愤^(guò)多,在振蕩運(yùn)動(dòng)過(guò)程中會(huì)溢出,造成黏度測(cè)量結(jié)果的誤差.在進(jìn)行高溫金屬熔體黏度測(cè)量時(shí),金屬熔體體積會(huì)發(fā)生膨脹,所以體積也因此不宜過(guò)大.保溫時(shí)間的長(zhǎng)短同樣對(duì)黏度有較大的影響,保溫初期,合金熔體黏度隨保溫時(shí)間變化會(huì)有所波動(dòng),但隨著保溫時(shí)間延長(zhǎng),合金熔體黏度波動(dòng)逐漸減小.為了得到盡可能準(zhǔn)確的黏度值,保溫時(shí)間應(yīng)控制在30 min以上[19].造成黏度隨保溫時(shí)間變化的原因可以從兩個(gè)角度來(lái)考慮:首先在保溫初期控制溫度和爐內(nèi)的溫度是不同步的,隨著保溫時(shí)間的增加,由于熱輻射的影響,內(nèi)層熔體的溫度逐漸與坩堝套的溫度趨于一致,黏度值趨于穩(wěn)定.其次在保溫初期,熔體經(jīng)過(guò)前期溫度的變化,原子運(yùn)動(dòng)相對(duì)比較劇烈,合金熔體的成分相對(duì)而言不均勻,合金熔體需要一定的時(shí)間達(dá)到均勻態(tài),這就導(dǎo)致測(cè)量結(jié)果隨保溫時(shí)間略有波動(dòng).

3 高溫金屬熔體的黏度

3.1 純組元單質(zhì)的黏度

目前,對(duì)金屬單質(zhì)高溫金屬熔體的黏度研究主要集中在較低熔點(diǎn)、不易氧化的元素種類上,對(duì)稀土基易氧化金屬以及Zr等高熔點(diǎn)金屬的黏度研究開(kāi)展得很少.這主要是受到實(shí)驗(yàn)條件(如真空度、最大溫度測(cè)量范圍)的限制.一般認(rèn)為,高溫金屬熔體的黏度在熔點(diǎn)以上隨著溫度的降低,黏度值不斷增大,并符合Arrhenius公式,其中η0是前指因子,E是流體團(tuán)簇的激活能,k為玻爾茲曼常數(shù),T為溫度.但相當(dāng)一部分實(shí)驗(yàn)結(jié)果表明在整個(gè)溫度范圍內(nèi),高溫金屬熔體的黏度隨溫度的變化并不是單調(diào)的,在某些特定的溫度附近,黏度存在突變現(xiàn)象,并在新變化的基礎(chǔ)上繼續(xù)符合Arrhenius公式.越來(lái)越多的研究結(jié)果顯示:高溫金屬熔體的黏度突變是一普遍的現(xiàn)象,而這一現(xiàn)象與熔體中團(tuán)簇的變化密切相關(guān).孫民華等[22]報(bào)道了Al熔體的黏度,研究了Al熔體黏度隨溫度的變化規(guī)律,發(fā)現(xiàn)在升溫過(guò)程中熔體黏度值在780?C和950?C左右發(fā)生突變;在降溫過(guò)程中,黏度的突變發(fā)生在930?C與750?C,如圖3所示.孫民華等通過(guò)Al熔體氫含量的測(cè)定以及分子動(dòng)力學(xué)模擬方法,進(jìn)一步探索了Al熔體液態(tài)微觀結(jié)構(gòu)與熔體黏度的內(nèi)在聯(lián)系.含氫量的測(cè)定結(jié)果表明:氫含量隨溫度變化在780?C發(fā)生突變.通過(guò)對(duì)液態(tài)Al的分子動(dòng)力學(xué)模擬,發(fā)現(xiàn)Al的第一近鄰配位原子的排布方式隨溫度的變化在780?C左右與950?C左右也存在突變.這似乎可以說(shuō)明含氫量以及配位數(shù)的變化是導(dǎo)致黏度突變的原因.至于在升溫和降溫過(guò)程中發(fā)生黏度突變的溫度不同,降溫過(guò)程黏度突變點(diǎn)存在的滯后性,孫民華等認(rèn)為是由于熔煉后的母合金存在微觀不均勻性,升溫過(guò)程消除了母合金的熱歷史所致.

為了更加具體地體現(xiàn)黏度在不同溫度段的Arrhenius關(guān)系,將黏度取對(duì)數(shù)后,并對(duì)溫度進(jìn)行適當(dāng)?shù)淖冃?得到圖4.可以看出,其黏度隨溫度的變化可以用三個(gè)Arrhenius公式來(lái)擬合.

圖3 (a)加熱時(shí)Al熔體的黏度隨溫度的變化;(b)降溫時(shí)Al熔體的黏度隨溫度的變化[22]Fig.3.(a)Variation cure of Al melt viscosity with temperature during heating;(b)variation cure of Al melt viscosity with temperature during cooling[22].

圖4 (a)升溫時(shí)Al熔體的Arrhenius曲線logη-1/(10?3T);(b)降溫時(shí)Al熔體的Arrhenius曲線logη-1/(10?3T)Fig.4.(a)Arrhenius plot logη-1/(10?3T)of Al viscosity during heating;(b)arrhenius plot logη-1/(10?3T)of Al viscosity during cooling.

Wang等[23]報(bào)道了Sb,In,Sn,Bi高溫金屬熔體的黏度,并通過(guò)差示掃描量熱法(DSC)掃描進(jìn)行了分析.結(jié)果顯示在黏度突變點(diǎn)附近,DSC曲線有明顯的放熱峰或吸熱峰的存在,表明在黏度突變點(diǎn)附近高溫金屬熔體發(fā)生了熱效應(yīng),此熱效應(yīng)是液態(tài)相變所致.Arrhenius曲線不連續(xù)點(diǎn)被認(rèn)為是相變點(diǎn),相變點(diǎn)把Arrhenius曲線分成幾個(gè)不同的溫度區(qū)域.在每一個(gè)溫度區(qū)域范圍內(nèi),公式中與黏度計(jì)算有關(guān)的系數(shù)η0和E不依懶于溫度,只依賴于在某一特定的溫度范圍內(nèi)液體的結(jié)構(gòu)和彌散分布的流體團(tuán)簇的尺寸.也就是說(shuō)在黏度突變點(diǎn)前后的η0和E是不一樣的,在沒(méi)發(fā)生黏度突變之前這一溫度范圍內(nèi)η0和E是一樣的,不隨溫度的變化而變化.圖5為Sb,In,Sn,Bi高溫金屬熔體的Arrhenius曲線,可以看出Sb,In的黏度突變現(xiàn)象并沒(méi)有Sn,Bi那么明顯,并且Sb,Sn,Bi只有一個(gè)黏度突變點(diǎn),分別在670?C,400?C,480?C發(fā)生突變,而In存在兩個(gè)黏度突變點(diǎn),在620?C,300?C發(fā)生突變.

圖5 (a)Sb,(b)In,(c)Sn,(d)Bi Arrhenius曲線logη-1/(10?3T)[23]Fig.5.(a)Sb,(b)In,(c)Sn,(d)Bi Arrhenius plot logη-1/(10?3T)[23].

有關(guān)Pb[24?26],Ni[27?29],Ag[27,30,31]黏度,不同文獻(xiàn)中報(bào)道的黏度值相差比較大.這極有可能是由于黏度的測(cè)量比較復(fù)雜、高溫下熔體極易被氧化、不可控雜質(zhì)濃度的影響、數(shù)據(jù)處理的方法或者實(shí)驗(yàn)條件的不同導(dǎo)致的.

3.2 二元合金體系的黏度

由于存在不同種類原子之間的相互作用,合金熔體的結(jié)構(gòu)往往要比金屬單質(zhì)熔體結(jié)構(gòu)復(fù)雜得多.合金在凝固的過(guò)程中除了存在能量起伏和結(jié)構(gòu)起伏外,還存在著成分起伏.可以認(rèn)為,合金有著比金屬單質(zhì)更為復(fù)雜的物化性質(zhì)以及更廣泛的應(yīng)用領(lǐng)域.目前文獻(xiàn)報(bào)道的二元合金體系的黏度值相對(duì)來(lái)說(shuō)比金屬單質(zhì)要多.Zhao等[32]報(bào)道了不同組分下的CuZr二元合金體系,如圖6所示.實(shí)驗(yàn)結(jié)果表明:不同組分下CuZr合金在1400 K左右均存在明顯的黏度突變,說(shuō)明液-液相變的發(fā)生;隨著Cu原子分?jǐn)?shù)的不斷增加,黏度突變發(fā)生時(shí)對(duì)應(yīng)的突變黏度值是不相同的,Cu60Zr40黏度突變時(shí)對(duì)應(yīng)的突變黏度值最小;不同組分下CuZr合金發(fā)生黏度突變的溫度范圍不同,且Cu62Zr38黏度突變最緩慢(黏度突變?cè)谝粋€(gè)很寬的溫度范圍內(nèi)完成).高溫金屬熔體在黏度突變范圍以上被定義為高溫轉(zhuǎn)變區(qū)(HTZ),黏度突變范圍以下被定義為低溫轉(zhuǎn)變區(qū)(LTZ).為了說(shuō)明CuZr熔體的動(dòng)力學(xué)變化,作者引用了過(guò)熱熔體脆性的概念M,分別計(jì)算了HTZ和LTZ的M值,將MHTZ與MLTZ的比值定義為轉(zhuǎn)變強(qiáng)度系數(shù)F.計(jì)算結(jié)果表明:MHTZ的值總是小于MLTZ的值,且F的值越大,玻璃形成能力越強(qiáng).從這一結(jié)果可以看出,高溫金屬熔體黏度的變化情況極有可能成為衡量金屬玻璃形成能力的一個(gè)重要判據(jù).而這個(gè)工作也直接體現(xiàn)了金屬玻璃制備過(guò)程中存在的液-固相關(guān)性.

圖6 (a)Cu48Zr52,(b)Cu50Zr50,(c)Cu52Zr48,(d)Cu54Zr46,(e)Cu56Zr44,(f)Cu58Zr42,(g)Cu60Zr40,(h)Cu62Zr38合金Arrhenius曲線logη-1/(10?3T)[32]Fig.6.(a)Cu48Zr52,(b)Cu50Zr50,(c)Cu52Zr48,(d)Cu54Zr46,(e)Cu56Zr44,(f)Cu58Zr42,(g)Cu60Zr40,(h)Cu62Zr38melt Arrhenius polt logη-1/(10?3T)[32]

Zhao等[33]報(bào)道了Cu60Sn40黏度隨溫度變化的規(guī)律,并結(jié)合Cu-Sn合金相圖分析了黏度非連續(xù)變化的原因.通過(guò)CuSn合金相圖可以發(fā)現(xiàn),化學(xué)計(jì)量數(shù)為Cu3Sn的成分在固態(tài)Cu-Sn合金中有一個(gè)很寬的存在范圍,在CuSn高溫金屬熔體中,存在著Cu3Sn團(tuán)簇結(jié)構(gòu).Cu60Sn40中的液-液相變,即黏度突變與Sn原子的聚集成鍵以及Cu3Sn團(tuán)簇結(jié)構(gòu)的形成有關(guān).對(duì)于大多數(shù)物質(zhì),蒸發(fā)焓要比熔化焓大數(shù)倍,這就意味著Cu60Sn40熔體在剛開(kāi)始熔化時(shí)只有部分的原子鍵破壞,剩余的大量原子以Cu3Sn團(tuán)簇和Sn團(tuán)簇的形式存在于高溫金屬熔體中.當(dāng)溫度超過(guò)900?C時(shí),Cu3Sn團(tuán)簇中的Cu—Sn鍵以及Sn團(tuán)簇中的Sn—Sn鍵被破壞,熔體傾向于均勻態(tài),黏度也因此發(fā)生突變.同理在降溫過(guò)程中,有相反的過(guò)程發(fā)生,Cu3Sn團(tuán)簇和Sn團(tuán)簇在900?C附近生成,并隨著溫度的逐漸降低結(jié)合更加牢固,黏度再次出現(xiàn)突變.這一可逆過(guò)程可以表示為:L(Cu3Sn)+L(Sn)?L(Cu60Sn40).更多的研究結(jié)果證明了相圖與高溫金屬熔體黏度突變可能存在的聯(lián)系:有多種中間化合物生成的相圖,其高溫金屬熔體的黏度突變程度往往比較明顯[32,34];流體團(tuán)簇的激活能的大小與液相線溫度存在正相關(guān)性[35];在熔點(diǎn)或液相線溫度以上相同的過(guò)熱度條件下,有金屬間化合物存在的區(qū)域黏度值較大[36].由于三元及以上相圖比較復(fù)雜,目前高溫合金熔體和相圖的關(guān)系研究主要集中在二元合金.

圖7 Cu60Sn40合金Arrhenius曲線logη-1/(10?3T)[33]Fig.7.Cu60Sn40melt Arrhenius plot logη-1/(10?3T)[33].

圖8[15]是Ga80Cr20在不同磁場(chǎng)強(qiáng)度下黏度的Arrhenius曲線.在磁場(chǎng)強(qiáng)度分別為0,650,935 Gs時(shí),Ga80Cr20均存在黏度突變,且黏度突變的溫度范圍一致.由圖8可以看出,磁場(chǎng)強(qiáng)度在使黏度升高的同時(shí)不會(huì)使高溫金屬熔體黏度的突變行為消失.這說(shuō)明黏度突變是金屬玻璃高溫金屬熔體本質(zhì)的體現(xiàn),不隨外界條件(比如磁場(chǎng)強(qiáng)度、升降溫速過(guò)程)的改變而改變.Xu等[15]進(jìn)一步定量描述了外加磁場(chǎng)與黏度之間的關(guān)系:黏度與磁場(chǎng)強(qiáng)度的平方成正比.二者之間的比例系數(shù)與電阻率有關(guān),且隨著熔體的溫度的升高,電阻率逐漸增大,這樣就可以通過(guò)對(duì)電阻率的測(cè)量來(lái)預(yù)測(cè)黏度在磁場(chǎng)作用下的變化情況.

圖8 不同磁場(chǎng)強(qiáng)度下的Ga80Cr20Arrhenius曲線logη-1/(10?3T)[15]Fig.8.Ga80Cr20melt Arrhenius plot logη-1/(10?3T)in different magnetic field intensity[15].

Konstantinova等[37]在Cu-Al合金熔體的運(yùn)動(dòng)學(xué)黏度一文中,探究了Cu-Al合金熔體的運(yùn)動(dòng)學(xué)黏度隨溫度和成分的變化規(guī)律.圖9所示為1100?C CuAl合金隨著Al含量的增加運(yùn)動(dòng)學(xué)黏度的變化情況.從圖9可以看出在Cu70Al30附近運(yùn)動(dòng)學(xué)黏度有一個(gè)平緩的極值,在化學(xué)計(jì)量數(shù)為CuAl3附近運(yùn)動(dòng)學(xué)黏度有一個(gè)尖銳的極值.根據(jù)Arrhenius公式計(jì)算出的不同組分下的η0和E發(fā)現(xiàn),在Cu70Al30附近有最大的E值,而在化學(xué)計(jì)量數(shù)為CuAl3附近有著最大η0值,這有可能是導(dǎo)致Cu-Al合金在圖9相應(yīng)的成分下出現(xiàn)極值的原因.

圖9 Cu-Al合金體系等溫運(yùn)動(dòng)學(xué)黏度[37]Fig.9.The isotherms of kinematic viscosity of melts of Cu-Al system[37].

3.3 多元合金體系的黏度

根據(jù)Inoue提出的制備塊體非晶的三原則[38],我們已經(jīng)知道制備大塊非晶通常需具備的三個(gè)經(jīng)驗(yàn)條件:1)合金由三種或以上組元組成;2)主要組元的混合焓為大的負(fù)值,即?H<0;3)組元之間的原子尺寸比要大于12%.多元合金體系的黏度往往比單質(zhì)或二元合金體系黏度值要大.Zhou等[34]在報(bào)道CuZr基熔體的液—液相變時(shí),測(cè)量了Cu46Zr46Al8和Cu48Zr48Al4的黏度.實(shí)驗(yàn)結(jié)果表明,Cu46Zr46Al8和Cu48Zr48Al4均在1350 K左右存在明顯的黏度突變,DSC掃描結(jié)果表明在1350 K左右有明顯的放熱峰存在,與黏度突變的溫度區(qū)間較好地符合,兩者均證明了液-液相變的發(fā)生.為了進(jìn)一步探究黏度突變的動(dòng)力學(xué)本質(zhì),他們用經(jīng)驗(yàn)公式[23,39]計(jì)算了高溫轉(zhuǎn)變區(qū)和低溫轉(zhuǎn)變區(qū)流體團(tuán)簇的尺寸,發(fā)現(xiàn)Cu46Zr46Al8在高溫下團(tuán)簇尺寸為7.62 nm,在低溫下的團(tuán)簇為6.12 nm;Cu48Zr48Al4在高溫下團(tuán)簇尺寸為9.58 nm,在低溫下的團(tuán)簇為7.69 nm.流體團(tuán)簇的尺寸隨著溫度的降低,逐漸減小,通過(guò)流體團(tuán)簇尺寸的變化說(shuō)明,在降溫的過(guò)程中,團(tuán)簇結(jié)構(gòu)變化是從高密度有序態(tài)變?yōu)榈兔芏葻o(wú)序態(tài)的一個(gè)過(guò)程.

郭海東[19]報(bào)道了Co68.25Fe4.25Si15.25B12.25熔體的黏度[19],探究了Co-Fe-Si-B合金升溫和降溫過(guò)程中黏度值的差異.實(shí)驗(yàn)發(fā)現(xiàn)升溫和降溫過(guò)程中黏度數(shù)據(jù)有著明顯不同,如圖11所示.在升溫過(guò)程中,Co-Fe-Si-B合金熔體是一個(gè)逐漸變化的過(guò)程,當(dāng)溫度達(dá)到1240?C時(shí),合金Co-Fe-Si-B熔體由低溫結(jié)構(gòu)轉(zhuǎn)變?yōu)楦邷亟Y(jié)構(gòu);在冷卻過(guò)程中,當(dāng)冷卻到1240?C時(shí)低溫結(jié)構(gòu)又再次恢復(fù).這說(shuō)明液態(tài)結(jié)構(gòu)轉(zhuǎn)變的可逆性.在降溫過(guò)程中,其黏度值明顯比升溫過(guò)程同一溫度下的黏度值低,這表明Co-Fe-Si-B熔體的高溫結(jié)構(gòu)具有一定的穩(wěn)定性,結(jié)構(gòu)轉(zhuǎn)變的進(jìn)行速度可能要比外界條件(如溫度)的變化速度小得多.

圖10 (a)降溫過(guò)程中Cu48Zr48Al4(紅色曲線)和Cu46Zr46Al8(黑色曲線)熔體的平衡黏度;(b)從1623 K以0.167 K/s速率降溫的DSC曲線[34]Fig.10.(a)The equilibrium viscosity(η)upon cooling for Cu48Zr48Al4(red one)and Cu46Zr46Al8(black one)melts;(b)corresponding DSC outputs in arbitrary unit(AU)upon cooling from 1623 K with a 0.167 K/s rate[34].

圖11 (a)Co68.25Fe4.25Si15.25B12.25熔體Arrhenius曲線logη-1/(10?3T);(b)升溫和降溫時(shí)Co68.25Fe4.25Si15.25B12.25熔體黏度變化曲線[19]Fig.11. (a)Co68.25Fe4.25Si15.25B12.25melt Arrhenius plot logη-1/(10?3T);(b)variation cure of viscosity of Co68.25Fe4.25Si15.25B12.25melt with temperature during heating or cooling[19].

鐵基非晶的玻璃形成能力比較差,這就限制了鐵基塊體非晶的制備與應(yīng)用.Zheng等[40]把鐵基非晶低的玻璃形成能力歸因于高溫金屬熔體的動(dòng)力學(xué)轉(zhuǎn)變(通過(guò)黏度的變化來(lái)體現(xiàn)),從流體團(tuán)簇和過(guò)冷液體脆性角度出發(fā),證實(shí)動(dòng)力學(xué)行為的轉(zhuǎn)折過(guò)程使鐵基合金過(guò)冷液體變得非常不穩(wěn)定,從而導(dǎo)致了鐵基合金低的玻璃形成能力.此外研究還發(fā)現(xiàn),動(dòng)力學(xué)轉(zhuǎn)變的程度對(duì)應(yīng)著鐵基合金不同的玻璃形成能力.如圖12所示:A1,A2,A3分別代表著x=0,1,2時(shí)(Fe75B15Si10)100?xNbx(x=0,1,2)的組分,同樣B1,B2,B3分別代表著x=0,1,2時(shí)Fe72?xCuxB20Si4Nb4(x=0,0.2,0.6)的組分. 從圖12可以看出隨著溫度的降低黏度逐漸增大,然而斜率卻發(fā)生了輕微的變化,黑色箭頭所指的溫度即為動(dòng)力學(xué)轉(zhuǎn)變的溫度.分別計(jì)算了動(dòng)力學(xué)轉(zhuǎn)變前后HTZ和LTZ的η0和E,發(fā)現(xiàn)HTZ有較大的E值,但η0比LTZ要小,而LTZ有著較小的E值和大的η0.η0和E的變化證明了熔體動(dòng)力學(xué)性質(zhì)的變化.該變化的發(fā)生對(duì)應(yīng)鐵基熔體中高溫有序態(tài)向低溫?zé)o序態(tài)的變化.

關(guān)于多元合金體系黏度的報(bào)道還有很多,如Zhang等[41]報(bào)道了Al-Cu-Mg-Si體系黏度;Jia等[42]報(bào)道了Al-(Ni)-Yb體系黏度;Gancarz和Gasior[43]報(bào)道Sn-Zn-Na熔體的黏度;Liu等[44]報(bào)道了Fe-Ni-Co三元合金體系的黏度.這些研究大都集中在通過(guò)黏度突變探究高溫金屬熔體的液-液相變等方面,但是關(guān)于黏度突變的本質(zhì)鮮有報(bào)道,液體結(jié)構(gòu)仍然無(wú)法被很好解釋.最近,Xiong等[45?50]通過(guò)實(shí)驗(yàn)和計(jì)算機(jī)模擬的手段研究了Al-Cu,Ag-Ge,Al-Ag,Ag-Ga合金原子結(jié)構(gòu)的演變過(guò)程,實(shí)驗(yàn)結(jié)果是在熔點(diǎn)以上不同溫度范圍內(nèi)通過(guò)高能X射線得到的結(jié)構(gòu)因子S(q)曲線,通過(guò)觀察S(q)曲線上峰值位置隨溫度的變化來(lái)觀察液體結(jié)構(gòu)的變化.研究發(fā)現(xiàn),S(q)曲線上峰值發(fā)生變化的溫度與G(r)曲線上峰值發(fā)生變化的溫度符合得很好.實(shí)驗(yàn)和計(jì)算機(jī)模擬的手段同時(shí)證明了液體結(jié)構(gòu)隨溫度的變化情況.

圖12 黏度隨溫度的變化 (a)A系列,(b)B系列[40]Fig.12.Variation cure of viscosity of with temperature.(a)Series A;(b)series B[40].

4 高溫金屬熔體黏度研究的前景與展望

利用高溫金屬熔體的黏度變化趨勢(shì)或規(guī)律,可以進(jìn)一步獲得過(guò)熱熔體的脆性等動(dòng)力學(xué)表征參數(shù),同時(shí)也是研究液-液相變現(xiàn)象、金屬玻璃形成液體的強(qiáng)脆轉(zhuǎn)變現(xiàn)象的重要依據(jù).作為一種由高溫金屬熔體經(jīng)快速冷卻凝固制備而成的先進(jìn)材料,金屬玻璃的制備更容易受到高溫金屬熔體黏度的直接影響.理解合金高溫金屬熔體的動(dòng)力學(xué)性質(zhì)變化規(guī)律,不僅有助于從液-固遺傳的角度理解金屬凝固過(guò)程的本質(zhì)、金屬晶化與非晶化的競(jìng)爭(zhēng)關(guān)系以及金屬材料固體性質(zhì)與其液體動(dòng)力學(xué)脆性的聯(lián)系,同時(shí)對(duì)生產(chǎn)實(shí)踐也具有重要的指導(dǎo)意義.如果把超快速凝固也作為一種先進(jìn)鑄造技術(shù)來(lái)看待,液-固遺傳性必然在金屬玻璃形成體系中同樣存在.而且,由于金屬玻璃的整體結(jié)構(gòu)更接近液體的結(jié)構(gòu),因此相比較傳統(tǒng)鑄造技術(shù)制備的金屬材料,金屬玻璃高溫金屬熔體的性質(zhì)應(yīng)該對(duì)其固體的性質(zhì)產(chǎn)生更為直接的影響.前期工作初步發(fā)現(xiàn)了高溫金屬熔體的黏度異常突變與非晶固體三階段焓異常弛豫行為的關(guān)系[51,52].具有明顯三階段焓異常弛豫行為的非晶合金,其高溫金屬熔體的黏度突變也非常明顯.而高溫金屬熔體黏度突變不明顯的玻璃形成體系,其固體焓弛豫行為具有單調(diào)性,僅僅在非晶條帶的起始弛豫激活能上可以觀察到三階段的變化.高溫金屬熔體的黏度突變程度還與玻璃形成能力有關(guān).已經(jīng)發(fā)現(xiàn),金屬過(guò)熱液體的脆性與玻璃形成能力具有反比關(guān)系.進(jìn)一步地,在發(fā)生液-液相變的金屬玻璃熔體體系中,其發(fā)生液-液相變的相變程度也與玻璃形成能力密切相關(guān)[32,40],但這種關(guān)系似乎更為復(fù)雜.在CuZr體系中,發(fā)生的黏度突變對(duì)應(yīng)明顯的放熱過(guò)程[32],而在鐵基玻璃形成體系中,其黏度的突變伴隨微弱的吸熱效應(yīng)[40].這一比較似乎可以說(shuō)明發(fā)生液-液相變后的CuZr熔體結(jié)構(gòu)更為穩(wěn)定,而鐵基則剛好相反.高溫金屬熔體黏度的異常變化同樣可以引起固體力學(xué)性質(zhì)和其他物理性質(zhì)的變化.通過(guò)對(duì)26種不同種類液體的數(shù)據(jù)分析,金屬玻璃液體的強(qiáng)脆轉(zhuǎn)變系數(shù)(即黏度隨溫度的變化率的改變)與非晶固體的泊松比成負(fù)的指數(shù)關(guān)系[53,54].目前,對(duì)金屬玻璃熔體的認(rèn)識(shí)還不夠深入,對(duì)其黏度突變伴隨的微觀結(jié)構(gòu)以及其他物性的變化還不明確.對(duì)于高溫金屬熔體黏度與液態(tài)金屬鍵的關(guān)系,國(guó)內(nèi)外研究較少.鍵型數(shù)的改變與黏度的改變有無(wú)因果關(guān)系,團(tuán)簇的哪些特征對(duì)黏度突變產(chǎn)生明顯的影響,仍缺乏定量與定性的分析.飛速發(fā)展的中子散射技術(shù)、高能X射線衍射技術(shù)等先進(jìn)實(shí)驗(yàn)手段的應(yīng)用,輔助以計(jì)算機(jī)模擬,以原子徑向密度分布函數(shù)、雙體分布函數(shù)、均方位移等多種結(jié)構(gòu)參數(shù)為橋梁,結(jié)合應(yīng)用鍵對(duì)分析技術(shù)等,使深入探索高溫金屬熔體黏度突變的微觀結(jié)構(gòu)來(lái)源成為可能.

另外,作為體現(xiàn)材料在不同溫度、壓力、成分下組織和結(jié)構(gòu)變化的依據(jù),相圖可以直觀體現(xiàn)不同成分下的合金在不同溫度范圍內(nèi)所處的組織狀態(tài).高溫金屬熔體的黏度變化是由于液體結(jié)構(gòu)變化所致,將相圖與高溫金屬熔體的黏度結(jié)合在一起,或許是探究高溫金屬熔體的黏度突變或液-液相變本質(zhì)的另一條出路[32,34?36].相圖與高溫金屬熔體黏度之間的內(nèi)在關(guān)系尚未明確,還需要進(jìn)行更為系統(tǒng)和深入的研究.以二元合金成分與黏度之間的關(guān)系為例,MolewynHoghes提出了一種理論[55],其表達(dá)式為

式中μA為合金的運(yùn)動(dòng)學(xué)黏度;x1,x2分別為兩種成分的摩爾分?jǐn)?shù);?μ為轉(zhuǎn)換能;μ1,μ2分別為兩種成分的運(yùn)動(dòng)學(xué)黏度.但用該理論計(jì)算的結(jié)果與實(shí)驗(yàn)結(jié)果相比較,存在著較大誤差.桂滿昌[56]通過(guò)對(duì)Al-Si合金的黏度研究發(fā)現(xiàn):在共晶點(diǎn)左邊亞共晶區(qū)域內(nèi),隨著含硅量的增加,黏度值下降;在共晶點(diǎn)右邊的過(guò)共晶區(qū)域,隨著含硅量的增加,黏度值增高.此外,在亞共晶鋁硅合金中,加入變質(zhì)元素Sr使合金的黏度增加,同時(shí)使黏度突變量減小,但并不改變黏度突變點(diǎn)的溫度;在過(guò)共晶鋁硅合金中加入變質(zhì)元素P后,使合金第一個(gè)黏度突變點(diǎn)溫度升高,并使合金黏度值增加.除此之外,其他學(xué)者也提出了一些有關(guān)二元合金成分與黏度之間關(guān)系的理論或半經(jīng)驗(yàn)的模型[57?61].但是這些模型往往只適用于某些特定的二元合金體系,當(dāng)換用其他體系時(shí),存在著比較大的誤差.由此可見(jiàn)黏度與高溫金屬熔體成分之間的關(guān)系比較復(fù)雜,對(duì)其關(guān)系的確立必須建立在大量準(zhǔn)確的黏度數(shù)據(jù)基礎(chǔ)之上.

[1]Han X F 2005M.S.Dissertation(Jinan:Shandong University)(in Chinese)[韓秀峰 2005碩士學(xué)位論文 (濟(jì)南:山東大學(xué))]

[2]Angell C A 1985J.Non-Cryst.Solids73 1

[3]Bian X F,Sun B A,Hu L N,Jia Y B 2005Phys.Lett.A335 61

[4]Meng Q G,Zhou J K,Zheng H X,Li J G 2006Scr.Mater.54 777

[5]Hu L N,Bian X F 2003Chin.Sci.Bull.48 2393(in Chinese)[胡麗娜,邊秀房 2003科學(xué)通報(bào) 48 2393]

[6]Hu L N,Zhang C Z,Yue Y Z,Bian X F 2010Chin.Sci.Bull.55 115(in Chinese)[胡麗娜,張春芝,岳遠(yuǎn)征,邊秀房2010科學(xué)通報(bào)55 115]

[7]Books R F,Dinsdale A T,Quested P N 2005Meas.Sci.Technol.16 354

[8]Dinsdale A T,Quested P N 2004J.Mater.Sci.39 7221

[9]Torklep K,Oye H A 1979J.Phys.E12 875

[10]SatoY,KamedaY,NagasawaT,SakamotoT,Moriguchi S,Yamamura T,Waseda Y 2003J.Cryst.Growth249 404

[11]Kehr M,Hoyer W,Egry I 2007Int.J.Thermophys.28 1017

[12]Nunes V M B,Santos F J V,de Castro C A N 1998Int.J.Thermophys.19 427

[13]Schenck H,Frohberg M G,Hoffmann K 1963Steel Res.Int.34 93

[14]Emadi D,Gruzleski J E,Toguri J M 1993Metall.Trans.B24 1055

[15]Xu Y P,Zhao X L,Yan T L 2017Chin.Phys.B26 036601

[16]Wu Y Q,Bian X F,Mao T,Li X L,Li T B,Wang C D 2006Phys.Lett.A361 265

[17]Sun C,Geng H,Liu J,Gneg H,Yang Z 2004Phys.Meas.1 16

[18]Wang C Z 2017M.S.Dissertation(Jinan:Shandong University)(in Chinese)[王春震 2017碩士學(xué)位論文 (濟(jì)南:山東大學(xué))]

[19]Guo H D 2008M.S.Dissertation(Harbin:Harbin Institute of Technology)(in Chinese)[郭海東 2008碩士學(xué)位論文(哈爾濱:哈爾濱工業(yè)大學(xué))]

[20]Sun C J,Geng H R,Zhang N,Teng X Y,Ji L L 2008Mater.Lett.62 73

[21]Mao T,Bian X F,Morioka S,Wu Y Q,Li X L,Lü X Q 2007Phys.Lett.366 155

[22]Sun M H,Geng H R,Bian X F,Liu Y 2000Acta Metal.Sin.36 1134(in Chinese)[孫民華,耿浩然,邊秀房,劉燕2000金屬學(xué)報(bào)36 1134]

[23]Wang L,Bian X F,Liu J T 2004Phys.Lett.A326 429

[24]Ofte D,Wittenberg L J 1963Trans.Metall.Soc.Aime.227 706

[25]Rothwell E 1961J.Inst.Metals90 389

[26]Gebhardt E,Kostlin K 1957Z.Metallkd.48 636

[27]Schenck H,Frohberg M G,Hoffmann K 1963Arch.Eisenhuettenw.34 93

[28]Cavalier G 1963Compt.Rend.256 1308

[29]Kaplun A B,Avaliani M 1977High Temp.15 259

[30]Nikolaev B,Vollmann J 1996J.Non-Cryst.Solids208 145

[31]Martin-Garin L,Martin-Garin R,Despre P 1978J.Less Common.Met.59 1

[32]Zhao X,Wang C Z,Zheng H J,Tian Z A,Hu L N 2017Phys.Chem.Chem.Phys.19 15962

[33]Zhao Y,Hou X X,Bian X F 2008Mater.Lett.62 3542

[34]Zhou C,Hu L N,Sun Q J,Bian X F,Yue Y Z 2013Appl.Phys.103 171904

[35]Ning S,Bian X F,Ren Z F 2010Phys.B:Condens.Matter405 3633

[36]Mao T,Bian X F,Xue X Y,Zhang Y N,Guo J,Sun B A 2007Phys.B:Phys.Condens.Matter387 1

[37]Konstantinova N Y,Popel’P S,Yagodin D A 2009High Temp.47 336

[38]Inoue A,Takeuchi A 2010Int.J.Appl.Glass Sci.1 273

[39]Wang L,Liu J T 2004Phys.Lett.A328 241

[40]Zheng H J,Hu L N,Zhao X,Wang C Z,Sun Q J,Wang T,Hui X D,Yue Y Z,Bian X F 2017J.Non-Cryst.Solids471 120

[41]Zhang F,Du Y,Liu S H,Jie W Q 2015Comput.Coupling Phase Diagrams Thermochem.49 79

[42]Jia R,Bian X F,Lu X Q,Song K K,Li X L 2010Phys.Mech.Astron.53 390

[43]Gancarz T,Gasior W 2016Fluid Phase Equilib.418 57

[44]Liu Y H,Lu X W,Bai C G,Lai P S,Wang J S 2015J.Ind.Eng.Chem.30 106

[45]Xiong L H,Lou H B,Wang X D,Debela T T,Cao Q P,Zhang D X,Wang S Y,Wang C Z,Jiang J Z 2014Acta Mater.68 1

[46]Xiong L H,Chen K,Ke F S,Lou H B,Yue G Q,Shen B,Dong F,Wang S Y,Chen L Y,Wang C Z,Ho K M,Wang X D,Lai L H,Xiao T Q,Jiang J Z 2015Acta Mater.92 109

[47]Xiong L H,Yoo H,Lou H B,Wang X D,Cao Q P,Zhang D X,Cao Q P,Zhang D X,Jian J Z,Xie H L,Xiao T Q,Jeon S,Lee G M 2015J.Phys.:Condens.Matter27 035102

[48]Xiong L H,Guo F M,Wang X D,Cao Q P,Zhang D X,Ren Y,Jiang J Z 2017J.Non-Cryst.Solids459 160

[49]Xiong L H,Wang X D,Cao Q P,Zhang D X,Xie H L,Xiao T Q,Jiang J Z 2017J.Phys.:Condens.Matter29 035101

[50]Su Y,Wang X D,Yu Q,Cao Q P,Ruett U,Zhang D X,Jiang J Z 2018J.Phys.:Condens.Matter30 015402

[51]Wang C W,Hu L N,Chen W,Tong X,Zhou C,Sun Q J,Hui X D,Yue Y Z 2014J.Phys.Chem.141 164507

[52]Hu L N,Zhou C,Zhang C Z,Yue Y Z 2013J.Phys.Chem.138 174508

[53]Sun Q J,Hu L N,Zhou C,Zheng H J,Yue Y Z 2015J.Phys.Chem.Lett.143 164504

[54]Sun Q J,Zhou C,Yue Y Z,Hu L N 2014J.Phys.Chem.Lett.5 1170

[55]Iida T,Roderick I L,1993The Properties of Liquid Metals(Oxford:University Press)pp147–199

[56]Gui M C 1994Ph.D.Dissertation(Harbin:Harbin Institute of Technology)(in Chinese)[桂滿昌 1994博士學(xué)位論文(哈爾濱:哈爾濱工業(yè)大學(xué))]

[57]Iidia T,Ueda M,Morita Z 1976Tetsu to Hagane62 1169

[58]Morita Z,Iida T,Ueda M 1997Inst.Phys.Conf.Ser.30 600

[59]Djemili B,Martin-Garin L,Hicter P 1980J.Phys.Colloq.C841 363

[60]Enskog D 1922Arkiv.Mth.Astron.Fys.16 16

[61]Tham M K,Gubbins K E 1971J.Chem.Phys.55 268