Atomistic simulations were performed on poly(N-vinyl-2-pyrrolidone) (PVP) and its blends with chitosan (CS) in different ratios using molecular mechanics (MM) and molecular dynamics (MD) simulations in three-dimensionally periodic and effective two-dimensionally periodic condensed phases. Four independent microstructures were generated to analyze their surface properties. The calculated surface-energy values for PVP compared quite well with the experimental data reported in the literature. The density profile was analyzed, and the structure of the films showed an interior region of the bulk density. Various components of the energetic interactions (torsional, van der Waals, etc.) were examined to gain deeper insight into the nature of regular and anomalous interactions between the bulk and the surface films. Surface energies of PVP/CS blends were computed by MD simulations using the bulk pressure-volume-temperature (PVT) parameters. Bulk properties such as the cohesive energy density (CED) and solubility parameter (delta) were calculated using MM and MD simulations in the NVT ensemble under periodic boundary conditions. The Flory equation of state was used to compute the thermal expansion coefficient as well as PVT parameters. These surface-energy values agreed well with the surface-energy data calculated using the Zisman equation, which were also in accordance with the experimental observations. The results from this study suggest that computer simulations would provide valuable information on polymers and polymer-blend surfaces.