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

Authors: Igor Pažanin; Marko Radulović
Journal: Mathematical Problems in Engineering
 (2018)

---

## 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.


---

## Body

## 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 Section2, 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 J=2. 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 (1)Ωϵ=x∈R3:x1∈R,x′=x2,x3∈Bϵ=ϵB, where (2)B=x′∈R2:x′<Ris the circular cross-section of the pipe with constant diameter R>0. We denote by (x1,(x2,x3))≡(x1,x′) the Cartesian coordinates, with x1 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:(3)∂uϵ∂t-ν+νrΔuϵ+uϵ·∇uϵ+∇pϵ=2νrrot⁡wϵ+fϵx∈Ωϵ,div⁡uϵ=0x∈Ωϵ,∂wϵ∂t-ca+cdΔwϵ+uϵ·∇wϵ-c0+cd-ca∇div⁡wϵ+4μrwϵ=2νrrot⁡uϵ+gϵx∈Ωϵ.Here uϵ(x1,x′,t)=(u1ϵ(x1,x′,t),u2ϵ(x1,x′,t),u3ϵ(x1,x′,t)) stands for the velocity field, pϵ(x1,x′,t) is the pressure, and wϵ(x′,x3,t)=(w1ϵ(x1,x′,t),w2ϵ(x1,x′,t),w3ϵ(x1,x′,t)) 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:(4)uϵ=0,wϵ=0on∂Ωϵ,vϵx,0=aϵx,wϵx,0=bϵx,0,along with the flux condition (5)∫Bϵu1ϵx1,x′,tdx′=Ft. The positive constants are the Newtonian viscosity ν and the microrotation viscosity νr, while c0, ca, and cd represent the coefficients of the angular viscosities. To simplify the notation, we abbreviate μ=ν+νr,α=ca+cd,β=c0+cd-ca,a=2νr. The external sources of linear and angular momentum are given by the functions fϵ=(f1ϵ,f2ϵ,f3ϵ) and gϵ=(g1ϵ,g2ϵ,g3ϵ), respectively.As indicated in the Introduction, we assume that the solution of the problem (3) has the micropolar Poiseuille form; namely,(6)uϵx,t=vϵx′,t,0,0,wϵx,t=0,w2ϵx′,t,w3ϵx′,t,pϵx,t=-qtx1+p0t,where p0(t) is an arbitrary function of t. We also assume that aϵ=(a1ϵ,0,0),bϵ=(0,b2ϵ,b3ϵ),fϵ=(f1ϵ,0,0),gϵ=(0,g2ϵ,g3ϵ), with a1ϵ=aϵ(x′), b2ϵ=b2ϵ(x′), b3ϵ=b3ϵ(x′), f1ϵ=fϵ(x′,t), g2ϵ(x′,t), and g3ϵ(x′,t) being independent of the longitudinal variable x1∈R.Substituting the micropolar Poiseuille solution (6) into system (3), we obtain the following initial boundary value problem for x′∈Bϵ:(7)∂pϵ∂x2=0,∂pϵ∂x3=0,∂vϵ∂t-μΔx′vϵ-qt=a∂w3ϵ∂x2-∂w2ϵ∂x3+fϵ,∂w2ϵ∂t-αΔx′w2ϵ-β∂∂x2∂w2ϵ∂x2+∂w3ϵ∂x3+2aw2ϵ=a∂vϵ∂x3+g2ϵ,∂w3ϵ∂t-αΔx′w3ϵ-β∂∂x3∂w2ϵ∂x2+∂w3ϵ∂x3+2aw3ϵ=-a∂vϵ∂x2+g3ϵ,where Δx′v=∂2v/∂x22+∂2v/∂x32. The corresponding boundary conditions read(8a)vϵ∂Bϵ=0,(8b)w2ϵ∂Bϵ=0,(8c)w3ϵ∂Bϵ=0,(8d)vϵx′,0=aϵx′,(8e)w1ϵx′,0=b2ϵx′,(8f)w2ϵx′,0=b3ϵx′,with given flux rate:(9)∫Bϵvϵx′,tdx′=Ft.Following the idea from [10–12], we represent the solution (vϵ(x′,t),w^ϵ(x′,t))=(vϵ(x′,t),w2ϵ(x′,t),w3ϵ(x′,t)) as (10)vϵx′,t=V1ϵx′,t+V2ϵx′,t,w^ϵx′,t=Waϵx′,t+Wbϵx′,t.Here (V1ϵ(x′,t),Waϵ(x′,t))=(V1ϵ(x′,t),W2,1ϵ(x′,t),W3,1ϵ(x′,t)) is the solution of the initial boundary value problem for the micropolar heat equation:(11)∂V1ϵ∂t-μΔx′V1ϵ=a∂W3,1ϵ∂x2-∂W2,1ϵ∂x3+fϵ,∂W2,1ϵ∂t-αΔx′W2,1ϵ-β∂∂x2∂W2,1ϵ∂x2+∂W3,1ϵ∂x3+2aW2,1ϵ=a∂V1ϵ∂x3+g2ϵ,∂W3,1ϵ∂t-αΔx′W3,1ϵ-β∂∂x3∂W2,1ϵ∂x2+∂W3,1ϵ∂x3+2aW3,1ϵ=-a∂V1ϵ∂x2+g3ϵ,V1ϵ∂Bϵ=0,W2,1ϵ∂Bϵ=0,W3,1ϵ∂Bϵ=0,V1ϵx′,0=aϵx′,W2,1ϵx′,0=b2ϵx′,W3,1ϵx′,0=b3ϵx′.On the other hand, (V2ϵ(x′,t),Wbϵ(x′,t),q(t))=(V2ϵ(x′,t),W2,2ϵ(x′,t),W3,2ϵ(x′,t),q(t)) denotes the solution of the micropolar inverse problem endowed with the homogeneous initial condition and the given flux rate on the cross-section Bϵ:(12)∂V2ϵ∂t-μΔx′V2ϵ=a∂W3,2ϵ∂x2-∂W2,2ϵ∂x3+qt,∂W2,2ϵ∂t-αΔx′W2,2ϵ-β∂∂x2∂W2,2ϵ∂x2+∂W3,2ϵ∂x3+2aW2,2ϵ=a∂V2ϵ∂x3,∂W3,2ϵ∂t-αΔx′W3,2ϵ-β∂∂x3∂W2,2ϵ∂x2+∂W3,2ϵ∂x3+2aW3,2ϵ=-a∂V2ϵ∂x2,V2ϵ∂Bϵ=0,W2,2ϵ∂Bϵ=0,W3,2ϵ∂Bϵ=0,V2ϵx′,0=0,W2,2ϵx′,0=0,W3,2ϵx′,0=0,∫BϵV2ϵx′,tdx′=Ht,where (13)Ht=Ft-∫BϵV1ϵx′,tdx′. From (9) we deduce(14)H0=0.Taking into account the problem under consideration, it is plausible to consider the following scaling with respect to small parameter ϵ: (15)aϵ,b2ϵ,b3ϵ~1,fϵ~ϵ-2,g2ϵ~ϵ-2,g3ϵ~ϵ-2,F~ϵ2. In view of that, we expand(16a)a~y′=a~0y′+ϵa~1y′+ϵ2a~2y′,(16b)b~2y′=b~20y′+ϵb~21y′+ϵ2b~22y′,(16c)b~3y′=b~30y′+ϵb~31+ϵ2b~32y′,(16d)f~y′,t=1ϵ2f~0y′,t+1ϵf~1y′,t+f~2y′,t,(16e)g~2y′,t=1ϵ2g~20y′,t+1ϵg~21y′,t+g~22y′,t,(16f)g~3y′,t=1ϵ2g~30y′,t+1ϵg~31y′,t+g~32y′,t,(16g)Ft=ϵ2F0t+ϵ3F1t+ϵ4F2t,and consider problem (7)–(9) with(17)aϵx=a~x′ϵ,b2ϵx=b~2x′ϵ,b3ϵx=b~3x′ϵ,fϵx,t=f~x′ϵ,t,g2ϵx,t=g~2x′ϵ,t,g3ϵx,t=g~3x′ϵ,tand F(t) 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 onB instead of Bϵ. We introduce the change of variables y′=x′/ϵ and obtain the following system of equations deduced from (11):(18)∂V~1y′,t∂t-μϵ2Δy′V~1y′,t=aϵ∂W~3,1y′,t∂y2-∂W~2,1y′,t∂y3+f~y′,t,∂W~2,1y′,t∂t-αϵ2Δy′W~2,1y′,t-βϵ2∂∂y2∂W~2,1y′,t∂y2+∂W~3,1y′,t∂y3+2aW~2,1y′,t=aϵ∂V~1y′,t∂y3+g~2y′,t,∂W~3,1y′,t∂t-αϵ2Δy′W~3,1y′,t-βϵ2∂∂y3∂W~2,1y′,t∂y2+∂W~3,1y′,t∂y3+2aW~3,1y′,t=-aϵ∂V~1y′,t∂y2+g~3y′,t,where V~1(y′,t)=V1ϵ(ϵy′,t),W~2,1(y′,t)=W2,1ϵ(ϵy′,t),W~3,1(y′,t)=W3,1ϵ(ϵy′,t). The boundary and initial conditions are the following:(19)V~1y′,t∂B=0,W~2,1y′,t∂B=0,W~3,1y′,t∂B=0,V~1y′,0=a~y′,W~2,1y′,0=b~2y′,W~3,1y′,0=b~3y′.We rewrite problem (18) as(20)-μΔy′V~1y′,t=aϵ∂W~3,1y′,t∂y2-∂W~2,1y′,t∂y3+ϵ2f~y′,t-ϵ2∂V~1y′,t∂t,-αΔy′W~2,1y′,t-β∂∂y2∂W~2,1y′,t∂y2+∂W~3,1y′,t∂y3=aϵ∂V~1y′,t∂y3+ϵ2g~2y′,t-ϵ2∂W~2,1y′,t∂t-2aϵ2W~2,1y′,t,-αΔy′W~3,1y′,t-β∂∂y3∂W~2,1y′,t∂y2+∂W~3,1y′,t∂y3=-aϵ∂V~1y′,t∂y2+ϵ2g~3y′,t-ϵ2∂W~3,1y′,t∂t-2aϵ2W~3,1y′,t.Now, we construct the formal asymptotic expansion of the solution (V~1apr(y′,t),W~2,1apr(y′,t),W~3,1apr(y′,t)) in powers of small parameter ϵ up to ϵ2 in the following way:(21)V~1apry′,t=V~10y′,t+ϵV~11y′,t+ϵ2V~12y′,t,W~2,1apry′,t=W~2,10y′,t+ϵW~2,11y′,t+ϵ2W~2,12y′,t,W~3,1apry′,t=W~3,10y′,t+ϵW~3,11y′,t+ϵ2W~3,12y′,t.Due to the small thickness of the pipe, it is reasonable to assume that the functions f~(y′,t), g~2(y′,t), g~3(y′,t) are independent of the cross-section variables y′=(y2,y3). 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 ϵ0=1, we get(22a)-μΔy′V~10y′,t=f~0t,(22b)-αΔy′W~2,10y′,t-β∂∂y2∂W~2,10y′,t∂y2+∂W~3,10y′,t∂y3=g~20t,(22c)-αΔy′W~3,10y′,t-β∂∂y3∂W~2,10y′,t∂y2+∂W~3,10y′,t∂y3=g~30t,(22d)V~10y′,t∂B=0,(22e)W~2,10y′,t∂B=0,(22f)W~3,10y′,t∂B=0.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(23)V~10y′,t=1μχ0y′f~0t,where χ0(y′) denotes the solution of the auxiliary problem posed on B:(24)-Δy′χ0y′=1,y′∈B,χ0y′∂B=0.Since we assumed the pipe has a circular cross-section, namely,(25)B=y′∈R2:y′<R,R>0,we can pass to polar coordinates (r,θ) and compute χ0 explicitly from (24):(26)χ0y′=14R2-y′2.The explicit expression for the zero-order velocity approximation V~10(y′,t) now reads(27)V~10y′,t=14μR2-y′2f~0t.Similarly, it can be straightforwardly verified that problem (22b), (22c) with the boundary conditions (22e)-(22f) for microrotation will be satisfied for(28a)W~2,10y′,t=22α+βχ0y′g~20t=122α+βR2-y′2g~20t,(28b)W~3,10y′,t=22α+βχ0y′g~30t=122α+βR2-y′2g~30t.


