董 帥， 劉立帥， 葉學(xué)民

(華北電力大學(xué) 能源動(dòng)力與機(jī)械工程學(xué)院，河北 保定 071003)

磁場(chǎng)作用下壁面導(dǎo)電率對(duì)導(dǎo)電流體流動(dòng)穩(wěn)定性的影響

董 帥， 劉立帥， 葉學(xué)民

(華北電力大學(xué) 能源動(dòng)力與機(jī)械工程學(xué)院，河北 保定 071003)

導(dǎo)電流體在外置靜態(tài)磁場(chǎng)作用下流動(dòng)時(shí)，其流動(dòng)特性會(huì)發(fā)生改變，進(jìn)而使流體流動(dòng)的穩(wěn)定性更難以預(yù)測(cè)。采用非正則模態(tài)穩(wěn)定性分析方法和數(shù)值模擬技術(shù)，研究了法向磁場(chǎng)作用下壁面導(dǎo)電率不同時(shí)平行平板內(nèi)導(dǎo)電流體的流動(dòng)穩(wěn)定性。通過(guò)建立并求解擾動(dòng)變量的原始和伴隨控制方程組，考察了平板導(dǎo)電和平板絕緣情形下初級(jí)擾動(dòng)的增長(zhǎng)特征及其空間分布形式。結(jié)果表明，平板導(dǎo)電(Ha≥7)和平板絕緣(Ha≥8)情形下最優(yōu)擾動(dòng)增長(zhǎng)倍數(shù)Gmax和局部雷諾數(shù)R的關(guān)系分別為Gmax≈5.34×10-4R2和Gmax≈6.12×10-4R2，相應(yīng)的最優(yōu)展向波數(shù)βopt與哈特曼數(shù)Ha成正比；平板導(dǎo)電時(shí)的擾動(dòng)增長(zhǎng)特征及其空間分布與平板絕緣情形下基本類似；但在平板導(dǎo)電條件下，磁場(chǎng)對(duì)擾動(dòng)的抑制效果更明顯。

穩(wěn)定性；法向磁場(chǎng)；壁面導(dǎo)電率；非正則模態(tài)；數(shù)值模擬

0 引 言

磁流體力學(xué)(Magnetohydrodynamics, 簡(jiǎn)稱MHD)是一門(mén)研究磁場(chǎng)和流場(chǎng)相互作用的學(xué)科，其研究結(jié)果可被廣泛應(yīng)用在材料制備過(guò)程[1]、航空飛行器的設(shè)計(jì)[2]及熱核聚變反應(yīng)堆包層設(shè)計(jì)和制造[3]等方面。導(dǎo)電流體(液態(tài)金屬、電解質(zhì)等)在磁場(chǎng)作用下流動(dòng)時(shí)會(huì)產(chǎn)生感生電流，該感生電流與外置磁場(chǎng)相互作用產(chǎn)生洛倫茲力，從而改變流體質(zhì)點(diǎn)的運(yùn)動(dòng)。在熱核聚變反應(yīng)堆的包層內(nèi)，液態(tài)金屬在強(qiáng)磁場(chǎng)環(huán)境下流動(dòng)，其流動(dòng)穩(wěn)定性會(huì)影響流體的傳熱傳質(zhì)過(guò)程，進(jìn)而影響包層的換熱效率及安全性。因此，研究磁場(chǎng)作用下液態(tài)金屬的流動(dòng)穩(wěn)定性對(duì)熱核聚變反應(yīng)堆的設(shè)計(jì)及運(yùn)行都有重要的理論意義。此外，磁流體流動(dòng)穩(wěn)定性的研究對(duì)晶體增長(zhǎng)[4]、翼型減震[5]、高分子聚合物擠壓成型[6]等方面也具有十分重要的意義。

由于導(dǎo)電流體高溫、不透明等條件的限制，流場(chǎng)結(jié)構(gòu)很難在實(shí)驗(yàn)過(guò)程中被直接觀測(cè)到，目前該領(lǐng)域的大部分工作，均集中于理論分析和數(shù)值模擬，從上個(gè)世紀(jì)30年代起至今，有大量學(xué)者對(duì)此進(jìn)行了研究。當(dāng)導(dǎo)電流體在外置法向靜態(tài)磁場(chǎng)作用下的平行平板內(nèi)流動(dòng)時(shí)，在平板的近壁面處會(huì)形成厚度很薄速度梯度很大的邊界層[7，8]，即哈特曼邊界層。研究者發(fā)現(xiàn)哈特曼邊界層的厚度δ和哈特曼數(shù)Ha有關(guān)，即：δ∝Ha-1。Lundquist[9]最早發(fā)現(xiàn)哈特曼邊界層的穩(wěn)定性和局部雷諾數(shù)R有關(guān)，即：R=Re/Ha，其中，Re為主流雷諾數(shù)。Krasnov等[10]通過(guò)非正則模態(tài)穩(wěn)定性分析方法和數(shù)值模擬在哈特曼邊界層穩(wěn)定性的研究中發(fā)現(xiàn)在法向磁場(chǎng)作用下，湍流轉(zhuǎn)變?yōu)閷恿鞯呐R界雷諾數(shù)Rc=350，這與Moresco和Alboussiere[11]經(jīng)實(shí)驗(yàn)獲得的臨界雷諾數(shù)Rc=380接近。Moreau[12]推導(dǎo)了磁場(chǎng)作用下平行平板內(nèi)的MHD方程，發(fā)現(xiàn)平板導(dǎo)電(σ=∞)情形下的主流速度比平板絕緣(σ=0)情形下小，且與Ha有關(guān)。Gerard-Varet[13]采用非正則模態(tài)穩(wěn)定性理論分析了單一無(wú)限大平板上的哈特曼邊界層的穩(wěn)定性；結(jié)果表明，在平板導(dǎo)電和平板絕緣情形下，其最優(yōu)擾動(dòng)為流向漩渦，最優(yōu)擾動(dòng)增長(zhǎng)倍數(shù)Gmax和局部雷諾數(shù)R的關(guān)系分別為Gmax≈4.93×10-4R2和Gmax≈5.65×10-4R2。Krasnov等[14]對(duì)展向磁場(chǎng)作用下，導(dǎo)電流體在無(wú)限大絕緣平行平板內(nèi)的流動(dòng)穩(wěn)定性開(kāi)展了研究，得到了最優(yōu)初始擾動(dòng)的空間分布及增長(zhǎng)率。Vorobev等[15]采用直接數(shù)值模擬分析了平行磁場(chǎng)作用下自由剪切層內(nèi)導(dǎo)電流體的流動(dòng)穩(wěn)定性問(wèn)題，通過(guò)計(jì)算獲得了擾動(dòng)增長(zhǎng)率及其空間分布。Hussain等[16]通過(guò)切比雪夫配置點(diǎn)譜方法研究了絕緣平行平板內(nèi)導(dǎo)電流體在外置展向磁場(chǎng)作用下流動(dòng)的線性穩(wěn)定性；結(jié)果表明，隨磁場(chǎng)強(qiáng)度增加，不穩(wěn)定區(qū)域變小。

近年來(lái)，國(guó)內(nèi)外學(xué)者對(duì)導(dǎo)電流體在矩形管道中的流動(dòng)特性和穩(wěn)定性做了較多研究，即考慮了側(cè)壁的影響。Priede等[17]研究了壁面導(dǎo)電率σ為有限值(0 <σ< ∞)時(shí)矩形管道內(nèi)液態(tài)金屬流動(dòng)的穩(wěn)定性，獲得了基本流的近似解。劉嬋等[18]研究了法向磁場(chǎng)作用下，液態(tài)金屬在上下壁面完全導(dǎo)電(σ=∞)、左右壁面導(dǎo)電率不同的方管中流動(dòng)的線性穩(wěn)定性，得到了擾動(dòng)的增長(zhǎng)率和擾動(dòng)速度在通道截面上的分布。

