含趨旋微生物納米流體在拉伸缸表面上的邊界層流動分析

許曉勤，黃惠

(1.福建船政交通職業(yè)學(xué)院汽車學(xué)院，福建 福州 350007； 2.福州大學(xué)機械工程及自動化學(xué)院，福建 福州 350108)

0 引言

眾所周知，水、 油、 乙二醇等傳統(tǒng)傳熱流體的低導(dǎo)熱系數(shù)是限制其傳熱性能的重要因素之一，在這些液體中添加固體納米顆?？商岣邔?dǎo)熱系數(shù).Choi和Eastman首次提出“納米流體”的概念，它由尺寸小于100 nm的金屬或金屬氧化物顆粒分散在基液中制成[1].納米流體在醫(yī)療、 化工、 交通、 電子、 制造及食品等領(lǐng)域得到廣泛應(yīng)用，是學(xué)者和工程師們研究的熱點問題之一.Sumanth等[2]研究納米流體羧基石墨烯對散熱器性能的影響，結(jié)果表明, 當(dāng)納米顆粒體積分數(shù)為0.02%，流量為5 L·min-1，進口溫度分別為40和50 ℃時，散熱器的效率分別提高27.38%和23.41%.Siricharoenpanich等[3]將納米流體用于電子元器件的冷卻系統(tǒng)，實驗表明, Ag/Fe3O4納米流體作為冷卻劑的散熱效率高于Ag納米流體或去離子水作為冷卻劑的散熱效率.Salari等[4]將納米流體用于食品工業(yè)的熱處理，發(fā)現(xiàn)與傳統(tǒng)流體相比，在極低固體顆粒濃度下使用納米流體可以顯著提高導(dǎo)熱系數(shù)和總傳熱系數(shù)，顯著縮短熱處理時間，且使用納米流體后，食品的顏色和pH值幾乎不變.值得一提的是，基液中的納米顆粒是不穩(wěn)定的，可能導(dǎo)致納米顆粒聚集，使導(dǎo)熱系數(shù)降低.

微生物是一切肉眼看不見或看不清的微小生物的統(tǒng)稱(如病毒、 細菌、 酵母菌、 單細胞藻類及原生動物等)，在自然界中廣泛存在[5].微生物體積小、 結(jié)構(gòu)簡單、 代謝能力強、 生長繁殖快、 容易培養(yǎng)、 容易變異、 有極強的環(huán)境適應(yīng)性，在環(huán)境保護、 農(nóng)業(yè)生產(chǎn)、 食品工業(yè)、 新能源開發(fā)、 醫(yī)療和制藥等方面均有廣泛的應(yīng)用[6].微生物具有的趨光或趨氧特性，使其自主向液體上層運動且聚集，形成上重下輕的不穩(wěn)定密度分層.當(dāng)這種密度分層達到一定程度時會引起流體宏觀上的流動，即生物對流[7].微生物在流體中同時受到重力矩和粘滯力矩的共同作用會產(chǎn)生一種回旋運動，這種現(xiàn)象稱為微生物的“趨旋性”.將趨旋微生物與納米流體混合，可以防止納米流體中的納米顆粒聚集，從而提高納米流體的穩(wěn)定性[8].納米流體中的生物對流還可以提高質(zhì)量傳遞和熱傳遞，為新型微流控制設(shè)備的設(shè)計開拓新思路[9].Chakraborty等[10]探討磁場和對流邊界條件對含趨旋微生物的納米流體在擴張板上流動的影響，結(jié)果表明，磁場的增強使擴張板表面的流體速度和Nusselt數(shù)減小，表面對流參數(shù)增大了微生物自遷移通量，而Peclet數(shù)的影響卻相反.Khan[11]研究含有納米顆粒和趨旋微生物的穩(wěn)態(tài)二階納米流體在薄膜中的生物對流現(xiàn)象，圖解分析各參數(shù)的影響.Kanta等[8]研究在非線性熱輻射、 化學(xué)反應(yīng)、 內(nèi)熱源和抽吸/注射效應(yīng)下，含有趨旋微生物的可變粘度納米流體在非線性拉伸板上的生物對流現(xiàn)象.

拉伸成型是把材料加工成平板、 空心柱體或特定形狀的一種工藝，材料在拉伸過程中往往需要嚴格的冷卻才能保證成品質(zhì)量.由于在金屬和工程領(lǐng)域的廣泛應(yīng)用，如金屬絲、 晶體生長、 紙張生產(chǎn)、 玻璃纖維、 聚合物擠出和食品制造等，流體通過拉伸表面的傳熱與流動現(xiàn)象成為學(xué)者們研究的熱點問題之一.Khan和Makinde[12]數(shù)值研究含趨旋微生物的水基納米流體在對流加熱拉伸板上的邊界層流動現(xiàn)象，結(jié)果表明，無量綱速度隨著浮力比和生物對流Rayleigh數(shù)的增大而減小，壁面無量綱溫度隨著對流參數(shù)的增大而增大，隨著浮力比的增大而減小.Akbar等[13]研究含趨旋微生物的納米流體通過延伸板的自然對流現(xiàn)象，探討磁場參數(shù)、 Brown運動參數(shù)和熱泳的影響.Hayat等[14]探討含趨旋微生物的磁性納米流體在拉伸缸上的流動現(xiàn)象，結(jié)果表明，缸體曲率參數(shù)越大，缸體附近溫度越低.

為解決納米顆粒聚集導(dǎo)致納米流體傳熱特性下降問題，本研究在納米流體中引入趨旋微生物，探討穩(wěn)態(tài)二維納米流體在拉伸缸表面上的邊界層流動現(xiàn)象，建立控制方程，采用打靶法和Runge-Kutta法進行數(shù)值求解，對不同物理參數(shù)下的速度場、 溫度場、 濃度場和微生物密度場的變化趨勢進行圖形化分析，并給出各物理參數(shù)的作用機制.研究結(jié)果有助于分析納米流體加入趨旋微生物后對其流動及傳熱傳質(zhì)特性的影響規(guī)律，為類似物理模型邊界層控制或?qū)嶒炁渲眠m合不同應(yīng)用場合的穩(wěn)定納米流體提供參考.

1 數(shù)學(xué)模型

考慮趨旋微生物的影響，研究納米流體在拉伸缸表面上的穩(wěn)態(tài)二維不可壓縮流動.物理模型和坐標系如圖1所示.缸體半徑為a，缸體以Vw(z)的速度向右拉伸，由于摩擦力作用，附著在拉伸缸表面的納米流體速度與拉伸缸表面一致(不考慮滑移效應(yīng))，無窮遠處納米流體速度為0(即靜止?fàn)顟B(tài))，拉伸缸浸泡在納米流體中，納米流體中加入趨旋微生物，趨旋微生物自主運動阻止納米顆粒聚集，使納米流體處于穩(wěn)定狀態(tài).

圖1所示的流動模型屬于軸對稱流動，假設(shè)缸體軸向為z軸，徑向為r軸.為了使大多數(shù)微生物存活并保持活性，選擇水基納米流體，納米顆粒為球形氧化鋁，假設(shè)納米流體穩(wěn)定且懸浮液稀釋，考慮納米流體在流動過程中遵守質(zhì)量守恒、 動量守恒、 能量守恒、 納米顆粒濃度守恒及微生物密度守恒等定律，在Buongiorno模型[15]基礎(chǔ)上建立控制方程[14]如下，

連續(xù)性方程：

(1)

r方向的動量方程：

(2)

z方向的動量方程：

(3)

能量方程：

(4)

濃度方程：

(5)

微生物密度方程：

(6)

邊界條件為：

1) 當(dāng)r=a時，

(7)

2) 當(dāng)r→∞時，

Vz→0，T→T∞，C→C∞，N→N∞

