Next Article in Journal
Challenges in the Fabrication of Biodegradable and Implantable Optical Fibers for Biomedical Applications
Next Article in Special Issue
Mechanical and Thermal Properties of Synthetic Polypropylene Fiber–Reinforced Renewable Oil Palm Shell Lightweight Concrete
Previous Article in Journal
High Magnesium and Sirolimus on Rabbit Vascular Cells—An In Vitro Proof of Concept
Previous Article in Special Issue
Investigation of the Residual Mechanical and Porosity Properties of Cement Mortar under Axial Stress during Heating
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Calculation of Dynamic Viscosity in Concentrated Cementitious Suspensions: Probabilistic Approximation and Bayesian Analysis

1
ETS de Ingenieros de Caminos, C. y P., Universidad de Castilla-La Mancha, Av. Camilo José Cela s/n, 13071 Ciudad Real, Spain
2
Real Academia de Ingeniería, Don Pedro 10, 28005 Madrid, Spain
3
Instituto de Cerámica y Vidrio (CSIC), C. Kelsen 5, Campus de Cantoblanco, 28049 Madrid, Spain
*
Author to whom correspondence should be addressed.
Materials 2021, 14(8), 1971; https://doi.org/10.3390/ma14081971
Submission received: 12 March 2021 / Revised: 5 April 2021 / Accepted: 11 April 2021 / Published: 14 April 2021
(This article belongs to the Special Issue Research of Mechanical Behavior of Cement and Concrete Composites)

Abstract

:
We present a new focus for the Krieger–Dougherty equation from a probabilistic point of view. This equation allows the calculation of dynamic viscosity in suspensions of various types, like cement paste and self-compacting mortar/concrete. The physical meaning of the parameters that intervene in the equation (maximum packing fraction of particles and intrinsic viscosity), together with the random nature associated with these systems, make the application of the Bayesian analysis desirable. This analysis permits the transformation of parametric-deterministic models into parametric-probabilistic models, which improves and enriches their results. The initial limitations of the Bayesian methods, due to their complexity, have been overcome by numerical methods (Markov Chain Monte Carlo and Gibbs Sampling) and the development of specific software (OpenBUGS). Here we use it to compute the probability density functions that intervene in the Krieger–Dougherty equation applied to the calculation of viscosity in several cement pastes, self-compacting mortars, and self-compacting concretes. The dynamic viscosity calculations made with the Bayesian distributions are significantly better than those made with the theoretical values.

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, η p , 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, w / c m [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:
η η 0 = 1 ϕ ϕ m [ η ] ϕ m
where,
  • η : dynamic viscosity of the suspension.
  • η 0 : dynamic viscosity of the fluid phase of the suspension.
  • ϕ : fraction in volume of the disperse solid phase of the suspension.
  • ϕ m : 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, ϕ m , 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, ϕ m 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 ϕ m 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 ϕ m , and equal to 2.5 for [ η ] (in monodisperse systems), and greater than 0.7 for ϕ m , and equal to 2.5 for [ η ] when polydispersion increases [1,13]. However, as the asymmetry of the particles increases, ϕ m 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 ( ϕ m 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 10 9 and 10 5 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 ϕ m 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 η 0 is not supposed to vary much), Equation (1) depends mainly on two parameters with physical significance: the maximum packing fraction of the particles, ϕ m , and the intrinsic viscosity of the system, [ η ] [6,10].
The maximum packing fraction of the disperse solid phase, ϕ m , is related to the particle size distribution and their shape [1,6]. Thus, in a system of monodisperse rigid spherical particles, the value of ϕ m 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 ϕ m . In polydisperse systems, the value of ϕ m increases with particle size variability as the space between them is filled more efficiently ( ϕ m ≈ 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 ϕ m 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, ( G , P ); the first, G , 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, P , 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:
f ( x 1 , x 2 , , x n ) = i = 1 n p ( x i | i ) ,
where i are the subsets of parents of the variables x i . 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, ϕ m 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:
η * = 1 ϕ p ϕ m p [ η ] p ϕ m p
where,
  • η * = η p η w : non-dimensional dynamic viscosity of the cement paste.
  • η p : dynamic viscosity of the cement paste.
  • η w : dynamic viscosity of the continuous fluid phase of the suspension (water).
  • ϕ p : fraction in volume of the disperse solid phase of the suspension (cementitious materials).
  • ϕ m p : maximum packing fraction of the particles of the disperse solid phase.
  • [ η ] p : 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 ϕ m p and [ η ] p 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, η p , 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:
η * N [ μ * , ν ]
μ * = 1 ϕ p ϕ m p [ η ] p ϕ m p
ϕ m p U [ ϕ m p min , ϕ m p max ]
η p U [ η ] p min , [ η ] p max
σ U [ σ min , σ max ]
where ν is an auxiliary variable depending on the standard deviation ( ν = σ 2 ).

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 η p , 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]:
η = 1 ϕ FA ϕ m FA [ η ] FA ϕ m FA
where,
  • η = η SCM η p : non-dimensional dynamic viscosity of self-compacting mortar.
  • η SCM : dynamic viscosity of self-compacting mortar.
  • η p : dynamic viscosity of cement paste.
  • ϕ FA : fraction in volume of the granular phase of the suspension (fine aggregate).
  • ϕ m FA : maximum packing fraction of particles of the granular phase.
  • [ η ] FA : 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:
η N [ μ , ν ]
μ = 1 ϕ FA ϕ m FA [ η ] FA ϕ m FA
ϕ m FA U [ ϕ m FA min , ϕ m FA max ]
η FA U [ η ] FA min , [ η ] FA max
σ U [ σ min , σ max ]

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]:
η = 1 ϕ FA ϕ m FA [ η ] FA ϕ m FA 1 ϕ CA ϕ m CA [ η ] CA ϕ m CA
where,
  • η = η SCC η p : non-dimensional dynamic viscosity of self-compacting concrete.
  • η SCC : dynamic viscosity of self-compacting concrete.
  • η p : dynamic viscosity of cement paste.
  • ϕ FA : fraction in volume of the finer granular phase of the suspension (fine aggregate).
  • ϕ m FA : maximum packing fraction of particles of the finer granular phase.
  • [ η ] FA : intrinsic viscosity of the system when adding the finer granular phase.
  • ϕ CA : fraction in volume of the thicker granular phase of the suspension (coarse aggregate).
  • ϕ m CA : maximum packing fraction of the thicker granular phase of the suspension (coarse aggregate).
  • [ η ] CA : 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:
η N [ μ , ν ]
μ = 1 ϕ FA ϕ m FA [ η ] FA ϕ m FA 1 ϕ CA ϕ m CA [ η ] CA ϕ m CA
ϕ m FA U [ ϕ m FA min , ϕ m FA max ]
ϕ m C A U [ ϕ m C A min , ϕ m C A max ]
η FA U [ η ] FA min , [ η ] FA max
η CA U [ η ] CA min , [ η ] CA max
σ U [ σ min , σ max ]

3. Materials and Methods

3.1. Cement Paste Suspensions

The cement pastes analyzed here belong to an experimental campaign recently published by De La Rosa et al. [5]. The cementitious suspensions were elaborated with two classes of Portland cement (c) and a mineral addition (ground granulated blast-furnace slag, GGBS). The particle size distribution of the cementitious materials (see Figure 1) was obtained by laser diffraction granulometry using an optical system Mastersizer 2000 (Malvern, UK). Two superplasticizer admixtures of a polymeric nature were used: one based on modified polycarboxylates in an aqueous base with a density equal to 1090 kg/m 3 and a dry solid residue of 40% (Sika ViscoCrete-20 HE), and the other of a poly-aryl-ether type with a density equal to 1058 kg/m 3 and a dry solid residue of 30% (BASF MasterEase-5025). The composition of the cement pastes in function of the different combinations of materials is the following:
  • 31 cement pastes with w / c m relationships equal to 0.35, 0.47, 0.53, and 0.63, each of which was composed of CEM I 52.5-SR, CEM II 32.5 B-L, 75% CEM I 52.5-SR, and 25% GGBS, 75% CEM II 32.5 B-L and 25% GGBS, respectively, and two relationships SP/ c m equal to 0.4 and 0.8% of a superplasticizer admixture based on polycarboxylates modified in an aqueous based (SikaViscoCrete-20 HE).
  • 8 cement pastes with w / c m relationships equal to 0.40, 0.47, 0.53, 0.63, made with CEM I 52.5-SR, and two relationships SP/ c m equal to 1.0 and 1.2%, with a superplasticizer admixture of a poly-aryl-ether type (BASF MasterEase-5025).
The mixture of the pastes was performed in the following manner: introduction of the cement materials and 90% of the water in the mixer and a 30 s rest. Then, mixing with a rotational velocity of the blades of 64 revolutions per minute (rpm) for 60 s. Subsequently, the mixing stopped, the material stuck to the sides of the recipient, the blades were scraped, and the superplasticizer was added with the remaining water. The mixing was restarted at a velocity of 92 rpm and after 90 s the process was stopped. This mixing protocol was adopted in accordance with the recommendations ASTM C305-99 and AASHTO T162-04 [36,37].
The dynamic viscosity of the cement pastes was measured using a rotational rheometer with a double cone-plate sensor Haake RS50 (Thermo Fischer), performing the tests at a constant temperature (25 C) with a water bath while controlling the shear rate and preventing the slippage between the suspension and the walls of the sensor. Measurements were performed twice with increasing and decreasing shear rate ramps from 0 to 600 s 1 for 3 min and dwelling time of 1 min at maximum rate. The samples were introduced in the rheometer and a first measurement was taken (considered as a pre-shear test). Subsequently, the sample was at rest for 60 s and then the test was repeated. The procedure was repeated again to have two measurements of η p for each sample (78 data).
The resulting flow curves are deemed that formed by the pairs of points shear stress-shear rate of the descent ramp. This permits determining the flow curve of cement suspensions (Figure 2). The slope of the adjusting straight line to the descending branch of the flow curve is η p , the value of which has been calculated in the range γ ˙ = 10–100 s 1 (note that the linear behavior of the descending branches extends beyond 100 s 1 , Figure 2). However, to make correct measurements in the typical range of application of cement pastes for fluid concrete applications (approximately γ ˙ = 100–200 s 1 ), when preparing the samples it is necessary to subject them to a higher shear rate to obtain a homogeneous and completely deflocculated mixture. Hence, the sample preparation shear rate range reached values of 600 s 1 .
The values of the descent branch of the flow curve were adjusted to a Bingham type rheological model in which the value of the slope of the descent branch is η p . These values of η p serve to adjust the parameters ϕ m p and [ η ] p of the Krieger–Dougherty equation, corresponding to cement pastes. Individual measurements η * (= η p / η w ) corresponding to each paste are included in Table A1 (Appendix A).

3.2. Self-Compacting Mortar/Concrete Suspensions

Next, eight self-compacting mortars from the research of Ouro et al. [38] were analyzed to evaluate the effect of adding the finest portion of a natural sand ( D m < 1.25 mm) on the parameters of the Krieger–Dougherty equation.
Besides, 17 self-compacting concretes from the research performed by Feys et al. [39], 16 self-compacting concretes from the research performed by Esmaeilkhanian et al. [40], and 9 self-compacting concretes from the doctoral thesis of Grünewald [41] were analyzed. The interest in using these self-compacting concretes lies in the nature and the size of their aggregates: crushed coarse aggregate was used (5–10 mm and 10–20 mm) by Feys et al.; round fine aggregate (0–2 mm) and crushed coarse aggregate (5–10 mm and 5–14 mm) were used by Esmaeilkhanian et al.; Grünewald employed both round fine (0.125–4 mm) and coarse (4–8 mm and 4–16 mm) aggregates.
The rheological measurement of self-compacting mortars was done using a parallel-plates rotational rheometer; every rheological measure of self-compacting concretes was done with rotational coaxial cylinder viscometers. However, none of the cited researches measured the value of the dynamic viscosity of the cement paste, so we have estimated these values from data of the scientific literature [6,7,42,43] according to Ghanbari et al. [6]. The composition of each of the mortars and concretes studied in this article is given in Appendix A.

4. Results and Discussion

The transformation of deterministic into probabilistic models was carried out using a Bayesian analysis methodology. For this purpose, the open-source software OpenBUGS was used which applies Markov Chain Monte Carlo and Gibbs Sampling to perform the simulations. In each model a total of 11,000 iterations were carried out to obtain the samples of the variables that can be interpreted as their density functions, which are the parameters of the deterministic models. We used all simulated samples but the first 1000, which belong to the burn-in stage. We did not experience any multimodal posterior density problem, which can appear in some cases, especially of extreme values or reliability analysis [44]. Besides, we did not consider noises since Bayesian methods deal with mixtures of the selected basic models instead of these models themselves, which provide them with more flexibility to reduce noise effects.

4.1. Ranges of the Parameters ϕ m i and [ η ] i to Be Used in the Bayesian Analysis of Cementitious Suspensions

4.1.1. Cement Paste Suspensions

From the rheometric tests performed on the cement pastes the values of the parameters ϕ m p and [ η ] p were calculated by means of the adjustment of the experimental results to the Krieger–Dougherty equation. Some of the adjustment curves may be seen in Figure 3. The results obtained from the statistical adjustment of the data are summarized in Table 1.
Previously, Struble et al. [1], and subsequently, Burgos-Montes et al. [3,45], carried out research on cement pastes of a different nature, that is, pastes with and without mineral additions and superplasticizer admixtures, in which they obtained good results from the adjustment of the experimental data to the Krieger–Dougherty equation. Struble et al. [1] concluded that in cement pastes the parameters ϕ m p and [ η ] p acquire adjustment values equal to 0.70 and 6.0, respectively. In turn, when the pastes are dispersed with a superplasticizer admixture, the value of ϕ m p varies between 0.64 and 0.80, while the values of [ η ] p range between 4.5 and 6.0, with average values of 0.70 and 5.0 for ϕ m p and [ η ] p , respectively (Table 2).
In light of the results obtained, it may be verified that the ranges of the values experimentally determined for ϕ m p and [ η ] p are very similar to those contributed by Struble et al. [1]. The differences that appear, basically in the lower value of the parameter ϕ m p , may be due to the physical–chemical interactions arising between the cement particles as a consequence of the chemical nature of the superplasticizer admixtures used in each case.

4.1.2. Self-Compacting Mortar/Concrete Suspensions

To establish the domains of validity of the parameters of the Krieger–Dougherty equation in each phase of self-compacting mortar/concrete, the next hypotheses have been followed. First, we have as reference the criterium established by Abo-Dhaheer et al. [7], who consider rigid spherical particles for which ϕ m increases with the addition of solid phases according to theoretical packing values ( ϕ m = 0.524 for powder phase, ϕ m = 0.63 for powder plus fine aggregate phase, and ϕ m = 0.74 for powder plus fine and coarse aggregate phase). Also, according to Hu et al. [27], ϕ m could de estimated under the same hypothesis as:
ϕ m = 1 0.45 d m D m 0.19
where d m and D m are the minimum and maximum dimension of the particles, respectively. Having this into account, we select the lower limit of ϕ m , equal to 0.55, for every phase of self-compacting concrete (powder, powder plus fine aggregate, powder plus fine and coarse aggregate).
To determine the upper limit of ϕ m FA , we consider the experimental maximum packing fraction measured by Grünewald for fine round aggregate [41]. Furthermore, the experimental maximum packing fraction measured for both fine and coarse round aggregates is 0.809 (40% fine aggregate, 60% coarse aggregate) [41]. Zentar et al. [46] indicate that ϕ m is lower with round aggregates ( ϕ m = 0.793) than with crushed aggregates ( ϕ m = 0.901). This is explained by the lubricating effect of the powder and cementitious phases, which is more important with crushed aggregates [46]. Taking into account this, and to establish an adequate value of reference, we take as upper limit of ϕ m = 0.894 (packing density of all solid particles). This value has been extracted from the investigation of Kwan et al. [47].
The morphology of the particles in any suspension is mainly controlled for the intrinsic viscosity [46]. Its value is lower for round aggregates than for crushed aggregates because these are closer to the spherical shape ( [ η ] = 2.5). The following equation, proposed by Pabst et al. [28], could be used to estimate [ η ] knowing the aspect ratio of the particles.
[ η ] = 2.5 + 0.123 D m d m 1 0.925
Also, there are other relationships between circularity of particle and intrinsic viscosity [30,31]. As the shape of the set of particles in concrete is broad, we have a lot of uncertainty regarding the maximum and minimum values of [ η ] that can be acquired in each phase, so we select the same interval in each one of them. Thus, we select a lower limit of [ η ] = 2.5 and an upper limit [ η ] = 9, taking into account these references.
Abo-Daheer et al. [7] consider the ϕ m [ η ] product (exponent of the Krieger–Dougherty equation) practically constant and equal on average to 1.9 in every phase, since both parameters have an inverse dependence on the shear rate [6,7]. This hypothesis is very close to other models similar to Krieger–Dougherty, such as that of Quemada (exponent equal to 2) [11]. However, in this parameter product, there is no information about the shape, morphology, particle size distribution, etc., of concrete raw materials.
The ranges of values for the parameters in the mortar phase can be extracted from this analysis. Actually, a lower value of ϕ m FA in the mortar phase can be calculated through Equation (23), considering monosize sand particles shape like spheres ( ϕ m = 0.55). The upper value of ϕ m FA is selected according to the experimental maximum packing fraction measured by Grünewald for round fine aggregate with a certain grade of polydispersion [41], which is equal to 0.717. Similarly, we select [ η ] FA within the 2.5 to 9 interval following the investigations of Choi and Szecsy [30,31].

4.2. Bayesian Analysis of the Krieger–Dougherty Equation in Cementitious Suspensions

4.2.1. Cement Paste Suspensions

Table 3 presents the intervals of the values detected for the parameters ϕ m p and [ η ] p , which appear to establish the validity domains in order to carry out the Bayesian analysis of the Krieger–Dougherty equation in cement pastes. The ranges in this work are calculated from the experimental measurements performed in [5].
The description of the model proposed from the Krieger–Dougherty equation in cement pastes and the domains of the definition of the parameters are shown below:
η * N [ μ * , ν ]
μ * = 1 ϕ p ϕ m p [ η ] p ϕ m p
ϕ m p U [ 0.510 , 0.830 ]
η p U [ 4.30 , 6.80 ]
σ U [ 0 , 400 ]
Figure 4 represents the graph of the Bayesian network, obtained for the model, which describes the hierarchical dependence structure of the collection of variables involved. As indicated, the graph of the model permits verifying that ϕ m p , [ η ] p and σ (standard deviation) are independent variables.
Table 4 summarizes the statistics calculated after performing the Bayesian analysis. The graphs of the density functions of the parameters of the Bayesian model are shown in Figure 5. They demonstrate the variability in the distribution of the values of the parameters ϕ m p and [ η ] p , which is indicative of the properties and the nature of these cement systems. Thus, the amplitude of the density function of the parameter ϕ m p describes the poly-dispersion insofar as the size of the particles within the system. Similarly, the characteristics of the density function of the parameter [ η ] p represents the variability of shapes of the particles present in the suspension. Figure 6a corresponds to the density function of the product of parameters ϕ m p [ η ] p (exponent of the Krieger–Dougherty equation). This function indicates the most likely value that the exponent of the equation in cement pastes can acquire. It reflects what is collected in the scientific literature [1], but now there is a probability density function and not a single value. Figure 6b shows the bivariate histogram of the parameters ϕ m p and [ η ] p .

4.2.2. Self-Compacting Mortar Suspensions

The Krieger–Dougherty model defined for the Bayesian analysis of self-compacting mortars (SCM) from the research of Ouro et al. [38] and the parameter definition domains are as follows:
η N [ μ , ν ]
μ = 1 ϕ FA ϕ m FA [ η ] FA ϕ m FA
ϕ m FA U [ 0.550 , 0.717 ]
η FA U [ 2.5 , 9.0 ]
σ U [ 0 , 400 ]
Figure 7 is the graph of the Bayesian network of the model. The statistics calculated from the parameters after performing the analysis are included in Table 5. In this case, mortars were manufactured with the finest portion of a natural sand ( D m ≤ 1.25 mm). The shape of the particles is very close to the sphere and we can consider this mortar phase like a monodisperse system, so the calculated values with the Bayesian analysis of the parameters ϕ m FA and [ η ] FA (Table 5) are very close to the theoretical values proposed by Abo-Daheer et al. (which are 0.63 and 2.5, respectively) [7]. Also, the exponent of the Krieger–Dougherty equation ϕ m FA [ η ] FA is 1.8, which is very approximate to the theoretical value proposed by Abo-Daheer et al. (which is 1.9) [7].
The graphics of density functions of the parameters are shown in Figure 8. Both ϕ m FA and, specially, [ η ] FA display a very clear peak of probability in a narrow range of values in the mortar. This confirms the hypothesis that the shape of the particles is very close to the sphere and we can consider this mortar phase as a monodisperse system. Thereby, this values are well known in the case of monodisperse system of spheres ( ϕ m FA 0.63, [ η ] FA 2.5).
To appreciate the improvement made with the Bayesian analysis better, we can use the mean values of the obtained parameters (Table 5) in the Krieger–Dougherty model to predict the actual experimental data of Ouro et al. [38]. Next, we compare them with the results obtained using the theoretical values proposed by Abo-Daheer et al. ( ϕ m FA = 0.63 for powder plus fine aggregate phase, ϕ m i [ η ] i = 1.9) [7]. Table 6 shows all the results and the errors obtained with both approaches, expressed as a percentage of the experimental values. The mean error using the theoretical values for the parameters is 43% (standard deviation 39%), which gets reduced to 33% (22%) when the mean values from the Bayesian analysis are used. It should be noted that these errors being high, they are reasonable values according to the uncertainties present in the system and similar to the acceptable error chosen by Ghanbari et al. [6] for their model (25%).

4.2.3. Self-Compacting Concrete Suspensions

The Krieger–Dougherty model defined for the Bayesian analysis of self-compacting concretes (SCC) from the research of Feys et al. [39], Esmaeilkhanian et al. [40], and Grünewald [41]; the parameter definition domains are as follows:
η N [ μ , ν ]
μ = 1 ϕ FA ϕ m FA [ η ] FA ϕ m FA 1 ϕ CA ϕ m CA [ η ] CA ϕ m CA
ϕ m FA U [ 0.550 , 0.717 ]
ϕ m CA U [ 0.550 , 0.894 ]
η FA U [ 2.5 , 9.0 ]
η CA U [ 2.5 , 9.0 ]
σ U [ 0 , 400 ]
The scheme of the Bayesian network of the model, which describes the hierarchy and independence of variables, is represented in Figure 9. The statistics calculated from the parameters of the model after performing the analysis are included in Table 7. The graphics of the functions of density of the parameters of every phase composing SCC and the density function of the product of parameters ϕ m i [ η ] i (exponent of the Krieger–Dougherty equation) are shown in Figure 10 (Feys et al. [39]), and Figure 11 (Esmaeilkhanian et al. [40]). Figure 12 represents the bivariate histogram of the parameters ϕ m i and [ η ] i in every phase of the research of Feys et al. [39] and Esmaeilkhanian et al. [40].
The sands of the research of Esmaeilkhanian et al. and Grünewald are round aggregates. The mean value of ϕ m FA is similar (≈0.68). However, this parameter has a lower value for the sand of Feys et al., which can indicate that is a crushed sand. Regarding [ η ] FA , the mean value calculated for the fine aggregate data of Feys et al. is higher than the value obtained in the data of Esmaeilkhanian et al. This could reinforce the hypothesis of crushed fine aggregate of Feys et al. The value obtained for [ η ] FA in Grünewald data is the highest (≈4.03). This could indicate that the shape of this fine aggregate has a lower circularity than the other sands [30,31].
The mean value of ϕ m CA is similar in the three coarse aggregates (≈0.74). We know that the self-compacting concretes of the research of Feys et al. and Esmaeilkhanian et al. are crushed aggregates whereas the coarse aggregate of Grünewald is rounded. This could be explained by the lubricating effect of the powder and cementitious phases, which is more important with crushed aggregates [46]. Regarding [ η ] CA the mean values obtained in every research are different and higher than 2.5. This means that the shape of the coarse aggregates moves away from the spherical shape, since the parameter [ η ] CA is related with the circularity of the aggregate [30,31]. This is plausible for the crushed coarse aggregates employed by Feys et al. and Esmaeilkhanian et al. Also, this is possible for elongated round aggregates. In this particular case, the investigation of Grünewald indicates that sand grains are much more rounded than coarse aggregates [41]. This fact is reflected in the Bayesian analysis of the data as well.
If we look at the non-parametric density functions, ϕ m FA displays a peak of probability for the fine aggregate phase of Esmaeilkhanian et al. whereas the density function is smoother for Feys et al. (Figure 10a,b and Figure 11a,b). This fact indicates that the maximum packing of the fine aggregate gets a value with high probability, which could be due to the use of fine round aggregate (the maximum packing value of spheres is well known). The plateau of values with similar probability for the fine aggregate of Feys et al. confirms that crushed fine aggregate was used in this research.
Respect to the coarse aggregate phase, in both cases non-parametric density functions of ϕ m CA and [ η ] CA (Figure 10c,d and Figure 11c,d) show a shape similar to a uniform density function with similar ranges of values and probability. Both investigations use crushed coarse aggregate; this indicates that there exists a variety of shapes of coarse aggregates leading to a significant dispersion in the values of ϕ m CA .
Figure 10e,f and Figure 11e,f show that the maximum probability of the exponent of the Krieger–Dougherty equation is, approximately, equal to the theoretical value of 1.9 in the fine granular phase. However, the width of the density function in the coarse aggregate phase indicates that the exponent could acquire an ample range of values with similar probability. Figure 12 shows the bivariate histogram of the parameters ϕ m i and [ η ] i in every phase composing the self-compacting concretes of Esmaeilkhanian et al. and Feys et al.
The improvement made with the Bayesian analysis can be visualized by carrying out the same process that has been done with mortars, applying the mean values of the parameters obtained in the Bayesian analysis (Table 7) of the Krieger–Dougherty model to the data of Feys et al. [39] (see the results in Table 8), Esmaeilkhanian et al. [40] (Table 9), and Grünewald [41] (Table 10), and compare them with the results obtained with the theoretical values proposed by Abo-Daheer et al. [7] ( ϕ m FA = 0.63 for powder plus fine aggregate phase; ϕ m CA = 0.74 for powder plus fine and coarse aggregate phase; ϕ m i [ η ] i = 1.9).
Again, the error given by the Krieger–Dougherty model in calculating the dynamic viscosity of the SCC is smaller in the three series of tests. Namely, the mean error goes from 77% (standard deviation of 8%) to 25% (24%) in the Feys et al. series [39], from 42% (24%) to 36% (26%) in the Esmaeilkhanian et al. tests [40] and from 71% (9%) to only 17% (22%) in the series reported by Grünewald [41]. In this last series of calculations, it is noteworthy that the error of most of the predicted values (8 out of 9) is smaller than 25%, i.e., the predictions are excellent according to the criterion established by Ghanbari et al. [6].

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, ϕ m i and their corresponding intrinsic viscosities [ η ] i ) 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 [ η ] FA 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 ϕ m FA 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 [ η ] FA . 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 [ η ] CA 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.

Author Contributions

Conceptualization, Á.D.L.R. and E.C.; methodology, all authors; simulations, Á.D.L.R.; validation, Á.D.L.R. and G.R.; resources: G.R and R.M.; writing—original draft preparation, Á.D.L.R.; writing—review and editing, all authors; supervision, G.R.; funding acquisition, G.R. All authors have read and agreed to the published version of the manuscript.

Funding

The authors appreciate the funding received from the Ministerio de Innovación, Ciencia y Universidades, Spain, through PID2019-110928RB-C31 and RTC-2017-6736-3 research projects, and from Junta de Comunidades de Castilla-La Mancha, Spain, through SBPLY/19/180501/000220 research project. Ángel De La Rosa also acknowledges the FPI BES-2016-077458 scholarship from the Ministerio de Innovación, Ciencia y Universidades, Spain.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Acknowledgments

We thank Cementos Portland Valderrivas, Cementos Tudela Veguín, SIKA, and BASF for the complimentary supply of materials.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
CIConfidence interval
dDiameter of the smallest particle in the system
DDiameter of the largest particle in the system
d m Minimum dimension of granular particle
D m Maximum dimension of granular particle
GGamma probability function
GGBSGround granulated blast-furnace slag
iNumber of nodes
nNumber of conditional probability densities
NNormal probability function
p ( x i | π i ) Probability of x i conditioned to π i
rpmNumber of revolutions per minute
SCCSelf-compacting concrete
SCMSelf-compacting mortar
SCSFRCSelf-compacting steel–fiber reinforced concrete
SP/ c m Superplasticizer–cementitious materials ratio
UUniform probability function
w / c m Water–cementitious materials ratio
ε * Residual error for non-dimensional dynamic viscosity of cement paste
ε Residual error for non-dimensional dynamic viscosity of self-compacting concrete
ε Residual error for non-dimensional dynamic viscosity of self-compacting mortar
γ ˙ Shear strain rate
η Dynamic viscosity
η SCC Dynamic viscosity of self-compacting concrete
η SCM Dynamic viscosity of self-compacting mortar
η p Dynamic viscosity of cement paste
η w Dynamic viscosity of water
η 0 Dynamic viscosity of the fluid phase in the suspension
η * Non-dimensional dynamic viscosity of cement paste
η Non-dimensional dynamic viscosity of self-compacting concrete
η Non-dimensional dynamic viscosity of self-compacting mortar
[ η ] Intrinsic viscosity
[ η ] CA Intrinsic viscosity of the coarse aggregate phase in self-compacting concrete
[ η ] FA Intrinsic viscosity of the fine aggregate phase in self-compacting mortar/concrete
[ η ] i Intrinsic viscosity of particles in i phase (i can be equal to p for cement paste,
FA for fine aggregate or CA for coarse aggregate)
[ η ] p Intrinsic viscosity of cement paste
μ * Mean value for non-dimensional dynamic viscosity of cement paste
μ Mean value for non-dimensional dynamic viscosity of self-compacting concrete
μ Mean value for non-dimensional dynamic viscosity of self– compacting mortar
ν = 1 σ 2 Auxiliary variable for the model of probability
π i Set of nodes X i in G
σ Standard deviation
ϕ Volume fraction of particles
ϕ CA Volume fraction of coarse aggregate
ϕ FA Volume fraction of fine aggregate
ϕ m Maximum packing fraction of particles
ϕ m CA Maximum packing fraction of particles in the coarse aggregate phase in SCC
ϕ m FA Maximum packing fraction of particles in the fine aggregate phase in SCM/SCC
ϕ m i Maximum packing fraction of particles in i phase (i can be equal to p for cement
paste, FA for fine aggregate or CA for coarse aggregate)
ϕ m p Maximum packing fraction of particles in cement paste
ϕ p Volume fraction of particles in cement paste
G Acyclic graph directed
P Associated joint probability density of all nodes
X Nodes or random variables
τ Shear stress

Appendix A. Data Annex

Appendix A.1. Experimental Dynamic Viscosities of the Cement Pastes

Table A1. Experimental dynamic viscosities of the cement pastes.
Table A1. Experimental dynamic viscosities of the cement pastes.
Paste η p [mPa s] ϕ p SP / cm η * Paste η p [mPa s] ϕ p SP / cm η *
CEM I 52.5 N-SR180.3360.419.3
380.3760.440.7CEM II 32.5 BL-II360.4140.438.6
510.4040.454.7130.3860.413.9
190.3360.420.4320.4140.834.3
370.3760.439.7230.3860.824.7
580.4040.462.2140.3450.815.0
260.3760.827.9250.3860.826.8
430.4040.846.1
200.3360.821.4 460.3990.449.3
300.3760.832.2 310.3710.433.2
470.4040.850.4 100.3320.410.7
130.3361.014.6 490.3990.452.5
320.3761.035.975% CEM I 52.5 N-SR340.3710.436.4
470.4041.052.8+ 25% GGBS400.3990.842.9
840.4431.094.3 270.3710.828.9
140.3361.015.7 110.3320.811.8
340.3761.038.2 390.3990.841.8
540.4041.060.6 290.3710.831.1
910.4431.0102.1
300.3761.233.7 340.4050.436.4
450.4041.250.575% CEM II 32.5 BL-II280.3770.430.0
890.4431.299.9+ 25% GGBS120.3380.412.9
150.3361.216.8 350.4050.837.5
480.4041.253.9 090.3380.89.6
860.4431.296.5

Appendix A.2. Dosages and Dynamic Viscosities of the Self-Compacting Mortars

Table A2. Dosages and dynamic viscosities of the self-compacting mortars in Ouro et al. [38].
Table A2. Dosages and dynamic viscosities of the self-compacting mortars in Ouro et al. [38].
NameCEM Portland + Fly Ash + Silica Fume [kg/m 3 ]Water [kg/m 3 ]SP [kg/m 3 ]Fine Agg. [kg/m 3 ] η SCM [Pa s] η p [Pa s]
N3, N48423357.968341.900.331
N7, N87923577.488181.100.286
N117843297.418861.660.313
N138533558.067741.410.313
N157383486.978930.940.268
E3, E46093536.329790.500.169
E75573815.789950.350.079
E135893716.119530.370.124

Appendix A.3. Dosages and Dynamic Viscosities of the Self-Compacting Concretes

Table A3. Dosages and dynamic viscosities of the self-compacting concretes in Feys et al. [39].
Table A3. Dosages and dynamic viscosities of the self-compacting concretes in Feys et al. [39].
NameCEM Portland + 8% SF (GUbSF)Water [kg/m 3 ]SP L + SP S [kg/m 3 ]Fine Agg. [kg/m 3 ]Coarse Agg. [kg/m 3 ] η SCC [Pa s] η p [Pa s]
SCC159918510.0 + 1.486374150.80.413
SCC260218310.0 + 1.785972642.60.422
SCC360518310.0 + 1.585672938.00.422
SCC459718610.0 + 1.687175241.40.413
SCC759417510.0 + 2.084973767.50.431
SCC864516515.3 + 2.484075059.00.458
SCC95581975.3 + 0.185772828.00.377
SCC1060118010.0 + 1.581676045.00.422
SCC1164419010.0 + 1.480071735.00.422
SCC1256217110.0 + 2.887977696.50.422
SCC1360317810.0 + 1.384572441.50.422
SCC1459817710.0 + 0.084373729.30.422
SCC1559718210.0 + 1.883872449.60.422
SCC1660217910.0 + 1.879778955.60.422
SCC1759617710.0 + 1.674383144.80.422
SCC1860217910.0 + 2.285173071.20.422
SCC1968115620.0 + 4.1855782155.80.485
Table A4. Dosages and dynamic viscosities of the self-compacting concretes in Esmaeilkhanian et al. [40].
Table A4. Dosages and dynamic viscosities of the self-compacting concretes in Esmaeilkhanian et al. [40].
NameCEM Port. (GU) [kg/m 3 ]Water [kg/m 3 ]SF [kg/m 3 ]SP 1 + SP 2 [kg/m 3 ]Fly Ash [kg/m 3 ]Fine Agg. [kg/m 3 ]Coarse Agg. [kg/m 3 ] η SCC [Pa s] η p [Pa s]
SCC1377165132.55 + 2.1513095071559.30.404
SCC2344181121.96 + 1.7211995071529.00.350
SCC3377165132.18 + 1.9213095071569.50.404
SCC438116502.45 + 21.5131950715124.00.413
SCC5320181282.13 + 1.8711695071525.00.340
SCC7315166111.88 + 1.6610999374962.00.350
SCC8365191.5132.18 + 1.9212691769125.00.350
SCC9308169.5111.92 + 2.02106101076028.00.331
SCC10341153.5121.92 + 2.95118101076072.00.395
SCC11374138133.05 + 3.601291010760128.00.449
SCC12341153.5122.80 + 2.85118101076071.00.395
SCC13341153.5122.90 + 2.95118101076067.00.395
SCC14308169.5112.12 + 1.73106100575635.00.331
SCC16308169.5112.12 + 2.18106100575637.00.331
SCC17308169.5112.12 + 2.18106100575639.00.331
SCC18308169.5112.12 + 2.18106100575930.00.331
Table A5. Dosages and dynamic viscosities of the self-compacting concretes in Grünewald [41].
Table A5. Dosages and dynamic viscosities of the self-compacting concretes in Grünewald [41].
NameCEM I 52.5 R [kg/m 3 ]CEM III 42.5 N [kg/m 3 ]Water [kg/m 3 ]SP LR + SP HR [kg/m 3 ]Fly Ash [kg/m 3 ]Fine Agg. [kg/m 3 ]Coarse Agg. [kg/m 3 ] η SCC [Pa s] η p [Pa s]
OS12491551722.58 + 1.5814291368269.20.404
OS22631491812.88 + 1.4417387665559.40.413
OS32491491712.59 + 2.12146108950887.90.413
OS42691431812.78 + 1.85173104548756.00.413
OS503351552.10 + 1.26168113452897.60.413
OS603521642.10 + 1.18192108950881.00.422
OS703671732.17 + 1.09217104548762.20.422
OS82281511812.68 + 1.49166110046771.30.395
OS92461641882.73 + 1.31180105844957.50.404

References

  1. Struble, L.; Sun, G. Viscosity of Portland cement pastes as a function of concentration. Adv. Cem. Based Mater. 1995, 2, 62–69. [Google Scholar] [CrossRef]
  2. Krieger, I.; Dougherthy, T. A mechanism for non-Newtonian flow in suspensions of rigid spheres. J. Rheol. 1959, 3, 137–152. [Google Scholar] [CrossRef]
  3. Burgos-Montes, O.; Alonso, M.; Puertas, F. Viscosity and water demand of limestone and fly ash–blended cement pastes in the presence of superplasticisers. Constr. Build. Mater. 2013, 48, 417–423. [Google Scholar] [CrossRef]
  4. De La Rosa, Á.; Poveda, E.; Ruiz, G.; Cifuentes, H. Proportioning of self-compacting steel-fiber reinforced concrete mixes based on target plastic viscosity and compressive strength: Mix-design procedure and experimental validation. Constr. Build. Mater. 2018, 189, 409–419. [Google Scholar] [CrossRef]
  5. De La Rosa, Á.; Poveda, E.; Ruiz, G.; Moreno, R.; Cifuentes, H.; Garijo, L. Determination of the plastic viscosity of superplasticized cement pastes through capillary viscometers. Constr. Build. Mater. 2020, 260, 119715. [Google Scholar] [CrossRef]
  6. Ghanbari, A.; Karihaloo, B.L. Prediction of the plastic viscosity of self-compacting steel fibre reinforced concrete. Cem. Concr. Res. 2009, 39, 1209–1216. [Google Scholar] [CrossRef]
  7. Abo Dhaheer, M.; Al-Rubaye, M.; Alyhya, W.S.; Karihaloo, B.; Kulasegaram, S. Proportioning of self–compacting concrete mixes based on target plastic viscosity and compressive strength: Part 1–Mix design procedure. J. Sustain. Cem. Based Mater. 2015, 5, 199–216. [Google Scholar] [CrossRef]
  8. De Kruif, C.; van Iersel, E.; Vrij, A.; Russel, W. Hard sphere colloidal dispersions: Viscosity as a function of shear rate and volume fraction. J. Chem. Phys. 1985, 83, 4717–4725. [Google Scholar] [CrossRef]
  9. Hinch, E.; Leal, L. The effect of Brownian motion on the rheological properties of a suspension of non-spherical particles. J. Fluid Mech. 1972, 52, 683–712. [Google Scholar] [CrossRef] [Green Version]
  10. Shewan, H.; Stokes, J. Analytically predicting the viscosity of hard sphere suspensions from the particle size distribution. J. Non-Newtonian Fluid Mech. 2015, 222, 72–81. [Google Scholar] [CrossRef]
  11. Quemada, D. Rheology of concentrated disperse systems and minimum energy–dissipation principle. 1 Viscosity–concentration relationship. Rheol. Acta 1977, 16, 82–94. [Google Scholar] [CrossRef]
  12. Mendoza, C.; Santamaria-Holek, I. The rheology of hard sphere suspensions at arbitrary volume fractions: An improved differential viscosity model. J. Chem. Phys. 2019, 130, 044904. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Barnes, H.; Hutton, J.; Walters, K. An Introduction to Rheology; Elsevier: Amsterdam, The Netherlands, 1989. [Google Scholar]
  14. Castillo, E.; Hadi, A.; Gutiérrez, J. Sistemas Expertos y Modelos de Redes Probabilísticas; Academia de Ingeniería: Madrid, Spain, 1996. [Google Scholar]
  15. Castillo, E.; Gutiérrez, J.; Hadi, A. Expert Systems and Probabilistic Network Models; Springer-Verlag: New York, NY, USA, 1997. [Google Scholar]
  16. Gelfand, A.E.; Smith, A.F.M. Sampling-based approaches to calculating marginal densities. J. Am. Stat. Assoc. 1990, 85, 398–409. [Google Scholar] [CrossRef]
  17. Lunn, D.; Spiegelhalter, D.; Thomas, A.; Best, N. The BUGS project: Evolution, critique and future directions (with discussion). Stat. Med. 2009, 28, 3049–3082. [Google Scholar] [CrossRef] [PubMed]
  18. Cowles, M. Applied Bayesian Statistics with R and OpenBUGS. Examples; Springer: New York, NY, USA, 2013. [Google Scholar]
  19. Castillo, E. Bayesian Methods with OpenBUGS. 2020. Available online: https://www.uclm.es/conocimiento/cursos/bayesian-methods-openbugs-en (accessed on 12 April 2021).
  20. Cavero, P. Self-compacting concrete used in the construction of the arch of the Almonte Viaduct. Hormigón y Acero 2020, 71, 87–97. [Google Scholar]
  21. Castillo, E.; Menéndez, J.; Sánchez-Cambronero, S. Predicting traffic flow using Bayesian networks. Transp. Res. Part B 2008, 42, 482–509. [Google Scholar] [CrossRef]
  22. Moreno, R. Reología de Suspensiones Cerámicas; Consejo Superior de Investigaciones Científicas: Madrid, Spain, 2005. [Google Scholar]
  23. Domènech, X.; Peral, J. Química Ambiental de Sistemas Terrestres; S.A. Reverté: Barcelona, Spain, 2006. [Google Scholar]
  24. Moreno, R. Better ceramics through colloid chemistry. J. Eur. Ceram. Soc. 2020, 40, 559–587. [Google Scholar] [CrossRef]
  25. Skalny, J. Materials Science of Concrete 1; U.S. Department of Energy Office of Scientific and Technical Information: Oak Ridge, TN, USA, 1989.
  26. Toutou, Z.; Roussel, N. Multi scale experimental study of concrete rheology: From water scale to gravel scale. Mater. Struct. 2006, 39, 189–199. [Google Scholar] [CrossRef]
  27. Hu, C.; de Larrard, F. The rheology of fresh high-performance concrete. Cem. Concr. Res. 1996, 26, 283–294. [Google Scholar] [CrossRef]
  28. Pabst, W.; Gregorova, E.; Berthold, C. Particle shape and suspension rheology of short-fiber systems. J. Eur. Ceram. Soc. 2006, 26, 149–160. [Google Scholar] [CrossRef]
  29. Brenner, H. Rheology of a dilute suspension of axisymmetric brownian particles. Int. J. Multiph. Flow 1974, 1, 195–341. [Google Scholar] [CrossRef]
  30. Choi, M. Numerical prediction on the effects of the coarse aggregate size to the pipe flow of pumped concrete. J. Adv. Concr. Technol. 2014, 12, 239–249. [Google Scholar] [CrossRef] [Green Version]
  31. Szecsy, R. Concrete Rheology. Ph.D. Thesis, University of Illinois at Urbana-Champaign, Urbana, IL, USA, 1997. [Google Scholar]
  32. Batchelor, G. The stress generated in a non–dilute suspension of elongated particles by pure straining motion. J. Fluid Mech. Digit. Arch. 1971, 46, 813–829. [Google Scholar] [CrossRef]
  33. Maron, S.; Pierce, P. Application of Ree-Eyring generalized flow theory to suspensions of spherical particles. J. Colloids Sci. 1956, 11, 80–95. [Google Scholar] [CrossRef]
  34. Castillo, E.; Menéndez, J.; Sánchez-Cambronero, S.; Calviño, A.; Sarabia, J. A hierarchical optimization problem: Estimating traffic flow using gamma random variables in a Bayesian context. Comput. Oper. Res. 2014, 41, 240–251. [Google Scholar] [CrossRef]
  35. OpenBUGS. 2009. Available online: www.openbugs.net (accessed on 12 April 2021).
  36. ASTM International. ASTM C305-99, Standard Method of Test for Mechanical Mixing of Hydraulic Cement Pastes and Mortars of Plastic Consistency; ASTM International: West Conshohocken, PA, USA, 2004. [Google Scholar]
  37. ASTM International. AASHTO T162-04, Standard Method of Test for Mechanical Mixing of Hydraulic Cement Pastes and Mortars of Plastic Consistency; ASTM International: West Conshohocken, PA, USA, 2004. [Google Scholar]
  38. Koura, B.O.; Hosseinpoor, M.; Yahia, A.; Kadri, E.; Kaci, A. A new proportioning approach of low and normal binder self-consolidating concrete based on the characteristics of fine mortar and granular skeleton. Constr. Build. Mater. 2020, 239, 117892. [Google Scholar] [CrossRef]
  39. Feys, D.; Khayat, K.; Pérez-Schell, A.; Khatib, R. Prediction of pumping pressure by means of new tribometer for highly-workable concrete. Cem. Concr. Compos. 2015, 57, 102–115. [Google Scholar] [CrossRef]
  40. Esmaeilkhanian, B.; Khayat, K.; Yahia, A.; Feys, D. Effects of mix design parameters and rheological properties on dynamic stability of self-consolidating concrete. Cem. Concr. Compos. 2014, 54, 21–28. [Google Scholar] [CrossRef]
  41. Grünewald, S. Performance-Based Design of Self-Compacting Fibre Reinforced Concrete. Ph.D. Thesis, Technische Universiteit Darmstadt, Delft, The Netherlands, 2004. [Google Scholar]
  42. Sun, Z.; Voigt, T.; Shah, S. Rheometric and ultrasonic investigations of viscoelastic properties of fresh portland cement pastes. Cem. Concr. Res. 2006, 36, 278–287. [Google Scholar] [CrossRef]
  43. Nehdi, M.; Rahman, M. Estimating rheological properties of cement pastes using various rheological models for different test geometry, gap and surface friction. Cem. Concr. Res. 2004, 34, 1993–2007. [Google Scholar] [CrossRef]
  44. Sundar, V.S.; Manohar, C.S. Updating reliability models of statically loaded instrumented structures. Struct. Saf. 2013, 41, 21–30. [Google Scholar] [CrossRef]
  45. Vikan, H. Rheology and Reactivity of Cementitious Binders with Plasticizers. Ph.D. Thesis, Department of Materials Science and Engineering, Norwegian University of Science and Technology, Trondheim, Norway, 2005. [Google Scholar]
  46. Zentar, R.; Bala, M.; Boustingorry, P. Influence of the granular mixture properties on the rheological properties of concrete: Yield stress determination using modified Chateau et al. model. Int. J. Archit. Civ. Constr. Sci. 2019, 13, 656–660. [Google Scholar]
  47. Kwan, A.; Li, L. Combined effects of water film, paste film and mortar film thicknesses on fresh properties of concrete. Constr. Build. Mater. 2014, 50, 598–608. [Google Scholar] [CrossRef]
Figure 1. Particle size distribution of the cementitious materials.
Figure 1. Particle size distribution of the cementitious materials.
Materials 14 01971 g001
Figure 2. Descending branches of the flow curves of cement pastes obtained with the rotational double cone-plate rheometer: (a) 75% CEM I 52.5 N-SR + 25% ground granulated blast-furnace slag (GGBS); (b) CEM I 52.5 N-SR.
Figure 2. Descending branches of the flow curves of cement pastes obtained with the rotational double cone-plate rheometer: (a) 75% CEM I 52.5 N-SR + 25% ground granulated blast-furnace slag (GGBS); (b) CEM I 52.5 N-SR.
Materials 14 01971 g002
Figure 3. Adjustment curves for the experimental measures to the Krieger–Dougherty equation: (a) all the pastes designed; (b) CEM I 52.5 N-SR.
Figure 3. Adjustment curves for the experimental measures to the Krieger–Dougherty equation: (a) all the pastes designed; (b) CEM I 52.5 N-SR.
Materials 14 01971 g003
Figure 4. Graph of the Bayesian network of the Krieger–Dougherty equation in cement pastes.
Figure 4. Graph of the Bayesian network of the Krieger–Dougherty equation in cement pastes.
Materials 14 01971 g004
Figure 5. Density functions of probability of the parameters of the cement paste suspensions: (a) [ η ] p ; (b) ϕ m p .
Figure 5. Density functions of probability of the parameters of the cement paste suspensions: (a) [ η ] p ; (b) ϕ m p .
Materials 14 01971 g005
Figure 6. (a) Density function of probability of the product of parameters ϕ m p [ η ] p (exponent of the Krieger–Dougherty equation); (b) bivariate histogram of ϕ m p and [ η ] p .
Figure 6. (a) Density function of probability of the product of parameters ϕ m p [ η ] p (exponent of the Krieger–Dougherty equation); (b) bivariate histogram of ϕ m p and [ η ] p .
Materials 14 01971 g006
Figure 7. Graph of the Bayesian network of the Krieger–Dougherty equation in self-compacting mortar.
Figure 7. Graph of the Bayesian network of the Krieger–Dougherty equation in self-compacting mortar.
Materials 14 01971 g007
Figure 8. Probability density functions of the parameters of the fine granular phase (FA) in self-compacting mortar: (a) ϕ m FA ; (b) [ η ] FA .
Figure 8. Probability density functions of the parameters of the fine granular phase (FA) in self-compacting mortar: (a) ϕ m FA ; (b) [ η ] FA .
Materials 14 01971 g008
Figure 9. Graph of the Bayesian network of the Krieger–Dougherty equation in self-compacting concretes.
Figure 9. Graph of the Bayesian network of the Krieger–Dougherty equation in self-compacting concretes.
Materials 14 01971 g009
Figure 10. Probability density functions of the parameters of the self-compacting concretes (SCC) phases in Feys et al. [39]: (a) ϕ m FA ; (b) [ η ] FA ; (c) ϕ m CA ; (d) [ η ] CA ; (e) [ η ] FA ϕ m FA ; and (f) [ η ] CA ϕ m CA .
Figure 10. Probability density functions of the parameters of the self-compacting concretes (SCC) phases in Feys et al. [39]: (a) ϕ m FA ; (b) [ η ] FA ; (c) ϕ m CA ; (d) [ η ] CA ; (e) [ η ] FA ϕ m FA ; and (f) [ η ] CA ϕ m CA .
Materials 14 01971 g010
Figure 11. Probability density functions of the parameters of the SCC phases in Esmaeilkhanian et al. [40]: (a) ϕ m FA ; (b) [ η ] FA ; (c) ϕ m CA ; (d) [ η ] CA ; (e) [ η ] FA ϕ m FA ; and (f) [ η ] CA ϕ m CA .
Figure 11. Probability density functions of the parameters of the SCC phases in Esmaeilkhanian et al. [40]: (a) ϕ m FA ; (b) [ η ] FA ; (c) ϕ m CA ; (d) [ η ] CA ; (e) [ η ] FA ϕ m FA ; and (f) [ η ] CA ϕ m CA .
Materials 14 01971 g011
Figure 12. Bivariate histograms in every SCC phase: (a) [ η ] FA and ϕ m FA in Feys et al. [39]; (b) [ η ] CA and ϕ m CA in Feys et al. [39]; (c) [ η ] FA and ϕ m FA in Esmaeilkhanian et al. [40]; and (d) [ η ] CA and ϕ m CA in Esmaeilkhanian et al. [40].
Figure 12. Bivariate histograms in every SCC phase: (a) [ η ] FA and ϕ m FA in Feys et al. [39]; (b) [ η ] CA and ϕ m CA in Feys et al. [39]; (c) [ η ] FA and ϕ m FA in Esmaeilkhanian et al. [40]; and (d) [ η ] CA and ϕ m CA in Esmaeilkhanian et al. [40].
Materials 14 01971 g012
Table 1. Values of ϕ m p and [ η ] p calculated in cement pastes [5]. Average range in rotational rheometry: 0–600 s 1 . CI stands for confidence interval.
Table 1. Values of ϕ m p and [ η ] p calculated in cement pastes [5]. Average range in rotational rheometry: 0–600 s 1 . CI stands for confidence interval.
Cementitious Material SP / cm η w [mPa s] ϕ m p 95% CI [ η ] p 95% CI R adj 2
CEM I 52.5 N-SR0.40.9330.680.43–0.936.64.4–7.70.980
0.80.9330.700.41–0.986.43.9–7.40.977
0.4, 0.80.9330.680.41–0.966.54.0–7.60.922
1.00.8910.750.56–0.946.85.5–7.70.981
1.20.8910.690.62–0.766.55.9–6.90.995
1.0, 1.20.8910.720.63–0.806.76.1–7.10.986
0.4, 0.8, 1.0, 1.20.9110.740.66–0.836.86.3–7.20.973
CEM II 32.5 BL-II0.40.9330.570.53–0.615.04.6–5.30.999
0.80.9330.830.54–1.006.24.8–6.80.993
0.4, 0.80.9330.660.52–0.795.54.5–6.20.979
75% CEM I 52.5 N-SR
+ 25% GGBS
0.40.9330.600.36–0.846.02.8–7.30.954
0.80.9330.830.50–1.006.54.8–7.30.991
0.4, 0.80.9330.630.41–0.856.13.9–7.20.917
75% CEM II 32.5 BL-II
+ 25% GGBS
0.40.9330.690.26–1.006.00.2–7.30.950
0.80.9330.510.45–0.564.33.6–4.90.987
0.4, 0.80.9330.580.42–0.735.13.4–6.10.909
All0.4, 0.8, 1.0, 1.20.9110.610.55–0.675.75.2–6.20.837
Table 2. Values of the parameters ϕ m p and [ η ] p in cement pastes calculated by Struble et al. [1]. Range of measure in rotational rheometer: 0–600 s 1 .
Table 2. Values of the parameters ϕ m p and [ η ] p in cement pastes calculated by Struble et al. [1]. Range of measure in rotational rheometer: 0–600 s 1 .
Cement Paste γ ˙ [ s 1 ] ϕ m p [ η ] p
Cement type I (dispersed)250.645.1
Cement type I (dispersed)5000.766.2
Cement type I (floculated)5000.646.3
White cement (dispersed)250.675.7
White cement (dispersed)5000.806.8
Cement type V (dispersed)Low limit0.704.7
Cement type V (dispersed)250.684.5
Cement type V (dispersed)5000.755.2
Table 3. Definition domains of the parameters ϕ m p and [ η ] p .
Table 3. Definition domains of the parameters ϕ m p and [ η ] p .
Investigation ϕ m p [ η ] p
Struble et al. [1]0.64–0.804.5–6.8
Present work & De La Rosa et al. [5]0.51–0.834.3–6.8
Table 4. Statistics of the parameters ϕ m p and [ η ] p .
Table 4. Statistics of the parameters ϕ m p and [ η ] p .
ParameterMeanStd. Dev.Percentage 2.5%MedianPercentage 97.5%
ϕ m p 0.8700.0310.8000.8790.899
[ η ] p 6.6510.1756.3336.6816.829
Table 5. Statistics of the parameters ϕ m FA and [ η ] FA for self-compacting mortar.
Table 5. Statistics of the parameters ϕ m FA and [ η ] FA for self-compacting mortar.
InvestigationParameterMeanStd. Dev.Percentage 2.5%MedianPercentage 97.5%
Ouro et al. [38] ϕ m FA 0.6550.0430.5650.6620.715
[ η ] FA 2.7530.1962.5102.7143.221
Table 6. Comparison between experimental values and model results for self-compacting mortar (Ouro et al. [38]).
Table 6. Comparison between experimental values and model results for self-compacting mortar (Ouro et al. [38]).
NameExperimental η SCM [Pa s]Calculated with Theoretical Values [7] η SCM [Pa s]Calculated with Bayesian Mean Values, η SCM [Pa s]Error with Theoretical Values [7] [%]Error with Bayesian Mean Values [%]
N3, N41.901.351.1729.038.6
N7, N81.101.120.982.211.3
N111.661.481.2610.524.1
N131.411.100.9721.931.3
N150.941.281.0836.015.3
E3, E40.501.060.87112.473.3
E70.350.490.4041.215.2
E130.370.720.5993.359.6
Table 7. Statistics of the parameters ϕ m FA , ϕ m CA , [ η ] FA , and [ η ] CA for self-compacting concretes.
Table 7. Statistics of the parameters ϕ m FA , ϕ m CA , [ η ] FA , and [ η ] CA for self-compacting concretes.
InvestigationParameterMeanStd. Dev.Percentage 2.5%MedianPercentage 97.5%
Feys et al. [39] ϕ m FA 0.6400.0470.5570.6430.713
[ η ] FA 3.6820.7432.5653.5975.206
ϕ m CA 0.7280.0990.5600.7310.886
[ η ] CA 5.6111.6072.7135.6728.402
Esmaeilkhanian et al. [40] ϕ m FA 0.6750.0310.6050.6820.716
[ η ] FA 3.1200.4252.5253.0524.009
ϕ m CA 0.7300.0990.5600.7340.887
[ η ] CA 4.7761.2912.6284.7517.151
Grünewald [41] ϕ m FA 0.6820.0290.6100.6890.716
[ η ] FA 4.0320.4252.9544.1024.679
ϕ m CA 0.7400.0970.5620.7490.887
[ η ] CA 5.2941.0943.1835.2417.741
Table 8. Comparison between experimental values and models, and error estimated for self-compacting concrete (Feys et al. [39]).
Table 8. Comparison between experimental values and models, and error estimated for self-compacting concrete (Feys et al. [39]).
NameExperimental η SCC [Pa s]Calculated with Theoretical Values [7] η SCC [Pa s]Calculated with Bayesian Mean Values, η SCC [Pa s]Error with Theoretical Values [7] [%]Error with Bayesian Mean Values [%]
SCC150.811.050.778.30.1
SCC242.611.149.874.117.0
SCC338.010.949.471.230.0
SCC441.411.452.972.527.8
SCC767.512.156.882.115.9
SCC859.011.553.580.49.3
SCC928.010.347.163.268.4
SCC1045.010.550.176.611.3
SCC1135.08.134.576.71.4
SCC1296.515.278.284.219.0
SCC1341.511.150.373.421.3
SCC1429.311.553.960.784.1
SCC1549.610.648.178.63.0
SCC1655.610.551.981.26.6
SCC1744.89.852.278.116.6
SCC1871.211.352.184.126.7
SCC1971.212.358.192.162.7
Table 9. Comparison between experimental values and models, and error estimated for self-compacting concrete (Esmaeilkhanian et al. [40]).
Table 9. Comparison between experimental values and models, and error estimated for self-compacting concrete (Esmaeilkhanian et al. [40]).
NameExperimental η SCC [Pa s]Calculated with Theoretical Values [7] η SCC [Pa s]Calculated with Bayesian Mean Values, η SCC [Pa s]Error with Theoretical Values [7] [%]Error with Bayesian Mean Values [%]
SCC159.319.534.467.242.0
SCC229.017.130.141.23.9
SCC369.519.634.671.850.2
SCC4124.016.128.987.076.7
SCC525.017.029.932.119.7
SCC762.031.651.849.016.5
SCC825.011.921.552.213.8
SCC928.032.752.616.887.9
SCC1072.038.061.547.214.6
SCC11128.041.767.867.447.1
SCC1271.037.660.847.114.3
SCC1367.037.460.744.19.4
SCC1435.031.851.59.047.1
SCC1637.031.651.214.538.3
SCC1739.031.651.218.931.2
SCC1830.031.751.55.771.6
Table 10. Comparison between experimental values, models, and error estimated for self-compacting concrete (Grünewald [41]).
Table 10. Comparison between experimental values, models, and error estimated for self-compacting concrete (Grünewald [41]).
NameExperimental η SCC [Pa s]Calculated with Theoretical Values [7] η SCC [Pa s]Calculated with Bayesian Mean Values, η SCC [Pa s]Error with Theoretical Values [7] [%]Error with Bayesian Mean Values [%]
OS169.216.174.676.87.9
OS259.411.449.080.917.4
OS387.926.091.870.54.5
OS456.015.854.571.82.7
OS597.649.3167.849.571.9
OS681.026.392.967.614.6
OS762.216.255.974.010.1
OS871.321.170.670.41.0
OS957.514.146.075.520.0
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

De La Rosa, Á.; Ruiz, G.; Castillo, E.; Moreno, R. Calculation of Dynamic Viscosity in Concentrated Cementitious Suspensions: Probabilistic Approximation and Bayesian Analysis. Materials 2021, 14, 1971. https://doi.org/10.3390/ma14081971

AMA Style

De La Rosa Á, Ruiz G, Castillo E, Moreno R. Calculation of Dynamic Viscosity in Concentrated Cementitious Suspensions: Probabilistic Approximation and Bayesian Analysis. Materials. 2021; 14(8):1971. https://doi.org/10.3390/ma14081971

Chicago/Turabian Style

De La Rosa, Ángel, Gonzalo Ruiz, Enrique Castillo, and Rodrigo Moreno. 2021. "Calculation of Dynamic Viscosity in Concentrated Cementitious Suspensions: Probabilistic Approximation and Bayesian Analysis" Materials 14, no. 8: 1971. https://doi.org/10.3390/ma14081971

APA Style

De La Rosa, Á., Ruiz, G., Castillo, E., & Moreno, R. (2021). Calculation of Dynamic Viscosity in Concentrated Cementitious Suspensions: Probabilistic Approximation and Bayesian Analysis. Materials, 14(8), 1971. https://doi.org/10.3390/ma14081971

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