A numerical method for tracking curve networks moving with curvature motion

Abstract A finite difference method is proposed to track curves whose normal velocity is given by their curvature and which meet at different types of junctions. The prototypical example is that of phase interfaces that meet at prescribed angles, although eutectic junctions and interactions through nonlocal effects are also considered. The method is based on a direct discretization of the underlying parabolic problem and boundary conditions. A linear stability analysis is presented for our scheme as well as computational studies that confirm the second order convergence to smooth solutions. After a singularity in the curve network where the solution is no longer smooth, we demonstrate "almost" second-order convergence. A numerical study of singularity types is done for the case of networks that meet at prescribed angles at triple junctions. Finally, different discretizations and methods for implicit time stepping are presented and compared.

[1]  C. S. Smith,et al.  Grains, Phases, and Interfaces an Interpretation of Microstructure , 1948 .

[2]  W. Mullins.  Two‐Dimensional Motion of Idealized Grain Boundaries , 1956 .

[3]  Stability of Difference Approximations to the Mixed Initial Boundary Value Problems for Parabolic Systems , 1971 .

[4]  John W. Cahn,et al.  Critical point wetting , 1977 .

[5]  Computer simulation of bidimensional grain growth , 1984 .

[6]  Wolfgang Hackbusch,et al.  Multi-grid methods and applications , 1985, Springer series in computational mathematics.

[7]  M. Grayson.  The heat equation shrinks embedded plane curves to round points , 1987 .

[8]  Short-range order in the arrangement of grains in two-dimensional polycrystals , 1987 .

[9]  C. V. Thompson,et al.  A two-dimensional computer simulation of capillarity-driven grain growth: Preliminary results , 1988 .

[10]  Katsuya Nakashima,et al.  Vertex models for two-dimensional grain growth , 1989 .

[11]  J. Strain.  A boundary integral approach to unstable solidification , 1989 .

[12]  James A. Glazier,et al.  Coarsening in the two-dimensional soap froth and the large-Q Potts model: A detailed comparison , 1990 .

[13]  J. Rubinstein,et al.  Nonlocal reaction−diffusion equations and nucleation , 1992 .

[14]  L. Bronsard,et al.  On three-phase boundary motion and the singular limit of a vector-valued Ginzburg-Landau equation , 1993 .

[15]  M. Glicksman,et al.  Topological rearrangements during 2D normal grain growth , 1993 .

[16]  L. Bronsard,et al.  Volume-preserving mean curvature flow as a limit of a nonlocal Ginzburg-Landau equation , 1997 .