紙張經(jīng)典滲吸模型相互關(guān)系仿真分析與理論推導(dǎo)

劉國(guó)棟， 成囡囡， 王　月, 馬　欣， 王　哲

(1.陜西科技大學(xué) 輕工與能源學(xué)院 陜西省造紙技術(shù)及特種紙品開發(fā)重點(diǎn)實(shí)驗(yàn)室, 陜西 西安　710021; 2.陜西青年職業(yè)學(xué)院, 陜西 西安　710068)

?

紙張經(jīng)典滲吸模型相互關(guān)系仿真分析與理論推導(dǎo)

劉國(guó)棟1， 成囡囡1， 王月2, 馬欣1， 王哲1

(1.陜西科技大學(xué) 輕工與能源學(xué)院 陜西省造紙技術(shù)及特種紙品開發(fā)重點(diǎn)實(shí)驗(yàn)室, 陜西 西安710021; 2.陜西青年職業(yè)學(xué)院, 陜西 西安710068)

摘要：針對(duì)紙張滲吸過程中不同毛細(xì)力學(xué)系統(tǒng)支配下所建立的經(jīng)典滲吸模型：純慣性力滲吸模型、Bosanquet滲吸模型和Lucas-Washburn滲吸模型.在統(tǒng)一的物理參數(shù)下，通過模型滲吸仿真分析和數(shù)學(xué)理論推導(dǎo)的方式，揭示了經(jīng)典滲吸模型之間存在的相互關(guān)系.同時(shí)，在滲吸分析的基礎(chǔ)上，討論了經(jīng)典滲吸模型應(yīng)用的特定情況，為紙張滲吸行為的精確描述提供切實(shí)的理論依據(jù).

關(guān)鍵詞：純慣性力滲吸模型； Bosanquet滲吸模型； Lucas-Washburn滲吸模型

0引言

印刷油墨等流體在紙張上的呈色表現(xiàn)是通過流體與紙張的相互作用來完成，該過程通常稱為流體(油墨)在紙張中的滲透[1]或者滲吸[2].在不考慮外在壓力的情況下，油墨流體在紙張中的傳輸過程主要表現(xiàn)為毛細(xì)力學(xué)系統(tǒng)的自發(fā)滲吸.所謂自發(fā)滲吸(簡(jiǎn)稱滲吸)是多孔介質(zhì)在毛細(xì)管力驅(qū)動(dòng)下自發(fā)地吸入某種潤(rùn)濕液體的過程[3].當(dāng)油墨流體與紙張發(fā)生作用時(shí)，在毛細(xì)力學(xué)系統(tǒng)的支配下，油墨的流相層(連接料等)滲吸到紙張中，油墨的顏料顆粒在紙張表面附著固化后形成圖文.因此，紙張滲吸特性的優(yōu)劣直接決定著滲吸過程中油墨的傳輸及滲吸后形成的印刷質(zhì)量.紙張滲吸行為描述時(shí)，通常將紙張介質(zhì)等效為毛細(xì)管結(jié)構(gòu)來研究紙張的滲吸特性，常見的經(jīng)典滲吸模型有Lucas-Washburn模型(簡(jiǎn)稱L-W模型)[4]，Bosanquet模型[5,6]和純慣性力模型[7]三種模型，這三種模型分別在不同的力學(xué)系統(tǒng)中描述流體在紙張等多孔介質(zhì)中的滲吸.在上述模型的應(yīng)用中，幾乎所有中文文獻(xiàn)都在利用L-W模型描述油墨等流體在紙張中的滲吸行為.但是Schoelkopf[8]和Ridgway等[9]在研究涂布紙涂層材料滲吸時(shí)，觀察到了L-W模型描述滲吸與實(shí)驗(yàn)滲吸的偏差，并分析了造成這種偏差的原因可能在于流體滲吸過程中的慣性力，提出了應(yīng)該利用考慮慣性力因素的Bosanquet模型來描述涂布紙涂層材料的滲吸.與此同時(shí)，Quere[10]通過實(shí)驗(yàn)觀察證實(shí)了流體在毛細(xì)管中的慣性流動(dòng)，并提出了純慣性力支配的滲吸模型.劉國(guó)棟[11]在研究涂層材料滲吸時(shí)，證實(shí)了慣性力在滲吸初期存在的積極作用，一定程度上證明了純慣性力模型在描述滲吸過程初期階段的可行性.

