Skip to main content

ORIGINAL RESEARCH article

Front. Phys., 15 July 2022
Sec. Interdisciplinary Physics

Characterization of the Induced Magnetic Field on Third-Grade Micropolar Fluid Flow Across an Exponentially Stretched Sheet

Aziz Ullah Awan
Aziz Ullah Awan1*Asia Ali AkbarAsia Ali Akbar1Haneen HamamHaneen Hamam2Fehmi GamaounFehmi Gamaoun3ElSyed M. Tag-ElDinElSyed M. Tag-ElDin4Amal AbdulrahmanAmal Abdulrahman5
  • 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:

.V=0,(1)
divT=J×B+ρDVDtρb.(2)

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:

T+pI=S,(3)
S=μB1+α1̂B2+α2̂B12+β1̂B3,(4)
+β2̂(B1B2+B2B1)+β3̂B1trcB12.

Here, p, I, T, and S are the pressure, identity tensor, Cauchy stress tensor, and the extra stress tensor, respectively. Furthermore, αk̂(k=1,2),βĵ(j=1,2,3) are metallic constants and Bi(i = 1, 2, 3) is the kinematic tensor defines as:

B1=LT+L,Bn=DBn−1Dt+Bn−1L+LTBn−1,n=2,3where,L=V.

Here, DDt represents the substantial derivative, and it is stated as:

D()DtV.()+()t.

In the case of third-grade fluid, the material moduli satisfy Clausius–Duhem inequality stated as:

β1̂=β2̂=0,α1̂0,μ0,β3̂0,(5)
α2̂+α1̂24μβ3̂.(6)

So, the stress tensor takes the following form.

T=μB1+α1̂B2+α2̂B12+β3̂B1trcB12.(7)

2.2 Constitutive Model for Micropolar Fluid

The field equations for the micropolar fluid following Papautsky et al. [32] are stated as follows:

.ρV+ρt=0,(8)
λ+κ+2μ.V+κ×Gκ+μ××V+ρf=ρV+p,(9)
ᾱ+γ+β̄.G+κ×V2κGγ××G=ρl+ρjG.(10)

Eqs. 810, 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 ᾱ, γ∗, and β̄ denote the swirl gradient viscosity coefficient. Also, (′) denotes the time derivative. For ᾱ=κ=β̄=γ=0 and in the absence of l and f, the micro-rotation vector G approaches zero, and Eq. 9 is simplified to the Navier–Stokes equations. For κ = 0, the gyration vector G and velocity vector V become detached, and the micro-movements have no influence on the entire motion of the fluid.

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.

FIGURE 1
www.frontiersin.org

FIGURE 1. Geometry of the flow.

Using Eq. 7 in Eq. of motion (Eq. 2) along with the boundary layer estimations [3336] in case of third-grade fluid, notably, inside the boundary layer px,2ux,ux, and u are O(1), v and y are O(δ), αĵρ(j=1,2) and ν be O(δ2), and βk̂ρ(k=1,2,3) being O(δ4), and the components of O(δ) are ignored (δ is boundary layer width). The equations for conservation of mass, magnetic field, and momentum are given as Eqs. 1112. Equations 1314 rise from the conservation of energy and concentration. Eq. 15 rises from micropolar boundary layer approximation by Papautsky et al. [32], and Eq. 16 resulted from the conservation of induced magnetic field. Hence, the mathematical model overseeing the considered flow problem is as follows:

vy+ux=0,H1x+H2y=0,(11)
vuy+uux=α1̂ρux2uy2+3uy2vy2+v3uy3+u3uxy2+2α2̂ρuy2vy2(12)
+6β3̂ρ2uy2uy2+kρNy+k+μρ2uy2σB02ρu,
vTy+uTx=1ρCpyTykT+k+μρCpuy21ρCpqry(13)
+τDBCyTy+Ty2DTT
vCy+uCx=DB2Cy2+DTT2Ty2,(14)
vNy+uNx=γρj2Ny2kjρuy+2N,(15)
vH1y+uH1x=H2uy+H1ux+ηo2H1y2.(16)

The concerned boundary conditions are as follows:

N=muy,u=Uwx=U0expx/l,T=kTy+Tf,
H2=0,DBCy+DTTTy=0,v=0,H1y=0,fory0,(17)
CC,H1HexTT,N0,u0,uy0fory.(18)

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:

γ=jμ+k2=μj1+K2,andj=2lμρU0expx/l,(19)

where K=kμ is the micropolar parameter. Temperature-based thermal conductivity is defined as follows:

kT=1+ϵTTTTfk.(20)

The Rosseland radiative heat flow (qr) along y-axis is defined as:

qr=4σ3kT4y,andT44T3T3T4.(21)

