Technology Journal of Energy JET Volume 12 (2019) p.p. 31-39 Issue 1, April 2019 Type of article 1.01 www.fe.um.si/en/jet.html ANALYSIS OF PIPELINE VIBRATION ANALIZA VIBRACIJ V CEVOVODIH Jurij Avsec1R, Urška Novosel1 Keywords: pipeline vibration, pipeline fluid flow, continuous vibration Abstract Vibrations occur in almost all energy systems. This article presents an analysis of vibrations in pipelines under the influence of fluxes of displaced persons with a pressure difference or with the help of electromagnetic forces. For this purpose, we analysed pipelines of different diameters and the flow of crude oil within the pipelines. Povzetek Vibracije se pojavljajo v skoraj vseh energetskih sistemih. V članku je prikazana analiza vibracij v cevovodih pod vplivom tokov gnanih s tlačno razliko ali s pomočjo elektromagnetnih sil. V ta namen smo analizirali cevovode različnih premerov. V cevovodih smo obravnavali tok surove nafte. Energy devices are affected by many mechanical influences, including mechanical fluctuations. The analysis (analytical and experimental) of mechanical oscillations of energy devices is in many cases of paramount importance. In the present article, we have studied pipeline vibrations for a simply supported pipeline with the help of continuous vibration theory, [1-4]. Fluid flow in pipelines and channels is driven due to the presence of the electric field, magnetic field, or pressure-driven flow, and some other effects, [1] (Fig. 1). Electrohydrodynamics (EHD), known as electrokinetics, is the theory of the fluid dynamics of electrically charged fluids. It is the study of the motions of ionized R Corresponding author: Corresponding author: Prof. Jurij Avsec, Ph. D., Tel.: +386-7-620-2217, Fax: +386-2-620-2222, Mailing address: Hočevarjev trg 1, 8270 Krško, Slovenia E-mail address: jurij.avsec@um.si, urska.novosel@um.si 1 University of Maribor, Faculty of Energy Technology, Laboratory for Thermomechanics, Applied Thermal Energy Technologies and Nanotechnologies, Hočevarjev trg 1, SI-8270 Krško, Slovenia 1 INTRODUCTION JET 31 JurijAvsec, UrskaNovosel JET Vol. 12 (2019) Issue 1 particles or molecules and their interactions with electric fields and the surrounding fluid. The fundamental concept for magnetohydrodynamics (MHD) is that magnetic fields can induce currents in a moving conductive fluid, which in turn creates forces on the fluid and changes the magnetic field itself. This paper develops a formulation of crude oil motion, in macro channels and mini-channels with pipeline fluid flow. A mathematical model is developed for channels with rectangular and circular cross-sections. Q 2 FLUID FLOW DUE TO PRESSURE AND MAGNETIC EFFECTS Consider electromagnetic flow in rectangular and circular channels (Fig. 2). The charged surface of a channel wall may attract ions of the opposite charge in the surrounding fluid. The general form of the momentum equation for electrohydrodynamic flow is: dp d2u . 0 = —— + ^—- + i Bz dx dz y z + pvVv = -Vp + V(y.Vv) + ixB, (2. 1) where the last term presents the electromagnetic force, and i and B refer to the current density and magnetic field strength, respectively. For steady-state flow in a channel at small Reynolds numbers, the transient and inertia terms can be neglected, so Eq. (2. 1) is simplified in the next equation: 0 = -Vp + V(p,Vv)++ixT. (2. 2) 32 JET Analysisof pipeline vibration Figure 2: Rectangular and circular micro channels Laminar flow During electrokinetic flow in channels, the charged surface of a channel wall may attract ions of the opposite charge in the surrounding fluid. Assuming the fluid velocity, magnetic field and current density are orthogonal, the reduced momentum equation becomes: 0 -- dP + dU + iß dx dz2 y z (2. 3) The terms represent pressure, viscous and electromagnetic forces in the liquid. Using Ohm's Law to express the current density in terms of fluid velocity: d2u „2 dp ß—F + ueB¡u -—, dz e z dx (2. 4) where: oe and Bz refer to the electrical conductivity and magnetic field strength. For fully developed flow in a channel, the pressure gradient becomes constant and independent of the magnetic field strength. In terms of the Hartman number, MH (MH = aBae / r ), d2u ( MH2,q | - dp dz2 ^ a2 J dx (2. 5) Applying the no-slip boundary conditions at z=0 and z=-2a, the analytical solution of Eq. (2. 5) becomes: JET 33 Jurij Avsec, UrskaNovosel JET Vol. 12 (2019) Issue 1 dp i I zMH dx ! —H Mh2M (1 + 2e2MH )Mh2^ a2 I dp) e2MH dx (1 + 2e2MH )Mh2^ The mean velocity within the channel becomes: 1 2a =— J u(z)dz 2a 0 a2 \^")(~MH + TanhM)) Mh3M (2. 6) Non-dimensionlizing this result (z*=z/a, u*=u/ub), we obtain the following equation: Mh -1 + (1 + 2e2MH ) (1 + 2e2MH ) (-Mh + Tanh (Mh )) (2. 7) Without electromagnetic effects, equations (5-6) transform into the following expressions: I dp ) u(z) = ^dx (-2az + z2), u, = — f u(z)di ' b 2a l 3ß dp dx ! * u = 2(2-z M (2. 8) For the circular channel without electromagnetic forces, the governing equation is: dp I d2u 1 du — = M\ —2 +-- dx { dr r dr (2. 9) Solving the differential equation subject to boundary conditions, 4V (r2 -R2 ) , u, = —1— fu2xrdr = v ' kR2 J dx 8ß u* = 2 - 2r2. (2. 10) u b b 2 34 JET Analysisof pipeline vibration If we wish to calculate the velocity profile for MHD flow in a circular channel, we have to solve the next differential equation: dp dx du 1 du —F +-- dr r dr Mh R2 (2. 11) The analytical solution of equation (2. 12) is slightly more complicated. We have obtained the next solution of a differential equation with the boundary conditions (u(R)=0, u'(0)=0): R2\dpï| -BesselI t[r ] = - dx J AAv 0,M ]+Bessell[0,-R- M^Bessell [0,M ] (2. 12) kR -ju2tfrdr ■ R2 [d-jBessell [2,M] M2^BesselI [0,M ] ' (2. 13) 1 b Turbulent flow In the case of turbulent flow, the logarithmic distribution of velocity that is the most suitable is, according to experimental work, the following: U K (^IfUL + B. (2. 14) U K ^ With the upper equation, we could calculate the mean velocity in the circular channel for turbulent flow [7,8]: 1 Rf 1 n (R - r)pu* - f - lny" ' 'r " + B kR20 [v V \ 2Krdr = 1 , i 2 , Rpu „„3 | —u I — In+ 2B--I. (2. 15) 2 ^ K K J In the presented model, we used the following values: k = 0.41 and B=5. ub = JET 35 JurijAvsec, UrskaNovosel JET Vol. 12 (2019) Issue 1 3 VIBRATION OF PIPELINES The equation of beam vibration containing flowing fluid could be written with the next equation, [9-13]: d4 z 2 d2z „ d2z d2z , , El — = -p2 — - 2 pv--m—. (3. 1) dx4 dx2 dxdt dt2 Equation (3. 1) is obtained with continuous vibration theory. In equation (3. 1) the v means the average fluid velocity, p the mass of fluid per unit length, E the modulus of elasticity, and I the moment of inertia, [10-12]. The first term is the force caused by the change in direction of the velocity of the fluid due to the curvature of the pipe. The second term is the force associated with the Coriolis acceleration, and the last term represents the force that is related to the vertical acceleration of the pipe. For the simply supported pipeline, we obtain the next solution of the differential equation, [912]: 32 pvm 3l «if Tifv 3l * (flf T (3. 2) 8 pvm In equation (3. 2), m means the angular velocity of the pipe relative to its length, and l the length between supports. With equation (3. 2), the relation between v and m is expressed. If we limit equation (3. 3) calculates the critical velocity at which static divergence occurs: Ef '3.3, V c 4 RESULTS AND VIBRATIONAL ANALYSIS Figure 3 shows analytical results for crude oil SAE 15W-40 flow at 200C through the pipeline with the next data: D=0.2m, ^=0.287 Pas, e=878.7 kg/m3. The fluid flow is also turbulent at small pressure differences. Figure 4 shows analytical results for crude oil flow through the pipeline with the following data: D=0.002m, ^=0.287 Pas, e=878.7 kg/m3. The fluid flow is also laminar at high-pressure differences. 36 JET Analysisof pipeline vibration ; (Pa/m] Figure 3: Menn pipe velocity and Reynolds number developed by the turbulent model u [m/s] 10000 SOOO 6000 4000 2000 dp/dx [Pa/m] Figure 4: Menn pipe velocity and Reynolds number developed by lnminnrfluid model Figure 5: Vibrational characteristics for steel pipeline containing crude oil (D=0.2m, ^=0.287 Pas, e=878.7 kg/m3, d=0.02m, L=10 m) JET 37 Jurij Avsec, UrskaNovosel JET Vol. 12 (2019) Issue 1 R Figure 6: Velocity profile with slip (6*=0.5) at Mh=5 (blue line), Mh=10 (sed line) and Mh=30 (gseen line)in circular minichannel without slip. The analysis shows the velocity profile of crude oil in the pipeline. The inner diameter of the pipeline is in the first case 0.2 m in the second case 0.002 m. In the pipeline with diameter 0.2 m, the fluid flow is turbulent; in the pipeline with the diameter of 0.002 m, the fluid flow is laminar. Even more interesting are results for the calculation of natural frequencies for the first and second mode vibration regarding pipelines with crude oil. As we see from Figure 5, the natural frequencies of the pipe decrease as the fluid flow increases; from the figures, the influence of the average velocity profile in macro and mini-pipelines is also evident. Figure 6hows velocity profile for MHD flow in the circular channel. 38 JET Analysisof pipeline vibration 5 CONCLUSION AND OUTLOOK The present article shows the vibrational analysis of pipes conveying crude oil SAE 15W-40 flow. The fluid flow was calculated for laminar and turbulent flow. The vibrational characteristics have been calculated with vibration continuous pipeline theory. References [1] V. Singhal, S. Garrimela, A. Raman: Microscnle pumping technologies fvr micrvchnnnel cvvling systems, Appl. Mech. Rev., Vol. 57: 191-222, 2015 [2] H. Chen: Vibration vf n Pipeline Containing Fluid Flvw with Elastic Support, MsC Thesis, Ohio University, 1991 [3] S. Rao: Vibration vf cvntinvus systems, Wiley, 2007 [4] G. F. Naterer, O. B. Adeyinka: Microfluidic Exergy Loss in n Non-Polarized Thermomngnetic Field, International Journal of Heat and Mass transfer, Vol. 48: 39453956, 2005 [5] G. D. Ngoma, E. Erchiqui: Heat Flux effects on Liquid Flow in n Microchnnnel, International Journal of Thermal Sciences, Vol. 46: 1076-1083, 2007 [6] H. P. Mazumdar, U. N. Ganguly, S. K. Vanketesan: Some Effects on Magnetic Field on the Flow of n Newtonian Fluid through n Circular Tube, Indian J. of Pure Applied Mathematics, 1996, Vol. 27: 519-524, 1996 [7] http://www.roymech.co.uk/Related/Fluids/Fluids_Viscosities.html [8] NPTEL: Principles of fluid mechanics, http://nptel.ac.in/courses/101103004/ [9] H. Dodds, H. L. Ruanyan: Effect of high velocity fluid flow on the bending vibrneivnAnnnd static divergence of n simply supported pipe, NASA TN D-2870, 1965 [10] H. Chen: Vibration of n pipeline containing fluid flow with elastic support, Master science thesis, Ohio University, 1991 [11] M. P. Paidoussis, S. J. Price, E. Langre: Fluid structure interactions, Cambridge Press, 2010 [12] R. D. Blevins: Flow-Induced Vibration, Van Nostrand Reinhold, 1990 JET 39