Journal of Applied Mathematics
Volume 2017 (2017), Article ID 1697135, 13 pages
https://doi.org/10.1155/2017/1697135

Research Article

Viscous Dissipation Effects on the Motion of Casson Fluid over an Upper Horizontal Thermally Stratified Melting Surface of a Paraboloid of Revolution: Boundary Layer Analysis

T. M. Ajayi, A. J. Omowaye, and I. L. Animasaun

Department of Mathematical Sciences, Federal University of Technology, Akure, Ondo State, Nigeria

Correspondence should be addressed to T. M. Ajayi; stillmetunde@gmail.com

Received 30 June 2016; Revised 20 October 2016; Accepted 6 November 2016; Published 4 January 2017

Academic Editor: Igor Andrianov

Copyright © 2017 T. M. Ajayi 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.

Abstract.

The problem of a non-Newtonian fluid flow past an upper surface of an object that is neither a perfect horizontal/vertical nor inclined/cone in which dissipation of energy is associated with temperature-dependent plastic dynamic viscosity is considered. An attempt has been made to focus on the case of two-dimensional Casson fluid flow over a horizontal melting surface embedded in a thermally stratified medium. Since the viscosity of the non-Newtonian fluid tends to take energy from the motion (kinetic energy) and transform it into internal energy, the viscous dissipation term is accommodated in the energy equation. Due to the existence of internal space-dependent heat source; plastic dynamic viscosity and thermal conductivity of the non-Newtonian fluid are assumed to vary linearly with temperature. Based on the boundary layer assumptions, suitable similarity variables are applied to nondimensionalized, parameterized and reduce the governing partial differential equations into a coupled ordinary differential equations. These equations along with the boundary conditions are solved numerically using the shooting method together with the Runge-Kutta technique. The effects of pertinent parameters are established. A significant increases in is guaranteed with when magnitude of is large. decreases with and .

1. Introduction

Within the last thirty years, the study of non-Newtonian fluid flow over a stretching surface has received significant attention due to its industrial applications. Such interest is fueled by its pertinent engineering applications in a number of fields as in spinning of filaments, continuous casting of metal, extrusion of polymers, crystal growing, glass fiber production, extrusion of plastic sheets, paper production, and process of condensation of metallic plates. Boundary layer analysis of fluid flow passing through a thick needle with variable diameter was investigated by Lee [1]. Historically, this can be referred to as the first report of flow adjacent to a surface with variable thickness where the effect of viscosity is highly significant. Thereafter, extensive studies were conducted on the boundary layer flows over a thin needle. Cebeci and Na [2] investigated the laminar free convection heat transfer from a needle. Ahmad et al. [3] examined the boundary layer flow over a moving thin needle with variable heat flux. The boundary layer flow over a stretching surface with variable thickness was analyzed by Fang et al. [4]. Recently, the flow of different fluids over an upper horizontal surface with variable thickness has been investigated extensively in [5–7]. This attracted Makinde and Animasaun [8, 9] to focus on the case of quartic autocatalysis kind of chemical reaction in the flow of an electrically conducting nanofluid containing gyrotactic-microorganism over an upper horizontal surface of a paraboloid of revolution in the presence and absence of thermophoresis and Brownian motion. In most cases, plastic dynamic viscosity of non-Newtonian Casson fluid tends to take energy away from the motion and transform it into internal energy. However, this area has been neglected.

In fluid mechanics, destruction of fluctuating velocity gradients due to viscous stresses is known as viscous dissipation. This partial irreversible processes is often referred to as transformation of kinetic energy into internal energy of the fluid (heating up the fluid due to viscosity since dissipation is high in the regions with large gradients). Pop [10] remarked that understanding the concept of energy dissipation and transport in nanoscale structures is of great importance for the design of energy-efficient circuits and energy-conversion systems. However, energy dissipation and transport of non-Newtonian fluid are also of importance to engineers and scientists. Motsumi and Makinde [11] examined the effects of viscous dissipation parameter (i.e., Eckert number), thermal diffusion, and thermal radiation on boundary layer flow of -water and -water nanofluids over a moving flat plate. In another theoretical study on combined effects of Newtonian heating and viscous dissipation parameter on boundary layer flow of copper and titania in water over stretchable wall, Makinde [12] reported an increase in the moving plate surface temperature and thermal boundary layer thickness. Depending on the admissible grouping of variables (parameterization), Eckert number and Brinkmann number () may be used to quantify viscous dissipation. In addition, unsteady mixed convection in the flow of air over a semi-infinite stretching sheet taking into account the effect of viscous dissipation was carried out by Abd El-Aziz [13]. Both at steady stage () and unsteady stage (), velocity and temperature of the flow increase with an increase in the magnitude of Eckert number. Recently, effects of viscous dissipation, Joule heating, and partial velocity slip on two-dimensional stagnation point flow were reported by Yasin et al. [14]. In another study conducted by Animasaun and Aluko [15], it is reported that when dynamic viscosity of air is assumed to vary linearly with temperature, normal negligible effect of Eckert number on velocity profiles will be noticed. Raju and Sandeep [16] focused on the motion of Casson fluid over a moving wedge with slip and observed a decrement in the temperature field with rising values of Eckert number. In the flow of non-Newtonian Casson fluid over an upper horizontal thermally stratified melting surface of a paraboloid of revolution, dissipation of smaller eddies due to molecular viscosity near the wall is significant. In the presence of constant magnetic field, electrically conducting Casson fluid flow over object that is neither a perfect horizontal/vertical nor inclined/cone is also an important issue.

