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).
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 readwith 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 deduceTaking into account the problem under consideration, it is plausible to consider the following scaling with respect to small parameter : In view of that, we expandand consider problem (7)–(9) withand 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) asNow, 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 getNote 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 takingwhere 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 readsSimilarly, 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 obtainThe 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 getLet us introduce , as the solutions of the following problems:Taking into account (25) and passing to the polar coordinates easily giveWe seek for the solution of system (30) in the formleading toNow we compute the corrector for the microrotation. Plugging expression (27) for into (29) yieldsSimilarly, 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 yieldTaking into account (27) and (36), we get the following problem for :We rewriteand introduce as the solution of the following problem:Passing to the polar coordinates providesWe seek for the solution of (39) in the formimplyingPlugging (28a)-(28b) and (34) into (37b) and (37c), we get the following system for the second-order microrotation corrector :After rewriting, we obtainTedious but straightforward calculation giveswherewhere , , , , , , , , ,