Research Article  Open Access
S. Abdul Gaffar, V. Ramachandra Prasad, Bhuvana Vijaya, O. Anwar Beg, "Mixed Convection Flow of Magnetic Viscoelastic Polymer from a Nonisothermal Wedge with Biot Number Effects", International Journal of Engineering Mathematics, vol. 2015, Article ID 287623, 15 pages, 2015. https://doi.org/10.1155/2015/287623
Mixed Convection Flow of Magnetic Viscoelastic Polymer from a Nonisothermal Wedge with Biot Number Effects
Abstract
Magnetic polymers are finding increasing applications in diverse fields of chemical and mechanical engineering. In this paper, we investigate the nonlinear steady boundary layer flow and heat transfer of such fluids from a nonisothermal wedge. The incompressible EyringPowell nonNewtonian fluid model is employed and a magnetohydrodynamic body force is included in the simulation. The transformed conservation equations are solved numerically subject to physically appropriate boundary conditions using a secondorder accurate implicit finite difference Keller Box technique. The numerical code is validated with previous studies. The influence of a number of emerging nondimensional parameters, namely, the EyringPowell rheological fluid parameter (), local nonNewtonian parameter based on length scale (), Prandtl number (Pr), Biot number (), pressure gradient parameter (), magnetic parameter (), mixed convection parameter (), and dimensionless tangential coordinate (), on velocity and temperature evolution in the boundary layer regime is examined in detail. Furthermore, the effects of these parameters on surface heat transfer rate and local skin friction are also investigated.
1. Introduction
The development of modern functional materials which can be manipulated using electromagnetic fields has stimulated great attention in polymer engineering sciences in recent years. These materials include ferromagnetic polymers [1], magnetotropic fluids [2], and liquid crystalline electrically conducting polymers [3]. Such fluids exhibit many complex characteristics including nonNewtonian behavior. For better understanding of the manufacture of such materials, magnetorheological fluid mechanics plays a central role. Although many different magnetic properties may arise in such materials, a fundamental methodology for their synthesis involves the application of transverse magnetic fields which alters momentum transfer and therefore also influences coupled processes such as heat transfer and mass transfer (species diffusion). Many constitutive material models have been developed to simulate the departure of such fluids from Newtonian viscous behavior. In these fluids, the constitutive relationship between stress and rate of strain is nonlinear in comparison to the NavierStokes equations which are generally good for Newtonian fluids. Most nonNewtonian models involve some form of modification to the momentum conservation equations. A comprehensive summary of such models which include the Maxwell model, powerlaw model, and also EyringPowell and Giesekus viscoelastic models, is provided in the treatise by Shaw [4]. Heat transfer is often also frequently present in polymer dynamics [5]. Furthermore, in many polymer fabrication flows, boundary layer phenomena arise. The convective boundary condition has therefore also attracted some interest and this is usually simulated via a Biot number in the wall thermal boundary condition. Recently, Ishak [6] discussed the similarity solutions for flow and heat transfer over a permeable surface with convective boundary condition. Aziz [7] provided a similarity solution for laminar thermal boundary layer over a flat surface with a convective surface boundary condition. Aziz [8] further studied hydrodynamic and thermal slip flow boundary layers with an isoflux thermal boundary condition. Buoyancy effects on thermal boundary layers subject to a convective surface boundary condition were examined by Makinde and Olanrewaju [9]. Gupta et al. [10] used a variational finite element to simulate mixed convectiveradiative nonNewtonian shrinking sheet flow with a convective boundary condition and Eringen’s micropolar material model. Swapna et al. [11] studied convective wall heating effects on hydromagnetic flow of a micropolar fluid. Makinde et al. [12] studied cross diffusion effects and Biot number influence on hydromagnetic Newtonian boundary layer flow with homogenous chemical reactions and MAPLE quadrature routines. Bég et al. [13] analyzed Biot number and buoyancy effects on magnetohydrodynamic thermal slip flows. Subhashini et al. [14] studied wall transpiration and cross diffusion effects on free convection boundary layers with a convective boundary condition. In many of these studies, magnetohydrodynamic effects were also considered generally via the introduction of a Lorentzian magnetic drag force and consideration of the Hartmann number. Relatively few studies however have considered magnetic field effects on EyringPowell viscoelastic polymer flows. This rheological model has certain advantages over the other nonNewtonian formulations, including simplicity, ease of computation, and physical robustness. Furthermore, it is deduced from kinetic theory of liquids rather than the empirical relation. Additionally, it correctly reduces to Newtonian behavior for low and high shear rates [15]. It is a threeconstant model which displays a nonzero bounded viscosity at both the upper and the lower limits. Distinct from the popular OstwaldDe Waele powerlaw model, the EyringPowell model does not demonstrate infinite effective viscosities for low shear rates. Investigations employing this model of relevance to polymeric transport processes include Hayat et al. [16] in radiative magnetic convection, Sirohi et al. [17] for wedge flows, and Adesanya and Gbadeyan [18] for channel flows. Very recently, Prasad et al. [19] studied the thermal convection flow in permeable materials with Biot number effects, observing that increasing rheological effect accelerates the flow and heats the boundary layer.
Very few of the above studies have considered FalknerSkan flows [20]. This family of boundary layer flows is associated with the twodimensional wedge configuration. NonNewtonian flows from wedge bodies arise in a number of chemical engineering systems which have been described in detail by Peddieson [21] employing the secondorder ReinerRivlin model. The mixed convection boundary layer flow from a heated wedge plate has also drawn some interest. The combined forced and free convection flow and heat transfer about a nonisothermal wedge subject to a nonuniform free stream velocity were first considered by Sparrow et al. [22]. Watanabe et al. [23] analyzed theoretically mixed convection flow over a perforated wedge with uniform suction or injection. Kafoussias and Nanousis [24] and Nanousis [25] studied the effect of suction or injection on MHD mixed convection flow past a wedge. Gorla [26] used a powerlaw model to study heat transfer in polymer flow past a wedge. Yih [27] evaluated radiation effects on mixed convection flow about an isothermal wedge embedded in a saturated porous medium. Rashidi et al. [28] developed homotopy solutions for third grade viscoelastic flow from a nonisothermal wedge. Chamkha et al. [29] presented computational solutions for MHD forced convection flow from a nonisothermal wedge in the presence of a heat source or sink with a finite difference method. Hsiao [30] reported on MHD convection of viscoelastic fluid past a porous wedge, observing that the elastic effect increases the local heat transfer coefficient and heat transfer rates at the wedge surface. Ishak et al. [31] obtained a selfsimilar solution for a moving wedge in a micropolar fluid. Ishak et al. [32] further studied numerically steady twodimensional laminar flow past a moving wedge in nonNewtonian fluid.
The objective of the present study is to investigate the laminar boundary layer flow and heat transfer of an EyringPowell nonNewtonian fluid from a nonisothermal wedge. Such a study has not appeared in the literature to the knowledge of the authors. The nondimensional equations with associated dimensionless boundary conditions are solved with the Keller implicit finite difference “Box” scheme [33]. The effects of the emerging thermophysical parameters, namely, the rheological parameters , Biot number , mixed convection parameter , pressure gradient parameter , magnetic parameter , and Prandtl number , on the velocity, temperature, local skin friction, and heat transfer rate (local Nusselt number) characteristics are studied. The present problem is relevant to the simulation of magnetized polymer materials fabrication processes.
2. NonNewtonian Constitutive EyringPowell Fluid Model
In the present study, a subclass of nonNewtonian fluids known as the EyringPowell fluid is employed. The Cauchy stress tensor for EyringPowell fluids [15] takes the hyperbolic form: where is dynamic viscosity and and are the rheological fluid parameters of the EyringPowell fluid model. Consider the secondorder approximation of the function as follows:When the EyringPowell formulation is introduced into the subsequent model, the momentum conservation equation is significantly modified with numerous velocity gradient terms including mixed derivatives. Strong nonlinearity therefore results in necessitating numerical solutions. The resulting boundary value problem is found to be well posed and permits an excellent mechanism for the assessment of rheological characteristics on the flow behaviour.
3. Mathematical Thermofluid Polymer Boundary Layer Model
Steady, laminar, and incompressible flow and heat transfer of an EyringPowell polymeric fluid from a nonisothermal horizontal wedge are considered, as illustrated in Figure 1.
The coordinate (streamwise) is measured along the wall of the wedge and coordinate (transverse) is directed normal to it. A uniform magnetic field of strength is applied parallel to the axis, that is, transverse to the wedge surface. Magnetic induction, Hall current, ion slip, and Maxwell displacement currents are neglected. Magnetic Reynolds number is small and the electric field due to polarization of charges is negligible. The gravitational acceleration acts downwards. We also assume that the Boussinesq approximation holds; that is, density variation is only experienced in the buoyancy term in the momentum equation. Initially, both the wedge and EyringPowell fluid are maintained at the same temperature. Instantaneously, they are raised to a temperature , the ambient temperature of the fluid which remains unchanged. In line with the approach of Yih [27] and introducing the boundary layer approximations, the equations for mass, momentum, and energy conservation can be written as follows: Here, and are the velocity components in  and directions, respectively, is the kinematic viscosity of EyringPowell fluid (polymer), is the coefficient of thermal expansion, is the thermal diffusivity, is the temperature, and is the density of the fluid. The EyringPowell fluid model therefore introduces a mixed derivative (second order, second degree) into the momentum boundary layer equation (4). The nonNewtonian effects feature in the shear terms only of (4) and not the convective (acceleration) terms. The fourth term on the right hand side of (4) represents the thermal buoyancy force on the flow field, with “+” and “−” signs referring to the buoyancyassisting and buoyancyopposing flow scenarios, and couples the velocity field (4) with the temperature field equation (5). The physical boundary conditions are given as follows:Here, is the free stream temperature, is the thermal conductivity, is the convective heat transfer coefficient, is the convective fluid temperature, is the free stream velocity, and is the Hartree pressure gradient parameter which corresponds to for a total angle Ω of the wedge. is a positive number. The stream function is defined by and , and, therefore, the continuity equation is automatically satisfied. In order to render the governing equations and the boundary conditions in dimensionless form, the following nondimensional quantities are introduced:All terms are defined in Nomenclature. In view of the transformation defined in (7), the boundary layer equations (4)–(6) are reduced to the following coupled, nonlinear, and dimensionless partial differential equations for momentum and energy for the domain:The transformed dimensionless boundary conditions are as follows:Here, primes denote the differentiation with respect to and is the Biot number. The wall thermal boundary condition in (10) corresponds to convective cooling. The skin friction coefficient (shear stress at the cylinder surface) and Nusselt number (heat transfer rate) can be defined using the transformations described above with the following expressions:The location, , corresponds to the vicinity of the lower stagnation point on the wedge. For this scenario, the model defined by (8) to (9) contracts to an ordinary differential boundary value problem: The general model is solved using a powerful and unconditionally stable finite difference technique introduced by Keller [33]. The Keller Box method has a secondorder accuracy with arbitrary spacing and attractive extrapolation features. It remains one of the most reliable numerical methods for nonlinear boundary layer flows.
4. Numerical Solution with Keller Box Implicit Method
The Keller Box implicit difference method is implemented to solve the nonlinear boundary value problem defined by (8) and (9) with boundary conditions (10). This technique, despite recent developments in other numerical methods, remains a powerful and very accurate approach for boundary layer flow equation systems which are generally parabolic in nature. It is unconditionally stable and achieves exceptional accuracy. An excellent summary of this technique is given in Keller [33]. Magnetohydrodynamics applications of Keller’s method are reviewed in Anwar Bég [34]. This method has also been applied successfully in many rheological flow problems in recent years. These include oblique micropolar stagnation flows [35], Walter’s B viscoelastic flows [36], Stokesian couple stress flows [37], hyperbolictangent convection flows from curved bodies [38], micropolar nanofluids [39], Jeffreys elastoviscous boundary layers [40], magnetic Williamson fluids [41], and Maxwell fluids [42]. The Keller Box discretization is fully coupled at each step which reflects the physics of parabolic systems, which are also fully coupled. Discrete calculus associated with the Keller Box scheme has also been shown to be fundamentally different from all other mimetic (physics capturing) numerical methods, as elaborated in Anwar Bég [34]. The Keller Box scheme comprises four stages:(1)Decomposition of the th order partial differential equation system to firstorder equations.(2)Finite Difference Discretization.(3)Quasilinearization of Nonlinear Keller Algebraic Equations.(4)Blocktridiagonal elimination solution of the Linearized Keller Algebraic Equations.The algebraic details for the present problem are provided in Appendix A. A typical mesh is given in Figure 2.
5. Numerical Results and Interpretation
Comprehensive solutions have been obtained and are presented in Figures 3–11 and Tables 1–5. The numerical problem comprises two independent variables , two dependent fluid dynamic variables , and eight thermophysical and body force control parameters, namely, , , , , , , , and . The following default parameter values, that is, , , , , , , , and , are prescribed (unless otherwise stated). Furthermore, the influence of streamwise (transverse) coordinate on heat transfer characteristics is also investigated.





(a)
(b)
(a)
(b)
(a)
(b)
(a)
(b)
(a)
(b)
(a)
(b)
(a)
(b)
(a)
(b)
(a)
(b)
6. Validation with Nakamura Tridiagonal Method (NTM)
To verify the accuracy of the Keller Box method (KBM) code, we have also solved the seventhorder nonlinear partial differential boundary value problem defined by the parabolic equations (8) and (9) with boundary conditions (10), with an alternative finite difference procedure, namely, the Nakamura tridiagonal method (NTM) [43]. NTM is also second order, accurate, and very easily implemented in boundary layer problems. It has proven to be particularly adept at resolving nonlinear and mixed derivatives which characterize nonNewtonian (e.g., viscoelastic) problems. Applications of NTM in rheological boundary layer flows include secondorder ReinerRivlin wedge convection in porous media [44], powerlaw convective flows in porous media [45], micropolar wall plume heat transfer [46], and nanofluid bioconvection in permeable materials [47] and magnetohydrodynamic micropolar enrobing flows [48]. Details of the technique with solutions are provided in Appendix B and Tables 6 and 7.


7. Discussion and Interpretation
In Figures 3–11, we note that asymptotically smooth profiles are always achieved in the free stream. An adequate value of infinity has therefore been prescribed in the computations. Figures 3(a) and 3(b) illustrate the effect of EyringPowell fluid parameter on the velocity and temperature distributions through the boundary layer regime. Velocity is significantly decreased with increasing at larger distance from the wedge surface. Conversely, temperature is consistently enhanced with increasing values of . The mathematical model reduces to the Newtonian viscous flow model as and . The momentum boundary layer equation in this case contracts to the familiar equation for Newtonian mixed convection from a plate; namely, . The thermal boundary layer equation (9) remains unchanged. In Figure 3(b), temperatures are clearly minimized for the Newtonian case () and maximized for the strongest nonNewtonian case . The nonNewtonian parameter (= ) is evidently inversely proportional to the viscosity and also to and (EyringPowell rheological fluid parameters). As is increased, the material viscosity therefore will decrease which would aid momentum development. Similarly, the parameters and are simultaneously reduced (they are strongly linked to viscosity) and the overall effect is a boost in material viscosity. This results in reduction in the magnitude which manifests in marked deceleration in the flow as shown in Figure 3(a). With decreased momentum diffusion rates, the velocity boundary layer thickness is decreased. Thermal diffusion is conversely enhanced and the fluid is heated (Figure 3(b)), resulting in a thickening in the thermal boundary layer.
Figures 4(a) and 4(b) depict the velocity and temperature distributions with increasing local nonNewtonian parameter, . Very little tangible effect is observed in Figure 4(a), although there is a very slight increase in velocity with increase in . Similarly, there is only a very slight depression in temperature magnitudes in Figure 4(b) with a rise in . This parameter is defined as and features in a single negative term in (9), namely, , unlike the rheological parameter which also appears in the shear term, . This parameter also features viscosity which is effectively decreased weakly and induces a slight acceleration in the flow. The momentum boundary layer thickness is therefore slightly decreased whereas since heat diffusion is decreased, thermal boundary layer thickness is therefore marginally reduced.
Figures 5(a) and 5(b) depict the evolution of velocity and temperature functions with a variation in Biot number, . Dimensionless velocity component (Figure 5(a)) is considerably enhanced with increasing . In Figure 5(b), an increase in Biot number is seen to considerably enhance temperatures throughout the boundary layer regime. For , that is, small Biot numbers, the regime is frequently designated as being “thermally simple” and there is a presence of more uniform temperature fields inside the boundary layer and the wedge solid surface. For , thermal fields are anticipated to be nonuniform within the solid body. The Biot number effectively furnishes a mechanism for comparing the conduction resistance within a solid body to the convection resistance external to that body (offered by the surrounding fluid) for heat transfer. We also note that a Biot number in excess of 0.1, as studied in Figures 5(a) and 5(b), corresponds to a “thermally thick” substance whereas Biot number less than 0.1 implies a “thermally thin” material. Since is inversely proportional to thermal conductivity (), as increases, thermal conductivity will be reduced at the cylinder surface and this will lead to a decrease in the rate of heat transfer from the boundary layer to within the cylinder, manifesting in a rise in temperature at the cylinder surface and in the body of the fluid; the maximum effect will be sustained at the surface, as witnessed in Figure 5(b). For a fixed wall convection coefficient and thermal conductivity, Biot number as defined in is also inversely proportional to the local Reynolds (free convection) number. However, as the Biot number increases, the local Reynolds number must decrease and this will induce the opposite effect, that is, accelerating the boundary layer flow, as shown in Figure 5(a).
Figures 6(a) and 6(b) depict the profiles for velocity and temperature for various values of pressure gradient parameter, . It is observed that increasing significantly accelerates the flow; that is, velocity increases. This trend is familiar from other studies, for example, Nanousis [25] and Gorla [26]. However, increasing is found to decelerate the temperature. Effectively, with greater pressure gradient effect, the momentum boundary layer thickness is reduced and thermal boundary layer thickness is decreased. In (9), with , the wedge flow becomes Blasius flow from a flat plate. With , we retrieve the case of flow in the vicinity of stagnation point flow on an infinite plate. While these two flow cases are of mathematical interest, in practical polymer processing, they are not. The generalized wedge case is more relevant for which , as studied in Figures 6(a) and 6(b).
Figures 7(a) and 7(b) depict the profiles for velocity and temperature for various values of mixed convection parameter, . It is found that an increase in increases the velocity. However, the temperature decreases significantly. The mixed convection parameter couples the momentum field with the energy field, via the buoyancy term, . Although this is a firstorder term, it has a profound influence on the interplay of velocity and temperature. With positive values, we have buoyancyaided flow and with negative values we have buoyancyopposed flow. Both are studied in Figures 7(a) and 7(b). Evidently when buoyancy aids the flow , velocities are increased, with the contrary effect for buoyancy opposition. The momentum boost achieved with assistive buoyancy force results however in a depletion of thermal energy which cools the boundary layer and decreases temperatures (Figure 7(b)). The reverse effect is generated with opposing buoyancy force which heats the boundary layer and increases thermal boundary layer thickness.
Figures 8(a) and 8(b) depict the profiles for velocity and temperature for various values of Prandtl number, . It is observed that an increase in the Prandtl number significantly decelerates the flow; that is, velocity decreases. Also increasing Prandtl number is found to decelerate the temperature. Prandtl numbers of 100 are more representative of weak polymer solutions and are associated with slower and cooler flows. Momentum boundary layer thickness is greater for such fluids as is thermal boundary layer thickness.
Figures 9(a) and 9(b) depict the profiles for velocity and temperature for various values of magnetic parameter, . It is observed that an increase in significantly decelerates the flow; that is, velocity decreases. Conversely, increasing is found to enhance the temperature. The parameter, , is a Hartmann number. It simulates the relative contribution of Lorentzian magnetohydrodynamic drag force relative to viscous hydrodynamic force. As is increased, greater opposition is generated to the flow past the wedge leading to deceleration. The supplementary work expended in dragging the polymer against the imposition of the transverse magnetic field creates heating in the polymer. This dissipation of heat leads to a temperature rise and thickening of thermal boundary layers in polymers. Such phenomena are documented extensively in magnetohydrodynamic studies, for example, in Makinde et al. [12] and Chamkha et al. [29].
Figures 10(a) and 10(b) show the influence of EyringPowell fluid parameter, , on dimensionless skin friction coefficient and heat transfer rate at the wedge surface. It is observed that the dimensionless skin friction is increased with the increase in ; that is, the boundary layer flow is accelerated with decreasing viscosity effects in the nonNewtonian regime. Conversely, the surface heat transfer rate is substantially decreased with increasing values. Decreasing viscosity of the fluid (induced by increasing the value) reduces thermal diffusion as compared with momentum diffusion. A decrease in heat transfer rate at the wall will imply that less heat is convected from the fluid regime to the wedge, thereby heating the boundary layer and enhancing temperatures.
Figures 11(a) and 11(b) illustrate the influence of the local nonNewtonian parameter, , on the dimensionless skin friction coefficient and heat transfer rate . The skin friction (Figure 11(a)) at the wedge surface is accentuated with increasing , however only for very large values of the transverse coordinate, . The flow is therefore strongly accelerated along the wedge surface far from the lower stagnation point. Heat transfer rate (local Nusselt number) is also enhanced with increasing , again at large values of , as computed in Figure 11(b).
In Table 1, we present the influence of the EyringPowell fluid parameter, , on the skin friction and heat transfer rate, along with a variation in and . With increasing , the skin friction is enhanced. The parameter is inversely proportional to the dynamic viscosity of the nonNewtonian fluid. There as is elevated, viscosity will be reduced and this will induce lower resistance to the flow at the surface of the wedge, that is, accelerating the flow leading to an escalation of shear stress. Furthermore, this trend is sustained at any Pr. However, an increase in Pr markedly reduces the shear stress magnitudes. Similarly increasing is observed to reduce heat transfer rates, again at all Pr’s, whereas Pr strongly accentuates heat transfer rates. Magnitudes of shear stress are always positive indicating that flow reversal (backflow) never arises. It is also observed that increasing mixed convection parameter, , increases both skin friction and heat transfer rates.
Table 2 shows results for the influence of the local nonNewtonian parameter (based on length scale ), that is, and the mixed convection parameter () along with the variation in on skin friction and heat transfer rate. Skin friction is generally decreased with increasing . However, heat transfer rate is found to be enhanced with increasing . and inspection of this definition shows that the effectively inverse proportionality of to kinematic viscosity () (with all other parameters being maintained constant) will generate a strong resistance to the flow leading to a deceleration, that is, drop in shear stresses. Conversely, the direct proportionality of to Reynolds number will imply that thermal buoyancy forces are enhanced as increases and this will cause a boost in heat transfer by convection from the wedge surface manifesting with the greater heat transfer rates observed in Table 2. This table also shows that, with an increase in , the skin friction and heat transfer rate are elevated.
Table 3 presents the results for the influence of the pressure gradient parameter () and the magnetic parameter () along with the variation in on skin friction and heat transfer rate. Skin friction and heat transfer rate are generally decreased with increasing . Flow deceleration is therefore confirmed with greater transverse magnetic field, associated with the escalation in Lorentzian magnetohydrodynamic impedance. However, skin friction and heat transfer rate are found to be enhanced with increasing .
Table 4 presents the Keller Box numerical solutions of the missing condition (in brackets) and skin friction for various values of and . It is found that skin friction is reduced with increasing values of . Furthermore, the skin friction is observed to be increased with a rise in the EyringPowell fluid parameter () for all values of the local nonNewtonian parameter ().
Table 5 presents the comparison results of for different values of which are found to be in good correlation with those of [49].
8. Conclusions
Numerical solutions have been presented for the nonsimilar, buoyancydriven boundary layer flow and heat transfer of an EyringPowell fluid external to a nonisothermal wedge. This model has been developed to analyse magnetic polymer FalknerSkan flows of relevance in chemical engineering processing. The Keller Box implicit secondorder accurate finite difference numerical scheme has been utilized to efficiently solve the transformed, dimensionless velocity and thermal boundary layer equations, subject to realistic boundary conditions. A comprehensive assessment of the effects of EyringPowell fluid parameter , local nonNewtonian parameter , Biot number , mixed convection parameter , pressure gradient parameter , magnetic body force parameter , and Prandtl number on thermofluid characteristics has been conducted. Significant flow deceleration has been demonstrated for increasing EyringPowell fluid parameter and negative mixed convection parameter , that is, buoyancy opposition and magnetic body force parameter with generally greater temperatures observed with stronger magnetic field. Excellent correlation with Nakamura tridiagonal method computations has also been achieved, verifying the validity of the present Keller Box solutions. Blasius and stagnation point flow cases have been retrieved as special cases of the generalized wedge flow model. Generally very stable and accurate solutions are obtained with the present finite difference code. The numerical code is able to solve nonlinear boundary layer equations very efficiently and therefore shows excellent promise in simulating transport phenomena in other nonNewtonian fluids. It is therefore presently being employed to study micropolar fluids and viscoplastic fluids which also represent other chemical engineering working fluids. The present study has also neglected time effects. Future simulations will also address transient polymeric boundary layer flows and will be presented soon.
Appendices
A. Keller Box Numerical Details
As mentioned earlier, the algebraic details for the Keller Box method are documented here.
Stage 1: Decomposition of thOrder Partial Differential Equation System to FirstOrder Equations. Equations (8) and (9) subject to the boundary conditions (10) are first cast as a multiple system of firstorder differential equations. New dependent variables are introduced:These denote the variables for velocity, temperature, and concentration, respectively. Now, (8) and (9) are solved as a set of fifthorder simultaneous differential equations:where primes denote differentiation with respect to the variable, . In terms of the dependent variables, the boundary conditions assume the following form:
Stage 2: Finite Difference Discretization. A twodimensional computational grid is imposed on the plane as depicted in Figure 2. The stepping process is defined bywhere is the spacing and is the spacing.
If denotes the value of any variable at , then the variables and derivatives of (A.2)–(A.6) at are replaced byThe finite difference approximation of (A.2)–(A.6) for the midpoint is as follows: Here, we have used the abbreviations:The boundary conditions are as follows:
Stage 3: Quasilinearization of Nonlinear Keller Algebraic Equations. Assuming , , , , to be known for , then (A.10) constitute a system of equations for the solution of unknowns , , , , , . This nonlinear system of algebraic equations is linearized by means of Newton’s method as explained in [33, 34].
Stage 4: BlockTridiagonal Elimination Solution of Linear Keller Algebraic Equations. The linearized system is solved by the blockelimination method, since it possesses a blocktridiagonal structure. The blocktridiagonal structure generated consists of block matrices. The complete linearized system is formulated as a block matrix system, where each element in the coefficient matrix is a matrix itself, and this system is solved using the efficient Keller Box method. The numerical results are strongly influenced by the number of mesh points in both directions. After some trials in the direction (radial coordinate), a larger number of mesh points are selected whereas in the direction (tangential coordinate) significantly less mesh points are utilized. has been set at 12 and this defines an adequately large value at which the prescribed boundary conditions are satisfied. is set at 3.0 for this flow domain. Mesh independence is achieved in the present computations. The numerical algorithm is executed in MATLAB on a PC. The method demonstrates excellent stability, convergence, and consistency, as elaborated by Keller [33].
B. Nakamura Tridiagonal Method (NTM) Verification
For twodimensional partial differential nonsimilar flows, NTM [43–48] executes the computations using both inner and outer loops during the iterative process, the former to advance the solution in the direction and the latter to advance it in the direction. The flow domain is discretized using an equispaced finite difference mesh in the directions. The partial derivatives for , with respect to , are evaluated by central difference approximations. A double iteration loop based on the method of successive substitution is employed. The finite difference discretized equations are solved as a linear secondorder boundary value problem on the domain. The momentum equation (8) is thirdorder whereas the heat equation (9) is a secondorder equation, and for (9), only a direct substitution is needed. Setting Equations (8) and (9) assume the following form.
Nakamura momentum equation is as follows: Nakamura heat (energy) equation is as follows:where , , and are the Nakamura matrix coefficients and are the Nakamura source terms containing a mixture of variables and derivatives associated with the variables (omitted for brevity). The Nakamura equations (B.2) and (B.3) are transformed to finite difference equations and these are orchestrated to form a tridiagonal system which is solved iteratively. Computations are performed in seconds on an SGI Octane dual processor machine. Tables 5 and 6 document the comparison of KBM and NTM solutions for the effects of pressure gradient , Hartmann parameter , streamwise coordinate , and EyringPowell rheological parameters . Excellent agreement is achieved. Confidence in the Keller Box solutions is therefore very high.
Nomenclature
:  Constant magnetic field intensity 
:  Skin friction coefficient 
:  Rheological fluid parameter 
:  Nondimensional stream function 
:  Local Grashof number 
:  Acceleration due to gravity 
:  Thermal conductivity of fluid 
:  Magnetic parameter 
:  Pressure gradient parameter 
:  Local Nusselt number 
:  Prandtl number 
:  Temperature of the fluid 
:  Nondimensional velocity components along the  and directions, respectively 
:  Velocity vector 
:  Streamwise coordinate 
:  Transverse coordinate. 
:  Thermal diffusivity 
:  Fluid parameter 
:  Nondimensional concentration 
:  Dimensionless radial coordinate 
:  Dynamic viscosity 
:  Kinematic viscosity 
:  Nondimensional temperature 
:  Density of nonNewtonian fluid (polymer) 
:  Dimensionless tangential coordinate 
:  Dimensionless stream function 
:  Total angle of the wedge 
:  Mixed convection parameter 
:  Rheological fluid parameter 
:  Biot number. 
:  Conditions at the wall (wedge surface) 
:  Free stream conditions. 
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
References
 G. Bossis, O. Volkova, S. Lacis, and A. Meunier, “Magnetorheology: fluids, structures and rheology,” in Ferrofluids: Magnetically Controllable Fluids and Their Applications, S. Odenbach, Ed., vol. 594 of Lecture Notes in Physics, pp. 202–230, Springer, Berlin, Germany, 2002. View at: Publisher Site  Google Scholar
 J. RodríguezLópez, H. C. Shum, L. Elvira, F. M. de Espinosa, and D. A. Weitz, “Fabrication and manipulation of polymeric magnetic particles with magnetorheological fluid,” Journal of Magnetism and Magnetic Materials, vol. 326, pp. 220–224, 2013. View at: Publisher Site  Google Scholar
 S.F. Fu and X.Z. Wang, “Effects of magnetic field on shear stress of nematic liquid crystalline polymer under simple shear flow,” Journal of Polymer Science B: Polymer Physics, vol. 48, no. 17, pp. 1919–1926, 2010. View at: Publisher Site  Google Scholar
 M. T. Shaw, Introduction to Polymer Rheology, Wiley, New York, NY, USA, 2012.
 S. Sato, K. Oka, and A. Murakami, “Heat transfer behavior of melting polymers in laminar flow field,” Polymer Engineering & Science, vol. 44, no. 3, pp. 423–432, 2004. View at: Publisher Site  Google Scholar
 A. Ishak, “Similarity solutions for flow and heat transfer over a permeable surface with convective boundary condition,” Applied Mathematics and Computation, vol. 217, no. 2, pp. 837–842, 2010. View at: Publisher Site  Google Scholar  MathSciNet
 A. Aziz, “A similarity solution for laminar thermal boundary layer over a flat plate with a convective surface boundary condition,” Communications in Nonlinear Science and Numerical Simulation, vol. 14, no. 4, pp. 1064–1068, 2009. View at: Publisher Site  Google Scholar
 A. Aziz, “Hydrodynamic and thermal slip flow boundary layers over a flat plate with constant heat flux boundary condition,” Communications in Nonlinear Science and Numerical Simulation, vol. 15, no. 3, pp. 573–580, 2010. View at: Publisher Site  Google Scholar
 O. D. Makinde and P. O. Olanrewaju, “Combined effects of internal heat generation and buoyancy force on boundary layer flow over a vertical plate with a convective surface boundary condition,” Canadian Journal of Chemical Engineering, vol. 90, no. 5, pp. 1289–1294, 2012. View at: Publisher Site  Google Scholar
 D. Gupta, L. Kumar, O. Anwar Bég, and B. Singh, “Finite element simulation of mixed convection flow of micropolar fluid over a shrinking sheet with thermal radiation,” Proceedings of the Institution of Mechanical Engineers Part E: Journal of Process Mechanical Engineering, vol. 228, no. 1, pp. 61–72, 2014. View at: Publisher Site  Google Scholar
 G. Swapna, L. Kumar, O. Anwar Bég, and B. Singh, “Effect of thermal radiation on the mixed convection flow of a magnetomicropolar fluid past a continuously moving plate with a convective surface boundary condition,” Journal of Engineering Thermophysics, In press. View at: Google Scholar
 O. D. Makinde, K. Zimba, and O. A. Bég, “Numerical study of chemicallyreacting hydromagnetic boundary layer flow with Soret/Dufour effects and a convective surface boundary condition,” International Journal of Thermal and Environmental Engineering, vol. 4, no. 1, pp. 89–98, 2012. View at: Publisher Site  Google Scholar
 O. A. Bég, M. J. Uddin, M. M. Rashidi, and N. Kavyani, “Doublediffusive radiative magnetic mixed convective slip flow with Biot and Richardson number effects,” Journal of Engineering Thermophysics, vol. 23, no. 2, pp. 79–97, 2014. View at: Publisher Site  Google Scholar
 S. V. Subhashini, N. Samuel, and I. Pop, “Doublediffusive convection from a permeable vertical surface under convective boundary condition,” International Communications in Heat and Mass Transfer, vol. 38, no. 9, pp. 1183–1188, 2011. View at: Publisher Site  Google Scholar
 R. E. Powell and H. Eyring, “Mechanisms for the relaxation theory of viscosity,” Nature, vol. 154, no. 3909, pp. 427–428, 1944. View at: Publisher Site  Google Scholar
 T. Hayat, S. Asad, M. Mustafa, and A. Alsaedi, “Radiation effects on the flow of powelleyring fluid past an unsteady inclined stretching sheet with nonuniform heat source/sink,” PLoS ONE, vol. 9, no. 7, Article ID e103214, 2014. View at: Publisher Site  Google Scholar
 V. Sirohi, M. G. Timol, and N. L. Kalathia, “Numerical treatment of PowellEyring fluid flow past a 90 degree wedge,” Regional Journal of Energy, Heat and Mass Transfer, vol. 6, no. 3, pp. 219–228, 1984. View at: Google Scholar
 O. Adesanya and J. A. Gbadeyan, “Adomian decomposition approach to steady viscoelastic fluid flow with slip through a planar channel,” International Journal of Nonlinear Science, vol. 9, pp. 86–94, 2010. View at: Google Scholar
 V. R. Prasad, S. A. Gaffar, E. K. Reddy, and O. A. Bég, “Computational study of nonnewtonian thermal convection from a vertical porous plate in a nondarcy porous medium with biot number effects,” Journal of Porous Media, vol. 17, no. 7, pp. 601–622, 2014. View at: Publisher Site  Google Scholar
 L. Rosenhead, Laminar Boundary Layers, Oxford University Press, Oxford, UK, 1963.
 J. Peddieson, “Wedge and cone flows of viscoelastic liquids,” AIChE Journal, vol. 19, no. 2, pp. 377–379, 1973. View at: Publisher Site  Google Scholar
 E. M. Sparrow, R. Eichhorn, and J. L. Gregg, “Combined forced and free convection in a boundary layer flow,” Physics of Fluids, vol. 2, pp. 319–328, 1959. View at: Publisher Site  Google Scholar  MathSciNet
 T. Watanabe, K. Funazaki, and H. Taniguchi, “Theoretical analysis on mixed convection boundary layer flow over a wedge with uniform suction or injection,” Acta Mechanica, vol. 105, no. 1–4, pp. 133–141, 1994. View at: Publisher Site  Google Scholar
 N. G. Kafoussias and N. D. Nanousis, “Magnetohydrodynamic laminar boundarylayer flow over a wedge with suction or injection,” Canadian Journal of Physics, vol. 75, no. 10, pp. 733–745, 1997. View at: Publisher Site  Google Scholar
 N. D. Nanousis, “Theoretical magnetohydrodynamic analysis of mixed convection boundarylayer flow over a wedge with uniform suction or injection,” Acta Mechanica, vol. 138, no. 12, pp. 21–30, 1999. View at: Publisher Site  Google Scholar
 R. S. R. Gorla, “Unsteady heat transfer in laminar nonNewtonian boundary layer over a wedg,” AIChE Journal, vol. 28, no. 1, pp. 56–60, 1982. View at: Publisher Site  Google Scholar
 K.A. Yih, “Radiation effect on mixed convection over an isothermal wedge in porous media: the entire regime,” Heat Transfer Engineering, vol. 22, no. 3, pp. 26–32, 2001. View at: Publisher Site  Google Scholar
 M. M. Rashidi, M. T. Rastegari, M. Asadi, and O. A. Bég, “A study of nonNewtonian flow and heat transfer over a nonisothermal wedge using the homotopy analysis method,” Chemical Engineering Communications, vol. 199, no. 2, pp. 231–256, 2012. View at: Publisher Site  Google Scholar
 A. J. Chamkha, M. Mujtaba, A. Quadri, and C. Issa, “Thermal radiation effects on MHD forced convection flow adjacent to a nonisothermal wedge in the presence of a heat source or sink,” Heat and Mass Transfer, vol. 39, no. 4, pp. 305–312, 2003. View at: Google Scholar
 K.L. Hsiao, “MHD mixed convection for viscoelastic fluid past a porous wedge,” International Journal of NonLinear Mechanics, vol. 46, no. 1, pp. 1–8, 2011. View at: Publisher Site  Google Scholar
 A. Ishak, R. Nazar, and I. Pop, “Moving wedge and flat plate in a micropolar fluid,” International Journal of Engineering Science, vol. 44, no. 1819, pp. 1225–1236, 2006. View at: Publisher Site  Google Scholar  MathSciNet
 A. Ishak, R. Nazar, and I. Pop, “Moving wedge and flat plate in a powerlaw fluid,” International Journal of NonLinear Mechanics, vol. 46, no. 8, pp. 1017–1021, 2011. View at: Publisher Site  Google Scholar
 H. B. Keller, “Numerical methods in boundarylayer theory,” Annual Review of Fluid Mechanics, vol. 10, pp. 417–433, 1978. View at: Publisher Site  Google Scholar
 O. Anwar Bég, “Numerical methods for multiphysical magnetohydrodynamics,” in New Developments in Hydrodynamics Research, chapter 1, pp. 1–112, Nova Science, New York, NY, USA, 2012. View at: Google Scholar
 Y. Y. Lok, I. Pop, and D. B. Ingham, “Oblique stagnation slip flow of a micropolar fluid,” Meccanica, vol. 45, no. 2, pp. 187–198, 2010. View at: Publisher Site  Google Scholar  MathSciNet
 T.B. Chang, A. Mehmood, O. A. Bég, M. Narahari, M. N. Islam, and F. Ameen, “Numerical study of transient free convective mass transfer in a WaltersB viscoelastic flow with wall suction,” Communications in Nonlinear Science and Numerical Simulation, vol. 16, no. 1, pp. 216–225, 2011. View at: Publisher Site  Google Scholar  MathSciNet
 D. Srinivasacharya and K. Kaladhar, “Mixed convection flow of couple stress fluid in a nonDarcy porous medium with Soret and Dufour effects,” Journal of Applied Science and Engineering, vol. 15, no. 4, pp. 415–422, 2012. View at: Google Scholar
 S. A. Gaffar, V. R. Prasad, and O. A. Bég, “Computational analysis of magnetohydrodynamic free convection flow and heat transfer of nonNewtonian tangent hyperbolic fluid from a horizontal circular cylinder with partial slip,” International Journal of Applied and Computational Mathematics, 2015. View at: Publisher Site  Google Scholar
 V. R. Prasad, S. A. Gaffar, and O. A. Bég, “Heat and mass transfer of nanofluid from horizontal cylinder to micropolar fluid,” Journal of Thermophysics and Heat Transfer, vol. 29, no. 1, pp. 127–139, 2015. View at: Publisher Site  Google Scholar
 V. R. Prasad, S. A. Gaffar, E. K. Reddy, and O. A. Bég, “Numerical study of nonnewtonian boundary layer flow of jeffreys fluid past a vertical porous plate in a nondarcy porous medium,” International Journal of Computational Methods in Engineering Science and Mechanics, vol. 15, no. 4, pp. 372–389, 2014. View at: Publisher Site  Google Scholar
 R. M. Darji and M. G. Timol, “On invariance analysis of MHD boundary layer equations for nonNewtonian Williamson fluids,” International Journal of Advances in Applied Mathematics and Mechanics, vol. 1, pp. 10–19, 2014. View at: Google Scholar
 V. Singh and S. Agarwal, “Flow and heat transfer of Maxwell fluid with variable viscosity and thermal conductivity over an exponentially stretching sheet,” American Journal of Fluid Dynamics, vol. 3, no. 4, pp. 87–95, 2013. View at: Publisher Site  Google Scholar
 S. Nakamura, Applied Numerical Methods and Software, PrenticeHall, New York, NY, USA, 1995.
 O. A. Bég, T. A. Bég, H. S. Takhar, and A. Raptis, “Mathematical and numerical modeling of nonNewtonian thermohydrodynamic flow in nonDarcy porous media,” International Journal of Fluid Mechanics Research, vol. 31, no. 1, pp. 1–12, 2004. View at: Publisher Site  Google Scholar  MathSciNet
 I. Pop and S. Nakamura, “Laminar boundary layer flow of powerlaw fluids over wavy surfaces,” Acta Mechanica, vol. 115, no. 1–4, pp. 55–65, 1996. View at: Publisher Site  Google Scholar  MathSciNet
 R. S. R. Gorla, A. Slaouti, and H. S. Takhar, “Free convection in micropolar fluids over a uniformly heated vertical plate,” International Journal of Numerical Methods for Heat & Fluid Flow, vol. 8, no. 56, pp. 504–518, 1998. View at: Publisher Site  Google Scholar
 O. A. Bég, V. R. Prasad, and B. Vasu, “Numerical study of mixed bioconvection in porous media saturated with nanofluid containing oxytactic microorganisms,” Journal of Mechanics in Medicine and Biology, vol. 13, no. 4, Article ID 1350067, 2013. View at: Publisher Site  Google Scholar
 O. Anwar Bég, J. Zueco, M. Norouzi, M. Davoodi, A. A. Joneidi, and A. F. Elsayed, “Network and Nakamura tridiagonal computational simulation of electricallyconducting biopolymer micromorphic transport phenomena,” Computers in Biology and Medicine, vol. 44, no. 1, pp. 44–56, 2014. View at: Publisher Site  Google Scholar
 N. A. Yacob, A. Ishak, and I. Pop, “FalknerSkan problem for a static or moving wedge in nanofluids,” International Journal of Thermal Sciences, vol. 50, no. 2, pp. 133–139, 2011. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2015 S. Abdul Gaffar et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.