The continuous-properties model has been developed to numerically solve problems of melting and solidification of substances with a mushy region. The new method is applicable to negligibly small mushy regions. Therefore, it can also be used to numerically calculate the behaviour of materials showing a discontinuous (first-order) phase transition. A stress number, comparable to the Reynolds number of fluid flow, can be introduced to quantify the difficulty of the computational task.