A numerical scheme based on a coordinate transform into stream function-velocity potential is proposed to solve heat and momentum transfer in porous media with phase change. A significant simplification of both computational domain and governing equations can be achieved by the transform. The dispersion term in the flow through porous media, which is important at the phase change interface, can be successfully incorporated into the numerical scheme without introducing any further computational complications.