A new method for the molecular dynamics simulation of bulk planar interfaces at high temperatures is presented. The method uses the basic Parrinallo-Rahman (constant stress) scheme modified, however, for the application to inhomogeneous systems. Since our computational cell contains only one interface with 2-d periodic boundary conditions we are able to study isolated interfaces all the way up to melting. The interaction between boundaries which may lead to their annihilation at higher temperatures, which is a problem when 3-d periodic boundaries are applied, is thus avoided. As an application, the method is used to study the stability of a grain boundary at high temperatures. Observations on a possible connection between grain boundary migration and “premelting” are discussed.