Unsteady Mixed Convection Flow along Symmetric Wedge with Variable Surface Temperature Embedded in a Porous Medium Saturated with a Nanofluid

Laminar two-dimensional unsteady mixed-convection boundary-layer flow of a viscous incompressible fluid past asymmetric wedge with variable surface temperature embedded in a porous medium saturated with a nanofluid has been studied. The employed mathematical model for the nanofluid takes into account the effects of Brownian motion and thermophoresis. The velocity in the potential flow is assumed to vary arbitrary with time. The non-Darcy effects including convective, boundary and inertial effects will be included in the analysis. The unsteadiness is due to the time-dependent free stream velocity. The governing boundary layer equations along with the boundary conditions are converted into dimensionless form by a non-similar transformation, and then resulting system of coupled non-linear partial differential equations are solved by perturbation solutions for small dimensionless time until the second order. Numerical solutions of the governing equations are obtained employing the implicit finite-difference scheme in combination with the quasi-linearization technique. To validating the method used, we compared our results with previous results in earlier papers on special cases of the problem and are found to be in agreement. Effects of various parameters on velocity, temperature and nanoparticle volume fraction profiles are graphically presented.

• Nanofluids
• Unsteady Flow
• Mixed Convection
• Boundary Layer
• Wedge Flow
• Finite Difference
• Non-Similarity Solutions
• Porous Media

1. Introduction

The study of mixed convection flow finds applications in several industrial and technical processes such as nuclear reactors cooled during emergency shutdown, solar central receivers exposed to winds, electronic devices cooled by fans, heat exchanges placed in a low-velocity environment, etc. The mixed convection flows become important when the buoyancy forces due to the temperature difference between the wall and the free stream becomes large.

The study and analysis of heat and mass transfer in porous media has been the subject of many investigations due to their frequent occurrence in industrial and technological applications. Examples of some applications include geothermal reservoirs, drying of porous solids, thermal insulations and many others. Smith [1] initiated the study of the unsteady incompressible forced convection boundary-layer flow past a semi-infinite wedge impulsively set into motion. This problem subsequently solved numerically by Nanbu [2] using the method proposed by Hall [3] and then that modified by Harris et al. [4]. This method solves the untransformed equations directly using an iterative procedure and by implicit finite difference method, which is well documented and widely used by Keller and Cebeci [5] and also by Hossain [6] [7] and Hossain et al. [8] for unsteady boundary layers. Watkins [9] has solved this problem numerically following a second order; he has also studied the unsteady heat transfer aspects of the semi-infinite wedge started impulsively from rest to include solutions of the energy equation. A new set of scaled coordinates introduced by Williams and Rhyne [10].

There is a large body of literature on unsteady, mixed convection, boundary- layer flows past bodies of different geometries. Ibrahim et al. [11] studied the problem of unsteady magneto hydrodynamic (MHD) micropolar fluid flow and heat transfer over a vertical porous plate through a porous medium in the presence of thermal and mass diffusion with constant heat source. The unsteady free convection flow in the stagnation-point region of a rotating sphere embedded in a porous medium has analyzed by Hassanien et al. [12]. Ibrahim [13] has investigated the problem of unsteady mixed conviction flow in the stagnation region of a three dimensional body embedded in a porous medium. Al-Harbi and Hassanien [14] have studied the problem of unsteady mixed convection flow in the stagnation region of a heated vertical plate embedded in a variable porosity medium with thermal dispersion effects. Also, Al-Harbi [15] has analyzed the problem of unsteady mixed convection flow on a rotating cone embedded in a porous medium saturated with a rotating fluid. Unsteady magneto hydrodynamic mixed convection flow of an electrically conducting nanofluid in a stagnation region of a rotating sphere is studied numerically through Sameh and Rashed [16].

The interested reader can find an excellent collection of papers on unsteady convective flow problems over heated bodies embedded in a fluid-saturated porous medium in the book by Pop and Ingham [17] and in the book by Nield and Bejan [18].

Nanofluids with enhanced thermal characteristics have widely been examined to improve the heat transfer performance of many engineering applications [19]. Most researchers argue that the addition of nanoparticles with relatively higher thermal conductivity to the base fluid results in an increase of the thermal performance of the resultant nanofluid [20] [21] [22] [23].

According to Yacob et al. [24], nanofluids are produced by dispersing the nanometer-scale solid particles into base liquids with low thermal conductivity such as water and ethylene glycol. Nanoparticles are usually made of metal, metal oxide, carbide, nitride and even immiscible nano-scale liquid droplets. Gorla et al. [25] presented a boundary layer analysis for the mixed convection past a vertical wedge in a porous medium saturated with a nanofluid. Hamad et al. [26] introduced a one parameter group to represent similarity reductions for the problem of magnetic field effects on free-convective nanofluid flow past a semi-infinite vertical flat plate following a nanofluid model proposed by Buongiorno [27]. Hamad [28] obtained the analytical solutions for convective flow and heat transfer of a viscous incompressible nanofluid past a semi-infinite vertical stretching sheet in the presence of magnetic field. Khan and Pop [29] obtained similar solutions depending on Prandtl number, Lewis number, Brownian motion number and thermophoresis number on the steady boundary layer flow, heat and mass over a stretching surface in its plane. Further, Abu-Nada and Chamkha [30] have presented the natural convection heat transfer characteristics in a differentially heated enclosure filled with a CuO-EG-water nanofluid for different variable thermal conductivity and variable viscosity models. For more information, see also Das et al. [31], and Kakaç and Pramuanjaroenkij [32]. Muthtamilselvan et al. [33] claimed that it is difficult to have a precise idea on how nanoparticle enhances the heat transfer characteristics of nanofluids. More recently, Hady et al. [34] studied a nonsimilar solution of steady forced convection boundary layer flow and heat transfer of a nanofluid past a stretching horizontal late sing one-phase model. On the other hand, the same authors [35] extended the problem of Cheng-Min- kowycz problem for natural convective boundary-layer flow in a porous medium saturated by a nanofluid to study the effect of the porosity of the plate. In addition, Hady et al. [36] studied the problem of the flow and heat transfer characteristics of a viscous nanofluid over a nonlinearly stretching sheet in the presence of thermal radiation and variable wall temperature. James et al. [37] have studied the effects of variable viscosity of nanofluid flow over a permeable wedge embedded in saturated non-Darcy porous medium with chemical reaction and thermal radiation. In existence of concerning magnetic field, heat together with mass transfer features on mixed convective copper-water nanofluid flow through inclined plate is investigated in surrounding porous medium together with viscous dissipation through Nasir et al. [38]. Chamkha et al. [39] introduced a mathematical model to accentuated the mixed bioconvective flow on a vertical wedge in a Darcy porous medium filled with a nanofluid containing both nanoparticles and gyrotactic microorganisms.