Here, Stefan–Boltzman constant and average absorption coefficient are expressed by σ∗ and k∗, respectively.

qry=16T3σ3k2Ty2.(22)

To simplify the analysis, the following appropriate similarity transformations are implemented.

ψx,y=2νlU0fηexpx2l,ψ1x,y=H0νU0gηexpx2l,ϕη=CCCfC,Nη=U0U02νlexp3x2lhη,θη=TTTfT,η=yU02νlexpx2l.(23)

Here, ψ and ψ1 are stream functions and can be represented as:

u=ψy=U0fηexpx/l,v=ψx=νU02lf+ηfexpx/2l,H1=ψ1y=H02lexpx/lgη,H2=ψ1x=H0ν4l2U0expx/2lg+ηg.(24)

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. 1216 are modified into the non-dimensional ODEs as stated:

1+Kf+3βff2+Kh+ff+α12ηffffiv9f2+6ff2f22Mfα2ηff+3f2=0,(25)
1Prθ+ϵθθ+ϵθ2+fθ+Ec1+Kf2+43Rdθ+Ntθ2+Nbθϕ=0,(26)
ϕ+NtNbθ+Lefϕ=0,(27)
1+K2h+fh3hf2KBf+2h=0,(28)
g+Prmfggf=0.(29)

The associated non-dimensional boundary constraints are given as:

Nbϕ0+Ntθ0=0,g0=0,f0=0,h0=mf0,h=0,ϕ=0,θ0=δθ0+1,θ=0,f=0,f0=1,f=0,g=1,g0=0.(30)

The non-dimensional parameters are defined as:

α1=α1̂U02ρνlexpx/l,α2=α2̂U0ρνlexpx/l,β=β3U03ρν2lexp3x/l,Pr=μCpk,Ec=U02exp2x/lCpTfT,Rd=4σT3μCpk,Nt=DTTfTτνT,Nb=DBτνCfC,Le=νDB,B=νlU0jexpx/l,Prm=νη0,δ=kU02νlexpx/2lM=σB02lρU0expx/l,K=kμ.

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:

Shx=xjwDBCfC,Cfx=τw12ρUw2,Nux=xqwkTfT,Ms=γUw2ρlNyy=0,(31)

where

jw=DBCyy=0,qw=16T3σ3k+kTyy=0,
τw=μ+kuy+α1̂v2uy2+2uyux+u2uyx+2β3uy3+kNy=0.

Utilizing the similarity transformation (Eq. 22), the aforementioned expressions in non-dimensional aspects are stated as follows:

12Rex1/2Cfx=1+Kf0+Kh0+βf03+α1f0f0+7f0f0,(32)
2XRex1/2Nux=1+43Rdθ0,(33)
2XRex1/2Shx=ϕ0,(34)
RexMs=1+K2h0,(35)

where X = x/l and Rex=lUw(x)ν is the local Reynolds number.

3 Numerical Procedure

Due to the extreme non-linearity, the coupled ODEs Eqs. 2529 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:

fη=y1;fη=y2;fη=y3;fη=y4;fivη=yy1;θη=y5;θη=y6;θη=yy2;ϕη=y7;ϕη=y8;ϕη=yy3;hη=y9;hη=y10;hη=yy4;gη=y11;gη=y12;gη=y13;gη=yy5.

Utilizing the aforementioned notations, the coupled ODEs are transformed into the following first-order ODEs:

yy1=α1y111+Ky42y2y2+y3y12My2+3βy3y3y4+Ky109α1+3α2y3y3α2+2α1y4xy3+6α1y4y2;,(36)
yy2=3Pr3+3ϵy5+4PrRd1y1y6+1+Ky3Ecy3+y6Nty6+Nby6,y8+ϵPry6y6;(37)
yy3=Ley1y8NtNbyy2,(38)
yy4=22+K12KB2y9+y3+3y2y9y1y10,(39)
yy5=Prmy1y13+y11y3.(40)

The concerned boundary constraints in MATLAB script are stated as:

y01;y021;y05δy061;y09+my03;yinf5;yinf2;Nby08+Nty06;yinf9;yinf7;y011;y013;yinf121;yinf3.(41)

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.

TABLE 1
www.frontiersin.org

TABLE 1. Comparative results of − f(0) for diverse values of M, setting all other parameters zero.

TABLE 2
www.frontiersin.org

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. 2529) 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.

FIGURE 2
www.frontiersin.org

FIGURE 2. Effect of α1 on the velocity curve.

FIGURE 3
www.frontiersin.org

FIGURE 3. Effect of α2 on the velocity curve.

FIGURE 4
www.frontiersin.org

FIGURE 4. Effect of β on the velocity curve.

FIGURE 5
www.frontiersin.org

FIGURE 5. Effect of M on the velocity curve.

