Two Semi-Implicit-Method for Pressure-Linked-Equations (SIMPLE) techniques have been performed to solve the dimensionless equations governing a dusty hybrid nanofluid flow in wavy enclosures contain a volumetric heat source. These techniques are applied to evaluate the pressure terms for both the hybrid nanofluid and the dusty particles based on the control volume solver. Two systems of equations are proposed to simulate the hybrid nanofluid phase and the dusty particles phase. In addition, a body-fitted method is applied to map the irregular domain into a rectangular domain, and an inverse map technique is used to present the obtained data inside the given wavy domain. The hybrid mixture, here, is consisting of water as a base fluid, and the nanoparticles are Al2O3 and Cu. The controlling parameters are the Rayleigh numbers RaE and Ral, the ratio of the densities of the mixture Ds, the dusty parameter αs, and the total nanoparticles volume fraction. The results revealed that the absolute values of the stream function are reduced by 54.5% when the heating modes are switched from (Ral/RaE) < 1 to (Ral /RaE) > 1. Also, the average Nusselt is enhanced by 5.2% at a = 0.9, 6.74% at a = 0.95, and 11.36% at a = 1.1 when the nanoparticles volume fraction is increased from 1% to 5%.