寧利中，張 迪，寧碧波，田偉利

(1.西安理工大學(xué) 水利水電學(xué)院，西安 710048； 2.嘉興學(xué)院 建筑工程學(xué)院， 浙江 嘉興 314001；3.上海大學(xué) 建筑系， 上海 200444)

0 引 言

從底部加熱的兩平板之間的自然對(duì)流，即Rayleigh-Benard對(duì)流，是研究流動(dòng)穩(wěn)定性和對(duì)流時(shí)空特性的模型之一。自1900年開(kāi)始研究Rayleigh-Benard對(duì)流以來(lái)，研究工作取得了很大進(jìn)展[1-2]。隨著計(jì)算機(jī)的出現(xiàn)，數(shù)值模擬已成為主要的研究手段之一。文獻(xiàn)[3-12]混合流體對(duì)流的流體力學(xué)基本方程組的數(shù)值模擬，研究了混合流體對(duì)流中的定常波，行波，局部行波形成，缺陷形成及其時(shí)空結(jié)構(gòu)與非線性動(dòng)力學(xué)特性等。對(duì)傳統(tǒng)的Rayleigh-Benard對(duì)流模型系統(tǒng)施加一個(gè)水平流動(dòng)，或者說(shuō)，對(duì)于水平流動(dòng)問(wèn)題從底部施加一個(gè)熱作用，這樣，形成一個(gè)新的對(duì)流模型系統(tǒng)，即具有水平流動(dòng)的Rayleigh-Benard對(duì)流問(wèn)題。文獻(xiàn)[13-18]通過(guò)具有水平流動(dòng)的單流體Rayleigh-Benard對(duì)流問(wèn)題的流體力學(xué)基本方程組的數(shù)值模擬，研究了單流體對(duì)流中邊界條件與進(jìn)口對(duì)對(duì)流斑圖的影響，局部行波特性及對(duì)流斑圖的分區(qū)等。文獻(xiàn)[19-24]通過(guò)具有水平流動(dòng)的混合流體Rayleigh-Benard對(duì)流問(wèn)題的流體力學(xué)基本方程組的數(shù)值模擬，研究了混合流體對(duì)流中周期變化的局部行波結(jié)構(gòu)，充分發(fā)展的行波結(jié)構(gòu)及水平流動(dòng)對(duì)對(duì)流斑圖的影響。這些工作都是基于流體力學(xué)基本方程組進(jìn)行的。擾動(dòng)方程組一般被使用在流動(dòng)穩(wěn)定性分析中獲得對(duì)流發(fā)生的臨界值?？墒牵墨I(xiàn)[25-29]通過(guò)混合流體Rayleigh-Benard對(duì)流問(wèn)題的流體力學(xué)擾動(dòng)方程組的數(shù)值模擬，獲得了混合流體對(duì)流中的擾動(dòng)的成長(zhǎng)，對(duì)流的行波結(jié)構(gòu)，局部行波結(jié)構(gòu)，對(duì)傳波結(jié)構(gòu)及其時(shí)空結(jié)構(gòu)與非線性動(dòng)力學(xué)特性等。文獻(xiàn)[30-34]在Rayleigh-Benard對(duì)流問(wèn)題研究中獲得了某些進(jìn)展。因此，為了用擾動(dòng)方程組研究具有水平流動(dòng)的Rayleigh-Benard對(duì)流問(wèn)題，有必要對(duì)具有水平流動(dòng)的Rayleigh-Benard對(duì)流問(wèn)題的擾動(dòng)方程組進(jìn)行推導(dǎo)與研究。從而推動(dòng)Rayleigh-Benard對(duì)流問(wèn)題的深入研究。

