二维不可压缩Navier-Stokes方程的七模类Lorenz方程组的动力学行为及其数值模拟

王贺元

（辽宁工业大学理学院，辽宁锦州121001）

二维不可压缩Navier-Stokes方程的七模类Lorenz方程组的动力学行为及其数值模拟

王贺元

（辽宁工业大学理学院，辽宁锦州121001）

本文研究了平面不可压缩的Navier-Stokes方程一个七模类Lorenz方程组的混沌行为问题.利用模式截断的方法，获得了一个七模类Lorenz方程组，证明了该方程组吸引子的存在性，并对其全局稳定性进行了分析和讨论.基于分岔图、最大李雅普诺夫指数、庞加莱截面、功率谱揭示了系统混沌行为的普适特征，仿真分析了系统动力学行为的演化过程.

Navier-Stokes方程；奇怪吸引子；李雅普诺夫函数

1　引言

流动现象是自然界及人类生产科研活动中最为常见的一种物理现象，流动稳定性是流动现象最为关键的问题.作为流动现象应普遍遵循的Navier-Stokes方程是一种典型的非线性偏微分方程，刻划着流体的运动规律，如大气运动、海洋流动、轴承润滑、透平机械内部流动等，研究它对人们认识和控制湍流至关重要.1963年美国气象学家E.Lorenz在研究大气对流时，首次给出了著名的Lorenz方程［9］.所采用的方法是对Navier-Stokes方程和热传导方程进行傅立叶级数展开，截取级数的前三项，得到三模的Lorenz系统.20世纪后期Valter Franceschini又在此方向上进一步扩展，多次和其他学者合作，将二维正方形区域T2=［0，2π］×［0，2π］上不可压缩的Navier-Stokes方程

（其中u为速度场函数；p为液体之间的压力；f为外力场函数，ν为动力粘性系数）进行傅立叶展开并截取其中的有限项，得出五模和七模或者任意模的非线性微分方程组（见文献［1-4］），讨论当雷诺数变化时方程组解的动力学行为.这种截断后来被扩展到三维空间，1988 年V.Franceschini，Inglese和Tebaldi在Commun.Mech.Phys.上发表了三维空间上的有关Navier-Stokes方程五模截断的文章［7］；1991年Franceschini和Zanasi在三维空间上对此方程傅立叶展开，进行七模截断后得到十四个非线性微分方程组成的方程组，随后又对这个复杂的方程组进行了详细的讨论［3］.国内王贺元等人选取不同的截断模式，并把这方面的研究扩展到磁流体，得到相应类Lorenz方程组并分析了系统的动力学行为［10，11］.Franceschini 在1981年给出的一个七模类Lorenz方程组［1］，讨论了这个七模模型定常解的线性稳定性，并对分歧行为进行了数值模拟.本文对此模型的动力学行为进行深入的分析和探讨，证明了该模型吸引子的存在性，并讨论了其全局稳定性，从而在理论上保证了数值模拟的有效性，并且数值模拟了分歧和混沌吸引子的发生过程.

2　七模类Lorenz方程组

下面对二维区域［0，2π］×［0，2π］上Navier-Stokes方程进行傅立叶展开.即对速度函数u，外力场函数f和流体之间的压力p进行如下傅立叶展开

其中K=（h1，h2）是波向量，K⊥=（h2，-h1），rK=rK（t）为时间t的函数.将（2.1）-（2.3）式代入到方程组（1.1），经过一系列运算得到如下形式的微分方程组［1，2］

其中L为波向量集合，并且满足若K∈L，则-K∈L.文献［1］取

其中

在ν=1时，分别令K为K1，K2，K3，K4，K5，K6，K7，代入到方程组（2.4）经大量计算，利用实条件作代换

为

这里xi=xi（t）（i=1，2，···，7）为谱展开系数.截得了七模非线性微分方程组的形式和Lorenz方程组相似，称其为类Lorenz系统.

3　平衡点及稳定性分析

由于系统（2.6）在平衡点处Jacobi矩阵与时间无关，故李雅普诺夫矩阵的特征指数就是Jacobi矩阵的特征值的实部［5，6］，它是刻画吸引子性质的重要指标，尤其对混沌吸引子更为重要.下面对类Lorenz方程组（2.6）线性化，然后根据各个平衡点的李雅普诺夫矩阵的特征指数的变化来讨论平衡点的稳定性.令

