Thank you for visiting nature.com. You are using a browser version with limited support for CSS. To obtain the best experience, we recommend you use a more up to date browser (or turn off compatibility mode in Internet Explorer). In the meantime, to ensure continued support, we are displaying the site without styles and JavaScript.
Scientific Reports volume 13, Article number: 22301 (2023 ) Cite this article Mortar Polycarboxylate Superplasticizer
![Mathematical models of apparent viscosity as a function of water–cement/binder ratio and superplasticizer in cement pastes | Scientific Reports](http://cdn.globalso.com/water-reducer/011.jpg)
The water–cement/binder ratio and the admixture of water-reducing agents strongly affect the rheological properties of cement pastes. This study develops mathematical models to predict the apparent viscosity of cement pastes with varying water-cement/binder ratios and polycarboxylate-based superplasticizer content by introducing the power law shear stress-shear strain relation of non-Newtonian fluids into the Navier–Stokes motion equations. The developed models are compared with the results of rheological experiments and verified for their accuracy in simulating the apparent viscosity of cement pastes. These models provide insight into the rheological behaviour of cement pastes and could have practical applications in the construction industry.
The rheology of fresh cement paste is closely linked to its workability, primarily determined by the water-cement ratio (w/c)1,2,3,4. After adding mineral admixtures to cementitious materials, the water-binder ratio (w/b) is used instead5. Numerous rheological models have been developed to characterize the rheological properties of cement pastes. For example, Einstein6 proposed a rigid sphere suspension model with a linear relationship between viscosity and solid volume fraction, closely related to w/c or w/b7. Other models with non-linear relationships, such as the power law (Krieger-Dougherty model), exponential law (Mooney model), and fraction forms (Eilers model, Robinson model, Quemada model, et al.)8,9,10, have been applied to describe non-linear results in cement pastes.
Several rheological models have been used to analyse the evolution of shear stress with shear rate and w/c in cement pastes or similar materials4,11,12,13,14,15,16,17. These models include the Bingham model which presents a linear relationship between shear stress and shear rate11,13, and the modified Bingham model which applies quadratic functions to discuss the non-linear parts of the shear stress-shear rate in fresh pastes12. Other models, such as the Casson model which adds a square root function of the shear rate to the Bingham model13, and the Ostwald model which uses a power law to obtain the whole constitutive relation between shear stress and shear rate14, have also been employed. The Herschel–Bulkey model combines the Bingham model and the Ostwald model15,16, while the Bingham–Papanastasiou model has a Papanastasiou’s function multiplier based on the Bingham model4. The Vipulanandan model takes a fraction function as the shear rate17. The parameters in these models are often adjusted based on experimental results from rheological tests with various shear protocols14, and they often change as w/c varies4.
Apparent viscosity, an important rheological index, is influenced by several factors, including hydration degree, thixotropy, and the dosage of additives and admixtures14. The experimental rheological results of Liu et al.18, shown in Fig. 1, demonstrate the relations between the apparent viscosity and the shear rate of cement pastes. These relations change with the w/b and the dosage of superplasticizer (SP). The cement pastes can show shear thickening behavior18, for cement paste with lower w/b and a high dosage of SP, as shown in Fig. 1b and c.
The relations between the apparent viscosity and the shear rate for pastes with different w/b18. (a) paste with w/b = 0.24, (b) paste with w/b = 0.20, (c) paste with w/b = 0.16, (d) paste with w/b = 0.32.
Most of the present models are developed based on the test results in rheological experiments to consider the influence of w/c, additives, and admixtures. The representatives of the models are Jones and Taylor’s model (six-parameter exponential form)19, Lapasin et al. model (linear function of w/c)20, Ivanov and Roshavelov model (twenty-parameter polynomial with w/c, SP, silica fume, tricalcium aluminate and sulfur trioxide)21, and Rosquo\(\ddot{\text{e}}\) t model (linear function of w/c and power law of the shear rate)22. Other empirical formulations consider the yield stress in the Bingham model14, the Herschel–Bulkey model14, the YODEL model23,24, and the Ma & Kawashima model25 as a fractional or exponential function of the cement volume fraction.
Different from the experimental or empirical method, the target of this paper is to give a theoretical explanation of the change in the viscosity of cement pastes considering w/c(w/b) and SP based on mathematical derivations. Firstly, the Ostwald model is introduced into the Navier–Stokes equations of solid mixing with liquid. Next, an ordinary differential equation is derived by the mean value theorem in a tiny domain. To solve the ordinary differential equation, the Bernstein polynomial approximation is applied in the whole domain to derive a concise mathematical model. The model is a four-parameter formulation of the shear stress-shear rate for cement pastes with one variable w/c. More importantly, by introducing the electrostatic repulsion and steric hindrance into the Navier–Stokes equations, the influence of the SP can be quantitively considered by a function of its dosage. Correspondingly, the mathematical model can be extended to present the shear stress changes with the shear rate, w/c, and SP, with the same four-parameter formulation. The developed model is verified by being compared with the experimental results from Rosquo\(\ddot{\text{e}}\) t et al.22, Jeong et al.4, Cyr et al.26, and Liu et al.18.
According to Kundu et al.27, the Navier–Stokes equations are
in which \(\frac{D}{Dt}\) is the material derivative, \(\rho\) is the density of the fluid, \(P\) is the liquid pressure, \(\eta\) is the apparent viscosity of the fluid, \({u}_{i}\) is the velocity in the xi (i = 1, 2, 3) direction, \({f}_{i}\) is the body force on a unit volume of the fluid. In the Cartesian coordinate system, O-x1x2x3, i (i = 1, 2, 3) corresponds to the x, y, and z axis respectively.
A simply unidirectional flow with constant pressure is considered, therefore the Newton’s law of viscosity is applicable as
where \({u}_{1}\) depends on the space variable y and the time variable t, \({P}_{0}\) is the initial pressure. By substituting Eq. (2) into Eq. (1), Eq. (3) is obtained
Considering the procedure of mixing the solid and the liquid to produce fresh cement paste28, Eq. (3) is changed into
in which \({\phi }^{c}\) , the volume fraction of liquid, can be taken as water in cement paste
here \(\phi\) , the volume fraction of solid, can be taken as the binder particles for cement with additives such as silica fume and fly ash. The solid phase serves as the structural phase in the paste and is assumed to sustain the body force of the liquid phase. By introducing the drag force, we have
where \({u}_{s}\) and \({m}_{s}\) are the velocity and the mass of the solid respectively. In unit volume, \({m}_{s}\) can be expressed as
\({\rho }_{s}\) is the density of the solid. Substituting Eqs. (5)–(7) into Eq. (4), it is obtained
In a tiny time increment \(\Delta t\) , the variables \(\phi\) and \({\rho }_{s}\) own little changes. Ignoring the variation of \(\phi\) and \({\rho }_{s}\) , Eq. (8) can be written as
where w/b is given as
Then, Eq. (9) is changed into
The shear rate \(\dot{\gamma }\) is
Substituting it into Eq. (11), we have
In pure cement pastes, the value of w/b is
According to Ostwald14, the power law shear stress-shear rate relation is
where K and n are the calculation parameters. Then the apparent viscosity can be expressed as
In fact, the apparent viscosity in Eq. (16) can be affected by many factors, e.g., w/c, cement components, mixing time, standing time before measurement, and the degree of hydration17,18. When the cement components, the mixing time, the standing time before measurement, and the degree of hydration are determined, the viscosity is a function of the w/c
which means the parameters K and n might be the function of w/c as
Substituting Eq. (16) into Eq. (13), we have
In a tiny space increment \(\Delta y\) , it can be approximated as
According to Eq. (10), the solid volume fraction \(\phi\) is
Substituting Eq. (21) into Eq. (20), it is obtained
In the \(\Delta y\) , \(\frac{\partial {\dot{\gamma }}^{n}}{\partial y}\) is approximately considered as a constant. It means K can be determined by w/c in a tiny domain \(({w}_{0}/{c}_{0},w/c)\) , where \({w}_{0}/{c}_{0}\) is an arbitrary known water cement ratio in the concrete mix proportion of interest and \(w/c\) is close to \({w}_{0}/{c}_{0}\) . Then, \(K\) can be given by the Lagrange mean value theorem in the tiny domain
where \(\theta \in \left(0, 1\right)\) . Substituting Eq. (23) into Eq. (22), it is derived
In the increment \(\Delta t\) , \(w/c\) owns little changes, the right side of Eq. (26) can be approximated as.
Then Eq. (26) can be changed as
Equation (28) can be simplified as
In the tiny neighborhood of the \({w}_{0}/{c}_{0}\) , the following approximation is given
The values of \(\frac{\partial {\dot{\gamma }}^{n}}{\partial y}\) and \(\frac{\partial u}{\partial t}\) are the constants when the shear rate keeps constant. Substituting Eq. (32) into Eq. (30), we have.
By dividing both sides of Eq. (33) by \({w}_{0}/{c}_{0}\) , Eq. (34) is obtained as
Solving the ordinary Eq. (34), we have
in which \({CC}_{\mathrm{I}}\) is the calculation parameter, \({CC}_{\mathrm{II}}\) and b are
In the whole domain, the solution (36) is considered as the basic function. The Bernstein polynomial approximation29 can be concisely constructed as
where D1, D2 and D3 are the calculation parameters, and the detailed derivation of Eq. (38) is shown in Appendix I (see the Supplementary Information Appendix I for details). Substituting Eq. (38) into Eq. (19), it is derived
Then, n can be expressed as
which is the implicit solution of n. The solution is complicated and therefore hard for engineering applications. For simplification, the n is expanded in the neighborhood of \(w/c={w}_{0}/{c}_{0}\) as
and the first-order approximation is
It is taken the parameters as
When the value of \({f}_{0}\) is
which means the parameter \({C}_{2}=0\) , and n is
Due to the arbitrariness of \({w}_{0}/{c}_{0}\) , the general form of n is presented as.
Substituting Eq. (38) and Eq. (49) into Eq. (16), the four-parameter constitutive model is given as.
Equation (50) describes the functional relation between the shear stress and the shear rate of cement pastes with varying \(w/c\) . It should be noticed that the \(w/c\) can be replaced with \(w/b\) based on the assumption in Eq. (14).
This section discusses the apparent viscosity model with the SP which has the chemical structure as shown in Fig. 2. Its specific density is 1.07 g/cm3, with the side chain length (average number of ethylene oxide units) of 53. The average molecular weight is 58.2 × 103 g/mol, and the polydispersity index is 2.018.
The chemical structure of the main component of the SP18.
The SP mainly disperses cement particles through electrostatic repulsion and steric hindrance after adsorption30,31. The electrostatic repulsion force \({F}_{e}\) is
where \(k\) is the Coulomb’s constant, \({q}_{1}\) and \({q}_{2}\) are the charges of ions, \(r\) is the distance between two ions. When the type of SP is given, the \(k\) , \({q}_{1}\) and \({q}_{2}\) are known. The distance \(r\) is related to the number or the concentration of ions in a certain space, which means the whole electrostatic repulsion force \({F}_{et}\) is
in which \({n}_{S}\) is the dosage of the SP in a fluid domain \(\Omega\) . The steric hindrance describes how the physical structure of SP affects its ability to react. Its force \({F}_{s}\) can be approximated as a constant for the SP. Introducing \({F}_{et}\) and \({F}_{s}\) into Eq. (4), it is obtained
According to Eq. (52), \({F}_{et}\) is related to the dosage of the SP, which means \({F}_{n}\) in Eq. (54) is the function of \({n}_{S}\) . Substituting Eqs. (6)–(7) into Eq. (53), we have
In a tiny domain \(\left(y,y+\Delta y\right)\) , Eq. (55) can be approximated as
which can be changed into
When the influence of the SP is ignored, \({F}_{n}\) equals zero. The apparent viscosity \(\eta\) in Eq. (57) is reduced to Eq. (13), which can be approximated as the four-parameter model from Eq. (50).
Based on Eq. (57), it is considered that the influence of SP comes from the additional item \({\left[\frac{w/b}{\rho }\frac{\partial }{\partial y}\left(\frac{\partial u}{\partial y}\right)\right]}^{-1}\frac{1}{\phi {\rho }_{s}}{F}_{n}\) , which also can affect the other parameters in Eq. (58). A modified item is added to illustrate the effect of SP in Eq. (58) and the expression of the apparent viscosity is constructed as
where \({f}_{3}\left(w/b,{n}_{S}\right)\) is the additional modified item, in \(\left(y,y+\Delta y\right)\) , \({f}_{3}\) is
According to the Bernstein first-order approximation (see the Appendix I for details), considering the effect of SP to \({D}_{1}\) and \({C}_{1}\) , the viscosity of whole domain in an arbitrarily known y position is approximated as
in which \({d}_{1}\) and \({d}_{3}\) are the calculation parameters. \({F}_{1}\left({n}_{S}\right)\) is the effective coefficient of the dosage of the SP, and \({F}_{3}\left(w/b,{n}_{S}\right)\) owns the same dimension of \({\left[\frac{w/b}{\rho }\frac{\partial }{\partial y}\left(\frac{\partial u}{\partial y}\right)\right]}^{-1}\frac{1}{\phi {\rho }_{s}}{F}_{n}\) . The detailed expressions of \({F}_{1}\) , \({F}_{2}\) and \({F}_{3}\) are shown in Appendix II (see the Supplementary Information Appendix II for details).
In this section, the developed models are verified by being compared with the examples of the experimental results in rheological tests with different w/c and SP of cement pastes.
Three examples are given from the experimental results of Rosquo\(\ddot{\text{e}}\) t et al.22, Cyr et al.26, and Jeong et al.4 to verify the developed model with w/c.
In the experiment of Rosquo\(\ddot{\text{e}}\) t et al.22, the Portland cement CEM I 52.5 PM ES CP2 is used. The Bogue compositions of the cement are listed as follows: \({\text{C}}_{3}{\text{S}}\) 63.30%, \({\mathrm{C}}_{2}\mathrm{S}\) 17.90%, \({\mathrm{C}}_{3}\mathrm{A}\) 4.74%, and \({\mathrm{C}}_{4}\mathrm{AF}\) 5.62%22.All measurements used a Rheomat 115 rheometer with an MS 145 coaxial cylinder spindle. A thermal control system maintains a constant temperature in the test. Tests were conducted at 0.5 min after mixing and at a temperature of 20 ± 1 \(^\circ \mathrm{C}\) . For each w/c, the measurements are conducted by increasing the shear rate from 23 s−1 to 1200 s−1. For each shear rate, the shear stress was measured after a minimum of 30 s without fluctuation22. Four experimental testing points with w/c = 0.5 were chosen as the collocation points to solve the calculation parameters in Eq. (50) as
Compared with the other results in the experiment, the results of the developed model are shown in Figs. 3 and 4 In the two figures, the values of the apparent viscosity plunge before the shear rate reaches 100 s−1. Then, the descending rates slow when the shear rate further increases from 200 s−1 to 1000 s−1. The maximum viscosity is about 0.325 Pa·s for w/c = 0.4. The calculated results agree with the measured22.
The comparison between model calculation and the experimental data when w/c = 0.4, 0.5, 0.7.
The comparison between model calculation and the experimental data when w/c = 0.45, 0.60, 1.00.
The second verification is conducted based on the measured results from Cyr et al.26. The Bogue compositions of the cement used are listed as follows: \({\mathrm{C}}_{3}\mathrm{S}\) 60%, \({\mathrm{C}}_{2}\mathrm{S}\) 13%, \({\mathrm{C}}_{3}\mathrm{A}\) 10%, \({\mathrm{C}}_{4}\mathrm{AF}\) 6%, Gypsum 5%, and others 6%26. The apparatus used was a modified Rotovisco RV2 (Haake) with a six-blade vane. The pastes were mixed for 8.5 min to obtain a good dispersion of the components26. When the influence of gypsum and other ingredients is ignored, the value of w/c equals that of w/b. There are four experimental points of the curve with w/b = 0.3 used to solve the parameters in Eq. (50) as
The calculated results of the developed model and the measured results are shown in Fig. 5. The curve with w/b = 0.3 shows that the shear stress increases from 0 to 379 Pa with the shear rate increasing. The calculated results match the measured results26.
The comparison between model calculation and the experimental data when w/b = 0.30, .035, 0.41, 0.48.
The third verification is based on the test results from the experiment of Jeong et al.4. The chemical compositions of the cement used are listed as follows: \(\mathrm{CaO}\) 60.84%, \({\mathrm{SiO}}_{2}\) 13.26%, \({\mathrm{Na}}_{2}\mathrm{O}\) 10.05%, \({\mathrm{SO}}_{3}\) 3.59%, \({\mathrm{Al}}_{2}{\mathrm{O}}_{3}\) 3.42%, \({\mathrm{Fe}}_{2}{\mathrm{O}}_{3}\) 3.14%, \(\mathrm{MgO}\) 2.35%, \({\mathrm{K}}_{2}\mathrm{O}\) 1.17%, and others 2.18% (all by the weight percentage of the cement)4. The matrix was mixed at low speed for 1 min, halting for 1.5 min, and then was mixed for another 0.5 min at high speed. The matrix was placed in the rheometer right after the mixing. The protocol of the rheological test consists of pre-shearing for 0.5 min at a shear rate of 500 s−1. The pre-shearing procedure avoids memory effects, e.g., the thixotropic effect4. Four experimental results of the shear stress-shear rate curve with w/c = 0.6 are chosen as the collocation points to solve the calculation parameters in Eq. (50) as
Compared with the other measured results in the experiment of Jeong et al.4, the calculated results of the developed model are presented in Figs. 6 and 7.
The comparison between model calculation and the experimental data when w/c = 0.3, 0.35, 0.40, 0.45.
The shear stress increases as the shear rate increases, which is applicable for pastes with w/c from 0.3 to 0.6, as shown in Figs. 6 and 7. The maximum shear stress is around 2600 Pa, measured in the paste with w/c = 0.3 at the shear rate of 1000 s−1. The difference can be observed between the calculated results of the developed model and the measured results at low shear rates, especially for the paste with w/c = 0.3. Two reasons result in the difference. First, the test results of cement paste with w/c = 0.6 were used to calibrate the four parameters of the model. Based on these parameters, the rheological behaviours of the other pastes were predicted and compared with the experimental results. For the cement paste with w/c = 0.3, its rheological behaviour is different from the paste with w/c = 0.6. When the w/c = 0.3, the paste behaves like a solid-like system before the shear stress reaches the static yield stress of the paste. After the shear stress exceeds the static yield stress, the paste starts to flow. For the paste with w/c = 0.6, the paste is nearly liquid-like and it starts flowing without a clear-cut yield stress4. The parameters are calibrated from the paste with w/c = 0.6 and lead to the variation for the paste with w/c = 0.3 at the low shear rate. As the shear rate increases, the flow resistance from the yield stress of the paste accounts for diminishingly and the variation reduces.
The second reason lies in the measurement protocol in Jeong’s research. In the measurement, fresh samples were left standing for 0.5 min after the preshear procedure. The 0.5 min is regarded as the minimum time interval to wait between the preshear and the successive measurement4. Then, the rotational speed of the rheometer starts to increase from zero. The paste changed from a flow state in the pre-shearing to a static state at the end of the standing time. The standing time enables the paste with w/c = 0.3 to recover to a solid-like system with static yield stress. Based on the abovementioned two points, the variation between test results and the proposed models can be observed for the paste with w/c = 0.3.
The test results from the experiment of Liu et al.18 in Fig. 1 were used for another verification. The chemical compositions of the cement used are listed as follows: \(\mathrm{CaO}\) 63.80%, \({\mathrm{SiO}}_{2}\) 19.41%, \({\mathrm{Al}}_{2}{\mathrm{O}}_{3}\) 4.33%, \({\mathrm{SO}}_{3}\) 3.89%, \({\mathrm{Fe}}_{2}{\mathrm{O}}_{3}\) 2.91%, \(\mathrm{MgO}\) 1.29%, \({\mathrm{Na}}_{2}\mathrm{O}\) 1.29%, \({\mathrm{K}}_{2}\mathrm{O}\) 0.68%, and \({\mathrm{TiO}}_{2}\) 0.28% (all by the weight percentage of the cement)18.
Four groups of mixed proportions are chosen with different w/b and SP, which are presented in Table 1.
The apparent viscosity, defined as the ratio of shear stress to the shear rate, was measured using a Brookfield R/S SST2000 rheometer with Spindle CC25. Since the shear rate during concrete pouring was from 10 s−1 to 20 s−1, the maximum shear rate for the pastes was set to 25 s−118. After placing the paste into the rheometer, the sample was left to equilibrate for 0.5 min and then sheared at a constant rate of 25 s−1 for 1 min. The collocation points are chosen from the curve with SP = 2 and w/b = 0.16, and the calculation parameters are solved as
Compared with the measured results, the calculated results of the present model are shown in Fig. 8. Figure 8 shows the changes in the viscosity of cement pastes with different w/b and SP. When w/b = 0.32, the viscosity monotonically decreases with SP = 0.5% (the minimum SP dosage). The other curves in Fig. 8a change from flat to monotonically increasing as the increase of the shear rate. The same phenomena are observed in Fig. 8b and c. In Fig. 8d, all curves monotonically increase as the shear rate increases. The maximum value of the viscosity is 8.448 Pa·s with w/b = 0.16, SP = 2.2%, and a shear rate of 25 s−1. The results of the developed model have the same trend of viscosity and agree with the experimental data.
The comparison between model calculation and the experimental data for different w/b and SP. (a) paste with w/b = 0.32, (b) paste with w/b = 0.24, (c) paste with w/b = 0.20, (d) paste with w/b = 0.16.
The models can be used to calculate the viscosity of cement pastes and give technique support for the pumping process of cement pastes, e.g., long-distance pumping32, oil-well cementing14, and cementitious pastes grouting33. In these applications, the viscosity of the cementitious materials is important. With the models proposed in the paper, the effect of adjusting the water-cement ratio and adding superplasticizer on the viscosity of the pastes can be quantified and serves these applications.
The paper presents a mathematical model for the apparent viscosity of cement pastes with varying water-cement/binder ratios and the influence of polycarboxylate-based superplasticizers (SP). The following conclusions can be drawn:
(1) An ordinary differential equation is developed considering the w/b and parameter K in the Ostwald model. It is derived based on the Navier–Stokes equations and the Ostwald model of shear stress-shear rate relations;
(2) The whole domain constitutive relation of a four-parameter formula of cement pastes in the rheological stage is approximated by the first-order Bernstein polynomial approximation with different water-cement/binder ratios;
(3) An approximate expression of the viscosity for one type of polycarboxylate-based superplasticizer is constructed considering the result of electrostatic repulsion and steric hindrance;
(4) The developed models are verified by several rheological experiments with different water-cement/binder ratios and dosages of superplasticizers.
It should be noted that the developed models are applicable during the initial mixing stage of cement pastes. Further studies are needed to consider the influence of other factors such as the degree of hydration and the addition of aggregates and other additives on the rheological behaviour of cement-based mixtures.
The datasets used and/or analysed during the current study are available from the corresponding author upon reasonable request.
Liquid velocity in i (i = 1, 2, 3) direction
Body force on unit volume
Chidiac, S. & Mahmoodzadeh, F. Plastic viscosity of fresh concrete–A critical review of predictions methods. Cem. Concr. Compos. 31, 535–544 (2009).
Marar, K. & Eren, Ö. Effect of cement content and water–cement ratio on fresh concrete properties without admixtures. Int. J. Phys. Sci. 6(24), 5752–5765 (2011).
Gupta, M. The effect of water to cement ratio on early age properties of cement paste and mortar by UPV. Int. J. Eng. Sci. 5, 550–556 (2016).
Jeong, J., Chuta, E., Ramézani, H. & Guillot, S. Rheological properties for fresh cement paste from colloidal suspension to the three-element Kelvin–Voigt model. Rheol. Acta 59, 47–61 (2020).
Aïtcin, P. C. in Science and Technology of Concrete Admixtures Ch. The importance of the water–cement and water–binder ratios, 3–13 (Woodhead Publishing, 2016).
Einstein, A. A new determination of molecular dimensions, ETH Zurich, (1905).
Justnes, H. & Vikan, H. Viscosity of cement slurries as a function of solids content. Annu. Trans. Nord. Rheol. Soc. 13, 75–82 (2005).
Struble, L. & Sun, G. K. Viscosity of Portland cement paste as a function of concentration. Adv. Cem. Based Mater. 2, 62–69 (1995).
Hodne, H., Galta, S. & Saasen, A. Rheological modelling of cementitious materials using the Quemada model. Cem. Concr. Res. 37, 543–550 (2007).
Skare, E. L. et al. Rheology modelling of cement paste with manufactured sand and silica fume: Comparing suspension models with artificial neural network predictions. Constr. Build. Mater. 317, 126114 (2022).
Bingham, E. C. Fluidity and Plasticity. Vol. 2 (McGraw-Hill, 1922).
Li, M., Han, J., Zhou, Y. & Yan, P. A rheological model for evaluating the behavior of shear thickening of highly flowable mortar. Molecules 26, 1011 (2021).
Article CAS PubMed PubMed Central Google Scholar
Peng, J., Deng, D., Liu, Z., Yuan, Q. & Ye, T. Rheological models for fresh cement asphalt paste. Constr. Build. Mater. 71, 254–262 (2014).
Tao, C., Kutchko, B. G., Rosenbaum, E. & Massoudi, M. A review of rheological modeling of cement slurry in oil well applications. Energies 13, 570 (2020).
Campos, R. & Maciel, G. Test protocol and rheological model influence on determining the rheological properties of cement pastes. J. Build. Eng. 44, 103206 (2021).
Rubio-Hernández, F. J. Rheological behavior of fresh cement pastes. Fluids 3, 106 (2018).
Vipulanandan, C. & Mohammed, A. Smart cement rheological and piezoresistive behavior for oil well applications. J. Pet. Sci. Eng. 135, 50–58 (2015).
Liu, J. et al. Influence of superplasticizer dosage on the viscosity of cement paste with low water-binder ratio. Constr. Build. Mater. 149, 359–366 (2017).
Jones, T. & Taylor, S. A mathematical model relating the flow curve of a cement paste to its water/cement ratio. Mag. Concr. Res. 29, 207–212 (1977).
Lapasin, R., Longo, V. & Rajgelj, S. Thixotropic behaviour of cement pastes. Cem. Concr. Res. 9, 309–318 (1979).
Ivanov, Y. & Roshavelov, T. in Rheology of Fresh Cement and Concrete: Proceedings of an International Conference, Liverpool (ed Banfill P. G. F) 23–26 (Taylor & Francis).
Rosquoët, F., Alexis, A., Khelidj, A. & Phelipot, A. Experimental study of cement grout: Rheological behavior and sedimentation. Cem. Concr. Res. 33, 713–722 (2003).
Flatt, R. J. & Bowen, P. Yodel: A yield stress model for suspensions. J. Am. Ceram. Soc. 89, 1244–1256 (2006).
Ukrainczyk, N., Thiedeitz, M., Kränkel, T., Koenders, E. & Gehlen, C. Modeling SAOS yield stress of cement suspensions: Microstructure-based computational approach.Materials 13, 2769 (2020).
Article ADS CAS PubMed PubMed Central Google Scholar
Ma, S. & Kawashima, S. A rheological approach to study the early-age hydration of oil well cement: Effect of temperature, pressure and nanoclay. Constr. Build. Mater. 215, 119–127 (2019).
Cyr, M., Legrand, C. & Mouret, M. Study of the shear thickening effect of superplasticizers on the rheological behaviour of cement pastes containing or not mineral additives. Cem. Concr. Res. 30, 1477–1483 (2000).
Kundu, P. K., Cohen, I. M. & Dowling, D. R. Fluid Mechanics 6th edn, Vol. 1 (Academic Press, 2015).
Derksen, J. Simulations of dense agitated solid–liquid suspensions—Effects of the distribution of particle sizes. Chem. Eng. Sci. 189, 56–64 (2018).
Farouki, R. T. The Bernstein polynomial basis: A centennial retrospective. Comput. Aided Geom. Des. 29, 379–419 (2012).
Faraday, M. V. Experimental researches in electricity. Philos. Trans. R. Soc. Lond. 122, 125–162 (1832).
Yoshioka, K., Sakai, E., Daimon, M. & Kitahara, A. Role of steric hindrance in the performance of superplasticizers for concrete. J. Am. Ceram. Soc. 80, 2667–2671 (1997).
Feys, D., De Schutter, G., Fataei, S., Martys, N. S. & Mechtcherine, V. Pumping of concrete: Understanding a common placement method with lots of challenges. Cem. Concr. Res. 154, 106720 (2022).
Anagnostopoulos, C. A. Effect of different superplasticisers on the physical and mechanical properties of cement grouts. Constr. Build. Mater. 50, 162–168 (2014).
The authors of this paper are indebted to the financial support from the National Key R&D Programme Intergovernmental Cooperation on International Science and Technology Innovation (No. 2021YFE0114100), the Sino-German Center for Research Promotion (No. GZ 1574), Science and Technology Commission of Shanghai Municipality (No. 21DZ1203505), and Science and Technology Project of Department of Transportation of Jiangxi Province (No. 2021C0008).
College of Civil Engineering, Tongji University, Shanghai, 200092, China
Yong Yuan & Xiaoyun Wang
College of Civil Engineering and Architecture, Jiaxing University, Zhejiang, 314001, China
Institute of Construction and Building Materials, Technical University of Darmstadt, 64287, Darmstadt, Germany
Peng Xiao & Eduardus Koenders
School of Aerospace Engineering and Applied Mechanics, Tongji University, Zhangwu Road 100, Shanghai, 200092, China
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
Y.Y.: Conceptualization, Methodology, Writing—original draft. X.W.: Methodology, Verification, Writing—review & editing. X.C.: Methodology, Verification, Visualization. P.X.: Writing—review & editing. E.K.: Project administration. Y.D.: Supervision.
The authors declare no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
Yuan, Y., Wang, X., Chen, X. et al. Mathematical models of apparent viscosity as a function of water–cement/binder ratio and superplasticizer in cement pastes. Sci Rep 13, 22301 (2023). https://doi.org/10.1038/s41598-023-48748-4
DOI: https://doi.org/10.1038/s41598-023-48748-4
Anyone you share the following link with will be able to read this content:
Sorry, a shareable link is not currently available for this article.
Provided by the Springer Nature SharedIt content-sharing initiative
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.
Scientific Reports (Sci Rep) ISSN 2045-2322 (online)
![Mathematical models of apparent viscosity as a function of water–cement/binder ratio and superplasticizer in cement pastes | Scientific Reports](http://cdn.globalso.com/water-reducer/44.jpg)
High Range Water Reducer Sign up for the Nature Briefing newsletter — what matters in science, free to your inbox daily.