### 3.2. First-Order Corrector

Now, we compute the first-order corrector(V~11(y′,t),W~2,11(y′,t),W~3,11(y′,t)). Inserting (21) and (16d)–(16f) into the system of equations (20), after collecting the terms of order ϵ, we obtain(29)-μΔy′V~11y′,t=a∂W~3,10y′,t∂y2-∂W~2,10y′,t∂y3+f~1t,-αΔy′W~2,11y′,t-β∂∂y2∂W~2,11y′,t∂y2+∂W~3,11y′,t∂y3=a∂V~10y′,t∂y3+g~21t,-αΔy′W~3,11y′,t-β∂∂y3∂W~2,11y′,t∂y2+∂W~3,11y′,t∂y3=-a∂V~10y′,t∂y2+g~31t,V~11y′,t∂B=0,W~2,11y′,t∂B=0,W~3,11y′,t∂B=0.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 (W~2,10(y′,t),W~3,10(y′,t)), we get(30)μΔy′V~11y′,t=ag~30t2α+βy2-g~20t2α+βy3-f~1t,V~11y′,t∂B=0.Let us introduce χi(y′),i=1,2, as the solutions of the following problems:(31)Δy′χiy′=yi+1,y′∈B,χiy′∂B=0.Taking into account (25) and passing to the polar coordinates easily give(32)χ1y′=18y′2-R2y2,χ2y′=18y′2-R2y3.We seek for the solution of system (30) in the form(33)V~11y′,t=1μag~30t2α+βχ1y′-ag~20t2α+βχ2y′+χ0y′f~1tleading to(34)V~11y′,t=18μy′2-R2ag~30t2α+βy2-ag~20t2α+βy3-2f~1t.Now we compute the corrector for the microrotation. Plugging expression (27) for V~10(y′,t) into (29) yields(35)αΔy′W~2,11y′,t+β∂∂y2∂W~2,11y′,t∂y2+∂W~3,11y′,t∂y3=af~0t2μy3-g~21t,αΔy′W~3,11y′,t+β∂∂y3∂W~2,11y′,t∂y2+∂W~3,11y′,t∂y3=-af~0t2μy2-g~31t,W~2,11y′,t∂B=0,W~3,11y′,t∂B=0.Similarly, as for V~11(y′,t), we obtain explicit expression for (W~2,11(y′,t),W~3,11(y′,t)); namely,(36)W~2,11y′,t=y′2-R2y3a16μαf~0t-122α+βg~21t,W~3,11y′,t=y′2-R2-y2a16μαf~0t-122α+βg~31t.


### 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 ϵ2 yield(37a)-μΔy′V~12y′,t=a∂W~3,11y′,t∂y2-∂W~2,11y′,t∂y3+f~2t-∂V~10y′,t∂t,(37b)-αΔy′W~2,12y′,t-β∂∂y2∂W~2,12y′,t∂y2+∂W~3,12y′,t∂y3=a∂V~11y′,t∂y3+g~22t-∂W~2,10y′,t∂t-2aW~2,10y′,t,(37c)-αΔy′W~3,12y′,t-β∂∂y3∂W~2,12y′,t∂y2+∂W~3,12y′,t∂y3=-a∂V~11y′,t∂y2+g~32t-∂W~3,10y′,t∂t-2aW~3,10y′,t,(37d)V~12y′,t∂B=0,(37e)W~2,12y′,t∂B=0,(37f)W~3,12y′,t∂B=0.Taking into account (27) and (36), we get the following problem for V~12(y′,t):(38)-μΔy′V~12y′,t=a-af~0t4μαy22+y32-R22-g~31t2α+βy2-g~21t2α+βy3+f~2t-R2-y22-y3214μddtf~0t,V~12y′,t∂B=0.We rewrite(39)-μΔy′V~12y′,t=-a2f~0t+αd/dtf~0t4μαy22+y32-ag~31t2α+βy2-ag~21t2α+βy3+f~2t+a2f~0t-2αd/dtf~0tR28μα,V~12y′,t∂B=0and introduce χ3(y′) as the solution of the following problem:(40)Δy′χ3y′=y22+y32,y′∈B,χ3y′∂B=0.Passing to the polar coordinates provides(41)χ3y′=116y′4-R4.We seek for the solution of (39) in the form(42)V~12y′,t=-1μ-a2f~0t+αd/dtf~0t4μαχ3y′+1μag~31t2α+βχ1y′+1μag~21t2α+βχ2y′+1μχ0y′f~2t+R2μa2f~0t-2αd/dtf~0t8μαχ0y′,implying(43)V~12y′,t=a2f~0t-αd/dtf~0t64μ2αy′4-R4+ag~31t8μ2α+βy2+ag~21t8μ2α+βy3y′2-R2+f~2t4μ+a2f~0t-2αd/dtf~0tR232μ2αR2-y′2.Plugging (28a)-(28b) and (34) into (37b) and (37c), we get the following system for the second-order microrotation corrector (W~2,12(y′,t),W~3,12(y′,t)):(44)-αΔy′W~2,12y′,t-β∂∂y2∂W~2,12y′,t∂y2+∂W~3,12y′,t∂y3=a-ag~20t8μ2α+βy22+3y32+ag~30t4μ2α+βy2y3+f~1t2μy3+aR2g~20t8μ2α+β+g~22t+122α+βy22+y32-R2ddtg~20t+a2α+βy22+y32-R2g~20t,-αΔy′W~3,12y′,t-β∂∂y3∂W~2,12y′,t∂y2+∂W~3,12y′,t∂y3=a-ag~30t8μ2α+β3y22+y32+ag~20t4μ2α+βy2y3+f~1t2μy2+aR2g~30t8μ2α+β+g~32t+122α+βy22+y32-R2ddtg~30t+a2α+βy22+y32-R2g~30t,W~2,12y′,t∂B=0,W~3,12y′,t∂B=0.After rewriting, we obtain(45)-αΔy′W~2,12y′,t-β∂∂y2∂W~2,12y′,t∂y2+∂W~3,12y′,t∂y3=-a2g~20t+4μd/dtg~20t+8aμg~20t8μ2α+βy22+-3a2g~20t+4μd/dtg~20t+8aμg~20t8μ2α+βy32+a2g~30t4μ2α+βy2y3+af~1t2μy3+a2R2g~20t+8μ2α+βg~22t-4R2μd/dtg~20t-8aμR2g~20t8μ2α+β,-αΔy′W~3,12y′,t-β∂∂y3∂W~2,12y′,t∂y2+∂W~3,12y′,t∂y3=-3a2g~30t+4μd/dtg~30t+8aμg~30t8μ2α+βy22+-a2g~30t+4μd/dtg~30t+8aμg~30t8μ2α+βy32+a2g~20t4μ2α+βy2y3+af~1t2μy2+a2R2g~30t+8μ2α+βg~32t-4R2μd/dtg~30t-8aμR2g~30t8μ2α+β,W~2,12y′,t∂B=0,W~3,12y′,t∂B=0.Tedious but straightforward calculation gives(46)W~2,12y′,t=y′2-R2C1y22+C2y2y3+C3y32+C4y2+C5y3+C6,W~3,12y′,t=y′2-R2C7y22+C8y2y3+C9y32+C10y2+C11y3+C12,where(47)C1=2α-12β2/12α+6βC3+A-3βH/12α+6β-14α-12β+12β2/12α+6β,C2=-4βC7+C9-C12α+6β,C3=B-3βH/12α+6β-14α-12β+12β2/12α+6β--2α-2β+12β2/12α+6βA-3βH/12α+6β-2α-2β+12β2/12α+6β2α-12β2/12α+6β+-14α-2β+12β2/12α+6β2,C4=0,C5=I+8α+2βC10-2β,C6=2αR2C1+C3+2βR2C1+βR2C8-E4α+2β,C7=F-3βC/12α+6β+2α+2β-12β2/12α+6βC9-14α-2β+12β2/12α+6β,C8=-4βC1+C3-H12α+6β,C9=G-3βC/12α+6β-14α-2β+12β2/12α+6β+2α-12β2/12α+6βF-3βC/12α+6β-2α+12β2/12α+6β2α+2β-12β2/12α+6β+-14α-12β+12β2/12α+6β-14α-2β+12β2/12α+6β,C10=2Dβ-I8α+2β8α+2β2-4β2,C11=0,C12=2αR2C7+C9+2βR2C9+βR2C2-J4α+2β,where A=-a2g~20(t)+4μd/dtg~20(t)+8aμg~20(t)/8μ(2α+β), B=-3a2g~20(t)+4μd/dtg~20(t)+8aμg~20(t)/8μ(2α+β), C=a2g~30(t)/4μ(2α+β), D=af~1(t)/2μ, E=a2R2g~20(t)+8μ(2α+β)g~22(t)-4R2μd/dtg~20(t)-8aμR2g~20(t)/8μ(2α+β), F=-3a2g~30(t)+4μd/dtg~30(t)+8aμg~30(t)/8μ(2α+β), G=-a2g~30(t)+4μd/dtg~30(t)+8aμg~30(t)/8μ(2α+β), H=a2g~20(t)/4μ(2α+β), I=af~1(t)/2μ, J=a2R2g~30(t)+8μ(2α+β)g~32(t)-4R2μd/dtg~30(t)-8aμR2g~30(t)/8μ(2α+β).


## 3.1. Zero-Order Approximation

Plugging (21) and (16d)–(16f) into system (20) and collecting the terms of order ϵ0=1, we get(22a)-μΔy′V~10y′,t=f~0t,(22b)-αΔy′W~2,10y′,t-β∂∂y2∂W~2,10y′,t∂y2+∂W~3,10y′,t∂y3=g~20t,(22c)-αΔy′W~3,10y′,t-β∂∂y3∂W~2,10y′,t∂y2+∂W~3,10y′,t∂y3=g~30t,(22d)V~10y′,t∂B=0,(22e)W~2,10y′,t∂B=0,(22f)W~3,10y′,t∂B=0.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(23)V~10y′,t=1μχ0y′f~0t,where χ0(y′) denotes the solution of the auxiliary problem posed on B:(24)-Δy′χ0y′=1,y′∈B,χ0y′∂B=0.Since we assumed the pipe has a circular cross-section, namely,(25)B=y′∈R2:y′<R,R>0,we can pass to polar coordinates (r,θ) and compute χ0 explicitly from (24):(26)χ0y′=14R2-y′2.The explicit expression for the zero-order velocity approximation V~10(y′,t) now reads(27)V~10y′,t=14μR2-y′2f~0t.Similarly, it can be straightforwardly verified that problem (22b), (22c) with the boundary conditions (22e)-(22f) for microrotation will be satisfied for(28a)W~2,10y′,t=22α+βχ0y′g~20t=122α+βR2-y′2g~20t,(28b)W~3,10y′,t=22α+βχ0y′g~30t=122α+βR2-y′2g~30t.


## 3.2. First-Order Corrector

Now, we compute the first-order corrector(V~11(y′,t),W~2,11(y′,t),W~3,11(y′,t)). Inserting (21) and (16d)–(16f) into the system of equations (20), after collecting the terms of order ϵ, we obtain(29)-μΔy′V~11y′,t=a∂W~3,10y′,t∂y2-∂W~2,10y′,t∂y3+f~1t,-αΔy′W~2,11y′,t-β∂∂y2∂W~2,11y′,t∂y2+∂W~3,11y′,t∂y3=a∂V~10y′,t∂y3+g~21t,-αΔy′W~3,11y′,t-β∂∂y3∂W~2,11y′,t∂y2+∂W~3,11y′,t∂y3=-a∂V~10y′,t∂y2+g~31t,V~11y′,t∂B=0,W~2,11y′,t∂B=0,W~3,11y′,t∂B=0.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 (W~2,10(y′,t),W~3,10(y′,t)), we get(30)μΔy′V~11y′,t=ag~30t2α+βy2-g~20t2α+βy3-f~1t,V~11y′,t∂B=0.Let us introduce χi(y′),i=1,2, as the solutions of the following problems:(31)Δy′χiy′=yi+1,y′∈B,χiy′∂B=0.Taking into account (25) and passing to the polar coordinates easily give(32)χ1y′=18y′2-R2y2,χ2y′=18y′2-R2y3.We seek for the solution of system (30) in the form(33)V~11y′,t=1μag~30t2α+βχ1y′-ag~20t2α+βχ2y′+χ0y′f~1tleading to(34)V~11y′,t=18μy′2-R2ag~30t2α+βy2-ag~20t2α+βy3-2f~1t.Now we compute the corrector for the microrotation. Plugging expression (27) for V~10(y′,t) into (29) yields(35)αΔy′W~2,11y′,t+β∂∂y2∂W~2,11y′,t∂y2+∂W~3,11y′,t∂y3=af~0t2μy3-g~21t,αΔy′W~3,11y′,t+β∂∂y3∂W~2,11y′,t∂y2+∂W~3,11y′,t∂y3=-af~0t2μy2-g~31t,W~2,11y′,t∂B=0,W~3,11y′,t∂B=0.Similarly, as for V~11(y′,t), we obtain explicit expression for (W~2,11(y′,t),W~3,11(y′,t)); namely,(36)W~2,11y′,t=y′2-R2y3a16μαf~0t-122α+βg~21t,W~3,11y′,t=y′2-R2-y2a16μαf~0t-122α+βg~31t.


## 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 ϵ2 yield(37a)-μΔy′V~12y′,t=a∂W~3,11y′,t∂y2-∂W~2,11y′,t∂y3+f~2t-∂V~10y′,t∂t,(37b)-αΔy′W~2,12y′,t-β∂∂y2∂W~2,12y′,t∂y2+∂W~3,12y′,t∂y3=a∂V~11y′,t∂y3+g~22t-∂W~2,10y′,t∂t-2aW~2,10y′,t,(37c)-αΔy′W~3,12y′,t-β∂∂y3∂W~2,12y′,t∂y2+∂W~3,12y′,t∂y3=-a∂V~11y′,t∂y2+g~32t-∂W~3,10y′,t∂t-2aW~3,10y′,t,(37d)V~12y′,t∂B=0,(37e)W~2,12y′,t∂B=0,(37f)W~3,12y′,t∂B=0.Taking into account (27) and (36), we get the following problem for V~12(y′,t):(38)-μΔy′V~12y′,t=a-af~0t4μαy22+y32-R22-g~31t2α+βy2-g~21t2α+βy3+f~2t-R2-y22-y3214μddtf~0t,V~12y′,t∂B=0.We rewrite(39)-μΔy′V~12y′,t=-a2f~0t+αd/dtf~0t4μαy22+y32-ag~31t2α+βy2-ag~21t2α+βy3+f~2t+a2f~0t-2αd/dtf~0tR28μα,V~12y′,t∂B=0and introduce χ3(y′) as the solution of the following problem:(40)Δy′χ3y′=y22+y32,y′∈B,χ3y′∂B=0.Passing to the polar coordinates provides(41)χ3y′=116y′4-R4.We seek for the solution of (39) in the form(42)V~12y′,t=-1μ-a2f~0t+αd/dtf~0t4μαχ3y′+1μag~31t2α+βχ1y′+1μag~21t2α+βχ2y′+1μχ0y′f~2t+R2μa2f~0t-2αd/dtf~0t8μαχ0y′,implying(43)V~12y′,t=a2f~0t-αd/dtf~0t64μ2αy′4-R4+ag~31t8μ2α+βy2+ag~21t8μ2α+βy3y′2-R2+f~2t4μ+a2f~0t-2αd/dtf~0tR232μ2αR2-y′2.Plugging (28a)-(28b) and (34) into (37b) and (37c), we get the following system for the second-order microrotation corrector (W~2,12(y′,t),W~3,12(y′,t)):(44)-αΔy′W~2,12y′,t-β∂∂y2∂W~2,12y′,t∂y2+∂W~3,12y′,t∂y3=a-ag~20t8μ2α+βy22+3y32+ag~30t4μ2α+βy2y3+f~1t2μy3+aR2g~20t8μ2α+β+g~22t+122α+βy22+y32-R2ddtg~20t+a2α+βy22+y32-R2g~20t,-αΔy′W~3,12y′,t-β∂∂y3∂W~2,12y′,t∂y2+∂W~3,12y′,t∂y3=a-ag~30t8μ2α+β3y22+y32+ag~20t4μ2α+βy2y3+f~1t2μy2+aR2g~30t8μ2α+β+g~32t+122α+βy22+y32-R2ddtg~30t+a2α+βy22+y32-R2g~30t,W~2,12y′,t∂B=0,W~3,12y′,t∂B=0.After rewriting, we obtain(45)-αΔy′W~2,12y′,t-β∂∂y2∂W~2,12y′,t∂y2+∂W~3,12y′,t∂y3=-a2g~20t+4μd/dtg~20t+8aμg~20t8μ2α+βy22+-3a2g~20t+4μd/dtg~20t+8aμg~20t8μ2α+βy32+a2g~30t4μ2α+βy2y3+af~1t2μy3+a2R2g~20t+8μ2α+βg~22t-4R2μd/dtg~20t-8aμR2g~20t8μ2α+β,-αΔy′W~3,12y′,t-β∂∂y3∂W~2,12y′,t∂y2+∂W~3,12y′,t∂y3=-3a2g~30t+4μd/dtg~30t+8aμg~30t8μ2α+βy22+-a2g~30t+4μd/dtg~30t+8aμg~30t8μ2α+βy32+a2g~20t4μ2α+βy2y3+af~1t2μy2+a2R2g~30t+8μ2α+βg~32t-4R2μd/dtg~30t-8aμR2g~30t8μ2α+β,W~2,12y′,t∂B=0,W~3,12y′,t∂B=0.Tedious but straightforward calculation gives(46)W~2,12y′,t=y′2-R2C1y22+C2y2y3+C3y32+C4y2+C5y3+C6,W~3,12y′,t=y′2-R2C7y22+C8y2y3+C9y32+C10y2+C11y3+C12,where(47)C1=2α-12β2/12α+6βC3+A-3βH/12α+6β-14α-12β+12β2/12α+6β,C2=-4βC7+C9-C12α+6β,C3=B-3βH/12α+6β-14α-12β+12β2/12α+6β--2α-2β+12β2/12α+6βA-3βH/12α+6β-2α-2β+12β2/12α+6β2α-12β2/12α+6β+-14α-2β+12β2/12α+6β2,C4=0,C5=I+8α+2βC10-2β,C6=2αR2C1+C3+2βR2C1+βR2C8-E4α+2β,C7=F-3βC/12α+6β+2α+2β-12β2/12α+6βC9-14α-2β+12β2/12α+6β,C8=-4βC1+C3-H12α+6β,C9=G-3βC/12α+6β-14α-2β+12β2/12α+6β+2α-12β2/12α+6βF-3βC/12α+6β-2α+12β2/12α+6β2α+2β-12β2/12α+6β+-14α-12β+12β2/12α+6β-14α-2β+12β2/12α+6β,C10=2Dβ-I8α+2β8α+2β2-4β2,C11=0,C12=2αR2C7+C9+2βR2C9+βR2C2-J4α+2β,where A=-a2g~20(t)+4μd/dtg~20(t)+8aμg~20(t)/8μ(2α+β), B=-3a2g~20(t)+4μd/dtg~20(t)+8aμg~20(t)/8μ(2α+β), C=a2g~30(t)/4μ(2α+β), D=af~1(t)/2μ, E=a2R2g~20(t)+8μ(2α+β)g~22(t)-4R2μd/dtg~20(t)-8aμR2g~20(t)/8μ(2α+β), F=-3a2g~30(t)+4μd/dtg~30(t)+8aμg~30(t)/8μ(2α+β), G=-a2g~30(t)+4μd/dtg~30(t)+8aμg~30(t)/8μ(2α+β), H=a2g~20(t)/4μ(2α+β), I=af~1(t)/2μ, J=a2R2g~30(t)+8μ(2α+β)g~32(t)-4R2μd/dtg~30(t)-8aμR2g~30(t)/8μ(2α+β).