The aim of the present paper is to study the unsteady mixed convection flow along a symmetric wedge embedded in a porous medium saturated with a nanofluid in the presence of first and second orders resistances, which to the best of our knowledge have not been investigated yet. Motivation to study mixed convection in porous media comes from the need to characterize the convective transport processes around deep geological repository for the disposal of high- level nuclear waste, e.g. spent fuel rods from nuclear reactors (see Lai [40] ).

The effect of the presence of the buoyancy forces and the isotropic solid matrix on the unsteady mixed convection flow along a symmetric wedge embedded in a porous medium saturated with a nanofluid are considered. In addition, the Brownian motion and the thermophoresis effects are considered. The unsteadiness in the flow field is caused by impulsively creating motion in the free stream and at the same time suddenly raising the surface temperature above its surroundings. The partial differential equations governing the flow and the heat transfer solved numerically using some different numerical methods as the finite difference scheme by Pereyra [41]. Particular cases of the present results compared with those of Scshadri et al. [42].

2. Governing Equations

Let us consider an unsteady mixed convection boundary layer flow of an incompressible fluid past along a symmetric wedge with variable surface temperature embedded in a porous medium saturated with a nanofluid as shown in Figure 1. The inviscid flow over the wedge develops instantaneously and its velocity is given by

u ¯ e ( x ¯ ) = U ∞ ( x ¯ L ) m           for     m ≤ 1 (1)

where L is a characteristic length and m is pressure gradient related in the included angle π β by m = β / ( 2 − β ) . It is clear that for negative values of m the solution becomes singular at x ¯ = 0 . While for m positive the solution can be defined for all values of x ¯ . It is assumed that the variable surface temperature of the wedge is ( T w > T ∞ ) where T ∞ is the ambient temperature of the fluid and

T ¯ is the temperature of the fluid. Under the above-mentioned assumptions, the boundary layer equations governing the flow can be expressed as follows:-

∂ u ¯ ∂ x ¯ + ∂ v ¯ ∂ y ¯ = 0 , (2)

∂ u ¯ ∂ t ¯ + u ¯ ∂ u ¯ ∂ x ¯ + v ¯ ∂ u ¯ ∂ y ¯ = u ¯ e d u ¯ e d x ¯ + υ ∂ 2 u ¯ ∂ y ¯ 2 + υ ε K ( u ¯ e − u ¯ ) + Γ ε 2 K 1 / 2 ( u ¯ e 2 − u ¯ 2 )       + [ ( 1 − φ ∞ ) β T ρ f ( T ¯ − T ∞ ) − ( ρ p − ρ f ) ( φ ¯ − φ ∞ ) / ρ f ] g (3)

∂ T ¯ ∂ t ¯ + u ¯ ∂ T ¯ ∂ x ¯ + v ¯ ∂ T ¯ ∂ y ¯ = α m ∂ 2 T ¯ ∂ y ¯ 2 + τ [ D B ∂ T ¯ ∂ y ¯ ∂ φ ¯ ∂ y ¯ + D T T ∞ ( ∂ T ¯ ∂ y ) 2 ] , (4)

∂ φ ¯ ∂ t ¯ + u ¯ ∂ φ ¯ ∂ x ¯ + v ¯ ∂ φ ¯ ∂ y ¯ = D B ∂ 2 φ ¯ ∂ y ¯ 2 + D T T ∞ ( ∂ 2 T ¯ ∂ y ¯ 2 ) , (5)

where u ¯ and v ¯ are the velocity components in is x ¯ and y ¯ coordinate, respectively. T ¯ is the temperature, φ ¯ is the nanoparticle volume fraction, g is the acceleration due to gravity, ρ f is the density of the base fluid and ν and β T are the viscosity and thermal conductivity of the fluid. While, ρ p is the density of the nanoparticles, ( ρ c ) f is the heat capacity of the fluid and ( ρ c ) p is the effective heat capacity of the nanoparticle material. α m = k m / ( ρ c ) f is the thermal diffusivity, τ = ( ρ c ) p / ( ρ c ) f is the ratio of the effective heat capacity of the nanoparticle material to the heat capacity of the fluid and k m is the effective thermal conductivity of the porous medium. K, ε and Γ are the permeability, the porosity of the porous medium and the empirical constant in the second- order resistance. The coefficients that appear in Equations ((4), (5)) are the Brownian diffusion coefficient D B , and the thermophoretic diffusion coefficient D T .

The boundary conditions for the present problem are

u ¯ = 0 ,           v ¯ = 0 ,         T ¯ = T w ( x ¯ ) = T ∞ + ( T r − T ∞ ) ( x ¯ L ) 2 m − 1 , φ ¯ = φ w ( x ¯ ) = φ ∞ + ( φ r − φ ∞ ) ( x ¯ L ) 2 m − 1 ,         at         y ¯ = 0 , u ¯ = u ¯ e ( x ¯ ) ,             T ¯ → T ∞ ,             φ ¯ → φ ∞                   as       y ¯ → ∞ . (6)

We introduce non dimensional dependent and independent variables according to,

x = x ¯ L ,             y = R e L 1 / 2 y ¯ L ,             u = u ¯ U ∞ ,             v = R e L 1 / 2 v ¯ U ∞ ,             u e = u ¯ e U ∞ , t = U ∞ L t ¯ ,               T = T ¯ − T ∞ T r − T ∞ ,           φ = φ ¯ − φ ∞ φ r − φ ∞ . (7)

The velocity over the wedge is now given by u e ( x ) = x m for m ≤ 1 .

The governing Equations (2)-(5) can be written as

∂ u ∂ x + ∂ v ∂ y = 0 , (8)

∂ u ∂ t + u ∂ u ∂ x + v ∂ v ∂ y = u e d u e d x + ∂ 2 u ∂ y 2 + ε D a L R e L ( u e − u ) + Γ ε 2 D a L 1 / 2 ( u e 2 − u 2 ) + λ ( T − N r φ ) , (9)

∂ T ∂ t + u ∂ T ∂ x + v ∂ T ∂ y = 1 P r ∂ 2 T ∂ y 2 + N o B ∂ T ∂ y ∂ φ ∂ y + N o T ( ∂ T ∂ y ) 2 , (10)

