Computer simulations play an important role in designing new polymers as well as in predicting properties of existing polymers. In this paper, the blend compatibility of poly(vinyl alcohol) (PVA) with poly(methyl methacrylate) (PMMA) was studied over the wide range of compositions allowed by the atomistic and mesoscopic simulation methods. The Flory-Huggins interaction parameter, chi, of the blends computed using the atomistic simulation confirmed the blend compatibility for compositions containing > 60 wt % of PVA. This observation was further supported by differential scanning calorimetric experiments. Solubility parameters of the polymers obtained from the simulation procedure were in good agreement with those of the literature data. Simulation results were further supported by the spectral and solution property measurements. From the atomistic simulations, chi versus concentration plots were constructed, which showed trends similar to those experimentally measured melting temperature versus concentration. The chi values for the blends, which satisfied the criteria of miscibility of two polymers by the atomistic simulation, agreed quite well with the solubility criteria related to order parameters calculated from the mesoscopic simulation. Kinetics of phase separation was examined via density profiles calculated using the MesoDyn approach for incompatible blends. The length and time scales spanned by these simulations were found to be relevant to the real application scales. The free energy computed in the mesoscopic simulation for blends reached equilibrium, particularly when the simulation was performed at a higher time step, indicating the stability of the blend system at certain compositions.