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 Al2O3-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 TC and TH, respectively. The top wall has a steady speed of u=u0. Pertinent non-dimensional parameters such as Reynolds number (Re=10,100,200,400), Hartmann number (Ha=0,10,20), Bingham number (Bn=0,2,5,10,50,100,200), nanoparticle volume fraction (ϕ=0,0.02,0.04), and Prandtl number (Pr=6.2) have been simulated numerically. The Richardson number Ri is calculated by combining the values of Re with a fixed value of Gr, 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 (Nu¯), total entropy generation (Es)t, and Bejan number (Beavg). 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 (Nu¯) drop when Ha and Bn are increased and grow when Re and ϕ are augmented. It is found that (Es)t is reduced by raising Ha, but (Es)t rises with the augmentation of ϕ and Re. It is also found that the ϕ and Re numbers have a positive sensitivity to the Nu¯, while the sensitivity of the Ha and Bn numbers is negative.