The present contribution provides a numerical treatment for the flow of unsteady natural convection of non-Newtonian Casson fluid model loaded with dusty particles over a vertical wavy plate. It is assumed that the fluid is incompressible and the dusty particles are considered to be spheres with the same size. Via the help of primitive variable formulation, the dimensionless boundary layer governing equations are reduced into a convenient coordinate system. The transformed resulting system of governing equations is tackled numerically employing the fully implicit finite difference method. The impact of various controlling physical parameters such as the amplitude of the wavy plate, the dimensionless time and fluid-particle interaction parameters on the velocity and temperature of both fluid and particle phase is analyzed. It is found that the magnitude of the velocity of both fluid and particle phase diminishes with increasing amplitude of the wavy plate parameter, while an opposite behavior is observed on temperature distributions as the amplitude of the wavy plate parameter rises. A significant impact is noticed for the heat transfer rate with enhancement values of fluid-particle interaction parameter. The numerical solutions are compared with the available data in the literature. Quantitative comparison illustrates good compatibility between the current and previous results. Keywords Dusty fluid • Wavy plate • Non-Newtonian Casson fluid • Natural convection • Finite difference List of symbolŝ a Amplitude of the wavy surface a Dimensionless amplitude of the wavy surface C Concentration of fluid C fx Local skin-friction coefficient c p Specific heat at constant pressure c s Specific heat of particle phase D Mass diffusivity of fluid D Mass concentration of the dust particle Ec Eckert number