## 4. Micropolar Inverse Problem

Now, we turn our attention to micropolar inverse problem (12). First, we notice that(48)Ht=∑j=02ϵj+2Fjt-∑j=02∫BϵϵjV~1jx′ϵ,tdx′=∑j=02ϵj+2Fjt-∑j=02∫Bϵj+2V~1jy′,tdy′=∑j=02ϵj+2Hjt.Denoting(49)Hjt=Fjt-∫BV~1jy′,tdy′,j=0,1,2,we have the compatibility condition (due to (14))(50)Hj0=0,j=0,1,2,or(51)Fj0=∫Ba~jy′dy′,j=0,1,2.Performing the change of variables y′=x′/ϵ and putting q(t)=1/ϵ2s(t), from (12) we deduce(52)-μΔy′V~2y′,t=aϵ∂W~3,2y′,t∂y2-∂W~2,2y′,t∂y3+st-ϵ2∂V~2∂ty′,t,-αΔy′W~2,2y′,t-β∂∂y2∂W~2,2y′,t∂y2+∂W~3,2y′,t∂y3=aϵ∂V~2y′,t∂y3-ϵ2∂W~2,2y′,t∂t-2aϵ2W~2,2y′,t,-αΔy′W~3,2y′,t-β∂∂y3∂W~2,2y′,t∂y2+∂W~3,2y′,t∂y3=-aϵ∂V~2y′,t∂y2-ϵ2∂W~3,2y′,t∂t-2aϵ2W~3,2y′,t,V~2y′,t∂B=0,W~2,2y′,t∂B=0,W~3,2y′,t∂B=0,V~2y′,0=0,W~2,2y′,0=0,W~3,2y′,0=0,∫BV~2y′,tdy′=1ϵ2Ht,where V~2(y′,t)=V2ϵ(ϵy′,t),W~2,2(y′,t)=W2,2ϵ(ϵy′,t),W~3,2(y′,t)=W3,2ϵ(ϵy′,t). In the sequel, we look for the asymptotic approximation in the form(53)V~2apry′,t=V~20y′,t+ϵV~21y′,t+ϵ2V~22y′,t,W~2,2apry′,t=W~2,20y′,t+ϵW~2,21y′,t+ϵ2W~2,22y′,t,W~3,2apry′,t=W~3,20y′,t+ϵW~3,21y′,t+ϵ2W~3,22y′,t,saprt=s0t+ϵs1t+ϵ2s2t.


### 4.1. Zero-Order Approximation

For the zero-order approximation(V~20(y′,t),W~2,20(y′,t),W~3,20(y′,t),s0(t)), we obtain the following problem:(54a)-μΔy′V~20y′,t=s0t,(54b)-αΔy′W~2,20y′,t-β∂∂y2∂W~2,20y′,t∂y2+∂W~3,20y′,t∂y3=0,(54c)-αΔy′W~3,20y′,t-β∂∂y3∂W~2,20y′,t∂y2+∂W~3,20y′,t∂y3=0,(54d)V~20y′,t∂B=0,(54e)W~2,20y′,t∂B=0,(54f)W~3,20y′,t∂B=0,(54g)∫BV~20y′,tdy′=H0t,where H0(t) is given by(55)H0t=F0t-∫BV~10y′,tdy′=F0t-R4π8μf~0t.Equation (54a) with the boundary condition (54d) and the flux rate (54g) can be solved by taking(56)V~20y′,t=14μR2-y′2s0t.Here we choose s0(t) to satisfy the flux condition:(57)s0t=κ0-1H0t,where(58)κ0=∫B14μR2-y′2=R4π8μ>0.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 (g~20(t),g~30(t))≡(0,0) (see Section 3.1). Thus, we conclude(59)W~2,20y′,t,W~3,20y′,t≡0,0.


### 4.2. First-Order Corrector

The first-order corrector(V~21(y′,t),W~2,21(y′,t),W~3,21(y′,t),s1(t)) is given by the following problem:(60a)-μΔy′V~21y′,t=a∂W~3,20y′,t∂y2-∂W~2,20y′,t∂y3+s1t,(60b)-αΔy′W~2,21y′,t-β∂∂y2∂W~2,21y′,t∂y2+∂W~3,21y′,t∂y3=a∂V~20y′,t∂y3,(60c)-αΔy′W~3,21y′,t-β∂∂y3∂W~2,21y′,t∂y2+∂W~3,21y′,t∂y3=-a∂V~20y′,t∂y2,(60d)V~21y′,t∂B=0,(60e)W~2,21y′,t∂B=0,(60f)W~3,21y′,t∂B=0,(60g)∫BV~21y′,tdy′=H1t,where H1(t) is given by(61)H1t=F1t-∫BV~11y′,tdy′=F1t-R4π8μf~1t.We notice that the system is not decoupled anymore. Inserting the expression for (V~20(y′,t),W~2,20(y′,t),W~3,20(y′,t)) given by (56) and (59), we obtain(62a)-μΔy′V~21y′,t=s1t,(62b)-αΔy′W~2,21y′,t-β∂∂y2∂W~2,21y′,t∂y2+∂W~3,21y′,t∂y3=-a2μy3s0t,(62c)-αΔy′W~3,21y′,t-β∂∂y3∂W~2,21y′,t∂y2+∂W~3,21y′,t∂y3=a2μy2s0t,(62d)V~21y′,t∂B=0,(62e)W~2,21y′,t∂B=0,(62f)W~3,21y′,t∂B=0,(62g)∫BV~21y′,tdy′=H1t.As above, (60a) with the boundary condition (60d) and flux rate (60g) can be solved by taking(63)V~21y′,t=14μR2-y′2s1t,with s1(t) obeying the flux condition:(64)s1t=κ1-1H1t,with(65)κ1=∫B14μR2-y′2=R4π8μ>0.One can easily verify that the system (62b), (62c) for the microrotation corrector with the boundary conditions (62e)-(62f) is satisfied for(66)W~2,21y′,t=y′2-R2aH0t16μακ0y3,W~3,21y′,t=y′2-R2-aH0t16μακ0y2.


### 4.3. Second-Order Corrector

To conclude the analysis for the micropolar inverse problem (52), we identify the terms of order ϵ2 to get the system for the second-order corrector (V~22(y′,t),W~2,22(y′,t),W~3,22(y′,t),s2(t)). It reads(67)-μΔy′V~22y′,t=a∂W~3,21y′,t∂y2-∂W~2,21y′,t∂y3+s2t-∂V~20y′,t∂t,-αΔy′W~2,22y′,t-β∂∂y2∂W~2,22y′,t∂y2+∂W~3,22y′,t∂y3=a∂V~21y′,t∂y3-∂W~2,20y′,t∂t-2aW~2,20y′,t,-αΔy′W~3,22y′,t-β∂∂y3∂W~2,22y′,t∂y2+∂W~3,22y′,t∂y3=-a∂V~21y′,t∂y2-∂W~3,20y′,t∂t-2aW~3,20y′,t,V~22y′,t∂B=0,W~2,22y′,t∂B=0,W~3,22y′,t∂B=0,∫∂BV~22y′,tdy′=H2t,where H2(t) is given by(68)H2t=F2t-∫BV~12y′,tdy′=F2t-πR6αd/dtf~0t-πR6a2f~0t96μ2α-8πμαR4f~2t+πR6a2f~0t-2αd/dtf~0t264μ3α.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: (69)-μΔy′V~22y′,t=-a2H0t4μακ0+d/dtH0t4μκ0y22+y32+s2t+R2a2H0t-2αR2H0t8μακ0,-αΔy′W~2,22y′,t-β∂∂y2∂W~2,22y′,t∂y2+∂W~3,22y′,t∂y3=-aH1t2μκ1y3,-αΔy′W~3,22y′,t-β∂∂y3∂W~2,22y′,t∂y2+∂W~3,22y′,t∂y3=aH1t2μκ1y2,V~22y′,t∂B=0,W~2,22y′,t∂B=0,W~3,22y′,t∂B=0,∫BV~22y′,tdy′=H2t.We can now explicitly compute the velocity corrector V~22(y′,t) as(70)V~22y′,t=-116μ-a2H0t4μακ0+d/dtH0t4μκ0y′4-R4-14μs2t+R2a2H0t-2αR2H0t8μακ0y′2-R2,where we choose s2(t) to satisfy the flux condition:(71)s2t=H2t+1/16μ-a2H0t/4μακ0+d/dtH0t/4μκ0-2πR6/3+1/4μR2a2H0t-2αR2H0t/8μακ0-R4π/2-1/4μ-R4π/2.Similarly, we can explicitly compute the corrector for the microrotation:(72)W~2,22y′,t=y′2-R2aH1t16μακ1y3,W~3,22y′,t=y′2-R2-aH1t16μακ1y2.


## 4.1. Zero-Order Approximation

For the zero-order approximation(V~20(y′,t),W~2,20(y′,t),W~3,20(y′,t),s0(t)), we obtain the following problem:(54a)-μΔy′V~20y′,t=s0t,(54b)-αΔy′W~2,20y′,t-β∂∂y2∂W~2,20y′,t∂y2+∂W~3,20y′,t∂y3=0,(54c)-αΔy′W~3,20y′,t-β∂∂y3∂W~2,20y′,t∂y2+∂W~3,20y′,t∂y3=0,(54d)V~20y′,t∂B=0,(54e)W~2,20y′,t∂B=0,(54f)W~3,20y′,t∂B=0,(54g)∫BV~20y′,tdy′=H0t,where H0(t) is given by(55)H0t=F0t-∫BV~10y′,tdy′=F0t-R4π8μf~0t.Equation (54a) with the boundary condition (54d) and the flux rate (54g) can be solved by taking(56)V~20y′,t=14μR2-y′2s0t.Here we choose s0(t) to satisfy the flux condition:(57)s0t=κ0-1H0t,where(58)κ0=∫B14μR2-y′2=R4π8μ>0.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 (g~20(t),g~30(t))≡(0,0) (see Section 3.1). Thus, we conclude(59)W~2,20y′,t,W~3,20y′,t≡0,0.


## 4.2. First-Order Corrector

The first-order corrector(V~21(y′,t),W~2,21(y′,t),W~3,21(y′,t),s1(t)) is given by the following problem:(60a)-μΔy′V~21y′,t=a∂W~3,20y′,t∂y2-∂W~2,20y′,t∂y3+s1t,(60b)-αΔy′W~2,21y′,t-β∂∂y2∂W~2,21y′,t∂y2+∂W~3,21y′,t∂y3=a∂V~20y′,t∂y3,(60c)-αΔy′W~3,21y′,t-β∂∂y3∂W~2,21y′,t∂y2+∂W~3,21y′,t∂y3=-a∂V~20y′,t∂y2,(60d)V~21y′,t∂B=0,(60e)W~2,21y′,t∂B=0,(60f)W~3,21y′,t∂B=0,(60g)∫BV~21y′,tdy′=H1t,where H1(t) is given by(61)H1t=F1t-∫BV~11y′,tdy′=F1t-R4π8μf~1t.We notice that the system is not decoupled anymore. Inserting the expression for (V~20(y′,t),W~2,20(y′,t),W~3,20(y′,t)) given by (56) and (59), we obtain(62a)-μΔy′V~21y′,t=s1t,(62b)-αΔy′W~2,21y′,t-β∂∂y2∂W~2,21y′,t∂y2+∂W~3,21y′,t∂y3=-a2μy3s0t,(62c)-αΔy′W~3,21y′,t-β∂∂y3∂W~2,21y′,t∂y2+∂W~3,21y′,t∂y3=a2μy2s0t,(62d)V~21y′,t∂B=0,(62e)W~2,21y′,t∂B=0,(62f)W~3,21y′,t∂B=0,(62g)∫BV~21y′,tdy′=H1t.As above, (60a) with the boundary condition (60d) and flux rate (60g) can be solved by taking(63)V~21y′,t=14μR2-y′2s1t,with s1(t) obeying the flux condition:(64)s1t=κ1-1H1t,with(65)κ1=∫B14μR2-y′2=R4π8μ>0.One can easily verify that the system (62b), (62c) for the microrotation corrector with the boundary conditions (62e)-(62f) is satisfied for(66)W~2,21y′,t=y′2-R2aH0t16μακ0y3,W~3,21y′,t=y′2-R2-aH0t16μακ0y2.


## 4.3. Second-Order Corrector

To conclude the analysis for the micropolar inverse problem (52), we identify the terms of order ϵ2 to get the system for the second-order corrector (V~22(y′,t),W~2,22(y′,t),W~3,22(y′,t),s2(t)). It reads(67)-μΔy′V~22y′,t=a∂W~3,21y′,t∂y2-∂W~2,21y′,t∂y3+s2t-∂V~20y′,t∂t,-αΔy′W~2,22y′,t-β∂∂y2∂W~2,22y′,t∂y2+∂W~3,22y′,t∂y3=a∂V~21y′,t∂y3-∂W~2,20y′,t∂t-2aW~2,20y′,t,-αΔy′W~3,22y′,t-β∂∂y3∂W~2,22y′,t∂y2+∂W~3,22y′,t∂y3=-a∂V~21y′,t∂y2-∂W~3,20y′,t∂t-2aW~3,20y′,t,V~22y′,t∂B=0,W~2,22y′,t∂B=0,W~3,22y′,t∂B=0,∫∂BV~22y′,tdy′=H2t,where H2(t) is given by(68)H2t=F2t-∫BV~12y′,tdy′=F2t-πR6αd/dtf~0t-πR6a2f~0t96μ2α-8πμαR4f~2t+πR6a2f~0t-2αd/dtf~0t264μ3α.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: (69)-μΔy′V~22y′,t=-a2H0t4μακ0+d/dtH0t4μκ0y22+y32+s2t+R2a2H0t-2αR2H0t8μακ0,-αΔy′W~2,22y′,t-β∂∂y2∂W~2,22y′,t∂y2+∂W~3,22y′,t∂y3=-aH1t2μκ1y3,-αΔy′W~3,22y′,t-β∂∂y3∂W~2,22y′,t∂y2+∂W~3,22y′,t∂y3=aH1t2μκ1y2,V~22y′,t∂B=0,W~2,22y′,t∂B=0,W~3,22y′,t∂B=0,∫BV~22y′,tdy′=H2t.We can now explicitly compute the velocity corrector V~22(y′,t) as(70)V~22y′,t=-116μ-a2H0t4μακ0+d/dtH0t4μκ0y′4-R4-14μs2t+R2a2H0t-2αR2H0t8μακ0y′2-R2,where we choose s2(t) to satisfy the flux condition:(71)s2t=H2t+1/16μ-a2H0t/4μακ0+d/dtH0t/4μκ0-2πR6/3+1/4μR2a2H0t-2αR2H0t/8μακ0-R4π/2-1/4μ-R4π/2.Similarly, we can explicitly compute the corrector for the microrotation:(72)W~2,22y′,t=y′2-R2aH1t16μακ1y3,W~3,22y′,t=y′2-R2-aH1t16μακ1y2.


## 5. Asymptotic Approximation

The asymptotic approximation related to the micropolar heat problem (11) has the following form:(73)V~1aprx′,t=V~10x′ϵ,t+ϵV~11x′ϵ,t+ϵ2V~12x′ϵ,t,W~2,1aprx′,t=W~2,10x′ϵ,t+ϵW~2,11x′ϵ,t+ϵ2W~2,12x′ϵ,t,W~3,1aprx′,t=W~3,10x′ϵ,t+ϵW~3,11x′ϵ,t+ϵ2W~3,12x′ϵ,t.Note we have calculated all the terms explicitly; that is, V~1jx′/ϵ,t,W~2,1jx′/ϵ,t,W~3,1jx′/ϵ,t,j=0,1,2, are given by (27) and (28a)-(28b), (34), (36), (43), and (46)-(47).For the micropolar inverse problem (12), we obtain the following:(74)V2aprx′,t=V~20x′ϵ,t+ϵV~21x′ϵ,t+ϵ2V~22x′ϵ,t,W2,2aprx′,t=W~2,20x′ϵ,t+ϵW~2,21x′ϵ,t+ϵ2W~2,22x′ϵ,t,W3,2aprx′,t=W~3,20x′ϵ,t+ϵW~3,21x′ϵ,t+ϵ2W~3,22x′ϵ,t,saprt=s0t+ϵs1t+ϵ2s2t,where (V~2j(x′/ϵ,t),W~2,2j(x′/ϵ,t),W~3,2j(x′/ϵ,t),sj(t)),j=0,1,2, are given by the explicit formulae (56)–(59), (63)–(66), and (70)–(72).Collecting the above results, we arrive at the asymptotic approximation(vapr(x′,t),w2apr(x′,t),w3apr(x′,t),qapr(t)) for problem (7)–(9):(75)vaprx′,t=V1aprx′,t+V2aprx′,t,w2aprx′,t=W~2,1aprx′,t+W~2,2aprx′,t,w3aprx′,t=W~3,1aprx′,t+W~3,2aprx′,t,qaprt=1ϵ2saprt. 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 t=0, 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 J) 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]):(76)ν=2.9×10-3,νr=2.32×10-4,a=2νr=4.64×10-4,μ=ν+νr,α=ca+cd=10-6,β=c0+cd-ca=10-6.The domain is circular with radius R=10. Furthermore, all the visualizations are done for fixed time t=10 and given flow rates F0(t)=t+50, F1(t)=t2, and F2(t)=t3. The linear and angular momentum force functions are neglected; namely, f~1(t)≡g~2(t)≡g~3(t)≡0. Therefore, it follows that (V1ϵ,2(x′,t),W2,1ϵ,2(x′,t),W3,1ϵ,2(x′,t))≡0 (see (27), (28a)-(28b), (34), (36), (43), and (46)-(47)) and Fj(t)≡Hj(t),j=0,1,2 (see (55), (61), and (68)). As a result, the asymptotic solution reads(77)vaprx′,t=V~20x′ϵ,t+ϵV~21x′ϵ,t+ϵ2V~22x′ϵ,t,w2aprx′,t=W~2,20x′ϵ,t+ϵW~2,21x′ϵ,t+ϵ2W~2,22x′ϵ,t,w3aprx′,t=W~3,20x′ϵ,t+ϵW~3,21x′ϵ,t+ϵ2W~3,22x′ϵ,t,where (V~2j(x′/ϵ,t),W~2,2j(x′/ϵ,t),W~3,2j(x′/ϵ,t)),j=0,1,2, 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 Figures2 and 3) as well as the asymptotic approximation for the velocity and microrotation, given by (77), for ϵ=0.1 (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 correctorV~22.Figure 3
Second-order microrotation correctorW~2,22 (a) and W~3,22 (b).
(a) (b)Figure 4
Velocity approximationvapr (a) and zero-order approximation V~20 (b).
(a) (b)Figure 5
Microrotation approximationw2apr (a) and w3apr (b).
(a) (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., ϵ=0.1). 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 ϵ2 in the asymptotic approximation (77).


### 6.2. Asymptotic Approximation

In Figure4 we compare our asymptotic approximation vapr and the zero-order approximation V~20. 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.


## 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., ϵ=0.1). 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 ϵ2 in the asymptotic approximation (77).


## 6.2. Asymptotic Approximation

In Figure4 we compare our asymptotic approximation vapr and the zero-order approximation V~20. 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.



---

*Source: 6759876-2018-04-11.xml*