# A Numerical Assessment of Time-Dependent Magneto-Convective Thermal-Material Transfer over a Vertical Permeable Plate

Authors: Mohammed Jahir Uddin; R. Nasrin
Journal: Journal of Applied Mathematics
 (2023)

---

## Abstract

The objective of this work is to investigate the influences of thermal radiation, heat generation, and buoyancy force on the time-dependent boundary layer (BL) flow across a vertical permeable plate. The fluid is unsteady, incompressible, viscous, and electrically insulating. The heat transfer mechanism happens due to free convection. The nondimensional partial differential equations of continuity, momentum, energy, and concentration are discussed using appropriate transformations. The impressions of thermal radiation and buoyancy forces are exposed in the energy and momentum equation, respectively. For numerical model, a set of nonlinear dimensionless partial differential equations can be solved using an explicit finite difference approach. The stability and convergence analyses are also established to complete the formulation of the model. The thermophysical effects of entering physical parameters on the flow, thermal, and material fields are analyzed. The variations in local and average skin friction, material, and heat transfer rates are also discussed for the physical interest. The analysis of the obtained findings is shown graphically, and relevant parameters pointedly prejudice the flow field. Studio Developer FORTRAN 6.2 and Tecplot 10.0 are applied to simulate the schematic model equations and graphical presentation numerically. The intensifying values of the magnetic field are affected decreasingly in the flow field. The temperature profiles decrease within the BL to increase the value of radiation parameters. The present study is on the consequences for petroleum engineering, agriculture engineering, extraction, purification processes, nuclear power plants, gas turbines, etc. To see the rationality of the present research, we compare these results and the results available in the literature with outstanding compatibility.


---

## Body

## 1. Introduction

MHD is the magnetic field that may produce currents in a fluid that magnetizes the fluid and alters the magnetic field. It is utilized in a variety of areas, including astrophysics, geophysics, sensors, earthquakes, engineering (such as plasma confinement, cooling of nuclear reactors, liquid-metal and polymer technology, and electromagnetic casting), biomedical engineering (like magnetic drug targeting, magnetic devices for cell separation, magnetic endoscopy, cancer tumor treatment, and cell death by hyperthermia), and meteorology, which is created by an alternating magnetic field. Many fields of research and engineering employ porous media, including filtration mechanics, soil mechanics, rock mechanics, groundwater hydrology and petroleum engineering, material science, geoscience, biology and biophysics, and metal forming. Heat convection in a porous medium serves the purposes of energy recovery, thermal energy storage, geothermal oil extraction, natural environment ventilation, central heating, and the flow of filtration devices. The simultaneous occurrence of heat and mass transfer in a moving fluid is crucial to designing chemical processing equipment, nuclear power plants, gas turbines, and thermal energy storage. The free convection flow often occurs in nature; fluid flows through porous media with magnetic fields are of particular interest and have attracted many scientists and engineers owing to their applications in science and industry.Ahmed and Alam [1] investigated the finite-difference solution of magnetohydrodynamics (MHD) mixed convection flow with chemical reaction and heat production. Ahmed et al. [2] studied naturally convective flow in water-copper nanofluid-filled circular and arc cavities. Olajuwon and Oahimire [3] investigated the effects of thermal radiation and Hall current on heat and mass transfer of unsteady MHD flow in saturated porous media containing a viscoelastic micropolar fluid. Das and Jana [4] studied natural convective magneto-nanofluid flow and thermal transmission across a vertically moving plate. Samiulhaq [5] investigated the influence of MHD free convection flow and a porous medium on thermal diffusion and ramping wall temperature. Reddy [6] observed unsteady MHD convection heat and mass transfer flow over a vertical semi-infinite porous plate with changing viscosity and thermal conductivity. Makinde and Aziz [7] examined the effects of MHD mixed convection flow across a vertical plate contained in a porous media and subjected to a transverse magnetic field. Reddy et al. [8] investigated the numerical solution of transient natural convection MHD radiation flow through a vertical cylinder contained in a porous medium with varying surface temperatures and concentrations.Ali et al. [9] conducted a numerical simulation of heat and mass transfer over a stretched wedge surface and identified the nanoparticle concentration solution for Brownian motion less than 0.2, thermophoresis less than 0.14, and Lewis number more than 1. Tarammim et al. [10] observed MHD free convection flow over a vertical plate. Alam et al. [11] obtained the numerical solution of coupled free-forced convection and mass flow through a vertical permeable medium with heat production and thermal transmission. Shamshuddin et al. [12] explored the effects of the thermal Péclet number, vortex viscosity, and Reynolds number on the two-dimensional flow of micropolar fluid through a channel owing to mixed convection. Ullah et al. [13] demonstrated magneto-convective flow over a vertical plate with radiation. Islam et al. [14] triumphed over the mathematical forming of a flow field with uniform/nonuniform thermal and magnetic strength in a half-moon-shaped region. The authors found that after the nondimension time of about 0.65, the numerical solution reached a steady-state phase from unsteadiness, and the rate of thermal transmission rises for rising buoyancy force whereas devalued with the elevated magnetic field. Magneto-convective flow over an inclined plate with Hall current was conquered by Alam et al. [15].Abdullameed et al. [16] revealed closed from solutions for unsteady MHD flow in a permeable medium with wall transpiration. Shamshuddin and Ghaffari [17] investigated radiative heat energy on Casson-type nanoliquid generated by a convectively heated porous medium in combination with thermophoresis and Brownian motion effects. Chamkha et al. [18] observed turbulent spontaneous convective power-law fluid flow across a vertical plate embedded in a non-Darcian porous medium with a homogenous chemical reaction.The innovation of this investigation is to explore the numerical solutions of magneto-convective heat and mass fluid flow over a vertical permeable plate with thermophysical parameters. The main novelty of this research is considering the influence of nine parameters and dimensionless numbers combinedly with the permeable plate. Also, it is required to study the inclusion of expressions on flow, thermal, and material profiles, as well as local and average Sherwood numbers, Nusselt number, and shear stress fluctuations in thermophysical and hydrodynamical parameters. In the first stage, the governing equations are established and mathematical analysis is given. The results and discussions including the graphs, tables, and comparisons are presented in the second stage. Finally, the significant outcomes of this investigation are summarized. Thereafter, future works and significant references related to this research work are included.


## 2. Formulation of Problem

The presentation of the two-dimensional unsteady fluid model with magneto-convective heat and mass transfer is analyzed in this research by adding thermal radiation and heat generation. According to Ali et al. [19], our research schedule is displayed in Figure 1.Figure 1
Flow chart.


### 2.1. Physical Model

Consider the trembling two-dimensional nonlinear MHD free convection flows of a viscoelastic, deformable, and electrically insulating fluid around a vertical plate embedded in porous media. This is happening while the thermal and material buoyancy effects are being encouraged. Figure2 depicts the entire flow mechanism structure with the appropriate coordinates. The x-axis is placed along the porous plate in a vertically upward orientation. Consequently, U is a function of y and t only, i.e., Uy,t. Let u and v represent the velocity components in the x and y directions, respectively. Initially, both the plate and fluid are the same at the temperature. Instantaneously, the plate temperature and species concentration are raised to Tw>T∞ and Cw>C∞, respectively.Figure 2
Physical configuration and coordinate system.Typically, a consistent magnetic fieldB=0,B0,0 is applied to the plate. All assumptions and approximations for the model flow are given as follows:
(i)
The induced magnetic field is insignificant compared to the applied magnetic field(ii)
The Joule heating effects are negligible in the energy equation(iii)
The magnetic Reynolds numbers and effective transversal magnetic field are insignificant(iv)
No applied electric force is expected(v)
The electric field is assumed to be zero(vi)
Viscous dissipation is used to quantify the temperature gradient of fluid particles within the surface of a flat plate(vii)
Flow in the medium is caused by the buoyancy force generated by the temperature gradient between the fluid and the porous plate(viii)
Under the Rosseland [20] approximation of radiative heat flow with Taylor series expansion, the thermal radiation is employed in a way that may be well estimatedUnder the above hypotheses and Boussinesq’s approximation and following Nasrin and Alim [21, 22], the equations of continuity, momentum, energy, and concentration can be put mathematically in the following form.
(1)∂u∂x+∂v∂y=0,(2)∂u∂t+u∂u∂x+v∂u∂y=ν∂2u∂y2+gβT−T∞+gβ∗C−C∞−σB02ρu−νκu,(3)∂T∂t+u∂T∂x+ν∂T∂y=κρCp.∂2T∂y2+νCp∂u∂y2−1ρCp∂qr∂y+Q0ρCPT−T∞,(4)∂C∂t+u∂C∂x+v∂C∂y=Dm∂2C∂y2.The corresponding boundary conditions, according to Ali et al. [23–25], are
(5)u=U0,v=0,T=Tw,C=Cwaty=0,(6)u⟶0,T⟶T∞,C⟶C∞asy⟶∞.Here,B0 is the constant magnetic field, ν is the kinematic viscosity, β is the coefficient of volumetric expansion, β∗ is the coefficient of expansion with concentration, σ is the electric conductivity, Q0 is the heat release per unit mass, k is the thermal conductivity, qr is the radiated heat flux, Dm is the coefficient of mass diffusivity, and U0 is the uniform velocity of the fluid, whereas the other symbols retain their ordinary meanings.


