用元胞自动机模型与研究二维激发介质中的非线性波.pdf

用元胞自动机模型研究二维激发介质中的非线性波 研究生：张立升 导 师：刘慕仁教授、孔令江教授 专 业：理论物理 方 向：计算物理 年级：2006 级 摘 要 激发介质广泛存在于物理系统、化学系统和生物系统中。人们对激发介质的研究已 经持续了多年，其应用领域也涉及到医学，表面科学等。在激发介质中存在的各种波形 中，以螺旋波的结构最为复杂，波形也最为丰富(如反向螺旋波、多臂螺旋波等) ，人们 对螺旋波的研究也最为深入。时至今日，激发介质中的非线性波仍然是非线性科学界研 究的热点之一。 本文以Greenberg and Hasting提出的元胞自动机激发介质模型为基础，对激发介 质中的行波、靶波和螺旋波进了具体的研究工作。主要工作如下： 一、数值模拟并分析了行波的波速u 与邻域半径R 和激发阈值K 的关系。结果表明， 对于单激发态下的行波，其传播速度随邻域半径增加而增加，随激发阈值增加而减小， 三者的具体的关系满足u=R +1-é K ù 。同时，我们分析得出了行波的激发条件，即 ê2R +1 ú ê ú é K ù é K ù。考虑到元胞的激发周期T 就是元胞的状态数n ，本文分析获得行波波 R+1- ³ ê2R+1ú ê2R+1ú ê ú ê ú 长与波速之间的关系 ， l = nu 。对于多激发态下的行波 ，其 激发条件 就是 æ K ö K é ù é ù。多激发态下的行波激发条件在E = 1 时就变为单激发态下的行波激 EçR+1- ÷³ ç ÷ ê2R+1ú ê2R+1ú ê ú ê ú è ø 发条件。在元胞激发态个数足够多的情况下( E ³ R ) ，多激发态下的行波激发条件就变 é K ù 为 £R 。显然，多激发态下的激发条件更宽松一些，这是因为元胞的激发态个数 ê2R +1 ú ê ú 增加导致系统的激发性变强；元胞的激发态个数的增加提升了系统的激发能力。 二、数值模拟和讨论了单激发态下靶波的性质。本文先模拟了靶波波速随不应态数、 激发阈值及邻域半径的变化。结果表明，靶波的不应态个数不影响波速；激发阈值变大 时，波速阶跃式减小。接着我们模拟和分析了靶波波长与波速的关系，所得到的关系与 单激发态下的行波波长与波速的关系相同。对于靶波的存在条件，我们分析得出结果： u 1 。 R 2 当此条件不满足时，靶波就先从曲率较小的地方破碎。计算机数值模拟验证了此结果。 三、对单激发态下螺旋波波头轨迹进行了计算机数值模拟，分别探讨了不应态数、 激发阈值对波头轨迹的影响。不应态个数越大，螺旋波的波头运动轨迹就越长；系统的 激发阈值越大，波头的运动轨迹越大。 四、以Greenberg-Hasting元胞自动机模型为基础，通过参考 Gil Bub等人的改进 - 3 - 的Greenberg-Hasting元胞自动机模型，并以Felipe等人的模型思想为指导，提出了一种 改进的元胞自动机激发介质模型。Felipe等人的模型考虑了激发介质中激发单元对邻域 内其他激发单元的作用，并将激发单元自身的激发能力也考虑了进来。采用本文的改进 模型，我们借助行波模拟和分析了：(一)系统各参数对系统激发性的影响：1、参数D 的 影响。在a = 0.01 ，n = 10 ，E = 2 ，R = 3 ，q = 0.5 ，当D 0.154 时，系统无法维持行 波的运行；当0.154D 0.32时，行波会在行进过程中很不稳定，虽然它不会消失，但行 波很容易断裂，产生激发碎片元，这些激发碎片元在演化过程中相互作用，最终使激发 介质陷入了湍流状态。当D 0.32时，系统的激发性会显