Get 20M+ Full-Text Papers For Less Than $1.50/day. Start a 14-Day Trial for You and Your Team.

Learn More →

Double stratification effects on heat and mass transfer in unsteady MHD nanofluid flow over a flat surface

Double stratification effects on heat and mass transfer in unsteady MHD nanofluid flow over a... mutukuwinnie@gmail.com Department The focus of this work is to theoretically investigate the problem of double stratifica- of Mathematics, Kenyatta tion on heat and mass transfer in an unsteady hydromagnetic boundary layer flow of a University, Nairobi, Kenya Full list of author information nanofluid over a flat surface. The model employed for the nanofluid transport equa- is available at the end of the tions incorporate the effects of Brownian motion and thermophoresis in the presence article of thermal and solutal stratification. The governing nonlinear partial differential equa- tions and their associated boundary conditions are initially transformed into dimen- sionless form by using similarity variables, before being solved numerically by employ- ing the Runge–Kutta–Fehlberg fourth-order method with shooting technique. The effects of different controlling parameters, viz. solutal and thermal stratification, Lewis number, thermophoresis, Brownian motion, magnetic field and unsteadiness on the fluid velocity, temperature, skin friction coefficient, the local Nusselt number, and the local Sherwood number are graphically depicted and quantitatively discussed in detail taking into account the practical applications of each profile. It is noted that thermal stratification reduces the fluid temperature, while the solutal stratification reduces the nanoparticle concentration. Keywords: Heat transfer, Mass transfer, Unsteady, MHD, Nanofluid, Double stratification Background Magnetohydrodynamics (MHD) boundary layer flow of electrically conducting fluids has diverse industrial and engineering applications in fields such as nuclear reactors, geother - mal engineering, liquid metals and plasma flows, petroleum industries, boundary layer control in aerodynamics and crystal growth. Several authors have studied the problem of MHD boundary layer flow, heat and mass transfer about different surface geometries in electrically conducting fluids. Makinde [ 1] numerically analysed the influence of mag - netic field on the steady heat and mass flow of an electrically conducting fluid by mixed convection along a semi-infinite vertical porous plate with constant heat flux taking into account Soret and Dufour effects. The effect of the Hall current on the MHD natural con - vection flow from a vertical permeable flat plate with a uniform heat flux in the presence of a transverse magnetic field was analysed by Saha et al. [ 2]. Rout et al. [3] investigated MHD heat and mass transfer of chemically reacting fluid flow over a moving vertical plate in the presence of a heat source with convective boundary condition. © The Author(s) 2017. This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. Mutuku and Makinde Asia Pac. J. Comput. Engin. (2017) 4:2 Page 2 of 16 Stratification is the layering of a fluid system due to temperature differences and varia - tions in concentration or the presence of different fluids with varying densities. In prac - tical situations where the heat and mass transfer mechanisms occur simultaneously, it is important to analyse the effect of double stratification (stratification of medium with respect to thermal and concentration fields) on the convective transport in the fluid. The effect of double stratification is important in applications involving fluids’ convective transport where heat and mass transfer run concurrently. It is encountered in several classical problems, for example, reservoir mixing, internal waves, shear flow instability, internal hydraulics, turbulence and jets, plumes and wakes. Stratified fluids are omni - present in nature and their occurrence is very common in heterogeneous fluid bodies, such as thermal stratification of reservoirs and oceans, salinity stratification in estuar - ies, rivers, groundwater reservoirs, and oceans, heterogeneous mixtures in industrial food manufacturing processing, density stratification of the atmosphere, and uncount - able similar examples [4]. For instance, density variation causes thermal stratification in reservoirs which leads to a reduction in the vertical mixing of oxygen to the point that bottom water becomes anoxic through the action of biological processes. Preventing, predicting, and solving such a reservoir problem, though dependent on other limno- logical factors, require an understanding of the dynamics of stratified fluids. MHD-con - vective heat transfer in thermally stratified fluid occurs in many industrial applications and is an important aspect in the study of heat transfer. Higher energy efficiency and increased system performance can be achieved with better thermal stratification. This has led to intensified research in mass and heat transfer in thermally stratified medium [5–7]. However, the aforementioned investigations analysed natural convection flow of convectional heat transfer fluids. Recently with the advent of nanofluids—a term coined by Choi [8] referring to a liq - uid containing a suspension of nanometre-sized solid particle (nanoparticles)—research has diversified to convective heat and mass transfer of nanofluids in doubly stratified medium. Although the effect of double stratification on heat and mass transfer in a fluid is significant, very little work has been reported in the literature [9–11]. On the other hand, studies on unsteady boundary layer flow, mass transfer, and convective heat transfer in a nanofluid are the topics of contemporary research areas in fluid science and engineering, owing to its novelty [12–17]. The mass and heat transfer in a doubly strati - fied medium in the presence of a magnetic field is not only fundamentally interesting but also finds wide range of natural, industrial, and engineering applications. These applica - tions include heat rejection into the environment such as lakes, rivers, and seas; thermal energy storage systems such as solar ponds; heat transfer from thermal sources such as condensers of power plants. Despite its numerous applications, it is evident from the literature herein that research has neglected mass and heat transfer in unsteady MHD nanofluid flow taking into account thermophoresis and Brownian motion in a doubly stratified medium. Hence, the motivation of this paper is to bridge this knowledge gap. This study aims to extend the recent work of Olanrewaju and Makinde [18] to incorporate unsteadiness, mag- netic field, and the effects of double stratification on boundary layer nanofluid flow. The shooting iteration technique coupled with fourth-order Runge–Kutta integration scheme is employed to numerically solve the model nonlinear system of equations of Mutuku and Makinde Asia Pac. J. Comput. Engin. (2017) 4:2 Page 3 of 16 this particular problem. Numerical results are displayed graphically and discussed for different values of the various dimensionless parameters controlling the flow regime. It is hoped that the results obtained will not only provide useful information for applications but also serve as a basis for studying other analogous systems that arise in engineering and industrial applications. Model formulation We consider heat and mass transfer in an unsteady flow of an incompressible electrically conducting viscous nanofluid past an infinite nonconducting flat plate in the presence of a transversely imposed unsteady magnetic field of strength B(t). The sketch of the physi - cal configuration and coordinate system is shown in Fig.  1. The induced magnetic field is neglected, which is justified for industrial application MHD flows at small magnetic Reynolds number. It is also assumed that the external electrical field is zero and that the electric field due to the polarization of charges is also negligible. The temperature and nanoparticle concentration at the plate surface and at the free stream are assumed to vary with time and axial distance along the plate. Under the above assumptions with the boundary layer approximations, the continuity, momentum, energy, and nanoparticle concentration equations are as follows: ∂u ∂v + = 0, (1) ∂x ∂y 2 2 ∂u ∂u ∂u ∂U ∂U ∂ u σ B(t) (u − U ) ∞ ∞ ∞ + u + v = + U + υ − , ∞ (2) ∂t ∂x ∂y ∂t ∂x ∂y ρ ∂T ∂T ∂T ∂T ∂T ∂ T υ ∂u ∞ ∞ + u + v = + U + α + ∂t ∂x ∂y ∂t ∂x ∂y c ∂y 2 2 σB(t) (u − U ) ∂C ∂T D ∂T ∞ T − + τ D + , ρc ∂y ∂y T ∂y p 0,∞ (3) Fig. 1 Physical configuration Mutuku and Makinde Asia Pac. J. Comput. Engin. (2017) 4:2 Page 4 of 16 2 2 ∂C ∂C ∂C ∂C ∂C ∂ C D ∂ T ∞ ∞ T + u + v = + U + D + , ∞ B (4) 2 2 ∂t ∂x ∂y ∂t ∂x ∂y T ∂y 0,∞ with the boundary conditions u = 0, v = 0, T = T (x, t), C = C (x, t) at y = 0, w w (5) u → U (x, t), C → C (x, t), T → T (x, t), as y →∞, ∞ ∞ ∞ where u and v are the velocity components along the x and y directions, respectively; t is time; T and C are temperature and nanoparticle volume fraction at the plate surface, w w respectively; U , T , and C are the free stream velocity, temperature, and nanopar- ∞ ∞ ∞ ticle volume fraction, respectively; ρ is the density; α is the thermal diffusivity; υ is the kinematic viscosity; c is the specific heat at constant pressure; σ is the electrical conduc - tivity; τ = (ρC) /(ρC) is the ratio of the effective heat capacitance of the nanoparticle to p f that of the base fluid; D is the Brownian motion coefficient; D is the thermophoretic B T coefficient; T is the constant free stream temperature; and C is the nanoparticle vol- 0,∞ ume fraction. Following Ibrahim and Shankar [19], we prescribe the free stream velocity and the stratified free stream and plate surface temperature and nanoparticles concentration together with the unsteady magnetic field as follows: ax bx cx U (x, t) = , T (x, t) = T + , T (x, t) = T + , ∞ w 0,∞ ∞ 0,∞ 1 − t 1 − t 1 − t (6) B mx nx B(t) = √ , C (x, t) = C + , C (x, t) = C + , w 0,∞ ∞ 0,∞ 1 − t 1 − t 1 − t −1 where λt  <  1, a, b, c, m, and n are positive constants and λ with dimension (time) is unsteadiness frequency parameter. Positive values of b, c, n, and m correspond to the assisting flows and their negative values correspond to opposing flows. In order to sim - plify the mathematical analysis of the problem, we introduce the following dimension- less variables: a aυ η = y, ψ = xf (η), υ(1 − t) (1 − t) T − T T − T cx ∞ 0,∞ θ(η) = = −   , (7) T − T T − T (1 − t) T − T w 0,∞ w 0,∞ w 0,∞ C − C C − C nx ∞ ∞ φ(η) = = − , C − C C − C (1 − t) C − C w 0,∞ w 0,∞ w 0,∞ where η is the similarity variable and ψ is the stream function defined as ∂ψ ∂ψ u = and v =− (8) ∂y ∂x Substituting Eqs. (7) and (8) into Eqs. (1) and (4), we obtain following system of simi- larity equations: ′′′ ′′ ′ ′ ′′ ′ f + ff − f − A f + f − 1 − Ha f − 1 + 1 = 0 (9) 2 Mutuku and Makinde Asia Pac. J. Comput. Engin. (2017) 4:2 Page 5 of 16 ηθ ′′ ′ ′ ′ θ − A Pr + θ + Pr(f θ + β − f θ − f β ) 1 1 (10) 2 2 ′′ ′ 2 ′ ′ ′ + Pr Ec[ f + Ha(f − 1) ]+ N θ φ + N θ = 0 b t ηφ N ′′ ′ ′ ′ ′′ φ − ALe + φ + Le(f φ + β − f φ − f β ) + θ = 0, 2 2 (11) 2 N subject to the boundary conditions f (0) = 0, f (0) = 0, θ(0) = 1 − β , φ(0) = 1 − β , 1 2 (12) (∞) = 1, θ(∞) = 0, φ(∞) = 0, where primes denote differentiation with respect to η and Ha is Hartmann number, Pr is Prandtl number, N is the Brownian motion parameter, N is the thermophoresis param- b t eter, Ec is Eckert number, Le is Lewis number, A is the unsteadiness parameter, β is the thermal stratification parameter, and β is the solutal stratification parameter. The above parameters are defined as follows: σB υ τD C − C τD T − T B w 0,∞ T w 0,∞ Ha = , Pr = , N = , N = , aρ α υ T υ o,∞ (13) U υ  c n Ec = , Le = , A = , β = , β = . 1 2 C T − T D a b m p w 0,∞ The significant physical quantities of engineering and industrial interest in this study are the skin friction coefficient C , Nusselt number Nu, and Sherwood number Sh defined, respectively, as τ xq xq w w m C = , Nu = , Sh = , (14) ρu k(T − T ) D (C − C ) w 0,∞ B w 0,∞ where τ , q , and q are the skin friction, surface heat flux, and the surface mass flux, w w m respectively, defined by ∂u ∂T ∂C τ = μ , q =−k , q =−D . w w m B (15) ∂y ∂y ∂y y=0 y=0 y=0 where µ and k are the dynamic viscosity and the thermal conductivity, respectively. Using Eqs. (15) into (14), we get 1/2 ′′ −1/2 ′ −1/2 ′ Re C = f (0), Re Nu =−θ (0), Re Sh = −φ (0), (16) x x x which are the local skin friction C , local Nusselt number Nu, and local Sherwood num- ber Sh, respectively, and Re  = U x/υ is the local Reynolds number. x ∞ Mutuku and Makinde Asia Pac. J. Comput. Engin. (2017) 4:2 Page 6 of 16 Numerical procedure The system of nonlinear differential Eqs.  (9)–(11) with the boundary conditions (12) has been solved numerically using shooting technique coupled with the fourth-order Runge–Kutta method and a modified version of the Newton–Raphson algorithm. Let ′ ′′ ′ ′ x = f , x = f , x = f , x = θ, x = θ , x = φ, x = φ . 1 2 3 4 5 6 7 (17) Then the set of higher order boundary value equations with their respective boundary conditions is reduced to a system of first-order differential equation with appropriate initial conditions, respectively, x = x , x = x , � � ′ 2 3 x =−x x + x + A x + x − 1 + Ha x − 1 − 1, ( ) 1 3 2 2 3 2 x = x , � � ηx x = A Pr + x − Pr(x x + β − x x − x β ) 4 1 5 1 2 4 2 1 � � �� (18) ′ ′ 2 − Pr Ec x + Ha x − 1 − N x x − N x , b 5 7 t 3 2 5 x = x , � � ηx x = A Le + x + Le(x x + β − x x − x β ) 6 1 7 2 2 6 2 2   � � ηx   A Pr + x − Pr(x x + β − x x − x β N 4 1 5 1 2 4 2 1 − , � � ��   N ′ ′ 2 − Pr Ec x + Ha x − 1 − N x x − N x b 5 7 t 3 2 5 subject to the initial conditions x (0) = 0, x (0) = 0, x (0) = s , x (0) = 1 − β , 1 2 3 1 4 1 (19) x (0) = s , x (0) = 1 − β , x = s , 5 2 6 2 7 3 By applying the shooting method, the unspecified initial conditions S , S , and S in 1 2 3 (19) are assumed and (18) integrated numerically as an initial valued problem to a given terminal point. The accuracy of the assumed missing initial conditions was checked by comparing the calculated value of the dependent variable at the terminal point with its given value there. If differences exist, improved values of the missing initial conditions are obtained and the process repeated. The method is programmed in MAPLE with the −7 step size of ∆η  = 0.001 selected to be satisfactory for a convergence criterion of 10 in nearly all cases. The maximum value of η to each group of parameters was determined when the values of the unknown boundary conditions at η = 0 do not change to success- −7 ful loop with error less than 10 . Results and discussion From the numerical computation process, the effects of the pertinent parameters, namely, unsteadiness parameter A, Lewis number Le, Brownian parameter N , ther- mophoresis parameter N , Hartmann number Ha, Eckert number Ec, thermal strati- fication parameter β , and solutal stratification parameter β on the fluid velocity, 1 2 temperature, nanoparticles volume as well as the skin friction coefficient, the local Mutuku and Makinde Asia Pac. J. Comput. Engin. (2017) 4:2 Page 7 of 16 Nusselt number, and the local Sherwood number are illustrated graphically (Figs.  2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14) and analysed quantitatively highlighting their industrial and engineering applications. In order to ascertain the validity of our numerical proce- dure, the special case of heat transfer in flow of a convectional fluid in the absence of buoyancy force, flow unsteadiness, magnetic field, and double stratification effect (i.e. Fig. 2 Eec ff ts of A and Ha on dimensionless velocity Fig. 3 Eec ff ts of Ec , Nt, and Nb on dimensionless temperature Mutuku and Makinde Asia Pac. J. Comput. Engin. (2017) 4:2 Page 8 of 16 Fig. 4 Eec ff ts of Ha on dimensionless temperature Fig. 5 Eec ff ts of A on dimensionless temperature A = β  = β  = Ha = φ = 0) over a surface can be compared with that of Aman et al. [20], 1 2 as shown in Table  1 and excellent agreement was established. This favourable compari - son lends confidence to the numerical results reported subsequently. Mutuku and Makinde Asia Pac. J. Comput. Engin. (2017) 4:2 Page 9 of 16 Fig. 6 Eec ff ts of β on dimensionless temperature Fig. 7 Eec ff ts of A and Ha on dimensionless nanoparticle concentration Velocity profiles Figure  2 shows the effect of unsteadiness and magnetic field strength on the nanoflu - ids’ velocity profiles. It can be seen that there is a gradual increase in the fluid veloc - ity, from zero at the plate surface to a maximum value as it approaches the free stream far away from the plate surface thus satisfying the boundary conditions. It is observed Mutuku and Makinde Asia Pac. J. Comput. Engin. (2017) 4:2 Page 10 of 16 Fig. 8 Eec ff ts of Le , Nt, and Nb on dimensionless nanoparticle concentration Fig. 9 Eec ff ts of β and Ec on dimensionless nanoparticle concentration that both the fluid velocity and the hydrodynamic boundary layer thickness decrease as the unsteadiness parameter (A) and Hartmann number (Ha) increases. It is well known that application of a transverse magnetic field orthogonal to the flow direction has a ten - dency to create a drag force known as Lorentz force which tends to resist the fluid flow and thus reducing its velocity. Mutuku and Makinde Asia Pac. J. Comput. Engin. (2017) 4:2 Page 11 of 16 Fig. 10 Eec ff ts of Ha and A on local skin friction Fig. 11 Eec ff ts of Ec , Ha, β , Nb, and Nt on local Nusselt number Temperature profiles Figures 3, 4, 5, 6 show the effects of various parameters on the temperature profile. Gen - erally, the temperature is the maximum at the plate surface but decreases exponentially to zero far away from the plate surface satisfying the free stream conditions. As shown in Mutuku and Makinde Asia Pac. J. Comput. Engin. (2017) 4:2 Page 12 of 16 Fig. 12 Eec ff ts of Le and A of local Nusselt number Fig. 13 Eec ff ts of Ha, Ec , β , Nt, and Le of local Sherwood number Fig. 3, an increase in Brownian parameter N , thermophoresis parameter N , and Eckert b t number Ec leads to an increase in both the fluid temperature and the thermal boundary layer thickness. This can be attributed to the additional heating due to viscous dissipa - tion, friction, and rapid collision of the nanoparticles as a result of Brownian motion and thermophoresis. Figure 4 shows that increasing Hartmann number Ha initially increases Mutuku and Makinde Asia Pac. J. Comput. Engin. (2017) 4:2 Page 13 of 16 Fig. 14 Eec ff ts of β and Nb of local Sherwood number Table 1 Computations showing comparison with Aman et al. [20] for different values of  Pr when A = β  = β  = Ha = φ = 0 1 2 ′′ ′ ′′ ′ Pr f (0) −θ (0) f (0) −θ (0) Amani et al. [20] Amani et al. [20] present present 0.1 2.50997 0.26815 2.50997 0.26815 0.3 2.50997 0.40338 2.50997 0.40338 0.7 2.50997 0.52965 2.50997 0.52965 6.2 2.50997 0.59014 2.50997 0.59014 the fluid temperature at the plate surface, but as the fluid moves from the plate surface the temperature decreases. The increase in temperature is due to ohmic heating as a result of the magnetic effect. On the other hand, as illustrated in Figures  5 and 6, there is a decrease in both the temperature and the thermal boundary layer thickness with an increase in thermal stratification parameter β and unsteadiness parameter A. When the thermal stratification effect is taken into consideration, the effective temperature differ - ence between the plate and the ambient fluid decreases; therefore, the temperature is reduced. Nanoparticle concentration profiles Figures 7, 8, 9 show the effects of various parameters on the nanoparticle concentration profile. The nanoparticle concentration is the highest at the plate surface but decreases exponentially to zero free stream value satisfying the boundary conditions. Unlike in the case of temperature, a decrease in the nanoparticle concentration is noted with an increase in unsteadiness parameter A, Lewis number Le, Brownian parameter N , b Mutuku and Makinde Asia Pac. J. Comput. Engin. (2017) 4:2 Page 14 of 16 thermophoresis parameter N , Hartmann number Ha, Eckert number Ec, and solutal stratification parameter β . Eec ff ts of parameter variations on C , Nu, Sh Figures 10, 11, 12, 13, 14 depict the various pertinent parameters at the plate surface for the skin friction coefficient, the local Nusselt number (rate of heat transfer), and the local Sherwood number with variation in the values of thermophysical parameters embed- ded in the problem. From Fig.  10, it can be seen that increasing both the unsteadiness parameter A and Hartmann number Ha leads to an increase in the skin friction coeffi - cient. This is as expected since both parameters causes an overshoot of the fluid towards the plate surface. From Fig.  11, it is noted increasing the Lewis number Le, Brownian parameter N , thermophoresis parameter N , Hartmann number Ha, Eckert number Ec, b t and thermal stratification parameter β leads to a decrease in the local Nusselt number which implies that there is reduced heat transfer rate at the plate surface. As seen earlier, the same parameters led to an increase to the fluid temperature; thus, it is expected that with increased temperature within the fluid the rate of heat transfer from the plate to the fluid will be reduced. On the other hand, an increase in the unsteadiness parameter A increases the local Nusselt number, as observed in Fig.  12. It is noted in Figs.  13 and 14 that increasing the Hartmann number Ha, Eckert number Ec, Lewis number Le, ther- mophoresis parameter N , and thermal stratification parameter β increases the Sher- t 1 wood number, while a decrease in the local Sherwood number is realized by increasing the Brownian parameter N and solutal stratification parameter β . It well known that b 2 at high Lewis number the nanoparticle concentration is low, since the concentration at the plate surface is higher than in the fluid, mass transfer from the plate to the fluid then occurs. Conclusions The combined effects of thermal stratification, solutal stratification, Brownian motion, and thermophoresis on an unsteady hydromagnetic boundary layer flow of a nanofluid over a flat surface are numerically investigated. Using a similarity variable, the governing nonlinear partial differential equations were obtained and transformed to a set of non - linear differential equations. Numerical results are obtained, then presented graphically for the velocity, temperature, and nanoparticle concentration profiles as well as skin fric - tion coefficient, reduced Nusselt number, and reduced Sherwood number. Our results reveal, among other conclusions, the following: • The fluid overshoots towards the plate surface with increasing A and Ha, thus reduc - ing both the fluid velocity and the hydrodynamic boundary layer thickness. • An increase in Ec, Nt, Nb, and Ha yields an increase in both the thermal boundary layer thicknesses and temperature, whereas the reverse is noted with increases in A and β . • A decrease in the nanoparticle concentration is noted with increase in A, Le, Nb, Nt, Ha, Ec, and β . • Increasing A and Ha increase the skin friction. Mutuku and Makinde Asia Pac. J. Comput. Engin. (2017) 4:2 Page 15 of 16 • A decrease in Nusselt number is realized with increased Ec, Ha, Nb, and β , whereas it increases with increasing values of A. • The local Sherwood number increases with increase in Ha, β , Le, Nt, and Ec but reduces with increased β and Nb. Abbreviations (u, v): velocity components in x- and y-direction; A: unsteadiness parameter (λ/a); B : magnetic field strength (T); C: nanoparticle volume fraction; C : skin friction coefficient (τ /U ρ); C : constant free stream volume fraction; D : Brownian diffusion coeffi - w ∞ 0,∞ B cient; D : thermophoresis diffusion coefficient; Ec: Eckert number ( U /C (T  − T )); T ∞ p w 0,∞ 2 2 Ha: magnetic parameter (σB B /αρ); k: thermal conductivity (W/mK); Le: Lewis num- ber (υ/D ); N : Brownian motion parameter (τD (C   −  C )/υ); N : thermophoresis m b B w 0,∞ t parameter (τD (T  − T )/T υ); Pr: Prandtl number (υ/α); Nu: Nusselt number (xq / T w 0,∞ 0,∞ w k(T  − T )); Sh: Sherwood number (xq /D (C  − C )); T: temperature of fluid (K) w 0,∞ m B w 0,∞ Greek symbols ψ: stream function; η: similarity variable; β : thermal stratification parameter; β : solutal 1 2 stratification parameter; υ: kinematic viscosity of fluid (kg/m  s); ρ: density of fluid (kg/ m ); φ: dimensionless concentration function; α: thermal diffusivity(W/m  k); σ: electrical conductivity of fluid (S/m); θ: dimensionless temperature; (c ) : heat capacity of fluid (J/ p f kg  K); (c ) : heat capacity of nanofluid; τ: parameter defined by (ρC) /(ρC) ; µ: dynamic p p p f viscosity of fluid (kg/m s) Subscripts ∞: condition at the free stream; w: condition at the surface Authors’ contributions Both authors were actively involved in the formulation of the problem, literature review, and simulations. Both authors read and approved the final manuscript. Author details 1 2 Department of Mathematics, Kenyatta University, Nairobi, Kenya. Faculty of Military Science, Stellenbosch University, Cape Town, South Africa. Acknowledgements None. Competing interests The authors declare that they have no competing interests. Funding There is no funding towards the publication of this manuscript. Ethics approval and consent to participate This study does not involve human participants, human data, or human tissue. Publisher’s Note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations. Received: 27 March 2017 Accepted: 3 May 2017 References 1. Makinde OD (2011) On MHD mixed convection with Soret and Dufour effects past a vertical plate embedded in a porous medium. LAAR 41(1):63–68 Mutuku and Makinde Asia Pac. J. Comput. Engin. (2017) 4:2 Page 16 of 16 2. Saha LK, Siddiqa S, Hossain MA (2011) Eec ff t of hall current on MHD natural convection flow from vertical perme - able flat plate with uniform surface heat flux. Appl Math Mech Engl 3(9):1127–1146 3. Rout BR, Panda SK, Panda S (2013) MHD heat and mass transfer of chemical reaction fluid flow over a mov- ing vertical plate in presence of heat source with convective surface boundary condition. Int J Chem Eng. doi:10.1155/2013/296834 4. Socolofsky AS, Jirka GH (2004) Stratified flow and buoyant mixing. Eng Lect Notes. Texas A&M University, College Station 5. Singh G, Sharma PR, Chamkha AJ (2010) Eec ff t of thermally stratified ambient fluid on MHD convective flow along a moving non-isothermal vertical plate. Int J Phys Sci 5(3):208–215 6. Saha SC, Hossain MA (2004) Natural convection flow with combined buoyancy effects due to thermal and mass diffusions in a thermally stratified media. NonLinear Anal 9(1):89–102 7. Mukhopadhyay S, Mondal C, Gorla SR (2012) Eec ff ts of thermal stratification on flow and heat transfer past a porous vertical stretching surface. Heat Mass Transf 48:915–921 8. Choi SUS (1995) Enhancing thermal conductivity of fluids with nanoparticle. In: Siginer DA, Wang HP (eds) Devel- opment and applications of non-Newtonian flow. The American Society of Mechanical Engineers, New York, pp 227–233 9. Srinivasacharya D, Upendar M (2013) Eec ff t of double stratification on MHD free convection in a micropolar fluid. J Egypt Math Soc. doi:10.1016/j.joems.2013.02.006 10. Narayana L, Murthy PA (2006) Free convective heat and mass transfer in a doubly stratified non-Darcy porous medium. J Heat Transf 128:1204–1212 11. Ibrahim W, Shankar B (2013) MHD boundary layer flow and heat transfer of a nanofluid past a permeable stretching sheet with velocity, thermal and solutal slip boundary conditions. J Comput Fluids 75:1–10 12. Makinde OD, Olanrewaju PO, Charles WM (2011) Unsteady convection with chemical reaction and radiative heat transfer past a flat porous plate moving through a binary mixture. Afr Mat 22(2011):65–78 13. Chamkha AJ, Rashady AM, Al-Meshaieiz E (2011) Melting effect on unsteady hydromagnetic flow of a nanofluid past a stretching sheet. Int J Chem Reactor Eng 9(113):1 14. Makinde OD (2012) Computational modelling of MHD unsteady flow and heat transfer towards a flat plate with Navier slip and Newtonian heating. Braz J Chem Eng 29(1):159–166 15. Khan MS, Karim I, Ali LE, Islam A (2012) Unsteady MHD free convection boundary layer flow of a nanofluid along a stretching sheet with thermal radiation and viscous dissipation effects. Int Nano Lett 2:24 16. Vajravelu K, Prasad KV, Ng C (2013) Unsteady convective boundary layer flow of a viscous fluid at a vertical surface with variable fluid properties. Nonlinear Anal Real 14:455–464 17. Makinde OD, Chinyoka T (2013) Numerical investigation of buoyancy effects on hydromagnetic unsteady flow through a porous channel with suction/injection. J Mech Sci Tech 27(5):1557–1568 18. Olanrewaju AM, Makinde OD (2013) On boundary layer stagnation point flow of a nanofluid over a permeable flat surface with Newtonian heating. Chem Eng Commun 200:836–852 19. Ibrahim W, Shanker B (2012) Unsteady MHD boundary layer flow and heat transfer due to stretching sheet in the presence of heat source or sink by quasi linearization technique. Int J Appl Math Mech 8(7):18–30 20. Aman F, Ishak A, Pop I (2011) Mixed convection boundary layer flow near stagnation point on vertical surface with slip. Appl Math Mech 32(12):1599–1606 Submit your next manuscript to BioMed Central and take full advantage of: • Convenient online submission • Thorough peer review • No space constraints or color figure charges • Immediate publication on acceptance • Inclusion in PubMed, CAS, Scopus and Google Scholar • Research which is freely available for redistribution Submit your manuscript at www.biomedcentral.com/submit http://www.deepdyve.com/assets/images/DeepDyve-Logo-lg.png Asia Pacific Journal on Computational Engineering Springer Journals

Double stratification effects on heat and mass transfer in unsteady MHD nanofluid flow over a flat surface

Loading next page...
 
/lp/springer-journals/double-stratification-effects-on-heat-and-mass-transfer-in-unsteady-x2jAkNVzlf
Publisher
Springer Journals
Copyright
Copyright © 2017 by The Author(s)
Subject
Engineering; Theoretical and Applied Mechanics; Computational Science and Engineering; Classical and Continuum Physics; Mathematical Applications in the Physical Sciences
eISSN
2196-1166
DOI
10.1186/s40540-017-0021-2
Publisher site
See Article on Publisher Site

Abstract

mutukuwinnie@gmail.com Department The focus of this work is to theoretically investigate the problem of double stratifica- of Mathematics, Kenyatta tion on heat and mass transfer in an unsteady hydromagnetic boundary layer flow of a University, Nairobi, Kenya Full list of author information nanofluid over a flat surface. The model employed for the nanofluid transport equa- is available at the end of the tions incorporate the effects of Brownian motion and thermophoresis in the presence article of thermal and solutal stratification. The governing nonlinear partial differential equa- tions and their associated boundary conditions are initially transformed into dimen- sionless form by using similarity variables, before being solved numerically by employ- ing the Runge–Kutta–Fehlberg fourth-order method with shooting technique. The effects of different controlling parameters, viz. solutal and thermal stratification, Lewis number, thermophoresis, Brownian motion, magnetic field and unsteadiness on the fluid velocity, temperature, skin friction coefficient, the local Nusselt number, and the local Sherwood number are graphically depicted and quantitatively discussed in detail taking into account the practical applications of each profile. It is noted that thermal stratification reduces the fluid temperature, while the solutal stratification reduces the nanoparticle concentration. Keywords: Heat transfer, Mass transfer, Unsteady, MHD, Nanofluid, Double stratification Background Magnetohydrodynamics (MHD) boundary layer flow of electrically conducting fluids has diverse industrial and engineering applications in fields such as nuclear reactors, geother - mal engineering, liquid metals and plasma flows, petroleum industries, boundary layer control in aerodynamics and crystal growth. Several authors have studied the problem of MHD boundary layer flow, heat and mass transfer about different surface geometries in electrically conducting fluids. Makinde [ 1] numerically analysed the influence of mag - netic field on the steady heat and mass flow of an electrically conducting fluid by mixed convection along a semi-infinite vertical porous plate with constant heat flux taking into account Soret and Dufour effects. The effect of the Hall current on the MHD natural con - vection flow from a vertical permeable flat plate with a uniform heat flux in the presence of a transverse magnetic field was analysed by Saha et al. [ 2]. Rout et al. [3] investigated MHD heat and mass transfer of chemically reacting fluid flow over a moving vertical plate in the presence of a heat source with convective boundary condition. © The Author(s) 2017. This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. Mutuku and Makinde Asia Pac. J. Comput. Engin. (2017) 4:2 Page 2 of 16 Stratification is the layering of a fluid system due to temperature differences and varia - tions in concentration or the presence of different fluids with varying densities. In prac - tical situations where the heat and mass transfer mechanisms occur simultaneously, it is important to analyse the effect of double stratification (stratification of medium with respect to thermal and concentration fields) on the convective transport in the fluid. The effect of double stratification is important in applications involving fluids’ convective transport where heat and mass transfer run concurrently. It is encountered in several classical problems, for example, reservoir mixing, internal waves, shear flow instability, internal hydraulics, turbulence and jets, plumes and wakes. Stratified fluids are omni - present in nature and their occurrence is very common in heterogeneous fluid bodies, such as thermal stratification of reservoirs and oceans, salinity stratification in estuar - ies, rivers, groundwater reservoirs, and oceans, heterogeneous mixtures in industrial food manufacturing processing, density stratification of the atmosphere, and uncount - able similar examples [4]. For instance, density variation causes thermal stratification in reservoirs which leads to a reduction in the vertical mixing of oxygen to the point that bottom water becomes anoxic through the action of biological processes. Preventing, predicting, and solving such a reservoir problem, though dependent on other limno- logical factors, require an understanding of the dynamics of stratified fluids. MHD-con - vective heat transfer in thermally stratified fluid occurs in many industrial applications and is an important aspect in the study of heat transfer. Higher energy efficiency and increased system performance can be achieved with better thermal stratification. This has led to intensified research in mass and heat transfer in thermally stratified medium [5–7]. However, the aforementioned investigations analysed natural convection flow of convectional heat transfer fluids. Recently with the advent of nanofluids—a term coined by Choi [8] referring to a liq - uid containing a suspension of nanometre-sized solid particle (nanoparticles)—research has diversified to convective heat and mass transfer of nanofluids in doubly stratified medium. Although the effect of double stratification on heat and mass transfer in a fluid is significant, very little work has been reported in the literature [9–11]. On the other hand, studies on unsteady boundary layer flow, mass transfer, and convective heat transfer in a nanofluid are the topics of contemporary research areas in fluid science and engineering, owing to its novelty [12–17]. The mass and heat transfer in a doubly strati - fied medium in the presence of a magnetic field is not only fundamentally interesting but also finds wide range of natural, industrial, and engineering applications. These applica - tions include heat rejection into the environment such as lakes, rivers, and seas; thermal energy storage systems such as solar ponds; heat transfer from thermal sources such as condensers of power plants. Despite its numerous applications, it is evident from the literature herein that research has neglected mass and heat transfer in unsteady MHD nanofluid flow taking into account thermophoresis and Brownian motion in a doubly stratified medium. Hence, the motivation of this paper is to bridge this knowledge gap. This study aims to extend the recent work of Olanrewaju and Makinde [18] to incorporate unsteadiness, mag- netic field, and the effects of double stratification on boundary layer nanofluid flow. The shooting iteration technique coupled with fourth-order Runge–Kutta integration scheme is employed to numerically solve the model nonlinear system of equations of Mutuku and Makinde Asia Pac. J. Comput. Engin. (2017) 4:2 Page 3 of 16 this particular problem. Numerical results are displayed graphically and discussed for different values of the various dimensionless parameters controlling the flow regime. It is hoped that the results obtained will not only provide useful information for applications but also serve as a basis for studying other analogous systems that arise in engineering and industrial applications. Model formulation We consider heat and mass transfer in an unsteady flow of an incompressible electrically conducting viscous nanofluid past an infinite nonconducting flat plate in the presence of a transversely imposed unsteady magnetic field of strength B(t). The sketch of the physi - cal configuration and coordinate system is shown in Fig.  1. The induced magnetic field is neglected, which is justified for industrial application MHD flows at small magnetic Reynolds number. It is also assumed that the external electrical field is zero and that the electric field due to the polarization of charges is also negligible. The temperature and nanoparticle concentration at the plate surface and at the free stream are assumed to vary with time and axial distance along the plate. Under the above assumptions with the boundary layer approximations, the continuity, momentum, energy, and nanoparticle concentration equations are as follows: ∂u ∂v + = 0, (1) ∂x ∂y 2 2 ∂u ∂u ∂u ∂U ∂U ∂ u σ B(t) (u − U ) ∞ ∞ ∞ + u + v = + U + υ − , ∞ (2) ∂t ∂x ∂y ∂t ∂x ∂y ρ ∂T ∂T ∂T ∂T ∂T ∂ T υ ∂u ∞ ∞ + u + v = + U + α + ∂t ∂x ∂y ∂t ∂x ∂y c ∂y 2 2 σB(t) (u − U ) ∂C ∂T D ∂T ∞ T − + τ D + , ρc ∂y ∂y T ∂y p 0,∞ (3) Fig. 1 Physical configuration Mutuku and Makinde Asia Pac. J. Comput. Engin. (2017) 4:2 Page 4 of 16 2 2 ∂C ∂C ∂C ∂C ∂C ∂ C D ∂ T ∞ ∞ T + u + v = + U + D + , ∞ B (4) 2 2 ∂t ∂x ∂y ∂t ∂x ∂y T ∂y 0,∞ with the boundary conditions u = 0, v = 0, T = T (x, t), C = C (x, t) at y = 0, w w (5) u → U (x, t), C → C (x, t), T → T (x, t), as y →∞, ∞ ∞ ∞ where u and v are the velocity components along the x and y directions, respectively; t is time; T and C are temperature and nanoparticle volume fraction at the plate surface, w w respectively; U , T , and C are the free stream velocity, temperature, and nanopar- ∞ ∞ ∞ ticle volume fraction, respectively; ρ is the density; α is the thermal diffusivity; υ is the kinematic viscosity; c is the specific heat at constant pressure; σ is the electrical conduc - tivity; τ = (ρC) /(ρC) is the ratio of the effective heat capacitance of the nanoparticle to p f that of the base fluid; D is the Brownian motion coefficient; D is the thermophoretic B T coefficient; T is the constant free stream temperature; and C is the nanoparticle vol- 0,∞ ume fraction. Following Ibrahim and Shankar [19], we prescribe the free stream velocity and the stratified free stream and plate surface temperature and nanoparticles concentration together with the unsteady magnetic field as follows: ax bx cx U (x, t) = , T (x, t) = T + , T (x, t) = T + , ∞ w 0,∞ ∞ 0,∞ 1 − t 1 − t 1 − t (6) B mx nx B(t) = √ , C (x, t) = C + , C (x, t) = C + , w 0,∞ ∞ 0,∞ 1 − t 1 − t 1 − t −1 where λt  <  1, a, b, c, m, and n are positive constants and λ with dimension (time) is unsteadiness frequency parameter. Positive values of b, c, n, and m correspond to the assisting flows and their negative values correspond to opposing flows. In order to sim - plify the mathematical analysis of the problem, we introduce the following dimension- less variables: a aυ η = y, ψ = xf (η), υ(1 − t) (1 − t) T − T T − T cx ∞ 0,∞ θ(η) = = −   , (7) T − T T − T (1 − t) T − T w 0,∞ w 0,∞ w 0,∞ C − C C − C nx ∞ ∞ φ(η) = = − , C − C C − C (1 − t) C − C w 0,∞ w 0,∞ w 0,∞ where η is the similarity variable and ψ is the stream function defined as ∂ψ ∂ψ u = and v =− (8) ∂y ∂x Substituting Eqs. (7) and (8) into Eqs. (1) and (4), we obtain following system of simi- larity equations: ′′′ ′′ ′ ′ ′′ ′ f + ff − f − A f + f − 1 − Ha f − 1 + 1 = 0 (9) 2 Mutuku and Makinde Asia Pac. J. Comput. Engin. (2017) 4:2 Page 5 of 16 ηθ ′′ ′ ′ ′ θ − A Pr + θ + Pr(f θ + β − f θ − f β ) 1 1 (10) 2 2 ′′ ′ 2 ′ ′ ′ + Pr Ec[ f + Ha(f − 1) ]+ N θ φ + N θ = 0 b t ηφ N ′′ ′ ′ ′ ′′ φ − ALe + φ + Le(f φ + β − f φ − f β ) + θ = 0, 2 2 (11) 2 N subject to the boundary conditions f (0) = 0, f (0) = 0, θ(0) = 1 − β , φ(0) = 1 − β , 1 2 (12) (∞) = 1, θ(∞) = 0, φ(∞) = 0, where primes denote differentiation with respect to η and Ha is Hartmann number, Pr is Prandtl number, N is the Brownian motion parameter, N is the thermophoresis param- b t eter, Ec is Eckert number, Le is Lewis number, A is the unsteadiness parameter, β is the thermal stratification parameter, and β is the solutal stratification parameter. The above parameters are defined as follows: σB υ τD C − C τD T − T B w 0,∞ T w 0,∞ Ha = , Pr = , N = , N = , aρ α υ T υ o,∞ (13) U υ  c n Ec = , Le = , A = , β = , β = . 1 2 C T − T D a b m p w 0,∞ The significant physical quantities of engineering and industrial interest in this study are the skin friction coefficient C , Nusselt number Nu, and Sherwood number Sh defined, respectively, as τ xq xq w w m C = , Nu = , Sh = , (14) ρu k(T − T ) D (C − C ) w 0,∞ B w 0,∞ where τ , q , and q are the skin friction, surface heat flux, and the surface mass flux, w w m respectively, defined by ∂u ∂T ∂C τ = μ , q =−k , q =−D . w w m B (15) ∂y ∂y ∂y y=0 y=0 y=0 where µ and k are the dynamic viscosity and the thermal conductivity, respectively. Using Eqs. (15) into (14), we get 1/2 ′′ −1/2 ′ −1/2 ′ Re C = f (0), Re Nu =−θ (0), Re Sh = −φ (0), (16) x x x which are the local skin friction C , local Nusselt number Nu, and local Sherwood num- ber Sh, respectively, and Re  = U x/υ is the local Reynolds number. x ∞ Mutuku and Makinde Asia Pac. J. Comput. Engin. (2017) 4:2 Page 6 of 16 Numerical procedure The system of nonlinear differential Eqs.  (9)–(11) with the boundary conditions (12) has been solved numerically using shooting technique coupled with the fourth-order Runge–Kutta method and a modified version of the Newton–Raphson algorithm. Let ′ ′′ ′ ′ x = f , x = f , x = f , x = θ, x = θ , x = φ, x = φ . 1 2 3 4 5 6 7 (17) Then the set of higher order boundary value equations with their respective boundary conditions is reduced to a system of first-order differential equation with appropriate initial conditions, respectively, x = x , x = x , � � ′ 2 3 x =−x x + x + A x + x − 1 + Ha x − 1 − 1, ( ) 1 3 2 2 3 2 x = x , � � ηx x = A Pr + x − Pr(x x + β − x x − x β ) 4 1 5 1 2 4 2 1 � � �� (18) ′ ′ 2 − Pr Ec x + Ha x − 1 − N x x − N x , b 5 7 t 3 2 5 x = x , � � ηx x = A Le + x + Le(x x + β − x x − x β ) 6 1 7 2 2 6 2 2   � � ηx   A Pr + x − Pr(x x + β − x x − x β N 4 1 5 1 2 4 2 1 − , � � ��   N ′ ′ 2 − Pr Ec x + Ha x − 1 − N x x − N x b 5 7 t 3 2 5 subject to the initial conditions x (0) = 0, x (0) = 0, x (0) = s , x (0) = 1 − β , 1 2 3 1 4 1 (19) x (0) = s , x (0) = 1 − β , x = s , 5 2 6 2 7 3 By applying the shooting method, the unspecified initial conditions S , S , and S in 1 2 3 (19) are assumed and (18) integrated numerically as an initial valued problem to a given terminal point. The accuracy of the assumed missing initial conditions was checked by comparing the calculated value of the dependent variable at the terminal point with its given value there. If differences exist, improved values of the missing initial conditions are obtained and the process repeated. The method is programmed in MAPLE with the −7 step size of ∆η  = 0.001 selected to be satisfactory for a convergence criterion of 10 in nearly all cases. The maximum value of η to each group of parameters was determined when the values of the unknown boundary conditions at η = 0 do not change to success- −7 ful loop with error less than 10 . Results and discussion From the numerical computation process, the effects of the pertinent parameters, namely, unsteadiness parameter A, Lewis number Le, Brownian parameter N , ther- mophoresis parameter N , Hartmann number Ha, Eckert number Ec, thermal strati- fication parameter β , and solutal stratification parameter β on the fluid velocity, 1 2 temperature, nanoparticles volume as well as the skin friction coefficient, the local Mutuku and Makinde Asia Pac. J. Comput. Engin. (2017) 4:2 Page 7 of 16 Nusselt number, and the local Sherwood number are illustrated graphically (Figs.  2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14) and analysed quantitatively highlighting their industrial and engineering applications. In order to ascertain the validity of our numerical proce- dure, the special case of heat transfer in flow of a convectional fluid in the absence of buoyancy force, flow unsteadiness, magnetic field, and double stratification effect (i.e. Fig. 2 Eec ff ts of A and Ha on dimensionless velocity Fig. 3 Eec ff ts of Ec , Nt, and Nb on dimensionless temperature Mutuku and Makinde Asia Pac. J. Comput. Engin. (2017) 4:2 Page 8 of 16 Fig. 4 Eec ff ts of Ha on dimensionless temperature Fig. 5 Eec ff ts of A on dimensionless temperature A = β  = β  = Ha = φ = 0) over a surface can be compared with that of Aman et al. [20], 1 2 as shown in Table  1 and excellent agreement was established. This favourable compari - son lends confidence to the numerical results reported subsequently. Mutuku and Makinde Asia Pac. J. Comput. Engin. (2017) 4:2 Page 9 of 16 Fig. 6 Eec ff ts of β on dimensionless temperature Fig. 7 Eec ff ts of A and Ha on dimensionless nanoparticle concentration Velocity profiles Figure  2 shows the effect of unsteadiness and magnetic field strength on the nanoflu - ids’ velocity profiles. It can be seen that there is a gradual increase in the fluid veloc - ity, from zero at the plate surface to a maximum value as it approaches the free stream far away from the plate surface thus satisfying the boundary conditions. It is observed Mutuku and Makinde Asia Pac. J. Comput. Engin. (2017) 4:2 Page 10 of 16 Fig. 8 Eec ff ts of Le , Nt, and Nb on dimensionless nanoparticle concentration Fig. 9 Eec ff ts of β and Ec on dimensionless nanoparticle concentration that both the fluid velocity and the hydrodynamic boundary layer thickness decrease as the unsteadiness parameter (A) and Hartmann number (Ha) increases. It is well known that application of a transverse magnetic field orthogonal to the flow direction has a ten - dency to create a drag force known as Lorentz force which tends to resist the fluid flow and thus reducing its velocity. Mutuku and Makinde Asia Pac. J. Comput. Engin. (2017) 4:2 Page 11 of 16 Fig. 10 Eec ff ts of Ha and A on local skin friction Fig. 11 Eec ff ts of Ec , Ha, β , Nb, and Nt on local Nusselt number Temperature profiles Figures 3, 4, 5, 6 show the effects of various parameters on the temperature profile. Gen - erally, the temperature is the maximum at the plate surface but decreases exponentially to zero far away from the plate surface satisfying the free stream conditions. As shown in Mutuku and Makinde Asia Pac. J. Comput. Engin. (2017) 4:2 Page 12 of 16 Fig. 12 Eec ff ts of Le and A of local Nusselt number Fig. 13 Eec ff ts of Ha, Ec , β , Nt, and Le of local Sherwood number Fig. 3, an increase in Brownian parameter N , thermophoresis parameter N , and Eckert b t number Ec leads to an increase in both the fluid temperature and the thermal boundary layer thickness. This can be attributed to the additional heating due to viscous dissipa - tion, friction, and rapid collision of the nanoparticles as a result of Brownian motion and thermophoresis. Figure 4 shows that increasing Hartmann number Ha initially increases Mutuku and Makinde Asia Pac. J. Comput. Engin. (2017) 4:2 Page 13 of 16 Fig. 14 Eec ff ts of β and Nb of local Sherwood number Table 1 Computations showing comparison with Aman et al. [20] for different values of  Pr when A = β  = β  = Ha = φ = 0 1 2 ′′ ′ ′′ ′ Pr f (0) −θ (0) f (0) −θ (0) Amani et al. [20] Amani et al. [20] present present 0.1 2.50997 0.26815 2.50997 0.26815 0.3 2.50997 0.40338 2.50997 0.40338 0.7 2.50997 0.52965 2.50997 0.52965 6.2 2.50997 0.59014 2.50997 0.59014 the fluid temperature at the plate surface, but as the fluid moves from the plate surface the temperature decreases. The increase in temperature is due to ohmic heating as a result of the magnetic effect. On the other hand, as illustrated in Figures  5 and 6, there is a decrease in both the temperature and the thermal boundary layer thickness with an increase in thermal stratification parameter β and unsteadiness parameter A. When the thermal stratification effect is taken into consideration, the effective temperature differ - ence between the plate and the ambient fluid decreases; therefore, the temperature is reduced. Nanoparticle concentration profiles Figures 7, 8, 9 show the effects of various parameters on the nanoparticle concentration profile. The nanoparticle concentration is the highest at the plate surface but decreases exponentially to zero free stream value satisfying the boundary conditions. Unlike in the case of temperature, a decrease in the nanoparticle concentration is noted with an increase in unsteadiness parameter A, Lewis number Le, Brownian parameter N , b Mutuku and Makinde Asia Pac. J. Comput. Engin. (2017) 4:2 Page 14 of 16 thermophoresis parameter N , Hartmann number Ha, Eckert number Ec, and solutal stratification parameter β . Eec ff ts of parameter variations on C , Nu, Sh Figures 10, 11, 12, 13, 14 depict the various pertinent parameters at the plate surface for the skin friction coefficient, the local Nusselt number (rate of heat transfer), and the local Sherwood number with variation in the values of thermophysical parameters embed- ded in the problem. From Fig.  10, it can be seen that increasing both the unsteadiness parameter A and Hartmann number Ha leads to an increase in the skin friction coeffi - cient. This is as expected since both parameters causes an overshoot of the fluid towards the plate surface. From Fig.  11, it is noted increasing the Lewis number Le, Brownian parameter N , thermophoresis parameter N , Hartmann number Ha, Eckert number Ec, b t and thermal stratification parameter β leads to a decrease in the local Nusselt number which implies that there is reduced heat transfer rate at the plate surface. As seen earlier, the same parameters led to an increase to the fluid temperature; thus, it is expected that with increased temperature within the fluid the rate of heat transfer from the plate to the fluid will be reduced. On the other hand, an increase in the unsteadiness parameter A increases the local Nusselt number, as observed in Fig.  12. It is noted in Figs.  13 and 14 that increasing the Hartmann number Ha, Eckert number Ec, Lewis number Le, ther- mophoresis parameter N , and thermal stratification parameter β increases the Sher- t 1 wood number, while a decrease in the local Sherwood number is realized by increasing the Brownian parameter N and solutal stratification parameter β . It well known that b 2 at high Lewis number the nanoparticle concentration is low, since the concentration at the plate surface is higher than in the fluid, mass transfer from the plate to the fluid then occurs. Conclusions The combined effects of thermal stratification, solutal stratification, Brownian motion, and thermophoresis on an unsteady hydromagnetic boundary layer flow of a nanofluid over a flat surface are numerically investigated. Using a similarity variable, the governing nonlinear partial differential equations were obtained and transformed to a set of non - linear differential equations. Numerical results are obtained, then presented graphically for the velocity, temperature, and nanoparticle concentration profiles as well as skin fric - tion coefficient, reduced Nusselt number, and reduced Sherwood number. Our results reveal, among other conclusions, the following: • The fluid overshoots towards the plate surface with increasing A and Ha, thus reduc - ing both the fluid velocity and the hydrodynamic boundary layer thickness. • An increase in Ec, Nt, Nb, and Ha yields an increase in both the thermal boundary layer thicknesses and temperature, whereas the reverse is noted with increases in A and β . • A decrease in the nanoparticle concentration is noted with increase in A, Le, Nb, Nt, Ha, Ec, and β . • Increasing A and Ha increase the skin friction. Mutuku and Makinde Asia Pac. J. Comput. Engin. (2017) 4:2 Page 15 of 16 • A decrease in Nusselt number is realized with increased Ec, Ha, Nb, and β , whereas it increases with increasing values of A. • The local Sherwood number increases with increase in Ha, β , Le, Nt, and Ec but reduces with increased β and Nb. Abbreviations (u, v): velocity components in x- and y-direction; A: unsteadiness parameter (λ/a); B : magnetic field strength (T); C: nanoparticle volume fraction; C : skin friction coefficient (τ /U ρ); C : constant free stream volume fraction; D : Brownian diffusion coeffi - w ∞ 0,∞ B cient; D : thermophoresis diffusion coefficient; Ec: Eckert number ( U /C (T  − T )); T ∞ p w 0,∞ 2 2 Ha: magnetic parameter (σB B /αρ); k: thermal conductivity (W/mK); Le: Lewis num- ber (υ/D ); N : Brownian motion parameter (τD (C   −  C )/υ); N : thermophoresis m b B w 0,∞ t parameter (τD (T  − T )/T υ); Pr: Prandtl number (υ/α); Nu: Nusselt number (xq / T w 0,∞ 0,∞ w k(T  − T )); Sh: Sherwood number (xq /D (C  − C )); T: temperature of fluid (K) w 0,∞ m B w 0,∞ Greek symbols ψ: stream function; η: similarity variable; β : thermal stratification parameter; β : solutal 1 2 stratification parameter; υ: kinematic viscosity of fluid (kg/m  s); ρ: density of fluid (kg/ m ); φ: dimensionless concentration function; α: thermal diffusivity(W/m  k); σ: electrical conductivity of fluid (S/m); θ: dimensionless temperature; (c ) : heat capacity of fluid (J/ p f kg  K); (c ) : heat capacity of nanofluid; τ: parameter defined by (ρC) /(ρC) ; µ: dynamic p p p f viscosity of fluid (kg/m s) Subscripts ∞: condition at the free stream; w: condition at the surface Authors’ contributions Both authors were actively involved in the formulation of the problem, literature review, and simulations. Both authors read and approved the final manuscript. Author details 1 2 Department of Mathematics, Kenyatta University, Nairobi, Kenya. Faculty of Military Science, Stellenbosch University, Cape Town, South Africa. Acknowledgements None. Competing interests The authors declare that they have no competing interests. Funding There is no funding towards the publication of this manuscript. Ethics approval and consent to participate This study does not involve human participants, human data, or human tissue. Publisher’s Note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations. Received: 27 March 2017 Accepted: 3 May 2017 References 1. Makinde OD (2011) On MHD mixed convection with Soret and Dufour effects past a vertical plate embedded in a porous medium. LAAR 41(1):63–68 Mutuku and Makinde Asia Pac. J. Comput. Engin. (2017) 4:2 Page 16 of 16 2. Saha LK, Siddiqa S, Hossain MA (2011) Eec ff t of hall current on MHD natural convection flow from vertical perme - able flat plate with uniform surface heat flux. Appl Math Mech Engl 3(9):1127–1146 3. Rout BR, Panda SK, Panda S (2013) MHD heat and mass transfer of chemical reaction fluid flow over a mov- ing vertical plate in presence of heat source with convective surface boundary condition. Int J Chem Eng. doi:10.1155/2013/296834 4. Socolofsky AS, Jirka GH (2004) Stratified flow and buoyant mixing. Eng Lect Notes. Texas A&M University, College Station 5. Singh G, Sharma PR, Chamkha AJ (2010) Eec ff t of thermally stratified ambient fluid on MHD convective flow along a moving non-isothermal vertical plate. Int J Phys Sci 5(3):208–215 6. Saha SC, Hossain MA (2004) Natural convection flow with combined buoyancy effects due to thermal and mass diffusions in a thermally stratified media. NonLinear Anal 9(1):89–102 7. Mukhopadhyay S, Mondal C, Gorla SR (2012) Eec ff ts of thermal stratification on flow and heat transfer past a porous vertical stretching surface. Heat Mass Transf 48:915–921 8. Choi SUS (1995) Enhancing thermal conductivity of fluids with nanoparticle. In: Siginer DA, Wang HP (eds) Devel- opment and applications of non-Newtonian flow. The American Society of Mechanical Engineers, New York, pp 227–233 9. Srinivasacharya D, Upendar M (2013) Eec ff t of double stratification on MHD free convection in a micropolar fluid. J Egypt Math Soc. doi:10.1016/j.joems.2013.02.006 10. Narayana L, Murthy PA (2006) Free convective heat and mass transfer in a doubly stratified non-Darcy porous medium. J Heat Transf 128:1204–1212 11. Ibrahim W, Shankar B (2013) MHD boundary layer flow and heat transfer of a nanofluid past a permeable stretching sheet with velocity, thermal and solutal slip boundary conditions. J Comput Fluids 75:1–10 12. Makinde OD, Olanrewaju PO, Charles WM (2011) Unsteady convection with chemical reaction and radiative heat transfer past a flat porous plate moving through a binary mixture. Afr Mat 22(2011):65–78 13. Chamkha AJ, Rashady AM, Al-Meshaieiz E (2011) Melting effect on unsteady hydromagnetic flow of a nanofluid past a stretching sheet. Int J Chem Reactor Eng 9(113):1 14. Makinde OD (2012) Computational modelling of MHD unsteady flow and heat transfer towards a flat plate with Navier slip and Newtonian heating. Braz J Chem Eng 29(1):159–166 15. Khan MS, Karim I, Ali LE, Islam A (2012) Unsteady MHD free convection boundary layer flow of a nanofluid along a stretching sheet with thermal radiation and viscous dissipation effects. Int Nano Lett 2:24 16. Vajravelu K, Prasad KV, Ng C (2013) Unsteady convective boundary layer flow of a viscous fluid at a vertical surface with variable fluid properties. Nonlinear Anal Real 14:455–464 17. Makinde OD, Chinyoka T (2013) Numerical investigation of buoyancy effects on hydromagnetic unsteady flow through a porous channel with suction/injection. J Mech Sci Tech 27(5):1557–1568 18. Olanrewaju AM, Makinde OD (2013) On boundary layer stagnation point flow of a nanofluid over a permeable flat surface with Newtonian heating. Chem Eng Commun 200:836–852 19. Ibrahim W, Shanker B (2012) Unsteady MHD boundary layer flow and heat transfer due to stretching sheet in the presence of heat source or sink by quasi linearization technique. Int J Appl Math Mech 8(7):18–30 20. Aman F, Ishak A, Pop I (2011) Mixed convection boundary layer flow near stagnation point on vertical surface with slip. Appl Math Mech 32(12):1599–1606 Submit your next manuscript to BioMed Central and take full advantage of: • Convenient online submission • Thorough peer review • No space constraints or color figure charges • Immediate publication on acceptance • Inclusion in PubMed, CAS, Scopus and Google Scholar • Research which is freely available for redistribution Submit your manuscript at www.biomedcentral.com/submit

Journal

Asia Pacific Journal on Computational EngineeringSpringer Journals

Published: May 12, 2017

References