Analysis of Eyring-Powell Fluid in Helical Screw Rheometer

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.


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, 2 The Scientific World Journal 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][16][17][18][19][20][21][22][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.

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, V is the velocity field, f is the body force per unit mass, and T is the Cauchy stress tensor expressed as where denotes the dynamic pressure, I the unit tensor, and S denotes the extra stress tensor. The constitutive equation for Eyring-Powell fluid is given by [24] as where is viscosity, C, B are material constants with dimensions second −1 and Pascal −1 , respectively, |A 1 | = √(1/2) tr(A 2 1 ), and A 1 is the Rivlin-Ericksen tensor defined as where ∇V is the velocity gradient.

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 ( where = − sin , = cos .
Dropping " * " onward (16)-(17), give = 0, at = 0, 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.
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 The Scientific World Journal 5 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, (44)

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 The Scientific World Journal are constant coefficients.
The shear force per unit width required to move the upper plate in -direction is (75) equation (75) gives * = − * 1 ; The shear force per unit width required to move the upper plate in -direction is where * = / , * = / are dimensionless shear forces and 1 = 1 /ℎ, 2 = 2 /ℎ are dimensionless lengths and 1 and 2 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 * = * sin + * cos , where * = / is dimensionless shear force.
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.  (88)

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 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 = 45 ∘ , which confirms the results given in [25]. The resultant velocity given in (71), reduces to the velocity profile in -direction when we take = 0 ∘ . When we put = 90 ∘ the resultant velocity given in (71) recovers the velocity profile in -direction.

Conclusion
The steady flow of an isothermal, homogeneous, and incompressible Eyring-Powell fluid is investigated in HSR. Using ADM, the expressions for the velocity profiles are calculated. Expressions for the shear stresses, shear stresses at barrel surface, shear forces exerted on the fluid, volume flow rates, and average velocity of the fluid are also calculated. Graphical representation is given for the velocity profiles. It is observed that the velocity field depends on the involved parameters. The increase in the value of non-Newtonian parameters and pressure gradients increases the flow of the fluid. It can be seen that the shear thinning effect of is larger than in the fluid. It is also observed that the net velocity of the fluid is due 14 The Scientific World Journal to the pressure gradient. It is also noticed that the resultant velocity attains its maximum value at = 45 ∘ . Thus, the profound conclusion is that the extrusion process depends on the involved parameters.