Wain R and Whitty J*
School of Engineering, UK
*Corresponding author: Whitty J, John Tyndall Institute, School of Engineering, Preston, PR1 2HR, UK
Submission: November 02, 2018; Published: November 19, 2018
ISSN: 2576-8816
Volume6 Issue5
A reasonably robust procedure is presented to evaluate each of the parameters of four non-Newtonian constitutive relationships for human blood, namely Cross, Carreau-Yasuda and modifications to Ostwald-de Waele and Sisko fluids. For each of the rheological models presented herein, a multidimensional optimization routine is used to find salient constitutive parameters from a compilation of digitized experimental data evident in the literature. It is shown for three of the models that, to conserve structural identification, the so-called low shear viscosity term can be set as it has little physiological relevance. The method presented herein is shown to minimize the square of the errors between the four suggested constitutive relationships and empirical data. It is shown that, for the data set investigated here, parameters which had previously been assumed to be fluid properties exhibit different values depending on the constitutive relationship selected.
Blood is a suspension of various solid cells in a liquid plasma, formed from a solution of proteins and minerals. The rheological properties of blood are influenced by individual characteristics of these components as well as factors such as haematocrit, amongst others [1,2]. Whilst plasma itself is a Newtonian fluid, addition of the blood cellular components alters its rheological behaviour towards a shear-thinning non-Newtonian fluid [1,2]. No single rheological model can capture all aspects of this complex fluid, and as such particular aspects can be modelled individually [2,3]. To this end, differing rheological models are applied when predicting flow characteristics depending on area of the circulatory system of interest. For example, blood is often approximated to a Newtonian fluid in arterial flows, though has also been commonly modelled using the Cross [4–7] and Carreau-Yasuda [5–11] models; their constitutive relationships being defined respectively by:
where, Φ(γ˙) is the so-called effective fluid viscocity, μ_{o} is the low shear viscosity, μ_{∞} (∼2-4cP) is the high shear viscosity, λ is the shear time constant and a is a shaping constant termed the Yasuda index. Whilst n∈(0,1) is the power index with values less than unity indicating shear thinning behavoiur, ergo values approaching unity approximate Newtonian behaviour. In this work we have used multivariate optimization routines resident in the SciLab data driven modelling open source software to approximate each of these. The values obtained being compared with the more traditional non- Newtonian constitute relations proposed by Oswald-de Waele [12] & Sisko [13] defined respectfully as:
where K is the viscosity consistency and μ_{l} is a low (non-zero) shear viscosity. To ensure dimensional consistency, each of these are presented here in slightly modified forms to those initially described. As such, they should be more correctly referred to as to as Modified Oswald [12] and Modified Sisko [13] fluids respectfully.
Parameters were approximated for each of the before mentioned rheological models using compiled experimental data from the literature [14-17]. This required a multivariate optimization of the square of the error function, viz:
Where μ_{m} is the mth measurement of the blood viscosity, Φ(γ˙) is constitutive relationship, i.e.: eqn. (1), (2), (3) or (4), and N is the total number of data points. For each of the rheological models we simultaneously fit the free parameter set to experimental data using a variation of the Nelder-Mead optimization algorithm (i.e. fminbnd routine^{1}) in SciLab open-source data driven modelling code [18]. To preserve structural identification of the Modified Sisko and Carreau-Yasuda fluid models the low shear viscosity terms (μ_{l} and μ0 respectively) were fixed to 150cP [8,11,19]. The literature maintaining that these terms have no physiological relevance so the calculations that follow are unaffected by this assumption. Cross fluid bounds were used as [1,300], [0.01,1] and [0.01,100] for λ, n and μ_{∞} respectively. For the Modified Sisko and Carreau-Yasuda models identical bounds
^{1}This being resident in the FOSSEE optimization toolbox developed at the University of Bombay:
were used for these parameters as for the Cross fluid. In the case of the latter, an additional bound [0.01,5] was used to represent the Yasuda index (a). In the case of the Oswald-de Waele fluid, the value of K in equation (3) is effectively set (to (150 − 2.3) ≈ 148cP) in order to render a structurally identifiable model with two free parameters. Refined accuracy being obtained (ipso facto verification) via application of an analogous constrained multivariate optimization algorithm. That is, the fmincon routine with initializations set to those previously evaluated from the fminbnd SciLab command. Each of the calculations allowed for 30000 interactions in a time of 16.7 minutes. A convergence tolerance of 10-10 being set for the initial fminbnd calculation with default values being left for the fmincon evaluation which followed.
In most of the cases investigated, the optimization routine reported that either an optimum solution was found or that a solution was obtained to acceptable accuracy. The noteable exception being in the case of the Carreau-Yasuda fluids which initially reported that the maximum number of iterations had been reached. In this case the bounds we altered using successive bisection until a solution to acceptable accuracy was obtained. These values were then used as initialization within the aforementioned SciLab fmincon routine thereby obtaining an optimized solution.
Table 1 compares these calculated parameters with those in the literature [1,3,4,20-23]. Those for Modified Sisko fluid are similar to the Ostwald-de Waele fluid, on the other hand equivalent parameters calculated for the Cross fluid are similar to those of Carreau-Yasuda model. The power index and high shear viscosity values are in good agreement with those quoted in the literature [4,24-26]. Particularly good agreement is demonstrated with respect to values obtained from the Cross and Carreau-Yasuda constitutive relationships and those reported in the literature. This said, analogous values obtained using identical optimization routines for the Modified Ostwald-de Waele and Modified Sisko fluids indicate significantly different values. However, similar values are obtained for the power index and shear time constant when comparing the Modified Sisko and Modified Ostwald-de Waele fluids.
Table 1:Rheological parameters to 2 s.f. obtained from fitting combined existing experimental data shown in Figure 1 to various non-Newtonian models. Experimental data is taken from a variety of sources.
*set value
Discrepancies in the power index calculations are of particular interest from a practical viewpoint. To date the power index has been thought of as a diversion from an idealized Newtonian fluid. The natural assumption being that this quantity is a property of the fluid, hence independent of constitutive relationship, these results certainly suggest that this is not the case. Implying that, at least in the current mathematical forms, the power index is not a property of the fluid.
Figure 1:A log-log plot demonstrating four non-Newtonian constitutive models fitted to collated experimental rheological data for blood adapted from the literature.
Figure 1 demonstrates that the empirical data are well predicted by each of the models to varying degrees; the primary source for these data being Ballyk et al. [20]. Crosses represent experimental blood rheological data collated from several sources i.e. Huang [14], Merrill [15], Cokelet [16] & Skalak [17]. Lines represent rheological models i.e. Ostwald [12], Cross [27], Carreau [28,29], and Modified Sisko [13], fitted to the compiled experimental data.
As expected each of the models converge at higher (imposed) apparent viscosity values. The models predicting the apparent viscosity of the blood at relatively lower shear rates γ ≤ 50s−1 very well. At higher shear rates the model predictions diverge. At the highest measured shear rate 950s−1 the Cross and Carreau-Yasuda models predict higher viscosity values of respectfully 3.77cP and 3.69cP, with the Sisko model predicting a 12% lower viscosity than that the measurement of 3.6cP. Moreover, Figure 1 demonstrates a deficiency with the Ostwald-de Waele model whereby at high shear rates the apparent viscosity is too low in comparison with the empirical data. Indeed, since this model is unbounded at both low and high shear limits unlike other shear-thinning models discussed, its use for evaluating blood is limited to a narrow range of shear rates [30], in this case probably between 1s−1 and 45s−1 (Table1).
A combination of optimization routines has been used to evaluate the constitutive non-Newtonian relationships from compilation of data evident in the literature, via minimization of appropriate least squares error. It was found that the values within each constitutive relationship vary depending on the data set used and non-Newtonian fluid model used. Specifically, this study has shown that:
A. Using these calculated optimized parameters most of the non-Newtonian models describe the empirical data rather well;
B. The Oswald-de Wael model is not suitable for predicating the effective blood viscosity at high shear rates;
C. Due to the similarity of the models the Sisko and Ostwaldde Waele similar constitutive non-Newtonian parameters are apparent; analogous similarities being evident regarding Cross and Carreau fluids;
D. The optimization routine predicts the Carreau-Yasuda and Cross models agree well about analogous parameters as well as values in the literature;
E. The non-Newtonian constitutive parameters (e.g. n,λ and μ_{∞}and a) may not be properties of blood.
We should also note that there are many more data sets available in the literature. Furthermore, much more advanced technology is available to measure blood viscosity. So, whilst the work herein does raise some interesting observations which have been properly verified, further work is of course required to validate these findings.
© 2018 Whitty J. This is an open access article distributed under the terms of the Creative Commons Attribution License , which permits unrestricted use, distribution, and build upon your work non-commercially.