# Simulation for Sludge Flocculation I: Brownian Dynamic Simulation for Perikinetic Flocculation of Charged Particle

Authors: Linshuang Liu; Guolu Yang; Minghui Yu
Journal: Mathematical Problems in Engineering
 (2012)

---

## Abstract

To investigate sludge drying process, a numerical simulation based on Brownian dynamic for the floc with uncharged and charged particles was conducted. The Langevin equation is used as dynamical equation for tracking each particle in a floc. An initial condition and periodic boundary condition which well conformed to reality is used for calculating the floc growth process. Each cell consists of 1000 primary particles with diameter 0.1~ 4 μm. Floc growth is related to the thermal force and the electrostatic force. The electrostatic force on a particle in the simulation cell is considered as the sum of electrostatic forces from other particles in the original cell and its replicate cells. It is assumed that flocs are charged with precharged primary particles in dispersion system by ionization. By the analysis of the simulation figures, on one hand, the effects of initial particle size and sludge density on floc smashing time, floc radius of gyration, and fractal dimension were discussed. On the other hand, the effects of ionization on floc smashing time and floc structure were presented. This study has important practical value in the high-turbidity water treatment, especially for sludge drying.


---

## Body

## 1. Introduction

There are numerous sediment-laden rivers in china. There are 42 rivers whose average annual sediment discharge is larger than 10 million. The sediment discharge of the Yellow River and the Yangtze River is 29.3% of the total top thirteen great rivers in the world. There are also many seasonal high-turbidity rivers on other regions of china. Therefore, the research of the high-turbidity water treatment has become a important research topic in china [1, 2].The effect of moisture content on the sludge volume in enormous. Moisture of sludge which taken out from the lakes and rivers was very high. This led a lot of work needed to be done before sludge can be used as resources. The moisture associated with flocs of sludge is: 70–75% free water, 20–25% floc water, and 1% each capillary and bound water. The free water can be removed by thickening. The floc water is trapped in the interstices of floc particles and is separable only by mechanical drying [3]. However, as the sludge flocculation mechanism is complex, it is still difficult to dry vast amount of sludge rapidly. As a method of dehydration, flocculation has been often used in drying engineering. Therefore, flocculation mechanism and its dynamcis simulation are very important to improve the sludge drying speed in the future studies [4].Sludge has many properties those differ from the coarse sediment and solution molecular. The floc of fine sediment particles occurs mostly in the flocculation process. In the sludge drying process, the particles grow to floc with electrochemical force [5] and adsorption [6]. These processes are very common in many hydrosol processes, especially for fine sediment flocculation. In the case of thermally stable particles, the Brownian motion should be considered as the dominant action. The aggregation characteristics are mainly affected by the morphology of aggregates [7]. As a fractal dimension can quantify the morphology of flocculation, it becomes a key variable for analyzing the flocculation dynamics. Ermak [8] simulated the flocculation process of charged particles in solution using Brownian dynamics. With more new applications for solid particles in solution recently, it is important to control the flocs size and its morphology in flocculation process. Xiong et al. [9] simulated the growth of game colloidal particle in the sectional colloid model. However, their results were obtained using the assumption of spherical particles. More work is necessary to better explain the realistic case. In this paper, we tried to simulate the motion of particles and analyse the effect of electrostatic force on the morphology in the sludge drying process. Therefore, we introduced a model that could describe the electrostatic force among charged sludge particles. This electrostatic simulation model was coupled with the Brownian dynamic simulation to deal with the flocculation of charged sludge particles. A comparison of floc growth process with different sludge density and particle diameters of charged and uncharged particles is present.


## 2. Theory

### 2.1. Langevin Equation

Brownian dynamics describe the dynamic behavior of a sediment particle whose mass and size are larger than those of the water molecules they are immersed in. These particles are subject to stochastic collisions with the water molecules, which leads to the random motion of the particles. It is assumed that a particle is in thermal equilibrium and surrounded by water molecules at a specific temperatureT. The particles in this system are subject to an external force Rn. In this case, the translational motion is described by the Langevin equation. The strict Langevin equation for the motion of a particle on the x,y,z coordinates is given by (2.1) together with several supplementary conditions [10, 11]:(2.1)d(mun)dt=Fn-mβun+Xn(n=x,y,z),(2.2)〈Xn〉=0,〈Xn2〉=2βmkbT,〈un(t)⋅un(0)〉=0(n=x,y,z).In (2.1), the left term of equation represents the Brownian particle’s momentum, and the Rn is systematic external force acting on the particle. In perikinetic flocculation, this term includes van der Waals force, electrostatic force, magnetic force, image force, and so forth. The other forces except for the electrostatic force can be negligible in the flocculation of charged sediment particles because the other forces are much smaller than the electrostatic force [5]. The specific calculation of the electrostatic force will be executed in Section 3. The second right term of equation is the frictional resistance of the fluid around the particle. β is the friction constant, and it’s the inverse of particle relaxation time. The last term Xn is an external force caused by the random collisions between the particle and the surrounding water molecule. The random force Xn is represented by the Gaussian random-distribution function. The mean and square mean of it are represented as (2.2).


### 2.2. Solution of the Langevin Equation