The study of electrically conducting fluid flow is of considerable interest in modern metallurgical. This can be traced to the fact that most fluids in this sector are electrically conducting fluid. Historically, the first report on the motion of electrically conducting fluid in the presence of magnetic field was presented by Rossow [17]. Thereafter, Alfvén [18] reported that if a conducting liquid is placed in a constant magnetic field, every motion of the liquid generates a force called electromotive force (e.m.f.) which produces electric currents. One of the most significant importance of these contributions are its applications in engineering problems such as MHD generators, plasma studies, nuclear reactors, and geothermal energy extractions. Soundalgekar and Murty [19] investigated heat transfer in MHD flow with pressure gradient, suction, and injection. It was observed that an increase in the magnetic field parameter leads to an increase in fluid velocity, skin friction, rate of heat transfer, and a fall in temperature. Rajeswari et al. [20] observed that, due to the uniform magnetic field and suction at the wall of the surface, the concentration of the fluid decreases with the increase in chemical reaction parameter. Ghosh et al. [21] found that an increase in inclination of the applied magnetic field opposes primary flow and also reduces Grashof numbers. Das [22] concluded that increasing magnetic field and thermal radiation leads to deceleration of velocity but reverse is the effect for the melting parameter when the solid surface and the free stream move in the same direction. Motsa and Animasaun [23] presented the behavior of unsteady non-Darcian magnetohydrodynamic fluid flow past an impulsively using bivariate spectral local linearization analysis. Koriko et al. [24] illustrated the dynamics of two-dimensional magnetohydrodynamics (MHD) free convective flow of micropolar fluid along a vertical porous surface embedded in a thermally stratified medium. It was concluded that velocity profiles and microrotation profiles are strongly influenced by the magnetic field in the boundary layer, which decreases with an increase in the magnitude of magnetic parameter. Recently, theoretical investigation of MHD natural convection flow in vertical microchannel formed by two electrically nonconducting infinite vertical parallel plates and effects of MHD mixed convection on the flow through vertical pipe with time periodic boundary condition was presented explicitly by Jha and Aina [25, 26].

Several processes involving melting heating transfer in non-Newtonian fluids have promising applications in thermal engineering, such as melting permafrost, oil extraction, magma solidification, and thermal insulation. As such, a lot of experimental and theoretical work has been conducted in the kinetics of heat transfer accompanied with melting or solidification effect. The process of melting of ice placed in a hot stream of air at a steady state was first reported by Roberts [27]. Historically, this report can be referred to as the pioneering analysis of the melting phenomenon. Another novel report on melting phenomenon during forced convection heat transfer when an iceberg drifts in warm sea water was presented by Tien and Yen [28]. From their investigation, they observed that melting at the interface results in a decrease in the Nusselt number. Epstein and Cho [29] discussed the laminar film condensation on a vertical surface. Much later, melting heat transfer in a nanofluid boundary layer on a stretching circular cylinder was examined by Gorla et al. [30]. In the study of the effect of radiation on MHD mixed convection flow from a vertical plate embedded in a saturated porous media with melting, Adegbie et al. [31] reported that the Nusselt number decreases with increase in melting parameter. Adegbie et al. [31] stated that the temperature of UCM fluid flow over a melting surface is an increasing function of variable thermal conductivity parameter. Omowaye and Animasaun [32] investigated boundary layer analysis in the flow of upper convected Maxwell fluid flow. Due to the fact that temperature at the wall is zero, classical temperature-dependent viscosity and thermal conductivity linear models were modified to suit the case of both melting heat transfer and thermal stratification. In another study of micropolar fluid flow in the presence of temperature-dependent and space heat source, the analysis of the case where vortex viscosity is a constant function of temperature was reported in [33].

Within the past two decades, the effects of temperature-dependent viscosity on the fluid flow have become more important to engineers dealings with geothermal systems, crude oil extraction, and machinery lubrication. Due to friction and internal heat generated between two layers of fluid, viscosity and thermal conductivity of fluid substance may be affected by temperature; for more details see Batchelor [34], Lai and Kulacki [35], and Abd El-Aziz [36]. Proper consideration of this fact in the study on inherent irreversibility in a variable viscosity Couette flow by Makinde and Maserumule [37], numerical investigation of micropolar fluid flow over a nonlinear stretching sheet taking into account the effects of a temperature-dependent viscosity by Rahman et al. [38], effects of MHD on Casson fluid flow in the presence of Cattaneo-Christov heat flux by Malik et al. [39], fluid flow through a pipe with variation in viscosity by Makinde [40], Casson fluid flow within boundary layer over an exponentially stretching surface embedded in a thermally stratified medium by Animasaun [41], steady fully developed natural convection flow in a vertical annular microchannel having temperature-dependent viscosity in the presence of velocity slip and temperature jump at the annular microchannel surfaces by Jha et al. [42] have enhanced the body of knowledge on fluid flow, boundary layer analysis, and heat/mass transfer. In a recent experiment, Alam et al. [43] concluded that thermal boundary layer decreases with an increasing temperature-dependent viscosity. Hayat et al. [44] discussed the effect of variable thermal conductivity on the mixed convective flow over a porous medium stretching surface. In the article, the kinematics viscosity of Casson fluid was considered as a function depending on plastic dynamic viscosity, density, and Casson parameter and hence reported that increase in the magnitude of temperature-dependent viscosity parameter leads to an increase in fluid’s velocity. The above literature review shows that there exists no published article on the effects of viscous dissipation in the flow of non-Newtonian Casson fluid over a upper horizontal thermally stratified melting surface of a paraboloid of revolution.

2. Formulation of the Problem

Consider a steady, incompressible, laminar flow of an electrically conducting non-Newtonian (Casson) fluid over a melting surface on upper horizontal paraboloid of revolution in the presence of viscous dissipation and thermal stratification. The -axis is taken in the direction of motion and -axis is normal to the flow as shown in Figures 1(a) and 1(b). A uniform magnetic field of strength is applied normal to the flow. The induced magnetic field due to the motion of an electrically conducting Casson fluid is assumed to be so small; hence it is neglected. The stretching velocity is and the wall is assumed impermeable. It is further assumed that the immediate fluid layer adjacent to the surface is specified as , where . Using Boussinesq approximation, the difference in inertia is negligible but gravity is sufficiently strong to make the specific weight appreciably different between any two layers of Casson fluid on the surface; hence the rate flow on this kind of surface is referred to as free convection. The body force term suitable to induce the flow over a surface which is neither a perfect horizontal/vertical nor inclined/cone is

Following the theory stated in Casson [45] and boundary layer assumptions, the rheological equation for an isotropic Casson fluid flow together with heat transfer is of the form

Suitable boundary conditions governing the flow along upper horizontal surface of a paraboloid of revolution are

The formulation of the second term in boundary (5) states that the heat conducted to the melting surface on paraboloid of revolution is equal to the combination of heat of melting and the sensible heat required to raise the solid temperature to its stratified melting temperature (for details, see [29]). For lubricating fluids, heat generated by the internal friction and the corresponding increase in temperature affects the viscosity of the fluid; hence, it may not be realistic to constant function of plastic dynamic viscosity knowing fully well that space-dependent internal heat source is significant. In order to account for this variation, it is valid to consider the modified mathematical models of both the temperature-dependent viscosity and thermal conductivity models proposed in [46] and adopted in [31, 32] as

Meanwhile, the models are still in good agreement with experimental data of Batchelor [34]. It is worth mentioning that the first and second terms in (7) are valid and reliable since in this study, whereas thermal stratification () at the melting wall and at the free stream is defined as

Using similarity variable for temperature in (7) to simplify we obtain suitable temperature difference in flow past thermally stratified horizontal melting surface of paraboloid of revolution as

From thermal stratification models in (9), the following relations can be easily obtained

where is known as reference temperature. A significant difference between and can be easily obtained from (10). This can be traced to the fact that the linear stratification occurs at all points of “” on the wall and at all points of “” at the free stream as . In view of this, it is valid to define temperature-dependent viscous parameter by considering the second term in (10) since . Mathematically, the ratio of the first two terms in (10) can thus produce dimensionless thermal stratification parameter () as

The stream function and similarity variable are of the form

It is important to note that the stream function automatically satisfies continuity (2). The nonlinear partial differential equations (3) and (4) are reduced to the following nonlinear coupled ordinary differential equations

(a)

(b)

Figure 1: (a) The coordinate system of Casson fluid flow over upper horizontal thermally stratified melting surface of a paraboloid of revolution. (b) Graphical illustration of fluid domain and conversion of domain from to .

In (13) and (14), melting parameter , Prandtl number , magnetic parameter , Eckert number , temperature-dependent thermal conductivity parameter , space-dependent internal heat source parameter , skin friction coefficient , Nusselt number , and buoyancy parameter depending on volumetric-expansion coefficient due to temperature are defined as

