Next Article in Journal
Teaching and Learning Floating and Sinking: Didactic Transformation in a Density-Based Approach
Next Article in Special Issue
Low-Noise Synthetic Turbulence Tailored to Lateral Periodic Boundary Conditions
Previous Article in Journal
Oscillating Magnetohydrodynamic Stokes Flow between Porous Plates with Spatiotemporally Periodic Reabsorption
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Precise Method to Estimate the Herschel-Bulkley Parameters from Pipe Rheometer Measurements

1
Total, R&D, 4029 Stavanger, Norway
2
Energy Department, Norwegian Research Centre, 4021 Stavanger, Norway
*
Author to whom correspondence should be addressed.
Fluids 2021, 6(4), 157; https://doi.org/10.3390/fluids6040157
Submission received: 24 March 2021 / Revised: 6 April 2021 / Accepted: 13 April 2021 / Published: 14 April 2021
(This article belongs to the Special Issue Complex Fluids and Flows: Algorithms and Applications)

Abstract

:
Accurate characterization of the rheological behavior of non-Newtonian fluids is critical in a wide range of industries as it governs process efficiency, safety, and end-product quality. When the rheological behavior of fluid may vary substantially over a relatively short period of time, it is desirable to measure its viscous properties on a more continuous basis than relying on spot measurements made with a viscometer on a few samples. An attractive solution for inline rheological measurements is to measure pressure gradients while circulating fluid at different bulk velocities in a circular pipe. Yet, extracting the rheological model parameters may be challenging as measurement uncertainty may influence the precision of the model fitting. In this paper, we present a method to calibrate the Herschel-Bulkley rheological model to a series of differential pressure measurements made at variable bulk velocities using a combination of physics-based equations and nonlinear optimization. Experimental validation of the method is conducted on non-Newtonian shear-thinning fluid based on aqueous solutions of polymers and the results are compared to those obtained with a scientific rheometer. It is found that using a physics-based method to estimate the parameters contributes to reducing prediction errors, especially at low flow rates. With the tested polymeric fluid, the proportion difference between the estimated Herschel-Bulkley parameters and those obtained using the scientific rheometer are −24% for the yield stress, 0.26% for the consistency index, and 0.30% for the flow behavior index. Finally, the computation requires limited resources, and the algorithm can be implemented on low-power devices such as an embedded single-board computer or a mobile device.

1. Introduction

Accurate characterization of non-Newtonian fluids is a critical issue in a wide range of domains from chemical to physical processing since it governs the process efficiency, safety, and end-product quality. Usually, the rheological properties are determined by sampling the fluid in the production system and measuring its characteristics with an off-line rheometer. However, the characterizations made off-line may be unsatisfactory when the viscous properties of the fluid change continuously. With non-Newtonian fluid, the in-line measurement of viscous properties is much more complicated than for Newtonian fluids where simple methods such as the tuning fork technique, for example, can be implemented even in very harsh conditions [1]. To meet the requirements of the industry for inline measurement of the viscous properties of non-Newtonian fluids, several methods have been developed. Most of them rely on the measurement of the fluid radial velocity profile losses along a known section of pipe. Different techniques were developed to estimate the velocity profile such as magnetic resonance imaging [2], laser velocimetry [3], ultrasound velocity profiling [4], and electrical resistance sensing [5]. Each technique comes with its own limitations due to the physical principles of the sensor, specific environmental constraints, cost, and/or safety. For instance, methods based on the measurement of the velocity profile are not adapted to investigate fluids that are opaque or which structure results in the diffraction of the probing signal that is sent by the radiating source, like for example with emulsions or fluids containing many solid particles.
An alternate approach to the one based on velocimetry consists of measuring the pressure gradient for various bulk fluid velocities while circulating the fluid in pipes [6]. This method has the advantage of working with fluids that are completely opaque to or that cause a high degree of diffraction of the radiations from possible illuminating sources. The inconvenience of the method is that measurements shall be made with multiple flow conditions in order to acquire enough information for calibrating a chosen rheological model. The choice of a rheological model depends on the type of fluid being measured. In turn, the minimum number of flow conditions is a function of the chosen rheological model. The simplest non-Newtonian rheological models have two parameters like the Bingham plastic or the power law models, while more realistic models have three parameters such as the Herschel-Bulkley [7], Robertson-Stiff [8], or Heinz-Casson models [9]. More complex models such as the ones from Carreaux [10] and Quemada [11] depend on four parameters.
The Herschel-Bulkley model is often used to describe the rheological behavior of non-Newtonian fluids, simply because it is a generalization of both the Bingham and power-law models. Therefore, we have chosen to use that model to describe the rheological behavior of the fluid being measured in the pipe rheometer. From differential pressure measurements, we derive a calculation method for the Hershel-Bulkley parameters using a combination of physics-based equations and non-linear optimization. The present method can accurately estimate the fluid rheologic characteristics based on pressure loss measurements and three hypotheses: the fluid rheological behavior can be approximated by the Herschel-Bulkley model, measurements are made in laminar flow conditions and slippage at the wall is negligible.

2. Herschel-Bulkley Flow in a Pipe

The Herschel-Bulkley (H-B) model, also referred to as the yield power law (YPL) model, is one of the many relations used to describe how certain viscous fluids behave; especially fluids that exhibit yield stress and for which the shear stress tends to behave like a power law at high shear rate.
The model relates the shear stress τ to the shear rate γ ˙ as follow [7]:
τ = τ γ + K γ ˙ n ,
where τ γ is the yield stress, K is the consistency index and n is the flow behavior index. Note that the model reduces to the Bingham plastic rheological behavior when n = 1 , i.e., τ = τ γ + K γ ˙ and the power law model when the yield stress is equal to zero, i.e., τ = K γ ˙ n .
For parameter identification purposes, a linearized expression of Equation (1) is often used [12]:
ln τ τ γ = ln K + n ln γ ˙ .

2.1. Practical Estimation of the Herschel-Bulkley Parameters

The Herschel-Bulkley parameters for a particular fluid are usually found experimentally by measuring several pairs of data τ ,   γ ˙ with a rheometer. Then, using an optimization method such as a direct search or a Newton-like method, the corresponding model is fitted, and its parameters are estimated. However, as shown by Kelessidis et al. [12], these methods often give a negative yield stress, τ γ , which is not physically meaningful. To overcome this challenge, different methods were developed which use the linearized model (Equation (2)) and require a prior estimation of τ γ . To fulfill this requirement, multiple approaches were tried: trial-and-error [13], graphical extrapolation of the rheogram [14], or golden section search [12], involving constrained optimization algorithms when needed. It is also possible to use machine learning algorithms, more specifically genetic algorithms, to directly determine the H-B parameters of a fluid, as shown in [15].
More recently, Mullineux showed in [14] that the resulting value of parameters K and γ ˙ were very sensitive to the initial estimation of τ γ . The problem is that the estimation of τ γ is heavily influenced by measurements at the lowest flow rates where uncertainty is the highest. Mullineux showed that for a fluid described by the following H-B parameters τ γ = 5   P a , K = 4   P a . s n and n = 0.35 , an error of 1 Pa on the estimation of τ γ leads to an error of 12.5% on K and 6% on n . Thus, he described a method to estimate the flow behavior index, n, independently from the determination of the two other parameters, therefore, allowing to estimate K and τ γ by simple linear regression. This method is especially interesting as it does not require extrapolation on uncertain data and returns the global optimum value for n.
When looking for the best H-B coefficients to model a fluid, an optimization step is usually required, and the method relies on the minimization of the least square error:
S = i = 1 m τ γ + K γ ˙ i n τ i 2 ,
where, τ i ;   γ ˙ i is one of the m measured pairs of shear stress and shear rate. In [14], the author showed that in the case of Herschel-Bulkley fluids, minimizing S is equivalent to solving the following equation:
F n = m γ ˙ i n γ ˙ i n ln γ ˙ i γ ˙ i n γ ˙ i 2 n γ ˙ i n ln γ ˙ i τ i γ ˙ i n τ i γ ˙ i n τ i ln γ ˙ i = 0 .
It has been shown in [14] as well that F n has only one root for n   0 ; + which is the n that minimizes the function described by Equation (3). n is found using a simple root search method such as the bisection method, K and τ γ are then determined by performing a linear regression based on Equation (1) on the measured data. The advantage of this method is that it allows us to find the H-B parameters corresponding to the global minimum of S at a very low computational cost.

2.2. Laboratory Flow Loop

To test the described method of the present study, experiments have been conducted in a laboratory flow loop. A positive displacement pump circulates fluid stored in a thermo-regulated tank with internal circulation. First, the fluid flows into a pipe with an internal diameter of 25 mm and then into a glass tube with an internal diameter of 15.5 mm. Along the glass tube, there are three differential pressure sensors. The distances between the ports of the differential pressure sensors are 0.209 m, 0.212 m, and 0.206 m, from upstream to downstream, respectively. The first sensor is placed 1.955 m away from the change in diameter from 25 mm to 15.5 mm. The second sensor is positioned 1.415 m apart from the first one and the last sensor is separated from the previous by 1.46 m. A schematic representation can be found in figure 8 of [16], reproduced here as Figure 1. There is also a pressure-reducing valve placed 1.5 m upstream from the change in diameter. The purpose of this valve is to ensure that the fluid is fully sheared before it continues its journey inside the flow loop. In the case of non-thixotropic fluids, one sensor would have been enough, but since three series of measurements were available, we decided to use all of them to increase the amount of data.
All the experiments were conducted at an average controlled temperature of 27.9 °C. The temperature variations during the experiments have not exceeded ±0.6 °C. The first measurements have been made with an aqueous solution of Carbopol 980 having a density of 997 kg/m3 at 25 °C. This is a non-Newtonian fluid showing no signs of thixotropy and is well described by the Herschel-Bulkley model.

2.3. Rheometer and Data Correction

