Mathematical Problems in Engineering
Volume 2018 (2018), Article ID 6759876, 16 pages
https://doi.org/10.1155/2018/6759876

Research Article

Asymptotic Approximation of the Nonsteady Micropolar Fluid Flow through a Circular Pipe

Igor Pažanin and Marko Radulović

Department of Mathematics, Faculty of Science, University of Zagreb, Zagreb, Croatia

Correspondence should be addressed to Igor Pažanin; pazanin@math.hr

Received 27 September 2017; Revised 2 January 2018; Accepted 29 January 2018; Published 11 April 2018

Academic Editor: Michael Vynnycky

Copyright © 2018 Igor Pažanin and Marko Radulović. 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.

We study the nonsteady flow of a micropolar fluid through a thin cylindrical pipe. The asymptotic behaviour of the flow is found via asymptotic analysis with respect to the small parameter , representing the pipe’s thickness. The asymptotic approximation is derived in the form of the explicit formulae for the fluid velocity and microrotation. We also provide the numerical examples in order to visualize the effects of the micropolar nature of the fluid. The illustrations indicate the influence of the micropolarity on the effective flow of the fluid in the whole domain. In particular, those effects are most clearly observed for the velocity approximation near the boundary of the domain.

1. Introduction

The theory of micropolar fluids was first introduced in 60s in the papers by Aero et al. [1] and Eringen [2]. Since then, it has been the case of study and interest of both the engineering and mathematical community. The importance of the model lays in the fact that it takes into account the microstructure of the fluid. The fluid particles with complex shapes of numerous fluids (e.g., liquid crystals, muddy fluids, polymeric suspensions, animal blood, and even water at small scale) exhibit microscopical effects (rotation, shrinking) which cannot be captured by the classical Navier–Stokes model. In order to overcome this limitation, that is, to describe the rotation of the particles, independently of the fluid as a whole, we introduce a new vector field, the angular velocity field of the particles. As a consequence, we obtain a coupled system of equations, with four new viscosities introduced. This non-Newtonian mathematical model obtained in such way describes the behaviour of various real fluids better than the classical Newtonian model. The monograph [3] provides a detailed overview of the mathematical theory underlying the micropolar fluid model.

Due to its practical importance in industrial and engineering applications, micropolar fluid flows have been extensively investigated, in particular in the last decade (see, e.g., [4–9]). In the present paper, we aim to study a time-dependent flow of a micropolar fluid through a thin pipe with constant circular cross-section. Traditionally, in order to describe the geometry of such pipes (naturally appearing in numerous applications), a small parameter is introduced, representing the pipe’s thickness. We start from the assumption that the solution of the governing problem has the so-called micropolar Poiseuille form, and then, following recent results obtained for nonsteady Newtonian flows [10–12], we separate the equations by linearity. Employing an asymptotic expansion of the solution in powers of , we simultaneously solve the obtained equations. Assuming that the external forces depend only on the time variable, we manage to analytically construct the asymptotic solution of the flow in the explicit form. It should be emphasized that deriving the explicit expressions for the fluid velocity and microrotation could be particularly important with regard to numerical simulations.

The paper is organized as follows: in Section 2, we describe the pipe’s geometry and the governing micropolar system of equations. Supposing the micropolar Poiseuille form of the solution, we split the problem into the micropolar heat and micropolar inverse problem. In Section 3, we study the micropolar heat problem and construct the asymptotic expansion of the solution up to the order . In Section 4, we do the same for the micropolar inverse problem, where we additionally have to deal with the pressure. In Section 5, we collect the obtained results and write the asymptotic approximation of the considered problem. Finally, in Section 6, we present some numerical examples and visualize the asymptotic solution.

To conclude, let us provide some more bibliographic remarks on the subject. Steady-state flow of a micropolar fluid through thin domains has been extensively discussed throughout the literature. In particular, this work is a natural continuation of the first author’s previous work on the stationary pipe flows; see [13–16]. However, to our knowledge, the results on time-dependent flow have been reported only for simplified settings in which the microrotation is taken to be a scalar function. We refer the reader to [17–20] which merit careful reading. In the present paper we do not impose such constraint leading to a more realistic flow. For that reason, we believe that the obtained result could prove useful in the engineering practice dealing with pipe flows of micropolar fluid. Last but not least, it should be mentioned that theoretical error analysis providing the justification of the asymptotic approximation formally derived in this paper will be brought in the forthcoming paper by the authors of the present paper.

2. Micropolar Equations

Let be a small positive parameter. We consider a thin straight pipe

where

is the circular cross-section of the pipe with constant diameter . We denote by the Cartesian coordinates, with being the direction coinciding with the axis of the pipe (see Figure 1). Let us introduce the following system of equations describing a nonsteady micropolar fluid flow:

Here stands for the velocity field, is the pressure, and is the angular velocity of rotation of the fluid particles (the microrotation field).

Figure 1: The considered flow problem.

To complete the problem, we impose the following boundary and initial conditions:

along with the flux condition

The positive constants are the Newtonian viscosity and the microrotation viscosity , while , , and represent the coefficients of the angular viscosities. To simplify the notation, we abbreviate . The external sources of linear and angular momentum are given by the functions and , respectively.

As indicated in the Introduction, we assume that the solution of the problem (3) has the micropolar Poiseuille form; namely,

where is an arbitrary function of . We also assume that , with , , , , , and being independent of the longitudinal variable .

Substituting the micropolar Poiseuille solution (6) into system (3), we obtain the following initial boundary value problem for :

where . The corresponding boundary conditions read

with given flux rate:

Following the idea from [10–12], we represent the solution as

Here is the solution of the initial boundary value problem for the micropolar heat equation:

On the other hand, denotes the solution of the micropolar inverse problem endowed with the homogeneous initial condition and the given flux rate on the cross-section :

where

From (9) we deduce

Taking into account the problem under consideration, it is plausible to consider the following scaling with respect to small parameter :

In view of that, we expand

and consider problem (7)–(9) with

and defined by (16g).

3. Micropolar Heat Problem

To perform the asymptotic analysis, we first need to rescale the domain, that is, to write the problem on instead of . We introduce the change of variables and obtain the following system of equations deduced from (11):

where . The boundary and initial conditions are the following:

We rewrite problem (18) as

Now, we construct the formal asymptotic expansion of the solution in powers of small parameter up to in the following way:

Due to the small thickness of the pipe, it is reasonable to assume that the functions , , are independent of the cross-section variables . Consequently, we are going to be in position to explicitly compute both zero-order approximation and the correctors.

3.1. Zero-Order Approximation

Plugging (21) and (16d)–(16f) into system (20) and collecting the terms of order , we get

Note that the problems for the velocity and microrotation are decoupled at this particular stage. Equation (22a) with the boundary condition (22d) can be solved by taking

where denotes the solution of the auxiliary problem posed on :

Since we assumed the pipe has a circular cross-section, namely,

we can pass to polar coordinates and compute explicitly from (24):

The explicit expression for the zero-order velocity approximation now reads

Similarly, it can be straightforwardly verified that problem (22b), (22c) with the boundary conditions (22e)-(22f) for microrotation will be satisfied for

3.2. First-Order Corrector

Now, we compute the first-order corrector . Inserting (21) and (16d)–(16f) into the system of equations (20), after collecting the terms of order , we obtain

The system is not decoupled anymore, so the effects of the microstructure on the fluid velocity will be present. Using expressions (28a) and (28b) for the microrotation zero-order approximation , we get

Let us introduce , as the solutions of the following problems:

Taking into account (25) and passing to the polar coordinates easily give

We seek for the solution of system (30) in the form

leading to

Now we compute the corrector for the microrotation. Plugging expression (27) for into (29) yields

Similarly, as for , we obtain explicit expression for ; namely,

3.3. Second-Order Corrector

To capture the effects of the time derivative as well, we have to look for the second-order corrector. Substituting (21) and (16d)–(16f) into system (20) and collecting the terms of order yield

Taking into account (27) and (36), we get the following problem for :

We rewrite

and introduce as the solution of the following problem:

Passing to the polar coordinates provides

We seek for the solution of (39) in the form

implying

Plugging (28a)-(28b) and (34) into (37b) and (37c), we get the following system for the second-order microrotation corrector :

After rewriting, we obtain

Tedious but straightforward calculation gives

where

where , , , , , , , , , .

4. Micropolar Inverse Problem

Now, we turn our attention to micropolar inverse problem (12). First, we notice that

Denoting

we have the compatibility condition (due to (14))

or

Performing the change of variables and putting , from (12) we deduce

where . In the sequel, we look for the asymptotic approximation in the form

4.1. Zero-Order Approximation

For the zero-order approximation , we obtain the following problem:

where is given by

Equation (54a) with the boundary condition (54d) and the flux rate (54g) can be solved by taking

Here we choose to satisfy the flux condition:

where

We now look at the system (54b), (54c), (54e)-(54f) and observe it has the same structure as the system (22b), (22c), (22e)-(22f) for the zero-order approximation of the micropolar heat problem with (see Section 3.1). Thus, we conclude

4.2. First-Order Corrector

The first-order corrector is given by the following problem:

where is given by

We notice that the system is not decoupled anymore. Inserting the expression for given by (56) and (59), we obtain

As above, (60a) with the boundary condition (60d) and flux rate (60g) can be solved by taking

with obeying the flux condition:

with

One can easily verify that the system (62b), (62c) for the microrotation corrector with the boundary conditions (62e)-(62f) is satisfied for

4.3. Second-Order Corrector

To conclude the analysis for the micropolar inverse problem (52), we identify the terms of order to get the system for the second-order corrector . It reads

where is given by

Employing the expressions for the zero-order approximation given by (56) and (59) and for the first-order corrector given by (63) and (66), from (67) we get the following:

We can now explicitly compute the velocity corrector as

where we choose to satisfy the flux condition:

Similarly, we can explicitly compute the corrector for the microrotation:

5. Asymptotic Approximation

The asymptotic approximation related to the micropolar heat problem (11) has the following form:

Note we have calculated all the terms explicitly; that is, , are given by (27) and (28a)-(28b), (34), (36), (43), and (46)-(47).

For the micropolar inverse problem (12), we obtain the following:

where , are given by the explicit formulae (56)–(59), (63)–(66), and (70)–(72).

Collecting the above results, we arrive at the asymptotic approximation for problem (7)–(9):

It should be noted that the above asymptotic approximation was computed to satisfy (7), the boundary conditions (8a)–(8c), and the flux condition (9). The initial conditions (8d)–(8f) were not taken into account in the process due to the fact that the time derivative appears only in the system for the second-order corrector (namely, as the time derivative of the zero-order approximation for the velocity and microrotation which is a known function on the right-hand side). Thus, taking into account the initial conditions while computing the correctors would yield overdetermined systems for the unknown terms in the asymptotic expansion. This essentially means that a boundary-layer-in-time phenomena appears; that is, near , we have some influence of the initial conditions that cannot be captured by the regular expansion. We can fix that by introducing the appropriate boundary-layer correctors, as proposed in [10] for Newtonian flow. However, it can be proved that those correctors would have the exponential decay towards zero; that is, it would not affect the approximation outside the boundary layers in time. More precisely, it would only serve for the convergence proof, namely, to derive satisfactory error estimates in the appropriate norms. Detailed analysis of the boundary layers together with the rigorous justification of the complete asymptotic expansion (up to a general order ) will be brought in the forthcoming paper by the authors of the present paper.

6. Numerical Example

Finally, in this section, we aim to visually present the derived asymptotic solution for problem (7)–(9). We take the following values for the viscosity constants (see e.g., [21, 22]):

The domain is circular with radius . Furthermore, all the visualizations are done for fixed time and given flow rates , , and . The linear and angular momentum force functions are neglected; namely, . Therefore, it follows that (see (27), (28a)-(28b), (34), (36), (43), and (46)-(47)) and (see (55), (61), and (68)). As a result, the asymptotic solution reads

where , have been explicitly computed in (56)–(59), (63)–(66), and (70)–(72).

In the following, we depict the second-order correctors for the velocity and microrotation (see Figures 2 and 3) as well as the asymptotic approximation for the velocity and microrotation, given by (77), for (see Figures 4 and 5). We omit the visualizations for the zero-order approximations and first-order correctors since they are of the classical Poiseuille form (for the velocity) and of the same form as the second-order corrector for the microrotation.

Figure 2: Second-order velocity corrector .

(a)

(b)

Figure 3: Second-order microrotation corrector (a) and (b).

(a)

(b)

Figure 4: Velocity approximation (a) and zero-order approximation (b).

(a)

(b)

Figure 5: Microrotation approximation (a) and (b).

6.1. Second-Order Correctors

We notice that the velocity second-order corrector is of shape and scale that will affect our asymptotic approximation if is not too small (e.g., ). This has been visually verified in Figure 4. The microrotation second-order corrector does not affect the asymptotic approximation in a significant way as the first-order corrector is of the same shape and similar scale. Note that the second-order correctors are scaled with the small parameter in the asymptotic approximation (77).

6.2. Asymptotic Approximation

In Figure 4 we compare our asymptotic approximation and the zero-order approximation . The second-order corrector for the velocity affects the approximation in the whole domain, with a clear impact near the boundary of the domain, correcting the Poiseuille zero-order approximation, and we can clearly observe the effect visually. The microrotation approximation is the scaled sum of the first- and second-order corrector, thus having the same form as the first- and second-order correctors.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

The authors of this work have been supported by the Croatian Science Foundation (Scientific Project 3955: Mathematical Modeling and Numerical Simulations of Processes in Thin or Porous Domains).

References

E. L. Aero, A. N. Bulygin, and E. V. Kuvshinskii, “Asymmetric hydromechanics,” Journal of Applied Mathematics and Mechanics, vol. 29, no. 2, pp. 333–346, 1965.

A. C. Eringen, “Theory of micropolar fluids,” Journal of Applied Mathematics and Mechanics, vol. 16, pp. 1–18, 1966.

G. Lukaszewicz, Micropolar Fluids: Theory and Application, Birkhäuser, Boston, Mass, USA, 1999.

C.-L. Chang, “Numerical simulation of micropolar fluid flow along a flat plate with wall conduction and buoyancy effects,” Journal of Physics D: Applied Physics, vol. 39, no. 6, pp. 1132–1140, 2006.

A. Ishak, Y. Y. Lok, and I. Pop, “Stagnation point flow over a shrinking sheet in micropolar fluid,” Chemical Engineering Communications, vol. 197, no. 11, pp. 1417–1427, 2010.

I. Abdullah and N. Amin, “A micropolar fluid model of blood flow through a tapered artery with a stenosis,” Mathematical Methods in the Applied Sciences, vol. 33, no. 16, pp. 1910–1923, 2010.

D. Dupuy, G. P. Panasenko, and R. Stavre, “Asymptotic methods for micropolar fluids in a tube structure,” Mathematical Models and Methods in Applied Sciences, vol. 14, no. 5, pp. 735–758, 2004.

I. Pažanin, “On the micropolar flow in a circular pipe: the effects of the viscosity coefficients,” Theoretical & Applied Mechanics Letters, vol. 1, Article ID 062004, 2011.

E. Marušić-Paloka, I. Pažanin, and M. Radulović, “Flow of a micropolar fluid through a channel with small boundary perturbation,” Zeitschrift fur Naturforschung A, vol. 71, no. 7, pp. 607–619, 2016.

G. Panasenko and K. Pileckas, “Asymptotic analysis of the nonsteady viscous flow with a given flow rate in a thin pipe,” Applicable Analysis: An International Journal, vol. 91, no. 3, pp. 559–574, 2012.

G. Panasenko and K. Pileckas, “Asymptotic analysis of the non-steady Navier-Stokes equations in a tube structure. I. The case without boundary-layer-in-time,” Nonlinear Analysis: Theory, Methods & Applications, vol. 122, pp. 125–168, 2015.

G. Panasenko and K. Pileckas, “Asymptotic analysis of the non-steady Navier-Stokes equations in a tube structure. II. General case,” Nonlinear Analysis: Theory, Methods & Applications, vol. 125, pp. 582–607, 2015.

I. Pažanin, “Effective flow of micropolar fluid through a thin or long pipe,” Mathematical Problems in Engineering, vol. 2011, Article ID 127070, 18 pages, 2011.

I. Pažanin, “Asymptotic behavior of micropolar fluid flow through a curved pipe,” Acta Applicandae Mathematicae, vol. 116, no. 1, pp. 1–25, 2011.

I. Pažanin, “Investigation of micropolar fluid flow in a helical pipe via asymptotic analysis,” Communications in Nonlinear Science and Numerical Simulation, vol. 18, no. 3, pp. 528–540, 2013.

M. Beneš and I. Pažanin, “Effective flow of incompressible micropolar fluid through a system of thin pipes,” Acta Applicandae Mathematicae, vol. 143, pp. 29–43, 2016.

A. A. El-Bary, “Exponential solution of a problem of two-dimensional motion of micropolar fluid in a half-plane,” Applied Mathematics and Computation, vol. 165, no. 1, pp. 81–93, 2005.

N. A. Khan and M. Jamil, “Analytic solution for creping flow of an unsteady micropolar fluid,” International Journal of Applied Mathematics and Mechanics, vol. 5, pp. 39–47, 2009.

A. Zaman, N. Ali, and O. Anwar Bég, “Numerical simulation of unsteady micropolar hemodynamics in a tapered catheterized artery with a combination of stenosis and aneurysm,” Medical & Biological Engineering & Computing, vol. 54, no. 9, pp. 1423–1436, 2016.

S. H. M. Saleh, N. M. Arifin, R. Nazar, and I. Pop, “Unsteady micropolar fluid over a permeable curved stretching shrinking surface,” Mathematical Problems in Engineering, vol. 2017, Article ID 3085249, 13 pages, 2017.

I. Papautsky, J. Brazzle, T. Ameel, and A. B. Frazier, “Laminar fluid behavior in microchannels using micropolar fluid theory,” Sensors and Actuators A: Physical, vol. 73, no. 1-2, pp. 101–108, 1999.

K. P. Madasu and M. K. Gurdatta, “Steady rotation of micropolar fluid sphere in concentric spherical container,” Procedia Engineering, vol. 127, pp. 469–475, 2015.