The Solution of a Coupled Nonlinear System Arising in a Three-Dimensional Rotating Flow Using Spline Method

The behavior of the non-linear-coupled systems arising in axially symmetric hydromagnetics flow between two horizontal plates in a rotating system is analyzed, where the lower is a stretching sheet and upper is a porous solid plate. The equations of conservation of mass and momentum are transformed to a system of coupled nonlinear ordinary differential equations. These equations for the velocity field are solved numerically by using quintic spline collocation method. To solve the nonlinear equation, quasilinearization technique has been used. The numerical results are presented through graphs, in which the effects of viscosity, through flow, magnetic flux, and rotational velocity on velocity field are discussed.


Introduction
In fluid mechanics, the problems associated with the flow that occurs due to the rotation of a single disk or that between two rotating disks have been found of interest of many researchers.Flows between finite disks were studied by Dijkstra and van Heijst 1 , Adams and Szeri 2 and Szeri et al. 3 .Berker 4 showed that when the two disks are rotating with the same angular speed, there exists a one parameter family of solutions all but one of which is not rotationally symmetric.This result has been extended by Parter and Rajagopal 5 , to disks rotating with differing angular speeds; they prove that the rotationally symmetric solutions are never isolated when considered within the full scope of the Navier-Stokes equations.The numerical study of the asymmetric flow has been carried out by Lai et al. 6,7 .Recently, hydromagnetic flow and heat transfer problems have become more important industrially.In view of these, Chakrabarti and Gupta 8 studied the hydromagnetic flow and heat transfer in a fluid, initially at rest and at uniform temperature, over a stretching International Journal of Mathematics and Mathematical Sciences sheet at a different uniform temperature.Banerjee 9 studied the effect of rotation on the hydromagnetic flow between two parallel plates where the upper plate is porous and solid, and the lower plate is a stretching sheet.
In this paper, we analyze the behavior of the solution of the nonlinear coupled systems arising in axially symmetric hydromagnetic flow between two horizontal plates in a rotating system, where the lower plate is a stretching sheet.The governing coupled ordinary differential equations are solved by quintic spline collocation method.
In Section 2, the mathematical model of the problem given by Vajravelu and Kumar 10 is presented.The quintic spline collocation method is explained in Section 3. The results are displayed in graphical manner in Section 4, and the discussion of results is drawn in Section 5.

Formulation of the Problem
We consider the steady flow of an electrically conducting fluid between two horizontal parallel plates when the fluid and the plates rotate in unison about an axis normal to the plates with an angular velocity Ω.A Cartesian coordinate system is considered in such a way that the x-axis is along the plate, the y-axis is perpendicular to it, and the z-axis is normal to the xy-plane as shown in Figure 1.
The origin is located at the centre of the channel, and the plates are located at y −h and h.The lower plate is is stretched by introducing two equal and opposite forces so that the position of the point 0, −h, 0 remains unchanged.A uniform magnetic flux with density B o is acting along y-axis about which the system is rotating.The upper plate is subjected to a constant wall injection with a velocity V o .The equations of motion in a rotating frame of reference are where u, υ, and ω denote the fluid velocity components along the x, y, and z directions, ν is the kinematics coefficient of viscosity, ρ is the fluid density, and p * is the modified fluid pressure.The velocity is independent of z and that all derivatives towards z therefore do not appear in the equations of motion.This leads to the absence of ∂p * /∂z in 2.3 which implies that there is a net cross-flow along the z-axis.

2.5
where a prime denotes differentiation with respect to η. Substituting 2.6 in 2.1 to 2.4 , we have where , the rotation parameter.

2.10
Equation 2.7 with the help of 2.8 can be written as where A is a constant.

International Journal of Mathematics and Mathematical Sciences
Differentiation of 2.11 with respect to η gives Thus we solve the following nonlinear system numerically for several sets of values of the parameters: subject to the boundary conditions where λ υ 0 /Eh a parameter depending on the y component of velocity at the upper plate.

Quintic Spline Collocation Method
The fifth degree spline is used to find numerical solutions to the boundary value problems discussed in 2.13 together with 2.14 .A detailed description of spline functions generated by subdivision is given by de Boor 11 .
Consider equally spaced knots of a partition π: Let S 5 π be the space of continuously differentiable, piecewise, Quintic polynomials on π.That is, S 5 π is the space of Quintic polynomials on π.The Quintic spline is given by Bickley 12 and by G. Micula and S. Micula 13 where the power function x − x k is defined as Consider a fourth-order linear boundary value problem of the form y iv x p x y x q x y x r x y x t x y x m x , a ≤ x ≤ b, 3.3 subject to the boundary conditions where y x , p x , q x , r x , t x , and m x are continuous functions defined in the interval x ∈ a, b ; η 0 , η 1 , η 2 , η 3 are finite real constants.Let 3.1 be an approximate solution of 3.3 , where a 0 , b 0 , c 0 , d 0 , e 0 , F 0 , F 1 , . . ., F n−1 are real coefficients to be determined.
Let x 0 , x 1 , . . ., x n be n 1 grid points in the interval a, b , so that It is required that the approximate solution 3.1 satisfies the differential equation at the point x x i .Putting 3.1 with its successive derivatives in 3.3 , we obtain the collocation equations as follows: From boundary conditions,

