Signiﬁcance of Temperature-Dependent Density on Dissipative and Reactive Flows of Nanoﬂuid along Magnetically Driven Sheet and Applications in Machining and Lubrications

: Nanoﬂuid lubrication and machining are challenging and signiﬁcant tasks in manufacturing industries that are used to control the removal of a material from a surface by using a cutting tool. The introduction of a nanoﬂuid to the cutting zone provides cooling, lubricating, and chip-cleaning beneﬁts that improve machining productivity. A nanoﬂuid is a cutting ﬂuid that is able to remove excessive friction and heat generation. Chemical reactions and temperature-dependent density are essential in the thermal behavior of a nanoﬂuid. The present study presents a careful inspection of the chemical reactions, temperature-dependent density, viscous dissipation, and thermophoresis during the heat and mass transfer of a nanoﬂuid along a magnetically driven sheet. The physical attitude of viscous dissipation and the chemical reaction improvement rate in magneto-nanoﬂuid ﬂow is the primary focus of the present research. By applying the proper transformation, nonlinear partial differential expressions are introduced to the structure of the ordinary differential framework. The ﬂow equations are simpliﬁed into nonlinear differential equations, and these equations are then computationally resolved via an efﬁcient computational technique known as the Keller box technique. Flow factors like the Eckert number, reaction rate, density parameter, magnetic force parameter, ther-mophoretic number, buoyancy number, and Prandtl parameter governing the velocity, temperature distribution, and concentration distribution are evaluated prominently via tables and graphs. The novelty of the current study is in computing a heat transfer assessment of the magneto-nanoﬂuid ﬂow with chemical reactions and temperature-dependent density to remove excessive friction and heating in cutting zones. Nanoﬂuids play signiﬁcant roles in minimum quantity lubrication (MQL), enhanced oil recovery (EOR), drilling, brake oil


