We study the interaction between heterogeneously charged surfaces in an electrolyte solution by employing classical Density Functional Theory (cDFT) and Monte Carlo simulations. We observe a consistent behavior between cDFT and Monte Carlo simulations regarding force curves and two-dimensional density profiles. Armed with the validated cDFT, we explore the system's behavior under parameters challenging to simulate directly .Our findings include impacts of domain size, domain charge, domain charge configuration, and bulk electrolyte concentration on the osmotic pressure. Remarkably, the force curve is more sensitive to the domain size for asymmetric configuration than symmetry configuration; the bulk concentration weakly influences the force curve independent of the system configurations.