Next Article in Journal
MHD Free Convection and Entropy Generation in a Corrugated Cavity Filled with a Porous Medium Saturated with Nanofluids
Previous Article in Journal
Generalized Distance-Based Entropy and Dimension Root Entropy for Simplified Neutrosophic Sets
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Simple Explicit Expression for the Flocculation Dynamics Modeling of Cohesive Sediment Based on Entropy Considerations

College of Water Sciences, Beijing Normal University, Xinjiekouwai Street 19, Beijing 100875, China
Entropy 2018, 20(11), 845; https://doi.org/10.3390/e20110845
Submission received: 13 September 2018 / Revised: 25 October 2018 / Accepted: 2 November 2018 / Published: 4 November 2018

Abstract

:
The flocculation of cohesive sediment plays an important role in affecting morphological changes to coastal areas, to dredging operations in navigational canals, to sediment siltation in reservoirs and lakes, and to the variation of water quality in estuarine waters. Many studies have been conducted recently to formulate a turbulence-induced flocculation model (described by a characteristic floc size with respect to flocculation time) of cohesive sediment by virtue of theoretical analysis, numerical modeling, and/or experimental observation. However, a probability study to formulate the flocculation model is still lacking in the literature. The present study, therefore, aims to derive an explicit expression for the flocculation of cohesive sediment in a turbulent fluid environment based on two common entropy theories: Shannon entropy and Tsallis entropy. This study derives an explicit expression for the characteristic floc size, assumed to be a random variable, as a function of flocculation time by maximizing the entropy function subject to the constraint equation using a hypothesis regarding the cumulative distribution function of floc size. It was found that both the Shannon entropy and the Tsallis entropy theories lead to the same expression. Furthermore, the derived expression was tested with experimental data from the literature and the results were compared with those of existing deterministic models, showing that it has good agreement with the experimental data and that it has a better prediction accuracy for the logarithmic growth pattern of data in comparison to the other models, whereas, for the sigmoid growth pattern of experimental data, the model of Keyvani and Strom or Son and Hsu model could be the better choice for floc size prediction. Finally, the maximum capacity of floc size growth, a key parameter incorporated into this expression, was found to exhibit an empirical power relationship with the flow shear rate.

1. Introduction

Cohesive sediment, which is different from non-cohesive sediments such as sand, gravel, and cobbles, is a mixture of water, fine-grained sediments, such as silt, clay, and organic matter of diverse natures [1,2]. When cohesive sediment particles are transported in rivers, reservoirs, lakes, estuarines, and coastal waters, they continually flocculate to form flocs of different sizes due to small-scale, particle-particle interactions. On the other hand, some fragile and loose flocs may break into small flocs and/or primary particles (floc breakage or floc disaggregation) due to the flow shear [3,4,5]. Flocs are totally different from primary sediment particles in terms of their larger sizes, lower excess density, and higher settling velocity in water [2,6]. Studying cohesive sediment flocculation in a turbulent flow environment is essential because it plays an important role in affecting the morphological changes to coastal areas, dredging operations in navigational canals, and sediment siltation in reservoirs and lakes [7,8]. Since some pollutants (such as heavy metals) and nutrients are absorbed on the surfaces of cohesive sediment particles due to the electrochemical attraction of clay particles and/or organic matter contained in the sediment, the flocculation of cohesive sediment is also a vital element in investigating the variation of water quality and ecosystem function in some waters such as lakes and estuarine and coastal waters, which contain an abundance of cohesive sediment [9,10].
The turbulence-induced flocculation of cohesive sediment and other particles (such as polystyrene/latex particles) in a fluid environment has been investigated by several researchers in many research fields, including chemical and environment engineering, oceanography, and river and estuarine mechanics (e.g., References [11,12,13,14,15,16,17,18,19,20,21,22,23]). Most studies regarding turbulence-induced particle flocculation have focused on two main aspects: (1) the floc properties (mainly characterized by floc size or floc structure) at the steady or equilibrium states, in which the property parameters reach constant values; and (2) temporal variations of the size distribution and the structural and morphological properties of the flocs (commonly characterized by different fractal dimensions of the flocs) during the flocculation/aggregation process.
Some experimental works have been performed to investigate the median value of the size distribution of the flocs at the steady state of flocculation with respect to various flow shear conditions (e.g., References [8,14,24,25,26,27]. These studies reported that the median floc size decreases as the flow shear stress increases. Furthermore, a power relationship function was commonly adopted to describe this dependence: the median size = c * (flow shear parameter) γ , where c and γ are two positive constants. The c is the floc strength, which strongly depends on the method used to measure the floc size, while the γ is the stable floc size exponent depending on the breakage mechanisms (erosion or fracture) for flocs smaller or larger than the smallest eddy (i.e., Kolmogorov microscale) in the turbulent flow [28,29]. Some studies have focused on the structural and morphological properties of the flocs at the steady state of flocculation with respect to various flow shear conditions (e.g., References [6,17,25,30]). For example, Stone and Krishnappan [30] showed that particle boundaries become more convoluted and the shape of larger particles are more irregular at higher levels of flow shear stress, whereas Zhu et al. [6] reported that with increasing flow shear rates, the flocs become less elongated and their boundary lines become tighter and more regular.
The time evolution of the size distribution of the flocs during the turbulence-induced flocculation has been investigated by some researchers using experimental observational techniques or numerical modeling methods (e.g., References [14,16,24,31,32]). Some experimental works have reported a typical trend of floc size with respect to flocculation time: the median value of the size distribution of the floc population grows rapidly with time at the beginning of the flocculation experiment. This is because the flow shear increases the collision and adhesion between primary particles, producing some large flocs in the system [24,25,31]. However, as flocculation progresses further, the rapidity with which the median floc size increases with time begins to decline (that is, the floc size experiences a slowly increasing process with flocculation time) because large flocs possess a fragile and loose structure and are susceptible to breakage induced by the flow shear [12,14]. Finally, the median floc size reaches a steady or equilibrium state after a long flocculation time as a result of a dynamic balance between the floc growth and the floc breakage, which are both caused by the flow shear [16,33]. On the other hand, some authors have developed theoretical and/or numerical models to describe the temporal evolution of the size distribution of a floc population, among which a simplified Lagrangian flocculation model is the focus of some works [1,3,34,35]. The earliest form of this model was developed by Winterwerp [3] with a linear combination of the formulations for the floc aggregation and the floc breakage process due to the flow turbulence, with a constant fractal dimension of flocs during flocculation. This model primarily tracks the time evolution of a characteristic floc size (commonly the median value of the size distribution of flocs) during flocculation, and its prediction accuracy is further improved by some authors by virtue of modifying some key parameters that have been incorporated into the model [1,34]. Additionally, there have also been some experimental studies to demonstrate the temporal variations of structural and geometric properties of flocs during the turbulence-induced flocculation process (e.g., References [17,25,36]).
This study focuses on the temporal evolution of the characteristic floc size during flocculation. To the best of my knowledge, most of the studies regarding turbulence-induced flocculation have been performed theoretically, with deterministic approaches. However, a probability approach to investigate the flocculation process is still lacking in the literature. For the last two decades, the probabilistic treatment of hydraulic engineering problems based on entropy theory has gained the attention of some researchers [37]. For example, entropy-based studies have been conducted by many researchers to investigate the velocity distribution (e.g., References [38,39,40,41]), sediment concentration [42,43,44], and shear stress distribution [45,46] in open channels. Recent work on the application of the entropy theory in a classical hydraulic engineering problem can be found in the study of Singh et al. [37]. It should be noted that the entropy-based probability method might also be an easy and applicable tool for predictions in some hydraulic problems, even though the traditional deterministic approaches have provided some physical insights.
This study attempts to derive an entropy-based expression for the temporal evolution of the characteristic floc size during flocculation by using two known entropy theories: Shannon entropy and Tsallis entropy. Section 2 and Section 3 derive the characteristic floc size as a function of flocculation time using these two entropy theories. The derived expression is tested against experimental observation data in Section 4, and Section 5 contains a comparison of the expression with some developed deterministic models, as well as a discussion of the maximum capacity of floc size growth, a key parameter that has been incorporated into the expression. Finally, Section 6 presents the concluding remarks.

2. Shannon Entropy Theory for Flocculation Expression

The determination of the flocculation process of cohesive sediment using the Shannon entropy theory entails the following steps: (1) the definition of the Shannon entropy; (2) the specification of constraints; (3) the maximization of entropy; (4) the determination of the Lagrange multiplier; (5) the hypothesis regarding cumulative probability distribution; and (6) the derivation of the flocculation expression.

2.1. Definition of Shannon Entropy

Considering entropy as a measure of information and, therefore, of uncertainty, Shannon [47] formulated what is referred to as the Shannon entropy theory. The Shannon entropy quantitatively measures the mean uncertainty associated with a probability distribution of a random variable. Consider the characteristic floc size D , during the flocculation process as a continuous random variable [48,49]. The objective of this study is to derive the characteristic floc size as a function of flocculation time, t . For simplicity, we let the excess floc size, D ^ ( t ) , be defined as D ^ ( t ) = D D ( t ) , where D is the steady or equilibrium state value of the floc size. Therefore, at the beginning of flocculation ( t = 0 ), there is D ^ ( 0 ) = D D 0 , where D 0 is the median size of primary particles, whereas after a long time of flocculation ( t ), there is a relation: D ^ ( t ) = 0 . Thus, the excess floc size D ^ ( t ) will vary from 0 to D ^ ( 0 ) . For the probability density function of the excess floc size D ^ ( t ) , f ( D ^ ) , the Shannon entropy, denoted by H S ( D ^ ) , can be expressed in the general form as
H S ( D ^ ) = 0 D D 0 f ( D ^ ) ( ln f ( D ^ ) ) d D ^  
Theoretically, the Shannon entropy is at a maximum when the probability density function is uniform within its limits. Equation (1) expresses a measure of uncertainty of f ( D ^ ) or the average information content of sample D ^ .

2.2. Specification of Constraint

The total probability law must be satisfied for the probability density function f ( D ^ ) . Therefore, the constraint equation can be written as
0 D D 0 f ( D ^ ) d D ^ = 1  

2.3. Maximization of Entropy

To derive the specific form of f ( D ^ ) , we adopted the principle of maximum entropy developed by Jaynes [50,51,52]. This principle states that the least biased probability of D ^ , f ( D ^ ) , will be the one that will maximize H ( D ^ ) given by Equation (1), subject to the given information on D ^ expressed as a constraint equation. Such a probability distribution is yielded by the maximization of the Shannon entropy. To that end, the method of the Euler–Lagrange calculus of variation is used [37]. The Lagrangian function L can be written as follows:
L = f ( D ^ ) ln f ( D ^ ) + ( 1 λ 0 ) f ( D ^ )  
where λ 0 is the zeroth Lagrange multiplier.
Differentiating Equation (3) with respect to f ( D ^ ) and equating the derivative to zero, the probability density function f ( D ^ ) of the excess floc size D ^ is given as
f ( D ^ ) = exp ( λ 0 )  
Therefore, in the cumulative distribution function (CDF), f ( D ^ ) of D ^ is obtained by using Equation (4) as follows:
F ( D ^ ) = P ( d D ^ ) = exp ( λ 0 ) D ^  
Both the probability density function and the cumulative distribution function depend on the value of the zeroth Lagrange multiplier λ 0 .

2.4. Determination of the Lagrange Multiplier

Inserting Equation (4) into the constraint equation (Equation (2)) leads to the following relation:
0 D D 0 exp ( λ 0 ) d D ^ = 1 f ( D ^ ) = 1 D D 0  
The combination of Equations (4) and (6) gives the following:
λ 0 = ln ( D D 0 )  
The value of the Lagrange multiplier λ 0 can be obtained as long as the values of D 0 and D are known from the observational data.

2.5. Hypothesis on the Cumulative Distribution Function

To derive the temporal evolution of the excess floc size D ^ in the real (space) domain, an equation connecting the probability domain to the space domain is required [37]; therefore, a hypothesis on the CDF of the excess floc size D ^ is made so that the hypothesized CDF can reflect the characteristic of D ^ .
Consider a simple flocculation element as shown in Figure 1. At the beginning of flocculation, some primary particles collide and adhere in small flocs due to the eddy motion of the turbulent flow, and the floc size D ( t ) increases significantly [18,19,25]. Whereas after a certain flocculation time, those formed fragile and loose flocs easily undergo a breakup due to the flow shear; therefore, the floc size growth begins to decline [12,14]. Let the floc size exiting the flocculation element be denoted as D which approximately equals the steady state of the floc size. The flocculation element will have a maximum capacity of floc size growth, denoted by S (its unit should be in m*s). If we define the cumulative floc size growth as J (its unit should also be in m*s), then 0 J S for the flocculation element. The continuity equation for the flocculation element, as shown in Figure 1, can be expressed as
d J d t = D D ( t ) ;   or   J ( t ) = D t 0 t D ( t ) d t  
It is hypothesized that the cumulative distribution function F ( D ) of the floc size can be defined as the ratio of the cumulative floc size growth to the maximum capacity of the floc size growth or maximum potential floc size growth, S :
F ( D ) = J S  
Here S has the same units as J . In Equation (9), it is implied that all of the values of the cumulative floc size growth are equally likely. A similar hypothesis has been employed by Chiu [53] and Kumbhakar and Ghoshal [41] for deriving a one-dimensional velocity distribution in open channels, by Chiu et al. [42] and Kumbhakar et al. [44] for deriving the sediment concentration profiles, and by Khozani and Bonakdari [54] for deriving the shear stress distribution in open channels. As Singh [55] showed, even if the above hypothesis is not strictly valid, it will not greatly influence the results because it merely allows the entropy theory to lead to the equation for floc size growth that is desired.
The differentiation of Equation (9) gives
d F ( D ) d D = f ( D ) = 1 S d J d D  

2.6. Derivation of the Flocculation Process

Combining Equations (6), (8), and (10) yields
1 D D 0 d D d t = 1 S ( D D )  
Integrating Equation (11) and using the initial condition: D = D 0 at t = 0 , we obtain
D ( t ) = D ( D D 0 ) exp ( D D 0 S t )  
Substituting Equation (12) into Equation (8), the cumulative floc size growth J is written as
J = S [ 1 exp ( D D 0 S t ) ]  
Inserting Equation (13) in Equation (9) yields the cumulative distribution function F ( D ) of the floc size as
F ( D ) = exp ( D D 0 S t )  
Finally, the entropy of the probability distribution of the floc size H S ( D ) is obtained by substituting Equation (6) into Equation (1) as follows:
H S ( D ) = ln ( D D 0 )  
Equation (15) states that the uncertainty of the floc size depends on the initial floc size value D 0 and the steady state value D .

3. Tsallis Entropy Theory for the Flocculation Model

The application of the Tsallis entropy theory into the derivation of the floc size as a function of flocculation time D ( t ) contains the same procedure as the Shannon entropy.
If the floc size D ( t ) is considered as a continuous random variable with a probability function defined as f ( D ) , another entropy function that has been termed as Tsallis entropy, H T ( D ) , which was proposed by Tsallis [56] as a generalized form of the Shannon entropy, can be written as follows:
H T ( D ) = 1 m 1 { 1 D 0 D [ f ( D ) ] m d D }  
where m is a real number not equal to 1. The Tsallis entropy is a non-extensive entropy that reduces to the Shannon entropy if the exponent m 1 in Equation (16). For any m , it takes its maximum value in the case of equiprobability, and this entropy function reaches its maximum (concave function) if m 0 and its minimum (convex function) if m 0 for a certain value of m [57].
Similar to the Shannon entropy method, the constraint equation that f ( D ) must satisfy becomes
D 0 D f ( D ) d D = 1  
Using the principle of maximum entropy, the Lagrangian function L for the Tsallis entropy can be written as follows:
L = 1 m 1 { 1 D 0 D [ f ( D ) ] m d D } + λ 0 [ D 0 D f ( D ) d D 1 ]  
where λ 0 is the zeroth Lagrange multiplier. Differentiating Equation (18) with respect to f ( D ) and equating the derivative to zero, the probability density function f ( D ) of the floc size is obtained as
f ( D ) = [ m 1 m ( 1 m 1 + λ 0 ) ] 1 m 1  
Substituting Equation (19) into Equation (17), we get
f ( D ) = 1 D D 0 ,   and   λ 0   =   m m 1 ( 1 D D 0 ) m 1 1 m 1  
Combining Equations (8), (10), and (20) yields: 1 / ( D D 0 ) d D / d t = ( D D ) / S , which is Equation (11). In the same way as the Shannon entropy method, by integrating this equation and using the initial condition D = D 0 at t = 0 , we obtain the function of D ( t ) as Equation (12). Similarly, the cumulative distribution function F ( D ) is also derived as Equation (14).
Inserting Equation (20) into Equation (16), we obtain the Tsallis entropy of the probability density function of the floc size as
H T ( D ) = 1 m 1 [ ( D D 0 ) ( D D 0 ) 1 m ] ,
which depends on three parameters: the initial floc size value D 0 , the steady state value D , and the parameter m .
It can be seen that both the Shannon entropy and the Tsallis entropy produce the same analytical expression (Equation (12)) that describes the temporal evolution of floc size during the turbulence-induced flocculation process. We also need to point out that the proposed model (Equation (12)) refers to a monodisperse distribution system, as the heterodisperse characteristic of aggregates were not considered at the presented modeling.

4. Results

Thirty-three experimental data sets regarding the floc size with respect to flocculation time in the published literature were collected to test the validity of the entropy-based expression (Equation (12)) in this study. Table 1 presents the information on these collected experimental data. The first column number is the experimental data. The second column introduces the particle material: some adopted the sediment material, whereas some used the polystyrene/latex material, and the third column presents the apparatus for generating the turbulent flocculation environment. In the fourth column, ϕ is the particle volumetric concentration (it is equal to the volume of the primary particle divided by the volume of the particle-liquid mixture). G is the flow shear rate (its unit is 1/s), defined as ε / ν , where ε is the turbulent dissipation rate of the turbulent flow and ν is the kinematic viscosity of the fluid, as adopted by many studies [12,14,16,25,29]. The measured size of the primary particles and the floc size at the steady state of flocculation are shown in the fifth and sixth columns, respectively, and the data source is identified in the last column. The criteria used to take experiments from the literature for modeling validation are that the selected data sets cover different flocculation materials (sediment or polystyrene/latex particle), different flocculation environments (Couette-flow system or baffled stirred tank), and various flow shear conditions (the low flow turbulent condition, for example, G = 0.45, 0.75, 0.96, 2.4 s 1 ; the moderate flow turbulent condition, for example, G = 19.4, 25, 37, 50 s 1 ; the strong turbulent condition, for example, G = 100, 150, 20, 246 s 1 ).
To evaluate the performance of the derived entropy-based flocculation expression with experimental observation data and some deterministic models, an error analysis is performed by computing the correlation coefficient R 2 between the modeled and the observed data, the relative bias (RBIAS) between the modeled and the observed data, defined as RBIAS = 1 N i = 1 N | m i o i o i | , and the root-mean-square error (RMSE), defined as RMSE = 1 N i = 1 N ( m i o i ) 2 , where m and o are the modeled and observed points, respectively, and N is the number of observed points. The goodness of fit increases as the R 2 value increases and both the RBIAS and RMSE values decrease.
Figure 2 shows the comparison of the proposed entropy-based model with the collected experimental data. Table 2 presents the comparison results. From the third, fourth, and fifth columns, it can be seen that there is a high R 2 value and low RBIAS and RMSE values for each case. Additionally, the entropy function values estimated by Equations (15) and (21) are also presented in the last columns. These results indicate that the proposed entropy-based model shows a good agreement with the experimental data.

5. Discussion

5.1. Comparison with the Deterministic Model

To further test the prediction accuracy of the proposed entropy-based expression, we compare it with some deterministic models. There have been three main simplified Lagrangian flocculation models: Winterwerp [3], Son and Hsu [1], and Son and Hsu [34]. Table 3 list these models. c and ρ s are the mass concentration and the density of primary particles, respectively; k A is the dimensionless coefficient for floc aggregation; k B is the dimensionless coefficient for floc breakage; d f and F are the fractal dimension and yield strength of the flocs, respectively; μ is the dynamic viscosity of the fluid; α and β are two coefficients; and B is a coefficient representing the cohesive force between the primary particles.
In the study of Son and Hsu [34], a comparison among the Winterwerp model, the Son and Hsu [1] model and the Son and Hsu [34] model were conducted with the experimental results of Burban et al. [58] and an experimental data set from Biggs and Lant [14]. To simplify the problem, we attempted to compare the proposed entropy-based model (Equation (12)) with the Winterwerp model, the Son and Hsu (2008) model, and the Son and Hsu (2009) model for the experimental data in this study. Table 4 presents the calculated R 2 , RBIAS, and RMSE values for these real cases. It can be observed that the proposed model has the highest R 2 value and the lowest RBIAS and RMSE values in comparison with the other three models for all of the real cases. For the case of ϕ = 1.66 × 10−3 from Burban et al. [58], the proposed model has the highest R 2 value and the lowest RBIAS value compared with the other models, whereas the model of Son and Hsu (2008) yields a lower RMSE value than the proposed model, which may be because of the very limited experimental data. Hence, this study shows the potential of the Shannon entropy together with the principle of maximum entropy to predict the temporal evolution of floc size during flocculation.
In the work of Keyvani and Strom [33], the effects of seven cycles of high and low turbulent shear on mud floc growth pattern and equilibrium size were investigated through a laboratory study. The measured temporal variation of the mean floc size in each of the seven cycles was modelled using the Winterwerp model after calibrating the collision and breakup efficiency coefficients for each cycle in their paper. Here we attempted to compare the entropy-based expression (Equation (12)) and the model of Keyvani and Strom for these measured data, as shown in Figure 3. Note that the horizontal axis in these figures refers to the logarithmic coordinates. Most of the measured data show a sigmoid growth pattern of floc size with flocculation time. For the sigmoid growth curve of floc size with flocculation time, it could be found from Figure 3 that the entropy-based expression (Equation (12)) did not exhibit a satisfactory fitting result for these measured data (especially in Figure 3d,e), however, the model of Keyvani and Strom presents a typical sigmoid growth property and shows a better prediction accuracy in comparison to Equation (12). In Figure 2, most of the measured data show a logarithmic growth pattern of floc size with flocculation time, and the entropy-based expression has provided a good fitting result for them. This is because this expression is based on the assumption that the cumulative distribution of floc size can be defined as the ratio of the cumulative floc size growth to the maximum potential floc size growth, and this leads to the logarithmic result. It could be concluded that the entropy-based expression developed in this study can fit well for the logarithmic growth pattern of floc size, whereas for the sigmoid growth pattern of floc size, the model of Keyvani and the Strom or Son and Hsu (2009) model could be the better choice for floc size prediction. This limitation of the developed entropy-based expression could be worthy of further investigation in future research.

5.2. Estimation of the Key Parameter

The key parameter that was incorporated into the expression (Equation (12)) is the capacity for floc size growth in the flocculation system S . Here, we compare the fitted values of S with the flow shear rate G for all experimental data (T1–T33), as shown in Figure 4, except for T1, T2, T7–T10, and T23–T25. The reason for these exceptions is that both studies of Burban et al. [58] (that is, T1–T2) and Spicer and Pratsinis [59] (that is, T7–T10) adopted the fixed G values, and the value of the absolute temperature T is not available in the study of Selomulya et al. [60] (that is, T23–T25); therefore, the calculated S does not have the same unit as the other real cases. It should be noted that in the study of Stone and Krishnappan [30] (that is T26–T28), only the bed shear stress τ b in the turbulent-generating equipment is provided, and the flow shear rate G cannot be calculated; thus, Figure 5 shows the calculated S values with respect to different bed shear stresses for T26–T28.
It can be seen from Figure 4 that there is a fitting relation between S and G as follows: S = 10 6 G 0.844 , with a very high coefficient of determination R 2 reaching 0.93. From Figure 5, there is a similar fitting relation between S and τ b as follows: S = 9426.8 τ b 1.517 , with a coefficient of determination R 2 reaching 0.9937. This implies that as the flow shear condition intensifies, the capacity for floc size growth in the flocculation system decreases. This is because the floc breakage caused by the increasing flow shear plays an increasingly important role in the flocculation process.
Substituting the aforementioned mathematical relation into Equation (12) leads to the final expression for floc size D as a function of flocculation time t as follows:
D ( t ) = D ( D D 0 ) exp [ 10 - 6 ( D D 0 ) G 0.844 t ]  
for the collected experimental data, except for Stone and Krishnappan [30], whereas for the experimental data of Stone and Krishnappan [30], the expression has the following form: D ( t ) = D ( D D 0 ) exp [ 10 - 4 ( D D 0 ) τ b 1.517 t ] .
We estimated the floc size D using Equation (22) and the aforementioned other expression, and compared them with the collected experimental data, as shown in Figure 6. It can be observed that there is a very high coefficient of determination R 2 between the estimated values and the observed ones for most data sets. Only four data sets (T3, T11, T12, and T16) have an R 2 value smaller than 0.90, and the reason may be that there is a data scattering perhaps due to the experimental measurement operation. These results imply that the proposed entropy-based expression has a good prediction ability for the temporal evolution of floc size during the turbulence-induced flocculation.
Similar to the Camp number P (this is equal to the product of flow shear rate and flocculation time, that is, P = G t ) defined in the research field of wastewater treatment [28], we define a new parameter P as P = G 0.844 t , and Equation (22) becomes D ( t ) = D ( D D 0 ) exp [ 10 - 6 ( D D 0 ) G 0.844 t ] = D ( D D 0 ) exp [ 10 - 6 ( D D 0 ) P ] . The floc size is a monotonic increasing function of the new parameter P : it firstly undergoes a rapidly increasing process, and then a slowly increasing process before approaching a final steady state. When D 0 D , the floc size approximately reaches the steady state if P is larger than 3 × 10 6 / D . It is easy and applicable to adopt Equation (22) to predict the floc size during flocculation as long as the values of D 0 , D , and G are given. Equation (22) provides a new method for flocculation dynamic modeling based on entropy considerations. It has a simple mathematical form, contains fewer parameter inputs compared with other existing deterministic models, and avoids an iteration calculation required in other models for the floc size estimation. This equation contains the effect of the flow shear on the floc breakage (as shown by Figure 4). However, some physical properties present in existing deterministic models are not incorporated into this equation. For example, the geometric structure of floc plays a role in both dynamic processes of floc growth and floc breakage. In the existing deterministic models, the fractal dimension of floc d f has been adopted to describe this property. However, the entropy-based expression does not contain this parameter.

6. Concluding Remarks

The following concluding remarks can be made from this study:
  • A simple explicit expression that describes the temporal evolution of the characteristic floc size during turbulence-induced flocculation was derived based on the entropy theory.
  • Both the Shannon entropy theory and the Tsallis entropy theory lead to the same expression for the function of floc size with respect to flocculation time.
  • The entropy-based expression was tested against the experimental data in the literature, and a good agreement was found.
  • The entropy-based expression was compared with other deterministic models, and it was found that the expression shows a better prediction accuracy for the logarithmic growth pattern of experimental data in comparison to the other models, whereas, for the sigmoid growth pattern of data, the model of Keyvani and Strom or the Son and Hsu (2009) model could be the better choice for floc size prediction.
  • The maximum capacity of floc size growth, a key parameter that was incorporated into the expression, exhibits an empirical power-law relation with the flow shear rate. As the flow shear condition intensifies, the capacity for floc size growth in the flocculation system decreases. This is because the floc breakage caused by the increasing flow shear plays an increasingly important role in the flocculation process.

Funding

This research was funded by the National Natural Science Foundation of China (51509004).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Son, M.; Hsu, T.J. Flocculation model of cohesive sediment using variable fractal dimension. Environ. Fluid Mech. 2008, 8, 55–71. [Google Scholar] [CrossRef]
  2. Pejrup, M.; Mikkelsen, O.A. Factors controlling the field settling velocity of cohesive sediment in estuaries. Estuar. Coast Shelf Sci. 2010, 87, 177–185. [Google Scholar] [CrossRef]
  3. Winterwerp, J.C. A simple model for turbulence induced flocculation of cohesive sediment. J. Hydraul. Res. 1998, 36, 309–326. [Google Scholar] [CrossRef]
  4. Xu, F.; Wang, D.P.; Riemer, N. Modeling flocculation processes of fine-grained particles using a size-resolved method: Comparison with published laboratory experiments. Cont. Shelf Res. 2008, 28, 2668–2677. [Google Scholar] [CrossRef]
  5. Dyer, K. Sediment processes in estuaries: Future research requirements. J. Geophys. Res. Oceans 1978–2012 1989, 94, 14327–14339. [Google Scholar] [CrossRef]
  6. Zhu, Z.; Yu, J.; Wang, H.; Dou, J.; Wang, C. Fractal dimension of cohesive sediment flocs at steady state under seven shear flow conditions. Water 2015, 7, 4385–4408. [Google Scholar] [CrossRef]
  7. Maggi, F. The settling velocity of mineral, biomineral, and biological particles and aggregates in water. J. Geophys. Res. Ocean 2013, 118, 2118–2132. [Google Scholar] [CrossRef] [Green Version]
  8. Guo, C.; He, Q.; van prooijen, B.C.; Guo, L.C.; Manning, A.J.; Bass, S. Investigation of flocculation dynamics under changing hydrodynamic forcing on an intertidal mudflat. Mar. Geol. 2018, 395, 120–132. [Google Scholar] [CrossRef]
  9. Maggi, F. Biological flocculation of suspended particles in nutrient-rich aqueous ecosystems. J. Hydrol. 2009, 376, 116–125. [Google Scholar] [CrossRef]
  10. Fang, H.W.; Lai, H.J.; Cheng, W.; Huang, L.; He, G.J. Modeling sediment transport with an integrated view of the biofilm effects. Water Resour. Res. 2017, 53, 7536–7557. [Google Scholar] [CrossRef]
  11. Parker, D.; Kaufman, W.; Jenkins, D. Floc breakup in turbulent flocculation processes. J. Sanit. Eng. Div. 1972, 98, 79–99. [Google Scholar]
  12. Serra, T.; Colomer, J.; Casamitjana, X. Aggregation and breakup of particles in a shear flow. J. Colloid Interface Sci. 1997, 187, 466–473. [Google Scholar] [CrossRef]
  13. Lu, S.; Ding, Y.; Guo, J. Kinetics of fine particle aggregation in turbulence. Adv. Colloid Inter. Sci. 1998, 78, 197–235. [Google Scholar] [CrossRef]
  14. Biggs, C.; Lant, P. Activated sludge flocculation: On-line determination of floc size and the effect of shear. Water Res. 2000, 34, 2542–2550. [Google Scholar] [CrossRef]
  15. van Leussen, W. Macroflocs, fine-grained sediment transport, and their longitudinal variations. Ocean Dyn. 2011, 61, 387–401. [Google Scholar] [CrossRef]
  16. Shen, X.; Maa, J.P.Y. Modeling floc size distribution of suspended cohesive sediments using quadrature method of moments. Mar. Geol. 2015, 359, 106–119. [Google Scholar] [CrossRef]
  17. Moruzzi, R.B.; de Oliveira, A.L.; da Conceicao, F.T.; Gregory, J.; Campos, L.C. Fractal dimensions of large aggregates under different flocculation conditions. Sci. Total Environ. 2017, 609, 807–814. [Google Scholar] [CrossRef] [PubMed]
  18. Winterwerp, J.C.; Manning, A.J.; Martens, C.; de Mulder, T.; Vanlede, J. A heuristic formula for turbulence-induced flocculation of cohesive sediment. Estuar. Coast Shelf Sci. 2006, 68, 195–207. [Google Scholar] [CrossRef]
  19. Xu, F.; Wang, D.; Riemer, N. An idealized model study of flocculation on sediment trapping in an estuarine turbidity maximum. Cont. Shelf Res. 2010, 30, 1314–1323. [Google Scholar] [CrossRef]
  20. Verney, R.; Lafite, R.; Brun-cottan, J.C.; Le, P. Behaviour of a floc population during a tidal cycle: Laboratory experiments and numerical modelling. Cont. Shelf Res. 2011, 31, S64–S83. [Google Scholar] [CrossRef]
  21. Strom, K.; Keyvani, A. Flocculation in a decaying shear field and its implications for mud removal in near-field river mouth discharges. J. Geophys. Res. Oceans 2016, 121, 2142–2162. [Google Scholar] [CrossRef]
  22. Sherwood, C.R.; Aretxabaleta, A.L.; Harris, C.K.; Rinehimer, J.P.; Verney, R.; Ferre, B. Cohesive and mixed sediment in the regional ocean modeling system (ROMS v3.6) implemented in the coupled ocean atmosphere wave sediment-transport modeling system (COAWST r1179). Geosci. Model Develop. 2018, 11, 1849–1871. [Google Scholar] [CrossRef]
  23. Kuprenas, R.; Tran, D.; Strom, K. A shear-limited flocculation model for dynamically predicting average floc size. J. Geophys. Res. Oceans 2018, 123, 6736–6752. [Google Scholar] [CrossRef]
  24. Oles, V. Shear-induced aggregation and breakup of polystyrene latex particles. J. Colloid Interface Sci. 1992, 154, 351–358. [Google Scholar] [CrossRef]
  25. Bubakova, P.; Pivokonsky, M.; Filip, P. Effect of shear rate on aggregate size and structure in the process of aggregation and at steady state. Powder Technol. 2013, 235, 540–549. [Google Scholar] [CrossRef]
  26. Hopkins, D.C.; Ducoste, J.J. Characterizing flocculation under heterogeneous turbulence. J. Colloid Interface Sci. 2003, 264, 184–194. [Google Scholar] [CrossRef]
  27. Manning, A.; Dyer, K. A laboratory examination of floc characteristics with regard to turbulent shearing. Mar. Geol. 1999, 160, 147–170. [Google Scholar] [CrossRef]
  28. Thomas, D.; Judd, S.; Fawcett, N. Flocculation modelling: a review. Water Res. 1999, 33, 1579–1592. [Google Scholar] [CrossRef]
  29. Zhu, Z. Theory on orthokinetic flocculation of cohesive sediment: a review. J. Geosci. Environ. Prot. 2014, 2, 24–31. [Google Scholar] [CrossRef]
  30. Stone, M.; Krishnappan, B. Floc morphology and size distributions of cohesive sediment in steady-state flow. Water Res. 2003, 37, 2739–2747. [Google Scholar] [CrossRef]
  31. Serra, T.; Casamitjana, X. Structure of the aggregates during the process of aggregation and breakup under a shear flow. J. Colloid Interface Sci. 1998, 206, 505–511. [Google Scholar] [CrossRef] [PubMed]
  32. Maggi, F.; Mietta, F.; Winterwerp, J.C. Effect of variable fractal dimension on the floc size distribution of suspended cohesive sediment. J. Hydrol. 2007, 343, 43–55. [Google Scholar] [CrossRef]
  33. Keyvani, A.; Strom, K. Influence of cycles of high and low turbulent shear on the growth rate and equilibrium size of mud flocs. Mar. Geol. 2014, 354, 1–14. [Google Scholar] [CrossRef]
  34. Son, M.; Hsu, T.J. The effect of variable yield strength and variable fractal dimension on flocculation of cohesive sediment. Water Res. 2009, 43, 3582–3592. [Google Scholar] [CrossRef] [PubMed]
  35. Xu, C.; Dong, P. A dynamic model for coastal mud flocs with distributed fractal dimension. J. Coast. Res. 2017, 33, 218–225. [Google Scholar] [CrossRef]
  36. Guerin, L.; Saudejaud, C.C.; Line, A.; Frances, C. Dynamics of aggregate size and shape properties under sequenced flocculation in a turbulent Taylor-Couette reactor. J. Colloid Interface Sci. 2017, 491, 167–178. [Google Scholar] [CrossRef] [PubMed]
  37. Singh, V.P.; Sivakumar, B.; Cui, H.J. Tsallis entropy theory for modelling in water engineering: A review. Entropy 2017, 19, 641–666. [Google Scholar] [CrossRef]
  38. Chiu, C.L.; Said, C.A.A. Maximum and mean velocities and entropy in open-channel flow. J. Hydraul. Eng. 1995, 121, 26–35. [Google Scholar] [CrossRef]
  39. Cui, H.; Singh, V.P. One dimensional velocity distribution in open channels using Tsallis entropy. J. Hydrol. Eng. 2014, 19, 290–298. [Google Scholar] [CrossRef]
  40. Luo, H.; Singh, V.P. Entropy theory for two-dimensional velocity distribution. J. Hydrol. Eng. 2011, 16, 303–315. [Google Scholar] [CrossRef]
  41. Kumbhakar, M.; Ghoshal, K. One-dimensional velocity distribution in open channels using Renyi entropy. Stochastic. Environ. Res. Risk Assess. 2017, 31, 949–959. [Google Scholar] [CrossRef]
  42. Chiu, C.L.; Jin, W.; Chen, Y.C. Mathematical models of distribution of sediment concentration. J. Hydraul. Eng. 2000, 1, 16–23. [Google Scholar] [CrossRef]
  43. Cui, H.; Singh, V.P. Suspended sediment concentration in open channels using Tsallis entropy. J. Hydrol. Eng. 2013, 19, 966–977. [Google Scholar] [CrossRef]
  44. Kumbhakar, M.; Ghoshal, K.; Singh, V.P. Derivation of Rouse equation for sediment concentration using Shannon entropy. Physics A 2017, 465, 494–499. [Google Scholar] [CrossRef]
  45. Sterling, M.; Knight, D. An attempt at using the entropy approach to predict the transverse distribution of boundary shear stress in open channel flow. Stochastic. Environ. Res. Risk Assess. 2002, 16, 127–142. [Google Scholar] [CrossRef]
  46. Bonakdari, H.; Sheikh, Z.; Tooshmalani, M. Comparison between Shannon and Tsallis entropies for prediction of shear stress distribution in open channels. Stochastic. Environ. Res. Risk Assess. 2015, 29, 1–11. [Google Scholar] [CrossRef]
  47. Shannon, C.E. A mathematical theory of communications, I and II. Bell Syst. Tech. J. 1948, 27, 379–423. [Google Scholar] [CrossRef]
  48. Maggi, F. Stochastic flocculation of cohesive sediment: analysis of floc mobility within the floc size spectrum. Water Resour. Res. 2008, 440, 168–182. [Google Scholar] [CrossRef]
  49. Shin, H.J.; Son, M.W.; Lee, G.H. Stochastic flocculation model for cohesive sediment suspended in water. Water 2015, 47, 2527–2541. [Google Scholar] [CrossRef]
  50. Jaynes, E.T. Information theory and statistical mechanics I. Phys. Rev. 1957, 106, 620–630. [Google Scholar] [CrossRef]
  51. Jaynes, E.T. Information theory and statistical mechanics II. Phys. Rev. 1957, 108, 171–190. [Google Scholar] [CrossRef]
  52. Jaynes, E.T. On the rationale of maximum entropy methods. Proc. IEEE 1982, 70, 939–952. [Google Scholar] [CrossRef]
  53. Chiu, C.L. Entropy and probability concepts in hydraulics. J. Hydraul. Eng. 1987, 113, 583–600. [Google Scholar] [CrossRef]
  54. Khozani, Z.S.; Bonakdari, H. Formulating the shear stress distribution in circular open channels based on the Renyi entropy. Physics A 2018, 490, 114–126. [Google Scholar] [CrossRef]
  55. Singh, V.P. A Shannon entropy-based general derivation of infiltration equations. Trans. ASABE 2011, 54, 123–129. [Google Scholar] [CrossRef]
  56. Tsallis, C. Possible generalization of Boltzmann-Gibbs statistics. J. Stat. Phys. 1988, 52, 479–487. [Google Scholar] [CrossRef]
  57. Singh, V.P. Tsallis entropy theory for derivation of infiltration equations. Trans. ASABE 2010, 53, 447–463. [Google Scholar] [CrossRef]
  58. Burban, P.Y.; Lick, W.; Lick, J. The flocculation of fine—Grained sediments in estuarine waters. J. Geophys. Res. Ocean 1989, 94, 8323–8330. [Google Scholar] [CrossRef]
  59. Spicer, P.T.; Pratsinis, S.E. Shear induced flocculation: the evolution of floc structure and the shape of the size distribution at steady state. Water Res. 1996, 30, 1048–1056. [Google Scholar] [CrossRef]
  60. Selomulya, C.; Bushell, G.; Amal, R.; Waite, T.D. Aggregation mechanisms of rates of different particle sizes in a controlled shear environment. Langmuir 2002, 18, 1974–1984. [Google Scholar] [CrossRef]
  61. Colomer, J.; Peters, F.; Marrasé, C. Experimental analysis of coagulation of particles under low-shear flow. Water Res. 2005, 39, 2994–3000. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Figure 1. The flocculation element with floc size growth D ( t ) = floc size entering the flocculation element, D ( t ) = floc size exiting the element, and S = the capacity of floc size growth.
Figure 1. The flocculation element with floc size growth D ( t ) = floc size entering the flocculation element, D ( t ) = floc size exiting the element, and S = the capacity of floc size growth.
Entropy 20 00845 g001
Figure 2. The comparison of the proposed entropy-based expression (Equation (12)) with thirty-three experimental data sets from the literature. In each figure, the magenta circles denote the measured data and the black line represents the proposed expression. For the cases of T23, T24, and T25 from Selomulya et al. [60], the horizontal axis is not the flocculation time t but a quantity t * T 0.45 , where T is the absolute temperature in the flocculation time (the unit is Kelvin), and we cannot obtain the value of the flocculation time from their paper.
Figure 2. The comparison of the proposed entropy-based expression (Equation (12)) with thirty-three experimental data sets from the literature. In each figure, the magenta circles denote the measured data and the black line represents the proposed expression. For the cases of T23, T24, and T25 from Selomulya et al. [60], the horizontal axis is not the flocculation time t but a quantity t * T 0.45 , where T is the absolute temperature in the flocculation time (the unit is Kelvin), and we cannot obtain the value of the flocculation time from their paper.
Entropy 20 00845 g002aEntropy 20 00845 g002bEntropy 20 00845 g002c
Figure 3. The comparison of the entropy-based expression with the model of Keyvani and Strom for measured data in each of the cycles: (a) ps1, (b) ps2, (c) ps3, (d) ps4, (e) ps5, (f) ps6, and (g) ps7 in the work of Keyvani and Strom [33] (ps was referred to as a “prior shear” case, corresponding to the cycle order of the high and low turbulent shear, in their paper). The blue circle denotes the measured data, the black line represents the entropy-based expression, and the red line shows the model of Keyvani and Strom.
Figure 3. The comparison of the entropy-based expression with the model of Keyvani and Strom for measured data in each of the cycles: (a) ps1, (b) ps2, (c) ps3, (d) ps4, (e) ps5, (f) ps6, and (g) ps7 in the work of Keyvani and Strom [33] (ps was referred to as a “prior shear” case, corresponding to the cycle order of the high and low turbulent shear, in their paper). The blue circle denotes the measured data, the black line represents the entropy-based expression, and the red line shows the model of Keyvani and Strom.
Entropy 20 00845 g003aEntropy 20 00845 g003b
Figure 4. The fitted parameter value S in the proposed entropy-based model with respect to different flow shear rate G values for the collected experimental data in the log-log space (a) and the normal space (b).
Figure 4. The fitted parameter value S in the proposed entropy-based model with respect to different flow shear rate G values for the collected experimental data in the log-log space (a) and the normal space (b).
Entropy 20 00845 g004
Figure 5. The fitted parameter value S in the proposed entropy-based model with respect to the different bed shear stress τ b values for the experimental data from Stone and Krishnappan [30].
Figure 5. The fitted parameter value S in the proposed entropy-based model with respect to the different bed shear stress τ b values for the experimental data from Stone and Krishnappan [30].
Entropy 20 00845 g005
Figure 6. The comparison between the observed floc size D and the estimated floc size D using Equation (22) for the experimental data from (a) Oles [24], (b) Serra et al. [12], (c) Serra and Casamitjana [31], (d) Biggs and Lant [14], (e) Colomer et al. [61], and (f) Stone and Krishnappan [30].
Figure 6. The comparison between the observed floc size D and the estimated floc size D using Equation (22) for the experimental data from (a) Oles [24], (b) Serra et al. [12], (c) Serra and Casamitjana [31], (d) Biggs and Lant [14], (e) Colomer et al. [61], and (f) Stone and Krishnappan [30].
Entropy 20 00845 g006aEntropy 20 00845 g006b
Table 1. The information on the collected experimental data in the literature.
Table 1. The information on the collected experimental data in the literature.
Experimental Data NumberExperimental MaterialTurbulence-Generating EnvironmentFlow Shear Condition D 0   ( μ m ) D   ( μ m ) Data Source
T1Detroit river sedimentCouette-flow chamber ϕ = 1.04 × 10−4; G = 200 s 1 487Burban et al. [58]
T2 ϕ = 1.66 × 10−3; G = 200 s 1 425.21
T3Polystyrene latexCouette-flow system formed by two cylinders ϕ = 5 × 10−5; G = 75 s 1 2.1739.54Oles [24]
T4 ϕ = 5 × 10−5; G = 100 s 1 2.1736.65
T5 ϕ = 5 × 10−5; G = 125 s 1 2.1726.52
T6 ϕ = 5 × 10−5; G = 150 s 1 2.1714.47
T7Polystyrene particleBaffled stirred tank ϕ = 2.10 × 10−5; G = 63 s−1
Alum concentration: 4.3 mg/L
0.8713.54Spicer and Pratsinis [59]
T8 ϕ = 2.10 × 10−5; G = 63 s−1
Alum concentration: 10.7 mg/L
0.8741.90
T9 ϕ = 2.10 × 10−5; G = 63 s−1
Alum concentration: 32 mg/L
0.8784.20
T10 ϕ = 2.10 × 10−5; G = 95 s−1
Alum concentration: 32 mg/L
0.8767.01
T11Latex particleCouette-flow system ϕ = 2.5 × 10−5; G = 25 s 1 246.06Serra et al. [12]
T12 ϕ = 2.5 × 10−5; G = 50 s 1 238.84
T13 ϕ = 2.5 × 10−5; G = 90 s 1 230
T14 ϕ = 2.5 × 10−5; G = 135 s 1 219.87
T15 ϕ = 2.5 × 10−5; G = 195 s 1 211.74
T16Latex particleCouette-flow system ϕ = 5 × 10−5; G = 25 s 1 241.36Serra and Casamitjana [31]
T17 ϕ = 5 × 10−5; G = 32 s 1 237.73
T18 ϕ = 5 × 10−5; G = 50 s 1 235.23
T19Activated sludgeBaffled batch vessel ϕ = 5 × 10−2; G = 19.4 s 1 15 ***121.27Biggs and Lant [20]
T20 ϕ = 5 × 10−2; G = 37 s 1 15 ***100.56
T21 ϕ = 5 × 10−2; G = 113 s 1 15 ***58.66
T22 ϕ = 5 × 10−2; G = 346 s 1 15 ***24.14
T23Polystyrene latex particleCouette-flow system ϕ = 3.76 × 10−5; G = 64 s−10.8170.94Selomulya et al. [60]
T24 ϕ = 3.76 × 10−5; G = 100 s−10.8167.76
T25 ϕ = 3.76 × 10−5; G = 246 s−10.8138.07
T26Hay river sediment, CanadaAnnular flumeBed shear stress = 0.123 Pa19.1128.97Stone and Krishnappan [30]
T27Bed shear stress = 0.212 Pa19.1178.1
T28Bed shear stress = 0.323 Pa19.1161.84
T29Polystyrene latex particleFlask shaking table ϕ = 2 × 10−5; G = 0.45 s 1 2.17.88Colomer et al. [61]
T30 ϕ = 2 × 10−5; G = 0.75 s 1 2.19.34
T31 ϕ = 2 × 10−5; G = 0.96 s 1 2.19.05
T32 ϕ = 2 × 10−5; G = 1.41 s 1 2.19.68
T33 ϕ = 2 × 10−5; G = 2.4 s 1 2.110.42
The “***” symbol indicated that the measured size by Biggs and Lant [14] at the beginning of the flocculation experiment is actually the floc size of 15 microns rather than the size of the primary particle (the primary particle size is actually 4 microns).
Table 2. The comparison results of the proposed entropy-based expression with the collected experimental data in the literature.
Table 2. The comparison results of the proposed entropy-based expression with the collected experimental data in the literature.
Experimental Data NumberData SourceFitting ResultEntropy Function
R 2   RBIASRMSE H S ( D )   H T ( D )
Assume m = 2
T1Burban et al. [58]0.9750.0544.1704.419 82.988
T20.9950.0230.6403.054 21.163
T3Oles [24]0.9440.2133.1343.621 37.343
T40.9480.1602.3413.540 34.451
T50.9890.0800.9603.193 24.309
T60.9820.0440.5122.510 12.219
T7Spicer and Pratsinis [59]0.9620.0761.0532.539 12.591
T80.9640.0693.2803.714 41.006
T90.9990.0141.5114.423 83.318
T100.9780.0394.0384.192 66.125
T11Serra et al. [12]0.9810.1182.4453.786 44.037
T120.9620.1213.0283.607 36.813
T130.9760.0451.2613.332 27.964
T140.9580.0440.9542.883 17.814
T150.8500.0761.0782.276 9.637
T16Serra and Casamitjana [31]0.8990.1213.6063.673 39.335
T170.9520.1042.6003.576 35.702
T180.9560.0722.0193.503 33.200
T19Biggs and Lant [14]0.9800.0273.4034.666 106.261
T200.9670.0374.1264.449 85.548
T210.9600.0362.0873.776 43.637
T220.9720.0170.5212.213 9.031
T23Selomulya et al. [60]0.8450.1247.6074.250 70.116
T240.8990.0413.6234.204 66.935
T250.9790.0191.1063.618 37.233
T26Stone and Krishnappan [30]0.8870.08513.3044.699 109.861
T270.9740.0358.3515.069 158.994
T280.9840.0235.9884.961 142.733
T29Colomer et al. [61]0.9930.0210.1891.754 5.607
T300.9920.0350.2861.980 7.102
T310.9930.0380.3041.939 6.806
T320.9940.0210.2202.026 7.448
T330.9880.0320.3502.119
Table 3. The three simplified Lagrangian flocculation models.
Table 3. The three simplified Lagrangian flocculation models.
Model NameFormulation
Winterwerp model d D d t = c ρ s k A d f G D 0 F 3 D 4 F k B D 0 d f μ F G 1.5 D 2 ( D D 0 )
Son and Hsu (2008) model d D d t = G D 0 β β ln D D 0 + 1 [ c 3 ρ s k A D 0 d f 3 D 4 d f β k B 3 D 0 μ G F D 2 β ( D D 0 ) ]
Son and Hsu (2009) model d D d t = G D 0 β β ln D D 0 + 1 [ c 3 ρ s k A D 0 d f 3 D 4 d f β k B 3 μ G B D 0 d f 3 1 D 1 β + 3 d f 3 ( D D 0 ) ]
Table 4. The comparison of the present entropy-based model with the deterministic models for the experimental data.
Table 4. The comparison of the present entropy-based model with the deterministic models for the experimental data.
ReferencesExperimental ConditionsFitting Effect
The Present Model Winterwerp Model Son and Hsu (2008) Model Son and Hsu (2009) Model
R 2   RBIASNRMSE R 2   RBIASNRMSE R 2   RBIASNRMSE R 2   RBIASNRMSE
Burban et al. [58] ϕ = 1.04 × 10−4; G = 200 s 1 0.980.0544.1700.830.28219.8600.860.25517.5870.900.19012.942
ϕ = 1.66 × 10−3; G = 200 s 1 0.990.0230.6400.970.0371.0830.970.0260.4240.980.0361.053
Biggs and Lant [14] ϕ = 5 × 10−2; G = 19.4 s 1 0.980.0273.4030.890.0537.2180.900.0597.9170.900.0678.889

Share and Cite

MDPI and ACS Style

Zhu, Z. A Simple Explicit Expression for the Flocculation Dynamics Modeling of Cohesive Sediment Based on Entropy Considerations. Entropy 2018, 20, 845. https://doi.org/10.3390/e20110845

AMA Style

Zhu Z. A Simple Explicit Expression for the Flocculation Dynamics Modeling of Cohesive Sediment Based on Entropy Considerations. Entropy. 2018; 20(11):845. https://doi.org/10.3390/e20110845

Chicago/Turabian Style

Zhu, Zhongfan. 2018. "A Simple Explicit Expression for the Flocculation Dynamics Modeling of Cohesive Sediment Based on Entropy Considerations" Entropy 20, no. 11: 845. https://doi.org/10.3390/e20110845

APA Style

Zhu, Z. (2018). A Simple Explicit Expression for the Flocculation Dynamics Modeling of Cohesive Sediment Based on Entropy Considerations. Entropy, 20(11), 845. https://doi.org/10.3390/e20110845

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop