旋转的Rayleigh-Bénard问题Lorenz模型的动力学行为及数值仿真

王贺元, 陈相霆

(沈阳师范大学 数学与系统科学学院, 沈阳 110034)

1963年,Lorenz[1]将截谱方法作用于Rayleigh-Bénard,在流模型上推导出了三维Lorenz方程。后来,国内外的众多学者研究了Lorenz模型的高维数推广、Lorenz吸引子的存在性及性态[2-3]、Couette-Taylor流的力学机理等[4-10]。Bhattacharjee等[11]和张银[12]先后对旋转Rayleigh-Bénard问题的四维Lorenz模型进行了讨论。本文对他们的模型进行了进一步讨论,进行了全局稳定性分析并对系统的动力学行为进行了数值仿真。

1 数学模型

旋转的Rayleigh-Bénard问题无量纲化后的扰动方程可通过如下的偏微分方程来描述[13]:

(1)

边界条件为自由边界∂zu|z=0,1=∂zv|z=0,1=w|z=0,1=θ|z=0,1=0。式(1)中:u=(u,v,w)表示速度的扰动;θ表示温度的扰动;p表示压强;R是Rayleigh数;T是Taylor数,物理意义是旋转的大小;k=(0,0,1);Pr是Prandtl数。

对方程组(1)化简后展开,选取模式后代入式(1)得到方程组(2)[13]:

(2)

方程组(2)中的X与对流的强度成比例,Y与上下层流体间的温度差成比例,Z与垂直温差的非线性强度成比例,G与流体涡旋强度成比例。

2 全局稳定性分析

2.1 全局吸引子的存在性

对于方程(2)作代换xx,yy,zz+r+Pr后作运算取u(t)=(X(t),G(t),Y(t),Z(t)),H=R4,由上述变量表示可知,参数r,Pr,b均为正数。

2.2 全局稳定性分析

下面利用李雅普诺夫函数法对系统(2)进行全局稳定性分析。构造李雅普诺夫函数

V(x,g,y,z)=13x2+13g2+5y2+5(z-63)2=K>0

很明显,当K是常数时,上式表示一个四维椭球面,把这个椭球面所包围的单连通区域记做E,K越大,椭球面越大。在方程组(2)的混沌区(Pr=5,b=8/3,r=50,τ=1)求V的导数:

显然,下式可以视为一个四维椭球面,记为U:

由李雅普诺夫定理[14]的分析可知,E外的轨线都会进入E内。由此可知,E是这个Rayleigh-Bénard系统的捕捉区。虽然该系统的平衡点都不稳定,但是系统依然具有全局稳定性。