Abstract

The diffusion equation is discretized in spacial direction and transformed into the ordinary differential equations. The ordinary differential equations are solved by Lie group method and the explicit Runge-Kutta method. Numerical results showed that Lie group method is more stable than the corresponding explicit Runge-Kutta method.

1. Introduction

The subject of solving stiff system has been well studied in the literature, including linearly implicit methods, semi-implicit methods, time-splitting methods, projection methods, multiscale methods, integrating factor methods, the exponential time difference (ETD) methods, and Lie group methods [112]. Lie group methods have been developed by Munthe-Kaas, Iserles, Nørsett, and their collaborators. In general, if a Lie group acts freely and transitively on a manifold , then a differential equation on uniquely determines a differential equation on . Thus, replacing the original differential on by the equivalent differential equation on Lie algebra allows one to apply appropriate Lie group integrators. There are many Lie group methods, such as Mangus method, RKMK method, and Fer method. Lie group methods can preserve the numerical solutions of the differential equations on the same manifolds, which have good stability and the same accuracy as the classical numerical methods [1, 2, 11, 13, 14]. In the paper, the exponential integrator method of the general nonlinear differential equations is proposed based on the idea of Lie group methods on manifolds and applied to the diffusion equations.

The paper is organized as follows: in Section 2, the Lie group method and the exponential time difference method are introduced. The exponential integrator method for the general nonlinear differential equations is proposed. In Section 3, the quasilinear diffusion equation is solved by the explicit exponential integrator method and the explicit Runge-Kutta methods. In Section 4, the Allen-Cahn equation is solved by the explicit exponential integrator method and the explicit Runge-Kutta methods. At last, some conclusions are obtained.

2. Exponential Integrator

Suppose that is a manifold. There exist a Lie algebra with a Lie bracket , a (left) Lie algebra action , and a function such that the equation for can be written as where is a homomorphism. The left algebra action assigns to the arguments and the solution of at time with the initial condition and some boundary conditions. There is freedom to choose the function , the corresponding differential equation in the Lie algebra Lie algebra is chosen as functions corresponding to simpler vector fields on which are easy to exponentiate. The correspondence between the Lie algebra and the vector fields on will be one to one. The Runge-Kutta method is applied to solve (4). It is the well-known RKMK method [1].

Let and be the coefficients of an -stage, order classical Runge-Kutta method, and let . The algorithm integrates (1) from to

Set .

For End where , and . The number denotes the order of the underlying Runge-Kutta method and is the th Bernoulli number [2].

A system of ordinary differential equations is often obtained where . The ETD method can be described in the context of solving (8). Integrating the equation over a single time step from to , we can get Equation (9) is exact, and the various ETD schemes come from the approximation to the integral. Suppose is the approximation to , and the first order ETD method is given by Higher order ETD schemes can also be found [7, 9, 10]. Some preliminary numerical studies conducted by a group of material scientists at the Penn State University have indicated that the higher order ETD schemes can be several orders of magnitudes faster than the low-order semi-implicit methods in some simulations of microstructure evolution [8].

In essence, for nonlinear time dependent equations, the ETD scheme provides a systematic coupling of the explicit treatment of nonlinearities and the implicit and possibly exact integration of the stiff linear part of the equation, while achieving high accuracy and maintaining good stability. The exponential integrator of (8) is constructed based on the idea of Lie group methods on manifolds. In general, (8) can be written as where . In addition, should be locally Lipschitz function (to ensure existence and uniqueness). The classical fourth order explicit RK scheme of (11) is

According to a classical result of Hausdorff, the solution of (11) is is the solution of the initial value problem where , and . Taking as the corresponding number and applying the Runge-Kutta method to (14), we can get the numerical solution of (14). So the numerical solution of (11) is obtained. The corresponding order exponential integrator of (11) is for and it is explicit provided that the underlying RK scheme is explicit. is the order truncation error of the . An explicit fourth order exponential integrator of (11) is as follows: As the Lie group method on manifolds, the exponential integrator method of the general nonlinear differential equation has the same accuracy as the corresponding Runge-Kutta methods. In the subsequent sections, the fourth order explicit exponential integrator (EEI) method and the corresponding explicit Runge-Kutta (ERK) method are applied to the diffusion equations.

3. Numerical Experiments I: The Quasilinear Convection Diffusion Equation

To the quasilinear convection diffusion equation the equation is the one-dimensional quasilinear parabolic differential equation, which is known as Burgers’ equation, where is a constant representing the kinematic viscosity of the fluid. Burgers’ equation first appeared in a paper by Bateman and he gave a special solution of (18). In remarkable series of papers from 1939 to 1965, Burgers’ investigated various aspects of turbulence and he used it as a model in studies of turbulence. Cole studied the general properties of Burgers’ equation and outlined some of its various applications [1517]. With the initial condition , , and , the exact solution of (18) is and is the Bessel function. Recently some numerical methods of Burgers’ equation have been proposed [1517].

Burgers’ equation can be transformed into the nonlinear ordinary differential equation by discretization in the spacial direction. The spacial derivative of Burger’ equation can be discretized as follows: . . So the variable can satisfy Supposing that , (20) can be written as where , , , and ,   . EEI method and classical ERK method are applied to (21).

In Table 1, we compare the stability property of ERK method and EEI methods with . The convergent steps with different spacial step are shown in Table 1. The convergent step lengths are obtained by computing three hundred steps to decide whether the step length is convergent. From Table 1, we can see that EEI method has better stability than ERK method. In Table 2, we compare the global errors of ERK method and EEI method with and . The global error is defined as , where is the exact solution at . They have the same accuracy. In Figure 1, the numerical solutions of Burgers’ equation at are obtained by ERK method and the EEI method with , , and . Figure 1(a) was obtained by ERK method and Figure 1(b) was obtained by EEI method. From Figure 1, we can get that they have the same accuracy.

4. Numerical Experiments II: The Allen-Cahn Equation

To the Allen-Cahn equation with and initial condition . Suppose that , and (23) can be discretized into where and . EEI method and classical ERK method are applied to (24).

In Table 3, we compare the stability property of ERR method and EEI method by solving (24). When is different, the convergent step with different spacial step are shown in Table 3. The convergent step lengths are obtained by computing three hundred steps to decide whether the step length is convergent. From Table 3, we can see that the EEI method has better stability than ERK method. In Table 4, we compare the errors of the numerical solution at with obtained by ERK method and the numerical solution obtained by EEI method. The error is up to at different times; the error is . We can conclude that they have the same accuracy. In Figure 2, the numerical solutions of the Allen-Cahn equation in were obtained by ERK method and EEI method with , , and . Figure 2(a) was obtained by ERK method. Figure 2(b) was obtained by EEI method. From Figure 2, we can get that the numerical results are the same. We can also conclude that the two methods have the same accuracy.

5. Conclusions

In this paper, the EEI method based on the idea of the Lie group method was proposed. The EEI method and the corresponding ERK method were applied to Burgers’ equation and Allen-Cahn equation. Numerical results showed that the EEI method has better stability than the corresponding ERK methods. The two methods have the same accuracy. It is obvious that the explicit EER method is better than the ERK method in computing some stiff ordinary differential equations.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgment

This work was supported by the National Natural Science Foundation of China (11161017 and 10905022).