Plasma-facing components in tokamak-type fusion reactors are subjected to intense heat loads during plasma disruptions. The influence of high heat fluxes on the depth of heat-affected zones of pure beryllium metal and beryllium containing very low levels of surface active impurities is studied by using a two-dimensional transient computer model that solves the equations of motion and energy. Results are presented for a range of energy densities and disruption times. Under certain conditions, impurities, through their effect on surface tension, create convective flows and hence influence the flow intensities and the resulting depths of the beryllium molten layers during plasma disruptions. The calculated depths of the molten layers are also compared with other mathematical models that are based on the assumption that heat is transported through the material by conduction only.