# Magnetohydrodynamic Fluctuating Free Convection Flow of Second-Grade Fluid Flow in a Porous Medium

Authors: Farhad Ali; Momin Khan; Madeha Gohar
Journal: Mathematical Problems in Engineering
 (2021)

---

## Abstract

The article’s objective is to study the magnetohydrodynamic fluctuating free convection flow of incompressible electrically conducting viscoelastic fluid in a porous medium in the presence of a pressure gradient. The flow is between two parallel plates; the lower plate is at rest, while the upper plate moves with a free stream velocityUt. A uniform magnetic field of strength B0 is applied transversely to the fluid motion. The dimensional governing equations and physical initial and boundary conditions are nondimensionalized by inserting appropriate dimensionless variables. Furthermore, Lighthill’s method is used to obtain exact analytical solutions for the velocity, temperature, and concentration distributions. Finally, the influences of inserted parameters are studied through plots and discussed physically. The numerical results are calculated and presented in tabular form for skin friction. It is worth mentioning that when the viscoelastic parameter α1 tends to zero, the obtained general solution is reduced to Newtonian fluid, which shows the validity and correctness of our obtained exact solutions.


---

## Body

## 1. Introduction

The understanding of flow characteristics of non-Newtonian fluid has become a subject of interest due to its widespread applications in industries and engineering, and it has been studied by Chhabra and Richardson [1]. Such interested motivations of a large number of applications have been studied in different fields such as geophysics, biology, chemistry, and petroleum industries [2]. The non-Newtonian fluids have been considered as more appropriate fluids than Newtonian fluids in industrial and technological applications. Because of particular applications of non-Newtonian fluid, it has become an increasingly applied topic of central research. These fluids fall into three classes known as integral, differential, and rate type. The differential sort of non-Newtonian fluids called viscoelastic type fluids has distinct characteristics compared to the viscous materials due to their higher-order governing equations compared to Navier-Stokes equations. Specific examples of non-Newtonian fluids are tomato sauce, polymers, and blood. The analysis of non-Newtonian fluids flow is significant in the field of engineering and oil and gas purification and filtration process. Researchers have reported on viscoelastic fluid flow applications using a porous medium [3–6].MHD has a vital role in the field of engineering and industrial level. The magnetohydrodynamics fluid flow has great interest due to its importance in MHD generator, designing the cooling system with liquid metal, MHD pumps accelerators, and flow meter. Characteristics of magnetic field effects have significant applications in science and engineering, such as MHD pumps, geothermal energy extraction, MHD generators, and many more. A useful review of MHD theory and its applications has been found in Moreau’s book [7]. MHD is the study of the magnetic properties and behavior of electrically conducting fluids. Examples of such magnetofluids include plasmas, liquid metals, salt water, and electrolytes. The word “magnetohydrodynamics” was derived from the word “magneto,” meaning magnetic field, the word “hydro,” meaning water, and the word “dynamics,” meaning movement. The field of MHD was initiated by Hannes Alfvén [8], for which he received the Nobel Prize in Physics in 1970. The fundamental concept behind MHD is that magnetic fields can induce currents in a moving conductive fluid, which in turn polarizes the fluid and reciprocally changes the magnetic field itself. The set of equations that describe MHD are a combination of the Navier-Stokes equations of fluid dynamics and Maxwell’s equations of electromagnetism. These differential equations have been solved simultaneously, either analytically or numerically. MHD applications have been seen in plasma, which is strongly collisional. The collisions’ time scale is shorter than the other characteristic times in the system, and the particle distributions are therefore close to Maxwellian. Interest in length scales much longer than the ion skin depth and Larmor radius perpendicular to the field, long enough along the field to ignore Landau damping, and time scales much longer than the ion gyration time. The importance of significant contributions of MHD in the boundary layer flow has increased among the researchers day by day. Many researchers have examined the MHD flow through a porous medium. The applications of viscoelastic fluid MHD flow have been studied in [9–13]. Makinde and Mhone [11] have discussed the unsteady MHD flow with slip boundary conditions. In another work, Makinde and Osalusi [14] have written a brief discussion on the MHD steady flow with slip and porous boundaries with uniform suction velocity. Makinde and Aziz [15] have investigated the MHD flow considering a porous medium. Barik et al. [16] have debated on the influence of MHD oscillatory flow through short effect. Sivaraj and Kumar [17] have investigated the impact of chemical reactions with a varying concentration on MHD flow. Singh [18] has derived an exact solution for MHD flow in a channel taking porous medium. That work is comprehensively discussed by Mehmood and Ali [19] considering slip condition. The perturbation technique has been used by Kumar et al. [20] to solve the similar problem of Mehmood and Ali [19] who considered the influence of slip boundary conditions. Farhad et al. [21] investigated unsteady magnetohydrodynamic oscillatory flow in a porous channel with heat and mass transfer. Baoku et al. [22] have discussed the effect of thermal radiation on MHD flow in a porous medium. Motsa and Shateyi [23] discussed numerically mixed convection magnetohydrodynamic flow under the effect of ohmic heating. Bakr [24] and Prakash et al. [25] have discussed the impact of unsteady MHD mixed convective flow over a moving porous plate. The nonlinear MHD flow with heat and mass transfer of an incompressible, viscoelastic fluid on a stretching surface and stratification effects have been investigated by Kandasamy et al. [26]. Ishak et al. [27] examined two-dimensional MHD stagnation point flow towards stretching sheet with variable surface temperature. Akbar et al. [28] have analyzed the MHD flow of viscoelastic fluids taking a porous medium. The study of MHD free convection flow of viscoelastic fluid through a porous wedge has been studied by Rashidi et al. [29]. Jang and Lee [30] investigated the micro pump study, where pumping mechanism was based on the MHD principles and electrically conducting fluids in electric and magnetic fields.The sum of all shear stresses of the second order is called second-grade fluids. Non-Newtonian fluids can be further categorized into three types on the shear bases: differential-type fluid, rate-type fluid, and integral-type fluid. The second-grade fluid is a subtype of the differential type of fluids. Erdogan and Irmak [31] have studied the transitory flow of second-grade fluid and reported that pressure is the cause of the fluid’s spontaneous velocity. An analytical solution for the velocity of second-grade fluid has been derived by Fetecau and zierep [32]. Erdogan and Irmak [33] discussed an analytical solution for the governing equation of velocity of a second-grade fluid in three-dimensional flows. Nazar et al. [34] obtained a few closed-form results using the Laplace transform method for oscillatory second-grade fluid flow. This work was further extended by Ali et al. [35] from nonporous hydrodynamic to porous hydromagnetic second-grade fluid flow to obtain closed-form solutions using the same technique. Unsteady MHD flow of second-grade fluid over an oscillating vertical plate in a porous medium with heat and mass transfer has been investigated by Farhad et al. [36]. The homotopy analysis method (HAM) has been used by Hayat et al. [37] for viscoelastic fluid flow to the derived analytical solution for velocity and temperature profile. Anwar et al. [38] have discussed the laminar flow of second grade with a horizontal cylinder. The steady magnetohydrodynamic flow of viscoelastic fluid is investigated by Eldabe et al. [39]. Mustafa [40] has considered a transverse magnetic field through a porous medium of second-grade fluid in a tube. Shah and Khan [41] have discussed the convective flow of second-grade fluid to obtain exact solutions for heat and velocity by Laplace transforms technique.The velocity of the fluid above the boundary layer or free stream velocity is defined as the speed at which the air is infinitely far from the aircraft moving when the object is kept as referential called free stream velocity. Free stream velocity is always greater than the boundary layer velocity because free stream velocity starts from the upper boundary of the fluid’s boundary layer velocity. In (1954), Lighthill has studied the laminar boundary layer, time-dependent free stream oscillation. An unsteady mixed convection flow over a moving vertical plate in a vertical free stream with combined effects of boundary force and thermal diffusion has been studied by Patil et al. [42]. A uniform free stream with a fixed or moving surface has been investigated by Ishak et al. [43]. The nonsimilar boundary layer flow of viscous fluid over a moving surface in a rotating fluid with free stream velocity has been explained by Takhar et al. [44]. The effect of horizontal oscillations of the free stream velocity on pitching airfoil oscillations using computational fluid dynamics has been studied by Gharali and Johnson [45]. Alassar et al. [46] investigated the viscous flow over a sphere with fluctuations in the free stream velocity. The effect of large temperature with the determination of skin friction for a laminar boundary layer flow with variable free stream velocity has been examined by Levy et al. [47]. Kelly [48] has studied the steady viscous flow near a stagnation point following a change free stream velocity with a finite instant of time in two dimensions. Mei and Adrian [49] have investigated unsteady flow over a stationary with small fluctuation in the free stream velocity at a finite Reynolds number. Rott et al. [50] have presented a complete solution of the laminar boundary layer flow of small fluctuations of the free stream velocity.Bearing in mind the above literature review, we found a gap that no one has considered the free stream fluctuating velocity on the boundary in their studies. Therefore, to fill this gap, we considered in the present work free convective, unsteady MHD flow of incompressible, electrically conducting, second-grade fluid flow, flowing through a porous medium. Furthermore, free stream velocity and pressure gradient are taken into account. Lighthill’s [51] method has been applied to obtain analytical solutions for the velocity, temperature, and concentration distributions. The effect of various embedded parameters on the skin friction is shown in tabular form. The graphs are sketched on the velocity, temperature, and concentration profiles to show the effect of several embedded parameters.


