Numerical Solutions for Laminar Boundary Layer Nanofluid Flow along with a Moving Cylinder with Heat Generation, Thermal Radiation, and Slip Parameter

Agbaje, Titilayo Morenike; Makanda, Gilbert

doi:https://doi.org/10.1155/2021/8288534

Abstract and Applied Analysis

On this page

Abstract Introduction Results and Discussion Conclusion Data Availability Conflicts of Interest Acknowledgments References Copyright Related Articles

Research Article | Open Access

Volume 2021 | Article ID 8288534 | https://doi.org/10.1155/2021/8288534

Numerical Solutions for Laminar Boundary Layer Nanofluid Flow along with a Moving Cylinder with Heat Generation, Thermal Radiation, and Slip Parameter

Titilayo Morenike Agbaje¹and Gilbert Makanda¹

Academic Editor: Fasma Diele

Received15 Apr 2021

Accepted05 Nov 2021

Published01 Dec 2021

Abstract

The investigation of the numerical solution of the laminar boundary layer flow along with a moving cylinder with heat generation, thermal radiation, and surface slip effect is carried out. The fluid mathematical model developed from the Navier-Stokes equations resulted in a system of partial differential equations which were then solved by the multidomain bivariate spectral quasilinearization method (MD-BSQLM). The results show that increasing the velocity slip factor results in an enhanced increase in velocity and temperature profiles. Increasing the heat generation parameter increases temperature profiles; increasing the radiation parameter and the Eckert numbers both increase the temperature profiles. The concentration profiles decrease with increasing radial coordinate. Increasing the Brownian motion and the thermophoresis parameter both destabilizes the concentration profiles. Increasing the Schmidt number reduces temperature profiles. The effect of increasing selected parameters: the velocity slip, Brownian motion, and the radiation parameter on all residual errors show that these errors do not deteriorate. This shows that the MD-BSQLM is very accurate and robust. The method was compared with similar results in the literature and was found to be in excellent agreement.

1. Introduction

The boundary layer flow on heat and mass transfer has been overmoving, and stretching surfaces have been studied and remained an active area in the past decade. This is because it has numerous applications in areas such as hot rolling, processes of polymer extrusions, wire drawing, extrusions in aerodynamic plastic sheets, process of condensation in metallic plates during cooling, and many other applications. According to Poply et al. [1], the study of flows in cylinders is considered two-dimensional when the cylinder radius is much larger than the boundary layer thickness. For lean and thin cylinders, the two dimensions may be of the same order; in this case, the flow is referred to as axisymmetric rather than two-dimensional. These dimensions affect velocity, temperature, and concentration profiles which in turn affect the skin friction coefficient.

In light of the importance of laminar boundary layer fluid flow, many researchers have considered several flow geometries with different boundary conditions, using many different techniques to solve similar models. These include the research carried out by Shateyi and Marewo [2] who used the successive linearization method (SRM) to solve a problem on laminar boundary layer flow and heat transfer in nonlinear differential equations; they also considered stretching cylinder, porous media, and thermal conductivity. In their investigation, they observed that the curvature significantly affects temperature and velocity fields. Also, both the skin friction coefficient and local Nusselt number increase as the curvature increases. Rangi and Naseem [3] used the Keller-box technique to solve the equations describing boundary layer flow of heat transfer with nonconstant thermal conductivity along with a stretching cylinder. Their results also showed that the cylinder curvature affects temperature, velocity, and skin friction fields. In the results, thermal conductivity and curvature aid heat transfer and reduction in fluid viscosity aids the convectional heat transfer rate. Numerical solutions concerning boundary layer flow, heat transfer along a stretching cylinder in porous media were obtained by Mukhopadhyay [4] using the shooting method. The results of this study found that increasing permeability parameter results in a decrease in velocity profiles. The skin friction and rate of heat transfer are much less on the flat plate than on the cylinder surface. Some other laminar boundary layer fluid models by different researchers including Elbashbeshy et al. [5], Lin and Shih [6], and Ali and Alabdulkarem [7] among others investigated the Casson fluid flow on a stretching surface based on an exponential model; they also applied theoretical analysis using lie groups in MHD fluid flow. In these works, a fourth-order Runge-Kutta method was used.

The study of the fluid flow and the physical properties of nanofluids has been widely studied in the past few years due to the wide applications of these fluids in cancer therapy, fuel cells, electronics, and pharmaceutical processes, just to mention a few. The concept of nanofluids was introduced by Choi [8], where the introduction for the proposal of nanoparticles suspended in base fluids such as ethylene glycol, oil, and water. A nonhomogeneous equation with two components for nanofluids was introduced by Buongiorno [9], where the seven slip mechanisms were proposed. The results explained effects in thermophoresis Brownian motion in nanofluids. Alamri et al. [10] used the homotopy analysis method to investigate fluid flow in nanofluids in the presence of second-order slip and Stefan blowing effects in a tube (Poiseuille). In this investigation, the results show that the retarding effects of Stefan blowing is observed for temperature and velocity profiles, the opposite effect was noticed for the case of particle concentrations. A more enhanced response in the field of velocity was observed in the slip of the second order than in the slip of the first order. Nadeem et al. [11] used the homotopy analysis method in the solutions for nonlinear differential equations in which the oblique stagnation point flow was considered for a Casson nanofluid flow with stretching surface and heat transfer. The Runge-Kutta method of the fourth order together with the shooting technique was used by Khan et al. [12] to find a numerical solution of the Maxwell nanofluid stagnation point fluid flow over a stretching sheet with slip conditions and chemical reaction effects. Their results showed that the skin friction coefficient is inversely proportional to slip parameter but an opposite is noticed in the case of fluid relaxation parameter.

Dhanai et al. [13] investigated MHD boundary layer fluid flow on multiple solutions of and heat transfer power-law nanofluids, viscous dissipation, and permeable shrinking/stretching effects in which the shooting method was used. In this study, viscous dissipation was significant and the Brownian motion was neglected in the case of heat transfer. The results of the analysis of the mass and heat transfer in nanofluid flow in Casson fluid flow between parallel plates with Hall current effects were obtained by Shah et al. [14] using an optimal and numerical method. The results showed that Hall currents decrease conductivity, which then increase the velocity and temperature of the fluid. Nadeem and Khan [15] obtained dual solutions of inclined stagnation point nanofluid flow with MHD over an oscillatory shrinking/stretching sheet using the fourth-order Runge-Kutta method. The results obtained indicated that dual solutions occur in both cases that is the shrinking and stretching cases. Furthermore, the obtained lower solution branch shows the same behaviour for the coefficient of skin friction in the shrinking case. On the contrary, the solution in the upper branch has perfect behaviour in both the shrinking and stretching cases. More recent studies of nanofluid flow include among others Kamal et al. [16], Besthapu et al. [17], Sadiq et al. [18], Nayak et al. [19], Yousif et al. [20], Salawu and Ogunseye [21], Kumar and Srinivas [22], Sreedevi and Reddy [23, 24], Gireesha et al. [25, 26], Archana et al. [27], and Awais et al. [28].

To the best of the authors’ knowledge, the study of laminar boundary layer nanofluid flow along a fixed or moving cylinder with heat generation, radiation parameter, and slip parameter has not been widely studied. This work is aimed at studying the effects of thermal radiation, heat generation, slip effects, and other important parameters discussed in Section 4. The model of Lin and Shih [6] is modified into a laminar boundary layer nanofluid model. The fluid flow includes the effects of thermophoretic forces, thermal radiation, heat generation, magnetic field, and Brownian motion. We will focus more on the nanoparticle slip boundary condition, thermophoretic force, and Brownian motion.

The solution of the partial differential equations is obtained by the multidomain bivariate spectral quasilinearization method (MD-BSQLM). This method uses more than one domain in its collocation process; it uses the quasilinearization method based on the Newton–Raphson technique previously used by Bellman and Kalaba [29]. The obtained system of equations is used together into multiple subintervals using the spectral method based on the Chebyshev and Lagrange interpolation polynomials. These are used on the linearized nonlinear equations of PDEs independently in both directions of time and space. This approach of the MD-BSQLM gives very accurate solutions. The method yields more accurate solutions when compared to other methods such as bivariate spectral quasilinearization [30], bivariate spectral homotopy analysis method [31], and bivariate spectral relaxation method (BSRM) [32], among others.

This paper mostly focuses on the robustness and accuracy of the MD-BSQLM and in finding solutions of problems of high complexity. The robustness and accuracy of the numerical method are obtained by analysis and computations of solutions of the momentum, heat, and mass transfer systems of equations. Obtained solutions are compared to those found in the literature, and an acceptable agreement is observed. The effect of different physical parameters on the temperature, velocity, and concentration fields is discussed in tabular and graphical forms in Section 4. The results show the different behaviours that occur when these parameters are changed. The paper will be arranged as follows: the mathematical formulation is in Section 2, the method of solution is in Section 3, and the results will be discussed in Section 4 and conclusions in Section 5.

2. Mathematical Formulation of the Problem

An axisymmetric, steady boundary layer flow of an incompressible viscous fluid along with a cylinder in motion as shown in Figure 1 is considered. The cylinder is moving with a constant in a . The , , , and are considered. The model will consider the motion of the cylinder in both directions. The mathematical model is presented as subject to the boundary conditions

To solve this problem, we introduce the following nondimensional groups and variables: where is the stream function which is defined as which satisfies the continuity equation (1). Substituting equation (9) in equations (2)–(5) gives subject to the boundary conditions: where prime defines derivatives with respect to , magnetic parameter is , Prandtl number is , thermal radiation parameter is , Eckert number is , , nondimensional heat generation parameter is , Brownian motion parameter is , thermophoresis parameter is , Schmidt number is , and nondimensional velocity slip parameter is .

3. Multidomain Bivariate Spectral Quasilinearization Method

In this section, we describe how to apply the multidomain bivariate spectral quasilinearization method (MD-BSQLM) for solving the system of coupled PDEs that are highly coupled as shown in equations (13)–(15). The solution is decomposed into a large interval with small subdivisions. The solution is evaluated at each time step at the end of each subinterval. The multidomain approach is applied in the direction. The system of equations is first linearized using the quasilinearization (QLM) as indicated in [29, 33]. The QLM approach makes use of the Taylor series approximation to linearize the system of equations. In this method, we assume that the value of the function on the current and previous iterations is small. Applying the QLM on (13)–(15) gives the following: where

Now, let , where and the domain is decomposed into nonoverlapping intervals as

The PDEs are solved separately at each of the subintervals . Solutions of the th interval are obtained by using solutions from the th interval. This process is implemented as follows: the boundary conditions are imposed

An appropriate initial guess is imposed to start the multidomain iteration. This gives the solution in the first interval and satisfies the boundary condition (16). The initial conditions at the next interval are prescribed by the conditions that ensure continuity as follows:

The domains in the and directions have transformed the domain for computation, at each subinterval using the linear transformation

where is a number that approximates infinity in . The Chebyshev-Gauss-Lobatto nodes defined in [34, 35] are used for collocation as follows:

where and are all the points in - and -directions used for collocation, respectively.

If the solutions , and are approximated at each subinterval by using the bivariate Lagrange interpolation polynomial of the form where the Lagrange cardinal polynomial functions and are defined as with

At the Chebyshev-Gauss-Lobatto points for , the derivative of , , and with respect to is evaluated as

where is the differentiation matrix of size as indicated in Trefethen [34]. The matrices , , and are defined as where denotes transpose of the matrix.

The derivative of order for , , and w.r.t are approximated as follows:

The derivatives of , , and are evaluated at each point for as where are the entries of the differentiation matrix in the th subinterval. Substituting equations (32)–(37) into equations (22)–(24), we have

The solution at of each subinterval is obtained from the solution at the previous level, redenoting and indices; equations (42)–(44) can be written as

For brevity, equations (45)–(47) can be written as where

4. Results and Discussion

A numerical solution for the boundary layer nanofluid flow along a moving cylinder with thermal radiation, heat generation, and slip effects is studied. The solutions for the system of partial differential equations are presented. The bivariate spectral quasilinearization method (MD-BSQLM) is used to solve the coupled system. To validate the numerical method, a comparison with results in [36] was made, and the results are in acceptable agreement. The finite difference-based method used in [36] shows a slightly more rigorous approach than the one applied in this study.

We present important and interesting results only. The results obtained are presented in Table 1 and Figures 2–18. The discussion of the effect of the Prandtl number and the magnetic parameter is solved for the purpose of validating results obtained in the literature as they have been widely reported. The most interesting and important results include the effect of the Brownian motion , the thermophoresis parameter , heat generation , velocity slip parameter , the temperature exponent , and the radial coordinate . The values , , and give the set of equations (13)–(15) to those of [36].

Figures 2 and 3 show the effect of varying the velocity slip parameter on the velocity and temperature profiles. Increasing result in increasing both the velocity and temperature profiles. Increasing the velocity slip parameter causes the fluid layer in contact with the cylinder surface to slip thereby causing the velocity at the surface to increase. The velocity decrease with increasing boundary layer. When fluid particles move faster on the cylinder surface, this causes more fluid particle interaction causing the temperature to increase with increasing .

Figures 4 and 5 depict the effect of changing the magnetic parameter on velocity profiles and changing the heat generation parameter on the temperature profiles. Increasing the magnetic parameter result in the decrease in velocity profiles, this result agrees with many other results displayed in the literature in [20, 37] and [21]. If the heat generation parameter is increased, this then results in the increase in the temperature profiles; this result has been shown in many other different studies. This shows a strong agreement in the results when the MD-BSQLM is used.

Figure 6 shows the effect of increasing the radiation parameter on the temperature profiles. If the radiation parameter is increased, this has an effect of increasing temperature profiles. A point of “inflexion” is noticed close to the surface of the cylinder. When radiation is applied, it causes a rise in temperature causing the mass particles close to the surface with lower temperature to rise. They then become heated as they move further from the surface causing the temperature to increase further away from this point.

Figures 7 and 8 show the effect of varying the Eckert number and the radial coordinate . Increasing the Eckert number results in the increase in temperature profiles. The trends in the increase in temperature fields are more pronounced at close to the cylinder surface. This is caused by the presence of velocity slip and radiation which increase temperature at the surface. The gradient of the velocity profile at the surface (heat transfer coefficient) can be seen increasing from negative to positive. Increasing the radial coordinate shows a decrease in the concentration profiles. This shows that as we move outward from the cylinder surface, the mass particle concentration decreases. This is consistent with the general settling of mass particles to the cylinder surface due to gravity. The settling is less enhanced due to the fluid flow which tends to carry the particles along with the flow regime.

Figures 9 and 10 show the effect of varying the thermophoresis parameter and the the Brownian motion parameter and on the concentration profiles. If the Brownian motion parameter is increased, this results in decreasing the concentration profile at close to the surface. A reverse effect is noticed further away from the cylinder surface. Increasing the Brownian motion parameter causes the mass particles to move more rapidly in random motion this destabilizes the particle motion. A more enhanced point of inflexion is noticed. Increasing the thermophoresis parameter results in decreasing the concentration profiles at close to the cylinder surface.

In Figure 11, it is shown that if the Schmidt number is increased, this result in the decrease in the concentration profiles. The result is consistent with those in the literature.

The effect of changing parameters of the coupled partial differential equations affects the stability of the numerical solution. Sometimes when parameters exceed certain ranges, the numerical method fails. In this section, we study the effect of changing various parameters on the stability of the numerical solution. This is shown in the graphs from Figures 12–18. The effect of varying the velocity slip on the residual error of the solution of , the Brownian motion on the residual error of , , and , and on , , and . Figure 12 shows the effect of increasing the velocity slip on the residual error of the solution of ; the error does not deteriorate with increasing iterations. The sudden drop in residual error in less than five iterations shows the accuracy of the MD BSQLM. Figure 13 shows the effect of increasing the Brownian motion parameter on the residual error of ; the residual error becomes very small with larger values of the Brownian motion parameter. The error becomes more stable with increasing iterations. In Figures 14 and 15, the errors and improve with increasing Brownian motion parameter. At low values of the parameter convergence is delayed, happens after five iterations. In both cases, the error does not deteriorate with increasing iterations.

In Figures 16–18, it is shown that increasing the radiation parameter does not make the errors , and deteriorate with increasing iterations. Increasing the radiation parameter result in improving the stability of the errors of the MD-BSQLM.

5. Conclusion

In this paper, we have studied the numerical solution of laminar boundary layer flow along with a moving cylinder with thermal radiation, heat generation, and surface slip effect. The coupled partial differential equations of the mass and heat transfer along with a cylinder in motion are solved by the multidomain bivariate spectral quasilinearization method (MD BSQLM). The most important results show that increasing the velocity slip factor results in an enhanced increase in velocity and temperature profiles. Increasing the heat generation parameter increases temperature profiles; if the radiation parameter and the Eckert number are both increased, this would then result in the increase temperature trends. The concentration fields decrease with increasing radial coordinate. Increasing the thermophoresis parameter and Brownian motion both destabilizes the concentration profiles. Increasing the Schmidt number reduce temperature fields. The effect of increasing selected parameters: the velocity slip, Brownian motion, and radiation parameter on all residual errors show that these errors do not deteriorate. This shows that the multidomain bivariate quasilinearization method (MD BSQLM) is very accurate and robust. The method is also easier to implement than the traditional finite difference methods as shown in Section 3. The method also shows very low residual errors only after five iterations. The method adds value to the way in which differential equations arising from fluid dynamics are solved.

Nomenclature

:	Velocity component in direction (m)
:	Velocity component in direction (m)
:	Solute ambient concentration (kg)
:	Surface solute concentration (kg)
:	Fluid ambient temperature ()
:	Fluid temperature at cylinder surface ()
:	Radial coordinate
:	Kinematic viscosity ()
:	Electrical conductivity ()
:	Fluid density ()
:	Magnetic field ()
:	Specific heat capacity (
:	Boltzmann constant ()
:	Mean absorption coefficient ()
:	Thermal diffusivity ()
:	Dynamic viscosity ()
:	Nondimensional ratio of specific heat capacities
:	Brownian diffusion coefficient
:	Thermophoretic diffusion coefficient ()
:	Heat generation constant (m)
:	Radius of cylinder ()
:	Velocity slip factor (m)
:	Free stream velocity ()
:	Surface velocity ()
:	Magnetic parameter (m)
:	Thermal radiation parameter (m)
:	Prandtl number
:	Temperature exponent
:	Eckert number
:	Brownian motion parameter
:	Thermophoresis parameter
:	Nondimensinal heat generation parameter
:	Schmidt number
:	Nondimensional velocity slip parameter

Data Availability

No data was used in this study; calculations were performed from code simulations in Matlab.

Conflicts of Interest

The authors declare that there is no conflict of interest regarding the publication of this article.

Acknowledgments

This work is supported by the National Research Foundation (NRF) of South Africa (Grant Number 129490), the DSI-NRF Centre of Excellence in Mathematical and Statistical Sciences (CoE-MaSS), South Africa, and the Central University of Technology, South Africa.

References

V. Poply, P. Singh, and K. K. Chaudhary, “Analysis of laminar boundary layer flow along a stretching cylinder in the presence of thermal radiation,” WSEAS Transactions on Fluid Mechanics, vol. 8, no. 4, pp. 159–164, 2013.
View at: Google Scholar
S. Shateyi and G. T. Marewo, “A new numerical approach for the laminar boundary layer flow and heat transfer along a stretching cylinder embedded in a porous medium with variable thermal conductivity,” Journal of Applied Mathematics, vol. 2013, Article ID 576453, 7 pages, 2013.
View at: Publisher Site | Google Scholar
R. R. Rangi and N. Ahmad, “Boundary layer flow past a stretching cylinder and heat transfer with variable thermal conductivity,” Applied Mathematics, vol. 3, no. 3, pp. 205–209, 2012.
View at: Publisher Site | Google Scholar
S. Mukhopadhyay, “Analysis of boundary layer flow and heat transfer along a stretching cylinder in a porous medium,” International Scholarly Research Network, vol. 2012, article 704984, pp. 1–7, 2012.
View at: Publisher Site | Google Scholar
E. M. A. Elbashbeshy, T. G. Emam, M. S. El-Azab, and K. M. Abdelgaber, “Laminar boundary layer flow along a stretching cylinder embedded in a porous medium,” International Journal of the Physical Sciences, vol. 7, no. 24, pp. 3067–3072, 2012.
View at: Publisher Site | Google Scholar
H. T. Lin and Y. P. Shih, “Laminar boundary layer heat transfer along a statitic and moving cylinders,” Journal of the Chinese Institute of Engineers, vol. 3, no. 1, pp. 73–79, 1980.
View at: Google Scholar
M. Ali and A. Alabdulkarem, “Laminar mixed convection boundary layer flow induced by a permeable surface stretched with prescribed skin friction boundary conditions,” Advances in Mechanical Engineering, vol. 7, no. 9, 2015.
View at: Publisher Site | Google Scholar
S. U. S. Choi, Enhancing thermal conductivity of fluids with nanoparticles, vol. 231, ASME-Publications-Federation, 1995.
J. Buongiorno, “Convective transport in nanofluids,” Journal of Heat Transfer, vol. 128, no. 3, pp. 240–250, 2006.
View at: Publisher Site | Google Scholar
S. Z. Alamri, R. Ellahi, N. Shehzad, and A. Zeeshan, “Convective radiative plane Poiseuille flow of nanofluid through porous medium with slip: an application of Stefan blowing,” Journal of Molecular Liquids, vol. 273, pp. 292–304, 2019.
View at: Publisher Site | Google Scholar
S. Nadeem, R. Mehmood, and N. S. Akbar, “Oblique stagnation point flow of a Casson-nano fluid towards a stretching surface with heat transfer,” Journal of Computational and Theoretical Nanoscience, vol. 11, no. 6, pp. 1422–1432, 2014.
View at: Publisher Site | Google Scholar
M. Khan, M. Y. Malik, T. Salahuddin, and F. Khan, “Generalized diffusion effects on Maxwell nanofluid stagnation point flow over a stretchable sheet with slip conditions and chemical reaction,” Journal of the Brazilian Society of Mechanical Sciences and Engineering, vol. 41, no. 3, p. 138, 2019.
View at: Publisher Site | Google Scholar
R. Dhanai, P. Rana, and L. Kumar, “Multiple solutions of MHD boundary layer flow and heat transfer behavior of nanofluids induced by a power-law stretching/shrinking permeable sheet with viscous dissipation,” Powder Technology, vol. 273, pp. 62–70, 2015.
View at: Publisher Site | Google Scholar
S. Shah, S. Islam, H. Ayaz, and S. Khan, “Radiative heat and mass transfer analysis of micropolar nanofluid flow of Casson fluid between two rotating parallel plates with effects of Hall current,” Journal of Heat Transfer, vol. 141, no. 2, pp. 022401–022413, 2019.
View at: Publisher Site | Google Scholar
S. Nadeem, M. Riaz khan, and A. U. Khan, “MHD oblique stagnation point flow of nanofluid over an oscillatory stretching/shrinking sheet: existence of dual solutions,” Physica Scripta, vol. 94, no. 7, article 075204, 2019.
View at: Publisher Site | Google Scholar
F. Kamal, K. Zaimi, A. Ishak, and I. Pop, “Stability analysis of MHD stagnation-point flow towards a permeable stretching/shrinking sheet in a nanofluid with chemical reactions effect,” Sains Malaysiana, vol. 48, no. 1, pp. 243–250, 2019.
View at: Publisher Site | Google Scholar
P. Besthapu, R. U. Haq, S. Bandari, and Q. M. Al-Mdallal, “Thermal radiation and slip effects on MHD stagnation point flow of non-Newtonian nanofluid over a convective stretching surface,” Neural Computing and Applications, vol. 31, no. 1, pp. 207–217, 2019.
View at: Publisher Site | Google Scholar
M. A. Sadiq, A. U. Khan, S. Saleem, and S. Nadeem, “Numerical simulation of oscillatory oblique stagnation point flow of a magneto micropolar nanofluid,” RSC Advances, vol. 9, no. 9, pp. 4751–4764, 2019.
View at: Publisher Site | Google Scholar
M. K. Nayak, T. M. Agbaje, S. Mondal, P. Sibanda, and P. G. L. Leach, “Thermodynamic effect in Darchy–Forchheimer nanofluid flow of a single-wall carbon nanotube/multi-wall carbon nanotube suspension due to a stretching/shrinking rotating disk: Buongiorno two-phase model,” Journal of Engineering Mathematics, vol. 120, no. 1, pp. 43–65, 2020.
View at: Publisher Site | Google Scholar
M. A. Yousif, H. F. Ismael, T. Abbas, and R. Ellahi, “Numerical study of momentum and heat transfer of MHD Carreau nanofluid over an exponentially stretched plate with internal heat source/sink and radiation,” Heat Transfer Research, vol. 50, no. 7, pp. 649–658, 2019.
View at: Publisher Site | Google Scholar
S. O. Salawu and H. A. Ogunseye, “Entropy generation of a radiative hydromagnetic Powell-Eyring chemical reaction nanofluid with variable conductivity and electric field loading,” Results in Engineering, vol. 5, article 100072, 2020.
View at: Publisher Site | Google Scholar
B. Kumar and S. Srinivas, “Unsteady hydromagnetic flow of Eyring-Powell nanofluid over an inclined permeable stretching sheet with joule heating and thermal radiation,” Journal of Applied and Computational Mechanics, vol. 6, no. 2, pp. 259–270, 2020.
View at: Google Scholar
P. Sreedevi and P. S. Reddy, “Effect of SWCNTs and MWCNTs Maxwell MHD nanofluid flow between two stretchable rotating disks under convective boundary conditions,” Heat Transfer–Asian Research, vol. 48, no. 8, pp. 4105–4132, 2019.
View at: Publisher Site | Google Scholar
P. Sreedevi, P. S. Reddy, and M. Sheremet, “A comparative study of _Al2O3and TiO2nanofluid flow over a wedge with non-linear thermal radiation,” International Journal of Numerical Methods for Heat & Fluid Flow, vol. 30, no. 3, pp. 1291–1317, 2019.
View at: Publisher Site | Google Scholar
B. J. Gireesha, M. Archana, B. Mahanthesh, and B. C. Prasannakumara, “Exploration of activation energy and binary chemical reaction effects on nano Casson fluid flow with thermal and exponential space-based heat source,” Multidiscipline Modeling in Materials and Structures, vol. 15, no. 1, pp. 227–245, 2019.
View at: Publisher Site | Google Scholar
B. J. Gireesha, M. Umeshaiah, B. C. Prasannakumara, N. S. Shashikumar, and M. Archana, “Impact of nonlinear thermal radiation on magnetohydrodynamic three dimensional boundary layer flow of Jeffrey nanofluid over a nonlinearly permeable stretching sheet,” Physica A: Statistical Mechanics and its Applications, vol. 549, article 124051, 2020.
View at: Publisher Site | Google Scholar
M. Archana, B. J. Gireesha, and B. C. Prasannakumara, “Triple diffusive flow of Casson nanofluid with buoyancy forces and nonlinear thermal radiation over a horizontal plate,” Archives of Thermodynamics, vol. 40, no. 1, pp. 49–69, 2019.
View at: Google Scholar
M. Awais, S. Bilal, M. Y. Malik, and Khalil-ur-Rehman, “Numerical analysis of magnetohydrodynamic Navier’s slip visco nanofluid flow induced by rotating disk with heat source/sink,” Communications in Theoretical Physics, vol. 71, no. 9, p. 1075, 2019.
View at: Publisher Site | Google Scholar
R. E. Bellman and R. E. Kalaba, Quasilinearization and nonlinear boundary-value problems, Elsevier Publishing Company, New York, 1965.
S. S. Motsa, “On the bivariate spectral homotopy analysis method approach for solving nonlinear evolution partial differential equations,” Abstract and Applied Analysis, vol. 2014, Article ID 350529, 8 pages, 2014.
View at: Publisher Site | Google Scholar
S. S. Motsa, V. M. Magagula, and P. Sibanda, “A bivariate Chebyshev spectral collocation quasilinearization method for nonlinear evolution parabolic equations,” The Scientific World Journal, vol. 2014, no. 13, Article ID 581987, 13 pages, 2014.
View at: Publisher Site | Google Scholar
V. M. Magagula, S. S. Motsa, P. Sibanda, and P. G. Dlamini, “On a bivariate spectral relaxation method for unsteady magneto-hydrodynamic flow in porous media,” Springerplus, vol. 5, no. 1, pp. 1–15, 2016.
View at: Publisher Site | Google Scholar
S. S. Motsa, P. G. Dlamini, and M. Khumalo, “Spectral relaxation method and spectral quasilinearization method for solving unsteady boundary layer flow problems,” Advances in Mathematical Physics, vol. 2014, Article ID 341964, 12 pages, 2014.
View at: Publisher Site | Google Scholar
L. N. Trefethen, Spectral Methods in MATLAB, SIAM, Philadelphia, 2000.
C. Canuto, M. Y. Hussaini, A. Quarteroni, and T. A. Zang, Spectral Methods in Fluid Dynamics, Springer-Verlag, Berlin, 1988.
A. Ishak, R. Nazar, and I. Pop, “Local similarity solutions for laminar boundary layer flow along a moving cylinder in a parallel stream,” in Computer Mathematics. ASCM 2007, D. Kapur, Ed., vol. 5081 of Lecture Notes in Computer Science, Springer, Berlin, Heidelberg, 2008.
View at: Publisher Site | Google Scholar
G. S. Seth, B. Kumar, and R. Nandkeolyar, “MHD mixed convection stagnation point flow of a micropolar nanofluid adjacent to stretching sheet: a revised model with successive linearization method,” Journal of Nanofluids, vol. 8, no. 3, pp. 620–630, 2019.
View at: Publisher Site | Google Scholar

Copyright

Copyright © 2021 Titilayo Morenike Agbaje and Gilbert Makanda. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

PDF Download Citation

Download other formats

Order printed copies

Views

616

Downloads

917

Citations