本文討論了混合流體對(duì)流的流體力學(xué)方程組和擾動(dòng)方程組?；诓夹聊箍私?，給出了考慮Soret效應(yīng)和Dufour效應(yīng)的混合流體對(duì)流的基本方程組。對(duì)于混合液體在忽略Dufour效應(yīng)的情況下，推導(dǎo)了具有水平流動(dòng)的混合流體對(duì)流的擾動(dòng)方程組。應(yīng)用擾動(dòng)方程組可以計(jì)算具有水平流動(dòng)的混合流體對(duì)流的穩(wěn)定性和對(duì)流時(shí)空結(jié)構(gòu)特性。如果忽略擾動(dòng)方程組中二階以上的高階項(xiàng)，方程組可簡(jiǎn)化成線性穩(wěn)定分析的擾動(dòng)方程組。

1 流體力學(xué)基本方程組

考慮在均勻重力場(chǎng)下從下面加熱的流體混合物水平層的對(duì)流問(wèn)題(圖1)。當(dāng)瑞利數(shù)超過(guò)臨界值時(shí)，發(fā)生對(duì)流運(yùn)動(dòng)。假設(shè)系統(tǒng)是一個(gè)忽略了滾動(dòng)軸向變化的二維對(duì)流滾動(dòng)。坐標(biāo)的原點(diǎn)位于腔體中底板與左端壁以直角相交的點(diǎn)上，x，z表示垂直于滾動(dòng)軸線的笛卡爾坐標(biāo)，z軸朝上為正，x軸朝右為正。在布辛涅斯克近似的框架內(nèi)，描述系統(tǒng)的基本方程可以寫(xiě)成：

▽·U=0

(1)

(2)

(3)

(4)

式中▽為梯度算子；▽2為L(zhǎng)aplace算子；U(U,0,W)，T，C，p*，ρ，ν,g，Cp,τ分別為速度矢量，溫度，濃度，壓強(qiáng)，密度，運(yùn)動(dòng)黏性系數(shù)，重力加速度，比熱，時(shí)間，下標(biāo)0表示相應(yīng)物理量的參考值；Jc為物質(zhì)通量；q為熱通量。Jc和q可根據(jù)Onsager的線性非平衡熱力學(xué)理論確定[35-36]，即：

(5)

(6)

圖1 對(duì)流示意圖Fig.1 Convection schemetic diagram

式中κT，λ，D分別為與Soret效應(yīng)有關(guān)的系數(shù)、導(dǎo)熱系數(shù)和濃度擴(kuò)散系數(shù)；μ為混合物的局部化學(xué)勢(shì)。筆者簡(jiǎn)單地描述Soret效應(yīng)和Dufour效應(yīng)。它們發(fā)生在一個(gè)力學(xué)平衡的非等溫流體系統(tǒng)中。對(duì)于具有擴(kuò)散系數(shù)D>0的混合物的流體層，當(dāng)在溶液中幾乎瞬間建立穩(wěn)定的溫度梯度時(shí)，熱擴(kuò)散產(chǎn)生組分的濃度梯度，從而使該組分的擴(kuò)散通量與熱擴(kuò)散方向相反。一段時(shí)間后，熱擴(kuò)散和反擴(kuò)散將完全平衡，整個(gè)系統(tǒng)將達(dá)到一個(gè)穩(wěn)定狀態(tài)。因此，在這個(gè)系統(tǒng)中，濃度梯度是由溫度梯度決定的。這種現(xiàn)象稱為Soret效應(yīng)。它已經(jīng)反映在物質(zhì)通量Jc的表達(dá)式中，即式(5)右側(cè)的第二項(xiàng)。Soret效應(yīng)的逆現(xiàn)象是Dufour效應(yīng)，它是濃度梯度導(dǎo)致的溫度梯度的形成。它已經(jīng)反映在熱通量q的表達(dá)式中，即式(6)右側(cè)的第二項(xiàng)。

考慮溫度T，濃度C離開(kāi)參考值偏差很小，質(zhì)量密度的狀態(tài)方程為:

ρ=ρ0[1-α(T-T0)-β(C-C0)]

(7)

式中恒壓和恒濃度下的熱膨脹系數(shù)α和恒溫和恒壓下的溶質(zhì)膨脹系數(shù)β定義為:

重力場(chǎng)g平行于z軸，方向向下，g=-gez，ez為z軸方向的單位矢量。代物質(zhì)通量Jc，熱通量q和密度ρ到方程(1)～(4)，得：

