The numerical investigation of magneto-hydrodynamic (MHD) mixed convection flow and entropy formation of non-Newtonian Bingham fluid in a lid-driven wavy square cavity filled with nanofluid was investigated by the finite volume method (FVM). The numerical data-based temperature and nanoparticle size-dependent correlations for the AlO-water nanofluids are used here. The physical model is a two-dimensional wavy square cavity with thermally adiabatic horizontal boundaries, while the right and left vertical walls maintain a temperature of and , respectively. The top wall has a steady speed of . Pertinent non-dimensional parameters such as Reynolds number (), Hartmann number (), Bingham number (), nanoparticle volume fraction (), and Prandtl number () have been simulated numerically. The Richardson number is calculated by combining the values of with a fixed value of , which is the governing factor for the mixed convective flow. Using the Response Surface Methodology (RSM) method, the correlation equations are obtained using the input parameters for the average Nusselt number (), total entropy generation (, and Bejan number (). The interactive effects of the pertinent parameters on the heat transfer rate are presented by plotting the response surfaces and the contours obtained from the RSM. The sensitivity of the output response to the input parameters is also tested. According to the findings, the mean Nusselt numbers () drop when and are increased and grow when and are augmented. It is found that is reduced by raising , but rises with the augmentation of and . It is also found that the and numbers have a positive sensitivity to the , while the sensitivity of the and numbers is negative.