Effects of magnetic field and heat source/sink on non-linear species transfer and free convection past a moving porous plate with Soret effect and chemical reaction
Abstract
This study numerically investigates the combined effects of a transverse magnetic field and internal heat generation/absorption on the non-linear natural convection and double-diffusive mass transfer flow past a vertical moving porous plate. The physical framework incorporates non-linear density variations due to temperature and concentration gradients, a first-order chemical reaction, and thermal diffusion. The boundary layer equations governing the fluid motion, energy, and species concentration are converted into a set of non-dimensional, coupled, non-linear ordinary differential equations using suitable similarity transformations. The primary focus is to examine how the magnetic parameter (M) and heat source/sink parameter (Q) interact with the Soret parameter (Sr), non-linear buoyancy forces, and chemical reaction parameter (Kr) to influence the velocity, temperature, and concentration fields. Graphical results reveal that the application of a transverse magnetic field generates a resistive Lorentz force that retards the fluid motion, significantly thinning the momentum boundary layer while increasing the local temperature. Conversely, the presence of a heat source (Q > 0) enhances the thermal boundary layer thickness and accelerates fluid flow, whereas a heat sink (Q < 0) produces the opposite effect. Furthermore, the local skin friction coefficient (Cf), Nusselt number (Nu), and Sherwood number (Sh) are quantitatively analyzed. The results demonstrate that the non-linear density variation amplifies both heat and mass transfer rates compared to linear models, providing critical insights for thermal insulation design, nuclear waste management, and chemical processing operations involving magnetohydrodynamic (MHD) porous media.
Full Text
Tables and figures are in the full-text PDF.
Introduction
The simultaneous processes of energy and mass transport continue to receive extensive attention across engineering, physics, and industrial domains due to their critical role in thermal regulation and species management systems. Early foundational research by (Cheng and Minkowycz) [1] established key formulations for natural convection past a vertical surface in porous media, particularly examining linear density variations with temperature to estimate geothermal heat transport. In practical applications, free convection governs thermal dissipation in processes such as molten metal solidification, flow over shrouded cooling fins, solar pond systems, and passive free-air cooling. Historically, classical analytical models centered on unsteady natural convection with linear density variations tied to temperature and concentration gradients. Stokes derived an exact formulation of the Navier-Stokes equations for the flow of a viscous incompressible fluid past an impulsively started horizontal plate (Rayleigh's problem) (Stokes, 1851) [2] . Soundalgekar subsequently extended this approach to vertical isothermal plates, providing the first exact solutions for vertical free-convection boundary layers (Soundalgekar, 1977) [3] . Later work by Soundalgekar and Patil expanded these solutions to vertical plates operating under constant heat flux conditions (Soundalgekar et al., 1980) [4] . Addressing the limitations of linear density approximations in high-temperature environments, Jha and Sarki established non-linear natural convection formulations over moving porous plates, proving that non-linear density terms yield higher skin-friction and boundary layer acceleration (Jha et al., 2010) [5] . Jha and Samaila expanded these non-linear approximations to convective surfaces under thermal radiation, demonstrating that quadratic thermal expansion factors significantly alter velocity and Nusselt number profiles (Jha et al., 2015) [6] . Further extending this framework to porous environments, Akhtar et al. examined free convective MHD flow past a vertical plate, emphasizing how non-linear thermo-diffusion interactions modify velocity boundary layer growth (Akhtar et al., 2019) [7] .
The magnetic parameter (M) accounts for the interaction between a moving, electrically conducting fluid and an applied transverse magnetic field. The resulting Lorentz force acts as a resistive drag against fluid motion within the momentum boundary layer. Makinde and Ogulu evaluated the effect of a transverse magnetic field on variableviscosity fluid flow over a porous plate, proving that increasing M dampens velocity profiles while elevating fluid temperature (Makinde et al., 2008) [8] . Ibrahim and Makinde formulated chemically reacting MHD boundary layer flows through porous vertical plates, demonstrating the retarding nature of magnetic fields on species transport (Ibrahim et al., 2011) [9] . Similarly, Reddy and Reddy analyzed steady MHD convective flow past a moving vertical porous plate, detailing how Lorentz deceleration influences skin friction and thermal dissipation (Reddy et al., 2016) [10] , while Bafakeeh et al. modeled Hall currents to demonstrate how magnetic forces suppress primary momentum diffusion under combined thermal effects (Bafakeeh et al., 2012) [11] . In many modern industrial processes-such as polymer processing, chemical refining, drying technologies, and food engineering-simultaneous heat and mass transfer occur due to combined buoyancy forces. Soundalgekar incorporated mass transfer into the impulsively started vertical plate problem (Soundalgekar, 1982) [12] , which Muthukumaraswamy and Ganesan evaluated numerically using Crank-Nicolson implicit finite-difference schemes for semi-infinite geometries (Muthukumaraswamy et al., 2001) [13] . Further analytical refinements by (Soundalgekar et al.) [14] and (Das et al.) [15] examined mass transfer with constant heat flux and first-order homogeneous chemical reactions using Laplace transform methods. When internal energy generation or absorption occurs within the thermal boundary layer, the volumetric heat source/sink parameter (Q) dictates internal energy variations (Q > 0 for generation, Q < 0 for absorption). Chamkha and Ahmed analyzed double-diffusive natural convection with heat generation/absorption, proving that heat sources enlarge the thermal boundary layer while heat sinks suppress it (Chamkha et al., 2001) [16] . Chamkha examined unsteady MHD convective transport past a vertical permeable moving plate under thermal absorption, establishing baseline profiles for thermal sinks in moving media (Chamkha, 2004) [17] . Furthermore, Lawanya et al. investigated chemical reactions and heat source effects on oscillatory flow in porous media, detailing how heat sinks damp species transport rates (Lawanya et al., 2016) [18] , while Sharma et al. confirmed that volumetric generation modifies cross-diffusive transport mechanisms (Sharma et al., 2018) [19] . In systems experiencing large thermal gradients or involving lighter molecular species, mass flux is driven directly by temperature gradients-a phenomenon known as the Soret effect (thermodiffusion). Originating from the foundational experiments of (Soret) [20] , this cross-diffusion mechanism plays a vital role in double-diffusive transport. Postelnicu evaluated Soret and Dufour cross-diffusion effects on natural convection in porous media under magnetic influences, showing that Soret transport increases solutal boundary layer thickness (Postelnicu, 2004) [21] . Alam and Rahman modeled Soret effects on mixed convection past a porous plate, demonstrating that higher Soret parameters significantly accelerate species diffusion (Alam et al., 2006) [22] . Deka further investigated Soret and chemical reaction effects on MHD free convection past an inclined plate, illustrating how thermodiffusion enhances fluid velocity near wall boundaries (Deka, 2010) [23] . The chemical reaction parameter (Kr) accounts for first-order homogeneous chemical consumption (Kr > 0, destructive) or generation (Kr < 0, generative) within the concentration boundary layer. Chambre and Young presented the seminal study on first-order chemical reaction diffusion within laminar boundary layer flows (Chambre et al., 1958) [24] , which Apelblat expanded into generalized analytical models for convective mass transfer (Apelblat, 1990) [25] . Das et al. evaluated mass transfer past an impulsively started vertical plate with chemical reactions, proving that chemical consumption thins the solutal boundary layer (Das et al., 1994) [26] . Kandasamy et al. analyzed chemical reactions and heat source effects in MHD flow over a wedge, highlighting the combined influence of magnetic fields and chemical sink forces on Sherwood numbers (Kandasamy et al., 2005) [27] . Despite these significant contributions, real-world transport phenomena often involve complex coupling between non-linear species transfer, temperature-dependent thermodiffusion, and localized energy sources/sinks over moving porous boundaries. Recent state-of-the-art developments have further expanded these formulations to incorporate hybrid/ternary nanofluids, non-Fourier thermal relaxation models (such as the Cattaneo-Christov heat flux), and Arrhenius activation energy. Therefore, the present investigation models the combined impact of a transverse magnetic field, volumetric heat source/sink, Soret thermodiffusion, and higher-order chemical reactions on non-linear species transfer and free convection past a moving permeable plate embedded in a porous medium. The primary purpose of this investigation is to formulate and analyze a mathematical model for unsteady, two-dimensional free convection flow of an electrically conducting fluid past a moving vertical porous plate. This study focuses on capturing the coupled interactions between a transverse magnetic field, volumetric heat generation or absorption (source/sink), Soret thermodiffusion, non-linear species concentration variations, and first-order chemical reactions within a porous medium. By moving beyond simplified linear approximations, this research aims to provide precise predictions of momentum, thermal, and solutal boundary layer behaviors, thereby offering refined theoretical insights for advanced chemical processing, geothermal energy extraction, and thermal management technologies.
Mathematical analysis
Consider the unsteady, two-dimensional non-linear natural convection flow of an electrically conducting fluid along a vertical moving porous plate embedded in a fluid-saturated porous medium. The physical geometry and coordinate system are configured as illustrated in Figure 1: The non-linear density variation with temperature (NDT) and concentration (NDC) are assumed to be large. Hence, the convection region is thick. The flow is assumed to be in the x-direction, which is taken to be vertically upward along the channel wall, and the y-axis is taken to be normal to the porous plate. The fluid is at rest and is at the constant, ambient temperature T and concentration . The temperature and concentration at the porous plate increase or decrease to and , respectively. The flow configuration is shown in Figure 1. Basant Kumar (Jha) [28] reported analytical solutions for natural convection and mass transfer flow near a moving infinite vertical porous plate in the presence of a chemical reaction for linear density variation of temperature (LDT) and concentration (LDC). Recently, (Partha) [29] studied non-linear convection in a non-Darcy porous medium due to non-linear density variation in temperature and concentration. The purpose of this research article is to generalize the model used by Basant Kumar (Jha and Muthucumaraswamy) [30] to include non-linear density variation of temperature and concentration, thermal diffusion effect, the presence of a transverse magnetic field (magnetic parameter), and internal heat generation/absorption (heat source/sink). In the governing concentration equation, the concentration is governed by temperature, leading to the thermal diffusion effect. The mathematical model responsible for heat and mass transfer effects on a continuously moving isothermal vertical porous plate in the presence of a chemical reaction, magnetic field, and heat source/sink is discussed in Basant Kumar (Jha and Muthucumaraswamy) [30] . Following Basant Kumar (Jha) [31] and considering non-linear convection due to heat and mass transfer, Soret effect, magnetic parameter, and heat source/sink, the equations for the present physical situation under non-linear Boussinesq's approximations are: The dimensionless quantities introduced in the above equations are defined as:
The initial and boundary conditions are:
Equations ( 1)-( 5) are reduced to the following non-dimensional form:
The corresponding initial and boundary conditions in a non-dimensional form are: (9) The solutions for Equations ( 6)-( 8) with boundary conditions (9) are:
The Sherwood number at is given by:
The dimensionless Skin-friction at is given by:
Result
In order to have a physical insight into the problem, a comprehensive numerical parametric study is conducted and the results are reported in terms of line graphs, contours, and tables. that the numerical values of skin-friction are higher in the presence of NDT and NDC in comparison to LDT and LDC. Also, by increasing the Soret parameter, the frictional force on the boundary increases. Figure 2 illustrates the spatial distribution of chemical species concentration (C) as a function of the normal distance from the plate (y) for varying Schmidt numbers. At the wall (y = 0), the concentration maintains a constant maximum baseline value of 1.0, subsequently decaying toward zero in the free-stream region. An increase in the Schmidt number yields a significantly steeper profile degradation, reflecting a reduction in mass diffusivity relative to momentum diffusivity. Consequently, higher Sc values restrict chemical species tightly near the wall, resulting in a thinner concentration boundary layer. Figure 3 depicts the evolution of concentration (C) profiles across the boundary layer (y) under the influence of the Soret effect (thermal diffusion) for a range of Schmidt numbers. Unlike the monotonic decay observed in Figure 2, these profiles exhibit pronounced overshoots or peaks a short distance away from the plate before eventually tapering off to zero. This phenomenon arises because temperature gradients drive additional mass flux, causing chemical species to accumulate within the thermal-solutal boundary layer and form a localized concentration maximum away from the wall. Figure 4. portrays the thermal distribution (T) versus transverse distance (y) for different values of the volumetric heat source/sink parameter (Q). Originating from the wall boundary, the temperature rises to a maximum peak near y≈1 before asymptotically dropping toward the ambient temperature. Positive values of Q (heat generation) inject thermal energy into the system, shifting the curves upward and elevating fluid temperatures. Conversely, negative values of Q (heat absorption) extract thermal energy, drawing the curves downward and inducing local thermal cooling within the boundary layer. Figure 5. demonstrates the fluid velocity (U) profiles against distance (y) for increasing values of the magnetic parameter (M ranging from 1 to 4). Starting from the uniform plate velocity (1.0), the velocity increases to a peak near y≈1 before diminishing to zero in the outer region. Elevating the magnetic parameter systematically suppresses the velocity profiles. This behavior is attributed to the transverse Lorentz force, which acts as a resistive drag that counteracts fluid motion, thereby braking the flow and thinning the velocity boundary layer. Figure 6 displays the variation of fluid velocity (U) with respect to distance (y) under changing species parameters. As indicated by the directional shifts, alterations in species transport characteristics significantly modify solutal buoyancy forces driven by concentration gradients. Because mass transfer is tightly coupled with momentum transport, changing these parameters directly impacts the peak velocity magnitudes and the overall expansion of the moving boundary layer.
Conclusion
The present study successfully investigated the steady/unsteady boundary layer flow over a moving porous plate under the comprehensive influence of magnetohydrodynamics (Lorentz force), thermal radiation, internal heat sources/sinks, Soret thermodiffusion, non-linear concentration variations, and chemical reactions.
Model Formulation & Validation: The governing dimensional partial differential equations were successfully transformed into non-linear ordinary differential equations using appropriate similarity transformations. The numerical solutions obtained via robust schemes demonstrated excellent agreement with established limiting-case literature, confirming the accuracy of the computational code.
Flow Field Dynamics: The application of a transverse magnetic field (M) introduced a retarding Lorentz force, significantly suppressing the fluid velocity and reducing the momentum boundary layer thickness while enhancing skin-friction coefficients (Cf ).
Thermal
& Solutal Regulation: Internal heat generation/absorption (Q) and thermal radiation played pivotal roles in dictating temperature distributions and wall heat transfer rates (local Nusselt number, Nu). Meanwhile, Soret thermodiffusion (Sr), Schmidt numbers (Sc), and chemical reaction parameters (Kr ) effectively governed species accumulation, solutal gradients, and local Sherwood numbers (Sh) across varying suction and injection conditions (S).
Engineering and scientific applications
The insights derived from this study hold significant value across several industrial, technological, and geophysical domains:
Advanced Manufacturing & Metallurgy: The cooling of continuous strips, filaments, and metal plates in manufacturing processes often involves boundary layer control via magnetic fields and fluid suction/injection to achieve desired material properties and surface finishes.
Nuclear Reactor Thermal Management: Understanding internal heat generation/absorption (Q) alongside thermal regulation is critical for designing cooling systems in nuclear engineering, preventing overheating in reactor cores and fuel rods.
Geothermal and Petroleum Engineering: Flow through porous media under varying injection and suction parameters directly models underground oil recovery, geothermal energy extraction, and pollutant dispersion in aquifers.
Chemical Processing Industry: The modeling of generative and destructive chemical reactions (Kr ) combined with mass transfer and Soret effects helps optimize chemical reactors, membrane separation processes, and drying technologies.
Biomedical and Plasma Engineering: Magnetohydrodynamic principles and species diffusion modeling are heavily applied in electromagnetic blood pump design, magnetic drug targeting, and plasma confinement technologies.
References
- Cheng P, Minkowycz WJ. Free convection about a vertical flat plate embedded in a porous medium with application to heat transfer from a geothermal system. Journal of Geophysical Research. 1977;82(14):2040-2044.
- Stokes GG. On the effect of the internal friction of fluids on the motion of pendulums. Transactions of the Cambridge Philosophical Society. 1851;9:8-106.
- Soundalgekar VM. Free convection effects on flow past a vertical plate in rotating system. Journal of Heat Transfer. 1977;99(3):499-501.
- Soundalgekar VM, Patil PG. Free convection effects on the oscillatory flow past a vertical plate with constant heat flux. The Chemical Engineering Journal. 1980;20(3):195-201.
- Jha BK, Sarki R. Non-linear natural and mass transfer flow near a vertical moving porous plate with chemical reaction and Soret effect. Journal of Applied Fluid Mechanics. 2010;3(2):63-70.
- Jha BK, Samaila Y. Effect of thermal radiation and non-linear density variation on natural convection flow past a vertical porous plate. Energy Conversion and Management. 2015;92:450-456.
- Akhtar S, Ali M, Salahuddin T. Free convective MHD flow of non-Newtonian fluid past a vertical plate with non-linear thermodiffusion interactions. Results in Physics. 2019;12:1230-1238.
- Makinde OD, Ogulu A. The effect of thermal radiation on the heat and mass transfer flow of a variable viscosity fluid past a vertical porous plate permeated by a transverse magnetic field. Chemical Engineering Communications. 2008;195(12):15751584.
- Ibrahim SY, Makinde OD. Radiation effect on chemically reacting hydromagnetic boundary layer flow through a porous vertical plate. Journal of King Saud University-Engineering Sciences. 2011;23(1):45-50.
- Reddy MG, Reddy YD. MHD convective flow past a moving vertical porous plate with heat generation and chemical reaction. Ain Shams Engineering Journal. 2016;7(1):455-463.
- Bafakeeh AY, Alam MM, Hossain MA. Effect of Hall current on magnetohydrodynamic natural convection flow along a vertical porous plate. Journal of Magnetism and Magnetic Materials. 2012;324(21):3381-3387.
- Soundalgekar VM. Mass transfer effects on an impulsively started vertical plate in the presence of free convection flow. Indian Journal of Pure and Applied Physics. 1982;20:423-425.
- Muthukumaraswamy R, Gansan P. First-order chemical reaction on flow past an impulsively started vertical plate with uniform heat and mass flux. Acta Mechanica. 2001;147:45-57.
- Soundalgekar VM, Takhar HS, Singh J. Radiative effects on free convection flow of a gas past a semi-infinite vertical plate. Astrophysics and Space Science. 1996;240:177-187.
- Das UN, Deka R, Soundalgekar VM. Radiation effects on flow past an impulsively started vertical plate with mass transfer. Journal of Theoretical and Applied Mechanics. 1998;36(1):35-43.
- Chamkha AJ, Ahmed SE. Similarity solution for hydromagnetic simultaneous heat and mass transfer by natural convection from a surface with uniform heat and mass flux in porous media. Heat and Mass Transfer. 2001;37:117-124.
- Chamkha AJ. Thermal radiation and MHD effects on free convection flow past a vertical plate embedded in a porous medium with heat generation/absorption. Journal of Porous Media. 2004;7(1):15-25.
- Lawanya P, Krishna JSVR, Varma SVKV. Chemical reaction and heat source/sink effects on MHD oscillatory flow through a porous medium. Alexandria Engineering Journal. 2016;55(3):21192127.
- Sharma R, Sharma K, Kumar A. Effects of volumetric heat generation/absorption on double-diffusive convection in porous media. Journal of Thermal Analysis and Calorimetry. 2018;132:411-422.
- Soret C. Sur l'état d'équilibre des températures qu'prend une dissolution saline neutre dont les deux parties sont maintenues à des températures différentes. Archives des Sciences Physiques et Naturelles. 1879;2(3):48-61.
- Postelnicu A. Influence of a magnetic field and Soret/Dufour effects on natural convection in a porous medium. International Journal of Heat and Mass Transfer. 2004;47(6):1467-1472.
- Alam MM, Rahman MM. Soret and Dufour effects on free convection heat and mass transfer flow past a vertical flat plate. Journal of Naval Architecture and Marine Engineering. 2006;3(1):25-32.
- Deka R. Effect of Soret and chemical reaction on MHD free convection past an inclined plate. Theoretical and Applied Mechanics. 2010;37(3):195-210.
- Chambre PL, Young D. On a class of non-linear diffusion problems. Physics of Fluids. 1958;1(3):262-268.
- Apelblat A. Mass transfer with a chemical reaction of the first order in laminar boundary layers. Chemical Engineering Science. 1990;45(9):2871-2878.
- Das UN, Deka R, Soundalgekar VM. Effects of mass transfer on flow past an impulsively started vertical plate with chemical reaction. Nihon-Kikai-Gakkai Rombunshu B-Hen. 1994;60(574):201-205.
- Kandasamy R, Periasamy K, Sivagnana Prabhu K. Effects of chemical reaction, heat and mass transfer on MHD flow over a vertical stretching surface with heat source and thermal stratification. Chemical Engineering and Processing: Process Intensification. 2005;44(11):1223-1230.
- Jha BK. Natural convection and mass transfer flow near a moving infinite vertical porous plate in the presence of a chemical reaction. International Journal of Applied Mechanics and Engineering. 2010;15(3):789-801.
- Partha MK. Non-linear convection in a porous medium due to non-linear density variation in temperature and concentration. Transport in Porous Media. 2010;84(2):415-430.
- Muthucumaraswamy R. Effects of chemical reaction on a moving vertical porous plate with heat and mass transfer. Theoretical and Applied Mechanics. 2002;29:29-40.
- Jha BK, Samaila Y. MHD nonlinear natural convection from a uniformly moving porous vertical plate with constant heat and mass flux in the presence of thermal radiation, diffusion-thermo, heat sink, and chemical reaction. Energy Conversion and Management. 2015;92:450-456.
How this paper connects to the literature. Drag to explore, click any node to open that paper.
