Multi-dimensional hydrodynamic code SOVA for interfacial ?ows:Application to the thermal layer effect. Eulerian, three-dimensional, numerical code, which conserves mass, momentum and energy si-multaneously both in the Lagrangian and remap steps, is developed. The use of special form of linearviscosity provides a weaker time step restriction as compared with the Courant condition. The code isdesigned to investigate the multi-material problems, including dusty ?ows. The performance of the code isillustrated by the modeling of shock wave interaction with a dusty thermal layer.