Structure-preserving numerical scheme for a generalized area-preserving crystalline curvature flow

The presented numerical scheme preserves variational structure of a generalized area-preserving crystalline curvature flow. The scheme is based on an iteration and a projection method. Several numerical examples will verify that the enclosed area is preserved in numerical computation with high accuracy in the sense of double precision. Numerical computations realize theoretical convexification results starting from almost convex polygon, and are extended to a general setting starting from nonconvex polygon.