Homotopy semi-numerical modeling of non-Newtonian nanofluid transport external to multiple geometries using a revised Buongiorno Model Homotopy Semi-Numerical Modeling of Non-Newtonian Nanoﬂuid Transport External to Multiple Geometries Using a Revised Buongiorno Model

Title Homotopy semi-numerical modeling of non-Newtonian nanofluid transport external to multiple geometries using a revised Buongiorno Model Abstract: A semi-analytical solution for the convection of a power-law nanoﬂuid external to three di ﬀ erent geometries (i.e., cone, wedge and plate), subject to convective boundary condition is presented. A revised Buongiorno model is employed for the nanoﬂuid transport over the various geometries with variable wall temperature and nanoparticle concentration conditions (non-isothermal and non-iso-solutal). Wall transpiration is included. The dimensional governing equations comprising the conservation of mass, momentum, energy and nanoparticle volume fraction are transformed to dimensionless form using appropriate transformations. The transformed equations are solved using a robust semi-analytical power series method known as the Homotopy analysis method (HAM). The convergence and validation of the series solutions is considered in detail. The variation of order of the approximation and computational time with respect to residual errors for temperature for the di ﬀ erent geometries is also elaborated. The inﬂuence of thermophysical parameters such as wall temperature parameter, wall concentration parameter for nanoﬂuid, Biot number, thermophoresis parameter, Brownian motion parameter and suction / blowing parameter on the velocity, temperature and nanoparticle volume fraction is visualized graphically and tabulated. The impact of these parameters on the engineering design functions, e.g., coe ﬃ cient of skin fraction factor, Nusselt number and Sherwood number is also shown in tabular form. The outcomes are compared with the existing results from the literature to validate the study. It is found that thermal and solute Grashof numbers both signiﬁcantly enhance the ﬂow velocity whereas they suppress the temperature and nanoparticle volume fraction for the three di ﬀ erent conﬁgurations, i.e., cone, wedge and plate. Furthermore, the thermal and concentration boundary layers are more dramatically modiﬁed for the wedge case, as compared to the plate and cone. This study has substantial applications in polymer engineering coating processes, ﬁber technology and nanoscale materials processing systems. to previous published works. A detailed parametric study is conducted of the inﬂuence of geometry and key nanoscale parameters on heat, mass and momentum characteristics with deep physical elaboration. The present investigation furnishes important benchmarks for future simulations of relevance to nanoscale polymeric coating analysis.


Introduction
The Convective flow over different geometrical configurations (e.g., wedge, cone, plate) arises in a diverse spectrum of engineering systems including polymer engineering, fiber technology, geothermal transport in the vicinity of rock formations, thermal-aerodynamics of high-speed flight, nuclear cooling systems, spray deposition techniques, coating and surface treatment, etc. Numerous studies have been communicated in this regard using both analytical and computational methods for multiple geometries and for both Newtonian and non-Newtonian fluids. Roy [1] studied the natural convection from vertical cone at high Prandtl number. Vajravelu and Nayfeh [2] inspected the effect of heat generation or absorption on the flow and heat transfer along the surface of a cone under magnetohydrodynamics (MHD) effects. Natural convection from a wedge and cone subject to mixed thermal boundary conditions has been addressed by Ramanaiah and Malaryizhi [3]. Watanabe et al. [4] have analyzed the response of mixed convection boundary layer flow over a wedge to the effect of surface transpiration, i.e., suction/blowing. Fang and Lee [5] considered the flow of rarefied gas free stream over a moving flat plate. Xu and Liao [6] discussed the flow and heat transfer of an Ostwald-De Waele power-law fluid over a stretching flat plate. Natural convection heat and mass transfer from a vertical cone in a porous medium with variable wall concentration was analyzed by Chen [7]. Gorla et al. [8] described the influence of mixed convection on the flow of nanofluid over a vertical porous wedge. Rahman et al. [9] studied flow of nanofluid from a wedge with convective boundary condition under the influence of magnetic field and in presence of heat generation/absorption. Chamkha and Rashad [10] studied the influence of chemical reaction and cross-diffusion effects on the MHD mixed convection flow from rotating vertical cone. A detailed study of heat transfer of third grade viscoelastic fluid flow over non-isothermal wedge using homotopy analysis method was presented by Rashidi et al. [11]. Recently, the flow of non-Newtonian fluids (with different constitutive models) over different geometries has attracted the focus of several researchers owing to its relevance to, for example, thermal polymer coating fabrication processes. Vasu and Kumar [12] examined the transient behaviour of laminar natural convective flow of a nanofluid over both a vertical cone and plate. Vasu and Gorla [13] have examined the mixed convection flow of Al 2 O 3 -water based nanofluid. Vasu et al. [14] used the finite difference method to simulate the influence of thermophoresis and heat sink/source on double-diffusive convection of a short memory viscoelastic fluid with thermal radiative flux effects. Very recently, Vasu et al. [15] reported on numerical solutions for transient mixed convection flow of a nanofluid in the forward stagnation region of a spinning sphere under the nonlinear Boussinesq approximation. Ray et al. [16] studied electrically-conducting viscoplastic nanofluid bioconvection thin film transport phenomena from a time-dependent extending sheet. Sreenivasulu et al. [17] used Lie algebra and computational solvers to investigate the radiative heat transfer and slip effects on oblique hydromagnetic flow of a tangent hyperbolic (shear-thinning) fluid containing carbon nanotubes.
It is well-known that a wide variety of fluids encountered in chemical and process mechanical engineering applications do not satisfy Newton's law of viscosity. These liquids are generally termed non-Newtonian or rheological fluids [18]. These include polymers, gels, colloids, lubricants and nano-coatings. Owing to the diverse characteristics of such liquids at different shear rates, a single relationship between shear stress and shear rate cannot be used to categorize all rheological behaviour. The classical Navier-Stokes equation is inadequate for simulating many of the characteristics of non-Newtonian liquids including normal stress differences, spurt, thixotropy, memory, relaxation, retardation, yield, etc. Therefore, various models have been developed to more precisely predict the behaviour of non-Newtonian fluid flow [19,20]. An elegant but simple constitutive equation for modeling viscosity variation in non-Newtonian fluids is the Ostwald-de-Waale model [21] which has considerable practical significance and in fact applies to a broad range of liquids in medicine, materials processing, petro-chemical engineering, etc. Acrivos et al. [22] and Schowalter [23] were among the earliest researchers to construct similarity solutions for boundary layer flows of power-law fluid. Many investigations in the past few decades have subsequently extended these earlier studies to more complex scenarios. Gorla and Kumari [24] obtained non-similar solutions for mixed convection in power-law fluid along a vertical plate embedded in an isotropic Darcian porous medium. Ikbal et al. [25] deployed the power-law model to compute the unsteady blood flow through stenosed artery in the presence of magnetic field. Mahapatra et al. [26] explored the heat transfer in power-law fluid magnetohydrodynamic boundary layer flow with lateral mass flux effects. Qi et al. [27] used the Lattice Boltzmann method to compute wake effects on particle and power-law fluid flow interaction. Vasu [28] studied viscous heating and Ohmic heating effects on natural magneto-convection flow of power-law nanofluid film along an inclined surface. Bég et al. [29] computed the Von Karman swirling flows of power-law liquids in isotropic, sparsely-packed non-linear porous media.
Improving the thermal conductivity of conventional fluids such as water, refrigerant, ethyl glycol and engine oil, has emerged as a key area of intensive activity with the objective of enhancing the heat transfer characteristics (thermal conductivity, heat transfer coefficient, etc.). Nanofluids [30][31][32] are nanoscale-engineered convectional fluids comprising metallic or non-metallic nano-particles which have significantly improved thermal performance in numerous systems. These include but are not restricted to cooling of electrical appliances, heat-exchanging devices, nuclear reactor technology, solar water heating, biomedical science, diesel generators and aerospace/mechanical coatings. The mathematical models which have been developed to describe the transport phenomena in nanofluids are generally either homogeneous [33][34][35] or non-homogeneous models [36][37][38]. In the homogeneous (single component) model, thermophysical properties of the base fluid are modified with the correlations of nanoparticles. The non-homogeneous (two-component) model is more elaborate and was introduced by Buongiorno [31]. It features multiple slip mechanisms including Inertia, Gravity, Magnus effect, Diffusiophoresis, Thermophoresis, Brownian motion and Fluid drainage. However, the Buongiorno model (which focuses on Brownian dynamics and thermophoretic body forces as the principal nanoscale mechanisms) has been modified to consider the influence of the distribution of nanoparticle volume fraction at a boundary. This "revised model" has been extensively implemented in computational nanoscale studies in recent years. Kuznetsov and Nield [39] applied the revised model in natural convection flow of nanofluid from a vertical plate. Malvandi et al. [40] investigated the effect of Brownian motion and thermophoresis on hydromagnetic flow inside circular micro-channel filled with alumina/water nanofluid. Moshizi et al. [41] have conducted a theoretical study of alumina-water nanofluid flow inside a concentric pipe under the influence of heat absorption/generation. Kameswaran et al. [42] have considered the modified Buongiorno model to study the mixed convection from a wavy surface embedded in porous medium. Bég et al. [43] investigated the hydrodynamic, thermal and mass slip effects in transient asymmetric bio-convective nanofluid flow in a porous micro-channel with deformable walls using MAPLE numerical quadrature.
To the authors' knowledge, the vast majority of nanofluid convection studies in external boundary layer flows have been restricted to a single geometry. Additionally, there is no attempt to simultaneously consider the transport characteristics of non-Newtonian fluids from multiple geometries with a convective boundary condition and revised Buongiorno nanofluid model. This is the main focus of the present article. In this regard, the wedge, plate and cone geometrical configurations are selected which feature extensively in polymer processing systems. A detailed mathematical formulation of the flow of nanofluid over the cone, wedge and plate by including convective boundary condition and modified Buongiorno nanofluid model is presented. The normalized nonlinear boundary value problems are solved using the homotopy analysis method (HAM). Aspects of the convergence of the homotopy series solution are elaborated. Extensive validation of the present study is conducted with regard to previous published works. A detailed parametric study is conducted of the influence of geometry and key nanoscale parameters on heat, mass and momentum characteristics with deep physical elaboration. The present investigation furnishes important benchmarks for future simulations of relevance to nanoscale polymeric coating analysis.

Mathematical Modeling
Consider the steady, incompressible and laminar flow of a non-Newtonian power-law nanofluid over three different geometrical surfaces, i.e., vertical cone, wedge and plate. Convective conditions are imposed at the boundary surface of the geometries. The effects of Brownian movement and thermophoresis are incorporated. The surface of the geometries is assumed to be stretched with velocity u W = xν l 2 and the surfaces are assumed to be transpiration (suction/injection) velocity v W . The temperature and concentration near the surface are prescribed as T w = T ∞ + ax r 1 and C w = C ∞ + bx r 2 , where r 1 and r 2 are the wall temperature and nanoparticle concentration parameters. r 1 = 0 implies constant temperature at the wall (iso-thermal) and r 2 = 0 corresponds to constant wall nanoparticle concentration (iso-solutal). The suspension of nanoparticles is dilute and stable with no chemical reaction and negligible radiative heat transfer. The physical models for the three geometries are shown in Figure 1. T ∞ and C ∞ are ambient temperature and nanoparticle concentration respectively. The two dimensional co-ordinate axes are orientated such that the x-axis coincides with the wall (boundary) of each geometry which is prescribed convective condition. The y-axis is orthogonal to the wall as shown in Figure 1. ϕ is the semi-vertex cone angle, r is the radius of the cone and Ω is the vertex angle of the wedge. Dilute nanofluids are considered with spherical, homogenously distributed nano-particles.
Dilute nanofluids As per the quoted restraints, the general conservation equations which govern the flow over all three geometries can be expressed as: The prescribed boundary conditions at the geometry wall and in the free stream: where u and v are velocity component along x and y directions. u w = xν l 2 . ρ f and ν = µ ρ represent nanofluid density and kinematic viscosity respectively (µ is dynamic viscosity and ρ is fluid density), k is nanofluid thermal conductivity, T is temperature, n is power index, D B is Brownian diffusion coefficient, D T is thermophoretic diffusion coefficient, T ∞ is ambient temperature, C ∞ is ambient concentration and τ is nanoparticle to base fluid heat capacity ratio. Apparent viscosity for power-law Introducing the following transformations Equations (2)-(4) thereby reduce to the following non-dimensional form (Equation (1) is automatically satisfied): 1 Pr 1 Sc The wall and free stream boundary conditions (5) reduced to: Parameters appear in Equations (7)-(9): Pr = α µC p is Prandtl number, Nb = denotes Brownian motion parameters and thermophoresis parameter, represent thermal Grashof numbers and solutal (nanoparticle species buoyancy) Grashof numbers respectively. Several gradient functions are of interest in materials processing (coating) engineering design. These are the friction factor C f (dimensionless surface shear stress function), local Nusselt number Nu (dimensionless wall heat transfer gradient) and local Sherwood number Sh (dimensionless nanoparticle wall mass transfer gradient) and these take the following definitions, respectively:

Solution Using Homotopy Analysis Method
Homotopy Analysis Method (HAM) has been used to solve the model defined by Equations (7)-(9) subject to the boundary conditions (10). HAM is a semi-analytical method developed by Liao [44] and has been implemented to solve many fluid dynamics, applied mathematics and multi-physical mechanics problems. This method does not rely upon any physical small/large parameter and hence provide convenient way to guarantee the series solution by using a special parameter known as convergence control parameter. Hashim et al. [45] obtained series solutions to systems of first and second order partial differential equations (PDEs) by applying HAM. Recently, Sheikholeslami et al. [46] used HAM to analyze the melting effect in nanofluid flow. Mabood et al. [47] implemented HAM to evaluate the MHD flow over exponentially stretching sheet with radiation effect. Sarvanthi and Gorla [48] scrutinized the influence of heat sink or source on chemically-reacting Maxwell elastic-viscous nanofluid over an exponential stretching surface with convective boundary conditions. Hayat et al. [49] examined the non-Fourier thermal convection of carbon nanotube fluids over a curved stretching surface. Other applications of HAM in biological fluid dynamics and nanoscale transport include [50][51][52]. To extract the solutions of Equations (7)-(9) subject to the boundary conditions (10) using HAM, we have considered initial guesses f 0 , θ 0 and φ 0 of f , θ and φ in the following form: The linear operators are selected as follows: with the following properties: are the arbitrary constants. If p ∈ [0, 1] is the embedding parameter, h 1 , h 2 and h 3 are convergence control parameters, then we construct the following zeroth-order deformation equations as: The appropriate forms of the boundary conditions are now: The nonlinear operators N 1 , N 2 and N 3 are defined based on Equations (7)-(9) as below: In Equation (17a), ε = n − 1. When p = 0 and p = 1, we obtain: The mth-order deformation equation can be obtained by differentiating the zeroth order deformation Equation 15a,d m times and dividing by m! This leads to: The associated mth order deformation boundary conditions are: Here: Therefore, as p increases from 0 to 1 then f (η; p), θ(η; p) and φ(η; p) vary from initial approximations to the exact solutions of the original nonlinear differential equations.

Convergence of Homotopy Series Solution
The homotopy solutions (21a)-(21c) for Equations (7)-(9) subject to the boundary conditions (10) is obtained using initial guesses f 0 , θ 0 and φ 0 given in Equation (12a-12c) and linear operators L 1 ( f ), L 2 (θ) and L 3 (φ) given in Equation (13a-13c) with the suitable values for the non-zero control parameters h 1 , h 2 and h 3 that have been obtained by plotting the h-curves shown in Figure 2. From Figure 2, it is seen that the valid regions of h 1 is (−0.1, −0.36) and of h 2 and h 3 is about (−0.1, −0.55). Here, h 1 = h 2 = h 3 = −0.35 is considered for the present study. Computations are performed using Mathematica 9 software with the thermophysical parameters: n = 3, Gr = Gc = 0.5,  The convergence of homotopy series solution is given Table 1 and concludes that 16th order of approximation is appropriate for computations since negligible variations are observed for higher orders of approximation. 16th order of approximation is taken for the computation throughout the study. The present study is compared with the existing results of Hassanien et al. [53] for the validation of the HAM computations as shown in Table 2. Computational time and error analysis with respect to different residual errors is explained in Table 3. Here, the obtained data elaborates that the computational time depends on the complexity of the geometries.
The convergence of homotopy series solution is given Table 1 and concludes that 16th order of approximation is appropriate for computations since negligible variations are observed for higher orders of approximation. 16th order of approximation is taken for the computation throughout the study. The present study is compared with the existing results of Hassanien et al. [53] for the validation of the HAM computations as shown in Table 2. Computational time and error analysis with respect to different residual errors is explained in Table 3. Here, the obtained data elaborates that the computational time depends on the complexity of the geometries. It is clear from Table 3 that computations require greater compilation times for the cone compared with the wedge and plate to achieve the same accuracy of order 10 −7 . With the same set of parameters, the order of approximation as well as the computational time for the plate is generally less. The computations have been performed in system processor: Intel(R) Core (TM) i5-5200U CPU@ 2.20 GHz and system type: 64-bit MS Windows 10 operating system. The symbolic software Mathematica has been used to compute the results.

Results and Discussion
The set of transformed coupled nonlinear ordinary differential Equations (7)-(9) subject to boundary condition (10) are solved very efficiently using HAM. To investigate the evolution in velocity, temperature and volume fraction profiles, we set values of the parameters, i.e., Prandtl number, Schmidt number, thermophoresis and Brownian motion parameter and suction/injection parameteras n = 3, Gr = Gc = 0.5, Pr = 2, S = 0.1, The main focus of this article is to analyze the effect of the revised Buongiorno model on the convective flow of non-Newtonian power-law fluid over different geometries. With the aid of tables, the skin factor, Nusselt number and Sherwood number are properly explained for pertinent parameters. The half-cone angle ϕ is taken 45 • . Figures 3-5 show the different response for Newtonian and non-Newtonian nanofluids via the velocity, temperature and nanoparticle concentration profiles for the cone, wedge and plate. It is clearly seen from Figure 3 that the velocity of both Newtonian and non-Newtonian fluid over the plate exceeds that computed for the wedge and plate. Additionally, the power-law nanofluid suppresses the velocity as compared to Newtonian nanofluid for all geometries, i.e., it thickens the momentum boundary layer for all three geometries. The temperature and nanoparticle profiles for Newtonian nanofluid shown in Figures 4 and 5 are higher than for non-Newtonian fluid near to the surface of the geometries. The contrary behaviour is observed away from the surface of the plate and wedge. Since, the velocity for power-law flow reduces because of higher viscosity when compared with Newtonian fluid flow. Due to this reason, there is clogging of nanoparticle near the surface of the geometries which results in increasing the nanoparticle volume fraction.
spherical and homogenously distributed in the base fluid). With greater Nb values smaller nanoparticles are present and this intensifies the thermal conduction heat transfer from the particles to the surrounding fluid ( Figure 31). However, stronger Brownian motion inhibits the diffusion of nanoparticles since smaller nanoparticles are less successful in migrating through the base fluid and are more susceptible to ballistic collisions, as noted by Shukla and Dhir [54]. Physically excessive concentrations of nanoparticles (higher volume fractions) are counter-productive in nano-coating design and intermediate sized nanoparticles have been shown to disperse more homogenously. Thermal and nanoparticle concentration boundary layer thicknesses for the wedge exceed those achieved for the cone but are less than those for the plate.        6-8 elucidate the effect of thermal Grashof number Gr on the velocity, temperature and nano-particles volume fraction (concentration) for the cone, wedge and plate. A rise in value of Grashof number enhances the velocity whereas it depletes the temperature and nanoparticle volume fraction. Therefore, hydrodynamic boundary layer thickness is reduced whereas thermal and nanoparticle concentration boundary layer thicknesses are elevated with thermal buoyancy. Increasing values of Grashof number correspond to a decrease in viscosity of the nanofluids which assist in momentum diffusion and flow acceleration. This boosts the velocity magnitudes for all geometries. It is also evident that the temperature and nanoparticle volume fraction are substantially greater for the wedge geometry as compared to the plate and cone geometries.                     Figures 15-17 visualize the effect of wall temperature (isothermal and non-isothermal effect) on velocity, temperature and nanoparticle volume fraction. Due to increase in wall temperature, the boundary layer behaviour is influenced strongly. Velocity is lower and momentum boundary layer thickness is greater for the wedge when compared with the cone whereas momentum boundary layer thickness is lower than for the plate. Velocity, temperature and nanoparticle volume fraction profiles show higher values for constant wall temperature (isothermal case) compared with the non-isothermal case.          Figure 22 indicates that the heat transfer behaviour is highly dependent on the Prandtl number, a most important parameter in thermal convection. The temperature of the fluid decreases monotonically with the increasing Pr, i.e., thermal boundary layer thickness is suppressed. Clearly different responses are computed for the different geometries indicating that external geometry is a key player in controlling thermo-fluid characteristics in nanofluid boundary layer flows. Figure 20. Effect of wall concentration parameter on nanoparticle volume fraction.           Figure 27 demonstrates that the momentum boundary layer thickness lessens with increasing value of Schmidt number since the flow is strongly accelerated. Figure 27 shows that an increase in Schmidt number results in a decrease in nanoparticle concentration values. This reduction in nano-particles concentration is due to the change in Brownian diffusion coefficient. Enlarging Schmidt number corresponds to weaker Brownian diffusion coefficient and inhibited diffusion of nanoparticles through the boundary layer regime. This also results in depletion in nanoparticle concentration boundary layer thickness.       Figures 28 and 29 depicts the effect of thermophoresis parameter on temperature and nanoparticle volume fraction respectively for different geometries. Higher thermophoresis implies greater migration of hot nanoparticles in the direction of a decreasing temperature gradient which encourages nanoparticle diffusion in the boundary layer. Thermophoretic forces exerted on the nanoparticles are in the opposite direction to the actual temperature gradient. This effectively results in a thickening in the nanoparticle concentration boundary layer thickness for all the geometries studied. Due to affect of convective boundary condition and modified Buongiorno model, the behaviour of temperature and nanoparticle volume fraction is different by varying Nt. Temperature field is a decreasing function of thermophoresis parameter Nt whereas nanoparticle volume fraction is an increasing function of Nt (see Figure 29). Thermal boundary layer thickness is therefore decreased with stronger thermophoresis. Temperature and nanoparticle volume fraction (concentration) profiles for the wedge exceed those for the plate and cone geometries.      Figure 31). However, stronger Brownian motion inhibits the diffusion of nanoparticles since smaller nanoparticles are less successful in migrating through the base fluid and are more susceptible to ballistic collisions, as noted by Shukla and Dhir [54]. Physically excessive concentrations of nanoparticles (higher volume fractions) are counter-productive in nano-coating design and intermediate sized nanoparticles have been shown to disperse more homogenously. Thermal and nanoparticle concentration boundary layer thicknesses for the wedge exceed those achieved for the cone but are less than those for the plate.   Tables 4-6 displays the effect of parameters on skin factor, Nusselt number and Sherwood number for cone, wedge and plate respectively. In all the geometries, each afore-mentioned parameter exhibited the same behaviour. It can be observed from the tables that a rise in Brownian motion, Schmidt number and wall concentration parameter, reduces the Nusselt number while the other parameters enhance Nusselt number (heat transfer rate). Thermal and solutal Grashof numbers enhance the heat and mass transfer rates. Increasing wall temperature 1 r and concentration parameter 2 r reduce the skin friction coefficient. However, friction factor and reduced mass transfer rate defined by Equation (11) are increased whereas heat transfer rate is reduced with an increase in Schmidt number.  Tables 4-6 displays the effect of parameters on skin factor, Nusselt number and Sherwood number for cone, wedge and plate respectively. In all the geometries, each afore-mentioned parameter exhibited the same behaviour. It can be observed from the tables that a rise in Brownian motion, Schmidt number and wall concentration parameter, reduces the Nusselt number while the other parameters enhance Nusselt number (heat transfer rate). Thermal and solutal Grashof numbers enhance the heat and mass transfer rates. Increasing wall temperature r 1 and concentration parameter r 2 reduce the skin friction coefficient. However, friction factor and reduced mass transfer rate defined by Equation (11) are increased whereas heat transfer rate is reduced with an increase in Schmidt number.

Conclusions
The present study investigates the external incompressible, steady-state boundary layer flow, heat and mass transfer behaviour in non-Newtonian power-law dilute nanofluid flow from three different over geometries (cone, wedge and plate) in the presence of variable wall thermal and solutal characteristics and wall transpiration. The effect of convective boundary conditions is considered and the revised Buongiorno nanofluid model is also deployed with focus on thermophoretic body forces and Brownian motion. The conservation equations are transformed with appropriate similarity variables and the resulting non-dimensional, nonlinear boundary value problems are solved using robust homotopy analysis method. The major findings of the study can be summarized as follows: • Momentum boundary layer thickness is more significantly modified for the plate compared to the cone and wedge whereas thermal and concentration boundary layer is more significantly altered for the wedge geometry.

•
Due to the boundary condition defined by modified Buongiorno model, the mass transfer rate (Sherwood number) decreases with increase in thermophoresis parameter and increases with increase in Brownian motion parameter for all the geometries.

•
Brownian motion exerts a stronger influence on mass transfer rates (Sherwood numbers) for all the geometries when compared to heat transfer rates (Nusselt numbers).

•
Increasing thermal Grashof number and solutal (nanoparticle) Grashof number reduce temperature and enhance wall heat transfer rates (Nusselt numbers).

•
Increasing thermal Grashof number and solutal (nanoparticle) Grashof number both elevate the skin friction factor for all geometries considered (cone, wedge and plate).

•
The Nusselt number for the cone is in excess of that for either a plate or wedge.

•
The convective boundary condition parameter, i.e., Biot number, controls the thermal and concentration boundary layer significantly and improves the heat transfer rates (Nusselt numbers) and in particular achieves high magnitudes for the wedge (Falkner-Skan case).

•
Increasing non-isothermal behaviour (rising value of wall temperature parameter) and non-iso-solutal effect (greater wall concentration parameter) magnifies the heat transfer and mass transfer rates (i.e., Nusselt and Sherwood numbers) for all geometries. • Thermal and solute Grashof number enhances the momentum boundary layer and suppresses the thermal and concentration boundary layer for all geometries. Both thermal and species buoyancy force therefore increase the rate of heat and nanoparticle mass transfer to the wall.

•
The dominant effect of increasing Schmidt number is to lessen the heat transfer rate and enhance skin friction and nanoparticle mass transfer rate for all geometries. • Increasing Schmidt number, Brownian motion and non-iso-solutal wall parameter deplete the Nusselt numbers, i.e., reduce heat transfer rates at the walls of all the geometries studied.
The present investigation has revealed some interesting insights into nanofluid transport characteristics in boundary layer flows from multiple geometries. However, it has been restricted to Fourier thermal conduction models. Future studies will consider thermal relaxation effects which may be simulated with non-Fourier models [55,56] and will be communicated imminently.