The heat equation to generate planar grids

An algorithm for the generation of quadrilateral grids on planar domains is presented. This algorithm is given by an iterative procedure, where, starting from an initial grid on the domain under consideration, the coordinates of the grid vertices are iteratively adjusted by using a local discrete variational approach. This procedure resembles the explicit difference scheme for a perturbed heat equation, where the perturbation can be dropped for convex domains. Experimental results on benchmark domains are presented, and show an interesting behavior of the proposed method.