An algorithm for the computation of design variables in layer crystallisation processes is presented, based on an algebraic solution of the moving boundary equation for heat transfer, as well as on an expression for the concentration profile. The algorithm was applied to the discontinuous crystallisation of a NaCl-H2O solution around the wall of an internally cooled cylindrical tube. Good agreement between simulation results and laboratory experiments was observed. After input of experimentally measurable values of the distribution coefficient, the algorithm enables the prediction of the following process variables along the crystallisation time: layer thickness, temperature, and composition profiles across the layer, as well as the end-values of solid and liquid mass and concentrations. Multistage processes can be studied by the series association of crystallisation cycles, enabling the determination of the optimal number of stages that meets a specific design criterion.