渗流模型的计算机模拟_毕业设计论文.doc

PAGE 毕业设计（论文）任务书 题目：渗流模型的计算机模拟 一、原始依据 1、论文的工作基础： 1逾渗是统计物理中的基础理论，在广泛的体系中得到应用。然而，除了极少的规则点阵，大多数结构的逾渗阈值都是依赖计算机模拟获得的。不规则结构的逾渗研究尤其薄弱。近年来，R.Ziff提出了计算机模拟新算法，并应用于正方形点阵和立方体点阵；该算法的模拟计算效率较高。 本课题组在前期工作中提出了一种新的计算机模型：RCP-LV模型，用于模拟多晶材料等无序不规则胞状结构；该模型的数据结构完整易用。这种结构的逾渗问题，还没有研究过。 2研究条件：电脑：cpu：P4 2.0GHz 内存：2；应用软件：Matlab6.5 3工作目的： 研究RCP-LV模型的基本问题。具体为 （1）设计统计发生逾渗时（估计）的逾渗值的方法和程序，用matlab编写程序并进行调试。 （2）对这些逾渗值进行进一步处理，利用统计学关系将这些统计结果经过推导和计算获得正则下的逾渗值。 （3）探讨提高Pc精度的方法。 （4）统计逾渗的最大集团，并进行结果分析。 二、参考文献 [1] R.泽仑．非晶态固体物理学[M]．北京：北京大学出版社，1988．155—217． [2]M. E. J. Newman, R. M. Ziff. Efficient Monte Carlo algorithm and high-precision results for percolation [J]．Physical review letters. 2000, 85(19): 1—4. [3]M. E. J. Newman, R. M. Ziff. Fast Monte Carlo algorithm for site or bond percolation [J]. Physical review E, 2001, 64(016706): 1—16. [4] Geoffrey Grimmett. Percolation [M]. Beijing: World Publishing Corporation，2000, 6(1), 1992 . 6-10． [5] G.R. Jerauld, L.E Scriven and H.T. Davis. Percolation and conduction on the 3D voronoi and regular networks: a second case study in topological disorder[J].J. Phys. C: Solid State Phys. 1984, 17(22): 1—10. [6]范智刚，吴裕功，赵选贺．Laguerre－Voronoi图软件包的设计和实现[J]．天津大报，2003，36(6): 1—2. [7] Scott Kirpatrick. Percolation and Conduction [J]. Reviews of modern physics，1973, 45(4): 1—13. 三、设计（研究）内容和要求 在本课题组的系列lognormal分布的Voronoi图上研究逾渗相关问题。在已经计算出的lognormal分布的Voronoi图上利用已有的数据结构，设计研究逾渗所需的数据结构，完成查找和连接等算法模块，对逾渗现象作模拟研究。具体目标参数如下所示。 主要指标和技术参数： 1、统计集团大小，比如发生逾渗时最大的集团大小，以及平均集团大小。 2、设计计算（估计）Pc的方法和程序，用matlab编写程序并进行调试。 3、尽量将Pc的精度提高（目前Pc值在各种模型即使是正方形模拟都没有精确值）。 具体要求：1、了解逾渗模型的基础理论 2、在本课题组的Voronoi图上，利用matlab建模进行仿真运算。 3、在本课题组的Voronoi图上，通过计算机实验，统计并通过推导和计算估计Pc值，研究提高其精度的方法。 4、统计最大集团的性质。 指导教师（签字） 年 月 日 审题小组组长（签字） 年 月 日 摘 要 本文采用了一种新的模型——LV（Laguerre-Voronoi）图——来进行逾渗试验，由于还没有相关的报道，这使得这次试验更有意义。 本文中，本文首先介绍了关于逾渗模型的一些概念，比如逾渗的定义，逾渗值，平均集团大小，以及各种应用。基于ziff的算法本试验给出了适于LV模型的算法以及计算机程序。在PC（cpu：P4 2.0G 内存1G）上完成整个模拟试验。 具体步骤是首先统计在不同规模下不同的逾渗值的大量数据（本实验中每个规模计算了104次），然后利用这些数据画出Rl（n）图，并记录下函数Rl（n）的数据，