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 velocity . 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 tends to zero, the obtained general solution is reduced to Newtonian fluid, which shows the validity and correctness of our obtained exact solutions.
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 layer is the free stream velocity away from the boundary layer. We assume that the medium’s lower boundary is at constant temperature and concentration, , while the ambient fluid has temperature and concentration level uniformly. Geometry of the flow is shown in Figure 1.
The fluid flow is along the x-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]
Initial and boundary conditions are given as
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, represents thermal conductivity of the fluid, is the radiative heat flux, “a” shows the width of the channel, is the thermal diffusion ratio, and is the frequency of oscillation. The temperature of the plates is which produces radiation heat transfer. The radiation flux of heat is assumed aswhere is the radiation absorption coefficient.
For dimensional analysis, the dimensionless variables are
For the free stream velocity, we convert the velocity of the fluid with free stream velocity all functions being considered as zero except and, thus, we get the following equation:
Putting equation (7) in equation (1), we get
After rearranging the terms, we have
Taking some standard terms and then arranging, we getwhere shows boundary layer velocity and is used for the free stream velocity.
After nondimensionalization, equations (2), (3), and (10) become
Initial and boundary conditions are
After nondimensionalization, we get some nondimensional parameters, such aswhere dimensionless viscoelastic parameter , the magnetic and permeability parameters , , and , and 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:
By using the assumed solution and solving equations (11)–(13), we get
The nonharmonic part and boundary conditions for the heat equation are
Solution of the harmonic part is
The harmonic part and conditions for the heat equation are
Solution of the nonharmonic part is
Putting the values of equations (21) and (23) in the assumed solution (17), we get
Now, we will solve the concentration equation:
The nonharmonic part and boundary conditions for concentration equation are
Solution of the nonharmonic part is
The harmonic part and conditions for concentration equation are
Solution of the harmonic part is
Putting equations (28) and (30) in the assumed solution (18), we get
Now solve the momentum equation:
Putting the values of equations (25) and (32) in (33), we have
The nonharmonic part and its conditions are
The solution of nonharmonic part is
The harmonic part and its conditions are
Solution of the harmonic part of the momentum equation is
Assume solution for the momentum equation:where m1 is a constant:
We have solved heat, concentration, and momentum equations and obtained the following results:
Using equations (43) and (44) in equation (40), we get
Now putting the values in equations (26) and (31), we can obtain the following results:
Furthermore, it has been noticed from equation (32) that the solution can be reduced into a Newtonian fluid by taking , 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. Figures 2–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 shows the influence of the magnetic parameter on the velocity profile; by increasing , 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 , which causes the drag forces or frictional force, resisting the flow and reducing the velocity. The influence ofillustrates that the greater value of increases the velocity.
Figure 3 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.
Figure 4 is sketched to demonstrate the behavior of mass Grashof number on the velocity profile. It is investigating that = 0 corresponds to the buoyancy force and > 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 .
Figure 5 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 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.
Figure 6 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 Figure 7 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 on the velocity shows that the magnetic parameter is stronger than the Soret number.
Figure 8 demonstrates the influence of Reynolds number on the velocity profile. It is noticed from the velocity profile that increasing the value of 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 , it retards velocity of the fluid.
Figure 9 shows the influence of the radiation parameter on the velocity profile. It is shown in Figure 9 that the velocity of the fluid increases by increasing radiation parameter , which consequently enhances thermal conductivity , which retards the velocity of the fluid. It means that velocity reduces due to high thermal radiation.
Figures 10 and 11 depict the effect of various values of heat and concentration distribution . Figure 10 shows that the greater the amount of thermal radiation is, the more the range of the temperature rises; it means that the thermal boundary layer thickness increases with increasing. Still, the concentration has opposite influence to the heat profile, which is shown in Figure 11.
Soret number’s influence is shown in Figure 12, 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 Table 1. It has been noticed that by increasing magnetic parameter that the skin friction rises. The skin friction reduces when the amount of mass Grashof number rises, and permeability parameter , Soret number , 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 and permeability parameters 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 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 of , , and .(3)When the value of magnetic parameter 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)Increasing the velocity of the fluid flow also increases the fluid velocity by increasing the radiation parameter .
Abbreviations
| : | Boundary layer velocity of the fluid |
| : | The velocity of the dust particle |
| : | Free stream velocity |
| : | The temperature of the fluid |
| : | The temperature of the particle |
| : | Fluid density |
| : | kinematic viscosity |
| : | Second-grade parameter |
| : | Stock’s resistance coefficient |
| : | The number density of the dust particle, which is assumed to be constant |
| : | Electrical conductivity |
| : | Applied magnetic field |
| : | Gravitational acceleration |
| : | Coefficient of thermal expansion |
| : | Coefficient of volumetric expansion |
| : | Thermal conductivity of the base fluid |
| : | Specific heat capacity |
| : | Radiative heat flux |
| : | Frequency of oscillation |
| : | Radiation absorption coefficient |
| : | Schmidt number |
| : | The density of the dust particle |
| : | Temperature relaxation time |
| : | Temperature of wall |
| : | Ambient temperature |
| : | Specific heat capacity of the dust particle |
| : | The average mass of the dust particle |
| : | Nondimensional parameter |
| : | Peclet number |
| Gr: | Grashof number |
| : | Reynolds number |
| : | Particle concentration parameter |
| : | Heat absorption coefficient |
| : | Dusty fluid parameter |
| : | Dusty fluid parameter |
| : | Nondimensional temperature relaxation time parameter |
| : | Nondimensional second-grade parameter |
| : | Mass diffusivity |
| : | Width channel |
| : | Thermal diffusion ratio |
| : | Soret number. |
Data Availability
All the data are available in the manuscript.
Conflicts of Interest
The authors declare that there are no conflicts of interest.