Abstract

This study aims to bridge the gap by conducting a numerical analysis of Maxwell fluid behaviour on a perpendicular plate within a porous medium, considering both chemical reaction and heat generation. The investigation also encompasses the study of energy and mass transfer within magnetohydrodynamic (MHD) Maxwell fluids. We utilise a transformation technique employing similarity variables to address the challenge posed by the nonlinear partial differential equations (PDEs). These transformed equations are subsequently solved via the bvp4c solver in MATLAB. The obtained results exhibit a high degree of agreement with the previously published work. The study systematically explores the influence of chemical reaction, energy generation, and Deborah number parameters on temperature and velocity, as well as concentration, presenting the outcomes graphically. In addition, we calculate local Sherwood numbers, Nusselt numbers, and skin friction coefficients to assess the impact of chemical reactions. Our findings notably indicate that Sherwood numbers and skin friction coefficients increase with higher levels of chemical reaction, while local Nusselt numbers decrease as chemical reactions become more pronounced. By studying Maxwell fluid flow over a perpendicular plate with chemical reactions, this research contributes to optimizing processes, enhancing product quality, and providing deeper insights into the behaviour of complex fluids in real-world scenarios.

1. Introduction

Chemical reactions are fundamental processes in various natural and industrial systems, transforming substances into new products. These reactions often involve the issue or absorption of energy through heat. When these chemical reactions occur in fluid flow over a porous vertical plate, they give rise to a complex interplay between the fluid dynamics and the heat generation or absorption effects. In this context, the Maxwell fluid flow refers to the flow of a viscoelastic fluid that exhibits non-Newtonian behaviour. Maxwell fluids are characterized by a time-dependent stress response to deformation, making their flow behaviour distinct from Newtonian fluids such as water or air. When such a Maxwell fluid interacts with a porous vertical plate, it introduces unique challenges and opportunities for understanding the system’s dynamics. The generation of heat due to chemical reactions adds complexity to the Maxwell fluid flow over the permeable vertical plate. This heat generation can have significant implications for the system’s overall behaviour. Depending on the nature and rate of the chemical reactions, the generated heat can either enhance or hinder fluid flow and transfer. Understanding these effects is crucial in various engineering and industrial applications, including chemical reactors, heat exchangers, and environmental processes. This study aims to investigate the intricate relationship among chemical reactions, Maxwell fluid flow, and energy generation over a porous vertical plate. By analyzing the flow patterns, temperature distribution, and heat transfer characteristics, we can gain insights into how these factors interact and influence each other. This research contributes to our fundamental understanding of fluid dynamics and heat transfer and has practical implications for optimizing processes and systems where these phenomena are encountered.

Kashif et al. [1] undertook a comprehensive analysis focusing on the dynamic interactions of pressure gradients and heat transfer within the context of a Maxwell fluid. Their investigation centred on the boundary layer phenomena occurring over a stretched sheet. Liu et al. [2] delved into the intricacies of energy transfer in a Maxwell fluid flow situated on a moving plate. Their research featured applying a time-fractional constitutive model to unravel this context’s complex heat transfer dynamics. Imran et al. [3] embarked on a study examining the magnetohydrodynamic (MHD) effects within a Maxwell fluid boundary layer flow across an exponentially infinite perpendicular surface. Their work introduced novel insights by incorporating slip conditions and Newtonian heating mechanisms at the boundary. Sudarmozhi et al. [4] directed their research efforts towards investigating radiation and heat transfer phenomena embedded within the dynamics of a Maxwell fluid flowing over an inclined permeable plate. This work shed light on the interplay of these complex physical processes. Na et al. [5] took a closer look at a Maxwell fluid flow confined between vertical plates, considering the influence of thermal flux and the presence of damped shear forces. Their study provided a comprehensive perspective on the behaviour of Maxwell fluids in this specific geometry. Riaz and Iftikha [6] thoroughly examined heat transfer flows involving magnetohydrodynamic (MHD) Maxwell fluids. They explored the comparative aspects of local and nonlocal differential operators, offering valuable insights into the heat transfer behaviour of these complex fluids. Rasa and Ullah [7] delved into the fascinating domain of MHD Maxwell fluid heat transfer. Their research featured a comparative analysis utilizing Caputo and Caputo–Fabrizio derivatives, contributing to a deeper understanding of the fluid’s thermal dynamics. Riaz et al. [8] focused on the intriguing power-law kernel characteristics inherent in magnetohydrodynamic (MHD) Maxwell fluids. Their study included the investigation of transport phenomena while considering ramping boundary conditions and specific function-based solutions. Hanif et al. [9] employed computational methods to ascertain the fractional behaviour of Maxwell fluids within boundary layers, particularly by examining heat transfer flows. Their research introduced innovative techniques for characterizing these complex fluids at the boundary. Khan and Rasheed [10] explored fluid flow dynamics featuring fractional Maxwell fluids passing through a Forchheimer medium. Their study provided insights into this specific flow scenario’s unique challenges and behaviour. Shah et al. [11] ventured into the realm of free convection flows, specifically considering Prabhakar-like fractional Maxwell fluids. Their research extended to encompass the influence of generalized heat transmission, enriching our understanding of these complex fluid dynamics. Zhang et al. [12] investigated the behaviour of Maxwell fluids within natural convection flows, focusing on generalized heat transfer and the application of Newtonian heating mechanisms. Their work contributed to comprehending heat transfer in these intricate fluid systems. Al Hajri et al. [13] researched into the Maxwell–Cattaneo properties within convective energy transfer flows inside a square enclosure filled with a permeable medium. Their research offered valuable insights into the thermal dynamics of such complex systems.

