Numerical and experimental works have been conducted to develop a visualization technique for phase distribution in two-phase systems by an electrical impedance tomography technique, which reconstructs the resistivity distribution with the electrical responses that are determined by corresponding excitations. In the conventional Newton-Raphson method, the impedance distribution, i.e., the phase distribution, is usually reconstructed in fixed elements inside the system. In practical cases, such as the impedance imaging of a two-phase flow, this model might not be useful. For two-phase flow systems, the impedance of each phase does not change, but, instead, the phase boundary depends on the distribution of the dispersed phase. In the present study, a new image reconstruction algorithm employing an adaptive mesh regeneration technique was developed for the detection of a phase boundary. The feasibility of this method is tested using some phantom experiments.