Paper received: 27.10.2008 Paper accepted: 03.04.2009 Determination of Pressure Losses in Hydraulic Pipeline Systems by Considering Temperature and Pressure Vladimir Savic 1 - Darko Knezevic 2 - Darko Lovrec 3 - Mitar Jocanovic 1 - Velibor Karanovic 1 1 Faculty of Tehnical Sciences, University of Novi Sad, Serbia 2 Faculty of Mechanical Engineering, Univesity of Banja Luka, Bosnia and Herzegovina 3 Faculty of Mechanical Engineering, University of Maribor, Slovenia Generally accepted methods for calculating pressure losses within flat pipelines as presented in literature and used in praxis, are based on the Reynolds number, which considers the viscosity and density of fluid, internal pipe friction coefficient, pipe geometry, and oil circulation velocity. Such an approach contains serious inconsequentiality. Namely, the only nominal values for viscosity and density are considered in the calculation, which differs substantially from real conditions. This often leads to inaccurate calculations of pressure losses. A numerical model has been developed in the work presented in the paper. It takes into account the actual changes in density and viscosity under the current oil pressure and temperature in order to overcome the above weaknesses of standard calculation procedures. Such an approach is novel and provides new capacity for an accurate pressure drop analysis of advanced hydraulic systems. © 2009 Journal of Mechanical Engineering. All rights reserved. Keywords: hydraulic pipeline, pressure losses, temperature, presure 0 INTRODUCTION Energy dissipation due to the frictional resistance of fluid flow within a hydraulic pipeline system is represented by any decrease in pressure. The pressure loss is determined according to the well-known Darcy-Weisbach equation: L v2 Ap = X - p — (1) d 2 where: L - pipeline length, d - diameter of the pipe, v - the fluid flow velocity, p - density of the hydraulic fluid and 1 = l{^ Re, d j friction factor 1 which depends on the Reynolds number Re, and pipe roughness. In the cylindrical tubes, the friction factor is usually determined by using the Moody diagram (Fig. 1). This diagram depicts the Reynolds number against the friction factor, and describes the corresponding fluid regimes. The fluid flow through a pipe can be either smooth or rough, depending on flow conditions. Various factors determine the nature of flow. In principle, the flow can be stated as either laminar flow, i.e. steady, smooth flow, or turbulent flow, i.e. the flow is disturbed. From a practical point of view, the Reynolds number indicates if a flow is laminar or turbulent. In the case of turbulent flow, the friction factor 1 depends on the Reynolds number, as in the case of laminar flow as well as on the coefficient of the relative roughness - the relationship between the absolute roughness of the pipe's inner wall surface in contact with the fluid, and the diameter of the pipe (s/d), bearing in mind that the pipe can be either perform as a smooth, rough or somewhere in-between -liminal (transitional). Several types of empirical equations that can help to determine the friction factor can be found in literature. Some fairly useful ones are given below [1] and [2]. a) Laminar flow area at Re < 2000: 1=64/Re (2) b) Turbulent flow area within the smooth pipe of a hydraulic system (according to Blasius [2]) at Re > 4000: 1=0.3164/Re025 (3) c) Transitional flow type area within range 2000 < Re < 4000: 1 = 3.9 Re + 0.0242. (4) * Corr. Author's Address: University of Novi Sad, Faculty of Technical Sciences, Trg D. Obradovica 6, 21000 Novi Sad, Serbia, savicv@uns.ns.ac.yu Fig. 1. Moody diagram for the determination of flow regimes with regard to internal friction coefficient [2] d) The Reynolds number value can be estimated using the following equation: Re= vd/v (5) where v [m2/s] - kinematic viscosity, d [m] - pipe diameter, v [m/s] - average fluid flow velocity. Most of the authors [1] legitimately assume that all the hydraulic pipeline systems are smooth because: a) Pipes are produced by precise drawing with low absolute roughness values 0.01 to 0.002 mm (average value 0.006 mm). b) The applied Reynolds number is most often in the order of medium values, i.e. between 1000 and 20000 so the friction factor value is at the limit of what is considered smooth (the lowest curve in Moody's diagram represents smooth hydraulic pipes). Numerical verification of the aforementioned statement is based on the calculation within the following boundary values: a) mostly used hydraulic pipelines range in diameter from 10 to 50 mm, so that the relative roughness s/d is within the range of 0.0006 to 0.00012 which is at the lower end of the relative roughness, as depicted in Moody's diagram (see Fig. 1); b) usually fluid flow velocity v =1 to 6 m/s; c) common kinematic viscosity at operating temperature (16 to 54)-10"6 m2/s are borderline values of the optimum viscosity for the e.g. axis piston pumps. For the borderline areas, the Reynolds number falls within the range indicated in Table 1. Table 1. Dependency of Reynolds number of the pipeline and fluid properties Dimension Number value Pipe diameter [mm] 10 50 Flow speed [m/s] 1.0 6.0 1.0 6.0 Kinematic viscosity [mm2/s] 16 54 16 54 16 54 16 54 Reynolds number [ - ] 625 185 3750 1110 3125 925 18750 5550 When analyzing the values of the Reynolds number, two conclusions can be drawn. The Reynolds number will increase with diameters d > 50 mm, flov speed v > 6.0 m/s, and kinematic viscosity v < 16 mm2/s, whereas the Reynolds number will decrease with d < 10 mm, v <1.0 m/s and v < 54 mm2/s. It follows, that the friction factor of the mineral hydraulic oil-flow through hydraulic pipeline can be determined with satisfactory accuracy, using the Eqs. (2) to (4) and the pipeline of the hydraulic system may be treated as smooth. 1 INFLUENCE OF TEMPERATURE AND PRESSURE CHANGES ON PHYSICAL-CHEMICAL PROPERTIES OF HYDRAULIC FLUID The characteristic physical values of all hydraulic fluids which have a deciding influence on the loss of energy are density p and viscosity v of the fluid. Both are temperature (T) and pressure (p) dependent values: p = p(p,T), v = v(p, T). 2.1 Density of Fluid Density is defined as the mass of fluid per unit volume: P = V (6) where: m [kg] - mass, V [m ] - volume. The density of a hydraulic fluid is measured under atmospheric pressure at T = 288 K. Any increase in temperature results in increased fluid volume, whereas the fluid volume decreases according to any rise in pressure (and contrarily in the reverse order of temperature and pressure changes). Density changes as a consequence of temperature and pressure changes can be given in the form of a total differential: dp J P dT + I ^ V dT J p I dp or in a rewritten form: d p _ ,m 1 dp dT +- p K -dr (7) T from where results Op - volumetric temperature fluid expansion coefficient by constant pressure: 1 ( dp a p = - p V dT and Kt - isothermal compression module Kt = p dp d p jt (8) (9) 2.1.1 The Effect of Temperature on the Oil Density Using the law of mass conservation (m = pV = const.) according Eqs. (8) and (9), two equivalent models can be expressed by the following equations: 1 ( dV a„ =— I- p V V dT and Kt = -V dp dV (10) (11) Coefficient ap serves to define the relative change of fluid (volumetric) density. This is dependent on the mean molecular weight of the hydrocarbon. Figure 2 shows a diagram of coefficient change, as drawn-up on the bases of the tested numerical values [3]. Based on the experimentally verified Equation [4]: d p\ —t— I = -ap p = const dT J p pF it is possible to derive an equation which describes the relationship between the hydraulic oil density for any temperature, and the atmospheric pressure, using the values that apply to the density at T = 288 K [4]: px =P15 [1 -ai5 (T - 288)] (12) where the value of a15=0.0007 is evident on the diagram in Fig. 2. 7.6 7.2 g 6.8 L 6.4 a 6.0 5.6 20 30 40 50 60 70 80 90 100 T [°C] Fig. 2. Changes in coefficient ap underconstant pressure pa -"-p= 100 bar — p= 200 bar - -p= 300 bar -*-p= 400 bar -*-p= 500 bar m T 2.1.2 The Effect of Pressure on the Hydraulic Oil density As mentioned above, the fluid volume reduces with any increase in pressure, while its density becomes greater at the same time. Compression module can be expressed in terms of the previously-stated equations: KT =p (=-v (*l y d p)T y dV JT so that any alteration in density, as a function of pressure change, for both the constant values of compression module KT and constant temperature, is worked out by the formula below: Ap = p ^ (13) kt whereby it follows that pp =px +AP = Px +Px-Kp = Px ^ + ] . (14) Coefficient (KT) change diagram modeled on the bases of the tested numeric values [2] is presented in Fig. 3. -■-T=150C -•-T=30oC ■ T=50°C -*-T=70°C -«-T=100°C 10000 100 200 300 P [bar] 400 T= 15°C T= 30°C T= 50°C -T= 70°C T= 100°C 500 100 400 500 200 300 P [bar] Fig. 3. Change of isothermal (KT) and adiabatic (KS) compression module as a function of operating pressure and oil temperature 2.2 Viscosity Viscosity represents fluid resistance to shear or flow and is a measure of the adhesive/cohesive or frictional fluid's property. Used in hydraulics, there is a concept of dynamic viscosity, the value of which is determined by Newton's equation concerning parallel linear laminar flow. A ^ dx ( ) 2 where: F [N] - is side force; A [m ] - side force surface; dv/dx [1/s] - velocity of angular deformation. Dynamic viscosity is expressed in Pascal seconds [Pa-s], which corresponds to [kg/m-s]. In technical practice, the term kinematic viscosity is used, which is defined as the ratio between dynamic viscosity and fluid density: v=n p' (16) 2.2.1 Influence of Temperature Change on the Viscosity of the Hydraulic Oil The kinematic viscosity of industrial oils is measured under atmospheric pressure at 40 °C. It is classified into VG viscosity classes according to the ISO standard. Hydraulic fluids of class ISO VG 46 are most frequently used in praxis. They are followed by fluids classed as VG 32, and less commonly by VG 68 or 22. Several formulae, derived for determining viscosity as a function of temperature, can be found in relevant literature, e.g. [6]: U = b e -a-T, U = (b+TA )c ■ according to Reynolds (17) - according to Slotte (18) T c (v + a) = b d A - according to Walther (19) ju = a e(Ta c) - according to Vogel. (20) The coefficients a, b, c, d contained in the Eqs. (17) to (20) are affected by the properties of individual hydraulic fluids, with TA signifying the absolute temperature. The most accurate among them is Eq. (20) but the more frequently used one is Eq. (19), which also forms the basis for the ASTM viscose temperature diagram's design [5]. The measurement for the change of kinematic viscosity with temperature is commonly determined by the viscosity index (VI), a non-dimensional number. Oils of faster viscosity alteration are classified as HM oils 0 a b (VI = 90), and those of slower alteration as HV oils (VI = 140). Figure 4 shows a diagram of viscosity alteration as a function of temperature and in relation to HM oils 32, 46, and 68. The diagram represents the results of tests for oil's viscosity at 40 and 100 °C [5]. 150 100 80 60 40 30 20 \ V \ v ' s. \ S s \ HM 68 V 1 V V HM X 46 32 32 N § s Os. V, v S V o 0 10 20 30 40 50 60 70 80 90 100 Temperature T [oC] Fig. 4. The alteration of viscosity as a function of temperature for oils HM 32, 46 and 68 2.2.2 The Effect of Pressure on the Viscosity of Hydraulic Oils Oil viscosity increases with any concurrent increase in pressure. The speed of this change depends on the molecular structure of the oil, meaning that the speed will be higher where the oil viscosity is lower. Numerical value of kinematic viscosity alteration as a function of the change in pressure can be worked out in two different ways: a) directly by kinematic viscosity calculation: Vp = v (1 + k p) (21) where: k = 0.002 to 0.003 stands for the change of viscosity coefficient; b) indirectly, i.e. obtaining a dynamic viscosity value by first using the Barus equation: M = M0 eap (22) where: a = a(p,T) - coefficient of viscosity depending on pressure, ^ - value of dynamic viscosity under atmospheric pressure, and at an appropriate operating temperature (see Fig. 5). The kinematic viscosity calculation drawn from this model is done in relation to the calculated value of dynamic viscosity and the density at the operating temperature and pressure using the following formula: np vp = — = • Mo e ap Pp +Kt (23) which follows from Eqs. (14), (16) and (23). 3.2 3.0 2.8 2.6 2.4 2.2 2.0 1.8 1.6 1.4 1.2 p=pa p=500 bar 20 40 60 T [°C] 80 100 Fig. 5. The alteration of the coefficient a as a function of any change in temperature, and operating pressure ranging from pa to 500 bar for mineral oils with paraffin base structure [5] 3 NUMERICAL VALIDATION OF THE PROPOSED MODEL The correlation of temperature and pressure with changes in the key properties of fluids, which should be taken into account when calculating pressure losses within hydraulic systems, is presented in Chapters 1 and 2. It remains unclear whether change in density and viscosity substantially effects pressure loss in hydraulic systems, and whether they should be considered while designing a hydraulic system; and if so, how to set up a general model for calculating the exact loss of pressure. It only makes sense to answer the question with regard to a representative model of the pipeline complying with defined criteria, as follows: - pipeline diameters 15, 25 and 40 mm alternately; - pipeline length 1 m; - velocity of flow 6 m/s; - type of hydraulic mineral oil HM46, viscosity index 90 to 100; - operating temperature of the hydraulic oil 50 °C; - density of the hydraulic oil 875 kg/m3 at 15 °C; - operating pressure 100, 200, 300 and 400 bar; - the nature of the pressure change: adiabatic. 0 3.1 Step 1 - Calculation Model for Determining Oil Density at Operating Temperature 50 °C, and Atmospheric Pressure. The volumetric temperature expansion _4 coefficient equals (Fig. 2) ap=7.0-10 , which makes the fluid density at operating temperature 50 °C and under atmospheric pressure equal to: P50 = P15 [1 - 0.0007 (T - 288)]= 853.6 kg/m3 . 3.2 Step 2 - Calculation Model for Determining Oil Density under Operating Pressure The compression module for the adiabatic change of pressure (Fig. 3) equals: a) pi = 100 bar ^ KSi = 17400 bar; b) p2 = 200 bar ^ KS2 = 18200 bar; c) p3 = 300 bar ^ KS3 = 19200 bar; d) p 4 = 400 bar ^ KS4 = 20500 bar. Oil density under the defined pressures equals: P100 = Po I 1 + Pi K Ti = 858.5 kg/m3 P200 = 862.9 kg/m3 P300 = 866.9 kg/m P300 = 870.3 kg/m3 . The density change curve as a function of operating temperature and temperature, is shown in Figure 6. p[kg/m p[bar] Fig. 6. Density change as a function ofpressure change at temperature 50 °C 3.3 Step 3 - Calculation Model for Determining Viscosity at Operating Temperature Kinematic viscosity of oil HM 46 is 46T0"6 m2/s at 40 °C and 30-10"6 m2/s at 50 °C. 3.4 Step 4 - Calculation Model for Determining Viscosity under Operating Pressures Kinematic viscosity at 50°C and under atmospheric pressure equals 25-10-6 m2/s, whereas pressure coefficient is 0.0025. Therefore, viscosity of the determinate operating pressures equals: - pi = 100 bar: v100 =v (1 + k ■ p1) = 37.5 -10-6m/s2 ; - p2 = 200 bar; v200 = 45 ■ 10-6 m/s2 ; - p3 = 300 bar; v300 = 52.5 ■ 10-6 m/s2 ; - p4 = 400 bar; v400 = 60 ■ 10-6 m/s2 . 3.5 Step 5 - Calculation Model for Determining Pressure Loss a) Viscosity corrected in accordance with the level of operating temperature. The starting parameters: v50 = 30-10-6 m2/s; p= 875 kg/m3; v = 6 m/s; d = 25 mm; L = 1 m. b) Viscosity corrected in accordance with the level of operating temperature and operating pressure. The starting values: v = 6 m/s; d = 25 and 15 mm; L = 1 m. - pa ^ V0 = 30-10-6 m2/s; P0 = 875 kg/m3; - px = 100 bar ^ v100 = 37.5-10-6 m2/s; P100 = 858.5 kg/m3; - p2 = 200 bar ^ V200 = 45-10-6 m2/s; P200 = 862.9 kg/m3; - p3 = 300 bar ^ v300 = 52.5-10-6 m2/s; P300 = 866.5 kg/m3; - p4 = 400 bar ^ v300 = 60-10-6 m2/s; P400 = 870.3 kg/m3. An example of the calculation model for determining turbulent flow: step 5a • The Reynolds number equals: Re = — = 5000 . v • The friction coefficient equals: I- 03164 = 0.0376. Re 0.25 • The pressure decrease within tube network equals: L v2 A p = 1--p = 23686 Pa . d 2 An example of the calculation model for determining the flow during the liminal phase: step 5b Re = — = 3000 2 = 3.9-10-6 • Re + 0.0242 = 0.0359 . L v2 A p = 2--p = 36435 Pa . d 2 An example of the calculation model for determining laminar flow: step 5 c Re = WA = 1500 v 64 1 = — = 0.0427. Re L v2 A p = 1--p = 44548 Pa . d2 Summary of the overall calculation for both pipe diameters is shown in in Figure 7. 50000 ]a 40000 p < s s 30000 o er ru s 20000 s er pr 10000 0 ^ d=15 mm d=25 mm W/T/TTrf, d=40 mm 100 200 300 operating pressure p [bar] 400 Fig. 7. Pressure loss as a function of the operating pressure level and the pipe diameter 4 CONCLUSIONS Generally accepted methods for calculating pressure losses within flat pipelines, as presented in technical literature and used in praxis are based on the Reynolds number, which considers viscosity and density of fluid, internal pipe friction coefficient, pipe geometry, and oil circulation speed. Technical literature treats the aspect of density and viscosity alteration as a function of temperature, and pressure change as an isolated issue with only tangential links to the problems of designing and calculating a hydraulic system. Nominal values for viscosity (at 40 °C) and density (at 15 °C and normal atmospheric pressure), which differ substantially from real conditions, are considered. Consequently, the derived values cannot be regarded as highly accurate, leading to inaccurate calculation of pressure losses (in some cases up to 200%). The numerical model, which takes into account the actual changes of density and viscosity at the current oil pressure and temperature in order to overcome the above weakness, is suggested in this paper. Such an approach is novel and provides a new capacity for an accurate pressure drop analysis of advanced hydraulic systems. 5 REFERENCES [1] Akers A., Gassman M., Smith R. Hydraulic Power System Analysis, CRC/Taylor & Francis, 2006, ISBN 0-8247-9956-9. [2] Munson B.R., Young D. F., Okiishi T. H. Fundamentals of fluid mechanics, John Wiley & Sons, 2006. ISBN 978-0-47167582-2. [3] Frauenstein, M. Directives about hydraulic fluids for design engineers, Mobil Oil AG, Hamburg, 3. Edition, 1980. (in German). [4] Knezevic, D. Effects of gap geometry in hydraulic components on hydraulic system efficiency, Ph.D. Thesis, University of Novi Sad, 2007. (in Serbian). [5] Knezevic, D., Savic, V. Mathematical modelling of changing of dynamic viscosity asa function of temperature and pressure of mineral oils for hydraulic system, Facta Universitatis, Series: Mechanical Engineering, 2007, vol. 4, no. 1, 2006, p. 2734, Niš. [6] Stachowiak G.W., Batchelor A.W., Engineering Tribology, ButterworthHeinemann, 2nd edition, 2001, ISBN 0-75067304-4. 0