Advances in Mathematical Physics

Volume 2017 (2017), Article ID 9724381, 10 pages

https://doi.org/10.1155/2017/9724381

## Unsteady Helical Flows of a Size-Dependent Couple-Stress Fluid

^{1}Department of Computer Sciences, Air University Multan Campus, Multan, Pakistan^{2}Abdus Salam School of Mathematical Sciences, GC University, Lahore, Pakistan^{3}Department of Mathematics, University of Management and Technology, Sialkot Campus, Pakistan^{4}Department of Management Sciences, Air University Multan Campus, Multan, Pakistan

Correspondence should be addressed to Qammar Rubbab

Received 1 December 2016; Accepted 15 January 2017; Published 8 February 2017

Academic Editor: John D. Clayton

Copyright © 2017 Qammar Rubbab et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

The helical flows of couple-stress fluids in a straight circular cylinder are studied in the framework of the newly developed, fully determinate linear couple-stress theory. The fluid flow is generated by the helical motion of the cylinder with time-dependent velocity. Also, the couple-stress vector is given on the cylindrical surface and the nonslip condition is considered. Using the integral transform method, analytical solutions to the axial velocity, azimuthal velocity, nonsymmetric force-stress tensor, and couple-stress vector are obtained. The obtained solutions incorporate the characteristic material length scale, which is essential to understand the fluid behavior at microscales. If characteristic length of the couple-stress fluid is zero, the results to the classical fluid are recovered. The influence of the scale parameter on the fluid velocity, axial flow rate, force-stress tensor, and couple-stress vector is analyzed by numerical calculus and graphical illustrations. It is found that the small values of the scale parameter have a significant influence on the flow parameters.

#### 1. Introduction

The existence of couple stresses in continuum mechanics is a consequence of the discrete character of the matter at the finest scale. Also, the noncentral property of forces between elementary particles of matter leads to the appearing of couple stresses [1]. E. Cosserat and F. Cosserat [2] were the first to have incorporated the couple stresses by considering the oriented material point triads and independent microrotation in theory of continuum solids. Toupin [3], Mindlin and Tiersten [4], and Koiter [5] have developed theories with couple stresses, in which the rigid body motion of the infinitesimal element of matter at each point of the continuum is described by six degrees of freedom.

Based on these theories with couple stresses, Stokes [6] has developed the theory of couple-stress fluids in order to study the size-dependent behavior of flows. Stokes’ theory represents the simplest generalization of the classical theory of fluids, which considers the presence of couple stresses. All the above models have various inconsistencies. In the context of Stokes’ theory, the indeterminacy of the spherical part of couple-stress tensor is a major shortcoming. Another problem appears in the theory of Mindlin and Tiersten; namely, in their theory, the constitutive relation for the force-stress tensor contains the body couples. The indeterminate spherical part of couple-stress tensor is ignored in Stokes’ theory, without justifications. Also, the presence of body couples in the constitutive relations was neglected by Stokes in his theory.

Hadjesfandiari et al. [1, 7, 8] have solved these issues by using arguments based upon the energy equation along with kinematical considerations. They developed a fully consistent couple-stress theory for solids and fluids, which establishes the couple-stress tensor and mean curvature rate tensor as skew-symmetric energy conjugate measures. Their results could be useful for examining fluid flow problems influenced by the mechanics at small scales. Significant application areas for size-dependent fluid mechanics involve modeling of blood flows, lubrication problems, liquid crystals, and polymeric suspensions. This branch of fluid mechanics has attracted a growing interest from researchers in the field.

Bakhti and Azrar [9] studied the steady flow of a couple-stress fluid through constricted tapered artery under influence of a transverse magnetic field, moving catheter, and slip velocity. Solutions to velocity and shear stress are expressed with Bessel’s functions. Pralhad and Schultz [10] used the couple-stress fluid model for the study of the steady flow of blood through stenosed artery. Blood velocity, the resistance to flow, and shear stress distribution have been obtained. Verma et al. [11] have studied the blood flow in a stenosed tube considering the blood as a couple-stress fluid. Effects of slip velocity in stenosed tube were highlighted. Shenoy and Pai [12] have made the static analysis of a misaligned externally adjustable fluid-film bearing including turbulence and couple-stress effects in lubricants blended with polymer additives. Naduvinamani and Patil [13] obtained a numerical solution to finite modified Reynolds equation for couple-stress squeeze film lubrication of porous journal bearings.

The unsteady three-dimensional flow of couple-stress fluid over a stretched surface with mass transfer and chemical reaction was investigated by Hayat et al. [14]. Devakar and Iyengar [15] investigated the generalized Stokes’ problems for incompressible couple-stress fluids. The effects of kinetic helicity (velocity-vorticity correlation) on turbulent momentum transport were investigated by Yokoi et al. [16, 17]. Other interesting topics can be found in references [18–21].

We must mention that the models considered in the above articles are based on the theory developed by Mindlin, Tiersten, and Stokes.

