Next Article in Journal
Finding Debt Cycles: QUBO Formulations for the Maximum Weighted Cycle Problem Solved Using Quantum Annealing
Next Article in Special Issue
Multi-Step-Ahead Wind Speed Forecast Method Based on Outlier Correction, Optimized Decomposition, and DLinear Model
Previous Article in Journal
Parameter Tuning for WECC Generic Photovoltaic System Models Using Latin Hypercube Sampling and Pareto Optimality
Previous Article in Special Issue
Studies on Finite Element Analysis in Hydroforming of Nimonic 90 Sheet
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Magnetohydrodynamics Williamson Nanofluid Flow over an Exponentially Stretching Surface with a Chemical Reaction and Thermal Radiation

1
Department of Mathematics and Computational Sciences, University of Zimbabwe, Mt. Pleasant, Harare P.O. Box MP167, Zimbabwe
2
Department of Mathematical and Computational Sciences, University of Venda, P. Bag X5050, Thohoyandou 0950, South Africa
*
Author to whom correspondence should be addressed.
Mathematics 2023, 11(12), 2740; https://doi.org/10.3390/math11122740
Submission received: 29 April 2023 / Revised: 5 June 2023 / Accepted: 13 June 2023 / Published: 16 June 2023

Abstract

:
Presented in this current study is the numerical analysis of magnetohydrodynamics Williamson nanofluid flow over an exponentially stretching surface. The most important aspect of the investigation is that the effects of the magnetic field, chemical reaction and thermal radiation in the fluid flow are taken into account. The partial differential equations governing the present Williamson nanofluid flow problem were observed to be highly nonlinear and coupled. Suitable similarity transformations were used to transmute the coupled system of nonlinear partial differential equations governing the fluid flow into a linear system. The linear system was solved numerically using the spectral quasi-linearization method. The MATLAB bvp4c numerical technique and a comparison with existing results for the skin friction coefficient were used to confirm the appropriateness of the method in solving the current problem. The influence of some pertinent physical parameters on the fluid’s velocity, temperature and concentration profiles were displayed graphically. The effects of all the physical parameters on the skin friction coefficient, Nusselt number and Sherwood number were portrayed in a tabular form. It was noted that enhancing the thermal radiation parameter reduces the fluid’s temperature, Nusselt number and the skin friction coefficient, while the Sherwood number is improved.

1. Introduction

Non-Newtonian fluids occur most often in industrial and engineering applications. The rheological properties of the non-Newtonian fluids cannot be explained using the famous Navier–Stokes equations. As a consequence, a number of models have been used to describe the characteristics of non-Newtonian fluids. These models include the Ellis model [1], Carreaus model [2], power law model [3], Cross model [4] and Casson model [5], to mention but a few. One special type of non-Newtonian model is the Williamson model [6], which was proposed to describe the flow of pseudoplastic materials. The boundary layer flow of the pseudoplastic materials has found applications in bio-engineering, chemical and nuclear industry, material processing and geophysics.
In fluid dynamics, Sakiadis [7] was the first researcher to study the boundary layer flow over a continuous stretching surface. An inaguaral study of fluid flow of Blasius type past a stretching surface was initiated by Crane [8]. The study of fluid flow over a stretching sheet has been a subject of interest in recent years due to its significant importance in areas such as metallurgical processes, polymer extrusion, plastic films, metal spinning, etc. There are quite a number of studies that have been performed on fluid flow past a stretching sheet, i.e., [9,10,11,12,13], among others.
In thermal engineering, the enhancement of the thermal characteristics of heat transfer fluids is a priority. The thermal conductivity and heat transfer qualities of the base fluid can be improved by dispersing nanosized (1–100 nm) solid particles into the fluid. These nanoparticles are usually metals, carbon nanotubes, oxides or carbines. The enhancement of heat transfer in fluids as a result of dispersing ultra-fine particles was first reported by Masuda et al. [14]. The term ‘nanofluid’, a fluid that contains dispersed nanoparticles, was introduced by Choi and Eastman [15]. A significant number of studies have been carried out on nanofluids, which inlcude the works by Elboughdiri et al. [16], Ashraf et al. [17], Nabwey et al. [18], Selimefendigil et al. [19] and Lou et al. [20].
Alfven [21] was the first to study the magnetic properties and the characteristics of fluids that are electrical conductors. Typical examples of such magnetofluids include electrolytes, plasmas, salt water and liquid metals. There has been growing interest in studying the MHD Williamson nanofluid. Abbas et al. [22] investigated the effects of heat generation and viscous dissipation on an MHD Williamson nanofluid that flows past a linear stretching sheet in a porous medium. The characteristics of MHD flow and heat transfer of a Williamson nanofluid flowing past a stretching sheet were examined by Reddy et al. [23]. Shawky et al. [24] used the Runge–Kutta method to analyze the heat and mass transfer of magnetohydrodynamic Williamson nanofluid flowing over a stretching sheet. The influence of Joule heating, heat generation/absorption, thermal radiation and chemical reaction on the MHD Williamson nanofluid flow over a stretching sheet through a porous medium was investigated by Bouslimi et al. [25]. Other notable works on the MHD Williamson nanofluid are [26,27,28,29,30,31,32,33], among others.
This current study mainly focuses on the Williamson nanofluid flow past an exponentially stretching surface with a chemical reaction and thermal radiation. This study has many applications in engineering and industrial processes. The Williamson fluid model with a chemical reaction has applications in water and air pollution, atmospheric flows and in chemical engineering problems such as food processes. Thermal radiation has applications in processes such as drying and distribution of temperature and moisture over agricultural fields  [34]. Nadeen and Hussain [35] used the homotopy analysis method to explore heat transfer effects on Williamson nanofluid flow over a porous exponentially stretching sheet. The Runge–Kutta–Fehlberg method was used to study the MHD flow of a Williamson nanofluid flow over an exponentially stretching surface by Kumar et al. [34]. Two cases of heat transfer, PEST and PEHF, were investigated on an MHD Williamson nanofluid flow over an exponentially stretching surface by Ahmed and Akbar [36]. Temperature-dependent viscosity and thermal conductivity in a Williamson nanofluid flow over an exponentially stretching sheet were studied by Amjad et al. [37]. Li et al. [38] used MATLAB’s bvp4c package to analyze heat and mass transfer in MHD Williamson nanofluid flow over an exponentially porous stretching surface.
Based on the aforementioned studies, it can be noted that there are many studies that have been performed on the Williamson nanofluid flow past an exponentially stretching surface. The novelty of this current study is the addition of thermophoresis and Brownian motion effects in the momentum equation. Additionally, the effects of the magnetic field, thermal radiation, chemical reaction, heat source and injection/suction parameters are simultaneously investigated in this model. The highly non-linear partial differential equations that govern the Williamson nanofluid flow are reduced into non-linear ordinary differential equations using suitable similarity transformations and then solved using the spectral quasi-linearization method (SQLM), developed by Motsa et al. [39]. The effects of some chosen pertinent parameters on the fluid velocity, temperature, concentration, skin friction coefficient, heat transfer rate and mass transfer rate were displayed using graphs and tables. The numerical results obtained in this current research work were validated by comparing the present results with those from MATLAB’s bvp4c routine and those results from already-published work. A very good agreement was established.

2. Fluid Model

The Williamson fluid model is used to describe the rheological behaviour of pseudoplastic materials over a wide range of shear stresses and shear rates. The continuity and momentum equations of an incompressible Williamson model are given, respectively, by [40]:
d i v V = 0 ,
ρ f d V d t = d i v S + ρ f b ,
where  d d t  is the time derivative and  b  is the specific body force vector. The Cauchy stress tensor  S = p I + τ * , [41], where p is the pressure term and I the identity vector. The extra stress tensor is given by:
τ * = μ + μ 0 μ 1 Γ γ ˙ A 1 ,
where the respective viscosities at zero and infinity shear rate are  μ 0  and  μ , respectively. The terms  A 1  and  Γ  are the first Revlin–Ericksen tensor and time constant, respectively. Additionally:
γ ˙ = π 2 , π = trace ( A 1 2 ) .
Choosing  μ = 0  and  Γ γ ˙ < 1  and applying the Binomial expansion, we have the extra stress tensor  τ * = μ 0 [ 1 + Γ γ ˙ ] A 1 .

3. Mathematical Analysis

Investigated in this current study is a two-dimensional flow of a steady incompressible Williamson nanofluid over a sheet that stretches exponentially. In this flow problem, the coordinate system is chosen in such a way that the x axis is along the stretching sheet and the y axis is measured normal to the sheet. At  y = 0 , the sheet is assumed to be stretching with a velocity  U w = U 0 e x l . The variable magnetic field  B ( x ) = B 0 e x 2 l  ( B 0  is a constant magnetic field) is applied perpendicular to the direction of flow. At the sheet, the fluid has temperature  T w = T + T 0 e x 2 l  and nanoparticle fraction  C w = C + C 0 e x 2 l . The ambient values of temperature and nanoparticle fraction, far way from the sheet, are denoted by  T  and  C , respectively. Figure 1 displays the schematic flow diagram and the coordinate system of the problem. Assuming that there is no pressure gradient and applying boundary layer approximations, the continuity, momentum and energy equations governing the flow are given by [25,40]:
u x + v y = 0 ,
u u x + v u y = ν y u y + Γ 2 u y 2 + g β T ( T T ) + g β C ( C C ) σ B 2 ρ f u , u T x + v T y = α 2 T y 2 + Q ( ρ c p ) f ( T T ) + τ D B T y C y + D T T T y 2
1 ( ρ c p ) f q r y ,
u C x + v C y = D B 2 C y 2 + D T T 2 T y 2 K ( C C ) ,
where u and v are the fluid velocity components in the x and y directions, respectively,  ν  is the kinematic viscosity of the fluid, g is the acceleration due to gravity,  β T  is the thermal expansion coefficient,  β C  is the concentration expansion coefficient,  σ  is the electrical conductivity,  α  is thermal diffusivity,  Q ( x ) ( = Q 0 e x l )  is the variable heat source,  ρ f  is the fluid density,  D B  is the Brownian diffusion coefficient,  D T  is the thermophoresis coefficient,  K ( x ) ( = K 0 e x l )  is the chemical reaction parameter and  τ = ( ρ c p ) p ( ρ c p ) f  is the ratio of the effective heat capacity of the nanoparticle material and heat capacity of the fluid.
The energy Equation (5) can be simplified by using the Rosseland approximation [42], which states that the radiative heat flux:
q r = 4 σ * 3 k * T 4 y ,
where  σ *  is the Stefan–Boltzmann constant and  k *  is the mean absorption coefficient. Assuming that the temperature differences within the flow are so small, the linear Taylor series expansion of  T 4  about  T  gives  T 4 4 T T 3 3 T 3  such that:
q r y = 16 σ * T 3 3 k * 2 T y 2 .
Using Equation (7) in Equation (5) gives:
u T x + v T y = α + 16 σ * T 3 3 ( ρ c p ) f k * 2 T y 2 + Q ( ρ c p ) f ( T T ) + τ D B T y C y + D T T T y 2 ,
The suitable boundary conditions for the system of Equations (3)–(6) are:
u = U w = U 0 e x l , v = γ ( x ) , where γ ( x ) = V 0 e x 2 l , T = T w = T + T 0 e x 2 l , C = C w = C + C 0 e x 2 l , at y = 0 , u 0 , T T , C C , as y .
The similarity transformations which are used to solve the governing equations are defined as follows [43]:
u = U 0 e x l f ( η ) , v = ν U 0 2 l e x 2 l [ f ( η ) + η f ( η ) ] , η = U 0 2 ν l y e x 2 l , T = T + T 0 e x 2 l θ ( η ) , C = C + C 0 e x 2 l ϕ ( η )
Using similarity transformations Equation (10), the continuity Equation (3) is identically satisfied and Equations (4)–(6) take the following form:
f + f f 2 f 2 + λ f f M 2 f + 2 G T θ + 2 G C ϕ = 0 ,
1 + 4 3 R d θ + P r ( f θ f θ + N b ϕ θ + N t θ 2 + δ θ ) = 0 ,
ϕ + S c ( f ϕ f ϕ K r ϕ ) + N t N b θ = 0 ,
subject to boundary conditions:
f ( 0 ) = S , f ( 0 ) = 1 , θ ( 0 ) = 1 , ϕ ( 0 ) = 1 , f ( ) 0 , θ ( ) 0 , ϕ ( ) 0 .
where  λ ( = Γ U 0 3 e 3 x l ν l )  is the Williamson fluid parameter,  M 2 ( = 2 l σ B 0 2 ρ U 0 )  is the magnetic field parameter,  G T ( = g l B T T 0 U 0 2 )  is the thermal Grashof number,  G C ( = g l B C C 0 U 0 2 )  is the mass Grashof number,  P r ( = ν α )  is the Prandlt number,  R d ( = 4 σ * T 3 k * κ )  is the radiation parameter,  N b ( = τ D B ν ( C w C ) )  is the Brownian motion parameter,  N t ( = τ D T ν T ( T w T ) )  is the thermophoresis parameter,  δ ( = 2 l Q 0 ρ c p U 0 )  is the heat generation parameter,  S c ( = ν D B )  is the Schmidt number,  S ( = V 0 2 l ν U 0 )  is the suction  ( S < 0 )  or the injection  ( S > 0 )  parameter and  K r ( = 2 l K 0 U 0 )  is the chemical reaction parameter.
The skin friction coefficient  ( c f ) , the local Nusselt number  ( N u x )  and the local Sherwood number  ( S h x )  are the physical quantities of engineering significance discussed in this study. Following the work by Ahmed and Akbar [36]:
c f = 1 ρ U w 2 μ u y + Γ 2 u y 2 y = 0 , N u x = 2 l ( T w T ) e x 2 l T y y = 0 , S h x = 2 l ( C w C ) e x 2 l C y y = 0 .
Using similarity transformations in Equation (10), the following dimensionless forms are obtained:
2 R e x c f = f ( 0 ) + λ 2 ( f ( 0 ) ) 2 , N u x R e x = θ ( 0 ) , S h x R e x = ϕ ( 0 ) ,
where  R e x = U w l ν  is the Reynolds number.

4. Method of Solution

In this study, the spectral quasi-linearization method is used to seek the numerical solution of the coupled system of Equations (11)–(13) subject to boundary conditions Equation (18). A Newton–Raphson-based quasi-linearization method [44], which uses first-order Taylor series expansion, is used to linearize the non-linear terms. Denote the respective solutions of Equations (11)–(13) at iteration level s by  f s θ s  and  ϕ s , respectively. Assuming that the difference between solutions at iteration level s and  s + 1  are sufficiently close, quasi-linearization gives the following iterative sequence of linear differential equations:
a 0 , s f s + 1 + a 1 , s f s + 1 + a 2 , s f s + 1 + a 3 , s f s + 1 + a 4 , s θ s + 1 + a 5 , s ϕ s + 1 = R 1 , s ,
b 0 , s θ s + 1 + b 1 , s θ s + 1 + b 2 , s θ s + 1 + b 3 , s f s + 1 + b 4 , s f s + 1 + b 5 , s ϕ s + 1 = R 2 , s ,
c 0 , s ϕ s + 1 + c 1 , s ϕ s + 1 + c 2 , s ϕ s + 1 + c 3 , s f s + 1 + c 4 , s f s + 1 + c 5 , s θ s + 1 = R 3 , s ,
where the variable coefficients known at iteration level s are defined as:
a 0 , s = 1 + λ f s , a 1 , s = f s + λ f s , a 2 , s = 4 f s M , a 3 , s = f s , a 4 , s = 2 G T , a 5 , s = 2 G C , b 0 , s = 1 + 4 3 R d , b 1 , s = P r ( f s + N b ϕ s + 2 N t θ s ) , b 2 , s = P r ( f s δ ) , b 3 , s = P r θ s , b 4 , s = P r θ s , b 5 , s = P r N b θ s , c 0 , s = 1 , c 1 , s = S c f s , c 2 , s = S c ( f s + K r ) , c 3 , s = S c ϕ s , c 4 , s = S c ϕ s , c 5 , s = N t N b .
The boundary conditions given in Equation (18) are transformed to:
f s + 1 ( 0 ) = 1 , f s + 1 ( 0 ) = S , θ s + 1 ( 0 ) = 1 , ϕ s + 1 ( 0 ) = 1 , f s + 1 ( ) 0 , θ s + 1 ( ) 0 , θ s + 1 ( ) 0 .
The terms on the right hand side are:
R 1 , s = f s f s 2 f s 2 + λ f s f s , R 2 , s = P r ( f s θ s f s θ + N b θ s ϕ s + N t θ s 2 ) , R 3 , s = S c ( f s ϕ s f s ϕ s )
The unknown functions  f s + 1 θ s + 1  and  ϕ s + 1  are approximated using Chebyshev interpolating polynomials, such that the their derivatives evaluated at Gauss–Lobatto collocation points  η i = cos π i N ( i = 0 , 1 , 2 , , N )  are given by:
d n f s + 1 d η ( η i ) = k = 0 N D i k n f s + 1 ( η k ) = D n F s + 1 , d n θ s + 1 d η ( η i ) = k = 0 N D i k n θ s + 1 ( η k ) = D n Θ s + 1 , d n ϕ s + 1 d η ( η i ) = k = 0 N D i k n ϕ s + 1 ( η k ) = D n Φ s + 1 ,
where
  • D = 2 L D F s + 1 = [ f s + 1 ( η 0 ) , f s + 1 ( η 1 ) , , f s + 1 ( η N 1 ) , f s + 1 ( η N ) ] T ,
  • Θ s + 1 = [ θ s + 1 ( η 0 ) , θ s + 1 ( η 1 ) , , θ s + 1 ( η N 1 ) , θ s + 1 ( η N ) ] T ,
  • Φ s + 1 = [ ϕ s + 1 ( η 0 ) , ϕ s + 1 ( η 1 ) , , ϕ s + 1 ( η N 1 ) , ϕ s + 1 ( η N ) ] T .
Using derivatives in Equation (19) in the system of Equations (15)–(17) yields a system in vector matrix form:
A 11 A 12 A 13 A 21 A 22 A 23 A 31 A 32 A 33 F s + 1 Θ s + 1 Φ s + 1 = R 1 , s R 2 , s R 3 , s ,
where
A 11 = a 0 , s D 3 + a 1 , s D 2 + a 2 , s D + a 3 , s , A 12 = a 4 , s I , A 13 = a 5 , s I , A 21 = b 3 , s D + b 4 , s , A 22 = b 0 , s D 2 + b 1 , s D + b 2 , s , A 23 = b 5 , s D , A 31 = c 3 , s D + c 4 , s , A 32 = c 5 , s D 2 , A 33 = c 0 , s D 2 + c 1 , s D + c 2 , s .
where  I  is an  ( N + 1 ) × ( N + 1 )  identity matrix. The spectral boundary conditions are:
f s + 1 ( η N ) = S , k = 0 N D N , k f s + 1 ( η N ) = 1 , θ s + 1 ( η N ) = 1 , ϕ s + 1 ( η N ) = 1 , k = 0 N D 0 , k f s + 1 ( η 0 ) = 0 , θ s + 1 ( η 0 ) = 0 , ϕ s + 1 ( η 0 ) = 0 .
The numerical iteration of the SQLM, coded in MATLAB R2022b on an Intel(R) Core(TM) i5, is started by using the initial guesses that satisfy the boundary conditions Equation (18), given by:
f 0 ( η ) = 1 e η + S , θ 0 ( η ) = e η , ϕ 0 ( η ) = e η .

5. Results and Discussion

5.1. Validation of Results

To confirm the accuracy of the SQLM used in this study, the values of the skin friction  ( f ( 0 ) + λ 2 ( f ( 0 ) ) 2 )  are compared against the MATLAB bvp4c solver results and the homotopy analysis method results obtained by Nadeem and Hussain [40] and Amjad et al. [45]. Considering the values  G T = G C = 0 , Equation (11) reduces to the problem by Amjad et al. [45], which is given by:
f + f f 2 f 2 + λ f f M 2 f = 0 ,
subject to boundary conditions:
f ( 0 ) = S , f ( 0 ) = 1 , f ( ) 0 .
Using  N = 40  collocation points,  η = 5.0 , the MATLAB SQLM algorithm for solving Equation (20) involves iteratively solving the following recursive sequence:
D 0 , 0 D 0 , 1 D 0 , N 1 , D 0 , N A D N , 0 D N , 1 D N , N 1 D N , N 0 0 1 f s + 1 ( η 0 ) F s + 1 f s + 1 ( η N 1 ) f s + 1 ( η N ) = 0 R s 1 S ,
where  A = a 0 , s D 3 + a 1 , s D 2 + a 2 , s D + a 3 , s F s + 1 = [ f s + 1 ( η 1 ) , f s + 1 ( η 2 ) , , f s + 1 ( η N 3 ) , f s + 1 ( η N 2 ) ] T  and  R s = 2 F s F s + F s F s ) + λ F s F s . Performing 20 iterations, the results obtained for  ( f ( 0 ) + λ 2 ( f ( 0 ) ) 2 )  are displayed in Table 1.
Table 1 displays the computed values of the skin friction coefficient compared against the results by Amjad et al. [45] for different values of  λ , S and M. A good match of the results is observed. The accuracy of the SQLM was validated by a direct comparison with the reported results.
Considering  G T = G C = 0 = M = 0 , Equation (11) reduces to the problem by Nadeem and Hussain [40], which can be solved using MATLAB’s bvp4c solver by first using the following substitutions:
y ( 1 ) = f , y ( 2 ) = f , y ( 3 ) = f
f = y ( 3 ) = 2 ( y ( 2 ) ) 2 y ( 1 ) y ( 3 ) ( 1 + λ y ( 3 ) )
and the boundary conditions are given as  y a ( 1 ) + S , y a ( 2 ) 1 , y b ( 2 ) . The three first-order equations are coded in the MATLAB’s bvp4c solver with the function name “odeBVP”, and with “odeBc” handling the boundary conditions. Choosing the interval of integration to  [ 0 , 40 ] , the solutions from the function “bvp4c” are given by:
  • sol = bvp4c(@odeBVP, @odeBc, solinit, options).
A comparison of the SQLM skin friction coefficient values against those from MATLAB’s bvp4c routine and the results by Nadeem and Hussain [40] is shown in Table 2. A perfect agreement was observed.

5.2. Results

The SQLM algorithm was implemented using MATLAB R2022b software. For all the numerical results, unless otherwise stated, the default parameters considered are:  N = 60 M = 0.1 G T = 0.1 λ = 0.3 G C = 0.1 N t = 0.1 N b = 0.1 P r = 0.5 δ = 0.2 S = 0.0 K r = 0.1 S c = 1.0 . The convergence and accuracy of the spectral quasi-linearization method were verified using the solution error norms and residual errors, respectively. The solution error norms, defined as the difference between values of successive iterations are denoted by [46]:
E r r [ F ( η ) ] = | | F s + 1 ( η ) F s ( η ) | | , E r r [ Θ ( η ) ] = | | Θ s + 1 ( η ) Θ s ( η ) | | , E r r [ Φ ( η ) ] = | | Φ s + 1 ( η ) Φ s ( η ) | | .
The residual error gives a measure of the extent to which the SQLM solution approximates the true solution. The residual  L  norms are given by Alharbey et al. [47] as:
R e s ( f ) = | | f + f f 2 f 2 + λ f f M 2 f + 2 G T θ + 2 G C ϕ | | , R e s ( θ ) = | | ( 1 + ( 4 / 3 ) R d θ + P r ( f θ f θ + N b ϕ θ + N t θ 2 + δ θ ) | | , R e s ( ϕ ) = | | ϕ + S c ( f ϕ f ϕ K r ϕ ) + ( N t / N b ) θ | | ,
Figure 2a reveals that the SQLM converges after only five iterations with a solution based error of order ≈ 10 10 . Additionally, after three iterations, the SQLM achieves an accuracy of order ≈ 10 9 , as shown in Figure 2b.
Figure 3, Figure 4 and Figure 5 display the effects of the magnetic parameter  ( M ) , suction/injection parameter  ( S )  and the Williamson parameter  ( λ ) , respectively, on the fluid velocity profiles. Figure 3 shows that the fluid velocity is depressed as the magnetic parameter is increased. Physically, the fluid velocity drops due to the resistive Lorentz force, which is induced by the magnetic parameter. It is depicted in Figure 4 that the Williamson nanofluid velocity profiles are depressed when the suction parameter is increased. Additionally, it is revealed in Figure 5 that there is an inverse relationship between the fluid velocity distribution and the non-Newtonian Williamson parameter. An increase in the values of  λ  causes a decrease in the fluid movement and reduces the boundary layer thickness.
The influences of the Prandtl number (Pr), Brownian motion parameter ( N b ), thermal Grashof parameter ( G T ) and thermal radiation parameter ( R d ) on the Williamson nanofluid dimensionless temperature ( θ ) are depicted in Figure 6, Figure 7, Figure 8 and Figure 9, respectively. It is displayed in Figure 6 that the fluid temperature and thermal boundary layer are reduced as the Prandtl number increases. The Prandtl number can be viewed as the ratio of momentum to thermal boundary layers. Physically, a high Prandtl number means a small thermal boundary layer. It is revealed in Figure 7 that an increase in the values of the Brownian motion parameter increases the fluid temperature profile. An increase in the Brownian motion parameter results in an increased kinetic energy of the Williamson nanoparticles, and hence, a temperature increase. Figure 8 depicts that the fluid temperature is depressed as the thermal Grashof number is enhanced. Essentially, the Grashof number signifies the ratio of buoyancy to viscous forces. Increasing  G T  results in an addition of more thermal energy in the fluid molecules, which in turn increases the fluid local heat transfer rate. The thermal boundary layer is reduced, and hence, the temperature profiles decrease. It is shown in Figure 9 that temperature is an increasing function of the thermal radiation parameter. The effect of increasing the thermal radiation parameter is thickening the thermal boundary layer, and hence, the temperature profiles are increased.
Figure 10 portrays the influence of the Schmidt number (Sc) on the Williamson nanoparticle concentration.  S c  can be defined as the ratio of momentum diffusivity and mass diffusivity. High values of  S c  corresponds to a weaker solute diffusivity and the concentration distribution and solute boundary layer decrease as a consequence. The fluid dimensionless concentration profiles are depressed when the chemical reaction parameter is increased as seen in Figure 11. Physically, when the chemical reaction parameter is increased, quite a number of solute molecules will undergo chemical reaction, and hence, the reduction in the concentration. The influence of the mass Grashof number  ( G C )  on the concentration is displayed in Figure 12 G C  relates species buoyancy force to the viscous hydrodynamic force. Increasing  G C  causes an enhancement of the concentration gradient, which in turn boosts the buoyancy effect. A resulting induced flow will cause a decrease in concentration, and hence, a decrease in the concentration profile, as depicted in Figure 12.
Table 3 displays the effects of all the pertinent thermo-physical parameters involved in the current problem on the skin friction, Nusselt number and Sherwood number. The skin friction coefficient upsurges as the values of injection parameter, Prandtl number, magnetic parameter, Brownian motion parameter and Schmidt number are increased. The fluid flow is improved by the functioning magnetic field regarding the Williamson nanofluid, and thus, increases the surface friction. The increase of  2 R e C f  with an increasing Prandtl number is attributed to increased fluid momentum. The opposite trend is observed when the Williamson fluid parameter, mass Grashof number, thermal Grashof number, thermal radiation parameter, thermophoresis parameter, heat generation parameter and chemical reaction parameter are increased.
It is also noted that the mass Grashof number and the thermal Grashof number are the only parameters whose increment enhances the Nusselt number. Increasing the values of the Williamson fluid parameter, injection parameter, magnetic parameter, Prandtl number, thermal radiation, Brownian motion parameter, thermophoresis parameter, heat generation parameter, Schmidt number and chemical reaction parameter suppresses the heat transfer rate. It is noted that Nusselt number is a decreasing function of  N b  and  N t . Physically, the effects of both Brownian motion and themorphoresis effects move the Williamson nanoparticles away from the stretching sheet sheet, intensifying the the diffusion of the nanoparticles into the boundary layer, and hence, causing a decrease in the Nusselt number.
The Sherwood number is improved as the values of mass Grashof number, thermal Grashof number, thermal radiation parameter, Brownian motion parameter, heat generation parameter, Schmidt number and chemical reaction parameter are increased and depreciates as the Williamson fluid parameter, injection parameter, magnetic parameter, Prandtl number and thermophoresis parameter are increased. The Schmidt number is the relative effectiveness of the momentum and mass transport by diffusion in the hydrodynamic and species boundary layers. An increase in  S c  will result in an increase in the Sherwood number. An increased chemical reaction parameter means there will be more interaction of species concentration with the momentum boundary, and hence, an increase in the Sherwood number.

6. Conclusions

In this manuscript, the spectral quasi-linearization method was applied to numerically analyze the magnetoxhydrodynamics Williamson nanofluid flow over an exponentially stretching surface with chemical reaction and thermal radiation nanofluid flow. A comparison of the skin friction coefficient results obtained from MATLAB’s bvp4c solver and published work confirmed that the method is reliable for solving the current problem. The key findings from the study are as follows:
  • The dimensionless velocity  ( f ( η ) )  diminishes as the values of the magnetic parameter are increased from 0 to  1.2 ;
  • The dimensionless temperature  ( θ ( η ) )  is an increasing function of  0.1 N b 1.6  and  0.0 R d 0.9 ;
  • The dimensionless concentration  ( ϕ ( η ) )  decreases for  0.5 S c 0.8  and  0.1 K r 0.7 ;
  • The skin friction coefficient increases as  M ( 0.1 M 0.3 )  and  N b ( 0.1 N b 0.5 )  increase and depressed for increased values of  N t ( 0.1 N t 0.5 ) ;
  • The Nusselt number diminishes as  M ( 0.1 M 0.3 ) N b ( 0.1 N b 0.5 )  and  N t ( 0.1 N t 0.5 )  are increased;
  • The Sherwood number decreases as  M ( 0.1 M 0.3 )  and  N t ( 0.1 N t 0.37 )  increase and decreases as  N b ( 0.3 N b 0.7 )  increases.

Author Contributions

Conceptualization, H.M. and S.S.; Methodology, H.M.; Validation, S.S.; Formal analysis, H.M.; Investigation, S.S.; Writing—original draft, H.M. and S.S.; Writing—review & editing, H.M. and S.S.; Supervision, S.S.; Funding acquisition, S.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

Data is contained within the article.

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

  x , y Cartesian coordinates [m]
  u , v Velocity components in the x and y directions, respectively [m s 1 ]
  U 0 Reference velocity [m s 1 ]
  β T Thermal expansion coefficient
  β C Concentration expansion coefficient
  B 0 Magnetic field strength [NmA 1 ]
  C f Skin friction coefficient
  P r Prandtl number
MMagnetic parameter [Te]
TFluid temperature [T]
  C w Concentration of nanoparticles at the surface [mol m 3 ]
CConcentration of nanoparticles [mol m 3 ]
  U w Velocity at the wall [m s 1 ]
QHeat source
  K r Chemical reaction parameter [Ms 1 ]
  T 0 Reference temperature [K]
  C 0 Reference concentration [mol m 3 ]
  R d Thermal radiation parameter
  q r Radiative heat flux [J]
SSuction/injection parameter
  θ ( η ) Dimensionless temperature
  ϕ ( η ) Dimensionless concentration
  N t Thermophoretic parameter
  N u x Local Nusselt number
  S h x Local Sherwood number
  T w Surface temperature [K]
  T Ambient temperature [K]
fDimensionless stream function
gAcceleration due to gravity [m s 2 ]
  S c Schmidt number
  D B Brownian diffusion coefficient [m 2  s 1 ]
  D T Thermophoresis diffusion coefficient [m 2  s 1 ]
  μ Infinite viscosity [Nsm 2 ]
  ( ρ c p ) f Heat capacity of the nanofluid [Jm 3 K 1 ]
  R e x Reynolds number
  f ( η ) Velocity profile
  η Dimensionless similarity variable
  σ Electrical conductivity [Sm 1 ]
  Γ Positive time constant
  α Thermal diffusivity [m 2  s 1 ]
  ( ρ c p ) p Heat capacity of the nanoparticles [Jm 3  K 1 ]
  ν Kinematic viscosity [m 2  s 1 ]
  ρ f Fluid density [kg m 3 ]
  λ Williamson fluid parameter
  G r T Thermal Grashof number
  G r C Concentration Grashof number

References

  1. Matsuhisa, S.; Bird, R.B. Analytical and numerical solutions for laminar flow of the non-Newtonian ellis fluid. AIChE J. 1965, 11, 588–595. [Google Scholar] [CrossRef]
  2. Carreau, P.J. Rheological Equations from Molecular Network Theories. Trans. Soc. Rheol. 1972, 16, 99–127. [Google Scholar] [CrossRef]
  3. Ostwald, W. Ueber die rechnerische Darstellung des Strukturgebietes der Viskosität. Kolloid Z. 1929, 47, 176–187. [Google Scholar] [CrossRef]
  4. Cross, M.M. Rheology of non-Newtonian fluids: A new flow equation for pseudoplastic systems. J. Colloid Sci. 1965, 20, 417–437. [Google Scholar] [CrossRef]
  5. Casson, N. A Flow Equation for Pigment-Oil Suspensions of the Printing Ink Type. In Rheology of Disperse Systems; Mill, C.C., Ed.; Pergamon Press: Oxford, UK, 1959; pp. 84–104. [Google Scholar]
  6. Williamson, R.V. The Flow of Pseudoplastic Materials. Ind. Eng. Chem. 1929, 21, 1108–1111. [Google Scholar] [CrossRef]
  7. Sakiadis, B.C. Boundary-layer behavior on continuous solid surfaces: I. Boundary-layer equations for two-dimensional and axisymmetric flow. AIChE J. 1961, 7, 26–28. [Google Scholar] [CrossRef]
  8. Crane, L.J. Flow past a stretching plate. Z. Angew. Math. Phys. ZAMP 1970, 21, 645–647. [Google Scholar] [CrossRef]
  9. Khan, S.; Selim, M.M.; Khan, A.; Ullah, A.; Abdeljawad, T.; Ikramullah; Ayaz, M.; Mashwani, W.K. On the Analysis of the Non-Newtonian Fluid Flow Past a Stretching/Shrinking Permeable Surface with Heat and Mass Transfer. Coatings 2021, 11, 566. [Google Scholar] [CrossRef]
  10. Zeb, H.; Bhatti, S.; Khan, U.; Wahab, H.A.; Mohamed, A.; Khan, I. Impact of Homogeneous-Heterogeneous Reactions on Flow of Non-Newtonian Ferrofluid over a Stretching Sheet. J. Nanomater. 2022, 2022, 2501263. [Google Scholar] [CrossRef]
  11. Mahabaleshwar, U.S.; Maranna, T.; Sofos, F. Analytical investigation of an incompressible viscous laminar Casson fluid flow past a stretching/shrinking sheet. Sci. Rep. 2022, 12, 18404. [Google Scholar] [CrossRef]
  12. Abbas, N.; Nadeem, S.; Shatanawi, W. Effects of radiation and heat generation for non-Newtonian fluid flow over slendering stretching sheet: Numerically. J. Appl. Math. Mech. 2022, 103, e202100299. [Google Scholar] [CrossRef]
  13. Akbar, N.S.; Al-Zubaidi, A.; Saleem, S.; Alsallami, S.A.M. Variable fluid properties analysis for thermally laminated 3-dimensional magnetohydrodynamic non-Newtonian nanofluid over a stretching sheet. Sci. Rep. 2023, 13, 3231. [Google Scholar] [CrossRef]
  14. Masuda, H.; Ebata, A.; Teramae, K.; Hishinuma, N. Alteration of Thermal Conductivity and Viscosity of Liquid by Dispersing Ultra-Fine Particles. Dispersion of Al2O3, SiO2 and TiO2 Ultra-Fine Particles. Netsu Bussei 1993, 7, 227–233. [Google Scholar] [CrossRef] [Green Version]
  15. Choi, S.U.S.; Eastman, J.A. Enhancing thermal conductivity of fluids with nanoparticle. In Proceedings of the 1995 International Mechanical Engineering Congress and Exhibition, San Francisco, CA, USA, 12–17 November 1995. [Google Scholar]
  16. Elboughdiri, N.; Ghernaout, D.; Muhammad, T.; Alshehri, A.; Sadat, R.; Ali, M.R.; Wakif, A. Towards a novel EMHD dissipative stagnation point flow model for radiating copper-based ethylene glycol nanofluids: An unsteady two-dimensional homogeneous second-grade flow case study. Case Stud. Therm. Eng. 2023, 45, 102914. [Google Scholar] [CrossRef]
  17. Ashraf, M.Z.; Rehman, S.U.; Farid, S.; Hussein, A.K.; Ali, B.; Shah, N.A.; Weera, W. Insight into Significance of Bioconvection on MHD Tangent Hyperbolic Nanofluid Flow of Irregular Thickness across a Slender Elastic Surface. Mathematics 2022, 10, 2592. [Google Scholar] [CrossRef]
  18. Nabwey, H.A.; Rahbar, F.; Armaghani, T.; Rashad, A.M.; Chamkha, A.J. A Comprehensive Review of Non-Newtonian Nanofluid Heat Transfer. Symmetry 2023, 15, 362. [Google Scholar] [CrossRef]
  19. Selimefendigil, F.; Şenol, G.; Öztop, H.F.; Abu-Hamdeh, N.H. A Review on Non-Newtonian Nanofluid Applications for Convection in Cavities under Magnetic Field. Symmetry 2022, 15, 41. [Google Scholar] [CrossRef]
  20. Lou, Q.; Ali, B.; Rehman, S.U.; Habib, D.; Abdal, S.; Shah, N.A.; Chung, J.D. Micropolar Dusty Fluid: Coriolis Force Effects on Dynamics of MHD Rotating Fluid When Lorentz Force Is Significant. Mathematics 2022, 10, 2630. [Google Scholar] [CrossRef]
  21. Alfven, H. Existence of Electromagnetic-Hydrodynamic Waves. Nature 1942, 150, 405–406. [Google Scholar] [CrossRef]
  22. Abbas, A.; Jeelani, M.B.; Alnahdi, A.S.; Ilyas, A. MHD Williamson Nanofluid Fluid Flow and Heat Transfer Past a Non-Linear Stretching Sheet Implanted in a Porous Medium: Effects of Heat Generation and Viscous Dissipation. Processes 2022, 10, 1221. [Google Scholar] [CrossRef]
  23. Reddy C, S.; Naikoti, K.; Rashidi, M.M. MHD flow and heat transfer characteristics of Williamson nanofluid over a stretching sheet with variable thickness and variable thermal conductivity. Trans. A. Razmadze Math. Inst. 2017, 171, 195–211. [Google Scholar] [CrossRef]
  24. Shawky, H.M.; Eldabe, N.T.M.; Kamel, K.A.; Abd-Aziz, E.A. MHD flow with heat and mass transfer of Williamson nanofluid over stretching sheet through porous medium. Microsyst. Technol. 2019, 25, 1155–1169. [Google Scholar] [CrossRef]
  25. Bouslimi, J.; Omri, M.; Mohamed, R.A.; Mahmoud, K.H.; Abo-Dahab, S.M.; Soliman, M.S. MHD Williamson nanofluid flow over a stretching sheet through a porous medium under effects of joule heating, nonlinear thermal radiation, heat generation/absorption, and chemical reaction. Adv. Math. Phys. 2021, 2021, 9950993. [Google Scholar] [CrossRef]
  26. Hayat, T.; Bashir, G.; Waqas, M.; Alsaedi, A. MHD 2D flow of Williamson nanofluid over a nonlinear variable thicked surface with melting heat transfer. J. Mol. Liq. 2016, 223, 836–844. [Google Scholar] [CrossRef]
  27. Ibrahim, W.; Negera, M. The Investigation of MHD Williamson Nanofluid over Stretching Cylinder with the Effect of Activation Energy. Adv. Math. Phys. 2020, 2020, 9523630. [Google Scholar] [CrossRef]
  28. Eswara Rao, M.; Siva Sankari, M.; Nagalakshmi, C.; Rajkumar, S. On the Role of Bioconvection and Activation Energy for MHD-Stretched Flow of Williamson and Casson Nanofluid Transportation across a Porous Medium Past a Permeable Sheet. J. Nanomater. 2023, 2023, 3995808. [Google Scholar] [CrossRef]
  29. Asjad, M.I.; Zahid, M.; Ali, B.; Jarad, F. Unsteady MHD Williamson Fluid Flow with the Effect of Bioconvection over Permeable Stretching Sheet. Math. Probl. Eng. 2022, 2022, 7980267. [Google Scholar] [CrossRef]
  30. Wang, F.; Asjad, M.I.; Rehman, S.U.; Ali, B.; Hussain, S.; Gia, T.N.; Muhammad, T. MHD Williamson Nanofluid Flow over a Slender Elastic Sheet of Irregular Thickness in the Presence of Bioconvection. Nanomaterials 2021, 11, 2297. [Google Scholar] [CrossRef]
  31. Ahmed, K.; Akbar, T.; Muhammad, T.; Alghamdi, M. Heat transfer characteristics of MHD flow of Williamson nanofluid over an exponential permeable stretching curved surface with variable thermal conductivity. Case Stud. Therm. Eng. 2021, 28, 101544. [Google Scholar] [CrossRef]
  32. Patil, V.S.; Humane, P.P.; Patil, A.B. MHD Williamson nanofluid flow past a permeable stretching sheet with thermal radiation and chemical reaction. Int. J. Model. Simul. 2023, 43, 185–199. [Google Scholar] [CrossRef]
  33. Khan, M.; Malik, M.Y.; Salahuddin, T.; Hussian, A. Heat and mass transfer of Williamson nanofluid flow yield by an inclined Lorentz force over a nonlinear stretching sheet. Results Phys. 2018, 8, 862–868. [Google Scholar] [CrossRef]
  34. Kumar, P.B.S.; Gireesha, B.J.; Gorla, R.S.R.; Mahanthesh, B. Magnetohydrodynamic Flow of Williamson Nanofluid Due to an Exponentially Stretching Surface in the Presence of Thermal Radiation and Chemical Reaction. J. Nanofluids 2017, 6, 264–272. [Google Scholar] [CrossRef]
  35. Nadeem, S.; Hussain, S.T.; Lee, C. Flow of a Williamson fluid over a stretching sheet. Braz. J. Chem. Eng. 2013, 30, 619–625. [Google Scholar] [CrossRef] [Green Version]
  36. Ahmed, K.; Akbar, T. Numerical investigation of magnetohydrodynamics Williamson nanofluid flow over an exponentially stretching surface. Adv. Mech. Eng. 2021, 13, 168781402110198. [Google Scholar] [CrossRef]
  37. Amjad, M.; Ahmed, I.; Ahmed, K.; Alqarni, M.S.; Akbar, T.; Muhammad, T. Numerical Solution of Magnetized Williamson Nanofluid Flow over an Exponentially Stretching Permeable Surface with Temperature Dependent Viscosity and Thermal Conductivity. Nanomaterials 2022, 12, 3661. [Google Scholar] [CrossRef]
  38. Li, Y.-X.; Alshbool, M.H.; Lv, Y.-P.; Khan, I.; Riaz Khan, M.; Issakhov, A. Heat and mass transfer in MHD Williamson nanofluid flow over an exponentially porous stretching surface. Case Stud. Therm. Eng. 2021, 26, 100975. [Google Scholar] [CrossRef]
  39. Motsa, S.S.; Dlamini, P.G.; Khumalo, M. Spectral Relaxation Method and Spectral Quasilinearization Method for Solving Unsteady Boundary Layer Flow Problems. Adv. Math. Phys. 2014, 2014, 341964. [Google Scholar] [CrossRef]
  40. Nadeem, S.; Hussain, S.T. Heat transfer analysis of Williamson fluid over exponentially stretching surface. Appl. Math. Mech. 2014, 35, 489–502. [Google Scholar] [CrossRef]
  41. Dapra, I.; Scarpi, G. Perturbation solution for pulsatile flow of a non-Newtonian Williamson fluid in a rock fracture. Int. J. Rock Mech. Min. Sci. 2007, 44, 271–278. [Google Scholar] [CrossRef]
  42. Rosseland, S. Theoretical Astrophysics; Oxford University Press: Oxford, UK, 1936. [Google Scholar]
  43. Seini, Y.I.; Makinde, O.D. MHD Boundary Layer Flow due to Exponential Stretching Surface with Radiation and Chemical Reaction. Math. Probl. Eng. 2013, 2013, 163614. [Google Scholar] [CrossRef] [Green Version]
  44. Bellman, R.; Kalaba, R. Quasilinearization and Nonlinear Boundary-Value Problems; American Elsevier Publishing Company: New York, NY, USA, 1965. [Google Scholar]
  45. Amjad, M.; Ahmed, K.; Akbar, T.; Muhammad, T.; Ahmed, I.; Alshomrani, A.S. Numerical investigation of double diffusion heat flux model in Williamson nanofluid over an exponentially stretching surface with variable thermal conductivity. Case Stud. Therm. Eng. 2022, 36, 102231. [Google Scholar] [CrossRef]
  46. Motsa, S. On the New Bivariate Local Linearisation Method for Solving Coupled Partial Differential Equations in Some Applications of Unsteady Fluid Flows with Heat and Mass Transfer. In Mass Transfer—Advancement in Process Modelling; InTech: Vienna, Austria, 2015. [Google Scholar] [CrossRef] [Green Version]
  47. Alharbey, R.A.; Mondal, H.; Behl, R. Spectral Quasi-Linearization Method for Non-Darcy Porous Medium with Convective Boundary Condition. Entropy 2019, 21, 838. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Schematic flow diagram and coordinate system.
Figure 1. Schematic flow diagram and coordinate system.
Mathematics 11 02740 g001
Figure 2. Error graphs of  f ( η ) θ ( η )  and  ϕ ( η ) .
Figure 2. Error graphs of  f ( η ) θ ( η )  and  ϕ ( η ) .
Mathematics 11 02740 g002
Figure 3. Influence of M on the nanofluid velocity.
Figure 3. Influence of M on the nanofluid velocity.
Mathematics 11 02740 g003
Figure 4. Influence of S on the nanofluid velocity.
Figure 4. Influence of S on the nanofluid velocity.
Mathematics 11 02740 g004
Figure 5. Influence of  λ  on fluid velocity.
Figure 5. Influence of  λ  on fluid velocity.
Mathematics 11 02740 g005
Figure 6. Influence of  P r  on fluid temperature.
Figure 6. Influence of  P r  on fluid temperature.
Mathematics 11 02740 g006
Figure 7. Influence of  N b  on fluid temperature.
Figure 7. Influence of  N b  on fluid temperature.
Mathematics 11 02740 g007
Figure 8. Influence of  G T  on fluid temperature.
Figure 8. Influence of  G T  on fluid temperature.
Mathematics 11 02740 g008
Figure 9. Influence of  R d  on fluid temperature.
Figure 9. Influence of  R d  on fluid temperature.
Mathematics 11 02740 g009
Figure 10. Influence of  S c  on fluid concentration.
Figure 10. Influence of  S c  on fluid concentration.
Mathematics 11 02740 g010
Figure 11. Influence of  K r  on fluid concentration.
Figure 11. Influence of  K r  on fluid concentration.
Mathematics 11 02740 g011
Figure 12. Influence of  G c  on fluid concentration.
Figure 12. Influence of  G c  on fluid concentration.
Mathematics 11 02740 g012
Table 1. Table of present values of  2 R e C f  compared against published results for varying values of  λ , S and M.
Table 1. Table of present values of  2 R e C f  compared against published results for varying values of  λ , S and M.
  ( f ( 0 ) + λ 2 ( f ( 0 ) ) 2 )
  λ SMAmjad et al. [45]SQLM
0.10.22.01.7542131.754213105760364
0.20.22.01.6786751.678675073146794
0.30.22.01.5798271.578533717157394
0.10.12.01.7992491.799249869955796
0.10.22.01.7542131.754213105760364
0.10.32.01.7104891.710489423953702
0.10.20.11.2015561.201559983439274
0.10.20.21.2372231.237224345281889
0.10.20.31.2718161.271816653083256
Table 2. Table of present values of  2 R e C f  compared against MATLAB’s bvp4c results for selected values of  λ  and  S = 0.1 .
Table 2. Table of present values of  2 R e C f  compared against MATLAB’s bvp4c results for selected values of  λ  and  S = 0.1 .
  ( f ( 0 ) + λ 2 ( f ( 0 ) ) 2 )
  λ Nadeem and Hussain [40]MATLAB’s bvp4cSQLM
0.01.329301.3293027360627211.329308462412963
0.11.298011.2980170712948071.298022829158239
0.21.263101.2631037960983371.263109548733657
0.31.222761.2227766171144271.222781418705920
Table 3. The numerical values of the skin friction coefficient, Nusselt number and Sherwood number for all the thermo-physical parameters.
Table 3. The numerical values of the skin friction coefficient, Nusselt number and Sherwood number for all the thermo-physical parameters.
  λ SMPr   G C   G T   R d   N b   N t   δ ScKr   2 Re C f   Nu x Re x   Sh x Re x
0.10.22.00.50.10.10.10.50.50.11.00.11.7767670.3082820.700530
0.3 1.7064940.2954630.687057
0.9 1.6089960.1695150.539204
0.10.1 1.7300040.3239070.736618
0.2 1.7767680.3082820.700530
0.3 1.8249880.2940090.665454
0.20.1 1.1742190.4341670.794308
0.2 1.2131710.4260830.787275
0.3 1.2510270.4181690.780543
2.00.1 1.7719390.1472230.810761
0.2 1.7732170.1915190.780377
0.3 1.7744610.2330310.751949
0.50.3 1.6355860.3472270.724385
0.6 1.4305630.3913370.757335
0.9 1.2323560.4246970.786362
0.10.3 1.4742630.3684640.735159
0.6 1.2460020.4194820.777179
0.9 1.0279070.4546210.811112
0.10.3 1.4709570.2736910.724148
0.5 1.4684210.2488890.741101
0.7 1.4664120.2302550.753849
0.10.1 1.3904810.3163430.032976
0.3 1.4616520.3123010.593189
0.5 1.4742630.3082820.700530
0.50.1 1.4913100.3583200.807764
0.3 1.4826460.3245440.753061
0.5 1.4742630.3082820.700530
0.10.0 1.4773660.3825480.658192
0.1 1.4742630.3082820.700530
0.2 1.4702400.1987530.756672
0.10.6 1.4465750.3223040.458305
0.8 1.4618640.3142010.587691
1.0 1.4742630.3082820.700530
1.00.21.4797270.3037870.778334
0.41.4885520.2979580.906974
0.61.4955330.2941691.015232
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Muzara, H.; Shateyi, S. Magnetohydrodynamics Williamson Nanofluid Flow over an Exponentially Stretching Surface with a Chemical Reaction and Thermal Radiation. Mathematics 2023, 11, 2740. https://doi.org/10.3390/math11122740

AMA Style

Muzara H, Shateyi S. Magnetohydrodynamics Williamson Nanofluid Flow over an Exponentially Stretching Surface with a Chemical Reaction and Thermal Radiation. Mathematics. 2023; 11(12):2740. https://doi.org/10.3390/math11122740

Chicago/Turabian Style

Muzara, Hillary, and Stanford Shateyi. 2023. "Magnetohydrodynamics Williamson Nanofluid Flow over an Exponentially Stretching Surface with a Chemical Reaction and Thermal Radiation" Mathematics 11, no. 12: 2740. https://doi.org/10.3390/math11122740

APA Style

Muzara, H., & Shateyi, S. (2023). Magnetohydrodynamics Williamson Nanofluid Flow over an Exponentially Stretching Surface with a Chemical Reaction and Thermal Radiation. Mathematics, 11(12), 2740. https://doi.org/10.3390/math11122740

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop