A general framework for high accuracy solutions to energy gradient flows from material science models

SCAIM Seminar
January 29, 2013 8:30 pm

Speaker:  Brian Wetton, Department of Mathematics, University of British Columbia

URL for Speaker:  http://www.math.ubc.ca/~wetton/

Location:  ESB 4133

Intended Audience:  Public

A computational framework is presented for materials science models that come from energy gradient flows that lead to the evolution of structure involving two or more phases. The models are considered in periodic cells and standard Fourier spectral discretization in space is used. Implicit time stepping is used, and the resulting implicit systems are solved iteratively with the preconditioned conjugate gradient method. The dependence of the condition number of the preconditioned system on the size of the time step and the order parameter in the model (that represents the scaled width of transition layers between phases) is investigated. The framework is easily extended to higher order derivative models, higher dimensional settings, and vector problems. Several examples of its application are demonstrated, including a sixth order problem in three dimensions. Higher order time stepping and a GPU implementation is described briefly. A comparison to time-stepping with operator splitting (into convex and concave parts) is done.