Application of Successive Linearisation Method to Squeezing Flow with Bifurcation

This paper employs the computational approach known as successive linearizationmethod (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.


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 crosssectional 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][4][5][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][10][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][18][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.the plates.Take a Cartesian coordinate system (, ) where  lies in the streamwise direction and  is the distance measured in the transverse direction.Let  and V be the velocity components in the directions of  and  increasing, respectively.It is assumed that the two plates are at  = ± 0 √(1 − ), where  0 is the position at time  = 0 as shown in Figure 1.

Mathematical Formulation
When  is positive, the two plates are squeezed symmetrically until they touch at  = 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][7][8], the two-dimensional governing equation of motion in terms of vorticity () and stream function (Ψ) formulation is given as with We introduce the following transformations: ( Substituting ( 3) into ( 1) and (2), we obtain where  =  2 0 /2] is the local Reynolds number ( > 0 represents squeezing and  < 0 represents separation).The wall skin friction is given by 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 and we obtain 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.

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 () can be expanded as where   are unknown functions.The solutions of   , ( = 1, 2, . ..) are obtained recursively by solving the linear part of the equation that results from substituting (10) in the governing equations (4) using  0 () as an initial approximation.The linearisation technique is based on the assumption that   becomes increasingly smaller as  becomes large; that is, lim The initial approximation  0 () must be chosen in such a way that it satisfies the boundary conditions ( 5) and ( 6).An appropriate initial guess is where  is an arbitrary constant which when varied results in multiple solutions.Substituting (10) in the governing equations and neglecting nonlinear terms in   ,    ,    , and    give where Starting from the initial approximation,  0 , the subsequent solutions for   ,  ≥ 1, are obtained iteratively by solving (13) subject the the boundary conditions Once each solution for   ( ≥ 1) has been obtained, the approximate solutions for () are obtained as where  is the order of SLM approximation.It is worth noting that the coefficient parameters and the right hand side of (13) for  = 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 where  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 The derivative of   at the collocation points is represented as where D = ((2/)D)  and D is the Chebyshev spectral differentiation matrix (see, e.g., [20,21]).Substituting ( 17)-( 19) in ( 13) results in the matrix equation In the above definitions, a  ( = 1, 2, 3, 4) are diagonal matrices of size ( + 1) × ( + 1).After modifying the matrix system (20) to incorporate boundary conditions, the solution is obtained as

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   (0) and   (1) at selected values of .By using different values of the constant  in the initial approximation (see (12)), it was found that both the SLM and bvp4c give multiple solutions when  < 0. Two solutions, called lower branch and upper branch, were identified when −3.45 ≤  ≤ −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   (1) decreases when  is increased in the lower branch and the opposite effect is observed when  is increased in the upper branch.
Figure 2 illustrates a slice of bifurcation diagram in both planes.For symmetrical squeezing of the plates; that is,  > 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 ( < 0); this is the upper solution branch.A turning point exists between the primary and secondary solution branches at  = −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.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.,  < 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 ( > 0) for the only solution branch in this region.Both the normal and axial velocity components decrease with an increase in plates squeezing.

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

Figure 4 :
Figure 4: Comparison between the 10th order SLM approximate solution (circles) and the bvp4c numerical results for the normal velocity profile () and longitudinal velocity profile   () for different values of .

Table 1 :
The tenth order SLM approximation for   (0) and   (1) at the two branches of solutions for different values of .