∂ φ ∂ t + u ∂ φ ∂ x + v ∂ φ ∂ y = 1 S c ∂ 2 φ ∂ y 2 + 1 S c N o T N o B ∂ 2 T ∂ y 2 , (11)

where D a L = K / L 2 is the Darcy number, λ = G r L / R e L 2 is the mixed convection parameter, R e L = U ∞ L / ν is the Reynolds number, G r L = g β ( T r − T ∞ ) L 3 / ν 2 is the Grashof number, N r = ( ρ p − ρ f ) ( φ r − φ ∞ ) ρ f ( 1 − φ ∞ ) ( T r − T ∞ ) is the nanofluid buoyancy ratio parameter, N o T = τ ν ( D T T ∞ ) ( T r − T ∞ ) is the thermophores is parameter, N o B = τ ν D B ( φ r − φ ∞ ) is the Brownian motion parameter, P r = α m / ν is the Pran- dtl number and S c = ν / B B is the Schmidt number.

For t ≥ 0 and m ≤ 1 the boundary conditions then may be written as:

u = 0 ,         v = 0 ,           T = x 2 m − 1 ,       φ = x 2 m − 1       at   y = 0 , u = u e ( x ) = x m ,           T → 0 ,       φ → 0     as     y → ∞ (12)

The number of independent variables in the governing Equations (8)-(11) can be reduced from three to two by applying the following transformations as

η = x ( m − 1 ) / 2 ( 2 ξ ) − 1 / 2   y ,           ψ = x ( m + 1 ) / 2 ( 2 ξ ) 1 / 2 f ( η , ξ ) , T = x 2 m − 1 θ ( η , ξ ) ,         φ = x 2 m − 1 ϕ ( η , ξ ) ,             ξ = 1 − e − τ ,       τ = x m − 1 t , (13)

where η is a non-dimensional similarity variable and ψ is the stream function, which is defined in the usual way, namely u = ∂ ψ / ∂ y and v = − ∂ ψ / ∂ y .

Substituting the transformations (13) into Equations (9)-(11), we obtained the following transformed equations for the momentum and thermal boundary layer equations:

f ‴ + η ( 1 − ξ ) f ″ + [ ( m + 1 ) ξ − ( m − 1 ) ( 1 − ξ ) ln ( 1 − ξ ) ] f f ″   + 2 ξ [ γ x ( 1 − f ′ ) + ( m + Δ x ) ( 1 − f ′ 2 ) ] = 2 ( 1 − m ) ξ ( 1 − ξ ) ln ( 1 − ξ ) ( f ′ ∂ f ′ ∂ ξ − f ″ ∂ f ∂ ξ )       + 2 ξ ( 1 − ξ ) ∂ f ′ ∂ ξ − 2 λ ξ ( θ − N r ϕ ) , (14)

1 P r θ ″ + η ( 1 − ξ ) θ ′ + [ ( m + 1 ) ξ − ( m − 1 ) ( 1 − ξ ) ln ( 1 − ξ ) ] f θ ′   − 2 ( 2 m − 1 ) ξ θ f ′ = 2 ( 1 − m ) ξ ( 1 − ξ ) ln ( 1 − ξ ) ( f ′ ∂ θ ∂ ξ − θ ′ ∂ f ∂ ξ )       + 2 ξ ( 1 − ξ ) ∂ θ ∂ ξ − N B θ ′ ϕ ′ − N T θ ′ 2 (15)

1 S c ϕ ″ + η ( 1 − ξ ) ϕ ′ + [ ( m + 1 ) ξ − ( m − 1 ) ( 1 − ξ ) ln ( 1 − ξ ) ] f ϕ ′   − 2 ( 2 m − 1 ) ξ ϕ f ′ = 2 ( 1 − m ) ξ ( 1 − ξ ) ln ( 1 − ξ ) ( f ′ ∂ ϕ ∂ ξ − ϕ ′ ∂ f ∂ ξ ) + 2 ξ ( 1 − ξ ) ∂ ϕ ∂ ξ − N T S c N B θ ″ (16)

where γ x = D a − 1 R e L − 1 x 1 − m is the local permeability parameter, Δ x = F x is the inertia coefficient parameter (see Chamkha [23] ). The modified thermophoresis and Brownian motion parameters are defined as:

N T = N o T x 2 m − 1 = τ ν ( D T T ∞ ) ( T r − T ∞ ) x 2 m − 1 = τ ν ( D T T ∞ ) ( T w − T ∞ ) , N B = N o B x 2 m − 1 = τ ν D B ( φ r − φ ∞ ) x 2 m − 1 = τ ν D B ( φ w − φ ∞ ) (17)

The boundary conditions to be satisfied by the Equations (14)-(16) are

f ( 0 , ξ ) = f ′ ( 0 , ξ ) = 0 ,           θ ( 0 , ξ ) = 1 ,         ϕ ( 0 , ξ ) = 1 , f ′ ( ∞ , ξ ) = 1 ,           θ ( ∞ , ξ ) = 0 ,           ϕ ( ∞ , ξ ) = 0 , (18)

In the above equations, prime denotes differentiation of the functions with respect to η only.

Now to find the numerical solutions we can get the easiest form by using the transformations ξ = 1 − e − τ from Equations (14)-(16) applicable for 0 ≤ τ ≤ ∞ .

f ‴ + η e − τ f ″ + [ ( m + 1 ) ( 1 − e − τ ) + ( m − 1 ) τ e − τ ] f f ″   + 2 ( 1 − e − τ ) [ γ x ( 1 − f ′ ) + ( m + Δ x ) ( 1 − f ′ 2 ) ] = 2 ( m − 1 ) τ ( 1 − e − τ ) ( f ′ ∂ f ′ ∂ τ − f ″ ∂ f ∂ τ ) + 2 ( 1 − e − τ ) ( ∂ f ′ ∂ τ )       − 2 λ ( 1 − e − τ ) ( θ − N r ϕ ) (19)

1 P r θ ″ + η e − τ θ ′ + [ ( m + 1 ) ( 1 − e − τ ) + ( m − 1 ) τ e − τ ] f θ ′   − 2 ( 2 m − 1 ) ( 1 − e − τ ) θ f ′ = 2 ( m − 1 ) τ ( 1 − e − τ ) ( f ′ ∂ θ ∂ τ − θ ′ ∂ f ∂ τ ) + 2 ( 1 − e − τ ) ∂ θ ∂ τ − N B θ ′ ϕ ′ − N T θ ′ 2 (20)

1 S c ϕ ″ + η e − τ ϕ ′ + [ ( m + 1 ) ( 1 − e − τ ) + ( m − 1 ) τ e − τ ] f ϕ ′   − 2 ( 2 m − 1 ) ( 1 − e − τ ) ϕ f ′ = 2 ( m − 1 ) τ ( 1 − e − τ ) ( f ′ ∂ ϕ ∂ τ − ϕ ′ ∂ f ∂ τ ) + 2 ( 1 − e − τ ) ∂ ϕ ∂ τ − N T S c N B θ ″ (21)

The boundary conditions to be satisfied by the above equations are

f ( 0 , τ ) = f ′ ( 0 , τ ) = 0 ,       θ ( 0 , τ ) = 1 ,         ϕ ( 0 , τ ) = 1 , f ′ ( ∞ , τ ) = 1 ,       θ ( ∞ , τ ) = 0 ,         θ ( ∞ , τ ) = 0. (22)

In practical applications, two quantities of physical interest are to be determined, such as, surface shear stress and the rate of heat and mass transfer at the surface. These may obtained in terms of the skin friction coefficient (wall shear stress) C f , local Nusselt number N u x and the local Sherwood number S h x , which are defined by:

C f = μ ρ f U e 2 ( x ¯ ) ( ∂ u ¯ ∂ y ¯ ) y = 0 ,         N u x = x ¯ T w − T ∞ ( ∂ T ¯ ∂ y ¯ ) y = 0 , S h x = x ¯ φ w − φ ∞ ( ∂ φ ¯ ∂ y ¯ ) y = 0 (23)

By introducing the non-dimensional variables (7) and the transformation (13), the skin friction coefficient, C f , the local Nusselt number, N u x and the local Sherwood number, S h x can now defined by:

C f R e x 1 / 2 = [ 2 ( 1 − e − τ ) ] − 1 / 2 f ″ ( 0 , τ ) , N u x R e x − 1 / 2 = − [ 2 ( 1 − e − τ ) ] − 1 / 2 θ ′ ( 0 , τ ) , S h x R e x − 1 / 2 = − [ 2 ( 1 − e − τ ) ] − 1 / 2 ϕ ′ ( 0 , τ ) . (24)

Perturbation solutions method for small time τ ≪ 1

For small τ Equations (19)-(21) become

f ‴ + η f ″ + 2 m τ f f ″ + 2 τ [ γ x ( 1 − f ′ ) + ( m + Δ x ) ( 1 − f ′ 2 ) + λ ( θ − N r ϕ ) ] = 2 ( m − 1 ) τ 2 [ f ′ ∂ f ′ ∂ τ − f ″ ∂ f ∂ τ ] + 2 τ [ ∂ f ′ ∂ τ ] (25)

1 P r θ ″ + η θ ′ + 2 m τ f θ ′ − 2 ( 2 m − 1 ) τ θ f ′ + N B θ ′ ϕ ′ + N T θ ′ 2 = 2 ( 1 − m ) τ 2 [ f ′ ∂ θ ∂ τ − θ ′ ∂ f ∂ τ ] + 2 τ ∂ θ ∂ τ (26)

1 S c ϕ ″ + η ϕ ′ + 2 m τ f ϕ ′ − 2 ( 2 m − 1 ) τ ϕ f ′ + N T S c N B θ ″ = 2 ( 1 − m ) τ 2 [ f ′ ∂ ϕ ∂ τ − ϕ ′ ∂ f ∂ τ ] + 2 τ ∂ ϕ ∂ τ (27)

and the corresponding boundary conditions (22) become

f ( 0 , τ ) = f ′ ( 0 , τ ) = 0 ,       θ ( 0 , τ ) = 1 ,       ϕ ( 0 , τ ) = 1 , f ′ ( ∞ , τ ) = 1 ,       θ ( ∞ , τ ) = 0 ,         ϕ ( ∞ , τ ) = 0 (28)

The resulting system of Equations (25)-(27), along with boundary conditions (28) is solved using the perturbation technique. Now, the non-dimensional stream function and the temperature functions can be written as:

f ( η , τ ) = f 0 ( η ) + τ f 1 ( η ) + τ 2 f 2 ( η ) + ⋯ (29)

θ ( η , τ ) = θ 0 ( η ) + τ θ 1 ( η ) + τ 2 θ 2 ( η ) + ⋯ (30)

ϕ ( η , τ ) = ϕ 0 ( η ) + τ ϕ 1 ( η ) + τ 2 ϕ 2 ( η ) + ⋯ (31)

Substituting Equations (29)-(31) into Equations (25)-(27) and equating the various coefficients of power of τ to zero (here we collect terms up to the second power of O ( τ 2 ) , we can get the following sets of ordinary differential equations:

Zero order:

f ‴ 0 + η f ″ 0 = 0 (32)

1 P r θ ″ 0 + η θ ′ 0 + N B θ ′ 0 ϕ ′ 0 + N T θ ′ 0 2 = 0 , (33)

1 S c ϕ ″ 0 + η ϕ ′ 0 + N T S c N B θ ″ 0 = 0 , (34)

with the corresponding boundary conditions:

f 0 ( 0 ) = f ′ 0 ( 0 ) = 0 ,         θ 0 ( 0 ) = ϕ 0 ( 0 ) = 1 ,         f ′ 0 ( ∞ ) = 1 ,         θ 0 ( ∞ ) = ϕ 0 ( ∞ ) = 0 (35)

First order:

f ‴ 1 + η f ″ 1 − 2 f ′ 1 = − 2 m ( 1 − f ′ 0 2 + f 0 f ″ 0 ) + 2 ( γ x f ′ 0 + Δ x f ″ 0 2 ) − 2 λ ( θ 0 − N r ϕ 0 ) , (36)

1 P r θ ″ 1 + η θ ′ 1 − 2 θ 1 + N B θ ′ 1 ϕ ′ 1 + N T θ ′ 1 2 = − 2 m f 0 θ ′ 0 + 2 ( 2 m − 1 ) θ 0 f ′ 0 , (37)

1 S c ϕ ″ 1 + η ϕ ′ 1 − 2 ϕ 1 + N T S c N B θ ″ 1 = − 2 m f 0 ϕ ′ 0 + 2 ( 2 m − 1 ) ϕ 0 f ′ 0 , (38)

with the corresponding boundary conditions:

f 1 ( 0 ) = f ′ 1 ( 0 ) = 0 ,         θ 1 ( 0 ) = ϕ 1 ( 0 ) = 0 ,         f ′ 1 ( ∞ ) = 0 ,         θ 1 ( ∞ ) = ϕ 1 ( ∞ ) = 0 (39)

Second order:

f ‴ 2 + η f ″ 2 − 4 f ′ 2 = 2 [ ( 1 − 2 m ) f 1 f ″ 0 + ( 3 m − 1 ) f ′ 1 f ′ 0 − m f 0 f ″ 1 − λ ( θ 1 − N r ϕ 1 ) ] + 2 γ x f ′ 1 + 4 Δ x f ′ 0 f ′ 1 (40)

1 P r θ ″ 2 + η θ ′ 2 − 4 θ 2 + N B θ ′ 2 ϕ ′ 2 + N T θ ′ 2 2 = 2 [ ( 1 − 2 m ) f 1 θ ′ 0 + ( 3 m − 2 ) θ 1 f ′ 0 − m θ ′ 1 f 0 + ( 2 m − 1 ) θ 0 f ′ 1 ] (41)

1 S c ϕ ″ 2 + η ϕ ′ 2 − 4 ϕ 2 + N T S c N B θ ″ 2 = 2 [ ( 1 − 2 m ) f 1 ϕ ′ 0 + ( 3 m − 2 ) ϕ 1 f ′ 0 − m ϕ ′ 1 f 0 + ( 2 m − 1 ) ϕ 0 f ′ 1 ] (42)

With the corresponding boundary conditions:

f 2 ( 0 ) = f ′ 2 ( 0 ) = 0 ,       θ 2 ( 0 ) = 0 ,       θ 2 ( 0 ) = 0 ,       f ′ 2 ( ∞ ) = 0 ,       θ 2 ( ∞ ) = 0 ,       ϕ 2 ( ∞ ) = 0 (43)

where primes denote differentiation with respect to η . Knowing the values of f ″ 0 ( 0 ) , f ″ 1 ( 0 ) , f ″ 2 ( 0 ) , θ ′ 0 ( 0 ) , θ ′ 1 ( 0 ) , θ ′ 2 ( 0 ) , ϕ ′ 0 ( 0 ) , ϕ ′ 1 ( 0 ) and ϕ ′ 2 ( 0 ) from the solutions of Equations (32)-(43), we get the values of the skin friction coefficient C f , local Nusselt number N u x and the local Sherwood number S h x , from the following expressions:

C f R e x 1 / 2 = ( 2 τ ) − 1 / 2 ( f ″ 0 ( 0 ) + τ f ″ 1 ( 0 ) + τ 2 f ″ 2 ( 0 ) ) , N u x R e x − 1 / 2 = − ( 2 τ ) − 1 / 2 ( θ ′ 0 ( 0 ) + τ θ ′ 1 ( 0 ) + τ 2 θ ′ 2 ( 0 ) ) , S h x R e x − 1 / 2 = ( 2 τ ) − 1 / 2 ( ϕ ′ 0 ( 0 ) +   τ ϕ ′ 1 ( 0 ) + τ 2 ϕ ′ 2 ( 0 ) ) . (44)

3. Results and Discussion

The sets of ordinary differential Equations (32)-(34) are solved successively by giving appropriate initial guess values for f ″ i ( 0 ) , θ ′ i ( 0 ) , ϕ ′ i ( 0 ) , i = 0 , 1 , 2 to match the values with the corresponding boundary conditions at f ′ i ( ∞ ) , θ i ( ∞ ) , ϕ i ( ∞ ) , i = 0 , 1 , 2 . The numerical values of the coefficients of skin friction, C f R e x 1 / 2 and the Nusselt number N u x R e x − 1 / 2 for smaller values of τ ∈ [ 0 , 1 ] obtained by perturbation method while P r = 1.0 and m = 0.2 are shown in Table 1 in the absence of the effect of the buoyancy force ( λ = 0 ), thermophoresis parameter ( N T = 0 ) and Brownian motion parameter ( N B = 0 ). These values compared with that of Hossain et al. [8] and Harris et al. [4]. The comparison shows excellent agreement between these two results.

Numerical results presented for some representative values of the governing parameters govern this problem. In order to see the physical insight, the numerical values of velocity f ′ ( η ) , temperature θ ( η ) , and nanoparticle volume fraction ϕ ( η ) within the boundary layer computed for different parameters as unsteadiness parameter τ , mixed convection parameter λ , nanofluid buoyancy ratio parameter N_r, thermophoresis parameter N_T, Brownian motion parameter N_B, first resistant parameter γ , second resistant parameter Δ , Prandtl number Pr and Schmidt number Sc. In addition, numerical results obtained to discuss the effects of the governing parameters on the skin friction coefficient, local Nussel number and local Shrewood number and displayed in tabular and graphical forms.

Figures 2-4 displayed the effects of nanonfluid parameters, thermophoresis parameter N_T, Brownian motion parameter N_B and buoyancy ratio parameter N_r. From these figures, it is observed that the velocity and the temperature profiles

Table 1

Comparison of the present numerical results of C f R e x 1 / 2 and N u x R e x − 1 / 2 for different values of small time while Pr = 1.0 and m = 0.2

τ	Harris et al. [4]		Hossian et al. [11]		Present result
	C f R e x 1 / 2	N u x R e x − 1 / 2	C f R e x 1 / 2	N u x R e x − 1 / 2	C f R e x 1 / 2	N u x R e x − 1 / 2
0.01	5.65797	5.64360	5.66886	5.64525	5.65769	5.63468
0.10	1.83491	1.78946	1.83693	1.76140	1.83312	1.75951
0.20	1.33334	1.26902	1.33189	1.22357	1.33090	1.22543
0.40	0.99341	0.90234	0.99548	0.83937	0.99038	0.83739
1.0	0.72370	0.57925	0.72665	0.46370	0.72116	0.46104

increase with an increase in both N_T and N_B. However, for the nanoparticle volume fraction profiles there is a crossing over point where the volume fraction profile decreases before that point and slightly increases after that (Figure 2 and Figure 3). Hence, the heat transfer enhancement is due to collision of high thermal energy with lower energy particles. From Figure 4, it is observed that with the increase in the buoyancy ratio parameter N_r, the velocity profile decreases and both the temperature and volume fraction profiles slightly increase.

The effects of the mixed convection parameter λ on the non-dimensional velocity, temperature and nanoparticle volume fraction are illustrated in Figure 5.

From this figure, it is observed that the velocity increases as λ increases, however; both the temperature and nanoparticle fraction profiles slightly decrease with increasing values of λ .

Figure 6 & Figure 7 exhibit the effects of first resistant parameter γ and second resistant parameter Δ of the porous medium on velocity, temperature

and nanoparticle volume fraction. An increase in γ or Δ leads to an increase in fluid velocity and a decrease on both temperature and nanoparticle volume fraction. This is due the fact that an increase in γ or Δ implies that there is a decrease in the resistance of the porous medium which tends to accelerate the fluid velocity in the boundary layer region.

Figure 8 displays the effect of Prandtl number Pr on velocity, temperature and nanoparticle volume fraction profiles. It is evident from the figure that, as Pr increases, velocity profiles have a small increasing. On the other hand, for both temperature and nanoparticle volume fraction profiles, there special points called the “crossing over points” and the temperature and the volume fraction profiles have completely conflicting behavior before and after these points for which the profiles at a fixed value of η decreases before that point and increases after it. Prandtl number signifies the ratio of momentum diffusivity to thermal diffusivity. In heat transfer problems, the Prandtl number Pr controls the relative thickening of the momentum and thermal boundary layers. Hence, Prandtl number is small, heat diffuses quickly compared to the velocity (momentum), which means that for liquid metals, the thickness of the thermal boundary layer is much bigger than the momentum boundary layer. Fluids with lower Prandtl number have higher thermal conductivities (and thicker thermal boundary layer structures) so that heat can diffuse from the sheet faster than for higherPr fluids (thinner boundary layer). So, Prandtl number can be used to increase the rate of cooling in conducting flows.

Figure 9 depicts the effect of Schimdt number Sc on velocity, temperature and species concentration profiles. It is noticed that effect of increasing value ofSc is

to decrease velocity, temperature and species concentration profiles. This is consistent with the fact that, increase in Sc means decrease of molecular diffusivity those results in decease of concentration boundary layer. Hence species concentration is higher for small values of Sc and lower for large value of Sc.

Figures 10-17 are presented to illustrate the variations of the local rate of shear stress, the local rate of heat transfer and the local rate of mass transfer for different values of the governing parameters.

Figure 10 display the effects of the thermophoresis parameter N_T on the local skin friction coefficient C f , local Nusselt number N u x and the local Sherwood number S h x . It can be seen that the thermophoresis parameter N_T appears in the thermal and the concentration boundary layer equations. As we note that, it is coupled with the temperature function and plays a strong role in determining the diffusion of heat and nanoparticles concentration in the boundary layer. Thus, an increase in the value of the thermophoresis parameter N_T leads to a decrease in the local Nusselt number and increases in both the local skin friction coefficient and the local Sherwood number.

Figure 11 display the effects of the Brownian motion parameter N_B on the local skin friction coefficient, local Nusselt number and the local Sherwood number. In nanofluid systems, owing to the size of the nanoparticles, Brownian motion takes place, and this can enhance the heat transfer properties. This is due to the fact that the Brownian diffusion promotes heat conduction. The nanoparticles increase the surface area for heat transfer. A nanofluid is a two-phase fluid where the nanoparticles move randomly and increase the energy exchange rates. However, the Brownian motion reduces nanoparticles diffusion. The increase in the local Sherwood number as N_B changes is relatively small. Therefore, as indicated before, increasing the value of the Brownian motion parameter N_B causes

increases in the velocity and temperature profiles while its volume fraction profile decreases. This yields reductions in both the local Nusselt number and the local Sherwood number and enhancement in the local skin-friction coefficient.

The effects of the buoyancy ratio parameter N_r are illustrated in Figure 12. It is observed from this figure that an increase in the value of N_r tends to decrease all of the local skinfriction coefficient, the local Nusselt number and the local Sherwood number.

The effects of mixed convection parameter λ on skin friction coefficient, local Nusselt number and local Sherwood number are are shown in Figure 13. It is observed that as λ increases all the skin friction coefficient, local Nusselt number and local Sherwood number increases. The physical reason is that the positive buoyancy force ( λ > 0 ) implies favorable pressure gradient, and the fluid gets accelerated, which results in thinner momentum boundary layer. Consequently, the local skin friction is also increased at all times.

The effects of both the first resistance parameter γ and the second resistance parameter Δ on the local skin-friction coefficient, the local Nusselt number and the local Sherwood number are illustrated in Figure 13 and Figure 14, respectively. From these figures, it is observed that both γ and Δ have the same behavior. From these figures, we conclude that both the first and second resistances enhance all of the local rate of shear stress (the local skin-friction coefficient), the local rate of heat transfer (local Nusselt number) and the local rate of mass transfer (local Sherwood number).

The effects of Prandtl number Pr and Schmidt number Scon the behaviors of the local skin-friction coefficient, local Nusselt number and the local Sherwood number are illustrated on Figure 16 and Figure 17, respectively. It is observed from these figures that increasing the value of Pr orSc produces increases in the local rate of shear stress and decreases in both of the heat transfer rate and the mass transfer rate. This is associated with the decreases in the nanoparticle boundary layer as Sc increases as discussed earlier.

Finally, from Figures 10-17, it is observed that as the unsteadiness parameter τ increases the local rate of shear stress increases while, both of the local rate of heat transfer and the local rate of mass transfer decreases.

4. Conclusions

In the present work, the problem of unsteady mixed convection flow along a sym- metric wedge embedded in a porous medium saturated with a nanofluidis studied theoretically. The resulting system of nonlinear partial differential equations is treated using the Sparrow-Quack-Boerner local non-similarity numerical method. The obtained system is solved numerically using an efficient numerical shooting technique with a fourth-fifth-order Runge-Kutta method scheme (MATLAB pac- kage). The solutions for the flow and the heat and mass transfer characteristics are evaluated numerically for various values of the governing parameters, namely the unsteadiness parameter τ , mixed convection parameter λ , nanofluid buo- yancy ratio parameter N_r, thermophoresis parameter N_T, Brownian motion parameter N_B, first resistance parameter γ , second resistance parameter Δ , Pran- dtl number Pr and the Schmidt number Sc. The following are brief summary conclusions drawn from the analysis:

1) The thickness of the momentum boundary layer slightly decreases with an increase in the nanofluid buoyancy ratio parameterN_r. However, it increases with the increase of all other parameters.

2) The thickness of the thermal boundary layer increases with an increase in both of N_T and N_B parameters while it slightly increases with an increase inN_r parameter.

3) The nanoparticle volume fraction boundary-layer thickness decreases obviously with an increase in both N_T andN_B parameters and slightly increases with an increase in theN_r parameter.

4) The buoyancy parameter λ enhances obviously the momentum boundary- layer thickness and slightly reduce both of the thermal and the nanoparticle volume fraction boundary-layer thicknesses.

5) Both the first and the second resistance parameters γ and Δ enhance the momentum boundary-layer thickness and reduce both of the thermal and the nanoparticle volume fraction boundary-layer thicknesses.

6) The magnitude of the skin-friction coefficient f ″ ( 0 ) decreases with increasing values of the nanofluid buoyancy ratio parameter N_r and increases with the increase of all other parameters.

7) The local Nusselt number decreases with all the nanofluid parametersN_T,N_B andN_r. While, it decreases with increases in the Prandtl number Pr and the Schmidt numberSc.

8) The local Sherwood number increases as theN_T parameter increases. However, it decreases as either ofN_B,N_r,Pr or Sc increases.

9) The unsteadiness parameter τ and the resistance parameters γ and Δ enhance the local skin-friction coefficient, local Nusselt number and the local Sher- wood number.

Acknowledgements

The authors would like to thank Institute of Scientific Research and Revival of Islamic Heritage at the Umm Al-Qura University (Project ID No. 43405010) for the financial support.

Conflicts of Interest

The authors declare no conflicts of interest regarding the publication of this paper.

References

1. Smith, S.H. (1967) The Impulsive Motion of a Wedge in a Viscous Fluid. Zeitschrift für Angewandte Mathematik und Physik, 18, 508-522. https://doi.org/10.1007/BF01601721

2. Nanbu, K. (1969) Unsteady Falkner-Skan Flow. Zeitschrift für Angewandte Mathematik und Physik, 22, 1167-1172. https://doi.org/10.1007/BF01590885

3. Hall, M.G. (1969) The Boundary-Layer over an Impulsively Started Flat Plate. Proceedings of the Royal Society of London, Series A, 310, 401-414. https://doi.org/10.1098/rspa.1969.0083

4. Harris, D.S., Ingham, D.B. and Pop, I. (2002) Unsteady Heat Transfer in Impulsive Falkner-Skan Flows: Constant Wall Temperature Case. European Journal of Mechanics-B/Fluids, 21, 447-468. https://doi.org/10.1016/S0997-7546(02)01193-7

5. Keller, H.B. and Cebeci, T. (1971) Accurate Numerical Methods for Boundary-Layer Flows. Part 1. Two-Dimensional Laminar Flows. Proceedings of the 2nd International Conference on Numerical Methods in Fluid Dynamics, Berkeley, 15-19 September 1970, 92-100. https://doi.org/10.1007/3-540-05407-3_15

6. Hossain, M.A. (1992) Simultaneous Heat and Mass Transfer on Oscillatory Free Convection Boundary-Layer Flow. International Journal of Energy Research, 12, 205-216. https://doi.org/10.1002/er.4440120203

7. Hossain, M.A. (1992) Effect of Transpiration on Combined Heat and Mass Transfer in Mixed Convection along a Vertical Plate. International Journal of Energy Research, 16, 761-769. https://doi.org/10.1002/er.4440160812

8. Hossain, M.A., Banu, N. and Nakayama, A. (1994) Non-Darcy Forced Convection Boundary-Layer Flow over a Wedge Embedded in a Saturated Porous Medium. Numerical Heat Transfer, Part A: Applications, 26, 399-414.https://doi.org/10.1080/10407789408956000

9. Watkins, C.B. (1976) Unsteady Heat Transfer in Impulsive Falkner-Skan Flows. International Journal of Heat and Mass Transfer, 19, 395-403.https://doi.org/10.1016/0017-9310(76)90095-8

10. Williams, J.C. and Rhyne, T.B. (1980) Boundary-Layer Development on a Wedge Impulsively Set into Motion. SIAM Journal on Applied Mathematics, 38, 215-224.https://doi.org/10.1137/0138019

11. Ibrahim, F.S., Hassanien, I.A. and Bakr, A.A. (2004) Unsteady MHD Micropolar Fluid Flow and Heat Transfer over a Vertical Porous Plate through Porous Medium in the Presence of Thermal and Mass Diffusion with Constant Heat Source. Canadian Journal of Physics, 82, 775-790. https://doi.org/10.1139/p04-021

12. Hassanien, I.A., Ibrahim, F.S. and Omer, G.M. (2004) Unsteady Free Convection Flow in the Stagnation-Point Region of a Rotating Sphere Embedded in a Porous Medium. Mechanics and Mechanical Engineering, 7, 89-98.

13. Ibrahim, F.S. (2008) Unsteady Mixed Conviction Flow in the Stagnation Region of a Three Dimensional Body Embedded in a Porous Medium. Nonlinear Analysis: Modeling and Control, 13, 31-46. https://doi.org/10.15388/NA.2008.13.1.14587

14. Al-Harbi, S.M. and Hassanien, I.A. (2011) Unsteady Mixed Convection Flow in the Stagnation Region of a Heated Vertical Plate Embedded in a Variable Porosity Medium with Thermal Dispersion Effects. In: Dos Santos Bernardes, M.A., Ed., Developments in Heat Transfer, IntechOpen, London.

15. Al-Harbi, S.M. (2011) Self-Similar Solution of Heat and Mass Transfer of Unsteady Mixed Convection Flow on a Rotating Cone Embedded in a Porous Medium Saturated with a Rotating Fluid. Applied Mathematics, 2, 1196-1203.http://dx.doi.org/10.4236/am.2011.210166

16. Sameh, E. A, and Rashed, Z.Z., 2018, Unsteady MHD Mixed Convection Flow with Slip of a Nanofluid in the Stagnation Region of an Impulsively Rotating Sphere with Effects of Thermal Radiation and Convective Boundary Conditions. World Journal of Mechanics, 8, 137-160. https://doi.org/10.4236/wjm.2018.85011

17. Pop, I. and Ingham, D.B. (2001) Convective Heat Transfer: Mathematical and Computational Modeling of Viscous Fluids and Porous Media. Pergamon, Oxford.

18. Nield, D.A. and Bejan, A. (2005) Convection in Porous Media. Springer-Verlag, New York.

19. Choi, S.U.S. (1995) Enhancing Thermal Conductivity of Fluid with Nanoparticles, Developments and Applications of Non-Newtonian FED-Vol. 231/MD-Vol. 66, ASME, New York, 99-105.

20. Khanafer, K., Vafai, K. and Lightstone, M. (2003) Buoyancy-Driven Heat Transfer Enhancement in a Two-Dimensional Enclosure Utilizing Nanofluids. International Journal of Heat and Mass Transfer, 46, 3639-3653. https://doi.org/10.1016/S0017-9310(03)00156-X

21. Nnanna, A.G.A. (2007) Experimental Model of Temperature-Driven Nanofluid. Journal of Heat Transfer, 129, 697-704. https://doi.org/10.1115/1.2717239

22. Oztop, H.F. and Abu-Nada, E. (2008) Numerical Study of Natural Convection in Partially Heated Rectangular Enclosures Filled with Nanofluids. International Journal of Heat and Fluid Flow, 29, 1326-1336.https://doi.org/10.1016/j.ijheatfluidflow.2008.04.009

23. Aminossadati, S.M. and Ghasemi, B. (2009) Natural Convection Cooling of a Localised Heat Source at the Bottom of a Nanofluid-Filled Enclosure. European Journal of Mechanics-B/Fluids, 28, 630-640. https://doi.org/10.1016/j.euromechflu.2009.05.006

24. Yacob, N., Ishak A. and Pop I. (2011) Falkner-Skan Problem for a Static or Moving Wedge in Nanofluids. International Journal of Thermal Sciences, 50, 133-139.https://doi.org/10.1016/j.ijthermalsci.2010.10.008

25. Gorla, R.S.R., Chamkha, A.J. and Rashad, A.M. (2011) Mixed Convective Boundary Layer Flow over a Vertical Wedge Embedded in a Porous Medium Saturated with a Nanofluid: Natural Convection Dominated Regime. Nanoscale Research Letters, 6, Article No. 207. https://doi.org/10.1186/1556-276X-6-207

26. Hamad, M., Pop, I. and Ismail, A. (2011) Magnetic Field Effects on Free Convection Flow of a Nanofluid Past a Semi-Infinite Vertical Flat Plate. Nonlinear Analysis: Real World Applications, 12, 1338-1346. https://doi.org/10.1016/j.nonrwa.2010.09.014

27. Buongiorno, J. (2006) Convective Transport in Nanofluids. Journal of Heat Transfer, 128, 240-250. https://doi.org/10.1115/1.2150834

28. Hamad, M. (2011) Analytical Solution of Natural Convection Flow of a Nanofluid over a Linearly Stretching Sheet in the Presence of Magnetic Field. International Communications in Heat and Mass Transfer, 38, 487-492.https://doi.org/10.1016/j.icheatmasstransfer.2010.12.042

29. Khan, W. and Pop, I. (2010) Boundary-Layer Flow of a Nanofluid Past a Stretching Sheet. International Journal of Heat and Mass Transfer, 53, 2477-2483.https://doi.org/10.1016/j.ijheatmasstransfer.2010.01.032

30. Abu-Nada, E. and Chamkha, A. (2010) Effect of Nanofluid Variable Properties on Natural Convection in Enclosures Filled with a CuO-EG-Water Nanofluid. International Journal of Thermal Sciences, 49, 2339-2352.https://doi.org/10.1016/j.ijthermalsci.2010.07.006

31. Das, S., Choi, S., Yu, W. and Pradeep, T. (2007) Nanofluids: Science and Technology. Wiley, New Jersey. https://doi.org/10.1002/9780470180693

32. Kakaç, S. and Pramuanjaroenkij, A. (2009) Review of Convective Heat Transfer Enhancement with Nanofluids. International Journal of Heat and Mass Transfer, 52, 3187-3196. https://doi.org/10.1016/j.ijheatmasstransfer.2009.02.006

33. Muthtamilselvan, M., Kandaswamy, P. and Lee, J. (2010) Heat Transfer Enhancement of Copper-Water Nanofluids in a Lid-Driven Enclosure. Communications in Nonlinear Science and Numerical Simulation, 15, 1501-1510.https://doi.org/10.1016/j.cnsns.2009.06.015

34. Hady, F.M., Ibrahim, F.S., El-Hawary, H.M.H. and Abdelhady, A.M. (2012) Forced Convection Flow of Nanofluids Past Power Law Stretching Horizontal Plates. Applied Mathematics, 3, 121-126. http://dx.doi.org/10.4236/am.2012.32019

35. Hady, F.M., Ibrahim, F.S., El-Hawary, H.M.H. and Abdelhady, A.M. (2012) Effect of Suction/Injection on Natural Convective Boundary-Layer Flow of a Nanofluid Past a Vertical Porous Plate through a Porous Medium. Journal of Modern Methods in Numerical Mathematics, 3, 53-63.

36. Hady, F.M., Ibrahim, F.S., Abdel-Gaied, S.M. and Eid, M.R. (2012) Radiation Effect on Viscous Flow of a Nanofluid and Heat Transfer over a Nonlinearly Stretching Sheet. Nanoscale Research Letters, 7, Article No. 229.https://doi.org/10.1186/1556-276X-7-229

37. James, M., Mureithi, E.W. and Kuznetsov, D. (2015) Effects of Variable Viscosity of Nanofluid Flow over a Permeable Wedge Embedded in Saturated Porous Medium with Chemical Reaction and Thermal Radiation. International Journal of Advances in Applied Mathematics and Mechanics, 2, 101-118.

38. Uddin, M., Alim, M. and Rahman, M. (2019) MHD Effects on Mixed Convective Nanofluid Flow with Viscous Dissipation in Surrounding Porous Medium. Journal of Applied Mathematics and Physics, 7, 968-982.https://doi.org/10.4236/jamp.2019.74065

39. Chamkha, A.J., Nabwey, H.A., Abdelrahman, Z.M.A. and Rashad, A.M. (2019) Mixed Bioconvective Flow over a Wedge in Porous Media Drenched with a Nanofluid. Journal of Nanofluids, 9, 24-35. https://doi.org/10.1166/jon.2020.1728

40. Lai, F.C. (2000) Mixed Convection in Saturated Porous Media. In: Vafai, K., Ed., Handbook of Porous Media, Marcel Dekker, New York, 605-661.

41. Pereyra, V. (1978) PASVA3, an Adaptive Finite Difference Fortran Program for First Order Non-Linear Boundary Value Problems. In: Childs, B., Scott, M., Daniel, J.W., Denman, E., Nelson, P., Eds., Codes for Boundary-Value Problems in Ordinary Differential Equations, Vol. 76, Springer, Berlin, Heidelberg.https://doi.org/10.1007/3-540-09554-3_4

42. Scshadri, R., Srccshylan, N. and Nath, G. (2002) Unsteady Mixed Convection Flow in the Stagnation Region of a Heated Vertical Plate Due to Impulsive Motion. International Journal of Heat and Mass Transfer, 45, 1345-1352.https://doi.org/10.1016/S0017-9310(01)00228-9