▽·U=0

(8)

(9)

(10)

(11)

2 無(wú)因次流體力學(xué)基本方程組

對(duì)方程(8)～(11)中的所有物理量用腔體高度d，熱擴(kuò)散率κ，參考密度ρ03個(gè)基本量無(wú)因次化。選用熱擴(kuò)散率κ同樣包含時(shí)間的因素。選用熱擴(kuò)散率κ更好的反映了熱對(duì)流中的熱擴(kuò)散特性?？紤]Soret效應(yīng)和Dufour效應(yīng)的無(wú)因次形式的混合流體對(duì)流基本方程組可表示為:

▽·δU=0

(12)

(13)

(14)

(15)

為了計(jì)算方便，濃度流δζ被定義為：

δζ=δC-ψδT

(16)

無(wú)因次形式的混合流體對(duì)流基本方程組可表示為：

▽·δU=0

(17)

(18)

(19)

(20)

方程(17)～(20)是考慮Soret效應(yīng)和Dufour效應(yīng)的無(wú)因次形式的混合流體對(duì)流基本方程組。Dufour效應(yīng)發(fā)生在混合氣體中，因此，對(duì)于混合氣體，Dufour效應(yīng)項(xiàng)必須保留，但對(duì)于混合液體，Dufour效應(yīng)項(xiàng)可以忽略，Soret效應(yīng)項(xiàng)必須保留。下面討論混合液體的對(duì)流問(wèn)題，忽略Dufour效應(yīng)項(xiàng)后，無(wú)因次形式的混合流體對(duì)流基本方程組可表示為：

▽·δU=0

(21)

(22)

(23)

(24)

方程(21)～(24)是考慮Soret效應(yīng)的無(wú)因次形式的混合流體對(duì)流基本方程組。利用流體力學(xué)基本方程(21)～(23)，文獻(xiàn)[3-12]模擬了混合流體對(duì)流的時(shí)空結(jié)構(gòu)和動(dòng)力學(xué)特性，文獻(xiàn)[19-24]模擬了具有水平流動(dòng)的混合流體對(duì)流的時(shí)空結(jié)構(gòu)和動(dòng)力學(xué)特性。在方程(21)～(24)簡(jiǎn)化成單流體對(duì)流方程組后，文獻(xiàn)[13-18]模擬了具有水平流動(dòng)的單流體對(duì)流的時(shí)空結(jié)構(gòu)和動(dòng)力學(xué)特性。

3 流體力學(xué)基本方程的傳導(dǎo)解

考慮具有水平流動(dòng)的混合流體對(duì)流系統(tǒng)，對(duì)流未發(fā)生時(shí)方程(21)～(24)的傳導(dǎo)解為：

δUcond=4RePrz(1-z)ex

(25)

δTcond=0.5-z

(26)

δCcond=ψ(0.5-z)

(27)

δζcond=0

(28)

(29)

4 流體力學(xué)擾動(dòng)方程組

4.1 具有水平流動(dòng)的混合流體對(duì)流擾動(dòng)方程組

如果對(duì)流發(fā)生，離開(kāi)傳導(dǎo)解的擾動(dòng)量可表示為：

u=δU-δUcond

(30)

θ=δT-δTcond

(31)

ξ=δC-δCcond

(32)

p′=p-pcond

(33)

η=δζ-δζcond

(34)

式中u=(u,0,w)，代擾動(dòng)量定義式(30)～式(34)到方程(21)～(24)，具有水平流動(dòng)的混合流體對(duì)流擾動(dòng)方程組可表示為:

▽·u=0

(35)

(36)

(37)

(38)

(39)

式中δUcond=RePrz(1-z)。擾動(dòng)方程(35)～(39)是非線性的，可以直接用于具有水平流動(dòng)的混合流體對(duì)流的穩(wěn)定性，擾動(dòng)成長(zhǎng)，對(duì)流時(shí)空結(jié)構(gòu)與動(dòng)力學(xué)特性的分析。如果忽略擾動(dòng)方程組中二階以上的高階項(xiàng)，即非線性項(xiàng)，可得具有水平流動(dòng)的混合流體對(duì)流的線性穩(wěn)定分析所依據(jù)的擾動(dòng)方程組。