where is the shear stress (skin friction) between Casson fluid and upper surface of horizontal paraboloid of revolution and is the heat flux at all points on the surface

In order to nondimensionalize the boundary conditions (5) and (6), it is pertinent to note that the minimum value of is not the starting point of the slot. This implies that all the conditions in (5) are not imposed at . As shown in Figures 1(a) and 1(b), it is obvious that it may not be realistic to say that at all points on upper horizontal melting surface of a paraboloid of revolution. Hence, it is not valid to set in similarity variable . Upon using the minimum value of which accurately corresponds to minimum value of similarity variable

This implies that, at the surface , the boundary condition suitable to scale the boundary layer flow is . The boundary condition becomes

Moreover, dimensionless governing equations (13) and (14) are depending on while the boundary conditions (18) and (19) are functions and/or derivatives depending on . In order to transform the domain from to it is valid to adopt and ; for more details, see Figure 1(b). Considering the fact that Prandtl number is strongly dependent on plastic dynamic viscosity and thermal conductivity, and it is assumed that both properties vary linearly with temperature; hence for more accurate analysis of boundary layer as suggested in [47, 48], the in (14) is

Equation (18) reveals that Prandlt number at free stream is denoted as . The final dimensionless governing equation (coupled system of nonlinear ordinary differential equation) is

Due to the fact that , the influence of temperature-dependent thermal conductivity on heat conduction during melting process diminishes. Dimensionless boundary condition reduces to

Upon substituting the similarity variables (12) and models of physical quantities (i.e., and ) at the wall into (16) we obtain

3. Numerical Solution

Numerical solutions of the boundary valued problem (21)–(24) are obtained using classical Runge-Kutta method with shooting techniques and MATLAB package (bvp5c). The boundary value problem cannot be solved on an infinite interval and it would be impractical to solve it for even a very large finite interval; hence at infinity is . Using the method of superposition by Na [49], the boundary value problem of ODE has been reduced to a system of five simultaneous equations of first order (IVP) for five unknowns following the method of superposition. In order to integrate the corresponding IVP, the values of and are required. However, such values do not exist after the nondimensionalization of the boundary conditions (5) and (6). The suitable guess values for and are chosen and then integration is carried out. The calculated values of and at infinity () are compared with the given boundary conditions in (24) and the estimated values and are adjusted to give a better approximation for the solution. Series of values for and are considered and applied with fourth-order classical Runge-Kutta method using step size . The above procedure is repeated until asymptotically converged results are obtained within a tolerance level of . It is very important to remark that setting , all profiles are compatible with the boundary layer theory and asymptotically satisfies the conditions at free stream as suggested by Pantokratoras [50]. It is worth mentioning that there exist no related published articles that can be used to validate the accuracy of the numerical results. In view of this, (21)–(24) can easily be solved using ODE solvers such as MATLAB’s bvp5c as explained in Kierzenka and Shampine [51] and Gökhan [52].

3.1. Verification of the Results

In order to verify the accuracy of the present analysis, the results of classical Runge-Kutta together with shooting have been compared with that of bvp5c solution for the limiting case when , , , , , , , , , , , and at various values of within the range . As shown in Table 1, the comparison in the above case is found to be in good agreement. This good agreement is an encouragement for further study of the effects of other parameters.

Table 1: Validation of numerical technique: comparison between the solutions of classical Runge-Kutta together with shooting (RK4SM) and MATLAB solver bvp5c for the limiting case.

(RK4SM)
(bvp5c)
(RK4SM)
(bvp5c)

0

0.5

1

4. Results and Discussion

The numerical computations have been carried out for various values of major parameters using the numerical scheme discussed in the previous section. This section presents the effect of different embedded physical parameters on the flow. Following Mustafa et al. [53], the ratio of momentum diffusivity to thermal diffusivity is considered to be unity (i.e., ) due to the fact that Casson fluid flow under consideration possesses substantial yield stress.

4.1. Influence of Velocity Index Parameter and Eckert Number

