Long Time Numerical Simulations for Phase-Field Problems Using $p$-Adaptive Spectral Deferred Correction Methods

2015 ◽  
Vol 37 (1) ◽  
pp. A271-A294 ◽  
Author(s):  
Xinlong Feng ◽  
Tao Tang ◽  
Jiang Yang
Keyword(s):  
Numerical Simulations ◽  
Phase Field ◽  
Deferred Correction ◽  
Long Time ◽  
Spectral Deferred Correction

IMPLEMENTING ARBITRARILY HIGH-ORDER SYMPLECTIC METHODS VIA KRYLOV DEFERRED CORRECTION TECHNIQUE

2010 ◽  
Vol 01 (02) ◽  
pp. 277-301 ◽  
Author(s):  
QUANDONG FENG ◽  
JINGFANG HUANG ◽  
NINGMING NIE ◽  
ZAIJIU SHANG ◽  
YIFA TANG
Keyword(s):  
Hamiltonian Systems ◽  
Gaussian Quadrature ◽  
Numerical Procedure ◽  
Original System ◽  
Algebraic Equations ◽  
Initial Value ◽  
Deferred Correction ◽  
Correction Technique ◽  
Long Time ◽  
Spectral Deferred Correction

In this paper, an efficient numerical procedure is presented to implement the Gaussian Runge–Kutta (GRK) methods (also called Gauss methods). The GRK technique first discretizes each marching step of the initial value problem using collocation formulations based on Gaussian quadrature. As is well known, it preserves the geometric structures of Hamiltonian systems. Existing analysis shows that the GRK discretization with s nodes is of order 2s, A-stable, B-stable, symplectic and symmetric, and hence "optimal" for solving initial value problems of general ordinary differential equations (ODEs). However, as the unknowns at different collocation points are coupled in the discretized system, direct solution of the resulting algebraic equations is in general inefficient. Instead, we use the Krylov deferred correction (KDC) method in which the spectral deferred correction (SDC) scheme is applied as a preconditioner to decouple the original system, and the resulting preconditioned nonlinear system is solved efficiently using Newton–Krylov schemes such as Newton–GMRES method. The KDC accelerated GRK methods have been applied to several Hamiltonian systems and preliminary numerical results are presented to show the accuracy, stability, and efficiency features of these methods for different accuracy requirements in short- and long-time simulations.

A Spectral Deferred Correction Method Applied to the Shallow Water Equations on a Sphere

2013 ◽  
Vol 141 (10) ◽  
pp. 3435-3449 ◽  
Author(s):  
Jun Jia ◽  
Judith C. Hill ◽  
Katherine J. Evans ◽  
George I. Fann ◽  
Mark A. Taylor
Keyword(s):  
Shallow Water ◽  
Shallow Water Equations ◽  
Correction Method ◽  
Higher Order ◽  
High Order ◽  
Benchmark Problems ◽  
Time Step ◽  
Deferred Correction ◽  
Long Time ◽  
Spectral Deferred Correction

Abstract Although significant gains have been made in achieving high-order spatial accuracy in global climate modeling, less attention has been given to the impact imposed by low-order temporal discretizations. For long-time simulations, the error accumulation can be significant, indicating a need for higher-order temporal accuracy. A spectral deferred correction (SDC) method is demonstrated of even order, with second- to eighth-order accuracy and A-stability for the temporal discretization of the shallow water equations within the spectral-element High-Order Methods Modeling Environment (HOMME). Because this method is stable and of high order, larger time-step sizes can be taken while still yielding accurate long-time simulations. The spectral deferred correction method has been tested on a suite of popular benchmark problems for the shallow water equations, and when compared to the explicit leapfrog, five-stage Runge–Kutta, and fully implicit (FI) second-order backward differentiation formula (BDF2) time-integration methods, it achieves higher accuracy for the same or larger time-step sizes. One of the benchmark problems, the linear advection of a Gaussian bell height anomaly, is extended to run for longer time periods to mimic climate-length simulations, and the leapfrog integration method exhibited visible degradation for climate length simulations whereas the second-order and higher methods did not. When integrated with higher-order SDC methods, a suite of shallow water test problems is able to replicate the test with better accuracy.

scholarly journals Well-posedness and long-time behavior for the modified phase-field crystal equation

2014 ◽  
Vol 24 (14) ◽  
pp. 2743-2783 ◽  
Author(s):  
Maurizio Grasselli ◽  
Hao Wu
Keyword(s):  
Phase Field ◽  
Time Derivative ◽  
Exponential Attractor ◽  
Well Posedness ◽  
Time Behavior ◽  
Long Time Behavior ◽  
Long Time ◽  
Global Existence And Uniqueness ◽  
And Diffusion ◽  
Phase Field Crystal

We consider a modification of the so-called phase-field crystal (PFC) equation introduced by K. R. Elder et al. This variant has recently been proposed by P. Stefanovic et al. to distinguish between elastic relaxation and diffusion time scales. It consists of adding an inertial term (i.e. a second-order time derivative) into the PFC equation. The mathematical analysis of the resulting equation is more challenging with respect to the PFC equation, even at the well-posedness level. Moreover, its solutions do not regularize in finite time as in the case of PFC equation. Here we analyze the modified PFC (MPFC) equation endowed with periodic boundary conditions. We first prove the global existence and uniqueness of a solution with initial data in a bounded energy space. This solution satisfies some uniform dissipative estimates which allow us to study the long-time behavior of the corresponding dynamical system. In particular, we establish the existence of the global attractor as well as an exponential attractor. Then we demonstrate that any trajectory originating from the bounded energy phase space converges to a single equilibrium. This is done by means of a suitable version of the Łojasiewicz–Simon inequality. An estimate on the convergence rate is also given.

Efficient numerical algorithms for the phase field crystal equation

2021 ◽  
pp. 2150042
Author(s):  
Qingqu Zhuang ◽  
Shuying Zhai ◽  
Zhifeng Weng
Keyword(s):  
Free Energy ◽  
Lower Bound ◽  
Numerical Simulations ◽  
Lagrange Multiplier ◽  
Phase Field ◽  
Numerical Algorithms ◽  
Energy Potential ◽  
Energy Stable ◽  
2D And 3D ◽  
Phase Field Crystal

In this paper, based on the Lagrange Multiplier approach in time and the Fourier-spectral scheme for space, we propose efficient numerical algorithms to solve the phase field crystal equation. The numerical schemes are unconditionally energy stable based on the original energy and do not need the lower bound hypothesis of the nonlinear free energy potential. The unconditional energy stability of the three semi-discrete schemes is proven. Several numerical simulations in 2D and 3D are demonstrated to verify the accuracy and efficiency of our proposed schemes.

Sign in / Sign up

Export Citation Format

Share Document