The complexities of magnetohydrodynamic-free convective flow in porous media (see for instance [14] and references therein) have piqued the interest of scientists across various domains, spanning from applications in flow meters, pumps, generators, and plasma jet engines to accelerators, metallurgical processes, chemical industry operations, material processing, industrial power engineering, and nuclear engineering. Mass transfer, which entails the movement of substances within a system due to differences in concentration, plays a pivotal role in many biological processes, such as respiration, blood oxygenation, renal function, osmosis, and absorption of nutrients and medications. Researchers have undertaken various investigations (see for instance [15, 16] and references therein) to unravel the intricacies of such dynamic systems. Reddy et al. [17] explored the behaviour of Casson fluid flowing across an oscillating permeable plate, systematically accounting for the influences of rotation and diffusion thermos on magnetohydrodynamic heat-generating mixed convection flow. Dadheech et al. [18] delved into the unique attributes of the Williamson fluid flow across a perpendicular plate, explicitly considering the impacts of Cattaneo–Christov heat flux.

Notably, the interplay between double-diffusive processes and chemical reactions adds complexity to these studies. The rates of double-diffusive phenomena can be substantially altered due to ongoing chemical reactions, which, in turn, depend on the reaction type—whether it is a mixed or identical reaction. The reaction order, characterized as first order when the reaction rate is directly proportional to the concentration, assumes paramount significance in understanding these dynamic systems. The natural environment rarely contains pure substances such as pristine water or unadulterated air. Consequently, comprehending the intricacies of chemical reactions is pivotal for improving numerous chemical technologies, including the synthesis of polymers, food processing, and fabricating ceramic or glass products. Technicians and scientists across various fields engage in processes where physical and chemical transformations convert raw materials into finished goods, underpinning key industries such as ammonia, sulfuric acid from sulphide ores, and a wide array of petrochemicals derived from petroleum feedstocks. The confluence of chemical reactions with heat and mass transfer has been the subject of rigorous investigation by researchers. Hussain et al. [19] explored natural convective heat transfer with the concurrent impact of hall current and chemical reactions, examining the impact of heat absorption past an accelerated moving plate in a rotating system. Saidulu and Reddy [20] delved into the intricacies of energy and mass transfer in hydromagnetic micropolar flow over a sheet, considering viscous dissipation along with chemical reactions. Sedia et al. [21] meticulously studied the convective heat-mass transport of a generalized Maxwell fluid under the impact of exponential heating and chemical reactions. Iftikhar et al. [22] investigated energy and mass transfer within magnetohydrodynamic Maxwell fluid dynamics across an infinite perpendicular plate.