因此本文針對(duì)經(jīng)典滲吸模型優(yōu)勢(shì)及不足，在統(tǒng)一滲吸物理參數(shù)下，對(duì)比分析及驗(yàn)證上述經(jīng)典滲吸模型在描述紙張滲吸行為的具體表現(xiàn)及可行性，為紙張滲吸特性的準(zhǔn)確描述提供可靠的理論依據(jù).

1表征紙張滲吸行為的經(jīng)典模型

1.1純慣性力支配的滲吸模型

當(dāng)流體與紙張接觸的初期，由于滲吸流體質(zhì)量很小，對(duì)應(yīng)的自身重力和粘性阻力也很小，在滲吸力學(xué)系統(tǒng)中，可將自身重力和粘性阻力忽略不計(jì)，因此起主要支配作用的是毛細(xì)力和慣性力.針對(duì)此種物理情形，滲吸模型[10]的表達(dá)如(1)式所示.

(1)

通過求解非線性微分方程，可得紙張滲吸等效毛細(xì)管滲吸的表達(dá)式，如式(2)所示：

(2)

其中h為滲吸的深度，θ為接觸角，γ為流體表面張力，R為毛細(xì)管等效半徑，ρ為流體的密度.從上式(2)可以看出，滲吸深度隨滲吸時(shí)間是線性變化的關(guān)系.

1.2Bosanquet滲吸模型

在滲吸過程中，隨著滲吸流體的進(jìn)一步上升，吸收量不斷增加，與之對(duì)應(yīng)的自身重力和粘性力也變得越來越明顯，與毛細(xì)驅(qū)動(dòng)力和慣性力一起共同作用滲吸過程.考慮到在該階段初期的自身重力仍較小，可以忽略不計(jì).因此Bosanquet[5]提出了忽略流體重力因素的滲吸模型(η為流體粘度)，如式(3)所示.

(3)

求解(3)可得

(4)

其中

從上述公式(4)中可以看出，Bosanquet滲吸模型所描述滲吸量與滲吸時(shí)間并非成簡(jiǎn)單的線性關(guān)系，其計(jì)算過程略顯復(fù)雜.

1.3Lucas-Washburn滲吸模型

隨著滲吸的進(jìn)行，流體粘性阻力不斷增大，與毛細(xì)驅(qū)動(dòng)壓力的相互抵消作用也越來越明顯，凈驅(qū)動(dòng)力越來越小，同時(shí)流體流動(dòng)的慣性力也越來越小.因此在忽略流體自身重力和慣性效應(yīng)下，受粘性力和毛細(xì)力共同支配Lucas-Washburn滲吸模型，如公式(5)所示(其中η為流體粘度).

(5)

對(duì)比分析上述三種經(jīng)典滲吸模型的建立過程及表達(dá)式，可以看出上述模型分別建立在不同的力學(xué)支配系統(tǒng)下，相關(guān)文獻(xiàn)[12-14],也分別在各自的物理情形和參數(shù)下，描述了上述滲吸模型在紙張滲吸表征中的應(yīng)用.但是究竟紙張的滲吸過程確切適合哪種模型及受何種力學(xué)系統(tǒng)支配卻一直未見相關(guān)報(bào)道，仍需進(jìn)一步分析對(duì)比.

2經(jīng)典滲吸模型相互關(guān)系的數(shù)學(xué)意義證明

對(duì)于Bosanquet滲吸方程而言，其表現(xiàn)形式正如公式(4)所示：

在初始的滲吸過程中，當(dāng)at<<1時(shí)，利用泰勒級(jí)數(shù)展開公式可得[15]

