Projection-free approximation of flows of harmonic maps with quadratic constraint accuracy and variable step sizes

TL;DR

This paper introduces a projection-free, energy-stable numerical method for simulating harmonic map flows into spheres, achieving second-order accuracy and efficient convergence with variable step sizes.

Contribution

It presents a novel linearly implicit, projection-free scheme that maintains energy stability and improves accuracy near singularities, with adaptive step sizing.

Findings

Achieves second-order accuracy under certain regularity conditions

Maintains unconditional energy stability

Outperforms traditional Euler and BDF methods in experiments

Abstract

We construct and analyze a projection-free linearly implicit method for the approximation of flows of harmonic maps into spheres. The proposed method is unconditionally energy stable and, under a sharp discrete regularity condition, achieves second-order accuracy with respect to the constraint violation. Furthermore, the method accommodates variable step sizes to speed up the convergence to stationary points and to improve the accuracy of the numerical solutions near singularities, without affecting the unconditional energy stability and the constraint violation property. We illustrate the accuracy in approximating the unit-length constraint and the performance of the method through a series of numerical experiments, and compare it with the linearly implicit Euler and two-step BDF methods.