Introduction and Literature Review
Nanofluids were initially created to improve the performance of minimum quantity lubrication (MQL) and overflow cooling.The particles of manufactured nanomaterials suspended in base fluids are known as nanofluids.MQL is the greatest choice as a base fluid for nanotechnologies because this fluid is currently popular for its high efficiency compared to traditional overflow cooling fluids.The idea of nanofluid lubrication was established as a result of developing nanotechnologies.Using nanofluid substituting for frequently employed cutting fluids, a procedure known as nano-lubrication, can efficiently cool and lubricate cutting regions during manufacturing operations.A substance's specified dimension, shape, and surface polish are obtained with the elimination of additional substances in the form of tiny pieces during cutting, which is now thought to be the best possible mechanical technique.The results of the present research can improve our awareness of nanofluid EOR procedures through the angle of oil-based substances, nanoparticles, and interfacial dynamics.They can also help us to figure out how nanofluids perform in multiphase fundamental movement modeling and assessments.For the majority of machining tools, an increased efficiency is anticipated to be mostly based on the significantly and substantially strong temperature-sensitive thermal conduction of nanofluids.Cutting tools are required to eliminate additional substances by making direct frictional contact.Due to its significance across multiple manufacturing and scientific procedures, including liquid film vaporization, the aerodynamic protrusion of polymeric spreadsheets, the formation of crystals, the cooling of solid papers, the design of manufacturing instruments and several heat exchange systems, and glass and polymeric sectors, the examination of the flow within a drastically stretched sheet is of significant importance.Khan et al. [1] used computational methods to examine the effects of solar radiation, heat production, and chemical processes on a stationary flow separation through a wedge-like structure in electromagnetic water.Magneto-nanofluid motion caused by an elongated surface was discussed by Anjum et al. [2].In the context of chemical reactions and radiation temperature factors, a researcher [3] explored the magnetic flow separation of tiny liquids and heat transmission through a dynamic stretchable surface.Zhang et al. [4] looked at the impacts of a plant-based oil's thickness, friction, and molecular composition on its ability to cool and its lubricated qualities.The radiative heat transmission of nanofluids across a wedge-shaped structure that is receptive to initial-order chemical reactions, temperature generation absorption, and vacuum forces was investigated by Kasmani et al. [5].Doshi et al. [6] discussed the difficulties involved in using a nanofluid as a cutting fluid for machine work.The impact of Darcy's homogenous chemical processes on higher-order nanofluids was examined by Jagadha et al. [7].
During either the initial or subsequent oil extraction operations, reserved oil can often be recovered via enhanced oil recovery (EOR).Chemical EOR has shown to achieve excessive oil extraction with minimal operating expenses when compared with traditional techniques, and nanofluids have drawn a lot of interest due to their beneficial characteristics of affordability, excessive oil extraction, and broad applicability.Additionally, a number of research findings have concentrated on the function of nanoparticles in the nanofluid EOR procedure.The wetting capability of materials (surface/particle micromechanics), the stiffness at the water-oil interface (massive fuel/molecule micromechanics), and the stress at the mechanically disjoining interface (heavy fuel/molecules/particle micromechanics) can all be affected by nanofluids.Moreover, they may reduce thickness (micromechanics of massive oil nanoparticles).Furthermore, nanomaterials are capable of sticking to the surfaces of water and oil, changing the surface tension and causing an emulsion to develop.Rao et al. [8] calculated the rates of heat transmission for a basic liquid and a mixture with nanoscale integration, and an analysis of finite elements was used to generate the heating ranges of a cutting instrument.Afify [9] quantitatively explored the effects of numerous slips with viscosity dissipations on the flow separation and heat transmission of unconventional nanoparticles across an extensible sheet.Hasin et al. [10] combined several different kinds of nanomaterials with a supporting liquid to enhance the machining and cutting abilities of ordinary plant-based oil.Srikant et al. [11] looked into various methods to substitute the deployment of cutting liquids for minimal lubrication.The significance of nanotechnology in machining procedures was studied by Kadirgama [12].The topic was covered in many sections, including the cutter effect, the roughness on the surface, tool durability, and tool corrosion.
The use of injecting liquids, nano-foam engineering, and nano-emulsion manufacturing is crucial for the enhanced oil recovery (EOR) procedure.The nanoparticle absorbing reaction serves as the foundation for the nanofluid EOR system.When a nanomaterial attaches to a metallic substrate, the wetness of the material changes from a remaining oil-wet state to moist with water.Additionally, a nanomaterial is capable of attaching to the outermost layers of water and oil, changing the surface tension and causing emulsions to develop.There is a greater need for brake devices with better and more effective thermal dissipation methods and characteristics to decrease drag forces to enhance automotive aerodynamic efficiency.Nanofluids are crucial for minimizing resistance and maximizing the transmission of heat to resolve this problem.The current mechanism is crucial for dual-phase flows, thermal exchange optimized performance, fuel cell development, and geothermal power installations.Entropy formation in a streaming liquid has been the topic of many investigations in the field of science.Li and Ali [13] provided a thorough discussion of the many facets of nanofluids, including their characterization, movement, and energy transfer procedures.The approaches for characterizing the several characteristics of nanofluids, such as their thermal conductivity and the variables affecting it, were presented.Sharma et al. [14] summarized several significant peer-reviewed papers on the use of nanofluids in cutting, tapping, crushing, and spinning operations.To improve thermal and heat transmission, a nanofluid is used as a cutting fluid in machining.For enhancing efficiency, micro-lubrication methods and method modification play similar roles [15].During machining, manufacturing heat is an essential occurrence.Rifat et al. [16] employed various cooling systems and tool durability to eliminate the heat-related negative consequences of skin burning.Sidik et al. [17] discussed the most current development and uses of nanofluids in machining operations.They reviewed numerous traditional and modern cooling strategies used in machining and to improve thermal conductivity and studied the characteristics of nanofluids.Amin et al. [18] examined the impact of nanofluid flow behavior in machining procedures.An experimental study on the effectiveness of environmentally friendly plant-based oil nanofluids in rotating movement was conducted by Pasam et al. [19].During machining, the cutting levels, cutting pressures, tool wear, and surface roughness were determined for stable cutting circumstances.For all interatomic connections, Lautenschlaeger [20] employed the Lennard-Jones shortened and shifting formulation.For the accuracy and justification of results, the numerical outcomes of heat transfer were compared with the existing literature in [21][22][23].The effects of temperaturedependent density or variable density on the heat transfer assessment of viscous and nanofluid flow along various geometries was reported in [24][25][26][27][28][29][30][31][32][33][34][35][36][37][38][39][40] and published in wellreputed journals.Dolatabadi et al. [41] performed a heat transfer improvement of engine oil by assuming the temperature-dependent viscosity.
The originality of the current work is to resolve problems in high-temperature systems.Most studies are found to deal with a small temperature difference between the surface and ambient nanofluid.However, circumstances arise where this temperature difference is high.In this situation, the density is assumed as an exponential function of temperature, and the nanofluid is assumed to be electrically conducting in the sense of ionized nanofluid due to the high operating temperature.The nanofluid density is assumed to be reduced exponentially with temperature and the relation for this mechanism.The magnetic field is applied on the surface to decrease excessive heat through tool chips and tool works during the machining procedure in the cutting region.The nanofluid works as a cutting fluid to lubricate the surface to prevent extreme temperatures.The introduction of an electrically conducting nanofluid to the cutting zone provides cooling, lubricating, and chip-cleaning benefits that improve machining productivity.Several scientists explored the temperature-dependent density flows in various fluid geometries ([2,3,9,24-40]).The novelty of the present analysis is to produce cooling in cutting tools in the presence of nanofluid lubrication.The impact of a chemical reaction, viscous dissipation, and temperature-dependent density on a magneto nanofluid along a stretching surface is explored.The thermal and concentration boundary layer of a nanofluid is examined along a stretching sheet numerically.The present work is original and significant because it satisfies the given boundary conditions.For the originality of the current research, the numerical findings of the heat and mass transfer are compared with the published research.