(6)

將式(6)代入到(4)可得

(7)

其中h1是起始位置的滲吸深度，h2是任意時(shí)刻t的滲吸深度.

(8)

當(dāng)滲吸的起始位置x1=0時(shí)，公式(8)將變?yōu)?/p>

(9)

公式(9)與純慣性力階段描述的滲吸模型是一致的，因此在初期的滲吸過程中(即當(dāng)at<<1時(shí))，純慣性力滲吸和Bosanquet滲吸是等效的.

另外，將公式(4)進(jìn)行如下的變形,

(10)

由于

(11)

因此

即滲吸過程中，當(dāng)物理參數(shù)：

(12)

對(duì)于一些高密度和低粘度的流體，公式(12)所描述的物理情形很容易達(dá)到.而這種流體常常用于油墨中的連接料，如正十六烷、礦物油、乙二醇、亞麻籽油和異三十烷等，都屬于高密度和低粘度的油墨流體.因此從式(10)～(12)可以看出，Bosanquet滲吸方程就演變?yōu)長(zhǎng)ucas-Washburn方程，即在滲吸過程中Bosanquet滲吸行為會(huì)逐漸轉(zhuǎn)換成Lucas-Washburn滲吸.

3經(jīng)典滲吸模型的對(duì)比分析

結(jié)合上述推導(dǎo)的經(jīng)典模型的相互關(guān)系，進(jìn)一步確定紙張滲吸過程中，基于何種滲吸模型及毛細(xì)力學(xué)支配系統(tǒng)，需要在統(tǒng)一的物理參數(shù)下對(duì)上述模型進(jìn)行對(duì)比仿真分析.由于正十六烷和礦物油作為常見的水基油墨和脂性油墨的載體已經(jīng)得到了普遍的應(yīng)用，本文以正十六烷(Sigma-Aldrich Chemie GmbH，D-89555，斯坦海姆，德國(guó))和礦物油(PKWF 4/7 af new，Haltermann，德國(guó))作為滲吸流體，具體的物理參數(shù)見表1所示.紙張介質(zhì)的毛細(xì)管等效半徑R分別設(shè)定為1μm，500μm，1 mm三種.紙張等多孔材料等效毛細(xì)管對(duì)正十六烷和礦物油的潤(rùn)濕接觸角為0[16].

表1　滲吸流體的物理性能參數(shù)[8]

基于Matlab計(jì)算軟件平臺(tái)，對(duì)上述模型進(jìn)行編程實(shí)現(xiàn)，為了更清楚的對(duì)比分析滲吸模型存在的相互關(guān)系，滲吸過程分別用較短的滲吸時(shí)間軸和較長(zhǎng)的滲吸時(shí)間軸來表示，較短的時(shí)間軸可以清楚看出滲吸初期的變化，較長(zhǎng)的時(shí)間軸可以看出滲吸后期的變化，如圖1～12所示.其中，圖1、圖3、圖5、圖7、圖9、圖11為滲吸模型滲吸初期的對(duì)比變化結(jié)果，圖2、圖4、圖6、圖8、圖10、圖12為滲吸模型后期的對(duì)比結(jié)果.

(1)毛細(xì)管等效半徑R=1μm時(shí)的滲吸模型仿真結(jié)果見圖1～4.

圖1　滲吸初期對(duì)比圖(礦物油，R=1 μm)

圖2　滲吸后期對(duì)比圖(礦物油，R=1 μm)

圖3　滲吸初期對(duì)比圖(正十六烷R=1 μm)

圖4　滲吸后期對(duì)比圖(正十六烷，R=1 μm)

(2)毛細(xì)管等效半徑R=500μm時(shí)的滲吸模型仿真結(jié)果見圖5～8.

圖5　滲吸初期對(duì)比圖(礦物油，R=500 μm)

圖6　滲吸后期對(duì)比圖(礦物油，R=500 μm)

圖7　滲吸初期對(duì)比圖(正十六烷，R=500 μm)

圖8　滲吸后期對(duì)比圖(正十六烷，R=500 μm)

(3)毛細(xì)管等效半徑R=1 mm時(shí)的滲吸模型仿真結(jié)果見圖9～12.

圖9　滲吸初期對(duì)比圖(礦物油，R=1 mm)

圖10　滲吸后期對(duì)比圖(礦物油，R=1 mm)

圖11　滲吸初期對(duì)比圖(正十六烷，R=1 mm)

圖12　滲吸后期對(duì)比圖(正十六烷，R=1 mm)

從圖1到圖12可以看出，三種不同力學(xué)系統(tǒng)支配下建立的滲吸模型，的確存在著一定相互關(guān)系.在滲吸的初期階段，當(dāng)R=1μm且滲吸流體為礦物油時(shí)，Bosanquet滲吸機(jī)制與純慣性力滲吸是重合的(如圖1所示)，說明了在Bosanquet滲吸模型的初期，發(fā)生了純慣性力支配的滲吸行為，其模型力學(xué)系統(tǒng)中流體的粘性力并沒有起到實(shí)際的作用.同時(shí)也說明了該階段滲吸量(滲吸體積或者滲吸高度)與滲吸時(shí)間標(biāo)度是線性變化的.

隨著滲吸的進(jìn)一步發(fā)生，發(fā)現(xiàn)Bosanquet滲吸行為又與Lucas-Washburn滲吸重合的(圖2).說明了隨著滲吸的進(jìn)行，Bosanquet滲吸模型中所強(qiáng)調(diào)的慣性力被消散，Bosanquet滲吸逐漸轉(zhuǎn)化成以毛細(xì)力和粘性力支配的Lucas-Washburn滲吸，也說明了Bosanquet滲吸模型的后期，慣性力沒有起到實(shí)際的作用，是可以被忽略的.也說明了本文第二部分推導(dǎo)中，Bosanquet方程和Lucas-Washburn方程等效時(shí)，是發(fā)生在滲吸的后期階段.

當(dāng)R=500μm和1 mm時(shí)，滲吸情況的具體表現(xiàn)跟R=1μm時(shí)完全相同，證明了滲吸初期慣性力起到的積極作用；隨著滲吸的進(jìn)行，滲吸支配力也慢慢過渡到毛細(xì)力和粘性力為主的力學(xué)系統(tǒng)中.對(duì)于滲吸流體正十六烷，獲得的滲吸仿真結(jié)果也同樣與礦物油的結(jié)果完全一致(筆者也驗(yàn)證了在其他不同等效毛細(xì)管參數(shù)下模型的對(duì)比分析結(jié)果，與上述的情況及變化趨勢(shì)完全一致，限于篇幅，這里不在重復(fù)羅列相關(guān)變化圖表).

另外，隨著毛細(xì)管直徑的增大(由1μm到500μm，再到1 mm)，Bosanquet機(jī)制階段與純慣性力的重合時(shí)間會(huì)逐漸增加或者延長(zhǎng)，對(duì)于礦物油而言，重合的時(shí)間終點(diǎn)分別為0.8×10-8s，0.004 s，0.012 s；對(duì)于正十六烷重合的時(shí)間終點(diǎn)分別為1.4×10-8s，0.005 s，0.015 s.即隨著紙張介質(zhì)的等效毛細(xì)管半徑越大，Bosanquet機(jī)制與純慣性力機(jī)制所描述滲吸過程的一致性越長(zhǎng).

對(duì)于滲吸的后期階段，Bosanquet機(jī)制階段和Lucas-Washburn滲吸階段重合時(shí)間也會(huì)隨著等效毛細(xì)管直徑的增大而漸漸推遲，即紙張介質(zhì)的等效毛細(xì)管半徑越小，Bosanquet滲吸與Lucas-Washburn滲吸重合的時(shí)間越早，滲吸過程的一致性越長(zhǎng).