3.7
Using the power function x − x k in the above equations, a system of n 5 linear equations in n 5 unknowns a 0 , b 0 , c 0 , d 0 , e 0 , F 0 , F 1 , . . ., F n−1 is thus obtained.This system can be written in matrix-vector form as follows: where The coefficient matrix A is an upper triangular Hessenberg matrix with a single lower subdiagonal, principal and upper diagonal having nonzero elements.Because of this nature of matrix A, the determination of the required quantities becomes simple and consumes less time.The values of these constants ultimately yield the quintic spline s x in 3.1 .
In case of nonlinear boundary value problem, the equations can be converted into linear form using quasilinearization method Bellman and Kalaba 14 , and hence this method can be used as iterative method.The procedure to obtain a spline approximation of y i i 0, 1, 2, . . ., j, where j denotes the number of iteration by an iterative method starts with fitting a curve satisfying the end conditions and this curve is designated as y i .We obtain the successive iterations y i 's with the help of an algorithm described as above till desired accuracy.

Quintic Spline Solution
We solve following nonlinear system numerically for several sets of values of the parameters: subject to the boundary conditions f λ, f 0, g 0 at η 1.

4.3
The spline collocation method is used to solve the differentiation system 4.1 to 4.3 .Equation 4.2 is a linear equation of order two, whereas 4.1 is a nonlinear equation of order four.For solving nonlinear equation by spline collocation method, we require a linear form of differentiation equation.The quasilinearization technique transforms 4.1 into linearized form as with boundary conditions

4.5
The linear equation 4.2 can be written as with boundary conditions We use the quintic spline method to find solutions of 4.4 of fourth order with the conditions 4.5 , whereas the cubic spline method is used for finding solutions of second-order equation 4.6 and 4.7 .

International Journal of Mathematics and Mathematical Sciences
The quintic spline given by is an approximate solution of the problem given by 4.4 and 4.5 .Substituting s η in both of above equations, we obtain the collocation as First two boundary conditions in 4.5 immediately give and other two boundary conditions are 1 120

4.11
International Journal of Mathematics and Mathematical Sciences 9 First of all, initial assumptions regarding f i , f i , f i , f i are necessary.Let a curve aη 3 bη 2 cη d be fitted through the points η −1 and η 1 satisfying the boundary conditions 4.5 .This requirement yields a In 4.6 and 4.7 , the following set of equations is obtained: International Journal of Mathematics and Mathematical Sciences and the last equation is A straight line g η aη b can be fitted to the boundary conditions for g.A line g η 0 satisfied the boundary conditions 4.7 .Therefore, initial values of g and g are also zero.
For N 10, R 0.1 and several sets of values of the parameters λ, M 2 and K 2 , the velocity components f, f , and g are obtained in graphical form.The behavior of velocity components are shown in Figures 2, 3, and 4. The results obtained by the spline collocation method are compared with the graphical solutions obtained by Vajravelu and Kumar 10 , which are tested by comparing the analytical solutions for small value of R with the exact solutions f and g.The numerical results thus obtained by the spline method to find f, f , and g are in close agreement with the available data.In Figures 2, 3, and 4, η equals a dimensionless injection parameter, M 2 is a magnetic parameter, and K 2 is a rotational parameter.
With improper initial guess, convergence is not observed every time.As discussed by Vajravelu and Kumar 10 , this is due to inherent instability of boundary-value problem.Because of this, integration from the starting point of the domain may produce rapidly increasing solutions, which may occasionally lead to overflow before the end-point is reached.

Discussion of the Results
The behavior of the nondimensional velocity component f for R 0.1 is shown in Figure 2.
From Figure 2, it is evident that the Lorenz force decreases the velocity component f see curves II and III , but the value of f increases with an increase in the value of the parameter λ as seen in curves III and IV.Further, in Figure 2, curves IV and V show that f decreases steadily for small K 2 , whereas for large value of K 2 , f increases near the lower plate and decreases near the upper plate.Figure 3 describes the behavior of f η , which is proportional to velocity along parallel plates, for several sets of values of the parameters only λ, M 2 and K 2 the rotation parameter are varied while R is kept constant at 0.1.From Figure 3, curves III and IV, it is evident that the velocity component f increases with the parameter λ υ 0 /Eh, y component of velocity at the upper plate , and the increase is maximized near the stretching plate.However, the effect of the Lorenz force magnetic parameter M 2 of f is to increase it near the upper plate and to decrease it near the lower stretching plate that is displayed from the curves II and III of Figure 3. Further, the effect of small K 2 is to increase f at the upper plate and to decrease it near the lower plate is evident from curves I and II.But the effect of large K 2 is to increase f near the plates and to decrease it at the centre of the channel.Also, the rotation of the channel brings humps near the plates, indicating the occurrence of a boundary layer near the plates.Furthermore, for large K 2 the coriolis force and the magnetic field that act against the pressure gradient cause reversal of the flow.

International Journal of Mathematics and Mathematical Sciences
The transverse velocity g increases as the parameter λ increases, which is shown in curves III and IV of Figure 4.But this is quite opposite to the phenomenon with the magnetic parameter M 2 .However, the rotation parameter K 2 increases g and the maximum of g occurs near the stretching sheet for large K 2 .
From above figures, it is observed that for large K 2 convergence in the results is not guaranteed.This fact is also observed by Vajravelu and Kumar 10 .
We conclude from the above numerical experience that the spline collocation method can be successfully applied to solve the fourth-order nonlinear coupled equation, which governs the hydromagnetic fluid flow.This widens the field of applicability to the higherorder coupled differential equations.This encourages the applications to other types of problems.
Another useful conclusion is that the selection of the domain is not restricted to positive interval only.That is we can successfully apply the above method for the negative intervals as the domain of the problem.