doi:10.1016/j.amc.2005.01.050
Copyright © 2005 Elsevier Inc. All rights reserved.
The Fourier spectral method for the Cahn–Hilliard equation
Department of Mathematics, Zhejiang University, Hangzhou 310028, PR China
Available online 3 March 2005.
References and further reading may be available for this article. To view references and further reading you must
purchase this article.
Abstract
In this paper, a Fourier spectral method for numerically solving Cahn–Hilliard equation with periodic boundary conditions is developed. We establish their semi-discrete and fully discrete schemes that inherit the energy dissipation property and mass conservation property from the associated continuous problem. we prove existence and uniqueness of the numerical solution and derive the optimal error bounds. we perform some numerical experiments which confirm our results.
Keywords: Cahn–Hilliard equation; Spectral method
Fig. 1. (a) and (b) show the evolution from t = 0 to t = 5. From these plots, we can easily see that the total mass remains indeed constant since the initial data has null mean value
. In numerical experiments, we observed that the solution u(t) tends to,and remains at for a long time, patterns that nearly piecewise constant, the steady-state solution is very close to a piecewise constant function. In [11], Elliott and French also observed the same results. In this experiment, the patterns hardly changed after t = 3.
Fig. 2. (a) shows the time dependency of energy functional (1.6) of numerical solutions in Fig. 1. As shown in the energy functional dissipation property, the energy functional
of numerical solutions theoretically decreases as time pass (see lk and lk). (a) indicates that the energy functional
of numerical solutions decreases and agree with the dissipation property. (b) shows the time dependency of mass of numerical solutions in Fig. 1. The mass of numerical solutions is theoretically independent of time and this is shown in lk and lk. Looking at (b), the mass of numerical solution is conservation quite well and it agrees with the conservation property also.
Table 1.
