On the Perturbation of the Three-Dimensional Stokes Flow of Micropolar Fluids by a Constant Uniform Magnetic Field in a Circular Cylinder

Modern engineering technology involves the micropolar magnetohydrodynamic flow of magnetic fluids. Here, we consider a colloidal suspension of non-conductive ferromagnetic material, which consists of small spherical particles that behave as rigid magnetic dipoles, in a carrier liquid of approximately zero conductivity and low-Reynolds number properties. The interaction of a 3D constant uniform magnetic field with the three-dimensional steady creeping motion Stokes flow of a viscous incompressible micropolar fluid in a circular cylinder is investigated, where the magnetization of the ferrofluid has been taken into account and the magnetic Stokes partial differential equations have been presented. Our goal is to apply the proper boundary conditions, so as to obtain the flow fields in a closed analytical form via the potential representation theory, and to study several characteristics of the flow. In view of this aim, we make use of an improved new complete and unique differential representation of magnetic Stokes flow, valid for nonaxisymmetric geometries, which provides the velocity and total pressure fields in terms of easyto-find potentials. We use these results to simulate the creeping flow of a magnetic fluid inside a circular duct and to obtain the flow fields associated with this kind of flow.


Introduction
Many chemical, biochemical, and other industrial or biological processes employ solid or soft matter in the form of small ferromagnetic particles, which are embedded in a Newtonian fluid and react in the presence of a magnetic field.Several engineering applications with high technical complexity that cover the large specific area that is offered by such systems have an inherent interest of physical and mathematical nature 1 .Those systems involve the pure hydrodynamic motion of the aforementioned aggregates of particles, where the application of a magnetic field perturbs the flow and the coexistence of liquid and magnetic properties provides us with useful information on the physical problem that has to be solved in each or biological systems that we have already described consider pipe flow problems, which can be modeled by the Stokes flow 7 approximation in the presence of a magnetic field, which involve the creeping motion of magnetic fluids in low Reynolds number 7 .Therefore, the aforementioned partial differential equations reduce to a simpler shape that is similar to the well-known Stokes equations 7 , the so-called magnetic Stokes equations, which are presented in every detail in 6 .Within the frame of the already known potential representation theory developed in early years 8-10 , a new complete and unique differential representation of magnetic Stokes flow has been constructed 6 , valid for nonaxisymmetric geometries, which provides the velocity and total pressure fields in terms of easy-to-find potentials, via an analytical fashion.This representation considers steady state, neglecting the current distribution in the fluid and approximating the magnetization by its equilibrium expression more specifically by the constant saturation expression for extremely strong magnetic fields at very low temperatures .These two last assumptions arise from the small size or the low concentration of the ferromagnetic particles neglecting magnetorestrictive effects and from the fact that, for creeping flow, very low velocities are considered.Such an analytical approach of creeping flow Stokes flow under the action of a magnetic field, arbitrarily orientated, in three-dimensional interior and exterior flow fields, appears for the first time, and, as it will be shown later, it is the key to our work.
A very large part of the applications in the technological field and in general biomagnetic processes involves the micropolar flow of magnetic fluids inside pipes or ducts and their reaction to the presence of magnetic fields.An approach of one-dimensional case appears in 1 by Berkovski and Bashtovoy, where it is examined under various strict conditions the theoretical analysis of a steady laminar flow in a cylinder under the action of an axial or transverse magnetic field, depending only on the axial direction of the flow.Similar simplified solutions like the problem of a plane Couette flow have been also developed by Rosensweig in 4 , while an interesting analytical work by Verma 11 on ferrohydrodynamics deals with a two-dimensional flow in Cartesian coordinates, introducing the well-known Stokes stream function 2 .Although the study of the reaction of micropolar fluids when they are disturbed by several types of magnetic fields is of great interest, it still remains useful to report papers concerning the flow of fluids with micropolar properties inside pipes or around bodies, where some works retain and others omit considering creeping flow the convection terms in the Navier-Stokes equations.Recently, a pure analytical solution was obtained by Sastry and Mohan Rao 12 , who discussed the micropolar fluid flow arising due to oscillations of a plane, when the system is subject to uniform rotation, whilst a hybrid paper of Calmelet-Eluhu and Majumdar 13 investigates the internal flow of a micropolar fluid inside a circular cylinder, which is subject to longitudinal and torsional oscillations, where the obtained analytical results are followed by numerical analysis.Under this aspect, Weng et al. 14 employed the theory of micropolar fluids in order to study the stability problem of flow between two concentric rotating cylinders, providing analytical and numerical results as well.On the other hand, Stokes flow in micropolar fluids appears frequently in the recent years, providing interesting information about the creeping flow of magnetic fluids.Indeed, Faltas and Saad 15 examined the axisymmetric Stokes flow of a sphere bisected by a free surface bounding a semi-infinite micropolar fluid, and three years later Sherif et al. 16 presented the Stokes axisymmetrical flow caused by a sphere translating in a micropolar fluid perpendicular to a plane wall at an arbitrary position from the wall.In addition, the same year, Moosaie and Atefi 17 provided an analytical solution for the creeping flow of a micropolar fluid past a rotating circular cylinder of infinite length in spanwise direction.Nevertheless, those analytical or semianalytical works introduce several physical or mathematical assumptions and conditions in order to overcome certain difficulties arisen by the particular analytical investigation.Those constraints multiply rapidly when someone has to face the problem of the interaction of micropolar fluids with magnetic fields and then the difficulties increase exponentially.Hence, the use of purely numerical methods was obligatory.Within this frame, two important works by Papadopoulos and Tzirtzilakis 18 as well as Khashan and Haik 19 are involved with the biomagnetic fluid dynamics under the action of a magnetic field.Moreover, in the aim of investigating the reaction of liquid metals to magnetic fields, several authors have taken into account the electric conductivity of such liquids.Under this aspect, two works by Sellers and Walker 20 and Witkowski and Walker 21 have used the principles of magnetohydrodynamics in order to construct numerical solutions for the flow of liquid metals, taking into account the electric current density.Recently, Papadopoulos et al. 22 studied numerically a first approach of the new theoretical model presented in 6 by examining the magnetohydrodynamic flow of a micropolar magnetic fluid in a straight square duct.Under this point of view, Vafeas et al. 23 provided analytical and numerical results about the flow of a micropolar fluid under the effect of a line dipole.However, all of these methods involve the use of elaborate computer codes for each case under consideration.But there is always room and need for strong analytical methods, which incorporate the appropriate geometrical and physical characteristics with the minimum of the simplifications.It is to this end that Stokes flow serves as the platform of the necessary assumptions that have to be made in order to keep consistency between the mathematics and the physics in real flow situations.During the last century, several authors have employed the creeping flow 7 in order to solve the corresponding physical boundary value problems.For example, such a large area of applications concerns the well-known particle-in-cell models for the Stokes flow through relatively homogeneous swarm of particles inorganic, organic, biological .Those systems provide a relatively simple platform for the analytical or semianalytical solution of heat or mass transport problems.Two reference papers by Dassios et al. 24 and Vafeas and Dassios 25 on this matter involve the solution of such problems in axisymmetric spheroidal and in 3D ellipsoidal coordinates, respectively.More specific, in 25 the analytical solution for the flow fields has been obtained with the aim of the Papkovich-Neuber differential representation 8, 10 , which was the motivation for the construction of a new differential representation for magnetic Stokes flow 6 .
Until nowadays, one can rarely retrieve reports of works in the literature that capture in an analytical fashion the influence of magnetic fields to flows of micropolar or electrically conductive fluids, where the complete anisotropy of the three-dimensional space must be seriously taken into account.For example, Martin Witkowski et al. 26 treated the nonaxisymmetric flow of an electrically conducting liquid in an insulating cylinder with a spatially uniform, transverse, rotating magnetic field.Here, in this paper, we obtain the analytical solution to the fully three-dimensional Stokes flow problem of a micropolar fluid with low Reynolds number properties neglect of the convection terms in the Navier-Stokes equations inside a circular cylinder, where the flow is perturbed by a constant vector uniform magnetic field that is arbitrarily orientated.In order to obtain the flow fields, we use a differential representation, which is based on the one developed in 6 , and it provides the velocity and the total pressure of the fluid in terms of easy-to-find potentials that contain the magnetic field.The representation used here considers the same assumptions made in 6 , namely, steady state flow and no current distribution in the fluid.The main improvement of the present representation compared to the one in 6 is that it approximates the magnetization by the more general equilibrium expression which contains the measure of the general applied magnetic field and not by the constant saturation expression.Thus, it has increased accuracy and more general applicability without being more complicated than the one in 6 .An outline of the differential representation that we employ as well as the additional steps needed for the present improved representation of the flow fields is presented in a later section.
Hence, under the assumption of low Reynolds number and small size or low concentration of ferromagnetic particles in the micropolar fluid, the general differential solution provides the velocity and the total pressure fields via potentials that satisfy easy-handled equations, which contain the applied magnetic field.For a constant magnetic field in the three dimensions as in our case, those equations reduce to Laplace equations for each potential function.Consequently, the unknown potentials are given via infinite series expansions of internal harmonic eigenfunctions in circular cylindrical coordinates 27-29 , since we face an interior flow problem inside a circular cylinder, and, then, the 3D flow fields are provided via closed forms of full series expansions through the differential representation.That way, the generality of the potentials is inherited to the magnetic flow fields, via a complete set of unknown constant coefficients that have to be calculated explicitly.In order to achieve that, we consider a circular cylinder of a finite length for our problem, and we supplement the general solution with the appropriate boundary conditions as a fair approximation for the creeping flow of a micropolar fluid inside a finite circular cylinder in the presence of a constant magnetic field.These conditions are nonslip condition on the wall of the cylinder, an imposed axial velocity of a known general form with no transversal components at the entrance of the cylinder, and cancellation of the axial derivatives of the transversal components of the velocity as well as zero axial stresses at the exit of the cylinder.We point out that the last condition has been introduced by Gresho 30 and Ganesh 31 .In order to secure consistency with the physical requirements, those conditions are supplemented with the demand of conservation of mass at the edge of the finite cylinder Bruneau and Fabrie 32 .All the boundary conditions reported here are extensively used in many papers, where their basic characteristics are concentrated within the references 30-32 .Applying the conditions to the general differential solutions for the flow fields, we perform many tedious and long calculations in order to evaluate the unknown constant coefficients and obtain the velocity and the total pressure in a closed analytical form of expansions of infinite series, in terms of the applied 3D constant magnetic field, of the interior circular cylindrical eigensolutions, and of the certain hydrodynamic or magnetic parameters.In addition, we accomplish to reveal the effect of the applied magnetic field on the magnetic Stokes flow fields and on the viscosity of the fluid, where its increase is clearly shown by a factor, which is the ratio of the apparent viscosity of the flow over the hydrodynamic viscosity of the fluid.The aforementioned closedtype solution is general and valid for the nonaxisymmetric circular cylindrical geometry, and all the details as well as the difficulties arisen during the calculation process are discussed extensively in the corresponding section of the paper.Here, we add that our calculations were made within the classical mathematical analysis framework using reliable bibliography 27-29, 33 , while a mathematical technique based on a theory for the completeness of the differential representations by Eubanks and Sternberg 34 also used in 25 was imposed.
The analytical section of the present paper is followed by the application of the obtained solution to the computation of the velocity and total pressure fields of creeping flow in a circular duct.The results include plots that depict the development of the flow as the magnetic fluid moves downstream of the duct, under the effect of the imposed magnetic field.We also prove that the additional friction losses depend upon the effective viscosity of the flow 35 ; hence, we present plots that depict the variation of the additional viscosity due to the magnetic field versus several characteristic parameters of the flow, such as the magnetic field magnitude and the concentration of particles inside the ferrofluid.
Finally, the necessary mathematical material from the theory of harmonic functions in circular cylindrical coordinates and the Bessel functions as well as some useful formulae associated with trigonometric and hyperbolic functions is collected in the appendix, along with the presentation of various vector identities.

