Adhesion of cell membranes involves multi-scale phenomena, ranging from specific molecular binding at Angstrom scale all the way up to membrane deformations and phase separation at micrometer scale. Consequently, theory and simulations of cell membrane adhesion require multi-scale modeling and suitable approximations that capture the essential physics of these phenomena. Here, we present a mesoscale model for membrane adhesion which we have employed in a series of our recent studies. This model quantifies, in particular, how nanoscale lipid clusters physically affect and respond to the intercellular receptor-ligand binding that mediates membrane adhesion. The goal of this Chapter is to present all details and subtleties of the mean-field theory and Monte Carlo simulations of this mesoscale model, which can be used to further explore physical phenomena related to cell membrane adhesion.