Abstract

In this article, the effects of chemical reaction and Joule heating on MHD generalized Couette flow between two vertical porous plates with induced magnetic field and Newtonian heating/cooling have been investigated. The mathematical model used for the MHD generalized Couette flow takes into account the effect of viscous dissipation. The system of nonlinear partial differential equations governing the flow was solved numerically using the finite difference method. The resulting numerical schemes are simulated in MATLAB to obtain the profiles of the flow variables such as velocity, induced magnetic field, temperature, and concentration profiles graphically. Also, the effects of the flow parameters on the skin-friction coefficient, Nusselt number, and Sherwood number are obtained and discussed numerically through tabular forms. The findings show that an increase in the chemical reaction parameter leads to a decrease in the concentration profiles. Also, increase in the Joule heating parameter and heat generation parameter leads to an increase in the temperature profiles. Induced magnetic field profiles increase with an increase in Reynold’s number. The findings of this study are important due to its application in developing a variety of chemical technologies, including polymer manufacturing, MHD pumps, food processing, chemical catalytic reactors, astronomy, MHD flow meters, and lubrication.

1. Introduction

Magnetohydrodynamic (MHD) generalized Couette flow continues to be of interest to researchers due to its various range of practical applications in lubrication, manufacturing processes, MHD power generator, polymer processing, astrophysical fluid dynamics, plasma aerodynamics, and geophysical fluid dynamics. The study of magnetohydrodynamic generalized Couette flows of different non-Newtonian fluids has been attracted by various researchers in fluid dynamics due to their applications in engineering and industry. The laminar flow of a viscous fluid between two parallel surfaces, one of which travels tangentially with respect to the other, is known as the Couette flow. This flow can be motivated by a pressure gradient that is present in the flow direction. Several scholars have investigated the magnetohydrodynamic generalized Couette flow between two parallel vertical plates. Khan et al. [1] analyzed the second law to determine how Newtonian heating affects the Couette flow of a viscoelastic dusty fluid and the transfer of heat in a rotating frame. Umavathi et al. [2] studied the Poiseuille–Couette flow with heat transfer in the inclined channel. They concluded that increases in the Grashof number, angle of inclination, and height ratio increase velocity, while increases in the Hartmann number, viscosity, and conductivity ratios decrease the velocity profiles. Bég et al. [3] analyzed the oblique magnetic field in a rotating highly permeable medium. Barikbin et al. [4] analyzed non-Newtonian fluid for MHD Couette flow using the Ritz–Galerkin method. Sreekala and Reddy [5] studied the effect of inclined magnetic field with the steady MHD Couette flow. Joseph et al. [6] investigated heat transfer in inclined magnetic field with unsteady MHD Couette flow between two infinitely parallel porous plates. The temperature-dependent transient MHD Couette flow and heat transfer of dusty fluid were studied by Mosayebidorcheh et al. [7]. Jha et al. [8] investigated thermal radiation with unsteady MHD free convective Couette flow. Ngiangia and Okechukwu [9] discussed the impact of variable electroconductivity and radiation. They indicated that increases in electroconductivity, Prandtl number, Reynolds number, and Grashof number result in an increase in velocity distribution, while increases in the magnetic field result in a drop in velocity. Raju et al. [10] discussed the influence of diffusion thermo and thermal diffusion on a natural convection Couette flow using the finite element method. Ali et al. [11] investigated convective cooling of the nanofluids in a rotating system. Job and Gunakala [12] studied the thermal radiation for vertical permeable plates using Galerkin’s finite element method. Also, Ali et al. [13] analyzed the Couette flow of a maxwell fluid for three dimensional with periodic injection/suction. Hussain et al. [14] analyzed instability of the MHD Couette flow of an electrically conducting fluid. Anyanwu et al. [15] discussed the influence of radioactive and a constant pressure gradient on an unsteady MHD Couette flow.

Chemical reaction can be categorized as homogeneous or heterogeneous reaction which can happen in a single phase or at the surface. The effects of chemical reaction depend on the order of the reaction rate for both heterogeneous and homogeneous reaction. The study of chemical reaction processes is beneficial in developing a variety of chemical technologies, including polymer manufacturing, MHD pumps, MHD flow meter, chemical catalytic reactors, and food processing. The chemical reactions have numerous applications in many sectors of science, engineering, and technology. The two-step exothermic chemical reaction of generalized Couette hydromagnetic flow with thermal criticality and entropy generation was presented by Kareem and Gbadeyan [16]. Jha and Sarki [17] analyzed the effects of Dufour and chemical reaction on a moving vertical porous plate. Soumya et al. [18] studied the effects of chemical reaction, variable thermal conductivity, and nanoparticles injection on the MHD three-dimensional free convective Couette flow. Reddy et al. [19] discussed the effects of chemical reaction and thermal radiation on an unsteady MHD Couette flow through a porous medium with periodic wall temperature. The study reveals that increasing effect of chemical reaction parameter leads to decrease in the concentration profiles. Increasing value of magnetic parameter leads to decrease in the velocity profiles. Ajibade and Umar [20] analyzed the effect of chemical reaction and radiation absorption on the unsteady MHD free convective Couette flow in a vertical channel. Apelblat [21] discussed analytically the effects of first order chemical reaction with mass transfer. Prakash et al. [22] analyzed the effects of chemical reaction and thermo-diffusion on the MHD free convective Couette flow. Abbaszadeh et al. [23] investigated the effects of two step exothermic reaction on generalized Couette hydromagnetic flow using the direct meshless local Petrov–Galerkin method. Chen and Arce [24] discussed the mathematical and numerical approach for convective-diffusive mass transfer with chemical reaction in the Couette flow. Kareem and Gbadeyan [16] studied entropy generation and thermal criticality of generalized Couette hydromagnetic flow of two-step exothermic chemical reaction in a channel. The study concluded that low-dissipation rates can promote irreversibility processes and low-heat source terms can reduce thermal explosion.

An electric current flows through a resistance, it can be said to “Joule heating” the object being passed through. Joule heating is important due to its various applications including portable fan heater, bulb that emits light, clothes iron, and the hot flow produced by a hairdryer. Sarkar [25] investigated the MHD Couette–Poiseuille flow with laminar forced convection and viscous and Joule dissipations. Zueco et al. [26] studied the effects of Joule heating, Hall and ion slip in a non-Darcian porous media channel with a nonlinear transient hydromagnetic partially ionized the dissipative Couette flow. Ramesh [27] analyzed the effects of Joule heating and viscous dissipation on a Jeffrey fluid’s Couette and Poiseuille flows by taking into account the slip boundary conditions. Bég et al. [28] analyzed numerically the effects of Hall current, Joule heating, viscous flow, and Ion slip on magnetohydrodynamic Hartmann–Couette flow and heat transfer in a Darcian channel. Saleel et al. [29] discussed the microfluidics flow control using the electro-osmotic Couette flow and Joule heating. Shamshuddin et al. [30] investigated the effects of Joule heating and viscous dissipation in non-Fourier MHD squeezing flow, as well as heat and mass transfer between Riga plates with thermal radiation. Shamshuddin and Satya Narayana [31] studied the effect of viscous dissipation and Joule heating combined with a Cattaneo–Christov heat flux on MHD flow past a Riga plate.

Induced magnetic on the hydromagnetic generalized Couette flow is beneficial in MHD pumps and MHD flow meters. The effects of induced magnetic field have numerous applications in many sectors of science, engineering, and technology. The magnetic field is modified by the induced magnetic field, which also creates its own magnetic field in the fluid. The flow of the fluid through the magnetic field also generates mechanical forces that change the fluid’s motion. Consequently, it is necessary to include the influence of the induced magnetic field in the hydromagnetic equations in a number of physical conditions. In the metallurgical sector, magnetic fields are employed for stirring, heating liquid metals, and levitating. Sarveshanand and Singh [32] discussed the effects of induced magnetic field on magnetohydrodynamic free convection between vertical parallel porous plates. The study concluded that increasing suction parameter leads to decrease the induced current density and velocity field and while it has increasing effect on the induced magnetic field. Mng’ang’a et al. [33] investigated the effects of induced magnetic field in the direction of the fluid flow on hydromagnetic surface driven flow. Kwanza and Balakiyema [34] analyzed the magnetic induction on MHD free convective flow past an infinite vertical porous plate. The study reveals that the increase of the Prandtl number and the suction velocity lead to a decrease in the temperature, induced magnetic field and velocity profiles in the boundary layer region whereas increase in the magnetic field parameter and Grashof number leads to increased induced magnetic field profiles. Singh et al. [35] discuss the effects of induced magnetic field on hydromagnetic free convection. Akbar et al. [36] investigated the effects of induced magnetic field on interaction of nanoparticles for the peristaltic flow in an asymmetric channel. Sher Akbar et al. [37] discuss the effects of induced magnetic field and heat flux in a permeable channel.

The novelty of this work is to investigate the combined effects of chemical reaction and Joule heating on the MHD generalized Couette flow between two parallel vertical porous plates with induced magnetic field and Newtonian heating/cooling by taking into the account the effects of viscous dissipation. The system of nonlinear partial differential equation of the governing equations was obtained using the finite difference method and implemented in MATLAB. The effects of various dimensionless parameters on velocity profiles, temperature profiles, concentration profiles, and induced magnetic field profiles are discussed in detail. The results of this problem can be helpful in various science and engineering. The results of this problem can be helpful in various devices subject to significant variations in the gravitational force, its application on heat exchanger designs, wire and glass fiber drawings, and its application in nuclear engineering in connection with reactor cooling. In line with the aforementioned objectives and aim, this study proposes answers to the following research questions.(i)How does the flow variables such as temperature, velocity, concentration, and induced magnetic field profiles varies with Prandtl number, Schmidt number, Grashof numbers for heat and mass transfer, Eckert number, Joule heating parameter, heat generation parameter, Reynold’s number, magnetic Prandtl number, suction/injection parameter, permeability parameter, chemical reaction parameter, and magnetic parameter?(ii)What are the variations in the Sherwood number, Nusselt number, and skin friction coefficient due to the enrichment in Joule heating parameter, heat generation parameter, Reynold’s number, Prandtl number, Schmidt number, Grashof numbers for heat and mass transfer, Eckert number, magnetic Prandtl number, suction/injection parameter, magnetic parameter, chemical reaction parameter, and permeability parameter?

2. Mathematical Formulation

Let us consider the bidimensional unsteady laminar flow of incompressible, electrically conducting, and viscous between two infinite parallel vertical porous plates with magnetic field strength vector , constant suction velocity , and a constant magnetic field applied perpendicular to the plates as depicted in Figure 1. Fluid flow between two parallel vertical porous plates with infinite lengths at and .

The fluid is flowing through a porous medium that is assumed to obey Darcy’s law. Electrically conducting fluid produce induced currents then produce an induced magnetic field in and direction, respectively, that opposite the applied field, resisting the conductor’s motion.

At initial concentration and temperature at time , the fluid is initially stationary, as are the plates at and . When , the porous plate at begins to move impulsively in its own plane at a constant velocity of , and its concentration and temperature rise to and , respectively. In contrast, the other porous plate, which is located at a distance of from it, remains fixed and maintains its concentration and temperature at and . The concentration, induced magnetic profile, temperature, and velocity are functions of and since the porous plates have an unlimited length. In the direction, the fluid is suctioned with constant velocity .

The unsteady MHD Couette flow between two vertical porous plates is governed by the following nonlinear partial differential equations. Using Boussinesq’s approximation, the governing equations for the current physical system in dimensional form are as follow [8]:-

Continuity equation

Momentum equation

Energy Equation

Concentration equation

Equation of magnetic induction

The corresponding initial and boundary conditions are as follows[8]:-

To nondimensionalize the governing equations, we use the following dimensionless quantities for the present hydromagnetic problem

Using the dimensionless variables defined in equations (2)–(6) and (8) in nondimensional forms are as follows:-where represents Joule heating parameter, represents magnetic parameter, represents injection/suction parameter, represents the permeability parameter, represents the Reynolds number, represents Eckert number, represents chemical reaction parameter, represents Schmidt number, represents Prandtl number, represents internal heat generation parameter, represents magnetic Prandtl number, represents Grashof number for heat transfer, and represents Grashof number for mass transfer.

The corresponding initial and boundary conditions are as follows:

3. Numerical Solution

The governing (9)–(13) are nonlinear partial differential equations, therefore cannot be solved analytically. The numerical method for nonlinear partial differential equations of momentum, concentration, energy, and induced given in (9)–(13) are solved using the finite difference method subject to the initial and boundary condition (14). The condition for time stability is the Courant–Friedrichs–Lewy or CFL condition which depends on space and time discretization. A mesh is fixed at and to ensure the stability and convergence. The transport equations (9)–(13) at the grid point (i, j) are expressed in difference form using Taylor’s series expansion. Thus, the values of V, C, Hx, Hy, and T at grid point t = 0 are known; hence the velocity, induced magnetic field, concentration, and temperature fields have been solved at time using the known values of the previous time for all . These processes are repeated till the required solution of velocity, induced magnetic field, concentration, and temperature field is converge.

The forward time central in space is used as the numerical scheme, in finite difference form can be approximated as

Using equation (15) into equations (9)–(12), in finite difference form the governing equations becomes

Subject to the following initial and boundary conditions

For practical engineering application and design, the quantities of interest are Sherwood number (Sh), skin friction coefficient , and Nusselt number (Nu) which are defined aswhere , , and

Substituting equation (8) into equation (18), we obtain

Both temperature, velocity, concentration, and induced magnetic field profiles were obtained and utilized to compute numerical values of the Sherwood number, skin friction coefficient, and Nusselt number in equation (19).

4. Results and Discussion

To investigate the physical significance of the modeled MHD generalized Couette flow, the influence of various physical parameters involved in the problem such as Joule heating parameter, heat generation parameter, Reynold’s number, Prandtl number, Schmidt number, Grashof numbers for heat and mass transfer, Eckert number, magnetic Prandtl number, suction/injection parameter, magnetic parameter, chemical reaction parameter, and permeability parameter on the nondimensional temperature, velocity, concentration, and induced magnetic field profiles are evaluated numerically and executed graphically and discussed to get the physical insight into the problem, whilst the numerical values for skin friction coefficient, Nusselt, and Sherwood numbers are represented in the tabular form. The numerical solutions have been conducted by adopting the default values of various involved physical parameters in the problem such as , and until otherwise specified particularly.

Figure 2 depicts the effects of Grashof number for mass transfer on velocity profiles. Graphically, it is observed that the velocity of the fluid increases with an increase in the Grashof number for mass transfer , since Grashof number for mass transfer represents the ratio of the species buoyancy force to the viscous force. Physically, increasing Grashof number leads to a decrease in the viscosity of the fluid which results to a decrease in the viscous force which leads to an increase in the species Bouyant force and thus increase the velocity of the fluid.

Figure 3 depicts the effects of Grashof number for heat transfer on velocity profiles. Graphically, it is observed that increasing the values of Grashof number for heat transfer leads to an increase in the velocity profiles. By definition, Grashof number for heat transfer represents the ratio of the thermal buoyancy force to the viscous force. Physically, increasing the Grashof number reduces the viscous force and increase thermal buoyancy force thus leads to an increase in the velocity profiles.

The effect of magnetic parameter on velocity profiles is shown in Figure 4. Graphically, it is observed that increasing the values of the magnetic parameter leads to a decrease in the velocity profiles, since magnetic parameter represents the ratio of electromagnetic force to an inertia force. Physically, increasing magnetic parameter leads to an increase in the electromagnetic force which caused by the magnetic field that control the fluid flow characteristics. The presence of a magnetic field in an electrically conducting fluid enhance the Lorentz force which acts against to the fluid flow when the magnetic field is applied normal to the direction of the fluid flow which leads to a slowly the motion fluid and thus decrease the velocity of the fluid.

Figure 5 illustrates the effects of permeability parameter on velocity profiles. Graphically, it is observed that increasing the values of the permeability parameter leads to a decrease in the velocity profiles. Physically, increasing the permeability parameter which opposes the flow and leads to enhance deceleration of the flow. Therefore, with an increase in the permeability parameter causes the resistance to the fluid motion which resulted to decrease the velocity boundary layer and thus the velocity profiles decrease.

The effects of suction/injection parameter on velocity profiles are shown in Figure 6. Graphically, it is observed that increasing the values of the suction parameter leads to a decrease in the velocity of the fluid. The injection parameter increases with a decrease in the velocity of the fluid. Physically, increasing suction parameter implies the removal of the fluid from the flow system which causes destabilization the velocity boundary layer leads to a decrease in the motion of the fluid and thus decrease the velocity profiles.

Figure 7 exhibits the effects of the Reynolds number on velocity profiles. Graphically, it is observed that increasing the values of the Reynolds number leads to an increase in the velocity profiles, since the Reynolds number represents the ratio of inertia force to the viscous force. Physically, increasing the Reynolds number reduces the viscous forces in the fluid and thus result to an increase in the velocity profiles.

