Research Article | Open Access

A. M. Siddiqui, T. Haroon, M. Zeb, "Analysis of Eyring-Powell Fluid in Helical Screw Rheometer", *The Scientific World Journal*, vol. 2014, Article ID 143968, 14 pages, 2014. https://doi.org/10.1155/2014/143968

# Analysis of Eyring-Powell Fluid in Helical Screw Rheometer

**Academic Editor:**Y. Liu

#### Abstract

This paper aims to study the flow of an incompressible, isothermal Eyring-Powell fluid in a helical screw rheometer. The complicated geometry of the helical screw rheometer is simplified by “unwrapping or flattening” the channel, lands, and the outside rotating barrel, assuming the width of the channel is larger as compared to the depth. The developed second order nonlinear differential equations are solved by using Adomian decomposition method. Analytical expressions are obtained for the velocity profiles, shear stresses, shear at wall, force exerted on fluid, volume flow rates, and average velocity. The effect of non-Newtonian parameters, pressure gradients, and flight angle on the velocity profiles is noticed with the help of graphical representation. The observation confirmed the vital role of involved parameters during the extrusion process.

#### 1. Introduction

The study of the rheological characteristics of different fluids is essential in the process of processing to obtain the desired quality and shape of the products. During processing, noticeable physical and chemical changes can occur [1]. For measuring the rheological properties of fluids in industries, mostly in the food industry, the available instruments are different types of viscometers. All these viscometers have their own advantages and limitations, they do not measure the fundamental physical parameters absolutely and they are empirical in nature [2, 3]. Kraynik et al. [4] have developed an instrument alternative to viscometers, called helical screw rheometer (HSR), for use in coal liquefaction processing that could characterize fluid suspensions accurately and consistently.

In order to optimize the processing and improve the quality of production, the literature on the classical extrusion theory contains the work of Carley et al. [5], Mohr and Mallouk [6], Booy [7], and Bird et al. [8].

Tamura et al. [9] have tried successfully the preceding analysis in the geometry of helical screw rheometer for Newtonian fluid and Power law fluid. In recent years, the study of non-Newtonian fluids have attracted many researchers. This is mostly due to their wide use in the food industry, chemical process industry, construction engineering, power engineering, petroleum production, and commercial and technological applications. Examples of non-Newtonian fluids are industrial materials, such as polymer melts, paints, gels, rubbers, soaps, inks, oils, concrete, ketchup, pastes, suspensions, slurries, and biological liquids such as blood and foodstuffs. The rheological knowledge of such fluids is of special importance owing to its application to many industrial problems. Considerable efforts have been made towards understanding their flows. Different models are available to characterize the non-Newtonian behavior of fluids out of which the fluid based on Eyring-Powell model is chosen for simplicity in the present work.

The basic governing equations for non-Newtonian fluids motion are highly nonlinear differential equations having no general solution, and only a limited number of exact solutions have been established for particular problems. Therefore, these problems should be treated by using some numerical or analytical methods. The analytical study of such type of nonlinear problems is important not only because of their technological significance, but also due to the interesting mathematical features presented by the governing differential equations of the flow. Apart from numerical methods, several analytical techniques such as the regular perturbation technique [10], the homotopy analysis method (HAM), the homotopy perturbation method (HPM) [11, 12], the variational iteration method (VIM) [13, 14], and the Adomian decomposition method (ADM) are mostly in use to overcome nonlinearity and get solutions. In this paper, we aim to apply the iterative technique ADM, which was introduced and developed by George Adomian and well addressed in the literature. ADM has recently received ample attention in the area of series solutions. A considerable amount of research work has been invested in the application of this method to a wide class of linear, nonlinear, partial differential equations, and integral equations [15–23].

This work considered the steady flow of an incompressible, isothermal, and homogeneous Eyring-Powell fluid in helical screw rheometer (HSR). Using Adomian decomposition method (ADM), analytical solutions are obtained for the governing equations in the geometry under consideration. Expressions of the velocity profiles, shear stresses, shear stresses at wall, forces exerted on fluid, volume flow rates, and average velocity are also calculated. The effect of involved dimensionless parameters on flow profiles are investigated through graphs and are discussed.