Furthermore, Hanif [23] applied advanced numerical techniques such as the Cattaneo–Friedrich model and the crank-Nicolson method to analyze the behaviour of upper-convected Maxwell fluid [24], passing over a perpendicular plate. Abbas et al. [25] observed the computational implications of the Darcy–Forchheimer [26] relation, reduced gravity, and the influence of an externally applied magnetic field on the radiative fluid flow and heat transfer on a sphere. In a separate study, Abbas et al. [27] analyzed MHD Williamson nanofluid flow and heat transfer, investigating its interaction with a nonlinear stretching sheet implanted in a permeable medium. Munir et al. [28] focused on a numerical study involving entropy generation in the stagnation point flow of Oldroyd-B nanofluid. Cui et al. [29] also delved into a numerical investigation of the thermal transport properties of the three-dimensional bioconvective nanofluid flow.

1.1. Novelty of the Subject, Motivation, and Its Application

Maxwell fluids refer to a class of non-Newtonian fluids that exhibit viscoelastic behaviour. When studying the flow of Maxwell fluids on a perpendicular plate in the presence of a chemical reaction, several unique aspects and applications arise, contributing to the novelty of this subject. Maxwell fluids display viscoelastic properties, exhibiting viscous and elastic responses to applied stresses. This behaviour differs significantly from Newtonian fluids, which exhibit viscosity without elasticity. Understanding the flow characteristics of such fluids on a perpendicular plate provides insights into how these materials deform and flow under various conditions. The flow of fluids over a perpendicular plate is a classical problem in fluid mechanics. Investigating the Maxwell fluid flow over a vertical plate involves examining how the fluid’s viscoelastic properties influence boundary layer formation, velocity profiles, and heat/mass transfer characteristics. The behaviour differs from Newtonian fluids due to the attendance of elastic effects. Introducing chemical reactions into the Maxwell fluid flow scenario adds another layer of complexity. The interaction between the fluid flow and chemical reactions can significantly alter flow properties, such as viscosity and density gradients. These alterations can, in turn, affect the overall flow pattern and boundary layer development. Studying the Maxwell fluid flow with a chemical reaction on a perpendicular plate is crucial in understanding heat and mass transfer processes. The chemical reaction can influence heat transfer and mass between the fluid and the plate, affecting the overall system’s behaviour and efficiency. Understanding the Maxwell fluid flow on a perpendicular plate with chemical reactions finds applications in various fields. These include chemical engineering (for designing reactors or separation processes), material sciences (in understanding and processing viscoelastic materials), and environmental engineering (in modelling pollutant transport and reactions in fluids). Researchers use mathematical models and simulations to predict and analyze the behaviour of the Maxwell fluid flow with chemical reactions. These models involve complex equations representing fluid dynamics, heat transfer, mass transfer, and chemical kinetics, providing insights that help optimise processes and design efficient systems. The motivation behind studying the Maxwell fluid flow on a vertical plate with chemical reactions lies in its potential to advance multiple fields, improve industrial processes, deepen scientific understanding, and contribute to developing innovative materials and technologies.

1.2. Novelty of the Proposed Problem

These collective research endeavours significantly deepen our comprehension of the intricate interplay among fluid dynamics, energy transfer, and mass transfer, along with chemical reactions across diverse physical and technological environments. While prior studies have offered valuable insights into different facets of fluid dynamics, they have not thoroughly addressed the intricate interconnection of viscoelastic fluid flow, magnetohydrodynamics (MHD) over a perpendicular plate, combined radiation effects, heat generation, and chemical interactions influencing energy and mass transfer through a porous medium. No authors have explored the Maxwell fluid flow on a perpendicular plate while considering this combination of parameters. Recognizing this research gap, our study’s primary aim is to fill this void in knowledge and conduct a comprehensive investigation into this multifaceted scenario. To confront this intricate problem, we transform the governing partial differential equations into more manageable ordinary differential equations (ODEs). These nonlinear ODEs are then solved via the robust MATLAB inbuilt software tool, the bvp4c solver. Our analysis thoroughly explores various critical parameters, encompassing temperature, velocity, and concentration contours and evaluating the local Sherwood number, skin friction coefficient, and Nusselt number. Our research seeks to illuminate the intricate dynamics in this complex fluid flow scenario through a systematic study and analysis of these governing parameters. This investigation not only enriches our fundamental understanding of fluid dynamics but also carries the potential to advance various practical applications in engineering, environmental science, and industrial processes. This research emphasizes and aims to address the following questions:(i)What impact does fluid flow have on the behaviour of skin friction, Nusselt number, and Sherwood number concerning the chemical reaction parameter?(ii)How does fluid flow influence the velocity profile, temperature profile, and concentration contour concerning the chemical reaction parameter?(iii)How does fluid flow influence the velocity, temperature, and concentration outline concerning the heat generation parameter?

