#### Abstract

Boundary layer flow of the Newtonian fluid that is caused by the vibration of inner sphere while the outer sphere is at rest is calculated. Vishik-Lyusternik (Nayfeh refers to this method as the method of composite expansions) method is employed to construct an asymptotic expansion of the solution of the Navier-Stokes equations in the limit of high-frequency vibrations for Reynolds number of . The effect of the Stokes drift of fluid particles is also considered.

#### 1. Introduction

In the present paper we consider boundary layer flow due to the vibrating sphere in a viscous fluid. Boundary layer flow for Newtonian and non-Newtonian fluid’s is studied by many researchers [1, 2]. The interesting effect of viscosity is the flow produced by high-frequency oscillations of a boundary where the flow oscillations do not average to zero, such steady flow is referred to as “steady streaming” or “acoustic streaming” or “nonlinear streaming” [3, 4].

Mathematical modeling of the problem involves two parameters: the inverse of Reynolds number and the inverse Strouhal number , given by where is the amplitude of the velocity of the oscillating sphere, is the radius of the sphere, is the angular frequency of the oscillations, and is the kinematic viscosity of the fluid. Parameter measures the ratio of the amplitude of the displacement of the oscillating sphere to its radius and is assumed to be small, that is, . An asymptotic solution of Navier-Stokes equations in the limit is obtained using the Vishik-Lyusternik method.

Steady secondary flow was investigated numerically in [5] for . The study was carried out to investigate the steady streaming near the rotating cylinder. In [6] steady streaming between two cylinders was studied experimentally for at different amplitude of oscillations. The study was carried out to investigate the behavior of outer boundary on the inner circulations. The unsteady flow around rotating and oscillating cylinder in a viscous incompressible micropolar fluid is studied numerically by [7] for and obtained a good agreement with the experiments. Stream lines pattern was observed between two oscillating walls by Thomas et al. [8]. The study is carried out to show that less force is required by the flat plate to move the fluid rather than wavy wall. Steady flow between two cylinders was studied by [9] where the inner cylinder performs transverse oscillation and outer cylinder performs rotatory oscillations for .

Wang [10] discussed the steady streaming due to oscillations of a sphere for using the method of inner and outer expansion. However his solution is incomplete, which had been understood and explained by Riley [11]. Riley discussed the steady streaming due to a sphere fixed in an oscillating fluid for using method of matched asymptotic expansions. He discussed two cases: (i) for and (ii) for , where in the second case choice of is arbitrary such that . Dohara [12] studied steady streaming for and obtained good results with the experimental results. Recently the experimental work of Kotas et al. [13–15] produced the visualization of steady streaming due to the oscillating spheroids for moderate Reynolds number and small amplitude of oscillations. In recent paper [13], numerical results for steady streaming are presented for and compared with the experimental results of Kotas et al.

In the present paper, steady streaming between two spheres is studied. The results are obtained for . Section 2 describes the mathematical formulation of the problem. Section 3 briefly explains the asymptotic procedure. To construct an asymptotic expansion of the solution of the Navier-Stokes equations, we apply the Vishik-Lyusternik method [16, 17]. This method has an essential advantage that it does not require the procedure of matching inner and outer expansion as in the case of matched asymptotic expansions. The Vishik-Lyusternik method had been used to study viscous boundary layer at a fixed impermeable boundary by Chudov [17]. Recently it has been applied to viscous boundary layers in high Reynolds number flows through a fixed domain with an inlet and an outlet [18] and to viscous flows in a half-plane produced by tangential vibrations on its boundary [19]. Vishik-Lyusternik method has successfully applied to the steady streaming due to a vibrating cylinder and by vibrating wavy wall [20–22]. In Section 4 asymptotic equations are solved. In Section 5 the steady Eulerian velocity is corrected through Stokes drift. This produces the steady Lagrangian velocity which is important because: (i) it is the Lagrangian velocity that is observed in experiments; (ii) it is the Lagrangian velocity that is invariant under the change of the frame of reference from the one fixed in the oscillating sphere to the one fixed in space. Section 6 contains the discussion of the results.

#### 2. Formulation of the Problem

We consider a three-dimensional flow of a viscous incompressible fluid between two spheres with radii and () produced by small translational vibrations of an inner sphere about the axis of the outer sphere which is fixed in space. Let be Cartesian coordinates in space, and let be the position of the center of the inner sphere at time . We assume that is oscillating in with angular frequency and period and has zero mean value, that is, The motion of the fluid is governed by the two-dimensional Navier-Stokes equations: The velocity of the fluid satisfies the standard no-slip condition on the surfaces of the spheres which is: Let , , , and be the characteristic scales for time, length, velocity, and pressure, respectively. Using the dimensionless variables we rewrite (3) in the form where are the inverse Strouhal number and the dimensionless viscosity (the inverse Reynolds number). Equations (4) become Here dots denote differentiation with respect to and is a given function which describes the motion of the inner sphere.

We are interested in the asymptotic behavior of periodic solution of (6) in the high-frequency limit .

It is convenient to introduce parameter such that . Then (6) becomes Boundary condition (8) takes the form We assume that the flow is axisymmetric (in spherical coordinates , axisymmetric flow means , where is the azimuthal angle and , ) in which the flow quantities do not depend upon the azimuthal angle . Thus the time-dependent boundary of the sphere can be described in the parametric form by the equations where is the parameter on the sphere boundary. Now the boundary condition on the sphere can be written as follows: Using the assumption that is small, we expand and in Taylor’s series at . This yields Note that each term on the left side of (13) is evaluated at the averaged position of the sphere.

In spherical polar coordinates , with origin at the axis of the outer sphere, the axisymmetric form of (9) takes the form where and are the velocity components along and directions and subscripts “,” “,” and “” denote partial derivatives, and The boundary condition on the outer sphere becomes where . The boundary condition (13) on the inner sphere takes the form Here Equation (14) is to be solved according to the boundary conditions (16)–(18).

#### 3. Asymptotic Expansion

We seek a solution of (14)–(18) in the form Here and are the boundary layer variables. Functions , , and represent a regular expansion of the solution in power series in (an outer solution), and , , and and , , and correspond to boundary layer corrections to this regular expansion. Superscripts “” and “” correspond to the boundary layers at the inner and outer spheres, respectively. We assume that the boundary layer part of the expansion rapidly decays outside the thin boundary layers, that is,

##### 3.1. Regular Part of the Expansion

Substituting (20)–(22) in (14), we get the following equations: where and for . Explicit expression of and for is given in Appendix A.

##### 3.2. Boundary Layer Expansion

*Boundary Layer at the Inner Sphere.* Let us assume that
Then we substitute (20)–(22) and (30) into (14) and take into account that , , and () satisfy (24)–(29). Then we make the change of variables , expand every function of in Taylor’s series at , and collect terms of equal powers in . This produces the following sequence of equations:
for . Functions , and are defined in terms of , , and . Explicit expression of these functions for is given in Appendix A.

*Boundary Layer at the Outer Sphere.* Let us assume that
Then we substitute (20)–(22) and (37) into (14) and take into account that , , and () satisfy (24)–(29). Then we make the change of variables , expand every function of in Taylor’s series at and collect terms of equal powers in . This produces the following sequence of equations:
for . Functions , and are defined in terms of , , and . Explicit expression of these functions for is given in Appendix A. We require that in all orders the boundary layer corrections to the outer solution rapidly decay outside boundary layers, that is,
for .

Before we try to solve these equations, it is convenient to discuss boundary conditions at and .

##### 3.3. Boundary Conditions

Substituting (20)–(22), (30), and (37) in (16)–(18) and collecting terms of equal powers in , we obtain the following boundary conditions: for . Functions and are defined in terms of , , and . Explicit expression of these functions for is given in Appendix A.

#### 4. Analysis of the Asymptotic Equations

##### 4.1. Leading Order Equations

Taking the derivative of (26) (for ) with respect to and substituting (24)-(25) we get The boundary conditions for the above equation can be obtained from (46) and (48) and with the relation . The boundary conditions are The solution of (54) together with (55) and (26) gives and , the (leading order) oscillating outer flow which are

*Inner Sphere.* Consider now the leading order equations (31)–(33). The condition of decay at infinity (in the boundary layer variable ) for and (32) has a consequence that . Hence, we have the standard heat equation
Boundary condition for at follows from (47):
The solution of (57) must also satisfy the decay condition (44).

Substitution of (56) in the boundary condition (58) yields Equation (57) subject to the boundary conditions (59) and (44) can be solved by standard methods. The solution is given by where On averaging (57), we obtain

The only solution of this equation that satisfies the decay condition at infinity and the boundary condition (which follows from averaging the boundary condition (59)) is zero solution. Thus,* in the leading order the boundary layer at the inner sphere is purely oscillatory*. This fact implies that the boundary condition for at (i.e., obtained by averaging the condition (47)) is . Similarly, averaging the condition (46) yields . Thus, we have
The normal velocity is determined from (33):
Here the constant of integration is chosen so as to guarantee that decays as . Evidently, , given by (64), does not satisfy the boundary condition . Now we recall that the correction to is , and, therefore, gives us the boundary condition for the next approximation of the outer solution. Indeed, according to (50) (For ), we must have

*Outer Sphere.* Consider now (38)–(40). It follows from (39) and the condition of decay at infinity (in the boundary layer variable ) for that . Hence, we have
Boundary condition for at follows from (49):
The solution of (66) that satisfies (67) and the decay condition (45) are given by
As mentioned before, the radial velocity is determined from the incompressibility condition (40):
Here again the constant of integration is chosen so as to guarantee that decays as . gives us the boundary condition for the next approximation of the outer solution:
The same arguments as in the case of the inner sphere lead to a conclusion that* in the leading order the boundary layer at the outer sphere is purely oscillatory*. This fact implies that the boundary condition for at (i.e., obtained by averaging the condition (49)) is . Similarly, averaging the condition (48) yields . Hence,

*Averaged Outer Flow.* Averaging (27) and (28) (for ) together with boundary conditions (63) and (71) implies that , that is,* there is no steady streaming in the leading order of the expansion*. First order solution.

##### 4.2. Higher Order Solutions

Using the same procedure as for the leading order equation we get the steady streaming at higher orders of . Equations (20)-(21) can be written as follows: To rewrite our asymptotic expansion in terms of the stream function, we get where is such that and for and where is defined as Taking the average of (73), we get where Similarly we get where where

#### 5. Stokes Drift

It is well-known that in oscillatory flows the observed averaged Lagrangian velocity differs from the Eulerian velocity by the term known as Stokes drift. The velocity observed in the experiments is the velocity of fluid particles, that is, the Lagrangian velocity. Our asymptotic expansion for the averaged Eulerian velocity has the form It is shown in Appendix B that the Lagrangian velocity of fluid particles is given by where is the Stokes drift velocity of the fluid particles. Comparing (80) with (82) and (81) with (83), we observe that the Stokes drift eliminates (i) from (80) and (ii) the term in (83). It also results in the additional term in the expansion of azimuthal velocity. Thus, the steady boundary layer at the inner cylinder disappears when we take account of the Stokes drift. This is a consequence of the fact that the steady Lagrangian velocity rather than the steady Eulerian velocity is invariant to the change of reference frame. In the reference frame fixed in the inner cylinder, we would have no steady boundary layer at the inner cylinder. Further calculations with the help of the known formula for , , and show that can be written as follows: where

*Asymptotic Expansion for Stream Function.* To rewrite our asymptotic expansion of terms of the averaged stream function, we get
In the last formula, and are obtained from (85):
is the stream function for the third-order Lagrangian velocity:
where

#### 6. Results Discussion

Equation (77) together with (87) and (88) represent nonzero terms in the asymptotic expansion of the stream function for the Lagrangian velocity. Let us first discuss the domain of applicability of formula (86). Steady streaming due to the vibration of the sphere is calculated. Our asymptotic expansion is formally valid for and for . Also we may expect that it will be valid for all and such that the contribution of the term to the right side of (86) is smaller than the contribution of the term. It is convenient to rewrite (86) in the form where Consider now the following quantity: which measures the magnitude of the second nonzero term relative to the first term. We expect that our theory will work for all , such that , and the smaller the is the better the theory should work.

In (86) if we take we get the results of steady streaming due to single sphere in infinite viscous incompressible fluid. The results of single vibrating sphere in an infinite fluid are Wang [10] studied the steady streaming when the sphere is at rest and fluid is oscillating at infinity for . However, the averaged stream function was incomplete because of the absence of the terms and which appears at in the expansion of stream function. Riley [11] studied the flow produced by a fixed sphere placed in an oscillating fluid using matched asymptotic expansion method and obtained a uniformly valid expansion of the stream function under the same assumption as in the present paper. He did not compute the , which appears in the expansion of stream function at . In order to obtain the invariant velocity field one should consider the Lagrangian velocity which is different from the Eulerian velocity by Stokes drift velocity of the fluid particles. In order to support the argument, here the results for steady streaming are also presented for the case when the sphere is at rest and fluid is oscillating at infinity using the Vishik-Lyusternik method. In Eulerian coordinates system the steady streaming is given as follows: After computing the Stokes drift the results of (93) are achieved.

Typical stream line pattern due to vibrating sphere is shown in Figure 1 for , , and . Graph of streamlines produced near the inner sphere produced by the formula (86) are shown in Figure 2 for and with outer sphere at . In Figure 3 two sets of circulations with opposite direction can be seen in each quadrant for , , and , respectively. Figure 3 shows that with the increase in the thickness of the inner circulation also increases. Similarly if we fix and change , it is observed that with the increase of the thickness of inner circulation also increases and outer circulations become thinner. Figure 3 shows the streamlines for and , respectively. The center of the circulation lies on the line at angle . It is observed that as and increase the stagnation point moves away from the surface of inner sphere. Figure 2 is in good comparison with the experimental results of Kotas et al. [13–15]. A recent paper [23] claims that distance of stagnation point to the surface of the sphere along the axis of oscillation is different from the distance of stagnation point to the surface of sphere along perpendicular axis. According to our results this does not happen. In Figure 4 the distance of stagnation points UPS (upper stagnation point) and LSP (left stagnation point) from the surface of the sphere is the same. Typical profiles of the radial velocity are shown in Figure 5 for and , respectively.

#### 7. Conclusion

(I)Vishik-Lyusternik method is successfully applied to study the steady streaming between two spheres. In comparison with the method of matched asymptotic expansions, it has two advantages: (i) it does not require the procedure of matching the inner and outer expansions and (ii) the boundary layer part of the expansion satisfies the condition of decay at infinity (in boundary layers) in all the orders of the expansion, which is not the case in the method of matched asymptotic expansion.(II)Stokes drift is calculated not for the outer flow but for the entire flow domain including the boundary layers.(III)We also calculated the averaged lagrangian velocity which remains invariant with the change of reference frame and is observed in experiment.(IV)The distance of stagnation points from the surface of sphere in each quadrant remains the same.(V)Qualitatively the graphs of the streamlines are very much similar to the experimental results which are done for ; our theory is valid for the small .

#### Appendix

#### A. Explicit Expression for the Right Hand Side of (27)–(29), (34)–(36), (41)–(43), and (50)-(51)

Functions , , and (29) for in (27)-(28) are given by Functions , , and for in (34)–(36) are given by Functions , , and for in (41)–(43) are given by