# Suction/Injection Effects on the Swirling Flow of a Reiner-Rivlin Fluid near a Rough Surface

Authors: Bikash Sahoo; Sébastien Poncet; Fotini Labropulu
Journal: Journal of Fluids
 (2015)

---

## Abstract

The similarity equations for the Bödewadt flow of a non-Newtonian Reiner-Rivlin fluid, subject to uniform suction/injection, are solved numerically. The conventional no-slip boundary conditions are replaced by corresponding partial slip boundary conditions, owing to the roughness of the infinite stationary disk. The combined effects of surface slip (λ), suction/injection velocity (W), and cross-viscous parameter (L) on the momentum boundary layer are studied in detail. It is interesting to find that suction dominates the oscillations in the velocity profiles and decreases the boundary layer thickness significantly. On the other hand, injection has opposite effects on the velocity profiles and the boundary layer thickness.


---

## Body

## 1. Introduction

The problem of Newtonian and non-Newtonian swirling flows near a rotating or stationary disk has occupied a central position in the field of fluid mechanics due firstly to the fact that similarity solutions to the Navier-Stokes equations may be found in some idealized infinite configurations and secondly to its industrial and technical applications in rotating machinery (centrifugal pumps, turbines, or computer storage devices), chemical engineering (spinning disk reactors, crystal growth processes, or rheometers), or oceanography among other things.Recently, Sahoo [1] and Sahoo and Poncet [2] have obtained numerical solution to similarity equations arising due to steady revolving flow (known as Bödewadt flow [3]) of a non-Newtonian Reiner-Rivlin fluid near an infinite rough stationary disk. In this short note, the flow problem studied by Sahoo and Poncet [2] has been reconsidered, including uniform suction/injection at the surface of the stationary disk. Knowledge of the flow structure close to a porous disk is of practical significance with regard to problems of lubrication of porous bearings or gaseous diffusion among other things. Suction or injection at the surface of a porous disk is also commonly used in chemical engineering to increase the electrochemical reaction time during electrolytic processes [4] or for control purpose as, under given conditions, it delays the transition to turbulence [5]. There are only few attempts in the literature to consider suction/injection at the disk surface. Kelson and Desseaux [6] revisited the von Kármán flow problem over a rotating disk including mass transfer through the disk. Attia [7] extended their work by considering the unsteady flow over an infinite rotating disk with uniform suction and injection and heat transfer effects. Ashraf et al. [8] considered the flow of a micropolar fluid between two stationary disks with constant injection velocity at the surface of one disk. Domairry and Aziz [9] investigated by a homotopy perturbation method the MHD flow between two parallel stationary disks with suction or injection through one of the two disks.None of these previous works considered the Bödewadt flow of a non-Newtonian fluid over a stationary rough disk with mass transfer through it. This paper is then an endeavour to fill this gap. A second-order finite difference method has been adopted here to solve the resulting system of fully coupled and highly nonlinear similarity equations arising due to Reiner-Rivlin swirling flow over an infinite stationary porous disk. The objective is to check if suction or injection is an effective way to reduce the chances of separation of the boundary layer.


## 2. Formulation of the Problem

One considers an incompressible non-Newtonian Reiner-Rivlin fluid, whose constitutive equation is given by(1)
T
=
-
p
I
+
ϕ
1
D
+
ϕ
2
D
2
,
where D
=
(
1
/
2
)
[
∇
v
+
(
∇
v
)
T
] [10]. The response functions ϕ
1 and ϕ
2 are functions of the scalar invariants (
tr
⁡
D
)
2, tr
⁡
(
D
2
), and det
⁡
D [11]. The fluid occupies the space z
>
0 over an infinite stationary disk, which coincides with z
=
0. The motion is due to the rotation of the fluid like a rigid body with constant angular velocity Ω at large distance from the stationary disk. Let v
=
(
u
,
v
,
w
) be the fluid velocity vector in a (
r
,
θ
,
z
) stationary reference frame (see [2] for a schematic view of the flow configuration). Using the von Kármán transformations [12],
(2)
u
=
r
Ω
F
(
ζ
)
,
v
=
r
Ω
G
(
ζ
)
,
w
=
Ω
ν
H
(
ζ
)
,
z
=
ν
Ω
ζ
,
p
ρ
=
-
ν
Ω
P
(
ζ
)
+
1
2
Ω
2
r
2
,
where ν is the fluid kinematic viscosity. By considering the usual boundary layer approximations, the equations of continuity and motion take the following forms [1, 2]:
(3)
d
H
d
ζ
+
2
F
=
0
,
(4)
d
2
F
d
ζ
2
-
H
d
F
d
ζ
-
F
2
+
G
2
-
1
2
L
d
F
d
ζ
2
-
3
d
G
d
ζ
2
-
2
F
d
2
F
d
ζ
2
=
1
,
(5)
d
2
G
d
ζ
2
-
H
d
G
d
ζ
-
2
F
G
+
L
d
F
d
ζ
d
G
d
ζ
+
F
d
2
G
d
ζ
2
=
0
,
(6)
d
2
H
d
ζ
2
-
H
d
H
d
ζ
-
7
2
L
d
H
d
ζ
d
2
H
d
ζ
2
+
d
P
d
ζ
=
0
,
where L
=
ϕ
2
Ω
/
ϕ
1 corresponds to the non-Newtonian cross-viscous parameter. The above system of equations has to be solved subject to following partial slip boundary conditions:
(7)
F
(
0
)
=
λ
[
F
′
(
0
)
-
L
F
0
F
′
(
0
)
]
,
G
(
0
)
=
η
[
G
′
(
0
)
-
2
L
F
0
G
′
(
0
)
]
,
H
(
0
)
=
W
,
F
(
∞
)
⟶
0
,
G
(
∞
)
⟶
1
,
where λ and η are nondimensional slip coefficients and W
=
W
0
/
Ω
ν is the uniform suction (W
<
0) or injection (W
>
0) velocity.The expression of the nondimensional moment coefficientC
m is given by
(8)
C
m
=
-
π
[
G
′
(
0
)
-
2
L
F
0
G
′
(
0
)
]
R
,
where R
=
R
2
Ω
/
ν is the Reynolds number based on the disk radius and the fluid velocity far from the disk surface. C
m represents the torque required to maintain the disk at rest.


## 3. Numerical Solution of the Problem

In this section, we will present briefly the finite difference method that has been used to solve the system of coupled, nonlinear equations (3)–(5) subject to slip boundary conditions (7). It is customary to mention that similar scheme has been used by Sahoo et al. [13] to solve the Bödewadt flow problem for a viscous fluid with Navier's slip boundary conditions. In this problem, as H
0
≠
0, a slightly modified scheme has been used in order to get diagonally dominant matrix while using generalized Gauss-Seidel method. The semi-infinite integration domain [
0
,
∞
) is replaced by a finite domain [
0
,
ζ
∞
). In practice, ζ
∞ should be chosen to be sufficiently large so that the numerical solution closely approximates the terminal boundary conditions and takes into account the asymptotical behavior far from the disk (see Appendix 1 in [14] for the asymptotical behavior of the solutions for large ζ). One approximates the functions and their derivatives by their finite difference counterparts and eventually solves a sequence of linear systems as explained below.(1)
One first solves(9)
1
+
L
F
k
F
′′
-
H
k
F
′
=
F
k
2
-
G
k
2
+
1
2
L
F
′
k
2
-
3
G
′
k
2
+
1
using mixed boundary conditions (7) and calls the solution of (9) as F
~
(
k
+
1
), with k being the iteration index. To obtain convergence, one defines F
(
k
+
1
) by the smoothing formula:
(10)
F
(
k
+
1
)
=
α
1
F
~
(
k
+
1
)
+
(
1
-
α
1
)
F
~
k
,
0
≤
α
1
≤
1
.
(2)
Then, one solves(11)
1
+
L
F
k
+
1
G
′′
+
L
F
′
k
+
1
-
H
k
G
′
=
2
F
(
k
+
1
)
G
(
k
)
using derivative boundary conditions (7) and calls the solution of (11) as G
~
(
k
+
1
). To obtain convergence, one defines G
(
k
+
1
) by the following smoothing formula:
(12)
G
(
k
+
1
)
=
α
2
G
~
(
k
+
1
)
+
(
1
-
α
2
)
G
~
k
,
0
≤
α
2
≤
1
.
(3)
In this step, one solves(13)
H
′
=
-
2
F
(
k
+
1
)
and calls the solution asH
~
(
k
+
1
). To obtain convergence, one defines H
(
k
+
1
) by the smoothing formula:
(14)
H
(
k
+
1
)
=
α
3
H
~
(
k
+
1
)
+
(
1
-
α
3
)
H
~
k
,
0
≤
α
3
≤
1
.
(4)
The iterations start with suitable initial guessesF
(
0
)
,
G
(
0
)
,
H
(
0
)
,
F
′
(
0
), and G
′
(
0
), borrowed from the work by Sahoo and Poncet [2]. If F
(
k
+
1
)
,
F
(
k
)
,
G
(
k
+
1
)
,
G
(
k
), and H
(
k
+
1
)
,
H
(
k
) are close enough to each other, iterations are stopped; otherwise one sets k
=
k
+
1 and goes to step (
1
).In order to solve the above system of equations by finite difference method, a uniform grid in 0
≤
ζ
≤
ζ
∞ is introduced by dividing it into n equal parts with a mesh size h equal to 0.01. One approximates the derivatives by their finite difference counterparts using second-order schemes as follows:
(15)
F
′
(
ζ
i
)
=
F
i
+
1
-
F
i
-
1
2
h
,
F
′′
(
ζ
i
)
=
F
i
+
1
-
2
F
i
+
F
i
-
1
h
2
,
i
=
1,2
,
…
,
n
-
1
.
In order to obtain a diagonally dominant linear algebraic system for (9) and (11), F
′ and G
′ are discretized by backward difference approximations as H
i
(
k
+
1
)
>
0 for Bödewadt flow. Finally, (13) is discretized by the central difference approximation. The above algebraic system of equations is solved by generalized Gauss-Seidel method [15]. The convergence of the generalized Gauss-Seidel method for the above diagonally dominant system of equations is quite fast. About 19–21 iterations are necessary to achieve an accuracy of 1
0
-
6. The Fortran 90 code was compiled and run using one of the NIT Rourkela high-end Linux servers. The typical time per iteration for a given mesh distribution (ζ
∞
=
20 and h
=
0.01) and L
=
λ
=
1 is 32.9 seconds for W
=
0 and increases up to 45.4 seconds for W
=
1.


## 4. Results and Discussions

The effects of slip (λ) and cross-viscous parameter (L) on the momentum boundary layer have been already precisely discussed in [2]. This short communication focuses only on the effects of the suction/injection velocity on the momentum boundary layer for fixed values of L
=
1 and λ
=
η
=
1.The velocity profiles for the Bödewadt problem exhibit oscillations unlike von Kármán flow. The oscillations occurring in the boundary layer when the fluid rotates near a stationary disk can be explained in the following manner. The radial inflow, induced in the vicinity of the stationary disk, tends to conserve angular momentum and thus to increase the tangential velocity when the local radius is decreased. For an overshoot, radial convection of the angular momentum near the disk must be strong enough to more than balance the dissipation of angular momentum caused by the wall shear. This inward radial convection of surplus angular momentum is possible as long as the distribution of circulation in the outer flow increases with increasing radius. A local overshoot in the tangential velocity increases the centrifugal force locally, which tends to induce a radial outflow. This radial outflow convects an angular momentum defect to force an undershoot in the tangential velocity profile, and the above process is repeated to yield oscillatory approach to infinity. It is interesting to observe that the oscillations in the three velocity components reduce asW changes its sign from positive (injection) to negative (suction), as seen in Figures 1 to 3. The boundary layer thickness decreases significantly when suction is applied. In fact, there is no oscillation in the velocity profiles for W
=
-
1. From Figure 1, it is clear that the intensity of the back flow near the disk surface decreases as W changes sign from positive to negative. Suction decreases the amplitude of the oscillations, whereas injection enhances them. To summarize, suction has a stabilizing effect on the velocity profiles, while injection destabilizes the flow. The influence of W on the components F and G confirms the previous results of Attia [7], in the case of a rotating disk in a porous medium. On the contrary, the effect of W on H has an opposite behavior. The axial velocity component H is relatively constant whatever the distance from the disk is when suction is applied. For future comparisons, the variations of H
∞, F
(
0
), and G
(
0
) for different combinations of the flow parameters are provided in Table 1. The velocity gradients close to the disk are more interesting. Suction tends to diminish them, which means that the shear stresses may be reduced by decreasing W. From Figure 2, it can be seen that the boundary layer thickness is also reduced when W decreases. With injection W
=
1, the fluid is pushed towards a larger distance from the disk, which induces a thickening of the boundary layer. The injection velocity W may be then used to adjust both the shear stress and the boundary layer thickness (and as a consequence the velocity profiles) in given engineering applications.Table 1
Variations ofH
∞, F
(
0
), and G
(
0
) with different flow parameters.
L
λ(=η)
W
H
∞
F
(
0
)
G
(
0
)
0.0
1.0
1.0
2.070589
−0.393258
0.510032
1.0
1.899414
−0.261875
0.409946
2.0
1.800522
−0.195825
0.345059
1.0
1.0
1.0
1.899414
−0.261875
0.409946
2.0
1.522734
−0.091550
0.311701
3.0
1.372901
−0.048644
0.239897
1.0
1.0
−1.0
−0.760271
−0.066407
0.593022
0.0
0.561328
−0.140015
0.497071
1.0
1.899414
−0.261876
0.409946Figure 1
Variation ofF with W.Figure 2
Variation ofG with W.Figure 3
Variation ofH with W.Finally, Figure4 shows the variation of the moment coefficient C
m with the suction/injection velocity W for three different sets of values of L and λ. Note that the Reynolds number has been fixed to R
=
1. The moment coefficient C
m remains negative for all values of W. It is interesting to observe that W has a significant effect on the moment coefficient for a viscous (L
=
0) Bödewadt flow with no-slip (λ
=
0) boundary conditions. In this case, the torque required to maintain the disk at rest is high when there is suction at the disk surface and it decreases in magnitude as W changes its sign from negative to positive. The rate of decrease of C
m as W changes sign from negative to positive for the other two cases (L
=
0
,
λ
=
1
;
L
=
1
,
λ
=
1) is not significant. It is also observed that, for the Reiner-Rivlin fluid (L
≠
0), C
m decreases in magnitude when W varies in the range [
-
1,1
] up to a critical value W
*. For L
=
1 and λ
=
1, W
*
≃
0.75.Figure 4
Variation ofC
m with W.


## 5. Conclusions

In this short communication, we have investigated the effects of suction and injection on the momentum boundary layer arising due to the swirling flow of a non-Newtonian Reiner-Rivlin fluid over an infinite rough stationary disk. A second-order finite difference method has been adopted to solve the resulting system of fully coupled and highly nonlinear similarity equations. It is observed that suction suppresses the oscillations in the velocity profiles, whereas injection enhances it. The boundary layer thickness decreases as the suction increases. Injection has an opposite effect on it.



---

*Source: 253504-2015-01-05.xml*