FIGURE 6
www.frontiersin.org

FIGURE 6. Effect of K on the velocity curve.

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.

FIGURE 7
www.frontiersin.org

FIGURE 7. Impact of K on the θ(η) curve.

FIGURE 8
www.frontiersin.org

FIGURE 8. Impact of Pr on the θ(η) curve.

FIGURE 9
www.frontiersin.org

FIGURE 9. Impact of ϵ on the θ(η) curve.

FIGURE 10
www.frontiersin.org

FIGURE 10. Impact of Ec on the θ(η) curve.

FIGURE 11
www.frontiersin.org

FIGURE 11. Impact of Rd on the θ(η) curve.

FIGURE 12
www.frontiersin.org

FIGURE 12. Impact of Nt on the θ(η) curve.

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.

FIGURE 13
www.frontiersin.org

FIGURE 13. Effect of Nb on the ϕ(η) profile.

FIGURE 14
www.frontiersin.org

FIGURE 14. Effect of Nt on the ϕ(η) profile.

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.

FIGURE 15
www.frontiersin.org

FIGURE 15. Effect of Le on the ϕ(η) profile.

FIGURE 16
www.frontiersin.org

FIGURE 16. Effect of K on the ϕ(η) profile.

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.

FIGURE 17
www.frontiersin.org

FIGURE 17. Effect of K on the h(η) profile.

FIGURE 18
www.frontiersin.org

FIGURE 18. Effect of B on the h(η) profile.

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 (η).

FIGURE 19
www.frontiersin.org

FIGURE 19. Influence of Prm on induced magnetic velocity.

4.6 Numerical Results

The influence of all the non-dimensional coefficients used in the considered article on the skin friction parameter (Rex1/2Cfx), Nusselt number. (Rex1/2Nux), Sherwood number. (Rex1/2Shx), and coupled wall stress (Rex Ms) are presented in Tables 3, 4. Table 3 probes the significance of involved coefficients on the skin friction. The values of the coefficients other than used in Table 3 are taken to be fixed Pr = 0.7, ϵ = 0.3, Ec = 0.1, K = 0.2, Rd = 0.3, Nt = 0.3, Nb = 1.0, and Le = 1.0. By increasing the value of material parameters α1 and α2, the skin friction coefficient reduces. It is because of the reason that the viscosity near the surface of sheet is enhanced with the progression in material constant, which causes the decline in Rex1/2Cfx values. As the rise in shear thickening coefficient β, improves the boundary layer thickness that give rise to the skin friction coefficient. The increase in micropolar parameter K, give rise to the value of Rex1/2Cfx, but the opposite behavior is seen for couple stress Rex Ms. A reduction is caused by a rise in the magnetic field constant M in both Rex1/2Cfx and Rex Ms values. By increasing the micro-inertia density B, the values of Rex1/2Cfx get enhanced but reverse behavior is seen for couple stress Rex Ms.

TABLE 3
www.frontiersin.org

TABLE 3. Numerical findings of RexMs and Rex1/2Cfx for various values of parameters

TABLE 4
www.frontiersin.org

TABLE 4. Numerical findings of Shx(Rex)0.5 and Nux(Rex)0.5 for diverse values of parameters.

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 Rex1/2Shx. As ϵ is enhanced, the Nusslet number decreases while the Sherwood number increases. The increase in the value of Ec causes the Nusselt number to decrease while the Sherwood number increases. The boosting values of micropolar parameter K enhances the Nusselt number but declines the Sherwood number. Higher values of the Rd parameter improve convective heat and matter transmission, which ultimately increases the average Nusselt and Sherwood numbers values.

The higher values of thermophoresis coefficient Nt causes a decline in both the NuxRex1/2 and ShxRex1/2 values. The increment in the Nb values results in the enhancement of Sherwood number because the mass transfer rate from higher to lower concentration area get increased. The growth in the values of Lewis number Le improves the thermal diffusion and declines the mass conductivity. Thus, the values of NuxRex1/2 decrease and that of ShxRex1/2 increase.

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

CrossRef Full Text | Google Scholar

2. Andersson HI. MHD Flow of a Viscoelastic Fluid Past a Stretching Surface. Acta Mech (1992) 95(1):227–30. doi:10.1007/bf01170814

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

24. Eringen A. Theory of Micropolar Fluids. Indiana Univ Math J (1966) 16(1):1–18. doi:10.1512/iumj.1967.16.16001

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

39. Elbashbeshy E. Heat Transfer over an Exponentially Stretching Continuous Surface with Suction. Arch Mech (2001) 53(6):643–51.

Google Scholar

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

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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, Pakistan

Reviewed by:

Poulomi De, Vellore Institute of Technology (VIT), India
Ghassan 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

Disclaimer: 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.