Numerical modelling of the operation of a two-phase thermosyphon