Mathematical Formulation and the New Differential Representation for Magnetic Stokes Flow
In order to formulate our boundary value problem, we consider a finite circular cylinder of radius α and of length L, and, under steady state conditions with no time dependence, we examine the interior creeping flow very small Reynolds number 7 of an incompressible constant mass density ρ and viscous constant dynamic viscosity η micropolar fluid inside the smooth bounded three-dimensional cylindrical domain V R 3 that is confined by the cylinder's prescribed dimensions.In what follows, every field will be written in terms of the position vector r x 1 x 1 x 2 x 2 x 3 x 3 expressed via the Cartesian basis x i , i 1, 2, 3, in Cartesian coordinates x 1 , x 2 , x 3 , where this dependence will be omitted for convenience in writing.This coordinate system is placed in such a way so as the x 1 -axis follows the direction of the axis of symmetry of the finite circular cylinder, while the other two axes lie on the perpendicular plane.On the other hand, its connection with the appropriate r, ϕ, z circular cylinder coordinate system is provided in the appendix via relation A.12 .Micropolar fluids 1 are characterized by colloidal suspension of very small solid ferromagnetic particles, which are considered as spherical of radius r p and of density ρ p due to their minor size and are finely divided in a continuous liquid medium that is considered nonconductive without any Lorentz force density 1, 6 in our case.This last hypothesis conforms to physical reality, since the volumetric concentration of the ferromagnetic particles remains in low levels, and, then, the current distribution is negligible.The ferrofluid consists of the carrier fluid and the suspension of magnetic particles, which behave as rigid magnetic dipoles.The particles do not interact to form agglomerations due to the Brownian motion, which is responsible for the stability of the ferrofluid.Hence, the application of an external magnetic field will prevent the rotation of each particle, increasing the effective viscosity of the fluid, and will change the total pressure.Under the aforementioned assumptions, the governing micropolar hydrodynamic equations of our physical problem relate the velocity field v and the total pressure field P p ρgx 3 with the applied magnetic field H of measure H |H|. It is noted that p is the thermodynamic pressure, g −g x 3 defines the acceleration of the gravity of measure g, and ρgx 3 refers to a hydrostatic pressure force, which corresponds to a height of reference in the x 3 -direction.Moreover, the induced magnetic field is taken approximately equal to zero, which is true in many applications.The governing equations are the Stokes magnetohydrodynamic equation for the creeping motion of magnetic fluids 1, 4, 6 and the continuity equation which secures the incompressibility of the ferrofluid.They are both expressed in terms of the well-known differential operators ∇ and Δ see also A.15 .Here, τ S r 2 p ρ p /15η 0 is the relaxation time of particle rotation η 0 corresponding to the rotational viscosity , I 8πr 5  p ρ p n/15 is the sum of moments of inertia of the spherical particles per unit volume n being the number of particles per unit volume , τ B 4πηr 3 p /KT is the relaxation time of Brownian rotation K being the Boltzmann's constant and T denoting the temperature, which we consider as constant for the present isothermal problem , and μ 0 is the magnetic permeability of the free space.Furthermore, as it is clearly explained in 6 , for creeping flow, as in our case, small velocities are obtained, and then the magnetization M of measure M |M| becomes approximately colinear with H. Their connection is provided via the equilibrium magnetization M 0 and in terms of the Langevin function L ξ by the equation where m specifies the magnetic moment of a particle, while M 0 M 0 H .It is obvious that, in some cases, we have to take into consideration the magnetization of the particle itself M p , which is connected with the saturation magnetization M s by the relation where φ is the volumetric concentration of particles.On the other hand, as ξ → 0, then L ξ → 0, whilst as ξ → ∞, it is easily verified that L ξ → 1.Hence, M 0 can sometimes be approached by the constant saturation expression for extremely strong magnetic fields at very low temperatures.The vorticity of the fluid Ω of measure Ω |Ω| is expressed as