Several methods have been used to integrate the Langevin equation [12]. In this study, we use the Monte Carlo method. Integration of (2.1) (2.2) from t tot′ yields the equation:(2.3)un-u0ne-βt-Fn0mβ⋅(t-e-βt)=1m∫0te-β(t-t′)Xn(t0+t′)dt′(n=x,y,z),rn-rn0un0β⋅(1-e-βt)-Fn0mβ[t-1β(1-e-βt)]=1mβ∫0t[1-e-β(t-t′)]Xn(t0+t′)dt′.
Equation (2.3) cannot be used directly to calculate the velocity and position changes due to the random nature of X(t). We calculate the statistical properties of the two integrals involving X(t) with one method derived by Luo et al. [11](2.4)un=u0ne-βt+F0nmβ⋅(t-e-βt)+B1n,rn=rn0+un0β⋅(1-e-βt)+Fn0mβ[t-1β(1-e-βt)]+B2n(n=x,y,z),
where B1n,B2n are random function of time chosen from the bivariate Gaussian distribution with properties(2.5)〈B1n〉=〈B2n〉=0,〈B1n2〉=3kTm(1-e-2βt),〈B2n2〉=kTmβ2(2βt-3+4e-βt-e-2βt),〈B1n⋅B2n〉=0(i=1,2,n=x,y,z).


#### 2.2.1. Floc Structure

The morphological state of the flocs is defined as the fractal number by Forrest and Witten [13]. They confirmed that there was a power law relation between the characteristic length and the number of primary particles in each floc. Fractal structure of floc has been researched by many researches [14, 15]. Generally, the radius of gyration of the flocculation is used as the characteristic length for the floc. The power law relation is expressed as(2.6)Rf=1n∑1n(rj-1n∑1nrj),2
where n is the number of primary particles in an floc, Rf is the radius of gyration.


## 2.1. Langevin Equation

Brownian dynamics describe the dynamic behavior of a sediment particle whose mass and size are larger than those of the water molecules they are immersed in. These particles are subject to stochastic collisions with the water molecules, which leads to the random motion of the particles. It is assumed that a particle is in thermal equilibrium and surrounded by water molecules at a specific temperatureT. The particles in this system are subject to an external force Rn. In this case, the translational motion is described by the Langevin equation. The strict Langevin equation for the motion of a particle on the x,y,z coordinates is given by (2.1) together with several supplementary conditions [10, 11]:(2.1)d(mun)dt=Fn-mβun+Xn(n=x,y,z),(2.2)〈Xn〉=0,〈Xn2〉=2βmkbT,〈un(t)⋅un(0)〉=0(n=x,y,z).In (2.1), the left term of equation represents the Brownian particle’s momentum, and the Rn is systematic external force acting on the particle. In perikinetic flocculation, this term includes van der Waals force, electrostatic force, magnetic force, image force, and so forth. The other forces except for the electrostatic force can be negligible in the flocculation of charged sediment particles because the other forces are much smaller than the electrostatic force [5]. The specific calculation of the electrostatic force will be executed in Section 3. The second right term of equation is the frictional resistance of the fluid around the particle. β is the friction constant, and it’s the inverse of particle relaxation time. The last term Xn is an external force caused by the random collisions between the particle and the surrounding water molecule. The random force Xn is represented by the Gaussian random-distribution function. The mean and square mean of it are represented as (2.2).


## 2.2. Solution of the Langevin Equation

Several methods have been used to integrate the Langevin equation [12]. In this study, we use the Monte Carlo method. Integration of (2.1) (2.2) from t tot′ yields the equation:(2.3)un-u0ne-βt-Fn0mβ⋅(t-e-βt)=1m∫0te-β(t-t′)Xn(t0+t′)dt′(n=x,y,z),rn-rn0un0β⋅(1-e-βt)-Fn0mβ[t-1β(1-e-βt)]=1mβ∫0t[1-e-β(t-t′)]Xn(t0+t′)dt′.
Equation (2.3) cannot be used directly to calculate the velocity and position changes due to the random nature of X(t). We calculate the statistical properties of the two integrals involving X(t) with one method derived by Luo et al. [11](2.4)un=u0ne-βt+F0nmβ⋅(t-e-βt)+B1n,rn=rn0+un0β⋅(1-e-βt)+Fn0mβ[t-1β(1-e-βt)]+B2n(n=x,y,z),
where B1n,B2n are random function of time chosen from the bivariate Gaussian distribution with properties(2.5)〈B1n〉=〈B2n〉=0,〈B1n2〉=3kTm(1-e-2βt),〈B2n2〉=kTmβ2(2βt-3+4e-βt-e-2βt),〈B1n⋅B2n〉=0(i=1,2,n=x,y,z).


### 2.2.1. Floc Structure

The morphological state of the flocs is defined as the fractal number by Forrest and Witten [13]. They confirmed that there was a power law relation between the characteristic length and the number of primary particles in each floc. Fractal structure of floc has been researched by many researches [14, 15]. Generally, the radius of gyration of the flocculation is used as the characteristic length for the floc. The power law relation is expressed as(2.6)Rf=1n∑1n(rj-1n∑1nrj),2
where n is the number of primary particles in an floc, Rf is the radius of gyration.


## 2.2.1. Floc Structure

The morphological state of the flocs is defined as the fractal number by Forrest and Witten [13]. They confirmed that there was a power law relation between the characteristic length and the number of primary particles in each floc. Fractal structure of floc has been researched by many researches [14, 15]. Generally, the radius of gyration of the flocculation is used as the characteristic length for the floc. The power law relation is expressed as(2.6)Rf=1n∑1n(rj-1n∑1nrj),2
where n is the number of primary particles in an floc, Rf is the radius of gyration.