The paper is organized as follows. Section 2 contains the basic equations governing the motion of the fluid. In Section 3 the problem under consideration is formulated. Section 4 is devoted to the analytical solutions of the flow profiles, shear stresses, shear stresses at wall, forces exerted on fluid, volume flow rates, and average velocity. Section 5 is related to the discussion about the effect of the involved parameters on the motion of the fluid. Appropriate conclusions are drawn in Section 6.

#### 2. Basic Equations

The equations of conservation of mass and momentum for an incompressible fluid are where is the constant fluid density, denotes the material time derivative, is the velocity field, is the body force per unit mass, and is the Cauchy stress tensor expressed as where denotes the dynamic pressure, the unit tensor, and denotes the extra stress tensor.

The constitutive equation for Eyring-Powell fluid is given by [24] as where is viscosity, , are material constants with dimensions and , respectively, , and is the Rivlin-Ericksen tensor defined as where is the velocity gradient.

#### 3. Problem Formulation

Consider the steady flow of an isothermal, incompressible, and homogeneous Eyring-Powell fluid in helical screw rheometer (HSR). The complicated geometry of HSR is simplified in such a way that the curvature of the screw channel is ignored, unrolled, and laid out on a flat surface. The barrel surface is also flattened. Assume that the screw surface, the lower plate, is stationary and the barrel surface, the upper plate, is moving across the top of the channel with velocity at an angle to the direction of the channel (Figure 1). The phenomenon is same as the barrel held stationary and the screw rotates. The geometry is approximated as a shallow infinite channel by assuming that the width of the channel is large compared with the depth ; edge effects in the fluid at the land are ignored. The coordinate axes are positioned in such a way that the -axis is perpendicular to the flight walls, -axis is normal to the barrel surface, and -axis is in down channel direction. The liquid wets all the surfaces and moves by the shear stresses produced by the relative movement of the barrel and channel. No leakage of the fluid occurs across the flights. For simplicity, the velocity of the barrel relative to the channel is decomposed into two components (see Figure 1): along -axis and along -axis [6]. The associated boundary conditions can be taken as seen below (Figure 1) where The geometry of the problem suggests that the velocity profile and extra stress tensor are

Using (8) in (5) and (4), we get the nonzero components of the extra stress tensor given as

Using (8), (1) is identically satisfied and the substitution of the velocity profile of (8) and (9) in (2), in the absence of body forces results in

Equation (11) shows that , since the right sides of (10) and (12) are functions of alone and ; this implies and . Maclaurin series expansion of the inverse sine hyperbolic function in (10) and (12), when neglecting higher powers as , is

By simplifying (13), we obtained

Introducing the dimensionless parameters in (6) and (14), we get where and are dimensionless non-Newtonian parameters. Equation (6), in dimensionless form reduces to

Dropping “*” onward (16)-(17), give

Equations (18) and (19) are second order nonlinear ordinary differential equations, with boundary conditions (20); the exact solutions seem to be difficult. In the following section we use the Adomian decomposition method to obtain the approximate solutions.

#### 4. Solution of the Problem

The Adomian decomposition method (ADM) [15–23] describes that the operator form (18) and (19) can be written as in where is the differential operator taken as the highest order derivative to avoid difficult integrations. Since is invertible, this implies that exists.

On applying to both sides of (21), we obtained where , , , and are arbitrary constants of integration and can be determined by using boundary conditions. According to the procedure of ADM, and can be written in component form as

Using (23) in (22) results in Since the nonlinear terms can be explored in the form of the Adomian polynomials, say and , which yield (24), in the form and the boundary conditions (20) will take the form

From the recursive relation in (26)–(29), we can identify the zeroth order problems as with boundary conditions

The remaining order problems are in the following form: with the boundary conditions

From (25) we can calculate the components of the Adomian polynomials, say and , as the remaining components of the Adomian polynomials can be generated easily.

The ADM solutions to (26) and (27) with the boundary conditions (28) and (29) will be the sum of all order solutions; that is,

##### 4.1. Zeroth Order Solution

Zeroth order solutions of (18)–(20) can be calculated from the relations given in (30)–(33), which are equations (45) and (46) give the solutions for the linearly viscous fluid by assuming .

##### 4.2. First Order Solution

Equations (34)–(37) give the first order problems as with the boundary conditions

Using (38) and (41) in (47)–(50), we get where are constant coefficients.

##### 4.3. Second Order Solution

Equations (34)–(37) give the second order problems as with the boundary conditions

Using (39) and (42) in (54)–(57), we obtained where are constant coefficients in (58) and (59).

##### 4.4. Third Order Solution

Equations (34)–(37) give the third order problems as with the boundary conditions

Using (40) and (43) in (61)–(64), we get where are constant coefficients in (65) and (66).

##### 4.5. Velocity Profiles

###### 4.5.1. Velocity Profile in -Direction

Equations (45), (51), (58), and (65) give the ADM solution for the velocity profile up to order three in the transverse plane:

###### 4.5.2. Velocity Profile in -Direction

Equations (46), (52), (59), and (66) give the ADM solution for the velocity profile up to order three in the down channel direction:

###### 4.5.3. Velocity in the Direction of the Axis of Screw

The velocity in the direction of the axis of the screw at any depth in the channel can be computed from (68) and (69) as Equation (71) represents the resultant velocity of the flow and shows that the forward velocity at a point in the channel only depends on the pressure gradients and .

##### 4.6. Shear Stresses

Using (68)-(69) in (9) we obtain where is the nondimensional shear stress.

The shears exerted by the fluid on the wall at are where and are nondimensional shears at wall.

The shear force per unit width required to move the upper plate in -direction is equation (75) gives

The shear force per unit width required to move the upper plate in -direction is equation (77) gives where , are dimensionless shear forces and , are dimensionless lengths and and are lengths of the channel in and -directions. The net shear force per unit width in the direction of the axis of the screw can be computed from (76) and (78) as where is dimensionless shear force.

##### 4.7. Volume Flow Rates

Volume flow rate in -direction per unit width is where is dimensionless volume flow rate in -direction; using (68), (80) gives

Volume flow rate in -direction per unit width is where is the dimensionless volume flow rate in -direction; using (69), (82) gives

Equation (71) gives the resultant volume flow rate forward in the screw channel, which is the product of the velocity and cross-sectional area integrated from the root of the screw to the barrel surface. Consider where is dimensionless volume flow rate in the direction of the axis of the screw and is the number of parallel flights in a multiflight screw. Using (71), (84) gives equation (85) can be written as

##### 4.8. Average Velocity

The average velocity in the direction of the axis of the screw is where is nondimensional average velocity. Using (71) in (87), we get

#### 5. Results and Discussion

In the present work we have considered the steady flow of an incompressible, isothermal, and homogeneous Eyring-Powell fluid in HSR (see Figure 1). By using ADM, solutions are obtained for velocity profiles in - and -directions and also in the direction of the axis of the screw . Expressions for the shear stresses ( and ), shear stresses at barrel surface, forces exerted on fluid, volume flow rates, and average velocity are also calculated. Here we discussed the effect of non-Newtonian parameters , flight angle , and pressure gradients and on the velocity profiles with the help of graphical representation. From Figures 2, 3, and 8, we can observe the behavior of velocity profiles against . Figure 2 is sketched for , back flow is seen toward the barrel surface after some points in the channel height which show that the fluid circulates inside the confined channel; thus the velocity in -direction helps in the process of mixing during processing. In Figure 3 we observe that with the increase in value of the velocity increases and helps to move the fluid in the forward direction in the channel. The resultant velocity is shown in Figure 8, which resembles the Poiseuille flow in the channel. Due to the fluid moves toward the die. It is worthwhile to note that the shear thinning occurs with the increase in value of . The velocity profiles for the Newtonian case are retrieved for [6]. Figures 4, 5, and 9 are plotted to notice the effect of on velocities , , and . It is seen that velocity profiles are in the same pattern for and . It is also noticed that both non-Newtonian parameters depict the shear thinning effects in the fluid. However, graphical representation shows that shear thinning effects of are larger than , as the increase in velocity profiles is observed to be larger for . It is noticed that during the extrusion process, thinning/thickening of the fluid can be controlled with the proper choice of and . Thus, both the parameters play a vital role in the process of processing.

Moreover, Figures 6 and 10 are sketched for the velocity profiles and for different values of . It is found that the rise in pressure gradient increases speed of flow. Figures 7 and 11 are plotted for the velocity profiles and for different values of . It is seen that the increase in the value of increases speed of flow. Figure 12 is plotted for different values of . It is observed that the resultant velocity attains its maximum value at , which confirms the results given in [25]. The resultant velocity given in (71), reduces to the velocity profile in -direction when we take . When we put the resultant velocity given in (71) recovers the velocity profile in -direction.