2.5
The quantity Ωτ B stands for the crucial characteristic that controls the nature of the micropolar flow.For instance, Stokes micropolar flow is confirmed by the condition Ωτ B 1 that reflects small velocities.For convenience to our calculations, we define the dimensionless parameter χ H as which is the ratio of the total viscosity of the micropolar fluid in the presence of magnetic particles that can respond to a magnetic field over the viscosity of the fluid in the absence of magnetic particles.The term Mathematical Problems in Engineering appearing in 2.1 is the additional effective viscosity caused by the interaction between the magnetic particles and the applied magnetic field.Inserting definitions 2.6 into the momentum equation for creeping magnetic flow 2.1 , we obtain the equivalent form ηχ H Δv ∇P − h H ∇ ln H .

2.8
The function h H , which is defined in 2.6 , stands for the effect of the Brownian motion of the particles on the total pressure.If there are no magnetic particles n φ 0 or if there is no magnetic field H 0 , then from 2.3 we have M 0 0, which means that h H 0 and, consequently, χ H 1 as it is revealed from relation 2.6 .Then, the pair of relationships 2.8 and 2.2 reduces to the already known Stokes equations for the hydrodynamic creeping flow 7 .
Our first purpose is focused on the analytical treatment of the coupled equations of momentum 2.8 and continuity 2.2 .Under this aspect we search for analytical solutions of these equations in the forms of differential representations 8-10 , which provide us with the velocity field v and the total pressure field P in terms of differential operators that act on particular potentials and in terms of the applied magnetic field H. Within this frame, in publication 6 , an extended paragraph with a proved theorem was devoted to that matter in the case where the equilibrium magnetization is taken nearly equal to the constant saturation magnetization M 0 ∼ M s nm .Nevertheless, in this paper it is our goal to extend this theorem in the case where the equilibrium magnetization is provided by its most general expression 2.3 , which contains the measure of the general applied magnetic field.In view of that, we present the following improved differential representation for magnetic Stokes flow, which provides the velocity and the total pressure of the fluid in terms of easy-to-find potentials that contain the magnetic field, in the more general case, where M 0 is given by 2.3 ; that is, where P 0 is a constant reference pressure and Φ, Ψ are potentials that satisfy respectively.The representation 2.9 -2.11 is the new, more complete, general solution that we present for the first time in this paper, which provides us with the magnetic flow fields of 2.8 and 2.2 in an analytical fashion.If there are no magnetic particles n φ 0 or if there is no magnetic field H 0 , we obtain M 0 0, h H h H 0, χ H 1, and our differential representation 2.9 -2.11 reduces to the already known general solution for Stokes flow 10 .The proof of completeness and of uniqueness of the differential solution 2.9 -2.11 follows the same steps provided in 6 by making the change h 6 H → h H − h H . Thus, in order to avoid repetition we choose to omit the proof and refer to 6 for further analysis and elaboration.However, it is not difficult to show that solution 2.9 and 2.10 with 2.11 satisfy 2.8 and 2.2 .Indeed, with the aid of identity A.1 we transform the momentum equation 2.8 to

2.12
Utilizing identity A.11 and substituting the velocity 2.9 and the total pressure 2.10 into 2.2 and 2.12 , then the last two are immediately satisfied, taking into account the relationships for the potentials 2.11 .
Our final task is to find an analytical form for h H , whose gradient appears in 2.10 .In order to achieve this, we proceed as follows.We make use of 2.3 , 2.6 , and 2.10 to calculate without loss of generality, since the arbitrary constant of the integration can be embodied inside the arbitrary constant pressure of reference P 0 from 2.10 .Expression 2.14 provides us with an analytical expression for h H .It is also easy to obtain from definition 2.6 with the aid of 2.3 a corresponding relation for h H , which is

