Skip to main content
U.S. Department of Energy
Office of Scientific and Technical Information

Advanced operator splitting-based semi-implicit spectral method to solve the binary phase-field crystal equations with variable coefficients

Journal Article · · Journal of Computational Physics
 [1];  [1]; ;  [2]; ;  [1]
  1. Brunel Centre for Advanced Solidification Technology, Brunel University, Uxbridge, Middlesex UB8 3PH (United Kingdom)
  2. Research Institute for Solid State Physics and Optics, P.O. Box 49, H-1525 Budapest (Hungary)
We present an efficient method to solve numerically the equations of dissipative dynamics of the binary phase-field crystal model proposed by Elder et al. [K.R. Elder, M. Katakowski, M. Haataja, M. Grant, Phys. Rev. B 75 (2007) 064107] characterized by variable coefficients. Using the operator splitting method, the problem has been decomposed into sub-problems that can be solved more efficiently. A combination of non-trivial splitting with spectral semi-implicit solution leads to sets of algebraic equations of diagonal matrix form. Extensive testing of the method has been carried out to find the optimum balance among errors associated with time integration, spatial discretization, and splitting. We show that our method speeds up the computations by orders of magnitude relative to the conventional explicit finite difference scheme, while the costs of the pointwise implicit solution per timestep remains low. Also we show that due to its numerical dissipation, finite differencing can not compete with spectral differencing in terms of accuracy. In addition, we demonstrate that our method can efficiently be parallelized for distributed memory systems, where an excellent scalability with the number of CPUs is observed.
OSTI ID:
21167758
Journal Information:
Journal of Computational Physics, Journal Name: Journal of Computational Physics Journal Issue: 5 Vol. 228; ISSN JCTPAH; ISSN 0021-9991
Country of Publication:
United States
Language:
English

Similar Records

The incorporation of the semi-implicit linear equations into Newton's method to solve radiation transfer equations
Journal Article · Mon Sep 10 00:00:00 EDT 2007 · Journal of Computational Physics · OSTI ID:21028268

An unconditionally stable, time-implicit algorithm for solving the one-dimensional Vlasov–Poisson system
Journal Article · Mon Mar 07 19:00:00 EST 2022 · Journal of Plasma Physics · OSTI ID:1893545