To reduce the overhead involved with the solution of the coupled system of physical and mesh PDEs, we propose a family of parallel solvers based on domain decomposition approaches. In this short talk I will review theoretical results obtained for the convergence of Schwarz methods for the steady state equidistribution principle at both the continuous and discrete levels and for a particular MMPDE solved by Rothe’s method.