2. Mathematical Formulation

Consideration is given to a non-Newtonian viscoelastic fluid that is laminar, steady, and two-dimensional. Consideration is given to porous media exposed to magnetic fields. It is presumed that the y-axis is over the plate and the x-axis is perpendicular to it. Let u and stand for the horizontal velocity in the x direction and the vertical velocity in the y direction, respectively. C represents the fluid concentration, and T represents the temperature. The convective boundary conditions obstruct the local temperature at the wall. A magnetic field with a constant strength of B0 is provided perpendicular to the plate and along the y-axis of the vertical porous plate. The continuous concentration and temperature of the vertical plate are also supposed to be and , whereas those of the surrounding fluid are T and C. The chemical reaction has been added to the flow, so the Boussinesq approximation says that all thermophysical properties are probably in the variable of the linear impetus expression. The physical mode of the problem is offered in Figure 1.

Upper-convected Maxwell fluid convinces are as follows:where is the stress tensor, is the zero shear rate viscosity, represents the symmetric part of the velocity gradient tensor, and is fluid’s relaxation time.

Upper-convected time derivative is as follows:where is the velocity tensor and is the material time derivative

“The Cartesian form of the governing equations, along with the boundary layer calculations as mentioned by Reference [30], is presented as follows.”

The equation of continuity of the governing equation is as follows:

In manufacturing and technological applications, non-Newtonian fluids are preferable over Newtonian fluids. The governing equations are complicated, and getting solutions to the following equations is more challenging. Numerous non-Newtonian fluid models have been proposed and are currently being researched due to the complexity of fluids. So, with the occurrence of Maxwell fluid, the equation of conservation of linear momentum is as follows:

The equation of conservation of energy in the occurrence of heat generation (Q) and thermal radiation is as follows:

Here, qr characterizes the radiative heat flux and Cp characterizes the specific heat; if the thermal conductivity of the fluid k and temperature of the fluid have a linear relation, then according to Roseland [31], qr is calculated as follows.

; here, it can become close using series development about the ambient temperature of the fluid represented as T by .

If we disregard higher-order terms, let us inscribe from the above equations .

The equation of conservation of mass with chemical reaction occurrence is as follows:

The given problem’s boundary conditions for the dimensional governing equations are as follows:

Using the similarity transformation, the nonlinear PDEs are changed into ODEs.

Dimensionless numbers lower the variables used to characterize them, minimizing the experimental data required to establish correlations between physical phenomena and scalable systems.

Flow characteristics typically expressed with the occurrence of stream function are given here.

PDEs are converted into ODEs using the above quantities.

The dimensionless momentum equation in the presence of dimensionless porous and MHD effect is as follows:

The dimensionless energy equation in the occurrence of radiation and heat generation is as follows:

In the chemical and water industries, double diffusion with higher-order chemical processes utilizing a porous medium has several applications. In addition to chemical reactions, industrial processes such as filtration, cooling, distillation, and drying rely on transmitting heat and mass. The relationship between the concentration and the rate of a chemical reaction is explained by the order of the reaction.

Nondimensional numbers help determine the flow properties of a fluid.

The porosity is typically given by the ratio of void space volume to the total volume of the material.

The Maxwell fluid model involves relaxation time (τ) and viscosity (μ). The nondimensional form of the Maxwell fluid parameter can be expressed as a ratio of relaxation time to characteristic time or any relevant length scale .

In the context of heat generation in a system, the nondimensional representation depends on the specific design and the governing equations involved. However, we can typically define nondimensional heat generation as the ratio of heat generated to some characteristic value or scale in the system.

The nondimensional representation of chemical reactions can vary based on the specific reaction rate and other factors.

The appropriate corresponding dimensionless boundary conditions of the current problem are as follows:

2.1. Essential Quantities of Engineering Interest
2.1.1. Coefficient of Skin Friction

Coefficient of skin friction is as follows:

2.1.2. Local Nusselt Number

Local Nusselt number is as follows:

2.1.3. Local Sherwood Number

Local Sherwood number is as follows:

Here, is the mass flux.

By using the dimensionless quantities listed above,

3. Methodology and Results of the Problem

To identify the answers to the dimensionless equations with boundary conditions, the shooting technique known as BVP4c in MATLAB software is utilized in order to obtain the numerical solution which is shown in Figures 210. Results can be summarized as graphs for the nondimensional ordinary differential equations. Computational results are displayed for each parameter through temperature, velocity, and concentration contours. We investigate the skin friction coefficient, the Nusselt number, and the Sherwood number about fluid flow parameters (see Tables 13). This work aims to determine how a chemical reaction affects the steady state with energy and mass transport of an electrically conducting fluid that can either generate or absorb heat. “If not specified differently, the study assumes default or unchanged parameters.” M= 1.0, Pr = 6.2, Rd = 0.5, Nr = 1.0, Ra = 0.5, λ= 1.0, β= 0.2, Le = 5, Q= 0.01, and Kr = 0.5 are taken by the standard literature. These parameter values are taken from the classic literature. Thus, unless otherwise noted on the relevant graph, all graphs correspond to these values.

3.1. Comparison of the Results to an Earlier Publication

To confirm the accuracy of our solution derived using the bvp4c MATLAB solver, we compared our recent results to those from a previous study conducted by Nield [32]. The significant alignment between our findings and those of the earlier research inspires us to examine the impacts of different controlling parameters on the behaviour of Maxwell fluid flowing over a porous perpendicular plate. The validation of our code against published research outcomes is presented in Table 4.

3.2. Velocity Profile

Figure 2 illustrates how the velocity profile is affected by the Maxwell fluid flow. As the value of the Maxwell fluid parameter (Deborah number) surges, so does the velocity profile. The increase in the Maxwell fluid parameter, also known as the Deborah number, leads to a corresponding increase in the velocity contour on a vertical plate in the Maxwell fluid flow. This is due to the physical mechanism of the fluid’s viscoelastic properties. Maxwell fluids exhibit both viscous and elastic behaviour. As the Deborah number increases, it signifies a higher dominance of flexible effects over viscous effects in the fluid flow. This higher elasticity allows the fluid to deform and respond more readily to applied forces, resulting in increased velocity profiles on the vertical plate. A higher Deborah number indicates that the liquid can more effectively transmit applied forces, causing more significant deformation and, consequently, higher velocity profiles along the vertical plate in the Maxwell fluid flow.

Figure 3 demonstrates the effect on the velocity outline when augmenting the chemical reaction (Kr) on a vertical plate in Maxwell fluid flow, resulting in an increased velocity outline. Elevating the Kr on the sheer scale during Maxwell fluid flow leads to an observed enhancement in the velocity profile. The physical mechanism behind the increase in the velocity profile due to an elevated chemical reaction parameter in Maxwell fluid flow on a vertical plate involves several factors related to the chemical reactions occurring within the fluid. Chemical reactions often involve the release or absorption of energy. This change in power can affect the temperature and, consequently, the fluid’s viscosity. In some cases, an increase in Kr can result in a decrease in density, facilitating faster flow and, hence, an increased velocity profile. Chemical reactions may involve species transport or generation, altering the mass transport properties of the fluid. This could lead to changes in density or composition, impacting the fluid’s flow behaviour and influencing the velocity profile. Chemical reactions can modify the rheological properties of the liquid, such as viscosity and elasticity. Alterations in these properties affect how the fluid responds to applied forces, potentially leading to changes in velocity profiles. Reactions that generate or consume heat can influence the temperature distribution within the fluid. Temperature changes can, in turn, affect the fluid’s viscosity and flow characteristics, controlling the velocity outline. Overall, the interaction between chemical reactions and fluid properties such as viscosity, density, temperature, and rheology directly impacts the flow behaviour and subsequently influences the velocity profile on the vertical plate in Maxwell fluid flow.

Figure 4 provides insight into how heat generation impacts the velocity profile, showcasing a direct correlation between increased heat generation and an augmented velocity profile. The physical mechanism underlying this phenomenon involves several interconnected processes within the fluid. As heat is generated, the fluid’s temperature rises, leading to thermal expansion. This expansion decreases the fluid density near the heated surface, creating a less dense layer that tends to increase due to buoyancy. The warmed fluid, being less viscous, ascends upward due to buoyancy forces. This upward movement of the heated fluid establishes a vertical flow component, contributing to an overall increase in the velocity profile. Heat generation induces a heightened temperature gradient in the liquid. This amplified gradient encourages convection, characterized by the bulk movement of the fluid. As the buoyant-heated fluid rises, it draws in denser fluid from below, fostering a convective flow that supplements the overall fluid velocity. Higher temperatures often lead to a reduction in fluid viscosity. With lower density, the fluid encounters less resistance to flow, enabling it to move more swiftly. This diminished resistance further bolsters the increase in fluid velocity.

3.3. Temperature Outline

Figure 5 visually elucidates the relationship between heightened chemical reaction values and their influence on the temperature profile in our study, incorporating the underlying physical mechanisms at play within the fluid. Chemical reactions typically produce heat as a byproduct. Increasing the chemical reaction values results in more significant heat generation within the system. The additional heat energy generated raises the overall temperature of the fluid. This thermal energy transfer heats the surrounding liquid near the reaction site. The heated fluid conducts heat to adjacent areas, fostering a widespread increase in temperature across the fluid. This conduction mechanism influences the entire temperature profile. Higher chemical reaction values intensify fluid motion and mixing. This increased agitation aids in spreading elevated temperatures throughout a broader region within the fluid. Chemical reactions create temperature gradients in the liquid. These gradients induce fluid motion through natural convection, where warmer fluid ascends while more excellent fluid descends. This circulation alters the temperature contour within the system.

Figure 6 depicts an escalation in the heat generation parameter directly correlating with a heightened temperature profile. This relationship is primarily governed by principles rooted in heat transfer and thermodynamics. The heat generation parameter signifies a heat source within the system, typically represented by various physical processes such as chemical reactions, electrical resistance, or mechanical work. Its increase results in a proportionate rise in the rate of energy generation within the system. The temperature profile denotes a material or design’s spatial and temporal temperature variation. It represents temperature distributions across different points within the system. Augmenting the heat generation parameter amplifies the rate at which heat is generated within the system. This increase causes a localized temperature rise near the heat source as more heat accumulates.

Consequently, the temperature gradient steepens within the system, indicating a more significant temperature disparity between regions close to the heat source and those farther away. Continued heat generation perpetuates its impact across the entire system, leading to an overall elevation in the temperature profile. Consequently, both the local temperature near the heat source and the temperature in the whole system increase due to the escalated heat generation.

3.4. Concentration Profile

Figure 7 illustrates the concentration outline as the value of the chemical reaction increases, revealing a declining concentration trend. The decrease in the concentration observed in the concentration outline with an increase in the value of the Kr, as depicted in Figure 7, can be attributed to several physical mechanisms associated with chemical reactions within the system. Elevating the chemical reaction parameter typically implies a higher rate of chemical reactions within the system. As reactions progress more rapidly, reactants are consumed faster, reducing their concentration over time. Chemical reactions often involve the consumption of reactants to produce products. As the reaction rate increases due to higher parameter values, more reactants get converted into products, consequently depleting the concentration of reactants in the system. Specific reactions might be reversible and tend towards reaching an equilibrium state. Increasing the chemical reaction parameter could shift the equilibrium towards products, thereby reducing the concentration of reactants and favouring product formation. Higher chemical reaction rates may intensify fluid motion and mixing within the system. This increased mixing can disperse the reactants more uniformly, potentially leading to lower localized concentrations as the reactants become more distributed. Accelerated chemical reactions induce faster diffusion rates of reactants, causing them to spread more rapidly within the system. This enhanced diffusion could contribute to a decrease in the concentration profile depicted in Figure 7.

3.5. Engineering Quantities

Figures 810 portray the profiles of skin friction, Nusselt number, and Sherwood number, respectively, concerning the escalating value of the chemical reaction parameter. Notably, when the chemical reaction parameter rises, the skin friction and Sherwood number profiles demonstrate an upward trend, whereas the Nusselt number profile exhibits a declining trend in isolation. The varying behaviour observed in the skin friction, Nusselt number, and Sherwood number profiles concerning the rise in the chemical reaction parameter can be attributed to distinct physical mechanisms associated with heat and mass transfer, as well as fluid dynamics.

3.5.1. Skin Friction

Physical mechanism: increasing the chemical reaction parameter often leads to elevated fluid motion and mixing due to higher reaction rates. This augmented fluid movement tends to increase the frictional forces exerted on the surface, increasing the skin friction profile. The heightened chemical reactions can intensify the momentum transfer between the fluid and the surface, leading to an enhanced skin friction profile.

3.5.2. Nusselt Number

Physical mechanism: the decrease in the Nusselt number profile can be attributed to the complex interplay between heat transfer mechanisms. Increasing chemical reactions might alter the flow patterns and energy transfer characteristics to suppress the convective heat transfer rate while potentially enhancing the conductive energy transfer. As a result, the convective heat transfer rate decreases despite the increased fluid motion, leading to a decline in the Nusselt number profile.

3.5.3. Sherwood Number

Physical mechanism: the surge in the Sherwood number profile is typically associated with mass transfer enhancement due to intensified fluid motion caused by higher chemical reaction rates. This increased fluid movement and mixing facilitate more excellent mass transport or diffusion of substances from the surface into the fluid, leading to an elevated Sherwood number profile.

4. Conclusion

This study fills a gap by numerically analyzing Maxwell fluid behaviour on a porous perpendicular plate, considering chemical reactions, heat generation, and magnetohydrodynamic effects. Using MATLAB’s bvp4c solver and transformation techniques, we tackle the challenges posed by nonlinear equations, exploring their impact on temperature, velocity, and concentration. In addition, we evaluate Sherwood numbers, Nusselt numbers, and skin friction coefficients to assess chemical reaction effects. Our detailed investigation sheds light on the intricate relationship among Maxwell fluid flow, chemical reactions, and heat generation on a vertical plate. These insights are crucial for applications in chemical processes, heat exchangers, and environmental engineering, offering opportunities for process optimization and further research in similar fluid flow scenarios. Key findings from our study are outlined as follows:(1)We observed that changes in the chemical reaction parameters and heat generation significantly affect the velocity profile. Higher chemical reaction rates and heat generation led to a surge in fluid velocity near the plate.(2)The temperature profile exhibited variations corresponding to the levels of chemical reactions and heat generation. Elevated chemical reaction rates and increased heat generation caused higher temperatures in the fluid near the plate.(3)Our study also considered the concentration profile, which reflects how chemical reactions influence fluid substance distribution. The concentration profiles demonstrated changes in response to varying chemical reaction rates and heat generation levels.(4)We analyzed the skin friction outline, indicating the fluid flow resistance along the plate. Chemical reactions and heat generation influenced skin friction, with higher values of these parameters leading to altered frictional forces.(5)The Nusselt number profile, a dimensionless energy transfer coefficient, showed significant variations with changes in chemical reactions and heat generation. This parameter is crucial for understanding heat transfer rates.(6)We also considered the Sherwood number profile related to mass transfer rates. Chemical reactions and heat generation were observed to impact the mass transfer rate, resulting in variations in the Sherwood number profile.

4.1. Future Work

There is a potential to expand the scope of this paper by introducing additional parameters or altering the geometry, alongside incorporating a more extensive range of fluid variables.

Nomenclature

C:Concentration of the fluid
:Temperature on the wall (K)
C:Concentration on the free stream
M:Magnetic field
β:Deborah number
T:Temperature of the fluid (K)
Pr:Prandtl number
, :Velocity of the fluid (ms−1)
B0:Magnetic component
:Dimensional chemical reaction
Cf:Skin friction
θ:Dimensionless temperature
DB:Mass diffusion
f:Dimensionless stream function
D:Diffusion term
K:Permeability of a porous medium
T:Temperature on the free stream
Kr:Dimensionless chemical reaction parameter
λ1:Maxwell fluid parameter (s)
Rd:Dimensionless radiation
Cp:Specific heat at constant pressure (JK−1·kg−1)
υ:Kinematic viscosity (m2s−1)
βc:Coefficient of mass expansion (K−1)
:Dimensional heat generation
ρ:Density of the fluid (kg−3)
Q:Dimensionless heat generation
Rex:Local Reynolds number
Ra:Rayleigh number
:Stefan–Boltzmann constant
λ:Porosity dimensionless parameter
βT:Coefficient of thermal expansion (K−1)
ϕ:Dimensionless concentration
η, ψ:Similarity variables (m2s−1)
μ:Dynamic viscosities (kg·m−1·s−1)
k:Thermal conductivity (Wm−1k−1)
:Gravitation (ms−2)
ρcp:Heat capacity (Jkg−1·K−1).

Data Availability

No data were required to perform this research.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

The authors extend their appreciation to the Deanship of Scientific Research at King Khalid University for funding this work through large group Research Project under grant no. RGP2/340/44.