# Swirling Flow in a Permeable Tube at Slowly Expanding and Contracting Wall

Authors: Sufian Munawar
Journal: Mathematical Problems in Engineering
 (2017)

---

## Abstract

The swirling flow inside a circular elastic tube with expanding and contracting permeable wall in the presence of a uniform magnetic field is studied analytically. The tube is also assumed to be rotating around its axis with an angular velocity. The governing equations for this multidimensional flow are reduced to nonlinear differential equations with similarity transformations. An analytic series solution is obtained by homotopy analysis method (HAM). The effects of physical parameters on various flow characteristics, such as the velocity, skin friction, and pressure variation, have been analysed briefly. The impact of surface expansion/contraction and rotation has been investigated on the internal boundary-layer flow inside the tube of uniform cross-section.


---

## Body

## 1. Introduction

The laminar boundary-layer flows inside elastic pipes or tubes have received great attention by numerous researchers since they are widely used in various engineering and biological applications. The examples include the binary gas diffusion, transport of biological liquids through contracting or expanding vessels, pulsating porous diaphragms, urine flow in urethras, natural transpiration and cooling, and the regression of the burning surface in solid rocket motors. Uchida and Aoki [1] initiated the study on the flow inside a tube with deformable wall and one closed end. Ohki [2] discussed the laminar incompressible flow in a semi-infinite porous pipe with expanding and contracting radius in axial direction. A theoretical investigation was performed by Goto and Uchida [3] to study the suction/injection effects in an expanding/contracting pipe. Bujurke et al. [4] investigated the same flow problem by calculating a series solution. Recently, Si et al. [5] calculated the multiple solutions for the contracting or expanding porous pipe at large suction Reynolds number using singular perturbation method. Srinivas et al. [6] obtained analytic series solution for the thermal-diffusion and diffusion-thermo effects in a channel with moveable permeable walls. Currently, the effect of mass transfer and chemical reaction inside a porous tube with expanding or contracting wall was addressed by Srinivas et al. [7]. The effect of MHD on the same type of flow phenomenon was studied by Srinivas et al. [8]. The flow inside a twisted pipe was investigated by O’Dea and Waters [9] and calculated the solute uptake inside a twisting pipe. Makinde [10] presented a mathematical model describing the flow inside a collapsible tube and obtained an analytic series solution. The effect of variable viscosity and viscous dissipation on the flow in a moving pipe was discussed by Makinde [11]. Xinhui et al. [12] investigated the slip effects in the micropolar fluid inside a porous channel with expanding/contracting walls. Another type of three-dimensional flows in polar coordinates is the boundary-layer flow over cylinder. In this regard, Sprague and Weidman [13] investigated the boundary-layer flow over a twisted cylinder and an asymptotic solution for large range of the Reynolds number was obtained. Keeping this fact in mind the boundary-layer flow over cylinder can also be considered three-dimensional if the surface of pipe moves in two lateral directions. Recently, Munawar et al. [14] made a complete thermodynamical analysis of a three-dimensional viscous flow over a stretching and rotating cylinder. Fang and Yao [15] and Fang et al. [16] discussed other sorts of the three-dimensional boundary-layer flows over cylinder which was rotating, expanding, and stretching with respect to its length.This analysis extends the idea by embarking on the impact of swirling on the flow inside a contracting/expanding semi-infinite porous tube. Since most of the biological fluids are electrically conducting, therefore a uniform magnetic field is taken along the radial direction. An approximate series solution is obtained by using homotopy analysis method (HAM). The convergence and accuracy of the solution are illustrated by plotting convergence tables, plotting theh-curves, and calculating the squared residual errors. To validate our HAM solution a comparison has also been made with the results reported by Li et al. [17] using numerical and analytic techniques.


## 2. Mathematical Formulation

Consider an axisymmetric flow of an incompressible viscous fluid inside a semi-infinite horizontally placed expanding and contracting tube. The tube is rotating with a time-dependent angular velocityω(t). The radius of tube R(t) is a continuous function of time and the surface of tube is assumed to be uniformly permeable. One end of the tube is kept closed and the wall of the tube expands and contracts with a time-dependent rate R˙(t). A uniform magnetic field with intensity B0 is applied along the radial direction. Assuming small magnetic Reynolds number, the induced magnetic field is assumed to be negligible. The x-axis is measured along the axis of tube and the r-axis is perpendicular to the surface (see Figure 1). Under these assumptions the governing equations of motion are given by(1)∂u∂x+1r∂rv∂r=0,∂u∂t+u∂u∂x+v∂u∂r=-1ρ∂p∂x+ν∂2u∂x2+1r∂∂rr∂u∂r-σB02ρu,∂v∂t+u∂v∂x+v∂v∂r=-1ρ∂p∂r+ν∂2v∂x2+1r∂∂rr∂v∂r-vr2,∂w∂t+u∂w∂x+v∂w∂r+vwr=ν∂2w∂r2+1r∂w∂r-wr2+∂2w∂x2-σB02ρw.Figure 1
Schematic diagram and the coordinate system.Subject to the boundary conditions(2)ux,R=0,vx,R=-V=-V0R˙,wx,R=ωR,∂u∂rr=0=0,vx,0=0,wx,0=0,where ρ is the density, p is the pressure, ν is the kinematic viscosity, B0 is the magnetic field, V0=V/R˙ is the measure of wall permeability (V0>0 corresponds to the injection and V0<0 to suction), and u, v, and w are the velocity components in x-, r-, and θ-directions, respectively.To normalize the governing equations (1)-(2) the following dimensionless quantities are used:(3)ψ=νxfη,η=rR,p=ρν2R2P,u=νxR2ηf′,v=-νRηf,w=ωRg.Substituting (3) into (1)-(2), we get the following set of nonlinear differential equations:(4)1ηddηηddη1ηf′+αf′′+f′η-M2ηf′-1η2f′2-ff′′+ff′η=1ξ∂P∂ξ,(5)-ddη1ηf′-αf′+1η3f2-1η2ff′+Reω2g=dPdη,(6)ddη1ηddηηg-ΩReω+M2g-αg-ηg′+1ηfg′+1η2fg=0.Subject to the boundary conditions(7)limη→0ddη1ηf′=0,f0=0,g0=0,f′1=0,f1=αV0=Re,g1=1,where prime denotes the differentiation with respect to η, α=RR˙/ν is the wall expansion parameter, such that its positive values indicate the tube contraction and negative values indicate expansion, M2=σB02R2/μ is the Hartmann number, Reω=ωR2/ν is the rotation Reynolds number, Re=RV/ν is the injection Reynolds number, Ω=ω˙/ω2 is the rotation parameter, and ξ=x/R is the dimensionless axial coordinate.Equations (4) and (5) can be further simplified by cross differentiation in order to obtain a fourth-order nonlinear differential equation given by(8)ddη1ηddηηddη1ηf′=f′η2f′′-f′η-αf′′′+f′′η-f′η2-fη2f′′′-3f′′η+3f′η2+M2f′′η-f′η2.The axial and tangential shear stressesτrz and τrθ have the forms(9)τxr=μ∂u∂rr=R=μνξR2f′′1,τθr=μ∂w∂r-wrr=R=μωRg′1-μω.The total shear stressτw at the surface of tube is(10)τw=τxr2+τθr2,and the skin friction coefficient is(11)Cf=τwρωR2=ξReω2f′′12+Reω2ξ2g′1-12.


## 3. Solution of the Problem

In order to obtain the solution of the nonlinear ordinary differential equations (6) and (8) along with the boundary conditions (7) the homotopy analysis method (HAM) is used for which the following initial guesses are chosen:(12)f0η=αV02-η2η2,g0η=η,and the auxiliary linear operators are(13)L1≡ddη1ηddηηH,L2g=Hηg,H≡ddη1ηddη.Since the HAM is well known and is broadly used by the community of fluid dynamics, therefore it seems appropriate to conceal the procedural details of the method. However, the readers are referred to Liao [18] for details about procedure of the HAM solution. We here directly write the mth-order deformation equation (for m≥1) given by(14)L1fmη-χmfm-1η=ℏ1R1mη,L2gmη-χmgm-1η=ℏ2R2mη.Subject to the boundary conditions(15)fm0=limη→0Hfm=gm0=fm1=fm′1=gm1=0,where(16)R1mη=HηHfm-1+αfm-1′′′+fm-1′′η-fm-1′η2-M2ηfm-1′′-fm-1′η+1η2∑k=0m-1fm-1-kfk′′′-3fk′′η+3fk′η2-1η2∑k=0m-1fm-1-k′fk′′-fk′η,R2mη=Hηgm-1-ΩReω+M2gm-1-αgm-1-ηgm-1′+1η∑k=0m-1fm-1-kgk′+1ηgk,χm=0,form≤1,1,form>1.For a particular value of “m”, (14) are ordinary differential equations which can easily be solved by a standard method. The final solution can be written in the form of an infinite series of functions:(17)fη=f0η+∑m=1∞fmη,gη=g0η+∑m=1∞gmη,where fm(η) and gm(η) are the solutions of mth (m≥1) order deformation equations.