縱觀整個(gè)仿真結(jié)果，從另外一個(gè)方面也說明了，Bosanquet滲吸機(jī)制包括了純慣性力階段和Lucas-Washburn滲吸行為.另外，當(dāng)?shù)刃?xì)管的直徑越小時(shí)，純慣性力階段在滲吸過程中的作用就越短.對(duì)于特定的等效毛細(xì)管直徑，純慣性力的滲吸就可以忽略不計(jì)，這也是有時(shí)在特定的材料和實(shí)驗(yàn)條件下，Lucas-Washburn滲吸能達(dá)到與實(shí)驗(yàn)結(jié)果比較匹配的重要原因之一.

4結(jié)論

通過模型的相互關(guān)系的數(shù)學(xué)推導(dǎo)和仿真分析可以看出，雖然經(jīng)典的純慣性力滲吸模型、Bosanquet滲吸模型及Lucas-Washburn滲吸模型，是基于不同的毛細(xì)力學(xué)支配系統(tǒng)建立的，但是，他們?cè)谝欢ㄎ锢砬樾蜗?，存在著確定的相關(guān)關(guān)系.

(1)Bosanquet滲吸模型包括了純慣性力階段的滲吸行為和Lucas-Washburn滲吸行為，即純慣性力滲吸和Lucas-Washburn滲吸是Bosanquet描述滲吸過程的兩個(gè)特殊階段，他們分別描述了在Bosanquet滲吸機(jī)制下滲吸的初始和后期滲吸的兩個(gè)重要階段.

(2)當(dāng)紙張介質(zhì)的等效毛細(xì)管半徑大時(shí)，Bosanquet機(jī)制與純慣性力機(jī)制所描述滲吸過程的一致性越長(zhǎng)，在強(qiáng)調(diào)初期滲吸的情況下，表達(dá)形式簡(jiǎn)單的純慣性力滲吸模型可以代替復(fù)雜Bosanquet滲吸模型.

(3)當(dāng)紙張介質(zhì)的等效毛細(xì)管半徑小時(shí)，Bosanquet滲吸與Lucas-Washburn滲吸重合的時(shí)間越早，滲吸過程的一致性越長(zhǎng)，在忽略短暫慣性效應(yīng)的情況下，表達(dá)形式簡(jiǎn)單的Lucas-Washburn滲吸可以代替復(fù)雜Bosanquet滲吸模型.

(4)當(dāng)需要描述整個(gè)階段的滲吸時(shí)，Bosanquet滲吸是三種模型中最佳的表征方案.Bosanquet滲吸完整的包括了滲吸過程存在的慣性力、毛細(xì)力和粘性力，但是由于Bosanquet滲吸模型的表達(dá)式及求解過程較其他模型復(fù)雜得多，一定程度上限制了該模型的實(shí)際應(yīng)用.

參考文獻(xiàn)

[1] Li Y.,He B.Characterization of ink pigment penetration and distribution related to surface topography of paper using confocal laser scanning microscopy[J].BioResources,2011,6(3):2 690-2 702.

[2] Schoelkopf J.,Gane P.,Ridgway C.,et al.Imbibition behaviour of offset inks:I gravimetric determination of oil imbibition rate into pigmented coating structures[J].Tappi Journal,2003,2(6):9-13.

[3] 蔡建超，郁伯銘．多孔介質(zhì)自發(fā)滲吸研究進(jìn)展[J].力學(xué)進(jìn)展,2012，42(6)：735-753．

[4] Washburn E.W.The dynamics of fluid flow[J].Phys.Rev.,1921,17:273-283.

[5] Bosanquet C.H.On the flow of liquids into capillary tubes[J].Philos.Mag,1923,45(267):525-531.

[6] Szekely J.,Neumann A.W.,Chuang Y.K.The rate of capillary penetraion and the applicability of the washburn equation[J].J.Colloid Interface Sci.,1971,35(2):273-278.

[7] Quere D. Inertial capillarity[J].Europhys Lett.,1997,39(5):533-538.

