1. Introduction
The rheological behavior of a cement paste depends principally on the contents of solid particles and their capacity to form flocs as a consequence of the particle interactions, which can be modulated by adding adequate dispersing agents to improve the dispersion state. As the particle contents increases, the value of dynamic viscosity,
, also increases [
1]. If these systems are considered homogenous suspensions, the Krieger–Dougherty equation [
1,
2] gives the dynamic viscosity from the volume fraction of cement particles, which may be obtained from the water/cementitious materials relationship,
[
3]. The interest and application of the Krieger–Dougherty equation in these suspensions are due, from a theoretical point of view, to its robustness and the fact that its parameters have a physical significance [
1,
4,
5,
6]. The Krieger–Dougherty equation [
1,
2,
6] is defined as follows:
where,
: dynamic viscosity of the suspension.
: dynamic viscosity of the fluid phase of the suspension.
: fraction in volume of the disperse solid phase of the suspension.
: maximum packing fraction of particles in the disperse solid phase.
: intrinsic viscosity, which depends on the specific viscosity (ratio between the viscosity of the suspension and the dispersion liquid) and the volume fraction of solids.
Much research has been performed on the rheology of the various types of cement suspensions in which the suitability of this equation for the calculation of dynamic viscosity has been verified [
2,
4,
5,
7,
8]. In addition, a comparison has been made of the values calculated with the experimental measures performed with rheometers, thereby obtaining good results with cement pastes with and without superplasticizer admixtures [
1,
3]. The Krieger–Dougherty equation [
1,
2,
6] depends on two parameters: the maximum packing fraction,
, and the intrinsic viscosity,
. The first is a measurement of the maximum packing that may be reached in a particle system, i.e., the maximum concentration of particles that can be added while maintaining some flowability. It depends on the dispersion conditions, on the shape of the particles and, especially, on their distribution and size [
3,
6]. The second measures the individual effect of the particles on the viscosity and is secondary to their form [
1]. Moreover, both parameters are affected by the shear rate,
, applied in the system: as
increases,
tends to increase while
shows the opposite effect [
6]. Nevertheless, it is certain that the behavior of
is very dependent on the slenderness of the particles [
1,
9].
Generally, in cement suspensions both
as well as
may be determined from the adjustment of a collection of experimental measurements, thereby assuming the hypothesis of sphericity and rigidity of the particles [
10,
11,
12] and the formation of monodisperse or polydisperse systems that adapt geometric packing of known values (cubic, random, or hexagonal compact packing) [
7]. In such cases, it is common practice to consider values ranging between 0.6–0.7 for
, and equal to 2.5 for
(in monodisperse systems), and greater than 0.7 for
, and equal to 2.5 for
when polydispersion increases [
1,
13]. However, as the asymmetry of the particles increases,
may reach values below 0.6,
between 3 and 5 when the particles are sharp and approximately equal, and between 4 and 10 in particles with more acicular shapes [
1,
13]. Also, the trend in the behavior of both parameters when shearing the system suggests that the product of the two, which appears in the exponent of the Krieger–Dougherty equation, remains practically constant in any case [
7,
8]. All the mentioned values for the parameters follow from the hypothesis on the shape and type of arrangement of the particles, rather than from experiments. Besides, the variations of these parameters due to the random shape and arrangement of the solid phases are not known, which makes the application of the Bayesian analysis to know their probability distributions desirable.
Since a physical system is not deterministic, but random, the models that simulate its behavior should be probabilistic, if a good approximation to reality is sought. Therefore, the application of parametric-deterministic models should be supplemented with parametric-probabilistic models, which signifies added value within the field of modeling in engineering practice. For this reason, it is adequate to identify the sources of randomness that are associated with these systems, observe them, and perform trials and collect data to establish processes for the construction of probabilistic models [
14,
15]. The transformation of parametric-deterministic models into parametric-probabilistic models is usually performed using Bayesian analysis methods. Due to their complexity, these methods were limited a short time ago to the conjugate family of functions, which were the only ones for which expressions of posterior distributions functions could be easily determined. Nevertheless, the appearance of Markov Chain Monte Carlo methods and the Gibbs Sampling has enabled the simulation of the posterior distribution directly, thereby generalizing its application [
16,
17,
18,
19].
Within the scope of the study of the rheology of systems that are based on cement materials, and more specifically, on suspensions of cement paste, the colloidal nature of the particles that intervene, and the physical–chemical interactions that are adjusted as a consequence of the inclusion of admixtures, make these suspensions especially attractive for the rheology study from a probabilistic point of view, and they may be analyzed through Bayesian methods. In the same manner, more complex cement systems, like self-compacting mortar and self-compacting concrete, show the same tendency insofar as the relationship between viscosity and particle content, even with some differences associated with the presence of larger size solids, such as aggregates.
These analyses are of great engineering relevance [
20], since self-compactibility, passing ability, pumpability, etc. are properties that depend on the viscosity and, in general, on the rheological properties of the mixture. In turn, these properties depend on the composition of the suspension, particularly water (its hardness and contents, especially metallic ions), and the geometry and compactness of cementitious and aggregate materials in suspension. For these reasons, the principal objective of this research consists of the transformation of the Krieger–Dougherty equation into a parametric-probabilistic model using Bayesian analysis, and to apply it to cement pastes, self-compacting mortars, and self-compacting concretes. We also want to use the new model to compute the probability density functions that intervene in the Krieger–Dougherty equation applied to these cementitious suspensions.
The article is structured as follows: Firstly, we define the characteristics of the Bayesian analysis and its suitability for application to the Krieger–Dougherty equation. Next, the paper describes the procedure performed with the methodology. Subsequently, there is a description of the experimental campaign undertaken and the scientific literature data used, which is followed by the results and its discussion. Lastly, the conclusions drawn from the research are explained.
2. Probabilistic Model and Bayesian Analysis of the Krieger–Dougherty Equation
The models of probabilistic networks are an appropriate methodology for dealing with problems in the engineering practice [
15] since the reality is random, not deterministic. When we deal with multivariate random variables, our aim is to obtain the multivariate density or probability function, because if these functions are known we can answer any probability question about this variable, for example, the marginal densities of any subset of variables, the regression expressions, any conditional distributions, etc. However, the definition of a multivariate random variable is neither trivial nor easy. Some definitions do not guarantee the existence of multivariate distributions. For example, the conditionals of
x given
y, and
y given
x can be incompatible. The best way of defining multivariate distributions is by means of a Bayesian network, because they always guarantee the existence of the joint multivariate density and, in addition, the definition of this density is made locally, that is, in small pieces that always guarantee this existence.
The objective of the equation of Krieger–Dougherty [
2] (Equation (
1)) is the analytical calculation and prediction of the viscosity as a function of the volume fraction of solids in suspensions of different types. Within systems based on cement materials, this equation has been used to calculate the viscosity of cement pastes [
1] and self-compacting concrete [
7]. In dealing with an equation in which two parameters (
and [
]) intervene, they may be adjusted to the experimental data [
10], which exhibits variety in the values that they may acquire and that depend on the characteristics of the suspensions. The physical-chemical nature of the materials and the interactions that are produced among them, inherently exhibit a collection of responses of a random nature. Thus, the parameters controlling such responses may be deemed random variables that follow density functions of probability, in lieu of being defined with a single value.
Therefore, the equation of Krieger–Dougherty may be expressed in probabilistic terms, thereby obtaining information on the variability of the estimation of
. If the matter is dealt with by using a classic methodology,
may be treated as a random variable of a parametric family, thereby reducing the problem to an estimation of equation parameters. Nevertheless, if a Bayesian methodology is applied, a series of parametric distribution families are considered, and their parameters are treated as random variables [
21].
2.1. Sources of Randomness in Cementitious Suspensions
A colloidal suspension is a system of two or more phases formed by a fluid dispersing medium and dispersed particles with diameters between
and
m [
22,
23]. They vary from large molecules, such as polymer chains of superplasticizer admixture, to small particles, such as cement materials and mineral fillers [
22,
24]. The shape and the size distribution of the particles, together with the surface chemistry and the interaction forces among them and with the dispersing medium, determine the properties of the suspension [
22]. The forces of interaction (van der Waals forces, electrostatic repulsion forces, steric hindrance, and Brownian forces) dominate on the inertials and the gravitationals at this scale [
22]. This constitutes one of the sources of randomness in suspensions of cement paste together with the random nature itself of the Brownian movement [
1,
9]. Another of the sources arises from the variety of shapes (pseudo-spherical, angular, elongated etc.) and sizes of the particles of the cement materials [
1,
25]. The shear rate,
, to which the system is subjected must also be taken into account. It has a direct influence on the packing and interaction among the particles, which is related to the values that the parameters
and
adopt [
1,
6,
7,
8].
Regarding self-compacting mortar (SCM) and self-compacting concrete (SCC), they may be deemed as systems composed of a solid granular phase of one or various sizes, respectively, that is found in suspension in a continuous fluid viscous phase, such as cement paste [
6,
7,
26]. The cement paste has a pronounced random nature as has been explained before. The rheological behavior of SCM and SCC is conditioned by the interactions between the aggregates and by the physical–chemical properties of the cement paste [
26]. Thus, performing a multi-scale approach seems appropriate since the different phases are defined by the maximum size of their particles [
26], and each exhibits a characteristic rheology.
If the hypothesis of considering SCM and SCC granular skeleton as a collection of rigid non-colloidal polydisperse spheres is proposed with respect to the distribution of the size of the particles, the viscosity of this system may be determined using the Krieger–Dougherty equation, Equation (
1). As the fluid phase is always water (and thus
is not supposed to vary much), Equation (
1) depends mainly on two parameters with physical significance: the maximum packing fraction of the particles,
, and the intrinsic viscosity of the system, [
] [
6,
10].
The maximum packing fraction of the disperse solid phase,
, is related to the particle size distribution and their shape [
1,
6]. Thus, in a system of monodisperse rigid spherical particles, the value of
is approximately equal to 0.648, regardless of the size of the sphere. Under this hypothesis, Hu et al. [
27] proposed an equation that can be used as a first approximation to the value of
. In polydisperse systems, the value of
increases with particle size variability as the space between them is filled more efficiently (
≈ 0.744) [
6,
7].
The intrinsic viscosity,
, depends on the individual effect of particles and their shape on viscosity [
1,
6]. It takes a value equal to 2.5 for rigid spherical geometries [
6]; when the particles deviate of this shape,
must be modified [
1,
13,
28,
29,
30]. The expression suggested by Pabst et al. [
28] can be useful to estimate
. Anyway, the correlation between particle shape and
is fundamentally complicated because, on the one hand, it is normally assumed that all particles have a similar shape and, on the other hand, the fit of
in the Krieger–Dougherty equation to the experimental measurements is subject to statistical and systematic errors [
28].
As the aggregates move away from the spherical shape, other values of
must be used to simulate the actual shape of the particles. Szecsy [
30,
31] established a relationship between
and the circularity of the particle, defined as the relationship between the area and the perimeter of the particle using digital image processing techniques. It is of utmost interest to analyze what occurs with the various concentrations of solids, especially when their values are near the maximum packing fraction, a point that corresponds to the volumetric fraction in which a strong increase of the dynamic viscosity occurs [
10]. In the same manner that occurs with cement paste, the variety of shapes and sizes of granular particles is a real fact that introduces a source of randomness related to aggregates to the system.
When the system is diluted the viscosity depends principally on the concentration of particles. Its value does not vary too much when
increases since the particles are not close to each other and the hydrodynamic interactions may be disregarded [
6,
32]. With high concentration of aggregates, which is what happens in SCC, the hydrodynamic interactions that occur between them as a consequence of the shear applied to the suspension [
6] must be taken into account in addition to the distribution of shapes and sizes. The parameters
and
of the Krieger–Dougherty equation depend on
, and its product remains practically constant when the aggregates are deemed rigid spheres (it takes an average value equal to 1.9 [
6] or 2 [
10,
33]). Thus, this energy introduces a source of randomness to the system as occurs with the systems comprised of only cement paste.
2.2. Bayesian Approach
The Bayesian networks are used in a multitude of disciplines and practical problems [
15,
21], in which the analysis and interpretation of data are important in taking decisions [
16,
17,
18,
19]. A Bayesian network consists of two elements, (
); the first,
, is a directed acyclic graph, that has the variables as nodes, and the links, the direct dependencies among the variables. Due to the directed acyclic graph, the variables are ordered, and each variable has no parents or a set of parents, normally a few, that are the variables on which there is a direct dependency. The acyclic graph permits answering questions such as: has a subset of variables
A new information on another subset
B, when a third subset
C is already known? This is a very interesting and useful property. Once we have the graph, which defines the qualitative information on the network, we need to quantify probabilities, and this is done by the second element,
, which is a set of conditional probabilities, one per variable, that gives the conditional probabilities of the variables conditioned on their parents. From this set we can obtain the joint density by multiplying all of the conditional probabilities, that is, we have:
where
are the subsets of parents of the variables
. Thus, defining a multivariate density means defining the directed acyclic graph, and defining each of the conditional probabilities, one per variable. However, these conditional probabilities are local, because normally involve a small number of variables, that is, easy to be defined. In Bayesian methods, the parameters of these conditional probabilities are assumed random, and then, they are converted into random variables, and incorporated into the Bayesian network. It is clear that, apart from some particular cases, the calculus of probabilities is complicated and normally involves integrals, that lead to very complex problems that have no analytical solution. One way to avoid these calculations is by using simulations. We simulate a very high number of samples and use the sample of the variable, marginal, or conditional, we are interested in, and we inform the user by means of a very large sample, instead of an analytical expression, that in practice is equivalent. The Markov Chain Monte Carlo methods optimize the sampling procedure using some asymptotic properties of the simulation procedures that reduce the sample size drastically. This is the method that is used in this paper.
The Bayesian methods deal with parametric families of distributions, the parameters of which are considered random variables [
19,
21]. These models do not choose a model of the initial family of distributions, but rather a linear convex combination of the different models of the selected family. This aspect is very important since it permits increasing the collection of models and enables the sample to indicate which is the most appropriate [
18,
19].
The Bayesian approach of a probabilistic model consists of the following steps [
34]:
Selection of the likelihood family.
Selection of the prior distribution of the parameters. It deals with a very important step in the methodology since the results for small samples strongly depend on it. The selection of this prior information may be done in the following manner:
- −
Using an imaginary sample. For that, an expert is consulted for the purpose of providing a virtual sample of a certain size as the most representative to their prior knowledge.
- −
Using previous non-updated information available in the specialized scientific literature.
- −
Based on data obtained in our own experiments.
Obtaining data from the sample.
Calculation of the posterior distribution.
Combining the posterior with the likelihood to obtain the predictive distribution, which is the one we used.
2.3. Formulation of the Probabilistic Model and Bayesian Analysis of the Krieger–Dougherty Equation
The proposed objective consists of converting the Krieger–Dougherty equation into a parametric model and performing a Bayesian analysis of it. In the first place, it will deal with the case of homogenous suspensions of cement paste. Below, the same procedure is performed for the case of self-compacting mortar and self-compacting concrete, thereby considering them as heterogeneous suspensions of granular particles within a homogeneous fluid, which is the cement paste.
Using the Krieger–Dougherty Equation (
1), dynamic viscosity of cementitious systems can be calculated; two main parameters intervene,
and
, which are going to be dealt with as random variables in the Bayesian analysis. In order for the proposed model to provide reasonable results, it is necessary to have prior adequate information, which may be obtained from the experimental data or consultation with the scientific literature. The quality of the information is fundamental, especially when sufficient data are not available [
34]. The Bayesian network of the model, which is to be created, will take into account the random nature of the average value of the dynamic viscosity as well as the variability of the parameters that intervene in the model.
The transformation of the parametric-deterministic model into a parametric-probabilistic model and the Bayesian analysis has been undertaken in this work using the open-source code, OpenBUGS [
35]. It involves a Bayesian inference program that uses Markov Chain Monte Carlo methods and the Gibbs Sampling as a basis (the Gibbs Sampling is a particular case of a simulation algorithm of a Markov Chain). These methods successively simulate the density function that has been proposed, which does not necessarily have to be similar to the posterior density function. Each value generated only depends on the value simulated previously (thus the denomination of Markov Chain). Besides, OpenBUGS permits the production of the graph or diagram of the Bayesian network of the model in question. The program simulates the posterior distribution of the parameters of a model, which requires the definition thereof, the incorporation of a collection of data and beginning values, the latter of which are an important aspect in the analysis of the quality of the simulations performed [
18]. Insofar as results, the program provides the graph of the hierarchical dependence structure between variables and parameters, the functions of the posterior density of the parameters, and a collection of statistics of the probabilistic model.
2.3.1. Cement Paste Suspensions
The first step is the definition of the Krieger–Dougherty model [
2] in a dimensionless format:
where,
: non-dimensional dynamic viscosity of the cement paste.
: dynamic viscosity of the cement paste.
: dynamic viscosity of the continuous fluid phase of the suspension (water).
: fraction in volume of the disperse solid phase of the suspension (cementitious materials).
: maximum packing fraction of the particles of the disperse solid phase.
: intrinsic viscosity of the cement paste.
To perform the conversion of the initial deterministic model into a probabilistic model and perform the Bayesian analysis, the variables
and
are considered independent random variables, which belong to a family of uniform probability density functions defined within a domain (minimum and maximum values of the validity interval). In addition, it is assumed that the Krieger–Dougherty equation provides the mean of the values of the random variable,
, which follows a normal density function (with mean
and standard deviation
). The observed residue,
, also follows a normal density function and
is to be defined by means of a non-informative uniform density function. Therefore, the syntax with which the model is to be defined in statistical notation is the following:
where
is an auxiliary variable depending on the standard deviation (
).
2.3.2. Self-Compacting Mortar Suspensions
Again, the first step is the definition of the Krieger–Dougherty model [
2] for self-compacting mortar suspensions in a dimensionless format. From
, and through the application of a micromechanical model, the viscosity of any cementitious paste with granular phases in suspension (self-compacting mortar and self-compacting concrete) can be calculated considering it as a two-phase suspension of particles within a viscous fluid [
6,
7]. The addition of each solid phase produces an increase of the viscosity of the fluid phase, which can be calculated by successively applying the Krieger–Dougherty equation when including each one of the granular phases. For the case of self-compacting mortar suspensions, the viscous fluid is the cement paste and the fine aggregate constitutes the only solid phase, thus the model can be expressed through [
7]:
where,
: non-dimensional dynamic viscosity of self-compacting mortar.
: dynamic viscosity of self-compacting mortar.
: dynamic viscosity of cement paste.
: fraction in volume of the granular phase of the suspension (fine aggregate).
: maximum packing fraction of particles of the granular phase.
: intrinsic viscosity of the system when adding the granular phase.
As occurred with the cement pastes, each one of the parameters of Equation (
9) are considered random variables that follow a uniform density function of probability within a domain of validity (with minimum and maximum values defined for each parameter). Likewise, Equation (
9) determines the mean of the dynamic viscosity of self-compacting mortar, which follows a normal density function of probability (average,
, standard deviation,
). The residual,
, also belongs to a normal family and comprises a density function of the uniform type. Thereby the syntax of the model in statistical notation is defined in this way:
2.3.3. Self-Compacting Concrete Suspensions
Similarly to self-compacting mortar, self-compacting concrete can be considered as a two-phase suspension of granular particles, in which the addition of each solid phase produces an increase of the dynamic viscosity. Again, applying successively the Krieger–Dougherty equation when including each one of the solid phases (fine aggregate and coarse aggregate) the model is defined as follows [
7]:
where,
: non-dimensional dynamic viscosity of self-compacting concrete.
: dynamic viscosity of self-compacting concrete.
: dynamic viscosity of cement paste.
: fraction in volume of the finer granular phase of the suspension (fine aggregate).
: maximum packing fraction of particles of the finer granular phase.
: intrinsic viscosity of the system when adding the finer granular phase.
: fraction in volume of the thicker granular phase of the suspension (coarse aggregate).
: maximum packing fraction of the thicker granular phase of the suspension (coarse aggregate).
: intrinsic viscosity of the system when adding the coarser granular phase.
Similarly to cement pastes and self-compacting mortars, the parameters of the Krieger–Dougherty equation in self-compacting concretes are treated as random variables with a uniform density function of probability within a domain of validity (with minimum and maximum values defined for each parameter). Thus, the Krieger–Dougherty equation gives the mean of the dynamic viscosity of self-compacting concrete, which follows normal density function of probability (average,
, standard deviation,
). Again, the residual,
, belongs to a normal family and comprises a density function of the uniform type. The syntax of the model in a statistical format is shown below:
5. Conclusions
We carried out the transformation of the Krieger–Dougherty equation into a probabilistic model using a Bayesian analysis methodology. The open-source software OpenBUGS was used, which employs Markov Chain Monte Carlo and Gibbs Sampling to perform the simulations to obtain the samples of the variables that can be interpreted as their density functions, which are the parameters of the deterministic models. The methodology has been applied to cement pastes, self-compacting mortars, and self-compacting concretes. The density functions of the parameters (maximum packing fractions of the phases involved, and their corresponding intrinsic viscosities ) are conditioned by the Bayesian graph, which describes the hierarchy and independence of variables, by the limits of the initial uniform distributions and by the limits of the final distribution. The analysis does not attribute the resulting distribution to a single cause (for example, the variations in the shape of the aggregate), but to all those that can physically condition the values of the parameters.
In particular, the Bayesian method applied to the cement pastes in De La Rosa et al. [
5] confirms that the theoretical values give reasonable results and, for the first time, calculates the distribution function of the parameters of the Krieger–Dougherty equation. The results obtained with the self-compacting mortars of Ouro et al. [
38] also confirm that theoretical values give a good approximation, and the abrupt shape of the distribution function for
indicates that the used sand was round and spherical. Additionally, the error in the dynamic viscosity predictions using the mean values of the distribution curves (33%) is less than the error with the theoretical values (43%) and the standard deviation is also reduced.
We also applied the Bayesian methodology to three series of self-compacting concretes reported by Feys et al. [
39], Esmaeilkhanian et al. [
40], and Grünewald [
41]. The Bayesian results detect that the sand used by Feys et al. [
39] was not round but crushed, since the maximum packing fraction
is smaller than in the other cases (0.64 versus 0.68) and the distribution is smoother. Similarly, the methodology reveals that the sand used by Grünewald [
41] was less spherical than the others due to the high value obtained for
. The study also reveals that the coarse aggregates used in the three SCC series must have low sphericity due to the high values obtained for
in all of them. Regarding the exponents of the Krieger–Dougherty equation, the one for fine aggregates is closer to the theoretical value (1.9) for the fine aggregate in the three series, whereas the uniform distribution for the exponents corresponding to the coarse aggregates reveals, again, the disparity in the shape of the crushed particles. Finally, the dynamic viscosity predictions made with the mean values of the Bayesian distributions were significantly better than those made with the theoretical values. The error diminished from 77% to 25% in the Feys et al. [
39] series, from 42% to 36% in the Esmaeilkhanian et al. [
40] SCCs and from 71% to 17% in the Grünewald [
41] concretes.