对F（X，Re）关于X求导数得到如下李雅普诺夫矩阵

由F（X，Re）=0求出（2.6）式的平衡点，下面根据Liapunov矩阵的特征指数的变化情况讨论各平衡点的稳定性（具体参见文献［1］）.

其中σ=±1，此时平衡点（1）变得不稳定，两个新平衡点（2）是稳定的.

3）当Re＞R2时，所有平衡点（1），（2）都是不稳定的.

4　吸引子的存在性和全局稳定性分析

耗散动力系统的混沌行为是由于存在着一个复杂的吸引子而引起的［8］，而这个吸引子就是系统的所有轨道当时间趋于无穷时收敛到的集合，可能是一个分形或康托集或康托集和一个区间的乘积.很自然地这个“吸引子”就成为数学上用来描述观察到的不稳定流的对象，它的复杂结构就是导致观察到的混沌现象的原因.因此，研究吸引子的存在性和数值模拟就成为一个重要的问题.下面就来证明系统（2.6）的吸引子存在性.

取H=R7，u（t）=（x1，...，x7），对Navier-Stokes方程的七模类Lorenz方程组（2.6）作如下运算

得

因此有

因此有

故有

非线性系统具有全局稳定性时，其轨线所收敛的单连通闭区域称为系统的捕捉区.只要能证明捕捉区的存在，不论其中的定常解是否稳定，系统均具有全局稳定性.而研究系统的全局稳定性主要借助于李雅普诺夫函数方法［5，6］.李雅普诺夫函数方法的基本思想是构造一个函数，然后利用它的性质和这个函数沿方程（2.6）的轨线方向的全导数的性质以确定（2.6）式平衡点的稳定性，以确定系统的捕捉区.

对系统（2.6）构造李雅普诺夫函数为令V（x1，x2，x3，x4，x5，x6，x7）=K，很明显，当K是一正常数时，上式表示H上的一球面，记为E.求V的导数，并利用（2.6）式

于是若把K取得充分大，E即可包围C.这样从式（4.1）式可知在C外面＜0，由李雅普诺夫定理［5］的分析得知E外（2.6）式的解轨线都将进入E内.可见E就是类Lorenz系统（2.6）的捕捉区.虽然这时系统平衡点（1），（2）都不稳定，但系统仍具有全局稳定性：系统最终要收缩到捕捉区内，而区内又无收点，因此系统只能在区内不停的振荡.于是轨线最终要在捕捉区内形成一个不变集合，这就是所谓的吸引子.人们称混沌运动这种具有独特性质和结构的吸引子为奇怪吸引子.它是整体稳定性和局部不稳定性一对矛盾的结合体.其具体形式如何呢？下面就来数值模拟系统（2.6）的奇怪吸引子.

5　数值模拟

随着雷诺数的增大，Lorenz方程组（2.6）的稳定性发生了变化，出现了Hopf分岔和混沌等非线性现象.下面就来详细数值模拟系统（2.6）从分岔到混沌的全过程.

1）当Re＜R2=30.2123时，系统（2.6）的新平衡点是稳定的，解轨线为螺旋线（如图1，2）.

图1：Re=15.60

图2：Re=20.04

2）通过数值计算得方程组（2.6）在Re=R2时平衡点p±处的李雅普诺夫矩阵的一对复共轭特征值穿越虚轴，其实部由负变正，因而系统（2.6）发生了Hopf分岔.从不稳定平衡点p±分叉出闭轨线.如图3-6.

图3：Re=30.2122

图4：Re=30.2124

图5：Re=30.2123

图6：Re=30.2123

3）当Re=71.31时，平衡点p±处分叉出闭轨线开始不稳定，分叉处环面，如图7，8.

4）当Re进一步增大时出现了滞后现象（各种吸引子共存），如图9-13.当Re=248.23时系统发生混沌，出现奇怪吸引子，如下图14-16分别给出了不同雷诺数时奇怪吸引子的大体状态.通过数值计算表明系统在高雷诺数下一直是混沌状态，这与文献［1］的结论是一致的.

图7：Re=71.4

图8：Re=72.0

图9：Re=72.5

图10：Re=73.8

图11：Re=202.4

图12：Re=220.54

图13：Re=230.24

图14：Re=248.23

图15：Re=249.44

图16：Re=255.64

图17：分叉图

图18：最大李雅普诺夫指数

5）图17，18分别给出了系统分岔图和最大李雅普诺夫指数，从分岔图17表明：当Re＜71.31时，系统是稳定的，当Re=71.31时，系统开始不稳定，分叉处环面，之后系统出现滞后现象，当Re=248.23时系统发生混沌，出现奇怪吸引子，算到Re=1000系统始终是混沌状态，在高雷诺数下系统处于湍流状态，这一点也与lorenz系统有明显的区别.图18中给出的最大李雅普诺夫指数与分岔图17是相符的.

6）图19给出了系统的庞加莱截面（Re=252.41），图20给出了系统的功率谱（Re= 255.24），它们均表明了系统的混沌运动特征.

图19：Re=252.41

图20：Re=255.24

6　结论

本文首先给出了七模类Lorenz方程组的推导过程，对此方程组线性化稳定性分析进行了简单介绍.然后证明了此方程组全局吸引子的存在性，并对其全局稳定性进行了分析和讨论，最后数值模拟了雷诺数变化时系统经由不变环面的失稳到达混沌的过程，运用分岔图、最大李雅普诺夫指数、庞加莱截面和功率谱揭示了系统混沌行为的普适特征.

［1］Valter Franceschini，Claudio Tebaldi.A seven-modes truncation of the plane incompressible Navier-Stokes equations［J］.J.Stat.Phys.，1981，25（3）：397-417.

［2］Carlo Boldrighini，Valter Franceschini.A five-dimensional truncation of the plane incompressible Navier-Stokes equations［J］.Commun.Math.Phys.，1979，64：159-170.

［3］Franceschini V，Zanasi R.Three-dimensional Navier-Stokes equations trancated on a torus［J］.Nonl.，1992，4：189-209.

［4］Valter Franceschini，Claudio Tebaldi.Breaking and disappearance of tori［J］.Commun.Math.Phys.，1984，94：317-329.

［5］刘秉正，彭建华.非线性动力学［M］.北京：高等教育出版社，2004，406-415.

［6］谢应齐.非线性动力学数学方法［M］.北京：气象出版社，2001，9-17.

［7］Franceschini V，Inglese G，Tebaldi C.A five-mode truncation of the Navier-Stokes equations on a three-dimensional torus［J］.Commun.Mech.Phys，1988，64：35-40.

［8］Chorin A，Marsden J，Smtle S.Mubalence seminar lecture notes in mathematics［M］.Berlin：Springer-Verlag，1988.

［9］Hilborn R C.Chaos and nonlinear dynamics［M］.Oxford：Oxford Univ.Press，1994．

［10］王贺元，姜悦岭，平面不可压缩Navier-Stokes方程新五模类Lorenz方程组的混沌行为［J］.数学杂志，2010，30（2）：269-272.

［11］高焱，磁流体动力学截断方程组的动力学行为研究［J］.数学杂志，2013，33（4）：671-678.

［12］李开泰，马逸尘.数理方程HILBERT空间方法（下）［M］.西安：西安交通大学出版社，1992.

THE DYNAMICAL BEHAVIOR AND THE NUMERICAL SIMULATION OF THE SEVEN-MODE TRUNCATION SYSTEM OF THE PLANE INCOMPRESSIBLE NAVIER-STOKES EQUATIONS

WANG He-yuan
（School of Sciences，Liaoning University of Technology，Jinzhou 121001，China）

The chaotic behavior of seven-mode Lorenz-like system for the plane incompressible Navier-Stokes equations is studied.By mode truncation，a seven-mode Lorenz equations is obtained.The existence of the attractor of the equations is proved，and the global stability of the equations is discussed.Based on numerical simulation results of bifurcation diagram，Lyapunov exponent spectrum，Poincare section and power spectrum of the system，general features of the system are revealed.The whole process，which shows a chaos behavior with the changing of Reynolds number，is simulated numerically.

Navier-Stokes equations；strange attractor；Liapunov function

MR（2010）主题分类号：65P20；65P40O175.14；O241.81

A

0255-7797（2016）05-1067-10

2014-07-13接收日期：2014-12-18

辽宁省教育厅科研基金（L2013248）；锦州市科学技术基金（13A1D32）资助；国家自然科学基金（11572146）.

王贺元（1963-），男，辽宁黑山，教授，主要研究方向：非线性系统的动力学行为及仿真.

2010 MR Subject Classification：65P20；65P40