有限扩散集团凝聚模型第二讲： DLCA模型的形态分析-CSDN博客

本文链接：https://blog.csdn.net/u013631121/article/details/76407376

上回书说道有限扩散集团絮凝模型（DLCA），是通过走碰粘沉的方式，让胶体粒子运动相碰，最终粘贴到一块，沉降到水底的过程。我们通过MATALB，已经实现了该粒子的运动状态模拟，然而，我们每次模拟的结果是不一样的，为啥呢，因为布朗运动是一种无规则的运动，所以每次模拟粒子走的位子都不是一样的，所以到最后导致生成的图形也不同。

既然每次都不同，我们咋知道那次好那次坏？那么需要进行分析——对粒子的形态。分析方法有好多好多，那么我先列举一下方式，然后一一的讲解。

形态分析的方式有：回转半径分析、孔隙率分析、回转分形维度、粒度分形维度。

（1）回转半径：首先讲一下回转半径，如图1所示，回转半径就是某胶体粒子与中心粒子的最大距离，这里的中心一般取所有的粒子坐标的平均值，举个例子，比如有三个粒子，它们的坐标分别是(1,2,1)、(2,3,2)、(3,4,3)，那么他们的中心就是((1+2+3)/3,(2+3+4)/3,(1+2+3)/3)=(2,3,2)，也就是各轴坐标之和的平均值。然后我们在计算，每个坐标距离(2,3,2)的距离，d1=√{(1-2)^2+(2-3)^2+(1-2)^2}=√3、d2=0、d3=√3，于是我们取距离的最大值√3作为回转半径。

图1

回转半径这个词不知道是谁瞎编的，这词很容易把人给误导成理论力学中的回转半径，非常的不好。我认为这不就是粒子之间的距离吗，最后得到的最大结果就好像是一个圆里面的半径，所以应该叫半径是才是最机智的。

下面说一下我们如何用它做统计呢，能统计什么？其实说白了，就是粒子碰撞到一块，变成一个絮凝体堆，当这些东西碰撞的差不多了，不怎么再动了，我们就数数这个有多少个粒子在这个堆里面，然后在计算一下它的回转半径，然后我们将每一个堆的粒子数目作为横坐标，将回转半径作为纵坐标，给他拟合出一条曲线，之后我们就再也没必要模拟了，想知道粒子数目为某数的絮凝体的半径是多少，我们直接就拿公式计算即可。

图2

（2）孔隙率：孔隙率指的是这些粒子粘附到一块了，大家看图1，里面的各个黑点虽然在一块了，但是之间还是有缝隙，这些缝隙就是孔隙率，孔隙率如果很大，那么水就能流入里面，这东西就会变重，然后沉降到水底就快。

孔隙率是怎么算的呢，孔隙率S=1-(N*r^3/R^3)，N值得是絮凝体中的粒子个数，r为粒子的半径，R为絮凝体的回转半径。以（1）中的粒子为例，絮凝体的回转半径是√3，每个粒子都假设半径是0.5，粒子的个数是3，那么孔隙率就是1-3*0.5^3/√3^3=0.928=92.8%。这说明孔隙率很大。你再看那三个点坐标，确实距离挺大的。

我们怎么统计比较好呢，有两种，第一种是随便选一个絮凝体，观察孔隙率随着时间的变化，但有一个缺点，就是不能掌控全局。第二种是将所有絮凝体的孔隙率求一个平均值，然后再看它随着时间的变化，这个虽然能掌控全局，但好像是以偏概全。我觉得应该选择几个特殊的絮凝体做单独，并且也把全局的絮凝体做平均来统计最好。

下面我放送一张第二种方式的图图，表名孔隙率是开始越来越大，后来就稳定在40~45%之间。

图3

（3）回转分形维度：回转分型维度说的其实就是回转半径与粒子数之间的关系，刚才说过它们能够拟合成曲线，但是看上面的图2，实在是太混乱了，我们怎么办，当然是用逻辑斯蒂回归来让他像一条线啊，啥叫逻辑斯蒂？其实就是我们说的log，我们将粒子数目与回转半径分别取自然对数（ln），并作为纵坐标与横坐标，拟合出一条直线，ln(N)=Kln(R)+b，K就是回转分形维度，它代表了回转半径随着半径的增长，粒子数变化的速度。接着我们就看着回转分形维度随着时间的变化是什么样子的，其实就是说絮凝体内粒子数去留变化的平均速度是怎么样子的。下面放一张图图，说明开始速度比较波动，后来速度稳定到一个值。

图4

（4）粒度分形维度：这个概念与上面的概念相似，就是把上面的公式ln(D)=Kln(R)+b，中的N给变成D，D值指的是回转半径小于等于R的群体个数，然后K指的就是分形维度。其意思是就是粒子数目的分布情况的变化速度。接着我们就看着K随着时间的变化是什么样子的，其实就是说絮凝体平均有多少个，下面放一张图图，可见维度不断降低，即絮凝体的数目在减少（因为碰撞）。