### 2.2. Mathematical Formulation

The subsequent usual transformations are introduced to make the nondimensional system of partial differential equations (PDEs) with boundary conditions (BCs).(7)u=U0U,v=VU0,Y=yU0ν,X=xU0ν,η=tU02v,T=T∞+Tw−T∞T¯,C=C∞+Cw−C∞C¯.Replacing the above transformation into equations (1), (2), (3), and (4) and consistent BCs (5) and applying some mathematical manipulation, the following nonlinear PDEs in expressions of dimensionless variables are
(8)∂U∂X+∂V∂Y=0,∂U∂η+U∂U∂X+V∂U∂Y=∂2U∂Y2+GrT¯+GmC¯−MU−NU,∂T¯∂η+U∂T¯∂X+V∂T↼∂Y=1Pr∂2T↼∂Y2+Ec∂U∂Y2+R∂2T↼∂Y2+HT↼,∂C¯∂η+U∂C¯∂X+V∂C¯∂Y=1Sc∂2C¯∂Y2.The BCs are(9)U=1,V=0,T¯=1,C¯=1atY=0,U=0,T¯=0,C¯=0asY⟶∞,where Grashof number Gr=νgβTw−T∞/U03, solutal Grashof number Gm=νgβ∗Cw−C∞/U03, magnetic parameter M=σνB02/ρU02, porosity parameter N=ν2/κU02, Prandtl number Pr=v/α, Eckert number Ec=U02/CpTw−T∞, radiative parameter R=16σ∗T∞3/3ρCpκ0Tw−T∞, Schmidt number Sc=v/Dm, and heat generation parameter H=Q0ν/U02ρCp.


### 2.3. Declaration of Physical Exigent Quantities in Engineering Inquisitiveness

The local skin friction coefficient and Sherwood number are the main physical curiosities in the engineering field. Typically, the shear stress at the plate is denoted by the skin friction coefficient. In compliance with Hossain et al. [26, 27] and Das et al. [28, 29], the nondimensional form of the local shear stress at the plate in the x direction is defined by τxL=μ0∂u/∂y and average shear stress in the x direction τxA=μ0∫∂u/∂ydx. Sherwood number is well known for the local mass transfer rate.The local Sherwood number is signified byShL=−μ0∂C/∂y, while the average Sherwood number is ShA=−μ0∫∂C/∂ydx. The Nusselt number is a well-known measure of the local heat transmission rate. The local Nusselt number is provided by the denoted NuL=−μ0∂T/∂y, and the average Nusselt number is inscribed as NxA=−μ0∫∂T/∂ydx.


## 2.1. Physical Model

Consider the trembling two-dimensional nonlinear MHD free convection flows of a viscoelastic, deformable, and electrically insulating fluid around a vertical plate embedded in porous media. This is happening while the thermal and material buoyancy effects are being encouraged. Figure2 depicts the entire flow mechanism structure with the appropriate coordinates. The x-axis is placed along the porous plate in a vertically upward orientation. Consequently, U is a function of y and t only, i.e., Uy,t. Let u and v represent the velocity components in the x and y directions, respectively. Initially, both the plate and fluid are the same at the temperature. Instantaneously, the plate temperature and species concentration are raised to Tw>T∞ and Cw>C∞, respectively.Figure 2
Physical configuration and coordinate system.Typically, a consistent magnetic fieldB=0,B0,0 is applied to the plate. All assumptions and approximations for the model flow are given as follows:
(i)
The induced magnetic field is insignificant compared to the applied magnetic field(ii)
The Joule heating effects are negligible in the energy equation(iii)
The magnetic Reynolds numbers and effective transversal magnetic field are insignificant(iv)
No applied electric force is expected(v)
The electric field is assumed to be zero(vi)
Viscous dissipation is used to quantify the temperature gradient of fluid particles within the surface of a flat plate(vii)
Flow in the medium is caused by the buoyancy force generated by the temperature gradient between the fluid and the porous plate(viii)
Under the Rosseland [20] approximation of radiative heat flow with Taylor series expansion, the thermal radiation is employed in a way that may be well estimatedUnder the above hypotheses and Boussinesq’s approximation and following Nasrin and Alim [21, 22], the equations of continuity, momentum, energy, and concentration can be put mathematically in the following form.
(1)∂u∂x+∂v∂y=0,(2)∂u∂t+u∂u∂x+v∂u∂y=ν∂2u∂y2+gβT−T∞+gβ∗C−C∞−σB02ρu−νκu,(3)∂T∂t+u∂T∂x+ν∂T∂y=κρCp.∂2T∂y2+νCp∂u∂y2−1ρCp∂qr∂y+Q0ρCPT−T∞,(4)∂C∂t+u∂C∂x+v∂C∂y=Dm∂2C∂y2.The corresponding boundary conditions, according to Ali et al. [23–25], are
(5)u=U0,v=0,T=Tw,C=Cwaty=0,(6)u⟶0,T⟶T∞,C⟶C∞asy⟶∞.Here,B0 is the constant magnetic field, ν is the kinematic viscosity, β is the coefficient of volumetric expansion, β∗ is the coefficient of expansion with concentration, σ is the electric conductivity, Q0 is the heat release per unit mass, k is the thermal conductivity, qr is the radiated heat flux, Dm is the coefficient of mass diffusivity, and U0 is the uniform velocity of the fluid, whereas the other symbols retain their ordinary meanings.


## 2.2. Mathematical Formulation

The subsequent usual transformations are introduced to make the nondimensional system of partial differential equations (PDEs) with boundary conditions (BCs).(7)u=U0U,v=VU0,Y=yU0ν,X=xU0ν,η=tU02v,T=T∞+Tw−T∞T¯,C=C∞+Cw−C∞C¯.Replacing the above transformation into equations (1), (2), (3), and (4) and consistent BCs (5) and applying some mathematical manipulation, the following nonlinear PDEs in expressions of dimensionless variables are
(8)∂U∂X+∂V∂Y=0,∂U∂η+U∂U∂X+V∂U∂Y=∂2U∂Y2+GrT¯+GmC¯−MU−NU,∂T¯∂η+U∂T¯∂X+V∂T↼∂Y=1Pr∂2T↼∂Y2+Ec∂U∂Y2+R∂2T↼∂Y2+HT↼,∂C¯∂η+U∂C¯∂X+V∂C¯∂Y=1Sc∂2C¯∂Y2.The BCs are(9)U=1,V=0,T¯=1,C¯=1atY=0,U=0,T¯=0,C¯=0asY⟶∞,where Grashof number Gr=νgβTw−T∞/U03, solutal Grashof number Gm=νgβ∗Cw−C∞/U03, magnetic parameter M=σνB02/ρU02, porosity parameter N=ν2/κU02, Prandtl number Pr=v/α, Eckert number Ec=U02/CpTw−T∞, radiative parameter R=16σ∗T∞3/3ρCpκ0Tw−T∞, Schmidt number Sc=v/Dm, and heat generation parameter H=Q0ν/U02ρCp.


## 2.3. Declaration of Physical Exigent Quantities in Engineering Inquisitiveness

The local skin friction coefficient and Sherwood number are the main physical curiosities in the engineering field. Typically, the shear stress at the plate is denoted by the skin friction coefficient. In compliance with Hossain et al. [26, 27] and Das et al. [28, 29], the nondimensional form of the local shear stress at the plate in the x direction is defined by τxL=μ0∂u/∂y and average shear stress in the x direction τxA=μ0∫∂u/∂ydx. Sherwood number is well known for the local mass transfer rate.The local Sherwood number is signified byShL=−μ0∂C/∂y, while the average Sherwood number is ShA=−μ0∫∂C/∂ydx. The Nusselt number is a well-known measure of the local heat transmission rate. The local Nusselt number is provided by the denoted NuL=−μ0∂T/∂y, and the average Nusselt number is inscribed as NxA=−μ0∫∂T/∂ydx.


## 3. Computational Technique