### 3.1. Convergence of the Series Solution

In order to verify that the series (17) is the approximate solutions of the differential equations (6) and (8) along with boundary conditions (7) it is necessary to show the convergence of the solution series. As mentioned by Liao [18] the convergence of the solution series strongly depends upon the auxiliary parameters ℏ1 and ℏ2 once the initial guesses and the linear operators have been selected. A rough estimate for the permissible regions of the values of ℏ1 and ℏ2 can be obtained by plotting the so-called ℏ-curves. In Figure 2 the ℏ-curves are plotted by keeping the other parameters fixed. From the figure it can be seen that the interval of allowed values of ℏ1 and ℏ2 is between −1.7 and −0.2 (roughly). However, to obtain the optimal values of ℏ1 and ℏ2 it is convenient to calculate the squared residual error. In Table 1 we have calculated the squared residual error at different orders of approximation for fixed values of the other parameters. From the table it is clear that as the order of approximation increases the squared residual errors decrease and are confined in an admissible range. The corresponding optimal values for ℏ1 and ℏ2 can be obtained by minimizing the squared residual error using the command “NMinimize” of the software Mathematica. In Table 1 we also calculate the HAM solution at different orders and corresponding optimal values of ℏ1 and ℏ2. Clearly, there are no corrections after 10th order of approximation up to 8 decimal places. This proves the convergence and accuracy of our series solution. A comparison of the current analytic solution is given with the numerical and asymptotic solution reported by Li et al. [17] in Table 2. The table shows a good agreement up to 4 decimal places between our results and the numerical results provided in Li et al. [17].Table 1
Convergence and accuracy of HAM solution at different orders of approximation whenα=1, M=0.5, V0=0.5, Ω=1, ξ=1, and N=20.
Order of approximationm f ′ ′(1) g ′(1) Residual errorE1m for (8) Residual errorE2m for (6) Optimalℏ1 Optimalℏ2 2 −3.60185569 0.98963266 5.756 × 10 - 2 3.832 × 10 - 3 −0.897 −0.916 5 −3.60529192 0.98947505 1.006 × 10 - 4 3.256 × 10 - 6 −0.948 −0.992 10 −3.60529631 0.98947481 1.774 × 10 - 9 4.968 × 10 - 11 −0.932 −0.997 15 −3.60529631 0.98947481 3.615 × 10 - 11 1.024 × 10 - 12 −0.845 −0.832 20 −3.60529631 0.98947481 3.561 × 10 - 10 4.223 × 10 - 11 −0.689 −0.865 25 −3.60529631 0.98947481 1.808 × 10 - 9 1.753 × 10 - 11 −0.580 −0.621 30 −3.60529631 0.98947481 3.981 × 10 - 9 3.781 × 10 - 11 −0.478 −0.518Table 2
Comparison of current analytic solution by HAM with the solution reported by Li et al. [17] for f′′(1)/4Re at α=0.05.
Re Asymptotic solution by Li et al. [17] Numerical solution by Li et al. [17] Current HAM solution 0.592 −2.0841 −2.0694 −2.0694 0.297 −2.0339 −2.0299 −2.0299 0.041 −1.9903 −1.9903 −1.9903 −0.052 1.9745 1.9744 1.9744 −0.141 1.9593 1.9584 1.9584 −0.422 1.9115 1.9015 1.9015Figure 2
ℏ-curves at 30th order of approximation.


## 3.1. Convergence of the Series Solution

In order to verify that the series (17) is the approximate solutions of the differential equations (6) and (8) along with boundary conditions (7) it is necessary to show the convergence of the solution series. As mentioned by Liao [18] the convergence of the solution series strongly depends upon the auxiliary parameters ℏ1 and ℏ2 once the initial guesses and the linear operators have been selected. A rough estimate for the permissible regions of the values of ℏ1 and ℏ2 can be obtained by plotting the so-called ℏ-curves. In Figure 2 the ℏ-curves are plotted by keeping the other parameters fixed. From the figure it can be seen that the interval of allowed values of ℏ1 and ℏ2 is between −1.7 and −0.2 (roughly). However, to obtain the optimal values of ℏ1 and ℏ2 it is convenient to calculate the squared residual error. In Table 1 we have calculated the squared residual error at different orders of approximation for fixed values of the other parameters. From the table it is clear that as the order of approximation increases the squared residual errors decrease and are confined in an admissible range. The corresponding optimal values for ℏ1 and ℏ2 can be obtained by minimizing the squared residual error using the command “NMinimize” of the software Mathematica. In Table 1 we also calculate the HAM solution at different orders and corresponding optimal values of ℏ1 and ℏ2. Clearly, there are no corrections after 10th order of approximation up to 8 decimal places. This proves the convergence and accuracy of our series solution. A comparison of the current analytic solution is given with the numerical and asymptotic solution reported by Li et al. [17] in Table 2. The table shows a good agreement up to 4 decimal places between our results and the numerical results provided in Li et al. [17].Table 1
Convergence and accuracy of HAM solution at different orders of approximation whenα=1, M=0.5, V0=0.5, Ω=1, ξ=1, and N=20.
Order of approximationm f ′ ′(1) g ′(1) Residual errorE1m for (8) Residual errorE2m for (6) Optimalℏ1 Optimalℏ2 2 −3.60185569 0.98963266 5.756 × 10 - 2 3.832 × 10 - 3 −0.897 −0.916 5 −3.60529192 0.98947505 1.006 × 10 - 4 3.256 × 10 - 6 −0.948 −0.992 10 −3.60529631 0.98947481 1.774 × 10 - 9 4.968 × 10 - 11 −0.932 −0.997 15 −3.60529631 0.98947481 3.615 × 10 - 11 1.024 × 10 - 12 −0.845 −0.832 20 −3.60529631 0.98947481 3.561 × 10 - 10 4.223 × 10 - 11 −0.689 −0.865 25 −3.60529631 0.98947481 1.808 × 10 - 9 1.753 × 10 - 11 −0.580 −0.621 30 −3.60529631 0.98947481 3.981 × 10 - 9 3.781 × 10 - 11 −0.478 −0.518Table 2
Comparison of current analytic solution by HAM with the solution reported by Li et al. [17] for f′′(1)/4Re at α=0.05.
Re Asymptotic solution by Li et al. [17] Numerical solution by Li et al. [17] Current HAM solution 0.592 −2.0841 −2.0694 −2.0694 0.297 −2.0339 −2.0299 −2.0299 0.041 −1.9903 −1.9903 −1.9903 −0.052 1.9745 1.9744 1.9744 −0.141 1.9593 1.9584 1.9584 −0.422 1.9115 1.9015 1.9015Figure 2
ℏ-curves at 30th order of approximation.


## 4. Flow Characteristics

In this section we analyse the flow characteristics of the problem by plotting the graphs of velocity profile, wall skin friction, and the pressure variations in the tube for various values of the parameters. The axial velocity profile versusη for various values of wall expansion parameter α can be seen in Figure 3. From the figure it is clear that as contraction increases the velocity profile also increases. It is further observed that the fluid velocity decreases as the tube wall shifts from contraction to expansion state. It can be noticed that as the tube either expands or contracts, the magnitude of velocity increases. However, for increasing expansion rate the backward flow accelerates, whereas, for increasing contraction, the forward flow exaggerates. Figure 4 has been prepared to express the variation in the tangential velocity of fluid induced by the expansion and contraction of tube wall. This figure shows that for increasing expansion rate the tangential velocity profile decreases, while for increasing the contraction rate the tangential velocity increases.Figure 3
Effect of wall expansion parameterα on the axial velocity.Figure 4
Effect of wall expansion parameterα on the tangential velocity.Figures5 and 6 are plotted to discuss the effects of the Hartmann number M on the axial and tangential velocity profiles, respectively. Figure 5 indicates that the axial velocity is large in case of hydrodynamic flow as compared to the magnetohydrodynamic flow near the tube centre. However, a reverse behaviour can be seen near the tube boundary. Here axial velocity is going to be augmented when large values of M are taken into contemplation. Figure 6 reveals that the tangential velocity of fluid improves by increasing the values of M throughout the tube. In this case magnetic field acts as an opposing force which resists the tangential velocity.Figure 5
Effect of the Hartmann numberM on the axial velocity.Figure 6
Effect of the Hartmann numberM on the tangential velocity.Figure7 shows the combined effects of wall expansion/contraction parameter α and the Hartmann number M on the wall skin friction coefficient. A close look at this Figure informs that, for the state of tube expansion, the magnetic field works as a supporting force to enhance the skin friction coefficient, while, for tube contraction state, the magnetic field becomes a resistive force. Whether the tube is in expansion state or in contraction state, the magnitude of skin friction always increases as M increases. The value α=0 corresponds to the rest state (when the tube neither contracts nor expands). In the rest state of tube, the value of skin friction is zero which indicates no flow across the tube. The influence of permeability parameter V0 on the axial and tangential velocities is presented by Figures 8 and 9. These figures are portrayed when the tube is observed to be in contracting state (α=2). Figure 8 renders that the axial velocity of fluid improves as one moves from the suction regime to the injection. Moreover, the variations in velocity are more noticeable near the tube centre. Figure 9 demonstrates that the tangential velocity of the fluid flowing through the tube increases by increasing injection and of course the tangential velocity profile decreases in the case when the suction phenomenon becomes strong.Figure 7
The skin friction coefficient againstα for various values of M.Figure 8
Effect of the suction/injection parameterV0 on the axial velocity.Figure 9
Effect of the parameterV0 on the tangential velocity.The skin friction coefficient for various values of permeability parameterV0 against α has been plotted in Figure 10. Clearly, two different variation trends are noticed for various values of V0. From the figure it is clear that for the case of tube contraction the skin friction becomes a decreasing function of V0, whereas the skin friction coefficient is going to be enhanced by increasing V0 in expansion state. This means that when tube contracts the presence of strong injection strengthens the negative skin friction or in other words, injection causes the forward flow to be enhanced when cylinder contracts. The consequence of rotational Reynolds number Reω on the tangential velocity has been captured in Figure 11 when the tube is in contracting state. This figure depicts that the fluid tangential velocity decreases when large values of Reω are taken into account. This conduct of fluid velocity is of course due to the presence of large inertial forces which resist the fluid tangential velocity.Figure 10
The skin friction coefficient againstα for various values of V0.Figure 11
Effect of the rotational Reynolds numberReω on the tangential velocity.The variations in the pressure versusη for different emerging parameters have been shown through Figures 12–15. In this regard, the selected parameters of concern are α, M, V0, and Reω. Figure 12 notifies that the pressure gradient increases and becomes adverse near the tube walls when tube is going to be relaxed back from the expansion state. Also in different phases of expansion and contraction states the pressure gradient shows miscellaneous variations throughout the tube. A careful glance at this figure reveals that the pressure variation curve becomes more concave when tube is going to be contracted more. At α=2 it is noticed that the pressure gradient becomes favourable near the centre of tube while near the wall it becomes more adverse. In Figure 13, it is observed that the pressure gradient becomes more favourable for large values of M in the region 0≤η<0.86 (roughly), although, in the region 0.86<η≤1, the magnetic force accelerates the adverse pressure gradient. Figure 14 reports that the variation in the values of permeability parameter V0 has miscellaneous effects on the pressure gradient. In this figure, it is noted that the curve for the pressure gradient changes from linear to parabolic shape as we alter from suction to injection situation. For the case of remarkable injection effects, the favourable pressure gradient increases near the tube centre when compared to its behaviour near the tube wall. Figure 15 explains that the adverse pressure gradient increases for large values of Reω throughout the tube. The variation in pressure accumulates near the tube centre and becomes more noticeable near the tube boundary.Figure 12
Pressure variation in the tube for different values of the parameterα.Figure 13
Pressure variation in the tube for different values of the Hartmann numberM.Figure 14
Effect of the suction/injection parameterV0 on pressure variation in the tube.Figure 15
Effect of the rotation Reynolds numberReω on pressure variation in the tube.


## 5. Conclusion

The magnetohydrodynamic rotating flow inside a uniformly porous tube with contracting/expanding wall is investigated. The analysis is based on the analytic technique HAM with the help of which we obtained an approximate series solution. The solution series is uniformly valid for all range of the involved parameters. The physical results categorically indicate the following findings:(i)
The magnitude of skin friction increases as the expansion/contraction rates increase.(ii)
The contraction of tube accelerates the torsional velocity of fluid particles while expansion of tube depreciates the rotating velocity of fluid.(iii)
The axial velocity has two variation trends as magnetic field intensifies; near the centre, the velocity is a decreasing function of the Hartmann number and near the wall the velocity is an increasing function of magnetic field.(iv)
Also the magnitude of skin friction increases for both expansion and contraction as the magnetic field becomes stronger.(v)
The pressure variation is large near the wall of pipe and vanishes near the centre of pipe.(vi)
The large contraction rate results in the favourable pressure gradient near the centre of the tube while at the wall adverse pressure gradient magnifies.(vii)
The augmentation of magnetic field results in the increasing of favourable pressure gradient in tube centre while it increases the adverse pressure gradient near the wall.(viii)
An increase in the Reynolds numberReω is to increase the adverse pressure gradient.



---

*Source: 6904161-2017-01-19.xml*