Lie Group Analysis of Unsteady Flow of Kerosene / Cobalt Ferroﬂuid Past A Radiated Stretching Surface with Navier Slip and Convective Heating

: In this work, we identiﬁed the characteristics of unsteady magnetohydrodynamic (MHD) ﬂow of ferroﬂuid past a radiated stretching surface. Cobalt–kerosene ferroﬂuid is considered and the impacts of Navier slip and convective heating are additionally considered. The mathematical model which describes the problem was built from some partial di ﬀ erential equations and then converted to self-similar equations with the assistance of the Lie group method; after that, the mathematical model was solved numerically with the aid of Runge–Kutta–Fehlberg method. Graphical representations were used to exemplify the impact of inﬂuential parameters on dimensionless velocity and temperature proﬁles; the obtained results for the skin friction coe ﬃ cient and Nusselt number were also examined graphically. It was demonstrated that the magnetic ﬁeld, Navier slip, and solid volume fraction of ferroparticles tended to reduce the dimensionless velocity, while the radiation parameter and Biot number had no e ﬀ ects on the dimensionless velocity. Moreover, the magnetic ﬁeld and solid volume fraction increase skin friction whereas Navier slip reduces the skin friction. Furthermore, the Navier slip and magnetic ﬁeld reduce the Nusselt number, whereas solid volume fraction of ferroparticles, convective heating, and radiation parameters help in increasing the Nusselt number.


Introduction
Flow and convective heat transfer through a stretching surface play an essential role in research due to their presence in many engineering and industrial applications. Many authors have emphasized this and the details are found in [1][2][3]. To overcome the poor thermal conductivity and increase the other thermophysical properties of the conventional fluids, nanoparticles were suspended in a base fluid. These nanoparticles are called nanofluids and can be generated from diverse operations or chemical deposition mechanisms. Enchantment in the surface area and the rate of heat transfer occurred and many improvements have recently been performed for this issue [4]. This scheme of nanofluid is processed by integrating the pure fluid and classical equations of mass. Many investigations ofnanofluid flow can be found in [5,6].

Problem Formulation
In the current research, it is assumed that a 2D unsteady magneto-forced convective flow of ferrofluid past a radiate stretchable surface with impacts on Navier slip and convective heating are additionally considered. In this work, Cobalt is considered and is treated as a base nanoparticle, with kerosene as a base ferrofluid. The stretchable surface switches on from a fine slot, which is positioned at the starting point of a 2D coordinate system (x, y). At this point, the x-axis is considered all along the stretching direction of the sheet, having stretched velocity U w = ax, which is applied vertically to the sheet externally. A constant magnetic strength B 0 is applied normal to the sheet. The mathematical model describing the system is (see Chamkha [29]) Subjected to the corresponding boundary conditions (see [30][31][32][33][34]): where t, u and v are the time and velocity components along the x and y axes and T is the temperature in the fluid phase. ρ f f stands for the density. µ f f stands for viscosity. β f f Stands for the ferrofluid volumetric thermal expansion coefficient. σ ff stands for electrical conductivity. α f f = k f f /(ρC p ) f f stands for the thermal diffusivity of the ferrofluid. L stands for the slip coefficient, which represents Navier slip, and h f stands for the heat transfer coefficient. T f stands for the uniform temperature of the stretchable surface. k ff stands for the thermal conductivity of ferrofluid. (ρC p ) f f stands for the specific heat of the ferrofluid at a constant pressure. The radiative heat flux qr is approached according to the Rosseland approximation (see [35,36]): where β R and σ 1 stand for the mean absorption coefficient and the Stefan-Boltzmann constant. As carried out by Raptis [35], the fluid-phase temperature variations within the flow are approached to be adequately tiny so that T 4 may be obvious as a linear function of temperature. This is created by extending T 4 in a Taylor series on the free-stream temperature T ∞ and removing higher-order terms to yield By applying Equations (5) and (6) in the last term of Equation (3), we obtain In the current investigation, the following thermophysical relations are applied [37]; Here, χ is nanoparticle volume fraction. Table 1 represents the thermophysical properties of ferrofluid. Table 1. Thermophysical properties of kerosene, water and cobalt [37].

Property
Kerosene Water Cobalt In this stage, the expressions for u, v, and θ will be defined as: Substituting Equations (7)-(9) into Equations (1)-(4), we obtain where

