Signiﬁcance of Thermophoretic Particle Deposition, Arrhenius Activation Energy and Chemical Reaction on the Dynamics of Wall Jet Nanoﬂuid Flow Subject to Lorentz Forces

: The need for effective heating and cooling systems in the automotive, chemical, and aerospace industries is driving a rapid proliferation of heat-transfer technology. In recent times, GO (Graphene Oxide) has been emerging as one of the most promising nanoparticles because of its uninterrupted behavior of electrical conductivity even at a minimum carrier concentration. Due to this incentive, the behavior of jet ﬂow with heat and mass transfer features of electrically conducting based kerosene oil (KO) ﬂuid dispensed by graphene nanoparticles was studied. In addition, the activation energy, irregular heat source/sink, thermophoretic particle deposition, and chemical reaction are also provoked. In order to provide numerical results, the boundary value problem of fourth-order (bvp4c) solver was used. The graphs were used to illustrate the effects of relevant parameters on the ﬂuid ﬂow, heat, and mass transfer rates. The incorporation of graphene nanoparticles signiﬁcantly improves heat conductivity. Additionally, the nanoparticle volume fraction augments the temperature and concentration proﬁle while the velocity proﬁle declines. Moreover, the temperature enhances due to the heat source, whilst the contrary behavior is observed in the presence of the heat sink. Furthermore, the shear stress increases up to 12.3%, the Nusselt number increases up to 0.119%, and the Sherwood number increases up to 0.006% due to the presence of nanoﬂuid. Finally, we can conclude that the latest work will be useful for thermal cooling systems, including cooling for engines and generators, nuclear systems, aviation refrigeration systems, and other systems.


Introduction
The past century has witnessed a revolution in the industries by employing extraordinary applications of nanofluids including electrical cooling, power generation, solar energy, building technology, and plant processing. Nanofluid offers novel alternate substances to utilize in several fields [1,2], one of them is the improvement of lubricants [3,4]. Graphene Oxide (GO) is among the most promising nanoparticles given its constant electrical conductivity even at low carrier concentrations. GO has an extremely high electron mobility when compared to other available particles. External magnetic fields have the power to control the heat-transfer and flow properties as well as the thermal and physical properties of magnetic nanofluids. This is highly beneficial for managing processes in metallurgy. A simple method used to create a nanocomposite of Co 3 O 4 nanoparticles embedded in conducting graphene as an improved anode material for battery cells was disclosed by Wu et al. [5]. They discovered that the Co 3 O 4 /graphene nanocomposite exhibits superior rechargeable performance with a significant amount of reversible capacity. Graphene-based particles have garnered a great deal of interest for usage in energy storage. Pumera [6] outlined the experimental as well as a theoretical exploration of using graphene-based applications such as storage systems of hydrogen, and lithium batteries. Ahmad et al. [7] present one of the significant uses of graphene in relation to electrical conductivity. Khan et al. [8] explored the magnetic impact on the nanofluid flow through the stretchable cylinder and found that the magnetic and thin film parameters decelerate the velocity and accelerate the temperature distribution. The non-Newtonian (Casson and Williamson) non-steady flow and thermal properties of electrically conducting water-based thin fluid film dispensed with graphene nanoparticles via a stretchy surface were considered by Ahmad et al. [9]. They observed that the low resistance for heat conduction is provided by the greater surface area of the film thickness due to graphene nanoparticles. Recently, Khan et al. [10] investigated the buoyancy flow of water-based graphene nanoparticles through a vertical surface with radiation impact and performed a stability analysis to find a stable solution.
Thermophoresis particle deposition in a fluid flow is important in several technical processes, including the protection of nuclear reactors, air cleaners, heat exchangers, and burners for powdered coal. The thermophoresis phenomenon is caused by several categories of particles acting differently when exposed to a temperature difference. Tiny particles dispersed in an NIG (non-isothermal gas). will reach a velocity during thermophoresis. This process greatly increases the deposition velocity of tiny particles in the direction of cooling temperature, while it has no effect on large particles. Thermophoresis allows small particles to mend on a cold surface. Alam et al. [11] inspected 2D (two-dimensional) steady magneto flow across an inclined semi-infinite plate including thermophoresis and erratic suction. They found that the suction parameter uplifts the thermal and concentration boundary layer thicknesses. The mixed convection flow through a vertical sheet with a magnetic effect and thermophoresis effect was examined by Damseh et al. [12]. They observed that the velocity of fluid decreases due to magnetic numbers. Rahman [13] demonstrated the thermophoretic particle deposition and enforced magnetic influence on the flow of nanofluid through the rotating system. It has been found that the rate at which thermophoretic particles deposit is greatly influenced by thermal diffusion, slip mechanism, magnetic field, diffusion-thermo, and radiation. Recently, numerous authors [14][15][16][17][18] have also reported the impact of thermophoretic particle deposition under diverse circumstances.
Numerous chemicals, mechanical, and manufacturing processes depend heavily on chemical changes. A large amount of oxygen known as activation energy is needed to control a process. The energy needed to initiate a chemical reaction is often referred to as activation energy. The fundamental idea behind this topic was created in 1889 by Svante Arrhenius. Shafique et al. [19] observed that the chemical reaction in connection to activation energy has a significant impact on the fluid-flow behavior. They discovered that the solute content in a binary mixture correlates with the activation energy and the rotation parameter. Khan et al. [20] explored the impact of activation energy on the radiative 3D (three-dimensional) flow induced by a cross nanofluid with binary chemical reaction by employing the Buongiorno model. It can be seen from the generated graphical data that Brownian motion and thermophoresis have a major impact on the heat-transfer process. The effects of Joule heating and activation energy with binary chemical change incorporating the MHD (magneto-hydrodynamics) Jeffrey fluid towards a contoured compressed surface were reviewed by Khan and Alzahrani [21]. They revealed that the curvature parameter enhances the entropy generation. Najafabadi et al. [22] studied the effects of radiation, conduction, and convection heat transfer in a movable fin with even velocity and changeable thermal conductivity. Rana et al. [23] utilized the non-Newtonian model to investigate the blood flow incorporated nanofluid through bio-convection with activation energy and chemical reaction. Recently, Reza-E-Rabbi et al. [24] investigated the radiative magnetism flow of a Casson nanofluid past a stretchy periphery with activation energy and chemical reaction. They observed that the temperature is elevated due to Brownian motion, thermophoretic, and radiation effects.
The migration of gaseous particles from hot to cold locations is known as thermophoresis. This process informs the system that micron-sized particles are prohibited in nonisothermal gas and move more quickly along the path of reducing thermal gradients. The significance of theoretical and experimental understanding of thermophoretic particle deposition cannot be understated due to its wide range of thermal engineering and industrial applications. The distribution of temperature and the deposition of particles may change as a result of thermogenesis' possible effects. These include heat exchangers, air cleaners, filtration, powdered coal burners, environmental and atmospheric pollution, ventilation systems in buildings, and safety for nuclear reactors. Previous explorations on the thermophoresis problem were presented by Goren [25], Talbot et al. [26], Mills et al. [27], and Batchelor and Shen [28]. Tsai [29] provided a straightforward method to examine the thermophoresis impact on the flow past a flat surface, which was inspired by earlier investigations. In this study, the thermophoresis is calculated by using a first-order chemical reaction and suction-like convection term. The hydromagnetic flow was then studied by Chamkha and Issa [30], and the porous medium impact was considered by Chamkha and Pop [31]. Alam et al. [32] additionally investigated the flow on an adiabatic porous sheet. Their studies revealed that the Nusselt number for the injection of fluid is sensitive due to the increasing inclination of the angle. Moreover, the changeable properties of the fluid were presented by Das et al. [33]. Notably, the fluid flow through a shrinkable sheet incorporating thermophoresis was observed by Sinha and Misra [34] and Zaib and Shafie [35]. In recent times, Jyothi et al. [36] explored the flow flanked by parallel plates with thermophoresis whereas Shah et al. [37] examined the flow of a second-grade fluid across a flat plate induced by erratic fluid properties.
The temperature of the entering jet should vary from that of the static fluid in industrial applications of jets such as melting plastic and metal sheets, drying fabrics, chilling microelectronic equipment, etc. As a result, the thermal boundary equation needs to be solved. It appears that Tetervin [38] presented the first solution for the conventional laminar wall jet. A famous non-trivial exact solution with exponential decline at infinity (e-jet) was obtained by Glauert [39]. Riley [40] investigated the thermal features of the preheated radial jet. Schwarz and Caswell [41] analyzed the characteristics of heat transfer in connecting to the Glauert wall jet flow, in which they presented the closed form solution including the impact of steady heat-flux and temperature. An approximate solution to the time-dependent energy equation was provided by Gorla [42] in the context of Glauert's wall jet flow. It has been shown that the features of heat transfer can be reduced when a variation in surface temperature is imposed. Aly and Pop [43] inspected the hybrid nanofluid induced by wall jet flow with thermal energy and found dual solutions. It was discovered that only suction makes the lower-branch and the upper-branch solutions conceivable.
The literature indicates that no effort has been made to investigate jet flow induced by oil-based graphene nanofluid including thermophoretic particle deposition, activation energy with chemical reaction, and erratic heat sink/source. Therefore, the novelty of the present work is to inspect the jet flow of GO nanofluid by including the thermophoretic particle deposition, activation energy with chemical reaction, and erratic heat source/sink which will close this research gap. The leading flow equations are turned into a set of ODEs (ordinary differential equations). The well-known bvp4c is used to solve the resulting equations. The objective is to improve the regular fluid's thermal conductivity for a better heat-transport phenomenon through portable density and heat exchangers employed in contemporary thermal procedures. The effects of all relevant parameters on the velocity, concentration, and temperature were graphically displayed and explained.

Mathematical Background of the Problem
Consider the phenomenon of the Glauert model which is schematically highlighted in the form of Figure 1. In this model, a 2D MHD laminar wall jet heat and mass transfer flow of nanofluid blown from a thin slit on the upper section of a static permeable surface with influences of an irregular heat sink/source was explored. Further, the analysis of the given problem also comprised the significant effect of thermophoretic particle deposition and Arrhenius activation energy with chemical reaction. It is assumed that the coordinates of the xand yaxes are measured along the wall and orthogonal to it, respectively. The uniform wall temperature and ambient temperature are denoted by T w and T ∞ , respectively. Further, assume that C = C w = 0 and C ∞ represent the respective uniform wall and ambient concentration of the fluid. However, the irregular heat sink/source term is taken in the posited energy equation which is symbolically denoted by Q and later it is explained in detail. Additionally, the magnetic field [44,45] is taken as B = B 0 x −3/4 to obtain the similarity equations. Meanwhile, B 0 denotes the magnetic field strength. These suppositions allow the leading equations to be formulated as [42,43]: subject to the BCs (boundary conditions) are: Lubricants 2022, 10, x FOR PEER REVIEW 4 of 22 The literature indicates that no effort has been made to investigate jet flow induced by oil-based graphene nanofluid including thermophoretic particle deposition, activation energy with chemical reaction, and erratic heat sink/source. Therefore, the novelty of the present work is to inspect the jet flow of GO nanofluid by including the thermophoretic particle deposition, activation energy with chemical reaction, and erratic heat source/sink which will close this research gap. The leading flow equations are turned into a set of ODEs (ordinary differential equations). The well-known bvp4c is used to solve the resulting equations. The objective is to improve the regular fluid's thermal conductivity for a better heat-transport phenomenon through portable density and heat exchangers employed in contemporary thermal procedures. The effects of all relevant parameters on the velocity, concentration, and temperature were graphically displayed and explained.

Mathematical Background of the Problem
Consider the phenomenon of the Glauert model which is schematically highlighted in the form of Figure 1. In this model, a 2D MHD laminar wall jet heat and mass transfer flow of nanofluid blown from a thin slit on the upper section of a static permeable surface with influences of an irregular heat sink/source was explored. Further, the analysis of the given problem also comprised the significant effect of thermophoretic particle deposition and Arrhenius activation energy with chemical reaction. It is assumed that the coordinates of the x − and y − axes are measured along the wall and orthogonal to it, respectively.
The uniform wall temperature and ambient temperature are denoted by w T and T ∞ , respectively. Further, assume that 0 w C C = = and C ∞ represent the respective uniform wall and ambient concentration of the fluid. However, the irregular heat sink/source term is taken in the posited energy equation which is symbolically denoted by ''' Q and later it is explained in detail. Additionally, the magnetic field [44,45] is taken as = to obtain the similarity equations. Meanwhile, 0 B denotes the magnetic field strength.
These suppositions allow the leading equations to be formulated as [42,43]:  Here, u and v are the components of velocity along the respective xand yaxes while T and C are the temperature and concentration of the nanofluid. It is further assumed that the mass transpiration velocity is v w (x) = − √ α f x −3/4 f wa where υ f signifies the respective kinematic viscosity and f wa the mass injection/suction parameter. Furthermore, D f identifies the mass diffusion coefficient, E a the activation energy, and k 2 r (x) the variable chemical reaction rate which is delineated later mathematically. Additionally, the final term Q on the RHS (right hand side) of Equation (4) demonstrates the role of the erratic heat source/sink, which can be illustrated as [see [46]]: where u r (x) = 4x −1/2 is the reference velocity which is occupied by Raees et al. [47] and the temperature-dependent heat sink/source and the exponentially decaying space coefficients are distinguished by B * a and A * a , respectively. Consequently, the positive values of A * a and B * a correspond to the phenomenon of a heat generation source factor while the negative values of A * a and B * a correspond to the phenomenon of the heat generation sink factor. In the aforementioned stated equations, the mathematical representations exercised for the nanofluid (NF) are heat capacity ρc p n f , thermal conductivity k n f , electrical conductivity σ n f , dynamic viscosity µ n f , and the requisite posited density ρ n f . The NF model expressions or correlations thermo-physically can be written as follows: ρc p n f Hence, the aforementioned leading equations exemplify the thermo-physical properties of the kerosene-oil-based graphene oxide nanofluid, where ϕ indicates the volume fraction of the nanoparticles. Additionally, the subscripts snp and f are the nanoparticles and the base working fluid (KO). Here, the shape factor is taken as a sphere whose value is 3. Generally, the data of the graphene oxide (GO) nanoparticles along with the regular-based fluid (KO) are provided in Table 1.

Similarity Variables
To ease the examination of the considered wall jet flow phenomenon, we stated the subsequent self-similarity variables according to Glauert [39] which are as follows: where ψ is the stream function and signified as u = ∂ψ/∂y and v = −∂ψ/∂x. Thus, the velocity component is calculated more simply in the form as follows:

Momentum Similarity Equation
For the momentum similarity equations, it is necessary to exercise the similarity variables of Equation (11) early along with the derived velocity components of Equation (12) into the governing Equations (1) and (2) which holds the continuity equation directly while the momentum equation is reduced to the form as: is called the magnetic parameter.

Energy Similarity Equation
It is better to write the state of Equation (3) in a more simplified form, therefore, we incorporate the expression of the non-uniform heat source/sink Q term to obtain the form as: In addition, substituting the respective similarity transformations (11) into the Equation (14) reduces to the following requisite posited dimensionless form as follows: where Pr = υ f α f is called the Prandtl number.

Concentration Similarity Equation
Before moving towards the working procedure of the concentration similarity equation, it is helpful to first define the logical term in Equation (4). The second-last term on the RHS of Equation (4) represents the thermophoretic particle deposition which is symbolically denoted by V TPD and mathematically expressed as: where T a is a reference temperature and k a is the thermophoretic coefficient with a range of values from 0.2 to 1.2 as designated by Batchelor and Shen [28]. Moreover, to obtain a similarity solution of the concentration equation after the use of variables, here, we express the variable reaction rate as k 2 r (x) = k 2 0 x −3/2 . Additionally, by plugging this equation and Equation (16) along with the self-similarity Equation (11) into the acquired leading governing Equation (4). Hence, one obtains the form as: where are called the Schmidt number, the chemical reaction rate, the temperature difference parameter, the activation energy parameter, and the thermophoretic parameter, respectively.

Transformed Boundary Conditions
According to the use of similarity Equation (11) in Equation (5), we obtain Here, f wa signifies the constant mass suction/injection factor with ( f wa > 0) for the mass suction and ( f wa < 0) for the mass injection.

Gradients
In the present study, the following key gradients (the skin friction coefficient C f , the Nusselt number (Nu x ), and the Sherwood number (Sh x )) are expressed as: where k n f and µ n f are the corresponding thermal conductivity and viscosity of the nanofluid. Now, employing the self-similarity transformations in the above stated Equation (19), we obtain where Re x = u r x υ f is called the local Reynolds number.

Numerical Methodology
This portion of the examination describes the complete working procedure (WP) as well as the reliability, confirmation, or validation of the given scheme. Initially for the WP, it is better to exercise the set or structure of the acquired bounded ordinary differential Equations (13), (15) and (17), along with the boundary restrictions (18). These equations were computationally simulated using the bvp4c solver, which relies on the 3-stage Lobatto IIIA algorithm or formula. In addition, there are several techniques such as RKF-45 (Runge-Kutta-Fehlberg fourth fifth order shooting method), shooting, and Keller box. However, the bvp4c is simple to use and has better accuracy. In the bvp4c technique, the structure of higher-order ODEs is altered into first-order ODEs to apply this technique. Let the variables be: and the changed initial conditions are In the present study, tolerance ε = 10 −6 is applied throughout the computation to meet the convergence criterion. Further, the proposed simulation technique needs initial or preliminary estimates to satisfy the boundary condition (23). In other geometry problems, the initial assumption for the first solution is fairly straightforward, but in the case of the wall jet nanofluid flow problem, it is a little more challenging. As a response, it is preferable to select the values of the other constant parameters that are pertinent to obtaining the necessary results quickly or easily. The computation range of numerical integration is ξ max = 25, which is sufficient for the graphical results to meet the requirement of boundary asymptotically.

Numerical Authentication of the Scheme
This aspect of the work illustrates the crucial aspect of the acknowledged accuracy, dependability, authenticity, or validation of the regarded numerical code. Therefore, it is preferable to create a graph or build a numerical table so that we can compare the new work with the published work for some limiting cases to verify the given code. Here, in this study, a possible Table 2 was prepared for the numerical values of shear stress with varying ϕ when f w = 0, and M a = 0. From the numerical table, it can be shown that the single current output of shear stress matched with the previously published works of Glauert [39] and Waini et al. [49] extremely accurately. Henceforth, the recent and earlier research are very well matched, therefore, we are confident that the wall jet flow problem will be solved using the suggested technique.

Analysis of the Results
This section breakdown was focused on interpreting the whole set of data for the dynamics of the wall jet heat and mass transfer characteristics flow of the kerosene-oil-based graphene oxide nanoparticles due to the influence of many unique parameters. The given problem comprised different control influential parameters such as the solid nanoparticle volume fraction ϕ, the magnetic parameter M a , the thermophoretic parameter Σ a , the radiation parameter R d , the activation parameter E a , the chemical reaction rate parameter β a , the dimensionless temperature difference parameter δ a , and the mass transpiration parameter f wa . The effects of these parameters on the velocity profile F (ξ), temperature profile G(ξ), concentration profile S(ξ), shear stress µ n f /µ f F (0), heat-transfer rate − k n f /k f + (4/3)R d G (0), and mass transfer rate S (0), are shown in Figures 2-14 as well as in Tables 3-5. For numeric calculations, the fixed values of the distinguished constraints were the following: ϕ = 0.025, m = 1, M a = 0.05, Σ a = 0.5, f wa = 0.5, Sc a = 10, β a = 0.5, δ a = 0.5, E a = 0.5, and R d = 2.0. In all graphs, the outcomes were presented in the form of solid gray, blue, and red colors.                                               Tables   Table 3 shows the quantitative outcome of µ n f /µ f F (0) for the single branch solutions with variations in ϕ, M a , and f wa . From the outcomes, it is observed that the µ n f /µ f F (0) escalates with superior values of ϕ,M a , and f wa while it shrinks with f wa . Moreover, the magnitude of the largest and smallest shear stress values was computationally noticed due to the mass suction and injection parameter, respectively. Furthermore, the shear stress increases up to 12.3% due to the presence of nanofluid. On the other hand, the numerical values of the heat and mass transfer rates of the wall jet kerosene-oil-based graphene oxide nanoparticles were created or quantitatively constrained in Tables 4 and 5, respectively. From Table 4, it is seen that the heat-transfer rate heightens with the higher value of ϕ, and the generation of the heat sink parameter A * a , B * a < 0 but is constant throughout the domain with the precise changed value of the temperature index parameter m. Meanwhile, the impact of the heat source parameter A * a , B * a > 0 decreases the heat-transfer rate. Furthermore, the numerical outputs of the mass transfer rate for the various influence parameters were highlighted in Table 5. Here, the mass transfer rate continuously declines with the superior values of ϕ, β a , δ a , and Σ a but augments with a higher Sc a . In addition, the largest mass transfer rate values were observed for the values of the thermophoretic parameter. It is worth mentioning that the Nusselt number increases up to 0.119% and the Sherwood number increases up to 0.006% due to the presence of nanofluid. Whereas the Nusselt number decreases up to 1.328% in the case of suction and increases up to 1.188% due to injection. On the other hand, the Sherwood number decreases up to 0.182%, 0.021%, and 0.022% due to Σ a , δ a , and β a , respectively. Whilst the Sherwood number increases up to 33.65% due to the Schmidt number.

Graphical Interpretation of the Velocity Profiles
Figures 2-4 exemplify the dimensionless velocity profile F (ξ) of the water-based graphene oxide nanoparticles for the solutions due to the influence of ϕ, M a , and f wa . Due to the cumulative effect of the influencing parameters, it is inferred from the results that the graphical behavior of the velocity profiles asymptotically converges and satisfies the boundary requirements or constraints. To explain more explicitly, the outcomes signpost in Figure 2 that the velocity field curves initially decays and then slightly improves with the growing value of ϕ. In other words, the wall jet velocity profile behaves distinctly versus the specific range of the pseudo-similarity variable ξ owing to the higher impressions of ϕ. In addition, the momentum boundary layer thickness decreases with a higher value of ϕ. Generally, the solid nanoparticle volume fraction improves the sign or evidence of viscosity, as a result, the motion of the wall jet flow reduces. On the other hand, the velocity field curves initially shrink and hasten with the superior values of M a and f wa as shown in Figures 3 and 4, respectively. Later, the graphical outcomes in both profiles behave inversely in the approximate ranges of (5 ≤ ξ ≤ 11) and 11 ≤ ξ ≤ 22 with the larger values of M a and f wa . Generally, the larger impact of magnetic force creates the drag-type force called the Lorentz force. This force reaction can stop the motion of the wall jet flow, as a response, the profile of velocity decreases. Furthermore, the gap for the curves with the higher impacts of f wa is slightly improved compared to the curves drawn for the value of M a .

Interpretation of the Temperature Profiles
The stimulus of the solid volume fraction of nanoparticles ϕ and the thermophoretic particle deposition parameter Σ a on the temperature G(ξ) of the water-based graphene oxide nanomaterials flow are respectively highlighted in Figures 5 and 6. Figure 5 shows that the thermal boundary layer thickness and the temperature profile curves monotonically rise with the higher impacts of ϕ. Generally, the behavior was spotted due to the wellknown facts. The higher impacts of a nanoparticle volume fraction creats a larger heat which significantly boosted up the values of the thermal conductivity. Therefore, the advanced thermal conductivity boosted the profile of temperature as well as the thickness of the thermal boundary layer (TTBL). Conversely, the temperature and the TTBL diminished with the superior effects of the thermophoretic parameter as shown graphically in Figure 6. Moreover, the thermophoretic parameter produces a much lower temperature because of the minimum thermal conductivity. Figures 7 and 8 demonstrate the temperature profile curves with the influence of the heat generation source and sink parameters, respectively. The outcomes reveal that temperature profile curves and the TTBL augmented with the higher impact of internal heat generation source parameter,A * a , B * a > 0, while it shrinks due to the generation of heat sink constraint A * a , B * a < 0. Physically, by adding more heat to the wall jet flow system, the heat source's presence ultimately boosts the fluid temperature. As a result, the temperature profile rises. While, in terms of physics, the heat sink factor that captivates heat from the temperature boundary layer causes the temperature profile to drop. Moreover, the TTBL rises and declines with the internal heat source and sink influences, respectively.

Interpretation of the Concentration Profiles
Figures 9-14 illustrate the dimensionless concentration profile S(ξ) of the keroseneoil-based graphene oxide nanoparticles flow for the branch of single outcomes with the influence of A * a , B * a , ϕ, Sc a , Σ a , and E a . Since the relevant boundary conditions are observed to be satisfied in every graph, the behavior of the solutions asymptotically converges when the many notable parameters are varied. Figures 9 and 10 represent the influences of the erratic heat source/sink factors on the concentration profile, whereas the concentration declines with the heat source constraint but accelerates with the internal heat sink factor. Moreover, the impact of such a parameter behaves oppositely as compared to the temperature profile. The impact of ϕ on the concentration profile curves is shown in Figure 11, where the behavior of the profile is growing. This behavior is seen due to the higher impressions of nanoparticles because they cause a significant improvement to the kerosene-oil-based graphene nanoparticles that augment the concentration and thickness of the concentration boundary layer. Moreover, the effect of Sc a and Σ a on the concentration profile is depicted in Figures 12 and 13, respectively. With the superior values of Sc a and Σ a , the concentration and TCBL (thickness of the concentration boundary layer) enriches. This behavior is typically brought on by larger repercussions of the parameter Σ a , which might enhance the thermophoretic coefficient k a , and, consequently, improve the concentration profile as well as the TCBL. Finally, the E a impact on the concentration field is graphically shown in Figure 14. Since the outcome is changed due to this term exp(−E a /(1 + δ a G))S, it is a better tool showing how the activation energy affects the nanoparticle concentration profile. Hence, with the larger E a , the TCBL and profile of nanoparticle concentration upsurge.

Conclusions
The aim of this paper was to inspect the effect of thermophoretic particle deposition on 2D wall jet heat and mass transfer flow that carries kerosene oil-based graphene oxide nanoparticles with thermal energy and mass transpiration velocity. Further, the Arrhenius activation energy and chemical reaction along with Lorentz forces were also added to the problem. The Glauert model transformations were used to renovate the governing equations into a similarity form of ODEs, which were then numerically solved using the bvp4c method. The key findings of the current investigation are:

•
With the larger value of the solid nanoparticle volume fractions, the velocity of the wall jet flow profile is moderate, while the temperature and concentration profile curves are augmented due to the higher number of nanoparticles, consequences which generate greater heat and significantly improved the thermal conductivity.

•
For the advanced value of the thermophoretic parameter, the temperature profile shrinks but the profile of the concentration is boosted which might enhance the thermophoretic coefficient.

•
The internal heat source factor escalates the temperature profile distribution but decelerates the concentration profiles while the influence of the internal heat sink factor is the opposite. By adding more heat to the wall jet flow system the heat source ultimately boosts the fluid temperature, whereas the heat sink factor that captivates heat from the temperature boundary layer causes the temperature profile to drop. • Initially, the velocity profile enriches and then declines with the advanced value of the mass suction parameter.

•
The shear stress progresses with the larger value of the magnetic parameter due to the Lorentz forces while it reduces magnitude-wise with the mass injection parameter.

•
The solid nanoparticle volume fractions escalate the heat transfer while reducing the rate of mass transfer.

•
The shear stress increases up to 12.3%, the Nusselt number increases up to 0.119%, and the Sherwood number increases up to 0.006% due to the presence of nanofluid.

•
The Nusselt number decreases up to 1.328% in the case of suction and increases up to 1.188% due to injection.

•
The Sherwood number decreases up to 0.182%, 0.021%, and 0.022% due to Σ a , δ a , and β a , respectively.  Derivative with respect to ξ