MHD flow and effects of hall current joule heating and heat source across a chemically reacting and
* Author to whom correspondence should be addressed.
Sigma Journal of Engineering and Natural Sciences 2025, Vol. 43, Issue 4, pp. 1179-1196; doi.org/10.14744/sigma.2024.00110
Abstract
Keywords: bvp4c; Chemical Reaction; Heat Transfer; Mass Transfer; Porous Medium; Radiation
Introduction
Magnetohydrodynamics (MHD) is an interdisciplinary discipline that investigates the behavior of electrically conducting fluids such as plasmas and liquid metals in the presence of magnetic fields. It investigates the complex interactions between fluid motion and magnetic forces,
including phenomena such as electric current production and magnetic field induction. MHD possesses applications in a variety of domains, including fusion research, astronomy, geophysics, material processing, and aerospace engineering. MHD, which combines fluid dynamics and electromagnetic principles, provides theoretical insights
*Corresponding author. *E-mail address: mth2291004_ankur@cottonuniversity.ac.in This paper was recommended for publication in revised form by Editor-in-Chief Ahmet Selim Dalkilic Published by Yıldız Technical University Press, İstanbul, Turkey © Author. This is an open access article under the CC BY-NC license (http://creativecommons.org/licenses/by-nc/4.0/).
Sigma J Eng Nat Sci, Vol. 43, No. 4, pp. 1179−1196, August, 2025
into complicated fluid behavior as well as practical solutions for technological breakthroughs that have the potential to revolutionize energy generation, space exploration, and industrial operations. Hall current refers to the flow of electric charge caused by the Hall effect, which occurs perpendicular to both an electric and magnetic field. When a conducting material carrying an electric current is exposed to a perpendicular magnetic field, the Lorentz force operates on the moving charges, deflecting them to one side of the conductor. This deflection causes charge to accumulate on one side, generating an electric field that resists the passage of charge. As a result, a secondary electric current, known as the Hall current, runs perpendicular to both the initial current and magnetic field. The Hall current is essential for a variety of applications, including magnetic field measurement [1], semiconductor device characterization [2], and the functioning of certain types of sensors [3] and electronic equipment [4]. In magnetohydrodynamics (MHD), joule heating is the process of creating heat within a conducting fluid (such as plasma or liquid metal) owing to the dissipation of electrical currents in the presence of a magnetic field. When an electrically conducting fluid flows through a magnetic field, it induces electric currents inside itself. These currents encounter resistance as they pass through the fluid, resulting in the conversion of electrical energy into heat. This behavior is most noticeable in plasmas, where the high temperatures and magnetic fields found in fusion experiments or astrophysical settings may produce considerable Joule heating effects. Joule heating in MHD systems can have a significant impact on energy transmission, plasma stability, and overall system performance. Understanding and managing joule heating is critical for improving MHD-based technologies, including fusion reactors, magnetohydrodynamic generators, and space propulsion systems. Efforts to regulate joule heating in MHD systems frequently include measures for reducing electrical resistance, controlling current densities, and mitigating thermal impacts to improve system performance and dependability. Ali et al. [5] explored the effects of Hall current on a vertically stretching sheet. The impact of radiation on MHD Newtonian fluid over an exponentially stretched sheet was studied by Kameswaran et al. [6]. The flow and heat transfer along a surface that is stretched exponentially and buried in porous material with a fluctuating surface heat flux were studied by Mandal and Mukhopadhyay [7]. The sequel of varying thermal conductivity and non-uniform heat source/sink on MHD flow and heat transfer across an exponentially stretched sheet in a porous medium saturated with Maxwell fluid were studied by Singh and Agarwal [8]. In a porous media with fluctuating viscosity, Khidir and Sibanda [9] probed MHD mixed convective flow, heat and mass transfer over an exponentially stretched surface, and the Soret and Dufour effects. Jaber [10] investigated on the
sequel of Hall current and changing viscosity in presence of radiation via stretched sheet. Gorla et al. [11], Hayat et al. [12], Eid [13] have investigated on heat source/sink, velocity slip condition and heat generation/absorption impacts on exponentially stretching sheet respectively. Srinivasacharya and Jagadeeshwar [14] studied the consequences of Joule heating and Hall current on an exponentially stretched sheet by usage of Chebyshev pseudo spectral method. Govindarajan and Vijayalakshmi [15] examined the impact of Hall current in a dusty fluid across a saturated permeable sheet. Veera Krishna et al. [16] investigated the sequel of radiation and Hall current when a vertical channel filled with porous media experiences an irregular flow. While accounting for the Hall current, Veera Krishna and Chamkha [17] investigated the sequel of heat generation/absorption and thermo-diffusion on an unstable second-grade fluid flow near an infinite vertical plate via porous media. Sarma and Pandit [18] studied the impacts of Hall current, Soret and rotation collectively via infinitely accelerated vertical plate. Rajput and Kanaujia [19] probed the behaviour of Hall current over a rotating inclined plate. Again, Rajput and Kumar [20] researched on the radiation and Hall current behaviour through a temperature-varying inclined plate. Reddy and Sreedevi [21] delved into the study of nanofluid flow via an inclined plate in occurrence of porous medium and chemical reaction. Goud et al. [22] probed the outcomes of Joule heating over a vertical plate with an inclined permeable surface when a chemical reaction is occurring. Sandhya et al. [23] probed the heat and mass transfer flow across an inclined plate in presence of radiation and chemical reaction. Vijayaragavan et al. [24] worked on the comprehensive study of MHD Casson fluid flow across an tilted plate. Dwivedi and Singh [25-26] indulged themselves into the study of cylindrical flow considering the sequel of Hall current. Jeevitha et al. [27] studied flow behaviour via a rotating vertical cone in occurrence of porous medium. In occurrence of Joule heating, Tasnim et al. [28] examined the sequel of MHD and entropy creation on a square enclosure filled with nanofluid and equipped with several heat-generating devices. Recently, several researchers [29-37] have considered the effects of hall current on the various fluid flow geometry. The transmission of heat and mass along an exponentially expanding surface has been the subject of previous research, but the combined sequel of slip and Hall currents in the occurrence of radiation and heat creation via a porous medium has received less attention. In the occurrence of radiation and heat production, together with Hall currents, this article concentrates on a novel study of mixed convection slip flow, heat, and mass transfer over a porous exponentially stretched sheet under suction or injection. We have considered a viscous Newtonian fluid in presence of porous medium using the combination of Brinkman equation with the Navier- Stokes equation
Sigma J Eng Nat Sci, Vol. 43, No. 4, pp. 1179−1196, August, 2025
[38,39]. This model of porous medium is known as the Darcy-Brinkman model of porous medium. This investigation has several applications in the fields of engineering, heat exchangers, mixing tanks, fluidized bed reactors [40], polycrystalline semiconductors [41], sensors [3], semiconductor detectors [2], and different electronic equipment [4].
Mathematical Formulation
Consider an incompressible viscous fluid with temperature Tw(x) and concentration Cw(x) flowing continuously over an exponentially stretched sheet via electrical conductivity. The z-axis corresponds with the leading edge of the sheet, the y-axis is perpendicular to the sheet in the outward direction of flow, and the positive x-axis is taken along the sheet in the direction of flow. The distance x from the slit causes an exponential shift in the sheet velocity. It is assumed that temperature and concentration of the ambient medium are T∞ and C∞. An intense magnetic field with strength B(x) is applied in the y direction, taking into account the effects of Hall current. The induced magnetic field is disregarded when a very low magnetic Reynolds number is assumed. This presumption is correct since liquid metals and partially ionized fluids which are frequently used in industrial applications have very low magnetic Reynolds numbers. The cross flow in the z-direction brought on by Hall current›s occurrence makes the flow three dimensional [42]. The concentration and temperature of the plate are denoted by C and T, respectively, whereas the velocity vector is (u, v, w). Figure 1 exhibits the geometrical representation of the flow. The following equations govern the current flow problem in the case of joule heating under the previously mentioned assumptions, together with Boussinesq and boundary layer approximations [14]:
(5) The relevant boundary conditions [14] associated with the surface are: (6) (7) (8) is the magnetic field and B0 is the Here is the stretching constant magnetic field, velocity, where U0 is the reference velocity, is the special velocity at the wall, where V0 is the initial is the velocity slip factor, strength of suction, where in case of no slip, N = 0 qr is described as follows
where σ1 indicates Stefan-Boltzmann constant and k represents Rosseland mean absorption coefficient. We presume that the distinction in internal flow temperature is suitably minimal, as shown by the Taylor series expansion about, while ignoring the terms with higher order as,
Now we introduce the similarity transformation [14] and the nondimensional quantities as follows:
Sigma J Eng Nat Sci, Vol. 43, No. 4, pp. 1179−1196, August, 2025
Using the above similarity transformations and nondimensional quantities the PDEqs. (2)-(5) are transformed into the following nonlinear ODEs as follows: (9)
(12) The boundary conditions (6)-(8) are transformed to the following form: (13) (14) where S < 0 respectively,
Materials And Methods
Boundary value problems, either linear or nonlinear, are evaluated using MATLAB solution bvp4c. It can accurately calculate y(x) for each x in [a,b] while taking into consideration the stage-by-stage boundary limitations. This method changes the boundary conditions at infinity to those at a point when it is reasonable to attempt to solve the current issue. The three-stage Lobatto IIIa formula is executed as a finite difference code by bvp4c. Shampine et al. [43] suggested this Lobato IIIa approach. This collocation polynomial divides the interval into subintervals based on a mesh of points, providing uniform fourth-order correct solutions across the interval. The solver makes an estimate of the numerical error on each subinterval. If the answer doesn›t satisfy the tolerance requirements, the mesh is adjusted, and the procedure is repeated. However, in addition to the initial mesh points, an initial assumption about the solution at the grid points is unavoidably needed. By adjusting the step size as it goes along and offering a preliminary estimate at the beginning mesh points, the algorithm produces specified precision. Equations (9), (10) and (11) are transformed as follows in order to apply the finite difference-based solver bvp4c:
is the suction/injection as S > 0 or represents the velocity slip
parameter. The expression of the various engineering terms, such as the wall shear stress in x and z- direction, heat and mass transfer rates, respectively are given as follows:
(15) (20) The local skin friction coefficient along x-direction Cx, local skin friction coefficient along z-direction Cz, local Nusselt number Nu and local Sherwood number Sh are defined as follows: (16) Using the nondimensional quantities the following results are obtained:
(21) The initial and boundary conditions (10) are transformed as follows: (22) (23)
To satisfy the convergence requirements, the accuracy is measured up to the sixth decimal place. Shampine et al. [43] gave the detailed explanation of this scheme. The algorithm of this scheme in given by Figure 2 [44].
Sigma J Eng Nat Sci, Vol. 43, No. 4, pp. 1179−1196, August, 2025
S, λ, ζ, Bu, δ, R, H, h, Sp, Sc, Ec and K for varying values of Pr.
Results And Discussion
Figure 2. Algorithm of bvp4c. [From Wahid et al. [44], with permission from Springer.]
Validation Table 1 provides a verified example of the correctness and precision of the methodology used in this work. In this instance, the bvp4c approach is applied to compare the outcomes generated by [45] and [14] for f ''(0) and f(∞), disregarding the values of S, λ, ζ and H are shown in Table 1. The approach is dependable since the current study›s calculated results are consistent with the earlier findings. Additionally, the values of −θ'(0) for various values of Pr show good agreement with the present results when comparing the equations discussed by Ishak [46] with the current equations. The findings for −θ'(0) are presumably compared using the bvp4c approach, ignoring the values of
Under the aforementioned considerations, the problem is numerically solved, and graphs from Figures. 3 to 33 illustrates the behaviors of tangential flow, cross flow, temperature and concentration for parameters, such as the porosity parameter (Sp), mixed convection parameter (ζ), radiation parameter (R), Magnetic parameter (H), Eckert number (Ec), Hall parameter (h), Bouyancy ratio parameter (Bu), Prandtl number (Pr), heat source (δ), Schmidt number (Sc), and chemical reaction parameter (K). The following standard values for the parameters are taken from [14] as Pr = 0.71, Ec = 0.5, H = 1.0, h = 2.0, Bu = 0.5, ζ = 1, R = 1, Sp = 0.4, Sc = 0.22, K = 1, δ = 1, S = 0.5, and λ = 1.0. Furthermore, variations in skin friction (τ), nusselt number (Nu) and sherwood number (Sh), under various conditions in relation to the different flow parameters are explored and illustrated in Tables 3 and 4. Tangential Velocity Profile (f'(η)) The impact of Bu on f '(η) is exhibited in Figure 3. It can be perceived that as Bu rises, f '(η) also rises near the sheet, whereas it shows reverse trend as we move away from sheet. The impact of H on f '(η) is exhibited in Figure 4. It exhibits that when the value of H escalates, the tangential velocity decreases near the sheet. A homogeneous magnetic field normal to the direction of flow is applied to produce the Lorentz force. The fluid›s velocity in the boundary layer tends to be slowed down by this force. As a result, the tangential velocity decreases as H increases. Also, it is depicted that the tangential velocity escalates with the rise in H away from the sheet. Figure 5 shows how the Hall parameter (h) effects f '(η). It is perceived from Figure 5 that as h escalates, so does the tangential velocity near the sheet. Effective conductivity decreases with increasing h, and this lowers the magnetic
Table 1. Comparison of f ''(0) and f(∞) with earlier results Magyari and Keller [45]
Sigma J Eng Nat Sci, Vol. 43, No. 4, pp. 1179−1196, August, 2025
damping force acting on f '(η). As a result f '(η) rises near the sheet, whereas it decreases as we move away from the sheet. Figure 6 shows the effect on f '(η) for higher values of Pr. f '(η) decreases as a result of the escalating values of Pr lowering the momentum boundary layer size. Prandtl number indicates how the thermal conductivity is related to the fluid viscosity at reference and the specific heat at constant pressure. Since the viscous forces are too strong to
keep the fluid particles securely at constant pressure, f '(η) decreases as Pr increases. The effects of Sp on f'(η) are seen in Figure 7. It is discovered that increasing Sp values decreases the velocity profiles and increases the thickness of the velocity boundary layer near the sheet, whereas it shows reverse trend as we move away from the sheet. As the flow resistance escalates in occurrence of porous medium, f'(η) reduces. It is evident
Sigma J Eng Nat Sci, Vol. 43, No. 4, pp. 1179−1196, August, 2025
Figure 7. f '(η) vs Sp. from Figure 8 that as the values of ζ grows, so does the tangential velocity near the sheet, whereas it shows reverse trend as we move away from the sheet. This is a result of positive values of ζ create a pressure gradient that is beneficial and enhances the boundary layer fluid flow near the sheet. Cross Flow Velocity Profile (g(η)) The impact of Bu on g(η) is exhibited in Figure 9. It can be observed that as Bu rises, g(η) decreases away from
Figure 8. f '(η) vs ζ. the sheet. g(η) results from the application of a high magnetic field, g(η) does not exist in absence of magnetic field (H=0). However, when magnetic field increases, the Hall effect causes a cross flow to automatically form transversely. It also shows that as Bu rises, g(η) escalates near the sheet. From Fig. 10, the effect of H on g(η) is explored. g(η), is produced by a strong applied magnetic field. Figure 10 shows that the g(η) increases together with the magnetic field value H. The z-axis direction is where the cross flow
Sigma J Eng Nat Sci, Vol. 43, No. 4, pp. 1179−1196, August, 2025
Figure 11. g(η) vs h. g(η) is induced by the Hall effect. Since the strong magnetic field is applied along the y-axis, which causes the Hall parameter h, the Hall force is parallel to the z-axis and flow rises in the direction of the Hall force. As illustrated in Figure 11, g(η) rises as h increases. Furthermore, with extremely high values of h, the magnetic force terms assume zero, such that the g(η) profiles trend to their classical hydrodynamic values as h approaches infinity. The influence of Pr on g(η) is shown in Figure 12.
Figure 12. g(η) vs Pr. The link between the product of fluid viscosity and specific heat and thermal conductivity is established by the Prandtl number under constant pressure using a reference point. g(η) decreases with rising values of Pr because under constant pressure, the viscous forces are too strong to retain the fluid particles (atoms/molecules). Figure 13 shows the impact of Sp on the g(η). It is found that when Sp values grow, the VBL thickens and the g(η) drops. g(η) decreases when the flow resistance rises due to
Sigma J Eng Nat Sci, Vol. 43, No. 4, pp. 1179−1196, August, 2025
Figure 15. θ(η) vs Bu. the presence of porous media. Figure 14 makes it clear that g(η), drops as ζ increases away from the sheet. This is due to pressure gradient produced by positive values of ζ, and this gradient lowers g(η) away from the sheet. Whereas, it shows reverse trend near the sheet. Temperature Profile (θ(η)) The impact of Bu on θ(η) is exhibited in Figure 15. It can be observed that as Bu rises, the fluid’s temperature decreases. In contrast, the temperature escalates with the
Figure 16. θ(η) w.r.t H. rise in H as shown in Figure 16. This is due to the fact that when an electrically charged fluid is exposed to a transverse magnetic field, it gives rise to a resistive force, known as the Lorentz force. By increasing the friction between its layers, this force causes the fluid to encounter resistance. As a result, the temperature in the boundary layer rises. Figure 17 depicts the sequel of Hall parameter h on θ(η). It is seen that with the rise in h, the fluid’s temperature decreases. The sequel of Pr on θ(η) is depicted in Figure 18. It is evident from the graph that temperature θ(η) decreases
Sigma J Eng Nat Sci, Vol. 43, No. 4, pp. 1179−1196, August, 2025
Figure 19. θ(η) vs Ec. as Prandtl number Pr increases. Fluids with high Prandtl number have comparatively poor thermal conductivity, which lowers the fluid’s temperature by decreasing conduction and TBL size. The temperature distribution for various Ec values is shown in Figure 19. The temperature distribution θ(η) escalates as Ec rises. This is consistent with the idea that energy is deposited in the fluid zone due to the dissipation caused by viscosity and elastic deformation. As the quantities of thermal radiation parameter R grow, Figure 20
Figure 20. θ(η) vs R. shows that θ(η) increases, leading to an increase in the size of TBL. This is because an escalation in radiation parameter R results in a drop in the mean Roseland absorption co-efficient k, which raises the temperature profile. Radiation usually happens at high temperatures. Figure 21 exhibits the sequel of δ on θ(η). It can be noticed that as δ escalates, the fluid’s temperature escalates, as a result thickens the TBL. Figure 22. depicts the sequel of Sc on θ(η). It exhibits that as Sc increases, the dimensionless temperature θ(η) rises as well, causing the TBL to narrow.
Sigma J Eng Nat Sci, Vol. 43, No. 4, pp. 1179−1196, August, 2025
Figure 23. θ(η) vs ζ. Figure 24. θ(η) vs K. It results as an increase in molecular diffusivity causes the fluid molecule’s kinetic energy to rise as well, resulting in an increase in average speed and range of motion. The fluid›s temperature escalates as a result. It can be observed from Figure 23 that when ζ values increase, the temperature drops. This can be attributed to the fluid›s acceleration caused by positive values of ζ, which lowers the TBL. From Figure 24 it can be noticed that as K increases, the fluid’s temperature escalates.
Concentration Profile (ϕ(η)) The impact of Bu on ϕ(η). is exhibited in Figure 25. It can be perceived that as Bu rises, the fluid’s concentration declines. Figure 26 shows the effect of H on ϕ(η). The concentration field ϕ(η) rises with escalating values of H, indicating that the magnetic field regulates the fluid flow. The strong magnetic field induces the Lorentz force, which regulates the fluid flow. Figure 27 depicts the effect of h on ϕ(η).ϕ(η) is not favoured by the Hall effect, meaning that concentration
Sigma J Eng Nat Sci, Vol. 43, No. 4, pp. 1179−1196, August, 2025
decreases as Hall parameter values increase. A crossflow in the z-direction is generated by the Hall force, which somewhat but not entirely supports concentration. Due to the CBL thickening, the concentration field in Figure 28 exhibits an increasing trend against Pr. The momentum diffusivity to thermal diffusivity ratio is known as Pr (Prandtl number). In actuality, Prandtl number strengthens the viscous forces, which causes the fluid›s viscosity to naturally rise and improve concentration.
Figure 29 makes it clear that when the thermal radiation parameter N is increased, the species concentration decreases. Reduced molecular activity results from a drop in fluid temperature caused by an increase in thermal radiation. Reduced fluid temperature decreases the thermal flux›s force on the concentration flux. The molecular concentration is lowered as a result. This suggests that thermal radiation has a tendency to detard the fluid›s species concentration. Figure 30 demonstrates how the non-dimensional
Sigma J Eng Nat Sci, Vol. 43, No. 4, pp. 1179−1196, August, 2025
ϕ(η) falls with changing Sc. The ratio of the viscosity and CBL thicknesses is characterized by Sc (Schmidt number), which is a crucial parameter in the mass transfer process. The CBL is greatly diminished by rise in Sc. A tiny induced flow moves along the sheet surface as a result of the solutal buoyancy effect, which is caused by this drop in solute concentration.
Figure 32. ϕ(η) vs ζ. Figure 31 illustrates, the effect of δ on ϕ(η). It can be noticed that as δ escalates, ϕ(η) escalates. It can be observed from Figure 32 that when ζ values increase, the concentration drops. This can be attributed to the fluid›s acceleration caused by positive values of ζ, which lowers the CBL. The impact of K on ϕ(η) is depicted in Figure 33. Concentration distributions are shown to decrease as the chemical reaction intensifies. Chemical reactions occur physically with several disruptions. This as a result leads to high molecular mobility, which diminishes the concentration distributions of fluid flow by increasing the transport phenomena. Table 3 shows the sequel of various fluid flow parameters on the coefficient of skin friction along x-axis f ''(0) It can be noticed that the magnitude of f ''(0) escalates with the escalating values of H, Sp, Pr, Sc and K . Whereas, decreases with the rise in Bu, h, δ, ζ, Ec and R. Table 3 also shows the sequel of various fluid flow parameters on the coefficient of skin friction along z-axis g'(0). It can be noticed that magnitude of g'(0) accelerates with the rise in Bu, H, ζ, Ec, R and δ. Whereas, decreases with the rise in h, Sp, Pr, Sc and K. Table 4 shows the sequel of various fluid flow parameters on the rate of heat transfer θ'(0). It can be contemplated that the magnitude of θ'(0) escalates with the rise in Bu, h, ζ and Pr. Whereas, decreases with the rise in H, Ec, Sp, R, Sc, δ and K. Table 4 also shows the sequel of various fluid flow parameters on rate of mass transfer ϕ'(0) It is demonstrated that magnitude of ϕ'(0) escalates with the rise in Bu, h, ζ, Ec, R, Pr, δ and K. Whereas, decreases with the rise in H, Sp and Pr.
Sigma J Eng Nat Sci, Vol. 43, No. 4, pp. 1179−1196, August, 2025
Table 3. Variation of skin friction coefficients f ''(0) and g'(0) for the various parameters Bu
Conclusion
A viscous, incompressible, electrically conducting fluid is transferred through an exponentially stretched sheet submerged in a porous material by MHD mixed convection flow, which transfers mass and heat. The impacts of Hall current, Heat source, Joule heating, radiation and chemical reaction are investigated. The bvp4c technique was utilized to provide numerical discussions for the governing equations. A thorough collection of graphics is provided, together with a discussion of how these variables depend on various physical parameters, for the fluid tangential, cross flow velocity, fluid temperature and species concentration.
The following are noteworthy findings: The tangential velocity escalates with the rise in Bu, h and ζ. Whereas, it detards with the rise in H, Pr and Sp. The cross flow velocity escalates with the rise in H, h and Pr. Whereas, it detards with the rise in Bu, ζ and Sp. The temperature profile escalates with the rise in H, Ec, R, δ, Sc and K. Whereas, it detards with the rise in Bu, h, Pr and ζ. The concentration profile escalates with the rise in H and Pr. Whereas, it detards with the rise in Bu, h, R, Sc, δ, ζ and K. The magnitude of local skin friction coefficient along x- direction escalates with the rise in H, Sp, Pr, Sc and K.
Sigma J Eng Nat Sci, Vol. 43, No. 4, pp. 1179−1196, August, 2025
Table 4. Variation of rate of heat transfer and mass transfer coefficients i.e., θ'(0) and ϕ'(0) for the various parameters Bu
Whereas, it decreases with the rise in Bu, h, δ, ζ, Ec and R. · The magnitude of local skin friction coefficient along zdirection escalates with the rise in Bu, ζ, h, Ec, Sp, R and δ. Whereas, it decreases with the rise in H, Pr, Sc and K. · The magnitude of rate of heat transfer increases with the rise in Bu, h, ζ and Pr. Whereas, it decreases with the rise in H, Ec, Sp, R, Sc, δ and K. · The magnitude of rate of mass transfer increases as Bu, h, ζ, Ec, R, Pr, δ and K escalates. Whereas, it decreases with the rise in H, Sp and Pr. Future work on the problem’s irregular behavior and its relationship to thermophoresis and activation energy should
be the main focus of the present inquiry. Furthermore, the study of entropy optimization in the magnetohydrodynamic flow of viscous liquids is of tremendous theoretical and technical significance. Metallic coating, crystal growth, electromagnetic pumps, power generators, MHD accelerators, and reactor cooling are just a few examples of the numerous industrial and technological applications where significant effects of hall current, thermal radiation, heat source, and chemical reaction parameters can be found on flow profiles. The production of paper, suspensions, heat exchanger technology, material processing, drying, and surface evaporation of water bodies are only a few of the possible uses for the study’s findings. Furthermore, we have taken into account a
Sigma J Eng Nat Sci, Vol. 43, No. 4, pp. 1179−1196, August, 2025
conventional fluid instead of a hybrid nanofluid or a nanofluid that has a higher mass transfer and heat efficiency than the fluid under consideration. Since we have left them openended for our research’s future direction.
Nomenclature
Cp K* f g Ec qr k1 k H Pr Sp R Gr D Sc Q h Re Nu Sh K Bu t T g1 C u v w
Specific heat at constant pressure, (J/kgK) Rate of first order chemical reaction Dimensionless stream function Cross flow velocity Eckert number radiative heat flux, (W/m2) Permeability of porous medium, (m2) Absorption coefficient Magnetic parameter Prandlt number Porosity parameter Radiation parameter Grashoff number mass diffusivity Schmidt number Heat source Hall parameter Reynolds number Nusselt number Sherwood number Chemical reaction parameter Buoyancy ratio parameter Dimensionless time, (K) fluid's temperature, (K) acceleration due to gravity, Species concentration, (kg/m3) fluid's velocity through x-axis, (m/s) fluid's velocity through y-axis, (m/s) fluid's velocity through z-axis, (m/s)
Greek Symbols ρ fluid’s density, (kg/m3) μ fluid's viscosity, (Pa s) σ Electrical conductivity, (1/Ωm) η Dimensionless similarity variable ν Kinematic viscosity, (m2 /s) κ Thermal conductivity, (W/mK) α Thermal diffusivity, (m2 /s) Coefficient of thermal expansion βt Coefficient of spatial expansion βc ψ Stream function τ Skin friction θ Dimensionless temperature δ Heat souce parameter ζ mixed convection parameter ϕ Dimensionless concentration Superscript ' Derivative w.r.t η
Subscript w Properties at the plate ∞ Free stream condition Abbreviation VBL Velocity boundary layer TBL Temperature boundary layer CBL Concentration boundary layer
Acknowlegdement
The Department of Mathematics, Cotton University was the site of this research. All of the reviewers’ insightful recommendations and remarks are also appreciated by the authors.
Data Availability Statement
The authors confirm that the data that supports the findings of this study are available within the article. Raw data that support the finding of this study are available from the corresponding author, upon reasonable request.
Conflict Of Interest
The author declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.
Ethics
There are no ethical issues with the publication of this manuscript.
References
- Kameswaran PK, Narayana M, Sibanda P, Makanda G. On radiation effects on hydromagnetic Newtonian MHD flow past a rotating inclined plate with heat liquid flow due to an exponential stretching sheet. transfer and mass diffusion. Int J Sci Res Math Stat Bound Value Probl 2012;2012:1–16. [CrossRef] Sci 2018;5:253–258.
- Mandal IC, Mukhopadhyay S. Heat transfer anal- ysis for fluid flow over an exponentially stretch- on MHD flow past an inclined plate with variable ing porous sheet with surface heat flux in porous wall temperature and mass diffusion in the presence medium. Ain Shams Eng J 2013;4:103–110. [CrossRef] of Hall current. Appl Appl Math 2018;13:31. [CrossRef]
- Singh V, Agarwal S. MHD flow and heat transfer for Maxwell fluid over an exponentially stretching and mass transfer flow of nanofluid through porous sheet with variable thermal conductivity in porous media over inclined plate with chemical reaction. medium. Therm Sci 2014;18:599–615. [CrossRef] Multidiscip Model Mater Struct 2021;17:317–336.
- Jaber KK. Effect of Hall currents and variable fluid [CrossRef] properties on MHD flow past stretching vertical [22] Goud BS, Bindu P, Srilatha P, Krishna YH. The Joule plate by the presence of radiation. J Appl Math Phys heating effect on MHD natural convective fluid flow 2014;2:888–902. [CrossRef] in a permeable medium over a semi-infinite inclined
- Khidir A, Sibanda P. Effect of temperature-depen- vertical plate in the presence of the chemical reaction. dent viscosity on MHD mixed convective flow from IOP Conf Ser Mater Sci Eng 2020;993:012111. [CrossRef] an exponentially stretching surface in porous media [23] Sandhya A, Reddy GR, Deekshitulu GVSR. Steady with cross-diffusion. Spec Top Rev Porous Media on MHD heat and mass transfer flow of an inclined 2014;5:157–170. [CrossRef] porous plate in the presence of radiation and chem-
- Gorla RSR, Gireesha BJ, Singh B. MHD flow and ical reaction. J Phys Conf Ser 2019;1344:012002. heat transfer of dusty nanofluid embedded in porous [CrossRef] medium over an exponentially stretching sheet. J [24] Bharathi V, Vijayaragavan R, Prakash J. Heat and Nanofluids 2015;4:449–460. [CrossRef] mass transfer effect of a magnetohydrodynamic
- Hayat T, Shafiq A, Alsaedi A, Shahzad SA. Unsteady Casson fluid flow in the presence of inclined plate. MHD flow over exponentially stretching sheet with Indian J Pure Appl Phys 2021;59:28–39. slip conditions. Appl Math Mech 2016;37:193–208. [25] Dwivedi N, Singh AK. Transient free convective [CrossRef] hydromagnetic flow in an infinite vertical cylinder
- Eid MR. Chemical reaction effect on MHD bound- with Hall current and heat source/sink. Heat Transf ary-layer flow of two-phase nanofluid model over 2020;49:4091–4108. [CrossRef] an exponentially stretching sheet with a heat gener- [26] Dwivedi N, Kumar Singh A. Influence of Hall cur- ation. J Mol Liq 2016;220:718–725. [CrossRef] rent on hydromagnetic natural convective flow
- Srinivasacharya D, Jagadeeshwar P. MHD flow between two vertical concentric cylinders in pres- with Hall current and Joule heating effects over ence of heat source/sink. Heat Transf 2020;49:1402– an exponentially stretching sheet. Nonlinear Eng 1417. [CrossRef] 2017;6:101–114. [CrossRef] [27] Jeevitha S, Chitra M, Rushi Kumar B. MHD flow in
- Govindarajan A, Vijayalakshmi R. Hall current a rotating vertical cone through a porous medium. effect in MHD oscillatory couple stress dusty fluid Heat Transf 2023;52:2165–2185. [CrossRef] through an inclined saturated permeable channel. [28] Tasnim S, Mitra A, Saha H, Islam MQ, Saha S. MHD Int J Eng Technol 2018;7:801–805. [CrossRef] conjugate natural convection and entropy genera-
- VeeraKrishna M, Subba Reddy G, Chamkha AJ. tion of a nanofluid filled square enclosure with mul- Hall effects on unsteady MHD oscillatory free con- tiple heat-generating elements in the presence of vective flow of second grade fluid through porous Joule heating. Results Eng 2023;17:100993. [CrossRef] medium between two vertical plates. Phys Fluids [29] Zainal NA, Nazar R, Naganthran K, Pop I. Viscous 2018;30:023106. [CrossRef] dissipation and MHD hybrid nanofluid flow towards
- VeeraKrishna M, Chamkha AJ. Hall effects on an exponentially stretching/shrinking surface. unsteady MHD flow of second grade fluid through Neural Comput Appl 2021;33:11285–11295. [CrossRef] porous medium with ramped wall temperature [30] Abbas S, Nazar M, Nisa ZU, Amjad M, Din SME, and ramped surface concentration. Phys Fluids Alanzi AM. Heat and mass transfer analysis of MHD 2018;30:053101. [CrossRef] Jeffrey fluid over a vertical plate with CPC fractional
- Sarma D, Pandit KK. Effects of Hall current, rota- derivative. Symmetry 2022;14:2491. [CrossRef] tion and Soret effects on MHD free convection heat [31] Swarnalathamma BV, Babu DP, Krishna MV. Hall and mass transfer flow past an accelerated vertical current effects on MHD flow of chemically reacting plate through a porous medium. Ain Shams Eng J fluid through a porous medium with heat source. 2018;9:631–646. [CrossRef] Heat Transf 2022;51:4123–4142. [CrossRef] 1196 Sigma J Eng Nat Sci, Vol. 43, No. 4, pp. 1179−1196, August, 2025
- Sheri SR, Peesu M, Mamidi Narsimha R. Hall cur- rent, chemical reaction, and radiation results on media. New York: Springer; 1991. p. 52. [CrossRef] transient magnetohydrodynamic flow past an [40] Shaw S, Guerrero D. Exploring the utilization of inclined plate: FEM. Heat Transf 2022;51:1876– Newtonian fluids in heat transfer applications. J 1899. [CrossRef] Fluid Flow 2023;10:120–130. [CrossRef]
- Pati AK, Misra A, Mishra SK, Mishra S, Sahu R, Panda S. Computational modelling of heat and mass semiconductors. Thin Solid Films 1978;52:379–395. transfer optimization in copper water nanofluid [CrossRef] flow with nanoparticle ionization. JP J Heat Mass Transf 2023;31:1–18. [CrossRef] [42] Cramer KR, Pai SI. Magnetofluid dynamics for engi-
- Bishnoi J, Kumar S. Hall current induced MHD con- neers and applied physicists. New York: McGraw- vection of a nanofluid subjected to the heat and salt Hill Book Company; 1973. [CrossRef] diffusion. Fluid Dyn 2023;58:1169–1182. [CrossRef] [43] Shampine LF, Kierzenka J, Reichelt MW. Solving
- Venthan SM, Kumar PS, Kumar SS, Sudarsan S, boundary value problems for ordinary differential Rangasamy G. A computational study of the impact equations in MATLAB with bvp4c. Tutorial Notes of fluid flow characteristics on convective heat transfer 2000;1–27. with Hall current using the MHD non-Newtonian fluid [44] Wahid NS, Arifin NM, Khashi’ie NS, Pop I, Bachok model. Chem Eng Res Des 2024;203:789–799. [CrossRef] N, Hafidzuddin MEH. Unsteady mixed convective
- Ali A, Khan HS, Noor I, Pasha AA, Irshad K, Al stagnation point flow of hybrid nanofluid in porous Mesfer MK, et al. Hall effects and Cattaneo–Christov medium. Neural Comput Appl 2022;34:14699– heat flux on MHD flow of hybrid nanofluid over a 14715. [CrossRef] varying thickness stretching surface. Mod Phys Lett
- Magyari E, Keller B. Heat and mass transfer in the B 2024;38:2450130. [CrossRef]
- Elsaid EM, Abdel‐Wahed MS. Thermal evalua- boundary layers on an exponentially stretching con- tion of MHD boundary‐layer flow of hybridity tinuous surface. J Phys D Appl Phys 1999;32:577– nanofluid via a 3D sinusoidal cylinder. ZAMM 585. [CrossRef] 2024;104:e202300186. [CrossRef] [46] Ishak A. MHD boundary layer flow due to an expo-
- Nield DA, Bejan A. Convection in porous media. nentially stretching sheet with radiation effect. Sains 2nd ed. New York: Springer; 1992. p. 12. [CrossRef] Malays 2011;40:391–395.
Share and Cite
SARMA, A.K.; SARMA, D.; MUDOI, S. MHD flow and effects of hall current joule heating and heat source across a chemically reacting and. Sigma Journal of Engineering and Natural Sciences 2025, Vol. 43, pp. 1179-1196. https://doi.org/10.14744/sigma.2024.00110

