This work aims to offer a mathematical model for two-phase flow that investigates the interaction of Casson nanofluid and dust particles across a stretching surface. MHD Darcy–Forchheimer porous medium and Fourier’s law through Cattaneo–Christove thermal flux are also considered. The governing equations for the two phases model are partial differential equations later transmuted into ordinary ones via similarity transforms. The Runge–Kutta method with the shooting tool is utilized numerically to solve the boundary layer equations computed in MATLAB to obtain numerical results for various pertinent parameters. The numerical outcomes of momentum, temperature, and concentration distribution are visible for both phases. The results of the skin friction, heat transfer coefficients, and the Sherwood number are also visible in the graphs. Furthermore, by comparing the current findings to the existing literature, the validity of the results is confirmed and found to be in good agreement. The fluid velocity is reduced against increasing strength of Casson fluid parameter, enhanced the fluid phase and dust phase fluid temperature. The temperature declines against the growing values of the relaxation time parameter in both phases. Dusty fluids are used in various engineering and manufacturing sectors, including petroleum transportation, car smoke emissions, power plant pipes, and caustic granules in mining.