Research Article | Open Access

Gossaye Aliy Adem, "Analytic Treatment for Electrical MHD Non-Newtonian Fluid Flow over a Stretching Sheet through a Porous Medium", *Advances in Mathematical Physics*, vol. 2020, Article ID 8879264, 14 pages, 2020. https://doi.org/10.1155/2020/8879264

# Analytic Treatment for Electrical MHD Non-Newtonian Fluid Flow over a Stretching Sheet through a Porous Medium

**Academic Editor:**John D. Clayton

#### Abstract

In this study, an attempt has been made to investigate the mass and heat transfer effects in a BLF through a porous medium of an electrically conducting viscoelastic fluid subject to a transverse magnetic field in the existence of an external electric field, heat source/sink, and chemical reaction. It has been considered the effects of the electric field, viscous and Joule dissipations, radiation, and internal heat generation/absorption. Closed-form solutions for the boundary layer equations of viscoelastic, second-grade, and Walters’ fluid models are considered. The method of the solution includes similarity transformation. The converted equations of thermal and mass transport are calculated using the optimal homotopy asymptotic method (OHAM). The solutions of the temperature field for both prescribed surface temperature (PST) and prescribed surface heat flux (PHF) are found. It is vital to remark that the interaction of the magnetic field is found to be counterproductive in enhancing velocity and concentration distribution, whereas the presence of chemical reaction, as well as a porous matrix with moderate values of the magnetic parameter, reduces the temperature and concentration fields at all points of the flow domain.

#### 1. Introduction

The significance of fluid flow over a stretching surface can be professed for its ever-increasing inexorable applications in industries and in present-day technology. The applications of the stretching sheet problem are such as polymer sheet extrusion from a dye, drawing, thinning and annealing of copper wires, glass fiber and paper production, and the cooling of a metallic plate in a cooling bath. The production of these sheets needs that the melt issues from a slit and is stretched to get the anticipated thickness. The final product depends on the rate of cooling in the process and the process of stretching. Sakadis [1] was the first to investigate the BL manners on a continuous solid surface moving with constant speed. Crane [2] was the first to attain a stylish analytical solution to the BL equations for the problem of steady 2D flow through a stretching surface in an inactive incompressible fluid. The work of Sakadis and Crane was extended and studied by several authors (Madaki et al. [3], Zheng et al. [4], Dessie and Kishan [5], and Pal [6]).

The rising need for chemical reaction and hydrometallurgical industries needs the study of heat and mass transfer with chemical reaction. There are several transport processes that are governed by the combined action of buoyancy forces due to both thermal and mass diffusion in the presence of chemical reactions and electric field effects. These processes occur in the nuclear reactor safety and combustion systems, solar collectors, metallurgical and chemical engineering, etc.

In real situations, most of the fluids used in engineering applications are almost non-Newtonian, specifically of the viscoelastic type than the viscous type. Khan [7], Madhu and Kishan [8], Babaelahi et al. [9], Raftari and Vajravelu [10], Goyal and Bhargava [11], and Madhu and Kishan [12] have investigated the flow of viscoelastic fluid flow under different conditions. Recently, Aliy and Kishan [13] examined the effect of electric field on MHD viscoelastic nanofluid flow and heat transfer over a stretching sheet with convective boundary conditions. They have considered the local nonsimilarity method and OHAM to solve the resulting equation. Jafar et al. [14] studied the MHD boundary layer flow of a viscoelastic fluid past a nonlinearly stretching sheet in the presence of a viscous dissipation effect.

Konda et al. [15] investigated the effect of nonuniform heat source/sink on MHD BLF and melting heat transfer of Williamson nanofluid in a porous medium. Hayat et al. [16] studied the effects of the chemical reaction of an unsteady 3D flow of the couple stress fluid over a stretching surface. Gireesha et al. [17] have studied the BLF and heat transfer of a dusty fluid flow over a stretching sheet in the existence of nonuniform heat source/sink and radiation. Parsa et al. [18] investigated the MHD boundary-layer flow over a stretching surface with internal heat generation or absorption. The thermal radiation effect was studied by different researchers (Weidman [19], Abel and Nandeppanavar [20], and Madaki et al. [21]).

OHAM is an approximate analytical technique that is directed forward to employ on a different type of problem, and the presence of any small or large parameters is not important. Marinca et al. [22] initially introduced the basic concept of this method in 2008. OHAM reduces the extent of the computational domain. It is a reliable analytical technique and has already been successfully applied to various nonlinear coupled differential equations occurring in science, engineering, and other fields of study. Marinca and Herisanu applied OHAM on the different problems (see, [23–25]). Many researchers applied OHAM to study fluid flow problems [26–29]. Recently, Gossaye and Kishan [13, 30, 31] applied the OHAM to the problem arisen from electrical MHD non-Newtonian fluid flows over a stretching sheet, and Khan et al. [32] used OHAM for Lane-Emden and Emden-Fowler initial and boundary value problems.

The novelty of the present study is starred due to the following aspects. (i)Inclusion of electric field is justified since, for a remarkable enhancement in the value of the electric field, the struggle among fluid particles rises, and hence, the Lorentz force inclines to speed up the body forces, and it leads to escalating in the flow velocity and thicker momentum boundary layer. And also, we observed that the concentration of the fluid is affected by a large amount of an electric field.(ii)The dimensionless ordinary differential equations with the corresponding boundary conditions are solved by using the approximate analytic technique OHAM.

The aim of this paper is to investigate the effect of electric field and thermal radiation on MHD viscoelastic fluid flow over a stretching surface through a porous medium with a heat source/sink. The strongly nonlinear differential equations are solved using the approximate analytic method OHAM. The effects of different governing parameters on velocity, temperature, and concentration profiles are displayed graphically and argued extensively. Moreover, the skin-friction coefficient is also exhibited using a table.

#### 2. Formulation of the Problem

Consider a steady 2D boundary layer flow of an electrically conducting, viscoelastic fluid past a stretching surface embedded in a porous medium, the flow being confined to. Two equal and opposite forces are applied along the axis, and hence, the surface is stretched, keeping the origin fixed. Assuming that a uniform magnetic field of strength is applied along that produces the magnetic effect in the direction.

The Lorentz force is given by , where , the Joule current can be expressed as , stands for the electrical conductivity, denotes the velocity of the fluid, and and are the transverse magnetic field and electric field vectors, correspondingly. The magnetic field and electric field are applied perpendicular to the flow, and hence, the Reynolds number is selected small. The induced magnetic field is insignificant for small magnetic Reynolds number because it is smaller when compared to the applied magnetic field. Under the usual BL conventions, the equations of continuity, momentum, energy, and concentration for the flow of viscoelastic fluid are:

The approximation of Rosseland for thermal radiation gives . It is assumed that the temperature variation within the flow is such that may be expanded in a Taylor series. Expanding about and neglecting the higher-order terms, we have

Taking the above equation and substituting it into equation (3), one can get

The boundary conditions are

#### 3. Solution of the Flow Field

Equations (1) and (2) disclose self-similar solutions of the type given below where represents the dimensionless stream function, denotes the similarity variable, represents the local Reynolds number, and symbolizes the stream function satisfying the continuity equation (1). Substituting equation (8) in equation (2), we get where is the viscoelastic parameter, is the magnetic parameter, is the electric field parameter, and is the nondimensional permeability parameter.

The BCs are: where is the suction/injection parameter, and and represent suction and injection, respectively.

Now, we apply OHAM to the nonlinear ODE (9) with the BCs (10) under the following assumption

where the parameter , is an auxiliary function which is different from zero, and are constants.

And also, assume the following where and are the nonlinear and linear operators, respectively. Hence, the OHAM family of equations is given by

Using the BCs (10) and comparing the coefficients of the same powers of -terms, we can get the hypothesis below:

Equating the coefficients of the zero-order equation , we obtain

Equating the coefficients of the first-order equation , we obtain

Equating the second-order equation , we obtain

After solving the ODEs (14) and (15) with the corresponding boundary conditions, we obtain

We cannot write the third term here because of its being large. Hence, the solution in terms of the parameters can be written as:

For equation (19), we define its residual equation in the form given below

The unknown convergence parameters will be optimally well-known using the conditions given below

Taking particular cases when , the convergence parameters are as follows

Hence, the approximate analytical solution is given as

The shear stress at the wall is defined as

The nondimensional form of skin friction coefficient at the well is

#### 4. Heat Transfer Analysis

##### 4.1. Case I: Prescribed Surface Temperature (PST)

In the prescribed surface temperature case, introducing nondimensional quantities as given below and using equation (7) equation (6) become with the BCs given below

Now let us apply the OHAM to the nonlinear ODE (27) using the following assumption where and are the nonlinear and linear operators, respectively. Hence, the OHAM family of equations is given by

Using the BCs (7) and comparing the coefficients of the same powers of -terms in the above equation, one can have the following:

Equating the coefficients of the zero-order equation , we obtain

Equating the coefficients of the first-order equation , we get

Equating the coefficients of the second-order equation , we get

.and .

After solving the ODEs (31), (34), and (36) with the corresponding boundary conditions, we obtain

We cannot write the third-term here because of it being large. Hence, the solution in terms of the parameters can be written as:

For equation (36), we define its residual equation in the form

The unknown convergence parameters can be optimally obtained from the following expressed conditions given below

In the particular case when , one can obtain the results of the convergence parameters. Hence, the values are calculated and obtained below

After substituting all the parameters, we get

The local Nusselt number for PST case is given by

##### 4.2. Case II: Prescribed Heat Flux (PHF)

In the prescribed heat flux case, introducing the similarity variable and using equation (7), equation (6) becomes with the boundary conditions

Now, let us apply the OHAM to the nonlinear ODE (42) using the following assumption where and are the nonlinear and linear operators, respectively. Hence, the OHAM family of equations is given by

Using the BCs (7) and comparing the coefficients of the same powers of -terms in the above equation, one can have the following:

Equating the coefficients of the zero-order equation , we obtain

Equating the coefficients of the first-order equation , we get

Equating the coefficients of the second-order equation , we get

and .

After solving the ODEs (46), (47), and (48) with the corresponding boundary conditions, we obtain

We cannot write the third term here because of it being large. Hence, the solution in terms of the parameters can be written as:

For equation (50), we define its residual equation in the form

The unknown convergence parameters can be optimally obtained from the following expressed conditions given below

In a particular case when , one can obtain the results of the convergence parameters. Hence, the values are calculated and obtained below

After substituting all the parameters, we get

The local Nusselt number for the PHF case is given by

#### 5. Mass Transfer Analysis

Introducing the similarity variables and using equation (7) in equation (4), we obtain the following with boundary condition

Now, let us apply the OHAM to the nonlinear ODE (56) using the following assumption

where and stands for the linear and nonlinear operators. Consequently, the OHAM family of equations is written as

Using the BCs (7) and comparing the coefficients of the same powers of -terms, one can have the following:

Equating the coefficients of the zero-order equation , one can obtain

Equating the coefficients of the first-order equation , we get

Equating the coefficients of the second-order equation , we get

After solving the ODEs (60), (61), and (63) with the corresponding boundary conditions, we obtain

We cannot write the third term here because of it being large. Hence, the solution in terms of the parameters can be written as:

For equation (65), we define its residual equation in the form

The unknown convergence parameters can be optimally obtained from the following expressed conditions given below

In the particular case when , one can obtain the results of the convergence parameters. Hence, the values are calculated as

After substituting all the parameters, we get the approximate analytic solution as given below

The local Sherwood number can be expressed as

#### 6. Result and Discussion

In the course of the discussion, the following aspects are highlighted: (i)Effect of electric field and permeability of the medium on flow characteristics(ii)Effect of diffusion species as well as the first-order chemical reaction(iii)Relative response of two viscoelastic models to temperature and velocity distribution in the existence of uniform porous matrix(iv)The dimensionless ODEs (9), (27), (36), and (56) with the corresponding boundary conditions are solved by using the approximate analytic technique OHAM

It is important to note that represent second grade, Walters , and viscous fluids, respectively.

Figure 1 displays the velocity distribution of second-grade fluid for suction (), injection (), and impermeable plate () in the presence () and absence () of the porous matrix. The porous medium () lessens the primary velocity at all points as a result of resistive force accessible by the porous medium which results in thinning of BL. Further, it is motivating to note that the suction at the plate decreases the velocity. Consequently, it is decided that the combined effect of porous matrix and suction is encouraging the thinning of BL which favors the stability of the flow.

Figure 2 displays the velocity profiles for Walters viscous () and second-grade flow (). It is perceived that velocity achieves low value in the case of Walters’ model, in the presence of a porous medium. Viscoelastic flows are disposed to instabilities as a result of nonlinearity in the constitutive equations. Instabilities are mainly driven by the fluid normal stresses (elasticity) or by the nature of the boundary conditions. Therefore, the elastic property of Walters’ flow model in combination with the permeability of the porous medium reduces the BL thickness and hence decreases the instability.

The variation of the velocity field for various values of is illustrated in Figure 3. One can easily understand from the figure that the velocity significantly declines with a rise in the values of in the absence of an electric field (). It is clear that the magnetic field depends on the Lorentz force, which is stronger for a greater magnetic field. Due to the nonexistence of an electric field, the Lorentz force raises the frictional force, which acts as a delaying force that opposes the viscoelastic fluid flow.

The effect of on is portrayed in Figure 4. As the magnitudes of boost up, the velocity BL rises nearer to the stretching plate. For remarkable enhancement in the value of , the struggle among the fluid particles rises, and hence, the Lorentz force inclines to speed up the body forces, and it leads to escalating in the flow velocity and thicker momentum boundary layer. The electric field has the same impact for the case of injection and suction as demonstrated in the figure.

Figures 5(a) and 5(b) display the temperature distribution in cases of PST and PHF cases, respectively, without suction/ injection. The effect of elasticity of the fluid subject to the present study on temperature distribution is opposite to that of the viscoelastic parameter that is the temperature declines at all points, and this is further contributed by the Walters model. Both the cases of PHF and PST show a similar effect.

**(a)**

**(b)**

Figures 6(a) and 6(b) demonstrate the effects of the Eckert number for second-grade fluid in the presence of suction/injection in both PST and PHF cases. The dissipation of energy of the flow pattern is measured by the Eckert number. It is seen that an increase in increases the temperature and hence increases the thermal BL thickness. This leads to a reduction in the rate of heat transfer from the plate surface. The same effect is observed in the case of injection and suction.

**(a)**

**(b)**

Figures 7(a) and 7(b) describe the effect of the radiation parameter in the case of second-grade fluid. It is perceived that an increase in increases the temperature of the fluid layer and the processes get accelerated in both cases of suction and injection. The rise in temperature with an increase in radiation parameter causes growth in temperature gradient of the wall in both PST and PHF cases.

**(a)**

**(b)**

Figures 8(a) and 8(b) display the outcome of the internal heat generation/absorption parameter on the temperature distribution (PST) and (PHF) in the case of second-grade fluid. This displays that a rise in heat source strength () increases the temperature. This is because of the generation of the heat in thermal BL which causes the temperature to escalate. In the same way, the heat sink causes temperature absorption resulting in a decline in temperature. The outcome holds good for both PHF and PST. The role of a porous matrix is to accelerate/decelerate the process in case of source/sink, respectively.

**(a)**

**(b)**