The mathematical simulation of crystallization process is based on numerical solution of the non-stationary Stephan problem, which is a labor-consuming procedure. A parallel algorithm for the Stephen problem solution within one time layer based on the physical processes decomposition scheme has been proposed. Such approach allows up to 5 times speedup, which is especially important for multi-parameterized simulation and research.