[8] Schoelkopf J.,Gane P.,Ridgway C.,et al.Practical observation of deviation from Lucas-washburn scaling in porous media[J].Colloid Surf.A.Physicochem.Eng.Asp.,2002,206:445-454.

[9] Ridgway C.J.,Gane P.A.C.,Schoelkopf J.Modified calcium carbonate coatings with rapid absorption and extensive liquid uptake capacity[J].Colloid Surf.A.Physicochem.Eng.Asp.,2004,236:91-102.

[10] Quere D.Wetting and roughness[J].Annu.Rev.Mater.Res.，2008,38:71-99.

[11] 劉國(guó)棟．造紙涂層材料自滲吸行為機(jī)制及理論建模研究[D]．西安：陜西科技大學(xué)，2014.

[12] Liu G.,Zhang M.,Ridgway C.,et al.Pore wall rugosity:The role of extended wetting contact line length during spontaneous liquid imbibition in porous media[J].Colloid Surf.A.Physicochem.Eng.Asp.,2014,443:286-295.

[13] Poulin N.,Tangup P.A.Numerical and physical modeling of the permeability of paper to CMC and coating liquids[J].The Canadian Journal of Chemical Engineering,1997,75:949-955.

[14] Liu G.,Zhang M.,Ridgway C.,et al.Spontaneous inertial lmbibition in porous media using a fractal representation of pore wall rugosity[J].Transp Porous Med,2014,104(1):231-251.

[15] Schoelkopf J.,Gane P.A.C.,Ridgway C.J.,et al, Influence of inertia on liquid absorption into paper coating structures[J].Nordic Pulp and Paper Research J.,2000,15(5):422-430.

[16] E.Chibowski,L.Holysz.On the use of Washburn′s equation for contact angle determination[J].J.Adhes.Sci.Technol.,1997,11(10):1 289-1 301.

【責(zé)任編輯：蔣亞儒】

The exploration of interaction relationship of classical imbibition models of paper-based media

LIU Guo-dong1, CHENG Nan-nan1, WANG Yue2， MA Xin1, WANG Zhe1

(1.College of Light Industry and Energy, Shaanxi Province Key Laboratory of Papermaking Technology and Specialty Paper, Shaanxi University of Science & Technology, Xi′an 710021, China; 2.Shaanxi Youth Vocational College, Xi′an 710068, China )

Abstract:Aiming at the different imbibition models of inertial,Bosanquet and Lucas-Washburn when charactering the imbibition process of paper-based materials.The interaction relationship between these models is discovered by means of specific imbibition analysis and theoretical mathematic method based on the uniform physical parameters and conditions. On the basis of the imbibition analysis of these models,the specific and reasonable application situations are discussed further eventually,which can provide the theoretical evidences to describe the imbibition performance of paper-based media.

Key words:inertial imbibition model； Bosanquet imbibition model； Lucas-Washburn imbibition model

中圖分類號(hào)：TS727+.3

文獻(xiàn)標(biāo)志碼：A

文章編號(hào)：1000-5811(2016)02-0035-06

作者簡(jiǎn)介:劉國(guó)棟(1981-)，男，陜西乾縣人，副教授，博士，研究方向：印刷油墨流體與紙張相互作用的關(guān)系

基金項(xiàng)目：國(guó)家自然科學(xué)基金項(xiàng)目(51402180); 陜西省教育廳專項(xiàng)科研計(jì)劃項(xiàng)目(14JK1108); 陜西省協(xié)同創(chuàng)新計(jì)劃項(xiàng)目(2015XT-64); 陜西省留學(xué)人員科技活動(dòng)項(xiàng)目([2014]1059);國(guó)家級(jí)大學(xué)生創(chuàng)新創(chuàng)業(yè)訓(xùn)練計(jì)劃項(xiàng)目(201510708143) ; 陜西科技大學(xué)博士科研啟動(dòng)基金項(xiàng)目( BJ15-11)

收稿日期:2015-12-01