Lie Group Framework
Obtaining the solutions of the PDEs (partial differential equations) (10)-(12) governing the investigation understudy is equivalent to satisfying the constant solutions of these equations under a special continuous one-parameter group. The proposed technique is to search for a transformation group from the primary collection of one parameter scaling transformation. The facilitated form of Lie group framework, namely, the scaling group of transformations ∆ (see [38][39][40][41][42][43][44][45][46]), will be presented here: where κ 1 , κ 2 , κ 3 , κ 4 , and κ 5 are transformation parameters and ε is a small parameter whose interrelationship will be determined by our investigation. Equation (13) may be scrutinized as a point transformation, which transfers the coordinates (x, y, t, Ψ, θ) to ( x , y , t , Ψ, θ). Substituting transformations Equation (13) in Equations (10)-(12), we obtain; The following relations should be determined to reserve the system to be constant: These relations give and the one-parameter group of transformations can be obtained as Developing by Taylor's technique in powers of ε, we obtain: where Γ 1 , Γ 2 , and Γ 3 are arbitrary functions which should be determined by its equations. η and τ are the similarity variable and dimensionless time.
To avert the fluid properties manifesting explicitly in the coefficients of the above equations, determining mass balance in Equation (1), with keeping generality, we have dropped three different convenient arbitrary constants based on the transformations performed previously by Nabwey [25] and Chamkha [29] as follows: As a consequence, we find with the assistance of these formulations in Equation (22). Equations (10)-(12) are characterized as subject to the following boundary conditions: where Ha = B 0 Rd = 4σ 1 T 3 ∞ /k f β R stands for the radiation parameter. δ = Lµ f /2 υ f /a stands for the velocity slip parameter.
The local skin-drag coefficient and local Nusselt number can be written respectively, as

Numerical Method
Following [47,48], Equations (23) and (24), subject to (25), will be solved using the local similarity method, where the first derivatives with respect to τ are neglected and the Equations (23) and (24) with boundary conditions (25) can be re-written as Ξ 2 Pr The boundary conditions (25) remain the same. These ordinary differential equations with the boundary conditions (25) can be solved numerically by applying the Runge-Kutta-Fehlberg method (RKF7 45). Following [47,48], for the local non-similarity solution, now we hold all the terms by assuming the new auxiliary functions F(τ, η), and Θ(τ, η), which are defined by Thus, Equations (23) and (24) can be expressed as The ODEs (31)-(33) subject to (25) and (34) were solved numerically by employing the Runge-Kutta-Fehlberg technique (RKF45) using MAPLE-19 software (MAPLE 2019.0, Maplesoft, Waterloo, ON, Canada). This method is generally known as one of the most excellent methods available for obtaining the solutions of nonlinear differential equations and provides more accurate results. The step size was selected. For the similarity variable η max , Equations (25) and (34), were replaced as The selection of η max = 5 guarantees that all numerical solutions approached the asymptotic values properly.

Results and Discussions
In this study, we investigated the unsteady magento-flow and heat transfer of Cobalt-kerosene ferrofluid past a stretchable surface. The influences of several key parameters on the dimensionless velocity f (τ, η), temperature θ(τ, η), skin friction C f (τ, 0), and Nusselt number Nu(τ, 0) are examined. The Lie group method is employed to reduce partial differential equations and local similar and non-similar models are solved employing the RK-45 technique.
The effects of the magnetic field Ha and dimensionless time τ on the velocity are symbolized in Figure 1a and on the dimensionless temperature in Figure 1b, respectively. As the time increases, the velocity at the surface rises. The magnetic field generates Lorentz strength on the fluid particles, which resist the fluid and reduce the fluid velocity, as shown in Figure 1a. Consequently, the velocity boundary layer thickness decreases. Due to the decline in velocity, the temperature increases. In the thermal boundary layer, the temperature declines to the ambient temperature. The thermal boundary layer thickness reduces with an enlargement of the dimensionless time, as exhibited in Figure 1b. The influence of the solid volume fraction of nanoparticles χ and Navierslip δ on the velocity and temperature is depicted in Figure 2a,b when τ = 0.5. In the absence of slip, the velocity is found to be higher for the pure regular fluid. At the surface, the velocity decreases with the increase in the slip and solid volume fraction, as shown in Figure 2a. No appreciable impact of χ could be observed at the surface as well as within the velocity boundary layer. The velocity boundary layer thickness enlarges with δ, which enlarges the thermal resistance and reduces the heat transfer rate; see Figure 2b. The variation of the dimensionless temperature with the solid volume fraction χ is depicted in Figure 2b. In the absence of slip, the temperature is lower at the wall and intensifies with δ. As expected, the temperature at the wall is higher for the regular fluid and dwindles with an intensify in the solid volume fraction χ. This is due to the higher thermal conductivity of Cobalt nanoparticles. With the addition of nanoparticles, the thermal conductivity of the ferrofluid increases and the heat transfer rate is enhanced.   Figure 3a,b presents the effects of radiation parameter Rd and Biot number Bi on the velocity and temperature curves. It is important to note that equations of momentum and energy are independent of each other. The momentum equation and the velocity boundary conditions are independent of the radiation parameter Rd and convective heating parameter Bi. Therefore, there is no influence of these parameters on the velocity, which is obvious from Figure 3a. On the other side, the surface temperature increases significantly with a strengthen in both Rd and Bi. As a result, the thermal boundary layer thickness is boosted, with an increase in both parameters, as depicted in Figure 3b. The radiation parameter Rd reveals an enhancement in radiative heat, which improves the thermal state of fluid, causing its surface temperature to increase. Similarly, as the convective heating parameter increases and tends to infinity, the convective boundary condition changes to an isothermal boundary condition. The variations in skin friction and Nusselt number with the magnetic field Ha are depicted in Figures 4 and 5 for different values of the velocity slip δ and the solid volume fraction χ at τ = 1 and τ = 2, respectively. In the presence of magnetic strength, a Lorentz force is generated which resists the fluid and reduces the velocity curve. Therefore, the skin friction enhances with Ha, as shown in Figures  4 and 5a at different dimensionless times. As expected, the skin friction increases with dimensionless time τ. In the absence of velocity slip δ, the velocity curves are higher at the surface and decline with an increment in slip parameter δ. Consequently, the skin friction declines with the boosting of slip parameter δ. For the pure regular fluid, the skin friction is lower and increases with a rise in the solid volume fraction χ. This is due to an evolution in the ferrofluid density with the increased volume fraction of cobalt nanoparticles. Figures 4 and 5b illustrate the variation of Nusselt number with the magnetic field Ha and the volume fraction of ferroparticles χ at different dimensionless times. Like skin friction, Nusselt number also increases with dimensionless time. Due to Lorentz force, the dimensionless velocity decreases and, as a result, the Nusselt number is reduced with an increasing magnetic field. Similarly, the velocity decreases due to an intension in the slip and the Nusselt number reduces. The thermal conductivity of ferroparticles increases with an increase in the volume fraction of ferroparticles. Consequently, the Nusselt number increases with increasing χ.   Figure 6a, b presents the comparison of Nusselt numbers for kerosene oil and water for the same parameters. Due to the smaller Prandtl number Pr for water, the Nusselt numbers are found to be lower than kerosene. The Prandtl number Pr compares the rate of thermal diffusion in comparison to the rate of momentum diffusion. The higher the Prandtl number Pr, the higher the Nusselt number will be. It is also noticed that an increase in the Biot number Bi and radiation parameter Rd leads to an increase in the Nusselt number. These Nusselt numbers also become greater with increasing dimensionless time.

Conclusions
In this study, application of the scaling group of transformations to the unsteady magneto-flow of ferrofluid past a stretching surface was employed. The impacts of Navier slip, radiation and solid volume fraction of ferroparticles, as well as convective heating, were also investigated. From this study, it was concluded that: • Employing the Lie group framework, the symmetries of the partial differential equations are presented exclusively in this investigation, these equations are reduced to self-similar equations utilizing translational and scaling symmetries. Numerical solutions for scaling symmetry are obtained applying the Runge-Kutta-Fehlberg method.

•
The magnetic field, Navier slip, and solid volume fraction of ferroparticles tend to reduce the dimensionless velocity.

•
The radiation parameter and Biot number have no effects on the dimensionless velocity.

•
Magnetic field, radiation, Biot number, and Navier slip increase the surface temperature, whereas the solid volume fraction of ferroparticles reduces the surface temperature.

•
The magnetic field, dimensionless time and solid volume fraction increase skin friction, whereas Navier slip reduces the skin friction.

•
The magnetic field and Navier slip reduce the Nusselt number, whereas solid volume fraction of ferroparticles, convective heating, and radiation parameters help in increasing the Nusselt number.

•
The Nusselt number for kerosene oil is higher than for water.