Photovoltaic Cells and Systems: 102-117 SQU Journal for Science, 16 (2011) © 2011 Sultan Qaboos University 102 Dynamical Properties of Omani Crude Oils for Flow Through a Vertical Annulus and a Cylindrical Pipe Sayyadul Arafin* and S.M. Mujibur Rahman Department of Physics, College of Science, Sultan Qaboos University, P.O.Box 36, Al-Khoud, 123 Muscat, Sultanate of Oman, *Email: sayfin@squ.edu.om. الخواص الديناميكية للنفط العماني الخام عند سريانه عبر أنبوب حلقي أو أسطواني سيد العارفين و الشيخ مجيب الرحمن عتماد الكثافة واللزوجة على درجة الحرارة لعينات من النفط الخام ابالتحقيق في مدى دراسة اية هذه القمنا في بد :خصمل تم استخدام بيانات القياسات في دراسة الخواص الديناميكية لهذه ثم تلفة في سلطنة عمان. مخ هيدروكربونيةجمعت من حقول تحت تأثير الجاذبية والضغوط المختلفة عند درجات حرارة مختلفة. لقد قمنا بنمذجة سريان العينات الهيدروكربونيةالمواقع التقريب النيوتوني لسريان المائع الطبقي وذلك خدامباستالمختلفة لخام النفط عبر أنبوب رأسي )أ( حلقي )ب( أسطواني إليجاد خواص سريان هذه العينات . تم تطوير برنامج حاسوبي بحيث يمكن إعتماد سريان المائع على درجة الحرارة من . 0222وأعلى من 0222أقل من الفصل بين وضع السريان الطبقي ووضع السريان المضطرب بالنسبة ألعداد رينولد ومعدل سريان الكتلة وقوة اللزوجة على السطح الجامد ليست بالجديدة اعتمدت للتشكيل الجانبي للسرعة ماذج التيالن مختلفة؛ AIPعينات نفط خام عمانية مع قيم الستخدام)المستحدثة( لكن الحسابات المقدمة في هذا البحث تهدف بالتحديد يناميكية لهذه العينات المحددة في إطار التقريب النيوتوني . يمكن أن تكون قى النتائج المحسوبة الضوء على الخواص الدتونل ذاث فائدة في أغراض مختلفت مثل استخراج النفط الخام ونقله عبر للخواصالخواص الفيزيائية المقاسة والحسابات الالحقة األنابيب. ABSTRACT: We have initially investigated the temperature dependence of density and viscosity of a number of crude oils, collected from various hydrocarbon reservoirs in Oman. The measured data are then utilized to investigate the flow dynamics of these hydrocarbon fluids under gravity and applied pressures at various temperatures. We have modeled the flow of the various crude oil samples through a vertical (a) annulus and (b) cylindrical pipe - all treated within the Newtonian fluid flow approximation of a laminar flow - to investigate the flow properties of these samples. A computer program is developed so that the temperature dependence of the fluid flow distinctly separates the laminar mode from a turbulent mode with respect to Reynolds numbers within the ranges Re<2000 and Re>2000. The adopted models of the velocity profiles, mass rate of flow and viscous force on the solid surface are not novel, but the present calculations aim to specifically use the various Omani crude oil samples with various AIP values; the calculated results shed some light on the dynamics of these specific samples within Newtonian approximation. The measured physical properties and the subsequent calculations of the relevant dynamical properties might be useful for various purposes e.g. extraction and transportation of crude oils through pipes. KEYWORDS: Dynamics; Flow properties; Hydrocarbon fluids; Cylinder and annulus. DYNAMICAL PROPERTIES OF OMANI CRUDE OILS 301 1. Introduction emperature dependence of thermo-physical properties, e.g., densities and viscosities of hydrocarbon fluids plays an important role in many fields of petroleum industries including enhanced oil recovery, oil purification, and transportation of produced fluids (Ahmed, 2000). When producing heavy oils, the high viscosity manifests as one of the impediments to recovering these oils from the oil rigs. The temperature dependence of the transport properties influences the relevant flow dynamics, which in turn affect the flow mechanism (Li et al. 2004). We present here a model study of the temperature dependence of certain dynamic properties such as velocity distribution, mass rate of flow and viscous force on the solid surface through a vertical annulus and a cylinder. The present investigation primarily aims at understanding the profile of a number of dynamic properties of Omani crude oils through a vertical annulus at various temperatures and applied pressures. We may mention that fluid flow through an annulus has wide applications in various branches of science and technology including nuclear reactor engineering, oil and gas production, aero-engines, turbo-machinery, chemical engineering, steam generators and heat exchangers. Knowledge of the dynamics of energy flows either in the form of heat transfer or fluid flows through a vertical annulus or normal cylinder under atmospheric and pressurized conditions, can help us understand the re-flooding phenomenon during the emergency cooling in a water-cooled reactor (Shiotsu and Hama, 2000). It is particularly important in the petroleum industry because annular flow of liquid and gas in a pipeline segregates the material of lighter molecular weight by restricting its flow down the center of the pipe while allowing the material of heavier molecular weight to form a thin film and flow along the pipe wall. The lighter mass fluid or gas can also be in the form of a mist or colloidal suspension known as an emulsion. The interface between the flowing materials may not be entirely precise, and can involve gas and liquid mixtures. It has been suggested (Prada and Bannwart, 2001) that the use of a core-annular flow pattern may be made attractive as an artificial lift method in heavy oil wells by inducing the flow pattern by the lateral injection of relatively small quantities of water in order to get a lubricated oil core along the pipe. Concurrently we have also investigated the flow properties of a number of Omani crude oil samples through a vertical cylinder. The research work on vertical upward core flow is nonetheless scanty; particularly, the theme of research on flow dynamics with Omani crudes is novel and in its infancy (Arafin et al. 2011). However, the works of Shertok (1975), Bai et al. (1996) and, Ho and Li (1994) on various samples are worth mentioning here. The prime interest in annular flow is to study the parameters relevant to the transport of fluids through straight pipes whether horizontal or vertical. However, presently, interest in annular flow goes beyond straight pipes to include intersections in pipe line networks, such as T-junctions, where phase segregation is likely to exist (Adechy and Issa, 2004). We have used here the Newtonian fluid approximation relevant only to laminar and linear fluids flows. Certainly crude oils are far from Newtonian fluids. It is relevant to mention here that the Newtonian approximation is highly limited to homogeneous, isotropic and non-compressible systems where stress-velocity gradient relations are non-linear. Above all, any non-Newtonian fluid has a non-zero term on the right hand side of the momentum-balance equation. Consideration of non-Newtonian approximation will certainly necessitate the non-linear form of the Navier-Stokes equation. Thus it is apparent that the present approximation yields qualitative results only. Thus for a non-laminar and nonlinear case, the situation will be certainly complex, and that will be considered in a future endeavor. In the present investigation, initially, the densities of the samples were measured by using the Anton Paar density meter (DMA 5000) within temperatures from 20°C to 70°C, at steps of 5°C. Subsequently, the kinematic and dynamic viscosities for all the samples were measured within this range of temperatures using a state-of-the-art Cannon-Fenske viscometer. Finally, the density and viscosity data were utilized to calculate the velocity distribution, maximum and average velocities of the bulk samples, and their mass flow rate under various thermal conditions. The calculations were constrained by prefixing the Reynolds number, Re, at 2000, ensuring a laminar flow as suggested by Bird et al. (2002). T http://www.springerlink.com/content/7581u35430540303/fulltext.html#CR1 http://www.springerlink.com/content/7581u35430540303/fulltext.html#CR6 http://www.springerlink.com/content/7581u35430540303/fulltext.html#CR3 SAYYADUL ARAFIN and S.M. MUJIBUR RAHMAN. 301 The layout of the paper is as follows: In Section 2 we describe the experimental details for measuring the densities and viscosities of the samples followed by a display of graphical results for temperature dependence of density and viscosity, and we end with fitting formulas. Basic model-based formulation for the dynamical quantities is presented in Section 3. Details of calculated results are presented and discussed in Section 4 followed by some concluding remarks in Section 5. 2. Experimental measurements In this section, we describe the experimental details of measuring the densities and viscosities at various temperatures for all the five samples. 2.1 Density measurements at various temperatures We briefly introduce here the classification of the crude oil samples that were procured from different oil fields in Oman. The two samples (Oman Export and Receive Line) comprising mixtures of crude oils were collected from the Petroleum Development of Oman (PDO). The other three samples were collected from Erad field, Mabruk field, and Zal-41 field located in different regions in Oman. We have used the American Petroleum Institute (API) oil gravity number to classify our oil samples: 0 141.5 131.54 API    (1) Here, 0  is the density (g cm−3) measured at temperature 15.6°C and at atmospheric pressure. The API numbers usually vary from 5 for very heavy oils to nearly 100 for light condensates (Batzle and Wang, 1992). Clearly, our samples may be classified into three types. The heavy type is represented by the Erad (API = 19.19) sample; to the intermediate type belong the Oman Export (API = 33.41) and Receive Line samples (API = 34.13), while the light type consists of Mabruk (API = 39.47) and Zal-41 (API = 40.89) samples. The densities of the samples were measured by using the Anton Paar density meter (DMA5000) as shown in Figure 1. The unit consists of a U-shaped oscillating tube, a system for electronic excitation, frequency counting, and a display. The injected sample volume is kept constant and vibrated. The density is calculated based on a measurement of the sample oscillation period and temperature. The temperature was controlled to ± 0.01°C during the measurement using a built-in thermostat. By measuring the damping of the U-tube’s oscillation caused by the viscosity of the filled-in sample, the instrument automatically corrects the viscosity- related errors. The Anton Paar density meter is calibrated to measure density to an accuracy of ± 5 × 10 −3 kg m −3 . This device was used to measure the density in the range of temperature varying from 20°C to 70°C, through temperature increments of 5°C. From the practical point of view, density data of crude oil as a function of temperature provide important information, which is useful for various industrial applications ranging from exploration to refining and transportation. The density data of crude oil were plotted in Figure 2. These data can be adequately represented by the equation:  ( ) r r T m T T    (2) where ρr is the density at 20°C, m is the slope of the density versus temperature curve (dρ/dT), T is the temperature (°C), and Tr is the room temperature (20°C). http://www.springerlink.com/content/7581u35430540303/fulltext.html#CR2 http://www.springerlink.com/content/7581u35430540303/fulltext.html#Fig1 http://www.springerlink.com/content/7581u35430540303/fulltext.html#Fig2 DYNAMICAL PROPERTIES OF OMANI CRUDE OILS 301 Figure 1. Anton Paar density meter (DMA 5000). 20 30 40 50 60 70 780 800 820 840 860 880 900 920 940 R 2 =0.99988 Erad R 2 =1 Oman R 2 =1 Receive R 2 =0.9999 Mabruk D e n s it y ,  ( k g m -3 ) Temperature, T( o C) Erad Oman Receive Mabruk Zal-41 R 2 =1 Zal - 41 Figure 2. Density versus temperature for various collected crude oils. The error in density, Δρ, is 0.005 kg m -3 .The equations for density as a function of temperature obtained by a fitting procedure are written in the form of Eq. 2 as follows: Erad sample: ( ) 933 . 008 0.638 ( ) r T T T    (3a) Oman Export sample: ( ) 868 . 980 0.699 ( ) r T T T    (3b) Receive Line sample: ( ) 851 . 426 0.701( ) r T T T    (3c) Mabruk sample: ( ) 824 .645 0.683( ) r T T T    (3d) Zal_41 sample: ( ) 817 .772 0.686( ) r T T T    (3e) SAYYADUL ARAFIN and S.M. MUJIBUR RAHMAN. 301 2.2 Viscosity measurements at various temperatures Two types of Cannon-Fenske viscometers, type 350 and 300 with calibration constant 5 × 10 −6 m 2  s −1 and 2.5 × 10 −6 m 2  s −1 respectively, were used to measure the viscosity. A Cannon-Fenske viscometer type 350, was used to measure viscosity of the Erad sample, the heaviest crude among the five samples. For the rest of the crude oil samples, a Cannon-Fenske viscometer type 300 was used. The working principle of the viscometer is based on the fact that the average velocity of steady flow in a round tube depends inversely on viscosity. The viscometer determines the kinematic viscosity by timing the fluid flow through a capillary tube as it passes between two etched lines on the glass wall. In our experiment, the viscometer was inserted into a constant temperature bath whose temperature was controlled to ± 0.1°C using a Haake D8 circulation thermostat. To establish the efflux time, the thermostat was set at a desired temperature and the sample liquid was allowed to fall freely down past an upper mark and the time taken for the meniscus to pass the lower mark was measured. The kinematic viscosity (ν k) was determined by multiplying the measured transit time of the fluid column in seconds with the calibration constant. The error of the viscosity measurement is < 0.35%. The error in the measurement of viscosity is too small to be shown in the plot (Figure 3a - b). From the Figures 3a and 3b, it is clear that the kinematic viscosity decreases exponentially (Arrhenius type) as the temperature increases; this correctly reflects the effect of temperature on the crude oil viscosity. From the plot, it is also evident that crude oil samples with lower API values have higher kinematic viscosity. An identical behavior of the viscosities of various Omani crude oils has been reported previously (George et al. 2006). The equations of kinematic viscosity as a function of temperature were obtained by exponential fit of the data of various crude oil samples. These are listed as follows: Erad: ( /11.98) ( ) 106.12 21384.40 T k T e    (4a) Oman Export: ( / 23.72) ( ) 3.21 47.46 T k T e    (4b) Receive Line: ( / 23.88) ( ) 4.66 36.44 T k T e    (4c) Mabruk: ( /17.95) ( ) 3.27 26.70 T k T e    (4d) Zal-41: ( / 28.6) ( ) 2.05 16.19 T k T e    (4e) The dynamic viscosity, μ(Pa s), is determined from the product of kinematic viscosity, ν, and density, ρ (μ = νρ). It is noted that the viscosity of a fluid is highly temperature-dependent. The viscosity of the crude oil decreases exponentially as the temperature increases from 20°C to 70°C. The present results are consistent with the correlation approach (Naseri et al. 2005) for prediction of crude oil viscosity. 3. Formulation The present study applies only to steady-state flow relevant to Newtonian fluids. By steady state flow it is understood that the flow conditions at each point in the stream do not change with time. For this approximation, the momentum balance equation is (Bird et al. 2002): 0 in out j j dp dp F dt dt    (5) where pin and pout are the momentum entering into and going out of the system respectively. ΣFj is the sum of the forces acting on the system. As shown in Figure 4, we focus our attention on a region of length L, http://www.springerlink.com/content/7581u35430540303/fulltext.html#Fig3 http://www.springerlink.com/content/7581u35430540303/fulltext.html#CR5 http://www.springerlink.com/content/7581u35430540303/fulltext.html#CR7 DYNAMICAL PROPERTIES OF OMANI CRUDE OILS 301 sufficiently far from the ends of the wall so that the entrance and exit disturbances are not included in L. This ensures that in this region the velocity component Vz does not depend on z. Figure 3. Kinematic viscosity versus temperature for the: (3a) Erad sample and (3b) Oman Export, Receive Line, Mabruk, and Zal-41 samples. Zal -41 Sample Erad Sample (low API) SAYYADUL ARAFIN and S.M. MUJIBUR RAHMAN. 301 3.1 Flow through a vertical annulus The flow of fluids in an annulus (Figure 4) is encountered frequently in physics, chemistry, biology, and engineering. The laminar flow of fluid in an annulus may be analyzed by means of momentum balance described in the previous section. By considering an incompressible fluid flowing in steady state in the annular region between two coaxial circular cylinders of radii κR and R (Figure 4), the equation for velocity profile of Bird et al. (2002) can be written with slight rearrangement in the form: Figure 4. Flow through an annulus. The various forces acting on the top and bottom surface of the annulus.   2 2 2 0 1 1 ln 4 ln L out z out out P P R r r v L R R                            (6) where κ is the ratio of the inner to the outer radius (κ = rin/Rout), Po is the atmospheric pressure, and r is the distance measured from the outer surface of the inner cylinder to the inner surface of the outer cylinder. The quantity PL represents the combined effect of pressure, pL and the gravitational term (ρgL). μ and ρ are viscosity and density of the fluid respectively. With the help of figure 4 it can be shown that the term (P0 - PL ) is the net pressure, ΔP which can be replaced by ΔP = Δp = pa- ρgL. pa is the applied pressure as shown in figure 4. Since ρ is a function of temperature, the net pressure, Δp, is a function of temperature as well. The fluid will remain static if the net pressure is equal to zero, which means that the applied pressure, pa at the bottom of the annulus will exactly balance ρgL. Since our prime objective is to determine the temperature dependence of flow DYNAMICAL PROPERTIES OF OMANI CRUDE OILS 301 properties we express all the relevant equations as a function of temperature through the temperature dependence of μ and ρ. Therefore Eq. 6 can be written in terms of net pressure as:   2 2 2 1 , 1 ln 4 ln out z out out p R r r v r T L R R                            (7) The second term in the square bracket may be termed as the pipe characteristic function, which at any point along r depends solely on the pipe dimension. Eq. 7 works only under the conditions that (i) the fluid is of constant density ρ (incompressible), (ii) the flow is laminar, (iii) the annulus length (L) is very large and (iv) there are no end effects. In fact, at the tube entrance and exit, the flow will not necessarily be parallel everywhere, so the tube surface effect will be ignored (Bird et al. 2002). Since, we need to know the maximum velocity, the average velocity and the shear stress at a surface we have to know the range of the applied pressure for the flow to be in laminar at a given radius, length, viscosit y and density. In this study we have assumed the flow to be laminar by restricting the Reynolds number, Re to 2000 and constraining Re by the equation (Bird et al. 2002):  Re 2 1out zR v        (8) where, z v  is the average velocity of the velocity profile, z v given by:       2 8 out z p T R v T B T L     (9) B may be termed the annular constant. It is obtained from the expression:    2 2 1 1 ln B               (10) Combining equations (8) and (9) it can be shown that:         2 3 4 Re 1 out TL p T R B T       (11) Because μ 2 (T) decreases more rapidly with increasing temperature than ρ does, Δp will therefore decrease with increasing temperature. Substituting Eq. 11 in Eq. 7 and Eq. 9, the expressions for z v and  z v in terms of Reynolds number, Re, reduce respectively to:          2 2 1Re , 1 ln 1 ln z out out out T r r v r T R T B R R                          (12) SAYYADUL ARAFIN and S.M. MUJIBUR RAHMAN. 330 with         Re 2 1 z out T v T T R        (13) Both z v and z v  will decrease with increasing temperature because the viscosity decreases more rapidly with increasing temperature than density. Although Re depends on temperature through ρ and μ, as suggested by Eq. 8, we have kept it fixed at 2000 for all temperatures. This will force the average velocity to adjust according to Eq. 8. By doing so no generality is lost with regard to the condition for laminar flow. This strategy makes the computation simple and less cumbersome. The mass rate of flow, w is given by (Bird et al. 2002):        2 21out zw T R T v T      (14) Substituting Eq. 13 in Eq. 14 one can get:      1 Re 2 out R T w T     (15) Again, the mass rate of flow, w will decrease with increasing temperature because of the presence of μ in the numerator of Eq. 14. By keeping Re fixed at 2000, the net pressure, Δp, z v , z v  and w have been calculated from Eq. 10, Eq. 11, Eq. 12 and Eq. 14 respectively as a function of temperature, T. The viscous force exerted by the fluid on the walls of the annulus, Fz, is given by (Bird et al. 2002):    2 2( ) 1z outF T R p T    (16) Since Δp decreases sharply with temperature (Eq. 11) at a fixed Re, Fz is expected to decrease with increasing temperature sharply as well. 3.2 Flow through a vertical cylinder The equation of the velocity profile, z v obtained for a laminar flow through a cylindrical tube (Bird et al. 2002) is expressed as:   22 ( ) 1 4 o L z P P R r v r L R             (17) where R is the radius of the cylinder. The variable r is the distance measured from the center to the inner wall of the cylinder. The temperature dependence of z v for a cylinder is:   22 ( , ) 1 ( ) 4 z p T R r v r T T L R                        (18) where     0 L ap P P p T gL     . The Reynolds number, Re for the flow in a cylinder is given (Bird et al. 2002) by the equation:       Re( ) 2 z T T R v T T      (19) DYNAMICAL PROPERTIES OF OMANI CRUDE OILS 333 with,       2 8 z p T R v T T L     (20) Figure 5. Flow in a vertical cylinder of length L and radius R. The various forces acting on the top and bottom faces of the cylinder are shown by arrows. As in the case of annular flow, we have kept the Reynolds number, Re, fixed at 2000 for the study of laminar flow in a cylinder as well. Following the same derivation procedure as in an annulus, one can show, in terms of Re, that the net pressure, Δp, average velocity, z v  and velocity profile, z v , for a cylinder can be written as:       2 3 4 Re TL p T R T     (21) ( ) Re ( ) ( ) 2 z T v T T R      (22) and 2 Re ( ) ( ) 1 ( ) z T r v T R T R              (23) And, finally the mass rate of flow,      4 ( ) 8 ( ) a p T gL R T w T T L       can be expressed as SAYYADUL ARAFIN and S.M. MUJIBUR RAHMAN. 331 Re ( ) ( ) 2 R w T T   (24) The viscous force exerted by the fluid on the wall of the cylinder, Fz, is given by (Bird et al. 2002):     2 z F T R p T  (25) As in the case of an annulus, Fz will decrease with increasing temperature in a similar fashion. 4. Results and discussions 4.1 Vertical annulus The calculated results on dynamical properties of the two extreme samples, namely Zal-41 [lightest] and Erad [heaviest] samples through an annulus are presented and discussed in this section. The outer radius (Rout) and inner radius (Rin) of the annulus are 0.5 m and 0.1 m respectively (rin/ Rout = κ = 0.2 ) and length, L=50 m. -0.6 -0.4 -0.2 0.0 0.2 0.4 0.6 -0.06 -0.04 -0.02 0.00 0.02 0.04 0.06 70 o C 70 o C 70 o C 70 o C 20 o C 20 o C 20 o C V e lo c it y p ro fi le , v z ( m s -1 ) Distance, r(m) 20 o C Zal-41 sample Re=2000 (a) -0.6 -0.4 -0.2 0.0 0.2 0.4 0.6 -2.0 -1.5 -1.0 -0.5 0.0 0.5 1.0 1.5 40 o C 40 o C 40 o C40 o C 30 o C30 o C 30 o C 30 o C 20 o C20 o C 20 o C V e lo c it y p r o fi le , v z ( m s - 1 ) Distance r(m) 20 o C Erad sample Re=2000 (b) Figure 6. The velocity profiles for the (6a) flow of light (Zal-41 sample) and (6b) heavy (Erad sample) crude oils through an annulus at various temperatures. Since the velocities [Eqs.7 and 9] within the bulk of the fluid depend on the viscosity ( )T , which decreases as the temperature increases, it seems that the profile of ( ) z v T as shown in Figure 6 is in contradiction with these equations. But the infused constraint owing to using a fixed Reynolds number Re = 2000 in Eq. 12 and Eq. 13, turns back the velocity profile; this is now reflected by the joint effect of )(/)( TT  which essentially decreases with temperature. DYNAMICAL PROPERTIES OF OMANI CRUDE OILS 331 20 30 40 50 60 70 0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 Re=100 Re=500 Re=1000 N e t p r e s s u r e ,  p ( P a ) Temperature, T( o c) Zal-41 Annulus Re=2000 (a) 20 30 40 50 60 70 0 200 400 600 800 1000 Re=100 Re=500 Re=1000 N e t p r e s s u r e ,  p ( P a ) Temperature, T o C) Re=2000 Erad sample Annulus (b) Figure 7. The net pressure as function of temperature at various Reynolds numbers, Re for (7a) Zal-41 sample and (7b) Erad sample. 20 30 40 50 60 70 0 2 4 6 8 10 12 14 16 18 20 Re=100 Re=500 Re=1000 M a s s r a te o f fl o w , w ( K g m s - 1 ) Temperature, T( o C) Zal-41 Annulus Re=2000 (a) 20 30 40 50 60 70 0 100 200 300 400 500 600 700 800 Re=100 Re=500 Re=1000 M a s s r a te o f fl o w , w ( K g m s - 1 ) Temperature, T( o C) Erad Annulus Re=2000 (b) Figure 8. The mass rate of flow w through an annulus as a function of temperature at various Reynolds number, Re, values for (8a) light (Zal-41) and (8b) heavy (Erad) crude oils. The profile of the curves in Figure7 can be explained in terms of Eq. 11 which clearly demonstrates that for a fixed Reynolds number Re, the net pressure ( )p T will decay because of the presence of the factor 2 ( ) / ( )T T  which decreases as T increases keeping also in mind that ( )T decreases more rapidly than ).(T The rapid increase of viscosity values at the limit of the low temperature range is also truly reflected Zal-41Sample Annulus Erad Sample Cylinder Zal-41 Annulus Erad Annulus SAYYADUL ARAFIN and S.M. MUJIBUR RAHMAN. 331 accordingly in the ( )p T curves in the figure. Concurrently it is also relevant to mention that the results ideally reflect the varying viscous nature of the lightest (Zal-41) and heaviest (Erad) samples through a rapid variation of the profile of ( )p T at this limit. The profiles of )(Tp and w(T) display some similar features; this is primarily because of the presence of the viscosity at the numerators in Eq. 11 and Eq. 15. 20 30 40 50 60 70 0.0 0.1 0.2 0.3 0.4 0.5 0.6 V is c o u s f o rc e o n t h e s o li d s u rf a c e , F z ( N ) Temperature, T( o C) Zal-41 sample (light) API=40.89 (a) Re=2000 20 30 40 50 60 70 0 100 200 300 400 500 600 700 800 V is c o u s f o rc e o n t h e s o li d s u rf a c e , F z ( N ) Temperature, T( o C) Erad sample (heavy) API=19.19 Re=2000 (b) Figure 9. The viscous force, Fz, on the walls of the annulus as a function of temperature at Reynolds number 2000 for: (9a) light (Zal-41) and (9b) heavy (Erad) crude oils. Owing to the presence of the net pressure ( )p T in the viscous force Fz, Figure 9 follows the profile of ( )p T as shown in figure 7. -0.10 -0.05 0.00 0.05 0.10 -0.30 -0.25 -0.20 -0.15 -0.10 -0.05 0.00 0.05 0.10 0.15 0.20 0.25 70 o C 70 o C 20 o C V e lo c it y p r o fi le , v z ( m s - 1 ) Distance, r(m) 20 o C Zal-41 sample Re=2000 (a) -0.10 -0.05 0.00 0.05 0.10 -8 -6 -4 -2 0 2 4 6 8 70 o C 70 o C 40 o C 40 o C 30 o C 30 o C 20 o C V e lo c it y p r o fi le , v z ( m s - 1 ) Distance r(m) Erad sample 20 o C Re=2000 (b) Figure 10. The velocity profile at various temperatures as a function of the distance, r from the center of the cylinder towards its surface for (10a) Zal-41 sample and (10b) Erad sample. Zal-41Sample (light) API = 40.89 Re = 2000 Erad Sample (heavy) API = 19.19 Re = 2000 DYNAMICAL PROPERTIES OF OMANI CRUDE OILS 331 4.2 Vertical cylinder The calculated results on dynamical properties of two extreme samples, namely Zal-41 [lightest] and Erad [heaviest] samples through cylindrical tubes are presented and discussed in this section. The length and radius of the cylinder are 30 m and 0.1m respectively. 20 30 40 50 60 70 0 5 10 15 20 25 30 Re=100 Re=500 Re=1000 N e t p r e s s u r e ,  p ( P a ) Temperature, T o C Zal-41 Cylinder Re=2000 (a) 20 30 40 50 60 70 0 5000 10000 15000 20000 25000 30000 35000 40000 Re=100 Re=500 Re=1000 N e t p r e s s u r e ,  p ( P a ) Temperature, T( o C) Re=2000 Erad sample Cylinder (b) Figure 11. The net pressure on the fluid samples in a cylinder as function of temperature at various Reynolds number Re values for: (11a) a Zal-41 sample and (11b) an Erad sample. Figures 10 (a-b) show plots of the net pressure values as a function of temperature for upward laminar flow (Re=2000) of Zal-41 and Erad crude oil samples through a cylindrical pipe. The pipe is 30 m long and has a radius of 10 cm. The net pressure is quite high at low temperatures and decreases sharply up to about 40 o C after which it slowly reduces to zero at higher temperatures. The heavy crude oil (Erad sample) which has low API value (19.19) shows a sharp decrease in net pressure in comparison with the light crude oil (Zal-41 sample) having a high API value of 40.89. The shape of the curves can be predicted from equation (11) which involves the viscosity term, μ 2 in its numerator indicating that this term decreases more rapidly than the density, ρ in the denominator. The zero net pressure means that the applied pressure, pa should equal the pressure (ρgL) due to the column of liquid in the cylinder. Since ρ decreases with increasing temperature, the applied pressure should decrease as well in order to balance the liquid column pressure for no flow through the pipe. pa in general will decrease with increasing temperature for a fixed Reynolds number. Owing to the presence of the net pressure ( )p T in the viscous force Fz, Figure 13 follows the profile of ( )p T as shown in Figure 11. 5. Conclusions We have investigated a few dynamical properties of some Omani crude oils through an annulus and a cylindrical pipe of fixed dimensions, in terms of the measured density and viscosity of these samples at various temperatures. The underlying methodology involves a simple model of fluid dynamics based on Newtonian approximation applied to laminar flow. Zal -41 Sample Cylinder Erad Sample Cylinder SAYYADUL ARAFIN and S.M. MUJIBUR RAHMAN. 331 The measured density and viscosity respectively follow the equations:  ( ) o r t m t t    and ( ) t E o t        within the selected range of temperatures. Subsequently we have studied the dynamical properties of these samples through an annulus and a finite-dimension cylindrical tube. 20 30 40 50 60 70 0.0 0.5 1.0 1.5 2.0 2.5 3.0 Re=100 Re=500 Re=1000 M a s s r a te o f fl o w , w ( K g s - 1 ) Temperature, T( o C) Re=2000 Zal-41 Cylinder (a) 20 30 40 50 60 70 0 20 40 60 80 100 120 Re=100 Re=500 Re=1000 M a s s r a te o f fl o w ,w ( k g m s - 1 ) Temperature, T( o C) Erad sample Cylinder Re=2000 (b) Figure 12. The mass rate of flow w through a cylindrical pipe as a function of temperature at various Reynolds number Re values for: (12a) a light Zal-41 sample and (12b) a heavy Erad sample. 20 30 40 50 60 70 0 5 10 15 20 25 V is c o u s f o r c e , F z ( N ) Temperature, T( o C) Zal-41 sample (light) API = 40.89 Re = 2000 (a) 20 30 40 50 60 70 0 5000 10000 15000 20000 25000 30000 V is c o u s f o r c e , F z ( N ) Temperature, T( o C) Erad sample(heavy) API=19.19 -500 (b) Figure 13. The viscous force Fz on the walls of the cylinder as a function of temperature at Reynolds number 2000 for: (13a) light (Zal-41) and (13b) heavy (Erad) crude oils. Zal -41 Sample Cylinder Erad Sample Cylinder Erad Sample (heavy) API = 19.19 Zal-41 Sample (light) API = 40.89 Re = 2000 DYNAMICAL PROPERTIES OF OMANI CRUDE OILS 331 The investigation of the dynamical properties using simplified models in conjunction with the measured viscosity data sheds some light on the Omani crude oil properties. However, we would like to make the following critical remarks in light of the presented findings: (a) The concept of momentum balance used in the formalism is merely an approximation applicable for ideal Newtonian liquid. Any non-Newtonian fluid has indeed a non-zero term on the right hand side of the momentum-balance equation. Consideration of non-Newtonian approximation will certainly necessitate the usage of the non-linear form of the Navier-Stokes equation. Thus it is apparent that the present approximation yields qualitative results only. (b) The dynamic properties reported in the article have been studied by using simplified models and geometry that have neglected wall and end effects, coarseness and non-linearity. These effects should be looked at using more realistic models and geometry within the realm of a non-Newtonian fluid model. (c) The Omani crude oils investigated have been treated as isotropic fluids without any inherent phase separation or segregation among components. These effects might be relevant to any investigation of dynamic properties. Presently we are working on a project considering the above points. 6. References ADECHY, D. and ISSA, R.I., 2004. Modeling of Annular Flow Through Pipes and T-junctions. Computers and Fluids, 33: 289-313. AHMED, T. 2000. Reservoir Engineering Handbook. Gulf Publishing, Houston,TX. ARAFIN, S., AL-HABSI, N. and RAHMAN, S.M.M., 2011. Transport Properties and Model-based Dynamical Properties of Omani Crude Oils. Arab J. Geosci. [Springer], DOI 10.1007/s/12517-011-0301-z. BAI, R., KELKAR, K., and JOSEPH, D.D., 1996. Direct Simulation of Interfacial Waves in a High-viscosity- Ratio and Axisymmetric Core-Annular Flow, J. Fluid Mech, 327: 1-34. BATZLE, M. and WANG, Z. 1992. Seismic Properties of Pore Fluids. Geophysics 57: 1396-1408. BIRD, R.B., STEWART, W.E. and LIGHTFOOT, E.N., 2002. Transport Phenomena. John Wiley & Sons, New York. GEORGE, A.K., ARAFIN, S., SINGH, R.N. and CARBONI, C., 2006. A Correlation Between Surface, Transport and Thermo-elastic Properties of Liquid Hydrocarbon: An Experimental Investigation. J. Phys: Condens. Matter 18: 3691-3701. HO, W.S. and Li, N.N. 1994. Core-annular Flow of Liquid Membrane Emulsion. AIChE J., 40: 1961-1968. LI H.-B., WU C. and ZHENG. Y-G., 2004. Plane Flow Model of Non-Newtonian Turbulent Stratified Flow in Wells and Pipes. J. Pet. Sci. & Eng. 44: 223-229. NASERI, A., NIKAZAR, M., MOUSAVI, S.A. and DEGHANI, A.A.M., 2005. A Correlation Approach for Prediction of Crude Oil Viscosities. J. Pet. Sci. & Eng. 47: 163-174. PRADA, J.W.V. and BANNWART, A.C., 2001. Modeling of Vertical Core-annular Flows and Application to Heavy Oil Production. Transaction of the ASME, 123: 194-199. SHERTOK, J.T. 1975. Velocity Profiles in Core-annular Flow using a Laser-Doppler Velocimeter. Ph.D. Thesis, Princeton University, Department of Chemical Engineering, Princeton, NJ. SHIOTSU, M. and HAMA, K., 2000. Film Boiling Heat Transfer from a Vertical Cylinder in Forced Flow of Liquids Under Saturated and Sub-cooled Conditions at Pressures. Nucl. Engg. Design, 200: 23-38. Received: 30 October 2011 Accepted: 20 November 2011