This paper employs the computational approach known as successive linearization method (SLM) to tackle a fourth order nonlinear differential equation modelling the transient flow of an incompressible viscous fluid between two parallel plates produced by a simple wall motion. Numerical and graphical results obtained show excellent agreement with the earlier results reported in the literature. We obtain solution branches as
well as a turning point in the flow field accurately. A comparison with numerical results generated using the inbuilt MATLAB boundary value solver, bvp4c, demonstrates that the SLM approach is a very efficient technique for tackling highly nonlinear differential equations of the type discussed in this paper.
1. Introduction
Studies related to transient flows produced by a simple wall motion have been of interest for several years due to its practical importance in understanding several engineering and physiological flow problems. For instance, the entire conduits in human body are flexible and also collapsible. That is, when the external pressure exceeds the internal pressure, the cross-sectional area can be significantly reduced, if not fully diminished. The cross-section may eventually return to its original shape when the external pressure is reduced, and, consequently, normal internal fluid flow can be restored [1]. Other applications can be found in unsteady loading, which is met frequently in many hydrodynamical machines and apparatus [2]. In the light of these applications, squeezing flow in a channel has been studied by many authors; mention may be made of research studies [3–6]. This problem admits similarity variable [7, 8], thereby reducing the unsteady Navier-Stokes equations to a parameter dependent fourth order nonlinear ordinary differential equation for the similarity function.
Generally speaking, nonlinear problems and their solutions provide an insight into inherently complex physical process in the system. The nonlinear nature of the model equations in most cases precludes its exact solution. Several approximation techniques have been developed to tackle this problem such as the homotopy analysis method [9–11], homotopy perturbation method [12, 13], spectral homotopy analysis method [14, 15], and variational iteration method [16]. In this paper, we employ the successive linearisation method [17–19] to tackle a fourth order nonlinear boundary value problem that governs the squeezing flow problem between parallel plates. In this work, we assess the applicability of the SLM approach in solving nonlinear problems with bifurcations. Such problems are very difficult to resolve numerically near the bifurcation point. Numerical and graphical results obtained using the new SLM approach are validated through comparison with numerical results generated using the inbuilt MATLAB boundary value solver, bvp4c, for different values of the governing physical parameters. In following sections, the problem is formulated, analysed, and discussed.
2. Mathematical Formulation
Consider a transient flow of an incompressible viscous fluid between parallel plates driven by the normal motion of the plates. Take a Cartesian coordinate system (x,y) where x lies in the streamwise direction and y is the distance measured in the transverse direction. Let u and v be the velocity components in the directions of x and y increasing, respectively. It is assumed that the two plates are at y=±a0(1-αt), where a0 is the position at time t=0 as shown in Figure 1.
Schematic diagram of the problem.
When α is positive, the two plates are squeezed symmetrically until they touch at t=1/α. Negative values of α represent the symmetrical separation of the plates. The length of the plates is assumed to be much larger than the gap width at any time such that the end effects could be neglected. Following [4, 6–8], the two-dimensional governing equation of motion in terms of vorticity (ω) and stream function (Ψ) formulation is given as
(1)∂ω∂t+∂(ω,Ψ)∂(x,y)=ν∇ω,ω=-∇2Ψ,
with
(2)∂Ψ∂y=0,∂Ψ∂x=-dadt,ony=a(t),∂2Ψ∂y2=0,∂Ψ∂x=0,ony=0.
We introduce the following transformations:
(3)η=ya01-αt,Ψ=αa0xF(η)2(1-αt),ω=-αx2a0((1-αt))3d2Fdη2.
Substituting (3) into (1) and (2), we obtain
(4)d4Fdη4=R(d2Fdη2dFdη-Fd3Fdη3+ηd3Fdη3+3d2Fdη2),(5)dFdη=0,F=1,onη=1,(6)d2Fdη2=0,F=0,onη=0,
where R=a02α/2ν is the local Reynolds number (R>0 represents squeezing and R<0 represents separation). The wall skin friction is given by
(7)τw=-μ∂u∂y=-αμx2a0((1-αt))3d2Fdη2,atη=1,
where μ is the dynamic coefficient of viscosity. From the axial component of the Navier-Stokes equations, the pressure drop in the longitudinal direction can be obtained. Let
(8)∂P∂x=μαxA2a02(1-αt)2,
and we obtain
(9)A=d3Fdη3-R[(dFdη)2-Fd2Fdη2+ηd2Fdη2+2dFdη].
In the following section, (4)–(6) will be solved using successive linearization method and other important flow properties like the skin friction and pressure drop will be determined.
3. Successive Linearisation Method (SLM) Approach
The proposed linearisation method of solution, hereinafter referred to as the successive linearisation method (SLM), is based on the assumption that the unknown function F(η) can be expanded as
(10)F(η)=Fi(η)+∑m=0i-1Fm(η),i=1,2,3,…,
where Fi are unknown functions. The solutions of Fm, (m=1,2,…) are obtained recursively by solving the linear part of the equation that results from substituting (10) in the governing equations (4) using F0(η) as an initial approximation. The linearisation technique is based on the assumption that Fi becomes increasingly smaller as i becomes large; that is,
(11)limi→∞Fi=0.
The initial approximation F0(η) must be chosen in such a way that it satisfies the boundary conditions (5) and (6). An appropriate initial guess is
(12)F0(η)=(e2+32)η-(3e2+12)η3+eη4,
where e is an arbitrary constant which when varied results in multiple solutions. Substituting (10) in the governing equations and neglecting nonlinear terms in Fi, Fi′, Fi′′, and Fi′′′ give
(13)Fiiv+a1Fi′′′+a2Fi′′+a3Fi′+a4Fi=ri-1,
where
(14)a1=R(∑m=0i-1Fm-η),a2=-R(3+∑m=0i-1Fm′),a3=-R∑m=0i-1Fm′′,a4=R∑m=0i-1Fm′′′,ri-1=-∑m=0i-1Fmiv+R(∑m=0i-1Fm′∑m=0i-1Fm′′-∑m=0i-1Fm∑m=0i-1Fm′′′fffffffffffffffffffffffff+η∑m=0i-1Fm′′′+3∑m=0i-1Fm′′).
Starting from the initial approximation, F0, the subsequent solutions for Fm, m≥1, are obtained iteratively by solving (13) subject the the boundary conditions
(15)Fi′=0,Fi=0,onη=1,Fi′′=0,Fi=0,onη=0.
Once each solution for Fi(i≥1) has been obtained, the approximate solutions for F(η) are obtained as
(16)F(η)≈∑m=0MFm(η),
where M is the order of SLM approximation. It is worth noting that the coefficient parameters and the right hand side of (13) for i=1,2,3,…, are known (from previous iterations). Thus, system (13) can easily be solved using numerical methods such as finite differences, finite elements, Runge-Kutta based shooting methods, or collocation methods. In this work, (13) is solved using the Chebyshev spectral collocation method. This method is based on approximating the unknown functions by the Chebyshev interpolating polynomials in such a way that they are collocated at the Gauss-Lobatto points defined as
(17)zj=cosπjN,j=0,1,…,N,
where N is the number of collocation points used (see, e.g., [20, 21]). In order to implement the method, the physical region [0,1] is transformed into the region [-1,1] using the mapping
(18)η=z+12,-1≤z≤1.
The derivative of Fi at the collocation points is represented as
(19)dsFidηs=∑k=0NDkjsFi(zk),j=0,1,…,N,
where D=((2/L)𝒟)s and 𝒟 is the Chebyshev spectral differentiation matrix (see e.g., [20, 21]). Substituting (17)–(19) in (13) results in the matrix equation
(20)Ai-1Fi=Ri-1,
in which Ai-1 is a (N+1)×(N+1) square matrix and Fi and Ri-1 are (N+1)×1 column vectors defined by
(21)Fi=[Fi(z0),Fi(z1),…,Fi(zN-1),Fi(zN)]T,Ri-1=[ri-1(z0),ri-1(z1),…,ri-1(zN-1),ri-1(zN)]T,Ai-1=D4+a1D3+a2D2+a3D+a4.
In the above definitions, ak(k=1,2,3,4) are diagonal matrices of size (N+1)×(N+1). After modifying the matrix system (20) to incorporate boundary conditions, the solution is obtained as
(22)Fi=Ai-1-1Ri-1.
4. Results and Discussion
In this section, we present the results for the solution of the governing nonlinear boundary value problem (4). To check the accuracy of the successive linearisation method (SLM), comparison is made with numerical solutions obtained using the MATLAB routine bvp4c. The MATLAB routine bvp4c is based on an adaptive Lobatto quadrature scheme [22, 23].
Table 1 gives a comparison between the 10th order SLM approximate results and the bvp4c numerical results for F′(0) and F′′(1) at selected values of R. By using different values of the constant e in the initial approximation (see (12)), it was found that both the SLM and bvp4c give multiple solutions when R<0. Two solutions, called lower branch and upper branch, were identified when -3.45≤R≤-3.25. We observe that the SLM results are in very good agreement with the bvp4c results for both branches of the solutions. Table 1 indicates that the skin friction F′′(1) decreases when R is increased in the lower branch and the opposite effect is observed when R is increased in the upper branch.
The tenth order SLM approximation for F′(0) and F′′(1) at the two branches of solutions for different values of R.
Lower branch
Upper branch
R
SLM
bvp4c
SLM
bvp4c
F′(0)
−3.45
3.39748
3.39748
5.45425
5.45425
−3.40
3.12576
3.12576
6.27172
6.27172
−3.35
2.94758
2.94758
7.05936
7.05936
−3.30
2.81355
2.81355
7.88138
7.88138
−3.25
2.70596
2.70596
8.77609
8.77609
F′′(1)
−3.45
7.28689
7.28689
19.95709
19.95709
−3.40
5.80900
5.80900
25.70501
25.70502
−3.35
4.86274
4.86274
31.60366
31.60366
−3.30
4.16178
4.16178
38.12272
38.12272
−3.25
3.60511
3.60511
45.62417
45.62417
Figure 2 illustrates a slice of bifurcation diagram in both planes. For symmetrical squeezing of the plates; that is, R>0, only one solution branch exists; this can be regarded as the lower solution branch. Another solution branch was identified in addition to the lower solution branch when the plates were symmetrically separated (R<0); this is the upper solution branch. A turning point exists between the primary and secondary solution branches at R=-3.495. This bifurcation result obtained using SLM is in perfect agreement with the one reported by Makinde et al. [5] using Hermite-Padé approximation technique. Moreover, this turning point is very significant with respect to application; it represents the symmetrical separation limit of plates during flow process.
Slice of the bifurcation diagrams for F′(0) and F′′(1) against R.
Figure 3 depicts both the fluid normal and axial velocity components during plate separation. It is interesting to note from the lower solution branch that both the fluid normal and axial velocity components increase with an increase in the plates separation (i.e., R<0) whereas the trend is opposite for upper solution branch. Meanwhile, we observe the possibility of flow reversal near the plates with increasing plate separation. Figure 4 displays the fluid normal and axial velocity components during plate squeezing (R>0) for the only solution branch in this region. Both the normal and axial velocity components decrease with an increase in plates squeezing.
Two branches of the 10th order SLM approximate solution for normal velocity profile F(η) and longitudinal velocity profile F′(η) for different values of R.
Comparison between the 10th order SLM approximate solution (circles) and the bvp4c numerical results for the normal velocity profile F(η) and longitudinal velocity profile F′(η) for different values of R.
5. Conclusion
In this work, we employed a very powerful new linearisation technique, known as the successive linearisation method (SLM), to solve a fourth order nonlinear differential equation modelling the transient flow of an incompressible viscous fluid between two parallel plates produced by a simple wall motion. The SLM results for the governing flow parameters were compared with results obtained using MATLAB’s bvp4c function and excellent agreement was observed. Using the SLM, it was also shown that the governing problem admits multiple solutions when R<0. The ability of the SLM to generate multiple solutions makes it superior to most numerical methods which are only capable of generating one solution of nonlinear equations. Another significant advantage of the SLM is that its implementation does not depend on small parameters unlike other traditional perturbation methods. The study confirms that the proposed SLM approach converges rapidly to the solution of the original nonlinear problem and can be used to solve many other nonlinear equations arising in fluid mechanics and nonlinear science in general.
Conflict of Interests
The authors declare that there is no conflict of interests.
BertramC. D.Unstable equilibrium behaviour in collapsible tubes198619161692-s2.0-002257902610.1016/0021-9290(86)90109-0JacksonJ. D.A study of squeezing flow19631111481522-s2.0-000136082910.1007/BF03184719BertramC. D.PedleyT. J.A mathematical model of unsteady collapsible tube behaviour198215139502-s2.0-002007725710.1016/0021-9290(82)90033-1BujurkeN. M.AcharP. K.PaiN. P.Computer extended series for squeezing flow between plates1995162-31731872-s2.0-002916366410.1016/0169-5983(94)00058-8MakindeO. D.MotsumiT. G.RamolloM. P.Squeezing flow between parallel plates: a bifurcation study2002928194MR1960522WangC. Y.Squeezing of fluid between two plates197643457958310.1115/1.3423935ZBL0345.76017IshizawaS.The unsteady flow between two parallel discs with arbitrary varying gap width196693553355010.1299/jsme1958.9.533MakindeO. D.Fluid dynamics of parallel plates viscometer: a case study of methods of series summation200326440541710.2989/16073600309486071MR2046146ZBL1152.76358AbbasbandyS.HayatT.Solution of the MHD Falkner-Skan flow by homotopy analysis method2009149-103591359810.1016/j.cnsns.2009.01.030MR2509950ZBL1221.76133LiaoS.20032Chapman & Hall/CRC PressModern Mechanics and MathematicsMR2058313LiaoS.-J.A uniformly valid analytic solution of two-dimensional viscous flow over a semi-infinite flat plate199938510112810.1017/S0022112099004292MR1690937ZBL0931.76017HeJ.-H.Homotopy perturbation technique19991783-425726210.1016/S0045-7825(99)00018-3MR1711041ZBL0956.70017HeJ.-H.A coupling method of a homotopy technique and a perturbation technique for non-linear problems2000351374310.1016/S0020-7462(98)00085-7MR1723761MotsaS. S.SibandaP.ShateyiS.A new spectral-homotopy analysis method for solving a nonlinear second order BVP20101592293230210.1016/j.cnsns.2009.09.019MR2602713ZBL1222.65090MotsaS. S.SibandaP.AwadF. G.ShateyiS.A new spectral-homotopy analysis method for the MHD Jeffery-Hamel problem20103971219122510.1016/j.compfluid.2010.03.004MR2645383ZBL1242.76363HeJ.-H.Variational iteration method—a kind of non-linear analytical technique: some examples19993446997082-s2.0-000009267310.1016/S0020-7462(98)00048-1ZBL05137891MotsaS. S.ShateyiS.A successive linearization method approach to solve Lane-Emden type of equations201220121428070210.1155/2012/280702MR2992435ZBL1264.65104MakukulaZ. G.SibandaP.MotsaS. S.ShateyiS.On new numerical techniques for the MHD flow past a shrinking sheet with heat and mass transfer in the presence of a chemical reaction201120111948921710.1155/2011/489217MR2846134ZBL1235.76108MotsaS. S.ShateyiS.Successive linearisation analysis of unsteady heat and mass transfer from a stretching surface embedded in a Porous medium with suction/injection and thermal radiation effects20119051323133510.1002/cjce.20640CanutoC.HussainiM. Y.QuarteroniA.ZangT. A.1988Berlin, GermanySpringerMR917480TrefethenL. N.200010Philadelphia, Pa, USASIAMSoftware, Environments, and Tools10.1137/1.9780898719598MR1776072KierzenkaJ.ShampineL. F.A BVP solver based on residual control and the MATLAB PSE200127329931610.1145/502800.502801MR1918120ShampineL. F.GladwellI.ThompsonS.2003Cambridge, UKCambridge University PressMR1985643