絮凝体的DLA分形仿真模拟及其分形维数的计算方法论文.pdfVIP

絮凝体的DLA分形仿真模拟及其分形维数的计算方法珠 金鹏康 王晓昌 郭坤 (西安建筑科技大学西北水资源与环境生态教育部重点实验室，西安，710055) 摘要：论文对絮凝体的成长过程，运用有限扩散凝聚(DLA)模型，进行了二维仿真模拟。对仿真模拟所得到的分 形图像，分别采用密度函数法、回转半径法、图像分析法进行分析计算，得到的分形维数均随模拟絮凝体的尺寸增 加而减小。三种计算方法的结果显示，密度函数法和回转半径法所得到的分形维数基本相当，但图像分析法所得到 的分形维数较小，约为密度函数法和回转半径法所得到的分形维数的0．8左右，其原因可能是图像分析过程中计算 机对图像的识别误差引起的。以DLA模型模拟所得到的絮凝体内部的孔隙率随絮凝体尺寸的增加而增加。絮凝体中 孔隙率的增加是絮凝体密度减小、结构松散的主要原因，也是絮凝体分形维数降低的主要因素。 关键词：絮凝体；有限扩散凝聚模型；仿真模拟；分形维数 ～般情况下，絮凝体的结构具有分形构造 特征[1’2】，絮凝体的分形形态直接影响着水处 理固液分离效果，如果絮凝体较为密实、形态 结构紧凑，则易于沉降。对水中粒子的絮凝过 程进行仿真模拟是目前探索絮凝体内部构造 特征的重要途径之一，得到的模拟絮凝体可以 为改善絮凝体构造，提高絮凝体密度以及水处 理固液分离提供理论依据，其中通常被采用的 模型是有限扩散凝聚(DLA)模型。另一方面， 描述絮凝体结构密实的方法是测定其分形维 数，分形维数越大，结构越密实。分形维数的 计算方法较多，不同的分析方法间可能会有较 大的出入，但迄今为止的研究尚未有关于不同 的分形维数分析方法的比较，本文对絮凝体的 成长过程，运用有限扩散凝聚(DLA)模型， 进行了二维仿真模拟，对模拟絮凝体分别采用 密度函数法、回转半径法、图像分析法进行了 分形维数计算，评价了絮凝体成长过程中结构 的变化特征以及三种方法的分析特点。 1 DLA模型 在1981年，T．A．Witten和L．M．Sander 提出了DLA模型[3】，即有限扩散凝聚模型 (Diflusion—Limited Aggregation)，是用来研 究悬浮在溶液或大气中的金属粉末、煤灰和烟 尘等微粒的无规扩散凝聚过程。DLA模型一 般指在有限空间内有一个中心粒子内，其它粒 子在该有限空间内的不断的自由运动，这些粒 子以随机行走的方式与该中心粒子结合，逐渐 图1 DLA模拟算法结构框图 形成某一粒径的颗粒物，由于微粒的运动受到 。旗金项II：同家自然科学甚金20547001) 106 有限空间的限制，其运动方式会受到空间边壁的影响。 2仿真计算方法 LxL的方形点阵，在该二维平面中央设置一个固定的粒子，也称之为“种子”。开始时，在平面的边 缘上随机地产生一个粒予，该粒子以布朗运动方式在平面上作无规运动。假定粒子的运动有两种可 能：一种是与“种子”相碰并附着于“种子”上，形成粒子簇；再随机地产生一个粒子，继续作无规运 动。另一种是运动到平面的边缘上，令其停止运动并消失，或者根据周期性边界条件继续作无规运 动。将上述过程重复下去，直至最后在中央形成一个树枝状的凝聚体，其算法结构框图如图1。 3分形维数的计算方法 目前，对于DLA生长模型所形成的分形图像主要采用密度相关函数法和回转半径法两种方法。 (1)密度相关函数法： 密度相关函数法一般认为分形维数服从以下规律： C(r)一r一“ (1) 式中C纠为密度相关函数，它与围绕一个给定点的空间中密度分布成正比。a是一个大于零小 (d为欧氏维数)。C纠和r可以通过模拟过程中的数据给出。 (2)回转半径法 回转半径法认为模拟絮凝体的分形维数服从以下规律： Ⅳ(r)．rD (2) 式中r为囤转半径，N纠为其内的粒子数，D为分形维数。 回转半径指粒子在凝聚过程中得到的模拟絮凝体与粒子结合所能形成的最大半径，显然，随着 模拟絮凝体尺寸的增大，回转半径也不断增加，相当于不断增长的同心圆半径(如图2所示)。因此， 可以由Ⅳ纠与r的双对数坐标曲线图的斜