Figure 8 presents the effects of chemical reaction parameter on concentration profiles. Graphically, it is observed that increasing the values of the chemical reaction parameter leads to a decrease in the fluid concentration. The positive values of indicate the destructive type of the chemical reaction. The influence of a destructive chemical reaction parameter caused a decrease in the concentration diffusion species. Physically, increasing chemical reaction for destructive results to molecular motion become higher, which enhances the transport phenomenon, thus decreases the concentration profiles.

Figure 9 illustrates the effects of suction/injection on concentration profiles. Graphically, it is observed that increasing the values of the suction leads to a decrease in the concentration profiles. Physically, increasing suction/injection parameter results in the thinning of the concentration boundary layer of the fluid flow which leads to a decrease in the concentration profiles.

The effects of Schmidt number on concentration Profiles are shown in Figure 10. Graphically, it is observed that the concentration profiles decrease as the Schmidt number increases. By the definition, the Schmidt number represents the ratio of kinematic air viscosity to the mass diffusivity. Physically, increasing Schmidt number reduces the mass diffusivity, which results to a decrease in the concentration profiles.

The effects of magnetic Prandtl number on induced magnetic field profiles are shown in Figures 11 and 12. Graphically, it is observed that increasing the values of the magnetic Prandtl number leads to a decrease in the induced magnetic field profiles. Physically, increasing the magnetic Prandtl number reduces the magnetic diffusivity, which leads to a decrease in the induced magnetic field by the motion of a conducting medium and thus the induced magnetic field profiles decrease.

The effects of suction/injection parameter on induced magnetic profiles are shown in Figures 13 and 14. Graphically, it is observed that increasing the values of the suction parameter leads to a decrease in the induced magnetic field profiles. Increasing suction/injection parameter leads to a decrease in the fluid the flow system and thus decreases the velocity of the fluid which reduces the interaction between the fluid and magnetic field and thus decreases the induced magnetic field profiles.

The effects of the Reynolds number on an induced magnetic field profiles along direction are shown in Figure 15. Graphically, it is observed that increasing the values of the Reynolds number leads to an increase in the induced magnetic field. By the definition, the Reynolds number represents the ratio of inertia force to the viscous forces. Physically, increasing the Reynolds number reduces the viscous force which leads to an increase of interaction between the fluid and magnetic field and thus increases induced magnetic field profiles.

The effects of the Reynolds number on induced magnetic field profiles along direction are shown in Figure 16. Graphically, it is observed that increase in the Reynolds number leads to a decrease in the induced magnetic field. The Reynolds number represents the ratio of inertia force to the viscous forces. Physically, increasing the Reynolds number leads to a decrease of interaction between the fluid and magnetic field and thus decreases induced magnetic field profiles.

Figure 17 demonstrates the effects of suction/injection parameter on temperature profiles. Graphically, it is observed that increase in the suction/injection parameter leads to an decrease in the fluid temperature. Physically, increasing suction/injection parameter results in the thinning of the thermal boundary layer of the fluid flow which leads to decrease the temperature profiles.

The effects of the Eckert number on temperature profiles are shown in Figure 18. Graphically, it is observed increasing the values of the Eckert number leads to an increase in the temperature profile. By the definition, the Eckert number represents the ratio of the fluid flow kinetic energy to the enthalpy of the fluid. Physically, increasing the Eckert number reduces the enthalpy of the fluid and thus increases the temperature of the fluid. The motion of the fluid increases as an increase of the Eckert number which transforms into kinetic energy which lead to an increase in the temperature profiles.

The effects of the Prandtl number on temperature profiles are shown in Figure 19. Graphically, it is observed that the temperature of the fluid decreases as the values of the Prandtl number increases. By the definition, the Prandtl number represents the ratio of viscous force/momentum diffusivity to the thermal diffusivity. Physically, increasing the Prandtl number reduces the thermal diffusivity and increase viscosity of the fluid which leads to an increase of temperature of the fluid.

Figure 20 illustrates the effects of varying heat generation parameter on temperature profiles. Graphically, it is observed that temperature profile is greatly influenced by heat generation parameter . Physically, increasing heat generation parameter leads to an increase in the temperature of the fluid. By increasing heat generation parameter leads to an increase in kinetic energy which leads to an increase in the thermal boundary layer and thus results to an increase in the temperature profiles.

The impacts of varying Joule heating parameter on temperature profiles are shown Figure 21. Graphically, it is noticed that the temperature profiles increase with an increase in the Joule heating parameter . Physically, increasing the Joule heating parameter leads to an increase in the friction, and thus the mechanical energy converted into thermal energy which results to an increase in the temperature of the fluid and thus temperature profiles increase.

Table 1 shows the effects of the Prandtl number, Eckert number, heat generation parameter, and Joule heating parameter on the rate of heat transfer (Nusselt number) at the stationary and moving plates. It is observed from the table that the Nusselt number at the stationary plate decreases with an increase in the Prandtl number whereas increase with an increase in Joule heating parameter, heat generation parameter, and Eckert number. The Nusselt number at the moving plate increases with an increase in Joule heating parameter, Eckert number, and heat generation parameter whereas decrease with an increase in the Prandtl number.

Table 2 shows the effects of magnetic parameter, Prandtl number, Reynolds number, heat generation parameter, and Joule heating parameter on the wall shear stress (skin friction coefficient) at the moving and stationary plates. It is observed from the table that the wall shear stress at the stationary plate increases with an increase in the Prandtl number, magnetic parameter, and Reynolds number whereas decrease with an increase in heat generation parameter and Joule heating parameter. The skin friction coefficient at the moving plate increases with an increase in the Prandtl number, magnetic parameter, and Reynolds number whereas decrease with an increase in heat generation parameter and Joule heating parameter.

Table 3 shows the effects of Schmidt number and chemical reaction parameter on the Sherwood number at the moving and stationary plates. It is observed from the table that the rates of mass transfer (Sherwood number) at the moving and stationary plates decrease with an increase in the chemical reaction parameter and Schmidt number.

5. Concluding Remarks

We examine the effects of chemical reaction and Joule heating on the MHD generalized Couette flow between two vertical porous plates with an induced magnetic field and Newtonian heating/cooling. The system of the governing equations of partial differential equations of the model is solved numerically using the finite difference method. The numerical results for the temperature, velocity, concentration, and induced magnetic field profiles are presented graphically for the pertinent parameters. The numerical values for skin friction coefficient, Nusselt, and Sherwood numbers are represented in the tabular form. The major findings of the present analysis can be summarized in the following points:(i)The velocity profiles increase with rising the values of , , and while increasing and have an opposite effect on the velocity profiles(ii)The temperature of the fluid increase with an increasing the values of , , and while increasing and leads to a decrease in the temperature fluid(iii)The rates of heat transfer at the moving and stationary plates increases with an increase in the , , and while decrease with an increase in (iv)The skin friction coefficient (wall shear stress) at the stationary and moving plates increases with an increase in , , and while decrease with an increase in the and (v)The rate of mass transfer at the moving plate decreases with an increase in the and . However, at the stationary plate decreases with an increase in and (vi)An increasing , , and decrease the concentration profiles(vii)An induced magnetic field profiles decrease with rising the values of , , and

Future research can be conducted on the MHD generalized Couette flow between two parallel vertical porous plates in presence of an inclined variable magnetic field.

Nomenclature

:Acceleration due to gravity
:Uniform applied magnetic field along x-axis
:Nondimensional velocity component along y-direction
:Coefficient of thermal conductivity
:Time
:Dimensionless time
:Induced magnetic field along y-direction
:Induced magnetic field along x-direction
:Current density
:Magnetic field strength
:Velocity components along y axis
:Concentration susceptibility (m/mole)
:Mean fluid temperature (k)
:Temperature
:Species concentration
:velocity of the fluid
:Specific heat at constant pressure
:Temperature of the fluid in the moving plate
:Temperature of the fluid in the stationary plate
:Concentration of the fluid in the moving plate
:Concentration of the fluid in the stationary plate
:Mass diffusivity/chemical molecular diffusivity
:Thermal diffusion ratio
:Injection/suction velocity
:Distance between two parallel plates
:Nondimensional induced magnetic field
:Constant velocity of the moving plate
:Viscosity
:Fluid density
:Gradient operator
:Coefficient of volume expansion for heat transfer
:Coefficient of thermal expansion for concentration gradient
:Electrical conductivity
:Time and distance intervals
:Nondimensional temperature
:Nondimensional concentration
:Kinematic viscosity
:Magnetic permeability
:Internal heating generation parameter/Newtonian heating/cooling parameter
:Chemical reaction parameter
:Reaction coefficient
:Joule heating parameter.

Data Availability

The data used to support the findings of this study are included within the article.

Conflicts of Interest

The author declare no conflicts of interest.

Acknowledgments

The author wish to express their very sincere thanks to anonymous reviewers for their valuable comments and suggestions.