Mathematical Formulation and Flow Shape
The research topic is based on a computational analysis of the effects of variable density, chemical reaction, and viscous dissipation on the heat transfer and nanoparticle characteristics of a nanofluid over a stretching sheet under thermophoretic and Brownian motion.The current problem is converted into a set of partially differentiated conditions through stream-function modifications and then transformed into an ordinary form.The constructed problem is resolved through the integration of the Keller box technique with the finite difference method (FDM) procedure.The logical outcomes for the physical characteristics under consideration will be represented visually and in tabulated form.
Consider the potentially magnetically driven two-dimensional nanomaterial shown in Figure 1, where u and v define the velocity distributions in x and y directions, respectively.The axes x and y are the directions along the surface and normal to the stretching surface.The temperature is expressed by T, the stream-free temperature is denoted by T ∞ , the fluid thermal conductivity is indicated by κ, and the particular heat is expressed by C p .The thickness is marked by ρ, the density of the particles is marked by ρ p , the kinematics viscosity of the fluid is marked by ν = µ/ρ, the acceleration due to gravitation is marked by g, and the fluid density is marked by ρ f .The mathematical models for the electromagnetic force, momentum, energy, and continuum are as follows [1][2][3]: The boundary conditions of the current methodology are The required boundary parameters for both the temperature and elements of velocity are described in the section above.P signifies pressure, ρ is for the base fluid density, α is for the heating diffusion, ν is for the kinematic viscosity, a is for a certain constant, D B is for the Brownian distribution, D T is for the thermophoretic distribution, and τ = (ρc)p/(ρc) f is for the proportion of the efficient heat ability of the nanoparticle substance to the fluid heat ability, where ρ appears for density and c is for the factor of heat conductance.and  = ()/() is for the proportion of the efficient heat ability of the nanoparticle substance to the fluid heat ability, where  appears for density and  is for the factor of heat conductance.

Stream Factors and Similarity Analysis
The stream factors and similarity factors for the conversion of partial equations into ordinary equations are expressed below [5,7,9].
The stream factors for concentration, velocity, and temperature distribution are expressed below.

Stream Factors and Similarity Analysis
The stream factors and similarity factors for the conversion of partial equations into ordinary equations are expressed below [5,7,9].
The stream factors for concentration, velocity, and temperature distribution are expressed below.
The convenient terms of the modified ordinary equations are evaluated in Equations ( 9)- (11). Here, is designated as the Brownian parameter, θ is represented as the non-dimensional temperature, Le = v D B is defined as the Lewis parameter, η is indicated as the similarity factor, Pr = ν α is denoted as the Prandtl number, σ = k r 2 a is the reaction rate, and ν = µ ρ is denoted as the kinematical viscosity.The appropriate boundary values are described in Equation (12).

Numerical Procedure and Method
The nonlinear ordinary differentiated form is reduced in the algebraic system of equations by offering the Keller box method (KBM) under defined boundaries.For smooth coding through MATLAB, the Newton-Raphson Scheme (NRC) is applied to alter the system of Equations ( 9)- (12) in the matrix notation.

