Abstract
A numerical method and software were developed to predict non-stationary conjugated conductive heat transfer, melting and possible evaporation of materials under high energy fluxes impinging onto solid body surface (plasma jet, arc spot, laser or electron beam), and also subsequent cooling and solidification of the melted substrate layer. In the numerical procedure, the finite-element method was employed. The processes of interest can have different characteristic time and spatial scales, which in addition can suffer drastic changes at heat flux densities q ?[108 ;1014 ] W/m2. An advanced procedure was developed to enable dynamic adaptive triangulation of domains involved in the current numerical solution and characterizing the different phase states (liquid or solid) of the materials. This procedure, belonging to the class of the frontal algorithms, allows one to break a solution domain into triangles based only on the domain boundaries. The model applications of the developed simulation software are illustrated.