2024年4月9日发(作者：裔念寒)
哈特曼边界层的初级稳定性分析
董帅;林殿吉;吕玉坤
【摘 要】导电流体在法向外置磁场的作用下,在贴近壁面处会形成哈特曼边界层.哈
特曼边界层的稳定性研究对电磁冶金过程和热核聚变反应冷却系统等相关设备的设
计和运行都有着十分重要的意义.本文采用非正则模态稳定性分析方法,对两无限大
绝缘平行平板内导电流体流动的稳定性进行了研究.通过在时间上迭代求解扰动变
量的控制方程组和伴随控制方程组,获得了在磁场作用下初级扰动的增长情况及其
空间分布形式,分析了磁场强度对最优扰动增长倍数Gmax、最优展向波数βopt和
最优时刻topt的影响,并考察了上下两个哈特曼边界层之间的相互作用.结果表明,
最优初始扰动的空间分布形式为沿着流场方向的漩涡,关于法向方向对称或者反对
称.当哈特曼数Ha较大时(Ha＞10),对称漩涡和反对称漩涡形式的初始扰动增长倍
数基本相等;上下两个哈特曼边界层可以认为是彼此独立的,不会相互影响,此时最优
扰动增长倍数Gmax与局部雷诺数R的平方成正比,相应的最优展向波数βopt和
最优时刻topt均正比于哈特曼数Ha.当哈特曼数Ha较小时(Ha＜10),反对称漩涡
形式的初始扰动更为不稳定,其增长倍数大于对称漩涡的增长倍数,且上下两个边界
层之间存在着一定的相互作用,并对整个流场的稳定性产生一定的影响.
【期刊名称】《力学学报》
【年(卷),期】2016(048)002
【总页数】9页(P327-335)
【关键词】哈特曼边界层;初级扰动;稳定性分析;非正则模态;数值模拟
【作 者】董帅;林殿吉;吕玉坤
【作者单位】华北电力大学能源动力与机械工程学院,河北保定071003;华北电力
大学能源动力与机械工程学院,河北保定071003;华北电力大学能源动力与机械工
程学院,河北保定071003
【正文语种】中 文
【中图分类】O361.5
2015–05–18收稿，2015–11–27录用,2016–01–05网络版发表.
1)国家自然科学基金青年科学基金(11302076)，河北省自然科学基金
(A2014502047)和中央高校基本科研业务费(2014MS111)资助项目.
2)董帅，讲师，主要研究方向：磁流体力学、流动稳定性分析.E-mail：
********************.cn,*****************
引用格式：董帅，林殿吉，吕玉坤.哈特曼边界层的初级稳定性分析.力学学
报,2016,48(2):327-335
Dong Shuai,Lin Dianji,Lu¨ primary stability analysis of Hartmann
boundary e Journal of Theoretical and Applied
Mechanics,2016,48(2):327-335
冶金行业中双金属复层合金的制备、高分子聚合物的挤压成型、气动挤压塑料片、
吹制玻璃等加工工艺中均会涉及到磁流体边界层内的流动问题[1-3].在这些领域中，
导电流体被置于电磁场环境下，通过施加电磁力可以控制导电流体的运动状态[4-
5].此外，在机翼上合理布置电磁场及相关设备，可以有效改变绕流边界层的结构，
从而达到减震、减阻等效果[6-7].
导电流体在法向外置磁场的作用下，在贴近壁面处会形成速度梯度很大、厚度很薄
的边界层，即哈特曼边界层[8-9].哈特曼边界层的稳定性研究对电磁制动[10]、电
磁搅拌[11]、晶体增长[12]、热核聚变反应冷却系统[13]等相关设备的设计和制造
都有着十分重要的意义.Lundquist[14]最早提出了判断哈特曼边界层是否稳定(层
流还是湍流)的参数——局部雷诺数R
式中，Re为主流雷诺数，Ha为哈特曼数.Murgatroyd[15]对哈特曼边界层的稳定
性进行了实验研究，发现局部临界雷诺数Rc位于150～250的范围内，哈特曼边
界层内流体开始从层流向紊流转捩.Moresco和Alboussi´ere[16]实验得到的临界
雷诺数Rc约为380，但是Lingwood和Alboussi`ere[17]采用正则模态线性稳定
性分析得到的临界雷诺数Rc均明显高于实验结果.最近，Krasnov等[18]采用直接
数值模拟的方法，通过在流场中施加随机背景噪声扰动来观测流场变化，得到了与
实验值相近的结果，临界雷诺数Rc约为350.基于直接数值模拟结果，他们提出了
如下转捩过程假说：(1)流向漩涡和哈特曼层流流动相互作用形成流向条纹结构；
(2)流向条纹结构失稳，进而过渡到湍流.
Gerard-Varet[19]通过非正则模态稳定性的分析方法对单一无限大平板上的哈特
曼边界层进行研究，认为在壁面绝缘情况下，其最优扰动为流向漩涡，最优扰动增
长倍数Gmax和局部雷诺数R呈平方关系，即Gmax≈5.65×10−和
Castets[20]获得的结果与之类似.
Krasnov等[21]随后对大哈特曼数Ha下，导电流体在方形截面管道内的流动进行
了直接数值模拟.当哈特曼数Ha从0变化到400时，流动从湍流逐渐过渡到层流，
哈特曼边界层厚度不断减小，上下两个哈特曼边界层之间的相互作用也逐渐减弱.
目前，对哈特曼边界层稳定性的研究主要针对较大哈特曼数Ha下的磁流体流动，
对小哈特曼数Ha (Ha < 10)下流动的研究还不充分.本文对两块平行的无限大绝缘
平板内导电流体流动的稳定性进行分析，重点考察初级扰动的增长情况，以及哈特
曼数Ha和雷诺数Re的影响.平板流内上下两个哈特曼边界层之间可能存在着一定
的相互作用，从而对整个流场的稳定性产生影响，这在以前的研究中并未引起足够
的重视，本文将就此展开研究和讨论.
1.1物理模型
本研究考察的物理模型为不可压缩导电流体在两块无限大的绝缘平行平板之间的流
动，在垂直流场的方向施加一个均匀恒定的静态磁场.其中，两平板的间距为d，
磁场强度为B，层流流场中心线上的速度为U.平板内的流体在压力作用下流动，
其体积流量保持恒定.
1.2控制方程组
本文是对低磁雷诺数下导电流体展开的研究，即忽略了流场对外界磁场的影响，只
考虑外界磁场对流场的影响[22].其无量纲控制方程组和边界条件为
其中，e = (0,0,1)为单位矢量，x，y，z分别代表了流向，展向和法向方向.u为速
度矢量，p为压力，φ为电势，u,v和w分别代表了流向、展向、法向速度.进行
无量纲化时，特征速度取哈特曼层流流动中心线上的速度U，特征长度L取平板
间距的一半，即L = d/2；磁场强度特征值取B，相应的电势特征值为σUB，时间
和压力的特征值分别为L/U和ρU2.式(2)中出现了两个重要的无量纲参数，哈特曼
数和主流雷诺数
本文研究的基本流是哈特曼层流流动，流场中的任一物理量ψ可以分解为
其中，Ψ(z,t)表示基本流场的物理量，ψp(x,y,z,t)表示扰动变量.将式(3)代入式(2)
中，减去基本流的控制方程组，并忽略关于扰动的高阶项，即可获得关于扰动的线
性方程组.
本研究主要采用非正则模态稳定性分析方法[23-27]，考察扰动的增长模式及其获
得的最大峰值，为了对扰动的增长情况进行分析，将扰动项展开成如下形式
其中，α为流向波数，β为展向波数.为了定量的描述扰动的增长，选取变量E(t)代
表扰动动能随时间的变化情况
其中，上标“*”代表变量的共轭复数.初始扰动的增长或者衰减可由下式表示
式中，T是所要考察的时刻，0为初始时刻，G称为扰动增长倍数.
当哈特曼数Ha、局部雷诺数R、流向波数α、展向波数β和时刻T分别取某一定
值时，扰动的分布形式仍有若干种可能，这些扰动都是控制方程组的解，目的是要
从中找到某一特定的扰动分布，在T时刻获得扰动增长倍数G的最大值.这一特定
的扰动分布被称之为最优扰动.
1.3求解方法
要获得G因子的最大值，数学上可以采用拉格朗日乘子法求解[28-29].即引入扰动
变量的伴随场，通过耦合求解关于扰动的原始控制方程组和相应的关于伴随场变量
的控制方程组，来获取最优扰动的空间分布形式和增长倍数.为了求解T时刻对应
的最优初始扰动ˆu(z,0)，可以采用循环迭代的数值算法[23]，如图1所示.图1中
顶标“～”代表变量的伴随场，顶标“ˆ”代表变量的原始场.首先由初始时刻的扰
动分布ˆu(z,0)，在时间轴上正向迭代求解控制方程组，即可获得在T时刻的扰动
分布ˆu(z,T)；该扰动分布经由适当变换可作为伴随场变量的初始值˜u(z,T)，然后在
时间轴上逆向迭代求解伴随场变量满足的控制方程组，可得˜u(z,0)，并用该扰动分
布去更新初始扰动ˆu(z,0).如此循环迭代下去，直至˜u(z,0) 和ˆu(z,0)满足一定的收
敛标准，最后需要再次求解控制方程组，计算最优扰动ˆu(z,0)在T时刻的扰动动
能，并代入式(6)中求解扰动增长倍数G.
本文采用了高精度的伪谱方法求解无量纲方程组，在法向方向上采用切比雪夫配置
点布置方式.这种网格节点布置方式的特点是中间稀疏，两边密集，从而可以较好
地求解近壁面边界层内的流动.当哈特曼数较小时，磁场效果不明显，哈特曼边界
层厚度较大，采用较少的网格节点(N = 64)可以精确求解控制方程组.当哈特曼数
较大时，哈特曼边界层厚度变薄，需要相应增加节点个数(N = 128)以满足计算精
度要求.
对于一组确定的局部雷诺数R和哈特曼数Ha，某一初始扰动能够在某一流向波数
α、展向波数β和时刻T下，得到所有扰动增长倍数G的最大值，称为最优扰动
增长倍数Gmax，对应的流向波数α、展向波数β和时刻T分别称为最优流向波
数αopt、最优展向波数βopt和最优时刻topt.对剪切流的研究发现，无论是平行
平板流还是圆管中的流动，最优扰动的初始分布形式均为沿着流场方向的漩涡(α=
0)，并且通常为关于法向反对称或者对称的漩涡[2631].所以本文选取初始扰动为
沿着流场方向的对称漩涡或者反对称漩涡两种情况进行分析，并设定流向波数α=
0，只考虑了展向波数β的变化情况.
图2和图3分别表示局部雷诺数R = 300，哈特曼数Ha = 5和15时，不同展向
波数β下、不同扰动增长倍数G和时刻t之间的关系.当哈特曼数Ha = 5时，相
同的展向波数β下，反对称漩涡在相应时刻的扰动增长倍数G均大于对称漩涡的
扰动增长倍数G，即反对称漩涡获得了更大的扰动能量，但反对称漩涡达到最优扰
动增长倍数Gmax的时刻大于对称漩涡.在哈特曼数Ha = 15时，相同的展向波数
β下，初始对称漩涡和反对称漩涡取得了相等的最优扰动增长倍数Gmax，说明在
较大哈特曼数Ha时，两种不同结构的初始扰动其增长情况基本相同.
图4和图5分别表示局部雷诺数R = 300，哈特曼数Ha = 5和15时，最优扰动
增长倍数Gmax和展向波数β随时间的变化情况.当哈特曼数Ha = 5时，反对称
漩涡扰动动能达到最大值的时刻要大于对称漩涡，展向波数β要小于对称漩涡，
这说明磁场对反对称漩涡的抑制效果更加明显.但是必须明确的是，反对称漩涡的
最大扰动动能要大于对称漩涡.汇总结果如表1所示.当哈特曼数Ha = 15时，磁场
对两种漩涡的影响相同，其最优扰动增长倍数Gmax和展向波数β随时间的变化
基本相同，并在相同的时刻和展向波数β下取得最优扰动.两种情况下，随着展向
波数β和时间t的增加，扰动增长倍数G均逐渐增大，直到G因子在某一时刻t
和某一展向波数β下获取峰值，即Gmax，相应的取得最优展向波数βopt和最优
时刻topt.
图6和图7分别显示局部雷诺数R = 300，哈特曼数Ha = 5和15时，最优初级
扰动在初始时刻的空间分布.哈特曼数Ha = 5时，对称漩涡和反对称漩涡扰动的幅
值是不同的；而哈特曼数Ha = 15时，对称漩涡和反对称漩涡的幅值相同，并且
无论流向、展向还是法向速度，在流场中心区域均保持在零值上下，说明最优扰动
主要分布在边界层附近，当哈特曼数越大(Ha = 15)时，这种趋势更加明显.从图中
还可以看出，当哈特曼数Ha较小时，上下两块平板上的边界层之间有着一定程度
的相互作用；而当哈特曼数较大时，两个边界层之间基本上不会互相影响.
图8显示了初始扰动为对称漩涡和反对称漩涡两种情况下，局部雷诺数R =
300,500,700,1000时，最优扰动增长倍数Gmax随哈特曼数Ha的变化规律.在局
部雷诺数R取值不同时，最优扰动增长倍数Gmax随哈特曼数Ha的主要变化趋
势都是相同的，即哈特曼数Ha < 8时，最优扰动增长倍数Gmax随着哈特曼数
Ha的增大而增大，当哈特曼数Ha > 8时，最优扰动增长倍数Gmax基本不变.
在局部雷诺数R相同，哈特曼数Ha < 8情况下，反对称漩涡的最优扰动增长倍数
Gmax均大于对称漩涡.说明这种形式的扰动更为不稳定，而上下两个边界层之间
的影响也是不可以忽略的.在哈特曼数Ha = 5附近，反对称漩涡的最优.最动扰动
增长倍数Gmax出现了先增大后减小的特征.这表明哈特曼数Ha < 5时，磁场作
用较弱，扰动增长主要受主流雷诺数的制约，即Gmax∝Re2；当哈特曼数Ha =
5时，最优扰动增长倍数Gmax达到最大值；当哈特曼数Ha > 5时，磁场强度的
增加，扰动增长主要受局部雷诺数的制约，即Gmax∝R2.
图9所示为不同哈特曼数Ha、不同漩涡结构时，最优扰动增长倍数Gmax和局部
雷诺数R之间的关系，并与Gerard-Varet[19]的结果进行了对比.在哈特曼数Ha
= 1,3,5,9时，对称和反对称漩涡的Gmax-R曲线出现了明显的差别，反对称漩涡
的最优扰动增长倍数Gmax均大于对称漩涡的值.哈特曼数Ha = 10,15,20时，对
称和反对称漩涡的Gmax-R曲线全部重合，即不仅哈特曼数Ha的大小不影响最
优扰动增长倍数Gmax的大小，漩涡结构对最优扰动增长倍数Gmax的影响也非
常有限.本文研究表明，在壁面绝缘情况下，最优扰动增长倍数Gmax和局部雷诺
数R呈平方关系，Gmax≈6.12×10−4R2，该结果与Gerard-Varet[19]的结果
Gmax≈5.65×10−4R2一致.因此，在哈特曼数Ha较大时，哈特曼边界层很薄，
平行平板间的上、下两个边界层可以认为是相互独立的，不会互相影响.但是对于
小哈特曼数，上下两平板间的哈特曼边界层距离较近，会互相影响.
图10显示了不同漩涡结构时，最优展向波数βopt和哈特曼数Ha之间的关系.哈
特曼数Ha≥10时，最优展向波数βopt和哈特曼数Ha呈线性增长关系；哈特曼
数Ha < 10时，最优展向波数βopt和哈特曼数Ha呈指数增长关系.哈特曼数Ha
< 10时，在相同的哈特曼数Ha下，反对称漩涡在更小的最优展向波数βopt下
达到了最优扰动增长倍数Gmax.
图11表明不同局部雷诺数R、不同漩涡结构时，最优时刻topt和哈特曼数Ha之
间的关系.哈特曼数Ha < 10时，最优时刻topt和哈特曼数Ha二者增长关系不明
显；哈特曼数Ha≥10时，最优时刻topt和哈特曼数Ha呈线性增长关系，且局
部雷诺数R越大，斜率越大.哈特曼数Ha < 10时，反对称漩涡对磁场更加敏感，
其最优时刻topt大于对称漩涡的最优时刻topt.哈特曼数Ha≥10时，磁场对两种
漩涡的抑制效果相同.
(1)当哈特曼数Ha较大时，磁场可以有效抑制扰动的增长，无论是对称漩涡还是
反对称漩涡，其效果相同.在哈特曼数Ha较小的时候，对称漩涡和反对称漩涡对
磁场的敏感性是不同的，反对称漩涡需更长的时间取得最优扰动增长倍数，说明磁
场抑制了其增长速度，但是其取得了比对称漩涡更大的扰动增长倍数，即这种形式
的扰动更为不稳定.
(2)哈特曼数Ha > 10时，最优扰动增长倍数Gmax与局部雷诺数R的平方成正比，
与哈特曼数Ha大小无关，对应的最优展向波数βopt和最优时刻topt均正比于
哈特曼数Ha的一次方.
(3)哈特曼数Ha > 10时，上下两边界层之间无相互影响；哈特曼数Ha < 10时，
上下两边界层存在一定的相互作用，并且对整个流场的稳定性产生一定的影响.
【相关文献】
1 Mukhopadhyay S,Layek GC,Samad of MHD boundary layer fl w over a heated
stretching sheet with variable ational Journal of Heat and Mass
Transfer,2005,48(21):4460-4466
2 Sharma R,Bhargava R,Bhargava P.A numerical solution of unsteady MHD convection
heat and mass transfer past a semi-infinit vertical porous moving plate using element free
Galerkin ational Materials Science,2010,48(3):537-543
3 Rao JS,Sankar cal simulation of MHD effects on convective heat transfer
characteristics of fl w of liquid metal in annular Engineering and
Design,2011,86(2):183-191
4 Molokov S,Moreau R,Moffatt ohydrodynamics-Historical Evolution and
Yorker:Springer,2007:231-233
5 Moreau York:Springer,1990:1-4
6 周本谋,范宝春,陈志华等.流体边界层上电磁力的控制效应研究.力学学报,2004,36(4):472-478
(Zhou Benmou,Fan Baochun,Chen Zhihua,et control effects of electromagnetic
force in the boundary e Journal of Theoretical and Applied
Mechanics,2004,36(4):472-478 (in Chinese))
7陈耀慧,范宝春,周本谋等.翼型绕流的电磁力控制.力学学报,2008,40(1):121-127 (Chen
Yaohui,Fan Baochun,Zhou Benmou,et o-magnetic control of hydrofoil
e Journal of Theoretical and Applied Mechanics,2008,40(1):121-127 (in
Chinese))
8 Hartmann -Dynamics I,Theory of the laminar fl w of an electrically conductive liquid
in a homogeneous magnetic field Kongelige Danske Videnskabernes Selskab,Matematisk-
Fysiske Meddelelser,1937,15(6):1-28
9 Hartmann J，Lazarus -dynamics II:Experimental investigation on the fl w of mercury
in a homogeneous magnetic field Kongelige Danske Videnskabernes Selskab,Matematisk-
Fysiske Meddelelser,1937,15(7):1-45
10 Ha MY,Lee HG,Seong cal simulation of threedimensional fl w,heat
transfer,and solidificatio of steel in continuous casting mold with electromagnetic
l of Materials Processing Technology,2003,133(3):322-339
11 Fujisaki -Mold electromagnetic stirring in continuous Transactions on
Industry Applications,2001,37(4):1098-1104
12 Dold P,Benz ng magnetic fields Fluid fl w and crystal growth
ss in Crystal Growth and Characterization of Materials,1999,38(1):39-58
13 Smolentsev S,Moreau R,Buhler L,et thermoflui issues of liquid-metal
blankets:Phenomena and Engineering and Design,2010,85(7):1196-1205
14 Lundquist s in f¨ur Fysik,1952,5(4):297-347
15 Murgatroyd ments on magnetohydrodynamic channel fl ophical
Magazine,1953,44(359):1348-1354
16 Moresco P,Alboussi´ere mental study of the instability of the Hartmann
l of Fluid Mechanics,2004,504(10):167-181
17 Lingwood RJ,Alboussi´ere the stability of the Hartmann s of
Fluids,1999,11(8):2058-2068
18 Krasnov D,Zienicke E,Zikanov O,et cal study of the instability of the Hartmann
l of Fluid Mechanics,2004,504(10):183-211
19 Gerard-Varet icatio of small perturbations in a Hartmann s of
Fluids,2002,14(4):1458-1467
20 Airiau C,Castets the amplificatio of small disturbances in a channel fl w with a
normal magnetic field Physics of Fluids,2004,16(8):2991-3005
21 Krasnov D,Zikanov O,Boeck cal studyof magnetohydrodynamic duct fl w at
high Reynolds and Hartmann l of Fluid Mechanics,2012,704(2):421-446
22 Davidson Introduction to dge:Cambridge
University Press,2001:117-119
23 Dong S,Krasnov D,Boeck ary energy growth and turbulence suppression in
conducting channel fl w with streamwise magnetic field Physics of
Fluids,2012,24(7):074101
24 李春曦,陈朋强,叶学民.含活性剂液滴在倾斜粗糙壁面上的铺展稳定性.物理学
报,2015,64(1):014702 (Li Chunxi,Chen Pengqiang,Ye ity of surfactant-laden
droplet spreading over an inclined heterogeneous Physica
Sinica,2015,64(1):014702 (in Chinese))
25 Luchini ds-number-independent instability of the boundary layer over a fla
surface:optimal l of Fluid Mechanics,2000,404(2):298-309
26 Butler KM,Farrell -dimensional optimal perturbations in viscousshearfl
sof Fluids A:Fluid Dynamics,1992,4(8):1637-1650
27 Andersson P,Berggren M,Henningson l disturbances and bypass transition
in boundary s of Fluids,1999,11(1):134-150
28 Schmid PJ,Henningson ity and Transition in Shear
Yorker:Springer,2001:515-519
29 Reddy SC,Schmid PJ,Baggett JS,et stability of streamwise streaks and transition
thresholds in plane channel fl l of Fluid Mechanics,1998,365:269-303
30 Zikanov the instability of pipe Poiseuille fl s of Fluids,1996,8(11):2923-
2932
31 Gustavsson growth of three-dimensional disturbances in plane Poiseuille fl
l of Fluid Mechanics,1991,224:241-260
2024年4月9日发(作者：裔念寒)
哈特曼边界层的初级稳定性分析
董帅;林殿吉;吕玉坤
【摘 要】导电流体在法向外置磁场的作用下,在贴近壁面处会形成哈特曼边界层.哈
特曼边界层的稳定性研究对电磁冶金过程和热核聚变反应冷却系统等相关设备的设
计和运行都有着十分重要的意义.本文采用非正则模态稳定性分析方法,对两无限大
绝缘平行平板内导电流体流动的稳定性进行了研究.通过在时间上迭代求解扰动变
量的控制方程组和伴随控制方程组,获得了在磁场作用下初级扰动的增长情况及其
空间分布形式,分析了磁场强度对最优扰动增长倍数Gmax、最优展向波数βopt和
最优时刻topt的影响,并考察了上下两个哈特曼边界层之间的相互作用.结果表明,
最优初始扰动的空间分布形式为沿着流场方向的漩涡,关于法向方向对称或者反对
称.当哈特曼数Ha较大时(Ha＞10),对称漩涡和反对称漩涡形式的初始扰动增长倍
数基本相等;上下两个哈特曼边界层可以认为是彼此独立的,不会相互影响,此时最优
扰动增长倍数Gmax与局部雷诺数R的平方成正比,相应的最优展向波数βopt和
最优时刻topt均正比于哈特曼数Ha.当哈特曼数Ha较小时(Ha＜10),反对称漩涡
形式的初始扰动更为不稳定,其增长倍数大于对称漩涡的增长倍数,且上下两个边界
层之间存在着一定的相互作用,并对整个流场的稳定性产生一定的影响.
【期刊名称】《力学学报》
【年(卷),期】2016(048)002
【总页数】9页(P327-335)
【关键词】哈特曼边界层;初级扰动;稳定性分析;非正则模态;数值模拟
【作 者】董帅;林殿吉;吕玉坤
【作者单位】华北电力大学能源动力与机械工程学院,河北保定071003;华北电力
大学能源动力与机械工程学院,河北保定071003;华北电力大学能源动力与机械工
程学院,河北保定071003
【正文语种】中 文
【中图分类】O361.5
2015–05–18收稿，2015–11–27录用,2016–01–05网络版发表.
1)国家自然科学基金青年科学基金(11302076)，河北省自然科学基金
(A2014502047)和中央高校基本科研业务费(2014MS111)资助项目.
2)董帅，讲师，主要研究方向：磁流体力学、流动稳定性分析.E-mail：
********************.cn,*****************
引用格式：董帅，林殿吉，吕玉坤.哈特曼边界层的初级稳定性分析.力学学
报,2016,48(2):327-335
Dong Shuai,Lin Dianji,Lu¨ primary stability analysis of Hartmann
boundary e Journal of Theoretical and Applied
Mechanics,2016,48(2):327-335
冶金行业中双金属复层合金的制备、高分子聚合物的挤压成型、气动挤压塑料片、
吹制玻璃等加工工艺中均会涉及到磁流体边界层内的流动问题[1-3].在这些领域中，
导电流体被置于电磁场环境下，通过施加电磁力可以控制导电流体的运动状态[4-
5].此外，在机翼上合理布置电磁场及相关设备，可以有效改变绕流边界层的结构，
从而达到减震、减阻等效果[6-7].
导电流体在法向外置磁场的作用下，在贴近壁面处会形成速度梯度很大、厚度很薄
的边界层，即哈特曼边界层[8-9].哈特曼边界层的稳定性研究对电磁制动[10]、电
磁搅拌[11]、晶体增长[12]、热核聚变反应冷却系统[13]等相关设备的设计和制造
都有着十分重要的意义.Lundquist[14]最早提出了判断哈特曼边界层是否稳定(层
流还是湍流)的参数——局部雷诺数R
式中，Re为主流雷诺数，Ha为哈特曼数.Murgatroyd[15]对哈特曼边界层的稳定
性进行了实验研究，发现局部临界雷诺数Rc位于150～250的范围内，哈特曼边
界层内流体开始从层流向紊流转捩.Moresco和Alboussi´ere[16]实验得到的临界
雷诺数Rc约为380，但是Lingwood和Alboussi`ere[17]采用正则模态线性稳定
性分析得到的临界雷诺数Rc均明显高于实验结果.最近，Krasnov等[18]采用直接
数值模拟的方法，通过在流场中施加随机背景噪声扰动来观测流场变化，得到了与
实验值相近的结果，临界雷诺数Rc约为350.基于直接数值模拟结果，他们提出了
如下转捩过程假说：(1)流向漩涡和哈特曼层流流动相互作用形成流向条纹结构；
(2)流向条纹结构失稳，进而过渡到湍流.
Gerard-Varet[19]通过非正则模态稳定性的分析方法对单一无限大平板上的哈特
曼边界层进行研究，认为在壁面绝缘情况下，其最优扰动为流向漩涡，最优扰动增
长倍数Gmax和局部雷诺数R呈平方关系，即Gmax≈5.65×10−和
Castets[20]获得的结果与之类似.
Krasnov等[21]随后对大哈特曼数Ha下，导电流体在方形截面管道内的流动进行
了直接数值模拟.当哈特曼数Ha从0变化到400时，流动从湍流逐渐过渡到层流，
哈特曼边界层厚度不断减小，上下两个哈特曼边界层之间的相互作用也逐渐减弱.
目前，对哈特曼边界层稳定性的研究主要针对较大哈特曼数Ha下的磁流体流动，
对小哈特曼数Ha (Ha < 10)下流动的研究还不充分.本文对两块平行的无限大绝缘
平板内导电流体流动的稳定性进行分析，重点考察初级扰动的增长情况，以及哈特
曼数Ha和雷诺数Re的影响.平板流内上下两个哈特曼边界层之间可能存在着一定
的相互作用，从而对整个流场的稳定性产生影响，这在以前的研究中并未引起足够
的重视，本文将就此展开研究和讨论.
1.1物理模型
本研究考察的物理模型为不可压缩导电流体在两块无限大的绝缘平行平板之间的流
动，在垂直流场的方向施加一个均匀恒定的静态磁场.其中，两平板的间距为d，
磁场强度为B，层流流场中心线上的速度为U.平板内的流体在压力作用下流动，
其体积流量保持恒定.
1.2控制方程组
本文是对低磁雷诺数下导电流体展开的研究，即忽略了流场对外界磁场的影响，只
考虑外界磁场对流场的影响[22].其无量纲控制方程组和边界条件为
其中，e = (0,0,1)为单位矢量，x，y，z分别代表了流向，展向和法向方向.u为速
度矢量，p为压力，φ为电势，u,v和w分别代表了流向、展向、法向速度.进行
无量纲化时，特征速度取哈特曼层流流动中心线上的速度U，特征长度L取平板
间距的一半，即L = d/2；磁场强度特征值取B，相应的电势特征值为σUB，时间
和压力的特征值分别为L/U和ρU2.式(2)中出现了两个重要的无量纲参数，哈特曼
数和主流雷诺数
本文研究的基本流是哈特曼层流流动，流场中的任一物理量ψ可以分解为
其中，Ψ(z,t)表示基本流场的物理量，ψp(x,y,z,t)表示扰动变量.将式(3)代入式(2)
中，减去基本流的控制方程组，并忽略关于扰动的高阶项，即可获得关于扰动的线
性方程组.
本研究主要采用非正则模态稳定性分析方法[23-27]，考察扰动的增长模式及其获
得的最大峰值，为了对扰动的增长情况进行分析，将扰动项展开成如下形式
其中，α为流向波数，β为展向波数.为了定量的描述扰动的增长，选取变量E(t)代
表扰动动能随时间的变化情况
其中，上标“*”代表变量的共轭复数.初始扰动的增长或者衰减可由下式表示
式中，T是所要考察的时刻，0为初始时刻，G称为扰动增长倍数.
当哈特曼数Ha、局部雷诺数R、流向波数α、展向波数β和时刻T分别取某一定
值时，扰动的分布形式仍有若干种可能，这些扰动都是控制方程组的解，目的是要
从中找到某一特定的扰动分布，在T时刻获得扰动增长倍数G的最大值.这一特定
的扰动分布被称之为最优扰动.
1.3求解方法
要获得G因子的最大值，数学上可以采用拉格朗日乘子法求解[28-29].即引入扰动
变量的伴随场，通过耦合求解关于扰动的原始控制方程组和相应的关于伴随场变量
的控制方程组，来获取最优扰动的空间分布形式和增长倍数.为了求解T时刻对应
的最优初始扰动ˆu(z,0)，可以采用循环迭代的数值算法[23]，如图1所示.图1中
顶标“～”代表变量的伴随场，顶标“ˆ”代表变量的原始场.首先由初始时刻的扰
动分布ˆu(z,0)，在时间轴上正向迭代求解控制方程组，即可获得在T时刻的扰动
分布ˆu(z,T)；该扰动分布经由适当变换可作为伴随场变量的初始值˜u(z,T)，然后在
时间轴上逆向迭代求解伴随场变量满足的控制方程组，可得˜u(z,0)，并用该扰动分
布去更新初始扰动ˆu(z,0).如此循环迭代下去，直至˜u(z,0) 和ˆu(z,0)满足一定的收
敛标准，最后需要再次求解控制方程组，计算最优扰动ˆu(z,0)在T时刻的扰动动
能，并代入式(6)中求解扰动增长倍数G.
本文采用了高精度的伪谱方法求解无量纲方程组，在法向方向上采用切比雪夫配置
点布置方式.这种网格节点布置方式的特点是中间稀疏，两边密集，从而可以较好
地求解近壁面边界层内的流动.当哈特曼数较小时，磁场效果不明显，哈特曼边界
层厚度较大，采用较少的网格节点(N = 64)可以精确求解控制方程组.当哈特曼数
较大时，哈特曼边界层厚度变薄，需要相应增加节点个数(N = 128)以满足计算精
度要求.
对于一组确定的局部雷诺数R和哈特曼数Ha，某一初始扰动能够在某一流向波数
α、展向波数β和时刻T下，得到所有扰动增长倍数G的最大值，称为最优扰动
增长倍数Gmax，对应的流向波数α、展向波数β和时刻T分别称为最优流向波
数αopt、最优展向波数βopt和最优时刻topt.对剪切流的研究发现，无论是平行
平板流还是圆管中的流动，最优扰动的初始分布形式均为沿着流场方向的漩涡(α=
0)，并且通常为关于法向反对称或者对称的漩涡[2631].所以本文选取初始扰动为
沿着流场方向的对称漩涡或者反对称漩涡两种情况进行分析，并设定流向波数α=
0，只考虑了展向波数β的变化情况.
图2和图3分别表示局部雷诺数R = 300，哈特曼数Ha = 5和15时，不同展向
波数β下、不同扰动增长倍数G和时刻t之间的关系.当哈特曼数Ha = 5时，相
同的展向波数β下，反对称漩涡在相应时刻的扰动增长倍数G均大于对称漩涡的
扰动增长倍数G，即反对称漩涡获得了更大的扰动能量，但反对称漩涡达到最优扰
动增长倍数Gmax的时刻大于对称漩涡.在哈特曼数Ha = 15时，相同的展向波数
β下，初始对称漩涡和反对称漩涡取得了相等的最优扰动增长倍数Gmax，说明在
较大哈特曼数Ha时，两种不同结构的初始扰动其增长情况基本相同.
图4和图5分别表示局部雷诺数R = 300，哈特曼数Ha = 5和15时，最优扰动
增长倍数Gmax和展向波数β随时间的变化情况.当哈特曼数Ha = 5时，反对称
漩涡扰动动能达到最大值的时刻要大于对称漩涡，展向波数β要小于对称漩涡，
这说明磁场对反对称漩涡的抑制效果更加明显.但是必须明确的是，反对称漩涡的
最大扰动动能要大于对称漩涡.汇总结果如表1所示.当哈特曼数Ha = 15时，磁场
对两种漩涡的影响相同，其最优扰动增长倍数Gmax和展向波数β随时间的变化
基本相同，并在相同的时刻和展向波数β下取得最优扰动.两种情况下，随着展向
波数β和时间t的增加，扰动增长倍数G均逐渐增大，直到G因子在某一时刻t
和某一展向波数β下获取峰值，即Gmax，相应的取得最优展向波数βopt和最优
时刻topt.
图6和图7分别显示局部雷诺数R = 300，哈特曼数Ha = 5和15时，最优初级
扰动在初始时刻的空间分布.哈特曼数Ha = 5时，对称漩涡和反对称漩涡扰动的幅
值是不同的；而哈特曼数Ha = 15时，对称漩涡和反对称漩涡的幅值相同，并且
无论流向、展向还是法向速度，在流场中心区域均保持在零值上下，说明最优扰动
主要分布在边界层附近，当哈特曼数越大(Ha = 15)时，这种趋势更加明显.从图中
还可以看出，当哈特曼数Ha较小时，上下两块平板上的边界层之间有着一定程度
的相互作用；而当哈特曼数较大时，两个边界层之间基本上不会互相影响.
图8显示了初始扰动为对称漩涡和反对称漩涡两种情况下，局部雷诺数R =
300,500,700,1000时，最优扰动增长倍数Gmax随哈特曼数Ha的变化规律.在局
部雷诺数R取值不同时，最优扰动增长倍数Gmax随哈特曼数Ha的主要变化趋
势都是相同的，即哈特曼数Ha < 8时，最优扰动增长倍数Gmax随着哈特曼数
Ha的增大而增大，当哈特曼数Ha > 8时，最优扰动增长倍数Gmax基本不变.
在局部雷诺数R相同，哈特曼数Ha < 8情况下，反对称漩涡的最优扰动增长倍数
Gmax均大于对称漩涡.说明这种形式的扰动更为不稳定，而上下两个边界层之间
的影响也是不可以忽略的.在哈特曼数Ha = 5附近，反对称漩涡的最优.最动扰动
增长倍数Gmax出现了先增大后减小的特征.这表明哈特曼数Ha < 5时，磁场作
用较弱，扰动增长主要受主流雷诺数的制约，即Gmax∝Re2；当哈特曼数Ha =
5时，最优扰动增长倍数Gmax达到最大值；当哈特曼数Ha > 5时，磁场强度的
增加，扰动增长主要受局部雷诺数的制约，即Gmax∝R2.
图9所示为不同哈特曼数Ha、不同漩涡结构时，最优扰动增长倍数Gmax和局部
雷诺数R之间的关系，并与Gerard-Varet[19]的结果进行了对比.在哈特曼数Ha
= 1,3,5,9时，对称和反对称漩涡的Gmax-R曲线出现了明显的差别，反对称漩涡
的最优扰动增长倍数Gmax均大于对称漩涡的值.哈特曼数Ha = 10,15,20时，对
称和反对称漩涡的Gmax-R曲线全部重合，即不仅哈特曼数Ha的大小不影响最
优扰动增长倍数Gmax的大小，漩涡结构对最优扰动增长倍数Gmax的影响也非
常有限.本文研究表明，在壁面绝缘情况下，最优扰动增长倍数Gmax和局部雷诺
数R呈平方关系，Gmax≈6.12×10−4R2，该结果与Gerard-Varet[19]的结果
Gmax≈5.65×10−4R2一致.因此，在哈特曼数Ha较大时，哈特曼边界层很薄，
平行平板间的上、下两个边界层可以认为是相互独立的，不会互相影响.但是对于
小哈特曼数，上下两平板间的哈特曼边界层距离较近，会互相影响.
图10显示了不同漩涡结构时，最优展向波数βopt和哈特曼数Ha之间的关系.哈
特曼数Ha≥10时，最优展向波数βopt和哈特曼数Ha呈线性增长关系；哈特曼
数Ha < 10时，最优展向波数βopt和哈特曼数Ha呈指数增长关系.哈特曼数Ha
< 10时，在相同的哈特曼数Ha下，反对称漩涡在更小的最优展向波数βopt下
达到了最优扰动增长倍数Gmax.
图11表明不同局部雷诺数R、不同漩涡结构时，最优时刻topt和哈特曼数Ha之
间的关系.哈特曼数Ha < 10时，最优时刻topt和哈特曼数Ha二者增长关系不明
显；哈特曼数Ha≥10时，最优时刻topt和哈特曼数Ha呈线性增长关系，且局
部雷诺数R越大，斜率越大.哈特曼数Ha < 10时，反对称漩涡对磁场更加敏感，
其最优时刻topt大于对称漩涡的最优时刻topt.哈特曼数Ha≥10时，磁场对两种
漩涡的抑制效果相同.
(1)当哈特曼数Ha较大时，磁场可以有效抑制扰动的增长，无论是对称漩涡还是
反对称漩涡，其效果相同.在哈特曼数Ha较小的时候，对称漩涡和反对称漩涡对
磁场的敏感性是不同的，反对称漩涡需更长的时间取得最优扰动增长倍数，说明磁
场抑制了其增长速度，但是其取得了比对称漩涡更大的扰动增长倍数，即这种形式
的扰动更为不稳定.
(2)哈特曼数Ha > 10时，最优扰动增长倍数Gmax与局部雷诺数R的平方成正比，
与哈特曼数Ha大小无关，对应的最优展向波数βopt和最优时刻topt均正比于
哈特曼数Ha的一次方.
(3)哈特曼数Ha > 10时，上下两边界层之间无相互影响；哈特曼数Ha < 10时，
上下两边界层存在一定的相互作用，并且对整个流场的稳定性产生一定的影响.
【相关文献】
1 Mukhopadhyay S,Layek GC,Samad of MHD boundary layer fl w over a heated
stretching sheet with variable ational Journal of Heat and Mass
Transfer,2005,48(21):4460-4466
2 Sharma R,Bhargava R,Bhargava P.A numerical solution of unsteady MHD convection
heat and mass transfer past a semi-infinit vertical porous moving plate using element free
Galerkin ational Materials Science,2010,48(3):537-543
3 Rao JS,Sankar cal simulation of MHD effects on convective heat transfer
characteristics of fl w of liquid metal in annular Engineering and
Design,2011,86(2):183-191
4 Molokov S,Moreau R,Moffatt ohydrodynamics-Historical Evolution and
Yorker:Springer,2007:231-233
5 Moreau York:Springer,1990:1-4
6 周本谋,范宝春,陈志华等.流体边界层上电磁力的控制效应研究.力学学报,2004,36(4):472-478
(Zhou Benmou,Fan Baochun,Chen Zhihua,et control effects of electromagnetic
force in the boundary e Journal of Theoretical and Applied
Mechanics,2004,36(4):472-478 (in Chinese))
7陈耀慧,范宝春,周本谋等.翼型绕流的电磁力控制.力学学报,2008,40(1):121-127 (Chen
Yaohui,Fan Baochun,Zhou Benmou,et o-magnetic control of hydrofoil
e Journal of Theoretical and Applied Mechanics,2008,40(1):121-127 (in
Chinese))
8 Hartmann -Dynamics I,Theory of the laminar fl w of an electrically conductive liquid
in a homogeneous magnetic field Kongelige Danske Videnskabernes Selskab,Matematisk-
Fysiske Meddelelser,1937,15(6):1-28
9 Hartmann J，Lazarus -dynamics II:Experimental investigation on the fl w of mercury
in a homogeneous magnetic field Kongelige Danske Videnskabernes Selskab,Matematisk-
Fysiske Meddelelser,1937,15(7):1-45
10 Ha MY,Lee HG,Seong cal simulation of threedimensional fl w,heat
transfer,and solidificatio of steel in continuous casting mold with electromagnetic
l of Materials Processing Technology,2003,133(3):322-339
11 Fujisaki -Mold electromagnetic stirring in continuous Transactions on
Industry Applications,2001,37(4):1098-1104
12 Dold P,Benz ng magnetic fields Fluid fl w and crystal growth
ss in Crystal Growth and Characterization of Materials,1999,38(1):39-58
13 Smolentsev S,Moreau R,Buhler L,et thermoflui issues of liquid-metal
blankets:Phenomena and Engineering and Design,2010,85(7):1196-1205
14 Lundquist s in f¨ur Fysik,1952,5(4):297-347
15 Murgatroyd ments on magnetohydrodynamic channel fl ophical
Magazine,1953,44(359):1348-1354
16 Moresco P,Alboussi´ere mental study of the instability of the Hartmann
l of Fluid Mechanics,2004,504(10):167-181
17 Lingwood RJ,Alboussi´ere the stability of the Hartmann s of
Fluids,1999,11(8):2058-2068
18 Krasnov D,Zienicke E,Zikanov O,et cal study of the instability of the Hartmann
l of Fluid Mechanics,2004,504(10):183-211
19 Gerard-Varet icatio of small perturbations in a Hartmann s of
Fluids,2002,14(4):1458-1467
20 Airiau C,Castets the amplificatio of small disturbances in a channel fl w with a
normal magnetic field Physics of Fluids,2004,16(8):2991-3005
21 Krasnov D,Zikanov O,Boeck cal studyof magnetohydrodynamic duct fl w at
high Reynolds and Hartmann l of Fluid Mechanics,2012,704(2):421-446
22 Davidson Introduction to dge:Cambridge
University Press,2001:117-119
23 Dong S,Krasnov D,Boeck ary energy growth and turbulence suppression in
conducting channel fl w with streamwise magnetic field Physics of
Fluids,2012,24(7):074101
24 李春曦,陈朋强,叶学民.含活性剂液滴在倾斜粗糙壁面上的铺展稳定性.物理学
报,2015,64(1):014702 (Li Chunxi,Chen Pengqiang,Ye ity of surfactant-laden
droplet spreading over an inclined heterogeneous Physica
Sinica,2015,64(1):014702 (in Chinese))
25 Luchini ds-number-independent instability of the boundary layer over a fla
surface:optimal l of Fluid Mechanics,2000,404(2):298-309
26 Butler KM,Farrell -dimensional optimal perturbations in viscousshearfl
sof Fluids A:Fluid Dynamics,1992,4(8):1637-1650
27 Andersson P,Berggren M,Henningson l disturbances and bypass transition
in boundary s of Fluids,1999,11(1):134-150
28 Schmid PJ,Henningson ity and Transition in Shear
Yorker:Springer,2001:515-519
29 Reddy SC,Schmid PJ,Baggett JS,et stability of streamwise streaks and transition
thresholds in plane channel fl l of Fluid Mechanics,1998,365:269-303
30 Zikanov the instability of pipe Poiseuille fl s of Fluids,1996,8(11):2923-
2932
31 Gustavsson growth of three-dimensional disturbances in plane Poiseuille fl
l of Fluid Mechanics,1991,224:241-260