This is Matlab code, with some C subroutines (mex files), that implements grain boundary motion of networks with arbitrary surface tensions and fixed mobilities, in two space dimensions. Note carefully that the mobility of each interface = 1 / its surface tension, i.e.
The algorithm is based on the following paper:
Esedoglu, S.; Otto, F. Threshold dynamics for networks with arbitrary surface tensions. Communications on Pure and Applied Mathematics. 68:5 (2015), pp. 808-864.
The research that led to this software was supported by the National Science Foundation grant DMS-074333.
|Version 1||November 15, 2015.|
After downloading the linked zipped directory, unzip it. For each C source file contained in the directory, issue the command
from within Matlab, where
in Matlab. Here is an example of how to call the program in Matlab:
This will generate randomly (via Voronoi construction) 100 initial grains in three dimensions on a