Fractional Model and Exact Solutions of Convection Flow of an Incompressible Viscous Fluid under the Newtonian Heating and Mass Diffusion

In this paper, we consider a natural convection flow of an incompressible viscous fluid subject to Newtonian heating and constant mass diffusion. %e proposed model has been described by the Caputo fractional operator. %e used derivative is compatible with physical initial and boundaries conditions. %e exact analytical solutions of the proposed model have been provided using the Laplace transform method. %e obtained solutions are expressed using some special functions as the Gaussian error function, Mittag–Leffler function, Wright function, and G-function. %e influences of the order of the fractional operator, parameters used in modeling the considered fluid, Nusselt number, and Sherwood number have been analyzed and discussed. %e physical interpretations of the influences of the parameters of our fluid model have been presented and analyzed as well. We use the graphical representations of the exact solutions of the model to support the findings of the paper.


Introduction
e main challenge during these decades is modeling realworld phenomena using fractional operators. is challenge is motivated by the memory effect observed in fractional operators and the generalization of the differential equations to noninteger order derivatives. Many explanations and tentatives to explain the utilities of fractional calculus have been addressed in the literature. e readers are advised to take a look at the following application papers in biology [1,2], science and engineering [3,4], physics [5][6][7], mathematical physics [8][9][10], and many others domains [11,12]. Fractional calculus is also a field with many controversies and discussions on the validity of fractional operators. Note that there exist many fractional operators in the literature: fractional operators with singularities as the Riemann-Liouville derivative and the Caputo fractional derivative [13,14]. ere also exist many generalizations of the previous operators; see [15]. Fractional operators without singularities are known as the fractional operator with exponential kernel proposed by Caputo-Fabrizio and the fractional operator with Mittag-Leffler kernel provided by Atangana and Baleanu [16][17][18]. All cited operators have their advantages and their conveniences reported in the literature [16,19]. e present paper focuses on applying the Caputo derivative to fluids models. e first reason for using this fractional operator is the memory effect, and the second reason is that this operator is more realistic because the derivative of a constant is zero, contrary to the Riemann-Liouville operator in which the derivative of a constant is not zero.
In the literature, there exist many papers in the same directions. We recall some of them, which will permit us to see the findings addressed in the literature and the advantages of our present paper regarding the existing literature. Khan et al. [19] model heat and mass transfer of secondgrade fluids over a vertical plate using derivative with exponential kernel and derivative with Mittag-Leffler kernel.
ey have used the Laplace transforms to get the analytical solutions of the proposed model. Samiulhaq et al. [20] studied the magnetic field influence on unsteady free convection flow of a second-grade fluid near an infinite vertical flat plate with ramped wall temperature embedded in a porous medium. e Laplace transform method has been used in the investigations related to the exact solutions of the model considered in the paper. Shah et al. [21] proposed an investigation related to the influence of magnetic field on double convection problem of viscous fluid described by fractional-order derivative over an exponentially moving vertical plate. Zafar et al. [22] proposed the exact solutions of the model described by fractional operator obeying the rotational flow of some fractional Maxwell fluid. Imran et al. [3] proposed the study on the semianalytic solutions of a viscous fluid with old and new definitions of fractional derivatives. Hussanan et al. [23] investigated the soret effects on unsteady magnetohydrodynamic mixed convection heatand-mass-transfer flow in a porous medium with Newtonian heating. Reyaz et al. [24] studied modeling with Caputo operator the MHD Casson fluid flow over an oscillating plate with thermal radiation. Narahari [25] presented the effects of thermal radiation and mass diffusion on free convection flow near a vertical plate with Newtonian heating. Sheikh et al. [26] proposed a comparative study between the fractionalorder derivatives without singularities to the convective flow of a generalized Casson fluid. Ali et al. [9] applied the fractional operation with the exponential kernel to MHD free convection flow of generalized Walters' B fluid model. e authors proposed a fractional model and tried to get solutions using the applications of the Laplace transforms. Abro [27] proposed using fractional operator investigation in analytic investigation of thermodiffusion process on free convection flow. Shen et al. [28] used Fourier sine transform to propose the analytical solutions of the Rayleigh-Stokes problem for a heated generalized second-grade fluid described by the fractional derivative. Vieru et al. [29] first proposed an investigation related to the analytical solutions of the fractional model described free convection flow near a vertical plate with Newtonian heating and mass diffusion using Laplace transform. For works related to the unsteady flow of generalized Casson fluid with fractional derivative, the readers are recommended to take a look at [6]. For papers in the same directions of research as addressed in our paper, see in [30][31][32][33].
is paper's question is to determine the exact solutions of the fluid model using the Laplace transforms. e advantages of this paper are as follows: at first, the solutions will be expressed using well-known functions as Gaussian error function, Mittag-Leffler function, Wright function, and G-function [29]. us, the graphical representations of the dynamics of the exact solutions of the proposed model will be easy to be obtained via MATLAB manipulations. Another advantage of this paper is that some new initial conditions not previously considered in the literature will be used and studied in detail. e influence of the order of the fractional operator will be studied and analyzed. e physical interpretations of the influence of the Caputo derivative will be proposed. e same procedure will be done with the parameters of the model. In other words, the influences of the model parameters in the dynamics will be focused on, analyzed, and interpreted. In our modeling, the influence of the Nusselt number and the Sherwood number will be presented and discussed, too. e fact that the solutions can be obtained via the resolution of the secondorder differential equation is also a significant advantage of the Laplace transform method. e present paper is structured as follows: In Section 2, the fractional operators and the special functions used in this paper have been recalled. In Section 3, the natural convection flow of an incompressible viscous fluid subject to Newtonian heating and constant mass diffusion model has been proposed using Caputo derivative. In the modeling, velocity, temperature, and concentrations of the fluid have been considered and analyzed for the considered model. In Section 4, the solutions procedure has been proposed via the Laplace transform. In Section 5, we discuss the main findings of the paper via graphical representations. In Section 6, we give the conclusion and the final remarks and discuss the main findings of the present paper and its advantages in the existing literature.

