Prediction of the Rheological Properties of Fresh Cementitious Suspensions Considering Microstructural Parameters

Supplementary cementitious materials (SCMs) are commonly used to partially replace cements. Although it is necessary to investigate the rheological properties of the individual supplementary cementitious materials (SCMs) for understanding complex rheological behaviors of the blended mixes, the study on the investigation of rheological properties of various SCMs such as fly ash, blast-furnace slag, and silica fume, according to various solid volume fractions and prediction models is fairly limited. This study investigated the rheological properties of non-blended cementitious suspensions with Portland cement (PC), fly ash (FA), blast-furnace slag (BS), and silica fume (SF) materials in the experiments and predicted using YODEL (Yield stress mODEL) and Krieger–Dougherty’s (K–D’s) equation. Experiments were designed with various solid volume fractions (ϕ) from 0.28 to 0.44, and the rheological properties of all cementitious suspensions were noted to increase with increasing ϕ, showing an improved flowability at low ϕ. YODEL, derived from the first principles considering particle-size distributions, interparticle forces and microstructural parameters predicted the yield stress. The YODEL predictions were consistent with the experiments with a positive correlation coefficient of above 0.96. K–D’s equation with the maximum particle fractions and intrinsic viscosity as key parameters predicted the plastic viscosity. The K–D’s equation predictions match up with the experiments with a positive correlation coefficient of above 0.94. Both models showed more quantitative predictions without any fitting parameters and could be applied to any multimodal powder suspensions.


Introduction
The rheological properties of fresh cementitious suspensions are of great interest for many applications in the construction industry [1][2][3][4]. Advances in the construction industry have led to the design of new and special mixtures such as self-compacting, ultra-highperformance, and engineered cementitious composites that can significantly improve the mechanical properties and impermeability, and reduce intrusion of aggressive agents [5][6][7][8][9][10]. These improved characteristics could be hindered by the inadequate rheological properties, which facilitates the placement process such as pumping, casting, and molding [11,12]. Sustainable cementitious mixtures with recycled aggregates or construction waste fines also have difficulty in obtaining proper rheological properties [13,14]. In general, the cementitious mixtures incorporate multi-scale fine particles, admixtures, fibers, superplasticizers, etc., making it difficult to experimentally investigate their rheological properties [15,16]. In addition, there are several factors that influence the rheological properties such as material properties, interparticle forces due to dispersion forces and electrostatic interactions, steric forces from adsorbed polymers, hydrodynamic reactions, and crowding factors [1][2][3][4]. Such complex behaviors make the experimental work more complicated and are in need of adopting multi-scale modeling approaches to predict and understand the rheological properties with minimal experimentations. widely applied to general cementitious paste, mortar and concrete, but the applications were undertaken only in macroscopic approach without consideration for the influence of individual material. This shortage limits the extendibility of application of the models to various kinds of cementitious mixtures. Therefore, in this study, it was intended to determine the parameters of the models suitable for suspensions with different kinds of cementitious materials, widening the application.
Moreover, the development of special cementitious mixtures over the last three decades requires a high binder volume, mainly achieved by partially replacing Portland cement (PC) with supplementary cementitious materials (SCMs) [31,32]. SCMs such as fly ash (FA), blast-furnace slag (BS), and silica fume (SF) are commonly used due to their ability to partially replace cement, making a more effective binder [33,34]. In addition to this, SCMs are also used for the purposes of improving strength, particle packing efficiency, and durability, sometimes for the purpose of reducing permeability, alkali-silica reaction, heat evolution during hydration [35][36][37]. To achieve the above intended positive effects, a mixture designed with SCMs should carefully consider the workability as well as the rheological properties. The partial replacement of SCMs can result in complex rheological behavior influenced by the physical properties as well as by the type and replacement ratio. The physical properties including specific surface area (SSA), packing density, particle shape and particle size influence the rheological behavior either by increasing or decreasing the flowability characteristics. Studies were focused on the influence of SCMs to rheology of cement pastes that have been designed as blended mixtures such as binary, tertiary, and quaternary mixtures [38][39][40][41]. However, these complex rheological behaviors are difficult to understand with the blended mixes due to uncontrolled fluctuations in the particle's properties. This can be achieved by investigating the rheological properties of the individual SCMs as non-blended mixes. Recently, a study designed non-blended cementitious suspensions with individual SCMs that considered the influence of inter-particle distances on the rheological properties [42]. Therefore, more detailed study to investigate the rheological behavior of non-blended mixes could help to better understand each individual rheological behavior and the factors influencing the rheological properties more precisely. For the purpose, this study investigated the rheological properties of non-blended cementitious suspensions with PC, FA, BS, and SF materials in the experiments and predicted using YODEL (Yield stress mODEL) and Krieger-Dougherty's (K-D's) equation. Based on the rheological information of non-blended mixes, the rheological properties of any blended mixtures with SCMs can be predicted and controlled to obtain more effective mixtures with adequate rheological properties.

YODEL
Flatt and Bowen [1,2] proposed YODEL, a yield stress model for multimodal suspensions based on the microstructural considerations of colloidal particle interactions. This equation was derived in consideration of the interparticle forces, suspension microstructure, and particle-size distribution. This also includes the physical parameters such as particle size, geometrical maximum packing fraction, percolation threshold, and minimum separation distance at contact. Based on the YODEL, the yield stress function (τ 0 ) was expressed as in Equation (1), where, m 1 indicates the predetermined factor that considers interparticle forces, particle size, and particle-size distribution as expressed in Equation (2). φ represents the solid volume fraction of the cementitious suspensions. φ 0 represents the percolation volume fraction that depends on the interaction between Brownian motion (dispersive) and colloidal attractive forces between particles. φ m represents the maximum packing fraction of cementitious particles.
where, R v,50 was the median volume radius, was derived from the particle-size distribution, and a * indicates the average (characteristic) radius of curvature at particle contacts, which was introduced to describe the dependence of m 1 on interparticle force. G max was the maximum attractive interparticle force normalized by the radius of curvature at the contact points as expressed in Equation (3), where, A 0 and H represents the Hamaker's constant and minimum separation distance. In addition, the effect of particle-size distribution was considered by the function f * σ,∆ and expressed in Equation (4), where, b i is the normalized particle radii by normalizing the particle radii a i by the mean volume radius R v, 50 . φ k indicates the volume fraction of particles of size b k in the size interval k. A s A c and S a,1 were derived from the geometrical model of Suzuki et al. [43], as shown in Equations (5) and (6). ∆v k,1 was a geometrical term that accounts for a change in the maximum packing fraction induced by each pair of undispersed particles of sizes a k and a l , which is expressed as shown in Equation (7). u k,k indicates the normalized factor as expressed in Equation (8), with the assumption of enclosing sphere model for the effective volume fraction of solids. In Flatt and Bowen [1,2], three geometrical models (truncated cone model without the particle portion, truncated cone model with the particle volume, and the enclosing sphere model without the particle volume) for effective volume fraction of solids were introduced to account for an increase in the effective volume of solids from undispersed particles, which is used in calculating ∆v k,1 and u k,k . Given that explaining the choice of a geometrical model is somewhat arbitrary, they suggested that the enclosing sphere model provided a better result in the data analysis, even though the impact of this choice has a less pronounced effect on size distribution. Therefore the equations of ∆v k,1 and u k,k based on the enclosing sphere model were adopted for this study.

Krieger-Dougherty's (K-D) Equation
The K-D's equation [27] predicts the plastic viscosity from the volume fraction of cementitious particles as expressed in Equation (9), where, µ represents the plastic viscosity of the suspensions; µ 0 represents the plastic viscosity of the fluid phase of the suspensions; φ indicates the solid volume fraction of the cementitious suspensions; φ m indicates the maximum packing fraction of cementitious particles, and; [η] represents the intrinsic viscosity which mainly depends on the shape of particles.

Materials
The raw materials used in this study consisted of: PC, FA, BS, and SF materials. The physical properties and the mean particle sizes (D mean ) of the PC, FA, BS, and SF are shown in Table 1. A scanning electron microscope (SEM, S-4300, Hitachi, Tokyo, Japan) was used to observe the micro-morphology of cementitious materials as shown in Figure 1. In the SEM images, the PC and BS consist of mostly angular particles, FA consists of spherical particles, and SF consists of an agglomerated spheroidal particle. The chemical compositions of PC, FA, BS, and SF were obtained using an X-ray fluorescence spectrometer (PW2400, Philips, Amsterdam, The Netherland), of which the main oxides are presented in Table 2. The cementitious materials predominantly consist of major oxides that includes SiO 2 , Al 2 O 3 , CaO, and Fe 2 O 3 , and a small amount of minor oxide such as MgO, respectively. A laser scattering particle-size distribution analyzer (LS 13 320, Beckman Coulter, Brea, CA, USA) was used to measure the particle-size distributions (PSDs) of PC, FA, BS, and SF particles as shown in Figure 2.

Mixture Formulations and Procedures
Four groups of non-blended mixtures were prepared using PC, FA, BS, and SF materials with different water-to-solid volume ratios. The water-to-solid volume ratios were quantified as the volumetric ratios, as the rheological properties generally governed by volume ratio rather than weight ratio. All the mixtures were proportioned with different water-to-solid volume ratios (w/s)v varied from 1.25, 1.50, 1.75, 2.00, 2.25, to 2.50, which

Mixture Formulations and Procedures
Four groups of non-blended mixtures were prepared using PC, FA, BS, and SF materials with different water-to-solid volume ratios. The water-to-solid volume ratios were quantified as the volumetric ratios, as the rheological properties generally governed by volume ratio rather than weight ratio. All the mixtures were proportioned with different water-to-solid volume ratios (w/s) v varied from 1.25, 1.50, 1.75, 2.00, 2.25, to 2.50, which correspond to the mass-based water-to-solid volume ratios (w/s) m between 0.40 to 1.00. The mix ratios that includes (w/s) v , and the corresponding (w/s) m ratios, and φ are shown in Table 3. For an individual cementitious material, six non-blended mixtures were prepared with different (w/s) v ratios, which totally consist of twenty-four mixtures. All the materials were preconditioned to a constant temperature (20 ± 3 • C) for 24 h inside the laboratory to minimize the temperature difference. To produce the mixtures, the solid ingredients were mixed in a dry state for 2 min in a Hobart-mixer. The water was slowly added to the dry mixture and mixed for about 3 min. Then, the sides of the mixer container were scraped, and the mixing continued until a homogeneous mixture was obtained.

Rheological Measurements
Immediately after the mixing, the mixtures were subjected to rheological tests using a Brookfield DV2T rheometer(AMETEC Brookfield, Middleborough, MA, USA) equipped with mortar type spindle. The rheometer measured the shear stress response to the applied strain rate, and the data was recorded using a PC connected to the rheometer. Figure 3 shows the spindle type, experimental set-up, and applied shear rates. This study adopted the Bingham model, which includes yield stress and viscosity as shown in Equation (10), where, τ 0 and µ indicate the yield stress and plastic viscosity of the cementitious suspensions. The torque required to rotate the spindle in the mixtures was determined, and the shear rate was increased from 0 to 33.15 s −1 and decreased from 33.15 to 0 s −1 . The shear rate at each step was maintained for 10 s to detect a stable shear stress. The curves at the decreasing rate were more consistent and a linear regression was performed to determine the plastic viscosity and yield stress with the slope and intercept of the regression analysis line plotted through shear stress against shear rate. For each mixture sample, three repetitive tests were performed, and the average yield stress and plastic viscosity were obtained.

Yield Stress in the Experiments with Different
From the rheological tests, the yield stresses of cementitious suspensions were obtained using a linear regression plotted through shear stress on the shear rate. Three tests were performed for each mixture and the average values with deviations were also obtained. Figure 4 shows the obtained yield stress of cementitious suspensions with different . For all the cementitious suspensions, the yield stresses were noted to increase with increasing , showing an improved flowability with the low . A non-linear relationship between and yield stress was established using the relationship = , as shown in the figure. By using the relationship, the maximum yield stress ( ), which can be theoretically obtained at = 1 and the increase factor ( ) with increasing ( ) were estimated. The maximum yield stresses of the PC, FA, BS, and SF suspensions were 3147.1 Pa, 113.98 Pa, 247.49 Pa, and 5463.9 Pa, respectively. The increase factors with increasing of the PC, FA, BS, and SF suspensions were 5.90, 4.50, 4.16, and 3.79, respectively. The above relationship can be used to predict the yield stress of the cementitious suspensions with different . Moreover, when comparing the yield stress of the PC with the FA and BS suspensions at the same , the FA and BS suspensions were observed to provide low yield-stress values, while the SF suspensions showed the opposite tendency. The yield stress of the FA suspensions reduced by 75% to 90% compared to PC suspensions, indicating that the FA suspensions were more fluidic mixtures. These fluidic observations can be governed by the spherical shape of FA particles and the filling effects that increase the physical separation between particles, which reduces the yield stress [44,45]. Similar to From the rheological tests, the yield stresses of cementitious suspensions were obtained using a linear regression plotted through shear stress on the shear rate. Three tests were performed for each mixture and the average values with deviations were also obtained. Figure 4 shows the obtained yield stress of cementitious suspensions with different φ. For all the cementitious suspensions, the yield stresses were noted to increase with increasing φ, showing an improved flowability with the low φ. A non-linear relationship between φ and yield stress was established using the relationship y = ax b , as shown in the figure. By using the relationship, the maximum yield stress (a), which can be theoretically obtained at φ = 1 and the increase factor (b) with increasing φ (x) were estimated. The maximum yield stresses of the PC, FA, BS, and SF suspensions were 3147.1 Pa, 113.98 Pa, 247.49 Pa, and 5463.9 Pa, respectively. The increase factors with increasing φ of the PC, FA, BS, and SF suspensions were 5.90, 4.50, 4.16, and 3.79, respectively. The above relationship can be used to predict the yield stress of the cementitious suspensions with different φ. Moreover, when comparing the yield stress of the PC with the FA and BS suspensions at the same φ, the FA and BS suspensions were observed to provide low yield-stress values, while the SF suspensions showed the opposite tendency. The yield stress of the FA suspensions reduced by 75% to 90% compared to PC suspensions, indicating that the FA suspensions were more fluidic mixtures. These fluidic observations can be governed by the spherical shape of FA particles and the filling effects that increase the physical separation between particles, which reduces the yield stress [44,45]. Similar to the FA suspensions, the BS suspensions have lower yield-stress values than the PC suspensions, i.e., reduced by 28% to 70%. This behavior may be associated with the smooth surface of BS particles, lowchemical effect, and the micro-filling effect [38,39]. In contrast, the yield stress of the SF suspensions increased significantly by 782% to 2372% compared to the PC suspensions. This can be explained by the dense packing and the fine size of SF particles, which results in the reduction in particle spacing and an increase in the number of direct contact points between particles, thus making it difficult for the particles to slide on each other [add]. The increase in thixotropy with SF is also unfavorable to the initial flow [35]. In addition, SF suspension shows relatively higher variation from the regression curve. It seems to be due to the difficulty in obtaining good dispersion of the fine particles with consistency because of the higher possibility of agglomeration for finer particles. The same tendency can be found in the results of plastic viscosity.

Calculation of YODEL Parameters
To obtain quantitative predictions using YODEL, the first step was to measure the PSDs of the cementitious materials. The PSDs of the cementitious materials were measured using the laser scattering particle-size distribution analyzer and are displayed in Figure 2. From the PSDs, the mean particle sizes of the PC, FA, BS, and SF were 21.58 µm, 30.20 µm, 15.85 µm, and 6.31 µm, respectively. In general, the cementitious materials exhibit a quasi-log-normal size distribution, and the cumulative distribution must fit well with the log-normal distributions, so that the ,∆ * values can be estimated using directly measured PSDs [1,2]. The log-normal size distributions of the cementitious materials were obtained and shown in Figure 5. The obtained log-normal distributions were fitted well to the cumulative distributions for all the cementitious materials. Therefore, the ,∆ * values were estimated as a function of the normal PSDs. The mean and median particle sizes of the cementitious materials, and the estimated ,∆ * and , (calculated using enclosing sphere model) values are shown in Table 4. The values were then calculated using 0 and . Firstly, 0 were estimated by using Equation (11) for the cementitious particles of same kind, where, 1 and 3 indicates the Hamaker's constant of cementitious particles and the suspending fluid (Water = 4.38 × 10 −20 J). The 1 was estimated by using the rule of mixtures as shown in Equation (12),

Calculation of YODEL Parameters
To obtain quantitative predictions using YODEL, the first step was to measure the PSDs of the cementitious materials. The PSDs of the cementitious materials were measured using the laser scattering particle-size distribution analyzer and are displayed in Figure 2. From the PSDs, the mean particle sizes of the PC, FA, BS, and SF were 21.58 µm, 30.20 µm, 15.85 µm, and 6.31 µm, respectively. In general, the cementitious materials exhibit a quasi-log-normal size distribution, and the cumulative distribution must fit well with the log-normal distributions, so that the f * σ,∆ values can be estimated using directly measured PSDs [1,2]. The log-normal size distributions of the cementitious materials were obtained and shown in Figure 5. The obtained log-normal distributions were fitted well to the cumulative distributions for all the cementitious materials. Therefore, the f * σ values were estimated as a function of the normal PSDs. The mean and median particle sizes of the cementitious materials, and the estimated f * σ,∆ and u k,k (calculated using enclosing sphere model) values are shown in Table 4. The G max values were then calculated using A 0 and H. Firstly, A 0 were estimated by using Equation (11) for the cementitious particles of same kind, where, A 1 and A 3 indicates the Hamaker's constant of cementitious particles and the suspending fluid (Water = 4.38 × 10 −20 J). The A 1 was estimated by using the rule of mixtures as shown in Equation (12), where, m 1 , m 2 , and m n indicates the mass (%) of each chemical component of cementitious materials; m indicates the total mass (%) of cementitious materials; Ah 1 , Ah 2 , and Ah n represents the corresponding Hamaker's constant of each chemical component. The chemical components of the cementitious materials measured through X-ray fluorescence spectroscopy and their corresponding Hamaker's constants [46,47] are given in Table 5.
As shown in the table, the estimated A 0 of the PC, FA, BS, and SF suspensions were 1.628 × 10 −20 J, 1.449 × 10 −20 J, 1.638 × 10 −20 J, and 6.885 × 10 −21 J, respectively. In the cementitious suspensions, the H lies in the order of 2 nm for flocculated systems and 10 nm for fully dispersed systems [4]. This study investigated the rheological properties of cementitious suspensions without any chemical additives such as superplasticizers, and therefore the suspensions were not totally dispersed systems. The H values were adjusted to fit the experimental yield-stress values and the obtained H of the PC, FA, BS, and SF suspensions were 2.5 nm, 7.0 nm, 6.1 nm, and 0.5 nm, respectively. Then, m 1 depending on the interparticle forces was calculated, assuming the contact points between particles as the average fixed radius of curvature a * [1,2]. The interparticle forces would not be proportional to the particle size of a material, but largely independent of the particle size and a * is therefore determined for each material regardless of particle size. The a * of alumina particles was estimated as 20 nm in Flatt and Bowen [1,2], and Roussel et al. [4] suggested that the a * of cement particles was 20 times larger than alumina particles. Considering the lack of informative previous studies about a * for all the materials and the limited experiments for this study, a constant a * was assumed for all the cementitious materials used for simplicity of analysis. Therefore, the a * was fixed as 300 nm for all the cementitious particles. The φ m was calculated using de Larrard's compressive packing model (CPM) [48], and the calculated φ m of the PC, FA, BS, and SF particles were 0.59, 0.60, 0.59, and 0.69, respectively. In addition, the φ 0 values considered for the cementitious materials were ranged from 0.024 to 0.04 and was assumed as 0.03 for all the suspensions [4]. The a * , φ 0 , and φ m are given in Table 6.  used for simplicity of analysis. Therefore, the * was fixed as 300 nm for all the cementitious particles. The was calculated using de Larrard's compressive packing model (CPM) [48], and the calculated of the PC, FA, BS, and SF particles were 0.59, 0.60, 0.59, and 0.69, respectively. In addition, the 0 values considered for the cementitious materials were ranged from 0.024 to 0.04 and was assumed as 0.03 for all the suspensions [4]. The * , 0 , and are given in Table 6.

Prediction of Yield Stress Using YODEL
By using the above estimated microstructural parameters required by YODEL, the yield stresses of the cementitious suspensions were predicted and fitted with the experimental data. Figure 6 shows the yield stress of cementitious suspensions obtained through the experiments and predicted through YODEL. Similar to the experimental data, the yield stresses of all the cementitious suspensions increased with the increase in φ in the YODEL predictions. The correlation coefficient (r), and the regression analysis with 95% confidence intervals (i.e., maximum, and minimum) were estimated between the experiments and YODEL predictions, as shown in the figure. In the PC suspensions, the YODEL predicted the yield stress with a mean standard deviation to the experimental data of approximately 1.14. In addition, the r between the experiment and prediction was 0.98, respectively. The average standard deviation and r values between the experiment and prediction of FA suspensions were 0.11 and 0.99, respectively. For the BS suspensions, the average standard deviation and r values between the experiment and prediction were 0.58 and 0.96, respectively. In addition, the average standard deviation and r values between the experiment and prediction of SF suspensions were 16 and 0.96, respectively. From the above results, the YODEL predicted the yield stress of cementitious suspensions with positive r above 0.96. As YODEL adopts microstructural approaches to predict yield stress as a function of PSDs, interparticle forces, and packing fractions, a more quantitative predictive capacity becomes apparent and was certainly confirmed in the above results. The PSDs plays a very prominent role, as it is the main input parameter. The YODEL also designed without introducing any fitting parameters and introduced CPM model to estimate φ m . Therefore, an appropriate input of PSDs and φ m was very essential for YODEL to obtain quantitative predictions.
Materials 2022, 15, x FOR PEER REVIEW 12 of 18 a function of PSDs, interparticle forces, and packing fractions, a more quantitative predictive capacity becomes apparent and was certainly confirmed in the above results. The PSDs plays a very prominent role, as it is the main input parameter. The YODEL also designed without introducing any fitting parameters and introduced CPM model to estimate . Therefore, an appropriate input of PSDs and was very essential for YODEL to obtain quantitative predictions.  Figure 7 shows the obtained plastic viscosity of cementitious suspensions with different . For all the cementitious suspensions, the plastic viscosity were noted to increase with increasing , showing an improved flowability with the low . A non-linear relationship between and plastic viscosity was established using the exponential function ( = ). By using the relationship, the initial plastic viscosity ( ) when the was 0, and the increase factor ( ) with increasing was estimated. By using the relationship, the plastic viscosity of cementitious suspensions with different can be predicted. The initial plastic viscosity was the measure of water, which means that there was no solid particles. The initial plastic viscosity of the PC, FA, BS, and SF suspensions were predicted as 0.0001 Pa.s, 0.00002 Pa.s, 0.0001 Pa.s, and 0.0196 Pa.s, which were relatively equivalent  Figure 7 shows the obtained plastic viscosity of cementitious suspensions with different φ. For all the cementitious suspensions, the plastic viscosity were noted to increase with increasing φ, showing an improved flowability with the low φ. A non-linear relationship between φ and plastic viscosity was established using the exponential function (y = ae −bx ). By using the relationship, the initial plastic viscosity (a) when the φ was 0, and the increase factor (b) with increasing φ was estimated. By using the relationship, the plastic viscosity of cementitious suspensions with different φ can be predicted. The initial plastic viscosity was the measure of water, which means that there was no solid particles. The initial plastic viscosity of the PC, FA, BS, and SF suspensions were predicted as 0.0001 Pa.s, 0.00002 Pa.s, 0.0001 Pa.s, and 0.0196 Pa.s, which were relatively equivalent to those of water (0.01 Pa.s at 20 • C, decreases as the temperature increases) [49] respectively. The low initial plastic viscosity of cementitious suspensions predicted from the experimental results may be due to the chemical reactions and initial hydration process that immediately raises the temperature when mixed with water. The increase factor with increasing φ of the PC, FA, BS, and SF suspensions were 21.98, 20.62, 20.38, and 11.56, respectively. Similar to the yield stress, the FA and BS suspensions were observed to provide lower plastic viscosity values than PC suspensions, whereas an opposite behavior was observed in SF suspensions. The plastic viscosity of the FA suspensions reduced by 78% to 86% compared to PC suspensions, indicating that the FA suspensions were more fluidic mixtures. These fluidic observations were explained above due to shape and filling effects [44,45]. The BS suspensions also showed a lower plastic viscosity value than the PC suspensions, reduced by 8% to 44%. It was due to the smooth surface of BS particles, low-chemical and micro-filling effects [38,39]. On contrary, the plastic viscosity of the SF suspensions increased significantly by 95% to 830% compared to the PC suspensions, caused by the fine size and dense packing. As explained previously, those properties make it difficult for the particles to slide over each other because of increased contacts between particles. SF suspension therefore presented much higher viscosity than PC, FA, and BS suspensions with a constant solid volume fraction. to those of water (0.01 Pa.s at 20 °C, decreases as the temperature increases) [49] respectively. The low initial plastic viscosity of cementitious suspensions predicted from the experimental results may be due to the chemical reactions and initial hydration process that immediately raises the temperature when mixed with water. The increase factor with increasing of the PC, FA, BS, and SF suspensions were 21.98, 20.62, 20.38, and 11.56, respectively. Similar to the yield stress, the FA and BS suspensions were observed to provide lower plastic viscosity values than PC suspensions, whereas an opposite behavior was observed in SF suspensions. The plastic viscosity of the FA suspensions reduced by 78% to 86% compared to PC suspensions, indicating that the FA suspensions were more fluidic mixtures. These fluidic observations were explained above due to shape and filling effects [44,45]. The BS suspensions also showed a lower plastic viscosity value than the PC suspensions, reduced by 8% to 44%. It was due to the smooth surface of BS particles, lowchemical and micro-filling effects [38,39]. On contrary, the plastic viscosity of the SF suspensions increased significantly by 95% to 830% compared to the PC suspensions, caused by the fine size and dense packing. As explained previously, those properties make it difficult for the particles to slide over each other because of increased contacts between particles. SF suspension therefore presented much higher viscosity than PC, FA, and BS suspensions with a constant solid volume fraction.

Calculation of K-D's Equation Parameters
The K-D equation mainly depends on two parameters: The for all the cementitious materials were calculated using de Larrard's CPM model [48]. The calculated of PC, FA, BS, and SF particles were 0.59, 0.60, 0.59, and 0.69, respectively. When the ranging between 0.6-0.7, the [ ] was assumed as 2.5 for monodisperse and polydisperse systems, 3 to 5 when the particles were sharp and angular, and between 4 to 10 when the particles consist of acicular, rods and fiber shapes [50]. The [ ] depends on the individual effect of particles and their shape [28,[50][51][52]. The [ ] can be chosen as 2.5 for rigid spherical geometries and must be modified if the particle shapes were found to be deviated [50][51][52]. An expression suggested by Pabst et al. [53] can be useful to estimate [ ], but the correlation between particle shape and [ ] was fundamentally complicated [50][51][52]. However, [ ] value generally assumes that all particles have a similar shape and most of the studies adopted [ ] to be adjusted to fit the experimental measurements [50][51][52]. Therefore, this study also adopted the [ ] to fit the K-D's equation by adjusting the

Calculation of K-D's Equation Parameters
The K-D equation mainly depends on two parameters: φ m and [η]. The φ m for all the cementitious materials were calculated using de Larrard's CPM model [48]. The calculated φ m of PC, FA, BS, and SF particles were 0.59, 0.60, 0.59, and 0.69, respectively. When the φ m ranging between 0.6-0.7, the [η] was assumed as 2.5 for monodisperse and polydisperse systems, 3 to 5 when the particles were sharp and angular, and between 4 to 10 when the particles consist of acicular, rods and fiber shapes [50]. The [η] depends on the individual effect of particles and their shape [28,[50][51][52]. The [η] can be chosen as 2.5 for rigid spherical geometries and must be modified if the particle shapes were found to be deviated [50][51][52]. An expression suggested by Pabst et al. [53] can be useful to estimate [η], but the correlation between particle shape and [η] was fundamentally complicated [50][51][52]. However, [η] value generally assumes that all particles have a similar shape and most of the studies adopted [η] to be adjusted to fit the experimental measurements [50][51][52]. Therefore, this study also adopted the [η] to fit the K-D's equation by adjusting the values to fit the experimental results. The adopted [η] values of PC, FA, BS, and SF particles were 9.0, 7.0, 8.5, and 13.0, respectively. Table 7 shows the estimated φ m and [η].

Conclusions
This study primarily investigated the rheological properties of cementitious suspensions, designed as non-blended mixtures with different φ. Based on the experimental and analytical results, the following conclusions can be drawn: (1) In the experiments, the rheological properties of all the cementitious suspensions were noted to increase with increasing φ, showing an improved flowability with the low φ.
Compared to the PC suspensions, the FA and BS suspensions showed an improved flowability, which was mainly due to the particle characteristics. The shape of FA was spherical, which facilitates isolation and dispersion, thus providing more flowable mixtures. In the case of BS, the improved flowability can be due to its smooth surface, less-chemical activity, micro-filling effect, and large surface area. However, the SF suspensions showed an opposite behavior, which was less flowable. As the SF consists of more fine particles, it could reduce particle spacing and increases direct contact points between particles, which makes it difficult for the particles to slide each other.