Operation of the fuel cycle of a thermonuclear fusion reactor naturally leads to accumulation of surplus protium, but in some cases it can also lead to accumulation of surplus deuterium. Both surplus protium and deuterium have to be separated, detritiated, and discharged to the environment, normally passing a final detritiation stage based on either the liquid phase catalytic exchange or water distillation process. The concept of a multicolumn cryogenic distillation (CD) system capable of discharging a time-varying surplus of deuterium is presented in this paper. A model of a CD column based on a UV (internal energy U – volume V) flash formulation and equation of state (EOS) thermodynamic model for hydrogen isotopologue mixtures is also presented at the principal step to a comprehensive model of the isotope separation system. Although fundamental for constant volume systems, the UV formulation of the thermodynamic state has not been widely used in transient simulations; in particular, for distillation dynamics modeling, other approaches are much more common. At the same time, in helium cryogenics the UV formulation has gained wide usage in large-scale dynamic simulations. It is known from the literature that a UV formulation of the distillation problem is very challenging for a numerically stable implementation. To cope with this situation, we present our findings on the sources of numerical instabilities and approaches.