綜上所述，壁面導(dǎo)電率對(duì)導(dǎo)電流體的流動(dòng)穩(wěn)定性有較大影響，研究此問(wèn)題對(duì)進(jìn)一步指導(dǎo)磁流體邊界層穩(wěn)定性在工業(yè)生產(chǎn)中的應(yīng)用有重要意義。為此，本文探討法向磁場(chǎng)作用下壁面導(dǎo)電和絕緣時(shí)平行平板內(nèi)導(dǎo)電流體的流動(dòng)穩(wěn)定性，分析兩種情形下擾動(dòng)的增長(zhǎng)特征及其空間分布，重點(diǎn)考察哈特曼數(shù)Ha、雷諾數(shù)Re的影響。

1 數(shù)值方法

1.1 物理模型

研究對(duì)象為不可壓縮導(dǎo)電流體在兩塊無(wú)限大平行平板內(nèi)的流動(dòng)，在流場(chǎng)中施加一個(gè)穩(wěn)恒均勻的靜態(tài)法向磁場(chǎng)。平板間距為d，流體的導(dǎo)電率、密度和運(yùn)動(dòng)粘度分別為σ0、ρ和υ，外置磁場(chǎng)強(qiáng)度為B，層流流場(chǎng)中心線上的速度為U。平板內(nèi)的流體靠壓力來(lái)驅(qū)動(dòng)，體積流量保持定值。

1.2 控制方程組

本文針對(duì)導(dǎo)電流體流動(dòng)處于低磁雷諾數(shù)下開(kāi)展的研究，即只考慮外界磁場(chǎng)對(duì)流場(chǎng)的作用，而忽略流場(chǎng)對(duì)磁場(chǎng)的影響，其相應(yīng)的無(wú)量綱控制方程組和邊界條件為

(1)

式中：p為壓力；φ為電勢(shì)；u為速度矢量，u、υ、w分別為流向x、展向y、法向z上的速度；e=(0,0,1)為單位矢量。對(duì)方程組和邊界條件進(jìn)行無(wú)量綱化時(shí)，取哈特曼層流流動(dòng)中心線上的速度U為特征速度，平板間距的一半L為特征長(zhǎng)度，即L=d/2；磁場(chǎng)強(qiáng)度特征值為B，相應(yīng)的電勢(shì)特征值為σ0UB，壓力和時(shí)間的特征值分別為ρU2和L/U。式(1)中有兩個(gè)重要的無(wú)量綱參數(shù)：表征洛倫茲力和粘性力相對(duì)大小的哈特曼數(shù)Ha：

表征慣性力和粘性力相對(duì)大小的主流雷諾數(shù)Re：

1.3 線性分析

將控制方程組的解分解為基本流和擾動(dòng)，因此，流場(chǎng)中的物理量可表示為

(2)

將式(2)代入式(1)，聯(lián)立基本流的控制方程組，并進(jìn)行線性化處理，可得關(guān)于擾動(dòng)的線性方程組：

(3)

為分析擾動(dòng)的增長(zhǎng)特征及其獲得的最大值，采用非正則模態(tài)穩(wěn)定性分析方法[19-21]。為此，將擾動(dòng)項(xiàng)展開(kāi)成如下形式：

(4)

式中：α和β分別為流向波數(shù)和展向波數(shù)。為定量描述擾動(dòng)的增長(zhǎng)特征，選取變量E(t)表示擾動(dòng)動(dòng)能隨時(shí)間的變化：

(5)

式中：上標(biāo)“*”代表變量的共軛復(fù)數(shù)。初始擾動(dòng)的衰減或增長(zhǎng)倍數(shù)由下式表示：

(6)

式中：T為所要考察的時(shí)刻；0是初始時(shí)刻；G稱為擾動(dòng)增長(zhǎng)倍數(shù)。

對(duì)于一組確定的哈特曼數(shù)Ha和局部雷諾數(shù)R，某一初始擾動(dòng)會(huì)在某一展向波數(shù)β、流向波數(shù)α、和時(shí)刻t下，達(dá)到所有擾動(dòng)增長(zhǎng)倍數(shù)G的最大值，即最優(yōu)擾動(dòng)增長(zhǎng)倍數(shù)Gmax，對(duì)應(yīng)的展向波數(shù)、流向波數(shù)、和時(shí)刻分別稱為最優(yōu)展向波數(shù)βopt、最優(yōu)流向波數(shù)αopt和最優(yōu)時(shí)刻topt。在剪切流研究中發(fā)現(xiàn)，無(wú)論是平行平板流[22]，還是邊界層流動(dòng)[23]，或圓管中的流動(dòng)[24]，初始最優(yōu)擾動(dòng)的類型均為沿著流場(chǎng)方向的漩渦(α=0)，且一般為關(guān)于法向?qū)ΨQ或者反對(duì)稱的漩渦。本文對(duì)擾動(dòng)類型為對(duì)稱漩渦的情況進(jìn)行分析，并設(shè)定流向波數(shù)α=0，只考慮展向波數(shù)β的變化。

1.4 模型驗(yàn)證

本文采用高精度的切比雪夫配置點(diǎn)譜方法，網(wǎng)格節(jié)點(diǎn)布置的特點(diǎn)是兩邊密集、中間稀疏，從而可以更好地求解哈特曼邊界層內(nèi)的流場(chǎng)結(jié)構(gòu)。首先，經(jīng)計(jì)算得到Ha=0下，Re=5000時(shí)平板流中的最優(yōu)擾動(dòng)增長(zhǎng)倍數(shù)Gmax=4 897.18，對(duì)應(yīng)的流向波數(shù)α=0，最優(yōu)展向波數(shù)βopt=2.044，這與Butler 和 Farrell[25]計(jì)算的結(jié)果吻合；并通過(guò)大量計(jì)算，驗(yàn)證了網(wǎng)格無(wú)關(guān)性。網(wǎng)格數(shù)目N不同時(shí)的計(jì)算結(jié)果如表1所示，當(dāng)Ha數(shù)增加時(shí)，需相應(yīng)增加網(wǎng)格節(jié)點(diǎn)個(gè)數(shù)，以提高計(jì)算精度。

表1Ha=0、5、10，Re=5 000，網(wǎng)格數(shù)目N不同時(shí)最優(yōu)擾動(dòng)增長(zhǎng)倍數(shù)Gmax

Tab.1 The optimal amplification of perturbationGmaxas function of differentNatHa=0、5、10 andRe=5 000

NGmaxHa=0Ha=5Ha=10324897.73516.84133.80644897.18516.79133.771284897.17516.77133.762564897.17516.77133.75

2 結(jié)果與討論

圖1給出了壁面導(dǎo)電率和哈特曼數(shù)Ha不同時(shí)的速度分布，u0為層流流場(chǎng)中流向速度分量?？梢园l(fā)現(xiàn)，平板導(dǎo)電和平板絕緣情形下無(wú)量綱速度分布規(guī)律一致，表明壁面導(dǎo)電率對(duì)兩種情形下無(wú)量綱速度分布沒(méi)有影響。當(dāng)哈特曼數(shù)Ha=0，即無(wú)磁場(chǎng)作用下，其流動(dòng)為平板泊肅葉流。隨Ha提高，中心區(qū)域速度分布趨于平緩，哈特曼層中速度梯度增大，這是因?yàn)槭┘哟艌?chǎng)后，流場(chǎng)中心區(qū)域感生電流的方向與磁場(chǎng)方向垂直，產(chǎn)生了阻礙流體流動(dòng)的洛倫茲力。而在近壁面處，感生電流會(huì)形成閉合回路，其電流方向與中心處的相反，從而產(chǎn)生推動(dòng)流體流動(dòng)的洛倫茲力。

圖1 壁面導(dǎo)電率和哈特曼數(shù)Ha不同時(shí)的速度分布Fig.1 The velocity distribution with different wall conductivity and Hartmann numbers

圖2為局部雷諾數(shù)R=300、哈特曼數(shù)Ha=5、展向波數(shù)β和壁面導(dǎo)電率不同時(shí)，擾動(dòng)增長(zhǎng)倍數(shù)G隨t的演化特征。由圖可知，當(dāng)β不變時(shí)，平板導(dǎo)電和平板絕緣情形下G隨t的演化規(guī)律基本一致，均呈單峰曲線變化；但在平板導(dǎo)電條件下，G的峰值小于平板絕緣情形，且向右移動(dòng)，表明平板導(dǎo)電時(shí)，G達(dá)到峰值時(shí)的t大于平板絕緣情形。

圖2 R=300，Ha=5，展向波數(shù)β和壁面導(dǎo)電率不同時(shí)G隨t的演化圖Fig.2 Evolution curves of the amplification of perturbation G as function of time t with different β and wall conductivity at R=300, Ha=5

圖3為R=300、Ha=5，壁面導(dǎo)電率不同時(shí)，擾動(dòng)增長(zhǎng)倍數(shù)G和β隨t的演化過(guò)程。隨t持續(xù)，擾動(dòng)增長(zhǎng)倍數(shù)G逐漸增大，直到G在某一時(shí)刻t和某一展向波數(shù)β下得到最大值，即Gmax，此時(shí)對(duì)應(yīng)最優(yōu)展向波數(shù)βopt和最優(yōu)時(shí)刻topt。與平板導(dǎo)電情形相比可以發(fā)現(xiàn)，平板絕緣時(shí)的Gmax在較小的topt和βopt下達(dá)到峰值，且Gmax的峰值高于平板導(dǎo)電情形，表明平板導(dǎo)電時(shí)，磁場(chǎng)對(duì)擾動(dòng)的抑制效果更明顯。由圖3(b)還可看出，同一t下，平板導(dǎo)電時(shí)的β大于平板絕緣情形。

圖3 R=300，Ha=5，壁面導(dǎo)電率不同時(shí)(a)Gmax隨t的演化圖(b) β隨t的演化圖Fig.3 Evolution curves of the optimal amplification of perturbation Gmax (a) and spanwise wavenumber (b) as function of time t with different wall conductivity at R=300, Ha=5

圖4和圖5為R=500，Ha=5和15，壁面導(dǎo)電率不同時(shí)最優(yōu)初級(jí)擾動(dòng)在初始時(shí)刻的空間分布。由圖可知，當(dāng)平板導(dǎo)電和平板絕緣時(shí)，最優(yōu)初級(jí)擾動(dòng)在初始時(shí)刻的空間分布類似，均為流向漩渦，且隨磁場(chǎng)強(qiáng)度增加，漩渦的尺寸變小，對(duì)應(yīng)的展向波數(shù)增加。流向漩渦和主流相互作用，從而產(chǎn)生流向條紋結(jié)構(gòu)，進(jìn)而導(dǎo)致主流失穩(wěn)，發(fā)生旁路轉(zhuǎn)捩。在強(qiáng)磁場(chǎng)作用下，流向漩渦受到強(qiáng)烈抑制作用，其尺度和作用范圍都減小，相應(yīng)的擾動(dòng)增長(zhǎng)倍數(shù)減小，流場(chǎng)也更為穩(wěn)定。兩種情形的主要區(qū)別體現(xiàn)在哈特曼邊界層附近平板導(dǎo)電時(shí)的各個(gè)擾動(dòng)分量略高于平板絕緣情形。由圖還可看出，兩種情形下法向和展向速度分量要比流向速度分量大一個(gè)數(shù)量級(jí)，表明初始時(shí)刻的擾動(dòng)動(dòng)能主要來(lái)自于法向和展向；當(dāng)Ha=15時(shí)，兩種情形下的流向、展向及法向速度在中心區(qū)域都保持在零附近，表明隨Ha增加，哈特曼邊界層厚度變薄，中心區(qū)域范圍擴(kuò)大，也意味著上、下兩個(gè)哈特曼邊界層的流動(dòng)行為基本上互不影響。

圖4 R=500，Ha=5，壁面導(dǎo)電率不同時(shí)最優(yōu)初級(jí)擾動(dòng)在初始時(shí)刻的空間分布Fig.4 Distribution of optimal primary perturbation at initial time with different wall conductivity at R=500, Ha=5

圖5 R=500，Ha=15，壁面導(dǎo)電率不同時(shí)最優(yōu)初級(jí)擾動(dòng)在初始時(shí)刻的空間分布Fig.5 Distribution of optimal primary perturbation at initial time with different wall conductivity at R=500,Ha=15

圖6給出了壁面導(dǎo)電率不同，局部雷諾數(shù)R=300、500、700、1 000時(shí)，Gmax和Ha的關(guān)系。由圖可知，當(dāng)R不同時(shí)，平板導(dǎo)電和平板絕緣情形下的Gmax隨Ha變化趨勢(shì)相同，即隨Ha提高Gmax增大，并在達(dá)到峰值后趨于平緩。在同一Ha下，隨R增加Gmax增大，且平板導(dǎo)電時(shí)的Gmax小于平板絕緣情形，表明平板導(dǎo)電時(shí)，磁場(chǎng)對(duì)擾動(dòng)的抑制作用較強(qiáng)。

圖6 R和壁面導(dǎo)電率不同時(shí)Gmax和Ha的關(guān)系Fig.6 The optimal amplification of perturbation Gmax as function of Ha with different R and wall conductivity

圖7為壁面導(dǎo)電率和Ha不同時(shí)Gmax和R的關(guān)系。由圖可知，當(dāng)Ha不同時(shí)，平板導(dǎo)電和平板絕緣情形下的Gmax隨R演化規(guī)律相同，即隨R提高Gmax增大。在同一R下，平板導(dǎo)電時(shí)的Gmax小于平板絕緣情形，由圖7(b)可看出，當(dāng)Ha較大時(shí)，Ha不會(huì)對(duì)Gmax產(chǎn)生影響。依據(jù)上述結(jié)果可得，平板導(dǎo)電情形下當(dāng)Ha≥7時(shí)，最優(yōu)擾動(dòng)增長(zhǎng)倍數(shù)Gmax和局部雷諾數(shù)R呈平方關(guān)系，即Gmax≈5.34×10-4R2；平板絕緣情形下當(dāng)Ha≥8時(shí)，最優(yōu)擾動(dòng)增長(zhǎng)倍數(shù)Gmax和局部雷諾數(shù)R的關(guān)系為Gmax≈6.12×10-4R2，這兩種情形下得到的結(jié)果與Gerard-Varet[13]的結(jié)果基本吻合，表明Ha較大時(shí)，上、下兩個(gè)哈特曼邊界層基本上各自獨(dú)立，其流動(dòng)穩(wěn)定性與單一無(wú)限大平板情形類似。

圖7 壁面導(dǎo)電率和Ha不同時(shí)Gmax和R的關(guān)系Fig.7 The optimal amplification of perturbation Gmax as function of R with different Ha and wall conductivity

圖8為壁面導(dǎo)電率不同時(shí)，最優(yōu)展向波數(shù)βopt隨哈特曼數(shù)Ha的變化。如圖所示，平板導(dǎo)電和平板絕緣情形下βopt隨Ha變化的規(guī)律基本一致，隨Ha提高βopt增大，其主要區(qū)別是平板導(dǎo)電時(shí)βopt大于平板絕緣情形。由擬合直線1可看出，平板導(dǎo)電情形下當(dāng)Ha≥7時(shí)βopt和Ha成正比，即βopt=0.845Ha；平板絕緣情形下當(dāng)Ha≥8時(shí)，βopt和Ha的關(guān)系為βopt=0.76Ha。

圖8 壁面導(dǎo)電率不同時(shí)βopt和Ha的關(guān)系Fig.8 The optimal spanwise wavenumber βopt as function of Ha with different wall conductivity

3 結(jié) 論

(1) 平板導(dǎo)電情形下當(dāng)Ha≥7時(shí)，最優(yōu)擾動(dòng)增長(zhǎng)倍數(shù)Gmax和局部雷諾數(shù)R呈平方關(guān)系，即Gmax≈5.34×10-4R2，相應(yīng)的最優(yōu)展向波數(shù)βopt與哈特曼數(shù)Ha的關(guān)系為βopt=0.845Ha；平板絕緣情形下當(dāng)Ha≥8時(shí)，最優(yōu)擾動(dòng)增長(zhǎng)倍數(shù)Gmax和局部雷諾數(shù)R的關(guān)系為Gmax≈6.12×10-4R2，相應(yīng)的最優(yōu)展向波數(shù)βopt與哈特曼數(shù)Ha的關(guān)系為βopt=0.76Ha。

(2)當(dāng)Ha較大時(shí)，上、下兩個(gè)哈特曼邊界層基本上各自獨(dú)立，其流動(dòng)穩(wěn)定性與單一無(wú)限大平板情形類似；哈特曼邊界層附近，平板導(dǎo)電時(shí)的各個(gè)擾動(dòng)分量略高于平板絕緣情形。

(3)平板導(dǎo)電時(shí)的擾動(dòng)增長(zhǎng)趨勢(shì)與平板絕緣情形基本一致；但在平板導(dǎo)電條件下，磁場(chǎng)對(duì)擾動(dòng)的抑制效果更明顯。

(4)平板導(dǎo)電和平板絕緣時(shí)，最優(yōu)初級(jí)擾動(dòng)在初始時(shí)刻的空間分布類似，均為沿流向的漩渦；隨磁場(chǎng)強(qiáng)度增加，漩渦尺寸減小，對(duì)應(yīng)的展向波數(shù)增加；受到磁場(chǎng)的限制作用，流向漩渦位于壁面處的哈特曼邊界層附近。

[1] FUJISAKI K. In-mold electromagnetic stirring in continuous casting[J].IEEE Transactions on Industry Applications,2001,37(4):1098-1104.

[2] 劉宗凱, 周本謀, 劉會(huì)星, 等. 電磁流體表面推進(jìn)機(jī)理與效果分析[J]. 物理學(xué)報(bào), 2011, 60(8): 084701.

[3] MULLER U, BUHLER L. Magnetofluiddynamics in channels and containers [M]. New York: Springer, 2001.

[4] 杜嬌嬌,李國(guó)建, 王強(qiáng), 等. 強(qiáng)磁場(chǎng)下不同晶粒尺寸 Fe 薄膜生長(zhǎng)模式演變及其對(duì)磁性能的影響[J]. 金屬學(xué)報(bào), 2015, 51(7): 799-806.

[5] 陳耀慧, 栗保明, 潘緒超, 等. 電磁力控制翼型繞流分離的增升減阻效率研究[J]. 力學(xué)學(xué)報(bào), 2015, 47(3): 414-421.

[6] RAO J S, SANKAR H. Numerical simulation of MHD effects on convective heat transfer characteristics of flow of liquid metal in annular tube [J]. Fusion Engineering and Design, 2011, 86(2): 183-191.

[7] HARTMANN J. Hg-dynamics i: Theory of the laminar flow of an electrically conductive liquid in a homogeneous magnetic field [J]. Danske Videnslabernes Selskab, Math. Fys. Medd, 1937, 15(6): 1-28.

[8] HARTMANN J, LAZARUS F. Hg-dynamics ii: Experimental investigation on the flow of mercury in a homogeneous magnetic field [J]. Danske Videnslabernes Selskab, Math. Fys. Medd, 1937, 15(7): 1-45.

[9] LUNDQUIST S. Studies in magneto-hydrodynamics [J]. Arkiv for Fysik, 1952, 5(4): 297-347.

[10] KRASNOV D S, ZIENICKE E, ZIKANOV O, et al. Numerical study of the instability of the Hartmann layer [J]. Journal of Fluid Mechanics, 2004, 504: 183-211.

[11] MORESCO P, ALBOUSSIERE T. Experimental study of the instability of the Hartmann layer [J]. Journal of Fluid Mechanics, 2004, 504: 167-181.

[12] MOREAU R. Magnetohydrodynamics [M]. Berlin: Springer Netherlands, 1990.

[13] GERARD V D. Amplification of small perturbations in a Hartmann layer [J]. Physics of Fluids, 2002, 14(4): 1458-1467.

[14] KRASNOV D, ROSSI M, ZIKANOV O, et al. Optimal growth and transition to turbulence in channel flow with spanwise magnetic field [J]. Journal of Fluid Mechanics, 2008, 596: 73-101.

[15] VOROBEV A, ZIKANOV O. Instability and transition to turbulence in a free shear layer affected by a parallel magnetic field [J]. Journal of Fluid Mechanics, 2007, 574: 131-154.

[16] HUSSAIN Z, 劉嬋, 倪明玖, 等. 平行平板間導(dǎo)電流體的磁流體動(dòng)力學(xué)穩(wěn)定性[J]. 中國(guó)科學(xué)院大學(xué)學(xué)報(bào), 2013, 30(3): 304-310.

[17] PRIEDE J, ARLT T, BüHLER L. Linear stability of magnetohydrodynamic flow in a square duct with thin conducting walls [J]. Journal of Fluid Mechanics, 2016, 788: 129-146.

[18] 劉嬋, 張年梅, 倪明玖. 側(cè)壁非對(duì)稱的方管 MHD 流動(dòng)的線性穩(wěn)定性分析[J]. 中國(guó)科學(xué)院大學(xué)學(xué)報(bào), 2016, 33(1): 42-49.

[19] 劉梅, 劉玉梅, 王松嶺. 非平整加熱斜面上液膜流動(dòng)穩(wěn)定性研究[J]. 華北電力大學(xué)學(xué)報(bào)(自然科學(xué)版), 2013, 40(5): 19-23.

[20] 裴建軍. 非平整基底上含不溶性活性劑溶液的流動(dòng)特性[D]. 保定：華北電力大學(xué), 2014.

[21] 葉學(xué)民, 閻維平. 液體薄膜流穩(wěn)定性和破斷特性的研究進(jìn)展[J]. 華北電力大學(xué)學(xué)報(bào)(自然科學(xué)版), 2006, 33(6): 63-67.

[22] REDDY S C, SCHMID P J, BAGGETT J S, et al. On the stability of streamwise streaks and transition thresholds in plane channel flows [J]. Journal of Fluid Mechanics, 1998, 365: 269-303.

[23] ANDERSSON P, BERGGREN M, HENNINGSON D S. Optimal disturbances and bypass transition in boundary layers [J]. Physics of Fluids, 1999, 11(1): 134-150.

[24] ZIKANOV O. On the instability of pipe Poiseuille flow [J]. Physics of Fluids, 1996, 8(11): 2923-2932.

[25] BUTLER K M, FARRELL B F. Three-dimensional optimal perturbations in viscous shear flow [J]. Physics of Fluids A: Fluid Dynamics, 1992, 4(8): 1637-1650.

Effect of Wall Conductivity on Stability of Conducting Fluid Flow Under Magnetic Field

DONG Shuai, LIU Lishuai, YE Xuemin

(School of Energy Power and Mechanical Engineering, North China Electric Power University, Baoding 071003, China)

When conducting fluid is passing through an outside static magnetic field, the characteristics of the flow will be different, what's more, the stability of the flow is more difficult to predict. Non-normal modal stability analysis and numerical simulation are used in this paper to study the stability of conducting fluid within a parallel plate with different wall conductivity in a vertical magnetic field. The growth rate and spatial distributions of primary perturbations are calculated respectively in conducting plate and insulating plate by solving iteratively the direct and adjoint governing equations of disturbance variables. The results show that the relationship between increased times of optimal disturbanceGmaxand local Reynolds numberRare different in conducting plate and insulating plate and the relationships can be expressed as follows: Gmax≈5.34×10-4Rand Gmax≈6.12×10-4R2. The corresponding optimal span-wise wave numberβoptis proportional to the Hartmann numberHain some range. The growth of disturbances and spatial distributions of primary perturbations in the conducting plates is similar to the one in the insulating plates. When the plates are conductive, the suppression effect of the magnetic field on the disturbances is more significant.

stability; vertical magnetic field; wall conductivity; non-normal modal; numerical simulation

10.3969/j.ISSN.1007-2691.2017.01.16

2016-03-26.

國(guó)家自然科學(xué)基金資助項(xiàng)目(11302076);河北省自然科學(xué)基金資助項(xiàng)目(A2014502047);中央高?；究蒲袠I(yè)務(wù)費(fèi)資助項(xiàng)目(2014MS111).

O 361.5

A

1007-2691(2017)01-0103-08

董帥(1982-)，男，講師，主要從事磁流體力學(xué)和流動(dòng)穩(wěn)定性分析的研究工作。