To confirm the results obtained with the pipe rheometer, a scientific rheometer was used (Physica MCR 301, manufactured by Anton Paar, headquarters: Graz, Austria). In a concentric rheometer configuration, the rotational speed is converted to a shear rate at the wall and the shear stress at the wall is deduced from the measured torque. However, the relationship between rotational speed and the shear rate at the wall is constant only for Newtonian fluids. For non-Newtonian fluids, the relationship is more complex and requires the utilization of all the measurements of shear rate and shear stresses in order to be estimated [17]. Similarly, the ratio of the torque to the shear stress is constant only for a Newtonian fluid. For non-Newtonian ones, corrections for end-effects shall be applied [18].
For this paper, we have applied a correction method based on Tikhonov regularization as described in [19]. Figure 2 shows the rheometer measurements and their corrected values for an aqueous solution of Carbopol 980 (Lubrizol, headquarters: Wickliffe, OH, USA).
After applying the correction on rheometer data and using the method from Mullineux for Herschel-Bulkley parameter identification, the following values for this fluid are found: n R h e o = 0.6389 ,   K R h e o = 0.2717   Pa . s n ,     τ γ ,   r h e o = 1.198   Pa .

3. Pipe Rheometer

The objective is to develop an algorithm able to estimate the H-B parameters of a fluid flowing in a pipe using only measurements of pressure loss along the pipe and the flow rate. Vajargah et al. [20,21] developed a similar method to estimate the rheology of drilling fluid during an oil well drilling operation, using only pressure loss measurements scattered along the drill string. In [21] they developed a pipe rheometer to measure the drilling mud characteristics on the rig surface. However, the authors used generic nonlinear optimization techniques to fit the rheological models. Other devices were developed for in-line fluid rheology measurements like a helicoidal pipe rheometer [22] to save space, or on-line dedicated sensors [23]. However, because of the high complexity brought by the geometry or the nature of the sensors, data was interpreted using a machine learning algorithm that needs to be trained and maintained thus raising the complexity of the apparatus and reducing reliability in an operational environment.

3.1. Background

Given the reasons evoked in Section 2.1, it is advisable to use Mullineux’s H-B parameters identification method instead of a regular nonlinear optimization technique. Thus, we need to calculate pairs of τ ,   γ ˙ from our measured Q ,   d P d L , where Q is the volumetric flowrate and d P d L is the pressure gradient.
In laminar flow, the wall shear stress τ w is directly proportional to the pressure gradient d P d L :
τ w = R 2 d P d L ,
where R is the inner pipe radius.
To calculate the wall shear rate, we can use the Weissenberg-Rabinowitsch-Mooney (WRM) relation for a pipe developed in [24]. The WRM equation is mainly based on assumptions on the flow instead of the fluid and therefore can be applied to non-Newtonian fluids:
γ ˙ w = 3 4 γ ˙ N , w + 1 4 γ ˙ N , w d ln γ ˙ N , w d ln τ w ,
with γ ˙ N , w = 4 Q π R 3 being the wall shear rate in the case of a Newtonian fluid, and called apparent shear rate at the wall.
The main difficulty in calculating γ ˙ w is that we need to estimate the quantity d ln γ ˙ N , w d ln τ w . Most of the time this quantity is evaluated by plotting ln γ ˙ N , w as a function of ln τ w , applying a linear regression and using the slope as the derivative. However, this assumption is only true if τ γ = 0 or in other words if we have a Power Law fluid. When a linear fit is not enough, it is generally advised to perform a polynomial fit [25]. However, in practice, a least-square polynomial fit tends to be inaccurate. Indeed, when collecting a lot of measurements, minor imprecisions of sensors or external perturbations produce a cloud of point which creates ambiguity when fitting a polynomial on the measurements. In addition, the polynomial fit is chosen for its ease of use and not because it corresponds to physical reality. In order to get a more accurate regression, we will derive a physics-based formula for d ln γ ˙ N , w d ln τ w :
First, let us express d ln γ ˙ N , w d ln τ w as a function of the flow-rate Q:
d ln γ ˙ N , w d ln τ w = d ln 4 v R d ln τ w = d ln 4 Q π R 3 d ln τ w = d ln Q + ln 4 π R 3 d ln τ w = d ln Q d ln R 2 d P d L ,
Using the equation developed by Kelessidis et al. in [12], we have:
Q = π n K 1 / n τ w τ γ 1 + 1 n τ w R 3 × τ w τ γ 2 1 + 3 n + 2 τ γ τ w τ γ 1 + 2 n + τ γ 2 1 + n
Let´s express ln Q from Equation (8):
ln Q = ln π n K 1 / n + n + 1 n ln τ w τ γ 3 ln τ w + 3 ln R + ln 1 + n 1 + 2 n τ w τ γ 2 + 1 + n 1 + 3 n 2 τ γ τ w τ γ + 1 + 2 n 1 + 3 n τ γ 2 ln 1 + n 1 + 2 n 1 + 3 n .
By regrouping the constants, we obtain:
ln Q = n + 1 n ln τ w τ γ 3 ln τ w + ln C 2 τ w 2 + C 1 τ w + C 0 + A ,
with A = ln π n K 1 / n + 3 ln R ln 1 + n 1 + 2 n 1 + 3 n , C 0 = 6 n 4 n 2 τ γ 2 , C 1 = 2 n τ γ 1 + n and C 2 = 1 + n 1 + 2 n .
Let us now express ln Q as a function of ln τ w :
ln Q ln τ w = n + 1 n ln e ln τ w τ γ 3 ln e ln τ w + ln C 2 e 2 ln τ w + C 1 e ln τ w + C 0 + A ,
Finally, we obtain:
d ln γ ˙ N , w d ln τ w = d ln Q d ln τ w = n + 1 n τ w τ w τ γ C 1 τ w + 2 C 0 C 2 τ w 2 + C 1 τ w + C 0 1 .

3.2. Experimental Determination of d ln γ ˙ N , w d ln τ w

As stated in Section 3.1, in order to compute accurately the derivative quantity (Equation (12)) it is better to avoid relying on curve fitting, however, a prior estimation of the H-B parameters is required in that case. The solution used here is to perform a nonlinear least square curve fitting with the Levenberg-Marquardt method [26] to fit Equation (8) on the measured pair data Q ,   d P d L and using Equation (5) to calculate τ w from   d P d L . The Levenberg-Marquardt algorithm has been developed to solve non-linear least square curve fitting. It is a more robust method than more classical methods such as the Gauss-Newton and gradient descent algorithms, but the method can only find a local minimum and not necessarily the global one. Therefore it is an important point to choose values of the initial parameters that are close to the optimum. Consequently, it is very important to find a physically-based method to estimate those initial parameters.
An estimation of n is obtained by calculating the limit of ln Q when τ w goes toward infinity:
lim ln τ w ln Q = lim τ w ln Q = lim τ w 1 n ln τ w = + .
The trend for τ w of ln Q is a line of slope equal to 1 n so we can do a linear regression on the last part of the curve ln(Q) as a function of ln( τ w ) , read the slope and use it to estimate n (see Figure 3). This behavior is easily explained by the fact that Herschel-Bulkley fluids are power-law fluids with a yield stress, thus, when τ w gets large, the influence of the yield stress disappears, and the behavior of the fluid can be described with a power law.
An initial estimation of the yield stress, τ γ , is obtained from the wall shear stress at a low flow rate. In Figure 4a, it is visible that there is a range of d P d L 0 ; 1150   P a . m 1 where either the fluid flows or it does not dependent on the measurements. This is due to the gelling property of the fluid: when left at rest, the fluid acts like an elastic solid. After a gelling period, there is a stress overshoot while increasing the shear strain. This behavior has been commonly observed on Carbopol fluids at different concentrations [27]. This stress overshoot depends on the chemical composition of the fluid but also the resting time. τ γ is always smaller than the stress overshoot, therefore it is between the smallest value of the wall shear stress measured under viscous flow and the highest value corresponding to gel breaking. It is possible to obtain a rough estimation of τ γ by taking a weighted average of the shear stress between the first flowing and the last non-flowing τ w (shown in Figure 4b).
τ γ is estimated using a weighted average of the results shown in Figure 4b:
τ γ = 1 m i = 1 m Q > 0 | τ w , i ε 2 τ w < τ w , i + ε 2 τ w , i ,
where m is the number of bins created for the histogram, ε is the shear stress span corresponding to a bin, and i the index to run through the list of bins. The term Q > 0 | τ w , i ε 2 τ w < τ w , i + ε 2 is the probability that there was a fluid movement within the conditions defined by the current bin, in other words when the wall shear stress belongs to a given range of values. In addition, it is required to exclude data points whose wall shear stress is above the last non-flowing wall shear stress because it is known that the yield stress value is below the stress overshoot value. With the above-described method, we estimate an initial guess for the yield stress of τ γ = 1.19   Pa , which is an acceptable value when compared to the one obtained with the scientific rheometer measurements, i.e., τ γ ,   r h e o = 1.20   Pa .
By rearranging Equation (8), we obtain an expression of the consistency index, K , that depends on Q ,   d P d L ,   n and τ γ :
K = π n Q τ w τ γ 1 + 1 n τ w R 3 × τ w τ γ 2 1 + 3 n + 2 τ γ τ w τ γ 1 + 2 n + τ γ 2 1 + n n ,
The consistency index is then estimated by randomly sampling the measurements points in the set of pairs Q ,   d P d L and applying Equation (15) using the previously estimated values for n and τ γ in order to determine K . The median value of all the calculations is taken in order to reduce the sensitivity of the estimation to outliers (see Figure 5).
The estimated Herschel-Bulkley parameters are now used as initial values for the Levenberg-Marquart optimization algorithm. The objective is to obtain an estimation for the Herschel-Bulkley parameters that are sufficiently precise to calculate d ln γ ˙ N , w d ln τ w using Equation (12). The estimated parameters using the Levenberg-Marquart optimization are: τ γ ,   L M = 1.03   P a   , K L M = 0.256   P a . s n , n L M = 0.649 .
Now, combining Equation (12) with the Weissenberg-Rabinowitsch-Mooney relation, i.e., Equation (6), it is possible to calculate the wall shear rate. Then by applying the calculated pairs τ w ,   γ ˙ w   , the method described by Mullineux, (see Section 2.3), the Herschel-Bulkley parameters are obtained based on the flow-loop measurements: τ γ ,   M = 0.909   P a   , K M = 0.272   P a . s n , n M = 0.641 . It is interesting to note that the consistency and flow behavior indices are closer to the ones obtained with the scientific rheometer measurements compared to those obtained with the Levenberg-Marquart optimization, K M K R h e o K R h e o 0.26 % , K L M K R h e o K R h e o 5.85 % , n M n R h e o n R h e o 0.30 % , n L M n R h e o n R h e o 1.66 % while the new estimation of yield stress deviates more than the value obtained with the Levenberg-Marquart optimization, τ γ ,   M τ γ ,   r h e o τ γ ,   r h e o 24.12 % ,
τ γ ,   L M τ γ ,   r h e o τ γ ,   r h e o 14.02 % .

3.3. Results

Figure 6 presents the measured and calculated pressure losses as a function of the flow rate. The pressure losses curves obtained with the Levenberg-Marquart and Mullineux method are so close that they appear superimposed. Looking at the Herschel-Bulkley parameters produced by each method, only the yield stress varies with a substantial amount (approximatively −24% for the Mullineux method, and −12% for the Levenberg-Marquart optimization method, based on the pressure gradient measurements compared to the value obtained with the scientific rheometer). This result is to be expected, knowing that the measurements at low flow rate are scattered over a large range of pressure loss values. In addition, the low yield stress of this sample makes the determination of τ γ even more sensitive to numerical and experimental noise. The resulting curves presented in Figure 6 are still in agreement with the measurements. One can notice that at a low flow rate the curves produced from the Levenberg-Marquart optimization and Mullineux methods are closer to the mean value of the measurements than the one obtained based on the Herschel-Bulkley parameters calibrated with the scientific rheometer measurements.
In Figure 7, the estimation made using Equation (12) shows a hyperbolic shape when the estimations made with a polynomial fit of the measured data show an almost flat curve. It is important to note that the hyperbolic shape is amplified with an increase of the yield stress τ γ value of the fluid. Therefore it shows that it is not possible to estimate d ln γ ˙ N , w d ln τ w using a polynomial fit of ln γ ˙ N , w ln τ w as it will result in errors while predicting the rheological behavior of the fluid at low shear stress (and therefore at low flow rate). Using a proper estimation of d ln γ ˙ N , w d ln τ w becomes more and more critical as the yield stress of the fluid increases.

4. Discussion

As shown in Figure 6, the parameters obtained with this algorithm are very close to those obtained with the scientific rheometer after correcting the measurements for non-Newtonian effects. However, it is important to mention that such results can only be obtained if the algorithm is applied on points that follow the Herschel-Bulkley behavior: it is required to remove outliers, remove measurements made at transitional and turbulent flow regimes and remove low flow rates data where gelling occurs. The method used here was to analyze the impact of the removal of data from the dataset. When removing the first extremum data, the value of the calculated parameters undergoes high amplitude variation, and the more outliers were removed, the more stable the parameters became. It was decided to cut the dataset where the variation of the dataset has the least influence on the result.
At the moment, the described algorithm is limited to circular pipes, but it can be adapted to other pipe geometry. Another interesting axis of research would be to extend this work to transitional and turbulent flow.
The calculations proposed in Section 3 do not require heavy computations and can be implemented on a single board computer (SBC) or a programmable logic controller (PLC). It is therefore possible to build an automatic non-Newtonian rheometer based on differential pressure measurements in a pipe that can provide continuous estimations of the Herschel-Bulkley parameters with an uncertainty that is of the same order of magnitude as those obtained with a Scientific rheometer.

5. Conclusions

It has been shown that:
  • It is possible to estimate the Herschel-Bulkley rheological behavior parameters utilizing differential pressure measurements along a circular pipe made at different volumetric flow rates. The advantage of using pressure gradients to obtain information about the viscous properties of non-Newtonian fluid is that it does not put any specific requirements on the transparency of the fluid nor the possible negative side effects of diffractions when attempting to measure the fluid velocity field when the fluid contains large proportions of solid particles.
  • The method described by Mullineux [14] to calibrate the Herschel-Bulkley parameters based on rheometer measurements, i.e., a series of pairs of shear rate and shear stress at the wall could also be transposed to the context of calibrating the parameters of Herschel-Bulkley fluid utilizing a series of pairs of volumetric flow rate and pressure gradients. The method has the advantage of being also precise at a low shear rate which is not the case of the method based on a logarithm development of the difference of the shear and yield stresses [12].
  • The obtained precision of the calibrated parameters is of the same order of magnitude as the one obtained with a scientific rheometer.
  • The calibration method based on the method from Mullineux is simple enough to be implemented on real-time computer systems such as single-board computers or programmable logic controllers, therefore allowing the possibility to devise real-time equipment that can measure continuously the rheological behavior of non-Newtonian fluids that follow the Herschel-Bulkley rheological behavior.
The possibility to measure continuously the rheological behavior of Herschel-Bulkley fluid in a setup where settling of particles is unlikely and which can be operated completely automatically and with a precision of the same order of magnitude as the one from a scientific rheometer, is important for industries that utilize non-Newtonian fluids with rheological characteristics that are likely to change with time.

Author Contributions

Conceptualization, E.M. and E.C.; methodology, E.M.; software, E.M.; validation, E.M.; formal analysis, E.M.; investigation, E.M.; resources, E.M.; data curation, E.M.; writing—original draft preparation, E.M.; writing—review and editing, E.C.; visualization, E.M.; supervision, E.C.; project administration, E.C.; funding acquisition, E.C. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Norwegian Research Council, AkerBP, Equinor, Repsol Norge AS, TOTAL, Vår Energi, and Petrobras through the Demo2000 funding program.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable

Data Availability Statement

Not applicable.

Acknowledgments

The authors would like to thank Amare Leulseged from the Norwegian Research Centre for his valuable contribution in gathering the experimental measurements that have been used in this work.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

Nomenclature

d P d L pressure gradient [ML−2T−2] (Pa/m)
K consistency index [ML−1Tn−2] (Pa.sn)
K M consistency index obtained with the Mullineux method based on the flow-loop measurements [ML−1Tn−2] (Pa.sn)
K L M consistency index obtained with the Levenberg-Marquart optimization [ML−1Tn−2] (Pa.sn)
K R h e o consistency index obtained with the Mullineux method based on the scientific rheometer measurements [ML−1Tn−2] (Pa.sn)
m number of measurements
n flow behavior index (dimensionless)
n M flow behavior index obtained from the Mullineux method based on the flow-loop measurements (dimensionless)
n L M flow behavior index obtained with the Levenberg-Marquart optimization (dimensionless)
n R h e o flow behavior index obtained with the Mullineux method on the scientific rheometer measurements (dimensionless)
p number of bins in the non-flowing condition histogram
Q volumetric flowrate [L3T−1](m3/s)
R pipe radius [L](m)
S least square error [M2L−2T−4] (Pa2)
Greek letters:
γ ˙ shear rate [T−1](1/s)
γ ˙ i measured shear rate [T−1](1/s)
γ ˙ N , w wall shear for a Newtonian fluid [T−1](1/s)
ε shear stress span for a bin in the non-flowing condition histogram [ML−1T−2] (Pa)
τ shear stress [ML−1T−2] (Pa)
τ i measured shear stress [ML−1T−2] (Pa)
τ γ yield stress [ML−1T−2] (Pa)
τ γ ,     M yield stress obtained with the Mullineux method based on the flow-loop measurements [ML−1T−2] (Pa)
τ γ ,   L M yield stress obtained with the Levenberg-Marquart optimization [ML−1T−2] (Pa)
τ γ ,   r h e o yield stress obtained from scientific rheometer measurements [ML−1T−2] (Pa)

References

  1. Heinisch, M.; Abdallah, A.; Dufour, I.; Jakoby, B. Resonant Steel Tuning Forks for Precise Inline Viscosity and Mass Density Measurements in Harsh Environments. Procedia Eng. 2014, 87, 1139–1142. [Google Scholar] [CrossRef] [Green Version]
  2. Ovarlez, G.; Mahaut, F.; Bertrand, F.; Chateau, X. Flows and heterogeneities with a vane tool: Magnetic resonance imaging measurements. J. Rheol. 2011, 55, 197–223. [Google Scholar] [CrossRef] [Green Version]
  3. Molki, A.; Khezzar, L.; Goharzadeh, A. Measurement of fluid velocity development in laminar pipe flow using laser Doppler velocimetry. Eur. J. Phys. 2013, 34, 1127–1134. [Google Scholar] [CrossRef]
  4. Wiklund, J.; Shahram, I.; Stading, M. Methodology for in-line rheology by ultrasound Doppler velocity profiling and pressure difference techniques. Chem. Eng. Sci. 2007, 62, 4277–4293. [Google Scholar] [CrossRef]
  5. Machin, T.D.; Wei, H.-Y.; Greenwood, R.W.; Simmons, M.J.H. In-pipe rheology and mixing characterisation using electrical resistance sensing. Chem. Eng. Sci. 2018, 187, 327–341. [Google Scholar] [CrossRef]
  6. Carlsen, L.A.; Rolland, N.L.; Nygaard, G.; Time, R. Simultaneous Continuous Monitoring of the Drilling-Fluid Friction Factor and Density. SPE Drill. Completion 2013, 28, 34–44. [Google Scholar] [CrossRef]
  7. Herschel, W.H.; Bulkley, R. Konsistenzmessungen von Gummi-Benzollösungen. Kolloid Z. 1926, 39, 291–300. [Google Scholar] [CrossRef]
  8. Robertson, R.E.; Stiff, H.A., Jr. An Improved Mathematical Model for Relating Shear Stress to Shear Rate in Drilling Fluids and Cement Slurries. Soc. Pet. Eng. J. 1976, 16, 31–36. [Google Scholar] [CrossRef]
  9. Heinz, W. The Casson flow equation: Its validity for suspension of paints. Mater. PrüFung 1959, 1, 311. [Google Scholar]
  10. Sochi, T. Analytical solutions for the flow of Carreau and Cross fluids in circular pipes and thin slits. Rheol. Acta 2015, 54, 745–756. [Google Scholar] [CrossRef] [Green Version]
  11. Quemada, D. Rheological modelling of complex fluids. I. The concept of effective volume fraction revisited. Eur. Phys. J. AP 1998, 1, 119–127. [Google Scholar] [CrossRef]
  12. Kelessidis, V.C.; Maglione, R.; Tsamantaki, C.; Aspirtakis, Y. Optimal determination of rheological parameters for Herschel–Bulkley drilling fluids and impact on pressure drop, velocity profiles and penetration rates during drilling. J. Pet. Sci. Eng. 2006, 53, 203–224. [Google Scholar] [CrossRef]
  13. Turian, R.M.; Ma, T.W.; Hsu, F.L.G.; Sung, D.J. Characterization, settling, and rheology of concentrated fine particulate mineral slurries. Powder Technol. 1997, 93, 219–233. [Google Scholar] [CrossRef]
  14. Mullineux, G. Non-linear least squares fitting of coefficients in the Herschel–Bulkley model. Appl. Math. Model. 2008, 32, 2538–2551. [Google Scholar] [CrossRef]
  15. Rooki, R.; Ardejani, F.D.; Moradzadeh, A.; Mirzaei, H.; Kelessidis, V.; Maglione, R.; Norouzi, M. Optimal determination of rheological parameters for herschel-bulkley drilling fluids using genetic algorithms (GAs). Korea-Aust. Rheol. J. 2012, 24, 163–170. [Google Scholar] [CrossRef]
  16. Cayeux, E.; Leulseged, A. The Effect of Thixotropy on Pressure Losses in a Pipe. Energies 2020, 13, 6165. [Google Scholar] [CrossRef]
  17. Skadsem, H.J.; Saasen, A. Concentric cylinder viscometer flows of Herschel-Bulkley fluids. Appl. Rheol. 2019, 29, 173–181. [Google Scholar] [CrossRef] [Green Version]
  18. Lac, É.; Parry, A. Non-Newtonian end-effects in standard oilfield rheometers. J. Rheol. 2017, 61, 833–843. [Google Scholar] [CrossRef]
  19. Cayeux, E.; Leulseged, A. Characterization of the rheological behavior of drilling fluids. In Proceedings of the ASME 2020 39th International Conference on Ocean, Offshore and Arctic Engineering, Fort Lauderdale, FL, USA, 28 June–3 July 2020. [Google Scholar] [CrossRef]
  20. Vajargah, A.K.; van Oort, E. Determination of drilling fluid rheology under downhole conditions by using real-time distributed pressure data. J. Nat. Gas Sci. Eng. 2015, 24, 400–411. [Google Scholar] [CrossRef]
  21. Vajargah, A.K.; Sullivan, G.; van Oort, E. Automated fluid rheology and ECD management. In Proceedings of the SPE Deepwater Drilling and Completions Conference, Galveston, TX, USA, 14–15 September 2016. [Google Scholar] [CrossRef]
  22. Gul, S.; Erge, O.; van Oort, E. Frictional pressure losses of non-Newtonian fluids in helical pipes: Applications for automated rheology measurements. J. Nat. Gas Sci. Eng. 2020, 73, 103042. [Google Scholar] [CrossRef]
  23. Magalhaes, S.; Scheid, C.M.; Calcada, L.A.; Folsta, M.; Martins, A.L.; Marques deSa, C.H. Development of on-line sensors for automated measurement of drilling fluid properties. In Proceedings of the IADC/SPE Drilling Conference and Exhibition, Fort Worth, TX, USA, 4–6 March 2014. [Google Scholar] [CrossRef]
  24. Macosco, C.W. Rheology: Principles, Measurements and Applications; Wiley-VCH: Weinheim, Germany, 1994. [Google Scholar]
  25. Morrison, F.A. Understanding Rheology; Oxford University Press: New York, NY, USA, 2001. [Google Scholar]
  26. Levenberg, K. A Method for the Solution of Certain Non-Linear Problems in Least Squares. Q. Appl. Math. 1944, 2, 164–168. [Google Scholar] [CrossRef] [Green Version]
  27. Divoux, T.; Barentin, C.; Manneville, S. Stress overshoot in a simple yield stress fluid: An extensive study combining rheology and velocimetry. Soft Matter 2011, 7, 9335–9349. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Schematic description of the laboratory flow loop (courtesy of Cayeux and Leulseged in [16]).
Figure 1. Schematic description of the laboratory flow loop (courtesy of Cayeux and Leulseged in [16]).
Fluids 06 00157 g001
Figure 2. Comparison of the effect of utilizing the Newtonian-fluid shear rate conversion instead of the non-Newtonian shear rate conversion of the rheometer speed for a rheogram obtained with an aqueous solution of Carbopol 980 having a density of 997 kg/m3 at 25 °C.
Figure 2. Comparison of the effect of utilizing the Newtonian-fluid shear rate conversion instead of the non-Newtonian shear rate conversion of the rheometer speed for a rheogram obtained with an aqueous solution of Carbopol 980 having a density of 997 kg/m3 at 25 °C.
Fluids 06 00157 g002
Figure 3. Graph of ln γ ˙ N , w as a function of ln τ w , shown in red dots where γ ˙ N , w is expressed in reciprocal seconds and τ w in Pascals. The blue line is a linear regression of the 25% highest τ w data points τ w 11.2   P a . The slope is 1.67, therefore the estimated n = 0.60   , which is an acceptable initial approximation considering that the value obtained with the scientific rheometer was 0.64 .
Figure 3. Graph of ln γ ˙ N , w as a function of ln τ w , shown in red dots where γ ˙ N , w is expressed in reciprocal seconds and τ w in Pascals. The blue line is a linear regression of the 25% highest τ w data points τ w 11.2   P a . The slope is 1.67, therefore the estimated n = 0.60   , which is an acceptable initial approximation considering that the value obtained with the scientific rheometer was 0.64 .
Fluids 06 00157 g003
Figure 4. Estimations of τ γ looking at pressure losses at a low flow rate. (a) Pressure losses at a low flow rate are not steady since it depends on the flow history of the fluid. (b) This histogram shows the percentage of occurrences of non-flowing state at a given wall shear stress τ w up to the stress overshoot value.
Figure 4. Estimations of τ γ looking at pressure losses at a low flow rate. (a) Pressure losses at a low flow rate are not steady since it depends on the flow history of the fluid. (b) This histogram shows the percentage of occurrences of non-flowing state at a given wall shear stress τ w up to the stress overshoot value.
Fluids 06 00157 g004
Figure 5. Estimations of K using 10% of the measured Q ,   d P d L pair data. The final estimation is obtained by taking the median of the calculations. Here, K = 0.219   Pa . s n , which is an acceptable initial value considering that the consistency index measured with the scientific rheometer was 0.272   Pa . s n .
Figure 5. Estimations of K using 10% of the measured Q ,   d P d L pair data. The final estimation is obtained by taking the median of the calculations. Here, K = 0.219   Pa . s n , which is an acceptable initial value considering that the consistency index measured with the scientific rheometer was 0.272   Pa . s n .
Fluids 06 00157 g005
Figure 6. Measured and calculated pressure losses as a function of the flow rate. Curves named Rheometer, Levenberg-Marquart, and Mullineux are the result of a calculation using the H-B parameters estimated by the technique referring to their name and inverting Equation (8).
Figure 6. Measured and calculated pressure losses as a function of the flow rate. Curves named Rheometer, Levenberg-Marquart, and Mullineux are the result of a calculation using the H-B parameters estimated by the technique referring to their name and inverting Equation (8).
Fluids 06 00157 g006
Figure 7. Estimation of d ln γ ˙ N , w d ln τ w using Equation (12) and the method corresponding to the legend for the first three curves and using a polynomial fit of ln γ ˙ N , w ln τ w in the order mentioned in the legend. γ ˙ N , w is expressed in reciprocal seconds and τ w in Pascals.
Figure 7. Estimation of d ln γ ˙ N , w d ln τ w using Equation (12) and the method corresponding to the legend for the first three curves and using a polynomial fit of ln γ ˙ N , w ln τ w in the order mentioned in the legend. γ ˙ N , w is expressed in reciprocal seconds and τ w in Pascals.
Fluids 06 00157 g007
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Magnon, E.; Cayeux, E. Precise Method to Estimate the Herschel-Bulkley Parameters from Pipe Rheometer Measurements. Fluids 2021, 6, 157. https://doi.org/10.3390/fluids6040157

AMA Style

Magnon E, Cayeux E. Precise Method to Estimate the Herschel-Bulkley Parameters from Pipe Rheometer Measurements. Fluids. 2021; 6(4):157. https://doi.org/10.3390/fluids6040157

Chicago/Turabian Style

Magnon, Elie, and Eric Cayeux. 2021. "Precise Method to Estimate the Herschel-Bulkley Parameters from Pipe Rheometer Measurements" Fluids 6, no. 4: 157. https://doi.org/10.3390/fluids6040157

APA Style

Magnon, E., & Cayeux, E. (2021). Precise Method to Estimate the Herschel-Bulkley Parameters from Pipe Rheometer Measurements. Fluids, 6(4), 157. https://doi.org/10.3390/fluids6040157

Article Metrics

Back to TopTop