## 2. Mathematical Formulation

An incompressible, viscoelastic, second-grade fluid flow through a porous medium has been considered in the present work. The fluid is taken to be electrically conducting, and the magnetic field is transversely applied to the flow. As the magnetic Reynolds number is smaller, we are neglecting the induced magnetic field’s influence. The fluid velocity in the boundary layerU is the free stream velocity away from the boundary layer. We assume that the medium’s lower boundary is at constant temperature T0 and concentration, C0, while the ambient fluid has temperature Tw and concentration level Cw uniformly. Geometry of the flow is shown in Figure 1.Figure 1
Geometry of the flow.The fluid flow is along thex-axis, and the y-axis is normal to the flow direction. Three main equations govern the given flow regime: one is the momentum equation, the second is the heat equation, and the third one is concentration distribution, such as [21](1)ρ∂uy,t∂t=−∂p∂x+μ∂2uy,t∂y2+α1∂∂t∂2uy,t∂y2+gβTTy,t−Ty,t0+gβCCy,t−C0y,t−σB0uy,t−ϕ1k11+α1∂μ∂tuy,t,(2)ρCp∂Ty,t∂t=k∂2Ty,t∂y2−∂qr∂y,(3)∂Cy,t∂t=D∂2Cy,t∂y2+DKtTm∂2Ty,t∂y2.Initial and boundary conditions are given as(4)u0,t=0u1,t=1+εeiωtT0,t=T0T1,t=TwC0,t=C0C1,t=Cw.The velocity, heat, and concentration equations show density, dynamic viscosity, second-grade parameter, electrical conductivity, thermal expansion coefficient, concentration expansion coefficient, gravitational acceleration and heat capacity, porosity, and permeability parameters, respectively.D shows the mass diffusivity, k represents thermal conductivity of the fluid, qr is the radiative heat flux, “a” shows the width of the channel, KT is the thermal diffusion ratio, and ω is the frequency of oscillation. The temperature of the plates is TwandT0 which produces radiation heat transfer. The radiation flux of heat is assumed as(5)−∂qr∂y=4α02T−T0,where α0 is the radiation absorption coefficient.For dimensional analysis, the dimensionless variables are(6)u∗=uUo,y∗=ya,t∗=U0at,θ=T−T0Tw−T0,ϕ=C−C0Cw−C0.For the free stream velocity, we convert the velocity of the fluiduy,t with free stream velocity Ut all functions being considered as zero except Ut and, thus, we get the following equation:(7)−1ρ∂p∂x=∂U∂t+σB02U+μϕ1ρk11+α1∂μ∂tU.Putting equation (7) in equation (1), we get(8)∂uy,t∂t=∂Ut∂t+σρB02Ut+μϕ1ρk11+α1∂μ∂tUt+μ∂2uy,tρ∂y2+α1∂ρ∂t∂2uy,t∂y2+gρβTT−T0+gρβCC−C0−σρB0uy,t−ϕ1ρk11+α1∂μ∂tuy,t.After rearranging the terms, we have(9)∂uy,t∂t=∂Ut∂t+μ∂2uy,tρ∂y2+α1∂ρ∂t∂2uy,t∂y2+gρβTT−T0+gρβCC−C0−σρB0uy,t+σρB02Ut−ϕ1ρk11+α1∂μ∂tuy,t+μϕ1ρk11+α1∂μ∂tUt.Taking some standard terms and then arranging, we get(10)∂uy,t∂t=∂Ut∂t+μ∂2uy,tρ∂y2+α1∂ρ∂t∂2uy,t∂y2+gρβTT−T0+gρβCC−C0−σρB02uy,t−Ut−μϕ1ρk11+α1∂μ∂tuy,t−Ut,where u shows boundary layer velocity and U is used for the free stream velocity.After nondimensionalization, equations (2), (3), and (10) become(11)Pe∂θ∂t=∂2θ∂y2+N2θ,(12)1Sc∂ϕ∂t=∂2ϕ∂y2+Sr∂2θ∂y2,(13)a1∂u∂t−a1dUdy−∂2u∂y2−α∂∂t∂2u∂y2+Hu−HU−Grθ−Gmϕ=0.Initial and boundary conditions are(14)u0,t=0u1,t=1+εeiωtt>0,θ0,t=1θ1,t=0t>0,ϕ0,t=1ϕ1,t=0t>0.After nondimensionalization, we get some nondimensional parameters, such as(15)Re=ρaU0μ,α=α1U0μa,M=σB02a2μ,1k=a2k1,Gr=gβta2Tw−T0μU0,Gm=gβca2Cw−C0μU0,Sc=DaU0,Sr=KtTw−T0TmCw−C0,Pe=ρCpaU0k,a1=Re+αk,H=M+1k,N2=4α02a2k,where dimensionless viscoelastic parameter α, the magnetic and permeability parameters M, K, and Re,N,Pe,Gr,Gm,Sc, and Sr are Reynolds number, radiation, Peclet, thermal Grashof, and mass Grashof parameters, as well as Schmidt and Soret numbers, respectively.


## 3. Solution of the Model

We are solving equations (11)–(13) by Lighthill’s method, and equations (16)–(18) are the assumed solutions:(16)uy,t=u0y+u1yεeiωt+Oε2,(17)θy,t=θ0y+θ1yεeiωt+Oε2,(18)ϕy,t=ϕ0y+ϕ1yεeiωt+Oε2.By using the assumed solution and solving equations (11)–(13), we get(19)Pe∂θ∂t=∂2θ∂y2+N2θ.The nonharmonic part and boundary conditions for the heat equation are(20)∂2θ0y∂y2+N2θ0y=0,θ00=1,θ0a=0.Solution of the harmonic part is(21)θ0y=cosNy−a2sinNy.The harmonic part and conditions for the heat equation are(22)∂2θ1y∂y2+N2−iωPeθ1y=0,θ10=0,θ1a=1.Solution of the nonharmonic part is(23)θ1y=0.Putting the values of equations (21) and (23) in the assumed solution (17), we get(24)θy,t=θ0y+θ1yεeiωt+Oε2,(25)θy=cosNy−a2sinNy.Now, we will solve the concentration equation:(26)1Sc∂ϕ∂t=∂2ϕ∂y2+Sr∂2θ∂y2.The nonharmonic part and boundary conditions for concentration equation are(27)∂2ϕ0y∂y2+Sr∂2θ0y∂y2=0,ϕ00=1,ϕ0a=0.Solution of the nonharmonic part is(28)ϕ0y=a4+a3y−SrcosNy−a2sinNy.The harmonic part and conditions for concentration equation are(29)∂2ϕ1y∂y2−iωScϕ1y=−Sr∂2ϕ1y∂y2,ϕ10=0,ϕ1a=0.Solution of the harmonic part is(30)ϕ1y=0.Putting equations (28) and (30) in the assumed solution (18), we get(31)ϕy,t=ϕ0y+ϕ1yεeiωt+Oε2,(32)ϕy=a4+a3y−SrcosNy−a2sinNy.Now solve the momentum equation:(33)a1∂u∂t−a1dUdy−∂2u∂y2−α∂∂t∂2u∂y2+Hu−HU−Grθ−Gmϕ=0.Putting the values of equations (25) and (32) in (33), we have(34)a1∂u∂t−a1dUdy−∂2u∂y2−α∂∂t∂2u∂y2+Hu−HU−GrcosNy−a2sinNy−Gma4+a3y−SrcosNy−a2sinNy=0.The nonharmonic part and its conditions are(35)∂2u0y∂y2−Hu0y=−H+Gma4−Gr−SrGmcosNy−a2sinNy−a3Gmy,u00=0,u0a=1.The solution of nonharmonic part is(36)u0y=C1coshHy+C2sinhHy+H+Gma4+a3yH+Gr−SrGmH+N2cosNy−a2sinNy.The harmonic part and its conditions are(37)∂2u1y∂y2−m1u1y=−m1,u10=0,u1a=1.Solution of the harmonic part of the momentum equation is(38)u1y=−coshm1y+a5sinhm1y+1.Assume solution for the momentum equation:(39)uy,t=u0y+u1yεeiωt+Oε2,(40)uy,t=C1coshHy+C2sinhHy+H+Gma4+a3yH+Gr−SrGmH+N2cosNy−a2sinNy+εeiωt−coshm1y+a5sinhm1y+1,where m1 is a constant:(41)m1=H+iωa1+iωα.We have solved heat, concentration, and momentum equations and obtained the following results:(42)θ1y=0,ϕ1y=0.(43)θ0y=cosNy−a2sinNy,(44)ϕ0y=a4+a3y−SrcosNy−a2sinNy.Using equations (43) and (44) in equation (40), we get(45)a2=cosNsinN,a3=SrcosN−a2sinN−1,a4=1+Sr,C1=−H+Gma4+a3H−Gr−SrGmH+N2,C2=1sinhH1−C1coshH−H+Gma4+a3H−Gr−SrGmH+N2cosN−a2sinN,a5=coshm1sinhm1.Now putting the valuesθ0y,θ1yandϕ0y,ϕ1y in equations (26) and (31), we can obtain the following results:(46)θy,t=cosNy−a2sinNy,aty=0,ϕy,t=a4+a3y−SrcosNy−a2sinNy,uy,t=C1coshHy+C2sinhHy+H+Gma4+a3yH+Gr−SrGmH+N2cosNy−a2sinNy+εeiωt−coshm1y+a5sinhm1y+1.Furthermore, it has been noticed from equation (32) that the solution can be reduced into a Newtonian fluid by taking α=0, which shows the validity of the problem.


## 4. Graphical Results and Discussion

The exact solution for free convection flow of second-grade fluid has been investigated. The governing partial differential equations and corresponding boundary conditions have been converted into ordinary differential equations and are analytically solved. Figures2–9 are plotted to influence different physical parameters on the velocity of the fluid. Figure 10 is drawn for the parametric influence of radiation parameters for temperature profile. Figures 11 and 12 depict the effect of radiation parameter and Soret number in concentration profile. Table 1 illustrates a variation in skin friction for different parameters.Figure 2
Velocity profile for the various values ofM when Re = 0.2, Gr = 1.1, K = 0.6, Gm = 1.2, N = 2.2, Sr = 1, α = 0.71, ε = 0.001, and ω = π/2.Figure 3
Velocity graph for the several values ofGr when Re = 0.2, M = 1.1, K = 0.5, Gm = 1.1, N = 2, Sr = 1, α = 0.71, ε = 0.001, and ω = π/2.Figure 4
Velocity profile used for the several values of Gm whenRe = 5, Gr = 1.1, M = 2.1, K = 0.5, N = 2, Sr = 1, α = 0.71, ε = 0.001, and ω = π/2.Figure 5
Velocity graph for various values ofK when Re = 0.2, Gr = 2, M = 1.1, Gm = 1.1, N = 2, Sr = 1, α = 0.71, ε = 0.001, and ω = π/2.Figure 6
Velocity profile of several values forα when Re = 5, Gr = 1.1, M = 2.1, K = 0.52, Gm = 1.1, N = 0.001, Sr = 1.1,ε = 2, and ω = π/2.Figure 7
Velocity for the several values ofSr when Re = 2.2, Gr = 2.2, M = 1.1, K = 1, Gm = 1, N = 2, α = 7, ε = 0.001, and ω = π/2.Figure 8
Velocity graph for the several values ofRe when Gr = 1.1, M = 2.2, K = 0.51, Gm = 1, N = 2, Sr = 1, α = 0.71, ε = 0.001, and ω = π/2.Figure 9
Velocity profile for the several values ofNwhen Re = 0.2, Gr = 5, M = 2, K = 5, Gm = 0.1, Sr = 0.4, α = 0.3, ε = 0.001, and ω = π/2.Figure 10
Effect of temperature equation for the several values ofN when t = 1.01.Figure 11
The effect of concentration equation for the various values ofN when Sr = 0.71.Figure 12
Effect of concentration equation for the various values ofSr when N = 1.1.Table 1
The influence of different parameters on (Cf) skin friction.
MGrGmKNSrαεReωCfM = 110.51.50.02220.20.110.80.5233M = 220.51.50.02220.20.110.80.5314M = 330.51.50.02220.20.110.80.5394Gr = 0.50.70.51.50.02220.20.110.80.5208Gr = 0.70.70.71.50.02220.20.110.80.5216Gr = 0.90.70.91.50.02220.20.110.80.5221Gm = 0.50.70.50.50.02220.20.110.80.5800Gm = 0.60.70.50.60.02220.20.110.80.5740Gm = 0.70.70.50.70.02220.20.110.80.5681K = 0.010.70.51.50.01220.20.110.80.8157K = 0.020.70.51.50.02220.20.110.80.5208K = 0.030.70.51.50.03220.20.110.80.3664N = 20.70.51.50.02220.20.110.80.5208N = 40.70.51.50.02420.20.110.80.5204N = 60.70.51.50.02620.20.110.80.5200Sr = 10.70.51.50.02210.20.110.80.5504Sr = 20.70.51.50.02220.20.110.80.5208Sr = 30.70.51.50.02230.20.110.80.4912α = 0.20.70.51.50.02220.20.110.80.5208α = 0.40.70.51.50.02220.40.110.80.4769α = 0.60.70.51.50.02220.60.110.80.4327ε = 0.20.70.51.50.02220.20.210.81.1304ε = 0.40.70.51.50.02220.20.410.82.3496ε = 0.60.70.51.50.02220.20.610.83.5688Re = 10.70.51.50.02220.20.110.80.5208Re = 20.70.51.50.02220.20.120.80.5187Re = 30.70.51.50.02220.20.130.80.5166ω = 0.10.70.51.50.02220.20.110.10.6216ω = 0.30.70.51.50.02220.20.110.30.6086ω = 0.50.70.51.50.02220.20.110.50.5827Figure2 shows the influence of the magnetic parameter on the velocity profile; by increasing M, the fluid’s velocity under the boundary layer decreases, which indicates the production of Lorentz force. In our case, the velocity of the fluid increases due to free steam velocity. Free stream velocity is always greater than the boundary layer velocity because the fluid’s velocity increases. Physically it means that under the boundary layer velocity of the fluid decreases by increasing the magnetic parameter M, which causes the drag forces or frictional force, resisting the flow and reducing the velocity. The influence ofMillustrates that the greater value of M increases the velocity.Figure3 represents the velocity profile for the thermal Grashof number, increasing the velocity because of the enhancement of buoyancy force, which indicates that the thermal buoyancy force accelerates the velocity. Thermal Grashof number tells about the velocity effects of thermal buoyancy force to viscous force. It is observed that an increase in Gr leads to increase velocity due to enhancement in buoyancy force. Physically, Gr signifies the relative effect of the thermal buoyancy force to the viscous hydrodynamic force. Increase of Grashof number Gr means an increase of temperature gradients due to which the contribution from the buoyancy near the plate becomes significant and hence a short rise in the velocity near the plate is observed. Gr > 0 means heating of the fluid or cooling of the plate by natural convection. For the positive value of Gr, heat is conducted away from the boundary into the fluid which increases the temperature and thereby enhances the buoyancy force.Figure4 is sketched to demonstrate the behavior of mass Grashof number on the velocity profile. It is investigating that Gm = 0 corresponds to the buoyancy force and Gm > 0 corresponds to the strong magnetic field chilling problem faced in atomic engineering connected with the chilling of the nuclear reactor. It is to be noted that a more significant number of buoyancy forces create and thus speed up the fluid flow. Boundary layer thickness becomes large along “a,” where “a” is the boundary-level wideness and velocity increases Gm.Figure5 illustrates the permeability parameter’s impact, which shows that increasing permeability produces the drag force or resistance, decreasing the velocity. The fact is that by increasing the amount of K the porous medium increases, and the opposing force becomes the cause of decreasing the velocity of the flow. Physically, it has been recorded that increasing the value of K increases the resistance force which decreases the velocity.Figure6 illustrates the influence of the viscoelastic parameter α on the velocity profile, which indicates that when α increases, velocity also increases. The area of the boundary layer increases in the center point, where α increases. The viscoelastic parameter α near the boundary y = 0 is weaker and increases far away from the boundary. At α = 0, physically, the second-grade fluid has a faster velocity than Newtonian fluid because of the elastic properties of a second-grade fluid, which are not found in Newtonian fluids.The influence of the Soret number is depicted in Figure7 for the various values. It has been examined that the velocity of the fluid reduces by increasing the Soret number. The Soret effect is the ratio between temperature and concentration distribution, which causes decrease velocity of the fluid. The influence of Sr on the velocity shows that the magnetic parameter is stronger than the Soret number.Figure8 demonstrates the influence of Reynolds number Re on the velocity profile. It is noticed from the velocity profile that increasing the value of Re causes the decrease of velocity. Physically it means that Reynolds number increases the fluid density due to the fact that velocity of the fluid decreases. Reynolds number directly relates to the density; hence, by increasing R, it retards velocity of the fluid.Figure9 shows the influence of the radiation parameter N on the velocity profile. It is shown in Figure 9 that the velocity of the fluid increases by increasing radiation parameter N, which consequently enhances thermal conductivity k⟶∞, which retards the velocity of the fluid. It means that velocity reduces due to high thermal radiation.Figures10 and 11 depict the effect of various values of heat θ and concentration distribution ϕ. Figure 10 shows that the greater the amount of thermal radiationN is, the more the range of the temperature rises; it means that the thermal boundary layer thickness increases with increasingN. Still, the concentration has opposite influence to the heat profile, which is shown in Figure 11.Soret number’s influence is shown in Figure12, which demonstrates that it is linked to the mass transfer of the fluid flow, and Soret number is the ratio between temperature and concentration. The concentration profile shows that, by increasing Soret number, a considerable variation has been seen on the concentration profile.The influences of several parameters are discussed in Table1. It has been noticed that by increasing magnetic parameter M that the skin friction rises. The skin friction reduces when the amount of mass Grashof number Gm rises, and permeability parameter K, Soret number Sr, viscoelastic parameter, Reynolds number, and phase angle reduce skin friction, but M has an opposite effect on the skin, increasing skin friction. The influences of M and permeability parameters K are discussed.


## 5. Conclusion

In the current work, unsteady MHD fluctuating free convection flow of second-grade fluid through a porous medium in the pressure gradient presence is studied. The effects of different parameters have been shown graphically, and skin friction is evaluated in tabular form. The perturbation technique is used to obtain exact analytical solutions for velocity, heat, and mass equations. The graphical results are discussed in detail. The primary outcomes of the present study are summarized as follows:(1)
When the viscoelastic parameterα and mass Grashof Gm increase, the fluid’s velocity increases, while skin friction decreases by increasing both these parameters.(2)
Velocity and skin friction decrease for the large amounts ofRe, Sr, and K.(3)
When the value of magnetic parameterM increased, velocity and skin friction increased.(4)
The effect of phase angle on velocity and skin friction shows that by increasingω fluid velocity increases while the skin friction decreases.(5)
IncreasingGr the velocity of the fluid flow also increases the fluid velocity by increasing the radiation parameter N.



---

*Source: 6648281-2021-04-07.xml*