Next Article in Journal
Biodiesel Production through Acid Catalyst In Situ Reactive Extraction of Chlorella vulgaris Foamate
Previous Article in Journal
Analysis of Impact of Design Solutions of an Electric Machine with Permanent Magnets for Bearing Voltages with Inverter Power Supply
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Theory of the Vom Berg Rheological Model and Its Use in Cloud-Native Application

by
Rafał Wiśniowski
and
Grzegorz Orłowicz
*
Faculty of Drilling Oil and Gas, Department of Drilling and Geoengineering, AGH University of Science and Technology, 30 Mickiewicza Av., 30-059 Kraków, Poland
*
Author to whom correspondence should be addressed.
Submission received: 23 May 2022 / Revised: 14 June 2022 / Accepted: 17 June 2022 / Published: 20 June 2022

Abstract

:
Various technological fluids, such as drilling muds, drill-in fluids, fracturing fluids, spacers, washes and cement slurries are used in the wellbore drilling process. The fundamental issue, which needs to be addressed in order to become acquainted with the phenomena occurring during fluids flow through a circulatory system, is to establish mutual dependencies between a stream of fluid being pumped and flow resistances. The awareness of these dependencies enables the optimisation of hydraulic parameters in order to minimise costs and maximise drilling works safety. This article presents rheological models of drilling fluids and proposes the application of a new rheological model, not used in the drilling industry so far, namely the Vom Berg model. The model has been presented in other publications; however, there is an unsolved and unpublished problem of determining the effect of rheological parameters of the model on the value of resistance to laminar and turbulent flow. In this article, algorithms and Cloud-Native application enabling numerical determination of rheological properties of the Vom Berg fluid are presented. What is more, an algorithm for calculating pressure losses during the laminar flow of fluid in a pipe is provided. Taking an example from the industry, a practical application of the proposed calculation methodology is presented.

1. Introduction

A number of boreholes are made every year for the needs of the energy and environmental engineering sectors. When drilling oil wells and geothermal boreholes, it is possible to acquire fluid energy raw materials and store energy and, when using trenchless technologies (horizontal directional drilling, microtunneling direct pipe), the underground transport of utilities to long distances is ensured. During making boreholes, various technological fluids are used, such as drilling muds, drill-in fluids, fracturing fluids, spacers, washes and cement slurries. The knowledge of physicochemical phenomena occurring during drilling fluid flow enables rational designing of the technique and technology in order to guarantee maximum works safety and minimum investment expenditures.
The fundamental issue, which needs to be solved in order to learn about the phenomena occurring during drilling fluids flow through a circulatory system, is to establish mutual dependencies between a stream of fluid being pumped and flow resistances occurring when it is pumped.
Flow resistances in any circulatory system component depend on the following:
  • The type of flowing fluid;
  • Geometry and physical properties of the element in which flow takes place;
  • Flow regime.
The most important factors affecting the amount of pressure losses include physical properties and rheological parameters of the flowing fluid.
This article presents rheological models of drilling fluids and proposes the application of a new rheological model, not used in the drilling industry so far, namely the Vom Berg model. The model has been presented in other publications [1,2]; however, there is an unsolved and unpublished problem of determining the effect of rheological parameters of the model on the value of resistance to laminar and turbulent flow. The aim of this study was to present in detail the laminar flow of fluid described by the Vom Berg model.

2. Division of Drilling Fluids According to Rheological Properties

In most cases, fluids used in the drilling industry are Newtonian fluids, viscoplastic (Bingham plastics) and pseudoplastic (Ostwald de Waele).
An API RP 13D (American Petroleum Institute Recommended Practice 13D) standard, containing recommendations concerning rheological studies and hydraulic calculations of drilling fluids, recommends using one of the three models: Bingham, Ostwald de Waele (power law) or Herschel–Bulkley (yield power law) [3,4,5,6]. The form of models are presented in Equations (1)–(3):
  • Bingham model:
τ = τ y + η ( dv dr )
  • Ostwald de Waele model (power law):
τ = k ( dv dr ) n
  • Herschel–Bulkley model (yield power law):
τ = τ y + k ( dv dr ) n
The application of chemical agents enabling the modification of physical properties and rheological parameters of drilling muds and cement slurries causes those fluids to potentially demonstrate both the features of viscoelastic fluids and pseudoplastic fluids.
Thus, for this kind of fluid, other rheological models should be considered. One of the possible ones is the Vom Berg model in the form of [1,2,7]:
τ = τ y + D sinh 1 ( dv dr G )

3. Laminar Flow of Reostable Fluids through a Pipe

Finding dependencies between a stream of flowing fluid and flow resistances is possible after considering fluid flow in a straight horizontal pipe. From the stream of fluid, we can distinguish (Figure 1) a part of fluid volume in the form of a cylinder with the radius r and the length of ΔL.
In order to obtain the condition of balance of forces applied at separate mass of fluids, the sum of projections of external forces on the direction of the pipe axis has to equal zero.
Assuming that flow occurs to the right p1 > p2, the forces distribution is as follows:
  • Pressure p1 = p acts on cross Section 1 with an area of πr2,
  • Pressure p2 acts on cross Section 2 with an area of πr2, that pressure by number equals to pressure p reduced by pressure losses at the length ΔL: p2 = ( p dp dl Δ l ) ,
  • At the cylinder side with an area of 2πrΔl friction force occurs during flow, which is caused by shear stresses τ coming from the fluid.
The equilibrium equation of forces can be thus written in the form of
p π r 2 2 π r Δ l τ ( p dp dl Δ l ) π r 2 = 0
After transformation, we receive:
dp dl = 2 r τ
The above equation is a basic dependency enabling the establishment of correlations between a stream of fluid flowing through a pipe (regardless of its rheological properties) and flow resistances in the laminar regime.
The sought relationships for the Vom Berg model can be obtained by inserting to Equation (6) rheological Equation (4):
dp dl = 2 r ( τ y + D sinh 1 ( dv dr G ) )
From Equation (7), one can determine a shear rate gradient   γ ˙ = ( dv dr )   by means of the following relationship:
( dv dr ) = G sinh ( r 2 D dp dl τ y D )
After variables separation in Equation (8) and integration, the following is obtained:
v = G 2 D dp dl cosh ( r 2 D dp dl τ y D ) + C 1
The constant C1 is determined from the boundary conditions. It is assumed that at the pipe wall, for the value r = R, the fluid stream flow rate is v = 0.
C 1 = G 2 D dp dl cosh ( R 2 D dp dl τ y D )
By inserting the value C1 to Equation (9), we obtain the distribution of the Vom Berg fluid flow rate along the radius of the pipe in which flow takes place:
v = G 2 D dp dl [ cosh ( R 2 D dp dl τ y D ) cosh ( r 2 D dp dl τ y D ) ]
The minimum speed value is zero and occurs at the pipe wall. The maximum fluid flow rate is determined from the following condition: v ( r 0 ) = 0 v ( r 0 ) < 0 . This rate occurs at a distance r 0 = 2 τ y dl dp from the pipe centre and is determined on the basis of the following formula:
V max = G 2 D dp dl [ cosh ( R 2 D dp dl τ y D ) 1 ]
The total stream of the Vom Berg fluid is the sum of the fluid nucleus stream Q0 (a cylinder with the radius r0 moving with constant velocity v0 = vmax) and stream Q1 flowing around a hollow cylinder with the internal radius r0 and the external radius R.
Q = 0 R v ( r ) dF = 0 r 0 v 0 dF + r 0 R v ( r ) dF
Assuming unit cross-section dF = 2πrdr, the total stream of the Vom Berg fluid flowing through a pipe is determined on the basis of the following dependence:
Q = Q 0 + Q 1 = v 0 π r 0 2 + r 0 R v ( r ) dF
Q = v 0 π r 0 2 + r 0 R G 2 D dp dl [ cosh ( R 2 D dp dl τ y D ) cosh ( r 2 D dp dl τ y D ) ] 2 π rdr
By determining the integral in Equation (15) and taking into account the relationship r 0 = 2 τ y dl dp , we obtain:
Q = 16 π GD 3 ( dp dl ) 3 [ ( R 2 8 D 2 ( dp dl ) 2 + 1 ) cosh ( R 2 D dp dl τ y D ) R 2 D dp dl sinh ( R 2 D dp dl τ y D ) ( τ y ) 2 2 D 2 1 ]
The dependence of unit pressure losses ( dp dl ) on the flow stream for laminar flow (Q) of the Vom Berg fluid is determined by means of the implicit function theorem (16).
In order to solve Equation (16), the following function is created:
g ( x ) = 16 π GD 3 ( x ) 3 [ ( R 2 8 D 2 ( x ) 2 + 1 ) cosh ( R 2 D x τ y D ) R 2 D x sinh ( R 2 D x τ y D ) ( τ y ) 2 2 D 2 1 ] Q
and next, its zero is designated. As can be seen in Figure 2, the g(x) function has two zeros.
The condition of physical interpretation of solving Equation (17) is the value of its zero, ( x = dp dl ) , ensuring a positive shear rate value:
( dv dr ) = G sinh ( r 2 D dp dl τ y D ) > 0
For G > 0, a physical condition is obtained, which has to meet the root of Equation (17):
x > 2 τ y R
It is proposed to solve the relationship g(x) = 0 by means of numerical methods: the bisection method, the Newton–Raphson method, the regula falsi method, combined methods, etc. [8,9,10]. To this end, the Department of Drilling and Geoengineering, at the Faculty of Drilling, Oil and Gas, AGH University of Science and Technology in Kraków, developed a numerical program, using the bisection method. The algorithm of solving numerical Equation (17) is presented in Figure 3.
The boundaries of the interval of the solution being sought ( x sr = dp dl ) are suggested as follows: a = 2 τ y R , b = 106, the accuracy of its estimation epsilon = 0.00001. Total pressure losses are determined from the relationship ( p   = x sr · L ) .

4. Determination of Rheological Parameters and Resistances of Drilling Fluid Laminar Flow

The selection of the accurate rheological model has an impact on the accuracy of the calculation results received. Drilling fluid rheological properties are measured by means of rotational viscometers of the following types: Fann, Chan, Brookfield, Haake, and Ofite [11,12,13,14,15,16]. A fluid test consists in determining shear stresses relationships (τ) in a shear rate function   (   γ ˙ = ( dv dr ) ) . Drilling fluid (Bingham, pseudoplastic and Herschel–Bulkley) the rheological parameters are specified on the basis of standards [3,4] and the relationships which are developed and not accounted for in the standards [6,17].
In order to determine rheological parameters for a three-parameter Vom Berg fluid model, it is necessary to know the shear stress values for three shear rates ( γ ˙ LOW ,   γ ˙ MID ,   γ ˙ TOP ) .
The basic problem is to determine a real value of fluid shear rate (   γ ˙ MID ) . Apart from its rheological parameters, also flow conditions, such as a stream of flow volume and geometry of an element in which flow takes place, have an impact on the real value of shear rates of flowing drilling fluid. Precise values of real shear rates of fluid can be determined from Formula (6) only after establishing fluid rheological parameters.
For the laminar flow of a Newtonian fluid through a pipe, the shear rate is determined as follows [18]:
  γ ˙ = ( dv dr ) = 8 v d
This value is proposed to be assumed as the starting point for establishing the actual scope of shear rates for the Vom Berg fluid.
Therefore, it is suggested, before starting calculations, to determine the shear rate value (   γ ˙ )   from Formula (20). Next, for a given viscometer, by means of which dependences between shear stresses and shear rate are determined, the most approximated value of   γ ˙ MID 8 v d should be adopted. The value will define the shear rate interval, which should be considered in calculations of drilling fluid rheological parameters   γ ˙ MID ∈(   γ ˙ LOW   ;   γ ˙ TOP ). The value   γ ˙ LOW is the closest shear rate value smaller than   γ ˙ MID recorded on a viscometer. The value   γ ˙ TOP is the closest shear rate value larger than   γ ˙ MID recorded on a viscometer. For the assumed shear rate values (   γ ˙ LOW ,   γ ˙ MID ,   γ ˙ TOP ) and shear stresses corresponding to them ( τ LOW , τ MID , τ TOP ) , one may determine rheological parameters of fluid (τy, D and G) described with the Vom Berg model.
For the Vom Berg fluid model, rheological parameters: τy, D and G are determined by solving the following system of equations:
{ τ TOP = τ y + D sinh 1 (   γ ˙ TOP G ) τ MID = τ y + D sinh 1 (   γ ˙ MID G ) τ LOW = τ y + D sinh 1 (   γ ˙ LOW G )
The value of the G parameter is obtained from the numerical solution of the equation:
sinh 1 (   γ ˙ TOP G )   sin h 1 (   γ ˙ MID G ) sinh 1 (   γ ˙ TOP G ) sinh 1 (   γ ˙ LOW G ) = τ TOP τ MID τ TOP τ LOW
In order to solve Equation (22), it is assumed that x = G and the g(x) function is created:
g ( x ) = sinh 1 (   γ ˙ TOP x )   sin h 1 (   γ ˙ MID x ) sinh 1 (   γ ˙ TOP x ) sinh 1 (   γ ˙ LOW x ) τ TOP τ MID τ TOP τ LOW
and next its zero is designated.
The solution of the g(x) = 0 relationship is proposed by means of the numerical methods: the bisection method, the Newton–Raphson method, the regula falsi method, combined methods, etc. [8,9,19]. The algorithm of solving Equation (23) by means of the secant method is presented in Figure 4.
When analysing the existence of the solution of the Equation (23), we should analyse the g(x) function flow. Assuming:
g ( x ) = f ( x ) W
where:
f ( x ) = sinh 1 (   γ ˙ TOP x )   sin h 1 (   γ ˙ MID x ) sinh 1 (   γ ˙ TOP x ) sinh 1 (   γ ˙ LOW x ) ,   and   W = τ TOP τ MID τ MID τ LOW
one should state that the g(x) function is the f(x) function shifted in relation to the axis of dependent variables by the value W. Analysing the f(x) function flow for positive values x, it can be stated that it is a continuous, increasing and differentiable function. It also does not intersect the OX axis.
The smallest value is obtained for x→0: f ( x ) MIN = lim x 0 f ( x ) .
Whereas the highest value is obtained for x→∞: f ( x ) MAX = lim x f ( x ) .
The g(x) function will intersect the OX axis when W f ( x ) MIN and when W f ( x ) MAX .
Figure 5 shows the g(x) function flow depending on the W parameter value.
The g(x) function will have a root when the following condition is met:
lim x 0 f ( x ) < τ TOP τ MID τ TOP τ LOW < lim x f ( x )
Taking into account
f ( x ) MIN = lim x 0 f ( x ) = lim x 0 sinh 1 ( γ ˙ TOP x )     sinh 1 ( γ ˙ MID x ) sinh 1 ( γ ˙ TOP x )     sinh 1 ( γ ˙ 200 x ) = lim x 0 ln ( γ ˙ TOP + ( γ ˙ TOP ) 2 + x 2 γ ˙ MID + ( γ ˙ MID ) 2 + x 2 ) ln ( γ ˙ TOP + ( γ ˙ TOP ) 2 + x 2 γ ˙ LOW + ( γ ˙ LOW ) 2 + x 2 ) = ln ( γ ˙ TOP γ ˙ MID ) ln ( γ ˙ TOP γ ˙ LOW )
f ( x ) MAX = lim x f ( x ) = lim x ln ( γ ˙ TOP + ( γ ˙ TOP ) 2 + x 2 γ ˙ MID + ( γ ˙ MID ) 2 + x 2 ) ln ( γ ˙ TOP + ( γ ˙ TOP ) 2 + x 2 γ ˙ LOW + ( γ ˙ LOW ) 2 + x 2 ) = lim x ( ln ( γ ˙ TOP + ( γ ˙ TOP ) 2 + x 2 γ ˙ MID + ( γ ˙ MID ) 2 + x 2 ) ) ( ln ( γ ˙ TOP + ( γ ˙ TOP ) 2 + x 2 γ ˙ LOW + ( γ ˙ LOW ) 2 + x 2 ) )
= lim x ( γ ˙ LOW x + (   γ ˙ LOW x ) 2 + 1 ) (   γ ˙ MID x + (   γ ˙ MID x ) 2 + 1 ) ( (   γ ˙ LOW x ) 2 + 1 ) ( (   γ ˙ MID x ) 2 + 1 )
γ ˙ TOP ( ( γ ˙ TOP x ) 2 + 1 )   +   ( γ ˙ TOP x ) 2   γ ˙ MID ( ( γ ˙ MID x ) 2 + 1 )     ( γ ˙ MID x ) 2 γ ˙ TOP ( ( γ ˙ TOP x ) 2 + 1 )   +   ( γ ˙ TOP x ) 2   γ ˙ LOW ( ( γ ˙ LOW x ) 2 + 1 )     ( γ ˙ LOW x ) 2 = γ ˙ TOP γ ˙ MID γ ˙ TOP γ ˙ LOW
the condition of occurrence of the root of the g(x) function is obtained:
ln (   γ ˙ TOP   γ ˙ MID ) ln (   γ ˙ TOP γ LOW ) < τ TOP τ MID τ TOP τ LOW <   γ ˙ TOP   γ ˙ MID   γ ˙ TOP   γ ˙ LOW
Knowing the value of the G parameter (G = x), the other rheological parameters should be determined from the following relationship:
D = τ TOP   τ MID sinh 1 (   γ ˙ TOP G ) sinh 1 (   γ ˙ MID G )
τ y = τ TOP D sinh 1 (   γ ˙ TOP G )
The calculated rheological parameters of the Vom Berg model (τy, D and G) enable to determine unit flow resistances ( dp dl )   of a fluid stream (Q) flowing inside a pipe with the radius (R).
The following algorithm is suggested:
1.
Measurement, by means of a viscometer, shear stress dependencies for different drilling fluid shear rate values.
2.
Estimation of the approximated shear rate value (   γ ˙ START ) (occurring when fluid stream (Q) flows inside a pipe with the radius (R). Formula (20).
3.
Determination of the most approximated shear rate value (   γ ˙ MID   γ ˙ START ) , taking into account shear stress dependencies on shear rate, obtained from a given viscometer; establishment of the value (   γ ˙ LOW ,   γ ˙ TOP ) . Determination of the measured values ( τ LOW , τ MID , τ TOP ) .
4.
Calculation of rheological parameters of the fluid described by the Vom Berg model. Numerical solution of the Equation (22) and use of relationships (27) and (28).
5.
Calculation of unit flow resistances   ( dp dl ) . Numerical solution of Equation (17).
6.
Calculation of real shear rate (   γ ˙ START ) , characterising fluid flow, Equation (8).
7.
Evaluation of the accuracy of the assumptions made in Item 3 of the algorithm. Checking of the condition   γ ˙ START (   γ ˙ LOW ,   γ ˙ TOP ) . If this condition is met, calculations of unit flow resistances are considered completed.
If the condition   γ ˙ START (   γ ˙ LOW ,   γ ˙ TOP )   is not met, calculations are repeated beginning with Item 3, assuming the value of the real shear rate (   γ ˙ START ) determined in Item 6.

5. Application of Derived Equations in Cloud-Hosted Computer-Aided Design (CAD) Tool

The main purpose of the presented work was to derive and present mathematical formulas of laminar flow of the Vom Berg model. In order to present the formulas application possibilities to larger audience, the authors created a new computer-aided design tool. The tool is called the HDD Designer, and it is available on the Internet via web browsers. The address of the tools is www.HDDDesigner.com, accessed on 10 June 2022.
Rheology Sim is the first finished module in the HDD Designer, and it includes following features:
  • Rheological parameters calculation.
  • Rheological model specific parameters calculation for Newton, Bingham, Casson, Ostwald De Waele, Herschel–Bulkley, Vom Berg and Eyring models.
  • Regression parameters calculation.
  • Graphical and tabular result presentation.
As the next step, the authors plan to develop another module for calculations of flow resistances and pressure losses during laminar flow of fluid in a drill pipe and annulus for various rheological models.
The graphical user interface of the tool is presented on the Figure 6. The HDD Designer has been built with use of Angular (Graphical User Interface) and NET Core (application programming interface), it is hosted on Microsoft Azure Cloud. As it uses the most recent available technologies and cloud hosting, it is an easily expandable, robust, scalable and high availability application.

6. A practical Example of the Proposed Methodology Application

The Vom Berg model can be applied to various technological fluids, such as drilling muds, drill-in fluids, fracturing fluids, spacers, washes and cement slurries. In order to apply the Vom Berg model in practice, the unit flow resistances ( dp dl ) were calculated for the selected sealing slurry.
In order to exemplify the developed relationships, flow resistances of cement slurry pumped to casing with the external diameter of 6 5/8″ (the radius of the pipe interior R = 0.0736 m), with a stream of flow volume Q = 0.03 m3/s were calculated.
By means of a Fann 35 rotational viscometer (in the arrangement R1-B1 and a spring F1), rheological properties of cement slurry to be used in a borehole [15] were measured. The measurement results of a torsion angle of the outer cylinder (Φ), read for the specified rotational speeds (n) and calculated shear stress values (τ) and shear rate   (   γ ˙ ) , occurring during measurements, are provided in Table 1. A graphical comparison of popular rheological models’ fitness to measured values is presented on Figure 7. Table 2 shows the tabular comparisons of statistical parameters from popular rheological models. Based on both comparisons it is visible that the Vom Berg model achieved the highest fitness for the analysed example.
Next, on the basis of the developed methodology, the slurry rheological parameters and unit flow resistances were calculated, assuming the Vom Berg fluid model.
To this end, the approximated shear rate inside 6 5/8″ casing (   γ ˙ START 96   s 1 ) was estimated, using Equation (18). That value was used for estimating the value of   γ ˙ LOW ,   γ ˙ MID and   γ ˙ TOP . For those values, the slurry rheological parameters and unit flow resistances were calculated. The calculation of real shear rate   γ ˙ START = 101.6   s 1 and comparing it with its acceptable range   γ ˙ START (   γ ˙ LOW ,   γ ˙ TOP ) confirmed the accuracy of the assumptions made. The calculation results are spresented in Table 3.

7. Conclusions

In the case of fluids that demonstrate features of both viscoelastic and pseudoplastic fluids, we should consider other rheological models than those used so far in the drilling industry. One of the possible models to be applied is the Vom Berg model.
When describing the laminar flow of Vom Berg fluids through a pipe, we should use the relationships introduced in this article. The above-presented limitations and conditions of obtaining an unambiguous solution of Equation (17) allow for numerical calculation of resistances of the laminar flow of Vom Berg fluids. The proposed block diagram and the developed HDD Designer application should be used in engineering practice.
Rheological parameters of the Vom Berg model should be determined depending on the real shear rate, recorded during drilling fluid flow. To this end, the methodology suggested in this article should be applied.

Author Contributions

Conceptualisation, R.W. and G.O.; Data curation, R.W. and G.O.; Formal analysis, R.W. and G.O.; Funding acquisition, R.W.; Investigation, R.W. and G.O.; Methodology, R.W. and G.O.; Project administration, R.W.; Resources, R.W.; Software, G.O.; Supervision, R.W.; Validation, R.W. and G.O.; Visualisation, R.W. and G.O.; Writing—original draft, R.W. and G.O.; Writing—review and editing, G.O. All authors have read and agreed to the published version of the manuscript.

Funding

Research works realized at AGH-UST in Kraków—Faculty of Drilling, Oil and Gas. This investigation was supported under grant WWNiG AGH.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

a , b interval boundaries in the bisection method or the regula falsi method, [-];
C 1 constant;
d internal pipe diameter, [m];
D rheological parameter in the Vom Berg model, [Pa];
η plastic viscosity, [Pas];
G rheological parameter in the Vom Berg model, [s−1];
γ ˙ shear rate gradient, [s−1];
γ ˙ LOW lower measured value of a shear rate gradient, [s−1];
γ ˙ MID middle measured value of a shear rate gradient, [s−1];
γ ˙ START reference shear rate gradient, [s−1];
γ ˙ TOP upper measured value of a shear rate gradient, [s−1];
k coefficient of drilling fluid consistency, [Pasn];
n exponential index, [-];
p pressure, [Pa];
Q flow rate, [m3/s];
R internal radius, [m];
r distance from the pipe axis, [m];
r 0 distance from the axis of a pipe in which fluid flows with constant velocity v0, [m];
τ shear stress, [Pa];
τ y yield point, [Pa];
τ LOW lower measured value of shear stress, [Pa];
τ MID middle measured value of shear stress, [Pa];
τ TOP upper measured value of shear stress, [Pa];
v flow velocity, [m/s];
v max maximum flow velocity, [m/s];
v 0 constant fluid flow velocity at a distance from 0 to r0 from the pipe axis, [m/s];

References

  1. Vom Berg, W. Influence of specific surface and concentration of solids upon the fow behaviour of cement pastes. Mag. Concr. Res. 1979, 31, 211–216. [Google Scholar] [CrossRef]
  2. Wiśniowski, R.; Skrzypaszek, K.; Małachowski, T. Selection of a Suitable Rheological Model for Drilling Fluid Using Applied Numerical Methods. Energies 2020, 13, 3192. [Google Scholar] [CrossRef]
  3. API RP 13D. Rheology and Hydraulics of Oil-Well Drilling Fluids; Norm of American Petroleum Institute: Washington, DC, USA, 2006. [Google Scholar]
  4. BN-90/1785-01; Płuczka Wiertnicza. Metody Badań w Warunkach Polowych. Wydawnictwa Normalizacjne ”ALFA”: Warszawa, Poland, 1990.
  5. Wiśniowski, R. Metodyka Określania Modelu Reologicznego Cieczy Wiertniczej; Wiertnictwo, Nafta, Gaz R. 18/1; Wydawnictwo AGH: Kraków, Poland, 2001. [Google Scholar]
  6. Wiśniowski, R. Zastosowanie Modelu Herschela-Bulkleya w Hydraulice Płuczek Wiertniczych; Nowoczesne Techniki i Technologie Bezwykopowe, Zeszyt nr 2/2000; HDD Consulting: Kraków, Poland, 2000; pp. 20–28. [Google Scholar]
  7. Benyounes, K. Rheological behavior of cement-based grout with Algerian bentonite. SN Appl. Sci. 2019, 1, 1037. [Google Scholar] [CrossRef] [Green Version]
  8. Ahmad, A.G. Comparative Study of Bisection and Newton-Rhapson Methods of Root-Finding Problems. Int. J. Math. Trends Technol. (IJMTT) 2015, V19, 121–129. [Google Scholar] [CrossRef]
  9. Dalquist, G.; Biorck, A. Numerical Methods; Dower Publication Inc.: Mineula, NY, USA, 2003. [Google Scholar]
  10. Ehiwario, J.C.; Aghamie, S.O. Comparative Study of Bisection, Newton-Raphson and Secant Methods of Root-Finding Problems. IOSR J. Eng. (IOSRJEN) 2014, 4, 1–7. [Google Scholar]
  11. Catalogue of Brookfield Company. Brookfield Viscometers/Rheometers; Catalogue of Brookfield Company: Middleboro, MA, USA, 1999. [Google Scholar]
  12. Thermo Fisher Scientific Inc. Haake Series 1; Materials from Thermo Fisher Scientific; Thermo Fisher Scientific Inc.: Hanover, Germany, 2018. [Google Scholar]
  13. Materials of Rheology Solution Victoria Company. Haake, Material Characterisation & Testing Products; Materials of Rheology Solution Victoria Company: Truganina, Australia, 2019. [Google Scholar]
  14. Hoechst. Hoechst Polymer Fluids; Information Materials from Hoechst; Hoechst: Höchst, Germany, 1989. [Google Scholar]
  15. Fann Instrument Company. Model 35 Instruction Manual, Manual No. 208878; Fann Instrument Company: Houston, TX, USA, 2016. [Google Scholar]
  16. OFI Testing Equipment, Inc. Testing Equipment for Drilling Fluids · Cement · Core Analysis · Laboratory Supplies · Reagents, Ofite Catalog. 2017. Available online: https://hiltoninstruments.com/pdfUploads/OFITE%20Catalogue.pdf (accessed on 10 June 2022).
  17. Klotz, J.A.; Brigham, W.E. To Determine Herschel-Bulkley Coefficients. J. Pet. Technol. 1998, 50, 80–81. [Google Scholar] [CrossRef]
  18. Bourgoyne, A.T.; Milheim, K.K.; Chenevert, M.E.; Young, F.S. Applied Drilling Engineering. In SPE Textbook; SPE: Houston, TX, USA, 1986; Chapter 4.8. [Google Scholar]
  19. Chen, J.; Li, W. An improved exponential regula falsi method with quadratic convergence of both diameter and point for solving nonlinear equations. Appl. Numer. Math. 2007, 57, 80–88. [Google Scholar] [CrossRef]
Figure 1. Distribution of stresses during laminar flow of fluid through a pipe.
Figure 1. Distribution of stresses during laminar flow of fluid through a pipe.
Energies 15 04481 g001
Figure 2. g(x) function graph.
Figure 2. g(x) function graph.
Energies 15 04481 g002
Figure 3. Block diagram for the calculation of root of the g(x) function being in the range [a, b] by means of the bisection method.
Figure 3. Block diagram for the calculation of root of the g(x) function being in the range [a, b] by means of the bisection method.
Energies 15 04481 g003
Figure 4. Block diagram for the calculation of the root of the g(x) function being in the range [a, b] by means of the secant method.
Figure 4. Block diagram for the calculation of the root of the g(x) function being in the range [a, b] by means of the secant method.
Energies 15 04481 g004
Figure 5. g(x) function flow depending on the W value: (a) W = 0, g(x) = f(x); (b) f(x)MIN < W < f(x)MAX; (c) W < f(x)MIN; (d) W > f(x)MAX.
Figure 5. g(x) function flow depending on the W value: (a) W = 0, g(x) = f(x); (b) f(x)MIN < W < f(x)MAX; (c) W < f(x)MIN; (d) W > f(x)MAX.
Energies 15 04481 g005
Figure 6. Graphical user interface (GUI) of the HDD Designer. The HDD Designer is available on www.hdddesigner.com (accessed on 10 June 2022).
Figure 6. Graphical user interface (GUI) of the HDD Designer. The HDD Designer is available on www.hdddesigner.com (accessed on 10 June 2022).
Energies 15 04481 g006
Figure 7. Graphical comparison of rheological models’ fitness to the measurement.
Figure 7. Graphical comparison of rheological models’ fitness to the measurement.
Energies 15 04481 g007
Table 1. Results of laboratory measurements and calculations of shear stresses and shear rates for cement slurry tested on a 12 speed Fann viscometer 35A/SR-12, with the arrangement of cylinders R1-B1 and spring F1 [15].
Table 1. Results of laboratory measurements and calculations of shear stresses and shear rates for cement slurry tested on a 12 speed Fann viscometer 35A/SR-12, with the arrangement of cylinders R1-B1 and spring F1 [15].
Laboratory Measurement
Rotational speed n, [rot/min]0.91.836306090100180200300600
Torsion1224162838488291141237
Angle Φ, [-]
Calculated Values
Shear rate γ ˙ , [s−1]1.533.075.1110.2251.10102.20153.31170.34306.61340.68511.021022.04
Shear stresses τ, [Pa]0.511.021.022.048.1814.3119.4224.5341.9046.5072.05121.11
Table 2. Tabular comparison of rheological models’ fitness to the measurement.
Table 2. Tabular comparison of rheological models’ fitness to the measurement.
ModelNewtonBinghamCassonOstwal De WaeleHerschel BulkleyEyringVom Berg
Pearson
Correlation Coefficient
0.99470.9960.99650.99390.9980.9970.9985
Fisher Snedecor
Coefficient
929.581230.191425.92805.642556.011649.973379.23
Sum of squares156.17118.31102.19179.957.1888.3943.29
Table 3. Results of calculations with regard to rheological parameters of the tested cement slurry.
Table 3. Results of calculations with regard to rheological parameters of the tested cement slurry.
Approximated shear rate γ ˙ S T A R T , [s−1]96
Shear rate γ ˙ , [s−1]   γ ˙ LOW =   γ ˙ 30   γ ˙ MID =   γ ˙ 60   γ ˙ TOP =   γ ˙ 90
51.10102.20153.31
Shear stresses τ, [Pa] τ LOW = τ 30 τ MID = τ 60 τ TOP = τ 90
8.1814.3119.42
Rheological parameters of the Vom Berg modelτyDG
1.244818.3547132.16
Unit flow resistances
dp/dl [Pa/m]
387
Real shear rate γ ˙ S T A R T , [s−1]
γ ˙ S T A R T ( γ ˙ L O W , γ ˙ T O P )
101.6
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Wiśniowski, R.; Orłowicz, G. Theory of the Vom Berg Rheological Model and Its Use in Cloud-Native Application. Energies 2022, 15, 4481. https://0-doi-org.brum.beds.ac.uk/10.3390/en15124481

AMA Style

Wiśniowski R, Orłowicz G. Theory of the Vom Berg Rheological Model and Its Use in Cloud-Native Application. Energies. 2022; 15(12):4481. https://0-doi-org.brum.beds.ac.uk/10.3390/en15124481

Chicago/Turabian Style

Wiśniowski, Rafał, and Grzegorz Orłowicz. 2022. "Theory of the Vom Berg Rheological Model and Its Use in Cloud-Native Application" Energies 15, no. 12: 4481. https://0-doi-org.brum.beds.ac.uk/10.3390/en15124481

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop