A generalized MBO diffusion generated motion for orthogonal matrix-valued fields

11/03/2017
by   Braxton Osting, et al.
0

We consider the problem of finding stationary points of the Dirichlet energy for orthogonal matrix-valued fields. Following the Ginzburg-Landau approach, this energy is relaxed by penalizing the matrix-valued field when it does not take orthogonal matrix values. A generalization of the MBO diffusion generated motion is introduced that effectively finds local minimizers of this energy by iterating two steps until convergence. In the first step, as in the original method, the current matrix-valued field is evolved by the diffusion equation. In the second step, the field is pointwise reassigned to the closest orthogonal matrix, which can be computed via the singular value decomposition. We extend the Lyapunov function of Esedoglu and Otto to show that the method is non-increasing on iterates and hence, unconditionally stable. We also prove that spatially discretized iterates converge to a stationary solution in a finite number of iterations. The algorithm is implemented using the closest point method and non-uniform fast Fourier transform. We conclude with several numerical experiments on flat tori and closed surfaces, which, unsurprisingly, exhibit classical behavior from the Allen-Cahn and complex Ginzburg Landau equations, but also new phenomena.

READ FULL TEXT

page 17

page 18

page 20

page 21

page 24

page 25

page 27

page 29

research
08/25/2021

Energy stability of Strang splitting for vector-valued and matrix-valued Allen-Cahn equations

We consider the second-order in time Strang-splitting approximation for ...
research
01/24/2020

Electric Field Propagation Through Singular Value Decomposition

We demonstrate that the singular value decomposition algorithm in conjun...
research
12/19/2019

Equilibria of an anisotropic nonlocal interaction equation: Analysis and numerics

In this paper, we study the equilibria of an anisotropic, nonlocal aggre...
research
02/17/2022

On H1, H(curl) and H(sym Curl) finite elements for matrix-valued Curl problems

In this work we test the numerical behaviour of matrix-valued fields app...
research
01/05/2014

Schatten-p Quasi-Norm Regularized Matrix Optimization via Iterative Reweighted Singular Value Minimization

In this paper we study general Schatten-p quasi-norm (SPQN) regularized ...
research
05/10/2020

An efficient diffusion generated motion method for wetting dynamics

By using the Onsager variational principle as an approximation tool, we ...
research
11/08/2018

Orthogonal Trace-Sum Maximization: Applications, Local Algorithms, and Global Optimality

This paper studies a problem of maximizing the sum of traces of matrix q...

Please sign up or login with your details

Forgot password? Click here to reset