## 3. Numerical Simulation

### 3.1. Particle Charging

When fine sediment dispersed in water, they were charged due to ionization and sorption [16, 17]. In this simulation, we assumed that primary particles were charged with some additive ions to simulate the flocculation with different potential flocculent. The equilibrium and diffusive properties of charged spherical particles in solution have been studied by Xiangyan et al. [8], who simulated the motion of particles in solution obeying classical statistical mechanics. In this case, the charge quantity of a floc is the sum of elementary charges on the primary particles. The quantity of electricity on each particle is calculated by ionization intensity formula. In (3.1), Ci is particle molarity, I is ionization intensity, and Zi is quantivalency of particle. In our research, Ci was relate to sludge density and the span of I was 0.001~0.1.(3.1)I=12⋅∑1nCi⋅Zi2.


### 3.2. Initial Condition

To simulate the flocculation process, initial conditions for the position and the velocity of the primary particles have to be given in Figure1 and Figure 2, respectively. The positions of these particles are obtained with a homogeneous distribution in 3 dimensions. It can be assumed that the primary particles experience only thermal force in the initial state. Therefore, we set the initial velocities with a normal random number generator that is characterized by the Gaussian distribution of zero mean and kbT/m0 variance.Figure 1
Particle initial position alongx,y,z axis.Particle initial velocity alongx,y,z axis.
(a)
ux0(b)
uy0(c)
uz0The primary particles diameters were chosen from0.1~4μm [18, 19], and space between two particles L0 were calculated with sludge density to coincided to real liquid dispersion. The space is decreased with the increase of sludge density ρs as well as particle diameters ρp as shown in Figure 3.Figure 3
Space between particles with different sludge density.We calculated the trajectories of 1000 primary particles in a simulation cell. The cell size was determined by the sludge densityρs. Ten simulations for each initial condition were performed to obtain statistically meaningful average values. The particle diameter dp was value to 0.1, 0.5, 1.0, 2.0, and 4.0 um separately and ρs was value to 1100, 1075, 1050, 1025, and 1010 kg m−3. In this work, we have fixed the values of T, particles density(ρp) and Δt to be 293 K, 2650 kg m−3, and 1/50 s, respectively. These conditions were obtained from the experimental conditions in a flocculation process.


#### 3.2.1. Boundary Condition

In this simulation, 300 particles were taken as boundary smash condition to signify the termination of floc growing. The periodic boundary condition was used not only for maintaining the internal bond strength but also for calculating the electrostatic force. The electrostatic force of theith particle is an infinite summation of each force from the other charged particles in the cell. Therefore, this force includes the electrostatic dispersion effect mentioned in the paper of Ermak [8]. The total electrostatic force on particle i is given by(3.2)F0i=∑j=1,j≠i∞Keqiqjrij2r⃗ij.In this paper, the electrostatic force can be calculated directly from the gradient of the pair potential. Ermak [8] introduced a new technique for treating an infinite number of pair potential. This potential includes the interaction with the images of the original cell in the replicated cells.Though an infinite number of particles affect the electrostatic force of a given particle, the force has a finite value because it is symmetrical in a macroscopic sense. In our simulation, electric charges locate in the center of mass of the floc and the dielectric constant of the water is uniform. By applying this method to the present system, the modified electrostatic force by the infinite number of particles can be represented as:(3.3)F0i=-∑j=1,j≠iNKeq2[∇ψ1(rij)+∇ψ2(rij)⃗],
where N is the total number of flocs in the original, cell;q  is quantity of charges on the particles. ψ1 and ψ2 are expressed as(3.4)ψ1(r)=erfc⁡(πr/L)r-1L-EmL,(3.5)ψ2(r⃗)=1L∑1⃗*(erfc⁡(π|(r⃗/L)+l⃗|)|(r⃗/L)+l⃗|)+exp⁡(-πl2)cos⁡(2πl⃗(r⃗/L))πl2,
where L is the length of the side of cell; Em is the Madelung constant and its value, for a simple cubic lattice, is Em = 2.83729479 [20],l̅ is a vector with integer components, l is its magnitude, and * denotes exclusion of thel̅=0. To treat the equation of (3.5), we have used the optimized expansion in Cubic harmonics, which was first introduced by Xiangyan et al. [15]. Two basic matrix P,Q in (3.7) are used in coefficient matrix [k1,k2,k3,k4,k15,k6,k7,k8,k9,k10,k11,k12,k13]T to calculated the ∇ψ1 and ∇ψ2. In this method, 15 coefficients are listed in (3.8) by Hansen in order to optimize the (3.9)
(3.6)P=[x1x2x3]，Q=[x1x2x3],(3.7)[k1k2k3k4k5k6k7k8k9k10k11k12k13]=[PQ⋅(PQ)TPQ3⋅(PQ3)Texp(-πP⋅PT)|Q2|14a46a68a410a44b46b68b410b4b4b6b8b10a4a6a8a106c68c810c1012c12c6c8c10c128d810d1012d12d8d10d12]⋅[1P⋅PT(P⋅PT)2(P⋅PT)3(P⋅PT)4(P⋅PT)5(P⋅PT)6],(3.8)[a4b4c6a6b6c8d8a8b8c10d10a10b10c12d12]=[-3.0901509.8010671-22.6452014.6764131-19.994017-656.4341357.7471497.681783-237.90423440.65802-570.95501-110.95681312.8662996.932540133.22453].Consequently, the gradients of∇ψ1 and ∇ψ2 are represented as:(3.9)∇ψ1(rij)/n=1L2(2exp⁡(-πP⋅PT)P⋅PT+erfc⁡(πP⋅PT)(P⋅PT)3/2)xn⋅en̂,∇ψ2(rij)/n=-1L2(m11xn+m2xn+m3xn3+m4xn7)⋅en̂n=1,2,3,
where m1,m2,m3,m4 are described as(3.10)m1=-2k3k11k4k52,m2=2πk3k9-k3k6+2πk3k8k1-k3k7k1k52+4πk3k8k1+2πk3k11k4-k3k10k4k53+6k3k11k4+2πk3k13k2-k3k2k12k54+k2k3k13k2k55,m3=-4k3k8k52,m4=-8k3k13k54.


### 3.3. Calculation Procedure

We define that a particle collides with the floc if the distance between two is located within the radius of gyration of floc. The momentum and the charge quantity are conserved before and after the collision, and immediately the new-formed floc is in thermal equilibrium state. We assumed that the sticking probability was 1 after each collision. The flocs grow with the repetition of the above sequence. The program was written with Fortran 95 and the procedure is shown in Figure4.Figure 4
Flow chart for Brownian dynamic simulation.


## 3.1. Particle Charging

When fine sediment dispersed in water, they were charged due to ionization and sorption [16, 17]. In this simulation, we assumed that primary particles were charged with some additive ions to simulate the flocculation with different potential flocculent. The equilibrium and diffusive properties of charged spherical particles in solution have been studied by Xiangyan et al. [8], who simulated the motion of particles in solution obeying classical statistical mechanics. In this case, the charge quantity of a floc is the sum of elementary charges on the primary particles. The quantity of electricity on each particle is calculated by ionization intensity formula. In (3.1), Ci is particle molarity, I is ionization intensity, and Zi is quantivalency of particle. In our research, Ci was relate to sludge density and the span of I was 0.001~0.1.(3.1)I=12⋅∑1nCi⋅Zi2.


## 3.2. Initial Condition

To simulate the flocculation process, initial conditions for the position and the velocity of the primary particles have to be given in Figure1 and Figure 2, respectively. The positions of these particles are obtained with a homogeneous distribution in 3 dimensions. It can be assumed that the primary particles experience only thermal force in the initial state. Therefore, we set the initial velocities with a normal random number generator that is characterized by the Gaussian distribution of zero mean and kbT/m0 variance.Figure 1
Particle initial position alongx,y,z axis.Particle initial velocity alongx,y,z axis.
(a)
ux0(b)
uy0(c)
uz0The primary particles diameters were chosen from0.1~4μm [18, 19], and space between two particles L0 were calculated with sludge density to coincided to real liquid dispersion. The space is decreased with the increase of sludge density ρs as well as particle diameters ρp as shown in Figure 3.Figure 3
Space between particles with different sludge density.We calculated the trajectories of 1000 primary particles in a simulation cell. The cell size was determined by the sludge densityρs. Ten simulations for each initial condition were performed to obtain statistically meaningful average values. The particle diameter dp was value to 0.1, 0.5, 1.0, 2.0, and 4.0 um separately and ρs was value to 1100, 1075, 1050, 1025, and 1010 kg m−3. In this work, we have fixed the values of T, particles density(ρp) and Δt to be 293 K, 2650 kg m−3, and 1/50 s, respectively. These conditions were obtained from the experimental conditions in a flocculation process.


### 3.2.1. Boundary Condition

In this simulation, 300 particles were taken as boundary smash condition to signify the termination of floc growing. The periodic boundary condition was used not only for maintaining the internal bond strength but also for calculating the electrostatic force. The electrostatic force of theith particle is an infinite summation of each force from the other charged particles in the cell. Therefore, this force includes the electrostatic dispersion effect mentioned in the paper of Ermak [8]. The total electrostatic force on particle i is given by(3.2)F0i=∑j=1,j≠i∞Keqiqjrij2r⃗ij.In this paper, the electrostatic force can be calculated directly from the gradient of the pair potential. Ermak [8] introduced a new technique for treating an infinite number of pair potential. This potential includes the interaction with the images of the original cell in the replicated cells.Though an infinite number of particles affect the electrostatic force of a given particle, the force has a finite value because it is symmetrical in a macroscopic sense. In our simulation, electric charges locate in the center of mass of the floc and the dielectric constant of the water is uniform. By applying this method to the present system, the modified electrostatic force by the infinite number of particles can be represented as:(3.3)F0i=-∑j=1,j≠iNKeq2[∇ψ1(rij)+∇ψ2(rij)⃗],
where N is the total number of flocs in the original, cell;q  is quantity of charges on the particles. ψ1 and ψ2 are expressed as(3.4)ψ1(r)=erfc⁡(πr/L)r-1L-EmL,(3.5)ψ2(r⃗)=1L∑1⃗*(erfc⁡(π|(r⃗/L)+l⃗|)|(r⃗/L)+l⃗|)+exp⁡(-πl2)cos⁡(2πl⃗(r⃗/L))πl2,
where L is the length of the side of cell; Em is the Madelung constant and its value, for a simple cubic lattice, is Em = 2.83729479 [20],l̅ is a vector with integer components, l is its magnitude, and * denotes exclusion of thel̅=0. To treat the equation of (3.5), we have used the optimized expansion in Cubic harmonics, which was first introduced by Xiangyan et al. [15]. Two basic matrix P,Q in (3.7) are used in coefficient matrix [k1,k2,k3,k4,k15,k6,k7,k8,k9,k10,k11,k12,k13]T to calculated the ∇ψ1 and ∇ψ2. In this method, 15 coefficients are listed in (3.8) by Hansen in order to optimize the (3.9)
(3.6)P=[x1x2x3]，Q=[x1x2x3],(3.7)[k1k2k3k4k5k6k7k8k9k10k11k12k13]=[PQ⋅(PQ)TPQ3⋅(PQ3)Texp(-πP⋅PT)|Q2|14a46a68a410a44b46b68b410b4b4b6b8b10a4a6a8a106c68c810c1012c12c6c8c10c128d810d1012d12d8d10d12]⋅[1P⋅PT(P⋅PT)2(P⋅PT)3(P⋅PT)4(P⋅PT)5(P⋅PT)6],(3.8)[a4b4c6a6b6c8d8a8b8c10d10a10b10c12d12]=[-3.0901509.8010671-22.6452014.6764131-19.994017-656.4341357.7471497.681783-237.90423440.65802-570.95501-110.95681312.8662996.932540133.22453].Consequently, the gradients of∇ψ1 and ∇ψ2 are represented as:(3.9)∇ψ1(rij)/n=1L2(2exp⁡(-πP⋅PT)P⋅PT+erfc⁡(πP⋅PT)(P⋅PT)3/2)xn⋅en̂,∇ψ2(rij)/n=-1L2(m11xn+m2xn+m3xn3+m4xn7)⋅en̂n=1,2,3,
where m1,m2,m3,m4 are described as(3.10)m1=-2k3k11k4k52,m2=2πk3k9-k3k6+2πk3k8k1-k3k7k1k52+4πk3k8k1+2πk3k11k4-k3k10k4k53+6k3k11k4+2πk3k13k2-k3k2k12k54+k2k3k13k2k55,m3=-4k3k8k52,m4=-8k3k13k54.


## 3.2.1. Boundary Condition

In this simulation, 300 particles were taken as boundary smash condition to signify the termination of floc growing. The periodic boundary condition was used not only for maintaining the internal bond strength but also for calculating the electrostatic force. The electrostatic force of theith particle is an infinite summation of each force from the other charged particles in the cell. Therefore, this force includes the electrostatic dispersion effect mentioned in the paper of Ermak [8]. The total electrostatic force on particle i is given by(3.2)F0i=∑j=1,j≠i∞Keqiqjrij2r⃗ij.In this paper, the electrostatic force can be calculated directly from the gradient of the pair potential. Ermak [8] introduced a new technique for treating an infinite number of pair potential. This potential includes the interaction with the images of the original cell in the replicated cells.Though an infinite number of particles affect the electrostatic force of a given particle, the force has a finite value because it is symmetrical in a macroscopic sense. In our simulation, electric charges locate in the center of mass of the floc and the dielectric constant of the water is uniform. By applying this method to the present system, the modified electrostatic force by the infinite number of particles can be represented as:(3.3)F0i=-∑j=1,j≠iNKeq2[∇ψ1(rij)+∇ψ2(rij)⃗],
where N is the total number of flocs in the original, cell;q  is quantity of charges on the particles. ψ1 and ψ2 are expressed as(3.4)ψ1(r)=erfc⁡(πr/L)r-1L-EmL,(3.5)ψ2(r⃗)=1L∑1⃗*(erfc⁡(π|(r⃗/L)+l⃗|)|(r⃗/L)+l⃗|)+exp⁡(-πl2)cos⁡(2πl⃗(r⃗/L))πl2,
where L is the length of the side of cell; Em is the Madelung constant and its value, for a simple cubic lattice, is Em = 2.83729479 [20],l̅ is a vector with integer components, l is its magnitude, and * denotes exclusion of thel̅=0. To treat the equation of (3.5), we have used the optimized expansion in Cubic harmonics, which was first introduced by Xiangyan et al. [15]. Two basic matrix P,Q in (3.7) are used in coefficient matrix [k1,k2,k3,k4,k15,k6,k7,k8,k9,k10,k11,k12,k13]T to calculated the ∇ψ1 and ∇ψ2. In this method, 15 coefficients are listed in (3.8) by Hansen in order to optimize the (3.9)
(3.6)P=[x1x2x3]，Q=[x1x2x3],(3.7)[k1k2k3k4k5k6k7k8k9k10k11k12k13]=[PQ⋅(PQ)TPQ3⋅(PQ3)Texp(-πP⋅PT)|Q2|14a46a68a410a44b46b68b410b4b4b6b8b10a4a6a8a106c68c810c1012c12c6c8c10c128d810d1012d12d8d10d12]⋅[1P⋅PT(P⋅PT)2(P⋅PT)3(P⋅PT)4(P⋅PT)5(P⋅PT)6],(3.8)[a4b4c6a6b6c8d8a8b8c10d10a10b10c12d12]=[-3.0901509.8010671-22.6452014.6764131-19.994017-656.4341357.7471497.681783-237.90423440.65802-570.95501-110.95681312.8662996.932540133.22453].Consequently, the gradients of∇ψ1 and ∇ψ2 are represented as:(3.9)∇ψ1(rij)/n=1L2(2exp⁡(-πP⋅PT)P⋅PT+erfc⁡(πP⋅PT)(P⋅PT)3/2)xn⋅en̂,∇ψ2(rij)/n=-1L2(m11xn+m2xn+m3xn3+m4xn7)⋅en̂n=1,2,3,
where m1,m2,m3,m4 are described as(3.10)m1=-2k3k11k4k52,m2=2πk3k9-k3k6+2πk3k8k1-k3k7k1k52+4πk3k8k1+2πk3k11k4-k3k10k4k53+6k3k11k4+2πk3k13k2-k3k2k12k54+k2k3k13k2k55,m3=-4k3k8k52,m4=-8k3k13k54.


## 3.3. Calculation Procedure

We define that a particle collides with the floc if the distance between two is located within the radius of gyration of floc. The momentum and the charge quantity are conserved before and after the collision, and immediately the new-formed floc is in thermal equilibrium state. We assumed that the sticking probability was 1 after each collision. The flocs grow with the repetition of the above sequence. The program was written with Fortran 95 and the procedure is shown in Figure4.Figure 4
Flow chart for Brownian dynamic simulation.


## 4. Results

### 4.1. Floc Growth Process

#### 4.1.1. Number of Particles in Floc

A series of comparison ofdp from 0.1 to 4.0 um were performed to interpret floc growth process. We chose ρs= 1010, 1025, 1050, 1075, and 1100 kg m−3 for simulation in order to intergrate with practice. As shown in Figure 5, in the early stage of flocculation the primary particle was sustaining to the a initial flocculent time. The particle number begun to increase after the first stage till reaching the smashing condition. In the final stage of flocculation, the sequence of smashing time with different sludge density is similar to that of initial flocculent time. The initial flocculent time and smashing time are both increased as the floc diameter increased.Variety of particle number as floc growing.
(a)
dp=0.1(b)
dp=1um(c)
dp=4um


#### 4.1.2. The Smashing Time

The floc size can not grow up unlimited in the practice. It could be affect by particle diameter and particle-bonding mode. In this paper, the time that 300 particles conglutinated in one floc is smashing time. Though smashing time was related to random velocity and displacement, it increased with the decrease of sludge density. Whereas it increased with the increase of particle diameter as shown in Figures6 and 7. This tendency may be understood, on one hand as particle diameter presents a linear increase, the particle number decreased in a cube way, which leads particle impact probability to decrease. On the other hand, the decrease of sludge density causes particle impact probability to decrease likewise.Figure 6
Effects of sludge density on smashing time.Figure 7
Effects of particle diameter on smashing time.


#### 4.1.3. Radius of Gyration

The relation between the radius of gyration and particle number was not linear positive correlation. The main factors influencing it were particle number and fractal dimension. In this simulation, the radius of gyration is calculated by four variables, they were particle number andx-, y-, z-axis coordinate values. The radius of gyration values in Figure 8 was calculated by three duplicate computations with the initial particle size being equal to 2 um. It was shown in the Figure 8 that the gyration radius decreased with the increase of sludge density. The relation formulas of sludge density to gyration radius is Rf=-1441(ρs)3+4616(ρs)2-4927(ρs)+1752, and correlation coefficient is 0.992.Figure 8
Effects of sludge density on floc diameter.


#### 4.1.4. Fractal Dimension

The fractal dimension was calculated with gyration radius method. the mean values of fractal dimension of the three series were shown in Figure9. The relation formulas of sludge density to gyration radius is Df=1.3338±0.0451+(0.6049±0.1424)×(ρs-1000)0.4107±0.4181, and correlation coefficient is 0.991. The gyration radius is increased with the increase of sludge density. The reason is that with the increase of sediment density, the collision probability of sediment particles and floc increased, while the opportunities for floc to accept new particles increased, which finally made the floc more compact. Several authors [21, 22] have researched the effects of sludge density on floc fractal dimension with numerical simulation method. The results show that relation formula of sludge density to gyration radius with Diffusion Limited Cluster Aggregation model is Df=Df0+aΦc, which is similar to our conclusion. Xiangyan et al. [15] has studied the effect of sludge density on gyration radius in two-dimensional space and the results presented in their paper are also similar to our conclusion.Figure 9
Effects of sludge density on fractal dimension.


### 4.2. Morphology of Floc

The morphology of the floc is a very important factor to determine the dynamic behavior of flocs. The morphology of a floc is visualized for the case of uncharged particles in Figure10(a). This floc, consisting of 300 primary particles, is obtained from a simulation based on 1000 primary particles with edge length decided by sludge density. Its structure is quite open to the outer space and it is very similar to the structures of aggregates observed by other researchers [21–23].Morphology of floc.
(a)
Floc formed by uncharged particles(b)
Floc formed by charged particles,q=-2.0×10-14C


### 4.3. Effects of Electrostatic Force on Flocs

Figure11 shows a comparison of the simple Coulomb force and the modified electrostatic force. The Coulomb force was calculated just within only one cell. There is no difference if two particles are close together, there is a relatively large difference as the distance becomes larger. Because of the periodic boundary condition, the electrostatic force must be zero if two particles are located at center and the edge of the simulation cell. The modified electrostatic force satisfies this condition well. Therefore, our model for the electrostatic force can be regarded as a good model coinciding with the real system. The morphology is quantified with a fractal dimension. It is affected by the number of primary particles and the radius of gyration, as shown in Figure 10(b). In this condition, the primary particles have −2.0 × 10−14C quantity of electricity according to (3.1).Figure 11
Comparison between modified electrostatic force versus simple Coulomb force.Four serieses of growth process are simulated to present the effects of charges on smashing time. As shown in Figure12, the effect of electric quantity on floc is obvious. However, the effects of unipolar and polar electrical charges on floc are still needed to be researched in the further work. Figures 13(a)–13(c) showed the effect of charges on smashing time. The smashing time is increased with sludge density as well as electrical charges.Figure 12
Variety of particle number as floc growing.Effects of sludge density on smashing time.
(a)
q=-2×10-14C(b)
q=-3×10-14C(c)
q=-4×10-14C


## 4.1. Floc Growth Process

### 4.1.1. Number of Particles in Floc

A series of comparison ofdp from 0.1 to 4.0 um were performed to interpret floc growth process. We chose ρs= 1010, 1025, 1050, 1075, and 1100 kg m−3 for simulation in order to intergrate with practice. As shown in Figure 5, in the early stage of flocculation the primary particle was sustaining to the a initial flocculent time. The particle number begun to increase after the first stage till reaching the smashing condition. In the final stage of flocculation, the sequence of smashing time with different sludge density is similar to that of initial flocculent time. The initial flocculent time and smashing time are both increased as the floc diameter increased.Variety of particle number as floc growing.
(a)
dp=0.1(b)
dp=1um(c)
dp=4um


### 4.1.2. The Smashing Time

The floc size can not grow up unlimited in the practice. It could be affect by particle diameter and particle-bonding mode. In this paper, the time that 300 particles conglutinated in one floc is smashing time. Though smashing time was related to random velocity and displacement, it increased with the decrease of sludge density. Whereas it increased with the increase of particle diameter as shown in Figures6 and 7. This tendency may be understood, on one hand as particle diameter presents a linear increase, the particle number decreased in a cube way, which leads particle impact probability to decrease. On the other hand, the decrease of sludge density causes particle impact probability to decrease likewise.Figure 6
Effects of sludge density on smashing time.Figure 7
Effects of particle diameter on smashing time.


### 4.1.3. Radius of Gyration

The relation between the radius of gyration and particle number was not linear positive correlation. The main factors influencing it were particle number and fractal dimension. In this simulation, the radius of gyration is calculated by four variables, they were particle number andx-, y-, z-axis coordinate values. The radius of gyration values in Figure 8 was calculated by three duplicate computations with the initial particle size being equal to 2 um. It was shown in the Figure 8 that the gyration radius decreased with the increase of sludge density. The relation formulas of sludge density to gyration radius is Rf=-1441(ρs)3+4616(ρs)2-4927(ρs)+1752, and correlation coefficient is 0.992.Figure 8
Effects of sludge density on floc diameter.


### 4.1.4. Fractal Dimension

The fractal dimension was calculated with gyration radius method. the mean values of fractal dimension of the three series were shown in Figure9. The relation formulas of sludge density to gyration radius is Df=1.3338±0.0451+(0.6049±0.1424)×(ρs-1000)0.4107±0.4181, and correlation coefficient is 0.991. The gyration radius is increased with the increase of sludge density. The reason is that with the increase of sediment density, the collision probability of sediment particles and floc increased, while the opportunities for floc to accept new particles increased, which finally made the floc more compact. Several authors [21, 22] have researched the effects of sludge density on floc fractal dimension with numerical simulation method. The results show that relation formula of sludge density to gyration radius with Diffusion Limited Cluster Aggregation model is Df=Df0+aΦc, which is similar to our conclusion. Xiangyan et al. [15] has studied the effect of sludge density on gyration radius in two-dimensional space and the results presented in their paper are also similar to our conclusion.Figure 9
Effects of sludge density on fractal dimension.


## 4.1.1. Number of Particles in Floc

A series of comparison ofdp from 0.1 to 4.0 um were performed to interpret floc growth process. We chose ρs= 1010, 1025, 1050, 1075, and 1100 kg m−3 for simulation in order to intergrate with practice. As shown in Figure 5, in the early stage of flocculation the primary particle was sustaining to the a initial flocculent time. The particle number begun to increase after the first stage till reaching the smashing condition. In the final stage of flocculation, the sequence of smashing time with different sludge density is similar to that of initial flocculent time. The initial flocculent time and smashing time are both increased as the floc diameter increased.Variety of particle number as floc growing.
(a)
dp=0.1(b)
dp=1um(c)
dp=4um


## 4.1.2. The Smashing Time

The floc size can not grow up unlimited in the practice. It could be affect by particle diameter and particle-bonding mode. In this paper, the time that 300 particles conglutinated in one floc is smashing time. Though smashing time was related to random velocity and displacement, it increased with the decrease of sludge density. Whereas it increased with the increase of particle diameter as shown in Figures6 and 7. This tendency may be understood, on one hand as particle diameter presents a linear increase, the particle number decreased in a cube way, which leads particle impact probability to decrease. On the other hand, the decrease of sludge density causes particle impact probability to decrease likewise.Figure 6
Effects of sludge density on smashing time.Figure 7
Effects of particle diameter on smashing time.


## 4.1.3. Radius of Gyration

The relation between the radius of gyration and particle number was not linear positive correlation. The main factors influencing it were particle number and fractal dimension. In this simulation, the radius of gyration is calculated by four variables, they were particle number andx-, y-, z-axis coordinate values. The radius of gyration values in Figure 8 was calculated by three duplicate computations with the initial particle size being equal to 2 um. It was shown in the Figure 8 that the gyration radius decreased with the increase of sludge density. The relation formulas of sludge density to gyration radius is Rf=-1441(ρs)3+4616(ρs)2-4927(ρs)+1752, and correlation coefficient is 0.992.Figure 8
Effects of sludge density on floc diameter.


## 4.1.4. Fractal Dimension

The fractal dimension was calculated with gyration radius method. the mean values of fractal dimension of the three series were shown in Figure9. The relation formulas of sludge density to gyration radius is Df=1.3338±0.0451+(0.6049±0.1424)×(ρs-1000)0.4107±0.4181, and correlation coefficient is 0.991. The gyration radius is increased with the increase of sludge density. The reason is that with the increase of sediment density, the collision probability of sediment particles and floc increased, while the opportunities for floc to accept new particles increased, which finally made the floc more compact. Several authors [21, 22] have researched the effects of sludge density on floc fractal dimension with numerical simulation method. The results show that relation formula of sludge density to gyration radius with Diffusion Limited Cluster Aggregation model is Df=Df0+aΦc, which is similar to our conclusion. Xiangyan et al. [15] has studied the effect of sludge density on gyration radius in two-dimensional space and the results presented in their paper are also similar to our conclusion.Figure 9
Effects of sludge density on fractal dimension.


## 4.2. Morphology of Floc

The morphology of the floc is a very important factor to determine the dynamic behavior of flocs. The morphology of a floc is visualized for the case of uncharged particles in Figure10(a). This floc, consisting of 300 primary particles, is obtained from a simulation based on 1000 primary particles with edge length decided by sludge density. Its structure is quite open to the outer space and it is very similar to the structures of aggregates observed by other researchers [21–23].Morphology of floc.
(a)
Floc formed by uncharged particles(b)
Floc formed by charged particles,q=-2.0×10-14C


## 4.3. Effects of Electrostatic Force on Flocs

Figure11 shows a comparison of the simple Coulomb force and the modified electrostatic force. The Coulomb force was calculated just within only one cell. There is no difference if two particles are close together, there is a relatively large difference as the distance becomes larger. Because of the periodic boundary condition, the electrostatic force must be zero if two particles are located at center and the edge of the simulation cell. The modified electrostatic force satisfies this condition well. Therefore, our model for the electrostatic force can be regarded as a good model coinciding with the real system. The morphology is quantified with a fractal dimension. It is affected by the number of primary particles and the radius of gyration, as shown in Figure 10(b). In this condition, the primary particles have −2.0 × 10−14C quantity of electricity according to (3.1).Figure 11
Comparison between modified electrostatic force versus simple Coulomb force.Four serieses of growth process are simulated to present the effects of charges on smashing time. As shown in Figure12, the effect of electric quantity on floc is obvious. However, the effects of unipolar and polar electrical charges on floc are still needed to be researched in the further work. Figures 13(a)–13(c) showed the effect of charges on smashing time. The smashing time is increased with sludge density as well as electrical charges.Figure 12
Variety of particle number as floc growing.Effects of sludge density on smashing time.
(a)
q=-2×10-14C(b)
q=-3×10-14C(c)
q=-4×10-14C


## 5. Conclusions

To investigate sludge drying process, a numerical simulation for the floc of uncharged and charged particles was presented with a Brownian dynamic. The initial conditions for the particle position were contained by sludge density and velocity was given with Gaussian distribution. To treat the electrostatic force, a modified electrostatic force has been introduced.In a fixed volume of water, the particle number increases gradually as the time and sludge density increase. In the meantime its increasing rate increases with time. However, it decreases as the increase of particle diameter. Though the smashing time was related to random velocity and displacement, it increased as the decrease of sludge density, yet it increased as the increase of initial particle diameter.The gyration radius is increased with the increase of sludge density. The relation formula of sludge density to fractal dimension isRf=-1441(ρs)3+4616(ρs)2-4927(ρs)+1752, and correlation coefficient(R2) is 0.991. The fractal dimension was calculated with gyration radius method. The relation formulas of sludge density to gyration radius is is Df=1.3338±0.0451+(0.6049±0.1424)×(ρs-1000)0.4107±0.4181, and correlation coefficient is 0.990. The conclusion is similar to the results of previous studies.Effects of electric quantity on floc were obvious. The smashing time was increased with the increase of electric quantity. Besides, the polar electrical charges on primal particles make the morphology of floc more open. However, this simulation was carried out under the assumption of the high-moisture content. Therefore, our conclusions cannot be applied directly to other cases, where the moisture content is extremely low. To find out the particles motiom law in the low moisture content, it needs some additional models, and we will perform additional simulations and experiments to confirm it as future works.



---

*Source: 527384-2012-03-22.xml*