Preliminaries
is part is devoted to recalling certain fractional operators necessary for our present investigations. As known, there exist many fractional operators such as Riemann-Liouville integral, the Caputo fractional derivative, Riemann-Liouville derivative, derivative with the exponential kernel, derivative with Mittag-Leffler kernel. e use of fractional operators in modeling real-world problems is motivated by the memory effect. Here, we will investigate using the Caputo operator because it satisfies, in particular, the initial conditions which we planned to use. Note that the Riemann-Liouville derivative does not satisfy the initial conditions used in the present investigations. We give the following definitions.
We represent the Riemann-Liouville integral as the following form [13,14]: where the function u: [0, +∞[ ⟶ R and the order obeys the condition α > 0. e function Γ(. . .) denotes the Gamma Euler function. e Caputo fractional derivative as reported in the literature can be represented as the following form [12,15]: where the function u: [0, +∞[ ⟶ R and the order is taken in the interval α ∈ (0.1).
Without losing generalities, we recall the Riemann-Liouville derivative as reported in the literature with the following form [12,15]: with the condition that t > 0, the order of the fractional operator satisfies the condition α ∈ (0, 1), and Γ(. . .) represents the Gamma Euler function.
ere exists a recent fractional operator in fractional calculus. We mean the following definition of the Caputo-Fabrizio fractional derivative [18] of the function u: [0, +∞[ ⟶ R, of order α in the following term: with the following condition t > 0, the order of the fractional operator satisfies α ∈ (0, 1) and CF(.) is the normalization term and respects the condition CF(0) � CF(1) � 1. e integral associated with the Caputo-Fabrizio derivative is represented as the following form: We have the fractional derivative with Mittag-Leffler kernel and its associated integral, seen in the literature particularly in the following paper [17]. e following is the definition of the Atangana-Baleanu derivative of the function u: [0, +∞[ ⟶ R, of order α: with respect to the condition t > 0, the order of the fractional operator α ∈ (0, 1), and AB(.) denotes the normalization term and satisfies the condition AB(0) � AB(1) � 1. e integral associated with the Atangana-Baleanu fractional derivative is represented as follows: We finish this section by recalling the Laplace transform because this transformation is the key of the present paper. In the following line, we recall the Laplace transform of the Caputo derivative, for the function u, and we have the following: where the order α satisfies the relationship α ∈ (0, 1). e problem consists of establishing the analytical solutions of the fluid models that need some special functions. We mean the Mittag-Leffler function and the Wright function. ey will play an important role in expressing the form of analytical solutions. For the Mittag-Leffler function with three parameters [29], we have the following definition: where α > 0, β ∈ R, and x ∈ C. Note that when α � β � 1, we have the classical exponential function. We define the Wright function [29] with three parameters as follows: where the conditions σ ∈ (0, 1), β ∈ R, and x ∈ C are satisfied.

Modeling Dynamics of the Fractional Model
In this part, we propose the constructive equations of the model subject of investigation in this paper. We consider an unsteady free conventional flow of an incompressible viscous fluid over an infinite vertical plate with mass diffusion and chemical reaction [3]. We let the x-axis be obtained along with the plate in the vertical direction, and the y-axis is taken normal to the plate as presented in the literature [2,29]. In the present modeling, the use of the Caputo derivative in the newly added assumption. At the neighborhood of the initial time of the diffusion processes, the temperature is constant at T ∞ and the concentration on the plate is fixed to C ∞ while plate and fluid are at the rest. We suppose that at t � 0 + , the heat and the mass from the plate to the fluid is raised to a variable temperature and the concentration level near the plate is considered as . For the rest of the modeling, we use the Boussinesq rules. Furthermore, we consider the Caputo fractional operator. Fick's first and second laws in the context of fractional operators reported in the literature [21] will be applied to get the constructive equations. Note that inspired to the modeling described in the paper [21], the Linear momentum equation with Caputo operator can conventionally be represented by the following differential equation: where ϵ represents the shear stress. We use the Caputo derivative to express the shear stress as the following form: e next step is to introduce equations (12) into (11). We obtain the following fractional differential equation: Applying the Riemann-Liouville integral I 1− α to both sides of equation (13), we arrive at our first fractional differential equation described as follows: e Fick second law related to the thermal balance equation as reported in the literature [21] can be represented by the following form: where q is the thermal flux and obeys Fick's first law reported with the fractional operator and represented as follows: We replace equations (16) into (15); then, the temperature distribution considered in the present paper can be represented as follows: We repeat the previous procedure by applying the Caputo operator I 1− α to both sides of equation (17). We arrive at our second fractional differential equation satisfied by the temperature distribution and described as follows: We now give the constructive equation in which the concentration obeys in our modeling. e diffusion equation can be expressed using the Riemann-Liouville integral as follows: where J satisfies the molecular diffusion equation given by Fick's first law represented as follows: Replacing equation (20) into (19), then the concentration distribution in our present modeling satisfies the following equation: erefore, the constructive equations of the present work are assigned in the following equations: e initial and boundaries conditions are summarized in the following lines, and we have Furthermore, we consider that To simplify the previous model for helpful manipulations, we will try to introduce dimensionless variables. Here, in our modeling, we set the dimensionless variables [29] as follows: Note that the rest of the modeling follows [29]. Plugging equations (29) into (22) by removing the ″ * ″ and multiplying by 1/ρ, we have the first constructive equation of our fluid model related to the velocity described as follows: where Gr and Gm are defined in the next lines. Replacing equations (29) into (23) by removing the ″ * ″ and divising by ρC p , we have the second constructive equation of our fluid model related to the temperature described as follows: where Pr is defined in the next lines. Plugging (29) into (24) by removing the ″ * ″ , we have the third constructive equation of our fluid model related to the concentration described as follows: where Sc is defined in the next lines. e velocity, temperature, and concentration according to dimensionless variables at initial time t � 0 are summarized as follows: Note that we consider velocity v, temperature ϕ, and concentration θ satisfying the boundary conditions represented in the following relationships: In our previous constructive equations, we consider that the Prandtl number, thermal Grashof number, mass Grashof Journal of Mathematics number, and chemical reaction parameter are defined, respectively, as follows: To complete the modeling, we give tables of nomenclature where all parameters described in this paper have been summarized for more comprehension (Tables 1 and 2).

Dynamics Approaches
is section provides the analytical solutions of the equations established in our modeling section. e procedure to obtain the solutions starts with the second and the third equation and finishes with the first equation. is procedure is adopted because the first equation of the constructive equation (30) contains solutions of the second and last equations of models (31) and (32). In this section, the technique will use the Laplace transform and the Mittag-Leffler function, Gaussian error function, and Wright function. ese functions will play an important role in expressing the forms of the analytical solutions.

Dynamics Approaches of the Concentration.
We begin with the fractional diffusion equation with the reaction term expressed in equation (32) under the initial and boundaries conditions (33)-(35). See the following fractional differential equation: Its initial and boundary conditions are described in the following equations: As previously mentioned, the considered initial and boundary conditions are one of the novelties of the present paper. e solution approach will apply the Laplace transform to both sides of (38). Note that the Laplace transform of the fractional operator will also be used, and we obtain the following relationships: Utilizing the Laplace transform of the initial condition described by θ(0, s) � 1/s 2 and solving the second-order differential equation (41) with the classical mathematical method, the exact analytical solution is given as follows: We can observe that in (42), we have a product of two Laplace transform functions. Using convolution rule, the analytical solution of (38) follows the form described by the relationship where, in particular, function f is represented as follows: and function g is a function obtained via integration and the Wright function: It is also important to propose the analytical solution when the order of the fractional operator is α � 1. To be more clear, we repeat the previous procedure. It will permit some readers to do the same investigation when the integerorder derivative is used. Note that, in integer-order We utilize the Laplace transform of the initial condition; that is θ(0, s) � 1/s 2 . e second-order differential equation (46) with classical method has the following solution: e inverse of the Laplace transform according to the second variable can be represented by the following equation: We finish this part by focusing on a special case that can be obtained when the α � 1 and K � 0. K � 0 means the chemical reaction parameter is null. e procedure of the solution is similar to the integer-order version. After calculation, we have the inverse as follows: Applying the inverse of the Laplace transform of the previous function (49), we obtain the following relationship:

Dynamics
Approaches of the Temperature. In this part, we try to determine the analytical solution of the fractional diffusion equation presented in (31). To be more clear, we recall the fractional differential (31), the problem is defined as follows: and its initial and boundary conditions are described by the following equations: e procedure of getting the solution does not change. We apply the Laplace transform in equations (51) and (53) and get the following relationships: Note that the Laplace transform of the initial condition described in (53) is given by ϕ(0, s) � 1/s 2 . Combining it with the solution of the second-order differential equation (54), we obtain the following relation: e inverse of this Laplace transform can be written using the wright function described in the preliminary section. We have the following exact analytical solution: We continue with a special case. Let the order of the fractional operator α � 1. e Laplace transform is represented as the following forms: e Laplace transform of the initial condition described in (53) gives ϕ(0, s) � 1/s 2 . Combining it with the solution of 6 Journal of Mathematics the second-order differential equation (57), we obtain the following relation: Applying the inverse of the Laplace transform of the previous function (46), we obtain the following relationship:

Dynamics Approaches of the Velocity.
We propose the velocity dynamics in the model proposed in the modeling section. Here, the procedure does not change, and the Laplace transform will be applied. Let the fractional differential equation be defined by the following equation: and its initial and boundary conditions are described in the following equations: For the simplification in the computations, we suppose that κ � 1 + 1/β. e Laplace transform as described in the previous sections will be applied. e particularity of this present step is that the Laplace transforms in equations (60)-(62) will be used. We have the following relationship: In the resolution of this equation, we need to replace the terms in the second member by their values. Using the function in equations (42) and (55), we have to solve the following second-order differential equation: In the resolution of this equation, we need to get the homogenous solution and the particular solution. We first consider the equation defined by Its solution considering the Laplace transform of the initial condition can be represented as the following form:

v(x, s) � a(x, s) + b(x, s) + c(x, s) + d(x, s),
with the functions defined as c(x, s) � − Gm KκSc Note that the constant b 1 � KκSc/(κSc − 1). e analytical solution of the model considered in our paper will be calculated by the inverse of the Laplace transforms of the functions represented in equations (67) to (70). We begin with the first function. We need to invert the function a using the convolution rule and have the following: where, in particular, the functions f is represented by the following form: We continue with the inverse of the Laplace transform of the function b and adopt the same procedure as previously made: We continue now with the inverse of the function c. Here, we use the procedure adopted in the determination of the concentration. We have the following: We can observe that, in equation (74), according to the product of two Laplace transform functions, the analytical form follows the form described by where, in particular, function f 2 is represented by the following form: Journal of Mathematics We finish the first part of the inversion by the function d. We use the same procedure as in the previous parts and get the following analytical form: where, in particular, function f 2 is represented by the following form: and the function g is a function obtained via integration and the Wright function: We now try to solve the second-order differential equation defined by the following relationship: e solution of equation (80) is not difficult to be established and can be represented as the following form: where the functions m and n are defined as Applying the inverse of the Laplace transform and combining it with the convolution product, the analytical solution of equation (80) can be represented as follows: where, in particular, functions f 1 and f 2 are represented by the following forms: e analytical solution of the fractional differential equation described in equation (60) can be obtained by combining the exact analytical solution in equations (71)-(77) and equation (83). We have the following: We finish this part with a special case. In this part, the special case is when the Caputo derivative is equivalent to the integer-order derivative. So, we consider that α � 1. For more clarity, the process to get the solution is repeated with the new form of the Laplace transform, and we have the following: (88) e resolution of the second-order differential (87) gives the following analytical Laplace solution: where the functions are defined as follows:  κSc). e analytical solution of (60) is obtained via inverting the Laplace transform of equation (89). We begin with the first function. e inverse of the function a is given by the following formula: (96) e inverse of the Laplace transform of the function b is given by the following formula: (97) e inverse of the Laplace transform of the function c is given by the following representation: e inverse of the Laplace transform of the function d is given by the following representation: (99)

Journal of Mathematics
For the inverse of the Laplace transform of the function m and n, we use convolution product properties and get the following representations: .
In the discussion section, these solutions will be implemented via MATLAB, and the solution will be drawn. e advantages of the analytical solutions are useful and not difficult to be implemented via MATLAB. e expression of the Nusselt number is obtained after the inverse of the Laplace transform of a certain function, and we use the temperature distribution given by After calculations, the Nusselt number is represented by the following value: e formula of the Sherwood number is obtained after the inverse of the Laplace transform of a certain function, and we use the concentration distribution given by is number is obtained as previously mentioned using the inverse of the function represented in equation (103) and expressed via a special function G(.) known as the Lorenzo-Hartley function. We have the following value:

Results and Interpretations
is section represents the dynamics using the exact analytical solutions described in the previous section. We will also analyze the impact of the order of the fractional operator. e influence of the mass Grashof number, thermal Grashof number, Prandtl number, Schmidt number will also be discussed and physically interpreted. We will also analyze and discuss the Nusselt number and Sherwood number in the fractional context.

Dynamics of the Concentration.
is section represents the dynamics of the concentration in our model (30)-(32). In Figures 1(a), 1(b), 2(a), and 2(b), we consider Sc � 4 and Sc � 10, K � 0.06. In the first part, we focus on the influence of different orders of the Caputo derivative. We represent the concentration of the fluid in Figures 1(a), 1(b), 2(a), and 2(b) versus the variations of the order of the fractional operator.
We now consider analyzing the influence of the Schmidt number Sc. We work with different values of the order of the Caputo derivative, and the time is fixed to τ � 4 and K � 0.06. e concentrations are represented in Figures 3(a), 3(b), 4(a), and 4(b).
We notice with Figures 3(a), 3(b), 4(a), and 4(b), when the Schmidt number Sc increases the concentration of the considered fluid decreases as well. is fact can be explained physically by the fact that when the Schmidt number and the order of the Caputo derivative increase, they reduce the diffusivity of the fluid. us, the concentration of the fluid decreases.
e same behaviors can be observed with the chemical reaction parameter. erefore, this part is omitted.
We finish the analysis of the influence of the parameters by representing the dynamics of the concentration when the chemical reaction parameter k � 0. If we use the analytical solution presented in equation (50) in the graphical representations, we have Figures 5(a) and 5(b).
In Figures 5(a) and 5(b), we can observe that the influence of the order of the Caputo derivative and Schmidt number does not change.

Dynamics of the Temperature.
e objective of this section is to give graphical representations of the temperature distribution and analyze the influence of the order of the Caputo derivative α and the Prandtl number Pr. Here, we work with equation (56). In this section, we consider different orders of the Caputo derivative and different values of the Prandtl number. We have the following graphical representations (Figures 6(a), 6(b), 7(a), and 7(b)) for the dynamics of our model (30)-(32) considering the temperature.
From Figures 6(a), 6(b), and 7(a), we can observe that as the order of the Caputo derivative increases when the time is fixed to τ � 5, the temperature distribution increases as well. We notice the acceleration effect generated by the Caputo derivative on the temperature dynamics of the considered fluid model. From Figure 7(b), we observe that when the Prandtl number increases, the temperature distribution decreases. is behavior is normal because big values of the Prandtl number reduce the diffusivity of the considered fluid too. us, the temperature decreases with it increase.

Dynamics of the Velocity.
is section is interesting and we will analyze the impact of all the parameters of the considered model on the velocity dynamics. Using the solution described in equation (87) and the resolution procedure adopted to obtain the analytical form of the velocity, we can notice both the concentration and the temperature impact the velocity, which make this section more interesting. e influence of the mass Grashof number Gm, ermal Grashof number Gr, volumetric coefficient of thermal expansion β, Prandtl number Pr, and Schmidt number Sc should be analyzed. In our first case, fix the following values Gr � 5, Pr � 5, κ � 9.5, Gm � 5, K � 1.5, and τ � 2. We have Figures 8(a), 8(b), 9(a), and 9(b) when we work with different orders of the fractional operator.
Firstly, we observe that in Figures 8(a), 8(b), 9(a), and 9(b), when the order of the Caputo operator increases, the velocity of the considered fluid increases as well. is behavior is due to the accumulation of the memory, which causes an acceleration effect in the velocity. e second remark is that the Schmidt number Sc generates a decrease in the velocity when its value increases. e same behaviors are obtained also with the increase of the Prandtl number Pr. ese behaviors are already mentioned with temperature and concentration because these values reduce diffusivity. We focus on the influence of the mass Grashof number Gm.
Before we represent the graphics with the values previously fixed, the difference at this stage is that we take different values for the mass Grashof number Gm. We have velocities. is behavior can be explained physically by the following fact: the mass Grashof number influences the ratio of the species buoyancy force to the viscous hydrodynamic force and its increase generates an increase in the species buoyancy force, which in turn implies an increase in the velocities. e same behavior can be observed with the thermal Grashof number Gm (Figures 12(a), 12(b), 13(a), and 13(b)). Referring to Figures 12(a), 12(b), 13(a), and 13(b), we observe that the velocities increase when the thermal Grashof number Gr increases. ese dynamics are because the thermal Grashof number generates an increase in the buoyancy forces, which in turn causes an increase in the momentum boundary layer thickness and then the velocity of the fluid considered in our model. e comparative studies will be made versus the results in [21,29]. As considered in [29], the same model has been considered. For the temperature distribution, the constructive equations described for the temperature in the present paper and [29] are the same, but the initial and boundary conditions differ. e dynamics observed for the temperature in [29] and the present paper in Section 5.2 are similar (Figure 14(a)). e similarity can be explained by the fact that our present initial and boundary conditions are a particular case of the boundary and initial conditions considered in [29]. e advantage of our initial condition, in this case, is that it is in a polynomial form, and it is more realistic because it excludes negative value for the derivative of the temperature with respect to the state. In Section 5.1, the concentration has been considered and the constructive equation is the same as the constructive equation for the concentration in [29]. However, the initial conditions are totally different. Here, our initial and boundary condition depends on time, but in [29], they are constant; the consequences are that the dynamics in the figures differ; see Figure 14(b). Note that the difference is caused in this case by the difference in the initial condition. e same initial condition for initial and boundaries conditions have been considered for the velocity. We observe the dynamics of the velocities in [29], and findings in Section 5.3 have the same form, but the values are impacted because the concentrations differ in the paper (Figure 15(a)).
A similar comparison can be made with paper [6].

Conclusion
is paper has addressed and focused on the analytical solutions of the fluid model described by the Caputo derivative using the Laplace transform method. e special functions have been used to describe the analytical solutions and their graphical representations. We have analyzed the influence of the order of the Caputo derivative and find that it plays acceleration or retardation effects. ese behaviors are due to the accumulation of the memory effect described in the paper. We have also analyzed the influence of the Prandtl and Schmidt numbers in the dynamics of the concentration and the temperature. We find that these numbers generate a decrease in these distributions. We also find that Grashof numbers generate a decrease in the velocities due to the impact of these numbers on the buoyancy forces. In the future direction of research, the analysis adopted in this paper can be addressed with the same model but with different initial and boundary conditions, which play an important role in the form of the exact analytical solutions. e graphical representations have been used in this paper to support all the paper's findings.

Data Availability
No data were used to support this study.

Conflicts of Interest
e author declares no conflicts of interest.