At metalimnion level of thermal stratification (), when , , , , , , , , , , and , it is worth noting that both vertical and horizontal velocities decrease with ; see Figures 2 and 3. Combination of these practical meanings of velocity index parameter justifies the decrease in both vertical and horizontal velocities profiles we obtained. Meanwhile, temperature distribution within the flow increases with from a few distance after wall till free stream; see Figure 4. With an increase in the magnitude of , Figure 5 reveals that temperature gradient profile increases near the wall and decreases thereafter as . It is very important to notice that velocity profile when perfectly satisfies free stream condition asymptotically. This can be traced to the fact that and as , . In addition, graphical illustration of the function which describes the immediate layer of Casson fluid next to an upper horizontal surface of a paraboloid of revolution against (i.e., and ) at various values of less than 1 and stretching velocity at the wall when stretching rate (i.e., are presented in Figure 6. It is noticed that as increases within , the thickness of the paraboloid of revolution decreases but corresponding influence of stretching on the flow is an increasing function. The variation in local skin friction coefficient and local Nusselt number which is proportional to local heat transfer rate as stated in (25) as a function of viscous dissipation term and velocity index parameter is shown in Figure 7 and Table 2. Figure 7 shows that, at a fixed value of , decreases with Eckert number . At constant value of , unequal decrease in with is also observed. In addition, decreases with at various values of . Table 2 shows that increases highly significant when magnitude of is large.

Table 2: Variations in Nusselt number with Eckert number at various values of velocity index parameter when , , , , , , , , , , , and .

When
When
When
When

Figure 2: Effect of on .

Figure 3: Effect of on .

Figure 4: Effect of on .

Figure 5: Effect of on .

Figure 6: Graphical illustrations of (i) ; (ii) stretching wall velocity at various values of .

Figure 7: Variations in with at various values of .

4.2. Influence of Non-Newtonian Casson Parameter , Thermal Stratification Parameter , and Eckert Number

Figures 8–12 illustrate the effect of increasing the magnitude of non-Newtonian Casson parameter on all the five profiles when , , , , , , , , , , , and . With an increase in the magnitude of , it is observed that vertical velocity decreases, horizontal velocity decreases, temperature distribution increases only within the fluid domain (), and the temperature gradient increases only near the wall. Figure 10 reveals that a distinct significant increase in is guarantee with an initial increase in from to . Physically, increase in the magnitude of non-Newtonian Casson parameter ( implies sharp transition in the flow behavior from non-Newtonian fluid flow to Newtonian fluid flow. In view of this, resistance in the fluid flow is produced. It is worth mentioning that an increase in implies a decrease in yield stress of the Casson fluid and increase in the magnitude of plastic dynamic viscosity .

Figure 8: Effect of on .

Figure 9: Effect of on .

Figure 10: Effect of on .

Figure 11: Effect of on .

Figure 12: Effect of on .

It is observed that the present study complements related studies on Casson fluid flow with temperature-dependent plastic dynamic viscosity on nonmelting surface; see Figures and in [41], Figures and in [54], and Figure reported by Jasmine Benazir et al. [55]. The relationship between non-Newtonian Casson parameter and Eckert number is sought for and illustrated graphically in Figures 13 and 14. Within , there exists no significant difference in with . As shown in Figure 13, when the magnitude of , a distinct significant increase in is observed due to an increase in the magnitude of viscous dissipation parameter. At small magnitude of , local Nusselt number () which is proportional to local heat transfer rate is found to be decreasing with . At large value of , increases with ; see Figure 14. The simulation was further extended to unravels the relationship between non-Newtonian Casson parameter , thermal stratification parameter , and local skin friction coefficients when , , , , , , , , , , and . It is revealed in Figure 15 that increases with at epilimnion stage which is known as the highest and warmest layer (). In addition, a significant decrease in is observed with an increase in at hypolimnion stage which can be referred to as the coolest layer. Mathematically, when , this implies that and maximum wall temperature at the wall explains the increase in since increase in corresponds to a decrease in yield stress .

Figure 13: Variations in with at various values of .

Figure 14: Variations in with at various values of .

Figure 15: Variations in with at various values of .

5. Conclusion

The boundary layer analysis of non-Newtonian Casson fluid flow over a horizontal melting surface embedded in a thermally stratified medium in the presence of viscous dissipation internal space heat source has been investigated numerically. The effects of the velocity power index, melting parameter, temperature-dependent viscous parameter, Eckert number, thermal conductivity, and magnetic interaction parameter were examined. Conclusions of the present analysis are as follows:

(1)Increase in the magnitude of velocity index parameter leads to a decrease in velocity and increase in temperature due to combine practical influence of the parameter.(2) decreases with at various values of within the interval . decreases with and .(3)At various values of viscous dissipation within and due to the nature of the flow past melting surface, a valid non-Newtonian parameter falls within the interval . Equivalent Newtonian fluid flow is guaranteed for .(4)Local skin friction coefficient increases negligible with when magnitude of is small. A significant increases in is guaranteed with when magnitude of is large.(5)In the case of Casson fluid flow over an upper horizontal thermally stratified melting surface of a paraboloid of revolution, decrease in horizontal velocity is guaranteed with an increase and .(6)With an increase in the magnitude of , the influence of stretching velocity at the wall on horizontal and vertical velocities is stronger than that of which describes the immediate fluid’s layer next to upper horizontal surface of a paraboloid of revolution due to melting heat transfer.

An extension of the present study to the case of Williamson and Prandtl fluid flow over an upper horizontal thermally stratified melting surface of a paraboloid of revolution is hereby recommended. For suitable parametrization to achieve a comparative study, see [7].

Nomenclature

: | Velocity component in direction
: | Velocity component in direction
: | Coefficient related to stretching sheet
: | Magnetic field strength
: | Skin friction coefficient
: | Specific heat at constant pressure
: | Distance along the surface
: | Distance normal to the surface
: | Parameter related to stretching sheet
: | Local Nusselt number
: | Prandtl number
: | Heat transfer
: | Dimensional fluid temperature
: | Stretching velocity at the wall
: | Thermal property of the fluid
: | Eckert number
: | Dimensionless vertical velocity
: | Thermal conductivity of the Casson fluid
: | Velocity power index
: | Melting temperature at wall
: | Ambient Temperature
: | Local Reynolds number
: | Non-Newtonian Casson parameter
: | Wall thickness parameter
: | Plastic dynamic viscosity
: | Latent heat of the fluid
: | Melting parameter
: | Dimensionless temperature
: | Kinematic viscosity
: | Temperature-dependent viscosity parameter
: | Temperature-dependent thermal conductivity parameter
: | Density of Casson fluid
: | Electrical conductivity of the fluid
: | Stream function.

Competing Interests

The authors declare that they have no competing interests.

References

L. L. Lee, “Boundary layer over a thin needle,” Physics of Fluids, vol. 10, no. 4, pp. 820–822, 1967.

T. Cebeci and T. Y. Na, “Laminar free-convection heat transfer from a needle,” Physics of Fluids, vol. 12, no. 2, pp. 463–465, 1969.

S. Ahmad, R. Nazar, and L. Pop, “Mathematical modeling of boundary layer flow over a moving thin needle with variable heat flux,” in Advances in Numerical Methods, vol. 11 of Lecture Notes in Electrical Engineering, pp. 43–54, Springer, Berlin, Germany, 2009.

T. Fang, J. Zhang, and Y. Zhong, “Boundary layer flow over a stretching sheet with variable thickness,” Applied Mathematics and Computation, vol. 218, no. 13, pp. 7241–7252, 2012.

I. L. Animasaun, “47nm alumina–water nanofluid flow within boundary layer formed on upper horizontal surface of paraboloid of revolution in the presence of quartic autocatalysis chemical reaction,” Alexandria Engineering Journal, vol. 55, no. 3, pp. 2375–2389, 2016.

I. L. Animasaun and N. Sandeep, “Buoyancy induced model for the flow of 36 nm alumina-water nanofluid along upper horizontal surface of a paraboloid of revolution with variable thermal conductivity and viscosity,” Powder Technology, vol. 301, pp. 858–867, 2016.

O. Abegunrin, S. Okhuevbie, and I. Animasaun, “Comparison between the flow of two non-Newtonian fluids over an upper horizontal surface of paraboloid of revolution: boundary layer analysis,” Alexandria Engineering Journal, vol. 55, no. 3, pp. 1915–1929, 2016.

O. D. Makinde and I. L. Animasaun, “Bioconvection in MHD nanofluid flow with nonlinear thermal radiation and quartic autocatalysis chemical reaction past an upper surface of a paraboloid of revolution,” International Journal of Thermal Sciences, vol. 109, pp. 159–171, 2016.

O. Makinde and I. Animasaun, “Thermophoresis and Brownian motion effects on MHD bioconvection of nanofluid with nonlinear thermal radiation and quartic chemical reaction past an upper horizontal surface of a paraboloid of revolution,” Journal of Molecular Liquids, vol. 221, pp. 733–743, 2016.

E. Pop, “Energy dissipation and transport in nanoscale devices,” Nano Research, vol. 3, no. 3, pp. 147–169, 2010.

T. G. Motsumi and O. D. Makinde, “Effects of thermal radiation and viscous dissipation on boundary layer flow of nanofluids over a permeable moving flat plate,” Physica Scripta, vol. 86, no. 4, Article ID 045003, 2012.

O. D. Makinde, “Analysis of Sakiadis flow of nanofluids with viscous dissipation and Newtonian heating,” Applied Mathematics and Mechanics (English Edition), vol. 33, no. 12, pp. 1545–1554, 2012.

M. Abd El-Aziz, “Unsteady mixed convection heat transfer along a vertical stretching surface with variable viscosity and viscous dissipation,” Journal of the Egyptian Mathematical Society, vol. 22, no. 3, pp. 529–537, 2014.

M. H. M. Yasin, A. Ishak, and I. Pop, “MHD stagnation-point flow and heat transfer with effects of viscous dissipation, Joule heating and partial velocity slip,” Scientific Reports, vol. 5, Article ID 17848, 2015.

I. L. Animasaun and O. B. Aluko, “Analysis on variable fluid viscosity of non-darcian flow over a moving vertical plate in a porous medium with suction and viscous dissipation,” IOSR Journal of Engineering, vol. 4, no. 8, pp. 18–32, 2014.

C. S. K. Raju and N. Sandeep, “MHD slip flow of a dissipative Casson fluid over a moving geometry with heat source/sink: a numerical study,” Acta Astronautica, 2016.

V. J. Rossow, “On the flow of electrically conducting fluid over a flat plate in the presence of a transverse magnetic field,” NACA Technical Report 3071, 1957.

H. Alfvén, “Existence of electromagnetic-hydrodynamic waves,” Nature, vol. 150, no. 3805, pp. 405–406, 1942.

V. M. Soundalgekar and T. V. R. Murty, “Heat transfer in MHD flow with pressure gradient, suction and injection,” Journal of Engineering Mathematics, vol. 14, no. 2, pp. 155–159, 1980.

R. Rajeswari, B. Jothiram, and V. K. Nelson, “Chemical reaction, heat and mass transfer on nonlinear MHD boundary layer flow through a vertical porous surface in the presence of suction,” Applied Mathematical Sciences, vol. 3, no. 49–52, pp. 2469–2480, 2009.

S. K. Ghosh, O. A. Bég, and A. Aziz, “A mathematical model for magnetohydrodynamic convection flow in a rotating horizontal channel with inclined magnetic field, magnetic induction and hall current effects,” World Journal of Mechanics, vol. 1, no. 3, pp. 137–154, 2011.

K. Das, “Radiation and melting effects on MHD boundary layer flow over a moving surface,” Ain Shams Engineering Journal, vol. 5, no. 4, pp. 1207–1214, 2014.

S. S. Motsa and I. L. Animasaun, “A new numerical investigation of some thermo-physical properties on unsteady mhd non-darcian flow past an impulsively started vertical surface,” Thermal Science, vol. 19, supplement 1, pp. 249–258, 2015.

O. K. Koriko, T. Oreyeni, A. J. Omowaye, and I. L. Animasaun, “Homotopy analysis of MHD free convective micropolar fluid flow along a vertical surface embedded in non-darcian thermally-stratified medium,” Open Journal of Fluid Dynamics, vol. 6, no. 3, pp. 198–221, 2016.

B. K. Jha and B. Aina, “Role of induced magnetic field on MHD natural convection flow in vertical microchannel formed by two electrically non-conducting infinite vertical parallel plates,” Alexandria Engineering Journal, vol. 55, no. 3, pp. 2087–2097, 2016.

B. K. Jha and B. Aina, “MHD mixed convection flow in a vertical pipe with time periodic boundary condition: steady periodic regime,” International Journal of Fluid Mechanics Research, vol. 43, no. 4, pp. 350–367, 2016.

L. Roberts, “On the melting of a semi-infinite body of ice placed in a hot stream of air,” Journal of Fluid Mechanics, vol. 4, pp. 505–528, 1958.

C. Tien and Y. Yen, “The effect of melting on forced convection heat transfer,” Journal of Applied Meteorology, vol. 4, no. 4, pp. 523–527, 1965.

M. Epstein and D. H. Cho, “Melting heat transfer in steady laminar flow over a flat plate,” Journal of Heat Transfer, vol. 98, no. 3, pp. 531–533, 1976.

R. S. Gorla, A. J. Chamkha, and E. Al-Meshaiei, “Melting heat transfer in a nanofluid boundary layer on a stretching circular cylinder,” Journal of Naval Architecture and Marine Engineering, vol. 9, no. 1, pp. 1–10, 2012.

K. S. Adegbie, A. J. Omowaye, A. B. Disu, and I. L. Animasaun, “Heat and mass transfer of upper convected maxwell fluid flow with variable thermo-physical properties over a horizontal melting surface,” Applied Mathematics, vol. 6, no. 8, pp. 1362–1379, 2015.

A. J. Omowaye and I. L. Animasaun, “Upper-convected maxwell fluid flow with variable thermo-physical properties over a melting surface situated in hot environment subject to thermal stratification,” Journal of Applied Fluid Mechanics, vol. 9, no. 4, pp. 1777–1790, 2016.

I. L. Animasaun, “Melting heat and mass transfer in stagnation point micropolar fluid flow of temperature dependent fluid viscosity and thermal conductivity at constant vortex viscosity,” Journal of the Egyptian Mathematical Society, 2016.

G. K. Batchelor, An Introduction to Fluid Dynamics, Cambridge University Press, Cambridge, UK, 1987.

F. C. Lai and F. A. Kulacki, “The effect of variable viscosity on convective heat transfer along a vertical surface in a saturated porous medium,” International Journal of Heat and Mass Transfer, vol. 33, no. 5, pp. 1028–1031, 1990.

M. Abd El-Aziz, “Temperature dependent viscosity and thermal conductivity effects on combined heat and mass transfer in MHD three-dimensional flow over a stretching surface with Ohmic heating,” Meccanica, vol. 42, article 375, 2007.

O. D. Makinde and R. L. Maserumule, “Thermal criticality and entropy analysis for a variable viscosity Couette flow,” Physica Scripta, vol. 78, no. 1, Article ID 015402, 2008.

M. M. Rahman, M. A. Rahman, M. A. Samad, and M. S. Alam, “Heat transfer in a micropolar fluid along a non-linear stretching sheet with a temperature-dependent viscosity and variable surface temperature,” International Journal of Thermophysics, vol. 30, no. 5, pp. 1649–1670, 2009.

M. Y. Malik, M. Khan, T. Salahuddin, and I. Khan, “Variable viscosity and MHD flow in Casson fluid with Cattaneo-Christov heat flux model: using Keller box method,” Engineering Science and Technology, an International Journal, vol. 19, no. 4, pp. 1985–1992, 2016.

O. D. Makinde, “Heat and mass transfer in a pipe with moving surface: effects of viscosity variation and energy dissipation,” Quaestiones Mathematicae, vol. 24, no. 1, pp. 93–104, 2001.

I. L. Animasaun, “Casson fluid flow with variable viscosity and thermal conductivity along exponentially stretching sheet embedded in a thermally stratified medium with exponentially heat generation,” Journal of Heat and Mass Transfer Research, vol. 2, no. 2, pp. 63–78, 2015.

B. K. Jha, B. Aina, and Z. Rilwanu, “Steady fully developed natural convection flow in a vertical annular microchannel having temperature dependent viscosity: an exact solution,” Alexandria Engineering Journal, vol. 55, no. 2, pp. 951–958, 2016.

M. S. Alam, M. M. Rahman, and M. A. Sattar, “Transient magnetohydrodynamic free convective heat and mass transfer flow with thermophoresis past a radiate inclined permeable plate in the presence of variable chemical reaction and temperature dependent viscosity,” Nonlinear Analysis: Modelling and Control, vol. 14, no. 1, pp. 3–20, 2009.

T. Hayat, S. A. Shehzad, M. Qasim, and A. Alsaedi, “Mixed convection flow by a porous sheet with variable thermal conductivityand convective boundary condition,” Brazilian Journal of Chemical Engineering, vol. 31, no. 1, pp. 109–117, 2014.

N. Casson, “A flow equation for pigment oil-suspensions of the printing ink type,” in Rheology of Disperse Systems, C. C. Mill, Ed., p. 84, Pergamon Press, 1959.

I. L. Animasaun, “Double diffusive unsteady convective micropolar flow past a vertical porous plate moving through binary mixture using modified Boussinesq approximation,” Ain Shams Engineering Journal, vol. 7, no. 2, pp. 755–765, 2015.

M. M. Rahman and K. M. Salahuddin, “Study of hydromagnetic heat and mass transfer flow over an inclined heated surface with variable viscosity and electric conductivity,” Communications in Nonlinear Science and Numerical Simulation, vol. 15, no. 8, pp. 2073–2085, 2010.

A. Pantokratoras, “Laminar free-convection over a vertical isothermal plate with uniform blowing or suction in water with variable physical properties,” International Journal of Heat and Mass Transfer, vol. 45, no. 5, pp. 963–977, 2002.

T. Y. Na, Computational Methods in Engineering Boundary Value Problems, Academic Press, New York, NY, USA, 2009.

A. Pantokratoras, “A common error made in investigation of boundary layer flows,” Applied Mathematical Modelling, vol. 33, no. 1, pp. 413–422, 2009.

J. Kierzenka and L. F. Shampine, “A BVP solver based on residual control and the Maltab PSE,” ACM Transactions on Mathematical Software, vol. 27, no. 3, pp. 299–316, 2001.

F. S. Gökhan, Effect of the Guess Function & Continuation Method on the Run Time of MATLAB BVP Solvers, Edited by C. M. Ionescu, 2011.

M. Mustafa, T. Hayat, I. Pop, and A. Aziz, “Unsteady boundary layer flow of a Casson fluid due to an impulsively started moving flat plate,” Heat Transfer—Asian Research, vol. 40, no. 6, pp. 563–576, 2011.

I. L. Animasaun, “Effects of thermophoresis, variable viscosity and thermal conductivity on free convective heat and mass transfer of non-darcian MHD dissipative Casson fluid flow with suction and nth order of chemical reaction,” Journal of the Nigerian Mathematical Society, vol. 34, no. 1, pp. 11–31, 2015.

A. Jasmine Benazir, R. Sivaraj, and M. M. Rashidi, “Comparison between casson fluid flow in the presence of heat and mass transfer from a vertical cone and flat plate,” Journal of Heat Transfer, vol. 138, no. 11, Article ID 112005, 2016.