The steady two-dimensional laminar boundary layer that develops on a wedge or in flows where the surface is not aligned with the flow direction is described by the Falkner-Skan boundary layer, named after Falkner and Skan. It also serves as a classic example of flow over a flat plate subjected to a pressure gradient along its length, a situation commonly encountered in wind tunnel experiments. The Falkner-Skan [1] extends the Blasius boundary layer model, which applies to flat plates with no pressure gradient, by incorporating pressure gradient effects.

For the boundary layer flow of a homogeneous incompressible fluid of second grade past a wedge positioned symmetrically with respect to the flow direction, non-similar solutions are developed. It is addressed how the skin friction varies in relation to non-Newtonian factors. Rajagopal et al. [2] . The angle of the wedge is taken as απ. Due to the case of unsteady flow, the velocity potential of uniform flow is considered as $U\left(x,t\right)$. In the case of flow over a wedge, the velocity potential is assumed to be proportional to a power of distance along the wall. In this analysis, a constant $\tilde{C}$ is introduced that depends on the boundary layer thickness, and it has been treated as a dimensionless measure of the unsteadiness parameter. Regarding the two-dimensional flow, the coordinate system is considered as two-dimensional. Because of the no-slip boundary condition, the results of this solution demonstrate lower skin friction or shear stress, boundary-layer thickness, velocity thickness, and momentum thickness.

Non-Newtonian fluids can transport mass, heat, and momentum; they have grown in significance in industrial and engineering applications. Non-Newtonian fluids are defined as liquids without Newtonian viscosity. There are technical and industrial uses for non-Newtonian fluids. There are some liquids in nature without Newtonian viscosities. Applesauce, polymer mixes, colloidal and fermented mixtures, mud, edibles, paper paste, slush, paints, lubricants, clay glazes, slush, shampoos, and so on are examples of these fluids. It is impossible to forecast the rheological characteristics of non-Newtonian liquids using just shear and stress rates. As a result, the three main types of non-Newtonian liquids are differential, integral, and rate. Numerous investigations of differential fluid flows in various geometries have been conducted. Shear thinning, shear thickening, and a normal stress are characteristics of these liquids; however, the effects of relaxation and retardation durations are not well represented by the differential fluid flow model.

Because this topic has several applications in engineering processes, there has been a significant surge in the last few decades in the research of nano fluid flow across an elongating sheet. A fundamental process with a plethora of applications is heat transfer. The rate of heat transfer, which is used to calculate the amount of heat required by machinery and processes, is dependent on thermal conduction in running liquids. Liquids with nanoparticles added to them significantly enhance the properties of the basic fluids, resulting in nano-fluids. The unsteady flow and heat transfer of a dusty fluid have a wide range of applications in air conditioning, refrigeration, chemical processing, pumps, and nuclear reactors. Datta et al. [3] obtained the solution of unsteady heat transfer to pulsatile flow of a dusty viscous incompressible fluid in a channel. In a boundary-layer study, the coupled convection along a vertical non-isothermal wedge buried in a porous material that was saturated with fluid was studied by Kumari & Gorla [4] . In their study, Hossain et al. [5] investigated the forced flow of a viscous incompressible fluid in two dimensions via a horizontal wedge with consistent flow of surface heat, and also represented that the laminar two-dimensional unsteady mixed-convection boundary-layer flow of a viscous incompressible fluid past a sharp wedge has been studied by Hossain et al. [6] . Attia [7] investigated an unsteady MHD Coutte flow and heat transfer of dusty fluid. Rajput et al. [8] discussed the unsteady nonlinear mixed convective flow of nanofluid over a wedge. Watanabe [9] analyzed the behavior of heat transfer under forced environments where the convection flow is stimulated by a wedge towards suction and injection. This was later studied by Ishak et al. [10] , who evaluated the two-dimensional boundary layer flow of viscous fluid induced by a moving wedge. Mahanthesh et al. [11] investigated the influence of nonlinear convective transport on non-Newtonian Dusty fluid flow through a stretched sheet. The impact of thermal radiation of squeezed dusty fluid flow with heat transfer between parallel plates was studied by Abbas [12] . Chandrawat et al. [13] studied numerically the unsteady flow of two immiscible micropolar and dusty fluids through a horizontal plate. In all these investigations, the paper’s goal is to investigate an unstable laminar flow of an incompressible conducting viscous dusty fluid between two parallel plates that are endlessly non-conducting, with the porous medium enclosing the upper plate by Parul Szxena et al. [14] . Islam et al. [15] discussed about on the unsteady laminar flow of heat transferable dusty fluid between two parallel Riga plates. Attia et al. [16] prescribed the finite difference approach to solve the coupled unsteady power-law conducting fluid flow and the continuous dusty viscous fluid flow under the influence of a magnetic field. Last few years, some of the authors [17] - [19] have analyzed the phenomenon of dust particles of Newtonian and non-Newtonian fluids with or without heat transfer between parallel plates.

From the above discussion, it follows that no author has previously given any clear idea about the unsteady dusty fluid flow over a horizontal wedge. Our main investigation is the unsteady dusty fluid flow over a horizontal wedge. The behavior of the flow properties is discussed and presented graphically. Here we want to focus on the study on Numerical Study on the Dusty Fluid along the horizontal System, which is also a significant use in the various practical field. The flow characteristics behavior has been illustrated visually, and a brief discussion of their important applications in several real-world fields follows.

Consider an unsteady, viscous incompressible dusty fluid that flows along a semi-infinite horizontal wedge at $y=0$ moving with a constant velocity $U\left(x,t\right)$. Suppose the direction of flow be taken along the x-axis. In the case of flow over a horizontal wedge, the velocity potential is assumed to be proportional to a power of distance along the wall. Under this consideration, the potential flow velocity of the wedge can be written as

$U\left(x,t\right)=\frac{\upsilon x^m}{{\delta }^{m+1}}$,

where m is a constant, $\delta$ is the time-dependent length scale, which is taken to be $\delta =\delta \left(t\right)$. In Figure 1 , the total angle of the wedge is defined by $\Omega =\alpha \pi$ and

$\alpha =\frac{2m}{m+1}$ is the wedge angle parameter. The velocity $u=U\left(x,t\right)$ is considered at the outer edge of the boundary layer, and then outside of the boundary layer:

$-\frac{1}{\rho }\frac{\partial p}{\partial x}=\frac{\partial U}{\partial t}+U\frac{\partial U}{\partial x}$

which gives that,

$-\frac{1}{\rho }\frac{\partial p}{\partial x}=\frac{{\upsilon }^2}{{\delta }^{2m+2}}\left[m-\tilde{C}\left(m+1\right)\right]x^{2m-1}$

where $\tilde{C}=\frac{{\delta }^m}{\upsilon x^{m-1}}\frac{\partial \delta }{\partial t}$ is taken to be a constant and thus it can be treated as a dimensionless measure of the unsteadiness parameter. The physical configuration of the model is shown in Figure 1 .

The equations relevant to the unsteady two-dimensional problem are governed by the following system of non-linear partial differential equations under the competent boundary layer approximation are given as follows:

Continuity equation:

$\frac{\partial u}{\partial x}+\frac{\partial v}{\partial y}=0$ (1)

$\frac{\partial u_p}{\partial x}+\frac{\partial v_p}{\partial y}=0$ (2)

Momentum equation:

$\frac{\partial u}{\partial t}+u\frac{\partial u}{\partial x}+v\frac{\partial u}{\partial y}=\upsilon \frac{{\partial }^{2}u}{\partial y^2}-\frac{1}{\rho }KN\left(u-u_p\right)+\frac{{\upsilon }^2}{{\delta }^{2m+2}}\left[m-\tilde{C}\left(m+1\right)\right]x^{2m-1}$ (3)

$m_p\left(\frac{\partial u_p}{\partial t}+u_p\frac{\partial u_p}{\partial x}+v_p\frac{\partial u_p}{\partial y}\right)=KN\left(u-u_p\right)$ (4)

It provides the boundary condition for the horizontal wedges as follows:

$\begin{array}{l}u=u_p=0\text{at}y=0\\ u=u_p=U\left(x,t\right)\text{at}y\to \infty \end{array}\}$ (5)

The solution procedure is the systematic approach employed to solve a problem or address a challenge. The solution procedure serves as a systematic framework that guides the computation from the initial formulation to the final solution. From the concept of the above discussion, for explicit finite difference method has been used to solve Equations (6)-(9) subject to the boundary conditions given by (10). A trial-and-error approach has been used to determine suitable step sizes and relaxation parameters in order to ensure numerical stability and convergence of the explicit finite difference solution. To obtain the difference equations, the region of the flow is divided into a grid or mesh of lines parallel to $x$ and $y$ axes where $x$ axis is along the plate and $y$-axis is normal to the plate. Here, we consider that the plate of length $x_{\max }$ i.e. $x$ varies from 0 to 40 and regard $y_{\max }$ as corresponding to $y\to \infty$ i.e. $y$ varies from 0 to 8. The grid pair has been taken (m, n) = (100, 80).

$U{p}_{i,j}^{k+1}=U{p}_{i,j}^k+\frac{\tau }{G}\left(U_{i,j}^k-U{p}_{i,j}^k\right)-\frac{\tau }{h}U{p}_{i,j}^k\left(U{p}_{i+1,j}^k-U{p}_{i,j}^k\right)-\frac{\tau }{l}{V}_{i,j}^k\left(U{p}_{i,j}^k-U{p}_{i,j-1}^k\right)$

The subscripts $i$ and $j$ designate the grid points with $x$ and $y$ co-ordinates respectively and the superscript prepresents a value of time, $t=p\Delta \tau$ where $p=0,1,2,3,\cdots$ from the initial condition.

In this research, the governing equations are organized using standard transformation technology. The explicit finite difference method is used in the transformed governing equations. The MATLAB programming language is utilized to construct appropriate software for resolving these equations. The governing equations are then numerically solved. The necessary values for velocity and concentration profile in various grapes are displayed with varying values of relevant parameters.

The physical properties are graphically discussed for four different values of parameters, namely Wedge angle (α), particle mass parameter (G), fluid concentration parameter (R)and Reynolds number (R_e). Based on the graphical representation of results and discussion, some important findings are mentioned as follows:

1) The velocity for clear fluid and dust phase increases with the increase of wedge angle α, while it decreases with the increase of G, R, R_e.

2) The local shear stress for clear fluid and dust phase decreases with the increase of wedge angle α, while it increases with the increase of G, R, R_e.

3) The average shear stress for clear fluid and dust phase decreases with the increase of wedge angle α, while it increases with the increase of G, R, R_e.