RayleighBenard对流温度信号的混沌特性分析研究.pdfVIP

Rayleigh-Benard 对流温度信号的混沌特性分析 钟纪华 付强 张朝曦 解放军理工大学理学院流体力学研究中心，江苏南京，211101 摘要：本文首先介绍了相空间重构及关联维，最大lyapunov指数及Kolmogorov熵的计算方法， 随后通过对Rayleigh-Benard对流温度信号的数据分析，计算得到了关联维，最大lyapunov指 数及Kolmogorov熵，证实了这组湍流信号并不是无从控制的随机过程，而存在低维奇异吸 引子，并且这些特征指数能很好的表征湍流的结构。 关键词：Rayleigh-Benard对流 混沌 lyapunov指数 Kolmogorov熵 关联维 一、前言 湍流是普遍存在的流动现象，广泛存在于自然界和工程应用领域。对湍流的研究具有很 大的理论意义和实际工程价值，但目前，湍流作为一个整体问题（不管是它的物理机制还是 它的模型计算）仍然存在。湍流是一个典型的复杂系统，湍流流场各物理量是脉动的，它的 这种随机性是否是非线性方程的解在初始条件下的演化，目前还很难回答。但对湍流信号的 时间序列进行分析，分别计算其时间延迟参数，嵌入维数，lyapunov 指数和 Kolmogorov 熵 等非线性系统研究的重要参数，可以帮助更好的理解湍流的混沌特性，从而更好的揭示湍流 过程的动力学特性。 人们已发现许多类型的湍流都具有内在随机性（即混沌特性）。李昕等[1]通过对低层大 气湍流观测数据计算其吸引子维数与最大lyapunov指数检验了大气湍流的混沌特性，李睿劬 等[2]通过对平板边界层湍流的相空间重构及计算其关联维指出流体从层流通向湍流的转捩 过程与混沌动力学本质上有着相同的性质，并进一步指出奇异吸引子的关联维数值在转捩过 程中从1.04演化到1.81表明了两者在定量上的对应关系；陈国等[3]指出随湍流Re数的增大 关联维也随之增大，表明关联维与能量随时间耗散有定量上的联系；金宁德等人[4]的计算 也表明混沌关联维能较好的表征气液两相流流型的变化。 Keng-Hwee Chiam[5]通过数值模拟论述了Rayleigh-Benard对流（以下简称RBC ）的时空 混沌特性，用唯象的观点阐述了对流结构及扩散过程中具有的混沌特性。MUAGNVS INGI EINARSSON[6]讨论了RBC 圆柱边界变化（通过参数κ控制，κ= radius/depth）对对流 混沌特性的影响，进行了数字模拟并且计算了各种情况下的lyapunov指数和关联维。指出 RBC并不是以前所认为的空间弱相关的，它具有与其他湍流所不同的一些特性。那么，RBC 到底具有哪些特性呢？RBC 的混沌特性参数关联维、最大lyapunov指数，Kolmogorov熵与其 对流结构是否具有定量相对应的关系，本文将通过对RBC各个位置点数据进行混沌特性检 验，进行一些有益的探索。 二、理论介绍 由于混沌系统的内在本质的复杂性，构造完整的系统模型是十分困难的，实验的状态下 甚至我们无法知道真实的系统模型，往往只能测得系统的一个或几个分量的时间序列，这时 可以通过相空间重构的方法重构系统的相空间。根据 Tekens 的相空间重构理论，系统的任 一变量的演化是由与之相互作用的其他分量决定的，因此，这些相关分量的信息就隐藏在任 意分量的发展过程中。于是，只考虑一个分量，并将在某些固定时间延迟点上的观测值作为 第 147 页 新维来处理，从而就可以通过嵌入的方法构造出一个与原系统等价的相空间，并且可以在这 个空间中恢复原来的动力系统，并研究其吸引子的性质。已经证明，当嵌入维数 m 和时滞 г的选择恰当时，重构的相空间可以具有与实际的动力系统相同的几何性质和信息性质，具 有真实空间的所有特征。 重构相空间的基本方法为：设时间序列为{x(t)}，t=1,2…N,(N 为序列长度) 。引入一个滞后 时间τ ，将时间序列做顺移，构造一个m 维的状态向量X (n) ，n=1,2…M( M N =−(m −1)τ ) 即 X (n ) {x (n ),x (n =+τ),x (n +2τ), L,x[n +(m −1)τ]} （1） 我们用这些序列构成一个 m 维的空间 R 。 1. 嵌入维数和延迟时间 Grassberger 和Procaccia