特殊边界条件下的二维反应扩散体系的研究.pdfVIP

第37卷 增刊4 稀有金属材料与工程 V01．37，Suppl．4 2008生 10月 RAREMETALMATERIAI，SANDENGINEERING October2008 特殊边界条件下的二维反应扩散体系研究 蒋 密，马 平，聂瑞娟，王福仁，戴远东 (北京大学人工微结构和介观物理国家重点实验室，北京100871) 摘要：反应扩散作为一种普遍的现象，可以描述大量实际问题。由于边界条件在各种问题中的重要性，在模拟过程 中一般采用零流或周期性边界条件。本文讨论了边界也为波源情况下螺旋波的演变过程。利用Field．Koros．Noyes模型 对于二维情况下不同大小介质的数值模拟表明，随着参数的变化，介质中的波有不同的表现，并会伴随失稳现象。 关键词：螺旋波；Field-Koros-Noyes模型：边界条件 中图法分类号：051l+．9 文献标识码：A 文章编号：1002-185X(2008)S4—064．03 况’研究两对螺旋波在此介质中的演变。 1 引 言 2原理 反应扩散是自然界中的普遍现象，在物理学、化 学、生物学的广阔领域里都出现了各种各样的反应扩 鉴于本文是在一般情况下研究螺旋波的行为，因 散体系。其中可激发体系中的反应扩散过程【¨】可以表 此本文数值模拟采用的模型是具有代表性的 现出更为丰富的现象。从物理化学系统到生物组织， Field．Koros．Noyes模型 从一维到三维介质，实验和理论计算以及数值模拟都 (1) 譬：一1[“l—v)一(2qa_生+∥)型】+V2v 表明介质在不同的条件下可以表现出各种斑图的形 研 F l—W V十∥ 成。一种典型的可激介质中的斑图是螺旋波【10】。 望竺：’，一口兰+v2w (2) 研 1一w 自从20世纪70年代，Winfree首先在化学 其中，v是激发变量(快变量)，w是恢复变最(慢变 Belousov．Zhabotinsky(Bz)系统中观察到了螺旋波现 象，大量的有关螺旋波以及它的三维衍生物卷轴波的 量)。8一般是远小于1的数，在这种情况下，系统一 理论和实验研究随即展开【6，‘71。至今。螺旋波和卷轴波 般将成为可激发系统。由于可激发系统中两个变量变 的相关问题仍是科学界的热点，其中一个重要的原因 化的时间尺度有较大的差别，因此系统的演化将表现 是20世纪90年代初期在实验上观察到了之前设想了 出一些特殊的区别于其他反应扩散体系的性质，而螺 多年的心外膜表面的螺旋波现象，从而引起了心脏病 旋波是可激发系统中出现的典型斑图。g，a，∥，∥是系 学家的关注。由于长久以来心脏的纤维性颤动这种致 统行为的参数。 命的心脏病一直困扰着人类，针对它的研究层出不穷。 3方法 对它的动力学演变过程的研究是一个重要的方面，目 前对纤维性颤动的认识是：它是由心脏中的螺旋波失 本文采用欧拉法进行数值模拟，时间步长△ 稳引起的严重疾病【5】。虽然大量的生理学实验也正集 中于纤维性颤动的认识，但鉴于实验条件的限制和某 些实验方法的不足，至今仍未在实验中观察到心脏中 胪0．000227，g=0．00051。初始条件为 的三维卷轴波。因此在尽前实验技术相对落后于