# data[].attributes.headline
		
	
	
		
			
				
					Computational Methods For Multiphase flow

# data[].attributes.subject
		
	
	
		
			
				
					Computational Methods For Multiphase flow

# data[].attributes.body
Compilation of research papers by Cambridge University

(Downloaded)

# data[].attributes.headline
		
	
	
		
			
				
					The Journal of

					Computational Multiphase Flows 
				
			
		

# data[].attributes.subject
		
	
	
		
			
				
					The Journal of

					Computational Multiphase Flows 
				
			
		

# data[].attributes.body
Research paper by Yu-Shu Wu, Bitao Lai and Jennifer L. Miskimins

(Downloaded)

# data[].attributes.headline
FULL RESEARCH PAPERS AND MATERIALS ON MULTIPHASE FLOW

# data[].attributes.subject
FULL RESEARCH PAPERS AND MATERIALS ON MULTIPHASE FLOW

# data[].attributes.headline
		
	
	
		
			
				
					CFD Predictions of Drilling Fluid Velocity and Pressure Profiles in Laminar Helical Flow 
				
			
		

# data[].attributes.subject
		
	
	
		
			
				
					CFD Predictions of Drilling Fluid Velocity and Pressure Profiles in Laminar Helical Flow 
				
			
		

# data[].attributes.body
Research paper by F. A. R. Pereira, M. A. S. Barrozo and C. H. Ataíde

(Downloaded)

# data[].attributes.headline
		
	
	
		
			
				
					CFD Calculations of Cuttings Transport through
Drilling Annuli at Various Angles 
				
			
		

# data[].attributes.subject
		
	
	
		
			
				
					CFD Calculations of Cuttings Transport through
Drilling Annuli at Various Angles 
				
			
		

# data[].attributes.body
Research paper by Uduak Mme and Pål Skalle

(Downloaded)

# data[].attributes.headline
FULL RESEARCH PAPERS RELEVANT TO DRILLING FLUIDS

# data[].attributes.subject
FULL RESEARCH PAPERS RELEVANT TO DRILLING FLUIDS

# data[].attributes.headline
PAPERS RELEVANT TO WELLBORES (ABSTRACTS ONLY)

# data[].attributes.subject
PAPERS RELEVANT TO WELLBORES (ABSTRACTS ONLY)

# data[].attributes.headline
Heat Transfer and dissipation effects in the flow of a drilling fluid

# data[].attributes.subject
Heat Transfer and dissipation effects in the flow of a drilling fluid

# data[].attributes.body
Research paper by Wei-Tao Wu and Mehrdad Massoudi
(Downloaded)

# data[].attributes.attachment_link.readable_text
Heat Transfer and Dissipation Effects in the Flow of a Drilling Fluid

Wei-Tao Wu 1 and Mehrdad Massoudi 2 ,*

Department of Biomedical Engineering, Carnegie Mellon University, Pittsburgh, PA 15213, USA

U.S. Department of Energy, National Energy Technology Laboratory (NETL), 626 Cochrans Mill Road, P.O. Box 10940, Pittsburgh, PA 15236, USA

Correspondence: Tel.: +1-412-386-4975; Fax: +1-412-386-5870

Academic Editors: Goodarz Ahmadi and Bekir S. Yilbas

Received: 29 December 2015 / Accepted: 17 March 2016 / Published: 18 March 2016

drilling mud; drilling fluid; non-Newtonian fluids; suspension; rheology; variable viscosity; viscous dissipation

1. Introduction

Flow of mixtures composed of solid particles of various shapes and sizes dispersed in a fluid occurs in many industrial processes as well as in nature (see [ 1, 2]). Examples of such complex fluid mixtures are coal slurries and drilling fluids. Drilling fluids are generally water-based or oil/polymer-based. In the formulation of either the water-based or the oil-based drilling fluids several additives are added to maintain proper rheological properties. These additives, and the fact that drilling fluids are composed of mud, pieces of rocks, water, oil, bubble, etc., make the drilling fluid behave as complex non-Newtonian fluids. Flow of these fluids becomes even more complicated during the drilling operation as continuously drill cuts of different size, shape, concentration, pH, etc., are added into the drilling fluids. Amongst the most important properties of the drilling fluid one can name density, viscosity, yield stress, etc. When the drilling fluid is subjected to high temperature and high pressure then ionic mobility, reactivity, solubility etc., are increased and these, in turn, change the rheological properties of the drilling fluid. The drilling operation is quite often, as a temporary solution, suspended in order to change the drill bit, to clean the pipe or whenever some other difficulties are encountered. During this stage, the drilling fluid should have the capability of holding the rock cuttings in suspended form and prevent them from settling down in the well annulus (see [ 3]). Our main interest in this paper is to study the behavior of non-linear fluids which in many ways resemble and behave similarly to drilling fluids.

The success of the drilling operation depends, to a large extent, on variables such as the drill bit design, the selection of the drilling fluid, the climate conditions, the geophysical properties of the ground rock material, etc. [ 4]. The role of the drilling fluid is considered by many researchers to be a prominent one [ 5]. Hole cleaning is an important issue in the drilling operation, especially in cases of horizontal drilling, as mentioned by Caenn, et al. [ 6]. Siginer and Bakhtiyavou [ 7] indicate that the effect of the eccentricity of pipe on the flow of well bore fluids should also be studied. For deep drilling applications, in addition to external factors such as temperature and pressure, factors such as the permeability and deformation properties of the rocks and mudrock instability should be considered [ 8]. Hole cleaning is related to the ability of the drilling fluid to transport and suspend drilled cuttings. Hole cleaning depends on many variables including well bore inclination, cuttings slip velocity, flow regime, rotary speed of the drill pipe, fluid rheology, flow rate, rate of penetration, cuttings size and shape, wellbore geometry, cuttings density, cuttings agglomeration, the fluid density, the wellbore geometry, pipe rotation, etc. [ 9].

In most drilling applications, the Bingham or the Power-law models are used to represent the drilling fluid behavior (see [ 10]). However, many drilling fluids behave according to the modified Power-law or the Herschel-Buckley rheological model (see [ 10]). The shear rate is low in the annulus, and as it is well known, in the low shear-rate regime, that the Power-law model underestimates while the Bingham plastic model overestimates the frictional pressure drops. The Herschel-Buckley model perhaps presents a more adequate constitutive relation [ 11]. Slawomriski [ 12], indicated that even though the Bingham plastic model is often used in modeling the drilling fluids, some studies have shown that many of the drilling fluids are non-linear fluids with memory. Briscoe et al. [ 13] indicate that in a high pressure and high temperature environment, the yield stress and the plastic viscosity of the muds are also influenced by pressure and temperature. In many ways, a drilling fluid, when modeled as a (single component) suspension, behaves as a non-linear fluid similar to a slurry where the effects of shear rate and concentration are to be included in the constitutive relation for the stress tensor.

The objective of the present paper is to numerically study the effects of viscous dissipation in the Couette flow and heat transfer in such a fluid, and explore the effects of concentration and the shear rate and temperature-dependent viscosity, along with a variable thermal conductivity. We assume that the viscosity of this non-linear fluid depends on the shear rate and the volume fraction of the particles. The effects of the particle concentration (volume fraction) and its behavior is governed by a convection-diffusion equation. In the next section, the governing equations of motion are provided. Section 3 focuses on the constitutive relations for the stress tensor, the diffusive particle flux vector, and the heat flux vector. In Section 4, we describe the geometry of the problem and provide the derivation for the one-dimensional form of the governing equations, as well as the boundary conditions. In Section 5, we outline the numerical scheme we have used. In Section 6, we perform a parametric study using the different dimensionless numbers.

2. Governing Equations

If the drilling fluid is treated as a single component material, in the absence of any electro-magnetic effects, the governing equations of motion are the conservations of mass, linear momentum, angular momentum, concentration and the energy equation. If the drilling fluid is modeled as a multi-component material, then we need to provide the governing equations for all the components, and a multi-phase approach should be used; this requires not only constitutive relations for each component, but also for the interactions among the components [ 2, 14, 15]. In this paper, we assume that the drilling fluid can be treated as a single component non-homogeneous fluid. As a result, the governing equations are:

Conservation of mass

∂ρ ∂t+d i v ( ρ v)=0

where ρ is the density of the fluid, ∂/∂t is the partial derivative with respect to time, and v is the velocity vector. For an isochoric motion we have d i v v=0.

Conservation of linear momentum

ρ d v d t=d i v T+ρ b

where b is the body force vector, T is the Cauchy stress tensor, and d/d t is the total time derivative, given by d (.)/d t=∂ (.)/∂t+ [ g r a d (.)]v. The conservation of angular momentum indicates that in the absence of couple stresses the stress tensor is symmetric, that is, T= T T .

In suspension flows, the particle concentration is, in general, not constant and in many applications, it is necessary to have an additional equation, often called the convection-diffusion equation. Here we use the particle concentration equation for ϕ as discussed in [ 16], based on [ 17].

Conservation of concentration

∂ϕ ∂t+ v i ∂ϕ ∂ x i=−d i v N

The first term on the left-hand side denotes the rate of accumulation of particles, the second term denotes the convected particle flux (where ∂ϕ ∂ x i denotes the gradient of the concentration), and the term on the right side denotes diffusive particle flux. Following [ 17], the diffusive particle flux N is composed of fluxes related to the Brownian motion, the variation of interaction frequency and the viscosity.

Conservation of Energy:

ρ d ε d t=T:L−d i v q+ρ r 1

where ε is the specific internal energy, L is the gradient of velocity, q is the heat flux vector, r 1 is the specific radiant energy, and " :" designates the scalar product of two tensors. Thermodynamical considerations require the application of the second law of thermodynamics or the entropy inequality. The local form of the entropy inequality is given by (see [ 18], p. 130):

ρ η˙+d i v φ−ρ s≥0

where η ( x,t) is the specific entropy density, the dot implies the material time derivative, φ ( x,t) is the entropy flux, and s is the entropy supply density. If it is assumed that φ= q θ, and s= r θ, where θ is the absolute temperature, then Equation (5) reduces to the Clausius-Duhem inequality:

ρ η˙+d i v q θ−ρ r θ≥0

In this paper, we do not consider the consequences of the Clausius-Duhem inequality [ 18, 19, 20, 21]. In the next section, we will address the constitutive relations for the stress tensor, the heat flux vector q and the diffusive particle flux N .

3. Constitutive Equations

Drilling fluids are complex multi-component fluids; in general, a drilling fluid is composed of water, oil, pieces of rocks, sand, mud, etc., with bubbles and other chemicals added to this mixture for various reasons. Most researchers assume that a drilling fluid is a suspension whose rheological properties can be described using a non-linear (non-Newtonian) fluid. Based on the available experimental observations, drilling fluids exhibit characteristics similar to those of non-linear materials, such as colloidal suspensions, polymers, rubber, slag, etc. For complex materials, the main points of departure from linear behavior are: (1) The ability to shear thin or shear thicken; (2) The ability to creep; (3) The ability to relax stresses; (4) The presence of normal stress differences in simple shear flows; (5) The presence of yield stress; (6) Memory effects, etc. For many practical fluid engineering applications, the viscosity may be a function of time, shear rate, concentration, temperature, pressure, electric field, magnetic field, etc. Therefore, in general, μ=μ ( t,Π,θ,ϕ,p,E,B), where t is the time, Π is some measure of the shear rate, θ is temperature, ϕ is concentration, p is pressure, E is the electric field, and B is the magnetic field. For some materials or some applications, the dependence of one or more of these can be dropped (see [ 22]).

3.1. Stress Tensor

A general constitutive relation for stress tensor of the drilling fluid should at the very least be able to predict (or include) some type of yield stress and a viscous stress with shear-thinning capabilities, i.e., where the coefficient of viscosity depends on the shear rate. Thus, we assume:

where is the yield stress and is the viscous stress. In general, the yield stress can be obtained from experiments and for the viscous stresses a model is needed where the shear viscosity depends on volume fraction, temperature, pressure, chemical composition and the shear rate. In this paper, we will not include the yield stress (In our next study, we plan to model and study the yield stress part of the stress tensor.) in our model and we only consider the viscous stress which we assume to be given by a generalized Power-law fluid model:

where is the reference viscosity, is the maximum volume fraction at which the suspension exhibits fluid behavior, β is an experimentally determined coefficient, and are reference values, and α is a constant. Here, 1 is the identity tensor, p is the pressure, θ is temperature, ϕ is the volume fraction (The function ϕ is an independent kinematical variable called the volume distribution or volume fraction function (related to concentration) having the property 0≤ϕ ( x,t)≤<1. The function ϕ is represented as a continuous function of position and time; in reality, ϕ in such a system is either one or zero at any position and time, depending on whether one is pointing to a particle or to the void space (fluid) at that location. Now, ρ , is related to (density of pure fluid) and ϕ through ρ= ( 1−φ) ρ f.), t r is the trace operator, Π is an invariant of D where D is the symmetric part of the velocity gradient, μ and Π together form the effective viscosity, and m is a material parameter. When m < 0, the fluid is shear thinning, and if m > 0, the fluid is shear thickening. For the remainder of this paper, we drop the dependency of μ on p, while recognizing that in many drilling applications, especially in the case of deep-ocean drilling, the effect of pressure should be included (We should mention that the viscosity of drilling muds at high pressure and high temperatures, will also depend on pressure and temperature. For example, Briscoe, et al. [ 13] suggested a high shear-rate viscosity of the type, μ= μ o ( 1+ α 1 (ϕ)+ α 2 (p)p)exp ( E u+p ( B θ−C))/k θ),where is a factor related to the activation energy , the Boltzmann constant ( k ) etc., B and C are constants, ϕ is the volume fraction, and α 2 (p) is a compressibility term.). Note that in our previous study, Zhou et al. [ 3], we used a similar but simpler model where the effect of temperature on viscosity was ignored; also, no heat transfer was considered in that paper. In the current model, as shown in Equation (8) which is based on our previous studies, including Gupta and Massoudi [ 23], Miao et al. [ 24] and Miao and Massoudi [ 25], we have included an exponentially decaying function for the effects of the temperature. The basic form of the viscosity is based on Krieger's viscosity-concentration correlation ( μ= μ r ( 1− ϕ ϕ m) −1.82) [ 26]. In the model described in Equation (8), we have assumed that the apparent viscosity is not only a function of temperature but also of volume fraction, obeying the Einstein-Roscoe relation [ 27, 28] (if β=2.5). It should be mentioned that alternative ways of modeling non-homogeneous fluids exist (see Massoudi and Vaidya [ 29, 30]).

3.2. Heat Flux Vector

According to Fourier [ 31] (see also Winterton [ 32]) the heat flux vector depends on the temperature gradient, where

and k is the thermal conductivity of the material. For many complex materials, k is considered an effective or modified form of the thermal conductivity that depends on concentration, temperature, etc.. In fact, for anisotropic materials, k becomes a second order tensor. (For a discussion of the effective thermal conductivity concept in porous media and multiphase flows, see [ 33] (p. 129) and [ 34, 35, 36]). Jeffrey [ 37] derived an expression for the effective thermal conductivity which includes the second order effects in the volume fraction [ 38]:

where

where

where ω is the ratio of conductivity of the particle to that of the matrix, k the effective conductivity of the suspension, the conductivity of the matrix, and ϕ is the solid volume fraction [ 34]. More recently, Pabst [ 39] has derived a relationship for the effective thermal conductivity,

We will use Equations (12)-(16) in this paper. Miao et al. [ 40] used Equation (12) but in their study they assumed that the thermal conductivity is a constant. In this paper, we will consider the case where the thermal conductivity is a function of volume fraction.

3.3. Particle Flux

As noted in [ 16], the various particle fluxes in the concentration equation are due to different mechanisms such as the Brownian motion, sedimentation, particle interactions, etc. In this paper, we assume that the Brownian motion can be neglected and consequently the diffusive particle flux is modeled as:

where is the flux due to particle interactions and is the flux associated with spatial variations in the viscosity. Based on the proposal of [ 17], we assume:

where a is the characteristic particle length (e.g., radius), γ˙ is the local shear rate, given by Equation (21), μ is the effective viscosity, and K c, K μ are empirically determined coefficients.

In the next section, we will examine the Couette flow of a fluid modeled as a concentrated suspension where the constitutive relations are given by Equations (8), (13)-(16) and (18)-(20).

4. Couette Flow

For the fully developed Couette flow of a suspension, (shown in Figure 1), we assume:

(i)

(ii)

the particle flux due to the Brownian diffusion is neglected;

(iii)

the constitutive equation for the stress tensor is given by Equation (8), the constitutive relation for the particle flux is given by Equations (18)-(20), the constitutive relation for the heat flux is given by Equations (13)-(16).

(iv)

the velocity, the volume fraction, and the temperature profiles are of the form:

{ v=v (r) e θ ϕ=ϕ (r) θ=θ (r)

Also, the specific radiant energy, r 1, is assumed to be negligible. Then the energy equation becomes

μ r ( 1− ϕ ) −β e α ( θ 0−θ) | d v d r| m+2− 1 r 2 d d r ( r 2 k m ( 1+3 ξ ϕ+ ξ^ ϕ 2) d θ d r)=0

where the first term on the left-hand side is the viscous dissipation (see [ 41]). Equations (23), (24) and (26) are subject to the following boundary conditions:

{ a t r= R i:v= R i Ω=V;θ= θ i a t r= R o:v=0;θ= θ 0

where and are the inner and outer radii, R i Ωis the inner wall velocity, and and are the temperature of the inner and the outer walls. Here we assume the no-slip condition for the velocity at the outer wall and a constant velocity at the inner wall; ϕ is the value of the volume fraction, and its value depends to some extent on the average volume fraction ϕ a v g through:

ϕ a v g= 2 R o 2 − R i 2 ∫ R i R o ϕ r d r

We need to assume a value for ϕ a v g so that the solution of ϕ can be determined accordingly. If the calculated volume fraction distribution satisfies condition Equation (28), then the original assumption is valid; otherwise a new guess for ϕ must be made until condition Equation (28) is satisfied.

Here Re is the Reynolds number, given by:

R e= ρ V R μ r; and M=α ( θ i− θ 0)

For the current problem, due to the kinematics of the flow, that is, the fact that the flow is steady and fully developed, we could have eliminated the Re. However, we keep the form of the equation as it is for problems involving the pressure gradient term. Substituting Equation (29) into Equations (24) and (26) results in:

( ϕ 2 d d r¯+ ϕ d ϕ d r¯) + ϕ 2 1 ( 1− ϕ ϕ m) −β e −M θ¯ | d v¯ d r¯|m d d r¯ [ ( 1− ϕ ϕ m) −β e −M θ¯ | d v¯ d r¯|m]=0

R 4 ( 1− ϕ ) −β e −M θ¯ | d v¯ d r¯| m+2− 1 r¯2 d d r¯ ( r¯2 ( 1+3 ξ ϕ+ ξ^ ϕ 2) d θ¯ d r¯)=0

where R 4= μ r V 2 k m ( θ i− θ 0) ( V R o)m, and K c K μ is the ratio of the effect of varying interaction frequency to varying the viscosity on the non-uniform distribution of volume fraction. Substituting Equation (29) into Equations (27) and (28) gives:

{ a t r¯= R i R o: v¯=1; θ¯=1. a t r¯=1: v¯=0; θ¯=0;

ϕ a v g= 2 1− R i 2 ∫ R i R o 1 ϕ r d r

Now Equations (30), (32) and (33) can be solved numerically with the conditions and constraints of Equations (34) and (35). Note that R 4 is related to the Prandtl and Eckert numbers ( R 4=Pr E c).

5. Numerical Results

The dimensionless differential equations are solved using the MATLAB solver bvp4c, which is a collocation boundary value problem solver. The step size is automatically adjusted by the solver and the default relative tolerance for the maximum residue is 0.001. The constraint Equation (35) for ϕ a v g was achieved by applying the shooting method. Table 1 shows the designated values of the studied dimensionless numbers and parameters in current work.

5.1. Effect of m

Figure 2, Figure 3 and Figure 4 show the effect of m on the distribution of velocity, volume fraction and temperature. Recall that when m > 0 the fluid is shear thickening and when m < 0 the fluid is shear thinning. Therefore, we can see that as m increases, that is, as the fluid increasingly shear thickens, contributing to a more uniform apparent viscosity, the velocity profiles become more linear as shown in Figure 2. From Figure 3 we can see that the shear-dependent viscosity has a strong effect on the particle distribution. For the case when m = 1, the value of the volume fraction at the outer wall is about ten times that of the value of the volume fraction at the inner wall. For the range of parameters considered in this study we can see that the effect of m on the temperature distribution is very small.

5.2. Effect of

Recall that ϕ mrepresents the maximum packing of solid particles in a solid-fluid system. A higher maximum packing can be attributed to the non-uniform distribution of the particle sizes, the deformability of the particles and so on. From Figure 5 we can see that a smaller produces a more non-linear velocity profile. Figure 6 shows that as increases, the volume fraction distribution becomes more non-uniform. For the case of ϕ m=0.45, it can be seen that near the outer wall, the particle distribution is very uniform and the value of the volume fraction is close to 0.45. Furthermore, in that region, the fluid motion is negligible (see Figure 5). In other words, the particles seem to accumulate in the region near the outer wall. Similar to the parameter m, the effect of on temperature distribution is negligible (see Figure 7).

5.3. Effect of M

In the current model, M is the parameter controlling the effect of temperature on the apparent viscosity. A larger M represents a higher sensitivity of the viscosity on the temperature change. Therefore, from Figure 8, we can see that as M increases, the velocity curve becomes more non-linear, which can be attributed to a more non-linear behavior of the apparent viscosity. It should be emphasized that if we had considered a pressure-driven flow, then as M increases we would expect a higher flow rate, since with the same temperature, a larger M implies a lower apparent viscosity. Figure 8 shows that at a larger M, such as M = 5, a higher velocity is reached mainly in the region near the inner wall. From Figure 9 and Figure 10, we can see that the effect of M on the volume fraction and temperature distribution is minimal.

5.4. Effect of

Recall that is defined as μ r V 2 k m ( θ 1− θ 0) ( V H)m, which is a measure of viscous dissipation and is related to the Prandtl and the Eckert numbers. Figure 11, Figure 12 and Figure 13 show the impact of in this problem. From Figure 11 and Figure 12, it can be seen that the effect of on the velocity and volume fraction distribution is negligible. From Figure 13 it can be seen that increases if the temperature profiles become more non-linear. From the expression of , we know that a larger implies a stronger effect of the shear stress on temperature distribution; furthermore, from the velocity profile, we can see that near the inner wall, the velocity gradient is larger and therefore a larger has the tendency to increase the temperature gradient near the inner wall region.

5.5. Effect of ω

Recall that ω represents the ratio of the conductivity of the particles to the matrix; ω only appears in the energy equation. Therefore, as expected from Figure 14 and Figure 15, we can see that the effect of ω on the velocity and the volume fraction distribution is meagre. Figure 16 shows the effect of ω on the temperature profile. It can be seen that as ω increases, the temperature distribution becomes more non-linear. This may be attributed to the higher dependency of the thermal conductivity on the volume fraction as ω increases.

5.6. Effect of

Figure 17, Figure 18 and Figure 19 show the effect of . This ratio represents the effect of varying interaction frequency to varying viscosity on the non-uniform distribution of the volume fraction. From Figure 18, we can see that as decreases, the particle distribution becomes more uniform. That is, when the is small, a small particle non-uniformity can cause viscosity gradients large enough to satisfy Equation (32) [ 17]. Figure 17 and Figure 19 show that as increases, the velocity and the temperature curves become more linear, though the effect is small.

5.7. Effect of the Average Volume Fraction

From Figure 20, it can be seen that when increases, the velocity curves become more non-linear. This may be attributed to higher non-uniform distribution of volume fraction when is high. As Figure 21 shows, the effect of on the volume fraction distribution is strong. For higher values of , the volume fraction profile becomes more non-linear; that is, more particles stack near the outer wall. At the same time, we can see that the velocity does not change much in the region where the particles accumulate. That is as expected, since when the bulk volume fraction is approaching the maximum packing, then a very high apparent viscosity is reached. From Figure 22, we can see that the effect of on temperature profile is negligible.

6. Concluding Remarks

In this paper, we have studied the fully developed flow of a drilling fluid between two rotating cylinders. We have modeled the drilling fluid as a non-Newtonian suspension. The viscosity is assumed to depend not only on the volume fraction and the shear rate but also on the temperature. To study the variation of the volume fraction, a convection-diffusion equation with a concentration flux suggested by Phillips et al. [ 17] is used. The condition of the fully developed flow is considered and the governing equations are non-dimensionalized and numerically solved. The parametric study indicates that the flow field is greatly affected by the dependency of the viscosity on the shear rate and the temperature, and that the concentration flux parameters, K c/ K μ, and the solid-particle maximum packing, , play a significant role in the particle distribution. For future studies, we plan to consider the yield stress portion of the stress tensor, as well as look at unsteady multi-dimensional flow situations. Furthermore, the effects of viscous dissipation near the tip of the drilling bit is considered to be more significant than inside the flow domain, and therefore, we plan to study the torsional and longitudinal oscillation of a cylindrical rod immersed in a non-linear fluid.

Author Contributions

Both authors contributed to the writing of the paper. Most of the mathematical formulation was primarily done by Mehrdad Massoudi, part of the derivation and all of the numerical work were done by Wei-Tao Wu.

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

Symbol

Explanation

b
body force vector

ϕ

concentration

D
symmetric part of the velocity gradient

g

acceleration due to gravity

l
identity tensor

L
gradient of the velocity vector

m

Power-law index

t

time

T
Cauchy stress tensor

ρ

bulk density

μ r

reference viscosity

η

effective viscosity

θ

temperature

div

divergence operator

∇ (or grad)

gradient operator

References

1. Joseph, D.; Renardy, Y. Fundamentals of two-fluid dynamics. J. Fluid Mech. 1995, 282, 405-406. [Google Scholar]

2. Rajagopal, K.; Tao, L. Mechanics of Mixtures; Series on Advances in Mathematics for Applied Sciences; World Scientific: Singapore, 1995; Volume 35. [Google Scholar]

3. Gray, G.R.; Darley, H. Composition and Properties of Oilwell Drilling Fluids; Gulf Publishing Company: London, UK, 1981. [Google Scholar]

4. Rabia, H. Oilwell Drilling Engineering: Principles and Practice; Graham and Trotman Ltd: London, UK, 1985. [Google Scholar]

5. Valluri, S.G.; Miska, S.Z.; Yu, M.; Ahmed, R.M.; Takach, N. Experimental study of effective hole cleaning using "sweeps" in horizontal wellbores. In Proceedings of the SPE Annual Technical Conference and Exhibition, San Antonio, TX, USA, 24-27 September 2006.

6. Coussot, P. Mudflow Rheology and Dynamics; Balkema: Amsterdam, The Netherlands, 1997. [Google Scholar]

7. Yield-power Law Model More Accurately Predicts Mud Rheology. Available online: http://www.ogj.com/articles/print/volume-91/issue-34/in-this-issue/drilling/yield-power-law-model-more-accurately-predicts-mud-rheology.html (accessed on 17 March 2016).

8. Liu, I.-S. Continuum Mechanics; Springer: Berlin, Germany, 2002. [Google Scholar]

9. Ziegler, H. An Introduction to Thermomechanics; North-Holland Publishing Company: Amsterdam, The Netherlands, 1983. [Google Scholar]

10. Truesdell, C.; Noll, W. The Nonlinear Field Theories of Mechanics; Springer-Verlag: Berlin, Geramny, 1992. [Google Scholar]

11. Roscoe, R. Flow Properties of Disperse Systems; Hermans, J., Ed.; Interscience Publishers: New York, NY, USA, 1953; pp. 1-38. [Google Scholar]

12. Fourier, J. The Analytical Theory of Heat; Dover Publications: New York, NY, USA, 1955. [Google Scholar]

13. Kaviany, M. Principles of Heat Transfer in Porous Media; Springer: Berlin, Germany, 1995. [Google Scholar]

Table 1. Designated values of the dimensionless numbers and parameters.

m

ϕ m

M

R 4

−0.5, 0.0, 0.5, 1.0

0.45, 0.5, 0.68, 0.9

0.0, 0.5, 1.0, 5.0

0.1, 0.5, 2.0, 3.0

ω

K c/ K μ

ϕ a v g

--

0.01, 1.0, 100.0, 1000.0

0.4, 0.6, 0.8, 0.9

0.1, 0.3, 0.4, 0.6

--

# data[].attributes.attachment_link.title
Heat Transfer and Dissipation Effects in the Flow of a Drilling Fluid

# data[].attributes.attachment_link.description
In this paper we study the effects of dissipation in the Couette flow and heat transfer in a drilling fluid, and explore the effects of concentration and the shear-rate and temperature-dependent viscosity, along with a variable thermal conductivity. A brief discussion on the constitutive relations for the stress tensor, the diffusive particle flux vector, and the heat flux vector is presented.

# data[].attributes.headline
Study of Rock-Bit Interaction and Wellbore Deviation

# data[].attributes.subject
Study of Rock-Bit Interaction and Wellbore Deviation

# data[].attributes.headline
Wellbore Temperature during Mud Circulation and Cementing 

# data[].attributes.subject
Wellbore Temperature during Mud Circulation and Cementing 

# data[].attributes.headline
Wellbore Pressure and Temperature during Gas-Well Testing

# data[].attributes.subject
Wellbore Pressure and Temperature during Gas-Well Testing

# data[].attributes.headline
The Effect of Wellbore Curvature on Tubular Buckling and Lockup

# data[].attributes.subject
The Effect of Wellbore Curvature on Tubular Buckling and Lockup

# data[].attributes.headline
Temp. Variation in a Circulating Wellbore Fluid

# data[].attributes.subject
Temp. Variation in a Circulating Wellbore Fluid

# data[].attributes.headline
Microannulus Cracks Developed Around a Wellbore

# data[].attributes.subject
Microannulus Cracks Developed Around a Wellbore