Unconstrained Scheme for Geometrically Constrained Gradient Flows

In this paper, we study the approximation of gradient flows of harmonic maps, which serve as model problems for applications in micromagnetics, liquid crystals, and nonlinear plate bending.

Harmonic maps are vector fields that are critical points of the Dirichlet energy subject to the constraint that the vector field be unit length pointwise.

Most existing time-stepping schemes for gradient flows deal with the constraint by linearizing the unit length constraint at every step, which involves solving for the solution increment in the tangent space of the constraint.

These schemes lead to robust control over the violation of the constraint, but require solving degenerate saddle point systems at every step that may be difficult to precondition.

In this paper, we propose a scheme that first computes the unconstrained increment and then projects this increment pointwise onto the tangent space.

With an additional stabilization, this scheme is energy stable under mild step size restrictions and provides robust control of the unit length constraint violation.

Our new scheme only requires the solution of decoupled symmetric positive definite systems at every step, which translates to a large increase in computational efficiency.

We also propose a computable a posteriori criterion and a variable time-stepping procedure that guarantee the stability of the scheme.

We conclude with computational examples demonstrating the efficacy of the scheme, and present a computational extension of the scheme to nonlinear plate bending.