Department of Mathematics
 Search | Help | Login | pdf version | printable version

Math @ Duke



Publications [#329008] of Saulo Orizaga

Papers Published

  1. Glasner, K; Orizaga, S, Improving the accuracy of convexity splitting methods for gradient flow equations, Journal of Computational Physics, vol. 315 (June, 2016), pp. 52-64, Elsevier BV [doi]
    (last updated on 2019/10/13)

    © 2016 Elsevier Inc. This paper introduces numerical time discretization methods which significantly improve the accuracy of the convexity-splitting approach of Eyre (1998) [7], while retaining the same numerical cost and stability properties.A first order method is constructed by iteration of a semi-implicit method based upon decomposing the energy into convex and concave parts. A second order method is also presented based on backwards differentiation formulas. Several extrapolation procedures for iteration initialization are proposed. We show that, under broad circumstances, these methods have an energy decreasing property, leading to good numerical stability.The new schemes are tested using two evolution equations commonly used in materials science: the Cahn-Hilliard equation and the phase field crystal equation. We find that our methods can increase accuracy by many orders of magnitude in comparison to the original convexity-splitting algorithm. In addition, the optimal methods require little or no iteration, making their computation cost similar to the original algorithm.
ph: 919.660.2800
fax: 919.660.2821

Mathematics Department
Duke University, Box 90320
Durham, NC 27708-0320