A numerical code was developed for modeling phase change micro heat exchangers. The model covers evaporating and condensing annular flows as well as slug/elongated bubble flows. The model was validated against a set of experimental results from different laboratories and depicted good agreement with the data. The model takes puts a particular emphasis at optimizing the geometry and operating conditions of micro heat exchagers in a fast and reliable manner.