The nondimensional governing set of nonlinear PDEs can be explained numerically with relevant BCs with an explicit finite difference method (FDM). The flow zone is divided into a grid of lines corresponding to thex- and y-axes, with the x-axis measured along the plate and the y-axis parallel to the fluid flow, to get the difference equations. Due to this bifurcation, these time- and space-dependent finite difference equations are discretized. In this study the plate of the size Xmax (=100) is measured, that is, X fluctuates between 0 and 100 and supposed Ymax (=25) as corresponding to Y⟶∞, that is, it ranges between 0 and 25. Estimation n′=250 is used to conduct the computation for numerical methods based on random trials of meshes with varying values. Figure 3 shows a visualization of the development of this mesh space. There are m′=250 and n′=250 grid layout in the X and Y conducted separately and occupied as trails Δx=0.40≤x≤100 and ΔY=0.10≤Y≤25 with respect to the minimal time phase, Δη=0.005. The phrase η∞=η∞+Δη determines the value of η∞ throughout each iteration of the loop. The extreme value of η∞ is determined when the successful loop results in no change in the value of the unknown BCs at η=0. Here, U′, T′¯, and C′¯ signify the value of U, T¯, and C¯ at the culmination of a time phase individually.Figure 3
Explicit finite difference system.Exhausting the explicit FDM into the PDEs ((8) and (9)), we get
(10)Ui,j−Ui,j−1ΔX+Vi,j−Vi,j−1ΔY=0,(11)Ui,j′−Ui,jΔη+Ui,jUi,j−Ui−1,jΔX+Vi,jUi,j+1−Ui,jΔY=Ui,j+1−2Ui,j+Ui,j−1ΔY2+GrT¯i,j+GmC¯i,j−MUi,j−NUi,j,(12)Ui,j′=Ui,j+Δη−Ui,jUi,j−Ui−1,jΔX−Vi,jUi,j+1−Ui,jΔY+Ui,j+1−2Ui,j+Ui,j−1ΔY2+GrT¯i,j+GmC¯i,j−MUi,j−NUi,j,(13)T′¯i,j−T¯i,jΔη+Ui,jT¯i,j−T¯i−1,jΔX+Vi,jT¯i,j+1−T¯i,jΔY=1PrT¯i,j+1−2T¯i,j+T¯i,j−1ΔY2+EcUi,j+1−Ui,jΔY2+RT¯i,j+1−2T¯i,j+T¯i,j−1ΔY2+HT¯ij,(14)T′¯i,j=T¯i,j+Δη−UI,JT¯i,j−T¯i−1,jΔX−Vi,jT¯i,j+1−T¯i,jΔY+1PrT¯i,j+1−2T¯i,j+T¯i,j−1ΔY2+EcUi,j+1−Ui,jΔY2+RT¯i,j+1−2T¯i,j+T¯i,j−1ΔY2+HT¯ij,(15)C¯i,j′−C¯i,jΔη+Ui,jC¯i,j−C¯i−1,jΔX+Vi,jC¯i,j+1−C¯i,jΔY=1ScC¯i,j+1−C¯i,j+C¯i,j−1ΔY2C¯i,j′=C¯i,j+Δη−Ui,jC¯i,j−C¯i−1,jΔX−Vi,jC¯i,j+1−C¯i,jΔY+1ScC¯i,j+1−C¯i,j+C¯i,j−1ΔY2,with boundary conditions
(16)Ui,0n=1,Vi,0n=0,T¯i,0n=1,C¯i,0n=1,Ui,Ln=0,T¯i,Ln=0,C¯i,Ln=0, where L⟶∞.Here,i and j stand for the nodes with coordinates, T¯ is the temperature, and C¯ is concentration.


### 3.1. Stability and Convergence Analyses

Stability analysis is an essential tool that gives a good prediction of the convenience of a problem. It is helpful to analyze stability before we go to accurate computations. So, it represents that our numerical solutions can be trusted and is reliable. Ingredients for the entire stability computing process are the minor time phase and the mesh phase values. According to Ali et al. [30], the general expressions of the Fourier expansion for U,θ,ϕ an arbitrary time, say, η=0, are eiα¯x¯ and eiβ¯y¯, where i=−1.Then,(17)U:Ψηeiα−x−⋅eiβ−y−,T¯:θηeiα¯x¯eiβ¯y¯,C¯:φηeiα¯x¯eiβ¯y¯.And after the time step, we get(18)U:Ψ′ηeiα−x−⋅eiβ−y−,T−:θ′ηeiα−x−⋅eiβ−y−,C−:φ′ηeiα−x−⋅eiβ−y−.Substituting equations (17) and (18) into equations (10)–(16), the following equations are obtained upon simplification.
(19)Ψ′−ΨΔη+UΨ1−e−iα−ΔXΔX+VΨeiβ−ΔY−1ΔY=2ΨΔY⋅cosβ−1ΔY2+Grθ′+Gmϕ′−MΨ−NΨ,(20)Ψ′−Ψ+ΔηΔXUΨ1−e−iα−ΔX+ΔηΔYVΨeiβ−ΔY−1=2ΔηΔY2ΨΔYcosβ−1+GrΔηθ′+GmΔηϕ′−MΨ−NΨ,⇒Ψ′=Ψ1−ΔηΔXU1−e−α−iΔX−ΔηΔYVeiβ−ΔY−1+2ΔηΔY2ΔYcosβ−1−M+NΔη+GrΔηθ′+GmΔηϕ′,Ψ′=AΨ+Bθ′+Cϕ′,where
(21)A=1−ΔηΔxU1−e−iα−ΔX−ΔηΔYVeiβ−ΔY−1+2ΔηΔY2cosβ⋅ΔY−1−M+N⋅Δη,B=GrΔη,C=GmΔη,(22)θ′η−θηΔη+Uθη1−e−iα−ΔXΔX+Vθη·e−iβ̄ΔY−1ΔY=1Pr2θηcosβ·ΔY−1ΔY2+EcUΨηeiβ−ΔY−1ΔY2+R2θηcosβ·ΔY−1ΔY2+Hθη,(23)θ′η=θη1−ΔηΔxU1−e−iα−ΔX−ΔηΔYVeiβ−ΔY−1+1Pr2ΔηΔY2cosβ⋅ΔY−1+EcU⋅ΔηΔY2Ψηeiβ−ΔY−1+RΔηcosβ⋅ΔY−1ΔY2+HΔη,θ′η=Gθ+KΨ,where
(24)G=1−ΔηΔXU1−e−iα¯ΔX−ΔηΔYVeiβ−ΔY−1+1Pr2Δηcosβ⋅ΔY−1ΔY2,K=EcU⋅ΔηΔY2eiβ−ΔY−1,P=RΔηcosβ⋅ΔY−1ΔY2+HΔη,(25)ϕ′η−ϕηΔη+Uϕη1−e−iα¯ΔXΔX+Vϕη⋅eiβ¯ΔY−1ΔY=1Sc⋅2ϕηcosβ⋅ΔY−1ΔY2,ϕ′η=ϕη1−ΔηΔXU1−e−iα−ΔX−ΔηΔYUeiβ−ΔY−1+1Sc⋅2ΔηΔY2cosβ⋅ΔY−1,(26)ϕ′η=Qϕη,where
(27)θ=1−ΔηΔXU1−e−iα−ΔX−ΔηΔYVeiβ−ΔY−1+1Sc2ΔηΔY2cosβ⋅ΔY−1.Equations (19), (23), and (25) can be written as
(28)Ψ′=AΨ+BGθ+KΨ+CQϕ=A+KΨ+BGθ+CQϕ,(29)Ψ′=A1Ψ+B1θ+C1ϕ,where
(30)A1=A+K,B1=BG,C1=CQ,(31)θ′=Gθ+KΨ,and ϕ′=Qϕ.
(32)ϕ′=0⋅Ψ+0⋅θ+Qϕ.Equations (28), (31), and (32) can be written as
(33)ψ′θ′ϕ′=A1B1C1KG000Qψθϕ,i.e.,η′=Tη,where
(34)T=A1B1C1KG000Q,η=ψθϕ.To obtain the stability condition, it is necessary to realize the eigenvalues of the augmentation matrixT.For this reason, letB1⟶0,C1⟶0,K⟶0.So, matrixT can be expressed as
(35)T=A1000G000Q.Thus, the eigenvalues areλ1=A1,λ2=G,λ3=Q.For stability, each eigenvalueλ1,λ2,λ3 must not exceed in modulus.Hence, the stability condition is(36)A1≤1,G≤1,Q≤1.Let,a=uΔη/ΔX,b=Δη⋅V,c=Δη/ΔY2.Hence,A=1−a1−e−iα−ΔX−beiβ−ΔY−1+2ccosβΔY−1−M+NΔη,
(37)K=EcUcCiβ−ΔY−1.The coefficients ofa,b,c are real and nonnegative. So, the maximum modulus of A1,G,Q occurs when α¯ΔX=mπ and β¯ΔY=nπ where m and n are integers, and hence, A1,G,Q are real. When m and n are odd integers, the amounts of A1,G,Q are maximum. So we get
(38)A1=A+K=1−2a−2b−4c−2cEc.To satisfyA1≤1,G≤1,Q≤1, the greatest negative allowable values are
(39)A1=−1,G=−1,Q=−1.So, the converging restrictions for the FDM solution schemes are(40)1−2a−2b−4c−2cEc≤−1,2a+b+2c+cEc≤2,a+b+2c+cEc≤1,UΔηΔX+VΔηΔY+2ΔηΔY2+ΔηΔY2Ec≤1.Likewise, the second condition(41)UΔηΔX+VΔηΔY+1Pr2ΔηΔY2+R2ΔηΔY2≤1.The third condition(42)UΔηΔX+VΔηΔY+1Sc2ΔηΔY2≤1.The values∆η=0.005, Δx=0.4, and Δy=0.1, along with the initial circumstances that determine the convergence requirements of the method, are Pr≥1 and Sc≥0.22 with stable constraints materially.