(8)

式中：r、z分別表示坐標系的徑向和軸向坐標；Vr、Vz分別表示沿r軸和z軸方向的速度分量(m·s-1)；T、C、N分別表示納米流體的溫度(K)、 納米顆粒體積分數(shù)及微生物密度(kg·m-3)；Tw、Cw、Nw分別表示缸體表面納米流體的溫度(K)、 納米顆粒體積分數(shù)及微生物密度(kg·m-3)；T∞、C∞、N∞分別表示無窮遠處納米流體的溫度(K)、 納米顆粒體積分數(shù)及微生物密度(kg·m-3)； ν為運動粘度(m2·s-1)；α為納米流體的熱擴散系數(shù)(m2·s-1)；τ=(ρpcp)/(ρfcf)為納米顆粒熱容與基液熱容之比；ρ為流體密度(kg·m-3)；c為比熱容(J·(kg·K)-1)；DB為Brown擴散系數(shù)(m2·s-1)；DT為熱泳擴散系數(shù)(m2·s-1)；Dn為微生物擴散率(m2·s-1)；b為趨化常數(shù)； ΔC=Cw-C∞；Wc為細胞游動的最大速度(m·s-1)；a為缸體半徑(m)；U0為參考速度(m·s-1)；Vw為缸體拉伸速度(m·s-1)；l為特征長度(m).

引入如下無量綱變量：

(9)

式中：ψ為流函數(shù).

則：

(10)

將式(9)、 (10)代入式(1)～(6)得：

(1+2ηγ)f?+2γf″+ff″-f′2=0

(11)

(1+2ηγ)θ″+2γθ′+Pr(1+2ηγ)(NTθ′2+NBθ′φ′)+Prfθ′=0

(12)

(13)

(1+2ηγ)φ″+2γφ′-Pe[(1+2ηγ)(φ′φ′+φφ″+σ1φ″)+2γ(φφ′+σ1φ′)]+Scfφ′=0

(14)

式中： 曲率參數(shù)γ=[νl/(U0a2)]1/2，與缸體半徑a成反比； Prandtl數(shù)Pr=ν/α，表示動量傳遞及熱量傳遞效果之比； 熱泳參數(shù)NT=τDT(Tw-T∞)/(νT∞)，用于描述在溫度梯度不為0的流體中，納米顆粒向較冷區(qū)域運動的現(xiàn)象； Brown運動參數(shù)NB=τDB(Cw-C∞)/ν，用于描述液體分子不停做無規(guī)則運動，不斷隨機撞擊懸浮微粒的現(xiàn)象； Lewis數(shù)Le=ν/DB，定義為納米流體粘度與布朗擴散系數(shù)之比； 生物對流Peclet數(shù)Pe=bWc/Dn，指在受迫運動傳質(zhì)時對流傳遞與擴散傳遞之比； 生物對流常數(shù)σ1=N∞/(Nw-N∞)； 生物對流Schmidt數(shù)Sc=ν/Dn，定義為動量擴散率和微生物擴散率之比.

邊界條件(7)、 (8)變?yōu)椋?/p>

1) 當(dāng)η=0時，

f(0)=0，f′(0)=θ(0)=φ(0)=φ(0)=1

(15)

2) 當(dāng)η→∞時，

f′(∞)=θ(∞)=φ(∞)=φ(∞)=0

(16)

局部表面摩擦系數(shù)Cfz、 局部Nusselt數(shù)Nuz、 局部Sherwood數(shù)Shz和微生物的局部密度數(shù)Nnz分別定義如下：

(17)

由相似變量(9)得：

(18)

式中：Rez=zVw/ν是局部Reynolds數(shù).

2 數(shù)值方法

控制方程(11)～(14)和邊界條件(15)、 (16)構(gòu)成兩點邊值問題，由于是高階非線性微分方程，解析法求解困難，可采用數(shù)值方法求解.先通過打靶法將其轉(zhuǎn)化為非線性微分方程的初值問題，其基本思想是將邊值問題中未知的初始值設(shè)為未知參數(shù)，從而將此問題看作是帶有未知參數(shù)的初值問題，通過對參數(shù)賦值后，使邊值問題轉(zhuǎn)化為初值問題.采用Runge-Kutta法求解初值問題，由于一次賦值通常不能得到滿足精度要求的解，通過牛頓迭代法可逐步調(diào)整初值參數(shù)，直到初值問題的解滿足邊值約束，即類似一個不斷調(diào)整射中靶心的打靶過程.Runge-Kutta法是一種在工程上應(yīng)用廣泛的高精度單步算法，它來源于泰勒公式，采用斜率近似表達微分，在積分區(qū)間多預(yù)算出幾個點的斜率，然后進行加權(quán)平均，用做下一點的依據(jù)，從而構(gòu)造出精度更高的數(shù)值積分計算方法.為了驗證打靶法和Runge-Kutta法的有效性，將Pr取不同值時-θ′(0)的計算結(jié)果與文獻[13, 16-17]的結(jié)果進行對比，結(jié)果吻合，如表1所示.

表1 Pr取不同值時的-θ′(0)(當(dāng)γ=NT=NB=0時)

3 實驗結(jié)果與討論

圖2給出曲率參數(shù)γ對速度f′(η)、 溫度θ(η)、 體積分數(shù)φ(η)和微生物密度φ(η)的影響.由圖2可知，在靠近缸體表面處(η值較小時)，速度、 溫度、 體積分數(shù)和微生物密度均隨曲率參數(shù)的增大而減??； 而在離缸體較遠處(η值較大時)，速度、 溫度、 體積分數(shù)和微生物密度均隨著曲率參數(shù)的增大而增大.因為在缸體表面附近的粘性力所產(chǎn)生的阻力比遠離缸體的粘性力所產(chǎn)生的阻力大得多，因此隨著曲率參數(shù)γ的增大，缸體半徑減小，對流體運動的阻力減小，在遠離缸體時流體速度增大.另外，缸體半徑減小，附著在缸體表面上的納米顆粒和微生物數(shù)量亦減小，而這些顆粒通過熱傳導(dǎo)進行傳熱，因此缸體表面附近的體積分數(shù)、 溫度和微生物密度均減小.

圖3給出Brown運動參數(shù)NB和熱泳參數(shù)NT對溫度θ(η)的影響.由圖3可知，Brown運動參數(shù)NB和熱泳參數(shù)NT對溫度的影響類似，NB越大，納米流體Brown運動越劇烈，使對流換熱強度增強，因此溫度越高； 從物理上講，對于較高的熱泳參數(shù)值NT，熱泳力增大，促進納米顆粒從較熱區(qū)域向較冷區(qū)域移動，導(dǎo)致熱邊界層厚度和溫度增加.納米流體的Brown運動和熱泳現(xiàn)象大大提高了傳統(tǒng)流體的傳熱特性.另外，溫度升高還會進一步降低納米流體粘度，改善納米流體流動特性，納米流體良好的散熱效果可保證最終產(chǎn)品質(zhì)量.

圖4揭示Brown運動參數(shù)NB和熱泳參數(shù)NT對納米顆粒體積分數(shù)φ(η)的影響.由圖4(a)可知，NB越大，體積分數(shù)越小.Brown運動參數(shù)NB是描述納米顆粒在基液中隨機運動劇烈程度的物理量，NB越大，顆粒無規(guī)則運動速度越快，濃度越低.NT的影響恰好相反，如圖4(b)所示.從物理上講，熱泳力是溫度梯度產(chǎn)生的，當(dāng)熱泳力增加時，納米顆粒從熱表面向外移動，使納米顆粒體積分數(shù)上升.

