Applications of semi-implicit Fourier-spectral method to phase field equations

Long-qing Chen, Jie Shen

Research output: Contribution to journalArticle

617 Citations (Scopus)

Abstract

An efficient and accurate numerical method is implemented for solving the time-dependent Ginzburg-Landau equation and the Cahn-Hilliard equation. The time variable is discretized by using semi-implicit schemes which allow much larger time step sizes than explicit schemes; the space variables are discretized by using a Fourier-spectral method whose convergence rate is exponential in contrast to second order by a usual finite-difference method. We have applied our method to predict the equilibrium profiles of an order parameter across a stationary planar interface and the velocity of a moving interface by solving the time-dependent Ginzburg-Landau equation, and compared the accuracy and efficiency of our results with those obtained by others. We demonstrate that, for a specified accuracy of 0.5%, the speedup of using semi-implicit Fourier-spectral method, when compared with the explicit finite-difference schemes, is at least two orders of magnitude in two dimensions, and close to three orders of magnitude in three dimensions. The method is shown to be particularly powerful for systems in which the morphologies and microstructures are dominated by long-range elastic interactions.

Original languageEnglish (US)
Pages (from-to)147-158
Number of pages12
JournalComputer Physics Communications
Volume108
Issue number2-3
StatePublished - Feb 1 1998

Fingerprint

spectral methods
Finite difference method
Numerical methods
Microstructure
Landau-Ginzburg equations
microstructure
profiles
interactions

All Science Journal Classification (ASJC) codes

  • Computer Science Applications
  • Physics and Astronomy(all)

Cite this

@article{0afec671b2734aa8bfe63597bca585ca,
title = "Applications of semi-implicit Fourier-spectral method to phase field equations",
abstract = "An efficient and accurate numerical method is implemented for solving the time-dependent Ginzburg-Landau equation and the Cahn-Hilliard equation. The time variable is discretized by using semi-implicit schemes which allow much larger time step sizes than explicit schemes; the space variables are discretized by using a Fourier-spectral method whose convergence rate is exponential in contrast to second order by a usual finite-difference method. We have applied our method to predict the equilibrium profiles of an order parameter across a stationary planar interface and the velocity of a moving interface by solving the time-dependent Ginzburg-Landau equation, and compared the accuracy and efficiency of our results with those obtained by others. We demonstrate that, for a specified accuracy of 0.5{\%}, the speedup of using semi-implicit Fourier-spectral method, when compared with the explicit finite-difference schemes, is at least two orders of magnitude in two dimensions, and close to three orders of magnitude in three dimensions. The method is shown to be particularly powerful for systems in which the morphologies and microstructures are dominated by long-range elastic interactions.",
author = "Long-qing Chen and Jie Shen",
year = "1998",
month = "2",
day = "1",
language = "English (US)",
volume = "108",
pages = "147--158",
journal = "Computer Physics Communications",
issn = "0010-4655",
publisher = "Elsevier",
number = "2-3",

}

Applications of semi-implicit Fourier-spectral method to phase field equations. / Chen, Long-qing; Shen, Jie.

In: Computer Physics Communications, Vol. 108, No. 2-3, 01.02.1998, p. 147-158.

Research output: Contribution to journalArticle

TY - JOUR

T1 - Applications of semi-implicit Fourier-spectral method to phase field equations

AU - Chen, Long-qing

AU - Shen, Jie

PY - 1998/2/1

Y1 - 1998/2/1

N2 - An efficient and accurate numerical method is implemented for solving the time-dependent Ginzburg-Landau equation and the Cahn-Hilliard equation. The time variable is discretized by using semi-implicit schemes which allow much larger time step sizes than explicit schemes; the space variables are discretized by using a Fourier-spectral method whose convergence rate is exponential in contrast to second order by a usual finite-difference method. We have applied our method to predict the equilibrium profiles of an order parameter across a stationary planar interface and the velocity of a moving interface by solving the time-dependent Ginzburg-Landau equation, and compared the accuracy and efficiency of our results with those obtained by others. We demonstrate that, for a specified accuracy of 0.5%, the speedup of using semi-implicit Fourier-spectral method, when compared with the explicit finite-difference schemes, is at least two orders of magnitude in two dimensions, and close to three orders of magnitude in three dimensions. The method is shown to be particularly powerful for systems in which the morphologies and microstructures are dominated by long-range elastic interactions.

AB - An efficient and accurate numerical method is implemented for solving the time-dependent Ginzburg-Landau equation and the Cahn-Hilliard equation. The time variable is discretized by using semi-implicit schemes which allow much larger time step sizes than explicit schemes; the space variables are discretized by using a Fourier-spectral method whose convergence rate is exponential in contrast to second order by a usual finite-difference method. We have applied our method to predict the equilibrium profiles of an order parameter across a stationary planar interface and the velocity of a moving interface by solving the time-dependent Ginzburg-Landau equation, and compared the accuracy and efficiency of our results with those obtained by others. We demonstrate that, for a specified accuracy of 0.5%, the speedup of using semi-implicit Fourier-spectral method, when compared with the explicit finite-difference schemes, is at least two orders of magnitude in two dimensions, and close to three orders of magnitude in three dimensions. The method is shown to be particularly powerful for systems in which the morphologies and microstructures are dominated by long-range elastic interactions.

UR - http://www.scopus.com/inward/record.url?scp=0001152970&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=0001152970&partnerID=8YFLogxK

M3 - Article

VL - 108

SP - 147

EP - 158

JO - Computer Physics Communications

JF - Computer Physics Communications

SN - 0010-4655

IS - 2-3

ER -