- 1Department of Mathematics, University of the Punjab, Lahore, Pakistan
- 2Mathematics Department, Umm Al-Qura University, Makkah, Saudi Arabia
- 3Department of Mechanical Engineering, College of Engineering, King Khalid University, Abha, Saudi Arabia
- 4Faculty of Engineering and Technology, Future University in Egypt New Cairo, New Cairo, Egypt
- 5Department of Industrial Engineering, College of Engineering, King Khalid University, Abha, Saudi Arabia
The current research article discusses the two-dimensional, laminar, steady, and incompressible third-grade viscoelastic micropolar fluid flow along with thermal radiation caused by an exponentially stretched sheet. The primary goal of this extensive study is to improve thermal transportation. Thermophoresis and Brownian motion are two key causes of nanoparticle migration in nanofluids, and their impacts on the thermophysical properties of nanofluids are significant. Micropolar fluids are investigated due to their micro-motions that are significant in convective thermal and mass transport polymer formation, nanotechnology, and electronics. The consequences of third-grade fluid parameters, thermophoresis and Brownian motion, induced magnetic field, micro-polarity, and micro-inertia density on the stream of an electrically conductive fluid are analyzed. A homogeneous magnetic field is supplied perpendicularly to the surface, and the liquid is believed to be electrically conducting. As the flow has a significant magnetic Reynolds number, the contribution of the evoked magnetic field is properly accounted in the governing equations. A mathematical model in the form of partial differential equations (PDEs) is built under certain assumptions. By invoking the suitable similarity transformation, the non-linear PDEs are modified into dimensionless coupled ordinary differential equations (ODEs). The MATLAB numerical technique bvp4c is employed to settle the subsequent ODEs together with the boundary constraints. The consequences of numerous physical parameters on the non-dimensional concentration, temperature, micropolar, velocity, and induced magnetic field profiles are portrayed in graphs. It is found that the concentration boundary layer, thermal boundary layer, and micropolar boundary layer thickness decelerate with the increment in the micro-polarity of the fluid.
1 Introduction
Magnetohydrodynamics (MHD) is a branch of fluid mechanics that studies the movement of an electrically conductive liquid in the existence of the magnetic field. Alfven [1] was the first who found that the movement of conducting fluid between magnetic field lines generates potential differences, that, in turn induces the flow of electric currents. The magnetic fields coupled with these electric currents alter the magnetic force that generates them. Alternatively, the fluid stream modifies the system’s electromagnetic structure. The propagation of electric current through a magnetic field, on the other hand, is coupled with a body force, known as the Lorentz force, which affects fluid motion. Scientists and researchers in the field of fluid dynamics have identified its use in metal dispersion, mining, fusion reactors, targeted drug delivery, MHD-based laser beam scanning, construction of MHD pumps, MHD generators, and MHD flow meters. Andersson [2] probed the MHD stream of a viscoelastic liquid past over an expanding sheet and showed that external magnetic field affects the flow similar as viscoelasticity. Hayat and Abbas [3] investigated the radiation impacts on the MHD stream of a viscoelastic liquid across a permeable surface. Nadeem and Hussain [4] explored the shrinking solutions in the existence of MHD stream of a viscous liquid toward a non-linear contracting surface. The effects of an exponentially stretched surface on MHD boundary layer stream and heat transmission characteristics in a permeable medium have been investigated by Ahmad et al. [5]. Hayat et al. [6] discussed the MHD boundary layer stream of a viscous fluid due to non-linear extending cylinder with thermal slip stratification and radiation. Sohail et al. [7] probed the boundary layer stream of the steady MHD Carreau liquid with bioconvection across the heated disk. They determined that the intensity of the magnetic field reduces fluid’s particle velocity, although mass transfer and fluid temperature increase with increasing magnetic field values. Riaz et al. [8] probed the MHD and entropy formation modeling of nanoliquid with three-dimensional (3-D) peristaltic cylindrical confinement. Employing molecular dynamical modeling, AbdulHussein et al. [9] investigated the boiling mechanism of several fluids in micro-channels in the existence of an exterior electromagnetic field.
The study of fluid flow across an expanding sheet is crucial for a variety of applications including, extrusion, glass blowing, cord depiction, copper spiraling, strengthening and tinning of copper wires, glass blowing, warm progressing, thermal conductivity of heat sinks, and melts of high molecular weight polymers. Hussain et al. [10] explored the influences on the boundary layer stream of a micropolar liquid flowing on the stretched surface. Abbas et al. [11] inspected the micropolar fluid flow behavior across the stretched sheet due to a variety of significant applications. Awan et al. have investigated different types of flows over the stretching surface and porous media, with various physical implications [12, 13]. Nadeem and Khan [14] studied the rotating Maxwell nanofluid flow between linear and exponential expanding sheets. Riaz et al. [15] exposed the bioconvection mechanism in the stream of magnetically polarized Williamson nanoparticles along with activating energy and heat source/sink.
Many scientists from all over the world are eager to learn more about non-Newtonian fluid flow. The reason for such a motivation in the study of these fluids is as a result of their use in industries and technologies including suspension fabrication, detergent and paint production, skincare creams, polymer production, spinning of metal etc. The non-Newtonian fluids in comparison to the Newtonian fluids are rheological in structure having non-linear correlation between shearing stress and velocity gradient. Non-Newtonian liquids are classified into three kinds depending upon their characteristics: rate type, differential type, and integral type fluids. The differential type liquid models have found to be well known among them. The second-grade fluid is the most basic subclass of these viscoelastic models, capturing typical stress variations but not anticipating shear thinning/thickening processes. In contrast, the third-grade liquid model can predict both ordinary stress and shear thinning/expanding processes. Abbasbandy et al. [16] investigated the both exact and series approaches for third-grade fluid by using thin film. Hayat et al. [17] probed the rotating stream of third-grade liquid among two permeable sheets by applying MHD effects. Hayat et al. [18] probed the MHD nanofluid stream of second-grade fluid caused by a non-linear stretched surface. In the vicinity of nanoparticles, the non-transient motions of a third-grade fluid driven by a pressure sort die are investigated by Mahanthesh and Joseph [19]. The characteristics of an applied magnetism and entropy formation on Jeffrey nanoliquid in between an annulus region of two small non-concentric pipes was disclosed by Riaz et al. [20]. Mondal et al. [21] discussed nanoliquid stream across a permeable vertical plate along with interior heat production and non-linear thermal radiations. Sangeetha and De [22] analyzed the bioconvection in nanoliquid stream with viscous losses and Ohmic heating. Riaz et al. [23] presented a comparative analysis of entropy assessment on a (3-D) wavy stream of Eyring–Powell nanoliquid.
Micropolar fluids have a microstructure and are coupled to fluids with a non-symmetrical stress tensor. Because of the liquid particles’ localized composition and micro-motions, micropolar liquids display particular microscopic features. They characterize fluids composed of stiff, arbitrarily oriented, or cylindrical particles dispersed in a viscous material, where fluid’s particle distortion is neglected. Eringen [24] was the first to investigate the hypothesis of micropolar fluids. In this hypothesis, the continuum is defined as the set of systematic particles that have not only momentum but also a sub-structure. In other words, each material volume element is made up of micro volume components that can translate and twist independently of each other. Gorla et al. [25] scrutinized the impacts of buoyancy on driven convection in an axially symmetric stagnation stream of micropolar fluids across a vertically placed cylinder. Rehman and Sattar [26] explained the MHD convective stream of a micropolar liquid over the consistently moving permeable sheet. The temperature boundary layer stream caused by a linearly stretched surface submerged in a constant density micropolar liquid together with radiation effects is probed by [27]. Gaffar et al. [28] examined the free convective boundary layer stream of viscoelastic third-grade micropolar liquid passing over a vertically positioned isothermal cone. Ali et al. [29] revealed the significance of MHD on the micropolar nanoliquid stream across an expanding surface along with radiation and heat stratification influences. Jiang et al. [30] conducted the numerical assessment of the passive usage of phase-change elements and the active usage of nanoliquid inside a rectangular channel.
To the extent of the writer’s insights, no research has been conducted to investigate the magnetohydrodynamics (MHD), and heat transmission effects of a viscous, micropolar, and third-grade fluid passed across an exponentially stretched sheet. The current research is presented to analyze the consequences of third-grade fluid parameters, thermophoresis and Brownian motion, induced magnetic field, micro-polarity, and micro-inertia density on the stream of an electrically conductive fluid. The primary purpose of this extensive study is to improve thermal transportation subject to the existence of micro-rotations of tiny nanoparticles. The mathematical model is built under considered flow assumptions and Buongiorno model. The similarity transformation is a technique that is commonly used for solving various flow problems, in which the system of partial differential equations (PDEs) is modified into ordinary differential equations (ODEs) to solve them analytically or numerically. The current model is simplified by applying similarity transformation. The resultant ODEs are numerically tackled in MATLAB using the bvp4c algorithm, and the numerical tables are constructed to ensure the validity of results. The obtained findings have significant engineering and technological applications.
2 Mathematical Analysis and Flow Geometry
2.1 Constitutive Model for Third-Grade Fluid
For an incompressible fluid, the equations of motion and continuity are as follows:
Here, ρ, V, T, b, J, and B are density of the fluid, velocity vector, Cauchy stress tensor, body force, electric current, and external magnetic field, respectively. Following Rajagopal and Fosdick [31], the third-grade fluid’s stress tensor is as follows:
Here, p, I, T, and S are the pressure, identity tensor, Cauchy stress tensor, and the extra stress tensor, respectively. Furthermore,
Here,
In the case of third-grade fluid, the material moduli satisfy Clausius–Duhem inequality stated as:
So, the stress tensor takes the following form.
2.2 Constitutive Model for Micropolar Fluid
The field equations for the micropolar fluid following Papautsky et al. [32] are stated as follows:
Eqs. 8–10, respectively, denote conservation of mass, linear momentum, and angular momentum. Furthermore, ρ represents the micropolar fluid density, G is the (gyration) micro-rotation vector, V depicts velocity vector, j represents micro-inertia, p is the pressure, l denotes body couple per unit mass vector , λ∗ reflects Eringen second-degree viscosity parameter, f depicts the body strength per unit mass vector, κ represents coefficient of vortex viscosity, μ denotes dynamic viscosity, and
2.3 Problem Formulation
In this section, the steady, laminar, two-dimensional, and incompressible third-grade viscoelastic micropolar fluid stream along with thermal radiation caused by an exponentially stretched surface is reported. The x-axis is marked parallel to the sheet, while the y-axis is directed orthogonal to the sheet. The uniform magnetic field is supplied to the surface in a normal direction, and the fluid is believed to be electrically conductive. The sheet and the third-grade micropolar liquid are both originally maintained at the identical temperature. The temperature is elevated to Tf, and Tf > T∞ where T∞ is the ambient temperature, which remains constant as (see Figure 1). The concentration at the sheet is taken as Cf, while the ambient concentration is C∞ (Cf > C∞). This is supposed that the flow’s magnetic Reynolds number is not low in magnitude, and thus the produced magnetic field is not negligible. The evoked magnetic field H2 is also assumed to be supplied along the x-axis. The parallel component H1 of the induced magnetic field tends to He(x) in the free stream flow, and the slope of the magnetic field approaches zero along y-axis.
Using Eq. 7 in Eq. of motion (Eq. 2) along with the boundary layer estimations [33–36] in case of third-grade fluid, notably, inside the boundary layer
The concerned boundary conditions are as follows:
Here, u, v are the velocity coefficients in x and y directions, ρ depicts the liquid’s density , k denotes the vertex viscosity, μ represents dynamic viscosity, N denotes the micropolar fluid’s angular velocity, T represents the temperature, U0 depicts the reference velocity, C depicts the concentration, B0 denotes the uniform magnetic field, η0 denotes the magnetic diffusivity, σ denotes the electrical conductance, Cp represents specific heat, and τ* reveals the quotient of the latent heat of the nanoparticles to the latent heat of the base liquid. Furthermore, DB reflects the coefficient of Brownian motion, DT denotes the coefficient of thermophoresis diffusion, and the micro-inertial density is depicted by j. The viscosity of the spin gradient γ* is given as follows:
where
The Rosseland radiative heat flow (qr) along y-axis is defined as:
Here, Stefan–Boltzman constant and average absorption coefficient are expressed by σ∗ and k∗, respectively.
To simplify the analysis, the following appropriate similarity transformations are implemented.
Here, ψ and ψ1 are stream functions and can be represented as:
where η denotes the dimensionless variable, and f′(η), θ(η), ϕ(η), h(η), and g(η) are the dimensionless velocity curve, temperature curve, concentration curve, micropolar curve, and induced magnetic field curve, respectively. Implying Eqs. 23, 24, the continuity equation and magnetic flux equation holds true, and the Eqs. 12–16 are modified into the non-dimensional ODEs as stated:
The associated non-dimensional boundary constraints are given as:
The non-dimensional parameters are defined as:
Here, α1, α2 are the non-dimensional viscoelastic coefficients, β is the third-grade fluid constant, M depicts the magnetic parameter, δ denotes the thermal slip parameter , Ec represents the Eckert no, Rd is the radiation constant, Nt represents the thermophoresis constant, K denotes the micropolar fluid coefficient, Pr stands for Prandtl number, Nb indicates the Brownian diffusivity coefficient, Le denotes the Lewis no, B is the micro-inertia density coefficient, and Prm denotes the magnetic Prandtl number. Physical variables of interest, such as local couple stress (Ms), Sherwood number (Shx), Nusselt number (Nux), and skin friction coefficient (Cfx) are defined as follows:
where
Utilizing the similarity transformation (Eq. 22), the aforementioned expressions in non-dimensional aspects are stated as follows:
where X = x/l and
3 Numerical Procedure
Due to the extreme non-linearity, the coupled ODEs Eqs. 25–29 along with the boundary constraints (Eq. 30) are unable to solve analytically. To address non-linear boundary value problems, different numerical techniques are applied in MATLAB. The bvp4c approach is a useful technique to solve such problems numerically. The underlying partial differential equations (PDEs) are transfigured into ODEs utilizing the similarity analysis. The resulting system of ODEs is resolved numerically employing the built-in bvp4c technique in MATLAB. The solution strategy is described as follows:
Utilizing the aforementioned notations, the coupled ODEs are transformed into the following first-order ODEs:
The concerned boundary constraints in MATLAB script are stated as:
The MATLAB code is processed to get the numerical solutions. In order to determine the validity of numerical technique, we have analyzed our outcomes for limiting cases with already existing research work. The outcomes are reported to be in excellent concordance as shown in Tables 1, 2. It is also probed that for β = α2 = 0, the mathematical model reduces to the second-grade fluid model, and for β = α1 = α2 = 0, we get the equations of motion for the classical viscous liquid.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/964653/fphy-10-964653-HTML-r3/image_m/fphy-10-964653-t002.jpg)
TABLE 2. Comparative outcomes of − θ′(0) for diverse values of Pr, setting all other parameters zero.
4 Computational Results
In this article, the micropolar, third-grade, nanofluid flow past across an exponentially expanding surface is considered. The mathematical model is built by analyzing fluid flow assumptions. The system of dimensionless ODEs (Eqs. 25–29) are solved numerically, utilizing the bvp4c strategy in MATLAB in addition to the boundary constraints (Eq. 30), and the impacts of numerous physical variables are explored. The significance of these coefficients on the temperature curve θ(η), micro-rotation curve h(η), velocity curve f′(η), produced magnetic field profile g(η), and concentration profile ϕ(η) is highlighted via tables and figures.
4.1 Velocity Profile f′(η)
Figure 2 depicts the tendency of velocity distribution with the rising values of material constant α1. It is clear that velocity and momentum boundary layer thickness increases for the higher values of viscoelastic constant α1. Physically, the viscosity of the liquid is inversely proportional to the material constant α1 because of that when the stress is applied, it reduces the strain and strengthens the elastic effects between the adjacent layers, and hence the velocity profile enhances. Figure 3 denotes the impact of second material constant α2 on the velocity profile. The velocity curve is the decreasing function for larger values of α2. This parameter causes shear thickening of the fluid and an increase in resistance that reduces the boundary layer flow, and originates a decrement in the size of the momentum boundary layer width. Figure 4 exhibits the ascending behavior of velocity with the augmentation of the third-grade fluid coefficient β. This coefficient is inversely proportionate to the square of the liquid viscosity. Alternatively, higher β readings indicate superior third-grade material characteristics (higher liquid elasticity) and smaller fluid viscosity. This causes the boundary layer stream to accelerate, resulting in higher f′(η) values. As this value is raised, the fluid needs less stress to flow, encouraging flow acceleration. In Figure 5, the influence of the magnetic field coefficient on velocity profile is examined. It is analyzed that the velocity curve declines monotonically with the augmented values of M ,and the velocity diminishes far away from the sheet. This process assists in controlling the size of the boundary layer. This is because of the fact that the existence of magnetism in an electrically conductive liquid generates a force known as the Lorentz force that operates opposite to the flow direction and forces velocity profile to decline. The velocity diminishes as the retardation to the flow enhances. In Figure 6 the influence of micropolar fluid parameter on the structure of velocity distribution is noticed. The increase in vortex viscosity strongly accelerates the fluid flow. The results show that the momentum transfer layer-by-layer is significantly affected by the rise in viscosity caused by the micro-rotation of the molecules. The micro-elements rotate more strongly, which helps to accelerate the liquid motion in the boundary layer. As a result, linear momentum diffusion is enhanced by micro-polarity, which explains why suspension fluids have thinner boundary layers than regular fluids.
4.2 Temperature Profile θ(η)
Figure 7 explains that temperature in the boundary layer increases with progressing micropolar coefficient, that is, higher K values. The regime is greatly heated, and the size of the temperature boundary layer enhances. The vortex’s enhanced viscosity promotes thermal diffusion and serves as a rotator. This boosts up the capacity of thermal diffusion within the fluid’s regime from the micro to the macro level and rapidly carries heat from the sheet boundary into the liquid body with greater intensity. Figure 8 depicts the relationship between the temperature profile and the Prandtl number Pr. The fraction of the momentum diffusion coefficient to heat diffusion coefficient is called the Prandtl number. This is examined that the temperature curve declines with increasing the Prandtl number. The higher values of Prandtl number affect the thermal diffusion. Prandtl number is inversely proportionate to the thermal diffusion; hence, higher Pr values indicate lower thermal diffusion, leading in lower temperature and a weaker temperature boundary layer. When the Prandtl number increases, rate of thermal conductivity gets lower. Consequently, heat is dissipated more quickly, and hence, the temperature boundary layer width and temperature of the liquid both diminish. As a result, the Prandtl number is employed to enhance the cooling tendency of fluids. The temperature-dependent thermal conductivity parameter’s effect over the temperature curve is explored in Figure 9. This is evident that the temperature curve is increasing with the augmentation of ϵ. The consequence of various Eckert number values over the temperature curve is explored in Figure 10. Ec is known as the quotient of the kinetic energy and enthalpy of heat transfer. As the Eckert number (Ec) increases, the liquid’s temperature rises. The fluctuation in the thermal curve θ(η) is examined in Figure 11 in respect to distinct amounts of the radiation parameter Rd. It is identified that a higher thermal radiation constant is related to a greater temperature and a wider temperature boundary layer. As the temperature field grows, the enhanced radiation transfers a significant amount of heat to the fluid. Physically, the involved radiation produces more heat which raises the fluid’s temperature. The influence of thermophoresis coefficient over the thermal distribution is probed in Figure 12. Thermophoresis is a phenomenon that occurs in mixture of sub-micron sized particles, in which the different particles respond differently to the force of a thermal gradient. The particle’s velocity is known as the thermophoretic velocity, and the stress exerted on the dispersed particles because of the temperature difference is known as the thermophoretic force. The thermophoretic force increases as the value of Nt grows, causing the temperature to increase.
4.3 Concentration Profile ϕ(η)
The influence of Brownian diffusivity and thermophoresis coefficients over the concentration distribution ϕ(η) is portrayed in Figures 13, 14. According to these graphs, nanofluid constants have opposite impacts on concentration distributions. Particularly, as the thermophoresis constant rises, the width of the concentration boundary layer enhances; nevertheless, as the Brownian motion constant increases, the ϕ(η) values decrease. The thermophoresis constant amplifies the thermophoretic force, resulting in the transfer of nanoparticles from warm to cool locations and an increment in nanoparticle volume. Thermophoresis has numerous uses, including radioactive particle deposition in nuclear reactors, silicon thin film deposition, and aerosol technologies. Furthermore, the progressing amount of the Brownian motion coefficient reduces the micro-mixing of nanoparticles into the fluid’s zone, which diminishes the boundary layer thickness of concentration distribution.
The behavior of ϕ(η) for the numerous values of Lewis number (Le) is examined in Figure 15. The Lewis number is stated as the rate of heat-to-mass diffusion coefficient. It is utilized to express the flow of liquid, in which heat and momentum transfer occur simultaneously. The concentration distribution becomes steeper when Lewis number is increased. The higher values of ‘Le’ imply the lesser values of mass diffusivity ‘DB,’ which causes a weaker penetration depth for the concentration boundary layer. In Figure 16, the concentration curve for the various amounts of the micropolar coefficient is explored. However, when K increases, the micro-rotation velocity ϕ(η) drops down within a shorter range in the concentration boundary layer.
4.4 Micropolar Profile h(η)
Figures 17, 18 are portrayed for the different values of micro-inertia density coefficient B and micropolar coefficient K. This is clear that the micropolar (angular) speed of the sub-micron sized particles decreases with the increase in the micro-polarity of the fluid. Similarly, the micro-inertia density coefficient B serves as a restricting force for the micro-rotation profile h(η). So, for the progressing amounts of B, the micropolar boundary layer width declines. This is possible because micropolar fluids provide a high barrier to fluid motion.
4.5 Induced Magnetic Field Profile h(η)
The influence of magnetic Prandtl number Prm over induced magnetic field curve g(η) is exhibited in Figure 19 , where a rise in Prm leads to an improvement in the produced magnetic field distribution. This is because the magnetic diffusivity over the boundary layer’s surface decreased while the fluid viscous dispersion rate enhanced. Magnetic Prandtl no.(Prm) is a non-dimensional quantity in Magnetohydrodynamics that estimates the ratio of momentum diffusion coefficient (ν) to magnetic diffusion coefficient (η).
4.6 Numerical Results
The influence of all the non-dimensional coefficients used in the considered article on the skin friction parameter
Table 4 depicts the influence of different parameters on the Nusselt number and Sherwood number. The other parameters are taken to be fixed β = 0.5, α1 = 0.1, B = 0.5, α2 = 0.5, K = 0.2, and M = 0.3. The increment in the values of Prandtl number causes a decline in the thermal diffusivity, and hence resists the rise in the heat transmission rate at the boundary. Boosting the Prandtl number raises the average Nusselt number at the heated surface, while a reverse effect is examined for the
The higher values of thermophoresis coefficient Nt causes a decline in both the
5 Conclusion
The third-grade micropolar fluid flowing over an exponentially stretched sheet is probed in this research article. The underlying PDEs are transfigured into a system of ODEs using the appropriate similarity analysis, and the subsequent ODEs are settled in MATLAB using the bvp4c technique. The graphical explanation for the velocity curve f′(η), concentration curve ϕ(η), micropolar curve h(η), temperature curve θ(η), and induced magnetic field curve g(η) is portrayed via graphs. The key points under the aforementioned study are:
• Nanofluid constants have reverse effects over the concentration distribution. Concentration distribution ϕ(η) continues to increase for the greater values of Nt. It is due to the fact that the thermophoretic force increases causing an increase in concentration. On the other hand, an opposing trend is seen for Nb.
• Fluid’s velocity enhances with increasing the viscoelastic parameter α1, micropolar parameter K, and third-grade fluid coefficient β. The reverse trend is analyzed for magnetic field parameter M and α2.
• The thermal boundary layer width declines with the increment in the micropolar coefficient K and Prandtl number Pr. But the temperature gradient show ascending behavior for higher values of ϵ, radiation parameter Rd, Eckert number Ec, and thermophoresis constant Nt.
• The micropolar distribution h(η) tends to decline with increasing micropolar parameter K and micro-rotation constant B. With increasing the magnetic Prandtl number Prm, the induced magnetic field’s velocity remains going up.
• The skin friction coefficient remains diminishing for the larger values of α1, α2, and M, while, reverse behavior is observed for β, K, and B values.
• Couple stress coefficient continues to decline with rising the parameters α2, K, M, and B. The ascending trend is noted for the higher values of α1 and β.
• Nusselt number continuous to increase with boosting the parameters Pr , K, and Rd, while, the decreasing trend is examined for the parameters ϵ, Ec, Nt, and Le.
• Sherwood number depreciates for higher values of Pr , K, and Nt, and grows up for the higher values of ϵ, Ec, Rd, Nb, and Le.
Data Availability Statement
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author.
Author Contributions
AUA: methodology and conceptualization. AAA: writing—original draft and validation. FG: funding acquisition, software, and revision. EMT-E: validation and funding acquisition. AmA: software and funding acquisition. HH contributed in revision and validation.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors, and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Acknowledgments
The authors would like to extend their appreciation to the Deanship of Scientific Research at King Khalid University, Saudi Arabia for funding this work through the Research Group Program under Grant No. RGP. 2/12/43. We will be thankful to you for your kind favor.
References
1. Alfvén H. Existence of Electromagnetic-Hydrodynamic Waves. Nature (1942) 150(3805):405–6. doi:10.1038/150405d0
2. Andersson HI. MHD Flow of a Viscoelastic Fluid Past a Stretching Surface. Acta Mech (1992) 95(1):227–30. doi:10.1007/bf01170814
3. Abbas Z, Hayat T. Radiation Effects on MHD Flow in a Porous Space. Int J Heat Mass Transf (2008) 51(5-6):1024–33. doi:10.1016/j.ijheatmasstransfer.2007.05.031
4. Nadeem S, Hussain A. MHD Flow of a Viscous Fluid on a Nonlinear Porous Shrinking Sheet with Homotopy Analysis Method. Appl Math Mech.-Engl Ed (2009) 30(12):1569–78. doi:10.1007/s10483-009-1208-6
5. Ahmad I, Sajid M, Awan W, Rafique M, Aziz W, Ahmed M, et al. MHD Flow of a Viscous Fluid over an Exponentially Stretching Sheet in a Porous Medium. J Appl Maths (2014) 2014:1–8. doi:10.1155/2014/256761
6. Khan MI, Tamoor M, Hayat T, Alsaedi A. MHD Boundary Layer thermal Slip Flow by Nonlinearly Stretching cylinder with Suction/blowing and Radiation. Results Phys (2017) 7:1207–11. doi:10.1016/j.rinp.2017.03.009
7. Sohail M, Nazir U, Chu YM, Al-Kouz W, Thounthong P. Bioconvection Phenomenon for the Boundary Layer Flow of Magnetohydrodynamic Carreau Liquid over a Heated Disk. Sci Iran (2021) 28(3):1896–907. doi:10.24200/SCI.2021.53970.3518
8. Riaz A, Abbas T, Zeeshan A, Doranehgard MH. Entropy Generation and MHD Analysis of a Nanofluid with Peristaltic Three Dimensional Cylindrical Enclosures. Hff (2021) 31:2698–714. doi:10.1108/HFF-11-2020-0704
9. AbdulHussein WA, Abed AM, Mohammed DB, Smaisim GF, Baghaei S. Investigation of Boiling Process of Different Fluids in Microchannels and Nanochannels in the Presence of External Electric Field and External Magnetic Field Using Molecular Dynamics Simulation. Case Stud Therm Eng (2022) 35:102105. doi:10.1016/j.csite.2022.102105
10. Hussain M, Ashraf M, Nadeem S, Khan M. Radiation Effects on the thermal Boundary Layer Flow of a Micropolar Fluid towards a Permeable Stretching Sheet. J Franklin Inst (2013) 350(1):194–210. doi:10.1016/j.jfranklin.2012.07.005
11. Abbas N, Saleem S, Nadeem S, Alderremy AA, Khan AU. On Stagnation point Flow of a Micro Polar Nanofluid Past a Circular cylinder with Velocity and thermal Slip. Results Phys (2018) 9:1224–32. doi:10.1016/j.rinp.2018.04.017
12. Awan AU, Abid S, Ullah N, Nadeem S. Magnetohydrodynamic Oblique Stagnation point Flow of Second Grade Fluid over an Oscillatory Stretching Surface. Results Phys (2020) 18:103233. doi:10.1016/j.rinp.2020.103233
13. Awan AU, Abid S, Abbas N. Theoretical Study of Unsteady Oblique Stagnation point Based Jeffrey Nanofluid Flow over an Oscillatory Stretching Sheet. Adv Mech Eng (2020) 12(11):1687814020971881. doi:10.1177/1687814020971881
14. Khan MN, Nadeem S. A Comparative Study between Linear and Exponential Stretching Sheet with Double Stratification of a Rotating Maxwell Nanofluid Flow. Surf Inter (2021) 22:100886. doi:10.1016/j.surfin.2020.100886
15. Abdelmalek Z, Khan SU, Waqas H, Riaz A, Khan IA, Tlili I. A Mathematical Model for Bioconvection Flow of Williamson Nanofluid over a Stretching cylinder Featuring Variable thermal Conductivity, Activation Energy and Second-Order Slip. J Therm Anal Calorim (2021) 144(1):205–17. doi:10.1007/s10973-020-09450-z
16. Abbasbandy S, Hayat T, Mahomed FM, Ellahi R. On Comparison of Exact and Series Solutions for Thin Film Flow of a Third-Grade Fluid. Int J Numer Meth Fluids (2009) 61(9):987–94. doi:10.1002/fld.1994
17. Hayat T, Naz R, Alsaedi A, Rashidi MM. Hydromagnetic Rotating Flow of Third Grade Fluid. Appl Math Mech.-Engl Ed (2013) 34(12):1481–94. doi:10.1007/s10483-013-1761-7
18. Hayat T, Aziz A, Muhammad T, Ahmad B. On Magnetohydrodynamic Flow of Second Grade Nanofluid over a Nonlinear Stretching Sheet. J Magnetism Magn Mater (2016) 408:99–106. doi:10.1016/j.jmmm.2016.02.017
19. Mahanthesh B, Joseph TV. Dynamics of magneto-nano third-grade fluid with Brownian motion and thermophoresis effects in the pressure type die. J Nanofluids (2019) 8(4):870–5. doi:10.1166/jon.2019.1642
20. Riaz A, Ellahi R, Sait SM, Muhammad T. Magnetized Jeffrey Nanofluid with Energy Loss in between an Annular Part of Two Micro Non-concentric Pipes. Energy Sourc A: Recovery, Utilization, Environ Effects (2020) 2020:1–20. doi:10.1080/15567036.2020.1798568
21. Mondal H, De P, Goqo S, Sibanda P. A Numerical Study of Nanofluid Flow over a Porous Vertical Plate with Internal Heat Generation and Nonlinear thermal Radiation. J Por Media (2020) 23(6):517–29. doi:10.1615/JPorMedia.2020026624
22. Sangeetha E, De P. Bioconvection in Nanofluid Flow Embedded in Non-darcy Porous Medium with Viscous Dissipation and Ohmic Heating. J Por Media (2021) 24(1):15–23. doi:10.1615/JPorMedia.2020036165
23. Riaz A, Zeeshan A, Bhatti MM. Entropy Analysis on a Three-Dimensional Wavy Flow of Eyring-Powell Nanofluid: A Comparative Study. Math Probl Eng (2021) 2021:1–14. doi:10.1155/2021/6672158
24. Eringen A. Theory of Micropolar Fluids. Indiana Univ Math J (1966) 16(1):1–18. doi:10.1512/iumj.1967.16.16001
25. Subba Reddy Gorla R, Mansour MA, Mohammedien AA. Combined Convection in an Axisymmetric Stagnation Flow of Micropolar Fluid. Int J Numer Methods Heat Fluid Flow (1996) 6(4):47–55. doi:10.1108/09615539610123441
26. Rahman MM, Sattar MA. Magnetohydrodynamic Convective Flow of a Micropolar Fluid Past a Continuously Moving Vertical Porous Plate in the Presence of Heat Generation/absorption. J Heat Transfer (2006) 128(2):142–52. doi:10.1115/1.2136918
27. Ishak A. Thermal Boundary Layer Flow over a Stretching Sheet in a Micropolar Fluid with Radiation Effect. Meccanica (2010) 45(3):367–73. doi:10.1007/s11012-009-9257-4
28. Gaffar SA, Bég OA, Prasad VR. Mathematical Modeling of Natural Convection in a Third-Grade Viscoelastic Micropolar Fluid from an Isothermal Inverted Cone. Iran J Sci Technol Trans Mech Eng (2020) 44(2):383–402. doi:10.1007/s40997-018-0262-x
29. Ali L, Liu X, Ali B, Mujeed S, Abdal S, Mutahir A. The Impact of Nanoparticles Due to Applied Magnetic Dipole in Micropolar Fluid Flow Using the Finite Element Method. Symmetry (2020) 12(4):520. doi:10.3390/sym12040520
30. Jiang Y, Smaisim GF, Mahmoud MZ, Li Z, Aybar HŞ, Abed AM. Simultaneous Numerical Investigation of the Passive Use of Phase-Change Materials and the Active Use of a Nanofluid inside a Rectangular Duct in the thermal Management of Lithium-Ion Batteries. J Power Sourc (2022) 541:231610. doi:10.1016/j.jpowsour.2022.231610
31. Fosdick RL, Rajagopal KR. Thermodynamics and Stability of Fluids of Third Grade. Proc R Soc Lond A (1980) 369(1738):351–77. doi:10.1098/rspa.1980.0005
32. Papautsky I, Brazzle J, Ameel T, Frazier AB. Laminar Fluid Behavior in Microchannels Using Micropolar Fluid Theory. Sens Actuator A Phys (1999) 73(1-2):101–8. doi:10.1016/S0924-4247(98)00261-1
33. Sajid M, Hayat T, Asghar S. Non-similar Analytic Solution for MHD Flow and Heat Transfer in a Third-Order Fluid over a Stretching Sheet. Int J Heat Mass Transfer (2007) 50(9-10):1723–36. doi:10.1016/j.ijheatmasstransfer.2006.10.011
34. Pakdemirli M. The Boundary Layer Equations of Third-Grade Fluids. Int J Non-Linear Mech (1992) 27(5):785–93. doi:10.1016/0020-7462(92)90034-5
35. Loganathan K, Mohana K, Mohanraj M, Sakthivel P, Rajan S. Impact of Third-Grade Nanofluid Flow across a Convective Surface in the Presence of Inclined Lorentz Force: an Approach to Entropy Optimization. J Therm Anal Calorim (2021) 144(5):1935–47. doi:10.1007/s10973-020-09751-3
36. Bég OA, Bakier AY, Prasad VR, Zueco J, Ghosh SK. Nonsimilar, Laminar, Steady, Electrically-Conducting Forced Convection Liquid Metal Boundary Layer Flow with Induced Magnetic Field Effects. Int J Therm Sci (2009) 48(8):1596–606. doi:10.1016/j.ijthermalsci.2008.12.007
37. Chaudhary S, Singh S, Chaudhary S. Thermal Radiation Effects on MHD Boundary Layer Flow over an Exponentially Stretching Surface. Appl Maths (2015) 06(02):295–303. doi:10.4236/am.2015.62027
38. Hayat T, Khan MI, Waqas M, Alsaedi A, Yasmeen T. Diffusion of Chemically Reactive Species in Third Grade Fluid Flow over an Exponentially Stretching Sheet Considering Magnetic Field Effects. Chin J Chem Eng (2017) 25(3):257–63. doi:10.1016/j.cjche.2016.06.008
39. Elbashbeshy E. Heat Transfer over an Exponentially Stretching Continuous Surface with Suction. Arch Mech (2001) 53(6):643–51.
40. Ali B, Naqvi RA, Hussain D, Aldossary OM, Hussain S. Magnetic Rotating Flow of a Hybrid Nano-Materials Ag-MoS2 and Go-MoS2 in C2H6O2-H2o Hybrid Base Fluid over an Extending Surface Involving Activation Energy: FE Simulation. Mathematics (2020) 8(10):1730. doi:10.3390/math8101730
41. Kumari M, Nath G. Transient Rotating Flow over a Moving Surface with a Magnetic Field. Int J Heat Mass Transfer (2005) 48(14):2878–85. doi:10.1016/j.ijheatmasstransfer.2005.02.009
Nomenclature
v Velocity in y-direction (m/s)
u Velocity in x-direction (m/s)
Cp Specific heat (J/kg.K)
C∞ Ambient concentration (kg/m3)
T∞ Ambient temperature (K)
Rd Radiation parameter (–)
k Vortex viscosity (Pa.s)
qr Radiative heat flux (W/m2)
p Pressure (N/m2)
Re Reynolds number (–)
T Temperature (K)
N Angular velocity (1/s)
Nt Thermophoresis diffusion parameter (–)
Ec Eckert number (–)
Prm Magnetic Prandtl number (–)
Cfx Skin friction coefficient (–)
Greek Symbols
τ* Ratio of latent heat capacities (–)
σ* Stepan–Boltzmann constant (W/m2.K4)
μ Viscosity (kgm−1s−1)
θ(η) Dimensionless temperature (–)
ηo Magnetic diffusivity (m2/s)
ρ Density (kg.m−3)
η Dimensionless variable
j Micro-inertia density (m2)
α1, α2, β Dimensionless fluid parameters (–)
(x, y) Cartesian coordinates (m)
Tf Temperature at wall (K)
C Fluid’s concentration (kg/m3)
Uo Reference velocity (m/s)
B Micro-inertia density coefficient (–)
DB Brownian diffusivity coefficient (m2/s)
DT Thermophoresis diffusivity coefficient (m2/s)
Bo Applied magnetic field (kg/s2.A)
Pr Prandtl number (–)
Le Lewis number (–)
K micro-polarity parameter (–)
Nb Brownian diffusion parameter (–)
M Magnetic parameter
H1, H2 Induced magnetism components (A/m)
Cf Concentration at sheet (kg/m3)
Nux Heat transfer coefficient (–)
f(η) Non-dimensional velocity (–)
ϕ(η) Non-dimensional concentration curve (–)
δ Thermal slip parameter (–)
g(η) Dimensionless induced field (–)
k* Average absorption coefficient (1/m)
h(η) Dimensionless micro-rotation curve (–)
γ* Fluid’s spin gradient viscosity (kg.m/s)
σ Electrical conductance (A2.s3/kg.m3)
Ψ(x, y) Stream function (m/s)
Keywords: magnetohydrodynamics, micropolar third-grade fluid, stretching sheet, Buongiorno model, bvp4c technique
Citation: Awan AU, Akbar AA, Hamam H, Gamaoun F, Tag-ElDin EM and Abdulrahman A (2022) Characterization of the Induced Magnetic Field on Third-Grade Micropolar Fluid Flow Across an Exponentially Stretched Sheet. Front. Phys. 10:964653. doi: 10.3389/fphy.2022.964653
Received: 08 June 2022; Accepted: 20 June 2022;
Published: 15 July 2022.
Edited by:
Arshad Riaz, University of Education Lahore, PakistanReviewed by:
Poulomi De, Vellore Institute of Technology (VIT), IndiaGhassan F. Smaisim, University of Kufa, Iraq
Copyright © 2022 Awan, Akbar, Hamam, Gamaoun, Tag-ElDin and Abdulrahman. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Aziz Ullah Awan, aziz.math@pu.edu.pk