In the present paper we consider the consistent theory of couple-stress fluids elaborated by Hadjesfandiari and his coworkers [1, 8] and, we study the helical flows of couple-stress fluids within a straight circular cylinder under general boundary conditions. In the studied problem, the fluid motion is generated by the helical motion of the cylindrical surface with the time-dependent axial and azimuthal velocities and by the time-dependent couple-stress vector on the cylinder surface. The nonslip conditions are, also, considered. Using suitable nondimensional variables, we determine analytical solutions to the axial and azimuthal velocities by means of the integral transform method. The axial angular velocity and the flow rate are also obtained. Components of the nonsymmetric force-stress tensor and the couple-stress vector are determined from the constitutive relations and velocity field. Obviously, if the characteristic length of couple-stress fluids is zero, we recover results for the classical fluid. The obtained analytical solutions are used in order to perform numerical calculations using the Mathcad software for particular external loadings. The results are graphically presented. It is found that all flow parameters are influenced by the fluid scale parameter. The significant influence is obtained for small values of the scale parameter.

#### 2. Statement of the Problem

We consider a homogeneous, incompressible, viscous couple-stress fluid flowing inside a straight circular cylinder of radius . The cylindrical coordinate system has the -axis identical with the cylinder axis. The governing equations of the couple-stress fluid are as follows [1, 8]:

(i) Continuity equation:where is the fluid velocity.

(ii) Equation of linear momentum (the body forces are neglected):where is the fluid density, is the dynamic viscosity, is the viscosity coefficient of couple-stress fluid, and is the thermodynamic pressure.

(iii) The constitutive equations:

The nonsymmetric force-stress tensor:whererepresent the strain rate tensor and the angular velocity tensor, respectively.

The polar couple-stress vector:

In this paper, we consider that the velocity field and pressure are functions of the form

Equation (1) is identically satisfied and (2)–(6) become

We consider the following initial-boundary conditions:

Functions , , , and are piecewise continuous functions on ; for every , they have exponential order at infinity and .

We define the characteristic material lengthwhich is absent in classical fluid mechanics but is fundamental for couple-stress fluids.

Introducing the nondimensional variablesinto (8)–(10) and dropping the star notation, we obtain the following nondimensional problem:

In the end of this section, we give two lemmas regarding some properties to the operators from (16) and (17). These properties will be used in order to find solutions of the above problem.

Lemma 1. *If are the positive roots of Bessel function , , and , then,*(i)* *(ii)* *

*Lemma 2. If are the positive roots of Bessel function , , and , then,(i) (ii) *

*The demonstration of relations (26)–(29) is made easy using integration by parts and properties of Bessel functions [22, 23]. It is noted that is zero-order finite Hankel transform of function and is the finite Hankel transform of the first order of function , respectively.*

*3. Solution of the Problem*

*3. Solution of the Problem*

*In order to find the solution of problem (15)–(25), we use the Laplace transform with respect to the variable time and finite Hankel transform with respect the radial coordinate [24, 25].*

*3.1. Velocity Field*

*3.1. Velocity Field*

*Applying the Laplace and finite Hankel transforms to (16) and (17), using the initial and boundary conditions (22)–(25) and Lemmas 1 and 2, we obtain the following transformed equations:where , , , , , and are the Laplace transforms of functions , , , , , and , respectively. Equations (30) can be written in the suitable forms*

*Now, using the integralsand applying the inverse Hankel and Laplace transforms, we obtain closed forms to azimuthal velocity and axial velocities as*

*In the above relations we have used the notation .*

*By introducing notations(33) are written as*

*It is observed that functions and have the property .*

*The axial angular velocity (the spin vector) is given byand the flow rate is*

*3.2. Force-Stress Tensor and Couple-Stress Vector*

*3.2. Force-Stress Tensor and Couple-Stress Vector*

*Replacing (35) in (18)–(21) and performing calculations, we get the following expressions for the force-stresses and couple-stress vector:*

*It is easy to see that, for (the ordinary Newtonian fluid), the force-stress tensor becomes a symmetric tensor and the couple-stress vector is zero.*

*3.3. Particular Case (Constant Velocity and Couple Stress on the Boundary)*

*3.3. Particular Case (Constant Velocity and Couple Stress on the Boundary)*

*Let us consider the following boundary conditions:where , , , and are constants and is the Heaviside unit step function. In this case, the derivatives of functions given by (39) are being the Dirac distribution. Functions and are written in simpler forms, as*

*4. Numerical Results and Discussion*

*4. Numerical Results and Discussion*

*Unsteady helical flows in the consistent theory of couple-stress fluids were considered, under general boundary conditions. By using suitable nondimensional variables, the governing flow equations are obtained in the dimensionless form. It is important to note that these equations contain as parameter the dimensionless scale flow parameter , defined as the square of the rate between the characteristic material length and the radius of circular cylinder. As a result from the studied particular problems, the scale parameter has a significant influence on the fluid behavior. Obviously, if the scale parameter equals zero, results corresponding to the Newtonian fluid are obtained. Solutions for fluid velocity, nonsymmetric force-stress tensor, and couple-stress vector were obtained using integral transforms method (Laplace transform with respect to the time variable and finite Hankel transform with respect to the radial coordinate). The axial angular velocity and the flow rate are also determined.*

*The obtained solutions contain in their expressions the positive roots of the Bessel functions and , denoted by and , roots which were generated by means of Mathcad subroutine “.” In the numerical simulations, we have used , values for which the numerical approximation accuracy is very good.*

*In our study, the azimuthal velocity, axial velocity, and the components of the couple-stress vector are given on the cylinder surface as arbitrary functions of the time ; therefore, the obtained solutions can generate solutions to various problems with practical applications.*

*The numerical results for Figures 1–3 and for Table 1 were generated under conditions .*