(ProQuest: ... denotes non-US-ASCII text omitted.)
Recommended by Alexander P. Seyranian
Department of Engineering Sciences, University of Patras, 265 04 Patras, Greece
Received 3 June 2010; Accepted 25 January 2011
1. 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 case. It is to this end that magnetic fluids [1], such as water, hydrocarbon, ester, fluorocarbon, and others, comprise a novel class of engineering materials, where their specific magnetic and hydrodynamic feature promotes heat and mass transfer, reaction rates, blood or other biological flows, and so forth, under the action of magnetic fields. More specifically, magnetic fluids can be considered as colloidal suspensions of spherical particles in a conductive liquid, where the assumption of complete isotropic shape of the magnetic particles is valid due to their small size. In this case, the so-called ferromagnetic particles follow the Brownian motion and behave as rigid magnetic dipoles. Thus, the application of an external magnetic field, apart from the creation of an induced magnetic field of minor significance, will prevent the rotation of each particle, increasing the effective viscosity of the fluid and will cause the appearance of an additional magnetic pressure. The combination of the hydrodynamic flow [2] with the application of an arbitrarily orientated magnetic field allows us to consider the fluid as micropolar [1, 3, 4] and provides us with the appropriate tools for developing and solving such boundary value problems of high mathematical and technical complexity. The most general consideration secures the consistency of many boundary value problems with the principles of both ferrohydrodynamics and magnetohydrodynamics, by including, respectively, both magnetization and electrical conductivity of the fluid within the partial differential equations, which are under analytical or numerical investigation [1-4]. Specifically, ferrohydrodynamics is concerned with the mechanics of the motion of the so-called micropolar fluids that is influenced by strong forces of magnetic polarization, introducing a magnetic stress tensor [5] into the momentum equation. On the other hand, magnetohydrodynamics deals with the current distribution of the electrically conductive carrier liquids. It is worth mentioning that nowadays a very big part of the scientific research is dedicated to the biomagnetic fluid dynamics, which combines the existence of the magnetization of the biofluid with its conductivity; see, for instance, blood flow problems. The system of hydrodynamic equations that governs the magnetic fluids flow [1-4] is constituted by the equations of continuity, of momentum, of Maxwell, of energy (when we deal with the thermomechanics of a magnetic fluid [3]), of angular momentum, and of magnetization, which are all coupled with each other. Those partial differential equations are expressed in terms of the velocity field, the total pressure field, the variation of the temperature, the magnetic field, the conduction current density, and the several specific hydrodynamic or magnetic parameters, which characterize the ferrofluid and the motion itself. The appropriate boundary conditions, depending on the particular physical problem, are adequate for the completeness of a well-posed boundary value problem.
Recently, a general three-dimensional theoretical model that conforms to physical reality and at the same time permits the analytical investigation of the aforementioned partial differential equations has been developed by Hatzikonstantinou and Vafeas [6]. In this novel model the authors focused both on ferrohydrodynamic and on magnetohydrodynamic problems of technological or biomechanical interest, which can be applied to various micropolar magnetic fluids, by constructing purely analytical solutions and using the minimum of the necessary assumptions for their model. This work [6] involves the study of the incompressible flow of a Newtonian carrier magnetic liquid, which contains a small concentration of magnetic particles under the effect of an arbitrarily orientated applied magnetic field. There, it has been presented the development of a new dyadic, as well as vector, expression for the momentum equation, which incorporates an explicit theoretical expression for the extra viscosity. The additional magnetic pressure generated by the effect of the magnetic field and a general three-dimensional analytical expression for the magnetization and the Lorentz forces are also included. However, many of the industrial 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 closed-type 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.
2. 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(...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=x1x...1 +x2x...2 +x3x...3 expressed via the Cartesian basis x...i , i=1,2,3 , in Cartesian coordinates (x1 ,x2 ,x3 ) , where this dependence will be omitted for convenience in writing. This coordinate system is placed in such a way so as the x1 -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,[straight phi],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 rp 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+ρgx3 with the applied magnetic field H of measure H=|H| . It is noted that p is the thermodynamic pressure, g=-gx...3 defines the acceleration of the gravity of measure g , and ρgx3 refers to a hydrostatic pressure force, which corresponds to a height of reference in the x3 -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] [figure omitted; refer to PDF] and the continuity equation [figure omitted; refer to PDF] 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 =rp2ρp /15η0 is the relaxation time of particle rotation (η0 corresponding to the rotational viscosity), I=8πrp5ρ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πηrp3 /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 M0 and in terms of the Langevin function L(ξ) by the equation [figure omitted; refer to PDF] where m specifies the magnetic moment of a particle, while M0 =M0 (H) . It is obvious that, in some cases, we have to take into consideration the magnetization of the particle itself Mp , which is connected with the saturation magnetization Ms by the relation [figure omitted; refer to PDF] where [varphi] is the volumetric concentration of particles. On the other hand, as ξ[arrow right]0 , then L(ξ)[arrow right]0 , whilst as ξ[arrow right]+∞ , it is easily verified that L(ξ)[arrow right]1 . Hence, M0 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 [figure omitted; refer to PDF] 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 [figure omitted; refer to PDF] 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 [figure omitted; refer to PDF] 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 [figure omitted; refer to PDF] 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=[varphi]=0 ) or if there is no magnetic field (H=0 ), then from (2.3) we have M0 =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 (M0 [congruent with]Ms =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 M0 is given by (2.3); that is, [figure omitted; refer to PDF] [figure omitted; refer to PDF] where P0 is a constant reference pressure and Φ, Ψ are potentials that satisfy [figure omitted; refer to PDF] 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=[varphi]=0 ) or if there is no magnetic field (H=0 ), we obtain M0 =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)[arrow right][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 [figure omitted; refer to PDF] 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 [figure omitted; refer to PDF] or [figure omitted; refer to PDF] without loss of generality, since the arbitrary constant of the integration can be embodied inside the arbitrary constant pressure of reference P0 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 [figure omitted; refer to PDF] Combining relationships (2.14), (2.15) and after simple analytical manipulations, we take [figure omitted; refer to PDF] 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 M0 , which is provided in relation (2.3). It is noted that, when M0 [congruent with]Ms =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)=μ0Ms 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 M0 is approximated by the constant value of the saturation magnetization Ms , then the last term of (2.12) becomes h(H)∇(ln M0 )=μ0M0 H∇(ln M0 )=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: [figure omitted; refer to PDF] 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 [figure omitted; refer to PDF] which leads us to the corresponding momentum equation for magnetic Stokes flow of [6] for very high values of ξ , that is, [figure omitted; refer to PDF]
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 [figure omitted; refer to PDF] where this situation provides us with the constant equilibrium magnetization [figure omitted; refer to PDF] and the constant dimensionless parameter [figure omitted; refer to PDF] Therefore, applying (2.21) and (2.22) into the magnetic Stokes equations (2.12) and (2.2), we recover the corresponding differential equations [figure omitted; refer to PDF] since hc ∇(ln M0,c )=∇hc =0 and [figure omitted; refer to PDF] The particular case of the application of a constant three-dimensional magnetic field leads via relation (2.22) to ∇[ln χc ]=0 and ∇h...(Hc )=0[implies]h...(Hc )≡h...c . Therefore, our improved differential representation (2.9)-(2.11) of the solutions of (2.23) and (2.24) takes the form [figure omitted; refer to PDF] [figure omitted; refer to PDF] where Pc ≡P0 +[hc -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, [figure omitted; refer to PDF] 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: [figure omitted; refer to PDF] [figure omitted; refer to PDF] [figure omitted; refer to PDF] 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,[straight phi]) , 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[congruent with]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, [figure omitted; refer to PDF]
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 [figure omitted; refer to PDF] respectively, where -A≡dPP /dz<0 is the constant pressure gradient at the z -direction of the magnetic Poiseuille flow (vP ,PP ) . Immediate integration gives [figure omitted; refer to PDF] We point out that the constant of the integration in order to find PP 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: [figure omitted; refer to PDF] where the pair of flow fields (vP ,PP ) refers to the Poiseuille flow, while the general velocity field vg and the general total pressure field Pg must satisfy the differential representation (2.25), (2.26) with (2.27) of the solutions of (2.23) and (2.24), [figure omitted; refer to PDF] in view of (2.20)-(2.22) for χc . 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 (vP ,PP ) given in (2.32) and (2.33). Thus, decomposition (2.34) holds true, and the boundary condition (2.32) is substituted by [figure omitted; refer to PDF] a condition that must be satisfied automatically when the 3D flow fields (2.34) are calculated.
Our goal is to solve the aforementioned boundary value problem (2.23)-(2.24) with the boundary conditions (2.28)-(2.30) supplemented with (2.31), using the decomposition (2.34) with (2.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.
3. 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 [figure omitted; refer to PDF] in terms of the functions (χc provided by (2.22)) [figure omitted; refer to PDF] 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 Nn (μr) must be excluded from any general harmonic expansion of the type (A.30). The constant parameter μ∈... 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 [figure omitted; refer to PDF] where, for n≥0, μ∈... , and i=1,2,3,4 , the coefficients enμ,(i) =anμ,(i)x...1 +bnμ,(i)x...2 +cnμ,(i)x...3 and dnμ,(i) 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 H -functions, usually named as eigenfunctions, for every kind i=1,2,3,4 are [figure omitted; refer to PDF] [figure omitted; refer to PDF] [figure omitted; refer to PDF] [figure omitted; refer to PDF] 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 [figure omitted; refer to PDF] for the velocity field, while for the total pressure field we obtain [figure omitted; refer to PDF] in the cylinder's prescribed dimensions V(...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 enμ,(i) and dnμ,(i) for n≥0 , μ∈... , 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 Δ(∇Hnμ,(i) )=∇(ΔHnμ,(i) )=0, n≥0, μ∈... , i=1,2,3,4 , which means that ∇Hnμ,(i) belong to the subspace produced by Hnμ,(i) , 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..., [straight phi]..., z... (see also (A.13)). Hence, in view of (A.13), (A.15) and the basic relations for the trigonometric functions (A.25)-(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 [figure omitted; refer to PDF] or [figure omitted; refer to PDF] for every n≥0 and μ∈... , where the prime at the Bessel functions denotes derivation with respect to the variable μr∈... . In the same way we calculate the rest three gradients, that is, [figure omitted; refer to PDF] for every n≥0 and μ∈... , where it is obvious that Δ(∇Hnμ,(i) )=∇(ΔHnμ,(i) )=0 . Here, we must add that, in order to avoid negative values of n within the above relationships, we impose [figure omitted; refer to PDF] Expressions (3.11) and (3.12) can be collected in the general formulae [figure omitted; refer to PDF] for every n≥0 , while ai =i+(-1)i+1 , bi =i+2ci for i=1,2,3,4 and ci =1, i=1,2 or ci =-1, i=3,4 . Relation (A.12) and Cartesian definition of enμ,(i) result in the product [figure omitted; refer to PDF] which by virtue of (3.14) incorporates with (3.8) to the following velocity field: [figure omitted; refer to PDF] At this point we manipulate properly the H -functions Hn+1μ,(i) , Hn-1μ,(i) , Hn+1μ,(bi ) , and Hn-1μ,(bi ) by a certain readjustment of the index n at the series to rewrite the velocity (3.16) as [figure omitted; refer to PDF] 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 anμ,(i) , bnμ,(i) , cnμ,(i) , and dnμ,(i) for n≥0, μ∈... , and i=1,2,3,4 , which are grouped via the simple relations [figure omitted; refer to PDF] [figure omitted; refer to PDF] [figure omitted; refer to PDF] [figure omitted; refer to PDF]
Similarly, formulae (3.14) and the group of constant coefficients (3.18)-(3.21) are interrelated, and the total pressure (3.9) becomes [figure omitted; refer to PDF] or since x...i ·x...j =δij (δij being the Kronecker delta) and with a proper readjustment of the index n at the series for the H -functions Hn+1μ,(i) ,Hn-1μ,(i) ,Hn+1μ,(bi ) , and Hn-1μ,(bi ) as previously mentioned, relation (3.22) yields [figure omitted; refer to PDF] for the total 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 anμ,(i) ,bnμ,(i) ,cnμ,(i) , and dnμ,(i) for n≥0 and i=1,2,3,4 , as well as the parameter μ∈... . During this process, we must keep in mind that ai =i+(-1)i+1 , bi =i+2ci for i=1,2,3,4 and ci =1, i=1,2 or ci =-1, i=3,4 .
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 [figure omitted; refer to PDF] since r...·z...=r...·x...1 =0 , while [straight phi]...·v=0 at z=0 renders [figure omitted; refer to PDF] since [straight phi]...·z...=[straight phi]...·x...1 =0 . We notice that both conditions (3.24) and (3.25) contain the H -eigenfunctions 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 sinn[straight phi] , n≥1 , and via the definition for the Bessel functions (A.19), we obtain the following relations for the constant coefficients: [figure omitted; refer to PDF] while for n=0 we can admit without loss of generality (sin0[straight phi]=0 ) that d0μ,(1) =b0μ,(1) =c0μ,(3) =0 , μ∈... . On the other hand, an easy manipulation of (3.26) reveals the nihilism of the constant coefficients dnμ,(1) ,bnμ,(1) , and cnμ,(3) for n≥1 and μ∈... . Recapitulating, [figure omitted; refer to PDF] whilst orthogonality of cos n[straight phi] , n≥0 , on (3.24) is provided with similar relations to (3.26), leading to [figure omitted; refer to PDF] 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)) [figure omitted; refer to PDF] and via (3.18)-(3.21) [figure omitted; refer to PDF] 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) and (3.23), taking into account the constant parameters ai , bi , ci for i=1,2,3,4 with the definitions for the H -functions (3.4)-(3.7), and we recover the velocity field [figure omitted; refer to PDF] while for the total pressure field we arrive at [figure omitted; refer to PDF] in terms of the constant coefficients (3.18)-(3.21) and the factor (2.22), where u(r), p(z) are given via relation (3.2).
Next, we 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 [figure omitted; refer to PDF] where we have used the fact that ∂r.../∂z=∂[straight phi].../∂z=0 (see (A.13)), while for reasons of convenience to our forthcoming manipulations we incorporate the two relations of (3.33) into [figure omitted; refer to PDF] 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 [figure omitted; refer to PDF] 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 sinn[straight phi] , n≥0 and a certain readjustment of the index n at the series infer the following formulae: [figure omitted; refer to PDF] Orthogonality arguments on the functions cos n[straight phi] , n≥0 with similar readjustment of n imply [figure omitted; refer to PDF] 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 [figure omitted; refer to PDF] 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 [figure omitted; refer to PDF] 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 [figure omitted; refer to PDF] [figure omitted; refer to PDF] respectively, where in view of definitions (3.19)-(3.20) and after some trivial processing, expressions (3.40) and (3.41) take the form [figure omitted; refer to PDF] [figure omitted; refer to PDF] 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 [figure omitted; refer to PDF] 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 [figure omitted; refer to PDF] 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 Jn (μr)sinn[straight phi]cosh (μz) and Jn (μr)cos n[straight phi]cosh (μz) for every n≥0 and μ∈... . Since they appear nowhere else in the velocity or the total pressure, we can group them as [figure omitted; refer to PDF] [figure omitted; refer to PDF] 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 enμ,(1) , enμ,(3) for n≥0 and μ∈... provided by (3.46), (3.47), respectively, is furnished by [figure omitted; refer to PDF] while the total pressure field (3.32) with the aid of relations (3.45) becomes [figure omitted; refer to PDF] whereas the A 's 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 anμ,(1) ,anμ,(2) ,anμ,(3) ,anμ,(4) ,enμ,(1) , and enμ,(3) for every n≥0 , μ∈... 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[congruent with]p , in order to take [figure omitted; refer to PDF] 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: [figure omitted; refer to PDF] Orthogonality arguments with respect to the eigenfunctions Jn (μr)sinn[straight phi] and Jn (μr)cos n[straight phi] for n≥0 and μ∈... yield [figure omitted; refer to PDF] [figure omitted; refer to PDF] for n≥0 and μ∈... , where we have expressed the e 's constant coefficients via the a 's constant coefficients into the flow fields (3.48) and (3.49).
Consequently, we are left with four sets of constant coefficients; those are anμ,(1) ,anμ,(2) ,anμ,(3) , and anμ,(4) for every n≥0 , μ∈... , which must be determined from the continuity restraint (2.31) and the second part of the inlet boundary condition (2.29). The constant parameter μ∈... 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 [figure omitted; refer to PDF] whereas λμ ∈... is the new parameter to be calculated for μ∈... , 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 [figure omitted; refer to PDF] for the velocity field and [figure omitted; refer to PDF] 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..., [straight phi]..., 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 [figure omitted; refer to PDF] where by definition [figure omitted; refer to PDF] for λμ ∈... and μ∈... . 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.58), we conclude that [figure omitted; refer to PDF] where we have used the Bessel ordinary differential equation (A.18) to evaluate [figure omitted; refer to PDF] Simple orthogonality arguments of the products μJn (μr)(anμ,(1) sinn[straight phi]+anμ,(3) cos n[straight phi]) for every n≥0 and μ∈... on (3.59) provides us with the constant parameter λμ as [figure omitted; refer to PDF] Actually, the flow fields retain only the parameter μ , and their equivalent form with respect to (3.61) is [figure omitted; refer to PDF] from velocity field (3.55) and [figure omitted; refer to PDF] from total pressure field (3.56), respectively, where the prime appearing at gμ[variant prime] (μz) denotes derivation with respect to the variable μz∈... . The new functions of expansions (3.62) and (3.63) are defined as [figure omitted; refer to PDF] where [figure omitted; refer to PDF] 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 , [figure omitted; refer to PDF] and provided that the remaining unknown constant coefficients anμ,(1) and anμ,(3) for n≥0 , as well as the unknown parameter μ∈... , are bounded, the velocity (3.62) and total pressure (3.63) asymptotically lead to [figure omitted; refer to PDF] 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, [figure omitted; refer to PDF] 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, [figure omitted; refer to PDF] a condition that provides us with the value of the parameter μ as discrete eigenvalues via [figure omitted; refer to PDF] where ρnm is the m -root (m≥1 ) of order n≥0 of the Bessel functions (Jn (ρnm )=0 ). Consequently, formula (3.70) determines the parameter μ∈... , 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 μnm -symbolism given now by (3.70), where μnm is the m -parameter (m≥1 ) of order n≥0 . Introducing the new constant coefficients that have to be calculated as [figure omitted; refer to PDF] and instead of (3.18) their equivalent forms [figure omitted; refer to PDF] for every n≥0, m≥1 , and μ≡μnm ∈... , the velocity field (3.62) is given via the eigenvalue-type expansion [figure omitted; refer to PDF] while the total pressure field (3.63), under the aforementioned readjustment, gives rise to [figure omitted; refer to PDF] 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: [figure omitted; refer to PDF] where [figure omitted; refer to PDF] Here, μnm ,χc and u(r), p(z) satisfy relations (3.70), (2.22), and (3.2), respectively, while the remaining constant coefficients anm and bnm 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,[straight phi]) for z=0 at the velocity (3.73) (note that x...1 =z... ) and use of (3.70) result in [figure omitted; refer to PDF] By virtue of the orthogonality relation (A.32) and the well-known orthogonality of the trigonometric functions, relation (3.77) reveals that [figure omitted; refer to PDF] for n≥0 and m≥1 , where [straight epsilon]n =2 for n≥1 , [straight epsilon]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[arrow right]+∞ . 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 b0m for m≥1 and the magnetic flow becomes axisymmetric. Moreover, in this case, if v(r)=u(r) , then (3.78) shows that anm =bnm =0 for n≥0 and m≥1 ; thus, the magnetic flow (3.73) becomes the magnetic Poiseuille flow v=u(r)z... .
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..., [straight phi]... , 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 [figure omitted; refer to PDF] [figure omitted; refer to PDF] where μnm =ρnm /α (ρnm 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) [figure omitted; refer to PDF] and of the total pressure field (3.74) [figure omitted; refer to PDF]
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,[straight phi]) 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 [figure omitted; refer to PDF] assume the expansions (3.81) and (3.82), respectively, where anm , bnm for n≥0 and m≥1 have been calculated through relations (3.78). The functions hnm (μnm z),gnm (μnm z),jn,±m (μnm r,μnm z) , and fnm (n[straight phi]) for n≥0 and m≥1 are given by (3.75), (3.76), (3.79), and (3.80), respectively, while μnm =ρnm /α , χc and u(r), p(z) satisfy relationships (3.70), (2.22), and (3.2), respectively. Since ρnm is the m -root (m≥1 ) of order n≥0 of the Bessel functions, that is, Jn (ρnm )=0 , then obviously jn,±m (ρnm ,μnm z)=0 . Finally, any information needed beyond this section, which concerns the Bessel functions Jn (μnm r) or the trigonometric and hyperbolic functions, is collected into the appendix.
4. 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 [figure omitted; refer to PDF] where Λ is the cross section at the inlet. Consequently, the mass flow rate at the inlet of the duct is [figure omitted; refer to PDF] If we consider the limiting case of a semi-infinite duct with L[arrow right]+∞ 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 [figure omitted; refer to PDF] which yields [figure omitted; refer to PDF] In order for our results to have general applicability, we introduce the nondimensional variables [figure omitted; refer to PDF] 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 [figure omitted; refer to PDF] which is very close to a uniform inlet velocity. Figure 1 shows contour plots of the dimensionless axial velocity component u...z ≡z...·v... , the dimensionless radial velocity component 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.
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.
(a) [figure omitted; refer to PDF]
(b) [figure omitted; refer to PDF]
(c) [figure omitted; refer to PDF]
(d) [figure omitted; refer to PDF]
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 cross-section 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.
Figure 2: Diagrams of the dimensionless axial velocity distribution u...z along the radial direction for different axial lengths z... .
[figure omitted; refer to PDF]
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 [figure omitted; refer to PDF] 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...[straight phi] ≡[straight phi]...·v... , and the dimensionless total pressure P... , at several downstream locations.
Contour plots of the dimensionless axial velocity component u...z , the dimensionless radial velocity component u...r , the dimensionless angular velocity component u...[straight phi] , and the dimensionless total pressure P... at several downstream locations.
(a) [figure omitted; refer to PDF]
(b) [figure omitted; refer to PDF]
(c) [figure omitted; refer to PDF]
(d) [figure omitted; refer to PDF]
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]: [figure omitted; refer to PDF] where Dh =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 =-αdPP /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 [figure omitted; refer to PDF] 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 fRe=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 [varphi] , which is defined in relation (2.4).
Figure 4: Variation of parameter χc with the magnetic field H and the volumetric concentration of magnetic particles [varphi] , where the plot corresponds to fully developed flow of a micropolar fluid.
[figure omitted; refer to PDF]
The plots correspond to constant temperature T=300 K, magnetic dipole strength m=2×10-18 A· m2 , viscosity η=10-3 kg/ms, and small magnetic particles of radius rp <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 rp . This can be seen from the relation (2.22) after the substitution of τS ,τB ,I and of M0,c from (2.21). From the plot of Figure 4 it is observed that, independently of the concentration [varphi] , the parameter χc undergoes a steep increase in the interval 103 <H<105 (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 [varphi] increases.
5. 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 three-dimensional 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.
[1] B. Berkovski, V. Bashtovoy Magnetic Fluids and Applications Handbook , Begell House, New York, NY, USA, 1996.
[2] R. W. Fox, A. T. McDonald Introduction to Fluid Mechanics , John Wiley & Sons, New York, NY, USA, 1985.
[3] V. G. Bashotovoy, B. M. Berkovski, A. N. Vislovich Introduction to Thermomechanics of Magnetic Fluids , Hemisphere, New York, NY, USA, 1988.
[4] R. E. Rosensweig Ferrohydrodynamics , Dover, New York, NY, USA, 1997.
[5] L. D. Landau, E. M. Lifshitz, L. P. Pitaevskii Electrodynamics of Continuous Media , Pergamon Press, New York, NY, USA, 1984.
[6] P. M. Hatzikonstantinou, P. Vafeas, "A general theoretical model for the magnetohydrodynamic flow of micropolar magnetic fluids. Application to Stokes flow," Mathematical Methods in the Applied Sciences , vol. 33, no. 2, pp. 233-248, 2010.
[7] J. Happel, H. Brenner Low Reynolds Number Hydrodynamics with Special Applications to Particulate Media , Prentice-Hall, Englewood Cliffs, NJ, USA, 1986.
[8] H. Neuber, "Ein neuer ansatz zur lösung räumblicher probleme der elastizitätstheorie," Journal of Applied Mathematics and Mechanics , vol. 14, pp. 203-212, 1934.
[9] P. M. Naghdi, C. S. Hsu, "On a representation of displacements in linear elasticity in terms of three stress functions," Journal of Mathematical Mechanics , vol. 10, pp. 233-245, 1961.
[10] X. S. Xu, M. Z. Wang, "General complete solutions of the equations of spatial and axisymmetric Stokes flow," The Quarterly Journal of Mechanics and Applied Mathematics , vol. 44, no. 4, pp. 537-548, 1991.
[11] P. D. S. Verma, "On a class of exact solutions to magnetic fluid equations," Journal of Magnetism and Magnetic Materials , vol. 39, no. 1-2, pp. 137-141, 1983.
[12] V. U. K. Sastry, V. R. Mohan Rao, "Micropolar fluid flow due to an oscillating plane subject to rotation," Acta Mechanica , vol. 33, no. 1-2, pp. 45-53, 1979.
[13] C. Calmelet-Eluhu, D. R. Majumdar, "Flow of a micropolar fluid through a circular cylinder subject to longitudinal and torsional oscillations," Mathematical and Computer Modelling , vol. 27, no. 8, pp. 69-78, 1998.
[14] H. C. Weng, C.-K. Chen, M.-H. Chang, "Stability of micropolar fluid flow between concentric rotating cylinders," Journal of Fluid Mechanics , vol. 631, pp. 343-362, 2009.
[15] M. S. Faltas, E. I. Saad, "Stokes flow with slip caused by the axisymmetric motion of a sphere bisected by a free surface bounding a semi-infinite micropolar fluid," International Journal of Engineering Science , vol. 43, no. 11-12, pp. 953-976, 2005.
[16] H. H. Sherif, M. S. Faltas, E. I. Saad, "Slip at the surface of a sphere translating perpendicular to a plane wall in micropolar fluid," Zeitschrift für Angewandte Mathematik und Physik , vol. 59, no. 2, pp. 293-312, 2008.
[17] A. Moosaie, G. Atefi, "Slow motion of a rotating circular cylinder through a micropolar fluid," Archives of Mechanics , vol. 60, no. 3, pp. 199-220, 2008.
[18] P. K. Papadopoulos, E. E. Tzirtzilakis, "Biomagnetic flow in a curved square duct under the influence of an applied magnetic field," Physics of Fluids , vol. 16, no. 8, pp. 2952-2962, 2004.
[19] S. A. Khashan, Y. Haik, "Numerical simulation of biomagnetic fluid downstream an eccentric stenotic orifice," Physics of Fluids , vol. 18, no. 11, 2006.
[20] C. C. Sellers, J. S. Walker, "Liquid-metal flow in an electrically insulated rectangular duct with a non-uniform magnetic field," International Journal of Engineering Science , vol. 37, no. 5, pp. 541-552, 1999.
[21] L. M. Witkowski, J. S. Walker, "Numerical solutions for the liquid-metal flow in a rotating cylinder with a weak transverse magnetic field," Fluid Dynamics Research , vol. 30, no. 3, pp. 127-137, 2002.
[22] P. K. Papadopoulos, P. Vafeas, P. M. Hatzikonstantinou, "Numerical study of a new model for the magnetohydrodynamic flow of micropolar magnetic fluids in straight square ducts," in Proceedings of the 6th International Conference on Engineering Computational Technology, pp. 1-19, Civil-Comp Press, 2008.
[23] P. Vafeas, P. K. Papadopoulos, P. M. Hatzikonstantinou, "Micropolar flow under the effect of a magnetic dipole," in Proceedings of the 6th International Conference of Computational Methods in Sciences and Engineering, vol. 1148, of American Institute of Physics Proceedings, pp. 566-570, 2008.
[24] G. Dassios, M. Hadjinicolaou, A. C. Payatakes, "Generalized eigenfunctions and complete semiseparable solutions for Stokes flow in spheroidal coordinates," Quarterly of Applied Mathematics , vol. 52, no. 1, pp. 157-191, 1994.
[25] P. Vafeas, G. Dassios, "Stokes flow in ellipsoidal geometry," Journal of Mathematical Physics , vol. 47, no. 9, 2006.
[26] L. Martin Witkowski, J. S. Walker, P. Marty, "Nonaxisymmetric flow in a finite-length cylinder with a rotating magnetic field," Physics of Fluids , vol. 11, no. 7, pp. 1821-1826, 1999.
[27] P. Moon, E. Spencer Field Theory Handbook , Springer, Berlin, Germany, 1971.
[28] P. M. Morse, H. Feshbach Methods of Theoretical Physics, Vol I, II , McGraw-Hill, New York, NY, USA, 1953.
[29] E. W. Hobson The Theory of Spherical and Ellipsoidal Harmonics , pp. xi+500, Chelsea, New York, NY, USA, 1965.
[30] P. M. Gresho, "On the theory of semi-implicit projection methods for viscous incompressible flow and its implementation via a finite element method that also introduces a nearly consistent mass matrix. I. Theory," International Journal for Numerical Methods in Fluids , vol. 11, no. 5, pp. 587-620, 1990.
[31] R. K. Ganesh, "Evaluation of the drag force by integrating the energy dissipation rate in Stokes flow for 2D domains using the FEM," International Journal for Numerical Methods in Fluids , vol. 13, no. 5, pp. 557-578, 1991.
[32] C. H. Bruneau, P. Fabrie, "Effective downstream boundary conditions for incompressible Navier-Stokes equations," International Journal for Numerical Methods in Fluids , vol. 19, no. 8, pp. 693-705, 1994.
[33] G. B. Arfken, H. J. Weber Mathematical Methods for Physicists , pp. xiv+1112, Elsevier Academic, New York, NY, USA, 2005., 5th.
[34] R. A. Eubanks, E. Sternberg, "On the completeness of the Boussinesq-Papkovich stress functions," Journal of Rational Mechanics and Analysis , vol. 5, pp. 735-746, 1956.
[35] W. M. Rohsenow, J. P. Hartnett, Y. I. Cho Handbook of Heat Transfer , McGraw-Hill, New York, NY, USA, 1998.
Appendix
Mathematical Material
In the interest of making this paper complete and independent it is necessary to provide some useful information concerning important mathematical tools, which are adequate for both the analytical and the numerical part of the present paper.
Differential Identities
Consequently, we begin by providing important identities [29], which were used extensively during our analytical calculations. Let u,v and f,g denote two scalar and two vector fields, respectively. Then, the basic vector differential identities involve the action of the gradient and the Laplace operators on certain expressions, that is, [figure omitted; refer to PDF] [figure omitted; refer to PDF] [figure omitted; refer to PDF] [figure omitted; refer to PDF] [figure omitted; refer to PDF] [figure omitted; refer to PDF] [figure omitted; refer to PDF] [figure omitted; refer to PDF] [figure omitted; refer to PDF] [figure omitted; refer to PDF] [figure omitted; refer to PDF] where the symbols "[ecedil]7; ", ": ", and "T " denote juxtaposition, double inner product, and transposition, respectively. The proofs of identities (A.1)-(A.11) are based on classical analysis by expanding the vectors f and g in Cartesian coordinates.
Circular Cylindrical Coordinate System
In terms of the three variables r∈[0,+∞) , [straight phi]∈[0,2π) , and z∈(-∞,+∞) , we define the implemented to our work circular cylindrical coordinate system [27] via the Cartesian basis as follows: [figure omitted; refer to PDF] where r is the position vector, while the coordinate vectors of this system r..., [straight phi]..., z... assume the forms [figure omitted; refer to PDF] or their inverse [figure omitted; refer to PDF] The gradient ∇ and the Laplacian Δ are provided by the expressions [figure omitted; refer to PDF] and then the forthcoming relationships are easily verified using the vector relations (A.12), (A.13), and (A.15), that is, [figure omitted; refer to PDF] as well as [figure omitted; refer to PDF] 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 x1 -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,[straight phi],z) right-handed system.
Bessel and Neumann Special Functions
In the sequel, we define the special functions used in this project. The Bessel functions Jn (x) and the Neumann functions Nn (x) [29] of order n≥0 are linear independent solutions of the Bessel differential equation [figure omitted; refer to PDF] for n=0,1,2,... , which is valid for every x∈... . These functions are defined as follows: [figure omitted; refer to PDF] [figure omitted; refer to PDF] for n≥0 and ν∈... , where the Bessel functions Jν (x) for real ν are furnished by the formula [figure omitted; refer to PDF] via the well-known Gamma function Γ(κ+ν+1) [29] for every κ≥0 and ν∈... . In our case, where ν=n=0,1,2,... , it is easily confirmed that J-n (x)=(-1)nJn (x) and N-n (x)=(-1)nNn (x) for x∈... . The limiting cases as x[arrow right]0 , where Jn (x) converge, or as x[arrow right]+∞ , where Nn (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 [figure omitted; refer to PDF] [figure omitted; refer to PDF] defined as Zn (x)≡Jn (x) or Zn (x)≡Nn (x) for n≥0 and x∈... . 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∈... , with μ∈... 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 μnm ∈... for n≥0 and m≥1 , then x is substituted by x=μr≡μnm r∈... , whereas μ≡μnm is the m -parameter (m≥1 ) of order n≥0 . Thus, both the Bessel and the Neumann functions satisfy the orthogonality relation [figure omitted; refer to PDF] for n≥0 and m,m[variant prime] ≥1 , with δmm[variant prime] being the Kronecker delta and a,b∈... being constants.
Trigonometric and Hyperbolic Functions
As far as the trigonometric functions sinn[straight phi] and cos n[straight phi] are concerned for n≥0 , which are used in this paper, the following expressions hold true: [figure omitted; refer to PDF] [figure omitted; refer to PDF] [figure omitted; refer to PDF] [figure omitted; refer to PDF] where [straight phi]∈[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 sinn[straight phi] and cos n[straight phi] 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, [figure omitted; refer to PDF] 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,[straight phi],z)=0 ) in the circular cylindrical coordinate system is written as [figure omitted; refer to PDF] where anμ , bnμ , cnμ , dnμ , enμ , and fnμ 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 α∈... , where the boundary conditions require nihilism of the Bessel functions on the surface of the cylinder, then the regularity of the solution demands Zn (μr)≡Jn (μr) and the boundary conditions provide us with the discrete values of the parameter as μ≡μnm =ρnm /α , where ρnm ∈... is the m -root (m≥1 ) of order n≥0 of the Bessel functions and the eigenvalue now μnm is the m -parameter (m≥1 ) of order n≥0 . Hence, for this particular problem of determining parameter μ , expansion (A.30) becomes [figure omitted; refer to PDF] where anm , bnm , cnm , dnm , enm , and fnm for n≥0 , m≥1 are the new constants, while we have a=0 , b=α and relation (A.24) takes the form [figure omitted; refer to PDF] 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 Jn (μr)≡Jn (μnm r) survive inside the expansion (A.30) or (A.31), and, thus, this is reflected to the vanishing of the corresponding constants, that is, bnμ =0 for n≥0 or bnm =0 for n≥0 and m≥1 , respectively.
You have requested "on-the-fly" machine translation of selected content from our databases. This functionality is provided solely for your convenience and is in no way intended to replace human translation. Show full disclaimer
Neither ProQuest nor its licensors make any representations or warranties with respect to the translations. The translations are automatically generated "AS IS" and "AS AVAILABLE" and are not retained in our systems. PROQUEST AND ITS LICENSORS SPECIFICALLY DISCLAIM ANY AND ALL EXPRESS OR IMPLIED WARRANTIES, INCLUDING WITHOUT LIMITATION, ANY WARRANTIES FOR AVAILABILITY, ACCURACY, TIMELINESS, COMPLETENESS, NON-INFRINGMENT, MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE. Your use of the translations is subject to all use restrictions contained in your Electronic Products License Agreement and by using the translation functionality you agree to forgo any and all claims against ProQuest or its licensors for your use of the translation functionality and any output derived there from. Hide full disclaimer
Copyright © 2011 Panayiotis Vafeas et al. Panayiotis Vafeas 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
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 non-axisymmetric geometries, which provides the velocity and total pressure fields in terms of easy-to-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.
You have requested "on-the-fly" machine translation of selected content from our databases. This functionality is provided solely for your convenience and is in no way intended to replace human translation. Show full disclaimer
Neither ProQuest nor its licensors make any representations or warranties with respect to the translations. The translations are automatically generated "AS IS" and "AS AVAILABLE" and are not retained in our systems. PROQUEST AND ITS LICENSORS SPECIFICALLY DISCLAIM ANY AND ALL EXPRESS OR IMPLIED WARRANTIES, INCLUDING WITHOUT LIMITATION, ANY WARRANTIES FOR AVAILABILITY, ACCURACY, TIMELINESS, COMPLETENESS, NON-INFRINGMENT, MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE. Your use of the translations is subject to all use restrictions contained in your Electronic Products License Agreement and by using the translation functionality you agree to forgo any and all claims against ProQuest or its licensors for your use of the translation functionality and any output derived there from. Hide full disclaimer