Next Article in Journal
Comparison of Splitting Methods for Deterministic/Stochastic Gross–Pitaevskii Equation
Previous Article in Journal
2D Hybrid Steady-State Magnetic Field Model for Linear Induction Motors
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Overlapping Multi-Domain Spectral Method for Conjugate Problems of Conduction and MHD Free Convection Flow of Nanofluids over Flat Plates

by
Musawenkhosi Mkhatshwa
1,*,
Sandile Motsa
1,2 and
Precious Sibanda
1
1
School of Mathematics, Statistics and Computer Science, University of Kwazulu Natal, Scottsville, Pietermaritzburg 3209, South Africa
2
Department of Mathematics, Faculty of Science and Engineering, University of Eswatini, Private Bag 4, Kwaluseni, Eswatini
*
Author to whom correspondence should be addressed.
Math. Comput. Appl. 2019, 24(3), 75; https://doi.org/10.3390/mca24030075
Submission received: 16 June 2019 / Revised: 17 July 2019 / Accepted: 24 July 2019 / Published: 5 August 2019

Abstract

:
An efficient overlapping multi-domain spectral method is used in the analysis of conjugate problems of heat conduction in solid walls coupled with laminar magnetohydrodynamic (MHD) free convective boundary layer flow of copper (Cu) water and silver (Ag) water nanofluids over vertical and horizontal flat plates. The combined effects of heat generation and thermal radiation on the flow has been analyzed by imposing a magnetic field along the direction of the flow to control the motion of electrically conducting fluid in nanoscale systems. We have assumed that the nanoparticle volume fraction at the wall may be actively controlled. The dimensionless flow equations are solved numerically using an overlapping multi-domain bivariate spectral quasilinearisation method. The effects of relevant parameters on the fluid properties are shown graphically and discussed in detail. Furthermore, the variations of the skin friction coefficient, surface temperature and the rate of heat transfer are shown in graphs and tables. The findings show that the surface temperature is enhanced due to the presence of nanoparticles in the base fluid and the inclusion of the thermal radiation, heat generation and transverse magnetic field in the system. An increase in the nanoparticle volume fraction, heat generation, thermal radiation, and magnetic field parameter enhances the nanofluid velocity and temperature while reducing the heat transfer rate. The results also indicate that the Ag–water nanofluid has higher skin friction and surface temperature than the Cu–water nanofluid, while the opposite behaviour is observed in the case of the rate of heat transfer. The computed numerical results are compared with previously published results and found to be in good agreement.

1. Introduction

Conjugate heat transfer (CHT) is the interaction between the conduction and the buoyancy forced flow of fluid along a solid surface. In numerous applications, the effect of conduction within the solid wall is significant and thus must be taken into account. Such applications include heat exchangers, heaters, nuclear reactors, and pipe insulation systems. In these applications, the analysis of CHT mechanisms, the coupling of the conduction in the solid body and the convection in the fluid surrounding is important [1]. CHT problems, in which the coupled heat transfer processes between conduction and convection mechanisms are considered simultaneously, have been studied by several researchers in the case of Newtonian fluids. For example, Miyamoto et al. [2] reviewed the early theoretical and experimental work of conjugate free convection including the methods and the principal results in the previously obtained solutions of conjugate problems. Miyamoto et al. [2] considered CHT problems of free convection from a vertical plate with a uniform temperature or a uniform heat flux on the outside surface of the plate. Sparrow and Chyu [3] studied CHT problems for a vertical fin with forced convection. Merkin and Pop [4] analyzed CHT over a vertical flat plate using an efficient finite-difference scheme. Pop et al. [5] presented a detailed numerical study of the conjugate mixed convection flow along a vertical flat plate. Luna et al. [6] investigated CHT across a thin horizontal wall separating two fluids at different temperatures numerically and asymptotically. Vasquez and Bula [7] studied the CHT process in cooling a horizontal plate in a steady state condition. Hajmohammadi and Nourazar [8] investigated conjugate forced convection heat transfer from a good conducting horizontal plate with temperature-dependent thermal conductivity. The horizontal plate was heated with uniform heat flux at the lower surface and cooled at the upper surface under laminar forced convection flow. Findings showed that, for a good conducting plate with a finite thickness, the distribution of the conjugate heat flux at the upper surface is significantly affected by the plate thickness. Yu and Lin [9] analyzed conjugate free convection over a vertical and horizontal plate using Keller’s finite-difference method. They proposed the new conjugate parameters and novel dimensionless coordinates to solve the conjugate free convection problem on vertical and horizontal plates. Hsiao [10] analysed the conjugate problems of conduction in solid and free convection in fluid flow using a novel improved formula. The flow equations were solved numerically using the finite difference, Runge–Kutta and Shooting method.
Many studies have been performed on magnetic field and heat generation effects on magnetohydrodynamic (MHD)-conjugate heat transfer. Azim and Chowdhury [11] investigated MHD-conjugate free convection from an isothermal horizontal circular cylinder with Joule heating and heat generation in the presence of a magnetic field. Azim et al. [12] studied the problem of steady CHT through an electrically-conducting fluid for a vertical flat plate with a transverse uniform magnetic field. Kaya [13] investigated mixed convection heat transfer about a thin vertical plate with magneto and CHT effects in a porous medium. Kaya [14] studied the effect of CHT on MHD mixed convection about a vertical slender hollow cylinder. Mamun et al. [15] studied the effects of conduction and viscous dissipation on natural convection flow of an incompressible, viscous and electrically conducting fluid with a transverse magnetic field. Mamun et al. [16] investigated the magnetic field, viscous dissipation and heat generation effects on natural convection flow of incompressible, viscous and electrically conducting fluid along a vertical flat plate with conduction. Hosain and Azim [17] studied the effects of viscous dissipation and heat generation on MHD conjugate free convection flow from an isothermal horizontal circular cylinder when the magnetic field was applied.
In the studies mentioned above, the fluid was assumed to be regular. However, traditional fluids such as water, oil and ethylene glycol might not have enough thermal conductivity to provide the desired efficiency. A good way to overcome this limitation is to add some solid nanoparticles with high thermal conductivity to the fluid. The resulting fluid is a suspension of the solid nanoparticles in the base fluid, which is called nanofluid. The thermal conductivities of nanofluids are believed to be greater than those of the base fluid due to the high thermal conductivity of the nanoparticles. Numerous investigations have been done on the effect of nanoparticles on thermal performance. For example, Choi et al. [18] experimentally studied the effective thermal conductivity of a nano-solid-liquid mixture. Their results revealed that the dispersion of a small amount (<1% by volume) of carbon nanotubes in a liquid increases its thermal conductivity remarkably (nearly 200 % ). Nanofluids have many applications in heat transfer such as microelectronics, fuel cells, Pharmaceutical processes, Hybrid-powered engines, engine cooling vehicles, domestic refrigerator, heat exchanger, nuclear reactor coolant, space technology, and boiler flue gas temperature reduction. Nanoparticles can exist in a variety of types such as metals, metal oxides, carbides, and carbon. The most common types of nanofluids available commercially include Aluminium oxide ( Al 2 O 3 ), Titanium oxide ( Ti O 2 ), Copper (Cu), and Silver (Ag)–water nanofluid [19].
The study of CHT in nanofluids has attracted the interest of many researchers. Jafarian et al. [20] studied CHT in MHD mixed convective flows of nanofluid about a vertical slender hollow Cylinder embedded in a porous medium. Nimmagadda and Venkatasubbaiah [21] analyzed CHT in a micro-channel using novel hybrid nanofluids ( Al 2 O 3 + Ag / Water ). Patrulescu and Grosan [22] studied CHT in a vertical channel filled with a nanofluid adjacent to a heat generating solid domain. Zahan and Alim [23] investigated the problem of developing laminar CHT of copper water nanofluid in a rectangular enclosure. Malvandi et al. [24] studied fluid flow and heat transfer of nanofluids over a flat plate with conjugate heat transfer by including the fluid effects of thermal resistance of the plate in the formulation. Zahan et al. [25] also studied the problem of MHD conjugate natural convection flow in a rectangular frame filled with a copper water nanofluid. Alsabery et al. [26] investigated the conjugate natural convection of A l 2 O 3 -water nanofluid in a square cavity using Buongiorno’s two-phase model. Amongst their findings, they reported that, when the heat conduction is dominated, the heat transfer is increased with the increment of the nanoparticles volume fraction.
High temperatures are required to perform many engineering processes. Nuclear power plants, gas turbines, missiles, satellites, different types of equipment for aircraft, to name a few, can be included in such processes. Accordingly, radiation heat transfer knowledge is very important to design relevant devices. Furthermore, radiation has a significant effect on MHD flow and heat transfer characteristics from an industrial point of view. Industrial applications of thermal radiation include polymer technology, food production, engineering and spinning of fibers and advanced energy conversion in heat transfer at high temperatures. The effect of thermal radiation on MHD convection flow has been investigated by many researchers in the case of regular fluids and nanofluids. Takhar et al. [27] studied the effect of radiation on natural convection flow and heat transfer for a semi-infinite vertical plate with the transverse magnetic field. Emad [28] investigated free convection heat transfer characteristics of an electrically conducting fluid along an isothermal sheet with a transverse magnetic field. In this analysis, the simultaneous effects of buoyancy and radiation with internal heat generation or absorption were considered over the linearly stretched sheet taking into account a uniform free stream of constant velocity and temperature. El-Naby et al. [29] investigated natural convection unsteady flow over a semi-finite vertical plate with variable temperature, radiation, and transverse magnetic field. Ali et al. [30] investigated thermal radiation effects on the time-independent hydromagnetic forced convective flow of an electrically conducting and heat generating-absorbing fluid over a non-isothermal wedge. Mbeledogu et al. [31] obtained the perturbation solutions of the problem formed by the simultaneous action of buoyancy and transverse magnetic field on free convection flow of compressible Boussinesq fluid past a moving vertical plate. The viscosity and thermal conductivity of the fluid were a function of temperature and the radiative flux was confirmed using the Rosseland approximation. Ali et al. [32] analyzed the effect of thermal radiation and heat generation on viscous Joule heating MHD-conjugate heat transfer along a vertical flat plate. Their results showed that thermal radiation, viscous Joule heating and internal heat generation in the presence of conduction effects have a significant effect on MHD natural convection flow and thermal fields. In the case of nanofluids, Elazem et al. [33] considered the effect of radiation on the steady MHD flow and heat transfer of C u –water and A g –water nanofluids flow over a stretching sheet. Raju et al. [34] investigated the influence of the magnetic field, radiation, and non-uniform heat source/sink on Cu–Ethyline glycol and Ag–Ethyline glycol nanofluids flow over a moving vertical plate in a porous medium. The results from these studies reveal that, as thermal radiation increases, the rate of energy transported to the fluid increases, consequently an increase in temperature occurs.
The thermal radiation effect on the MHD-conjugate flow of nanofluids over flat plates with internal heat generation can be important in many industrial and theoretical applications. However, the literature review shows that no significant study investigated the combined effects of thermal radiation and heat generation on MHD-conjugate heat transfer flow on natural convection in a nanofluid filled enclosure. The objective of this study is to extend the work of Yu and Lin [9] by analyzing the conjugate heat transfer in MHD free convective flow of Cu–water and Ag–water nanofluids along the vertical and horizontal plates with internal heat generation and thermal radiation. This study is theoretical and can have practical significance in designing and operation of plate heat exchangers. It is worth mentioning that the problem considered has applications in industries such as flat fins and cooling of electronic boards due to the inclusion of nanoparticles. The flow is subject to a uniform magnetic field imposed along the direction of the flow. We further demonstrate the application of efficient overlapping multi-domain bivariate spectral quasilinearisation method in solving a nonlinear system of partial differential equations (PDEs) modeling CHT problems. This method is more accurate than the non-overlapping multi-domain bivariate spectral quasilinearisation method (MD-BSQLM) [35]. The non-overlapping MD-BSQLM applies the multi-domain technique only in the time interval. However, the method considered in the present work applies the multi-domain technique in both space and time intervals. In addition to that, the method uses the overlapping multi-domain technique in the space interval. The overlapping grid strategy can improve the accuracy of spectral collocation based methods. The accuracy improvement is achieved through making the coefficient matrix in the matrix equation (resulting from the collocation process) less dense. This means that the coefficient matrix will be sparse. The sparsity of matrices caused by overlapping sub-domains can help to minimize the storage of large matrices and make it easy to perform matrix-vector multiplications. This is because there will be a lot of multiplication by zero which reduces the computational time and enables the matrices to be stored efficiently. Since the method combines the bivariate spectral quasilinearisation method [36], non-overlapping and overlapping multi-domain technique, for reference purposes, we shall refer to the method as the overlapping multi-domain bivariate spectral quasilinearisation method (OMD-BSQLM). The use of spectral collocation-based methods such as the OMD-BSQLM for solving systems of PDEs can be a most promising tool in the study of conjugate heat transfer problems. From the literature review, several studies [2,9] concluded that it is very difficult to obtain analytical solutions of conjugate heat transfer problems due to the matching conditions at the solid–fluid interface. These studies proposed the use of numerical methods such as finite difference schemes as the most promising procedures for performing this matching. However, the finite difference methods have a lot of limitations when compared to spectral methods. Spectral methods are highly accurate and more efficient than traditional methods such as the finite difference methods [37]. When applied to problems with smooth solutions, they use few grid points and require minimal computational time to generate accurate solutions, thus they are better than traditional methods. The spectral method algorithm is easy to implement in scientific computing software. To establish the accuracy of the OMD-BSQLM, certain limiting solutions of the flow equations are studied.

2. Mathematical Formulation

Let us consider the viscous, steady, incompressible, electrically conducting and free convection flow of nanofluid over a vertical flat plate and a horizontal flat plate of finite length l and thickness b . The thickness of the plates is assumed to be smaller than the length. The base fluid is water and nanoparticles (Cu and Ag) are in thermal equilibrium with no slip between them. The thermophysical properties of the base fluid and different nanoparticles are shown in Table 1. It is assumed that the left side of the vertical plate and the lower side of the horizontal plate are maintained at the constant temperature T b , such that T b > T , where T is the temperature of the ambient nanofluid. Heat is transferred by conduction from the outside surface of the solid plates coupled with the free convection in the nanofluid, while the axial heat conduction in the plates is neglected. A uniform magnetic field B ( x ) is imposed along the direction of flow. The applied transverse magnetic field can be chosen in its special form as B ( x ) = B 0 α f 1 / 2 x 1 , where B 0 is the steady strength of the magnetic field towards the y-axis. It is assumed that the induced magnetic field and the external electric field are negligible. Thermal radiation and internal heat generation terms are included in the energy equation. The geometry and coordinate system for the vertical and horizontal flat plates are shown in Figure 1. With the above assumptions, equations of the conservation of mass, momentum, and energy are given by
u x + v y = 0 ,
u u x + v u y = 1 ρ n f p x + ν n f 2 u y 2 + g β n f ( T T ) sin φ σ n f B 2 ( x ) ρ n f u ,
0 = 1 ρ n f p y + g β n f ( T T ) cos φ ,
u T x + v T y = k n f ( ρ C p ) n f 2 T y 2 + Q 0 ( ρ C p ) n f ( T T ) 1 ( ρ C p ) n f q r y ,
where u and v are the velocity components in the x- and y- directions, p is the pressure, g is the acceleration due to gravity, T is the fluid temperature near the plate, q r is the radiative heat flux, Q 0 is the rate of heat generation, σ n f is the electrical conductivity, ν n f is the kinematic viscosity, μ n f is the dynamic viscosity, ρ n f is the effective density, α n f is the thermal diffusivity, k n f is the effective thermal conductivity, β n f is the thermal expansion coefficient and ( ρ C p ) n f is the heat capacity of the nanofluid. The term Q 0 ( T T ) represents the amount of heat generated or absorbed per unit volume, where Q 0 is a constant which may be either positive for a heat sink or negative for a heat source. The radiative heat flux q r with Rosseland approximation has the form q r = 4 σ * 3 k * T 4 y , where σ * is the Stefan–Boltzmann constant and k * is the mean absorption coefficient. The temperature differences within the flow are assumed to be sufficiently small such that T 4 may be expressed as a linear function of temperature. Expanding T 4 using Taylor series and neglecting higher order terms yields T 4 4 T 3 T 3 T 4 . The nanofluid constants are defined as [38,39,40]
ν n f = μ n f ρ n f , μ n f = μ f ( 1 ϕ ) 2 . 5 , k n f k f = ( k s + 2 k f ) 2 ϕ ( k f k s ) ( k s + 2 k f ) + ϕ ( k f k s ) , σ n f σ f = 1 + 3 σ s σ f 1 ϕ σ s σ f + 2 σ s σ f 1 ϕ , ρ n f = ( 1 ϕ ) ρ f + ϕ ρ s , ( ρ C p ) n f = ( 1 ϕ ) ( ρ C p ) f + ϕ ( ρ C p ) s , β n f = ( 1 ϕ ) β f + ϕ β s , α n f = k n f ( ρ C p ) n f ,
where ϕ is the solid volume fraction of nanoparticles, β is the thermal expansion, subscripts f , s and n f denote fluid, solid and nanofluid, respectively. For the vertical plate, the angle φ is π / 2 and p / x , and p / y are both equal to zero. For the horizontal plate, the angle φ to the horizontal is equal to zero. In formulating Equations (1)–(4), viscous dissipation and compression work have been neglected. Moreover, the physical properties of the fluid are assumed to be constant except for the density variation that induces a buoyancy force. The boundary conditions for Equations (1)–(4) are
u = 0 , v = 0 , at y = 0 ,
u 0 , p 0 , T T as y .

2.1. Dimensionless Equations for the Vertical Plate

Yu and Lin [9] developed the following non-dimensional variables:
ψ ( x , y ) = α f λ f ( ξ , η ) , η ( x , y ) = ( y / x ) λ , ξ ( x ) = 1 + σ R a t / ( σ R a h ) 4 / 5 1 , θ ( ξ , η ) = T T T b T ξ 1 ,
where R a t = g β ( T b T ) x 3 / α f ν is the Rayleigh number, R a h = g β ( q h x / k r ) x 3 / α f ν is the Reyleigh number for a plate with constant wall flux q h = k s ( T b T ) / b , λ = ( σ R a t ) 1 + ( σ R a h ) 4 / 5 1 / 4 , σ = P r / ( 1 + P r ) , ψ ( x , y ) is the stream function defined by u = ψ / y and v = ψ / x , η ( x , y ) and ξ ( x ) are the dimensionless coordinates, f ( ξ , η ) is the dimensionless stream function and θ ( ξ , η ) is the dimensionless temperature. By using Equation (8), Equations (1)–(4) along with the boundary conditions (6) and (7) are reduced to the following two-point boundary value problem:
P r f + ϕ 1 16 ξ 20 f f 6 ξ 10 f 2 M 2 ϕ 2 f + ϕ 3 ( 1 + P r ) θ = ϕ 1 5 ξ ( 1 ξ ) f f ξ f f ξ ,
1 + k f k n f R d θ + ϕ 4 16 ξ 20 f θ 1 ξ 5 f θ + k f k n f Q ξ θ = ϕ 4 5 ξ ( 1 ξ ) f θ ξ θ f ξ ,
where M 2 = σ f B 0 2 ρ f λ 2 is the magnetic field parameter, P r = ν f α f is the Prandtl number, R d = 16 σ * T 3 3 k * k f is the radiation parameter and Q = Q 0 x 2 k f λ 2 is the heat generation parameter. The nanoparticle volume fractions ϕ 1 , ϕ 2 , ϕ 3 and ϕ 4 depend on the thermal properties of the nanofluid and are defined as
ϕ 1 = [ 1 ϕ ] 2 . 5 1 ϕ + ϕ ρ s ρ f , ϕ 2 = 1 + 3 ( σ s / σ f 1 ) ϕ ( σ s / σ f + 2 ) ( σ s / σ f 1 ) ϕ 1 ( 1 ϕ ) + ϕ ( ρ s ρ f ) , ϕ 3 = ( 1 ϕ ) + ϕ ( β s / β f ) , ϕ 4 = k s + 2 k f + ϕ ( k f k s ) k s + 2 k f 2 ϕ ( k f k s ) ( 1 ϕ ) + ϕ ( ρ C p ) s ( ρ C p ) f .
The corresponding boundary conditions in dimensionless form are
f ( ξ , 0 ) = 0 , f ( ξ , 0 ) = 0 , ξ θ ( ξ , 0 ) ( 1 ξ ) 5 / 4 θ ( ξ , 0 ) = 1 , f ( ξ , ) = 0 , θ ( ξ , ) = 0 .

2.2. Dimensionless Equations for the Horizontal Plate

Using the following non-dimensional variables in [9], namely
ψ ( x , y ) = α f λ f ( ξ , η ) , η ( x , y ) = ( y / x ) λ , ξ ( x ) = 1 + σ R a t / ( σ R a h ) 5 / 6 1 . 5 ,
θ ( ξ , η ) = T T T b T ξ 1 , ω ( ξ , η ) = σ p x 2 / ρ α f ν λ 4 ,
where λ = ( σ R a t ) 1 + ( σ R a h ) 5 / 6 1 . 5 and ω ( ξ , η ) is the dimensionless pressure. Equations (1)–(4) along with their boundary conditions (6) and (7) for the horizontal plate are reduced to
P r f + ϕ 1 10 ξ 15 f f 5 2 ξ 15 f 2 ϕ 2 M 2 f + ( 1 ϕ ) 2 . 5 15 ( 1 + P r ) [ ( 5 + ξ ) η ω ( 10 4 ξ ) ω ]
= ϕ 1 3 ξ ( 1 ξ ) f f ξ f f ξ + ( 1 + P r ) ω ξ ,
ω = θ ,
1 + k f k n f R d θ + ϕ 4 10 ξ 15 f θ 1 ξ 3 f θ + k f k n f Q ξ θ = ϕ 4 3 ξ ( 1 ξ ) f θ ξ θ f ξ .
The corresponding boundary conditions in dimensionless form are
f ( ξ , 0 ) = 0 , f ( ξ , 0 ) = 0 , ξ θ ( ξ , 0 ) ( 1 ξ ) 6 / 5 θ ( ξ , 0 ) = 1 ,
f ( ξ , ) = 0 , θ ( ξ , ) = 0 , ω ( ξ , ) = 0 .

3. Solution Procedure

In this section, we describe the application of the OMD-BSQLM to find numerical solutions of the transformed nonlinear PDEs. The method uses the overlapping multi-domain technique, Chebyshev–Gauss–Lobatto grid points [41,42], and the quasilinearisation method [43], together with spectral collocation on approximate functions defined as bivariate Lagrange interpolation polynomials. The multi-domain approach divides the time interval into non-overlapping sub-intervals and the space interval into overlapping sub-intervals. The quasilinearisation technique helps to linearise the nonlinear PDEs. The spectral collocation method is applied independently both in space and time variables in the linearized equations. In order to apply the OMD-BSQLM, the time interval ξ [ 0 , ξ F ] is decomposed into q non-overlapping sub-intervals defined as
J υ = ( ξ υ 1 , ξ υ ) , υ = 1 , 2 , 3 , , q , with 0 = ξ 0 < ξ 1 < ξ 2 < < ξ q 1 < ξ q = ξ F ,
For the semi-finite space domain [ 0 , ) , a truncated grid [ 0 , η ] is used. We choose a finite value of η that is large enough such that the flow properties at η resemble those at . The truncated space interval [ 0 , η ] is decomposed into p overlapping sub-intervals of length L , denoted by
I μ = [ η 0 μ , η N η μ ] , μ = 1 , 2 , 3 , , p ,
where each I μ interval is further discretized into N η + 1 collocation points. Without loss of generality, we will consider that each subinterval has the same length given by
L = η p + 1 2 ( 1 p ) ( 1 cos π N η )
for the overlap to be possible, and the same number of collocation points ( N η + 1 ) is used in each subinterval. In the domain decomposition scheme, we use overlapping subintervals I μ , where the first two points of the interval I μ + 1 coincide with the last two points of the interval I μ , that is, η 0 1 = 0 , η N η p = η , η N η 1 μ = η 0 μ + 1 and η N η μ = η 1 μ + 1 . The non-overlapping and overlapping multi-domain grids are shown in Figure 2 and Figure 3, respectively.

3.1. Numerical Solution for the Vertical Plate

Applying the quasilinearisation method in each sub-interval to Equations (9) and (10) gives the following system of linear PDEs:
α 1 , 3 , r ( 1 , μ , υ ) 3 f r + 1 ( μ , υ ) η 3 + α 1 , 2 , r ( 1 , μ , υ ) 2 f r + 1 ( μ , υ ) η 2 + α 1 , 1 , r ( 1 , μ , υ ) f r + 1 ( μ , υ ) η + α 1 , 0 , r ( 1 , μ , υ ) f r + 1 ( μ , υ ) + α 2 , 0 , r ( 1 , μ , υ ) θ r + 1 ( μ , υ )
+ γ 1 , r ( 1 , μ , υ ) ξ f r + 1 ( μ , υ ) η + β 1 , r ( 1 , μ , υ ) f r + 1 ( μ , υ ) ξ = R 1 , r ( μ , υ ) , α 2 , 2 , r ( 2 , μ , υ ) 2 θ r + 1 ( μ , υ ) η 2 + α 2 , 1 , r ( 2 , μ , υ ) θ r + 1 ( μ , υ ) η + α 2 , 0 , r ( 2 , μ , υ ) θ r + 1 ( μ , υ ) + α 1 , 1 , r ( 2 , μ , υ ) f r + 1 ( μ , υ ) η + α 1 , 0 , r ( 2 , μ , υ ) f r + 1 ( μ , υ )
+ β 2 , r ( 2 , μ , υ ) θ r + 1 ( μ , υ ) ξ + β 1 , r ( 2 , μ , υ ) f r + 1 ( μ , υ ) ξ = R 2 , r ( μ , υ ) ,
where the variable coefficients are given by
α 1 , 3 , r ( 1 , μ , υ ) = P r , α 1 , 2 , r ( 1 , μ , υ ) = ϕ 1 ( 16 ξ ) 20 f r ( μ , υ ) + ϕ 1 ξ ( 1 ξ ) 5 f r ( μ , υ ) ξ , α 1 , 0 , r ( 1 , μ , υ ) = ϕ 1 ( 16 ξ ) 20 2 f r ( μ , υ ) η 2 , α 1 , 1 , r ( 1 , μ , υ ) = ϕ 1 ( 6 ξ ) 5 f r ( μ , υ ) η M ϕ 1 ϕ 2 ϕ 1 ξ ( 1 ξ ) 5 ξ f r ( μ , υ ) η , α 2 , 0 , r ( 1 , μ , υ ) = ϕ 1 ϕ 3 ( 1 + P r ) , α 2 , 2 , r ( 2 , μ , υ ) = 1 + 1 + k f k n f R d , α 2 , 0 , r ( 2 , μ , υ ) = k k k n f Q ξ ϕ 4 ( 1 ξ ) 5 f r ( μ , υ ) η , α 2 , 1 , r ( 2 , μ , υ ) = ϕ 4 ( 16 ξ ) 20 f r ( μ , υ ) + ϕ 4 ξ ( 1 ξ ) 5 f r ( μ , υ ) ξ , α 1 , 1 , r ( 2 , μ , υ ) = ϕ 4 ( 1 ξ ) 5 θ r ( μ , υ ) ϕ 4 ξ ( 1 ξ ) 5 θ r ( μ , υ ) ξ , α 1 , 0 , r ( 2 , μ , υ ) = ϕ 4 ( 16 ξ ) 20 θ r ( μ , υ ) η , γ 1 , r ( 1 , μ , υ ) = ϕ 1 ξ ( 1 ξ ) 5 f r ( μ , υ ) η , β 1 , r ( 1 , μ , υ ) = ϕ 1 ξ ( 1 ξ ) 5 2 f r ( μ , υ ) η 2 , β 2 , r ( 2 , μ , υ ) = ϕ 4 ξ ( 1 ξ ) 5 f r ( μ , υ ) η , β 1 , r ( 2 , μ , υ ) = ϕ 4 ξ ( 1 ξ ) 5 θ r ( μ , υ ) η , R 1 , r ( μ , υ ) = ϕ 1 ( 16 ξ ) 20 f r ( μ , υ ) 2 f r ( μ , υ ) η 2 ϕ 1 ( 6 ξ ) 10 f r ( μ , υ ) η 2 ϕ 1 ξ ( 1 ξ ) 5 f r ( μ , υ ) η ξ f r ( μ , υ ) η + ϕ 1 ξ ( 1 ξ ) 5 2 f r ( μ , υ ) η 2 f r ( μ , υ ) ξ , R 2 , r ( μ , υ ) = ϕ 4 ( 16 ξ ) 20 f r ( μ , υ ) θ r ( μ , υ ) η ϕ 4 ( 1 ξ ) 5 θ r ( μ , υ ) f r ( μ , υ ) η ϕ 4 ξ ( 1 ξ ) 5 f r ( μ , υ ) η θ r ( μ , υ ) ξ + ϕ 4 ξ ( 1 ξ ) 5 θ r ( μ , υ ) η f r ( μ , υ ) ξ ,
subject to boundary conditions
f r + 1 ( μ , υ ) ( ξ , 0 ) = 0 , f r + 1 ( μ , υ ) η ( ξ , 0 ) = 0 , ξ θ r + 1 ( μ , υ ) ( ξ , 0 ) ( 1 ξ ) 5 / 4 θ r + 1 ( μ , υ ) η ( ξ , 0 ) = 1 , f r + 1 ( μ , υ ) η ( ξ , ) = 0 , θ r + 1 ( μ , υ ) ( ξ , ) = 0 .
The constants r and r + 1 denote previous and current iterations, respectively. The system of linear PDEs (23) and (24) is discretized using the spectral collocation method in both η and ξ directions. Before applying the spectral method on the sub-intervals, the time interval ξ [ ξ υ 1 , ξ υ ] is transformed to τ [ 1 , 1 ] using the linear transformation
ξ j υ = 1 2 ( ξ υ ξ υ 1 ) τ j + 1 2 ( ξ υ + ξ υ 1 ) , τ j = cos π j N ξ ,
and the space region η [ η 0 μ , η N η μ ] is transformed to z [ 1 , 1 ] using the linear transformation
η i μ = L 2 ( z i + 1 ) , z i = cos π i N η .
We assume that, at each sub-interval, the required solution, say f ( η , ξ ) , can be approximated by a bivariate Lagrange interpolation polynomial of the form
f ( μ , υ ) ( η , ξ ) i = 0 N η j = 0 N ξ f ( μ , υ ) ( z i , τ j ) L i ( z ) L j ( τ ) ,
for μ = 1 , 2 , 3 , , p and υ = 1 , 2 , 3 , , q . The bivariate interpolation polynomial interpolates f ( μ , υ ) ( z , τ ) at selected points ( z i , τ j ) in both z and τ directions, for i = 0 , 1 , 2 , , N η and j = 0 , 1 , 2 , , N ξ . The function L i ( z ) and L j ( τ ) are the well known characteristic Lagrange cardinal polynomial based on the Chebyshev–Gauss–Lobatto points. The required solution for θ ( η , ξ ) can be approximated in a similar manner. The solution procedure requires that the derivatives of L i ( z ) and L j ( τ ) with respect to z and τ , respectively be evaluated at the Chebyshev–Gauss–Lobatto grid points. The derivatives of f ( μ , υ ) ( η , ξ ) with respect to η and ξ at the Chebyshev–Gauss–Lobatto points ( z k , τ i ) , are computed as
f ( μ , υ ) η ( z k , τ i ) = ω N η j = 0 N ξ f ( μ , υ ) ( z ω , τ j ) d L ω ( z k ) d z L j ( τ i )
= ω = 0 N η D k , ω ( μ ) f ( μ , υ ) ( z ω , τ i ) = D ( μ ) F i ( μ , υ ) , f ( μ , υ ) ξ ( z k , τ i ) = ω N η j = 0 N ξ f ( μ , υ ) ( z ω , τ j ) L ω ( z k ) d L j ( τ i ) d τ
= j = 0 N ξ d i , j f ( μ , υ ) ( z k , τ j ) = j = 0 N ξ d i , j F j ( μ , υ ) ,
where d i , j = d L j ( τ i ) d τ is the ith and jth entry of the standard first derivative Chebyshev–Gauss–Lobatto based differentiation matrix d = [ d i , j ] , for i , j = 0 , 1 , 2 , 3 , , N ξ , of size ( N ξ + 1 ) × ( N ξ + 1 ) , D k , ω ( μ ) = 2 η N η μ η 0 μ D k , ω with D k , ω = d L ω ( z k ) d z being the kth and ω th entries of the standard first derivative Chebyshev–Gauss–Lobatto differentiation matrix of size ( M + 1 ) × ( M + 1 ) , where M = N η + ( N η 1 ) ( p 1 ) is the total number of collocation points over a single domain [ 1 , 1 ] . In general, to find an sth order derivative with respect to η , we have
s f ( μ , υ ) η s ( z k , τ i ) = ω = 0 N η D k , ω ( μ ) s f ( μ , υ ) ( z ω , τ i ) = D ( μ ) s F i ( μ , υ ) .
The vector F i ( μ , υ ) is defined as
F i ( μ , υ ) = f ( μ , υ ) ( z 0 ( μ ) , τ i ( υ ) ) , f ( μ , υ ) ( z 1 ( μ ) , τ i ( υ ) ) , f ( μ , υ ) ( z 2 ( μ ) , τ i ( υ ) ) , , f ( μ , υ ) ( z N η ( μ ) , τ i ( υ ) ) T ,
where T denotes the matrix transpose. The derivatives d and D ( μ ) s are scaled by multiplying by the factors Λ = 2 ξ υ ξ υ 1 and Ω s = 2 η N η μ η 0 μ s = 2 L s , respectively. The space and time derivatives of θ at each sub-interval can be transformed to discrete matrix form in a similar manner. Applying the spectral collocation method by evaluating Equations (23) and (24) at the collocation points and making use of the derivative matrices as well as incorporating the initial condition which corresponds to ξ N ξ = 1 gives
A 1 , 1 ( μ , υ ) F i , r + 1 ( μ , υ ) + A 1 , 2 ( μ , υ ) Θ i , r + 1 ( μ , υ ) + γ 1 , r ( 1 , μ , υ ) j = 0 N ξ 1 d i , j D ( μ ) F j ( μ , υ ) + β 1 , r ( 1 , μ , υ ) j = 0 N ξ 1 d i , j F j ( μ , υ ) = K 1 , i ( μ , υ ) ,
A 2 , 1 ( μ , υ ) F i , r + 1 ( μ , υ ) + A 2 , 2 ( μ , υ ) Θ i , r + 1 ( μ , υ ) + β 1 , r ( 2 , μ , υ ) j = 0 N ξ 1 d i , j F j ( μ , υ ) + β 2 , r ( 2 , μ , υ ) j = 0 N ξ 1 d i , j Θ j ( μ , υ ) = K 2 , i ( μ , υ ) ,
where
K 1 , i ( μ , υ ) = R 1 , i ( μ , υ ) γ 1 , r ( 1 , μ , υ ) d i , N ζ D ( μ ) F N ξ ( μ , υ ) β 1 , r ( 1 , μ , υ ) d i , N ξ F N ξ ( μ , υ ) , K 2 , i ( μ , υ ) = R 2 , i ( μ , υ ) β 1 , r ( 2 , μ , υ ) d i , N ζ F N ξ ( μ , υ ) β 2 , r ( 2 , μ , υ ) d i , N ξ Θ N ξ ( μ , υ )
For i = 0 , 1 , 2 , N ξ , Equations (34) and (35) form an N ξ ( M + 1 ) × N ξ ( M + 1 ) matrix system
Mca 24 00075 i001
where
A i , i ( 1 , 1 , p , υ ) = α 1 , 3 , r ( 1 , μ , υ ) D ( μ ) 3 + α 1 , 2 , r ( 1 , μ , υ ) D ( μ ) 2 + α 1 , 1 , r ( 1 , μ , υ ) D ( μ ) + α 1 , 0 , r ( 1 , μ , υ ) + β 1 , r ( 1 , μ , υ ) d i , i I + γ 1 , r ( 1 , μ , υ ) d i , i D ( μ ) , A i , i ( 1 , 2 , p , υ ) = α 2 , 0 , r ( 1 , μ , υ ) I , A i , i ( 2 , 1 , p , υ ) = α 1 , 1 , r ( 2 , μ , υ ) D ( μ ) + α 1 , 0 , r ( 2 , μ , υ ) + β 1 , r ( 2 , μ , υ ) d i , i I , A i , i ( 2 , 2 , p , υ ) = α 2 , 2 , r ( 2 , μ , υ ) D ( μ ) 2 + α 2 , 1 , r ( 2 , μ , υ ) D ( μ ) + α 2 , 0 , r ( 2 , μ , υ ) + β 2 , r ( 2 , μ , υ ) d i , i I , when i = j
and
A i , j ( 1 , 1 , p , υ ) = β 1 , r ( 1 , μ , υ ) d i , j I + γ 1 , r ( 1 , μ , υ ) d i , j D ( μ ) , A i , j ( 1 , 2 , p , υ ) = 0 , A i , j ( 2 , 1 , p , υ ) = β 1 , r ( 2 , μ , υ ) d i , j I , A i , j ( 2 , 2 , p , υ ) = β 2 , r ( 2 , μ , υ ) d i , j I , when i j .
The vectors F i , r + 1 ( μ , υ ) , and Θ i , r + 1 ( μ , υ ) denote the values of f and θ approximated at the collocation points, and I is the standard ( M + 1 ) × ( M + 1 ) identity matrix. Starting from suitable initial guesses, the numerical solution for f ( η , ξ ) and θ ( η , ξ ) are obtained by solving matrix Equation (36) iteratively for r = 1 , 2 , , σ , where σ is the number of iterations to be used.

3.2. Numerical Solution for the Horizontal Plate

Applying a quasilinearisation method in each subinterval to the system of nonlinear PDEs (15)–(17) gives the following system of linear PDEs:
α 1 , 3 , r ( 1 , μ , υ ) 3 f r + 1 ( μ , υ ) η 3 + α 1 , 2 , r ( 1 , μ , υ ) 2 f r + 1 ( μ , υ ) η 2 + α 1 , 1 , r ( 1 , μ , υ ) f r + 1 ( μ , υ ) η + α 1 , 0 , r ( 1 , μ , υ ) f r + 1 ( μ , υ ) + α 2 , 1 , r ( 1 , μ , υ ) ω r + 1 ( μ , υ ) η ,
+ α 2 , 0 , r ( 1 , μ , υ ) ω r + 1 ( μ , υ ) + γ 1 , 1 , r ( 1 , μ , υ ) ξ f r + 1 ( μ , υ ) η + β 1 , 1 , r ( 1 , μ , υ ) f r + 1 ( μ , υ ) ξ + β 2 , 1 , r ( 1 , μ , υ ) ω r + 1 ( μ , υ ) ξ = R 1 , r ( μ , υ ) ,
α 2 , 1 , r ( 2 , μ , υ ) ω r + 1 ( μ , υ ) η + α 3 , 0 , r ( 2 , μ , υ ) θ r + 1 ( μ , υ ) = 0 , α 3 , 2 , r ( 3 , μ , υ ) 2 θ r + 1 ( μ , υ ) η 2 + α 3 , 1 , r ( 3 , μ , υ ) θ r + 1 ( μ , υ ) η + α 3 , 0 , r ( 3 , μ , υ ) θ r + 1 ( μ , υ ) + α 1 , 1 , r ( 3 , μ , υ ) f r + 1 ( μ , υ ) η + α 1 , 0 , r ( 3 , μ , υ ) f r + 1 ( μ , υ )
+ β 3 , 1 , r ( 3 , μ , υ ) θ r + 1 ( μ , υ ) ξ + β 1 , 1 , r ( 3 , μ , υ ) f r + 1 ( μ , υ ) ξ = R 3 , r ( μ , υ ) ,
where the variable coefficients are given by
α 1 , 3 , r ( 1 , μ , υ ) = P r , α 1 , 2 , r ( 1 , μ , υ ) = ϕ 1 ( 10 ξ ) 15 f r ( μ , υ ) + ϕ 1 ξ ( 1 ξ ) 3 f r ( μ , υ ) ξ , α 1 , 0 , r ( 1 , μ , υ ) = ϕ 1 ( 10 ξ ) 15 2 f r ( μ , υ ) η 2 , α 1 , 1 , r ( 1 , μ , υ ) = ϕ 1 ( 10 4 ξ ) 15 f r ( μ , υ ) η M 2 ϕ 1 ϕ 2 ϕ 1 ξ ( 1 ξ ) 3 ξ f r ( μ , υ ) η , α 2 , 1 , r ( 1 , μ , υ ) = ( 1 ϕ ) 2 . 5 15 ( 1 + P r ) ( 10 4 ξ ) , ( 1 ϕ ) 2 . 5 15 ( 1 + P r ) ( 5 + ξ ) η , α 2 , 1 , r ( 2 , μ , υ ) = 1 , α 3 , 0 , r ( 2 , μ , υ ) = 1 , α 3 , 2 , r ( 3 , μ , υ ) = 1 + 1 + k f k n f R d , α 3 , 1 , r ( 3 , μ , υ ) = ϕ 4 ( 10 ξ ) 15 f r ( μ , υ ) + ϕ 4 ξ ( 1 ξ ) 3 f r ( μ , υ ) ξ , α 3 , 0 , r ( 3 , μ , υ ) = k f k n f Q ξ ϕ 4 ( 1 ξ ) 3 f r ( μ , υ ) η , α 1 , 1 , r ( 3 , μ , υ ) = ϕ 4 ( 1 ξ ) 3 θ r ( μ , υ ) ϕ 4 ξ ( 1 ξ ) 3 θ r ( μ , υ ) ξ , α 1 , 0 , r ( 3 , μ , υ ) = ϕ 4 ( 10 ξ ) 15 θ r ( μ , υ ) η , γ 1 , 1 , r ( 1 , μ , υ ) = ϕ 1 ξ ( 1 ξ ) 3 f r ( μ , υ ) η , β 1 , 1 , r ( 1 , μ , υ ) = ϕ 1 ξ ( 1 ξ ) 3 2 f r ( μ , υ ) η 2 , β 2 , 1 , r ( 1 , μ , υ ) = ϕ 1 ξ ( 1 ξ ) 3 ( 1 + P r ) , β 3 , 1 , r ( 3 , μ , υ ) = ϕ 4 ξ ( 1 ξ ) 3 f r ( μ , υ ) η , β 1 , 1 , r ( 3 , μ , υ ) = ϕ 4 ξ ( 1 ξ ) 5 θ r ( μ , υ ) η , R 1 , r ( μ , υ ) = ϕ 1 ( 10 ξ ) 15 f r ( μ , υ ) 2 f r ( μ , υ ) η 2 ϕ 1 ( 5 2 ξ ) 15 f r ( μ , υ ) η 2 ϕ 1 ξ ( 1 ξ ) 3 f r ( μ , υ ) η ξ f r ( μ , υ ) η + ϕ 1 ξ ( 1 ξ ) 3 2 f r ( μ , υ ) η 2 f r ( μ , υ ) ξ , R 2 , r ( μ , υ ) = 0 , R 3 , r ( μ , υ ) = ϕ 4 ( 10 ξ ) 15 f r ( μ , υ ) θ r ( μ , υ ) η ϕ 4 ( 1 ξ ) 3 θ r ( μ , υ ) f r ( μ , υ ) η ϕ 4 ξ ( 1 ξ ) 5 f r ( μ , υ ) η θ r ( μ , υ ) ξ + ϕ 4 ξ ( 1 ξ ) 5 θ r ( μ , υ ) η f r ( μ , υ ) ξ ,
subject to boundary conditions
f r + 1 ( μ , υ ) ( ξ , 0 ) = 0 , f r + 1 ( μ , υ ) η ( ξ , 0 ) = 0 , ξ θ r + 1 ( μ , υ ) ( ξ , 0 ) ( 1 ξ ) 6 / 5 θ r + 1 ( μ , υ ) η ( ξ , 0 ) = 1 , f r + 1 ( μ , υ ) η ( ξ , ) = 0 , ω r + 1 ( μ , υ ) ( ξ , ) = 0 θ r + 1 ( μ , υ ) ( ξ , ) = 0 .
We apply the Chebyshev spectral collocation method that uses bivariate Lagrange interpolation polynomials as basic functions as in the vertical plate. Thus, evaluating Equations (39) and (41) at the collocation points and making use of the derivative matrices as well as incorporating the initial conditions which corresponds to ξ N ξ , we obtain
A 1 , 1 ( μ , υ ) F i , r + 1 ( μ , υ ) + A 1 , 2 ( μ , υ ) Ω i , r + 1 ( μ , υ ) + A 1 , 3 ( μ , υ ) Θ i , r + 1 ( μ , υ ) + γ 1 , r ( 1 , 1 , μ , υ ) j = 0 N ξ d i , j D ( μ ) F j ( μ , υ ) + β 1 , 1 , r ( 1 , μ , υ ) j = 0 N ξ d i , j F j ( μ , υ ) ,
+ β 2 , 1 , r ( 1 , μ , υ ) j = 0 N ξ d i , j Ω j ( μ , υ ) = R 1 , i ( μ , υ ) ,
A 2 , 1 ( μ , υ ) F i , r + 1 ( μ , υ ) + A 2 , 2 ( μ , υ ) Ω i , r + 1 ( μ , υ ) + A 2 , 3 ( μ , υ ) Θ i , r + 1 ( μ , υ ) = R 2 , i ( μ , υ ) ,
A 3 , 1 ( μ , υ ) F i . r + 1 ( μ , υ ) + A 3 , 2 ( μ , υ ) Ω i , r + 1 ( μ , υ ) + A 3 , 3 ( μ , υ ) Θ i , r + 1 ( μ , υ ) + β 1 , 1 , r ( 3 , μ , υ ) j = 0 N ξ d i , j F j ( μ , υ ) + β 3 , 1 , r ( 3 , μ , υ ) j = 0 N ξ d i , j Θ j ( μ , υ ) = R 3 , i ( μ , υ ) .
The vectors F i , r + 1 ( μ , υ ) , Ω i , r + 1 ( μ , υ ) and Θ i , r + 1 ( μ , υ ) denote the values of f , ω and θ approximated at the collocation points. Imposing boundary conditions for i = 0 , 1 , 2 , 3 , , N ξ 1 , Equations (43)–(45) can be expressed as a matrix system of size N ξ ( M + 1 ) × N ξ ( M + 1 ) as in the previous subsection.

4. Results and Discussion

The transformed nonlinear PDEs for the vertical and horizontal plates were solved numerically using the OMD-BSQLM for Cu–water and Ag–water nanofluids. Numerical computations are carried out using P r = 0.7 [9], M = 0.5 and Q = 0.01 [12,16]. However, the parametric values of the radiation parameter and nanoparticle volume fraction were chosen as R d = 0.6 and ϕ = 0.3 . All of these values are treated the same in the entire study except the varied values in respective figures. The space domain η was truncated to η = 15 . The numerical results were generated using N ξ = 5 , N η = 20 collocation points. The number of sub-intervals in both space and time are taken as p = q = 5 . In order to obtain a clear understanding of the physics of the problem, a parametric study was undertaken to determine the impact of the different physical parameters on the fluid properties and flow characteristics.
To determine the accuracy of our numerical results, the local skin friction coefficient and the surface temperature are compared with the non-overlapping MD-BSQLM and published results by Yi and Lin [9] in Table 2. The table gives a comparison of the OMD-BSQLM results when ξ = M = ϕ = 0 for different values of the Prandtl number P r . It is observed that, for increasing values of the Prandtl number, the results are in good agreement with values in the literature and those obtained using the non-overlapping MD-BSQLM. Hence, the use of the present method is justified. It is also noted that the OMD-BSQLM can give accurate results with a minimal number of grid points compared to the non-overlapping MD-BSQLM. Table 3 presents results for the local skin friction f ( ξ , 0 ) , surface temperature θ ( ξ , 0 ) and heat transfer rate θ ( ξ , 0 ) for varying values of the dimensionless streamwise coordinate ξ and different nanofluids. The table shows clearly that the skin friction, interfacial temperature, and heat transfer rate decrease with increasing values of ξ . This is due to the increase of the momentum boundary layer thickness and thermal boundary layer thickness.
Figure 4, Figure 5, Figure 6 and Figure 7 depict the effects of nanoparticle volume fraction, thermal radiation, heat generation and magnetic field parameter on the velocity profiles for both Ag and Cu nanofluids. It is observed from the figures that Ag–water nanofluid shows better enhancement in the velocity profiles than Cu–water nanofluid. This is because the viscosity of the Ag–water nanofluid is higher compared to that of Cu–water nanofluid. The effect of using different types of nanofluids is more significant in the vertical plate than on the horizontal plate. Figure 4 shows the influence of the magnetic parameter on the dimensionless velocity. It is noted that the velocity is higher near the wall and lower far from the wall for hydrodynamic flows ( M = 0 ). The opposite trend is observed for the hydromagnetic flows ( M 0 ). Moreover, increasing the magnetic parameter reduces the velocity distribution near the wall. The magnetic parameter is known to represent the Lorentz force that opposes the flow. The peak velocity decreases with the increasing values of the magnetic parameter due to the retarding effect in the boundary layer region. As a result, the separation of the boundary layer occurs earlier since the momentum boundary layer becomes thick. These findings concur with results reported by Mamun et al. [12] and Azim et al. [16] in regular fluids.
Figure 5 shows the effect of nanoparticle volume fraction on the velocity profiles. It is seen that the flow velocity increases around the vertical and horizontal plates with an increase in the nanoparticle volume fraction. For both the vertical and horizontal plates, it is clear that the flow velocity is significantly low for the conventional fluid ( ϕ = 0 ) than for nanofluids ( ϕ 0 ) . As expected, for the conventional fluid, there is no change in velocity profiles for both plates. However, as the volume fraction of nanoparticles increases, the velocity distribution also increases. This is due to an increase in the momentum boundary layer thickness which is attributed to adding nanoparticles to the base fluid. The nanoparticles enhance the velocity profiles due to the higher thermal conductivity of nanofluids. For the horizontal plate, we also observe that, near the wall, the momentum boundary layer thickness decreases as the volume fraction of silver particles increases and away from the wall, the boundary layer thickness increases. Figure 6 is presented to show the effect of the heat generation parameter on the dimensionless velocity. It is observed that more heat is generated within the boundary as the heat generation parameter increases and, consequently, the fluid velocity increases as well. The increase in velocity is consistent with the physical consequence as the internal energy generation resulted from the heat generation increases the buoyancy forces, which in turn enhance more flow along both the vertical and horizontal plates. The effect of the thermal radiation parameter on the velocity distribution is shown in Figure 7 for both Ag–water and Cu–water nanofluids. We observe that the velocity increases within the boundary layer thickness as the thermal radiation parameter increases. Radiation accelerates the fluid motion, thus enhancing the velocity of nanofluids.
Figure 8, Figure 9, Figure 10 and Figure 11 show the influence of thermal radiation, nanoparticle volume fraction, heat generation and magnetic field parameter on the temperature profiles for both nanofluids. The temperature distribution in the case of Ag–water nanofluid is relatively higher than in the case of Cu–water nanofluid. This is because the thermal conductivity of silver nanoparticles is higher than that of Copper nanoparticles. The effect of using different types of nanofluids is more clear on the horizontal plate than for the vertical plate. Figure 8 depicts the influence of the magnetic field parameter on the temperature profiles. The figure shows that the magnetic field enhances the thickness of the thermal boundary layer, thus increasing the temperature profiles. The effect of nanoparticle volume fraction on the temperature profiles is shown in Figure 9. For the vertical and horizontal plates, the thermal boundary layer thickness is enhanced when the nanoparticle volume fraction increases. Physically, increasing the nanoparticle volume fraction causes an increase in the thermal conductivity of the nanofluid, which in turn enhances the boundary layer thickness and an augmentation in the temperature profiles. Similar results were reported by Shahzad et al. [44]. It is worth mentioning that the temperature is significantly higher in the case of nanofluids than in the regular fluid ( ϕ = 0 ). This is due to the presence of high conductive silver and copper nanoparticles.
Figure 10 illustrates the effect of the heat generation parameter on the temperature profiles. It is seen that the thermal boundary layer is enhanced when the heat generation parameter increases. The energy resulted from internal heat generation increases the temperature of the fluid within the boundary and increases the motion of the fluid. The influence of thermal radiation on the temperature profiles is shown in Figure 11. The figure depicts that an increase in the thermal radiation parameter improves the temperature profiles. As the temperature increases with increasing radiation parameter, the thickness of the thermal boundary layer is enhanced. The larger values of the amount of k * k n f 4 σ * T 3 in the radiation parameter indicate dominance in the thermal radiation over conduction. Thus, there is a large amount of radiative heat energy being poured into the system. The fluid within the boundary layer absorbs imitated heat from the heated plate because of the radiation effect. The radiated heat ultimately increases the temperature of the fluid. The greater values of thermal radiation parameter generate higher temperature and, consequently, the fluid motion is accelerated. Similar results were obtained by Ali et al. [32] in the case of regular fluid.
Figure 12, Figure 13, Figure 14 and Figure 15 present the variation of the local skin friction at different values of the thermal radiation, heat generation, nanoparticle volume fraction and magnetic parameter for the Ag–water and Cu–water nanofluids. The skin friction coefficient is higher in the Ag–water nanofluid compared to the Cu–water nanofluid. Hence, the Ag–water nanofluid gives a high drag force in opposition to the flow compared to the Cu–water nanofluid. For the vertical plate, the skin friction is higher for the Ag–water nanofluid throughout the surface. However, for the horizontal plate, the skin friction is higher for a Ag–water nanofluid close to the wall and higher for a Cu–water nanofluid far from the wall. The effect of the magnetic parameter on the skin friction coefficient is shown in Figure 12. The figure shows that, when the magnetic parameter increases, the skin friction coefficient decreases. The magnetic force that opposes the flow decreases the shear stress at the wall, thus reducing the skin friction coefficient.
The behaviour of the skin friction coefficient against the streamwise coordinate ξ for different values of the nanoparticle volume fraction is plotted in Figure 13. The figure shows that an increase in the nanoparticle volume fraction causes a decrease in the skin friction at the plates. In Figure 14, the impact of the heat generation parameter on the local skin friction is exhibited. The figure reflects that the skin friction factor increases with increasing heat generation parameter. As mentioned earlier, increasing the heat generation parameter accelerates the flow and generates greater buoyancy force and thus increases the skin friction coefficient. The influence of the thermal radiation parameter on the skin friction coefficient is shown in Figure 15. The increase in the fluid motion due to thermal radiation enhances the shear stress at the wall which in turn causes an increase in the skin friction coefficient.
Figure 16, Figure 17, Figure 18 and Figure 19 show the effects of the thermal radiation, heat generation, nanoparticle volume fraction and magnetic field parameter on the rate of heat transfer for both Ag and Cu nanofluids. The rate of heat transfer is observed to be higher in the Cu–water nanofluid than in the Ag–water nanofluid for both vertical and horizontal plates. An increase in the magnetic parameter reduces the rate of heat transfer as seen in Figure 16. The increasing magnetic field parameter enhances the thermal boundary layer thickness and consequently the heat transfer rate decreases due to an increase in the magnetic field strength. In addition, the rate of heat transfer depends on the gradient of temperature and, as the temperature gradient decreases with increasing values of the magnetic parameter, the heat transfer rate also decreases.
Figure 17 depicts the impact of the nanoparticle volume fraction on the skin friction for the different nanofluids. Increasing the nanoparticle volume fraction enhances the thermal conductivity of the nanofluids, which reduces the thermal boundary layer thickness and the temperature gradient at the wall as observed from the figure. The influence of the heat generation on the heat transfer rate is depicted in Figure 18. The figure shows that the heat transfer rate decreases with increasing heat generation parameter. Since higher values of the heat generation parameter create a hot layer of fluid near the surface which results in the temperature of the fluid to exceed the surface temperature, accordingly, the rate of heat transfer from the surface decreases. Figure 19 depicts that increasing values of the thermal radiation parameter enhances the fluid interfacial temperature, which in turn makes the flow of the heat rate slower.
Figure 20, Figure 21, Figure 22 and Figure 23 display the influence of the thermal radiation, heat generation, nanoparticle volume fraction and magnetic field parameter on the surface temperature for the Ag–water and Cu–water nanofluids. The surface temperature is noted to be higher in the Ag–water nanofluid than in the Cu–water nanofluid for both the vertical and horizontal plates. In Figure 20, we observe that, when the magnetic field is applied in the system, the surface temperature is enhanced for both the vertical and horizontal plates. As the magnetic field increases, the surface temperature is enhanced. The interaction between the magnetic field and the fluid motion increases the temperature of the fluid within the boundary layer which in turn increases the thermal boundary layer thickness as well as the surface temperature. Figure 21 shows that adding nanoparticles to the fluid enhances the surface temperature since the surface temperature increases with increasing nanoparticle volume fraction. The surface temperature increases with increasing values of the heat generation parameter as observed in Figure 22. This is because the temperature within the boundary layer increases for increasing heat generation parameter and thus enhances the surface temperature. Figure 23 shows that increasing the thermal radiation parameter also enhances the surface temperature.

5. Conclusions

The multi-domain bivariate spectral quasilinearisation method was used to analyze conjugate heat transfer in MHD free convection flow of copper water and silver water nanofluids over vertical and horizontal plates. The comparison with previously published results was performed and the results were in good agreement. The effects of nanofluids, heat generation, thermal radiation, nanoparticle volume fraction and magnetic field parameter on the fluid properties and flow characteristics were discussed appropriately with numerical computations. The results obtained in the present study can have practical importance in various problems such as ablation or perspiration cooling problems. The results of the paper are of engineering interest where heat transfer processes are controlled in polymer processing and nuclear reactor cooling systems. and designing and operation of plate heat exchangers. From the obtained results and discussion, the following conclusions can be drawn:
  • The Ag–water nanofluid has higher velocity and temperature profiles, skin friction coefficient, and surface temperature than the Cu–water nanofluid. However, the reverse is true for the rate of heat transfer.
  • Heat generation, thermal radiation, nanoparticle volume fraction and magnetic field parameter enhance the velocity of the nanofluid far from the wall. However, an increase in the magnetic field parameter significantly decreases the velocity of the nanofluid near the wall.
  • Increasing the heat generation, thermal radiation, nanoparticle volume fraction and magnetic field parameter improves the temperature distribution and the surface temperature while reducing the rate of heat transfer.
  • The overlapping multi-domain bivariate spectral quasilinearisation method holds great potential for solving highly nonlinear conjugate heat transfer problems since the method gives accurate results using a minimal number of grid points.

Author Contributions

Conceptualization, M.M.; Data curation, M.M.; Methodology, M.M. and S.M.; Supervision, S.M. and P.S.; Validation, S.M. and P.S.; Writing—original draft, M.M.; Writing—review and editing, M.M., S.M. and P.S.

Funding

This research received no external funding

Acknowledgments

The authors would like to thank the anonymous reviewers for their valuable comments and suggestions to improve the quality of the paper.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
B ( x ) External uniform magnetic field
B 0 Magnetic strength
pPressure
R a Rayleigh number
gGravitational acceleration
kThermal conductivity (W/m K)
C p Specific heat capacity
TFluid temperature (K or C )
q h Heat flux
fDimensionless stream function
( u , v ) Velocity component in Cartesian coordinate
T b Constant temperature
T Ambient temperature
Q 0 Rate of heat generation
q r Radiative heat flux
MMagnetic field parameter
P r Prandtl number
R d Radiation parameter
QHeat generation parameter
Greek Symbols
η Scaled boundary layer coordinate
ξ Streamwise coordinate
σ Electrical conductivity (S m 1 )
α Thermal diffusivity m 2 s 1
μ Dynamic viscosity kg m 1 s 1
θ Dimensionless temperature
ϕ Nanoparticle volume fraction parameter
ψ Stream function m 2 s 1
ρ Density of the fluid ( Kg/m 3 )
β Thermal expansion coefficient
ν Kinematic viscosity m 2 s 1
Subscripts
n f Nanofluid phase
fFluid phase
sSolid phase

References

  1. Chang, C.L. Buoyancy and wall conduction effects on forced convection of micropolar fluid flow along a vertical slender hollow circular cylinder. Int. J. Heat Mass Transf. 2006, 49, 4932–4942. [Google Scholar] [CrossRef]
  2. Miyamoto, M.; Sumikawa, J.; Akiyoshi, T.; Nakamura, T. Effects of axial heat conduction in a vertical flat plate on free convection heat transfer. Int. J. Heat Mass Transf. 1980, 23, 1545–1553. [Google Scholar] [CrossRef]
  3. Sparrow, E.M.; Chyu, M.K. Conjugated forced convection-conduction analysis of heat transfer in a plate fin. ASME J. Heat Transf. 1982, 104, 204–206. [Google Scholar] [CrossRef]
  4. Merkin, J.H.; Pop, I. A note on the free convection boundary layer on a horizontal circular cylinder with constant heat flux. Wärme Stoffübertrag. 1988, 22, 79–81. [Google Scholar] [CrossRef]
  5. Pop, I.; Lesnic, D.; Ingham, D.B. Conjugate mixed convection on a vertical surface in a porous medium. Int. J. Heat Mass Transf. 1995, 38, 1517–1525. [Google Scholar] [CrossRef]
  6. Luna, E.; Trevino, C.; Higuera, F.J. Conjugate natural convection heat transfer between two fluids seperated by a horizontal wall: Steady-state analysis. Heat Mass Transf. 1996, 31, 353–358. [Google Scholar] [CrossRef]
  7. Vasquez, R.S.; Bula, A.J. Uncoupling the Conjugate Heat Transfer Problem in a Horizontal Plate Under the Influence of a Laminar Flow, Heat Transfer. In Proceedings of the ASME 2003 International Mechanical Engineering Congress and Exposition, Washington, DC, USA, 15–21 November 2003; Volume 1, pp. 11–18. [Google Scholar]
  8. Hajmohammadi, M.R.; Nourazar, S.S. Conjugate Forced Convection Heat Transfer From a Heated Flat Plate of Finite Thickness and Temperature-Dependent Thermal Conductivity. Heat Transf. Eng. 2014, 35, 863–874. [Google Scholar] [CrossRef]
  9. Yu, W.S.; Lin, H.T. Conjugate problems of conduction and free convection on vertical and horizontal flat plates. Int. J. Heat Mass Transf. 1993, 36, 1303–1313. [Google Scholar] [CrossRef]
  10. Hsiao, K.L. Conjugate heat transfer for free convection along vertical plate fin. J. Therm. Sci. 2010, 19, 337–345. [Google Scholar] [CrossRef]
  11. Azim, N.H.M.; Chowdhury, M.K. MHD-conjugate free convection from an isothermal horizontal circular cylinder with Joule heating and heat generation. J. Comput. Meth. Phys. 2013, 2013, 1–11. [Google Scholar] [CrossRef]
  12. Azim, M.; Mamun, A.; Rahman, M. Viscous Joule heating MHD-conjugate heat transfer for a vertical flat plate in the presence of heat generation. Int. J. Commun. Heat. Mass. Transf. 2010, 37, 666–674. [Google Scholar] [CrossRef]
  13. Kaya, A. Effects of Conjugate Heat Transfer on Steady MHD Mixed Convective Heat Transfer Flow over a Thin Vertical Plate Embedded in a Porous Medium with High Porosity. Math. Prob. Eng. 2012, 2012, 1–19. [Google Scholar] [CrossRef] [Green Version]
  14. Kaya, A. The effect of conjugate heat transfer on MHD mixed convection about a vertical slender hollow cylinder. Commun. Nonlinear Sci. Numer. Simul. 2012, 16, 1905–1916. [Google Scholar] [CrossRef]
  15. Al-Mamun, A.; Azim, N.H.M.A.; Maleque, M.A. Combined effect of conduction and viscous dissipation on MHD free convection flow along a vertical flat plate. J. Naval Archit. Mar. Eng. 1970, 4, 87–98. [Google Scholar] [CrossRef]
  16. Mamun, A.A.; Chowdhury, Z.R.; Azim, M.A.; Molla, M.M. MHD-conjugate heat transfer analysis for a vertical flat plate in the presence of viscous dissipation and heat generation. Int. Commun. Heat Mass Transf. 2008, 35, 1275–1280. [Google Scholar] [CrossRef]
  17. Azim, N.H.M.A. Effects of Viscous Dissipation and Heat Generation on MHD Conjugate Free Convection Flow from an Isothermal Horizontal Circular Cylinder. SOP Trans. Appl. Phys. 2014, 2014, 1–11. [Google Scholar] [CrossRef]
  18. Choi, S.U.S.; Zhang, Z.G.; Yu, W.; Lockwood, F.E.; Grulke, E.A. Anomalous thermal conductivity enhancement in nanotube suspensions. Appl. Phys. Lett. 2001, 79, 2252–2254. [Google Scholar] [CrossRef]
  19. Malvandi, A.; Hedayati, F.; Nobari, M.R.H. An HAM analysis of stagnation-point flow of nanofluid over a porous stretching sheet with heat generation. J. Appl. FluidMech. 2014, 7, 135–145. [Google Scholar]
  20. Jafarian, B.; Hajipour, M.; Khademi, R. Conjugate Heat Transfer of MHD non-Darcy Mixed Convection Flow of a Nanofluid over a Vertical Slender Hollow Cylinder Embedded in Porous Media. Trans. Phenom. Nano Micro Scales 2016, 4, 1–10. [Google Scholar]
  21. Nimmagadda, R.; Venkatasubbaiah, K. Conjugate heat transfer analysis of micro-channel using novel hybrid nanofluids (Al2O3 + Ag/Water). Eur. J. Mech. B Fluids. 2015, 52, 19–27. [Google Scholar] [CrossRef]
  22. Patrulescu, F.; Grosan, T. Conjugate heat transfer in a vertical channel filled with a nanofluid adjacent to a heat generating solid domain. Rev. Anal. Numer. Theor. Approx. 2010, 36, 141–149. [Google Scholar]
  23. Zahan, I.; Alim, M.A. Effect of Conjugate heat transfer on flow of nanofluid in a rectangular enclosure. Int. J. Heat Tech. 2008, 36, 397–405. [Google Scholar] [CrossRef]
  24. Malvandi, A.; Hedayati, F.; Ganji, D.D. Fluid and Heat Transfer of Nanofluids over a flat plate with conjugate heat transfer. Trans. Phenom. Nano. Micro Scales 2014, 2, 108–117. [Google Scholar]
  25. Zahan, I.; Nasrin, R.; Alim, M.A. MHD effect on conjugate heat transfer in a nanofluid filled rectangular enclosure. Int. J. Petrochem. Sci. Eng. 2018, 3, 114–123. [Google Scholar] [CrossRef]
  26. Alsabery, A.I.; Sheremet, M.A.; Chamkha, A.J.; Hashim, I. Conjugate natural convection of Al2O3-water nanofluid in a square cavity with a concentric solid insert using Buongiorno’s two-phase model. Int. J. Mech. Sci. 2018, 136, 200–219. [Google Scholar] [CrossRef]
  27. Takher, H.S.; Gorla, R.S.R.; Soundalgekar, V.M. Short communication radiation effects on MHD free convection flow of a gas past a semi-finite vertical plate. J. Numer. Meth. Heat Fluid Flow. 1996, 6, 77–83. [Google Scholar] [CrossRef]
  28. AboEldahab, E.M. Radiation effects on heat transfer in an electrically conducting fluid at a stretching surface with a uniform free stream. J. Phys. D Appl. Phys. 2000, 33, 3180. [Google Scholar] [CrossRef]
  29. El-Naby, A.; Elbarbary, E.M.E.; Abdelazem, N.Y. Finite difference solution of radiation effects on MHD unsteady free convection flow over vertical plate with variable surface temperature. J. Appl. Math. 2003, 2, 65–86. [Google Scholar] [CrossRef]
  30. Chamkha, A.J.; Mujtaba, M.; Quadri, A.; Issa, C. Thermal radiation effects on MHD forced convection flow adjacent to a non-isothermal wedge in the presence of a heat source/sink. Heat Mass Transf. 2003, 3, 305–312. [Google Scholar] [CrossRef]
  31. Mbeledogu, I.U.; Amakiri, A.R.C.; Ogulu, A. Unsteady MHD free convective flow of a compressible fluid past a moving vertical plate in the presence of radiative heat transfer. Int. J. Heat Mass Transf. 2007, 50, 1668–1674. [Google Scholar] [CrossRef]
  32. Ali, M.M.; Mamun, A.A.; Maleque, M.A. Radiation and heat generation effects on viscous Joule heating MHD-conjugate heat transfer for a vertical flat plate. Can. J. Phys. 2014, 92, 509–521. [Google Scholar] [CrossRef]
  33. Elazem, N.Y.A.; Ebaid, A.; Aly, E.H. Radiation effects of MHD on Cu-water and Ag-water nanofluids flow over a stretching sheet: Numerical study. J. Appl. Comput. Math. 2015, 4, 235. [Google Scholar] [CrossRef]
  34. Raju, C.S.K.; Babu, M.J.; Sandeep, N.; Krishna, P.M. Influence of non-uniform heat source/sink on MHD nanofluid flow over a moving vertical plate in porous medium. Chem. Process Eng. Res. 2015, 6, 31–42. [Google Scholar]
  35. Magagula, V.M.; Motsa, S.S.; Sibanda, P. A Multi-domain Bivariate Pseudospectral Method for Evolution Equations. Int. J. Comput. Meth. 2017, 14, 1750041. [Google Scholar] [CrossRef]
  36. Motsa, S.S.; Magagula, V.M.; Sibanda, P. A Bivariate Chebyshev Spectral Collocation Quasilinearization Method for Nonlinear Evolution Parabolic Equations. Sci. World J. 2014, 2014, 1–13. [Google Scholar] [CrossRef] [Green Version]
  37. Boyd, J.P. Chebyshev and Fourier Spectral Methods; Springer: Berlin, Germany; New York, NY, USA, 1989. [Google Scholar]
  38. Hayat, T.; Kiran, A.; Imtiaz, M.; Alsaedi, A. Hydromagnetic mixed convection flow of copper and silver water nanofluids due to a curved stretching sheet. Results Phys. 2016, 6, 904–910. [Google Scholar] [CrossRef] [Green Version]
  39. Raza, J.; Rohni, A.M.; Omar, Z. MHD flow and heat transfer of Cu-water nanofluid in a semi porous channel with stretching walls. Int. J. Heat Mass Transf. 2016, 103, 336–340. [Google Scholar] [CrossRef]
  40. Prasad, P.D.; Kumar, R.V.M.S.S.; Varma, S.V.K. Heat and mass transfer analysis for the MHD flow of nanofluid with radiation absorption. Ain Shaims Eng. J. 2018, 9, 801–813. [Google Scholar] [CrossRef] [Green Version]
  41. Canuto, C.; Hussaini, M.Y.; Quarteroni, A.; Zang, T.A. Spectral Methods in Fluid Dynamics; Springer: Berlin, Germany, 1988. [Google Scholar]
  42. Trefethen, L.N. Spectral Methods in Matlab; SIAM: Philadelphia, PA, USA, 2000. [Google Scholar]
  43. Bellman, R.E.; Kalaba, R.E. Quasilinearization and Nonlinear Boundary Value Problems; RAND Corporation: Santa Monica, CA, USA, 1965. [Google Scholar]
  44. Shahzad, F.; Sagheer, M.; Hussain, S. Numerical simulation of MHD Jeffrey nanofluid flow and heat transfer over a stretching sheet considering Joule heating and viscous dissipation. AIP Adv. 2018, 8, 065316. [Google Scholar] [CrossRef]
Figure 1. Physical model and coordinate system. (a) vertical plate; (b) horizontal plate.
Figure 1. Physical model and coordinate system. (a) vertical plate; (b) horizontal plate.
Mca 24 00075 g001
Figure 2. Non-overlapping grid ( ξ -domain).
Figure 2. Non-overlapping grid ( ξ -domain).
Mca 24 00075 g002
Figure 3. Overlapping grid ( η -domain).
Figure 3. Overlapping grid ( η -domain).
Mca 24 00075 g003
Figure 4. Velocity profiles for various values M. (a) vertical plate; (b) horizontal plate.
Figure 4. Velocity profiles for various values M. (a) vertical plate; (b) horizontal plate.
Mca 24 00075 g004
Figure 5. Velocity profiles for various values ϕ . (a) vertical plate; (b) horizontal plate.
Figure 5. Velocity profiles for various values ϕ . (a) vertical plate; (b) horizontal plate.
Mca 24 00075 g005
Figure 6. Velocity profiles for various values Q. (a) vertical plate; (b) horizontal plate.
Figure 6. Velocity profiles for various values Q. (a) vertical plate; (b) horizontal plate.
Mca 24 00075 g006
Figure 7. Velocity profiles for various values R d . (a) vertical plate; (b) horizontal plate.
Figure 7. Velocity profiles for various values R d . (a) vertical plate; (b) horizontal plate.
Mca 24 00075 g007
Figure 8. Temperature profiles for various values of M. (a) vertical plate; (b) horizontal plate.
Figure 8. Temperature profiles for various values of M. (a) vertical plate; (b) horizontal plate.
Mca 24 00075 g008
Figure 9. Temperature profiles for various values ϕ . (a) vertical plate; (b) horizontal plate.
Figure 9. Temperature profiles for various values ϕ . (a) vertical plate; (b) horizontal plate.
Mca 24 00075 g009
Figure 10. Temperature profiles for various values of Q. (a) vertical plate; (b) horizontal plate.
Figure 10. Temperature profiles for various values of Q. (a) vertical plate; (b) horizontal plate.
Mca 24 00075 g010
Figure 11. Temperature profiles for various values R d . (a) vertical plate; (b) horizontal plate.
Figure 11. Temperature profiles for various values R d . (a) vertical plate; (b) horizontal plate.
Mca 24 00075 g011
Figure 12. Skin friction coefficient for various values of M. (a) vertical plate; (b) horizontal plate.
Figure 12. Skin friction coefficient for various values of M. (a) vertical plate; (b) horizontal plate.
Mca 24 00075 g012
Figure 13. Skin friction coefficient for various values of ϕ . (a) vertical plate; (b) horizontal plate.
Figure 13. Skin friction coefficient for various values of ϕ . (a) vertical plate; (b) horizontal plate.
Mca 24 00075 g013
Figure 14. Skin friction for various values of Q. (a) vertical plate; (b) horizontal plate.
Figure 14. Skin friction for various values of Q. (a) vertical plate; (b) horizontal plate.
Mca 24 00075 g014
Figure 15. Skin friction coefficient for various values of R d . (a) vertical plate; (b) horizontal plate.
Figure 15. Skin friction coefficient for various values of R d . (a) vertical plate; (b) horizontal plate.
Mca 24 00075 g015
Figure 16. Heat transfer rate for various values of M. (a) vertical plate; (b) horizontal plate.
Figure 16. Heat transfer rate for various values of M. (a) vertical plate; (b) horizontal plate.
Mca 24 00075 g016
Figure 17. Heat transfer rate for various values of ϕ . (a) vertical plate; (b) horizontal plate.
Figure 17. Heat transfer rate for various values of ϕ . (a) vertical plate; (b) horizontal plate.
Mca 24 00075 g017
Figure 18. Heat transfer rate for various values of Q. (a) vertical plate; (b) horizontal plate.
Figure 18. Heat transfer rate for various values of Q. (a) vertical plate; (b) horizontal plate.
Mca 24 00075 g018
Figure 19. Heat transfer rate for various values of R d . (a) vertical plate; (b) horizontal plate.
Figure 19. Heat transfer rate for various values of R d . (a) vertical plate; (b) horizontal plate.
Mca 24 00075 g019
Figure 20. Surface temperature for various values of M. (a) vertical plate; (b) horizontal plate.
Figure 20. Surface temperature for various values of M. (a) vertical plate; (b) horizontal plate.
Mca 24 00075 g020
Figure 21. Surface temperature for various values of ϕ . (a) vertical plate; (b) horizontal plate.
Figure 21. Surface temperature for various values of ϕ . (a) vertical plate; (b) horizontal plate.
Mca 24 00075 g021
Figure 22. Surface temperature for various values of Q. (a) vertical plate; (b) horizontal plate.
Figure 22. Surface temperature for various values of Q. (a) vertical plate; (b) horizontal plate.
Mca 24 00075 g022
Figure 23. Surface temperature for various values of R d . (a) vertical plate; (b) horizontal plate.
Figure 23. Surface temperature for various values of R d . (a) vertical plate; (b) horizontal plate.
Mca 24 00075 g023
Table 1. Thermophysical properties of the base fluid and the nanoparticles [38].
Table 1. Thermophysical properties of the base fluid and the nanoparticles [38].
Base FluidNanoparticles
Physical PropertiesWaterCopper (Cu)Silver (Ag)
C p ( J / kgK ) 4179385235
ρ ( Kg / m 3 ) 997.1893310,500
k ( W / mK ) 0.613401429
σ ( S m 1 ) 0.05 5.96 × 10 7 6.3 × 10 7
β × 10 5 ( K 1 ) 21 1.67 1.89
Table 2. Comparison of the OMD-BSQLM results with MD-BSQLM, Yi and Lin [9] for f ( 0 , 0 ) and θ ( 0 , 0 ) at different values of P r when ξ = M = ϕ = 0 .
Table 2. Comparison of the OMD-BSQLM results with MD-BSQLM, Yi and Lin [9] for f ( 0 , 0 ) and θ ( 0 , 0 ) at different values of P r when ξ = M = ϕ = 0 .
Yi and Lin [9] MD-BSQLM OMD-BSQLM
η Pr f ( 0 , 0 ) θ ( 0 , 0 ) f ( 0 , 0 ) θ ( 0 , 0 ) N η f ( 0 , 0 ) θ ( 0 , 0 ) N η
Vertical plate
120.00154.7451.334554.74635211.334435610054.74635211.334435620
120.0116.9291.375916.92955161.375856210016.92955161.375856220
120.15.25021.48241.25023421.48239991001.25023421.482399920
150.72.31231.61322.31234801.61291661002.31234801.612916620
1571.57481.65201.57435191.65189401001.57435191.651894020
Horizontal plate
120.00147.1661.225847.20486731.225770310047.20486731.225770320
120.0114.5491.272014.55012641.272014910014.55012641.272014920
120.14.54241.39444.54233691.39437241004.54233691.394372420
150.72.02051.55832.07573561.55304461002.07573561.553044620
1571.36221.64101.36185151.64134641001.36185151.641346420
Table 3. OMD-BSQLM results for the skin friction coefficient, heat transfer rate and surface temperature at different values of ξ when P r = 0.7 , ϕ = 0.3 , M = 0.5 , Q = 0.01 and R d = 0.6 .
Table 3. OMD-BSQLM results for the skin friction coefficient, heat transfer rate and surface temperature at different values of ξ when P r = 0.7 , ϕ = 0.3 , M = 0.5 , Q = 0.01 and R d = 0.6 .
Vertical Plate
Cu -WaterNanofluid Ag -WaterNanofluid
ξ f ( ξ , 0 ) θ ( ξ , 0 ) θ ( ξ , 0 ) f ( ξ , 0 ) θ ( ξ , 0 ) θ ( ξ , 0 )
0.13.15021970.88862842.21025383.49012090.88597062.2335515
0.22.97831420.78369472.03530933.29547050.77911752.0526246
0.32.80552750.68745671.86610873.10064930.68166611.8784674
0.42.63541580.60138141.70607582.90965750.59498331.7145224
0.52.47110010.52600861.55768122.72592920.51946601.5631829
0.62.31490740.46102861.42223832.55194930.45465411.4256179
0.72.16818870.40552141.29994882.38908150.39949901.3018590
0.82.03129070.35823741.19010802.23757140.35265781.1910409
0.91.90347390.31779381.09125462.09648770.31269131.0915734
11.78085200.28241921.00000001.96151490.27781501.0000000
Horizontal plate
0.12.09299520.87682852.27309372.17091630.87037432.3299696
0.21.91475770.76031632.09147841.96511280.74902662.1346661
0.31.74314560.65372721.91299131.76961980.63931351.9443078
0.41.58576710.55950961.74224461.59368130.54357801.7638210
0.51.44926740.47886341.58312521.44505720.46273001.5971702
0.61.33814260.41163241.43819611.32864580.39621841.4467514
0.71.25421620.35657451.30845631.24598790.34241281.3132268
0.81.19708740.31186101.19349241.19600980.29918301.1957896
0.91.16563470.27559511.09179011.17703210.26446151.0925707
11.16742800.24643221.00000001.20052700.23706781.0000000

Share and Cite

MDPI and ACS Style

Mkhatshwa, M.; Motsa, S.; Sibanda, P. Overlapping Multi-Domain Spectral Method for Conjugate Problems of Conduction and MHD Free Convection Flow of Nanofluids over Flat Plates. Math. Comput. Appl. 2019, 24, 75. https://doi.org/10.3390/mca24030075

AMA Style

Mkhatshwa M, Motsa S, Sibanda P. Overlapping Multi-Domain Spectral Method for Conjugate Problems of Conduction and MHD Free Convection Flow of Nanofluids over Flat Plates. Mathematical and Computational Applications. 2019; 24(3):75. https://doi.org/10.3390/mca24030075

Chicago/Turabian Style

Mkhatshwa, Musawenkhosi, Sandile Motsa, and Precious Sibanda. 2019. "Overlapping Multi-Domain Spectral Method for Conjugate Problems of Conduction and MHD Free Convection Flow of Nanofluids over Flat Plates" Mathematical and Computational Applications 24, no. 3: 75. https://doi.org/10.3390/mca24030075

APA Style

Mkhatshwa, M., Motsa, S., & Sibanda, P. (2019). Overlapping Multi-Domain Spectral Method for Conjugate Problems of Conduction and MHD Free Convection Flow of Nanofluids over Flat Plates. Mathematical and Computational Applications, 24(3), 75. https://doi.org/10.3390/mca24030075

Article Metrics

Back to TopTop