2.15
Combining relationships 2.14 , 2.15 and after simple analytical manipulations, we take where ξ contains the measure of the applied magnetic field via 2.3 .Thus, the total pressure 2.10 is calculated, and, once the potentials Φ and Ψ are obtained from 2.11 , the flow fields 2.9 and 2.10 are known.Of course, the continuity equation 2.2 is immediately satisfied.
The basic difference of the present representation compared to that of 6 is that the equilibrium magnetization is provided by its most general expression 2.3 .This difference is expressed by the additional term h H which appears in the pressure 2.10 and by the fact that both h H and h H are calculated for the magnetization M 0 , which is provided in relation 2.3 .It is noted that, when M 0 ∼ M s nm from the definition of h H in 2.10 , it is obtained that h H const. ≡ 0 without loss of generality and that h H h 6 H μ 0 M s H μ 0 nmH nKT ξ.Consequently, representation 2.9 -2.11 reduces to the already published general solution 6 .It is also useful to note that if the equilibrium magnetization M 0 is approximated by the constant value of the saturation magnetization M s , then the last term of 2.12 becomes h H ∇ ln M 0 μ 0 M 0 H∇ ln M 0 0, and consequently 2.12 coincides with the momentum equation for magnetic Stokes flow presented in 6 .Finally, in order to verify the consistency of result 2.16 with the corresponding one from 6 , we make use of 2.10 to rewrite the momentum equation 2.12 as follows: 2.17 Then, we follow an easy limiting process, where with the aid of the novel result 2.16 and relation h 6 H nKT ξ for high values of ξ , we obtain which leads us to the corresponding momentum equation for magnetic Stokes flow of 6 for very high values of ξ, that is, In this paper we investigate the special case, where the main three-dimensional flow is perturbed by a known constant vector uniform magnetic field, arbitrarily orientated in the 3D space, of the form where this situation provides us with the constant equilibrium magnetization and the constant dimensionless parameter Therefore, applying 2.21 and 2.22 into the magnetic Stokes equations 2.12 and 2.2 , we recover the corresponding differential equations The particular case of the application of a constant three-dimensional magnetic field leads via relation 2.22 to ∇ ln χ c 0 and ∇h H c 0 ⇒ h H c ≡ h c .Therefore, our improved differential representation 2.9 -2.11 of the solutions of 2.23 and 2.24 takes the form where P c ≡ P 0 h c − h c specifies the new characteristic reference pressure that is arbitrarily chosen, since the total pressure enters the momentum equation 2.23 under the gradient action.As a consequence of 2.22 , relations 2.11 reduce to Laplace's equations for the potentials Φ and Ψ, that is, respectively.Consequently, we have to find two harmonic potentials and calculate the magnetic flow fields from the differential representation 2.25 and 2.26 in circular cylindrical geometry, whenever the boundary conditions of the corresponding physical problem are known.
Our boundary value problem must be supplemented by the appropriate boundary conditions fixed to the precisely defined boundaries of the circular cylinder of radius α and of finite length L. Those conditions are the following: where the first one defines the nonslip boundary condition on the wall of the cylinder at r α, the second one refers to the entrance of the cylinder at z 0 demanding the transversal components of the velocity to vanish and specifying a known imposed axial velocity v r, ϕ , while the third one cancels the axial derivatives of the transversal components of the velocity and demands that the axial stresses are set to zero see also 30, 31 at the exit of the cylinder at z L.Although the second part of boundary condition 2.30 comprises the thermodynamic pressure, it is to our convenience to consider small height variations, so as to neglect the hydrostatic pressure force and, consequently, to set the total pressure from 2.26 approximately equal to p, meaning that P ∼ p.In addition, in order to secure consistency with the physical requirements of our problem, these conditions are supplemented with the demand of conservation of mass at the edge of the finite cylinder see also 32 , that is, In the limiting case, where the length of the cylinder tends to infinity, we demand that the velocity field at the outlet of the duct obtains a fully developed profile lim respectively, where −A ≡ dP P /dz < 0 is the constant pressure gradient at the z-direction of the magnetic Poiseuille flow v P , P P .Immediate integration gives We point out that the constant of the integration in order to find P P has been taken equal to AL in order to be in accordance with 2.32 and the same time keep consistency at infinity.Moreover, the effect of the magnetic field 2.20 is inherited to the velocity profile 2.32 , which carries the parameter χ c in the absence of particles or magnetic field, χ c 1, and we obtain the hydrodynamic Poiseuille flow .On the other hand, the nature of the asymptotic condition 2.32 determines the character of the flow fields v and P .In simple words, condition 2.32 allows us to decompose the flow fields as follows:

The Magnetic 3D Flow Fields in Circular Cylindrical Coordinates
According to the analysis described above, the required 3D flow fields 2.34 of our problem, in view of 2.32 , 2.33 , and 2.34 , assume the form in terms of the functions χ c provided by 2.22 Since our case involves an interior flow problem, it is imposed the use of regular solutions on the axis of symmetry of the circular cylinder r 0 , which means that the Neumann functions N n μr must be excluded from any general harmonic expansion of the type A.30 .The constant parameter μ ∈ R from the separation of variables of the Laplace equation in our system see the appendix for that matter will be determined from boundary condition 2.28 in the end of our analytical procedure.However, until then, we will introduce the symbol " μ • • • ", which denotes integration if μ takes continuous values or summation in the case where μ is a parameter with discrete values.Consequently, the complete representation of the regular potentials Φ and Ψ, which belong to the kernel space of Δ, is where, for n ≥ 0, μ ∈ R, and i 1, 2, 3, 4, the coefficients e μ, i n x 3 and d μ, i n denote the unknown vector and the scalar constant coefficients of the harmonic potentials Φ and Ψ, respectively.On the other hand, within the frame of A.30 the Hfunctions, usually named as eigenfunctions, for every kind i 1, 2, 3, 4 are where n ≥ 0, while all the needed information for the Bessel, the trigonometric, and the hyperbolic functions is summarized in the appendix.Inserting the potentials 3.3 in the flow fields 3.1 and by extensive use of identities A.2 , A.3 , and A.5 , we derive the relation Mathematical Problems in Engineering for the velocity field, while for the total pressure field we obtain in the cylinder's prescribed dimensions V R 3 .That way, the generality of the potentials is inherited to the magnetic flow fields 3.8 and 3.9 , via the set of unknown constant coefficients e μ, i n and d μ, i n for n ≥ 0, μ ∈ R, and i 1, 2, 3, 4, which have to be calculated explicitly from the proper conditions 2.28 -2.31 .Additionally, the parameter μ will be evaluated from the same conditions.
Since the vector character of the vector harmonic potential Φ is reflected upon the corresponding constant coefficients, which are written in Cartesian coordinates, we are obliged to work in the Cartesian system.Thereupon, before we proceed to the boundary conditions, it is necessary, for our convenience in calculations, to evaluate the gradient of the harmonic H-functions 3.4 -3.7 that appear in the flow fields 3.8 and 3.9 .Obviously Δ ∇H , and it is feasible to be written as a function of them in Cartesian coordinates.Therefore, we will work on Cartesian coordinates as far as the vector character of the velocity 3.8 is concerned, and using the transformation A.14 we will return to the circular cylindrical basis r, ϕ, z see also A. 13 .Hence, in view of A.13 , A.15 and the basic relations for the trigonometric functions A. 25 for every n ≥ 0 and μ ∈ R, where the prime at the Bessel functions denotes derivation with respect to the variable μr ∈ R. In the same way we calculate the rest three gradients, that is, 3.12 for every n ≥ 0 and μ ∈ R, where it is obvious that Δ ∇H Here, we must add that, in order to avoid negative values of n within the above relationships, we impose 3.13 Expressions 3.11 and 3.12 can be collected in the general formulae for every n ≥ 0, while which by virtue of 3.14 incorporates with 3.8 to the following velocity field:

Mathematical Problems in Engineering
At this point we manipulate properly the H-functions H

and H μ, b i
n−1 by a certain readjustment of the index n at the series to rewrite the velocity 3.16 as 3.17 where we note that z x 1 , the parabolic velocity profile u r is given by 3.2 , and the new constant coefficients into the velocity field 3.17 are provided as a function of a and d   μ, i  n for n ≥ 0, μ ∈ R, and i 1, 2, 3, 4, which are grouped via the simple relations

3.21
Similarly, formulae 3.14 and the group of constant coefficients 3.18 -3.21 are interrelated, and the total pressure 3.9 becomes for n ≥ 0 and i 1, 2, 3, 4, as well as the parameter μ ∈ R.During this process, we must keep in mind that We begin with the first part of the inlet boundary condition 2.29 at the entrance of the circular cylinder z 0. Applying this condition on the velocity field 3.17 , condition r • v 0 at z 0 results in We notice that both conditions 3.24 and 3.25 contain the Heigenfunctions of kind i 1, 3 relations 3.4 , 3.6 , where those survive for z 0 as it is revealed from A.29 .Hence, we utilize definitions 3.4 and 3.6 , relations A.13 or A.14 , as well as 3.19 -3.21 , in order to handle firstly condition 3.24 with extensive use of the recurrence relations of the trigonometric functions A.25 -A.28 .Therefore, after long and tedious calculations on 3.24 with orthogonality arguments of the trigonometric functions sin nϕ, n ≥ 1, and via the definition for the Bessel functions A.19 , we obtain the following relations for the constant coefficients: 0 for n ≥ 0, μ ∈ R.

3.28
Since we have finished with condition 3.24 , we deal, now, with condition 3.25 , and we observe that results 3.27 , 3.28 correlate with 3.18 -3.21 so as to satisfy automatically this condition 3.25 .Hence, we completed our analysis with the first part of the inlet boundary condition 2.29 , which implied see 3.27 and 3.28 and via 3.18 -3.21

3.30
It is worth mentioning that the very same results are obtained in an easier way by using the velocity field from 3.8 and not from 3.17 .Nevertheless, it was a worthwhile way to crosscheck the relations 3.29 .At this point, before we proceed to the next boundary conditions, we will write down the flow fields with their form until now.Under this aim, we insert relationships 3.29 , 3.30 into the flow fields 3.17 where we have used the fact that ∂ r/∂z ∂ ϕ/∂z 0 see A.13 , while for reasons of convenience to our forthcoming manipulations we incorporate the two relations of 3.33 into 0 for z L.

3.34
For the respective calculations based on 3.34 , the upper index will refer to x 2 , whilst the lower one will correspond to x 3 .Within this consideration, we apply the boundary conditions 3.34 to the velocity field 3.31 , and since z • x 2 z • x 3 0 with z x 1 , this act gives rise to Mathematical Problems in Engineering which are two similar relations at z L, and the case n 0 is taken by the general case for n ≥ 1 into 3.35 .We take into account the recurrence relations of the trigonometric functions A.25 -A.28 , we use the Bessel recurrence relation A.23 , and, after some calculations via A.19 , orthogonality of sin nϕ, n ≥ 0 and a certain readjustment of the index n at the series infer the following formulae: for n ≥ 0, μ ∈ R.

3.36
Orthogonality arguments on the functions cos nϕ, n ≥ 0 with similar readjustment of n imply

3.37
where 3.36 and 3.37 are four double-equality relationships.Recall that the upper and the lower indices reflect the vanishing of x 2 and x 3 transversal components of the velocity field at z L, respectively.
In order to proceed we have to simplify the velocity field 3.31 and the total pressure field 3.32 in terms of the expressions 3.36 and 3.37 for the constant coefficients.In view of that, we are primarily involved with the velocity 3.31 , and we smartly manipulate certain group of terms of the x 2 and of the x 3 components with respect to A.27 , A.28 , as well as with a particular change of n, to obtain × J n μr sin nϕ sinh μz ,

3.38
where we have used the two equalities of relation 3.36 and the recurrence formula A.23 .Similarly, another group of terms of the velocity 3.31 shows that where, here we omitted common steps like those followed in 3.38 , we made use of A.25 , A.26 and of relationship 3.37 .Thus, we accomplished to simplify the x 2 -component and the x 3 -component of the velocity field 3.31 .In addition, some further analysis on the second equalities of 3.36 and 3.37 gives

3.44
On the other hand, the total pressure 3.32 can also be simplified by using again definitions 3.19 -3.20 and considering 3.42 , 3.43 to evaluate

3.45
Recapitulating the effect of the boundary condition 3.34 , via 3.36 and 3.37 , on the flow fields 3.31 and 3.32 , we work as follows.In terms of 3.38 , 3.39 , 3.44 , and 3.45 , the only b's, c's, and d's constant coefficients appear inside the x 1 -component of the velocity field 3.31 , associated with each other as coefficients of the eigenfunctions J n μr sin nϕ cosh μz and J n μr cos nϕ cosh μz for every n ≥ 0 and μ ∈ R. Since they appear nowhere else in the velocity or the total pressure, we can group them as  , and e μ, 3 n for every n ≥ 0, μ ∈ R from the remaining boundary conditions.
Continuing, we handle the second part of the outlet boundary condition 2.30 , and, as we mentioned earlier, we impose no height variations by setting P ∼ p, in order to take from which, with the contribution of the z-component of the velocity field 3.48 and the total pressure 3.49 note that x 1 z, while expressions 3.2 denote p L 0 and ∂u r /∂z 0 , we obtain for z L the following condition: for every n ≥ 0, μ ∈ R, which must be determined from the continuity restraint 2.31 and the second part of the inlet boundary condition 2.29 .The constant parameter μ ∈ R will be evaluated from the nonslip boundary condition 2.28 in the end of this analysis.Afterwards we will verify that the limit 2.32 is valid.Although we deal with a well-posed boundary value problem, we observe that we have four conditions in order to calculate four constant coefficients and also parameter μ.Hence, one constant coefficient will be undetermined.However, in view of a detailed analysis on the completeness of general solutions explained in 34 and also applied in 25 , the differential representation used here 2.25 and 2.26 with 2.27 , as well as the corresponding differential solutions 8-10 , offers certain degrees of freedom to our solutions and permit the interrelation of the harmonic potentials Φ and Ψ between each other in such a way so as to avoid any kind of indeterminacies to particular analytical solutions.This interrelation does not affect the generality of the analytical solutions see, e.g., 25 , and it is applied to the constant coefficients.Therefore, we adopt the theory of 34 , and we make a mathematical technique by choosing to express the a's constant coefficients of kind i 2, 4 in terms of those of kind i 1, 3, via the relationships of proportionality n,± , for n ≥ 0, μ ∈ R, 3.54 whereas λ μ ∈ R is the new parameter to be calculated for μ ∈ R, while the second equality of 3.54 is a consequence of definitions 3.18 .Relationship 3.54 does not constrain the generality of our analytical method and at the same time will provide us with a closed unique solution.By substitution of relations 3.53 and 3.54 into the flow fields 3.48 , 3.49 , we can easily show that

3.55
for the velocity field and

3.56
for the total pressure field, respectively.Recall that the A's constant coefficients are provided via the a's constant coefficients from 3.18 and the factor χ c is given through 2.22 , while u r and p z assume the relations 3.2 .
In order to impose the continuity restriction 2.31 at z L, it would provide us with analytical convenience if we could write the velocity field 3.55 in pure cylindrical coordinates, meaning with the cylindrical unit vectors r, ϕ, z.In view of that, it is obvious that x 1 z, while substituting the unit vectors in Cartesian coordinates x 2 and x 3 from A.14 , using 3.18 , A.22 , A.23 , A.25 -A.28 and for one more time rearranging the index n inside the series, we perform some analytical steps to express the Cartesian part of 3.55 in cylindrical coordinates as cos nϕ for every n ≥ 0 and μ ∈ R on 3.59 provides us with the constant parameter λ μ as λ μ μL − 2 −1 2χ c 1 sinh 2μL 2χ c 1 cosh 2 μL − 1 with sinh 2μL 2 sinh μL cosh μL for μ ∈ R.

3.61
Actually, the flow fields retain only the parameter μ, and their equivalent form with respect to 3.61 is 3.62 from velocity field 3.55 and sin nϕ a μ, 3 n cos nϕ , 3.63 from total pressure field 3.56 , respectively, where the prime appearing at g μ μz denotes derivation with respect to the variable μz ∈ R. The new functions of expansions 3.62 and 3.63 are defined as

3.65
Here, we must recall that the A's constant coefficients are given through the a's constant coefficients via 3.18 and the factor χ c is provided by 2.22 , while u r and p z assume the well-known expressions 3.2 .
One can easily manipulate relationships 3.64 and 3.65 to prove that, for z L, for n ≥ 0, as well as the unknown parameter μ ∈ R, are bounded, the velocity 3.62 and total pressure 3.63 asymptotically lead to lim

3.67
Hence, we obtain the predictable parabolic profile of the Poiseuille flow 3.2 , and, consequently, the asymptotic limit 2.32 holds true.Our next step involves the determination of the constant parameter μ, which will be calculated by the nonslip condition 2.28 , where the velocity 3.62 vanishes v 0 at r α, that is, sin nϕ x 3 0,

3.68
Mathematical Problems in Engineering which is a three-component condition, and each component must be set to nil at r α.In the aim of satisfying all the three conditions from 3.68 at the same time, we easily observe that it is required the nihilism of the Bessel functions on the surface of the cylinder, that is, a condition that provides us with the value of the parameter μ as discrete eigenvalues via where ρ m n is the m-root m ≥ 1 of order n ≥ 0 of the Bessel functions J n ρ m n 0 .Consequently, formula 3.70 determines the parameter μ ∈ R, which is rather an expected result for such kind of problems.Since the constant parameter μ takes discrete values, it is obvious that the defined symbol " μ • • • " used inside the fields 3.62 , 3.63 must be substituted by the standard series symbol " ∞ m 1 • • • ".Thus, the orthogonality relation for the Bessel functions A.32 can be applied in our forthcoming calculations at which we will retain the μ m n -symbolism given now by 3.70 , where μ m n is the m-parameter m ≥ 1 of order n ≥ 0. Introducing the new constant coefficients that have to be calculated as and instead of 3.18 their equivalent forms for every n ≥ 0, m ≥ 1, and μ ≡ μ m n ∈ R, the velocity field 3.62 is given via the eigenvaluetype expansion

3.73
while the total pressure field 3.63 , under the aforementioned readjustment, gives rise to

3.74
The functions 3.64 and 3.65 , which are involved into the flow fields 3.73 and 3.74 , are now provided by the following formulae:

3.76
Here, μ m n , χ c and u r , p z satisfy relations 3.70 , 2.22 , and 3.2 , respectively, while the remaining constant coefficients a m n and b m n for n ≥ 0 and m ≥ 1 will be determined from the second part of boundary condition 2.29 , which is the last one to deal with.
Straightforward application of the imposed entrance velocity z • v v r, ϕ for z 0 at the velocity 3.73 note that x 1 z and use of 3. In the sense of making this work more complete, we follow some easy analytical steps, identical to those followed for the production of relationship 3.57 , with extensive use of formulae A.25 -A.28 and with a proper rearranging of the index n ≥ 0, in order to write the velocity 3.73 in pure circular cylindrical coordinates with unit vectors r, ϕ, and z.In view of this aim, we substitute the Cartesian unit vectors x i , i 1, 2, 3, from A.14 into 3.73 , and we utilize expressions 3.72 in terms of the functions

3.82
Recapitulating, we completed our analytical method with the application of all the required boundary conditions of the particular physical problem that we had to solve.Hence, given a known velocity field v r, ϕ at the entrance of the circular pipe, the velocity, and the total pressure of a magnetic fluid perturbed by the three-dimensional constant magnetic field 2.20 within the prescribed dimensions n is the m-root m ≥ 1 of order n ≥ 0 of the Bessel functions, that is, J n ρ m n 0, then obviously j m n,± ρ m n , μ m n z 0. Finally, any information needed beyond this section, which concerns the Bessel functions J n μ m n r or the trigonometric and hyperbolic functions, is collected into the appendix.

Numerical Results and Discussion
The problem under consideration is the creeping flow inside a straight circular tube under the influence of a constant vector uniform magnetic field arbitrarily orientated.The radius of the pipe is α, and its length is L. The mean axial velocity of the flow at the inlet of the pipe is given by the relation where Λ is the cross section at the inlet.Consequently, the mass flow rate at the inlet of the duct is ṁ ρUΛ.

4.2
If we consider the limiting case of a semi-infinite duct with L → ∞ and apply mass conservation at the inlet and the outlet of the duct, then it is obtained that the mass rate at which the fluid enters the cylinder must be equal to the mass rate at the exit.Taking into account the analysis presented in 2.32 , we obtain a relation that expresses the axial pressure gradient −A in terms of the mean axial velocity at the inlet U, through the fully developed velocity profile u r of 3.2 as ṁ ρUΛ ρ Λ u r rdr dϕ Λ rdr dϕ

4.4
In order for our results to have general applicability, we introduce the nondimensional variables

4.5
Thus, we can apply solution 3.81 and 3.82 to the problem of developing flow of a micropolar fluid in the entrance region of a circular duct to obtain the corresponding flow field.In all the cases that we examined in the summations of 3.81 and 3.82 , we used the first 18 terms.
In the first case the inlet velocity profile was chosen to be v r, ϕ 1 − r 5 , 4.6 which is very close to a uniform inlet velocity.u r ≡ r • v, and the dimensionless total pressure P .The plots are taken in a longitudinal-radial section of the duct, as shown in Figure 1.
It is observed that, as the flow moves downstream, it becomes fully developed.The radial velocity vanishes, the pressure drop becomes constant and uniform over the crosssection of the duct, and the axial velocity obtains a parabolic form.This last observation is clearer in the plots of Figure 2 that depict the development of the axial velocity in the longitudinal direction.
This figure shows the dimensionless axial velocity distribution u z along the radial direction for different axial lengths.It is seen that the axial velocity has initially the inlet profile 4.6 , which gradually transforms to the parabolic profile.It must be noted, however, that the inlet profile is slightly distorted in the area close to the center at r 0, where we see a small reduction of the axial velocity although the inlet profile 4.6 is nondecreasing.This is attributed to the function of the inlet profile, which changes very rapidly with r.Since the solution that we present is able to predict 3D flow without any symmetries, we chose to examine a flow with a nonsymmetrical profile.Thus, we imposed an inlet velocity profile given by the function v r, ϕ 1 − r 5 sin πr sin ϕ.

4.7
Figure 3 shows contours of the dimensionless axial velocity component u z , the dimensionless radial velocity component u r , the dimensionless angular velocity component u ϕ ≡ ϕ • v, and the dimensionless total pressure P , at several downstream locations.A general remark is that, as the flow enters the pipe, there is intense movement of the fluid in the cross-sectional plane, which gradually leads to a fully developed velocity.The pressure distribution, which is initially highly nonuniform, changes as the flow moves downstream, and it eventually reaches a uniform profile.
Another interesting characteristic of micropolar flow is the increase of the viscosity, and consequently of the friction losses, due to the magnetic field.In order to investigate these effects, we define the Reynolds number Re and the Fanning friction factor f as in 35 : where D h 2α is the hydraulic diameter and τ w is wall shear stress.For fully developed flow in a circular duct, the wall shear stress is related to the axial pressure drop via the relation 2τ w −αdP P /dz αA.Applying this relation to 4.8 and using 4.4 , we obtain the relation for the friction factor f product with the Reynolds number f Re 16χ c , 4.9 which is valid for fully developed magnetic flow in a circular duct.If there are no magnetic particles or if there is no magnetic field, then χ c 1 and 4.9 gives the well-known relation f Re 16.Thus, in order to assess the additional friction losses due to the magnetic field, we examine the parameter χ c .
Figure 4 depicts the dependence of χ c on the magnetic field H and on the volumetric concentration of magnetic particles φ, which is defined in relation 2. 4 .
The plots correspond to constant temperature T 300 K, magnetic dipole strength m 2 × 10 −18 A•m 2 , viscosity η 10 −3 kg/ms, and small magnetic particles of radius r p < 10 −7 m so that the assumptions that were made for the derivation of the Stokes equations are valid.It is noted that the parameter χ c does not depend on the particle radius r p .This can be seen from the relation 2.22 after the substitution of τ S , τ B , I and of M 0,c from 2.21 .From the plot of Figure 4 it is observed that, independently of the concentration φ, the parameter χ c undergoes a steep increase in the interval 10 3 < H < 10 5 A/m , while it remains steady for smaller or higher values of the magnetic field.Moreover, it is observed that the parameter χ c increases in a linear manner as the concentration φ increases.

Conclusions
Many applications of practical interest in modern engineering technology deal with the micropolar magnetohydrodynamic flow of magnetic fluids in the presence of magnetic fields.In this paper we examined how a 3D constant uniform magnetic field perturbs a threedimensional steady creeping motion Stokes flow of a viscous incompressible micropolar fluid of approximately zero conductivity in a circular cylinder of a finite length.The ferrofluid was considered to be a colloidal suspension of nonconductive ferromagnetic material of very small spherical particles, which follow the Brownian motion and behave as rigid magnetic dipoles, while the magnetization of the carrier liquid had been taken into account from its general equilibrium expression.We used a special case of a general three-dimensional theoretical model that governs the micropolar hydrodynamic flow in such liquids, which was based on the reduction of the partial differential equations to a simpler shape that is similar to low-Reynolds number flow and are called magnetic Stokes equations.The interaction of the applied magnetic field with the ferromagnetic particles has been inherited within those equations through the additional effective viscosity and the magnetic pressure of the fluid in terms of both its hydrodynamic and magnetic properties.
We used an improved new complete and unique differential representation of magnetic Stokes flow, drawn from the potential representation theory, which is valid for nonaxisymmetric micropolar flows and provides in an analytical fashion the velocity and total pressure fields in terms of easy-to-find potentials.Therein, we applied the proper boundary conditions according to the physical requirements on the boundaries, and we managed in a cumbersome and tedious way to obtain the three-dimensional magnetic Stokes flow fields in a closed analytical form of expansions of infinite series, in terms of the applied 3D constant magnetic field, of the interior circular cylindrical eigensolutions, and of the certain hydrodynamic or magnetic parameters.
In order to study the characteristics of the 3D flow and compute the velocity and total pressure fields associated with creeping magnetic flow, we proceeded to a necessary numerical implementation of our results in a circular duct.The development of the flow, as the magnetic fluid moves downstream of the duct under the effect of the imposed magnetic field, has been presented graphically for various imposed velocities at the inlet, while the variation of the additional viscosity due to the magnetic field versus the magnetic field magnitude and the concentration of particles inside the ferrofluid has been depicted through specific plots.
or their inverse x 1 z, x 2 cos ϕ r − sin ϕ ϕ, x 3 sin ϕ r cos ϕ ϕ. A.14 The gradient ∇ and the Laplacian Δ are provided by the expressions ∂ ∂r and then the forthcoming relationships are easily verified using the vector relations A.12 , A.13 , and A.15 , that is, as well as ∇ ⊗ r I 3 i 1 x i ⊗ x i r ⊗ r ϕ ⊗ ϕ z ⊗ z, A.17 whereas I stands for the unit dyadic and 0 refers to the zero dyadic 28 .Our boundary value problem is adjusted to the type of circular cylindrical geometry introduced here, where the x 1 -axis is the axis of symmetry of an infinite circular cylinder and the other two axes are located properly so as to obtain a r, ϕ, z right-handed system.

Bessel and Neumann Special Functions
In the sequel, we define the special functions used in this project.The Bessel functions J n x and the Neumann functions N n x 29 of order n ≥ 0 are linear independent solutions of the Bessel differential equation x 2 y x xy x x 2 − n 2 y x 0, A.18 for n 0, 1, 2, . .., which is valid for every x ∈ R.These functions are defined as follows: via the well-known Gamma function Γ κ ν 1 29 for every κ ≥ 0 and ν ∈ R. In our case, where ν n 0, 1, 2, . .., it is easily confirmed that J −n x −1 n J n x and N −n x −1 n N n x for x ∈ R. The limiting cases as x → 0, where J n x converge, or as x → ∞, where N n x are regular, are discussed in many standard references such as 29 , and it is not worth to mention them here.The Bessel and the Neumann functions satisfy the important recurrence relations

Trigonometric and Hyperbolic Functions
As far as the trigonometric functions sin nϕ and cos nϕ are concerned for n ≥ 0, which are used in this paper, the following expressions hold true: where ϕ ∈ 0, 2π stands for the azimuthal angle, taken for the first period of the trigonometric circle, while the orthogonality here is obvious for the functions sin nϕ and cos nϕ for n ≥ 0. On the other hand, the hyperbolic functions introduced in our paper are well known and can be taken also by their exponential form, that is, where z ∈ −∞, ∞ , whilst the parameter μ, as we mentioned previously, comes from the boundary conditions of a specific boundary value problem.

Harmonic Functions in Circular Cylindrical Geometry
Then, in view of our physical problem, every harmonic function w Δw r, ϕ, z 0 in the circular cylindrical coordinate system is written as w r, ϕ, z for n ≥ 0 are constants, whilst the symbol " μ • • • " refers to the constant parameter μ and denotes integration if μ takes continuous values or summation in the case where μ is a parameter with discrete values.At this point we must make the following basic remark.Here, as we pointed earlier, the parameter μ comes from the method of separation of variables for the Laplace equation in circular cylindrical coordinates, where only the particular physical problem and the boundary conditions can dictate its value.For example, as in our case, when one has to face an interior problem inside a circular cylinder of radius α ∈ R, where the boundary conditions require nihilism of the Bessel functions on the surface of the cylinder, then the regularity of the solution demands Z n μr ≡ J n μr and the boundary conditions provide us with the discrete values of the parameter as μ ≡ μ m n ρ m n /α, where ρ m n ∈ R is the m-root m ≥ 1 of order n ≥ 0 of the Bessel functions and the eigenvalue now μ m n is the m-parameter m ≥ 1 of order n ≥ 0. Hence, for this particular problem of determining parameter μ, expansion A.

3 . 26 while 27 Mathematical
for n 0 we can admit without loss of generality sin 0ϕ 0 that d μ ∈ R. On the other hand, an easy manipulation of 3.26 reveals the nihilism of the constant coefficients d n ≥ 0, μ ∈ R,3.Problems in Engineeringwhilst orthogonality of cos nϕ, n ≥ 0, on 3.24 is provided with similar relations to 3 and 3.23 , taking into account the constant parameters a i , b i , c i for i 1, 2, 3, 4 with the definitions for the H-functions 3.4 -3.7 , and we recover the velocity field .41 respectively, where in view of definitions 3.19 -3.20 and after some trivial processing, expressions 3.40 and 3.41 take the form b μ, 2 .43 respectively.The last equations 3.42 , 3.43 provide us with the simplification of the third component of the velocity 3.31 .Hence, proceeding to the x 1 -component of 3.31 , for once more we utilize the relations A.25 -A.28 , A.23 , and we apply a certain readjustment of the index n so as to write the orthogonality relation A.32 and the well-known orthogonality of the trigonometric functions, relation 3.77 reveals that n ≥ 0 and m ≥ 1, where ε n 2 for n ≥ 1, ε 0 1 and u r is given in 3.2 .It is trivial to check the nice behavior of the constant coefficients 3.78 as L → ∞.Relations 3.78 coincide with 3.72 and provide us with the final form of the flow fields 3.73 and 3.74 .Two basic special cases are the following.If the imposed velocity is only a function of the variable r, that is, v r , then from relations 3.78 the only constant coefficients that survive are the b m 0 for m ≥ 1 and the magnetic flow becomes axisymmetric.Moreover, in this case, if v r u r , then 3.78 shows that a n ≥ 0 and m ≥ 1; thus, the magnetic flow 3.73 becomes the magnetic Poiseuille flow v u r z.
sin nϕ b m n cos nϕ for n ≥ 0, m ≥ 1ρ m n being the m-root m ≥ 1 of order n ≥ 0 of the Bessel functions .Thus, we obtain the equivalent form of the velocity field 3.73 v r, ϕ, z u r

3 .83 assume the expansions 3 .
81 and 3.82 , respectively, where a m n , b m n for n ≥ 0 and m ≥ 1 have been calculated through relations 3.78 .The functions h m n μ m n z , g m n μ m n z , j m n,± μ m n r, μ m n z , and f m n nϕ for n ≥ 0 and m ≥ 1 are given by 3.75 , 3.76 , 3.79 , and 3.80 , respectively, while μ m n ρ m n /α, χ c and u r , p z satisfy relationships 3.70 , 2.22 , and 3.2 , respectively.Since ρ m

Figure 1 Figure 1 :
Figure 1: Contour plots of the dimensionless axial velocity component u z , the dimensionless radial velocity component u r , and the dimensionless total pressure P over a longitudinal-radial section of the duct.

Figure 2 :
Figure 2: Diagrams of the dimensionless axial velocity distribution u z along the radial direction for different axial lengths z.

Figure 3 :
Figure 3: Contour plots of the dimensionless axial velocity component u z , the dimensionless radial velocity component u r , the dimensionless angular velocity component u ϕ , and the dimensionless total pressure P at several downstream locations.

Figure 4 :
Figure 4: Variation of parameter χ c with the magnetic field H and the volumetric concentration of magnetic particles φ, where the plot corresponds to fully developed flow of a micropolar fluid.

with Γ κ ν 1 ∞ 0 e
n ≥ 0 and ν ∈ R, where the Bessel functions J ν x for real ν are furnished by the formula −t t κ ν dt for x ∈ R, A.21 1 x Z n 1 x for x ∈ R, A.23 defined as Z n x ≡ J n x or Z n x ≡ N n x for n ≥ 0 and x ∈ R. The Bessel differential equation A.18 and its solutions A.19 -A.21 appear in several physical boundary value problems such as ours, where under the same notation the variable x is written as x μr ∈ R, with μ ∈ R being a constant parameter, which comes from the method of separation of variables of a partial differential equation in the circular cylindrical coordinate system.This parameter is specified by the boundary conditions of the particular physical problem that has to be satisfied in each case.Upon that procedure, if the parameter μ becomes an eigenvalue μ m n ∈ R for n ≥ 0 and m ≥ 1, then x is substituted by x μr ≡ μ m n r ∈ R, whereas μ ≡ μ m n is the m-parameter m ≥ 1 of order n ≥ 0. Thus, both the Bessel and the Neumann functions satisfy the orthogonality relation n ≥ 0 and m, m ≥ 1, with δ mm being the Kronecker delta and a, b ∈ R being constants.

2 J n 1 ρ m n 2
f m n for n ≥ 0, m ≥ 1 are the new constants, while we have a 0, b α and relation A.24 takes the form , n ≥ 0, m, m ≥ 1, A.32 40 Mathematical Problems in Engineering which is the orthogonality expression frequently utilized within our calculations.In addition, this situation of an internal problem requires the use of regular solutions on the axis of symmetry of the circular cylinder r 0 , which means that only the Bessel functions J n μr ≡ J n μ m n r survive inside the expansion A.30 or A.31 , and, thus, this is reflected to the vanishing of the corresponding constants, that is, b μ n 0 for n ≥ 0 or b m n 0 for n ≥ 0 and m ≥ 1, respectively.
Straightforward calculations within the frame of identities A.3 and A.10 confirm the satisfaction of the magnetic Stokes equations 2.23 and 2.24 by the Poiseuille flow fields v P , P P given in 2.32 and 2.33 .Thus, decomposition 2.34 holds true, and the boundary condition 2.32 is substituted by lim .32 or 2.36 , 2.33 and the differential representation 2.35 in circular cylindrical coordinates, in order to construct the three-dimensional flow fields v, P in a closed 3D analytical form, accompanied by a particular numerical implementation of the results.
g , P P P P g , 2.34 where the pair of flow fields v P , P P refers to the Poiseuille flow, while the general velocity field v g and the general total pressure field P g must satisfy the differential representation 2.25 , 2.26 with 2.27 of the solutions of 2.23 and 2.24 , v g Φ − 1 2 ∇ r • Φ Ψ , P g −ηχ c ∇ • Φ P c with ΔΦ 0, ΔΨ 0, 2.35 in view of 2.20 -2.22 for χ c .
3,4, which means that ∇H -A.28 , utilizing the recurrence relations for the Bessel functions A.22 , A.23 , and in terms of the H-functions 3.4 -3.7 , for i 1, we are led to n μr cos nϕ cosh μz zμJ n μr sin nϕ sinh μz μJ n μr sin nϕ sinh μz x 1 μ cosh μz × x 2 cos ϕ sin nϕJ n μr − sin ϕ cos nϕ 22or since x i • x j δ ij δ ij being the Kronecker delta and with a proper readjustment of the index pressure field, where the corresponding pressure p z of the parabolic velocity profile is provided by 3.2 .The velocity 3.17 and the total pressure 3.23 of the micropolar fluid are ready to accept the boundary conditions 2.28 -2.30 and condition 2.31 of our physical problem, in order to determine the unknown constant coefficients a refer to the first part of the outlet boundary condition 2.30 , which forces us to set the axial derivatives of the transversal components of the velocity to nil.Since the velocity field 3.31 is written in Cartesian coordinates with z x 1 , we replace this part of 2.30 by the proper identical relationships .47 without loss of generality.We insert the simplified relations 3.38 , 3.39 , and 3.44 into the velocity field 3.31 , where with the aim of the new constant coefficients e n μr 2μ 1 μL tanh μL sinh μz − μz cosh μz constant coefficients are given through the a's constant coefficients via 3.18 , the factor χ c satisfies 2.22 , while u r and p z are given by relations 3.2 .Hence, in order to complete our analytical method and obtain the final closed form of the flow fields 3.48 and 3.49 , we have to calculate the rest of the constant coefficients, that is, the a tanh μL − λ μ μ z − L sinh μz − μz cosh μz , 3.58 for λ μ ∈ R and μ ∈ R. By virtue of condition 2.31 in circular cylindrical coordinates see the gradient from A.15 on the velocity 3.55 , with the aim of 3.57 and definition 3.