Results and Discussion
The relationship between the chemical reaction, heat density, and viscous dissipation in a laminar, two-dimensional, mixed convective flow of a magneto nanofluid over a stretching sheet is shown in the present investigation.The fluid density is assumed as an exponential component of temperature in the current model.By using the variable density, the solutions are valid across any significant temperature variation.The boundary layer equations are developed for the present mechanism under defined boundary conditions.It is important to note that the dimensionless governing mathematical equations are additionally simplified to a group of non-similar calculations with the support of the stream function structure and are numerically incorporated through the Keller box approach to compare the findings with those of previous investigations.The numerical findings from this technique are discussed above geometrically with acceptable agreement.The governing flow factors like the Eckert number, density parameter, magnetic force parameter, thermophoretic number, buoyancy number, and Prandtl parameter on the velocity, temperature distribution, and concentration distribution are evaluated prominently via tables and graphs.
The impact of the reaction rate σ on the velocity, temperature distribution, and concentration distribution is depicted in Figure 2a-c.The velocity graph increases as the reaction rate increases with good variations.Similarly, the temperature attains maximum magnitude as the reaction rate increases in the presence of entropy generation.The mass distribution decreases as the reaction rate decreases.The prominent variations in each plot are noted.The consequences of the thermal density parameter n on the velocity, temperature distribution, and concentration distribution are plotted in Figure 3a-c in the presence of entropy generation.The prominent variations in the velocity distribution are examined for each n.The velocity of a water-based fluid increases for a lower density but decreases as n increases.The certain amplitude in the temperature distribution is illustrated as n = 0 but decreases as the density parameter increases.Similarly, the concentration distribution increases as the density of fluid decreases for Pr = 7.0.The minimum concentration is evaluated for the maximum density in the presence of entropy generation.Physically, the buoyant pressure increases for higher values of the density factor, which leads to growth density variations with temperatures.In Figure 4a, the impact of the Eckert number Ec on the velocity distribution yields an enhancement with feasible asymptotic behavior with n = 0.4.The temperature of the water-based fluid increases as the Eckert number increases with prominent variations in Figure 4b.Physically, it is valid because the entropy generation acts like a heat source to produce the ability in the temperature of the fluid.The concentration φ is increased for a lower entropy generation and decreases as the Eckert number is enhanced in Figure 4c.It can be observed that each physical profile yields prominent variations.The influence of the Prandtl number on the velocity distribution, concentration, and temperature distribution is evaluated in Figure 5a-c for Pr = 0.7, 1.0, 3.0, 5.0, and 7.0 with some fixed quantities.The Prandtl number has a significant impact on the flow rate through the plotted profiles.The velocity of the fluid boosts for a lower Pr and gradually decreases as the Pr increases.The prominent enhancement in the temperature of the fluid is assessed for the maximum Pr but decreases as the Pr decreases.Moreover, a reasonable change in the concentration distribution is evaluated for each Pr with entropy generation.According to scientific logic, this is valid because the entropy growth significantly boosts the thermal energy of the fluid, which increases the fluid's temperature.Additionally, the influence of the magneto force ξ is drafted in Figure 6a-c  because the entropy generation acts like a heat source to produce the ability in the temperature of the fluid.The concentration  is increased for a lower entropy generation and decreases as the Eckert number is enhanced in Figure 4c.It can be observed that each physical profile yields prominent variations.The influence of the Prandtl number on the velocity distribution, concentration, and temperature distribution is evaluated in Figure 5a-c for Pr = 0.7, 1.0, 3.0, 5.0, and 7.0 with some fixed quantities.The Prandtl number has a significant impact on the flow rate through the plotted profiles.The velocity of the fluid boosts for a lower Pr and gradually decreases as the Pr increases.The prominent enhancement in the temperature of the fluid is assessed for the maximum Pr but decreases as the Pr decreases.Moreover, a reasonable change in the concentration distribution is evaluated for each Pr with entropy generation.According to scientific logic, this is valid because the entropy growth significantly boosts the thermal energy of the fluid, which increases the fluid's temperature.Additionally, the influence of the magneto force  is drafted in Figure 6a-c         In addition, the impact of the thermophoretic number  on the physical features of heat and mass transmission across the heated surface is assessed in Figure 8a,b.By increasing , the heat and mass transfer rates are gradually enhanced with prominent changes in the presence of thermal density and entropy generation.On the other hand, minimum values of  decrease the heat and mass transmission rates with a higher Pr because the heating conductivity of a liquid with a greater Prandtl quantity is often lesser, which decreases the conductions and improves the amount of heat ′ transmission at the surface's interface.In Figure 9a,b, the impact of the Brownian motion parameter  on the mass and heat transmission is evaluated in the presence of thermal density and entropy generation.It can be examined that the dimensionless heat ′ transmission declined for the maximum  with Pr =7.0.However, the dimensionless mass transfer is increased for the maximum  under buoyant and magnetic attractions.The combination of the density influence and the magnetic attraction serves as a preventing power.The thermal efficiency can be controlled by this preventing power, which has a variety of uses, including the magnetic covering of cables and metallic material and the development of magnetohydrodynamic energy.In addition, the impact of the thermophoretic number Nt on the physical features of heat and mass transmission across the heated surface is assessed in Figure 8a,b.By increasing Nt, the heat and mass transfer rates are gradually enhanced with prominent changes in the presence of thermal density and entropy generation.On the other hand, minimum values of Nt decrease the heat and mass transmission rates with a higher Pr because the heating conductivity of a liquid with a greater Prandtl quantity is often lesser, which decreases the conductions and improves the amount of heat θ transmission at the surface's interface.In Figure 9a,b, the impact of the Brownian motion parameter Nb on the mass and heat transmission is evaluated in the presence of thermal density and entropy generation.It can be examined that the dimensionless heat θ transmission declined for the maximum Nb with Pr = 7.0.However, the dimensionless mass transfer is increased for the maximum Nb under buoyant and magnetic attractions.The combination of the density influence and the magnetic attraction serves as a preventing power.The thermal efficiency can be controlled by this preventing power, which has a variety of uses, including the magnetic covering of cables and metallic material and the development of magnetohydrodynamic energy.
entropy generation.It can be examined that the dimensionless heat ′ transmission declined for the maximum  with Pr =7.0.However, the dimensionless mass transfer is increased for the maximum  under buoyant and magnetic attractions.The combination of the density influence and the magnetic attraction serves as a preventing power.The thermal efficiency can be controlled by this preventing power, which has a variety of uses, including the magnetic covering of cables and metallic material and the development of magnetohydrodynamic energy.In Tables 1 and 2, the computational and numerical calculations of the dimensionless skin friction, heat rate, and mass transfer are computed for diverse factors of the buoyancy parameter  and the magneto force parameter  through the stretched surface with temperature-density and entropy effects.The enhancement in friction and Sherwood coefficients is evaluated for the maximum buoyancy  = 5.0, but the heat transfer is decreased.From a physical standpoint, it was anticipated that the stress variations in the nanofluid nanoparticles would be enhanced by the larger buoyant attraction.The impact of the magneto force  on the numerical features of heat rate, skin friction, and mass transmission is shown in Table 2.The maximum outcomes of each physical property are determined for higher values of the magneto force.The dimensionless heat and mass rate decreases as the magnetic force decreases.Tables 3 and 4 show the comparable results of the Prandtl number as well as the Lewis variable with the existing works in the literature, including Ibrahim [21], Khan et al. [22], and Dawar et al. [23], with prominent agreement.In Tables 1 and 2, the computational and numerical calculations of the dimensionless skin friction, heat rate, and mass transfer are computed for diverse factors of the buoyancy parameter λ and the magneto force parameter ξ through the stretched surface with temperature-density and entropy effects.The enhancement in friction and Sherwood coefficients is evaluated for the maximum buoyancy λ = 5.0, but the heat transfer is decreased.From a physical standpoint, it was anticipated that the stress variations in the nanofluid nanoparticles would be enhanced by the larger buoyant attraction.The impact of the magneto force ξ on the numerical features of heat rate, skin friction, and mass transmission is shown in Table 2.The maximum outcomes of each physical property are determined for higher values of the magneto force.The dimensionless heat and mass rate decreases as the magnetic force decreases.Tables 3 and 4 show the comparable results of the Prandtl number as well as the Lewis variable with the existing works in the literature, including Ibrahim [21], Khan et al. [22], and Dawar et al. [23], with prominent agreement.

Conclusions
The relationship between the chemical reaction, heat density, and viscous dissipation in a laminar, two-dimensional, mixed convective chemically reactive flow of a magneto nanofluid over a stretching sheet is shown in the present investigation.The fluid density is assumed to be an exponential component of the temperature rather than a linear component as in the Boussinesq model.Because of this, the solution to the issue remains valid across any significant temperature variations, and the outcomes are more precise.The boundary layer equations were developed for the present mechanism under defined boundary conditions.It is important to note that the dimensionless governing mathematical equations were additionally simplified to a group of non-similar calculations with the support of the stream function structure and were numerically incorporated through the Keller box approach to compare the findings with those of previous investigations.The numerical findings obtained from the technique discussed above were geometrically evaluated, and it was discovered that they are in acceptable agreement.The effects of governing flow factors like the Eckert number, density parameter, magnetic force parameter, thermophoretic number, buoyancy number, and Prandtl parameter on the velocity, temperature distribution, and concentration distribution were evaluated prominently via tables and graphs.The prominent key outcomes are described as follows:

•
The certain amplitude in the temperature distribution is illustrated as n = 0 but decreases as the density parameter increases.Similarly, the concentration distribution increases as the density of fluid decreases for Pr = 7.0.

•
The temperature of the water-based fluid increases as the Eckert number increases with prominent variations.Physically, it is valid because the entropy generation acts like a heat source to produce the ability in the temperature of the fluid.

•
The prominent enhancement in the temperature of the fluid is assessed for the maximum Pr but decreases as the Pr decreases.Moreover, the reasonable change in the concentration distribution is evaluated for each Pr with entropy generation.

•
The heat and mass transfer rates are gradually enhanced with a prominent change as Nt increases under thermal densities and viscous dissipations.

•
It is found that the dimensionless Nusselt coefficient is decreased for the maximum Nb under Pr = 7.0.However, the dimensionless mass transfer is increased for the maximum Nb in the presence of buoyancy and magnetic forces.
with the entropy generation and density variations.It is assessed that the velocity of the fluid and temperature decreases as ξ increases.The prominent enhancements in both the temperature and velocity are depicted for lower values of ξ.On the other hand, the concentration distribution enhances as ξ is reduced.It is normally understood because a boost in the Prandtl value causes a reduction in the fluid's thermal conduction, which, in turn, reduces the amount of frictional force within the viscosity layers.The numerical results of the velocity f (η), temperature θ, and concentration φ are drafted in Figure7a-cfor various choices of the Brownian motion Nb under the impacts of viscous dissipation and the magnetic field across the stretched sheet for a high temperature difference.
with the entropy generation and density variations.It is assessed that the velocity of the fluid and temperature decreases as  increases.The prominent enhancements in both the temperature and velocity are depicted for lower values of .On the other hand, the concentration distribution enhances as  is reduced.It is normally understood because a boost in the Prandtl value causes a reduction in the fluid's thermal conduction, which, in turn, reduces the amount of frictional force within the viscosity layers.The numerical results of the velocity  ′ (), temperature , and concentration  are drafted in Figure7a-cfor various choices of the Brownian motion  under the impacts of viscous dissipation and the magnetic field across the stretched sheet for a high temperature difference.

Figure 8 .
(a,b) Graphical behavior of Nusselt and Sherwood values for .

Figure 9 .
Figure 9. (a,b) Graphical behavior of Nusselt and Sherwood values for Nb.

Table 1 .
The computational and numerical calculations for different choices of

Table 2 .
The computational and numerical calculations for different choices of

Table 1 .
The computational and numerical calculations for different choices of λ.

Table 2 .
The computational and numerical calculations for different choices of ξ.
Author Contributions: Conceptualization, Z.U. and A.H.; methodology, Z.U.; software, Z.U.; validation, Z.U.; formal analysis, M.S.A.; investigation, N.H.A.; resources, S.S.; data curation, M.S.A.; writing-original draft preparation, Z.U.; writing-review and editing, Z.U. and A.H.; visualization, S.S.; supervision, Z.U.; project administration, S.S.; funding acquisition, Z.U. and A.H. All authors have read and agreed to the published version of the manuscript.This research was supported by the Princess Nourah bint Abdulrahman University Researchers Supporting Project (PNURSP2023R259), Princess Nourah bint Abdulrahman University, Riyadh, Saudi Arabia.This work is supported by The University of Lahore under HEC policy.This research was supported by the Princess Nourah bint Abdulrahman University Researchers Supporting Project (PNURSP2023R259), Princess Nourah bint Abdulrahman University, Riyadh, Saudi Arabia.The authors declare no conflict of interest.
Funding:Data Availability Statement: Data sharing is not applicable to this article.Acknowledgments: