An explicit fully discrete finite element method, which satisfies the nonconvex side constraint at every node, is developed for approximating the p-harmonic flow for p ∈ (1,∞). Convergence of the method is established under certain conditions on the domain and mesh. Computational examples are presented to demonstrate finite-time blow-ups and qualitative geometric changes of weak solutions of th...
