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

In this article, the e ﬀ ect of electromagnetic force with the e ﬀ ect of thermal radiation on the Williamson nano ﬂ uid on a stretching surface through a porous medium was studied considering the e ﬀ ect of both heat generation/absorption and Joule heating. On the other hand, the e ﬀ ect of Brownian motion and thermophoresis coe ﬃ cients was considered. The system of nonlinear partial di ﬀ erential equations governing the study of ﬂ uid ﬂ ow has transformed into a system of ordinary di ﬀ erential equations using similarity transformations and nondimensional variables which were subsequently solved numerically by using the Rung-Kutta fourth-order method with shooting technique. Moreover, the e ﬀ ect of the resulting physical parameters on the distributions of velocity, temperature, and concentration of nanoparticles has been studied by using graphical forms with an interest in providing physical meanings to each parameter. Finally, special diagrams were made to explain the study of the e ﬀ ect of some physical parameters on the skin friction coe ﬃ cient and the local Nusselt number; these results led to reinforcement in the values of the skin friction coe ﬃ cient for the increased values of the magnetic ﬁ eld and the Darcy number while the e ﬀ ect on the local Nusselt number by thermal radiation as well as the heat generation/absorption coe ﬃ cients became negative.


Introduction
In recent years, the study of non-Newtonian fluids has received the attention of researchers in the field of hydrodynamics around the world due to the enormous scientific developments in their applications. The Williamson fluid is one of the most important non-Newtonian fluids characterized by less viscosity with an increase in the rate of shear stress and very similar in its properties of polymeric solutions, for example. In another meaning, in the Williamson fluid model, the effective viscosity should lessen indefinitely with the rising shear rate, which is nothing but infinite viscosity at stationary and nil viscosity as the shear rate tends to infinity. A model of Williamson has been discovered by Williamson [1] in 1929, while Subbarayudu et al. [2] investigated the assessment of time-dependent flow of Williamson fluid with radiative blood flow against a wedge. On the other hand, Lyubimova et al. [3] analyzed the stability of quasiequilibrium states and supercritical regimes of thermal vibrational convection of Williamson fluid in zero gravity conditions, but multiple solutions for MHD transient flow of Williamson nanofluids with convective heat transport were studied by Hashim et al. [4]. Hamid et al. [5] discussed an investigation of thermal and solutal stratification effects on mixed convection flow and heat transfer of Williamson nanofluid.
Nanofluids are a modern class of fluids discovered by Choi [6]. Researchers in the field of studying fluids around the world have paid attention to this type of fluids, which plays an important role in many modern technological applications that serve life. On the other hand, nanofluids enter many industries, for example, the satellite industry. Nanofluids also have an important role in the medical fields, for example, also in the use of nanoparticles of gold in the treatment of cancerous tumors, as well as the manufacture of microscopic bombs that are used to eliminate cancerous tumors as well. Sajadifar et al. [7] examined fluid flow and heat transfer of non-Newtonian nanofluid in a microtube considering slip velocity and temperature jump boundary conditions, while Khan and Pop [8] investigated boundarylayer flow of a nanofluid past a stretching sheet, and at the same time, Farooq et al. [9] analyzed MHD stagnation point flow of viscoelastic nanofluid with nonlinear radiation effects. Very recently, Ahmed et al. [10] explained magnetohydrodynamic Maxwell nanofluid flow over a stretching surface through a porous medium: effects of nonlinear thermal radiation, convective boundary conditions, and heat generation/absorption. Alshomrani [11] studied generalized Fourier's and Fick's laws in bioconvection flow of magnetized Burgers nanofluid utilizing motile microorganisms. Alshomrani [12] investigated numerical investigation for bioconvection flow of viscoelastic nanofluid with magnetic dipole and motile microorganisms. Loganathan and Rajan [13] excogited an entropy approach of Williamson nanofluid flow with Joule heating and zero nanoparticle mass flux. It was found from this study that the velocity distribution decreases with increasing the values of the magnetic field parameter while the temperature distribution increases with increasing the values of this parameter. Khan et al. [14] illustrated evaluating the characteristics of the magnetic dipole for shearthinning Williamson nanofluid with thermal radiation. Khan et al. [15] represented changes in viscosity of Williamson nanofluid flow due to thermal and solutal stratification. Hayat et al. [16] explained mixed convective three-dimensional flow of Williamson nanofluid subject to chemical reaction. Ramzan et al. [17] analyzed Darcy-Forchheimer 3D Williamson nanofluid flow with generalized Fourier and Fick's laws in a stratified medium. It was observed from this study that the velocity distribution was decreased by the influence of the magnetic field parameter and also by the influence of the Williamson fluid parameter.
Electromagnetic fluid dynamics (magnetohydrodynamic fluid) is the study of fluids that have the characteristic of electromagnetic conductivity such as plasma fluids, metal fluids, and salt water, and the first to start studying this field is Alfvén [18] in 1942. Mohamed et al. [19] studied MHD Jeffrey nanofluid flow over a stretching sheet through a porous medium in the presence of nonlinear thermal radiation and heat generation/absorption. Also, Chandrashekar et al. [20] examined a discontinuous Galerkin method for a twodimensional reduced resistive MHD model. A numerical study of MHD mixed convection under volumetric heat source in a vertical square duct with wall effects investigated by Liu et al. [21]. Tassone et al. [22] excogitated MHD mixed convection flow in the WCLL: heat transfer analysis and cooling system optimization. Motsa et al. [23] illuminated the spectral relaxation method and spectral quasilinearization. Khan et al. [24] studied peristaltic transport of a Jeffrey fluid with variable viscosity through a porous medium in an asymmetric channel. Ellahi [25] elucidated a study on the convergence of series solutions of non-Newtonian thirdgrade fluid with variable viscosity by means of the homotopy analysis method. Umar et al. [26] investigated numerical treatment for the three-dimensional Eyring-Powell fluid flow over a stretching sheet with velocity slip and activation energy. Ramesh et al. [27] explained heat transfer in MHD dusty boundary layer flow over an inclined stretching sheet with nonuniform heat source/sink.
Joule heating (also referred to as resistive or ohmic heating) describes the process where the energy of an electric current is converted into heat as it flows through a resistance. The explanation for this is that when electrical current flows through solid or liquid materials that conduct electricity, electrical energy is converted into heat energy through the occurrence of resistances inside the conductor. In this case, free electrons transfer energy by collisions. Ramzan et al. [28] studied radiative and Joule heating effects in the MHD flow of a micropolar fluid with partial slip and convective boundary condition. Mohamed et al. [29] discussed the thermal radiation and MHD effects on free convective flow of a polar fluid through a porous medium in the presence of internal heat generation and chemical reaction. Also, Tetbirt et al. [30] discussed a numerical study of magnetic effect on the velocity distribution field in a macro-/microscale of a micropolar and viscous fluid in a vertical channel. At some time, Ghadikolaei et al. [31] investigated numerical study on magnetohydrodynamic CNT-water nanofluids as a micropolar dusty fluid influenced by nonlinear thermal radiation and joule heating effect. It is worth noting that Gireesha et al. [32] investigated hall effects on dusty nanofluid twophase transient flow past a stretching sheet using the KVL model. It was also reported that Aghanajafi et al. [33] explained numerical simulation of laminar forced convection of water-CuO nanofluid inside a triangular duct. But Hussain et al. [34] examined the effects of viscous dissipation and Joule heating on MHD Sisko nanofluid over a stretching cylinder. Moreover, free convective heat and mass transfer of MHD non-Newtonian nanofluids over a cone in the presence of a nonuniform heat source/sink was studied by Raju et al. [35]. Hayat et al. [36] excogitated radiative flow of Jeffrey fluid in a porous medium with power law heat flux and heat source. Meanwhile, Hartnett and Kostic [37] analyzed heat transfer to Newtonian and non-Newtonian fluids in rectangular ducts.
The main objective of the present work is to study the flow of a non-Newtonian Williamson fluid that contains nanoparticles on a stretching sheet through a porous medium under the influences of the magnetic field, nonlinear thermal radiation, and Joule heating in the presence of heat generation/absorption and chemical reaction on the distributions of 2 Advances in Mathematical Physics velocity, temperature, and concentration of nanoparticles taking into account studying effects of the Brownian motion coefficient and thermophoresis coefficient. On the other hand, the similarity transformations and nondimensional variables were used in converting the system of partial differential equations governing the movement of flow into a system of ordinary differential equations, which were solved by using the Runge-Kutta numerical method with a shooting technique. The graphs were used to study the effects of all physical parameters on the distributions of velocity, temperature, and concentration of nanoparticles.

Formulation of the Problem
In the beginning, the study of Williamson nanofluid flow can be shown in several steps: the first step is the fluid flow being steady and incompressible, the second step is that the flowing process is two-dimensional flow on a stretching surface through a porous medium, and the third step is the plate is stretched along the x-axis with a velocity U w ðxÞ = Bx, where B is the stretching rate. The basic equations for the balance of mass, momentum, energy, and nanoparticle volume fraction of the flow problem can be expressed in vector form as follows.

Continuity Equation.
The continuity equation for conservation of mass becomes where ∇ is the differential operator and q ! is the flow velocity vector.

Conservation of Momentum Equation.
The Navier-Stokes equation for the balance of linear momentum is given by where ρ f is the density of the nanofluid, t is the time, S is the Cauchy stress tensor, μ is the coefficient of dynamic viscosity, is the Lorentz force produced by the interaction of the applied magnetic field with velocity of fluid, J is the external magnetic field, σ is the electrical conductivity, and K is the permeability of the porous medium.

Conservation of Energy Equation.
In the absence of viscous dissipation effects, the conservation of energy for heat transfer is given by the constitutive equation as where c P is the specific heat at constant pressure; ðρc P Þ f and ðρc P Þ p are the specific heat of the nanofluid and the nanoparticles, respectively; T is the temperature of the nanofluid; k is the thermal conductivity; D B is the Brownian diffusion coefficient; D T is the thermophoresis diffusion coefficient; T ∞ is the ambient temperature; q r is the radiative heat flux; J h = F ! e · q ! is the Joule heating, and Q is the uniform volumetric heat generation/absorption.

Conservation of Nanoparticle Concentration Equation.
For a homogeneous chemical reaction [38], the concentration equation of nanoparticle volume fraction becomes where C is the nanoparticle volume fraction and the rates R * > 0 and R * < 0 denote destructive and constructive reaction rates, respectively. For the Williamson fluid model, Cauchy stress tensor S is defined as [26] S = −P I + τ, ð5Þ where P is pressure, I is the unit tensor, τ is the extra stress tensor, μ 0 is limiting viscosity at zero shear rate, μ ∞ is limiting viscosity at the infinite shear rate, Γ > 0 is a time constant, A 1 is the first Rivlin-Erickson tensor, and γ : is defined as follows: Here, it is considered the case for μ ∞ = 0 and Γ γ : < 1; thus, Equation (6) can be written as follows: The right-hand side of Equation (8) is done by using binomial expansion. Making use of Equations (5) and (8) into Equations (1) to (4), the two-dimensional boundary layer equations governing the flow can be written as follows [39]:

Advances in Mathematical Physics
where u is the velocity components in the direction of the x -axis and v is the velocity components in the direction of the y-axis; in addition, α represent the thermal diffusivity, ρ refer to density of the fluid, and ν is the kinematic viscosity of fluid; also, T and T ∞ , respectively, are fluid temperature and ambient fluid temperature. A uniform magnetic field of strength B 0 is applied in the transverse direction of the flow; due to the small magnetic Reynolds number, it is not necessary to introduce the effect of the induced magnetic field. It should be noted that D B , D T , σ, c P , and R * are, respectively, the Brownian diffusion coefficient, the thermophoresis diffusion coefficient, the electrical conductivity, the specific heat at constant pressure, and destructive and constructive reaction rates. Here, C and C ∞ are the concentration of nanoparticles and ambient nanoparticle concentration. It is worth noting that the effect of nonlinear thermal radiation, heat source/sink, and Joule heating is taken in the energy equation. The viscous dissipation is assumed to be negligibly small in the energy equation. The effects of the homogenous chemical reaction are taken in the concentration equation. The boundary conditions for the present investigation problem can be written in the formula: where U w ðxÞ = Bx represent the stretching surface velocity and B > 0 is the stretching rate, and the radiation heat flux q r can be written as follows: when substituting by Equation (16) into Equation (11), the energy equation takes the following form: where σ * and k * are the Stefan-Boltzmann constant and mean absorption coefficient; the similarity transformations and nondimensional variables can be written as follows: where η is the similarity variable and f ðηÞ, θðηÞ, and ϕðηÞ, respectively, are the dimensionless stream function, temperature, and concentration of nanoparticles; the similarity transformations and the nondimensional variables (18) were used to convert the boundary layer governing partial differential Equations (10), (12), and (17) to a set of ordinary differential equations taking the mathematical formulas: where the boundary conditions (13) and (14) can be written after converting as new follows: with the knowledge that M = σB 0 2 /ρB is the magnetic field parameter, R = 16σ * T 3 ∞ /3kk * is the nonlinear thermal radiation parameter, Pr = ν/α is the Prandtl number, Advances in Mathematical Physics Expression for the local Nusselt number Nu x and the skin friction coefficient C f are defined as And the reader will note that the dimensionless mass flux represented by a Sherwood number Sh x is now identically zero and q w and τ w are the heat flux and the shear stress along the stretching surface, respectively, on the following mathematical formulas: When substituting Equations (26) and (27) into Equations (24) and (25), the mathematical nondimensional formulas of Equations (24) and (25) are ffiffiffiffiffiffiffi where Re x = xU w ðxÞ/ν is the local Reynolds number based on the stretching velocity U w ðxÞ. It should be noted that in the absence of a magnetic field, porous medium, nonlinear thermal radiation, heat generation/absorption, and finally the chemical reaction, the work is due to Nadeem and Hussain [39]. Nt

Numerical Solution
The Runge-Kutta numerical method is one of the most important numerical analysis methods that are used in solving a system of ordinary differential equations through which the differential equation is reduced degree to make it easier and simpler to solve. On the other hand, there are different formulas for the solution by the numerical method of Runge-Kutta method; for example, there is the Runge-Kutta method of the fourth order and also the Runge-Kutta method of the fifth order, and the most used method is the Runge-Kutta method of the fourth order because it gives accurate results and is easy to use, and its derivation depends on the Euler method. The numerical solutions of the current study are listed in several steps. In the first step, the system of partial differential Equations (10), (12), and (17) for the current study has been converted into a system of ordinary differential Equations (19)-(21) using the similarity transformation and nondimensional variables (18); also, the boundary conditions (13) and (14) have been converted into the new forms (22) and (23). In the second step, the appropriate numerical method used in solving the new system of prior equations is the Runge-Kutta fourth-order with shooting technique, and by the way, this method was used to reduce the order of differential equations. In the third step, the program used to solve the system of ordinary differential equations is MATLAB, and the step size Δη = 0:001 is used to obtain the numerical solution with η max = 10; also, it used the bvp4c function in solving these equations. Differential Equations (19)- (21) can be rewritten in the simplest form as follows: In the last step, Equations (30)- (32) have been written inside the MATLAB program on the following formulas: Advances in Mathematical Physics With the boundary conditions: Taking into account both: To ensure the accuracy and correctness of the numerical solutions of the current study, a numerical comparison was made between the numerical values and results of the current study with the numerical values and results of the work published by Khan and Pop [8] in Table 1. The great convergence between the two studies was noted, which gives high credibility to the current study.

Results and Discussion
After converting the system of partial differential equations ruling to study the flow of fluid into a system of ordinary differential equations, it has a set of important parameters that we list in the following order M, Da, and λ which are the magnetic field, the Darcy number, and non-Newtonian Williamson parameters, respectively, while R, θ w , and S represent the nonlinear thermal radiation, the ratio temperature, and the heat generation/absorption parameters, respectively; also, Pr is the Prandtl number and Nt is the thermophoresis parameter, while Ec, γ, Le, and Nb are the Eckert number, the chemical reaction, the Lewis number, and Brownian motion parameters, respectively. The effect of all the previous physical parameters on the velocity, temperature, and concentration of nanoparticle distributions has been studied by making graphical figures that clarify this and by showing the physical meanings of each parameter and its importance in this study, and we can list the results of this study in detail.

Velocity Distributions.
It is known that the rate at which the fluid flows onto a specific surface is entirely determined by the distribution of velocity; therefore, velocity distributions play an important role in studying the changes that occur in the flow rate of the fluid. On the other hand, when affecting the flowing fluid with external forces, the behavior of the fluid's movement is changing, and one of these types of external forces is the magnetic field. Figure 1 shows the effect of the strength of the magnetic field parameter M on the fluid velocity distribution f ′ ðηÞ. It is noticeable that the result of this effect is negative in the sense that the relationship between the magnetic field and the velocity distribution is inverse; the increase in the values of the magnetic field means a decrease in the velocity of the fluid. Physically, when influencing a moving fluid by the magnetic field, the fluid particles are stimulated, which creates a kind of counter force that slows and reduces the fluid's motion; moreover, this force is perpendicular to the velocity vector on the one hand and also perpendicular to the magnetic field vector on the other hand which is originally a resistance force called Lorentz force, while highlighting that the increase in the magnetic field reduces the thickness of the boundary layer. If you want to know the effect of the Darcy number Da on the velocity distribution, you should look at Figure 2. It becomes clear to you that the velocity distribution decreases under the influence of the large values that the Darcy number Da takes. Physically, the porous medium is the medium that contains a group of small voids called pores that are interrupted by the fluid when it moves on this medium. On the other hand, the permeability of the fluid through the porous medium is related to its porosity.
Also, the greater the values of the Darcy number Da, the greater the resistance of the porous medium to the movement of the fluid on it on the one hand, in addition to the fluid viscosity on the other hand, which leads to a decrease in the fluid velocity. Moving to Figure 3, it reviews the effect of the non-Newtonian Williamson coefficient λ on the velocity distribution, and it clears that this effect is negative in the sense that the increase in the values of this coefficient is followed by a decrease in velocity of the fluid, which makes its movement slow.

Temperature Distributions.
The change in the fluid temperature plays an important role in the fluid's behavior in addition to its effect also on the particles inside the fluid. Now, we will study the effect of some physical parameters on the fluid temperature distribution θðηÞ and note the changes taking place; for example, Figures 4 and 5 illustrate the effects of the nonlinear thermal radiation parameter R   Advances in Mathematical Physics and the ratio temperature parameter θ w on the temperature distribution θðηÞ; it was noticed that the large values taken by the nonlinear thermal radiation parameter R and ratio temperature parameter θ w work to enhance the heat transfer inside the fluid and work to stimulate it. Physically, the positive effect of the nonlinear thermal radiation parameter R on the temperature of the nanofluid leads to three things. Firstly, the heat transfer of all distances of the boundary layer is supported, or in other words, it raises the temperature of the boundary layer regularly. Secondly, it makes the nanoparticles inside the fluid gain thermal energy, which improves the transfer and thermal diffusion within the fluid due to the thermal conductivity of the nanoparticles.
Thirdly, it works to enhance the thermal transfer methods of the nanofluid, which is the thermal transfer method by conduction and the method of thermal transfer by load. On the other hand, the effect of the ratio temperature increases the boundary layer temperature and increases its thickness. When talking about the effect of the heat generation/absorption parameter S on the fluid temperature distribution θðηÞ, as shown in Figure 6, there is a significant enhancement in the distribution of the fluid temperature due to the enhancement in the values of parameter S.
Physically, in the case of the phenomenon of heat generation, there is an enhancement in the transfer and thermal spread of the fluids, which raises temperature of fluid and also raises the temperature and thickness of the boundary layer with the increase in the values of the heat source type ðS > 0Þ, while the exact opposite occurs in the case of the heat absorption type ðS < 0Þ, but when moving upward from the state of heat absorption type to the state of heat generation type, there is an improvement in the rate of thermal diffusion and the thickness of the boundary layer. Figure 7 displays the impact of Prandtl number Pr on the temperature distribution θðηÞ; when this parameter takes a large value, the fluid temperature distribution is decreased. Physically, Prandtl number Pr is defined as the ratio of momentum diffusivity to thermal diffusivity; on the other hand, there are two types of Prandtl number Pr: the first type is Pr > 1; this means that the rate of momentum diffusion dominates the rate of thermal diffusion, and the second type is Pr < 1 which urges the opposite, and Figure 7 shows the case in which Pr > 1 such that the thermal diffusion within the fluid becomes small; thus, the fluid temperature and thickness of the boundary layer become decreasing. The distribution of the fluid temperature θðηÞ changes with the effect of the thermophoresis parameter Nt values; as shown in Figure 8, the direct relationship between them refers to the enhancement in the values of the thermophoresis parameter Nt leading to an     Advances in Mathematical Physics increase in the fluid temperature distribution θðηÞ. Physically, there are two different sources: the first is hot and the other is cold; the particles close to the hot medium absorb heat energy, which makes them move from the hot medium to the cold medium and raises its temperature by means of convection; this process is called the process of thermal potential difference; this is exactly what happens inside the fluid when it is exposed to a heat source. The fluid particles absorb heat energy that makes it spread throughout the fluid, raising its temperature as well as enhancing the thickness of the boundary layer of the fluid. Figure 9 studies the effect of the Eckert number Ec on the temperature distribution θðηÞ, so it was noted that the large values that take by this parameter affect the temperature distribution making this effect positive in the sense that it enhances the effect of temperature. Finally, the influence of the coefficient of the magnetic field M as well as the effect of the Darcy number Da on the temperature distribution θðηÞ is shown graphically in Figures 10 and 11, respectively; it is very clear from the study of the two figures that the fluid temperature and the boundary layer thickness increase with increasing values of these two parameters.     Advances in Mathematical Physics the effect is very positive. Physically, firstly, the effect of the magnetic field M on the fluid works to reduce the movement of the fluid accompanied slowly in the movement of the nanofluid particles, which makes its concentration large. Secondly, the effect of the Darcy number Da on the fluid movement gives a great opportunity for collision of the nanoparticles with the pores of the porous medium, which causes it to accumulate inside the pores and the porous medium itself, and accordingly, the concentration of nanoparticles becomes increasing. On the other hand, Figures 14 and 15 show the effects of the nonlinear thermal radiation parameter R and the ratio temperature parameter θ w on the concentration of nanoparticle distribution ϕðηÞ within the fluid. It was noted that the reinforcement in the values of these two parameters works to strengthen the temperature of the nanoparticles and thickness of the boundary layer of the fluid, as it works on a significant increase in the spread of the temperature in all parts of the fluid. Physically, nanoparticle fluid molecules are greatly affected by the radiation source that generates heat, which causes fluid particles to interact positively with this heat, which increases the mechanical energy of the nanoparticles, and this actually leads to their stimulation and increased concentration. The concentration of nanoparticles gives a good or positive impression in the event that it is affected by the heat generation/absorption parameter S because the mechanism of action of this parameter is very similar to the mechanism of nonlinear thermal radiation parameter R only in the case of heat generation; on the other side, the heat generation works to stimulate the temperature of nanoparticles and thus increases the concentration of nanoparticles ϕðηÞ and adds to the positive effect on the

10
Advances in Mathematical Physics boundary layer in terms of thickness, while in the case of heat absorption, the opposite occurs, and this is shown in Figure 16. As for Figure 17, it shows the negative effect of Prandtl number Pr on the concentration of nanoparticle distribution ϕðηÞ because in the case where Pr > 1, thermal diffusion equipment is very small compared to the momentum diffusion rate, and accordingly, the concentration of nanoparticles becomes small that is the physical meaning. Figures 18 and 19 show the effects of both the Eckert number Ec and the thermophoresis parameter Nt on the distribution of the concentration of nanoparticles ϕðηÞ within the fluid; it was found that the increase in this parameters gives a large concentration of the nanofluid particles, at a time when the effects of both Lewis number Le and chemical reaction parameter γ on the concentration of nanoparticles ϕðηÞ inside the fluid were negative; this is illustrated in Figures 20 and 21. Finally, Figure 22 clears the effect of the Brownian motion parameter Nb on the distribution of the concentration of nanoparticles ϕðηÞ within the fluid. It is noticeable that the distribution of the concentration of nanoparticles ϕðηÞ decreases in the event that this parameter is enhanced. Physically, the Brownian motion of the nanofluid is a random motion of the nanoparticles present in this fluid in the absence of any external effects so that the movement of these particles is free in all directions, but in the case of a large thermal diffusion resulting from the effect of thermal radiation and the coefficient of thermophoresis, it becomes that the random movement of the nanoparticles is limited.

Conclusion
In this article, the two-dimensional electromagnetic flow of the Williamson nanofluid has been studied on a stretching sheet through a porous medium; on the other hand, the governing partial differential equations were converted into a system of ordinary differential equations using the similarity transformations and the nondimensional variable under the effects of nonlinear thermal radiation, heat generation/absorption, chemical reaction, Brownian motion parameter, thermophoresis parameter, magnetic field, Darcy number, Joule heating, Lewis number, Prandtl number, and Williamson non-Newtonian parameter on distribution of velocity, temperature, and concentration of nanoparticles; the most important points drawn from this study are the following:  for values of M, Da, and λ parameters when R = 1:2, Nb = 0:5, Nt = 0:5, Pr = 5, Le = 2, S = −0:1, Ec = 0:4, θ w = 1:2, and γ = 0:6. Thermal conductivity k * :

M
Mean absorption coefficient Le: Lewis number