### 3.2. Questing for Suitable Mesh

After repeated experiments with various mesh sizes, the ideal mesh size is strategically significant for the flow model inside the boundary layer flow. The acquaintance with the graphs’ smoothness and convergence leads us to conclude that the mesh sizen=250 can be measured for numerical designs.


### 3.3. Questing Steady State

The uninterrupted testing of various parameter values attains the steady-state scenario for the flow model equations. The steady-state controls for magnetic field on shear stressY=0.3, Y=0.5, Y=1.5, and Y=2.5 are presented in Figure 4, with M=9.5, Gr=3.5, Gm=4.5, Pr=1, Ec=0.1, Sc=0.22, N=1, R=1, and H=0.3. Observations show that the consequence of the magnetic field on shear stress expresses minor fluctuation at Y=0.3 and Y=0.5 while Y=1 and Y=2 exhibit almost ineffectual variations as an acquisitive steady state for the revealing fields. The steady-state performance is exposed in (Table 1).Figure 4
Validations of steady state for magnetic field on shear stress atY=0.3, Y=0.5, Y=1.5, and Y=2.5.Table 1
Validations of steady state for magnetic field on shear stress for Figure4 at Y=0.3, Y=0.5, Y=1.5, and Y=2.5.
Average skin frictionParameterUnsteadyFluctuation (unstable) %Steady stateFluctuation (unstable) %Y=0.3Y=0.5Y=1.5Y=2.5M=9.50.00008287112723340.00008752836583880.0004657240.00008753045403860.00008753045403860M=10.50.00008671601972310.00008752994472168.13925E−050.00011325149534970.00011325149534970M=11.50.00014564763114320.00014690928219350.0001261650.00014690999523740.00014690999523740


### 3.4. Verifications of the Consequence of Parameter

The magnitudes of the Grashof number on concentration profile are given different values for radiation and in the absence of radiation which assures the accuracy of the effect of the parameter. This experimentation is demonstrated through Figures5(a) and 5(b), with M=9.5, Gr=3.5, Gm=4.5, Pr=1, Ec=0.1, Sc=0.22, N=1, R=1, H=0.3 and M=9.5, Gr=3.5, Gm=4.5, Pr=1, Ec=0.1, Sc=0.22, N=1, R=0, H=0.3. In the presence of radiation/absence of radiation, curve-to-curve fluctuation on concentration profiles is shown in (Table 2).Figure 5
Justifications of the parameter (radiation) on concentration profiles for Grashof number (a) in the presence of radiation and (b) in the absence of radiation.
(a)(b)Table 2
Justifications of the parameter (radiation) for Figures5(a) and 5(b) at Y=1.
ConcentrationRadiationAbsence of radiationParameterNumerical valueIncrease/decrease in %Numerical valueIncrease/decrease in %Effect of variation in %Gr=3.50.31471711400.5337356329Gr=4.50.2613359995.33811152% decrease0.49168854954.20470834% decrease1.13340318 (%) decreaseGr=5.50.20020052796.11354709% decrease0.42426365616.74248934% decrease0.62894225% decrease


### 3.5. Verifications of Codes with the Reported Results

The published result of Deepthi et al. [31] and the present results show excellent potential. The results of the pictorial matching are shown in Figures 6(a) and 6(b). The effects of heat generation on the temperature inside the boundary layer are carried out in this comparison.Figure 6
The effects of heat generation on temperature profiles in (a) Deepthi et al. [29] and (b) the present study.
(a)(b)


## 3.1. Stability and Convergence Analyses

Stability analysis is an essential tool that gives a good prediction of the convenience of a problem. It is helpful to analyze stability before we go to accurate computations. So, it represents that our numerical solutions can be trusted and is reliable. Ingredients for the entire stability computing process are the minor time phase and the mesh phase values. According to Ali et al. [30], the general expressions of the Fourier expansion for U,θ,ϕ an arbitrary time, say, η=0, are eiα¯x¯ and eiβ¯y¯, where i=−1.Then,(17)U:Ψηeiα−x−⋅eiβ−y−,T¯:θηeiα¯x¯eiβ¯y¯,C¯:φηeiα¯x¯eiβ¯y¯.And after the time step, we get(18)U:Ψ′ηeiα−x−⋅eiβ−y−,T−:θ′ηeiα−x−⋅eiβ−y−,C−:φ′ηeiα−x−⋅eiβ−y−.Substituting equations (17) and (18) into equations (10)–(16), the following equations are obtained upon simplification.
(19)Ψ′−ΨΔη+UΨ1−e−iα−ΔXΔX+VΨeiβ−ΔY−1ΔY=2ΨΔY⋅cosβ−1ΔY2+Grθ′+Gmϕ′−MΨ−NΨ,(20)Ψ′−Ψ+ΔηΔXUΨ1−e−iα−ΔX+ΔηΔYVΨeiβ−ΔY−1=2ΔηΔY2ΨΔYcosβ−1+GrΔηθ′+GmΔηϕ′−MΨ−NΨ,⇒Ψ′=Ψ1−ΔηΔXU1−e−α−iΔX−ΔηΔYVeiβ−ΔY−1+2ΔηΔY2ΔYcosβ−1−M+NΔη+GrΔηθ′+GmΔηϕ′,Ψ′=AΨ+Bθ′+Cϕ′,where
(21)A=1−ΔηΔxU1−e−iα−ΔX−ΔηΔYVeiβ−ΔY−1+2ΔηΔY2cosβ⋅ΔY−1−M+N⋅Δη,B=GrΔη,C=GmΔη,(22)θ′η−θηΔη+Uθη1−e−iα−ΔXΔX+Vθη·e−iβ̄ΔY−1ΔY=1Pr2θηcosβ·ΔY−1ΔY2+EcUΨηeiβ−ΔY−1ΔY2+R2θηcosβ·ΔY−1ΔY2+Hθη,(23)θ′η=θη1−ΔηΔxU1−e−iα−ΔX−ΔηΔYVeiβ−ΔY−1+1Pr2ΔηΔY2cosβ⋅ΔY−1+EcU⋅ΔηΔY2Ψηeiβ−ΔY−1+RΔηcosβ⋅ΔY−1ΔY2+HΔη,θ′η=Gθ+KΨ,where
(24)G=1−ΔηΔXU1−e−iα¯ΔX−ΔηΔYVeiβ−ΔY−1+1Pr2Δηcosβ⋅ΔY−1ΔY2,K=EcU⋅ΔηΔY2eiβ−ΔY−1,P=RΔηcosβ⋅ΔY−1ΔY2+HΔη,(25)ϕ′η−ϕηΔη+Uϕη1−e−iα¯ΔXΔX+Vϕη⋅eiβ¯ΔY−1ΔY=1Sc⋅2ϕηcosβ⋅ΔY−1ΔY2,ϕ′η=ϕη1−ΔηΔXU1−e−iα−ΔX−ΔηΔYUeiβ−ΔY−1+1Sc⋅2ΔηΔY2cosβ⋅ΔY−1,(26)ϕ′η=Qϕη,where
(27)θ=1−ΔηΔXU1−e−iα−ΔX−ΔηΔYVeiβ−ΔY−1+1Sc2ΔηΔY2cosβ⋅ΔY−1.Equations (19), (23), and (25) can be written as
(28)Ψ′=AΨ+BGθ+KΨ+CQϕ=A+KΨ+BGθ+CQϕ,(29)Ψ′=A1Ψ+B1θ+C1ϕ,where
(30)A1=A+K,B1=BG,C1=CQ,(31)θ′=Gθ+KΨ,and ϕ′=Qϕ.
(32)ϕ′=0⋅Ψ+0⋅θ+Qϕ.Equations (28), (31), and (32) can be written as
(33)ψ′θ′ϕ′=A1B1C1KG000Qψθϕ,i.e.,η′=Tη,where
(34)T=A1B1C1KG000Q,η=ψθϕ.To obtain the stability condition, it is necessary to realize the eigenvalues of the augmentation matrixT.For this reason, letB1⟶0,C1⟶0,K⟶0.So, matrixT can be expressed as
(35)T=A1000G000Q.Thus, the eigenvalues areλ1=A1,λ2=G,λ3=Q.For stability, each eigenvalueλ1,λ2,λ3 must not exceed in modulus.Hence, the stability condition is(36)A1≤1,G≤1,Q≤1.Let,a=uΔη/ΔX,b=Δη⋅V,c=Δη/ΔY2.Hence,A=1−a1−e−iα−ΔX−beiβ−ΔY−1+2ccosβΔY−1−M+NΔη,
(37)K=EcUcCiβ−ΔY−1.The coefficients ofa,b,c are real and nonnegative. So, the maximum modulus of A1,G,Q occurs when α¯ΔX=mπ and β¯ΔY=nπ where m and n are integers, and hence, A1,G,Q are real. When m and n are odd integers, the amounts of A1,G,Q are maximum. So we get
(38)A1=A+K=1−2a−2b−4c−2cEc.To satisfyA1≤1,G≤1,Q≤1, the greatest negative allowable values are
(39)A1=−1,G=−1,Q=−1.So, the converging restrictions for the FDM solution schemes are(40)1−2a−2b−4c−2cEc≤−1,2a+b+2c+cEc≤2,a+b+2c+cEc≤1,UΔηΔX+VΔηΔY+2ΔηΔY2+ΔηΔY2Ec≤1.Likewise, the second condition(41)UΔηΔX+VΔηΔY+1Pr2ΔηΔY2+R2ΔηΔY2≤1.The third condition(42)UΔηΔX+VΔηΔY+1Sc2ΔηΔY2≤1.The values∆η=0.005, Δx=0.4, and Δy=0.1, along with the initial circumstances that determine the convergence requirements of the method, are Pr≥1 and Sc≥0.22 with stable constraints materially.


## 3.2. Questing for Suitable Mesh

After repeated experiments with various mesh sizes, the ideal mesh size is strategically significant for the flow model inside the boundary layer flow. The acquaintance with the graphs’ smoothness and convergence leads us to conclude that the mesh sizen=250 can be measured for numerical designs.


## 3.3. Questing Steady State

The uninterrupted testing of various parameter values attains the steady-state scenario for the flow model equations. The steady-state controls for magnetic field on shear stressY=0.3, Y=0.5, Y=1.5, and Y=2.5 are presented in Figure 4, with M=9.5, Gr=3.5, Gm=4.5, Pr=1, Ec=0.1, Sc=0.22, N=1, R=1, and H=0.3. Observations show that the consequence of the magnetic field on shear stress expresses minor fluctuation at Y=0.3 and Y=0.5 while Y=1 and Y=2 exhibit almost ineffectual variations as an acquisitive steady state for the revealing fields. The steady-state performance is exposed in (Table 1).Figure 4
Validations of steady state for magnetic field on shear stress atY=0.3, Y=0.5, Y=1.5, and Y=2.5.Table 1
Validations of steady state for magnetic field on shear stress for Figure4 at Y=0.3, Y=0.5, Y=1.5, and Y=2.5.
Average skin frictionParameterUnsteadyFluctuation (unstable) %Steady stateFluctuation (unstable) %Y=0.3Y=0.5Y=1.5Y=2.5M=9.50.00008287112723340.00008752836583880.0004657240.00008753045403860.00008753045403860M=10.50.00008671601972310.00008752994472168.13925E−050.00011325149534970.00011325149534970M=11.50.00014564763114320.00014690928219350.0001261650.00014690999523740.00014690999523740


## 3.4. Verifications of the Consequence of Parameter

The magnitudes of the Grashof number on concentration profile are given different values for radiation and in the absence of radiation which assures the accuracy of the effect of the parameter. This experimentation is demonstrated through Figures5(a) and 5(b), with M=9.5, Gr=3.5, Gm=4.5, Pr=1, Ec=0.1, Sc=0.22, N=1, R=1, H=0.3 and M=9.5, Gr=3.5, Gm=4.5, Pr=1, Ec=0.1, Sc=0.22, N=1, R=0, H=0.3. In the presence of radiation/absence of radiation, curve-to-curve fluctuation on concentration profiles is shown in (Table 2).Figure 5
Justifications of the parameter (radiation) on concentration profiles for Grashof number (a) in the presence of radiation and (b) in the absence of radiation.
(a)(b)Table 2
Justifications of the parameter (radiation) for Figures5(a) and 5(b) at Y=1.
ConcentrationRadiationAbsence of radiationParameterNumerical valueIncrease/decrease in %Numerical valueIncrease/decrease in %Effect of variation in %Gr=3.50.31471711400.5337356329Gr=4.50.2613359995.33811152% decrease0.49168854954.20470834% decrease1.13340318 (%) decreaseGr=5.50.20020052796.11354709% decrease0.42426365616.74248934% decrease0.62894225% decrease


## 3.5. Verifications of Codes with the Reported Results

The published result of Deepthi et al. [31] and the present results show excellent potential. The results of the pictorial matching are shown in Figures 6(a) and 6(b). The effects of heat generation on the temperature inside the boundary layer are carried out in this comparison.Figure 6
The effects of heat generation on temperature profiles in (a) Deepthi et al. [29] and (b) the present study.
(a)(b)


## 4. Results with Explications

To attain the physical implication of the problems, we have plotted flow, thermal, material, Sherwood number, Nusselt number, and skin friction for distinct values of physical fields like magnetic parameter (M), Prandtl number (Pr), Eckert number (Ec), porosity (N), and Schmidt number (Sc). The following default parameter values are assumed for computations in the present study: M=9.5, Pr=1, Ec=0.1, Sc=0.22, and N=1. To be realistic, the numerical values of the Prandtl number are chosen as 1, 1.38, and 7.2, which correspond to salt water, ammonia, and seawater (at 20°C), respectively. Also, the typical values of Schmidt number are preferred for hydrogen (0.22), ethyl ether (1.66), and toluene (1.84) at 25°C and 1 atm, which denotes the diffusion of the most prevalent chemical species. The numerical values for the remaining parameters are determined at random. Grashof number (Gr=0.5), solutal Grashof number (Gm=4.5), radiation (R=1), and heat generation (H=0.3) are kept constant.


### 4.1. Influence of Magnetic Field

Figures7(a)–7(c) show flow, thermal, and concentration profile trends for various values of the magnetic parameters while maintaining the remaining parameters constant. Figure 7(a) demonstrates that the velocity decreased for increasing values of M (Table 3). It is a known fact that a transverse magnetic field in proximity to an electrically conducting fluid will cause an increase in the body force known as the Lorentz force. This force tends to resist the flow of the fluid and slow down its motion in the BL region. In consequence, tangential velocity falls. In Figure 7(b), it is apparent that the temperature of the fluid rises in 0≤Y≤1.2approx. and then falls as the value of M increases (Table 4). Figure 7(c) depicts that the concentration profiles decrease for enhancement of M (Table 5), indicating that it slows the rate of heat convection in the flow; hence, the flow temperature seems to rise, and the concentration gradient decreases. Figures 7(d) and 7(g) show that local and average Nusselt number declines for growing values of M. Figures 7(e) and 7(h) show that the local and average Sherwood number upsurges for rising values of M.Figure 7
Influence of magnetic field on (a) flow, (b) thermal, (c) material fields, (d)NuL, (e) ShL, (f) τxL, (g) NuA, (h) ShA (i) τxA, (j) flow fields, (k) τxL, and (l) τxA in 3D plot with Gr=3.5, Gm=4.5, Pr=1, Ec=0.1, Sc=0.22, N=1, R=1, and H=0.3.
(a)(b)(c)(d)(e)(f)(g)(h)(i)(j)(k)(l)Table 3
Variability of curve whenY=1 in Figure 7(a).
VelocityValue of parameterNumerical valueIncrease/decrease in %M=9.50.2696958184M=10.50.25263792281.70578956% decreaseM=11.50.23844483491.4193087% decreaseTable 4
Variability of the curve whenY=1 in Figure 7(b).
TemperatureValue of parameterNumerical valueIncrease/decrease in %M=9.50.40298757M=10.50.38619664311.10608041% decreaseM=11.50.3751358391.10608041% decreaseTable 5
Variability of the curve whenY=1 in Figure 7(c).
ConcentrationValue of parameterNumerical valueIncrease/decrease in %M=9.50.3147171140M=10.50.34433361892.96165049% increaseM=11.50.36980471012.54710912% increaseFigures7(f) and 7(i) exhibit that the local and average shear stress is the upsurge for increasing values of M, since an enhancement in portion causes an accumulation in the fluids’ thermal conductivity. In the 3D graph, Figure 7(j) exemplifies that the flow profiles decrease with growing values of the magnetic fields. Figures 7(k) and 7(l) show that for rising values of M, the local and average shear stress in 3D graphs increases.


### 4.2. Influence of Prandtl Number

Figure8(a) illustrates that the velocity profiles increase 0≤Y≤.5 and then decline for increasing values of Pr. The reason is that lesser value Pr is equivalent to the rise in the thermal conductivity of the fluid, and it also falls very quickly for higher values of Pr which creates high viscosity in the boundary layer. Figure 8(b) demonstrates that the thermal profile increases in the region 0≤Y≤.8 and then falls for growing values of Pr. Temperature gradients for influenced acceptable quantities of Pr=1, 1.38, and 7.2 which are significant to salt water, ammonia, and seawater (at 20°C), respectively, and they associate with salt water and ammonia. Also, we see that increasing Pr lessens thermal diffusivity; heat transmission capacity substantially reduces as thermal diffusivity lowers, resulting in a decline in the temperature distribution. The increased Pr is proportional to reduced thermal conductivity. As a result, concentration profile rise shown in Figure 8(c). Figures 8(d) and 8(g) reveal that local and average Nusselt number declines for rising values of Pr (Tables 6 and 7), which creates the reduction in thickness of the thermal boundary layer, as well as an overall lessening of the average temperature drop inside the boundary layer.Figure 8
Influence of Prandtl number on (a) flow, (b) thermal, (c) material fields, (d)NuL, (e) ShL, (f) τxL, (g) NuA, (h) ShA (i) τxA, and (j) ShL and ShA comparatively in a 3D plot with M=9.5, Gr=3.5, Gm=4.5, Ec=0.1, Sc=0.22, N=1, R=1, and H=0.3.
(a)(b)(c)(d)(e)(f)(g)(h)(i)(j)Table 6
Variability of the curve whenY=1 in Figure 8(c).
Local Nusselt numberValue of parameterNumerical valueIncrease/decrease in %Pr=1−0.0149532239884138Pr=1.38−0.01763824746012690.706585124% decreasePr=7.2−0.02874221466481690.190789815% decreaseTable 7
Variability of the curve whenY=1 in Figure 8(g).
Average Nusselt numberValue of parameterNumerical valueIncrease/decrease in %Pr=1−0.0149532239884138Pr=1.38−0.01763824746012690.706585124% decreasePr=7.2−0.02874221466481690.190789815% decreaseThe local and average Sherwood number losses (0≤Y≤2.5) and then upsurges for mounting values of Pr, which are shown in Figures 8(e) and 8(h) (Tables 8 and 9). Heat distributes gradually, and the thermal boundary layer is lower for a more significant Prandtl number, whereas the reverse scenario is shown at lower Pr values. Figures 8(f) and 8(i) reveal that local and average shear stress decreases for upgrading Pr values. Enlarging Pr, thermal diffusivity diminishes during fluid velocity augmentations. A 3D graph comparison of the local and average Sherwood number of Pr is shown in Figure 8(j).Table 8
Variability of the curve whenY=0.5 in Figure 8(e).
Local Sherwood numberValue of parameterNumerical valueIncrease/decrease in %Pr=10.0000157634385687Pr=1.380.00006680995284110.013433293% increasePr=7.20.00009517795115240.000487423% increaseTable 9
Variability of the curve whenY=0.5 in Figure 8(h).
Average Sherwood numberValue of parameterNumerical valueIncrease/decrease in %Pr=1−0.0149532239884138Pr=1.38−0.01763824746012690.706585124% decreasePr=7.2−0.02874221466481690.190789815% decrease


### 4.3. Influence of Porosity

Figure9(a) illustrates that the flow profiles rise for growing values of the porosity parameter N. It arises due to the rise permeability of the medium, and the fluid motion is subsequently accelerated. Figure 9(b) illustrates the effect of porosity parameters N on temperature, which leads to the rise in temperature profiles (η<14), and the temperature profiles decrease due to enhancement of N. Also, it is seen that the porosity parameter strongly affects velocity profiles compared to temperature profiles. Figure 9(c) shows that the material profiles rise with rising values of N. Because the porous medium slows down the motion of the fluid, as a result, mass transfer rises. Figures 9(d) and 9(g) reveal that local and average Nusselt number declines for increasing values of N. Figures 9(e) and 9(h) expose local and average Sherwood number reductions for growing values of N. Figures 9(f) and 9(i) reveal local and average shear stress decrease for rising values of N (Tables 10 and 11). It is true that when the permeability parameter’s values increase, the tightness of the porous medium’s resistance to flow diminishes, which in turn causes the fluid’s velocity to increase. Figure 9(j) exemplifies that the flow profiles upsurge with rising values of the porosity parameter in the 3D graph.Figure 9
Influence of porosity on (a) flow, (b) thermal, (c) material fields, (d)NuL, (e) ShL, (f) τxL, (g) NuA, (h) ShA (i) τxA, and (j) flow fields in the 3D plot with M=9.5, Gr=3.5, Gm=4.5, Pr=1, Ec=0.1, Sc=0.22, R=1, and H=0.3.
(a)(b)(c)(d)(e)(f)(g)(h)(i)(j)Table 10
Variability of the curve whenY=1 in Figure 9(f).
Local skin frictionValue of parameterNumerical valueIncrease/decrease in %N=10.0043904874473810N=20.00483940448611970.044891704% decreaseN=30.00538563309237360.054622861% decreaseTable 11
Variability of the curve whenY=1 in Figure 9(i).
Average skin frictionValue of parameterNumerical valueIncrease/decrease in %N=10.0043904874473810N=20.00483940448611970.044891704% decreaseN=30.00538563309237360.054622861% decrease


### 4.4. Influence of Eckert Number

Figure10(a) shows that the flow profiles decrease (η>2) owing to rises in Ec reducing the driving force to heat transfer to the kinetic energy of the flow; that is, there is no effect of Ec on velocity profiles near the plate (Table 12). Figure 10(b) exposes thermal profile declines for rising values of Ec. The accumulation of heat energy in liquid happens due to the fact that frictional heating arises. Consequently, the fluid flow rate increases, and the heat transmission rate falls. As a result, the thermal profiles decrease. The effects of Ec on material profiles are observed in Figure 10(c). The concentration profiles increase with growing values of Ec. It allows energy to be stored in the fluid, which creates frictional heating; for this reason, the concentration profiles increase. Figures 10(d) and 10(g) reveal that local and average Nusselt numbers rise for increasing values of Ec. In Figures 10(e), 10(h), 10(f), and 10(i), the increase in the viscous dissipation term (Eckert number (Ec)) is perceived to decline the local and average Sherwood number as well as shear stress. Materially, the rise in Eckert number causes the heat energy inside the fluid layer to increase. Due to the drop in flow regimes, the heat energy formed by the Eckert number is not conserved, and there is no assurance of frictional heating.Figure 10
Influence of Eckert number on (a) flow, (b) thermal, (c) material fields, (d)NuL, (e) ShL, (f) τxL, (g) NuA, (h) ShA, and (i) τxA with M=9.5, Gr=3.5, Gm=4.5, Pr=1, Sc=0.22, N=1, R=1, and H=0.3.
(a)(b)(c)(d)(e)(f)(g)(h)(i)Table 12
Variability of the curve whenY=1 in Figure 10(a).
VelocityValue of parameterNumerical valueIncrease/decrease in %Ec=0.10.269695818Ec=0.20.19307398838.3109152% decreaseEc=0.30.1365558828.25905385% decrease


### 4.5. Influence of Schmidt Number

To be accurate, the typical values of Schmidt number (Sc) are preferred for hydrogen (Sc=0.22), ethyl ether (Sc=1.66), and toluene (Sc=1.84) at 25°C and 1 atm, which denotes chemical species dispersing in the most prevalent aspect. Figures 11(a)–11(c) exhibit characteristic behavior of flow, thermal, and material profiles for rising quantities of Sc while maintaining the other parameters constant. It can be shown in Figures 11(a) and 11(b) that the flow and thermal fields increase owing to the enhancement of Sc. Schmidt number measures the relative efficiency of momentum and mass transmission through diffusions in the hydrodynamic and material boundary layer (BL). Figure 11(c) depicts the influence of the Schmidt number on concentration profiles and see that the material profile upsurges with rising values of Sc (Table 13). It is a fact that the rising effect tends to decline the boundary layer and concentrations buoyancy effects. In this case, the flow of fluid diminishes, and the rate of concentration increases. Figures 11(d) and 11(g) reveal that local and average Nusselt numbers rise significantly for growing Sc values. Figures 11(e) and 11(h) illustrate the local and average Sherwood number rise for increasing values of Sc. Otherwise, the rate of mass transmission at the plate drops for high mass diffusivity. Figures 11(f) and 11(i) reveal local and average shear stress decreases for upgrading Sc values. The impact represents the ratio between momentum and mass diffusivity. The flow and material profiles are seen to decline as Sc rises. As a result, the material buoyancy effects lessen, which lowers the fluid velocity. Figure 11(j) shows the consequence of the Schmidt number on thermal profiles in 3D graphs.Figure 11
Influence of Schmidt number on (a) flow, (b) thermal, (c) material fields, (d)NuL, (e) ShL, (f) τxL, (g) NuA, (h) ShA (i) τxA, and (j) thermal fields in the 3D plot with M=9.5, Gr=3.5, Gm=4.5, Pr=1, Ec=0.1, N=1, R=1, and H=0.3.
(a)(b)(c)(d)(e)(f)(g)(h)(i)(j)Table 13
Variability of the curve whenY=1 in Figure 11(c).
ConcentrationValue of parameterNumerical valueIncrease/decrease in %Sc=0.220.0111813480Sc=1.660.01415115871.649894833% increaseSc=1.840.240881905015.74519072% increase


### 4.6. Comparison

The local Nusselt number acts on the enlargement of heat flow at any point, and the average Nusselt number is significant at all points. So, the average Nusselt number is equal to the local Nusselt number in these circumstances for the consequences of several factors. The local Sherwood number acts on the development of mass transfer at any point, and the average Sherwood number is substantial at all points. So, the average Sherwood number is equivalent to the local Sherwood number for the effects of numerous parameters. In these cases, local and average skin friction is shown to have identical behaviors for the impacts of rising parameters. The consequent table (Table14) presents a qualitative comparison between the described work and several researchers from the literature studies.Table 14
Comparison of the current findings qualitatively to other published results.
(a)ParametersRao et al. [32]Ali et al. [33]Present resultsUTCUTCUTCMLes——LesLesEscLesPrLesLesEscLesLesEscLesLes/EscEscNEsc——EscEscEscEscEsc(b)ParametersSambath et al. [34]Present resultsτxLShLNuLτxLShLNuLMLesEscEscLesPrLesLesLesLesLes/EscLesEcLesLesEscNEscEscEscEscLesScLesEscLesEscEscLes means lessening; Esc means escalation.This research is aimed at finding the numerical solution for an MHD heat and mass transfer flow that is viscous, electrically conducting, incompressible, radiating, and flows over a permeable plate. The BL heat absorption, radiation, chemical reaction, and Soret effect-based flow of MHD fluid through a semi-infinite vertical permeable moving plate were examined by Rao et al. [32]. Ali et al. [33] investigated MHD heat and mass transfer portents using magnetic field stimulation along a vertical plate. The effect of Sherwood number, Nusselt number, and skin friction on MHD radiative mass transfer flow from an inclined vertical plate with a heat source and sink was investigated by Sambath et al. [34].


## 4.1. Influence of Magnetic Field

Figures7(a)–7(c) show flow, thermal, and concentration profile trends for various values of the magnetic parameters while maintaining the remaining parameters constant. Figure 7(a) demonstrates that the velocity decreased for increasing values of M (Table 3). It is a known fact that a transverse magnetic field in proximity to an electrically conducting fluid will cause an increase in the body force known as the Lorentz force. This force tends to resist the flow of the fluid and slow down its motion in the BL region. In consequence, tangential velocity falls. In Figure 7(b), it is apparent that the temperature of the fluid rises in 0≤Y≤1.2approx. and then falls as the value of M increases (Table 4). Figure 7(c) depicts that the concentration profiles decrease for enhancement of M (Table 5), indicating that it slows the rate of heat convection in the flow; hence, the flow temperature seems to rise, and the concentration gradient decreases. Figures 7(d) and 7(g) show that local and average Nusselt number declines for growing values of M. Figures 7(e) and 7(h) show that the local and average Sherwood number upsurges for rising values of M.Figure 7
Influence of magnetic field on (a) flow, (b) thermal, (c) material fields, (d)NuL, (e) ShL, (f) τxL, (g) NuA, (h) ShA (i) τxA, (j) flow fields, (k) τxL, and (l) τxA in 3D plot with Gr=3.5, Gm=4.5, Pr=1, Ec=0.1, Sc=0.22, N=1, R=1, and H=0.3.
(a)(b)(c)(d)(e)(f)(g)(h)(i)(j)(k)(l)Table 3
Variability of curve whenY=1 in Figure 7(a).
VelocityValue of parameterNumerical valueIncrease/decrease in %M=9.50.2696958184M=10.50.25263792281.70578956% decreaseM=11.50.23844483491.4193087% decreaseTable 4
Variability of the curve whenY=1 in Figure 7(b).
TemperatureValue of parameterNumerical valueIncrease/decrease in %M=9.50.40298757M=10.50.38619664311.10608041% decreaseM=11.50.3751358391.10608041% decreaseTable 5
Variability of the curve whenY=1 in Figure 7(c).
ConcentrationValue of parameterNumerical valueIncrease/decrease in %M=9.50.3147171140M=10.50.34433361892.96165049% increaseM=11.50.36980471012.54710912% increaseFigures7(f) and 7(i) exhibit that the local and average shear stress is the upsurge for increasing values of M, since an enhancement in portion causes an accumulation in the fluids’ thermal conductivity. In the 3D graph, Figure 7(j) exemplifies that the flow profiles decrease with growing values of the magnetic fields. Figures 7(k) and 7(l) show that for rising values of M, the local and average shear stress in 3D graphs increases.


## 4.2. Influence of Prandtl Number

Figure8(a) illustrates that the velocity profiles increase 0≤Y≤.5 and then decline for increasing values of Pr. The reason is that lesser value Pr is equivalent to the rise in the thermal conductivity of the fluid, and it also falls very quickly for higher values of Pr which creates high viscosity in the boundary layer. Figure 8(b) demonstrates that the thermal profile increases in the region 0≤Y≤.8 and then falls for growing values of Pr. Temperature gradients for influenced acceptable quantities of Pr=1, 1.38, and 7.2 which are significant to salt water, ammonia, and seawater (at 20°C), respectively, and they associate with salt water and ammonia. Also, we see that increasing Pr lessens thermal diffusivity; heat transmission capacity substantially reduces as thermal diffusivity lowers, resulting in a decline in the temperature distribution. The increased Pr is proportional to reduced thermal conductivity. As a result, concentration profile rise shown in Figure 8(c). Figures 8(d) and 8(g) reveal that local and average Nusselt number declines for rising values of Pr (Tables 6 and 7), which creates the reduction in thickness of the thermal boundary layer, as well as an overall lessening of the average temperature drop inside the boundary layer.Figure 8
Influence of Prandtl number on (a) flow, (b) thermal, (c) material fields, (d)NuL, (e) ShL, (f) τxL, (g) NuA, (h) ShA (i) τxA, and (j) ShL and ShA comparatively in a 3D plot with M=9.5, Gr=3.5, Gm=4.5, Ec=0.1, Sc=0.22, N=1, R=1, and H=0.3.
(a)(b)(c)(d)(e)(f)(g)(h)(i)(j)Table 6
Variability of the curve whenY=1 in Figure 8(c).
Local Nusselt numberValue of parameterNumerical valueIncrease/decrease in %Pr=1−0.0149532239884138Pr=1.38−0.01763824746012690.706585124% decreasePr=7.2−0.02874221466481690.190789815% decreaseTable 7
Variability of the curve whenY=1 in Figure 8(g).
Average Nusselt numberValue of parameterNumerical valueIncrease/decrease in %Pr=1−0.0149532239884138Pr=1.38−0.01763824746012690.706585124% decreasePr=7.2−0.02874221466481690.190789815% decreaseThe local and average Sherwood number losses (0≤Y≤2.5) and then upsurges for mounting values of Pr, which are shown in Figures 8(e) and 8(h) (Tables 8 and 9). Heat distributes gradually, and the thermal boundary layer is lower for a more significant Prandtl number, whereas the reverse scenario is shown at lower Pr values. Figures 8(f) and 8(i) reveal that local and average shear stress decreases for upgrading Pr values. Enlarging Pr, thermal diffusivity diminishes during fluid velocity augmentations. A 3D graph comparison of the local and average Sherwood number of Pr is shown in Figure 8(j).Table 8
Variability of the curve whenY=0.5 in Figure 8(e).
Local Sherwood numberValue of parameterNumerical valueIncrease/decrease in %Pr=10.0000157634385687Pr=1.380.00006680995284110.013433293% increasePr=7.20.00009517795115240.000487423% increaseTable 9
Variability of the curve whenY=0.5 in Figure 8(h).
Average Sherwood numberValue of parameterNumerical valueIncrease/decrease in %Pr=1−0.0149532239884138Pr=1.38−0.01763824746012690.706585124% decreasePr=7.2−0.02874221466481690.190789815% decrease


## 4.3. Influence of Porosity