4.2 混合流體對(duì)流擾動(dòng)方程組

考慮方程(35)～(39)中的水平流動(dòng)δUcond=0，可得混合流體Rayleigh-Benard對(duì)流的擾動(dòng)方程組，即：

▽·u=0

(40)

(41)

(42)

(43)

(44)

方程(40)～(44)與文獻(xiàn)[25-29]在混合流體Rayleigh-Benard對(duì)流時(shí)空結(jié)構(gòu)與動(dòng)力學(xué)特性分析中采用的擾動(dòng)方程組一致。只是無(wú)因次化時(shí)本文采用了d，κ，ρ03個(gè)基本量，文獻(xiàn)[25-29]采用d，ν，ρ03個(gè)基本量。由于基本量選擇的不同，方程的形式稍有不同而已。

4.3 具有水平流動(dòng)的單流體對(duì)流擾動(dòng)方程組

忽略方程(39)，考慮方程(35)～(38)中的分離比ψ=0，擾動(dòng)濃度流η=0，可得具有水平流動(dòng)的單流體Rayleigh-Benard對(duì)流的擾動(dòng)方程組，即：

▽·u=0

(45)

(46)

(47)

(48)

4.4 單流體對(duì)流擾動(dòng)方程組

忽略方程(38)，考慮方程(35)～(38)中的分離比ψ=0，擾動(dòng)濃度流η=0，水平流動(dòng)δUcond=0，可得單流體Rayleigh-Benard對(duì)流的擾動(dòng)方程組，即：

▽·u=0

(49)

(50)

(51)

(52)

4.5 水平流動(dòng)問(wèn)題的擾動(dòng)方程組

忽略方程(38)～(39)，考慮方程(35)～(37)中的分離比ψ=0，擾動(dòng)濃度流η=0，擾動(dòng)溫度θ=0，可得單流體水平流動(dòng)問(wèn)題的擾動(dòng)方程組，即：

▽·u=0

(53)

(54)

(55)

方程(53)～(55)與文獻(xiàn)[38]采用的單流體水平流動(dòng)問(wèn)題的擾動(dòng)方程組一致，只是無(wú)因次化時(shí)采用的基本量不同，方程的形式稍有不同而已。

上述4.2-4.5節(jié)中的擾動(dòng)方程組是非線性的，可以直接用于對(duì)流穩(wěn)定性和對(duì)流特性的分析。如果忽略二階以上的高階項(xiàng)，即非線性項(xiàng)，就可以獲得線性穩(wěn)定分析所依據(jù)的擾動(dòng)方程組。

5 結(jié) 論

本文討論了具有水平流動(dòng)的混合流體對(duì)流的流體力學(xué)方程組和擾動(dòng)方程組?；诓夹聊箍私疲o出了考慮Soret效應(yīng)和Dufour效應(yīng)的混合流體對(duì)流的基本方程組，給出了基本方程組的無(wú)對(duì)流運(yùn)動(dòng)的傳導(dǎo)解。對(duì)于混合液體在忽略Dufour效應(yīng)的情況下，推導(dǎo)了具有水平流動(dòng)的混合流體對(duì)流的擾動(dòng)方程組。根據(jù)不同的簡(jiǎn)化條件，前述擾動(dòng)方程組可以應(yīng)用于具有水平流動(dòng)的單流體對(duì)流問(wèn)題，混合流體對(duì)流問(wèn)題，單流體對(duì)流問(wèn)題，或者單流體水平流動(dòng)問(wèn)題等情況。本文建議的擾動(dòng)方程組可以應(yīng)用于相應(yīng)條件的對(duì)流穩(wěn)定性分析及對(duì)流時(shí)空結(jié)構(gòu)與動(dòng)力學(xué)特性的模擬。