圖5給出Brown運動參數(shù)NB和熱泳參數(shù)NT對微生物密度φ(η)的影響.由圖5可知，Brown運動參數(shù)NB和熱泳參數(shù)NT對微生物密度的影響趨勢與對濃度的影響類似.從物理角度來看，NB越大，納米顆粒在基液中隨機運動越劇烈，納米顆粒之間、 納米顆粒與水分子之間相互碰撞，微生物運動時所受阻力增大，導(dǎo)致微生物密度下降.NT越大，熱泳力越大，促進微生物生長，因此微生物密度上升.

圖6給出生物對流Peclet數(shù)(Pe)、 Schmidt數(shù)(Sc)和生物對流常數(shù)(σ1)對微生物密度φ(η)的影響.由圖6(a)可知，Pe對φ(η)影響較大，Pe越大，φ(η)越小.從Pe的表達式(Pe=bWc/Dn)可以看出，Pe與細胞游動最大速度Wc和趨化常數(shù)b成正比，與微生物擴散率Dn成反比，Pe越大，微生物平均游動速度越快，使微生物密度及其邊界層厚度降低.隨著Sc(Sc=ν/Dn)的增加，動量擴散率與微生物擴散率比值越大，微生物的生長速度受到抑制，微生物密度下降，見圖6(b).由圖6(c)可知，生物對流常數(shù)σ1(σ1=N∞/(Nw-N∞))越大，無量綱微生物密度φ(η)(φ(η)=(N-N∞)/(Nw-N∞))越小.

圖7給出了曲率參數(shù)γ對局部表面摩擦系數(shù)f″(0)、 局部Nusselt數(shù)-θ′(0)、 局部Sherwood數(shù)-φ′(0)和運動微生物的局部密度數(shù)-φ′(0)的影響，其他物理參數(shù)的影響見表2.由圖7可知，|f″(0)|、 -θ′(0)、 -φ′(0)和-φ′(0)均隨著γ的增大而增大.意味著曲率參數(shù)越大，靠近壁面納米流體流動速度越慢，而熱擴散率、 濃度擴散率、 運動微生物遷移率均越大.由表2可知，當(dāng)Pe=Sc=1，γ由0增大到1時，|f″(0)|、 -θ′(0)、 -φ′(0)、 -φ′(0)分別增大35.3%、 57.9%、 48.0%、 33.3%；Pe和Sc對f″(0)、 -θ′(0)和-φ′(0)的影響可以忽略不計，但對-φ′(0)的影響較大，當(dāng)γ=Sc=1，Pe由1增大到3時，-φ′(0)增大156%；當(dāng)γ=Pe=1，Sc由1增大到3時，-φ′(0)增大15.2%.

表2 其他物理參數(shù)下的流動及傳熱傳質(zhì)特性

4 結(jié)語

由于趨旋微生物對流能避免納米顆粒團聚，解決納米流體不穩(wěn)定問題，可進一步促進高導(dǎo)熱納米流體在工程中的應(yīng)用，因此研究含趨旋微生物的穩(wěn)態(tài)二維納米流體在可拉伸缸表面上的邊界層流動現(xiàn)象，探討各物理參數(shù)對速度場、 溫度場、 濃度場和微生物密度場的影響，得出以下3個主要結(jié)論.

1) 在靠近缸體表面處，納米流體速度、 溫度、 體積分數(shù)和微生物密度均隨著曲率參數(shù)的增大而減小，而在遠離缸體表面時變化情況卻相反；

2)Brown運動參數(shù)和熱泳參數(shù)對納米流體溫度的影響類似，而對體積分數(shù)的影響卻相反；

3) 微生物密度隨著Brown運動參數(shù)、 Peclet數(shù)、 Schmidt 數(shù)和生物對流常數(shù)的增大而減小，而隨著熱泳參數(shù)的增大而增大.

以上結(jié)論可為類似物理模型邊界層控制或?qū)嶒炁渲眠m合不同應(yīng)用場合的穩(wěn)定納米流體提供參考，進一步推進納米流體在工程中的應(yīng)用.