Figure9(a) illustrates that the flow profiles rise for growing values of the porosity parameter N. It arises due to the rise permeability of the medium, and the fluid motion is subsequently accelerated. Figure 9(b) illustrates the effect of porosity parameters N on temperature, which leads to the rise in temperature profiles (η<14), and the temperature profiles decrease due to enhancement of N. Also, it is seen that the porosity parameter strongly affects velocity profiles compared to temperature profiles. Figure 9(c) shows that the material profiles rise with rising values of N. Because the porous medium slows down the motion of the fluid, as a result, mass transfer rises. Figures 9(d) and 9(g) reveal that local and average Nusselt number declines for increasing values of N. Figures 9(e) and 9(h) expose local and average Sherwood number reductions for growing values of N. Figures 9(f) and 9(i) reveal local and average shear stress decrease for rising values of N (Tables 10 and 11). It is true that when the permeability parameter’s values increase, the tightness of the porous medium’s resistance to flow diminishes, which in turn causes the fluid’s velocity to increase. Figure 9(j) exemplifies that the flow profiles upsurge with rising values of the porosity parameter in the 3D graph.Figure 9
Influence of porosity on (a) flow, (b) thermal, (c) material fields, (d)NuL, (e) ShL, (f) τxL, (g) NuA, (h) ShA (i) τxA, and (j) flow fields in the 3D plot with M=9.5, Gr=3.5, Gm=4.5, Pr=1, Ec=0.1, Sc=0.22, R=1, and H=0.3.
(a)(b)(c)(d)(e)(f)(g)(h)(i)(j)Table 10
Variability of the curve whenY=1 in Figure 9(f).
Local skin frictionValue of parameterNumerical valueIncrease/decrease in %N=10.0043904874473810N=20.00483940448611970.044891704% decreaseN=30.00538563309237360.054622861% decreaseTable 11
Variability of the curve whenY=1 in Figure 9(i).
Average skin frictionValue of parameterNumerical valueIncrease/decrease in %N=10.0043904874473810N=20.00483940448611970.044891704% decreaseN=30.00538563309237360.054622861% decrease


## 4.4. Influence of Eckert Number

Figure10(a) shows that the flow profiles decrease (η>2) owing to rises in Ec reducing the driving force to heat transfer to the kinetic energy of the flow; that is, there is no effect of Ec on velocity profiles near the plate (Table 12). Figure 10(b) exposes thermal profile declines for rising values of Ec. The accumulation of heat energy in liquid happens due to the fact that frictional heating arises. Consequently, the fluid flow rate increases, and the heat transmission rate falls. As a result, the thermal profiles decrease. The effects of Ec on material profiles are observed in Figure 10(c). The concentration profiles increase with growing values of Ec. It allows energy to be stored in the fluid, which creates frictional heating; for this reason, the concentration profiles increase. Figures 10(d) and 10(g) reveal that local and average Nusselt numbers rise for increasing values of Ec. In Figures 10(e), 10(h), 10(f), and 10(i), the increase in the viscous dissipation term (Eckert number (Ec)) is perceived to decline the local and average Sherwood number as well as shear stress. Materially, the rise in Eckert number causes the heat energy inside the fluid layer to increase. Due to the drop in flow regimes, the heat energy formed by the Eckert number is not conserved, and there is no assurance of frictional heating.Figure 10
Influence of Eckert number on (a) flow, (b) thermal, (c) material fields, (d)NuL, (e) ShL, (f) τxL, (g) NuA, (h) ShA, and (i) τxA with M=9.5, Gr=3.5, Gm=4.5, Pr=1, Sc=0.22, N=1, R=1, and H=0.3.
(a)(b)(c)(d)(e)(f)(g)(h)(i)Table 12
Variability of the curve whenY=1 in Figure 10(a).
VelocityValue of parameterNumerical valueIncrease/decrease in %Ec=0.10.269695818Ec=0.20.19307398838.3109152% decreaseEc=0.30.1365558828.25905385% decrease


## 4.5. Influence of Schmidt Number

To be accurate, the typical values of Schmidt number (Sc) are preferred for hydrogen (Sc=0.22), ethyl ether (Sc=1.66), and toluene (Sc=1.84) at 25°C and 1 atm, which denotes chemical species dispersing in the most prevalent aspect. Figures 11(a)–11(c) exhibit characteristic behavior of flow, thermal, and material profiles for rising quantities of Sc while maintaining the other parameters constant. It can be shown in Figures 11(a) and 11(b) that the flow and thermal fields increase owing to the enhancement of Sc. Schmidt number measures the relative efficiency of momentum and mass transmission through diffusions in the hydrodynamic and material boundary layer (BL). Figure 11(c) depicts the influence of the Schmidt number on concentration profiles and see that the material profile upsurges with rising values of Sc (Table 13). It is a fact that the rising effect tends to decline the boundary layer and concentrations buoyancy effects. In this case, the flow of fluid diminishes, and the rate of concentration increases. Figures 11(d) and 11(g) reveal that local and average Nusselt numbers rise significantly for growing Sc values. Figures 11(e) and 11(h) illustrate the local and average Sherwood number rise for increasing values of Sc. Otherwise, the rate of mass transmission at the plate drops for high mass diffusivity. Figures 11(f) and 11(i) reveal local and average shear stress decreases for upgrading Sc values. The impact represents the ratio between momentum and mass diffusivity. The flow and material profiles are seen to decline as Sc rises. As a result, the material buoyancy effects lessen, which lowers the fluid velocity. Figure 11(j) shows the consequence of the Schmidt number on thermal profiles in 3D graphs.Figure 11
Influence of Schmidt number on (a) flow, (b) thermal, (c) material fields, (d)NuL, (e) ShL, (f) τxL, (g) NuA, (h) ShA (i) τxA, and (j) thermal fields in the 3D plot with M=9.5, Gr=3.5, Gm=4.5, Pr=1, Ec=0.1, N=1, R=1, and H=0.3.
(a)(b)(c)(d)(e)(f)(g)(h)(i)(j)Table 13
Variability of the curve whenY=1 in Figure 11(c).
ConcentrationValue of parameterNumerical valueIncrease/decrease in %Sc=0.220.0111813480Sc=1.660.01415115871.649894833% increaseSc=1.840.240881905015.74519072% increase


## 4.6. Comparison

The local Nusselt number acts on the enlargement of heat flow at any point, and the average Nusselt number is significant at all points. So, the average Nusselt number is equal to the local Nusselt number in these circumstances for the consequences of several factors. The local Sherwood number acts on the development of mass transfer at any point, and the average Sherwood number is substantial at all points. So, the average Sherwood number is equivalent to the local Sherwood number for the effects of numerous parameters. In these cases, local and average skin friction is shown to have identical behaviors for the impacts of rising parameters. The consequent table (Table14) presents a qualitative comparison between the described work and several researchers from the literature studies.Table 14
Comparison of the current findings qualitatively to other published results.
(a)ParametersRao et al. [32]Ali et al. [33]Present resultsUTCUTCUTCMLes——LesLesEscLesPrLesLesEscLesLesEscLesLes/EscEscNEsc——EscEscEscEscEsc(b)ParametersSambath et al. [34]Present resultsτxLShLNuLτxLShLNuLMLesEscEscLesPrLesLesLesLesLes/EscLesEcLesLesEscNEscEscEscEscLesScLesEscLesEscEscLes means lessening; Esc means escalation.This research is aimed at finding the numerical solution for an MHD heat and mass transfer flow that is viscous, electrically conducting, incompressible, radiating, and flows over a permeable plate. The BL heat absorption, radiation, chemical reaction, and Soret effect-based flow of MHD fluid through a semi-infinite vertical permeable moving plate were examined by Rao et al. [32]. Ali et al. [33] investigated MHD heat and mass transfer portents using magnetic field stimulation along a vertical plate. The effect of Sherwood number, Nusselt number, and skin friction on MHD radiative mass transfer flow from an inclined vertical plate with a heat source and sink was investigated by Sambath et al. [34].


## 5. Concluding Remarks

The numerical solutions of MHD free convection heat and mass transfer fluid flow over a vertical permeable plate with thermal radiation and heat generation are analyzed depending on the tabular and graphical implications drawn previously. The study dispensed with an FDM explaining approach with time-field contemplation, i.e., for the inconsistent results. This analysis brings out the following conclusions with engineering inquisitiveness:(i)
The velocity profile decreases for rising values ofM, Pr, and Ec and increases with rising values of N and Sc(ii)
The temperature profile upsurges owing to the various values ofM, Pr, N, and Sc and decreases for Ec(iii)
The concentration profile increases with risingPr, Ec, and Sc values and decreases for M and N(iv)
The local and average skin friction rises for rising values ofM and N. Conversely, it declines with growing values of Pr, Ec, and Sc(v)
The local and average Nusselt number increases owing to the different values ofEc and Sc. In contrast, it decreases with rising M, Pr, and N values(vi)
The local and average Sherwood number increases due to the differentM, N, and Sc values. In contrast, the local and average Sherwood number decreases for increasing values of Pr and Ec


## 6. Future Research

Numerical analysis of time-dependent oriented magneto-convective thermal-material transport across a vertical permeable plate is the subject of our research. The following recommendations can be drawn in future research with applicability from this research:(i)
Considering the influence of Hall current on the inclined plate with the flow direction of the flown flow(ii)
Considering the influence of Hall current on the flat plate(iii)
Considering the pointing vector when electric and magnetic fields are present(iv)
Considering the rotating magnetic field without the radiative parameter



---

*Source: 9977857-2023-03-02.xml*