Abstract. Multifrequency Electrical Impedance Tomography (MFEIT) reconstructs the distribution of conductivity by exploiting the dependence of tissue conductivity on frequency. MFEIT can be performed on a single instance of data, making it promising for applications such as stroke and cancer imaging, where it is not possible to obtain a 'baseline' measurement of healthy tissue. A nonlinear MFEIT algorithm able to reconstruct the volume fraction distribution of tissue rather than conductivities has been developed previously. For each volume, the fraction of a certain tissue should be either 1 or 0; this implies that the sharp changes of the fractions, representing the boundaries of tissue, contain all the relevant information However, these boundaries are blurred by traditional regularisation methods employing 2 l norm. The Total Variation (TV) regularisation can overcome this problem, but it is difficult to solve due to its non-differentiability. As the fraction must be between 0 and 1, this imposes a constraint on the MFEIT method based on the fraction model. Therefore, a constrained optimisation method capable of dealing with non-differentiable problems is required. We propose a new constrained TV regularised method, to solve the fraction reconstruction problem, based on the Primal and Dual Interior Point Method (PDIPM) method. The noise performance of the new MFEIT method is analysed using simulations on a 2D cylindrical mesh. Convergence performance is also analysed, through experiments using a cylindrical tank. Finally, simulations on an anatomically realistic head-shaped mesh are demonstrated. The proposed MFEIT method with TV regularisation shows higher spatial resolution, particularly at the edges of the perturbation, and stronger noise robustness, and its image noise and shape error are 20% to 30% lower than the traditional fraction method.