格子玻尔兹曼方法的研究激发介质中的非线性波.pdf

格子玻尔兹曼方法研究激发介质中的非线性波 研究生：陈 超 导 师：邓敏艺副教授 孔令江教授 专 业：理论物理 研究方向：计算物理 年级：2007 级 摘 要 激发介质中非线性波的性质是当前研究热点之一。在很多实际的系统中螺旋波及其 破碎是有害的，因此对非线性波的性质研究有着重大的现实意义。随着格子玻尔兹曼方法 的发展和成熟，并凭借其精度较高、程序代码简单、稳定性好等优点，成为了研究非线性 波性质的有效工具。 本文采用格子玻尔兹曼方法，对 Selkov 反应扩散系统中的非线性波进行研究。首先 建立了反应扩散方程的九速正方格子模型，由该模型的格子玻尔兹曼方程推导得到反应扩 散方程。接着进行计算机数值模拟，模拟时将二维系统划分为 300×300 个格子，采用无 流边界条件。主要的数值模拟工作由以下三部分组成: 一、 选取分岔参数a 0.76 ，b 0.02 ，χ 0.1 ，κ 5.0 ,DX DY 0.1，数值模拟结 果显示系统具有可激发性；在相同参数下，从不同的初始态出发，系统可演化到达不同的 非线性波；在不同的参数下，从相同的初始态出发，系统可演化到不同的三种状态：螺旋 波、时空混沌态和均匀定态。发现Selkov反应扩散系统中螺旋波的失稳是Doppler失稳。 对多种参数组合下的系统演化行为进行了模拟，得到了参数在一定范围内反映系统状态的 相图。 二、 根 据 格 子 玻 尔 兹 曼 方 法 理 论 定 义 了 系 统 内 能 函 数 ， 其 表 达 式 为 1 300 300 8 2 E ∑∑∑ek f k 。数值模拟结果表明：在均匀态下激发介质中的内能随分岔参数a 2 i 1 j 1 k 0 增大而线性增加，在螺旋波态下激发介质中的内能随参数a 的增大反而以指数形式递减。 在相同的系统参数下，系统的状态分别为行波、靶波、螺旋波时，系统内能随时间作小幅 的周期变化，且变化幅度不相同，经处理后发现不同波态的平均内能相差不大，这是因为 在参数相同时（即系统提供能量大致不变）由不同的初始条件得到不同的稳定波态，它们 的内能应该比较接近。通过计算机数值计算螺旋波失稳前后系统的内能还发现：螺旋波失 稳前后系统的内能急剧降低，这说明螺旋波失稳原因在于系统提供的能量不足。 三、 根据九速正方格子模型，定义相应的格子玻尔兹曼熵函数 H ，定义式为 300 300 2 4 3 4 4 3 H ( f (ln f −1) +ln( ) f +ln(6) f +f (ln f −1) +ln( )f ) 。对相 ∑∑∑∑ lk lk ∑ 1k ∑ 2k 0 0 0 i 1 j 1 l 1 k 1 2 k 1 k 1 8 I 同参数不同波态的系统以及螺旋波失稳前后的熵进行计算机数值模拟，结果显示：相同参 数下各稳定波态的熵随时间作小幅的周期变化，且变化幅度不相同，经处理发现不同稳定 波态的熵并不相同，说明各稳定波态的有序度不相同。螺旋波失稳后系统的熵突然增大， 这结果说明系统处于自发的相变过程时，系统总是向着熵增大的方向发展。 由于非线性系统的复杂性，激发介质中非线性波的性质还有许多问题有待深入研究， 本文结束时在总结全文工作的基础上对此作了分析和展望。 关键词：可激发介质，格子玻尔兹曼方法，螺旋波，内能，熵 II Lattice Boltzmann Method For Nonlinear Waves In