扩散过程的元胞自动机模拟

第23卷第1期2011年3月河南工程学院学报(自然科学版)J OURNA L O F HENAN I N ST ITUTE OF ENG I N EER I NG V o l 23,N o 1M ar .2011扩散过程的元胞自动机模拟李延升1,2,侯珂珂1,张保林2(1.许昌学院化学化工学院,河南许昌461000; 2.郑州大学化工与能源学院,河南郑州450007)摘 要:以元胞自动机为研究方法,通过编写M atlab 程序,以演化示意图的形式,模拟了圆域内的扩散过程,给出了累积演化率曲线.同时,用微积分的方法,解析了该扩散过程,并结合实例与元胞自动机方法对比.研究证明,元胞自动机可以生动形象地模拟扩散过程,且与微积分得到的结论一致.关键词:扩散;元胞自动机;扩散方程;模拟中图分类号:TQ019 文献标识码:A 文章编号:1674-330X (2011)01-0001-04收稿日期:2011-02-11作者简介:李延升(1971-),男,山东章丘人,讲师,博士,主要从事化工传递过程与精细化工的教学研究工作.通讯作者:张保林(1947-),男,河南西平人,教授,博士生导师,主要从事化工传递过程与精细化工的教学研究工作.扩散过程是重要的质量传递方式之一,元胞自动机则是近年来新兴的仿真模拟方法.但是,将两者结合,即用元胞自动机研究扩散过程的文献却相对较少.有关此类的报道,多是和反应过程相关联且研究的侧重点放在后者[1],很少见到用元胞自动机专门研究扩散过程的报道.另外,在有关的报道中,元胞自动机模型的参数往往过多且相互之间的关系复杂,分析时需要综合多门学科的理论,这限制了它的实际应用.作为一种新兴的研究手段,元胞自动机的意图是以极其简单的规则解释或模拟复杂的现象,而有关的报道多数违背了这一意图,无法体现元胞自动机的优越性.本文即从此方面入手,将元胞自动机方法与数学分析方法相对比,给出了生动形象的扩散过程动态画面,为进一步的研究奠定了基础.1 元胞自动机简介元胞自动机有时也被称为细胞自动机、点格自动机、分子自动机或单元自动机,它是现代计算机之父N eum ann 及其追随者提出的想法.20世纪末21世纪初,Stephen 将这种带有强烈的纯游戏色彩的原始想法从学术上加以分类整理,最终使之上升到了科学方法论[2].元胞自动机是一时间和空间都离散的动力系统,散布在规则格网(Lattice Grid)中的每一元胞(Ce ll)取有限的离散状态,遵循同样的演化规则,依据确定的局部规则同步更新,大量元胞通过简单的相互作用而构成动态系统的演化.不同于一般的动力学模型,元胞自动机不是由严格定义的物理方程或函数确定,而是由一系列模型构造的规则构成,凡是满足这些规则的模型都可以算作元胞自动机模型.因此,元胞自动机是一类模型的总称,或者说是一个方法框架[1-2].元胞自动机最基本的组成为元胞、元胞空间、邻居及规则这4部分.简单来讲,元胞自动机可以视为由一个元胞空间和定义于该空间的变换函数所组成[1-3].2 扩散过程元胞自动机模型的提出同普通的元胞自动机一样,本文的元胞自动机也是采用等间隔的点作为元胞.不失一般性,以二维的圆形区域为离散域,以处于中心的某个元胞为圆心,取适当长度的半径,在圆域内的元胞即为该元胞自动机的元胞,正好在圆上的交叉点也看作圆域内的元胞,见图1.以纵横方向上均匀分布的点为元胞,圆域内外的元胞区分非常明确.为了让画面清晰,该图的元胞较少,在实际应用时,元胞数量要比该图多得多.河南工程学院学报(自然科学版)2011年图1 模拟扩散过程的元胞自动机离散域及元胞Fig .1 The d iscrete region and cells of cellu l ar au to ma ta u sed for si m u l ate th e d iffu si on p rocess现以元胞自动机分析某一典型的扩散过程.该扩散过程的假设是:(1)参与扩散的物质为A 和溶剂M.(2)扩散过程在某一圆形区域及其外界进行.(3)初始时刻,物质A 以溶液的形式,均匀地分布在该圆形区域内.随后,A 即开始由内到外扩散.显然,在圆形区域中,越靠近圆心的位置,A 的浓度越大.(4)圆形区域内同时还有物质B ,该物质不参与扩散过程,其浓度均一且不发生变化.(5)圆域外的区域视为无穷大.根据这一假设,如果设该圆形区域的半径为a,某一元胞所在的位置处半径为r ,那么在r a 处,A 的浓度为0.根据扩散过程的特点,元胞的设定及演化规则是:(1)以元胞0、1、2分别代表物质M 、A 、B ,再用k 代表演化次数.显然,k 和时间t 对应,也是变量.(2)用元胞1的随机行走代表物质A 的扩散.具体规定是:第k 次演化时,某元胞1的上、下、前、后、左上、左下、右上、右下的8个邻居中,若有n 个为元胞0(n 1),则在第k +1次演化时,该元胞1与这n 个元胞0的任意一个交换位置.若该交换完毕后,元胞1的位置已经在圆域的边界处,则该元胞1演化为元胞0.(3)元胞2不发生演化.(4)根据前文的假设(5),圆域外界只有元胞0,无元胞1和元胞2.3 结果及讨论3.1 演化结果用一个实例分析.设初始时刻,元胞自动机中,元胞0、1、2所占的量分别为40%、20%、40%,且均匀分布,元胞自动机的离散域是直径为100的圆,用M atlab 编写程序得到的演化图见图2,该图可模拟扩散过程且较为生动形象.图2 模拟扩散过程的元胞自动机演化示意图Fig .2 Sk etch maps si m u lating the d iffusion process w ith ce ll u lar au tom ata定义:(1)演化次数为k 时,圆域内已经扩散到外界的元胞1的数量占初始时刻元胞1的数量的比例为累计扩散率.(2)累积演化的元胞1的数量占初始时刻数量的99%时,对应的演化次数k 为ke,k e 即为扩散 2第1期李延升,等:扩散过程的元胞自动机模拟的终点.元胞1的累计扩散率曲线k e /t-k 见图3.从图2与图3可以看出,随着k 的增大,累计扩散率曲线由陡变缓,说明扩散速率由大逐渐变小.图3 元胞自动机模拟的累积扩散量曲线以及用微积分求得的累积扩散量率曲线F ig .3 T he accu m u l ate diffuse a m ount curves obta i ned w ith cell u lar au to m ata and the calcu lous res p ective l y改变元胞自动机的参数,如元胞的数量、比例等,运行程序,得到的累计扩散率曲线基本不变,说明该曲线只受演化规则的影响,与元胞自动机的参数无关.3.2 微积分分析再用微积分的方法分析这一扩散过程,并与元胞自动机方法对比.该扩散过程实际上是球形扩散体系在球的大圆截面上的扩散过程.为了更符合实际情况,以球坐标分析.取球心为坐标原点,建立球坐标.设球体内,物质A 的浓度为c A ,则c A 是空间半径变量r 和时间变量t 的函数,角度方向上无扩散分量.物质A 的传递方程为c A t =D AB 1r 2 r (r 2 c A r ),(1)初始条件:c A =c 0,(2)边界条件:c A (l ,t)=0,(3)以上三式中,c A 与元胞1的数量对应;r 与元胞自动机某处的位置对应;t 为时间变量,与演化次数k 对应;D AB 为物质A 在物质B 中的扩散系数,与元胞1的随机行走的情况对应.这三式组成的偏微分方程的解为c A (r ,t)=2 c 0 r !∀n=1(-1)n+1sin (n r / )n exp (-n 2 2D AB2t),r #0,n =1,2,3,∃2c 0!∀n=1(-1)n+1exp (-n 2 2D AB2t),r =0,n =1,2,3,∃(4)由此,可求得物质A 的扩散速率N A (t)和累计扩散率x 1:N A (t)=-D AB c A r |r= =2D AB c 0 !∀n =1exp (-n 2 2D AB2t),n =1,2,3,∃(5)x 1=%t 04 2N A (t)dt /(43 3c 0)=6 2!∀n =11n 2[1-exp (-D AB n 2 2t 2)],n =1,2,3,∃(6)为了进一步阐述的需要,兹列举一个实例,该实例来自文献[6],实例为球形的晕海宁胶囊的扩散控制的释放系统,其药物扩散体系符合本文的假设.胶囊的直径为3.26mm ,晕海宁在控释凝胶层中的扩散系数D AB 为3.0&10-7m 2/s ,由此得x 1=6 2!∀n=11n 2[1-exp (-2.78&10-5n 2t)].(6a)3河南工程学院学报(自然科学版)2011年用M a tlab辅助计算,取n=1000,即x1∋62!1000n=11n2[1-exp(-2.78&10-5n2t)],(6b)则截断误差r1000的范围为r1000<62!∀n=11n2-62!1000n=11n2=1-0.9994=6&10-4.截断误差已经足够小,得到的累积扩散率曲线x1/te-t亦示于图3中,便于对比.图中,te为累积释放率达到99%时的t值,与元胞自动机的ke对应.从图3可以看出,两条曲线非常接近.实际上,任意取一D AB的值,并作x1/te-t曲线,结果都是如此,从而说明用元胞自动机的研究结论与用微积分研究的结论是一致的.4 结 论以元胞自动机作为方法和手段,通过编写M atlab程序,可以形象、逼真地模拟扩散过程.元胞自动机的核心是,各参数与扩散现象紧密相关,用元胞的演化模拟扩散的过程.元胞自动机模拟的结果与微积分的结论一致,并由实例为佐证.用元胞自动机模拟扩散过程,虽然具有规则简单、模拟效果较直观等优点,但是所用元胞自动机的演化规则和重要参数,基本上是根据释放的过程而作的人为规定,可以说没有完全脱离经验模型的框架和模式.当然,目前绝大多数元胞自动机模拟的情况也都是这样,这也是元胞自动机尚未被广泛应用的重要原因之一.如何通过理论分析得到演化规则和演化参数,使该模拟过程上升到较高的理论层次,属于进一步研究的内容.参考文献:[1] Basti en Chopa rd,M e ichel D roz.物理系统的元胞自动机模拟[M].祝玉学,赵学龙,译.北京:清华大学出版社,2003:1-47.[2] Stephen W o lfra m.A new k i nd o f sc i ence[M].W o lfram m edia,2002:1-50.[3] 曹伟.元胞自动机与计算机模拟[J].丹东纺专学报,2005,12(2):1-4.[4] Crank J.The m a t he m atics o f diff usion[M].O x f o rd:O x f o rd university press,1975:89-103.[5] 李延升,张保林,张雪梅,等.膜控型缓控释肥料养分释放模型的研究[J].江苏农业学报,2009,25(5):1033-1038.[6] W e lty J R,W icks C E,W ilson R E.动量、热量和质量传递原理[M].马紫峰,吴卫生,译.北京:化学工业出版社,2005:338-341.The Sim ulation of D iffusion P rocess Based on C ell u lar A uto m at aLI Yansheng1,2,HOU K eke1,ZHANG Baoli n2(1.School of Che m istry&Che m icalE ngineering,Xuchang C ollege,X uchang461000,China;2.School of Che m icalE ng ineering&Ener gy Sources,Zhengzhou University,Zhengzhou450007,China) Abstrac t:T ak i ng ce ll ular auto m a ta as research me t hod and m atl ab prog ram as t he m eans,the diff usion process i n round c ircu m scri p tion was si m ulated in the for m o f evo l u tion s ketch m aps.M ean w hile,by m eans of ca lcu l us,t h i s diffus i on process was ana l y zed associati ng w ith ex a m ple and comparing w ith cell u l ar auto m ata me t hod.The conc l usi on w as t hat the ce ll u lar auto m a ta a re a ble to si m u l a te the diff usion process viv i dly.T he resu lt by m eans o f cell u l ar au t om ata and t hat by m eans o f ca lcu l us ma tches each o t her.K ey word s:d iffus i on;ce ll u l a r auto m a ta;diff usion equation;si m u lati on4。